Diversity of CFTR variants across ancestries characterized using 454,727 UK biobank whole exome sequences.

BACKGROUND: Limited understanding of the diversity of variants in the cystic fibrosis transmembrane conductance regulator (CFTR) gene across ancestries hampers efforts to advance molecular diagnosis of cystic fibrosis (CF). The consequences pose a risk of delayed diagnoses and subsequently worsened health outcomes for patients. Therefore, characterizing the spectrum of CFTR variants across ancestries is critical for revolutionizing molecular diagnoses of CF. METHODS: We analyzed 454,727 UK Biobank (UKBB) whole-exome sequences to characterize the diversity of CFTR variants across ancestries. Using the PanUKBB classification, the participants were assigned into six major groups: African (AFR), American/American Admixed (AMR), Central South Asia (CSA), East Asian (EAS), European (EUR), and Middle East (MID). We segregated ancestry-specific CFTR variants, including those that are CF-causing or clinically relevant. The ages of certain CF-causing variants were determined and analyzed for selective pressure effects, and curated phenotype analysis was performed for participants with clinically relevant CFTR genotypes. RESULTS: We detected over 4000 CFTR variants, including novel ancestry-specific variants, across six ancestries. Europeans had the most unique CFTR variants [n = 2212], while the American group had the least unique variants [n = 23]. F508del was the most prevalent CF-causing variant found in all ancestries, except in EAS, where V520F was the most prevalent. Common EAS variants such as 3600G > A, V456A, and V520, which appeared approximately 270, 215, and 338 generations ago, respectively, did not show evidence of selective pressure. Sixteen participants had two CF-causing variants, with two being diagnosed with CF. We found 154 participants harboring a CF-causing and varying clinical consequences (VCC) variant. Phenotype analysis performed for participants with multiple clinically relevant variants returned significant associations with CF and its pulmonary phenotypes [Bonferroni-adjusted p < 0.05]. CONCLUSIONS: We leveraged the UKBB database to comprehensively characterize the broad spectrum of CFTR variants across ancestries. The detection of over 4000 CFTR variants, including several ancestry-specific and uncharacterized CFTR variants, warrants the need for further characterization of their functional and clinical relevance. Overall, the presentation of classical CF phenotypes seen in non-CF diagnosed participants with more than one CF-causing variant indicates that they may benefit from current CFTR modulator therapies.

Functional characterization of a single nucleotide polymorphism associated with Alzheimer's disease in a hiPSC-based neuron model.

Neurodegenerative diseases encompass a group of debilitating conditions resulting from progressive nerve cell death. Of these, Alzheimer's disease (AD) occurs most frequently, but is currently incurable and has limited treatment success. Late onset AD, the most common form, is highly heritable but is caused by a combination of non-genetic risk factors and many low-effect genetic variants whose disease-causing mechanisms remain unclear. By mining the FinnGen study database of phenome-wide association studies, we identified a rare variant, rs148726219, enriched in the Finnish population that is associated with AD risk and dementia, and appears to have arisen on a common haplotype with older AD-associated variants such as rs429358. The rs148726219 variant lies in an overlapping intron of the FosB proto-oncogene (FOSB) and ERCC excision repair 1 (ERCC1) genes. To understand the impact of this SNP on disease phenotypes, we performed CRISPR/Cas9 editing in a human induced pluripotent stem cell (hiPSC) line to generate isogenic clones harboring heterozygous and homozygous alleles of rs148726219. hiPSC clones differentiated into induced excitatory neurons (iNs) did not exhibit detectable molecular or morphological variation in differentiation potential compared to isogenic controls. However, global transcriptome analysis showed differential regulation of nearby genes and upregulation of several biological pathways related to neuronal function, particularly synaptogenesis and calcium signaling, specifically in mature iNs harboring rs148726219 homozygous and heterozygous alleles. Functional differences in iN circuit maturation as measured by calcium imaging were observed across genotypes. Edited mature iNs also displayed downregulation of unfolded protein response and cell death pathways. This study implicates a phenotypic impact of rs148726219 in the context of mature neurons, consistent with its identification in late onset AD, and underscores a hiPSC-based experimental model to functionalize GWAS-identified variants.

