Ligand binding and protein dynamics: a fluorescence depolarization study of aspartate transcarbamylase from Escherichia coli

The chromatin remodeling enzyme Chd4 regulates genome architecture in the mouse brain

The development and function of the brain require tight control of gene expression. Genome architecture is thought to play a critical regulatory role in gene expression, but the mechanisms governing genome architecture in the brain in vivo remain poorly understood. Here, we report that conditional knockout of the chromatin remodeling enzyme Chd4 in granule neurons of the mouse cerebellum increases accessibility of gene regulatory sites genome-wide in vivo. Conditional knockout of Chd4 promotes recruitment of the architectural protein complex cohesin preferentially to gene enhancers in granule neurons in vivo. Importantly, in vivo profiling of genome architecture reveals that conditional knockout of Chd4 strengthens interactions among developmentally repressed contact domains as well as genomic loops in a manner that tightly correlates with increased accessibility, enhancer activity, and cohesin occupancy at these sites. Collectively, our findings define a role for chromatin remodeling in the control of genome architecture organization in the mammalian brain.

Sequence-based genetic mapping of Cynodon dactylon Pers. reveals new insights into genome evolution in Poaceae

Bermudagrass (Cynodon dactylon Pers.) is an important warm-season perennial used extensively for turf, forage, soil conservation and remediation worldwide. However, limited genomic information has hindered the application of molecular tools towards understanding genome evolution and in breeding new cultivars. We genotype a first-generation selfed population derived from the tetraploid (4x = 36) ‘A12359’ using genotyping-by-sequencing. A high-density genetic map of 18 linkage groups (LGs) is constructed with 3,544 markers. Comparative genomic analyses reveal that each of nine homeologous LG pairs of C. dactylon corresponds to one of the first nine chromosomes of Oropetium thomaeum. Two nested paleo-ancestor chromosome fusions (ρ6-ρ9-ρ6, ρ2-ρ10-ρ2) may have resulted in a 12-to-10 chromosome reduction. A segmental dissemination of the paleo-chromosome ρ12 (ρ1-ρ12-ρ1, ρ6-ρ12-ρ6) leads to the 10-to-9 chromosome reduction in C. dactylon genome. The genetic map will assist in an ongoing whole genome sequence assembly and facilitate marker-assisted selection (MAS) in developing new cultivars.

Tilin Fang et al. study the genome of Bermudagrass (Cynodon dactylon Pers.). They use genotyping-by-sequencing and provide a genetic map with a 10-fold increase in genetic marker density. Comparative genomics analyses reveal chromosome rearrangements. Their work will contribute to whole genome assembly efforts which will be beneficial for developing new cultivars.

Select amino acids in DGCR8 are essential for the UGU-pri-miRNA interaction and processing

Microprocessor, composed of DROSHA and DGCR8, processes primary microRNAs (pri-miRNAs) in miRNA biogenesis. Its cleavage efficiency and accuracy are enhanced because DGCR8 interacts with the apical UGU motif of pri-miRNAs. However, the mechanism and influence of DGCR8–UGU interaction on cellular miRNA expression are still elusive. In this study, we demonstrated that Rhed (i.e., the RNA-binding heme domain, amino acids 285–478) of DGCR8 interacts with UGU. In addition, we identified three amino acids 461–463 in Rhed, which are critical for the UGU interaction and essential for Microprocessor to accurately and efficiently process UGU-pri-miRNAs in vitro and UGU-miRNA expression in human cells. Furthermore, we found that within the DGCR8 dimer, the amino acids 461–463 from one monomer are capable of discriminating between UGU- and noUGU-pri-miRNAs. Our findings improve the current understanding of the substrate-recognizing mechanism of DGCR8 and implicate the roles of this recognition in differentiating miRNA expression in human cells.

Thi Lieu Dang et al. study the mechanisms for the interaction between DGCR8 and the apical UGU motif of pri-miRNAs. They demonstrate that three amino acids in the Rhed domain of DGCR8 are critical to recognize and interact with UGU and to process pri-miRNAs. They further show amino acids 461–463 in one of the DGCR8 dimer are necessary to distinguish UGU- and noUGU-primiRNAs.

Laboratory culture of the California Sea Firefly Vargula tsujii (Ostracoda: Cypridinidae): Developing a model system for the evolution of marine bioluminescence

Bioluminescence, or the production of light by living organisms via chemical reaction, is widespread across Metazoa. Laboratory culture of bioluminescent organisms from diverse taxonomic groups is important for determining the biosynthetic pathways of bioluminescent substrates, which may lead to new tools for biotechnology and biomedicine. Some bioluminescent groups may be cultured, including some cnidarians, ctenophores, and brittle stars, but those use luminescent substrates (luciferins) obtained from their diets, and therefore are not informative for determination of the biosynthetic pathways of the luciferins. Other groups, including terrestrial fireflies, do synthesize their own luciferin, but culturing them is difficult and the biosynthetic pathway for firefly luciferin remains unclear. An additional independent origin of endogenous bioluminescence is found within ostracods from the family Cypridinidae, which use their luminescence for defense and, in Caribbean species, for courtship displays. Here, we report the first complete life cycle of a luminous ostracod (Vargula tsujii Kornicker & Baker, 1977, the California Sea Firefly) in the laboratory. We also describe the late-stage embryogenesis of Vargula tsujii and discuss the size classes of instar development. We find embryogenesis in V. tsujii ranges from 25–38 days, and this species appears to have five instar stages, consistent with ontogeny in other cypridinid lineages. We estimate a complete life cycle at 3–4 months. We also present the first complete mitochondrial genome for Vargula tsujii. Bringing a luminous ostracod into laboratory culture sets the stage for many potential avenues of study, including learning the biosynthetic pathway of cypridinid luciferin and genomic manipulation of an autogenic bioluminescent system.

Identification of distinct loci for de novo DNA methylation by DNMT3A and DNMT3B during mammalian development

De novo establishment of DNA methylation is accomplished by DNMT3A and DNMT3B. Here, we analyze de novo DNA methylation in mouse embryonic fibroblasts (2i-MEFs) derived from DNA-hypomethylated 2i/L ES cells with genetic ablation of Dnmt3a or Dnmt3b. We identify 355 and 333 uniquely unmethylated genes in Dnmt3a and Dnmt3b knockout (KO) 2i-MEFs, respectively. We find that Dnmt3a is exclusively required for de novo methylation at both TSS regions and gene bodies of Polycomb group (PcG) target developmental genes, while Dnmt3b has a dominant role on the X chromosome. Consistent with this, tissue-specific DNA methylation at PcG target genes is substantially reduced in Dnmt3a KO embryos. Finally, we find that human patients with DNMT3 mutations exhibit reduced DNA methylation at regions that are hypomethylated in Dnmt3 KO 2i-MEFs. In conclusion, here we report a set of unique de novo DNA methylation target sites for both DNMT3 enzymes during mammalian development that overlap with hypomethylated sites in human patients.

Measuring and interpreting transposable element expression

