Gene Discovery for Complex Traits: Lessons from Africa.

The genetics of African populations reveals an otherwise "missing layer" of human variation that arose between 100,000 and 5 million years ago. Both the vast number of these ancient variants and the selective pressures they survived yield insights into genes responsible for complex traits in all populations.

Spatial and temporal mapping of de novo mutations in schizophrenia to a fetal prefrontal cortical network.

Genes disrupted in schizophrenia may be revealed by de novo mutations in affected persons from otherwise healthy families. Furthermore, during normal brain development, genes are expressed in patterns specific to developmental stage and neuroanatomical structure. We identified de novo mutations in persons with schizophrenia and then mapped the responsible genes onto transcriptome profiles of normal human brain tissues from age 13 weeks gestation to adulthood. In the dorsolateral and ventrolateral prefrontal cortex during fetal development, genes harboring damaging de novo mutations in schizophrenia formed a network significantly enriched for transcriptional coexpression and protein interaction. The 50 genes in the network function in neuronal migration, synaptic transmission, signaling, transcriptional regulation, and transport. These results suggest that disruptions of fetal prefrontal cortical neurogenesis are critical to the pathophysiology of schizophrenia. These results also support the feasibility of integrating genomic and transcriptome analyses to map critical neurodevelopmental processes in time and space in the brain.

Long-read DNA and cDNA sequencing identify cancer-predisposing deep intronic variation in tumor-suppressor genes.

The vast majority of deeply intronic genomic variants are benign, but some extremely rare or private deep intronic variants lead to exonification of intronic sequence with abnormal transcriptional consequences. Damaging variants of this class are likely underreported as causes of disease for several reasons: Most clinical DNA and RNA testing does not include full intronic sequences; many of these variants lie in complex repetitive regions that cannot be aligned from short-read whole-genome sequence; and, until recently, consequences of deep intronic variants were not accurately predicted by in silico tools. We evaluated the frequency and consequences of rare deep intronic variants for families severely affected with breast, ovarian, pancreatic, and/or metastatic prostate cancer, but with no causal variant identified by any previous genomic or cDNA-based approach. For 10 tumor-suppressor genes, we used multiplexed adaptive sampling long-read DNA sequencing and cDNA sequencing, based on patient-derived DNA and RNA, to systematically evaluate deep intronic variation. We identified all variants across the full genomic loci of targeted genes, applied the in silico tools SpliceAI and Pangolin to predict variants of functional consequence, and then carried out long-read cDNA sequencing to identify aberrant transcripts. For eight of the 120 (6%) previously unsolved families, rare deep intronic variants in BRCA1, PALB2, and ATM create intronic pseudoexons that are spliced into transcripts, leading to premature truncations. These results suggest that long-read DNA and cDNA sequencing can be integrated into variant discovery, with strategies for accurately characterizing pathogenic variants.

Host genotype-specific therapies can optimize the inflammatory response to mycobacterial infections.

Susceptibility to tuberculosis is historically ascribed to an inadequate immune response that fails to control infecting mycobacteria. In zebrafish, we find that susceptibility to Mycobacterium marinum can result from either inadequate or excessive acute inflammation. Modulation of the leukotriene A(4) hydrolase (LTA4H) locus, which controls the balance of pro- and anti-inflammatory eicosanoids, reveals two distinct molecular routes to mycobacterial susceptibility converging on dysregulated TNF levels: inadequate inflammation caused by excess lipoxins and hyperinflammation driven by excess leukotriene B(4). We identify therapies that specifically target each of these extremes. In humans, we identify a single nucleotide polymorphism in the LTA4H promoter that regulates its transcriptional activity. In tuberculous meningitis, the polymorphism is associated with inflammatory cell recruitment, patient survival and response to adjunctive anti-inflammatory therapy. Together, our findings suggest that host-directed therapies tailored to patient LTA4H genotypes may counter detrimental effects of either extreme of inflammation.

Genetic heterogeneity in human disease.

Strong evidence suggests that rare mutations of severe effect are responsible for a substantial portion of complex human disease. Evolutionary forces generate vast genetic heterogeneity in human illness by introducing many new variants in each generation. Current sequencing technologies offer the possibility of finding rare disease-causing mutations and the genes that harbor them.

