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.

Autosomal dominant in cis D4Z4 repeat array duplication alleles in facioscapulohumeral dystrophy.

Facioscapulohumeral dystrophy (FSHD) has a unique genetic aetiology resulting in partial chromatin relaxation of the D4Z4 macrosatellite repeat array on 4qter. This D4Z4 chromatin relaxation facilitates inappropriate expression of the transcription factor DUX4 in skeletal muscle. DUX4 is encoded by a retrogene that is embedded within the distal region of the D4Z4 repeat array. In the European population, the D4Z4 repeat array is usually organized in a single array that ranges between 8 and 100 units. D4Z4 chromatin relaxation and DUX4 derepression in FSHD is most often caused by repeat array contraction to 1-10 units (FSHD1) or by a digenic mechanism requiring pathogenic variants in a D4Z4 chromatin repressor like SMCHD1, combined with a repeat array between 8 and 20 units (FSHD2). With a prevalence of 1.5% in the European population, in cis duplications of the D4Z4 repeat array, where two adjacent D4Z4 arrays are interrupted by a spacer sequence, are relatively common but their relationship to FSHD is not well understood. In cis duplication alleles were shown to be pathogenic in FSHD2 patients; however, there is inconsistent evidence for the necessity of an SMCHD1 mutation for disease development. To explore the pathogenic nature of these alleles we compared in cis duplication alleles in FSHD patients with or without pathogenic SMCHD1 variant. For both groups we showed duplication-allele-specific DUX4 expression. We studied these alleles in detail using pulsed-field gel electrophoresis-based Southern blotting and molecular combing, emphasizing the challenges in the characterization of these rearrangements. Nanopore sequencing was instrumental to study the composition and methylation of the duplicated D4Z4 repeat arrays and to identify the breakpoints and the spacer sequence between the arrays. By comparing the composition of the D4Z4 repeat array of in cis duplication alleles in both groups, we found that specific combinations of proximal and distal repeat array sizes determine their pathogenicity. Supported by our algorithm to predict pathogenicity, diagnostic laboratories should now be furnished to accurately interpret these in cis D4Z4 repeat array duplications, alleles that can easily be missed in routine settings.

Incipient genome erosion and metabolic streamlining for antibiotic production in a defensive symbiont.

Genome erosion is a frequently observed result of relaxed selection in insect nutritional symbionts, but it has rarely been studied in defensive mutualisms. Solitary beewolf wasps harbor an actinobacterial symbiont of the genus Streptomyces that provides protection to the developing offspring against pathogenic microorganisms. Here, we characterized the genomic architecture and functional gene content of this culturable symbiont using genomics, transcriptomics, and proteomics in combination with in vitro assays. Despite retaining a large linear chromosome (7.3 Mb), the wasp symbiont accumulated frameshift mutations in more than a third of its protein-coding genes, indicative of incipient genome erosion. Although many of the frameshifted genes were still expressed, the encoded proteins were not detected, indicating post-transcriptional regulation. Most pseudogenization events affected accessory genes, regulators, and transporters, but "Streptomyces philanthi" also experienced mutations in central metabolic pathways, resulting in auxotrophies for biotin, proline, and arginine that were confirmed experimentally in axenic culture. In contrast to the strong A+T bias in the genomes of most obligate symbionts, we observed a significant G+C enrichment in regions likely experiencing reduced selection. Differential expression analyses revealed that-compared to in vitro symbiont cultures-"S. philanthi" in beewolf antennae showed overexpression of genes for antibiotic biosynthesis, the uptake of host-provided nutrients and the metabolism of building blocks required for antibiotic production. Our results show unusual traits in the early stage of genome erosion in a defensive symbiont and suggest tight integration of host-symbiont metabolic pathways that effectively grants the host control over the antimicrobial activity of its bacterial partner.

Core transcriptional networks in Williams syndrome: IGF1-PI3K-AKT-mTOR, MAPK and actin signaling at the synapse echo autism.

