Mouse genomic associations with in vitro sensitivity to simulated space radiation.

Exposure to ionizing radiation is considered by NASA to be a major health hazard for deep space exploration missions. Ionizing radiation sensitivity is modulated by both genomic and environmental factors. Understanding their contributions is crucial for designing experiments in model organisms, evaluating the risk of deep space (i.e. high-linear energy transfer, or LET, particle) radiation exposure in astronauts, and also selecting therapeutic irradiation regimes for cancer patients. We identified single nucleotide polymorphisms in 15 strains of mice, including 10 collaborative cross model strains and 5 founder strains, associated with spontaneous and ionizing radiation-induced in vitro DNA damage quantified based on immunofluorescent tumor protein p53 binding protein (53BP1) positive nuclear foci. Statistical analysis suggested an association with pathways primarily related to cellular signaling, metabolism, tumorigenesis and nervous system damage. We observed different genomic associations in early (4 and 8 h) responses to different LET radiation, while later (24 hour) DNA damage responses showed a stronger overlap across all LETs. Furthermore, a subset of pathways was associated with spontaneous DNA damage, suggesting 53BP1 positive foci as a potential biomarker for DNA integrity in mouse models. Our results suggest several mouse strains as new models to further study the impact of ionizing radiation and validate the identified genetic loci. We also highlight the importance of future human in vitro studies to refine the association of genes and pathways with the DNA damage response to ionizing radiation and identify targets for space travel countermeasures.

Complete genomic and epigenetic maps of human centromeres.

Existing human genome assemblies have almost entirely excluded repetitive sequences within and near centromeres, limiting our understanding of their organization, evolution, and functions, which include facilitating proper chromosome segregation. Now, a complete, telomere-to-telomere human genome assembly (T2T-CHM13) has enabled us to comprehensively characterize pericentromeric and centromeric repeats, which constitute 6.2% of the genome (189.9 megabases). Detailed maps of these regions revealed multimegabase structural rearrangements, including in active centromeric repeat arrays. Analysis of centromere-associated sequences uncovered a strong relationship between the position of the centromere and the evolution of the surrounding DNA through layered repeat expansions. Furthermore, comparisons of chromosome X centromeres across a diverse panel of individuals illuminated high degrees of structural, epigenetic, and sequence variation in these complex and rapidly evolving regions.

An integrated host-microbiome response to atrazine exposure mediates toxicity in Drosophila.

The gut microbiome produces vitamins, nutrients, and neurotransmitters, and helps to modulate the host immune system-and also plays a major role in the metabolism of many exogenous compounds, including drugs and chemical toxicants. However, the extent to which specific microbial species or communities modulate hazard upon exposure to chemicals remains largely opaque. Focusing on the effects of collateral dietary exposure to the widely used herbicide atrazine, we applied integrated omics and phenotypic screening to assess the role of the gut microbiome in modulating host resilience in Drosophila melanogaster. Transcriptional and metabolic responses to these compounds are sex-specific and depend strongly on the presence of the commensal microbiome. Sequencing the genomes of all abundant microbes in the fly gut revealed an enzymatic pathway responsible for atrazine detoxification unique to Acetobacter tropicalis. We find that Acetobacter tropicalis alone, in gnotobiotic animals, is sufficient to rescue increased atrazine toxicity to wild-type, conventionally reared levels. This work points toward the derivation of biotic strategies to improve host resilience to environmental chemical exposures, and illustrates the power of integrative omics to identify pathways responsible for adverse health outcomes.

Correction to: Genetic and metabolic links between the murine microbiome and memory.

An amendment to this paper has been published and can be accessed via the original article.

Genetic and metabolic links between the murine microbiome and memory.