The lta4h locus modulates susceptibility to mycobacterial infection in zebrafish and humans.

Exposure to Mycobacterium tuberculosis produces varied early outcomes, ranging from resistance to infection to progressive disease. Here we report results from a forward genetic screen in zebrafish larvae that identify multiple mutant classes with distinct patterns of innate susceptibility to Mycobacterium marinum. A hypersusceptible mutant maps to the lta4h locus encoding leukotriene A(4) hydrolase, which catalyzes the final step in the synthesis of leukotriene B(4) (LTB(4)), a potent chemoattractant and proinflammatory eicosanoid. lta4h mutations confer hypersusceptibility independent of LTB(4) reduction, by redirecting eicosanoid substrates to anti-inflammatory lipoxins. The resultant anti-inflammatory state permits increased mycobacterial proliferation by limiting production of tumor necrosis factor. In humans, we find that protection from both tuberculosis and multibacillary leprosy is associated with heterozygosity for LTA4H polymorphisms that have previously been correlated with differential LTB(4) production. Our results suggest conserved roles for balanced eicosanoid production in vertebrate resistance to mycobacterial infection.

2020 William Allan Award address: genetics as a way of thinking-cultural inheritance from our teachers.

This article is based on the address given by the author at the 2020 virtual meeting of The American Society of Human Genetics (ASHG) on October 26, 2020. The video of the original address can be found at the ASHG website.

Targeted long-read sequencing identifies missing disease-causing variation.

Despite widespread clinical genetic testing, many individuals with suspected genetic conditions lack a precise diagnosis, limiting their opportunity to take advantage of state-of-the-art treatments. In some cases, testing reveals difficult-to-evaluate structural differences, candidate variants that do not fully explain the phenotype, single pathogenic variants in recessive disorders, or no variants in genes of interest. Thus, there is a need for better tools to identify a precise genetic diagnosis in individuals when conventional testing approaches have been exhausted. We performed targeted long-read sequencing (T-LRS) using adaptive sampling on the Oxford Nanopore platform on 40 individuals, 10 of whom lacked a complete molecular diagnosis. We computationally targeted up to 151 Mbp of sequence per individual and searched for pathogenic substitutions, structural variants, and methylation differences using a single data source. We detected all genomic aberrations-including single-nucleotide variants, copy number changes, repeat expansions, and methylation differences-identified by prior clinical testing. In 8/8 individuals with complex structural rearrangements, T-LRS enabled more precise resolution of the mutation, leading to changes in clinical management in one case. In ten individuals with suspected Mendelian conditions lacking a precise genetic diagnosis, T-LRS identified pathogenic or likely pathogenic variants in six and variants of uncertain significance in two others. T-LRS accurately identifies pathogenic structural variants, resolves complex rearrangements, and identifies Mendelian variants not detected by other technologies. T-LRS represents an efficient and cost-effective strategy to evaluate high-priority genes and regions or complex clinical testing results.

Clinically Complex LRBA Deficiency Due to a Founder Allele in the Georgian Jewish Population.

Pathogenic variants in LRBA, encoding the LPS Responsive Beige-Like Anchor (LRBA) protein, are responsible for recessive, early-onset hypogammaglobulinemia, severe multi-organ autoimmunity, and lymphoproliferation, with increased risk for malignancy. LRBA deficiency has a wide clinical spectrum with variable age of onset and disease severity. Three apparently unrelated patients with LRBA deficiency, of Georgian Jewish descent, were homozygous for LRBA c.6640C > T, p.R2214*, leading to a stop upstream of the LRBA BEACH domain. Despite carrying the same LRBA genotype, the three patients differed in clinical course: the first patient was asymptomatic until age 25 years; the second presented with failure to thrive at age 3 months; and the third presented at age 7 years with immune cytopenias and severe infections. Two of the patients developed malignancies: the first patient was diagnosed with recurrent Hodgkin's disease at age 36 years, and the second patient developed aggressive gastric cancer at age 15 years. Among Georgian Jews, the carrier frequency of the LRBA p.R2214* allele was 1.6% (4 of 236 Georgian Jewish controls). The allele was absent from other populations. Haplotype analysis showed a shared origin of the mutation. These three patients revealed a pathogenic LRBA founder allele in the Georgian Jewish population, support the diverse and complex clinical spectrum of LRBA deficiency, and support the possibility that LRBA deficiency predisposes to malignancy.

