Integrating gene annotation with orthology inference at scale.

Annotating coding genes and inferring orthologs are two classical challenges in genomics and evolutionary biology that have traditionally been approached separately, limiting scalability. We present TOGA (Tool to infer Orthologs from Genome Alignments), a method that integrates structural gene annotation and orthology inference. TOGA implements a different paradigm to infer orthologous loci, improves ortholog detection and annotation of conserved genes compared with state-of-the-art methods, and handles even highly fragmented assemblies. TOGA scales to hundreds of genomes, which we demonstrate by applying it to 488 placental mammal and 501 bird assemblies, creating the largest comparative gene resources so far. Additionally, TOGA detects gene losses, enables selection screens, and automatically provides a superior measure of mammalian genome quality. TOGA is a powerful and scalable method to annotate and compare genes in the genomic era.

Contradictory Phylogenetic Signals in the Laurasiatheria Anomaly Zone.

Relationships among laurasiatherian clades represent one of the most highly disputed topics in mammalian phylogeny. In this study, we attempt to disentangle laurasiatherian interordinal relationships using two independent genome-level approaches: (1) quantifying retrotransposon presence/absence patterns, and (2) comparisons of exon datasets at the levels of nucleotides and amino acids. The two approaches revealed contradictory phylogenetic signals, possibly due to a high level of ancestral incomplete lineage sorting. The positions of Eulipotyphla and Chiroptera as the first and second earliest divergences were consistent across the approaches. However, the phylogenetic relationships of Perissodactyla, Cetartiodactyla, and Ferae, were contradictory. While retrotransposon insertion analyses suggest a clade with Cetartiodactyla and Ferae, the exon dataset favoured Cetartiodactyla and Perissodactyla. Future analyses of hitherto unsampled laurasiatherian lineages and synergistic analyses of retrotransposon insertions, exon and conserved intron/intergenic sequences might unravel the conflicting patterns of relationships in this major mammalian clade.

Gene losses in the common vampire bat illuminate molecular adaptations to blood feeding.

Vampire bats are the only mammals that feed exclusively on blood. To uncover genomic changes associated with this dietary adaptation, we generated a haplotype-resolved genome of the common vampire bat and screened 27 bat species for genes that were specifically lost in the vampire bat lineage. We found previously unknown gene losses that relate to reduced insulin secretion (FFAR1 and SLC30A8), limited glycogen stores (PPP1R3E), and a unique gastric physiology (CTSE). Other gene losses likely reflect the biased nutrient composition (ERN2 and CTRL) and distinct pathogen diversity of blood (RNASE7) and predict the complete lack of cone-based vision in these strictly nocturnal bats (PDE6H and PDE6C). Notably, REP15 loss likely helped vampire bats adapt to high dietary iron levels by enhancing iron excretion, and the loss of CYP39A1 could have contributed to their exceptional cognitive abilities. These findings enhance our understanding of vampire bat biology and the genomic underpinnings of adaptations to blood feeding.

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.

A highly contiguous genome assembly of the bat hawkmoth Hyles vespertilio (Lepidoptera: Sphingidae).

BACKGROUND: Adapted to different ecological niches, moth species belonging to the Hyles genus exhibit a spectacular diversity of larval color patterns. These species diverged ∼7.5 million years ago, making this rather young genus an interesting system to study a wide range of questions including the process of speciation, ecological adaptation, and adaptive radiation. RESULTS: Here we present a high-quality genome assembly of the bat hawkmoth Hyles vespertilio, the first reference genome of a member of the Hyles genus. We generated 51× Pacific Biosciences long reads with an average read length of 8.9 kb. Pacific Biosciences reads longer than 4 kb were assembled into contigs, resulting in a 651.4-Mb assembly consisting of 530 contigs with an N50 value of 7.5 Mb. The circular mitochondrial contig has a length of 15,303 bp. The H. vespertilio genome is very repeat-rich and exhibits a higher repeat content (50.3%) than other Bombycoidea species such as Bombyx mori (45.7%) and Manduca sexta (27.5%). We developed a comprehensive gene annotation workflow to obtain consensus gene models from different evidence including gene projections, protein homology, transcriptome data, and ab initio predictions. The resulting gene annotation is highly complete with 94.5% of BUSCO genes being completely present, which is higher than the BUSCO completeness of the B. mori (92.2%) and M. sexta (90%) annotations. CONCLUSIONS: Our gene annotation strategy has general applicability to other genomes, and the H. vespertilio genome provides a valuable molecular resource to study a range of questions in this genus, including phylogeny, incomplete lineage sorting, speciation, and hybridization. A genome browser displaying the genome, alignments, and annotations is available at https://genome-public.pks.mpg.de/cgi-bin/hgTracks?db=HLhylVes1.