BACKGROUND: Recent evidence has linked the gut microbiome to host behavior via the gut-brain axis [1-3]; however, the underlying mechanisms remain unexplored. Here, we determined the links between host genetics, the gut microbiome and memory using the genetically defined Collaborative Cross (CC) mouse cohort, complemented with microbiome and metabolomic analyses in conventional and germ-free (GF) mice. RESULTS: A genome-wide association analysis (GWAS) identified 715 of 76,080 single-nucleotide polymorphisms (SNPs) that were significantly associated with short-term memory using the passive avoidance model. The identified SNPs were enriched in genes known to be involved in learning and memory functions. By 16S rRNA gene sequencing of the gut microbial community in the same CC cohort, we identified specific microorganisms that were significantly correlated with longer latencies in our retention test, including a positive correlation with Lactobacillus. Inoculation of GF mice with individual species of Lactobacillus (L. reuteri F275, L. plantarum BDGP2 or L. brevis BDGP6) resulted in significantly improved memory compared to uninoculated or E. coli DH10B inoculated controls. Untargeted metabolomics analysis revealed significantly higher levels of several metabolites, including lactate, in the stools of Lactobacillus-colonized mice, when compared to GF control mice. Moreover, we demonstrate that dietary lactate treatment alone boosted memory in conventional mice. Mechanistically, we show that both inoculation with Lactobacillus or lactate treatment significantly increased the levels of the neurotransmitter, gamma-aminobutyric acid (GABA), in the hippocampus of the mice. CONCLUSION: Together, this study provides new evidence for a link between Lactobacillus and memory and our results open possible new avenues for treating memory impairment disorders using specific gut microbial inoculants and/or metabolites. Video Abstract.

Complete Genome Sequence of Agrobacterium sp. Strain 33MFTa1.1, Isolated from Thlaspi arvense Roots.

Agrobacterium sp. strain 33MFTa1.1 was isolated for functional host-microbe interaction studies from the Thlaspi arvense root-associated microbiome. The complete genome is comprised of a circular chromosome of 2,771,937 bp, a linear chromosome of 2,068,443 bp, and a plasmid of 496,948 bp, with G+C contents of 59%, 59%, and 58%, respectively.

Haplotypes spanning centromeric regions reveal persistence of large blocks of archaic DNA.

Despite critical roles in chromosome segregation and disease, the repetitive structure and vast size of centromeres and their surrounding heterochromatic regions impede studies of genomic variation. Here we report the identification of large-scale haplotypes (cenhaps) in humans that span the centromere-proximal regions of all metacentric chromosomes, including the arrays of highly repeated α-satellites on which centromeres form. Cenhaps reveal deep diversity, including entire introgressed Neanderthal centromeres and equally ancient lineages among Africans. These centromere-spanning haplotypes contain variants, including large differences in α-satellite DNA content, which may influence the fidelity and bias of chromosome transmission. The discovery of cenhaps creates new opportunities to investigate their contribution to phenotypic variation, especially in meiosis and mitosis, as well as to more incisively model the unexpectedly rich evolution of these challenging genomic regions.

Rhizobacteria Mediate the Phytotoxicity of a Range of Biorefinery-Relevant Compounds.

Advances in engineering biology have expanded the list of renewable compounds that can be produced at scale via biological routes from plant biomass. In most cases, these chemical products have not been evaluated for effects on biological systems, defined in the present study as bioactivity, that may be relevant to their manufacture. For sustainable chemical and fuel production, the industry needs to transition from fossil to renewable carbon sources, resulting in unprecedented expansion in the production and environmental distribution of chemicals used in biomanufacturing. Further, although some chemicals have been assessed for mammalian toxicity, environmental and agricultural hazards are largely unknown. We assessed 6 compounds that are representative of the emerging biofuel and bioproduct manufacturing process for their effect on model plants (Arabidopsis thaliana, Sorghum bicolor) and show that several alter plant seedling physiology at submillimolar concentrations. However, these responses change in the presence of individual bacterial species from the A. thaliana root microbiome. We identified 2 individual microbes that change the effect of chemical treatment on root architecture and a pooled microbial community with different effects relative to its constituents individually. The present study indicates that screening industrial chemicals for bioactivity on model organisms in the presence of their microbiomes is important for biologically and ecologically relevant risk analyses. Environ Toxicol Chem 2019;38:1911-1922. © 2019 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals, Inc. on behalf of SETAC.

