Trans-ethnic and Ancestry-Specific Blood-Cell Genetics in 746,667 Individuals from 5 Global Populations.

Most loci identified by GWASs have been found in populations of European ancestry (EUR). In trans-ethnic meta-analyses for 15 hematological traits in 746,667 participants, including 184,535 non-EUR individuals, we identified 5,552 trait-variant associations at p < 5 × 10-9, including 71 novel associations not found in EUR populations. We also identified 28 additional novel variants in ancestry-specific, non-EUR meta-analyses, including an IL7 missense variant in South Asians associated with lymphocyte count in vivo and IL-7 secretion levels in vitro. Fine-mapping prioritized variants annotated as functional and generated 95% credible sets that were 30% smaller when using the trans-ethnic as opposed to the EUR-only results. We explored the clinical significance and predictive value of trans-ethnic variants in multiple populations and compared genetic architecture and the effect of natural selection on these blood phenotypes between populations. Altogether, our results for hematological traits highlight the value of a more global representation of populations in genetic studies.

The Polygenic and Monogenic Basis of Blood Traits and Diseases.

Blood cells play essential roles in human health, underpinning physiological processes such as immunity, oxygen transport, and clotting, which when perturbed cause a significant global health burden. Here we integrate data from UK Biobank and a large-scale international collaborative effort, including data for 563,085 European ancestry participants, and discover 5,106 new genetic variants independently associated with 29 blood cell phenotypes covering a range of variation impacting hematopoiesis. We holistically characterize the genetic architecture of hematopoiesis, assess the relevance of the omnigenic model to blood cell phenotypes, delineate relevant hematopoietic cell states influenced by regulatory genetic variants and gene networks, identify novel splice-altering variants mediating the associations, and assess the polygenic prediction potential for blood traits and clinical disorders at the interface of complex and Mendelian genetics. These results show the power of large-scale blood cell trait GWAS to interrogate clinically meaningful variants across a wide allelic spectrum of human variation.

High-throughput splicing assays identify missense and silent splice-disruptive POU1F1 variants underlying pituitary hormone deficiency.

Pituitary hormone deficiency occurs in ∼1:4,000 live births. Approximately 3% of the cases are due to mutations in the alpha isoform of POU1F1, a pituitary-specific transcriptional activator. We found four separate heterozygous missense variants in unrelated individuals with hypopituitarism that were predicted to affect a minor isoform, POU1F1 beta, which can act as a transcriptional repressor. These variants retain repressor activity, but they shift splicing to favor the expression of the beta isoform, resulting in dominant-negative loss of function. Using a high-throughput splicing reporter assay, we tested 1,070 single-nucleotide variants in POU1F1. We identified 96 splice-disruptive variants, including 14 synonymous variants. In separate cohorts, we found two additional synonymous variants nominated by this screen that co-segregate with hypopituitarism. This study underlines the importance of evaluating the impact of variants on splicing and provides a catalog for interpretation of variants of unknown significance in POU1F1.

Burden Testing of Rare Variants Identified through Exome Sequencing via Publicly Available Control Data.

The genetic causes of many Mendelian disorders remain undefined. Factors such as lack of large multiplex families, locus heterogeneity, and incomplete penetrance hamper these efforts for many disorders. Previous work suggests that gene-based burden testing-where the aggregate burden of rare, protein-altering variants in each gene is compared between case and control subjects-might overcome some of these limitations. The increasing availability of large-scale public sequencing databases such as Genome Aggregation Database (gnomAD) can enable burden testing using these databases as controls, obviating the need for additional control sequencing for each study. However, there exist various challenges with using public databases as controls, including lack of individual-level data, differences in ancestry, and differences in sequencing platforms and data processing. To illustrate the approach of using public data as controls, we analyzed whole-exome sequencing data from 393 individuals with idiopathic hypogonadotropic hypogonadism (IHH), a rare disorder with significant locus heterogeneity and incomplete penetrance against control subjects from gnomAD (n = 123,136). We leveraged presumably benign synonymous variants to calibrate our approach. Through iterative analyses, we systematically addressed and overcame various sources of artifact that can arise when using public control data. In particular, we introduce an approach for highly adaptable variant quality filtering that leads to well-calibrated results. Our approach "re-discovered" genes previously implicated in IHH (FGFR1, TACR3, GNRHR). Furthermore, we identified a significant burden in TYRO3, a gene implicated in hypogonadotropic hypogonadism in mice. Finally, we developed a user-friendly software package TRAPD (Test Rare vAriants with Public Data) for performing gene-based burden testing against public databases.

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.

