Conflicting Evolutionary Histories of the Mitochondrial and Nuclear Genomes in New World Myotis Bats

A comprehensive sampling of mitogenomes shows the utility to infer phylogeny of termites (Blattodea: Termitoidae)

The mitogenome sequence data have been widely used in inferring the phylogeny of insects. In this study, we determined the complete mitogenome for Macrotermes sp. (Termitidae, Macrotermitinae) using next-generation sequencing. Macrotermes sp. possesses a typical insect mitogenome, displaying an identical gene order and gene content to other existing termite mitogenomes. We present the first prediction of the secondary structure of ribosomal RNA genes in termites. The rRNA secondary structures of Macrotermes sp. exhibit similarities to closely related insects and also feature distinctive characteristics in their helical structures. Together with 321 published mitogenomes of termites as ingroups and 8 cockroach mitogenomes as outgroups, we compiled the most comprehensive mitogenome sequence matrix for Termitoidae to date. Phylogenetic analyses were conducted using datasets employing different data coding strategies and various inference methods. Robust relationships were recovered at the family or subfamily level, demonstrating the utility of comprehensive mitogenome sampling in resolving termite phylogenies. The results supported the monophyly of Termitoidae, and consistent relationships within this group were observed across different analyses. Mastotermitidae was consistently recovered as the sister group to all other termite families. The families Hodotermitidae, Stolotermitidae, and Archotermopsidae formed the second diverging clade, followed by the Kalotermitidae. The Neoisoptera was consistently supported with strong node support, with Stylotermitidae being sister to the remaining families. Rhinotermitidae was found to be non-monophyletic, and Serritermitidae nested within the basal clades of Rhinotermitidae and was sister to Psammotermitinae. Overall, our phylogenetic results are largely consistent with earlier mitogenome studies.

Karyotypic stasis and swarming influenced the evolution of viral tolerance in a species-rich bat radiation

The emergence of COVID-19 and severe acute respiratory syndrome (SARS) has prioritized understanding bats' viral tolerance. Myotis bats are exceptionally species rich and have evolved viral tolerance. They also exhibit swarming, a cryptic behavior where large, multi-species assemblages gather for mating, which has been hypothesized to promote interspecific hybridization. To resolve the coevolution of genome architecture and their unusual antiviral tolerance, we undertook a phylogenomic analysis of 60 Old World Myotis genomes. We demonstrate an extensive history of introgressive hybridization that has replaced the species phylogeny across 17%-93% of the genome except for pericentromeric regions of macrochromosomes. Introgression tracts were enriched on microchromosome regions containing key antiviral pathway genes overexpressed during viral challenge experiments. Together, these results suggest that the unusual Myotis karyotype may have evolved to selectively position immune-related genes in high recombining genomic regions prone to introgression of divergent alleles, including a diversity of interleukin loci responsible for the release of pro-inflammatory cytokines.

Mito-nuclear discordance within Anthozoa, with notes on unique properties of their mitochondrial genomes

Whole mitochondrial genomes are often used in phylogenetic reconstruction. However, discordant patterns in species relationships between mitochondrial and nuclear phylogenies are commonly observed. Within Anthozoa (Phylum Cnidaria), mitochondrial (mt)-nuclear discordance has not yet been examined using a large and comparable dataset. Here, we used data obtained from target-capture enrichment sequencing to assemble and annotate mt genomes and reconstruct phylogenies for comparisons to phylogenies inferred from hundreds of nuclear loci obtained from the same samples. The datasets comprised 108 hexacorals and 94 octocorals representing all orders and > 50% of extant families. Results indicated rampant discordance between datasets at every taxonomic level. This discordance is not attributable to substitution saturation, but rather likely caused by introgressive hybridization and unique properties of mt genomes, including slow rates of evolution driven by strong purifying selection and substitution rate variation. Strong purifying selection across the mt genomes caution their use in analyses that rely on assumptions of neutrality. Furthermore, unique properties of the mt genomes were noted, including genome rearrangements and the presence of nad5 introns. Specifically, we note the presence of the homing endonuclease in ceriantharians. This large dataset of mitochondrial genomes further demonstrates the utility of off-target reads generated from target-capture data for mt genome assembly and adds to the growing knowledge of anthozoan evolution.

Population Genomic Evidence for the Diversification of Bellamya aeruginosa in Different River Systems in China