Diverse tumour susceptibility in Collaborative Cross mice: identification of a new mouse model for human gastric tumourigenesis.

OBJECTIVE: The Collaborative Cross (CC) is a mouse population model with diverse and reproducible genetic backgrounds used to identify novel disease models and genes that contribute to human disease. Since spontaneous tumour susceptibility in CC mice remains unexplored, we assessed tumour incidence and spectrum. DESIGN: We monitored 293 mice from 18 CC strains for tumour development. Genetic association analysis and RNA sequencing were used to identify susceptibility loci and candidate genes. We analysed genomes of patients with gastric cancer to evaluate the relevance of genes identified in the CC mouse model and measured the expression levels of ISG15 by immunohistochemical staining using a gastric adenocarcinoma tissue microarray. Association of gene expression with overall survival (OS) was assessed by Kaplan-Meier analysis. RESULTS: CC mice displayed a wide range in the incidence and types of spontaneous tumours. More than 40% of CC036 mice developed gastric tumours within 1 year. Genetic association analysis identified Nfκb1 as a candidate susceptibility gene, while RNA sequencing analysis of non-tumour gastric tissues from CC036 mice showed significantly higher expression of inflammatory response genes. In human gastric cancers, the majority of human orthologues of the 166 mouse genes were preferentially altered by amplification or deletion and were significantly associated with OS. Higher expression of the CC036 inflammatory response gene signature is associated with poor OS. Finally, ISG15 protein is elevated in gastric adenocarcinomas and correlated with shortened patient OS. CONCLUSIONS: CC strains exhibit tremendous variation in tumour susceptibility, and we present CC036 as a spontaneous laboratory mouse model for studying human gastric tumourigenesis.

No difference in 4-nitroquinoline induced tumorigenesis between germ-free and colonized mice.

Variations in oral bacterial communities have been linked to oral cancer suggesting that the oral microbiome is an etiological factor that can influence oral cancer development. The 4-nitroquinoline 1-oxide (4-NQO)-induced murine oral and esophageal cancer model is frequently used to assess the effects of preventive and/or therapeutic agents. We used this model to assess the impact of the microbiome on tumorigenesis using axenic (germ-free) and conventionally housed mice. Increased toxicity was observed in germ-free mice, however, no difference in tumor incidence, multiplicity, and size was observed. Transcriptional profiling of liver tissue from germ-free and conventionally housed mice identified 254 differentially expressed genes including ten cytochrome p450 enzymes, the largest family of phase-1 drug metabolizing enzymes in the liver. Gene ontology revealed that differentially expressed genes were enriched for liver steatosis, inflammation, and oxidative stress in livers of germ-free mice. Our observations emphasize the importance of the microbiome in mediating chemical toxicity at least in part by altering host gene expression. Studies on the role of the microbiome in chemical-induced cancer using germ-free animal models should consider the potential difference in dose due to the microbiome-mediated changes in host metabolizing capacity, which might influence the ability to draw conclusions especially for tumor induction models that are dose dependent.

Short-term early exposure to thirdhand cigarette smoke increases lung cancer incidence in mice.

Exposure to thirdhand smoke (THS) is a recently described health concern that arises in many indoor environments. However, the carcinogenic potential of THS, a critical consideration in risk assessment, remains untested. Here we investigated the effects of short-term early exposure to THS on lung carcinogenesis in A/J mice. Forty weeks after THS exposure from 4 to 7 weeks of age, the mice had increased incidence of lung adenocarcinoma, tumor size and, multiplicity, compared with controls. In vitro studies using cultured human lung cancer cells showed that THS exposure induced DNA double-strand breaks and increased cell proliferation and colony formation. RNA sequencing analysis revealed that THS exposure induced endoplasmic reticulum stress and activated p53 signaling. Activation of the p53 pathway was confirmed by an increase in its targets p21 and BAX. These data indicate that early exposure to THS is associated with increased lung cancer risk.

