scholarly journals Microsatellites as Agents of Adaptive Change: An RNA-Seq-Based Comparative Study of Transcriptomes from Five Helianthus Species

Symmetry ◽  
2021 ◽  
Vol 13 (6) ◽  
pp. 933
Author(s):  
Chathurani Ranathunge ◽  
Sreepriya Pramod ◽  
Sébastien Renaut ◽  
Gregory L. Wheeler ◽  
Andy D. Perkins ◽  
...  

Mutations that provide environment-dependent selective advantages drive adaptive divergence among species. Many phenotypic differences among related species are more likely to result from gene expression divergence rather than from non-synonymous mutations. In this regard, cis-regulatory mutations play an important part in generating functionally significant variation. Some proposed mechanisms that explore the role of cis-regulatory mutations in gene expression divergence involve microsatellites. Microsatellites exhibit high mutation rates achieved through symmetric or asymmetric mutation processes and are abundant in both coding and non-coding regions in positions that could influence gene function and products. Here we tested the hypothesis that microsatellites contribute to gene expression divergence among species with 50 individuals from five closely related Helianthus species using an RNA-seq approach. Differential expression analyses of the transcriptomes revealed that genes containing microsatellites in non-coding regions (UTRs and introns) are more likely to be differentially expressed among species when compared to genes with microsatellites in the coding regions and transcripts lacking microsatellites. We detected a greater proportion of shared microsatellites in 5′UTRs and coding regions compared to 3′UTRs and non-coding transcripts among Helianthus spp. Furthermore, allele frequency differences measured by pairwise FST at single nucleotide polymorphisms (SNPs), indicate greater genetic divergence in transcripts containing microsatellites compared to those lacking microsatellites. A gene ontology (GO) analysis revealed that microsatellite-containing differentially expressed genes are significantly enriched for GO terms associated with regulation of transcription and transcription factor activity. Collectively, our study provides compelling evidence to support the role of microsatellites in gene expression divergence.

Author(s):  
Chathurani Ranathunge ◽  
Sreepriya Pramod ◽  
Sébastien Renaut ◽  
Gregory Wheeler ◽  
Andy Perkins ◽  
...  

Mutations that provide environment dependent selective advantages drive adaptive divergence among species. Many phenotypic differences among related species are more likely to result from gene expression divergence rather than from non-synonymous mutations. In this regard, cis-regulatory mutations play an important part in generating functionally significant variation. Some proposed mechanisms that explore the role of cis-regulatory mutations in gene expression divergence involve microsatellites. Microsatellites exhibit high mutation rates and are abundant in both coding and non-coding regions and could influence gene function and products. Here we tested the hypothesis that microsatellites contribute to gene expression divergence among species with 50 individuals from nine closely related Helianthus species using an RNA-seq approach. Differential expression analyses of the transcriptomes revealed that genes containing microsatellites in non-coding regions (UTRs and introns) are more likely to be differentially expressed among species when compared to genes with microsatellites in the coding regions and transcripts lacking microsatellites. We detected a greater proportion of shared microsatellites in 5’UTRs and coding regions compared to 3’UTRs and non-coding transcripts among Helianthus spp. Further, allele frequency differences measured by pairwise FST at single nucleotide polymorphisms (SNPs), indicate greater genetic divergence in transcripts containing microsatellites compared to those lacking microsatellites. A gene ontology (GO) analysis revealed that microsatellite-containing differentially expressed genes are significantly enriched for GO terms associated with regulation of transcription and transcription factor activity. Collectively, our study provides compelling evidence to support the role of microsatellites in gene expression divergence.


2021 ◽  
Vol 12 ◽  
Author(s):  
Jie Liu ◽  
Yan Zhou ◽  
Xin Hu ◽  
Jingchao Yang ◽  
Qiuxia Lei ◽  
...  

