oFlowSeq: a quantitative approach to identify protein coding mutations affecting cell type enrichment using mosaic CRISPR-Cas9 edited cerebral organoids.

Cerebral organoids are comprised of diverse cell types found in the developing human brain, and can be leveraged in the identification of critical cell types perturbed by genetic risk variants in common, neuropsychiatric disorders. There is great interest in developing high-throughput technologies to associate genetic variants with cell types. Here, we describe a high-throughput, quantitative approach (oFlowSeq) by utilizing CRISPR-Cas9, FACS sorting, and next-generation sequencing. Using oFlowSeq, we found that deleterious mutations in autism-associated gene KCTD13 resulted in increased proportions of Nestin+ cells and decreased proportions of TRA-1-60+ cells within mosaic cerebral organoids. We further identified that a locus-wide CRISPR-Cas9 survey of another 18 genes in the 16p11.2 locus resulted in most genes with > 2% maximum editing efficiencies for short and long indels, suggesting a high feasibility for an unbiased, locus-wide experiment using oFlowSeq. Our approach presents a novel method to identify genotype-to-cell type imbalances in an unbiased, high-throughput, quantitative manner.

The Genetic Landscape of Diamond-Blackfan Anemia.

Diamond-Blackfan anemia (DBA) is a rare bone marrow failure disorder that affects 7 out of 1,000,000 live births and has been associated with mutations in components of the ribosome. In order to characterize the genetic landscape of this heterogeneous disorder, we recruited a cohort of 472 individuals with a clinical diagnosis of DBA and performed whole-exome sequencing (WES). We identified relevant rare and predicted damaging mutations for 78% of individuals. The majority of mutations were singletons, absent from population databases, predicted to cause loss of function, and located in 1 of 19 previously reported ribosomal protein (RP)-encoding genes. Using exon coverage estimates, we identified and validated 31 deletions in RP genes. We also observed an enrichment for extended splice site mutations and validated their diverse effects using RNA sequencing in cell lines obtained from individuals with DBA. Leveraging the size of our cohort, we observed robust genotype-phenotype associations with congenital abnormalities and treatment outcomes. We further identified rare mutations in seven previously unreported RP genes that may cause DBA, as well as several distinct disorders that appear to phenocopy DBA, including nine individuals with biallelic CECR1 mutations that result in deficiency of ADA2. However, no new genes were identified at exome-wide significance, suggesting that there are no unidentified genes containing mutations readily identified by WES that explain >5% of DBA-affected case subjects. Overall, this report should inform not only clinical practice for DBA-affected individuals, but also the design and analysis of rare variant studies for heterogeneous Mendelian disorders.

An enhanced CRISPR repressor for targeted mammalian gene regulation.

The RNA-guided endonuclease Cas9 can be converted into a programmable transcriptional repressor, but inefficiencies in target-gene silencing have limited its utility. Here we describe an improved Cas9 repressor based on the C-terminal fusion of a rationally designed bipartite repressor domain, KRAB-MeCP2, to nuclease-dead Cas9. We demonstrate the system's superiority in silencing coding and noncoding genes, simultaneously repressing a series of target genes, improving the results of single and dual guide RNA library screens, and enabling new architectures of synthetic genetic circuits.

Genome-wide Analysis of Body Proportion Classifies Height-Associated Variants by Mechanism of Action and Implicates Genes Important for Skeletal Development.

Human height is a composite measurement, reflecting the sum of leg, spine, and head lengths. Many common variants influence total height, but the effects of these or other variants on the components of height (body proportion) remain largely unknown. We studied sitting height ratio (SHR), the ratio of sitting height to total height, to identify such effects in 3,545 African Americans and 21,590 individuals of European ancestry. We found that SHR is heritable: 26% and 39% of the total variance of SHR can be explained by common variants in European and African Americans, respectively, and global European admixture is negatively correlated with SHR in African Americans (r(2) ≈ 0.03). Six regions reached genome-wide significance (p < 5 × 10(-8)) for association with SHR and overlapped biological candidate genes, including TBX2 and IGFBP3. We found that 130 of 670 height-associated variants are nominally associated (p < 0.05) with SHR, more than expected by chance (p = 5 × 10(-40)). At these 130 loci, the height-increasing alleles are associated with either a decrease (71 loci) or increase (59 loci) in SHR, suggesting that different height loci disproportionally affect either leg length or spine/head length. Pathway analyses via DEPICT revealed that height loci affecting SHR, and especially those affecting leg length, show enrichment of different biological pathways (e.g., bone/cartilage/growth plate pathways) than do loci with no effect on SHR (e.g., embryonic development). These results highlight the value of using a pair of related but orthogonal phenotypes, in this case SHR with height, as a prism to dissect the biology underlying genetic associations in polygenic traits and diseases.