Temporal trends and yield of clinical diagnostic genetic testing in adult neurology.

While genetics evaluation is increasingly utilized in adult neurology patients, its usage and efficacy are not well characterized. Here, we report our experience with 1461 consecutive patients evaluated in an adult neurogenetics clinic at a large academic medical center between January 2015 and March 2020. Of the 1461 patients evaluated, 1215 patients were referred for the purposes of identifying a genetic diagnosis for an undiagnosed condition, 90.5% of whom underwent genetic testing. The modalities of genetic testing utilized varied across referral diagnostic categories, including a range of utilization of whole exome sequencing (WES) as an initial test in 13.9% of neuromuscular patients to 52.9% in white matter disorder patients. The usage of WES increased over time, from 7.7% of initial testing in 2015 to a peak of 27.3% in 2019. Overall, genetic testing yielded a causal genetic diagnosis in 30.7% of patients. This yield was higher in certain referring diagnosis categories, such as neuromuscular (39.0%) and epilepsy (29.8%). Our study demonstrates that evaluation at an adult neurogenetics referral center can yield diagnoses in a substantial fraction of patients. Additional research will be needed to determine optimal genetic testing strategies and cost effectiveness of adult neurogenetics evaluation.

TMEM106B Effect on cognition in Parkinson disease and frontotemporal dementia.

