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.

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.

Detection of elusive DNA copy-number variations in hereditary disease and cancer through the use of noncoding and off-target sequencing reads

Copy-number variants (CNVs) play a substantial role in the molecular pathogenesis of hereditary disease and cancer, as well as in normal human interindividual variation. However, they are still rather difficult to identify in mainstream sequencing projects, especially involving exome sequencing, because they often occur in DNA regions that are not targeted for analysis. To overcome this problem, we developed OFF-PEAK, a user-friendly CNV detection tool that builds on a denoising approach and the use of "off-target" DNA reads, which are usually discarded by sequencing pipelines. We benchmarked OFF-PEAK on data from targeted sequencing of 96 cancer samples, as well as 130 exomes of individuals with inherited retinal disease from three different populations. For both sets of data, OFF-PEAK demonstrated excellent performance (>95% sensitivity and >80% specificity vs. experimental validation) in detecting CNVs from in silico data alone, indicating its immediate applicability to molecular diagnosis and genetic research.

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.

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.

An immunogenomic exome landscape of triple positive primary antiphospholipid patients

Primary antiphospholipid syndrome is characterized by thrombosis and autoantibodies directed against phospholipids or associated proteins. The genetic etiology of PAPS remains unknown. We enrolled 21 patients with thromboembolic events associated to lupus anticoagulant, anticardiolipin and anti β2 glycoprotein1 autoantibodies. We performed whole exome sequencing and a systematic variant-based analysis in genes associated with thrombosis, in candidate genes previously associated with APS or inborn errors of immunity. Data were compared to public databases and to a control cohort of 873 non-autoimmune patients. Variants were identified following a state-of-the-art pipeline. Enrichment analysis was performed by comparing with the control cohort. We found an absence of significant HLA bias and genetic heterogeneity in these patients, including when testing combinations of rare variants in genes encoding for proteins involved in thrombosis and of variants in genes linked with inborn errors of immunity. These results provide evidence of genetic heterogeneity in PAPS, even in a homogenous series of triple positive patients. At the individual scale, a combination of variants may participate to the breakdown of B cell tolerance and to the vessel damage.

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.

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.

Diversity of CFTR variants across ancestries characterized using 454,727 UK biobank whole exome sequences

BACKGROUND

Limited understanding of the diversity of variants in the cystic fibrosis transmembrane conductance regulator (CFTR) gene across ancestries hampers efforts to advance molecular diagnosis of cystic fibrosis (CF). The consequences pose a risk of delayed diagnoses and subsequently worsened health outcomes for patients. Therefore, characterizing the spectrum of CFTR variants across ancestries is critical for revolutionizing molecular diagnoses of CF.

METHODS

We analyzed 454,727 UK Biobank (UKBB) whole-exome sequences to characterize the diversity of CFTR variants across ancestries. Using the PanUKBB classification, the participants were assigned into six major groups: African (AFR), American/American Admixed (AMR), Central South Asia (CSA), East Asian (EAS), European (EUR), and Middle East (MID). We segregated ancestry-specific CFTR variants, including those that are CF-causing or clinically relevant. The ages of certain CF-causing variants were determined and analyzed for selective pressure effects, and curated phenotype analysis was performed for participants with clinically relevant CFTR genotypes.

RESULTS

We detected over 4000 CFTR variants, including novel ancestry-specific variants, across six ancestries. Europeans had the most unique CFTR variants [n = 2212], while the American group had the least unique variants [n = 23]. F508del was the most prevalent CF-causing variant found in all ancestries, except in EAS, where V520F was the most prevalent. Common EAS variants such as 3600G > A, V456A, and V520, which appeared approximately 270, 215, and 338 generations ago, respectively, did not show evidence of selective pressure. Sixteen participants had two CF-causing variants, with two being diagnosed with CF. We found 154 participants harboring a CF-causing and varying clinical consequences (VCC) variant. Phenotype analysis performed for participants with multiple clinically relevant variants returned significant associations with CF and its pulmonary phenotypes [Bonferroni-adjusted p < 0.05].

CONCLUSIONS

We leveraged the UKBB database to comprehensively characterize the broad spectrum of CFTR variants across ancestries. The detection of over 4000 CFTR variants, including several ancestry-specific and uncharacterized CFTR variants, warrants the need for further characterization of their functional and clinical relevance. Overall, the presentation of classical CF phenotypes seen in non-CF diagnosed participants with more than one CF-causing variant indicates that they may benefit from current CFTR modulator therapies.

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.

Validated WGS and WES protocols proved saliva-derived gDNA as an equivalent to blood-derived gDNA for clinical and population genomic analyses

BACKGROUND

Whole exome sequencing (WES) and whole genome sequencing (WGS) have become standard methods in human clinical diagnostics as well as in population genomics (POPGEN). Blood-derived genomic DNA (gDNA) is routinely used in the clinical environment. Conversely, many POPGEN studies and commercial tests benefit from easy saliva sampling. Here, we evaluated the quality of variant call sets and the level of genotype concordance of single nucleotide variants (SNVs) and small insertions and deletions (indels) for WES and WGS using paired blood- and saliva-derived gDNA isolates employing genomic reference-based validated protocols.

METHODS

The genomic reference standard Coriell NA12878 was repeatedly analyzed using optimized WES and WGS protocols, and data calls were compared with the truth dataset published by the Genome in a Bottle Consortium. gDNA was extracted from the paired blood and saliva samples of 10 participants and processed using the same protocols. A comparison of paired blood-saliva call sets was performed in the context of WGS and WES genomic reference-based technical validation results.

RESULTS

The quality pattern of called variants obtained from genomic-reference-based technical replicates correlates with data calls of paired blood-saliva-derived samples in all levels of tested examinations despite a higher rate of non-human contamination found in the saliva samples. The F1 score of 10 blood-to-saliva-derived comparisons ranged between 0.8030-0.9998 for SNVs and between 0.8883-0.9991 for small-indels in the case of the WGS protocol, and between 0.8643-0.999 for SNVs and between 0.7781-1.000 for small-indels in the case of the WES protocol.

CONCLUSION

Saliva may be considered an equivalent material to blood for genetic analysis for both WGS and WES under strict protocol conditions. The accuracy of sequencing metrics and variant-detection accuracy is not affected by choosing saliva as the gDNA source instead of blood but much more significantly by the genomic context, variant types, and the sequencing technology used.

Genetic Factors and Long-term Treatment-Related Neurocognitive Deficits, Anxiety, and Depression in Childhood Leukemia Survivors: An Exome-Wide Association Study

BACKGROUND

An increased risk of neurocognitive deficits, anxiety, and depression has been reported in childhood cancer survivors.

METHODS

We analyzed associations of neurocognitive deficits, as well as anxiety and depression, with common and rare genetic variants derived from whole-exome sequencing data of acute lymphoblastic leukemia (ALL) survivors from the PETALE cohort. In addition, significant associations were assessed using stratified and multivariable analyses. Next, top-ranking common associations were analyzed in an independent SJLIFE replication cohort of ALL survivors.

RESULTS

