Exome sequencing reveals genetic heterogeneity and clinically actionable findings in children with cerebral palsy.

Cerebral palsy (CP) is the most common motor disability in children. To ascertain the role of major genetic variants in the etiology of CP, we conducted exome sequencing on a large-scale cohort with clinical manifestations of CP. The study cohort comprised 505 girls and 1,073 boys. Utilizing the current gold standard in genetic diagnostics, 387 of these 1,578 children (24.5%) received genetic diagnoses. We identified 412 pathogenic and likely pathogenic (P/LP) variants across 219 genes associated with neurodevelopmental disorders, and 59 P/LP copy number variants. The genetic diagnostic rate of children with CP labeled at birth with perinatal asphyxia was higher than the rate in children without asphyxia (P = 0.0033). Also, 33 children with CP manifestations (8.5%, 33 of 387) had findings that were clinically actionable. These results highlight the need for early genetic testing in children with CP, especially those with risk factors like perinatal asphyxia, to enable evidence-based medical decision-making.

Tumor Antigens beyond the Human Exome.

With the advent of immunotherapeutics, a new era in the combat against cancer has begun. Particularly promising are neo-epitope-targeted therapies as the expression of neo-antigens is tumor-specific. In turn, this allows the selective targeting and killing of cancer cells whilst healthy cells remain largely unaffected. So far, many advances have been made in the development of treatment options which are tailored to the individual neo-epitope repertoire. The next big step is the achievement of efficacious "off-the-shelf" immunotherapies. For this, shared neo-epitopes propose an optimal target. Given the tremendous potential, a thorough understanding of the underlying mechanisms which lead to the formation of neo-antigens is of fundamental importance. Here, we review the various processes which result in the formation of neo-epitopes. Broadly, the origin of neo-epitopes can be categorized into three groups: canonical, noncanonical, and viral neo-epitopes. For the canonical neo-antigens that arise in direct consequence of somatic mutations, we summarize past and recent findings. Beyond that, our main focus is put on the discussion of noncanonical and viral neo-epitopes as we believe that targeting those provides an encouraging perspective to shape the future of cancer immunotherapeutics.

Whole-exome-based single nucleotide variants and copy number analysis for prenatal diagnosis of compound heterozygosity of SMPD4.

BACKGROUND: Biallelic loss-of-function variants in SMPD4 cause a rare and severe neurodevelopmental disorder. These variants have been identified in a group of children with neurodevelopmental disorders with microcephaly, arthrogryposis, and structural brain anomalies. SMPD4 encodes a sphingomyelinase that hydrolyzes sphingomyelin into ceramide at neutral pH and can thereby affect membrane lipid homeostasis. SMPD4 localizes to the membranes of the endoplasmic reticulum and nuclear envelope and interacts with nuclear pore complexes. MATERIALS AND METHODS: For the efficient prenatal diagnosis of rare and undiagnosed diseases, the parallel detection of copy number variants (CNVs) and single nucleotide variants using whole-exome analysis is required. A physical examination of the parents was performed. Karyotype and whole-exome analysis were performed for the fetus and the parents. RESULTS: A fetus with microcephaly and arthrogryposis; biallelic null variants (c.387-1G>A; Chr2[GRCh38]: g.130142742_130202459del) were detected by whole-exome sequencing (WES). We have reported for the first time the biallelic loss-of-function mutations in SMPD4 in patients born to unrelated parents in China. CONCLUSION: WES could replace chromosomal microarray analysis and copy number variation sequencing as a more cost-effective genetic test for detecting CNVs and diagnosing highly heterogeneous conditions.

Exome sequencing implicates ancestry-related Mendelian variation at SYNE1 in childhood-onset essential hypertension.