Gene networks for disorders of social behavior provide the mechanisms critical for identifying therapeutic targets and biomarkers. Large behavioral phenotypic effects of small human deletions make the positive sociality of Williams syndrome (WS) ideal for determining transcriptional networks for social dysfunction currently based on DNA variations for disorders such as autistic spectrum disorder (ASD) and schizophrenia (SCHZ). Consensus on WS networks has been elusive due to the need for larger cohort size, sensitive genome-wide detection and analytic tools. We report a core set of WS network perturbations in a cohort of 58 individuals (34 with typical, 6 atypical deletions and 18 controls). Genome-wide exon-level expression arrays robustly detected changes in differentially expressed gene (DEG) transcripts from WS deleted genes that ranked in the top 11 of 12 122 transcripts, validated by quantitative reverse transcription PCR, RNASeq and western blots. WS DEG's were strictly dosed in the full but not the atypical deletions that revealed a breakpoint position effect on non-deleted CLIP2, a caveat for current phenotypic mapping based on copy number variants. Network analyses tested the top WS DEG's role in the dendritic spine, employing GeneMANIA to harmonize WS DEGs with comparable query gene-sets. The results indicate perturbed actin cytoskeletal signaling analogous to the excitatory dendritic spines. Independent protein-protein interaction analyses of top WS DEGs generated a 100-node graph annotated topologically revealing three interacting pathways, MAPK, IGF1-PI3K-AKT-mTOR/insulin and actin signaling at the synapse. The results indicate striking similarity of WS transcriptional networks to genome-wide association study-based ASD and SCHZ risk suggesting common network dysfunction for these disorders of divergent sociality.

Intron mutations and early transcription termination in Duchenne and Becker muscular dystrophy.

DMD pathogenic variants for Duchenne and Becker muscular dystrophy are detectable with high sensitivity by standard clinical exome analyses of genomic DNA. However, up to 7% of DMD mutations are deep intronic and analysis of muscle-derived RNA is an important diagnostic step for patients who have negative genomic testing but abnormal dystrophin expression in muscle. In this study, muscle biopsies were evaluated from 19 patients with clinical features of a dystrophinopathy, but negative clinical DMD mutation analysis. Reverse transcription-polymerase chain reaction or high-throughput RNA sequencing methods identified 19 mutations with one of three pathogenic pseudoexon types: deep intronic point mutations, deletions or insertions, and translocations. In association with point mutations creating intronic splice acceptor sites, we observed the first examples of DMD pseudo 3'-terminal exon mutations causing high efficiency transcription termination within introns. This connection between splicing and premature transcription termination is reminiscent of U1 snRNP-mediating telescripting in sustaining RNA polymerase II elongation across large genes, such as DMD. We propose a novel classification of three distinct types of mutations identifiable by muscle RNA analysis, each of which differ in potential treatment approaches. Recognition and appropriate characterization may lead to therapies directed toward full-length dystrophin expression for some patients.

Long-range genomic regulators of THBS1 and LTBP4 modify disease severity in duchenne muscular dystrophy.

OBJECTIVE: Duchenne muscular dystrophy (DMD) is a severe X-linked recessive disease caused by loss-of-function dystrophin (DMD) mutations in boys, who typically suffer loss of ambulation by age 12. Previously, we reported that coding variants in latent transforming growth factor beta (TGFß)-binding protein 4 (LTBP4) were associated with reduced TGFß signaling and prolonged ambulation (p = 1.0 × 10-3 ) in DMD patients; this result was subsequently replicated by other groups. In this study, we evaluated whether additional DMD modifier genes are observed using whole-genome association in the original cohort. METHODS: We performed a genome-wide association study (GWAS) for single-nucleotide polymorphisms (SNPs) influencing loss of ambulation (LOA) in the same cohort of 253 DMD patients used to detect the candidate association with LTBP4 coding variants. Gene expression and chromatin interaction databases were used to fine-map association signals above the threshold for genome-wide significance. RESULTS: Despite the small sample size, two loci associated with prolonged ambulation met genome-wide significance and were tagged by rs2725797 (chr15, p = 6.6 × 10-9 ) and rs710160 (chr19, p = 4.7 × 10-8 ). Gene expression and chromatin interaction data indicated that the latter SNP tags regulatory variants of LTBP4, whereas the former SNP tags regulatory variants of thrombospondin-1 (THBS1): an activator of TGFß signaling by direct binding to LTBP4 and an inhibitor of proangiogenic nitric oxide signaling. INTERPRETATION: Together with previous evidence implicating LTBP4, the THBS1 modifier locus emphasizes the role that common regulatory variants in gene interaction networks can play in mitigating disease progression in muscular dystrophy. Ann Neurol 2018;84:234-245.

