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.

FinnGen provides genetic insights from a well-phenotyped isolated population

Population isolates such as those in Finland benefit genetic research because deleterious alleles are often concentrated on a small number of low-frequency variants (0.1% ≤ minor allele frequency < 5%). These variants survived the founding bottleneck rather than being distributed over a large number of ultrarare variants. Although this effect is well established in Mendelian genetics, its value in common disease genetics is less explored^1,2. FinnGen aims to study the genome and national health register data of 500,000 Finnish individuals. Given the relatively high median age of participants (63 years) and the substantial fraction of hospital-based recruitment, FinnGen is enriched for disease end points. Here we analyse data from 224,737 participants from FinnGen and study 15 diseases that have previously been investigated in large genome-wide association studies (GWASs). We also include meta-analyses of biobank data from Estonia and the United Kingdom. We identified 30 new associations, primarily low-frequency variants, enriched in the Finnish population. A GWAS of 1,932 diseases also identified 2,733 genome-wide significant associations (893 phenome-wide significant (PWS), P < 2.6 × 10^-11) at 2,496 (771 PWS) independent loci with 807 (247 PWS) end points. Among these, fine-mapping implicated 148 (73 PWS) coding variants associated with 83 (42 PWS) end points. Moreover, 91 (47 PWS) had an allele frequency of <5% in non-Finnish European individuals, of which 62 (32 PWS) were enriched by more than twofold in Finland. These findings demonstrate the power of bottlenecked populations to find entry points into the biology of common diseases through low-frequency, high impact 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.

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.

Comparative RNA editing in autistic and neurotypical cerebella

Adenosine-to-inosine (A-to-I) RNA editing is a neurodevelopmentally regulated epigenetic modification shown to modulate complex behavior in animals. Little is known about human A-to-I editing, but it is thought to constitute one of many molecular mechanisms connecting environmental stimuli and behavioral outputs. Thus, comprehensive exploration of A-to-I RNA editing in human brains may shed light on gene-environment interactions underlying complex behavior in health and disease. Synaptic function is a main target of A-to-I editing, which can selectively recode key amino acids in synaptic genes, directly altering synaptic strength and duration in response to environmental signals. Here, we performed a high-resolution survey of synaptic A-to-I RNA editing in a human population, and examined how it varies in autism, a neurodevelopmental disorder in which synaptic abnormalities are a common finding. Using ultra-deep (>1000 × ) sequencing, we quantified the levels of A-to-I editing of 10 synaptic genes in postmortem cerebella from 14 neurotypical and 11 autistic individuals. A high dynamic range of editing levels was detected across individuals and editing sites, from 99.6% to below detection limits. In most sites, the extreme ends of the population editing distributions were individuals with autism. Editing was correlated with isoform usage, clusters of correlated sites were identified, and differential editing patterns examined. Finally, a dysfunctional form of the editing enzyme adenosine deaminase acting on RNA B1 was found more commonly in postmortem cerebella from individuals with autism. These results provide a population-level, high-resolution view of A-to-I RNA editing in human cerebella and suggest that A-to-I editing of synaptic genes may be informative for assessing the epigenetic risk for autism.

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.

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.

Gene expression profiling of skeletal muscles treated with a soluble activin type IIB receptor

Inhibition of the myostatin signaling pathway is emerging as a promising therapeutic means to treat muscle wasting and degenerative disorders. Activin type IIB receptor (ActRIIB) is the putative myostatin receptor, and a soluble activin receptor (ActRIIB-Fc) has been demonstrated to potently inhibit a subset of transforming growth factor (TGF)-β family members including myostatin. To determine reliable and valid biomarkers for ActRIIB-Fc treatment, we assessed gene expression profiles for quadriceps muscles from mice treated with ActRIIB-Fc compared with mice genetically lacking myostatin and control mice. Expression of 134 genes was significantly altered in mice treated with ActRIIB-Fc over a 2-wk period relative to control mice (fold change > 1.5, P < 0.001), whereas the number of significantly altered genes in mice treated for 2 days was 38, demonstrating a time-dependent response to ActRIIB-Fc in overall muscle gene expression. The number of significantly altered genes in Mstn(-/-) mice relative to control mice was substantially higher (360), but for most of these genes the expression levels in the 2-wk treated mice were closer to the levels in the Mstn(-/-) mice than in control mice (P < 10⁻³⁰). Expression levels of 30 selected genes were further validated with quantitative real-time polymerase chain reaction (qPCR), and a correlation of ≥ 0.89 was observed between the fold changes from the microarray analysis and the qPCR analysis. These data suggest that treatment with ActRIIB-Fc results in overlapping but distinct gene expression signatures compared with myostatin genetic mutation. Differentially expressed genes identified in this study can be used as potential biomarkers for ActRIIB-Fc treatment, which is currently in clinical trials as a therapeutic agent for muscle wasting and degenerative disorders.

