Metabolic phenotyping of marine heterotrophs on refactored media reveals diverse metabolic adaptations and lifestyle strategies

Microbial communities, through their metabolism, drive carbon cycling in marine environments. These complex communities are composed of many different microorganisms including heterotrophic bacteria, each with its own nutritional needs and metabolic capabilities. Yet, models of ecosystem processes typically treat heterotrophic bacteria as a "black box", which does not resolve metabolic heterogeneity nor address ecologically important processes such as the successive modification of different types of organic matter. Here we directly address the heterogeneity of metabolism by characterizing the carbon source utilization preferences of 63 heterotrophic bacteria representative of several major marine clades. By systematically growing these bacteria on 10 media containing specific subsets of carbon sources found in marine biomass, we obtained a phenotypic fingerprint that we used to explore the relationship between metabolic preferences and phylogenetic or genomic features. At the class level, these bacteria display broadly conserved patterns of preference for different carbon sources. Despite these broad taxonomic trends, growth profiles correlate poorly with phylogenetic distance or genome-wide gene content. However, metabolic preferences are strongly predicted by a handful of key enzymes that preferentially belong to a few enriched metabolic pathways, such as those involved in glyoxylate metabolism and biofilm formation. We find that enriched pathways point to enzymes directly involved in the metabolism of the corresponding carbon source and suggest potential associations between metabolic preferences and other ecologically-relevant traits. The availability of systematic phenotypes across multiple synthetic media constitutes a valuable resource for future quantitative modeling efforts and systematic studies of inter-species interactions.

Genomic diversification of long polynucleotide fragments is a signature of emerging SARS-CoV-2 variants of concern

