Monoallelic and biallelic mutations in RELN underlie a graded series of neurodevelopmental disorders

Reelin, a large extracellular protein, plays several critical roles in brain development and function. It is encoded by RELN, first identified as the gene disrupted in the reeler mouse, a classic neurological mutant exhibiting ataxia, tremors and a 'reeling' gait. In humans, biallelic variants in RELN have been associated with a recessive lissencephaly variant with cerebellar hypoplasia, which matches well with the homozygous mouse mutant that has abnormal cortical structure, small hippocampi and severe cerebellar hypoplasia. Despite the large size of the gene, only 11 individuals with RELN-related lissencephaly with cerebellar hypoplasia from six families have previously been reported. Heterozygous carriers in these families were briefly reported as unaffected, although putative loss-of-function variants are practically absent in the population (probability of loss of function intolerance = 1). Here we present data on seven individuals from four families with biallelic and 13 individuals from seven families with monoallelic (heterozygous) variants of RELN and frontotemporal or temporal-predominant lissencephaly variant. Some individuals with monoallelic variants have moderate frontotemporal lissencephaly, but with normal cerebellar structure and intellectual disability with severe behavioural dysfunction. However, one adult had abnormal MRI with normal intelligence and neurological profile. Thorough literature analysis supports a causal role for monoallelic RELN variants in four seemingly distinct phenotypes including frontotemporal lissencephaly, epilepsy, autism and probably schizophrenia. Notably, we observed a significantly higher proportion of loss-of-function variants in the biallelic compared to the monoallelic cohort, where the variant spectrum included missense and splice-site variants. We assessed the impact of two canonical splice-site variants observed as biallelic or monoallelic variants in individuals with moderately affected or normal cerebellum and demonstrated exon skipping causing in-frame loss of 46 or 52 amino acids in the central RELN domain. Previously reported functional studies demonstrated severe reduction in overall RELN secretion caused by heterozygous missense variants p.Cys539Arg and p.Arg3207Cys associated with lissencephaly suggesting a dominant-negative effect. We conclude that biallelic variants resulting in complete absence of RELN expression are associated with a consistent and severe phenotype that includes cerebellar hypoplasia. However, reduced expression of RELN remains sufficient to maintain nearly normal cerebellar structure. Monoallelic variants are associated with incomplete penetrance and variable expressivity even within the same family and may have dominant-negative effects. Reduced RELN secretion in heterozygous individuals affects only cortical structure whereas the cerebellum remains intact. Our data expand the spectrum of RELN-related neurodevelopmental disorders ranging from lethal brain malformations to adult phenotypes with normal brain imaging.

Anti-tumor activity of a T-helper 1 multiantigen vaccine in a murine model of prostate cancer

Prostate cancer is one of the few malignancies that includes vaccination as a treatment modality. Elements of an effective cancer vaccine should include the ability to elicit a Type I T-cell response and target multiple antigenic proteins expressed early in the disease. Using existing gene datasets encompassing normal prostate tissue and tumors with Gleason Score ≤ 6 and ≥ 8, 10 genes were identified that were upregulated and conserved in prostate cancer regardless of the aggressiveness of disease. These genes encoded proteins also expressed in prostatic intraepithelial neoplasia. Putative Class II epitopes derived from these proteins were predicted by a combination of algorithms and, using human peripheral blood, epitopes which selectively elicited IFN-γ or IL-10 dominant antigen specific cytokine secretion were determined. Th1 selective epitopes were identified for eight antigens. Epitopes from three antigens elicited Th1 dominant immunity in mice; PSMA, HPN, and AMACR. Each single antigen vaccine demonstrated significant anti-tumor activity inhibiting growth of implanted Myc-Cap cells after immunization as compared to control. Immunization with the combination of antigens, however, was superior to each alone in controlling tumor growth. When vaccination occurred simultaneously to tumor implant, multiantigen immunized mice had significantly smaller tumors than controls (p = 0.002) and a significantly improved overall survival (p = 0.0006). This multiantigen vaccine shows anti-tumor activity in a murine model of prostate cancer.

Profiling PI3K-AKT-MTOR variants in focal brain malformations reveals new insights for diagnostic care

Focal malformations of cortical development including focal cortical dysplasia, hemimegalencephaly and megalencephaly, are a spectrum of neurodevelopmental disorders associated with brain overgrowth, cellular and architectural dysplasia, intractable epilepsy, autism and intellectual disability. Importantly, focal cortical dysplasia is the most common cause of focal intractable paediatric epilepsy. Gain and loss of function variants in the PI3K-AKT-MTOR pathway have been identified in this spectrum, with variable levels of mosaicism and tissue distribution. In this study, we performed deep molecular profiling of common PI3K-AKT-MTOR pathway variants in surgically resected tissues using droplet digital polymerase chain reaction (ddPCR), combined with analysis of key phenotype data. A total of 159 samples, including 124 brain tissue samples, were collected from 58 children with focal malformations of cortical development. We designed an ultra-sensitive and highly targeted molecular diagnostic panel using ddPCR for six mutational hotspots in three PI3K-AKT-MTOR pathway genes, namely PIK3CA (p.E542K, p.E545K, p.H1047R), AKT3 (p.E17K) and MTOR (p.S2215F, p.S2215Y). We quantified the level of mosaicism across all samples and correlated genotypes with key clinical, neuroimaging and histopathological data. Pathogenic variants were identified in 17 individuals, with an overall molecular solve rate of 29.31%. Variant allele fractions ranged from 0.14 to 22.67% across all mutation-positive samples. Our data show that pathogenic MTOR variants are mostly associated with focal cortical dysplasia, whereas pathogenic PIK3CA variants are more frequent in hemimegalencephaly. Further, the presence of one of these hotspot mutations correlated with earlier onset of epilepsy. However, levels of mosaicism did not correlate with the severity of the cortical malformation by neuroimaging or histopathology. Importantly, we could not identify these mutational hotspots in other types of surgically resected epileptic lesions (e.g. polymicrogyria or mesial temporal sclerosis) suggesting that PI3K-AKT-MTOR mutations are specifically causal in the focal cortical dysplasia-hemimegalencephaly spectrum. Finally, our data suggest that ultra-sensitive molecular profiling of the most common PI3K-AKT-MTOR mutations by targeted sequencing droplet digital polymerase chain reaction is an effective molecular approach for these disorders with a good diagnostic yield when paired with neuroimaging and histopathology.

Machine Learning Prediction of Non-Coding Variant Impact in Human Retinal cis-Regulatory Elements

Purpose

Prior studies have demonstrated the significance of specific cis-regulatory variants in retinal disease; however, determining the functional impact of regulatory variants remains a major challenge. In this study, we utilized a machine learning approach, trained on epigenomic data from the adult human retina, to systematically quantify the predicted impact of cis-regulatory variants.

Methods

We used human retinal DNA accessibility data (ATAC-seq) to determine a set of 18.9k high-confidence, putative cis-regulatory elements. Eighty percent of these elements were used to train a machine learning model utilizing a gapped k-mer support vector machine-based approach. In silico saturation mutagenesis and variant scoring was applied to predict the functional impact of all potential single nucleotide variants within cis-regulatory elements. Impact scores were tested in a 20% hold-out dataset and compared to allele population frequency, phylogenetic conservation, transcription factor (TF) binding motifs, and existing massively parallel reporter assay data.

Results