Disruption of an AP-2alpha binding site in an IRF6 enhancer is associated with cleft lip

Previously we have shown that nonsyndromic cleft lip with or without cleft palate (NSCL/P) is strongly associated with SNPs in IRF6 (interferon regulatory factor 6). Here, we use multispecies sequence comparisons to identify a common SNP (rs642961, G>A) in a newly identified IRF6 enhancer. The A allele is significantly overtransmitted (P = 1 x 10(-11)) in families with NSCL/P, in particular those with cleft lip but not cleft palate. Further, there is a dosage effect of the A allele, with a relative risk for cleft lip of 1.68 for the AG genotype and 2.40 for the AA genotype. EMSA and ChIP assays demonstrate that the risk allele disrupts the binding site of transcription factor AP-2alpha and expression analysis in the mouse localizes the enhancer activity to craniofacial and limb structures. Our findings place IRF6 and AP-2alpha in the same developmental pathway and identify a high-frequency variant in a regulatory element contributing substantially to a common, complex disorder.

Genetic variants in IRF6 and the risk of facial clefts: single-marker and haplotype-based analyses in a population-based case-control study of facial clefts in Norway

Mutations in the gene encoding interferon regulatory factor 6 (IRF6) underlie a common form of syndromic clefting known as Van der Woude syndrome. Lip pits and missing teeth are the only additional features distinguishing the syndrome from isolated clefts. Van der Woude syndrome, therefore, provides an excellent model for studying the isolated forms of clefting. From a population-based case-control study of facial clefts in Norway (1996-2001), we selected 377 cleft lip with or without cleft palate (CL/P), 196 cleft palate only (CPO), and 763 control infant-parent triads for analysis. We genotyped six single nucleotide polymorphisms within the IRF6 locus and estimated the relative risks (RR) conferred on the child by alleles and haplotypes of the child and of the mother. On the whole, there were strong statistical associations with CL/P but not CPO in our data. In single-marker analyses, mothers with a double-dose of the 'a'-allele at rs4844880 had an increased risk of having a child with CL/P (RR=1.85, 95% confidence interval: 1.04-3.25; P=0.036). An RR of 0.38 (95% confidence interval: 0.16-0.92; P=0.031) was obtained when the child carried a single-dose of the 'a'-allele at rs2235371 (the p.V274I polymorphism). The P-value for the overall test was <0.001. In haplotype analyses, several of the fetal and maternal haplotype relative risks were statistically significant individually but were not strong enough to show up on the overall test (P=0.113). Taken together, these findings further support a role for IRF6 variants in clefting of the lip and provide specific risk estimates in a Norwegian population.

GLI2 mutations in four Brazilian patients: how wide is the phenotypic spectrum?

We report four patients with GLI2 mutations together with their associated phenotypes: (1) holoprosencephaly-like phenotype, (2) anophthalmia, branchial arch anomalies, and CNS abnormalities, (3) heminasal aplasia and orbital anomalies, and (4) lobar holoprosencephaly. This diversity of phenotypes expands our understanding. Findings include not only (1) holoprosencephaly or a holoprosencephaly-like phenotype, but also (2) heminasal aplasia with orbital anomalies, and (3) branchial arch anomalies of the type seen in hemifacial microsomia with anophthalmia and in oculoauriculofrontonasal syndrome. Finally, this is the first report of a double mutation involving GLI2 and PTCH in the same patient.

