Using de novo assembly to identify structural variation of eight complex immune system gene regions.

Driven by the necessity to survive environmental pathogens, the human immune system has evolved exceptional diversity and plasticity, to which several factors contribute including inheritable structural polymorphism of the underlying genes. Characterizing this variation is challenging due to the complexity of these loci, which contain extensive regions of paralogy, segmental duplication and high copy-number repeats, but recent progress in long-read sequencing and optical mapping techniques suggests this problem may now be tractable. Here we assess this by using long-read sequencing platforms from PacBio and Oxford Nanopore, supplemented with short-read sequencing and Bionano optical mapping, to sequence DNA extracted from CD14+ monocytes and peripheral blood mononuclear cells from a single European individual identified as HV31. We use this data to build a de novo assembly of eight genomic regions encoding four key components of the immune system, namely the human leukocyte antigen, immunoglobulins, T cell receptors, and killer-cell immunoglobulin-like receptors. Validation of our assembly using k-mer based and alignment approaches suggests that it has high accuracy, with estimated base-level error rates below 1 in 10 kb, although we identify a small number of remaining structural errors. We use the assembly to identify heterozygous and homozygous structural variation in comparison to GRCh38. Despite analyzing only a single individual, we find multiple large structural variants affecting core genes at all three immunoglobulin regions and at two of the three T cell receptor regions. Several of these variants are not accurately callable using current algorithms, implying that further methodological improvements are needed. Our results demonstrate that assessing haplotype variation in these regions is possible given sufficiently accurate long-read and associated data. Continued reductions in the cost of these technologies will enable application of these methods to larger samples and provide a broader catalogue of germline structural variation at these loci, an important step toward making these regions accessible to large-scale genetic association studies.

Genetic and functional insights into CDA-I prevalence and pathogenesis.

BACKGROUND: Congenital dyserythropoietic anaemia type I (CDA-I) is a hereditary anaemia caused by biallelic mutations in the widely expressed genes CDAN1 and C15orf41. Little is understood about either protein and it is unclear in which cellular pathways they participate. METHODS: Genetic analysis of a cohort of patients with CDA-I identifies novel pathogenic variants in both known causative genes. We analyse the mutation distribution and the predicted structural positioning of amino acids affected in Codanin-1, the protein encoded by CDAN1. Using western blotting, immunoprecipitation and immunofluorescence, we determine the effect of particular mutations on both proteins and interrogate protein interaction, stability and subcellular localisation. RESULTS: We identify six novel CDAN1 mutations and one novel mutation in C15orf41 and uncover evidence of further genetic heterogeneity in CDA-I. Additionally, population genetics suggests that CDA-I is more common than currently predicted. Mutations are enriched in six clusters in Codanin-1 and tend to affect buried residues. Many missense and in-frame mutations do not destabilise the entire protein. Rather C15orf41 relies on Codanin-1 for stability and both proteins, which are enriched in the nucleolus, interact to form an obligate complex in cells. CONCLUSION: Stability and interaction data suggest that C15orf41 may be the key determinant of CDA-I and offer insight into the mechanism underlying this disease. Both proteins share a common pathway likely to be present in a wide variety of cell types; however, nucleolar enrichment may provide a clue as to the erythroid specific nature of CDA-I. The surprisingly high predicted incidence of CDA-I suggests that better ascertainment would lead to improved patient care.

The chromatin remodelling factor ATRX suppresses R-loops in transcribed telomeric repeats.

ATRX is a chromatin remodelling factor found at a wide range of tandemly repeated sequences including telomeres (TTAGGG)n ATRX mutations are found in nearly all tumours that maintain their telomeres via the alternative lengthening of telomere (ALT) pathway, and ATRX is known to suppress this pathway. Here, we show that recruitment of ATRX to telomeric repeats depends on repeat number, orientation and, critically, on repeat transcription. Importantly, the transcribed telomeric repeats form RNA-DNA hybrids (R-loops) whose abundance correlates with the recruitment of ATRX Here, we show loss of ATRX is also associated with increased R-loop formation. Our data suggest that the presence of ATRX at telomeres may have a central role in suppressing deleterious DNA secondary structures that form at transcribed telomeric repeats, and this may account for the increased DNA damage, stalling of replication and homology-directed repair previously observed upon loss of ATRX function.

MicroRNAs of the miR-290-295 Family Maintain Bivalency in Mouse Embryonic Stem Cells.

Numerous developmentally regulated genes in mouse embryonic stem cells (ESCs) are marked by both active (H3K4me3)- and polycomb group (PcG)-mediated repressive (H3K27me3) histone modifications. This bivalent state is thought to be important for transcriptional poising, but the mechanisms that regulate bivalent genes and the bivalent state remain incompletely understood. Examining the contribution of microRNAs (miRNAs) to the regulation of bivalent genes, we found that the miRNA biogenesis enzyme DICER was required for the binding of the PRC2 core components EZH2 and SUZ12, and for the presence of the PRC2-mediated histone modification H3K27me3 at many bivalent genes. Genes that lost bivalency were preferentially upregulated at the mRNA and protein levels. Finally, reconstituting Dicer-deficient ESCs with ESC miRNAs restored bivalent gene repression and PRC2 binding at formerly bivalent genes. Therefore, miRNAs regulate bivalent genes and the bivalent state itself.