Longitudinal comparative transcriptomics reveals unique mechanisms underlying extended healthspan in bats.

Bats are the longest-lived mammals, given their body size. However, the underlying molecular mechanisms of their extended healthspans are poorly understood. To address this question we carried out an eight-year longitudinal study of ageing in long-lived bats (Myotis myotis). We deep-sequenced ~1.7 trillion base pairs of RNA from 150 blood samples collected from known aged bats to ascertain the age-related transcriptomic shifts and potential microRNA-directed regulation that occurred. We also compared ageing transcriptomic profiles between bats and other mammals by analysis of 298 longitudinal RNA sequencing datasets. Bats did not show the same transcriptomic changes with age as commonly observed in humans and other mammals, but rather exhibited a unique, age-related gene expression pattern associated with DNA repair, autophagy, immunity and tumour suppression that may drive their extended healthspans. We show that bats have naturally evolved transcriptomic signatures that are known to extend lifespan in model organisms, and identify novel genes not yet implicated in healthy ageing. We further show that bats' longevity profiles are partially regulated by microRNA, thus providing novel regulatory targets and pathways for future ageing intervention studies. These results further disentangle the ageing process by highlighting which ageing pathways contribute most to healthy ageing in mammals.

Recurrent loss of HMGCS2 shows that ketogenesis is not essential for the evolution of large mammalian brains.

Apart from glucose, fatty acid-derived ketone bodies provide metabolic energy for the brain during fasting and neonatal development. We investigated the evolution of HMGCS2, the key enzyme required for ketone body biosynthesis (ketogenesis). Unexpectedly, we found that three mammalian lineages, comprising cetaceans (dolphins and whales), elephants and mastodons, and Old World fruit bats have lost this gene. Remarkably, many of these species have exceptionally large brains and signs of intelligent behavior. While fruit bats are sensitive to starvation, cetaceans and elephants can still withstand periods of fasting. This suggests that alternative strategies to fuel large brains during fasting evolved repeatedly and reveals flexibility in mammalian energy metabolism. Furthermore, we show that HMGCS2 loss preceded brain size expansion in toothed whales and elephants. Thus, while ketogenesis was likely important for brain size expansion in modern humans, ketogenesis is not a universal precondition for the evolution of large mammalian brains.

Population level mitogenomics of long-lived bats reveals dynamic heteroplasmy and challenges the Free Radical Theory of Ageing.

Bats are the only mammals capable of true, powered flight, which drives an extremely high metabolic rate. The "Free Radical Theory of Ageing" (FTRA) posits that a high metabolic rate causes mitochondrial heteroplasmy and the progressive ageing phenotype. Contrary to this, bats are the longest-lived order of mammals given their small size and high metabolic rate. To investigate if bats exhibit increased mitochondrial heteroplasmy with age, we performed targeted, deep sequencing of mitogenomes and measured point heteroplasmy in wild, long lived Myotis myotis. Blood was sampled from 195 individuals, aged between <1 and at 6+ years old, and whole mitochondria deep-sequenced, with a subset sampled over multiple years. The majority of heteroplasmies were at a low frequency and were transitions. Oxidative mutations were present in only a small number of individuals, suggesting local oxidative stress events. Cohort data showed no significant increase in heteroplasmy with age, while longitudinal data from recaptured individuals showed heteroplasmy is dynamic, and does not increase uniformly over time. We show that bats do not suffer from the predicted, inevitable increase in heteroplasmy as posited by the FRTA, instead heteroplasmy was found to be dynamic, questioning its presumed role as a primary driver of ageing.

Growing old, yet staying young: The role of telomeres in bats' exceptional longevity.

Understanding aging is a grand challenge in biology. Exceptionally long-lived animals have mechanisms that underpin extreme longevity. Telomeres are protective nucleotide repeats on chromosome tips that shorten with cell division, potentially limiting life span. Bats are the longest-lived mammals for their size, but it is unknown whether their telomeres shorten. Using >60 years of cumulative mark-recapture field data, we show that telomeres shorten with age in Rhinolophus ferrumequinum and Miniopterus schreibersii, but not in the bat genus with greatest longevity, Myotis. As in humans, telomerase is not expressed in Myotis myotis blood or fibroblasts. Selection tests on telomere maintenance genes show that ATM and SETX, which repair and prevent DNA damage, potentially mediate telomere dynamics in Myotis bats. Twenty-one telomere maintenance genes are differentially expressed in Myotis, of which 14 are enriched for DNA repair, and 5 for alternative telomere-lengthening mechanisms. We demonstrate how telomeres, telomerase, and DNA repair genes have contributed to the evolution of exceptional longevity in Myotis bats, advancing our understanding of healthy aging.

