Substitutions Are Boring: Some Arguments about Parallel Mutations and High Mutation Rates.

Extant genomes are largely shaped by global transposition, copy-number fluctuation, and rearrangement of DNA sequences rather than by substitutions of single nucleotides. Although many of these large-scale mutations have low probabilities and are unlikely to repeat, others are recurrent or predictable in their effects, leading to stereotyped genome architectures and genetic variation in both eukaryotes and prokaryotes. Such recurrent, parallel mutation modes can profoundly shape the paths taken by evolution and undermine common models of evolutionary genetics. Similar patterns are also evident at the smaller scales of individual genes or short sequences. The scale and extent of this 'non-substitution' variation has recently come into focus through the advent of new genomic technologies; however, it is still not widely considered in genotype-phenotype association studies. In this review we identify common features of these disparate mutational phenomena and comment on the importance and interpretation of these mutational patterns.

Litter Commensal Bacteria Can Limit the Horizontal Gene Transfer of Antimicrobial Resistance to Salmonella in Chickens.

Fostering a "balanced" gut microbiome through the administration of beneficial microbes that can competitively exclude pathogens has gained a lot of attention and use in human and animal medicine. However, little is known about how microbes affect the horizontal gene transfer of antimicrobial resistance (AMR). To shed more light on this question, we challenged neonatal broiler chicks raised on reused broiler chicken litter-a complex environment made up of decomposing pine shavings, feces, uric acid, feathers, and feed-with Salmonella enterica serovar Heidelberg (S. Heidelberg), a model pathogen. Neonatal chicks challenged with S. Heidelberg and raised on reused litter were more resistant to S. Heidelberg cecal colonization than chicks grown on fresh litter. Furthermore, chicks grown on reused litter were at a lower risk of colonization with S. Heidelberg strains that encoded AMR on IncI1 plasmids. We used 16S rRNA gene sequencing and shotgun metagenomics to show that the major difference between chicks grown on fresh litter and those grown on reused litter was the microbiome harbored in the litter and ceca. The microbiome of reused litter samples was more uniform and enriched in functional pathways related to the biosynthesis of organic and antimicrobial molecules than that in fresh litter samples. We found that Escherichia coli was the main reservoir of plasmids encoding AMR and that the IncI1 plasmid was maintained at a significantly lower copy per cell in reused litter compared to fresh litter. These findings support the notion that commensal bacteria play an integral role in the horizontal transfer of plasmids encoding AMR to pathogens like Salmonella. IMPORTANCE Antimicrobial resistance spread is a worldwide health challenge, stemming in large part from the ability of microorganisms to share their genetic material through horizontal gene transfer. To address this issue, many countries and international organizations have adopted a One Health approach to curtail the proliferation of antimicrobial-resistant bacteria. This includes the removal and reduction of antibiotics used in food animal production and the development of alternatives to antibiotics. However, there is still a significant knowledge gap in our understanding of how resistance spreads in the absence of antibiotic selection and the role commensal bacteria play in reducing antibiotic resistance transfer. In this study, we show that commensal bacteria play a key role in reducing the horizontal gene transfer of antibiotic resistance to Salmonella, provide the identity of the bacterial species that potentially perform this function in broiler chickens, and also postulate the mechanism involved.

Massive variation of short tandem repeats with functional consequences across strains of Arabidopsis thaliana.

Short tandem repeat (STR) mutations may comprise more than half of the mutations in eukaryotic coding DNA, yet STR variation is rarely examined as a contributor to complex traits. We assessed this contribution across a collection of 96 strains of Arabidopsis thaliana, genotyping 2046 STR loci each, using highly parallel STR sequencing with molecular inversion probes. We found that 95% of examined STRs are polymorphic, with a median of six alleles per STR across these strains. STR expansions (large copy number increases) are found in most strains, several of which have evident functional effects. These include three of six intronic STR expansions we found to be associated with intron retention. Coding STRs were depleted of variation relative to noncoding STRs, and we detected a total of 56 coding STRs (11%) showing low variation consistent with the action of purifying selection. In contrast, some STRs show hypervariable patterns consistent with diversifying selection. Finally, we detected 133 novel STR-phenotype associations under stringent criteria, most of which could not be detected with SNPs alone, and validated some with follow-up experiments. Our results support the conclusion that STRs constitute a large, unascertained reservoir of functionally relevant genomic variation.

