PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “Bioinformatics analysis”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25Linked to original sources

Global Genomic Surveillance.

Global genomic surveillance has emerged as a foundational pillar of public health in the twenty-first century, enabling real-time tracking of pathogen evolution and informing outbreak response. This chapter examines the strategic architecture of global genomic surveillance, focusing on its application to arboviruses such as chikungunya virus (CHIKV). It explores the integration of genomic data with epidemiological, clinical, and environmental information within a One Health framework, while addressing critical challenges in governance, equity, and interoperability. The discussion covers the entire genomic surveillance workflow, from sample collection and sequencing to bioinformatic analysis and phylogenetic inference, and highlights the transformative role of artificial intelligence (AI) in predictive surveillance. By analyzing global initiatives, operational barriers, and emerging technologies, this chapter underscores the necessity of sustainable, equitable, and interoperable genomic systems to proactively address current and future infectious disease threats.

Humans↗

Local RNA target structure influences siRNA efficacy: systematic analysis of intentionally designed binding regions.

Contradictory reports in the literature have emphasised either the sequence of small interfering RNAs (siRNA) or the structure of their target molecules to be the major determinant of the efficiency of RNA interference (RNAi) approaches. In the present study, we analyse systematically the contributions of these parameters to siRNA activity by using deliberately designed mRNA constructs. The siRNA target sites were included in well-defined structural elements rendering them either highly accessible or completely involved in stable base-pairing. Furthermore, complementary sequence elements and various hairpins with different stem lengths and designs were used as target sites. Only one of the strands of the siRNA duplex was found to be capable of silencing via its respective target site, indicating that thermodynamic characteristics intrinsic to the siRNA strands are a basic determinant of siRNA activity. A significant obstruction of gene silencing by the same siRNA, however, was observed to be caused by structural features of the substrate RNA. Bioinformatic analysis of the mRNA structures suggests a direct correlation between the extent of gene-knockdown and the local free energy in the target region. Our findings indicate that, although a favourable siRNA sequence is a necessary prerequisite for efficient RNAi, complex target structures may limit the applicability even of carefully chosen siRNAs.

Animals↗

Transcriptome analysis of frontal cortex in alcohol-preferring and nonpreferring rats.

Although it is widely accepted that alcohol abuse and alcoholism have a significant genetic component of risk, the identities of the genes themselves remain obscure. To illuminate such potential genetic contributions, DNA macroarrays were used to probe for differences in normative cortical gene expression between rat strains genetically selected for alcohol self-administration preference, AA (Alko, alcohol) and P (Indiana, preferring), or avoidance, ANA (Alko, nonalcohol) and NP (Indiana, nonpreferring). Among 1,176 genes studied, six demonstrated confirmable, differential expression following comparison of ethanol-naive AA and ANA rats. Specifically, the mRNA level for metabotropic glutamate receptor 3 (mGluR3) was down-regulated in the AA vs. ANA lines. In contrast, calcium channel subunit alpha2delta1 (cacna2d1), vesicle-associated membrane protein 2 (VAMP2), syntaxin 1 (both syntaxin 1a and 1b; STX1a and STX1b), and syntaxin binding protein (MUNC-18) mRNAs were found to be increased in frontal cortex following comparison of AA with ANA animals. Bioinformatic analysis of these molecular targets showed that mGluR3 and cacna2d1 fall within chromosomal locations reported to be alcohol-related by the Collaborative Study on the Genetics of Alcoholism (COGA) as well as quantitative trait loci (QTL) studies. To determine further whether these differences were strain specific, the above-mentioned genes were compared in ethanol-preferring (P) and -nonpreferring (NP) selected lines. VAMP2 was the only gene that displayed statistically different mRNA levels in a comparison of P and NP rats. In conclusion, the altered cortical gene expression illuminated here would have the effect of altering neurotransmitter release in AA rats (compared with ANA rats). Such alterations, however, might not be a universal characteristic of all animal models of alcohol abuse and will also require further investigation in post-mortem human samples.