Clarifying the genetic structure can facilitate the understanding of a species evolution history. It is crucial for the management of germplasm resources and providing useful guidance for effective selective breeding. Bellamya is an economically and ecologically important freshwater snail for fish, birds and even humans. Population genetic structures of the Bellamya species, however, were unknown in previous studies. Population genomics approaches with tens to hundreds of thousands of single nucleotide polymorphisms (SNPs) make it possible to detect previously unidentified structures. The population genomic study of seven populations of B. aeruginosa across three river systems (Yellow River, Yangtze River and Pearl River) in China was conducted by SLAF-seq. SLAF-seq obtained a total of 4737 polymorphisms SLAF-tags and 25,999 high-consistency genome-wide SNPs. The population genetic structure showed a clear division among populations from the Yellow River basin (YH and WL) and the Pearl River basin (QSH and LB), as well as population YC from the Yangtze River basin using the SNPs data. However, there existed no distinct population structure using the mitochondrial DNA (mtDNA). The anthropogenic translocation from the Yangtze River basin to the Pearl River basin and the passive dispersion from the Yangtze River basin to the Yellow River basin by flooding have weakened the phylogeographic pattern of B. aeruginosa. The divergence of B. aeruginosa in the three river systems suggests that the anthropogenic dispersal for aquaculture and breeding requires serious consideration of the population structure for the preservation of genetic diversity and effective utilization of germplasm resources.

Whole-genome Analysis Reveals Contrasting Relationships Among Nuclear and Mitochondrial Genomes Between Three Sympatric Bat Species

Understanding mechanisms involved in speciation can be challenging, especially when hybridization or introgression blurs species boundaries. In bats, resolving relationships of some closely related groups has proved difficult due subtle interspecific variation both in morphometrics and molecular data sets. The endemic South American Histiotus bats, currently considered a subgenus of Eptesicus, harbor unresolved phylogenetic relationships and of those is a trio consisting of two closely related species: Eptesicus (Histiotus) macrotus and Eptesicus (Histiotus) montanus, and their relationship with a third, Eptesicus (Histiotus) magellanicus. The three sympatric species bear marked resemblance to each other, but can be differentiated morphologically. Furthermore, previous studies have been unable to differentiate the species from each other at a molecular level. In order to disentangle the phylogenetic relationships of these species, we examined the differentiation patterns and evolutionary history of the three Eptesicus (H.) species at the whole-genome level. The nuclear DNA statistics between the species suggest strong gene flow and recent hybridization between E. (H.) montanus and E. (H.) macrotus, whereas E. (H.) magellanicus shows a higher degree of isolation. In contrast, mitochondrial DNA shows a closer relationship between E. (H.) magellanicus and E. (H.) montanus. Opposing patterns in mtDNA and nuclear markers are often due to differences in dispersal, and here it could be both as a result of isolation in refugia during the last glacial maximum and female philopatry and male-biased dispersal. In conclusion, this study shows the importance of both the nuclear and mitochondrial DNA in resolving phylogenetic relationships and species histories.

Systematics of the New World bats Eptesicus and Histiotus suggest trans-marine dispersal followed by Neotropical cryptic diversification

Biodiversity can be boosted by colonization of new habitats such as remote islands and separated continents. Molecular studies have suggested that recently evolved organisms probably colonized already separated continents by dispersal, either via land bridge connections or crossing the ocean. Here we test the on-land and trans-marine dispersal hypotheses by evaluating possibilities of colonization routes over the Bering land bridge and across the Atlantic Ocean in the cosmopolitan bat genus Eptesicus (Chiroptera, Vespertilionidae). Previous molecular studies have found New World Eptesicus more closely related to Histiotus, a Neotropical endemic lineage with enlarged ears, than to Old World Eptesicus. However, phylogenetic relationships within the New World group remained unresolved and their evolutionary history was unclear. Here we studied the systematics of New World Eptesicus and Histiotus using extensive taxonomic and geographic sampling, and genomic data from thousands of ultra-conserved elements (UCEs). We estimated phylogenetic trees using concatenation and multispecies coalescent. All analyses supported four major New World clades and a novel topology where E. fuscus and Histiotus are sister clades that together diverged from two sister clades of Neotropical Eptesicus. Intra-clade divergence suggested cryptic diversity that has been concealed by morphological features, especially in the Neotropics where taxonomic re-evaluations are warranted. Molecular dating estimated that Old World and New World clades diverged around 17 million years ago followed by radiation of major New World clades in the mid-Miocene, when climatic changes might have facilitated global dispersal and radiation events. Biogeographic ancestral reconstruction supported the Neotropical origin of the New World clades, suggesting a trans-Atlantic colonization route from North Africa to the northern Neotropics. We highlight that trans-marine dispersal may be more prevalent than currently acknowledged and may be an important first step to global biodiversification.

Old specimens for old branches: Assessing effects of sample age in resolving a rapid Neotropical radiation of squirrels

