Genome-wide data reveal cryptic diversity and genetic introgression in an Oriental cynopterine fruit bat radiation

Cryptic diversity and phylogeographic patterns of Mediodactylus species in the Eastern Mediterranean region

Cryptic diversity poses a great obstacle in our attempts to assess the current biodiversity crisis and may hamper conservation efforts. The gekkonid genus Mediodactylus, a well-known case of hidden species and genetic diversity, has been taxonomically reclassified several times during the last decade. Focusing on the Mediterranean populations, a recent study within the M. kotschyi species complex using classic mtDNA/nuDNA markers suggested the existence of five distinct species, some being endemic and some possibly threatened, yet their relationships have not been fully resolved. Here, we generated genome-wide SNPs (using ddRADseq) and applied molecular species delimitation approaches and population genomic analyses to further disentangle these relationships. Τhe most extensive nuclear dataset, so far, encompassing 2,360 loci and ∼ 699,000 bp from across the genome of Mediodactylus gecko, enabled us to resolve previously obscure phylogenetic relationships among the five, recently elevated, Mediodactylus species and to support the hypothesis that the taxon includes several new, undescribed species. Population genomic analyses within each of the proposed species showed strong genetic structure and high levels of genetic differentiation among populations.

When colors mislead: Genomics and bioacoustics prompt re-classification of Asian flycatcher radiation (Aves: Niltavinae)

Traditional classification of many animals, including birds, has been highly dependent on external morphological characters like plumage coloration. However, both bioacoustics and genetic or genomic data have revolutionized our understanding of the relationships of certain lineages and led to sweeping taxonomic re-organizations. In this study, we present a case of erroneous delimitation of genus boundaries in the species-rich flycatcher subfamily Niltavinae. Genera within this subfamily have historically been delineated based on blue versus brown male body plumage until recent studies based on a few mitochondrial and nuclear loci unearthed several cases of generic misclassification. Here we use extensive bioacoustic data from 43 species and genomic data from 28 species for a fundamental reclassification of species in the Niltavinae. Our study reveals that song is an important trait to classify these birds even at the genus level, whereas plumage traits exhibit ample convergence and have led to numerous historic misattributions. Our taxonomic re-organization leads to new biogeographic limits of major genera, such that the genus Cyornis now only extends as far east as the islands of Sulawesi, Sula, and Banggai, whereas Eumyias is redefined to extend far beyond Wallace's Line to the islands of Seram and Timor. Our conclusions advise against an over-reliance on morphological traits and underscore the importance of integrative datasets.

Cryptic diversity begets challenges and opportunities in biodiversity research

How many species of life are there on Earth? This is a question that we want to know but cannot yet answer. Some scholars speculate that the number of species may reach 2.2 billion when considering cryptic diversity and that each morphology-based insect species may contain an average of 3.1 cryptic species. With nearly two million described species, such high estimates of cryptic diversity would suggest that cryptic species are widespread. The development of molecular species delimitation has led to the discovery of a large number of cryptic species, and cryptic biodiversity has gradually entered our field of vision and attracted more attention. This paper introduces the concept of cryptic species, how they evolve, and methods by which they may be discovered and confirmed, and provides theoretical and methodological guidance for the study of hidden species. A workflow of how to confirm cryptic species is provided. In addition, the importance and reliability of multi-evidence-based integrated taxonomy are reaffirmed as a way to better standardize decision-making processes. Special focus on cryptic diversity and increased funding for taxonomy is needed to ensure that cryptic species in hyperdiverse groups are discoverable and described. An increased focus on cryptic species in the future will naturally arise as more difficult groups are studied, and thereby, we may finally better understand the rules governing the evolution and maintenance of cryptic biodiversity.

Next Generation Sequencing Revolutionizes Organismal Biology Research in Bats

