Deciphering D4Z4 CpG methylation gradients in fascioscapulohumeral muscular dystrophy using nanopore sequencing.

Fascioscapulohumeral muscular dystrophy (FSHD) is caused by a unique genetic mechanism that relies on contraction and hypomethylation of the D4Z4 macrosatellite array on the Chromosome 4q telomere allowing ectopic expression of the DUX4 gene in skeletal muscle. Genetic analysis is difficult because of the large size and repetitive nature of the array, a nearly identical array on the 10q telomere, and the presence of divergent D4Z4 arrays scattered throughout the genome. Here, we combine nanopore long-read sequencing with Cas9-targeted enrichment of 4q and 10q D4Z4 arrays for comprehensive genetic analysis including determination of the length of the 4q and 10q D4Z4 arrays with base-pair resolution. In the same assay, we differentiate 4q from 10q telomeric sequences, determine A/B haplotype, identify paralogous D4Z4 sequences elsewhere in the genome, and estimate methylation for all CpGs in the array. Asymmetric, length-dependent methylation gradients were observed in the 4q and 10q D4Z4 arrays that reach a hypermethylation point at approximately 10 D4Z4 repeat units, consistent with the known threshold of pathogenic D4Z4 contractions. High resolution analysis of individual D4Z4 repeat methylation revealed areas of low methylation near the CTCF/insulator region and areas of high methylation immediately preceding the DUX4 transcriptional start site. Within the DUX4 exons, we observed a waxing/waning methylation pattern with a 180-nucleotide periodicity, consistent with phased nucleosomes. Targeted nanopore sequencing complements recently developed molecular combing and optical mapping approaches to genetic analysis for FSHD by adding precision of the length measurement, base-pair resolution sequencing, and quantitative methylation analysis.

Evidence-based path to newborn screening for Duchenne muscular dystrophy.

OBJECTIVE: Creatine kinase (CK) levels are increased on dried blood spots in newborns related to the birthing process. As a marker for newborn screening, CK in Duchenne muscular dystrophy (DMD) results in false-positive testing. In this report, we introduce a 2-tier system using the dried blood spot to first assess CK with follow-up DMD gene testing. METHODS: A fluorometric assay based upon the enzymatic transphosphorylation of adenosine diphosphate to adenosine triphosphate was used to measure CK activity. Preliminary studies established a population-based range of CK in newborns using 30,547 deidentified anonymous dried blood spot samples. Mutation analysis used genomic DNA extracted from the dried blood spot followed by whole genome amplification with assessment of single-/multiexon deletions/duplications in the DMD gene using multiplex ligation-dependent probe amplification. RESULTS: DMD gene mutations (all exonic deletions) were found in 6 of 37,649 newborn male subjects, all of whom had CK levels>2,000U/l. In 3 newborns with CK>2,000U/l in whom DMD gene abnormalities were not found, we identified limb-girdle muscular dystrophy gene mutations affecting DYSF, SGCB, and FKRP. INTERPRETATION: A 2-tier system of analysis for newborn screening for DMD has been established. This path for newborn screening fits our health care system, minimizes false-positive testing, and uses predetermined levels of CK on dried blood spots to predict DMD gene mutations.

Deciphering D4Z4 CpG methylation gradients in fascioscapulohumeral muscular dystrophy using nanopore sequencing.

Fascioscapulohumeral muscular dystrophy (FSHD) is caused by a unique genetic mechanism that relies on contraction and hypomethylation of the D4Z4 macrosatellite array on the chromosome 4q telomere allowing ectopic expression of the DUX4 gene in skeletal muscle. Genetic analysis is difficult due to the large size and repetitive nature of the array, a nearly identical array on the 10q telomere, and the presence of divergent D4Z4 arrays scattered throughout the genome. Here, we combine nanopore long-read sequencing with Cas9-targeted enrichment of 4q and 10q D4Z4 arrays for comprehensive genetic analysis including determination of the length of the 4q and 10q D4Z4 arrays with base-pair resolution. In the same assay, we differentiate 4q from 10q telomeric sequences, determine A/B haplotype, identify paralogous D4Z4 sequences elsewhere in the genome, and estimate methylation for all CpGs in the array. Asymmetric, length-dependent methylation gradients were observed in the 4q and 10q D4Z4 arrays that reach a hypermethylation point at approximately 10 D4Z4 repeat units, consistent with the known threshold of pathogenic D4Z4 contractions. High resolution analysis of individual D4Z4 repeat methylation revealed areas of low methylation near the CTCF/insulator region and areas of high methylation immediately preceding the DUX4 transcriptional start site. Within the DUX4 exons, we observed a waxing/waning methylation pattern with a 180-nucleotide periodicity, consistent with phased nucleosomes. Targeted nanopore sequencing complements recently developed molecular combing and optical mapping approaches to genetic analysis for FSHD by adding precision of the length measurement, base-pair resolution sequencing, and quantitative methylation analysis.