Ultraconserved Elements (UCEs) have been useful to resolve challenging phylogenies of non-model clades, unpuzzling long-conflicted relationships in key branches of the Tree of Life at both deep and shallow levels. UCEs are often reliably recovered from historical samples, unlocking a vast number of preserved natural history specimens for analysis. However, the extent to which sample age and preservation method impact UCE recovery as well as downstream inferences remains unclear. Furthermore, there is an ongoing debate on how to curate, filter, and properly analyze UCE data when locus recovery is uneven across sample age and quality. In the present study we address these questions with an empirical dataset composed of over 3800 UCE loci from 219 historical and modern samples of Sciuridae, a globally distributed and ecologically important family of rodents. We provide a genome-scale phylogeny of two squirrel subfamilies (Sciurillinae and Sciurinae: Sciurini) and investigate their placement within Sciuridae. For historical specimens, recovery of UCE loci and mean length per locus were inversely related to sample age; deeper sequencing improved the number of UCE loci recovered but not locus length. Most of our phylogenetic inferences-performed on six datasets with alternative data-filtering strategies, and using three distinct optimality criteria-resulted in distinct topologies. Datasets containing more loci (40% and 50% taxa representativeness matrices) yielded more concordant topologies and higher support values than strictly filtered datasets (60% matrices) particularly with IQ-Tree and SVDquartets, while filtering based on information content provided better topological resolution for inferences with the coalescent gene-tree based approach in ASTRAL-III. We resolved deep relationships in Sciuridae (including among the five currently recognized subfamilies) and relationships among the deepest branches of Sciurini, but conflicting relationships remain at both genus- and species-levels for the rapid Neotropical tree squirrel radiation. Our results suggest that phylogenomic consensus can be difficult and heavily influenced by the age of available samples and the filtering steps used to optimize dataset properties.

SINE-Based Phylogenomics Reveal Extensive Introgression and Incomplete Lineage Sorting in Myotis

Using presence/absence data from over 10,000 Ves SINE insertions, we reconstructed a phylogeny for 11 Myotis species. With nearly one-third of individual Ves gene trees discordant with the overall species tree, phylogenetic conflict appears to be rampant in this genus. From the observed conflict, we infer that ILS is likely a major contributor to the discordance. Much of the discordance can be attributed to the hypothesized split between the Old World and New World Myotis clades and with the first radiation of Myotis within the New World. Quartet asymmetry tests reveal signs of introgression between Old and New World taxa that may have persisted until approximately 8 MYA. Our introgression tests also revealed evidence of both historic and more recent, perhaps even contemporary, gene flow among Myotis species of the New World. Our findings suggest that hybridization likely played an important role in the evolutionary history of Myotis and may still be happening in areas of sympatry. Despite limitations arising from extreme discordance, our SINE-based phylogeny better resolved deeper relationships (particularly the positioning of M. brandtii) and was able to identify potential introgression pathways among the Myotis species sampled.

Variation in the Myotis keaysi complex (Chiroptera, Vespertilionidae), with description of a new species from Ecuador

The genus Myotis comprises a diverse group of vespertilionid bats with worldwide distribution. Twenty-eight Neotropical species are currently recognized. In this study, we evaluate molecular and morphological variation in the M. keaysi complex, a high elevation clade of Neotropical myotine bats characterized by complex taxonomy and high morphological variation. A phylogeny inferred with cytochrome-b sequences recovered two clades composed of samples traditionally assigned to M. keaysi, with 9% of genetic divergence between them. These clades were also suggested as putative distinct species by molecular species delimitation methods. Qualitative and quantitative morphological analyses indicated a phenotypic discontinuity between specimens from central Andes (including the holotype of M. keaysi) and western lowlands of Ecuador, showing strong congruence between molecular and morphological approaches. We therefore describe a new species for the Tumbes-Choco-Magdalena region, documenting their external and cranial diagnostic characters by comparing them with other Neotropical species. In addition, we provide an emended diagnosis for our new concept of M. keaysi.

A new species of Myotis (Chiroptera, Vespertilionidae) from Uruguay

Abstract

The genus Myotis comprises a diverse group of vesper bats with worldwide distribution. Twenty-eight neotropical species are currently recognized. Based on a morphological approach, we describe a new species of Myotis from the Uruguayan Pampas grasslands, an ecoregion under high anthropogenic pressure with a largely unknown bat fauna. Qualitative and quantitative morphological analyses support the recognition of the new species and we present a set of external and cranial diagnostic characters by comparing them with other neotropical Myotis species. The new species reassembles Myotis riparius, but can be distinguished by a set of qualitative and quantitative morphological traits, including its clearly bicolored dorsal fur, tricolored ventral fur, a pelage on the dorsal surface of uropatagium, sagittal crest lower, braincase lower in lateral view and overall smaller size.

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.

Comment on “Population genetics reveal Myotis keenii (Keen’s myotis) and Myotis evotis (long-eared myotis) to be a single species”

