myTAI: evolutionary transcriptomics with R

Heat stress reprograms herbivory-induced defense responses in potato plants

Climate change is predicted to increase the occurrence of extreme weather events such as heatwaves, which may thereby impact the outcome of plant-herbivore interactions. While elevated temperature is known to directly affect herbivore growth, it remains largely unclear if it indirectly influences herbivore performance by affecting the host plant they feed on. In this study, we investigated how transient exposure to high temperature influences plant herbivory-induced defenses at the transcript and metabolic level. To this end, we studied the interaction between potato (Solanum tuberosum) plants and the larvae of the potato tuber moth (Phthorimaea operculella) under different temperature regimes. We found that P. operculella larvae grew heavier on leaves co-stressed by high temperature and insect herbivory than on leaves pre-stressed by herbivory alone. We also observed that high temperature treatments altered phylotranscriptomic patterns upon herbivory, which changed from an evolutionary hourglass pattern, in which transcriptomic responses at early and late time points after elicitation are more variable than the ones in the middle, to a vase pattern. Specifically, transcripts of many herbivory-induced genes in the early and late defense stage were suppressed by HT treatment, whereas those in the intermediate stage peaked earlier. Additionally, we observed that high temperature impaired the induction of jasmonates and defense compounds upon herbivory. Moreover, using jasmonate-reduced (JA-reduced, irAOC) and -elevated (JA-Ile-elevated, irCYP94B3s) potato plants, we showed that high temperature suppresses JA signaling mediated plant-induced defense to herbivore attack. Thus, our study provides evidences on how temperature reprograms plant-induced defense to herbivores.

Conserved and specific gene expression patterns in the embryonic development of tardigrades

Tardigrades, commonly known as water bears, are enigmatic organisms characterized by their remarkable resilience to extreme environments despite their simple and compact body structure. To date, there is still much to understand about their evolutionary and developmental features contributing to their special body plan and abilities. This research provides preliminary insights on the conserved and specific gene expression patterns during embryonic development of water bears, focusing on the species Hypsibius exemplaris. The developmental dynamic expression analysis of the genes with various evolutionary age grades indicated that the mid-conserved stage of H. exemplaris corresponds to the period of ganglia and midgut development, with the late embryonic stage showing a transition from non-conserved to conserved state. Additionally, a comparison with Drosophila melanogaster highlighted the absence of certain pathway nodes in development-related pathways, such as Maml and Hairless, which are respectively the transcriptional co-activator and co-repressor of NOTCH regulated genes. We also employed Weighted Gene Co-expression Network Analysis (WGCNA) to investigate the expression patterns of tardigrade-specific genes during embryo development. Our findings indicated that the module containing the highest proportion of tardigrade-specific genes (TSGs) exhibits high expression levels before the mid-conserved stage, potentially playing a role in glutathione and lipid metabolism. These functions may be associated to the ecdysone synthesis and storage cell formation, which is unique to tardigrades.

An evolutionary timeline of the oxytocin signaling pathway

Oxytocin is a neuropeptide associated with both psychological and somatic processes like parturition and social bonding. Although oxytocin homologs have been identified in many species, the evolutionary timeline of the entire oxytocin signaling gene pathway has yet to be described. Using protein sequence similarity searches, microsynteny, and phylostratigraphy, we assigned the genes supporting the oxytocin pathway to different phylostrata based on when we found they likely arose in evolution. We show that the majority (64%) of genes in the pathway are 'modern'. Most of the modern genes evolved around the emergence of vertebrates or jawed vertebrates (540 - 530 million years ago, 'mya'), including OXTR, OXT and CD38. Of those, 45% were under positive selection at some point during vertebrate evolution. We also found that 18% of the genes in the oxytocin pathway are 'ancient', meaning their emergence dates back to cellular organisms and opisthokonta (3500-1100 mya). The remaining genes (18%) that evolved after ancient and before modern genes were classified as 'medium-aged'. Functional analyses revealed that, in humans, medium-aged oxytocin pathway genes are highly expressed in contractile organs, while modern genes in the oxytocin pathway are primarily expressed in the brain and muscle tissue.

Genome assemblies of 11 bamboo species highlight diversification induced by dynamic subgenome dominance

