scholarly journals Ecological divergence in sympatry causes gene misregulation in hybrids

2019 ◽  
Author(s):  
Joseph A. McGirr ◽  
Christopher H. Martin

AbstractEcological speciation occurs when reproductive isolation evolves as a byproduct of adaptive divergence between populations. However, it is unknown whether divergent ecological selection on gene regulation can directly cause reproductive isolation. Selection favoring regulatory divergence between species could result in gene misregulation in F1 hybrids and ultimately lower hybrid fitness. We combined 58 resequenced genomes with 124 transcriptomes to test this hypothesis in a young, sympatric radiation of Cyprinodon pupfishes endemic to San Salvador Island, Bahamas, which consists of a dietary generalist and two novel trophic specialists – a molluscivore and a scale-eater. We found more differential gene expression between closely related sympatric specialists than between allopatric generalist populations separated by 1000 km. Intriguingly, 9.6% of genes that were differentially expressed between sympatric species were also misregulated in their F1 hybrids. Consistent with divergent ecological selection causing misregulation, a subset of these genes were in highly differentiated genomic regions and enriched for functions important for trophic specialization, including head, muscle, and brain development. These regions also included genes that showed evidence of hard selective sweeps and were significantly associated with oral jaw length – the most rapidly diversifying skeletal trait in this radiation. Our results indicate that divergent ecological selection in sympatry can cause hybrid gene misregulation which may act as a primary reproductive barrier between nascent species.SignificanceIt is unknown whether the same genes that regulate ecological traits can simultaneously contribute to reproductive barriers between species. We measured gene expression in two trophic specialist species of Cyprinodon pupfishes that rapidly diverged from a generalist ancestor. We found genes differentially expressed between species that also showed extreme expression levels in their hybrid offspring. Many of these genes showed signs of selection and have putative effects on the development of traits that are important for ecological specialization. This suggests that genetic variants contributing to adaptive trait divergence between parental species negatively interact to cause hybrid gene misregulation, potentially producing unfit hybrids. Such loci may be important barriers to gene flow during the early stages of speciation, even in sympatry.

eLife ◽  
2019 ◽  
Vol 8 ◽  
Author(s):  
Jukka-Pekka Verta ◽  
Felicity C Jones

Regulation of gene expression is thought to play a major role in adaptation, but the relative importance of cis- and trans- regulatory mechanisms in the early stages of adaptive divergence is unclear. Using RNAseq of threespine stickleback fish gill tissue from four independent marine-freshwater ecotype pairs and their F1 hybrids, we show that cis-acting (allele-specific) regulation consistently predominates gene expression divergence. Genes showing parallel marine-freshwater expression divergence are found near to adaptive genomic regions, show signatures of natural selection around their transcription start sites and are enriched for cis-regulatory control. For genes with parallel increased expression among freshwater fish, the quantitative degree of cis- and trans-regulation is also highly correlated across populations, suggesting a shared genetic basis. Compared to other forms of regulation, cis-regulation tends to show greater additivity and stability across different genetic and environmental contexts, making it a fertile substrate for the early stages of adaptive evolution.


2018 ◽  
Author(s):  
Jukka-Pekka Verta ◽  
Felicity C. Jones

AbstractRegulation of gene expression is thought to play a major role in adaptation but the relative importance of cis- and trans-regulatory mechanisms in the early stages of adaptive divergence is unclear. Using RNAseq of threespine stickleback fish gill tissue from four independent marine-freshwater ecotype pairs and their F1 hybrids, we show that cis-acting (allele-specific) regulation consistently predominates gene expression divergence. Genes showing parallel marine-freshwater expression divergence are found near to adaptive genomic regions, show signatures of natural selection around their transcription start sites and are enriched for cis-regulatory control. For genes with parallel increased expression among freshwater fish, the quantitative degree of cis- and trans-regulation is also highly correlated across populations, suggesting a shared genetic basis. Compared to other forms of regulation, cis-regulation tends to show greater additivity and stability across different genetic and environmental contexts, making it a fertile substrate for the early stages of adaptive evolution.


2019 ◽  
Vol 59 (4) ◽  
pp. 900-911 ◽  
Author(s):  
M Tobler ◽  
N Barts ◽  
R Greenway