Evolutionary assembly patterns of prokaryotic genomes.

Evolutionary innovation must occur in the context of some genomic background, which limits available evolutionary paths. For example, protein evolution by sequence substitution is constrained by epistasis between residues. In prokaryotes, evolutionary innovation frequently happens by macrogenomic events such as horizontal gene transfer (HGT). Previous work has suggested that HGT can be influenced by ancestral genomic content, yet the extent of such gene-level constraints has not yet been systematically characterized. Here, we evaluated the evolutionary impact of such constraints in prokaryotes, using probabilistic ancestral reconstructions from 634 extant prokaryotic genomes and a novel framework for detecting evolutionary constraints on HGT events. We identified 8228 directional dependencies between genes and demonstrated that many such dependencies reflect known functional relationships, including for example, evolutionary dependencies of the photosynthetic enzyme RuBisCO. Modeling all dependencies as a network, we adapted an approach from graph theory to establish chronological precedence in the acquisition of different genomic functions. Specifically, we demonstrated that specific functions tend to be gained sequentially, suggesting that evolution in prokaryotes is governed by functional assembly patterns. Finally, we showed that these dependencies are universal rather than clade-specific and are often sufficient for predicting whether or not a given ancestral genome will acquire specific genes. Combined, our results indicate that evolutionary innovation via HGT is profoundly constrained by epistasis and historical contingency, similar to the evolution of proteins and phenotypic characters, and suggest that the emergence of specific metabolic and pathological phenotypes in prokaryotes can be predictable from current genomes.

MIPSTR: a method for multiplex genotyping of germline and somatic STR variation across many individuals.

Short tandem repeats (STRs) are highly mutable genetic elements that often reside in regulatory and coding DNA. The cumulative evidence of genetic studies on individual STRs suggests that STR variation profoundly affects phenotype and contributes to trait heritability. Despite recent advances in sequencing technology, STR variation has remained largely inaccessible across many individuals compared to single nucleotide variation or copy number variation. STR genotyping with short-read sequence data is confounded by (1) the difficulty of uniquely mapping short, low-complexity reads; and (2) the high rate of STR amplification stutter. Here, we present MIPSTR, a robust, scalable, and affordable method that addresses these challenges. MIPSTR uses targeted capture of STR loci by single-molecule Molecular Inversion Probes (smMIPs) and a unique mapping strategy. Targeted capture and our mapping strategy resolve the first challenge; the use of single molecule information resolves the second challenge. Unlike previous methods, MIPSTR is capable of distinguishing technical error due to amplification stutter from somatic STR mutations. In proof-of-principle experiments, we use MIPSTR to determine germline STR genotypes for 102 STR loci with high accuracy across diverse populations of the plant A. thaliana. We show that putatively functional STRs may be identified by deviation from predicted STR variation and by association with quantitative phenotypes. Using DNA mixing experiments and a mutant deficient in DNA repair, we demonstrate that MIPSTR can detect low-frequency somatic STR variants. MIPSTR is applicable to any organism with a high-quality reference genome and is scalable to genotyping many thousands of STR loci in thousands of individuals.

An E3 ubiquitin ligase prevents ectopic localization of the centromeric histone H3 variant via the centromere targeting domain.

Proper centromere function is critical to maintain genomic stability and to prevent aneuploidy, a hallmark of tumors and birth defects. A conserved feature of all eukaryotic centromeres is an essential histone H3 variant called CENP-A that requires a centromere targeting domain (CATD) for its localization. Although proteolysis prevents CENP-A from mislocalizing to euchromatin, regulatory factors have not been identified. Here, we identify an E3 ubiquitin ligase called Psh1 that leads to the degradation of Cse4, the budding yeast CENP-A homolog. Cse4 overexpression is toxic to psh1Δ cells and results in euchromatic localization. Strikingly, the Cse4 CATD is a key regulator of its stability and helps Psh1 discriminate Cse4 from histone H3. Taken together, we propose that the CATD has a previously unknown role in maintaining the exclusive localization of Cse4 by preventing its mislocalization to euchromatin via Psh1-mediated degradation.