Childhood-onset essential hypertension (COEH) is an uncommon form of hypertension that manifests in childhood or adolescence and, in the United States, disproportionately affects children of African ancestry. The etiology of COEH is unknown, but its childhood onset, low prevalence, high heritability, and skewed ancestral demography suggest the potential to identify rare genetic variation segregating in a Mendelian manner among affected individuals and thereby implicate genes important to disease pathogenesis. However, no COEH genes have been reported to date. Here, we identify recessive segregation of rare and putatively damaging missense variation in the spectrin domain of spectrin repeat containing nuclear envelope protein 1 (SYNE1), a cardiovascular candidate gene, in 3 of 16 families with early-onset COEH without an antecedent family history. By leveraging exome sequence data from an additional 48 COEH families, 1,700 in-house trios, and publicly available data sets, we demonstrate that compound heterozygous SYNE1 variation in these COEH individuals occurred more often than expected by chance and that this class of biallelic rare variation was significantly enriched among individuals of African genetic ancestry. Using in vitro shRNA knockdown of SYNE1, we show that reduced SYNE1 expression resulted in a substantial decrease in the elasticity of smooth muscle vascular cells that could be rescued by pharmacological inhibition of the downstream RhoA/Rho-associated protein kinase pathway. These results provide insights into the molecular genetics and underlying pathophysiology of COEH and suggest a role for precision therapeutics in the future.

Exome sequencing in four families with neurodevelopmental disorders: genotype-phenotype correlation and identification of novel disease-causing variants in VPS13B and RELN.

Neurodevelopmental disorders (NDDs) are a clinically and genetically heterogeneous group of early-onset pediatric disorders that affect the structure and/or function of the central or peripheral nervous system. Achieving a precise molecular diagnosis for NDDs may be challenging due to the diverse genetic underpinnings and clinical variability. In the current study, we investigated the underlying genetic cause(s) of NDDs in four unrelated Pakistani families. Using exome sequencing (ES) as a diagnostic approach, we identified disease-causing variants in established NDD-associated genes in all families, including one hitherto unreported variant in RELN and three recurrent variants in VPS13B, DEGS1, and SPG11. Overall, our study highlights the potential of ES as a tool for clinical diagnosis.

Identifying proteomic risk factors for cancer using prospective and exome analyses of 1463 circulating proteins and risk of 19 cancers in the UK Biobank.

The availability of protein measurements and whole exome sequence data in the UK Biobank enables investigation of potential observational and genetic protein-cancer risk associations. We investigated associations of 1463 plasma proteins with incidence of 19 cancers and 9 cancer subsites in UK Biobank participants (average 12 years follow-up). Emerging protein-cancer associations were further explored using two genetic approaches, cis-pQTL and exome-wide protein genetic scores (exGS). We identify 618 protein-cancer associations, of which 107 persist for cases diagnosed more than seven years after blood draw, 29 of 618 were associated in genetic analyses, and four had support from long time-to-diagnosis ( > 7 years) and both cis-pQTL and exGS analyses: CD74 and TNFRSF1B with NHL, ADAM8 with leukemia, and SFTPA2 with lung cancer. We present multiple blood protein-cancer risk associations, including many detectable more than seven years before cancer diagnosis and that had concordant evidence from genetic analyses, suggesting a possible role in cancer development.

Decoding triancestral origins, archaic introgression, and natural selection in the Japanese population by whole-genome sequencing.

We generated Japanese Encyclopedia of Whole-Genome/Exome Sequencing Library (JEWEL), a high-depth whole-genome sequencing dataset comprising 3256 individuals from across Japan. Analysis of JEWEL revealed genetic characteristics of the Japanese population that were not discernible using microarray data. First, rare variant-based analysis revealed an unprecedented fine-scale genetic structure. Together with population genetics analysis, the present-day Japanese can be decomposed into three ancestral components. Second, we identified unreported loss-of-function (LoF) variants and observed that for specific genes, LoF variants appeared to be restricted to a more limited set of transcripts than would be expected by chance, with PTPRD as a notable example. Third, we identified 44 archaic segments linked to complex traits, including a Denisovan-derived segment at NKX6-1 associated with type 2 diabetes. Most of these segments are specific to East Asians. Fourth, we identified candidate genetic loci under recent natural selection. Overall, our work provided insights into genetic characteristics of the Japanese population.

Exome-wide analysis implicates rare protein-altering variants in human handedness.