Transposable elements (TEs) are insertional mutagens that contribute greatly to the plasticity of eukaryotic genomes, influencing the evolution and adaptation of species as well as physiology or disease in individuals. Measuring TE expression helps to understand not only when and where TE mobilization can occur but also how this process alters gene expression, chromatin accessibility or cellular signalling pathways. Although genome-wide gene expression assays such as RNA sequencing include transposon-derived transcripts, most computational analytical tools discard or misinterpret TE-derived reads. Emerging approaches are improving the identification of expressed TE loci and helping to discriminate TE transcripts that permit TE mobilization from chimeric gene-TE transcripts or pervasive transcription. Here we review the main challenges associated with the detection of TE expression, including mappability, insertional and internal sequence polymorphisms, and the diversity of the TE transcriptional landscape, as well as the different experimental and computational strategies to solve them.

Functional shell matrix proteins tentatively identified by asymmetric snail shell morphology

Molluscan shell matrix proteins (SMPs) are essential in biomineralization. Here, we identify potentially important SMPs by exploiting the asymmetric shell growth in snail, Lymnaea stagnalis. Asymmetric shells require bilaterally asymmetric expression of SMP genes. We examined expression levels of 35,951 transcripts expressed in the left and right sides of mantle tissue of the pond snail, Lymnaea stagnalis. This transcriptome dataset was used to identify 207 SMPs by LC-MS/MS. 32 of the 207 SMP genes show asymmetric expression patterns, which were further verified for 4 of the 32 SMPs using quantitative PCR analysis. Among asymmetrically expressed SMPs in dextral snails, those that are more highly expressed on the left side than the right side are 3 times more abundant than those that are more highly expressed on the right than the left, suggesting potentially inhibitory roles of SMPs in shell formation. The 32 SMPs thus identified have distinctive features, such as conserved domains and low complexity regions, which may be essential in biomineralization.

Pseudomonas aeruginosa core metabolism exerts a widespread growth-independent control on virulence

To assess the role of core metabolism genes in bacterial virulence - independently of their effect on growth - we correlated the genome, the transcriptome and the pathogenicity in flies and mice of 30 fully sequenced Pseudomonas strains. Gene presence correlates robustly with pathogenicity differences among all Pseudomonas species, but not among the P. aeruginosa strains. However, gene expression differences are evident between highly and lowly pathogenic P. aeruginosa strains in multiple virulence factors and a few metabolism genes. Moreover, 16.5%, a noticeable fraction of the core metabolism genes of P. aeruginosa strain PA14 (compared to 8.5% of the non-metabolic genes tested), appear necessary for full virulence when mutated. Most of these virulence-defective core metabolism mutants are compromised in at least one key virulence mechanism independently of auxotrophy. A pathway level analysis of PA14 core metabolism, uncovers beta-oxidation and the biosynthesis of amino-acids, succinate, citramalate, and chorismate to be important for full virulence. Strikingly, the relative expression among P. aeruginosa strains of genes belonging in these metabolic pathways is indicative of their pathogenicity. Thus, P. aeruginosa strain-to-strain virulence variation, remains largely obscure at the genome level, but can be dissected at the pathway level via functional transcriptomics of core metabolism.

EXOSC9 depletion attenuates P-body formation, stress resistance, and tumorigenicity of cancer cells

Cancer cells adapt to various stress conditions by optimizing gene expression profiles via transcriptional and translational regulation. However, whether and how EXOSC9, a component of the RNA exosome complex, regulates adaptation to stress conditions and tumorigenicity in cancer cells remain unclear. Here, we examined the effects of EXOSC9 depletion on cancer cell growth under various stress conditions. EXOSC9 depletion attenuated growth and survival under various stress conditions in cancer cells. Interestingly, this also decreased the number of P-bodies, which are messenger ribonucleoprotein particles (mRNPs) required for stress adaptation. Meanwhile, EXOSC2/EXOSC4 depletion also attenuated P-body formation and stress resistance with decreased EXOSC9 protein. EXOSC9-mediated stress resistance and P-body formation were found to depend on the intact RNA-binding motif of this protein. Further, RNA-seq analyses identified 343 EXOSC9-target genes, among which, APOBEC3G contributed to defects in stress resistance and P-body formation in MDA-MB-231 cells. Finally, EXOSC9 also promoted xenografted tumor growth of MDA-MB-231 cells in an intact RNA-binding motif-dependent manner. Database analyses further showed that higher EXOSC9 activity, estimated based on the expression of 343 target genes, was correlated with poorer prognosis in some cancer patients. Thus, drugs targeting activity of the RNA exosome complex or EXOSC9 might be useful for cancer treatment.

Three-dimensional genome restructuring across timescales of activity-induced neuronal gene expression

Neuronal activation induces rapid transcription of immediate early genes (IEGs) and longer-term chromatin remodeling around secondary response genes (SRGs). Here, we use high-resolution chromosome-conformation-capture carbon-copy sequencing (5C-seq) to elucidate the extent to which long-range chromatin loops are altered during short- and long-term changes in neural activity. We find that more than 10% of loops surrounding select IEGs, SRGs, and synaptic genes are induced de novo during cortical neuron activation. IEGs Fos and Arc connect to activity-dependent enhancers via singular short-range loops that form within 20 min after stimulation, prior to peak messenger RNA levels. By contrast, the SRG Bdnf engages in both pre-existing and activity-inducible loops that form within 1-6 h. We also show that common single-nucleotide variants that are associated with autism and schizophrenia are colocalized with distinct classes of activity-dependent, looped enhancers. Our data link architectural complexity to transcriptional kinetics and reveal the rapid timescale by which higher-order chromatin architecture reconfigures during neuronal stimulation.

Targeted reprogramming of H3K27me3 resets epigenetic memory in plant paternal chromatin

Epigenetic marks are reprogrammed in the gametes to reset genomic potential in the next generation. In mammals, paternal chromatin is extensively reprogrammed through the global erasure of DNA methylation and the exchange of histones with protamines1,2. Precisely how the paternal epigenome is reprogrammed in flowering plants has remained unclear since DNA is not demethylated and histones are retained in sperm3,4. Here, we describe a multi-layered mechanism by which H3K27me3 is globally lost from histone-based sperm chromatin in Arabidopsis. This mechanism involves the silencing of H3K27me3 writers, activity of H3K27me3 erasers and deposition of a sperm-specific histone, H3.10 (ref. 5), which we show is immune to lysine 27 methylation. The loss of H3K27me3 facilitates the transcription of genes essential for spermatogenesis and pre-configures sperm with a chromatin state that forecasts gene expression in the next generation. Thus, plants have evolved a specific mechanism to simultaneously differentiate male gametes and reprogram the paternal epigenome.

NFI transcription factors provide chromatin access to maintain stem cell identity while preventing unintended lineage fate choices

Tissue homeostasis and regeneration rely on resident stem cells (SCs), whose behaviour is regulated through niche-dependent crosstalk. The mechanisms underlying SC identity are still unfolding. Here, using spatiotemporal gene ablation in murine hair follicles, we uncover a critical role for the transcription factors (TFs) nuclear factor IB (NFIB) and IX (NFIX) in maintaining SC identity. Without NFI TFs, SCs lose their hair-regenerating capability, and produce skin bearing striking resemblance to irreversible human alopecia, which also displays reduced NFIs. Through single-cell transcriptomics, ATAC-Seq and ChIP-Seq profiling, we expose a key role for NFIB and NFIX in governing super-enhancer maintenance of the key hair follicle SC-specific TF genes. When NFIB and NFIX are genetically removed, the stemness epigenetic landscape is lost. Super-enhancers driving SC identity are decommissioned, while unwanted lineages are de-repressed ectopically. Together, our findings expose NFIB and NFIX as crucial rheostats of tissue homeostasis, functioning to safeguard the SC epigenome from a breach in lineage confinement that otherwise triggers irreversible tissue degeneration.

