Evolutionary rate covariation is a reliable predictor of co-functional interactions but not necessarily physical interactions

Co-functional proteins tend to have rates of evolution that covary over time. This correlation between evolutionary rates can be measured over the branches of a phylogenetic tree through methods such as evolutionary rate covariation (ERC), and then used to construct gene networks by the identification of proteins with functional interactions. The cause of this correlation has been hypothesized to result from both compensatory coevolution at physical interfaces and nonphysical forces such as shared changes in selective pressure. This study explores whether coevolution due to compensatory mutations has a measurable effect on the ERC signal. We examined the difference in ERC signal between physically interacting protein domains within complexes compared to domains of the same proteins that do not physically interact. We found no generalizable relationship between physical interaction and high ERC, although a few complexes ranked physical interactions higher than nonphysical interactions. Therefore, we conclude that coevolution due to physical interaction is weak, but present in the signal captured by ERC, and we hypothesize that the stronger signal instead comes from selective pressures on the protein as a whole and maintenance of the general function.

The First Genome-Wide Mildew Locus O Genes Characterization in the Lamiaceae Plant Family

Powdery mildew (PM) is a widespread plant disease that causes significant economic losses in thousands crops of temperate climates, including Lamiaceae species. Multiple scientific studies describe a peculiar form of PM-resistance associated at the inactivation of specific members of the Mildew Locus O (MLO) gene family, referred to as mlo-resistance. The characterization of Lamiaceae MLO genes, at the genomic level, would be a first step toward their potential use in breeding programs. We carried out a genome-wide characterization of the MLO gene family in 11 Lamiaceae species, providing a manual curated catalog of 324 MLO proteins. Evolutionary history and phylogenetic relationships were studied through maximum likelihood analysis and motif patter reconstruction. Our approach highlighted seven different clades diversified starting from an ancestral MLO domain pattern organized in 18 highly conserved motifs. In addition, 74 Lamiaceae putative PM susceptibility genes, clustering in clade V, were identified. Finally, we performed a codon-based evolutionary analysis, revealing a general high level of purifying selection in the eleven Lamiaceae MLO gene families, and the occurrence of few regions under diversifying selection in candidate susceptibility factors. The results of this work may help to address further biological questions concerning MLOs involved in PM susceptibility. In follow-up studies, it could be investigated whether the silencing or loss-of-function mutations in one or more of these candidate genes may lead to PM resistance.

Functional inference of long non-coding RNAs through exploration of highly conserved regions

Background: Long non-coding RNAs (lncRNAs), which are generally less functionally characterized or less annotated, evolve more rapidly than mRNAs and substantially possess fewer sequence conservation patterns than protein-coding genes across divergent species. People assume that the functional inference could be conducted on the evolutionarily conserved long non-coding RNAs as they are most likely to be functional. In the past decades, substantial progress has been made in discussions on the evolutionary conservation of non-coding genomic regions from multiple perspectives. However, understanding their conservation and the functions associated with sequence conservation in relation to further corresponding phenotypic variability or disorders still remains incomplete. Results: Accordingly, we determined a highly conserved region (HCR) to verify the sequence conservation among long non-coding RNAs and systematically profiled homologous long non-coding RNA clusters in humans and mice based on the detection of highly conserved regions. Moreover, according to homolog clustering, we explored the potential function inference via highly conserved regions on representative long non-coding RNAs. On lncRNA XACT, we investigated the potential functional competence between XACT and lncRNA XIST by recruiting miRNA-29a, regulating the downstream target genes. In addition, on lncRNA LINC00461, we examined the interaction relationship between LINC00461 and SND1. This interaction or association may be perturbed during the progression of glioma. In addition, we have constructed a website with user-friendly web interfaces for searching, analyzing, and downloading to present the homologous clusters of humans and mice. Conclusion: Collectively, homolog clustering via the highly conserved region definition and detection on long non-coding RNAs, as well as the functional explorations on representative sequences in our research, would provide new evidence for the potential function of long non-coding RNAs. Our results on the remarkable roles of long non-coding RNAs would presumably provide a new theoretical basis and candidate diagnostic indicators for tumors.

Studying Plant-Insect Interactions through the Analyses of the Diversity, Composition, and Functional Inference of Their Bacteriomes

