The dynamics of alternative pathways to compensatory substitution

Phylogenies of the 16S rRNA gene and its hypervariable regions lack concordance with core genome phylogenies

BACKGROUND

The 16S rRNA gene is used extensively in bacterial phylogenetics, in species delineation, and now widely in microbiome studies. However, the gene suffers from intragenomic heterogeneity, and reports of recombination and an unreliable phylogenetic signal are accumulating. Here, we compare core gene phylogenies to phylogenies constructed using core gene concatenations to estimate the strength of signal for the 16S rRNA gene, its hypervariable regions, and all core genes at the intra- and inter-genus levels. Specifically, we perform four intra-genus analyses (Clostridium, n = 65; Legionella, n = 47; Staphylococcus, n = 36; and Campylobacter, n = 17) and one inter-genus analysis [41 core genera of the human gut microbiome (31 families, 17 orders, and 12 classes), n = 82].

RESULTS

At both taxonomic levels, the 16S rRNA gene was recombinant and subject to horizontal gene transfer. At the intra-genus level, the gene showed one of the lowest levels of concordance with the core genome phylogeny (50.7% average). Concordance for hypervariable regions was lower still, with entropy masking providing little to no benefit. A major factor influencing concordance was SNP count, which showed a positive logarithmic association. Using this relationship, we determined that 690 ± 110 SNPs were required for 80% concordance (average 16S rRNA gene SNP count was 254). We also found a wide range in 16S-23S-5S rRNA operon copy number among genomes (1-27). At the inter-genus level, concordance for the whole 16S rRNA gene was markedly higher (73.8% - 10th out of 49 loci); however, the most concordant hypervariable regions (V4, V3-V4, and V1-V2) ranked in the third quartile (62.5 to 60.0%).

CONCLUSIONS

Ramifications of a poor phylogenetic performance for the 16S rRNA gene are far reaching. For example, in addition to incorrect species/strain delineation and phylogenetic inference, it has the potential to confound community diversity metrics if phylogenetic information is incorporated - for example, with popular approaches such as Faith's phylogenetic diversity and UniFrac. Our results highlight the problematic nature of these approaches and their use (along with entropy masking) is discouraged. Lastly, the wide range in 16S rRNA gene copy number among genomes also has a strong potential to confound diversity metrics. Video Abstract.

Compensatory Mutations in GI and ZTL May Modulate Temperature Compensation in the Circadian Clock

Circadian systems share the three properties of entrainment, free-running period, and temperature compensation (TC). TC ensures nearly the same period over a broad range of physiologically relevant temperatures; however, the mechanisms behind TC remain poorly understood. Here, we identify single point mutations in two key elements of the Arabidopsis circadian clock, GIGANTEA (GI) and ZEITLUPE (ZTL), which likely act as compensatory substitutions to establish a remarkably constant free-running period over a wide range of temperatures. Using near-isogenic lines generated from the introgression of the Cape Verde Islands (Cvi) alleles of GI and ZTL into the Landsberg erecta (Ler) background, we show how longer periods in the Cvi background at higher temperatures correlate with a difference in strength of the GI/ZTL interaction. Pairwise interaction testing of all GI/ZTL allelic combinations shows similar affinities for isogenic alleles at 22°C, but very poor interaction between GI (Cvi) and ZTL (Cvi) at higher temperature. In vivo, this would result in lower ZTL levels at high temperatures leading to longer periods in the Cvi background. Mismatched allelic combinations result in extremely strong or weak GI/ZTL interactions, indicating how the corresponding natural variants likely became fixed through epistatic selection. Additionally, molecular characterization of GI (Cvi) reveals a novel functional motif that can modulate the GI/ZTL interaction as well as nucleocytoplasmic partitioning. Taken together, these results identify a plausible temperature-dependent molecular mechanism, which contributes to the robustness of TC through natural variation in GI and ZTL alleles found on the Cape Verde Islands.

Reconstruction of ancestral RNA sequences under multiple structural constraints

BACKGROUND

Secondary structures form the scaffold of multiple sequence alignment of non-coding RNA (ncRNA) families. An accurate reconstruction of ancestral ncRNAs must use this structural signal. However, the inference of ancestors of a single ncRNA family with a single consensus structure may bias the results towards sequences with high affinity to this structure, which are far from the true ancestors.

METHODS

In this paper, we introduce achARNement, a maximum parsimony approach that, given two alignments of homologous ncRNA families with consensus secondary structures and a phylogenetic tree, simultaneously calculates ancestral RNA sequences for these two families.

RESULTS

We test our methodology on simulated data sets, and show that achARNement outperforms classical maximum parsimony approaches in terms of accuracy, but also reduces by several orders of magnitude the number of candidate sequences. To conclude this study, we apply our algorithms on the Glm clan and the FinP-traJ clan from the Rfam database.

CONCLUSIONS

Our results show that our methods reconstruct small sets of high-quality candidate ancestors with better agreement to the two target structures than with classical approaches. Our program is freely available at: http://csb.cs.mcgill.ca/acharnement .