Genetic exchange and hybridization appear common among the western long-eared bats from North America. Multiple sources of evidence indicate that lineages within this group are evolving independently, despite genetic exchange. However, evidence of gene flow raises questions about the species-level status of some lineages. C.L. Lausen et al. (2019. Can. J. Zool. 97(3): 267–279) proposed that Myotis evotis (H. Allen, 1864) (long-eared myotis) and Myotis keenii (Merriam, 1895) (Keen’s myotis) are one species, not two. This conclusion is based on analyses of cytochrome b and microsatellite data suggesting gene flow between these taxa. Microsatellites are not reliable markers for identifying species because homoplasy can be a major confounding factor, which appears to be true in this case. We reanalyzed the dataset of C.L. Lausen et al. (2019) and show that it is not reliable to distinguish between gene flow or homoplasy, and that these data do not support the conclusion that M. evotis and M. keenii represent a single species. Previous morphological and genomic studies indicate that these are separate species despite previous genetic exchange between them. Failing to recognize that gene flow can occur between independently evolving lineages is counterproductive for conservation because it can lead to neglect of important independent lineages, and likewise failing to use proper tools to delimit species is counterproductive to efforts to quantify biodiversity and design conservation strategies.

First record of disk-footed bat Eudiscopus denticulus (Chiroptera, Vespertilionidae) from China and resolution of phylogenetic position of the genus

The disk-footed bat Eudiscopus denticulus(Osgood, 1932) is a rare species in Southeast Asia. During two chiropteran surveys in the summer of 1981 and 2019, eight and three small Myotis-like bats with distinct disk-like hindfeet were collected from Yunnan Province, China, respectively. External, craniodental, and phylogenetic evidence confirmed these specimens as E. denticulus, representing a new genus in China. The complete mitochondrial genome consistently showed robust support for E. denticulus as a basal lineage within Myotinae. The coding patterns and characteristics of its mitochondrial genome were similar to that of other published genomes fromMyotis. The echolocation signals of the newly collected individuals were analyzed. The potential distribution range of Eudiscopus in Southeast Asia inferred using the MaxEnt model indicated its potential occurrence along the southern border region of Yunnan, China.

Molecular and morphological revision of small Myotinae from the Himalayas shed new light on the poorly known genus Submyotodon (Chiroptera: Vespertilionidae)

The systematics status of the constituent species of the M. mystacinus morphogroup in the Himalayan region has long been marred by uncertainty. Lack of integrative studies combining morphological and genetic data from specimens recently collected in this region has hampered our understanding of cryptic variations in this complex taxonomic group. To address this issue, new material from the Himalayan region of India and Nepal was obtained and vouchered specimens in the holdings of various museums were also re-examined. As comparative material, a large series of relevant specimens from South and Southeast Asia were also included in this revision. Using a combination of multivariate analysis of craniodental characters and molecular reconstructions, we critically evaluated the systematic position of the small Myotinae in the Himalayas. We establish that M. nipalensis forms a very distinct lineage (which also includes the recently described M. annatessae) and refute previous taxonomic suggestions that it is related to M. davidii. Our study also conclusively proved the common occurrence of the poorly known genus Submyotodon in the Himalayan region (Afghanistan, Pakistan, India, Nepal and China) and evidenced species-level divergences within that genus. Submyotodon species share nyctalodont or semi-nyctalodont lower molar configuration with few other small and unrelated Myotinae from Asia suggesting that these unusual dental characters are homoplasious in this subfamily. We also noticed a very confused taxonomic situation associated with many DNA sequences of Asian Myotis deposited in public repositories and call for possibilities of better data curation.

Heterogeneous taxonomic resolution of cytochrome b gene identification of bats from Argentina: Implications for field studies

Bats are among the most diverse, widespread, and abundant mammals. In Argentina, 67 species of bats have been recorded, belonging to 5 families and 29 genera. These high levels of biodiversity are likely to complicate identification at fieldwork, especially between closely related species, where external morphology-based approaches are the only immediate means for a priori species assignment. The use of molecular markers can enhance species identification, and acquires particular relevance in capture-release studies. In this study, we discuss the extent of the use of the mitochondrial cytochrome b gene for species identification, comparing external morphology identification with a molecular phylogenetic classification based on this marker, under the light of current bat systematics. We analyzed 33 samples collected in an eco-epidemiological survey in the province of Santa Fe (Argentina). We further sequenced 27 museum vouchers to test the accuracy of cytochrome b -based phylogenies in taxonomic identification of bats occurring in the Pampean/Chacoan regions of Argentina. The cytochrome b gene was successfully amplified in all Molossid and Vespertilionid species except for Eptesicus, for which we designed a new reverse primer. The resulting Bayesian phylogeny was congruent with current systematics. Cytochrome b proved useful for species-level delimitation in non-conflicting genera (Eumops, Dasypterus, Molossops) and has infrageneric resolution in more complex lineages (Eptesicus, Myotis, Molossus). We discuss four sources of incongruence that may act separately or in combination: 1) molecular processes, 2) biology, 3) limitations in identification, and 4) errors in the current taxonomy. The present study confirms the general applicability of cytochrome b -based phylogenies in eco-epidemiological studies, but its resolution and reliability depend mainly, but not solely, on the level of genetic differentiation within each bat genus.

Target-capture phylogenomics provide insights on gene and species tree discordances in Old World treefrogs (Anura: Rhacophoridae)

