Longitudinal multi-omics of host-microbe dynamics in prediabetes.

Type 2 diabetes mellitus (T2D) is a growing health problem, but little is known about its early disease stages, its effects on biological processes or the transition to clinical T2D. To understand the earliest stages of T2D better, we obtained samples from 106 healthy individuals and individuals with prediabetes over approximately four years and performed deep profiling of transcriptomes, metabolomes, cytokines, and proteomes, as well as changes in the microbiome. This rich longitudinal data set revealed many insights: first, healthy profiles are distinct among individuals while displaying diverse patterns of intra- and/or inter-personal variability. Second, extensive host and microbial changes occur during respiratory viral infections and immunization, and immunization triggers potentially protective responses that are distinct from responses to respiratory viral infections. Moreover, during respiratory viral infections, insulin-resistant participants respond differently than insulin-sensitive participants. Third, global co-association analyses among the thousands of profiled molecules reveal specific host-microbe interactions that differ between insulin-resistant and insulin-sensitive individuals. Last, we identified early personal molecular signatures in one individual that preceded the onset of T2D, including the inflammation markers interleukin-1 receptor agonist (IL-1RA) and high-sensitivity C-reactive protein (CRP) paired with xenobiotic-induced immune signalling. Our study reveals insights into pathways and responses that differ between glucose-dysregulated and healthy individuals during health and disease and provides an open-access data resource to enable further research into healthy, prediabetic and T2D states.

A serum-induced gene signature in hepatocytes is associated with pediatric nonalcoholic fatty liver disease.

OBJECTIVE: Pediatric nonalcoholic fatty liver disease (NAFLD) is a growing problem, but its underlying mechanisms are poorly understood. We used transcriptomic reporter cell assays to investigate differences in transcriptional signatures induced in hepatocyte reporter cells by the sera of children with and without NAFLD. METHODS: We studied serum samples from 45 children with NAFLD and 28 children without NAFLD. The sera were used to induce gene expression in cultured HepaRG cells and RNA-sequencing was used to determine gene expression. Computational techniques were used to compare gene expression patterns. RESULTS: Sera from children with NAFLD induced the expression of 195 genes that were significantly differentially expressed in hepatocytes compared to controls with obesity. NAFLD was associated with increased expression of genes promoting inflammation, collagen synthesis, and extracellular matrix remodeling. Additionally, there was lower expression of genes involved in endobiotic and xenobiotic metabolism, and downregulation of peroxisome function, oxidative phosphorylation, and xenobiotic, bile acid, and fatty acid metabolism. A 13-gene signature, including upregulation of TREM1 and MMP1 and downregulation of CYP2C9, was consistently associated with all diagnostic categories of pediatric NAFLD. CONCLUSION: The extracellular milieu of sera from children with NAFLD induced specific gene profiles distinguishable by a hepatocyte reporter system. Circulating factors may contribute to inflammation and extracellular matrix remodeling and impair xenobiotic and endobiotic metabolism in pediatric NAFLD.

Improved DNA Extraction and Amplification Strategy for 16S rRNA Gene Amplicon-Based Microbiome Studies.

