Missense and truncating variants in CHD5 in a dominant neurodevelopmental disorder with intellectual disability, behavioral disturbances, and epilepsy.

Located in the critical 1p36 microdeletion region, the chromodomain helicase DNA-binding protein 5 (CHD5) gene encodes a subunit of the nucleosome remodeling and deacetylation (NuRD) complex required for neuronal development. Pathogenic variants in six of nine chromodomain (CHD) genes cause autosomal dominant neurodevelopmental disorders, while CHD5-related disorders are still unknown. Thanks to GeneMatcher and international collaborations, we assembled a cohort of 16 unrelated individuals harboring heterozygous CHD5 variants, all identified by exome sequencing. Twelve patients had de novo CHD5 variants, including ten missense and two splice site variants. Three familial cases had nonsense or missense variants segregating with speech delay, learning disabilities, and/or craniosynostosis. One patient carried a frameshift variant of unknown inheritance due to unavailability of the father. The most common clinical features included language deficits (81%), behavioral symptoms (69%), intellectual disability (64%), epilepsy (62%), and motor delay (56%). Epilepsy types were variable, with West syndrome observed in three patients, generalized tonic-clonic seizures in two, and other subtypes observed in one individual each. Our findings suggest that, in line with other CHD-related disorders, heterozygous CHD5 variants are associated with a variable neurodevelopmental syndrome that includes intellectual disability with speech delay, epilepsy, and behavioral problems as main features.

wg-blimp: an end-to-end analysis pipeline for whole genome bisulfite sequencing data.

BACKGROUND: Analysing whole genome bisulfite sequencing datasets is a data-intensive task that requires comprehensive and reproducible workflows to generate valid results. While many algorithms have been developed for tasks such as alignment, comprehensive end-to-end pipelines are still sparse. Furthermore, previous pipelines lack features or show technical deficiencies, thus impeding analyses. RESULTS: We developed wg-blimp (whole genome bisulfite sequencing methylation analysis pipeline) as an end-to-end pipeline to ease whole genome bisulfite sequencing data analysis. It integrates established algorithms for alignment, quality control, methylation calling, detection of differentially methylated regions, and methylome segmentation, requiring only a reference genome and raw sequencing data as input. Comparing wg-blimp to previous end-to-end pipelines reveals similar setups for common sequence processing tasks, but shows differences for post-alignment analyses. We improve on previous pipelines by providing a more comprehensive analysis workflow as well as an interactive user interface. To demonstrate wg-blimp's ability to produce correct results we used it to call differentially methylated regions for two publicly available datasets. We were able to replicate 112 of 114 previously published regions, and found results to be consistent with previous findings. We further applied wg-blimp to a publicly available sample of embryonic stem cells to showcase methylome segmentation. As expected, unmethylated regions were in close proximity of transcription start sites. Segmentation results were consistent with previous analyses, despite different reference genomes and sequencing techniques. CONCLUSIONS: wg-blimp provides a comprehensive analysis pipeline for whole genome bisulfite sequencing data as well as a user interface for simplified result inspection. We demonstrated its applicability by analysing multiple publicly available datasets. Thus, wg-blimp is a relevant alternative to previous analysis pipelines and may facilitate future epigenetic research.

Identification and characterisation of pathogenic and non-pathogenic FGF14 repeat expansions.

Repeat expansions in FGF14 cause autosomal dominant late-onset cerebellar ataxia (SCA27B) with estimated pathogenic thresholds of 250 (incomplete penetrance) and 300 AAG repeats (full penetrance), but the sequence of pathogenic and non-pathogenic expansions remains unexplored. Here, we demonstrate that STRling and ExpansionHunter accurately detect FGF14 expansions from short-read genome data using outlier approaches. By combining long-range PCR and nanopore sequencing in 169 patients with cerebellar ataxia and 802 controls, we compare FGF14 expansion alleles, including interruptions and flanking regions. Uninterrupted AAG expansions are significantly enriched in patients with ataxia from a lower threshold (180-200 repeats) than previously reported based on expansion size alone. Conversely, AAGGAG hexameric expansions are equally frequent in patients and controls. Distinct 5' flanking regions, interruptions and pre-repeat sequences correlate with repeat size. Furthermore, pure AAG (pathogenic) and AAGGAG (non-pathogenic) repeats form different secondary structures. Regardless of expansion size, SCA27B is a recognizable clinical entity characterized by frequent episodic ataxia and downbeat nystagmus, similar to the presentation observed in a family with a previously unreported nonsense variant (SCA27A). Overall, this study suggests that SCA27B is a major overlooked cause of adult-onset ataxia, accounting for 23-31% of unsolved patients. We strongly recommend re-evaluating pathogenic thresholds and integrating expansion sequencing into the molecular diagnostic process.