The overdue promise of short tandem repeat variation for heritability.

Short tandem repeat (STR) variation has been proposed as a major explanatory factor in the heritability of complex traits in humans and model organisms. However, we still struggle to incorporate STR variation into genotype-phenotype maps. We review here the promise of STRs in contributing to complex trait heritability and highlight the challenges that STRs pose due to their repetitive nature. We argue that STR variants are more likely than single-nucleotide variants to have epistatic interactions, reiterate the need for targeted assays to genotype STRs accurately, and call for more appropriate statistical methods in detecting STR-phenotype associations. Lastly, we suggest that somatic STR variation within individuals may serve as a read-out of disease susceptibility, and is thus potentially a valuable covariate for future association studies.

Genome-scale co-evolutionary inference identifies functions and clients of bacterial Hsp90.

The molecular chaperone Hsp90 is essential in eukaryotes, in which it facilitates the folding of developmental regulators and signal transduction proteins known as Hsp90 clients. In contrast, Hsp90 is not essential in bacteria, and a broad characterization of its molecular and organismal function is lacking. To enable such characterization, we used a genome-scale phylogenetic analysis to identify genes that co-evolve with bacterial Hsp90. We find that genes whose gain and loss were coordinated with Hsp90 throughout bacterial evolution tended to function in flagellar assembly, chemotaxis, and bacterial secretion, suggesting that Hsp90 may aid assembly of protein complexes. To add to the limited set of known bacterial Hsp90 clients, we further developed a statistical method to predict putative clients. We validated our predictions by demonstrating that the flagellar protein FliN and the chemotaxis kinase CheA behaved as Hsp90 clients in Escherichia coli, confirming the predicted role of Hsp90 in chemotaxis and flagellar assembly. Furthermore, normal Hsp90 function is important for wild-type motility and/or chemotaxis in E. coli. This novel function of bacterial Hsp90 agreed with our subsequent finding that Hsp90 is associated with a preference for multiple habitats and may therefore face a complex selection regime. Taken together, our results reveal previously unknown functions of bacterial Hsp90 and open avenues for future experimental exploration by implicating Hsp90 in the assembly of membrane protein complexes and adaptation to novel environments.

Background-dependent effects of polyglutamine variation in the Arabidopsis thaliana gene ELF3.

Tandem repeats (TRs) have extremely high mutation rates and are often considered to be neutrally evolving DNA. However, in coding regions, TR copy number mutations can significantly affect phenotype and may facilitate rapid adaptation to new environments. In several human genes, TR copy number mutations that expand polyglutamine (polyQ) tracts beyond a certain threshold cause incurable neurodegenerative diseases. PolyQ-containing proteins exist at a considerable frequency in eukaryotes, yet the phenotypic consequences of natural variation in polyQ tracts that are not associated with disease remain largely unknown. Here, we use Arabidopsis thaliana to dissect the phenotypic consequences of natural variation in the polyQ tract encoded by EARLY FLOWERING 3 (ELF3), a key developmental gene. Changing ELF3 polyQ tract length affected complex ELF3-dependent phenotypes in a striking and nonlinear manner. Some natural ELF3 polyQ variants phenocopied elf3 loss-of-function mutants in a common reference background, although they are functional in their native genetic backgrounds. To test the existence of background-specific modifiers, we compared the phenotypic effects of ELF3 polyQ variants between two divergent backgrounds, Col and Ws, and found dramatic differences. In fact, the Col-ELF3 allele, encoding the shortest known ELF3 polyQ tract, was haploinsufficient in Ws × Col F(1) hybrids. Our data support a model in which variable polyQ tracts drive adaptation to internal genetic environments.

High-throughput identification of viral termini and packaging mechanisms in virome datasets using PhageTermVirome.

