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.

Maximum depth sequencing reveals an ON/OFF replication slippage switch and apparent in vivo selection for bifidobacterial pilus expression.

The human gut microbiome, of which the genus Bifidobacterium is a prevalent and abundant member, is thought to sustain and enhance human health. Several surface-exposed structures, including so-called sortase-dependent pili, represent important bifidobacterial gut colonization factors. Here we show that expression of two sortase-dependent pilus clusters of the prototype Bifidobacterium breve UCC2003 depends on replication slippage at an intragenic G-tract, equivalents of which are present in various members of the Bifidobacterium genus. The nature and extent of this slippage is modulated by the host environment. Involvement of such sortase-dependent pilus clusters in microbe-host interactions, including bacterial attachment to the gut epithelial cells, has been shown previously and is corroborated here for one case. Using a Maximum Depth Sequencing strategy aimed at excluding PCR and sequencing errors introduced by DNA polymerase reagents, specific G-tract sequences in B. breve UCC2003 reveal a range of G-tract lengths whose plasticity within the population is functionally utilized. Interestingly, replication slippage is shown to be modulated under in vivo conditions in a murine model. This in vivo modulation causes an enrichment of a G-tract length which appears to allow biosynthesis of these sortase-dependent pili. This work provides the first example of productive replication slippage influenced by in vivo conditions. It highlights the potential for microdiversity generation in "beneficial" gut commensals.

Distinguishing activated T regulatory cell and T conventional cells by single-cell technologies.

Resting conventional T cells (Tconv) can be distinguished from T regulatory cells (Treg) by the canonical markers FOXP3, CD25 and CD127. However, the expression of these proteins alters after T-cell activation leading to overlap between Tconv and Treg. The objective of this study was to distinguish resting and antigen-responsive T effector (Tconv) and Treg using single-cell technologies. CD4+ Treg and Tconv cells were stimulated with antigen and responsive and non-responsive populations processed for targeted and non-targeted single-cell RNAseq. Machine learning was used to generate a limited set of genes that could distinguish responding and non-responding Treg and Tconv cells and which was used for single-cell multiplex qPCR and to design a flow cytometry panel. Targeted scRNAseq clearly distinguished the four-cell populations. A minimal set of 27 genes was identified by machine learning algorithms to provide discrimination of the four populations at >95% accuracy. In all, 15 of the genes were validated to be differentially expressed by single-cell multiplex qPCR. Discrimination of responding Treg from responding Tconv could be achieved by a flow cytometry strategy that included staining for CD25, CD127, FOXP3, IKZF2, ITGA4, and the novel marker TRIM which was strongly expressed in Tconv and weakly expressed in both responding and non-responding Treg. A minimal set of genes was identified that discriminates responding and non-responding CD4+ Treg and Tconv cells and, which have identified TRIM as a marker to distinguish Treg by flow cytometry.

Interspecies transcriptomics identify genes that underlie disproportionate foot growth in jerboas.

Despite the great diversity of vertebrate limb proportion and our deep understanding of the genetic mechanisms that drive skeletal elongation, little is known about how individual bones reach different lengths in any species. Here, we directly compare the transcriptomes of homologous growth cartilages of the mouse (Mus musculus) and bipedal jerboa (Jaculus jaculus), the latter of which has "mouse-like" arms but extremely long metatarsals of the feet. Intersecting gene-expression differences in metatarsals and forearms of the two species revealed that about 10% of orthologous genes are associated with the disproportionately rapid elongation of neonatal jerboa feet. These include genes and enriched pathways not previously associated with endochondral elongation as well as those that might diversify skeletal proportion in addition to their known requirements for bone growth throughout the skeleton. We also identified transcription regulators that might act as "nodes" for sweeping differences in genome expression between species. Among these, Shox2, which is necessary for proximal limb elongation, has gained expression in jerboa metatarsals where it has not been detected in other vertebrates. We show that Shox2 is sufficient to increase mouse distal limb length, and a nearby putative cis-regulatory region is preferentially accessible in jerboa metatarsals. In addition to mechanisms that might directly promote growth, we found evidence that jerboa foot elongation may occur in part by de-repressing latent growth potential. The genes and pathways that we identified here provide a framework to understand the modular genetic control of skeletal growth and the remarkable malleability of vertebrate limb proportion.

Functional and metabolic fitness of human CD4+ T lymphocytes during metabolic stress.

