Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “RNA sequencing analysis”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

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

At least 469 records · Page 26Linked to original sources

Independent segregation of the VP4 and the VP7 genes in bovine rotaviruses as confirmed by VP4 sequence analysis of G8 and G10 bovine rotavirus strains.

We determined the complete nucleotide sequences of the VP4 genes of five bovine rotavirus strains (A5, 61A, A44, B223 and KK3). The deduced VP4 amino acid sequence of strain A5 with G8 serotype specificity is 91.1% identical to that of strain NCDV of the serotype G6, and strain 61A with G10 serotype specificity has a VP4 amino acid sequence (95.2% identity) similar to that of the UK strain of the G6 serotype. In contrast, the VP4 amino acid sequences of strains A44, B223 and KK3 of the G10 serotype, isolated in Thailand, the U.S.A. and Japan respectively, have very similar sequences to each other, but less similarity (50 to 60%) to other group A rotavirus strains reported so far. Their VP4 genes are 2352 nucleotides in length and encode 772 amino acids, four amino acids fewer than the VP4 proteins of other animal rotaviruses. Thus, the presence of three different VP4 types (P types) and the independent segregation of G serotype and P type were postulated in bovine rotaviruses by VP4 sequence analysis.

Amino Acid Sequence↗

Characterization of neutralization epitopes on the VP7 surface protein of serotype G11 porcine rotaviruses.

Rotavirus strain A253, isolated from the faeces of a diarrhoeic piglet in Venezuela, was classified as serotype G11 by cross-neutralization studies and by comparison of the deduced amino acid sequence of the VP7 surface protein. The epitopes involved in neutralization of the two G11 porcine rotavirus strains A253 and YM were analysed using neutralization-resistant mutants selected with seven neutralizing monoclonal antibodies (MAbs), monotype-specific (M-) MAbs and serotype-specific (S-) MAbs, produced against VP7 of strain A253. Cross-neutralization tests and sequence analysis of the escape mutants selected from strains A253 and YM indicated the presence of two antigenic sites, one common to both M-MAbs and S-MAbs in region A (positions 87, 91 and 96) and the other defined by one S-MAb in region C (position 223). All A253 variants selected with M-MAbs and two S-MAbs, although having different amino acid substitutions, had a change at amino acid position 87, whereas YM variants involved residues 91 and 96, part of the same antigenic site. Compared to strain A253, the YM stain presents an amino acid substitution at position 87 and was not recognized by M-MAbs. These results suggest that in the VP7 of G11 serotype specificity, the amino acid at position 87 is an important component of a neutralization site associated with region A and the intraserotypic variation between strains A253 and YM may account for the selection of mutations at different positions by a single MAb.

Amino Acid Sequence↗

Systematic evaluation of metatranscriptomic differential gene expression in silico, in vitro, and in vivo enables elucidation of inter-species cross-feeding.

Metatranscriptomic (MTX) sequencing quantifies gene expression from the collective genomes of microbial communities (microbiomes), enabling assessment of functional activity rather than functional potential. While differential expression testing is instrumental to RNA-sequencing analysis, current metatranscriptomic approaches have been benchmarked only on simulated data and not under real operating conditions, resulting in a lack of standard practices. Here, we evaluate the performance of statistical differential expression methods on both simulated datasets and data collected from real bacterial 'mock communities' designed for this purpose. We assess the robustness of individual methods to organisms' low relative abundance, differential abundance, low prevalence, and transcription rate changes, showing that no existing methods perform adequately across all confounding conditions. We then apply the same approaches to metatranscriptomic datasets generated from gnotobiotic mice colonized with defined consortia of human bacterial strains and show that the method nominated by our mock community comparisons successfully inferred cross-feeding dynamics which were validated in vitro. We conclude that MTX method benchmarking on real, not simulated, datasets can and should optimize model implementation, enabling inference and validation of cross-feeding and other inter-species and host-microbe dynamics from in vivo studies.

Journal Article↗

Analytical description of finite size effects for RNA secondary structures.

The ensemble of RNA secondary structures of uniform sequences is studied analytically. We calculate the partition function for very long sequences and discuss how the crossover length, beyond which asymptotic scaling laws apply, depends on thermodynamic parameters. For realistic choices of parameters this length can be much longer than natural RNA molecules. This has to be taken into account when applying asymptotic theory to interpret experiments or numerical results.

Base Pairing↗