Polyploidy (genome duplication) is a pivotal force in evolution. However, the interactions between parental genomes in a polyploid nucleus, frequently involving subgenome dominance, are poorly understood. Here we showcase analyses of a bamboo system (Poaceae: Bambusoideae) comprising a series of lineages from diploid (herbaceous) to tetraploid and hexaploid (woody), with 11 chromosome-level de novo genome assemblies and 476 transcriptome samples. We find that woody bamboo subgenomes exhibit stunning karyotype stability, with parallel subgenome dominance in the two tetraploid clades and a gradual shift of dominance in the hexaploid clade. Allopolyploidization and subgenome dominance have shaped the evolution of tree-like lignified culms, rapid growth and synchronous flowering characteristic of woody bamboos as large grasses. Our work provides insights into genome dominance in a remarkable polyploid system, including its dependence on genomic context and its ability to switch which subgenomes are dominant over evolutionary time.

Performance of a Shotgun Prediction Model for Colorectal Cancer When Using 16S rRNA Sequencing Data

Colorectal cancer (CRC), the third most common cancer globally, has shown links to disturbed gut microbiota. While significant efforts have been made to establish a microbial signature indicative of CRC using shotgun metagenomic sequencing, the challenge lies in validating this signature with 16S ribosomal RNA (16S) gene sequencing. The primary obstacle is reconciling the differing outputs of these two methodologies, which often lead to divergent statistical models and conclusions. In this study, we introduce an algorithm designed to bridge this gap by mapping shotgun-derived taxa to their 16S counterparts. This mapping enables us to assess the predictive performance of a shotgun-based microbiome signature using 16S data. Our results demonstrate a reduction in performance when applying the 16S-mapped taxa in the shotgun prediction model, though it retains statistical significance. This suggests that while an exact match between shotgun and 16S data may not yet be feasible, our approach provides a viable method for comparative analysis and validation in the context of CRC-associated microbiome research.

oggmap: a Python package to extract gene ages per orthogroup and link them with single-cell RNA data

SUMMARY

For model species, single-cell RNA-based cell atlases are available. A good cell atlas includes all major stages in a species' ontogeny, and soon, they will be standard even for nonmodel species. Here, we propose a Python package called oggmap, which allows for the easy extraction of an orthomap (gene ages per orthogroup) for any given query species from OrthoFinder and other gene family data resources, like homologous groups from eggNOG or PLAZA. oggmap provides extracted gene ages for more than thousand eukaryotic species which can be further used to calculate gene age-weighted expression data from scRNA sequencing objects using the Python Scanpy toolkit. Not limited to one transcriptome evolutionary index, oggmap can visualize the individual gene category (e.g. age class, nucleotide diversity bin) and their corresponding expression profiles to investigate scRNA-based cell type assignments in an evolutionary context.

AVAILABILITY AND IMPLEMENTATION

oggmap source code is available at https://github.com/kullrich/oggmap, documentation is available at https://oggmap.readthedocs.io/en/latest/. oggmap can be installed via PyPi or directly used via a docker container.

Discovery of putative long non-coding RNAs expressed in the eyes of Astyanax mexicanus (Actinopterygii: Characidae)

Astyanax mexicanus is a well-known model species, that has two morphotypes, cavefish, from subterranean rivers and surface fish, from surface rivers. They are morphologically distinct due to many troglomorphic traits in the cavefish, such as the absence of eyes. Most studies on A. mexicanus are focused on eye development and protein-coding genes involved in the process. However, lncRNAs did not get the same attention and very little is known about them. This study aimed to fill this knowledge gap, identifying, describing, classifying, and annotating lncRNAs expressed in the embryo's eye tissue of cavefish and surface fish. To do so, we constructed a concise workflow to assemble and evaluate transcriptomes, annotate protein-coding genes, ncRNAs families, predict the coding potential, identify putative lncRNAs, map them and predict interactions. This approach resulted in the identification of 33,069 and 19,493 putative lncRNAs respectively mapped in cavefish and surface fish. Thousands of these lncRNAs were annotated and identified as conserved in human and several species of fish. Hundreds of them were validated in silico, through ESTs. We identified lncRNAs associated with genes related to eye development. This is the case of a few lncRNAs associated with sox2, which we suggest being isomorphs of the SOX2-OT, a lncRNA that can regulate the expression of sox2. This work is one of the first studies to focus on the description of lncRNAs in A. mexicanus, highlighting several lncRNA targets and opening an important precedent for future studies focusing on lncRNAs expressed in A. mexicanus.