Viruses that infect bacteria (phages) are increasingly recognized for their importance in diverse ecosystems but identifying and annotating them in large-scale sequence datasets is still challenging. Although efficient scalable virus identification tools are emerging, defining the exact ends (termini) of phage genomes is still particularly difficult. The proper identification of termini is crucial, as it helps in characterizing the packaging mechanism of bacteriophages and provides information on various aspects of phage biology. Here, we introduce PhageTermVirome (PTV) as a tool for the easy and rapid high-throughput determination of phage termini and packaging mechanisms using modern large-scale metagenomics datasets. We successfully tested the PTV algorithm on a mock virome dataset and then used it on two real virome datasets to achieve the rapid identification of more than 100 phage termini and packaging mechanisms, with just a few hours of computing time. Because PTV allows the identification of free fully formed viral particles (by recognition of termini present only in encapsidated DNA), it can also complement other virus identification softwares to predict the true viral origin of contigs in viral metagenomics datasets. PTV is a novel and unique tool for high-throughput characterization of phage genomes, including phage termini identification and characterization of genome packaging mechanisms. This software should help researchers better visualize, map and study the virosphere. PTV is freely available for downloading and installation at https://gitlab.pasteur.fr/vlegrand/ptv .

Horizontal Gene Transfer Is the Main Driver of Antimicrobial Resistance in Broiler Chicks Infected with Salmonella enterica Serovar Heidelberg.

The overuse and misuse of antibiotics in clinical settings and in food production have been linked to the increased prevalence and spread of antimicrobial resistance (AR). Consequently, public health and consumer concerns have resulted in a remarkable reduction in antibiotics used for food animal production. However, there are no data on the effectiveness of antibiotic removal in reducing AR shared through horizontal gene transfer (HGT). In this study, we used neonatal broiler chicks and Salmonella enterica serovar Heidelberg, a model food pathogen, to test if chicks raised antibiotic free harbor transferable AR. We challenged chicks with an antibiotic-susceptible S. Heidelberg strain using various routes of inoculation and determined if S. Heidelberg isolates recovered carried plasmids conferring AR. We used antimicrobial susceptibility testing and whole-genome sequencing (WGS) to show that chicks grown without antibiotics harbored an antimicrobial resistant S. Heidelberg population at 14 days after challenge and chicks challenged orally acquired AR at a higher rate than chicks inoculated via the cloaca. Using 16S rRNA gene sequencing, we found that S. Heidelberg infection perturbed the microbiota of broiler chicks, and we used metagenomics and WGS to confirm that a commensal Escherichia coli population was the main reservoir of an IncI1 plasmid acquired by S. Heidelberg. The carriage of this IncI1 plasmid posed no fitness cost to S. Heidelberg but increased its fitness when exposed to acidic pH in vitro. These results suggest that HGT of plasmids carrying AR shaped the evolution of S. Heidelberg and that antibiotic use reduction alone is insufficient to limit antibiotic resistance transfer from commensal bacteria to Salmonella enterica. IMPORTANCE The reported increase in antibiotic-resistant bacteria in humans has resulted in a major shift away from antibiotic use in food animal production. This shift has been driven by the assumption that removing antibiotics will select for antibiotic susceptible bacterial taxa, which in turn will allow the currently available antibiotic arsenal to be more effective. This change in practice has highlighted new questions that need to be answered to assess the effectiveness of antibiotic removal in reducing the spread of antibiotic resistance bacteria. This research demonstrates that antibiotic-susceptible Salmonella enterica serovar Heidelberg strains can acquire multidrug resistance from commensal bacteria present in the gut of neonatal broiler chicks, even in the absence of antibiotic selection. We demonstrate that exposure to acidic pH drove the horizontal transfer of antimicrobial resistance plasmids and suggest that simply removing antibiotics from food animal production might not be sufficient to limit the spread of antimicrobial resistance.

High-Quality Assemblies for Three Invasive Social Wasps from the Vespula Genus.

