Genome-wide meta-analysis identifies BARX1 and EML4-MTA3 as new loci associated with infantile hypertrophic pyloric stenosis.

Infantile hypertrophic pyloric stenosis (IHPS) is a disorder of young infants with a population incidence of ∼2/1000 live births, caused by hypertrophy of the pyloric sphincter smooth muscle. Reported genetic loci associated with IHPS explain only a minor proportion of IHPS risk. To identify new risk loci, we carried out a genome-wide meta-analysis on 1395 surgery-confirmed cases and 4438 controls, with replication in a set of 2427 cases and 2524 controls. We identified and replicated six independent genomic loci associated with IHPS risk at genome wide significance (P < 5 × 10-8), including novel associations with two single nucleotide polymorphisms (SNPs). One of these SNPs, rs6736913 [odds ratio (OR) = 2.32; P = 3.0 × 10-15], is a low frequency missense variant in EML4 at 2p21. The second SNP, rs1933683 (OR = 1.34; P = 3.1 × 10-9) is 1 kb downstream of BARX1 at 9q22.32, an essential gene for stomach formation in embryogenesis. Using the genome-wide complex trait analysis method, we estimated the IHPS SNP heritability to be 30%, and using the linkage disequilibrium score regression method, we found support for a previously reported genetic correlation of IHPS with lipid metabolism. By combining the largest collection of IHPS cases to date (3822 cases), with results generalized across populations of different ancestry, we elucidate novel mechanistic avenues of IHPS disease architecture.

Rare variants in dynein heavy chain genes in two individuals with situs inversus and developmental dyslexia: a case report.

BACKGROUND: Developmental dyslexia (DD) is a neurodevelopmental learning disorder with high heritability. A number of candidate susceptibility genes have been identified, some of which are linked to the function of the cilium, an organelle regulating left-right asymmetry development in the embryo. Furthermore, it has been suggested that disrupted left-right asymmetry of the brain may play a role in neurodevelopmental disorders such as DD. However, it is unknown whether there is a common genetic cause to DD and laterality defects or ciliopathies. CASE PRESENTATION: Here, we studied two individuals with co-occurring situs inversus (SI) and DD using whole genome sequencing to identify genetic variants of importance for DD and SI. Individual 1 had primary ciliary dyskinesia (PCD), a rare, autosomal recessive disorder with oto-sino-pulmonary phenotype and SI. We identified two rare nonsynonymous variants in the dynein axonemal heavy chain 5 gene (DNAH5): a previously reported variant c.7502G > C; p.(R2501P), and a novel variant c.12043 T > G; p.(Y4015D). Both variants are predicted to be damaging. Ultrastructural analysis of the cilia revealed a lack of outer dynein arms and normal inner dynein arms. MRI of the brain revealed no significant abnormalities. Individual 2 had non-syndromic SI and DD. In individual 2, one rare variant (c.9110A > G;p.(H3037R)) in the dynein axonemal heavy chain 11 gene (DNAH11), coding for another component of the outer dynein arm, was identified. CONCLUSIONS: We identified the likely genetic cause of SI and PCD in one individual, and a possibly significant heterozygosity in the other, both involving dynein genes. Given the present evidence, it is unclear if the identified variants also predispose to DD and further studies into the association between laterality, ciliopathies and DD are needed.

CTNND2-a candidate gene for reading problems and mild intellectual disability.

BACKGROUND: Cytogenetically visible chromosomal translocations are highly informative as they can pinpoint strong effect genes even in complex genetic disorders. METHODS AND RESULTS: Here, we report a mother and daughter, both with borderline intelligence and learning problems within the dyslexia spectrum, and two apparently balanced reciprocal translocations: t(1;8)(p22;q24) and t(5;18)(p15;q11). By low coverage mate-pair whole-genome sequencing, we were able to pinpoint the genomic breakpoints to 2 kb intervals. By direct sequencing, we then located the chromosome 5p breakpoint to intron 9 of CTNND2. An additional case with a 163 kb microdeletion exclusively involving CTNND2 was identified with genome-wide array comparative genomic hybridisation. This microdeletion at 5p15.2 is also present in mosaic state in the patient's mother but absent from the healthy siblings. We then investigated the effect of CTNND2 polymorphisms on normal variability and identified a polymorphism (rs2561622) with significant effect on phonological ability and white matter volume in the left frontal lobe, close to cortical regions previously associated with phonological processing. Finally, given the potential role of CTNND2 in neuron motility, we used morpholino knockdown in zebrafish embryos to assess its effects on neuronal migration in vivo. Analysis of the zebrafish forebrain revealed a subpopulation of neurons misplaced between the diencephalon and telencephalon. CONCLUSIONS: Taken together, our human genetic and in vivo data suggest that defective migration of subpopulations of neuronal cells due to haploinsufficiency of CTNND2 contribute to the cognitive dysfunction in our patients.