Human CD4+ T cells are essential mediators of immune responses. By altering the mitochondrial and metabolic states, we defined metabolic requirements of human CD4+ T cells for in vitro activation, expansion, and effector function. T-cell activation and proliferation were reduced by inhibiting oxidative phosphorylation, whereas early cytokine production was maintained by either OXPHOS or glycolytic activity. Glucose deprivation in the presence of mild mitochondrial stress markedly reduced all three T-cell functions, contrasting the exposure to resveratrol, an antioxidant and sirtuin-1 activator, which specifically inhibited cytokine production and T-cell proliferation, but not T-cell activation. Conditions that inhibited T-cell activation were associated with the down-regulation of 2',5'-oligoadenylate synthetase genes via interferon response pathways. Our findings indicate that T-cell function is grossly impaired by stressors combined with nutrient deprivation, suggesting that correcting nutrient availability, metabolic stress, and/or the function of T cells in these conditions will improve the efficacy of T-cell-based therapies.

Convergent Losses of TLR5 Suggest Altered Extracellular Flagellin Detection in Four Mammalian Lineages.

Toll-like receptors (TLRs) play an important role for the innate immune system by detecting pathogen-associated molecular patterns. TLR5 encodes the major extracellular receptor for bacterial flagellin and frequently evolves under positive selection, consistent with coevolutionary arms races between the host and pathogens. Furthermore, TLR5 is inactivated in several vertebrates and a TLR5 stop codon polymorphism is widespread in human populations. Here, we analyzed the genomes of 120 mammals and discovered that TLR5 is convergently lost in four independent lineages, comprising guinea pigs, Yangtze river dolphin, pinnipeds, and pangolins. Validated inactivating mutations, absence of protein-coding transcript expression, and relaxed selection on the TLR5 remnants confirm these losses. PCR analysis further confirmed the loss of TLR5 in the pinniped stem lineage. Finally, we show that TLR11, encoding a second extracellular flagellin receptor, is also absent in these four lineages. Independent losses of TLR5 and TLR11 suggest that a major pathway for detecting flagellated bacteria is not essential for different mammals and predicts an impaired capacity to sense extracellular flagellin.

Losses of human disease-associated genes in placental mammals.

We systematically investigate whether losses of human disease-associated genes occurred in other mammals during evolution. We first show that genes lost in any of 62 non-human mammals generally have a lower degree of pleiotropy, and are highly depleted in essential and disease-associated genes. Despite this under-representation, we discovered multiple genes implicated in human disease that are truly lost in non-human mammals. In most cases, traits resembling human disease symptoms are present but not deleterious in gene-loss species, exemplified by losses of genes causing human eye or teeth disorders in poor-vision or enamel-less mammals. We also found widespread losses of PCSK9 and CETP genes, where loss-of-function mutations in humans protect from atherosclerosis. Unexpectedly, we discovered losses of disease genes (TYMP, TBX22, ABCG5, ABCG8, MEFV, CTSE) where deleterious phenotypes do not manifest in the respective species. A remarkable example is the uric acid-degrading enzyme UOX, which we found to be inactivated in elephants and manatees. While UOX loss in hominoids led to high serum uric acid levels and a predisposition for gout, elephants and manatees exhibit low uric acid levels, suggesting alternative ways of metabolizing uric acid. Together, our results highlight numerous mammals that are 'natural knockouts' of human disease genes.

Author Correction: A genomics approach reveals insights into the importance of gene losses for mammalian adaptations.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Gene Expression-Based Identification of Antigen-Responsive CD8+ T Cells on a Single-Cell Level.

CD8+ T cells are important effectors of adaptive immunity against pathogens, tumors, and self antigens. Here, we asked how human cognate antigen-responsive CD8+ T cells and their receptors could be identified in unselected single-cell gene expression data. Single-cell RNA sequencing and qPCR of dye-labeled antigen-specific cells identified large gene sets that were congruently up- or downregulated in virus-responsive CD8+ T cells under different antigen presentation conditions. Combined expression of TNFRSF9, XCL1, XCL2, and CRTAM was the most distinct marker of virus-responsive cells on a single-cell level. Using transcriptomic data, we developed a machine learning-based classifier that provides sensitive and specific detection of virus-responsive CD8+ T cells from unselected populations. Gene response profiles of CD8+ T cells specific for the autoantigen islet-specific glucose-6-phosphatase catalytic subunit-related protein differed markedly from virus-specific cells. These findings provide single-cell gene expression parameters for comprehensive identification of rare antigen-responsive cells and T cell receptors.