From head to rootlet: comparative transcriptomic analysis of a rhizocephalan barnacle Peltogaster reticulata (Crustacea: Rhizocephala)

Background: Rhizocephalan barnacles stand out in the diverse world of metazoan parasites. The body of a rhizocephalan female is modified beyond revealing any recognizable morphological features, consisting of the interna, a system of rootlets, and the externa, a sac-like reproductive body. Moreover, rhizocephalans have an outstanding ability to control their hosts, literally turning them into "zombies". Despite all these amazing traits, there are no genomic or transcriptomic data about any Rhizocephala. Methods: We collected transcriptomes from four body parts of an adult female rhizocephalan Peltogaster reticulata: the externa, and the main, growing, and thoracic parts of the interna. We used all prepared data for the de novo assembly of the reference transcriptome. Next, a set of encoded proteins was determined, the expression levels of protein-coding genes in different parts of the parasite's body were calculated and lists of enriched bioprocesses were identified. We also in silico identified and analyzed sets of potential excretory / secretory proteins. Finally, we applied phylostratigraphy and evolutionary transcriptomics approaches to our data.  Results: The assembled reference transcriptome included transcripts of 12,620 protein-coding genes and was the first for any rhizocephalan. Based on the results obtained, the spatial heterogeneity of protein-coding gene expression in different regions of the adult female body of P. reticulata was established. The results of both transcriptomic analysis and histological studies indicated the presence of germ-like cells in the lumen of the interna. The potential molecular basis of the interaction between the nervous system of the host and the parasite's interna was also determined. Given the prolonged expression of development-associated genes, we suggest that rhizocephalans "got stuck in their metamorphosis", even at the reproductive stage. Conclusions: The results of the first comparative transcriptomic analysis for Rhizocephala not only clarified but also expanded the existing ideas about the biology of these extraordinary parasites.

Integrated lipidomic and transcriptomic analyses reveal the mechanism of large yellow croaker roe phospholipids on lipid metabolism in normal-diet mice

Large yellow croaker roe phospholipids (LYCRPLs) could regulate the accumulation of triglycerides and blood lipid levels. However, there exists little research on the mechanism of LYCRPLs on lipid metabolism in normal-diet mice. In this work, the mice on a normal diet were given low-dose, medium-dose, and high-dose LYCRPLs by intragastric administration for 6 weeks. At the same time, the physiological and biochemical indicators of the mice were determined, and the histomorphological observation of the liver and epididymal fat was carried out. In addition, we examined the gene expression and lipid metabolites in the liver of mice using transcriptomic and lipidomic and performed a correlation analysis. The results showed that LYCRPLs regulated the lipid metabolism of normal-diet mice by affecting the expression of the glycerolipid metabolism pathway, insulin resistance pathway, and cholesterol metabolism pathway. This study not only elucidated the main pathway by which LYCRPLs regulate lipid metabolism, but also laid a foundation for exploring LYCRPLs as functional food supplements.

Oxytocin receptor expression patterns in the human brain across development

Oxytocin plays a vital role in social behavior and homeostatic processes, with animal models indicating that oxytocin receptor (OXTR) expression patterns in the brain influence behavior and physiology. However, the developmental trajectory of OXTR gene expression is unclear. By analyzing gene expression data in human post-mortem brain samples, from the prenatal period to late adulthood, we demonstrate distinct patterns of OXTR gene expression in the developing brain, with increasing OXTR expression along the course of the prenatal period culminating in a peak during early childhood. This early life OXTR expression peak pattern appears slightly earlier in a comparative macaque sample, which is consistent with the relative immaturity of the human brain during early life compared to macaques. We also show that a network of genes with strong spatiotemporal couplings with OXTR is enriched in several psychiatric illness and body composition phenotypes. Taken together, these results demonstrate that oxytocin signaling plays an important role in a diverse set of psychological and somatic processes across the lifespan.

From head to rootlet: comparative transcriptomic analysis of a rhizocephalan barnacle Peltogaster reticulata (Crustacea: Rhizocephala)