Highly transmissible or immuno-evasive SARS-CoV-2 variants have intermittently emerged and outcompeted previously circulating strains, resulting in repeated COVID-19 surges, reinfections, and breakthrough infections in vaccinated individuals. With over 5 million SARS-CoV-2 genomes sequenced globally over the last 2 years, there is unprecedented data to decipher how competitive viral evolution results in the emergence of fitter SARS-CoV-2 variants. Much attention has been directed to studying how specific mutations in the Spike protein impact its binding to the ACE2 receptor or viral neutralization by antibodies, but there is limited knowledge of genomic signatures shared primarily by dominant variants. Here we introduce a methodology to quantify the genome-wide distinctiveness of polynucleotide fragments of various lengths (3- to 240-mers) that constitute SARS-CoV-2 lineage genomes. Compared to standard phylogenetic distance metrics and overall mutational load, the quantification of distinctive 9-mer polynucleotides provides a higher resolution of separation between variants of concern (Reference = 89, IQR: 65-108; Alpha = 166, IQR: 150-182; Beta 130, IQR: 113-147; Gamma = 165, IQR: 152-180; Delta = 234, IQR: 216-253; and Omicron = 294, IQR: 287-315). The similar scoring of the Alpha and Gamma variants by our methodology is consistent with these strains emerging at approximately the same time and circulating in distinct geographical regions as dominant strains. Furthermore, evaluation of genomic distinctiveness for 1,363 lineages annotated in GISAID highlights that polynucleotide diversity has increased over time (R2 = 0.37) and that VOCs show high distinctiveness compared to non-VOC contemporary lineages. To facilitate similar real-time assessments on the competitive fitness potential of future variants, we are launching a freely accessible resource for infusing pandemic preparedness with genomic inference ("GENI" — https://academia.nferx.com/GENI). This study demonstrates the value of characterizing new SARS-CoV-2 variants by their genome-wide polynucleotide distinctiveness and emphasizes the need to go beyond a narrow set of mutations at known functionally salient sites.

Reconstructing mitochondrial genomes from ancient DNA through iterative mapping: an evaluation of software, parameters, and bait reference

(1) Within evolutionary biology, mitochondrial genomes (mitogenomes) provide useful insights at both population and species level. Several approaches are available to assemble mitogenomes. However, most are not suitable for divergent, extinct species, due to the requirement of a reference mitogenome from a conspecific or close relative, and relatively high-quality DNA. (2) Iterative mapping can overcome the lack of a close reference sequence, and has been applied to an array of extinct species. Despite its widespread use, the accuracy of the reconstructed assemblies are yet to be comprehensively assessed. Here, we investigated the influence of mapping software (BWA or MITObim), parameters, and bait reference phylogenetic distance on the accuracy of the reconstructed assembly using two simulated datasets: (i) spotted hyena and various mammalian bait references, and (ii) southern cassowary and various avian bait references. Specifically, we assessed the accuracy of results through pairwise distance (PWD) to the reference conspecific mitogenome, number of incorrectly inserted base pairs (bp), and total length of the reconstructed assembly. (3) We found large discrepancies in the accuracy of reconstructed assemblies using different mapping software, parameters, and bait references. PWD to the reference conspecific mitogenome, which reflected the level of incorrect base calls, was consistently higher with BWA than MITObim. The same was observed for the number of incorrectly inserted bp. In contrast, the total sequence length was lower. Overall, the most accurate results were obtained with MITObim using mismatch values of 3 or 5, and the phylogenetically closest bait reference sequence. Accuracy could be further improved by combining results from multiple bait references. (4) We present the first comprehensive investigation of how mapping software, parameters, and bait reference influence mitogenome reconstruction from ancient DNA through iterative mapping. Our study provides information on how mitogenomes are best reconstructed from divergent, short-read data. By obtaining the most accurate reconstruction possible, one can be more confident as to the reliability of downstream analyses, and the evolutionary inferences made from them.

Intact rDNA arrays of Potentilla-origin detected in Erythronium nucleus suggest recent eudicot-to-monocot horizontal transfer

The occurrence of horizontal gene transfer (HGT) in Eukarya is increasingly gaining recognition. Nuclear-to-nuclear jump of DNA between plant species at high phylogenetic distance and devoid of intimate association (e.g., parasitism) is still scarcely reported. Within eukaryotes, components of ribosomal DNA (rDNA) multigene family have been found to be horizontally transferred in protists, fungi and grasses. However, in neither case HGT occurred between phylogenetic families, nor the transferred rDNA remained tandemly arrayed and transcriptionally active in the recipient organism. This study aimed to characterize an alien eudicot-type of 45S nuclear rDNA, assumingly transferred horizontally to the genome of monocot European Erythronium (Liliaceae). Genome skimming coupled by PacBio HiFi sequencing of a BAC clone were applied to determine DNA sequence of the alien rDNA. A clear phylogenetic signal traced the origin of the alien rDNA of Erythronium back to the Argentea clade of Potentilla (Rosaceae) and deemed the transfer to have occurred in the common ancestor of E. dens-canis and E. caucasicum. Though being discontinuous, transferred rDNA preserved its general tandemly arrayed feature in the host organism. Southern blotting, molecular cytogenetics, and sequencing of a BAC clone derived from flow-sorted nuclei indicated integration of the alien rDNA into the recipient's nuclear genome. Unprecedently, dicot-type alien rDNA was found to be transcribed in the monocot Erythronium albeit much less efficiently than the native counterpart. This study adds a new example to the growing list of naturally transgenic plants while holding the scientific community continually in suspense about the mode of DNA transfer.

Darwin’s naturalization conundrum disentangled: the role of phylogenetic relatedness depends on the invasion stages

Darwin’s naturalization hypothesis predicts successful invaders to be distantly related to native species, whereas his pre-adaptation hypothesis predicts the opposite. It has been suggested that depending on the invasion stage (i.e. introduction, naturalization, and invasiveness), both hypotheses, now known as Darwin’s naturalization conundrum, could hold true. We tested this by analysing whether the likelihood of introduction for cultivation as well as subsequent stages of naturalization and invasion of species alien to Southern Africa are correlated with their phylogenetic distance to the native flora of this region. While species were more likely to be introduced for cultivation if they are distantly related to the native flora, the probability of subsequent naturalization was higher for species closely related to the native flora. Furthermore, the probability of becoming invasive was higher for naturalized species distantly related to the native flora. These results were consistent across three different metrics of phylogenetic distance. Our study reveals that the relationship between phylogenetic distance to the native flora and success of an alien species depends on the invasion stage.

Detailed Evolutionary Analyses of the F Gene in the Respiratory Syncytial Virus Subgroup A

We performed evolution, phylodynamics, and reinfection-related antigenicity analyses of respiratory syncytial virus subgroup A (RSV-A) fusion (F) gene in globally collected strains (1465 strains) using authentic bioinformatics methods. The time-scaled evolutionary tree using the Bayesian Markov chain Monte Carlo method estimated that a common ancestor of the RSV-A, RSV-B, and bovine-RSV diverged at around 450 years ago, and RSV-A and RSV-B diverged around 250 years ago. Finally, the RSV-A F gene formed eight genotypes (GA1‑GA7 and NA1) over the last 80 years. Phylodynamics of RSV-A F gene, including all genotype strains, increased twice in the 1990s and 2010s, while patterns of each RSV-A genotype were different. Phylogenetic distance analysis suggested that the genetic distances of the strains were relatively short (less than 0.05). No positive selection sites were estimated, while many negative selection sites were found. Moreover, the F protein 3D structure mapping and conformational epitope analysis implied that the conformational epitopes did not correspond to the neutralizing antibody binding sites of the F protein. These results suggested that the RSV-A F gene is relatively conserved, and mismatches between conformational epitopes and neutralizing antibody binding sites of the F protein are responsible for the virus reinfection.

Correlation between acoustic divergence and phylogenetic distance in soniferous European gobiids (Gobiidae; Gobius lineage)

In fish, species identity can be encoded by sounds, which have been thoroughly investigated in European gobiids (Gobiidae, Gobius lineage). Recent evolutionary studies suggest that deterministic and/or stochastic forces could generate acoustic differences among related animal species, though this has not been investigated in any teleost group to date. In the present comparative study, we analysed the sounds from nine soniferous gobiids and quantitatively assessed their acoustic variability. Our interspecific acoustic study, incorporating for the first time the representative acoustic signals from the majority of soniferous gobiids, suggested that their sounds are truly species-specific (92% of sounds correctly classified into exact species) and each taxon possesses a unique set of spectro-temporal variables. In addition, we reconstructed phylogenetic relationships from a concatenated molecular dataset consisting of multiple molecular markers to track the evolution of acoustic signals in soniferous gobiids. The results of this study indicated that the genus Padogobius is polyphyletic, since P. nigricans was nested within the Ponto-Caspian clade, while the congeneric P. bonelli turned out to be a sister taxon to the remaining investigated soniferous species. Lastly, by extracting the acoustic and genetic distance matrices, sound variability and genetic distance were correlated for the first time to assess whether sound evolution follows a similar phylogenetic pattern. The positive correlation between the sound variability and genetic distance obtained here emphasizes that certain acoustic features from representative sounds could carry the phylogenetic signal in soniferous gobiids. Our study was the first attempt to evaluate the mutual relationship between acoustic variation and genetic divergence in any teleost fish.

Enterococcus faecium are associated with the modification of gut microbiota and shrimp post-larvae survival

Background Probiotics are widely used to promote host health. Compared to mammals and terrestrial invertebrates, little is known the role of probiotics in aquatic invertebrates. In this study, eighteen tanks with eight hundred of shrimp post-larvae individuals each were randomly grouped into three groups, one is shrimps administered with E. faecium as probiotic (Tre) and others are shrimps without probiotic-treatment (CK1: blank control, CK2: medium control). We investigated the correlations between a kind of commercial Enterococcus faecium (E. faecium) powder and microbiota composition with function potentials in shrimp post-larvae gut. Results We sequenced the 16S rRNA gene (V4) of gut samples to assess diversity and composition of the shrimp gut microbiome and used differential abundance and Tax4Fun2 analyses to identify the differences of taxonomy and predicted function between different treatment groups. The ingested probiotic bacteria (E. faecium) were tracked in gut microbiota of Tre and the shrimps here showed the best growth performance especially in survival ratio (SR). The distribution of SR across samples was similar to that in PCoA plot based on Bray-Curits and two subgroups generated (SL: SR < 70%, SH: SR ≥ 70%). The gut microbiota structure and predicted function were correlated with both treatment and SR, and SR was a far more important factor driving taxonomic and functional differences than treatment. Both Tre and SH showed a low and uneven community species and shorted phylogenetic distance. We detected a shift in composition profile at phylum and genus level and further identified ten OTUs as relevant taxa that both closely associated with treatment and SR. The partial least squares path model further supported the important role of relevant taxa related to shrimp survival ratio. Conclusions Overall, we found gut microbiota correlated to both shrimp survival and ingested probiotic bacteria (E. faecium). These correlations should not be dismissed without merit and will uncover a promising strategy for developing novel probiotics through certain consortium of gut microbiota.

The evolution of ovary-specific gene expression in Hawaiian Drosophilidae

As detailed data on gene expression become accessible from more species, we have an opportunity to test the extent to which our understanding of developmental genetics from model organisms helps predict expression patterns across species. Central to this is the question: how much variation in gene expression do we expect to observe between species? Here we provide an answer by comparing RNAseq data between twelve species of Hawaiian Drosophilidae flies, focusing on gene expression differences between the ovary and other tissues. We show that there exists a cohort of ovary-specific genes that is stable across species, and that largely corresponds to described expression patterns from laboratory model Drosophila species. However, our results also show that, as phylogenetic distance increases, variation between species overwhelms variation between tissues. Using ancestral state reconstruction of expression, we describe the distribution of evolutionary changes in tissue-biased expression profiles, and use this to identify gains and losses of ovarian expression across these twelve species. We then use this distribution to calculate the correlation in expression evolution between genes, and demonstrate that genes with known interactions in D. melanogaster are significantly more correlated in their evolution than genes with no or unknown interactions. Finally, we use this correlation matrix to infer new networks of genes that have similar evolutionary trajectories, and we provide these as a dataset of novel testable hypotheses about genetic roles and interactions.

Phylosymbiosis in the Rhizosphere Microbiome Extends to Nitrogen Cycle Functional Potential

Most plants rely on specialized root-associated microbes to obtain essential nitrogen (N), yet not much is known about the evolutionary history of the rhizosphere–plant interaction. We conducted a common garden experiment to investigate the plant root–rhizosphere microbiome association using chloridoid grasses sampled from around the world and grown from seed in a greenhouse. We sought to test whether plants that are more closely related phylogenetically have more similar root bacterial microbiomes than plants that are more distantly related. Using metagenome sequencing, we found that there is a conserved core and a variable rhizosphere bacterial microbiome across the chloridoid grasses. Additionally, phylogenetic distance among the host plant species was correlated with bacterial community composition, suggesting the plant hosts prefer specific bacterial lineages. The functional potential for N utilization across microbiomes fluctuated extensively and mirrored variation in the microbial community composition across host plants. Variation in the bacterial potential for N fixation was strongly affected by the host plants’ phylogeny, whereas variation in N recycling, nitrification, and denitrification was unaffected. This study highlights the evolutionary linkage between the N fixation traits of the microbial community and the plant host and suggests that not all functional traits are equally important for plant–microbe associations.