An excess of risk-increasing low-frequency variants can be a signal of polygenic inheritance in complex diseases.

In most complex diseases, much of the heritability remains unaccounted for by common variants. It has been postulated that lower-frequency variants contribute to the remaining heritability. Here, we describe a method to test for polygenic inheritance from lower-frequency variants by using GWAS summary association statistics. We explored scenarios with many causal low-frequency variants and showed that there is more power to detect risk variants than to detect protective variants, resulting in an increase in the ratio of detected risk to protective variants (R/P ratio). Such an excess can also occur if risk variants are present and kept at lower frequencies because of negative selection. The R/P ratio can be falsely elevated because of reasons unrelated to polygenic inheritance, such as uneven sample sizes or asymmetric population stratification, so precautions to correct for these confounders are essential. We tested our method on published GWAS results and observed a strong signal in some diseases (schizophrenia and type 2 diabetes) but not others. We also explored the shared genetic component in overlapping phenotypes related to inflammatory bowel disease (Crohn disease [CD] and ulcerative colitis [UC]) and diabetic nephropathy (macroalbuminuria and end-stage renal disease [ESRD]). Although the signal was still present when both CD and UC were jointly analyzed, the signal was lost when macroalbuminuria and ESRD were jointly analyzed, suggesting that these phenotypes should best be studied separately. Thus, our method may also help guide the design of future genetic studies of various traits and diseases.

A novel test for recessive contributions to complex diseases implicates Bardet-Biedl syndrome gene BBS10 in idiopathic type 2 diabetes and obesity.

Rare-variant association studies in common, complex diseases are customarily conducted under an additive risk model in both single-variant and burden testing. Here, we describe a method to improve detection of rare recessive variants in complex diseases termed RAFT (recessive-allele-frequency-based test). We found that RAFT outperforms existing approaches when the variant influences disease risk in a recessive manner on simulated data. We then applied our method to 1,791 Finnish individuals with type 2 diabetes (T2D) and 2,657 matched control subjects. In BBS10, we discovered a rare variant (c.1189A>G [p.Ile397Val]; rs202042386) that confers risk of T2D in a recessive state (p = 1.38 × 10(-6)) and would be missed by conventional methods. Testing of this variant in an established in vivo zebrafish model confirmed the variant to be pathogenic. Taken together, these data suggest that RAFT can effectively reveal rare recessive contributions to complex diseases overlooked by conventional association tests.

Distribution and medical impact of loss-of-function variants in the Finnish founder population.

Exome sequencing studies in complex diseases are challenged by the allelic heterogeneity, large number and modest effect sizes of associated variants on disease risk and the presence of large numbers of neutral variants, even in phenotypically relevant genes. Isolated populations with recent bottlenecks offer advantages for studying rare variants in complex diseases as they have deleterious variants that are present at higher frequencies as well as a substantial reduction in rare neutral variation. To explore the potential of the Finnish founder population for studying low-frequency (0.5-5%) variants in complex diseases, we compared exome sequence data on 3,000 Finns to the same number of non-Finnish Europeans and discovered that, despite having fewer variable sites overall, the average Finn has more low-frequency loss-of-function variants and complete gene knockouts. We then used several well-characterized Finnish population cohorts to study the phenotypic effects of 83 enriched loss-of-function variants across 60 phenotypes in 36,262 Finns. Using a deep set of quantitative traits collected on these cohorts, we show 5 associations (p<5×10⁻8) including splice variants in LPA that lowered plasma lipoprotein(a) levels (P = 1.5×10⁻¹¹7). Through accessing the national medical records of these participants, we evaluate the LPA finding via Mendelian randomization and confirm that these splice variants confer protection from cardiovascular disease (OR = 0.84, P = 3×10⁻4), demonstrating for the first time the correlation between very low levels of LPA in humans with potential therapeutic implications for cardiovascular diseases. More generally, this study articulates substantial advantages for studying the role of rare variation in complex phenotypes in founder populations like the Finns and by combining a unique population genetic history with data from large population cohorts and centralized research access to National Health Registers.

Ataxia, dementia, and hypogonadotropism caused by disordered ubiquitination.