Next-generation sequencing technology has driven the rapid advancement of human microbiome studies by enabling community-level sequence profiling of microbiomes. Although all microbiome sequencing methods depend on recovering the DNA from a sample as a first critical step, lysis methods can be a major determinant of microbiome profile bias. Gentle enzyme-based DNA preparation methods preserve DNA quality but can bias the results by failing to open difficult-to-lyse bacteria. Mechanical methods like bead beating can also bias DNA recovery because the mechanical energy required to break tougher cell walls may shear the DNA of the more easily lysed microbes, and shearing can vary depending on the time and intensity of beating, influencing reproducibility. We introduce a non-mechanical, non-enzymatic, novel rapid microbial DNA extraction procedure suitable for 16S rRNA gene-based microbiome profiling applications that eliminates bead beating. The simultaneous application of alkaline, heat, and detergent ('Rapid' protocol) to milligram quantity samples provided consistent representation across the population of difficult and easily lysed bacteria equal to or better than existing protocols, producing sufficient high-quality DNA for full-length 16S rRNA gene PCR. The novel 'Rapid' method was evaluated using mock bacterial communities containing both difficult and easily lysed bacteria. Human fecal sample testing compared the novel Rapid method with a standard Human Microbiome Project (HMP) protocol for samples from lung cancer patients and controls. DNA recovered from both methods was analyzed using 16S rRNA gene sequencing of the V1V3 and V4 regions on the Illumina platform and the V1V9 region on the PacBio platform. Our findings indicate that the 'Rapid' protocol consistently yielded higher levels of Firmicutes species, which reflected the profile of the bacterial community structure more accurately, which was confirmed by mock community evaluation. The novel 'Rapid' DNA lysis protocol reduces population bias common to bead beating and enzymatic lysis methods, presenting opportunities for improved microbial community profiling, combined with the reduction in sample input to 10 milligrams or less, and it enables rapid transfer and simultaneous lysis of 96 samples in a standard plate format. This results in a 20-fold reduction in sample handling time and an overall 2-fold time advantage when compared to widely used commercial methods. We conclude that the novel 'Rapid' DNA extraction protocol offers a reliable alternative for preparing fecal specimens for 16S rRNA gene amplicon sequencing.

Metagenomics and chemotherapy-induced nausea: A roadmap for future research.

Uncontrolled chemotherapy-induced nausea and vomiting can reduce patients' quality of life and may result in premature discontinuation of chemotherapy. Although nausea and vomiting are commonly grouped together, research has shown that antiemetics are clinically effective against chemotherapy-induced vomiting (CIV) but less so against chemotherapy-induced nausea (CIN). Nausea remains a problem for up to 68% of patients who are prescribed guideline-consistent antiemetics. Despite the high prevalence of CIN, relatively little is known regarding its etiology independent of CIV. This review summarizes a metagenomics approach to the study and treatment of CIN with the goal of encouraging future research. Metagenomics focuses on genetic risk factors and encompasses both human (ie, host) and gut microbial genetic variation. Little work to date has focused on metagenomics as a putative biological mechanism of CIN. Metagenomics has the potential to be a powerful tool in advancing scientific understanding of CIN by identifying new biological pathways and intervention targets. The investigation of metagenomics in the context of well-established demographic, clinical, and patient-reported risk factors may help to identify patients at risk and facilitate the prevention and management of CIN.

Enteric defensins are essential regulators of intestinal microbial ecology.

Antimicrobial peptides are important effectors of innate immunity throughout the plant and animal kingdoms. In the mammalian small intestine, Paneth cell alpha-defensins are antimicrobial peptides that contribute to host defense against enteric pathogens. To determine if alpha-defensins also govern intestinal microbial ecology, we analyzed the intestinal microbiota of mice expressing a human alpha-defensin gene (DEFA5) and in mice lacking an enzyme required for the processing of mouse alpha-defensins. In these complementary models, we detected significant alpha-defensin-dependent changes in microbiota composition, but not in total bacterial numbers. Furthermore, DEFA5-expressing mice had striking losses of segmented filamentous bacteria and fewer interleukin 17 (IL-17)-producing lamina propria T cells. Our data ascribe a new homeostatic role to alpha-defensins in regulating the makeup of the commensal microbiota.

HCoDES reveals chromosomal DNA end structures with single-nucleotide resolution.

The structure of broken DNA ends is a critical determinant of the pathway used for DNA double-strand break (DSB) repair. Here, we develop an approach involving the hairpin capture of DNA end structures (HCoDES), which elucidates chromosomal DNA end structures at single-nucleotide resolution. HCoDES defines structures of physiologic DSBs generated by the RAG endonuclease, as well as those generated by nucleases widely used for genome editing. Analysis of G1 phase cells deficient in H2AX or 53BP1 reveals DNA ends that are frequently resected to form long single-stranded overhangs that can be repaired by mutagenic pathways. In addition to 3' overhangs, many of these DNA ends unexpectedly form long 5' single-stranded overhangs. The divergence in DNA end structures resolved by HCoDES suggests that H2AX and 53BP1 may have distinct activities in end protection. Thus, the high-resolution end structures obtained by HCoDES identify features of DNA end processing during DSB repair.