The advent of next generation sequencing technologies (NGS) has greatly accelerated our understanding of critical aspects of organismal biology from non-model organisms. Bats form a particularly interesting group in this regard, as genomic data have helped unearth a vast spectrum of idiosyncrasies in bat genomes associated with bat biology, physiology, and evolution. Bats are important bioindicators and are keystone species to many eco-systems. They often live in proximity to humans and are frequently associated with emerging infectious diseases, including the COVID-19 pandemic. Nearly four dozen bat genomes have been published to date, ranging from drafts to chromosomal level assemblies. Genomic investigations in bats have also become critical towards our understanding of disease biology and host-pathogen coevolution. In addition to whole genome sequencing, low coverage genomic data like reduced representation libraries, resequencing data, etc. have contributed significantly towards our understanding of the evolution of natural populations, and their responses to climatic and anthropogenic perturbations. In this review, we discuss how genomic data have enhanced our understanding of physiological adaptations in bats (particularly related to ageing, immunity, diet, etc.), pathogen discovery, and host pathogen co-evolution. In comparison, the application of NGS towards population genomics, conservation, biodiversity assessment, and functional genomics has been appreciably slower. We reviewed the current areas of focus, identifying emerging topical research directions and providing a roadmap for future genomic studies in bats.

Novel de Novo Genome of Cynopterus brachyotis Reveals Evolutionarily Abrupt Shifts in Gene Family Composition across Fruit Bats

Major novel physiological or phenotypic adaptations often require accompanying modifications at the genic level. Conversely, the detection of considerable contractions and/or expansions of gene families can be an indicator of fundamental but unrecognized physiological change. We sequenced a novel fruit bat genome (Cynopterus brachyotis) and adopted a comparative approach to reconstruct the evolution of fruit bats, mapping contractions and expansions of gene families along their evolutionary history. Despite a radical change in life history as compared with other bats (e.g., loss of echolocation, large size, and frugivory), fruit bats have undergone surprisingly limited change in their genic composition, perhaps apart from a potentially novel gene family expansion relating to telomere protection and longevity. In sharp contrast, within fruit bats, the new Cynopterus genome bears the signal of unusual gene loss and gene family contraction, despite its similar morphology and lifestyle to two other major fruit bat lineages. Most missing genes are regulatory, immune-related, and olfactory in nature, illustrating the diversity of genomic strategies employed by bats to contend with responses to viral infection and olfactory requirements. Our results underscore that significant fluctuations in gene family composition are not always associated with obvious examples of novel physiological and phenotypic adaptations but may often relate to less-obvious shifts in immune strategies.

Current and Forthcoming Approaches for Benchmarking Genetic and Genomic Diversity

The current attrition of biodiversity extends beyond loss of species and unique populations to steady loss of a vast genomic diversity that remains largely undescribed. Yet the accelerating development of new techniques allows us to survey entire genomes ever faster and cheaper, to obtain robust samples from a diversity of sources including degraded DNA and residual DNA in the environment, and to address conservation efforts in new and innovative ways. Here we review recent studies that highlight the importance of carefully considering where to prioritize collection of genetic samples (e.g., organisms in rapidly changing landscapes or along edges of geographic ranges) and what samples to collect and archive (e.g., from individuals of little-known subspecies or populations, even of species not currently considered endangered). Those decisions will provide the sample infrastructure to detect the disappearance of certain genotypes or gene complexes, increases in inbreeding levels, and loss of genomic diversity as environmental conditions change. Obtaining samples from currently endangered, protected, and rare species can be particularly difficult, thus we also focus on studies that use new, non-invasive ways of obtaining genomic samples and analyzing them in these cases where other sampling options are highly constrained. Finally, biological collections archiving such samples face an inherent contradiction: their main goal is to preserve biological material in good shape so it can be used for scientific research for centuries to come, yet the technologies that can make use of such materials are advancing faster than collections can change their standardized practices. Thus, we also discuss current and potential new practices in biological collections that might bolster their usefulness for future biodiversity conservation research.

The wall lizards of the Balkan peninsula: Tackling questions at the interface of phylogenomics and population genomics

