Deleterious, protein-altering variants in the transcriptional coregulator ZMYM3 in 27 individuals with a neurodevelopmental delay phenotype.

Neurodevelopmental disorders (NDDs) result from highly penetrant variation in hundreds of different genes, some of which have not yet been identified. Using the MatchMaker Exchange, we assembled a cohort of 27 individuals with rare, protein-altering variation in the transcriptional coregulator ZMYM3, located on the X chromosome. Most (n = 24) individuals were males, 17 of which have a maternally inherited variant; six individuals (4 male, 2 female) harbor de novo variants. Overlapping features included developmental delay, intellectual disability, behavioral abnormalities, and a specific facial gestalt in a subset of males. Variants in almost all individuals (n = 26) are missense, including six that recurrently affect two residues. Four unrelated probands were identified with inherited variation affecting Arg441, a site at which variation has been previously seen in NDD-affected siblings, and two individuals have de novo variation resulting in p.Arg1294Cys (c.3880C>T). All variants affect evolutionarily conserved sites, and most are predicted to damage protein structure or function. ZMYM3 is relatively intolerant to variation in the general population, is widely expressed across human tissues, and encodes a component of the KDM1A-RCOR1 chromatin-modifying complex. ChIP-seq experiments on one variant, p.Arg1274Trp, indicate dramatically reduced genomic occupancy, supporting a hypomorphic effect. While we are unable to perform statistical evaluations to definitively support a causative role for variation in ZMYM3, the totality of the evidence, including 27 affected individuals, recurrent variation at two codons, overlapping phenotypic features, protein-modeling data, evolutionary constraint, and experimentally confirmed functional effects strongly support ZMYM3 as an NDD-associated gene.

DNMT1 modulates gene expression without its catalytic activity partially through its interactions with histone-modifying enzymes.

While DNA methyltransferase1 (DNMT1) is classically known for its functions as a maintenance methyltransferase enzyme, additional roles for DNMT1 in gene expression are not as clearly understood. Several groups have shown that deletion of the catalytic domain from DNMT1 does not abolish repressive activity of the protein against a reporter gene. In our studies, we examine the repressor function of catalytically inactive DNMT1 at endogenous genes. First, potential DNMT1 target genes were identified by searching for genes up-regulated in HCT116 colon cancer cells genetically disrupted for DNMT1 (DNMT1(-/-) hypomorph cells). Next, the requirement for DNMT1 activity for repression of these genes was assessed by stably restoring expression of wild-type or catalytically inactive DNMT1. Both wild-type and mutant proteins are able to occupy the promoters and repress the expression of a set of target genes, and induce, at these promoters, both the depletion of active histone marks and the recruitment of a H3K4 demethylase, KDM1A/LSD1. Together, our findings show that there are genes for which DNMT1 acts as a transcriptional repressor independent from its methyltransferase function and that this repressive function may invoke a role for a scaffolding function of the protein at target genes.

An enhancer RNA recruits KMT2A to regulate transcription of Myb.

The Myb proto-oncogene encodes the transcription factor c-MYB, which is critical for hematopoiesis. Distant enhancers of Myb form a hub of interactions with the Myb promoter. We identified a long non-coding RNA (Myrlin) originating from the -81-kb murine Myb enhancer. Myrlin and Myb are coordinately regulated during erythroid differentiation. Myrlin TSS deletion using CRISPR-Cas9 reduced Myrlin and Myb expression and LDB1 complex occupancy at the Myb enhancers, compromising enhancer contacts and reducing RNA Pol II occupancy in the locus. In contrast, CRISPRi silencing of Myrlin left LDB1 and the Myb enhancer hub unperturbed, although Myrlin and Myb expressions were downregulated, decoupling transcription and chromatin looping. Myrlin interacts with the KMT2A/MLL1 complex. Myrlin CRISPRi compromised KMT2A occupancy in the Myb locus, decreasing CDK9 and RNA Pol II binding and resulting in Pol II pausing in the Myb first exon/intron. Thus, Myrlin directly participates in activating Myb transcription by recruiting KMT2A.

An enhancer RNA recruits MLL1 to regulate transcription of Myb.

The Myb proto-oncogene encodes the transcription factor c-MYB, which is critical for hematopoiesis. Distant enhancers of Myb form a hub of interactions with the Myb promoter. We identified a long non-coding RNA (Myrlin) originating from the -81 kb murine Myb enhancer. Myrlin and Myb are coordinately regulated during erythroid differentiation. Myrlin TSS deletion using CRISPR/Cas9 reduced Myrlin and Myb expression and LDB1 complex occupancy at the Myb enhancers, compromising enhancer contacts and reducing RNA Pol II occupancy in the locus. In contrast, CRISPRi silencing of Myrlin left LDB1 and the Myb enhancer hub unperturbed, although Myrlin and Myb expression were downregulated, decoupling transcription and chromatin looping. Myrlin interacts with the MLL1 complex. Myrlin CRISPRi compromised MLL1 occupancy in the Myb locus, decreasing CDK9 and RNA Pol II binding and resulting in Pol II pausing in the Myb first exon/intron. Thus, Myrlin directly participates in activating Myb transcription by recruiting MLL1.