As with many other trophic interactions, the interchange of microorganisms between plants and their herbivorous insects is unavoidable. To test the hypothesis that the composition and diversity of the insect bacteriome are driven by the bacteriome of the plant, the bacteriomes of both the plant Datura inoxia and its specialist insect Lema daturaphila were characterised using 16S sRNA gene amplicon sequencing. Specifically, the bacteriomes associated with seeds, leaves, eggs, guts, and frass were described and compared. Then, the functions of the most abundant bacterial lineages found in the samples were inferred. Finally, the patterns of co-abundance among both bacteriomes were determined following a multilayer network approach. In accordance with our hypothesis, most genera were shared between plants and insects, but their abundances differed significantly within the samples collected. In the insect tissues, the most abundant genera were Pseudomonas (24.64%) in the eggs, Serratia (88.46%) in the gut, and Pseudomonas (36.27%) in the frass. In contrast, the most abundant ones in the plant were Serratia (40%) in seeds, Serratia (67%) in foliar endophytes, and Hymenobacter (12.85%) in foliar epiphytes. Indeed, PERMANOVA analysis showed that the composition of the bacteriomes was clustered by sample type (F = 9.36, p < 0.001). Functional inferences relevant to the interaction showed that in the plant samples, the category of Biosynthesis of secondary metabolites was significantly abundant (1.4%). In turn, the category of Xenobiotics degradation and metabolism was significantly present (2.5%) in the insect samples. Finally, the phyla Proteobacteria and Actinobacteriota showed a pattern of co-abundance in the insect but not in the plant, suggesting that the co-abundance and not the presence−absence patterns might be more important when studying ecological interactions.

Inferring microbiota functions from taxonomic genes: a review

Deciphering microbiota functions is crucial to predict ecosystem sustainability in response to global change. High-throughput sequencing at the individual or community level has revolutionized our understanding of microbial ecology, leading to the big data era and improving our ability to link microbial diversity with microbial functions. Recent advances in bioinformatics have been key for developing functional prediction tools based on DNA metabarcoding data and using taxonomic gene information. This cheaper approach in every aspect serves as an alternative to shotgun sequencing. Although these tools are increasingly used by ecologists, an objective evaluation of their modularity, portability, and robustness is lacking. Here, we reviewed 100 scientific papers on functional inference and ecological trait assignment to rank the advantages, specificities, and drawbacks of these tools, using a scientific benchmarking. To date, inference tools have been mainly devoted to bacterial functions, and ecological trait assignment tools, to fungal functions. A major limitation is the lack of reference genomes-compared with the human microbiota-especially for complex ecosystems such as soils. Finally, we explore applied research prospects. These tools are promising and already provide relevant information on ecosystem functioning, but standardized indicators and corresponding repositories are still lacking that would enable them to be used for operational diagnosis.

Updates on metaQuantome Software for Quantitative Metaproteomics

metaQuantome is a software suite that enables the quantitative analysis, statistical evaluation. and visualization of mass-spectrometry-based metaproteomics data. In the latest update of this software, we have provided several extensions, including a step-by-step training guide, the ability to perform statistical analysis on samples from multiple conditions, and a comparative analysis of metatranscriptomics data. The training module, accessed via the Galaxy Training Network, will help users to use the suite effectively both for functional as well as for taxonomic analysis. We extend the ability of metaQuantome to now perform multi-data-point quantitative and statistical analyses so that studies with measurements across multiple conditions, such as time-course studies, can be analyzed. With an eye on the multiomics analysis of microbial communities, we have also initiated the use of metaQuantome statistical and visualization tools on outputs from metatranscriptomics data, which complements the metagenomic and metaproteomic analyses already available. For this, we have developed a tool named MT2MQ ("metatranscriptomics to metaQuantome"), which takes in outputs from the ASaiM metatranscriptomics workflow and transforms them so that the data can be used as an input for comparative statistical analysis and visualization via metaQuantome. We believe that these improvements to metaQuantome will facilitate the use of the software for quantitative metaproteomics and metatranscriptomics and will enable multipoint data analysis. These improvements will take us a step toward integrative multiomic microbiome analysis so as to understand dynamic taxonomic and functional responses of these complex systems in a variety of biological contexts. The updated metaQuantome and MT2MQ are open-source software and are available via the Galaxy Toolshed and GitHub.