Abstract Mitochondria have been known to be involved in speciation through the generation of Dobzhansky–Muller incompatibilities, where functionally neutral co-evolution between mitochondrial and nuclear genomes can cause dysfunction when alleles are recombined in hybrids. We propose that adaptive mitochondrial divergence between populations can not only produce intrinsic (Dobzhansky–Muller) incompatibilities, but could also contribute to reproductive isolation through natural and sexual selection against migrants, post-mating prezygotic isolation, as well as by causing extrinsic reductions in hybrid fitness. We describe how these reproductive isolating barriers can potentially arise through adaptive divergence of mitochondrial function in the absence of mito-nuclear coevolution, a departure from more established views. While a role for mitochondria in the speciation process appears promising, we also highlight critical gaps of knowledge: (1) many systems with a potential for mitochondrially-mediated reproductive isolation lack crucial evidence directly linking reproductive isolation and mitochondrial function; (2) it often remains to be seen if mitochondrial barriers are a driver or a consequence of reproductive isolation; (3) the presence of substantial gene flow in the presence of mito-nuclear incompatibilities raises questions whether such incompatibilities are strong enough to drive speciation to completion; and (4) it remains to be tested how mitochondrial effects on reproductive isolation compare when multiple mechanisms of reproductive isolation coincide. We hope this perspective and the proposed research plans help to inform future studies of mitochondrial adaptation in a manner that links genotypic changes to phenotypic adaptations, fitness, and reproductive isolation in natural systems, helping to clarify the importance of mitochondria in the formation and maintenance of biological diversity.


2017 ◽  
Author(s):  
Joseph A. McGirr ◽  
Christopher H. Martin

AbstractParallel evolution of gene expression commonly underlies convergent niche specialization, but parallel changes in expression could also underlie divergent specialization. We investigated divergence in gene expression and whole-genome genetic variation across three sympatric Cyprinodon pupfishes endemic to San Salvador Island, Bahamas. This recent radiation consists of a generalist and two derived specialists adapted to novel niches – a ‘scale-eater’ and a ‘snail-eater.’ We sampled total mRNA from all three species at two early developmental stages and compared gene expression with whole-genome genetic differentiation among all three species in 42 resequenced genomes. 80% of genes that were differentially expressed between snail-eaters and generalists were up or downregulated in the same direction between scale-eaters and generalists; however, there were no fixed variants shared between species underlying these parallel changes in expression. Genes showing parallel evolution of expression were enriched for effects on metabolic processes, whereas genes showing divergent expression were enriched for effects on cranial skeleton development and pigment biosynthesis, reflecting the most divergent phenotypes observed between specialist species. Our findings reveal that even divergent niche specialists may exhibit convergent adaptation to higher trophic levels through shared genetic pathways. This counterintuitive result suggests that parallel evolution in gene expression can accompany divergent ecological speciation during adaptive radiation.Impact SummaryAdaptations that result in unique forms of ecological specialization are central to research in evolutionary biology, yet little is known about their molecular foundations. We combined transcriptome sequencing with whole-genome divergence scans to study the molecular evolution of two specialist Cyprinodon pupfish species – a ‘scale-eater’ and a ‘snail-eater’ – that rapidly diverged from a sympatric generalist ancestor within the last 10,000 years. While parallel evolution of gene expression driving convergent niche specialization seems common, we present, to our knowledge, the first example of significant parallel changes in expression coinciding with divergent niche specialization. 80% of genes that were differentially expressed between snail-eaters and generalists showed the same direction of expression in scale-eaters relative to generalists. Furthermore, parallel evolution in expression seem to be controlled by unique genetic variants in each specialist species. Genes showing parallel changes in expression were enriched for metabolic processes that may facilitate adaptation to a higher trophic level, while genes showing divergent expression likely shape the striking morphological differences between specialists. These findings contribute to a more nuanced understanding of convergent adaptations that arise during speciation, and highlight how species can evolve similar expression profiles adapted to divergent niches.


2020 ◽  
Vol 26 (29) ◽  
pp. 3619-3630
Author(s):  
Saumya Choudhary ◽  
Dibyabhaba Pradhan ◽  
Noor S. Khan ◽  
Harpreet Singh ◽  
George Thomas ◽  
...  