We generated a model that distinguishes between human retinal regulatory elements and negative test sequences with 95% accuracy. Among a hold-out test set of 3.7k human retinal CREs, all possible single nucleotide variants were scored. Variants with negative impact scores correlated with higher phylogenetic conservation of the reference allele, disruption of predicted TF binding motifs, and massively parallel reporter expression.

Conclusions

We demonstrated the utility of human retinal epigenomic data to train a machine learning model for the purpose of predicting the impact of non-coding regulatory sequence variants. Our model accurately scored sequences and predicted putative transcription factor binding motifs. This approach has the potential to expedite the characterization of pathogenic non-coding sequence variants in the context of unexplained retinal disease.

Translational Relevance

This workflow and resulting dataset serve as a promising genomic tool to facilitate the clinical prioritization of functionally disruptive non-coding mutations in the retina.

Cell-specific cis-regulatory elements and mechanisms of non-coding genetic disease in human retina and retinal organoids

Cis-regulatory elements (CREs) play a critical role in the development and disease-states of all human cell types. In the retina, CREs have been implicated in several inherited disorders. To better characterize human retinal CREs, we performed single-nucleus assay for transposase-accessible chromatin sequencing (snATAC-seq) and single-nucleus RNA sequencing (snRNA-seq) on the developing and adult human retina and on induced pluripotent stem cell (iPSC)-derived retinal organoids. These analyses identified developmentally dynamic, cell-class-specific CREs, enriched transcription-factor-binding motifs, and putative target genes. CREs in the retina and organoids are highly correlated at the single-cell level, and this supports the use of organoids as a model for studying disease-associated CREs. As a proof of concept, we disrupted a disease-associated CRE at 5q14.3, confirming its principal target gene as the miR-9-2 primary transcript and demonstrating its role in neurogenesis and gene regulation in mature glia. This study provides a resource for characterizing human retinal CREs and showcases organoids as a model to study the function of CREs that influence development and disease.

Activated interleukin-7 receptor signaling drives B-cell acute lymphoblastic leukemia in mice

Philadelphia chromosome-like acute lymphoblastic leukemia (Ph-like ALL) is a high-risk subtype of B-ALL often associated with genetic variants that alter cytokine receptor signaling, including mutations in the interleukin-7 receptor (IL7R). To investigate whether IL7R variants are leukemia-initiating, we built mouse models expressing activated Il7r (aIL7R). B-cell intrinsic aIL7R mice developed spontaneous B-ALL, demonstrating sufficiency of Il7r activating mutations in leukemogenesis. Concomitant introduction of a knock-out allele in the associated adapter protein Lnk (encoded by Sh2b3) or a dominant-negative variant of the transcription factor Ikaros (Ikzf1) increased disease penetrance. The resulting murine leukemias displayed monoclonality and recurrent somatic Kras mutations and efficiently engrafted into immunocompetent mice. Phosphoproteomic analyses of aIL7R leukemic cells revealed constitutive Stat5 signaling and B cell receptor (BCR)-like signaling despite the absence of surface pre-BCR. Finally, in vitro treatment of aIL7R leukemic B-cells with Jak, mTOR, or Syk inhibitors blocked growth, confirming that each pathway is active in this mouse model of IL7R-driven B-ALL.

The Autoimmune Risk R262W Variant of the Adaptor SH2B3 Improves Survival in Sepsis

The single-nucleotide polymorphism (SNP) rs3184504 is broadly associated with increased risk for multiple autoimmune and cardiovascular diseases. Although the allele is uniquely enriched in European descent, the mechanism for the widespread selective sweep is not clear. In this study, we find the rs3184504*T allele had a strong association with reduced mortality in a human sepsis cohort. The rs3184504*T allele associates with a loss-of-function amino acid change (p.R262W) in the adaptor protein SH2B3, a likely causal variant. To better understand the role of SH2B3 in sepsis, we used mouse modeling and challenged SH2B3-deficient mice with a polymicrobial cecal-ligation puncture (CLP) procedure. We found SH2B3 deficiency improved survival and morbidity with less organ damage and earlier bacterial clearance compared with control mice. The peritoneal infiltrating cells exhibited augmented phagocytosis in Sh2b3 ^-/- mice with enriched recruitment of Ly6C^hi inflammatory monocytes despite equivalent or reduced chemokine expression. Rapid cycling of monocytes and progenitors occurred uniquely in the Sh2b3 ^-/- mice following CLP, suggesting augmented myelopoiesis. To model the hypomorphic autoimmune risk allele, we created a novel knockin mouse harboring a similar point mutation in the murine pleckstrin homology domain of SH2B3. At baseline, phenotypic changes suggested a hypomorphic allele. In the CLP model, homozygous knockin mice displayed improved mortality and morbidity compared with wild-type or heterozygous mice. Collectively, these data suggest that hypomorphic SH2B3 improves the sepsis response and that balancing selection likely contributed to the relative frequency of the autoimmune risk variant.

Spiny mice activate unique transcriptional programs after severe kidney injury regenerating organ function without fibrosis

Fibrosis-driven solid organ failure is an enormous burden on global health. Spiny mice (Acomys) are terrestrial mammals that can regenerate severe skin wounds without scars to avoid predation. Whether spiny mice also regenerate internal organ injuries is unknown. Here, we show that despite equivalent acute obstructive or ischemic kidney injury, spiny mice fully regenerate nephron structure and organ function without fibrosis, whereas C57Bl/6 or CD1 mice progress to complete organ failure with extensive renal fibrosis. Two mechanisms for vertebrate regeneration have been proposed that emphasize either extrinsic (pro-regenerative macrophages) or intrinsic (surviving cells of the organ itself) controls. Comparative transcriptome analysis revealed that the Acomys genome appears poised at the time of injury to initiate regeneration by surviving kidney cells, whereas macrophage accumulation was not detected until about day 7. Thus, we provide evidence for rapid activation of a gene expression signature for regenerative wound healing in the spiny mouse kidney.

Gene regulatory networks controlling temporal patterning, neurogenesis, and cell-fate specification in mammalian retina

Gene regulatory networks (GRNs), consisting of transcription factors and their target sites, control neurogenesis and cell-fate specification in the developing central nervous system. In this study, we use integrated single-cell RNA and single-cell ATAC sequencing (scATAC-seq) analysis in developing mouse and human retina to identify multiple interconnected, evolutionarily conserved GRNs composed of cell-type-specific transcription factors that both activate genes within their own network and inhibit genes in other networks. These GRNs control temporal patterning in primary progenitors, regulate transition from primary to neurogenic progenitors, and drive specification of each major retinal cell type. We confirm that NFI transcription factors selectively activate expression of genes promoting late-stage temporal identity in primary retinal progenitors and identify other transcription factors that regulate rod photoreceptor specification in postnatal retina. This study inventories cis- and trans-acting factors that control retinal development and can guide cell-based therapies aimed at replacing retinal neurons lost to disease.

Spatial and cell type transcriptional landscape of human cerebellar development

The human neonatal cerebellum is one-fourth of its adult size yet contains the blueprint required to integrate environmental cues with developing motor, cognitive and emotional skills into adulthood. Although mature cerebellar neuroanatomy is well studied, understanding of its developmental origins is limited. In this study, we systematically mapped the molecular, cellular and spatial composition of human fetal cerebellum by combining laser capture microscopy and SPLiT-seq single-nucleus transcriptomics. We profiled functionally distinct regions and gene expression dynamics within cell types and across development. The resulting cell atlas demonstrates that the molecular organization of the cerebellar anlage recapitulates cytoarchitecturally distinct regions and developmentally transient cell types that are distinct from the mouse cerebellum. By mapping genes dominant for pediatric and adult neurological disorders onto our dataset, we identify relevant cell types underlying disease mechanisms. These data provide a resource for probing the cellular basis of human cerebellar development and disease.