Within-host evolution of Burkholderia pseudomallei in four cases of acute melioidosis.

Little is currently known about bacterial pathogen evolution and adaptation within the host during acute infection. Previous studies of Burkholderia pseudomallei, the etiologic agent of melioidosis, have shown that this opportunistic pathogen mutates rapidly both in vitro and in vivo at tandemly repeated loci, making this organism a relevant model for studying short-term evolution. In the current study, B. pseudomallei isolates cultured from multiple body sites from four Thai patients with disseminated melioidosis were subjected to fine-scale genotyping using multilocus variable-number tandem repeat analysis (MLVA). In order to understand and model the in vivo variable-number tandem repeat (VNTR) mutational process, we characterized the patterns and rates of mutations in vitro through parallel serial passage experiments of B. pseudomallei. Despite the short period of infection, substantial divergence from the putative founder genotype was observed in all four melioidosis cases. This study presents a paradigm for examining bacterial evolution over the short timescale of an acute infection. Further studies are required to determine whether the mutational process leads to phenotypic alterations that impact upon bacterial fitness in vivo. Our findings have important implications for future sampling strategies, since colonies in a single clinical sample may be genetically heterogeneous, and organisms in a culture taken late in the infective process may have undergone considerable genetic change compared with the founder inoculum.

Phylogeographic reconstruction of a bacterial species with high levels of lateral gene transfer.

BACKGROUND: Phylogeographic reconstruction of some bacterial populations is hindered by low diversity coupled with high levels of lateral gene transfer. A comparison of recombination levels and diversity at seven housekeeping genes for eleven bacterial species, most of which are commonly cited as having high levels of lateral gene transfer shows that the relative contributions of homologous recombination versus mutation for Burkholderia pseudomallei is over two times higher than for Streptococcus pneumoniae and is thus the highest value yet reported in bacteria. Despite the potential for homologous recombination to increase diversity, B. pseudomallei exhibits a relative lack of diversity at these loci. In these situations, whole genome genotyping of orthologous shared single nucleotide polymorphism loci, discovered using next generation sequencing technologies, can provide very large data sets capable of estimating core phylogenetic relationships. We compared and searched 43 whole genome sequences of B. pseudomallei and its closest relatives for single nucleotide polymorphisms in orthologous shared regions to use in phylogenetic reconstruction. RESULTS: Bayesian phylogenetic analyses of >14,000 single nucleotide polymorphisms yielded completely resolved trees for these 43 strains with high levels of statistical support. These results enable a better understanding of a separate analysis of population differentiation among >1,700 B. pseudomallei isolates as defined by sequence data from seven housekeeping genes. We analyzed this larger data set for population structure and allele sharing that can be attributed to lateral gene transfer. Our results suggest that despite an almost panmictic population, we can detect two distinct populations of B. pseudomallei that conform to biogeographic patterns found in many plant and animal species. That is, separation along Wallace's Line, a biogeographic boundary between Southeast Asia and Australia. CONCLUSION: We describe an Australian origin for B. pseudomallei, characterized by a single introduction event into Southeast Asia during a recent glacial period, and variable levels of lateral gene transfer within populations. These patterns provide insights into mechanisms of genetic diversification in B. pseudomallei and its closest relatives, and provide a framework for integrating the traditionally separate fields of population genetics and phylogenetics for other bacterial species with high levels of lateral gene transfer.

Identification of melioidosis outbreak by multilocus variable number tandem repeat analysis.

Endemic melioidosis is caused by genetically diverse Burkholderia pseudomallei strains. However, clonal outbreaks (multiple cases caused by 1 strain) have occurred, such as from contaminated potable water. B. pseudomallei is designated a group B bioterrorism agent, which necessitates rapidly recognizing point-source outbreaks. Pulsed-field gel electrophoresis (PFGE) and multilocus sequence typing (MLST) can identify genetically related isolates, but results take several days to obtain. We developed a simplified 4-locus multilocus variable number tandem repeat analysis (MLVA-4) for rapid typing and compared results with PFGE and MLST for a large number of well-characterized B. pseudomallei isolates. MLVA-4 compared favorably with MLST and PFGE for the same isolates; it discriminated between 65 multilocus sequence types and showed relatedness between epidemiologically linked isolates from outbreak clusters and between isolates from individual patients. MLVA-4 can establish or refute that a clonal outbreak of melioidosis has occurred within 8 hours of receipt of bacterial strains.