The developmental complexity of muscle arises from elaborate gene regulation. Long non-coding RNAs (lncRNAs) play critical roles in muscle development through the regulation of transcription and post-transcriptional gene expression. In chickens, previous studies have focused on the lncRNA profile during the embryonic periods, but there are no studies that explore the profile from the embryonic to post-hatching period. Here, we reconstructed 14,793 lncRNA transcripts and identified 2,858 differentially expressed lncRNA transcripts and 4,282 mRNAs from 12-day embryos (E12), 17-day embryos (E17), 1-day post-hatch chicks (D1), 14-day post-hatch chicks (D14), 56-day post-hatch chicks (D56), and 98-day post-hatch chicks (D98), based on our published RNA-seq datasets. We performed co-expression analysis for the differentially expressed lncRNAs and mRNAs, using STEM, and identified two profiles with opposite expression trends: profile 4 with a downregulated pattern and profile 21 with an upregulated pattern. The cis- and trans-regulatory interactions between the lncRNAs and mRNAs were predicted within each profile. Functional analysis of the lncRNA targets showed that lncRNAs in profile 4 contributed to the cell proliferation process, while lncRNAs in profile 21 were mainly involved in metabolism. Our work highlights the lncRNA profiles involved in the development of chicken breast muscle and provides a foundation for further experiments on the role of lncRNAs in the regulation of muscle development.


2018 ◽  
Author(s):  
Ana Llopart ◽  
Evgeny Brud ◽  
Nikale Pettie ◽  
Josep M. Comeron

ABSTRACTInteractions among divergent elements of transcriptional networks from different species can lead to misexpression in hybrids through regulatory incompatibilities, some with the potential to generate sterility. Genes with male-biased expression tend to be overrepresented among genes misexpressed in hybrid males. While the possible contribution of faster-male evolution to this misexpression has been explored, the role of the hemizygous X chromosome (i.e., the dominance theory for transcriptomes) remains yet to be determined. Here we study genome-wide patterns of gene expression in females and males of Drosophila yakuba and D. santomea and their hybrids. We used attached-X stocks to specifically test the dominance theory, and we uncovered a significant contribution of recessive alleles on the X chromosome to hybrid misexpression. Our analysis of gene expression patterns suggests that there is a contribution of weakly deleterious regulatory mutations to gene expression divergence in the sex towards which the expression is biased. In the opposite sex (e.g., genes with female-biased expression analyzed in male transcriptomes), we detect stronger selective constraints on gene expression divergence. Although genes with high degree of male-biased expression show a clear signal of faster-X evolution for gene expression divergence, we also detected slower-X evolution of gene expression in other gene classes (e.g. female-biased genes) that is mediated by significant decreases of cis- and trans-regulatory divergence. The distinct behavior of X-linked genes with high degree of male-biased expression is consistent with these genes experiencing a higher incidence of positively selected regulatory mutations than their autosomal counterparts. We propose that both dominance theory and faster-X evolution of gene expression may be major contributors to hybrid misexpression and possibly the large X-effect in these species.


2018 ◽  
Vol 27 (5) ◽  
pp. 1188-1199 ◽  
Author(s):  
Chathurani Ranathunge ◽  
Gregory L. Wheeler ◽  
Melody E. Chimahusky ◽  
Meaghan M. Kennedy ◽  
Jesse I. Morrison ◽  
...  

2021 ◽  
Vol 80 (Suppl 1) ◽  
pp. 12.2-12
Author(s):  
I. Muller ◽  
M. Verhoeven ◽  
H. Gosselt ◽  
M. Lin ◽  
T. De Jong ◽  
...  

