1
|
Wu J, Duan C, Lan L, Chen W, Mao X. Sex Differences in Cochlear Transcriptomes in Horseshoe Bats. Animals (Basel) 2024; 14:1177. [PMID: 38672325 PMCID: PMC11047584 DOI: 10.3390/ani14081177] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2024] [Revised: 04/11/2024] [Accepted: 04/11/2024] [Indexed: 04/28/2024] Open
Abstract
Sexual dimorphism of calls is common in animals, whereas studies on the molecular basis underlying this phenotypic variation are still scarce. In this study, we used comparative transcriptomics of cochlea to investigate the sex-related difference in gene expression and alternative splicing in four Rhinolophus taxa. Based on 31 cochlear transcriptomes, we performed differential gene expression (DGE) and alternative splicing (AS) analyses between the sexes in each taxon. Consistent with the degree of difference in the echolocation pulse frequency between the sexes across the four taxa, we identified the largest number of differentially expressed genes (DEGs) and alternatively spliced genes (ASGs) in R. sinicus. However, we also detected multiple DEGs and ASGs in taxa without sexual differences in echolocation pulse frequency, suggesting that these genes might be related to other parameters of echolocation pulse rather than the frequency component. Some DEGs and ASGs are related to hearing loss or deafness genes in human or mice and they can be considered to be candidates associated with the sexual differences of echolocation pulse in bats. We also detected more than the expected overlap of DEGs and ASGs in two taxa. Overall, our current study supports the important roles of both DGE and AS in generating or maintaining sexual differences in animals.
Collapse
Affiliation(s)
| | | | | | | | - Xiuguang Mao
- School of Ecological and Environmental Sciences, East China Normal University, Shanghai 200062, China; (J.W.); (C.D.); (L.L.); (W.C.)
| |
Collapse
|
2
|
Mitonuclear mismatch alters nuclear gene expression in naturally introgressed Rhinolophus bats. Front Zool 2021; 18:42. [PMID: 34488775 PMCID: PMC8419968 DOI: 10.1186/s12983-021-00424-x] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2021] [Accepted: 07/20/2021] [Indexed: 01/23/2023] Open
Abstract
Background Mitochondrial function involves the interplay between mitochondrial and nuclear genomes. Such mitonuclear interactions can be disrupted by the introgression of mitochondrial DNA between taxa or divergent populations. Previous studies of several model systems (e.g. Drosophila) indicate that the disruption of mitonuclear interactions, termed mitonuclear mismatch, can alter nuclear gene expression, yet few studies have focused on natural populations. Results Here we study a naturally introgressed population in the secondary contact zone of two subspecies of the intermediate horseshoe bat (Rhinolophus affinis), in which individuals possess either mitonuclear matched or mismatched genotypes. We generated transcriptome data for six tissue types from five mitonuclear matched and five mismatched individuals. Our results revealed strong tissue-specific effects of mitonuclear mismatch on nuclear gene expression with the largest effect seen in pectoral muscle. Moreover, consistent with the hypothesis that genes associated with the response to oxidative stress may be upregulated in mitonuclear mismatched individuals, we identified several such gene candidates, including DNASE1L3, GPx3 and HSPB6 in muscle, and ISG15 and IFI6 in heart. Conclusion Our study reveals how mitonuclear mismatch arising from introgression in natural populations is likely to have fitness consequences. Underlying the processes that maintain mitonuclear discordance is a step forward to understand the role of mitonuclear interactions in population divergence and speciation. Supplementary Information The online version contains supplementary material available at 10.1186/s12983-021-00424-x.
Collapse
|
3
|
Chen W, Mao X. Extensive alternative splicing triggered by mitonuclear mismatch in naturally introgressed Rhinolophus bats. Ecol Evol 2021; 11:12003-12010. [PMID: 34522356 PMCID: PMC8427577 DOI: 10.1002/ece3.7966] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2021] [Revised: 07/07/2021] [Accepted: 07/16/2021] [Indexed: 11/16/2022] Open
Abstract
Mitochondrial function needs strong interactions of mitochondrial and nuclear (mitonuclear) genomes, which can be disrupted by mitonuclear mismatch due to mitochondrial DNA (mtDNA) introgression between two formerly isolated populations or taxa. This mitonuclear disruption may cause severe cellular stress in mismatched individuals. Gene expression changes and alternative splicing (AS) are two important transcriptional regulations to respond to environmental or cellular stresses. We previously identified a naturally introgressed population in the intermediate horseshoe bat (Rhinolophus affinis). Individuals from this population belong to R. a. himalayanus and share almost identical nuclear genetic background; however, some of them had mtDNA from another subspecies (R. a. macrurus). With this unique natural system, we examined gene expression changes in six tissues between five mitonuclear mismatched and five matched individuals. A small number of differentially expressed genes (DEGs) were identified, and functional enrichment analysis revealed that most DEGs were related to immune response although some may be involved in response to oxidative stress. In contrast, we identified extensive AS events and alternatively spliced genes (ASGs) between mismatched and matched individuals. Functional enrichment analysis revealed that multiple ASGs were directly or indirectly associated with energy production in mitochondria which is vital for survival of organism. To our knowledge, this is the first study to examine the role of AS in responding to cellular stress caused by mitonuclear mismatch in natural populations. Our results suggest that AS may play a more important role than gene expression regulation in responding to severe environmental or cellular stresses.
Collapse
Affiliation(s)
- Wenli Chen
- School of Ecological and Environmental SciencesEast China Normal UniversityShanghaiChina
| | - Xiuguang Mao
- School of Ecological and Environmental SciencesEast China Normal UniversityShanghaiChina
- Institute of Eco‐Chongming (IEC)East China Normal UniversityShanghaiChina
| |
Collapse
|
4
|
Genetic Introgression and Morphological Variation in Naked-Back Bats (Chiroptera: Mormoopidae: Pteronotus Species) along Their Contact Zone in Central America. DIVERSITY 2021. [DOI: 10.3390/d13050194] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Two sibling bare-backed bat species (Pteronotus fulvus and P. gymnonotus) have been traditionally differentiated by their size. However, intermediate specimens between the two species have been found in sympatric populations along southern Mexico and it has been suggested that they may be the outcome of a hybridization process between the two species. We used one mitochondrial (COI), three nuclear markers (PRKCL, STAT5A and RAG2) and 13 microsatellites to explore the evolutionary relationships between these two species and elucidate whether the intermediate morphotypes correspond to hybrid individuals. These markers have been analyzed in sympatric and allopatric populations of the two species plus the closely related species Pteronotus davyi. We confirmed the species-level differentiation of the three lineages (P. fulvus, P. davyi and P. gymnonotus), but the phylogenetic hypotheses suggested by the nuclear and mitochondrial markers were discordant. We confirm that the discordance between markers is due to genetic introgression through the mitochondrial capture of P. fulvus in P. gymnonotus populations. Such introgression was found in all P. gymnonotus specimens across its sympatric distribution range (Mexico to Costa Rica) and is related to expansion/retraction species distribution pulses associated with changes in forest distribution during the Quaternary climate cycles. Microsatellite analyses showed contemporary genetic contact between the two sympatric species and 3.0% of the samples studied were identified as hybrids. In conclusion, we found a historical and asymmetric genetic introgression (through mitochondrial capture) of P. fulvus into P. gymnonotus in Mexico and Central America and a limited contemporary gene exchange between the two species. However, no relationship was found between hybridization and the intermediate-sized specimens from southern Mexico, which might likely result from a clinal variation with latitude. These results confirm the need for caution when using forearm size to identify these species in the field and when differentiating them in the laboratory based on mitochondrial DNA alone.
Collapse
|
5
|
Wells HL, Letko M, Lasso G, Ssebide B, Nziza J, Byarugaba DK, Navarrete-Macias I, Liang E, Cranfield M, Han BA, Tingley MW, Diuk-Wasser M, Goldstein T, Johnson CK, Mazet JAK, Chandran K, Munster VJ, Gilardi K, Anthony SJ. The evolutionary history of ACE2 usage within the coronavirus subgenus Sarbecovirus. Virus Evol 2021; 7:veab007. [PMID: 33754082 PMCID: PMC7928622 DOI: 10.1093/ve/veab007] [Citation(s) in RCA: 43] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
Severe acute respiratory syndrome coronavirus 1 (SARS-CoV-1) and SARS-CoV-2 are not phylogenetically closely related; however, both use the angiotensin-converting enzyme 2 (ACE2) receptor in humans for cell entry. This is not a universal sarbecovirus trait; for example, many known sarbecoviruses related to SARS-CoV-1 have two deletions in the receptor binding domain of the spike protein that render them incapable of using human ACE2. Here, we report three sequences of a novel sarbecovirus from Rwanda and Uganda that are phylogenetically intermediate to SARS-CoV-1 and SARS-CoV-2 and demonstrate via in vitro studies that they are also unable to utilize human ACE2. Furthermore, we show that the observed pattern of ACE2 usage among sarbecoviruses is best explained by recombination not of SARS-CoV-2, but of SARS-CoV-1 and its relatives. We show that the lineage that includes SARS-CoV-2 is most likely the ancestral ACE2-using lineage, and that recombination with at least one virus from this group conferred ACE2 usage to the lineage including SARS-CoV-1 at some time in the past. We argue that alternative scenarios such as convergent evolution are much less parsimonious; we show that biogeography and patterns of host tropism support the plausibility of a recombination scenario, and we propose a competitive release hypothesis to explain how this recombination event could have occurred and why it is evolutionarily advantageous. The findings provide important insights into the natural history of ACE2 usage for both SARS-CoV-1 and SARS-CoV-2 and a greater understanding of the evolutionary mechanisms that shape zoonotic potential of coronaviruses. This study also underscores the need for increased surveillance for sarbecoviruses in southwestern China, where most ACE2-using viruses have been found to date, as well as other regions such as Africa, where these viruses have only recently been discovered.
Collapse
Affiliation(s)
- H L Wells
- Department of Ecology, Evolution, and Environmental Biology, Columbia University, 1200 Amsterdam Ave, New York, NY 10027, USA
| | - M Letko
- Laboratory of Virology, Division of Intramural Research, National Institute of Allergy and Infectious Diseases, National Institutes of Health, 903 S. 4th St, Hamilton, MT 59840, USA.,Paul G. Allen School for Global Animal Health, Washington State University, 1155 College Ave, Pullman, WA 99164, USA
| | - G Lasso
- Department of Microbiology and Immunology, Albert Einstein College of Medicine, 1300 Morris Park Ave, Bronx, NY 10462, USA
| | - B Ssebide
- Gorilla Doctors, c/o MGVP, Inc., 1089 Veterinary Medicine Drive, Davis, CA 95616, USA
| | - J Nziza
- Gorilla Doctors, c/o MGVP, Inc., 1089 Veterinary Medicine Drive, Davis, CA 95616, USA
| | - D K Byarugaba
- Makerere University Walter Reed Project, Plot 42, Nakasero Road, Kampala, Uganda.,Makerere University, College of Veterinary Medicine, Living Stone Road, Kampala, Uganda
| | - I Navarrete-Macias
- Center for Infection and Immunity, Mailman School of Public Health, Columbia University, 722 W 168th St, New York, NY 10032, USA
| | - E Liang
- Center for Infection and Immunity, Mailman School of Public Health, Columbia University, 722 W 168th St, New York, NY 10032, USA
| | - M Cranfield
- One Health Institute and Karen C. Drayer Wildlife Health Center, School of Veterinary Medicine, University of California Davis, 1089 Veterinary Medicine Drive, Davis, CA 95616, USA.,Department of Microbiology and Immunology, University of North Carolina School of Medicine, 125 Mason Farm Road, Chapel Hill, NC 27599, USA
| | - B A Han
- Cary Institute of Ecosystem Studies, 2801 Sharon Turnpike, Millbrook, NY 12545, USA
| | - M W Tingley
- Department of Ecology and Evolutionary Biology, University of California Los Angeles, 612 Charles E. Young Drive South, Los Angeles, CA 90095, USA
| | - M Diuk-Wasser
- Department of Ecology, Evolution, and Environmental Biology, Columbia University, 1200 Amsterdam Ave, New York, NY 10027, USA
| | - T Goldstein
- One Health Institute and Karen C. Drayer Wildlife Health Center, School of Veterinary Medicine, University of California Davis, 1089 Veterinary Medicine Drive, Davis, CA 95616, USA
| | - C K Johnson
- One Health Institute and Karen C. Drayer Wildlife Health Center, School of Veterinary Medicine, University of California Davis, 1089 Veterinary Medicine Drive, Davis, CA 95616, USA
| | - J A K Mazet
- One Health Institute and Karen C. Drayer Wildlife Health Center, School of Veterinary Medicine, University of California Davis, 1089 Veterinary Medicine Drive, Davis, CA 95616, USA
| | - K Chandran
- Department of Microbiology and Immunology, Albert Einstein College of Medicine, 1300 Morris Park Ave, Bronx, NY 10462, USA
| | - V J Munster
- Paul G. Allen School for Global Animal Health, Washington State University, 1155 College Ave, Pullman, WA 99164, USA
| | - K Gilardi
- Makerere University Walter Reed Project, Plot 42, Nakasero Road, Kampala, Uganda.,One Health Institute and Karen C. Drayer Wildlife Health Center, School of Veterinary Medicine, University of California Davis, 1089 Veterinary Medicine Drive, Davis, CA 95616, USA
| | - S J Anthony
- Department of Pathology, Microbiology, and Immunology, School of Veterinary Medicine, University of California Davis, One Shields Avenue, Davis, CA 95616, USA
| |
Collapse
|
6
|
Wells H, Letko M, Lasso G, Ssebide B, Nziza J, Byarugaba D, Navarrete-Macias I, Liang E, Cranfield M, Han B, Tingley M, Diuk-Wasser M, Goldstein T, Johnson C, Mazet J, Chandran K, Munster V, Gilardi K, Anthony S. The evolutionary history of ACE2 usage within the coronavirus subgenus Sarbecovirus. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2021:2020.07.07.190546. [PMID: 32676605 PMCID: PMC7359528 DOI: 10.1101/2020.07.07.190546] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
SARS-CoV-1 and SARS-CoV-2 are not phylogenetically closely related; however, both use the ACE2 receptor in humans for cell entry. This is not a universal sarbecovirus trait; for example, many known sarbecoviruses related to SARS-CoV-1 have two deletions in the receptor binding domain of the spike protein that render them incapable of using human ACE2. Here, we report three sequences of a novel sarbecovirus from Rwanda and Uganda which are phylogenetically intermediate to SARS-CoV-1 and SARS-CoV-2 and demonstrate via in vitro studies that they are also unable to utilize human ACE2. Furthermore, we show that the observed pattern of ACE2 usage among sarbecoviruses is best explained by recombination not of SARS-CoV-2, but of SARS-CoV-1 and its relatives. We show that the lineage that includes SARS-CoV-2 is most likely the ancestral ACE2-using lineage, and that recombination with at least one virus from this group conferred ACE2 usage to the lineage including SARS-CoV-1 at some time in the past. We argue that alternative scenarios such as convergent evolution are much less parsimonious; we show that biogeography and patterns of host tropism support the plausibility of a recombination scenario; and we propose a competitive release hypothesis to explain how this recombination event could have occurred and why it is evolutionarily advantageous. The findings provide important insights into the natural history of ACE2 usage for both SARS-CoV-1 and SARS-CoV-2, and a greater understanding of the evolutionary mechanisms that shape zoonotic potential of coronaviruses. This study also underscores the need for increased surveillance for sarbecoviruses in southwestern China, where most ACE2-using viruses have been found to date, as well as other regions such as Africa, where these viruses have only recently been discovered.
Collapse
Affiliation(s)
- H.L Wells
- Department of Ecology, Evolution, and Environmental Biology, Columbia University, New York, NY, USA
| | - M Letko
- Laboratory of Virology, Division of Intramural Research, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Hamilton, MT, USA
- Paul G. Allen School for Global Animal Health, Washington State University, Pullman, WA, USA
| | - G Lasso
- Department of Microbiology and Immunology, Albert Einstein College of Medicine, New York, NY 10461, USA
| | - B Ssebide
- Gorilla Doctors, c/o MGVP, Inc., Davis, California, USA
| | - J Nziza
- Gorilla Doctors, c/o MGVP, Inc., Davis, California, USA
| | - D.K Byarugaba
- Makerere University Walter Reed Project, Kampala, Uganda
- Makerere University, College of Veterinary Medicine, Kampala, Uganda
| | - I Navarrete-Macias
- Center for Infection and Immunity, Mailman School of Public Health, Columbia University, New York, NY, USA
| | - E Liang
- Center for Infection and Immunity, Mailman School of Public Health, Columbia University, New York, NY, USA
| | - M Cranfield
- One Health Institute and Karen C. Drayer Wildlife Health Center, School of Veterinary Medicine, University of California Davis, California, USA
- Department of Microbiology and Immunology, University of North Carolina School of Medicine, Chapel Hill, North Carolina, USA
| | - B.A Han
- Cary Institute of Ecosystem Studies, Millbrook, New York, USA
| | - M.W Tingley
- Department of Ecology and Evolutionary Biology, University of California Los Angeles, Los Angeles, CA, USA
| | - M Diuk-Wasser
- Laboratory of Virology, Division of Intramural Research, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Hamilton, MT, USA
| | - T Goldstein
- One Health Institute and Karen C. Drayer Wildlife Health Center, School of Veterinary Medicine, University of California Davis, California, USA
| | - C.K Johnson
- One Health Institute and Karen C. Drayer Wildlife Health Center, School of Veterinary Medicine, University of California Davis, California, USA
| | - J Mazet
- One Health Institute and Karen C. Drayer Wildlife Health Center, School of Veterinary Medicine, University of California Davis, California, USA
| | - K Chandran
- Gorilla Doctors, c/o MGVP, Inc., Davis, California, USA
| | - V.J Munster
- Paul G. Allen School for Global Animal Health, Washington State University, Pullman, WA, USA
| | - K Gilardi
- Makerere University Walter Reed Project, Kampala, Uganda
- One Health Institute and Karen C. Drayer Wildlife Health Center, School of Veterinary Medicine, University of California Davis, California, USA
| | - S.J Anthony
- Department of Ecology, Evolution, and Environmental Biology, Columbia University, New York, NY, USA
- Department of Pathology, Microbiology, and Immunology, School of Veterinary Medicine, University of California Davis, California, USA
| |
Collapse
|
7
|
Ma L, Sun H, Mao X. Transcriptome sequencing of cochleae from constant-frequency and frequency-modulated echolocating bats. Sci Data 2020; 7:341. [PMID: 33051458 PMCID: PMC7554033 DOI: 10.1038/s41597-020-00686-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2020] [Accepted: 09/14/2020] [Indexed: 02/06/2023] Open
Abstract
Echolocating bats are fascinating for their ability to ‘see’ the world in the darkness. Ultrahigh frequency hearing is essential for echolocation. In this study we collected cochlear tissues from constant-frequency (CF) bats (two subspecies of Rhinolophus affinis, Rhinolophidae) and frequency-modulated (FM) bats (Myotis ricketti, Vespertilionidae) and applied PacBio single-molecule real-time isoform sequencing (Iso-seq) technology to generate the full-length (FL) transcriptomes for the three taxa. In total of 10103, 9676 and 10504 non-redundant FL transcripts for R. a. hainanus, R. a. himalayanus and Myotis ricketti were obtained respectively. These data present a comprehensive list of transcripts involved in ultrahigh frequency hearing of echolocating bats including 26342 FL transcripts, 24833 of which are annotated by public databases. No further comparative analyses were performed on the current data in this study. This data can be reused to quantify gene or transcript expression, assess the level of alternative splicing, identify novel transcripts and improve genome annotation of bat species. Measurement(s) | cochlea • transcriptome • sequence feature annotation | Technology Type(s) | isoform sequencing • sequence annotation | Factor Type(s) | species | Sample Characteristic - Organism | Rhinolophus affinis • Myotis ricketti |
Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12999614
Collapse
Affiliation(s)
- Lu Ma
- School of Ecological and Environmental Sciences, East China Normal University, Shanghai, 200062, China.,Changsha Central Hospital, University of South China, Changsha, 410011, China
| | - Haijian Sun
- School of Ecological and Environmental Sciences, East China Normal University, Shanghai, 200062, China
| | - Xiuguang Mao
- School of Ecological and Environmental Sciences, East China Normal University, Shanghai, 200062, China. .,Institute of Eco-Chongming (IEC), East China Normal University, Shanghai, 200062, China.
| |
Collapse
|
8
|
Sun H, Chen W, Wang J, Zhang L, Rossiter SJ, Mao X. Echolocation call frequency variation in horseshoe bats: molecular basis revealed by comparative transcriptomics. Proc Biol Sci 2020; 287:20200875. [PMID: 32900318 DOI: 10.1098/rspb.2020.0875] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023] Open
Abstract
Recently diverged taxa with contrasting phenotypes offer opportunities for unravelling the genetic basis of phenotypic variation in nature. Horseshoe bats are a speciose group that exhibit a derived form of high-duty cycle echolocation in which the inner ear is finely tuned to echoes of the narrowband call frequency. Here, by focusing on three recently diverged subspecies of the intermediate horseshoe bat (Rhinolophus affinis) that display divergent echolocation call frequencies, we aim to identify candidate loci putatively involved in hearing frequency variation. We used de novo transcriptome sequencing of two mainland taxa (himalayanus and macrurus) and one island taxon (hainanus) to compare expression profiles of thousands of genes. By comparing taxa with divergent call frequencies (around 15 kHz difference), we identified 252 differentially expressed genes, of which six have been shown to be involved in hearing or deafness in human/mouse. To obtain further validation of these results, we applied quantitative reverse transcription-PCR to the candidate gene FBXL15 and found a broad association between the level of expression and call frequency across taxa. The genes identified here represent strong candidate loci associated with hearing frequency variation in bats.
Collapse
Affiliation(s)
- Haijian Sun
- School of Ecological and Environmental Sciences, East China Normal University, Shanghai 200062, People's Republic of China
| | - Wenli Chen
- School of Ecological and Environmental Sciences, East China Normal University, Shanghai 200062, People's Republic of China
| | - Jiaying Wang
- Institute of Estuarine and Coastal Research, East China Normal University, Shanghai 200062, People's Republic of China
| | - Libiao Zhang
- Guangdong Key Laboratory of Animal Conservation and Resource Utilization, Guangdong Public Laboratory of Wild Animal Conservation and Utilization, Guangdong Institute of Applied Biological Resources, Guangzhou 510260, People's Republic of China
| | - Stephen J Rossiter
- School of Biological and Chemical Sciences, Queen Mary University of London, London E1 4NS, UK
| | - Xiuguang Mao
- School of Ecological and Environmental Sciences, East China Normal University, Shanghai 200062, People's Republic of China.,Institute of Eco-Chongming (IEC), East China Normal University, Shanghai 200062, People's Republic of China
| |
Collapse
|
9
|
Mao X, Rossiter SJ. Genome-wide data reveal discordant mitonuclear introgression in the intermediate horseshoe bat (Rhinolophus affinis). Mol Phylogenet Evol 2020; 150:106886. [PMID: 32534185 DOI: 10.1016/j.ympev.2020.106886] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2019] [Revised: 06/01/2020] [Accepted: 06/03/2020] [Indexed: 10/24/2022]
Abstract
Closely related taxa often exhibit mitonuclear discordance attributed to introgression of mitochondrial DNA (mtDNA), yet few studies have considered the underlying causes of mtDNA introgression. Here we test for demographic versus adaptive processes as explanations for mtDNA introgression in three subspecies of the intermediate horseshoe bat (Rhinolophus affinis). We generated sequences of 1692 nuclear genes and 13 mitochondrial protein-coding genes for 48 individuals. Phylogenetic reconstructions based on 320 exon sequences and 2217 single nucleotide polymorphisms (SNPs) both revealed conflicts between the species tree and mtDNA tree. These results, together with geographic patterns of mitonuclear discordance, and shared identical or near-identical mtDNA sequences, suggest extensive introgression of mtDNA between the two parapatric mainland subspecies. Under demographic hypotheses, we would also expect to uncover traces of ncDNA introgression, however, population structure and gene flow analyses revealed little nuclear admixture. Furthermore, we found inconsistent estimates of the timing of population expansion and that of the most recent common ancestor for the clade containing introgressed haplotypes. Without a clear demographic explanation, we also examined whether introgression likely arises from adaptation. We found that five mtDNA genes contained fixed amino acid differences between introgressed and non-introgressed individuals, including putative positive selection found in one codon, although this did not show introgression. While our evidence for rejecting demographic hypotheses is arguably stronger than that for rejecting adaptation, we find no definitive support for either explanation. Future efforts will focus on larger-scale resequencing to decipher the underlying causes of discordant mitonuclear introgression in this system.
Collapse
Affiliation(s)
- Xiuguang Mao
- School of Ecological and Environmental Sciences, East China Normal University, Shanghai 200062, China; Institute of Eco-Chongming (IEC), East China Normal University, Shanghai 200062, China.
| | - Stephen J Rossiter
- School of Biological and Chemical Sciences, Queen Mary University of London, London E1 4NS, UK.
| |
Collapse
|
10
|
Mao X, Tsagkogeorga G, Thong VD, Rossiter SJ. Resolving evolutionary relationships among six closely related taxa of the horseshoe bats (Rhinolophus) with targeted resequencing data. Mol Phylogenet Evol 2019; 139:106551. [DOI: 10.1016/j.ympev.2019.106551] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2019] [Revised: 06/11/2019] [Accepted: 07/01/2019] [Indexed: 01/04/2023]
|
11
|
Liu T, Sun K, Csorba G, Zhang K, Zhang L, Zhao H, Jin L, Thong VD, Xiao Y, Feng J. Species delimitation and evolutionary reconstruction within an integrative taxonomic framework: A case study on Rhinolophus macrotis complex (Chiroptera: Rhinolophidae). Mol Phylogenet Evol 2019; 139:106544. [PMID: 31252069 DOI: 10.1016/j.ympev.2019.106544] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2018] [Revised: 06/21/2019] [Accepted: 06/24/2019] [Indexed: 10/26/2022]
Abstract
Species delimitation and evolutionary reconstruction remain challenging for non-model species that have experienced reticulate evolution and exhibit conflicting patterns of differentiation among multiple lines of evidence, such as mitochondrial and nuclear data and phenotypes. Here, we applied an integrative taxonomic approach to a case study of Rhinolophus macrotis complex, whose taxonomic status remains controversial, to provide insight into the systematics and evolutionary history of these species. By integrating traditional genetic markers with different modes of inheritance, genome-wide SNPs as well as phenotypic characteristics, we clarified the presence of three closely related species, R. episcopus, R. siamensis, and R. osgoodi, within this complex, and proposed a new taxonomic treatment for R. osgoodi. Our results suggested that hybridization and introgression are the main causes of low mtDNA divergence in these species. Combined with the demographic inference, we deduced that glacial-interglacial cycles drove geographic isolation and secondary contacts of these species, then promoted hybridization and lineage fusion among them, finally resulting in a reticulate evolutionary pattern. Overall, our study highlights the importance of combining multiple types of data to delimit species, especially those with conserved morphology, and to reveal the sophisticated processes of speciation.
Collapse
Affiliation(s)
- Tong Liu
- Jilin Provincial Key Laboratory of Animal Resource Conservation and Utilization, Northeast Normal University, Changchun, China
| | - Keping Sun
- Jilin Provincial Key Laboratory of Animal Resource Conservation and Utilization, Northeast Normal University, Changchun, China; Key Laboratory of Vegetation Ecology, Ministry of Education, Changchun, China.
| | - Gábor Csorba
- Department of Zoology, Hungarian Natural History Museum, Budapest, Hungary
| | - Kangkang Zhang
- Jilin Provincial Key Laboratory of Animal Resource Conservation and Utilization, Northeast Normal University, Changchun, China
| | - Lin Zhang
- Jilin Provincial Key Laboratory of Animal Resource Conservation and Utilization, Northeast Normal University, Changchun, China
| | - Hanbo Zhao
- Jilin Provincial Key Laboratory of Animal Resource Conservation and Utilization, Northeast Normal University, Changchun, China
| | - Longru Jin
- Jilin Provincial Key Laboratory of Animal Resource Conservation and Utilization, Northeast Normal University, Changchun, China
| | - Vu Dinh Thong
- Graduate University of Science and Technology, Vietnam Academy of Science and Technology, Hanoi, Viet Nam; Institute of Ecology and Biological Resources, Vietnam Academy of Science and Technology, Hanoi, Viet Nam
| | - Yanhong Xiao
- Jilin Provincial Key Laboratory of Animal Resource Conservation and Utilization, Northeast Normal University, Changchun, China
| | - Jiang Feng
- Jilin Provincial Key Laboratory of Animal Resource Conservation and Utilization, Northeast Normal University, Changchun, China; College of Life Science, Jilin Agricultural University, Changchun, China.
| |
Collapse
|
12
|
Nieto-Lugilde M, Werner O, McDaniel SF, Koutecký P, Kučera J, Rizk SM, Ros RM. Peripatric speciation associated with genome expansion and female-biased sex ratios in the moss genus Ceratodon. AMERICAN JOURNAL OF BOTANY 2018; 105:1009-1020. [PMID: 29957852 DOI: 10.1002/ajb2.1107] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/21/2017] [Accepted: 03/21/2018] [Indexed: 06/08/2023]
Abstract
PREMISE OF THE STUDY A period of allopatry is widely believed to be essential for the evolution of reproductive isolation. However, strict allopatry may be difficult to achieve in some cosmopolitan, spore-dispersed groups, like mosses. We examined the genetic and genome size diversity in Mediterranean populations of the moss Ceratodon purpureus s.l. to evaluate the role of allopatry and ploidy change in population divergence. METHODS We sampled populations of the genus Ceratodon from mountainous areas and lowlands of the Mediterranean region, and from Western and Central Europe. We performed phylogenetic and coalescent analyses on sequences from five nuclear introns and a chloroplast locus to reconstruct their evolutionary history. We also estimated genome size using flow cytometry (employing propidium iodide) and determined the sex of samples using a sex-linked PCR marker. KEY RESULTS Two well-differentiated clades were resolved, discriminating two homogeneous groups: the widespread C. purpureus and a local group mostly restricted to the mountains in Southern Spain. The latter also possessed a genome size 25% larger than the widespread C. purpureus, and the samples of this group consist entirely of females. We also found hybrids, and some of them had a genome size equivalent to the sum of the C. purpureus and Spanish genome, suggesting that they arose by allopolyploidy. CONCLUSIONS These data suggest that a new species of Ceratodon arose via peripatric speciation, potentially involving a genome size change and a strong female-biased sex ratio. The new species has hybridized in the past with C. purpureus.
Collapse
Affiliation(s)
- Marta Nieto-Lugilde
- Departamento de Biología Vegetal, Facultad de Biología, Universidad de Murcia, Campus de Espinardo, 30100, Murcia, Spain
| | - Olaf Werner
- Departamento de Biología Vegetal, Facultad de Biología, Universidad de Murcia, Campus de Espinardo, 30100, Murcia, Spain
| | - Stuart F McDaniel
- Biology Department, University of Florida, Gainesville, Florida, 32611, USA
| | - Petr Koutecký
- Faculty of Science, University of South Bohemia, Branišovská 1760, CZ-370 05, České Budějovice, Czech Republic
| | - Jan Kučera
- Faculty of Science, University of South Bohemia, Branišovská 1760, CZ-370 05, České Budějovice, Czech Republic
| | - Samah Mohamed Rizk
- Genetics Department, Faculty of Agriculture, Ain Shams University, 68 Hadayek Shubra, 11241, Cairo, Egypt
| | - Rosa M Ros
- Departamento de Biología Vegetal, Facultad de Biología, Universidad de Murcia, Campus de Espinardo, 30100, Murcia, Spain
| |
Collapse
|
13
|
Mao X, Zhang S, Rossiter SJ. Differential introgression suggests candidate beneficial and barrier loci between two parapatric subspecies of Pearson's horseshoe bat Rhinolophus pearsoni. Curr Zool 2016; 62:405-412. [PMID: 29491929 PMCID: PMC5829442 DOI: 10.1093/cz/zow017] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2015] [Accepted: 10/05/2015] [Indexed: 11/13/2022] Open
Abstract
Observations that rates of introgression between taxa can vary across loci are
increasingly common. Here, we test for differential locus-wise introgression in 2
parapatric subspecies of Pearson′s horseshoe bat (Rhinolophus pearsoni
chinensis and R. p. pearsoni). To efficiently identify
putative speciation genes and/or beneficial genes in our current system, we used a
candidate gene approach by including loci from X chromosome that are suggested to be more
likely involved in reproductive isolation in other organisms and loci underlying hearing
that have been suggested to spread across the hybrid zone in another congeneric species.
Phylogenetic and coalescent analyses were performed at 2 X-linked, 4 hearing genes, as
well as 2 other autosomal loci individually. Likelihood ratio tests could not reject the
model of zero gene flow at 2 X-linked and 2 autosomal genes. In contrast, gene flow was
supported at 3 of 4 hearing genes. While this introgression could be adaptive, we cannot
rule out stochastic processes. Our results highlight the utility of the candidate gene
approach in searching for speciation genes and/or beneficial genes across the species
boundary in natural populations.
Collapse
Affiliation(s)
- Xiuguang Mao
- Institute of Estuarine and Coastal Research, East China Normal University, Shanghai 200062, China and.,School of Biological and Chemical Sciences, Queen Mary University of London, London E1 4NS, UK
| | - Shuyi Zhang
- Institute of Estuarine and Coastal Research, East China Normal University, Shanghai 200062, China and
| | - Stephen J Rossiter
- School of Biological and Chemical Sciences, Queen Mary University of London, London E1 4NS, UK
| |
Collapse
|
14
|
Dool SE, Puechmaille SJ, Foley NM, Allegrini B, Bastian A, Mutumi GL, Maluleke TG, Odendaal LJ, Teeling EC, Jacobs DS. Nuclear introns outperform mitochondrial DNA in inter-specific phylogenetic reconstruction: Lessons from horseshoe bats (Rhinolophidae: Chiroptera). Mol Phylogenet Evol 2016; 97:196-212. [PMID: 26826601 DOI: 10.1016/j.ympev.2016.01.003] [Citation(s) in RCA: 61] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2015] [Revised: 01/07/2016] [Accepted: 01/08/2016] [Indexed: 01/22/2023]
Abstract
Despite many studies illustrating the perils of utilising mitochondrial DNA in phylogenetic studies, it remains one of the most widely used genetic markers for this purpose. Over the last decade, nuclear introns have been proposed as alternative markers for phylogenetic reconstruction. However, the resolution capabilities of mtDNA and nuclear introns have rarely been quantified and compared. In the current study we generated a novel ∼5kb dataset comprising six nuclear introns and a mtDNA fragment. We assessed the relative resolution capabilities of the six intronic fragments with respect to each other, when used in various combinations together, and when compared to the traditionally used mtDNA. We focused on a major clade in the horseshoe bat family (Afro-Palaearctic clade; Rhinolophidae) as our case study. This old, widely distributed and speciose group contains a high level of conserved morphology. This morphological stasis renders the reconstruction of the phylogeny of this group with traditional morphological characters complex. We sampled multiple individuals per species to represent their geographic distributions as best as possible (122 individuals, 24 species, 68 localities). We reconstructed the species phylogeny using several complementary methods (partitioned Maximum Likelihood and Bayesian and Bayesian multispecies-coalescent) and made inferences based on consensus across these methods. We computed pairwise comparisons based on Robinson-Foulds tree distance metric between all Bayesian topologies generated (27,000) for every gene(s) and visualised the tree space using multidimensional scaling (MDS) plots. Using our supported species phylogeny we estimated the ancestral state of key traits of interest within this group, e.g. echolocation peak frequency which has been implicated in speciation. Our results revealed many potential cryptic species within this group, even in taxa where this was not suspected a priori and also found evidence for mtDNA introgression. We demonstrated that by using just two introns one can recover a better supported species tree than when using the mtDNA alone, despite the shorter overall length of the combined introns. Additionally, when combining any single intron with mtDNA, we showed that the result is highly similar to the mtDNA gene tree and far from the true species tree and therefore this approach should be avoided. We caution against the indiscriminate use of mtDNA in phylogenetic studies and advocate for pilot studies to select nuclear introns. The selection of marker type and number is a crucial step that is best based on critical examination of preliminary or previously published data. Based on our findings and previous publications, we recommend the following markers to recover phylogenetic relationships between recently diverged taxa (<20 My) in bats and other mammals: ACOX2, COPS7A, BGN, ROGDI and STAT5A.
Collapse
Affiliation(s)
- Serena E Dool
- Department of Biological Sciences, Animal Evolution and Systematics Group, University of Cape Town, Cape Town, South Africa; Zoological Institute and Museum, University of Greifswald, Soldmann-Straße 14, D-17487 Greifswald, Germany.
| | - Sebastien J Puechmaille
- Zoological Institute and Museum, University of Greifswald, Soldmann-Straße 14, D-17487 Greifswald, Germany; Midi-Pyrénées bat group (CREN-GCMP), Toulouse, France; School of Biology and Environmental Science, University College Dublin, Belfield, Dublin 4, Ireland.
| | - Nicole M Foley
- School of Biology and Environmental Science, University College Dublin, Belfield, Dublin 4, Ireland.
| | | | - Anna Bastian
- Department of Biological Sciences, Animal Evolution and Systematics Group, University of Cape Town, Cape Town, South Africa.
| | - Gregory L Mutumi
- Department of Biological Sciences, Animal Evolution and Systematics Group, University of Cape Town, Cape Town, South Africa.
| | - Tinyiko G Maluleke
- Department of Biological Sciences, Animal Evolution and Systematics Group, University of Cape Town, Cape Town, South Africa.
| | - Lizelle J Odendaal
- Department of Biological Sciences, Animal Evolution and Systematics Group, University of Cape Town, Cape Town, South Africa.
| | - Emma C Teeling
- School of Biology and Environmental Science, University College Dublin, Belfield, Dublin 4, Ireland.
| | - David S Jacobs
- Department of Biological Sciences, Animal Evolution and Systematics Group, University of Cape Town, Cape Town, South Africa.
| |
Collapse
|