Drosophila Histone Demethylase KDM4A Has Enzymatic and Non-enzymatic Roles in Controlling Heterochromatin Integrity.

Eukaryotic genomes are broadly divided between gene-rich euchromatin and the highly repetitive heterochromatin domain, which is enriched for proteins critical for genome stability and transcriptional silencing. This study shows that Drosophila KDM4A (dKDM4A), previously characterized as a euchromatic histone H3 K36 demethylase and transcriptional regulator, predominantly localizes to heterochromatin and regulates heterochromatin position-effect variegation (PEV), organization of repetitive DNAs, and DNA repair. We demonstrate that dKDM4A demethylase activity is dispensable for PEV. In contrast, dKDM4A enzymatic activity is required to relocate heterochromatic double-strand breaks outside the domain, as well as for organismal survival when DNA repair is compromised. Finally, DNA damage triggers dKDM4A-dependent changes in the levels of H3K56me3, suggesting that dKDM4A demethylates this heterochromatic mark to facilitate repair. We conclude that dKDM4A, in addition to its previously characterized role in euchromatin, utilizes both enzymatic and structural mechanisms to regulate heterochromatin organization and functions.

Influence of early life exposure, host genetics and diet on the mouse gut microbiome and metabolome.

Although the gut microbiome plays important roles in host physiology, health and disease1, we lack understanding of the complex interplay between host genetics and early life environment on the microbial and metabolic composition of the gut. We used the genetically diverse Collaborative Cross mouse system2 to discover that early life history impacts the microbiome composition, whereas dietary changes have only a moderate effect. By contrast, the gut metabolome was shaped mostly by diet, with specific non-dietary metabolites explained by microbial metabolism. Quantitative trait analysis identified mouse genetic trait loci (QTL) that impact the abundances of specific microbes. Human orthologues of genes in the mouse QTL are implicated in gastrointestinal cancer. Additionally, genes located in mouse QTL for Lactobacillales abundance are implicated in arthritis, rheumatic disease and diabetes. Furthermore, Lactobacillales abundance was predictive of higher host T-helper cell counts, suggesting an important link between Lactobacillales and host adaptive immunity.

Identification of genetic factors that modify motor performance and body weight using Collaborative Cross mice.

Evidence has emerged that suggests a link between motor deficits, obesity and many neurological disorders. However, the contributing genetic risk factors are poorly understood. Here we used the Collaborative Cross (CC), a large panel of newly inbred mice that captures 90% of the known variation among laboratory mice, to identify the genetic loci controlling rotarod performance and its relationship with body weight in a cohort of 365 mice across 16 CC strains. Body weight and rotarod performance varied widely across CC strains and were significantly negatively correlated. Genetic linkage analysis identified 14 loci that were associated with body weight. However, 45 loci affected rotarod performance, seven of which were also associated with body weight, suggesting a strong link at the genetic level. Lastly, we show that genes identified in this study overlap significantly with those related to neurological disorders and obesity found in human GWA studies. In conclusion, our results provide a genetic framework for studies of the connection between body weight, the central nervous system and behavior.

Comparative analysis of metazoan chromatin organization.