Wall lizards of the genus Podarcis (Sauria, Lacertidae) are the predominant reptile group in southern Europe, including 24 recognized species. Mitochondrial DNA data have shown that, with the exception of P. muralis, the Podarcis species distributed in the Balkan peninsula form a species group that is further sub-divided into two subgroups: the one of "P. tauricus" consisting of P. tauricus, P. milensis, P. gaigeae, and P. melisellensis, and the other of "P. erhardii" comprising P. erhardii, P. levendis, P. cretensis, and P. peloponnesiacus. In an attempt to explore the Balkan Podarcis phylogenomic relationships, assess the levels of genetic structure and to re-evaluate the number of extant species, we employed phylogenomic and admixture approaches on ddRADseq (double digested Restriction site Associated DNA sequencing) genomic data. With this efficient Next Generation Sequencing approach, we were able to obtain a large number of genomic loci randomly distributed throughout the genome and use them to resolve the previously obscure phylogenetic relationships among the different Podarcis species distributed in the Balkans. The obtained phylogenomic relationships support the monophyly of both aforementioned subgroups and revealed several divergent lineages within each subgroup, stressing the need for taxonomic re-evaluation of Podarcis' species in Balkans. The phylogenomic trees and the species delimitation analyses confirmed all recently recognized species (P. levendis, P. cretensis, and P. ionicus) and showed the presence of at least two more species, one in P. erhardii and the other in P. peloponnesiacus.

Gene Flow in Volant Vertebrates: Species Biology, Ecology and Climate Change

Gene flow, the exchange of genetic material between populations is an important biological process, which shapes and maintains biodiversity. The successful movement of individuals between populations depends on multiple factors determined by species biology and the environment. One of the most important factors regulating gene flow is the ability to move, and flight allows individuals to easily move across geographical barriers. Volant vertebrates are found on some of the remotest islands and contribute significantly to the biodiversity and ecosystem. The availability of next-generation sequencing data for non-model animals has substantially improved our understanding of gene flow and its consequences, allowing us to look at fine-scale patterns. However, most of our understanding regarding gene flow comes from the temperate regions and the Neotropics. The lack of studies from species-rich Asia is striking. In this review, we outline the importance of gene flow and the factors affecting gene flow, especially for volant vertebrates. We especially discuss research studies from tropical biomes of South and Southeast Asia, highlight the lacuna in literature and provide an outline for future studies in this species-rich region.

Population genomics of two congeneric Palaearctic shorebirds reveals differential impacts of Quaternary climate oscillations across habitats types

Intracontinental biotic divisions across the vast Palaearctic region are not well-characterized. Past research has revealed patterns ranging from a lack of population structure to deep divergences along varied lines of separation. Here we compared biogeographic patterns of two Palaearctic shorebirds with different habitat preferences, Whimbrel (Numenius phaeopus) and Eurasian curlew (N. arquata). Using genome-wide markers from populations across the Palaearctic, we applied a multitude of population genomic and phylogenomic approaches to elucidate population structure. Most importantly, we tested for isolation by distance and visualized barriers and corridors to gene flow. We found shallow Palaearctic population structure in subpolar bog and tundra-breeding whimbrels, consistent with other species breeding at a similarly high latitude, indicating connectivity across the tundra belt, both presently and during southward shifts in periods of global cooling. In contrast, the temperate grassland-breeding Eurasian curlew emerged in three distinct clades corresponding to glacial refugia. Barriers to gene flow coincided with areas of topographic relief in the central Palaearctic for whimbrels and further east for Eurasian curlews. Our findings highlight the interplay of historic and ecological factors in influencing present-day population structure of Palaearctic biota.

Conservation genomics in the fight to help the recovery of the critically endangered Siamese crocodile Crocodylus siamensis