MYT1 role in the microtia-craniofacial microsomia spectrum

Background

Craniofacial microsomia (CFM), also known as the oculo‐auriculo‐vertebral spectrum, comprises a variable phenotype with the most common features including microtia and mandibular hypoplasia on one or both sides, in addition to lateral oral clefts, epibulbar dermoids, cardiac, vertebral, and renal abnormalities. The etiology of CFM is largely unknown. The MYT1 gene has been reported as a candidate based in mutations found in three unrelated individuals. Additional patients with mutations in this gene are required to establish its causality. We present two individuals with CFM that have rare variants in MYT1 contributing to better understand the genotype and phenotype associated with mutations in this gene.

Methods/Results

We conducted genetic analysis using whole‐exome and ‐genome sequencing in 128 trios with CFM. Two novel MYT1 mutations were identified in two participants. Sanger sequencing was used to confirm these mutations.

Conclusion

We identified two additional individuals with CFM who carry rare variants in MYT1, further supporting the presumptive role of this gene in the CFM spectrum.

Transcriptome data of temporal and cingulate cortex in the Rett syndrome brain

Rett syndrome is an X-linked neurodevelopmental disorder caused by mutation in the methyl-CpG-binding protein 2 gene (MECP2) in the majority of cases. We describe an RNA sequencing dataset of postmortem brain tissue samples from four females clinically diagnosed with Rett syndrome and four age-matched female donors. The dataset contains 16 transcriptomes, including two brain regions, temporal and cingulate cortex, for each individual. We compared our dataset with published transcriptomic analyses of postmortem brain tissue from Rett syndrome and found consistent gene expression alterations among regions of the cerebral cortex. Our data provide a valuable resource to explore the biology of the human brain in Rett syndrome.

Mapping the cis-regulatory architecture of the human retina reveals noncoding genetic variation in disease

The interplay of transcription factors and cis-regulatory elements (CREs) orchestrates the dynamic and diverse genetic programs that assemble the human central nervous system (CNS) during development and maintain its function throughout life. Genetic variation within CREs plays a central role in phenotypic variation in complex traits including the risk of developing disease. We took advantage of the retina, a well-characterized region of the CNS known to be affected by pathogenic variants in CREs, to establish a roadmap for characterizing regulatory variation in the human CNS. This comprehensive analysis of tissue-specific regulatory elements, transcription factor binding, and gene expression programs in three regions of the human visual system (retina, macula, and retinal pigment epithelium/choroid) reveals features of regulatory element evolution that shape tissue-specific gene expression programs and defines regulatory elements with the potential to contribute to Mendelian and complex disorders of human vision.

Genotype correlates with clinical severity in PIK3CA-associated lymphatic malformations

Lymphatic malformations (LMs) are congenital, nonneoplastic vascular malformations associated with postzygotic activating PIK3CA mutations. The mutation spectrum within LMs is narrow, with the majority having 1 of 3 hotspot mutations. Despite this relative genetic homogeneity, clinical presentations differ dramatically. We used molecular inversion probes and droplet digital polymerase chain reaction to perform deep, targeted sequencing of PIK3CA in 271 affected and unaffected tissue samples from 81 individuals with isolated LMs and retrospectively collected clinical data. Pathogenic PIK3CA mutations were identified in affected LM tissue in 64 individuals (79%) with isolated LMs, with variant allele fractions (VAFs) ranging from 0.1% to 13%. Initial analyses revealed no correlation between VAF and phenotype variables. Recognizing that different mutations activate PI3K to varying degrees, we developed a metric, the genotype-adjusted VAF (GVAF), to account for differences in mutation strength, and found significantly higher GVAFs in LMs with more severe clinical characteristics including orofacial location or microcystic structure. In addition to providing insight into LM pathogenesis, we believe GVAF may have broad applicability for genotype-phenotype analyses in mosaic disorders.

Spatiotemporal expansion of primary progenitor zones in the developing human cerebellum

We present histological and molecular analyses of the developing human cerebellum from 30 days after conception to 9 months after birth. Differences in developmental patterns between humans and mice include spatiotemporal expansion of both ventricular and rhombic lip primary progenitor zones to include subventricular zones containing basal progenitors. The human rhombic lip persists longer through cerebellar development than in the mouse and undergoes morphological changes to form a progenitor pool in the posterior lobule, which is not seen in other organisms, not even in the nonhuman primate the macaque. Disruptions in human rhombic lip development are associated with posterior cerebellar vermis hypoplasia and Dandy-Walker malformation. The presence of these species-specific neural progenitor populations refines our insight into human cerebellar developmental disorders.

Redefining the Etiologic Landscape of Cerebellar Malformations

Cerebellar malformations are diverse congenital anomalies frequently associated with developmental disability. Although genetic and prenatal non-genetic causes have been described, no systematic analysis has been performed. Here, we present a large-exome sequencing study of Dandy-Walker malformation (DWM) and cerebellar hypoplasia (CBLH). We performed exome sequencing in 282 individuals from 100 families with DWM or CBLH, and we established a molecular diagnosis in 36 of 100 families, with a significantly higher yield for CBLH (51%) than for DWM (16%). The 41 variants impact 27 neurodevelopmental-disorder-associated genes, thus demonstrating that CBLH and DWM are often features of monogenic neurodevelopmental disorders. Though only seven monogenic causes (19%) were identified in more than one individual, neuroimaging review of 131 additional individuals confirmed cerebellar abnormalities in 23 of 27 genetic disorders (85%). Prenatal risk factors were frequently found among individuals without a genetic diagnosis (30 of 64 individuals [47%]). Single-cell RNA sequencing of prenatal human cerebellar tissue revealed gene enrichment in neuronal and vascular cell types; this suggests that defective vasculogenesis may disrupt cerebellar development. Further, de novo gain-of-function variants in PDGFRB, a tyrosine kinase receptor essential for vascular progenitor signaling, were associated with CBLH, and this discovery links genetic and non-genetic etiologies. Our results suggest that genetic defects impact specific cerebellar cell types and implicate abnormal vascular development as a mechanism for cerebellar malformations. We also confirmed a major contribution for non-genetic prenatal factors in individuals with cerebellar abnormalities, substantially influencing diagnostic evaluation and counseling regarding recurrence risk and prognosis.

The Effect of Mouse Strain, Sex, and Carcinogen Dose on Toxicity and the Development of Lung Dysplasia and Squamous Cell Carcinomas in Mice

In order to translate new treatments to the clinic, it is necessary to use animal models that closely recapitulate human disease. Lung cancer develops after extended exposure to carcinogens. It has one of the highest mutation rates of all cancer and is highly heterogenic. Topical treatment with N-nitrosotris-(2-chloroethyl)urea (NTCU) induces lung squamous cell carcinoma (SCC) with nonsynonymous mutation rates similar to those reported for human non-small cell lung cancer. However, NTCU induces lung cancer with variable efficacy and toxicity depending on the mouse strain. A detailed characterization of the NTCU model is needed. We have compared the effect of three different NTCU doses (20, 30, and 40 mmol/L) in female and male of NIH Swiss, Black Swiss, and FVB mice on tumor incidence, survival, and toxicity. The main findings in this study are (1) NIH Swiss mice present with a higher incidence of SCC and lower mortality compared with Black Swiss and FVB mice; (2) 30 mmol/L NTCU dose induces SCC at the same rate and incidence as the 40 mmol/L dose with lower mortality; (3) female mice present higher grade and incidence of preinvasive lesions and SCC compared with males; (4) NTCU-induced transformation is principally within the respiratory system; and (5) NTCU treatment does not affect the ability to elicit a specific adaptive immune response. This study provides a reference point for experimental designs to evaluate either preventive or therapeutic treatments for lung SCC, including immunotherapies, before initiating human clinical trials.