Medical sequencing of candidate genes for nonsyndromic cleft lip and palate

Nonsyndromic or isolated cleft lip with or without cleft palate (CL/P) occurs in wide geographic distribution with an average birth prevalence of 1/700. We used direct sequencing as an approach to study candidate genes for CL/P. We report here the results of sequencing on 20 candidate genes for clefts in 184 cases with CL/P selected with an emphasis on severity and positive family history. Genes were selected based on expression patterns, animal models, and/or role in known human clefting syndromes. For seven genes with identified coding mutations that are potentially etiologic, we performed linkage disequilibrium studies as well in 501 family triads (affected child/mother/father). The recently reported MSX1 P147Q mutation was also studied in an additional 1,098 cleft cases. Selected missense mutations were screened in 1,064 controls from unrelated individuals on the Centre d'Étude du Polymorphisme Humain (CEPH) diversity cell line panel. Our aggregate data suggest that point mutations in these candidate genes are likely to contribute to 6% of isolated clefts, particularly those with more severe phenotypes (bilateral cleft of the lip with cleft palate). Additional cases, possibly due to microdeletions or isodisomy, were also detected and may contribute to clefts as well. Sequence analysis alone suggests that point mutations in FOXE1, GLI2, JAG2, LHX8, MSX1, MSX2, SATB2, SKI, SPRY2, and TBX10 may be rare causes of isolated cleft lip with or without cleft palate, and the linkage disequilibrium data support a larger, as yet unspecified, role for variants in or near MSX2, JAG2, and SKI. This study also illustrates the need to test large numbers of controls to distinguish rare polymorphic variants and prioritize functional studies for rare point mutations.

Nonsyndromic or isolated cleft lip with or without cleft palate (CL/P) is a birth defect with wide geographic distribution, occurring with an average frequency of 1/700 live births. Treatment can be provided, but it will involve medical, surgical, dental, and psychological personnel. Several different genes have been implicated in different cases. Here the researchers report the results of sequencing 20 different genes in 184 CL/P cases selected with an emphasis on more severe cases and cases with a positive family history for CL/P. Genes were selected based on previous work done by others and by the researchers' group. The authors' results suggest that point mutations in these candidate genes are likely to contribute to about 5% of CL/P, and particularly those with more severe phenotypes (bilateral cleft of the lip with cleft palate). This study also illustrates the need to test large numbers of controls to distinguish rare polymorphic variants and allow investigators to focus functional studies on the rare point mutations that seem to be disease-causing, so that researchers might better understand the mechanisms that play a role in CL/P.

Variation in the HLA-G promoter region influences miscarriage rates

The HLA-G gene is primarily expressed in placental cells that invade the maternal decidua during pregnancy. This gene encodes multiple isoforms that fulfill a variety of functions at the maternal-fetal interface throughout gestation. Recently, a null allele for the most abundant HLA-G isoform was associated with recurrent miscarriage in two independent studies, suggesting that reduced levels of the HLA-G1 protein may compromise successful pregnancy. We initiated the present study to determine whether other polymorphisms that could affect expression levels of HLA-G were associated with fetal loss in women participating in a 15-year prospective study of pregnancy outcome. We genotyped these subjects for 18 single-nucleotide polymorphisms in the 1,300 bp upstream of exon 1, 13 of which were identified as part of this study, as well as for an insertion/deletion (in/del) polymorphism in the 3' untranslated region. The 18 SNPs defined eight unique haplotypes. One polymorphism, -725C/G, was associated with fetal loss, with an increased risk for miscarriage in couples in which both partners carried the -725G allele, compared with couples not carrying this allele (odds ratio 2.76, 95% confidence interval 1.08-7.09; P=.035). Further, the G at nucleotide -725 creates a CpG dinucleotide, and we demonstrate that this CpG site is methylated on -725G alleles. Overall, this study identified extraordinary levels of variation in the 5'-upstream regulatory region of HLA-G and provides evidence for an association between a promoter-region SNP and fetal loss rates, further attesting to the novel features and critical role of this gene in pregnancy.