Endangered species are often characterized by low genetic diversity and it is imperative for conservation efforts to incorporate the knowledge obtained from genetic studies for effective management. However, despite the promise of technological advances in sequencing, application of genome-wide data to endangered populations remains uncommon. In the present study we pursued a holistic conservation-genomic approach to inform a field-based management programme of a Critically Endangered species, the Siamese crocodile Crocodylus siamensis. Using thousands of single nucleotide polymorphisms from throughout the genome, we revealed signals of introgression from two other crocodile species within our sample of both wild and captive-bred Siamese crocodiles from Cambodia. Our genetic screening of the Siamese crocodiles resulted in the subsequent re-introduction of 12 individuals into the wild as well as the selection of four individuals for captive breeding programmes. Comparison of intraspecific genetic diversity revealed an alarmingly low contemporary effective population size in the wild (<50) with evidence of a recent bottleneck around Tonle Sap Lake. We also projected a probable future extinction in the wild (within fewer than five generations) in this population in the absence of re-introduction efforts. However, an increase in the number of potential breeders through re-introductions, including the one resulting from this project, could counter this trend. Our results have been implemented in ongoing re-introduction and captive breeding programmes, with major implications for the conservation management of Siamese crocodiles, and provide a blueprint for the rescue effort of other "terminally ill" populations of critically endangered species.

Population genetics of fruit bat reservoir informs the dynamics, distribution and diversity of Nipah virus

The structure and connectivity of wildlife host populations may influence zoonotic disease dynamics, evolution and therefore spillover risk to people. Fruit bats in the genus Pteropus, or flying foxes, are the primary natural reservoir for henipaviruses-a group of emerging paramyxoviruses that threaten livestock and public health. In Bangladesh, Pteropus medius is the reservoir for Nipah virus-and viral spillover has led to human fatalities nearly every year since 2001. Here, we use mitochondrial DNA and nuclear microsatellite markers to measure the population structure, demographic history and phylogeography of P. medius in Bangladesh. We combine this with a phylogeographic analysis of all known Nipah virus sequences and strains currently available to better inform the dynamics, distribution and evolutionary history of Nipah virus. We show that P. medius is primarily panmictic, but combined analysis of microsatellite and morphological data shows evidence for differentiation of two populations in eastern Bangladesh, corresponding to a divergent strain of Nipah virus also found in bats from eastern Bangladesh. Our demographic analyses indicate that a large, expanding population of flying foxes has existed in Bangladesh since the Late Pleistocene, coinciding with human population expansion in South Asia, suggesting repeated historical spillover of Nipah virus likely occurred. We present the first evidence of mitochondrial introgression, or hybridization, between P. medius and flying fox species found in South-East Asia (P. vampyrus and P. hypomelanus), which may help to explain the distribution of Nipah virus strains across the region.

Last exit before the brink: Conservation genomics of the Cambodian population of the critically endangered southern river terrapin

The southern river terrapin, Batagur affinis is one of the world's 25 most endangered freshwater turtle species. The major portion of the global population is currently found in peninsular Malaysia, with the only remnant Indochinese population in southern Cambodia. For more than a decade, wild nests in this remnant Cambodian population have been fenced and hatchlings reared in captivity. Here we amplified 10 microsatellite markers from all 136 captive individuals, obtained 2,658 presumably unlinked and neutral single nucleotide polymorphisms from 72 samples with ddRAD-seq, and amplified 784 bp of mtDNA from 50 samples. Our results reveal that the last Indochinese population comprised only four kinship groups as of 2012, with all offspring sired from <10 individuals in the wild. We demonstrate an obvious decrease in genetic contributions of breeders in the wild from 2006-2012 and identify high-value breeders instrumental for ex-situ management of the contemporary genetic stock of the species.

Mitochondrial DNA identified bat species in northeast India: electrocution mortality and biodiversity loss