Diagnostic yield of chromosomal microarray and trio whole exome sequencing in cryptogenic cerebral palsy.

OBJECTIVE: To determine the yield of genetic diagnoses using chromosomal microarray (CMA) and trio whole exome sequencing (WES), separately and combined, among patients with cryptogenic cerebral palsy (CP). METHODS: Trio WES of patients with prior CMA analysis for cryptogenic CP, defined as disabling, non-progressive motor symptoms beginning before the age of 3 years without known cause. RESULTS: Given both CMA analysis and trio WES, clinically significant genetic findings were identified for 58% of patients (26 of 45). Diagnoses were eight large CNVs detected by CMA and 18 point mutations detected by trio WES. None had more than one severe mutation. Approximately half of events (14 of 26) were de novo. Yield was significantly higher in patients with CP with comorbidities (69%, 22 of 32) than in those with pure motor CP (31%, 4 of 13; p=0.02). Among patients with genetic diagnoses, CNVs were more frequent than point mutations among patients with congenital anomalies (OR 7.8, 95% CI 1.2 to 52.4) or major dysmorphic features (OR 10.5, 95% CI 1.4 to 73.7). Clinically significant mutations were identified in 18 different genes: 14 with known involvement in CP-related disorders and 4 responsible for other neurodevelopmental conditions. Three possible new candidate genes for CP were ARGEF10, RTF1 and TAOK3. CONCLUSIONS: Cryptogenic CP is genetically highly heterogeneous. Genomic analysis has a high yield and is warranted in all these patients. Trio WES has higher yield than CMA, except in patients with congenital anomalies or major dysmorphic features, but these methods are complementary. Patients with negative results with one approach should also be tested by the other.

Genomic analysis of inherited hearing loss in the Palestinian population.

The genetic characterization of a common phenotype for an entire population reveals both the causes of that phenotype for that place and the power of family-based, population-wide genomic analysis for gene and mutation discovery. We characterized the genetics of hearing loss throughout the Palestinian population, enrolling 2,198 participants from 491 families from all parts of the West Bank and Gaza. In Palestinian families with no prior history of hearing loss, we estimate that 56% of hearing loss is genetic and 44% is not genetic. For the great majority (87%) of families with inherited hearing loss, panel-based genomic DNA sequencing, followed by segregation analysis of large kindreds and transcriptional analysis of participant RNA, enabled identification of the causal genes and mutations, including at distant noncoding sites. Genetic heterogeneity of hearing loss was striking with respect to both genes and alleles: The 337 solved families harbored 143 different mutations in 48 different genes. For one in four solved families, a transcription-altering mutation was the responsible allele. Many of these mutations were cryptic, either exonic alterations of splice enhancers or silencers or deeply intronic events. Experimentally calibrated in silico analysis of transcriptional effects yielded inferences of high confidence for effects on splicing even of mutations in genes not expressed in accessible tissue. Most (58%) of all hearing loss in the population was attributable to consanguinity. Given the ongoing decline in consanguineous marriage, inherited hearing loss will likely be much rarer in the next generation.

Molecular diagnosis of childhood immune dysregulation, polyendocrinopathy, and enteropathy, and implications for clinical management.