Handedness is a manifestation of brain hemispheric specialization. Left-handedness occurs at increased rates in neurodevelopmental disorders. Genome-wide association studies have identified common genetic effects on handedness or brain asymmetry, which mostly involve variants outside protein-coding regions and may affect gene expression. Implicated genes include several that encode tubulins (microtubule components) or microtubule-associated proteins. Here we examine whether left-handedness is also influenced by rare coding variants (frequencies ≤ 1%), using exome data from 38,043 left-handed and 313,271 right-handed individuals from the UK Biobank. The beta-tubulin gene TUBB4B shows exome-wide significant association, with a rate of rare coding variants 2.7 times higher in left-handers than right-handers. The TUBB4B variants are mostly heterozygous missense changes, but include two frameshifts found only in left-handers. Other TUBB4B variants have been linked to sensorineural and/or ciliopathic disorders, but not the variants found here. Among genes previously implicated in autism or schizophrenia by exome screening, DSCAM and FOXP1 show evidence for rare coding variant association with left-handedness. The exome-wide heritability of left-handedness due to rare coding variants was 0.91%. This study reveals a role for rare, protein-altering variants in left-handedness, providing further evidence for the involvement of microtubules and disorder-relevant genes.

Genetic associations of protein-coding variants in venous thromboembolism.

Previous genetic studies of venous thromboembolism (VTE) have been largely limited to common variants, leaving the genetic determinants relatively incomplete. We performed an exome-wide association study of VTE among 14,723 cases and 334,315 controls. Fourteen known and four novel genes (SRSF6, PHPT1, CGN, and MAP3K2) were identified through protein-coding variants, with broad replication in the FinnGen cohort. Most genes we discovered exhibited the potential to predict future VTE events in longitudinal analysis. Notably, we provide evidence for the additive contribution of rare coding variants to known genome-wide polygenic risk in shaping VTE risk. The identified genes were enriched in pathways affecting coagulation and platelet activation, along with liver-specific expression. The pleiotropic effects of these genes indicated the potential involvement of coagulation factors, blood cell traits, liver function, and immunometabolic processes in VTE pathogenesis. In conclusion, our study unveils the valuable contribution of protein-coding variants in VTE etiology and sheds new light on its risk stratification.

Comprehensive Molecular Profiling of NPM1-Mutated Acute Myeloid Leukemia Using RNAseq Approach.

Acute myeloid leukemia (AML) is a complex hematologic malignancy with high morbidity and mortality. Nucleophosmin 1 (NPM1) mutations occur in approximately 30% of AML cases, and NPM1-mutated AML is classified as a distinct entity. NPM1-mutated AML patients without additional genetic abnormalities have a favorable prognosis. Despite this, 30-50% of them experience relapse. This study aimed to investigate the potential of total RNAseq in improving the characterization of NPM1-mutated AML patients. We explored genetic variations independently of myeloid stratification, revealing a complex molecular scenario. We showed that total RNAseq enables the uncovering of different genetic alterations and clonal subtypes, allowing for a comprehensive evaluation of the real expression of exome transcripts in leukemic clones and the identification of aberrant fusion transcripts. This characterization may enhance understanding and guide improved treatment strategies for NPM1mut AML patients, contributing to better outcomes. Our findings underscore the complexity of NPM1-mutated AML, supporting the incorporation of advanced technologies for precise risk stratification and personalized therapeutic strategies. The study provides a foundation for future investigations into the clinical implications of identified genetic variations and highlights the importance of evolving diagnostic approaches in leukemia management.

Identification of skewed X chromosome inactivation using exome and transcriptome sequencing in patients with suspected rare genetic disease.

BACKGROUND: X-chromosome inactivation (XCI) is an epigenetic process that occurs during early development in mammalian females by randomly silencing one of two copies of the X chromosome in each cell. The preferential inactivation of either the maternal or paternal copy of the X chromosome in a majority of cells results in a skewed or non-random pattern of X inactivation and is observed in over 25% of adult females. Identifying skewed X inactivation is of clinical significance in patients with suspected rare genetic diseases due to the possibility of biased expression of disease-causing genes present on the active X chromosome. The current clinical test for the detection of skewed XCI relies on the methylation status of the methylation-sensitive restriction enzyme (Hpall) binding site present in proximity of short tandem polymorphic repeats on the androgen receptor (AR) gene. This approach using one locus results in uninformative or inconclusive data for 10-20% of tests. Further, recent studies have shown inconsistency between methylation of the AR locus and the state of inactivation of the X chromosome. Herein, we develop a method for estimating X inactivation status, using exome and transcriptome sequencing data derived from blood in 227 female samples. We built a reference model for evaluation of XCI in 135 females from the GTEx consortium. We tested and validated the model on 11 female individuals with different types of undiagnosed rare genetic disorders who were clinically tested for X-skew using the AR gene assay and compared results to our outlier-based analysis technique. RESULTS: In comparison to the AR clinical test for identification of X inactivation, our method was concordant with the AR method in 9 samples, discordant in 1, and provided a measure of X inactivation in 1 sample with uninformative clinical results. We applied this method on an additional 81 females presenting to the clinic with phenotypes consistent with different hereditary disorders without a known genetic diagnosis. CONCLUSIONS: This study presents the use of transcriptome and exome sequencing data to provide an accurate and complete estimation of X-inactivation and skew status in a cohort of female patients with different types of suspected rare genetic disease.

