Evidence for extensive horizontal gene transfer from the draft genome of a tardigrade

Antioxidant Defense in the Toughest Animals on the Earth: Its Contribution to the Extreme Resistance of Tardigrades

Tardigrades are unique among animals in their resistance to dehydration, mainly due to anhydrobiosis and tun formation. They are also very resistant to high-energy radiation, low and high temperatures, low and high pressure, and various chemical agents, Interestingly, they are resistant to ionizing radiation both in the hydrated and dehydrated states to a similar extent. They are able to survive in the cosmic space. Apparently, many mechanisms contribute to the resistance of tardigrades to harmful factors, including the presence of trehalose (though not common to all tardigrades), heat shock proteins, late embryogenesis-abundant proteins, tardigrade-unique proteins, DNA repair proteins, proteins directly protecting DNA (Dsup and TDR1), and efficient antioxidant system. Antioxidant enzymes and small-molecular-weight antioxidants are an important element in the tardigrade resistance. The levels and activities of many antioxidant proteins is elevated by anhydrobiosis and UV radiation; one explanation for their induction during dehydration is provided by the theory of "preparation for oxidative stress", which occurs during rehydration. Genes coding for some antioxidant proteins are expanded in tardigrades; some genes (especially those coding for catalases) were hypothesized to be of bacterial origin, acquired by horizontal gene transfer. An interesting antioxidant protein found in tardigrades is the new Mn-dependent peroxidase.

How to Survive without Water: A Short Lesson on the Desiccation Tolerance of Budding Yeast

Water is essential to all life on earth. It is a major component that makes up living organisms and plays a vital role in multiple biological processes. It provides a medium for chemical and enzymatic reactions in the cell and is a major player in osmoregulation and the maintenance of cell turgidity. Despite this, many organisms, called anhydrobiotes, are capable of surviving under extremely dehydrated conditions. Less is known about how anhydrobiotes adapt and survive under desiccation stress. Studies have shown that morphological and physiological changes occur in anhydrobiotes in response to desiccation stress. Certain disaccharides and proteins, including heat shock proteins, intrinsically disordered proteins, and hydrophilins, play important roles in the desiccation tolerance of anhydrobiotes. In this review, we summarize the recent findings of desiccation tolerance in the budding yeast Saccharomyces cerevisiae. We also propose that the yeast under desiccation could be used as a model to study neurodegenerative disorders.

Crossroads of assembling a moss genome: navigating contaminants and horizontal gene transfer in the moss Physcomitrellopsis africana

The first chromosome-scale reference genome of the rare narrow-endemic African moss Physcomitrellopsis africana (P. africana) is presented here. Assembled from 73 × Oxford Nanopore Technologies (ONT) long reads and 163 × Beijing Genomics Institute (BGI)-seq short reads, the 414 Mb reference comprises 26 chromosomes and 22,925 protein-coding genes [Benchmarking Universal Single-Copy Ortholog (BUSCO) scores: C:94.8% (D:13.9%)]. This genome holds 2 genes that withstood rigorous filtration of microbial contaminants, have no homolog in other land plants, and are thus interpreted as resulting from 2 unique horizontal gene transfers (HGTs) from microbes. Further, P. africana shares 176 of the 273 published HGT candidates identified in Physcomitrium patens (P. patens), but lacks 98 of these, highlighting that perhaps as many as 91 genes were acquired in P. patens in the last 40 million years following its divergence from its common ancestor with P. africana. These observations suggest rather continuous gene gains via HGT followed by potential losses during the diversification of the Funariaceae. Our findings showcase both dynamic flux in plant HGTs over evolutionarily "short" timescales, alongside enduring impacts of successful integrations, like those still functionally maintained in extant P. africana. Furthermore, this study describes the informatic processes employed to distinguish contaminants from candidate HGT events.

Are Viruses Taxonomic Units? A Protein Domain and Loop-Centric Phylogenomic Assessment

Virus taxonomy uses a Linnaean-like subsumption hierarchy to classify viruses into taxonomic units at species and higher rank levels. Virus species are considered monophyletic groups of mobile genetic elements (MGEs) often delimited by the phylogenetic analysis of aligned genomic or metagenomic sequences. Taxonomic units are assumed to be independent organizational, functional and evolutionary units that follow a 'natural history' rationale. Here, I use phylogenomic and other arguments to show that viruses are not self-standing genetically-driven systems acting as evolutionary units. Instead, they are crucial components of holobionts, which are units of biological organization that dynamically integrate the genetics, epigenetic, physiological and functional properties of their co-evolving members. Remarkably, phylogenomic analyses show that viruses share protein domains and loops with cells throughout history via massive processes of reticulate evolution, helping spread evolutionary innovations across a wider taxonomic spectrum. Thus, viruses are not merely MGEs or microbes. Instead, their genomes and proteomes conduct cellularly integrated processes akin to those cataloged by the GO Consortium. This prompts the generation of compositional hierarchies that replace the 'is-a-kind-of' by a 'is-a-part-of' logic to better describe the mereology of integrated cellular and viral makeup. My analysis demands a new paradigm that integrates virus taxonomy into a modern evolutionarily centered taxonomy of organisms.

Pinpointing the microbiota of tardigrades: What is really there?

Microbiota are considered significant in the biology of tardigrades, yet their diversity and distribution remain largely unexplored. This is partly due to the methodological challenges associated with studying the microbiota of small organisms that inhabit microbe-rich environments. In our study, we characterized the microbiota of 31 species of cultured tardigrades using 16S rRNA amplicon sequencing. We employed various sample preparation strategies and multiple types of controls and estimated the number of microbes in samples using synthetic DNA spike-ins. We also reanalysed data from previous tardigrade microbiome studies. Our findings suggest that the microbial communities of cultured tardigrades are predominantly composed of bacterial genotypes originating from food, medium, or reagents. Despite numerous experiments, we found it challenging to identify strains that were enriched in certain tardigrades, which would have indicated likely symbiotic associations. Putative tardigrade-associated microbes rarely constituted more than 20% of the datasets, although some matched symbionts identified in other studies. We also uncovered serious contamination issues in previous tardigrade microbiome studies, casting doubt on some of their conclusions. We concluded that tardigrades are not universally dependent on specialized microbes. Our work underscores the need for rigorous safeguards in studies of the microbiota of microscopic organisms and serves as a cautionary tale for studies involving samples with low microbiome abundance.

Horizontal gene transfer in eukaryotes: aligning theory with data

Horizontal gene transfer (HGT), or lateral gene transfer, is the non-sexual movement of genetic information between genomes. It has played a pronounced part in bacterial and archaeal evolution, but its role in eukaryotes is less clear. Behaviours unique to eukaryotic cells - phagocytosis and endosymbiosis - have been proposed to increase the frequency of HGT, but nuclear genomes encode fewer HGTs than bacteria and archaea. Here, I review the existing theory in the context of the growing body of data on HGT in eukaryotes, which suggests that any increased chance of acquiring new genes through phagocytosis and endosymbiosis is offset by a reduced need for these genes in eukaryotes, because selection in most eukaryotes operates on variation not readily generated by HGT.