Genome function is dynamically regulated in part by chromatin, which consists of the histones, non-histone proteins and RNA molecules that package DNA. Studies in Caenorhabditis elegans and Drosophila melanogaster have contributed substantially to our understanding of molecular mechanisms of genome function in humans, and have revealed conservation of chromatin components and mechanisms. Nevertheless, the three organisms have markedly different genome sizes, chromosome architecture and gene organization. On human and fly chromosomes, for example, pericentric heterochromatin flanks single centromeres, whereas worm chromosomes have dispersed heterochromatin-like regions enriched in the distal chromosomal 'arms', and centromeres distributed along their lengths. To systematically investigate chromatin organization and associated gene regulation across species, we generated and analysed a large collection of genome-wide chromatin data sets from cell lines and developmental stages in worm, fly and human. Here we present over 800 new data sets from our ENCODE and modENCODE consortia, bringing the total to over 1,400. Comparison of combinatorial patterns of histone modifications, nuclear lamina-associated domains, organization of large-scale topological domains, chromatin environment at promoters and enhancers, nucleosome positioning, and DNA replication patterns reveals many conserved features of chromatin organization among the three organisms. We also find notable differences in the composition and locations of repressive chromatin. These data sets and analyses provide a rich resource for comparative and species-specific investigations of chromatin composition, organization and function.

Nucleosomes shape DNA polymorphism and divergence.

An estimated 80% of genomic DNA in eukaryotes is packaged as nucleosomes, which, together with the remaining interstitial linker regions, generate higher order chromatin structures [1]. Nucleosome sequences isolated from diverse organisms exhibit ∼10 bp periodic variations in AA, TT and GC dinucleotide frequencies. These sequence elements generate intrinsically curved DNA and help establish the histone-DNA interface. We investigated an important unanswered question concerning the interplay between chromatin organization and genome evolution: do the DNA sequence preferences inherent to the highly conserved histone core exert detectable natural selection on genomic divergence and polymorphism? To address this hypothesis, we isolated nucleosomal DNA sequences from Drosophila melanogaster embryos and examined the underlying genomic variation within and between species. We found that divergence along the D. melanogaster lineage is periodic across nucleosome regions with base changes following preferred nucleotides, providing new evidence for systematic evolutionary forces in the generation and maintenance of nucleosome-associated dinucleotide periodicities. Further, Single Nucleotide Polymorphism (SNP) frequency spectra show striking periodicities across nucleosomal regions, paralleling divergence patterns. Preferred alleles occur at higher frequencies in natural populations, consistent with a central role for natural selection. These patterns are stronger for nucleosomes in introns than in intergenic regions, suggesting selection is stronger in transcribed regions where nucleosomes undergo more displacement, remodeling and functional modification. In addition, we observe a large-scale (∼180 bp) periodic enrichment of AA/TT dinucleotides associated with nucleosome occupancy, while GC dinucleotide frequency peaks in linker regions. Divergence and polymorphism data also support a role for natural selection in the generation and maintenance of these super-nucleosomal patterns. Our results demonstrate that nucleosome-associated sequence periodicities are under selective pressure, implying that structural interactions between nucleosomes and DNA sequence shape sequence evolution, particularly in introns.

An interferon signature identified by RNA-sequencing of mammary tissues varies across the estrous cycle and is predictive of metastasis-free survival.

The concept that a breast cancer patient's menstrual stage at the time of tumor surgery influences risk of metastases remains controversial. The scarcity of comprehensive molecular studies of menstrual stage-dependent fluctuations in the breast provides little insight in this observation. To gain a deeper understanding of the biological changes in mammary tissue and blood during the menstrual cycle and to determine the influence of environmental exposures, such as low-dose ionizing radiation (LDIR), we used the mouse to characterize estrous-cycle variations in mammary gene transcripts by RNA-sequencing, peripheral white blood cell (WBC) counts and plasma cytokine levels. We identified an estrous-variable and hormone-dependent gene cluster enriched for Type-1 interferon genes. Cox regression identified a 117-gene signature of interferon-associated genes, which correlated with lower frequencies of metastasis in breast cancer patients. LDIR (10cGy) exposure had no detectable effect on mammary transcripts. However, peripheral WBC counts varied across the estrous cycle and LDIR exposure reduced lymphocyte counts and cytokine levels in tumor-susceptible mice. Our finding of variations in mammary Type-1 interferon and immune functions across the estrous cycle provides a mechanism by which timing of breast tumor surgery during the menstrual cycle may have clinical relevance to a patient's risk for distant metastases.