De novo variants in the non-coding spliceosomal snRNA gene RNU4-2 are a frequent cause of syndromic neurodevelopmental disorders.

Around 60% of individuals with neurodevelopmental disorders (NDD) remain undiagnosed after comprehensive genetic testing, primarily of protein-coding genes1. Increasingly, large genome-sequenced cohorts are improving our ability to discover new diagnoses in the non-coding genome. Here, we identify the non-coding RNA RNU4-2 as a novel syndromic NDD gene. RNU4-2 encodes the U4 small nuclear RNA (snRNA), which is a critical component of the U4/U6.U5 tri-snRNP complex of the major spliceosome2. We identify an 18 bp region of RNU4-2 mapping to two structural elements in the U4/U6 snRNA duplex (the T-loop and Stem III) that is severely depleted of variation in the general population, but in which we identify heterozygous variants in 119 individuals with NDD. The vast majority of individuals (77.3%) have the same highly recurrent single base-pair insertion (n.64_65insT). We estimate that variants in this region explain 0.41% of individuals with NDD. We demonstrate that RNU4-2 is highly expressed in the developing human brain, in contrast to its contiguous counterpart RNU4-1 and other U4 homologs, supporting RNU4-2's role as the primary U4 transcript in the brain. Overall, this work underscores the importance of non-coding genes in rare disorders. It will provide a diagnosis to thousands of individuals with NDD worldwide and pave the way for the development of effective treatments for these individuals.

Correction: Association of Rare Genetic Variants in Opioid Receptors with Tourette Syndrome.

[This corrects the article DOI: 10.5334/tohm.464.].

Early-set POMC methylation variability is accompanied by increased risk for obesity and is addressable by MC4R agonist treatment.

Increasing evidence points toward epigenetic variants as a risk factor for developing obesity. We analyzed DNA methylation of the POMC (pro-opiomelanocortin) gene, which is pivotal for satiety regulation. We identified sex-specific and nongenetically determined POMC hypermethylation associated with a 1.4-fold (confidence interval, 1.03 to 2.04) increased individual risk of developing obesity. To investigate the early embryonic establishment of POMC methylation states, we established a human embryonic stem cell (hESC) model. Here, hESCs (WA01) were transferred into a naïve state, which was associated with a reduction of DNA methylation. Naïve hESCs were differentiated via a formative state into POMC-expressing hypothalamic neurons, which was accompanied by re-establishment of DNA methylation patterning. We observed that reduced POMC gene expression was associated with increased POMC methylation in POMC-expressing neurons. On the basis of these findings, we treated POMC-hypermethylated obese individuals (n = 5) with an MC4R agonist and observed a body weight reduction of 4.66 ± 2.16% (means ± SD) over a mean treatment duration of 38.4 ± 26.0 weeks. In summary, we identified an epigenetic obesity risk variant at the POMC gene fulfilling the criteria for a metastable epiallele established in early embryonic development that may be addressable by MC4R agonist treatment to reduce body weight.

Genome-wide case-only analysis of gene-gene interactions with known Parkinson's disease risk variants reveals link between LRRK2 and SYT10.