Significant associations were identified in the entire discovery cohort (N = 229) between the AK8 gene and changes in neurocognitive function, whereas PTPRZ1, MUC16, TNRC6C-AS1 were associated with anxiety. Following stratification according to sex, the ZNF382 gene was linked to a neurocognitive deficit in males, whereas APOL2 and C6orf165 were associated with anxiety and EXO5 with depression. Following stratification according to prognostic risk groups, the modulatory effect of rare variants on depression was additionally found in the CYP2W1 and PCMTD1 genes. In the replication SJLIFE cohort (N = 688), the male-specific association in the ZNF382 gene was not significant; however, a P value<0.05 was observed when the entire SJLIFE cohort was analyzed. ZNF382 was significant in males in the combined cohorts as shown by meta-analyses as well as the depression-associated gene EXO5.

CONCLUSIONS

Further research is needed to confirm whether the current findings, along with other known risk factors, may be valuable in identifying patients at increased risk of these long-term complications.

IMPACT

Our results suggest that specific genes may be related to increased neuropsychological consequences.

Unraveling the genetic architecture of congenital vertebral malformation with reference to the developing spine

Congenital vertebral malformation, affecting 0.13-0.50 per 1000 live births, has an immense locus heterogeneity and complex genetic architecture. In this study, we analyze exome/genome sequencing data from 873 probands with congenital vertebral malformation and 3794 control individuals. Clinical interpretation identifies Mendelian etiologies in 12.0% of the probands and reveals a muscle-related disease mechanism. Gene-based burden test of ultra-rare variants identifies risk genes with large effect sizes (ITPR2, TBX6, TPO, H6PD, and SEC24B). To further investigate the biological relevance of the genetic association signals, we perform single-nucleus RNAseq on human embryonic spines. The burden test signals are enriched in the notochord at early developmental stages and myoblast/myocytes at late stages, highlighting their critical roles in the developing spine. Our work provides insights into the developmental biology of the human spine and the pathogenesis of spine malformation.

Development of a coding SNP panel for tracking the origin of whole-exome sequencing samples

Whole-exome sequencing (WES) is widely used to diagnose complex genetic diseases and rare conditions. The implementation of a robust and effective quality control system for sample identification and tracking throughout the WES process is essential. We established a multiplex panel that included 22 coding single-nucleotide polymorphism (cSNP) loci. The personal identification and paternity identification abilities of the panel were evaluated, and a preliminary validation of the practical feasibility of the panel was conducted in a clinical WES case. These results indicate that the cSNP panel could be a useful tool for sample tracking in WES.

In-Depth Genomic Analysis: The New Challenge in Congenital Heart Disease

The use of next-generation sequencing has provided new insights into the causes and mechanisms of congenital heart disease (CHD). Examinations of the whole exome sequence have detected detrimental gene variations modifying single or contiguous nucleotides, which are characterised as pathogenic based on statistical assessments of families and correlations with congenital heart disease, elevated expression during heart development, and reductions in harmful protein-coding mutations in the general population. Patients with CHD and extracardiac abnormalities are enriched for gene classes meeting these criteria, supporting a common set of pathways in the organogenesis of CHDs. Single-cell transcriptomics data have revealed the expression of genes associated with CHD in specific cell types, and emerging evidence suggests that genetic mutations disrupt multicellular genes essential for cardiogenesis. Metrics and units are being tracked in whole-genome sequencing studies.

Prenatal diagnosis of polycystic kidney caused by biallelic hypomorphic variants in the PKD1 gene

Heterozygous loss-of-function variants in the PKD1 gene are commonly associated with adult-onset autosomal dominant polycystic kidney disease (ADPKD), where the formation of renal cysts depends on the dosage of the PKD1 gene. Biallelic null PKD1 variants are not viable, but biallelic hypomorphic variants could lead to early-onset PKD. We report a non-consanguineous Chinese family with recurrent fetal polycystic kidney and negative findings in the coding region of the PKHD1 gene or chromosomal microarray analysis. Trio exome analysis revealed compound heterozygous variants of uncertain significance in the PKD1 gene in the index pregnancy: a novel paternally inherited c.7863 + 5G > C and a maternally inherited c.9739C > T, p.(Arg3247Cys). Segregation analysis through long-range PCR followed by nested PCR and Sanger sequencing confirmed another affected fetus had both variants, while the other two normal siblings and the parents carried either variant. Thus, these two variants, both of which were hypomorphic as opposed to null variants, co-segregated with prenatal onset polycystic kidney disease in this family. Functional studies are needed to further determine the impact of these two variants. Our findings highlight the biallelic inheritance of hypomorphic PKD1 variants causing prenatal onset polycystic kidney disease, which provides a better understanding of phenotype-genotype correlation and valuable information for reproductive counseling.

Should secondary pharmacogenomic variants be actively screened and reported when diagnostic genome-wide sequencing is performed in a child?

This white paper was prepared by the Global Alliance for Genomics and Health Regulatory and Ethics Work Stream's Pediatric Task Team to review and provide perspective with respect to ethical, legal, and social issues regarding the return of secondary pharmacogenomic variants in children who have a serious disease or developmental disorder and are undergoing exome or genome sequencing to identify a genetic cause of their condition. We discuss actively searching for and reporting pharmacogenetic/genomic variants in pediatric patients, different methods of returning secondary pharmacogenomic findings to the patient/parents and/or treating clinicians, maintaining these data in the patient's health record over time, decision supports to assist using pharmacogenetic results in future treatment decisions, and sharing information in public databases to improve the clinical interpretation of pharmacogenetic variants identified in other children. We conclude by presenting a series of points to consider for clinicians and policymakers regarding whether, and under what circumstances, routine screening and return of pharmacogenomic variants unrelated to the indications for testing is appropriate in children who are undergoing genome-wide sequencing to assist in the diagnosis of a suspected genetic disease.

Mobile element insertions in rare diseases: a comparative benchmark and reanalysis of 60,000 exome samples

Mobile element insertions (MEIs) are a known cause of genetic disease but have been underexplored due to technical limitations of genetic testing methods. Various bioinformatic tools have been developed to identify MEIs in Next Generation Sequencing data. However, most tools have been developed specifically for genome sequencing (GS) data rather than exome sequencing (ES) data, which remains more widely used for routine diagnostic testing. In this study, we benchmarked six MEI detection tools (ERVcaller, MELT, Mobster, SCRAMble, TEMP2 and xTea) on ES data and on GS data from publicly available genomic samples (HG002, NA12878). For all the tools we evaluated sensitivity and precision of different filtering strategies. Results show that there were substantial differences in tool performance between ES and GS data. MELT performed best with ES data and its combination with SCRAMble increased substantially the detection rate of MEIs. By applying both tools to 10,890 ES samples from Solve-RD and 52,624 samples from Radboudumc we were able to diagnose 10 patients who had remained undiagnosed by conventional ES analysis until now. Our study shows that MELT and SCRAMble can be used reliably to identify clinically relevant MEIs in ES data. This may lead to an additional diagnosis for 1 in 3000 to 4000 patients in routine clinical ES.