BACKGROUND: Most patients with childhood-onset immune dysregulation, polyendocrinopathy, and enteropathy have no genetic diagnosis for their illness. These patients may undergo empirical immunosuppressive treatment with highly variable outcomes. OBJECTIVE: We sought to determine the genetic basis of disease in patients referred with Immune dysregulation, polyendocrinopathy, enteropathy, X-linked-like (IPEX-like) disease, but with no mutation in FOXP3; then to assess consequences of genetic diagnoses for clinical management. METHODS: Genomic DNA was sequenced using a panel of 462 genes implicated in inborn errors of immunity. Candidate mutations were characterized by genomic, transcriptional, and (for some) protein analysis. RESULTS: Of 123 patients with FOXP3-negative IPEX-like disease, 48 (39%) carried damaging germline mutations in 1 of the following 27 genes: AIRE, BACH2, BCL11B, CARD11, CARD14, CTLA4, IRF2BP2, ITCH, JAK1, KMT2D, LRBA, MYO5B, NFKB1, NLRC4, POLA1, POMP, RAG1, SH2D1A, SKIV2L, STAT1, STAT3, TNFAIP3, TNFRSF6/FAS, TNRSF13B/TACI, TOM1, TTC37, and XIAP. Many of these genes had not been previously associated with an IPEX-like diagnosis. For 42 of the 48 patients with genetic diagnoses, knowing the critical gene could have altered therapeutic management, including recommendations for targeted treatments and for or against hematopoietic cell transplantation. CONCLUSIONS: Many childhood disorders now bundled as "IPEX-like" disease are caused by individually rare, severe mutations in immune regulation genes. Most genetic diagnoses of these conditions yield clinically actionable findings. Barriers are lack of testing or lack of repeat testing if older technologies failed to provide a diagnosis.

Mutational spectrum of breast cancer susceptibility genes among women ascertained in a cancer risk clinic in Northeast Brazil.

PURPOSE: There is a paucity of data on the spectrum and prevalence of pathogenic variants among women of African ancestry in the Northeast region of Brazil. METHODS: We performed BROCA panel sequencing to identify inherited loss-of-function variants in breast cancer susceptibility genes among 292 Brazilian women referred to a single institution cancer risk assessment program. RESULTS: The study included a convenient cohort of 173 women with invasive breast cancer (cases) and 119 women who were cancer-free at the time of ascertainment. The majority of the women self-reported as African-descended (67% for cases and 90.8% for unaffected volunteers). Thirty-seven pathogenic variants were found in 36 (20.8%) patients. While the spectrum of pathogenic variants was heterogeneous, the majority (70.3%) of the pathogenic variants were detected in high-risk genes BRCA1, BRCA2, PALB2, and TP53. Pathogenic variants were also found in the ATM, BARD1, BRIP1, FAM175A, FANCM, NBN, and SLX4 genes in 6.4% of the affected women. Four recurrent pathogenic variants were detected in 11 patients of African ancestry. Only one unaffected woman had a pathogenic variant in the RAD51C gene. Different risk assessment models examined performed well in predicting risk of carrying germline loss-of-function variants in BRCA1 and/or BRCA2 in breast cancer cases. CONCLUSION: The high prevalence and heterogenous spectrum of pathogenic variants identified among self-reported African descendants in Northeast Brazil is consistent with studies in other African ancestry populations with a high burden of aggressive young onset breast cancer. It underscores the need to integrate comprehensive cancer risk assessment and genomic testing in the management of newly diagnosed Black women with breast cancer across the African Diaspora, enabling improved cancer control in admixed underserved and understudied populations.

CRISPR-Cas9/long-read sequencing approach to identify cryptic mutations in BRCA1 and other tumour suppressor genes.

Current clinical approaches for mutation discovery are based on short sequence reads (100-300 bp) of exons and flanking splice sites targeted by multigene panels or whole exomes. Short-read sequencing is highly accurate for detection of single nucleotide variants, small indels and simple copy number differences but is of limited use for identifying complex insertions and deletions and other structural rearrangements. We used CRISPR-Cas9 to excise complete BRCA1 and BRCA2 genomic regions from lymphoblast cells of patients with breast cancer, then sequenced these regions with long reads (>10 000 bp) to fully characterise all non-coding regions for structural variation. In a family severely affected with early-onset bilateral breast cancer and with negative (normal) results by gene panel and exome sequencing, we identified an intronic SINE-VNTR-Alu retrotransposon insertion that led to the creation of a pseudoexon in the BRCA1 message and introduced a premature truncation. This combination of CRISPR-Cas9 excision and long-read sequencing reveals a class of complex, damaging and otherwise cryptic mutations that may be particularly frequent in tumour suppressor genes replete with intronic repeats.

Characterization of splice-altering mutations in inherited predisposition to cancer.