Population-Based Prevalence of Myotonic Dystrophy Type 1 Using Genetic Analysis of Statewide Blood Screening Program.

OBJECTIVE: To determine whether the genetic prevalence of the CTG expansion in the DMPK gene associated with myotonic dystrophy type 1 (DM1) in an unbiased cohort is higher than previously reported population estimates, ranging from 5 to 20 per 100,000 individuals. METHODS: This study used a cross-sectional cohort of deidentified dried blood spots from the newborn screening program in the state of New York, taken from consecutive births from 2013 to 2014. Blood spots were screened for the CTG repeat expansion in the DMPK gene using triplet-repeat primed PCR and melt curve analysis. Melt curve morphology was assessed by 4 blinded reviewers to identify samples with possible CTG expansion. Expansion of the CTG repeat was validated by PCR fragment sizing using capillary electrophoresis for samples classified as positive or premutation to confirm the result. Prevalence was calculated as the number of samples with CTG repeat size ≥50 repeats compared to the overall cohort. RESULTS: Of 50,382 consecutive births, there were 24 with a CTG repeat expansion ≥50, consistent with a diagnosis of DM1. This represents a significantly higher DM1 prevalence of 4.76 per 10,000 births (95% confidence interval 2.86-6.67) or 1 in every 2,100 births. There were an additional 96 samples (19.1 per 10,000 or 1 in 525 births) with a CTG expansion in the DMPK gene in the premutation range (CTG)35-49. CONCLUSION: The prevalence of individuals with CTG repeat expansions in DMPK is up to 5 times higher than previous reported estimates. This suggests that DM1, with multisystemic manifestations, is likely underdiagnosed in practice.

High throughput screening for expanded CTG repeats in myotonic dystrophy type 1 using melt curve analysis.

BACKGROUND: Myotonic dystrophy type 1 (DM1) is caused by CTG repeat expansions in the DMPK gene and is the most common form of muscular dystrophy. Patients can have long delays from onset to diagnosis, since clinical signs and symptoms are often nonspecific and overlapping with other disorders. Clinical genetic testing by Southern blot or triplet-primed PCR (TP-PCR) is technically challenging and cost prohibitive for population surveys. METHODS: Here, we present a high throughput, low-cost screening tool for CTG repeat expansions using TP-PCR followed by high resolution melt curve analysis with saturating concentrations of SYBR GreenER dye. RESULTS: We determined that multimodal melt profiles from the TP-PCR assay are a proxy for amplicon length stoichiometry. In a screen of 10,097 newborn blood spots, melt profile analysis accurately reflected the tri-modal distribution of common alleles from 5 to 35 CTG repeats, and identified the premutation and full expansion alleles. CONCLUSION: We demonstrate that robust detection of expanded CTG repeats in a single tube can be achieved from samples derived from specimens with minimal template DNA such as dried blood spots (DBS). This technique is readily adaptable to large-scale testing programs such as population studies and newborn screening programs.

Genome degeneration and adaptation in a nascent stage of symbiosis.

Symbiotic associations between animals and microbes are ubiquitous in nature, with an estimated 15% of all insect species harboring intracellular bacterial symbionts. Most bacterial symbionts share many genomic features including small genomes, nucleotide composition bias, high coding density, and a paucity of mobile DNA, consistent with long-term host association. In this study, we focus on the early stages of genome degeneration in a recently derived insect-bacterial mutualistic intracellular association. We present the complete genome sequence and annotation of Sitophilus oryzae primary endosymbiont (SOPE). We also present the finished genome sequence and annotation of strain HS, a close free-living relative of SOPE and other insect symbionts of the Sodalis-allied clade, whose gene inventory is expected to closely resemble the putative ancestor of this group. Structural, functional, and evolutionary analyses indicate that SOPE has undergone extensive adaptation toward an insect-associated lifestyle in a very short time period. The genome of SOPE is large in size when compared with many ancient bacterial symbionts; however, almost half of the protein-coding genes in SOPE are pseudogenes. There is also evidence for relaxed selection on the remaining intact protein-coding genes. Comparative analyses of the whole-genome sequence of strain HS and SOPE highlight numerous genomic rearrangements, duplications, and deletions facilitated by a recent expansion of insertions sequence elements, some of which appear to have catalyzed adaptive changes. Functional metabolic predictions suggest that SOPE has lost the ability to synthesize several essential amino acids and vitamins. Analyses of the bacterial cell envelope and genes encoding secretion systems suggest that these structures and elements have become simplified in the transition to a mutualistic association.