Evaluation of the diagnostic accuracy of exome sequencing and its impact on diagnostic thinking for patients with rare disease in a publicly funded health care system: A prospective cohort study

PURPOSE

To evaluate the diagnostic utility of publicly funded clinical exome sequencing (ES) for patients with suspected rare genetic diseases.

METHODS

We prospectively enrolled 297 probands who met eligibility criteria and received ES across 5 sites in Ontario, Canada, and extracted data from medical records and clinician surveys. Using the Fryback and Thornbury Efficacy Framework, we assessed diagnostic accuracy by examining laboratory interpretation of results and assessed diagnostic thinking by examining the clinical interpretation of results and whether clinical-molecular diagnoses would have been achieved via alternative hypothetical molecular tests.

RESULTS

Laboratories reported 105 molecular diagnoses and 165 uncertain results in known and novel genes. Of these, clinicians interpreted 102 of 105 (97%) molecular diagnoses and 6 of 165 (4%) uncertain results as clinical-molecular diagnoses. The 108 clinical-molecular diagnoses were in 104 families (35% diagnostic yield). Each eligibility criteria resulted in diagnostic yields of 30% to 40%, and higher yields were achieved when >2 eligibility criteria were met (up to 45%). Hypothetical tests would have identified 61% of clinical-molecular diagnoses.

CONCLUSION

We demonstrate robustness in eligibility criteria and high clinical validity of laboratory results from ES testing. The importance of ES was highlighted by the potential 40% of patients that would have gone undiagnosed without this test.

Genetic architecture and biology of youth-onset type 2 diabetes

The prevalence of youth-onset type 2 diabetes (T2D) and childhood obesity has been rising steadily1, producing a growing public health concern1 that disproportionately affects minority groups2. The genetic basis of youth-onset T2D and its relationship to other forms of diabetes are unclear3. Here we report a detailed genetic characterization of youth-onset T2D by analysing exome sequences and common variant associations for 3,005 individuals with youth-onset T2D and 9,777 adult control participants matched for ancestry, including both males and females. We identify monogenic diabetes variants in 2.4% of individuals and three exome-wide significant (P < 2.6 × 10-6) gene-level associations (HNF1A, MC4R, ATXN2L). Furthermore, we report rare variant association enrichments within 25 gene sets related to obesity, monogenic diabetes and β-cell function. Many youth-onset T2D associations are shared with adult-onset T2D, but genetic risk factors of all frequencies-and rare variants in particular-are enriched within youth-onset T2D cases (5.0-fold increase in the rare variant and 3.4-fold increase in common variant genetic liability relative to adult-onset cases). The clinical presentation of participants with youth-onset T2D is influenced in part by the frequency of genetic risk factors within each individual. These findings portray youth-onset T2D as a heterogeneous disease situated on a spectrum between monogenic diabetes and adult-onset T2D.

An integrative framework for clinical diagnosis and knowledge discovery from exome sequencing data

Non-silent single nucleotide genetic variants, like nonsense changes and insertion-deletion variants, that affect protein function and length substantially are prevalent and are frequently misclassified. The low sensitivity and specificity of existing variant effect predictors for nonsense and indel variations restrict their use in clinical applications. We propose the Pathogenic Mutation Prediction (PMPred) method to predict the pathogenicity of single nucleotide variations, which impair protein function by prematurely terminating a protein's elongation during its synthesis. The prediction starts by monitoring functional effects (Gene Ontology annotation changes) of the change in sequence, using an existing ensemble machine learning model (UniGOPred). This, in turn, reveals the mutations that significantly deviate functionally from the wild-type sequence. We have identified novel harmful mutations in patient data and present them as motivating case studies. We also show that our method has increased sensitivity and specificity compared to state-of-the-art, especially in single nucleotide variations that produce large functional changes in the final protein. As further validation, we have done a comparative docking study on such a variation that is misclassified by existing methods and, using the altered binding affinities, show how PMPred can correctly predict the pathogenicity when other tools miss it. PMPred is freely accessible as a web service at https://pmpred.kansil.org/, and the related code is available at https://github.com/kansil/PMPred.

DNA methylome, R-loop and clinical exome profiling of patients with sporadic amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disorder characterized by the death of motor neurons, the aetiology of which is essentially unknown. Here, we present an integrative epigenomic study in blood samples from seven clinically characterised sporadic ALS patients to elucidate molecular factors associated with the disease. We used clinical exome sequencing (CES) to study DNA variants, DNA-RNA hybrid immunoprecipitation sequencing (DRIP-seq) to assess R-loop distribution, and reduced representation bisulfite sequencing (RRBS) to examine DNA methylation changes. The above datasets were combined to create a comprehensive repository of genetic and epigenetic changes associated with the ALS cases studied. This repository is well-suited to unveil new correlations within individual patients and across the entire patient cohort. The molecular attributes described here are expected to guide further mechanistic studies on ALS, shedding light on the underlying genetic causes and facilitating the development of new epigenetic therapies to combat this life-threatening disease.

Human whole-exome genotype data for Alzheimer's disease

The heterogeneity of the whole-exome sequencing (WES) data generation methods present a challenge to a joint analysis. Here we present a bioinformatics strategy for joint-calling 20,504 WES samples collected across nine studies and sequenced using ten capture kits in fourteen sequencing centers in the Alzheimer's Disease Sequencing Project. The joint-genotype called variant-called format (VCF) file contains only positions within the union of capture kits. The VCF was then processed specifically to account for the batch effects arising from the use of different capture kits from different studies. We identified 8.2 million autosomal variants. 96.82% of the variants are high-quality, and are located in 28,579 Ensembl transcripts. 41% of the variants are intronic and 1.8% of the variants are with CADD > 30, indicating they are of high predicted pathogenicity. Here we show our new strategy can generate high-quality data from processing these diversely generated WES samples. The improved ability to combine data sequenced in different batches benefits the whole genomics research community.

TG468: a text graph convolutional network for predicting clinical response to immune checkpoint inhibitor therapy

Enhancing cancer treatment efficacy remains a significant challenge in human health. Immunotherapy has witnessed considerable success in recent years as a treatment for tumors. However, due to the heterogeneity of diseases, only a fraction of patients exhibit a positive response to immune checkpoint inhibitor (ICI) therapy. Various single-gene-based biomarkers and tumor mutational burden (TMB) have been proposed for predicting clinical responses to ICI; however, their predictive ability is limited. We propose the utilization of the Text Graph Convolutional Network (GCN) method to comprehensively assess the impact of multiple genes, aiming to improve the predictive capability for ICI response. We developed TG468, a Text GCN model framing drug response prediction as a text classification task. By combining natural language processing (NLP) and graph neural network techniques, TG468 effectively handles sparse and high-dimensional exome sequencing data. As a result, TG468 can distinguish survival time for patients who received ICI therapy and outperforms single gene biomarkers, TMB and some classical machine learning models. Additionally, TG468's prediction results facilitate the identification of immune status differences among specific patient types in the Cancer Genome Atlas dataset, providing a rationale for the model's predictions. Our approach represents a pioneering use of a GCN model to analyze exome data in patients undergoing ICI therapy and offers inspiration for future research using NLP technology to analyze exome sequencing data.