Lung eosinophilia induced by house dust mites or ovalbumin is modulated by nicotinic receptor α7 and inhibited by cigarette smoke.

Eosinophilia (EOS) is an important component of airway inflammation and hyperresponsiveness in allergic reactions including those leading to asthma. Although cigarette smoking (CS) is a significant contributor to long-term adverse outcomes in these lung disorders, there are also the curious reports of its ability to produce acute suppression of inflammatory responses including EOS through poorly understood mechanisms. One possibility is that proinflammatory processes are suppressed by nicotine in CS acting through nicotinic receptor α7 (α7). Here we addressed the role of α7 in modulating EOS with two mouse models of an allergic response: house dust mites (HDM; Dermatophagoides sp.) and ovalbumin (OVA). The influence of α7 on EOS was experimentally resolved in wild-type mice or in mice in which a point mutation of the α7 receptor (α7E260A:G) selectively restricts normal signaling of cellular responses. RNA analysis of alveolar macrophages and the distal lung epithelium indicates that normal α7 function robustly impacts gene expression in the epithelium to HDM and OVA but to different degrees. Notable was allergen-specific α7 modulation of Ccl11 and Ccl24 (eotaxins) expression, which was enhanced in HDM but suppressed in OVA EOS. CS suppressed EOS induced by both OVA and HDM, as well as the inflammatory genes involved, regardless of α7 genotype. These results suggest that EOS in response to HDM or OVA is through signaling pathways that are modulated in a cell-specific manner by α7 and are distinct from CS suppression.

Complex signatures of natural selection at GYPA.

The human MN blood group antigens are isoforms of glycophorin A (GPA) encoded by the gene, GYPA, and are the most abundant erythrocyte sialoglycoproteins. The distribution of MN antigens has been widely studied in human populations yet the evolutionary and/or demographic factors affecting population variation remain elusive. While the primary function of GPA is yet to be discovered, it serves as the major binding site for the 175-kD erythrocyte-binding antigen (EB-175) of the malarial parasite, Plasmodium falciparum, a major selective pressure in recent human history. More specifically, exon two of GYPA encodes the receptor-binding ligand to which P. falciparum binds. Accordingly, there has been keen interest in understanding what impact, if any, natural selection has had on the distribution of variation in GYPA and exon two in particular. To this end, we resequenced GYPA in individuals sampled from both P. falciparum endemic (sub-Saharan Africa and South India) and non-endemic (Europe and East Asia) regions of the world. Observed patterns of variation suggest that GYPA has been subject to balancing selection in populations living in malaria endemic areas and in Europeans, but no such evidence was found in samples from East Asia, Oceania, and the Americas. These results are consistent with malaria acting as a selective pressure on GYPA, but also suggest that another selective force has resulted in a similar pattern of variation in Europeans. Accordingly, GYPA has perhaps a more complex evolutionary history, wherein on a global scale, spatially varying selective pressures have governed its natural history.

Patterns of host gene expression associated with harboring a foregut microbial community.