The constitution and functional prediction of the microbiota in necrotizing enterocolitis with a gestational age of over 28 weeks

BACKGROUND

To explore the features and function of gut microbiota in necrotizing enterocolitis patients over 28 gestational age weeks through a case-control study.

METHODS

Fecal samples from patients with NEC over 28 gestational week age and matched control cases were collected. DNA of the fecal samples was extracted for 16 s rRNA sequencing to estimate the composition of the microbiota. Functional inference analyses were conducted through PICRUSt based on the sequencing raw data.

RESULTS

There was no significant difference in the total diversity of microbiota between the fecal samples from the patients with NEC and the controls (P = .40). Propionibacterium was more abundant in the NEC cases than in the controls. Conversely, Lactobacillus, Phascolarctobacterium, and Streptococcus_salivarius were found to be more plentiful in the controls through LEfSe analysis. Functional inference analysis revealed that the xenobiotic biodegradation and metabolic activity was lower in the NEC cases than in the controls (P < .05).

CONCLUSION

The NEC cohort with a gestational age of over 28 weeks has a different pattern of microbiota compared with the controls. Functional inference analysis indicated that the potential function of the microbiota may also differ between these groups.

The chordate proteome history database

The chordate proteome history database (http://ioda.univ-provence.fr) comprises some 20,000 evolutionary analyses of proteins from chordate species. Our main objective was to characterize and study the evolutionary histories of the chordate proteome, and in particular to detect genomic events and automatic functional searches. Firstly, phylogenetic analyses based on high quality multiple sequence alignments and a robust phylogenetic pipeline were performed for the whole protein and for each individual domain. Novel approaches were developed to identify orthologs/paralogs, and predict gene duplication/gain/loss events and the occurrence of new protein architectures (domain gains, losses and shuffling). These important genetic events were localized on the phylogenetic trees and on the genomic sequence. Secondly, the phylogenetic trees were enhanced by the creation of phylogroups, whereby groups of orthologous sequences created using OrthoMCL were corrected based on the phylogenetic trees; gene family size and gene gain/loss in a given lineage could be deduced from the phylogroups. For each ortholog group obtained from the phylogenetic or the phylogroup analysis, functional information and expression data can be retrieved. Database searches can be performed easily using biological objects: protein identifier, keyword or domain, but can also be based on events, eg, domain exchange events can be retrieved. To our knowledge, this is the first database that links group clustering, phylogeny and automatic functional searches along with the detection of important events occurring during genome evolution, such as the appearance of a new domain architecture.

Nonparametric Residue Analysis of Dynamic PET Data With Application to Cerebral FDG Studies in Normals

Kinetic analysis is used to extract metabolic information from dynamic positron emission tomography (PET) uptake data. The theory of indicator dilutions, developed in the seminal work of Meier and Zierler (1954), provides a probabilistic framework for representation of PET tracer uptake data in terms of a convolution between an arterial input function and a tissue residue. The residue is a scaled survival function associated with tracer residence in the tissue. Nonparametric inference for the residue, a deconvolution problem, provides a novel approach to kinetic analysis-critically one that is not reliant on specific compartmental modeling assumptions. A practical computational technique based on regularized cubic B-spline approximation of the residence time distribution is proposed. Nonparametric residue analysis allows formal statistical evaluation of specific parametric models to be considered. This analysis needs to properly account for the increased flexibility of the nonparametric estimator. The methodology is illustrated using data from a series of cerebral studies with PET and fluorodeoxyglucose (FDG) in normal subjects. Comparisons are made between key functionals of the residue, tracer flux, flow, etc., resulting from a parametric (the standard two-compartment of Phelps et al. 1979) and a nonparametric analysis. Strong statistical evidence against the compartment model is found. Primarily these differences relate to the representation of the early temporal structure of the tracer residence-largely a function of the vascular supply network. There are convincing physiological arguments against the representations implied by the compartmental approach but this is the first time that a rigorous statistical confirmation using PET data has been reported. The compartmental analysis produces suspect values for flow but, notably, the impact on the metabolic flux, though statistically significant, is limited to deviations on the order of 3%-4%. The general advantage of the nonparametric residue analysis is the ability to provide a valid kinetic quantitation in the context of studies where there may be heterogeneity or other uncertainty about the accuracy of a compartmental model approximation of the tissue residue.