Background: Psoriasis is a chronic immune mediated skin disorder with global prevalence of 0.2- 11.4%. Despite rare mortality, the severity of the disease could be understood by the accompanying comorbidities, that has even led to psychological problems among several patients. The cause and the disease mechanism still remain elusive. Objective: To identify potential therapeutic targets and affecting pathways for better insight of the disease pathogenesis. Method: The gene expression profile GSE13355 and GSE14905 were retrieved from NCBI, Gene Expression Omnibus database. The GEO profiles were integrated and the DEGs of lesional and non-lesional psoriasis skin were identified using the affy package in R software. The Kyoto Encyclopaedia of Genes and Genomes pathways of the DEGs were analyzed using clusterProfiler. Cytoscape, V3.7.1 was utilized to construct protein interaction network and analyze the interactome map of candidate proteins encoded in DEGs. Functionally relevant clusters were detected through Cytohubba and MCODE. Results: A total of 1013 genes were differentially expressed in lesional skin of which 557 were upregulated and 456 were downregulated. Seven dysregulated genes were extracted in non-lesional skin. The disease gene network of these DEGs revealed 75 newly identified differentially expressed gene that might have a role in development and progression of the disease. GO analysis revealed keratinocyte differentiation and positive regulation of cytokine production to be the most enriched biological process and molecular function. Cytokines -cytokine receptor was the most enriched pathways. Among 1013 identified DEGs in lesional group, 36 DEGs were found to have altered genetic signature including IL1B and STAT3 which are also reported as hub genes. CCNB1, CCNA2, CDK1, IL1B, CXCL8, MKI 67, ESR1, UBE2C, STAT1 and STAT3 were top 10 hub gene. Conclusion: The hub genes, genomic altered DEGs and other newly identified differentially dysregulated genes would improve our understanding of psoriasis pathogenesis, moreover, the hub genes could be explored as potential therapeutic targets for psoriasis.


2020 ◽  
Vol 15 ◽  
Author(s):  
Chen-An Tsai ◽  
James J. Chen

Background: Gene set enrichment analyses (GSEA) provide a useful and powerful approach to identify differentially expressed gene sets with prior biological knowledge. Several GSEA algorithms have been proposed to perform enrichment analyses on groups of genes. However, many of these algorithms have focused on identification of differentially expressed gene sets in a given phenotype. Objective: In this paper, we propose a gene set analytic framework, Gene Set Correlation Analysis (GSCoA), that simultaneously measures within and between gene sets variation to identify sets of genes enriched for differential expression and highly co-related pathways. Methods: We apply co-inertia analysis to the comparisons of cross-gene sets in gene expression data to measure the costructure of expression profiles in pairs of gene sets. Co-inertia analysis (CIA) is one multivariate method to identify trends or co-relationships in multiple datasets, which contain the same samples. The objective of CIA is to seek ordinations (dimension reduction diagrams) of two gene sets such that the square covariance between the projections of the gene sets on successive axes is maximized. Simulation studies illustrate that CIA offers superior performance in identifying corelationships between gene sets in all simulation settings when compared to correlation-based gene set methods. Result and Conclusion: We also combine between-gene set CIA and GSEA to discover the relationships between gene sets significantly associated with phenotypes. In addition, we provide a graphical technique for visualizing and simultaneously exploring the associations of between and within gene sets and their interaction and network. We then demonstrate integration of within and between gene sets variation using CIA and GSEA, applied to the p53 gene expression data using the c2 curated gene sets. Ultimately, the GSCoA approach provides an attractive tool for identification and visualization of novel associations between pairs of gene sets by integrating co-relationships between gene sets into gene set analysis.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Ying Xie ◽  
Xiaofeng Hang ◽  
Wensheng Xu ◽  
Jing Gu ◽  
Yuanjing Zhang ◽  
...  

Abstract Background Most of the biological functions of circular RNAs (circRNAs) and the potential underlying mechanisms in hepatocellular carcinoma (HCC) have not yet been discovered. Methods In this study, using circRNA expression data from HCC tumor tissues and adjacent tissues from the Gene Expression Omnibus database, we identified out differentially expressed circRNAs and verified them by qRT-PCT. Functional experiments were performed to evaluate the effects of circFAM13B in HCC in vitro and in vivo. Results We found that circFAM13B was the most significantly differentially expressed circRNA in HCC tissue. Subsequently, in vitro and in vivo studies also demonstrated that circFAM13B promoted the proliferation of HCC. Further studies revealed that circFAM13B, a sponge of miR-212, is involved in the regulation of E2F5 gene expression by competitively binding to miR-212, inhibits the activation of the P53 signalling pathway, and promotes the proliferation of HCC cells. Conclusions Our findings revealed the mechanism underlying the regulatory role played by circFAM13B, miR-212 and E2F5 in HCC. This study provides a new theoretical basis and novel target for the clinical prevention and treatment of HCC.


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 22 (1) ◽  
Author(s):  
Rowan AlEjielat ◽  
Anas Khaleel ◽  
Amneh H. Tarkhan

