1
|
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
|
2
|
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
|
3
|
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
|
4
|
Tong C. Convergent genomics and Arctic adaptation of ruminants. Proc Biol Sci 2024; 291:20232448. [PMID: 38166424 PMCID: PMC10762444 DOI: 10.1098/rspb.2023.2448] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2023] [Accepted: 12/04/2023] [Indexed: 01/04/2024] Open
Affiliation(s)
- Chao Tong
- School of Life Sciences, Arizona State University, Tempe, AZ 85287, USA
| |
Collapse
|
5
|
Square TA, Mackey EJ, Sundaram S, Weksler NC, Chen ZZ, Narayanan SN, Miller CT. Modulation of tooth regeneration through opposing responses to Wnt and BMP signals in teleosts. Development 2023; 150:dev202168. [PMID: 38059590 PMCID: PMC10730089 DOI: 10.1242/dev.202168] [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: 07/13/2023] [Accepted: 11/02/2023] [Indexed: 12/08/2023]
Abstract
Most vertebrate species undergo tooth replacement throughout adult life. This process is marked by the shedding of existing teeth and the regeneration of tooth organs. However, little is known about the genetic circuitry regulating tooth replacement. Here, we tested whether fish orthologs of genes known to regulate mammalian hair regeneration have effects on tooth replacement. Using two fish species that demonstrate distinct modes of tooth regeneration, threespine stickleback (Gasterosteus aculeatus) and zebrafish (Danio rerio), we found that transgenic overexpression of four different genes changed tooth replacement rates in the direction predicted by a hair regeneration model: Wnt10a and Grem2a increased tooth replacement rate, whereas Bmp6 and Dkk2 strongly inhibited tooth formation. Thus, similar to known roles in hair regeneration, Wnt and BMP signals promote and inhibit regeneration, respectively. Regulation of total tooth number was separable from regulation of replacement rates. RNA sequencing of stickleback dental tissue showed that Bmp6 overexpression resulted in an upregulation of Wnt inhibitors. Together, these data support a model in which different epithelial organs, such as teeth and hair, share genetic circuitry driving organ regeneration.
Collapse
Affiliation(s)
- Tyler A. Square
- Department of Molecular & Cell Biology, University of California, Berkeley, CA 94720, USA
| | - Emma J. Mackey
- Department of Molecular & Cell Biology, University of California, Berkeley, CA 94720, USA
| | - Shivani Sundaram
- Department of Molecular & Cell Biology, University of California, Berkeley, CA 94720, USA
| | - Naama C. Weksler
- Department of Molecular & Cell Biology, University of California, Berkeley, CA 94720, USA
| | - Zoe Z. Chen
- Department of Molecular & Cell Biology, University of California, Berkeley, CA 94720, USA
| | - Sujanya N. Narayanan
- Department of Molecular & Cell Biology, University of California, Berkeley, CA 94720, USA
| | - Craig T. Miller
- Department of Molecular & Cell Biology, University of California, Berkeley, CA 94720, USA
| |
Collapse
|
6
|
Rahi ML, Mather PB, de Bello Cioffi M, Ezaz T, Hurwood DA. Genomic Basis of Freshwater Adaptation in the Palaemonid Prawn Genus Macrobrachium: Convergent Evolution Following Multiple Independent Colonization Events. J Mol Evol 2023; 91:976-989. [PMID: 38010517 DOI: 10.1007/s00239-023-10149-6] [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: 10/14/2022] [Accepted: 11/14/2023] [Indexed: 11/29/2023]
Abstract
Adaptation to different salinity environments can enhance morphological and genomic divergence between related aquatic taxa. Species of prawns in the genus Macrobrachium naturally inhabit different osmotic niches and possess distinctive lifecycle traits associated with salinity tolerance. This study was conducted to investigate the patterns of adaptive genomic divergence during freshwater colonization in 34 Macrobrachium species collected from four continents; Australia, Asia, North and South America. Genotyping-by-sequencing (GBS) technique identified 5018 loci containing 82,636 single nucleotide polymorphisms (SNPs) that were used to reconstruct a phylogenomic tree. An additional phylogeny was reconstructed based on 43 candidate genes, previously identified as being potentially associated with freshwater adaptation. Comparison of the two phylogenetic trees revealed contrasting topologies. The GBS tree indicated multiple independent continent-specific invasions into freshwater by Macrobrachium lineages following common marine ancestry, as species with abbreviated larval development (ALD), i.e., species having a full freshwater life history, appeared reciprocally monophyletic within each continent. In contrast, the candidate gene tree showed convergent evolution for all ALD species worldwide, forming a single, well-supported clade. This latter pattern is likely the result of common evolutionary pressures selecting key mutations favored in continental freshwater habitats Results suggest that following multiple independent invasions into continental freshwaters at different evolutionary timescales, Macrobrachium taxa experienced adaptive genomic divergence, and in particular, convergence in the same genomic regions with parallel shifts in specific conserved phenotypic traits, such as evolution of larger eggs with abbreviated larval developmental.
Collapse
Affiliation(s)
- Md Lifat Rahi
- Fisheries and Marine Resource Technology Discipline, Khulna University, Khulna, Bangladesh.
| | - Peter B Mather
- Faculty of Science, Queensland University of Technology (QUT), Brisbane, QLD, 4001, Australia
| | - Marcelo de Bello Cioffi
- Department of Genetics and Evolution, Federal University of Sao Carlos, São Carlos, SP, Brazil
| | - Tariq Ezaz
- Institute for Applied Ecology (IAE), University of Canberra (UC), Canberra, ACT, 2617, Australia
| | - David A Hurwood
- Faculty of Science, Queensland University of Technology (QUT), Brisbane, QLD, 4001, Australia
| |
Collapse
|
7
|
Li C, Huang W, Han X, Zhao G, Zhang W, He W, Nie B, Chen X, Zhang T, Bai W, Zhang X, He J, Zhao C, Fernie AR, Tschaplinski TJ, Yang X, Yan S, Wang L. Diel dynamics of multi-omics in elkhorn fern provide new insights into weak CAM photosynthesis. PLANT COMMUNICATIONS 2023; 4:100594. [PMID: 36960529 PMCID: PMC10504562 DOI: 10.1016/j.xplc.2023.100594] [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: 10/27/2022] [Revised: 02/19/2023] [Accepted: 03/20/2023] [Indexed: 05/29/2023]
Abstract
Crassulacean acid metabolism (CAM) has high water-use efficiency (WUE) and is widely recognized to have evolved from C3 photosynthesis. Different plant lineages have convergently evolved CAM, but the molecular mechanism that underlies C3-to-CAM evolution remains to be clarified. Platycerium bifurcatum (elkhorn fern) provides an opportunity to study the molecular changes underlying the transition from C3 to CAM photosynthesis because both modes of photosynthesis occur in this species, with sporotrophophyll leaves (SLs) and cover leaves (CLs) performing C3 and weak CAM photosynthesis, respectively. Here, we report that the physiological and biochemical attributes of CAM in weak CAM-performing CLs differed from those in strong CAM species. We investigated the diel dynamics of the metabolome, proteome, and transcriptome in these dimorphic leaves within the same genetic background and under identical environmental conditions. We found that multi-omic diel dynamics in P. bifurcatum exhibit both tissue and diel effects. Our analysis revealed temporal rewiring of biochemistry relevant to the energy-producing pathway (TCA cycle), CAM pathway, and stomatal movement in CLs compared with SLs. We also confirmed that PHOSPHOENOLPYRUVATE CARBOXYLASE KINASE (PPCK) exhibits convergence in gene expression among highly divergent CAM lineages. Gene regulatory network analysis identified candidate transcription factors regulating the CAM pathway and stomatal movement. Taken together, our results provide new insights into weak CAM photosynthesis and new avenues for CAM bioengineering.
Collapse
Affiliation(s)
- Cheng Li
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
| | - Wenjie Huang
- Guangdong Key Laboratory for Crop Germplasm Resources Preservation and Utilization, Agro-biological Gene Research Center, Guangdong Academy of Agricultural Sciences, Guangzhou, China
| | - Xiaoxu Han
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
| | - Guohua Zhao
- Key Laboratory of Southern Subtropical Plant Diversity, Fairy Lake Botanical Garden, Shenzhen & Chinese Academy of Science, Shenzhen, China
| | - Wenyang Zhang
- Guangdong Key Laboratory for Crop Germplasm Resources Preservation and Utilization, Agro-biological Gene Research Center, Guangdong Academy of Agricultural Sciences, Guangzhou, China
| | - Weijun He
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
| | - Bao Nie
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
| | - Xufeng Chen
- Guangdong Key Laboratory for Crop Germplasm Resources Preservation and Utilization, Agro-biological Gene Research Center, Guangdong Academy of Agricultural Sciences, Guangzhou, China
| | - Taijie Zhang
- Institute of Plant Protection, Guangdong Academy of Agricultural Sciences, Key Laboratory of Green Prevention and Control on Fruits and Vegetables in South China, Ministry of Agriculture and Rural Affairs, Guangdong Provincial Key Laboratory of High Technology for Plant Protection, Guangzhou, China
| | - Wenhui Bai
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
| | - Xiaopeng Zhang
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
| | - Jingjing He
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
| | - Cheng Zhao
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
| | - Alisdair R Fernie
- Max Planck Institute of Molecular Plant Physiology, Am Muhlenberg 1, 14476 Potsdam-Golm, Germany
| | - Timothy J Tschaplinski
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA; The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA
| | - Xiaohan Yang
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA; The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA.
| | - Shijuan Yan
- Guangdong Key Laboratory for Crop Germplasm Resources Preservation and Utilization, Agro-biological Gene Research Center, Guangdong Academy of Agricultural Sciences, Guangzhou, China.
| | - Li Wang
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China; Kunpeng Institute of Modern Agriculture at Foshan, Chinese Academy of Agricultural Sciences, Foshan, China.
| |
Collapse
|
8
|
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
|
9
|
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
|
10
|
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
|
11
|
Stepaniak MD, Square TA, Miller CT. Evolved Bmp6 enhancer alleles drive spatial shifts in gene expression during tooth development in sticklebacks. Genetics 2021; 219:6374454. [PMID: 34849839 PMCID: PMC8664583 DOI: 10.1093/genetics/iyab151] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2021] [Accepted: 08/31/2021] [Indexed: 11/13/2022] Open
Abstract
Mutations in enhancers have been shown to often underlie natural variation but the evolved differences in enhancer activity can be difficult to identify in vivo. Threespine sticklebacks (Gasterosteus aculeatus) are a robust system for studying enhancer evolution due to abundant natural genetic variation, a diversity of evolved phenotypes between ancestral marine and derived freshwater forms, and the tractability of transgenic techniques. Previous work identified a series of polymorphisms within an intronic enhancer of the Bone morphogenetic protein 6 (Bmp6) gene that are associated with evolved tooth gain, a derived increase in freshwater tooth number that arises late in development. Here, we use a bicistronic reporter construct containing a genetic insulator and a pair of reciprocal two-color transgenic reporter lines to compare enhancer activity of marine and freshwater alleles of this enhancer. In older fish, the two alleles drive partially overlapping expression in both mesenchyme and epithelium of developing teeth, but the freshwater enhancer drives a reduced mesenchymal domain and a larger epithelial domain relative to the marine enhancer. In younger fish, these spatial shifts in enhancer activity are less pronounced. Comparing Bmp6 expression by in situ hybridization in developing teeth of marine and freshwater fish reveals similar evolved spatial shifts in gene expression. Together, these data support a model in which the polymorphisms within this enhancer underlie evolved tooth gain by shifting the spatial expression of Bmp6 during tooth development, and provide a general strategy to identify spatial differences in enhancer activity in vivo.
Collapse
Affiliation(s)
- Mark D Stepaniak
- Department of Molecular and Cell Biology, University of California, Berkeley, Berkeley, CA 94720, USA
| | - Tyler A Square
- Department of Molecular and Cell Biology, University of California, Berkeley, Berkeley, CA 94720, USA
| | - Craig T Miller
- Department of Molecular and Cell Biology, University of California, Berkeley, Berkeley, CA 94720, USA
| |
Collapse
|
12
|
Lyu MJA, Gowik U, Kelly S, Covshoff S, Hibberd JM, Sage RF, Ludwig M, Wong GKS, Westhoff P, Zhu XG. The coordination of major events in C 4 photosynthesis evolution in the genus Flaveria. Sci Rep 2021; 11:15618. [PMID: 34341365 PMCID: PMC8329263 DOI: 10.1038/s41598-021-93381-8] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2021] [Accepted: 05/31/2021] [Indexed: 12/13/2022] Open
Abstract
C4 photosynthesis is a remarkable complex trait, elucidations of the evolutionary trajectory of C4 photosynthesis from its ancestral C3 pathway can help us better understand the generic principles of the evolution of complex traits and guide the engineering of C3 crops for higher yields. Here, we used the genus Flaveria that contains C3, C3-C4, C4-like and C4 species as a system to study the evolution of C4 photosynthesis. We first mapped transcript abundance, protein sequence and morphological features onto the phylogenetic tree of the genus Flaveria, and calculated the evolutionary correlation of different features; we then predicted the relative changes of ancestral nodes of those features to illustrate the major events during the evolution of C4 photosynthesis. We found that gene expression and protein sequence showed consistent modification patterns in the phylogenetic tree. High correlation coefficients ranging from 0.46 to 0.9 among gene expression, protein sequence and morphology were observed. The greatest modification of those different features consistently occurred at the transition between C3-C4 species and C4-like species. Our results show highly coordinated changes in gene expression, protein sequence and morphological features, which support evolutionary major events during the evolution of C4 metabolism.
Collapse
Affiliation(s)
- Ming-Ju Amy Lyu
- National Key Laboratory of Plant Molecular Genetics, CAS Center for Excellence in Molecular Plant Sciences, Institute of Plant Physiology and Ecology, Chinese Academy of Sciences, Shanghai, China
| | - Udo Gowik
- Institute of Plant Molecular and Developmental Biology, Heinrich-Heine-University, Dusseldorf, Germany
| | - Steve Kelly
- Department of Plant Sciences, University of Oxford, Oxford, UK
| | - Sarah Covshoff
- Department of Plant Sciences, University of Cambridge, Cambridge, UK
| | - Julian M Hibberd
- Department of Plant Sciences, University of Cambridge, Cambridge, UK
| | - Rowan F Sage
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, Canada
| | - Martha Ludwig
- School of Molecular Sciences, University of Western Australia, Crawley, WA, Australia
| | - Gane Ka-Shu Wong
- BGI-Shenzhen, Beishan Industrial Zone, Yantian District, Shenzhen, 518083, China
- Department of Medicine and Department of Biological Sciences, The University of Alberta, Edmonton, AB, T6G 2E1, Canada
| | - Peter Westhoff
- Institute of Plant Molecular and Developmental Biology, Heinrich-Heine-University, Dusseldorf, Germany
| | - Xin-Guang Zhu
- National Key Laboratory of Plant Molecular Genetics, CAS Center for Excellence in Molecular Plant Sciences, Institute of Plant Physiology and Ecology, Chinese Academy of Sciences, Shanghai, China.
| |
Collapse
|
13
|
Waters JM, McCulloch GA. Reinventing the wheel? Reassessing the roles of gene flow, sorting and convergence in repeated evolution. Mol Ecol 2021; 30:4162-4172. [PMID: 34133810 DOI: 10.1111/mec.16018] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2021] [Revised: 06/08/2021] [Accepted: 06/11/2021] [Indexed: 12/31/2022]
Abstract
Biologists have long been intrigued by apparently predictable and repetitive evolutionary trajectories inferred across a variety of lineages and systems. In recent years, high-throughput sequencing analyses have started to transform our understanding of such repetitive shifts. While researchers have traditionally categorized such shifts as either "convergent" or "parallel," based on relatedness of the lineages involved, emerging genomic insights provide an opportunity to better describe the actual evolutionary mechanisms at play. A synthesis of recent genomic analyses confirms that convergence is the predominant driver of repetitive evolution among species, whereas repeated sorting of standing variation is the major driver of repeated shifts within species. However, emerging data reveal numerous notable exceptions to these expectations, with recent examples of de novo mutations underpinning convergent shifts among even very closely related lineages, while repetitive sorting processes have occurred among even deeply divergent taxa, sometimes via introgression. A number of very recent analyses have found evidence for both processes occurring on different scales within taxa. We suggest that the relative importance of convergent versus sorting processes depends on the interplay between gene flow among populations, and phylogenetic relatedness of the lineages involved.
Collapse
|
14
|
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
|
15
|
Karagic N, Schneider RF, Meyer A, Hulsey CD. A Genomic Cluster Containing Novel and Conserved Genes is Associated with Cichlid Fish Dental Developmental Convergence. Mol Biol Evol 2021; 37:3165-3174. [PMID: 32579214 DOI: 10.1093/molbev/msaa153] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
The two toothed jaws of cichlid fishes provide textbook examples of convergent evolution. Tooth phenotypes such as enlarged molar-like teeth used to process hard-shelled mollusks have evolved numerous times independently during cichlid diversification. Although the ecological benefit of molar-like teeth to crush prey is known, it is unclear whether the same molecular mechanisms underlie these convergent traits. To identify genes involved in the evolution and development of enlarged cichlid teeth, we performed RNA-seq on the serially homologous-toothed oral and pharyngeal jaws as well as the fourth toothless gill arch of Astatoreochromis alluaudi. We identified 27 genes that are highly upregulated on both tooth-bearing jaws compared with the toothless gill arch. Most of these genes have never been reported to play a role in tooth formation. Two of these genes (unk, rpfA) are not found in other vertebrate genomes but are present in all cichlid genomes. They also cluster genomically with two other highly expressed tooth genes (odam, scpp5) that exhibit conserved expression during vertebrate odontogenesis. Unk and rpfA were confirmed via in situ hybridization to be expressed in developing teeth of Astatotilapia burtoni. We then examined expression of the cluster's four genes in six evolutionarily independent and phylogenetically disparate cichlid species pairs each with a large- and a small-toothed species. Odam and unk commonly and scpp5 and rpfA always showed higher expression in larger toothed cichlid jaws. Convergent trophic adaptations across cichlid diversity are associated with the repeated developmental deployment of this genomic cluster containing conserved and novel cichlid-specific genes.
Collapse
Affiliation(s)
- Nidal Karagic
- Department for Zoology and Evolutionary Biology, University of Konstanz, Konstanz, Germany
| | - Ralf F Schneider
- Department for Zoology and Evolutionary Biology, University of Konstanz, Konstanz, Germany
| | - Axel Meyer
- Department for Zoology and Evolutionary Biology, University of Konstanz, Konstanz, Germany
| | - C Darrin Hulsey
- Department for Zoology and Evolutionary Biology, University of Konstanz, Konstanz, Germany
| |
Collapse
|
16
|
Square TA, Sundaram S, Mackey EJ, Miller CT. Distinct tooth regeneration systems deploy a conserved battery of genes. EvoDevo 2021; 12:4. [PMID: 33766133 PMCID: PMC7995769 DOI: 10.1186/s13227-021-00172-3] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2020] [Accepted: 01/13/2021] [Indexed: 01/01/2023] Open
Abstract
BACKGROUND Vertebrate teeth exhibit a wide range of regenerative systems. Many species, including most mammals, reptiles, and amphibians, form replacement teeth at a histologically distinct location called the successional dental lamina, while other species do not employ such a system. Notably, a 'lamina-less' tooth replacement condition is found in a paraphyletic array of ray-finned fishes, such as stickleback, trout, cod, medaka, and bichir. Furthermore, the position, renewal potential, and latency times appear to vary drastically across different vertebrate tooth regeneration systems. The progenitor cells underlying tooth regeneration thus present highly divergent arrangements and potentials. Given the spectrum of regeneration systems present in vertebrates, it is unclear if morphologically divergent tooth regeneration systems deploy an overlapping battery of genes in their naïve dental tissues. RESULTS In the present work, we aimed to determine whether or not tooth progenitor epithelia could be composed of a conserved cell type between vertebrate dentitions with divergent regeneration systems. To address this question, we compared the pharyngeal tooth regeneration processes in two ray-finned fishes: zebrafish (Danio rerio) and threespine stickleback (Gasterosteus aculeatus). These two teleost species diverged approximately 250 million years ago and demonstrate some stark differences in dental morphology and regeneration. Here, we find that the naïve successional dental lamina in zebrafish expresses a battery of nine genes (bmpr1aa, bmp6, cd34, gli1, igfbp5a, lgr4, lgr6, nfatc1, and pitx2), while active Wnt signaling and Lef1 expression occur during early morphogenesis stages of tooth development. We also find that, despite the absence of a histologically distinct successional dental lamina in stickleback tooth fields, the same battery of nine genes (Bmpr1a, Bmp6, CD34, Gli1, Igfbp5a, Lgr4, Lgr6, Nfatc1, and Pitx2) are expressed in the basalmost endodermal cell layer, which is the region most closely associated with replacement tooth germs. Like zebrafish, stickleback replacement tooth germs additionally express Lef1 and exhibit active Wnt signaling. Thus, two fish systems that either have an organized successional dental lamina (zebrafish) or lack a morphologically distinct successional dental lamina (sticklebacks) deploy similar genetic programs during tooth regeneration. CONCLUSIONS We propose that the expression domains described here delineate a highly conserved "successional dental epithelium" (SDE). Furthermore, a set of orthologous genes is known to mark hair follicle epithelial stem cells in mice, suggesting that regenerative systems in other epithelial appendages may utilize a related epithelial progenitor cell type, despite the highly derived nature of the resulting functional organs.
Collapse
Affiliation(s)
- Tyler A Square
- Department of Molecular & Cell Biology, University of California, Berkeley, USA.
| | - Shivani Sundaram
- Department of Molecular & Cell Biology, University of California, Berkeley, USA
| | - Emma J Mackey
- Department of Molecular & Cell Biology, University of California, Berkeley, USA
| | - Craig T Miller
- Department of Molecular & Cell Biology, University of California, Berkeley, USA.
| |
Collapse
|
17
|
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
|
18
|
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
|
19
|
Wang Y, Zhao Y, Wang Y, Li Z, Guo B, Merilä J. Population transcriptomics reveals weak parallel genetic basis in repeated marine and freshwater divergence in nine‐spined sticklebacks. Mol Ecol 2020; 29:1642-1656. [DOI: 10.1111/mec.15435] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2020] [Revised: 03/21/2020] [Accepted: 03/30/2020] [Indexed: 12/20/2022]
Affiliation(s)
- Yingnan Wang
- Key Laboratory of Zoological Systematics and Evolution Institute of Zoology Chinese Academy of Sciences Beijing China
- University of Chinese Academy of Sciences Beijing China
| | - Yongxin Zhao
- Key Laboratory of Zoological Systematics and Evolution Institute of Zoology Chinese Academy of Sciences Beijing China
| | - Yu Wang
- Key Laboratory of Zoological Systematics and Evolution Institute of Zoology Chinese Academy of Sciences Beijing China
- University of Chinese Academy of Sciences Beijing China
| | - Zitong Li
- Ecological Genetics Research Unit, Organismal and Evolutionary Biology Research Programme Faculty of Biological and Environmental Sciences University of Helsinki Helsinki Finland
| | - Baocheng Guo
- Key Laboratory of Zoological Systematics and Evolution Institute of Zoology Chinese Academy of Sciences Beijing China
- University of Chinese Academy of Sciences Beijing China
- Center for Excellence in Animal Evolution and Genetics Chinese Academy of Sciences Kunming China
| | - Juha Merilä
- Ecological Genetics Research Unit, Organismal and Evolutionary Biology Research Programme Faculty of Biological and Environmental Sciences University of Helsinki Helsinki Finland
| |
Collapse
|
20
|
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
|
21
|
Tan MP, Wong LL, Razali SA, Afiqah-Aleng N, Mohd Nor SA, Sung YY, Van de Peer Y, Sorgeloos P, Danish-Daniel M. Applications of Next-Generation Sequencing Technologies and Computational Tools in Molecular Evolution and Aquatic Animals Conservation Studies: A Short Review. Evol Bioinform Online 2019; 15:1176934319892284. [PMID: 31839703 PMCID: PMC6896124 DOI: 10.1177/1176934319892284] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2019] [Accepted: 11/12/2019] [Indexed: 12/21/2022] Open
Abstract
Aquatic ecosystems that form major biodiversity hotspots are critically threatened due to environmental and anthropogenic stressors. We believe that, in this genomic era, computational methods can be applied to promote aquatic biodiversity conservation by addressing questions related to the evolutionary history of aquatic organisms at the molecular level. However, huge amounts of genomics data generated can only be discerned through the use of bioinformatics. Here, we examine the applications of next-generation sequencing technologies and bioinformatics tools to study the molecular evolution of aquatic animals and discuss the current challenges and future perspectives of using bioinformatics toward aquatic animal conservation efforts.
Collapse
Affiliation(s)
- Min Pau Tan
- Institute of Marine Biotechnology, Universiti Malaysia Terengganu, Kuala Nerus, Malaysia.,Faculty of Fisheries and Food Science, Universiti Malaysia Terengganu, Kuala Nerus, Malaysia
| | - Li Lian Wong
- Institute of Marine Biotechnology, Universiti Malaysia Terengganu, Kuala Nerus, Malaysia.,Institute of Tropical Aquaculture, Universiti Malaysia Terengganu, Kuala Nerus, Malaysia
| | - Siti Aisyah Razali
- Institute of Marine Biotechnology, Universiti Malaysia Terengganu, Kuala Nerus, Malaysia
| | - Nor Afiqah-Aleng
- Institute of Marine Biotechnology, Universiti Malaysia Terengganu, Kuala Nerus, Malaysia
| | - Siti Azizah Mohd Nor
- Institute of Marine Biotechnology, Universiti Malaysia Terengganu, Kuala Nerus, Malaysia
| | - Yeong Yik Sung
- Institute of Marine Biotechnology, Universiti Malaysia Terengganu, Kuala Nerus, Malaysia
| | - Yves Van de Peer
- Institute of Marine Biotechnology, Universiti Malaysia Terengganu, Kuala Nerus, Malaysia.,Center for Plant Systems Biology, VIB, Ghent, Belgium.,Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium
| | - Patrick Sorgeloos
- Institute of Marine Biotechnology, Universiti Malaysia Terengganu, Kuala Nerus, Malaysia.,Laboratory of Aquaculture & Artemia Reference Center, Faculty of Bioscience Engineering, Ghent University, Ghent, Belgium
| | - Muhd Danish-Daniel
- Institute of Marine Biotechnology, Universiti Malaysia Terengganu, Kuala Nerus, Malaysia.,Faculty of Fisheries and Food Science, Universiti Malaysia Terengganu, Kuala Nerus, Malaysia
| |
Collapse
|
22
|
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
|
23
|
Kitano J, Ishikawa A, Kusakabe M. Parallel transcriptome evolution in stream threespine sticklebacks. Dev Growth Differ 2018; 61:104-113. [DOI: 10.1111/dgd.12576] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2018] [Revised: 10/05/2018] [Accepted: 10/05/2018] [Indexed: 12/29/2022]
Affiliation(s)
- Jun Kitano
- Division of Ecological Genetics; National Institute of Genetics; Mishima, Shizuoka Japan
| | - Asano Ishikawa
- Division of Ecological Genetics; National Institute of Genetics; Mishima, Shizuoka Japan
| | - Makoto Kusakabe
- Department of Biological Science; Faculty of Science; Shizuoka University; Surugaku, Shizuoka Japan
| |
Collapse
|