Diagnostic genome sequencing improves diagnostic yield: a prospective single-centre study in 1000 patients with inherited eye diseases

PURPOSE

Genome sequencing (GS) is expected to reduce the diagnostic gap in rare disease genetics. We aimed to evaluate a scalable framework for genome-based analyses 'beyond the exome' in regular care of patients with inherited retinal degeneration (IRD) or inherited optic neuropathy (ION).

METHODS

PCR-free short-read GS was performed on 1000 consecutive probands with IRD/ION in routine diagnostics. Complementary whole-blood RNA-sequencing (RNA-seq) was done in a subset of 74 patients. An open-source bioinformatics analysis pipeline was optimised for structural variant (SV) calling and combined RNA/DNA variation interpretation.

RESULTS

A definite genetic diagnosis was established in 57.4% of cases. For another 16.7%, variants of uncertain significance were identified in known IRD/ION genes, while the underlying genetic cause remained unresolved in 25.9%. SVs or alterations in non-coding genomic regions made up for 12.7% of the observed variants. The RNA-seq studies supported the classification of two unclear variants.

CONCLUSION

GS is feasible in clinical practice and reliably identifies causal variants in a substantial proportion of individuals. GS extends the diagnostic yield to rare non-coding variants and enables precise determination of SVs. The added diagnostic value of RNA-seq is limited by low expression levels of the major IRD disease genes in blood.

Genomics and tumor microenvironment of breast mucoepidermoid carcinoma based on whole-exome and RNA sequencing

Mammary mucoepidermoid carcinoma (MEC) is a rare entity. The molecular characteristics of breast MEC have not been fully investigated due to its rarity. We performed a retrospective study among 1000 patients with breast carcinomas and identified four cases of breast MEC. Clinical and demographic data were collected. Immunohistochemistry panels which were used to diagnose salivary gland MEC and breast carcinomas were also performed. MAML2 rearrangements were detected by FISH and fusion partners were identified by RNA sequencing. Whole-exome sequencing (WES) was used to reveal the genomes of these four breast MEC. Then, the biological functions and features of breast MEC were further compared with those of invasive breast carcinomas and salivary gland MEC.According to Ellis and Auclair's methods, these four breast MEC could be classified as low-grade breast MEC. All the patients were alive, and disease-free survival (PFS) ranged from 20 months to 67 months. Among these four breast MEC, two cases were triple-negative, and the other two cases were found to be ER positive, with one also showing HER2 equivocal by immunohistochemical staining, but no amplification in FISH. FISH analysis confirmed the presence of the MAML2 translocation in three of four tumors, and CRTC1-MAML2 fusion was confirmed in two of them by RNA-sequencing. The average coverage size of WES for the tumor mutation burden estimation was 32 Mb. MUC4, RP1L1 and QRICH2 mutations were identified in at least three tumors, and these mutation also existed in breast invasive carcinoma databases (TCGA, Cell 2015; TCGA, Nature 2012). The results showed that there were many genes in breast MEC overlapping with the breast invasive carcinoma databases mentioned above, range from 5 to 63 genes (median:21 genes). Next, we assessed immune cell infiltration levels in these tumors. In all these tumors, M2 macrophages and plasma cell were in the high infiltration group. Our breast MEC showed different results from the salivary gland MEC, whose plasma cells were in the low infiltration group. Overall, we first analyzed the genomics and tumor microenvironment of breast mucoepidermoid carcinoma and proposed our hypothesis that although MECs arising in the breast resemble their salivary gland counterparts phenotypically, our findings indicate that breast MECs probably resemble invasive breast carcinomas at the genetic level and immune cell infiltration levels. More cases and in deep research need to be done to further understand this rare carcinoma.

Whole-Exome Screening and Analysis of Signaling Pathways in Multiple Endocrine Neoplasia Type 1 Patients with Different Outcomes: Insights into Cellular Mechanisms and Possible Functional Implications

Multiple endocrine neoplasia type 1 (MEN1) is a syndrome characterized by tumors in multiple organs. Although being a dominantly inherited monogenic disease, disease phenotypes are unpredictable and differ even among members of the same family. There is growing evidence for the role of modifier genes in the alteration of the course of this disease. However, genome-wide screening data are still lacking. In our study, we addressed the different outcomes of the disease, focusing on pituitary and adrenocortical tumors. By means of exome sequencing we identified the affected signaling pathways that segregated with those symptoms. Most significantly, we identified damaging alterations in numerous structural genes responsible for cell adhesion and migration. Additionally, in the case of pituitary tumors, genes related to neuronal function, survival, and morphogenesis were repeatedly identified, while in patients with adrenocortical tumors, TLR10, which is involved in the regulation of the innate immunity, was commonly modified. Our data show that using exome screening, it is possible to find signatures which correlate with the given clinical MEN1 outcomes, providing evidence that studies addressing modifier effects in MEN1 are reasonable.

Comprehensive review and expanding the genetic landscape of Cornelia-de-Lange spectrum: insights from novel mutations and skin biopsy in exome-negative cases

BACKGROUND

Cornelia de Lange Syndrome (CdLS) is a rare genetic disorder characterized by a range of physical, cognitive, and behavioral abnormalities. This study aimed to perform a comprehensive review of the literature on CdLS and investigate two cases of CdLS with distinct phenotypes that underwent WES to aid in their diagnosis.

METHODS

We conducted a comprehensive review of the literature on CdLS along with performing whole-exome sequencing on two CdLS patients with distinct phenotypes, followed by Sanger sequencing validation and in-silico analysis.

RESULTS

The first case exhibited a classic CdLS phenotype, but the initial WES analysis of blood-derived DNA failed to identify any mutations in CdLS-related genes. However, a subsequent WES analysis of skin-derived DNA revealed a novel heterozygous mutation in the NIPBL gene (NM_133433.4:c.6534_6535del, p.Met2178Ilefs*8). The second case was presented with a non-classic CdLS phenotype, and WES analysis of blood-derived DNA identified a heterozygous missense variant in the SMC1A gene (NM_006306.4:c.2320G>A, p.Asp774Asn).

CONCLUSIONS

The study shows the importance of considering mosaicism in classic CdLS cases and the value of WES for identifying genetic defects. These findings contribute to our understanding of CdLS genetics and underscore the need for comprehensive genetic testing to enhance the diagnosis and management of CdLS patients.

Using multi-scale genomics to associate poorly annotated genes with rare diseases

BACKGROUND

Next-generation sequencing (NGS) has significantly transformed the landscape of identifying disease-causing genes associated with genetic disorders. However, a substantial portion of sequenced patients remains undiagnosed. This may be attributed not only to the challenges posed by harder-to-detect variants, such as non-coding and structural variations but also to the existence of variants in genes not previously associated with the patient's clinical phenotype. This study introduces EvORanker, an algorithm that integrates unbiased data from 1,028 eukaryotic genomes to link mutated genes to clinical phenotypes.

METHODS

