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 Ly6Chi 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.

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.

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.

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.

Transcriptomic Signatures of Single-Suture Craniosynostosis Phenotypes.

Craniosynostosis is a birth defect where calvarial sutures close prematurely, as part of a genetic syndrome or independently, with unknown cause. This study aimed to identify differences in gene expression in primary calvarial cell lines derived from patients with four phenotypes of single-suture craniosynostosis, compared to controls. Calvarial bone samples (N = 388 cases/85 controls) were collected from clinical sites during reconstructive skull surgery. Primary cell lines were then derived from the tissue and used for RNA sequencing. Linear models were fit to estimate covariate adjusted associations between gene expression and four phenotypes of single-suture craniosynostosis (lambdoid, metopic, sagittal, and coronal), compared to controls. Sex-stratified analysis was also performed for each phenotype. Differentially expressed genes (DEGs) included 72 genes associated with coronal, 90 genes associated with sagittal, 103 genes associated with metopic, and 33 genes associated with lambdoid craniosynostosis. The sex-stratified analysis revealed more DEGs in males (98) than females (4). There were 16 DEGs that were homeobox (HOX) genes. Three TFs (SUZ12, EZH2, AR) significantly regulated expression of DEGs in one or more phenotypes. Pathway analysis identified four KEGG pathways associated with at least one phenotype of craniosynostosis. Together, this work suggests unique molecular mechanisms related to craniosynostosis phenotype and fetal sex.

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.

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.

Delineating the molecular and phenotypic spectrum of the SETD1B-related syndrome.

PURPOSE: Pathogenic variants in SETD1B have been associated with a syndromic neurodevelopmental disorder including intellectual disability, language delay, and seizures. To date, clinical features have been described for 11 patients with (likely) pathogenic SETD1B sequence variants. This study aims to further delineate the spectrum of the SETD1B-related syndrome based on characterizing an expanded patient cohort. METHODS: We perform an in-depth clinical characterization of a cohort of 36 unpublished individuals with SETD1B sequence variants, describing their molecular and phenotypic spectrum. Selected variants were functionally tested using in vitro and genome-wide methylation assays. RESULTS: Our data present evidence for a loss-of-function mechanism of SETD1B variants, resulting in a core clinical phenotype of global developmental delay, language delay including regression, intellectual disability, autism and other behavioral issues, and variable epilepsy phenotypes. Developmental delay appeared to precede seizure onset, suggesting SETD1B dysfunction impacts physiological neurodevelopment even in the absence of epileptic activity. Males are significantly overrepresented and more severely affected, and we speculate that sex-linked traits could affect susceptibility to penetrance and the clinical spectrum of SETD1B variants. CONCLUSION: Insights from this extensive cohort will facilitate the counseling regarding the molecular and phenotypic landscape of newly diagnosed patients with the SETD1B-related syndrome.

Evidence of disrupted rhombic lip development in the pathogenesis of Dandy-Walker malformation.

Dandy-Walker malformation (DWM) and Cerebellar vermis hypoplasia (CVH) are commonly recognized human cerebellar malformations diagnosed following ultrasound and antenatal or postnatal MRI. Specific radiological criteria are used to distinguish them, yet little is known about their differential developmental disease mechanisms. We acquired prenatal cases diagnosed as DWM and CVH and studied cerebellar morphobiometry followed by histological and immunohistochemical analyses. This was supplemented by laser capture microdissection and RNA-sequencing of the cerebellar rhombic lip, a transient progenitor zone, to assess the altered transcriptome of DWM vs control samples. Our radiological findings confirm that the cases studied fall within the accepted biometric range of DWM. Our histopathological analysis points to reduced foliation and inferior vermian hypoplasia as common features in all examined DWM cases. We also find that the rhombic lip, a dorsal stem cell zone that drives the growth and maintenance of the posterior vermis is specifically disrupted in DWM, with reduced proliferation and self-renewal of the progenitor pool, and altered vasculature, all confirmed by transcriptomics analysis. We propose a unified model for the developmental pathogenesis of DWM. We hypothesize that rhombic lip development is disrupted through either aberrant vascularization and/or direct insult which causes reduced proliferation and failed expansion of the rhombic lip progenitor pool leading to disproportionate hypoplasia and dysplasia of the inferior vermis. Timing of insult to the developing rhombic lip (before or after 14 PCW) dictates the extent of hypoplasia and distinguishes DWM from CVH.

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.

Analysis of 17 genes detects mutations in 81% of 811 patients with lissencephaly.

PURPOSE: To estimate diagnostic yield and genotype-phenotype correlations in a cohort of 811 patients with lissencephaly or subcortical band heterotopia. METHODS: We collected DNA from 756 children with lissencephaly over 30 years. Many were tested for deletion 17p13.3 and mutations of LIS1, DCX, and ARX, but few other genes. Among those tested, 216 remained unsolved and were tested by a targeted panel of 17 genes (ACTB, ACTG1, ARX, CRADD, DCX, LIS1, TUBA1A, TUBA8, TUBB2B, TUBB, TUBB3, TUBG1, KIF2A, KIF5C, DYNC1H1, RELN, and VLDLR) or by whole-exome sequencing. Fifty-five patients studied at another institution were added as a validation cohort. RESULTS: The overall mutation frequency in the entire cohort was 81%. LIS1 accounted for 40% of patients, followed by DCX (23%), TUBA1A (5%), and DYNC1H1 (3%). Other genes accounted for 1% or less of patients. Nineteen percent remained unsolved, which suggests that several additional genes remain to be discovered. The majority of unsolved patients had posterior pachygyria, subcortical band heterotopia, or mild frontal pachygyria. CONCLUSION: The brain-imaging pattern correlates with mutations in single lissencephaly-associated genes, as well as in biological pathways. We propose the first LIS classification system based on the underlying molecular mechanisms.

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.

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.

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.

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.

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.