Prioritization of oligogenic variant combinations in whole exomes.

MOTIVATION: Whole exome sequencing (WES) has emerged as a powerful tool for genetic research, enabling the collection of a tremendous amount of data about human genetic variation. However, properly identifying which variants are causative of a genetic disease remains an important challenge, often due to the number of variants that need to be screened. Expanding the screening to combinations of variants in two or more genes, as would be required under the oligogenic inheritance model, simply blows this problem out of proportion. RESULTS: We present here the High-throughput oligogenic prioritizer (Hop), a novel prioritization method that uses direct oligogenic information at the variant, gene and gene pair level to detect digenic variant combinations in WES data. This method leverages information from a knowledge graph, together with specialized pathogenicity predictions in order to effectively rank variant combinations based on how likely they are to explain the patient's phenotype. The performance of Hop is evaluated in cross-validation on 36 120 synthetic exomes for training and 14 280 additional synthetic exomes for independent testing. Whereas the known pathogenic variant combinations are found in the top 20 in approximately 60% of the cross-validation exomes, 71% are found in the same ranking range when considering the independent set. These results provide a significant improvement over alternative approaches that depend simply on a monogenic assessment of pathogenicity, including early attempts for digenic ranking using monogenic pathogenicity scores. AVAILABILITY AND IMPLEMENTATION: Hop is available at https://github.com/oligogenic/HOP.

Clustering of predicted loss-of-function variants in genes linked with monogenic disease can explain incomplete penetrance.

BACKGROUND: Genetic variants that severely alter protein products (e.g. nonsense, frameshift) are often associated with disease. For some genes, these predicted loss-of-function variants (pLoFs) are observed throughout the gene, whilst in others, they occur only at specific locations. We hypothesised that, for genes linked with monogenic diseases that display incomplete penetrance, pLoF variants present in apparently unaffected individuals may be limited to regions where pLoFs are tolerated. To test this, we investigated whether pLoF location could explain instances of incomplete penetrance of variants expected to be pathogenic for Mendelian conditions. METHODS: We used exome sequence data in 454,773 individuals in the UK Biobank (UKB) to investigate the locations of pLoFs in a population cohort. We counted numbers of unique pLoF, missense, and synonymous variants in UKB in each quintile of the coding sequence (CDS) of all protein-coding genes and clustered the variants using Gaussian mixture models. We limited the analyses to genes with ≥ 5 variants of each type (16,473 genes). We compared the locations of pLoFs in UKB with all theoretically possible pLoFs in a transcript, and pathogenic pLoFs from ClinVar, and performed simulations to estimate the false-positive rate of non-uniformly distributed variants. RESULTS: For most genes, all variant classes fell into clusters representing broadly uniform variant distributions, but genes in which haploinsufficiency causes developmental disorders were less likely to have uniform pLoF distribution than other genes (P < 2.2 × 10-6). We identified a number of genes, including ARID1B and GATA6, where pLoF variants in the first quarter of the CDS were rescued by the presence of an alternative translation start site and should not be reported as pathogenic. For other genes, such as ODC1, pLoFs were located approximately uniformly across the gene, but pathogenic pLoFs were clustered only at the end, consistent with a gain-of-function disease mechanism. CONCLUSIONS: Our results suggest the potential benefits of localised constraint metrics and that the location of pLoF variants should be considered when interpreting variants.

Exome copy number variant detection, analysis, and classification in a large cohort of families with undiagnosed rare genetic disease.