Alterations in subgingival microbiota during full-fixed appliance orthodontic treatment-A prospective study.

BACKGROUND: Full-fixed appliance orthodontic treatment (commonly called braces) increases plaque accumulation and the risk of gingivitis and periodontitis. However, little consensus exists on changes to subgingival microbiota and specific periodontopathogens during treatment with braces. Prior studies have been hampered by selection biases due to dependence on culture conditions, candidate-based PCR and shallow sequencing methods. OBJECTIVE: The objective was to provide the first longitudinal, culture-free and deep-sequence profiling of subgingival bacteria in subjects during early stages of full-fixed orthodontic treatment. METHODS: We performed 16S rRNA next-generation sequencing (NGS) on 168 subgingival samples collected at 4 distinct mandibular tooth sites per subject before (0 weeks) and during (6 and 12 weeks) orthodontic intervention in 9 experimental and 5 control subjects not undergoing treatment. RESULTS: Overall, we noted that orthodontic intervention led to increased microbial richness, accompanied by an increased incidence of localized gingivitis/mild periodontitis in subjects requiring orthodontic treatment compared to controls, as well as significant baseline variations in subgingival microbiomes in all subjects. Moreover, we confirmed individual- and site-dependent microbiome variability (in particular, the lingual site harboured higher microbiome diversity than buccal sites) that orthodontic bands may lead to more prolonged shifts in microbial changes compared to brackets, and evidence of adaptive enrichment of consensus bacteria with orthodontic intervention (12 novel, consensus bacterial species were identified). CONCLUSION: Our study, along with evolving global profiling methods and data analyses, builds a strong foundation for further analyses of subgingival microbiomes during full-fixed orthodontic treatment.

Host genetic control of gut microbiome composition.

The gut microbiome plays a significant role in health and disease, and there is mounting evidence indicating that the microbial composition is regulated in part by host genetics. Heritability estimates for microbial abundance in mice and humans range from (0.05-0.45), indicating that 5-45% of inter-individual variation can be explained by genetics. Through twin studies, genetic association studies, systems genetics, and genome-wide association studies (GWAS), hundreds of specific host genetic loci have been shown to associate with the abundance of discrete gut microbes. Using genetically engineered knock-out mice, at least 30 specific genes have now been validated as having specific effects on the microbiome. The relationships among of host genetics, microbiome composition, and abundance, and disease is now beginning to be unraveled through experiments designed to test causality. The genetic control of disease and its relationship to the microbiome can manifest in multiple ways. First, a genetic variant may directly cause the disease phenotype, resulting in an altered microbiome as a consequence of the disease phenotype. Second, a genetic variant may alter gene expression in the host, which in turn alters the microbiome, producing the disease phenotype. Finally, the genetic variant may alter the microbiome directly, which can result in the disease phenotype. In order to understand the processes that underlie the onset and progression of certain diseases, future research must take into account the relationship among host genetics, microbiome, and disease phenotype, and the resources needed to study these relationships.

Microbiome Signatures Associated With Steatohepatitis and Moderate to Severe Fibrosis in Children With Nonalcoholic Fatty Liver Disease.