A calcineurin-Hoxb13 axis regulates growth mode of mammalian cardiomyocytes

A major factor in the progression to heart failure in humans is the inability of the adult heart to repair itself after injury. We recently demonstrated that the early postnatal mammalian heart is capable of regeneration following injury through proliferation of preexisting cardiomyocytes1,2 and that Meis1, a three amino acid loop extension (TALE) family homeodomain transcription factor, translocates to cardiomyocyte nuclei shortly after birth and mediates postnatal cell cycle arrest3. Here we report that Hoxb13 acts as a cofactor of Meis1 in postnatal cardiomyocytes. Cardiomyocyte-specific deletion of Hoxb13 can extend the postnatal window of cardiomyocyte proliferation and reactivate the cardiomyocyte cell cycle in the adult heart. Moreover, adult Meis1-Hoxb13 double-knockout hearts display widespread cardiomyocyte mitosis, sarcomere disassembly and improved left ventricular systolic function following myocardial infarction, as demonstrated by echocardiography and magnetic resonance imaging. Chromatin immunoprecipitation with sequencing demonstrates that Meis1 and Hoxb13 act cooperatively to regulate cardiomyocyte maturation and cell cycle. Finally, we show that the calcium-activated protein phosphatase calcineurin dephosphorylates Hoxb13 at serine-204, resulting in its nuclear localization and cell cycle arrest. These results demonstrate that Meis1 and Hoxb13 act cooperatively to regulate cardiomyocyte maturation and proliferation and provide mechanistic insights into the link between hyperplastic and hypertrophic growth of cardiomyocytes.

Genomic prediction applied to multiple traits and environments in second season maize hybrids

Genomic selection has become a reality in plant breeding programs with the reduction in genotyping costs. Especially in maize breeding programs, it emerges as a promising tool for predicting hybrid performance. The dynamics of a commercial breeding program involve the evaluation of several traits simultaneously in a large set of target environments. Therefore, multi-trait multi-environment (MTME) genomic prediction models can leverage these datasets by exploring the correlation between traits and Genotype-by-Environment (G×E) interaction. Herein, we assess predictive abilities of univariate and multivariate genomic prediction models in a maize breeding program. To this end, we used data from 415 maize hybrids evaluated in 4 years of second season field trials for the traits grain yield, number of ears, and grain moisture. Genotypes of these hybrids were inferred in silico based on their parental inbred lines using single nucleotide polymorphisms (SNPs) markers obtained via genotyping-by-sequencing (GBS). Because genotypic information was available for only 257 hybrids, we used the genomic and pedigree relationship matrices to obtain the H matrix for all 415 hybrids. Our results demonstrated that in the single-environment context the use of multi-trait models was always superior in comparison to their univariate counterparts. Besides that, although MTME models were not particularly successful in predicting hybrid performance in untested years, they improved the ability to predict the performance of hybrids that had not been evaluated in any environment. However, the computational requirements of this kind of model could represent a limitation to its practical implementation and further investigation is necessary.

A phylogenetically novel cyanobacterium most closely related to Gloeobacter

Clues to the evolutionary steps producing innovations in oxygenic photosynthesis may be preserved in the genomes of organisms phylogenetically placed between non-photosynthetic Vampirovibrionia (formerly Melainabacteria) and the thylakoid-containing Cyanobacteria. However, only two species with published genomes are known to occupy this phylogenetic space, both within the genus Gloeobacter. Here, we describe nearly complete, metagenome-assembled genomes (MAGs) of an uncultured organism phylogenetically placed near Gloeobacter, for which we propose the name Candidatus Aurora vandensis {Au’ro.ra. L. fem. n. aurora, the goddess of the dawn in Roman mythology; van.de’nsis. N.L. fem. adj. vandensis of Lake Vanda, Antarctica}. The MAG of A. vandensis contains homologs of most genes necessary for oxygenic photosynthesis including key reaction center proteins. Many accessory subunits associated with the photosystems in other species either are missing from the MAG or are poorly conserved. The MAG also lacks homologs of genes associated with the pigments phycocyanoerethrin, phycoeretherin and several structural parts of the phycobilisome. Additional characterization of this organism is expected to inform models of the evolution of oxygenic photosynthesis.

Metagenome-assembled genomes reveal unique metabolic adaptations of a basal marine Thaumarchaeota lineage

Thaumarchaeota constitute an abundant and ubiquitous phylum of Archaea that play critical roles in the global nitrogen and carbon cycles. Most well-characterized members of the phylum are chemolithoautotrophic ammonia-oxidizing archaea (AOA), which comprise up to 5 and 20% of the total single-celled life in soil and marine systems, respectively. Using two high-quality metagenome-assembled genomes (MAGs), here we describe a divergent marine thaumarchaeal clade that is devoid of the ammonia-oxidation machinery and the AOA-specific carbon-fixation pathway. Phylogenomic analyses placed these genomes within the uncultivated and largely understudied marine pSL12-like thaumarchaeal clade. The predominant mode of nutrient acquisition appears to be aerobic heterotrophy, evidenced by the presence of respiratory complexes and various organic carbon degradation pathways. Both genomes encoded several pyrroloquinoline quinone (PQQ)-dependent alcohol dehydrogenases, as well as a form III RuBisCO. Metabolic reconstructions suggest anaplerotic CO₂ assimilation mediated by RuBisCO, which may be linked to the central carbon metabolism. We conclude that these genomes represent a hitherto unrecognized evolutionary link between predominantly anaerobic basal thaumarchaeal lineages and mesophilic marine AOA, with important implications for diversification within the phylum Thaumarchaeota.

BMP4 resets mouse epiblast stem cells to naive pluripotency through ZBTB7A/B-mediated chromatin remodelling

BMP4 regulates a plethora of developmental processes, including the dorsal-ventral axis and neural patterning. Here, we report that BMP4 reconfigures the nuclear architecture during the primed-to-naive transition (PNT). We first established a BMP4-driven PNT and show that BMP4 orchestrates the chromatin accessibility dynamics during PNT. Among the loci opened early by BMP4, we identified Zbtb7a and Zbtb7b (Zbtb7a/b) as targets that drive PNT. ZBTB7A/B in turn facilitate the opening of naive pluripotent chromatin loci and the activation of nearby genes. Mechanistically, ZBTB7A not only binds to chromatin loci near to the genes that are activated, but also strategically occupies those that are silenced, consistent with a role of BMP4 in both activating and suppressing gene expression during PNT at the chromatin level. Our results reveal a previously unknown function of BMP4 in regulating nuclear architecture and link its targets ZBTB7A/B to chromatin remodelling and pluripotent fate control.

Genome diversification in globally distributed novel marine Proteobacteria is linked to environmental adaptation