Targeted high-throughput sequencing of candidate genes for chronic obstructive pulmonary disease.

BACKGROUND: Reduced lung function in patients with chronic obstructive pulmonary disease (COPD) is likely due to both environmental and genetic factors. We report here a targeted high-throughput DNA sequencing approach to identify new and previously known genetic variants in a set of candidate genes for COPD. METHODS: Exons in 22 genes implicated in lung development as well as 61 genes and 10 genomic regions previously associated with COPD were sequenced using individual DNA samples from 68 cases with moderate or severe COPD and 66 controls matched for age, gender and smoking. Cases and controls were selected from the Obstructive Lung Disease in Northern Sweden (OLIN) studies. RESULTS: In total, 37 genetic variants showed association with COPD (p < 0.05, uncorrected). Several variants previously discovered to be associated with COPD from genetic genome-wide analysis studies were replicated using our sample. Two high-risk variants were followed-up for functional characterization in a large eQTL mapping study of 1,111 human lung specimens. The C allele of a synonymous variant, rs8040868, predicting a p.(S45=) in the gene for cholinergic receptor nicotinic alpha 3 (CHRNA3) was associated with COPD (p = 8.8 x 10-3). This association remained (p = 0.003 and OR = 1.4, 95 % CI 1.1-1.7) when analysing all available cases and controls in OLIN (n = 1,534). The rs8040868 variant is in linkage disequilibrium with rs16969968 previously associated with COPD and altered expression of the CHRNA5 gene. A follow-up analysis for detection of expression quantitative trait loci revealed that rs8040868-C was found to be significantly associated with a decreased expression of the nearby gene cholinergic receptor, nicotinic, alpha 5 (CHRNA5) in lung tissue. CONCLUSION: Our data replicate previous result suggesting CHRNA5 as a candidate gene for COPD and rs8040868 as a risk variant for the development of COPD in the Swedish population.

DCDC2 polymorphism is associated with left temporoparietal gray and white matter structures during development.

Three genes, DYX1C1, DCDC2, and KIAA0319, have been previously associated with dyslexia, neuronal migration, and ciliary function. Three polymorphisms within these genes, rs3743204 (DYX1C1), rs793842 (DCDC2), and rs6935076 (KIAA0319) have also been linked to normal variability of left temporoparietal white matter volume connecting the middle temporal cortex to the angular and supramarginal gyri. Here, we assessed whether these polymorphisms are also related to the cortical thickness of the associated regions during childhood development using a longitudinal dataset of 76 randomly selected children and young adults who were scanned up to three times each, 2 years apart. rs793842 in DCDC2 was significantly associated with the thickness of left angular and supramarginal gyri as well as the left lateral occipital cortex. The cortex was significantly thicker for T-allele carriers, who also had lower white matter volume and lower reading comprehension scores. There was a negative correlation between white matter volume and cortical thickness, but only white matter volume predicted reading comprehension 2 years after scanning. These results show how normal variability in reading comprehension is related to gene, white matter volume, and cortical thickness in the inferior parietal lobe. Possibly, the variability of gray and white matter structures could both be related to the role of DCDC2 in ciliary function, which affects both neuronal migration and axonal outgrowth.

Mutation in CEP63 co-segregating with developmental dyslexia in a Swedish family.