Background: Rhizocephalan barnacles stand out in the diverse world of metazoan parasites. The body of a rhizocephalan female is modified beyond revealing any recognizable morphological features, consisting of the interna, a system of rootlets, and the externa, a sac-like reproductive body. Moreover, rhizocephalans have an outstanding ability to control their hosts, literally turning them into "zombies". Despite all these amazing traits, there are no genomic or transcriptomic data about any Rhizocephala. Methods: We collected transcriptomes from four body parts of an adult female rhizocephalan Peltogaster reticulata: the externa, and the main, growing, and thoracic parts of the interna. We used all prepared data for the de novo assembly of the reference transcriptome. Next, a set of encoded proteins was determined, the expression levels of protein-coding genes in different parts of the parasite's body were calculated and lists of enriched bioprocesses were identified. We also in silico identified and analyzed sets of potential excretory / secretory proteins. Finally, we applied phylostratigraphy and evolutionary transcriptomics approaches to our data.  Results: The assembled reference transcriptome included transcripts of 12,620 protein-coding genes and was the first for any rhizocephalan. Based on the results obtained, the spatial heterogeneity of protein-coding gene expression in different regions of the adult female body of P. reticulata was established. The results of both transcriptomic analysis and histological studies indicated the presence of germ-like cells in the lumen of the interna. The potential molecular basis of the interaction between the nervous system of the host and the parasite's interna was also determined. Given the prolonged expression of development-associated genes, we suggest that rhizocephalans "got stuck in their metamorphosis", even at the reproductive stage. Conclusions: The results of the first comparative transcriptomic analysis for Rhizocephala not only clarified but also expanded the existing ideas about the biology of these extraordinary parasites.

Large-scale k-mer-based analysis of the informational properties of genomes, comparative genomics and taxonomy

Information theoretic approaches are ubiquitous and effective in a wide variety of bioinformatics applications. In comparative genomics, alignment-free methods, based on short DNA words, or k-mers, are particularly powerful. We evaluated the utility of varying k-mer lengths for genome comparisons by analyzing their sequence space coverage of 5805 genomes in the KEGG GENOME database. In subsequent analyses on four k-mer lengths spanning the relevant range (11, 21, 31, 41), hierarchical clustering of 1634 genus-level representative genomes using pairwise 21- and 31-mer Jaccard similarities best recapitulated a phylogenetic/taxonomic tree of life with clear boundaries for superkingdom domains and high subtree similarity for named taxons at lower levels (family through phylum). By analyzing ~14.2M prokaryotic genome comparisons by their lowest-common-ancestor taxon levels, we detected many potential misclassification errors in a curated database, further demonstrating the need for wide-scale adoption of quantitative taxonomic classifications based on whole-genome similarity.

New Genes Interacted with Recent Whole Genome Duplicates in the Fast Stem Growth of Bamboos

As drivers of evolutionary innovations, new genes allow organisms to explore new niches. However, clear examples of this process remain scarce. Bamboos, the unique grass lineage diversifying into the forest, have evolved with a key innovation of fast growth of woody stem, reaching up to 1 m/day. Here, we identify 1,622 bamboo-specific orphan genes that appeared in recent 46 million years, and 19 of them evolved from noncoding ancestral sequences with entire de novo origination process reconstructed. The new genes evolved gradually in exon−intron structure, protein length, expression specificity, and evolutionary constraint. These new genes, whether or not from de novo origination, are dominantly expressed in the rapidly developing shoots, and make transcriptomes of shoots the youngest among various bamboo tissues, rather than reproductive tissue in other plants. Additionally, the particularity of bamboo shoots has also been shaped by recent whole-genome duplicates (WGDs), which evolved divergent expression patterns from ancestral states. New genes and WGDs have been evolutionarily recruited into coexpression networks to underline fast-growing trait of bamboo shoot. Our study highlights the importance of interactions between new genes and genome duplicates in generating morphological innovation.

Relative qPCR to quantify colonization of plant roots by arbuscular mycorrhizal fungi

