Whole-exome sequencing in syndromic craniosynostosis increases diagnostic yield and identifies candidate genes in osteogenic signaling pathways.

Craniosynostosis (CS) is a common congenital anomaly defined by premature fusion of one or more cranial sutures. Syndromic CS involves additional organ anomalies or neurocognitive deficits and accounts for 25%-30% of the cases. In a recent population-based study by our group, 84% of the syndromic CS cases had a genetically verified diagnosis after targeted analyses. A number of different genetic causes were detected, confirming that syndromic CS is highly heterogeneous. In this study, we performed whole-exome sequencing of 10 children and parents from the same cohort where previous genetic results were negative. We detected pathogenic, or likely pathogenic, variants in four additional genes (NFIA, EXTL3, POLR2A, and FOXP2) associated with rare conditions. In two of these (POLR2A and FOXP2), CS has not previously been reported. We further detected a rare predicted damaging variant in SH3BP4, which has not previously been related to human disease. All findings were clustered in genes involved in the pathways of osteogenesis and suture patency. We conclude that whole-exome sequencing expands the list of genes associated with syndromic CS, and provides new candidate genes in osteogenic signaling pathways.

Genetic risk variants for autoimmune diseases that influence gene expression in thymus.

Genome-wide association studies (GWAS) have boosted our knowledge of genetic risk variants in autoimmune diseases (AIDs). Most risk variants are located within or near genes with immunological functions, and the majority is found to be non-coding, pointing towards a regulatory role. In this study, we performed a cis expression quantitative trait locus (eQTL) screen restricted to 353 AID associated risk variants selected from the GWAS catalog to investigate whether these single nucleotide polymorphisms (SNPs) influence gene expression in thymus. Genotypes were obtained by Immunochip (Ichip) and tested against expression of surrounding genes (±1 Mb) in human thymic tissue (n = 42). We identified eight significant eQTLs located within seven genetic regions (FCRL3, RNASET2, C2orf74, NPIPB8, SIRPG, SYS1 and AJ006998.2) where the expression was associated with AID risk SNPs at a study-wide level of significance (P < 2.7 × 10-5). In NPIPB8 and AJ006998.2, the eQTL signals appeared to be thymus-specific. Furthermore, many AID risk SNPs from GWAS have been subsequently fine-mapped in recent Ichip projects, and fine-mapped AID SNPs overlapped with the thymic eQTLs within RNASET2 and SIRPG Finally, in all the eQTL regions, except C2orf74, SNPs underlying the thymic eQTLs were predicted to interfere with transcription factors important in T cell development. Our study therefore reveals autoimmune risk variants that act as eQTLs in thymus, and suggest that thymic gene regulation may play a functional role at some AID risk loci.

Biochemical and genetic characterization of an unusual mild PEX3-related Zellweger spectrum disorder.

Patients with PEX3 mutations usually present with a severe form of Zellweger spectrum disorder with death in the first year of life. Whole exome sequencing in adult siblings with intellectual disability revealed a homozygous variant in PEX3 that abolishes the normal splice site. A cryptic acceptor splice site is activated and an in-frame transcript with a deletion is produced. This transcript translates into a protein with residual activity explaining the relatively mild peroxisomal abnormalities and clinical phenotype.

Qualitative and quantitative analysis of FBN1 mRNA from 16 patients with Marfan Syndrome.

BACKGROUND: Pathogenic mutations in FBN1, encoding the glycoprotein, fibrillin-1, cause Marfan syndrome (MFS) and related connective tissue disorders. In the present study, qualitative and quantitative effects of 16 mutations, identified in FBN1 in MFS patients with systematically described phenotypes, were investigated in vitro. METHODS: Qualitative analysis was performed with reverse transcription-PCR (RT-PCR) and gel electrophoresis, and quantitative analysis to determine the FBN1 mRNA levels in fibroblasts from the 16 patients with MFS was performed with real-time PCR. RESULTS: Qualitative analysis documented that the mutations c.4817-2delA and c.A4925G led to aberrant FBN1 mRNA splicing leading to in frame deletion of exon 39 and in exon 39, respectively. No difference in the mean FBN1 mRNA level was observed between the entire group of cases and controls, nor between the group of patients with missense mutations and controls. The mean expression levels associated with premature termination codon (PTC) and splice site mutations were significantly lower than the levels in patients with missense mutations. A high level of FBN1 mRNA in the patient with the missense mutation c.G2447T did not segregate with the mutation in three of his first degree relatives. No association was indicated between the FBN1 transcript level and specific phenotypic manifestations. CONCLUSIONS: Abnormal FBN1 transcripts were indicated in fibroblasts from patients with the splice site mutation c.4817-2delA and the missense mutation c.A4925G. While the mean FBN1 mRNA expression level in fibroblasts from patients with splice site and PTC mutations were lower than the mean level in patients with missense mutations and controls, inter-individual variability was high. The observation that high level of FBN1 mRNA in the patient with the missense mutation c.G2447T did not segregate with the mutation in the family suggests that variable expression of the normal FBN1 allele may contribute to explain the variability in FBN1 mRNA level.

STIM1 R304W causes muscle degeneration and impaired platelet activation in mice.

STIM1 and ORAI1 regulate store-operated Ca2+ entry (SOCE) in most cell types, and mutations in these proteins have deleterious and diverse effects. We established a mouse line expressing the STIM1 R304 W gain-of-function mutation causing Stormorken syndrome to explore effects on organ and cell physiology. While STIM1 R304 W was lethal in the homozygous state, surviving mice presented with reduced growth, skeletal muscle degeneration, and reduced exercise endurance. Variable STIM1 expression levels between tissues directly impacted cellular SOCE capacity. In contrast to patients with Stormorken syndrome, STIM1 was downregulated in fibroblasts from Stim1R304W/R304W mice, which maintained SOCE despite constitutive protein activity. In studies using foetal liver chimeras, STIM1 protein was undetectable in homozygous megakaryocytes and platelets, resulting in impaired platelet activation and absent SOCE. These data indicate that downregulation of STIM1 R304 W effectively opposes the gain-of-function phenotype associated with this mutation, and highlight the importance of STIM1 in skeletal muscle development and integrity.

A comprehensive screen for SNP associations on chromosome region 5q31-33 in Swedish/Norwegian celiac disease families.

Celiac disease (CD) is a gluten-induced enteropathy, which results from the interplay between environmental and genetic factors. There is a strong human leukocyte antigen (HLA) association with the disease, and HLA-DQ alleles represent a major genetic risk factor. In addition to HLA-DQ, non-HLA genes appear to be crucial for CD development. Chromosomal region 5q31-33 has demonstrated linkage with CD in several genome-wide studies, including in our Swedish/Norwegian cohort. In a European meta-analysis 5q31-33 was the only region that reached a genome-wide level of significance except for the HLA region. To identify the genetic variant(s) responsible for this linkage signal, we performed a comprehensive single nucleotide polymorphism (SNP) association screen in 97 Swedish/Norwegian multiplex families who demonstrate linkage to the region. We selected tag SNPs from a 16 Mb region representing the 95% confidence interval of the linkage peak. A total of 1,404 SNPs were used for the association analysis. We identified several regions with SNPs demonstrating moderate single- or multipoint associations. However, the isolated association signals appeared insufficient to account for the linkage signal seen in our cohort. Collective effects of multiple risk genes within the region, incomplete genetic coverage or effects related to copy number variation are possible explanations for our findings.