Developmental dyslexia is the most common learning disorder in children. Problems in reading and writing are likely due to a complex interaction of genetic and environmental factors, resulting in reduced power of studies of the genetic factors underlying developmental dyslexia. Our approach in the current study was to perform exome sequencing of affected and unaffected individuals within an extended pedigree with a familial form of developmental dyslexia. We identified a two-base mutation, causing a p.R229L amino acid substitution in the centrosomal protein 63 kDa (CEP63), co-segregating with developmental dyslexia in this pedigree. This mutation is novel, and predicted to be highly damaging for the function of the protein. 3D modelling suggested a distinct conformational change caused by the mutation. CEP63 is localised to the centrosome in eukaryotic cells and is required for maintaining normal centriole duplication and control of cell cycle progression. We found that a common polymorphism in the CEP63 gene had a significant association with brain white matter volume. The brain regions were partly overlapping with the previously reported region influenced by polymorphisms in the dyslexia susceptibility genes DYX1C1 and KIAA0319. We hypothesise that CEP63 is particularly important for brain development and might control the proliferation and migration of cells when those two events need to be highly coordinated.

Variant Profiling of Candidate Genes in Pancreatic Ductal Adenocarcinoma.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) has a poor prognosis. Variant profiling is crucial for developing personalized treatment and elucidating the etiology of this disease. METHODS: Patients with PDAC undergoing surgery from 2007 to 2012 (n = 73) were followed from diagnosis until death or the end of the study. We applied an anchored multiplex PCR (AMP)-based next-generation sequencing (NGS) method to a panel of 65 selected genes and assessed analytical performance by sequencing a quantitative multiplex DNA reference standard. In clinical PDAC samples, detection of low-level KRAS (Kirsten rat sarcoma viral oncogene homolog) mutations was validated by allele-specific PCR and digital PCR. We compared overall survival of patients according to KRAS mutation status by log-rank test and applied logistic regression to evaluate the association between smoking and tumor variant types. RESULTS: The AMP-based NGS method could detect variants with allele frequencies as low as 1% given sufficient sequencing depth (>1500×). Low-frequency KRAS G12 mutations (allele frequency 1%-5%) were all confirmed by allele-specific PCR and digital PCR. The most prevalent genetic alterations were in KRAS (78% of patients), TP53 (tumor protein p53) (25%), and SMAD4 (SMAD family member 4) (8%). Overall survival in T3-stage PDAC patients differed among KRAS mutation subtypes (P = 0.019). Transversion variants were more common in ever-smokers than in never-smokers (odds ratio 5.7; 95% CI 1.2-27.8). CONCLUSIONS: The AMP-based NGS method is applicable for profiling tumor variants. Using this approach, we demonstrated that in PDAC patients, KRAS mutant subtype G12V is associated with poorer survival, and that transversion variants are more common among smokers.

Polymorphisms in DCDC2 and S100B associate with developmental dyslexia.

Genetic studies of complex traits have become increasingly successful as progress is made in next-generation sequencing. We aimed at discovering single nucleotide variation present in known and new candidate genes for developmental dyslexia: CYP19A1, DCDC2, DIP2A, DYX1C1, GCFC2 (also known as C2orf3), KIAA0319, MRPL19, PCNT, PRMT2, ROBO1 and S100B. We used next-generation sequencing to identify single-nucleotide polymorphisms in the exons of these 11 genes in pools of 100 DNA samples of Finnish individuals with developmental dyslexia. Subsequent individual genotyping of those 100 individuals, and additional cases and controls from the Finnish and German populations, validated 92 out of 111 different single-nucleotide variants. A nonsynonymous polymorphism in DCDC2 (corrected P = 0.002) and a noncoding variant in S100B (corrected P = 0.016) showed a significant association with spelling performance in families of German origin. No significant association was found for the variants neither in the Finnish case-control sample set nor in the Finnish family sample set. Our findings further strengthen the role of DCDC2 and implicate S100B, in the biology of reading and spelling.

Polymorphisms in the dopamine receptor 2 gene region influence improvements during working memory training in children and adolescents.