Social wasps of the genus Vespula have spread to nearly all landmasses worldwide and have become significant pests in their introduced ranges, affecting economies and biodiversity. Comprehensive genome assemblies and annotations for these species are required to develop the next generation of control strategies and monitor existing chemical control. We sequenced and annotated the genomes of the common wasp (Vespula vulgaris), German wasp (Vespula germanica), and the western yellowjacket (Vespula pensylvanica). Our chromosome-level Vespula assemblies each contain 176-179 Mb of total sequence assembled into 25 scaffolds, with 10-200 unanchored scaffolds, and 16,566-18,948 genes. We annotated gene sets relevant to the applied management of invasive wasp populations, including genes associated with spermatogenesis and development, pesticide resistance, olfactory receptors, immunity and venom. These genomes provide evidence for active DNA methylation in Vespidae and tandem duplications of venom genes. Our genomic resources will contribute to the development of next-generation control strategies, and monitoring potential resistance to chemical control.

Linking the resistome and plasmidome to the microbiome.

The rapid spread of antibiotic resistance among bacterial pathogens is a serious human health threat. While a range of environments have been identified as reservoirs of antibiotic resistance genes (ARGs), we lack understanding of the origins of these ARGs and their spread from environment to clinic. This is partly due to our inability to identify the natural bacterial hosts of ARGs and the mobile genetic elements that mediate this spread, such as plasmids and integrons. Here we demonstrate that the in vivo proximity-ligation method Hi-C can reconstruct a known plasmid-host association from a wastewater community, and identify the in situ host range of ARGs, plasmids, and integrons by physically linking them to their host chromosomes. Hi-C detected both previously known and novel associations between ARGs, mobile genetic elements and host genomes, thus validating this method. We showed that IncQ plasmids and class 1 integrons had the broadest host range in this wastewater, and identified bacteria belonging to Moraxellaceae, Bacteroides, and Prevotella, and especially Aeromonadaceae as the most likely reservoirs of ARGs in this community. A better identification of the natural carriers of ARGs will aid the development of strategies to limit resistance spread to pathogens.

Degradation of Recalcitrant Polyurethane and Xenobiotic Additives by a Selected Landfill Microbial Community and Its Biodegradative Potential Revealed by Proximity Ligation-Based Metagenomic Analysis.

Polyurethanes (PU) are the sixth most produced plastics with around 18-million tons in 2016, but since they are not recyclable, they are burned or landfilled, generating damage to human health and ecosystems. To elucidate the mechanisms that landfill microbial communities perform to attack recalcitrant PU plastics, we studied the degradative activity of a mixed microbial culture, selected from a municipal landfill by its capability to grow in a water PU dispersion (WPUD) as the only carbon source, as a model for the BP8 landfill microbial community. The WPUD contains a polyether-polyurethane-acrylate (PE-PU-A) copolymer and xenobiotic additives (N-methylpyrrolidone, isopropanol and glycol ethers). To identify the changes that the BP8 microbial community culture generates to the WPUD additives and copolymer, we performed chemical and physical analyses of the biodegradation process during 25 days of cultivation. These analyses included Nuclear magnetic resonance, Fourier transform infrared spectroscopy, Thermogravimetry, Differential scanning calorimetry, Gel permeation chromatography, and Gas chromatography coupled to mass spectrometry techniques. Moreover, for revealing the BP8 community structure and its genetically encoded potential biodegradative capability we also performed a proximity ligation-based metagenomic analysis. The additives present in the WPUD were consumed early whereas the copolymer was cleaved throughout the 25-days of incubation. The analysis of the biodegradation process and the identified biodegradation products showed that BP8 cleaves esters, C-C, and the recalcitrant aromatic urethanes and ether groups by hydrolytic and oxidative mechanisms, both in the soft and the hard segments of the copolymer. The proximity ligation-based metagenomic analysis allowed the reconstruction of five genomes, three of them from novel species. In the metagenome, genes encoding known enzymes, and putative enzymes and metabolic pathways accounting for the biodegradative activity of the BP8 community over the additives and PE-PU-A copolymer were identified. This is the first study revealing the genetically encoded potential biodegradative capability of a microbial community selected from a landfill, that thrives within a WPUD system and shows potential for bioremediation of polyurethane- and xenobiotic additives-contamitated sites.