Genome-scale data have greatly facilitated the resolution of recalcitrant nodes that Sanger-based datasets have been unable to resolve. However, phylogenomic studies continue to use traditional methods such as bootstrapping to estimate branch support; and high bootstrap values are still interpreted as providing strong support for the correct topology. Furthermore, relatively little attention has been given to assessing discordances between gene and species trees, and the underlying processes that produce phylogenetic conflict. We generated novel genomic datasets to characterize and determine the causes of discordance in Old World treefrogs (Family: Rhacophoridae)-a group that is fraught with conflicting and poorly supported topologies among major clades. Additionally, a suite of data filtering strategies and analytical methods were applied to assess their impact on phylogenetic inference. We showed that incomplete lineage sorting was detected at all nodes that exhibited high levels of discordance. Those nodes were also associated with extremely short internal branches. We also clearly demonstrate that bootstrap values do not reflect uncertainty or confidence for the correct topology and, hence, should not be used as a measure of branch support in phylogenomic datasets. Overall, we showed that phylogenetic discordances in Old World treefrogs resulted from incomplete lineage sorting and that species tree inference can be improved using a multi-faceted, total-evidence approach, which uses the most amount of data and considers results from different analytical methods and datasets.

Parallel ddRAD and Genome Skimming Analyses Reveal a Radiative and Reticulate Evolutionary History of the Temperate Bamboos

Rapid evolutionary radiations are among the most challenging phylogenetic problems, wherein different types of data (e.g., morphology and molecular) or genetic markers (e.g., nuclear and organelle) often yield inconsistent results. The tribe Arundinarieae, that is, the temperate bamboos, is a clade of tetraploid originated 22 Ma and subsequently radiated in East Asia. Previous studies of Arundinarieae have found conflicting relationships and/or low support. Here, we obtain nuclear markers from ddRAD data for 213 Arundinarieae taxa and parallel sampling of chloroplast genomes from genome skimming for 147 taxa. We first assess the feasibility of using ddRAD-seq data for phylogenetic estimates of paleopolyploid and rapidly radiated lineages, optimize clustering thresholds, and analysis workflow for orthology identification. Reference-based ddRAD data assembly approaches perform well and yield strongly supported relationships that are generally concordant with morphology-based taxonomy. We recover five major lineages, two of which are notable (the pachymorph and leptomorph lineages), in that they correspond with distinct rhizome morphologies. By contrast, the phylogeny from chloroplast genomes differed significantly. Based on multiple lines of evidence, the ddRAD tree is favored as the best species tree estimation for temperate bamboos. Using a time-calibrated ddRAD tree, we find that Arundinarieae diversified rapidly around the mid-Miocene corresponding with intensification of the East Asian monsoon and the evolution of key innovations including the leptomorph rhizomes. Our results provide a highly resolved phylogeny of Arundinarieae, shed new light on the radiation and reticulate evolutionary history of this tribe, and provide an empirical example for the study of recalcitrant plant radiations. [Arundinarieae; ddRAD; paleopolyploid; genome skimming; rapid diversification; incongruence.]

The published complete mitochondrial genome of Eptesicus serotinus is a chimera of Vespertilio sinensis and Hypsugo alaschanicus (Mammalia: Chiroptera)

The mitogenome of Eptesicus serotinus (Serotine bat) was published in 2013 with GenBank accession number KF111725 and NCBI Reference Sequence number NC_022474. This sequence was placed with Vespertilio sinensis (Asian parti-colored bat) in a COI gene tree but with Hypsugo alashanicus (Alashanian pipistrelle) in a cytochrome b gene tree. Direct comparison of mitogenomes showed that 92.4% of this mitogenome is similar to Vespertilio sinensis, 5.9% to Hypsugo alaschanicus, and that 1.6% of the mitogenome could not be attributed to either species, or any other species. This mitogenome has been re-used in at least 17 phylogenies. Our findings suggest that mitogenomes are best verified with multiple gene trees, followed by direct comparison of sequences. We conclude that greater vigilance is warranted to ensure that problematic sequences do not enter the scientific record, and are not re-used in subsequent studies.

Larger, unfiltered datasets are more effective at resolving phylogenetic conflict: Introns, exons, and UCEs resolve ambiguities in Golden-backed frogs (Anura: Ranidae; genus Hylarana)