BACKGROUND: The combination of ataxia and hypogonadism was first described more than a century ago, but its genetic basis has remained elusive. METHODS: We performed whole-exome sequencing in a patient with ataxia and hypogonadotropic hypogonadism, followed by targeted sequencing of candidate genes in similarly affected patients. Neurologic and reproductive endocrine phenotypes were characterized in detail. The effects of sequence variants and the presence of an epistatic interaction were tested in a zebrafish model. RESULTS: Digenic homozygous mutations in RNF216 and OTUD4, which encode a ubiquitin E3 ligase and a deubiquitinase, respectively, were found in three affected siblings in a consanguineous family. Additional screening identified compound heterozygous truncating mutations in RNF216 in an unrelated patient and single heterozygous deleterious mutations in four other patients. Knockdown of rnf216 or otud4 in zebrafish embryos induced defects in the eye, optic tectum, and cerebellum; combinatorial suppression of both genes exacerbated these phenotypes, which were rescued by nonmutant, but not mutant, human RNF216 or OTUD4 messenger RNA. All patients had progressive ataxia and dementia. Neuronal loss was observed in cerebellar pathways and the hippocampus; surviving hippocampal neurons contained ubiquitin-immunoreactive intranuclear inclusions. Defects were detected at the hypothalamic and pituitary levels of the reproductive endocrine axis. CONCLUSIONS: The syndrome of hypogonadotropic hypogonadism, ataxia, and dementia can be caused by inactivating mutations in RNF216 or by the combination of mutations in RNF216 and OTUD4. These findings link disordered ubiquitination to neurodegeneration and reproductive dysfunction and highlight the power of whole-exome sequencing in combination with functional studies to unveil genetic interactions that cause disease. (Funded by the National Institutes of Health and others.).

A pharmacogenetic study implicates SLC9a9 in multiple sclerosis disease activity.

OBJECTIVE: A proportion of multiple sclerosis (MS) patients experience disease activity despite treatment. The early identification of the most effective drug is critical to impact long-term outcome and to move toward a personalized approach. The aim of the present study is to identify biomarkers for further clinical development and to yield insights into the pathophysiology of disease activity. METHODS: We performed a genome-wide association study in interferon-ß (IFNß)-treated MS patients followed by validation in 3 independent cohorts. The role of the validated variant was examined in several RNA data sets, and the function of the presumed target gene was explored using an RNA interference approach in primary T cells in vitro. RESULTS: We found an association between rs9828519(G) and nonresponse to IFNß (pdiscovery = 4.43 × 10(-8)) and confirmed it in a meta-analysis across 3 replication data sets (preplication = 7.78 × 10(-4)). Only 1 gene is found in the linkage disequilibrium block containing rs9828519: SLC9A9. Exploring the function of this gene, we see that SLC9A9 mRNA expression is diminished in MS subjects who are more likely to have relapses. Moreover, SLC9A9 knockdown in T cells in vitro leads an increase in expression of IFNγ, which is a proinflammatory molecule. INTERPRETATION: This study identifies and validates the role of rs9828519, an intronic variant in SLC9A9, in IFNß-treated subjects, demonstrating a successful pharmacogenetic screen in MS. Functional characterization suggests that SLC9A9, an Na(+) -H(+) exchanger found in endosomes, appears to influence the differentiation of T cells to a proinflammatory fate and may have a broader role in MS disease activity, outside of IFNß treatment.

Integrated model of de novo and inherited genetic variants yields greater power to identify risk genes.

De novo mutations affect risk for many diseases and disorders, especially those with early-onset. An example is autism spectrum disorders (ASD). Four recent whole-exome sequencing (WES) studies of ASD families revealed a handful of novel risk genes, based on independent de novo loss-of-function (LoF) mutations falling in the same gene, and found that de novo LoF mutations occurred at a twofold higher rate than expected by chance. However successful these studies were, they used only a small fraction of the data, excluding other types of de novo mutations and inherited rare variants. Moreover, such analyses cannot readily incorporate data from case-control studies. An important research challenge in gene discovery, therefore, is to develop statistical methods that accommodate a broader class of rare variation. We develop methods that can incorporate WES data regarding de novo mutations, inherited variants present, and variants identified within cases and controls. TADA, for Transmission And De novo Association, integrates these data by a gene-based likelihood model involving parameters for allele frequencies and gene-specific penetrances. Inference is based on a Hierarchical Bayes strategy that borrows information across all genes to infer parameters that would be difficult to estimate for individual genes. In addition to theoretical development we validated TADA using realistic simulations mimicking rare, large-effect mutations affecting risk for ASD and show it has dramatically better power than other common methods of analysis. Thus TADA's integration of various kinds of WES data can be a highly effective means of identifying novel risk genes. Indeed, application of TADA to WES data from subjects with ASD and their families, as well as from a study of ASD subjects and controls, revealed several novel and promising ASD candidate genes with strong statistical support.