EvORanker utilizes clinical data, multi-scale phylogenetic profiling, and other omics data to prioritize disease-associated genes. It was evaluated on solved exomes and simulated genomes, compared with existing methods, and applied to 6260 knockout genes with mouse phenotypes lacking human associations. Additionally, EvORanker was made accessible as a user-friendly web tool.

RESULTS

In the analyzed exomic cohort, EvORanker accurately identified the "true" disease gene as the top candidate in 69% of cases and within the top 5 candidates in 95% of cases, consistent with results from the simulated dataset. Notably, EvORanker outperformed existing methods, particularly for poorly annotated genes. In the case of the 6260 knockout genes with mouse phenotypes, EvORanker linked 41% of these genes to observed human disease phenotypes. Furthermore, in two unsolved cases, EvORanker successfully identified DLGAP2 and LPCAT3 as disease candidates for previously uncharacterized genetic syndromes.

CONCLUSIONS

We highlight clade-based phylogenetic profiling as a powerful systematic approach for prioritizing potential disease genes. Our study showcases the efficacy of EvORanker in associating poorly annotated genes to disease phenotypes observed in patients. The EvORanker server is freely available at https://ccanavati.shinyapps.io/EvORanker/ .

HAT: de novo variant calling for highly accurate short-read and long-read sequencing data

MOTIVATION

de novo variants (DNVs) are variants that are present in offspring but not in their parents. DNVs are both important for examining mutation rates as well as in the identification of disease-related variation. While efforts have been made to call DNVs, calling of DNVs is still challenging from parent-child sequenced trio data. We developed Hare And Tortoise (HAT) as an automated DNV detection workflow for highly accurate short-read and long-read sequencing data. Reliable detection of DNVs is important for human genomics and HAT addresses this need.

RESULTS

HAT is a computational workflow that begins with aligned read data (i.e. CRAM or BAM) from a parent-child sequenced trio and outputs DNVs. HAT detects high-quality DNVs from Illumina short-read whole-exome sequencing, Illumina short-read whole-genome sequencing, and highly accurate PacBio HiFi long-read whole-genome sequencing data. The quality of these DNVs is high based on a series of quality metrics including number of DNVs per individual, percent of DNVs at CpG sites, and percent of DNVs phased to the paternal chromosome of origin.

AVAILABILITY AND IMPLEMENTATION

https://github.com/TNTurnerLab/HAT.

Carrier frequency estimation of pathogenic variants of autosomal recessive and X-linked recessive mendelian disorders using exome sequencing data in 1,642 Thais

BACKGROUND

People with autosomal recessive disorders often were born without awareness of the carrier status of their parents. The American College of Medical Genetics and Genomics (ACMG) recommends screening 113 genes known to cause autosomal recessive and X-linked conditions in couples seeking to learn about their risk of having children with these disorders to have an appropriate reproductive plan.

METHODS

We analyzed the exome sequencing data of 1,642 unrelated Thai individuals to identify the pathogenic variant (PV) frequencies in genes recommended by ACMG.

RESULTS

In the 113 ACMG-recommended genes, 165 PV and likely PVs in 60 genes of 559 exomes (34%, 559/1642) were identified. The carrier rate was increased to 39% when glucose-6-phosphate dehydrogenase (G6PD) was added. The carrier rate was still as high as 14.7% when thalassemia and hemoglobinopathies were excluded. In addition to thalassemia, hemoglobinopathies, and G6PD deficiency, carrier frequencies of > 1% were found for Gaucher disease, primary hyperoxaluria, Pendred syndrome, and Wilson disease. Nearly 2% of the couples were at risk of having offsprings with the tested autosomal recessive conditions.

CONCLUSIONS

Based on the study samples, the expanded carrier screening, which specifically targeted common autosomal recessive conditions in Thai individuals, will benefit clinical outcomes, regarding preconception/prenatal genetic carrier screening.

Phenotype-driven reanalysis reveals five novel pathogenic variants in 40 exome-negative families with Charcot-Marie-Tooth Disease

BACKGROUND

To identify genetic causes in 40 whole exome sequencing (WES)-negative Charcot-Marie-Tooth (CMT) families and provide a summary of the clinical and genetic features of the diagnosed patients.

METHODS

The clinical information and sequencing data of 40 WES-negative families out of 131 CMT families were collected, and phenotype-driven reanalysis was conducted using the Exomiser software.

RESULTS

The molecular diagnosis was regained in 4 families, increasing the overall diagnosis rate by 3.0%. One family with adolescent-onset pure CMT1 was diagnosed [POLR3B: c.2810G>A (p.R937Q)] due to the novel genotype-phenotype association. One infantile-onset, severe CMT1 family with deep sensory disturbance was diagnosed by screening the BAM file and harbored c.1174C>T (p.R392*) and 875_927delinsCTGCCCACTCTGCCCACTCTGCCCACTCTG (p.V292Afs53) of PRX. Two families were diagnosed due to characteristic phenotypes, including an infantile-onset ICMT family with renal dysfunction harboring c.213_233delinsGAGGAGCA (p.S72Rfs34) of INF2 and an adolescent-onset CMT2 family with optic atrophy harboring c.560C>T (p.P187L) and c.616A>G (p.K206E) of SLC25A46. According to the American College of Medical Genetics and Genomics (ACMG) guidelines, the variants of POLR3B and SLC25A46 were classified as likely pathogenic, and the variants of INF2 and PRX were pathogenic. All these variants were first reported worldwide except for p.R392* of PRX.

CONCLUSIONS

We identified five novel pathogenic variants in POLR3B, PRX, INF2, and SLC25A46, which broaden their phenotypic and genotypic spectrums. Regular phenotype-driven reanalysis is a powerful strategy for increasing the diagnostic yield of WES-negative CMT patients, and long-term follow-up and screening BAM files for contiguous deletion and missense variants are both essential for reanalysis.

Whole exome germline sequencing in early-onset prostate cancer patients: Genomic findings and clinical outcomes

BACKGROUND

Whole exome sequencing (WES) furthered our understanding of various tumors. We assessed the occurrence of germline likely pathogenic/pathogenic (LP/P) variants, disease features, and clinical outcomes in early-onset prostate cancer.

METHODS

This retrospective study (N = 134) included consecutive prostate cancer patients who donated blood samples for research purposes to the Kahn-Sagol-Maccabi biobank. Patients diagnosed at >65 years were excluded. Clinical characteristics were extracted from the medical records. Germline WES was performed with analysis reporting on oncogenetic, two immunogenic, and a secondary minimum list panels (121, 468, 76, and 59 genes, respectively).

RESULTS

Median age at diagnosis was 61 (range 46-65) years; 131 (98%) were diagnosed with local disease. The median follow-up time from diagnosis was 14 (range <1-25) years. Of the patients with local disease, 32 (24%) and 10 (8%) had biochemical and distant recurrences, respectively. Twenty-five patients (19%) had ≥1 additional cancer (excluding non-melanoma skin cancer), most frequently bladder (6), colorectal (5), and lymphoma (5). Seven (5%) deaths were reported, with only one related to prostate cancer. LP/P variants were identified in 8 patients (6%), all in genes from the oncogenetic panel: ATM, BRCA1 (in two patients), BRCA2 (in two patients), HOXB13, MUTYH, and MYH7. Of these eight patients, with a median follow-up of 7 years (range <1-15), two (25%) had biochemical recurrences, one had (12.5%) distant recurrence, and no deaths were reported.