Using FrogCap, a recently-developed sequence-capture protocol, we obtained >12,000 highly informative exons, introns, and ultraconserved elements (UCEs), which we used to illustrate variation in evolutionary histories of these classes of markers, and to resolve long-standing systematic problems in Southeast Asian Golden-backed frogs of the genus-complex Hylarana. We also performed a comprehensive suite of analyses to assess the relative performance of different genetic markers, data filtering strategies, tree inference methods, and different measures of branch support. To reduce gene tree estimation error, we filtered the data using different thresholds of taxon completeness (missing data) and parsimony informative sites (PIS). We then estimated species trees using concatenated datasets and Maximum Likelihood (IQ-TREE) in addition to summary (ASTRAL-III), distance-based (ASTRID), and site-based (SVDQuartets) multispecies coalescent methods. Topological congruence and branch support were examined using traditional bootstrap, local posterior probabilities, gene concordance factors, quartet frequencies, and quartet scores. Our results did not yield a single concordant topology. Instead, introns, exons, and UCEs clearly possessed different phylogenetic signals, resulting in conflicting, yet strongly-supported phylogenetic estimates. However, a combined analysis comprising the most informative introns, exons, and UCEs converged on a similar topology across all analyses, with the exception of SVDQuartets. Bootstrap values were consistently high despite high levels of incongruence and high proportions of gene trees supporting conflicting topologies. Although low bootstrap values did indicate low heuristic support, high bootstrap support did not necessarily reflect congruence or support for the correct topology. This study reiterates findings of some previous studies, which demonstrated that traditional bootstrap values can produce positively misleading measures of support in large phylogenomic datasets. We also showed a remarkably strong positive relationship between branch length and topological congruence across all datasets, implying that very short internodes remain a challenge to resolve, even with orders of magnitude more data than ever before. Overall, our results demonstrate that more data from unfiltered or combined datasets produced superior results. Although data filtering reduced gene tree incongruence, decreased amounts of data also biased phylogenetic estimation. A point of diminishing returns was evident, at which higher congruence (from more stringent filtering) at the expense of amount of data led to topological error as assessed by comparison to more complete datasets across different genomic markers. Additionally, we showed that applying a parameter-rich model to a partitioned analysis of concatenated data produces better results compared to unpartitioned, or even partitioned analysis using model selection. Despite some lingering uncertainties, a combined analysis of our genomic data and sequences supplemented from GenBank (on the basis of a few gene regions) revealed highly supported novel systematic arrangements. Based on these new findings, we transfer Amnirana nicobariensis into the genus Indosylvirana; and I. milleti and Hylarana celebensis to the genus Papurana. We also provisionally place H. attigua in the genus Papurana pending verification from positively identified (voucher substantiated) samples.

Patterns of mtDNA introgression suggest population replacement in Palaearctic whiskered bat species

Secondary contacts can play a major role in the evolutionary histories of species. Various taxa diverge in allopatry and later on come into secondary contact during range expansions. When they meet, their interactions and the extent of gene flow depend on the level of their ecological differentiation and the strength of their reproductive isolation. In this study, we present the multilocus phylogeography of two cryptic whiskered bat species, Myotis mystacinus and M. davidii, with a particular focus on their putative sympatric zone. Our findings suggest that M. mystacinus and M. davidii evolved in allopatry and came into secondary contact during range expansions. Individuals in the area of secondary contact, in Anatolia and the Balkans, have discordant population assignments based on the mitochondrial and the nuclear datasets. These observed patterns suggest that the local M. mystacinus populations hybridized with expanding M. davidii populations, which resulted in mitochondrial introgression from the former. In the introgression area, M. mystacinus individuals with concordant nuclear and mitochondrial genotypes were identified in relatively few locations, suggesting that the indigenous populations might have been largely replaced by invading M. davidii. Changing environmental conditions coupled with ecological competition is the likely reason for this replacement. Our study presents one possible example of a historical population replacement that was captured in phylogeographic patterns.

An Annotated Chromosome-Level Reference Genome of the Red-Eared Slider Turtle (Trachemys scripta elegans)

Among vertebrates, turtles have many unique characteristics providing biologists with opportunities to study novel evolutionary innovations and processes. We present here a high-quality, partially phased, and chromosome-level Red-Eared Slider (Trachemys scripta elegans, TSE) genome as a reference for future research on turtle and tetrapod evolution. This TSE assembly is 2.269 Gb in length, has one of the highest scaffold N50 and N90 values of any published turtle genome to date (N50 = 129.68 Mb and N90 = 19 Mb), and has a total of 28,415 annotated genes. We introduce synteny analyses using BUSCO single-copy orthologs, which reveal two chromosome fusion events accounting for differences in chromosome counts between emydids and other cryptodire turtles and reveal many fission/fusion events for birds, crocodiles, and snakes relative to TSE. This annotated chromosome-level genome will provide an important reference genome for future studies on turtle, vertebrate, and chromosome evolution.

NGS-Indel Coder: A pipeline to code indel characters in phylogenomic data with an example of its application in milkweeds (Asclepias)