Proteobacteria constitute one of the most diverse and abundant groups of microbes on Earth. In productive marine environments like deep-sea hydrothermal systems, Proteobacteria are implicated in autotrophy coupled to sulfur, methane, and hydrogen oxidation, sulfate reduction, and denitrification. Beyond chemoautotrophy, little is known about the ecological significance of poorly studied Proteobacteria lineages that are globally distributed and active in hydrothermal systems. Here we apply multi-omics to characterize 51 metagenome-assembled genomes from three hydrothermal vent plumes in the Pacific and Atlantic Oceans that are affiliated with nine Proteobacteria lineages. Metabolic analyses revealed these organisms to contain a diverse functional repertoire including chemolithotrophic ability to utilize sulfur and C₁ compounds, and chemoorganotrophic ability to utilize environment-derived fatty acids, aromatics, carbohydrates, and peptides. Comparative genomics with marine and terrestrial microbiomes suggests that lineage-associated functional traits could explain niche specificity. Our results shed light on the ecological functions and metabolic strategies of novel Proteobacteria in hydrothermal systems and beyond, and highlight the relationship between genome diversification and environmental adaptation.

Bifidobacteria shape host neural circuits during postnatal development by promoting synapse formation and microglial function

We hypothesized that early-life gut microbiota support the functional organization of neural circuitry in the brain via regulation of synaptic gene expression and modulation of microglial functionality. Germ-free mice were colonized as neonates with either a simplified human infant microbiota consortium consisting of four Bifidobacterium species, or with a complex, conventional murine microbiota. We examined the cerebellum, cortex, and hippocampus of both groups of colonized mice in addition to germ-free control mice. At postnatal day 4 (P4), conventionalized mice and Bifidobacterium-colonized mice exhibited decreased expression of synapse-promoting genes and increased markers indicative of reactive microglia in the cerebellum, cortex and hippocampus relative to germ-free mice. By P20, both conventional and Bifidobacterium-treated mice exhibited normal synaptic density and neuronal activity as measured by density of VGLUT2+ puncta and Purkinje cell firing rate respectively, in contrast to the increased synaptic density and decreased firing rate observed in germ-free mice. The conclusions from this study further reveal how bifidobacteria participate in establishing functional neural circuits. Collectively, these data indicate that neonatal microbial colonization of the gut elicits concomitant effects on the host CNS, which promote the homeostatic developmental balance of neural connections during the postnatal time period.

Genetic basis for virulence differences of various Cryptosporidium parvum carcinogenic isolates

Cryptosporidium parvum is known to cause life-threatening diarrhea in immunocompromised hosts and was also reported to be capable of inducing digestive adenocarcinoma in a rodent model. Interestingly, three carcinogenic isolates of C. parvum, called DID, TUM1 and CHR, obtained from fecal samples of naturally infected animals or humans, showed higher virulence than the commercially available C. parvum IOWA isolate in our animal model in terms of clinical manifestations, mortality rate and time of onset of neoplastic lesions. In order to discover the potential genetic basis of the differential virulence observed between C. parvum isolates and to contribute to the understanding of Cryptosporidium virulence, entire genomes of the isolates DID, TUM1 and CHR were sequenced then compared to the C. parvum IOWA reference genome. 125 common SNVs corresponding to 90 CDSs were found in the C. parvum genome that could explain this differential virulence. In particular variants in several membrane and secreted proteins were identified. Besides the genes already known to be involved in parasite virulence, this study identified potential new virulence factors whose functional characterization can be achieved through CRISPR/Cas9 technology applied to this parasite.

Targeted analysis of polymorphic loci from low-coverage shotgun sequence data allows accurate genotyping of HLA genes in historical human populations

The highly polymorphic human leukocyte antigen (HLA) plays a crucial role in adaptive immunity and is associated with various complex diseases. Accurate analysis of HLA genes using ancient DNA (aDNA) data is crucial for understanding their role in human adaptation to pathogens. Here, we describe the TARGT pipeline for targeted analysis of polymorphic loci from low-coverage shotgun sequence data. The pipeline was successfully applied to medieval aDNA samples and validated using both simulated aDNA and modern empirical sequence data from the 1000 Genomes Project. Thus the TARGT pipeline enables accurate analysis of HLA polymorphisms in historical (and modern) human populations.

Mismatched and wobble base pairs govern primary microRNA processing by human Microprocessor

MicroRNAs (miRNAs) are small RNAs that regulate gene expression. miRNAs are produced from primary miRNAs (pri-miRNAs), which are cleaved by Microprocessor. Microprocessor, therefore, plays a crucial role in determining the efficiency and precision of miRNA production, and thus the function of the final miRNA product. Here, we conducted high-throughput enzymatic assays to investigate the catalytic mechanism of Microprocessor cleaving randomized pri-miRNAs. We identified multiple mismatches and wobble base pairs in the upper stem of pri-miRNAs, which influence the efficiency and accuracy of their processing. The existence of these RNA elements helps to explain the alternative cleavage of Microprocessor for some human pri-miRNAs. We also demonstrated that miRNA biogenesis can be altered via modification of the RNA elements by RNA-editing events or single nucleotide polymorphisms (SNPs). These findings improve our understanding of pri-miRNA processing mechanisms and provide a foundation for interpreting differential miRNA expression due to RNA modifications and SNPs.

MicroRNA genes are transcribed to long primary transcripts called primary microRNAs, which are cleaved by Microprocessor. Here the authors employ high-throughput sequencing and Microprocessor assay to show that mismatches and wobble base pairs in primary microRNAs affect the accuracy and efficiency of Microprocessor processing.

Super-resolution imaging reveals the evolution of higher-order chromatin folding in early carcinogenesis

Genomic DNA is folded into a higher-order structure that regulates transcription and maintains genomic stability. Although progress has been made on understanding biochemical characteristics of epigenetic modifications in cancer, the in-situ higher-order folding of chromatin structure during malignant transformation remains largely unknown. Here, using optimized stochastic optical reconstruction microscopy (STORM) for pathological tissue (PathSTORM), we uncover a gradual decompaction and fragmentation of higher-order chromatin folding throughout all stages of carcinogenesis in multiple tumor types, and prior to tumor formation. Our integrated imaging, genomic, and transcriptomic analyses reveal functional consequences in enhanced transcription activities and impaired genomic stability. We also demonstrate the potential of imaging higher-order chromatin disruption to detect high-risk precursors that cannot be distinguished by conventional pathology. Taken together, our findings reveal gradual decompaction and fragmentation of higher-order chromatin structure as an enabling characteristic in early carcinogenesis to facilitate malignant transformation, which may improve cancer diagnosis, risk stratification, and prevention.

The stepwise assembly of the neonatal virome is modulated by breastfeeding

The gut of healthy human neonates is usually devoid of viruses at birth, but quickly becomes colonized, in some cases leading to gastrointestinal disorders^1–4. Here we report that viral community assembly in neonates takes place in distinct steps. Fluorescent staining of virus-like particles purified from infant meconium/early stool samples show few or no particles, but by one month of life particle numbers achieve 10⁹ per gram, and these numbers appear to persist through life^5–7. We investigated the origin of these viral populations using shotgun metagenomic sequencing of viral-enriched preparations and whole microbial communities, and followed up with targeted microbiological analyses. Results indicate that, early after birth, pioneer bacteria colonize the infant gut, and by one month prophage induced from these bacteria provide the predominant population of virus-like particles. By four months of life, identifiable viruses that replicate in human cells become more prominent. Multiple human viruses were more abundant in stool samples from babies exclusively fed formula versus those fed partially or fully on breast milk, paralleling reports that breast milk can be protective against viral infections^8–10. Phage populations also differed associated with breastfeeding. Evidently colonization of the infant gut is stepwise, first mainly by temperate bacteriophages induced from pioneer bacteria, and later by viruses that replicate in human cells, with the second phase modulated by breastfeeding.