Anaerobic fungi in the tortoise alimentary tract illuminate early stages of host-fungal symbiosis and Neocallimastigomycota evolution

Anaerobic gut fungi (AGF, Neocallimastigomycota) reside in the alimentary tract of herbivores. While their presence in mammals is well documented, evidence for their occurrence in non-mammalian hosts is currently sparse. Culture-independent surveys of AGF in tortoises identified a unique community, with three novel deep-branching genera representing >90% of sequences in most samples. Representatives of all genera were successfully isolated under strict anaerobic conditions. Transcriptomics-enabled phylogenomic and molecular dating analyses indicated an ancient, deep-branching position in the AGF tree for these genera, with an evolutionary divergence time estimate of 104-112 million years ago (Mya). Such estimates push the establishment of animal-Neocallimastigomycota symbiosis from the late to the early Cretaceous. Further, tortoise-associated isolates (T-AGF) exhibited limited capacity for plant polysaccharides metabolism and lacked genes encoding several carbohydrate-active enzyme (CAZyme) families. Finally, we demonstrate that the observed curtailed degradation capacities and reduced CAZyme repertoire is driven by the paucity of horizontal gene transfer (HGT) in T-AGF genomes, compared to their mammalian counterparts. This reduced capacity was reflected in an altered cellulosomal production capacity in T-AGF. Our findings provide insights into the phylogenetic diversity, ecological distribution, evolutionary history, evolution of fungal-host nutritional symbiosis, and dynamics of genes acquisition in Neocallimastigomycota.

The genome of the invasive and broadly polyphagous Diaprepes root weevil, Diaprepes abbreviatus (Coleoptera), reveals an arsenal of putative polysaccharide-degrading enzymes

The Diaprepes root weevil (DRW), Diaprepes abbreviatus, is a broadly polyphagous invasive pest of agriculture in the southern United States and the Caribbean. Its genome was sequenced, assembled, and annotated to study genomic correlates of specialized plant-feeding and invasiveness and to facilitate the development of new methods for DRW control. The 1.69 Gb D. abbreviatus genome assembly was distributed across 653 contigs, with an N50 of 7.8 Mb and the largest contig of 62 Mb. Most of the genome was comprised of repetitive sequences, with 66.17% in transposable elements, 5.75% in macrosatellites, and 2.06% in microsatellites. Most expected orthologous genes were present and fully assembled, with 99.5% of BUSCO genes present and 1.5% duplicated. One hundred and nine contigs (27.19 Mb) were identified as putative fragments of the X and Y sex chromosomes, and homology assessment with other beetle X chromosomes indicated a possible sex chromosome turnover event. Genome annotation identified 18,412 genes, including 43 putative horizontally transferred (HT) loci. Notably, 258 genes were identified from gene families known to encode plant cell wall degrading enzymes and invertases, including carbohydrate esterases, polysaccharide lyases, and glycoside hydrolases (GH). GH genes were unusually numerous, with 239 putative genes representing 19 GH families. Interestingly, several other beetle species with large numbers of GH genes are (like D. abbreviatus) successful invasive pests of agriculture or forestry.

ContScout: sensitive detection and removal of contamination from annotated genomes

Contamination of genomes is an increasingly recognized problem affecting several downstream applications, from comparative evolutionary genomics to metagenomics. Here we introduce ContScout, a precise tool for eliminating foreign sequences from annotated genomes. It achieves high specificity and sensitivity on synthetic benchmark data even when the contaminant is a closely related species, outperforms competing tools, and can distinguish horizontal gene transfer from contamination. A screen of 844 eukaryotic genomes for contamination identified bacteria as the most common source, followed by fungi and plants. Furthermore, we show that contaminants in ancestral genome reconstructions lead to erroneous early origins of genes and inflate gene loss rates, leading to a false notion of complex ancestral genomes. Taken together, we offer here a tool for sensitive removal of foreign proteins, identify and remove contaminants from diverse eukaryotic genomes and evaluate their impact on phylogenomic analyses.

GTax: improving de novo transcriptome assembly by removing foreign RNA contamination

The cost and complexity of generating a complete reference genome means that many organisms lack an annotated reference. An alternative is to use a de novo reference transcriptome. This technology is cost-effective but is susceptible to off-target RNA contamination. In this manuscript, we present GTax, a taxonomy-structured database of genomic sequences that can be used with BLAST to detect and remove foreign contamination in RNA sequencing samples before assembly. In addition, we use a de novo transcriptome assembly of Solanum lycopersicum (tomato) to demonstrate that removing foreign contamination in sequencing samples reduces the number of assembled chimeric transcripts.

Population Genomics of Pooled Samples: Unveiling Symbiont Infrapopulation Diversity and Host-Symbiont Coevolution

Microscopic symbionts represent crucial links in biological communities. However, they present technical challenges in high-throughput sequencing (HTS) studies due to their small size and minimal high-quality DNA yields, hindering our understanding of host-symbiont coevolution at microevolutionary and macroevolutionary scales. One approach to overcome those barriers is to pool multiple individuals from the same infrapopulation (i.e., individual host) and sequence them together (Pool-Seq), but individual-level information is then compromised. To simultaneously address both issues (i.e., minimal DNA yields and loss of individual-level information), we implemented a strategic Pool-Seq approach to assess variation in sequencing performance and categorize genetic diversity (single nucleotide polymorphisms (SNPs)) at both the individual-level and infrapopulation-level for microscopic feather mites. To do so, we collected feathers harboring mites (Proctophyllodidae: Amerodectes protonotaria) from four individual Prothonotary Warblers (Parulidae: Protonotaria citrea). From each of the four hosts (i.e., four mite infrapopulations), we conducted whole-genome sequencing on three extraction pools consisting of different numbers of mites (1 mite, 5 mites, and 20 mites). We found that samples containing pools of multiple mites had more sequencing reads map to the feather mite reference genome than did the samples containing only a single mite. Mite infrapopulations were primarily genetically structured by their associated individual hosts (not pool size) and the majority of SNPs were shared by all pools within an infrapopulation. Together, these results suggest that the patterns observed are driven by evolutionary processes occurring at the infrapopulation level and are not technical signals due to pool size. In total, despite the challenges presented by microscopic symbionts in HTS studies, this work highlights the value of both individual-level and infrapopulation-level sequencing toward our understanding of host-symbiont coevolution at multiple evolutionary scales.

HELIOS: High-speed sequence alignment in optics