Rec-I-DCM3: a fast algorithmic technique for reconstructing large phylogenetic trees.

Phylogenetic trees are commonly reconstructed based on hard optimization problems such as maximum parsimony (MP) and maximum likelihood (ML). Conventional MP heuristics for producing phylogenetic trees produce good solutions within reasonable time on small datasets (up to a few thousand sequences), while ML heuristics are limited to smaller datasets (up to a few hundred sequences). However, since MP (and presumably ML) is NP-hard, such approaches do not scale when applied to large datasets. In this paper, we present a new technique called Recursive-Iterative-DCM3 (Rec-I-DCM3), which belongs to our family of Disk-Covering Methods (DCMs). We tested this new technique on ten large biological datasets ranging from 1,322 to 13,921 sequences and obtained dramatic speedups as well as significant improvements in accuracy (better than 99.99%) in comparison to existing approaches. Thus, high-quality reconstructions can be obtained for datasets at least ten times larger than was previously possible.

Algorithms↗

JViz.Rna--a Java tool for RNA secondary structure visualization.

Many tools have been developed for visualization of RNA secondary structures using a variety of techniques and output formats. However, each tool is typically limited to one or two of the visualization models discussed in this paper, supports only a single file format, and is tied to a specific platform. In order for structure prediction researchers to better understand the results of their algorithms and to enable life science researchers to interpret RNA structure easily, it is helpful to provide them with a flexible and powerful tool.jViz.Rna is a multiplatform visualization tool capable of displaying RNA secondary structures encoded in a variety of file formats. The same structure can be viewed using any of the models supported, including linked graph, circle graph, dot plot, and classical structure. Also, the output is dynamic and can easily be further manipulated by the user. In addition, any of the drawings produced can be saved in either the EPS or PNG file formats enabling easy usage in publications and presentations.

Computer Graphics↗

Visualization of barrier tree sequences.

Dynamical models that explain the formation of spatial structures of RNA molecules have reached a complexity that requires novel visualization methods that help to analyze the validity of these models. Here, we focus on the visualization of so-called folding landscapes of a growing RNA molecule. Folding landscapes describe the energy of a molecule as a function of its spatial configuration; thus they are huge and high dimensional. Their most salient features, however, are encapsulated by their so-called barrier tree that reflects the local minima and their connecting saddle points. For each length of the growing RNA chain there exists a folding landscape. We visualize the sequence of folding landscapes by an animation of the corresponding barrier trees. To generate the animation, we adapt the foresight layout with tolerance algorithm for general dynamic graph layout problems. Since it is very general, we give a detailed description of each phase: constructing a supergraph for the trees, layout of that supergraph using a modified DoT algorithm, and presentation techniques for the final animation.

Base Sequence↗

Host-dependent serum specificity to the V3 domain of HIV-1.

The previous studies of anti-HIV-1 humoral immune response have found the stable individual structural bias in the antibody production that reflects the changes in the immune functional network caused by HIV-1. The correlation between antigen structure, in particular the V3 domain of HIV-1, and the serum specificity to the structure remains unclear. To clarify the role of host individual factors in the serum specificity, we used the slightly modified HIV-1 serotyping methodological approach for the sera of the patients infected with HIV-1 with a different level of genetic similarity. Substantial discrepancy between serum specificities and antigen structure was found. Patients infected with HIV-1 carrying similar and identical V3 sequences had significantly different serum specificities. The opposite situation was often observed as well. The results of the study suggest that the influence of the V3 loop structure in the development of humoral immune response to the V3 loop is substantially modified in a patient-specific manner.

Antibody Formation↗

SLA2 is Associated With Immune evasion and Exhaustion of CD8+ T Cells in Gastric Cancer.

The Src-like adaptor 2 (SLA2) functions as a negative regulator of T cell receptor signalling. However, its involvement in the tumour microenvironment (TME) of gastric cancer (GC) remains unexplored. In this study, we found that SLA2 expression was significantly elevated in GC tissues, and a high level of SLA2 was associated with poor prognosis in GC patients. Bioinformatics analyses revealed a close association between SLA2 and TME in GC. Single-cell RNA sequencing analysis indicated that SLA2 was significantly enriched in CD8+ T cells in GC tissues. Functional validation demonstrated that SLA2 overexpression contributed to the exhaustion of CD8+ T cells by suppressing their proliferation, upregulating the expression of exhaustion markers, reducing the secretion of effector cytokines (IFN-γ and TNF-α) and impairing cytotoxic function. SLA2 knockdown in in vitro-generated exhausted CD8 T cells significantly alleviated T cell exhaustion. Mechanistically, we found that inverse promoter methylation and active histone marks (H3K27ac, H3K4me3 and H3K4me1) may regulate SLA2 expression. Our findings suggest that SLA2 may modulate the TME and promote immune evasion via CD8+ T cell exhaustion in GC.

