Physician and informal care use explained by the Pediatric Quality of Life Inventory (PedsQL) in children with suspected genetic disorders.

PURPOSE: To examine associations between Pediatric Quality of Life Inventory (PedsQL) 4.0 Generic Core Scales and PedsQL Infant Scales with formal health care resource utilization (HCRU) and informal caregiver burden. METHODS: We studied a pediatric cohort of 837 patients (median age: 8.4 years) with suspected genetic disorders enrolled January 2019 through July 2021 in the NYCKidSeq program for diagnostic sequencing. Using linked ~ nine-month longitudinal survey and physician claims data collected through May 2022, we modeled the association between baseline PedsQL scores and post-baseline HCRU (median follow-up: 21.1 months) and informal care. We also assessed the longitudinal change in PedsQL scores with physician services using linear mixed-effects models. RESULTS: Lower PedsQL total and physical health scores were independently associated with increases in 18-month physician services, encounters, and weekly informal care. Comparing low vs. median total scores, increases were 10.6 services (95% CI: 1.0-24.6), 3.3 encounters (95% CI: 0.5-6.8), and $668 (95% CI: $350-965), respectively. For the psychosocial domain, higher scores were associated with decreased informal care. Based on adjusted linear mixed-effects modeling, every additional ten physician services was associated with diminished improvement in longitudinal PedsQL total score trajectories by 1.1 point (95% confidence interval: 0.6-1.6) on average. Similar trends were observed in the physical and psychosocial domains. CONCLUSION: PedsQL scores were independently associated with higher utilization of physician services and informal care. Moreover, longitudinal trajectories of PedsQL scores became less favorable with increased physician services. Adding PedsQL survey instruments to conventional measures for improved risk stratification should be evaluated in further research.

The Pediatric Quality of Life Inventory (PedsQL) is widely used to measure health-related quality of life in pediatric patients; however, few studies have examined whether the PedsQL is indicative of longitudinal outcomes of morbidity and health care needs. This study captures associations between PedsQL scores with utilization of physician and informal care in children with suspected genetic disorders. We demonstrate that lower PedsQL total and physical health scores are independently associated with greater utilization of physician services and informal care. Moreover, longitudinal trajectories of PedsQL scores become less favorable with increased physician services. Results can inform future applications of PedsQL instruments.

Genomic and phenotypic correlates of mosaic loss of chromosome Y in blood.

Mosaic loss of Y (mLOY) is the most common somatic chromosomal alteration detected in human blood. The presence of mLOY is associated with altered blood cell counts and increased risk of Alzheimer's disease, solid tumors, and other age-related diseases. We sought to gain a better understanding of genetic drivers and associated phenotypes of mLOY through analyses of whole genome sequencing of a large set of genetically diverse males from the Trans-Omics for Precision Medicine (TOPMed) program. This approach enabled us to identify differences in mLOY frequencies across populations defined by genetic similarity, revealing a higher frequency of mLOY in the European American (EA) ancestry group compared to those of Hispanic American (HA), African American (AA), and East Asian (EAS) ancestry. Further, we identified two genes ( CFHR1 and LRP6 ) that harbor multiple rare, putatively deleterious variants associated with mLOY susceptibility, show that subsets of human hematopoietic stem cells are enriched for activity of mLOY susceptibility variants, and that certain alleles on chromosome Y are more likely to be lost than others.

Managing differential performance of polygenic risk scores across groups: Real-world experience of the eMERGE Network.

The differential performance of polygenic risk scores (PRSs) by group is one of the major ethical barriers to their clinical use. It is also one of the main practical challenges for any implementation effort. The social repercussions of how people are grouped in PRS research must be considered in communications with research participants, including return of results. Here, we outline the decisions faced and choices made by a large multi-site clinical implementation study returning PRSs to diverse participants in handling this issue of differential performance. Our approach to managing the complexities associated with the differential performance of PRSs serves as a case study that can help future implementers of PRSs to plot an anticipatory course in response to this issue.

Beyond the Human Genome Project: The Age of Complete Human Genome Sequences and Pangenome References.