Merkel cell polyomavirus activates LSD1-mediated blockade of non-canonical BAF to regulate transformation and tumorigenesis

Merkel cell carcinoma (MCC), a neuroendocrine cancer of the skin, is caused by integration of Merkel cell polyomavirus (MCV) and persistent expression of Large T antigen (LT) and Small T antigen (ST). We report that ST in complex with MYCL and the EP400 complex activates expression of LSD1 (KDM1A), RCOR2, and INSM1 to repress gene expression by the lineage transcription factor ATOH1. LSD1 inhibition reduces growth of MCC in vitro and in vivo. Through a forward-genetics CRISPR-Cas9 screen, we identified an antagonistic relationship between LSD1 and the non-canonical BAF (ncBAF) chromatin remodeling complex. Changes in gene expression and chromatin accessibility caused by LSD1 inhibition could be partially rescued by BRD9 inhibition, revealing that LSD1 and ncBAF antagonistically regulate an overlapping set of genes. Our work provides mechanistic insight into the dependence of MCC on LSD1 and a tumor suppressor role for ncBAF in cancer.

RXRs control serous macrophage neonatal expansion and identity and contribute to ovarian cancer progression

Tissue-resident macrophages (TRMs) populate all tissues and play key roles in homeostasis, immunity and repair. TRMs express a molecular program that is mostly shaped by tissue cues. However, TRM identity and the mechanisms that maintain TRMs in tissues remain poorly understood. We recently found that serous-cavity TRMs (LPMs) are highly enriched in RXR transcripts and RXR-response elements. Here, we show that RXRs control mouse serous-macrophage identity by regulating chromatin accessibility and the transcriptional regulation of canonical macrophage genes. RXR deficiency impairs neonatal expansion of the LPM pool and reduces the survival of adult LPMs through excess lipid accumulation. We also find that peritoneal LPMs infiltrate early ovarian tumours and that RXR deletion diminishes LPM accumulation in tumours and strongly reduces ovarian tumour progression in mice. Our study reveals that RXR signalling controls the maintenance of the serous macrophage pool and that targeting peritoneal LPMs may improve ovarian cancer outcomes.

Macrophages can differentiate to perform homeostatic tissue-specific functions. Here the authors show that RXR signalling is critical for large peritoneal macrophage (LPM) expansion during neonatal life and LPM lipid metabolism and survival during adult homeostasis, and that ovarian cancer growth relies on RXR-dependent LPMs. 

Rapid morphological divergence following a human-mediated introduction: the role of drift and directional selection

Theory predicts that when populations are established by few individuals, random founder effects can facilitate rapid phenotypic divergence even in the absence of selective processes. However, empirical evidence from historically documented colonisations suggest that, in most cases, drift alone is not sufficient to explain the rate of morphological divergence. Here, using the human-mediated introduction of the silvereye (Zosterops lateralis) to French Polynesia, which represents a potentially extreme example of population founding, we reassess the potential for morphological shifts to arise via drift alone. Despite only 80 years of separation from their New Zealand ancestors, French Polynesian silvereyes displayed significant changes in body and bill size and shape, most of which could be accounted for by drift, without the need to invoke selection. However, signatures of selection at genes previously identified as candidates for bill size and body shape differences in a range of bird species, also suggests a role for selective processes in driving morphological shifts within this population. Twenty-four SNPs in our RAD-Seq dataset were also found to be strongly associated with phenotypic variation. Hence, even under population founding extremes, when it is difficult to reject drift as the sole mechanism based on rate tests of phenotypic shifts, the additional role of divergent natural selection in novel environments can be revealed at the level of the genome.

Gut microbiome diversity detected by high-coverage 16S and shotgun sequencing of paired stool and colon sample

The gut microbiome has a fundamental role in human health and disease. However, studying the complex structure and function of the gut microbiome using next generation sequencing is challenging and prone to reproducibility problems. Here, we obtained cross-sectional colon biopsies and faecal samples from nine participants in our COLSCREEN study and sequenced them in high coverage using Illumina pair-end shotgun (for faecal samples) and IonTorrent 16S (for paired feces and colon biopsies) technologies. The metagenomes consisted of between 47 and 92 million reads per sample and the targeted sequencing covered more than 300 k reads per sample across seven hypervariable regions of the 16S gene. Our data is freely available and coupled with code for the presented metagenomic analysis using up-to-date bioinformatics algorithms. These results will add up to the informed insights into designing comprehensive microbiome analysis and also provide data for further testing for unambiguous gut microbiome analysis.

Cell lineage-specific transcriptome analysis for interpreting cell fate specification of proembryos

In Arabidopsis, a zygote undergoes asymmetrical cell division that establishes the first two distinct cell types of early proembryos, apical and basal cells. However, the genome-wide transcriptional activities that guide divergence of apical and basal cell development remain unknown. Here, we present a comprehensive transcriptome analysis of apical and basal cell lineages, uncovering distinct molecular pathways during cell lineage specification. Selective deletion of inherited transcripts and specific de novo transcription contribute to the establishment of cell lineage-specific pathways for cell fate specification. Embryo-related pathways have been specifically activated in apical cell lineage since 1-cell embryo stage, but quick transcriptome remodeling toward suspensor-specific pathways are found in basal cell lineage. Furthermore, long noncoding RNAs and alternative splicing isoforms may be involved in cell lineage specification. This work also provides a valuable lineage-specific transcriptome resource to elucidate the molecular pathways for divergence of apical and basal cell lineages at genome-wide scale.

Asymmetric division of the Arabidopsis zygote produces apical and basal cells that mainly develop into embryo and suspensor, respectively. Here, Zhou et al. show that de novo transcription and selective RNA turnover establish distinct apical and basal transcriptomes as early as the 1-cell stage of embryo development.

Enhanced nucleotide chemistry and toehold nanotechnology reveals lncRNA spreading on chromatin

Understanding the targeting and spreading patterns of long non-coding RNAs (lncRNAs) on chromatin requires a technique that can detect both high-intensity binding sites and reveal genome-wide changes in spreading patterns with high precision and confidence. Here we determine lncRNA localization using biotinylated locked nucleic acid (LNA)-containing oligonucleotides with toehold architecture capable of hybridizing to target RNA through strand-exchange reaction. During hybridization, a protecting strand competitively displaces contaminating species, leading to highly specific RNA capture of individual RNAs. Analysis of Drosophila roX2 lncRNA using this approach revealed that heat shock, unlike the unfolded protein response, leads to reduced spreading of roX2 on the X chromosome, but surprisingly also to relocalization to sites on autosomes. Our results demonstrate that this improved hybridization capture approach can reveal previously uncharacterized changes in the targeting and spreading of lncRNAs on chromatin.

Uncovering diversity and metabolic spectrum of animals in dead zone sediments

Ocean deoxygenation driven by global warming and eutrophication is a primary concern for marine life. Resistant animals may be present in dead zone sediments, however there is lack of information on their diversity and metabolism. Here we combined geochemistry, microscopy, and RNA-seq for estimating taxonomy and functionality of micrometazoans along an oxygen gradient in the largest dead zone in the world. Nematodes are metabolically active at oxygen concentrations below 1.8 µmol L-1, and their diversity and community structure are different between low oxygen areas. This is likely due to toxic hydrogen sulfide and its potential to be oxidized by oxygen or nitrate. Zooplankton resting stages dominate the metazoan community, and these populations possibly use cytochrome c oxidase as an oxygen sensor to exit dormancy. Our study sheds light on mechanisms of animal adaptation to extreme environments. These biological resources can be essential for recolonization of dead zones when oxygen conditions improve.