Genes lost during the transition from land to water in cetaceans highlight genomic changes associated with aquatic adaptations.

The transition from land to water in whales and dolphins (cetaceans) was accompanied by remarkable adaptations. To reveal genomic changes that occurred during this transition, we screened for protein-coding genes that were inactivated in the ancestral cetacean lineage. We found 85 gene losses. Some of these were likely beneficial for cetaceans, for example, by reducing the risk of thrombus formation during diving (F12 and KLKB1), erroneous DNA damage repair (POLM), and oxidative stress-induced lung inflammation (MAP3K19). Additional gene losses may reflect other diving-related adaptations, such as enhanced vasoconstriction during the diving response (mediated by SLC6A18) and altered pulmonary surfactant composition (SEC14L3), while loss of SLC4A9 relates to a reduced need for saliva. Last, loss of melatonin synthesis and receptor genes (AANAT, ASMT, and MTNR1A/B) may have been a precondition for adopting unihemispheric sleep. Our findings suggest that some genes lost in ancestral cetaceans were likely involved in adapting to a fully aquatic lifestyle.

Coding Exon-Structure Aware Realigner (CESAR): Utilizing Genome Alignments for Comparative Gene Annotation.

Alignment-based gene identification methods utilize sequence conservation between orthologous protein-coding genes to annotate genes in newly sequenced genomes. CESAR is an approach that makes use of existing genome alignments to transfer genes from one genome to other aligned genomes, and thus generates comparative gene annotations. To accurately detect conserved exons that exhibit an intact reading frame and consensus splice sites, CESAR produces a new alignment between orthologous exons, taking information about the exon's reading frame and splice site positions into account. Furthermore, CESAR is able to detect most evolutionary splice site shifts, which helps to annotate exon boundaries at high precision. Here, we describe how to apply CESAR to generate comparative gene annotations for one or many species, and discuss the strengths and limitations of this approach. CESAR is available at https://github.com/hillerlab/CESAR2.0 .

Retapamulin-Assisted Ribosome Profiling Reveals the Alternative Bacterial Proteome.

The use of alternative translation initiation sites enables production of more than one protein from a single gene, thereby expanding the cellular proteome. Although several such examples have been serendipitously found in bacteria, genome-wide mapping of alternative translation start sites has been unattainable. We found that the antibiotic retapamulin specifically arrests initiating ribosomes at start codons of the genes. Retapamulin-enhanced Ribo-seq analysis (Ribo-RET) not only allowed mapping of conventional initiation sites at the beginning of the genes, but strikingly, it also revealed putative internal start sites in a number of Escherichia coli genes. Experiments demonstrated that the internal start codons can be recognized by the ribosomes and direct translation initiation in vitro and in vivo. Proteins, whose synthesis is initiated at internal in-frame and out-of-frame start sites, can be functionally important and contribute to the "alternative" bacterial proteome. The internal start sites may also play regulatory roles in gene expression.

Convergent gene losses illuminate metabolic and physiological changes in herbivores and carnivores.

The repeated evolution of dietary specialization represents a hallmark of mammalian ecology. To detect genomic changes that are associated with dietary adaptations, we performed a systematic screen for convergent gene losses associated with an obligate herbivorous or carnivorous diet in 31 placental mammals. For herbivores, our screen discovered the repeated loss of the triglyceride lipase inhibitor PNLIPRP1, suggesting enhanced triglyceride digestion efficiency. Furthermore, several herbivores lost the pancreatic exocytosis factor SYCN, providing an explanation for continuous pancreatic zymogen secretion in these species. For carnivores, we discovered the repeated loss of the hormone-receptor pair INSL5-RXFP4 that regulates appetite and glucose homeostasis, which likely relates to irregular feeding patterns and constant gluconeogenesis. Furthermore, reflecting the reduced need to metabolize plant-derived xenobiotics, several carnivores lost the xenobiotic receptors NR1I3 and NR1I2 Finally, the carnivore-associated loss of the gastrointestinal host defense gene NOX1 could be related to a reduced gut microbiome diversity. By revealing convergent gene losses associated with differences in dietary composition, feeding patterns, and gut microbiomes, our study contributes to understanding how similar dietary specializations evolved repeatedly in mammals.

