Bat pluripotent stem cells reveal unusual entanglement between host and viruses.

Bats are distinctive among mammals due to their ability to fly, use laryngeal echolocation, and tolerate viruses. However, there are currently no reliable cellular models for studying bat biology or their response to viral infections. Here, we created induced pluripotent stem cells (iPSCs) from two species of bats: the wild greater horseshoe bat (Rhinolophus ferrumequinum) and the greater mouse-eared bat (Myotis myotis). The iPSCs from both bat species showed similar characteristics and had a gene expression profile resembling that of cells attacked by viruses. They also had a high number of endogenous viral sequences, particularly retroviruses. These results suggest that bats have evolved mechanisms to tolerate a large load of viral sequences and may have a more intertwined relationship with viruses than previously thought. Further study of bat iPSCs and their differentiated progeny will provide insights into bat biology, virus host relationships, and the molecular basis of bats' special traits.

Six reference-quality genomes reveal evolution of bat adaptations.

Bats possess extraordinary adaptations, including flight, echolocation, extreme longevity and unique immunity. High-quality genomes are crucial for understanding the molecular basis and evolution of these traits. Here we incorporated long-read sequencing and state-of-the-art scaffolding protocols1 to generate, to our knowledge, the first reference-quality genomes of six bat species (Rhinolophus ferrumequinum, Rousettus aegyptiacus, Phyllostomus discolor, Myotis myotis, Pipistrellus kuhlii and Molossus molossus). We integrated gene projections from our 'Tool to infer Orthologs from Genome Alignments' (TOGA) software with de novo and homology gene predictions as well as short- and long-read transcriptomics to generate highly complete gene annotations. To resolve the phylogenetic position of bats within Laurasiatheria, we applied several phylogenetic methods to comprehensive sets of orthologous protein-coding and noncoding regions of the genome, and identified a basal origin for bats within Scrotifera. Our genome-wide screens revealed positive selection on hearing-related genes in the ancestral branch of bats, which is indicative of laryngeal echolocation being an ancestral trait in this clade. We found selection and loss of immunity-related genes (including pro-inflammatory NF-κB regulators) and expansions of anti-viral APOBEC3 genes, which highlights molecular mechanisms that may contribute to the exceptional immunity of bats. Genomic integrations of diverse viruses provide a genomic record of historical tolerance to viral infection in bats. Finally, we found and experimentally validated bat-specific variation in microRNAs, which may regulate bat-specific gene-expression programs. Our reference-quality bat genomes provide the resources required to uncover and validate the genomic basis of adaptations of bats, and stimulate new avenues of research that are directly relevant to human health and disease1.

Comparative Genome Microsynteny Illuminates the Fast Evolution of Nuclear Mitochondrial Segments (NUMTs) in Mammals.

The escape of DNA from mitochondria into the nuclear genome (nuclear mitochondrial DNA, NUMT) is an ongoing process. Although pervasively observed in eukaryotic genomes, their evolutionary trajectories in a mammal-wide context are poorly understood. The main challenge lies in the orthology assignment of NUMTs across species due to their fast evolution and chromosomal rearrangements over the past 200 million years. To address this issue, we systematically investigated the characteristics of NUMT insertions in 45 mammalian genomes and established a novel, synteny-based method to accurately predict orthologous NUMTs and ascertain their evolution across mammals. With a series of comparative analyses across taxa, we revealed that NUMTs may originate from nonrandom regions in mtDNA, are likely found in transposon-rich and intergenic regions, and unlikely code for functional proteins. Using our synteny-based approach, we leveraged 630 pairwise comparisons of genome-wide microsynteny and predicted the NUMT orthology relationships across 36 mammals. With the phylogenetic patterns of NUMT presence-and-absence across taxa, we constructed the ancestral state of NUMTs given the mammal tree using a coalescent method. We found support on the ancestral node of Fereuungulata within Laurasiatheria, whose subordinal relationships are still controversial. This study broadens our knowledge on NUMT insertion and evolution in mammalian genomes and highlights the merit of NUMTs as alternative genetic markers in phylogenetic inference.

