Search

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions in lipid metabolism. Large-scale whole-genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess more associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with measurement of blood lipids and lipoproteins (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare-variant aggregate association tests using the STAAR (variant-set test for association using annotation information) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare-coding variants in nearby protein-coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500-kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variation and rare protein-coding variation at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNAs.

Despite the growing number of genome-wide association studies (GWASs), it remains unclear to what extent gene-by-gene and gene-by-environment interactions influence complex traits in humans. The magnitude of genetic interactions in complex traits has been difficult to quantify because GWASs are generally underpowered to detect individual interactions of small effect. Here, we develop a method to test for genetic interactions that aggregates information across all trait-associated loci. Specifically, we test whether SNPs in regions of European ancestry shared between European American and admixed African American individuals have the same causal effect sizes. We hypothesize that in African Americans, the presence of genetic interactions will drive the causal effect sizes of SNPs in regions of European ancestry to be more similar to those of SNPs in regions of African ancestry. We apply our method to two traits: gene expression in 296 African Americans and 482 European Americans in the Multi-Ethnic Study of Atherosclerosis (MESA) and low-density lipoprotein cholesterol (LDL-C) in 74K African Americans and 296K European Americans in the Million Veteran Program (MVP). We find significant evidence for genetic interactions in our analysis of gene expression; for LDL-C, we observe a similar point estimate, although this is not significant, most likely due to lower statistical power. These results suggest that gene-by-gene or gene-by-environment interactions modify the effect sizes of causal variants in human complex traits.

Whole-genome sequencing (WGS), a powerful tool for detecting novel coding and non-coding disease-causing variants, has largely been applied to clinical diagnosis of inherited disorders. Here we leveraged WGS data in up to 62,653 ethnically diverse participants from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program and assessed statistical association of variants with seven red blood cell (RBC) quantitative traits. We discovered 14 single variant-RBC trait associations at 12 genomic loci, which have not been reported previously. Several of the RBC trait-variant associations (RPN1, ELL2, MIDN, HBB, HBA1, PIEZO1, and G6PD) were replicated in independent GWAS datasets imputed to the TOPMed reference panel. Most of these discovered variants are rare/low frequency, and several are observed disproportionately among non-European Ancestry (African, Hispanic/Latino, or East Asian) populations. We identified a 3 bp indel p.Lys2169del (g.88717175_88717177TCT[4]) (common only in the Ashkenazi Jewish population) of PIEZO1, a gene responsible for the Mendelian red cell disorder hereditary xerocytosis (MIM: 194380), associated with higher mean corpuscular hemoglobin concentration (MCHC). In stepwise conditional analysis and in gene-based rare variant aggregated association analysis, we identified several of the variants in HBB, HBA1, TMPRSS6, and G6PD that represent the carrier state for known coding, promoter, or splice site loss-of-function variants that cause inherited RBC disorders. Finally, we applied base and nuclease editing to demonstrate that the sentinel variant rs112097551 (nearest gene RPN1) acts through a cis-regulatory element that exerts long-range control of the gene RUVBL1 which is essential for hematopoiesis. Together, these results demonstrate the utility of WGS in ethnically diverse population-based samples and gene editing for expanding knowledge of the genetic architecture of quantitative hematologic traits and suggest a continuum between complex trait and Mendelian red cell disorders.

Cell division terminates with cytokinesis and cellular separation. Autosomal-recessive primary microcephaly (MCPH) is a neurodevelopmental disorder characterized by a reduction in brain and head size at birth in addition to non-progressive intellectual disability. MCPH is genetically heterogeneous, and 16 loci are known to be associated with loss-of-function mutations predominantly affecting centrosomal-associated proteins, but the multiple roles of centrosomes in cellular function has left questions about etiology. Here, we identified three families affected by homozygous missense mutations in CIT, encoding citron rho-interacting kinase (CIT), which has established roles in cytokinesis. All mutations caused substitution of conserved amino acid residues in the kinase domain and impaired kinase activity. Neural progenitors that were differentiated from induced pluripotent stem cells (iPSCs) derived from individuals with these mutations exhibited abnormal cytokinesis with delayed mitosis, multipolar spindles, and increased apoptosis, rescued by CRISPR/Cas9 genome editing. Our results highlight the importance of cytokinesis in the pathology of primary microcephaly.

Elevated low-density lipoprotein cholesterol (LDL-C) is a treatable, heritable risk factor for cardiovascular disease. Genome-wide association studies (GWASs) have identified 157 variants associated with lipid levels but are not well suited to assess the impact of rare and low-frequency variants. To determine whether rare or low-frequency coding variants are associated with LDL-C, we exome sequenced 2,005 individuals, including 554 individuals selected for extreme LDL-C (>98(th) or

The corpus callosum is the principal cerebral commissure connecting the right and left hemispheres. The development of the corpus callosum is under tight genetic control, as demonstrated by abnormalities in its development in more than 1,000 genetic syndromes. We recruited more than 25 families in which members affected with corpus callosum hypoplasia (CCH) lacked syndromic features and had consanguineous parents, suggesting recessive causes. Exome sequence analysis identified C12orf57 mutations at the initiator methionine codon in four different families. C12orf57 is ubiquitously expressed and encodes a poorly annotated 126 amino acid protein of unknown function. This protein is without significant paralogs but has been tightly conserved across evolution. Our data suggest that this conserved gene is required for development of the human corpus callosum.

A genome-wide association study is undertaken in a high-risk isolate for multiple sclerosis (MS) in order to identify previously unreported high-impact variants associated with MS. Results reveal that SNP (rs744166) within the STAT3 gene is associated with MS.

Bulk-tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, and context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from the blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell-type proportions, we demonstrate that cell-type iQTLs could be considered as proxies for cell-type-specific QTL effects, particularly for the most abundant cell type in the tissue. The interpretation of age iQTLs, however, warrants caution because the moderation effect of age on the genotype and molecular phenotype association could be mediated by changes in cell-type composition. Finally, we show that cell-type iQTLs contribute to cell-type-specific enrichment of diseases that, in combination with additional functional data, could guide future functional studies. Overall, this study highlights the use of iQTLs to gain insights into the context specificity of regulatory effects. To identify regulatory variants with plasticity in effect size, we performed interaction molecular quantitative trait loci (iQTL) mapping with cell-type abundance, age, sex, and smoking as the environmental factors. Our results highlight the usefulness of iQTLs for gaining insights into the context specificity of regulatory effects. [ABSTRACT FROM AUTHOR]

Finnish samples have been extensively utilized in studying single-gene disorders, where the founder effect has clearly aided in discovery, and more recently in genome-wide association studies of complex traits, where the founder effect has had less obvious impacts. As the field starts to explore rare variants’ contribution to polygenic traits, it is of great importance to characterize and confirm the Finnish founder effect in sequencing data and to assess its implications for rare-variant association studies. Here, we employ forward simulation, guided by empirical deep resequencing data, to model the genetic architecture of quantitative polygenic traits in both the general European and the Finnish populations simultaneously. We demonstrate that power of rare-variant association tests is higher in the Finnish population, especially when variants’ phenotypic effects are tightly coupled with fitness effects and therefore reflect a greater contribution of rarer variants. SKAT-O, variable-threshold tests, and single-variant tests are more powerful than other rare-variant methods in the Finnish population across a range of genetic models. We also compare the relative power and efficiency of exome array genotyping to those of high-coverage exome sequencing. At a fixed cost, less expensive genotyping strategies have far greater power than sequencing; in a fixed number of samples, however, genotyping arrays miss a substantial portion of genetic signals detected in sequencing, even in the Finnish founder population. As genetic studies probe sequence variation at greater depth in more diverse populations, our simulation approach provides a framework for evaluating various study designs for gene discovery. [Copyright &y& Elsevier]

Joubert syndrome and related disorders (JSRDs) are genetically heterogeneous and characterized by a distinctive mid-hindbrain malformation. Causative mutations lead to primary cilia dysfunction, which often results in variable involvement of other organs such as the liver, retina, and kidney. We identified predicted null mutations in CSPP1 in six individuals affected by classical JSRDs. CSPP1 encodes a protein localized to centrosomes and spindle poles, as well as to the primary cilium. Despite the known interaction between CSPP1 and nephronophthisis-associated proteins, none of the affected individuals in our cohort presented with kidney disease, and further, screening of a large cohort of individuals with nephronophthisis demonstrated no mutations. CSPP1 is broadly expressed in neural tissue, and its encoded protein localizes to the primary cilium in an in vitro model of human neurogenesis. Here, we show abrogated protein levels and ciliogenesis in affected fibroblasts. Our data thus suggest that CSPP1 is involved in neural-specific functions of primary cilia. [ABSTRACT FROM AUTHOR]

Cobblestone brain malformation (COB) is a neuronal migration disorder characterized by protrusions of neurons beyond the first cortical layer at the pial surface of the brain. It is usually seen in association with dystroglycanopathy types of congenital muscular dystrophies (CMDs) and ocular abnormalities termed muscle-eye-brain disease. Here we report homozygous deleterious mutations in LAMB1, encoding laminin subunit beta-1, in two families with autosomal-recessive COB. Affected individuals displayed a constellation of brain malformations including cortical gyral and white-matter signal abnormalities, severe cerebellar dysplasia, brainstem hypoplasia, and occipital encephalocele, but they had less apparent ocular or muscular abnormalities than are typically observed in COB. LAMB1 is localized to the pial basement membrane, suggesting that defective connection between radial glial cells and the pial surface mediated by LAMB1 leads to this malformation. [ABSTRACT FROM AUTHOR]

Rare sarcomere protein variants cause dominant hypertrophic and dilated cardiomyopathies. To evaluate whether allelic variants in eight sarcomere genes are associated with cardiac morphology and function in the community, we sequenced 3,600 individuals from the Framingham Heart Study (FHS) and Jackson Heart Study (JHS) cohorts. Out of the total, 11.2% of individuals had one or more rare nonsynonymous sarcomere variants. The prevalence of likely pathogenic sarcomere variants was 0.6%, twice the previous estimates; however, only four of the 22 individuals had clinical manifestations of hypertrophic cardiomyopathy. Rare sarcomere variants were associated with an increased risk for adverse cardiovascular events (hazard ratio: 2.3) in the FHS cohort, suggesting that cardiovascular risk assessment in the general population can benefit from rare variant analysis. [ABSTRACT FROM AUTHOR]

Delineation of structural variants (SVs) at sequence resolution in highly repetitive genomic regions has long been intractable. The sequence properties, origins, and functional effects of classes of genomic rearrangements such as ring chromosomes and Robertsonian translocations thus remain unknown. To resolve these complex structures, we leveraged several recent milestones in the field, including (1) the emergence of long-read sequencing, (2) the gapless telomere-to-telomere (T2T) assembly, and (3) a tool (BigClipper) to discover chromosomal rearrangements from long reads. We applied these technologies across 13 cases with ring chromosomes, Robertsonian translocations, and complex SVs that were unresolved by short reads, followed by validation using optical genome mapping (OGM). Our analyses resolved 10 of 13 cases, including a Robertsonian translocation and all ring chromosomes. Multiple breakpoints were localized to genomic regions previously recalcitrant to sequencing such as acrocentric p-arms, ribosomal DNA arrays, and telomeric repeats, and involved complex structures such as a deletion-inversion and interchromosomal dispersed duplications. We further performed methylation profiling from long-read data to discover phased differential methylation in a gene promoter proximal to a ring fusion, suggesting a long-range position effect (LRPE) with heterochromatin spreading. Breakpoint sequences suggested mechanisms of SV formation such as microhomology-mediated and non-homologous end-joining, as well as non-allelic homologous recombination. These methods provide some of the first glimpses into the sequence resolution of Robertsonian translocations and illuminate the structural diversity of ring chromosomes and complex chromosomal rearrangements with implications for genome biology, prediction of LRPEs from integrated multi-omics technologies, and molecular diagnostics in rare disease cases., Competing Interests: Declaration of interests The Talkowski laboratory receives reagents and/or research support from Illumina Inc and Microsoft Inc. Bionano Genomics provided data and analysis support for Optical Genome Mapping. J.H. and A.R.H. are employees of Bionano Genomics, Inc., (Copyright © 2024 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

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., (Copyright © 2018 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Dendritic spines represent the major site of neuronal activity in the brain; they serve as the receiving point for neurotransmitters and undergo rapid activity-dependent morphological changes that correlate with learning and memory. Using a combination of homozygosity mapping and next-generation sequencing in two consanguineous families affected by nonsyndromic autosomal-recessive intellectual disability, we identified truncating mutations in formin 2 (FMN2), encoding a protein that belongs to the formin family of actin cytoskeleton nucleation factors and is highly expressed in the maturing brain. We found that FMN2 localizes to punctae along dendrites and that germline inactivation of mouse Fmn2 resulted in animals with decreased spine density; such mice were previously demonstrated to have a conditioned fear-learning defect. Furthermore, patient neural cells derived from induced pluripotent stem cells showed correlated decreased synaptic density. Thus, FMN2 mutations link intellectual disability either directly or indirectly to the regulation of actin-mediated synaptic spine density., (Copyright © 2014 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Whole-exome sequencing (WES), which analyzes the coding sequence of most annotated genes in the human genome, is an ideal approach to studying fully penetrant autosomal-recessive diseases, and it has been very powerful in identifying disease-causing mutations even when enrollment of affected individuals is limited by reduced survival. In this study, we combined WES with homozygosity analysis of consanguineous pedigrees, which are informative even when a single affected individual is available, to identify genetic mutations responsible for Walker-Warburg syndrome (WWS), a genetically heterogeneous autosomal-recessive disorder that severely affects the development of the brain, eyes, and muscle. Mutations in seven genes are known to cause WWS and explain 50%-60% of cases, but multiple additional genes are expected to be mutated because unexplained cases show suggestive linkage to diverse loci. Using WES in consanguineous WWS-affected families, we found multiple deleterious mutations in GTDC2 (also known as AGO61). GTDC2's predicted role as an uncharacterized glycosyltransferase is consistent with the function of other genes that are known to be mutated in WWS and that are involved in the glycosylation of the transmembrane receptor dystroglycan. Therefore, to explore the role of GTDC2 loss of function during development, we used morpholino-mediated knockdown of its zebrafish ortholog, gtdc2. We found that gtdc2 knockdown in zebrafish replicates all WWS features (hydrocephalus, ocular defects, and muscular dystrophy), strongly suggesting that GTDC2 mutations cause WWS., (Copyright © 2012 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)