Northeast India with two biodiversity hotspots is recognized as a biodiversity-rich region. However, several extant animals including chiropterans are currently at jeopardy due to habitat loss, electrocution mortality, and other anthropogenic threats. This study examines the efficacy of mitochondrial Cytochrome b (mtCytb) sequences for species-level identification of five electrocuted bat specimens from Manipur state. The similarity search results in the global database, Kimura 2 parameter (K2P) genetic distances, and neighbor-joining (NJ) tree identified all bat specimens into two species, Cynopterus sphinx and Megaerops niphanae. The detection of M. niphanae is the first record of this mammal from the state. In comparison with other Pteropodidae species, the genetic distances clearly discriminate both C. sphinx (7.9–30.2%) and M. niphanae (12.2–25.7%). In addition, the combined tree analysis of present and earlier genetic information of C. sphinx suggested the presence of cryptic lineages and sympatric population in India. This similar approach with more sampling from a wide distribution area could assist the future genetics research on chiropterans and their precise conservation.

High cryptic species diversity is revealed by genome-wide polymorphisms in a wild relative of banana, Musa itinerans, and implications for its conservation in subtropical China

BACKGROUND

Species delimitation is a challenging but essential task in conservation biology. Morphologically similar species are sometimes difficult to recognize even after examination by experienced taxonomists. With the advent of molecular approaches in species delimitation, this hidden diversity has received much recent attention. In addition to DNA barcoding approaches, analytical tools based on the multi-species coalescence model (MSC) have been developed for species delimitation. Musa itinerans is widely distributed in subtropical Asia, and at least six varieties have been documented. However, the number of evolutionarily distinct lineages remains unknown.

RESULTS

Using genome resequencing data of five populations (making up four varieties), we examined genome-wide variation and found four varieties that were evolutionary significant units. A Bayesian Phylogenetics and Phylogeography (BP&P) analysis using 123 single copy nuclear genes support three speciation events of M. itinerans varieties with robust posterior speciation probabilities; However, a Bayes factor delimitation of species with genomic data (BFD*) analysis using 1201 unlinked single nucleotide polymorphisms gave decisive support for a five-lineage model. When reconciling divergence time estimates with a speciation time scale, a modified three-lineage model was consistent with that of BP&P, in which the speciation time of two varieties (M. itinerans var. itinerans and M. itinerans var. lechangensis) were dated to 26.2 kya and 10.7 kya, respectively. In contrast, other two varieties (M. itinerans var. chinensis and M. itinerans var. guangdongensis) diverged only 3.8 kya in the Anthropocene; this may be a consequence of genetic drift rather than a speciation event.

CONCLUSION

Our results showed that the M. itinerans species complex harbours high cryptic species diversity. We recommend that M. itinerans var. itinerans and M. itinerans var. lechangensis be elevated to subspecies status, and the extremely rare latter subspecies be given priority for conservation. We also recommend that the very recently diverged M. itinerans var. chinensis and M. itinerans var. guangdongensis should be merged under the subspecies M. itinerans var. chinensis. Finally, we speculate that species delimitation of recently diverged lineages may be more effective using genome-wide bi-allelic SNP markers with BFD* than by using unlinked loci and BP&P.

Single-nucleotide polymorphism discovery and panel characterization in the African forest elephant

The continuing decline in forest elephant (Loxodonta cyclotis) numbers due to poaching and habitat reduction is driving the search for new tools to inform management and conservation. For dense rainforest species, basic ecological data on populations and threats can be challenging and expensive to collect, impeding conservation action in the field. As such, genetic monitoring is being increasingly implemented to complement or replace more burdensome field techniques. Single-nucleotide polymorphisms (SNPs) are particularly cost-effective and informative markers that can be used for a range of practical applications, including population census, assessment of human impact on social and genetic structure, and investigation of the illegal wildlife trade. SNP resources for elephants are scarce, but next-generation sequencing provides the opportunity for rapid, inexpensive generation of SNP markers in nonmodel species. Here, we sourced forest elephant DNA from 23 samples collected from 10 locations within Gabon, Central Africa, and applied double-digest restriction-site-associated DNA (ddRAD) sequencing to discover 31,851 tags containing SNPs that were reduced to a set of 1,365 high-quality candidate SNP markers. A subset of 115 candidate SNPs was then selected for assay design and validation using 56 additional samples. Genotyping resulted in a high conversion rate (93%) and a low per allele error rate (0.07%). This study provides the first panel of 107 validated SNP markers for forest elephants. This resource presents great potential for new genetic tools to produce reliable data and underpin a step-change in conservation policies for this elusive species.