Alcohol Drinking↗

MDA-MB-435 cells are derived from M14 melanoma cells--a loss for breast cancer, but a boon for melanoma research.

BACKGROUND: The tissue of origin of the cell line MDA-MB-435 has been a matter of debate since analysis of DNA microarray data led Ross et al. (2000, Nat Genet 24(3):227-235) to suggest they might be of melanocyte origin due to their similarity to melanoma cell lines. We have previously shown that MDA-MB-435 cells maintained in multiple laboratories are of common origin to those used by Ross et al. and concluded that MDA-MB-435 cells are not a representative model for breast cancer. We could not determine, however, whether the melanoma-like properties of the MDA-MB-435 cell line are the result of misclassification or due to transdifferention to a melanoma-like phenotype. METHODS: We used karyotype, comparative genomic hybridization (CGH), and microsatalite polymorphism analyses, combined with bioinformatics analysis of gene expression and single nucleotide polymorphism (SNP) data, to test the hypothesis that the MDA-MB-435 cell line is derived from the melanoma cell line M14. RESULTS: We show that the MDA-MB-435 and M14 cell lines are essentially identical with respect to cytogenetic characteristics as well as gene expression patterns and that the minor differences found can be explained by phenotypic and genotypic clonal drift. CONCLUSIONS: All currently available stocks of MDA-MB-435 cells are derived from the M14 melanoma cell line and can no longer be considered a model of breast cancer. These cells are still a valuable system for the study of cancer metastasis and the extensive literature using these cells since 1982 represent a valuable new resource for the melanoma research community.

Breast Neoplasms↗

Expression analysis of a tyrosinase promoter sequence in zebrafish.

Sequence comparisons and functional analysis of the 5' upstream regions of tyrosinase genes have revealed the importance of cis-regulatory elements acting to control the spatiotemporal expression of tyrosinase in the melanocytes and retinal pigmented epithelium of developing embryos. To date there are no reports addressing the control of tyrosinase gene transcription in zebrafish, a vertebrate model organism of increasing importance. To exploit the tyrosinase gene as a marker in zebrafish we set out to clone its promoter and analyse its regulation during embryogenesis. Amplification of a zebrafish tyrosinase complementary DNA fragment by reverse transcriptase polymerase chain reaction allowed us to isolate and sequence a 1041 nt genomic DNA fragment that includes a transcription initiation site and 73 nt of the open reading frame. Bioinformatic analysis of this genomic sequence revealed five E-box motifs, including one CATGTG type E-box present in a putative initiation region. These are conserved positive regulatory elements in vertebrate tyrosinase promoters. We show that a region of 814 nt upstream from the translation start site of the zebrafish tyrosinase gene can drive expression in retinal pigmented epithelium in transiently transgenic zebrafish embryos but that its activity is not restricted to melanin-producing cells. This region is unable to drive transcription in human melanoma cell lines. Ectopic expression from this zebrafish tyrosinase promoter fragment is probably due to the absence of positive and negative cis-regulatory elements, such as a tyrosinase distal element, which is known to function as a pigment cell-specific enhancer.

Animals↗

Functional specificity of shuttling hnRNPs revealed by genome-wide analysis of their RNA binding profiles.