Mutations responsible for inherited disease may act by disrupting normal transcriptional splicing. Such mutations can be difficult to detect, and their effects difficult to characterize, because many lie deep within exons or introns where they may alter splice enhancers or silencers or introduce new splice acceptors or donors. Multiple mutation-specific and genome-wide approaches have been developed to evaluate these classes of mutations. We introduce a complementary experimental approach, cBROCA, which yields qualitative and quantitative assessments of the effects of genomic mutations on transcriptional splicing of tumor suppressor genes. cBROCA analysis is undertaken by deriving complementary DNA (cDNA) from puromycin-treated patient lymphoblasts, hybridizing the cDNA to the BROCA panel of tumor suppressor genes, and then multiplex sequencing to very high coverage. At each splice junction suggested by split sequencing reads, read depths of test and control samples are compared. Significant Z scores indicate altered transcripts, over and above naturally occurring minor transcripts, and comparisons of read depths indicate relative abundances of mutant and normal transcripts. BROCA analysis of genomic DNA suggested 120 rare mutations from 150 families with cancers of the breast, ovary, uterus, or colon, in >600 informative genotyped relatives. cBROCA analysis of their transcripts revealed a wide variety of consequences of abnormal splicing in tumor suppressor genes, including whole or partial exon skipping, exonification of intronic sequence, loss or gain of exonic and intronic splicing enhancers and silencers, complete intron retention, hypomorphic alleles, and combinations of these alterations. Combined with pedigree analysis, cBROCA sequencing contributes to understanding the clinical consequences of rare inherited mutations.

Inherited predisposition to malignant mesothelioma and overall survival following platinum chemotherapy.

Survival from malignant mesothelioma, particularly pleural mesothelioma, is very poor. For patients with breast, ovarian, or prostate cancers, overall survival is associated with increased sensitivity to platinum chemotherapy due to loss-of-function mutations in DNA repair genes. The goal of this project was to evaluate, in patients with malignant mesothelioma, the relationship between inherited loss-of-function mutations in DNA repair and other tumor suppressor genes and overall survival following platinum chemotherapy. Patients with histologically confirmed malignant mesothelioma were evaluated for inherited mutations in tumor suppressor genes. Survival was evaluated with respect to genotype and site of mesothelioma. Among 385 patients treated with platinum chemotherapy, median overall survival was significantly longer for patients with loss-of-function mutations in any of the targeted genes compared with patients with no such mutation (P = 0.0006). The effect of genotype was highly significant for patients with pleural mesothelioma (median survival 7.9 y versus 2.4 y, P = 0.0012), but not for patients with peritoneal mesothelioma (median survival 8.2 y versus 5.4 y, P = 0.47). Effect of patient genotype on overall survival, measured at 3 y, remained independently significant after adjusting for gender and age at diagnosis, two other known prognostic factors. Patients with pleural mesothelioma with inherited mutations in DNA repair and other tumor suppressor genes appear to particularly benefit from platinum chemotherapy compared with patients without inherited mutations. These patients may also benefit from other DNA repair targeted therapies such as poly-ADP ribose polymerase (PARP) inhibitors.

Inherited thrombocytopenia associated with mutation of UDP-galactose-4-epimerase (GALE).

Severe thrombocytopenia, characterized by dysplastic megakaryocytes and intracranial bleeding, was diagnosed in six individuals from a consanguineous kindred. Three of the individuals were successfully treated by bone marrow transplant. Whole-exome sequencing and homozygosity mapping of multiple family members, coupled with whole-genome sequencing to reveal shared non-coding variants, revealed one potentially functional variant segregating with thrombocytopenia under a recessive model: GALE p.R51W (c.C151T, NM_001127621). The mutation is extremely rare (allele frequency = 2.5 × 10-05), and the likelihood of the observed co-segregation occurring by chance is 1.2 × 10-06. GALE encodes UDP-galactose-4-epimerase, an enzyme of galactose metabolism and glycosylation responsible for two reversible reactions: interconversion of UDP-galactose with UDP-glucose and interconversion of UDP-N-acetylgalactosamine with UDP-N-acetylglucosamine. The mutation alters an amino acid residue that is conserved from yeast to humans. The variant protein has both significantly lower enzymatic activity for both interconversion reactions and highly significant thermal instability. Proper glycosylation is critical to normal hematopoiesis, in particular to megakaryocyte and platelet development, as reflected in the presence of thrombocytopenia in the context of congenital disorders of glycosylation. Mutations in GALE have not previously been associated with thrombocytopenia. Our results suggest that GALE p.R51W is inadequate for normal glycosylation and thereby may impair megakaryocyte and platelet development. If other mutations in GALE are shown to have similar consequences, this gene may be proven to play a critical role in hematopoiesis.