Targeted genome sequencing approaches allow characterization of evolutionary relationships using a considerable number of nuclear genes and informative characters. However, most phylogenomic analyses only utilize single nucleotide polymorphisms (SNPs). Studies at the species level, especially in groups that have recently radiated, often recover low amounts of phylogenetically informative variation in coding regions, and require non-coding sequences, which are richer in indels, to resolve gene trees. Here, NGS-Indel Coder, a pipeline to detect and omit false positive indels inferred from assemblies of short read sequence data, was developed to resolve the relationships among and within major clades of the American milkweeds (Asclepias), which are the result of a rapid and recent evolutionary radiation, and whose phylogeny has been difficult to resolve. This pipeline was applied to a Hyb-Seq data set of 768 loci including targeted exons and flanking intron regions from 33 milkweed species. Robust species tree inference was improved by excluding small alignment partitions (<100 bp) that increased gene tree ambiguity and incongruence. To further investigate the robustness of indel coding, data sets that included small and large indels were explored, and species trees derived from concatenated loci versus coalescent methods based on gene trees were compared. The phylogeny of Asclepias obtained using nuclear data was well resolved, and phylogenetic information from indels improved resolution of specific nodes. The Temperate North American, Mexican Highland, and Incarnatae clades were well supported as monophyletic. Asclepias coulteri, which has been considered part of the Sonoran Desert clade based on plastome analyses, was placed as sister to all the other milkweed species studied here, rather than as a member of that clade. Two groups within the Temperate North American and Mexican clades were not resolved, and the inferred relationships strongly conflicted when comparing results based on data sets that did or did not include indel characters. This new pipeline represents a step forward in making maximal use of the information content in phylogenomic data sets.

Phylogeography and niche modelling: reciprocal enlightenment

Phylogeography examines the spatial genetic structure of species. Environmental niche modelling (or ecological niche modelling; ENM) examines the environmental limits of a species’ ecological niche. These two fields have great potential to be used together. ENM can shed light on how phylogeographical patterns develop and help identify possible drivers of spatial structure that need to be further investigated. Specifically, ENM can be used to test for niche differentiation among clades, identify factors limiting individual clades and identify barriers and contact zones. It can also be used to test hypotheses regarding the effects of historical and future climate change on spatial genetic patterns by projecting niches using palaeoclimate or future climate data. Conversely, phylogeographical information can populate ENM with within-species genetic diversity. Where adaptive variation exists among clades within a species, modelling their niches separately can improve predictions of historical distribution patterns and future responses to climate change. Awareness of patterns of genetic diversity in niche modelling can also alert conservationists to the potential loss of genetically diverse areas in a species’ range. Here, we provide a simplistic overview of both fields, and focus on their potential for integration, encouraging researchers on both sides to take advantage of the opportunities available.

Population genetics revealMyotis keenii(Keen’s myotis) andMyotis evotis(long-eared myotis) to be a single species

Recognizing delineations of gene flow among groups of animals can be challenging but is necessary for conservation and management. Of particular importance is the identification of species boundaries. Several physical and genetic traits have been used with mixed success to distinguish Myotis keenii (Merriam, 1895) (Keen’s myotis) and Myotis evotis (H. Allen, 1864) (long-eared myotis), but it is unclear whether species distinction is biologically warranted. We generated 12–14 microsatellite locus genotypes for 275 long-eared Myotis representing four species — M. keenii, M. evotis, Myotis septentrionalis (Trouessart, 1897) (northern myotis), and Myotis thysanodes Miller, 1897 (fringed myotis) — from across northwestern North America and 23 Myotis lucifugus (Le Conte, 1831) (little brown myotis) as the outgroup. Population genetic analyses revealed four well-defined groups (species): M. septentrionalis, M. thysanodes, M. lucifugus, and a single group comprising M. keenii and M. evotis. We document high rates of gene flow within M. evotis/M. keenii. Cytochrome b gene (mtDNA) sequencing failed to resolve morphologically identifiable species. We highlight the importance of geographically thorough investigation of genetic connectivity (nuclear markers) when assessing taxonomic status of closely related groups. We document a morphometric cline within M. evotis/M. keenii that may in part explain earlier analyses that led to the description of the smaller bodied M. keenii (type locality: Haida Gwaii). We conclude that M. keenii does not qualify as a genetic or biological species.

Porous barriers? Assessment of gene flow within and among sympatric long-eared bat species

Species are the basic units for measuring biodiversity and for comprehending biological interactions. Yet, their delineation is often contentious, especially in groups that are both diverse and phenotypically conservative. Three cryptic species of long-eared bats, Plecotus auritus, P. austriacus, and P. macrobullaris, co-occur over extensive areas of Western Europe. The latter is a fairly recent discovery, questioning the overall diversity of the entire Plecotus complex. Yet, high morphological and acoustic similarities compromise the reliable identification of long-eared bats in the field. We postulate that such extensive phenotypic overlap, along with the recurrent observation of morphologically intermediate individuals, may hide rampant interspecific hybridization. Based on a geographic sampling centered on areas of sympatry in the Alps and Corsica, we assessed the level of reproductive isolation of these three Plecotus species with mitochondrial and nuclear markers, looking at both inter- and intraspecific genetic population structuring. No sign of hybridization was detected between these three species that appear well separated biologically. Genetic structuring of populations, however, reflected different species-specific responses to environmental connectivity, that is, to the presence of orographic or sea barriers. While the Alpine range and the Ligurian Sea coincided with sharp genetic discontinuities in P. macrobullaris and P. austriacus, the more ubiquitous P. auritus showed no significant population structuration. There were clear phylogeographic discrepancies between microsatellite and mitochondrial markers at the intraspecific level, however, which challenges the reliance on simple barcoding approaches for the delineation of sound conservation units.