SRSF7 maintains its homeostasis through the expression of Split-ORFs and nuclear body assembly

SRSF7 is an essential RNA-binding protein whose misexpression promotes cancer. Here, we describe how SRSF7 maintains its protein homeostasis in murine P19 cells using an intricate negative feedback mechanism. SRSF7 binding to its premessenger RNA promotes inclusion of a poison cassette exon and transcript degradation via nonsense-mediated decay (NMD). However, elevated SRSF7 levels inhibit NMD and promote translation of two protein halves, termed Split-ORFs, from the bicistronic SRSF7-PCE transcript. The first half acts as dominant-negative isoform suppressing poison cassette exon inclusion and instead promoting the retention of flanking introns containing repeated SRSF7 binding sites. Massive SRSF7 binding to these sites and its oligomerization promote the assembly of large nuclear bodies, which sequester SRSF7 transcripts at their transcription site, preventing their export and restoring normal SRSF7 protein levels. We further show that hundreds of human and mouse NMD targets, especially RNA-binding proteins, encode potential Split-ORFs, some of which are expressed under specific cellular conditions.

Endogenous murine microbiota member Faecalibaculum rodentium and its human homologue protect from intestinal tumour growth

The microbiota has been shown to promote intestinal tumourigenesis, but a possible anti-tumourigenic effect has also been postulated. Here, we demonstrate that changes in the microbiota and mucus composition are concomitant with tumourigenesis. We identified two anti-tumourigenic strains of the microbiota-Faecalibaculum rodentium and its human homologue, Holdemanella biformis-that are strongly under-represented during tumourigenesis. Reconstitution of ApcMin/+ or azoxymethane- and dextran sulfate sodium-treated mice with an isolate of F. rodentium (F. PB1) or its metabolic products reduced tumour growth. Both F. PB1 and H. biformis produced short-chain fatty acids that contributed to control protein acetylation and tumour cell proliferation by inhibiting calcineurin and NFATc3 activation in mouse and human settings. We have thus identified endogenous anti-tumourigenic bacterial strains with strong diagnostic, therapeutic and translational potential.

Integrating microarray-based spatial transcriptomics and single-cell RNA-seq reveals tissue architecture in pancreatic ductal adenocarcinomas

Single-cell RNA sequencing (scRNA-seq) enables the systematic identification of cell populations in a tissue, but characterizing their spatial organization remains challenging. We combine a microarray-based spatial transcriptomics method that reveals spatial patterns of gene expression using an array of spots, each capturing the transcriptomes of multiple adjacent cells, with scRNA-Seq generated from the same sample. To annotate the precise cellular composition of distinct tissue regions, we introduce a method for multimodal intersection analysis. Applying multimodal intersection analysis to primary pancreatic tumors, we find that subpopulations of ductal cells, macrophages, dendritic cells and cancer cells have spatially restricted enrichments, as well as distinct coenrichments with other cell types. Furthermore, we identify colocalization of inflammatory fibroblasts and cancer cells expressing a stress-response gene module. Our approach for mapping the architecture of scRNA-seq-defined subpopulations can be applied to reveal the interactions inherent to complex tissues.

New insight into the molecular mechanism of colour differentiation among floral segments in orchids

An unbalanced pigment distribution among the sepal and petal segments results in various colour patterns of orchid flowers. Here, we explored this type of mechanism of colour pattern formation in flowers of the Cattleya hybrid 'KOVA'. Our study showed that pigment accumulation displayed obvious spatiotemporal specificity in the flowers and was likely regulated by three R2R3-MYB transcription factors. Before flowering, RcPAP1 was specifically expressed in the epichile to activate the anthocyanin biosynthesis pathway, which caused substantial cyanin accumulation and resulted in a purple-red colour. After flowering, the expression of RcPAP2 resulted in a low level of cyanin accumulation in the perianths and a pale pink colour, whereas RcPCP1 was expressed only in the hypochile, where it promoted α-carotene and lutein accumulation and resulted in a yellow colour. Additionally, we propose that the spatiotemporal expression of different combinations of AP3- and AGL6-like genes might participate in KOVA flower colour pattern formation.

A comprehensive non-redundant gene catalog reveals extensive within-community intraspecies diversity in the human vagina

Analysis of metagenomic and metatranscriptomic data is complicated and typically requires extensive computational resources. Leveraging a curated reference database of genes encoded by members of the target microbiome can make these analyses more tractable. In this study, we assemble a comprehensive human vaginal non-redundant gene catalog (VIRGO) that includes 0.95 million non-redundant genes. The gene catalog is functionally and taxonomically annotated. We also construct a vaginal orthologous groups (VOG) from VIRGO. The gene-centric design of VIRGO and VOG provides an easily accessible tool to comprehensively characterize the structure and function of vaginal metagenome and metatranscriptome datasets. To highlight the utility of VIRGO, we analyze 1,507 additional vaginal metagenomes, and identify a high degree of intraspecies diversity within and across vaginal microbiota. VIRGO offers a convenient reference database and toolkit that will facilitate a more in-depth understanding of the role of vaginal microorganisms in women’s health and reproductive outcomes.

Reference databases are essential for studies on host-microbiota interactions. Here, the authors present the construction of VIRGO, a human vaginal non-redundant gene catalog, which represents a comprehensive resource for taxonomic and functional profiling of vaginal microbiomes from metagenomic and metatranscriptomic datasets.

Candidate silencer elements for the human and mouse genomes

The study of gene regulation is dominated by a focus on the control of gene activation or increase in the level of expression. Just as critical is the process of gene repression or silencing. Chromatin signatures have identified enhancers, however, genome-wide identification of silencers by computational or experimental approaches are lacking. Here, we first define uncharacterized cis-regulatory elements likely containing silencers and find that 41.5% of ~7500 tested elements show silencer activity using massively parallel reporter assay (MPRA). We trained a support vector machine classifier based on MPRA data to predict candidate silencers in over 100 human and mouse cell or tissue types. The predicted candidate silencers exhibit characteristics expected of silencers. Leveraging promoter-capture HiC data, we find that over 50% of silencers are interacting with gene promoters having very low to no expression. Our results suggest a general strategy for genome-wide identification and characterization of silencer elements.

Identification of silencer elements by computational or experimental approaches in a genome-wide manner is still challenging. Here authors define uncharacterized cis-regulatory elements (CREs) in human and mouse genomes likely containing silencer elements, and test them in cells using massively parallel reporter assays to identify silencer elements that showed silencer activity.

Regulatory sites for splicing in human basal ganglia are enriched for disease-relevant information

Genome-wide association studies have generated an increasing number of common genetic variants associated with neurological and psychiatric disease risk. An improved understanding of the genetic control of gene expression in human brain is vital considering this is the likely modus operandum for many causal variants. However, human brain sampling complexities limit the explanatory power of brain-related expression quantitative trait loci (eQTL) and allele-specific expression (ASE) signals. We address this, using paired genomic and transcriptomic data from putamen and substantia nigra from 117 human brains, interrogating regulation at different RNA processing stages and uncovering novel transcripts. We identify disease-relevant regulatory loci, find that splicing eQTLs are enriched for regulatory information of neuron-specific genes, that ASEs provide cell-specific regulatory information with evidence for cellular specificity, and that incomplete annotation of the brain transcriptome limits interpretation of risk loci for neuropsychiatric disease. This resource of regulatory data is accessible through our web server, http://braineacv2.inf.um.es/.