CONCLUSIONS

In this cohort of 134 early-onset prostate cancer patients, we identified germline LP/P variants in an oncogenetic panel in 6% of participants, with no unique clinical outcome.

Intestinal-Type Adenocarcinoma in Head and Neck: Dissecting Oncogenic Gene Alterations Through Whole Transcriptome and Exome Analysis

Adenocarcinomas of the nasal/paranasal sinuses are uncommon, but intestinal-type adenocarcinomas (ITACs) are important. Due to the rarity of these tumors, their molecular profile is not well known. To further investigate the molecular profile and find potential oncogenic drivers, we compared the whole transcriptome and exome of ITACs at different anatomic locations in the head and neck. Twenty-one head and neck adenocarcinomas were used in this study, divided into 10 sinonasal adenocarcinomas (SNT) and 11 extrasinonasal (T) head and neck adenocarcinomas according to anatomic location and histology. Tumor samples along with normal mucosa were microdissected from formalin-fixed, paraffin-embedded samples, and RNA and DNA were subjected to whole-transcriptome and -exome shotgun sequencing. Analysis of ITACs at sinonasal locations showed 410 subtype-specific differentially expressed (DE) genes and noncoding transcripts compared with the group of other anatomic locations, with 2909 subtype-specific DE genes. The groups shared 872 genes, with 17 highly different or opposing DE genes. Whole-exome mutation analysis revealed the gene MLL3 (KMT2C) to be exhibiting the most frequent loss-of-function mutations in all adenocarcinomas investigated. The results suggest that the head and neck ITACs investigated were mainly caused by loss-of-function mutations in MLL3 that disabled chromatin methylation and remodeling of all MLL3-targeted enhancers in the tumors. This changed the activity of multiple genes/gene clusters, supporting oncogenicity mostly via pathways of signaling, dedifferentiation, proliferation, migration, and immune and inflammatory deregulation, indicating a truly epigenetic event as the root cause for the heterogenous diversity of these enteric types of cancer. The data of this study form the basis for understanding cell fate determination and cellular homeostasis in the normal respiratory mucosa at different anatomic sites and show the contribution of different mucosal components to the etiology/molecular pathology of ITAC.

Whole exome sequencing reveals novel candidate variants for endometriosis utilizing multiple affected members in a single family

BACKGROUND

Endometriosis is an estrogen-dependent, chronic inflammatory disease that affects 10% of women during the reproductive ages. Despite the estimated 50% heritability for the condition, only 26% was associated with common genetic variants. Thus, necessity of identifying rare variants for the missing heritability is implicated in the literature. Therefore, our study aimed to identify novel rare genetic variants involved in the pathogenesis of endometriosis utilizing a family of multiple affected members.

METHODS

A family composed of four affected women along with their two unaffected mothers were recruited at a single gynecology and infertility clinic specialized in endometriosis. All patients presented with endometriomas, which was visualized by transvaginal ultrasonography. Two affected individuals had received laparoscopic endometrioma excision and therefore were diagnosed with recurrent disease. One mother had a history of endometrial serous adenocarcinoma (ESC) for which she underwent hysterectomy with bilateral oophorectomy. Three endometriosis cases were whole exome sequenced on Illumina NextSeq 550 platform with an average of 90% coverage. Candidate genes were confirmed by Sanger sequencing and followed-up with family segregation.

RESULTS

Novel rare variants were identified in TNFRSF1B (NM_001066.3: c.1072G>A, p.(Ala358Thr)) and GEN1 (NM_001130009.3: c.1574C>T, p.(Ser525Leu)) as possible genetic causes of endometriosis. A third novel rare variant was identified in CRABP1 (NM_004378.3:c.54G>C, p.(Glu18Asp)) only on the mother with ESC history and her daughters.

CONCLUSION

Novel candidate genetic variants that might contribute to endometriosis were suggested that need replication through independent cohorts or validation by functional studies. The family has also received genetic counseling and that the affected daughters are on clinical follow-up, accordingly.

Exome-informed formulations of food proteins enhance body growth and feed conversion efficiency in ad libitum-fed mice

Meeting requirements for dietary proteins, especially of essential amino acids (EAAs), is critical for the life-long health of living organisms. However, defining EAA targets for preparing biologically-matched nutrition that satisfies metabolic requirements for protein remains challenging. Previous research has shown the advantages of 'exome matching' in representing the specific requirement of dietary AAs, where the target dietary AA profile was derived from in silico translation of the genome of an organism, specifically responsible for protein expression (the 'exome'). However, past studies have assessed these effects in only one sex, for few parameters (body mass and composition), and have used purified diets in which protein is supplied as a mixture of individual AAs. Here, for the first time, we utilise a computational method to guide the formulation of custom protein blends and test if exome matching can be achieved at the intact protein level, through blending standard protein ingredients, ultimately leading to optimal growth, longevity and reproductive function. Mice were provided ad libitum (ad lib) access to one of the four iso-energetic protein-limited diets, two matched and two mis-matched to the mouse exome target, and fed at a fixed protein energy level of 6.2%. During or following 13-weeks of feeding, the food intake, body growth, composition and reproductive functions were measured. Compared to the two mis-matched diets, male and female animals on the exome-matched diet with protein digestibility correction applied, exhibited significantly improved growth rates and final body mass. The feed conversion efficiency in the same diet was also increased by 62% and 40% over the worst diets for males and females, respectively. Male, not female, exhibited higher accretion of lean body mass with the matched, digestibility-corrected diet. All reproductive function measures in both sexes were comparable among diets, with the exception of testicular daily sperm production in males, which was higher in the two matched diets versus the mis-matched diets. The results collectively demonstrate the pronounced advantages of exome-matching in supporting body growth and improving feed conversion efficiency in both sexes. However, the potential impact of this approach in enhancing fertility needs further investigation.

A comprehensive assessment of exome capture methods for RNA sequencing of formalin-fixed and paraffin-embedded samples

RNA-Seq analysis of Formalin-Fixed and Paraffin-Embedded (FFPE) samples has emerged as a highly effective approach and is increasingly being used in clinical research and drug development. However, the processing and storage of FFPE samples are known to cause extensive degradation of RNAs, which limits the discovery of gene expression or gene fusion-based biomarkers using RNA sequencing, particularly methods reliant on Poly(A) enrichment. Recently, researchers have developed an exome targeted RNA-Seq methodology that utilizes biotinylated oligonucleotide probes to enrich RNA transcripts of interest, which could overcome these limitations. Nevertheless, the standardization of this experimental framework, including probe designs, sample multiplexing, sequencing read length, and bioinformatic pipelines, remains an essential requirement. In this study, we conducted a comprehensive comparison of three main commercially available exome capture kits and evaluated key experimental parameters, to provide the overview of the advantages and limitations associated with the selection of library preparation protocols and sequencing platforms. The results provide valuable insights into the best practices for obtaining high-quality data from FFPE samples.

