Abstract Neurofibromatosis type 1 (NF1), a neuroectodermal disorder, is caused by germline mutations in the NF1 gene. NF1 affects approximately 1/3,000 individuals worldwide, with about 50% of cases representing de novo mutations. Although the NF1 gene was identified in 1990, the underlying gene mutations still remain undetected in a small but obdurate minority of NF1 patients. We postulated that in these patients, hitherto undetected pathogenic mutations might occur in regulatory elements far upstream of the NF1 gene. In an attempt to identify such remotely acting regulatory elements, we reasoned that some of them might reside within DNA sequences that (1) have the potential to interact at distance with the NF1 gene and (2) lie within a histone H3K27ac-enriched region, a characteristic of active enhancers. Combining Hi-C data, obtained by means of the chromosome conformation capture technique, with data on the location and level of histone H3K27ac enrichment upstream of the NF1 gene, we predicted in silico the presence of two remotely acting regulatory regions, located, respectively, approximately 600 kb and approximately 42 kb upstream of the NF1 gene. These regions were then sequenced in 47 NF1 patients in whom no mutations had been found in either the NF1 or SPRED1 gene regions. Five patients were found to harbour DNA sequence variants in the distal H3K27ac-enriched region. Although these variants are of uncertain pathological significance and still remain to be functionally characterized, this approach promises to be of general utility for the detection of mutations underlying other inherited disorders that may be caused by mutations in remotely acting regulatory elements.
The cytosine-guanine (CpG) dinucleotide has long been known to be a hotspot for pathological mutation in the human genome. This hypermutability is related to its role as the major site of cytosine methylation with the attendant risk of spontaneous deamination of 5-methylcytosine (5mC) to yield thymine. Cytosine methylation, however, also occurs in the context of CpNpG sites in the human genome, an unsurprising finding since the intrinsic symmetry of CpNpG renders it capable of supporting a semi-conservative model of replication of the methylation pattern. Recently, it has become clear that significant DNA methylation occurs in a CpHpG context (where H = A, C or T) in a variety of human somatic tissues. If we assume that CpHpG methylation also occurs in the germline, and that 5mC deamination can occur within a CpHpG context, then we might surmise that methylated CpHpG sites could also constitute mutation hotspots causing human genetic disease. To test this postulate, 54,625 missense and nonsense mutations from 2,113 genes causing inherited disease were retrieved from the Human Gene Mutation Database http://www.hgmd.org . Some 18.2 per cent of these pathological lesions were found to be C → T and G → A transitions located in CpG dinucleotides (compatible with a model of methylation-mediated deamination of 5mC), an approximately ten-fold higher proportion than would have been expected by chance alone. The corresponding proportion for the CpHpG trinucleotide was 9.9 per cent, an approximately two-fold higher proportion than would have been expected by chance. We therefore estimate that ~5 per cent of missense/nonsense mutations causing human inherited disease may be attributable to methylation-mediated deamination of 5mC within a CpHpG context.
Abstract Understanding structure--function relationships in the brain remains an important challenge in neuroscience. However, whilst structural brain networks are intrinsically directed, due to the prevalence of chemical synapses in the cortex, most studies in network neuroscience represent the brain as an undirected network. Here, we explore the role that directionality plays in shaping transition dynamics of functional brain states. Using a system of Hopfield neural elements with heterogeneous structural connectivity given by different species and parcellations (cat, Caenorhabditis elegans and two macaque networks), we investigate the effect of removing directionality of connections on brain capacity, which we quantify via its ability to store attractor states. In addition to determining large numbers of fixed-point attractor sets, we deploy the recently developed basin stability technique in order to assess the global stability of such brain states, which can be considered a proxy for network state robustness. Our study indicates that not only can directed network topology have a significant effect on the information capacity of connectome-based networks, but it can also impact significantly the domains of attraction of the aforementioned brain states. In particular, we find network modularity to be a key mechanism underlying the formation of neural activity patterns, and moreover, our results suggest that neglecting network directionality has the scope to eliminate states that correlate highly with the directed modular structure of the brain. A numerical analysis of the distribution of attractor states identified a small set of prototypical direction-dependent activity patterns that potentially constitute a `skeleton' of the non-stationary dynamics typically observed in the brain. This study thereby emphasizes the substantial role network directionality can have in shaping the brain's ability to both store and process information.
Translocations and gross deletions are important causes of both cancer and inherited disease. Such gene rearrangements are nonrandomly distributed in the human genome as a consequence of selection for growth advantage and/or the inherent potential of some DNA sequences to be frequently involved in breakage and recombination. Using the Gross Rearrangement Breakpoint Database [GRaBD; www.uwcm.ac.uk/uwcm/mg/grabd/grabd.html] (containing 397 germ-line and somatic DNA breakpoint junction sequences derived from 219 different rearrangements underlying human inherited disease and cancer), we have analyzed the sequence context of translocation and deletion breakpoints in a search for general characteristics that might have rendered these sequences prone to rearrangement. The oligonucleotide composition of breakpoint junctions and a set of reference sequences, matched for length and genomic location, were compared with respect to their nucleotide composition. Deletion breakpoints were found to be AT-rich whereas by comparison, translocation breakpoints were GC-rich. Alternating purine-pyrimidine sequences were found to be significantly over-represented in the vicinity of deletion breakpoints while polypyrimidine tracts were over-represented at translocation breakpoints. A number of recombination-associated motifs were found to be over-represented at translocation breakpoints (including DNA polymerase pause sites/frameshift hotspots, immunoglobulin heavy chain class switch sites, heptamer/nonamer V(D)J recombination signal sequences, translin binding sites, and the chi element) but, with the exception of the translin-binding site and immunoglobulin heavy chain class switch sites, none of these motifs were over-represented at deletion breakpoints. Alu sequences were found to span both breakpoints in seven cases of gross deletion that may thus be inferred to have arisen by homologous recombination. Our results are therefore consistent with a role for homologous unequal recombination in deletion mutagenesis and a role for nonhomologous recombination in the generation of translocations.
Nonallelic homologous recombination (NAHR) is one of the major mechanisms underlying copy number variation in the human genome. Although several disease-associated meiotic NAHR breakpoints have been analyzed in great detail, hotspots for mitotic NAHR are not well characterized. Type-2 NF1 microdeletions, which are predominantly of postzygotic origin, constitute a highly informative model with which to investigate the features of mitotic NAHR. Here, a custom-designed MLPA- and PCR-based approach was used to identify 23 novel NAHR-mediated type-2 NF1 deletions. Breakpoint analysis of these 23 type-2 deletions, together with 17 NAHR-mediated type-2 deletions identified previously, revealed that the breakpoints are nonuniformly distributed within the paralogous SUZ12 and SUZ12P sequences. Further, the analysis of this large group of type-2 deletions revealed breakpoint recurrence within short segments (ranging in size from 57 to 253-bp) as well as the existence of a novel NAHR hotspot of 1.9-kb (termed PRS4). This hotspot harbored 20% (8/40) of the type-2 deletion breakpoints and contains the 253-bp recurrent breakpoint region BR6 in which four independent type-2 deletion breakpoints were identified. Our findings indicate that a combination of an open chromatin conformation and short non-B DNA-forming repeats may predispose to recurrent mitotic NAHR events between SUZ12 and its pseudogene.
Mucolipidosis type III (MLIII) is an autosomal recessive disorder affecting lysosomal hydrolase trafficking. In a study of 10 patients from seven families with a clinical phenotype and enzymatic diagnosis of MLIII, six novel GNPTG gene mutations were identified. These included missense (p.T286M) and nonsense (p.W111X) mutations and a transition in the obligate AG-dinucleotide of the intron 8 acceptor splice site (c.610-2A>G). Three microdeletions were also identified, two of which (c.611delG and c.640_667del28) were located within the coding region whereas one (c.609+28_610-16del) was located entirely within intron 8. RT-PCR analysis of the c.610-2A>G transition demonstrated that the change altered splicing, leading to the production of two distinct aberrantly spliced forms, viz. the skipping of exon 9 (p.G204_K247del) or the retention of introns 8 and 9 (p.G204VfsX28). RT-PCR analysis, performed on a patient homozygous for the intronic deletion (c.609+28_610-16del), failed to detect any GNPTG RNA transcripts. To determine whether c.609+28_610-16del allele-derived transcripts were subject to nonsense-mediated mRNA decay (NMD), patient fibroblasts were incubated with the protein synthesis inhibitor anisomycin. An RT-PCR fragment retaining 43 bp of intron 8 was consistently detected suggesting that the 33-bp genomic deletion had elicited NMD. Quantitative real-time PCR and GNPTG western blot analysis confirmed that the homozygous microdeletion p.G204VfsX17 had elicited NMD resulting in failure to synthesize GNPTG protein. Analysis of the sequences surrounding the microdeletion breakpoints revealed either intrinsic repetitivity of the deleted region or short direct repeats adjacent to the breakpoint junctions. This is consistent with these repeats having mediated the microdeletions via replication slippage and supports the view that the mutational spectrum of the GNPTG gene is strongly influenced by the properties of the local DNA sequence environment.