The Human Genome Project was an enormous accomplishment, providing a foundation for countless explorations into the genetics and genomics of the human species. Yet for many years, the human genome reference sequence remained incomplete and lacked representation of human genetic diversity. Recently, two major advances have emerged to address these shortcomings: complete gap-free human genome sequences, such as the one developed by the Telomere-to-Telomere Consortium, and high-quality pangenomes, such as the one developed by the Human Pangenome Reference Consortium. Facilitated by advances in long-read DNA sequencing and genome assembly algorithms, complete human genome sequences resolve regions that have been historically difficult to sequence, including centromeres, telomeres, and segmental duplications. In parallel, pangenomes capture the extensive genetic diversity across populations worldwide. Together, these advances usher in a new era of genomics research, enhancing the accuracy of genomic analysis, paving the path for precision medicine, and contributing to deeper insights into human biology.

MUSSEL: Enhanced Bayesian polygenic risk prediction leveraging information across multiple ancestry groups.

Polygenic risk scores (PRSs) are now showing promising predictive performance on a wide variety of complex traits and diseases, but there exists a substantial performance gap across populations. We propose MUSSEL, a method for ancestry-specific polygenic prediction that borrows information in summary statistics from genome-wide association studies (GWASs) across multiple ancestry groups via Bayesian hierarchical modeling and ensemble learning. In our simulation studies and data analyses across four distinct studies, totaling 5.7 million participants with a substantial ancestral diversity, MUSSEL shows promising performance compared to alternatives. For example, MUSSEL has an average gain in prediction R2 across 11 continuous traits of 40.2% and 49.3% compared to PRS-CSx and CT-SLEB, respectively, in the African ancestry population. The best-performing method, however, varies by GWAS sample size, target ancestry, trait architecture, and linkage disequilibrium reference samples; thus, ultimately a combination of methods may be needed to generate the most robust PRSs across diverse populations.

Principles and methods for transferring polygenic risk scores across global populations.

Polygenic risk scores (PRSs) summarize the genetic predisposition of a complex human trait or disease and may become a valuable tool for advancing precision medicine. However, PRSs that are developed in populations of predominantly European genetic ancestries can increase health disparities due to poor predictive performance in individuals of diverse and complex genetic ancestries. We describe genetic and modifiable risk factors that limit the transferability of PRSs across populations and review the strengths and weaknesses of existing PRS construction methods for diverse ancestries. Developing PRSs that benefit global populations in research and clinical settings provides an opportunity for innovation and is essential for health equity.

Physician services and costs after disclosure of diagnostic sequencing results in the NYCKidSeq program.

PURPOSE: To better understand the effects of returning diagnostic sequencing results on clinical actions and economic outcomes for pediatric patients with suspected genetic disorders. METHODS: Longitudinal physician claims data after diagnostic sequencing were obtained for patients aged 0 to 21 years with neurologic, cardiac, and immunologic disorders with suspected genetic etiology. We assessed specialist consultation rates prompted by primary diagnostic results, as well as marginal effects on overall 18-month physician services and costs. RESULTS: We included data on 857 patients (median age: 9.6 years) with a median follow-up of 17.3 months after disclosure of diagnostic sequencing results. The likelihood of having ≥1 recommendation for specialist consultation in 155 patients with positive findings was high (72%) vs 23% in 443 patients with uncertain findings and 21% in 259 patients with negative findings (P < .001). Follow-through consultation occurred in 30%. Increases in 18-month physician services and costs following a positive finding diminished after multivariable adjustment. Also, no significant differences between those with uncertain and negative findings were demonstrated. CONCLUSION: Our study did not provide evidence for significant increases in downstream physician services and costs after returning positive or uncertain diagnostic sequencing findings. More large-scale longitudinal studies are needed to confirm these findings.

The NYCKidSeq randomized controlled trial: Impact of GUÍA digitally enhanced genetic results disclosure in diverse families.