Humans↗

The H3K27me3 reader GmLHP1 impairs Phytophthora sojae resistance by repressing ethylene precursor accumulation in soybean.

Phytophthora root rot, caused by Phytophthora sojae, is a devastating soilborne disease of soybean (Glycine max). However, the epigenetic regulation of soybean responses to P. sojae remains incompletely understood. Here, using genetic, molecular and biochemical approaches, we characterized the functions of LIKE HETEROCHROMATIN PROTEIN 1 (GmLHP1) and its downstream regulatory network. We demonstrated that GmLHP1, as a reader of H3K27me3, negatively regulates soybean resistance to P. sojae. GmLHP1 binds to H3K27me3 peptides in vitro and colocalizes with H3K27me3 marks genome-wide in vivo. The integrated chromatin immunoprecipitation sequencing and RNA sequencing analysis identified the ethylene biosynthesis pathway gene 1-AMINO-CYCLOPROPANE-1-CARBOXYLATE SYNTHASE 18 (GmACS18) as being enriched for H3K27me3 and bound by GmLHP1, leading to its transcriptional downregulation. Notably, GmLHP1 associates with the GmACS18 promoter by directly binding to AATTAA motifs and recognizing H3K27me3 marks. Moreover, GmACS18 enhances defense against P. sojae by accumulating the ethylene precursor 1-aminocyclopropane-1-carboxylic acid (ACC). Further analysis unveiled that recognition of H3K27me3 by GmLHP1 is essential for regulating soybean resistance to P. sojae through repressing GmACS18 transcription and decreasing ACC accumulation. Our findings reveal a novel epigenetic regulatory mechanism in which the H3K27me3 reader GmLHP1 blocks soybean resistance to P. sojae by repressing ethylene precursor ACC accumulation.

ACC accumulation↗

The Tartary Buckwheat FtMYB46-FtNRAMP3 Module Enhances Plant Lead and Cadmium Tolerance.

The presence of toxic heavy metals lead (Pb) and cadmium (Cd) in polluted soil damage crop production and consequently harms human and livestock health. Tartary buckwheat (Fagopyrum tataricum) is a potential model plant for heavy metal phytoremediation because of its valuable characteristics of high heavy metal tolerance and abundant biomass production. Here, we report that the Tartary buckwheat FtMYB46-FtNRAMP3 module enhances plant Pb and Cd tolerance. RNA sequencing analysis showed that Pb treatment specifically induced expression of FtNRAMP3, a member of the NRAMP (Natural Resistance-Associated Macrophage Protein) transporter gene family. Further cytological and biochemical analysis revealed that FtNRAMP3 was localised to the plasma membrane and significantly contributed to increased tolerance to Pb and Cd in yeast cells. Consistently, transgenic overexpression of FtNRAMP3 in Arabidopsis significantly increased plant tolerance to Pb and Cd applications, reducing Pb concentration but increasing Cd concentration in the overexpression transgenic plants. Subsequent yeast one-hybrid and electrophoretic mobility shift assays showed that the transcription factor FtMYB46 directly binds to the FtNRAMP3 promoter. Further, FtMYB46 promoted FtNRAMP3 expression and increased plant Pb and Cd tolerance. Overall, this study demonstrates the important role of the FtMYB46-FtNRAMP3 module and its potential value in the phytoremediation of Pb and Cd stress.

Cadmium↗

NAM and CUC3 boundary genes maintain shoot apical meristem viability and suppress the development of axillary shoot in rice seedlings.

