1
|
Segovia‐Ramírez MG, Ramírez‐Sánchez O, Decena Segarra LP, Rios‐Carlos H, Rovito SM. Determinants of genetic diversity in Neotropical salamanders (Plethodontidae: Bolitoglossini). Ecol Evol 2023; 13:e10707. [PMID: 38020701 PMCID: PMC10654480 DOI: 10.1002/ece3.10707] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2023] [Revised: 10/09/2023] [Accepted: 10/23/2023] [Indexed: 12/01/2023] Open
Abstract
Genetic diversity is the raw material of evolution, yet the reasons why it varies among species remain poorly understood. While studies at deeper phylogenetic scales point to the influence of life history traits on genetic diversity, it appears to be more affected by population size but less predictable at shallower scales. We used proxies for population size, mutation rate, direct selection, and linked selection to test factors affecting genetic diversity within a diverse assemblage of Neotropical salamanders, which vary widely for these traits. We estimated genetic diversity of noncoding loci using ddRADseq and coding loci using RNAseq for an assemblage of Neotropical salamanders distributed from northern Mexico to Costa Rica. Using ddRADseq loci, we found no significant association with genetic diversity, while for RNAseq data we found that environmental heterogeneity and proxies of population size predict a substantial portion of the variance in genetic diversity across species. Our results indicate that diversity of coding loci may be more predictable than that of noncoding loci, which appears to be mostly unpredictable at shallower phylogenetic scales. Our results suggest that coding loci may be more appropriate for genetic diversity estimates used in conservation planning because of the lack of any association between the variables we used and genetic diversity of noncoding loci.
Collapse
Affiliation(s)
| | - Obed Ramírez‐Sánchez
- Unidad de Genómica AvanzadaCentro de Investigación y de Estudios Avanzados del Instituto Politécnico NacionalIrapuatoMexico
| | - Louis Paul Decena Segarra
- Unidad de Genómica AvanzadaCentro de Investigación y de Estudios Avanzados del Instituto Politécnico NacionalIrapuatoMexico
| | - Hairo Rios‐Carlos
- Unidad de Genómica AvanzadaCentro de Investigación y de Estudios Avanzados del Instituto Politécnico NacionalIrapuatoMexico
| | - Sean M. Rovito
- Unidad de Genómica AvanzadaCentro de Investigación y de Estudios Avanzados del Instituto Politécnico NacionalIrapuatoMexico
| |
Collapse
|
2
|
Cunha RL, Nicastro KR, Zardi GI, Madeira C, McQuaid CD, Cox CJ, Castilho R. Comparative mitogenomic analyses and gene rearrangements reject the alleged polyphyly of a bivalve genus. PeerJ 2022; 10:e13953. [PMID: 36187748 PMCID: PMC9521344 DOI: 10.7717/peerj.13953] [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: 03/03/2022] [Accepted: 08/05/2022] [Indexed: 01/19/2023] Open
Abstract
Background The order and orientation of genes encoded by animal mitogenomes are typically conserved, although there is increasing evidence of multiple rearrangements among mollusks. The mitogenome from a Brazilian brown mussel (hereafter named B1) classified as Perna perna Linnaeus, 1758 and assembled from Illumina short-length reads revealed an unusual gene order very different from other congeneric species. Previous mitogenomic analyses based on the Brazilian specimen and other Mytilidae suggested the polyphyly of the genus Perna. Methods To confirm the proposed gene rearrangements, we sequenced a second Brazilian P. perna specimen using the "primer-walking" method and performed the assembly using as reference Perna canaliculus. This time-consuming sequencing method is highly effective when assessing gene order because it relies on sequentially-determined, overlapping fragments. We also sequenced the mitogenomes of eastern and southwestern South African P. perna lineages to analyze the existence of putative intraspecific gene order changes as the two lineages show overlapping distributions but do not exhibit a sister relationship. Results The three P. perna mitogenomes sequenced in this study exhibit the same gene order as the reference. CREx, a software that heuristically determines rearrangement scenarios, identified numerous gene order changes between B1 and our P. perna mitogenomes, rejecting the previously proposed gene order for the species. Our results validate the monophyly of the genus Perna and indicate a misidentification of B1.
Collapse
Affiliation(s)
- Regina L. Cunha
- Centre of Marine Sciences, CCMAR, University of Algarve, Campus de Gambelas, Faro, Algarve, Portugal
| | - Katy R. Nicastro
- Centre of Marine Sciences, CCMAR, University of Algarve, Campus de Gambelas, Faro, Algarve, Portugal,CNRS, Univ. Littoral Côte d’Opale, UMR 8187 – LOG – Laboratoire d’Océanologie et de Géosciences, Université de Lille, Lille, France,Department of Zoology and Entomology, Rhodes University, Grahamstown, South Africa
| | - Gerardo I. Zardi
- Department of Zoology and Entomology, Rhodes University, Grahamstown, South Africa
| | - Celine Madeira
- Centre of Marine Sciences, CCMAR, University of Algarve, Campus de Gambelas, Faro, Algarve, Portugal
| | | | - Cymon J. Cox
- Centre of Marine Sciences, CCMAR, University of Algarve, Campus de Gambelas, Faro, Algarve, Portugal
| | - Rita Castilho
- Centre of Marine Sciences, CCMAR, University of Algarve, Campus de Gambelas, Faro, Algarve, Portugal
| |
Collapse
|
3
|
Ivan J, Moritz C, Potter S, Bragg J, Turakulov R, Hua X. Temperature predicts the rate of molecular evolution in Australian Eugongylinae skinks. Evolution 2022; 76:252-261. [PMID: 34486736 DOI: 10.1111/evo.14342] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2021] [Revised: 07/12/2021] [Accepted: 08/15/2021] [Indexed: 01/21/2023]
Abstract
Temperature differences over time and space have been hypothesized to cause variation in the rate of molecular evolution of species, but empirical evidence is mixed. To further test this hypothesis, we utilized a large exon-capture sequence data of Australian Eugongylinae skinks, exemplifying a radiation of temperature-sensitive ectotherms spanning a large latitudinal gradient. The association between temperature (and other species traits) and long-term substitution rate was assessed based on 1268 sequenced exons of 44 species pairs from the Eugongylinae subfamily using regression analyses. Temperature is the strongest, positively correlated predictor of variation in substitution rate across the Australian Eugongylinae. It explains 45% of variation in synonymous substitution rate, and 11% after controlling for all the other factors. Synonymous substitution rate is also negatively associated with body size, with a 6% variation explained by body size after controlling for the effects of temperature. Other factors are not associated with synonymous substitution rate after controlling for temperature. Overall, this study points to temperature as a strong predictor of the molecular evolution rate in the Eugongylinae subfamily, and demonstrates the power of large-scale exonic data to identify correlates of the rate of molecular evolution.
Collapse
Affiliation(s)
- Jeremias Ivan
- Department of Bioinformatics, School of Life Sciences, Indonesia International Institute for Life Sciences, Jakarta, Indonesia
| | - Craig Moritz
- Division of Ecology and Evolution, Research School of Biology, Australian National University, Canberra, Australian Capital Territory, Australia
| | - Sally Potter
- Division of Ecology and Evolution, Research School of Biology, Australian National University, Canberra, Australian Capital Territory, Australia
| | - Jason Bragg
- Research Centre for Ecosystem Resilience, Australian Institute of Botanical Science, The Royal Botanic Garden Sydney, Sydney, New South Wales, Australia
| | - Rust Turakulov
- Laboratory of Pathology, National Cancer Institute, Bethesda, Maryland, United States
| | - Xia Hua
- Mathematical Sciences Institute, Australian National University, Canberra, Australian Capital Territory, Australia
| |
Collapse
|
4
|
Potter S, Bragg JG, Turakulov R, Eldridge MDB, Deakin J, Kirkpatrick M, Edwards RJ, Moritz C. Limited introgression between rock-wallabies with extensive chromosomal rearrangements. Mol Biol Evol 2021; 39:6448774. [PMID: 34865126 PMCID: PMC8788226 DOI: 10.1093/molbev/msab333] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Chromosome rearrangements can result in the rapid evolution of hybrid incompatibilities. Robertsonian fusions, particularly those with monobrachial homology, can drive reproductive isolation amongst recently diverged taxa. The recent radiation of rock-wallabies (genus Petrogale) is an important model to explore the role of Robertsonian fusions in speciation. Here, we pursue that goal using an extensive sampling of populations and genomes of Petrogale from north-eastern Australia. In contrast to previous assessments using mitochondrial DNA or nuclear microsatellite loci, genomic data are able to separate the most closely related species and to resolve their divergence histories. Both phylogenetic and population genetic analyses indicate introgression between two species that differ by a single Robertsonian fusion. Based on the available data, there is also evidence for introgression between two species which share complex chromosomal rearrangements. However, the remaining results show no consistent signature of introgression amongst species pairs and where evident, indicate generally low introgression overall. X-linked loci have elevated divergence compared with autosomal loci indicating a potential role for genic evolution to produce reproductive isolation in concert with chromosome change. Our results highlight the value of genome scale data in evaluating the role of Robertsonian fusions and structural variation in divergence, speciation, and patterns of molecular evolution.
Collapse
Affiliation(s)
- Sally Potter
- Division of Ecology and Evolution, Research School of Biology, The Australian National University, Acton, ACT, Australia.,Australian Museum Research Institute, Australian Museum, Sydney, NSW, Australia
| | - Jason G Bragg
- National Herbarium of New South Wales, The Royal Botanical Gardens and Domain Trust, Sydney, NSW, Australia
| | - Rustamzhon Turakulov
- Division of Ecology and Evolution, Research School of Biology, The Australian National University, Acton, ACT, Australia
| | - Mark D B Eldridge
- Australian Museum Research Institute, Australian Museum, Sydney, NSW, Australia
| | - Janine Deakin
- Institute for Applied Ecology, University of Canberra, Bruce, ACT, Australia
| | - Mark Kirkpatrick
- Department of Integrative Biology, University of Texas, Austin, TX, United States of America
| | - Richard J Edwards
- School of Biotechnology and Biomolecular Sciences, University of New South Wales, Sydney, NSW, 2052, Australia
| | - Craig Moritz
- Division of Ecology and Evolution, Research School of Biology, The Australian National University, Acton, ACT, Australia
| |
Collapse
|
5
|
Medina R, Wogan GOU, Bi K, Termignoni-García F, Bernal MH, Jaramillo-Correa JP, Wang IJ, Vázquez-Domínguez E. Phenotypic and genomic diversification with isolation by environment along elevational gradients in a neotropical treefrog. Mol Ecol 2021; 30:4062-4076. [PMID: 34160853 DOI: 10.1111/mec.16035] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2020] [Revised: 06/10/2021] [Accepted: 06/16/2021] [Indexed: 01/03/2023]
Abstract
Understanding how geographic and environmental heterogeneity drive local patterns of genetic variation is a major goal of ecological genomics and a key question in evolutionary biology. The tropical Andes and inter-Andean valleys are shaped by markedly heterogeneous landscapes, where species experience strong selective processes. We examined genome-wide SNP data together with behavioural and ecological traits (mating calls and body size) known to contribute to genetic isolation in anurans in the banana tree-dwelling frog, Boana platanera, distributed across an environmental gradient in Central Colombia (northern South America). Here, we analysed the relationships between environmentally (temperature and precipitation) associated genetic and phenotypic differentiation and the potential drivers of isolation by environment along an elevation gradient. We identified candidate SNPs associated with temperature and body size, which follow a clinal pattern of genome-wide differentiation tightly coupled with phenotypic variation: as elevation increases, B. platanera exhibits larger body size and longer call duration with more pulses but lower pulse rate and frequency. Thus, the environmental landscape has rendered a scenario where isolation by environment and candidate loci show concordance with phenotypic divergence in this tropical frog along an elevation gradient in the Colombian Andes. Our study sets the basis for evaluating the role of temperature in the genetic structure and local adaptation in tropical treefrogs and its putative effect on life cycle (embryos, tadpoles, adults) along elevation gradients.
Collapse
Affiliation(s)
- Ricardo Medina
- Departamento de Ecología de la Biodiversidad, Instituto de Ecología, Universidad Nacional Autónoma de México, Coyoacán, Ciudad de México, México.,Posgrado en Ciencias Biológicas, Universidad Nacional Autónoma de México, Coyoacán, Ciudad de México, México.,Grupo de Herpetología, Eco-Fisiología & Etología, Departamento de Biología, Universidad del Tolima, Altos de Santa Helena, Ibagué, Colombia
| | - Guinevere O U Wogan
- Department of Environmental Science, Policy, and Management, College of Natural Resources, University of California, Berkeley, California, USA.,Department of Integrative Biology, Oklahoma State University, Oklahoma, USA
| | - Ke Bi
- Museum of Vertebrate Zoology, University of California, Berkeley, California, USA.,Computational Genomics Resource Laboratory (CGRL, California Institute for Quantitative Biosciences (QB3, University of California, Berkeley, California, USA
| | - Flavia Termignoni-García
- Department of Organismic and Evolutionary Biology (OEB, Harvard University, Cambridge, Massachusetts, USA
| | - Manuel Hernando Bernal
- Grupo de Herpetología, Eco-Fisiología & Etología, Departamento de Biología, Universidad del Tolima, Altos de Santa Helena, Ibagué, Colombia
| | - Juan P Jaramillo-Correa
- Departamento de Ecología de la Biodiversidad, Instituto de Ecología, Universidad Nacional Autónoma de México, Coyoacán, Ciudad de México, México
| | - Ian J Wang
- Department of Environmental Science, Policy, and Management, College of Natural Resources, University of California, Berkeley, California, USA
| | - Ella Vázquez-Domínguez
- Departamento de Ecología de la Biodiversidad, Instituto de Ecología, Universidad Nacional Autónoma de México, Coyoacán, Ciudad de México, México
| |
Collapse
|
6
|
Ferris KG, Chavez AS, Suzuki TA, Beckman EJ, Phifer-Rixey M, Bi K, Nachman MW. The genomics of rapid climatic adaptation and parallel evolution in North American house mice. PLoS Genet 2021; 17:e1009495. [PMID: 33914747 PMCID: PMC8084166 DOI: 10.1371/journal.pgen.1009495] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2020] [Accepted: 03/17/2021] [Indexed: 12/23/2022] Open
Abstract
Parallel changes in genotype and phenotype in response to similar selection pressures in different populations provide compelling evidence of adaptation. House mice (Mus musculus domesticus) have recently colonized North America and are found in a wide range of environments. Here we measure phenotypic and genotypic differentiation among house mice from five populations sampled across 21° of latitude in western North America, and we compare our results to a parallel latitudinal cline in eastern North America. First, we show that mice are genetically differentiated between transects, indicating that they have independently colonized similar environments in eastern and western North America. Next, we find genetically-based differences in body weight and nest building behavior between mice from the ends of the western transect which mirror differences seen in the eastern transect, demonstrating parallel phenotypic change. We then conduct genome-wide scans for selection and a genome-wide association study to identify targets of selection and candidate genes for body weight. We find some genomic signatures that are unique to each transect, indicating population-specific responses to selection. However, there is significant overlap between genes under selection in eastern and western house mouse transects, providing evidence of parallel genetic evolution in response to similar selection pressures across North America.
Collapse
Affiliation(s)
- Kathleen G. Ferris
- Museum of Vertebrate Zoology and Department of Integrative Biology, University of California Berkeley, Berkeley, California, United States of America
| | - Andreas S. Chavez
- Museum of Vertebrate Zoology and Department of Integrative Biology, University of California Berkeley, Berkeley, California, United States of America
| | - Taichi A. Suzuki
- Museum of Vertebrate Zoology and Department of Integrative Biology, University of California Berkeley, Berkeley, California, United States of America
| | - Elizabeth J. Beckman
- Museum of Vertebrate Zoology and Department of Integrative Biology, University of California Berkeley, Berkeley, California, United States of America
| | - Megan Phifer-Rixey
- Museum of Vertebrate Zoology and Department of Integrative Biology, University of California Berkeley, Berkeley, California, United States of America
| | - Ke Bi
- 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
|
7
|
Byrne AQ, Richards-Zawacki CL, Voyles J, Bi K, Ibáñez R, Rosenblum EB. Whole exome sequencing identifies the potential for genetic rescue in iconic and critically endangered Panamanian harlequin frogs. GLOBAL CHANGE BIOLOGY 2021; 27:50-70. [PMID: 33150627 DOI: 10.1111/gcb.15405] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/18/2020] [Revised: 09/24/2020] [Accepted: 10/02/2020] [Indexed: 06/11/2023]
Abstract
Avoiding extinction in a rapidly changing environment often relies on a species' ability to quickly adapt in the face of extreme selective pressures. In Panamá, two closely related harlequin frog species (Atelopus varius and Atelopus zeteki) are threatened with extinction due to the fungal pathogen Batrachochytrium dendrobatidis (Bd). Once thought to be nearly extirpated from Panamá, A. varius have recently been rediscovered in multiple localities across their historical range; however, A. zeteki are possibly extinct in the wild. By leveraging a unique collection of 186 Atelopus tissue samples collected before and after the Bd outbreak in Panama, we describe the genetics of persistence for these species on the brink of extinction. We sequenced the transcriptome and developed an exome-capture assay to sequence the coding regions of the Atelopus genome. Using these genetic data, we evaluate the population genetic structure of historical A. varius and A. zeteki populations, describe changes in genetic diversity over time, assess the relationship between contemporary and historical individuals, and test the hypothesis that some A. varius populations have rapidly evolved to resist or tolerate Bd infection. We found a significant decrease in genetic diversity in contemporary (compared to historical) A. varius populations. We did not find strong evidence of directional allele frequency change or selection for Bd resistance genes, but we uncovered a set of candidate genes that warrant further study. Additionally, we found preliminary evidence of recent migration and gene flow in one of the largest persisting A. varius populations in Panamá, suggesting the potential for genetic rescue in this system. Finally, we propose that previous conservation units should be modified, as clear genetic breaks do not exist beyond the local population level. Our data lay the groundwork for genetically informed conservation and advance our understanding of how imperiled species might be rescued from extinction.
Collapse
Affiliation(s)
- Allison Q Byrne
- Department of Environmental Science, Policy, and Management, University of California Berkeley, Berkeley, CA, USA
- Museum of Vertebrate Zoology, University of California Berkeley, Berkeley, CA, USA
| | | | - Jamie Voyles
- Department of Biology, University of Nevada Reno, Reno, NV, USA
| | - Ke Bi
- Museum of Vertebrate Zoology, University of California Berkeley, Berkeley, CA, USA
| | - Roberto Ibáñez
- Smithsonian Tropical Research Institute, Panamá, República de Panamá
- Sistema Nacional de Investigación, SENACYT, Clayton, Panamá, República de Panamá
| | - Erica Bree Rosenblum
- Department of Environmental Science, Policy, and Management, University of California Berkeley, Berkeley, CA, USA
- Museum of Vertebrate Zoology, University of California Berkeley, Berkeley, CA, USA
| |
Collapse
|
8
|
Gwee CY, Garg KM, Chattopadhyay B, Sadanandan KR, Prawiradilaga DM, Irestedt M, Lei F, Bloch LM, Lee JGH, Irham M, Haryoko T, Soh MCK, Peh KSH, Rowe KMC, Ferasyi TR, Wu S, Wogan GOU, Bowie RCK, Rheindt FE. Phylogenomics of white-eyes, a 'great speciator', reveals Indonesian archipelago as the center of lineage diversity. eLife 2020; 9:e62765. [PMID: 33350381 PMCID: PMC7775107 DOI: 10.7554/elife.62765] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2020] [Accepted: 12/21/2020] [Indexed: 01/09/2023] Open
Abstract
Archipelagoes serve as important 'natural laboratories' which facilitate the study of island radiations and contribute to the understanding of evolutionary processes. The white-eye genus Zosterops is a classical example of a 'great speciator', comprising c. 100 species from across the Old World, most of them insular. We achieved an extensive geographic DNA sampling of Zosterops by using historical specimens and recently collected samples. Using over 700 genome-wide loci in conjunction with coalescent species tree methods and gene flow detection approaches, we untangled the reticulated evolutionary history of Zosterops, which comprises three main clades centered in Indo-Africa, Asia, and Australasia, respectively. Genetic introgression between species permeates the Zosterops phylogeny, regardless of how distantly related species are. Crucially, we identified the Indonesian archipelago, and specifically Borneo, as the major center of diversity and the only area where all three main clades overlap, attesting to the evolutionary importance of this region.
Collapse
Affiliation(s)
- Chyi Yin Gwee
- National University of Singapore, Department of Biological SciencesSingaporeSingapore
| | - Kritika M Garg
- National University of Singapore, Department of Biological SciencesSingaporeSingapore
| | - Balaji Chattopadhyay
- National University of Singapore, Department of Biological SciencesSingaporeSingapore
| | - Keren R Sadanandan
- National University of Singapore, Department of Biological SciencesSingaporeSingapore
- Max Planck Institute for OrnithologySeewiesenGermany
| | - Dewi M Prawiradilaga
- Division of Zoology, Research Center for Biology, Indonesian Institute of Sciences (LIPI), Cibinong Science CenterCibinongIndonesia
| | - Martin Irestedt
- Department of Bioinformatics and Genetics, Swedish Museum of Natural HistoryStockholmSweden
| | - Fumin Lei
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of SciencesBeijingChina
- Center for Excellence in Animal Evolution and Genetics, Chinese Academy of SciencesKunmingChina
| | - Luke M Bloch
- Museum of Vertebrate Zoology and Department of Integrative Biology, University of California, BerkeleyBerkeleyUnited States
| | | | - Mohammad Irham
- Division of Zoology, Research Center for Biology, Indonesian Institute of Sciences (LIPI), Cibinong Science CenterCibinongIndonesia
| | - Tri Haryoko
- Division of Zoology, Research Center for Biology, Indonesian Institute of Sciences (LIPI), Cibinong Science CenterCibinongIndonesia
| | - Malcolm CK Soh
- University of Western Australia, School of Biological SciencesPerthAustralia
| | - Kelvin S-H Peh
- University of Southampton, School of Biological Sciences, UniversitySouthamptonUnited Kingdom
| | - Karen MC Rowe
- Sciences Department, Museums VictoriaMelbourneAustralia
| | - Teuku Reza Ferasyi
- Faculty of Veterinary Medicine, Universitas Syiah KualaDarussalamIndonesia
- Jiangsu Key Laboratory of Phylogenomics and Comparative Genomics, School of Life Sciences, Jiangsu Normal UniversityXuzhouChina
| | - Shaoyuan Wu
- Department of Biochemistry and Molecular Biology, 2011 Collaborative Innovation Center of Tianjin for Medical Epigenetics, Tianjin Key Laboratory of Medical Epigenetics, School of Basic Medical Sciences, Tianjin Medical UniversityTianjinChina
- Center for Tropical Veterinary Studies – One Health Collaboration Center, Universitas Syiah KualaDarussalamIndonesia
| | - Guinevere OU Wogan
- Museum of Vertebrate Zoology and Department of Environmental Science, Policy, and Management, University of California, BerkeleyBerkeleyUnited States
| | - Rauri CK Bowie
- Museum of Vertebrate Zoology and Department of Integrative Biology, University of California, BerkeleyBerkeleyUnited States
| | - Frank E Rheindt
- National University of Singapore, Department of Biological SciencesSingaporeSingapore
| |
Collapse
|
9
|
Meléndez-Rosa J, Bi K, Lacey EA. Mating system is correlated with immunogenetic diversity in sympatric species of Peromyscine mice. PLoS One 2020; 15:e0236084. [PMID: 32701975 PMCID: PMC7377423 DOI: 10.1371/journal.pone.0236084] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2019] [Accepted: 06/29/2020] [Indexed: 11/18/2022] Open
Abstract
The number of reproductive partners per individual varies markedly across animal mating systems. This variation may be an important determinant of patterns of immunogenetic diversity, particularly at Major Histocompatibility Complex (MHC) Class I and II loci. To compare immunogenetic variation in taxa with markedly different mating systems, we used RNAseq-generated data to quantify genotypic diversity in three species of Peromyscine rodents: the monogamous California mouse (Peromyscus californicus) and the polygynandrous deer mouse (P. maniculatus) and brush mouse (P. boylii). By sampling populations of these species from multiple localities in California, we were able to conduct replicated analyses of the relationship between mating system and immunogenetic variation. Across the four localities sampled, diversity at MHC Class I and II genes was consistently higher in the two polygynandrous species. We found no evidence that sampling location (i.e., variation in habitat conditions) contributed to observed differences in MHC variation among populations or species. Collectively, our data indicate that immunogenetic variation in Peromyscine mice is associated with reproductive behavior, rather than geographic locality or habitat type. The consistently greater variability detected in the polygynandrous species examined suggests that balancing selection imposed by behaviorally-mediated pathogen exposure is important in maintaining variation at MHC genes in these animals.
Collapse
Affiliation(s)
- Jesyka Meléndez-Rosa
- Department of Integrative Biology, University of California, Berkeley, California, United States of America
- Museum of Vertebrate Zoology, University of California, Berkeley, California, United States of America
| | - Ke Bi
- Museum of Vertebrate Zoology, University of California, Berkeley, California, United States of America
- Computational Genomics Resource Laboratory, University of California, Berkeley, California, United States of America
| | - Eileen A. Lacey
- Department of Integrative Biology, University of California, Berkeley, California, United States of America
- Museum of Vertebrate Zoology, University of California, Berkeley, California, United States of America
| |
Collapse
|
10
|
Hanf ZR, Chavez AS. A Comprehensive Multi-Omic Approach Reveals a Relatively Simple Venom in a Diet Generalist, the Northern Short-Tailed Shrew, Blarina brevicauda. Genome Biol Evol 2020; 12:1148-1166. [PMID: 32520994 PMCID: PMC7486961 DOI: 10.1093/gbe/evaa115] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 06/05/2020] [Indexed: 12/15/2022] Open
Abstract
Animals that use venom to feed on a wide diversity of prey may evolve a complex mixture of toxins to target a variety of physiological processes and prey-defense mechanisms. Blarina brevicauda, the northern short-tailed shrew, is one of few venomous mammals, and is also known to eat evolutionarily divergent prey. Despite their complex diet, earlier proteomic and transcriptomic studies of this shrew's venom have only identified two venom proteins. Here, we investigated with comprehensive molecular approaches whether B. brevicauda venom is more complex than previously understood. We generated de novo assemblies of a B. brevicauda genome and submaxillary-gland transcriptome, as well as sequenced the salivary proteome. Our findings show that B. brevicauda's venom composition is simple relative to their broad diet and is likely limited to seven proteins from six gene families. Additionally, we explored expression levels and rate of evolution of these venom genes and the origins of key duplications that led to toxin neofunctionalization. We also found three proteins that may be involved in endogenous self-defense. The possible synergism of the toxins suggests that vertebrate prey may be the main target of the venom. Further functional assays for all venom proteins on both vertebrate and invertebrate prey would provide further insight into the ecological relevance of venom in this species.
Collapse
Affiliation(s)
- Zachery R Hanf
- Department of Evolution, Ecology, and Organismal Biology, The Ohio State University
| | - Andreas S Chavez
- Department of Evolution, Ecology, and Organismal Biology, The Ohio State University
- Translational Data Analytics Institute, The Ohio State University
| |
Collapse
|
11
|
Denney DA, Jameel MI, Bemmels JB, Rochford ME, Anderson JT. Small spaces, big impacts: contributions of micro-environmental variation to population persistence under climate change. AOB PLANTS 2020; 12:plaa005. [PMID: 32211145 PMCID: PMC7082537 DOI: 10.1093/aobpla/plaa005] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/28/2019] [Accepted: 02/06/2020] [Indexed: 05/05/2023]
Abstract
Individuals within natural populations can experience very different abiotic and biotic conditions across small spatial scales owing to microtopography and other micro-environmental gradients. Ecological and evolutionary studies often ignore the effects of micro-environment on plant population and community dynamics. Here, we explore the extent to which fine-grained variation in abiotic and biotic conditions contributes to within-population variation in trait expression and genetic diversity in natural plant populations. Furthermore, we consider whether benign microhabitats could buffer local populations of some plant species from abiotic stresses imposed by rapid anthropogenic climate change. If microrefugia sustain local populations and communities in the short term, other eco-evolutionary processes, such as gene flow and adaptation, could enhance population stability in the longer term. We caution, however, that local populations may still decline in size as they contract into rare microhabitats and microrefugia. We encourage future research that explicitly examines the role of the micro-environment in maintaining genetic variation within local populations, favouring the evolution of phenotypic plasticity at local scales and enhancing population persistence under global change.
Collapse
Affiliation(s)
- Derek A Denney
- Department of Plant Biology, University of Georgia, Athens, GA, USA
| | - M Inam Jameel
- Department of Genetics, University of Georgia, Athens, GA, USA
| | - Jordan B Bemmels
- Department of Genetics, University of Georgia, Athens, GA, USA
- Department of Biological Sciences, University of Toronto Scarborough, Toronto, ON, Canada
| | - Mia E Rochford
- Department of Plant Biology, University of Georgia, Athens, GA, USA
| | - Jill T Anderson
- Department of Genetics, University of Georgia, Athens, GA, USA
| |
Collapse
|
12
|
Ceschin DG, Pires NS, Mardirosian MN, Lascano CI, Venturino A. The Rhinella arenarum transcriptome: de novo assembly, annotation and gene prediction. Sci Rep 2020; 10:1053. [PMID: 31974515 PMCID: PMC6978513 DOI: 10.1038/s41598-020-57961-4] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2019] [Accepted: 01/09/2020] [Indexed: 12/13/2022] Open
Abstract
The common toad Rhinella arenarum is widely distributed in Argentina, where it is utilised as an autochthonous model in ecotoxicological research and environmental toxicology. However, the lack of a reference genome makes molecular assays and gene expression studies difficult to carry out on this non-model species. To address this issue, we performed a genome-wide transcriptome analysis on R. arenarum larvae through massive RNA sequencing, followed by de novo assembly, annotation, and gene prediction. We obtained 57,407 well-annotated transcripts representing 99.4% of transcriptome completeness (available at http://rhinella.uncoma.edu.ar). We also defined a set of 52,800 high-confidence lncRNA transcripts and demonstrated the reliability of the transcriptome data to perform phylogenetic analysis. Our comprehensive transcriptome analysis of R. arenarum represents a valuable resource to perform functional genomic studies and to identify potential molecular biomarkers in ecotoxicological research.
Collapse
Affiliation(s)
- Danilo Guillermo Ceschin
- Centro de Investigaciones en Toxicología Ambiental y Agrobiotecnología del Comahue (CITAAC), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET)-Universidad Nacional del Comahue, Buenos Aires 1400, Neuquén, CP 8300, Neuquén, Argentina.
- Laboratorio de Bioinformática Traslacional, Centro de Investigaciones en Medicina Traslacional Severo Amuchástegui, Instituto Universitario de Ciencias Biomédicas de Córdoba. Av. Naciones Unidas 420, CP 5000, Córdoba, Argentina.
| | - Natalia Susana Pires
- Centro de Investigaciones en Toxicología Ambiental y Agrobiotecnología del Comahue (CITAAC), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET)-Universidad Nacional del Comahue, Buenos Aires 1400, Neuquén, CP 8300, Neuquén, Argentina
| | - Mariana Noelia Mardirosian
- Centro de Investigaciones en Toxicología Ambiental y Agrobiotecnología del Comahue (CITAAC), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET)-Universidad Nacional del Comahue, Buenos Aires 1400, Neuquén, CP 8300, Neuquén, Argentina
| | - Cecilia Inés Lascano
- Centro de Investigaciones en Toxicología Ambiental y Agrobiotecnología del Comahue (CITAAC), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET)-Universidad Nacional del Comahue, Buenos Aires 1400, Neuquén, CP 8300, Neuquén, Argentina
| | - Andrés Venturino
- Centro de Investigaciones en Toxicología Ambiental y Agrobiotecnología del Comahue (CITAAC), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET)-Universidad Nacional del Comahue, Buenos Aires 1400, Neuquén, CP 8300, Neuquén, Argentina
| |
Collapse
|
13
|
Portik DM, Bell RC, Blackburn DC, Bauer AM, Barratt CD, Branch WR, Burger M, Channing A, Colston TJ, Conradie W, Dehling JM, Drewes RC, Ernst R, Greenbaum E, Gvoždík V, Harvey J, Hillers A, Hirschfeld M, Jongsma GFM, Kielgast J, Kouete MT, Lawson LP, Leaché AD, Loader SP, Lötters S, Meijden AVD, Menegon M, Müller S, Nagy ZT, Ofori-Boateng C, Ohler A, Papenfuss TJ, Rößler D, Sinsch U, Rödel MO, Veith M, Vindum J, Zassi-Boulou AG, McGuire JA. Sexual Dichromatism Drives Diversification within a Major Radiation of African Amphibians. Syst Biol 2020; 68:859-875. [PMID: 31140573 DOI: 10.1093/sysbio/syz023] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2018] [Revised: 02/15/2019] [Accepted: 04/09/2019] [Indexed: 01/11/2023] Open
Abstract
Theory predicts that sexually dimorphic traits under strong sexual selection, particularly those involved with intersexual signaling, can accelerate speciation and produce bursts of diversification. Sexual dichromatism (sexual dimorphism in color) is widely used as a proxy for sexual selection and is associated with rapid diversification in several animal groups, yet studies using phylogenetic comparative methods to explicitly test for an association between sexual dichromatism and diversification have produced conflicting results. Sexual dichromatism is rare in frogs, but it is both striking and prevalent in African reed frogs, a major component of the diverse frog radiation termed Afrobatrachia. In contrast to most other vertebrates, reed frogs display female-biased dichromatism in which females undergo color transformation, often resulting in more ornate coloration in females than in males. We produce a robust phylogeny of Afrobatrachia to investigate the evolutionary origins of sexual dichromatism in this radiation and examine whether the presence of dichromatism is associated with increased rates of net diversification. We find that sexual dichromatism evolved once within hyperoliids and was followed by numerous independent reversals to monochromatism. We detect significant diversification rate heterogeneity in Afrobatrachia and find that sexually dichromatic lineages have double the average net diversification rate of monochromatic lineages. By conducting trait simulations on our empirical phylogeny, we demonstrate that our inference of trait-dependent diversification is robust. Although sexual dichromatism in hyperoliid frogs is linked to their rapid diversification and supports macroevolutionary predictions of speciation by sexual selection, the function of dichromatism in reed frogs remains unclear. We propose that reed frogs are a compelling system for studying the roles of natural and sexual selection on the evolution of sexual dichromatism across micro- and macroevolutionary timescales.
Collapse
Affiliation(s)
- Daniel M Portik
- Museum of Vertebrate Zoology, University of California, Berkeley, CA 94720, USA.,Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, AZ 85721, USA
| | - Rayna C Bell
- Museum of Vertebrate Zoology, University of California, Berkeley, CA 94720, USA.,Department of Vertebrate Zoology, National Museum of Natural History, Smithsonian Institution, Washington, DC 20560-0162, USA
| | - David C Blackburn
- Florida Museum of Natural History, University of Florida, Gainesville, FL 32611, USA
| | - Aaron M Bauer
- Department of Biology, Villanova University, 800 Lancaster Avenue, Villanova, PA 19085, USA
| | - Christopher D Barratt
- Department of Environmental Sciences, University of Basel, Basel 4056, Switzerland.,German Centre for Integrative Biodiversity Research (iDiv) Halle-Jena-Leipzig, Leipzig 0413, Germany.,Max Planck Institute for Evolutionary Anthropology, Leipzig 0413, Germany
| | - William R Branch
- Port Elizabeth Museum, P.O. Box 11347, Humewood 6013, South Africa.,Department of Zoology, Nelson Mandela Metropolitan University, P.O. Box 77000, Port Elizabeth 6031, South Africa
| | - Marius Burger
- African Amphibian Conservation Research Group, Unit for Environmental Sciences and Management, North-West University, Potchefstroom 2520, South Africa.,Flora Fauna & Man, Ecological Services Ltd. Tortola, British Virgin, Island
| | - Alan Channing
- Unit for Environmental Sciences and Management, North-West University, Potchefstroom 2520, South Africa
| | - Timothy J Colston
- Department of Biological Sciences, Florida State University, Tallahassee, FL 32306, USA.,Zoological Natural History Museum, Addis Ababa University, Arat Kilo, Addis Ababa, Ethiopia
| | - Werner Conradie
- Port Elizabeth Museum, P.O. Box 11347, Humewood 6013, South Africa.,School of Natural Resource Management, Nelson Mandela University, George Campus, George 6530, South Africa
| | - J Maximilian Dehling
- Department of Biology, Institute of Sciences, University of Koblenz-Landau, Universitätsstr. 1, D-56070 Koblenz, Germany
| | - Robert C Drewes
- California Academy of Sciences, San Francisco, CA 94118, USA
| | - Raffael Ernst
- Museum of Zoology, Senckenberg Natural History Collections Dresden, Königsbrücker Landstr. 159, Dresden 01109, Germany.,Department of Ecology, Technische Universität Berlin, Rothenburgstr. 12, Berlin 12165, Germany
| | - Eli Greenbaum
- Department of Biological Sciences, University of Texas at El Paso, El Paso, TX 79968, USA
| | - Václav Gvoždík
- The Czech Academy of Sciences, Institute of Vertebrate Biology, Brno, Czech Republic.,Department of Zoology, National Museum, Prague, Czech Republic
| | | | - Annika Hillers
- Museum für Naturkunde, Leibniz Institute for Evolution and Biodiversity Science, Biodiversity Dynamics, Invalidenstr. 43, Berlin 10115, Germany.,Across the River - A Transboundary Peace Park for Sierra Leone and Liberia, The Royal Society for the Protection of Birds, 164 Dama Road, Kenema, Sierra Leone
| | - Mareike Hirschfeld
- Museum für Naturkunde, Leibniz Institute for Evolution and Biodiversity Science, Biodiversity Dynamics, Invalidenstr. 43, Berlin 10115, Germany
| | - Gregory F M Jongsma
- Florida Museum of Natural History, University of Florida, Gainesville, FL 32611, USA
| | - Jos Kielgast
- Natural History Museum of Denmark, University of Copenhagen, Universitetsparken 15, Copenhagen 2100, Denmark
| | - Marcel T Kouete
- Florida Museum of Natural History, University of Florida, Gainesville, FL 32611, USA
| | - Lucinda P Lawson
- Department of Biological Sciences, University of Cincinnati, 614 Rieveschl Hall, Cincinnati, OH 45220, USA.,Life Sciences, Field Museum of Natural History, 1400 S. Lake Shore Dr., Chicago, IL 60605, USA
| | - Adam D Leaché
- Department of Biology, Burke Museum of Natural History and Culture, University of Washington, Seattle, WA, USA
| | - Simon P Loader
- Life Sciences Department, Natural History Museum, London SW7 5BD, UK
| | - Stefan Lötters
- Biogeography Department, Trier University, Universitätsring 15, Trier 54296, Germany
| | - Arie Van Der Meijden
- CIBIO Research Centre in Biodiversity and Genetic Resources, InBIO, Universidade do Porto, Campus Agrario de Vairão, Rua Padre Armando Quintas, No. 7, 4485-661 Vairão, Vila do Conde, Portugal
| | - Michele Menegon
- Tropical Biodiversity Section, Science Museum of Trento, Corso del lavoro e della Scienza 3, Trento 38122, Italy
| | - Susanne Müller
- Biogeography Department, Trier University, Universitätsring 15, Trier 54296, Germany
| | - Zoltán T Nagy
- Royal Belgian Institute of Natural Sciences, OD Taxonomy and Phylogeny, Rue Vautier 29, B-1000 Brussels, Belgium
| | | | - Annemarie Ohler
- Département Origines et Evolution, Muséum National d'Histoire Naturelle, UMR 7205 ISYEB, 25 rue Cuvier, Paris 75005, France
| | | | - Daniela Rößler
- Biogeography Department, Trier University, Universitätsring 15, Trier 54296, Germany
| | - Ulrich Sinsch
- Department of Biology, Institute of Sciences, University of Koblenz-Landau, Universitätsstr. 1, D-56070 Koblenz, Germany
| | - Mark-Oliver Rödel
- Museum für Naturkunde, Leibniz Institute for Evolution and Biodiversity Science, Biodiversity Dynamics, Invalidenstr. 43, Berlin 10115, Germany
| | - Michael Veith
- Biogeography Department, Trier University, Universitätsring 15, Trier 54296, Germany
| | - Jens Vindum
- California Academy of Sciences, San Francisco, CA 94118, USA
| | - Ange-Ghislain Zassi-Boulou
- Institut National de Recherche en Sciences Exactes et Naturelles, Brazzaville BP 2400, République du Congo
| | - Jimmy A McGuire
- Museum of Vertebrate Zoology, University of California, Berkeley, CA 94720, USA
| |
Collapse
|
14
|
Wogan GOU, Yuan ML, Mahler DL, Wang IJ. Genome-wide epigenetic isolation by environment in a widespread Anolis lizard. Mol Ecol 2019; 29:40-55. [PMID: 31710739 DOI: 10.1111/mec.15301] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2019] [Revised: 10/26/2019] [Accepted: 11/08/2019] [Indexed: 12/12/2022]
Abstract
Epigenetic changes can provide a pathway for organisms to respond to local environmental conditions by influencing gene expression. However, we still know little about the spatial distribution of epigenetic variation in natural systems, how it relates to the distribution of genetic variation and the environmental structure of the landscape, and the processes that generate and maintain it. Studies examining spatial patterns of genetic and epigenetic variation can provide valuable insights into how ecological and population processes contribute to epigenetic divergence across heterogeneous landscapes. Here, we perform a comparative analysis of spatial genetic and epigenetic variation based on 8,459 single nucleotide polymorphisms (SNPs) and 8,580 single methylation variants (SMVs) from eight populations of the Puerto Rican crested anole, Anolis cristatellus, an abundant lizard in the adaptive radiations of anoles on the Greater Antilles that occupies a diverse range of habitats. Using generalized dissimilarity modelling and multiple matrix regression, we found that genome-wide epigenetic differentiation is strongly correlated with environmental divergence, even after controlling for the underlying genetic structure. We also detected significant associations between key environmental variables and 96 SMVs, including 42 located in promoter regions or gene bodies. Our results suggest an environmental basis for population-level epigenetic differentiation in this system and contribute to better understanding how environmental gradients structure epigenetic variation in nature.
Collapse
Affiliation(s)
- Guinevere O U Wogan
- Department of Environmental Science, Policy, and Management, College of Natural Resources, University of California, Berkeley, CA, USA
| | - Michael L Yuan
- Department of Environmental Science, Policy, and Management, College of Natural Resources, University of California, Berkeley, CA, USA
| | - D Luke Mahler
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, ON, Canada
| | - Ian J Wang
- Department of Environmental Science, Policy, and Management, College of Natural Resources, University of California, Berkeley, CA, USA
| |
Collapse
|
15
|
Lim MCW, Witt CC, Graham CH, Dávalos LM. Parallel Molecular Evolution in Pathways, Genes, and Sites in High-Elevation Hummingbirds Revealed by Comparative Transcriptomics. Genome Biol Evol 2019; 11:1552-1572. [PMID: 31028697 PMCID: PMC6553502 DOI: 10.1093/gbe/evz101] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/12/2019] [Indexed: 12/13/2022] Open
Abstract
High-elevation organisms experience shared environmental challenges that include low oxygen availability, cold temperatures, and intense ultraviolet radiation. Consequently, repeated evolution of the same genetic mechanisms may occur across high-elevation taxa. To test this prediction, we investigated the extent to which the same biochemical pathways, genes, or sites were subject to parallel molecular evolution for 12 Andean hummingbird species (family: Trochilidae) representing several independent transitions to high elevation across the phylogeny. Across high-elevation species, we discovered parallel evolution for several pathways and genes with evidence of positive selection. In particular, positively selected genes were frequently part of cellular respiration, metabolism, or cell death pathways. To further examine the role of elevation in our analyses, we compared results for low- and high-elevation species and tested different thresholds for defining elevation categories. In analyses with different elevation thresholds, positively selected genes reflected similar functions and pathways, even though there were almost no specific genes in common. For example, EPAS1 (HIF2α), which has been implicated in high-elevation adaptation in other vertebrates, shows a signature of positive selection when high-elevation is defined broadly (>1,500 m), but not when defined narrowly (>2,500 m). Although a few biochemical pathways and genes change predictably as part of hummingbird adaptation to high-elevation conditions, independent lineages have rarely adapted via the same substitutions.
Collapse
Affiliation(s)
- Marisa C W Lim
- Department of Ecology and Evolution, Stony Brook University
| | - Christopher C Witt
- Museum of Southwestern Biology and Department of Biology, University of New Mexico
| | - Catherine H Graham
- Department of Ecology and Evolution, Stony Brook University.,Swiss Federal Research Institute (WSL), Birmensdorf, Switzerland
| | - Liliana M Dávalos
- Department of Ecology and Evolution, Stony Brook University.,Consortium for Inter-Disciplinary Environmental Research, Stony Brook University
| |
Collapse
|
16
|
Skipwith PL, Bi K, Oliver PM. Relicts and radiations: Phylogenomics of an Australasian lizard clade with east Gondwanan origins (Gekkota: Diplodactyloidea). Mol Phylogenet Evol 2019; 140:106589. [DOI: 10.1016/j.ympev.2019.106589] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2019] [Revised: 08/06/2019] [Accepted: 08/12/2019] [Indexed: 10/26/2022]
|
17
|
Goodman KR, Prost S, Bi K, Brewer MS, Gillespie RG. Host and geography together drive early adaptive radiation of Hawaiian planthoppers. Mol Ecol 2019; 28:4513-4528. [PMID: 31484218 DOI: 10.1111/mec.15231] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2018] [Revised: 08/19/2019] [Accepted: 08/27/2019] [Indexed: 11/30/2022]
Abstract
The interactions between insects and their plant host have been implicated in driving diversification of both players. Early arguments highlighted the role of ecological opportunity, with the idea that insects "escape and radiate" on new hosts, with subsequent hypotheses focusing on the interplay between host shifting and host tracking, coupled with isolation and fusion, in generating diversity. Because it is rarely possible to capture the initial stages of diversification, it is particularly difficult to ascertain the relative roles of geographic isolation versus host shifts in initiating the process. The current study examines genetic diversity between populations and hosts within a single species of endemic Hawaiian planthopper, Nesosydne umbratica (Hemiptera, Delphacidae). Given that the species was known as a host generalist occupying unrelated hosts, Clermontia (Campanulaceae) and Pipturus (Urticaceae), we set out to determine the relative importance of geography and host in structuring populations in the early stages of differentiation on the youngest islands of the Hawaiian chain. Results from extensive exon capture data showed that N. umbratica is highly structured, both by geography, with discrete populations on each volcano, and by host plant, with parallel radiations on Clermontia and Pipturus leading to extensive co-occurrence. The marked genetic structure suggests that populations can readily become established on novel hosts provided opportunity; subsequent adaptation allows monopolization of the new host. The results support the role of geographic isolation in structuring populations and with host shifts occurring as discrete events that facilitate subsequent parallel geographic range expansion.
Collapse
Affiliation(s)
- Kari Roesch Goodman
- Department of Environmental Science, Policy, and Management, University of California, Berkeley, CA, USA
| | - Stefan Prost
- Department of Integrative Biology, University of California, Berkeley, CA, USA.,LOEWE-Centre for Translational Biodiversity Genomics, Senckenberg Research Institute, Frankfurt/Main, Germany
| | - Ke Bi
- Computational Genomics Resource Laboratory (CGRL), California Institute for Quantitative Biosciences (QB3), University of California, Berkeley, CA, USA.,Ancestry, San Francisco, CA, USA.,Museum of Vertebrate Zoology, University of California, Berkeley, CA, USA
| | - Michael S Brewer
- Department of Biology, East Carolina University, Greenville, NC, USA
| | - Rosemary G Gillespie
- Department of Environmental Science, Policy, and Management, University of California, Berkeley, CA, USA
| |
Collapse
|
18
|
Blom MPK, Matzke NJ, Bragg JG, Arida E, Austin CC, Backlin AR, Carretero MA, Fisher RN, Glaw F, Hathaway SA, Iskandar DT, McGuire JA, Karin BR, Reilly SB, Rittmeyer EN, Rocha S, Sanchez M, Stubbs AL, Vences M, Moritz C. Habitat preference modulates trans-oceanic dispersal in a terrestrial vertebrate. Proc Biol Sci 2019; 286:20182575. [PMID: 31161911 DOI: 10.1098/rspb.2018.2575] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
The importance of long-distance dispersal (LDD) in shaping geographical distributions has been debated since the nineteenth century. In terrestrial vertebrates, LDD events across large water bodies are considered highly improbable, but organismal traits affecting dispersal capacity are generally not taken into account. Here, we focus on a recent lizard radiation and combine a summary-coalescent species tree based on 1225 exons with a probabilistic model that links dispersal capacity to an evolving trait, to investigate whether ecological specialization has influenced the probability of trans-oceanic dispersal. Cryptoblepharus species that occur in coastal habitats have on average dispersed 13 to 14 times more frequently than non-coastal species and coastal specialization has, therefore, led to an extraordinarily widespread distribution that includes multiple continents and distant island archipelagoes. Furthermore, their presence across the Pacific substantially predates the age of human colonization and we can explicitly reject the possibility that these patterns are solely shaped by human-mediated dispersal. Overall, by combining new analytical methods with a comprehensive phylogenomic dataset, we use a quantitative framework to show how coastal specialization can influence dispersal capacity and eventually shape geographical distributions at a macroevolutionary scale.
Collapse
Affiliation(s)
- Mozes P K Blom
- 1 Research School of Biology, The Australian National University , Canberra , Australia.,2 Museum für Naturkunde, Leibniz Institut für Evolutions- und Biodiversitätsforschung , Berlin , Germany
| | - Nicholas J Matzke
- 1 Research School of Biology, The Australian National University , Canberra , Australia.,3 School of Biological Sciences, University of Auckland , Auckland , New Zealand
| | - Jason G Bragg
- 1 Research School of Biology, The Australian National University , Canberra , Australia
| | - Evy Arida
- 4 Research Center for Biology, The Indonesian Institute of Sciences , Cibinong , Indonesia
| | | | - Adam R Backlin
- 6 U.S. Geological Survey, Western Ecological Research Center , Santa Ana, CA , USA
| | | | - Robert N Fisher
- 8 U.S. Geological Survey, Western Ecological Research Center , San Diego, CA , USA
| | - Frank Glaw
- 9 Department of Herpetology, Zoologische Staatssamlung Münich , Munich , Germany
| | - Stacie A Hathaway
- 8 U.S. Geological Survey, Western Ecological Research Center , San Diego, CA , USA
| | - Djoko T Iskandar
- 10 School of Life Sciences and Technology, Institut Teknologi , Bandung , Indonesia
| | - Jimmy A McGuire
- 11 Museum of Vertebrate Zoology and Department of Integrative Biology, University of California Berkeley , Berkeley, CA , USA
| | - Benjamin R Karin
- 11 Museum of Vertebrate Zoology and Department of Integrative Biology, University of California Berkeley , Berkeley, CA , USA
| | - Sean B Reilly
- 11 Museum of Vertebrate Zoology and Department of Integrative Biology, University of California Berkeley , Berkeley, CA , USA
| | - Eric N Rittmeyer
- 1 Research School of Biology, The Australian National University , Canberra , Australia.,5 Museum of Natural Science, Louisiana State University , Baton Rouge, LA , USA
| | - Sara Rocha
- 12 Department of Biochemistry, Genetics and Immunology & Biomedical Research Center (CINBIO), University of Vigo , Vigo , Spain
| | | | - Alexander L Stubbs
- 11 Museum of Vertebrate Zoology and Department of Integrative Biology, University of California Berkeley , Berkeley, CA , USA
| | - Miguel Vences
- 14 Zoological Institute, Technische Universität Braunschweig , Braunschweig , Germany
| | - Craig Moritz
- 1 Research School of Biology, The Australian National University , Canberra , Australia
| |
Collapse
|
19
|
Lim MCW, Witt CC, Graham CH, Dávalos LM. Divergent Fine-Scale Recombination Landscapes between a Freshwater and Marine Population of Threespine Stickleback Fish. Genome Biol Evol 2019; 11:1573-1585. [PMID: 31028697 PMCID: PMC6553502 DOI: 10.1093/gbe/evz090] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 04/17/2019] [Indexed: 12/27/2022] Open
Abstract
Meiotic recombination is a highly conserved process that has profound effects on genome evolution. At a fine-scale, recombination rates can vary drastically across genomes, often localized into small recombination "hotspots" with highly elevated rates, surrounded by regions with little recombination. In most species studied, the location of hotspots within genomes is highly conserved across broad evolutionary timescales. The main exception to this pattern is in mammals, where hotspot location can evolve rapidly among closely related species and even among populations within a species. Hotspot position in mammals is controlled by the gene, Prdm9, whereas in species with conserved hotspots, a functional Prdm9 is typically absent. Due to a limited number of species where recombination rates have been estimated at a fine-scale, it remains unclear whether hotspot conservation is always associated with the absence of a functional Prdm9. Threespine stickleback fish (Gasterosteus aculeatus) are an excellent model to examine the evolution of recombination over short evolutionary timescales. Using a linkage disequilibrium-based approach, we found recombination rates indeed varied at a fine-scale across the genome, with many regions organized into narrow hotspots. Hotspots had highly divergent landscapes between stickleback populations, where only ∼15% of these hotspots were shared. Our results indicate that fine-scale recombination rates may be diverging between closely related populations of threespine stickleback fish. Interestingly, we found only a weak association of a PRDM9 binding motif within hotspots, which suggests that threespine stickleback fish may possess a novel mechanism for targeting recombination hotspots at a fine-scale.
Collapse
Affiliation(s)
- Marisa C W Lim
- Department of Ecology and Evolution, Stony Brook University
| | - Christopher C Witt
- Museum of Southwestern Biology and Department of Biology, University of New Mexico
| | - Catherine H Graham
- Department of Ecology and Evolution, Stony Brook University
- Swiss Federal Research Institute (WSL), Birmensdorf, Switzerland
| | - Liliana M Dávalos
- Department of Ecology and Evolution, Stony Brook University
- Consortium for Inter-Disciplinary Environmental Research, Stony Brook University
| |
Collapse
|
20
|
Meléndez‐Rosa J, Bi K, Lacey EA. Differential gene expression in relation to mating system in Peromyscine rodents. Ecol Evol 2019; 9:5975-5990. [PMID: 31161013 PMCID: PMC6540711 DOI: 10.1002/ece3.5181] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2019] [Revised: 03/15/2019] [Accepted: 03/29/2019] [Indexed: 12/25/2022] Open
Abstract
Behaviors that increase an individual's exposure to pathogens are expected to have important effects on immunoactivity. Because sexual reproduction typically requires close contact among conspecifics, mating systems provide an ideal opportunity to study the immunogenetic correlates of behaviors with high versus low risks of pathogen exposure. Despite logical links between polygynandrous mating behavior, increased pathogen exposure, and greater immunoactivity, these relationships have seldom been examined in nonhuman vertebrates. To explore interactions among these variables in a different lineage of mammals, we used RNAseq to study the gene expression profiles of liver tissue-a highly immunoactive organ-from sympatric populations of the monogamous California mouse (Peromyscus californicus) and two polygynandrous congeners (P. maniculatus and P. boylii). Differential expression and co-expression analyses revealed distinct patterns of gene activity among species, with much of this variation associated with differences in mating system. This tendency was particularly pronounced for MHC genes, with multiple MHC Class I genes being upregulated in the two polygynandrous species, as expected if exposure to sexually transmitted pathogens varies with mating system. Our results underscore the role of mating behavior in influencing patterns of gene expression and highlight the use of emerging transcriptomic tools in behavioral studies of free-living animals.
Collapse
Affiliation(s)
- Jesyka Meléndez‐Rosa
- Department of Integrative BiologyUniversity of CaliforniaBerkeleyCalifornia
- Museum of Vertebrate ZoologyUniversity of CaliforniaBerkeleyCalifornia
| | - Ke Bi
- Museum of Vertebrate ZoologyUniversity of CaliforniaBerkeleyCalifornia
- Computational Genomics Resource LaboratoryUniversity of CaliforniaBerkeleyCalifornia
| | - Eileen A. Lacey
- Department of Integrative BiologyUniversity of CaliforniaBerkeleyCalifornia
- Museum of Vertebrate ZoologyUniversity of CaliforniaBerkeleyCalifornia
| |
Collapse
|
21
|
Bi K, Linderoth T, Singhal S, Vanderpool D, Patton JL, Nielsen R, Moritz C, Good JM. Temporal genomic contrasts reveal rapid evolutionary responses in an alpine mammal during recent climate change. PLoS Genet 2019; 15:e1008119. [PMID: 31050681 PMCID: PMC6519841 DOI: 10.1371/journal.pgen.1008119] [Citation(s) in RCA: 50] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2018] [Revised: 05/15/2019] [Accepted: 04/01/2019] [Indexed: 12/12/2022] Open
Abstract
Many species have experienced dramatic changes in their abundance and distribution during recent climate change, but it is often unclear whether such ecological responses are accompanied by evolutionary change. We used targeted exon sequencing of 294 museum specimens (160 historic, 134 modern) to generate independent temporal genomic contrasts spanning a century of climate change (1911-2012) for two co-distributed chipmunk species: an endemic alpine specialist (Tamias alpinus) undergoing severe range contraction and a stable mid-elevation species (T. speciosus). Using a novel analytical approach, we reconstructed the demographic histories of these populations and tested for evidence of recent positive directional selection. Only the retracting species showed substantial population genetic fragmentation through time and this was coupled with positive selection and substantial shifts in allele frequencies at a gene, Alox15, involved in regulation of inflammation and response to hypoxia. However, these rapid population and gene-level responses were not detected in an analogous temporal contrast from another area where T. alpinus has also undergone severe range contraction. Collectively, these results highlight that evolutionary responses may be variable and context dependent across populations, even when they show seemingly synchronous ecological shifts. Our results demonstrate that temporal genomic contrasts can be used to detect very recent evolutionary responses within and among contemporary populations, even in the face of complex demographic changes. Given the wealth of specimens archived in natural history museums, comparative analyses of temporal population genomic data have the potential to improve our understanding of recent and ongoing evolutionary responses to rapidly changing environments.
Collapse
Affiliation(s)
- Ke Bi
- Museum of Vertebrate Zoology, University of California, Berkeley, California, United States of America
- Computational Genomics Resource Laboratory (CGRL), California Institute for Quantitative Biosciences (QB3), University of California, Berkeley, California, United States of America
| | - Tyler Linderoth
- Museum of Vertebrate Zoology, University of California, Berkeley, California, United States of America
- Department of Integrative Biology, University of California, Berkeley, California, United States of America
| | - Sonal Singhal
- Museum of Vertebrate Zoology, University of California, Berkeley, California, United States of America
- Department of Integrative Biology, University of California, Berkeley, California, United States of America
| | - Dan Vanderpool
- Division of Biological Sciences, University of Montana, Missoula, Montana, United States of America
| | - James L. Patton
- Museum of Vertebrate Zoology, University of California, Berkeley, California, United States of America
- Department of Integrative Biology, University of California, Berkeley, California, United States of America
| | - Rasmus Nielsen
- Department of Integrative Biology, University of California, Berkeley, California, United States of America
| | - Craig Moritz
- Museum of Vertebrate Zoology, University of California, Berkeley, California, United States of America
- Department of Integrative Biology, University of California, Berkeley, California, United States of America
- Research School of Biology and Centre for Biodiversity Analysis, Australian National University, Canberra, ACT, Australia
| | - Jeffrey M. Good
- Division of Biological Sciences, University of Montana, Missoula, Montana, United States of America
- Wildlife Biology Program, University of Montana, Missoula, MT, United States of America
| |
Collapse
|
22
|
Wei C, Li M, Qin J, Xu Y, Zhang Y, Wang H. Transcriptome analysis reveals the effects of grafting on sweetpotato scions during the full blooming stages. Genes Genomics 2019; 41:895-907. [PMID: 31030407 DOI: 10.1007/s13258-019-00823-9] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2018] [Accepted: 04/20/2019] [Indexed: 01/21/2023]
Abstract
BACKGROUND Sweetpotato (Ipomoea batatas) is a hexaploid plant and generally most genotypes do not flower at all in sub-tropics. Heterografting was carried out between sweetpotato cultivar 'Xushu 18' and Japanese morning glory (Ipomoea nil). With sweetpotato as 'scion' and I. nil as 'rootstock', sweetpotato was induced flowering in the autumn. However, little is known about the molecular mechanisms underlying sweetpotato responses to grafting, especially during the full blooming stages. OBJECTIVES To investigate the poorly understood molecular responses underlying the grafting-induced phenotypic processes in sweetpotato at full anthesis. METHODS In this study, to explore the transcriptome diversity and complexity of sweetpotato, PacBio Iso-Seq and Illumina RNA-seq analysis were combined to obtain full-length transcripts and to profile the changes in gene expression of five tissues: scion flowers (SF), scion leaves (SL), scion stems (SS), own-rooted leaves (OL) and own-rooted stems (OS). RESULTS A total of 138,151 transcripts were generated with an average length of 2255 bp, and more than 72% (100,396) of the transcripts were full-length. During full blooming, to examine the difference in gene expression of sweetpotato under grafting and natural growth conditions, 7905, 7795 and 15,707 differentially expressed genes were detected in pairwise comparisons of OS versus SS, OL versus SL and SL versus SF, respectively. Moreover, differential transcription of genes associated with anthocyanin biosynthesis, light pathway and photosynthesis, ethylene signal transduction pathway was observed in scion responses to grafting. CONCLUSION Our study is useful in understanding the molecular basis of grafting-induced flowering in grafted sweetpotatoes, and will lay a foundation for further research on sweetpotato breeding in the future.
Collapse
Affiliation(s)
- Changhe Wei
- Key Laboratory of Bio-resources and Eco-environment, Ministry of Education, Sichuan Key Laboratory of Molecular Biology and Biotechnology, College of Life Sciences, Sichuan University, Chengdu, 610064, China
| | - Ming Li
- Key Laboratory of Bio-resources and Eco-environment, Ministry of Education, Sichuan Key Laboratory of Molecular Biology and Biotechnology, College of Life Sciences, Sichuan University, Chengdu, 610064, China.,Institute of Biotechnology and Nuclear Technology, Sichuan Academy of Agricultural Sciences, Chengdu, 610061, China
| | - Jia Qin
- Key Laboratory of Bio-resources and Eco-environment, Ministry of Education, Sichuan Key Laboratory of Molecular Biology and Biotechnology, College of Life Sciences, Sichuan University, Chengdu, 610064, China
| | - Yunfan Xu
- Key Laboratory of Bio-resources and Eco-environment, Ministry of Education, Sichuan Key Laboratory of Molecular Biology and Biotechnology, College of Life Sciences, Sichuan University, Chengdu, 610064, China
| | - Yizheng Zhang
- Key Laboratory of Bio-resources and Eco-environment, Ministry of Education, Sichuan Key Laboratory of Molecular Biology and Biotechnology, College of Life Sciences, Sichuan University, Chengdu, 610064, China
| | - Haiyan Wang
- Key Laboratory of Bio-resources and Eco-environment, Ministry of Education, Sichuan Key Laboratory of Molecular Biology and Biotechnology, College of Life Sciences, Sichuan University, Chengdu, 610064, China.
| |
Collapse
|
23
|
Frias-Soler RC, Villarín Pildaín L, Hotz-Wagenblatt A, Kolibius J, Bairlein F, Wink M. De novo annotation of the transcriptome of the Northern Wheatear ( Oenanthe oenanthe). PeerJ 2018; 6:e5860. [PMID: 30498627 PMCID: PMC6251345 DOI: 10.7717/peerj.5860] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2017] [Accepted: 10/02/2018] [Indexed: 11/20/2022] Open
Abstract
We have sequenced a partial transcriptome of the Northern Wheatear (Oenanthe oenanthe), a species with one of the longest migrations on Earth. The transcriptome was constructed de novo using RNA-Seq sequence data from the pooled mRNA of six different tissues: brain, muscle, intestine, liver, adipose tissue and skin. The samples came from nine captive-bred wheatears collected at three different stages of the endogenous autumn migratory period: (1) lean birds prior the onset of migration, (2) during the fattening stage and (3) individuals at their migratory body mass plateau, when they have almost doubled their lean body mass. The sample structure used to build up the transcriptome of the Northern Wheatears concerning tissue composition and time guarantees the future survey of the regulatory genes involved in the development of the migratory phenotype. Through the pre-migratory period, birds accomplish outstanding physical and behavioural changes that involve all organ systems. Nevertheless, the molecular mechanisms through which birds synchronize and control hyperphagia, fattening, restlessness increase, immunity boosting and tuning the muscles for such endurance flight are still largely unknown. The use of RNA-Seq has emerged as a powerful tool to analyse complex traits on a broad scale, and we believe it can help to characterize the migratory phenotype of wheatears at an unprecedented level. The primary challenge to conduct quantitative transcriptomic studies in non-model species is the availability of a reference transcriptome, which we have constructed and described in this paper. The cDNA was sequenced by pyrosequencing using the Genome Sequencer Roche GS FLX System; with single paired-end reads of about 400 bp. We estimate the total number of genes at 15,640, of which 67% could be annotated using Turkey and Zebra Finch genomes, or protein sequence information from SwissProt and NCBI databases. With our study, we have made a first step towards understanding the migratory phenotype regarding gene expression of a species that has become a model to study birds long-distance migrations.
Collapse
Affiliation(s)
- Roberto Carlos Frias-Soler
- Institute of Pharmacy and Molecular Biotechnology, Heidelberg University, Heidelberg, Baden Württemberg, Germany.,Institute of Avian Research, Wilhelmshaven, Germany
| | - Lilian Villarín Pildaín
- Institute of Pharmacy and Molecular Biotechnology, Heidelberg University, Heidelberg, Baden Württemberg, Germany
| | - Agnes Hotz-Wagenblatt
- Bioinformatics Group, Core Facility Genomics and Proteomics, German Cancer Research Center, Heidelberg University, Heidelberg, Baden Württemberg, Germany
| | - Jonas Kolibius
- Institute of Pharmacy and Molecular Biotechnology, Heidelberg University, Heidelberg, Baden Württemberg, Germany
| | | | - Michael Wink
- Institute of Pharmacy and Molecular Biotechnology, Heidelberg University, Heidelberg, Baden Württemberg, Germany
| |
Collapse
|
24
|
Corl A, Bi K, Luke C, Challa AS, Stern AJ, Sinervo B, Nielsen R. The Genetic Basis of Adaptation following Plastic Changes in Coloration in a Novel Environment. Curr Biol 2018; 28:2970-2977.e7. [PMID: 30197088 DOI: 10.1016/j.cub.2018.06.075] [Citation(s) in RCA: 62] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2018] [Revised: 06/20/2018] [Accepted: 06/29/2018] [Indexed: 02/07/2023]
Abstract
Phenotypic plasticity has been hypothesized to precede and facilitate adaptation to novel environments [1-8], but examples of plasticity preceding adaptation in wild populations are rare (but see [9, 10]). We studied a population of side-blotched lizards, Uta stansburiana, living on a lava flow that formed 22,500 years ago [11] to understand the origin of their novel melanic phenotype that makes them cryptic on the black lava. We found that lizards living on and off of the lava flow exhibited phenotypic plasticity in coloration but also appeared to have heritable differences in pigmentation. We sequenced the exomes of 104 individuals and identified two known regulators of melanin production, PREP and PRKAR1A, which had markedly increased levels of divergence between lizards living on and off the lava flow. The derived variants in PREP and PRKAR1A were only found in the lava population and were associated with increased pigmentation levels in an experimental cohort of hatchling lizards. Simulations suggest that the derived variants in the PREP and PRKAR1A genes arose recently and were under strong positive selection in the lava population. Overall, our results suggest that ancestral plasticity for coloration facilitated initial survival in the lava environment and was followed by genetic changes that modified the phenotype in the direction of the induced plastic response, possibly through de novo mutations. These observations provide a detailed example supporting the hypothesis that plasticity aids in the initial colonization of a novel habitat, with natural selection subsequently refining the phenotype with genetic adaptations to the new environment. VIDEO ABSTRACT.
Collapse
Affiliation(s)
- Ammon Corl
- Museum of Vertebrate Zoology, University of California, Berkeley, 3101 Valley Life Sciences Building, Berkeley, CA 94720, USA; Department of Ecology and Evolutionary Biology, University of California, Santa Cruz, Santa Cruz, CA 95064, USA.
| | - Ke Bi
- Museum of Vertebrate Zoology, University of California, Berkeley, 3101 Valley Life Sciences Building, Berkeley, CA 94720, USA; Computational Genomics Resource Laboratory (CGRL), California Institute for Quantitative Biosciences (QB3), University of California, Berkeley, Berkeley, CA 94720, USA
| | - Claudia Luke
- Center for Environmental Inquiry, Sonoma State University, 1801 East Cotati Ave, Rohnert Park, CA 94928, USA
| | - Akshara Sree Challa
- Museum of Vertebrate Zoology, University of California, Berkeley, 3101 Valley Life Sciences Building, Berkeley, CA 94720, USA
| | - Aaron James Stern
- Department of Integrative Biology, University of California, Berkeley, 3040 Valley Life Sciences Building # 3140, Berkeley, CA 94720-3140, USA; Graduate Group in Computational Biology, University of California, Berkeley, Berkeley, CA 94720, USA
| | - Barry Sinervo
- Department of Ecology and Evolutionary Biology, University of California, Santa Cruz, Santa Cruz, CA 95064, USA
| | - Rasmus Nielsen
- Department of Integrative Biology, University of California, Berkeley, 3040 Valley Life Sciences Building # 3140, Berkeley, CA 94720-3140, USA
| |
Collapse
|
25
|
Phifer-Rixey M, Bi K, Ferris KG, Sheehan MJ, Lin D, Mack KL, Keeble SM, Suzuki TA, Good JM, Nachman MW. The genomic basis of environmental adaptation in house mice. PLoS Genet 2018; 14:e1007672. [PMID: 30248095 PMCID: PMC6171964 DOI: 10.1371/journal.pgen.1007672] [Citation(s) in RCA: 41] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2018] [Revised: 10/04/2018] [Accepted: 08/30/2018] [Indexed: 01/18/2023] Open
Abstract
House mice (Mus musculus) arrived in the Americas only recently in association with European colonization (~400-600 generations), but have spread rapidly and show evidence of local adaptation. Here, we take advantage of this genetic model system to investigate the genomic basis of environmental adaptation in house mice. First, we documented clinal patterns of phenotypic variation in 50 wild-caught mice from a latitudinal transect in Eastern North America. Next, we found that progeny of mice from different latitudes, raised in a common laboratory environment, displayed differences in a number of complex traits related to fitness. Consistent with Bergmann's rule, mice from higher latitudes were larger and fatter than mice from lower latitudes. They also built bigger nests and differed in aspects of blood chemistry related to metabolism. Then, combining exomic, genomic, and transcriptomic data, we identified specific candidate genes underlying adaptive variation. In particular, we defined a short list of genes with cis-eQTL that were identified as candidates in exomic and genomic analyses, all of which have known ties to phenotypes that vary among the studied populations. Thus, wild mice and the newly developed strains represent a valuable resource for future study of the links between genetic variation, phenotypic variation, and climate.
Collapse
Affiliation(s)
- Megan Phifer-Rixey
- Department of Biology, Monmouth University, West Long Branch, New Jersey, United States of America
- Department of Integrative Biology and Museum of Vertebrate Zoology, University of California, Berkeley, Berkeley, California, United States of America
| | - Ke Bi
- Department of Integrative Biology and Museum of Vertebrate Zoology, University of California, Berkeley, Berkeley, California, United States of America
- Computational Genomics Resource Laboratory, California Institute for Quantitative Biosciences, University of California, Berkeley, Berkeley, California, United States of America
| | - Kathleen G. Ferris
- Department of Integrative Biology and Museum of Vertebrate Zoology, University of California, Berkeley, Berkeley, California, United States of America
| | - Michael J. Sheehan
- Department of Integrative Biology and Museum of Vertebrate Zoology, University of California, Berkeley, Berkeley, California, United States of America
- Department of Neurobiology and Behavior, Cornell University, Ithaca, New York, United States of America
| | - Dana Lin
- Department of Integrative Biology and Museum of Vertebrate Zoology, University of California, Berkeley, Berkeley, California, United States of America
| | - Katya L. Mack
- Department of Integrative Biology and Museum of Vertebrate Zoology, University of California, Berkeley, Berkeley, California, United States of America
| | - Sara M. Keeble
- Division of Biological Sciences, University of Montana, Missoula, Missoula, Montana, United States of America
- Department of Molecular and Computational Biology, University of Southern California, Los Angeles, Los Angeles, California, United States of America
| | - Taichi A. Suzuki
- Department of Integrative Biology and Museum of Vertebrate Zoology, University of California, Berkeley, Berkeley, California, United States of America
| | - Jeffrey M. Good
- Division of Biological Sciences, University of Montana, Missoula, Missoula, Montana, United States of America
| | - Michael W. Nachman
- Department of Integrative Biology and Museum of Vertebrate Zoology, University of California, Berkeley, Berkeley, California, United States of America
| |
Collapse
|
26
|
Almeida AMR, Piñeyro-Nelson A, Yockteng RB, Specht CD. Comparative analysis of whole flower transcriptomes in the Zingiberales. PeerJ 2018; 6:e5490. [PMID: 30155368 PMCID: PMC6110254 DOI: 10.7717/peerj.5490] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2017] [Accepted: 07/30/2018] [Indexed: 01/14/2023] Open
Abstract
The advancement of next generation sequencing technologies (NGS) has revolutionized our ability to generate large quantities of data at a genomic scale. Despite great challenges, these new sequencing technologies have empowered scientists to explore various relevant biological questions on non-model organisms, even in the absence of a complete sequenced reference genome. Here, we analyzed whole flower transcriptome libraries from exemplar species across the monocot order Zingiberales, using a comparative approach in order to gain insight into the evolution of the molecular mechanisms underlying flower development in the group. We identified 4,153 coding genes shared by all floral transcriptomes analyzed, and 1,748 genes that are only retrieved in the Zingiberales. We also identified 666 genes that are unique to the ginger lineage, and 2,001 that are only found in the banana group, while in the outgroup species Dichorisandra thyrsiflora J.C. Mikan (Commelinaceae) we retrieved 2,686 unique genes. It is possible that some of these genes underlie lineage-specific molecular mechanisms of floral diversification. We further discuss the nature of these lineage-specific datasets, emphasizing conserved and unique molecular processes with special emphasis in the Zingiberales. We also briefly discuss the strengths and shortcomings of de novo assembly for the study of developmental processes across divergent taxa from a particular order. Although this comparison is based exclusively on coding genes, with particular emphasis in transcription factors, we believe that the careful study of other regulatory mechanisms, such as non-coding RNAs, might reveal new levels of complexity, which were not explored in this work.
Collapse
Affiliation(s)
- Ana Maria R Almeida
- Department of Biological Sciences, California State University, Hayward, Hayward, CA, United States of America
| | - Alma Piñeyro-Nelson
- Department of Food and Animal Production, Autonomous Metropolitan University, Xochimilco, Mexico City, DF, Mexico
| | - Roxana B Yockteng
- Centro de Investigaciones Tibaitatá, Corporación Colombiana de Investigación Agropecuaria (AGROSAVIA), Tibaitatá, Colombia.,Institut de Systématique, Evolution, Biodiversité-UMR-CNRS, National Museum of Natural History, Paris, France
| | - Chelsea D Specht
- School of Integrative Plant Sciences, Section of Plant Biology and the L.H. Bailey Hortorium, Cornell University, Ithaca, NY, United States of America
| |
Collapse
|
27
|
MacManes MD. The Oyster River Protocol: a multi-assembler and kmer approach for de novo transcriptome assembly. PeerJ 2018; 6:e5428. [PMID: 30083482 PMCID: PMC6078068 DOI: 10.7717/peerj.5428] [Citation(s) in RCA: 60] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2017] [Accepted: 07/21/2018] [Indexed: 11/24/2022] Open
Abstract
Characterizing transcriptomes in non-model organisms has resulted in a massive increase in our understanding of biological phenomena. This boon, largely made possible via high-throughput sequencing, means that studies of functional, evolutionary, and population genomics are now being done by hundreds or even thousands of labs around the world. For many, these studies begin with a de novo transcriptome assembly, which is a technically complicated process involving several discrete steps. The Oyster River Protocol (ORP), described here, implements a standardized and benchmarked set of bioinformatic processes, resulting in an assembly with enhanced qualities over other standard assembly methods. Specifically, ORP produced assemblies have higher Detonate and TransRate scores and mapping rates, which is largely a product of the fact that it leverages a multi-assembler and kmer assembly process, thereby bypassing the shortcomings of any one approach. These improvements are important, as previously unassembled transcripts are included in ORP assemblies, resulting in a significant enhancement of the power of downstream analysis. Further, as part of this study, I show that assembly quality is unrelated with the number of reads generated, above 30 million reads. Code Availability: The version controlled open-source code is available at https://github.com/macmanes-lab/Oyster_River_Protocol. Instructions for software installation and use, and other details are available at http://oyster-river-protocol.rtfd.org/.
Collapse
Affiliation(s)
- Matthew D MacManes
- Department of Molecular, Cellular and Biomedical Sciences, University of New Hampshire, Durham, NH, USA
| |
Collapse
|
28
|
Ashman LG, Bragg JG, Doughty P, Hutchinson MN, Bank S, Matzke NJ, Oliver P, Moritz C. Diversification across biomes in a continental lizard radiation. Evolution 2018; 72:1553-1569. [PMID: 29972238 DOI: 10.1111/evo.13541] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2017] [Revised: 06/15/2018] [Accepted: 06/15/2018] [Indexed: 12/23/2022]
Abstract
Ecological opportunity is a powerful driver of evolutionary diversification, and predicts rapid lineage and phenotypic diversification following colonization of competitor-free habitats. Alternatively, topographic or environmental heterogeneity could be key to generating and sustaining diversity. We explore these hypotheses in a widespread lineage of Australian lizards: the Gehyra variegata group. This clade occurs across two biomes: the Australian monsoonal tropics (AMT), where it overlaps a separate, larger bodied clade of Gehyra and is largely restricted to rocks; and in the larger Australian arid zone (AAZ) where it has no congeners and occupies trees and rocks. New phylogenomic data and coalescent analyses of AAZ taxa resolve lineages and their relationships and reveal high diversity in the western AAZ (Pilbara region). The AMT and AAZ radiations represent separate radiations with no difference in speciation rates. Most taxa occur on rocks, with small geographic ranges relative to widespread generalist taxa across the vast central AAZ. Rock-dwelling and generalist taxa differ morphologically, but only the lineage-poor central AAZ taxa have accelerated evolution. This accords with increasing evidence that lineage and morphological diversity are poorly correlated, and suggests environmental heterogeneity and refugial dynamics have been more important than ecological release in elevating lineage diversity.
Collapse
Affiliation(s)
- L G Ashman
- Research School of Biology, Australian National University, Canberra, ACT 0200, Australia
| | - J G Bragg
- Research School of Biology, Australian National University, Canberra, ACT 0200, Australia
- Royal Botanic Garden, Sydney, NSW 2000, Australia
| | - P Doughty
- Department of Terrestrial Zoology, Western Australian Museum, Perth, WA 6016, Australia
| | - M N Hutchinson
- South Australian Museum, Adelaide, SA 5000, Australia
- School of Biological Sciences, University of Adelaide, Adelaide, SA 5005, Australia
- School of Biological Sciences, Flinders University, Adelaide, SA 5042, Australia
| | - S Bank
- Research School of Biology, Australian National University, Canberra, ACT 0200, Australia
- Johann-Friedrich-Blumenbach Institute for Zoology and Anthropology, University of Göttingen, Göttingen 37073, Germany
| | - N J Matzke
- Research School of Biology, Australian National University, Canberra, ACT 0200, Australia
- School of Biological Sciences, University of Auckland, Auckland 1142, New Zealand
| | - P Oliver
- Research School of Biology, Australian National University, Canberra, ACT 0200, Australia
- Environmental Futures Research Institute, Griffith University, Nathan, QLD 4111, Australia
- Biodiversity and Geosciences Program, Queensland Museum, Brisbane, QLD 4101, Australia
| | - C Moritz
- Research School of Biology, Australian National University, Canberra, ACT 0200, Australia
| |
Collapse
|
29
|
Evaluating the Performance of De Novo Assembly Methods for Venom-Gland Transcriptomics. Toxins (Basel) 2018; 10:toxins10060249. [PMID: 29921759 PMCID: PMC6024825 DOI: 10.3390/toxins10060249] [Citation(s) in RCA: 47] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2018] [Revised: 06/14/2018] [Accepted: 06/15/2018] [Indexed: 11/17/2022] Open
Abstract
Venom-gland transcriptomics is a key tool in the study of the evolution, ecology, function, and pharmacology of animal venoms. In particular, gene-expression variation and coding sequences gained through transcriptomics provide key information for explaining functional venom variation over both ecological and evolutionary timescales. The accuracy and usefulness of inferences made through transcriptomics, however, is limited by the accuracy of the transcriptome assembly, which is a bioinformatic problem with several possible solutions. Several methods have been employed to assemble venom-gland transcriptomes, with the Trinity assembler being the most commonly applied among them. Although previous evidence of variation in performance among assembly software exists, particularly regarding recovery of difficult-to-assemble multigene families such as snake venom metalloproteinases, much work to date still employs a single assembly method. We evaluated the performance of several commonly used de novo assembly methods for the recovery of both nontoxin transcripts and complete, high-quality venom-gene transcripts across eleven snake and four scorpion transcriptomes. We varied k-mer sizes used by some assemblers to evaluate the impact of k-mer length on transcript recovery. We showed that the recovery of nontoxin transcripts and toxin transcripts is best accomplished through different assembly software, with SDT at smaller k-mer lengths and Trinity being best for nontoxin recovery and a combination of SeqMan NGen and a seed-and-extend approach implemented in Extender as the best means of recovering a complete set of toxin transcripts. In particular, Extender was the only means tested capable of assembling multiple isoforms of the diverse snake venom metalloproteinase family, while traditional approaches such as Trinity recovered at most one metalloproteinase transcript. Our work demonstrated that traditional metrics of assembly performance are not predictive of performance in the recovery of complete and high quality toxin genes. Instead, effective venom-gland transcriptomic studies should combine and quality-filter the results of several assemblers with varying algorithmic strategies.
Collapse
|
30
|
Cotoras DD, Bi K, Brewer MS, Lindberg DR, Prost S, Gillespie RG. Co-occurrence of ecologically similar species of Hawaiian spiders reveals critical early phase of adaptive radiation. BMC Evol Biol 2018; 18:100. [PMID: 29921226 PMCID: PMC6009049 DOI: 10.1186/s12862-018-1209-y] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2017] [Accepted: 06/05/2018] [Indexed: 01/15/2023] Open
Abstract
Background The processes through which populations originate and diversify ecologically in the initial stages of adaptive radiation are little understood because we lack information on critical steps of early divergence. A key question is, at what point do closely related species interact, setting the stage for competition and ecological specialization? The Hawaiian Islands provide an ideal system to explore the early stages of adaptive radiation because the islands span ages from 0.5–5 Mya. Hawaiian spiders in the genus Tetragnatha have undergone adaptive radiation, with one lineage (“spiny legs”) showing four different ecomorphs (green, maroon, large brown, small brown); one representative of each ecomorph is generally found at any site on the older islands. Given that the early stages of adaptive radiation are characterized by allopatric divergence between populations of the same ecomorph, the question is, what are the steps towards subsequent co-occurrence of different ecomorphs? Using a transcriptome-based exon capture approach, we focus on early divergence among close relatives of the green ecomorph to understand processes associated with co-occurrence within the same ecomorph at the early stages of adaptive radiation. Results The major outcomes from the current study are first that closely related species within the same green ecomorph of spiny leg Tetragnatha co-occur on the same single volcano on East Maui, and second that there is no evidence of genetic admixture between these ecologically equivalent species. Further, that multiple genetic lineages exist on a single volcano on Maui suggests that there are no inherent dispersal barriers and that the observed limited distribution of taxa reflects competitive exclusion. Conclusions The observation of co-occurrence of ecologically equivalent species on the young volcano of Maui provides a missing link in the process of adaptive radiation between the point when recently divergent species of the same ecomorph occur in allopatry, to the point where different ecomorphs co-occur at a site, as found throughout the older islands. More importantly, the ability of close relatives of the same ecomorph to interact, without admixture, may provide the conditions necessary for ecological divergence and independent evolution of ecomorphs associated with adaptive radiation. Electronic supplementary material The online version of this article (10.1186/s12862-018-1209-y) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Darko D Cotoras
- Department of Integrative Biology, University of California, 3060 Valley Life Sciences Building, Berkeley, CA, 94720-3140, USA. .,Department of Ecology & Evolutionary Biology, University of California Santa Cruz, Santa Cruz, CA, 95064, USA. .,Department of Entomology / Center for Comparative Genomics, California Academy of Sciences, San Francisco, CA, 94118, USA.
| | - Ke Bi
- Museum of Vertebrate Zoology, University of California, 3101 Valley Life Sciences Building, Berkeley, CA, 94720-3160, USA.,Computational Genomics Resource Laboratory (CGRL), California Institute for Quantitative Biosciences (QB3), University of California, Berkeley, CA, 94720-3102, USA
| | - Michael S Brewer
- Department of Biology, East Carolina University, 1000 E 5th St, Greenville, NC, 27858-4353, USA
| | - David R Lindberg
- Department of Integrative Biology, University of California, 3060 Valley Life Sciences Building, Berkeley, CA, 94720-3140, USA.,Museum of Paleontology, University of California, 1101 Valley Life Sciences Building, Berkeley, CA, 94720, USA
| | - Stefan Prost
- Department of Integrative Biology, University of California, 3060 Valley Life Sciences Building, Berkeley, CA, 94720-3140, USA.,Department of Biology, Stanford University, Stanford, CA, 94305-5020, USA
| | - Rosemary G Gillespie
- Department of Environmental Science, University of California, 130 Mulford Hall, Berkeley, CA, 94720-3114, USA
| |
Collapse
|
31
|
Lin D, Bi K, Conroy CJ, Lacey EA, Schraiber JG, Bowie RCK. Mito-nuclear discordance across a recent contact zone for California voles. Ecol Evol 2018; 8:6226-6241. [PMID: 29988439 PMCID: PMC6024151 DOI: 10.1002/ece3.4129] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2017] [Revised: 03/16/2018] [Accepted: 03/29/2018] [Indexed: 11/17/2022] Open
Abstract
To examine the processes that maintain genetic diversity among closely related taxa, we investigated the dynamics of introgression across a contact zone between two lineages of California voles (Microtus californicus). We tested the prediction that introgression of nuclear loci would be greater than that for mitochondrial loci, assuming ongoing gene flow across the contact zone. We also predicted that genomic markers would show a mosaic pattern of differentiation across this zone, consistent with genomes that are semi-permeable. Using mitochondrial cytochrome b sequences and genome-wide loci developed via ddRAD-seq, we analyzed genetic variation for 10 vole populations distributed along the central California coast; this transect included populations from within the distributions of both parental lineages as well as the putative contact zone. Our analyses revealed that (1) the two lineages examined are relatively young, having diverged ca. 8.5-54 kya, (2) voles from the contact zone in Santa Barbara County did not include F1 or early generation backcrossed individuals, and (3) there appeared to be little to no recurrent gene flow across the contact zone. Introgression patterns for mitochondrial and nuclear markers were not concordant; only mitochondrial markers revealed evidence of introgression, putatively due to historical hybridization. These differences in genetic signatures are intriguing given that the contact zone occurs in a region of continuous vole habitat, with no evidence of past or present physical barriers. Future studies that examine specific isolating mechanisms, such as microhabitat use and mate choice, will facilitate our understanding of how genetic boundaries are maintained in this system.
Collapse
Affiliation(s)
- Dana Lin
- Museum of Vertebrate ZoologyUniversity of California, BerkeleyBerkeleyCalifornia
- Department of Integrative BiologyUniversity of California, BerkeleyBerkeleyCalifornia
| | - Ke Bi
- Computational Genomics Resource LaboratoryCalifornia Institute for Quantitative BiosciencesUniversity of California, BerkeleyBerkeleyCalifornia
| | - Christopher J. Conroy
- Department of Integrative BiologyUniversity of California, BerkeleyBerkeleyCalifornia
| | - Eileen A. Lacey
- Museum of Vertebrate ZoologyUniversity of California, BerkeleyBerkeleyCalifornia
- Department of Integrative BiologyUniversity of California, BerkeleyBerkeleyCalifornia
| | - Joshua G. Schraiber
- Department of BiologyCenter for Computational Genetics and GenomicsTemple UniversityPhiladelphiaPennsylvania
- Institute for Genomics and Evolutionary MedicineTemple UniversityPhiladelphiaPennsylvania
| | - Rauri C. K. Bowie
- Museum of Vertebrate ZoologyUniversity of California, BerkeleyBerkeleyCalifornia
- Department of Integrative BiologyUniversity of California, BerkeleyBerkeleyCalifornia
| |
Collapse
|
32
|
Singhal S, Hoskin CJ, Couper P, Potter S, Moritz C. A Framework for Resolving Cryptic Species: A Case Study from the Lizards of the Australian Wet Tropics. Syst Biol 2018; 67:1061-1075. [DOI: 10.1093/sysbio/syy026] [Citation(s) in RCA: 51] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2017] [Accepted: 03/27/2018] [Indexed: 12/19/2022] Open
Affiliation(s)
- Sonal Singhal
- Museum of Zoology and Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI 48109, USA
- Department of Biology, California State University—Dominguez Hills, Carson, CA 90747, USA
| | - Conrad J Hoskin
- College of Science & Engineering, James Cook University, Townsville, Queensland 4811, Australia
| | - Patrick Couper
- Biodiversity Program, Queensland Museum, South Brisbane, Queensland 4101, Australia
| | - Sally Potter
- Division of Ecology and Evolution, Research School of Biology and Centre for Biodiversity Analysis, Australian National University, Acton, ACT 2601, Australia
| | - Craig Moritz
- Division of Ecology and Evolution, Research School of Biology and Centre for Biodiversity Analysis, Australian National University, Acton, ACT 2601, Australia
| |
Collapse
|
33
|
Bragg JG, Potter S, Afonso Silva AC, Hoskin CJ, Bai BYH, Moritz C. Phylogenomics of a rapid radiation: the Australian rainbow skinks. BMC Evol Biol 2018; 18:15. [PMID: 29402211 PMCID: PMC5800007 DOI: 10.1186/s12862-018-1130-4] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2017] [Accepted: 01/25/2018] [Indexed: 12/13/2022] Open
Abstract
Background The application of target capture with next-generation sequencing now enables phylogenomic analyses of rapidly radiating clades of species. But such analyses are complicated by extensive incomplete lineage sorting, demanding the use of methods that consider this process explicitly, such as the multispecies coalescent (MSC) model. However, the MSC makes strong assumptions about divergence history and population structure, and when using the full Bayesian implementation, current computational limits mean that relatively few loci and samples can be analysed for even modest sized radiations. We explore these issues through analyses of an extensive (> 1000 loci) dataset for the Australian rainbow skinks. This group consists of 3 genera and 41 described species, which likely diversified rapidly in Australia during the mid-late Miocene to occupy rainforest, woodland, and rocky habitats with corresponding diversity of morphology and breeding colouration. Previous phylogenetic analyses of this group have revealed short inter-nodes and high discordance among loci, limiting the resolution of inferred trees. A further complication is that many species have deep phylogeographic structure – this poses the question of how to sample individuals within species for analyses using the MSC. Results Phylogenies obtained using concatenation and summary coalescent species tree approaches to the full dataset are well resolved with generally consistent topology, including for previously intractable relationships near the base of the clade. As expected, branch lengths at the tips are substantially overestimated using concatenation. Comparisons of different strategies for sampling haplotypes for full Bayesian MSC analyses (for one clade and using smaller sets of loci) revealed, unexpectedly, that combining haplotypes across divergent phylogeographic lineages yielded consistent species trees. Conclusions This study of more than 1000 loci provides a strongly-supported estimate of the phylogeny of the Australian rainbow skinks, which will inform future research on the evolution and taxonomy of this group. Our analyses suggest that species tree estimation with the MSC can be quite robust to violation of the assumption that the individuals representing a taxon are sampled from a panmictic population. Electronic supplementary material The online version of this article (10.1186/s12862-018-1130-4) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Jason G Bragg
- Research School of Biology and Centre for Biodiversity Analysis, Australian National University, Canberra, Australia. .,Herbarium of NSW, Royal Botanic Gardens & Domain Trust, Sydney, Australia.
| | - Sally Potter
- Research School of Biology and Centre for Biodiversity Analysis, Australian National University, Canberra, Australia
| | - Ana C Afonso Silva
- Research School of Biology and Centre for Biodiversity Analysis, Australian National University, Canberra, Australia.,cE3c - Centre for Ecology, Evolution and Environmental Changes, Faculdade de Ciências, Universidade de Lisboa, Campo Grande, 1749-016, Lisbon, Portugal
| | - Conrad J Hoskin
- College of Science & Engineering, James Cook University, Qld, Townsville, 4811, Australia
| | - Benjamin Y H Bai
- Research School of Biology and Centre for Biodiversity Analysis, Australian National University, Canberra, Australia.,Present address: Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridgeshire, CB10 1SA, UK
| | - Craig Moritz
- Research School of Biology and Centre for Biodiversity Analysis, Australian National University, Canberra, Australia
| |
Collapse
|
34
|
Blom MPK, Bragg JG, Potter S, Moritz C. Accounting for Uncertainty in Gene Tree Estimation: Summary-Coalescent Species Tree Inference in a Challenging Radiation of Australian Lizards. Syst Biol 2018; 66:352-366. [PMID: 28039387 DOI: 10.1093/sysbio/syw089] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2015] [Accepted: 09/27/2016] [Indexed: 11/12/2022] Open
Abstract
Accurate gene tree inference is an important aspect of species tree estimation in a summary-coalescent framework. Yet, in empirical studies, inferred gene trees differ in accuracy due to stochastic variation in phylogenetic signal between targeted loci. Empiricists should, therefore, examine the consistency of species tree inference, while accounting for the observed heterogeneity in gene tree resolution of phylogenomic data sets. Here, we assess the impact of gene tree estimation error on summary-coalescent species tree inference by screening ${\sim}2000$ exonic loci based on gene tree resolution prior to phylogenetic inference. We focus on a phylogenetically challenging radiation of Australian lizards (genus Cryptoblepharus, Scincidae) and explore effects on topology and support. We identify a well-supported topology based on all loci and find that a relatively small number of high-resolution gene trees can be sufficient to converge on the same topology. Adding gene trees with decreasing resolution produced a generally consistent topology, and increased confidence for specific bipartitions that were poorly supported when using a small number of informative loci. This corroborates coalescent-based simulation studies that have highlighted the need for a large number of loci to confidently resolve challenging relationships and refutes the notion that low-resolution gene trees introduce phylogenetic noise. Further, our study also highlights the value of quantifying changes in nodal support across locus subsets of increasing size (but decreasing gene tree resolution). Such detailed analyses can reveal anomalous fluctuations in support at some nodes, suggesting the possibility of model violation. By characterizing the heterogeneity in phylogenetic signal among loci, we can account for uncertainty in gene tree estimation and assess its effect on the consistency of the species tree estimate. We suggest that the evaluation of gene tree resolution should be incorporated in the analysis of empirical phylogenomic data sets. This will ultimately increase our confidence in species tree estimation using summary-coalescent methods and enable us to exploit genomic data for phylogenetic inference. [Coalescence; concatenation; Cryptoblepharus; exon capture; gene tree; phylogenomics; species tree.].
Collapse
Affiliation(s)
- Mozes P K Blom
- Research School of Biology, Australian National University, Canberra ACT 0200, Australia
| | - Jason G Bragg
- Research School of Biology, Australian National University, Canberra ACT 0200, Australia
| | - Sally Potter
- Research School of Biology, Australian National University, Canberra ACT 0200, Australia
| | - Craig Moritz
- Research School of Biology, Australian National University, Canberra ACT 0200, Australia
| |
Collapse
|
35
|
Potter S, Xue AT, Bragg JG, Rosauer DF, Roycroft EJ, Moritz C. Pleistocene climatic changes drive diversification across a tropical savanna. Mol Ecol 2017; 27:520-532. [DOI: 10.1111/mec.14441] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2017] [Revised: 10/20/2017] [Accepted: 10/31/2017] [Indexed: 01/06/2023]
Affiliation(s)
- Sally Potter
- Research School of Biology The Australian National University Acton ACT Australia
- Centre for Biodiversity Analysis Acton ACT Australia
| | - Alexander T. Xue
- Department of Biology City University of New York New York NY USA
- Department of Genetics Rutgers University Piscataway NJ USA
| | - Jason G. Bragg
- Research School of Biology The Australian National University Acton ACT Australia
- Centre for Biodiversity Analysis Acton ACT Australia
| | - Dan F. Rosauer
- Research School of Biology The Australian National University Acton ACT Australia
- Centre for Biodiversity Analysis Acton ACT Australia
| | - Emily J. Roycroft
- School of Biosciences The University of Melbourne Parkville Vic. Australia
- Sciences Department Museums Victoria Melbourne Vic. Australia
| | - Craig Moritz
- Research School of Biology The Australian National University Acton ACT Australia
- Centre for Biodiversity Analysis Acton ACT Australia
| |
Collapse
|
36
|
Blom MPK, Horner P, Moritz C. Convergence across a continent: adaptive diversification in a recent radiation of Australian lizards. Proc Biol Sci 2017; 283:rspb.2016.0181. [PMID: 27306048 DOI: 10.1098/rspb.2016.0181] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2016] [Accepted: 05/19/2016] [Indexed: 12/11/2022] Open
Abstract
Recent radiations are important to evolutionary biologists, because they provide an opportunity to study the mechanisms that link micro- and macroevolution. The role of ecological speciation during adaptive radiation has been intensively studied, but radiations can arise from a diversity of evolutionary processes; in particular, on large continental landmasses where allopatric speciation might frequently precede ecological differentiation. It is therefore important to establish a phylogenetic and ecological framework for recent continental-scale radiations that are species-rich and ecologically diverse. Here, we use a genomic (approx. 1 200 loci, exon capture) approach to fit branch lengths on a summary-coalescent species tree and generate a time-calibrated phylogeny for a recent and ecologically diverse radiation of Australian scincid lizards; the genus Cryptoblepharus We then combine the phylogeny with a comprehensive phenotypic dataset for over 800 individuals across the 26 species, and use comparative methods to test whether habitat specialization can explain current patterns of phenotypic variation in ecologically relevant traits. We find significant differences in morphology between species that occur in distinct environments and convergence in ecomorphology with repeated habitat shifts across the continent. These results suggest that isolated analogous habitats have provided parallel ecological opportunity and have repeatedly promoted adaptive diversification. By contrast, speciation processes within the same habitat have resulted in distinct lineages with relatively limited morphological variation. Overall, our study illustrates how alternative diversification processes might have jointly stimulated species proliferation across the continent and generated a remarkably diverse group of Australian lizards.
Collapse
Affiliation(s)
- Mozes P K Blom
- Research School of Biology, The Australian National University, Canberra ACT 0200, Australia
| | - Paul Horner
- Museum and Art Gallery of the Northern Territory, GPO Box 4646, Darwin NT 0801, Australia
| | - Craig Moritz
- Research School of Biology, The Australian National University, Canberra ACT 0200, Australia
| |
Collapse
|
37
|
Chery JG, Sass C, Specht CD. Development of single-copy nuclear intron markers for species-level phylogenetics: Case study with Paullinieae (Sapindaceae). APPLICATIONS IN PLANT SCIENCES 2017; 5:apps1700051. [PMID: 28989824 PMCID: PMC5628028 DOI: 10.3732/apps.1700051] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/11/2017] [Accepted: 08/07/2017] [Indexed: 06/07/2023]
Abstract
PREMISE OF THE STUDY We developed a bioinformatic pipeline that leverages a publicly available genome and published transcriptomes to design primers in conserved coding sequences flanking targeted introns of single-copy nuclear loci. Paullinieae (Sapindaceae) is used to demonstrate the pipeline. METHODS AND RESULTS Transcriptome reads phylogenetically closer to the lineage of interest are aligned to the closest genome. Single-nucleotide polymorphisms are called, generating a "pseudoreference" closer to the lineage of interest. Several filters are applied to meet the criteria of single-copy nuclear loci with introns of a desired size. Primers are designed in conserved coding sequences flanking introns. Using this pipeline, we developed nine single-copy nuclear intron markers for Paullinieae. CONCLUSIONS This pipeline is highly flexible and can be used for any group with available genomic and transcriptomic resources. This pipeline led to the development of nine variable markers for phylogenetic study without generating sequence data de novo.
Collapse
Affiliation(s)
- Joyce G. Chery
- Department of Integrative Biology, University of California, Berkeley, 3040 Valley Life Sciences Building #3140, Berkeley, California 94720 USA
- University and Jepson Herbaria, University of California, Berkeley, 1001 Valley Life Sciences Building #2465, Berkeley, California 94720 USA
| | - Chodon Sass
- University and Jepson Herbaria, University of California, Berkeley, 1001 Valley Life Sciences Building #2465, Berkeley, California 94720 USA
- Department of Plant and Microbial Biology, University of California, Berkeley, 111 Koshland Hall, Berkeley, California 94720 USA
| | - Chelsea D. Specht
- School of Integrative Plant Sciences, Section of Plant Biology, Cornell University, 412 Mann Library Building, Ithaca, New York 14853 USA
| |
Collapse
|
38
|
Afonso Silva AC, Bragg JG, Potter S, Fernandes C, Coelho MM, Moritz C. Tropical specialist vs. climate generalist: Diversification and demographic history of sister species of
Carlia
skinks from northwestern Australia. Mol Ecol 2017; 26:4045-4058. [DOI: 10.1111/mec.14185] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2016] [Revised: 04/29/2016] [Accepted: 05/02/2017] [Indexed: 01/05/2023]
Affiliation(s)
- Ana C. Afonso Silva
- Research School of Biology and Centre for Biodiversity Analysis Australian National University Acton ACT Australia
- cE3c ‐ Centre for Ecology, Evolution and Environmental Changes Departamento de Biologia Animal Faculdade de Ciências Universidade de Lisboa Lisboa Portugal
| | - Jason G. Bragg
- Research School of Biology and Centre for Biodiversity Analysis Australian National University Acton ACT Australia
- Royal Botanic Garden Sydney NSW Australia
| | - Sally Potter
- Research School of Biology and Centre for Biodiversity Analysis Australian National University Acton ACT Australia
| | - Carlos Fernandes
- cE3c ‐ Centre for Ecology, Evolution and Environmental Changes Departamento de Biologia Animal Faculdade de Ciências Universidade de Lisboa Lisboa Portugal
| | - Maria Manuela Coelho
- cE3c ‐ Centre for Ecology, Evolution and Environmental Changes Departamento de Biologia Animal Faculdade de Ciências Universidade de Lisboa Lisboa Portugal
| | - Craig Moritz
- Research School of Biology and Centre for Biodiversity Analysis Australian National University Acton ACT Australia
| |
Collapse
|
39
|
Singhal S, Bi K. History cleans up messes: The impact of time in driving divergence and introgression in a tropical suture zone. Evolution 2017; 71:1888-1899. [DOI: 10.1111/evo.13278] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2017] [Accepted: 05/04/2017] [Indexed: 02/01/2023]
Affiliation(s)
- Sonal Singhal
- Department of Ecology and Evolutionary Biology University of Michigan, 830 North University Ann Arbor Michigan 48109
- Museum of Zoology University of Michigan, 1109 Geddes Avenue Ann Arbor Michigan 48109
| | - Ke Bi
- Museum of Vertebrate Zoology University of California, Berkeley, 3101 Valley Life Sciences Building Berkeley California 94720
- Computational Genomics Resource Laboratory, California Institute for Quantitative Biosciences University of California Berkeley California 94720
| |
Collapse
|
40
|
Pimsler ML, Jackson JM, Lozier JD. Population genomics reveals a candidate gene involved in bumble bee pigmentation. Ecol Evol 2017; 7:3406-3413. [PMID: 28515876 PMCID: PMC5433978 DOI: 10.1002/ece3.2935] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2016] [Revised: 02/27/2017] [Accepted: 03/03/2017] [Indexed: 01/03/2023] Open
Abstract
Variation in bumble bee color patterns is well‐documented within and between species. Identifying the genetic mechanisms underlying such variation may be useful in revealing evolutionary forces shaping rapid phenotypic diversification. The widespread North American species Bombus bifarius exhibits regional variation in abdominal color forms, ranging from red‐banded to black‐banded phenotypes and including geographically and phenotypically intermediate forms. Identifying genomic regions linked to this variation has been complicated by strong, near species level, genome‐wide differentiation between red‐ and black‐banded forms. Here, we instead focus on the closely related black‐banded and intermediate forms that both belong to the subspecies B. bifarius nearcticus. We analyze an RNA sequencing (RNAseq) data set and identify a cluster of single nucleotide polymorphisms (SNPs) within one gene, Xanthine dehydrogenase/oxidase‐like, that exhibit highly unusual differentiation compared to the rest of the sequenced genome. Homologs of this gene contribute to pigmentation in other insects, and results thus represent a strong candidate for investigating the genetic basis of pigment variation in B. bifarius and other bumble bee mimicry complexes.
Collapse
Affiliation(s)
- Meaghan L Pimsler
- Department of Biological Sciences University of Alabama Tuscaloosa AL USA
| | - Jason M Jackson
- Department of Biological Sciences University of Alabama Tuscaloosa AL USA
| | - Jeffrey D Lozier
- Department of Biological Sciences University of Alabama Tuscaloosa AL USA
| |
Collapse
|
41
|
Bragg JG, Potter S, Bi K, Catullo R, Donnellan SC, Eldridge MDB, Joseph L, Keogh JS, Oliver P, Rowe KC, Moritz C. Resources for phylogenomic analyses of Australian terrestrial vertebrates. Mol Ecol Resour 2016; 17:869-876. [DOI: 10.1111/1755-0998.12633] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2016] [Revised: 10/27/2016] [Accepted: 11/02/2016] [Indexed: 12/01/2022]
Affiliation(s)
- Jason G. Bragg
- Research School of Biology and Centre for Biodiversity Analysis; The Australian National University; Canberra ACT 0200 Australia
| | - Sally Potter
- Research School of Biology and Centre for Biodiversity Analysis; The Australian National University; Canberra ACT 0200 Australia
| | - Ke Bi
- Computational Genomics Resource Laboratory (CGRL); California Institute for Quantitative Biosciences (QB3); University of California; Berkeley CA 94720 USA
| | - Renee Catullo
- Biological Sciences; Macquarie University; Sydney NSW 2109 Australia
| | | | - Mark D. B. Eldridge
- Australian Museum Research Institute, Australian Museum; 1 William St Sydney NSW 2010 Australia
| | - Leo Joseph
- The Commonwealth Scientific and Industrial Research Organization; National Research Collections Australia; GPO Box 1700 Canberra ACT 2601 Australia
| | - J. Scott Keogh
- Research School of Biology and Centre for Biodiversity Analysis; The Australian National University; Canberra ACT 0200 Australia
| | - Paul Oliver
- Research School of Biology and Centre for Biodiversity Analysis; The Australian National University; Canberra ACT 0200 Australia
| | | | - Craig Moritz
- Research School of Biology and Centre for Biodiversity Analysis; The Australian National University; Canberra ACT 0200 Australia
| |
Collapse
|
42
|
Portik DM, Smith LL, Bi K. An evaluation of transcriptome‐based exon capture for frog phylogenomics across multiple scales of divergence (Class: Amphibia, Order: Anura). Mol Ecol Resour 2016; 16:1069-83. [DOI: 10.1111/1755-0998.12541] [Citation(s) in RCA: 73] [Impact Index Per Article: 9.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2015] [Revised: 04/28/2016] [Accepted: 05/06/2016] [Indexed: 11/26/2022]
Affiliation(s)
- Daniel M. Portik
- Museum of Vertebrate Zoology University of California Berkeley CA 94720 USA
- Department of Integrative Biology University of California Valley Life Sciences Building Berkeley CA 94720 USA
- Department of Biology The University of Texas at Arlington Arlington TX 76019 USA
| | - Lydia L. Smith
- Museum of Vertebrate Zoology University of California Berkeley CA 94720 USA
- Department of Integrative Biology University of California Valley Life Sciences Building Berkeley CA 94720 USA
| | - Ke Bi
- Museum of Vertebrate Zoology University of California Berkeley CA 94720 USA
- Computational Genomics Resource Laboratory (CGRL) California Institute for Quantitative Biosciences (QB3) University of California Berkeley CA 94720 USA
| |
Collapse
|
43
|
Gerchen JF, Reichert SJ, Röhr JT, Dieterich C, Kloas W, Stöck M. A Single Transcriptome of a Green Toad (Bufo viridis) Yields Candidate Genes for Sex Determination and -Differentiation and Non-Anonymous Population Genetic Markers. PLoS One 2016; 11:e0156419. [PMID: 27232626 PMCID: PMC4883742 DOI: 10.1371/journal.pone.0156419] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2016] [Accepted: 05/15/2016] [Indexed: 12/13/2022] Open
Abstract
Large genome size, including immense repetitive and non-coding fractions, still present challenges for capacity, bioinformatics and thus affordability of whole genome sequencing in most amphibians. Here, we test the performance of a single transcriptome to understand whether it can provide a cost-efficient resource for species with large unknown genomes. Using RNA from six different tissues from a single Palearctic green toad (Bufo viridis) specimen and Hiseq2000, we obtained 22,5 Mio reads and publish >100,000 unigene sequences. To evaluate efficacy and quality, we first use this data to identify green toad specific candidate genes, known from other vertebrates for their role in sex determination and differentiation. Of a list of 37 genes, the transcriptome yielded 32 (87%), many of which providing the first such data for this non-model anuran species. However, for many of these genes, only fragments could be retrieved. In order to allow also applications to population genetics, we further used the transcriptome for the targeted development of 21 non-anonymous microsatellites and tested them in genetic families and backcrosses. Eleven markers were specifically developed to be located on the B. viridis sex chromosomes; for eight markers we can indeed demonstrate sex-specific transmission in genetic families. Depending on phylogenetic distance, several markers, which are sex-linked in green toads, show high cross-amplification success across the anuran phylogeny, involving nine systematic anuran families. Our data support the view that single transcriptome sequencing (based on multiple tissues) provides a reliable genomic resource and cost-efficient method for non-model amphibian species with large genome size and, despite limitations, should be considered as long as genome sequencing remains unaffordable for most species.
Collapse
Affiliation(s)
- Jörn F Gerchen
- Department of Ecophysiology and Aquaculture, Leibniz-Institute of Freshwater Ecology and Inland Fisheries (IGB), Berlin, Germany
| | - Samuel J Reichert
- Department of Ecophysiology and Aquaculture, Leibniz-Institute of Freshwater Ecology and Inland Fisheries (IGB), Berlin, Germany
| | - Johannes T Röhr
- Leibniz Institute for Research on Evolution and Biodiversity, Berlin, Germany.,Department of Mathematics and Computer Science, Freie Universität Berlin, Berlin, Germany
| | | | - Werner Kloas
- Department of Ecophysiology and Aquaculture, Leibniz-Institute of Freshwater Ecology and Inland Fisheries (IGB), Berlin, Germany
| | - Matthias Stöck
- Department of Ecophysiology and Aquaculture, Leibniz-Institute of Freshwater Ecology and Inland Fisheries (IGB), Berlin, Germany
| |
Collapse
|
44
|
Potter S, Bragg JG, Peter BM, Bi K, Moritz C. Phylogenomics at the tips: inferring lineages and their demographic history in a tropical lizard, Carlia amax. Mol Ecol 2016; 25:1367-80. [DOI: 10.1111/mec.13546] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2015] [Revised: 12/06/2015] [Accepted: 01/07/2016] [Indexed: 01/04/2023]
Affiliation(s)
- Sally Potter
- Research School of Biology; The Australian National University; Acton Act 2601 Australia
- Centre for Biodiversity Analysis; Acton Act 2601 Australia
| | - Jason G. Bragg
- Research School of Biology; The Australian National University; Acton Act 2601 Australia
- Centre for Biodiversity Analysis; Acton Act 2601 Australia
| | - Benjamin M. Peter
- Department of Integrative Biology; University of California; Berkeley CA 94720-3102 USA
- Department of Human Genetics; University of Chicago; Chicago IL 60637 USA
| | - Ke Bi
- Computational Genomics Resource Laboratory (CGRL); California Institute for Quantitative Biosciences (QB3); University of California; Berkeley CA 94720 USA
| | - Craig Moritz
- Research School of Biology; The Australian National University; Acton Act 2601 Australia
- Centre for Biodiversity Analysis; Acton Act 2601 Australia
| |
Collapse
|
45
|
Lozier JD, Jackson JM, Dillon ME, Strange JP. Population genomics of divergence among extreme and intermediate color forms in a polymorphic insect. Ecol Evol 2016; 6:1075-91. [PMID: 26811748 PMCID: PMC4722823 DOI: 10.1002/ece3.1928] [Citation(s) in RCA: 26] [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/11/2015] [Revised: 12/08/2015] [Accepted: 12/09/2015] [Indexed: 01/01/2023] Open
Abstract
Geographic variation in insect coloration is among the most intriguing examples of rapid phenotypic evolution and provides opportunities to study mechanisms of phenotypic change and diversification in closely related lineages. The bumble bee Bombus bifarius comprises two geographically disparate color groups characterized by red‐banded and black‐banded abdominal pigmentation, but with a range of spatially and phenotypically intermediate populations across western North America. Microsatellite analyses have revealed that B. bifarius in the USA are structured into two major groups concordant with geography and color pattern, but also suggest ongoing gene flow among regional populations. In this study, we better resolve the relationships among major color groups to better understand evolutionary mechanisms promoting and maintaining such polymorphism. We analyze >90,000 and >25,000 single‐nucleotide polymorphisms derived from transcriptome (RNAseq) and double digest restriction site associated DNA sequencing (ddRAD), respectively, in representative samples from spatial and color pattern extremes in B. bifarius as well as phenotypic and geographic intermediates. Both ddRAD and RNAseq data illustrate substantial genome‐wide differentiation of the red‐banded (eastern) color form from both black‐banded (western) and intermediate (central) phenotypes and negligible differentiation among the latter populations, with no obvious admixture among bees from the two major lineages. Results thus indicate much stronger background differentiation among B. bifarius lineages than expected, highlighting potential challenges for revealing loci underlying color polymorphism from population genetic data alone. These findings will have significance for resolving taxonomic confusion in this species and in future efforts to investigate color‐pattern evolution in B. bifarius and other polymorphic bumble bee species.
Collapse
Affiliation(s)
- Jeffrey D Lozier
- Department of Biological Sciences University of Alabama Tuscaloosa Alabama
| | - Jason M Jackson
- Department of Biological Sciences University of Alabama Tuscaloosa Alabama
| | - Michael E Dillon
- Department of Zoology & Physiology and Program in Ecology University of Wyoming Laramie Wyoming
| | - James P Strange
- USDA-ARS Pollinating Insect Research Unit Utah State University Logan Utah
| |
Collapse
|
46
|
Sass C, Iles WJD, Barrett CF, Smith SY, Specht CD. Revisiting the Zingiberales: using multiplexed exon capture to resolve ancient and recent phylogenetic splits in a charismatic plant lineage. PeerJ 2016; 4:e1584. [PMID: 26819846 PMCID: PMC4727956 DOI: 10.7717/peerj.1584] [Citation(s) in RCA: 43] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2015] [Accepted: 12/18/2015] [Indexed: 12/15/2022] Open
Abstract
The Zingiberales are an iconic order of monocotyledonous plants comprising eight families with distinctive and diverse floral morphologies and representing an important ecological element of tropical and subtropical forests. While the eight families are demonstrated to be monophyletic, phylogenetic relationships among these families remain unresolved. Neither combined morphological and molecular studies nor recent attempts to resolve family relationships using sequence data from whole plastomes has resulted in a well-supported, family-level phylogenetic hypothesis of relationships. Here we approach this challenge by leveraging the complete genome of one member of the order, Musa acuminata, together with transcriptome information from each of the other seven families to design a set of nuclear loci that can be enriched from highly divergent taxa with a single array-based capture of indexed genomic DNA. A total of 494 exons from 418 nuclear genes were captured for 53 ingroup taxa. The entire plastid genome was also captured for the same 53 taxa. Of the total genes captured, 308 nuclear and 68 plastid genes were used for phylogenetic estimation. The concatenated plastid and nuclear dataset supports the position of Musaceae as sister to the remaining seven families. Moreover, the combined dataset recovers known intra- and inter-family phylogenetic relationships with generally high bootstrap support. This is a flexible and cost effective method that gives the broader plant biology community a tool for generating phylogenomic scale sequence data in non-model systems at varying evolutionary depths.
Collapse
Affiliation(s)
- Chodon Sass
- Department of Plant and Microbial Biology, Department of Integrative Biology and the University and Jepson Herbaria, University of California, Berkeley , Berkeley, CA , United States
| | - William J D Iles
- Department of Plant and Microbial Biology, Department of Integrative Biology and the University and Jepson Herbaria, University of California, Berkeley , Berkeley, CA , United States
| | - Craig F Barrett
- Division of Plant and Soil Sciences, West Virginia University, Morgantown, WV, United States; Department of Biology, California State University, Los Angeles, Los Angeles, CA, United States
| | - Selena Y Smith
- Department of Earth & Environmental Sciences and the Museum of Paleontology, University of Michigan , Ann Arbor, MI , United States
| | - Chelsea D Specht
- Department of Plant and Microbial Biology, Department of Integrative Biology and the University and Jepson Herbaria, University of California, Berkeley , Berkeley, CA , United States
| |
Collapse
|
47
|
Kono N, Nakamura H, Ito Y, Tomita M, Arakawa K. Evaluation of the impact of RNA preservation methods of spiders for de novo transcriptome assembly. Mol Ecol Resour 2015; 16:662-72. [PMID: 26561354 DOI: 10.1111/1755-0998.12485] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2015] [Revised: 10/26/2015] [Accepted: 10/29/2015] [Indexed: 11/30/2022]
Abstract
With advances in high-throughput sequencing technologies, de novo transcriptome sequencing and assembly has become a cost-effective method to obtain comprehensive genetic information of a species of interest, especially in nonmodel species with large genomes such as spiders. However, high-quality RNA is essential for successful sequencing, and sample preservation conditions require careful consideration for the effective storage of field-collected samples. To this end, we report a streamlined feasibility study of various storage conditions and their effects on de novo transcriptome assembly results. The storage parameters considered include temperatures ranging from room temperature to -80°C; preservatives, including ethanol, RNAlater, TRIzol and RNAlater-ICE; and sample submersion states. As a result, intact RNA was extracted and assembly was successful when samples were preserved at low temperatures regardless of the type of preservative used. The assemblies as well as the gene expression profiles were shown to be robust to RNA degradation, when 30 million 150-bp paired-end reads are obtained. The parameters for sample storage, RNA extraction, library preparation, sequencing and in silico assembly considered in this work provide a guideline for the study of field-collected samples of spiders.
Collapse
Affiliation(s)
- Nobuaki Kono
- Institute for Advanced Biosciences, Keio University, Mizukami 246-2, Kakuganji, Tsuruoka, Yamagata, 997-0052, Japan
| | - Hiroyuki Nakamura
- Spiber Inc., Mizukami 234-1, Kakuganji, Tsuruoka, Yamagata, 997-0052, Japan
| | - Yusuke Ito
- Spiber Inc., Mizukami 234-1, Kakuganji, Tsuruoka, Yamagata, 997-0052, Japan
| | - Masaru Tomita
- Institute for Advanced Biosciences, Keio University, Mizukami 246-2, Kakuganji, Tsuruoka, Yamagata, 997-0052, Japan
| | - Kazuharu Arakawa
- Institute for Advanced Biosciences, Keio University, Mizukami 246-2, Kakuganji, Tsuruoka, Yamagata, 997-0052, Japan
| |
Collapse
|
48
|
Liu H, Lamm MS, Rutherford K, Black MA, Godwin JR, Gemmell NJ. Large-scale transcriptome sequencing reveals novel expression patterns for key sex-related genes in a sex-changing fish. Biol Sex Differ 2015; 6:26. [PMID: 26613014 PMCID: PMC4660848 DOI: 10.1186/s13293-015-0044-8] [Citation(s) in RCA: 67] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 09/11/2015] [Accepted: 11/09/2015] [Indexed: 12/25/2022] Open
Abstract
Background Teleost fishes exhibit remarkably diverse and plastic sexual developmental patterns. One of the most astonishing is the rapid socially controlled female-to-male (protogynous) sex change observed in bluehead wrasses (Thalassoma bifasciatum). Such functional sex change is widespread in marine fishes, including species of commercial importance, yet its underlying molecular basis remains poorly explored. Methods RNA sequencing was performed to characterize the transcriptomic profiles and identify genes exhibiting sex-biased expression in the brain (forebrain and midbrain) and gonads of bluehead wrasses. Functional annotation and enrichment analysis were carried out for the sex-biased genes in the gonad to detect global differences in gene products and genetic pathways between males and females. Results Here we report the first transcriptomic analysis for a protogynous fish. Expression comparison between males and females reveals a large set of genes with sex-biased expression in the gonad, but relatively few such sex-biased genes in the brain. Functional annotation and enrichment analysis suggested that ovaries are mainly enriched for metabolic processes and testes for signal transduction, particularly receptors of neurotransmitters and steroid hormones. When compared to other species, many genes previously implicated in male sex determination and differentiation pathways showed conservation in their gonadal expression patterns in bluehead wrasses. However, some critical female-pathway genes (e.g., rspo1 and wnt4b) exhibited unanticipated expression patterns. In the brain, gene expression patterns suggest that local neurosteroid production and signaling likely contribute to the sex differences observed. Conclusions Expression patterns of key sex-related genes suggest that sex-changing fish predominantly use an evolutionarily conserved genetic toolkit, but that subtle variability in the standard sex-determination regulatory network likely contributes to sexual plasticity in these fish. This study not only provides the first molecular data on a system ideally suited to explore the molecular basis of sexual plasticity and tissue re-engineering, but also sheds some light on the evolution of diverse sex determination and differentiation systems. Electronic supplementary material The online version of this article (doi:10.1186/s13293-015-0044-8) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Hui Liu
- Department of Anatomy, University of Otago, Dunedin, New Zealand
| | - Melissa S Lamm
- Department of Biological Sciences, North Carolina State University, Raleigh, NC USA ; W.M. Keck Center for Behavioral Biology, North Carolina State University, Raleigh, NC USA
| | - Kim Rutherford
- Department of Anatomy, University of Otago, Dunedin, New Zealand
| | - Michael A Black
- Department of Biochemistry, University of Otago, Dunedin, New Zealand
| | - John R Godwin
- Department of Biological Sciences, North Carolina State University, Raleigh, NC USA ; W.M. Keck Center for Behavioral Biology, North Carolina State University, Raleigh, NC USA
| | - Neil J Gemmell
- Department of Anatomy, University of Otago, Dunedin, New Zealand
| |
Collapse
|
49
|
De Novo Assembly of the Pea (Pisum sativum L.) Nodule Transcriptome. Int J Genomics 2015; 2015:695947. [PMID: 26688806 PMCID: PMC4672141 DOI: 10.1155/2015/695947] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2015] [Revised: 09/28/2015] [Accepted: 10/25/2015] [Indexed: 11/17/2022] Open
Abstract
The large size and complexity of the garden pea (Pisum sativum L.) genome hamper its sequencing and the discovery of pea gene resources. Although transcriptome sequencing provides extensive information about expressed genes, some tissue-specific transcripts can only be identified from particular organs under appropriate conditions. In this study, we performed RNA sequencing of polyadenylated transcripts from young pea nodules and root tips on an Illumina GAIIx system, followed by de novo transcriptome assembly using the Trinity program. We obtained more than 58,000 and 37,000 contigs from "Nodules" and "Root Tips" assemblies, respectively. The quality of the assemblies was assessed by comparison with pea expressed sequence tags and transcriptome sequencing project data available from NCBI website. The "Nodules" assembly was compared with the "Root Tips" assembly and with pea transcriptome sequencing data from projects indicating tissue specificity. As a result, approximately 13,000 nodule-specific contigs were found and annotated by alignment to known plant protein-coding sequences and by Gene Ontology searching. Of these, 581 sequences were found to possess full CDSs and could thus be considered as novel nodule-specific transcripts of pea. The information about pea nodule-specific gene sequences can be applied for gene-based markers creation, polymorphism studies, and real-time PCR.
Collapse
|
50
|
Hykin SM, Bi K, McGuire JA. Fixing Formalin: A Method to Recover Genomic-Scale DNA Sequence Data from Formalin-Fixed Museum Specimens Using High-Throughput Sequencing. PLoS One 2015; 10:e0141579. [PMID: 26505622 PMCID: PMC4623518 DOI: 10.1371/journal.pone.0141579] [Citation(s) in RCA: 90] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2015] [Accepted: 10/09/2015] [Indexed: 01/23/2023] Open
Abstract
For 150 years or more, specimens were routinely collected and deposited in natural history collections without preserving fresh tissue samples for genetic analysis. In the case of most herpetological specimens (i.e. amphibians and reptiles), attempts to extract and sequence DNA from formalin-fixed, ethanol-preserved specimens—particularly for use in phylogenetic analyses—has been laborious and largely ineffective due to the highly fragmented nature of the DNA. As a result, tens of thousands of specimens in herpetological collections have not been available for sequence-based phylogenetic studies. Massively parallel High-Throughput Sequencing methods and the associated bioinformatics, however, are particularly suited to recovering meaningful genetic markers from severely degraded/fragmented DNA sequences such as DNA damaged by formalin-fixation. In this study, we compared previously published DNA extraction methods on three tissue types subsampled from formalin-fixed specimens of Anolis carolinensis, followed by sequencing. Sufficient quality DNA was recovered from liver tissue, making this technique minimally destructive to museum specimens. Sequencing was only successful for the more recently collected specimen (collected ~30 ybp). We suspect this could be due either to the conditions of preservation and/or the amount of tissue used for extraction purposes. For the successfully sequenced sample, we found a high rate of base misincorporation. After rigorous trimming, we successfully mapped 27.93% of the cleaned reads to the reference genome, were able to reconstruct the complete mitochondrial genome, and recovered an accurate phylogenetic placement for our specimen. We conclude that the amount of DNA available, which can vary depending on specimen age and preservation conditions, will determine if sequencing will be successful. The technique described here will greatly improve the value of museum collections by making many formalin-fixed specimens available for genetic analysis.
Collapse
Affiliation(s)
- Sarah M. Hykin
- Department of Integrative Biology, 3101 Valley Life Sciences Building, University of California, Berkeley, California, United States of America
- Museum of Vertebrate Zoology, 3101 Valley Life Sciences Building, University of California, Berkeley, California, United States of America
- * E-mail: (SMH); (KB)
| | - Ke Bi
- Computational Genomics Resource Laboratory (CGRL), California Institute for Quantitative Biosciences (QB3), University of California, Berkeley, California, United States of America
- Museum of Vertebrate Zoology, 3101 Valley Life Sciences Building, University of California, Berkeley, California, United States of America
- * E-mail: (SMH); (KB)
| | - Jimmy A. McGuire
- Department of Integrative Biology, 3101 Valley Life Sciences Building, University of California, Berkeley, California, United States of America
- Museum of Vertebrate Zoology, 3101 Valley Life Sciences Building, University of California, Berkeley, California, United States of America
| |
Collapse
|