MACF1 Mutations Encoding Highly Conserved Zinc-Binding Residues of the GAR Domain Cause Defects in Neuronal Migration and Axon Guidance

To date, mutations in 15 actin- or microtubule-associated genes have been associated with the cortical malformation lissencephaly and variable brainstem hypoplasia. During a multicenter review, we recognized a rare lissencephaly variant with a complex brainstem malformation in three unrelated children. We searched our large brain-malformation databases and found another five children with this malformation (as well as one with a less severe variant), analyzed available whole-exome or -genome sequencing data, and tested ciliogenesis in two affected individuals. The brain malformation comprised posterior predominant lissencephaly and midline crossing defects consisting of absent anterior commissure and a striking W-shaped brainstem malformation caused by small or absent pontine crossing fibers. We discovered heterozygous de novo missense variants or an in-frame deletion involving highly conserved zinc-binding residues within the GAR domain of MACF1 in the first eight subjects. We studied cilium formation and found a higher proportion of mutant cells with short cilia than of control cells with short cilia. A ninth child had similar lissencephaly but only subtle brainstem dysplasia associated with a heterozygous de novo missense variant in the spectrin repeat domain of MACF1. Thus, we report variants of the microtubule-binding GAR domain of MACF1 as the cause of a distinctive and most likely pathognomonic brain malformation. A gain-of-function or dominant-negative mechanism appears likely given that many heterozygous mutations leading to protein truncation are included in the ExAC Browser. However, three de novo variants in MACF1 have been observed in large schizophrenia cohorts.

PIK3CA-associated developmental disorders exhibit distinct classes of mutations with variable expression and tissue distribution

Mosaicism is increasingly recognized as a cause of developmental disorders with the advent of next-generation sequencing (NGS). Mosaic mutations of PIK3CA have been associated with the widest spectrum of phenotypes associated with overgrowth and vascular malformations. We performed targeted NGS using 2 independent deep-coverage methods that utilize molecular inversion probes and amplicon sequencing in a cohort of 241 samples from 181 individuals with brain and/or body overgrowth. We identified PIK3CA mutations in 60 individuals. Several other individuals (n = 12) were identified separately to have mutations in PIK3CA by clinical targeted-panel testing (n = 6), whole-exome sequencing (n = 5), or Sanger sequencing (n = 1). Based on the clinical and molecular features, this cohort segregated into three distinct groups: (a) severe focal overgrowth due to low-level but highly activating (hotspot) mutations, (b) predominantly brain overgrowth and less severe somatic overgrowth due to less-activating mutations, and (c) intermediate phenotypes (capillary malformations with overgrowth) with intermediately activating mutations. Sixteen of 29 PIK3CA mutations were novel. We also identified constitutional PIK3CA mutations in 10 patients. Our molecular data, combined with review of the literature, show that PIK3CA-related overgrowth disorders comprise a discontinuous spectrum of disorders that correlate with the severity and distribution of mutations.

The clinical and molecular spectrum of PIK3CA-related developmental disorders are correlated with types of mutations, tissue distributions, and levels of mosaicism with the clinical phenotype.

Bi-allelic mutations of CCDC88C are a rare cause of severe congenital hydrocephalus

Congenital or infantile hydrocephalus is caused by genetic and non-genetic factors and is highly heterogeneous in etiology. In recent studies, a limited number of genetic causes of hydrocephalus have been identified. To date, recessive mutations in the CCDC88C gene have been identified as a cause of non-syndromic congenital hydrocephalus in three reported families. Here, we report the fourth known family with two affected individuals with congenital hydrocephalus due to a homozygous mutation in the CCDC88C gene identified by whole exome sequencing. Our two newly described children, as well as the previously published ones, all shared several features including severe infantile-onset hydrocephalus, mild to severe intellectual delay, varying degrees of motor delay, and infantile onset seizures. All identified homozygous mutations in CCDC88C abolish the PDZ binding site necessary for proper CCDC88C protein function in the Wnt signaling pathway. Our report further establishes CCDC88C as one of the few known recessive causes of severe prenatal-onset hydrocephalus. Recognition of this syndrome has important diagnostic and genetic implications for families identified in the future.

Single suture craniosynostosis: Identification of rare variants in genes associated with syndromic forms

We report RNA-Sequencing results on a cohort of patients with single suture craniosynostosis and demonstrate significant enrichment of heterozygous, rare, and damaging variants among key craniosynostosis-related genes. Genetic burden analysis identified a significant increase in damaging variants in ATR, EFNA4, ERF, MEGF8, SCARF2, and TGFBR2. Of 391 participants, 15% were found to have damaging and potentially causal variants in 29 genes. We observed transmission in 96% of the affected individuals, and thus penetrance, epigenetics, and oligogenic factors need to be considered when recommending genetic testing in patients with nonsyndromic craniosynostosis.

Mutations of AKT3 are associated with a wide spectrum of developmental disorders including extreme megalencephaly

Mutations of genes within the phosphatidylinositol-3-kinase (PI3K)-AKT-MTOR pathway are well known causes of brain overgrowth (megalencephaly) as well as segmental cortical dysplasia (such as hemimegalencephaly, focal cortical dysplasia and polymicrogyria). Mutations of the AKT3 gene have been reported in a few individuals with brain malformations, to date. Therefore, our understanding regarding the clinical and molecular spectrum associated with mutations of this critical gene is limited, with no clear genotype-phenotype correlations. We sought to further delineate this spectrum, study levels of mosaicism and identify genotype-phenotype correlations of AKT3-related disorders. We performed targeted sequencing of AKT3 on individuals with these phenotypes by molecular inversion probes and/or Sanger sequencing to determine the type and level of mosaicism of mutations. We analysed all clinical and brain imaging data of mutation-positive individuals including neuropathological analysis in one instance. We performed ex vivo kinase assays on AKT3 engineered with the patient mutations and examined the phospholipid binding profile of pleckstrin homology domain localizing mutations. We identified 14 new individuals with AKT3 mutations with several phenotypes dependent on the type of mutation and level of mosaicism. Our comprehensive clinical characterization, and review of all previously published patients, broadly segregates individuals with AKT3 mutations into two groups: patients with highly asymmetric cortical dysplasia caused by the common p.E17K mutation, and patients with constitutional AKT3 mutations exhibiting more variable phenotypes including bilateral cortical malformations, polymicrogyria, periventricular nodular heterotopia and diffuse megalencephaly without cortical dysplasia. All mutations increased kinase activity, and pleckstrin homology domain mutants exhibited enhanced phospholipid binding. Overall, our study shows that activating mutations of the critical AKT3 gene are associated with a wide spectrum of brain involvement ranging from focal or segmental brain malformations (such as hemimegalencephaly and polymicrogyria) predominantly due to mosaic AKT3 mutations, to diffuse bilateral cortical malformations, megalencephaly and heterotopia due to constitutional AKT3 mutations. We also provide the first detailed neuropathological examination of a child with extreme megalencephaly due to a constitutional AKT3 mutation. This child has one of the largest documented paediatric brain sizes, to our knowledge. Finally, our data show that constitutional AKT3 mutations are associated with megalencephaly, with or without autism, similar to PTEN-related disorders. Recognition of this broad clinical and molecular spectrum of AKT3 mutations is important for providing early diagnosis and appropriate management of affected individuals, and will facilitate targeted design of future human clinical trials using PI3K-AKT pathway inhibitors.