Abstract Background Ankylosing spondylitis (AS) is a rare inflammatory disorder affecting the spinal joints. Although we know some of the genetic factors that are associated with the disease, the molecular basis of this illness has not yet been fully elucidated, and the genes involved in AS pathogenesis have not been entirely identified. The current study aimed at constructing a gene network that may serve as an AS gene signature and biomarker, both of which will help in disease diagnosis and the identification of therapeutic targets. Previously published gene expression profiles of 16 AS patients and 16 gender- and age-matched controls that were profiled on the Illumina HumanHT-12 V3.0 Expression BeadChip platform were mined. Patients were Portuguese, 21 to 64 years old, were diagnosed based on the modified New York criteria, and had Bath Ankylosing Spondylitis Disease Activity Index scores > 4 and Bath Ankylosing Spondylitis Functional Index scores > 4. All patients were receiving only NSAIDs and/or sulphasalazine. Functional enrichment and pathway analysis were performed to create an interaction network of differentially expressed genes. Results ITM2A, ICOS, VSIG10L, CD59, TRAC, and CTLA-4 were among the significantly differentially expressed genes in AS, but the most significantly downregulated genes were the HLA-DRB6, HLA-DRB5, HLA-DRB4, HLA-DRB3, HLA-DRB1, HLA-DQB1, ITM2A, and CTLA-4 genes. The genes in this study were mostly associated with the regulation of the immune system processes, parts of cell membrane, and signaling related to T cell receptor and antigen receptor, in addition to some overlaps related to the IL2 STAT signaling, as well as the androgen response. The most significantly over-represented pathways in the data set were associated with the “RUNX1 and FOXP3 which control the development of regulatory T lymphocytes (Tregs)” and the “GABA receptor activation” pathways. Conclusions Comprehensive gene analysis of differentially expressed genes in AS reveals a significant gene network that is involved in a multitude of important immune and inflammatory pathways. These pathways and networks might serve as biomarkers for AS and can potentially help in diagnosing the disease and identifying future targets for treatment.


2021 ◽  
Vol 22 (14) ◽  
pp. 7624
Author(s):  
Mohammad Saeed ◽  
Alejandro Ibáñez-Costa ◽  
Alejandra María Patiño-Trives ◽  
Laura Muñoz-Barrera ◽  
Eduardo Collantes Estévez ◽  
...  

Objectives: This study employed genetic and functional analyses using OASIS meta-analysis of multiple existing GWAS and gene-expression datasets to identify novel SLE genes. Methods: Four hundred and ten genes were mapped using SNIPPER to 30 SLE GWAS loci and investigated for expression in three SLE GEO-datasets and the Cordoba GSE50395-dataset. Blood eQTL for significant SNPs in SLE loci and STRING for functional pathways of differentially expressed genes were used. Confirmatory qPCR on SLE monocytes was performed. The entire 12p11 locus was investigated for genetic association using two additional GWAS. Expression of 150 genes at this locus was assessed. Based on this significance, qPCRs for DNM1L and KRAS were performed. Results: Fifty genes were differentially expressed in at least two SLE GEO-datasets, with all probes directionally aligned. DDX11, an RNA helicase involved in genome stability, was downregulated in both GEO and Cordoba datasets. The most significant SNP, rs3741869 in OASIS locus 12p11.21, containing DDX11, was a cis-eQTL regulating DDX11 expression. DDX11 was found repressed. The entire 12p11 locus showed three association peaks. Gene expression in GEO datasets identified DNM1L and KRAS, besides DDX11. Confirmatory qPCR validated DNM1L as an SLE susceptibility gene. DDX11, DNM1L and KRAS interact with each other and multiple known SLE genes including STAT1/STAT4 and major components of IFN-dependent gene expression, and are responsible for signal transduction of cytokines, hormones, and growth-factors, deregulation of which is involved in SLE-development. Conclusion: A genomic convergence approach with OASIS analysis of multiple GWAS and expression datasets identified DDX11 and DNM1L as novel SLE-genes, the expression of which is altered in monocytes from SLE patients. This study lays the foundation for understanding the pathogenic involvement of DDX11 and DNM1L in SLE by identifying them using a systems-biology approach, while the 12p11 locus harboring these genes was previously missed by four independent GWAS.


Sign in / Sign up

Export Citation Format

Share Document