BACKGROUND & AIMS: The intestinal microbiome might affect the development and severity of nonalcoholic fatty liver disease (NAFLD). We analyzed microbiomes of children with and without NAFLD. METHODS: We performed a prospective, observational, cross-sectional study of 87 children (age range, 8-17 years) with biopsy-proven NAFLD and 37 children with obesity without NAFLD (controls). Fecal samples were collected and microbiome composition and functions were assessed using 16S ribosomal RNA amplicon sequencing and metagenomic shotgun sequencing. Microbial taxa were identified using zero-inflated negative binomial modeling. Genes contributing to bacterial pathways were identified using gene set enrichment analysis. RESULTS: Fecal microbiomes of children with NAFLD had lower α-diversity than those of control children (3.32 vs 3.52, P = .016). Fecal microbiomes from children with nonalcoholic steatohepatitis (NASH) had the lowest α-diversity (control, 3.52; NAFLD, 3.36; borderline NASH, 3.37; NASH, 2.97; P = .001). High abundance of Prevotella copri was associated with more severe fibrosis (P = .036). Genes for lipopolysaccharide biosynthesis were enriched in microbiomes from children with NASH (P < .001). Classification and regression tree model with level of alanine aminotransferase and relative abundance of the lipopolysaccharide pathway gene encoding 3-deoxy-d-manno-octulosonate 8-phosphate-phosphatase identified patients with NASH with an area under the receiver operating characteristic curve value of 0.92. Genes involved in flagellar assembly were enriched in the fecal microbiomes of patients with moderate to severe fibrosis (P < .001). Classification and regression tree models based on level of alanine aminotransferase and abundance of genes encoding flagellar biosynthesis protein had good accuracy for identifying case children with moderate to severe fibrosis (area under the receiver operating characteristic curve, 0.87). CONCLUSIONS: In an analysis of fecal microbiomes of children with NAFLD, we associated NAFLD and NASH with intestinal dysbiosis. NAFLD and its severity were associated with greater abundance of genes encoding inflammatory bacterial products. Alterations to the intestinal microbiome might contribute to the pathogenesis of NAFLD and be used as markers of disease or severity.

MSC: a metagenomic sequence classification algorithm.

MOTIVATION: Metagenomics is the study of genetic materials directly sampled from natural habitats. It has the potential to reveal previously hidden diversity of microscopic life largely due to the existence of highly parallel and low-cost next-generation sequencing technology. Conventional approaches align metagenomic reads onto known reference genomes to identify microbes in the sample. Since such a collection of reference genomes is very large, the approach often needs high-end computing machines with large memory which is not often available to researchers. Alternative approaches follow an alignment-free methodology where the presence of a microbe is predicted using the information about the unique k-mers present in the microbial genomes. However, such approaches suffer from high false positives due to trading off the value of k with the computational resources. In this article, we propose a highly efficient metagenomic sequence classification (MSC) algorithm that is a hybrid of both approaches. Instead of aligning reads to the full genomes, MSC aligns reads onto a set of carefully chosen, shorter and highly discriminating model sequences built from the unique k-mers of each of the reference sequences. RESULTS: Microbiome researchers are generally interested in two objectives of a taxonomic classifier: (i) to detect prevalence, i.e. the taxa present in a sample, and (ii) to estimate their relative abundances. MSC is primarily designed to detect prevalence and experimental results show that MSC is indeed a more effective and efficient algorithm compared to the other state-of-the-art algorithms in terms of accuracy, memory and runtime. Moreover, MSC outputs an approximate estimate of the abundances. AVAILABILITY AND IMPLEMENTATION: The implementations are freely available for non-commercial purposes. They can be downloaded from https://drive.google.com/open?id=1XirkAamkQ3ltWvI1W1igYQFusp9DHtVl.

Rapid growth is a dominant predictor of hepcidin suppression and declining ferritin in Gambian infants.

Iron deficiency and iron deficiency anemia are highly prevalent in low-income countries, especially among young children. Hepcidin is the major regulator of systemic iron homeostasis. It controls dietary iron absorption, dictates whether absorbed iron is made available in circulation for erythropoiesis and other iron-demanding processes, and predicts response to oral iron supplementation. Understanding how hepcidin is itself regulated is therefore important, especially in young children. We investigated how changes in iron-related parameters, inflammation and infection status, seasonality, and growth influenced plasma hepcidin and ferritin concentrations during infancy using longitudinal data from two birth cohorts of infants in rural Gambia (n=114 and n=193). This setting is characterized by extreme seasonality, prevalent childhood anemia, undernutrition, and frequent infection. Plasma was collected from infants at birth and at regular intervals, up to 12 months of age. Hepcidin, ferritin and plasma iron concentrations declined markedly during infancy, with reciprocal increases in soluble transferrin receptor and transferrin concentrations, indicating declining iron stores and increasing tissue iron demand. In cross-sectional analyses at 5 and 12 months of age, we identified expected relationships of hepcidin with iron and inflammatory markers, but also observed significant negative associations between hepcidin and antecedent weight gain. Correspondingly, longitudinal fixed effects modeling demonstrated weight gain to be the most notable dynamic predictor of decreasing hepcidin and ferritin through infancy across both cohorts. Infants who grow rapidly in this setting are at particular risk of depletion of iron stores, but since hepcidin concentrations decrease with weight gain, they may also be the most responsive to oral iron interventions.

