1
|
Leclercq J, Torres-Paz J, Policarpo M, Agnès F, Rétaux S. Evolution of the regulation of developmental gene expression in blind Mexican cavefish. Development 2024; 151:dev202610. [PMID: 39007346 DOI: 10.1242/dev.202610] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2023] [Accepted: 07/08/2024] [Indexed: 07/16/2024]
Abstract
Developmental evolution and diversification of morphology can arise through changes in the regulation of gene expression or protein-coding sequence. To unravel mechanisms underlying early developmental evolution in cavefish of the species Astyanax mexicanus, we compared transcriptomes of surface-dwelling and blind cave-adapted morphs at the end of gastrulation. Twenty percent of the transcriptome was differentially expressed. Allelic expression ratios in cave X surface hybrids showed that cis-regulatory changes are the quasi-exclusive contributors to inter-morph variations in gene expression. Among a list of 108 genes with change at the cis-regulatory level, we explored the control of expression of rx3, which is a master eye gene. We discovered that cellular rx3 levels are cis-regulated in a cell-autonomous manner, whereas rx3 domain size depends on non-autonomous Wnt and Bmp signalling. These results highlight how uncoupled mechanisms and regulatory modules control developmental gene expression and shape morphological changes. Finally, a transcriptome-wide search for fixed coding mutations and differential exon use suggested that variations in coding sequence have a minor contribution. Thus, during early embryogenesis, changes in gene expression regulation are the main drivers of cavefish developmental evolution.
Collapse
Affiliation(s)
- Julien Leclercq
- Paris-Saclay Institute of Neuroscience, CNRS and University Paris-Saclay, 91400 Saclay, France
| | - Jorge Torres-Paz
- Paris-Saclay Institute of Neuroscience, CNRS and University Paris-Saclay, 91400 Saclay, France
| | - Maxime Policarpo
- Paris-Saclay Institute of Neuroscience, CNRS and University Paris-Saclay, 91400 Saclay, France
| | - François Agnès
- Paris-Saclay Institute of Neuroscience, CNRS and University Paris-Saclay, 91400 Saclay, France
| | - Sylvie Rétaux
- Paris-Saclay Institute of Neuroscience, CNRS and University Paris-Saclay, 91400 Saclay, France
| |
Collapse
|
2
|
Luo M, Hu J. Alternative splicing in parallel evolution and the evolutionary potential in sticklebacks. J Anim Ecol 2024. [PMID: 39056271 DOI: 10.1111/1365-2656.14157] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2023] [Accepted: 06/28/2024] [Indexed: 07/28/2024]
Abstract
Repeatability of adaptation to similar environments provides opportunity to evaluate the predictability of natural selection. While many studies have investigated gene expression differences between populations adapted to contrasting environments, the role of post-transcriptional processes such as alternative splicing has rarely been evaluated in the context of parallel adaptation. To address the aforementioned knowledge gap, we reanalysed transcriptomic data from three pairs of threespine stickleback (Gasterosteus aculeatus) ecotypes adapted to marine or freshwater environment. First, we identified genes with repeated expression or splicing divergence across ecotype pairs, and compared the genetic architecture and biological processes between parallelly expressed and parallelly spliced loci. Second, we analysed the extent to which parallel adaptation was reflected at gene expression and alternative splicing levels. Finally, we tested how the two axes of transcriptional variation differed in their potential for evolutionary change. Although both repeated differential splicing and differential expression across ecotype pairs showed tendency for parallel divergence, the degree of parallelism was lower for splicing than expression. Furthermore, parallel divergences in splicing and expression were likely to be associated with distinct cis-regulatory genetic variants and functionally unique set of genes. Finally, we found that parallelly spliced genes showed higher nucleotide diversity than parallelly expressed genes, indicating splicing is less susceptible to genetic variation erosion during parallel adaptation. Our results provide novel insight into the role of splicing in parallel adaptation, and underscore the contribution of splicing to the evolutionary potential of wild populations under environmental change.
Collapse
Affiliation(s)
- Man Luo
- Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, Institute of Biodiversity Science, Center for Evolutionary Biology, School of Life Sciences, Fudan University, Shanghai, China
| | - Juntao Hu
- Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, Institute of Biodiversity Science, Center for Evolutionary Biology, School of Life Sciences, Fudan University, Shanghai, China
| |
Collapse
|
3
|
Okude G, Yamasaki YY, Toyoda A, Mori S, Kitano J. Genome-wide analysis of histone modifications can contribute to the identification of candidate cis-regulatory regions in the threespine stickleback fish. BMC Genomics 2024; 25:685. [PMID: 38992624 PMCID: PMC11241946 DOI: 10.1186/s12864-024-10602-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2024] [Accepted: 07/08/2024] [Indexed: 07/13/2024] Open
Abstract
BACKGROUND Cis-regulatory mutations often underlie phenotypic evolution. However, because identifying the locations of promoters and enhancers in non-coding regions is challenging, we have fewer examples of identified causative cis-regulatory mutations that underlie naturally occurring phenotypic variations than of causative amino acid-altering mutations. Because cis-regulatory elements have epigenetic marks of specific histone modifications, we can detect cis-regulatory elements by mapping and analyzing them. Here, we investigated histone modifications and chromatin accessibility with cleavage under targets and tagmentation (CUT&Tag) and assay for transposase-accessible chromatin-sequencing (ATAC-seq). RESULTS Using the threespine stickleback (Gasterosteus aculeatus) as a model, we confirmed that the genes for which nearby regions showed active marks, such as H3K4me1, H3K4me3, and high chromatin accessibility, were highly expressed. In contrast, the expression levels of genes for which nearby regions showed repressive marks, such as H3K27me3, were reduced, suggesting that our chromatin analysis protocols overall worked well. Genomic regions with peaks of histone modifications showed higher nucleotide diversity within and between populations. By comparing gene expression in the gills of the marine and stream ecotypes, we identified several insertions and deletions (indels) with transposable element fragments in the candidate cis-regulatory regions. CONCLUSIONS Thus, mapping and analyzing histone modifications can help identify cis-regulatory elements and accelerate the identification of causative mutations in the non-coding regions underlying naturally occurring phenotypic variations.
Collapse
Affiliation(s)
- Genta Okude
- Ecological Genetics Laboratory, National Institute of Genetics, Yata 1111, Mishima, Shizuoka, 411-8540, Japan.
| | - Yo Y Yamasaki
- Ecological Genetics Laboratory, National Institute of Genetics, Yata 1111, Mishima, Shizuoka, 411-8540, Japan
| | - Atsushi Toyoda
- Comparative Genetics Laboratory, National Institute of Genetics, Mishima, Shizuoka, Japan
| | - Seiichi Mori
- Faculty of Economics, Gifu-Kyoritsu University, Ogaki, Gifu, Japan
| | - Jun Kitano
- Ecological Genetics Laboratory, National Institute of Genetics, Yata 1111, Mishima, Shizuoka, 411-8540, Japan.
| |
Collapse
|
4
|
Venu V, Harjunmaa E, Dreau A, Brady S, Absher D, Kingsley DM, Jones FC. Fine-scale contemporary recombination variation and its fitness consequences in adaptively diverging stickleback fish. Nat Ecol Evol 2024; 8:1337-1352. [PMID: 38839849 PMCID: PMC11239493 DOI: 10.1038/s41559-024-02434-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2023] [Accepted: 05/02/2024] [Indexed: 06/07/2024]
Abstract
Despite deep evolutionary conservation, recombination rates vary greatly across the genome and among individuals, sexes and populations. Yet the impact of this variation on adaptively diverging populations is not well understood. Here we characterized fine-scale recombination landscapes in an adaptively divergent pair of marine and freshwater populations of threespine stickleback from River Tyne, Scotland. Through whole-genome sequencing of large nuclear families, we identified the genomic locations of almost 50,000 crossovers and built recombination maps for marine, freshwater and hybrid individuals at a resolution of 3.8 kb. We used these maps to quantify the factors driving variation in recombination rates. We found strong heterochiasmy between sexes but also differences in recombination rates among ecotypes. Hybrids showed evidence of significant recombination suppression in overall map length and in individual loci. Recombination rates were lower not only within individual marine-freshwater-adaptive loci, but also between loci on the same chromosome, suggesting selection on linked gene 'cassettes'. Through temporal sampling along a natural hybrid zone, we found that recombinants showed traits associated with reduced fitness. Our results support predictions that divergence in cis-acting recombination modifiers, whose functions are disrupted in hybrids, may play an important role in maintaining differences among adaptively diverging populations.
Collapse
Affiliation(s)
- Vrinda Venu
- Friedrich Miescher Laboratory of the Max Planck Society, Tübingen, Germany.
- Los Alamos National Laboratory, New Mexico, NM, USA.
| | - Enni Harjunmaa
- Friedrich Miescher Laboratory of the Max Planck Society, Tübingen, Germany
- CeGAT GmbH, Tübingen, Germany
| | - Andreea Dreau
- Friedrich Miescher Laboratory of the Max Planck Society, Tübingen, Germany
- Evotec SE 'Campus Curie', Toulouse, France
| | - Shannon Brady
- Deptartment of Developmental Biology, Stanford University School of Medicine, Stanford, CA, USA
| | - Devin Absher
- HudsonAlpha Institute for Biotechnology, Huntsville, AL, USA
| | - David M Kingsley
- Deptartment of Developmental Biology, Stanford University School of Medicine, Stanford, CA, USA
- Howard Hughes Medical Institute, Chevy Chase, MD, USA
| | - Felicity C Jones
- Friedrich Miescher Laboratory of the Max Planck Society, Tübingen, Germany.
- Groningen Institute for Evolutionary Life Sciences (GELIFES), University of Groningen, Groningen, the Netherlands.
| |
Collapse
|
5
|
Singh PP, Reeves GA, Contrepois K, Papsdorf K, Miklas JW, Ellenberger M, Hu CK, Snyder MP, Brunet A. Evolution of diapause in the African turquoise killifish by remodeling the ancient gene regulatory landscape. Cell 2024; 187:3338-3356.e30. [PMID: 38810644 DOI: 10.1016/j.cell.2024.04.048] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2021] [Revised: 11/30/2023] [Accepted: 04/30/2024] [Indexed: 05/31/2024]
Abstract
Suspended animation states allow organisms to survive extreme environments. The African turquoise killifish has evolved diapause as a form of suspended development to survive a complete drought. However, the mechanisms underlying the evolution of extreme survival states are unknown. To understand diapause evolution, we performed integrative multi-omics (gene expression, chromatin accessibility, and lipidomics) in the embryos of multiple killifish species. We find that diapause evolved by a recent remodeling of regulatory elements at very ancient gene duplicates (paralogs) present in all vertebrates. CRISPR-Cas9-based perturbations identify the transcription factors REST/NRSF and FOXOs as critical for the diapause gene expression program, including genes involved in lipid metabolism. Indeed, diapause shows a distinct lipid profile, with an increase in triglycerides with very-long-chain fatty acids. Our work suggests a mechanism for the evolution of complex adaptations and offers strategies to promote long-term survival by activating suspended animation programs in other species.
Collapse
Affiliation(s)
| | - G Adam Reeves
- Department of Genetics, Stanford University, Stanford, CA, USA
| | - Kévin Contrepois
- Department of Genetics, Stanford University, Stanford, CA, USA; Stanford Cardiovascular Institute, Stanford University, Stanford, CA, USA
| | | | - Jason W Miklas
- Department of Genetics, Stanford University, Stanford, CA, USA
| | | | - Chi-Kuo Hu
- Department of Genetics, Stanford University, Stanford, CA, USA
| | - Michael P Snyder
- Department of Genetics, Stanford University, Stanford, CA, USA; Stanford Cardiovascular Institute, Stanford University, Stanford, CA, USA; Stanford Diabetes Research Center, Stanford University, Stanford, CA, USA
| | - Anne Brunet
- Department of Genetics, Stanford University, Stanford, CA, USA; Glenn Center for the Biology of Aging, Stanford University, Stanford, CA, USA; Wu Tsai Neurosciences Institute, Stanford University, Stanford, CA, USA; Chan Zuckerberg Biohub, San Francisco, San Francisco, CA, USA.
| |
Collapse
|
6
|
Innes PA, Goebl AM, Smith CCR, Rosenberger K, Kane NC. Gene expression and alternative splicing contribute to adaptive divergence of ecotypes. Heredity (Edinb) 2024; 132:120-132. [PMID: 38071268 PMCID: PMC10924094 DOI: 10.1038/s41437-023-00665-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2023] [Revised: 11/23/2023] [Accepted: 11/24/2023] [Indexed: 03/10/2024] Open
Abstract
Regulation of gene expression is a critical link between genotype and phenotype explaining substantial heritable variation within species. However, we are only beginning to understand the ways that specific gene regulatory mechanisms contribute to adaptive divergence of populations. In plants, the post-transcriptional regulatory mechanism of alternative splicing (AS) plays an important role in both development and abiotic stress response, making it a compelling potential target of natural selection. AS allows organisms to generate multiple different transcripts/proteins from a single gene and thus may provide a source of evolutionary novelty. Here, we examine whether variation in alternative splicing and gene expression levels might contribute to adaptation and incipient speciation of dune-adapted prairie sunflowers in Great Sand Dunes National Park, Colorado, USA. We conducted a common garden experiment to assess transcriptomic variation among ecotypes and analyzed differential expression, differential splicing, and gene coexpression. We show that individual genes are strongly differentiated for both transcript level and alternative isoform proportions, even when grown in a common environment, and that gene coexpression networks are disrupted between ecotypes. Furthermore, we examined how genome-wide patterns of sequence divergence correspond to divergence in transcript levels and isoform proportions and find evidence for both cis and trans-regulation. Together, our results emphasize that alternative splicing has been an underappreciated mechanism providing source material for natural selection at short evolutionary time scales.
Collapse
Affiliation(s)
- Peter A Innes
- Ecology and Evolutionary Biology Department, University of Colorado, Boulder, CO, USA.
| | - April M Goebl
- Ecology and Evolutionary Biology Department, University of Colorado, Boulder, CO, USA
- Research and Conservation Department, Denver Botanic Gardens, Denver, CO, USA
| | - Chris C R Smith
- Ecology and Evolutionary Biology Department, University of Colorado, Boulder, CO, USA
- Institute of Ecology and Evolution, University of Oregon, Eugene, OR, USA
| | - Kaylee Rosenberger
- Ecology and Evolutionary Biology Department, University of Colorado, Boulder, CO, USA
| | - Nolan C Kane
- Ecology and Evolutionary Biology Department, University of Colorado, Boulder, CO, USA
| |
Collapse
|
7
|
Benson CW, Sheltra MR, Huff DR. The genome of Salmacisia buchloëana, the parasitic puppet master pulling strings of sexual phenotypic monstrosities in buffalograss. G3 (BETHESDA, MD.) 2024; 14:jkad238. [PMID: 37847611 PMCID: PMC10849329 DOI: 10.1093/g3journal/jkad238] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/28/2023] [Revised: 10/04/2023] [Accepted: 10/06/2023] [Indexed: 10/19/2023]
Abstract
To complete its parasitic lifecycle, Salmacisia buchloëana, a biotrophic fungus, manipulates reproductive organ development, meristem determinacy, and resource allocation in its dioecious plant host, buffalograss (Bouteloua dactyloides; Poaceae). To gain insight into S. buchloëana's ability to manipulate its host, we sequenced and assembled the 20.1 Mb genome of S. buchloëana into 22 chromosome-level pseudomolecules. Phylogenetic analysis suggests that S. buchloëana is nested within the genus Tilletia and diverged from Tilletia caries and Tilletia walkeri ∼40 MYA. We find that S. buchloëana contains a novel chromosome arm with no syntenic relationship to other publicly available Tilletia genomes, and that genes on the novel arm are upregulated upon infection, suggesting that this unique chromosomal segment may have played a critical role in S. buchloëana's evolution and host specificity. Salmacisia buchloëana has one of the largest fractions of serine peptidases (1.53% of the proteome) and one of the highest GC contents (62.3%) in all classified fungi. Analysis of codon base composition indicated that GC content is controlled more by selective constraints than directional mutation, and that S. buchloëana has a unique bias for the serine codon UCG. Finally, we identify 3 inteins within the S. buchloëana genome, 2 of which are located in a gene often used in fungal taxonomy. The genomic and transcriptomic resources generated here will aid plant pathologists and breeders by providing insight into the extracellular components contributing to sex determination in dioecious grasses.
Collapse
Affiliation(s)
- Christopher W Benson
- Department of Plant Science, Pennsylvania State University, University Park, PA 16801, USA
- Intercollegiate Graduate Degree Program in Plant Biology, Pennsylvania State University, University Park, PA 16801, USA
| | - Matthew R Sheltra
- Department of Plant Science, Pennsylvania State University, University Park, PA 16801, USA
- Intercollegiate Graduate Degree Program in Plant Biology, Pennsylvania State University, University Park, PA 16801, USA
| | - David R Huff
- Department of Plant Science, Pennsylvania State University, University Park, PA 16801, USA
| |
Collapse
|
8
|
Durkin SM, Ballinger MA, Nachman MW. Tissue-specific and cis-regulatory changes underlie parallel, adaptive gene expression evolution in house mice. PLoS Genet 2024; 20:e1010892. [PMID: 38306396 PMCID: PMC10866503 DOI: 10.1371/journal.pgen.1010892] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2023] [Revised: 02/14/2024] [Accepted: 01/22/2024] [Indexed: 02/04/2024] Open
Abstract
Changes in gene regulation have long been appreciated as a driving force of adaptive evolution, however the relative contributions of cis- and trans-acting changes to gene regulation over short evolutionary timescales remain unclear. Instances of recent, parallel phenotypic evolution provide an opportunity to assess whether parallel patterns are seen at the level of gene expression, and to assess the relative contribution of cis- and trans- changes to gene regulation in the early stages of divergence. Here, we studied gene expression in liver and brown adipose tissue in two wild-derived strains of house mice that independently adapted to cold, northern environments, and we compared them to a strain of house mice from a warm, tropical environment. To investigate gene regulatory evolution, we studied expression in parents and allele-specific expression in F1 hybrids of crosses between warm-adapted and cold-adapted strains. First, we found that the different cold-adapted mice showed both unique and shared changes in expression, but that the proportion of shared changes (i.e. parallelism) was greater than expected by chance. Second, we discovered that expression evolution occurred largely at tissue-specific and cis-regulated genes, and that these genes were over-represented in parallel cases of evolution. Finally, we integrated the expression data with scans for selection in natural populations and found substantial parallelism in the two northern populations for genes under selection. Furthermore, selection outliers were associated with cis-regulated genes more than expected by chance; cis-regulated genes under selection influenced phenotypes such as body size, immune functioning, and activity level. These results demonstrate that parallel patterns of gene expression in mice that have independently adapted to cold environments are driven largely by tissue-specific and cis-regulatory changes, providing insight into the mechanisms of adaptive gene regulatory evolution at the earliest stages of divergence.
Collapse
Affiliation(s)
- Sylvia M. Durkin
- Museum of Vertebrate Zoology and Department of Integrative Biology, University of California, Berkeley, Berkeley, California, United States of America
| | - Mallory A. Ballinger
- Museum of Vertebrate Zoology and Department of Integrative Biology, University of California, Berkeley, Berkeley, California, United States of America
| | - Michael W. Nachman
- Museum of Vertebrate Zoology and Department of Integrative Biology, University of California, Berkeley, Berkeley, California, United States of America
| |
Collapse
|
9
|
Mañes-García J, Marco-Ferreres R, Beccari L. Shaping gene expression and its evolution by chromatin architecture and enhancer activity. Curr Top Dev Biol 2024; 159:406-437. [PMID: 38729683 DOI: 10.1016/bs.ctdb.2024.01.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/12/2024]
Abstract
Transcriptional regulation plays a pivotal role in orchestrating the intricate genetic programs governing embryonic development. The expression of developmental genes relies on the combined activity of several cis-regulatory elements (CREs), such as enhancers and silencers, which can be located at long linear distances from the genes that they regulate and that interact with them through establishment of chromatin loops. Mutations affecting their activity or interaction with their target genes can lead to developmental disorders and are thought to have importantly contributed to the evolution of the animal body plan. The income of next-generation-sequencing approaches has allowed identifying over a million of sequences with putative regulatory potential in the human genome. Characterizing their function and establishing gene-CREs maps is essential to decode the logic governing developmental gene expression and is one of the major challenges of the post-genomic era. Chromatin 3D organization plays an essential role in determining how CREs specifically contact their target genes while avoiding deleterious off-target interactions. Our understanding of these aspects has greatly advanced with the income of chromatin conformation capture techniques and fluorescence microscopy approaches to visualize the organization of DNA elements in the nucleus. Here we will summarize relevant aspects of how the interplay between CRE activity and chromatin 3D organization regulates developmental gene expression and how it relates to pathological conditions and the evolution of animal body plan.
Collapse
Affiliation(s)
| | | | - Leonardo Beccari
- Centro de Biología Molecular Severo Ochoa, CSIC-UAM, Madrid, Spain.
| |
Collapse
|
10
|
Kautt AF, Chen J, Lewarch CL, Hu C, Turner K, Lassance JM, Baier F, Bedford NL, Bendesky A, Hoekstra HE. Evolution of gene expression across brain regions in behaviourally divergent deer mice. Mol Ecol 2024:e17270. [PMID: 38263608 DOI: 10.1111/mec.17270] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2023] [Revised: 01/03/2024] [Accepted: 01/08/2024] [Indexed: 01/25/2024]
Abstract
The evolution of innate behaviours is ultimately due to genetic variation likely acting in the nervous system. Gene regulation may be particularly important because it can evolve in a modular brain-region specific fashion through the concerted action of cis- and trans-regulatory changes. Here, to investigate transcriptional variation and its regulatory basis across the brain, we perform RNA sequencing (RNA-Seq) on ten brain subregions in two sister species of deer mice (Peromyscus maniculatus and P. polionotus)-which differ in a range of innate behaviours, including their social system-and their F1 hybrids. We find that most of the variation in gene expression distinguishes subregions, followed by species. Interspecific differential expression (DE) is pervasive (52-59% of expressed genes), whereas the number of DE genes between sexes is modest overall (~3%). Interestingly, the identity of DE genes varies considerably across brain regions. Much of this modularity is due to cis-regulatory divergence, and while 43% of genes were consistently assigned to the same gene regulatory class across subregions (e.g. conserved, cis- or trans-regulatory divergence), a similar number were assigned to two or more different gene regulatory classes. Together, these results highlight the modularity of gene expression differences and divergence in the brain, which may be key to explain how the evolution of brain gene expression can contribute to the astonishing diversity of animal behaviours.
Collapse
Affiliation(s)
- Andreas F Kautt
- Department of Organismic & Evolutionary Biology, Department of Molecular & Cellular Biology, Center for Brain Science, Harvard University, Cambridge, Massachusetts, USA
| | - Jenny Chen
- Department of Organismic & Evolutionary Biology, Department of Molecular & Cellular Biology, Center for Brain Science, Harvard University, Cambridge, Massachusetts, USA
| | - Caitlin L Lewarch
- Department of Organismic & Evolutionary Biology, Department of Molecular & Cellular Biology, Center for Brain Science, Harvard University, Cambridge, Massachusetts, USA
| | - Caroline Hu
- Department of Organismic & Evolutionary Biology, Department of Molecular & Cellular Biology, Center for Brain Science, Harvard University, Cambridge, Massachusetts, USA
| | - Kyle Turner
- Department of Organismic & Evolutionary Biology, Department of Molecular & Cellular Biology, Center for Brain Science, Harvard University, Cambridge, Massachusetts, USA
| | - Jean-Marc Lassance
- Department of Organismic & Evolutionary Biology, Department of Molecular & Cellular Biology, Center for Brain Science, Harvard University, Cambridge, Massachusetts, USA
| | - Felix Baier
- Department of Organismic & Evolutionary Biology, Department of Molecular & Cellular Biology, Center for Brain Science, Harvard University, Cambridge, Massachusetts, USA
| | - Nicole L Bedford
- Department of Organismic & Evolutionary Biology, Department of Molecular & Cellular Biology, Center for Brain Science, Harvard University, Cambridge, Massachusetts, USA
| | - Andres Bendesky
- Department of Organismic & Evolutionary Biology, Department of Molecular & Cellular Biology, Center for Brain Science, Harvard University, Cambridge, Massachusetts, USA
| | - Hopi E Hoekstra
- Department of Organismic & Evolutionary Biology, Department of Molecular & Cellular Biology, Center for Brain Science, Harvard University, Cambridge, Massachusetts, USA
| |
Collapse
|
11
|
Tsouris A, Brach G, Schacherer J, Hou J. Non-additive genetic components contribute significantly to population-wide gene expression variation. CELL GENOMICS 2024; 4:100459. [PMID: 38190102 PMCID: PMC10794783 DOI: 10.1016/j.xgen.2023.100459] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/31/2023] [Revised: 09/19/2023] [Accepted: 11/09/2023] [Indexed: 01/09/2024]
Abstract
Gene expression variation, an essential step between genotype and phenotype, is collectively controlled by local (cis) and distant (trans) regulatory changes. Nevertheless, how these regulatory elements differentially influence gene expression variation remains unclear. Here, we bridge this gap by analyzing the transcriptomes of a large diallel panel consisting of 323 unique hybrids originating from genetically divergent Saccharomyces cerevisiae isolates. Our analysis across 5,087 transcript abundance traits showed that non-additive components account for 36% of the gene expression variance on average. By comparing allele-specific read counts in parent-hybrid trios, we found that trans-regulatory changes underlie the majority of gene expression variation in the population. Remarkably, most cis-regulatory variations are also exaggerated or attenuated by additional trans effects. Overall, we showed that the transcriptome is globally buffered at the genetic level mainly due to trans-regulatory variation in the population.
Collapse
Affiliation(s)
- Andreas Tsouris
- Université de Strasbourg, CNRS, GMGM UMR, 7156 Strasbourg, France
| | - Gauthier Brach
- Université de Strasbourg, CNRS, GMGM UMR, 7156 Strasbourg, France
| | - Joseph Schacherer
- Université de Strasbourg, CNRS, GMGM UMR, 7156 Strasbourg, France; Institut Universitaire de France (IUF), Paris, France.
| | - Jing Hou
- Université de Strasbourg, CNRS, GMGM UMR, 7156 Strasbourg, France.
| |
Collapse
|
12
|
Mirchandani CD, Shultz AJ, Thomas GWC, Smith SJ, Baylis M, Arnold B, Corbett-Detig R, Enbody E, Sackton TB. A Fast, Reproducible, High-throughput Variant Calling Workflow for Population Genomics. Mol Biol Evol 2024; 41:msad270. [PMID: 38069903 PMCID: PMC10764099 DOI: 10.1093/molbev/msad270] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2023] [Revised: 10/27/2023] [Accepted: 11/22/2023] [Indexed: 01/05/2024] Open
Abstract
The increasing availability of genomic resequencing data sets and high-quality reference genomes across the tree of life present exciting opportunities for comparative population genomic studies. However, substantial challenges prevent the simple reuse of data across different studies and species, arising from variability in variant calling pipelines, data quality, and the need for computationally intensive reanalysis. Here, we present snpArcher, a flexible and highly efficient workflow designed for the analysis of genomic resequencing data in nonmodel organisms. snpArcher provides a standardized variant calling pipeline and includes modules for variant quality control, data visualization, variant filtering, and other downstream analyses. Implemented in Snakemake, snpArcher is user-friendly, reproducible, and designed to be compatible with high-performance computing clusters and cloud environments. To demonstrate the flexibility of this pipeline, we applied snpArcher to 26 public resequencing data sets from nonmammalian vertebrates. These variant data sets are hosted publicly to enable future comparative population genomic analyses. With its extensibility and the availability of public data sets, snpArcher will contribute to a broader understanding of genetic variation across species by facilitating the rapid use and reuse of large genomic data sets.
Collapse
Affiliation(s)
- Cade D Mirchandani
- Department of Biomolecular Engineering, University of California Santa Cruz, Santa Cruz, CA 95064, USA
- Genomics Institute, University of California Santa Cruz, Santa Cruz, CA 95064, USA
| | - Allison J Shultz
- Ornithology Department, Natural History Museum of Los Angeles County, Los Angeles, CA 90007, USA
| | | | - Sara J Smith
- Informatics Group, Harvard University, Cambridge, MA, USA
- Biology, Mount Royal University, Calgary, AB T3E 6K6, Canada
| | - Mara Baylis
- Department of Biomolecular Engineering, University of California Santa Cruz, Santa Cruz, CA 95064, USA
- Genomics Institute, University of California Santa Cruz, Santa Cruz, CA 95064, USA
| | - Brian Arnold
- Department of Ecology and Evolutionary Biology, Princeton University, Princeton, NJ, USA
- Center for Statistics and Machine Learning, Princeton University, Princeton, NJ, USA
| | - Russ Corbett-Detig
- Department of Biomolecular Engineering, University of California Santa Cruz, Santa Cruz, CA 95064, USA
- Genomics Institute, University of California Santa Cruz, Santa Cruz, CA 95064, USA
| | - Erik Enbody
- Department of Biomolecular Engineering, University of California Santa Cruz, Santa Cruz, CA 95064, USA
| | | |
Collapse
|
13
|
Fulgione D, Maselli V, Rivieccio E, Aceto S, Salvemini M, Buglione M. Evolutionary Plasticity in Insular Lizard, Adapting over Reproduction, Metabolism, and Color Variation. BIOLOGY 2023; 12:1478. [PMID: 38132304 PMCID: PMC10740616 DOI: 10.3390/biology12121478] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/24/2023] [Revised: 11/22/2023] [Accepted: 11/27/2023] [Indexed: 12/23/2023]
Abstract
The Italian wall lizard (Podarcis siculus) living on islets exhibits a melanic skin coloration and a suite of adaptive traits lacking in nearby mainland populations. On islets, the unpredictable environmental conditions and highly fluctuating population densities are believed to have produced reversed island syndrome (RIS). Several physiological, behavioral, and life-history changes based on the RIS could result from positive selection on increased activity of melanocortins. We hypothesize that phenotypes on islets are the product of a plastic variation depending on the regulation of specific genes. Focusing on control systems that determine the insular-adapted phenotype, we demonstrated that reproductive markers, involved in the hypothalamus-hypophysis-gonadal axis, and metabolism markers, flags for hypophysis-melanocortin receptors, are all up-regulated in island lizards under the RIS. This behavior, combined with the observed limited variation in the mitochondrial genome, agrees with the hypothesis that plasticity enables populations to persist in novel environmental conditions and that over time, natural selection will "fine-tune" the population to the environment by modifying the phenotype under selection. We believe that analysis of the transcriptome and the single gene expression, such that all the variations observed in the island populations, can be useful to shed light on evolutionary plasticity as a process affecting animals' populations in general.
Collapse
Affiliation(s)
- Domenico Fulgione
- Department of Biology, University of Naples Federico II, 80126 Naples, Italy; (D.F.); (S.A.); (M.S.); (M.B.)
| | - Valeria Maselli
- Department of Biology, University of Naples Federico II, 80126 Naples, Italy; (D.F.); (S.A.); (M.S.); (M.B.)
| | - Eleonora Rivieccio
- Department of Humanities Studies, University of Naples Federico II, 80138 Naples, Italy;
| | - Serena Aceto
- Department of Biology, University of Naples Federico II, 80126 Naples, Italy; (D.F.); (S.A.); (M.S.); (M.B.)
| | - Marco Salvemini
- Department of Biology, University of Naples Federico II, 80126 Naples, Italy; (D.F.); (S.A.); (M.S.); (M.B.)
| | - Maria Buglione
- Department of Biology, University of Naples Federico II, 80126 Naples, Italy; (D.F.); (S.A.); (M.S.); (M.B.)
| |
Collapse
|
14
|
Ballinger MA, Mack KL, Durkin SM, Riddell EA, Nachman MW. Environmentally robust cis-regulatory changes underlie rapid climatic adaptation. Proc Natl Acad Sci U S A 2023; 120:e2214614120. [PMID: 37725649 PMCID: PMC10523592 DOI: 10.1073/pnas.2214614120] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2022] [Accepted: 07/26/2023] [Indexed: 09/21/2023] Open
Abstract
Changes in gene expression are thought to play a major role in adaptive evolution. While it is known that gene expression is highly sensitive to the environment, very few studies have determined the influence of genetic and environmental effects on adaptive gene expression differences in natural populations. Here, we utilize allele-specific expression to characterize cis and trans gene regulatory divergence in temperate and tropical house mice in two metabolic tissues under two thermal conditions. First, we show that gene expression divergence is pervasive between populations and across thermal conditions, with roughly 5 to 10% of genes exhibiting genotype-by-environment interactions. Second, we found that most expression divergence was due to cis-regulatory changes that were stable across temperatures. In contrast, patterns of expression plasticity were largely attributable to trans-effects, which showed greater sensitivity to temperature. Nonetheless, we found a small subset of temperature-dependent cis-regulatory changes, thereby identifying loci underlying expression plasticity. Finally, we performed scans for selection in wild house mice to identify genomic signatures of rapid adaptation. Genomic outliers were enriched in genes with evidence for cis-regulatory divergence. Notably, these genes were associated with phenotypes that affected body weight and metabolism, suggesting that cis-regulatory changes are a possible mechanism for adaptive body size evolution between populations. Our results show that gene expression plasticity, largely controlled in trans, may facilitate the colonization of new environments, but that evolved changes in gene expression are largely controlled in cis, illustrating the genetic and nongenetic mechanisms underlying the establishment of populations in new environments.
Collapse
Affiliation(s)
- Mallory A. Ballinger
- Museum of Vertebrate Zoology, University of California, Berkeley, CA94720
- Department of Integrative Biology, University of California, Berkeley, CA94720
- Department of Biology, Utah State University, Logan, UT84322
| | - Katya L. Mack
- Department of Biology, Stanford University, Stanford, CA94305
| | - Sylvia M. Durkin
- Museum of Vertebrate Zoology, University of California, Berkeley, CA94720
- Department of Integrative Biology, University of California, Berkeley, CA94720
| | - Eric A. Riddell
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, Ames, IA50011
| | - Michael W. Nachman
- Museum of Vertebrate Zoology, University of California, Berkeley, CA94720
- Department of Integrative Biology, University of California, Berkeley, CA94720
| |
Collapse
|
15
|
Podobnik M, Singh AP, Fu Z, Dooley CM, Frohnhöfer HG, Firlej M, Stednitz SJ, Elhabashy H, Weyand S, Weir JR, Lu J, Nüsslein-Volhard C, Irion U. kcnj13 regulates pigment cell shapes in zebrafish and has diverged by cis-regulatory evolution between Danio species. Development 2023; 150:dev201627. [PMID: 37530080 PMCID: PMC10482006 DOI: 10.1242/dev.201627] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2023] [Accepted: 07/21/2023] [Indexed: 08/03/2023]
Abstract
Teleost fish of the genus Danio are excellent models to study the genetic and cellular bases of pigment pattern variation in vertebrates. The two sister species Danio rerio and Danio aesculapii show divergent patterns of horizontal stripes and vertical bars that are partly caused by the divergence of the potassium channel gene kcnj13. Here, we show that kcnj13 is required only in melanophores for interactions with xanthophores and iridophores, which cause location-specific pigment cell shapes and thereby influence colour pattern and contrast in D. rerio. Cis-regulatory rather than protein coding changes underlie kcnj13 divergence between the two Danio species. Our results suggest that homotypic and heterotypic interactions between the pigment cells and their shapes diverged between species by quantitative changes in kcnj13 expression during pigment pattern diversification.
Collapse
Affiliation(s)
- Marco Podobnik
- Max Planck Institute for Biology, 72076 Tübingen, Germany
| | - Ajeet P. Singh
- Chemical Biology and Therapeutics, Novartis Institutes for BioMedical Research, Cambridge, MA 02139, USA
| | - Zhenqiang Fu
- School of Marine Sciences, Sun Yat-sen University, Zhuhai 519082, China
| | - Christopher M. Dooley
- Department of Genetics, Max Planck Institute for Heart and Lung Research, 61231 Bad Nauheim, Germany
| | | | - Magdalena Firlej
- Friedrich Miescher Laboratory of the Max Planck Society, 72076 Tübingen, Germany
| | - Sarah J. Stednitz
- Department of Anatomy & Physiology, University of Melbourne, Victoria, 3010, Melbourne, Australia
| | - Hadeer Elhabashy
- Department of Protein Evolution, Max Planck Institute for Biology, 72076 Tübingen, Germany
- Institute for Bioinformatics and Medical Informatics, University of Tübingen, 72076 Tübingen, Germany
- Department of Computer Science, University of Tübingen, 72076 Tübingen, Germany
| | - Simone Weyand
- Department of Biochemistry, University of Cambridge, Cambridge, CB2 1QW, UK
| | - John R. Weir
- Friedrich Miescher Laboratory of the Max Planck Society, 72076 Tübingen, Germany
| | - Jianguo Lu
- School of Marine Sciences, Sun Yat-sen University, Zhuhai 519082, China
| | | | - Uwe Irion
- Max Planck Institute for Biology, 72076 Tübingen, Germany
| |
Collapse
|
16
|
Hao Y, Song G, Zhang YE, Zhai W, Jia C, Ji Y, Tang S, Lv H, Qu Y, Lei F. Divergent contributions of coding and noncoding sequences to initial high-altitude adaptation in passerine birds endemic to the Qinghai-Tibet Plateau. Mol Ecol 2023; 32:3524-3540. [PMID: 37000417 DOI: 10.1111/mec.16942] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2022] [Revised: 02/27/2023] [Accepted: 03/27/2023] [Indexed: 04/01/2023]
Abstract
Early events in the evolution of an ancestral lineage can shape the adaptive patterns of descendant species, but the evolutionary mechanisms driving initial adaptation from an ancestor remain largely unexplored. High-altitude adaptations have been extensively explored from the viewpoint of protein-coding genes; however, the contribution of noncoding regions remains relatively neglected. Here, we integrate genomic and transcriptomic data to investigate adaptive evolution in the ancestor of three high-altitude snowfinch species endemic to the Qinghai-Tibet Plateau. Our genome-wide scan for adaptation in the snowfinch ancestor identifies strong adaptation signals in functions of development and metabolism for the coding genes, but in functions of the nervous system development for noncoding regions. This pattern is exclusive to the snowfinch ancestor compared to a control ancestral lineage subject to weak selection. Changes in noncoding regions in the snowfinch ancestor, especially those nearest to coding genes, may be disproportionately associated with the differential expression of genes in the brain tissue compared to other tissues. Extensive gene expression in the brain tissue can be further altered via genetic regulatory networks of transcription factors harbouring potential accelerated regulatory regions (e.g., the development-related transcription factor YEATS4). Altogether, our study provides new evidence concerning how coding and noncoding sequences work through decoupled pathways in initial adaptation to the selective pressure of high-altitude environments. The analysis highlights the idea that noncoding sequences may be promising elements in facilitating the rapid evolution and adaptation to high altitudes.
Collapse
Affiliation(s)
- Yan Hao
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
| | - Gang Song
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
| | - Yong E Zhang
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
- College of Life Science, University of Chinese Academy of Sciences, Beijing, China
- Center for Excellence in Animal Evolution and Genetics, Chinese Academy of Sciences, Kunming, China
| | - Weiwei Zhai
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
- Center for Excellence in Animal Evolution and Genetics, Chinese Academy of Sciences, Kunming, China
| | - Chenxi Jia
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
| | - Yanzhu Ji
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
| | - Shiyu Tang
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
- College of Life Science, University of Chinese Academy of Sciences, Beijing, China
| | - Hongrui Lv
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
- College of Life Science, University of Chinese Academy of Sciences, Beijing, China
| | - Yanhua Qu
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
| | - Fumin Lei
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
- College of Life Science, University of Chinese Academy of Sciences, Beijing, China
- Center for Excellence in Animal Evolution and Genetics, Chinese Academy of Sciences, Kunming, China
| |
Collapse
|
17
|
Cutter AD. Speciation and development. Evol Dev 2023; 25:289-327. [PMID: 37545126 DOI: 10.1111/ede.12454] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2023] [Revised: 06/13/2023] [Accepted: 07/20/2023] [Indexed: 08/08/2023]
Abstract
Understanding general principles about the origin of species remains one of the foundational challenges in evolutionary biology. The genomic divergence between groups of individuals can spawn hybrid inviability and hybrid sterility, which presents a tantalizing developmental problem. Divergent developmental programs may yield either conserved or divergent phenotypes relative to ancestral traits, both of which can be responsible for reproductive isolation during the speciation process. The genetic mechanisms of developmental evolution involve cis- and trans-acting gene regulatory change, protein-protein interactions, genetic network structures, dosage, and epigenetic regulation, all of which also have roots in population genetic and molecular evolutionary processes. Toward the goal of demystifying Darwin's "mystery of mysteries," this review integrates microevolutionary concepts of genetic change with principles of organismal development, establishing explicit links between population genetic process and developmental mechanisms in the production of macroevolutionary pattern. This integration aims to establish a more unified view of speciation that binds process and mechanism.
Collapse
Affiliation(s)
- Asher D Cutter
- Department of Ecology & Evolutionary Biology, University of Toronto, Toronto, Ontario, Canada
| |
Collapse
|
18
|
Mack KL, Square TA, Zhao B, Miller CT, Fraser HB. Evolution of Spatial and Temporal cis-Regulatory Divergence in Sticklebacks. Mol Biol Evol 2023; 40:7048494. [PMID: 36805962 PMCID: PMC10015619 DOI: 10.1093/molbev/msad034] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2022] [Revised: 02/02/2023] [Accepted: 02/08/2023] [Indexed: 02/22/2023] Open
Abstract
Cis-regulatory changes are thought to play a major role in adaptation. Threespine sticklebacks have repeatedly colonized freshwater habitats in the Northern Hemisphere, where they have evolved a suite of phenotypes that distinguish them from marine populations, including changes in physiology, behavior, and morphology. To understand the role of gene regulatory evolution in adaptive divergence, here we investigate cis-regulatory changes in gene expression between marine and freshwater ecotypes through allele-specific expression (ASE) in F1 hybrids. Surveying seven ecologically relevant tissues, including three sampled across two developmental stages, we identified cis-regulatory divergence affecting a third of genes, nearly half of which were tissue-specific. Next, we compared allele-specific expression in dental tissues at two timepoints to characterize cis-regulatory changes during development between marine and freshwater fish. Applying a genome-wide test for selection on cis-regulatory changes, we find evidence for lineage-specific selection on several processes between ecotypes, including the Wnt signaling pathway in dental tissues. Finally, we show that genes with ASE, particularly those that are tissue-specific, are strongly enriched in genomic regions of repeated marine-freshwater divergence, supporting an important role for these cis-regulatory differences in parallel adaptive evolution of sticklebacks to freshwater habitats. Altogether, our results provide insight into the cis-regulatory landscape of divergence between stickleback ecotypes across tissues and during development, and support a fundamental role for tissue-specific cis-regulatory changes in rapid adaptation to new environments.
Collapse
Affiliation(s)
- Katya L Mack
- Department of Biology, Stanford University, Stanford, CA
| | - Tyler A Square
- Department of Molecular and Cell Biology, University of California, Berkeley, CA
| | - Bin Zhao
- Department of Biology, Stanford University, Stanford, CA
| | - Craig T Miller
- Department of Molecular and Cell Biology, University of California, Berkeley, CA
| | | |
Collapse
|
19
|
Fuhrmann N, Prakash C, Kaiser TS. Polygenic adaptation from standing genetic variation allows rapid ecotype formation. eLife 2023; 12:82824. [PMID: 36852484 PMCID: PMC9977305 DOI: 10.7554/elife.82824] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2022] [Accepted: 02/07/2023] [Indexed: 03/01/2023] Open
Abstract
Adaptive ecotype formation can be the first step to speciation, but the genetic underpinnings of this process are poorly understood. Marine midges of the genus Clunio (Diptera) have recolonized Northern European shore areas after the last glaciation. In response to local tide conditions they have formed different ecotypes with respect to timing of adult emergence, oviposition behavior and larval habitat. Genomic analysis confirms the recent establishment of these ecotypes, reflected in massive haplotype sharing between ecotypes, irrespective of whether there is ongoing gene flow or geographic isolation. QTL mapping and genome screens reveal patterns of polygenic adaptation from standing genetic variation. Ecotype-associated loci prominently include circadian clock genes, as well as genes affecting sensory perception and nervous system development, hinting to a central role of these processes in ecotype formation. Our data show that adaptive ecotype formation can occur rapidly, with ongoing gene flow and largely based on a re-assortment of existing alleles.
Collapse
Affiliation(s)
- Nico Fuhrmann
- Max Planck Institute for Evolutionary BiologyPlönGermany
| | - Celine Prakash
- Max Planck Institute for Evolutionary BiologyPlönGermany
| | | |
Collapse
|
20
|
Martínez Sosa F, Pilot M. Molecular Mechanisms Underlying Vertebrate Adaptive Evolution: A Systematic Review. Genes (Basel) 2023; 14:416. [PMID: 36833343 PMCID: PMC9957108 DOI: 10.3390/genes14020416] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2022] [Revised: 01/24/2023] [Accepted: 02/01/2023] [Indexed: 02/08/2023] Open
Abstract
Adaptive evolution is a process in which variation that confers an evolutionary advantage in a specific environmental context arises and is propagated through a population. When investigating this process, researchers have mainly focused on describing advantageous phenotypes or putative advantageous genotypes. A recent increase in molecular data accessibility and technological advances has allowed researchers to go beyond description and to make inferences about the mechanisms underlying adaptive evolution. In this systematic review, we discuss articles from 2016 to 2022 that investigated or reviewed the molecular mechanisms underlying adaptive evolution in vertebrates in response to environmental variation. Regulatory elements within the genome and regulatory proteins involved in either gene expression or cellular pathways have been shown to play key roles in adaptive evolution in response to most of the discussed environmental factors. Gene losses were suggested to be associated with an adaptive response in some contexts. Future adaptive evolution research could benefit from more investigations focused on noncoding regions of the genome, gene regulation mechanisms, and gene losses potentially yielding advantageous phenotypes. Investigating how novel advantageous genotypes are conserved could also contribute to our knowledge of adaptive evolution.
Collapse
Affiliation(s)
| | - Małgorzata Pilot
- Museum and Institute of Zoology, Polish Academy of Sciences, 80-680 Gdańsk, Poland
- Faculty of Biology, University of Gdańsk, 80-308 Gdańsk, Poland
| |
Collapse
|
21
|
Wang C, Li A, Cong R, Qi H, Wang W, Zhang G, Li L. Cis- and Trans-variations of Stearoyl-CoA Desaturase Provide New Insights into the Mechanisms of Diverged Pattern of Phenotypic Plasticity for Temperature Adaptation in Two Congeneric Oyster Species. Mol Biol Evol 2023; 40:6994358. [PMID: 36661848 PMCID: PMC9949715 DOI: 10.1093/molbev/msad015] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2022] [Revised: 12/21/2022] [Accepted: 01/13/2023] [Indexed: 01/21/2023] Open
Abstract
The evolution of phenotypic plasticity plays an essential role in adaptive responses to climate change; however, its regulatory mechanisms in marine organisms which exhibit high phenotypic plasticity still remain poorly understood. The temperature-responsive trait oleic acid content and its major gene stearoyl-CoA desaturase (Scd) expression have diverged in two allopatric congeneric oyster species, cold-adapted Crassostrea gigas and warm-adapted Crassostrea angulata. In this study, genetic and molecular methods were used to characterize fatty acid desaturation and membrane fluidity regulated by oyster Scd. Sixteen causative single-nucleotide polymorphisms (SNPs) were identified in the promoter/cis-region of the Scd between wild C. gigas and C. angulata. Further functional experiments showed that an SNP (g.-333C [C. gigas allele] >T [C. angulata allele]) may influence Scd transcription by creating/disrupting the binding motif of the positive trans-factor Y-box factor in C. gigas/C. angulata, which mediates the higher/lower constitutive expression of Scd in C. gigas/C. angulata. Additionally, the positive trans-factor sterol-regulatory element-binding proteins (Srebp) were identified to specifically bind to the promoter of Scd in both species, and were downregulated during cold stress in C. gigas compared to upregulated in C. angulata. This partly explains the relatively lower environmental sensitivity (plasticity) of Scd in C. gigas. This study serves as an experimental case to reveal that both cis- and trans-variations shape the diverged pattern of phenotypic plasticity, which provides new insights into the formation of adaptive traits and the prediction of the adaptive potential of marine organisms to future climate change.
Collapse
Affiliation(s)
- Chaogang Wang
- CAS and Shandong Province Key Laboratory of Experimental Marine Biology, Center for Ocean Mega-Science, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China,Laboratory for Marine Biology and Biotechnology, Pilot National Laboratory for Marine Science and Technology, Qingdao, China,University of Chinese Academy of Sciences, Beijing, China
| | - Ao Li
- CAS and Shandong Province Key Laboratory of Experimental Marine Biology, Center for Ocean Mega-Science, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China,Laboratory for Marine Fisheries Science and Food Production Processes, Pilot National Laboratory for Marine Science and Technology, Qingdao, China,National and Local Joint Engineering Laboratory of Ecological Mariculture, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China
| | - Rihao Cong
- CAS and Shandong Province Key Laboratory of Experimental Marine Biology, Center for Ocean Mega-Science, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China,Laboratory for Marine Biology and Biotechnology, Pilot National Laboratory for Marine Science and Technology, Qingdao, China,National and Local Joint Engineering Laboratory of Ecological Mariculture, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China
| | - Haigang Qi
- CAS and Shandong Province Key Laboratory of Experimental Marine Biology, Center for Ocean Mega-Science, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China,Laboratory for Marine Biology and Biotechnology, Pilot National Laboratory for Marine Science and Technology, Qingdao, China,National and Local Joint Engineering Laboratory of Ecological Mariculture, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China
| | - Wei Wang
- CAS and Shandong Province Key Laboratory of Experimental Marine Biology, Center for Ocean Mega-Science, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China,Laboratory for Marine Biology and Biotechnology, Pilot National Laboratory for Marine Science and Technology, Qingdao, China,National and Local Joint Engineering Laboratory of Ecological Mariculture, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China
| | - Guofan Zhang
- CAS and Shandong Province Key Laboratory of Experimental Marine Biology, Center for Ocean Mega-Science, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China,University of Chinese Academy of Sciences, Beijing, China,Laboratory for Marine Fisheries Science and Food Production Processes, Pilot National Laboratory for Marine Science and Technology, Qingdao, China,National and Local Joint Engineering Laboratory of Ecological Mariculture, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China
| | - Li Li
- Corresponding author: E-mail:
| |
Collapse
|
22
|
Hu J, Barrett RDH. The role of plastic and evolved DNA methylation in parallel adaptation of threespine stickleback (Gasterosteus aculeatus). Mol Ecol 2022; 32:1581-1591. [PMID: 36560898 DOI: 10.1111/mec.16832] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2022] [Revised: 12/17/2022] [Accepted: 12/20/2022] [Indexed: 12/24/2022]
Abstract
Repeated phenotypic patterns among populations undergoing parallel evolution in similar environments provide support for the deterministic role of natural selection. Epigenetic modifications can mediate plastic and evolved phenotypic responses to environmental change and might make important contributions to parallel adaptation. While many studies have explored the genetic basis of repeated phenotypic divergence, the role of epigenetic processes during parallel adaptation remains unclear. The parallel evolution of freshwater ecotypes of threespine stickleback fish (Gasterosteus aculeatus) following colonization of thousands of lakes and streams from the ocean is a classic example of parallel phenotypic and genotypic adaptation. To investigate epigenetic modifications during parallel adaptation of threespine stickleback, we reanalysed three independent data sets that investigated DNA methylation variation between marine and freshwater ecotypes. Although we found widespread methylation differentiation between ecotypes, there was no significant tendency for CpG sites associated with repeated methylation differentiation across studies to be parallel versus nonparallel. To next investigate the role of plastic versus evolved changes in methylation during freshwater adaptation, we explored if CpG sites exhibiting methylation plasticity during salinity change were more likely to also show evolutionary divergence in methylation between ecotypes. The directions of divergence between ecotypes were generally in the opposite direction to those observed for plasticity when ecotypes were challenged with non-native salinity conditions, suggesting that most plastic responses are likely to be maladaptive during colonization of new environments. Finally, we found a greater number of CpG sites showing evolved changes when ancestral marine ecotypes are acclimated to freshwater environments, whereas plastic changes predominate when derived freshwater ecotypes transition back to their ancestral marine environments. These findings provide evidence for an epigenetic contribution to parallel adaptation and demonstrate the contrasting roles of plastic and evolved methylation differences during adaptation to new environments.
Collapse
Affiliation(s)
- Juntao Hu
- Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, Institute of Biodiversity Science, Center of Evolutionary Biology, School of Life Sciences, Fudan University, Shanghai, P. R. China
| | - Rowan D H Barrett
- Redpath Museum and Department of Biology, McGill University, Montreal, Quebec, Canada
| |
Collapse
|
23
|
Trans-driven variation in expression is common among detoxification genes in the extreme generalist herbivore Tetranychus urticae. PLoS Genet 2022; 18:e1010333. [DOI: 10.1371/journal.pgen.1010333] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2022] [Revised: 11/28/2022] [Accepted: 10/31/2022] [Indexed: 11/15/2022] Open
Abstract
The extreme adaptation potential of the generalist herbivore Tetranychus urticae (the two-spotted spider mite) to pesticides as well as diverse host plants has been associated with clade-specific gene expansions in known detoxifying enzyme families, and with extensive and rapid transcriptional responses. However, how this broad transcriptional potential is regulated remains largely unknown. Using a parental/F1 design in which four inbred strains were crossed to a common inbred strain, we assessed the genetic basis and inheritance of gene expression variation in T. urticae. Mirroring known phenotypic variation in the progenitor strains of the inbreds, we confirmed that the inbred strains we created were genetically distinct, varied markedly in pesticide resistance, and also captured variation in host plant fitness as is commonly observed in this species. By examining differences in gene expression between parents and allele-specific expression in F1s, we found that variation in RNA abundance was more often explained in trans as compared to cis, with the former associated with dominance in inheritance. Strikingly, in a gene ontology analysis, detoxification genes of the cytochrome P450 monooxygenase (CYP) family, as well as dioxygenases (DOGs) acquired from horizontal gene transfer from fungi, were specifically enriched at the extremes of trans-driven up- and downregulation. In particular, multiple CYPs and DOGs with broad substrate-specificities for pesticides or plant specialized compounds were exceptionally highly upregulated as a result of trans-regulatory variation, or in some cases synergism of cis and trans, in the most multi-pesticide resistant strains. Collectively, our findings highlight the potential importance of trans-driven expression variation in genes associated with xenobiotic metabolism and host plant use for rapid adaptation in T. urticae, and also suggests modular control of these genes, a regulatory architecture that might ameliorate negative pleiotropic effects.
Collapse
|
24
|
Carruthers M, Edgley DE, Saxon AD, Gabagambi NP, Shechonge A, Miska EA, Durbin R, Bridle JR, Turner GF, Genner MJ. Ecological Speciation Promoted by Divergent Regulation of Functional Genes Within African Cichlid Fishes. Mol Biol Evol 2022; 39:msac251. [PMID: 36376993 PMCID: PMC10101686 DOI: 10.1093/molbev/msac251] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022] Open
Abstract
Rapid ecological speciation along depth gradients has taken place repeatedly in freshwater fishes, yet molecular mechanisms facilitating such diversification are typically unclear. In Lake Masoko, an African crater lake, the cichlid Astatotilapia calliptera has diverged into shallow-littoral and deep-benthic ecomorphs with strikingly different jaw structures within the last 1,000 years. Using genome-wide transcriptome data, we explore two major regulatory transcriptional mechanisms, expression and splicing-QTL variants, and examine their contributions to differential gene expression underpinning functional phenotypes. We identified 7,550 genes with significant differential expression between ecomorphs, of which 5.4% were regulated by cis-regulatory expression QTLs, and 9.2% were regulated by cis-regulatory splicing QTLs. We also found strong signals of divergent selection on differentially expressed genes associated with craniofacial development. These results suggest that large-scale transcriptome modification plays an important role during early-stage speciation. We conclude that regulatory variants are important targets of selection driving ecologically relevant divergence in gene expression during adaptive diversification.
Collapse
Affiliation(s)
- Madeleine Carruthers
- School of Biological Sciences, University of Bristol,
Bristol BS8 1TQ, United
Kingdom
| | - Duncan E Edgley
- School of Biological Sciences, University of Bristol,
Bristol BS8 1TQ, United
Kingdom
| | - Andrew D Saxon
- School of Biological Sciences, University of Bristol,
Bristol BS8 1TQ, United
Kingdom
| | - Nestory P Gabagambi
- Tanzanian Fisheries Research Institute, Kyela Research
Centre, P.O. Box 98, Kyela, Mbeya, Tanzania
| | - Asilatu Shechonge
- Tanzanian Fisheries Research Institute, Dar es Salaam Research
Centre, P.O. Box 9750, Dar es Salaam, Tanzania
| | - Eric A Miska
- Wellcome/CRUK Gurdon Institute, University of Cambridge,
Cambridge CB2 1QN, United
Kingdom
- Department of Genetics, University of Cambridge,
Cambridge CB2 3EH, United
Kingdom
- Wellcome Sanger Institute, Wellcome Genome Campus,
Cambridge CB10 1SA, United Kingdom
| | - Richard Durbin
- Department of Genetics, University of Cambridge,
Cambridge CB2 3EH, United
Kingdom
- Wellcome Sanger Institute, Wellcome Genome Campus,
Cambridge CB10 1SA, United Kingdom
| | - Jon R Bridle
- School of Biological Sciences, University of Bristol,
Bristol BS8 1TQ, United
Kingdom
| | - George F Turner
- School of Natural Sciences, Bangor University,
Bangor, Wales LL57 2UW, United
Kingdom
| | - Martin J Genner
- School of Biological Sciences, University of Bristol,
Bristol BS8 1TQ, United
Kingdom
| |
Collapse
|
25
|
Ereful NC, Lalusin AG, Laurena AC. Assessing Loss of Regulatory Divergence, Genome-Transcriptome Incongruence, and Preferential Expression Switching in Abaca × Banana Backcrosses. Genes (Basel) 2022; 13:genes13081396. [PMID: 36011307 PMCID: PMC9407414 DOI: 10.3390/genes13081396] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2022] [Revised: 08/02/2022] [Accepted: 08/03/2022] [Indexed: 11/16/2022] Open
Abstract
The Musa textilis var. Abuab has high fiber quality (FQ) but is susceptible to abaca bunchy top virus (AbBTV); the Musa balbisiana var. Pacol has low FQ but is resistant against AbBTV. Their backcrosses (BC2 and BC3) possess both desirable traits. Analysis using RNA-seq showed that the regulatory divergence of Abuab and Pacol is largely explained by cis differences with 27.4% and 22.3% if we are to assess it using BC2 and BC3, respectively. Cis differences between the two genotypes are significantly reduced from BC2 to BC3 due to changes in genomic constitution. Trans, on the other hand, is robust to changes in allelic composition. All these are attributed to the loss of heterozygosity in BC3 relative to BC2. Further analysis showed that both backcrosses exhibited genome-wide preferential expression of Pacol- over Abuab-specific alleles, despite the wider genetic presence of the latter in the hybrids. The ratio of the two genotype-specific expressed transcripts and the ratio of their corresponding genetic make-up are significantly disproportionate, a phenomenon that we refer to here as "genome-transcriptome incongruence". We also observed preferential expression switching in which several genes prefer the Abuab- (or Pacol-) specific allele in BC2 but switched to the Pacol- (or Abuab-) specific allele in the BC3 genome.
Collapse
Affiliation(s)
- Nelzo C. Ereful
- Biochemistry Laboratory–Plant Physiology Laboratory, Institute of Plant Breeding, College of Agriculture and Food Science, University of the Philippines Los Baños, Laguna 4031, Philippines
- Philippine Genome Center for Agriculture, University of the Philippines Los Baños, Laguna 4031, Philippines
- Correspondence:
| | - Antonio G. Lalusin
- Institute of Crop Science, College of Agriculture and Food Science, University of the Philippines Los Baños, Laguna 4031, Philippines
| | - Antonio C. Laurena
- Philippine Genome Center for Agriculture, University of the Philippines Los Baños, Laguna 4031, Philippines
| |
Collapse
|
26
|
Ji M, Bui H, Ramirez RA, Clark RM. Concerted cis and trans effects underpin heightened defense gene expression in multi-herbivore-resistant maize lines. THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2022; 111:508-528. [PMID: 35575017 DOI: 10.1111/tpj.15812] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/19/2021] [Revised: 04/04/2022] [Accepted: 05/09/2022] [Indexed: 06/15/2023]
Abstract
In maize (Zea mays ssp. mays), agriculturally damaging herbivores include lepidopteran insects, such as the European corn borer (Ostrinia nubilalis), and distantly related arthropods, like the two-spotted spider mite (Tetranychus urticae). A small number of maize lines, including B96 and B75, are highly resistant to both herbivores, and B96 is also resistant to thrips. Using T. urticae as a representative pest that causes significant leaf tissue damage, we examined the gene expression responses of these lines to herbivory in comparison with each other and with the susceptible line B73. Upon herbivory, the most resistant line, B96, showed the strongest gene expression response, with a dramatic upregulation of genes associated with jasmonic acid biosynthesis and signaling, as well as the biosynthesis of specialized herbivore deterrent compounds, such as death acids and benzoxazinoids. Extending this work with allele-specific expression analyses in F1 hybrids, we inferred that the concerted upregulation of many defense genes, including the majority of benzoxazinoid biosynthetic genes in B96, as compared with B73, for the herbivore treatment, resulted from an assemblage of trans control and multiple cis effects acting with similar directionality on gene expression. Further, at the level of individual and potentially rate limiting genes in several major defense pathways, cis and trans effects acted in a reinforcing manner to result in exceptionally high expression in B96. Our study provides a comprehensive resource of cis elements for maize lines important in breeding efforts for herbivore resistance, and reveals potential genetic underpinnings of the origins of multi-herbivore resistance in plant populations.
Collapse
Affiliation(s)
- Meiyuan Ji
- School of Biological Sciences, University of Utah, 257 South 1400 East Rm 201, Salt Lake City, UT, 84112, USA
| | - Huyen Bui
- School of Biological Sciences, University of Utah, 257 South 1400 East Rm 201, Salt Lake City, UT, 84112, USA
| | - Ricardo A Ramirez
- Department of Biology, Utah State University, 5305 Old Main Hill, Logan, UT, 84332, USA
| | - Richard M Clark
- School of Biological Sciences, University of Utah, 257 South 1400 East Rm 201, Salt Lake City, UT, 84112, USA
- Henry Eyring Center for Cell and Genome Science, University of Utah, 1390 Presidents Circle, Salt Lake City, UT, 84112, USA
| |
Collapse
|
27
|
Huang Y, Shang R, Lu GA, Zeng W, Huang C, Zou C, Tang T. Spatiotemporal Regulation of a Single Adaptively Evolving Trans-Regulatory Element Contributes to Spermatogenetic Expression Divergence in Drosophila. Mol Biol Evol 2022; 39:6605656. [PMID: 35687719 PMCID: PMC9254010 DOI: 10.1093/molbev/msac127] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022] Open
Abstract
Due to extensive pleiotropy, trans-acting elements are often thought to be evolutionarily constrained. While the impact of trans-acting elements on gene expression evolution has been extensively studied, relatively little is understood about the contribution of a single trans regulator to interspecific expression and phenotypic divergence. Here, we disentangle the effects of genomic context and miR-983, an adaptively evolving young microRNA, on expression divergence between Drosophila melanogaster and D. simulans. We show miR-983 effects promote interspecific expression divergence in testis despite its antagonism with the often-predominant context effects. Single-cyst RNA-seq reveals that distinct sets of genes gain and lose miR-983 influence under disruptive or diversifying selection at different stages of spermatogenesis, potentially helping minimize antagonistic pleiotropy. At the round spermatid stage, the effects of miR-983 are weak and distributed, coincident with the transcriptome undergoing drastic expression changes. Knocking out miR-983 causes reduced sperm length with increased within-individual variation in D. melanogaster but not in D. simulans, and the D. melanogaster knockout also exhibits compromised sperm defense ability. Our results provide empirical evidence for the resolution of antagonistic pleiotropy and also have broad implications for the function and evolution of new trans regulators.
Collapse
Affiliation(s)
- Yumei Huang
- State Key Laboratory of Biocontrol, School of Life Sciences, Sun Yat-sen University, 510275 Guangzhou, Guangdong Province, China
| | - Rui Shang
- State Key Laboratory of Biocontrol, School of Life Sciences, Sun Yat-sen University, 510275 Guangzhou, Guangdong Province, China
| | - Guang-An Lu
- State Key Laboratory of Biocontrol, School of Life Sciences, Sun Yat-sen University, 510275 Guangzhou, Guangdong Province, China
| | - Weishun Zeng
- State Key Laboratory of Biocontrol, School of Life Sciences, Sun Yat-sen University, 510275 Guangzhou, Guangdong Province, China
| | - Chenglong Huang
- State Key Laboratory of Biocontrol, School of Life Sciences, Sun Yat-sen University, 510275 Guangzhou, Guangdong Province, China
| | - Chuangchao Zou
- State Key Laboratory of Biocontrol, School of Life Sciences, Sun Yat-sen University, 510275 Guangzhou, Guangdong Province, China
| | - Tian Tang
- State Key Laboratory of Biocontrol, School of Life Sciences, Sun Yat-sen University, 510275 Guangzhou, Guangdong Province, China
| |
Collapse
|
28
|
Li Q, Lindtke D, Rodríguez-Ramírez C, Kakioka R, Takahashi H, Toyoda A, Kitano J, Ehrlich RL, Chang Mell J, Yeaman S. Local Adaptation and the Evolution of Genome Architecture in Threespine Stickleback. Genome Biol Evol 2022; 14:6589818. [PMID: 35594844 PMCID: PMC9178229 DOI: 10.1093/gbe/evac075] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/16/2022] [Indexed: 12/11/2022] Open
Abstract
Theory predicts that local adaptation should favor the evolution of a concentrated genetic architecture, where the alleles driving adaptive divergence are tightly clustered on chromosomes. Adaptation to marine versus freshwater environments in threespine stickleback has resulted in an architecture that seems consistent with this prediction: divergence among populations is mainly driven by a few genomic regions harboring multiple quantitative trait loci for environmentally adapted traits, as well as candidate genes with well-established phenotypic effects. One theory for the evolution of these "genomic islands" is that rearrangements remodel the genome to bring causal loci into tight proximity, but this has not been studied explicitly. We tested this theory using synteny analysis to identify micro- and macro-rearrangements in the stickleback genome and assess their potential involvement in the evolution of genomic islands. To identify rearrangements, we conducted a de novo assembly of the closely related tubesnout (Aulorhyncus flavidus) genome and compared this to the genomes of threespine stickleback and two other closely related species. We found that small rearrangements, within-chromosome duplications, and lineage-specific genes (LSGs) were enriched around genomic islands, and that all three chromosomes harboring large genomic islands have experienced macro-rearrangements. We also found that duplicates and micro-rearrangements are 9.9× and 2.9× more likely to involve genes differentially expressed between marine and freshwater genotypes. While not conclusive, these results are consistent with the explanation that strong divergent selection on candidate genes drove the recruitment of rearrangements to yield clusters of locally adaptive loci.
Collapse
Affiliation(s)
- Qiushi Li
- Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Canada T2N 1N4
| | - Dorothea Lindtke
- Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Canada T2N 1N4
| | - Carlos Rodríguez-Ramírez
- Division of Evolutionary Ecology, Institute of Ecology and Evolution, University of Bern, Bern, Switzerland
| | - Ryo Kakioka
- Tropical Biosphere Research Center, University of the Ryukyus, Nishihara, Nakagami-gun, Okinawa 903-0213, Japan
| | - Hiroshi Takahashi
- National Fisheries University, 2-7-1 Nagata-honmachi, Shimonoseki, Yamaguchi 759-6595, Japan
| | - Atsushi Toyoda
- Comparative Genomics Laboratory, National Institute of Genetics, Mishima, Shizuoka 411-8540, Japan
| | - Jun Kitano
- Ecological Genetics Laboratory, National Institute of Genetics, Mishima, Shizuoka 411-8540, Japan
| | - Rachel L Ehrlich
- Department of Microbiology & Immunology, Drexel University College of Medicine, Philadelphia 19102, PA, USA
| | - Joshua Chang Mell
- Department of Microbiology & Immunology, Drexel University College of Medicine, Philadelphia 19102, PA, USA
| | - Sam Yeaman
- Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Canada T2N 1N4
| |
Collapse
|
29
|
Kim C, Wang X, Kültz D. Prediction and Experimental Validation of a New Salinity-Responsive Cis-Regulatory Element (CRE) in a Tilapia Cell Line. Life (Basel) 2022; 12:life12060787. [PMID: 35743818 PMCID: PMC9225295 DOI: 10.3390/life12060787] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2022] [Revised: 05/12/2022] [Accepted: 05/16/2022] [Indexed: 11/16/2022] Open
Abstract
Transcriptional regulation is a major mechanism by which organisms integrate gene x environment interactions. It can be achieved by coordinated interplay between cis-regulatory elements (CREs) and transcription factors (TFs). Euryhaline tilapia (Oreochromis mossambicus) tolerate a wide range of salinity and thus are an appropriate model to examine transcriptional regulatory mechanisms during salinity stress in fish. Quantitative proteomics in combination with the transcription inhibitor actinomycin D revealed 19 proteins that are transcriptionally upregulated by hyperosmolality in tilapia brain (OmB) cells. We searched the extended proximal promoter up to intron1 of each corresponding gene for common motifs using motif discovery tools. The top-ranked motif identified (STREME1) represents a binding site for the Forkhead box TF L1 (FoxL1). STREME1 function during hyperosmolality was experimentally validated by choosing two of the 19 genes, chloride intracellular channel 2 (clic2) and uridine phosphorylase 1 (upp1), that are enriched in STREME1 in their extended promoters. Transcriptional induction of these genes during hyperosmolality requires STREME1, as evidenced by motif mutagenesis. We conclude that STREME1 represents a new functional CRE that contributes to gene x environment interactions during salinity stress in tilapia. Moreover, our results indicate that FoxL1 family TFs are contribute to hyperosmotic induction of genes in euryhaline fish.
Collapse
Affiliation(s)
- Chanhee Kim
- Stress-Induced Evolution Laboratory, Department of Animal Sciences, University of California, Davis, CA 95616, USA;
| | - Xiaodan Wang
- Laboratory of Aquaculture Nutrition and Environmental Health, School of Life Sciences, East China Normal University, Shanghai 200241, China;
| | - Dietmar Kültz
- Stress-Induced Evolution Laboratory, Department of Animal Sciences, University of California, Davis, CA 95616, USA;
- Correspondence: ; Tel.: +1-530-752-2991
| |
Collapse
|
30
|
Ozerov M, Noreikiene K, Kahar S, Huss M, Huusko A, Kõiv T, Sepp M, López M, Gårdmark A, Gross R, Vasemägi A. Whole-genome sequencing illuminates multifaceted targets of selection to humic substances in Eurasian perch. Mol Ecol 2022; 31:2367-2383. [PMID: 35202502 PMCID: PMC9314028 DOI: 10.1111/mec.16409] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2021] [Revised: 02/10/2022] [Accepted: 02/17/2022] [Indexed: 11/30/2022]
Abstract
Extreme environments are inhospitable to the majority of species, but some organisms are able to survive in such hostile conditions due to evolutionary adaptations. For example, modern bony fishes have colonized various aquatic environments, including perpetually dark, hypoxic, hypersaline and toxic habitats. Eurasian perch (Perca fluviatilis) is among the few fish species of northern latitudes that is able to live in very acidic humic lakes. Such lakes represent almost "nocturnal" environments; they contain high levels of dissolved organic matter, which in addition to creating a challenging visual environment, also affects a large number of other habitat parameters and biotic interactions. To reveal the genomic targets of humic-associated selection, we performed whole-genome sequencing of perch originating from 16 humic and 16 clear-water lakes in northern Europe. We identified over 800,000 SNPs, of which >10,000 were identified as potential candidates under selection (associated with >3,000 genes) using multiple outlier approaches. Our findings suggest that adaptation to the humic environment may involve hundreds of regions scattered across the genome. Putative signals of adaptation were detected in genes and gene families with diverse functions, including organism development and ion transportation. The observed excess of variants under selection in regulatory regions highlights the importance of adaptive evolution via regulatory elements, rather than via protein sequence modification. Our study demonstrates the power of whole-genome analysis to illuminate multifaceted nature of humic adaptation and provides the foundation for further investigation of causal mutations underlying phenotypic traits of ecological and evolutionary importance.
Collapse
Affiliation(s)
- Mikhail Ozerov
- Department of Aquatic Resources, Institute of Freshwater Research, Swedish University of Agricultural Sciences, 17893, Drottningholm, Sweden.,Department of Biology, University of Turku, 20014, Turku, Finland.,Biodiversity Unit, University of Turku, 20014, Turku, Finland
| | - Kristina Noreikiene
- Chair of Aquaculture, Institute of Veterinary Medicine and Animal Sciences, Estonian University of Life Sciences, Kreutzwaldi 46, 51006, Tartu, Estonia
| | - Siim Kahar
- Chair of Aquaculture, Institute of Veterinary Medicine and Animal Sciences, Estonian University of Life Sciences, Kreutzwaldi 46, 51006, Tartu, Estonia
| | - Magnus Huss
- Swedish University of Agricultural Sciences, Department of Aquatic Resources, 74242, Öregrund, Sweden
| | - Ari Huusko
- Natural resources Institute Finland (Luke), 88300, Paltamo, Finland
| | - Toomas Kõiv
- Chair of Hydrobiology and Fishery, Institute of Agricultural and Environmental Sciences, Estonian University of Life Sciences, Kreutzwaldi 5, 51006, Tartu, Estonia
| | - Margot Sepp
- Chair of Hydrobiology and Fishery, Institute of Agricultural and Environmental Sciences, Estonian University of Life Sciences, Kreutzwaldi 5, 51006, Tartu, Estonia
| | - María López
- Department of Aquatic Resources, Institute of Freshwater Research, Swedish University of Agricultural Sciences, 17893, Drottningholm, Sweden
| | - Anna Gårdmark
- Swedish University of Agricultural Sciences, Department of Aquatic Resources, 74242, Öregrund, Sweden
| | - Riho Gross
- Chair of Aquaculture, Institute of Veterinary Medicine and Animal Sciences, Estonian University of Life Sciences, Kreutzwaldi 46, 51006, Tartu, Estonia
| | - Anti Vasemägi
- Department of Aquatic Resources, Institute of Freshwater Research, Swedish University of Agricultural Sciences, 17893, Drottningholm, Sweden.,Chair of Aquaculture, Institute of Veterinary Medicine and Animal Sciences, Estonian University of Life Sciences, Kreutzwaldi 46, 51006, Tartu, Estonia
| |
Collapse
|
31
|
cis-Regulatory changes in locomotor genes are associated with the evolution of burrowing behavior. Cell Rep 2022; 38:110360. [PMID: 35172153 DOI: 10.1016/j.celrep.2022.110360] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2021] [Revised: 11/24/2021] [Accepted: 01/19/2022] [Indexed: 12/13/2022] Open
Abstract
How evolution modifies complex, innate behaviors is largely unknown. Divergence in many morphological traits, and some behaviors, is linked to cis-regulatory changes in gene expression. Given this, we compare brain gene expression of two interfertile sister species of Peromyscus mice that show large and heritable differences in burrowing behavior. Species-level differential expression and allele-specific expression in F1 hybrids indicate a preponderance of cis-regulatory divergence, including many genes whose cis-regulation is affected by burrowing behavior. Genes related to locomotor coordination show the strongest signals of lineage-specific selection on burrowing-induced cis-regulatory changes. Furthermore, genetic markers closest to these candidate genes associate with variation in burrow shape in a genetic cross, suggesting an enrichment for loci affecting burrowing behavior near these candidate locomotor genes. Our results provide insight into how cis-regulated gene expression can depend on behavioral context and how this dynamic regulatory divergence between species may contribute to behavioral evolution.
Collapse
|
32
|
Pélissié B, Chen YH, Cohen ZP, Crossley MS, Hawthorne DJ, Izzo V, Schoville SD. Genome resequencing reveals rapid, repeated evolution in the Colorado potato beetle. Mol Biol Evol 2022; 39:6511499. [PMID: 35044459 PMCID: PMC8826761 DOI: 10.1093/molbev/msac016] [Citation(s) in RCA: 27] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Insecticide resistance and rapid pest evolution threatens food security and the development of sustainable agricultural practices, yet the evolutionary mechanisms that allow pests to rapidly adapt to control tactics remains unclear. Here we examine how a global super-pest, the Colorado potato beetle (CPB), Leptinotarsa decemlineata, rapidly evolves resistance to insecticides. Using whole genome resequencing and transcriptomic data focused on its ancestral and pest range in North America, we assess evidence for three, non-mutually exclusive models of rapid evolution: pervasive selection on novel mutations, rapid regulatory evolution, and repeated selection on standing genetic variation. Population genomic analysis demonstrates that CPB is geographically structured, even among recently established pest populations. Pest populations exhibit similar levels of nucleotide diversity, relative to non-pest populations, and show evidence of recent expansion. Genome scans provide clear signatures of repeated adaptation across CPB populations, with especially strong evidence of selection on insecticide resistance genes in different populations. Analyses of gene expression show that constitutive upregulation of candidate insecticide resistance genes drives distinctive population patterns. CPB evolves insecticide resistance repeatedly across agricultural regions, leveraging similar genetic pathways but different genes, demonstrating a polygenic trait architecture for insecticide resistance that can evolve from standing genetic variation. Despite expectations, we do not find support for strong selection on novel mutations, or rapid evolution from selection on regulatory genes. These results suggest that integrated pest management practices must mitigate the evolution of polygenic resistance phenotypes among local pest populations, in order to maintain the efficacy and sustainability of novel control techniques.
Collapse
Affiliation(s)
- Benjamin Pélissié
- Department of Entomology, University of Wisconsin-Madison, Madison, WI 53706, USA
| | - Yolanda H Chen
- Department of Plant and Soil Science, University of Vermont, Burlington, VT 05405, USA
| | - Zachary P Cohen
- Department of Entomology, University of Wisconsin-Madison, Madison, WI 53706, USA
| | - Michael S Crossley
- Department of Entomology, University of Wisconsin-Madison, Madison, WI 53706, USA
| | - David J Hawthorne
- Department of Entomology, University of Maryland, College Park, MD 20742, USA
| | - Victor Izzo
- Department of Plant and Soil Science, University of Vermont, Burlington, VT 05405, USA
| | - Sean D Schoville
- Department of Entomology, University of Wisconsin-Madison, Madison, WI 53706, USA
| |
Collapse
|
33
|
Verta JP, Jacobs A. The role of alternative splicing in adaptation and evolution. Trends Ecol Evol 2021; 37:299-308. [PMID: 34920907 DOI: 10.1016/j.tree.2021.11.010] [Citation(s) in RCA: 38] [Impact Index Per Article: 12.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2021] [Revised: 10/26/2021] [Accepted: 11/19/2021] [Indexed: 01/02/2023]
Abstract
Regulation of gene expression plays a central role in adaptive divergence and evolution. Although the role of gene regulation in microevolutionary processes is gaining wide acceptance, most studies have only investigated the evolution of transcript levels, ignoring the potentially significant role of transcript structures. We argue that variation in alternative splicing plays an important and widely unexplored role in adaptation (e.g., by increasing transcriptome and/or proteome diversity, or buffering potentially deleterious genetic variation). New studies increasingly highlight the potential for independent evolution in alternative splicing and transcript level, providing alternative paths for selection to act upon. We propose that alternative splicing and transcript levels can provide contrasting, nonredundant mechanisms of equal importance for adaptive diversification of gene function and regulation.
Collapse
Affiliation(s)
- Jukka-Pekka Verta
- Organismal and Evolutionary Biology Research Programme, University of Helsinki, Viikinkaari 9, 00790, Helsinki, Finland.
| | - Arne Jacobs
- Institute of Biodiversity, Animal Health, and Comparative Medicine, University of Glasgow, G12 8QQ, Glasgow, UK.
| |
Collapse
|
34
|
James ME, Wilkinson MJ, Bernal DM, Liu H, North HL, Engelstädter J, Ortiz-Barrientos D. Phenotypic and genotypic parallel evolution in parapatric ecotypes of Senecio. Evolution 2021; 75:3115-3131. [PMID: 34687472 PMCID: PMC9299460 DOI: 10.1111/evo.14387] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2021] [Revised: 10/08/2021] [Accepted: 10/12/2021] [Indexed: 12/11/2022]
Abstract
The independent and repeated adaptation of populations to similar environments often results in the evolution of similar forms. This phenomenon creates a strong correlation between phenotype and environment and is referred to as parallel evolution. However, we are still largely unaware of the dynamics of parallel evolution, as well as the interplay between phenotype and genotype within natural systems. Here, we examined phenotypic and genotypic parallel evolution in multiple parapatric Dune‐Headland coastal ecotypes of an Australian wildflower, Senecio lautus. We observed a clear trait‐environment association in the system, with all replicate populations having evolved along the same phenotypic evolutionary trajectory. Similar phenotypes have arisen via mutational changes occurring in different genes, although many share the same biological functions. Our results shed light on how replicated adaptation manifests at the phenotypic and genotypic levels within populations, and highlight S. lautus as one of the most striking cases of phenotypic parallel evolution in nature.
Collapse
Affiliation(s)
- Maddie E James
- School of Biological Sciences, The University of Queensland, St. Lucia, QLD, 4072, Australia.,Australian Research Council Centre of Excellence for Plant Success in Nature and Agriculture, The University of Queensland, St Lucia, QLD, 4072, Australia
| | - Melanie J Wilkinson
- School of Biological Sciences, The University of Queensland, St. Lucia, QLD, 4072, Australia.,Australian Research Council Centre of Excellence for Plant Success in Nature and Agriculture, The University of Queensland, St Lucia, QLD, 4072, Australia
| | - Diana M Bernal
- School of Biological Sciences, The University of Queensland, St. Lucia, QLD, 4072, Australia.,Current Address: Biousos Neotropicales S.A.S, Bogotá, Colombia
| | - Huanle Liu
- School of Biological Sciences, The University of Queensland, St. Lucia, QLD, 4072, Australia.,Current Address: Centre for Genomic Regulation, The Barcelona Institute of Science and Technology, Barcelona, 08003, Spain
| | - Henry L North
- School of Biological Sciences, The University of Queensland, St. Lucia, QLD, 4072, Australia.,Current Address: Department of Zoology, University of Cambridge, Cambridge, CB2 3EJ, United Kingdom
| | - Jan Engelstädter
- School of Biological Sciences, The University of Queensland, St. Lucia, QLD, 4072, Australia
| | - Daniel Ortiz-Barrientos
- School of Biological Sciences, The University of Queensland, St. Lucia, QLD, 4072, Australia.,Australian Research Council Centre of Excellence for Plant Success in Nature and Agriculture, The University of Queensland, St Lucia, QLD, 4072, Australia
| |
Collapse
|
35
|
Jacobs A, Elmer KR. Alternative splicing and gene expression play contrasting roles in the parallel phenotypic evolution of a salmonid fish. Mol Ecol 2021; 30:4955-4969. [PMID: 33502030 PMCID: PMC8653899 DOI: 10.1111/mec.15817] [Citation(s) in RCA: 36] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2020] [Revised: 01/06/2021] [Accepted: 01/18/2021] [Indexed: 12/25/2022]
Abstract
Understanding the contribution of different molecular processes to evolution and development is crucial for identifying the mechanisms of adaptation. Here, we used RNA-sequencing data to test the importance of alternative splicing and differential gene expression in a case of parallel adaptive evolution, the replicated postglacial divergence of the salmonid fish Arctic charr (Salvelinus alpinus) into sympatric benthic and pelagic ecotypes across multiple independent lakes. We found that genes differentially spliced between ecotypes were mostly not differentially expressed (<6% overlap) and were involved in different biological processes. Differentially spliced genes were primarily enriched for muscle development and functioning, while differentially expressed genes were involved in metabolism, immunity and growth. Furthermore, alternative splicing and gene expression were mostly controlled by independent cis-regulatory quantitative trait loci (<3.4% overlap). Cis-regulatory regions were associated with the parallel divergence in splicing (16.5% of intron clusters) and expression (6.7%-10.1% of differentially expressed genes), indicating shared regulatory variation across ecotype pairs. Contrary to theoretical expectation, we found that differentially spliced genes tended to be highly central in regulatory networks ("hub genes") and were annotated to significantly more gene ontology terms compared to nondifferentially spliced genes, consistent with a higher level of pleiotropy. Together, our results suggest that the concerted regulation of alternative splicing and differential gene expression through different regulatory regions leads to the divergence of complementary processes important for local adaptation. This provides novel insights into the importance of contrasting but putatively complementary molecular processes in rapid parallel adaptive evolution.
Collapse
Affiliation(s)
- Arne Jacobs
- Institute of Biodiversity, Animal Health and Comparative Medicine, College of Medical, Veterinary & Life SciencesUniversity of GlasgowGlasgowUK
- Department of Natural ResourcesCornell UniversityIthacaNYUSA
| | - Kathryn R. Elmer
- Institute of Biodiversity, Animal Health and Comparative Medicine, College of Medical, Veterinary & Life SciencesUniversity of GlasgowGlasgowUK
| |
Collapse
|
36
|
Abstract
The repeated adaptation of oceanic threespine sticklebacks to fresh water has made it a premier organism to study parallel evolution. These small fish have multiple distinct ecotypes that display a wide range of diverse phenotypic traits. Ecotypes are easily crossed in the laboratory, and families are large and develop quickly enough for quantitative trait locus analyses, positioning the threespine stickleback as a versatile model organism to address a wide range of biological questions. Extensive genomic resources, including linkage maps, a high-quality reference genome, and developmental genetics tools have led to insights into the genomic basis of adaptation and the identification of genomic changes controlling traits in vertebrates. Recently, threespine sticklebacks have been used as a model system to identify the genomic basis of highly complex traits, such as behavior and host-microbiome and host-parasite interactions. We review the latest findings and new avenues of research that have led the threespine stickleback to be considered a supermodel of evolutionary genomics.
Collapse
Affiliation(s)
- Kerry Reid
- Department of Ecology and Evolution, Stony Brook University, Stony Brook, New York 11794, USA;
| | - Michael A Bell
- University of California Museum of Paleontology, Berkeley, California 94720, USA
| | - Krishna R Veeramah
- Department of Ecology and Evolution, Stony Brook University, Stony Brook, New York 11794, USA;
| |
Collapse
|
37
|
Benowitz KM, Coleman JM, Allan CW, Matzkin LM. Contributions of cis- and trans-Regulatory Evolution to Transcriptomic Divergence across Populations in the Drosophila mojavensis Larval Brain. Genome Biol Evol 2021; 12:1407-1418. [PMID: 32653899 PMCID: PMC7495911 DOI: 10.1093/gbe/evaa145] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/06/2020] [Indexed: 12/22/2022] Open
Abstract
Natural selection on gene expression was originally predicted to result primarily in cis- rather than trans-regulatory evolution, due to the expectation of reduced pleiotropy. Despite this, numerous studies have ascribed recent evolutionary divergence in gene expression predominantly to trans-regulation. Performing RNA-seq on single isofemale lines from genetically distinct populations of the cactophilic fly Drosophila mojavensis and their F1 hybrids, we recapitulated this pattern in both larval brains and whole bodies. However, we demonstrate that improving the measurement of brain expression divergence between populations by using seven additional genotypes considerably reduces the estimate of trans-regulatory contributions to expression evolution. We argue that the finding of trans-regulatory predominance can result from biases due to environmental variation in expression or other sources of noise, and that cis-regulation is likely a greater contributor to transcriptional evolution across D. mojavensis populations. Lastly, we merge these lines of data to identify several previously hypothesized and intriguing novel candidate genes, and suggest that the integration of regulatory and population-level transcriptomic data can provide useful filters for the identification of potentially adaptive genes.
Collapse
Affiliation(s)
| | - Joshua M Coleman
- Department of Entomology, University of Arizona.,Department of Biological Sciences, University of Alabama in Huntsville
| | | | - Luciano M Matzkin
- Department of Entomology, University of Arizona.,Department of Ecology and Evolutionary Biology, University of Arizona.,BIO5 Institute, University of Arizona
| |
Collapse
|
38
|
Unraveling regulatory divergence, heterotic malleability, and allelic imbalance switching in rice due to drought stress. Sci Rep 2021; 11:13489. [PMID: 34188147 PMCID: PMC8241847 DOI: 10.1038/s41598-021-92938-x] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2021] [Accepted: 06/17/2021] [Indexed: 01/10/2023] Open
Abstract
The indica ecotypes, IR64, an elite drought-susceptible variety adapted to irrigated ecosystem, and Apo (IR55423-01 or NSIC RC9), a moderate drought-tolerant upland genotype together with their hybrid (IR64 × Apo) were exposed to non- and water-stress conditions. By sequencing (RNA-seq) these genotypes, we were able to map genes diverging in cis and/or trans factors. Under non-stress condition, cis dominantly explains (11.2%) regulatory differences, followed by trans (8.9%). Further analysis showed that water-limiting condition largely affects trans and cis + trans factors. On the molecular level, cis and/or trans regulatory divergence explains their genotypic differences and differential drought response. Between the two parental genotypes, Apo appears to exhibit more photosynthetic efficiency even under water-limiting condition and is ascribed to trans. Statistical analyses showed that regulatory divergence is significantly influenced by environmental conditions. Likewise, the mode of parental expression inheritance which drives heterosis (HET) is significantly affected by environmental conditions indicating the malleability of heterosis to external factors. Further analysis revealed that the HET class, dominance, was significantly enriched under water-stress condition. We also identified allelic imbalance switching in which several genes prefer IR64- (or Apo-) specific allele under non-stress condition but switched to Apo- (or IR64-) specific allele when exposed to water-stress condition.
Collapse
|
39
|
Huang Y, Lack JB, Hoppel GT, Pool JE. Parallel and Population-specific Gene Regulatory Evolution in Cold-Adapted Fly Populations. Genetics 2021; 218:6275754. [PMID: 33989401 PMCID: PMC8864734 DOI: 10.1093/genetics/iyab077] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2020] [Accepted: 05/10/2021] [Indexed: 11/15/2022] Open
Abstract
Changes in gene regulation at multiple levels may comprise an important share of the molecular changes underlying adaptive evolution in nature. However, few studies have assayed within- and between-population variation in gene regulatory traits at a transcriptomic scale, and therefore inferences about the characteristics of adaptive regulatory changes have been elusive. Here, we assess quantitative trait differentiation in gene expression levels and alternative splicing (intron usage) between three closely related pairs of natural populations of Drosophila melanogaster from contrasting thermal environments that reflect three separate instances of cold tolerance evolution. The cold-adapted populations were known to show population genetic evidence for parallel evolution at the SNP level, and here we find evidence for parallel expression evolution between them, with stronger parallelism at larval and adult stages than for pupae. We also implement a flexible method to estimate cis- vs trans-encoded contributions to expression or splicing differences at the adult stage. The apparent contributions of cis- vs trans-regulation to adaptive evolution vary substantially among population pairs. While two of three population pairs show a greater enrichment of cis-regulatory differences among adaptation candidates, trans-regulatory differences are more likely to be implicated in parallel expression changes between population pairs. Genes with significant cis-effects are enriched for signals of elevated genetic differentiation between cold- and warm-adapted populations, suggesting that they are potential targets of local adaptation. These findings expand our knowledge of adaptive gene regulatory evolution and our ability to make inferences about this important and widespread process.
Collapse
Affiliation(s)
- Yuheng Huang
- Laboratory of Genetics, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA.,Department of Ecology and Evolutionary Biology, University of California, Irvine, Irvine, CA 92697, USA
| | - Justin B Lack
- Laboratory of Genetics, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA.,Advanced Biomedical Computational Science, Frederick National Laboratory for Cancer Research, Frederick, MD 21701, USA
| | - Grant T Hoppel
- Laboratory of Genetics, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - John E Pool
- Laboratory of Genetics, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| |
Collapse
|
40
|
RNA sequencing describes both population structure and plasticity-selection dynamics in a non-model fish. BMC Genomics 2021; 22:273. [PMID: 33858341 PMCID: PMC8048188 DOI: 10.1186/s12864-021-07592-4] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2021] [Accepted: 04/05/2021] [Indexed: 01/03/2023] Open
Abstract
Background Messenger RNA sequencing is becoming more common in studies of non-model species and is most often used for gene expression-based investigations. However, the method holds potential for numerous other applications as well—including analyses of alternative splicing, population structure, and signatures of selection. To maximize the utility of mRNA data sets, distinct analyses may be combined such as by exploring dynamics between gene expression with signatures of selection in the context of population structure. Here, we compare two published data sets describing two populations of a minnow species endemic to the San Francisco Estuary (Sacramento splittail, Pogonichthys macrolepidotus): a microsatellite data set showing population structure, and an mRNA whole transcriptome data set obtained after the two populations were exposed to a salinity challenge. We compared measures of population structure and genetic variation using single nucleotide polymorphisms (SNPs) called from mRNA from the whole transcriptome sequencing study with those patterns determined from microsatellites. For investigating plasticity and evolution, intra- and inter-population transcriptome plasticity was investigated with differential gene expression, differential exon usage, and gene expression variation. Outlier SNP analysis was also performed on the mRNA data set and signatures of selection and phenotypic plasticity were investigated on an individual-gene basis. Results We found that mRNA sequencing revealed patterns of population structure consistent with those found with microsatellites, but with lower magnitudes of genetic variation and population differentiation consistent with widespread purifying selection expected when using mRNA. In addition, within individual genes, phenotypic plasticity or signatures of selection were found in almost mutual exclusion (except heatr6, nfu1, slc22a6, sya, and mmp13). Conclusions These results show that an mRNA sequencing data set may have multiple uses, including describing population structure and for investigating the mechanistic interplay of evolution and plasticity in adaptation. MRNA sequencing thus complements traditional sequencing methods used for population genetics, in addition to its utility for describing phenotypic plasticity. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-021-07592-4.
Collapse
|
41
|
Gillard GB, Grønvold L, Røsæg LL, Holen MM, Monsen Ø, Koop BF, Rondeau EB, Gundappa MK, Mendoza J, Macqueen DJ, Rohlfs RV, Sandve SR, Hvidsten TR. Comparative regulomics supports pervasive selection on gene dosage following whole genome duplication. Genome Biol 2021; 22:103. [PMID: 33849620 PMCID: PMC8042706 DOI: 10.1186/s13059-021-02323-0] [Citation(s) in RCA: 35] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2020] [Accepted: 03/23/2021] [Indexed: 12/13/2022] Open
Abstract
BACKGROUND Whole genome duplication (WGD) events have played a major role in eukaryotic genome evolution, but the consequence of these extreme events in adaptive genome evolution is still not well understood. To address this knowledge gap, we used a comparative phylogenetic model and transcriptomic data from seven species to infer selection on gene expression in duplicated genes (ohnologs) following the salmonid WGD 80-100 million years ago. RESULTS We find rare cases of tissue-specific expression evolution but pervasive expression evolution affecting many tissues, reflecting strong selection on maintenance of genome stability following genome doubling. Ohnolog expression levels have evolved mostly asymmetrically, by diverting one ohnolog copy down a path towards lower expression and possible pseudogenization. Loss of expression in one ohnolog is significantly associated with transposable element insertions in promoters and likely driven by selection on gene dosage including selection on stoichiometric balance. We also find symmetric expression shifts, and these are associated with genes under strong evolutionary constraints such as ribosome subunit genes. This possibly reflects selection operating to achieve a gene dose reduction while avoiding accumulation of "toxic mutations". Mechanistically, ohnolog regulatory divergence is dictated by the number of bound transcription factors in promoters, with transposable elements being one likely source of novel binding sites driving tissue-specific gains in expression. CONCLUSIONS Our results imply pervasive adaptive expression evolution following WGD to overcome the immediate challenges posed by genome doubling and to exploit the long-term genetic opportunities for novel phenotype evolution.
Collapse
Affiliation(s)
- Gareth B. Gillard
- Faculty of Chemistry, Biotechnology and Food Science, Norwegian University of Life Sciences, Ås, Norway
| | - Lars Grønvold
- Center for Integrative Genetics, Department of Animal and Aquacultural Sciences, Faculty of Biosciences, Norwegian University of Life Sciences, Ås, Norway
| | - Line L. Røsæg
- Center for Integrative Genetics, Department of Animal and Aquacultural Sciences, Faculty of Biosciences, Norwegian University of Life Sciences, Ås, Norway
| | - Matilde Mengkrog Holen
- Center for Integrative Genetics, Department of Animal and Aquacultural Sciences, Faculty of Biosciences, Norwegian University of Life Sciences, Ås, Norway
| | - Øystein Monsen
- Center for Integrative Genetics, Department of Animal and Aquacultural Sciences, Faculty of Biosciences, Norwegian University of Life Sciences, Ås, Norway
| | - Ben F. Koop
- Department of Biology, University of Victoria, Victoria, Canada
| | - Eric B. Rondeau
- Department of Biology, University of Victoria, Victoria, Canada
| | - Manu Kumar Gundappa
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, The University of Edinburgh, Midlothian, UK
| | - John Mendoza
- Department of Computer Science, San Francisco State University, San Francisco, USA
| | - Daniel J. Macqueen
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, The University of Edinburgh, Midlothian, UK
| | - Rori V. Rohlfs
- Department of Biology, San Francisco State University, San Francisco, USA
| | - Simen R. Sandve
- Center for Integrative Genetics, Department of Animal and Aquacultural Sciences, Faculty of Biosciences, Norwegian University of Life Sciences, Ås, Norway
| | - Torgeir R. Hvidsten
- Faculty of Chemistry, Biotechnology and Food Science, Norwegian University of Life Sciences, Ås, Norway
| |
Collapse
|
42
|
Hu J, Wuitchik SJS, Barry TN, Jamniczky HA, Rogers SM, Barrett RDH. Heritability of DNA methylation in threespine stickleback (Gasterosteus aculeatus). Genetics 2021; 217:1-15. [PMID: 33683369 PMCID: PMC8045681 DOI: 10.1093/genetics/iyab001] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2020] [Accepted: 12/30/2020] [Indexed: 12/13/2022] Open
Abstract
Epigenetic mechanisms underlying phenotypic change are hypothesized to contribute to population persistence and adaptation in the face of environmental change. To date, few studies have explored the heritability of intergenerationally stable methylation levels in natural populations, and little is known about the relative contribution of cis- and trans-regulatory changes to methylation variation. Here, we explore the heritability of DNA methylation, and conduct methylation quantitative trait loci (meQTLs) analysis to investigate the genetic architecture underlying methylation variation between marine and freshwater ecotypes of threespine stickleback (Gasterosteus aculeatus). We quantitatively measured genome-wide DNA methylation in fin tissue using reduced representation bisulfite sequencing of F1 and F2 crosses, and their marine and freshwater source populations. We identified cytosines (CpG sites) that exhibited stable methylation levels across generations. We found that additive genetic variance explained an average of 24-35% of the methylation variance, with a number of CpG sites possibly autonomous from genetic control. We also detected both cis- and trans-meQTLs, with only trans-meQTLs overlapping with previously identified genomic regions of high differentiation between marine and freshwater ecotypes. Finally, we identified the genetic architecture underlying two key CpG sites that were differentially methylated between ecotypes. These findings demonstrate a potential role for DNA methylation in facilitating adaptation to divergent environments and improve our understanding of the heritable basis of population epigenomic variation.
Collapse
Affiliation(s)
- Juntao Hu
- National Observation and Research Station for Yangtze Estuarine Wetland Ecosystems, and Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, School of Life Sciences, Institute of Biodiversity Science, Fudan University, Shanghai 200438, China
- Redpath Museum and Department of Biology, McGill University, Montreal, QC H3A 0C4, Canada
| | - Sara J S Wuitchik
- Informatics Group, Harvard University, Cambridge, MA 02138, USA
- Department of Biology, Boston University, Boston, MA 02215, USA
- Department of Biological Sciences, University of Calgary, Calgary, AB T2N 1N4, Canada
| | - Tegan N Barry
- Department of Biological Sciences, University of Calgary, Calgary, AB T2N 1N4, Canada
| | - Heather A Jamniczky
- Department of Cell Biology and Anatomy, McCaig Institute for Bone and Joint Health, University of Calgary, Calgary, AB T2N 1N4, Canada
| | - Sean M Rogers
- Department of Biological Sciences, University of Calgary, Calgary, AB T2N 1N4, Canada
| | - Rowan D H Barrett
- Redpath Museum and Department of Biology, McGill University, Montreal, QC H3A 0C4, Canada
| |
Collapse
|
43
|
Wos G, Bohutínská M, Nosková J, Mandáková T, Kolář F. Parallelism in gene expression between foothill and alpine ecotypes in Arabidopsis arenosa. THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2021; 105:1211-1224. [PMID: 33258160 DOI: 10.1111/tpj.15105] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/15/2020] [Revised: 11/13/2020] [Accepted: 11/26/2020] [Indexed: 06/12/2023]
Abstract
Parallel adaptation results from the independent evolution of similar traits between closely related lineages and allows us to test to what extent evolution is repeatable. Similar gene expression changes are often detected but the identity of genes shaped by parallel selection and the causes of expression parallelism remain largely unknown. By comparing genomes and transcriptomes of four distinct foothill-alpine population pairs across four treatments, we addressed the genetic underpinnings, plasticity and functional consequences of gene expression parallelism in alpine adaptation. Seeds of eight populations of Arabidopsis arenosa were raised under four treatments that differed in temperature and irradiance, factors varying strongly with elevation. Parallelism in differential gene expression between the foothill and alpine ecotypes was quantified by RNA-seq in leaves of young plants. By manipulating temperature and irradiance, we also tested for parallelism in plasticity (i.e., gene-environment interaction, GEI). In spite of global non-parallel patterns transcriptome wide, we found significant parallelism in gene expression at the level of individual loci with an over-representation of genes involved in biotic stress response. In addition, we demonstrated significant parallelism in GEI, indicating a shared differential response of the originally foothill versus alpine populations to environmental variation across mountain regions. A fraction of genes showing expression parallelism also encompassed parallel outliers for genomic differentiation, with greater enrichment of such variants in cis-regulatory elements in some mountain regions. In summary, our results suggest frequent evolutionary repeatability in gene expression changes associated with the colonization of a challenging environment that combines constitutive expression differences and plastic interaction with the surrounding environment.
Collapse
Affiliation(s)
- Guillaume Wos
- Department of Botany, Charles University, Prague, 128 01, Czech Republic
| | - Magdalena Bohutínská
- Department of Botany, Charles University, Prague, 128 01, Czech Republic
- Institute of Botany, The Czech Academy of Sciences, Průhonice, 252 43, Czech Republic
| | - Jana Nosková
- Department of Botany, Charles University, Prague, 128 01, Czech Republic
| | - Terezie Mandáková
- Central European Institute of Technology and Faculty of Science, Masaryk University, Brno, 625 00, Czech Republic
| | - Filip Kolář
- Department of Botany, Charles University, Prague, 128 01, Czech Republic
- Institute of Botany, The Czech Academy of Sciences, Průhonice, 252 43, Czech Republic
| |
Collapse
|
44
|
Yocca AE, Lu Z, Schmitz RJ, Freeling M, Edger PP. Evolution of Conserved Noncoding Sequences in Arabidopsis thaliana. Mol Biol Evol 2021; 38:2692-2703. [PMID: 33565589 PMCID: PMC8233505 DOI: 10.1093/molbev/msab042] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
Recent pangenome studies have revealed a large fraction of the gene content within a species exhibits presence-absence variation (PAV). However, coding regions alone provide an incomplete assessment of functional genomic sequence variation at the species level. Little to no attention has been paid to noncoding regulatory regions in pangenome studies, though these sequences directly modulate gene expression and phenotype. To uncover regulatory genetic variation, we generated chromosome-scale genome assemblies for thirty Arabidopsis thaliana accessions from multiple distinct habitats and characterized species level variation in Conserved Noncoding Sequences (CNS). Our analyses uncovered not only PAV and positional variation (PosV) but that diversity in CNS is nonrandom, with variants shared across different accessions. Using evolutionary analyses and chromatin accessibility data, we provide further evidence supporting roles for conserved and variable CNS in gene regulation. Additionally, our data suggests that transposable elements contribute to CNS variation. Characterizing species-level diversity in all functional genomic sequences may later uncover previously unknown mechanistic links between genotype and phenotype.
Collapse
Affiliation(s)
- Alan E Yocca
- Department of Plant Biology, Michigan State University, East Lansing, MI, USA.,Department of Horticulture, Michigan State University, East Lansing, MI, USA
| | - Zefu Lu
- Department of Genetics, University of Georgia, Athens, GA, USA
| | | | - Michael Freeling
- Department of Plant and Microbial Biology, University of California, Berkeley, CA, USA
| | - Patrick P Edger
- Department of Horticulture, Michigan State University, East Lansing, MI, USA.,Ecology, Evolutionary Biology and Behavior, Michigan State University, East Lansing, MI, USA
| |
Collapse
|
45
|
McGirr JA, Martin CH. Few Fixed Variants between Trophic Specialist Pupfish Species Reveal Candidate Cis-Regulatory Alleles Underlying Rapid Craniofacial Divergence. Mol Biol Evol 2021; 38:405-423. [PMID: 32877534 PMCID: PMC7826174 DOI: 10.1093/molbev/msaa218] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
Investigating closely related species that rapidly evolved divergent feeding morphology is a powerful approach to identify genetic variation underlying variation in complex traits. This can also lead to the discovery of novel candidate genes influencing natural and clinical variation in human craniofacial phenotypes. We combined whole-genome resequencing of 258 individuals with 50 transcriptomes to identify candidate cis-acting genetic variation underlying rapidly evolving craniofacial phenotypes within an adaptive radiation of Cyprinodon pupfishes. This radiation consists of a dietary generalist species and two derived trophic niche specialists-a molluscivore and a scale-eating species. Despite extensive morphological divergence, these species only diverged 10 kya and produce fertile hybrids in the laboratory. Out of 9.3 million genome-wide SNPs and 80,012 structural variants, we found very few alleles fixed between species-only 157 SNPs and 87 deletions. Comparing gene expression across 38 purebred F1 offspring sampled at three early developmental stages, we identified 17 fixed variants within 10 kb of 12 genes that were highly differentially expressed between species. By measuring allele-specific expression in F1 hybrids from multiple crosses, we found that the majority of expression divergence between species was explained by trans-regulatory mechanisms. We also found strong evidence for two cis-regulatory alleles affecting expression divergence of two genes with putative effects on skeletal development (dync2li1 and pycr3). These results suggest that SNPs and structural variants contribute to the evolution of novel traits and highlight the utility of the San Salvador Island pupfish system as an evolutionary model for craniofacial development.
Collapse
Affiliation(s)
- Joseph A McGirr
- Environmental Toxicology Department, University of California, Davis, CA
| | - Christopher H Martin
- Department of Integrative Biology and Museum of Vertebrate Zoology, University of California, Berkeley, CA
| |
Collapse
|
46
|
Thorstensen MJ, Jeffrey JD, Treberg JR, Watkinson DA, Enders EC, Jeffries KM. Genomic signals found using RNA sequencing show signatures of selection and subtle population differentiation in walleye ( Sander vitreus) in a large freshwater ecosystem. Ecol Evol 2020; 10:7173-7188. [PMID: 32760520 PMCID: PMC7391302 DOI: 10.1002/ece3.6418] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2020] [Revised: 05/04/2020] [Accepted: 05/07/2020] [Indexed: 12/29/2022] Open
Abstract
RNA sequencing is an effective approach for studying aquatic species yielding both physiological and genomic data. However, its population genetic applications are not well-characterized. We investigate this possible role for RNA sequencing for population genomics in Lake Winnipeg, Manitoba, Canada, walleye (Sander vitreus). Lake Winnipeg walleye represent the largest component of the second-largest freshwater fishery in Canada. In the present study, large female walleye were sampled via nonlethal gill biopsy over two years at three spawning sites representing a latitudinal gradient in the lake. Genetic variation from sequenced mRNA was analyzed for neutral and adaptive markers to investigate population structure and possible adaptive variation. We find low population divergence (F ST = 0.0095), possible northward gene flow, and outlier loci that vary latitudinally in transcripts associated with cell membrane proteins and cytoskeletal function. These results indicate that Lake Winnipeg walleye may be effectively managed as a single demographically connected metapopulation with contributing subpopulations and suggest genomic differences possibly underlying observed phenotypic differences. Despite its high cost relative to other genotyping methods, RNA sequencing data can yield physiological in addition to genetic information discussed here. We therefore argue that it is useful for addressing diverse molecular questions in the conservation of freshwater species.
Collapse
Affiliation(s)
| | | | - Jason R. Treberg
- Department of Biological SciencesUniversity of ManitobaWinnipegMBCanada
| | | | - Eva C. Enders
- Freshwater Institute, Fisheries and Oceans CanadaWinnipegMBCanada
| | - Ken M. Jeffries
- Department of Biological SciencesUniversity of ManitobaWinnipegMBCanada
| |
Collapse
|
47
|
Freedman AH, Gaspar JM, Sackton TB. Short paired-end reads trump long single-end reads for expression analysis. BMC Bioinformatics 2020; 21:149. [PMID: 32306895 PMCID: PMC7168855 DOI: 10.1186/s12859-020-3484-z] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2019] [Accepted: 04/07/2020] [Indexed: 01/02/2023] Open
Abstract
Background Typical experimental design advice for expression analyses using RNA-seq generally assumes that single-end reads provide robust gene-level expression estimates in a cost-effective manner, and that the additional benefits obtained from paired-end sequencing are not worth the additional cost. However, in many cases (e.g., with Illumina NextSeq and NovaSeq instruments), shorter paired-end reads and longer single-end reads can be generated for the same cost, and it is not obvious which strategy should be preferred. Using publicly available data, we test whether short-paired end reads can achieve more robust expression estimates and differential expression results than single-end reads of approximately the same total number of sequenced bases. Results At both the transcript and gene levels, 2 × 40 paired-end reads unequivocally provide expression estimates that are more highly correlated with 2 × 125 than 1 × 75 reads; in nearly all cases, those correlations are also greater than for 1 × 125, despite the greater total number of sequenced bases for the latter. Across an array of metrics, differential expression tests based upon 2 × 40 consistently outperform those using 1 × 75. Conclusion Researchers seeking a cost-effective approach for gene-level expression analysis should prefer short paired-end reads over a longer single-end strategy. Short paired-end reads will also give reasonably robust expression estimates and differential expression results at the isoform level.
Collapse
Affiliation(s)
- Adam H Freedman
- Informatics Group, Harvard University, Cambridge, MA, 02138, USA
| | - John M Gaspar
- Informatics Group, Harvard University, Cambridge, MA, 02138, USA
| | | |
Collapse
|
48
|
Verta JP, Jones F. mRNA Extraction from Gill Tissue for RNA-sequencing. Bio Protoc 2020; 10:e3539. [PMID: 33659513 PMCID: PMC7842710 DOI: 10.21769/bioprotoc.3539] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2019] [Revised: 01/07/2020] [Accepted: 01/15/2020] [Indexed: 11/10/2022] Open
Abstract
Adaptation is thought to proceed in part through spatial and temporal changes in gene expression. Fish species such as the threespine stickleback are powerful vertebrate models to study the genetic architecture of adaptive changes in gene expression since divergent adaptation to different environments is common, they are abundant and easy to study in the wild and lab, and have well-established genetic and genomic resources. Fish gills, due to their respiratory and osmoregulatory roles, show many physiological adaptations to local water chemistry, including differences in gene expression. However, obtaining high-quality RNA using popular column-based extraction methods can be challenging from small tissue samples high in cartilage and bone such as fish gills. Here, we describe a bead-based mRNA extraction and transcriptome RNA-seq protocol that does not use purification columns. The protocol can be readily scaled according to sample size for the purposes of diverse gene expression experiments using animal or plant tissue.
Collapse
Affiliation(s)
- Jukka-Pekka Verta
- Organismal and Evolutionary Biology Research Programme, University of Helsinki, Helsinki, Finland.,Friedrich Miescher Laboratory of the Max Planck Society, Tuebingen, Germany
| | - Felicity Jones
- Friedrich Miescher Laboratory of the Max Planck Society, Tuebingen, Germany
| |
Collapse
|
49
|
Huang Y, Feulner PGD, Eizaguirre C, Lenz TL, Bornberg-Bauer E, Milinski M, Reusch TBH, Chain FJJ. Genome-Wide Genotype-Expression Relationships Reveal Both Copy Number and Single Nucleotide Differentiation Contribute to Differential Gene Expression between Stickleback Ecotypes. Genome Biol Evol 2020; 11:2344-2359. [PMID: 31298693 PMCID: PMC6735750 DOI: 10.1093/gbe/evz148] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/10/2019] [Indexed: 12/11/2022] Open
Abstract
Repeated and independent emergence of trait divergence that matches habitat differences is a sign of parallel evolution by natural selection. Yet, the molecular underpinnings that are targeted by adaptive evolution often remain elusive. We investigate this question by combining genome-wide analyses of copy number variants (CNVs), single nucleotide polymorphisms (SNPs), and gene expression across four pairs of lake and river populations of the three-spined stickleback (Gasterosteus aculeatus). We tested whether CNVs that span entire genes and SNPs occurring in putative cis-regulatory regions contribute to gene expression differences between sticklebacks from lake and river origins. We found 135 gene CNVs that showed a significant positive association between gene copy number and gene expression, suggesting that CNVs result in dosage effects that can fuel phenotypic variation and serve as substrates for habitat-specific selection. Copy number differentiation between lake and river sticklebacks also contributed to expression differences of two immune-related genes in immune tissues, cathepsin A and GIMAP7. In addition, we identified SNPs in cis-regulatory regions (eSNPs) associated with the expression of 1,865 genes, including one eSNP upstream of a carboxypeptidase gene where both the SNP alleles differentiated and the gene was differentially expressed between lake and river populations. Our study highlights two types of mutations as important sources of genetic variation involved in the evolution of gene expression and in potentially facilitating repeated adaptation to novel environments.
Collapse
Affiliation(s)
- Yun Huang
- Department of Evolutionary Ecology, Max Planck Institute for Evolutionary Biology, Plön, Germany.,Biodiversity Research Center, Academia Sinica, Taipei, Taiwan, ROC
| | - Philine G D Feulner
- Department of Fish Ecology and Evolution, Centre of Ecology, Evolution and Biogeochemistry, EAWAG Swiss Federal Institute of Aquatic Science and Technology, Kastanienbaum, Switzerland.,Division of Aquatic Ecology and Evolution, Institute of Ecology and Evolution, University of Bern, Switzerland
| | - Christophe Eizaguirre
- School of Biological and Chemical Sciences, Queen Mary University of London, United Kingdom
| | - Tobias L Lenz
- Department of Evolutionary Ecology, Max Planck Institute for Evolutionary Biology, Plön, Germany
| | - Erich Bornberg-Bauer
- Evolutionary Bioinformatics, Institute for Evolution and Biodiversity, Westfälische Wilhelms University, Münster, Germany
| | - Manfred Milinski
- Department of Evolutionary Ecology, Max Planck Institute for Evolutionary Biology, Plön, Germany
| | - Thorsten B H Reusch
- Marine Evolutionary Ecology, GEOMAR Helmholtz Centre for Ocean Research Kiel, Germany
| | - Frédéric J J Chain
- Department of Biological Sciences, University of Massachusetts Lowell, USA
| |
Collapse
|
50
|
Verta JP, Jones FC. Predominance of cis-regulatory changes in parallel expression divergence of sticklebacks. eLife 2019; 8:e43785. [PMID: 31090544 PMCID: PMC6550882 DOI: 10.7554/elife.43785] [Citation(s) in RCA: 48] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2018] [Accepted: 05/01/2019] [Indexed: 12/15/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Jukka-Pekka Verta
- Friedrich Miescher Laboratory of the Max Planck SocietyMax-Planck-RingTübingenGermany
- Organismal and Evolutionary Biology Research ProgrammeUniversity of HelsinkiHelsinkiFinland
| | - Felicity C Jones
- Friedrich Miescher Laboratory of the Max Planck SocietyMax-Planck-RingTübingenGermany
| |
Collapse
|