In response to the imperfections of current sequence alignment methods, originated from the inherent serialism within their corresponding electrical systems, a few optical approaches for biological data comparison have been proposed recently. However, due to their low performance, raised from their inefficient coding scheme, this paper presents a novel all-optical high-throughput method for aligning DNA, RNA, and protein sequences, named HELIOS. The HELIOS method employs highly sophisticated operations to locate character matches, single or multiple mutations, and single or multiple indels within various biological sequences. On the other hand, the HELIOS optical architecture exploits high-speed processing and operational parallelism in optics, by adopting wavelength and polarization of optical beams. For evaluation, the functionality and accuracy of the HELIOS method are approved through behavioral and optical simulation studies, while its complexity and performance are estimated through analytical computation. The accuracy evaluations indicate that the HELIOS method achieves a precise pairwise alignment of two sequences, highly similar to those of Smith-Waterman, Needleman-Wunsch, BLAST, MUSCLE, ClustalW, ClustalΩ, T-Coffee, Kalign, and MAFFT. According to our performance evaluations, the HELIOS optical architecture outperforms all alternative electrical and optical algorithms in terms of processing time and memory requirement, relying on its highly sophisticated method and optical architecture. Moreover, the employed compact coding scheme highly escalates the number of input characters, and hence, it offers reduced time and space complexities, compared to the electrical and optical alternatives. It makes the HELIOS method and optical architecture highly applicable for biomedical applications.

Propagation, detection and correction of errors using the sequence database network

Nucleotide and protein sequences stored in public databases are the cornerstone of many bioinformatics analyses. The records containing these sequences are prone to a wide range of errors, including incorrect functional annotation, sequence contamination and taxonomic misclassification. One source of information that can help to detect errors are the strong interdependency between records. Novel sequences in one database draw their annotations from existing records, may generate new records in multiple other locations and will have varying degrees of similarity with existing records across a range of attributes. A network perspective of these relationships between sequence records, within and across databases, offers new opportunities to detect-or even correct-erroneous entries and more broadly to make inferences about record quality. Here, we describe this novel perspective of sequence database records as a rich network, which we call the sequence database network, and illustrate the opportunities this perspective offers for quantification of database quality and detection of spurious entries. We provide an overview of the relevant databases and describe how the interdependencies between sequence records across these databases can be exploited by network analyses. We review the process of sequence annotation and provide a classification of sources of error, highlighting propagation as a major source. We illustrate the value of a network perspective through three case studies that use network analysis to detect errors, and explore the quality and quantity of critical relationships that would inform such network analyses. This systematic description of a network perspective of sequence database records provides a novel direction to combat the proliferation of errors within these critical bioinformatics resources.

Trehalose and tardigrade CAHS proteins work synergistically to promote desiccation tolerance

Tardigrades are microscopic animals renowned for their ability to survive extreme desiccation. Unlike many desiccation-tolerant organisms that accumulate high levels of the disaccharide trehalose to protect themselves during drying, tardigrades accumulate little or undetectable levels. Using comparative metabolomics, we find that despite being enriched at low levels, trehalose is a key biomarker distinguishing hydration states of tardigrades. In vitro, naturally occurring stoichiometries of trehalose and CAHS proteins, intrinsically disordered proteins with known protective capabilities, were found to produce synergistic protective effects during desiccation. In vivo, this synergistic interaction is required for robust CAHS-mediated protection. This demonstrates that trehalose acts not only as a protectant, but also as a synergistic cosolute. Beyond desiccation tolerance, our study provides insights into how the solution environment tunes intrinsically disordered proteins’ functions, many of which are vital in biological contexts such as development and disease that are concomitant with large changes in intracellular chemistry.

The disaccharide trehalose is a synergistic cosolute and key biomarker of desiccation tolerance in tardigrades.

Computational analysis on two putative mitochondrial protein-coding genes from the Emydura subglobosa genome: A functional annotation approach

Rapid advancements in automated genomic technologies have uncovered many unique findings about the turtle genome and its associated features including olfactory gene expansions and duplications of toll-like receptors. However, despite the advent of large-scale sequencing, assembly, and annotation, about 40-50% of genes in eukaryotic genomes are left without functional annotation, severely limiting our knowledge of the biological information of genes. Additionally, these automated processes are prone to errors since draft genomes consist of several disconnected scaffolds whose order is unknown; erroneous draft assemblies may also be contaminated with foreign sequences and propagate to cause errors in annotation. Many of these automated annotations are thus incomplete and inaccurate, highlighting the need for functional annotation to link gene sequences to biological identity. In this study, we have functionally annotated two genes of the red-bellied short-neck turtle (Emydura subglobosa), a member of the relatively understudied pleurodire lineage of turtles. We improved upon initial ab initio gene predictions through homology-based evidence and generated refined consensus gene models. Through functional, localization, and structural analyses of the predicted proteins, we discovered conserved putative genes encoding mitochondrial proteins that play a role in C21-steroid hormone biosynthetic processes and fatty acid catabolism-both of which are distantly related by the tricarboxylic acid (TCA) cycle and share similar metabolic pathways. Overall, these findings further our knowledge about the genetic features underlying turtle physiology, morphology, and longevity, which have important implications for the treatment of human diseases and evolutionary studies.

Deciphering the Biological Enigma-Genomic Evolution Underlying Anhydrobiosis in the Phylum Tardigrada and the Chironomid Polypedilum vanderplanki

Anhydrobiosis, an ametabolic dehydrated state triggered by water loss, is observed in several invertebrate lineages. Anhydrobiotes revive when rehydrated, and seem not to suffer the ultimately lethal cell damage that results from severe loss of water in other organisms. Here, we review the biochemical and genomic evidence that has revealed the protectant molecules, repair systems, and maintenance pathways associated with anhydrobiosis. We then introduce two lineages in which anhydrobiosis has evolved independently: Tardigrada, where anhydrobiosis characterizes many species within the phylum, and the genus Polypedilum, where anhydrobiosis occurs in only two species. Finally, we discuss the complexity of the evolution of anhydrobiosis within invertebrates based on current knowledge, and propose perspectives to enhance the understanding of anhydrobiosis.

Tardigrades and their emergence as model organisms

Experimentally tractable organisms like C. elegans, Drosophila, zebrafish, and mouse are popular models for addressing diverse questions in biology. In 1997, two of the most valuable invertebrate model organisms to date-C. elegans and Drosophila-were found to be much more closely related to each other than expected. C. elegans and Drosophila belong to the nematodes and arthropods, respectively, and these two phyla and six other phyla make up a clade of molting animals referred to as the Ecdysozoa. The other ecdysozoan phyla could be valuable models for comparative biology, taking advantage of the rich and continual sources of research findings as well as tools from both C. elegans and Drosophila. But when the Ecdysozoa was first recognized, few tools were available for laboratory studies in any of these six other ecdysozoan phyla. In 1999 I began an effort to develop tools for studying one such phylum, the tardigrades. Here, I describe how the tardigrade species Hypsibius exemplaris and tardigrades more generally have emerged over the past two decades as valuable new models for answering diverse questions. To date, these questions have included how animal body plans evolve and how biological materials can survive some remarkably extreme conditions.

Contamination detection in genomic data: more is not enough

The decreasing cost of sequencing and concomitant augmentation of publicly available genomes have created an acute need for automated software to assess genomic contamination. During the last 6 years, 18 programs have been published, each with its own strengths and weaknesses. Deciding which tools to use becomes more and more difficult without an understanding of the underlying algorithms. We review these programs, benchmarking six of them, and present their main operating principles. This article is intended to guide researchers in the selection of appropriate tools for specific applications. Finally, we present future challenges in the developing field of contamination detection.