Beyond the exome: utility of long-read whole genome sequencing in exome-negative autosomal recessive diseases

BACKGROUND

Long-read whole genome sequencing (lrWGS) has the potential to address the technical limitations of exome sequencing in ways not possible by short-read WGS. However, its utility in autosomal recessive Mendelian diseases is largely unknown.

METHODS

In a cohort of 34 families in which the suspected autosomal recessive diseases remained undiagnosed by exome sequencing, lrWGS was performed on the Pacific Bioscience Sequel IIe platform.

RESULTS

Likely causal variants were identified in 13 (38%) of the cohort. These include (1) a homozygous splicing SV in TYMS as a novel candidate gene for lethal neonatal lactic acidosis, (2) a homozygous non-coding SV that we propose impacts STK25 expression and causes a novel neurodevelopmental disorder, (3) a compound heterozygous SV in RP1L1 with complex inheritance pattern in a family with inherited retinal disease, (4) homozygous deep intronic variants in LEMD2 and SNAP91 as novel candidate genes for neurodevelopmental disorders in two families, and (5) a promoter SNV in SLC4A4 causing non-syndromic band keratopathy. Surprisingly, we also encountered causal variants that could have been identified by short-read exome sequencing in 7 families. The latter highlight scenarios that are especially challenging at the interpretation level.

CONCLUSIONS

Our data highlight the continued need to address the interpretation challenges in parallel with efforts to improve the sequencing technology itself. We propose a path forward for the implementation of lrWGS sequencing in the setting of autosomal recessive diseases in a way that maximizes its utility.

A novel 268 kb deletion combined with a splicing variant in IL7R causes of severe combined immunodeficiency in a Chinese family: a case report

BACKGROUND

Severe combined immunodeficiency (SCID) is a group of fatal primary immunodeficiencies characterized by the severe impairment of T-cell differentiation. IL7R deficiency is a rare form of SCID that usually presents in the first months of life with severe and opportunistic infections, failure to thrive, and a high risk of mortality unless treated. Although recent improvements in early diagnosis have been achieved through newborn screening, few IL7R-related SCID patients had been reported in the Chinese population.

CASE PRESENTATION

Here, we retrospectively analyzed a case of SCID in a 5-month-old girl with symptoms, including severe T-cell depletion, recurrent fever, oral ulcers, pneumonia, hepatosplenomegaly, bone marrow hemophagocytosis, and bacterial and viral infections. Whole-exome sequencing (WES), quantitative PCR (qPCR), and chromosome microarray analysis (CMA) were performed to identify the patient's genetic etiology. We identified a 268 kb deletion and a splicing variant, c.221 + 1G > A, in the proband. These two variants of IL7R were inherited from the father and mother.

CONCLUSIONS

To our knowledge, this is the first report of whole IL7R gene deletion in combination with a pathogenic splicing variant in a patient with SCID. This deletion also expands the pathogenic variation spectrum of SCID caused by IL7R. The incorporation of exome-based copy number variant analysis makes WES a powerful molecular diagnostic technique for the clinical diagnosis of pediatric patients.

Benchmarking and improving the performance of variant-calling pipelines with RecallME

MOTIVATION

The steady increment of Whole Genome/Exome sequencing and the development of novel Next Generation Sequencing-based gene panels requires continuous testing and validation of variant calling (VC) pipelines and the detection of sequencing-related issues to be maintained up-to-date and feasible for the clinical settings. State of the art tools are reliable when used to compute standard performance metrics. However, the need for an automated software to discriminate between bioinformatic and sequencing issues and to optimize VC parameters remains unmet.

RESULTS

The aim of the current work is to present RecallME, a bioinformatic suite that tracks down difficult-to-detect variants as insertions and deletions in highly repetitive regions, thus providing the maximum reachable recall for both single nucleotide variants and small insertion and deletions and to precisely guide the user in the pipeline optimization process.

AVAILABILITY AND IMPLEMENTATION

Source code is freely available under MIT license at https://github.com/mazzalab-ieo/recallme. RecallME web application is available at https://translational-oncology-lab.shinyapps.io/recallme/. To use RecallME, users must obtain a license for ANNOVAR by themselves.

Exome analysis focusing on epilepsy-related genes in children and adults with sudden unexplained death

PURPOSE

Genetic studies in sudden infant death syndrome (SIDS) and sudden unexplained death (SUD) cohorts have indicated that cardiovascular diseases might have contributed to sudden unexpected death in 20-35 % of autopsy-negative cases. Sudden unexpected death can also occur in people with epilepsy, termed as sudden unexpected death in epilepsy (SUDEP). The pathophysiological mechanisms of SUDEP are not well understood, but are likely multifactorial, including seizure-induced hypoventilation and arrhythmias as well as genetic risk factors. The sudden death of some of the SIDS/SUD victims might also be explained by genetic epilepsy, therefore this study aimed to expand the post-mortem genetic analysis of SIDS/SUD cases to epilepsy-related genes.

METHODS

Existing whole-exome sequencing data from our 155 SIDS and 45 SUD cases were analyzed, with a focus on 365 epilepsy-related genes. Nine of the SUD victims had a known medical history of epilepsy, seizures or other underlying neurological conditions and were therefore classified as SUDEP cases.

RESULTS

In our SIDS and SUD cohorts, we found epilepsy-related pathogenic/likely pathogenic variants in the genes OPA1, RAI1, SCN3A, SCN5A and TSC2.

CONCLUSION

Post-mortem analysis of epilepsy-related genes identified potentially disease-causing variants that might have contributed to the sudden death events in our SIDS/SUD cases. However, the interpretation of identified variants remains challenging and often changes over time as more data is gathered. Overall, this study contributes insight in potentially pathophysiological epilepsy-related mechanisms in SIDS, SUD and SUDEP victims and underlines the importance of sensible counselling on the risk and preventive measures in genetic epilepsy.

A Novel Truncating Mutation in PAX1 Gene Causes Otofaciocervical Syndrome Without Immunodeficiency

Otofaciocervical syndrome (OTFCS) is a rare genetic disorder of both autosomal recessive and autosomal dominant patterns of inheritance. It is caused by biallelic or monoallelic mutations in PAX1 or EYA1 genes, respectively. Here, we report an OTFCS2 female patient of 1st consanguineous healthy parents. She manifested facial dysmorphism, hearing loss, intellectual disability (ID), and delayed language development (DLD) as the main clinical phenotype. The novel homozygous variant c.1212dup (p.Gly405Argfs*51) in the PAX1 gene was identified by whole exome sequencing (WES), and family segregation confirmed the heterozygous status of the mutation in the parents using the Sanger sequencing. The study recorded a novel PAX1 variant representing the sixth report of OTFCS2 worldwide and the first Egyptian study expanding the geographic area where the disorder was confined.

Solving the unsolved genetic epilepsies: Current and future perspectives