Assignment of virus and antimicrobial resistance genes to microbial hosts in a complex microbial community by combined long-read assembly and proximity ligation.

We describe a method that adds long-read sequencing to a mix of technologies used to assemble a highly complex cattle rumen microbial community, and provide a comparison to short read-based methods. Long-read alignments and Hi-C linkage between contigs support the identification of 188 novel virus-host associations and the determination of phage life cycle states in the rumen microbial community. The long-read assembly also identifies 94 antimicrobial resistance genes, compared to only seven alleles in the short-read assembly. We demonstrate novel techniques that work synergistically to improve characterization of biological features in a highly complex rumen microbial community.

Assembly of 913 microbial genomes from metagenomic sequencing of the cow rumen.

The cow rumen is adapted for the breakdown of plant material into energy and nutrients, a task largely performed by enzymes encoded by the rumen microbiome. Here we present 913 draft bacterial and archaeal genomes assembled from over 800 Gb of rumen metagenomic sequence data derived from 43 Scottish cattle, using both metagenomic binning and Hi-C-based proximity-guided assembly. Most of these genomes represent previously unsequenced strains and species. The draft genomes contain over 69,000 proteins predicted to be involved in carbohydrate metabolism, over 90% of which do not have a good match in public databases. Inclusion of the 913 genomes presented here improves metagenomic read classification by sevenfold against our own data, and by fivefold against other publicly available rumen datasets. Thus, our dataset substantially improves the coverage of rumen microbial genomes in the public databases and represents a valuable resource for biomass-degrading enzyme discovery and studies of the rumen microbiome.

Variability in a Short Tandem Repeat Mediates Complex Epistatic Interactions in Arabidopsis thaliana.

Short tandem repeats (STRs) are hypervariable genetic elements that occur frequently in coding regions. Their high mutation rate readily generates genetic variation, contributing to adaptive evolution and human diseases. We previously reported that natural ELF3 polyglutamine variants cause reciprocal genetic incompatibilities in two divergent Arabidopsis thaliana backgrounds. Here, we dissect the genetic architecture of this incompatibility, revealing as many as four loci putatively interacting with ELF3 We were able to specifically identify one such ELF3-interacting gene, LSH9 We further used a yeast two-hybrid strategy to identify proteins whose physical interactions with ELF3 were affected by polyglutamine tract length. We found two proteins for which this was the case, ELF4 and AtGLDP1. Using these two approaches, we identify specific genetic interactions and physical mechanisms by which the ELF3 polyglutamine tract may mediate the observed genetic incompatibilities. Our work elucidates how STR variation, which is generally underascertained in population-scale sequencing, can contribute to phenotypic variation. Furthermore, our results support our proposal that highly variable STR loci can contribute to the epistatic component of heritability.

PIF4 and ELF3 Act Independently in Arabidopsis thaliana Thermoresponsive Flowering.

Plants have evolved elaborate mechanisms controlling developmental responses to environmental stimuli. A particularly important stimulus is temperature. Previous work has identified the interplay of PIF4 and ELF3 as a central circuit underlying thermal responses in Arabidopsis thaliana. However, thermal responses vary widely among strains, possibly offering mechanistic insights into the wiring of this circuit. ELF3 contains a polyglutamine (polyQ) tract that is crucial for ELF3 function and varies in length across strains. Here, we use transgenic analysis to test the hypothesis that natural polyQ variation in ELF3 is associated with the observed natural variation in thermomorphogenesis. We found little evidence that the polyQ tract plays a specific role in thermal responses beyond modulating general ELF3 function. Instead, we made the serendipitous discovery that ELF3 plays a crucial, PIF4-independent role in thermoresponsive flowering under conditions more likely to reflect field conditions. We present evidence that ELF3 acts through the photoperiodic pathway, pointing to a previously unknown symmetry between low and high ambient temperature responses. Moreover, in analyzing two strain backgrounds with different thermal responses, we demonstrate that responses may be shifted rather than fundamentally rewired across strains. Our findings tie together disparate observations into a coherent framework in which multiple pathways converge in accelerating flowering in response to temperature, with some such pathways modulated by photoperiod.