Chiropterans Are a Hotspot for Horizontal Transfer of DNA Transposons in Mammalia.

Horizontal transfer of transposable elements (TEs) is an important mechanism contributing to genetic diversity and innovation. Bats (order Chiroptera) have repeatedly been shown to experience horizontal transfer of TEs at what appears to be a high rate compared with other mammals. We investigated the occurrence of horizontally transferred (HT) DNA transposons involving bats. We found over 200 putative HT elements within bats; 16 transposons were shared across distantly related mammalian clades, and 2 other elements were shared with a fish and two lizard species. Our results indicate that bats are a hotspot for horizontal transfer of DNA transposons. These events broadly coincide with the diversification of several bat clades, supporting the hypothesis that DNA transposon invasions have contributed to genetic diversification of bats.

Hibernation telomere dynamics in a shifting climate: insights from wild greater horseshoe bats.

Hibernation is linked with various hypotheses to explain the extended lifespan of hibernating mammals compared with their non-hibernating counterparts. Studies on telomeres, markers of ageing and somatic maintenance, suggest telomere shortening slows during hibernation, and lengthening may reflect self-maintenance with favourable conditions. Bats in temperate zones adjust body temperatures during winter torpor to conserve energy and exploit mild conditions for foraging. Climate change may impact the hibernation cycle of bats, but more research is needed regarding the role of telomeres in understanding their response to a changing climate. Here, relative telomere length (rTL) was measured in the long-lived greater horseshoe bat Rhinolophus ferrumequinum (n = 223 individuals) over three winters, considering climatic conditions. Cross-sectional analyses revealed between-individual variation in rTL with a strong year effect, likely linked to varying weather conditions and foraging success. Additionally, within-individual increases of rTL occurred in 51% of consecutive measurements, with evidence of increasing telomerase expression during hibernation in this species. These findings highlight the beneficial effects of hibernation on telomeres and potential consequences of changing climatic conditions for long-lived temperate bats. Understanding the interplay between hibernation, telomeres, and climate can provide insights into the adaptive capacity and survival of bat populations facing environmental challenges.

Broad host range of SARS-CoV-2 predicted by comparative and structural analysis of ACE2 in vertebrates.

The novel coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of COVID-19. The main receptor of SARS-CoV-2, angiotensin I converting enzyme 2 (ACE2), is now undergoing extensive scrutiny to understand the routes of transmission and sensitivity in different species. Here, we utilized a unique dataset of ACE2 sequences from 410 vertebrate species, including 252 mammals, to study the conservation of ACE2 and its potential to be used as a receptor by SARS-CoV-2. We designed a five-category binding score based on the conservation properties of 25 amino acids important for the binding between ACE2 and the SARS-CoV-2 spike protein. Only mammals fell into the medium to very high categories and only catarrhine primates into the very high category, suggesting that they are at high risk for SARS-CoV-2 infection. We employed a protein structural analysis to qualitatively assess whether amino acid changes at variable residues would be likely to disrupt ACE2/SARS-CoV-2 spike protein binding and found the number of predicted unfavorable changes significantly correlated with the binding score. Extending this analysis to human population data, we found only rare (frequency <0.001) variants in 10/25 binding sites. In addition, we found significant signals of selection and accelerated evolution in the ACE2 coding sequence across all mammals, and specific to the bat lineage. Our results, if confirmed by additional experimental data, may lead to the identification of intermediate host species for SARS-CoV-2, guide the selection of animal models of COVID-19, and assist the conservation of animals both in native habitats and in human care.

Taking flight: An ecological, evolutionary and genomic perspective on bat telomeres.