Examples of Extreme Survival: Tardigrade Genomics and Molecular Anhydrobiology

Tardigrades are ubiquitous meiofauna that are especially renowned for their exceptional extremotolerance to various adverse environments, including pressure, temperature, and even ionizing radiation. This is achieved through a reversible halt of metabolism triggered by desiccation, a phenomenon called anhydrobiosis. Recent establishment of genome resources for two tardigrades, Hypsibius exemplaris and Ramazzottius varieornatus, accelerated research to uncover the molecular mechanisms behind anhydrobiosis, leading to the discovery of many tardigrade-unique proteins. This review focuses on the history, methods, discoveries, and current state and challenges regarding tardigrade genomics, with an emphasis on molecular anhydrobiology. Remaining questions and future perspectives regarding prospective approaches to fully elucidate the molecular machinery of this complex phenomenon are discussed.

Accurate annotation of protein coding sequences with IDTAXA

The observed diversity of protein coding sequences continues to increase far more rapidly than knowledge of their functions, making classification algorithms essential for assigning a function to proteins using only their sequence. Most pipelines for annotating proteins rely on searches for homologous sequences in databases of previously annotated proteins using BLAST or HMMER. Here, we develop a new approach for classifying proteins into a taxonomy of functions and demonstrate its utility for genome annotation. Our algorithm, IDTAXA, was more accurate than BLAST or HMMER at assigning sequences to KEGG ortholog groups. Moreover, IDTAXA correctly avoided classifying sequences with novel functions to existing groups, which is a common error mode for classification approaches that rely on E-values as a proxy for confidence. We demonstrate IDTAXA's utility for annotating eukaryotic and prokaryotic genomes by assigning functions to proteins within a multi-level ontology and applied IDTAXA to detect genome contamination in eukaryotic genomes. Finally, we re-annotated 8604 microbial genomes with known antibiotic resistance phenotypes to discover two novel associations between proteins and antibiotic resistance. IDTAXA is available as a web tool (http://DECIPHER.codes/Classification.html) or as part of the open source DECIPHER R package from Bioconductor.

Endurance of extremely prolonged nutrient prevention across kingdoms of life

Numerous observations demonstrate that microorganisms can survive very long periods of nutrient deprivation and starvation. Moreover, it is evident that prolonged periods of starvation are a feature of many habitats, and many cells in all kingdoms of life are found in prolonged starvation conditions. Bacteria exhibit a range of responses to long-term starvation. These include genetic adaptations such as the long-term stationary phase and the growth advantage in stationary phase phenotypes characterized by mutations in stress-signaling genes and elevated mutation rates. Here, we suggest using the term "endurance of prolonged nutrient prevention" (EPNP phase), to describe this phase, which was also recently described in eukaryotes. Here, we review this literature and describe the current knowledge about the adaptations to very long-term starvation conditions in bacteria and eukaryotes, its conceptual and structural conservation across all kingdoms of life, and point out possible directions that merit further research.

Fields, biochemistry fast autooxidation of a Bis-Histidyl-ligated globin from the anhydrobiotic tardigrade, ramazzottius varieornatus, by molecular oxygen

Tardigrades, a phylum of meiofaunal organisms, exhibit extraordinary tolerance to various environmental conditions, including extreme temperatures (-272 to 151 °C) and exposure to ionizing radiation. Proteins from anhydrobiotic tardigrades with homology to known proteins from other organisms are new potential targets for structural genomics. Recently, we reported spectroscopic and structural characterization of a hexacoordinated hemoglobin (Kumaglobin [Kgb]) found in an anhydrobiotic tardigrade. In the absence of its exogenous ligand, Kgb displays hexacoordination with distal and proximal histidines. In this work, we analyzed binding of the molecular oxygen ligand following reduction of heme in Kgb using a pulse radiolysis technique. Radiolytically generated hydrated electrons (eaq-) reduced the heme iron of Kgb within 20 µs. Subsequently, ferrous heme reacted with O2 to form a ferrous-dioxygen intermediate with a second-order rate constant of 3.0 × 106 M-1 s-1. The intermediate was rapidly (within 0.1 s) autooxidized to the ferric form. Redox potential measurements revealed an E'0 of -400 mV (vs. SHE) in the ferric/ferrous couple. Our results suggest that Kgb may serve as a physiological generator of O2·- via redox signaling and/or electron transfer.

Streamlining data-intensive biology with workflow systems

As the scale of biological data generation has increased, the bottleneck of research has shifted from data generation to analysis. Researchers commonly need to build computational workflows that include multiple analytic tools and require incremental development as experimental insights demand tool and parameter modifications. These workflows can produce hundreds to thousands of intermediate files and results that must be integrated for biological insight. Data-centric workflow systems that internally manage computational resources, software, and conditional execution of analysis steps are reshaping the landscape of biological data analysis and empowering researchers to conduct reproducible analyses at scale. Adoption of these tools can facilitate and expedite robust data analysis, but knowledge of these techniques is still lacking. Here, we provide a series of strategies for leveraging workflow systems with structured project, data, and resource management to streamline large-scale biological analysis. We present these practices in the context of high-throughput sequencing data analysis, but the principles are broadly applicable to biologists working beyond this field.

Conservative route to genome compaction in a miniature annelid

The causes and consequences of genome reduction in animals are unclear because our understanding of this process mostly relies on lineages with often exceptionally high rates of evolution. Here, we decode the compact 73.8-megabase genome of Dimorphilus gyrociliatus, a meiobenthic segmented worm. The D. gyrociliatus genome retains traits classically associated with larger and slower-evolving genomes, such as an ordered, intact Hox cluster, a generally conserved developmental toolkit and traces of ancestral bilaterian linkage. Unlike some other animals with small genomes, the analysis of the D. gyrociliatus epigenome revealed canonical features of genome regulation, excluding the presence of operons and trans-splicing. Instead, the gene-dense D. gyrociliatus genome presents a divergent Myc pathway, a key physiological regulator of growth, proliferation and genome stability in animals. Altogether, our results uncover a conservative route to genome compaction in annelids, reminiscent of that observed in the vertebrate Takifugu rubripes.

This study reports the genome of the miniature segmented annelid Dimorphilus gyrociliatus and reveals no drastic changes in genome architecture and regulation, unlike other cases of genome miniaturization.

The biology of tardigrade disordered proteins in extreme stress tolerance

Abstract