BACKGROUND: Harboring foregut microbial communities is considered a key innovation that allows herbivorous mammals to colonize new ecological niches. However, the functions of these chambers have only been well studied at the molecular level in ruminants. Here, we investigate gene expression in the foregut chamber of herbivorous rodents and ask whether these gene expression patterns are consistent with results in ruminants. We compared gene expression in foregut tissues of two rodent species: Stephen's woodrat (Neotoma stephensi), which harbors a dense foregut microbial community, and the lab rat (Rattus norvegicus), which lacks such a community. RESULTS: We found that woodrats have higher abundances of transcripts associated with smooth muscle processes, specifically a higher expression of the smoothelin-like 1 gene, which may assist in contractile properties of this tissue to retain food material in the foregut chamber. The expression of genes associated with keratinization and cornification exhibited a complex pattern of differences between the two species, suggesting distinct molecular mechanisms. Lab rats exhibited higher abundances of transcripts associated with immune function, likely to inhibit microbial growth in the foregut of this species. CONCLUSIONS: Some of our results were consistent with previous findings in ruminants (high expression of facilitative glucose transporters, lower expression of B4galnt2), suggestive of possible convergent evolution, while other results were unclear, and perhaps represent novel host-microbe interactions in rodents. Overall, our results suggest that harboring a foregut microbiota is associated with changes to the functions and host-microbe interactions of the foregut tissues.

Clinical phenotypes as predictors of the outcome of skipping around DMD exon 45.

OBJECTIVE: Exon-skipping therapies aim to convert Duchenne muscular dystrophy (DMD) into less severe Becker muscular dystrophy (BMD) by altering pre-mRNA splicing to restore an open reading frame, allowing translation of an internally deleted and partially functional dystrophin protein. The most common single exon deletion-exon 45 (Δ45)-may theoretically be treated by skipping of either flanking exon (44 or 46). We sought to predict the impact of these by assessing the clinical severity in dystrophinopathy patients. METHODS: Phenotypic data including clinical diagnosis, age at wheelchair use, age at loss of ambulation, and presence of cardiomyopathy were analyzed from 41 dystrophinopathy patients containing equivalent in-frame deletions. RESULTS: As expected, deletions of either exons 45 to 47 (Δ45-47) or exons 45 to 48 (Δ45-48) result in BMD in 97% (36 of 37) of subjects. Unexpectedly, deletion of exons 45 to 46 (Δ45-46) is associated with the more severe DMD phenotype in 4 of 4 subjects despite an in-frame transcript. Notably, no patients with a deletion of exons 44 to 45 (Δ44-45) were found within the United Dystrophinopathy Project database, and this mutation has only been reported twice before, which suggests an ascertainment bias attributable to a very mild phenotype. INTERPRETATION: The observation that Δ45-46 patients have typical DMD suggests that the conformation of the resultant protein may result in protein instability or altered binding of critical partners. We conclude that in DMD patients with Δ45, skipping of exon 44 and multiexon skipping of exons 46 and 47 (or exons 46-48) are better potential therapies than skipping of exon 46 alone.

A novel human-infection-derived bacterium provides insights into the evolutionary origins of mutualistic insect-bacterial symbioses.

Despite extensive study, little is known about the origins of the mutualistic bacterial endosymbionts that inhabit approximately 10% of the world's insects. In this study, we characterized a novel opportunistic human pathogen, designated "strain HS," and found that it is a close relative of the insect endosymbiont Sodalis glossinidius. Our results indicate that ancestral relatives of strain HS have served as progenitors for the independent descent of Sodalis-allied endosymbionts found in several insect hosts. Comparative analyses indicate that the gene inventories of the insect endosymbionts were independently derived from a common ancestral template through a combination of irreversible degenerative changes. Our results provide compelling support for the notion that mutualists evolve from pathogenic progenitors. They also elucidate the role of degenerative evolutionary processes in shaping the gene inventories of symbiotic bacteria at a very early stage in these mutualistic associations.

LTBP4 genotype predicts age of ambulatory loss in Duchenne muscular dystrophy.