Copy number variants (CNVs) are significant contributors to the pathogenicity of rare genetic diseases and, with new innovative methods, can now reliably be identified from exome sequencing. Challenges still remain in accurate classification of CNV pathogenicity. CNV calling using GATK-gCNV was performed on exomes from a cohort of 6,633 families (15,759 individuals) with heterogeneous phenotypes and variable prior genetic testing collected at the Broad Institute Center for Mendelian Genomics of the Genomics Research to Elucidate the Genetics of Rare Diseases consortium and analyzed using the seqr platform. The addition of CNV detection to exome analysis identified causal CNVs for 171 families (2.6%). The estimated sizes of CNVs ranged from 293 bp to 80 Mb. The causal CNVs consisted of 140 deletions, 15 duplications, 3 suspected complex structural variants (SVs), 3 insertions, and 10 complex SVs, the latter two groups being identified by orthogonal confirmation methods. To classify CNV variant pathogenicity, we used the 2020 American College of Medical Genetics and Genomics/ClinGen CNV interpretation standards and developed additional criteria to evaluate allelic and functional data as well as variants on the X chromosome to further advance the framework. We interpreted 151 CNVs as likely pathogenic/pathogenic and 20 CNVs as high-interest variants of uncertain significance. Calling CNVs from existing exome data increases the diagnostic yield for individuals undiagnosed after standard testing approaches, providing a higher-resolution alternative to arrays at a fraction of the cost of genome sequencing. Our improvements to the classification approach advances the systematic framework to assess the pathogenicity of CNVs.

Focused Exome Sequencing Gives a High Diagnostic Yield in the Indian Subcontinent.

The genetically isolated yet heterogeneous and highly consanguineous Indian population has shown a higher prevalence of rare genetic disorders. However, there is a significant socioeconomic burden for genetic testing to be accessible to the general population. In the current study, we analyzed next-generation sequencing data generated through focused exome sequencing from individuals with different phenotypic manifestations referred for genetic testing to achieve a molecular diagnosis. Pathogenic or likely pathogenic variants are reported in 280 of 833 cases with a diagnostic yield of 33.6%. Homozygous sequence and copy number variants were found as positive diagnostic findings in 131 cases (15.7%) because of the high consanguinity in the Indian population. No relevant findings related to reported phenotype were identified in 6.2% of the cases. Patients referred for testing due to metabolic disorder and neuromuscular disorder had higher diagnostic yields. Carrier testing of asymptomatic individuals with a family history of the disease, through focused exome sequencing, achieved positive diagnosis in 54 of 118 cases tested. Copy number variants were also found in trans with single-nucleotide variants and mitochondrial variants in a few of the cases. The diagnostic yield and the findings from this study signify that a focused exome test is a good lower-cost alternative for whole-exome and whole-genome sequencing and as a first-tier approach to genetic testing.

Genetics Visit Uptake Among Individuals Receiving Clinically Actionable Genomic Screening Results.

