1
|
Myburgh AM, Barnes A, Henriques R, Daniels SR. Congruent patterns of cryptic cladogenesis revealed using RADseq and Sanger sequencing in a velvet worm species complex (Onychophora: Peripatopsidae: Peripatopsis sedgwicki). Mol Phylogenet Evol 2024; 198:108132. [PMID: 38909874 DOI: 10.1016/j.ympev.2024.108132] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2024] [Revised: 05/24/2024] [Accepted: 06/15/2024] [Indexed: 06/25/2024]
Abstract
In the present study, first generation DNA sequencing (mitochondrial cytochrome c oxidase subunit one, COI) and reduced-representative genomic RADseq data were used to understand the patterns and processes of diversification of the velvet worm, Peripatopsis sedgwicki species complex across its distribution range in South Africa. For the RADseq data, three datasets (two primary and one supplementary) were generated corresponding to 1,259-11,468 SNPs, in order to assess the diversity and phylogeography of the species complex. Tree topologies for the two primary datasets were inferred using maximum likelihood and Bayesian inferences methods. Phylogenetic analyses using the COI datasets retrieved four distinct, well-supported clades within the species complex. Five species delimitation methods applied to the COI data (ASAP, bPTP, bGMYC, STACEY and iBPP) all showed support for the distinction of the Fort Fordyce Nature Reserve specimens. In the main P. sedgwicki species complex, the species delimitation methods revealed a variable number of operational taxonomic units and overestimated the number of putative taxa. Divergence time estimates coupled with the geographic exclusivity of species and phylogeographic results suggest recent cladogenesis during the Plio/Pleistocene. The RADseq data were subjected to a principal components analysis and a discriminant analysis of principal components, under a maximum-likelihood framework. The latter results corroborate the four main clades observed using the COI data, however, applying additional filtering revealed additional diversity. The high overall congruence observed between the RADseq data and COI data suggest that first generation sequence data remain a cheap and effective method for evolutionary studies, although RADseq does provide a far greater resolution of contemporary temporo-spatial patterns.
Collapse
Affiliation(s)
- Angus Macgregor Myburgh
- Department of Botany and Zoology, Private Bag X1, Stellenbosch University, 7602, South Africa
| | - Aaron Barnes
- Department of Botany and Zoology, Private Bag X1, Stellenbosch University, 7602, South Africa
| | - Romina Henriques
- Department of Biochemistry, Genetics and Microbiology, University of Pretoria, South Africa
| | - Savel R Daniels
- Department of Botany and Zoology, Private Bag X1, Stellenbosch University, 7602, South Africa.
| |
Collapse
|
2
|
Peñafiel Loaiza N, Chafe AH, Moraes R M, Oleas NH, Roncal J. Genotyping-by-sequencing informs conservation of Andean palms sources of non-timber forest products. Evol Appl 2024; 17:e13765. [PMID: 39091352 PMCID: PMC11291087 DOI: 10.1111/eva.13765] [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: 04/06/2023] [Revised: 06/30/2024] [Accepted: 07/18/2024] [Indexed: 08/04/2024] Open
Abstract
Conservation and sustainable management of lineages providing non-timber forest products are imperative under the current global biodiversity loss. Most non-timber forest species, however, lack genomic studies that characterize their intraspecific variation and evolutionary history, which inform species' conservation practices. Contrary to many lineages in the Andean biodiversity hotspot that exhibit high diversification, the genus Parajubaea (Arecaceae) has only three species despite the genus' origin 22 million years ago. Two of the three palm species, P. torallyi and P. sunkha, are non-timber forest species endemic to the Andes of Bolivia and are listed as IUCN endangered. The third species, P. cocoides, is a vulnerable species with unknown wild populations. We investigated the evolutionary relationships of Parajubaea species and the genetic diversity and structure of wild Bolivian populations. Sequencing of five low-copy nuclear genes (3753 bp) challenged the hypothesis that P. cocoides is a cultigen that originated from the wild Bolivian species. We further obtained up to 15,134 de novo single-nucleotide polymorphism markers by genotyping-by-sequencing of 194 wild Parajubaea individuals. Our total DNA sequencing effort rejected the taxonomic separation of the two Bolivian species. As expected for narrow endemic species, we observed low genetic diversity, but no inbreeding signal. We found three genetic clusters shaped by geographic distance, which we use to propose three management units. Different percentages of missing genotypic data did not impact the genetic structure of populations. We use the management units to recommend in situ conservation by creating new protected areas, and ex situ conservation through seed collection.
Collapse
Affiliation(s)
- Nicolás Peñafiel Loaiza
- Department of BiologyMemorial University of NewfoundlandSt. John'sNewfoundland and LabradorCanada
- Present address:
Chone y BabahoyoLojaEcuador
| | - Abigail H. Chafe
- Department of BiologyMemorial University of NewfoundlandSt. John'sNewfoundland and LabradorCanada
| | - Mónica Moraes R
- Herbario Nacional de Bolivia, Instituto de EcologíaUniversidad Mayor de San AndrésLa PazBolivia
| | - Nora H. Oleas
- Centro de Investigación de la Biodiversidad y Cambio Climático – BioCamb e Ingeniería en Biodiversidad y Recursos Genéticos, Facultad de Ciencias de Medio AmbienteUniversidad IndoaméricaQuitoEcuador
| | - Julissa Roncal
- Department of BiologyMemorial University of NewfoundlandSt. John'sNewfoundland and LabradorCanada
| |
Collapse
|
3
|
Muharromah AF, Carvajal TM, Regilme MAF, Watanabe K. Fine-scale adaptive divergence and population genetic structure of Aedes aegypti in Metropolitan Manila, Philippines. Parasit Vectors 2024; 17:233. [PMID: 38769579 PMCID: PMC11107013 DOI: 10.1186/s13071-024-06300-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2024] [Accepted: 04/23/2024] [Indexed: 05/22/2024] Open
Abstract
BACKGROUND The adaptive divergence of Aedes aegypti populations to heterogeneous environments can be a driving force behind the recent expansion of their habitat distribution and outbreaks of dengue disease in urbanized areas. In this study, we investigated the population genomics of Ae. aegypti at a regional scale in Metropolitan Manila, Philippines. METHODS We used the Pool-Seq double digestion restriction-site association DNA sequencing (ddRAD-Seq) approach to generate a high number of single nucleotide polymorphisms (SNPs), with the aim to determine local adaptation and compare the population structure with 11 microsatellite markers. A total of 217 Ae. aegypti individuals from seven female and seven male populations collected from Metropolitan Manila were used in the assays. RESULTS We detected 65,473 SNPs across the populations, of which 76 were non-neutral SNPs. Of these non-neutral SNPs, the multivariate regression test associated 50 with eight landscape variables (e.g. open space, forest, etc.) and 29 with five climate variables (e.g. air temperature, humidity, etc.) (P-value range 0.005-0.045) in female and male populations separately. Male and female populations exhibited contrasting spatial divergence, with males exhibiting greater divergence than females, most likely reflecting the different dispersal abilities of male and female mosquitoes. In the comparative analysis of the same Ae. aegypti individuals, the pairwise FST values of 11 microsatellite markers were lower than those of the neutral SNPs, indicating that the neutral SNPs generated via pool ddRAD-Seq were more sensitive in terms of detecting genetic differences between populations at fine-spatial scales. CONCLUSIONS Overall, our study demonstrates the utility of pool ddRAD-Seq for examining genetic differences in Ae. aegypti populations in areas at fine-spatial scales that could inform vector control programs such as Wolbachia-infected mosquito mass-release programs. This in turn would provide information on mosquito population dispersal patterns and the potential barriers to mosquito movement within and around the release area. In addition, the potential of environmental adaptability observed in Ae. aegypti could help population control efforts.
Collapse
Affiliation(s)
- Atikah Fitria Muharromah
- Center for Marine Environmental Studies (CMES), Ehime University, Bunkyo-cho 3, Matsuyama, Ehime, 7908577, Japan
- Graduate School of Science and Engineering, Ehime University, Bunkyo-cho 3, Matsuyama, Ehime, 7908577, Japan
- Department of Tropical Biology, Faculty of Biology, Universitas Gadjah Mada, Yogyakarta, 55281, Indonesia
| | - Thaddeus M Carvajal
- Center for Marine Environmental Studies (CMES), Ehime University, Bunkyo-cho 3, Matsuyama, Ehime, 7908577, Japan
- Biological Control Research Unit, Center for Natural Sciences and Environmental Research, De La Salle University, 2401 Taft Avenue, 1004, Manila, Philippines
| | - Maria Angenica F Regilme
- Center for Marine Environmental Studies (CMES), Ehime University, Bunkyo-cho 3, Matsuyama, Ehime, 7908577, Japan
| | - Kozo Watanabe
- Center for Marine Environmental Studies (CMES), Ehime University, Bunkyo-cho 3, Matsuyama, Ehime, 7908577, Japan.
| |
Collapse
|
4
|
Rainey TA, Tryc EE, Nicholson KE. Comparing skin swabs, buccal swabs, and toe clips for amphibian genetic sampling, a case study with a small anuran ( Acris blanchardi). Biol Methods Protoc 2024; 9:bpae030. [PMID: 38818027 PMCID: PMC11139508 DOI: 10.1093/biomethods/bpae030] [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: 03/28/2024] [Revised: 04/26/2024] [Accepted: 05/14/2024] [Indexed: 06/01/2024] Open
Abstract
Multiple methods for collecting genetic samples from amphibians exist, each with their own implications for study design, animal welfare, and costs. Toe clipping is one common method, but there is ongoing debate regarding its potential detriment. Less invasive methods should be implemented, if efficacious, as amphibians are a particularly vulnerable vertebrate group. Skin and buccal swabbing are less invasive methods for genetic sampling, but the potential for contamination and a lower yield of DNA may exist. To compare these methods, we gathered skin swabs, buccal swabs, and toe clips from the same individuals of a relatively small anuran species, Blanchard's Cricket Frog (Acris blanchardi). We then compared DNA yield, DNA purity, amplification success rate, and genotypic data quality among sample types. We found toe clips and buccal swabs generated similar DNA yield and purity, with skin swabs yielding significantly less DNA of significantly lower purity than the other sample types. Amplification success rate was significantly higher using toe clips compared to the other sample types, though buccal swab samples amplified more readily than skin swabs. Genotypic data from toe clips and buccal swabs did not differ significantly in quality, but skin swab data quality was significantly lowest among sample types. Thus, skin swabbing could produce erroneous data in some situations, but buccal swabbing is likely an effective substitute to toe clipping, even for small species. Our results can help future researchers select which genetic sampling method might best suit their research needs.
Collapse
Affiliation(s)
- Travis A Rainey
- Department of Biology, Central Michigan University, Mount Pleasant, MI 48859, United States
| | - Emily E Tryc
- Department of Biology, Central Michigan University, Mount Pleasant, MI 48859, United States
| | - Kirsten E Nicholson
- Department of Biology, Central Michigan University, Mount Pleasant, MI 48859, United States
| |
Collapse
|
5
|
Koontz AC, Schumacher EK, Spence ES, Hoban SM. Ex situ conservation of two rare oak species using microsatellite and SNP markers. Evol Appl 2024; 17:e13650. [PMID: 38524684 PMCID: PMC10960078 DOI: 10.1111/eva.13650] [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: 10/03/2023] [Revised: 12/27/2023] [Accepted: 01/14/2024] [Indexed: 03/26/2024] Open
Abstract
Plant collections held by botanic gardens and arboreta are key components of ex situ conservation. Maintaining genetic diversity in such collections allows them to be used as resources for supplementing wild populations. However, most recommended minimum sample sizes for sufficient ex situ genetic diversity are based on microsatellite markers, and it remains unknown whether these sample sizes remain valid in light of more recently developed next-generation sequencing (NGS) approaches. To address this knowledge gap, we examine how ex situ conservation status and sampling recommendations differ when derived from microsatellites and single nucleotide polymorphisms (SNPs) in garden and wild samples of two threatened oak species. For Quercus acerifolia, SNPs show lower ex situ representation of wild allelic diversity and slightly lower minimum sample size estimates than microsatellites, while results for each marker are largely similar for Q. boyntonii. The application of missing data filters tends to lead to higher ex situ representation, while the impact of different SNP calling approaches is dependent on the species being analyzed. Measures of population differentiation within species are broadly similar between markers, but larger numbers of SNP loci allow for greater resolution of population structure and clearer assignment of ex situ individuals to wild source populations. Our results offer guidance for future ex situ conservation assessments utilizing SNP data, such as the application of missing data filters and the usage of a reference genome, and illustrate that both microsatellites and SNPs remain viable options for botanic gardens and arboreta seeking to ensure the genetic diversity of their collections.
Collapse
Affiliation(s)
| | | | - Emma S. Spence
- Morton ArboretumCenter for Tree ScienceLisleIllinoisUSA
- Cornell UniversityDepartment of Public and Ecosystem HealthIthacaNew YorkUSA
| | - Sean M. Hoban
- Morton ArboretumCenter for Tree ScienceLisleIllinoisUSA
| |
Collapse
|
6
|
Ambu J, Martínez-Solano Í, Suchan T, Hernandez A, Wielstra B, Crochet PA, Dufresnes C. Genomic phylogeography illuminates deep cyto-nuclear discordances in midwife toads (Alytes). Mol Phylogenet Evol 2023; 183:107783. [PMID: 37044190 DOI: 10.1016/j.ympev.2023.107783] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2022] [Revised: 03/28/2023] [Accepted: 04/06/2023] [Indexed: 04/14/2023]
Abstract
The advent of genomic methods allows us to revisit the evolutionary history of organismal groups for which robust phylogenies are still lacking, particularly in species complexes that frequently hybridize. In this study, we conduct RAD-sequencing (RAD-seq) analyses of midwife toads (genus Alytes), an iconic group of western Mediterranean amphibians famous for their parental care behavior, but equally infamous for the difficulties to reconstruct their evolutionary history. Through admixture and phylogenetic analyses of thousands of loci, we provide the most comprehensive phylogeographic framework for the A. obstetricans complex to date, as well as the first fully resolved phylogeny for the entire genus. As part of this effort, we carefully explore the influence of different sampling schemes and data filtering thresholds on tree reconstruction, showing that several, slightly different, yet robust topologies may be retrieved with small datasets obtained by stringent SNP calling parameters, especially when admixed individuals are included. In contrast, analyses of incomplete but larger datasets converged on the same phylogeny, irrespective of the reconstruction method used or the proportion of missing data. The Alytes tree features three Miocene-diverged clades corresponding to the proposed subgenera Ammoryctis (A. cisternasii), Baleaphryne (A. maurus, A. dickhilleni and A. muletensis), and Alytes (A. obstetricans complex). The latter consists of six evolutionary lineages, grouped into three clades of Pliocene origin, and currently delimited as two species: (1) A. almogavarii almogavarii and A. a. inigoi; (2) A. obstetricans obstetricans and A. o. pertinax; (3) A. o. boscai and an undescribed taxon (A. o. cf. boscai). These results contradict the mitochondrial tree, due to past mitochondrial captures in A. a. almogavarii (central Pyrenees) and A. o. boscai (central Iberia) by A. obstetricans ancestors during the Pleistocene. Patterns of admixture between subspecies appear far more extensive than previously assumed from microsatellites, causing nomenclatural uncertainties, and even underlying the reticulate evolution of one taxon (A. o. pertinax). All Ammoryctis and Baleaphryne species form shallow clades, so their taxonomy should remain stable. Amid the prevalence of cyto-nuclear discordance among terrestrial vertebrates and the usual lack of resolution of conventional nuclear markers, our study advocates for phylogeography based on next-generation sequencing, but also encourages properly exploring parameter space and sampling schemes when building and analyzing genomic datasets.
Collapse
Affiliation(s)
- Johanna Ambu
- LASER, College of Biology and the Environment, Nanjing Forestry University, Nanjing, People's Republic of China
| | - Íñigo Martínez-Solano
- Departamento de Biodiversidad y Biología Evolutiva, Museo Nacional de Ciencias Naturales, CSIC, Madrid, Spain
| | - Tomasz Suchan
- W. Szafer Institute of Botany, Polish Academy of Sciences, Kraków, Poland
| | - Axel Hernandez
- LASER, College of Biology and the Environment, Nanjing Forestry University, Nanjing, People's Republic of China
| | - Ben Wielstra
- Naturalis Biodiversity Center, P.O. Box 9517, 2300 RA Leiden, The Netherlands; Institute of Biology Leiden, Leiden University, P.O. Box 9505, 2300 RA Leiden, The Netherlands
| | | | - Christophe Dufresnes
- LASER, College of Biology and the Environment, Nanjing Forestry University, Nanjing, People's Republic of China
| |
Collapse
|
7
|
Mohd Rodzik FF, Sudirman NA, Teh CK, Ong AL, Heng HY, Yaakop S, Mohd-Assaad N, Ong-Abdullah M, Ata N, Amit S, Saragih B, Appleton DR, Kulaveerasingam H. Development of Nuclear DNA Markers for Applications in Genetic Diversity Study of Oil Palm-Pollinating Weevil Populations. INSECTS 2023; 14:157. [PMID: 36835726 PMCID: PMC9967927 DOI: 10.3390/insects14020157] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/10/2022] [Revised: 12/19/2022] [Accepted: 12/20/2022] [Indexed: 06/18/2023]
Abstract
The oil palm-pollinating weevil (Elaeidobius kamerunicus Faust) was introduced from Cameroon, West Africa, to Malaysia in 1981, and subsequently, to other oil palm-growing countries as well. This study aims to develop a set of robust E. kamerunicus-specific nuclear DNA markers to directly assess the genetic diversity of the weevil populations. A total of 19,148 SNP and 223,200 SSR were discovered from 48 weevils representing three origins (Peninsular Malaysia, Sabah, and Riau) using RAD tag sequencing. Subsequent filtering steps further reduced these to 1000 SNP and 120 SSR. The selected 220 SNP exhibited a polymorphism information content (PIC) of 0.2387 (±0.1280), and 8 SSR had the PIC of 0.5084 (±0.1928). These markers were found to show sufficient polymorphism, making it possible to assign 180 weevils into three major clusters from Ghana, Cameroon, and Southeast Asia (mainly in Malaysia and Indonesia). These DNA markers successfully confirmed the Cameroon origin of the Southeast Asian cluster. However, the presence of null alleles in the SSR markers, due to limited flexibility of the probe design on the short RAD tags, led to an underestimation of heterozygosity within the populations. Hence, the developed SNP markers turned out to be more efficient than the SSR markers in the genetic diversity assessment of the E. kamerunicus populations. The genetic information provides useful insight into developing guidelines for the genetic monitoring and conservation planning of E. kamerunicus.
Collapse
Affiliation(s)
- Fairuz Farhana Mohd Rodzik
- Sime Darby Plantation R&D Centre, Serdang 43400, Selangor, Malaysia
- Centre for Insect Systematics, Department of Biological Science and Biotechnology, Faculty of Science and Technology, Universiti Kebangsaan Malaysia (UKM), Bangi 43600, Selangor, Malaysia
| | - Nurshazwani Amalina Sudirman
- Sime Darby Plantation R&D Centre, Serdang 43400, Selangor, Malaysia
- Centre for Insect Systematics, Department of Biological Science and Biotechnology, Faculty of Science and Technology, Universiti Kebangsaan Malaysia (UKM), Bangi 43600, Selangor, Malaysia
| | - Chee-Keng Teh
- Sime Darby Plantation R&D Centre, Serdang 43400, Selangor, Malaysia
| | - Ai-Ling Ong
- Sime Darby Plantation R&D Centre, Serdang 43400, Selangor, Malaysia
| | - Huey-Ying Heng
- Sime Darby Plantation R&D Centre, Serdang 43400, Selangor, Malaysia
| | - Salmah Yaakop
- Centre for Insect Systematics, Department of Biological Science and Biotechnology, Faculty of Science and Technology, Universiti Kebangsaan Malaysia (UKM), Bangi 43600, Selangor, Malaysia
| | - Norfarhan Mohd-Assaad
- Department of Applied Physics, Faculty of Science and Technology, Universiti Kebangsaan Malaysia (UKM), Bangi 43600, Selangor, Malaysia
- Institute of Systems Biology (INBIOSIS), Universiti Kebangsaan Malaysia, Bangi 43600, Selangor, Malaysia
| | - Meilina Ong-Abdullah
- Advanced Biotechnology and Breeding Centre, Malaysian Palm Oil Board (MPOB), 6 Persiaran Institusi, Bandar Baru Bangi, Kajang 43000, Selangor, Malaysia
| | - Nabeel Ata
- Advanced Biotechnology and Breeding Centre, Malaysian Palm Oil Board (MPOB), 6 Persiaran Institusi, Bandar Baru Bangi, Kajang 43000, Selangor, Malaysia
| | - Samsudin Amit
- Sime Darby Plantation R&D Centre, Serdang 43400, Selangor, Malaysia
| | - Burhanuddin Saragih
- Sime Darby Plantation R&D Centre, Serdang 43400, Selangor, Malaysia
- Minamas Research Centre Pekanbaru, Jalan Baru Bakal, Tualang Timur, Kecamatan Tualang, Kabupaten Siak, Perawang 28772, Provinsi Riau, Indonesia
| | | | | |
Collapse
|
8
|
Dufresnes C, Dutoit L, Brelsford A, Goldstein-Witsenburg F, Clément L, López-Baucells A, Palmeirim J, Pavlinić I, Scaravelli D, Ševčík M, Christe P, Goudet J. Inferring genetic structure when there is little: population genetics versus genomics of the threatened bat Miniopterus schreibersii across Europe. Sci Rep 2023; 13:1523. [PMID: 36707640 PMCID: PMC9883447 DOI: 10.1038/s41598-023-27988-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2022] [Accepted: 01/11/2023] [Indexed: 01/28/2023] Open
Abstract
Despite their paramount importance in molecular ecology and conservation, genetic diversity and structure remain challenging to quantify with traditional genotyping methods. Next-generation sequencing holds great promises, but this has not been properly tested in highly mobile species. In this article, we compared microsatellite and RAD-sequencing (RAD-seq) analyses to investigate population structure in the declining bent-winged bat (Miniopterus schreibersii) across Europe. Both markers retrieved general patterns of weak range-wide differentiation, little sex-biased dispersal, and strong isolation by distance that associated with significant genetic structure between the three Mediterranean Peninsulas, which could have acted as glacial refugia. Microsatellites proved uninformative in individual-based analyses, but the resolution offered by genomic SNPs illuminated on regional substructures within several countries, with colonies sharing migrators of distinct ancestry without admixture. This finding is consistent with a marked philopatry and spatial partitioning between mating and rearing grounds in the species, which was suspected from marked-recaptured data. Our study advocates that genomic data are necessary to properly unveil the genetic footprints left by biogeographic processes and social organization in long-distant flyers, which are otherwise rapidly blurred by their high levels of gene flow.
Collapse
Affiliation(s)
- Christophe Dufresnes
- Laboratory for Amphibian Systematic and Evolutionary Research, College of Biology and the Environment, Nanjing Forestry University, Nanjing, People's Republic of China.
| | - Ludovic Dutoit
- Department of Ecology and Evolution, University of Lausanne, 1015, Lausanne, Switzerland.,Department of Zoology, University of Otago, Dunedin, New Zealand
| | - Alan Brelsford
- Department of Evolution, Ecology, and Organismal Biology, University of California Riverside, Riverside, CA, USA
| | | | - Laura Clément
- Department of Ecology and Evolution, University of Lausanne, 1015, Lausanne, Switzerland
| | - Adria López-Baucells
- Bat Research Area, Granollers Museum of Natural Sciences, Carrer Palaudaries 102, 08402, Granollers, Spain
| | - Jorge Palmeirim
- Department of Animal Biology, Centre for Ecology, Evolution and Environmental Change - cE3c, University of Lisbon, 1749-016, Lisbon, Portugal
| | - Igor Pavlinić
- Department of Zoology, Croatian Natural History Museum, Demetrova 1, 10000, Zagreb, Croatia
| | - Dino Scaravelli
- Department of Biological, Geological, and Environmental Sciences, University of Bologna, Via Selmi 3, 40126, Bologna, Italy
| | - Martin Ševčík
- Department of Zoology, Faculty of Science, Charles University in Prague, Viničná 7, 128 44, Prague 2, Czech Republic
| | - Philippe Christe
- Department of Ecology and Evolution, University of Lausanne, 1015, Lausanne, Switzerland.
| | - Jérôme Goudet
- Department of Ecology and Evolution, University of Lausanne, 1015, Lausanne, Switzerland.
| |
Collapse
|
9
|
Byerly PA, Chesser RT, Fleischer RC, McInerney N, Przelomska NAS, Leberg PL. Museum Genomics Provide Evidence for Persistent Genetic Differentiation in a Threatened Seabird Species in the Western Atlantic. Integr Comp Biol 2022; 62:1838-1848. [PMID: 35781565 DOI: 10.1093/icb/icac107] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2022] [Revised: 06/18/2022] [Accepted: 06/27/2022] [Indexed: 01/05/2023] Open
Abstract
Connectivity among wildlife populations facilitates exchange of genetic material between groups. Changes to historical connectivity patterns resulting from anthropogenic activities can therefore have negative consequences for genetic diversity, particularly for small or isolated populations. DNA obtained from museum specimens can enable direct comparison of temporal changes in connectivity among populations, which can aid in conservation planning and contribute to the understanding of population declines. However, museum DNA can be degraded and only available in low quantities, rendering it challenging for use in population genomic analyses. Applications of genomic methodologies such as targeted sequencing address this issue by enabling capture of shared variable sites, increasing quantity and quality of recovered genomic information. We used targeted sequencing of ultra-conserved Elements (UCEs) to evaluate potential changes in connectivity and genetic diversity of roseate terns (Sterna dougallii) with a breeding distribution in the northwestern Atlantic and the Caribbean. Both populations experienced range contractions and population declines due to anthropogenic activity in the 20th century, which has the potential to alter historical connectivity regimes. Instead, we found that the two populations were differentiated historically as well as contemporaneously, with little evidence of migration between them for either time period. We also found no evidence for temporal changes in genetic diversity, although these interpretations may have been limited due to sequencing artifacts caused by the degraded nature of the museum samples. Population structuring in migratory seabirds is typically reflective of low rates of divergence and high connectivity among geographically segregated subpopulations. Our contrasting results suggest the potential presence of ecological mechanisms driving population differentiation, and highlight the value of targeted sequencing on DNA derived from museum specimens to uncover long-term patterns of genetic differentiation in wildlife populations.
Collapse
Affiliation(s)
- Paige A Byerly
- University of Louisiana at Lafayette, 104 E University Ave, Lafayette, LA 70504, USA.,Smithsonian's National Zoo and Conservation Biology Institute, 3001 Connecticut Avenue, NW, Washington, DC 20008, USA
| | - R Terry Chesser
- Eastern Ecological Science Center, U.S. Geological Survey, 12100 Beech Forest Road, Laurel, MD 20708, USA.,National Museum of Natural History, 10th St. and Constitution Avenue, NW, Washington, DC 20560, USA
| | - Robert C Fleischer
- Smithsonian's National Zoo and Conservation Biology Institute, 3001 Connecticut Avenue, NW, Washington, DC 20008, USA
| | - Nancy McInerney
- Smithsonian's National Zoo and Conservation Biology Institute, 3001 Connecticut Avenue, NW, Washington, DC 20008, USA
| | - Natalia A S Przelomska
- National Museum of Natural History, 10th St. and Constitution Avenue, NW, Washington, DC 20560, USA.,Smithsonian's National Zoo and Conservation Biology Institute, 3001 Connecticut Avenue, NW, Washington, DC 20008, USA.,Royal Botanic Gardens, Kew, Richmond TW9 3AE, UK
| | - Paul L Leberg
- University of Louisiana at Lafayette, 104 E University Ave, Lafayette, LA 70504, USA
| |
Collapse
|
10
|
Scherer BP, Mast A. Red Mangrove Propagule Bacterial Communities Vary With Geographic, But Not Genetic Distance. MICROBIAL ECOLOGY 2022:10.1007/s00248-022-02147-w. [PMID: 36441249 DOI: 10.1007/s00248-022-02147-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/14/2022] [Accepted: 11/12/2022] [Indexed: 06/16/2023]
Abstract
Bacterial communities associated with plant propagules remain understudied, despite the opportunities that propagules represent as dispersal vectors for bacteria to new sites. These communities may be the product of a combination of environmental influence and inheritance from parent to offspring. The relative role of these mechanisms could have significant implications for our understanding of plant-microbe interactions. We studied the correlates of microbiome community similarities across an invasion front of red mangroves (Rhizophora mangle L.) in Florida, where the species is expanding northward. We collected georeferenced propagule samples from 110 individuals of red mangroves across 11 populations in Florida and used 16S rRNA gene (iTag) sequencing to describe their bacterial communities. We found no core community of bacterial amplicon sequence variants (ASVs) across the Florida range of red mangroves, though there were some ASVs shared among individuals within most populations. Populations differed significantly as measured by Bray-Curtis dissimilarity, but not Unifrac distance. We generated data from 6 microsatellite loci from 60 individuals across 9 of the 11 populations. Geographic distance was correlated with beta diversity, but genetic distance was not. We conclude that red mangrove propagule bacterial communities are likely influenced more by local environmental acquisition than by inheritance.
Collapse
Affiliation(s)
- Brendan P Scherer
- Department of Biological Science, Florida State University, King Life Sciences Building, 319 Stadium Drive, Tallahassee, Fl, 32304, USA.
| | - Austin Mast
- Department of Biological Science, Florida State University, King Life Sciences Building, 319 Stadium Drive, Tallahassee, Fl, 32304, USA
| |
Collapse
|
11
|
Cohen JI, Ruane LG. Conservation genetics of Phlox hirsuta, a serpentine endemic. CONSERV GENET 2022. [DOI: 10.1007/s10592-022-01478-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
12
|
Halsey MK, Stuhler JD, Bayona-Vásquez NJ, Platt RN, Goetze JR, Martin RE, Matocha KG, Bradley RD, Stevens RD, Ray DA. Comparison of genetic variation between rare and common congeners of Dipodomys with estimates of contemporary and historical effective population size. PLoS One 2022; 17:e0274554. [PMID: 36099283 PMCID: PMC9469943 DOI: 10.1371/journal.pone.0274554] [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: 01/27/2022] [Accepted: 08/31/2022] [Indexed: 11/18/2022] Open
Abstract
Species with low effective population sizes are at greater risk of extinction because of reduced genetic diversity. Such species are more vulnerable to chance events that decrease population sizes (e.g. demographic stochasticity). Dipodomys elator, (Texas kangaroo rat) is a kangaroo rat that is classified as threatened in Texas and field surveys from the past 50 years indicate that the distribution of this species has decreased. This suggests geographic range reductions that could have caused population fluctuations, potentially impacting effective population size. Conversely, the more common and widespread D. ordii (Ord’s kangaroo rat) is thought to exhibit relative geographic and demographic stability. We assessed the genetic variation of D. elator and D. ordii samples using 3RAD, a modified restriction site associated sequencing approach. We hypothesized that D. elator would show lower levels of nucleotide diversity, observed heterozygosity, and effective population size when compared to D. ordii. We were also interested in identifying population structure within contemporary samples of D. elator and detecting genetic variation between temporal samples to understand demographic dynamics. We analyzed up to 61,000 single nucleotide polymorphisms. We found that genetic variability and effective population size in contemporary D. elator populations is lower than that of D. ordii. There is slight, if any, population structure within contemporary D. elator samples, and we found low genetic differentiation between spatial or temporal historical samples. This indicates little change in nuclear genetic diversity over 30 years. Results suggest that genetic diversity of D. elator has remained stable despite reduced population size and/or abundance, which may indicate a metapopulation-like system, whose fluctuations might counteract species extinction.
Collapse
Affiliation(s)
- Michaela K. Halsey
- Department of Biological Sciences, Texas Tech University, Lubbock, Texas, United States of America
- Department of Natural Resources Management, Texas Tech University, Lubbock, Texas, United States of America
| | - John D. Stuhler
- Department of Natural Resources Management, Texas Tech University, Lubbock, Texas, United States of America
| | - Natalia J. Bayona-Vásquez
- Department of Environmental Health Science, University of Georgia, Athens, Georgia, United States of America
- Institute of Bioinformatics, University of Georgia, Athens, Georgia, United States of America
| | - Roy N. Platt
- Texas Biomedical Research Institute, San Antonio, Texas, United States of America
| | - Jim R. Goetze
- Natural Sciences Department, Laredo College, Laredo, Texas, United States of America
| | - Robert E. Martin
- Department of Biology, McMurry University, Abilene, Texas, United States of America
| | - Kenneth G. Matocha
- Department of Biology, South Arkansas Community College, El Dorado, Arkansas, United States of America
| | - Robert D. Bradley
- Department of Biological Sciences, Texas Tech University, Lubbock, Texas, United States of America
- Natural Science Research Laboratory, Museum of Texas Tech, Lubbock, Texas, United States of America
| | - Richard D. Stevens
- Department of Natural Resources Management, Texas Tech University, Lubbock, Texas, United States of America
- Natural Science Research Laboratory, Museum of Texas Tech, Lubbock, Texas, United States of America
| | - David A. Ray
- Department of Biological Sciences, Texas Tech University, Lubbock, Texas, United States of America
- * E-mail:
| |
Collapse
|
13
|
Vitek NS, McDaniel SF, Bloch JI. Microevolutionary variation in molar morphology of Onychomys leucogaster decoupled from genetic structure. Evolution 2022; 76:2032-2048. [PMID: 35872621 DOI: 10.1111/evo.14576] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2021] [Revised: 04/22/2022] [Accepted: 04/29/2022] [Indexed: 01/22/2023]
Abstract
In neutral models of quantitative trait evolution, both genetic and phenotypic divergence scale as random walks, producing a correlation between the two measures. However, complexity in the genotype-phenotype map may alter the correlation between genotypic and phenotypic divergence, even when both are evolving neutrally or nearly so. Understanding this correlation between phenotypic and genetic variation is critical for accurately interpreting the fossil record. This study compares the geographic structure and scaling of morphological variation of the shape of the first lower molar of 77 individuals of the northern grasshopper mouse Onychomys leucogaster to genome-wide SNP variation in the same sample. We found strong genetic structure but weak or absent morphological structure indicating that the scaling of each type of variation is decoupled from one another. Low PST values relative to FST values are consistent with a lack of morphological divergence in contrast to genetic divergence between groups. This lack of phenotypic structure and the presence of notable within-sample phenotypic variance are consistent with uniform selection or constraints on molar shape across a wide geographic and environmental range. Over time, this kind of decoupling may result in patterns of phenotypic stasis masking underlying genetic patterns.
Collapse
Affiliation(s)
- Natasha S Vitek
- Department of Biology, University of Florida, Gainesville, Florida, 32611.,Florida Museum of Natural History, University of Florida, Gainesville, Florida, 32611.,Department of Ecology and Evolution, Stony Brook University, Stony Brook, New York, 11794
| | - Stuart F McDaniel
- Department of Biology, University of Florida, Gainesville, Florida, 32611
| | - Jonathan I Bloch
- Florida Museum of Natural History, University of Florida, Gainesville, Florida, 32611
| |
Collapse
|
14
|
Aylward M, Sagar V, Natesh M, Ramakrishnan U. How methodological changes have influenced our understanding of population structure in threatened species: insights from tiger populations across India. Philos Trans R Soc Lond B Biol Sci 2022; 377:20200418. [PMID: 35430878 PMCID: PMC9014192 DOI: 10.1098/rstb.2020.0418] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
Unprecedented advances in sequencing technology in the past decade allow a better understanding of genetic variation and its partitioning in natural populations. Such inference is critical to conservation: to understand species biology and identify isolated populations. We review empirical population genetics studies of Endangered Bengal tigers within India, where 60–70% of wild tigers live. We assess how changes in marker type and sampling strategy have impacted inferences by reviewing past studies, and presenting three novel analyses including a single-nucleotide polymorphism (SNP) panel, genome-wide SNP markers, and a whole-mitochondrial genome network. At a broad spatial scale, less than 100 SNPs revealed the same patterns of population clustering as whole genomes (with the exception of one additional population sampled only in the SNP panel). Mitochondrial DNA indicates a strong structure between the northeast and other regions. Two studies with more populations sampled revealed further substructure within Central India. Overall, the comparison of studies with varied marker types and sample sets allows more rigorous inference of population structure. Yet sampling of some populations is limited across all studies, and these should be the focus of future sampling efforts. We discuss challenges in our understanding of population structure, and how to further address relevant questions in conservation genetics.
This article is part of the theme issue ‘Celebrating 50 years since Lewontin's apportionment of human diversity’.
Collapse
Affiliation(s)
- Megan Aylward
- National Centre for Biological Sciences, TIFR, Bangalore, India, 560065
| | - Vinay Sagar
- National Centre for Biological Sciences, TIFR, Bangalore, India, 560065
| | - Meghana Natesh
- Indian Institute of Science Education and Research, Tirupati, India, 517507
| | - Uma Ramakrishnan
- National Centre for Biological Sciences, TIFR, Bangalore, India, 560065
- Senior Fellow, DBT Wellcome Trust India Alliance, Hyderabad, Telangana, India, 500034
| |
Collapse
|
15
|
Mosher CM, Johnson CJ, Murray BW. Reduced genetic diversity associated with the northern expansion of an amphibian species with high habitat specialization,
Ascaphus truei
, resolved using two types of genetic markers. Ecol Evol 2022; 12:e8716. [PMID: 35342604 PMCID: PMC8931771 DOI: 10.1002/ece3.8716] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2021] [Revised: 01/18/2022] [Accepted: 01/24/2022] [Indexed: 11/13/2022] Open
Abstract
Reconstruction of historical relationships between geographic regions within a species’ range can indicate dispersal patterns and help predict future responses to shifts in climate. Ascaphus truei (coastal tailed frog) is an indicator species of the health of forests and perennial streams in the Coastal and Cascade Mountains of the Pacific Northwest of North America. We used two genetic techniques—microsatellite and genotype‐by‐sequencing (GBS)—to compare the within‐region genetic diversity of populations near the northern extent of the species’ range (British Columbia, Canada) to two geographic regions in British Columbia and two in Washington, USA, moving toward the core of the range. Allelic richness and heterozygosity declined substantially as latitude increased. The northernmost region had the lowest mean expected heterozygosities for both techniques (microsatellite, M = 0.20, SE = 0.080; GBS, M = 0.025, SE = 0.0010) and the southernmost region had the highest (microsatellite, M = 0.88, SE = 0.054; GBS, M = 0.20, SE = 0.0029). The northernmost regions (NC and MC) clustered together in population structure models for both genetic techniques. Our discovery of reduced diversity may have important conservation and management implications for population connectivity and the response of A. truei to climate change.
Collapse
Affiliation(s)
- Cherie M. Mosher
- University of Northern British Columbia Prince George British Columbia Canada
| | - Chris J. Johnson
- University of Northern British Columbia Prince George British Columbia Canada
| | - Brent W. Murray
- University of Northern British Columbia Prince George British Columbia Canada
| |
Collapse
|
16
|
Linan AG, Lowry Ii PP, Miller AJ, Schatz GE, Sevathian JC, Edwards CE. Interspecific hybridization and island colonization history, not rarity, most strongly affect the genetic diversity in a clade of Mascarene-endemic trees. J Hered 2022; 113:336-352. [PMID: 35192705 DOI: 10.1093/jhered/esac010] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2021] [Accepted: 02/21/2022] [Indexed: 11/14/2022] Open
Abstract
Many factors shape the genetic diversity of island-endemic trees, with important implications for conservation. Oceanic island-endemic lineages undergo an initial founding bottleneck during the colonization process and subsequently accumulate diversity following colonization. Moreover, many island endemics occur in small populations and are further threatened by anthropogenic factors that cause population declines, making them susceptible to losses in genetic diversity through genetic drift, inbreeding, and bottlenecks. However, life-history traits commonly found in trees, such as outcrossing mechanisms, long lifespans, and a propensity for interspecific hybridization, may help buffer against losses of genetic variation. To assess the relative importance of colonization history, rarity, and distribution in shaping genetic diversity of island-endemic trees, we conducted a comparative population genomic analysis of 13 species of Diospyros (Ebenaceae) endemic to the Mascarene Islands that differ in island colonization history, distribution, population size, and IUCN threat status. We genotyped 328 individuals across the islands using 2b-RADseq, compared genetic diversity both among and within species, and assessed patterns of genetic structure. Genetic diversity did not vary significantly by IUCN status, but we found that species that co-occur with others on the same intermediate-aged island (Mauritius) had much greater genetic diversity than those that occur solitarily on an island (Réunion and Rodrigues), likely because of greater interspecific hybridization among species with overlapping distributions and processes related to time since island colonization. Results presented here were used to determine priority localities for in situ and ex situ conservation efforts to maximize the genetic diversity of each Mascarene Diospyros species.
Collapse
Affiliation(s)
- Alexander G Linan
- Center for Conservation and Sustainable Development, Missouri Botanical Garden, 4344 Shaw Blvd., St. Louis, MO 63110, USA.,Department of Biology, Saint Louis University, Macelwane Hall, 3507 Laclede Avenue, St. Louis, MO 63103, USA
| | - Porter P Lowry Ii
- Africa and Madagascar Program, Missouri Botanical Garden, 4344 Shaw Blvd., St. Louis, MO 63110, USA.,Institut de Systématique, Évolution et Biodiversité (ISYEB), Muséum National d'Histoire Naturelle, Centre National de la Recherche Scientifique, Sorbonne Université, École Pratique des Hautes Études, Université des Antilles, C.P. 39, 57 rue Cuvier, 75005 Paris, France
| | - Allison J Miller
- Department of Biology, Saint Louis University, Macelwane Hall, 3507 Laclede Avenue, St. Louis, MO 63103, USA.,Donald Danforth Plant Science Center, 975 North Warson Road, St. Louis, MO 63132, USA
| | - George E Schatz
- Africa and Madagascar Program, Missouri Botanical Garden, 4344 Shaw Blvd., St. Louis, MO 63110, USA
| | - Jean Claude Sevathian
- Botanist, Sustainability Consultant and Landscape Care and Maintenance Service, 20 Bis Morcellement La Confiance, Beau Bassin 71504, Mauritius
| | - Christine E Edwards
- Center for Conservation and Sustainable Development, Missouri Botanical Garden, 4344 Shaw Blvd., St. Louis, MO 63110, USA
| |
Collapse
|
17
|
de Lima Ferreira P, Batista R, Andermann T, Groppo M, Bacon CD, Antonelli A. Target sequence capture of Barnadesioideae (Compositae) demonstrates the utility of low coverage loci in phylogenomic analyses. Mol Phylogenet Evol 2022; 169:107432. [DOI: 10.1016/j.ympev.2022.107432] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2021] [Revised: 12/21/2021] [Accepted: 01/14/2022] [Indexed: 11/26/2022]
|
18
|
Conservation genomics of an Australian cycad Cycas calcicola, and the Absence of Key Genotypes in Botanic Gardens. CONSERV GENET 2022. [DOI: 10.1007/s10592-022-01428-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
|
19
|
Han J, Munro JE, Kocoski A, Barry AE, Bahlo M. Population-level genome-wide STR discovery and validation for population structure and genetic diversity assessment of Plasmodium species. PLoS Genet 2022; 18:e1009604. [PMID: 35007277 PMCID: PMC8782505 DOI: 10.1371/journal.pgen.1009604] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2021] [Revised: 01/21/2022] [Accepted: 12/14/2021] [Indexed: 11/18/2022] Open
Abstract
Short tandem repeats (STRs) are highly informative genetic markers that have been used extensively in population genetics analysis. They are an important source of genetic diversity and can also have functional impact. Despite the availability of bioinformatic methods that permit large-scale genome-wide genotyping of STRs from whole genome sequencing data, they have not previously been applied to sequencing data from large collections of malaria parasite field samples. Here, we have genotyped STRs using HipSTR in more than 3,000 Plasmodium falciparum and 174 Plasmodium vivax published whole-genome sequence data from samples collected across the globe. High levels of noise and variability in the resultant callset necessitated the development of a novel method for quality control of STR genotype calls. A set of high-quality STR loci (6,768 from P. falciparum and 3,496 from P. vivax) were used to study Plasmodium genetic diversity, population structures and genomic signatures of selection and these were compared to genome-wide single nucleotide polymorphism (SNP) genotyping data. In addition, the genome-wide information about genetic variation and other characteristics of STRs in P. falciparum and P. vivax have been available in an interactive web-based R Shiny application PlasmoSTR (https://github.com/bahlolab/PlasmoSTR).
Collapse
Affiliation(s)
- Jiru Han
- Population Health and Immunity Division, The Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
- Department of Medical Biology, The University of Melbourne, Melbourne, Australia
| | - Jacob E. Munro
- Population Health and Immunity Division, The Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
- Department of Medical Biology, The University of Melbourne, Melbourne, Australia
| | - Anthony Kocoski
- Population Health and Immunity Division, The Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
- Department of Mathematics and Statistics, The University of Melbourne, Melbourne, Australia
| | - Alyssa E. Barry
- Population Health and Immunity Division, The Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
- Department of Medical Biology, The University of Melbourne, Melbourne, Australia
- Disease Elimination Program, Burnet Institute, Melbourne, Australia
- IMPACT Institute for Innovation in Mental and Physical Health and Clinical Translation, Deakin University, Geelong, Australia
| | - Melanie Bahlo
- Population Health and Immunity Division, The Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
- Department of Medical Biology, The University of Melbourne, Melbourne, Australia
- * E-mail:
| |
Collapse
|
20
|
Jerney J, Rengefors K, Nagai S, Krock B, Sjöqvist C, Suikkanen S, Kremp A. Seasonal genotype dynamics of a marine dinoflagellate: Pelagic populations are homogeneous and as diverse as benthic seed banks. Mol Ecol 2021; 31:512-528. [PMID: 34716943 PMCID: PMC9298838 DOI: 10.1111/mec.16257] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2021] [Revised: 10/19/2021] [Accepted: 10/25/2021] [Indexed: 11/27/2022]
Abstract
Genetic diversity is the basis for evolutionary adaptation and selection under changing environmental conditions. Phytoplankton populations are genotypically diverse, can become genetically differentiated within small spatiotemporal scales and many species form resting stages. Resting stage accumulations in sediments (seed banks) are expected to serve as reservoirs for genetic information, but so far their role in maintaining phytoplankton diversity and in evolution has remained unclear. In this study we used the toxic dinoflagellate Alexandrium ostenfeldii (Dinophyceae) as a model organism to investigate if (i) the benthic seed bank is more diverse than the pelagic population and (ii) the pelagic population is seasonally differentiated. Resting stages (benthic) and plankton (pelagic) samples were collected at a coastal bloom site in the Baltic Sea, followed by cell isolation and genotyping using microsatellite markers (MS) and restriction site associated DNA sequencing (RAD). High clonal diversity (98%–100%) combined with intermediate to low gene diversity (0.58–0.03, depending on the marker) was found. Surprisingly, the benthic and pelagic fractions of the population were equally diverse, and the pelagic fraction was temporally homogeneous, despite seasonal fluctuation of environmental selection pressures. The results of this study suggest that continuous benthic–pelagic coupling, combined with frequent sexual reproduction, as indicated by persistent linkage equilibrium, prevent the dominance of single clonal lineages in a dynamic environment. Both processes harmonize the pelagic with the benthic population and thus prevent seasonal population differentiation. At the same time, frequent sexual reproduction and benthic–pelagic coupling maintain high clonal diversity in both habitats.
Collapse
Affiliation(s)
- Jacqueline Jerney
- Tvärminne Zoological Station, University of Helsinki, Hanko, Finland.,Marine Research Center, Finnish Environment Institute, Helsinki, Finland
| | | | - Satoshi Nagai
- National Research Institute of Fisheries Science, Yokohama, Kanagawa, Japan
| | - Bernd Krock
- Alfred-Wegener-Institut Helmholtz-Zentrum für Polar- und Meeresforschung, Bremerhaven, Germany
| | - Conny Sjöqvist
- Faculty of Science and Engineering, Environmental and Marine Biology, Åbo Akademi University, Turku, Finland
| | - Sanna Suikkanen
- Marine Research Center, Finnish Environment Institute, Helsinki, Finland
| | - Anke Kremp
- Marine Research Center, Finnish Environment Institute, Helsinki, Finland
| |
Collapse
|
21
|
Ortiz D, Pekár S, Dianat M. Phylogenomics and loci dropout patterns of deeply diverged Zodarion ant-eating spiders suggest a high potential of RAD-seq for genus-level spider phylogenetics. Cladistics 2021; 38:320-334. [PMID: 34699083 DOI: 10.1111/cla.12493] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 10/02/2021] [Indexed: 11/28/2022] Open
Abstract
RAD sequencing yields large amounts of genome-wide data at a relatively low cost and without requiring previous taxon-specific information, making it ideal for evolutionary studies of highly diversified and neglected organisms. However, concerns about information decay with phylogenetic distance have discouraged its use for assessing supraspecific relationships. Here, using Double Digest Restriction Associated DNA (ddRAD) data, we perform the first deep-level approach to the phylogeny of Zodarion, a highly diversified spider genus. We explore the impact of loci and taxon filtering across concatenated and multispecies coalescent reconstruction methods and investigate the patterns of information dropout in reference to both the time of divergence and the mitochondrial divergence between taxa. We found that relaxed loci-filtering and nested taxon-filtering strategies maximized the amount of molecular information and improved phylogenetic inference. As expected, there was a clear pattern of allele dropout towards deeper time and mitochondrial divergences, but the phylogenetic signal remained strong throughout the phylogeny. Therefore, we inferred topologies that were almost fully resolved, highly supported, and noticeably congruent between setups and inference methods, which highlights overall inconsistency in the taxonomy of Zodarion. Because Zodarion appears to be among the oldest and most mitochondrially diversified spider genera, our results suggest that ddRAD data show high potential for inferring intra-generic relationships across spiders and probably also in other taxonomic groups.
Collapse
Affiliation(s)
- David Ortiz
- Department of Botany and Zoology, Faculty of Science, Masaryk University, Kotlářská 2, 611 37 Brno, Czechia
| | - Stano Pekár
- Department of Botany and Zoology, Faculty of Science, Masaryk University, Kotlářská 2, 611 37 Brno, Czechia
| | - Malahat Dianat
- Department of Botany and Zoology, Faculty of Science, Masaryk University, Kotlářská 2, 611 37 Brno, Czechia
| |
Collapse
|
22
|
Nazareno AG, Knowles LL. There Is No 'Rule of Thumb': Genomic Filter Settings for a Small Plant Population to Obtain Unbiased Gene Flow Estimates. FRONTIERS IN PLANT SCIENCE 2021; 12:677009. [PMID: 34721447 PMCID: PMC8551369 DOI: 10.3389/fpls.2021.677009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/07/2021] [Accepted: 06/16/2021] [Indexed: 06/13/2023]
Abstract
The application of high-density polymorphic single-nucleotide polymorphisms (SNP) markers derived from high-throughput sequencing methods has heralded plenty of biological questions about the linkages of processes operating at micro- and macroevolutionary scales. However, the effects of SNP filtering practices on population genetic inference have received much less attention. By performing sensitivity analyses, we empirically investigated how decisions about the percentage of missing data (MD) and the minor allele frequency (MAF) set in bioinformatic processing of genomic data affect direct (i.e., parentage analysis) and indirect (i.e., fine-scale spatial genetic structure - SGS) gene flow estimates. We focus specifically on these manifestations in small plant populations, and particularly, in the rare tropical plant species Dinizia jueirana-facao, where assumptions implicit to analytical procedures for accurate estimates of gene flow may not hold. Avoiding biases in dispersal estimates are essential given this species is facing extinction risks due to habitat loss, and so we also investigate the effects of forest fragmentation on the accuracy of dispersal estimates under different filtering criteria by testing for recent decrease in the scale of gene flow. Our sensitivity analyses demonstrate that gene flow estimates are robust to different setting of MAF (0.05-0.35) and MD (0-20%). Comparing the direct and indirect estimates of dispersal, we find that contemporary estimates of gene dispersal distance (σ r t = 41.8 m) was ∼ fourfold smaller than the historical estimates, supporting the hypothesis of a temporal shift in the scale of gene flow in D. jueirana-facao, which is consistent with predictions based on recent, dramatic forest fragmentation process. While we identified settings for filtering genomic data to avoid biases in gene flow estimates, we stress that there is no 'rule of thumb' for bioinformatic filtering and that relying on default program settings is not advisable. Instead, we suggest that the approach implemented here be applied independently in each separate empirical study to confirm appropriate settings to obtain unbiased population genetics estimates.
Collapse
Affiliation(s)
- Alison G. Nazareno
- Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI, United States
- Department of Genetics, Ecology and Evolution, Federal University of Minas Gerais, Belo Horizonte, Brazil
| | - L. Lacey Knowles
- Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI, United States
| |
Collapse
|
23
|
Suyama Y, Hirota SK, Matsuo A, Tsunamoto Y, Mitsuyuki C, Shimura A, Okano K. Complementary combination of multiplex high‐throughput
DNA
sequencing for molecular phylogeny. Ecol Res 2021. [DOI: 10.1111/1440-1703.12270] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Yoshihisa Suyama
- Graduate School of Agricultural Science Tohoku University Osaki Miyagi Japan
| | - Shun K. Hirota
- Graduate School of Agricultural Science Tohoku University Osaki Miyagi Japan
| | - Ayumi Matsuo
- Graduate School of Agricultural Science Tohoku University Osaki Miyagi Japan
| | - Yoshihiro Tsunamoto
- Graduate School of Agricultural Science Tohoku University Osaki Miyagi Japan
| | - Chika Mitsuyuki
- Graduate School of Agricultural Science Tohoku University Osaki Miyagi Japan
| | - Atsuki Shimura
- Graduate School of Agricultural Science Tohoku University Osaki Miyagi Japan
| | - Kunihiro Okano
- Department of Biological Environment Akita Prefectural University Akita Japan
| |
Collapse
|
24
|
Yi X, Latch EK. Nonrandom missing data can bias Principal Component Analysis inference of population genetic structure. Mol Ecol Resour 2021; 22:602-611. [PMID: 34463035 DOI: 10.1111/1755-0998.13498] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2020] [Revised: 08/17/2021] [Accepted: 08/26/2021] [Indexed: 11/29/2022]
Abstract
Population genetic studies in non-model systems increasingly use next-generation sequencing to obtain more loci, but such methods also generate more missing data that may affect downstream analyses. Here we focus on the principal component analysis (PCA) which has been widely used to explore and visualize population structure with mean-imputed missing data. We simulated data of different population models with various total missingness (1%, 10%, 20%) introduced either randomly or biased among individuals or populations. We found that individuals biased with missing data would be dragged away from their real population clusters to the origin of PCA plots, making them indistinguishable from true admixed individuals and potentially leading to misinterpreted population structure. We also generated empirical data of the big brown bat (Eptesicus fuscus) using restriction site-associated DNA sequencing (RADseq). We filtered three data sets with 19.12%, 9.87%, and 1.35% total missingness, all showing nonrandom missing data with biased individuals dragged towards the PCA origin, consistent with results from simulations. We highlight the importance of considering missing data effects on PCA in non-model systems where nonrandom missing data are common due to varying sample quality. To help detect missing data effects, we suggest to (1) plot PCA with a colour gradient showing per sample missingness, (2) interpret samples close to the PCA origin with extra caution, (3) explore filtering parameters with and without the missingness-biased samples, and (4) use complementary analyses (e.g., model-based methods) to cross-validate PCA results and help interpret population structure.
Collapse
Affiliation(s)
- Xueling Yi
- Behavioral and Molecular Ecology Research Group, Department of Biological Sciences, University of Wisconsin-Milwaukee, Milwaukee, Wisconsin, USA
| | - Emily K Latch
- Behavioral and Molecular Ecology Research Group, Department of Biological Sciences, University of Wisconsin-Milwaukee, Milwaukee, Wisconsin, USA
| |
Collapse
|
25
|
Leal BSS, Chaves CJN, Graciano VA, Boury C, Huacre LAP, Heuertz M, Palma-Silva C. Evidence of local adaptation despite strong drift in a Neotropical patchily distributed bromeliad. Heredity (Edinb) 2021; 127:203-218. [PMID: 33953353 PMCID: PMC8322333 DOI: 10.1038/s41437-021-00442-9] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2020] [Revised: 04/16/2021] [Accepted: 04/17/2021] [Indexed: 02/03/2023] Open
Abstract
Both genetic drift and divergent selection are predicted to be drivers of population differentiation across patchy habitats, but the extent to which these forces act on natural populations to shape traits is strongly affected by species' ecological features. In this study, we infer the genomic structure of Pitcairnia lanuginosa, a widespread herbaceous perennial plant with a patchy distribution. We sampled populations in the Brazilian Cerrado and the Central Andean Yungas and discovered and genotyped SNP markers using double-digest restriction-site associated DNA sequencing. In addition, we analyzed ecophysiological traits obtained from a common garden experiment and compared patterns of phenotypic and genetic divergence (PST-FST comparisons) in a subset of populations from the Cerrado. Our results from molecular analyses pointed to extremely low genetic diversity and a remarkable population differentiation, supporting a major role of genetic drift. Approximately 0.3% of genotyped SNPs were flagged as differentiation outliers by at least two distinct methods, and Bayesian generalized linear mixed models revealed a signature of isolation by environment in addition to isolation by distance for high-differentiation outlier SNPs among the Cerrado populations. PST-FST comparisons suggested divergent selection on two ecophysiological traits linked to drought tolerance. We showed that these traits vary among populations, although without any particular macro-spatial pattern, suggesting local adaptation to differences in micro-habitats. Our study shows that selection might be a relevant force, particularly for traits involved in drought stress, even for populations experiencing strong drift, which improves our knowledge on eco-evolutionary processes acting on non-continuously distributed species.
Collapse
Affiliation(s)
- Bárbara Simões Santos Leal
- grid.410543.70000 0001 2188 478XDepartamento de Ecologia, Instituto de Biociências, Universidade Estadual Paulista, Rio Claro, São Paulo Brazil
| | - Cleber Juliano Neves Chaves
- grid.410543.70000 0001 2188 478XDepartamento de Ecologia, Instituto de Biociências, Universidade Estadual Paulista, Rio Claro, São Paulo Brazil
| | - Vanessa Araujo Graciano
- grid.410543.70000 0001 2188 478XDepartamento de Ecologia, Instituto de Biociências, Universidade Estadual Paulista, Rio Claro, São Paulo Brazil
| | - Christophe Boury
- grid.412041.20000 0001 2106 639XINRAE, Univ. Bordeaux, Biogeco, Cestas France
| | - Luis Alberto Pillaca Huacre
- grid.10800.390000 0001 2107 4576Departamento de Ecología, Museo de Historia Natural de la Universidad Nacional Mayor de San Marcos, Lima, Peru
| | - Myriam Heuertz
- grid.412041.20000 0001 2106 639XINRAE, Univ. Bordeaux, Biogeco, Cestas France
| | - Clarisse Palma-Silva
- grid.410543.70000 0001 2188 478XDepartamento de Ecologia, Instituto de Biociências, Universidade Estadual Paulista, Rio Claro, São Paulo Brazil ,grid.411087.b0000 0001 0723 2494Departamento de Biologia Vegetal, Instituto de Biologia, Universidade Estadual de Campinas, Campinas, São Paulo, Brazil
| |
Collapse
|
26
|
How well do genetic markers inform about responses to intraspecific admixture? A comparative analysis of microsatellites and RADseq. BMC Genom Data 2021; 22:22. [PMID: 34182923 PMCID: PMC8237422 DOI: 10.1186/s12863-021-00974-3] [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: 02/17/2021] [Accepted: 05/20/2021] [Indexed: 11/21/2022] Open
Abstract
Background Fitness consequences of intraspecific genetic admixture can vary from positive to negative depending on the genetic composition of the populations and environmental conditions. Because admixture has potential to influence the success of management and conservation efforts, genetic similarity has been suggested to be used as a proxy to predict the outcome. Studies utilizing microsatellites (a neutral marker) to investigate associations between genetic distance and admixture effects show conflicting results. Marker types that yield information on genome-wide and/or adaptive variation might be more useful for predicting responses to inter-population hybridization. In this study we utilized published data for three populations of pike (Esox lucius) to investigate associations between offspring performance (hatching success) and parental genetic similarity in experimentally purebred and admixed families, based on neutral (microsatellites), genome-wide neutral (RADseq SNPs), and adaptive (SNPs under selection) markers. Results Estimated similarity varied among the markers, likely reflecting differences in their inherent properties, but was consistently higher in purebred than admixed families. A significant interaction between marker type and admixture treatment reflected that neutral SNPs yielded higher estimates than adaptive SNPs for admixed families whereas no difference was found for purebred families, which indicates that neutral similarity was not reflective of adaptive similarity. When all samples were pooled, no association between similarity and performance was found for any marker. For microsatellites, similarity was positively correlated with hatching success in purebred families, whereas no association was found in admixed families; however, the direction of the effect differed between the population combinations. Conclusions The results strengthen the notion that, as of today, there is no proxy that can reliably predicted the outcome of admixture. This emphasizes the need of further studies to advance knowledge that can shed light on how to safeguard against negative consequences of admixture, and thereby inform management and promote conservation of biological diversity. Supplementary Information The online version contains supplementary material available at 10.1186/s12863-021-00974-3.
Collapse
|
27
|
Chafin TK, Zbinden ZD, Douglas MR, Martin BT, Middaugh CR, Gray MC, Ballard JR, Douglas ME. Spatial population genetics in heavily managed species: Separating patterns of historical translocation from contemporary gene flow in white-tailed deer. Evol Appl 2021; 14:1673-1689. [PMID: 34178112 PMCID: PMC8210790 DOI: 10.1111/eva.13233] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2020] [Accepted: 03/10/2021] [Indexed: 01/16/2023] Open
Abstract
Approximately 100 years ago, unregulated harvest nearly eliminated white-tailed deer (Odocoileus virginianus) from eastern North America, which subsequently served to catalyze wildlife management as a national priority. An extensive stock-replenishment effort soon followed, with deer broadly translocated among states as a means of re-establishment. However, an unintended consequence was that natural patterns of gene flow became obscured and pretranslocation signatures of population structure were replaced. We applied cutting-edge molecular and biogeographic tools to disentangle genetic signatures of historical management from those reflecting spatially heterogeneous dispersal by evaluating 35,099 single nucleotide polymorphisms (SNPs) derived via reduced-representation genomic sequencing from 1143 deer sampled statewide in Arkansas. We then employed Simpson's diversity index to summarize ancestry assignments and visualize spatial genetic transitions. Using sub-sampled transects across these transitions, we tested clinal patterns across loci against theoretical expectations of their response under scenarios of re-colonization and restricted dispersal. Two salient results emerged: (A) Genetic signatures from historic translocations are demonstrably apparent; and (B) Geographic filters (major rivers; urban centers; highways) now act as inflection points for the distribution of this contemporary ancestry. These results yielded a statewide assessment of contemporary population structure in deer as driven by historic translocations as well as ongoing processes. In addition, the analytical framework employed herein to effectively decipher extant/historic drivers of deer distribution in Arkansas is also applicable for other biodiversity elements with similarly complex demographic histories.
Collapse
Affiliation(s)
- Tyler K. Chafin
- Department of Biological SciencesUniversity of ArkansasFayettevilleARUSA
- Present address:
Department of Ecology and Evolutionary BiologyUniversity of ColoradoBoulderCOUSA
| | - Zachery D. Zbinden
- Department of Biological SciencesUniversity of ArkansasFayettevilleARUSA
| | - Marlis R. Douglas
- Department of Biological SciencesUniversity of ArkansasFayettevilleARUSA
| | - Bradley T. Martin
- Department of Biological SciencesUniversity of ArkansasFayettevilleARUSA
| | | | - M. Cory Gray
- Research DivisionArkansas Game and Fish CommissionLittle RockARUSA
| | | | - Michael E. Douglas
- Department of Biological SciencesUniversity of ArkansasFayettevilleARUSA
| |
Collapse
|
28
|
Casanova A, Maroso F, Blanco A, Hermida M, Ríos N, García G, Manuzzi A, Zane L, Verissimo A, García-Marín JL, Bouza C, Vera M, Martínez P. Low impact of different SNP panels from two building-loci pipelines on RAD-Seq population genomic metrics: case study on five diverse aquatic species. BMC Genomics 2021; 22:150. [PMID: 33653268 PMCID: PMC7927381 DOI: 10.1186/s12864-021-07465-w] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2020] [Accepted: 02/22/2021] [Indexed: 02/06/2023] Open
Abstract
Background The irruption of Next-generation sequencing (NGS) and restriction site-associated DNA sequencing (RAD-seq) in the last decade has led to the identification of thousands of molecular markers and their genotyping for refined genomic screening. This approach has been especially useful for non-model organisms with limited genomic resources. Many building-loci pipelines have been developed to obtain robust single nucleotide polymorphism (SNPs) genotyping datasets using a de novo RAD-seq approach, i.e. without reference genomes. Here, the performances of two building-loci pipelines, STACKS 2 and Meyer’s 2b-RAD v2.1 pipeline, were compared using a diverse set of aquatic species representing different genomic and/or population structure scenarios. Two bivalve species (Manila clam and common edible cockle) and three fish species (brown trout, silver catfish and small-spotted catshark) were studied. Four SNP panels were evaluated in each species to test both different building-loci pipelines and criteria for SNP selection. Furthermore, for Manila clam and brown trout, a reference genome approach was used as control. Results Despite different outcomes were observed between pipelines and species with the diverse SNP calling and filtering steps tested, no remarkable differences were found on genetic diversity and differentiation within species with the SNP panels obtained with a de novo approach. The main differences were found in brown trout between the de novo and reference genome approaches. Genotyped vs missing data mismatches were the main genotyping difference detected between the two building-loci pipelines or between the de novo and reference genome comparisons. Conclusions Tested building-loci pipelines for selection of SNP panels seem to have low influence on population genetics inference across the diverse case-study scenarios here studied. However, preliminary trials with different bioinformatic pipelines are suggested to evaluate their influence on population parameters according with the specific goals of each study. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-021-07465-w.
Collapse
Affiliation(s)
- Adrián Casanova
- Department of Zoology, Genetics and Physical Anthropology, ACUIGEN group, Faculty of Veterinary, Universidade de Santiago de Compostela, Campus of Lugo, 27002, Lugo, Spain
| | - Francesco Maroso
- Department of Zoology, Genetics and Physical Anthropology, ACUIGEN group, Faculty of Veterinary, Universidade de Santiago de Compostela, Campus of Lugo, 27002, Lugo, Spain.,Present address: Dipartimento di Scienze della Vita e Biotecnologia (SVeB), Università degli Studi di Ferrara, via Luigi Borsari, 46 - 44121, Ferrara, Italy
| | - Andrés Blanco
- Department of Zoology, Genetics and Physical Anthropology, ACUIGEN group, Faculty of Veterinary, Universidade de Santiago de Compostela, Campus of Lugo, 27002, Lugo, Spain
| | - Miguel Hermida
- Department of Zoology, Genetics and Physical Anthropology, ACUIGEN group, Faculty of Veterinary, Universidade de Santiago de Compostela, Campus of Lugo, 27002, Lugo, Spain
| | - Néstor Ríos
- Sección Genética Evolutiva. Facultad de Ciencias, UdelaR, Iguá 4225, 11400, Montevideo, Uruguay
| | - Graciela García
- Sección Genética Evolutiva. Facultad de Ciencias, UdelaR, Iguá 4225, 11400, Montevideo, Uruguay
| | - Alice Manuzzi
- National Institute of Aquatic Resources, Technical University of Denmark, Vejlsøvej 39, 8600, Silkeborg, Denmark
| | - Lorenzo Zane
- Department of Biology, University of Padova, via U. Bassi 58/B, 35131, Padova, Italy.,Consorzio Nazionale Interuniversitario per le Scienze del Mare (CoNISMa), Piazzale Flaminio 9, 00196, Rome, Italy
| | - Ana Verissimo
- CIBIO - U.P. - Research Center for Biodiversity and Genetic Resources, Campus Agrário de Vairão, 4485-661, Vairão, Portugal.,Virginia Institute of Marine Science, College of William and Mary, Route 1208, Greate Road, Gloucester Point, VA, 23062, USA
| | - José-Luís García-Marín
- Laboratori d'Ictiologia Genètica, Departamento de Biología, Faculty of Sciences, University of Girona, Campus of Montilivi, ES-17071, Girona, Spain
| | - Carmen Bouza
- Department of Zoology, Genetics and Physical Anthropology, ACUIGEN group, Faculty of Veterinary, Universidade de Santiago de Compostela, Campus of Lugo, 27002, Lugo, Spain.,Instituto de Acuicultura, Universidade de Santiago de Compostela, 15705, Santiago de Compostela, Spain
| | - Manuel Vera
- Department of Zoology, Genetics and Physical Anthropology, ACUIGEN group, Faculty of Veterinary, Universidade de Santiago de Compostela, Campus of Lugo, 27002, Lugo, Spain. .,Instituto de Acuicultura, Universidade de Santiago de Compostela, 15705, Santiago de Compostela, Spain.
| | - Paulino Martínez
- Department of Zoology, Genetics and Physical Anthropology, ACUIGEN group, Faculty of Veterinary, Universidade de Santiago de Compostela, Campus of Lugo, 27002, Lugo, Spain.,Instituto de Acuicultura, Universidade de Santiago de Compostela, 15705, Santiago de Compostela, Spain
| |
Collapse
|
29
|
Cerca J, Maurstad MF, Rochette NC, Rivera‐Colón AG, Rayamajhi N, Catchen JM, Struck TH. Removing the bad apples: A simple bioinformatic method to improve loci‐recovery in de novo RADseq data for non‐model organisms. Methods Ecol Evol 2021. [DOI: 10.1111/2041-210x.13562] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- José Cerca
- Frontiers in Evolutionary Zoology Natural History MuseumUniversity of Oslo Oslo Norway
- Department of Environmental Science, Policy, and Management University of California Berkeley CA USA
- Department of Natural History NTNU University MuseumNorwegian University of Science and Technology Trondheim Norway
| | - Marius F. Maurstad
- Frontiers in Evolutionary Zoology Natural History MuseumUniversity of Oslo Oslo Norway
- Centre for Ecological and Evolutionary Synthesis University of Oslo Oslo Norway
| | - Nicolas C. Rochette
- Department of Evolution, Ecology, and Behavior University of Illinois at Urbana‐ChampaignUrbana‐Champaign IL USA
- Department of Ecology and Evolutionary Biology University of California Los Angeles CA USA
| | - Angel G. Rivera‐Colón
- Department of Evolution, Ecology, and Behavior University of Illinois at Urbana‐ChampaignUrbana‐Champaign IL USA
| | - Niraj Rayamajhi
- Department of Evolution, Ecology, and Behavior University of Illinois at Urbana‐ChampaignUrbana‐Champaign IL USA
| | - Julian M. Catchen
- Department of Evolution, Ecology, and Behavior University of Illinois at Urbana‐ChampaignUrbana‐Champaign IL USA
| | - Torsten H. Struck
- Frontiers in Evolutionary Zoology Natural History MuseumUniversity of Oslo Oslo Norway
| |
Collapse
|
30
|
Cerca J, Rivera-Colón AG, Ferreira MS, Ravinet M, Nowak MD, Catchen JM, Struck TH. Incomplete lineage sorting and ancient admixture, and speciation without morphological change in ghost-worm cryptic species. PeerJ 2021; 9:e10896. [PMID: 33614296 PMCID: PMC7879940 DOI: 10.7717/peerj.10896] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2020] [Accepted: 01/13/2021] [Indexed: 12/14/2022] Open
Abstract
Morphologically similar species, that is cryptic species, may be similar or quasi-similar owing to the deceleration of morphological evolution and stasis. While the factors underlying the deceleration of morphological evolution or stasis in cryptic species remain unknown, decades of research in the field of paleontology on punctuated equilibrium have originated clear hypotheses. Species are expected to remain morphologically identical in scenarios of shared genetic variation, such as hybridization and incomplete lineage sorting, or in scenarios where bottlenecks reduce genetic variation and constrain the evolution of morphology. Here, focusing on three morphologically similar Stygocapitella species, we employ a whole-genome amplification method (WGA) coupled with double-digestion restriction-site associated DNA sequencing (ddRAD) to reconstruct the evolutionary history of the species complex. We explore population structure, use population-level statistics to determine the degree of connectivity between populations and species, and determine the most likely demographic scenarios which generally reject for recent hybridization. We find that the combination of WGA and ddRAD allowed us to obtain genomic-level data from microscopic eukaryotes (∼1 millimetre) opening up opportunities for those working with population genomics and phylogenomics in such taxa. The three species share genetic variance, likely from incomplete lineage sorting and ancient admixture. We speculate that the degree of shared variation might underlie morphological similarity in the Atlantic species complex.
Collapse
Affiliation(s)
- José Cerca
- Department of Environmental Science, Policy, and Management, University of California, University of California, Berkeley, Berkeley, CA, United States of America
- Department of Natural History, NTNU University Museum, Norwegian University of Science and Technology, Trondheim, Norway
- Natural History Museum, University of Oslo, Oslo, Norway
| | - Angel G. Rivera-Colón
- Department of Evolution, Ecology, and Behavior, University of Illinois at Urbana-Champaign, Urbana Champaign, IL, United States of America
| | - Mafalda S. Ferreira
- Division of Biological Sciences, University of Montana, Missoula, MT, United States of America
- Departamento de Biologia, Universidade do Porto, Porto, Porto, Portugal
- CIBIO, Centro de Investigação em Biodiversidade e Recursos Genéticos, InBIO Laboratório Associado, Universidade do Porto, Porto, Porto, Portugal
| | - Mark Ravinet
- School of Life Sciences, University of Nottingham, Nottingham, United Kingdom
- Centre for Ecological and Evolutionary Synthesis, University of Oslo, Oslo, Norway
| | | | - Julian M. Catchen
- Department of Evolution, Ecology, and Behavior, University of Illinois at Urbana-Champaign, Urbana Champaign, IL, United States of America
| | | |
Collapse
|
31
|
Undin M, Lockhart PJ, Hills SFK, Castro I. Genetic Rescue and the Plight of Ponui Hybrids. FRONTIERS IN CONSERVATION SCIENCE 2021. [DOI: 10.3389/fcosc.2020.622191] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
Long-term sustainable and resilient populations is a key goal of conservation. How to best achieve this is controversial. There are, for instance, polarized views concerning the fitness and conservation value of hybrid populations founded through multi-origin translocations. A classic example concerns Apteryx (kiwi) in New Zealand. The A. mantelli of Ponui Island constitute a hybrid population where the birds are highly successful in their island habitat. A key dilemma for managers is understanding the reason for this success. Are the hybrid birds of Ponui Island of “no future conservation value” as recently asserted, or do they represent an outstanding example of genetic rescue and an important resource for future translocations? There has been a paradigm shift in scientific thinking concerning hybrids, but the ecological significance of admixed genomes remains difficult to assess. This limits what we can currently predict in conservation science. New understanding from genome science challenges the sufficiency of population genetic models to inform decision making and suggests instead that the contrasting outcomes of hybridization, “outbreeding depression” and “heterosis,” require understanding additional factors that modulate gene and protein expression and how these factors are influenced by the environment. We discuss these findings and the investigations that might help us to better understand the birds of Ponui, inform conservation management of kiwi and provide insight relevant for the future survival of Apteryx.
Collapse
|
32
|
Ortiz D, Pekár S, Bilat J, Alvarez N. Poor performance of DNA barcoding and the impact of RAD loci filtering on the species delimitation of an Iberian ant-eating spider. Mol Phylogenet Evol 2020; 154:106997. [PMID: 33164854 DOI: 10.1016/j.ympev.2020.106997] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2020] [Revised: 10/19/2020] [Accepted: 10/20/2020] [Indexed: 11/28/2022]
Abstract
Genomic data provide unprecedented power for species delimitation. However, current implementations are still time and resource consuming. In addition, bioinformatic processing is contentious and its impact on downstream analyses is insufficiently understood. Here we employ ddRAD sequencing and a thorough sampling for species delimitation in Zodarion styliferum, a widespread Iberian ant-eating spider. We explore the influence of the loci filtering strategy on the downstream phylogenetic analyses, genomic clustering and coalescent species delimitation. We also assess the accuracy of one mitochondrial (COI) and one nuclear (ITS) barcode for fast and inexpensive species delineation in the group. Our genomic data strongly support two morphologically cryptic but ecologically divergent lineages, mainly restricted to the central-eastern and western parts of the Iberian Peninsula, respectively. Larger matrices with more missing data showed increased genomic diversity, supporting that bioinformatic strategies to maximize matrix completion disproportionately exclude loci with the highest mutation rates. Moderate loci filtering gave the best results across analyses: although larger matrices returned concatenated phylogenies with higher support, middle-sized matrices performed better in genetic structure analyses. COI displayed high diversity and a conspicuous barcode gap, revealing 13 mitochondrial lineages. Mitonuclear discordance is consistent with ancestral isolation in multiple groups, probably in glacial refugia, followed by range expansion and secondary contact that produced genomic homogenization. Several apparently (unidirectionally) introgressed specimens further challenge the accuracy of species identification through mitochondrial barcodes in the group. Conversely, ITS failed to separate both lineages of Z. styliferum. This study shows an extreme case of mitonuclear discordance that highlights the limitations of single molecular barcodes for species delimitation, even in presence of distinct barcode gaps, and brings new light on the effects of parameterization on shallow-divergence studies using RAD data.
Collapse
Affiliation(s)
- David Ortiz
- Department of Botany and Zoology, Faculty of Science, Masaryk University, Brno, Czech Republic.
| | - Stano Pekár
- Department of Botany and Zoology, Faculty of Science, Masaryk University, Brno, Czech Republic
| | - Julia Bilat
- Geneva Natural History Museum, Geneva, Switzerland
| | - Nadir Alvarez
- Geneva Natural History Museum, Geneva, Switzerland; Department of Genetics & Evolution, University of Geneva, Geneva, Switzerland
| |
Collapse
|
33
|
Ledent A, Gauthier J, Pereira M, Overson R, Laenen B, Mardulyn P, Gradstein SR, de Haan M, Ballings P, Van der Beeten I, Zartman CE, Vanderpoorten A. What do tropical cryptogams reveal? Strong genetic structure in Amazonian bryophytes. THE NEW PHYTOLOGIST 2020; 228:640-650. [PMID: 32488881 DOI: 10.1111/nph.16720] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/31/2020] [Accepted: 05/25/2020] [Indexed: 06/11/2023]
Abstract
Lowland tropical bryophytes have been perceived as excellent dispersers. In such groups, the inverse isolation hypothesis proposes that spatial genetic structure is erased beyond the limits of short-distance dispersal. Here, we determine the influence of environmental variation and geographic barriers on the spatial genetic structure of a widely dispersed and phylogenetically independent sample of Amazonian bryophytes. Single nucleotide polymorphism data were produced from a restriction site-associated DNA sequencing protocol for 10 species and analyzed through F-statistics and Mantel tests. Neither isolation-by-environment nor the impact of geographic barriers were recovered from the analyses. However, significant isolation-by-distance patterns were observed for 8 out of the 10 investigated species beyond the scale of short-distance dispersal (> 1 km), offering evidence contrary to the inverse isolation hypothesis. Despite a cadre of life-history traits and distributional patterns suggesting that tropical bryophytes are highly vagile, our analyses reveal spatial genetic structures comparable to those documented for angiosperms, whose diaspores are orders of magnitude larger. Dispersal limitation for tropical bryophytes flies in the face of traditional assumptions regarding their dispersal potential, and suggests that the plight of this component of cryptic biodiversity is more dire than previously considered in light of accelerated forest fragmentation in the Amazon.
Collapse
Affiliation(s)
- Alice Ledent
- Institute of Botany, University of Liège, B22 Sart Tilman, Liège, 4000, Belgium
| | - Jérémy Gauthier
- Muséum d'Histoire Naturelle, route de Malagnou 1, Genève, 1208, Switzerland
| | - Martinha Pereira
- Department of Biodiversity, National Institute for Amazonian Research, Petrópolis, CEP 69060-001, Manaus, Amazonas, Brazil
| | - Rick Overson
- Global Institute of Sustainability, Arizona State University, Tempe, AZ, 85281, USA
| | - Benjamin Laenen
- Department of Ecology, Environment and Plant Sciences, Science for Life laboratory, Stockholm University, Stockholm, 106 91, Sweden
| | - Patrick Mardulyn
- Evolutionary Biology and Ecology, Université Libre de Bruxelles, Brussels, 1050, Belgium
| | - S Robbert Gradstein
- Institut de Systématique, Evolution, Biodiversité (ISYEB), Muséum national d'Histoire naturelle, CNRS, EPHE, Sorbonne Université, CP 39, 57 rue Cuvier, Paris, F-75005, France
| | - Myriam de Haan
- Research Department, Meise Botanic Garden, Meise, 1860, Belgium
| | - Petra Ballings
- Research Department, Meise Botanic Garden, Meise, 1860, Belgium
| | | | - Charles E Zartman
- Muséum d'Histoire Naturelle, route de Malagnou 1, Genève, 1208, Switzerland
| | - Alain Vanderpoorten
- Institute of Botany, University of Liège, B22 Sart Tilman, Liège, 4000, Belgium
| |
Collapse
|
34
|
Population genetic structure of the great star coral, Montastraea cavernosa, across the Cuban archipelago with comparisons between microsatellite and SNP markers. Sci Rep 2020; 10:15432. [PMID: 32963271 PMCID: PMC7508986 DOI: 10.1038/s41598-020-72112-5] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2020] [Accepted: 08/17/2020] [Indexed: 11/25/2022] Open
Abstract
Coral reef habitats surrounding Cuba include relatively healthy, well-developed shallow and mesophotic (30–150 m) scleractinian communities at the cross-currents of the Tropical Western Atlantic (TWA). However, Cuba’s coral communities are not immune to the declines observed throughout the TWA, and there is limited information available regarding genetic connectivity, diversity, and structure among these populations. This represents an immense gap in our understanding of coral ecology and population dynamics at both local and regional scales. To address this gap, we evaluated the population genetic structure of the coral Montastraea cavernosa across eight reef sites surrounding Cuba. Colonies were genotyped using nine microsatellite markers and > 9,000 single nucleotide polymorphism (SNP) markers generated using the 2bRAD approach to assess fine-scale genetic structure across these sites. Both the microsatellite and SNP analyses identified patterns of genetic differentiation among sample populations. While the microsatellite analyses did not identify significant genetic structure across the seven shallow M. cavernosa sampling sites, the SNP analyses revealed significant pairwise population differentiation, suggesting that differentiation is greater between eastern and western sites. This study provides insight into methodological differences between microsatellite and SNP markers including potential trade-offs between marker-specific biases, sample size, sequencing costs, and the ability to resolve subtle patterns of population genetic structure. Furthermore, this study suggests that locations in western Cuba may play important roles in this species’ regional metapopulation dynamics and therefore may merit incorporation into developing international management efforts in addition to the local management the sites receive.
Collapse
|
35
|
Camacho-Sanchez M, Velo-Antón G, Hanson JO, Veríssimo A, Martínez-Solano Í, Marques A, Moritz C, Carvalho SB. Comparative assessment of range-wide patterns of genetic diversity and structure with SNPs and microsatellites: A case study with Iberian amphibians. Ecol Evol 2020; 10:10353-10363. [PMID: 33072264 PMCID: PMC7548196 DOI: 10.1002/ece3.6670] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2020] [Accepted: 07/22/2020] [Indexed: 11/11/2022] Open
Abstract
Reduced representation genome sequencing has popularized the application of single nucleotide polymorphisms (SNPs) to address evolutionary and conservation questions in nonmodel organisms. Patterns of genetic structure and diversity based on SNPs often diverge from those obtained with microsatellites to different degrees, but few studies have explicitly compared their performance under similar sampling regimes in a shared analytical framework. We compared range‐wide patterns of genetic structure and diversity in two amphibians endemic to the Iberian Peninsula: Hyla molleri and Pelobates cultripes, based on microsatellite (18 and 14 loci) and SNP (15,412 and 33,140 loci) datasets of comparable sample size and spatial extent. Model‐based clustering analyses with STRUCTURE revealed minor differences in genetic structure between marker types, but inconsistent values of the optimal number of populations (K) inferred. SNPs yielded more repeatable and less admixed ancestries with increasing K compared to microsatellites. Genetic diversity was weakly correlated between marker types, with SNPs providing a better representation of southern refugia and of gradients of genetic diversity congruent with the demographic history of both species. Our results suggest that the larger number of loci in a SNP dataset can provide more reliable inferences of patterns of genetic structure and diversity than a typical microsatellite dataset, at least at the spatial and temporal scales investigated.
Collapse
Affiliation(s)
- Miguel Camacho-Sanchez
- CIBIO/InBIO Centro de Investigação em Biodiversidade e Recursos Genéticos da Universidade do Porto Vairão Portugal
| | - Guillermo Velo-Antón
- CIBIO/InBIO Centro de Investigação em Biodiversidade e Recursos Genéticos da Universidade do Porto Vairão Portugal
| | - Jeffrey O Hanson
- CIBIO/InBIO Centro de Investigação em Biodiversidade e Recursos Genéticos da Universidade do Porto Vairão Portugal
| | - Ana Veríssimo
- CIBIO/InBIO Centro de Investigação em Biodiversidade e Recursos Genéticos da Universidade do Porto Vairão Portugal
| | | | - Adam Marques
- CIBIO/InBIO Centro de Investigação em Biodiversidade e Recursos Genéticos da Universidade do Porto Vairão Portugal
| | - Craig Moritz
- Centre for Biodiversity Analysis and Research School of Biology The Australian National University Canberra ACT Australia
| | - Sílvia B Carvalho
- CIBIO/InBIO Centro de Investigação em Biodiversidade e Recursos Genéticos da Universidade do Porto Vairão Portugal
| |
Collapse
|
36
|
Muñoz-Ramírez CP, Barnes DKA, Cárdenas L, Meredith MP, Morley SA, Roman-Gonzalez A, Sands CJ, Scourse J, Brante A. Gene flow in the Antarctic bivalve Aequiyoldia eightsii (Jay, 1839) suggests a role for the Antarctic Peninsula Coastal Current in larval dispersal. ROYAL SOCIETY OPEN SCIENCE 2020; 7:200603. [PMID: 33047024 PMCID: PMC7540763 DOI: 10.1098/rsos.200603] [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: 04/10/2020] [Accepted: 08/21/2020] [Indexed: 05/12/2023]
Abstract
The Antarctic Circumpolar Current (ACC) dominates the open-ocean circulation of the Southern Ocean, and both isolates and connects the Southern Ocean biodiversity. However, the impact on biological processes of other Southern Ocean currents is less clear. Adjacent to the West Antarctic Peninsula (WAP), the ACC flows offshore in a northeastward direction, whereas the Antarctic Peninsula Coastal Current (APCC) follows a complex circulation pattern along the coast, with topographically influenced deflections depending on the area. Using genomic data, we estimated genetic structure and migration rates between populations of the benthic bivalve Aequiyoldia eightsii from the shallows of southern South America and the WAP to test the role of the ACC and the APCC in its dispersal. We found strong genetic structure across the ACC (between southern South America and Antarctica) and moderate structure between populations of the WAP. Migration rates along the WAP were consistent with the APCC being important for species dispersal. Along with supporting current knowledge about ocean circulation models at the WAP, migration from the tip of the Antarctic Peninsula to the Bellingshausen Sea highlights the complexities of Southern Ocean circulation. This study provides novel biological evidence of a role of the APCC as a driver of species dispersal and highlights the power of genomic data for aiding in the understanding of the influence of complex oceanographic processes in shaping the population structure of marine species.
Collapse
Affiliation(s)
- Carlos P. Muñoz-Ramírez
- Instituto de Entomología, Facultad de Ciencias Básicas, Universidad Metropolitana de Ciencias de la Educación, Santiago, Chile
- Facultad de Ciencias, Universidad Católica de la Santísima Concepción, Concepción, Chile
- Centro de Investigación en Biodiversidad y Ambientes Sustentables (CIBAS), Universidad Católica de la Santísima Concepción, Concepción, Chile
| | - David K. A. Barnes
- British Antarctic Survey, Natural Environment Research Council, Cambridge, UK
| | - Leyla Cárdenas
- Centro FONDAP de Investigación Dinámica de Ecosistemas Marinos de Altas Latitudes (IDEAL), Instituto de Ciencias Ambientales y Evolutivas, Facultad de Ciencias, UniversidadAustral de Chile, Valdivia, Chile
| | - Michael P. Meredith
- British Antarctic Survey, Natural Environment Research Council, Cambridge, UK
| | - Simon A. Morley
- British Antarctic Survey, Natural Environment Research Council, Cambridge, UK
| | | | - Chester J. Sands
- British Antarctic Survey, Natural Environment Research Council, Cambridge, UK
| | - James Scourse
- College of Life and Environmental Sciences, University of Exeter, Penryn, Cornwall TR10 9EZ, UK
| | - Antonio Brante
- Facultad de Ciencias, Universidad Católica de la Santísima Concepción, Concepción, Chile
- Centro de Investigación en Biodiversidad y Ambientes Sustentables (CIBAS), Universidad Católica de la Santísima Concepción, Concepción, Chile
| |
Collapse
|
37
|
Silva PIT, Silva-Junior OB, Resende LV, Sousa VA, Aguiar AV, Grattapaglia D. A 3K Axiom SNP array from a transcriptome-wide SNP resource sheds new light on the genetic diversity and structure of the iconic subtropical conifer tree Araucaria angustifolia (Bert.) Kuntze. PLoS One 2020; 15:e0230404. [PMID: 32866150 PMCID: PMC7458329 DOI: 10.1371/journal.pone.0230404] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2020] [Accepted: 08/05/2020] [Indexed: 12/30/2022] Open
Abstract
High-throughput SNP genotyping has become a precondition to move to higher precision and wider genome coverage genetic analysis of natural and breeding populations of non-model species. We developed a 44,318 annotated SNP catalog for Araucaria angustifolia, a grandiose subtropical conifer tree, one of the only two native Brazilian gymnosperms, critically endangered due to its valuable wood and seeds. Following transcriptome assembly and annotation, SNPs were discovered from RNA-seq and pooled RAD-seq data. From the SNP catalog, an Axiom® SNP array with 3,038 validated SNPs was developed and used to provide a comprehensive look at the genetic diversity and structure of 15 populations across the natural range of the species. RNA-seq was a far superior source of SNPs when compared to RAD-seq in terms of conversion rate to polymorphic markers on the array, likely due to the more efficient complexity reduction of the huge conifer genome. By matching microsatellite and SNP data on the same set of A. angustifolia individuals, we show that SNPs reflect more precisely the actual genome-wide patterns of genetic diversity and structure, challenging previous microsatellite-based assessments. Moreover, SNPs corroborated the known major north-south genetic cline, but allowed a more accurate attribution to regional versus among-population differentiation, indicating the potential to select ancestry-informative markers. The availability of a public, user-friendly 3K SNP array for A. angustifolia and a catalog of 44,318 SNPs predicted to provide ~29,000 informative SNPs across ~20,000 loci across the genome, will allow tackling still unsettled questions on its evolutionary history, toward a more comprehensive picture of the origin, past dynamics and future trend of the species' genetic resources. Additionally, but not less importantly, the SNP array described, unlocks the potential to adopt genomic prediction methods to accelerate the still very timid efforts of systematic tree breeding of A. angustifolia.
Collapse
Affiliation(s)
- Pedro Italo T. Silva
- Plant Genetics Laboratory, EMBRAPA Genetic Resources and Biotechnology, DF, Brasilia, Brazil
- University of Brasília, Cell Biology Department, Campus Universitário, DF, Brasília, Brazil
| | - Orzenil B. Silva-Junior
- Plant Genetics Laboratory, EMBRAPA Genetic Resources and Biotechnology, DF, Brasilia, Brazil
| | - Lucileide V. Resende
- Plant Genetics Laboratory, EMBRAPA Genetic Resources and Biotechnology, DF, Brasilia, Brazil
| | - Valderes A. Sousa
- Empresa Brasileira de Pesquisa Agropecuária–EMBRAPA Florestas, PR, Colombo, Brazil
| | - Ananda V. Aguiar
- Empresa Brasileira de Pesquisa Agropecuária–EMBRAPA Florestas, PR, Colombo, Brazil
| | - Dario Grattapaglia
- Plant Genetics Laboratory, EMBRAPA Genetic Resources and Biotechnology, DF, Brasilia, Brazil
- University of Brasília, Cell Biology Department, Campus Universitário, DF, Brasília, Brazil
- Graduate Program in Genomic Sciences, Universidade Católica de Brasília, Brasília, DF, Brazil
| |
Collapse
|
38
|
Bohling J. Evaluating the effect of reference genome divergence on the analysis of empirical RADseq datasets. Ecol Evol 2020; 10:7585-7601. [PMID: 32760550 PMCID: PMC7391306 DOI: 10.1002/ece3.6483] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2020] [Revised: 05/18/2020] [Accepted: 05/20/2020] [Indexed: 12/29/2022] Open
Abstract
The advent of high-throughput sequencing (HTS) has made genomic-level analyses feasible for nonmodel organisms. A critical step of many HTS pipelines involves aligning reads to a reference genome to identify variants. Despite recent initiatives, only a fraction of species has publically available reference genomes. Therefore, a common practice is to align reads to the genome of an organism related to the target species; however, this could affect read alignment and bias genotyping. In this study, I conducted an experiment using empirical RADseq datasets generated for two species of salmonids (Actinopterygii; Teleostei; Salmonidae) to address these questions. There are currently reference genomes for six salmonids of varying phylogenetic distance. I aligned the RADseq data to all six genomes and identified variants with several different genotypers, which were then fed into population genetic analyses. Increasing phylogenetic distance between target species and reference genome reduced the proportion of reads that successfully aligned and mapping quality. Reference genome also influenced the number of SNPs that were generated and depth at those SNPs, although the affect varied by genotyper. Inferences of population structure were mixed: increasing reference genome divergence reduced estimates of differentiation but similar patterns of population relationships were found across scenarios. These findings reveal how the choice of reference genome can influence the output of bioinformatic pipelines. It also emphasizes the need to identify best practices and guidelines for the burgeoning field of biodiversity genomics.
Collapse
Affiliation(s)
- Justin Bohling
- Abernathy Fish Technology Center US Fish and Wildlife Service Longview WA USA
| |
Collapse
|
39
|
Ackiss AS, Larson WA, Stott W. Genotyping-by-sequencing illuminates high levels of divergence among sympatric forms of coregonines in the Laurentian Great Lakes. Evol Appl 2020; 13:1037-1054. [PMID: 32431751 PMCID: PMC7232772 DOI: 10.1111/eva.12919] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2019] [Revised: 01/05/2020] [Accepted: 01/07/2020] [Indexed: 02/06/2023] Open
Abstract
Effective resource management depends on our ability to partition diversity into biologically meaningful units. Recent evolutionary divergence, however, can often lead to ambiguity in morphological and genetic differentiation, complicating the delineation of valid conservation units. Such is the case with the "coregonine problem," where recent postglacial radiations of coregonines into lacustrine habitats resulted in the evolution of numerous species flocks, often with ambiguous taxonomy. The application of genomics methods is beginning to shed light on this problem and the evolutionary mechanisms underlying divergence in these ecologically and economically important fishes. Here, we used restriction site-associated DNA (RAD) sequencing to examine genetic diversity and differentiation among sympatric forms in the Coregonus artedi complex in the Apostle Islands of Lake Superior, the largest lake in the Laurentian Great Lakes. Using 29,068 SNPs, we were able to clearly distinguish among the three most common forms for the first time, as well as identify putative hybrids and potentially misidentified specimens. Population assignment rates for these forms using our RAD data were 93%-100% with the only mis-assignments arising from putative hybrids, an improvement from 62% to 77% using microsatellites. Estimates of pairwise differentiation (F ST: 0.045-0.056) were large given the detection of hybrids, suggesting that reduced fitness of hybrid individuals may be a potential mechanism for the maintenance of differentiation. We also used a newly built C. artedi linkage map to look for islands of genetic divergence among forms and found widespread differentiation across the genome, a pattern indicative of long-term drift, suggesting that these forms have been reproductively isolated for a substantial amount of time. The results of this study provide valuable information that can be applied to develop well-informed management strategies and stress the importance of re-evaluating conservation units with genomic tools to ensure they accurately reflect species diversity.
Collapse
Affiliation(s)
- Amanda S. Ackiss
- Wisconsin Cooperative Fishery Research UnitCollege of Natural ResourcesUniversity of Wisconsin‐Stevens PointStevens PointWisconsin
| | - Wesley A. Larson
- U.S. Geological SurveyWisconsin Cooperative Fishery Research UnitCollege of Natural ResourcesUniversity of Wisconsin‐Stevens PointStevens PointWisconsin
| | - Wendylee Stott
- U.S. Geological SurveyGreat Lakes Science CenterAnn ArborMichigan
| |
Collapse
|
40
|
D'Aloia CC, Andrés JA, Bogdanowicz SM, McCune AR, Harrison RG, Buston PM. Unraveling hierarchical genetic structure in a marine metapopulation: A comparison of three high-throughput genotyping approaches. Mol Ecol 2020; 29:2189-2203. [PMID: 32147850 DOI: 10.1111/mec.15405] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2019] [Revised: 02/05/2020] [Accepted: 03/03/2020] [Indexed: 01/04/2023]
Abstract
Marine metapopulations often exhibit subtle population structure that can be difficult to detect. Given recent advances in high-throughput sequencing, an emerging question is whether various genotyping approaches, in concert with improved sampling designs, will substantially improve our understanding of genetic structure in the sea. To address this question, we explored hierarchical patterns of structure in the coral reef fish Elacatinus lori using a high-resolution approach with respect to both genetic and geographic sampling. Previously, we identified three putative E. lori populations within Belize using traditional genetic markers and sparse geographic sampling: barrier reef and Turneffe Atoll; Glover's Atoll; and Lighthouse Atoll. Here, we systematically sampled individuals at ~10 km intervals throughout these reefs (1,129 individuals from 35 sites) and sequenced all individuals at three sets of markers: 2,418 SNPs; 89 microsatellites; and 57 nonrepetitive nuclear loci. At broad spatial scales, the markers were consistent with each other and with previous findings. At finer spatial scales, there was new evidence of genetic substructure, but our three marker sets differed slightly in their ability to detect these patterns. Specifically, we found subtle structure between the barrier reef and Turneffe Atoll, with SNPs resolving this pattern most effectively. We also documented isolation by distance within the barrier reef. Sensitivity analyses revealed that the number of loci (and alleles) had a strong effect on the detection of structure for all three marker sets, particularly at small spatial scales. Taken together, these results illustrate empirically that high-throughput genotyping data can elucidate subtle genetic structure at previously-undetected scales in a dispersive marine fish.
Collapse
Affiliation(s)
- Cassidy C D'Aloia
- Department of Biological Sciences, University of New Brunswick, Saint John, NB, Canada
| | - Jose A Andrés
- Department of Ecology and Evolutionary Biology, Cornell University, Ithaca, NY, USA
| | - Steven M Bogdanowicz
- Department of Ecology and Evolutionary Biology, Cornell University, Ithaca, NY, USA
| | - Amy R McCune
- Department of Ecology and Evolutionary Biology, Cornell University, Ithaca, NY, USA
| | - Richard G Harrison
- Department of Ecology and Evolutionary Biology, Cornell University, Ithaca, NY, USA
| | - Peter M Buston
- Department of Biology and Marine Program, Boston University, Boston, MA, USA
| |
Collapse
|
41
|
Sunde J, Yıldırım Y, Tibblin P, Forsman A. Comparing the Performance of Microsatellites and RADseq in Population Genetic Studies: Analysis of Data for Pike ( Esox lucius) and a Synthesis of Previous Studies. Front Genet 2020; 11:218. [PMID: 32231687 PMCID: PMC7082332 DOI: 10.3389/fgene.2020.00218] [Citation(s) in RCA: 53] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2019] [Accepted: 02/24/2020] [Indexed: 01/06/2023] Open
Abstract
Population genetic studies reveal biodiversity patterns and inform about drivers of evolutionary differentiation and adaptation, including gene flow, drift and selection. This can advance our understanding and aid decision making regarding management and conservation efforts. Microsatellites have long been used in population genetic studies. Thanks to the development of newer techniques, sequencing approaches such as restriction site associated DNA sequencing (RADseq) are on their way to replace microsatellites for some applications. However, the performance of these two marker types in population genetics have rarely been systematically compared. We utilized three neutrally and adaptively differentiated populations of anadromous pike (Esox lucius) to assess the relative performance of microsatellites and RADseq with respect to resolution and conclusiveness of estimates of population differentiation and genetic structure. To this end, the same set of individuals (N = 64) were genotyped with both RADseq and microsatellite markers. To assess effects of sample size, the same subset of 10 randomly chosen individuals from each population (N = 30 in total) were also genotyped with both methods. Comparisons of estimated genetic diversity and structure showed that both markers were able to uncover genetic structuring. The full RADseq dataset provided the clearest detection of the finer scaled genetic structuring, and the other three datasets (full and subset microsatellite, and subset RADseq) provided comparable results. A search for outlier loci performed on the full SNP dataset pointed to signs of selection potentially associated with salinity and temperature, exemplifying the utility of RADseq to inform about the importance of different environmental factors. To evaluate whether performance differences between the markers are general or context specific, the results of previous studies that have investigated population structure using both marker types were synthesized. The synthesis revealed that RADseq performed as well as, or better than microsatellites in detecting genetic structuring in the included studies. The differences in the ability to detect population structure, both in the present and the previous studies, are likely explained by the higher number of loci typically utilized in RADseq compared to microsatellite analysis, as increasing the number of markers will (regardless of the marker type) increase power and allow for clearer detection and higher resolution of genetic structure.
Collapse
Affiliation(s)
- Johanna Sunde
- Department of Biology and Environmental Science, Centre for Ecology and Evolution in Microbial Model Systems, EEMiS, Linnaeus University, Kalmar, Sweden
| | | | | | | |
Collapse
|
42
|
Multiplex microsatellite PCR panels for the neotropical red mangrove, Rhizophora mangle: combining efforts towards a cost-effective and modifiable tool to better inform conservation and management. CONSERV GENET RESOUR 2020. [DOI: 10.1007/s12686-020-01138-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
Abstract
AbstractBetter-informed mangrove conservation and management practices are needed as the ecosystem services provided by these intertidal forests continue to be threatened by increasing anthropogenic pressures and climate change. Multiple layers of knowledge are required to achieve this goal, including insights into population genetics of mangrove species. Understanding the importance of population-genetic insights to conservation, multiple research groups have developed microsatellite loci for the widespread, neotropical red mangrove, Rhizophora mangle. However, although a wealth of genetic markers exist, empirical research is limited in the number of these loci employed. Here, we designed two multiplex PCR panels that combine seven novel loci developed for this work and eight previously-developed loci from three research groups to generate 15-locus genotypes, more than twice the average number of loci used in previous research, in only two PCR. We demonstrated utility in R. mangle from four sites across ~ 2500 km near this species’ northern latitudinal limits, and that these multiplex panels were better able to delineate populations than data subsets with numbers of loci comparable to previous research. We focus our discussion on how this tool is a more-informative, efficient (both in terms of time and resources), and easily-modifiable alternative to address many pressing conservation and management issues, such as the generation of baseline genetic data for areas not yet studied, better defining management units, and monitoring genetic effects of restoration projects. We also provide a quick protocol that outlines each step in this procedure to facilitate the use of this tool by others.
Collapse
|
43
|
Shi Y, Zhou Z, Liu B, Kong S, Chen B, Bai H, Li L, Pu F, Xu P. Construction of a High-Density Genetic Linkage Map and QTL Mapping for Growth-Related Traits in Takifugu bimaculatus. MARINE BIOTECHNOLOGY (NEW YORK, N.Y.) 2020; 22:130-144. [PMID: 31900733 DOI: 10.1007/s10126-019-09938-2] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/02/2019] [Accepted: 12/03/2019] [Indexed: 06/10/2023]
Abstract
Takifugu bimaculatus is a euryhaline species, distributed ranging from the southern Yellow Sea to the South China Sea. Their tolerance to a wide range of salinity and temperature, coupled with a desirable firm texture, makes T. bimaculatus a strong candidate for Takifugu aquaculture in subtropics areas. Due to the increasing demand in markets and emerging of the Takifugu aquaculture industry, close attention has been paid to improvement on the T. bimaculatus production. In aquaculture, the great effort has been put into marker-assisted selective breeding, and efficient improvement was realized. However, few genetic resources on T. bimaculatus are provided so far. Aiming at understanding the genetic basis underlying important economic growth traits, facilitating genetic improvement and enriching the genetic resource in T. bimaculatus, we constructed the first genetic linkage map for T. bimaculatus via double digestion restriction-site association DNA sequencing and conducted quantitative traits locus (QTL) mapping for growth-related traits. The map comprised 1976 single nucleotide polymorphism markers distributed on 22 linkage groups (LG), with a total genetic distance of 2039.74 cM. Based on the linkage map, a chromosome-level assembly was constructed whereby we carried out comparative genomics analysis, verifying the high accuracy on contigs ordering of the linkage map. On the other hand, 18 QTLs associated with growth traits were detected on LG6, LG7, LG8, LG10, LG20, and LG21 with phenotypical variance ranging from 15.1 to 56.4%. Candidate genes participating in cartilage development, fat accumulation, and other growth-related regulation activities were identified from these QTLs, including col11a1, foxa2, and thrap3. The linkage map provided a solid foundation for chromosomes assembly and refinement. QTLs reported here unraveled the genomic architecture of some growth traits, which will advance the investigation of aquaculture breeding efforts in T. bimaculatus.
Collapse
Affiliation(s)
- Yue Shi
- State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
| | - Zhixiong Zhou
- State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
| | - Bo Liu
- Fisheries Research Institute of Fujian, Xiamen, China
| | - Shengnan Kong
- State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
| | - Baohua Chen
- State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
| | - Huaqiang Bai
- State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
| | - Leibin Li
- Fisheries Research Institute of Fujian, Xiamen, China
| | - Fei Pu
- State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
| | - Peng Xu
- State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences, Xiamen University, Xiamen, China.
- Fujian Key Laboratory of Genetics and Breeding of Marine Organisms, Xiamen University, Xiamen, China.
- Shenzhen Research Institute of Xiamen University, Shenzhen, China.
| |
Collapse
|
44
|
Graham CF, Boreham DR, Manzon RG, Stott W, Wilson JY, Somers CM. How "simple" methodological decisions affect interpretation of population structure based on reduced representation library DNA sequencing: A case study using the lake whitefish. PLoS One 2020; 15:e0226608. [PMID: 31978053 PMCID: PMC6980518 DOI: 10.1371/journal.pone.0226608] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2019] [Accepted: 12/01/2019] [Indexed: 12/30/2022] Open
Abstract
Reduced representation (RRL) sequencing approaches (e.g., RADSeq, genotyping by sequencing) require decisions about how much to invest in genome coverage and sequencing depth, as well as choices of values for adjustable bioinformatics parameters. To empirically explore the importance of these “simple” methodological decisions, we generated two independent sequencing libraries for the same 142 individual lake whitefish (Coregonus clupeaformis) using a nextRAD RRL approach: (1) a larger number of loci at low sequencing depth based on a 9mer (library A); and (2) fewer loci at higher sequencing depth based on a 10mer (library B). The fish were selected from populations with different levels of expected genetic subdivision. Each library was analyzed using the STACKS pipeline followed by three types of population structure assessment (FST, DAPC and ADMIXTURE) with iterative increases in the stringency of sequencing depth and missing data requirements, as well as more specific a priori population maps. Library B was always able to resolve strong population differentiation in all three types of assessment regardless of the selected parameters, largely due to retention of more loci in analyses. In contrast, library A produced more variable results; increasing the minimum sequencing depth threshold (-m) resulted in a reduced number of retained loci, and therefore lost resolution at high -m values for FST and ADMIXTURE, but not DAPC. When detecting fine population differentiation, the population map influenced the number of loci and missing data, which generated artefacts in all downstream analyses tested. Similarly, when examining fine scale population subdivision, library B was robust to changing parameters but library A lost resolution depending on the parameter set. We used library B to examine actual subdivision in our study populations. All three types of analysis found complete subdivision among populations in Lake Huron, ON and Dore Lake, SK, Canada using 10,640 SNP loci. Weak population subdivision was detected in Lake Huron with fish from sites in the north-west, Search Bay, North Point and Hammond Bay, showing slight differentiation. Overall, we show that apparently simple decisions about library construction and bioinformatics parameters can have important impacts on the interpretation of population subdivision. Although potentially more costly on a per-locus basis, early investment in striking a balance between the number of loci and sequencing effort is well worth the reduced genomic coverage for population genetics studies. More conservative stringency settings on STACKS parameters lead to a final dataset that was more consistent and robust when examining both weak and strong population differentiation. Overall, we recommend that researchers approach “simple” methodological decisions with caution, especially when working on non-model species for the first time.
Collapse
Affiliation(s)
- Carly F. Graham
- Department of Biology, University of Regina, Regina, Saskatchewan, Canada
| | - Douglas R. Boreham
- Medical Sciences, Northern Ontario School of Medicine, Greater Sudbury, Ontario, Canada
| | - Richard G. Manzon
- Department of Biology, University of Regina, Regina, Saskatchewan, Canada
| | - Wendylee Stott
- Department of Fisheries and Wildlife, Michigan State University, East Lansing, Michigan, USA
| | - Joanna Y. Wilson
- Department of Biology, McMaster University, Hamilton, Ontario, Canada
| | | |
Collapse
|
45
|
Zucchi MI, Cordeiro EMG, Wu X, Lamana LM, Brown PJ, Manjunatha S, Viana JPG, Omoto C, Pinheiro JB, Clough SJ. Population Genomics of the Neotropical Brown Stink Bug, Euschistus heros: The Most Important Emerging Insect Pest to Soybean in Brazil. Front Genet 2019; 10:1035. [PMID: 31749834 PMCID: PMC6844245 DOI: 10.3389/fgene.2019.01035] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2018] [Accepted: 09/27/2019] [Indexed: 11/29/2022] Open
Abstract
Recent changes in soybean management like the adoption of transgenic crops and no-till farming, in addition to the expansion of cultivated areas into new virgin frontiers, are some of the hypotheses that can explain the rise of secondary pests, such as the Neotropical brown stink bug, Euschistus heros, in Brazil. To better access the risk of increased pests like E. heros and to determine probabilities for insecticide resistance spreading, it is necessary first to access the levels of the genetic diversity, how the genetic diversity is distributed, and how natural selection is acting upon the natural variation. Using the genotyping by sequencing (GBS) technique, we generated ~60,000 single-nucleotide polymorphisms (SNPs) distributed across the E. heros genome to answer some of those questions. The SNP data was used to investigate the pattern of genetic structure, hybridization and natural selection of this emerging pest. We found that E. heros populations presented similar levels of genetic diversity with slightly higher values at several central locations in Brazil. Our results also showed strong genetic structure separating northern and southern Brazilian regions (FST = 0.22; p-value = 0.000) with a very distinct hybrid zone at the central region. The analyses also suggest the possibility that GABA channels and odorant receptors might play a role in the process of natural selection. At least one marker was associated with soybean and beans crops, but no association between allele frequency and cotton was found. We discuss the implications of these findings in the management of emerging pests in agriculture, particularly in the context of large areas of monoculture such as soybean and cotton.
Collapse
Affiliation(s)
- Maria I Zucchi
- Institute of Biology, University of Campinas, Campinas, Brazil.,Agência Paulista de Tecnologia dos Agronegócios, Pólo Regional Centro-Sul, Piracicaba, Brazil
| | - Erick M G Cordeiro
- Department of Entomology, University of São Paulo-ESALQ/USP, Piracicaba, Brazil
| | - Xing Wu
- Department of Crop Sciences, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - Letícia Marise Lamana
- Department of Fitotecnia e Fitossanidade, Universidade Estadual de Ponta Grossa, Ponta Grossa, Brazil
| | - Patrick J Brown
- Department of Crop Sciences, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - Shilpa Manjunatha
- Department of Crop Sciences, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - João Paulo Gomes Viana
- Department of Crop Sciences, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - Celso Omoto
- Department of Entomology, University of São Paulo-ESALQ/USP, Piracicaba, Brazil
| | - José B Pinheiro
- Department of Genetics, University of São Paulo-ESALQ/USP, Piracicaba, Brazil
| | - Steven J Clough
- Department of Crop Sciences, University of Illinois at Urbana-Champaign, Urbana, IL, United States.,US Department of Agriculture, Agricultural Research Services, Soybean/Maize Germplasm, Pathology, and Genetics Research Unit, Urbana, IL, United States
| |
Collapse
|
46
|
Colicchio J, Monnahan PJ, Wessinger CA, Brown K, Kern JR, Kelly JK. Individualized mating system estimation using genomic data. Mol Ecol Resour 2019; 20:333-347. [DOI: 10.1111/1755-0998.13094] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2019] [Revised: 07/27/2019] [Accepted: 08/19/2019] [Indexed: 11/29/2022]
Affiliation(s)
- Jack Colicchio
- Department of Ecology and Evolutionary Biology University of Kansas Lawrence KS USA
| | - Patrick J. Monnahan
- Department of Ecology and Evolutionary Biology University of Kansas Lawrence KS USA
| | - Carolyn A. Wessinger
- Department of Ecology and Evolutionary Biology University of Kansas Lawrence KS USA
| | - Keely Brown
- Department of Ecology and Evolutionary Biology University of Kansas Lawrence KS USA
| | - James Russell Kern
- Department of Ecology and Evolutionary Biology University of Kansas Lawrence KS USA
| | - John K. Kelly
- Department of Ecology and Evolutionary Biology University of Kansas Lawrence KS USA
| |
Collapse
|
47
|
Zhou Z, Han K, Wu Y, Bai H, Ke Q, Pu F, Wang Y, Xu P. Genome-Wide Association Study of Growth and Body-Shape-Related Traits in Large Yellow Croaker (Larimichthys crocea) Using ddRAD Sequencing. MARINE BIOTECHNOLOGY (NEW YORK, N.Y.) 2019; 21:655-670. [PMID: 31332575 DOI: 10.1007/s10126-019-09910-0] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/23/2019] [Accepted: 06/26/2019] [Indexed: 06/10/2023]
Abstract
Large yellow croaker (Larimichthys crocea) is an economically important marine fish species of China. Due to overfishing and marine pollution, the wild stocks of this croaker have collapsed in the past decades. Meanwhile, the cultured croaker is facing the difficulties of reduced genetic diversity and low growth rate. To explore the molecular markers related to the growth traits of croaker and providing the related SNPs for the marker-assisted selection, we used double-digest restriction-site associated DNA (ddRAD) sequencing to dissect the genetic bases of growth traits in a cultured population and identify the SNPs that associated with important growth traits by GWAS. A total of 220 individuals were genotyped by ddRAD sequencing. After quality control, 27,227 SNPs were identified in 220 samples and used for GWAS analysis. We identified 13 genome-wide significant associated SNPs of growth traits on 8 chromosomes, and the beta P of these SNPs ranged from 0.01 to 0.86. Through the definition of candidate regions and gene annotation, candidate genes related to growth were identified, including important regulators such as fgf18, fgf1, nr3c1, cyp8b1, fabp2, cyp2r1, ppara, and ccm2l. We also identified SNPs and candidate genes that significantly associated with body shape, including bmp7, col1a1, col11a2, and col18a1, which are also economically important traits for large yellow croaker aquaculture. The results provided insights into the genetic basis of growth and body shape in large yellow croaker population and would provide reliable genetic markers for molecular marker-assisted selection in the future. Meanwhile, the result established a basis for our subsequent fine mapping and related gene study.
Collapse
Affiliation(s)
- Zhixiong Zhou
- State Key Laboratory of Large Yellow Croaker Breeding, Ningde Fufa Fisheries Company Limited, Ningde, 352130, China
- State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences, Xiamen University, Xiamen, 361102, China
| | - Kunhuang Han
- State Key Laboratory of Large Yellow Croaker Breeding, Ningde Fufa Fisheries Company Limited, Ningde, 352130, China
- Key Laboratory of Healthy Mariculture for the East China Sea, Fisheries College, Jimei University, Xiamen, 361021, China
| | - Yidi Wu
- State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences, Xiamen University, Xiamen, 361102, China
| | - Huaqiang Bai
- State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences, Xiamen University, Xiamen, 361102, China
| | - Qiaozhen Ke
- State Key Laboratory of Large Yellow Croaker Breeding, Ningde Fufa Fisheries Company Limited, Ningde, 352130, China
- State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences, Xiamen University, Xiamen, 361102, China
| | - Fei Pu
- State Key Laboratory of Large Yellow Croaker Breeding, Ningde Fufa Fisheries Company Limited, Ningde, 352130, China
- State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences, Xiamen University, Xiamen, 361102, China
| | - Yilei Wang
- State Key Laboratory of Large Yellow Croaker Breeding, Ningde Fufa Fisheries Company Limited, Ningde, 352130, China.
- Key Laboratory of Healthy Mariculture for the East China Sea, Fisheries College, Jimei University, Xiamen, 361021, China.
| | - Peng Xu
- State Key Laboratory of Large Yellow Croaker Breeding, Ningde Fufa Fisheries Company Limited, Ningde, 352130, China.
- State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences, Xiamen University, Xiamen, 361102, China.
- Laboratory for Marine Biology and Biotechnology, Qingdao National Laboratory for Marine Science and Technology, Qingdao, 266071, China.
| |
Collapse
|
48
|
Paetzold C, Wood KR, Eaton DAR, Wagner WL, Appelhans MS. Phylogeny of Hawaiian Melicope (Rutaceae): RAD-seq Resolves Species Relationships and Reveals Ancient Introgression. FRONTIERS IN PLANT SCIENCE 2019; 10:1074. [PMID: 31608076 PMCID: PMC6758601 DOI: 10.3389/fpls.2019.01074] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/24/2019] [Accepted: 08/07/2019] [Indexed: 05/11/2023]
Abstract
Hawaiian Melicope are one of the major adaptive radiations of the Hawaiian Islands comprising 54 endemic species. The lineage is monophyletic with an estimated crown age predating the rise of the current high islands. Phylogenetic inference based on Sanger sequencing has not been sufficient to resolve species or deeper level relationships. Here, we apply restriction site-associated DNA sequencing (RAD-seq) to the lineage to infer phylogenetic relationships. We employ Quartet Sampling to assess information content and statistical support, and to quantify discordance as well as partitioned ABBA-BABA tests to uncover evidence of introgression. Our new results drastically improved resolution of relationships within Hawaiian Melicope. The lineage is divided into five fully supported main clades, two of which correspond to morphologically circumscribed infrageneric groups. We provide evidence for both ancestral and current hybridization events. We confirm the necessity for a taxonomic revision of the Melicope section Pelea, as well as a re-evaluation of several species complexes by combining genomic and morphological data.
Collapse
Affiliation(s)
- Claudia Paetzold
- Department of Systematics, Biodiversity and Evolution of Plants (with Herbarium), University of Göttingen, Goettingen, Germany
| | - Kenneth R. Wood
- National Tropical Botanical Garden, Kalaheo, HI, United States
| | - Deren A. R. Eaton
- Department of Ecology, Evolution and Environmental Biology, Columbia University, New York, NY, USA
- Department of Ecology, Evolution, and Environmental Biology, Columbia University, New York, NY, United States
| | - Warren L. Wagner
- Department of Botany, Smithsonian Institution, Washington, DC, United States
| | - Marc S. Appelhans
- Department of Systematics, Biodiversity and Evolution of Plants (with Herbarium), University of Göttingen, Goettingen, Germany
- Department of Botany, Smithsonian Institution, Washington, DC, United States
| |
Collapse
|
49
|
Warschefsky EJ, von Wettberg EJB. Population genomic analysis of mango (Mangifera indica) suggests a complex history of domestication. THE NEW PHYTOLOGIST 2019; 222:2023-2037. [PMID: 30730057 DOI: 10.1111/nph.15731] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/01/2018] [Accepted: 01/21/2019] [Indexed: 05/22/2023]
Abstract
Humans have domesticated diverse species from across the plant kingdom, yet much of our foundational knowledge of domestication has come from studies investigating relatively few of the most important annual food crops. Here, we examine the impacts of domestication on genetic diversity in a tropical perennial fruit species, mango (Mangifera indica). We used restriction site associated DNA sequencing to generate genomic single nucleotide polymorphism (SNP) data from 106 mango cultivars from seven geographical regions along with 52 samples of closely related species and unidentified cultivars to identify centers of mango genetic diversity and examine how post-domestication dispersal shaped the geographical distribution of diversity. We identify two gene pools of cultivated mango, representing Indian and Southeast Asian germplasm. We found no significant genetic bottleneck associated with the introduction of mango into new regions of the world. By contrast, we show that mango populations in introduced regions have elevated levels of diversity. Our results suggest that mango has a more complex history of domestication than previously supposed, perhaps including multiple domestication events, hybridization and regional selection. Our work has direct implications for mango breeding and genebank management, and also builds on recent efforts to understand how woody perennial crops respond to domestication.
Collapse
Affiliation(s)
- Emily J Warschefsky
- Biological Sciences, Florida International University, 11200 SW 8th St., Miami, FL, 33199, USA
| | - Eric J B von Wettberg
- Biological Sciences, Florida International University, 11200 SW 8th St., Miami, FL, 33199, USA
- Plant and Soil Science, The University of Vermont, 63 Carrigan Drive, Burlington, VT, USA
| |
Collapse
|
50
|
Curto M, Winter S, Seiter A, Schmid L, Scheicher K, Barthel LMF, Plass J, Meimberg H. Application of a SSR-GBS marker system on investigation of European Hedgehog species and their hybrid zone dynamics. Ecol Evol 2019; 9:2814-2832. [PMID: 30891219 PMCID: PMC6405497 DOI: 10.1002/ece3.4960] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2018] [Revised: 12/10/2018] [Accepted: 12/31/2018] [Indexed: 01/02/2023] Open
Abstract
By applying second-generation sequencing technologies to microsatellite genotyping, sequence information is produced which can result in high-resolution population genetics analysis populations and increased replicability between runs and laboratories. In the present study, we establish an approach to study the genetic structure patterns of two European hedgehog species Erinaceaus europaeus and E. roumanicus. These species are usually associated with human settlements and are good models to study anthropogenic impacts on the genetic diversity of wild populations. The short sequence repeats genotyping by sequence (SSR-GBS) method presented uses amplicon sequences to determine genotypes for which allelic variants can be defined according to both length and single nucleotide polymorphisms (SNPs). To evaluate whether complete sequence information improved genetic structure definition, we compared this information with datasets based solely on length information. We identified a total of 42 markers which were successfully amplified in both species. Overall, genotyping based on complete sequence information resulted in a higher number of alleles, as well as greater genetic diversity and differentiation between species. Additionally, the structure patterns were slightly clearer with a division between both species and some potential hybrids. There was some degree of genetic structure within species, although only in E. roumanicus was this related to geographical distance. The statistically significant results obtained by SSR-GBS demonstrate that it is superior to electrophoresis-based methods for SSR genotyping. Moreover, the greater reproducibility and throughput with lower effort which can be obtained with SSR-GBS and the possibility to include degraded DNA into the analysis, allow for continued relevance of SSR markers during the genomic era.
Collapse
Affiliation(s)
- Manuel Curto
- Institute for Integrative Nature Conservation ResearchUniversity of Natural Resources and Life Sciences (BOKU)ViennaAustria
| | - Silvia Winter
- Institute for Integrative Nature Conservation ResearchUniversity of Natural Resources and Life Sciences (BOKU)ViennaAustria
- Division of Plant ProtectionUniversity of Natural Resources and Life Sciences (BOKU)ViennaAustria
| | - Anna Seiter
- Institute for Integrative Nature Conservation ResearchUniversity of Natural Resources and Life Sciences (BOKU)ViennaAustria
| | - Lukas Schmid
- Institute for Integrative Nature Conservation ResearchUniversity of Natural Resources and Life Sciences (BOKU)ViennaAustria
| | - Klaus Scheicher
- Institute of MathematicsUniversity of Natural Resources and Life Sciences (BOKU)ViennaAustria
| | - Leon M. F. Barthel
- Evolutionary EcologyLeibniz Institute for Zoo and Wildlife Research (IZW)BerlinGermany
| | - Jürgen Plass
- Biologiezentrum LinzOberösterreich LandesmuseumLinzAustria
| | - Harald Meimberg
- Institute for Integrative Nature Conservation ResearchUniversity of Natural Resources and Life Sciences (BOKU)ViennaAustria
| |
Collapse
|