The effects of one genetic factor upon Parkinson's disease (PD) risk may be modified by other genetic factors. Such gene-gene interaction (G×G) could explain some of the 'missing heritability' of PD and the reduced penetrance of known PD risk variants. Using the largest single nucleotide polymorphism (SNP) genotype data set currently available for PD (18,688 patients), provided by the International Parkinson's Disease Genomics Consortium, we studied G×G with a case-only (CO) design. To this end, we paired each of 90 SNPs previously reported to be associated with PD with one of 7.8 million quality-controlled SNPs from a genome-wide panel. Support of any putative G×G interactions found was sought by the analysis of independent genotype-phenotype and experimental data. A total of 116 significant pairwise SNP genotype associations were identified in PD cases, pointing towards G×G. The most prominent associations involved a region on chromosome 12q containing SNP rs76904798, which is a non-coding variant of the LRRK2 gene. It yielded the lowest interaction p-value overall with SNP rs1007709 in the promoter region of the SYT10 gene (interaction OR = 1.80, 95% CI: 1.65-1.95, p = 2.7 × 10-43). SNPs around SYT10 were also associated with the age-at-onset of PD in an independent cohort of carriers of LRRK2 mutation p.G2019S. Moreover, SYT10 gene expression during neuronal development was found to differ between cells from affected and non-affected p.G2019S carriers. G×G interaction on PD risk, involving the LRRK2 and SYT10 gene regions, is biologically plausible owing to the known link between PD and LRRK2, its involvement in neural plasticity, and the contribution of SYT10 to the exocytosis of secretory vesicles in neurons.

No evidence for intervention-associated DNA methylation changes in monocytes of patients with posttraumatic stress disorder.

DNA methylation patterns can be responsive to environmental influences. This observation has sparked interest in the potential for psychological interventions to influence epigenetic processes. Recent studies have observed correlations between DNA methylation changes and therapy outcome. However, most did not control for changes in cell composition. This study had two aims: first, we sought to replicate therapy-associated changes in DNA methylation of commonly assessed candidate genes in isolated monocytes from 60 female patients with post-traumatic stress disorder (PTSD). Our second, exploratory goal was to identify novel genomic regions with substantial pre-to-post intervention DNA methylation changes by performing whole-genome bisulfite sequencing (WGBS) in two patients with PTSD. Equivalence testing and Bayesian analyses provided evidence against physiologically meaningful intervention-associated DNA methylation changes in monocytes of PTSD patients in commonly investigated target genes (NR3C1, FKBP5, SLC6A4, OXTR). Furthermore, WGBS yielded only a limited set of candidate regions with suggestive evidence of differential DNA methylation pre- to post-therapy. These differential DNA methylation patterns did not prove replicable when investigated in the entire cohort. We conclude that there is no evidence for major, recurrent intervention-associated DNA methylation changes in the investigated genes in monocytes of patients with PTSD.

Systematic analysis and prediction of genes associated with monogenic disorders on human chromosome X.

Disease gene discovery on chromosome (chr) X is challenging owing to its unique modes of inheritance. We undertook a systematic analysis of human chrX genes. We observe a higher proportion of disorder-associated genes and an enrichment of genes involved in cognition, language, and seizures on chrX compared to autosomes. We analyze gene constraints, exon and promoter conservation, expression, and paralogues, and report 127 genes sharing one or more attributes with known chrX disorder genes. Using machine learning classifiers trained to distinguish disease-associated from dispensable genes, we classify 247 genes, including 115 of the 127, as having high probability of being disease-associated. We provide evidence of an excess of variants in predicted genes in existing databases. Finally, we report damaging variants in CDK16 and TRPC5 in patients with intellectual disability or autism spectrum disorders. This study predicts large-scale gene-disease associations that could be used for prioritization of X-linked pathogenic variants.

PHIP-associated Chung-Jansen syndrome: Report of 23 new individuals.