Disordered proteins have long been known to help mediate tolerance to different abiotic stresses including freezing, osmotic stress, high temperatures, and desiccation in a diverse set of organisms. Recently, three novel families of intrinsically disordered proteins were identified in tardigrades, microscopic animals capable of surviving a battery of environmental extremes. These three families include the Cytoplasmic-, Secreted-, and Mitochondrial- Abundant Heat Soluble (CAHS, SAHS, and MAHS) proteins, which are collectively termed Tardigrade Disordered Proteins (TDPs). At the level of sequence conservation TDPs are unique to tardigrades, and beyond their high degree of disorder the CAHS, SAHS, and MAHS families do not resemble one another. All three families are either highly expressed constitutively, or significantly enriched in response to desiccation. In vivo, ex vivo, and in vitro experiments indicate functional roles for members of each TDP family in mitigating cellular perturbations induced by various abiotic stresses. What is currently lacking is a comprehensive and holistic understanding of the fundamental mechanisms by which TDPs function, and the properties of TDPs that allow them to function via those mechanisms. A quantitative and systematic approach is needed to identify precisely what cellular damage TDPs work to prevent, what sequence features are important for these functions, and how those sequence features contribute to the underlying mechanisms of protection. Such an approach will inform us not only about these fascinating proteins, but will also provide insights into how the sequence of a disordered protein can dictate its functional, structural, and dynamic properties.

A Chromosome-Level Genome Assembly of Garlic (Allium sativum) Provides Insights into Genome Evolution and Allicin Biosynthesis

Garlic, an economically important vegetable, spice, and medicinal crop, produces highly enlarged bulbs and unique organosulfur compounds. Here, we report a chromosome-level genome assembly for garlic, with a total size of approximately 16.24 Gb, as well as the annotation of 57 561 predicted protein-coding genes, making garlic the first Allium species with a sequenced genome. Analysis of this garlic genome assembly reveals a recent burst of transposable elements, explaining the substantial expansion of the garlic genome. We examined the evolution of certain genes associated with the biosynthesis of allicin and inulin neoseries-type fructans, and provided new insights into the biosynthesis of these two compounds. Furthermore, a large-scale transcriptome was produced to characterize the expression patterns of garlic genes in different tissues and at various growth stages of enlarged bulbs. The reference genome and large-scale transcriptome data generated in this study provide valuable new resources for research on garlic biology and breeding.

Loss of intermediate regions of perpendicular body axes contributed to miniaturization of tardigrades

Tardigrades have a miniaturized body plan. Miniaturization in tardigrades is associated with the loss of several organ systems and an intermediate region of their anteroposterior (AP) axis. However, how miniaturization has affected tardigrade legs is unclear. In arthropods and in onychophorans, the leg gap genes are expressed in regionalized proximodistal (PD) patterns in the legs. Functional studies indicate that these genes regulate growth in their respective expression domains and establish PD identities, partly through mutually antagonistic regulatory interactions. Here, we investigated the expression patterns of tardigrade orthologs of the leg gap genes. Rather than being restricted to a proximal leg region, as in arthropods and onychophorans, we detected coexpression of orthologues of homothorax and extradenticle broadly across the legs of the first three trunk segments in the tardigrade Hypsibius exemplaris. We could not identify a dachshund orthologue in tardigrade genomes, a gene that is expressed in an intermediate region of developing legs in arthropods and onychophorans, suggesting that this gene was lost in the tardigrade lineage. We detected Distal-less expression broadly across all developing leg buds in H. exemplaris embryos, unlike in arthropods and onychophorans, in which it exhibits a distally restricted expression domain. The broad expression patterns of the remaining leg gap genes in H. exemplaris legs may reflect the loss of dachshund and the accompanying loss of an intermediate region of the legs in the tardigrade lineage. We propose that the loss of intermediate regions of both the AP and PD body axes contributed to miniaturization of Tardigrada.

Desiccation does not drastically increase the accessibility of exogenous DNA to nuclear genomes: evidence from the frequency of endosymbiotic DNA transfer

Background

Although horizontal gene transfer (HGT) is a widely accepted force in the evolution of prokaryotic genomes, its role in the evolution of eukaryotic genomes remains hotly debated. Some bdelloid rotifers that are resistant to extreme desiccation and radiation undergo a very high level of HGT, whereas in another desiccation-resistant invertebrate, the tardigrade, the pattern does not exist. Overall, the DNA double-strand breaks (DSBs) induced by prolonged desiccation have been postulated to open a gateway to the nuclear genome for exogenous DNA integration and thus to facilitate the HGT process, thereby enhancing the rate of endosymbiotic DNA transfer (EDT).

Results

We first surveyed the abundance of nuclear mitochondrial DNAs (NUMTs) and nuclear plastid DNAs (NUPTs) in five eukaryotes that are highly resistant to desiccation: the bdelloid rotifers Adineta vaga and Adineta ricciae, the tardigrade Ramazzottius varieornatus, and the resurrection plants Dorcoceras hygrometricum and Selaginella tamariscina. Excessive NUMTs or NUPTs were not detected. Furthermore, we compared 24 groups of desiccation-tolerant organisms with their relatively less desiccation-tolerant relatives but did not find a significant difference in NUMT/NUPT contents.

Conclusions

Desiccation may induce DSBs, but it is unlikely to dramatically increase the frequency of exogenous sequence integration in most eukaryotes. The capture of exogenous DNA sequences is possible only when DSBs are repaired through a subtype of non-homologous end joining, named alternative end joining (alt-EJ). Due to the deleterious effects of the resulting insertion mutations, alt-EJ is less frequently initiated than other mechanisms.

Horizontal Gene Transfer: From Evolutionary Flexibility to Disease Progression

Flexibility in the exchange of genetic material takes place between different organisms of the same or different species. This phenomenon is known to play a key role in the genetic, physiological, and ecological performance of the host. Exchange of genetic materials can cause both beneficial and/or adverse biological consequences. Horizontal gene transfer (HGT) or lateral gene transfer (LGT) as a general mechanism leads to biodiversity and biological innovations in nature. HGT mediators are one of the genetic engineering tools used for selective introduction of desired changes in the genome for gene/cell therapy purposes. HGT, however, is crucial in development, emergence, and recurrence of various human-related diseases, such as cancer, genetic-, metabolic-, and neurodegenerative disorders and can negatively affect the therapeutic outcome by promoting resistant forms or disrupting the performance of genome editing toolkits. Because of the importance of HGT and its vital physio- and pathological roles, here the variety of HGT mechanisms are reviewed, ranging from extracellular vesicles (EVs) and nanotubes in prokaryotes to cell-free DNA and apoptotic bodies in eukaryotes. Next, we argue that HGT plays a role both in the development of useful features and in pathological states associated with emerging and recurrent forms of the disease. A better understanding of the different HGT mediators and their genome-altering effects/potentials may pave the way for the development of more effective therapeutic and diagnostic regimes.

Accurate and complete genomes from metagenomes