Arbuscular mycorrhiza fungi (AMF) are beneficial soil fungi that can promote the growth of their host plants. Accurate quantification of AMF in plant roots is important because the level of colonization is often indicative of the activity of these fungi. Root colonization is traditionally measured with microscopy methods which visualize fungal structures inside roots. Microscopy methods are labor-intensive, and results depend on the observer. In this study, we present a relative qPCR method to quantify AMF in which we normalized the AMF qPCR signal relative to a plant gene. First, we validated the primer pair AMG1F and AM1 in silico, and we show that these primers cover most AMF species present in plant roots without amplifying host DNA. Next, we compared the relative qPCR method with traditional microscopy based on a greenhouse experiment with Petunia plants that ranged from very high to very low levels of AMF root colonization. Finally, by sequencing the qPCR amplicons with MiSeq, we experimentally confirmed that the primer pair excludes plant DNA while amplifying mostly AMF. Most importantly, our relative qPCR approach was capable of discriminating quantitative differences in AMF root colonization and it strongly correlated (Spearman Rho = 0.875) with quantifications by traditional microscopy. Finally, we provide a balanced discussion about the strengths and weaknesses of microscopy and qPCR methods. In conclusion, the tested approach of relative qPCR presents a reliable alternative method to quantify AMF root colonization that is less operator-dependent than traditional microscopy and offers scalability to high-throughput analyses.

Embryo-Like Features in Developing Bacillus subtilis Biofilms

Correspondence between evolution and development has been discussed for more than two centuries. Recent work reveals that phylogeny-ontogeny correlations are indeed present in developmental transcriptomes of eukaryotic clades with complex multicellularity. Nevertheless, it has been largely ignored that the pervasive presence of phylogeny-ontogeny correlations is a hallmark of development in eukaryotes. This perspective opens a possibility to look for similar parallelisms in biological settings where developmental logic and multicellular complexity are more obscure. For instance, it has been increasingly recognized that multicellular behavior underlies biofilm formation in bacteria. However, it remains unclear whether bacterial biofilm growth shares some basic principles with development in complex eukaryotes. Here we show that the ontogeny of growing Bacillus subtilis biofilms recapitulates phylogeny at the expression level. Using time-resolved transcriptome and proteome profiles, we found that biofilm ontogeny correlates with the evolutionary measures, in a way that evolutionary younger and more diverged genes were increasingly expressed toward later timepoints of biofilm growth. Molecular and morphological signatures also revealed that biofilm growth is highly regulated and organized into discrete ontogenetic stages, analogous to those of eukaryotic embryos. Together, this suggests that biofilm formation in Bacillus is a bona fide developmental process comparable to organismal development in animals, plants, and fungi. Given that most cells on Earth reside in the form of biofilms and that biofilms represent the oldest known fossils, we anticipate that the widely adopted vision of the first life as a single-cell and free-living organism needs rethinking.

Carbon metabolism and transcriptome in developmental paths differentiation of a homokaryotic Coprinopsis cinerea strain

The balance and interplay between sexual and asexual reproduction is one of the most intriguing mysteries in the study of fungi. The choice of developmental strategy reflects the ability of fungi to adapt to the changing environment. However, the evolution of developmental paths and the metabolic regulation during differentiation and morphogenesis are poorly understood. Here, an analysis was performed of carbohydrate metabolism and gene expression regulation during the early differentiation process from the vegetative mycelium, to the differentiated structures, fruiting body, oidia and sclerotia, of a homokaryotic fruiting Coprinopsis cinerea strain A43mutB43mut pab1-1 #326. Changes during morphogenesis and the evolution of developmental strategies were followed. Conversion between glucose and glycogen and between glucose and beta-glucan were the main carbon flows in the differentiation processes. Genes related to carbohydrate transport and metabolism were significantly differentially expressed among paths. Sclerotia displayed a set of specifically up-regulated genes that were enriched in the carbon metabolism and energy production and conversion processes. Evolutionary transcriptomic analysis of four developmental paths showed that all transcriptomes were under the purifying selection, and the more stressful the environment, the younger the transcriptome age. Oidiation has the lowest value of transcriptome age index (TAI) and transcriptome divergence index (TDI), while the fruiting process has the highest of both indexes. These findings provide new insights into the regulations of carbon metabolism and gene expressions during the early stages of fungal developmental paths differentiation, and improve our understanding of the evolutionary process of life history and reproductive strategy in fungi.