Loss of Enzymes in the Bile Acid Synthesis Pathway Explains Differences in Bile Composition among Mammals.

Bile acids are important for absorbing nutrients. Most mammals produce cholic and chenodeoxycholic bile acids. Here, we investigated genes in the bile acid synthesis pathway in four mammals that deviate from the usual mammalian bile composition. First, we show that naked-mole rats, elephants, and manatees repeatedly inactivated CYP8B1, an enzyme uniquely required for cholic acid synthesis, which explains the absence of cholic acid in these species. Second, no gene-inactivating mutations were found in any pathway gene in the rhinoceros, a species that lacks bile acids, indicating an evolutionarily recent change in its bile composition. Third, elephants and/or manatees that also lack bile acids altogether have lost additional nonessential enzymes (SLC27A5, ACOX2). Apart from uncovering genomic differences explaining deviations in bile composition, our analysis of bile acid enzymes in bile acid-lacking species suggests that essentiality prevents gene loss, while loss of pleiotropic genes is permitted if their other functions are compensated by functionally related proteins.

CD4+ T cell activation, function, and metabolism are inhibited by low concentrations of DMSO.

Dimethyl sulfoxide (DMSO) is a polar organic solvent used in a wide range of biological applications. DMSO is routinely used as a cryoprotectant for long-term cell freezing as well as to dissolve peptides and drugs for immune cell functional assays. Here, human CD4+ T cell activation, cytokine production, proliferation, and metabolism were investigated after stimulation in the presence of 0.01% to 1%, DMSO, representing concentrations commonly used in vitro. Surface expression of the activation markers CD69, CD25 and CD154 after polyclonal activation of CD4+ T cells was inhibited by 0.25% or higher concentrations of DMSO. The frequencies of IL-21+, IL-4+, and IL-22+ CD4+ T cells, following polyclonal activation were variably inhibited by DMSO at concentrations ranging from 0.25% to 1%, whereas IFNγ+ cells were unaffected. CD4+ T cell proliferation after anti-CD3 or antigen stimulation was inhibited by 0.5% DMSO and abolished by 1% DMSO. After polyclonal stimulation, glucose uptake was inhibited in the presence of 1% DMSO, but only minor effects on CD4+ T cell respiration were observed. Consistent with the immune effects, the gene expression of early signaling and activation pathways were inhibited in CD4+ T cells in the presence of 1% DMSO. Our study revealed that DMSO at concentrations generally used for in vitro studies of T cells impacts multiple features of T cell function. Therefore, we urge care when adding DMSO-containing preparations to T cell cultures.

Loss of RXFP2 and INSL3 genes in Afrotheria shows that testicular descent is the ancestral condition in placental mammals.

Descent of testes from a position near the kidneys into the lower abdomen or into the scrotum is an important developmental process that occurs in all placental mammals, with the exception of five afrotherian lineages. Since soft-tissue structures like testes are not preserved in the fossil record and since key parts of the placental mammal phylogeny remain controversial, it has been debated whether testicular descent is the ancestral or derived condition in placental mammals. To resolve this debate, we used genomic data of 71 mammalian species and analyzed the evolution of two key genes (relaxin/insulin-like family peptide receptor 2 [RXFP2] and insulin-like 3 [INSL3]) that induce the development of the gubernaculum, the ligament that is crucial for testicular descent. We show that both RXFP2 and INSL3 are lost or nonfunctional exclusively in four afrotherians (tenrec, cape elephant shrew, cape golden mole, and manatee) that completely lack testicular descent. The presence of remnants of once functional orthologs of both genes in these afrotherian species shows that these gene losses happened after the split from the placental mammal ancestor. These "molecular vestiges" provide strong evidence that testicular descent is the ancestral condition, irrespective of persisting phylogenetic discrepancies. Furthermore, the absence of shared gene-inactivating mutations and our estimates that the loss of RXFP2 happened at different time points strongly suggest that testicular descent was lost independently in Afrotheria. Our results provide a molecular mechanism that explains the loss of testicular descent in afrotherians and, more generally, highlight how molecular vestiges can provide insights into the evolution of soft-tissue characters.

A genomics approach reveals insights into the importance of gene losses for mammalian adaptations.