Genomic islands from five strains of Burkholderia pseudomallei.

BACKGROUND: Burkholderia pseudomallei is the etiologic agent of melioidosis, a significant cause of morbidity and mortality where this infection is endemic. Genomic differences among strains of B. pseudomallei are predicted to be one of the major causes of the diverse clinical manifestations observed among patients with melioidosis. The purpose of this study was to examine the role of genomic islands (GIs) as sources of genomic diversity in this species. RESULTS: We found that genomic islands (GIs) vary greatly among B. pseudomallei strains. We identified 71 distinct GIs from the genome sequences of five reference strains of B. pseudomallei: K96243, 1710b, 1106a, MSHR668, and MSHR305. The genomic positions of these GIs are not random, as many of them are associated with tRNA gene loci. In particular, the 3' end sequences of tRNA genes are predicted to be involved in the integration of GIs. We propose the term "tRNA-mediated site-specific recombination" (tRNA-SSR) for this mechanism. In addition, we provide a GI nomenclature that is based upon integration hotspots identified here or previously described. CONCLUSION: Our data suggest that acquisition of GIs is one of the major sources of genomic diversity within B. pseudomallei and the molecular mechanisms that facilitate horizontally-acquired GIs are common across multiple strains of B. pseudomallei. The differential presence of the 71 GIs across multiple strains demonstrates the importance of these mobile elements for shaping the genetic composition of individual strains and populations within this bacterial species.

A KDM5 Inhibitor Increases Global H3K4 Trimethylation Occupancy and Enhances the Biological Efficacy of 5-Aza-2'-Deoxycytidine.

The H3K4 demethylase KDM5B is amplified and overexpressed in luminal breast cancer, suggesting it might constitute a potential cancer therapy target. Here, we characterize, in breast cancer cells, the molecular effects of a recently developed small-molecule inhibitor of the KDM5 family of proteins (KDM5i), either alone or in combination with the DNA-demethylating agent 5-aza-2'-deoxycytidine (DAC). KDM5i treatment alone increased expression of a small number of genes, whereas combined treatment with DAC enhanced the effects of the latter for increasing expression of hundreds of DAC-responsive genes. ChIP-seq studies revealed that KDM5i resulted in the broadening of existing H3K4me3 peaks. Furthermore, cells treated with the drug combination exhibited increased promoter and gene body H3K4me3 occupancy at DAC-responsive genes compared with DAC alone. Importantly, treatment with either DAC or DAC+KDM5i induced a dramatic increase in H3K27ac at enhancers with an associated significant increase in target gene expression, suggesting a previously unappreciated effect of DAC on transcriptional regulation. KDM5i synergized with DAC to reduce the viability of luminal breast cancer cells in in vitro assays. Our study provides the first look into the molecular effects of a novel KDM5i compound and suggests that combinatorial inhibition along with DAC represents a new area to explore in translational epigenetics.Significance: This study offers a first look into the molecular effects of a novel KDM5 inhibitory compound, suggesting how its use in combination with DNA methylation inhibitors presents new opportunities to explore in translational cancer epigenetics. Cancer Res; 78(5); 1127-39. ©2017 AACR.

VNTR analysis of selected outbreaks of Burkholderia pseudomallei in Australia.

Molecular typing methods for Burkholderia pseudomallei have been successful at assigning isolates into epidemiologically related groups, but have not been able to detect differences and define evolutionary patterns within groups. Our variable number tandem repeat (VNTR) analysis of a set of 121 Australian B. pseudomallei isolates, 104 of which were associated with nine epidemiological groups, provides fine scale differentiation even among very closely related isolates. We used a Bayesian model based upon mutation accumulation patterns to define the close phylogenetic relationships within these epidemiological groups. Our results reveal that genetic diversity can exist within a very small geographic area, and that low levels of diversity can exist even within a single infection. These methods provide the ability to generate robust evolutionary hypotheses that enable tracking of B. pseudomallei in forensic and epidemiological outbreaks at fine phylogenetic scales.

Fine-scale genetic diversity among Burkholderia pseudomallei soil isolates in northeast Thailand.

Burkholderia pseudomallei soil isolates from northeast Thailand were genotyped using multiple-locus variable-number tandem repeat (VNTR) analysis (MLVA) and multilocus sequence typing (MLST). MLVA identified 19 genotypes within three clades, while MLST revealed two genotypes. These close genetic relationships imply a recent colonization followed by localized expansion, similar to what occurs in an outbreak situation.