Regulation of gene expression and splicing are thought to be tissue-specific. Here, the authors obtain genomic and transcriptomic data from putamen and substantia nigra of 117 neurologically healthy human brains and find that splicing eQTLs are enriched for neuron-specific regulatory information.

A novel DNA methylation motif identified in Bacillus pumilus BA06 and possible roles in the regulation of gene expression

Single-molecule real-time (SMRT) sequencing can be used to identify a wide variety of chemical modifications of the genome, such as methylation. Here, we applied this approach to identify N6-methyl-adenine (m6A) and N4-methyl-cytosine (m4C) modification in the genome of Bacillus pumilus BA06. A typical methylation recognition motif of the type I restriction-modification system (R-M), 5'-TC^m6AN8TTGG-3'/3'-AGTN8^m6AACC-5', was identified. We confirmed that this motif was a new type I methylation site using REBASE analysis and that it was recognized by a type I R-M system, Bpu6ORFCP, according to methylation sensitivity assays in vivo and vitro. Furthermore, we found that deletion of the R-M system Bpu6ORFCP induced transcriptional changes in many genes and led to increased gene expression in pathways related to ABC transporters, sulfur metabolism, ribosomes, cysteine and methionine metabolism and starch and sucrose metabolism, suggesting that the R-M system in B. pumilus BA06 has other significant biological functions beyond protecting the B. pumilus BA06 genome from foreign DNA.

Bacterial coexistence driven by motility and spatial competition

Elucidating elementary mechanisms that underlie bacterial diversity is central to ecology^1,2 and microbiome research³. Bacteria are known to coexist by metabolic specialization⁴, cooperation⁵ and cyclic warfare^6-8. Many species are also motile⁹, which is studied in terms of mechanism^10,11, benefit^12,13, strategy^14,15, evolution^16,17 and ecology^18,19. Indeed, bacteria often compete for nutrient patches that become available periodically or by random disturbances^2,20,21. However, the role of bacterial motility in coexistence remains unexplored experimentally. Here we show that-for mixed bacterial populations that colonize nutrient patches-either population outcompetes the other when low in relative abundance. This inversion of the competitive hierarchy is caused by active segregation and spatial exclusion within the patch: a small fast-moving population can outcompete a large fast-growing population by impeding its migration into the patch, while a small fast-growing population can outcompete a large fast-moving population by expelling it from the initial contact area. The resulting spatial segregation is lost for weak growth-migration trade-offs and a lack of virgin space, but is robust to population ratio, density and chemotactic ability, and is observed in both laboratory and wild strains. These findings show that motility differences and their trade-offs with growth are sufficient to promote diversity, and suggest previously undescribed roles for motility in niche formation and collective expulsion-containment strategies beyond individual search and survival.

ATM and PRDM9 regulate SPO11-bound recombination intermediates during meiosis

Meiotic recombination is initiated by SPO11-induced double-strand breaks (DSBs). In most mammals, the methyltransferase PRDM9 guides SPO11 targeting, and the ATM kinase controls meiotic DSB numbers. Following MRE11 nuclease removal of SPO11, the DSB is resected and loaded with DMC1 filaments for homolog invasion. Here, we demonstrate the direct detection of meiotic DSBs and resection using END-seq on mouse spermatocytes with low sample input. We find that DMC1 limits both minimum and maximum resection lengths, whereas 53BP1, BRCA1 and EXO1 play surprisingly minimal roles. Through enzymatic modifications to END-seq, we identify a SPO11-bound meiotic recombination intermediate (SPO11-RI) present at all hotspots. We propose that SPO11-RI forms because chromatin-bound PRDM9 asymmetrically blocks MRE11 from releasing SPO11. In Atm-/- spermatocytes, trapped SPO11 cleavage complexes accumulate due to defective MRE11 initiation of resection. Thus, in addition to governing SPO11 breakage, ATM and PRDM9 are critical local regulators of mammalian SPO11 processing.

Nucleosome destabilization by nuclear non-coding RNAs

In the nucleus, genomic DNA is wrapped around histone octamers to form nucleosomes. In principle, nucleosomes are substantial barriers to transcriptional activities. Nuclear non-coding RNAs (ncRNAs) are proposed to function in chromatin conformation modulation and transcriptional regulation. However, it remains unclear how ncRNAs affect the nucleosome structure. Eleanors are clusters of ncRNAs that accumulate around the estrogen receptor-α (ESR1) gene locus in long-term estrogen deprivation (LTED) breast cancer cells, and markedly enhance the transcription of the ESR1 gene. Here we detected nucleosome depletion around the transcription site of Eleanor2, the most highly expressed Eleanor in the LTED cells. We found that the purified Eleanor2 RNA fragment drastically destabilized the nucleosome in vitro. This activity was also exerted by other ncRNAs, but not by poly(U) RNA or DNA. The RNA-mediated nucleosome destabilization may be a common feature among natural nuclear RNAs, and may function in transcription regulation in chromatin.

The Eleanor cluster of non-coding RNAs is localised upstream of estrogen receptor-α (ESR1) gene locus in estrogen-deprived breast cancer cells. Fujita et al find that RNA fragments of Eleanor2 and of other non-coding RNAs are able to destabilise nucleosomes in vitro, suggesting a role in transcriptional regulation.

H3K4 methylation at active genes mitigates transcription-replication conflicts during replication stress

Transcription-replication conflicts (TRCs) occur when intensive transcriptional activity compromises replication fork stability, potentially leading to gene mutations. Transcription-deposited H3K4 methylation (H3K4me) is associated with regions that are susceptible to TRCs; however, the interplay between H3K4me and TRCs is unknown. Here we show that H3K4me aggravates TRC-induced replication failure in checkpoint-defective cells, and the presence of methylated H3K4 slows down ongoing replication. Both S-phase checkpoint activity and H3K4me are crucial for faithful DNA synthesis under replication stress, especially in highly transcribed regions where the presence of H3K4me is highest and TRCs most often occur. H3K4me mitigates TRCs by decelerating ongoing replication, analogous to how speed bumps slow down cars. These findings establish the concept that H3K4me defines the transcriptional status of a genomic region and defends the genome from TRC-mediated replication stress and instability.

Transcription-replication conflicts (TRC) can contribute to genome instability. Here the authors reveal that under replication stress H3K4 methylation can play a role in TRC prevention.

A combination strategy targeting enhancer plasticity exerts synergistic lethality against BETi-resistant leukemia cells

Primary and acquired drug resistance imposes a major threat to achieving optimized clinical outcomes during cancer treatment. Aberrant changes in epigenetic modifications are closely involved in drug resistance of tumor cells. Using BET inhibitor (BETi) resistant leukemia cells as a model system, we demonstrated herein that genome-wide enhancer remodeling played a pivotal role in driving therapeutic resistance via compensational re-expression of pro-survival genes. Capitalizing on the CRISPR interference technology, we identified the second intron of IncRNA, PVT1, as a unique bona fide gained enhancer that restored MYC transcription independent of BRD4 recruitment in leukemia. A combined BETi and CDK7 inhibitor treatment abolished MYC transcription by impeding RNAPII loading without affecting PVT1-mediated chromatin looping at the MYC locus in BETi-resistant leukemia cells. Together, our findings have established the feasibility of targeting enhancer plasticity to overcome drug resistance associated with epigenetic therapies.