OBJECTIVE: Duchenne muscular dystrophy (DMD) displays a clinical range that is not fully explained by the primary DMD mutations. Ltbp4, encoding latent transforming growth factor-ß binding protein 4, was previously discovered in a genome-wide scan as a modifier of murine muscular dystrophy. We sought to determine whether LTBP4 genotype influenced DMD severity in a large patient cohort. METHODS: We analyzed nonsynonymous single nucleotide polymorphisms (SNPs) from human LTBP4 in 254 nonambulatory subjects with known DMD mutations. These SNPs, V194I, T787A, T820A, and T1140M, form the VTTT and IAAM LTBP4 haplotypes. RESULTS: Individuals homozygous for the IAAM LTBP4 haplotype remained ambulatory significantly longer than those heterozygous or homozygous for the VTTT haplotype. Glucocorticoid-treated patients who were IAAM homozygotes lost ambulation at 12.5 ± 3.3 years compared to 10.7 ± 2.1 years for treated VTTT heterozygotes or homozygotes. IAAM fibroblasts exposed to transforming growth factor (TGF) ß displayed reduced phospho-SMAD signaling compared to VTTT fibroblasts, consistent with LTBP4' role as a regulator of TGFß. INTERPRETATION: LTBP4 haplotype influences age at loss of ambulation, and should be considered in the management of DMD patients.

Association of the CHRNA4 neuronal nicotinic receptor subunit gene with frequency of binge drinking in young adults.

BACKGROUND: Binge drinking is responsible for over half of all alcohol-related deaths and results in significant health and economic costs to individuals and society. Knowledge of genetic aspects of this behavior, particularly as it emerges in young adulthood, could lead to improved treatment and prevention programs. METHODS: We have focused on the association of variation in neuronal nicotinic receptor subunit genes (CHRNs) in a cohort of 702 Hispanic and non-Hispanic White young adults who are part of the Social and Emotional Contexts of Adolescent Smoking Patterns (SECASP) study. Fifty-five single nucleotide polymorphisms (SNPs) covering the variation in 5 CHRNs (CHRNA4, CHRNB2, CHRNA2, CHRNB3A6, and CHRNA5A3B4) were studied. RESULTS: Frequency of binge drinking and other correlated alcohol consumption measures were significantly associated with SNPs in CHRNA4 (p-values ranged from 0.0003 to 0.02), but not with SNPs in other CHRNs. This association was independent of smoking status in our cohort. CONCLUSIONS: Variants in CHRNA4 may contribute to risk of binge drinking in young adults in this cohort. Results will need to be confirmed in independent samples.

Effect of neuronal nicotinic acetylcholine receptor genes (CHRN) on longitudinal cigarettes per day in adolescents and young adults.

INTRODUCTION: Few studies have sought to identify specific genetic markers associated with cigarettes per day (CPD) during adolescence and young adulthood, the period of greatest vulnerability for the development of nicotine dependence. METHODS: We used a longitudinal design to investigate the effect of neuronal nicotinic acetylcholine receptor (CHRN) subunit genes on CPD from 15 to 21 years of age in young smokers of European descent (N = 439, 59% female). The number of CPD typically smoked during the previous 30 days was self-reported. Single nucleotide polymorphisms (SNPs) from CHRN genes were genotyped using DNA extracted from saliva samples collected at the 5-year assessment. Mixed-model analyses of SNP effects were computed across age at the time of assessment using log-transformed CPD as the phenotype. Data from the 1000 Genomes Project were used to clarify the architecture of CHRN genes to inform SNP selection and interpretation of results. RESULTS: CPD was associated with a CHRNB3A6 region tagged by rs2304297, with CHRNA5A3B4 haplotype C (tagged by rs569207), and with the CHRNA2 SNP rs2271920, ps < .004. The reliability of single-SNP associations was supported by the correspondence between a more extensive set of SNP signals and the underlying genetic architecture. The 3 signals identified in this study appear to make independent contributions to CPD, and their combined effect accounts for 5.5% of the variance in log-transformed CPD. CONCLUSIONS: Level of CPD during adolescence and young adulthood is associated with CHRNB3A6, CHRNA5A3B4, and CHRNA2.