OBJECTIVE: Common variants near TMEM106B associate with risk of developing frontotemporal dementia (FTD). Emerging evidence suggests a role for TMEM106B in neurodegenerative processes beyond FTD. We evaluate the effect of TMEM106B genotype on cognitive decline across multiple neurogenerative diseases. METHODS: We longitudinally followed 870 subjects with diagnoses of Parkinson disease (PD; n = 179), FTD (n = 179), Alzheimer disease (AD; n = 300), memory-predominant mild cognitive impairment (MCI; n = 75), or neurologically normal control subjects (NC; n = 137) at the University of Pennsylvania (UPenn). All participants had annual Mini-Mental State Examination (MMSE; median follow-up duration = 3.0 years) and were genotyped at TMEM106B index single nucleotide polymorphism rs1990622. Genotype effects on cognition were confirmed by extending analyses to additional cognitive instruments (Mattis Dementia Rating Scale-2 [DRS-2] and Montreal Cognitive Assessment [MoCA]) and to an international validation cohort (Parkinson's Progression Markers Initiative [PPMI], N = 371). RESULTS: The TMEM106B rs1990622T allele, linked to increased risk of FTD, associated with greater MMSE decline over time in PD subjects but not in AD or MCI subjects. For FTD subjects, rs1990622T associated with more rapid decrease in MMSE only under the minor-allele, rs1990622C , dominant model. Among PD patients, rs1990622T carriers from the UPenn cohort demonstrated more rapid longitudinal decline in DRS-2 scores. Finally, in the PPMI cohort, TMEM106B risk allele carriers demonstrated more rapid longitudinal decline in MoCA scores. INTERPRETATION: Irrespective of cognitive instrument or cohort assessed, TMEM106B acts as a genetic modifier for cognitive trajectory in PD. Our results implicate lysosomal dysfunction in the pathogenesis of cognitive decline in 2 different proteinopathies. ANN NEUROL 2019;85:801-811.

Genetics of rheumatoid arthritis contributes to biology and drug discovery.

A major challenge in human genetics is to devise a systematic strategy to integrate disease-associated variants with diverse genomic and biological data sets to provide insight into disease pathogenesis and guide drug discovery for complex traits such as rheumatoid arthritis (RA). Here we performed a genome-wide association study meta-analysis in a total of >100,000 subjects of European and Asian ancestries (29,880 RA cases and 73,758 controls), by evaluating ∼10 million single-nucleotide polymorphisms. We discovered 42 novel RA risk loci at a genome-wide level of significance, bringing the total to 101 (refs 2 - 4). We devised an in silico pipeline using established bioinformatics methods based on functional annotation, cis-acting expression quantitative trait loci and pathway analyses--as well as novel methods based on genetic overlap with human primary immunodeficiency, haematological cancer somatic mutations and knockout mouse phenotypes--to identify 98 biological candidate genes at these 101 risk loci. We demonstrate that these genes are the targets of approved therapies for RA, and further suggest that drugs approved for other indications may be repurposed for the treatment of RA. Together, this comprehensive genetic study sheds light on fundamental genes, pathways and cell types that contribute to RA pathogenesis, and provides empirical evidence that the genetics of RA can provide important information for drug discovery.

Comprehensive population-based genome sequencing provides insight into hematopoietic regulatory mechanisms.

Genetic variants affecting hematopoiesis can influence commonly measured blood cell traits. To identify factors that affect hematopoiesis, we performed association studies for blood cell traits in the population-based Estonian Biobank using high-coverage whole-genome sequencing (WGS) in 2,284 samples and SNP genotyping in an additional 14,904 samples. Using up to 7,134 samples with available phenotype data, our analyses identified 17 associations across 14 blood cell traits. Integration of WGS-based fine-mapping and complementary epigenomic datasets provided evidence for causal mechanisms at several loci, including at a previously undiscovered basophil count-associated locus near the master hematopoietic transcription factor CEBPA The fine-mapped variant at this basophil count association near CEBPA overlapped an enhancer active in common myeloid progenitors and influenced its activity. In situ perturbation of this enhancer by CRISPR/Cas9 mutagenesis in hematopoietic stem and progenitor cells demonstrated that it is necessary for and specifically regulates CEBPA expression during basophil differentiation. We additionally identified basophil count-associated variation at another more pleiotropic myeloid enhancer near GATA2, highlighting regulatory mechanisms for ordered expression of master hematopoietic regulators during lineage specification. Our study illustrates how population-based genetic studies can provide key insights into poorly understood cell differentiation processes of considerable physiologic relevance.

Determinants of Power in Gene-Based Burden Testing for Monogenic Disorders.

Whole-exome sequencing has enabled new approaches for discovering genes associated with monogenic disorders. One such approach is gene-based burden testing, in which the aggregate frequency of "qualifying variants" is compared between case and control subjects for each gene. Despite substantial successes of this approach, the genetic causes for many monogenic disorders remain unknown or only partially known. It is possible that particular genetic architectures lower rates of discovery, but the influence of these factors on power has not been rigorously evaluated. Here, we leverage large-scale exome-sequencing data to create an empirically based simulation framework to evaluate the impact of key parameters (background variation rates, locus heterogeneity, mode of inheritance, penetrance) on power in gene-based burden tests in the context of monogenic disorders. Our results demonstrate that across genes, there is a wide range in sample sizes needed to achieve power due to differences in the background rate of rare variants in each gene. Increasing locus heterogeneity results in rapid increases in sample sizes needed to achieve adequate power, particularly when individual genes contribute to less than 5% of cases under a dominant model. Interestingly, incomplete penetrance as low as 10% had little effect on power due to the low prevalence of monogenic disorders. Our results suggest that moderate incomplete penetrance is not an obstacle in this gene-based burden testing approach but that dominant disorders with high locus heterogeneity will require large sample sizes. Our simulations also provide guidance on sample size needs and inform study design under various genetic architectures.

Estimating yields of prenatal carrier screening and implications for design of expanded carrier screening panels.

PURPOSE: Prenatal genetic carrier screening can identify parents at risk of having a child affected by a recessive condition. However, the conditions/genes most appropriate for screening remain a matter of debate. Estimates of carrier rates across genes are needed to guide construction of carrier screening panels. METHOD: We leveraged an exome sequencing database (n = 123,136) to estimate carrier rates across six major ancestries for 415 genes associated with severe recessive conditions. RESULTS: We found that 32.6% (East Asian) to 62.9% (Ashkenazi Jewish) of individuals are variant carriers in at least one of the 415 genes. For couples, screening all 415 genes would identify 0.17-2.52% of couples as being at risk for having a child affected by one of these conditions. Screening just the 40 genes with carrier rate >1.0% would identify more than 76% of these at-risk couples. An ancestry-specific panel designed to capture genes with carrier rates >1.0% would include 5 to 28 genes, while a comparable panethnic panel would include 40 genes. CONCLUSION: Our work guides the design of carrier screening panels and provides data to assist in counseling prospective parents. Our results highlight a high cumulative carrier rate across genes, underscoring the need for careful selection of genes for screening.

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.

The Importance of Offering Exome or Genome Sequencing in Adult Neuromuscular Clinics.

Advances in gene-specific therapeutics for patients with neuromuscular disorders (NMDs) have brought increased attention to the importance of genetic diagnosis. Genetic testing practices vary among adult neuromuscular clinics, with multi-gene panel testing currently being the most common approach; follow-up testing using broad-based methods, such as exome or genome sequencing, is less consistently offered. Here, we use five case examples to illustrate the unique ability of broad-based testing to improve diagnostic yield, resulting in identification of SORD-neuropathy, HADHB-related disease, ATXN2-ALS, MECP2 related progressive gait decline and spasticity, and DNMT1-related cerebellar ataxia, deafness, narcolepsy, and hereditary sensory neuropathy type 1E. We describe in each case the technological advantages that enabled identification of the causal gene, and the resultant clinical and personal implications for the patient, demonstrating the importance of offering exome or genome sequencing to adults with NMDs.

Inferring compound heterozygosity from large-scale exome sequencing data.

Recessive diseases arise when both copies of a gene are impacted by a damaging genetic variant. When a patient carries two potentially causal variants in a gene, accurate diagnosis requires determining that these variants occur on different copies of the chromosome (that is, are in trans) rather than on the same copy (that is, in cis). However, current approaches for determining phase, beyond parental testing, are limited in clinical settings. Here we developed a strategy for inferring phase for rare variant pairs within genes, leveraging genotypes observed in the Genome Aggregation Database (v2, n = 125,748 exomes). Our approach estimates phase with 96% accuracy, both in trio data and in patients with Mendelian conditions and presumed causal compound heterozygous variants. We provide a public resource of phasing estimates for coding variants and counts per gene of rare variants in trans that can aid interpretation of rare co-occurring variants in the context of recessive disease.

Genetic testing in adults with neurologic disorders: indications, approach, and clinical impacts.

The role of genetic testing in neurologic clinical practice has increased dramatically in recent years, driven by research on genetic causes of neurologic disease and increased availability of genetic sequencing technology. Genetic testing is now indicated for adults with a wide range of common neurologic conditions. The potential clinical impacts of a genetic diagnosis are also rapidly expanding, with a growing list of gene-specific treatments and clinical trials, in addition to important implications for prognosis, surveillance, family planning, and diagnostic closure. The goals of this review are to provide practical guidance for clinicians about the role of genetics in their practice and to provide the neuroscience research community with a broad survey of current progress in this field. We aim to answer three questions for the neurologist in practice: Which of my patients need genetic testing? What testing should I order? And how will genetic testing help my patient? We focus on common neurologic disorders and presentations to the neurology clinic. For each condition, we review the most current guidelines and evidence regarding indications for genetic testing, expected diagnostic yield, and recommended testing approach. We also focus on clinical impacts of genetic diagnoses, highlighting a number of gene-specific therapies recently approved for clinical use, and a rapidly expanding landscape of gene-specific clinical trials, many using novel nucleotide-based therapeutic modalities like antisense oligonucleotides and gene transfer. We anticipate that more widespread use of genetic testing will help advance therapeutic development and improve the care, and outcomes, of patients with neurologic conditions.

Redefining the progeroid form of Ehlers-Danlos syndrome: report of the fourth patient with B4GALT7 deficiency and review of the literature.

Proteoglycans are a component of the extracellular matrix and are critical for cellular and tissue function. Mutations in proteoglycan components and enzymes involved in proteoglycan synthesis have been implicated in several growth disorders, with common features including short stature and skeletal dysplasia. For example, mutations in B4GALT7, a gene whose protein product catalyzes proteoglycan synthesis, have been associated with the rare progeroid variant of Ehlers-Danlos syndrome. Here, we conducted exome sequencing in a patient with a previously undiagnosed growth disorder and identified compound heterozygous mutations in B4GALT7. This patient is just the fourth individual with genetically confirmed progeroid variant of Ehlers-Danlos syndrome. The mutations include a previously characterized c.808C>T p.Arg270Cys substitution, and a novel c.122T>C p.Leu41Pro substitution. We demonstrate that the novel mutation caused decreased levels of the enzyme, supporting the pathogenicity of the mutation. Our report identifies a novel mutation in B4GALT7 causing the progeroid variant of Ehlers-Danlos syndrome and contributes an extensive phenotypic characterization of a patient with the syndrome. We also reviewed the previous literature in addition to the present patient, and conclude that the key features associated with B4GALT7 deficiency are short stature, developmental anomalies of the forearm bones and elbow, and bowing of the extremities, in addition to the classic features of Ehlers-Danlos syndrome. This report helps define the phenotype of the progeroid variant of Ehlers-Danlos syndrome and furthers our understanding of the effect of proteoglycan defects in growth disorders.

Polygenic burden of short tandem repeat expansions promote risk for Alzheimer's disease.

Studies of the genetics of Alzheimer's disease (AD) have largely focused on single nucleotide variants and short insertions/deletions. However, most of the disease heritability has yet to be uncovered, suggesting that there is substantial genetic risk conferred by other forms of genetic variation. There are over one million short tandem repeats (STRs) in the genome, and their link to AD risk has not been assessed. As pathogenic expansions of STR cause over 30 neurologic diseases, it is important to ascertain whether STRs may also be implicated in AD risk. Here, we genotyped 321,742 polymorphic STR tracts genome-wide using PCR-free whole genome sequencing data from 2,981 individuals (1,489 AD case and 1,492 control individuals). We implemented an approach to identify STR expansions as STRs with tract lengths that are outliers from the population. We then tested for differences in aggregate burden of expansions in case versus control individuals. AD patients had a 1.19-fold increase of STR expansions compared to healthy elderly controls (p=8.27×10-3, two-sided Mann Whitney test). Individuals carrying > 30 STR expansions had 3.62-fold higher odds of having AD and had more severe AD neuropathology. AD STR expansions were highly enriched within active promoters in post-mortem hippocampal brain tissues and particularly within SINE-VNTR-Alu (SVA) retrotransposons. Together, these results demonstrate that expanded STRs within active promoter regions of the genome promote risk of AD.

Inferring compound heterozygosity from large-scale exome sequencing data.

Recessive diseases arise when both the maternal and the paternal copies of a gene are impacted by a damaging genetic variant in the affected individual. When a patient carries two different potentially causal variants in a gene for a given disorder, accurate diagnosis requires determining that these two variants occur on different copies of the chromosome (i.e., are in trans) rather than on the same copy (i.e. in cis). However, current approaches for determining phase, beyond parental testing, are limited in clinical settings. We developed a strategy for inferring phase for rare variant pairs within genes, leveraging genotypes observed in exome sequencing data from the Genome Aggregation Database (gnomAD v2, n=125,748). When applied to trio data where phase can be determined by transmission, our approach estimates phase with 95.7% accuracy and remains accurate even for very rare variants (allele frequency < 1×10-4). We also correctly phase 95.9% of variant pairs in a set of 293 patients with Mendelian conditions carrying presumed causal compound heterozygous variants. We provide a public resource of phasing estimates from gnomAD, including phasing estimates for coding variants across the genome and counts per gene of rare variants in trans, that can aid interpretation of rare co-occurring variants in the context of recessive disease.