Using Somatic Mutations from Tumors to Classify Variants in Mismatch Repair Genes.

Present guidelines for classification of constitutional variants do not incorporate inferences from mutations seen in tumors, even when these are associated with a specific molecular phenotype. When somatic mutations and constitutional mutations lead to the same molecular phenotype, as for the mismatch repair genes, information from somatic mutations may enable interpretation of previously unclassified variants. To test this idea, we first estimated likelihoods that somatic variants in MLH1, MSH2, MSH6, and PMS2 drive microsatellite instability and characteristic IHC staining patterns by calculating likelihoods of high versus low normalized variant read fractions of 153 mutations known to be pathogenic versus those of 760 intronic passenger mutations from 174 paired tumor-normal samples. Mutations that explained the tumor mismatch repair phenotype had likelihood ratio for high variant read fraction of 1.56 (95% CI 1.42-1.71) at sites with no loss of heterozygosity and of 26.5 (95% CI 13.2-53.0) at sites with loss of heterozygosity. Next, we applied these ratios to 165 missense, synonymous, and splice variants observed in tumors, combining in a Bayesian analysis the likelihood ratio corresponding with the adjusted variant read fraction with pretest probabilities derived from published analyses and public databases. We suggest classifications for 86 of 165 variants: 7 benign, 31 likely benign, 22 likely pathogenic, and 26 pathogenic. These results illustrate that for mismatch repair genes, characterization of tumor mutations permits tumor mutation data to inform constitutional variant classification. We suggest modifications to incorporate molecular phenotype in future variant classification guidelines.

Essential Role of BRCA2 in Ovarian Development and Function.

The causes of ovarian dysgenesis remain incompletely understood. Two sisters with XX ovarian dysgenesis carried compound heterozygous truncating mutations in the BRCA2 gene that led to reduced BRCA2 protein levels and an impaired response to DNA damage, which resulted in chromosomal breakage and the failure of RAD51 to be recruited to double-stranded DNA breaks. The sisters also had microcephaly, and one sister was in long-term remission from leukemia, which had been diagnosed when she was 5 years old. Drosophila mutants that were null for an orthologue of BRCA2 were sterile, and gonadal dysgenesis was present in both sexes. These results revealed a new role for BRCA2 and highlight the importance to ovarian development of genes that are critical for recombination during meiosis. (Funded by the Israel Science Foundation and others.).

De novo mutation in RING1 with epigenetic effects on neurodevelopment.

RING1 is an E3-ubiquitin ligase that is involved in epigenetic control of transcription during development. It is a component of the polycomb repressive complex 1, and its role in that complex is to ubiquitylate histone H2A. In a 13-year-old girl with syndromic neurodevelopmental disabilities, we identified a de novo mutation, RING1 p.R95Q, which alters a conserved arginine residue in the catalytic RING domain. In vitro assays demonstrated that the mutant RING1 retains capacity to catalyze ubiquitin chain formation, but is defective in its ability to ubiquitylate histone H2A in nucleosomes. Consistent with this in vitro effect, cells of the patient showed decreased monoubiquitylation of histone H2A. We modeled the mutant RING1 in Caenorhabditis elegans by editing the comparable amino acid change into spat-3, the suggested RING1 ortholog. Animals with either the missense mutation or complete knockout of spat-3 were defective in monoubiquitylation of histone H2A and had defects in neuronal migration and axon guidance. Relevant to our patient, animals heterozygous for either the missense or knockout allele also showed neuronal defects. Our results support three conclusions: mutation of RING1 is the likely cause of a human neurodevelopmental syndrome, mutation of RING1 can disrupt histone H2A ubiquitylation without disrupting RING1 catalytic activity, and the comparable mutation in C. elegans spat-3 both recapitulates the effects on histone H2A ubiquitylation and leads to neurodevelopmental abnormalities. This role for RING1 adds to our understanding of the importance of aberrant epigenetic effects as causes of human neurodevelopmental disorders.