Importance: Screening unselected populations for clinically actionable genetic disease risk can improve ascertainment and facilitate risk management. Genetics visits may encourage at-risk individuals to perform recommended management, but little has been reported on genetics visit completion or factors associated with completion in genomic screening programs. Objective: To identify factors associated with postdisclosure genetics visits in a genomic screening cohort. Design, Setting, and Participants: This was a cohort study of biobank data in a health care system in central Pennsylvania. Participants' exome sequence data were reviewed for pathogenic or likely pathogenic (P/LP) results in all genes on the American College of Medical Genetics and Genomics Secondary Findings list. Clinically confirmed results were disclosed by phone and letter. Participants included adult MyCode biobank participants who received P/LP results between July 2015 and November 2019. Data were analyzed from May 2021 to March 2022. Exposure: Clinically confirmed P/LP result disclosed by phone or letter. Main Outcomes and Measures: Completion of genetics visit in which the result was discussed and variables associated with completion were assessed by electronic health record (EHR) review. Results: Among a total of 1160 participants (703 [60.6%] female; median [IQR] age, 57.0 [42.1-68.5] years), fewer than half of participants (551 of 1160 [47.5%]) completed a genetics visit. Younger age (odds ratio [OR] for age 18-40 years, 2.98; 95% CI, 1.40-6.53; OR for age 41-65 years, 2.36; 95% CI, 1.22-4.74; OR for age 66-80 years, 2.60; 95% CI, 1.41-4.98 vs age ≥81 years); female sex (OR, 1.49; 95% CI, 1.14-1.96); being married (OR, 1.74; 95% CI, 1.23-2.47) or divorced (OR, 1.80; 95% CI, 1.11-2.91); lower Charlson comorbidity index (OR for score of 0-2, 1.76; 95% CI, 1.16-2.68; OR for score of 3-4, 1.73; 95% CI, 1.18-2.54 vs score of ≥5); EHR patient portal use (OR, 1.42; 95% CI, 1.06-1.89); living closer to a genetics clinic (OR, 1.64; 95% CI, 1.14-2.36 for <8.9 miles vs >20.1 miles); successful results disclosure (OR for disclosure by genetic counselor, 16.32; 95% CI, 8.16-37.45; OR for disclosure by research assistant, 20.30; 95% CI, 10.25-46.31 vs unsuccessful phone disclosure); and having a hereditary cancer result (OR, 2.13; 95% CI, 1.28-3.58 vs other disease risk) were significantly associated with higher rates of genetics visit completion. Preference to follow up with primary care was the most common reported reason for declining a genetics visit (68 of 152 patients [44.7%]). Conclusions and Relevance: This cohort study of a biobank-based population genomic screening program suggests that targeted patient engagement, improving multidisciplinary coordination, and reducing barriers to follow-up care may be necessary for enhancing genetics visit uptake.

Implementation of whole-exome sequencing for pharmacogenomics profiling and exploring its potential clinical utilities.

Whole-exome sequencing (WES) is widely used in clinical settings; however, the exploration of its use in pharmacogenomic analysis remains limited. Our study compared the variant callings for 28 core absorption, distribution, metabolism and elimination genes by WES and array-based technology using clinical trials samples. The results revealed that WES had a positive predictive value of 0.71-0.92 and a sensitivity of single-nucleotide variants between 0.68 and 0.95, compared with array-based technology, for the variants in the commonly targeted regions of the WES and PhamacoScan™ assay. Besides the common variants detected by both assays, WES identified 200-300 exclusive variants per sample, totalling 55 annotated exclusive variants, including important modulators of metabolism such as rs2032582 (ABCB1) and rs72547527 (SULT1A1). This study highlights the potential clinical advantages of using WES to identify a wider range of genetic variations and enabling precision medicine.

Added value of whole-exome and RNA sequencing in advanced and refractory cancer patients with no molecular-based treatment recommendation based on a 90-gene panel.

INTRODUCTION: The objective was to determine the added value of comprehensive molecular profile by whole-exome and RNA sequencing (WES/RNA-Seq) in advanced and refractory cancer patients who had no molecular-based treatment recommendation (MBTR) based on a more limited targeted gene panel (TGP) plus array-based comparative genomic hybridization (aCGH). MATERIALS AND METHODS: In this retrospective analysis, we selected 50 patients previously included in the PROFILER trial (NCT01774409) for which no MBT could be recommended based on a targeted 90-gene panel and aCGH. For each patient, the frozen tumor sample mirroring the FFPE sample used for TGP/aCGH analysis were processed for WES and RNA-Seq. Data from TGP/aCGH were reanalyzed, and together with WES/RNA-Seq, findings were simultaneously discussed at a new molecular tumor board (MTB). RESULTS: After exclusion of variants of unknown significance, a total of 167 somatic molecular alterations were identified in 50 patients (median: 3 [1-10]). Out of these 167 relevant molecular alterations, 51 (31%) were common to both TGP/aCGH and WES/RNA-Seq, 19 (11%) were identified by the TGP/aCGH only and 97 (58%) were identified by WES/RNA-Seq only, including two fusion transcripts in two patients. A MBTR was provided in 4/50 (8%) patients using the information from TGP/aCGH versus 9/50 (18%) patients using WES/RNA-Seq findings. Three patients had similar recommendations based on TGP/aCGH and WES/RNA-Seq. CONCLUSIONS: In advanced and refractory cancer patients in whom no MBTR was recommended from TGP/aCGH, WES/RNA-Seq allowed to identify more alterations which may in turn, in a limited fraction of patients, lead to new MBTR.