A homozygous mutation in KCTD7 links neuronal ceroid lipofuscinosis to the ubiquitin-proteasome system.

Neuronal ceroid lipofuscinosis (NCL) is a genetically heterogeneous group of lysosomal diseases that collectively compose the most common Mendelian form of childhood-onset neurodegeneration. It is estimated that ∼8% of individuals diagnosed with NCL by conservative clinical and histopathologic criteria have been ruled out for mutations in the nine known NCL-associated genes, suggesting that additional genes remain unidentified. To further understand the genetic underpinnings of the NCLs, we performed whole-exome sequencing on DNA samples from a Mexican family affected by a molecularly undefined form of NCL characterized by infantile-onset progressive myoclonic epilepsy (PME), vision loss, cognitive and motor regression, premature death, and prominent NCL-type storage material. Using a recessive model to filter the identified variants, we found a single homozygous variant, c.550C>T in KCTD7, that causes a p.Arg184Cys missense change in potassium channel tetramerization domain-containing protein 7 (KCTD7) in the affected individuals. The mutation was predicted to be deleterious and was absent in over 6,000 controls. The identified variant altered the localization pattern of KCTD7 and abrogated interaction with cullin-3, a ubiquitin-ligase component and known KCTD7 interactor. Intriguingly, murine cerebellar cells derived from a juvenile NCL model (CLN3) showed enrichment of endogenous KCTD7. Whereas KCTD7 mutations have previously been linked to PME without lysosomal storage, this study clearly demonstrates that KCTD7 mutations also cause a rare, infantile-onset NCL subtype designated as CLN14.

Whole-exome sequencing and homozygosity analysis implicate depolarization-regulated neuronal genes in autism.

Although autism has a clear genetic component, the high genetic heterogeneity of the disorder has been a challenge for the identification of causative genes. We used homozygosity analysis to identify probands from nonconsanguineous families that showed evidence of distant shared ancestry, suggesting potentially recessive mutations. Whole-exome sequencing of 16 probands revealed validated homozygous, potentially pathogenic recessive mutations that segregated perfectly with disease in 4/16 families. The candidate genes (UBE3B, CLTCL1, NCKAP5L, ZNF18) encode proteins involved in proteolysis, GTPase-mediated signaling, cytoskeletal organization, and other pathways. Furthermore, neuronal depolarization regulated the transcription of these genes, suggesting potential activity-dependent roles in neurons. We present a multidimensional strategy for filtering whole-exome sequence data to find candidate recessive mutations in autism, which may have broader applicability to other complex, heterogeneous disorders.

Expression of ALS-PFN1 impairs vesicular degradation in iPSC-derived microglia.

Microglia play a pivotal role in neurodegenerative disease pathogenesis, but the mechanisms underlying microglia dysfunction and toxicity remain to be elucidated. To investigate the effect of neurodegenerative disease-linked genes on the intrinsic properties of microglia, we studied microglia-like cells derived from human induced pluripotent stem cells (iPSCs), termed iMGs, harboring mutations in profilin-1 (PFN1) that are causative for amyotrophic lateral sclerosis (ALS). ALS-PFN1 iMGs exhibited evidence of lipid dysmetabolism, autophagy dysregulation and deficient phagocytosis, a canonical microglia function. Mutant PFN1 also displayed enhanced binding affinity for PI3P, a critical signaling molecule involved in autophagic and endocytic processing. Our cumulative data implicate a gain-of-toxic function for mutant PFN1 within the autophagic and endo-lysosomal pathways, as administration of rapamycin rescued phagocytic dysfunction in ALS-PFN1 iMGs. These outcomes demonstrate the utility of iMGs for neurodegenerative disease research and implicate microglial vesicular degradation pathways in the pathogenesis of these disorders.

FACS-Based Sequencing Approach to Evaluate Cell Type to Genotype Associations Using Cerebral Organoids.

Recent technological developments have led to widespread applications of large-scale transcriptomics-based sequencing methods to identify genotype-to-cell type associations. Here we describe a fluorescence-activated cell sorting (FACS)-based sequencing method to utilize CRISPR/Cas9 edited mosaic cerebral organoids to identify or validate genotype-to-cell type associations. Our approach is high-throughput and quantitative and uses internal controls to enable comparisons of the results across different antibody markers and experiments.

Lack of ethnic diversity in single-cell transcriptomics hinders cell type detection and precision medicine inclusivity.

