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.

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.

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.

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.

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.

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.

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.

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.

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.