Association of MTOR Mutations With Developmental Brain Disorders, Including Megalencephaly, Focal Cortical Dysplasia, and Pigmentary Mosaicism

IMPORTANCE

Focal cortical dysplasia (FCD), hemimegalencephaly, and megalencephaly constitute a spectrum of malformations of cortical development with shared neuropathologic features. These disorders are associated with significant childhood morbidity and mortality.

OBJECTIVE

To identify the underlying molecular cause of FCD, hemimegalencephaly, and diffuse megalencephaly.

DESIGN, SETTING, AND PARTICIPANTS

Patients with FCD, hemimegalencephaly, or megalencephaly (mean age, 11.7 years; range, 2-32 years) were recruited from Pediatric Hospital A. Meyer, the University of Hong Kong, and Seattle Children's Research Institute from June 2012 to June 2014. Whole-exome sequencing (WES) was performed on 8 children with FCD or hemimegalencephaly using standard-depth (50-60X) sequencing in peripheral samples (blood, saliva, or skin) from the affected child and their parents and deep (150-180X) sequencing in affected brain tissue. Targeted sequencing and WES were used to screen 93 children with molecularly unexplained diffuse or focal brain overgrowth. Histopathologic and functional assays of phosphatidylinositol 3-kinase-AKT (serine/threonine kinase)-mammalian target of rapamycin (mTOR) pathway activity in resected brain tissue and cultured neurons were performed to validate mutations.

MAIN OUTCOMES AND MEASURES

Whole-exome sequencing and targeted sequencing identified variants associated with this spectrum of developmental brain disorders.

RESULTS

Low-level mosaic mutations of MTOR were identified in brain tissue in 4 children with FCD type 2a with alternative allele fractions ranging from 0.012 to 0.086. Intermediate-level mosaic mutation of MTOR (p.Thr1977Ile) was also identified in 3 unrelated children with diffuse megalencephaly and pigmentary mosaicism in skin. Finally, a constitutional de novo mutation of MTOR (p.Glu1799Lys) was identified in 3 unrelated children with diffuse megalencephaly and intellectual disability. Molecular and functional analysis in 2 children with FCD2a from whom multiple affected brain tissue samples were available revealed a mutation gradient with an epicenter in the most epileptogenic area. When expressed in cultured neurons, all MTOR mutations identified here drive constitutive activation of mTOR complex 1 and enlarged neuronal size.

CONCLUSIONS AND RELEVANCE

In this study, mutations of MTOR were associated with a spectrum of brain overgrowth phenotypes extending from FCD type 2a to diffuse megalencephaly, distinguished by different mutations and levels of mosaicism. These mutations may be sufficient to cause cellular hypertrophy in cultured neurons and may provide a demonstration of the pattern of mosaicism in brain and substantiate the link between mosaic mutations of MTOR and pigmentary mosaicism in skin.

Improvement of ENU Mutagenesis Efficiency Using Serial Injection and Mismatch Repair Deficiency Mice

ENU mutagenesis is a powerful method for generating novel lines of mice that are informative with respect to both fundamental biological processes and human disease. Rapid developments in genomic technology have made the task of identifying causal mutations by positional cloning remarkably efficient. One limitation of this approach remains the mutation frequency achievable using standard treatment protocols, which currently generate approximately 1–2 sequence changes per megabase when optimized. In this study we used two strategies to attempt to increase the number of mutations induced by ENU treatment. One approach employed mice carrying a mutation in the DNA repair enzyme Msh6. The second strategy involved injection of ENU to successive generations of mice. To evaluate the number of ENU-induced mutations, single mice or pooled samples were analyzed using whole exome sequencing. The results showed that there is considerable variability in the induced mutation frequency using these approaches, but an overall increase in ENU-induced variants from one generation to another was observed. The analysis of the mice deficient for Msh6 also showed an increase in the ENU-induced variants compared to the wild-type ENU-treated mice. However, in both cases the increase in ENU-induced mutation frequency was modest.

A novel rasopathy caused by recurrent de novo missense mutations in PPP1CB closely resembles Noonan syndrome with loose anagen hair

Noonan syndrome is a rasopathy caused by mutations in multiple genes encoding components of the RAS/MAPK pathway. Despite its variable phenotype, limited genotype-phenotype correlations exist. Noonan syndrome with loose anagen hair (NS-LAH) is characterized by its distinctive hair anomalies, developmental differences, and structural brain abnormalities and is caused by a single recurrent missense SHOC2 mutation. SHOC2 forms a complex with protein phosphatase 1 (PP1C). Protein phosphatases counterbalance kinases and control activation of signaling proteins, such as the mitogen-activated protein kinases of the RAS/MAPK pathway. Here we report four patients with de novo missense mutations in protein phosphatase one catalytic subunit beta (PPP1CB), sharing a recognizable phenotype. Three individuals had the recurrent PPP1CB c.146G>C, p.Pro49Arg mutation, the fourth had a c.166G>C, p.Ala56Pro change. All had relative or absolute macrocephaly, low-set and posteriorly angulated ears, and developmental delay. Slow growing and/or sparse hair and/or an unruly hair texture was present in all. Three individuals had feeding difficulties requiring feeding tubes. One of two males had cryptorchidism, another had pectus excavatum. Short stature was present in three. A female with the recurrent mutation had a Dandy-Walker malformation and optic nerve hypoplasia. Mild ventriculomegaly occurred in all, cerebellar tonsillar ectopia was seen in two and progressed to Chiari 1 malformation in one individual. Based on the combination of phenotypic findings and PPP1CB's effect on RAF dephosphorylation within the RAS/MAPK pathway, this novel condition can be considered a rasopathy, most similar to NS-LAH. Collectively, these mutations meet the standardized criteria for pathogenicity. © 2016 Wiley Periodicals, Inc.

Update on the ACTG1-associated Baraitser-Winter cerebrofrontofacial syndrome

Baraitser-Winter cerebrofrontofacial syndrome is caused by heterozygous missense mutations in one of the two ubiquitous cytoplasmic actin-encoding genes ACTB and ACTG1. Recently, we characterized the large cohort of 41 patients presenting with this condition. Our series contained 34 patients with mutations in ACTB and only nine with ACTG1 mutations. Here, we report on seven unrelated patients with six mutations in ACTG1-four novel and two previously reported. Only one of seven patients was clinically diagnosed with this disorder and underwent ACTB/ACTG1 targeted sequencing, four patients were screened as a part of the large lissencephaly cohort and two were tested with exome sequencing. Retrospectively, facial features were compatible with the diagnosis but significantly milder than previously reported in four patients, and non-specific in one. The pattern of malformations of cortical development was highly similar in four of six patients with available MRI images and encompassed frontal predominant pachygyria merging with the posterior predominant band heterotopia. Two remaining patients showed mild involvement consistent with bilaterally simplified gyration over the frontal lobes. Taken together, we expand the clinical spectrum of the ACTG1-associated Baraitser-Winter cerebrofrontofacial syndrome demonstrating the mild end of the facial and brain manifestations. © 2016 Wiley Periodicals, Inc.