Digital solutions are needed to support rapid increases in the application of genetic/genomic tests (GTs) in diverse clinical settings and patient populations. We developed GUÍA, a bilingual digital application that facilitates disclosure of GT results. The NYCKidSeq randomized controlled trial enrolled diverse children with neurologic, cardiac, and immunologic conditions who underwent GTs. The trial evaluated GUÍA's impact on understanding the GT results by randomizing families to results disclosure genetic counseling with GUÍA (intervention) or standard of care (SOC). Parents/legal guardians (participants) completed surveys at baseline, post-results disclosure, and 6 months later. Survey measures assessed the primary study outcomes of participants' perceived understanding of and confidence in explaining their child's GT results and the secondary outcome of objective understanding. The analysis included 551 diverse participants, 270 in the GUÍA arm and 281 in SOC. Participants in the GUÍA arm had significantly higher perceived understanding post-results (OR = 2.8, CI[1.004, 7.617], p = 0.049) and maintained higher objective understanding over time (OR = 1.1, CI[1.004, 1.127], p = 0.038) compared to SOC. There was no impact on perceived confidence. Hispanic/Latino(a) individuals in the GUÍA arm maintained higher perceived understanding (OR = 3.9, CI[1.603, 9.254], p = 0.003), confidence (OR = 2.7, CI[1.021, 7.277], p = 0.046), and objective understanding (OR = 1.1, CI[1.009, 1.212], p = 0.032) compared to SOC. This trial demonstrates that GUÍA positively impacts understanding of GT results in diverse parents of children with suspected genetic conditions and builds a case for utilizing GUÍA to deliver complex results. Continued development and evaluation of digital applications in diverse populations are critical for equitably scaling GT offerings in specialty clinics.

Estimating heritability explained by local ancestry and evaluating stratification bias in admixture mapping from summary statistics.

The heritability explained by local ancestry markers in an admixed population (hγ2) provides crucial insight into the genetic architecture of a complex disease or trait. Estimation of hγ2 can be susceptible to biases due to population structure in ancestral populations. Here, we present heritability estimation from admixture mapping summary statistics (HAMSTA), an approach that uses summary statistics from admixture mapping to infer heritability explained by local ancestry while adjusting for biases due to ancestral stratification. Through extensive simulations, we demonstrate that HAMSTA hγ2 estimates are approximately unbiased and are robust to ancestral stratification compared to existing approaches. In the presence of ancestral stratification, we show a HAMSTA-derived sampling scheme provides a calibrated family-wise error rate (FWER) of ∼5% for admixture mapping, unlike existing FWER estimation approaches. We apply HAMSTA to 20 quantitative phenotypes of up to 15,988 self-reported African American individuals in the Population Architecture using Genomics and Epidemiology (PAGE) study. We observe hˆγ2 in the 20 phenotypes range from 0.0025 to 0.033 (mean hˆγ2 = 0.012 ± 9.2 × 10-4), which translates to hˆ2 ranging from 0.062 to 0.85 (mean hˆ2 = 0.30 ± 0.023). Across these phenotypes we find little evidence of inflation due to ancestral population stratification in current admixture mapping studies (mean inflation factor of 0.99 ± 0.001). Overall, HAMSTA provides a fast and powerful approach to estimate genome-wide heritability and evaluate biases in test statistics of admixture mapping studies.

An Ethical Framework for Research Using Genetic Ancestry.

A wide range of research uses patterns of genetic variation to infer genetic similarity between individuals, typically referred to as genetic ancestry. This research includes inference of human demographic history, understanding the genetic architecture of traits, and predicting disease risk. Researchers are not just structuring an intellectual inquiry when using genetic ancestry, they are also creating analytical frameworks with broader societal ramifications. This essay presents an ethics framework in the spirit of virtue ethics for these researchers: rather than focus on rule following, the framework is designed to build researchers' capacities to react to the ethical dimensions of their work. The authors identify one overarching principle of intellectual freedom and responsibility, noting that freedom in all its guises comes with responsibility, and they identify and define four principles that collectively uphold researchers' intellectual responsibility: truthfulness, justice and fairness, anti-racism, and public beneficence. Researchers should bring their practices into alignment with these principles, and to aid this, the authors name three common ways research practices infringe these principles, suggest a step-by-step process for aligning research choices with the principles, provide rules of thumb for achieving alignment, and give a worked case. The essay concludes by identifying support needed by researchers to act in accord with the proposed framework.

