1
|
Setsuko S, Narita S, Tamaki I, Sugai K, Nagano AJ, Ujino-Ihara T, Kato H, Isagi Y. Adaptive radiation of the Callicarpa genus in the Bonin Islands revealed through double-digest restriction site-associated DNA sequencing analysis. Ecol Evol 2024; 14:e70216. [PMID: 39279792 PMCID: PMC11393766 DOI: 10.1002/ece3.70216] [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: 06/07/2024] [Revised: 08/06/2024] [Accepted: 08/09/2024] [Indexed: 09/18/2024] Open
Abstract
The Bonin Islands, comprised of the Mukojima, Chichijima, and Hahajima Islands, are known for their isolated and distinctive habitats, hosting a diverse array of endemic flora and fauna. In these islands, adaptive radiation has played a remarkable role in speciation, particularly evident in the Callicarpa genus that is represented by three species: Callicarpa parvifolia and C. glabra exclusive to the Chichijima Islands, and Callicarpa subpubescens, distributed across the entire Bonin Islands. Notably, C. subpubescens exhibits multiple ecotypes, differing in leaf hair density, flowering time, and tree size. In this study, we aimed to investigate species and ecotype diversification patterns, estimate divergence times, and explore cryptic species within Callicarpa in the Bonin Islands using phenotypic and genetic data (double-digest restriction site-associated DNA sequencing). Genetic analysis revealed that C. parvifolia and C. glabra both formed single, distinct genetic groups. Conversely, C. subpubescens consisted of six genetic groups corresponding to different ecotypes and regions, and a hybrid group resulting from the hybridization between two of these genetic groups. Population demography analysis focusing on six Chichijima and Hahajima Islands-based species/ecotypes indicated that all species and ecotypes except one ecotype diverged simultaneously around 73-77 kya. The star-shaped neighbor-net tree also suggests the simultaneous divergence of species and ecotypes. The species and ecotypes that simultaneously diverged adapted to dry environments and understory forests, suggesting that aridification may have contributed to this process of adaptive radiation. Moreover, leaf morphology, flowering time, and genetic analyses suggested the presence of two cryptic species and one hybrid species within C. subpubescens.
Collapse
Affiliation(s)
- Suzuki Setsuko
- Department of Forest Molecular Genetics and Biotechnology, Forestry and Forest Products Research Institute Forest Research and Management Organization Tsukuba Ibaraki Japan
| | - Satoshi Narita
- Graduate School of Agriculture Kyoto University Sakyo-ku Kyoto Japan
| | - Ichiro Tamaki
- Faculty of Applied Biological Sciences, Gifu Field Science Center Gifu University Gifu Gifu Japan
| | - Kyoko Sugai
- Institute of Agricultural and Life Sciences Academic Assembly, Shimane University Matsue Shimane Japan
| | - Atsushi J Nagano
- Faculty of Agriculture Ryukoku University Otsu Shiga Japan
- Institute for Advanced Biosciences Keio University Tsuruoka Yamagata Japan
| | - Tokuko Ujino-Ihara
- Department of Forest Molecular Genetics and Biotechnology, Forestry and Forest Products Research Institute Forest Research and Management Organization Tsukuba Ibaraki Japan
| | - Hidetoshi Kato
- Makino Herbarium Tokyo Metropolitan University Hachioji Tokyo Japan
| | - Yuji Isagi
- Graduate School of Agriculture Kyoto University Sakyo-ku Kyoto Japan
| |
Collapse
|
2
|
Meleshko O, Martin M, Flatberg K, Stenøien H, Korneliussen T, Szövényi P, Hassel K. Linked Selection and Gene Density Shape Genome-Wide Patterns of Diversification in Peatmosses. Evol Appl 2024; 17:e13767. [PMID: 39165607 PMCID: PMC11333200 DOI: 10.1111/eva.13767] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2023] [Revised: 07/22/2024] [Accepted: 07/29/2024] [Indexed: 08/22/2024] Open
Abstract
Genome evolution under speciation is poorly understood in nonmodel and nonvascular plants, such as bryophytes-the largest group of nonvascular land plants. Their genomes are structurally different from angiosperms and likely subjected to stronger linked selection pressure, which may have profound consequences on genome evolution in diversifying lineages, even more so when their genome architecture is conserved. We use the highly diverse, rapidly radiated group of peatmosses (Sphagnum) to characterize the processes affecting genome diversification in bryophytes. Using whole-genome sequencing data from populations of 12 species sampled at different phylogenetic and geographical scales, we describe high correlation of the genomic landscapes of differentiation, divergence, and diversity in Sphagnum. Coupled with evidence from the patterns of covariation among different measures of genetic diversity, phylogenetic discordance, and gene density, this provides strong support that peatmoss genome evolution has been shaped by the long-term effects of linked selection, constrained by distribution of selection targets in the genome. Thus, peatmosses join the growing number of animal and plant groups where functional features of the genome, such as gene density, and linked selection drive genome evolution along predetermined and highly similar routes in different species. Our findings demonstrate the great potential of bryophytes for studying the genomics of speciation and highlight the urgent need to expand the genomic resources in this remarkable group of plants.
Collapse
Affiliation(s)
- Olena Meleshko
- Department of Natural History, NTNU University MuseumNorwegian University of Science and TechnologyTrondheimNorway
| | - Michael D. Martin
- Department of Natural History, NTNU University MuseumNorwegian University of Science and TechnologyTrondheimNorway
| | - Kjell Ivar Flatberg
- Department of Natural History, NTNU University MuseumNorwegian University of Science and TechnologyTrondheimNorway
| | - Hans K. Stenøien
- Department of Natural History, NTNU University MuseumNorwegian University of Science and TechnologyTrondheimNorway
| | | | - Péter Szövényi
- Department of Systematic and Evolutionary Botany & Zurich‐Basel Plant Science CenterUniversity of ZurichZurichSwitzerland
| | - Kristian Hassel
- Department of Natural History, NTNU University MuseumNorwegian University of Science and TechnologyTrondheimNorway
| |
Collapse
|
3
|
Lichman V, Ozerov M, López ME, Noreikiene K, Kahar S, Pukk L, Burimski O, Gross R, Vasemägi A. Whole-genome analysis reveals phylogenetic and demographic history of Eurasian perch. JOURNAL OF FISH BIOLOGY 2024. [PMID: 38897597 DOI: 10.1111/jfb.15821] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/25/2023] [Revised: 04/19/2024] [Accepted: 05/14/2024] [Indexed: 06/21/2024]
Abstract
The contemporary diversity and distribution of species are shaped by their evolutionary and ecological history. This can be deciphered with the help of phylogenetic and demographic analysis methods, ideally combining and supplementing information from mitochondrial and nuclear genomes. In this study, we investigated the demographic history of Eurasian perch (Perca fluviatilis), a highly adaptable teleost with a distribution range across Eurasia. We combined whole-genome resequencing data with available genomic resources to analyse the phylogeny, phylogeography, and demographic history of P. fluviatilis populations from Europe and Siberia. We identified five highly diverged evolutionary mtDNA lineages, three of which show a strong signal of admixture in the Baltic Sea region. The estimated mean divergence time between these lineages ranged from 0.24 to 1.42 million years. Based on nuclear genomes, two distinct demographic trajectories were observed in European and Siberian samples reflecting contrasting demographic histories ca. 30,000-100,000 years before the present. A comparison of mtDNA and nuclear DNA evolutionary trees and AMOVA revealed concordances, as well as incongruences, between the two types of data, most likely reflecting recent postglacial colonization and hybridization events. Overall, our findings demonstrate the power and usefulness of genome-wide information for delineating historical processes that have shaped the genome of P. fluviatilis. We also highlight the added value of data-mining existing transcriptomic resources to complement novel sequence data, helping to shed light on putative glacial refugia and postglacial recolonization routes.
Collapse
Affiliation(s)
- Vitalii Lichman
- Institute of Veterinary Medicine and Animal Sciences, Chair of Aquaculture, Estonian University of Life Sciences, Tartu, Estonia
| | - Mikhail Ozerov
- Biodiversity Unit, University of Turku, Turku, Finland
- Department of Aquatic Resources, Institute of Freshwater Research, Swedish University of Agricultural Sciences, Drottningholm, Sweden
- Department of Biology, University of Turku, Turku, Finland
| | - María-Eugenia López
- Department of Aquatic Resources, Institute of Freshwater Research, Swedish University of Agricultural Sciences, Drottningholm, Sweden
| | - Kristina Noreikiene
- Institute of Veterinary Medicine and Animal Sciences, Chair of Aquaculture, Estonian University of Life Sciences, Tartu, Estonia
- Department of Botany and Genetics, Vilnius University, Vilnius, Lithuania
| | - Siim Kahar
- Institute of Veterinary Medicine and Animal Sciences, Chair of Aquaculture, Estonian University of Life Sciences, Tartu, Estonia
| | - Lilian Pukk
- Institute of Veterinary Medicine and Animal Sciences, Chair of Aquaculture, Estonian University of Life Sciences, Tartu, Estonia
| | - Oksana Burimski
- Institute of Veterinary Medicine and Animal Sciences, Chair of Aquaculture, Estonian University of Life Sciences, Tartu, Estonia
| | - Riho Gross
- Institute of Veterinary Medicine and Animal Sciences, Chair of Aquaculture, Estonian University of Life Sciences, Tartu, Estonia
| | - Anti Vasemägi
- Institute of Veterinary Medicine and Animal Sciences, Chair of Aquaculture, Estonian University of Life Sciences, Tartu, Estonia
- Department of Aquatic Resources, Institute of Freshwater Research, Swedish University of Agricultural Sciences, Drottningholm, Sweden
| |
Collapse
|
4
|
Kamimura Y, Nishikawa M, Yamasako J. DNA barcoding of Japanese earwig species (Insecta, Dermaptera), with sequence diversity analyses of three species of Anisolabididae. Biodivers Data J 2023; 11:e107001. [PMID: 38318508 PMCID: PMC10840520 DOI: 10.3897/bdj.11.e107001] [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: 05/26/2023] [Accepted: 09/19/2023] [Indexed: 02/07/2024] Open
Abstract
Dermaptera is a polyneopteran insect order that includes more than 2,000 described species, commonly known as earwigs, that mainly inhabit tropical, subtropical and warm temperate regions. Although 40 species have been found in Japan, their distribution and habitat preferences have remained ambiguous due to sample misidentification, particularly amongst immature specimens. To overcome this problem, we sequenced and analysed the DNA barcoding region of the mitochondrial cytochrome oxidase I gene (cox1) of dermapteran species recorded from Japan. Including publicly available data, 72.5% of known Japanese dermapteran species were subjected to molecular identification. We extensively sampled three wingless species of subfamily Anisolabidinae (Anisolabididae): Anisolabismaritima, Anisolabellamarginalis and Euborelliapallipes. Although these species exhibit similar habitat preferences as semi-synanthropes, A.maritima, a cosmopolitan species with the highest affinity to seashore, had significantly higher sequence diversity than the latter two species, which are considered endemic to East Asia. A similar trend was observed for (at least partly) winged cosmopolitan species of other families. Introgression with the congener Anisolabisseirokui is also suggested for A.maritima. Possible causes of the varying levels of sequence diversity are discussed.
Collapse
Affiliation(s)
| | - Masaru Nishikawa
- Ehime University, Matsuyama, JapanEhime UniversityMatsuyamaJapan
| | - Junsuke Yamasako
- National Agriculture and Food Research Organization, Tsukuba, JapanNational Agriculture and Food Research OrganizationTsukubaJapan
| |
Collapse
|
5
|
Dai M, He SL, Chen B, Li TJ. Phylogeny of Rhynchium and Its Related Genera (Hymenoptera: Eumeninae) Based on Universal Single-Copy Orthologs and Ultraconserved Elements. INSECTS 2023; 14:775. [PMID: 37754743 PMCID: PMC10532281 DOI: 10.3390/insects14090775] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/12/2023] [Revised: 09/15/2023] [Accepted: 09/18/2023] [Indexed: 09/28/2023]
Abstract
The subfamily Eumeninae is a large group of fierce predatory insects that prey mainly on the larvae of Lepidoptera pests. Because of the highly similar morphologies of the genus Rhynchium and its related genera in the subfamily, including Rhynchium Spinola, Allorhynchium van der Vecht, Anterhynchium de Saussure, Pararrhynchium de Saussure, it is essential to delineate their relationships. A previous phylogenetic analysis based on mitochondrial genomes suggested the inconsistent relationships of these genera under traditional classification based on morphological characters. In this study, we first used single-copy orthologs [USCO] and ultraconserved elements [UCE] extracted from 10 newly sequenced low-coverage whole genomes to resolve the phylogenetic relationships of the above genera. The newly sequenced genomes are 152.99 Mb to 211.49 Mb in size with high completeness (BUSCO complete: 91.5-95.6%) and G + C content (36.31-38.76%). Based on extracted 5811 USCOs and 2312 UCEs, the phylogenetic relationships of Rhynchium and its related genera were: ((Allorhynchium + Lissodynerus) + (Pararrhynchium + (Anterhynchium + (Dirhynchium + Rhynchium)))), which was consistent with the mitochondrial genome results. The results supported the genus Rhynchium as monophyletic, whereas Anterhynchium was recovered as paraphyletic, with Anterhynchium (Dirhynchium) as a sister to Rhynchium and hence deserving generic status; In addition, in the genus Pararrhynchium, P. septemfasciatus feanus and P. venkataramani were separated, not clustered on a branch, just as P. septemfasciatus feanus was not together with P. striatum based on mitochondrial genomes. Since Lissodynerus septemfasciatus, the type species of the genus Lissodynerus, was transferred to Pararrhynchium, it is considered that the genus Lissodynerus should be restituted as a valid genus, not a synonym of Pararrhynchium.
Collapse
Affiliation(s)
| | | | | | - Ting-Jing Li
- Chongqing Key Laboratory of Vector Insects, Institute of Entomology and Molecular Biology, College of Life Science, Chongqing Normal University, Chongqing 401331, China; (M.D.); (S.-L.H.); (B.C.)
| |
Collapse
|
6
|
Dong X, Zhang H, Zhu X, Wang K, Xue H, Ye Z, Zheng C, Bu W. Mitochondrial introgression and mito-nuclear discordance obscured the closely related species boundaries in Cletus Stål from China (Heteroptera: Coreidae). Mol Phylogenet Evol 2023; 184:107802. [PMID: 37221926 DOI: 10.1016/j.ympev.2023.107802] [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: 11/08/2022] [Revised: 04/07/2023] [Accepted: 04/18/2023] [Indexed: 05/25/2023]
Abstract
Accurate taxonomy and delimitation are of great importance for pest control strategies and management programs. Here, we focus on Cletus (Insecta: Hemiptera: Coreidae), which includes many crop pests. The species boundaries still conflict and only cytochrome c oxidase subunit I (COI) barcoding has been previously used for molecular studies. We generated new mitochondrial genome and nuclear genome-wide SNPs to explore the species boundaries of 46 Cletus samples from China using multiple species delimitation approaches. All results recovered a monophyly with high support, except for two closely related species in clade I - C. punctiger and C. graminis. Mitochondrial data demonstrated admixture in clade I, while genome-wide SNPs unambiguously identified two separate species, which were confirmed by morphological classification. Inconsistent nuclear and mitochondrial data indicated mito-nuclear discordance. Mitochondrial introgression is the most likely explanation, and more extensive sampling and more comprehensive data are needed to ascertain a pattern. Accurate species delimitation will shed light on species status; thus, an accurate taxonomy is of particular concern, as there is a pressing need to implement precise control of agricultural pests and to perform further research on diversification.
Collapse
Affiliation(s)
- Xue Dong
- Institute of Entomology, College of Life Sciences, Nankai University, Tianjin 300071, China
| | - Haiguang Zhang
- College of Life Science, Linyi University, Linyi 276000, China
| | - Xiuxiu Zhu
- Institute of Entomology, College of Life Sciences, Nankai University, Tianjin 300071, China
| | - Kaibin Wang
- Institute of Entomology, College of Life Sciences, Nankai University, Tianjin 300071, China
| | - Huaijun Xue
- Institute of Entomology, College of Life Sciences, Nankai University, Tianjin 300071, China
| | - Zhen Ye
- Institute of Entomology, College of Life Sciences, Nankai University, Tianjin 300071, China.
| | - Chenguang Zheng
- Institute of Entomology, College of Life Sciences, Nankai University, Tianjin 300071, China.
| | - Wenjun Bu
- Institute of Entomology, College of Life Sciences, Nankai University, Tianjin 300071, China.
| |
Collapse
|
7
|
Hibbins MS, Breithaupt LC, Hahn MW. Phylogenomic comparative methods: Accurate evolutionary inferences in the presence of gene tree discordance. Proc Natl Acad Sci U S A 2023; 120:e2220389120. [PMID: 37216509 PMCID: PMC10235958 DOI: 10.1073/pnas.2220389120] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2022] [Accepted: 04/24/2023] [Indexed: 05/24/2023] Open
Abstract
Phylogenetic comparative methods have long been a mainstay of evolutionary biology, allowing for the study of trait evolution across species while accounting for their common ancestry. These analyses typically assume a single, bifurcating phylogenetic tree describing the shared history among species. However, modern phylogenomic analyses have shown that genomes are often composed of mosaic histories that can disagree both with the species tree and with each other-so-called discordant gene trees. These gene trees describe shared histories that are not captured by the species tree, and therefore that are unaccounted for in classic comparative approaches. The application of standard comparative methods to species histories containing discordance leads to incorrect inferences about the timing, direction, and rate of evolution. Here, we develop two approaches for incorporating gene tree histories into comparative methods: one that constructs an updated phylogenetic variance-covariance matrix from gene trees, and another that applies Felsenstein's pruning algorithm over a set of gene trees to calculate trait histories and likelihoods. Using simulation, we demonstrate that our approaches generate much more accurate estimates of tree-wide rates of trait evolution than standard methods. We apply our methods to two clades of the wild tomato genus Solanum with varying rates of discordance, demonstrating the contribution of gene tree discordance to variation in a set of floral traits. Our approaches have the potential to be applied to a broad range of classic inference problems in phylogenetics, including ancestral state reconstruction and the inference of lineage-specific rate shifts.
Collapse
Affiliation(s)
- Mark S Hibbins
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, ON M5S 3B2, Canada
- Department of Biology, Indiana University, Bloomington, IN 47405
| | - Lara C Breithaupt
- Department of Biology, Indiana University, Bloomington, IN 47405
- Department of Computer Science, Duke University, Durham, NC 27710
| | - Matthew W Hahn
- Department of Biology, Indiana University, Bloomington, IN 47405
- Department of Computer Science, Indiana University, Bloomington, IN 47405
| |
Collapse
|
8
|
DeSalle R, Narechania A, Tessler M. Multiple Outgroups Can Cause Random Rooting in Phylogenomics. Mol Phylogenet Evol 2023; 184:107806. [PMID: 37172862 DOI: 10.1016/j.ympev.2023.107806] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2022] [Revised: 02/06/2023] [Accepted: 04/26/2023] [Indexed: 05/15/2023]
Abstract
Outgroup selection has been a major challenge since the rise of phylogenetics, and it has remained so in the phylogenomic era. Our goal here is to use large phylogenomic animal datasets to examine the impact of outgroup selection on the final topology. The results of our analyses further solidify the fact that distant outgroups can cause random rooting, and that this holds for concatenated and coalescent-based methods. The results also indicate that the standard practice of using multiple outgroups often causes random rooting. Most researchers go out of their way to get multiple outgroups, as this has been standard practice for decades. Based on our findings, this practice should stop. Instead, our results suggest that a single (most closely) related relative should be selected as the outgroup, unless all outgroups are roughly equally closely related to the ingroup.
Collapse
Affiliation(s)
- Rob DeSalle
- Institute for Comparative Genomics, American Museum of Natural History, New York, NY 10024, USA; Division of Invertebrate Zoology, American Museum of Natural History, New York, NY 10024, USA
| | - Apurva Narechania
- Institute for Comparative Genomics, American Museum of Natural History, New York, NY 10024, USA
| | - Michael Tessler
- Institute for Comparative Genomics, American Museum of Natural History, New York, NY 10024, USA; Division of Invertebrate Zoology, American Museum of Natural History, New York, NY 10024, USA; St. Francis College, Department of Biology, Brooklyn, NY 11201, USA
| |
Collapse
|
9
|
Phylogeography of Ramalina farinacea (Lichenized Fungi, Ascomycota) in the Mediterranean Basin, Europe, and Macaronesia. DIVERSITY 2023. [DOI: 10.3390/d15030310] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/25/2023]
Abstract
Ramalina farinacea is an epiphytic lichen-forming fungus with a broad geographic distribution, especially in the Northern Hemisphere. In the eighties of the last century, it was hypothesized that R. farinacea had originated in the Macaronesian–Mediterranean region, with the Canary Islands as its probable southernmost limit, and thereafter it would have increased its distribution area. In order to explore the phylogeography of this emblematic lichen, we analyzed 120 thalli of R. farinacea collected in 38 localities distributed in temperate and boreal Europe, the Western Mediterranean Basin, and several Macaronesian archipelagos in the Atlantic Ocean. Data from two nuclear markers (nrITS and uid70) of the mycobiont were obtained to calculate genetic diversity indices to infer the phylogenies and haplotype networks and to investigate population structure. In addition, dating analysis was conducted to provide a valuable hypothesis of the timing of the origin and diversification of R. farinacea and its close allies. Our results highlight that phylogenetic species circumscription in the “Ramalina farinacea group” is complex and suggests that incomplete lineage sorting is at the base of conflicting phylogenetic signals. The existence of a high number of haplotypes restricted to the Macaronesian region, together with the diversification of R. farinacea in the Pleistocene, suggests that this species and its closest relatives originated during relatively recent geological times and then expanded its range to higher latitudes. However, our data cannot rule out whether the species originated from the Macaronesian archipelagos exclusively or also from the Mediterranean Basin. In conclusion, the present work provides a valuable biogeographical hypothesis for disentangling the evolution of this epiphytic lichen in space and time.
Collapse
|
10
|
Coelho MAG, Pearson GA, Boavida JRH, Paulo D, Aurelle D, Arnaud‐Haond S, Gómez‐Gras D, Bensoussan N, López‐Sendino P, Cerrano C, Kipson S, Bakran‐Petricioli T, Ferretti E, Linares C, Garrabou J, Serrão EA, Ledoux J. Not out of the Mediterranean: Atlantic populations of the gorgonian Paramuricea clavata are a separate sister species under further lineage diversification. Ecol Evol 2023; 13:e9740. [PMID: 36789139 PMCID: PMC9912747 DOI: 10.1002/ece3.9740] [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: 11/28/2022] [Revised: 12/22/2022] [Accepted: 12/27/2022] [Indexed: 01/31/2023] Open
Abstract
The accurate delimitation of species boundaries in nonbilaterian marine taxa is notoriously difficult, with consequences for many studies in ecology and evolution. Anthozoans are a diverse group of key structural organisms worldwide, but the lack of reliable morphological characters and informative genetic markers hampers our ability to understand species diversification. We investigated population differentiation and species limits in Atlantic (Iberian Peninsula) and Mediterranean lineages of the octocoral genus Paramuricea previously identified as P. clavata. We used a diverse set of molecular markers (microsatellites, RNA-seq derived single-copy orthologues [SCO] and mt-mutS [mitochondrial barcode]) at 49 locations. Clear segregation of Atlantic and Mediterranean lineages was found with all markers. Species-tree estimations based on SCO strongly supported these two clades as distinct, recently diverged sister species with incomplete lineage sorting, P. cf. grayi and P. clavata, respectively. Furthermore, a second putative (or ongoing) speciation event was detected in the Atlantic between two P. cf. grayi color morphotypes (yellow and purple) using SCO and supported by microsatellites. While segregating P. cf. grayi lineages showed considerable geographic structure, dominating circalittoral communities in southern (yellow) and western (purple) Portugal, their occurrence in sympatry at some localities suggests a degree of reproductive isolation. Overall, our results show that previous molecular and morphological studies have underestimated species diversity in Paramuricea occurring in the Iberian Peninsula, which has important implications for conservation planning. Finally, our findings validate the usefulness of phylotranscriptomics for resolving evolutionary relationships in octocorals.
Collapse
Affiliation(s)
- Márcio A. G. Coelho
- Centre for Marine Sciences (CCMAR)University of AlgarveFaroPortugal
- MARE – Marine and Environmental Sciences CentreISPA‐Instituto UniversitárioLisboaPortugal
| | | | | | - Diogo Paulo
- Centre for Marine Sciences (CCMAR)University of AlgarveFaroPortugal
| | - Didier Aurelle
- Aix Marseille Univ., Université de Toulon, CNRS, IRD, MIOMarseilleFrance
- Institut de Systématique, Evolution, Biodiversité (ISYEB), Muséum National d'Histoire Naturelle, CNRSSorbonne UniversitéParisFrance
| | - Sophie Arnaud‐Haond
- MARBEC (Marine Biodiversity, Exploitation and Conservation)Univ. Montpellier, IFREMER, CNRS, IRDSète CedexFrance
| | - Daniel Gómez‐Gras
- Hawai‘i Institute of Marine BiologyUniversity of Hawai‘i at MānoaKaneoheHawaiiUSA
- Departament de Biologia Evolutiva, Ecologia i Ciències AmbientalsUniversitat de Barcelona (UB)BarcelonaSpain
- Institut de Recerca de la Biodiversitat (IRBio)Universitat de Barcelona (UB)BarcelonaSpain
| | - Nathaniel Bensoussan
- Aix Marseille Univ., Université de Toulon, CNRS, IRD, MIOMarseilleFrance
- Departament de Biologia MarinaInstitut de Ciències del Mar (CSIC)BarcelonaSpain
| | - Paula López‐Sendino
- Departament de Biologia MarinaInstitut de Ciències del Mar (CSIC)BarcelonaSpain
| | - Carlo Cerrano
- Dipartimento di Scienze della Vita e dell’Ambiente (DiSVA)Università Politecnica delle MarcheAnconaItaly
- Consorzio Nazionale Interuniversitario per le Scienze del Mare (CoNISMa)RomeItaly
- Stazione Zoologica Anton DohrnNaplesItaly
- Fano Marine CenterFanoItaly
| | - Silvija Kipson
- Department of Biology, Faculty of ScienceUniversity of ZagrebZagrebCroatia
- SEAFAN – Marine Research & ConsultancyZagrebCroatia
| | | | - Eliana Ferretti
- Studio Associato GAIA s.n.c.GenoaItaly
- Institute of Marine ScienceThe University of AucklandAucklandNew Zealand
| | - Cristina Linares
- Departament de Biologia Evolutiva, Ecologia i Ciències AmbientalsUniversitat de Barcelona (UB)BarcelonaSpain
- Institut de Recerca de la Biodiversitat (IRBio)Universitat de Barcelona (UB)BarcelonaSpain
| | - Joaquim Garrabou
- Aix Marseille Univ., Université de Toulon, CNRS, IRD, MIOMarseilleFrance
- Departament de Biologia MarinaInstitut de Ciències del Mar (CSIC)BarcelonaSpain
| | - Ester A. Serrão
- Centre for Marine Sciences (CCMAR)University of AlgarveFaroPortugal
- CIBIO/InBIO‐Centro de Investigação em Biodiversidade e Recursos GenéticosVairãoPortugal
| | - Jean‐Baptiste Ledoux
- CIIMAR/CIMAR, Centro Interdisciplinar de Investigação Marinha e AmbientalUniversidade do PortoPortoPortugal
| |
Collapse
|
11
|
Rivera D, Prates I, Caldwell JP, Rodrigues MT, Fujita MK. Testing assertions of widespread introgressive hybridization in a clade of neotropical toads with low mate selectivity (Rhinella granulosa species group). Heredity (Edinb) 2023; 130:14-21. [PMID: 36333595 DOI: 10.1038/s41437-022-00571-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2021] [Revised: 10/20/2022] [Accepted: 10/21/2022] [Indexed: 11/06/2022] Open
Abstract
Discordance between different genomic regions, often identified through multilocus sequencing of selected markers, presents particular difficulties in identifying historical processes which drive species diversity and boundaries. Mechanisms causing discordance, such as incomplete lineage sorting or introgression due to interspecific hybridization, are better identified based on population-level genomic datasets. In the toads of the Rhinella granulosa species group, patterns of mito-nuclear discordance and potential hybridization have been reported by several studies. However, these patterns were proposed based on few loci, such that alternative mechanisms behind gene-tree heterogeneity cannot be ruled out. Using genome-wide ddRADseq loci from a subset of species within this clade, we found only partial concordance between currently recognized species-level taxon boundaries and patterns of genetic structure. While most taxa within the R. granulosa group correspond to clades, genetic clustering analyses sometimes grouped distinct taxonomic units into a single cluster. Moreover, levels of admixture between inferred clusters were limited and restricted to a single taxon pair which is best explained by incomplete lineage sorting as opposed to introgressive hybridization, according to D-statistics results. These findings contradict previous assertions of widespread cryptic diversity and gene flow within the R. granulosa clade. Lastly, our analyses suggest that diversification events within the Rhinella granulosa group mostly dated back to the early Pliocene, being generally younger than species divergences in other closely related clades that present high levels of cross-species gene flow. This finding uniquely contradicts common assertions that this young clade of toads exhibits interspecific hybridization.
Collapse
Affiliation(s)
- Danielle Rivera
- Department of Applied Ecology, North Carolina State University, Raleigh, NC, USA.
- Amphibian and Reptile Diversity Research Center, University of Texas at Arlington, Arlington, TX, USA.
| | - Ivan Prates
- Department of Ecology and Evolutionary Biology and Museum of Zoology, University of Michigan, Ann Arbor, MI, USA
| | - Janalee P Caldwell
- Sam Noble Museum and Department of Biology, University of Oklahoma, Norman, OK, USA
| | - Miguel Trefaut Rodrigues
- Departamento de Zoologia, Instituto de Biociências, Universidade de São Paulo, São Paulo, SP, Brazil
| | - Matthew K Fujita
- Amphibian and Reptile Diversity Research Center, University of Texas at Arlington, Arlington, TX, USA
- Department of Biology and Amphibian and Reptile Diversity Research Center, University of Texas at Arlington, Arlington, TX, USA
| |
Collapse
|
12
|
Klymus KE, Hrabik RA, Thompson NL, Cornman RS. Genome resequencing clarifies phylogeny and reveals patterns of selection in the toxicogenomics model Pimephales promelas. PeerJ 2022; 10:e13954. [PMID: 36042859 PMCID: PMC9420404 DOI: 10.7717/peerj.13954] [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: 04/04/2022] [Accepted: 08/05/2022] [Indexed: 01/19/2023] Open
Abstract
Background The fathead minnow (Pimephales promelas) is a model species for toxicological research. A high-quality genome reference sequence is available, and genomic methods are increasingly used in toxicological studies of the species. However, phylogenetic relationships within the genus remain incompletely known and little population-genomic data are available for fathead minnow despite the potential effects of genetic background on toxicological responses. On the other hand, a wealth of extant samples is stored in museum collections that in principle allow fine-scale analysis of contemporary and historical genetic variation. Methods Here we use short-read shotgun resequencing to investigate sequence variation among and within Pimephales species. At the genus level, our objectives were to resolve phylogenetic relationships and identify genes with signatures of positive diversifying selection. At the species level, our objective was to evaluate the utility of archived-sample resequencing for detecting selective sweeps within fathead minnow, applied to a population introduced to the San Juan River of the southwestern United States sometime prior to 1950. Results We recovered well-supported but discordant phylogenetic topologies for nuclear and mitochondrial sequences that we hypothesize arose from mitochondrial transfer among species. The nuclear tree supported bluntnose minnow (P. notatus) as sister to fathead minnow, with the slim minnow (P. tenellus) and bullhead minnow (P. vigilax) more closely related to each other. Using multiple methods, we identified 11 genes that have diversified under positive selection within the genus. Within the San Juan River population, we identified selective-sweep regions overlapping several sets of related genes, including both genes that encode the giant sarcomere protein titin and the two genes encoding the MTORC1 complex, a key metabolic regulator. We also observed elevated polymorphism and reduced differentation among populations (FST) in genomic regions containing certain immune-gene clusters, similar to what has been reported in other taxa. Collectively, our data clarify evolutionary relationships and selective pressures within the genus and establish museum archives as a fruitful resource for characterizing genomic variation. We anticipate that large-scale resequencing will enable the detection of genetic variants associated with environmental toxicants such as heavy metals, high salinity, estrogens, and agrichemicals, which could be exploited as efficient biomarkers of exposure in natural populations.
Collapse
Affiliation(s)
- Katy E. Klymus
- U.S. Geological Survey, Columbia Ecological Research Center, Columbia, MO, USA
| | | | - Nathan L. Thompson
- U.S. Geological Survey, Columbia Ecological Research Center, Columbia, MO, USA
| | - Robert S. Cornman
- U.S. Geological Survey, Fort Collins Science Center, Fort Collins, CO, USA
| |
Collapse
|
13
|
H Tomasco I, Giorello FM, Boullosa N, Feijoo M, Lanzone C, Lessa EP. The contribution of incomplete lineage sorting and introgression to the evolutionary history of the fast-evolving genus Ctenomys (Rodentia, Ctenomyidae). Mol Phylogenet Evol 2022; 176:107593. [PMID: 35905819 DOI: 10.1016/j.ympev.2022.107593] [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: 03/24/2022] [Revised: 06/28/2022] [Accepted: 07/21/2022] [Indexed: 10/31/2022]
Abstract
Incomplete lineage sorting and introgression have been increasingly recognized as important processes involved in biological differentiation. Both incomplete lineage sorting and introgression result in incongruences between gene trees and species trees, consequently causing difficulties in phylogenetic reconstruction. This is particularly the case for rapid radiations, as short internodal distances and incomplete reproductive isolation increase the likelihood of both ILS and introgression. Estimation of the relative frequency of these processes requires assessments across many genomic regions. We use transcriptomics to test for introgression and estimate the frequency of incomplete lineage sorting in a set of three closely related and geographically adjacent South American tuco-tucos species (Ctenomys), a genus comprising 64 species resulting from recent, rapid radiation. After cleaning and filtering, 5764 orthologous genes strongly support paraphyly of C. pearsoni relative to C. brasiliensis (putatively represented by the population of Villa Serrana). In line with earlier phylogenetic work, the C. pearsoni - C. brasiliensis pair is closely related to C. torquatus, whereas C. rionegrensis is more distantly related to these three nominal species. Classical Patterson's D-statistic shows significant signals of introgression from C. torquatus into C. brasiliensis. However, a 5-taxon test shows no significant results. Incomplete lineage sorting was estimated to have involved about 9% of the loci, suggesting it represents an important process in the incipient diversification of tuco-tucos.
Collapse
Affiliation(s)
- Ivanna H Tomasco
- Departamento de Ecología y Evolución, Facultad de Ciencias, Universidad de la República. Iguá 4225. Montevideo, 11400. Uruguay.
| | - Facundo M Giorello
- Facundo M. Giorello. PDU Espacio de Biología Vegetal del Noreste, Centro Universitario de Tacuarembó (CUT), Universidad de la República, Ruta 5 km 386,200, 45000, Tacuarembó, Uruguay
| | - Nicolás Boullosa
- Departamento de Ecología y Evolución, Facultad de Ciencias, Universidad de la República. Iguá 4225. Montevideo, 11400. Uruguay
| | - Matías Feijoo
- Matías Feijoo. Departamento de Sistemas Agrarios y Paisajes Culturales, Centro Universitario Regional Este (CURE). Universidad de la República. Ruta 8 Km 281, Treinta y Tres, Uruguay
| | - Cecilia Lanzone
- Cecilia Lanzone. Laboratorio de Genética Evolutiva, IBS (CONICET-UNaM), FCEQyN, Félix de Azara 1553, Posadas,3300. Misiones, Argentina
| | - Enrique P Lessa
- Departamento de Ecología y Evolución, Facultad de Ciencias, Universidad de la República. Iguá 4225. Montevideo, 11400. Uruguay
| |
Collapse
|
14
|
The Impact of Fast Radiation on the Phylogeny of Bactrocera Fruit Flies as Revealed by Multiple Evolutionary Models and Mutation Rate-Calibrated Clock. INSECTS 2022; 13:insects13070603. [PMID: 35886779 PMCID: PMC9319077 DOI: 10.3390/insects13070603] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/29/2022] [Revised: 06/22/2022] [Accepted: 06/28/2022] [Indexed: 11/16/2022]
Abstract
Several true fruit flies (Tephritidae) cause major damage to agriculture worldwide. Among them, species of the genus Bactrocera are extensively studied to understand the traits associated with their invasiveness and ecology. Comparative approaches based on a reliable phylogenetic framework are particularly effective, but several nodes of the Bactrocera phylogeny are still controversial, especially concerning the reciprocal affinities of the two major pests B. dorsalis and B. tryoni. Here, we analyzed a newly assembled genomic-scaled dataset using different models of evolution to infer a phylogenomic backbone of ten representative Bactrocera species and two outgroups. We further provide the first genome-scaled inference of their divergence by calibrating the clock using fossil records and the spontaneous mutation rate. The results reveal a closer relationship of B. dorsalis with B. latifrons than to B. tryoni, contrary to what was previously supported by mitochondrial-based phylogenies. By employing coalescent-aware and heterogeneous evolutionary models, we show that this incongruence likely derives from a hitherto undetected systematic error, exacerbated by incomplete lineage sorting and possibly hybridization. This agrees with our clock analysis, which supports a rapid and recent radiation of the clade to which B. dorsalis, B. latifrons and B. tryoni belong. These results provide a new picture of Bactrocera phylogeny that can serve as the basis for future comparative analyses.
Collapse
|
15
|
Plastid phylogenomic analyses of the Selaginella sanguinolenta group (Selaginellaceae) reveal conflict signatures resulting from sequence types, outlier genes, and pervasive RNA editing. Mol Phylogenet Evol 2022; 173:107507. [PMID: 35589053 DOI: 10.1016/j.ympev.2022.107507] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2021] [Revised: 04/11/2022] [Accepted: 04/21/2022] [Indexed: 11/23/2022]
Abstract
Different from the generally conserved plastomes (plastid genomes) of most land plants, the Selaginellaceae plastomes exhibit dynamic structure, high GC content and high substitution rates. Previous plastome analyses identified strong conflict on several clades in Selaginella, however the factors causing the conflictions and the impact on the phylogenetic inference have not been sufficiently investigated. Here, we dissect the distribution of phylogenetic signals and conflicts in Selaginella sanguinolenta group, the plastome of which is DR (direct repeats) structure and with genome-wide RNA editing. We analyzed the data sets including 22 plastomes representing all species of the S. sanguinolenta group, covering the entire geographical distribution from the Himalayas to Siberia and the Russian Far East regions. We recovered four different topologies by applying multispecies coalescent (ASTRAL) and concatenation methods (IQ-TREE and RAxML) on four data sets of PC (protein-coding genes), NC (non-coding sequences), PCN (the concatenated PC and NC), and RC (predicted RNA editing sites "C" were corrected by "T"), respectively. Six monophyletic clades, S. nummularifolia clade, S. rossii clade, S. sajanensis clade, S. sanguinolenta I clade, S. sanguinolenta II clade, and S. sanguinolenta III clade, were consistently resolved and supported by the characteristics of GC content, RNA editing frequency, and gene content. However, the relationships among these clades varied across the four topologies. To explore the underlying causes of the uncertainty, we compared the phylogenetic signals of the four topologies. We identified that the sequence types (coding versus non-coding), outlier genes (genes with extremely high |ΔGLS| values), and C-to-U RNA editing frequency in the protein-coding genes were responsible for the unstable phylogenomic relationship. We further revealed a significant positive correlation between the |ΔGLS| values and the variation coefficient of the RNA editing number. Our results demonstrated that the coalescent method performed better than the concatenation method in overcoming the problems caused by outlier genes and extreme RNA editing events. Our study particularly focused on the importance of exploring the plastid phylogenomic conflicts and suggested conducting concatenated analyses cautiously when adopting organelle genome data.
Collapse
|
16
|
Feng S, Bai M, Rivas-González I, Li C, Liu S, Tong Y, Yang H, Chen G, Xie D, Sears KE, Franco LM, Gaitan-Espitia JD, Nespolo RF, Johnson WE, Yang H, Brandies PA, Hogg CJ, Belov K, Renfree MB, Helgen KM, Boomsma JJ, Schierup MH, Zhang G. Incomplete lineage sorting and phenotypic evolution in marsupials. Cell 2022; 185:1646-1660.e18. [PMID: 35447073 PMCID: PMC9200472 DOI: 10.1016/j.cell.2022.03.034] [Citation(s) in RCA: 31] [Impact Index Per Article: 15.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2020] [Revised: 12/22/2021] [Accepted: 03/21/2022] [Indexed: 12/19/2022]
Abstract
Incomplete lineage sorting (ILS) makes ancestral genetic polymorphisms persist during rapid speciation events, inducing incongruences between gene trees and species trees. ILS has complicated phylogenetic inference in many lineages, including hominids. However, we lack empirical evidence that ILS leads to incongruent phenotypic variation. Here, we performed phylogenomic analyses to show that the South American monito del monte is the sister lineage of all Australian marsupials, although over 31% of its genome is closer to the Diprotodontia than to other Australian groups due to ILS during ancient radiation. Pervasive conflicting phylogenetic signals across the whole genome are consistent with some of the morphological variation among extant marsupials. We detected hundreds of genes that experienced stochastic fixation during ILS, encoding the same amino acids in non-sister species. Using functional experiments, we confirm how ILS may have directly contributed to hemiplasy in morphological traits that were established during rapid marsupial speciation ca. 60 mya.
Collapse
Affiliation(s)
- Shaohong Feng
- BGI-Shenzhen, Shenzhen 518083, China; State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming 650223, China
| | - Ming Bai
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China; School of Agriculture, Ningxia University, Yinchuan 750021, China; College of Plant Protection, Hebei Agricultural University, Baoding 071001, China
| | | | - Cai Li
- School of Life Sciences, Sun Yat-sen University, Guangzhou 510275, China
| | | | - Yijie Tong
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China; State Key Laboratory of North China Crop Improvement and Regulation, Hebei Agricultural University, Baoding, Hebei 071001, China; Hainan Yazhou Bay Seed Lab, Building 1, No. 7 Yiju Road, Yazhou District, Sanya, Hainan 572024, China
| | - Haidong Yang
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China; Guangdong Key Laboratory of Animal Conservation and Resource Utilization, Guangdong Public Laboratory of Wild Animal Conservation and Utilization, Institute of Zoology, Guangdong Academy of Sciences, Guangzhou 510260, China
| | - Guangji Chen
- BGI-Shenzhen, Shenzhen 518083, China; College of Life Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Duo Xie
- Shenzhen Institute of Synthetic Biology, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen 518055, China
| | - Karen E Sears
- Department of Ecology and Evolutionary Biology, University of California at Los Angeles, Los Angeles, CA 90095, USA
| | - Lida M Franco
- Facultad de Ciencias Naturales y Matemáticas, Universidad de Ibagué, Carrera 22 Calle 67, Ibagué, Colombia
| | - Juan Diego Gaitan-Espitia
- The Swire Institute of Marine Science and School of Biological Sciences, The University of Hong Kong, Pokfulam, Hong Kong SAR, China
| | - Roberto F Nespolo
- Instituto de Ciencias Ambientales y Evolutivas, Facultad de Ciencias, Universidad Austral de Chile, Campus Isla Teja, Valdivia 5090000, Chile; Center of Applied Ecology and Sustainability (CAPES), Facultad de Ciencias Biológicas, Universidad Católica de Chile, Santiago 6513677, Chile; Millenium Institute for Integrative Biology (iBio), Santiago, Chile; Millennium Nucleus of Patagonian Limit of Life (LiLi), Valdivia, Chile
| | - Warren E Johnson
- Center for Species Survival, Smithsonian Conservation Biology Institute, National Zoological Park, 1500 Remont Road, Front Royal, VA 22630, USA; The Walter Reed Biosystematics Unit, Museum Support Center MRC-534, Smithsonian Institution, 4210 Silver Hill Rd., Suitland, MD 20746-2863, USA; Walter Reed Army Institute of Research, 503 Robert Grant Avenue, Silver Spring, MD 20910, USA
| | - Huanming Yang
- BGI-Shenzhen, Shenzhen 518083, China; James D. Watson Institute of Genome Sciences, Hangzhou 310058, China
| | - Parice A Brandies
- School of Life and Environmental Sciences, University of Sydney, NSW 2006, Australia
| | - Carolyn J Hogg
- School of Life and Environmental Sciences, University of Sydney, NSW 2006, Australia
| | - Katherine Belov
- School of Life and Environmental Sciences, University of Sydney, NSW 2006, Australia
| | - Marilyn B Renfree
- School of BioSciences, The University of Melbourne, Melbourne, VIC 3010, Australia
| | - Kristofer M Helgen
- Australian Museum Research Institute, Australian Museum, Sydney, NSW 2010, Australia; Australian Research Council Centre of Excellence for Australian Biodiversity and Heritage, University of New South Wales, Sydney, NSW 2052, Australia
| | - Jacobus J Boomsma
- Section for Ecology and Evolution, Department of Biology, Universitetsparken 15, University of Copenhagen, 2100 Copenhagen, Denmark
| | | | - Guojie Zhang
- BGI-Shenzhen, Shenzhen 518083, China; State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming 650223, China; Villum Centre for Biodiversity Genomics, Section for Ecology and Evolution, Department of Biology, Universitetsparken 15, University of Copenhagen, 2100 Copenhagen, Denmark; Center for Excellence in Animal Evolution and Genetics, Chinese Academy of Sciences, Kunming 650223, China.
| |
Collapse
|
17
|
Li F, Rane RV, Luria V, Xiong Z, Chen J, Li Z, Catullo RA, Griffin PC, Schiffer M, Pearce S, Lee SF, McElroy K, Stocker A, Shirriffs J, Cockerell F, Coppin C, Sgrò CM, Karger A, Cain JW, Weber JA, Santpere G, Kirschner MW, Hoffmann AA, Oakeshott JG, Zhang G. Phylogenomic analyses of the genus Drosophila reveals genomic signals of climate adaptation. Mol Ecol Resour 2022; 22:1559-1581. [PMID: 34839580 PMCID: PMC9299920 DOI: 10.1111/1755-0998.13561] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2020] [Accepted: 11/10/2021] [Indexed: 01/13/2023]
Abstract
Many Drosophila species differ widely in their distributions and climate niches, making them excellent subjects for evolutionary genomic studies. Here, we have developed a database of high-quality assemblies for 46 Drosophila species and one closely related Zaprionus. Fifteen of the genomes were newly sequenced, and 20 were improved with additional sequencing. New or improved annotations were generated for all 47 species, assisted by new transcriptomes for 19. Phylogenomic analyses of these data resolved several previously ambiguous relationships, especially in the melanogaster species group. However, it also revealed significant phylogenetic incongruence among genes, mainly in the form of incomplete lineage sorting in the subgenus Sophophora but also including asymmetric introgression in the subgenus Drosophila. Using the phylogeny as a framework and taking into account these incongruences, we then screened the data for genome-wide signals of adaptation to different climatic niches. First, phylostratigraphy revealed relatively high rates of recent novel gene gain in three temperate pseudoobscura and five desert-adapted cactophilic mulleri subgroup species. Second, we found differing ratios of nonsynonymous to synonymous substitutions in several hundred orthologues between climate generalists and specialists, with trends for significantly higher ratios for those in tropical and lower ratios for those in temperate-continental specialists respectively than those in the climate generalists. Finally, resequencing natural populations of 13 species revealed tropics-restricted species generally had smaller population sizes, lower genome diversity and more deleterious mutations than the more widespread species. We conclude that adaptation to different climates in the genus Drosophila has been associated with large-scale and multifaceted genomic changes.
Collapse
Affiliation(s)
- Fang Li
- BGI‐ShenzhenShenzhenChina
- Section for Ecology and EvolutionDepartment of BiologyUniversity of CopenhagenCopenhagenDenmark
| | - Rahul V. Rane
- Commonwealth Scientific and Industrial Research OrganisationActonACTAustralia
- Bio21 InstituteSchool of BioSciencesUniversity of MelbourneParkvilleVic.Australia
| | - Victor Luria
- Department of Systems BiologyHarvard Medical SchoolBostonMassachusettsUSA
| | - Zijun Xiong
- BGI‐ShenzhenShenzhenChina
- State Key Laboratory of Genetic Resources and EvolutionKunming Institute of ZoologyChinese Academy of Sciences (CAS)KunmingYunnanChina
- College of Life SciencesUniversity of Chinese Academy of SciencesBeijingChina
| | | | | | - Renee A. Catullo
- Commonwealth Scientific and Industrial Research OrganisationActonACTAustralia
- Division of Ecology and EvolutionCentre for Biodiversity AnalysisThe Australian National UniversityActonACTAustralia
| | - Philippa C. Griffin
- Bio21 InstituteSchool of BioSciencesUniversity of MelbourneParkvilleVic.Australia
| | - Michele Schiffer
- Bio21 InstituteSchool of BioSciencesUniversity of MelbourneParkvilleVic.Australia
- Daintree Rainforest ObservatoryJames Cook UniversityCape TribulationQldAustralia
| | - Stephen Pearce
- Commonwealth Scientific and Industrial Research OrganisationActonACTAustralia
| | - Siu Fai Lee
- Commonwealth Scientific and Industrial Research OrganisationActonACTAustralia
- Applied BioSciencesMacquarie UniversityNorth RydeNSWAustralia
| | - Kerensa McElroy
- Commonwealth Scientific and Industrial Research OrganisationActonACTAustralia
| | - Ann Stocker
- Bio21 InstituteSchool of BioSciencesUniversity of MelbourneParkvilleVic.Australia
| | - Jennifer Shirriffs
- Bio21 InstituteSchool of BioSciencesUniversity of MelbourneParkvilleVic.Australia
| | - Fiona Cockerell
- School of Biological SciencesMonash UniversityClaytonVic.Australia
| | - Chris Coppin
- Commonwealth Scientific and Industrial Research OrganisationActonACTAustralia
| | - Carla M. Sgrò
- School of Biological SciencesMonash UniversityClaytonVic.Australia
| | - Amir Karger
- IT ‐ Research ComputingHarvard Medical SchoolBostonMassachusettsUSA
| | - John W. Cain
- Department of MathematicsHarvard UniversityCambridgeMassachusettsUSA
| | - Jessica A. Weber
- Department of GeneticsHarvard Medical SchoolBostonMassachusettsUSA
| | - Gabriel Santpere
- Neurogenomics Group, Research Programme on Biomedical Informatics (GRIB)Department of Experimental and Health Sciences (DCEXS)Hospital del Mar Medical Research Institute (IMIM)Universitat Pompeu FabraBarcelonaCataloniaSpain
| | - Marc W. Kirschner
- Department of Systems BiologyHarvard Medical SchoolBostonMassachusettsUSA
| | - Ary A. Hoffmann
- Bio21 InstituteSchool of BioSciencesUniversity of MelbourneParkvilleVic.Australia
| | - John G. Oakeshott
- Commonwealth Scientific and Industrial Research OrganisationActonACTAustralia
- Applied BioSciencesMacquarie UniversityNorth RydeNSWAustralia
| | - Guojie Zhang
- BGI‐ShenzhenShenzhenChina
- Section for Ecology and EvolutionDepartment of BiologyUniversity of CopenhagenCopenhagenDenmark
- State Key Laboratory of Genetic Resources and EvolutionKunming Institute of ZoologyChinese Academy of Sciences (CAS)KunmingYunnanChina
- Center for Excellence in Animal Evolution and GeneticsChinese Academy of SciencesKunmingChina
| |
Collapse
|
18
|
Di-Nizo CB, Suárez-Villota EY, Silva MJJ. Species limits and recent diversification of Cerradomys (Sigmodontinae: Oryzomyini) during the Pleistocene. PeerJ 2022; 10:e13011. [PMID: 35480563 PMCID: PMC9037131 DOI: 10.7717/peerj.13011] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2021] [Accepted: 02/04/2022] [Indexed: 01/11/2023] Open
Abstract
Cerradomys is a genus of the tribe Oryzomyini with eight species currently recognized, and a controversial taxonomy. These species are mainly distributed in the South America dry diagonal, but some species extend into Atlantic Forest, reaching the coastal sandy plains known as Restingas. This study aimed to address species limits and patterns of diversification of Cerradomys species. For this purpose, we performed cytogenetic and molecular analyses (phylogeny, coalescent species delimitation, barcoding, and divergence times estimation) using multiple mitochondrial and nuclear markers on a comprehensive sampling, representing all nominal taxa reported so far. Chromosomal information was a robust marker recognizing eight Cerradomys species. Reciprocal monophyly was recovered for all the species, except for C. subflavus. These results together with coalescent analyses recovered eight species as the most congruent species delimitation scenario for the genus (mean C tax : 0.72). Divergence time estimates revealed that Cerradomys' diversification occurred about 1.32 million years ago (Mya) during the Pleistocene. Although our results conservatively support the eight Cerradomys species described so far, different lines of evidence suggest that C. langguthi and C. subflavus could potentially be species-complexes. We discussed this scenario in the light of multiple evolutionary processes within and between species and populations, since Cerradomys comprises a species group with recent diversification affected by Pleistocene climatic changes and by the complex biogeographic history of South America dry diagonal. This work supports that the diversity of Cerradomys is underestimated and reiterates that interdisciplinary approaches are mandatory to identify small rodent species properly, and to unhide cryptic species.
Collapse
Affiliation(s)
- Camilla B. Di-Nizo
- Laboratório de Ecologia e Evolução, Instituto Butantan, São Paulo, São Paulo, Brazil
- Leibniz Institute for the Analysis of Biodiversity Change, Zoological Research Museum A. Koenig, Bonn, Germany
| | - Elkin Y. Suárez-Villota
- Instituto de Ciencias Naturales, Facultad de Medicina Veterinaria y Agronomía, Universidad de las Américas, Concepción, Chile
| | - Maria José J. Silva
- Laboratório de Ecologia e Evolução, Instituto Butantan, São Paulo, São Paulo, Brazil
| |
Collapse
|
19
|
Sanderson MJ, Búrquez A, Copetti D, McMahon MM, Zeng Y, Wojciechowski MF. Origin and diversification of the saguaro cactus (Carnegiea gigantea): a within-species phylogenomic analysis. Syst Biol 2022; 71:1178-1194. [PMID: 35244183 DOI: 10.1093/sysbio/syac017] [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: 03/25/2020] [Revised: 02/18/2022] [Accepted: 02/25/2022] [Indexed: 11/14/2022] Open
Abstract
Reconstructing accurate historical relationships within a species poses numerous challenges, not least in many plant groups in which gene flow is high enough to extend well beyond species boundaries. Nonetheless, the extent of tree-like history within a species is an empirical question on which it is now possible to bring large amounts of genome sequence to bear. We assess phylogenetic structure across the geographic range of the saguaro cactus, an emblematic member of Cactaceae, a clade known for extensive hybridization and porous species boundaries. Using 200 Gb of whole genome resequencing data from 20 individuals sampled from 10 localities, we assembled two data sets comprising 150,000 biallelic single nucleotide polymorphisms (SNPs) from protein coding sequences. From these we inferred within-species trees and evaluated their significance and robustness using five qualitatively different inference methods. Despite the low sequence diversity, large census population sizes, and presence of wide-ranging pollen and seed dispersal agents, phylogenetic trees were well resolved and highly consistent across both data sets and all methods. We inferred that the most likely root, based on marginal likelihood comparisons, is to the east and south of the region of highest genetic diversity, which lies along the coast of the Gulf of California in Sonora, Mexico. Together with striking decreases in marginal likelihood found to the north, this supports hypotheses that saguaro's current range reflects post-glacial expansion from the refugia in the south of its range. We conclude with observations about practical and theoretical issues raised by phylogenomic data sets within species, in which SNP-based methods must be used rather than gene tree methods that are widely used when sequence divergence is higher. These include computational scalability, inference of gene flow, and proper assessment of statistical support in the presence of linkage effects.
Collapse
Affiliation(s)
- Michael J Sanderson
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, AZ 85721, USA
| | - Alberto Búrquez
- Instituto de Ecología, Unidad Hermosillo, Universidad Nacional Autónoma de México, Hermosillo, Sonora, Mexico
| | - Dario Copetti
- Arizona Genomics Institute, School of Plant Sciences, University of Arizona, Tucson, AZ, 85721 USA
| | | | - Yichao Zeng
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, AZ 85721, USA
| | | |
Collapse
|
20
|
Souza CS, Melo BF, M. T. Mattox G, Oliveira C. Phylogenomic analysis of the Neotropical fish subfamily Characinae using ultraconserved elements (Teleostei: Characidae). Mol Phylogenet Evol 2022; 171:107462. [DOI: 10.1016/j.ympev.2022.107462] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2021] [Revised: 03/10/2022] [Accepted: 03/13/2022] [Indexed: 11/16/2022]
|
21
|
Hagberg L, Celemín E, Irisarri I, Hawlitschek O, Bella JL, Mott T, Pereira RJ. Extensive introgression at late stages of species formation: Insights from grasshopper hybrid zones. Mol Ecol 2022; 31:2384-2399. [PMID: 35191134 DOI: 10.1111/mec.16406] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2021] [Revised: 02/02/2022] [Accepted: 02/14/2022] [Indexed: 11/30/2022]
Abstract
The process of species formation is characterised by the accumulation of multiple reproductive barriers. The evolution of hybrid male sterility, or Haldane's rule, typically characterises later stages of species formation, when reproductive isolation is strongest. Yet, understanding how quickly reproductive barriers evolve and their consequences for maintaining genetic boundaries between emerging species remains a challenging task because it requires studying taxa that hybridise in nature. Here, we address these questions using the meadow grasshopper Pseudochorthippus parallelus, where populations that show multiple reproductive barriers, including hybrid male sterility, hybridise in two natural hybrid zones. Using mitochondrial data, we infer that such populations have diverged some 100,000 years ago, at the beginning of the last glacial cycle in Europe. Nuclear data shows that contractions at multiple glacial refugia, and post-glacial expansions have facilitated genetic differentiation between lineages that today interact in hybrid zones. We find extensive introgression throughout the sampled species range, irrespective of current strength of reproductive isolation. Populations exhibiting hybrid male sterility in two hybrid zones show repeatable patterns of genomic differentiation, consistent with shared genomic constraints affecting ancestral divergence or with the role of those regions in reproductive isolation. Together, our results suggest that reproductive barriers that characterise late stages of species formation can evolve relatively quickly, particularly when associated with strong demographic changes. Moreover, we show that such barriers persist in the face of extensive gene flow, allowing future studies to identify associated genomic regions.
Collapse
Affiliation(s)
- Linda Hagberg
- Division of Evolutionary Biology, Faculty of Biology II, Ludwig-Maximilians-Universität München, Grosshaderner Strasse 2, 82152, Planegg-Martinsried, Germany
| | - Enrique Celemín
- Division of Evolutionary Biology, Faculty of Biology II, Ludwig-Maximilians-Universität München, Grosshaderner Strasse 2, 82152, Planegg-Martinsried, Germany.,Unit of Evolutionary Biology/Systematic Zoology, Institute of Biochemistry and Biology, Universität Potsdam, Karl-Liebknecht-Strasse 24-25, 14476, Potsdam, Germany
| | - Iker Irisarri
- University of Goettingen, Institute for Microbiology and Genetics, Department of Applied Bioinformatics, Goldschmidtstr. 1, 37077, Göttingen, Germany.,Campus Institute Data Science (CIDAS), Göttingen, Germany
| | - Oliver Hawlitschek
- Leibniz Institute for the Analysis of Biodiversity Change, Zoological Museum, Martin-Luther-King-Platz 3, 20146, Hamburg, Germany.,Zoologische Staatssammlung (SNSB-ZSM), Münchhausenstr. 21, 81247, Munich, Germany
| | - José L Bella
- Departamento de Biología (Genética), Facultad de Ciencias, Universidad Autónoma de Madrid, 28049, Madrid, Spain.,Centro de Investigación en Biodiversidad y Cambio Global (CIBC-UAM), Universidad Autónoma de Madrid, 28049, Madrid, Spain
| | - Tamí Mott
- Instituto de Ciências Biológicas e da Saúde, Universidade Federal de Alagoas, 57072-900, Maceió, Alagoas, Brazil
| | - Ricardo J Pereira
- Division of Evolutionary Biology, Faculty of Biology II, Ludwig-Maximilians-Universität München, Grosshaderner Strasse 2, 82152, Planegg-Martinsried, Germany
| |
Collapse
|
22
|
Yan Z, Smith ML, Du P, Hahn MW, Nakhleh L. Species Tree Inference Methods Intended to Deal with Incomplete Lineage Sorting Are Robust to the Presence of Paralogs. Syst Biol 2022; 71:367-381. [PMID: 34245291 PMCID: PMC8978208 DOI: 10.1093/sysbio/syab056] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2020] [Revised: 06/23/2021] [Accepted: 06/30/2021] [Indexed: 11/24/2022] Open
Abstract
Many recent phylogenetic methods have focused on accurately inferring species trees when there is gene tree discordance due to incomplete lineage sorting (ILS). For almost all of these methods, and for phylogenetic methods in general, the data for each locus are assumed to consist of orthologous, single-copy sequences. Loci that are present in more than a single copy in any of the studied genomes are excluded from the data. These steps greatly reduce the number of loci available for analysis. The question we seek to answer in this study is: what happens if one runs such species tree inference methods on data where paralogy is present, in addition to or without ILS being present? Through simulation studies and analyses of two large biological data sets, we show that running such methods on data with paralogs can still provide accurate results. We use multiple different methods, some of which are based directly on the multispecies coalescent model, and some of which have been proven to be statistically consistent under it. We also treat the paralogous loci in multiple ways: from explicitly denoting them as paralogs, to randomly selecting one copy per species. In all cases, the inferred species trees are as accurate as equivalent analyses using single-copy orthologs. Our results have significant implications for the use of ILS-aware phylogenomic analyses, demonstrating that they do not have to be restricted to single-copy loci. This will greatly increase the amount of data that can be used for phylogenetic inference.[Gene duplication and loss; incomplete lineage sorting; multispecies coalescent; orthology; paralogy.].
Collapse
Affiliation(s)
- Zhi Yan
- Department of Computer Science, Rice University,
6100 Main Street, Houston, TX 77005, USA
| | - Megan L Smith
- Department of Biology and Department of Computer Science,
Indiana University, 1001 East Third Street, Bloomington,
IN 47405, USA
| | - Peng Du
- Department of Computer Science, Rice University,
6100 Main Street, Houston, TX 77005, USA
| | - Matthew W Hahn
- Department of Biology and Department of Computer Science,
Indiana University, 1001 East Third Street, Bloomington,
IN 47405, USA
| | - Luay Nakhleh
- Department of Computer Science, Rice University,
6100 Main Street, Houston, TX 77005, USA
- Department of BioSciences, Rice University, 6100
Main Street, Houston, TX 77005, USA
| |
Collapse
|
23
|
Hibbins MS, Hahn MW. Phylogenomic approaches to detecting and characterizing introgression. Genetics 2022; 220:iyab173. [PMID: 34788444 PMCID: PMC9208645 DOI: 10.1093/genetics/iyab173] [Citation(s) in RCA: 51] [Impact Index Per Article: 25.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2021] [Accepted: 10/02/2021] [Indexed: 12/26/2022] Open
Abstract
Phylogenomics has revealed the remarkable frequency with which introgression occurs across the tree of life. These discoveries have been enabled by the rapid growth of methods designed to detect and characterize introgression from whole-genome sequencing data. A large class of phylogenomic methods makes use of data across species to infer and characterize introgression based on expectations from the multispecies coalescent. These methods range from simple tests, such as the D-statistic, to model-based approaches for inferring phylogenetic networks. Here, we provide a detailed overview of the various signals that different modes of introgression are expected leave in the genome, and how current methods are designed to detect them. We discuss the strengths and pitfalls of these approaches and identify areas for future development, highlighting the different signals of introgression, and the power of each method to detect them. We conclude with a discussion of current challenges in inferring introgression and how they could potentially be addressed.
Collapse
Affiliation(s)
- Mark S Hibbins
- Department of Biology, Indiana University, Bloomington, IN 47405, USA
| | - Matthew W Hahn
- Department of Biology, Indiana University, Bloomington, IN 47405, USA
- Department of Computer Science, Indiana University, Bloomington, IN 47405, USA
| |
Collapse
|
24
|
Borges R, Boussau B, Szöllősi GJ, Kosiol C. Nucleotide Usage Biases Distort Inferences of the Species Tree. Genome Biol Evol 2022; 14:6496956. [PMID: 34983052 PMCID: PMC8829901 DOI: 10.1093/gbe/evab290] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/27/2021] [Indexed: 12/15/2022] Open
Abstract
Despite the importance of natural selection in species’ evolutionary history, phylogenetic methods that take into account population-level processes typically ignore selection. The assumption of neutrality is often based on the idea that selection occurs at a minority of loci in the genome and is unlikely to compromise phylogenetic inferences significantly. However, genome-wide processes like GC-bias and some variation segregating at the coding regions are known to evolve in the nearly neutral range. As we are now using genome-wide data to estimate species trees, it is natural to ask whether weak but pervasive selection is likely to blur species tree inferences. We developed a polymorphism-aware phylogenetic model tailored for measuring signatures of nucleotide usage biases to test the impact of selection in the species tree. Our analyses indicate that although the inferred relationships among species are not significantly compromised, the genetic distances are systematically underestimated in a node-height-dependent manner: that is, the deeper nodes tend to be more underestimated than the shallow ones. Such biases have implications for molecular dating. We dated the evolutionary history of 30 worldwide fruit fly populations, and we found signatures of GC-bias considerably affecting the estimated divergence times (up to 23%) in the neutral model. Our findings call for the need to account for selection when quantifying divergence or dating species evolution.
Collapse
Affiliation(s)
- Rui Borges
- Institut für Populationsgenetik, Vetmeduni Vienna, Wien, Austria
| | - Bastien Boussau
- Université de Lyon, Université Claude Bernard Lyon 1, CNRS UMR 5558, LBBE, Villeurbanne, France
| | - Gergely J Szöllősi
- Department of Biological Physics, Eötvös University, Budapest , Hungary.,MTA-ELTE "Lendület" Evolutionary Genomics Research Group, Budapest, Hungary.,Evolutionary Systems Research Group, Centre for Ecological Research, Hungarian Academy of Sciences, Tihany, Hungary
| | - Carolin Kosiol
- Institut für Populationsgenetik, Vetmeduni Vienna, Wien, Austria.,Centre for Biological Diversity, University of St Andrews, St Andrews, United Kingdom
| |
Collapse
|
25
|
Hibbins MS, Hahn MW. The effects of introgression across thousands of quantitative traits revealed by gene expression in wild tomatoes. PLoS Genet 2021; 17:e1009892. [PMID: 34748547 PMCID: PMC8601620 DOI: 10.1371/journal.pgen.1009892] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2021] [Revised: 11/18/2021] [Accepted: 10/18/2021] [Indexed: 01/13/2023] Open
Abstract
It is now understood that introgression can serve as powerful evolutionary force, providing genetic variation that can shape the course of trait evolution. Introgression also induces a shared evolutionary history that is not captured by the species phylogeny, potentially complicating evolutionary analyses that use a species tree. Such analyses are often carried out on gene expression data across species, where the measurement of thousands of trait values allows for powerful inferences while controlling for shared phylogeny. Here, we present a Brownian motion model for quantitative trait evolution under the multispecies network coalescent framework, demonstrating that introgression can generate apparently convergent patterns of evolution when averaged across thousands of quantitative traits. We test our theoretical predictions using whole-transcriptome expression data from ovules in the wild tomato genus Solanum. Examining two sub-clades that both have evidence for post-speciation introgression, but that differ substantially in its magnitude, we find patterns of evolution that are consistent with histories of introgression in both the sign and magnitude of ovule gene expression. Additionally, in the sub-clade with a higher rate of introgression, we observe a correlation between local gene tree topology and expression similarity, implicating a role for introgressed cis-regulatory variation in generating these broad-scale patterns. Our results reveal a general role for introgression in shaping patterns of variation across many thousands of quantitative traits, and provide a framework for testing for these effects using simple model-informed predictions. It is now known from studying large genetic datasets that species often hybridize and cross with each other over many generations – a phenomenon known as introgression. Introgression introduces new genetic variation into a population, and this variation can cause traits to be shared among the introgressing species. When researchers study the evolution of trait variation among species, this source of trait sharing is rarely accounted for. Here, we present a statistical model of the effects of introgression on trait variation. This model predicts that, when averaged across many thousands of traits, introgressing species are consistently more similar than expected from standard approaches. Researchers studying gene expression often consider the expression of many thousands of genes, making this a case where the expected effects of introgression are likely to manifest. We tested our model prediction using ovule gene expression data from the wild tomato genus Solanum, in two groups of species with evidence of historical introgression. We found that patterns of expression similarity in both groups are consistent with their histories of introgression and the predictions from our model. Our results highlight the importance of accounting for introgression as a source of trait variation among species.
Collapse
Affiliation(s)
- Mark S. Hibbins
- Department of Biology, Indiana University, Bloomington, Indiana, United States of America
- * E-mail:
| | - Matthew W. Hahn
- Department of Biology, Indiana University, Bloomington, Indiana, United States of America
- Department of Computer Science, Indiana University, Bloomington, Indiana, United States of America
| |
Collapse
|
26
|
Yang W, Feiner N, Salvi D, Laakkonen H, Jablonski D, Pinho C, Carretero MA, Sacchi R, Zuffi MAL, Scali S, Plavos K, Pafilis P, Poulakakis N, Lymberakis P, Jandzik D, Schulte U, Aubret F, Badiane A, Perez I de Lanuza G, Abalos J, While GM, Uller T. Population genomics of wall lizards reflects the dynamic history of the Mediterranean Basin. Mol Biol Evol 2021; 39:6413643. [PMID: 34718699 PMCID: PMC8760935 DOI: 10.1093/molbev/msab311] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The Mediterranean Basin has experienced extensive change in geology and climate over the past six million years. Yet, the relative importance of key geological events for the distribution and genetic structure of the Mediterranean fauna remains poorly understood. Here, we use population genomic and phylogenomic analyses to establish the evolutionary history and genetic structure of common wall lizards (Podarcis muralis). This species is particularly informative because, in contrast to other Mediterranean lizards, it is widespread across the Iberian, Italian, and Balkan Peninsulas, and in extra-Mediterranean regions. We found strong support for six major lineages within P. muralis, which were largely discordant with the phylogenetic relationship of mitochondrial DNA. The most recent common ancestor of extant P. muralis was likely distributed in the Italian Peninsula, and experienced an “Out-of-Italy” expansion following the Messinian salinity crisis (∼5 Mya), resulting in the differentiation into the extant lineages on the Iberian, Italian, and Balkan Peninsulas. Introgression analysis revealed that both inter- and intraspecific gene flows have been pervasive throughout the evolutionary history of P. muralis. For example, the Southern Italy lineage has a hybrid origin, formed through admixture between the Central Italy lineage and an ancient lineage that was the sister to all other P. muralis. More recent genetic differentiation is associated with the onset of the Quaternary glaciations, which influenced population dynamics and genetic diversity of contemporary lineages. These results demonstrate the pervasive role of Mediterranean geology and climate for the evolutionary history and population genetic structure of extant species.
Collapse
Affiliation(s)
- Weizhao Yang
- Department of Biology, Lund University, Lund, 223 62, Sweden
| | - Nathalie Feiner
- Department of Biology, Lund University, Lund, 223 62, Sweden
| | - Daniele Salvi
- Department of Health, Life and Environmental Sciences, University of L'Aquila, 67100, Italy Coppito L'Aquila
| | - Hanna Laakkonen
- Department of Biology, Lund University, Lund, 223 62, Sweden
| | - Daniel Jablonski
- Department of Zoology, Comenius University, Ilkovičova 6, 842 15, Slovakia in Bratislava, Bratislava
| | - Catarina Pinho
- CIBIO/InBIO Research Centre in Biodiversity and Genetic Resources, Universidade do Porto, Campus de Vairão, 4485-661, Vairão, Portugal
| | - Miguel A Carretero
- CIBIO/InBIO Research Centre in Biodiversity and Genetic Resources, Universidade do Porto, Campus de Vairão, 4485-661, Vairão, Portugal.,Departamento de Biologia, Faculdade de Ciências da Universidade do Porto, R. Campo Alegre, s/n, Porto, 4169 - 007, Portugal
| | - Roberto Sacchi
- Department of Earth and Environmental Sciences, University of Pavia, Pavia, 27100, Italy
| | - Marco A L Zuffi
- Museum Natural History, University of Pisa, Pisa, 56011, Italy
| | - Stefano Scali
- Museum of Natural History of Milan, Milano, 20121, Italy
| | | | - Panayiotis Pafilis
- National & Kapodistrian University of Athens, School of Science, Faculty of Biology, Panepistimiopolis 15701, Athens, Greece
| | - Nikos Poulakakis
- Natural History Museum of Crete, School of Sciences and Engineering, University of Crete, Knossos Avenue, Heraklion, 71409, Greece.,Biology Department, School of Sciences and Engineering, University of Crete, Voutes University Campus, Heraklion, 70013, Greece.,Institute of Molecular Biology and Biotechnology (IMBB), Heraklion, 70013, GreeceFoundation for Research and Technology-Hellas (FORTH)
| | - Petros Lymberakis
- Natural History Museum of Crete, School of Sciences and Engineering, University of Crete, Knossos Avenue, Heraklion, 71409, Greece
| | - David Jandzik
- Department of Zoology, Comenius University, Ilkovičova 6, 842 15, Slovakia in Bratislava, Bratislava
| | - Ulrich Schulte
- Büro für Faunistische Gutachten-Dr. Ulrich Schulte, Kaiserstraße 2, Borgholzhausen, 33829, Germany
| | - Fabien Aubret
- Station d'Ecologie Théorique et Expérimentale, CNRS, 09200, Moulis, France.,School of Molecular and Life Sciences, Curtin University, WA, 6102, Australia
| | - Arnaud Badiane
- IMBE, Aix-Marseille Université, Avignon Université, CNRS, IRD, Marseille, France
| | - Guillem Perez I de Lanuza
- Institut Cavanilles de Biodiversitat i Biologia Evolutiva, Universitat de València, Valencia, APT. 22085, 46071, Spain
| | - Javier Abalos
- Institut Cavanilles de Biodiversitat i Biologia Evolutiva, Universitat de València, Valencia, APT. 22085, 46071, Spain
| | - Geoffrey M While
- School of Natural Sciences, University of Tasmania, Sandy Bay, Tasmania, 7005, Australia
| | - Tobias Uller
- Department of Biology, Lund University, Lund, 223 62, Sweden
| |
Collapse
|
27
|
Lin X, Jiang K, Liu W, Liu W, Bu W, Wang X, Mo L. Toward a global DNA barcode reference library of the intolerant nonbiting midge genus Rheocricotopus Brundin, 1956. Ecol Evol 2021; 11:12161-12172. [PMID: 34522368 PMCID: PMC8427567 DOI: 10.1002/ece3.7979] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2021] [Revised: 07/16/2021] [Accepted: 07/19/2021] [Indexed: 01/01/2023] Open
Abstract
Environmental DNA metabarcoding is becoming a predominant tool in biodiversity assessment, as this time- and cost-efficient tactics have the ability to increase monitoring accuracy. As a worldwide distributed genus, Rheocricotopus Brundin, 1956 still does not possess a complete and comprehensive global DNA barcode reference library for biodiversity monitoring. In the present study, we compiled a cytochrome c oxidase subunit 1 (COI) DNA barcode library of Rheocricotopus with 434 barcodes around the world, including 121 newly generated DNA barcodes of 32 morphospecies and 313 public barcodes. Automatic Barcode Gap Discovery (ABGD) was applied on the 434 COI barcodes to provide a comparison between the operational taxonomic units (OTU) number calculated from the Barcode Index Number (BIN) with the "Barcode Gap Analysis" and neighbor-joining (NJ) tree analysis. Consequently, these 434 COI barcodes were clustered into 78 BINs, including 42 new BINs. ABGD yielded 51 OTUs with a prior intraspecific divergence of Pmax = 7.17%, while NJ tree revealed 52 well-separated clades. Conservatively, 14 unknown species and one potential synonym were uncovered with reference to COI DNA barcodes. Besides, based on our ecological analysis, we discovered that annual mean temperature and annual precipitation could be considered as key factors associated with distribution of certain members from this genus. Our global DNA barcode reference library of Rheocricotopus provides one fundamental database for accurate species delimitation in Chironomidae taxonomy and facilitates the biodiversity monitoring of aquatic biota.
Collapse
Affiliation(s)
- Xiao‐Long Lin
- College of Life SciencesNankai UniversityTianjinChina
| | - Kun Jiang
- College of Life SciencesNankai UniversityTianjinChina
| | - Wen‐Bin Liu
- Tianjin Key Laboratory of Conservation and Utilization of Animal DiversityTianjin Normal UniversityTianjinChina
| | - Wei Liu
- College of Life SciencesNankai UniversityTianjinChina
| | - Wen‐Jun Bu
- College of Life SciencesNankai UniversityTianjinChina
| | - Xin‐Hua Wang
- College of Life SciencesNankai UniversityTianjinChina
| | - Lidong Mo
- Crowther LabInstitute of Integrative BiologyETH Zurich (Swiss Federal Institute of Technology)ZurichSwitzerland
| |
Collapse
|
28
|
Sun Q, Yang Y, Hao X, Xiao J, Liu J, Yuan X. Comparative Mitogenomic Analysis of Five Awl Skippers (Lepidoptera: Hesperiidae: Coeliadinae) and Their Phylogenetic Implications. INSECTS 2021; 12:insects12080757. [PMID: 34442323 PMCID: PMC8397065 DOI: 10.3390/insects12080757] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/07/2021] [Revised: 08/11/2021] [Accepted: 08/21/2021] [Indexed: 11/16/2022]
Abstract
Simple Summary The subfamily Coeliadinae (Lepidoptera: Hesperiidae) is a unique group of over 70 species in the butterfly family, and its mitochondrial genome data still needs to be supplemented. This study sequenced and analyzed five additional complete mitochondrial genomes of the Coeliadinae species (Hasora schoenherr, Burara miracula, B. oedipodea, B. harisa, and Badamia exclamationis) and compared them in detail with those of the other known skipper mitogenomes. All five of these mitogenomes have the typical lepidopteran mitogenome characteristics of 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and a non-coding region. Our results indicate that their structure, nucleotide composition, codon usage, secondary structure of tRNAs, and so on, are highly conserved. Expanded sampling and gene data from the GenBank, phylogenetic analyses using maximum likelihood, and Bayesian inference methods indicate that Coeliadinae is monophyletic. These results contribute toward refining the phylogeny. Abstract To determine the significance of mitochondrial genome characteristics in revealing phylogenetic relationships and to shed light on the molecular evolution of the Coeliadinae species, the complete mitochondrial genomes (mitogenomes) of five Coeliadinae species were newly sequenced and analyzed, including Hasora schoenherr, Burara miracula, B. oedipodea, B. harisa, and Badamia exclamationis. The results show that all five mitogenomes are double-strand circular DNA molecules, with lengths of 15,340 bp, 15,295 bp, 15,304 bp, 15,295 bp, and 15,289 bp, respectively, and contain the typical 37 genes and a control region. Most protein-coding genes (PCGs) begin with ATN, with 3 types of stop codons including TAA, TAG, and an incomplete codon T-; most of the genes terminate with TAA. All of the transfer RNA genes (tRNAs) present the typical cloverleaf secondary structure except for the trnS1. Several conserved structural elements are found in the AT-rich region. Phylogenetic analyses based on three datasets (PCGs, PRT, and 12PRT) and using maximum likelihood (ML) and Bayesian inference (BI) methods show strong support for the monophyly of Coeliadinae, and the relationships of the five species are (B. exclamationis + ((B. harisa + (B. oedipodea + B. miracula)) + H. schoenherr)).
Collapse
Affiliation(s)
- Qi Sun
- College of Life Sciences, Northwest A&F University, Yangling 712100, China; (Q.S.); (X.H.)
| | - Yumeng Yang
- College of Natural Resources and Environment, Northwest A&F University, Yangling 712100, China;
| | - Xiangyu Hao
- College of Life Sciences, Northwest A&F University, Yangling 712100, China; (Q.S.); (X.H.)
| | - Jintian Xiao
- Key Laboratory of Plant Protection Resources and Pest Management, Ministry of Education, College of Plant Protection, Northwest A&F University, Yangling 712100, China; (J.X.); (J.L.)
| | - Jiaqi Liu
- Key Laboratory of Plant Protection Resources and Pest Management, Ministry of Education, College of Plant Protection, Northwest A&F University, Yangling 712100, China; (J.X.); (J.L.)
| | - Xiangqun Yuan
- Key Laboratory of Plant Protection Resources and Pest Management, Ministry of Education, College of Plant Protection, Northwest A&F University, Yangling 712100, China; (J.X.); (J.L.)
- Correspondence: ; Tel.: +86-137-5998-5152
| |
Collapse
|
29
|
Genomic Variation Shaped by Environmental and Geographical Factors in Prairie Cordgrass Natural Populations Collected across Its Native Range in the USA. Genes (Basel) 2021; 12:genes12081240. [PMID: 34440416 PMCID: PMC8391649 DOI: 10.3390/genes12081240] [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: 07/06/2021] [Revised: 08/10/2021] [Accepted: 08/11/2021] [Indexed: 11/17/2022] Open
Abstract
Prairie cordgrass (Spartina pectinata Link) is a native perennial warm-season (C4) grass common in North American prairies. With its high biomass yield and abiotic stress tolerance, there is a high potential of developing prairie cordgrass for conservation practices and as a dedicated bioenergy crop for sustainable cellulosic biofuel production. However, as with many other undomesticated grass species, little information is known about the genetic diversity or population structure of prairie cordgrass natural populations as compared to their ecotypic and geographic adaptation in North America. In this study, we sampled and characterized a total of 96 prairie cordgrass natural populations with 9315 high quality SNPs from a genotyping-by-sequencing (GBS) approach. The natural populations were collected from putative remnant prairie sites throughout the Midwest and Eastern USA, which are the major habitats for prairie cordgrass. Partitioning of genetic variance using SNP marker data revealed significant variance among and within populations. Two potential gene pools were identified as being associated with ploidy levels, geographical separation, and climatic separation. Geographical factors such as longitude and altitude, and environmental factors such as annual temperature, annual precipitation, temperature of the warmest month, precipitation of the wettest month, precipitation of Spring, and precipitation of the wettest month are important in affecting the intraspecific distribution of prairie cordgrass. The divergence of prairie cordgrass natural populations also provides opportunities to increase breeding value of prairie cordgrass as a bioenergy and conservation crop.
Collapse
|
30
|
Zhao T, Zwaenepoel A, Xue JY, Kao SM, Li Z, Schranz ME, Van de Peer Y. Whole-genome microsynteny-based phylogeny of angiosperms. Nat Commun 2021; 12:3498. [PMID: 34108452 PMCID: PMC8190143 DOI: 10.1038/s41467-021-23665-0] [Citation(s) in RCA: 41] [Impact Index Per Article: 13.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2020] [Accepted: 05/10/2021] [Indexed: 02/05/2023] Open
Abstract
Plant genomes vary greatly in size, organization, and architecture. Such structural differences may be highly relevant for inference of genome evolution dynamics and phylogeny. Indeed, microsynteny-the conservation of local gene content and order-is recognized as a valuable source of phylogenetic information, but its use for the inference of large phylogenies has been limited. Here, by combining synteny network analysis, matrix representation, and maximum likelihood phylogenetic inference, we provide a way to reconstruct phylogenies based on microsynteny information. Both simulations and use of empirical data sets show our method to be accurate, consistent, and widely applicable. As an example, we focus on the analysis of a large-scale whole-genome data set for angiosperms, including more than 120 available high-quality genomes, representing more than 50 different plant families and 30 orders. Our 'microsynteny-based' tree is largely congruent with phylogenies proposed based on more traditional sequence alignment-based methods and current phylogenetic classifications but differs for some long-contested and controversial relationships. For instance, our synteny-based tree finds Vitales as early diverging eudicots, Saxifragales within superasterids, and magnoliids as sister to monocots. We discuss how synteny-based phylogenetic inference can complement traditional methods and could provide additional insights into some long-standing controversial phylogenetic relationships.
Collapse
Affiliation(s)
- Tao Zhao
- State Key Laboratory of Crop Stress Biology for Arid Areas/Shaanxi Key Laboratory of Apple, College of Horticulture, Northwest A&F University, Yangling, China.
- Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium.
- Center for Plant Systems Biology, VIB, Ghent, Belgium.
| | - Arthur Zwaenepoel
- Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium
- Center for Plant Systems Biology, VIB, Ghent, Belgium
| | - Jia-Yu Xue
- College of Horticulture, Academy for Advanced Interdisciplinary Studies, Nanjing Agricultural University, Nanjing, China
- Institute of Botany, Jiangsu Province and Chinese Academy of Sciences, Nanjing, China
| | - Shu-Min Kao
- Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium
- Center for Plant Systems Biology, VIB, Ghent, Belgium
| | - Zhen Li
- Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium
- Center for Plant Systems Biology, VIB, Ghent, Belgium
| | - M Eric Schranz
- Biosystematics Group, Wageningen University and Research, Wageningen, The Netherlands
| | - Yves Van de Peer
- Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium.
- Center for Plant Systems Biology, VIB, Ghent, Belgium.
- College of Horticulture, Academy for Advanced Interdisciplinary Studies, Nanjing Agricultural University, Nanjing, China.
- Center for Microbial Ecology and Genomics, Department of Biochemistry, Genetics and Microbiology, University of Pretoria, Pretoria, South Africa.
| |
Collapse
|
31
|
Hao X, Liu J, Chiba H, Xiao J, Yuan X. Complete mitochondrial genomes of three skippers in the tribe Aeromachini (Lepidoptera: Hesperiidae: Hesperiinae) and their phylogenetic implications. Ecol Evol 2021; 11:8381-8393. [PMID: 34188893 PMCID: PMC8216930 DOI: 10.1002/ece3.7666] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2020] [Revised: 04/19/2021] [Accepted: 04/28/2021] [Indexed: 11/09/2022] Open
Abstract
The mitochondrial genome is now widely used in the study of phylogenetics and molecular evolution due to its maternal inheritance, fast evolutionary rate, and highly conserved gene content. To explore the phylogenetic relationships of the tribe Aeromachini within the subfamily Hesperiinae at the mitochondrial genomic level, we sequenced and annotated the complete mitogenomes of 3 skippers: Ampittia virgata, Halpe nephele, and Onryza maga (new mitogenomes for 2 genera) with a total length of 15,333 bp, 15,291 bp, and 15,381 bp, respectively. The mitogenomes all contain 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (rRNAs), and a noncoding A + T-rich region and are consistent with other lepidopterans in gene order and type. In addition, we reconstructed the phylogenetic trees of Hesperiinae using maximum likelihood (ML) and Bayesian inference (BI) methods based on mitogenomic data. Results show that the tribe Aeromachini in this study robustly constitute a monophyletic group in the subfamily Hesperiinae, with the relationships Coeliadinae + (Euschemoninae + (Pyrginae + ((Eudaminae + Tagiadinae) + (Heteropterinae + ((Trapezitinae + Barcinae) + Hesperiinae))))). Moreover, our study supports the view that Apostictopterus fuliginosus and Barca bicolor should be placed out of the subfamily Hesperiinae.
Collapse
Affiliation(s)
- Xiangyu Hao
- College of Life SciencesNorthwest A&F UniversityYanglingChina
| | - Jiaqi Liu
- Key Laboratory of Plant Protection Resources and Pest ManagementMinistry of EducationEntomological MuseumCollege of Plant ProtectionNorthwest A&F UniversityYanglingChina
| | | | - Jintian Xiao
- Key Laboratory of Plant Protection Resources and Pest ManagementMinistry of EducationEntomological MuseumCollege of Plant ProtectionNorthwest A&F UniversityYanglingChina
| | - Xiangqun Yuan
- Key Laboratory of Plant Protection Resources and Pest ManagementMinistry of EducationEntomological MuseumCollege of Plant ProtectionNorthwest A&F UniversityYanglingChina
| |
Collapse
|
32
|
Fegies AC, Carmignotto AP, Perez MF, Guilardi MD, Lessinger AC. Molecular Phylogeny of Cryptonanus (Didelphidae: Thylamyini): Evidence for a recent and complex diversification in South American open biomes. Mol Phylogenet Evol 2021; 162:107213. [PMID: 34029717 DOI: 10.1016/j.ympev.2021.107213] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2020] [Revised: 05/07/2021] [Accepted: 05/19/2021] [Indexed: 11/30/2022]
Abstract
Systematic revisions of South American marsupials have contributed to our knowledge about genus and species diversity in the last decades, including studies of the most recently described genus Cryptonanus (Didelphidae), currently comprising four recognized species. Herein we provide the first phylogeny for these mouse opossums based on comprehensive sampling, including representatives from all nominal taxa, encompassing most of the geographic distribution of the genus while also extending its known range. The taxonomic status of Cryptonanus species was explored by analyses of multiple mitochondrial and nuclear DNA markers to assess phylogenetic relationships and to provide divergence time estimates, species delimitations and biogeographical hypotheses. Cryptonanus monophyly remained highly supported despite the inclusion of abundant new data from more than a hundred specimens, comprising 10 independent evolutionary lineages. Species-complexes within valid nominal taxa reveal higher species richness in the genus. Based on divergence estimates from a dated phylogeny, we suggest that Cryptonanus diversified along the Quaternary, with speciation events occurring well into the Pleistocene. The best supported biogeographical hypothesis endorses speciation by vicariance and subset speciation across open formations in shaping the evolutionary history of this didelphid genus, strongly associated with dry tropical landscapes of South America.
Collapse
Affiliation(s)
- Ana Cláudia Fegies
- Departamento de Engenharia Ambiental, Universidade Estadual Paulista (UNESP), Campus Sorocaba, Av. Três de Março 511, Sorocaba, São Paulo CEP 18087-180, Brazil
| | - Ana Paula Carmignotto
- Departamento de Biologia, Universidade Federal de São Carlos (UFSCar), Campus Sorocaba, Rodovia João Leme dos Santos km 110, Sorocaba, São Paulo CEP 18052-780, Brazil.
| | - Manolo Fernandez Perez
- Departamento de Genética e Evolução, Universidade Federal de São Carlos (UFSCar), Rodovia Washington Luís km 235, São Carlos, São Paulo CEP 13565-905, Brazil
| | - Mariana Dias Guilardi
- Laboratório de Ecologia e Evolução, Instituto Butantan, Avenida Vital Brazil 1500, São Paulo, São Paulo CEP 05503-000, Brazil
| | - Ana Cláudia Lessinger
- Departamento de Biologia, Universidade Federal de São Carlos (UFSCar), Campus Sorocaba, Rodovia João Leme dos Santos km 110, Sorocaba, São Paulo CEP 18052-780, Brazil.
| |
Collapse
|
33
|
Dong X, Liang Q, Li J, Feng P. Positive selection drives the evolution of a primate bitter taste receptor gene. Ecol Evol 2021; 11:5459-5467. [PMID: 34026020 PMCID: PMC8131804 DOI: 10.1002/ece3.7440] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2021] [Revised: 02/20/2021] [Accepted: 02/24/2021] [Indexed: 12/22/2022] Open
Abstract
Bitter taste perception is important in food choice of animals, and it is mediated by bitter taste receptor (T2R) containing three regions: extracellular regions (ECs), transmembrane regions (TMs), and intracellular regions (ICs). It is hypothesized that ECs, TMs, and ICs are under different selective pressures, with ECs being unstable while TMs and ICs being constrained. To test this hypothesis, we examined the selective pressures on one of the bitter taste receptor genes-T2R1 and its different areas from 35 primates and found that T2R1 was under neutral evolution but with some positively selected sites in it. Further analyses suggested that TMs, ICs, and the concatenated transmembrane region TM1237 were under purifying selection; in contrast, extracellular regions, the first and second extracellular loop (EL1, EL2), were subject to positive selection with several positively selected sites in them. Therefore, this study supported the above-mentioned hypothesis. The reason why EL1 and EL2 of T2R1 have positively selected sites is probably due to their participation in forming the cap-like structure involved in ligand binding. Positive selection acts as a driving force of the T2R1 functional differentiation and confers the ability to discern various bitter substances for primates, which could help them to adapt to the changing environment during the evolutionary course.
Collapse
Affiliation(s)
- Xiaoyan Dong
- Key Laboratory of Ecology of Rare and Endangered Species and Environmental Protection (Guangxi Normal University)Ministry of EducationGuilinChina
- Guangxi Key Laboratory of Rare and Endangered Animal EcologyGuangxi Normal UniversityGuilinChina
| | - Qiufang Liang
- Key Laboratory of Ecology of Rare and Endangered Species and Environmental Protection (Guangxi Normal University)Ministry of EducationGuilinChina
- Guangxi Key Laboratory of Rare and Endangered Animal EcologyGuangxi Normal UniversityGuilinChina
| | - Jiaping Li
- Key Laboratory of Ecology of Rare and Endangered Species and Environmental Protection (Guangxi Normal University)Ministry of EducationGuilinChina
- Guangxi Key Laboratory of Rare and Endangered Animal EcologyGuangxi Normal UniversityGuilinChina
| | - Ping Feng
- Key Laboratory of Ecology of Rare and Endangered Species and Environmental Protection (Guangxi Normal University)Ministry of EducationGuilinChina
- Guangxi Key Laboratory of Rare and Endangered Animal EcologyGuangxi Normal UniversityGuilinChina
| |
Collapse
|
34
|
Zhao F, Wu YW, Drew BT, Yao G, Chen YP, Cai J, Liu ED, Li B, Xiang CL. Systematic Placement of the Enigmatic Southeast Asian Genus Paralamium and an Updated Phylogeny of Tribe Pogostemoneae (Lamiaceae Subfamily Lamioideae). FRONTIERS IN PLANT SCIENCE 2021; 12:646133. [PMID: 33936133 PMCID: PMC8085563 DOI: 10.3389/fpls.2021.646133] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/25/2020] [Accepted: 03/22/2021] [Indexed: 06/12/2023]
Abstract
Paralamium (Lamiaceae) is a monotypic genus within the subfamily Lamioideae and has a sporadic distribution in subtropical mountains of southeast Asia. Although recent studies have greatly improved our understanding of generic relationships within Lamioideae, the second most species-rich subfamily of Lamiaceae, the systematic position of Paralamium within the subfamily remains unclear. In this study, we investigate the phylogenetic placement of the genus using three datasets: (1) a 69,276 bp plastome alignment of Lamiaceae; (2) a five chloroplast DNA region dataset of tribe Pogostemoneae, and (3) a nuclear ribosomal internal transcribed spacer region dataset of Pogostemoneae. These analyses demonstrate that Paralamium is a member of Pogostemoneae and sister to the monotypic genus Craniotome. In addition, generic-level phylogenetic relationships within Pogostemoneae are also discussed, and a dichotomous key for genera within Pogostemoneae is provided.
Collapse
Affiliation(s)
- Fei Zhao
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, China
| | - Yi-Wen Wu
- College of Life Sciences, Shaanxi Normal University, Xi’an, China
| | - Bryan T. Drew
- Department of Biology, University of Nebraska at Kearney, Kearney, NE, United States
| | - Gang Yao
- South China Limestone Plants Center, College of Forestry and Landscape Architecture, South China Agricultural University, Guangzhou, China
| | - Ya-Ping Chen
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, China
| | - Jie Cai
- Germplasm Bank of Wild Species in Southwest China, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, China
| | - En-De Liu
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, China
| | - Bo Li
- Research Centre of Ecological Sciences, College of Agronomy, Jiangxi Agricultural University, Nanchang, China
| | - Chun-Lei Xiang
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, China
| |
Collapse
|
35
|
Cserhati M. A tail of two pandas- whole genome k-mer signature analysis of the red panda (Ailurus fulgens) and the Giant panda (Ailuropoda melanoleuca). BMC Genomics 2021; 22:228. [PMID: 33794768 PMCID: PMC8015091 DOI: 10.1186/s12864-021-07531-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: 05/18/2020] [Accepted: 03/14/2021] [Indexed: 11/29/2022] Open
Abstract
Background The red panda (Ailurus fulgens) is a riddle of morphology, making it hard to tell whether it is an ursid, a procyonid, a mustelid, or a member of its own family. Previous genetic studies have given quite contradictory results as to its phylogenetic placement. Results A recently developed whole genome-based algorithm, the Whole Genome K-mer Signature algorithm was used to analyze the genomes of 28 species of Carnivora, including A. fulgens and several felid, ursid, mustelid, one mephitid species. This algorithm has the advantage of holistically using all the information in the genomes of these species. Being a genomics-based algorithm, it also reduces stochastic error to a minimum. Besides the whole genome, the mitochondrial DNA from 52 mustelids, mephitids, ursids, procyonids and A. fulgens were aligned to draw further phylogenetic inferences. The results from the whole genome study suggested that A. fulgens is a member of the mustelid clade (p = 9·10− 97). A. fulgens also separates from the mephitid Spilogala gracilis. The giant panda, Ailuropoda melanoleuca also clusters away from A. fulgens, together with other ursids (p = 1.2·10− 62). This could be due to the geographic isolation of A. fulgens from other mustelid species. However, results from the mitochondrial study as well as neighbor-joining methods based on the sequence identity matrix suggests that A. fulgens forms a monophyletic group. A Maximum Likelihood tree suggests that A. fulgens and Ursidae form a monophyletic group, although the bootstrap value is weak. Conclusions The main conclusion that we can draw from this study is that on a whole genome level A. fulgens possibly belongs to the mustelid clade, and not an ursid or a mephitid. This despite the fact that previously some researchers classified A. fulgens and A. melanoleuca as relatives. Since the genotype determines the phenotype, molecular-based classification takes precedence over morphological classifications. This affirms the results of some previous studies, which studied smaller portions of the genome. However, mitochondrial analyses based on neighbor-joining and maximum likelihood methods suggest otherwise. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-021-07531-3.
Collapse
Affiliation(s)
- Matyas Cserhati
- Independent Scholar, 2615C Muscatel Avenue, Rosemead, CA, 91770, USA.
| |
Collapse
|
36
|
Meleshko O, Martin MD, Korneliussen TS, Schröck C, Lamkowski P, Schmutz J, Healey A, Piatkowski BT, Shaw AJ, Weston DJ, Flatberg KI, Szövényi P, Hassel K, Stenøien HK. Extensive Genome-Wide Phylogenetic Discordance Is Due to Incomplete Lineage Sorting and Not Ongoing Introgression in a Rapidly Radiated Bryophyte Genus. Mol Biol Evol 2021; 38:2750-2766. [PMID: 33681996 PMCID: PMC8233498 DOI: 10.1093/molbev/msab063] [Citation(s) in RCA: 40] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
The relative importance of introgression for diversification has long been a highly disputed topic in speciation research and remains an open question despite the great attention it has received over the past decade. Gene flow leaves traces in the genome similar to those created by incomplete lineage sorting (ILS), and identification and quantification of gene flow in the presence of ILS is challenging and requires knowledge about the true phylogenetic relationship among the species. We use whole nuclear, plastid, and organellar genomes from 12 species in the rapidly radiated, ecologically diverse, actively hybridizing genus of peatmoss (Sphagnum) to reconstruct the species phylogeny and quantify introgression using a suite of phylogenomic methods. We found extensive phylogenetic discordance among nuclear and organellar phylogenies, as well as across the nuclear genome and the nodes in the species tree, best explained by extensive ILS following the rapid radiation of the genus rather than by postspeciation introgression. Our analyses support the idea of ancient introgression among the ancestral lineages followed by ILS, whereas recent gene flow among the species is highly restricted despite widespread interspecific hybridization known in the group. Our results contribute to phylogenomic understanding of how speciation proceeds in rapidly radiated, actively hybridizing species groups, and demonstrate that employing a combination of diverse phylogenomic methods can facilitate untangling complex phylogenetic patterns created by ILS and introgression.
Collapse
Affiliation(s)
- Olena Meleshko
- Department of Natural History, NTNU University Museum, Norwegian University of Science and Technology, Trondheim, Norway
| | - Michael D Martin
- Department of Natural History, NTNU University Museum, Norwegian University of Science and Technology, Trondheim, Norway
| | | | | | - Paul Lamkowski
- Institute of Botany and Landscape Ecology, University of Greifswald, Greifswald, Germany
| | - Jeremy Schmutz
- United States Department of Energy, Joint Genome Institute, Berkeley, CA, USA.,HudsonAlpha Institute for Biotechnology, Huntsville, AL, USA
| | - Adam Healey
- HudsonAlpha Institute for Biotechnology, Huntsville, AL, USA
| | | | | | - David J Weston
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, USA.,Climate Change Science Institute, Oak Ridge National Laboratory, Oak Ridge, TN, USA
| | - Kjell Ivar Flatberg
- Department of Natural History, NTNU University Museum, Norwegian University of Science and Technology, Trondheim, Norway
| | - Péter Szövényi
- Department of Systematic and Evolutionary Botany & Zurich-Basel Plant Science Center, University of Zurich, Zurich, Switzerland
| | - Kristian Hassel
- Department of Natural History, NTNU University Museum, Norwegian University of Science and Technology, Trondheim, Norway
| | - Hans K Stenøien
- Department of Natural History, NTNU University Museum, Norwegian University of Science and Technology, Trondheim, Norway
| |
Collapse
|
37
|
Hime PM, Lemmon AR, Lemmon ECM, Prendini E, Brown JM, Thomson RC, Kratovil JD, Noonan BP, Pyron RA, Peloso PLV, Kortyna ML, Keogh JS, Donnellan SC, Mueller RL, Raxworthy CJ, Kunte K, Ron SR, Das S, Gaitonde N, Green DM, Labisko J, Che J, Weisrock DW. Phylogenomics Reveals Ancient Gene Tree Discordance in the Amphibian Tree of Life. Syst Biol 2021; 70:49-66. [PMID: 32359157 PMCID: PMC7823230 DOI: 10.1093/sysbio/syaa034] [Citation(s) in RCA: 86] [Impact Index Per Article: 28.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2019] [Revised: 04/14/2020] [Accepted: 04/14/2020] [Indexed: 11/30/2022] Open
Abstract
Molecular phylogenies have yielded strong support for many parts of the amphibian Tree of Life, but poor support for the resolution of deeper nodes, including relationships among families and orders. To clarify these relationships, we provide a phylogenomic perspective on amphibian relationships by developing a taxon-specific Anchored Hybrid Enrichment protocol targeting hundreds of conserved exons which are effective across the class. After obtaining data from 220 loci for 286 species (representing 94% of the families and 44% of the genera), we estimate a phylogeny for extant amphibians and identify gene tree-species tree conflict across the deepest branches of the amphibian phylogeny. We perform locus-by-locus genealogical interrogation of alternative topological hypotheses for amphibian monophyly, focusing on interordinal relationships. We find that phylogenetic signal deep in the amphibian phylogeny varies greatly across loci in a manner that is consistent with incomplete lineage sorting in the ancestral lineage of extant amphibians. Our results overwhelmingly support amphibian monophyly and a sister relationship between frogs and salamanders, consistent with the Batrachia hypothesis. Species tree analyses converge on a small set of topological hypotheses for the relationships among extant amphibian families. These results clarify several contentious portions of the amphibian Tree of Life, which in conjunction with a set of vetted fossil calibrations, support a surprisingly younger timescale for crown and ordinal amphibian diversification than previously reported. More broadly, our study provides insight into the sources, magnitudes, and heterogeneity of support across loci in phylogenomic data sets.[AIC; Amphibia; Batrachia; Phylogeny; gene tree-species tree discordance; genomics; information theory.].
Collapse
Affiliation(s)
- Paul M Hime
- Biodiversity Institute, University of Kansas, Lawrence, KS 66045, USA
- Department of Biology, University of Kentucky, Lexington, KY 40506, USA
| | - Alan R Lemmon
- Department of Scientific Computing, Florida State University, Tallahassee, FL 32306, USA
| | | | - Elizabeth Prendini
- Division of Vertebrate Zoology: Herpetology, American Museum of Natural History, New York, NY 10024, USA
| | - Jeremy M Brown
- Department of Biological Sciences and Museum of Natural Science, Louisiana State University, Baton Rouge, LA 70803, USA
| | - Robert C Thomson
- School of Life Sciences, University of Hawai’i, Honolulu, HI 96822, USA
| | - Justin D Kratovil
- Department of Biology, University of Kentucky, Lexington, KY 40506, USA
- Department of Entomology, University of Kentucky, Lexington, KY 40546, USA
| | - Brice P Noonan
- Department of Biology, University of Mississippi, Oxford, MS 38677, USA
| | - R Alexander Pyron
- Department of Biological Sciences, The George Washington University, Washington, DC 20052, USA
| | - Pedro L V Peloso
- Division of Vertebrate Zoology: Herpetology, American Museum of Natural History, New York, NY 10024, USA
- Instituto de Ciências Biológicas, Universidade Federal do Pará, Belém, 66075-750, Brazil
| | - Michelle L Kortyna
- Department of Biological Science, Florida State University, Tallahassee, FL 32306, USA
| | - J Scott Keogh
- Division of Ecology and Evolution, Research School of Biology, The Australian National University, Canberra, 2601, Australia
| | - Stephen C Donnellan
- South Australian Museum, North Terrace, Adelaide 5000, Australia
- School of Biological Sciences, University of Adelaide, Adelaide 5005, Australia
| | | | - Christopher J Raxworthy
- Division of Vertebrate Zoology: Herpetology, American Museum of Natural History, New York, NY 10024, USA
| | - Krushnamegh Kunte
- National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bengaluru 560065, India
| | - Santiago R Ron
- Museo de Zoología, Escuela de Biología, Pontificia Universidad Católica del Ecuador, Quito, Ecuador
| | - Sandeep Das
- Forest Ecology and Biodiversity Conservation Division, Kerala Forest Research Institute, Peechi, Kerala 680653, India
| | - Nikhil Gaitonde
- National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bengaluru 560065, India
| | - David M Green
- Redpath Museum, McGill University, Montreal, Quebec H3A 0C4, Canada
| | - Jim Labisko
- The Durrell Institute of Conservation and Ecology, School of Anthropology and Conservation, The University of Kent, Canterbury, Kent, CT2 7NR, UK
- Island Biodiversity and Conservation Centre, University of Seychelles, PO Box 1348, Anse Royale, Mahé, Seychelles
| | - Jing Che
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Kunming 650223, China
- Center for Excellence in Animal Evolution and Genetics, Chinese Academy of Sciences, Kunming 650223, China
| | - David W Weisrock
- Department of Biology, University of Kentucky, Lexington, KY 40506, USA
| |
Collapse
|
38
|
Hibbins MS, Gibson MJS, Hahn MW. Determining the probability of hemiplasy in the presence of incomplete lineage sorting and introgression. eLife 2020; 9:e63753. [PMID: 33345772 PMCID: PMC7800383 DOI: 10.7554/elife.63753] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2020] [Accepted: 12/18/2020] [Indexed: 12/11/2022] Open
Abstract
The incongruence of character states with phylogenetic relationships is often interpreted as evidence of convergent evolution. However, trait evolution along discordant gene trees can also generate these incongruences - a phenomenon known as hemiplasy. Classic comparative methods do not account for discordance, resulting in incorrect inferences about the number, timing, and direction of trait transitions. Biological sources of discordance include incomplete lineage sorting (ILS) and introgression, but only ILS has received theoretical consideration in the context of hemiplasy. Here, we present a model that shows introgression makes hemiplasy more likely, such that methods that account for ILS alone will be conservative. We also present a method and software (HeIST) for making statistical inferences about the probability of hemiplasy and homoplasy in large datasets that contain both ILS and introgression. We apply our methods to two empirical datasets, finding that hemiplasy is likely to contribute to the observed trait incongruences in both.
Collapse
Affiliation(s)
- Mark S Hibbins
- Department of Biology, Indiana UniversityBloomingtonUnited States
| | | | - Matthew W Hahn
- Department of Biology, Indiana UniversityBloomingtonUnited States
- Department of Computer Science, Indiana UniversityBloomingtonUnited States
| |
Collapse
|
39
|
Hosegood J, Humble E, Ogden R, de Bruyn M, Creer S, Stevens GMW, Abudaya M, Bassos-Hull K, Bonfil R, Fernando D, Foote AD, Hipperson H, Jabado RW, Kaden J, Moazzam M, Peel LR, Pollett S, Ponzo A, Poortvliet M, Salah J, Senn H, Stewart JD, Wintner S, Carvalho G. Phylogenomics and species delimitation for effective conservation of manta and devil rays. Mol Ecol 2020; 29:4783-4796. [PMID: 33164287 DOI: 10.1111/mec.15683] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2020] [Revised: 09/25/2020] [Accepted: 09/30/2020] [Indexed: 02/06/2023]
Abstract
Practical biodiversity conservation relies on delineation of biologically meaningful units. Manta and devil rays (Mobulidae) are threatened worldwide, yet morphological similarities and a succession of recent taxonomic changes impede the development of an effective conservation strategy. Here, we generate genome-wide single nucleotide polymorphism (SNP) data from a geographically and taxonomically representative set of manta and devil ray samples to reconstruct phylogenetic relationships and evaluate species boundaries under the general lineage concept. We show that nominal species units supported by alternative data sources constitute independently evolving lineages, and find robust evidence for a putative new species of manta ray in the Gulf of Mexico. Additionally, we uncover substantial incomplete lineage sorting indicating that rapid speciation together with standing variation in ancestral populations has driven phylogenetic uncertainty within Mobulidae. Finally, we detect cryptic diversity in geographically distinct populations, demonstrating that management below the species level may be warranted in certain species. Overall, our study provides a framework for molecular genetic species delimitation that is relevant to wide-ranging taxa of conservation concern, and highlights the potential for genomic data to support effective management, conservation and law enforcement strategies.
Collapse
Affiliation(s)
- Jane Hosegood
- Molecular Ecology and Fisheries Genetics Laboratory, Bangor University, Bangor, UK.,The Manta Trust, Catemwood House, Norwood Lane, Corscombe, Dorset, DT2 0NT, UK.,NERC Biomolecular Analysis Facility, Department of Animal and Plant Sciences, University of Sheffield, Sheffield, UK
| | - Emily Humble
- The Manta Trust, Catemwood House, Norwood Lane, Corscombe, Dorset, DT2 0NT, UK.,Royal (Dick) School of Veterinary Studies and the Roslin Institute, University of Edinburgh, Edinburgh, UK
| | - Rob Ogden
- Royal (Dick) School of Veterinary Studies and the Roslin Institute, University of Edinburgh, Edinburgh, UK.,TRACE Wildlife Forensics Network, Edinburgh, UK
| | - Mark de Bruyn
- Molecular Ecology and Fisheries Genetics Laboratory, Bangor University, Bangor, UK.,School of Life and Environmental Sciences, The University of Sydney, Sydney, NSW, Australia
| | - Simon Creer
- Molecular Ecology and Fisheries Genetics Laboratory, Bangor University, Bangor, UK
| | - Guy M W Stevens
- The Manta Trust, Catemwood House, Norwood Lane, Corscombe, Dorset, DT2 0NT, UK
| | | | - Kim Bassos-Hull
- Mote Marine Laboratory, The Center for Shark Research, Sarasota, FL, USA
| | | | - Daniel Fernando
- The Manta Trust, Catemwood House, Norwood Lane, Corscombe, Dorset, DT2 0NT, UK.,Department of Biology and Environmental Science, Linnaeus University, Kalmar, Sweden.,Blue Resources Trust, Colombo, Sri Lanka
| | - Andrew D Foote
- Molecular Ecology and Fisheries Genetics Laboratory, Bangor University, Bangor, UK.,Department of Natural History, Norwegian University of Science and Technology (NTNU), University Museum, Trondheim, Norway
| | - Helen Hipperson
- NERC Biomolecular Analysis Facility, Department of Animal and Plant Sciences, University of Sheffield, Sheffield, UK
| | | | - Jennifer Kaden
- RZSS WildGenes Lab, Royal Zoological Society of Scotland, Edinburgh, UK
| | | | - Lauren R Peel
- The Manta Trust, Catemwood House, Norwood Lane, Corscombe, Dorset, DT2 0NT, UK.,School of Biological Sciences, University of Western Australia, Crawley, WA, Australia.,The Australian Institute of Marine Science, Crawley, WA, Australia.,Save Our Seas Foundation - D'Arros Research Centre, Geneva, Switzerland
| | - Stephen Pollett
- The Manta Trust, Catemwood House, Norwood Lane, Corscombe, Dorset, DT2 0NT, UK
| | - Alessandro Ponzo
- Large Marine Vertebrates Research Institute Philippines, Jagna, Philippines
| | | | - Jehad Salah
- Ministry of Agriculture Directorate General of Fisheries, Gaza City, Palestine
| | - Helen Senn
- RZSS WildGenes Lab, Royal Zoological Society of Scotland, Edinburgh, UK
| | - Joshua D Stewart
- The Manta Trust, Catemwood House, Norwood Lane, Corscombe, Dorset, DT2 0NT, UK
| | - Sabine Wintner
- KwaZulu-Natal Sharks Board, Umhlanga Rocks, South Africa.,School of Life Sciences, University of KwaZulu-Natal, Durban, South Africa
| | - Gary Carvalho
- Molecular Ecology and Fisheries Genetics Laboratory, Bangor University, Bangor, UK
| |
Collapse
|
40
|
A phylogenomic study of Steganinae fruit flies (Diptera: Drosophilidae): strong gene tree heterogeneity and evidence for monophyly. BMC Evol Biol 2020; 20:141. [PMID: 33138771 PMCID: PMC7607883 DOI: 10.1186/s12862-020-01703-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2020] [Accepted: 10/19/2020] [Indexed: 12/18/2022] Open
Abstract
BACKGROUND The Drosophilidae family is traditionally divided into two subfamilies: Drosophilinae and Steganinae. This division is based on morphological characters, and the two subfamilies have been treated as monophyletic in most of the literature, but some molecular phylogenies have suggested Steganinae to be paraphyletic. To test the paraphyletic-Steganinae hypothesis, here, we used genomic sequences of eight Drosophilidae (three Steganinae and five Drosophilinae) and two Ephydridae (outgroup) species and inferred the phylogeny for the group based on a dataset of 1,028 orthologous genes present in all species (> 1,000,000 bp). This dataset includes three genera that broke the monophyly of the subfamilies in previous works. To investigate possible biases introduced by small sample sizes and automatic gene annotation, we used the same methods to infer species trees from a set of 10 manually annotated genes that are commonly used in phylogenetics. RESULTS Most of the 1,028 gene trees depicted Steganinae as paraphyletic with distinct topologies, but the most common topology depicted it as monophyletic (43.7% of the gene trees). Despite the high levels of gene tree heterogeneity observed, species tree inference in ASTRAL, in PhyloNet, and with the concatenation approach strongly supported the monophyly of both subfamilies for the 1,028-gene dataset. However, when using the concatenation approach to infer a species tree from the smaller set of 10 genes, we recovered Steganinae as a paraphyletic group. The pattern of gene tree heterogeneity was asymmetrical and thus could not be explained solely by incomplete lineage sorting (ILS). CONCLUSIONS Steganinae was clearly a monophyletic group in the dataset that we analyzed. In addition to ILS, gene tree discordance was possibly the result of introgression, suggesting complex branching processes during the early evolution of Drosophilidae with short speciation intervals and gene flow. Our study highlights the importance of genomic data in elucidating contentious phylogenetic relationships and suggests that phylogenetic inference for drosophilids based on small molecular datasets should be performed cautiously. Finally, we suggest an approach for the correction and cleaning of BUSCO-derived genomic datasets that will be useful to other researchers planning to use this tool for phylogenomic studies.
Collapse
|
41
|
Gabaldón T. Patterns and impacts of nonvertical evolution in eukaryotes: a paradigm shift. Ann N Y Acad Sci 2020; 1476:78-92. [PMID: 32860228 PMCID: PMC7589212 DOI: 10.1111/nyas.14471] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2020] [Revised: 07/19/2020] [Accepted: 07/27/2020] [Indexed: 12/14/2022]
Abstract
Evolution of eukaryotic species and their genomes has been traditionally understood as a vertical process in which genetic material is transmitted from parents to offspring along a lineage, and in which genetic exchange is restricted within species boundaries. However, mounting evidence from comparative genomics indicates that this paradigm is often violated. Horizontal gene transfer and mating between diverged lineages blur species boundaries and challenge the reconstruction of evolutionary histories of species and their genomes. Nonvertical evolution might be more restricted in eukaryotes than in prokaryotes, yet it is not negligible and can be common in certain groups. Recognition of such processes brings about the need to incorporate this complexity into our models, as well as to conceptually reframe eukaryotic diversity and evolution. Here, I review the recent work from genomics studies that supports the effects of nonvertical modes of evolution including introgression, hybridization, and horizontal gene transfer in different eukaryotic groups. I then discuss emerging patterns and effects, illustrated by specific examples, that support the conclusion that nonvertical processes are often at the root of important evolutionary transitions and adaptations. I will argue that a paradigm shift is needed to naturally accommodate nonvertical processes in eukaryotic evolution.
Collapse
Affiliation(s)
- Toni Gabaldón
- Barcelona Supercomputing Centre (BCS-CNS), Barcelona, Spain.,Institute for Research in Biomedicine (IRB), The Barcelona Institute of Science and Technology (BIST), Barcelona, Spain.,Institució Catalana de Recerca i Estudis Avançats (ICREA), Barcelona, Spain
| |
Collapse
|
42
|
Yin J, Zhang C, Mirarab S. ASTRAL-MP: scaling ASTRAL to very large datasets using randomization and parallelization. Bioinformatics 2020; 35:3961-3969. [PMID: 30903685 DOI: 10.1093/bioinformatics/btz211] [Citation(s) in RCA: 34] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2018] [Revised: 03/12/2019] [Accepted: 03/21/2019] [Indexed: 01/11/2023] Open
Abstract
MOTIVATION Evolutionary histories can change from one part of the genome to another. The potential for discordance between the gene trees has motivated the development of summary methods that reconstruct a species tree from an input collection of gene trees. ASTRAL is a widely used summary method and has been able to scale to relatively large datasets. However, the size of genomic datasets is quickly growing. Despite its relative efficiency, the current single-threaded implementation of ASTRAL is falling behind the data growth trends is not able to analyze the largest available datasets in a reasonable time. RESULTS ASTRAL uses dynamic programing and is not trivially parallel. In this paper, we introduce ASTRAL-MP, the first version of ASTRAL that can exploit parallelism and also uses randomization techniques to speed up some of its steps. Importantly, ASTRAL-MP can take advantage of not just multiple CPU cores but also one or several graphics processing units (GPUs). The ASTRAL-MP code scales very well with increasing CPU cores, and its GPU version, implemented in OpenCL, can have up to 158× speedups compared to ASTRAL-III. Using GPUs and multiple cores, ASTRAL-MP is able to analyze datasets with 10 000 species or datasets with more than 100 000 genes in <2 days. AVAILABILITY AND IMPLEMENTATION ASTRAL-MP is available at https://github.com/smirarab/ASTRAL/tree/MP. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- John Yin
- Department of Mathematics, University of California at San Diego, La Jolla, CA, USA
| | - Chao Zhang
- Bioinformatics and Systems Biology, University of California at San Diego, La Jolla, CA, USA
| | - Siavash Mirarab
- Department of Electrical and Computer Engineering, University of California at San Diego, La Jolla, CA, USA
| |
Collapse
|
43
|
Moura AE, Shreves K, Pilot M, Andrews KR, Moore DM, Kishida T, Möller L, Natoli A, Gaspari S, McGowen M, Chen I, Gray H, Gore M, Culloch RM, Kiani MS, Willson MS, Bulushi A, Collins T, Baldwin R, Willson A, Minton G, Ponnampalam L, Hoelzel AR. Phylogenomics of the genus Tursiops and closely related Delphininae reveals extensive reticulation among lineages and provides inference about eco-evolutionary drivers. Mol Phylogenet Evol 2020; 146:106756. [DOI: 10.1016/j.ympev.2020.106756] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2019] [Revised: 01/02/2020] [Accepted: 01/28/2020] [Indexed: 12/30/2022]
|
44
|
Hsieh P, Vollger MR, Dang V, Porubsky D, Baker C, Cantsilieris S, Hoekzema K, Lewis AP, Munson KM, Sorensen M, Kronenberg ZN, Murali S, Nelson BJ, Chiatante G, Maggiolini FAM, Blanché H, Underwood JG, Antonacci F, Deleuze JF, Eichler EE. Adaptive archaic introgression of copy number variants and the discovery of previously unknown human genes. Science 2020; 366:366/6463/eaax2083. [PMID: 31624180 DOI: 10.1126/science.aax2083] [Citation(s) in RCA: 50] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2019] [Revised: 07/05/2019] [Accepted: 09/12/2019] [Indexed: 01/01/2023]
Abstract
Copy number variants (CNVs) are subject to stronger selective pressure than single-nucleotide variants, but their roles in archaic introgression and adaptation have not been systematically investigated. We show that stratified CNVs are significantly associated with signatures of positive selection in Melanesians and provide evidence for adaptive introgression of large CNVs at chromosomes 16p11.2 and 8p21.3 from Denisovans and Neanderthals, respectively. Using long-read sequence data, we reconstruct the structure and complex evolutionary history of these polymorphisms and show that both encode positively selected genes absent from most human populations. Our results collectively suggest that large CNVs originating in archaic hominins and introgressed into modern humans have played an important role in local population adaptation and represent an insufficiently studied source of large-scale genetic variation.
Collapse
Affiliation(s)
- PingHsun Hsieh
- Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA
| | - Mitchell R Vollger
- Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA
| | - Vy Dang
- Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA
| | - David Porubsky
- Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA
| | - Carl Baker
- Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA
| | - Stuart Cantsilieris
- Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA
| | - Kendra Hoekzema
- Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA
| | - Alexandra P Lewis
- Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA
| | - Katherine M Munson
- Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA
| | - Melanie Sorensen
- Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA
| | - Zev N Kronenberg
- Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA
| | - Shwetha Murali
- Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA.,Howard Hughes Medical Institute, University of Washington, Seattle, WA, USA
| | - Bradley J Nelson
- Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA
| | - Giorgia Chiatante
- Dipartimento di Biologia, Università degli Studi di Bari "Aldo Moro," Bari, Italy
| | | | - Hélène Blanché
- Fondation Jean Dausset-Centre d'Etude du Polymorphisme Humain, Paris, France
| | - Jason G Underwood
- Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA.,Pacific Biosciences (PacBio) of California, Inc., Menlo Park, CA, USA
| | - Francesca Antonacci
- Dipartimento di Biologia, Università degli Studi di Bari "Aldo Moro," Bari, Italy
| | | | - Evan E Eichler
- Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA. .,Howard Hughes Medical Institute, University of Washington, Seattle, WA, USA
| |
Collapse
|
45
|
Zenker MM, Specht A, Fonseca VG. Assessing insect biodiversity with automatic light traps in Brazil: Pearls and pitfalls of metabarcoding samples in preservative ethanol. Ecol Evol 2020; 10:2352-2366. [PMID: 32184986 PMCID: PMC7069332 DOI: 10.1002/ece3.6042] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2019] [Revised: 12/23/2019] [Accepted: 12/27/2019] [Indexed: 12/24/2022] Open
Abstract
Automated species identification based on data produced with metabarcoding offers an alternative for assessing biodiversity of bulk insect samples obtained with traps. We used a standard two-step PCR approach to amplify a 313 bp fragment of the barcoding region of the mitochondrial COI gene. The PCR products were sequenced on an Illumina MiSeq platform, and the OTUs production and taxonomic identifications were performed with a customized pipeline and database. The DNA used in the PCR procedures was extracted directly from the preservative ethanol of bulk insect samples obtained with automatic light traps in 12 sampling areas located in different biomes of Brazil, during wet and dry seasons. Agricultural field and forest edge habitats were collected for all sampling areas. A total of 119 insect OTUs and nine additional OTUs assigned to other arthropod taxa were obtained at a ≥97% sequence similarity level. The alpha and beta diversity analyses comparing biomes, habitats, and seasons were mostly inconclusive, except for a significant difference in beta diversity between biomes. In this study, we were able to metabarcode and HTS adult insects from their preservative medium. Notwithstanding, our results underrepresent the true magnitude of insect diversity expected from samples obtained with automatic light traps in Brazil. Although biological and technical factors might have impacted our results, measures to optimize and standardize eDNA HTS should be in place to improve taxonomic coverage of samples of unknown diversity and stored in suboptimal conditions, which is the case of most eDNA samples.
Collapse
Affiliation(s)
- Mauricio M. Zenker
- Zoological Research Museum Alexander KoenigBonnGermany
- Present address:
Rua Eulo MaroniSão PauloBrazil
| | - Alexandre Specht
- Embrapa CerradosBrasiliaBrazil
- Present address:
Embrapa CerradosPlanaltinaFederal DistrictBrazil
| | - Vera G. Fonseca
- Zoological Research Museum Alexander KoenigBonnGermany
- Present address:
Centre for EnvironmentFisheries and Aquaculture Science (Cefas)WeymouthUK
| |
Collapse
|
46
|
Koenen EJM, Ojeda DI, Steeves R, Migliore J, Bakker FT, Wieringa JJ, Kidner C, Hardy OJ, Pennington RT, Bruneau A, Hughes CE. Large-scale genomic sequence data resolve the deepest divergences in the legume phylogeny and support a near-simultaneous evolutionary origin of all six subfamilies. THE NEW PHYTOLOGIST 2020; 225:1355-1369. [PMID: 31665814 PMCID: PMC6972672 DOI: 10.1111/nph.16290] [Citation(s) in RCA: 66] [Impact Index Per Article: 16.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/28/2019] [Accepted: 09/14/2019] [Indexed: 05/02/2023]
Abstract
Phylogenomics is increasingly used to infer deep-branching relationships while revealing the complexity of evolutionary processes such as incomplete lineage sorting, hybridization/introgression and polyploidization. We investigate the deep-branching relationships among subfamilies of the Leguminosae (or Fabaceae), the third largest angiosperm family. Despite their ecological and economic importance, a robust phylogenetic framework for legumes based on genome-scale sequence data is lacking. We generated alignments of 72 chloroplast genes and 7621 homologous nuclear-encoded proteins, for 157 and 76 taxa, respectively. We analysed these with maximum likelihood, Bayesian inference, and a multispecies coalescent summary method, and evaluated support for alternative topologies across gene trees. We resolve the deepest divergences in the legume phylogeny despite lack of phylogenetic signal across all chloroplast genes and the majority of nuclear genes. Strongly supported conflict in the remainder of nuclear genes is suggestive of incomplete lineage sorting. All six subfamilies originated nearly simultaneously, suggesting that the prevailing view of some subfamilies as 'basal' or 'early-diverging' with respect to others should be abandoned, which has important implications for understanding the evolution of legume diversity and traits. Our study highlights the limits of phylogenetic resolution in relation to rapid successive speciation.
Collapse
Affiliation(s)
- Erik J. M. Koenen
- Department of Systematic and Evolutionary BotanyUniversity of ZurichZollikerstrasse 107CH‐8008ZurichSwitzerland
| | - Dario I. Ojeda
- Service Évolution Biologique et ÉcologieFaculté des SciencesUniversité Libre de BruxellesAvenue Franklin Roosevelt 501050BrusselsBelgium
- Norwegian Institute of Bioeconomy ResearchHøgskoleveien 81433ÅsNorway
| | - Royce Steeves
- Institut de Recherche en Biologie Végétale and Département de Sciences BiologiquesUniversité de Montréal4101 Sherbrooke St EMontrealQCH1X 2B2Canada
- Fisheries & Oceans CanadaGulf Fisheries Center343 Université AveMonctonNBE1C 5K4Canada
| | - Jérémy Migliore
- Service Évolution Biologique et ÉcologieFaculté des SciencesUniversité Libre de BruxellesAvenue Franklin Roosevelt 501050BrusselsBelgium
| | - Freek T. Bakker
- Biosystematics GroupWageningen UniversityDroevendaalsesteeg 16708 PBWageningenthe Netherlands
| | - Jan J. Wieringa
- Naturalis Biodiversity Center, LeidenDarwinweg 22333 CRLeidenthe Netherlands
| | - Catherine Kidner
- Royal Botanic Gardens Edinburgh20a Inverleith RowEdinburghEH3 5LRUK
- School of Biological SciencesUniversity of EdinburghKing's Buildings, Mayfield RdEdinburghEH9 3JUUK
| | - Olivier J. Hardy
- Service Évolution Biologique et ÉcologieFaculté des SciencesUniversité Libre de BruxellesAvenue Franklin Roosevelt 501050BrusselsBelgium
| | - R. Toby Pennington
- Royal Botanic Gardens Edinburgh20a Inverleith RowEdinburghEH3 5LRUK
- GeographyUniversity of ExeterAmory Building, Rennes DriveExeterEX4 4RJUK
| | - Anne Bruneau
- Institut de Recherche en Biologie Végétale and Département de Sciences BiologiquesUniversité de Montréal4101 Sherbrooke St EMontrealQCH1X 2B2Canada
| | - Colin E. Hughes
- Department of Systematic and Evolutionary BotanyUniversity of ZurichZollikerstrasse 107CH‐8008ZurichSwitzerland
| |
Collapse
|
47
|
Phylogeny and highland adaptation of Chinese species in Allium section Daghestanica (Amaryllidaceae) revealed by transcriptome sequencing. Mol Phylogenet Evol 2020; 146:106737. [PMID: 31982455 DOI: 10.1016/j.ympev.2020.106737] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2019] [Revised: 01/20/2020] [Accepted: 01/21/2020] [Indexed: 01/27/2023]
Abstract
Allium L. is one of the largest monocotyledonous genera with extensive distribution in the Northern Hemisphere. The fundamental phylogenies of Allium have been investigated using many morphological and molecular characters. However, the morphological characters may not agree with the molecular results in some Allium groups or sections (such as the Chinese Allium section Daghestanica), which may result in ambiguous species relationships and hinder further evolutionary and adaptive researches. Here, transcriptome sequences of the six Chinese endemics from Allium section Daghestanica were collected, with their single-copy genes (SCGs) were extracted. The interspecies relationships were analyzed using concatenation and coalescent methods. The branch-site model (BSM) was conducted to detect the positively selected genes (PSGs) in five highland species of this section. Based on 1644, 1281 and 1580 SCGs in flowers, leaves, and flowers-leaves combination respectively, a robust consistent and well-resolved phylogeny was generated from the concatenation method. Strong conflicts among individual gene trees were detected in the coalescent method, and morphological characters were incongruent with molecular relationships to some degree. Many PSGs were involved in responses of various stresses and stimuli (e.g. hypoxia, low temperature, aridity), DNA repair, metabolism, nutrient or energy intake, photosynthesis, and signal transduction. Our study revealed a clear interspecies relationship of Chinese endemics in Allium section Daghestanica and suggested that the discordance between morphological characters and molecular relationships might result from that the former are more susceptible to convergence compared with the latter. PSGs detected in our study may provide some insights into highland adaptation in Allium species.
Collapse
|
48
|
Skim-Sequencing Reveals the Likely Origin of the Enigmatic Endangered Sunflower Helianthus schweinitzii. Genes (Basel) 2019; 10:genes10121040. [PMID: 31847493 PMCID: PMC6947582 DOI: 10.3390/genes10121040] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2019] [Revised: 12/05/2019] [Accepted: 12/12/2019] [Indexed: 12/13/2022] Open
Abstract
Resolving the origin of endangered taxa is an essential component of conservation. This information can be used to guide efforts of bolstering genetic diversity, and also enables species recovery and future evolutionary studies. Here, we used low-coverage whole genome sequencing to clarify the origin of Helianthus schweinitzii, an endangered tetraploid sunflower that is endemic to the Piedmont Plateau in the eastern United States. We surveyed four accessions representing four populations of H. schweinitzii and 38 accessions of six purported parental species. Using de novo approaches, we assembled 87,004 bp of the chloroplast genome and 6770 bp of the nuclear 35S rDNA. Phylogenetic reconstructions based on the chloroplast genome revealed no reciprocal monophyly of taxa. In contrast, nuclear rDNA data strongly supported the currently accepted sections of the genus Helianthus. Information from combined cpDNA and rDNA provided evidence that H. schweinitzii is likely an allo-tetraploid that formed as a result of hybridization between the diploids Helianthus giganteus and Helianthus microcephalus.
Collapse
|
49
|
Choo A, Nguyen TNM, Ward CM, Chen IY, Sved J, Shearman D, Gilchrist AS, Crisp P, Baxter SW. Identification of Y-chromosome scaffolds of the Queensland fruit fly reveals a duplicated gyf gene paralogue common to many Bactrocera pest species. INSECT MOLECULAR BIOLOGY 2019; 28:873-886. [PMID: 31150140 DOI: 10.1111/imb.12602] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/19/2019] [Revised: 05/22/2019] [Accepted: 05/28/2019] [Indexed: 06/09/2023]
Abstract
Bactrocera tryoni (Queensland fruit fly) are polyphagous horticultural pests of eastern Australia. Heterogametic males contain a sex-determining Y-chromosome thought to be gene poor and repetitive. Here, we report 39 Y-chromosome scaffolds (~700 kb) from B. tryoni identified using genotype-by-sequencing data and whole-genome resequencing. Male diagnostic PCR assays validated eight Y-scaffolds, and one (Btry4096) contained a novel gene with five exons that encode a predicted 575 amino acid protein. The Y-gene, referred to as typo-gyf, is a truncated Y-chromosome paralogue of X-chromosome gene gyf (1773 aa). The Y-chromosome contained ~41 copies of typo-gyf, and expression occurred in male flies and embryos. Analysis of 13 tephritid transcriptomes confirmed typo-gyf expression in six additional Bactrocera species, including Bactrocera latifrons, Bactrocera dorsalis and Bactrocera zonata. Molecular dating estimated typo-gyf evolved within the past 8.02 million years (95% highest posterior density 10.56-5.52 million years), after the split with Bactrocera oleae. Phylogenetic analysis also highlighted complex evolutionary histories among several Bactrocera species, as discordant nuclear (116 genes) and mitochondrial (13 genes) topologies were observed. B. tryoni Y-sequences may provide useful sites for future transgene insertions, and typo-gyf could act as a Y-chromosome diagnostic marker for many Bactrocera species, although its function is unknown.
Collapse
Affiliation(s)
- Amanda Choo
- School of Biological Sciences, University of Adelaide, Adelaide, South Australia, Australia
| | - Thu N M Nguyen
- School of Biological Sciences, University of Adelaide, Adelaide, South Australia, Australia
| | - Christopher M Ward
- School of Biological Sciences, University of Adelaide, Adelaide, South Australia, Australia
| | - Isabel Y Chen
- School of Biological Sciences, University of Adelaide, Adelaide, South Australia, Australia
- South Australian Research and Development Institute, Adelaide, South Australia, Australia
| | - John Sved
- Evolution and Ecology Research Centre, University of New South Wales, Sydney, New South Wales, Australia
| | - Deborah Shearman
- Evolution and Ecology Research Centre, University of New South Wales, Sydney, New South Wales, Australia
| | - Anthony S Gilchrist
- Evolution and Ecology Research Centre, University of New South Wales, Sydney, New South Wales, Australia
| | - Peter Crisp
- South Australian Research and Development Institute, Adelaide, South Australia, Australia
| | - Simon W Baxter
- School of Biological Sciences, University of Adelaide, Adelaide, South Australia, Australia
| |
Collapse
|
50
|
Zhang ML, Li ML, Ayoola AO, Murphy RW, Wu DD, Shao Y. Conserved sequences identify the closest living relatives of primates. Zool Res 2019; 40:532-540. [PMID: 31393097 PMCID: PMC6822925 DOI: 10.24272/j.issn.2095-8137.2019.057] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Abstract
Elucidating the closest living relatives of extant primates is essential for fully understanding important biological processes related to the genomic and phenotypic evolution of primates, especially of humans. However, the phylogenetic placement of these primate relatives remains controversial, with three primary hypotheses currently espoused based on morphological and molecular evidence. In the present study, we used two algorithms to analyze differently partitioned genomic datasets consisting of 45.4 Mb of conserved non-coding elements and 393 kb of concatenated coding sequences to test these hypotheses. We assessed different genomic histories and compared with other molecular studies found solid support for colugos being the closest living relatives of primates. Our phylogeny showed Cercopithecinae to have low levels of nucleotide divergence, especially for Papionini, and gibbons to have a high rate of divergence. The MCMCtree comprehensively updated divergence dates of early evolution of Primatomorpha and Primates.
Collapse
Affiliation(s)
- Mei-Ling Zhang
- Department of Acute Infectious Diseases Control and Prevention, Yunnan Provincial Centre for Disease Control and Prevention, Kunming Yunnan 650022, China
| | - Ming-Li Li
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming Yunnan 650223, China.,Kunming College of Life Science, University of the Chinese Academy of Sciences, Kunming Yunnan 650223, China
| | - Adeola Oluwakemi Ayoola
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming Yunnan 650223, China.,Kunming College of Life Science, University of the Chinese Academy of Sciences, Kunming Yunnan 650223, China
| | - Robert W Murphy
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming Yunnan 650223, China.,Centre for Biodiversity and Conservation Biology, Royal Ontario Museum, Toronto M5S 2C6, Canada
| | - Dong-Dong Wu
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming Yunnan 650223, China, E-mail:
| | - Yong Shao
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming Yunnan 650223, China, E-mail:
| |
Collapse
|