phylostratr: a framework for phylostratigraphy

MOTIVATION

The goal of phylostratigraphy is to infer the evolutionary origin of each gene in an organism. This is done by searching for homologs within increasingly broad clades. The deepest clade that contains a homolog of the protein(s) encoded by a gene is that gene's phylostratum.

RESULTS

We have created a general R-based framework, phylostratr, to estimate the phylostratum of every gene in a species. The program fully automates analysis: selecting species for balanced representation, retrieving sequences, building databases, inferring phylostrata and returning diagnostics. Key diagnostics include: detection of genes with inferred homologs in old clades, but not intermediate ones; proteome quality assessments; false-positive diagnostics, and checks for missing organellar genomes. phylostratr allows extensive customization and systematic comparisons of the influence of analysis parameters or genomes on phylostrata inference. A user may: modify the automatically generated clade tree or use their own tree; provide custom sequences in place of those automatically retrieved from UniProt; replace BLAST with an alternative algorithm; or tailor the method and sensitivity of the homology inference classifier. We show the utility of phylostratr through case studies in Arabidopsis thaliana and Saccharomyces cerevisiae.

AVAILABILITY AND IMPLEMENTATION

Source code available at https://github.com/arendsee/phylostratr.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Evolutionary transcriptomics of metazoan biphasic life cycle supports a single intercalation origin of metazoan larvae

The transient larva-bearing biphasic life cycle is the hallmark of many metazoan phyla, but how metazoan larvae originated remains a major enigma in animal evolution. There are two hypotheses for larval origin. The 'larva-first' hypothesis suggests that the first metazoans were similar to extant larvae, with later evolution of the adult-added biphasic life cycle; the 'adult-first' hypothesis suggests that the first metazoans were adult forms, with the biphasic life cycle arising later via larval intercalation. Here, we investigate the evolutionary origin of primary larvae by conducting ontogenetic transcriptome profiling for Mollusca-the largest marine phylum characterized by a trochophore larval stage and highly variable adult forms. We reveal that trochophore larvae exhibit rapid transcriptome evolution with extraordinary incorporation of novel genes (potentially contributing to adult shell evolution), and that cell signalling/communication genes (for example, caveolin and innexin) are probably crucial for larval evolution. Transcriptome age analysis of eight metazoan species reveals the wide presence of young larval transcriptomes in both trochozoans and other major metazoan lineages, therefore arguing against the prevailing larva-first hypothesis. Our findings support an adult-first evolutionary scenario with a single metazoan larval intercalation, and suggest that the first appearance of proto-larva probably occurred after the divergence of direct-developing Ctenophora from a metazoan ancestor.

Gene expression variation in the brains of harvester ant foragers is associated with collective behavior

Natural selection on collective behavior acts on variation among colonies in behavior that is associated with reproductive success. In the red harvester ant (Pogonomyrmex barbatus), variation among colonies in the collective regulation of foraging in response to humidity is associated with colony reproductive success. We used RNA-seq to examine gene expression in the brains of foragers in a natural setting. We find that colonies differ in the expression of neurophysiologically-relevant genes in forager brains, and a fraction of these gene expression differences are associated with two colony traits: sensitivity of foraging activity to humidity, and forager brain dopamine to serotonin ratio. Loci that were correlated with colony behavioral differences were enriched in neurotransmitter receptor signaling & metabolic functions, tended to be more central to coexpression networks, and are evolving under higher protein-coding sequence constraint. Natural selection may shape colony foraging behavior through variation in gene expression.

Evidence for phylogenetically and catabolically diverse active diazotrophs in deep-sea sediment