Pervasive genetic hitchhiking and clonal interference in forty evolving yeast populations.

The dynamics of adaptation determine which mutations fix in a population, and hence how reproducible evolution will be. This is central to understanding the spectra of mutations recovered in the evolution of antibiotic resistance, the response of pathogens to immune selection, and the dynamics of cancer progression. In laboratory evolution experiments, demonstrably beneficial mutations are found repeatedly, but are often accompanied by other mutations with no obvious benefit. Here we use whole-genome whole-population sequencing to examine the dynamics of genome sequence evolution at high temporal resolution in 40 replicate Saccharomyces cerevisiae populations growing in rich medium for 1,000 generations. We find pervasive genetic hitchhiking: multiple mutations arise and move synchronously through the population as mutational 'cohorts'. Multiple clonal cohorts are often present simultaneously, competing with each other in the same population. Our results show that patterns of sequence evolution are driven by a balance between these chance effects of hitchhiking and interference, which increase stochastic variation in evolutionary outcomes, and the deterministic action of selection on individual mutations, which favours parallel evolutionary solutions in replicate populations.

Genomic variation landscape of the human gut microbiome.

Whereas large-scale efforts have rapidly advanced the understanding and practical impact of human genomic variation, the practical impact of variation is largely unexplored in the human microbiome. We therefore developed a framework for metagenomic variation analysis and applied it to 252 faecal metagenomes of 207 individuals from Europe and North America. Using 7.4 billion reads aligned to 101 reference species, we detected 10.3 million single nucleotide polymorphisms (SNPs), 107,991 short insertions/deletions, and 1,051 structural variants. The average ratio of non-synonymous to synonymous polymorphism rates of 0.11 was more variable between gut microbial species than across human hosts. Subjects sampled at varying time intervals exhibited individuality and temporal stability of SNP variation patterns, despite considerable composition changes of their gut microbiota. This indicates that individual-specific strains are not easily replaced and that an individual might have a unique metagenomic genotype, which may be exploitable for personalized diet or drug intake.

Identification of functional variants for cleft lip with or without cleft palate in or near PAX7, FGFR2, and NOG by targeted sequencing of GWAS loci.

Although genome-wide association studies (GWASs) for nonsyndromic orofacial clefts have identified multiple strongly associated regions, the causal variants are unknown. To address this, we selected 13 regions from GWASs and other studies, performed targeted sequencing in 1,409 Asian and European trios, and carried out a series of statistical and functional analyses. Within a cluster of strongly associated common variants near NOG, we found that one, rs227727, disrupts enhancer activity. We furthermore identified significant clusters of non-coding rare variants near NTN1 and NOG and found several rare coding variants likely to affect protein function, including four nonsense variants in ARHGAP29. We confirmed 48 de novo mutations and, based on best biological evidence available, chose two of these for functional assays. One mutation in PAX7 disrupted the DNA binding of the encoded transcription factor in an in vitro assay. The second, a non-coding mutation, disrupted the activity of a neural crest enhancer downstream of FGFR2 both in vitro and in vivo. This targeted sequencing study provides strong functional evidence implicating several specific variants as primary contributory risk alleles for nonsyndromic clefting in humans.

Diet during Pregnancy and Infancy and the Infant Intestinal Microbiome.