The genetic underpinnings of variable penetrance and expressivity of pathogenic mutations in cardiometabolic traits.

Over three percent of people carry a dominant pathogenic mutation, yet only a fraction of carriers develop disease (incomplete penetrance), and phenotypes from mutations in the same gene range from mild to severe (variable expressivity). Here, we investigate underlying mechanisms for this heterogeneity: variable variant effect sizes, carrier polygenic backgrounds, and modulation of carrier effect by genetic background (epistasis). We leveraged exomes and clinical phenotypes from the UK Biobank and the Mt. Sinai Bio Me Biobank to identify carriers of pathogenic variants affecting cardiometabolic traits. We employed recently developed methods to study these cohorts, observing strong statistical support and clinical translational potential for all three mechanisms of variable penetrance and expressivity. For example, scores from our recent model of variant pathogenicity were tightly correlated with phenotype amongst clinical variant carriers, they predicted effects of variants of unknown significance, and they distinguished gain- from loss-of-function variants. We also found that polygenic scores predicted phenotypes amongst pathogenic carriers and that epistatic effects can exceed main carrier effects by an order of magnitude.

Assessing efficiency of fine-mapping obesity-associated variants through leveraging ancestry architecture and functional annotation using PAGE and UKBB cohorts.

Inadequate representation of non-European ancestry populations in genome-wide association studies (GWAS) has limited opportunities to isolate functional variants. Fine-mapping in multi-ancestry populations should improve the efficiency of prioritizing variants for functional interrogation. To evaluate this hypothesis, we leveraged ancestry architecture to perform comparative GWAS and fine-mapping of obesity-related phenotypes in European ancestry populations from the UK Biobank (UKBB) and multi-ancestry samples from the Population Architecture for Genetic Epidemiology (PAGE) consortium with comparable sample sizes. In the investigated regions with genome-wide significant associations for obesity-related traits, fine-mapping in our ancestrally diverse sample led to 95% and 99% credible sets (CS) with fewer variants than in the European ancestry sample. Lead fine-mapped variants in PAGE regions had higher average coding scores, and higher average posterior probabilities for causality compared to UKBB. Importantly, 99% CS in PAGE loci contained strong expression quantitative trait loci (eQTLs) in adipose tissues or harbored more variants in tighter linkage disequilibrium (LD) with eQTLs. Leveraging ancestrally diverse populations with heterogeneous ancestry architectures, coupled with functional annotation, increased fine-mapping efficiency and performance, and reduced the set of candidate variants for consideration for future functional studies. Significant overlap in genetic causal variants across populations suggests generalizability of genetic mechanisms underpinning obesity-related traits across populations.

Admixture mapping of peripheral artery disease in a Dominican population reveals a putative risk locus on 2q35.