Cell division and differentiation within the shoot apical meristem (SAM) are essential for the morphogenesis of aboveground plant organs. This study reveals that the boundary genes OsNAM and OsCUC3 collaboratively maintain SAM activity. Loss of function in both OsNAM and OsCUC3 during the fourth leaf stage reduced SAM size, with the osnam oscuc3 mutant exhibiting abnormal leaf number and morphology. Furthermore, OsNAM and OsCUC3 inhibited the growth of axillary shoots. In the osnam oscuc3 mutant, the number of new leaves decreased, while buds in the coleoptile and the axil of the first leaf developed into tillers. Since OsNAM and OsCUC3 are involved in regulating both SAM activity and the growth of lateral shoots, we examined their expression patterns at the base of the main shoot. β-Glucuronidase (GUS) reporter activity and GFP reporter lines demonstrated that OsNAM and OsCUC3 have distinct expression patterns. Specifically, OsNAM was expressed throughout the SAM, whereas OsCUC3 was expressed only at the base of the SAM, with its expression gradually decreasing as seedlings develop. RNA sequencing analysis showed that the expression of genes related to leaf epidermal cell development, cell wall components, and hormonal signal transduction was altered in response to the loss of function of OsNAM and OsCUC3. Therefore, the boundary genes OsNAM and OsCUC3 not only inhibit the growth of axillary shoots but also regulate the development of aboveground organs, including leaf morphology and number, by maintaining the SAM activity in the main shoot.

Meristem↗

Defective EV-mediated transport of SHH alters neural fate specification in EPM1 epilepsy.

The extracellular milieu, including extracellular vesicles (EVs), plays a pivotal role in brain development. In this study, we sought to elucidate the pathogenesis of progressive myoclonus epilepsy type 1 (EPM1), a disease caused by mutations in the CSTB gene, using cerebral organoids (COs) derived from patient cells. The results demonstrate that EPM1 COs display increased electrophysiological activity and disrupted excitatory/inhibitory (E/I) balance. Single-cell RNA sequencing analysis of ventral EPM1-COs revealed an abnormal specification of progenitor fate, with a shift toward dorsal neuron identities. We demonstrated that this misspecification is driven by a functional alteration of the ventral signaling niche, resulting from impaired EV dynamics and altered protein cargo. Mechanistically, we identified Sonic Hedgehog (SHH) as a direct physical interactor of CSTB and demonstrated that CSTB deficiency leads to reduced SHH content and secretion. Our findings establish CSTB as a safeguard of ventral patterning and identify the CSTB-SHH-EV axis as a potential therapeutic target for mitigating the E/I imbalance associated with EPM1.

Hedgehog Proteins↗

RR06 activates transcription of spr1996 and cbpA in Streptococcus pneumoniae.

Streptococcus pneumoniae colonizes at the nasopharynx of humans and is able to disseminate and cause various infections. The hallmark of pneumococcal disease is rapid bacterial replication in different tissue sites leading to intense inflammation. The genetic basis of pneumococcal adaptation to different host niches remains sketchy. In this study, we investigated the regulatory effect of RR06, a response regulator protein, on gene expression of S. pneumoniae. Microarray and Northern blot analyses showed that RR06 is specifically required for transcription of spr1996 and cbpA. While the function of Spr1996 is unknown, CbpA has been well characterized as a surface-exposed protective antigen and a virulence factor of S. pneumoniae. A recombinant form of RR06 was able to bind to a 19-bp conserved sequence shared by the spr1996 and cbpA promoter regions. Furthermore, inactivation of rr06 resulted in loss of CbpA expression as detected by antibody staining and bacterial adhesion. CbpA expression was restored in trans by the intact rr06 gene. However, a mutant, RR06(D51A), with a point mutation in the aspartate residue at position 51 (a predicted major phosphorylation site) of RR06, completely abolished the CbpA expression, suggesting that RR06 phosphorylation is required for transcriptional activation of spr1996 and cbpA. Finally, inactivation of rr06 in additional pneumococcal strains also led to the loss of CbpA expression. These data implicate that RR06 activates the expression of spr1996 and cbpA in many other pneumococcal strains.

Amino Acid Sequence↗

Eight genes in region 5 that form an operon are essential for invasion of epithelial cells by Shigella flexneri 2a.