Over 20% of all living mammals are bats (order Chiroptera). Bats possess extraordinary adaptations including powered flight, laryngeal echolocation and a unique immune system that enables them to tolerate a diversity of viral infections without presenting clinical disease symptoms. They occupy multiple trophic niches and environments globally. Significant physiological and ecological diversity occurs across the order. Bats also exhibit extreme longevity given their body size with many species showing few signs of ageing. The molecular basis of this extended longevity has recently attracted attention. Telomere maintenance potentially underpins bats' extended healthspan, although functional studies are still required to validate the causative mechanisms. In this review, we detail the current knowledge on bat telomeres, telomerase expression, and how these relate to ecology, longevity and life-history strategies. Patterns of telomere shortening and telomerase expression vary across species, and comparative genomic analyses suggest that alternative telomere maintenance mechanisms evolved in the longest-lived bats. We discuss the unique challenges faced when working with populations of wild bats and highlight ways to advance the field including expanding long-term monitoring across species that display contrasting life-histories and occupy different environmental niches. We further review how new high quality, chromosome-level genome assemblies can enable us to uncover the molecular mechanisms governing telomere dynamics and how phylogenomic analyses can reveal the adaptive significance of telomere maintenance and variation in bats.

Large-scale genome sampling reveals unique immunity and metabolic adaptations in bats.

Comprising more than 1,400 species, bats possess adaptations unique among mammals including powered flight, unexpected longevity, and extraordinary immunity. Some of the molecular mechanisms underlying these unique adaptations includes DNA repair, metabolism and immunity. However, analyses have been limited to a few divergent lineages, reducing the scope of inferences on gene family evolution across the Order Chiroptera. We conducted an exhaustive comparative genomic study of 37 bat species, one generated in this study, encompassing a large number of lineages, with a particular emphasis on multi-gene family evolution across immune and metabolic genes. In agreement with previous analyses, we found lineage-specific expansions of the APOBEC3 and MHC-I gene families, and loss of the proinflammatory PYHIN gene family. We inferred more than 1,000 gene losses unique to bats, including genes involved in the regulation of inflammasome pathways such as epithelial defence receptors, the natural killer gene complex and the interferon-gamma induced pathway. Gene set enrichment analyses revealed genes lost in bats are involved in defence response against pathogen-associated molecular patterns and damage-associated molecular patterns. Gene family evolution and selection analyses indicate bats have evolved fundamental functional differences compared to other mammals in both innate and adaptive immune system, with the potential to enhance antiviral immune response while dampening inflammatory signalling. In addition, metabolic genes have experienced repeated expansions related to convergent shifts to plant-based diets. Our analyses support the hypothesis that, in tandem with flight, ancestral bats had evolved a unique set of immune adaptations whose functional implications remain to be explored.

As Blind as a Bat? Opsin Phylogenetics Illuminates the Evolution of Color Vision in Bats.

Through their unique use of sophisticated laryngeal echolocation bats are considered sensory specialists amongst mammals and represent an excellent model in which to explore sensory perception. Although several studies have shown that the evolution of vision is linked to ecological niche adaptation in other mammalian lineages, this has not yet been fully explored in bats. Recent molecular analysis of the opsin genes, which encode the photosensitive pigments underpinning color vision, have implicated high-duty cycle (HDC) echolocation and the adoption of cave roosting habits in the degeneration of color vision in bats. However, insufficient sampling of relevant taxa has hindered definitive testing of these hypotheses. To address this, novel sequence data was generated for the SWS1 and MWS/LWS opsin genes and combined with existing data to comprehensively sample species representing diverse echolocation types and niches (SWS1 n = 115; MWS/LWS n = 45). A combination of phylogenetic analysis, ancestral state reconstruction, and selective pressure analyses were used to reconstruct the evolution of these visual pigments in bats and revealed that although both genes are evolving under purifying selection in bats, MWS/LWS is highly conserved but SWS1 is highly variable. Spectral tuning analyses revealed that MWS/LWS opsin is tuned to a long wavelength, 555-560 nm in the bat ancestor and the majority of extant taxa. The presence of UV vision in bats is supported by our spectral tuning analysis, but phylogenetic analyses demonstrated that the SWS1 opsin gene has undergone pseudogenization in several lineages. We do not find support for a link between the evolution of HDC echolocation and the pseudogenization of the SWS1 gene in bats, instead we show the SWS1 opsin is functional in the HDC echolocator, Pteronotus parnellii. Pseudogenization of the SWS1 is correlated with cave roosting habits in the majority of pteropodid species. Together these results demonstrate that the loss of UV vision in bats is more widespread than was previously considered and further elucidate the role of ecological niche specialization in the evolution of vision in bats.

