1
|
Zou D, Huang S, Tian S, Kilunda FK, Murphy RW, Dahn HA, Zhou Y, Lee PS, Chen JM. Comparative genomics sheds new light on the convergent evolution of infrared vision in snakes. Proc Biol Sci 2024; 291:20240818. [PMID: 39043244 PMCID: PMC11265913 DOI: 10.1098/rspb.2024.0818] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2024] [Revised: 05/30/2024] [Accepted: 06/19/2024] [Indexed: 07/25/2024] Open
Abstract
Infrared vision is a highly specialized sensory system that evolved independently in three clades of snakes. Apparently, convergent evolution occurred in the transient receptor potential ankyrin 1 (TRPA1) proteins of infrared-sensing snakes. However, this gene can only explain how infrared signals are received, and not the transduction and processing of those signals. We sequenced the genome of Xenopeltis unicolor, a key outgroup species of pythons, and performed a genome-wide analysis of convergence between two clades of infrared-sensing snakes. Our results revealed pervasive molecular adaptation in pathways associated with neural development and other functions, with parallel selection on loci associated with trigeminal nerve structural organization. In addition, we found evidence of convergent amino acid substitutions in a set of genes, including TRPA1 and TRPM2. The analysis also identified convergent accelerated evolution in non-coding elements near 12 genes involved in facial nerve structural organization and optic nerve development. Thus, convergent evolution occurred across multiple dimensions of infrared vision in vipers and pythons, as well as amino acid substitutions, non-coding elements, genes and functions. These changes enabled independent groups of snakes to develop and use infrared vision.
Collapse
|
2
|
Lin G, Huang Z, Yue T, Chai J, Li Y, Yang H, Qin W, Yang G, Murphy RW, Zhang YP, Zhang Z, Zhou W, Luo J. Puzzle Hi-C: An accurate scaffolding software. PLoS One 2024; 19:e0298564. [PMID: 39008464 PMCID: PMC11249255 DOI: 10.1371/journal.pone.0298564] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2024] [Accepted: 06/25/2024] [Indexed: 07/17/2024] Open
Abstract
High-quality, chromosome-scale genomes are essential for genomic analyses. Analyses, including 3D genomics, epigenetics, and comparative genomics rely on a high-quality genome assembly, which is often accomplished with the assistance of Hi-C data. Curation of genomes reveal that current Hi-C-assisted scaffolding algorithms either generate ordering and orientation errors or fail to assemble high-quality chromosome-level scaffolds. Here, we offer the software Puzzle Hi-C, which uses Hi-C reads to accurately assign contigs or scaffolds to chromosomes. Puzzle Hi-C uses the triangle region instead of the square region to count interactions in a Hi-C heatmap. This strategy dramatically diminishes scaffolding interference caused by long-range interactions. This software also introduces a dynamic, triangle window strategy during assembly. Initially small, the window expands with interactions to produce more effective clustering. Puzzle Hi-C outperforms available scaffolding tools.
Collapse
|
3
|
Peng MS, Liu YH, Shen QK, Zhang XH, Dong J, Li JX, Zhao H, Zhang H, Zhang X, He Y, Shi H, Cui C, Ouzhuluobu, Wu TY, Liu SM, Gonggalanzi, Baimakangzhuo, Bai C, Duojizhuoma, Liu T, Dai SS, Murphy RW, Qi XB, Dong G, Su B, Zhang YP. Genetic and cultural adaptations underlie the establishment of dairy pastoralism in the Tibetan Plateau. BMC Biol 2023; 21:208. [PMID: 37798721 PMCID: PMC10557253 DOI: 10.1186/s12915-023-01707-x] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2023] [Accepted: 09/20/2023] [Indexed: 10/07/2023] Open
Abstract
BACKGROUND Domestication and introduction of dairy animals facilitated the permanent human occupation of the Tibetan Plateau. Yet the history of dairy pastoralism in the Tibetan Plateau remains poorly understood. Little is known how Tibetans adapted to milk and dairy products. RESULTS We integrated archeological evidence and genetic analysis to show the picture that the dairy ruminants, together with dogs, were introduced from West Eurasia into the Tibetan Plateau since ~ 3600 years ago. The genetic admixture between the exotic and indigenous dogs enriched the candidate lactase persistence (LP) allele 10974A > G of West Eurasian origin in Tibetan dogs. In vitro experiments demonstrate that - 13838G > A functions as a LP allele in Tibetans. Unlike multiple LP alleles presenting selective signatures in West Eurasians and South Asians, the de novo origin of Tibetan-specific LP allele - 13838G > A with low frequency (~ 6-7%) and absence of selection corresponds - 13910C > T in pastoralists across eastern Eurasia steppe. CONCLUSIONS Results depict a novel scenario of genetic and cultural adaptations to diet and expand current understanding of the establishment of dairy pastoralism in the Tibetan Plateau.
Collapse
|
4
|
Nguyen SN, Nguyen VDH, LE MVAN, Nguyen LT, Vo TDH, Vo BD, Che J, Murphy RW. A new snake of the genus Dendrelaphis Boulenger, 1890 (Squamata: Colubridae) from the coastal area of southern Vietnam. Zootaxa 2023; 5318:130-144. [PMID: 37518394 DOI: 10.11646/zootaxa.5318.1.6] [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/17/2023] [Indexed: 08/01/2023]
Abstract
We describe a new species of Dendrelaphis from the coastal area in southern Vietnam based on morphological data and nucleotide sequences from one partial mitochondrial gene (Cytb). Dendrelaphis binhi sp. nov. is diagnosed by the following morphological characters: body scale rows 13 at neck and midbody, 9-11 before vent; vertebral scale row feebly enlarged; ventrals 154-158 in males and 161-170 in females; subcaudals 95-106 in males and 95-102 in females; one supralabial touching the eye; 18 or 19 maxillary teeth; hemipenis spinose, reaching 10th or 11th SC, with a slender papilla; cloacal plate divided; dorsum dark brown; and a white stripe along the lower flank. The new species differs from its congeners by an uncorrected p-distance in Cytb sequences of at least 8.1%.
Collapse
|
5
|
Zhong H, Luo J, Tang L, Liao S, Lu Z, Lin G, Murphy RW, Liu L. Association filtering and generative adversarial networks for predicting lncRNA-associated disease. BMC Bioinformatics 2023; 24:234. [PMID: 37277721 DOI: 10.1186/s12859-023-05368-z] [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/06/2023] [Accepted: 05/29/2023] [Indexed: 06/07/2023] Open
Abstract
BACKGROUND Long non-coding RNA (lncRNA) closely associates with numerous biological processes, and with many diseases. Therefore, lncRNA-disease association prediction helps obtain relevant biological information and understand pathogenesis, and thus better diagnose preventable diseases. RESULTS Herein, we offer the LDAF_GAN method for predicting lncRNA-associated disease based on association filtering and generative adversarial networks. Experimentation used two types of data: lncRNA-disease associated data without lncRNA sequence features, and fused lncRNA sequence features. LDAF_GAN uses a generator and discriminator, and differs from the original GAN by the addition of a filtering operation and negative sampling. Filtering allows the generator output to filter out unassociated diseases before being fed into the discriminator. Thus, the results generated by the model focuses only on lncRNAs associated with disease. Negative sampling takes a portion of disease terms with 0 from the association matrix as negative samples, which are assumed to be unassociated with lncRNA. A regular term is added to the loss function to avoid producing a vector with all values of 1, which can fool the discriminator. Thus, the model requires that generated positive samples are close to 1, and negative samples are close to 0. The model achieved a superior fitting effect; LDAF_GAN had superior performance in predicting fivefold cross-validations on the two datasets with AUC values of 0.9265 and 0.9278, respectively. In the case study, LDAF_GAN predicted disease association for six lncRNAs-H19, MALAT1, XIST, ZFAS1, UCA1, and ZEB1-AS1-and with the top ten predictions of 100%, 80%, 90%, 90%, 100%, and 90%, respectively, which were reported by previous studies. CONCLUSION LDAF_GAN efficiently predicts the potential association of existing lncRNAs and the potential association of new lncRNAs with diseases. The results of fivefold cross-validation, tenfold cross-validation, and case studies suggest that the model has great predictive potential for lncRNA-disease association prediction.
Collapse
|
6
|
Liu Y, Pu Y, Chen S, Wang X, Murphy RW, Wang X, Liao R, Tang K, Yue B, Liu S. Revalidation and expanded description of Mustela aistoodonnivalis (Mustelidae: Carnivora) based on a multigene phylogeny and morphology. Ecol Evol 2023; 13:e9944. [PMID: 37082328 PMCID: PMC10111237 DOI: 10.1002/ece3.9944] [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: 09/23/2022] [Revised: 03/02/2023] [Accepted: 03/13/2023] [Indexed: 04/22/2023] Open
Abstract
The lacked-teeth pygmy weasel, Mustela aistoodonnivalis Wu & Kao, 1991, was originally described as being from Taibai Mountain and Zhashui county, Shaanxi, China. Subsequently, it was considered a subspecies or synonym of Mustela nivalis. In a faunal survey of northwestern Sichuan, eight specimens of M. aistoodonnivalis were collected. A molecular phylogenetic analysis of one mitochondrial and six nuclear genes clustered the specimens as a distinct clade and not with M. nivalis. Morphologically, the lack of the second lower molar differentiated them from M. nivalis, and genetic distances were typical of discrete species. These analyses confirmed that M. aistoodonnivalis is an independent species in the genus Mustela.
Collapse
|
7
|
Wu YH, Hou SB, Yuan ZY, Jiang K, Huang RY, Wang K, Liu Q, Yu ZB, Zhao HP, Zhang BL, Chen JM, Wang LJ, Stuart BL, Chambers EA, Wang YF, Gao W, Zou DH, Yan F, Zhao GG, Fu ZX, Wang SN, Jiang M, Zhang L, Ren JL, Wu YY, Zhang LY, Yang DC, Jin JQ, Yin TT, Li JT, Zhao WG, Murphy RW, Huang S, Guo P, Zhang YP, Che J. DNA barcoding of Chinese snakes reveals hidden diversity and conservation needs. Mol Ecol Resour 2023. [PMID: 36924341 DOI: 10.1111/1755-0998.13784] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2022] [Revised: 02/25/2023] [Accepted: 03/10/2023] [Indexed: 03/17/2023]
Abstract
DNA barcoding has greatly facilitated studies of taxonomy, biodiversity, biological conservation, and ecology. Here, we establish a reliable DNA barcoding library for Chinese snakes, unveiling hidden diversity with implications for taxonomy, and provide a standardized tool for conservation management. Our comprehensive study includes 1638 cytochrome c oxidase subunit I (COI) sequences from Chinese snakes that correspond to 17 families, 65 genera, 228 named species (80.6% of named species) and 36 candidate species. A barcode gap analysis reveals gaps, where all nearest neighbour distances exceed maximum intraspecific distances, in 217 named species and all candidate species. Three species-delimitation methods (ABGD, sGMYC, and sPTP) recover 320 operational taxonomic units (OTUs), of which 192 OTUs correspond to named and candidate species. Twenty-eight other named species share OTUs, such as Azemiops feae and A. kharini, Gloydius halys, G. shedaoensis, and G. intermedius, and Bungarus multicinctus and B. candidus, representing inconsistencies most probably caused by imperfect taxonomy, recent and rapid speciation, weak taxonomic signal, introgressive hybridization, and/or inadequate phylogenetic signal. In contrast, 43 species and candidate species assign to two or more OTUs due to having large intraspecific distances. If most OTUs detected in this study reflect valid species, including the 36 candidate species, then 30% more species would exist than are currently recognized. Several OTU divergences associate with known biogeographic barriers, such as the Taiwan Strait. In addition to facilitating future studies, this reliable and relatively comprehensive reference database will play an important role in the future monitoring, conservation, and management of Chinese snakes.
Collapse
|
8
|
Li Y, Wang S, Zhang Z, Luo J, Lin GL, Deng WD, Guo Z, Han FM, Wang LL, Li J, Wu SF, Liu HQ, He S, Murphy RW, Zhang ZJ, Cooper DN, Wu DD, Zhang YP. Large-Scale Chromosomal Changes Lead to Genome-Level Expression Alterations, Environmental Adaptation, and Speciation in the Gayal (Bos frontalis). Mol Biol Evol 2023; 40:6980758. [PMID: 36625089 PMCID: PMC9874039 DOI: 10.1093/molbev/msad006] [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: 08/09/2022] [Revised: 12/20/2022] [Accepted: 01/05/2023] [Indexed: 01/11/2023] Open
Abstract
Determining the functional consequences of karyotypic changes is invariably challenging because evolution tends to obscure many of its own footprints, such as accumulated mutations, recombination events, and demographic perturbations. Here, we describe the assembly of a chromosome-level reference genome of the gayal (Bos frontalis) thereby revealing the structure, at base-pair-level resolution, of a telo/acrocentric-to-telo/acrocentric Robertsonian translocation (2;28) (T/A-to-T/A rob[2;28]). The absence of any reduction in the recombination rate or genetic introgression within the fusion region of gayal served to challenge the long-standing view of a role for fusion-induced meiotic dysfunction in speciation. The disproportionate increase noted in the distant interactions across pro-chr2 and pro-chr28, and the change in open-chromatin accessibility following rob(2;28), may, however, have led to the various gene expression irregularities observed in the gayal. Indeed, we found that many muscle-related genes, located synthetically on pro-chr2 and pro-chr28, exhibited significant changes in expression. This, combined with genome-scale structural variants and expression alterations in genes involved in myofibril composition, may have driven the rapid sarcomere adaptation of gayal to its rugged mountain habitat. Our findings not only suggest that large-scale chromosomal changes can lead to alterations in genome-level expression, thereby promoting both adaptation and speciation, but also illuminate novel avenues for studying the relationship between karyotype evolution and speciation.
Collapse
|
9
|
Babarinde IA, Adeola AC, Djagoun CAMS, Nneji LM, Okeyoyin AO, Niba G, Wanzie NK, Oladipo OC, Adebambo AO, Bello SF, Ng'ang'a SI, Olaniyi WA, Okoro VMO, Adedeji BE, Olatunde O, Ayoola AO, Matouke MM, Wang YY, Sanke OJ, Oseni SO, Nwani CD, Murphy RW. Population structure and evolutionary history of the greater cane rat ( Thryonomys swinderianus) from the Guinean Forests of West Africa. Front Genet 2023; 14:1041103. [PMID: 36923796 PMCID: PMC10010571 DOI: 10.3389/fgene.2023.1041103] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2022] [Accepted: 02/07/2023] [Indexed: 03/02/2023] Open
Abstract
Grasscutter (Thryonomys swinderianus) is a large-body old world rodent found in sub-Saharan Africa. The body size and the unique taste of the meat of this major crop pest have made it a target of intense hunting and a potential consideration as a micro-livestock. However, there is insufficient knowledge on the genetic diversity of its populations across African Guinean forests. Herein, we investigated the genetic diversity, population structures and evolutionary history of seven Nigerian wild grasscutter populations together with individuals from Cameroon, Republic of Benin, and Ghana, using five mitochondrial fragments, including D-loop and cytochrome b (CYTB). D-loop haplotype diversity ranged from 0.571 (± 0.149) in Republic of Benin to 0.921 (± 0.013) in Ghana. Within Nigeria, the haplotype diversity ranged from 0.659 (± 0.059) in Cross River to 0.837 (± 0.075) in Ondo subpopulation. The fixation index (FST), haplotype frequency distribution and analysis of molecular variance revealed varying levels of population structures across populations. No significant signature of population contraction was detected in the grasscutter populations. Evolutionary analyses of CYTB suggests that South African population might have diverged from other populations about 6.1 (2.6-10.18, 95% CI) MYA. Taken together, this study reveals the population status and evolutionary history of grasscutter populations in the region.
Collapse
|
10
|
Liu S, Zhou C, Meng G, Wan T, Tang M, Yang C, Murphy RW, Fan Z, Liu Y, Zeng T, Zhao Y, Liu S. Evolution and diversification of Mountain voles (Rodentia: Cricetidae). Commun Biol 2022; 5:1417. [PMID: 36572770 PMCID: PMC9792541 DOI: 10.1038/s42003-022-04371-z] [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: 03/05/2022] [Accepted: 12/13/2022] [Indexed: 12/27/2022] Open
Abstract
The systematics of the Cricetid genus Neodon have long been fraught with uncertainty due to sampling issues and a lack of comprehensive datasets. To gain better insights into the phylogeny and evolution of Neodon, we systematically sampled Neodon across the Hengduan and Himalayan Mountains, which cover most of its range in China. Analyses of skulls, teeth, and bacular structures revealed 15 distinct patterns corresponding to 15 species of Neodon. In addition to morphological analyses, we generated a high-quality reference genome for the mountain vole and generated whole-genome sequencing data for 47 samples. Phylogenomic analyses supported the recognition of six new species, revealing a long-term underestimation of Neodon diversity. We further identified positively selected genes potentially related to high-elevation adaptation. Together, our results illuminate how climate change caused the plateau to become the centre of Neodon origin and diversification and how mountain voles have adapted to the hypoxic high-altitude plateau environment.
Collapse
|
11
|
French CM, Berezin CT, Overcast I, Méndez De La Cruz FR, Basu S, Martínez Bernal RL, Murphy RW, Hickerson MJ, Blair C. Forest cover and geographical distance influence fine-scale genetic structure of leaf-toed geckos in the tropical dry forests of western Mexico. Biol J Linn Soc Lond 2022. [DOI: 10.1093/biolinnean/blac118] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
Abstract
Abstract
The biodiversity within tropical dry forests (TDFs) is astounding and yet poorly catalogued due to inadequate sampling and the presence of cryptic species. In the Mexican TDF, endemic species are common, and the landscape has been continually altered by geological and anthropogenic changes. To understand how landscape and environmental variables have shaped the population structure of endemic species, we studied the recently described species of leaf-toed gecko, Phyllodactylus benedettii, in coastal western Mexico. Using double-digest restriction site-associated DNA sequencing data, we first explore population structure and estimate the number of ancestral populations. The results indicate a high degree of genetic structure with little admixture, and patterns corresponding to both latitudinal and altitudinal gradients. We find that genetic structure cannot be explained purely by geographical distance, and that ecological corridors may facilitate dispersal and gene flow. We then model the spatial distribution of P. benedettii in the TDF through time and find that the coastline has been climatically suitable for the species since the Last Glacial Maximum. Landscape genetic analyses suggest that the combination of isolation by distance (IBD) and isolation by resistance (IBR; forest cover) has influenced the spatial genetic structure of the species. Overall, our genomic data demonstrate fine-scale population structure in TDF habitat, a complex colonization history, and spatial patterns consistent with both IBD and other ecological factors. These results further highlight the Mexican TDF as a diversity hotspot and suggest that continued anthropogenic changes are likely to affect native fauna.
Collapse
|
12
|
LaFond J, Martin KR, Dahn H, Richmond JQ, Murphy RW, Rollinson N, Savage AE. Invasive Bullfrogs Maintain MHC Polymorphism Including Alleles Associated with Chytrid Fungal Infection. Integr Comp Biol 2022; 62:262-274. [PMID: 35588059 DOI: 10.1093/icb/icac044] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2022] [Revised: 04/19/2022] [Accepted: 05/16/2022] [Indexed: 11/14/2022] Open
Abstract
Maintenance of genetic diversity at adaptive loci may facilitate invasions by non-native species by allowing populations to adapt to novel environments, despite the loss of diversity at neutral loci that typically occurs during founder events. To evaluate this prediction, we compared genetic diversity at major histocompatibility complex (MHC) and cytochrome b (cytb) loci from 20 populations of the American bullfrog (Rana catesbeiana) across theinvasive and native ranges in North America and quantified the presence of the pathogen Batrachochytrium dendrobatidis (Bd). Compared to native populations, invasive populations had significantly higher Bd prevalence and intensity, significantly higher pairwise MHC and cytb FST, and significantly lower cytb diversity, but maintained similar levels of MHC diversity. The two most common MHC alleles (LiCA_B and Rapi_33) were associated with a significant decreased risk of Bd infection, and we detected positive selection acting on four peptide binding residues. Phylogenetic analysis suggested invasive populations likely arose from a single founding population in the American Midwest with a possible subsequent invasion in the northwest. Overall, our study suggests that the maintenance of diversity at adaptive loci may contribute to invasion success and highlights the importance of quantifying diversity at functional loci to assess the evolutionary potential of invasive populations.
Collapse
|
13
|
Dahn HA, Mountcastle J, Balacco J, Winkler S, Bista I, Schmitt AD, Pettersson OV, Formenti G, Oliver K, Smith M, Tan W, Kraus A, Mac S, Komoroske LM, Lama T, Crawford AJ, Murphy RW, Brown S, Scott AF, Morin PA, Jarvis ED, Fedrigo O. Benchmarking ultra-high molecular weight DNA preservation methods for long-read and long-range sequencing. Gigascience 2022; 11:6659719. [PMID: 35946988 PMCID: PMC9364683 DOI: 10.1093/gigascience/giac068] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2021] [Revised: 01/26/2022] [Accepted: 06/16/2022] [Indexed: 11/14/2022] Open
Abstract
BACKGROUND Studies in vertebrate genomics require sampling from a broad range of tissue types, taxa, and localities. Recent advancements in long-read and long-range genome sequencing have made it possible to produce high-quality chromosome-level genome assemblies for almost any organism. However, adequate tissue preservation for the requisite ultra-high molecular weight DNA (uHMW DNA) remains a major challenge. Here we present a comparative study of preservation methods for field and laboratory tissue sampling, across vertebrate classes and different tissue types. RESULTS We find that storage temperature was the strongest predictor of uHMW fragment lengths. While immediate flash-freezing remains the sample preservation gold standard, samples preserved in 95% EtOH or 20-25% DMSO-EDTA showed little fragment length degradation when stored at 4°C for 6 hours. Samples in 95% EtOH or 20-25% DMSO-EDTA kept at 4°C for 1 week after dissection still yielded adequate amounts of uHMW DNA for most applications. Tissue type was a significant predictor of total DNA yield but not fragment length. Preservation solution had a smaller but significant influence on both fragment length and DNA yield. CONCLUSION We provide sample preservation guidelines that ensure sufficient DNA integrity and amount required for use with long-read and long-range sequencing technologies across vertebrates. Our best practices generated the uHMW DNA needed for the high-quality reference genomes for phase 1 of the Vertebrate Genomes Project, whose ultimate mission is to generate chromosome-level reference genome assemblies of all ∼70,000 extant vertebrate species.
Collapse
|
14
|
Tu F, Zhao L, Murphy RW, Zhai X, Wang J. Phylogeography of the Hong Kong Whipping Frog (<i>Polypedates megacephalus</i>) throughout Hainan Island, China. RUSS J HERPETOL 2022. [DOI: 10.30906/1026-2296-2022-29-3-137-148] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
The Hong Kong Whipping Frog (Polypedates megacephalus) is widespread on Hainan Island and yet nothing is known about its history. To investigate this, we use sequences of a partial fragment of mitochondrial ND1 and its flanking region from 248 individuals taken from 30 localities. Analyses reveal low genetic diversity and commonly shared haplotypes. Montane areas exhibit a significantly greater genetic diversity than lowland areas (P 0.01). The origin of P. megacephalus on Hainan Island dates to 1.06 Ma in the Pleistocene. Population expansion events were associated with the transition from the last interglacial period to the last glacial maximum. The matrilineal history indicates dispersal from the montane to the lowland. Bidirectional dispersion occurs between northeastern and southwestern parts of the island. Low FCT values for the northeastern and the southwestern areas associate with rivers, whereas the highest values occur among two geographical groups. Thus, mountains and rivers appear to have effected dispersal. Pleistocene climatic cooling and the creation of a landbridge with mainland China owing to lowered sea levels best explains the occurrence of P. megacephalus on Hainan Island. The analyses provide insights into the effects of Pleistocene climatic fluctuations on an amphibian on Hainan Island.
Collapse
|
15
|
Gao W, Yu CX, Zhou WW, Zhang BL, Chambers EA, Dahn HA, Jin JQ, Murphy RW, Zhang YP, Che J. Species persistence with hybridization in toad-headed lizards driven by divergent selection and low recombination. Mol Biol Evol 2022; 39:6561330. [PMID: 35356979 PMCID: PMC9007161 DOI: 10.1093/molbev/msac064] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Speciation plays a central role in evolutionary studies, and particularly how reproductive isolation (RI) evolves. The origins and persistence of RI are distinct processes that require separate evaluations. Treating them separately clarifies the drivers of speciation and then it is possible to link the processes to understand large-scale patterns of diversity. Recent genomic studies have focused predominantly on how species or RI originate. However, we know little about how species persist in face of gene flow. Here, we evaluate a contact zone of two closely related toad-headed lizards (Phrynocephalus) using a chromosome-level genome assembly and population genomics. To some extent, recent asymmetric introgression from Phrynocephalus putjatai to P. vlangalii reduces their genomic differences. However, their highly divergent regions (HDRs) have heterogeneous distributions across the genomes. Functional gene annotation indicates that many genes within HDRs are involved in reproduction and RI. Compared with allopatric populations, contact areas exhibit recent divergent selection on the HDRs and a lower population recombination rate. Taken together, this implies that divergent selection and low genetic recombination help maintain RI. This study provides insights into the genomic mechanisms that drive RI and two species persistence in the face of gene flow during the late stage of speciation.
Collapse
|
16
|
Yuan ZY, Wu YK, Yan F, Murphy RW, Papenfuss TJ, Wake DB, Zhang YP, Che J. Comparative multi-locus assessment of modern Asian newts (<i>Cynops</i>, <i>Paramesotriton</i>, and <i>Pachytriton</i>: Salamandridae) in southern China suggests a shared biogeographic history. Zool Res 2022; 43:706-718. [PMID: 35927393 PMCID: PMC9486527 DOI: 10.24272/j.issn.2095-8137.2022.080] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|
17
|
Chai J, Lu CQ, Yi MR, Dai NH, Weng XD, Di MX, Peng Y, Tang Y, Shan QH, Wang K, Liu HZ, Zhao HP, Jin JQ, Cao RJ, Lu P, Luo LC, Murphy RW, Zhang YP, Che J. Discovery of a wild, genetically pure Chinese giant salamander creates new conservation opportunities. Zool Res 2022; 43:469-480. [PMID: 35514224 PMCID: PMC9113980 DOI: 10.24272/j.issn.2095-8137.2022.101] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Effective conservation of threatened biota relies on accurate assessments and scientific guidance. As an unfortunate example, Chinese giant salamanders (Andrias, CGS) remain critically endangered in nature. Misguided conservation efforts, e.g., commercial propagation and releasing of millions of likely non-indigenous or interspecific hybrids, have further compromised conservation initiatives. Limited information on wild populations of CGS poses a significant conservation challenge. Following 18-month long field monitoring, we now report the discovery of a wild population of CGS in a closed nature reserve in Jiangxi Province, China. Genomic assessments reveal its genetic distinctiveness and do not detect genetic admixture with other species. Based on morphological and molecular evidences, we describe this CGS as a new species Andrias jiangxiensissp. nov. This is the only known species of CGS today with a genetically pure, reproducing, in situ population. This discovery emphasizes the important role that closed nature reserves play in protecting species, and the necessity of integrating long-term field monitoring and genetic assessments. It sets a new pathway for discovering and conserving endangered species, especially for those biotas that are similarly being extirpated by anthropogenic translocations and overexploitation.
Collapse
|
18
|
Smart U, Ingrasci MJ, Sarker GC, Lalremsanga H, Murphy RW, Ota H, Tu MC, Shouche Y, Orlov NL, Smith EN. A comprehensive appraisal of evolutionary diversity in venomous Asian coralsnakes of the genus
Sinomicrurus
(Serpentes: Elapidae) using Bayesian coalescent inference and supervised machine learning. J ZOOL SYST EVOL RES 2021. [DOI: 10.1111/jzs.12547] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
19
|
Zhang T, Chen J, Zhang J, Guo YT, Zhou X, Li MW, Zheng ZZ, Zhang TZ, Murphy RW, Nevo E, Shi P. Phenotypic and genomic adaptations to the extremely high elevation in plateau zokor (Myospalax baileyi). Mol Ecol 2021; 30:5765-5779. [PMID: 34510615 DOI: 10.1111/mec.16174] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2021] [Revised: 08/07/2021] [Accepted: 08/31/2021] [Indexed: 12/31/2022]
Abstract
The evolutionary outcomes of high elevation adaptation have been extensively described. However, whether widely distributed high elevation endemic animals adopt uniform mechanisms during adaptation to different elevational environments remains unknown, especially with respect to extreme high elevation environments. To explore this, we analysed the phenotypic and genomic data of seven populations of plateau zokor (Myospalax baileyi) along elevations ranging from 2,700 to 4,300 m. Based on whole-genome sequencing data and demographic reconstruction of the evolutionary history, we show that two populations of plateau zokor living at elevations exceeding 3,700 m diverged from other populations nearly 10,000 years ago. Further, phenotypic comparisons reveal stress-dependent adaptation, as two populations living at elevations exceeding 3,700 m have elevated ratios of heart mass to body mass relative to other populations, and the highest population (4,300 m) displays alterations in erythrocytes. Correspondingly, genomic analysis of selective sweeps indicates that positive selection might contribute to the observed phenotypic alterations in these two extremely high elevation populations, with the adaptive cardiovascular phenotypes of both populations possibly evolving under the functional constrains of their common ancestral population. Taken together, phenotypic and genomic evidence demonstrates that heterogeneous stressors impact adaptations to extreme elevations and reveals stress-dependent and genetically constrained adaptation to hypoxia, collectively providing new insights into the high elevation adaptation.
Collapse
|
20
|
Shen QK, Peng MS, Adeola AC, Kui L, Duan S, Miao YW, Eltayeb NM, Lichoti JK, Otecko NO, Strillacci MG, Gorla E, Bagnato A, Charles OS, Sanke OJ, Dawuda PM, Okeyoyin AO, Musina J, Njoroge P, Agwanda B, Kusza S, Nanaei HA, Pedar R, Xu MM, Du Y, Nneji LM, Murphy RW, Wang MS, Esmailizadeh A, Dong Y, Ommeh SC, Zhang YP. Erratum to: Genomic analyses unveil helmeted guinea fowl (Numida meleagris) domestication in West Africa. Genome Biol Evol 2021; 13:6355925. [PMID: 34417609 PMCID: PMC8379371 DOI: 10.1093/gbe/evab174] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
|
21
|
Yan F, Nneji LM, Jin JQ, Yuan ZY, Chen JM, Mi X, Chen HM, Murphy RW, Che J. Multi-locus genetic analyses of Quasipaa from throughout its distribution. Mol Phylogenet Evol 2021; 163:107218. [PMID: 34082130 DOI: 10.1016/j.ympev.2021.107218] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2020] [Revised: 02/24/2021] [Accepted: 05/27/2021] [Indexed: 11/16/2022]
Abstract
Montane frogs of the genus Quasipaa Dubois, 1992 occur from southern China to Southeast Asia (Frost 2021). Analyses of mtDNA (Cytb) and nuDNA data (Rag1, Rag2, Rhod, Tyr) for samples from 93 localities throughout its distribution yield a phylogeny. Clades A and B occur in Southeast Asia, clade C in northern Yangtze River, China, clade D in southwestern China, and clades E and F in southeastern China. Results place Q. yei within monophyletic Quasipaa and identify two new species. Based on nuDNA data, the basal split of clade A and B indicates an Indochinese origin of Quasipaa. The west-east diversification of five species across South China (Q. spinosa, Q. exilispinosa, Q. jiulongensis, Q. shini, Q. boulengeri) corresponds to topographic terrains II and III of China. Divergence of species from southeastern China (Q. shini, Q. jiulongensis, Q. spinosa, Q. exilispinosa) and southwestern China (Q. boulengeri) dates to 15.30-16.56 Ma (million years ago). A principal component analysis (PCA) and t-test involving 19 bioclimatic variables identifies significantly different environmental conditions between the two regions. Species' distribution models (SDM) for Q. spinosa and Q. boulengeri identify the best areas to be eastern and western South China, respectively. Thus, environmental variation appears to have influenced the genetic divergence and distributions of Quasipaa in South China. Mito-nuclear discordance indicates that some individuals of Q. exilispinosa and Q. spinosa hybridized historically.
Collapse
|
22
|
Lam NQ, Thu TTA, Nguyen LT, Murphy RW, Nguyen SN. A new species of Ahaetulla Link, 1807 (Squamata: Colubridae) from the Mekong Delta, Vietnam. Zootaxa 2021; 4966:290304. [PMID: 34186616 DOI: 10.11646/zootaxa.4966.3.2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2021] [Indexed: 11/04/2022]
Abstract
A new vine snake, genus Ahaetulla, from Soc Trang Province, southern Vietnam is described based on morphological data and nucleotide sequences from COI and Cytb. Ahaetulla rufusoculara sp. nov. is diagnosed by the following morphological characters: bright red eyes; snout without dermal appendage; internasal separated from supralabial by loreal; body scale rows 15-15-13; ventrals 186190 in males and 182185 in females; subcaudals 143153 in males and 128 or 129 in females; 8 (rarely 9) supralabials, fourth and fifth or fourth to sixth entering orbit; 1214 maxillary teeth; hemipenis short, reaching 6th or 7th SC; cloacal plate divided; dorsum bright green; and yellow or white stripe along the lower flank. The new species differs from its congeners by an uncorrected p-distance in COI and Cytb sequences of at least 7.7% and 7.5%, respectively.
Collapse
|
23
|
Nguyen VDH, Nguyen LT, Orlov NL, Murphy RW, Nguyen SN. Gekko phuyenensis sp. nov. (Squamata: Gekkonidae), a new species from southern Vietnam. Zootaxa 2021; 4966:175186. [PMID: 34186624 DOI: 10.11646/zootaxa.4966.2.5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2021] [Indexed: 11/04/2022]
Abstract
A new species of the genus Gekko Laurenti, 1768 is described from Phu Yen Province, southern Vietnam based on morphological and molecular evidence. Gekko phuyenensis sp. nov. differs from other congeners by its moderate size, with a snout-vent length of 90.292.8 mm; nares in contact with rostral; internasal 1; interorbitals 35 or 36; ventrals 3032; dorsal tubercle rows 9 or 10; precloacal pores 12 or 13 in males and precloacal pits 13 in female; subdigital lamellae under first finger 1821, fourth finger 19 or 20, first toe 1820, and fourth toe 2123; subcaudals median enlarged; no webbing between toes; and the absence of tubercles on hind limbs. Molecular phylogenetic analyses of mitochondrial gene COI (681 bp) using the geckos from Phu Yen and its congeners in southern Vietnam resolve the new species as the sister lineage to G. cf. grossmanni from Khanh Hoa Province. Genetic distances between G. phuyenensis sp. nov. and its congeners in southern Vietnam range from 16.924.2%. The discovery of this new species raises the number of known species of Gekko in Vietnam to 14.
Collapse
|
24
|
Shen QK, Peng MS, Adeola AC, Kui L, Duan S, Miao YW, Eltayeb NM, Lichoti JK, Otecko NO, Strillacci MG, Gorla E, Bagnato A, Charles OS, Sanke OJ, Dawuda PM, Okeyoyin AO, Musina J, Njoroge P, Agwanda B, Kusza S, Nanaei HA, Pedar R, Xu MM, Du Y, Nneji LM, Murphy RW, Wang MS, Esmailizadeh A, Dong Y, Ommeh SC, Zhang YP. Genomic Analyses of Unveil Helmeted Guinea Fowl (Numida meleagris) Domestication in West Africa. Genome Biol Evol 2021; 13:6261762. [PMID: 34009300 PMCID: PMC8214406 DOI: 10.1093/gbe/evab090] [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] [Accepted: 04/30/2021] [Indexed: 12/22/2022] Open
Abstract
Domestication of the helmeted guinea fowl (HGF; Numida meleagris) in Africa remains elusive. Here we report a high-quality de novo genome assembly for domestic HGF generated by long- and short-reads sequencing together with optical and chromatin interaction mapping. Using this assembly as the reference, we performed population genomic analyses for newly sequenced whole-genomes for 129 birds from Africa, Asia, and Europe, including domestic animals (n = 89), wild progenitors (n = 34), and their closely related wild species (n = 6). Our results reveal domestication of HGF in West Africa around 1,300-5,500 years ago. Scanning for selective signals characterized the functional genes in behavior and locomotion changes involved in domestication of HGF. The pleiotropy and linkage in genes affecting plumage color and fertility were revealed in the recent breeding of Italian domestic HGF. In addition to presenting a missing piece to the jigsaw puzzle of domestication in poultry, our study provides valuable genetic resources for researchers and breeders to improve production in this species.
Collapse
|
25
|
Ramírez‐Reyes T, Flores‐Villela O, Piñero D, Lathrop A, Murphy RW. Genomic assessment of the
Phyllodactylus tuberculosus
complex (Reptilia: Phyllodactylidae) in America. ZOOL SCR 2021. [DOI: 10.1111/zsc.12492] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|