Conservation priorities for endangered Indian tigers through a genomic lens

Tigers have lost 93% of their historical range worldwide. India plays a vital role in the conservation of tigers since nearly 60% of all wild tigers are currently found here. However, as protected areas are small (<300 km² on average), with only a few individuals in each, many of them may not be independently viable. It is thus important to identify and conserve genetically connected populations, as well as to maintain connectivity within them. We collected samples from wild tigers (Panthera tigris tigris) across India and used genome-wide SNPs to infer genetic connectivity. We genotyped 10,184 SNPs from 38 individuals across 17 protected areas and identified three genetically distinct clusters (corresponding to northwest, southern and central India). The northwest cluster was isolated with low variation and high relatedness. The geographically large central cluster included tigers from central, northeastern and northern India, and had the highest variation. Most genetic diversity (62%) was shared among clusters, while unique variation was highest in the central cluster (8.5%) and lowest in the northwestern one (2%). We did not detect signatures of differential selection or local adaptation. We highlight that the northwest population requires conservation attention to ensure persistence of these tigers.

Identifying hybridization and admixture using SNPs: application of the DArTseq platform in phylogeographic research on vertebrates

Next-generation sequencing (NGS) approaches are increasingly being used to generate multi-locus data for phylogeographic and evolutionary genetics research. We detail the applicability of a restriction enzyme-mediated genome complexity reduction approach with subsequent NGS (DArTseq) in vertebrate study systems at different evolutionary and geographical scales. We present two case studies using SNP data from the DArTseq molecular marker platform. First, we used DArTseq in a large phylogeographic study of the agamid lizard Ctenophorus caudicinctus, including 91 individuals and spanning the geographical range of this species across arid Australia. A low-density DArTseq assay resulted in 28 960 SNPs, with low density referring to a comparably reduced set of identified and sequenced markers as a cost-effective approach. Second, we applied this approach to an evolutionary genetics study of a classic frog hybrid zone (Litoria ewingii-Litoria paraewingi) across 93 individuals, which resulted in 48 117 and 67 060 SNPs for a low- and high-density assay, respectively. We provide a docker-based workflow to facilitate data preparation and analysis, then analyse SNP data using multiple methods including Bayesian model-based clustering and conditional likelihood approaches. Based on comparison of results from the DArTseq platform and traditional molecular approaches, we conclude that DArTseq can be used successfully in vertebrates and will be of particular interest to researchers working at the interface between population genetics and phylogenetics, exploring species boundaries, gene exchange and hybridization.

The effects of Pleistocene climate change on biotic differentiation in a montane songbird clade from Wallacea

The role of Pleistocene Ice Age in tropical diversification is poorly understood, especially in archipelagos, in which glaciation-induced sea level fluctuations may lead to complicated changes in land distribution. To assess how Pleistocene land bridges may have facilitated gene flow in tropical archipelagos, we investigated patterns of diversification in the rarely-collected rusty-bellied fantail Rhipidura teysmanni (Passeriformes: Rhipiduridae) complex from Wallacea using a combination of bioacoustic traits and whole-genome sequencing methods (dd-RADSeq). We report a biogeographic leapfrog pattern in the vocalizations of these birds, and uncover deep genomic divergence among island populations despite the presence of intermittent land connections between some. We demonstrate how rare instances of genetic introgression have affected the evolution of this species complex, and document the presence of double introgressive mitochondrial sweeps, highlighting the dangers of using only mitochondrial DNA in evolutionary research. By applying different tree inference approaches, we demonstrate how concatenation methods can give inaccurate results when investigating divergence in closely-related taxa. Our study highlights high levels of cryptic avian diversity in poorly-explored Wallacea, elucidates complex patterns of Pleistocene climate-mediated diversification in an elusive montane songbird, and suggests that Pleistocene land bridges may have accounted for limited connectivity among montane Wallacean biota.