Suppression of the alternative lengthening of telomere pathway by the chromatin remodelling factor ATRX.

Fifteen per cent of cancers maintain telomere length independently of telomerase by the homologous recombination (HR)-associated alternative lengthening of telomeres (ALT) pathway. A unifying feature of these tumours are mutations in ATRX. Here we show that expression of ectopic ATRX triggers a suppression of the pathway and telomere shortening. Importantly ATRX-mediated ALT suppression is dependent on the histone chaperone DAXX. Re-expression of ATRX is associated with a reduction in replication fork stalling, a known trigger for HR and loss of MRN from telomeres. A G-quadruplex stabilizer partially reverses the effect of ATRX, inferring ATRX may normally facilitate replication through these sequences that, if they persist, promote ALT. We propose that defective telomere chromatinization through loss of ATRX promotes the persistence of aberrant DNA secondary structures, which in turn present a barrier to DNA replication, leading to replication fork stalling, collapse, HR and subsequent recombination-mediated telomere synthesis in ALT cancers.

ATRX dysfunction induces replication defects in primary mouse cells.

The chromatin remodeling protein ATRX, which targets tandem repetitive DNA, has been shown to be required for expression of the alpha globin genes, for proliferation of a variety of cellular progenitors, for chromosome congression and for the maintenance of telomeres. Mutations in ATRX have recently been identified in tumours which maintain their telomeres by a telomerase independent pathway involving homologous recombination thought to be triggered by DNA damage. It is as yet unknown whether there is a central underlying mechanism associated with ATRX dysfunction which can explain the numerous cellular phenomena observed. There is, however, growing evidence for its role in the replication of various repetitive DNA templates which are thought to have a propensity to form secondary structures. Using a mouse knockout model we demonstrate that ATRX plays a direct role in facilitating DNA replication. Ablation of ATRX alone, although leading to a DNA damage response at telomeres, is not sufficient to trigger the alternative lengthening of telomere pathway in mouse embryonic stem cells.

The ISWI and CHD1 chromatin remodelling activities influence ADH2 expression and chromatin organization.

Nucleosome remodelling complexes play a key role in gene activation in response to environmental changes by driving promoter chromatin to reach an accessible configuration. They also mediate genome-wide chromatin organization, although their role in processes other than activation-related chromatin remodelling are poorly understood. The Saccharomyces cerevisiae ADH2 gene represents an excellent model for understanding the role of chromatin structure and remodelling in gene regulation. Following glucose depletion, highly positioned promoter nucleosomes are destabilized leading to strictly regulated kinetics of transcriptional activation. Nevertheless, no chromatin remodelling activities responsible for establishing or remodelling ADH2 chromatin structure have been identified to date. Here we show that the absence of the Isw1 and Chd1 ATP-dependent chromatin remodelling activities delays the maximal expression of ADH2 without impairing the chromatin remodelling that occurs upon activation. Instead, a destabilized chromatin structure on the ADH2 coding and termination region is observed in the absence of Isw1 or Chd1 in repressing conditions. The specific Isw1 complex involved in this nucleosome repositioning is Isw1b because the deletion of Ioc2 and Ioc4, but not of Ioc3, causes the same phenotype as the deletion of Isw1. Moreover, the lack of Chd1 combined with the absence of Isw1 and Isw2 impairs nucleosome spacing along the ADH2 gene, and genome-wide in S. cerevisiae. Thus, the ISWI and Chd1 remodelling factors are not only involved in transcription-related chromatin remodelling, but also are required to maintain a specific chromatin configuration across the yeast genome.

Common chromatin architecture, common chromatin remodeling, and common transcription kinetics of Adr1-dependent genes in Saccharomyces cerevisiae.

The chromatin structure of several Saccharomyces cerevisiae ADR1-dependent genes was comparatively analyzed in vivo in order to evaluate the role of promoter architecture in transcriptional control. In repressing conditions (high glucose) a nucleosome particle always obstructs the TATA box, immediately adjacent to an upstream-located nucleosome-free region containing a cluster of Adr1 binding sites. Upon derepression the TATA box-containing nucleosome is destabilized according to a mechanism shared by all of the genes studied. The transcription factor Adr1 is always required for the observed chromatin remodeling. mRNA accumulation of all of the genes analyzed is strongly delayed in the absence of the acetyltransferase Gcn5 and is decreased in the presence of a temperature-sensitive Esa1 mutant. The results suggest that a defined promoter chromatin structure, controlled by DNA conformational features, is relevant for the activation of coregulated genes.