Peripheral artery disease (PAD) is a form of atherosclerotic cardiovascular disease, affecting ∼8 million Americans, and is known to have racial and ethnic disparities. PAD has been reported to have a significantly higher prevalence in African Americans (AAs) compared to non-Hispanic European Americans (EAs). Hispanic/Latinos (HLs) have been reported to have lower or similar rates of PAD compared to EAs, despite having a paradoxically high burden of PAD risk factors; however, recent work suggests prevalence may differ between sub-groups. Here, we examined a large cohort of diverse adults in the BioMe biobank in New York City. We observed the prevalence of PAD at 1.7% in EAs vs. 8.5% and 9.4% in AAs and HLs, respectively, and among HL sub-groups, the prevalence was found at 11.4% and 11.5% in Puerto Rican and Dominican populations, respectively. Follow-up analysis that adjusted for common risk factors demonstrated that Dominicans had the highest increased risk for PAD relative to EAs [OR = 3.15 (95% CI 2.33-4.25), p < 6.44 × 10-14]. To investigate whether genetic factors may explain this increased risk, we performed admixture mapping by testing the association between local ancestry and PAD in Dominican BioMe participants (N = 1,813) separately from European, African, and Native American (NAT) continental ancestry tracts. The top association with PAD was an NAT ancestry tract at chromosome 2q35 [OR = 1.96 (SE = 0.16), p < 2.75 × 10-05) with 22.6% vs. 12.9% PAD prevalence in heterozygous NAT tract carriers versus non-carriers, respectively. Fine-mapping at this locus implicated tag SNP rs78529201 located within a long intergenic non-coding RNA (lincRNA) LINC00607, a gene expression regulator of key genes related to thrombosis and extracellular remodeling of endothelial cells, suggesting a putative link of the 2q35 locus to PAD etiology. Efforts to reproduce the signal in other Hispanic cohorts were unsuccessful. In summary, we showed how leveraging health system data helped understand nuances of PAD risk across HL sub-groups and admixture mapping approaches elucidated a putative risk locus in a Dominican population.

Carriers of rare damaging CCR2 genetic variants are at lower risk of atherosclerotic disease.

Background: The CCL2/CCR2 axis governs monocyte trafficking and recruitment to atherosclerotic lesions. Human genetic analyses and population-based studies support an association between circulating CCL2 levels and atherosclerosis. Still, it remains unknown whether pharmacological targeting of CCR2, the main CCL2 receptor, would provide protection against human atherosclerotic disease. Methods: In whole-exome sequencing data from 454,775 UK Biobank participants (40-69 years), we identified predicted loss-of-function (LoF) or damaging missense (REVEL score >0.5) variants within the CCR2 gene. We prioritized variants associated with lower monocyte count (p<0.05) and tested associations with vascular risk factors and risk of atherosclerotic disease over a mean follow-up of 14 years. The results were replicated in a pooled cohort of three independent datasets (TOPMed, deCODE and Penn Medicine BioBank; total n=441,445) and the effect of the most frequent damaging variant was experimentally validated. Results: A total of 45 predicted LoF or damaging missense variants were identified in the CCR2 gene, 4 of which were also significantly associated with lower monocyte count, but not with other white blood cell counts. Heterozygous carriers of these variants were at a lower risk of a combined atherosclerosis outcome, showed a lower burden of atherosclerosis across four vascular beds, and were at a lower lifetime risk of coronary artery disease and myocardial infarction. There was no evidence of association with vascular risk factors including LDL-cholesterol, blood pressure, glycemic status, or C-reactive protein. Using a cAMP assay, we found that cells transfected with the most frequent CCR2 damaging variant (3:46358273:T:A, M249K, 547 carriers, frequency: 0.14%) show a decrease in signaling in response to CCL2. The associations of the M249K variant with myocardial infarction were consistent across cohorts (ORUKB: 0.62 95%CI: 0.39-0.96; ORexternal: 0.64 95%CI: 0.34-1.19; ORpooled: 0.64 95%CI: 0.450.90). In a phenome-wide association study, we found no evidence for higher risk of common infections or mortality among carriers of damaging CCR2 variants. Conclusions: Heterozygous carriers of damaging CCR2 variants have a lower burden of atherosclerosis and lower lifetime risk of myocardial infarction. In conjunction with previous evidence from experimental and epidemiological studies, our findings highlight the translational potential of CCR2-targeting as an atheroprotective approach.

The NYCKidSeq randomized controlled trial: Impact of GUÍA digitally enhanced genetic counseling in racially and ethnically diverse families.