The complete mitochondrial genome of the Greater Mouse-Eared bat, Myotis myotis (Chiroptera: Vespertilionidae).

In this study, we report the complete mitochondrial genome of the Greater Mouse-Eared Bat, Myotis myotis. The mitogenome is 17 213 bp with base composition A (34.2%), G (13%), C (22.4%), and T (30.5%). The genome shows conserved synteny with other mammalian mitogenomes, containing 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and 1 control region (D-Loop). The majority of the genes are encoded on the H-Strand, except for ND6 and eight tRNAs. All protein-coding genes start with the ATG start codon, except for ND2, ND3, and ND5 which begin with ATT or ATA. Seven protein-coding genes terminated in a canonical stop codon, TAA or TAG, five contain incomplete stop codons, T or TA. Cytochrome b terminates in the mitochondria specific stop codon AGA. This mitogenome provides a valuable resource for future studies of M. myotis and other bat and mammal species.

The complete mitochondrial genome of the Bechstein's bat, Myotis bechsteinii (Chiroptera, Vespertilionidae).

In this study, we present the complete mitochondrial genome of the Bechstein's bat, Myotis bechsteinii. The mitogenome is 17,151 bp in length and is AT-rich with base composition A (27.8%), C (22%), G (16.1%), and T (34.1%). The mitogenome shows conserved gene content and order similar with other mammalian mitogenomes, being composed of 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and one control region. The majority of genes are encoded on the H-Strand except for ND6 and 8 tRNAs as found in other bat species. The field identification of Myotis bechsteinii was confirmed by phylogenetic analyses using datasets comprising whole mitogenomes and COXI. This mitogenome is a resource for future studies of Myotis bats and other mammals.

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.

The complete mitochondrial genome of Kuhl's pipistrelle, Pipistrellus kuhlii (Chiroptera: Vespertilionidae).

The Kuhl's pipistrelle (Pipistrellus kuhlii) is a small, vespertilionid bat species, with a large range extending from the Iberian Peninsula into the Near East and the Arabian Peninsula. In this study, we determine for the first time the complete mitogenome of this species. The P. kuhlii mitogenome is 16,991 base pairs long with 37 genes and 1 control region, showing conserved gene content and order with other vertebrate mitogenomes. The length of the 22 tRNA genes ranges between 60 bp (tRNA-Ser) and 75 bp (tRNA-Leu). The D-loop region is 1553 bp long with low CG content (39.8%).

A nonlethal sampling method to obtain, generate and assemble whole blood transcriptomes from small, wild mammals.

The acquisition of tissue samples from wild populations is a constant challenge in conservation biology, especially for endangered species and protected species where nonlethal sampling is the only option. Whole blood has been suggested as a nonlethal sample type that contains a high percentage of bodywide and genomewide transcripts and therefore can be used to assess the transcriptional status of an individual, and to infer a high percentage of the genome. However, only limited quantities of blood can be nonlethally sampled from small species and it is not known if enough genetic material is contained in only a few drops of blood, which represents the upper limit of sample collection for some small species. In this study, we developed a nonlethal sampling method, the laboratory protocols and a bioinformatic pipeline to sequence and assemble the whole blood transcriptome, using Illumina RNA-Seq, from wild greater mouse-eared bats (Myotis myotis). For optimal results, both ribosomal and globin RNAs must be removed before library construction. Treatment of DNase is recommended but not required enabling the use of smaller amounts of starting RNA. A large proportion of protein-coding genes (61%) in the genome were expressed in the blood transcriptome, comparable to brain (65%), kidney (63%) and liver (58%) transcriptomes, and up to 99% of the mitogenome (excluding D-loop) was recovered in the RNA-Seq data. In conclusion, this nonlethal blood sampling method provides an opportunity for a genomewide transcriptomic study of small, endangered or critically protected species, without sacrificing any individuals.

Transgene- and locus-dependent imprinting reveals allele-specific chromosome conformations.

When positioned into the integrin α-6 gene, an Hoxd9lacZ reporter transgene displayed parental imprinting in mouse embryos. While the expression from the paternal allele was comparable with patterns seen for the same transgene when present at the neighboring HoxD locus, almost no signal was scored at this integration site when the transgene was inherited from the mother, although the Itga6 locus itself is not imprinted. The transgene exhibited maternal allele-specific DNA hypermethylation acquired during oogenesis, and its expression silencing was reversible on passage through the male germ line. Histone modifications also corresponded to profiles described at known imprinted loci. Chromosome conformation analyses revealed distinct chromatin microarchitectures, with a more compact structure characterizing the maternally inherited repressed allele. Such genetic analyses of well-characterized transgene insertions associated with a de novo-induced parental imprint may help us understand the molecular determinants of imprinting.