Mosaic Activating Mutations in FGFR1 Cause Encephalocraniocutaneous Lipomatosis

Encephalocraniocutaneous lipomatosis (ECCL) is a sporadic condition characterized by ocular, cutaneous, and central nervous system anomalies. Key clinical features include a well-demarcated hairless fatty nevus on the scalp, benign ocular tumors, and central nervous system lipomas. Seizures, spasticity, and intellectual disability can be present, although affected individuals without seizures and with normal intellect have also been reported. Given the patchy and asymmetric nature of the malformations, ECCL has been hypothesized to be due to a post-zygotic, mosaic mutation. Despite phenotypic overlap with several other disorders associated with mutations in the RAS-MAPK and PI3K-AKT pathways, the molecular etiology of ECCL remains unknown. Using exome sequencing of DNA from multiple affected tissues from five unrelated individuals with ECCL, we identified two mosaic mutations, c.1638C>A (p.Asn546Lys) and c.1966A>G (p.Lys656Glu) within the tyrosine kinase domain of FGFR1, in two affected individuals each. These two residues are the most commonly mutated residues in FGFR1 in human cancers and are associated primarily with CNS tumors. Targeted resequencing of FGFR1 in multiple tissues from an independent cohort of individuals with ECCL identified one additional individual with a c.1638C>A (p.Asn546Lys) mutation in FGFR1. Functional studies of ECCL fibroblast cell lines show increased levels of phosphorylated FGFRs and phosphorylated FRS2, a direct substrate of FGFR1, as well as constitutive activation of RAS-MAPK signaling. In addition to identifying the molecular etiology of ECCL, our results support the emerging overlap between mosaic developmental disorders and tumorigenesis.

Variant mapping and mutation discovery in inbred mice using next-generation sequencing

Background

The development of powerful new methods for DNA sequencing enable the discovery of sequence variants, their utilization for the mapping of mutant loci, and the identification of causal variants in a single step. We have applied this approach for the analysis of ENU-mutagenized mice maintained on an inbred background.

Results

We ascertained ENU-induced variants in four different phenotypically mutant lines. These were then used as informative markers for positional cloning of the mutated genes. We tested both whole genome (WGS) and whole exome (WES) datasets.

Conclusion

Both approaches were successful as a means to localize a region of homozygosity, as well as identifying mutations of candidate genes, which could be individually assessed. As expected, the WGS strategy was more reliable, since many more ENU-induced variants were ascertained.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2173-1) contains supplementary material, which is available to authorized users.

Characterisation of mutations of the phosphoinositide-3-kinase regulatory subunit, PIK3R2, in perisylvian polymicrogyria: a next-generation sequencing study

Background

Bilateral perisylvian polymicrogyria (BPP), the most common form of regional polymicrogyria, causes the congenital bilateral perisylvian syndrome, featuring oromotor dysfunction, cognitive impairment and epilepsy. BPP is etiologically heterogeneous, but only a few genetic causes have been reported. The aim of this study was to identify additional genetic etiologies of BPP and delineate their frequency in this patient population.

Methods

We performed child-parent (trio)-based whole exome sequencing (WES) on eight children with BPP. Following the identification of mosaic PIK3R2 mutations in two of these eight children, we performed targeted screening of PIK3R2 in a cohort of 118 children with BPP who were ascertained from 1980 until 2015 using two methods. First, we performed targeted sequencing of the entire PIK3R2 gene by single molecule molecular inversion probes (smMIPs) on 38 patients with BPP with normal-large head size. Second, we performed amplicon sequencing of the recurrent PIK3R2 mutation (p.Gly373Arg) on 80 children with various types of polymicrogyria including BPP. One additional patient underwent clinical WES independently, and was included in this study given the phenotypic similarity to our cohort. All patients included in this study were children (< 18 years of age) with polymicrogyria enrolled in our research program.

Findings

Using WES, we identified a mosaic mutation (p.Gly373Arg) in the regulatory subunit of the PI3K-AKT-MTOR pathway, PIK3R2, in two children with BPP. Of the 38 patients with BPP and normal-large head size who underwent targeted next generation sequencing by smMIPs, we identified constitutional and mosaic PIK3R2 mutations in 17 additional children. In parallel, one patient was found to have the recurrent PIK3R2 mutation by clinical WES. Seven patients had BPP alone, and 13 had BPP in association with features of the megalencephaly-polymicrogyria-polydactyly-hydrocephalus syndrome (MPPH). Nineteen patients had the same mutation (Gly373Arg), and one had a nearby missense mutation (p.Lys376Glu). Across the entire cohort, mutations were constitutional in 12 and mosaic in eight patients. Among mosaic patients, we observed substantial variation in alternate (mutant) allele levels ranging from 2·5% (10/377) to 36·7% (39/106) of reads, equivalent to 5–73·4% of cells analyzed. Levels of mosaicism varied from undetectable to 17·1% (37/216) of reads in blood-derived compared to 29·4% (2030/6889) to 43·3% (275/634) in saliva-derived DNA.

Interpretation

Constitutional and mosaic mutations in the PIK3R2 gene are associated with a spectrum of developmental brain disorders ranging from BPP with a normal head size to the megalencephaly-polymicrogyria-polydactyly-hydrocephalus syndrome. The phenotypic variability and low-level mosaicism challenging conventional molecular methods have important implications for genetic testing and counseling.

Osteoblast differentiation profiles define sex specific gene expression patterns in craniosynostosis

Single suture craniosynostosis (SSC) is the premature fusion of one calvarial suture and occurs in 1-1700-2500 live births. Congenital fusion of either the sagittal, metopic, or coronal sutures represents 95% of all cases of SSC. Sagittal and metopic synostosis have a male preponderance (3:1) while premature fusion of the coronal suture has a female preponderance (2:1). Although environmental and genetic factors contribute to SSC, the etiology of the majority of SSC cases remains unclear. In this study, 227 primary calvarial osteoblast cell lines from patients with coronal, metopic, or sagittal synostosis and unaffected controls were established and assayed for ALP activity and BrdU incorporation (n = 226) as respective measures of early stage osteoblast differentiation and proliferation. Primary osteoblast cell lines from individuals with sagittal synostosis demonstrated higher levels of ALP activity and reduced proliferation when compared to control lines. In order to address the sex differences in SSC types, the data was further stratified by sex. Osteoblasts from males and females with sagittal synostosis as well as males with metopic synostosis demonstrated higher levels of ALP activity when compared to sex matched controls, and males with sagittal or metopic synostosis demonstrated reduced levels of proliferation. In order to elucidate genes and pathways involved in these observed phenotypes, correlation analyses comparing ALP activity and proliferation to global gene expression was performed. Transcripts related to osteoblast differentiation were identified both differentially up and downregulated, correlated with ALP activity when compared to controls, and demonstrated a striking sex specific gene expression pattern. These data support that the dysregulation of osteoblast differentiation plays a role in the development of SSC and that genetic factors contribute to the observed sex related differences.

Heritable GATA2 mutations associated with familial AML-MDS: a case report and review of literature