AGILE: an assembled genome mining pipeline.

SUMMARY: A number of limiting factors mean that traditional genome annotation tools either fail or perform sub-optimally when trying to detect coding sequences in poor quality genome assemblies/genome reports. This means that potentially useful data is accessible only to those with specific skills and expertise in assembly and annotation. We present an Assembled-Genome mIning pipeLinE (AGILE) written in Perl that combines bioinformatics tools with a number of steps to overcome the limitations imposed by such assemblies when applied to highly fragmented genomes. Our methodology uses user-specified query genes from a closely related species to mine and annotate coding sequences that would traditionally be missed by standard annotation packages. Despite a focus on mammalian genomes, the generalized implementation means that it may be applied to any genome assembly, providing a means for non-specialists to gather gene sequences for downstream analyses. AVAILABILITY AND IMPLEMENTATION: Source code and associated files are available at: https://github.com/batlabucd/GenomeMining and https://bitbucket.org/BatlabUCD/genomemining/src. Singularity and Virtual Box images available at https://figshare.com/s/a0004bf93dc43484b0c0. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Drivers of longitudinal telomere dynamics in a long-lived bat species, Myotis myotis.

Age-related telomere shortening is considered a hallmark of the ageing process. However, a recent cross-sectional ageing study of relative telomere length (rTL) in bats failed to detect a relationship between rTL and age in the long-lived genus Myotis (M. myotis and M. bechsteinii), suggesting some other factors are responsible for driving telomere dynamics in these species. Here, we test if longitudinal rTL data show signatures of age-associated telomere attrition in M. myotis and differentiate which intrinsic or extrinsic factors are likely to drive telomere length dynamics. Using quantitative polymerase chain reaction, rTL was measured in 504 samples from a marked population, from Brittany, France, captured between 2013 and 2016. These represent 174 individuals with an age range of 0 to 7+ years. We find no significant relationship between rTL and age (p = .762), but demonstrate that within-individual rTL is highly variable from year to year. To investigate the heritability of rTL, a population pedigree (n = 1744) was constructed from genotype data generated from a 16-microsatellite multiplex, designed from an initial, low-coverage, Illumina genome for M. myotis. Heritability was estimated in a Bayesian, mixed model framework, and showed that little of the observed variance in rTL is heritable (h2  = 0.01-0.06). Rather, correlations of first differences, correlating yearly changes in telomere length and weather variables, demonstrate that, during the spring transition, average temperature, minimum temperature, rainfall and windspeed correlate with changes in longitudinal telomere dynamics. As such, rTL may represent a useful biomarker to quantify the physiological impact of various environmental stressors in bats.

The Birth and Death of Olfactory Receptor Gene Families in Mammalian Niche Adaptation.

The olfactory receptor (OR) gene families, which govern mammalian olfaction, have undergone extensive expansion and contraction through duplication and pseudogenization. Previous studies have shown that broadly defined environmental adaptations (e.g., terrestrial vs. aquatic) are correlated with the number of functional and non-functional OR genes retained. However, to date, no study has examined species-specific gene duplications in multiple phylogenetically divergent mammals to elucidate OR evolution and adaptation. Here, we identify the OR gene families driving adaptation to different ecological niches by mapping the fate of species-specific gene duplications in the OR repertoire of 94 diverse mammalian taxa, using molecular phylogenomic methods. We analyze >70,000 OR gene sequences mined from whole genomes, generated from novel amplicon sequencing data, and collated with data from previous studies, comprising one of the largest OR studies to date. For the first time, we demonstrate statistically significant patterns of OR species-specific gene duplications associated with the presence of a functioning vomeronasal organ. With respect to dietary niche, we uncover a novel link between a large number of duplications in OR family 5/8/9 and herbivory. Our results also highlight differences between social and solitary niches, indicating that a greater OR repertoire expansion may be associated with a solitary lifestyle. This study demonstrates the utility of species-specific duplications in elucidating gene family evolution, revealing how the OR repertoire has undergone expansion and contraction with respect to a number of ecological adaptations in mammals.