The 7-kb region 5 on the large 230-kb plasmid pMYSH6000 in Shigella flexneri 2a YSH6000 is one of the virulence-associated DNA segments required for the invasion of epithelial cells (C. Sasakawa, K. Kamata, T. Sakai, S. Makino, M. Yamada, N. Okada, and M. Yoshikawa, J. Bacteriol. 170:2480-2484, 1988). To elucidate the functional organization of region 5 and to determine the virulence-associated genes encoded by region 5, we performed insertion and deletion mutagenesis, DNA subcloning, and complete nucleotide sequencing of region 5 and found that region 5 contained 11 open reading frames (ORFs) named ORF-1 through ORF-11 which could be translated into proteins with molecular masses of 15.1, 47.5, 13.2, 33.0, 33.4, 24.2, 9.4, 28.5, 39.9, 9.1, and 10.4 kDa, respectively. Complementation tests of the 14 Tn5-induced noninvasive mutants of region 5 with the above plasmid constructs have indicated that region 5 consists of an operon and that ORF-2 through ORF-9, but not ORF-1, ORF-10, and ORF-11, are essential for invasion, and 7 of 8 ORFs (ORF-2 and ORF-4 through ORF-9) and presumably the remaining ORF (ORF-3) are required for secretion of the Ipa proteins. The transcriptional organization, as determined by a promoter-proving vector, S1 nuclease protection, and primer extension RNA sequencing analysis revealed that region 5 is transcribed from a promoter located 47 bp upstream of the 5' end of ORF-2 for the 47.5-kDa protein and that the promoter activity identified was regulated by the virB gene, the transcriptional activator on the 230-kb plasmid.

Amino Acid Sequence↗

Rnall: an efficient algorithm for predicting RNA local secondary structural landscape in genomes.

BACKGROUND: The information of RNA local secondary structures (LSSs) can help retrieve biologically important motifs and study functions of RNA molecules. Most of the current RNA secondary structure prediction tools are not suitable for RNA LSS prediction on the genome scale due to high computational complexity. METHODS: We developed a new computer package Rnall based on a dynamic programming technique, which scans an RNA sequence with a sliding window and extracts all RNA LSSs with sizes no larger than the window size using the nearest neighbor thermodynamic parameters. The worst case running time of Rnall is O(W(3)L), where W is the window size and L is the query sequence length. In practice we observed a running time of O(W(2)L). We further introduced the concept of energy landscape for illustrating RNA LSS, which may facilitate RNA motif mining on the genomic scale. RESULTS: Rnall shows better prediction accuracy than two other prediction tools Lfold and Quickfold. Rnall is also applied to scan for RNA LSSs in three genomes, and the prediction maps well with known RNA motifs. CONCLUSIONS: Rnall is designed for RNA LSS prediction and together with the energy landscape, it has unique features that could be used for RNA structural motif mining. Rnall is freely available for download at http://digbio.missouri.edu/~wanx/Rnall or http://www.sysbio.muohio.edu/Rnall.

Algorithms↗

The skeletal muscle of aged male mice exhibits sustained growth regulatory transcriptional profile following glucocorticoid exposure compared with young males.

Excess glucocorticoids induce skeletal muscle myopathy by changing gene expression. Advanced age augments glucocorticoid-mediated muscle phenotypes, yet the transcriptional responses underlying those augmented phenotypes are unclear. The purpose of this study was to define the glucocorticoid-responsive transcriptome in young and aged muscle following both acute and more prolonged glucocorticoid treatment. Young (4-mo-old) or aged (24-mo-old) male mice were administered either an acute injection of dexamethasone (DEX) or vehicle or daily DEX or vehicle injections for 7 days. Muscles were harvested 6.5 h after the final or only injection. The tibialis anterior (TA) was selected for RNA sequencing analysis as DEX treatment lowered TA mass specifically in aged males. In silico analyses identified enriched pathways and transcription factors predicted to regulate DEX-sensitive genes. Acute DEX altered similar numbers of genes in young (950) versus aged males (913), although aged males had greater magnitudes of fold change. After 7 days of DEX treatment, aged muscle exhibited more DEGs compared with acute exposure (1,196 vs. 913), whereas young muscle exhibited fewer DEGs than after acute exposure (599 vs. 950). In aged males, glucocorticoid-sensitive genes were consistently enriched for growth regulatory processes across both time points, a pattern that was not evident in young males. Despite those age-associated transcriptional differences, the transcription factors predicted to regulate the glucocorticoid-sensitive genes were similar in young and aged males. These data expand our understanding into how aging modifies the transcriptional response to excess glucocorticoids in skeletal muscle.NEW & NOTEWORTHY Glucocorticoids promote mass loss in certain muscles with advanced age but not at younger ages. In a muscle whose mass is lost in response to elevated glucocorticoids only in advanced age in males, we show that glucocorticoids initiate a unique and exaggerated transcriptional profile after both acute exposure to the hormone and after prolonged treatment that is consistent with muscle atrophy. These findings expand our understanding of the effect primary aging has on glucocorticoid-induced atrophy in males.

Animals↗