Background: Digital solutions are needed to support rapid increases in the application of genetic and genomic tests (GT) in diverse clinical settings and patient populations. We developed GUÍA, a bi-lingual web-based platform that facilitates disclosure of GT results. The NYCKidSeq randomized controlled trial evaluated GUÍA's impact on understanding of GT results. Methods: NYCKidSeq enrolled diverse children with neurologic, cardiac, and immunologic conditions who underwent GT. Families were randomized to genetic counseling with GUÍA (intervention) or standard of care (SOC) genetic counseling for results disclosure. Parents/legal guardians (participants) completed surveys at baseline, post-results disclosure, and 6-months later. Survey measures assessed the primary study outcomes of perceived understanding of and confidence in explaining their child's GT results and the secondary outcome of objective understanding. We used regression models to evaluate the association between the intervention and the study outcomes. Results: The analysis included 551 participants, 270 in the GUÍA arm and 281 in SOC. Participants' mean age was 41.1 years and 88.6% were mothers. Most participants were Hispanic/Latino(a) (46.3%), White/European American (24.5%), or Black/African American (15.8%). Participants in the GUÍA arm had significantly higher perceived understanding post-results (OR=2.8, CI[1.004,7.617], P=0.049) and maintained higher objective understanding over time (OR=1.1, CI[1.004, 1.127], P=0.038) compared to those in the SOC arm. There was no impact on perceived confidence. Hispanic/Latino(a) individuals in the GUÍA arm maintained higher perceived understanding (OR=3.9, CI[1.6, 9.3], P=0.003), confidence (OR=2.7, CI[1.021, 7.277], P=0.046), and objective understanding (OR=1.1, CI[1.009, 1.212], P=0.032) compared to SOC . Conclusions: This trial demonstrates that GUÍA positively impacts understanding of GT results in diverse parents of children with suspected genetic conditions. These findings build a case for utilizing GUÍA to deliver complex and often ambiguous genetic results. Continued development and evaluation of digital applications in diverse populations are critical for equitably scaling GT offerings in specialty clinics. Trial Registration: Clinicaltrials.gov identifier NCT03738098.

Disease risk and healthcare utilization among ancestrally diverse groups in the Los Angeles region.

An individual's disease risk is affected by the populations that they belong to, due to shared genetics and environmental factors. The study of fine-scale populations in clinical care is important for identifying and reducing health disparities and for developing personalized interventions. To assess patterns of clinical diagnoses and healthcare utilization by fine-scale populations, we leveraged genetic data and electronic medical records from 35,968 patients as part of the UCLA ATLAS Community Health Initiative. We defined clusters of individuals using identity by descent, a form of genetic relatedness that utilizes shared genomic segments arising due to a common ancestor. In total, we identified 376 clusters, including clusters with patients of Afro-Caribbean, Puerto Rican, Lebanese Christian, Iranian Jewish and Gujarati ancestry. Our analysis uncovered 1,218 significant associations between disease diagnoses and clusters and 124 significant associations with specialty visits. We also examined the distribution of pathogenic alleles and found 189 significant alleles at elevated frequency in particular clusters, including many that are not regularly included in population screening efforts. Overall, this work progresses the understanding of health in understudied communities and can provide the foundation for further study into health inequities.

Identification of copy number variants with genome sequencing: Clinical experiences from the NYCKidSeq program.

Copy number variations (CNVs) play a significant role in human disease. While chromosomal microarray has traditionally been the first-tier test for CNV detection, use of genome sequencing (GS) is increasing. We report the frequency of CNVs detected with GS in a diverse pediatric cohort from the NYCKidSeq program and highlight specific examples of its clinical impact. A total of 1052 children (0-21 years) with neurodevelopmental, cardiac, and/or immunodeficiency phenotypes received GS. Phenotype-driven analysis was used, resulting in 183 (17.4%) participants with a diagnostic result. CNVs accounted for 20.2% of participants with a diagnostic result (37/183) and ranged from 0.5 kb to 16 Mb. Of participants with a diagnostic result (n = 183) and phenotypes in more than one category, 5/17 (29.4%) were solved by a CNV finding, suggesting a high prevalence of diagnostic CNVs in participants with complex phenotypes. Thirteen participants with a diagnostic CNV (35.1%) had previously uninformative genetic testing, of which nine included a chromosomal microarray. This study demonstrates the benefits of GS for reliable detection of CNVs in a pediatric cohort with variable phenotypes.