Nab2, Npl3, and Nab4/Hrp1 are essential RNA binding proteins of the shuttling hnRNP class that are required for the efficient export of mRNA. To characterize the in vivo transcript specificity of these proteins, we identified their mRNA binding partners using a microarray-based assay. Each of the three proteins was coimmunoprecipitated with many different mRNA transcripts. Interestingly, each protein exhibits preferential associations with a distinct set of mRNAs. Notably, some of these appear to denote specific functional classes. For example, the ribosomal protein mRNAs and other highly expressed transcripts significantly favor association with Npl3 over Nab2, and Nab4/Hrp1 is strongly enriched with transcripts required for amino acid metabolism. Significantly, nab4 mutants showed a striking, desensitized growth phenotype when exposed to amino acid stress conditions suggesting a biological consequence to the associations we observed. Supporting the hypothesis that these proteins display transcript specificity, we identified a unique 7-nucleotide sequence overrepresented in the transcripts highly associated with Nab2 and Nab4/Hrp1 using the REDUCE algorithm. Validating our approach, our bioinformatics analysis correctly identified the known binding site for Nab4/Hrp1. These specialized associations of the hnRNP proteins of Saccharomyces cerevisiae suggest the opportunity to regulate the processing of particular transcripts between transcription and translation.

Genome, Fungal↗

Investigation of the catalytic triad of arylamine N-acetyltransferases: essential residues required for acetyl transfer to arylamines.

The NATs (arylamine N-acetyltransferases) are a well documented family of enzymes found in both prokaryotes and eukaryotes. NATs are responsible for the acetylation of a range of arylamine, arylhydrazine and hydrazine compounds. We present here an investigation into the catalytic triad of residues (Cys-His-Asp) and other structural features of NATs using a variety of methods, including site-directed mutagenesis, X-ray crystallography and bioinformatics analysis, in order to investigate whether each of the residues of the catalytic triad is essential for catalytic activity. The catalytic triad of residues, Cys-His-Asp, is a well defined motif present in several families of enzymes. We mutated each of the catalytic residues in turn to investigate the role they play in catalysis. We also mutated a key residue, Gly126, implicated in acetyl-CoA binding, to examine the effects on acetylation activity. In addition, we have solved the structure of a C70Q mutant of Mycobacterium smegmatis NAT to a resolution of 1.45 A (where 1 A=0.1 nm). This structure confirms that the mutated protein is correctly folded, and provides a structural model for an acetylated NAT intermediate. Our bioinformatics investigation analysed the extent of sequence conservation between all eukaryotic and prokaryotic NAT enzymes for which sequence data are available. This revealed several new sequences, not yet reported, of NAT paralogues. Together, these studies have provided insight into the fundamental core of NAT enzymes, and the regions where sequence differences account for the functional diversity of this family. We have confirmed that each of the three residues of the triad is essential for acetylation activity.

Amino Acid Substitution↗

Genome-wide detection of tissue-specific alternative splicing in the human transcriptome.

We have developed an automated method for discovering tissue-specific regulation of alternative splicing through a genome-wide analysis of expressed sequence tags (ESTs). Using this approach, we have identified 667 tissue-specific alternative splice forms of human genes. We validated our muscle-specific and brain-specific splice forms for known genes. A high fraction (8/10) were reported to have a matching tissue specificity by independent studies in the published literature. The number of tissue-specific alternative splice forms is highest in brain, while eye-retina, muscle, skin, testis and lymph have the greatest enrichment of tissue-specific splicing. Overall, 10-30% of human alternatively spliced genes in our data show evidence of tissue-specific splice forms. Seventy-eight percent of our tissue-specific alternative splices appear to be novel discoveries. We present bioinformatics analysis of several tissue-specific splice forms, including automated protein isoform sequence and domain prediction, showing how our data can provide valuable insights into gene function in different tissues. For example, we have discovered a novel kidney-specific alternative splice form of the WNK1 gene, which appears to specifically disrupt its N-terminal kinase domain and may play a role in PHAII hypertension. Our database greatly expands knowledge of tissue-specific alternative splicing and provides a comprehensive dataset for investigating its functional roles and regulation in different human tissues.

Alternative Splicing↗

Identification, expression pattern, and subcellular location of human RIP isoforms.