In 2016 and 2018, Chung, Jansen and others described a new syndrome caused by haploinsufficiency of PHIP (pleckstrin homology domain interacting protein, OMIM *612,870) and mainly characterized by developmental delay (DD), learning difficulties/intellectual disability (ID), behavioral abnormalities, facial dysmorphism and obesity (CHUJANS, OMIM #617991). So far, PHIP alterations appear to be a rare cause of DD/ID. "Omics" technologies such as exome sequencing or array analyses have led to the identification of distinct types of alterations of PHIP, including, truncating variants, missense substitutions, splice variants and large deletions encompassing portions of the gene or the entire gene as well as adjacent genomic regions. We collected clinical and genetic data of 23 individuals with PHIP-associated Chung-Jansen syndrome (CHUJANS) from all over Europe. Follow-up investigations (e.g. Sanger sequencing, qPCR or Fluorescence-in-situ-Hybridization) and segregation analysis showed either de novo occurrence or inheritance from an also (mildly) affected parent. In accordance with previously described patients, almost all individuals reported here show developmental delay (22/23), learning disability or ID (22/23), behavioral abnormalities (20/23), weight problems (13/23) and characteristic craniofacial features (i.e. large ears/earlobes, prominent eyebrows, anteverted nares and long philtrum (23/23)). To further investigate the facial gestalt of individuals with CHUJANS, we performed facial analysis using the GestaltMatcher approach. By this, we could establish that PHIP patients are indistinguishable based on the type of PHIP alteration (e.g. missense, loss-of-function, splice site) but show a significant difference to the average face of healthy individuals as well as to individuals with Prader-Willi syndrome (PWS, OMIM #176270) or with a CUL4B-alteration (Intellectual developmental disorder, X-linked, syndromic, Cabezas type, OMIM #300354). Our findings expand the mutational and clinical spectrum of CHUJANS. We discuss the molecular and clinical features in comparison to the published individuals. The fact that some variants were inherited from a mildly affected parent further illustrates the variability of the associated phenotype and outlines the importance of a thorough clinical evaluation combined with genetic analyses for accurate diagnosis and counselling.

The different clinical facets of SYN1-related neurodevelopmental disorders.

Synapsin-I (SYN1) is a presynaptic phosphoprotein crucial for synaptogenesis and synaptic plasticity. Pathogenic SYN1 variants are associated with variable X-linked neurodevelopmental disorders mainly affecting males. In this study, we expand on the clinical and molecular spectrum of the SYN1-related neurodevelopmental disorders by describing 31 novel individuals harboring 22 different SYN1 variants. We analyzed newly identified as well as previously reported individuals in order to define the frequency of key features associated with these disorders. Specifically, behavioral disturbances such as autism spectrum disorder or attention deficit hyperactivity disorder are observed in 91% of the individuals, epilepsy in 82%, intellectual disability in 77%, and developmental delay in 70%. Seizure types mainly include tonic-clonic or focal seizures with impaired awareness. The presence of reflex seizures is one of the most representative clinical manifestations related to SYN1. In more than half of the cases, seizures are triggered by contact with water, but other triggers are also frequently reported, including rubbing with a towel, fever, toothbrushing, fingernail clipping, falling asleep, and watching others showering or bathing. We additionally describe hyperpnea, emotion, lighting, using a stroboscope, digestive troubles, and defecation as possible triggers in individuals with SYN1 variants. The molecular spectrum of SYN1 variants is broad and encompasses truncating variants (frameshift, nonsense, splicing and start-loss variants) as well as non-truncating variants (missense substitutions and in-frame duplications). Genotype-phenotype correlation revealed that epileptic phenotypes are enriched in individuals with truncating variants. Furthermore, we could show for the first time that individuals with early seizures onset tend to present with severe-to-profound intellectual disability, hence highlighting the existence of an association between early seizure onset and more severe impairment of cognitive functions. Altogether, we present a detailed clinical description of the largest series of individuals with SYN1 variants reported so far and provide the first genotype-phenotype correlations for this gene. A timely molecular diagnosis and genetic counseling are cardinal for appropriate patient management and treatment.

Whole-genome methylation analysis of testicular germ cells from cryptozoospermic men points to recurrent and functionally relevant DNA methylation changes.

BACKGROUND: Several studies have reported an association between male infertility and aberrant sperm DNA methylation patterns, in particular in imprinted genes. In a recent investigation based on whole methylome and deep bisulfite sequencing, we have not found any evidence for such an association, but have demonstrated that somatic DNA contamination and genetic variation confound methylation studies in sperm of severely oligozoospermic men. To find out whether testicular germ cells (TGCs) of such patients might carry aberrant DNA methylation, we compared the TGC methylomes of four men with cryptozoospermia (CZ) and four men with obstructive azoospermia, who had normal spermatogenesis and served as controls (CTR). RESULTS: There was no difference in DNA methylation at the whole genome level or at imprinted regions between CZ and CTR samples. However, using stringent filters to identify group-specific methylation differences, we detected 271 differentially methylated regions (DMRs), 238 of which were hypermethylated in CZ (binominal test, p < 2.2 × 10-16). The DMRs were enriched for distal regulatory elements (p = 1.0 × 10-6) and associated with 132 genes, 61 of which are differentially expressed at various stages of spermatogenesis. Almost all of the 67 DMRs associated with the 61 genes (94%) are hypermethylated in CZ (63/67, p = 1.107 × 10-14). As judged by single-cell RNA sequencing, 13 DMR-associated genes, which are mainly expressed during meiosis and spermiogenesis, show a significantly different pattern of expression in CZ patients. In four of these genes, the promoter is hypermethylated in CZ men, which correlates with a lower expression level in these patients. In the other nine genes, eight of which downregulated in CZ, germ cell-specific enhancers may be affected. CONCLUSIONS: We found that impaired spermatogenesis is associated with DNA methylation changes in testicular germ cells at functionally relevant regions of the genome. We hypothesize that the described DNA methylation changes may reflect or contribute to premature abortion of spermatogenesis and therefore not appear in the mature, motile sperm.

Lasp1 regulates adherens junction dynamics and fibroblast transformation in destructive arthritis.

The LIM and SH3 domain protein 1 (Lasp1) was originally cloned from metastatic breast cancer and characterised as an adaptor molecule associated with tumourigenesis and cancer cell invasion. However, the regulation of Lasp1 and its function in the aggressive transformation of cells is unclear. Here we use integrative epigenomic profiling of invasive fibroblast-like synoviocytes (FLS) from patients with rheumatoid arthritis (RA) and from mouse models of the disease, to identify Lasp1 as an epigenomically co-modified region in chronic inflammatory arthritis and a functionally important binding partner of the Cadherin-11/ß-Catenin complex in zipper-like cell-to-cell contacts. In vitro, loss or blocking of Lasp1 alters pathological tissue formation, migratory behaviour and platelet-derived growth factor response of arthritic FLS. In arthritic human TNF transgenic mice, deletion of Lasp1 reduces arthritic joint destruction. Therefore, we show a function of Lasp1 in cellular junction formation and inflammatory tissue remodelling and identify Lasp1 as a potential target for treating inflammatory joint disorders associated with aggressive cellular transformation.

The sperm epigenome does not display recurrent epimutations in patients with severely impaired spermatogenesis.

BACKGROUND: In the past 15 years, numerous studies have described aberrant DNA methylation of imprinted genes (e.g. MEST and H19) in sperm of oligozoospermic men, but the prevalence and genomic extent of abnormal methylation patterns have remained unknown. RESULTS: Using deep bisulfite sequencing (DBS), we screened swim-up sperm samples from 40 normozoospermic and 93 patients diagnosed as oligoasthenoteratozoospermic, oligoteratozoospermic or oligozoospermic, which are termed OATs throughout the manuscript, for H19 and MEST methylation. Based on this screening, we defined three patient groups: normal controls (NC), abnormally methylated oligozoospermic (AMO; n = 7) and normally methylated oligozoospermic (NMO; n = 86). Whole-genome bisulfite sequencing (WGBS) of five NC and five AMO samples revealed abnormal methylation levels of all 50 imprinting control regions in each AMO sample. To investigate whether this finding reflected epigenetic germline mosaicism or the presence of residual somatic DNA, we made a genome-wide inventory of soma-germ cell-specific DNA methylation. We found that > 2000 germ cell-specific genes are promoter-methylated in blood and that AMO samples had abnormal methylation levels at these genes, consistent with the presence of somatic cell DNA. The comparison between the five NC and six NMO samples revealed 19 differentially methylated regions (DMRs), none of which could be validated in an independent cohort of 40 men. Previous studies reported a higher incidence of epimutations at single CpG sites in the CTCF-binding region 6 of H19 in infertile patients. DBS analysis of this locus, however, revealed an association between DNA methylation levels and genotype (rs2071094), but not fertility phenotype. CONCLUSIONS: Our results suggest that somatic DNA contamination and genetic variation confound methylation studies in sperm of infertile men. While we cannot exclude the existence of rare patients with slightly abnormal sperm methylation at non-recurrent CpG sites, the prevalence of aberrant methylation in swim-up purified sperm of infertile men has likely been overestimated, which is reassuring for patients undergoing assisted reproduction.

Evidence of methanesulfonate utilizers in the Sargasso Sea metagenome.

Methanesulfonate (MSA) is one of the products of the photo-oxidation of dimethylsulfide in the atmosphere. The genes responsible for the import of MSA into the cell (msm EFGH) and for its oxidation to formaldehyde (msm ABCD) have been previously sequenced from the soil bacterium Methylosulfonomonas methylovora str. M2 while genes for an MSA monooxygenase have been sequenced from marine bacterium Marinosulfonomonas methylotropha str. TR3. We performed a sequence-based screening of the Sargasso Sea metagenome for homologues of the MSA monooxygenase (MSAMO) and MSA import genes. Our search retrieved one scaffold bearing genes with high identity to the msm ABCD cluster plus two scaffolds bearing genes highly identical to the msm EFGH operon. We increased the available data by sequencing two metagenome plasmids, which revealed more msm genes. In these three cases synteny with the original msm operons was revealed. We also retrieved several singletons showing high identity to shorter segments of the msm clusters or individual msm genes. Furthermore, a characteristic 26-aa internal spacer of the MsmA Rieske-type motif was conserved. Our findings support the case for a significant role of MSA degraders in the marine sulfur cycle and seem to suggest that they may be prominent members of the methylotrophic community in surface ocean waters.

Cyanobacterial hydrogenases: diversity, regulation and applications.

Cyanobacteria may possess two distinct nickel-iron (NiFe)-hydrogenases: an uptake enzyme found in N(2)-fixing strains, and a bidirectional one present in both non-N(2)-fixing and N(2)-fixing strains. The uptake hydrogenase (encoded by hupSL) catalyzes the consumption of the H(2) produced during N(2) fixation, while the bidirectional enzyme (hoxEFUYH) probably plays a role in fermentation and/or acts as an electron valve during photosynthesis. hupSL constitute a transcriptional unit, and are essentially transcribed under N(2)-fixing conditions. The bidirectional hydrogenase consists of a hydrogenase and a diaphorase part, and the corresponding five hox genes are not always clustered or cotranscribed. The biosynthesis/maturation of NiFe-hydrogenases is highly complex, requiring several core proteins. In cyanobacteria, the genes that are thought to affect hydrogenases pleiotropically (hyp), as well as the genes presumably encoding the hydrogenase-specific endopeptidases (hupW and hoxW) have been identified and characterized. Furthermore, NtcA and LexA have been implicated in the transcriptional regulation of the uptake and the bidirectional enzyme respectively. Recently, the phylogenetic origin of cyanobacterial and algal hydrogenases was analyzed, and it was proposed that the current distribution in cyanobacteria reflects a differential loss of genes according to their ecological needs or constraints. In addition, the possibilities and challenges of cyanobacterial-based H(2) production are addressed.

Association of Rare Genetic Variants in Opioid Receptors with Tourette Syndrome.

Background: Genes involved in Tourette syndrome (TS) remain largely unknown. We aimed to identify genetic factors contributing to TS in a French cohort of 120 individuals using a combination of hypothesis-driven and exome-sequencing approaches. Methods: We first sequenced exons of SLITRK1-6 and HDC in the TS cohort and subsequently sequenced the exome of 12 individuals harboring rare variants in these genes to find additional rare variants contributing to the disorder under the hypothesis of oligogenic inheritance. We further screened three candidate genes (OPRK1, PCDH10, and NTSR2) preferentially expressed in the basal ganglia, and three additional genes involved in neurotensin and opioid signaling (OPRM1, NTS, and NTSR1), and compared variant frequencies in TS patients and 788 matched control individuals. We also investigated the impact of altering the expression of Oprk1 in zebrafish. Results: Thirteen ultrarare missense variants of SLITRK1-6 and HDC were identified in 12 patients. Exome sequencing in these patients revealed rare possibly deleterious variants in 3,041 genes, 54 of which were preferentially expressed in the basal ganglia. Comparison of variant frequencies altering selected candidate genes in TS and control individuals revealed an excess of potentially disrupting variants in OPRK1, encoding the opioid kappa receptor, in TS patients. Accordingly, we show that downregulation of the Oprk1 orthologue in zebrafish induces a hyperkinetic phenotype in early development. Discussion: These results support a heterogeneous and complex genetic etiology of TS, possibly involving rare variants altering the opioid pathway in some individuals, which could represent a novel therapeutic target in this disorder.

Unstable TTTTA/TTTCA expansions in MARCH6 are associated with Familial Adult Myoclonic Epilepsy type 3.

Familial Adult Myoclonic Epilepsy (FAME) is a genetically heterogeneous disorder characterized by cortical tremor and seizures. Intronic TTTTA/TTTCA repeat expansions in SAMD12 (FAME1) are the main cause of FAME in Asia. Using genome sequencing and repeat-primed PCR, we identify another site of this repeat expansion, in MARCH6 (FAME3) in four European families. Analysis of single DNA molecules with nanopore sequencing and molecular combing show that expansions range from 3.3 to 14 kb on average. However, we observe considerable variability in expansion length and structure, supporting the existence of multiple expansion configurations in blood cells and fibroblasts of the same individual. Moreover, the largest expansions are associated with micro-rearrangements occurring near the expansion in 20% of cells. This study provides further evidence that FAME is caused by intronic TTTTA/TTTCA expansions in distinct genes and reveals that expansions exhibit an unexpectedly high somatic instability that can ultimately result in genomic rearrangements.

Locus-Specific DNA Methylation Analysis by Targeted Deep Bisulfite Sequencing.

DNA methylation, i.e., the methylation of cytosine at carbon atom C5, has an important role in the regulation of gene expression. The methylation status of each cytosine in a specific genomic region can be determined by targeted deep bisulfite sequencing at single-molecule resolution. Here we describe the design of PCR primers that are used to amplify specific sequences from bisulfite-converted DNA, the preparation of sequencing libraries, the sequencing of these libraries on the MiSeq system, as well as the analysis of the sequence reads. Using appropriate software tools such as amplikyzer2, it is easy to analyze complex multiplexed samples with several regions of interest, to determine the mean methylation values of all CpG dinucleotides in a region or of each CpG dinucleotide across all or selected reads, and to compare these values between different samples and between different alleles within a sample.

Regions of common inter-individual DNA methylation differences in human monocytes: genetic basis and potential function.

BACKGROUND: There is increasing evidence for inter-individual methylation differences at CpG dinucleotides in the human genome, but the regional extent and function of these differences have not yet been studied in detail. For identifying regions of common methylation differences, we used whole genome bisulfite sequencing data of monocytes from five donors and a novel bioinformatic strategy. RESULTS: We identified 157 differentially methylated regions (DMRs) with four or more CpGs, almost none of which has been described before. The DMRs fall into different chromatin states, where methylation is inversely correlated with active, but not repressive histone marks. However, methylation is not correlated with the expression of associated genes. High-resolution single nucleotide polymorphism (SNP) genotyping of the five donors revealed evidence for a role of cis-acting genetic variation in establishing methylation patterns. To validate this finding in a larger cohort, we performed genome-wide association studies (GWAS) using SNP genotypes and 450k array methylation data from blood samples of 1128 individuals. Only 30/157 (19%) DMRs include at least one 450k CpG, which shows that these arrays miss a large proportion of DNA methylation variation. In most cases, the GWAS peak overlapped the CpG position, and these regions are enriched for CREB group, NF-1, Sp100 and CTCF binding motifs. In two cases, there was tentative evidence for a trans-effect by KRAB zinc finger proteins. CONCLUSIONS: Allele-specific DNA methylation occurs in discrete chromosomal regions and is driven by genetic variation in cis and trans, but in general has little effect on gene expression.