Molecular diagnostic yield of genome sequencing versus targeted gene panel testing in racially and ethnically diverse pediatric patients.

PURPOSE: Adoption of genome sequencing (GS) as a first-line test requires evaluation of its diagnostic yield. We evaluated the GS and targeted gene panel (TGP) testing in diverse pediatric patients (probands) with suspected genetic conditions. METHODS: Probands with neurologic, cardiac, or immunologic conditions were offered GS and TGP testing. Diagnostic yield was compared using a fully paired study design. RESULTS: A total of 645 probands (median age 9 years) underwent genetic testing, and 113 (17.5%) received a molecular diagnosis. Among 642 probands with both GS and TGP testing, GS yielded 106 (16.5%) and TGPs yielded 52 (8.1%) diagnoses (P < .001). Yield was greater for GS vs TGPs in Hispanic/Latino(a) (17.2% vs 9.5%, P < .001) and White/European American (19.8% vs 7.9%, P < .001) but not in Black/African American (11.5% vs 7.7%, P = .22) population groups by self-report. A higher rate of inconclusive results was seen in the Black/African American (63.8%) vs White/European American (47.6%; P = .01) population group. Most causal copy number variants (17 of 19) and mosaic variants (6 of 8) were detected only by GS. CONCLUSION: GS may yield up to twice as many diagnoses in pediatric patients compared with TGP testing but not yet across all population groups.

Poison exon annotations improve the yield of clinically relevant variants in genomic diagnostic testing.

PURPOSE: Neurodevelopmental disorders (NDDs) often result from rare genetic variation, but genomic testing yield for NDDs remains below 50%, suggesting that clinically relevant variants may be missed by standard analyses. Here, we analyze "poison exons" (PEs), which are evolutionarily conserved alternative exons often absent from standard gene annotations. Variants that alter PE inclusion can lead to loss of function and may be highly penetrant contributors to disease. METHODS: We curated published RNA sequencing data from developing mouse cortex to define 1937 conserved PE regions potentially relevant to NDDs, and we analyzed variants found by genome sequencing in multiple NDD cohorts. RESULTS: Across 2999 probands, we found 6 novel clinically relevant variants in PE regions. Five of these variants are in genes that are part of the sodium voltage-gated channel alpha subunit family (SCN1A, SCN2A, and SCN8A), which is associated with epilepsies. One variant is in SNRPB, associated with cerebrocostomandibular syndrome. These variants have moderate to high computational impact assessments, are absent from population variant databases, and in genes with gene-phenotype associations consistent with each probands reported features. CONCLUSION: With a very minimal increase in variant analysis burden (average of 0.77 variants per proband), annotation of PEs can improve diagnostic yield for NDDs and likely other congenital conditions.

TOPORS as a novel causal gene for Joubert syndrome.

Joubert syndrome (JBTS) is a Mendelian disorder of the primary cilium defined by the clinical triad of hypotonia, developmental delay, and a distinct cerebellar malformation called the molar tooth sign. JBTS is inherited in an autosomal recessive, autosomal dominant, or X-linked recessive manner. Though over 40 genes have been identified as causal for JBTS, molecular diagnosis is not made in 30%-40% of individuals who meet clinical criteria. TOPORS encodes topoisomerase I-binding arginine/serine-rich protein, and homozygosity for a TOPORS missense variant (c.29C > A; p.(Pro10Gln)) was identified in individuals with the ciliopathy oral-facial-digital syndrome in two families of Dominican descent. Here, we report an additional proband of Dominican ancestry with JBTS found by exome sequencing to be homozygous for the identical p.(Pro10Gln) TOPORS missense variant. Query of the Mount Sinai BioMe biobank, which includes 1880 individuals of Dominican ancestry, supports a high carrier frequency of the TOPORS p.(Pro10Gln) variant in individuals of Dominican descent. Our data nominates TOPORS as a novel causal gene for JBTS and suggests that TOPORS variants should be considered in the differential of ciliopathy-spectrum disease in individuals of Dominican ancestry.