OBJECTIVES: To determine the association between diet during pregnancy and infancy, including breastfeeding vs formula feeding, solid food introduction, and the infant intestinal microbiome. STUDY DESIGN: Infants participating in the Vitamin D Antenatal Asthma Reduction Trial were included in this study (n = 323). Maternal and infant diets were assessed by questionnaire. Infant stool samples were collected at age 3-6 months. Stool sequencing was performed using the Roche 454 platform. Analyses were stratified by race/ethnicity. RESULTS: Breastfeeding, compared with formula feeding, was independently associated with infant intestinal microbial diversity. Breastfeeding also had the most consistent associations with individual taxa that have been previously linked to early-life diet and health outcomes (eg, Bifidobacterium). Maternal diet during pregnancy and solid food introduction were less associated with the infant gut microbiome than breastfeeding status. We found evidence of a possible interaction between breastfeeding and child race/ethnicity on microbial composition. CONCLUSIONS: Breastfeeding vs formula feeding is the dietary factor that is most consistently independently associated with the infant intestinal microbiome. The relationship between breastfeeding status and intestinal microbiome composition varies by child race/ethnicity. Future studies will need to investigate factors, including genomic factors, which may influence the response of the microbiome to diet. TRIAL REGISTRATION: ClinicalTrials.gov: NCT00920621.

A prospective microbiome-wide association study of food sensitization and food allergy in early childhood.

BACKGROUND: Alterations in the intestinal microbiome are prospectively associated with the development of asthma; less is known regarding the role of microbiome alterations in food allergy development. METHODS: Intestinal microbiome samples were collected at age 3-6 months in children participating in the follow-up phase of an interventional trial of high-dose vitamin D given during pregnancy. At age 3, sensitization to foods (milk, egg, peanut, soy, wheat, walnut) was assessed. Food allergy was defined as caretaker report of healthcare provider-diagnosed allergy to the above foods prior to age 3 with evidence of IgE sensitization. Analysis was performed using Phyloseq and DESeq2; P-values were adjusted for multiple comparisons. RESULTS: Complete data were available for 225 children; there were 87 cases of food sensitization and 14 cases of food allergy. Microbial diversity measures did not differ between food sensitization and food allergy cases and controls. The genera Haemophilus (log2 fold change -2.15, P=.003), Dialister (log2 fold change -2.22, P=.009), Dorea (log2 fold change -1.65, P=.02), and Clostridium (log2 fold change -1.47, P=.002) were underrepresented among subjects with food sensitization. The genera Citrobacter (log2 fold change -3.41, P=.03), Oscillospira (log2 fold change -2.80, P=.03), Lactococcus (log2 fold change -3.19, P=.05), and Dorea (log2 fold change -3.00, P=.05) were underrepresented among subjects with food allergy. CONCLUSIONS: The temporal association between bacterial colonization and food sensitization and allergy suggests that the microbiome may have a causal role in the development of food allergy. Our findings have therapeutic implications for the prevention and treatment of food allergy.

Genomic approaches to studying the human microbiota.

The human body is colonized by a vast array of microbes, which form communities of bacteria, viruses and microbial eukaryotes that are specific to each anatomical environment. Every community must be studied as a whole because many organisms have never been cultured independently, and this poses formidable challenges. The advent of next-generation DNA sequencing has allowed more sophisticated analysis and sampling of these complex systems by culture-independent methods. These methods are revealing differences in community structure between anatomical sites, between individuals, and between healthy and diseased states, and are transforming our view of human biology.

A multi-megabase copy number gain causes maternal transmission ratio distortion on mouse chromosome 2.