Genomes are an integral component of the biological information about an organism; thus, the more complete the genome, the more informative it is. Historically, bacterial and archaeal genomes were reconstructed from pure (monoclonal) cultures, and the first reported sequences were manually curated to completion. However, the bottleneck imposed by the requirement for isolates precluded genomic insights for the vast majority of microbial life. Shotgun sequencing of microbial communities, referred to initially as community genomics and subsequently as genome-resolved metagenomics, can circumvent this limitation by obtaining metagenome-assembled genomes (MAGs); but gaps, local assembly errors, chimeras, and contamination by fragments from other genomes limit the value of these genomes. Here, we discuss genome curation to improve and, in some cases, achieve complete (circularized, no gaps) MAGs (CMAGs). To date, few CMAGs have been generated, although notably some are from very complex systems such as soil and sediment. Through analysis of about 7000 published complete bacterial isolate genomes, we verify the value of cumulative GC skew in combination with other metrics to establish bacterial genome sequence accuracy. The analysis of cumulative GC skew identified potential misassemblies in some reference genomes of isolated bacteria and the repeat sequences that likely gave rise to them. We discuss methods that could be implemented in bioinformatic approaches for curation to ensure that metabolic and evolutionary analyses can be based on very high-quality genomes.

Ultrastructural analysis of the dehydrated tardigrade Hypsibius exemplaris unveils an anhydrobiotic-specific architecture

Tardigrades can cope with adverse environmental conditions by turning into anhydrobiotes with a characteristic tun shape. Tun formation is an essential morphological adaptation for tardigrade entry into the anhydrobiotic state. The tun cell structure and ultrastructure have rarely been explored in tardigrades in general and never in Hypsibius exemplaris. We used transmission electron microscopy to compare cellular organization and ultrastructures between hydrated and anhydrobiotic H. exemplaris. Despite a globally similar cell organelle structure and a number of cells not significantly different between hydrated and desiccated tardigrades, reductions in the sizes of both cells and mitochondria were detected in dehydrated animals. Moreover, in anhydrobiotes, secretory active cells with a dense endoplasmic reticulum network were observed. Interestingly, these anhydrobiote-specific cells are in a close relationship with a specific extracellular structure surrounding each cell. It is possible that this rampart-like extracellular structure resulted from the accumulation of anhydrobiotic-specific material to protect the cells. Interestingly, after five hours of rehydration, the number of secretory cells decreased, and the specific extracellular structure began to disappear. Twenty-four hours after the beginning of rehydration, the cellular structure and ultrastructure were comparable to those observed in hydrated tardigrades.

Transposable Elements Cross Kingdom Boundaries and Contribute to Inflammation and Ageing: Somatic Acquisition of Foreign Transposable Elements as a Catalyst of Genome Instability, Epigenetic Dysregulation, Inflammation, Senescence, and Ageing

The de-repression of transposable elements (TEs) in mammalian genomes is thought to contribute to genome instability, inflammation, and ageing, yet is viewed as a cell-autonomous event. In contrast to mammalian cells, prokaryotes constantly exchange genetic material through TEs, crossing both cell and species barriers, contributing to rapid microbial evolution and diversity in complex communities such as the mammalian gut. Here, it is proposed that TEs released from prokaryotes in the microbiome or from pathogenic infections regularly cross the kingdom barrier to the somatic cells of their eukaryotic hosts. It is proposed this horizontal transfer of TEs from microbe to host is a stochastic, ongoing catalyst of genome destabilization, resulting in structural and epigenetic variations, and activation of well-evolved host defense mechanisms contributing to inflammation, senescence, and biological ageing. It is proposed that innate immunity pathways defend against the horizontal acquisition of microbial TEs, and that activation of this pathway during horizontal transposon transfer promotes chronic inflammation during ageing. Finally, it is suggested that horizontal acquisition of prokaryotic TEs into mammalian genomes has been masked and subsequently under-reported due to flaws in current sequencing pipelines, and new strategies to uncover these events are proposed.

Integrative description of bisexual Paramacrobiotus experimentalis sp. nov. (Macrobiotidae) from republic of Madagascar (Africa) with microbiome analysis

In a moss samples collected on Madagascar two populations of Paramacrobiotus experimentalis sp. nov. were found. Paramacrobiotus experimentalis sp. nov. with the presence of a microplacoid and areolatus type of eggs is similar to Pam. danielae, Pam. garynahi, Pam. hapukuensis, Pam. peteri, Pam. rioplatensis and Pam. savai, but it differs from them by some morphological and morphometric characters of the eggs. The p-distance between two COI haplotypes of Pam. experimentalis sp. nov. was 0.17%. In turn, the ranges of uncorrected genetic p-distances of all Paramacrobiotus species available in GenBank was from 18.27% (for Pam. lachowskae) to 25.26% (for Pam. arduus) with an average distance of 20.67%. We also found that Pam. experimentalis sp. nov. is bisexual. This observation was congruent on three levels: (i) morphological - specimen size dimorphism; (ii) structural (primary sexual characteristics) - females have an unpaired ovary while males have an unpaired testis and (iii) molecular - heterozygous and homozygous strains of the ITS-2 marker. Although symbiotic associations of hosts with bacteria (including endosymbiotic bacteria) are common in nature and these interactions exert various effects on the evolution, biology and reproductive ecology of hosts, there is still very little information on the bacterial community associated with tardigrades. To fill this gap and characterise the bacterial community of Pam. experimentalis sp. nov. populations and microbiome of its microhabitat, high throughput sequencing of the V3-V4 hypervariable regions in the bacterial 16S rRNA gene fragment was performed. The obtained 16S rRNA gene sequences ranged from 92,665 to 131,163. In total, 135 operational taxonomic units (OTUs) were identified across the rarefied dataset. Overall, both Pam. experimentalis sp. nov. populations were dominated by OTUs ascribed to the phylum Proteobacteria (89-92%) and Firmicutes (6-7%). In the case of samples from tardigrades' laboratory habitat, the most abundant bacterial phylum was Proteobacteria (51-90%) and Bacteroides (9-48%). In all compared microbiome profiles, only 16 of 137 OTUs were shared. We found also significant differences in beta diversity between the partly species-specific microbiome of Pam. experimentalis sp. nov. and its culturing environment. Two OTUs belonging to a putative bacterial endosymbiont were identified - Rickettsiales and Polynucleobacter. We also demonstrated that each bacterial community was rich in genes involved in membrane transport, amino acid metabolism, and carbohydrate metabolism.

Finding relationships among biological entities

Confusion over the concepts of “relationships” and “similarities” lies at the heart of many battles over the direction and intent of research projects. Here is a short story that demonstrates the difference between the two concepts: You look up at the clouds, and you begin to see the shape of a lion. The cloud has a tail, like a lion’s tale, and a fluffy head, like a lion’s mane. With a little imagination the mouth of the lion seems to roar down from the sky. You have succeeded in finding similarities between the cloud and a lion. If you look at a cloud and you imagine a tea kettle producing a head of steam and you recognize that the physical forces that create a cloud and the physical forces that produced steam from a heated kettle are the same, then you have found a relationship. Most popular classification algorithms operate by grouping together data objects that have similar properties or values. In so doing, they may miss finding the true relationships among objects. Traditionally, relationships among data objects are discovered by an intellectual process. In this chapter, we will discuss the scientific gains that come when we classify biological entities by relationships, not by their similarities.

On the evolutionary significance of horizontal gene transfers in plants