XRIP alpha was identified as an adapter protein involved in RAP nuclear import. Several homologs were reported in mammal EST analysis, but the expression pattern and genomic organization of hRIP isoforms were not clarified yet. We isolated nine isoforms of hRIP from a premade human fetal brain library. hRIP alpha is the longest isoform with 219 residues, containing a N-terminal arginine-rich basic region, followed by an acidic region and two C-terminal Zn finger-like structures. hRIP beta deletes one Zn-finger-like structure. Three hRIP alpha isoforms and four hRIP sigma isoforms express truncated proteins due to frame shift. hRIP gamma isoforms lost the C-terminal Zn-finger-like structure. hRIP delta isoforms only contain the N-terminal arginine-rich basic region and the core sequence of the acidic region. The genomic organization of hRIP was identified by bioinformatic analysis. hRIP, containing seven exons, is located at human chromosome 17p13. hRIP was expressed in all 16 detected human tissues with a similar pattern. All EGFP-hRIP fusion proteins were located at the nucleus in the HEK293 cell. The two-polar molecular structure of hRIP might be involved in the basic cell function, and plays a role in the alternative nuclear ingress.

Amino Acid Motifs↗

Proteomic analysis of normal human urinary proteins isolated by acetone precipitation or ultracentrifugation.

BACKGROUND: Proteomic techniques have recently become available for large-scale protein analysis. The utility of these techniques in identification of urinary proteins is poorly defined. We constructed a proteome map of normal human urine as a reference protein database by using two differential fractionated techniques to isolate the proteins. METHODS: Proteins were isolated from urine obtained from normal human volunteers by acetone precipitation or ultracentrifugation, separated by two-dimensional polyacrylamide gel electrophoresis (2D-PAGE) and identified by matrix-assisted laser desorption ionization-time-of-flight (MALDI-TOF) mass spectrometry followed by peptide mass fingerprinting. RESULTS: A total of 67 protein forms of 47 unique proteins were identified, including transporters, adhesion molecules, complement, chaperones, receptors, enzymes, serpins, cell signaling proteins and matrix proteins. Acetone precipitated more acidic and hydrophilic proteins, whereas ultracentrifugation fractionated more basic, hydrophobic, and membrane proteins. Bioinformatic analysis predicted glycosylation to be the most common explanation for multiple forms of the same protein. CONCLUSIONS: Combining two differential isolation techniques magnified protein identification from human urine. Proteomic analysis of urinary proteins is a promising tool to study renal physiology and pathophysiology and to determine biomarkers of renal disease.

Acetone↗

Genome-wide analysis of defense-responsive genes in bacterial blight resistance of rice mediated by the recessive R gene xa13.

Defense responses triggered by dominant and recessive disease resistance (R) genes are presumed to be regulated by different molecular mechanisms. In order to characterize the genes activated in defense responses against bacterial blight mediated by the recessive R gene xa13, two pathogen-induced subtraction cDNA libraries were constructed using the resistant rice line IRBB13--which carries xa13--and its susceptible, near-isogenic, parental line IR24. Clustering analysis of expressed sequence tags (ESTs) identified 702 unique expressed sequences as being involved in the defense responses triggered by xa13; 16% of these are new rice ESTs. These sequences define 702 genes, putatively encoding a wide range of products, including defense-responsive genes commonly involved in different host-pathogen interactions, genes that have not previously been reported to be associated with pathogen-induced defense responses, and genes (38%) with no homology to previously described functional genes. In addition, R-like genes putatively encoding nucleotide-binding site/leucine rich repeat (NBS-LRR) and LRR receptor kinase proteins were observed to be induced in the disease resistance activated by xa13. A total of 568 defense-responsive ESTs were mapped to 588 loci on the rice molecular linkage map through bioinformatic analysis. About 48% of the mapped ESTs co-localized with quantitative trait loci (QTLs) for resistance to various rice diseases, including bacterial blight, rice blast, sheath blight and yellow mottle virus. Furthermore, some defense-responsive sequences were conserved at similar locations on different chromosomes. These results reveal the complexity of xa13-mediated resistance. The information obtained in this study provides a large source of candidate genes for understanding the molecular bases of defense responses activated by recessive R genes and of quantitative disease resistance.