Diazotrophic microorganisms regulate marine productivity by alleviating nitrogen limitation. However, we know little about the identity and activity of diazotrophs in deep-sea sediments, a habitat covering nearly two-thirds of the planet. Here, we identify candidate diazotrophs from Pacific Ocean sediments collected at 2893 m water depth using ¹⁵N-DNA stable isotope probing and a novel pipeline for nifH sequence analysis. Together, these approaches detect an unexpectedly diverse assemblage of active diazotrophs, including members of the Acidobacteria, Firmicutes, Nitrospirae, Gammaproteobacteria, and Deltaproteobacteria. Deltaproteobacteria, predominately members of the Desulfobacterales and Desulfuromonadales, are the most abundant diazotrophs detected, and display the most microdiversity of associated nifH sequences. Some of the detected lineages, including those within the Acidobacteria, have not previously been shown to fix nitrogen. The diazotrophs appear catabolically diverse, with the potential for using oxygen, nitrogen, iron, sulfur, and carbon as terminal electron acceptors. Therefore, benthic diazotrophy may persist throughout a range of geochemical conditions and provide a stable source of fixed nitrogen over geologic timescales. Our results suggest that nitrogen-fixing communities in deep-sea sediments are phylogenetically and catabolically diverse, and open a new line of inquiry into the ecology and biogeochemical impacts of deep-sea microorganisms.

Reconstructing the Transcriptional Ontogeny of Maize and Sorghum Supports an Inverse Hourglass Model of Inflorescence Development

Assembling meaningful comparisons between species is a major limitation in studying the evolution of organismal form. To understand development in maize and sorghum, closely related species with architecturally distinct inflorescences, we collected RNA-seq profiles encompassing inflorescence body-plan specification in both species. We reconstructed molecular ontogenies from 40 B73 maize tassels and 47 BTx623 sorghum panicles and separated them into transcriptional stages. To discover new markers of inflorescence development, we used random forest machine learning to determine stage by RNA-seq. We used two descriptions of transcriptional conservation to identify hourglass-like stages during inflorescence development. Despite a relatively short 12 million years since their last common ancestor, we found maize and sorghum inflorescences are most different during their hourglass-like stages of development, following an inverse-hourglass model of development. We discuss whether agricultural selection may account for the rapid divergence signatures in these species and the observed separation of evolutionary pressure and developmental reprogramming.

Pervasive convergent evolution and extreme phenotypes define chaperone requirements of protein homeostasis

Maintaining protein homeostasis is an essential requirement for cell and organismal viability. An elaborate regulatory system within cells, the protein homeostasis network, safeguards that proteins are correctly folded and functional. At the heart of this regulatory system lies a class of specialized protein quality control enzymes called chaperones that are tasked with assisting proteins in their folding, avoiding aggregation and degradation. Failure and decline of protein homeostasis are directly associated with conditions of aging and aging-related neurodegeneration. However, it is not clear what tips the balance of protein homeostasis and leads to onset of aging and diseases. Here, using a comparative genomics approach we report general principles of maintaining protein homeostasis across the eukaryotic tree of life. Expanding a previous study of 16 eukaryotes to the quantitative analysis of 216 eukaryotic genomes, we find a strong correlation between the composition of eukaryotic chaperone networks and genome complexity that is distinct for different species kingdoms. Organisms with pronounced phenotypes clearly buck this trend. Northobranchius furzeri, the shortest-lived vertebrate and a widely used model for fragile protein homeostasis, is found to be chaperone limited while Heterocephalus glaber as the longest-lived rodent and thus an especially robust organism is characterized by above-average numbers of chaperones. Strikingly, the relative size of chaperone networks is found to generally correlate with longevity in Metazoa. Our results thus indicate that the balance in protein homeostasis may be a key variable in explaining organismal robustness.

Elucidating the endogenous synovial fluid proteome and peptidome of inflammatory arthritis using label-free mass spectrometry

Background

Inflammatory arthritis (IA) is an immunological disorder in which loss of immune tolerance to endogenous self-antigens perpetuates synovitis and eventual destruction of the underlying cartilage and bone. Pathological changes in the joint are expected to be represented by synovial fluid (SF) proteins and peptides. In the present study, a mass spectrometry-based approach was utilized for the identification of key protein and peptide mediators of IA.

Methods

Age-matched SF samples from 10 rheumatoid arthritis patients, 10 psoriatic arthritis patients and 10 cadaveric controls were subjected to an integrated proteomic and peptidomic protocol using liquid chromatography tandem mass spectrometry. Significant differentially abundant proteins and peptides were identified between cohorts according to the results of a Mann-Whitney U test coupled to the Benjamini-Hochberg correction for multiple hypothesis testing. Fold change ratios were computed for each protein and peptide according to their log-transformed extracted ion current. Pathway analysis and antimicrobial peptide (AMP) prediction were conducted to clarify the pathophysiological relevance of identified proteins and peptides to IA.