Significant departures from expected Mendelian inheritance ratios (transmission ratio distortion, TRD) are frequently observed in both experimental crosses and natural populations. TRD on mouse Chromosome (Chr) 2 has been reported in multiple experimental crosses, including the Collaborative Cross (CC). Among the eight CC founder inbred strains, we found that Chr 2 TRD was exclusive to females that were heterozygous for the WSB/EiJ allele within a 9.3 Mb region (Chr 2 76.9 - 86.2 Mb). A copy number gain of a 127 kb-long DNA segment (designated as responder to drive, R2d) emerged as the strongest candidate for the causative allele. We mapped R2d sequences to two loci within the candidate interval. R2d1 is located near the proximal boundary, and contains a single copy of R2d in all strains tested. R2d2 maps to a 900 kb interval, and the number of R2d copies varies from zero in classical strains (including the mouse reference genome) to more than 30 in wild-derived strains. Using real-time PCR assays for the copy number, we identified a mutation (R2d2WSBdel1) that eliminates the majority of the R2d2WSB copies without apparent alterations of the surrounding WSB/EiJ haplotype. In a three-generation pedigree segregating for R2d2WSBdel1, the mutation is transmitted to the progeny and Mendelian segregation is restored in females heterozygous for R2d2WSBdel1, thus providing direct evidence that the copy number gain is causal for maternal TRD. We found that transmission ratios in R2d2WSB heterozygous females vary between Mendelian segregation and complete distortion depending on the genetic background, and that TRD is under genetic control of unlinked distorter loci. Although the R2d2WSB transmission ratio was inversely correlated with average litter size, several independent lines of evidence support the contention that female meiotic drive is the cause of the distortion. We discuss the implications and potential applications of this novel meiotic drive system.

An Analysis of the Epidemic of Klebsiella pneumoniae Carbapenemase-Producing K. pneumoniae: Convergence of Two Evolutionary Mechanisms Creates the "Perfect Storm".

Background: Carbapenem resistance is a critical healthcare challenge worldwide. Particularly concerning is the widespread dissemination of Klebsiella pneumoniae carbapenemase (KPC). Klebsiella pneumoniae harboring blaKPC (KPC-Kpn) is endemic in many areas including the United States, where the epidemic was primarily mediated by the clonal dissemination of Kpn ST258. We postulated that the spread of blaKPC in other regions occurs by different and more complex mechanisms. To test this, we investigated the evolution and dynamics of spread of KPC-Kpn in Colombia, where KPC became rapidly endemic after emerging in 2005. Methods: We sequenced the genomes of 133 clinical isolates recovered from 24 tertiary care hospitals located in 10 cities throughout Colombia, between 2002 (before the emergence of KPC-Kpn) and 2014. Phylogenetic reconstructions and evolutionary mapping were performed to determine temporal and genetic associations between the isolates. Results: Our results indicate that the start of the epidemic was driven by horizontal dissemination of mobile genetic elements carrying blaKPC-2, followed by the introduction and subsequent spread of clonal group 258 (CG258) isolates containing blaKPC-3. Conclusions: The combination of 2 evolutionary mechanisms of KPC-Kpn within a challenged health system of a developing country created the "perfect storm" for sustained endemicity of these multidrug-resistant organisms in Colombia.

R2d2 Drives Selfish Sweeps in the House Mouse.

A selective sweep is the result of strong positive selection driving newly occurring or standing genetic variants to fixation, and can dramatically alter the pattern and distribution of allelic diversity in a population. Population-level sequencing data have enabled discoveries of selective sweeps associated with genes involved in recent adaptations in many species. In contrast, much debate but little evidence addresses whether "selfish" genes are capable of fixation-thereby leaving signatures identical to classical selective sweeps-despite being neutral or deleterious to organismal fitness. We previously described R2d2, a large copy-number variant that causes nonrandom segregation of mouse Chromosome 2 in females due to meiotic drive. Here we show population-genetic data consistent with a selfish sweep driven by alleles of R2d2 with high copy number (R2d2(HC)) in natural populations. We replicate this finding in multiple closed breeding populations from six outbred backgrounds segregating for R2d2 alleles. We find that R2d2(HC) rapidly increases in frequency, and in most cases becomes fixed in significantly fewer generations than can be explained by genetic drift. R2d2(HC) is also associated with significantly reduced litter sizes in heterozygous mothers, making it a true selfish allele. Our data provide direct evidence of populations actively undergoing selfish sweeps, and demonstrate that meiotic drive can rapidly alter the genomic landscape in favor of mutations with neutral or even negative effects on overall Darwinian fitness. Further study will reveal the incidence of selfish sweeps, and will elucidate the relative contributions of selfish genes, adaptation and genetic drift to evolution.