Chromosomes, Plant↗

[Cluster analysis and identify significance of novel genes related to glioma in molecule classification].

OBJECTIVE: To screen differentially expressed genes in the development of human glioma and establish molecular classification of glioma preliminary based on gene expression using cDNA microarray. METHODS: Brain specimens were obtained from 18 patients with glioma, 10 males and 8 females, aged 14 approximately 62 with an average age of 44.4. The total RNAs of these glioma specimens and 2 specimens of donated brain of normal adults were extracted. BioStarH140S microarray (including 8347 old genes and 5592 novel genes) were adopted and hybridized with probes which were prepared from the total RNAs. Differentially expressed genes between the normal tissues and glioma tissues were assayed after scanning cDNA microarray with ScanArray 4000. Northern hybridization, and in situ hybridization (ISH) were used to identify the functions of novel genes. Those differentially expressed genes were studied with Hierarchical method and molecular classification of glioma was preliminarily carried out. RESULTS: Among the 13 939 target genes, there were 1200 (8.61%) differentially expressed genes and 395 (2.83%) novel genes. 348 genes were up-regulated and 852 genes were down-regulated in glioma. The results of bioinformatical analysis, Northern hybridization and ISH revealed that those novel genes were highly associated with glioma. There were multiple genes which were relevance to classification by Hierarchical method, such as MAP gene, cytoskeleton and matrix motility genes, etc. Molecular classification of glioma with Hierarchical cluster was in accordance with pathology and revealed internal essence in tumorigenesis and development. CONCLUSION: Multiple genes play important roles in development of glioma. cDNA microarray technology is a powerful technique in screening for differentially expressed genes between two different kinds of tissues.

Adolescent↗

Genetic analysis of the HAMP domain of the Aer aerotaxis sensor localizes flavin adenine dinucleotide-binding determinants to the AS-2 helix.

HAMP domains are signal transduction domains typically located between the membrane anchor and cytoplasmic signaling domain of the proteins in which they occur. The prototypical structure consists of two helical amphipathic sequences (AS-1 and AS-2) connected by a region of undetermined structure. The Escherichia coli aerotaxis receptor, Aer, has a HAMP domain and a PAS domain with a flavin adenine dinucleotide (FAD) cofactor that senses the intracellular energy level. Previous studies reported mutations in the HAMP domain that abolished FAD binding to the PAS domain. In this study, using random and site-directed mutagenesis, we identified the distal helix, AS-2, as the component of the HAMP domain that stabilizes FAD binding. AS-2 in Aer is not amphipathic and is predicted to be buried. Mutations in the sequence coding for the contiguous proximal signaling domain altered signaling by Aer but did not affect FAD binding. The V264M residue replacement in this region resulted in an inverted response in which E. coli cells expressing the mutant Aer protein were repelled by oxygen. Bioinformatics analysis of aligned HAMP domains indicated that the proximal signaling domain is conserved in other HAMP domains that are not involved in chemotaxis or aerotaxis. Only one null mutation was found in the coding sequence for the HAMP AS-1 and connector regions, suggesting that these are not active signal transduction sites. We consider a model in which the signal from FAD is transmitted across a PAS-HAMP interface to AS-2 or the proximal signaling domain.

Adenylyl Cyclases↗

DNA microarray analysis of striatal gene expression in symptomatic transgenic Huntington's mice (R6/2) reveals neuroinflammation and insulin associations.