Identifying the genomic changes that underlie phenotypic adaptations is a key challenge in evolutionary biology and genomics. Loss of protein-coding genes is one type of genomic change with the potential to affect phenotypic evolution. Here, we develop a genomics approach to accurately detect gene losses and investigate their importance for adaptive evolution in mammals. We discover a number of gene losses that likely contributed to morphological, physiological, and metabolic adaptations in aquatic and flying mammals. These gene losses shed light on possible molecular and cellular mechanisms that underlie these adaptive phenotypes. In addition, we show that gene loss events that occur as a consequence of relaxed selection following adaptation provide novel insights into species' biology. Our results suggest that gene loss is an evolutionary mechanism for adaptation that may be more widespread than previously anticipated. Hence, investigating gene losses has great potential to reveal the genomic basis underlying macroevolutionary changes.

Transition to an Aquatic Habitat Permitted the Repeated Loss of the Pleiotropic KLK8 Gene in Mammals.

Kallikrein related peptidase 8 (KLK8; also called neuropsin) is a serine protease that plays distinct roles in the skin and hippocampus. In the skin, KLK8 influences keratinocyte proliferation and desquamation, and activates antimicrobial peptides in sweat. In the hippocampus, KLK8 affects memory acquisition. Here, we examined the evolution of KLK8 in mammals and discovered that, out of 70 placental mammals, KLK8 is exclusively lost in three independent fully-aquatic lineages, comprising dolphin, killer whale, minke whale, and manatee. In addition, while the sperm whale has an intact KLK8 reading frame, the gene evolves neutrally in this species. We suggest that the distinct functions of KLK8 likely became obsolete in the aquatic environment, leading to the subsequent loss of KLK8 in several fully-aquatic mammalian lineages. First, the cetacean and manatee skin lacks sweat glands as an adaptation to the aquatic environment, which likely made the epidermal function of KLK8 obsolete. Second, cetaceans and manatees exhibit a proportionally small hippocampus, which may have rendered the hippocampal functions of KLK8 obsolete. Together, our results shed light on the genomic changes that correlate with skin and neuroanatomical differences of aquatic mammals, and show that even pleiotropic genes can be lost during evolution if an environmental change nullifies the need for the different functions of such genes.

CESAR 2.0 substantially improves speed and accuracy of comparative gene annotation.

MOTIVATION: Homology-based gene prediction is a powerful concept to annotate newly sequenced genomes. We have previously demonstrated that whole genome alignments can be utilized for accurate comparative coding gene annotation. RESULTS: Here we present CESAR 2.0 that utilizes genome alignments to transfer coding gene annotations from one reference to many other aligned genomes. We show that CESAR 2.0 is 77 times faster and requires 31 times less memory compared to its predecessor. CESAR 2.0 substantially improves the ability to align splice sites that have shifted over larger distances, allowing for precise identification of the exon boundaries in the aligned genome. Finally, CESAR 2.0 supports entire genes, which enables the annotation of joined exons that arose by complete intron deletions. CESAR 2.0 can readily be applied to new genome alignments to annotate coding genes in many other genomes at improved accuracy and without necessitating large-computational resources. AVAILABILITY AND IMPLEMENTATION: Source code is freely available at https://github.com/hillerlab/CESAR2.0. CONTACT: hiller@mpi-cbg.de. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Increased alignment sensitivity improves the usage of genome alignments for comparative gene annotation.

Genome alignments provide a powerful basis to transfer gene annotations from a well-annotated reference genome to many other aligned genomes. The completeness of these annotations crucially depends on the sensitivity of the underlying genome alignment. Here, we investigated the impact of the genome alignment parameters and found that parameters with a higher sensitivity allow the detection of thousands of novel alignments between orthologous exons that have been missed before. In particular, comparisons between species separated by an evolutionary distance of >0.75 substitutions per neutral site, like human and other non-placental vertebrates, benefit from increased sensitivity. To systematically test if increased sensitivity improves comparative gene annotations, we built a multiple alignment of 144 vertebrate genomes and used this alignment to map human genes to the other 143 vertebrates with CESAR. We found that higher alignment sensitivity substantially improves the completeness of comparative gene annotations by adding on average 2382 and 7440 novel exons and 117 and 317 novel genes for mammalian and non-mammalian species, respectively. Our results suggest a more sensitive alignment strategy that should generally be used for genome alignments between distantly-related species. Our 144-vertebrate genome alignment and the comparative gene annotations (https://bds.mpi-cbg.de/hillerlab/144VertebrateAlignment_CESAR/) are a valuable resource for comparative genomics.