Characterization and potential mechanisms of highly antibiotic tolerant VBNC Escherichia coli induced by low level chlorination

Escherichia coli is an important pathogenic indicator in drinking water. Viable but non-culturable (VBNC) E. coli induced by low level chlorination was found to have higher antibiotic tolerance. The emerging of VBNC bacteria in drinking water systems is posing challenges to the control of bio-safety. It is necessary to study the underlying mechanisms of VBNC state E. coli under actual residual chlorine condition of drinking water pipe network. In this study, we investigated the changes of morphology and gene expressions that might present such state. The results indicated that the size of VBNC E. coli was not remarkably changed or recovered culturability under favorable environmental conditions. Results from transcriptomic analysis revealed that the regulated genes related to fimbrial-like adhesin protein, putative periplasmic pilin chaperone, regulators of the transcriptional regulation, antibiotic resistance genes and stress-induced genes, rendering VBNC cells more tolerant to adverse environmental conditions. In total of 16 genes were significantly up-regulated under the VBNC state, including three genes encoding toxic protein (ygeG, ibsD, shoB), indicating that VBNC E. coil was still a threat to human. The work is of great relevance in the context of better understanding this poorly understood physiological state.

The microbiota programs DNA methylation to control intestinal homeostasis and inflammation

Although much research has been done on the diversity of the gut microbiome, little is known about how it influences intestinal homeostasis under normal and pathogenic conditions. Epigenetic mechanisms have recently been suggested to operate at the interface between the microbiota and the intestinal epithelium. We performed whole-genome bisulfite sequencing on conventionally raised and germ-free mice, and discovered that exposure to commensal microbiota induced localized DNA methylation changes at regulatory elements, which are TET2/3-dependent. This culminated in the activation of a set of 'early sentinel' response genes to maintain intestinal homeostasis. Furthermore, we demonstrated that exposure to the microbiota in dextran sodium sulfate-induced acute inflammation results in profound DNA methylation and chromatin accessibility changes at regulatory elements, leading to alterations in gene expression programs enriched in colitis- and colon-cancer-associated functions. Finally, by employing genetic interventions, we show that microbiota-induced epigenetic programming is necessary for proper intestinal homeostasis in vivo.

The structural basis for cohesin-CTCF-anchored loops

Cohesin catalyses the folding of the genome into loops that are anchored by CTCF1. The molecular mechanism of how cohesin and CTCF structure the 3D genome has remained unclear. Here we show that a segment within the CTCF N terminus interacts with the SA2-SCC1 subunits of human cohesin. We report a crystal structure of SA2-SCC1 in complex with CTCF at a resolution of 2.7 Å, which reveals the molecular basis of the interaction. We demonstrate that this interaction is specifically required for CTCF-anchored loops and contributes to the positioning of cohesin at CTCF binding sites. A similar motif is present in a number of established and newly identified cohesin ligands, including the cohesin release factor WAPL2,3. Our data suggest that CTCF enables the formation of chromatin loops by protecting cohesin against loop release. These results provide fundamental insights into the molecular mechanism that enables the dynamic regulation of chromatin folding by cohesin and CTCF.

ARID1A influences HDAC1/BRD4 activity, intrinsic proliferative capacity and breast cancer treatment response

Using genome-wide clustered regularly interspaced short palindromic repeats (CRISPR) screens to understand endocrine drug resistance, we discovered ARID1A and other SWI/SNF complex components as the factors most critically required for response to two classes of estrogen receptor-alpha (ER) antagonists. In this context, SWI/SNF-specific gene deletion resulted in drug resistance. Unexpectedly, ARID1A was also the top candidate in regard to response to the bromodomain and extraterminal domain inhibitor JQ1, but in the opposite direction, with loss of ARID1A sensitizing breast cancer cells to bromodomain and extraterminal domain inhibition. We show that ARID1A is a repressor that binds chromatin at ER cis-regulatory elements. However, ARID1A elicits repressive activity in an enhancer-specific, but forkhead box A1-dependent and active, ER-independent manner. Deletion of ARID1A resulted in loss of histone deacetylase 1 binding, increased histone 4 lysine acetylation and subsequent BRD4-driven transcription and growth. ARID1A mutations are more frequent in treatment-resistant disease, and our findings provide mechanistic insight into this process while revealing rational treatment strategies for these patients.

Clades of huge phages from across Earth's ecosystems

Bacteriophages typically have small genomes1 and depend on their bacterial hosts for replication2. Here we sequenced DNA from diverse ecosystems and found hundreds of phage genomes with lengths of more than 200 kilobases (kb), including a genome of 735 kb, which is-to our knowledge-the largest phage genome to be described to date. Thirty-five genomes were manually curated to completion (circular and no gaps). Expanded genetic repertoires include diverse and previously undescribed CRISPR-Cas systems, transfer RNAs (tRNAs), tRNA synthetases, tRNA-modification enzymes, translation-initiation and elongation factors, and ribosomal proteins. The CRISPR-Cas systems of phages have the capacity to silence host transcription factors and translational genes, potentially as part of a larger interaction network that intercepts translation to redirect biosynthesis to phage-encoded functions. In addition, some phages may repurpose bacterial CRISPR-Cas systems to eliminate competing phages. We phylogenetically define the major clades of huge phages from human and other animal microbiomes, as well as from oceans, lakes, sediments, soils and the built environment. We conclude that the large gene inventories of huge phages reflect a conserved biological strategy, and that the phages are distributed across a broad bacterial host range and across Earth's ecosystems.

Evidence of increased hypoxia signaling in fetal liver from maternal nutrient restriction in mice

BACKGROUND

Intrauterine growth restriction (IUGR) is a pregnancy condition where fetal growth is reduced, and offspring from IUGR pregnancies are at increased risk for type II diabetes as adults. The liver is susceptible to fetal undernutrition experienced by IUGR infants and animal models of growth restriction. This study aimed to examine hepatic expression changes in a maternal nutrient restriction (MNR) mouse model of IUGR to understand fetal adaptations that influence adult metabolism.

METHODS

Liver samples of male offspring from MNR (70% of ad libitum starting at E6.5) or control pregnancies were obtained at E18.5 and differential expression was assessed by RNAseq and western blots.

RESULTS

Forty-nine differentially expressed (FDR < 0.1) transcripts were enriched in hypoxia-inducible pathways including Fkbp5 (1.6-fold change), Ccng2 (1.5-fold change), Pfkfb3 (1.5-fold change), Kdm3a (1.2-fold change), Btg2 (1.6-fold change), Vhl (1.3-fold change), and Hif-3a (1.3-fold change) (FDR < 0.1). Fkbp5, Pfkfb3, Kdm3a, and Hif-3a were confirmed by qPCR, but only HIF-2a (2.2-fold change, p = 0.002) and HIF-3a (1.3 p = 0.03) protein were significantly increased.

CONCLUSION

Although a moderate impact, these data support evidence of fetal adaptation to reduced nutrients by increased hypoxia signaling in the liver.