Huntington's disease (HD) is an inherited, progressive neurodegenerative disorder caused by CAG repeat expansion in the gene that codes for the protein huntingtin. The underlying neuropathological events leading to the selectivity of striatal neuronal loss are unknown. However, the huntingtin mutation interferes at several levels of normal cell function. The complexity of this disease makes microarray analysis an appealing technique to begin the identification of common pathways that may contribute to the pathology. In this study, striatal tissue was extracted for gene expression profiling from wild-type and symptomatic transgenic Huntington mice (R6/2) expressing part of the human Huntington's disease gene. We interrogated a 15 K high-density mouse EST array not previously used for HD and identified 170 significantly differentially expressed ESTs in symptomatic R6/2 mice. Of the 80 genes with known function, 9 genes had previously been identified as altered in HD. 71 known genes were associated with HD for the first time. The data obtained from this study confirm and extend previous observations using DNA microarray techniques on genetic models for HD, revealing novel changes in expression in a number of genes not previously associated with HD. Further bioinformatic analysis, using software to construct biological association maps, focused attention on proteins such as insulin and TH1-mediated cytokines, suggesting that they may be important regulators of affected genes. These results may provide insight into the regulation and interaction of genes that contribute to adaptive and pathological processes involved in HD.

Animals↗

The CRAL/TRIO and GOLD domain protein CGR-1 promotes induction of vulval cell fates in Caenorhabditis elegans and interacts genetically with the Ras signaling pathway.

Ras-mediated signaling is necessary for the induction of vulval cell fates during Caenorhabditis elegans development. We identified cgr-1 by screening for suppressors of the ectopic vulval cell fates caused by a gain-of-function mutation of the let-60 ras gene. Analysis of two cgr-1 loss-of-function mutations indicates that cgr-1 positively regulates induction of vulval cell fates. cgr-1 is likely to function at a step in the Ras signaling pathway that is downstream of let-60, which encodes Ras, and upstream of lin-1, which encodes a transcription factor, if these genes function in a linear signaling pathway. These genetic studies are also consistent with the model that cgr-1 functions in a parallel pathway that promotes vulval cell fates. Localized expression studies suggest that cgr-1 functions cell autonomously to affect vulval cell fates. cgr-1 also functions early in development, since cgr-1 is necessary for larval viability. CGR-1 contains a CRAL/TRIO domain likely to bind a small hydrophobic ligand and a GOLD domain that may mediate interactions with proteins. A bioinformatic analysis revealed that there is a conserved family of CRAL/TRIO and GOLD domain-containing proteins that includes members from vertebrates and Drosophila. The analysis of cgr-1 identifies a novel in vivo function for a member of this family and a potential new regulator of Ras-mediated signaling.

Amino Acid Sequence↗

Novel polyphenol oxidase mined from a metagenome expression library of bovine rumen: biochemical properties, structural analysis, and phylogenetic relationships.