Horizontal gene transfer (HGT) has long been seen as a crucial process in the evolution of prokaryotic species, but until recently it was thought to have little, if any, effect on the evolution of eukaryotic life forms. Detecting and describing HGT events in eukaryotes is difficult, making this phenomenon at times controversial. However, modern advances in genomics and bioinformatics have radically altered our view of HGT in eukaryotes, especially in plants. It now appears that HGT to and from plant lineages is more common than previously suspected. Importantly, the transfer of functional nuclear genes with adaptive significance has been reported in numerous taxa. Here we review several recent studies that have found evidence of the horizontal transfer of nuclear genes, and argue that HGT has undoubtedly had profound impacts on plant evolution as a whole.

Further insights in the Tardigrada microbiome: phylogenetic position and prevalence of infection of four new Alphaproteobacteria putative endosymbionts

Data from a previous study showed that microbiomes of six tardigrade species are species-specific and distinct from associated environmental microbes. We here performed a more in-depth analyses of those data, to identify and characterize new potential symbionts. The most abundant bacterial operational taxonomic units (OTUs) found in tardigrades are classified, and their prevalence in other environments is assessed using public databases. A subset of OTUs was selected for molecular phylogenetic analyses based on their affiliation with host-associated bacterial families in tardigrades. Almost 22.6% of the most abundant OTUs found do not match any sequence at 99% identity in the IMNGS database. These novel OTUs include four putative tardigrade endosymbionts from Alphaproteobacteria (Anaplasmataceae and Candidatus Tenuibacteraceae), which are characterized by 16S rRNA gene analysis and investigated for their infection rates in: Echiniscus trisetosus, Richtersisus coronifer and Macrobiotus macrocalix. These putative endosymbionts have an infection prevalence between 9.1% and 40.0%, and are, therefore, likely secondary symbionts, not essential for tardigrade survival and reproduction. Using fluorescence in situ hybridization (FISH), we detected bacteria on the cuticle and within the ovary of E. trisetosus, suggesting possible vertical transmission. This study highlights the great contribution in biodiversity discovery that neglected phyla can provide in microbiome and symbiosis studies.

Mechanisms and evolution of resistance to environmental extremes in animals

When animals are exposed to an extreme environmental stress, one of three possible outcomes takes place: the animal dies, the animal avoids the environmental stress and survives, or the animal tolerates the environmental stress and survives. This review is concerned with the third possibility, and will look at mechanisms that rare animals use to survive extreme environmental stresses including freezing, desiccation, intense heat, irradiation, and low-oxygen conditions (hypoxia). In addition, an increasing understanding of the molecular mechanisms involved in environmental stress tolerance allows us to speculate on how these tolerances arose. Uncovering the mechanisms of extreme environmental stress tolerance and how they evolve has broad implications for our understanding of the evolution of early life on this planet, colonization of new environments, and the search for novel forms of life both on Earth and elsewhere, as well as a number of agricultural and health-related applications.

Acoustomicrofluidic separation of tardigrades from raw cultures for sample preparation

Tardigrades are microscopic animals widely known for their ability to survive in extreme conditions. They are the focus of current research in the fields of taxonomy, biogeography, genomics, proteomics, development, space biology, evolution and ecology. Tardigrades, such as Hypsibius exemplaris, are being advocated as a next-generation model organism for genomic and developmental studies. The raw culture of H. exemplaris usually contains tardigrades themselves, their eggs, faeces and algal food. Experimentation with tardigrades often requires the demanding and laborious separation of tardigrades from raw samples to prepare pure and contamination-free tardigrade samples. In this paper, we propose a two-step acoustomicrofluidic separation method to isolate tardigrades from raw samples. In the first step, a passive microfluidic filter composed of an array of traps is used to remove large algal clusters in the raw sample. In the second step, a surface acoustic wave-based active microfluidic separation device is used to deflect tardigrades continuously from their original streamlines inside the microchannel and thus isolate them selectively from algae and eggs. The experimental results demonstrated the efficient separation of tardigrades, with a recovery rate of 96% and an impurity of 4% algae on average in a continuous, contactless, automated, rapid and biocompatible manner.

Comparative transcriptomics suggest unique molecular adaptations within tardigrade lineages

Background

Tardigrades are renowned for their ability to enter cryptobiosis (latent life) and endure extreme stress, including desiccation and freezing. Increased focus is on revealing molecular mechanisms underlying this tolerance. Here, we provide the first transcriptomes from the heterotardigrade Echiniscoides cf. sigismundi and the eutardigrade Richtersius cf. coronifer, and compare these with data from other tardigrades and six eukaryote models. Investigating 107 genes/gene families, our study provides a thorough analysis of tardigrade gene content with focus on stress tolerance.

Results

E. cf. sigismundi, a strong cryptobiont, apparently lacks expression of a number of stress related genes. Most conspicuous is the lack of transcripts from genes involved in classical Non-Homologous End Joining. Our analyses suggest that post-cryptobiotic survival in tardigrades could rely on high fidelity transcription-coupled DNA repair. Tardigrades seem to lack many peroxins, but they all have a comprehensive number of genes encoding proteins involved in antioxidant defense. The “tardigrade unique proteins” (CAHS, SAHS, MAHS, RvLEAM), seem to be missing in the heterotardigrade lineage, revealing that cryptobiosis in general cannot be attributed solely to these proteins. Our investigation further reveals a unique and highly expressed cold shock domain. We hypothesize that the cold shock protein acts as a RNA-chaperone involved in regulation of translation following freezing.

Conclusions

Our results show common gene family contractions and expansions within stress related gene pathways in tardigrades, but also indicate that evolutionary lineages have a high degree of divergence. Different taxa and lineages may exhibit unique physiological adaptations towards stress conditions involving possible unknown functional homologues and/or novel physiological and biochemical mechanisms. To further substantiate the current results genome assemblies coupled with transcriptome data and experimental investigations are needed from tardigrades belonging to different evolutionary lineages.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5912-x) contains supplementary material, which is available to authorized users.

Horizontal Gene Transfer as an Indispensable Driver for Evolution of Neocallimastigomycota into a Distinct Gut-Dwelling Fungal Lineage