Results

We determined that 144 proteins showed significant differential abundance between the IA and control SF proteomes, of which 11 protein candidates were selected for future follow-up studies. Similar analyses applied to our peptidomic data identified 15 peptide sequences, originating from 4 protein precursors, to have significant differential abundance in IA compared to the control SF peptidome. Pathway enrichment analysis of the IA SF peptidome along with AMP prediction suggests a possible mechanistic role of microbes in eliciting an immune response which drives the development of IA.

Conclusions

The discovery-phase data generated herein has provided a basis for the identification of candidates with the greatest potential to serve as novel serum biomarkers specific to inflammatory arthritides. Moreover, these findings facilitate the understanding of possible disease mechanisms specific to each subtype.

Gene Expression Does Not Support the Developmental Hourglass Model in Three Animals with Spiralian Development

It has been proposed that animals have a pattern of developmental evolution resembling an hourglass because the most conserved development stage—often called the phylotypic stage—is always in midembryonic development. Although the topic has been debated for decades, recent studies using molecular data such as RNA-seq gene expression data sets have largely supported the existence of periods of relative evolutionary conservation in middevelopment, consistent with the phylotypic stage and the hourglass concepts. However, so far this approach has only been applied to a limited number of taxa across the tree of life. Here, using established phylotranscriptomic approaches, we found a surprising reverse hourglass pattern in two molluscs and a polychaete annelid, representatives of the Spiralia, an understudied group that contains a large fraction of metazoan body plan diversity. These results suggest that spiralians have a divergent midembryonic stage, with more conserved early and late development, which is the inverse of the pattern seen in almost all other organisms where these phylotranscriptomic approaches have been reported. We discuss our findings in light of proposed reasons for the phylotypic stage and hourglass model in other systems.

Developmental Constraints on Genome Evolution in Four Bilaterian Model Species

Developmental constraints on genome evolution have been suggested to follow either an early conservation model or an "hourglass" model. Both models agree that late development strongly diverges between species, but debate on which developmental period is the most conserved. Here, based on a modified "Transcriptome Age Index" approach, that is, weighting trait measures by expression level, we analyzed the constraints acting on three evolutionary traits of protein coding genes (strength of purifying selection on protein sequences, phyletic age, and duplicability) in four species: Nematode worm Caenorhabditis elegans, fly Drosophila melanogaster, zebrafish Danio rerio, and mouse Mus musculus. In general, we found that both models can be supported by different genomic properties. Sequence evolution follows an hourglass model, but the evolution of phyletic age and of duplicability follow an early conservation model. Further analyses indicate that stronger purifying selection on sequences in the middle development are driven by temporal pleiotropy of these genes. In addition, we report evidence that expression in late development is enriched with retrogenes, which usually lack efficient regulatory elements. This implies that expression in late development could facilitate transcription of new genes, and provide opportunities for acquisition of function. Finally, in C. elegans, we suggest that dosage imbalance could be one of the main factors that cause depleted expression of high duplicability genes in early development.

A comparative transcriptional landscape of maize and sorghum obtained by single-molecule sequencing

Maize and sorghum are both important crops with similar overall plant architectures, but they have key differences, especially in regard to their inflorescences. To better understand these two organisms at the molecular level, we compared expression profiles of both protein-coding and noncoding transcripts in 11 matched tissues using single-molecule, long-read, deep RNA sequencing. This comparative analysis revealed large numbers of novel isoforms in both species. Evolutionarily young genes were likely to be generated in reproductive tissues and usually had fewer isoforms than old genes. We also observed similarities and differences in alternative splicing patterns and activities, both among tissues and between species. The maize subgenomes exhibited no bias in isoform generation; however, genes in the B genome were more highly expressed in pollen tissue, whereas genes in the A genome were more highly expressed in endosperm. We also identified a number of splicing events conserved between maize and sorghum. In addition, we generated comprehensive and high-resolution maps of poly(A) sites, revealing similarities and differences in mRNA cleavage between the two species. Overall, our results reveal considerable splicing and expression diversity between sorghum and maize, well beyond what was reported in previous studies, likely reflecting the differences in architecture between these two species.