RL5, a gene coding for a novel polyphenol oxidase, was identified through activity screening of a metagenome expression library from bovine rumen microflora. Characterization of the recombinant protein produced in Escherichia coli revealed a multipotent capacity to oxidize a wide range of substrates (syringaldazine > 2,6-dimethoxyphenol > veratryl alcohol > guaiacol > tetramethylbenzidine > 4-methoxybenzyl alcohol > 2,2'-azino-bis(3-ethylbenzthiazoline-6-sulfonic acid) (ABTS) >> phenol red) over an unusually broad range of pH from 3.5 to 9.0. Apparent Km and kcat values for ABTS, syringaldazine, and 2,6-dimetoxyphenol obtained from steady-state kinetic measurements performed at 40 degrees C, pH 4.5, yielded values of 26, 0.43, and 0.45 microm and 18, 660, and 1175 s(-1), respectively. The Km values for syringaldazine and 2,6-dimetoxyphenol are up to 5 times lower, and the kcat values up to 40 times higher, than values previously reported for this class of enzyme. RL5 is a 4-copper oxidase with oxidation potential values of 745, 400, and 500 mV versus normal hydrogen electrode for the T1, T2, and T3 copper sites. A three-dimensional model of RL5 and site-directed mutants were generated to identify the copper ligands. Bioinformatic analysis of the gene sequence and the sequences and contexts of neighboring genes suggested a tentative phylogenetic assignment to the genus Bacteroides. Kinetic, electrochemical, and EPR analyses provide unequivocal evidence that the hypothetical proteins from Bacteroides thetaiotaomicron and from E. coli, which are closely related to the deduced protein encoded by the RL5 gene, are also multicopper proteins with polyphenol oxidase activity. The present study shows that these three newly characterized enzymes form a new family of functional multicopper oxidases with laccase activity related to conserved hypothetical proteins harboring the domain of unknown function DUF152 and suggests that some other of these proteins may also be laccases.

Amino Acid Sequence↗

Metabolic reconstruction of sulfur assimilation in the extremophile Acidithiobacillus ferrooxidans based on genome analysis.

BACKGROUND: Acidithiobacillus ferrooxidans is a gamma-proteobacterium that lives at pH2 and obtains energy by the oxidation of sulfur and iron. It is used in the biomining industry for the recovery of metals and is one of the causative agents of acid mine drainage. Effective tools for the study of its genetics and physiology are not in widespread use and, despite considerable effort, an understanding of its unusual physiology remains at a rudimentary level. Nearly complete genome sequences of A. ferrooxidans are available from two public sources and we have exploited this information to reconstruct aspects of its sulfur metabolism. RESULTS: Two candidate mechanisms for sulfate uptake from the environment were detected but both belong to large paralogous families of membrane transporters and their identification remains tentative. Prospective genes, pathways and regulatory mechanisms were identified that are likely to be involved in the assimilation of sulfate into cysteine and in the formation of Fe-S centers. Genes and regulatory networks were also uncovered that may link sulfur assimilation with nitrogen fixation, hydrogen utilization and sulfur reduction. Potential pathways were identified for sulfation of extracellular metabolites that may possibly be involved in cellular attachment to pyrite, sulfur and other solid substrates. CONCLUSIONS: A bioinformatic analysis of the genome sequence of A. ferrooxidans has revealed candidate genes, metabolic process and control mechanisms potentially involved in aspects of sulfur metabolism. Metabolic modeling provides an important preliminary step in understanding the unusual physiology of this extremophile especially given the severe difficulties involved in its genetic manipulation and biochemical analysis.

Adenosine Phosphosulfate↗

Dual RNA isolation from blood: an optimized protocol for host and bacterial RNA purification for dual RNA-sequencing analysis in whole blood sepsis samples.

Dual RNA-sequencing (dual RNA-seq) holds significant promise for deciphering bacterial virulence mechanisms during systemic infections. However, its application in sepsis research is hindered by technical challenges, including a low bacterial burden in blood and limited sample volumes and RNA yield from vulnerable populations, such as neonates. We developed an optimized protocol [dual RNA isolation from blood (DRIB)] for simultaneous stabilization, isolation and purification of high-quality host leukocyte and bacterial RNA from low-volume whole blood samples (0.5 ml). This protocol is compatible with clinical sample collection workflows and high-throughput RNA sequencing. The feasibility of DRIB for dual RNA-seq was validated using a pilot cohort of clinical adult sepsis samples, enabling the investigation of host-bacterial gene expression during sepsis. The DRIB protocol yielded 2.10-6.91 µg of total RNA per clinical sample in our pilot cohort. Dual-species ribosomal RNA (rRNA) depletion and RNA-seq generated 16.6-24.8 million filtered reads per sample, with 63±7% of reads uniquely mapped to host or bacterial sequences. Host genes accounted for 51-68% (8.4-10.9 million) reads, while 0.5-6.7% (79,496-789,808 reads) mapped to bacterial genomes. Bioinformatic analysis revealed that both shared and individual transcriptional patterns were identified in host and bacterial responses, including pathways related to immune metabolism and metal-ion binding. Our optimized DRIB protocol and RNA-seq pipeline effectively captured both host and bacterial RNA transcription in clinical sepsis samples. Expanding this approach to larger cohorts and varying disease timepoints will provide crucial new insights into host-bacterial gene co-expression dynamics in sepsis progression and outcomes.

Humans↗