A 50-year-old woman was diagnosed with acute myeloid leukemia (AML). She has history of thrombocytopenia for 25 years and a significant family history of thrombocytopenia, affecting her mother, siblings and their children, as well as her own children. Both her mother and maternal aunt died from myelodysplastic syndrome (MDS). Additional genetic analysis was performed and identified two heterozygous missence mutations in the second zinc finger domain of GATA2 gene (p.Thr358Lys, and p.Leu359Val), occurring in cis on the same allele. Given the patient’s family history and clinical manifestation, this was interpreted as an acute myeloid leukemia with heritable GATA2 mutations associated with familial AML-MDS. Germline GATA2 mutations are involved in a group of complex syndromes with overlapping clinical features of immune deficiency, lymphedema and propensity to acute myeloid leukemia or myelodysplastic syndrome (AML-MDS). Here we reported a case of familial AML-MDS with two novel GATA2 mutations. This case illustrates the importance of recognizing the clinical features for this rare category of AML-MDS and performing the appropriate molecular testing. The diagnosis of heritable gene mutations associated familial AML-MDS has significant clinical implication for the patients and affected families. Clinical trials are available to further investigate the role of allogeneic hematopoietic stem cell transplant in managing these patients.

Support for the N-methyl-D-aspartate receptor hypofunction hypothesis of schizophrenia from exome sequencing in multiplex families

IMPORTANCE

Schizophrenia is a complex genetic disorder demonstrating considerable heritability. Genetic studies have implicated many different genes and pathways, but much of the genetic liability remains unaccounted for. Investigation of genetic forms of schizophrenia will lead to a better understanding of the underlying molecular pathways, which will then enable targeted approaches for disease prevention and treatment.

OBJECTIVE

To identify new genetic factors strongly predisposing to schizophrenia in families with multiple affected individuals with schizophrenia.

DESIGN

We performed genome-wide array comparative genomic hybridization, linkage analysis, and exome sequencing in multiplex families with schizophrenia.

SETTING

Probands and their family members were recruited from academic medical centers.

PARTICIPANTS

We intended to identify rare disease-causing mutations in 5 large families where schizophrenia transmission appears consistent with single-gene inheritance.

INTERVENTION

Array comparative genomic hybridization was used to identify copy number variants, while exome sequencing was used to identify variants shared in all affected individuals and linkage analysis was used to further filter shared variants of interest. Analysis of select variants was performed in cultured cells to assess their functional consequences.

MAIN OUTCOME MEASURES

Rare inherited disease-related genetic mutations.

RESULTS

No segregating rare copy number variants were detected by array comparative genomic hybridization. However, in all 5 families, exome sequencing detected rare protein-altering variants in 1 of 3 genes associated with the N -methyl-D-aspartate (NMDA) receptor. One pedigree shared a missense and frameshift substitution of GRM5, encoding the metabotropic glutamate receptor subtype 5 (mGluR5), which is coupled to the NMDA receptor and potentiates its signaling; the frameshift disrupts binding to the scaffolding protein tamalin and increases mGluR5 internalization. Another pedigree transmitted a missense substitution in PPEF2, encoding a calmodulin-binding protein phosphatase, which we show influences mGluR5 levels. Three pedigrees demonstrated different missense substitutions within LRP1B, encoding a low-density lipoprotein receptor-related protein tied to both the NMDA receptor and located in a chromosome 2q22 region previously strongly linked to schizophrenia.

CONCLUSIONS AND RELEVANCE

Exome sequencing of multiplex pedigrees uncovers new genes associated with risk for developing schizophrenia and suggests potential novel therapeutic targets.

Heritable GATA2 mutations associated with familial myelodysplastic syndrome and acute myeloid leukemia

We report the discovery of the GATA2 gene as a new myelodysplastic syndrome (MDS)/acute myeloid leukemia (AML) predisposition gene. We found the same, novel heterozygous c.1061C>T (p.Thr354Met) missense mutation in the GATA2 transcription factor gene segregating with the multigenerational transmission of MDS/AML in three families, and a GATA2 c.1063_1065delACA (p.Thr355del) mutation at an adjacent codon in a fourth MDS/AML family. The mutations reside within the second zinc finger of GATA2 which mediates DNA-binding and protein-protein interactions. We show differential effects of the mutants on transactivation of target genes, cellular differentiation, apoptosis and global gene expression. Identification of such predisposing genes to familial forms of MDS and AML is critical for more effective diagnosis and prognosis, counselling, selection of related bone marrow transplant donors, and development of therapies.

KLHDC8B in Hodgkin lymphoma and possibly twinning

A key feature of Hodgkin lymphoma is that the malignant cells are binucleated, as a consequence of failed cytokinesis. We recently ascertained a family in which multiple cases of Hodgkin lymphoma had occurred among individuals who inherited a balanced chromosomal translocation. We cloned the translocation breakpoints and found that it disrupted a previously uncharacterized gene, KLHDC8B, encoding a Kelch family protein whose deficiency impairs cytokinesis and leads to binucleated cells. In other families we found a rare single nucleotide polymorphism affecting mitotic translation of KLHDC8B that was associated with and linked to Hodgkin lymphoma. Interestingly, the index family demonstrated an unusual frequency of twins, and there is a previously reported association between Hodgkin lymphoma and twins. Here we review the unusual genetic features of Hodgkin lymphoma, including gender concordance among siblings, and genetically test the hypothesis that KLHDC8B may participate in twinning by disrupting cytokinesis through impediment of polar body separation from oocytes.

Combined analysis of three whole genome linkage scans for Ankylosing Spondylitis

OBJECTIVE

Ankylosing spondylitis (AS) is a debilitating chronic inflammatory condition with a high degree of familiality (lambda(s) = 82) and heritability (>90%) that primarily affects spinal and sacroiliac joints. Whole genome scans for linkage to AS phenotypes have been conducted, although results have been inconsistent between studies and all have had modest sample sizes. One potential solution to these issues is to combine data from multiple studies in a retrospective meta-analysis.

METHODS

The International Genetics of Ankylosing Spondylitis Consortium combined data from three whole genome linkage scans for AS (n = 3744 subjects) to determine chromosomal markers that show evidence of linkage with disease. Linkage markers typed in different centres were integrated into a consensus map to facilitate effective data pooling. We performed a weighted meta-analysis to combine the linkage results, and compared them with the three individual scans and a combined pooled scan.

RESULTS

In addition to the expected region surrounding the HLA-B27 gene on chromosome 6, we determined that several marker regions showed significant evidence of linkage with disease status. Regions on chromosome 10q and 16q achieved 'suggestive' evidence of linkage, and regions on chromosomes 1q, 3q, 5q, 6q, 9q, 17q and 19q showed at least nominal linkage in two or more scans and in the weighted meta-analysis. Regions previously associated with AS on chromosome 2q (the IL-1 gene cluster) and 22q (CYP2D6) exhibited nominal linkage in the meta-analysis, providing further statistical support for their involvement in susceptibility to AS.

CONCLUSION

These findings provide a useful guide for future studies aiming to identify the genes involved in this highly heritable condition.

Replication of association of IL1 gene complex members with ankylosing spondylitis in Taiwanese Chinese

OBJECTIVE

To test the association of interleukin 1 (IL1) gene family members with ankylosing spondylitis (AS), previously reported in Europid subjects, in an ethnically remote population.

METHODS

200 Taiwanese Chinese AS patients and 200 ethnically matched healthy controls were genotyped for five single nucleotide polymorphisms (SNPs) and the IL1RN.VNTR, markers previously associated with AS. Allele, genotype, and haplotype frequencies were compared between cases and controls.

RESULTS