Evolutionary History of the Asian Horned Frogs (Megophryinae): Integrative Approaches to Timetree Dating in the Absence of a Fossil Record.

Molecular dating studies typically need fossils to calibrate the analyses. Unfortunately, the fossil record is extremely poor or presently nonexistent for many species groups, rendering such dating analysis difficult. One such group is the Asian horned frogs (Megophryinae). Sampling all generic nomina, we combined a novel ∼5 kb dataset composed of four nuclear and three mitochondrial gene fragments to produce a robust phylogeny, with an extensive external morphological study to produce a working taxonomy for the group. Expanding the molecular dataset to include out-groups of fossil-represented ancestral anuran families, we compared the priorless RelTime dating method with the widely used prior-based Bayesian timetree method, MCMCtree, utilizing a novel combination of fossil priors for anuran phylogenetic dating. The phylogeny was then subjected to ancestral phylogeographic analyses, and dating estimates were compared with likely biogeographic vicariant events. Phylogenetic analyses demonstrated that previously proposed systematic hypotheses were incorrect due to the paraphyly of genera. Molecular phylogenetic, morphological, and timetree results support the recognition of Megophryinae as a single genus, Megophrys, with a subgenus level classification. Timetree results using RelTime better corresponded with the known fossil record for the out-group anuran tree. For the priorless in-group, it also outperformed MCMCtree when node date estimates were compared with likely influential historical biogeographic events, providing novel insights into the evolutionary history of this pan-Asian anuran group. Given a relatively small molecular dataset, and limited prior knowledge, this study demonstrates that the computationally rapid RelTime dating tool may outperform more popular and complex prior reliant timetree methodologies.

ExUTR: a novel pipeline for large-scale prediction of 3'-UTR sequences from NGS data.