Survival and growth of the anaerobic gut fungi (AGF; Neocallimastigomycota) in the herbivorous gut necessitate the possession of multiple abilities absent in other fungal lineages. We hypothesized that horizontal gene transfer (HGT) was instrumental in forging the evolution of AGF into a phylogenetically distinct gut-dwelling fungal lineage. The patterns of HGT were evaluated in the transcriptomes of 27 AGF strains, 22 of which were isolated and sequenced in this study, and 4 AGF genomes broadly covering the breadth of AGF diversity. We identified 277 distinct incidents of HGT in AGF transcriptomes, with subsequent gene duplication resulting in an HGT frequency of 2 to 3.5% in AGF genomes. The majority of HGT events were AGF specific (91.7%) and wide (70.8%), indicating their occurrence at early stages of AGF evolution. The acquired genes allowed AGF to expand their substrate utilization range, provided new venues for electron disposal, augmented their biosynthetic capabilities, and facilitated their adaptation to anaerobiosis. The majority of donors were anaerobic fermentative bacteria prevalent in the herbivorous gut. This study strongly indicates that HGT indispensably forged the evolution of AGF as a distinct fungal phylum and provides a unique example of the role of HGT in shaping the evolution of a high-rank taxonomic eukaryotic lineage.IMPORTANCE The anaerobic gut fungi (AGF) represent a distinct basal phylum lineage (Neocallimastigomycota) commonly encountered in the rumen and alimentary tracts of herbivores. Survival and growth of anaerobic gut fungi in these anaerobic, eutrophic, and prokaryote-dominated habitats necessitates the acquisition of several traits absent in other fungal lineages. We assess here the role of horizontal gene transfer as a relatively fast mechanism for trait acquisition by the Neocallimastigomycota postsequestration in the herbivorous gut. Analysis of 27 transcriptomes that represent the broad diversity of Neocallimastigomycota identified 277 distinct HGT events, with subsequent gene duplication resulting in an HGT frequency of 2 to 3.5% in AGF genomes. These HGT events have allowed AGF to survive in the herbivorous gut by expanding their substrate utilization range, augmenting their biosynthetic pathway, providing new routes for electron disposal by expanding fermentative capacities, and facilitating their adaptation to anaerobiosis. HGT in the AGF is also shown to be mainly a cross-kingdom affair, with the majority of donors belonging to the bacteria. This study represents a unique example of the role of HGT in shaping the evolution of a high-rank taxonomic eukaryotic lineage.

Protist Evolution: Stealing Genes to Gut It Out

The amount and evolutionary impact of horizontal gene transfer in eukaryotes remain contentious issues. A new phylogenomic study suggests that gene transfer from prokaryotes has contributed significantly to the adaptation and metabolic evolution of Blastocystis, the most widespread human gut eukaryotic parasite.

The genomes of polyextremophilic cyanidiales contain 1% horizontally transferred genes with diverse adaptive functions

The role and extent of horizontal gene transfer (HGT) in eukaryotes are hotly disputed topics that impact our understanding of the origin of metabolic processes and the role of organelles in cellular evolution. We addressed this issue by analyzing 10 novel Cyanidiales genomes and determined that 1% of their gene inventory is HGT-derived. Numerous HGT candidates share a close phylogenetic relationship with prokaryotes that live in similar habitats as the Cyanidiales and encode functions related to polyextremophily. HGT candidates differ from native genes in GC-content, number of splice sites, and gene expression. HGT candidates are more prone to loss, which may explain the absence of a eukaryotic pan-genome. Therefore, the lack of a pan-genome and cumulative effects fail to provide substantive arguments against our hypothesis of recurring HGT followed by differential loss in eukaryotes. The maintenance of 1% HGTs, even under selection for genome reduction, underlines the importance of non-endosymbiosis related foreign gene acquisition.

Human contamination in bacterial genomes has created thousands of spurious proteins

Contaminant sequences that appear in published genomes can cause numerous problems for downstream analyses, particularly for evolutionary studies and metagenomics projects. Our large-scale scan of complete and draft bacterial and archaeal genomes in the NCBI RefSeq database reveals that 2250 genomes are contaminated by human sequence. The contaminant sequences derive primarily from high-copy human repeat regions, which themselves are not adequately represented in the current human reference genome, GRCh38. The absence of the sequences from the human assembly offers a likely explanation for their presence in bacterial assemblies. In some cases, the contaminating contigs have been erroneously annotated as containing protein-coding sequences, which over time have propagated to create spurious protein “families” across multiple prokaryotic and eukaryotic genomes. As a result, 3437 spurious protein entries are currently present in the widely used nr and TrEMBL protein databases. We report here an extensive list of contaminant sequences in bacterial genome assemblies and the proteins associated with them. We found that nearly all contaminants occurred in small contigs in draft genomes, which suggests that filtering out small contigs from draft genome assemblies may mitigate the issue of contamination while still keeping nearly all of the genuine genomic sequences.

Inferring Ancient Relationships with Genomic Data: A Commentary on Current Practices

Contemporary phylogeneticists enjoy an embarrassment of riches, not only in the volumes of data now available, but also in the diversity of bioinformatic tools for handling these data. Here, I discuss a subset of these tools I consider well-suited to the task of inferring ancient relationships with coding sequence data in particular, encompassing data generation, orthology assignment, alignment and gene tree inference, supermatrix construction, and analysis under the best-fitting models applicable to large-scale datasets. Throughout, I compare and critique methods, considering both their theoretical principles and the details of their implementation, and offering practical tips on usage where appropriate. I also entertain different motivations for analyzing what are almost always originally DNA sequence data as codons, amino acids, and higher-order recodings. Although presented in a linear order, I see value in using the diversity of tools available to us to assess the sensitivity of clades of biological interest to different gene and taxon sets and analytical modes, which can be an indication of the presence of systematic error, of which a few forms remain poorly controlled by even the best available inference methods.

Protecting activity of desiccated enzymes

Protein-based biological drugs and many industrial enzymes are unstable, making them prohibitively expensive. Some can be stabilized by formulation with excipients, but most still require low temperature storage. In search of new, more robust excipients, we turned to the tardigrade, a microscopic animal that synthesizes cytosolic abundant heat soluble (CAHS) proteins to protect its cellular components during desiccation. We find that CAHS proteins protect the test enzymes lactate dehydrogenase and lipoprotein lipase against desiccation-, freezing-, and lyophilization-induced deactivation. Our data also show that a variety of globular and disordered protein controls, with no known link to desiccation tolerance, protect our test enzymes. Protection of lactate dehydrogenase correlates, albeit imperfectly, with the charge density of the protein additive, suggesting an approach to tune protection by modifying charge. Our results support the potential use of CAHS proteins as stabilizing excipients in formulations and suggest that other proteins may have similar potential.

Lateral transfers of large DNA fragments spread functional genes among grasses

A fundamental tenet of multicellular eukaryotic evolution is that vertical inheritance is paramount, with natural selection acting on genetic variants transferred from parents to offspring. This lineal process means that an organism's adaptive potential can be restricted by its evolutionary history, the amount of standing genetic variation, and its mutation rate. Lateral gene transfer (LGT) theoretically provides a mechanism to bypass many of these limitations, but the evolutionary importance and frequency of this process in multicellular eukaryotes, such as plants, remains debated. We address this issue by assembling a chromosome-level genome for the grass Alloteropsis semialata, a species surmised to exhibit two LGTs, and screen it for other grass-to-grass LGTs using genomic data from 146 other grass species. Through stringent phylogenomic analyses, we discovered 57 additional LGTs in the A. semialata nuclear genome, involving at least nine different donor species. The LGTs are clustered in 23 laterally acquired genomic fragments that are up to 170 kb long and have accumulated during the diversification of Alloteropsis. The majority of the 59 LGTs in A. semialata are expressed, and we show that they have added functions to the recipient genome. Functional LGTs were further detected in the genomes of five other grass species, demonstrating that this process is likely widespread in this globally important group of plants. LGT therefore appears to represent a potent evolutionary force capable of spreading functional genes among distantly related grass species.