Background:Tocilizumab (TCZ) is a monoclonal antibody that binds to the interleukin 6 receptor (IL-6R), inhibiting IL-6R signal transduction to downstream inflammatory mediators. TCZ has shown to be effective as monotherapy in early rheumatoid arthritis (RA) patients (1). However, approximately one third of patients inadequately respond to therapy and the biological mechanisms underlying lack of efficacy for TCZ remain elusive (1). Here we report gene expression differences, in both whole blood and peripheral blood mononuclear cells (PBMC) RNA samples between early RA patients, categorized by clinical TCZ response (reaching DAS28 < 3.2 at 6 months). These findings could lead to identification of predictive biomarkers for TCZ response and improve RA treatment strategies.Objectives:To identify potential baseline gene expression markers for TCZ response in early RA patients using an RNA-sequencing approach.Methods:Two cohorts of RA patients were included and blood was collected at baseline, before initiating TCZ treatment (8 mg/kg every 4 weeks, intravenously). DAS28-ESR scores were calculated at baseline and clinical response to TCZ was defined as DAS28 < 3.2 at 6 months of treatment. In the first cohort (n=21 patients, previously treated with DMARDs), RNA-sequencing (RNA-seq) was performed on baseline whole blood PAXgene RNA (Illumina TruSeq mRNA Stranded) and differential gene expression (DGE) profiles were measured between responders (n=14) and non-responders (n=7). For external replication, in a second cohort (n=95 therapy-naïve patients receiving TCZ monotherapy), RNA-seq was conducted on baseline PBMC RNA (SMARTer Stranded Total RNA-Seq Kit, Takara Bio) from the 2-year, multicenter, double-blind, placebo-controlled, randomized U-Act-Early trial (ClinicalTrials.gov identifier: NCT01034137) and DGE was analyzed between 84 responders and 11 non-responders.Results:Whole blood DGE analysis showed two significantly higher expressed genes in TCZ non-responders (False Discovery Rate, FDR < 0.05): urotensin 2 (UTS2) and caveolin-1 (CAV1). Subsequent analysis of U-Act-Early PBMC DGE showed nine differentially expressed genes (FDR < 0.05) of which expression in clinical TCZ non-responders was significantly higher for eight genes (MTCOP12, ZNF774, UTS2, SLC4A1, FECH, IFIT1B, AHSP, and SPTB) and significantly lower for one gene (TND2P28M). Both analyses were corrected for baseline DAS28-ESR, age and gender. Expression of UTS2, with a proposed function in regulatory T-cells (2), was significantly higher in TCZ non-responders in both cohorts. Furthermore, gene ontology enrichment analysis revealed no distinct gene ontology or IL-6 related pathway(s) that were significantly different between TCZ-responders and non-responders.Conclusion:Several genes are differentially expressed at baseline between responders and non-responders to TCZ therapy at 6 months. Most notably, UTS2 expression is significantly higher in TCZ non-responders in both whole blood as well as PBMC cohorts. UTS2 could be a promising target for further analyses as a potential predictive biomarker for TCZ response in RA patients in combination with clinical parameters (3).References:[1]Bijlsma JWJ, Welsing PMJ, Woodworth TG, et al. Early rheumatoid arthritis treated with tocilizumab, methotrexate, or their combination (U-Act-Early): a multicentre, randomised, double-blind, double-dummy, strategy trial. Lancet. 2016;388(10042):343-55.[2]Bhairavabhotla R, Kim YC, Glass DD, et al. Transcriptome profiling of human FoxP3+ regulatory T cells. Human Immunology. 2016;77(2):201-13.[3]Gosselt HR, Verhoeven MMA, Bulatovic-Calasan M, et al. Complex machine-learning algorithms and multivariable logistic regression on par in the prediction of insufficient clinical response to methotrexate in rheumatoid arthritis. Journal of Personalized Medicine. 2021;11(1).Disclosure of Interests:None declared


2021 ◽  
Vol 15 (1) ◽  
Author(s):  
Weitong Cui ◽  
Huaru Xue ◽  
Lei Wei ◽  
Jinghua Jin ◽  
Xuewen Tian ◽  
...  

Abstract Background RNA sequencing (RNA-Seq) has been widely applied in oncology for monitoring transcriptome changes. However, the emerging problem that high variation of gene expression levels caused by tumor heterogeneity may affect the reproducibility of differential expression (DE) results has rarely been studied. Here, we investigated the reproducibility of DE results for any given number of biological replicates between 3 and 24 and explored why a great many differentially expressed genes (DEGs) were not reproducible. Results Our findings demonstrate that poor reproducibility of DE results exists not only for small sample sizes, but also for relatively large sample sizes. Quite a few of the DEGs detected are specific to the samples in use, rather than genuinely differentially expressed under different conditions. Poor reproducibility of DE results is mainly caused by high variation of gene expression levels for the same gene in different samples. Even though biological variation may account for much of the high variation of gene expression levels, the effect of outlier count data also needs to be treated seriously, as outlier data severely interfere with DE analysis. Conclusions High heterogeneity exists not only in tumor tissue samples of each cancer type studied, but also in normal samples. High heterogeneity leads to poor reproducibility of DEGs, undermining generalization of differential expression results. Therefore, it is necessary to use large sample sizes (at least 10 if possible) in RNA-Seq experimental designs to reduce the impact of biological variability and DE results should be interpreted cautiously unless soundly validated.