Genomic variation in natural populations of Drosophila melanogaster.

This report of independent genome sequences of two natural populations of Drosophila melanogaster (37 from North America and 6 from Africa) provides unique insight into forces shaping genomic polymorphism and divergence. Evidence of interactions between natural selection and genetic linkage is abundant not only in centromere- and telomere-proximal regions, but also throughout the euchromatic arms. Linkage disequilibrium, which decays within 1 kbp, exhibits a strong bias toward coupling of the more frequent alleles and provides a high-resolution map of recombination rate. The juxtaposition of population genetics statistics in small genomic windows with gene structures and chromatin states yields a rich, high-resolution annotation, including the following: (1) 5'- and 3'-UTRs are enriched for regions of reduced polymorphism relative to lineage-specific divergence; (2) exons overlap with windows of excess relative polymorphism; (3) epigenetic marks associated with active transcription initiation sites overlap with regions of reduced relative polymorphism and relatively reduced estimates of the rate of recombination; (4) the rate of adaptive nonsynonymous fixation increases with the rate of crossing over per base pair; and (5) both duplications and deletions are enriched near origins of replication and their density correlates negatively with the rate of crossing over. Available demographic models of X and autosome descent cannot account for the increased divergence on the X and loss of diversity associated with the out-of-Africa migration. Comparison of the variation among these genomes to variation among genomes from D. simulans suggests that many targets of directional selection are shared between these species.

A large-scale screen for mutagen-sensitive loci in Drosophila.

In a screen for new DNA repair mutants, we tested 6275 Drosophila strains bearing homozygous mutagenized autosomes (obtained from C. Zuker) for hypersensitivity to methyl methanesulfonate (MMS) and nitrogen mustard (HN2). Testing of 2585 second-chromosome lines resulted in the recovery of 18 mutants, 8 of which were alleles of known genes. The remaining 10 second-chromosome mutants were solely sensitive to MMS and define 8 new mutagen-sensitive genes (mus212-mus219). Testing of 3690 third chromosomes led to the identification of 60 third-chromosome mutants, 44 of which were alleles of known genes. The remaining 16 mutants define 14 new mutagen-sensitive genes (mus314-mus327). We have initiated efforts to identify these genes at the molecular level and report here the first two identified. The HN2-sensitive mus322 mutant defines the Drosophila ortholog of the yeast snm1 gene, and the MMS- and HN2-sensitive mus301 mutant defines the Drosophila ortholog of the human HEL308 gene. We have also identified a second-chromosome mutant, mus215(ZIII-2059), that uniformly reduces the frequency of meiotic recombination to <3% of that observed in wild type and thus defines a function required for both DNA repair and meiotic recombination. At least one allele of each new gene identified in this study is available at the Bloomington Stock Center.

Obtaining maximal concatenated phylogenetic data sets from large sequence databases.

To improve the accuracy of tree reconstruction, phylogeneticists are extracting increasingly large multigene data sets from sequence databases. Determining whether a database contains at least k genes sampled from at least m species is an NP-complete problem. However, the skewed distribution of sequences in these databases permits all such data sets to be obtained in reasonable computing times even for large numbers of sequences. We developed an exact algorithm for obtaining the largest multigene data sets from a collection of sequences. The algorithm was then tested on a set of 100,000 protein sequences of green plants and used to identify the largest multigene ortholog data sets having at least 3 genes and 6 species. The distribution of sizes of these data sets forms a hollow curve, and the largest are surprisingly small, ranging from 62 genes by 6 species, to 3 genes by 65 species, with more symmetrical data sets of around 15 taxa by 15 genes. These upper bounds to sequence concatenation have important implications for building the tree of life from large sequence databases.