Many patients with epilepsy undergo exome or genome sequencing as part of a diagnostic workup; however, many remain genetically unsolved. There are various factors that account for negative results in exome/genome sequencing for patients with epilepsy: (1) the underlying cause is not genetic; (2) there is a complex polygenic explanation; (3) the illness is monogenic but the causative gene remains to be linked to a human disorder; (4) family segregation with reduced penetrance; (5) somatic mosaicism or the complexity of, for example, a structural rearrangement; or (6) limited knowledge or diagnostic tools that hinder the proper classification of a variant, resulting in its designation as a variant of unknown significance. The objective of this review is to outline some of the diagnostic options that lie beyond the exome/genome, and that might become clinically relevant within the foreseeable future. These options include: (1) re-analysis of older exome/genome data as knowledge increases or symptoms change; (2) looking for somatic mosaicism or long-read sequencing to detect low-complexity repeat variants or specific structural variants missed by traditional exome/genome sequencing; (3) exploration of the non-coding genome including disruption of topologically associated domains, long range non-coding RNA, or other regulatory elements; and finally (4) transcriptomics, DNA methylation signatures, and metabolomics as complementary diagnostic methods that may be used in the assessment of variants of unknown significance. Some of these tools are currently not integrated into standard diagnostic workup. However, it is reasonable to expect that they will become increasingly available and improve current diagnostic capabilities, thereby enabling precision diagnosis in patients who are currently undiagnosed.

Informing a value care model: lessons from an integrated adult neurogenomics clinic

BACKGROUND

Advances in genomics provide improved opportunities for diagnosis of complex neurogenetic disorders, yet the optimal approach to translate these benefits to the outpatient clinic is unclear.

AIMS

We retrospectively reviewed referral indications and outcomes of an integrated multidisciplinary team (MDT) clinic pathway for adults with suspected neurogenetic disorders. The associated cost implications were estimated.

METHODS

Consecutive patients who attended the neurogenomics clinic from January 2017 to April 2020 were included. The clinic comprised neurologists, clinical geneticists and genetic counsellors, who assessed each patient concurrently.

RESULTS

Ninety-nine new patients were referred spanning 45 different clinical diagnoses. Following MDT clinical assessment, 23% (23/99) of referral diagnoses were revised prior to molecular testing. Eighty-one patients (82%) underwent genetic testing, including 43 exome-based panels, 15 whole-genome sequencing, 14 single gene tests, 27 repeat-primed polymerase chain reaction testing and two chromosomal microarrays. Overall, 33/99 patients (33%) received a diagnosis, either a molecular diagnosis (n = 24, of which 22 were diagnostic and two were predictive) or a clinical diagnosis (n = 9). Of the clinical diagnosis cohort, five patients received a diagnosis without molecular testing and four patients whose negative testing (one diagnostic and three predictive) allowed exclusion of genetic differentials and, hence, confirmation of clinical diagnoses. The diagnostic rate following MDT and diagnostic testing was 30% (28/94), excluding the five predictive testing cases. MDT assessment aligned with eventual molecular diagnoses in 96% of cases. The estimated average costs were AU$1386 per patient undergoing MDT assessment and AU$4159 per diagnosis achieved.

CONCLUSIONS

We present an integrated multidisciplinary neurogenomics clinic pathway providing a diagnostic yield of 33% (30% excluding predictive testing cases), with costing implications. The relatively high diagnostic yield may be attributed to multidisciplinary input integrating accurate phenotyping of complex disorders and interpretation of genomic findings.

Skim exome capture genotyping in wheat

Next-generation sequencing (NGS) technology advancements continue to reduce the cost of high-throughput genome-wide genotyping for breeding and genetics research. Skim sequencing, which surveys the entire genome at low coverage, has become feasible for quantitative trait locus (QTL) mapping and genomic selection in various crops. However, the genome complexity of allopolyploid crops such as wheat (Triticum aestivum L.) still poses a significant challenge for genome-wide genotyping. Targeted sequencing of the protein-coding regions (i.e., exome) reduces sequencing costs compared to whole genome re-sequencing and can be used for marker discovery and genotyping. We developed a method called skim exome capture (SEC) that combines the strengths of these existing technologies and produces targeted genotyping data while decreasing the cost on a per-sample basis compared to traditional exome capture. Specifically, we fragmented genomic DNA using a tagmentation approach, then enriched those fragments for the low-copy genic portion of the genome using commercial wheat exome baits and multiplexed the sequencing at different levels to achieve desired coverage. We demonstrated that for a library of 48 samples, ∼7-8× target coverage was sufficient for high-quality variant detection. For higher multiplexing levels of 528 and 1056 samples per library, we achieved an average coverage of 0.76× and 0.32×, respectively. Combining these lower coverage SEC sequencing data with genotype imputation using a customized wheat practical haplotype graph database that we developed, we identified hundreds of thousands of high-quality genic variants across the genome. The SEC method can be used for high-resolution QTL mapping, genome-wide association studies, genomic selection, and other downstream applications.

A loss of function variant in AGPAT3 underlies intellectual disability and retinitis pigmentosa (IDRP) syndrome

Intellectual disability (ID) and retinal dystrophy (RD) are the frequently found features of multiple syndromes involving additional systemic manifestations. Here, we studied a family with four members presenting severe ID and retinitis pigmentosa (RP). Using genome wide genotyping and exome sequencing, we identified a nonsense variant c.747 C > A (p.Tyr249Ter) in exon 7 of AGPAT3 which co-segregates with the disease phenotype. Western blot analysis of overexpressed WT and mutant AGPAT3 in HEK293T cells showed the absence of AGPAT3, suggesting instability of the truncated protein. Knockdown of Agpat3 in the embryonic mouse brain caused marked deficits in neuronal migration, strongly suggesting that reduced expression of AGPAT3 affects neuronal function. Altogether, our data indicates that AGPAT3 activity is essential for neuronal functioning and loss of its activity probably causes intellectual disability and retinitis pigmentosa (IDRP) syndrome.

Targeted gene panel provides advantages over whole-exome sequencing for diagnosing obesity and diabetes mellitus

A small fraction of patients diagnosed with obesity or diabetes mellitus has an underlying monogenic cause. Here, we constructed a targeted gene panel consisting of 83 genes reported to be causative for monogenic obesity or diabetes. We performed this panel in 481 patients to detect causative variants and compared these results with whole-exome sequencing (WES) data available for 146 of these patients. The coverage of targeted gene panel sequencing was significantly higher than that of WES. The diagnostic yield in patients sequenced by the panel was 32.9% with subsequent WES leading to three additional diagnoses with two novel genes. In total, 178 variants in 83 genes were detected in 146 patients by targeted sequencing. Three of the 178 variants were missed by WES, although the WES-only approach had a similar diagnostic yield. For the 335 samples only receiving targeted sequencing, the diagnostic yield was 32.2%. In conclusion, taking into account the lower costs, shorter turnaround time, and higher quality of data, targeted sequencing is a more effective screening method for monogenic obesity and diabetes compared to WES. Therefore, this approach could be routinely established and used as a first-tier test in clinical practice for specific patients.