Position of glycine substitutions in the triple helix of COL6A1, COL6A2, and COL6A3 is correlated with severity and mode of inheritance in collagen VI myopathies.

Glycine substitutions in the conserved Gly-X-Y motif in the triple helical (TH) domain of collagen VI are the most commonly identified mutations in the collagen VI myopathies including Ullrich congenital muscular dystrophy, Bethlem myopathy, and intermediate (INT) phenotypes. We describe clinical and genetic characteristics of 97 individuals with glycine substitutions in the TH domain of COL6A1, COL6A2, or COL6A3 and add a review of 97 published cases, for a total of 194 cases. Clinical findings include severe, INT, and mild phenotypes even from patients with identical mutations. INT phenotypes were most common, accounting for almost half of patients, emphasizing the importance of INT phenotypes to the overall phenotypic spectrum. Glycine substitutions in the TH domain are heavily clustered in a short segment N-terminal to the 17th Gly-X-Y triplet, where they are acting as dominants. The most severe cases are clustered in an even smaller region including Gly-X-Y triplets 10-15, accounting for only 5% of the TH domain. Our findings suggest that clustering of glycine substitutions in the N-terminal region of collagen VI is not based on features of the primary sequence. We hypothesize that this region may represent a functional domain within the triple helix.

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.

A genome-wide association analysis of loss of ambulation in dystrophinopathy patients suggests multiple candidate modifiers of disease severity.

The major determinant of disease severity in Duchenne muscular dystrophy (DMD) or milder Becker muscular dystrophy (BMD) is whether the dystrophin gene (DMD) mutation truncates the mRNA reading frame or allows expression of a partially functional protein. However, even in the complete absence of dystrophin, variability in disease severity is observed, and candidate gene studies have implicated several genes as modifiers. Here we present the largest genome-wide search to date for loci influencing severity in N = 419 DMD patients. Availability of subjects for such studies is quite limited, leading to modest sample sizes, which present a challenge for GWAS design. We have therefore taken special steps to minimize heterogeneity within our dataset at the DMD locus itself, taking a novel approach to mutation classification to effectively exclude the possibility of residual dystrophin expression, and utilized statistical methods that are well adapted to smaller sample sizes, including the use of a novel linear regression-like residual for time to ambulatory loss and the application of evidential statistics for the GWAS approach. Finally, we applied an unbiased in silico pipeline, utilizing functional genomic datasets to explore the potential impact of the best supported SNPs. In all, we obtained eight SNPs (out of 1,385,356 total) with posterior probability of trait-marker association (PPLD) ≥ 0.4, representing six distinct loci. Our analysis prioritized likely non-coding SNP regulatory effects on six genes (ETAA1, PARD6G, GALNTL6, MAN1A1, ADAMTS19, and NCALD), each with plausibility as a DMD modifier. These results support both recurrent and potentially new pathways for intervention in the dystrophinopathies.

Independent evolution of bitter-taste sensitivity in humans and chimpanzees.

It was reported over 65 years ago that chimpanzees, like humans, vary in taste sensitivity to the bitter compound phenylthiocarbamide (PTC). This was suggested to be the result of a shared balanced polymorphism, defining the first, and now classic, example of the effects of balancing selection in great apes. In humans, variable PTC sensitivity is largely controlled by the segregation of two common alleles at the TAS2R38 locus, which encode receptor variants with different ligand affinities. Here we show that PTC taste sensitivity in chimpanzees is also controlled by two common alleles of TAS2R38; however, neither of these alleles is shared with humans. Instead, a mutation of the initiation codon results in the use of an alternative downstream start codon and production of a truncated receptor variant that fails to respond to PTC in vitro. Association testing of PTC sensitivity in a cohort of captive chimpanzees confirmed that chimpanzee TAS2R38 genotype accurately predicts taster status in vivo. Therefore, although Fisher et al.'s observations were accurate, their explanation was wrong. Humans and chimpanzees share variable taste sensitivity to bitter compounds mediated by PTC receptor variants, but the molecular basis of this variation has arisen twice, independently, in the two species.