Association of alleles and genotypes of the markers IL1F10.3, IL1RN.4, and IL1RN.VNTR was observed with AS (p<0.05). Haplotypes of pairs of these markers and of the markers IL1RN.6/1 and IL1RN.6/2 were also significantly associated with AS. The strongest associations observed were with the marker IL1RN.4, and with the two-marker haplotype IL1RN.4-IL1RN.VNTR (both p = 0.004). Strong linkage disequilibrium was observed between all marker pairs except those involving IL1B-511 (D' 0.4 to 0.9, p<0.01).

CONCLUSIONS

The IL1 gene cluster is associated with AS in Taiwanese Chinese. This finding provides strong statistical support that the previously observed association of this gene cluster with AS is a true positive finding. These authors contributed equally to the study.

The interleukin 1 gene cluster contains a major susceptibility locus for ankylosing spondylitis

Ankylosing spondylitis (AS) is a common and highly heritable inflammatory arthropathy. Although the gene HLA-B27 is almost essential for the inheritance of the condition, it alone is not sufficient to explain the pattern of familial recurrence of the disease. We have previously demonstrated suggestive linkage of AS to chromosome 2q13, a region containing the interleukin 1 (IL-1) family gene cluster, which includes several strong candidates for involvement in the disease. In the current study, we describe strong association and transmission of IL-1 family gene cluster single-nucleotide polymorphisms and haplotypes with AS.

Dissection of class III major histocompatibility complex haplotypes associated with rheumatoid arthritis

OBJECTIVE

We have previously identified a single-nucleotide polymorphism (SNP) haplotype involving the lymphotoxin alpha (LTA) and tumor necrosis factor (TNF) loci (termed haplotype LTA-TNF2) on chromosome 6 that shows differential association with rheumatoid arthritis (RA) on HLA-DRB1*0404 and *0401 haplotypes, suggesting the presence of additional non-HLA-DRB1 RA susceptibility genes on these haplotypes. To refine this association, we performed a case-control association study using both SNPs and microsatellite markers in haplotypes matched either for HLA-DRB1*0404 or for HLA-DRB1*0401.

METHODS

Fourteen SNPs lying between HLA-DRB1 and LTA were genotyped in 87 DRB1*04-positive families. High-density microsatellite typing was performed using 24 markers spanning 2,500 kb centered around the TNF gene in 305 DRB1*0401 or *0404 cases and 400 DRB1*0401 or *0404 controls. Single-marker, 2-marker, and 3-marker minihaplotypes were constructed and their frequencies compared between the DRB1*0401 and DRB1*0404 matched case and control haplotypes.

RESULTS

Marked preservation of major histocompatibility complex haplotypes was seen, with chromosomes carrying LTA-TNF2 and either DRB1*0401 or DRB1*0404 both carrying an identical SNP haplotype across the 1-Mb region between TNF and HLA-DRB1. Using microsatellite markers, we observed two 3-marker minihaplotypes that were significantly overrepresented in the DRB1*0404 case haplotypes (P = 0.00024 and P = 0.00097).

CONCLUSION

The presence of a single extended SNP haplotype between LTA-TNF2 and both DRB1*0401 and DRB1*0404 is evidence against this region harboring the genetic effects in linkage disequilibrium with LTA-TNF2. Two RA-associated haplotypes on the background of DRB1*0404 were identified in a 126-kb region surrounding and centromeric to the TNF locus.

Investigation of the role of ANKH in ankylosing spondylitis

OBJECTIVE

The ank/ank mouse develops a phenotype similar to ankylosing spondylitis (AS) in humans. ANKH, the human homolog of the mutated gene in the ank/ank mouse, has been implicated in familial autosomal-dominant chondrocalcinosis and autosomal-dominant craniometaphyseal dysplasia. This study was undertaken to investigate the role of ANKH in susceptibility to and clinical manifestations of AS.

METHODS

Sequence variants were identified by genomic sequencing of the 12 ANKH exons and their flanking splice sites in 48 AS patients; variants were then screened in 233 patients and 478 controls. Linkage to the ANKH locus was assessed in 185 affected-sibling-pair families.

RESULTS

Five single-nucleotide polymorphisms were identified within the coding region and flanking splice sites. No association between either susceptibility to AS or its clinical manifestations and these novel polymorphisms, or between disease susceptibility and 3 known promoter variants, was seen. No linkage between the ANKH locus and AS was observed. Multipoint exclusion mapping rejected the hypothesis of a locus of a magnitude lambda>/=1.4 (logarithm of odds score <-2) (equivalent to a genetic contribution of >10% to the AS sibling recurrence risk ratio) within this area contributing to AS.

CONCLUSION

These findings indicate that ANKH is not significantly involved in susceptibility to or clinical manifestations of AS.

Genetic testing for haemochromatosis in patients with chondrocalcinosis

Hereditary haemochromatosis (HH) is the most common lethal monogenic human disease, affecting roughly 1 in 300 white northern Europeans. Homozygosity for the C282Y polymorphism within the HFE gene causes more than 80% of cases, with compound heterozygosity of the C282Y and H63D polymorphism also increasing susceptibility to disease. The aim of this study was to determine the frequency of the C282Y and H63D polymorphisms in the disease, and to assess the risk of HH in heterozygotes for the C282Y polymorphism. 128 patients were recruited because of either radiographic chondrocalcinosis (at least bicompartmental knee disease or joints other than the knee involved) or CPPD pseudogout. Genotyping of the HFE C282Y and H63D mutations was performed using PCR/SSP and genotypes for the C282Y polymorphism confirmed by PCR/RFLP. Historical white European control data were used for comparison. Two previously undiagnosed C282Y homozygotes (1.6%), and 16 C282Y heterozygotes (12.5%), including four (3.1%) C282Y/H63D compound heterozygotes were identified. This represents a significant overrepresentation of C282Y homozygotes (relative risk 3.4, p=0.037), but the number of heterozygotes was not significantly increased. At a cost per test of pound1 for each subject, screening all patients with chondrocalcinosis using the above ascertainment criteria costs only pound64 for each case of haemochromatosis identified, clearly a highly cost effective test given the early mortality associated with untreated haemochromatosis. Routine screening for haemochromatosis in patients with appreciable chondrocalcinosis is recommended.

A search for type 1 diabetes susceptibility genes in families from the United Kingdom

Genetic analysis of a mouse model of major histocompatability complex (MHC)-associated autoimmune type 1 (insulin-dependent) diabetes mellitus (IDDM) has shown that the disease is caused by a combination of a major effect at the MHC and at least ten other susceptibility loci elsewhere in the genome. A genome-wide scan of 93 affected sibpair families (ASP) from the UK (UK93) indicated a similar genetic basis for human type 1 diabetes, with the major genetic component at the MHC locus (IDDM1) explaining 34% of the familial clustering of the disease (lambda(s)=2.5; refs 3,4). In the present report, we have analysed a further 263 multiplex families from the same population (UK263) to provide a total UK data set of 356 ASP families (UK356). Only four regions of the genome outside IDDM1/MHC, which was still the only major locus detected, were not excluded at lambda(s)=3 and lod=-2, of which two showed evidence of linkage: chromosome 10p13-p11 (maximum lod score (MLS)=4.7, P=3x10(-6), lambda(s)=1.56) and chromosome 16q22-16q24 (MLS=3.4, P=6.5x10(-5), lambda(s)=1.6). These and other novel regions, including chromosome 14q12-q21 and chromosome 19p13-19q13, could potentially harbour disease loci but confirmation and fine mapping cannot be pursued effectively using conventional linkage analysis. Instead, more powerful linkage disequilibrium-based and haplotype mapping approaches must be used; such data is already emerging for several type 1 diabetes loci detected initially by linkage.