Validity of the Diplostomoidea and Diplostomida (Digenea, Platyhelminthes) upheld in phylogenomic analysis

Higher systematics within the Digenea, Carus 1863 have been relatively stable since a phylogenetic analysis of partial nuclear ribosomal markers (rDNA) led to the erection of the Diplostomida Olson, Cribb, Tkach, Bray, and Littlewood, 2003. However, recent mitochondrial (mt) genome phylogenies suggest this order might be paraphyletic. These analyses show members of two diplostomidan superfamilies are more closely related to the Plagiorchiida La Rue, 1957 than to other members of the Diplostomida. A recent phylogeny based on partial cytochrome c oxidase I also indicates one of the groups implicated, the Diplostomoidea Poirier, 1886, is non-monophyletic. To determine if these results were robust to additional taxon sampling, we analyzed mt genomes from seven diplostomoids in three families. To choose between phylogenetic alternatives based on mt genomes and the prior rDNA-based topology, we analyzed hundreds of ultra-conserved genomic elements assembled from shotgun sequencing. The Diplostomida was paraphyletic in the mt genome phylogeny but supported in the ultra-conserved genomic element phylogeny. We speculate this mitonuclear discordance is related to ancient, rapid radiation in the Digenea. Both ultra-conserved genomic elements and mt genomes support the monophyly of the Diplostomoidea and show congruent relationships within it. The Cyathocotylidae Mühling, 1898 are early diverging descendants of a paraphyletic clade of Diplostomidae Poirier, 1886, in which are nested members of the Strigeidae Railliet, 1919; the results support prior suggestions that the Crassiphialinae Sudarikov, 1960 will rise to the family level. Morphological traits of diplostomoid metacercariae appear to be more useful for differentiating clades than those of adults. We describe a new species of Cotylurus Szidat, 1928, resurrect a species of Hysteromorpha Lutz, 1931, and find support for a species of Alaria Schrank, 1788 of contested validity. Complete rDNA operons from seven diplostomoid species are provided as a resource for future studies.

Systematics of Clupeiformes and testing for ecological limits on species richness in a trans-marine/freshwater clade

ABSTRACT Clupeiformes (herring, sardines, shad, anchovies and allies) are a globally distributed clade with nearly 400 marine, freshwater, and diadromous species. Although best known as filter feeding fishes that form large schools, this group occupies a diverse array of trophic guilds and habitats. Theory suggests that species richness in clades is modulated by ecological limits, which results in diversity-dependent clade growth, a pattern that most clades exhibit. As a trans-marine/freshwater clade that has undergone repeated transitions between marine and freshwaters, Clupeiformes are an excellent system for investigating the interplay between ecological diversity and macroevolutionary dynamics. In this study we review the systematics of Clupeiformes and explore discordance in phylogenetic relationships and divergence times between mitochondrial and nuclear loci. We then use comparative methods to test whether ecological limits regulate diversity in Clupeiformes. We find discordance in phylogenetic relationships at various taxonomic scales, but also considerable agreement between genomes. Our results suggest that trans-marine/freshwater clades are able to circumvent ecological limits on clade growth at regional, but not on local scales. Our study demonstrates that phylogenies are a critical link between ecology and macroevolutionary dynamics, and suggests habitat transitions can play a key role in shaping diversity patterns, particularly in the neotropics.

Bat Biology, Genomes, and the Bat1K Project: To Generate Chromosome-Level Genomes for All Living Bat Species

Bats are unique among mammals, possessing some of the rarest mammalian adaptations, including true self-powered flight, laryngeal echolocation, exceptional longevity, unique immunity, contracted genomes, and vocal learning. They provide key ecosystem services, pollinating tropical plants, dispersing seeds, and controlling insect pest populations, thus driving healthy ecosystems. They account for more than 20% of all living mammalian diversity, and their crown-group evolutionary history dates back to the Eocene. Despite their great numbers and diversity, many species are threatened and endangered. Here we announce Bat1K, an initiative to sequence the genomes of all living bat species (n∼1,300) to chromosome-level assembly. The Bat1K genome consortium unites bat biologists (>148 members as of writing), computational scientists, conservation organizations, genome technologists, and any interested individuals committed to a better understanding of the genetic and evolutionary mechanisms that underlie the unique adaptations of bats. Our aim is to catalog the unique genetic diversity present in all living bats to better understand the molecular basis of their unique adaptations; uncover their evolutionary history; link genotype with phenotype; and ultimately better understand, promote, and conserve bats. Here we review the unique adaptations of bats and highlight how chromosome-level genome assemblies can uncover the molecular basis of these traits. We present a novel sequencing and assembly strategy and review the striking societal and scientific benefits that will result from the Bat1K initiative.