Studying the effects of cognitive training can lead to finding better treatments, but it can also be a tool for investigating factors important for brain plasticity and acquisition of cognitive skills. In this study, we investigated how single-nucleotide polymorphisms (SNPs) and ratings of intrinsic motivation were associated to interindividual differences in improvement during working memory training. The study included 256 children aged 7-19 years who were genotyped for 13 SNPs within or near eight candidate genes previously implicated in learning: COMT, SLC6A3 (DAT1), DRD4, DRD2, PPP1R1B (DARPP32), MAOA, LMX1A, and BDNF. Ratings on the intrinsic motivation inventory were also available for 156 of these children. All participants performed at least 20 sessions of working memory training, and performance during the training was logged and used as the outcome variable. We found that two SNPs, rs1800497 and rs2283265, located near and within the dopamine receptor 2 (DRD2) gene, respectively, were significantly associated with improvements during training (p < .003 and p < .0004, respectively). Scores from a questionnaire regarding intrinsic motivation did not correlate with training outcome. However, we observed both the main effect of genotype at those two loci as well as the interaction between genotypes and ratings of intrinsic motivation (perceived competence). Both SNPs have previously been shown to affect DRD2 receptor density primarily in the BG. Our results suggest that genetic variation is accounting for some interindividual differences in how children acquire cognitive skills and that part of this effect is also seen on intrinsic motivation. Moreover, they suggest that dopamine D2 transmission in the BG is a key factor for cognitive plasticity.

Mutations in the gene encoding fibroblast growth factor 10 are associated with aplasia of lacrimal and salivary glands.

Autosomal dominant aplasia of lacrimal and salivary glands (ALSG; OMIM 180920 and OMIM 103420) is a rare condition characterized by irritable eyes and dryness of the mouth. We mapped ALSG to 5p13.2-5q13.1, which coincides with the gene fibroblast growth factor 10 (FGF10). In two extended pedigrees, we identified heterozygous mutations in FGF10 in all individuals with ALSG. Fgf10(+/-) mice have a phenotype similar to ALSG, providing a model for this disorder. We suggest that haploinsufficiency for FGF10 during a crucial stage of development results in ALSG.

Case report: a novel deep intronic splice-altering variant in DMD as a cause of Becker muscular dystrophy.