Systematic single-variant and gene-based association testing of thousands of phenotypes in 394,841 UK Biobank exomes.

Genome-wide association studies have successfully discovered thousands of common variants associated with human diseases and traits, but the landscape of rare variations in human disease has not been explored at scale. Exome-sequencing studies of population biobanks provide an opportunity to systematically evaluate the impact of rare coding variations across a wide range of phenotypes to discover genes and allelic series relevant to human health and disease. Here, we present results from systematic association analyses of 4,529 phenotypes using single-variant and gene tests of 394,841 individuals in the UK Biobank with exome-sequence data. We find that the discovery of genetic associations is tightly linked to frequency and is correlated with metrics of deleteriousness and natural selection. We highlight biological findings elucidated by these data and release the dataset as a public resource alongside the Genebass browser for rapidly exploring rare-variant association results.

Concordance between gene expression in peripheral whole blood and colonic tissue in children with inflammatory bowel disease.

BACKGROUND: Presenting features of inflammatory bowel disease (IBD) are non-specific. We hypothesized that mRNA profiles could (1) identify genes and pathways involved in disease pathogenesis; (2) identify a molecular signature that differentiates IBD from other conditions; (3) provide insight into systemic and colon-specific dysregulation through study of the concordance of the gene expression. METHODS: Children (8-18 years) were prospectively recruited at the time of diagnostic colonoscopy for possible IBD. We used transcriptome-wide mRNA profiling to study gene expression in colon biopsies and paired whole blood samples. Using blood mRNA measurements, we fit a regression model for disease state prediction that was validated in an independent test set of adult subjects (GSE3365). RESULTS: Ninety-eight children were recruited [39 Crohn's disease, 18 ulcerative colitis, 2 IBDU, 39 non-IBD]. There were 1,118 significantly differentially (IBD vs non-IBD) expressed genes in colon tissue, and 880 in blood. The direction of relative change in expression was concordant for 106/112 genes differentially expressed in both tissue types. The regression model from the blood mRNA measurements distinguished IBD vs non-IBD disease status in the independent test set with 80% accuracy using only 6 genes. The overlap of 5 immune and metabolic pathways in the two tissue types was significant (p<0.001). CONCLUSIONS: Blood and colon tissue from patients with IBD share a common transcriptional profile dominated by immune and metabolic pathways. Our results suggest that peripheral blood expression levels of as few as 6 genes (IL7R, UBB, TXNIP, S100A8, ALAS2, and SLC2A3) may distinguish patients with IBD from non-IBD.

Transgenic zebrafish model of DUX4 misexpression reveals a developmental role in FSHD pathogenesis.

Facioscapulohumeral dystrophy type 1 (FSHD-1) is the most common autosomal dominant form of muscular dystrophy with a prevalence of ∼1 in 8000 individuals. It is considered a late-onset form of muscular dystrophy and leads to asymmetric muscle weakness in the facial, scapular, trunk and lower extremities. The prevalent hypothesis on disease pathogenesis is explained by misexpression of a germ line, primate-specific transcription factor DUX4-fl (double homeobox 4, full-length isoform) linked to the chromosome 4q35. In vitro and in vivo studies have demonstrated that very low levels of DUX4-fl expression are sufficient to induce an apoptotic and/or lethal phenotype, and therefore modeling of the disease has proved challenging. In this study, we expand upon our previously established injection model of DUX4 misexpression in zebrafish and describe a DUX4-inducible transgenic zebrafish model that better recapitulates the expression pattern and late onset phenotype characteristic of FSHD patients. We show that an induced burst of DUX4 expression during early development results in the onset of FSHD-like phenotypes in adulthood, even when DUX4 is no longer detectable. We also utilize our injection model to study long-term consequences of DUX4 expression in those that fail to show a developmental phenotype. Herein, we introduce a hypothesis that DUX4 expression during developmental stages is sufficient to induce FSHD-like phenotypes in later adulthood. Our findings point to a developmental role of DUX4 misexpression in the pathogenesis of FSHD and should be factored into the design of future therapies.

CD82 Is a Marker for Prospective Isolation of Human Muscle Satellite Cells and Is Linked to Muscular Dystrophies.

Cell-surface markers for prospective isolation of stem cells from human skeletal muscle have been difficult to identify. Such markers would be powerful tools for studying satellite cell function during homeostasis and in pathogenesis of diseases such as muscular dystrophies. In this study, we show that the tetraspanin KAI/CD82 is an excellent marker for prospectively isolating stem cells from human fetal and adult skeletal muscle. Human CD82+ muscle cells robustly engraft into a mouse model of muscular dystrophy. shRNA knockdown of CD82 in myogenic cells reduces myoblast proliferation, suggesting it is functionally involved in muscle homeostasis. CD82 physically interacts with alpha7beta1 integrin (α7ß1-ITG) and with α-sarcoglycan, a member of the Dystrophin-Associated Glycoprotein Complex (DAPC), both of which have been linked to muscular dystrophies. Consistently, CD82 expression is decreased in Duchenne muscular dystrophy patients. Together, these findings suggest that CD82 function may be important for muscle stem cell function in muscular disorders.

A missense mutation in TFRC, encoding transferrin receptor 1, causes combined immunodeficiency.

Patients with a combined immunodeficiency characterized by normal numbers but impaired function of T and B cells had a homozygous p.Tyr20His substitution in transferrin receptor 1 (TfR1), encoded by TFRC. The substitution disrupts the TfR1 internalization motif, resulting in defective receptor endocytosis and markedly increased TfR1 expression on the cell surface. Iron citrate rescued the lymphocyte defects, and expression of wild-type but not mutant TfR1 rescued impaired transferrin uptake in patient-derived fibroblasts. Tfrc(Y20H/Y20H) mice recapitulated the immunological defects of patients. Despite the critical role of TfR1 in erythrocyte development and function, patients had only mild anemia and only slightly increased TfR1 expression in erythroid precursors. We show that STEAP3, a metalloreductase expressed in erythroblasts, associates with TfR1 and partially rescues transferrin uptake in patient-derived fibroblasts, suggesting that STEAP3 may provide an accessory TfR1 endocytosis signal that spares patients from severe anemia. These findings demonstrate the importance of TfR1 in adaptive immunity.

Congenital sideroblastic anemia due to mutations in the mitochondrial HSP70 homologue HSPA9.

The congenital sideroblastic anemias (CSAs) are relatively uncommon diseases characterized by defects in mitochondrial heme synthesis, iron-sulfur (Fe-S) cluster biogenesis, or protein synthesis. Here we demonstrate that mutations in HSPA9, a mitochondrial HSP70 homolog located in the chromosome 5q deletion syndrome 5q33 critical deletion interval and involved in mitochondrial Fe-S biogenesis, result in CSA inherited as an autosomal recessive trait. In a fraction of patients with just 1 severe loss-of-function allele, expression of the clinical phenotype is associated with a common coding single nucleotide polymorphism in trans that correlates with reduced messenger RNA expression and results in a pseudodominant pattern of inheritance.

Emerging preclinical animal models for FSHD.

Facioscapulohumeral dystrophy (FSHD) is a unique and complex genetic disease that is not entirely solved. Recent advances in the field have led to a consensus genetic premise for the disorder, enabling researchers to now pursue the design of preclinical models. In this review we explore all available FSHD models (DUX4-dependent and -independent) for their utility in therapeutic discovery and potential to yield novel disease insights. Owing to the complex nature of FSHD, there is currently no single model that accurately recapitulates the genetic and pathophysiological spectrum of the disorder. Existing models emphasize only specific aspects of the disease, highlighting the need for more collaborative research and novel paradigms to advance the translational research space of FSHD.

MicroRNA-486-dependent modulation of DOCK3/PTEN/AKT signaling pathways improves muscular dystrophy-associated symptoms.

Duchenne muscular dystrophy (DMD) is caused by mutations in the gene encoding dystrophin, which results in dysfunctional signaling pathways within muscle. Previously, we identified microRNA-486 (miR-486) as a muscle-enriched microRNA that is markedly reduced in the muscles of dystrophin-deficient mice (Dmdmdx-5Cv mice) and in DMD patient muscles. Here, we determined that muscle-specific transgenic overexpression of miR-486 in muscle of Dmdmdx-5Cv mice results in reduced serum creatine kinase levels, improved sarcolemmal integrity, fewer centralized myonuclei, increased myofiber size, and improved muscle physiology and performance. Additionally, we identified dedicator of cytokinesis 3 (DOCK3) as a miR-486 target in skeletal muscle and determined that DOCK3 expression is induced in dystrophic muscles. DOCK3 overexpression in human myotubes modulated PTEN/AKT signaling, which regulates muscle hypertrophy and growth, and induced apoptosis. Furthermore, several components of the PTEN/AKT pathway were markedly modulated by miR-486 in dystrophin-deficient muscle. Skeletal muscle-specific miR-486 overexpression in Dmdmdx-5Cv animals decreased levels of DOCK3, reduced PTEN expression, and subsequently increased levels of phosphorylated AKT, which resulted in an overall beneficial effect. Together, these studies demonstrate that stable overexpression of miR-486 ameliorates the disease progression of dystrophin-deficient skeletal muscle.

Human skeletal muscle xenograft as a new preclinical model for muscle disorders.

Development of novel therapeutics requires good animal models of disease. Disorders for which good animal models do not exist have very few drugs in development or clinical trial. Even where there are accepted, albeit imperfect models, the leap from promising preclinical drug results to positive clinical trials commonly fails, including in disorders of skeletal muscle. The main alternative model for early drug development, tissue culture, lacks both the architecture and, usually, the metabolic fidelity of the normal tissue in vivo. Herein, we demonstrate the feasibility and validity of human to mouse xenografts as a preclinical model of myopathy. Human skeletal muscle biopsies transplanted into the anterior tibial compartment of the hindlimbs of NOD-Rag1(null) IL2rγ(null) immunodeficient host mice regenerate new vascularized and innervated myofibers from human myogenic precursor cells. The grafts exhibit contractile and calcium release behavior, characteristic of functional muscle tissue. The validity of the human graft as a model of facioscapulohumeral muscular dystrophy is demonstrated in disease biomarker studies, showing that gene expression profiles of xenografts mirror those of the fresh donor biopsies. These findings illustrate the value of a new experimental model of muscle disease, the human muscle xenograft in mice, as a feasible and valid preclinical tool to better investigate the pathogenesis of human genetic myopathies and to more accurately predict their response to novel therapeutics.

An etiologic regulatory mutation in IRF6 with loss- and gain-of-function effects.

DNA variation in Interferon Regulatory Factor 6 (IRF6) causes Van der Woude syndrome (VWS), the most common syndromic form of cleft lip and palate (CLP). However, an etiologic variant in IRF6 has been found in only 70% of VWS families. To test whether DNA variants in regulatory elements cause VWS, we sequenced three conserved elements near IRF6 in 70 VWS families that lack an etiologic mutation within IRF6 exons. A rare mutation (350dupA) was found in a conserved IRF6 enhancer element (MCS9.7) in a Brazilian family. The 350dupA mutation abrogated the binding of p63 and E47 transcription factors to cis-overlapping motifs, and significantly disrupted enhancer activity in human cell cultures. Moreover, using a transgenic assay in mice, the 350dupA mutation disrupted the activation of MCS9.7 enhancer element and led to failure of lacZ expression in all head and neck pharyngeal arches. Interestingly, disruption of the p63 Motif1 and/or E47 binding sites by nucleotide substitution did not fully recapitulate the effect of the 350dupA mutation. Rather, we recognized that the 350dupA created a CAAAGT motif, a binding site for Lef1 protein. We showed that Lef1 binds to the mutated site and that overexpression of Lef1/ß-Catenin chimeric protein repressed MCS9.7-350dupA enhancer activity. In conclusion, our data strongly suggest that 350dupA variant is an etiologic mutation in VWS patients and disrupts enhancer activity by a loss- and gain-of-function mechanism, and thus support the rationale for additional screening for regulatory mutations in patients with CLP.

Exome sequencing identifies a novel SMCHD1 mutation in facioscapulohumeral muscular dystrophy 2.

FSHD2 is a rare form of facioscapulohumeral muscular dystrophy (FSHD) characterized by the absence of a contraction in the D4Z4 macrosatellite repeat region on chromosome 4q35 that is the hallmark of FSHD1. However, hypomethylation of this region is common to both subtypes. Recently, mutations in SMCHD1 combined with a permissive 4q35 allele were reported to cause FSHD2. We identified a novel p.Lys275del SMCHD1 mutation in a family affected with FSHD2 using whole-exome sequencing and linkage analysis. This mutation alters a highly conserved amino acid in the ATPase domain of SMCHD1. Subject III-11 is a male who developed asymmetrical muscle weakness characteristic of FSHD at 13 years. Physical examination revealed marked bilateral atrophy at biceps brachii, bilateral scapular winging, some asymmetrical weakness at tibialis anterior and peroneal muscles, and mild lower facial weakness. Biopsy of biceps brachii in subject II-5, the father of III-11, demonstrated lobulated fibers and dystrophic changes. Endomysial and perivascular inflammation was found, which has been reported in FSHD1 but not FSHD2. Given the previous report of SMCHD1 mutations in FSHD2 and the clinical presentations consistent with the FSHD phenotype, we conclude that the SMCHD1 mutation is the likely cause of the disease in this family.

The cell biology of disease: cellular and molecular mechanisms underlying muscular dystrophy.

The muscular dystrophies are a group of heterogeneous genetic diseases characterized by progressive degeneration and weakness of skeletal muscle. Since the discovery of the first muscular dystrophy gene encoding dystrophin, a large number of genes have been identified that are involved in various muscle-wasting and neuromuscular disorders. Human genetic studies complemented by animal model systems have substantially contributed to our understanding of the molecular pathomechanisms underlying muscle degeneration. Moreover, these studies have revealed distinct molecular and cellular mechanisms that link genetic mutations to diverse muscle wasting phenotypes.

Telomere position effect regulates DUX4 in human facioscapulohumeral muscular dystrophy.

Telomeres may regulate human disease by at least two independent mechanisms. First, replicative senescence occurs once short telomeres generate DNA-damage signals that produce a barrier to tumor progression. Second, telomere position effects (TPE) could change gene expression at intermediate telomere lengths in cultured human cells. Here we report that telomere length may contribute to the pathogenesis of facioscapulohumeral muscular dystrophy (FSHD). FSHD is a late-onset disease genetically residing only 25-60 kilobases from the end of chromosome 4q. We used a floxable telomerase to generate isogenic clones with different telomere lengths from affected patients and their unaffected siblings. DUX4, the primary candidate for FSHD pathogenesis, is upregulated over ten-fold in FSHD myoblasts and myotubes with short telomeres, and its expression is inversely proportional to telomere length. FSHD may be the first known human disease in which TPE contributes to age-related phenotype.

Transcriptional profiling in facioscapulohumeral muscular dystrophy to identify candidate biomarkers.

Facioscapulohumeral muscular dystrophy (FSHD) is a progressive neuromuscular disorder caused by contractions of repetitive elements within the macrosatellite D4Z4 on chromosome 4q35. The pathophysiology of FSHD is unknown and, as a result, there is currently no effective treatment available for this disease. To better understand the pathophysiology of FSHD and develop mRNA-based biomarkers of affected muscles, we compared global analysis of gene expression in two distinct muscles obtained from a large number of FSHD subjects and their unaffected first-degree relatives. Gene expression in two muscle types was analyzed using GeneChip Gene 1.0 ST arrays: biceps, which typically shows an early and severe disease involvement; and deltoid, which is relatively uninvolved. For both muscle types, the expression differences were mild: using relaxed cutoffs for differential expression (fold change ≥1.2; nominal P value <0.01), we identified 191 and 110 genes differentially expressed between affected and control samples of biceps and deltoid muscle tissues, respectively, with 29 genes in common. Controlling for a false-discovery rate of <0.25 reduced the number of differentially expressed genes in biceps to 188 and in deltoid to 7. Expression levels of 15 genes altered in this study were used as a "molecular signature" in a validation study of an additional 26 subjects and predicted them as FSHD or control with 90% accuracy based on biceps and 80% accuracy based on deltoids.

Facioscapulohumeral muscular dystrophy family studies of DUX4 expression: evidence for disease modifiers and a quantitative model of pathogenesis.

Facioscapulohumeral muscular dystrophy (FSHD), the most prevalent myopathy afflicting both children and adults, is predominantly associated with contractions in the 4q35-localized macrosatellite D4Z4 repeat array. Recent studies have proposed that FSHD pathology is caused by the misexpression of the DUX4 (double homeobox 4) gene resulting in production of a pathogenic protein, DUX4-FL, which has been detected in FSHD, but not in unaffected control myogenic cells and muscle tissue. Here, we report the analysis of DUX4 mRNA and protein expression in a much larger collection of myogenic cells and muscle biopsies derived from biceps and deltoid muscles of FSHD affected subjects and their unaffected first-degree relatives. We confirmed that stable DUX4-fl mRNA and protein were expressed in myogenic cells and muscle tissues derived from FSHD affected subjects, including several genetically diagnosed adult FSHD subjects yet to show clinical manifestations of the disease in the assayed muscles. In addition, we report DUX4-fl mRNA and protein expression in muscle biopsies and myogenic cells from genetically unaffected relatives of the FSHD subjects, although at a significantly lower frequency. These results establish that DUX4-fl expression per se is not sufficient for FSHD muscle pathology and indicate that quantitative modifiers of DUX4-fl expression and/or function and family genetic background are determinants of FSHD muscle disease progression.

A unique library of myogenic cells from facioscapulohumeral muscular dystrophy subjects and unaffected relatives: family, disease and cell function.

To explore possible mechanisms of pathology in facioscapulohumeral muscular dystrophy (FSHD), we generated a novel library of myogenic cells composed of paired cultures derived from FSHD subjects and unaffected first-degree relatives. We prepared cells from biopsies of both biceps and deltoid muscles obtained from each of 10 FSHD and 9 unaffected donors. We used this new collection to determine how family background and disease affected patterns of growth and differentiation, expression of a panel of candidate, and muscle-specific genes, and responses to exogenous stressors. We found that FSHD and unaffected cells had, on average, indistinguishable patterns of differentiation, gene expression, and dose-response curves to staurosporine, paraquat, hydrogen peroxide, and glutathione depletion. Differentiated FSHD and unaffected cultures were both more sensitive to glutathione depletion than proliferating cultures, but showed similar responses to paraquat, staurosporine, and peroxide. For stress responses, the sample size was sufficient to detect a 10% change in effect at the observed variability with a power of >99%. In contrast, for each of these properties, we found significant differences among cells from different cohorts, and these differences were independent of disease status, gender, or muscle biopsied. Thus, though none of the properties we examined could be used to reliably distinguish between FSHD and unaffected cells, family of origin was an important contributor to gene-expression patterns and stressor responses in cultures of both FSHD and unaffected myogenic cells.

Genetics of nonsyndromic orofacial clefts.

With an average worldwide prevalence of approximately 1.2/1000 live births, orofacial clefts are the most common craniofacial birth defects in humans. Like other complex disorders, these birth defects are thought to result from the complex interplay of multiple genes and environmental factors. Significant progress in the identification of underlying genes and pathways has benefited from large populations available for study, increased international collaboration, rapid advances in genotyping technology, and major improvements in analytic approaches. Here we review recent advances in genetic epidemiological approaches to complex traits and their applications to studies of nonsyndromic orofacial clefts. Our main aim is to bring together a discussion of new and previously identified candidate genes to create a more cohesive picture of interacting pathways that shape the human craniofacial region. In future directions, we highlight the need to search for copy number variants that affect gene dosage and rare variants that are possibly associated with a higher disease penetrance. In addition, sequencing of protein-coding regions in candidate genes and screening for genetic variation in noncoding regulatory elements will help advance this important area of research.