Perhaps one of the most revolutionary next generation sequencing technologies is single-cell (SC) transcriptomics, which was recognized by Nature in 2013 as the method of the year. SC-technologies delve deep into genomics at the single-cell level, revealing previously restricted, valuable information on the identity of single cells, particularly highlighting their heterogeneity. Understanding the cellular heterogeneity of complex tissue provides insight about the gene expression and regulation across different biological and environmental conditions. This vast heterogeneity of cells and their markers makes identifying populations and sub-clusters especially difficult, even more so in rare cell types limited by the absence of rare sub-population markers. One particularly overlooked challenge is the lack of adequate ethnic representation in single-cell data. As the availability of cell types and their markers grow exponentially through new discoveries, the need to study ethnically driven heterogeneity becomes more feasible, while offering the opportunity to further elaborate ethnicity-related heterogeneity. In this commentary, we will discuss this major single-cell limitation particularly focusing on the repercussions it has on disease research, therapeutic design, and precision medicine.

Reliable multiplex generation of pooled induced pluripotent stem cells.

Reprogramming somatic cells into pluripotent stem cells (iPSCs) enables the study of systems in vitro. To increase the throughput of reprogramming, we present induction of pluripotency from pooled cells (iPPC)-an efficient, scalable, and reliable reprogramming procedure. Using our deconvolution algorithm that employs pooled sequencing of single-nucleotide polymorphisms (SNPs), we accurately estimated individual donor proportions of the pooled iPSCs. With iPPC, we concurrently reprogrammed over one hundred donor lymphoblastoid cell lines (LCLs) into iPSCs and found strong correlations of individual donors' reprogramming ability across multiple experiments. Individual donors' reprogramming ability remains consistent across both same-day replicates and multiple experimental runs, and the expression of certain immunoglobulin precursor genes may impact reprogramming ability. The pooled iPSCs were also able to differentiate into cerebral organoids. Our procedure enables a multiplex framework of using pooled libraries of donor iPSCs for downstream research and investigation of in vitro phenotypes.

Expression of ALS-PFN1 impairs vesicular degradation in iPSC-derived microglia.

Microglia play a pivotal role in neurodegenerative disease pathogenesis, but the mechanisms underlying microglia dysfunction and toxicity remain to be fully elucidated. To investigate the effect of neurodegenerative disease-linked genes on the intrinsic properties of microglia, we studied microglia-like cells derived from human induced pluripotent stem cells (iPSCs), termed iMGs, harboring mutations in profilin-1 (PFN1) that are causative for amyotrophic lateral sclerosis (ALS). ALS-PFN1 iMGs exhibited lipid dysmetabolism and deficits in phagocytosis, a critical microglia function. Our cumulative data implicate an effect of ALS-linked PFN1 on the autophagy pathway, including enhanced binding of mutant PFN1 to the autophagy signaling molecule PI3P, as an underlying cause of defective phagocytosis in ALS-PFN1 iMGs. Indeed, phagocytic processing was restored in ALS-PFN1 iMGs with Rapamycin, an inducer of autophagic flux. These outcomes demonstrate the utility of iMGs for neurodegenerative disease research and highlight microglia vesicular degradation pathways as potential therapeutic targets for these disorders.

Modeling of mitochondrial genetic polymorphisms reveals induction of heteroplasmy by pleiotropic disease locus 10398A>G.

Mitochondrial (MT) dysfunction has been associated with several neurodegenerative diseases including Alzheimer's disease (AD). While MT-copy number differences have been implicated in AD, the effect of MT heteroplasmy on AD has not been well characterized. Here, we analyzed over 1800 whole genome sequencing data from four AD cohorts in seven different tissue types to determine the extent of MT heteroplasmy present. While MT heteroplasmy was present throughout the entire MT genome for blood samples, we detected MT heteroplasmy only within the MT control region for brain samples. We observed that an MT variant 10398A>G (rs2853826) was significantly associated with overall MT heteroplasmy in brain tissue while also being linked with the largest number of distinct disease phenotypes of all annotated MT variants in MitoMap. Using gene-expression data from our brain samples, our modeling discovered several gene networks involved in mitochondrial respiratory chain and Complex I function associated with 10398A>G. The variant was also found to be an expression quantitative trait loci (eQTL) for the gene MT-ND3. We further characterized the effect of 10398A>G by phenotyping a population of lymphoblastoid cell-lines (LCLs) with and without the variant allele. Examination of RNA sequence data from these LCLs reveal that 10398A>G was an eQTL for MT-ND4. We also observed in LCLs that 10398A>G was significantly associated with overall MT heteroplasmy within the MT control region, confirming the initial findings observed in post-mortem brain tissue. These results provide novel evidence linking MT SNPs with MT heteroplasmy and open novel avenues for the investigation of pathomechanisms that are driven by this pleiotropic disease associated loci.