Viruses ◽  
2021 ◽  
Vol 13 (2) ◽  
pp. 244 ◽  
Author(s):  
Antonio Victor Campos Coelho ◽  
Rossella Gratton ◽  
João Paulo Britto de Melo ◽  
José Leandro Andrade-Santos ◽  
Rafael Lima Guimarães ◽  
...  

HIV-1 infection elicits a complex dynamic of the expression various host genes. High throughput sequencing added an expressive amount of information regarding HIV-1 infections and pathogenesis. RNA sequencing (RNA-Seq) is currently the tool of choice to investigate gene expression in a several range of experimental setting. This study aims at performing a meta-analysis of RNA-Seq expression profiles in samples of HIV-1 infected CD4+ T cells compared to uninfected cells to assess consistently differentially expressed genes in the context of HIV-1 infection. We selected two studies (22 samples: 15 experimentally infected and 7 mock-infected). We found 208 differentially expressed genes in infected cells when compared to uninfected/mock-infected cells. This result had moderate overlap when compared to previous studies of HIV-1 infection transcriptomics, but we identified 64 genes already known to interact with HIV-1 according to the HIV-1 Human Interaction Database. A gene ontology (GO) analysis revealed enrichment of several pathways involved in immune response, cell adhesion, cell migration, inflammation, apoptosis, Wnt, Notch and ERK/MAPK signaling.


2012 ◽  
Vol 111 (suppl_1) ◽  
Author(s):  
Emma L Robinson ◽  
Syed Haider ◽  
Hillary Hei ◽  
Richard T Lee ◽  
Roger S Foo

Heart failure comprises of clinically distinct inciting causes but a consistent pattern of change in myocardial gene expression supports the hypothesis that unifying biochemical mechanisms underlie disease progression. The recent RNA-seq revolution has enabled whole transcriptome profiling, using deep-sequencing technologies. Up to 70% of the genome is now known to be transcribed into RNA, a significant proportion of which is long non-coding RNAs (lncRNAs), defined as polyribonucleotides of ≥200 nucleotides. This project aims to discover whether the myocardium expression of lncRNAs changes in the failing heart. Paired end RNA-seq from a 300-400bp library of ‘stretched’ mouse myocyte total RNA was carried out to generate 76-mer sequence reads. Mechanically stretching myocytes with equibiaxial stretch apparatus mimics pathological hypertrophy in the heart. Transcripts were assembled and aligned to reference genome mm9 (UCSC), abundance determined and differential expression of novel transcripts and alternative splice variants were compared with that of control (non-stretched) mouse myocytes. Five novel transcripts have been identified in our RNA-seq that are differentially expressed in stretched myocytes compared with non-stretched. These are regions of the genome that are currently unannotated and potentially are transcribed into non-coding RNAs. Roles of known lncRNAs include control of gene expression, either by direct interaction with complementary regions of the genome or association with chromatin remodelling complexes which act on the epigenome.Changes in expression of genes which contribute to the deterioration of the failing heart could be due to the actions of these novel lncRNAs, immediately suggesting a target for new pharmaceuticals. Changes in the expression of these novel transcripts will be validated in a larger sample size of stretched myocytes vs non-stretched myocytes as well as in the hearts of transverse aortic constriction (TAC) mice vs Sham (surgical procedure without the aortic banding). In vivo investigations will then be carried out, using siLNA antisense technology to silence novel lncRNAs in mice.


Sign in / Sign up

Export Citation Format

Share Document