BACKGROUND: The three prime untranslated region (3'-UTR) is known to play a pivotal role in modulating gene expression by determining the fate of mRNA. Many crucial developmental events, such as mammalian spermatogenesis, tissue patterning, sex determination and neurogenesis, rely heavily on post-transcriptional regulation by the 3'-UTR. However, 3'-UTR biology seems to be a relatively untapped field, with only limited tools and 3'-UTR resources available. To elucidate the regulatory mechanisms of the 3'-UTR on gene expression, firstly the 3'-UTR sequences must be identified. Current 3'-UTR mining tools, such as GETUTR, 3USS and UTRscan, all depend on a well-annotated reference genome or curated 3'-UTR sequences, which hinders their application on a myriad of non-model organisms where the genomes are not available. To address these issues, the establishment of an NGS-based, automated pipeline is urgently needed for genome-wide 3'-UTR prediction in the absence of reference genomes. RESULTS: Here, we propose ExUTR, a novel NGS-based pipeline to predict and retrieve 3'-UTR sequences from RNA-Seq experiments, particularly designed for non-model species lacking well-annotated genomes. This pipeline integrates cutting-edge bioinformatics tools, databases (Uniprot and UTRdb) and novel in-house Perl scripts, implementing a fully automated workflow. By taking transcriptome assemblies as inputs, this pipeline identifies 3'-UTR signals based primarily on the intrinsic features of transcripts, and outputs predicted 3'-UTR candidates together with associated annotations. In addition, ExUTR only requires minimal computational resources, which facilitates its implementation on a standard desktop computer with reasonable runtime, making it affordable to use for most laboratories. We also demonstrate the functionality and extensibility of this pipeline using publically available RNA-Seq data from both model and non-model species, and further validate the accuracy of predicted 3'-UTR using both well-characterized 3'-UTR resources and 3P-Seq data. CONCLUSIONS: ExUTR is a practical and powerful workflow that enables rapid genome-wide 3'-UTR discovery from NGS data. The candidates predicted through this pipeline will further advance the study of miRNA target prediction, cis elements in 3'-UTR and the evolution and biology of 3'-UTRs. Being independent of a well-annotated reference genome will dramatically expand its application to much broader research area, encompassing all species for which RNA-Seq is available.

Identification of the same polyomavirus species in different African horseshoe bat species is indicative of short-range host-switching events.

Polyomaviruses (PyVs) are considered to be highly host-specific in different mammalian species, with no well-supported evidence for host-switching events. We examined the species diversity and host specificity of PyVs in horseshoe bats (Rhinolophus spp.), a broadly distributed and highly speciose mammalian genus. We annotated six PyV genomes, comprising four new PyV species, based on pairwise identity within the large T antigen (LTAg) coding region. Phylogenetic comparisons revealed two instances of highly related PyV species, one in each of the Alphapolyomavirus and Betapolyomavirus genera, present in different horseshoe bat host species (Rhinolophus blasii and R. simulator), suggestive of short-range host-switching events. The two pairs of Rhinolophus PyVs in different horseshoe bat host species were 99.9 and 88.8 % identical with each other over their respective LTAg coding sequences and thus constitute the same virus species. To corroborate the species identification of the bat hosts, we analysed mitochondrial cytb and a large nuclear intron dataset derived from six independent and neutrally evolving loci for bat taxa of interest. Bayesian estimates of the ages of the most recent common ancestors suggested that the near-identical and more distantly related PyV species diverged approximately 9.1E4 (5E3-2.8E5) and 9.9E6 (4E6-18E6) years before the present, respectively, in contrast to the divergence times of the bat host species: 12.4E6 (10.4E6-15.4E6). Our findings provide evidence that short-range host-switching of PyVs is possible in horseshoe bats, suggesting that PyV transmission between closely related mammalian species can occur.

Blood miRNomes and transcriptomes reveal novel longevity mechanisms in the long-lived bat, Myotis myotis.

BACKGROUND: Chiroptera, the bats, are the only order of mammals capable of true self-powered flight. Bats exhibit a number of other exceptional traits such as echolocation, viral tolerance and, perhaps most puzzlingly, extreme longevity given their body size. Little is known about the molecular mechanisms driving their extended longevity particularly at the levels of gene expression and post-transcriptional regulation. To elucidate the molecular mechanisms that may underlie their unusual longevity, we have deep sequenced 246.5 million small RNA reads from whole blood of the long-lived greater mouse-eared bats, Myotis myotis, and conducted a series of genome-wide comparative analyses between bat and non-bat mammals (human, pig and cow) in both blood miRNomes and transcriptomes, for the first time. RESULTS: We identified 539 miRNA gene candidates from bats, of which 468 unique mature miRNA were obtained. More than half of these miRNA (65.1 %) were regarded as bat-specific, regulating genes involved in the immune, ageing and tumorigenesis pathways. We have also developed a stringent pipeline for genome-wide miRNome comparisons across species, and identified 37 orthologous miRNA groups shared with bat, human, pig and cow, 6 of which were differentially expressed. For bats, 3 out of 4 up-regulated miRNA (miR-101-3p, miR-16-5p, miR-143-3p) likely function as tumor suppressors against various kinds of cancers, while one down-regulated miRNA (miR-221-5p) acts as a tumorigenesis promoter in human breast and pancreatic cancers. Additionally, a genome-wide comparison of mRNA transcriptomes across species also revealed specific gene expression patterns in bats. 127 up-regulated genes were enriched mainly in mitotic cell cycle and DNA repair mechanisms, while 364 down-regulated genes were involved primarily in mitochondrial activity. CONCLUSIONS: Our comprehensive and integrative analyses revealed bat-specific and differentially expressed miRNA and mRNA that function in key longevity pathways, producing a distinct bat gene expression pattern. For the first time, we show that bats may possess unique regulatory mechanisms for resisting tumorigenesis, repairing cellular damage and preventing oxidative stresses, all of which likely contribute to the extraordinary lifespan of Myotis myotis.

How and why overcome the impediments to resolution: lessons from rhinolophid and hipposiderid bats.

The phylogenetic and taxonomic relationships among the Old World leaf-nosed bats (Hipposideridae) and the closely related horseshoe bats (Rhinolophidae) remain unresolved. In this study, we generated a novel approximately 10-kb molecular data set of 19 nuclear exon and intron gene fragments for 40 bat species to elucidate the phylogenetic relationships within the families Rhinolophidae and Hipposideridae. We estimated divergence times and explored potential reasons for any incongruent phylogenetic signal. We demonstrated the effects of outlier taxa and genes on phylogenetic reconstructions and compared the relative performance of intron and exon data to resolve phylogenetic relationships. Phylogenetic analyses produced a well-resolved phylogeny, supporting the familial status of Hipposideridae and demonstrated the paraphyly of the largest genus, Hipposideros. A fossil-calibrated timetree and biogeographical analyses estimated that Rhinolophidae and Hipposideridae diverged in Africa during the Eocene approximately 42 Ma. The phylogram, the timetree, and a unique retrotransposon insertion supported the elevation of the subtribe Rhinonycterina to family level and which is diagnosed herein. Comparative analysis of diversification rates showed that the speciose genera Rhinolophus and Hipposideros underwent diversification during the Mid-Miocene Climatic Optimum. The intron versus exon analyses demonstrated the improved nodal support provided by introns for our optimal tree, an important finding for large-scale phylogenomic studies, which typically rely on exon data alone. With the recent outbreak of Middle East respiratory syndrome, caused by a novel coronavirus, the study of these species is urgent as they are considered the natural reservoir for emergent severe acute respiratory syndrome (SARS)-like coronaviruses. It has been shown that host phylogeny is the primary factor that determines a virus's persistence, replicative ability, and can act as a predictor of new emerging disease. Therefore, this newly resolved phylogeny can be used to direct future assessments of viral diversity and to elucidate the origin and development of SARS-like coronaviruses in mammals.

Wing Membrane Biopsies for Bat Cytogenetics: Finding of 2n = 54 in Irish Rhinolophushipposideros (Rhinolophidae, Chiroptera, Mammalia) Supports Two Geographically Separated Chromosomal Variants in Europe.

In Europe, 2 different diploid chromosome numbers, 2n = 54 and 2n = 56, have been described in the lesser horseshoe bat (Rhinolophushipposideros). The eastern form with 2n = 56 extends from the Czech Republic to Greece. To date, specimens with 54 chromosomes have been reported only from Spain and Germany. This study expands the distributional area of the western variant to Ireland. Strikingly, this distribution of European chromosomal variants is in contrast to the available molecular data that indicate little genetic differentiation of R. hipposideros populations spanning Northwestern to Central Europe. Further, we have developed an optimized protocol for establishing fibroblast cell cultures, suitable for karyotype analyses, from 3-mm wing membrane biopsies. This is a useful technique for cytogenetic studies of endangered bat species, as this non-lethal sampling method imposes only minimum stress to the animal without lasting adverse effects and is routinely used to sample tissue probes for molecular genetic studies in bats.

Nuclear introns outperform mitochondrial DNA in inter-specific phylogenetic reconstruction: Lessons from horseshoe bats (Rhinolophidae: Chiroptera).

Despite many studies illustrating the perils of utilising mitochondrial DNA in phylogenetic studies, it remains one of the most widely used genetic markers for this purpose. Over the last decade, nuclear introns have been proposed as alternative markers for phylogenetic reconstruction. However, the resolution capabilities of mtDNA and nuclear introns have rarely been quantified and compared. In the current study we generated a novel ∼5kb dataset comprising six nuclear introns and a mtDNA fragment. We assessed the relative resolution capabilities of the six intronic fragments with respect to each other, when used in various combinations together, and when compared to the traditionally used mtDNA. We focused on a major clade in the horseshoe bat family (Afro-Palaearctic clade; Rhinolophidae) as our case study. This old, widely distributed and speciose group contains a high level of conserved morphology. This morphological stasis renders the reconstruction of the phylogeny of this group with traditional morphological characters complex. We sampled multiple individuals per species to represent their geographic distributions as best as possible (122 individuals, 24 species, 68 localities). We reconstructed the species phylogeny using several complementary methods (partitioned Maximum Likelihood and Bayesian and Bayesian multispecies-coalescent) and made inferences based on consensus across these methods. We computed pairwise comparisons based on Robinson-Foulds tree distance metric between all Bayesian topologies generated (27,000) for every gene(s) and visualised the tree space using multidimensional scaling (MDS) plots. Using our supported species phylogeny we estimated the ancestral state of key traits of interest within this group, e.g. echolocation peak frequency which has been implicated in speciation. Our results revealed many potential cryptic species within this group, even in taxa where this was not suspected a priori and also found evidence for mtDNA introgression. We demonstrated that by using just two introns one can recover a better supported species tree than when using the mtDNA alone, despite the shorter overall length of the combined introns. Additionally, when combining any single intron with mtDNA, we showed that the result is highly similar to the mtDNA gene tree and far from the true species tree and therefore this approach should be avoided. We caution against the indiscriminate use of mtDNA in phylogenetic studies and advocate for pilot studies to select nuclear introns. The selection of marker type and number is a crucial step that is best based on critical examination of preliminary or previously published data. Based on our findings and previous publications, we recommend the following markers to recover phylogenetic relationships between recently diverged taxa (<20 My) in bats and other mammals: ACOX2, COPS7A, BGN, ROGDI and STAT5A.

Patterns of orofacial clefting in the facial morphology of bats: a possible naturally occurring model of cleft palate.

A normal feature of the facial anatomy of many species of bat is the presence of bony discontinuities or clefts, which bear a remarkable similarity to orofacial clefts that occur in humans as a congenital pathology. These clefts occur in two forms: a midline cleft between the two premaxillae (analogous to the rare midline craniofacial clefts in humans) and bilateral paramedian clefts between the premaxilla and the maxillae (analogous to the typical cleft lip and palate in humans). Here, we describe the distribution of orofacial clefting across major bat clades, exploring the relationship of the different patterns of clefting to feeding mode, development of the vomeronasal organ, development of the nasolacrimal duct and mode of emission of the echolocation call in different bat groups. We also present the results of detailed radiographic and soft tissue dissections of representative examples of the two types of cleft. The midline cleft has arisen independently multiple times in bat phylogeny, whereas the paramedian cleft has arisen once and is a synapomorphy uniting the Rhinolophidae and Hipposideridae. In all cases examined, the bony cleft is filled in by a robust fibrous membrane, continuous with the periosteum of the margins of the cleft. In the paramedian clefts, this membrane splits to enclose the premaxilla but forms a loose fold laterally between the premaxilla and maxilla, allowing the premaxilla and nose-leaf to pivot dorsoventrally in the sagittal plane under the action of facial muscles attached to the nasal cartilages. It is possible that this is a specific adaptation for echolocation and/or aerial insectivory. Given the shared embryological location of orofacial clefts in bats and humans, it is likely that aspects of the developmental control networks that produce cleft lip and palate in humans may also be implicated in the formation of these clefts as a normal feature in some bats. A better understanding of craniofacial development in bats with and without clefts may therefore suggest avenues for research into abnormal craniofacial development in humans.