We present the case of a male patient who was ultimately diagnosed with Becker muscular dystrophy (BMD; MIM# 300376) after the onset of muscle weakness in his teens progressively led to significant walking difficulties in his twenties. A genetic diagnosis was pursued but initial investigation revealed no aberrations in the dystrophin gene (DMD), although immunohistochemistry and Western blot analysis suggested the diagnosis of dystrophinopathy. Eventually, after more than 10 years, an RNA analysis captured abnormal splicing where 154 nucleotides from intron 43 were inserted between exon 43 and 44 resulting in a frameshift and a premature stop codon. Normal splicing of the DMD gene was also observed. Additionally, a novel variant c.6291-13537A>G in DMD was confirmed in the genomic DNA of the patient. The predicted function of the variant aligns with the mRNA results. To conclude, we here demonstrate that mRNA analysis can guide the diagnosis of non-coding genetic variants in DMD.

Functional interaction of DYX1C1 with estrogen receptors suggests involvement of hormonal pathways in dyslexia.

Dyslexia, or specific reading disability, is the unexpected failure in learning to read and write when intelligence and senses are normal. One of the susceptibility genes, DYX1C1, has been implicated in neuronal migration, but little is known about its interactions and functions. As DYX1C1 was suggested to interact with the U-box protein CHIP (carboxy terminus of Hsc70-interacting protein), which also participates in the degradation of estrogen receptors alpha (ERalpha) and beta (ERbeta), we hypothesized that the effects of DYX1C1 might be at least in part mediated through the regulation of ERs. ERs have shown to be important in brain development and cognitive functions. Indeed, we show that DYX1C1 interacts with both ERs in the presence of 17beta-estradiol, as determined by co-localization, co-immunoprecipitation and proximity ligation assays. Protein levels of endogenous ERalpha or exogenous ERbeta were reduced upon over-expression of DYX1C1, resulting in decreased transcriptional responses to 17beta-estradiol. Furthermore, we detected in vivo complexes of DYX1C1 with ERalpha or ERbeta at endogenous levels along neurites of primary rat hippocampal neurons. Taken together, our data suggest that DYX1C1 is involved in the regulation of ERalpha and ERbeta, and may thus affect the brain development and regulate cognitive functions. These findings provide novel insights into the function of DYX1C1 and link neuronal migration and developmental dyslexia to the estrogen-signaling effects in the brain.

SNP variations in the 7q33 region containing DGKI are associated with dyslexia in the Finnish and German populations.

Four genes, DYX1C1, ROBO1, DCDC2 and KIAA0319 have been studied both genetically and functionally as candidate genes for developmental dyslexia, a common learning disability in children. The identification of novel genes is crucial to better understand the molecular pathways affected in dyslectic individuals. Here, we report results from a fine-mapping approach involving linkage and association analysis in Finnish and German dyslexic cohorts. We restrict a candidate region to 0.3 Mb on chromosome 7q33. This region harbours the gene diacylglycerol kinase, iota (DGKI) which contains overlapping haplotypes associated with dyslexia in both Finnish and German sample sets.

Alpha-cardiac actin mutations produce atrial septal defects.

Atrial septal defect (ASD) is one of the most frequent congenital heart defects (CHDs) with a variable phenotypic effect depending on the size of the septal shunt. We identified two pedigrees comprising 20 members segregating isolated autosomal dominant secundum ASD. By genetic mapping, we identified the gene-encoding alpha-cardiac actin (ACTC1), which is essential for cardiac contraction, as the likely candidate. A mutation screen of the coding regions of ACTC1 revealed a founder mutation predicting an M123V substitution in affected individuals of both pedigrees. Functional analysis of ACTC1 with an M123V substitution shows a reduced affinity for myosin, but with retained actomyosin motor properties. We also screened 408 sporadic patients with CHDs and identified a case with ASD and a 17-bp deletion in ACTC1 predicting a non-functional protein. Morpholino (MO) knockdown of ACTC1 in chick embryos produces delayed looping and reduced atrial septa, supporting a developmental role for this protein. The combined results indicate, for the first time, that ACTC1 mutations or reduced ACTC1 levels may lead to ASD without signs of cardiomyopathy.

Multiple epiphyseal dysplasia.

BACKGROUND: Multiple epiphyseal dysplasia (MED) is a common genetically and clinically heterogeneous skeletal dysplasia characterized by early-onset osteoarthritis, mainly in the hip and knee, and mild-to-moderate short stature. Here we report on a 6-generation MED family with 17 affected members. METHOD: The clinical and radiographic data on the 12 affected members still living were scrutinized. A structured inquiry comprising state of health and MED-related symptoms since birth up to the present time and the osteoarthritis outcome (KOOS) questionnaire were sent to all living family members with MED. The 5 known gene loci for autosomal dominant MED were analyzed for linkage, using fluorescence-labeled microsatellite markers. Linkage was ascertained with markers close to the COL9A2 gene, which was analyzed for mutations by sequencing. RESULTS: We identified an exon 3 donor splice mutation in the COL9A2 gene in all affected family members. Clinical, radiographic, and questionnaire data from affected family members suggested that MED caused by COL9A2 mutations starts in early childhood with knee pain accompanied by delayed ossification of femoral epiphyses. The disease then either stabilizes during puberty or progresses with additional joints becoming affected; joint surgery might be necessary. The progression of the disease also affects muscles, with increasing atrophy, resulting in muscle fatigue and pain. Muscular atrophy has not been reported earlier in cases with COL9A2 mutations. INTERPRETATION: In a patient with clinically suspected or verified MED, it is important to perform DNA-based analysis to identify a possible disease-causing mutation. This information can be used to carry out genetic risk assessment of other family members and to achieve an early and correct diagnosis in the children.

Evaluation of the ISL1 gene in the pathogenesis of bladder exstrophy in a Swedish cohort.

Bladder exstrophy is a congenital closure defect of the urinary bladder with a profound effect on morbidity. Although the malformation is usually sporadic, a genetic background is supported by an increased recurrence risk in relatives, higher concordance rates in monozygotic twins and several associated chromosomal aberrations. Recently, the ISL1 gene was presented as a candidate gene for bladder exstrophy and epispadias complex (BEEC) development in two different studies. In our study, we screened for genetic variants in the ISL1 gene in DNA from 125 Swedish patients using Sanger sequencing and array-CGH analysis. In addition, we evaluated ISL1 expression in RNA of human bladder during embryonic and fetal weeks 5-10 relative to that in lung tissue (week 9). In total, 21 single-nucleotide variants were identified, including a potentially novel missense variant, c.137C>G p.(Ala46Gly), substituting a conserved amino acid. This variant was inherited from an unaffected mother. No structural variants were identified. RNA sequencing revealed ISL1 mRNA expression during the critical time frame of human bladder development. In conclusion, we did not detect any known or likely pathogenic variants in the ISL1 gene in 125 Swedish BEEC patients, indicating that variation in the ISL1 gene is not a common genetic mechanism of BEEC development in the Swedish population.

Genome-wide association study of Hirschsprung disease detects a novel low-frequency variant at the RET locus.

Hirschsprung disease (HSCR) is a congenital disorder with a population incidence of ~1/5000 live births, defined by an absence of enteric ganglia along variable lengths of the colon. HSCR genome-wide association studies (GWAS) have found common associated variants at RET, SEMA3, and NRG1, but they still fail to explain all of its heritability. To enhance gene discovery, we performed a GWAS of 170 cases identified from the Danish nationwide pathology registry with 4717 controls, based on 6.2 million variants imputed from the haplotype reference consortium panel. We found a novel low-frequency variant (rs144432435), which, when conditioning on the lead RET single-nucleotide polymorphism (SNP), was of genome-wide significance in the discovery analysis. This conditional association signal was replicated in a Swedish HSCR cohort with discovery plus replication meta-analysis conditional odds ratio of 6.6 (P = 7.7 × 10-10; 322 cases and 4893 controls). The conditional signal was, however, not replicated in two HSCR cohorts from USA and Finland, leading to the hypothesis that rs144432435 tags a rare haplotype present in Denmark and Sweden. Using the genome-wide complex trait analysis method, we estimated the SNP heritability of HSCR to be 88%, close to estimates based on classical family studies. Moreover, by using Lasso (least absolute shrinkage and selection operator) regression we were able to construct a genetic HSCR predictor with a area under the receiver operator characteristics curve of 76% in an independent validation set. In conclusion, we combined the largest collection of sporadic Hirschsprung cases to date (586 cases) to further elucidate HSCR's genetic architecture.

Targeted disruption of the ribosomal protein S19 gene is lethal prior to implantation.

The ribosomal protein S19 (RPS19) is located in the small (40S) subunit and is one of 79 ribosomal proteins. The gene encoding RPS19 is mutated in approximately 25% of patients with Diamond-Blackfan anemia, which is a rare congenital erythroblastopenia. Affected individuals present with decreased numbers or the absence of erythroid precursors in the bone marrow, and associated malformations of various organs are common. We produced C57BL/6J mice with a targeted disruption of murine Rps19 to study its role in erythropoiesis and development. Mice homozygous for the disrupted Rps19 were not identified as early as the blastocyst stage, indicating a lethal effect. In contrast, mice heterozygous for the disrupted Rps19 allele have normal growth and organ development, including that of the hematopoietic system. Our findings indicate that zygotes which are Rps19(-/-) do not form blastocysts, whereas one normal Rps19 allele in C57BL/6J mice is sufficient to maintain normal ribosomal and possibly extraribosomal functions.

Human ROBO1 regulates white matter structure in corpus callosum.

The axon guidance receptor, Robo1, controls the pathfinding of callosal axons in mice. To determine whether the orthologous ROBO1 gene is involved in callosal development also in humans, we studied polymorphisms in the ROBO1 gene and variation in the white matter structure in the corpus callosum using both structural magnetic resonance imaging and diffusion tensor magnetic resonance imaging. We found that five polymorphisms in the regulatory region of ROBO1 were associated with white matter density in the posterior part of the corpus callosum pathways. One of the polymorphisms, rs7631357, was also significantly associated with the probability of connections to the parietal cortical regions. Our results demonstrate that human ROBO1 may be involved in the regulation of the structure and connectivity of posterior part of corpus callosum.