npj Genomic Medicine

Assessment of candidate high-grade serous ovarian carcinoma predisposition genes through integrated germline and tumour sequencing
Subramanian DN, Zethoven M, Pishas KI, Marinović ER, McInerny S, Rowley SM, Allan PE, Devereux L, Cheasley D, James PA and Campbell IG
High-grade serous ovarian carcinoma (HGSOC) has a significant hereditary component, only half of which is explained. Previously, we performed germline exome sequencing on BRCA1 and BRCA2-negative HGSOC patients, revealing three proposed and 43 novel candidate genes enriched with rare loss-of-function variants. For validation, we undertook case-control analyses using genomic data from disease-free controls. This confirms enrichment for nearly all previously identified genes. Additionally, one-hundred-and-eleven HGSOC tumours from variant carriers were sequenced alongside other complementary studies, seeking evidence of biallelic inactivation as supportive evidence. PALB2 and ATM validate as HGSOC predisposition genes, with 6/8 germline carrier tumours exhibiting biallelic inactivation accompanied by characteristic mutational signatures. Among candidate genes, only LLGL2 consistently shows biallelic inactivation and protein expression loss, supporting it as a novel HGSOC susceptibility gene. The remaining candidate genes fail to validate. Integrating case-control analyses with tumour sequencing is thus crucial for accurate gene discovery in familial cancer studies.
Whole genome sequencing completes the molecular genetic testing workflow of patients with Lynch syndrome
Horti-Oravecz K, Bozsik A, Pócza T, Vereczkey I, Strausz T, Tóth E, Sedlackova T, Rusnakova D, Szemes T, Likó I, Oláh E, Butz H, Patócs A, Papp J and Grolmusz VK
Multigene panel tests (MGPTs) revolutionized the diagnosis of Lynch syndrome (LS), however noncoding pathogenic variants (PVs) can only be detected by complementary methods including whole genome sequencing (WGS). Here we present a DNA-, RNA- and tumor tissue-based WGS prioritization workflow for patients with a suspicion of LS where MGPT detected no LS-related PV. Among the 100 enrolled patients, MGPT detected 28 simple PVs and an additional 3 complex PVs. Among the 69 MGPT-negative patients, the lack of somatic MLH1 promoter methylation in a patient with a distinguished MLH1 allelic imbalance selected this sample for WGS. This returned a germline deep intronic MLH1 variant, with further functional studies confirming its' pathogenicity. Interestingly, all three complex PVs and the MLH1 deep intronic PV were found to be recurrent at our center. Our straightforward and cost-effective prioritization workflow can optimally include WGS in the genetic diagnosis of LS.
Pre-T cell receptor-α immunodeficiency detected exclusively using whole genome sequencing
Merico D, Sharfe N, Dadi H, Thiruvahindrapuram B, de Rijke J, Dahi Z, Zarrei M, Al Ghamdi A, Al Shaqaq A, Vong L, Scherer SW and Roifman CM
Maturation of αβ lineage T cells in the thymus relies on the formation and cell surface expression of a pre-T cell receptor (TCR) complex, composed of TCRβ chain and pre-TCRα (pTCRα) chain heterodimers, giving rise to a diverse T cell repertoire. Genetic aberrations in key molecules involved in T cell development lead to profound T cell immunodeficiency. Definitive genetic diagnosis guides treatment choices and counseling. In this study, we describe the role of whole genome sequencing (WGS) in providing a definitive diagnosis for a child with T cell deficiency, where targeted panel sequencing of SCID genes and whole exome sequencing had failed. A novel homozygous 8kb deletion in PTCRA, encoding pTCRα, was identified. To date, use of WGS remains restricted and for many geographical regions, is clinically unavailable.
Adaptive evolution of SARS-CoV-2 during a persistent infection for 521 days in an immunocompromised patient
Schmidt H, Schick L, Podlech J, Renzaho A, Lieb B, Diederich S, Hankeln T, Plachter B and Kriege O
Immunocompromised patients struggle to adequately clear viral infections, offering the virus the opportunity to adapt to the immune system in the host. Here we present a case study of a patient undergoing allogeneic hematopoietic stem cell transplantation with a 521-day follow-up of a SARS-CoV-2 infection with the BF.7.21 variant. Virus samples from five time points were submitted to whole genome sequencing. Between the first detection of SARS-CoV-2 infection and its clearance, the patient's virus population acquired 34 amino acid substitutions and 8 deletions in coding regions. With 11 amino acid substitutions in the receptor binding domain of the virus' spike protein, substitutions were 15 times more abundant than expected for a random distribution in this highly functional region. Amongst them were the substitutions S:K417T, S:N440S, S:K444R, S:V445A, S:G446N, S:L452Q, S:N460K, and S:E484V at positions that are notorious for their resistance-mediating effects. The substitution patterns found indicate ongoing adaptive evolution.
Germline structural variant as the cause of Lynch Syndrome in a family from Ecuador
Llargués-Sistac G, Bonjoch L, Muñoz J, Domínguez-Rovira X, Ocaña T, Alvarez-Mora MI, Badenas C, Esteve-Codina A, Reyes-Silva C, Jaramillo-Koupermann G, Rodrigo MT, López-Prades S, Cuatrecasas M, Castells A, Balaguer F, Moreira L, Fernandez G and Castellví-Bel S
Colorectal cancer (CRC) is one of the most common cancers worldwide. Lynch Syndrome (LS) is the most common form of hereditary CRC and it is caused by germline defects in the DNA-mismatch repair (MMR) pathway. It is of extreme importance for affected LS patients and their relatives to identify the germline causative alteration to provide intensified surveillance to those at risk and allow early diagnosis and cancer prevention. Current approaches for LS molecular diagnosis typically involve screening of the MMR genes by targeted gene-panel sequencing and rearrangement screening. We report the identification and characterization of a novel germline structural variant encompassing 48.757 kb, involving the 3'-ends of the MLH1 and LRRFIP2 genes, as the cause of LS in a family of Ecuador. Whole-genome sequencing and transcriptomics allowed the identification of the genomic rearrangement and highlights the importance of the use of these additional approaches to achieve a comprehensive molecular diagnosis in some LS patients.
Exome sequencing of 18,994 ethnically diverse patients with suspected rare Mendelian disorders
Han H, Seo GH, Hyun SI, Kwon K, Ryu SW, Khang R, Lee E, Kim J, Song Y, Jeong WC, Han J, Kim DW, Yang S, Lee S, Jang S, Lee J and Lee H
We investigated the effectiveness of exome sequencing (ES) in diagnosing ethnically diverse patients with rare genetic disorders. A total of 18,994 patients referred to a single reference laboratory for ES between 2020 and 2022 were studied for the diagnostic rate and factors influencing the diagnostic rate. The overall diagnostic rate was 31.8%. Dermatological disorders, skeletal disorders, and neurodevelopmental disorders disease categories, early age-of-onset, presence of consanguinity, and the presence of parental sequencing data were found to be correlated with a higher diagnostic rate. Nearly 68K variants were identified in our dataset at a higher frequency than that observed in gnomAD 4.0. Of these, 507 variants could be classified as likely benign, representing 0.04% of non-benign variants in ClinVar (507/1,433,904) and 0.20% of the non-benign ClinVar variants observed at least once in our cohort (507/276,777). The overall diagnostic rate is comparable to that observed in other large cohort studies with less diverse ethnic backgrounds.
Hospital-wide access to genomic data advanced pediatric rare disease research and clinical outcomes
French CE, Andrews NC, Beggs AH, Boone PM, Brownstein CA, Chopra M, Chou J, Chung WK, D'Gama AM, Doan RN, Ebrahimi-Fakhari D, Goldstein RD, Irons M, Jacobsen C, Kenna M, Lee T, Madden JA, Majmundar AJ, Mann N, Morton SU, Poduri A, Randolph AG, Roberts AE, Roberts S, Sampson MG, Shao DD, Shao W, Sharma A, Shearer E, Shimamura A, Snapper SB, Srivastava S, Thiagarajah JR, Whitman MC, Wojcik MH, Rockowitz S and Sliz P
Boston Children's Hospital has established a genomic sequencing and analysis research initiative to improve clinical care for pediatric rare disease patients. Through the Children's Rare Disease Collaborative (CRDC), the hospital offers CLIA-grade exome and genome sequencing, along with other sequencing types, to patients enrolled in specialized rare disease research studies. The data, consented for broad research use, are harmonized and analyzed with CRDC-supported variant interpretation tools. Since its launch, 66 investigators representing 26 divisions and 45 phenotype-based cohorts have joined the CRDC. These studies enrolled 4653 families, with 35% of analyzed cases having a finding either confirmed or under further investigation. This accessible and harmonized genomics platform also supports additional institutional data collections, research and clinical, and now encompasses 13,800+ patients and their families. This has fostered new research projects and collaborations, increased genetic diagnoses and accelerated innovative research via integration of genomics research with clinical care.
The genetic landscape of autism spectrum disorder in an ancestrally diverse cohort
Gogate A, Kaur K, Khalil R, Bashtawi M, Morris MA, Goodspeed K, Evans P and Chahrour MH
Autism spectrum disorder (ASD) comprises neurodevelopmental disorders with wide variability in genetic causes and phenotypes, making it challenging to pinpoint causal genes. We performed whole exome sequencing on a modest, ancestrally diverse cohort of 195 families, including 754 individuals (222 with ASD), and identified 38,834 novel private variants. In 68 individuals with ASD (~30%), we identified 92 potentially pathogenic variants in 73 known genes, including BCORL1, CDKL5, CHAMP1, KAT6A, MECP2, and SETD1B. Additionally, we identified 158 potentially pathogenic variants in 120 candidate genes, including DLG3, GABRQ, KALRN, KCTD16, and SLC8A3. We also found 34 copy number variants in 31 individuals overlapping known ASD loci. Our work expands the catalog of ASD genetics by identifying hundreds of variants across diverse ancestral backgrounds, highlighting convergence on nervous system development and signal transduction. These findings provide insights into the genetic underpinnings of ASD and inform molecular diagnosis and potential therapeutic targets.
Long-read genome sequencing resolves complex genomic rearrangements in rare genetic syndromes
Showpnil IA, E Hernandez Gonzalez M, Ramadesikan S, Marhabaie M, Daley A, Dublin-Ryan L, Pastore MT, Gurusamy U, Hunter JM, Stone BS, Bartholomew DW, Manickam K, Miller AR, Wilson RK, Stottmann RW and Koboldt DC
Long-read sequencing can often overcome the deficiencies in routine microarray or short-read technologies in detecting complex genomic rearrangements. Here we used Pacific Biosciences circular consensus sequencing to resolve complex rearrangements in two patients with rare genetic anomalies. Copy number variants (CNVs) identified by clinical microarray -chr8p deletion and chr8q duplication in patient 1, and interstitial deletions of chr18q in patient 2-were suggestive of underlying rearrangements. Long-read genome sequencing not only confirmed these CNVs but also revealed their genomic structures. In patient 1, we resolved a novel recombinant chromosome 8 (Rec8)-like rearrangement with a 3.43 Mb chr8q terminal duplication that was linked to a 7.25-8.21 Mb chr8p terminal deletion. In patient 2, we uncovered a novel complex rearrangement involving a 1.17 Mb rearranged segment and four interstitial deletions ranging from 9 bp to 12.39 Mb. Our results underscore the diversity of clinically relevant structural rearrangements and the power of long-read sequencing in unraveling their nuanced architectures.
Analysis of exonic deletions in a large population study provides novel insights into NRXN1 pathology
Montalbano S, Krebs MD, Rosengren A, Vaez M, Hellberg KG, Mortensen PB, Børglum AD, Geschwind DH, , Raznahan A, Thompson WK, Helenius D, Werge T and Ingason A
The NRXN1 locus is a hotspot for non-recurrent copy number variants and exon-disrupting NRXN1 deletions have been associated with increased risk of neurodevelopmental disorders in case-control studies. However, corresponding population-based estimates of prevalence and disease-associated risk are currently lacking. Also, most studies have not differentiated between deletions affecting exons of different NRXN1 splice variants nor considered intronic deletions. We used the iPSYCH2015 case-cohort sample to obtain unbiased estimates of the prevalence of NRXN1 deletions and their associated risk of autism, schizophrenia, depression, and ADHD. Most exon-disrupting deletions affected exons specific to the alpha isoform, and almost half of the non-exonic deletions represented a previously reported segregating founder deletion. Carriage of exon-disrupting NRXN1 deletions was associated with a threefold and twofold increased risk of autism and ADHD, respectively, whereas no significantly increased risk of depression or schizophrenia was observed. Our results highlight the importance of using population-based samples in genetic association studies.
Efficient reinterpretation of rare disease cases using Exomiser
Vestito L, Jacobsen JOB, Walker S, Cipriani V, Harris NL, Haendel MA, Mungall CJ, Robinson P and Smedley D
Whole genome sequencing has transformed rare disease research; however, 50-80% of rare disease patients remain undiagnosed after such testing. Regular reanalysis can identify new diagnoses, especially in newly discovered disease-gene associations, but efficient tools are required to support clinical interpretation. Exomiser, a phenotype-driven variant prioritisation tool, fulfils this role; within the 100,000 Genomes Project (100kGP), diagnoses were identified after reanalysis in 463 (2%) of 24,015 unsolved patients after previous analysis for variants in known disease genes. However, extensive manual interpretation was required. This led us to develop a reanalysis strategy to efficiently reveal candidates from recent disease gene discoveries or newly designated pathogenic/likely pathogenic variants. Optimal settings to highlight new candidates from Exomiser reanalysis were identified with high recall (82%) and precision (88%) when including Exomiser's automated ACMG/AMP classifier, which correctly converted 92% of variants from unknown significance to pathogenic/likely pathogenic. In conclusion, Exomiser efficiently reinterprets previously unsolved cases.
Somatic mutation in autosomal dominant polycystic kidney disease revealed by deep sequencing human kidney cysts
Mallawaarachchi AC, Hort Y, Wedd L, Lo K, Senum S, Toumari M, Chen W, Utsiwegota M, Mawson J, Leslie S, Laurence J, Anderson L, Snelling P, Salomon R, Rangan GK, Furlong T, Shine J and Cowley MJ
Autosomal Dominant Polycystic Kidney Disease (ADPKD) results in progressive cysts that lead to kidney failure, and is caused by heterozygous germline variants in PKD1 or PKD2. Cyst pathogenesis is not definitively understood. Somatic second-hit mutations have been implicated in cyst pathogenesis, though technical sequencing challenges have limited investigation. We used unique molecular identifiers, high-depth massively parallel sequencing and custom analysis techniques to identify somatic second-hit mutations in 24 whole cysts from disparate regions of six human ADPKD kidneys, utilising replicate samples and orthogonal confirmation. Average depth of coverage of 1166 error-corrected reads for PKD1 and 539 reads for PKD2 was obtained. 58% (14/24) of cysts had a detectable PKD1 somatic variant, with 5/6 participants having at least one cyst with a somatic variant. We demonstrate that low-frequency somatic mutations are detectable in a proportion of cysts from end-stage ADPKD human kidneys. Further studies are required to understand the drivers of this somatic mutation.
Functional assessment of IDUA variants of uncertain significance identified by newborn screening
Yu SH, Kubaski F, Arno G, Phinney W, Wood TC, Flanagan-Steet H, Pollard LM and Steet R
With the expansion of newborn screening efforts for MPS disorders, the number of identified variants of uncertain significance in IDUA continues to increase. To better define functional consequences of identified IDUA variants, we developed a HEK293-based expression platform that can be used to determine the relative specific activity of variant α-iduronidases by combining a fluorescence-based activity assay and semi-quantitative western blotting. We employed the current platform to characterize over thirty different IDUA variants, including known benign and pathogenic variants, as well as multiple variants of uncertain significance identified through newborn screening. This analysis allowed the stratification of variant enzymes based on their relative specific activity, and uncovered distinct effects of the different variants on enzyme folding, processing, and stability. While relative specific activity serves as a useful first-level test for enzyme function, our observations reinforce the need for secondary analyses of enzyme function to fully assess variant pathogenicity.
Clinical and genetic characterization of patients with late onset Wilson's disease
Yang W, Yang Y, Wang H, Wang J and Zhang S
Wilson's disease (WD) typically manifests in children and young adults, with little knowledge of its late-onset forms. In this study, we performed a retrospective cohort study of 105 WD patients (99 index cases, 6 siblings) with an onset age ≥35 years. We compared 99 index late-onset patients with 1237 early-onset patients and analyzed the ATP7B variant penetrance referring to the Genome Aggregation Database (gnomAD). Sixty-two ATP7B variants were identified in the late-onset patients, among which A874V, V1106I, R919G, and T935M were correlated with late presentation of WD. Regarding gnomAD, V1106I and T935M exhibited significantly low penetrance, and there is a lack of patients carrying a genotype of V1106I/V1106I, R919G/R919G, T935M/T935M, V1106I/T935M, V1106I/R919G, or T935M/R919G. Our data revealed that patients carrying a combination of two late-onset variants may be overlooked due to atypical or lack of WD symptoms, which may provide valuable insights into the genetic basis and diagnosis of WD.
Author Correction: The genetic landscape of autism spectrum disorder in an ancestrally diverse cohort
Gogate A, Kaur K, Khalil R, Bashtawi M, Morris MA, Goodspeed K, Evans P and Chahrour MH
Comparative transcriptomic, epigenomic and immunological analyses identify drivers of disparity in high-grade serous ovarian cancer
Huang H, Keathley R, Kim U, Cardenas H, Xie P, Wei J, Lengyel E, Nephew KP, Zhao G, Fu Z, Barber EL, Kocherginsky M, Bae-Jump V, Zhang B and Matei D
Black women face the highest mortality-to-incidence ratio from high grade serous ovarian cancer (HGSOC). This study investigated biological differences in HGSOC tumors from Black vs. White women. HGSOC from 35 Black and 31 White patients were analyzed by Infinium Methyation-EPIC array and RNA sequencing. 191 CpG sites were differentially methylated (FDR < 0.05, β value change> 10%) and 277 genes were differentially expressed (FDR < 0.05). Gene Ontology identified enriched pathways related to DNA damage response, p53/apoptosis signaling, and cholesterol/lipid metabolism directly connected with genes like INSR, FOXA1 and FOXB1. INSR and FOXA1 knockdown enhanced cisplatin sensitivity and inhibited cell proliferation and colony formation. Tumors from Black patients were infiltrated by fewer CD4+ naïve and regulatory T-cells. Overall, differences in DNA methylation, transcriptomic profiles and immune cell infiltration were detected in tumors from Black vs. White patients. Further investigation is warranted into how these differences may affect treatment response and outcomes in Black women.
Implementing genomic newborn screening as an effective public health intervention: sidestepping the hype and criticism
Friedman JM
Genome-wide sequencing of the DNA that can be obtained from a newborn screening blood spot could provide predictions of thousands of genetic diseases that are not currently included in universal newborn screening. Most of the serious ethical, legal, privacy, and social concerns raised by genome sequencing of all infants can be avoided by implementing genomic newborn screening in accordance with widely-accepted public health criteria.
Coding and non-coding variants in the ciliopathy gene CFAP410 cause early-onset non-syndromic retinal degeneration
Sangermano R, Gupta P, Price C, Han J, Navarro J, Condroyer C, Place EM, Antonio A, Mukai S, Zanlonghi X, Sahel JA, DiTroia S, O'Heir E, Duncan JL, Pierce EA, Zeitz C, Audo I, Huckfeldt RM and Bujakowska KM
Inherited retinal degenerations are blinding genetic disorders characterized by high genetic and phenotypic heterogeneity. In this retrospective study, we describe sixteen families with early-onset non-syndromic retinal degenerations in which affected probands carried rare bi-allelic variants in CFAP410, a ciliary gene previously associated with recessive Jeune syndrome. We detected twelve variants, eight of which were novel, including c.373+91A>G, which led to aberrant splicing. To our knowledge this is the first likely pathogenic deep-intronic variant identified in this gene. Analysis of all reported and novel CFAP410 variants revealed no clear correlation between the severity of the CFAP410-associated phenotypes and the identified causal variants. This is supported by the fact that the frequently encountered missense variant p.(Arg73Pro), often found in syndromic cases, was also associated with non-syndromic retinal degeneration. This study expands the current knowledge of CFAP410-associated ciliopathy by enriching its mutational landscape and supports its association with non-syndromic retinal degeneration.
Germline sequence variation in cancer genes in Rwandan breast and prostate cancer cases
Manirakiza AV, Baichoo S, Uwineza A, Dukundane D, Uwinkindi F, Ngendahayo E, Rubagumya F, Muhawenimana E, Nsabimana N, Nzeyimana I, Maniragaba T, Ntirenganya F, Rurangwa E, Mugenzi P, Mutamuliza J, Runanira D, Niyibizi BA, Rugengamanzi E, Besada J, Nielsen SM, Bucknor B, Nussbaum RL, Koeller D, Andrews C, Mutesa L, Fadelu T and Rebbeck TR
Cancer genetic data from Sub-Saharan African (SSA) are limited. Patients with female breast (fBC), male breast (mBC), and prostate cancer (PC) in Rwanda underwent germline genetic testing and counseling. Demographic and disease-specific information was collected. A multi-cancer gene panel was used to identify germline Pathogenic Variants (PV) and Variants of Uncertain Significance (VUS). 400 patients (201 with BC and 199 with PC) were consented and recruited to the study. Data was available for 342 patients: 180 with BC (175 women and 5 men) and 162 men with PC. PV were observed in 18.3% fBC, 4.3% PC, and 20% mBC. BRCA2 was the most common PV. Among non-PV carriers, 65% had ≥1 VUS: 31.8% in PC and 33.6% in BC (female and male). Our findings highlight the need for germline genetic testing and counseling in cancer management in SSA.
Common protein-altering variant in GFAP is associated with white matter lesions in the older Japanese population
Furuta Y, Akiyama M, Hirabayashi N, Honda T, Shibata M, Ohara T, Hata J, Terao C, Momozawa Y, Tatewaki Y, Taki Y, Nakaji S, Maeda T, Ono K, Mimura M, Nakashima K, Iga JI, Takebayashi M, Ninomiya T and
The genetic architecture of white matter lesions (WMLs) in Asian populations has not been well-characterized. Here, we performed a genome-wide association study (GWAS) to identify loci associated with the WML volume. Brain MRI and DNA samples were collected from 9479 participants in the Japan Prospective Studies Collaboration for Aging and Dementia (JPSC-AD). The GWAS confirmed three known WML-associated loci (SH3PXD2A, GFAP, and TRIM47). The lead variant of GFAP was a common missense variant (p.D295N) in East Asians. Meta-GWAS using the publicly available summary statistics of UK Biobank identified one previously unreported locus 6q23.2 (SLC2A12). Integration with expression quantitative trait locus data implied the newly identified locus affects SLC2A12 expression. The effect sizes of 20 lead variants at the WML-associated loci were moderately correlated between JPSC-AD and UK Biobank. These results indicate that the alteration in GFAP protein caused by the common missense variant in East Asians influences the WML volume.
Implementing genomic medicine in clinical practice for adults with undiagnosed rare diseases
Ahn JH, Yoon JG, Cho J, Lee S, Kim S, Kim MJ, Kim SY, Lee ST, Chu K, Lee SK, Kim HJ, Youn J, Jang JH, Chae JH, Moon J and Cho JW
The global burden of undiagnosed diseases, particularly in adults, is rising due to their significant socioeconomic impact. To address this, we enrolled 232 adult probands with undiagnosed conditions, utilizing bioinformatics tools for genetic analysis. Alongside exome and genome sequencing, repeat-primed PCR and Cas9-mediated nanopore sequencing were applied to suspected short tandem repeat disorders. Probands were classified into probable genetic (n = 128) or uncertain (n = 104) origins. The study found genetic causes in 66 individuals (28.4%) and non-genetic causes in 12 (5.2%), with a longer diagnostic journey for those in the probable genetic group or with pediatric symptom onset, emphasizing the need for increased efforts in these populations. Genetic diagnoses facilitated effective surveillance, cascade screening, drug repurposing, and pregnancy planning. This study demonstrates that integrating sequencing technologies improves diagnostic accuracy, may shorten the time to diagnosis, and enhances personalized management for adults with undiagnosed diseases.