Genomic variant annotation and prioritization with ANNOVAR and wANNOVAR

Comparative analyses of Pleurotus pulmonarius mitochondrial genomes reveal two major lineages of mini oyster mushroom cultivars

Pleurotus pulmonarius, commonly known as the mini oyster mushroom, is highly esteemed for its crisp texture and umami flavor. Limited genetic diversity among P. pulmonarius cultivars raises concerns regarding its sustainable industrial production. To delve into the maternal genetic diversity of the principal P. pulmonarius cultivars, 36 cultivars and five wild isolates were subjected to de novo sequencing and assembly to generate high-quality mitogenome sequences. The P. pulmonarius mitogenomes had lengths ranging from 69,096 to 72,905 base pairs. The mitogenome sizes of P. pulmonarius and those of other mushroom species in the Pleurotus genus showed a significant positive correlation with the counts of LAGLIDAG and GIY-YIG homing endonucleases encoded by intronic open reading frames. A comparison of gene arrangements revealed an inversion of a fragment containing atp9-nad3-nad2 between P. pulmonarius and P. ostreatus. The mitogenomes of P. pulmonarius were clustered into three distinct clades, two of which were crowded with commercial cultivars. Clade I, all of which possess an inserted dpo gene, shared a maternal origin linked to an ancestral cultivar from Taiwan. Primers were designed to target the dpo gene, potentially safeguarding intellectual property rights. The wild isolates in Clade III exhibited more divergent mitogenomes, rendering them valuable for breeding.

Genomic and transcriptomic landscape of human gastrointestinal stromal tumors

Gastrointestinal stromal tumor (GISTs) are clinically heterogenous exhibiting varying degrees of disease aggressiveness in individual patients. We comprehensively describe the genomic and transcriptomic landscape of a cohort of 117 GISTs including 31 low-risk, 18 intermediate-risk, 29 high-risk, 34 metastatic and 5 neoadjuvant GISTs from 105 patients. GISTs have notably low tumor mutation burden but widespread copy number variations. Aggressive GISTs harbor remarkably more genomic aberrations than low-/intermediate-risk GISTs. Complex genomic alterations, chromothripsis and kataegis, occur selectively in aggressive GISTs. Despite the paucity of mutations, recurrent inactivating YLPM1 mutations are identified (10.3%, 7 of 68 patients), enriched in high-risk/metastatic GIST and functional study further demonstrates YLPM1 inactivation promotes GIST proliferation, growth and oxidative phosphorylation. Spatially and temporally separated GISTs from individual patients demonstrate complex tumor heterogeneity in metastatic GISTs. Finally, four prominent subtypes are proposed with different genomic features, expression profiles, immune characteristics, clinical characteristics and subtype-specific treatment strategies. This large-scale analysis depicts the landscape and provides further insights into GIST pathogenesis and precise treatment.

Identification of novel BCL11A variant in a patient with developmental delay and behavioural differences

BACKGROUND

The BCL11A gene is involved in disorders including intellectual disability syndrome (IDS), encodes a zinc finger protein highly expressed in haematopoietic and brain and acts as a transcriptional repressor of foetal haemoglobin (HbF). De novo variants in BCL11A have been associated with IDS, which is characterized by developmental delays, autism spectrum disorder (ASD) and speech and language delays. The reports of BCL11A gene variants are still limited worldwide, and additional cases are needed to expand the variant and phenotype spectrums.

METHODS

The patient is a 5-year-old girl who was hospitalized due to developmental delays. After analysing her clinical and pathological characterizations, whole-exome sequencing (WES) was performed for pathogenic genetic variants of developmental delay and behavioural differences. Candidate variant in BCL11A gene was identified and further validated by Sanger sequencing.

RESULTS

A heterozygous variant, c.1442delA (p.Glu481Glyfs*25), was identified in exon 4 of the BCL11A gene through WES. This variant results in a truncated expression of the encoded protein. This de novo variant was confirmed by Sanger sequencing. The language delay and behavioural differences were prominent in our patient.

CONCLUSION

Our finding demonstrates that BCL11A variant may cause developmental delay and behavioural differences, expanding the genetic spectrum of the BCL11A gene.

Genomic insights into the contribution of de novo lipogenesis to intramuscular fat deposition in chicken

INTRODUCTION

The proportion of animal based foods in daily diet of consumers is constantly increasing, with chicken being highly favored due to its high protein and low fat characteristics. The consumption of chicken around the world is steadily increasing. Intramuscular fat (IMF) is a key indicator affecting meat quality.

OBJECT

High IMF content can contribute to improve the quality of chicken meat. The regulatory mechanism of IMF deposition in chicken is poorly understood, so its complete elucidation is essential to improve chicken meat quality.

METHOD

Here, we performed whole genome resequencing on 516 yellow feather chickens and single-cell RNA sequencing on 3 63-day-old female JXY chickens. In addition, transcriptome sequencing techniques were also performed on breast muscle tissue of JXY chickens at different developmental stages. And 13C isotope tracing technique was applied.

RESULTS

In this study, a large-scale genetic analysis of an IMF-selected population and a control population identified fatty acid synthase (FASN) as a key gene for improving IMF content. Also, contrary to conventional view, de novo lipogenesis (DNL) was deemed to be an important contributor to IMF deposition. As expected, further analyses by isotope tracing and other techniques, confirmed that DNL mainly occurs in myocytes, contributing about 40% of the total fatty acids through the regulation of FASN, using the available FAs as substrates. Additionally, we also identified a relevant causal mutation in the FASN gene with effects on FA composition.

CONCLUSION

These findings contribute to the understanding of fat metabolism in muscle tissue of poultry, and provide the feasible strategy for the production of high-quality chicken meat.

Analysis of genetic structure and identification of important genes associated with muscle growth in Fujian Muscovy duck

Fujian Muscovy duck is a well-known meat waterfowl in Fujian Province due to its high meat production, superior breeding potential, and strong resistance. To fully explore the genetic characteristics of these advantages, Fujian black Muscovy duck and white Muscovy duck were used for whole-genome re-sequencing and transcriptome analyses. Population structure analysis showed significant differentiation between the two feather strains. Runs of homozygosity analysis indicated a stronger artificial influence on the black-feathered strain, with ROH island genes notably enriched in muscle tissue-related terms and pathways. Selective sweep and transcriptome analysis revealed a significant enrichment of genes linked to muscle tissue and muscle fiber-related terms and pathways. Key candidate genes identified, such as MEF2C, MYOZ2, and METTL21C, are believed to play crucial roles in meat production in Fujian Muscovy duck. This study offers a new perspective on improving meat production in Fujian Muscovy duck, which can benefit breeding strategies and production management.

Clinical relevance of protein-truncating variants of germline DNA repair genes in prostate cancer

BACKGROUND

Interpreting genetic variants remains a challenge in prostate cancer (PCa). Although many annotation tools are available for prioritizing causal variants, the clinical relevance of these variants is rarely studied.

METHODS

We collected a cohort study that included 274 PCa patients from June 2017 to December 2020 and sequenced 19 DNA damage repair (DDR) genes in these patients and explored the clinical consequence of these different approaches. We also examined all-cause and PCa-specific survival in DDR gene mutation carriers compared to non-carriers after androgen receptor (AR)-directed therapy.

RESULTS

We identified 13 variants from 19 DDR genes in a total of 14 (5.1%) patients who had at least one presumed pathogenic mutation using different annotation methods. Four variants were annotated as pathogenic, 11 variants were predicted as protein-truncating variants (PTVs), four variants received proxy-deleterious (Combined Annotation-Dependent Depletion scores of > 30), and only one variant was identified as a pathogenic variant or as having a functional effect by all three methods. PCa patients with PTVs were significantly associated with early onset, high cancer stage, and a worse response to AR-directed treatment. However, patients carrying a proxy-deleterious variant were only associated with a higher T (tumor) stage and N (node) stage than those without such a variant, but not associated with other clinical characteristics. In patients treated with AR-directed therapy, patients with a PTV showed an increased risk of all-cause death (adjusted hazard ratio (aHR) = 3.51, 95% confidence interval (CI): 1.06 ~ 11.56) and PCa-specific death (aHR = 4.49, 95% CI: 1.87 ~ 10.77) compared to non-PTV carriers after adjustment. We were unable to examine gene-specific risks due to the small number of patients.

CONCLUSIONS

PTVs may assist in guiding treatment and early screening in PCa, while population-specific data for pathogenic variants are still being amassed.

Identification of ubiquitin markers for survival and prognosis of ovarian cancer

Ovarian cancer (OC) is one of the most common malignancies and a leading cause of death among women worldwide. The ubiquitin pathway plays an important role in OC development. Using the single nucleotide polymorphism data obtained using the prevalence and dominance strategies, four ubiquitin marker genes were identified according to their expression levels: BARD1, BRCA2, FANCA, and BRCA1. Based on these four genes, a consensus clustering of OC expression data was performed. The significant differences in the survival analysis, ESTIMATE, and CIBERSORT results among the clusters indicated the pivotal role of these four genes in OC development. Of the ubiquitin-representative genes in each cluster, two ubiquitin genes, TOP2A and MYLIP, were identified in the survival risk model after univariate survival, Least Absolute Shrinkage and Selection Operator regression, and multivariate survival analyses. The reliability and robustness of both the training and validation data were confirmed by comparing the significant survival difference between high- and low-risk patients. We further explored the association between our risk model and clinical outcomes as well as predicted potentially interacting drugs. The co-expression network showed clear interactions among the four marker genes and two model genes and between high- and low-risk differentially expressed genes. Significantly enriched genes were found in pathways associated with ion channels, channel activity, and neuroactive ligand-receptor interactions. Therefore, suggesting the involvement of ubiquitin genes in the survival and development of OC through neurohumoral regulation. Our results will provide valuable reference or supplementary information for studies investigating OC diagnosis and therapies.

Whole exome sequencing identifies a novel mutation in Annexin A4 that is associated with recurrent spontaneous abortion

Background

Recurrent spontaneous abortion (RSA) is a multifactorial disease, the exact causes of which are still unknown. Environmental, maternal, and genetic factors have been shown to contribute to this condition. The aim of this study was to investigate the presence of mutations in the ANXA4 gene in patients with RSA.

Methods

Genomic DNA was extracted from 325 patients with RSA and 941 control women with a normal reproductive history for whole-exome sequencing (WES). The detected variants were annotated and filtered, and the pathogenicity of the variants was predicted through the SIFT online tool, functional enrichment analyses, Sanger sequencing validation, prediction of changes in protein structure, and evolutionary conservation analysis. Furthermore, plasmid construction, Western blotting, RT-qPCR, and cell migration, invasion and adhesion assays were used to detect the effects of ANXA4 mutations on protein function.

Results

An ANXA4 mutation (p.G8D) in 1 of the 325 samples from patients with RSA (RSA-219) was identified through WES. This mutation was not detected in 941 controls or included in public databases. Evolutionary conservation analysis revealed that the amino acid residue affected by the mutation (p.G8D) was highly conserved among 13 vertebrate species, and the SIFT program and structural modeling analysis predicted that this mutation was harmful. Furthermore, functional assays revealed that this mutation could inhibit cell migration, invasion and adhesion.

Conclusion

Our study suggests that an unreported novel ANXA4 mutation (p.G8D) plays an important role in the pathogenesis of RSA and may contribute to the genetic diagnosis of RSA.

Large-scale multi-omic analysis identifies noncoding somatic driver mutations and nominates ZFP36L2 as a driver gene for pancreatic ductal adenocarcinoma

Identification of somatic driver mutations in the noncoding genome remains challenging. To comprehensively characterize noncoding driver mutations for pancreatic ductal adenocarcinoma (PDAC), we first created genome-scale maps of accessible chromatin regions (ACRs) and histone modification marks (HMMs) in pancreatic cell lines and purified pancreatic acinar and duct cells. Integration with whole-genome mutation calls from 506 PDACs revealed 314 ACRs/HMMs significantly enriched with 3,614 noncoding somatic mutations (NCSMs). Functional assessment using massively parallel reporter assays (MPRA) identified 178 NCSMs impacting reporter activity (19.45% of those tested). Focused luciferase validation confirmed negative effects on gene regulatory activity for NCSMs near CDKN2A and ZFP36L2. For the latter, CRISPR interference (CRISPRi) further identified ZFP36L2 as a target gene (16.0 - 24.0% reduced expression, P = 0.023-0.0047) with disrupted KLF9 binding likely mediating the effect. Our integrative approach provides a catalog of potentially functional noncoding driver mutations and nominates ZFP36L2 as a PDAC driver gene.

Predicting functional outcome in ischemic stroke patients using genetic, environmental, and clinical factors: a machine learning analysis of population-based prospective cohort study

Ischemic stroke (IS) is a leading cause of adult disability that can severely compromise the quality of life for patients. Accurately predicting the IS functional outcome is crucial for precise risk stratification and effective therapeutic interventions. We developed a predictive model integrating genetic, environmental, and clinical factors using data from 7819 IS patients in the Third China National Stroke Registry. Employing an 80:20 split, we randomly divided the dataset into development and internal validation cohorts. The discrimination and calibration performance of models were evaluated using the area under the receiver operating characteristic curves (AUC) for discrimination and Brier score with calibration curve in the internal validation cohort. We conducted genome-wide association studies (GWAS) in the development cohort, identifying rs11109607 (ANKS1B) as the most significant variant associated with IS functional outcome. We employed principal component analysis to reduce dimensionality on the top 100 significant variants identified by the GWAS, incorporating them as genetic factors in the predictive model. We employed a machine learning algorithm capable of identifying nonlinear relationships to establish predictive models for IS patient functional outcome. The optimal model was the XGBoost model, which outperformed the logistic regression model (AUC 0.818 versus 0.756, P < .05) and significantly improved reclassification efficiency. Our study innovatively incorporated genetic, environmental, and clinical factors for predicting the IS functional outcome in East Asian populations, thereby offering novel insights into IS functional outcome.

Genetic variants associated with response to anti-CGRP monoclonal antibody therapy in a chronic migraine Han Chinese population

BACKGROUND

Anti-calcitonin gene-related peptide (CGRP) monoclonal antibodies have emerged as promising therapeutic options for the treatment of chronic migraine. However, treatment response varies considerably among individuals, suggesting a potential role for genetic factors. This study aimed to identify genetic variants affecting the efficacy of anti-CGRP monoclonal antibody therapy in chronic migraine among the Han Chinese population in Taiwan to enhance treatment precision and to understand the genetic architecture of migraine.

METHODS

We conducted a quantitative trait locus (QTL) association study in patients with chronic migraines from a tertiary medical center in Taiwan using the Taiwan Precision Medicine Array Chip. The patients received fremanezumab or galcanezumab for at least 12 weeks. Treatment efficacy was assessed based on the improvement rate in monthly migraine days. Genetic variants were identified, and their associations with treatment efficacy were examined through quantitative trait loci analysis, linkage disequilibrium studies, and functional annotations using the Gene Ontology database.

RESULTS

Six single nucleotide polymorphisms (SNPs) relative variants were significantly associated with anti-CGRP therapy response (p < 1 × 10- 7): rs116870564, rs75244870, rs56216870, rs12938101, rs74655790, and rs149540851. These variants are located in or near genes, including LRRC4C, ATAD2B, and OXR1, which are involved in neuronal development, DNA-dependent ATPase activity, and oxidation-reduction processes, respectively. The rs116870564 variant in LRRC4C showed the strongest association (β = -0.551, p = 6.65 × 10- 9). The functional impact of these variants is attributed to their regulatory effects on gene expression, which are influenced by intron splicing regulation, transcription factors, and changes in chromatin structure.

CONCLUSION

The identification of key genetic markers associated with response to anti-CGRP therapy emphasizes the importance of genetic variability in treatment efficacy. This could lead to more personalized chronic migraine management strategies and tailored therapeutic approaches based on individual genetic profiles. Further research in larger, diverse populations is warranted to validate these findings and refine our understanding of the role of CGRP in chronic migraine pathophysiology.

TRIAL REGISTRATION

Not applicable.

Whole-genome resequencing and RNA-seq analysis implicates GPR75 as a potential genetic basis related to retarded growth in South China carp (Cyprinus carpio rubrofuscus)

The south China carp (Cyprinus carpio rubrofuscus) is an indigenous and important fish species, widely cultured in south China. However, part of individuals experienced retarded growth, the genetic basis of which has yet to be elucidated. In this study, whole-genome resequencing of 35 fast-growing and 35 retarded-growing south China carp were conducted to identify promising genes associated with retarded growth. Twelve candidate SNPs were detected and annotated to the Gpr75 gene, which has been reported to be related with body weight through regulating insulin homeostasis. RNA-seq analysis of muscle suggested that differentially expressed genes were significantly enriched in the insulin signaling pathway. Additionally, the fasting serum insulin level was significantly lower while the blood glucose level was significantly higher in the retarded-growing group. Our preliminary study provides insights into the genetic basis underlying the retarded growth and may facilitate further genetic improvement of south China carp.

Variant Impact Predictor database (VIPdb), version 2: trends from three decades of genetic variant impact predictors

BACKGROUND

Variant interpretation is essential for identifying patients' disease-causing genetic variants amongst the millions detected in their genomes. Hundreds of Variant Impact Predictors (VIPs), also known as Variant Effect Predictors (VEPs), have been developed for this purpose, with a variety of methodologies and goals. To facilitate the exploration of available VIP options, we have created the Variant Impact Predictor database (VIPdb).

RESULTS

The Variant Impact Predictor database (VIPdb) version 2 presents a collection of VIPs developed over the past three decades, summarizing their characteristics, ClinGen calibrated scores, CAGI assessment results, publication details, access information, and citation patterns. We previously summarized 217 VIPs and their features in VIPdb in 2019. Building upon this foundation, we identified and categorized an additional 190 VIPs, resulting in a total of 407 VIPs in VIPdb version 2. The majority of the VIPs have the capacity to predict the impacts of single nucleotide variants and nonsynonymous variants. More VIPs tailored to predict the impacts of insertions and deletions have been developed since the 2010s. In contrast, relatively few VIPs are dedicated to the prediction of splicing, structural, synonymous, and regulatory variants. The increasing rate of citations to VIPs reflects the ongoing growth in their use, and the evolving trends in citations reveal development in the field and individual methods.

CONCLUSIONS

VIPdb version 2 summarizes 407 VIPs and their features, potentially facilitating VIP exploration for various variant interpretation applications. VIPdb is available at  https://genomeinterpretation.org/vipdb.

Clinical efficacy of tumor organoid-guided cancer therapy for locally advanced unresectable or metastatic breast cancer

Patient-derived organoids (PDOs) may facilitate treatment selection. This retrospective cohort study evaluated the feasibility and clinical benefit of using PDOs to guide personalized treatment in metastatic breast cancer (MBC). Patients diagnosed with MBC were recruited between January 2019 and August 2022. PDOs were established and the efficacy of customized drug panels was determined by measuring cell mortality after drug exposure. Patients receiving organoid-guided treatment (OGT) were matched 1:2 by nearest neighbor propensity scores with patients receiving treatment of physician's choice (TPC). The primary outcome was progression-free survival. Secondary outcomes included objective response rate and disease control rate. Targeted gene sequencing and pathway enrichment analysis were performed. Forty-six PDOs (46 of 51, 90.2%) were generated from 45 MBC patients. PDO drug screening showed an accuracy of 78.4% (95% CI 64.9%-91.9%) in predicting clinical responses. Thirty-six OGT patients were matched to 69 TPC patients. OGT was associated with prolonged median progression-free survival (11.0 months vs. 5.0 months; hazard ratio 0.53 [95% CI 0.33-0.85]; p = .01) and improved disease control (88.9% vs. 63.8%; odd ratio 4.26 [1.44-18.62]) compared with TPC. The objective response rate of both groups was similar. Pathway enrichment analysis in hormone receptor-positive, human epidermal growth factor receptor 2-negative patients demonstrated differentially modulated pathways implicated in DNA repair and transcriptional regulation in those with reduced response to capecitabine/gemcitabine, and pathways associated with cell cycle regulation in those with reduced response to palbociclib. Our study shows that PDO-based functional precision medicine is a feasible and effective strategy for MBC treatment optimization and customization.

Boolean Modeling of Biological Network Applied to Protein-Protein Interaction Network of Autism Patients

Cellular molecules interact with one another in a structured manner, defining a regulatory network topology that describes cellular mechanisms. Genetic mutations alter these networks' pathways, generating complex disorders such as autism spectrum disorder (ASD). Boolean models have assisted in understanding biological system dynamics since Kauffman's 1969 discovery, and various analytical tools for regulatory networks have been developed. This study examined the protein-protein interaction network created in our previous publication of four ASD patients using the SPIDDOR R package, a Boolean model-based method. The aim is to examine how patients' genetic variations in INTS6L, USP9X, RSK4, FGF5, FLNA, SUMF1, and IDS affect mTOR and Wnt cell signaling convergence. The Boolean network analysis revealed abnormal activation levels of essential proteins such as β-catenin, MTORC1, RPS6, eIF4E, Cadherin, and SMAD. These proteins affect gene expression, translation, cell adhesion, shape, and migration. Patients 1 and 2 showed consistent patterns of increased β-catenin activity and decreased MTORC1, RPS6, and eIF4E activity. However, patient 2 had an independent decrease in Cadherin and SMAD activity due to the FLNA mutation. Patients 3 and 4 have an abnormal activation of the mTOR pathway, which includes the MTORC1, RPS6, and eIF4E genes. The shared mTOR pathway behavior in these patients is explained by a shared mutation in two closely related proteins (SUMF1 and IDS). Diverse activities in β-catenin, MTORC1, RPS6, eIF4E, Cadherin, and SMAD contributed to the reported phenotype in these individuals. Furthermore, it unveiled the potential therapeutic options that could be suggested to these individuals.

Exome-wide association study identifies KDELR3 mutations in extreme myopia

Extreme myopia (EM), defined as a spherical equivalent (SE) ≤ -10.00 diopters (D), is one of the leading causes of sight impairment. Known EM-associated variants only explain limited risk and are inadequate for clinical decision-making. To discover risk genes, we performed a whole-exome sequencing (WES) on 449 EM individuals and 9606 controls. We find a significant excess of rare protein-truncating variants (PTVs) in EM cases, enriched in the retrograde vesicle-mediated transport pathway. Employing single-cell RNA-sequencing (scRNA-seq) and a single-cell polygenic burden score (scPBS), we pinpointed PI16 + /SFRP4+ fibroblasts as the most relevant cell type. We observed that KDELR3 is highly expressed in scleral fibroblast and involved in scleral extracellular matrix (ECM) organization. The zebrafish model revealed that kdelr3 downregulation leads to elongated ocular axial length and increased lens diameter. Together, our study provides insight into the genetics of EM in humans and highlights KDELR3's role in EM pathogenesis.

Identification of Schizophrenia Susceptibility Loci in the Urban Taiwanese Population

Background and Objectives: Genomic studies have identified several SNP loci associated with schizophrenia in East Asian populations. Environmental factors, particularly urbanization, play a significant role in schizophrenia development. This study aimed to identify schizophrenia susceptibility loci and characterize their biological functions and molecular pathways in Taiwanese urban Han individuals. Materials and Methods: Participants with schizophrenia were recruited from the Taiwan Precision Medicine Initiative at Tri-Service General Hospital. Genotype-phenotype association analysis was performed, with significant variants annotated and analyzed for functional relevance. Results: A total of 137 schizophrenia patients and 26,129 controls were enrolled. Ten significant variants (p < 1 × 10-5) and 15 expressed genes were identified, including rs1010840 (SOWAHC and RGPD6), rs11083963 (TRPM4), rs11619878 (LINC00355 and LINC01052), rs117010638 (AGBL1 and MIR548AP), rs1170702 (LINC01680 and LINC01720), rs12028521 (KAZN and PRDM2), rs12859097 (DMD), rs1556812 (ATP11A), rs78144262 (LINC00977), and rs9997349 (ENPEP). These variants and associated genes are involved in immune response, blood pressure regulation, muscle function, and the cytoskeleton. Conclusions: Identified variants and associated genes suggest a potential genetic predisposition to schizophrenia in the Taiwanese urban Han population, highlighting the importance of potential comorbidities, considering population-specific genetic and environmental interactions.

Selective advantage of mutant stem cells in human clonal hematopoiesis is associated with attenuated response to inflammation and aging

Clonal hematopoiesis (CH) arises when hematopoietic stem cells (HSCs) acquire mutations, most frequently in the DNMT3A and TET2 genes, conferring a competitive advantage through mechanisms that remain unclear. To gain insight into how CH mutations enable gradual clonal expansion, we used single-cell multi-omics with high-fidelity genotyping on human CH bone marrow (BM) samples. Most of the selective advantage of mutant cells occurs within HSCs. DNMT3A- and TET2-mutant clones expand further in early progenitors, while TET2 mutations accelerate myeloid maturation in a dose-dependent manner. Unexpectedly, both mutant and non-mutant HSCs from CH samples are enriched for inflammatory and aging transcriptomic signatures, compared with HSCs from non-CH samples, revealing a non-cell-autonomous effect. However, DNMT3A- and TET2-mutant HSCs have an attenuated inflammatory response relative to wild-type HSCs within the same sample. Our data support a model whereby CH clones are gradually selected because they are resistant to the deleterious impact of inflammation and aging.

Germline variants in patients diagnosed with pediatric soft tissue sarcoma

BACKGROUND

While soft tissue sarcomas affect younger patients, few studies have assessed the distribution of underlying pathogenic germline variants.

PATIENTS AND METHODS

We retrospectively identified all pediatric and young adult patients (0-22 years) at Haukeland University Hospital, Norway (1981-2019), through clinical and pathological records. We identified n = 46 eligible patients. From these 46 patients, adequate material representing normal tissue was available for n = 41 cases (n = 24 diagnosed with rhabdomyosarcoma, 9 with synovial sarcomas, 2 with Ewing sarcomas, and 6 without further classification), with matching tumor tissue for n = 40. Normal tissue samples were analyzed for germline pathogenic variants (PVs) by targeted sequencing of 360 cancer genes.

RESULTS

Out of the 41 analyzed cases, we found PVs or likely PVs in 7 (17%). These variants were found in TP53, MUTYH, FANCC, DICER1, FANCA, MYO3A, and MYO5B. Supporting the causality of these PVs, four cases revealed loss of heterozygosity (LOH) of the wild-type allele in the tumor tissue, one patient with a PV in DICER1 had a second somatic variant in DICER1, and a patient with a PV in TP53 had the altered allele amplified in the tumor. For three out of five with available family history, a history of other cancers in relatives was recorded. Among genes with variants of uncertain significance, CHD1L was of particular interest, revealing a stop-gain and a missense variant.

INTERPRETATION

A high fraction of young patients with soft tissue sarcoma harbor PVs. Among the genes affected, we substantiate a potential role of MYO5B and propose a potential role for MYO3A.

Whole-exome sequencing identifies high-confidence genes for tic disorders in a Chinese Han population

BACKGROUND

Tic disorder (TD) is a polygenic neurodevelopmental disorder with high susceptibility. However, identifying high-confidence risk genes has been challenging due to poor replication across multiple studies.

METHODS

Whole-exome sequencing was performed on 390 TD patients and 372 unaffected individuals in a Chinese Han population. Analysis of variance, burden analysis and in silico prediction were used to identify candidate genes for TD. To facilitate data analysis and to focus on high-confidence genes, we defined a panel of 160 genes as known causal or candidate TD genes from previous studies. Gene enrichment and protein-protein interaction analysis were utilized to detect potential novel TD risk genes.

RESULTS

Totally, 14 variants across 12 known TD candidate genes were considered potential susceptibility variants. Ten variants across 10 known TD candidate genes were identified as potential disease-causing variants. Burden analysis identified variants of 28 known genes were significantly excess in TD patients. In addition, 354 previously unproven TD genes are over-represented in patients. Genes enriched in the PI3K-Akt signaling, sphingolipid metabolism and serotonergic synaptic pathways, as well as those interacting with FN1, were considered potential new candidate genes for TD.

CONCLUSIONS

This is the largest WES study focusing on TD patients in a Chinese Han population. Several variants recurring in our cohort were identified as high-confidence risk loci for TD. Moreover, we provided potential new risk genes that may be prioritized for further investigation.

Immune dysregulation due to bi-allelic mutation of the actin remodeling protein DIAPH1

Children with severe inflammatory diseases are challenging to diagnose and treat, and the etiology of disease often remains unexplained. Here we present DIAPH1 deficiency as an unexpected genetic finding in a child with fatal inflammatory bowel disease who also displayed complex neurological and developmental phenotypes. Bi-allelic mutations of DIAPH1 were first described in patients with a severe neurological phenotype including microcephaly, intellectual disability, seizures, and blindness. Recent findings have expanded the clinical phenotype of DIAPH1 deficiency to include severe susceptibility to infections, placing this monogenic disease amongst the etiologies of inborn errors of immunity. Immune phenotypes in DIAPH1 deficiency are largely driven aberrant lymphocyte activation, particularly the failure to form an effective immune synapse in T cells. We present the case of a child with a novel homozygous deletion in DIAPH1, leading to a premature truncation in the Lasso domain of the protein. Unlike other cases of DIAPH1 deficiency, this patient did not have seizures or lung infections. Her major immune-related clinical symptoms were inflammation and enteropathy, diarrhea and failure to thrive. This patient did not show T or B cell lymphopenia but did have dramatically reduced naïve CD4+ and CD8+ T cells, expanded CD4-CD8- T cells, and elevated IgE. Similar to other cases of DIAPH1 deficiency, this patient had non-hematological phenotypes including microcephaly, developmental delay, and impaired vision. This patient's symptSoms of immune dysregulation were not successfully controlled and were ultimately fatal. This case expands the clinical spectrum of DIAPH1 deficiency and reveals that autoimmune or inflammatory enteropathy may be the most prominent immunological manifestation of disease.

Novel drug resistance mechanisms and drug targets in BRAF-mutated peritoneal metastasis from colorectal cancer

BACKGROUND

Patients with peritoneal metastasis from colorectal cancer (PM-CRC) have inferior prognosis and respond particularly poorly to chemotherapy. This study aims to identify the molecular explanation for the observed clinical behavior and suggest novel treatment strategies in PM-CRC.

METHODS

Tumor samples (230) from a Norwegian national cohort undergoing surgery and hyperthermic intraperitoneal chemotherapy (HIPEC) with mitomycin C (MMC) for PM-CRC were subjected to targeted DNA sequencing, and associations with clinical data were analyzed. mRNA sequencing was conducted on a subset of 30 samples to compare gene expression in tumors harboring BRAF or KRAS mutations and wild-type tumors.

RESULTS

BRAF mutations were detected in 27% of the patients, and the BRAF-mutated subgroup had inferior overall survival compared to wild-type cases (median 16 vs 36 months, respectively, p < 0.001). BRAF mutations were associated with RNF43/RSPO aberrations and low expression of negative Wnt regulators (ligand-dependent Wnt activation). Furthermore, BRAF mutations were associated with gene expression changes in transport solute carrier proteins (specifically SLC7A6) and drug metabolism enzymes (CES1 and CYP3A4) that could influence the efficacy of MMC and irinotecan, respectively. BRAF-mutated tumors additionally exhibited increased expression of members of the novel butyrophilin subfamily of immune checkpoint molecules (BTN1A1 and BTNL9).

CONCLUSIONS

BRAF mutations were frequently detected and were associated with particularly poor survival in this cohort, possibly related to ligand-dependent Wnt activation and altered drug transport and metabolism that could confer resistance to MMC and irinotecan. Drugs that target ligand-dependent Wnt activation or the BTN immune checkpoints could represent two novel therapy approaches.

Variant Impact Predictor database (VIPdb), version 2: Trends from 25 years of genetic variant impact predictors

Background

Variant interpretation is essential for identifying patients' disease-causing genetic variants amongst the millions detected in their genomes. Hundreds of Variant Impact Predictors (VIPs), also known as Variant Effect Predictors (VEPs), have been developed for this purpose, with a variety of methodologies and goals. To facilitate the exploration of available VIP options, we have created the Variant Impact Predictor database (VIPdb).

Results

The Variant Impact Predictor database (VIPdb) version 2 presents a collection of VIPs developed over the past 25 years, summarizing their characteristics, ClinGen calibrated scores, CAGI assessment results, publication details, access information, and citation patterns. We previously summarized 217 VIPs and their features in VIPdb in 2019. Building upon this foundation, we identified and categorized an additional 186 VIPs, resulting in a total of 403 VIPs in VIPdb version 2. The majority of the VIPs have the capacity to predict the impacts of single nucleotide variants and nonsynonymous variants. More VIPs tailored to predict the impacts of insertions and deletions have been developed since the 2010s. In contrast, relatively few VIPs are dedicated to the prediction of splicing, structural, synonymous, and regulatory variants. The increasing rate of citations to VIPs reflects the ongoing growth in their use, and the evolving trends in citations reveal development in the field and individual methods.

Conclusions

VIPdb version 2 summarizes 403 VIPs and their features, potentially facilitating VIP exploration for various variant interpretation applications.

Availability

VIPdb version 2 is available at https://genomeinterpretation.org/vipdb.

Genetic Characterization of 191 Probands with Inherited Retinal Dystrophy by Targeted NGS Analysis

Inherited retinal diseases (IRDs) represent a frequent cause of blindness in children and adults. As a consequence of the phenotype and genotype heterogeneity of the disease, it is difficult to have a specific diagnosis without molecular testing. To date, over 340 genes and loci have been associated with IRDs. We present the molecular finding of 191 individuals with IRD, analyzed by targeted next-generation sequencing (NGS). For 67 of them, we performed a family segregation study, considering a total of 126 relatives. A total of 359 variants were identified, 44 of which were novel. Genetic diagnostic yield was 41%. However, after stratifying the patients according to their clinical suspicion, diagnostic yield was higher for well-characterized diseases such as Stargardt disease (STGD), at 65%, and for congenital stationary night blindness 2 (CSNB2), at 64%. Diagnostic yield was higher in the patient group where family segregation analysis was possible (68%) and it was higher in younger (55%) than in older patients (33%). The results of this analysis demonstrated that targeted NGS is an effective method for establishing a molecular genetic diagnosis of IRDs. Furthermore, this study underlines the importance of segregation studies to understand the role of genetic variants with unknow pathogenic role.

Comprehensive variant analysis of phospholipase A2 superfamily genes in large Chinese Parkinson' s disease cohorts

To clarify the genetic role of phospholipase A2 (PLA2) genes in Parkinson's disease (PD), we performed a genetic association study in large Chinese population cohorts using next-generation sequencing. In this study, we analyzed both rare and common variants of 38 phospholipase A2 genes in two large cohorts. We detected 1558 and 1115 rare variants in these two cohorts, respectively. In both cohorts, we observed suggestive associations between specific subgroups and the risk of PD. At the single-gene level, several genes (PLA2G2D, PLA2G12A, PLA2G12B, PLA2G4F, PNPLA1, PNPLA3, PNPLA7, PLA2G7, PLA2G15, PLAAT5, and ABHD12) are suggestively associated with PD. Meanwhile, 364 and 2261 common variants were identified in two cohorts, respectively. Our study has expanded the genetic spectrum of the PLA2 family genes and suggested potential pathogenetic roles of PLA2 superfamily in PD.

Interplay Between Plasma Glycine and Branched-Chain Amino Acids Contributes to the Development of Hypertension and Coronary Heart Disease

BACKGROUND

Higher levels of plasma glycine are linked to a reduced risk, while increased levels of total branched-chain amino acids (tBCAAs) are associated with a higher risk of essential hypertension and coronary heart disease (CHD). As these metabolic components are interconnected, analyzing the tBCAAs/glycine ratio may help to understand their interplay in the pathogenesis of cardiovascular disease.

METHODS

The Cox regression approach was combined with the development of novel genetic tools for assessments of associations between plasma metabolomic data (glycine, tBCAAs, and tBCAAs/glycine ratio) from the UK Biobank and the development of hypertension and CHD. Genome-wide association study was performed on 186 523 White UK Biobank participants to identify new independent genetic instruments for the 2-sample Mendelian randomization analyses. P-gain statistic >10 identified instruments associated with tBCAAs/glycine ratio significantly stronger compared with individual amino acids. Outcomes of genome-wide association study on hypertension and CHD were derived from the UK Biobank (nonoverlapping sample), FinnGen, and CARDIoGRAMplusC4D.

RESULTS

The tBCAAs/glycine ratio was prospectively associated with a higher risk of developing hypertension and CHD (hazard ratio quintile Q5 versus Q1, 1.196 [95% CI, 1.109-1.289] and 1.226 [95% CI, 1.160-1.296], respectively). Mendelian randomization analysis demonstrated that tBCAAs/glycine ratio (P-gain >10) was a risk factor for hypertension (meta-analyzed inverse-variance weighted causal estimate 0.45 log odds ratio/SD (95% CI, 0.26-0.64) and CHD (0.48 [95% CI, 0.29-0.67]) with an absolute effect significantly larger compared with the effect of glycine (-0.06 [95% CI, -0.1 to -0.03] and -0.08 [95% CI, -0.11 to -0.05], respectively) or tBCAAs (0.22 [95% CI, 0.09-0.34] and 0.12 [95% CI, 0.01-0.24], respectively).

CONCLUSIONS

The total BCAAs/glycine ratio is a key element of the metabolic signature contributing to hypertension and CHD, which may reflect biological pathways shared by glycine and tBCAAs.

CytoSIP: an annotated structural atlas for interactions involving cytokines or cytokine receptors

Therapeutic agents targeting cytokine-cytokine receptor (CK-CKR) interactions lead to the disruption in cellular signaling and are effective in treating many diseases including tumors. However, a lack of universal and quick access to annotated structural surface regions on CK/CKR has limited the progress of a structure-driven approach in developing targeted macromolecular drugs and precision medicine therapeutics. Herein we develop CytoSIP (Single nucleotide polymorphisms (SNPs), Interface, and Phenotype), a rich internet application based on a database of atomic interactions around hotspots in experimentally determined CK/CKR structural complexes. CytoSIP contains: (1) SNPs on CK/CKR; (2) interactions involving CK/CKR domains, including CK/CKR interfaces, oligomeric interfaces, epitopes, or other drug targeting surfaces; and (3) diseases and phenotypes associated with CK/CKR or SNPs. The database framework introduces a unique tri-level SIP data model to bridge genetic variants (atomic level) to disease phenotypes (organism level) using protein structure (complexes) as an underlying framework (molecule level). Customized screening tools are implemented to retrieve relevant CK/CKR subset, which reduces the time and resources needed to interrogate large datasets involving CK/CKR surface hotspots and associated pathologies. CytoSIP portal is publicly accessible at https://CytoSIP.biocloud.top , facilitating the panoramic investigation of the context-dependent crosstalk between CK/CKR and the development of targeted therapeutic agents.

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

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

Comparative analysis of in-silico tools in identifying pathogenic variants in dominant inherited retinal diseases

Inherited retinal diseases (IRDs) are a group of rare genetic eye conditions that cause blindness. Despite progress in identifying genes associated with IRDs, improvements are necessary for classifying rare autosomal dominant (AD) disorders. AD diseases are highly heterogenous, with causal variants being restricted to specific amino acid changes within certain protein domains, making AD conditions difficult to classify. Here, we aim to determine the top-performing in-silico tools for predicting the pathogenicity of AD IRD variants. We annotated variants from ClinVar and benchmarked 39 variant classifier tools on IRD genes, split by inheritance pattern. Using area-under-the-curve (AUC) analysis, we determined the top-performing tools and defined thresholds for variant pathogenicity. Top-performing tools were assessed using genome sequencing on a cohort of participants with IRDs of unknown etiology. MutScore achieved the highest accuracy within AD genes, yielding an AUC of 0.969. When filtering for AD gain-of-function and dominant negative variants, BayesDel had the highest accuracy with an AUC of 0.997. Five participants with variants in NR2E3, RHO, GUCA1A, and GUCY2D were confirmed to have dominantly inherited disease based on pedigree, phenotype, and segregation analysis. We identified two uncharacterized variants in GUCA1A (c.428T>A, p.Ile143Thr) and RHO (c.631C>G, p.His211Asp) in three participants. Our findings support using a multi-classifier approach comprised of new missense classifier tools to identify pathogenic variants in participants with AD IRDs. Our results provide a foundation for improved genetic diagnosis for people with IRDs.

Regulatory elements in SEM1-DLX5-DLX6 (7q21.3) locus contribute to genetic control of coronal nonsyndromic craniosynostosis and bone density-related traits

Purpose

The etiopathogenesis of coronal nonsyndromic craniosynostosis (cNCS), a congenital condition defined by premature fusion of 1 or both coronal sutures, remains largely unknown.

Methods

We conducted the largest genome-wide association study of cNCS followed by replication, fine mapping, and functional validation of the most significant region using zebrafish animal model.

Results

Genome-wide association study identified 6 independent genome-wide-significant risk alleles, 4 on chromosome 7q21.3 SEM1-DLX5-DLX6 locus, and their combination conferred over 7-fold increased risk of cNCS. The top variants were replicated in an independent cohort and showed pleiotropic effects on brain and facial morphology and bone mineral density. Fine mapping of 7q21.3 identified a craniofacial transcriptional enhancer (eDlx36) within the linkage region of the top variant (rs4727341; odds ratio [95% confidence interval], 0.48[0.39-0.59]; P = 1.2E-12) that was located in SEM1 intron and enriched in 4 rare risk variants. In zebrafish, the activity of the transfected human eDlx36 enhancer was observed in the frontonasal prominence and calvaria during skull development and was reduced when the 4 rare risk variants were introduced into the sequence.

Conclusion

Our findings support a polygenic nature of cNCS risk and functional role of craniofacial enhancers in cNCS susceptibility with potential broader implications for bone health.

Temporal landscape and translational regulation of A-to-I RNA editing in mouse retina development

BACKGROUND

The significance of A-to-I RNA editing in nervous system development is widely recognized; however, its influence on retina development remains to be thoroughly understood.

RESULTS

In this study, we performed RNA sequencing and ribosome profiling experiments on developing mouse retinas to characterize the temporal landscape of A-to-I editing. Our findings revealed temporal changes in A-to-I editing, with distinct editing patterns observed across different developmental stages. Further analysis showed the interplay between A-to-I editing and alternative splicing, with A-to-I editing influencing splicing efficiency and the quantity of splicing events. A-to-I editing held the potential to enhance translation diversity, but this came at the expense of reduced translational efficiency. When coupled with splicing, it could produce a coordinated effect on gene translation.

CONCLUSIONS

Overall, this study presents a temporally resolved atlas of A-to-I editing, connecting its changes with the impact on alternative splicing and gene translation in retina development.

Negative Selection Allows DNA Mismatch Repair-Deficient Mouse Fibroblasts In Vitro to Tolerate High Levels of Somatic Mutations

Substantial numbers of somatic mutations have been found to accumulate with age in different human tissues. Clonal cellular amplification of some of these mutations can cause cancer and other diseases. However, it is as yet unclear if and to what extent an increased burden of random mutations can affect cellular function without clonal amplification. We tested this in cell culture, which avoids the limitation that an increased mutation burden in vivo typically leads to cancer. We performed single-cell whole-genome sequencing of primary fibroblasts from DNA mismatch repair (MMR) deficient Msh2-/- mice and littermate control animals after long-term passaging. Apart from analyzing somatic mutation burden we analyzed clonality, mutational signatures, and hotspots in the genome, characterizing the complete landscape of somatic mutagenesis in normal and MMR-deficient mouse primary fibroblasts during passaging. While growth rate of Msh2-/- fibroblasts was not significantly different from the controls, the number of de novo single-nucleotide variants (SNVs) increased linearly up until at least 30,000 SNVs per cell, with the frequency of small insertions and deletions (INDELs) plateauing in the Msh2-/- fibroblasts to about 10,000 INDELS per cell. We provide evidence for negative selection and large-scale mutation-driven population changes, including significant clonal expansion of preexisting mutations and widespread cell-strain-specific hotspots. Overall, our results provide evidence that increased somatic mutation burden drives significant cell evolutionary changes in a dynamic cell culture system without significant effects on growth. Since similar selection processes against mutations preventing organ and tissue dysfunction during aging are difficult to envision, these results suggest that increased somatic mutation burden can play a causal role in aging and diseases other than cancer.

Low-frequency and rare genetic variants associated with rheumatoid arthritis risk

Rheumatoid arthritis (RA) has an estimated heritability of nearly 50%, which is particularly high in seropositive RA. HLA alleles account for a large proportion of this heritability, in addition to many common single-nucleotide polymorphisms with smaller individual effects. Low-frequency and rare variants, such as those captured by next-generation sequencing, can also have a large role in heritability in some individuals. Rare variant discovery has informed the development of drugs such as inhibitors of PCSK9 and Janus kinases. Some 34 low-frequency and rare variants are currently associated with RA risk. One variant (19:10352442G>C in TYK2) was identified in five separate studies, and might therefore represent a promising therapeutic target. Following a set of best practices in future studies, including studying diverse populations, using large sample sizes, validating RA and serostatus, replicating findings, adjusting for other variants and performing functional assessment, could help to ensure the relevance of identified variants. Exciting opportunities are now on the horizon for genetics in RA, including larger datasets and consortia, whole-genome sequencing and direct applications of findings in the management, and especially treatment, of RA.

Pan-Genome Analysis of Wolbachia, Endosymbiont of Diaphorina citri, Reveals Independent Origin in Asia and North America

Wolbachia, a group of Gram-negative symbiotic bacteria, infects nematodes and a wide range of arthropods. Diaphorina citri Kuwayama, the vector of Candidatus Liberibacter asiaticus (CLas) that causes citrus greening disease, is naturally infected with Wolbachia (wDi). However, the interaction between wDi and D. citri remains poorly understood. In this study, we performed a pan-genome analysis using 65 wDi genomes to gain a comprehensive understanding of wDi. Based on average nucleotide identity (ANI) analysis, we classified the wDi strains into Asia and North America strains. The ANI analysis, principal coordinates analysis (PCoA), and phylogenetic tree analysis supported that the D. citri in Florida did not originate from China. Furthermore, we found that a significant number of core genes were associated with metabolic pathways. Pathways such as thiamine metabolism, type I secretion system, biotin transport, and phospholipid transport were highly conserved across all analyzed wDi genomes. The variation analysis between Asia and North America wDi showed that there were 39,625 single-nucleotide polymorphisms (SNPs), 2153 indels, 10 inversions, 29 translocations, 65 duplications, 10 SV-based insertions, and 4 SV-based deletions. The SV-based insertions and deletions involved genes encoding transposase, phage tail tube protein, ankyrin repeat (ANK) protein, and group II intron-encoded protein. Pan-genome analysis of wDi contributes to our understanding of the geographical population of wDi, the origin of hosts of D. citri, and the interaction between wDi and its host, thus facilitating the development of strategies to control the insects and huanglongbing (HLB).

Fine-mapping analysis including over 254,000 East Asian and European descendants identifies 136 putative colorectal cancer susceptibility genes

Genome-wide association studies (GWAS) have identified more than 200 common genetic variants independently associated with colorectal cancer (CRC) risk, but the causal variants and target genes are mostly unknown. We sought to fine-map all known CRC risk loci using GWAS data from 100,204 cases and 154,587 controls of East Asian and European ancestry. Our stepwise conditional analyses revealed 238 independent association signals of CRC risk, each with a set of credible causal variants (CCVs), of which 28 signals had a single CCV. Our cis-eQTL/mQTL and colocalization analyses using colorectal tissue-specific transcriptome and methylome data separately from 1299 and 321 individuals, along with functional genomic investigation, uncovered 136 putative CRC susceptibility genes, including 56 genes not previously reported. Analyses of single-cell RNA-seq data from colorectal tissues revealed 17 putative CRC susceptibility genes with distinct expression patterns in specific cell types. Analyses of whole exome sequencing data provided additional support for several target genes identified in this study as CRC susceptibility genes. Enrichment analyses of the 136 genes uncover pathways not previously linked to CRC risk. Our study substantially expanded association signals for CRC and provided additional insight into the biological mechanisms underlying CRC development.

Molecular characterization of the evolution of premalignant lesions in the upper aerodigestive tract

Introduction

Early relapse and development of metastatic disease are some of the primary reasons for the poor prognosis of patients with head and neck squamous cell carcinoma (HNSCC). HNSCC is a heterogeneous disease which may develop in large premalignant fields of genetically altered cells. Yet knowing which individuals will progress and develop clinically significant cancers during their lifetimes remains one of the most important challenges of reducing HNSCC morbidity and mortality. To further elucidate the molecular mechanisms, we performed a focused analysis of the genome and immune microenvironment from multiple, matched normal squamous tissue, premalignant lesions, as well as primary and recurrent tumors from seven patients with p16-negative HNSCC.

Methods

We performed targeted panel Next Generation Sequencing (161 genes) to analyze somatic variants from sequentially collected, matched formalin-fixed paraffin-embedded tissue (normal, premalignant, HNSCC) from two patients. These samples plus samples from five additional patients were analyzed with the Nanostring PanCancer Immune Panel. In addition, we performed shallow whole genome sequencing (0.5x coverage on average) on samples from three of these patients. Patients were, apart from one case, primarily treated with curative-intent surgery, and received subsequent adjuvant treatment, if indicated.

Results

The most frequently mutated genes were TP53 and NOTCH1. Other mutated genes included NOTCH3 and CDKN2A, among others. A significant number of mutations were private to dysplasia and invasive carcinoma, respectively, however, almost 20% were shared between them. Increasing genomic instability was observed when comparing histologically normal squamous mucosa with higher levels of dysplasia. High-grade dysplasia showed similarly rearranged genomes as invasive carcinoma. Pathways related to interferon alpha and gamma response were upregulated even in moderate dysplastic lesions with increasing expression in higher grades of dysplasia and carcinoma. SPINK5, a known tumor suppressor gene in HNSCC, was already downregulated in low-grade dysplastic lesions, indicating an early deactivation in the evolution of the disease.

Conclusion

Genomic alterations as well as aberrant immune gene expression can be observed early in the evolution of tumors of the upper aerodigestive tract, highlighting the potential for targeting early mechanisms of disease progression.

Predisposing deleterious variants in the cancer-associated human kinases in the global populations

Human kinases play essential and diverse roles in the cellular activities of maintaining homeostasis and growth. Genetic mutations in the genes encoding the kinases (or phosphotransferases) have been linked with various types of cancers. In this study, we cataloged mutations in 500 kinases genes in >65,000 individuals of global populations from the Human Genetic Diversity Project (HGDP) and ExAC databases, and assessed their potentially deleterious impact by using the in silico tools SIFT, Polyphen2, and CADD. The analysis highlighted 35 deleterious non-synonymous SNVs in the ExAC and 5 SNVs in the HGDP project. Notably, a higher number of deleterious mutations was observed in the Non-Finnish Europeans (26 SNVs), followed by the Africans (14 SNVs), East Asians (13 SNVs), and South Asians (12 SNVs). The gene set enrichment analysis highlighted NTRK1 and FGFR3 being most significantly enriched among the kinases. The gene expression analysis revealed over-expression of NTRK1 in liver cancer, whereas, FGFR3 was found over-expressed in lung, breast, and liver cancers compared to their expression in the respective normal tissues. Also, 13 potential drugs were identified that target the NTRK1 protein, whereas 6 potential drugs for the FGFR3 target were identified. Taken together, the study provides a framework for exploring the predisposing germline mutations in kinases to suggest the underlying pathogenic mechanisms in cancers. The potential drugs are also suggested for personalized cancer management.

TMBcalc: a computational pipeline for identifying pan-cancer Tumor Mutational Burden gene signatures

Background

In the precision medicine era, identifying predictive factors to select patients most likely to benefit from treatment with immunological agents is a crucial and open challenge in oncology.

Methods

This paper presents a pan-cancer analysis of Tumor Mutational Burden (TMB). We developed a novel computational pipeline, TMBcalc, to calculate the TMB. Our methodology can identify small and reliable gene signatures to estimate TMB from custom targeted-sequencing panels. For this purpose, our pipeline has been trained on top of 17 cancer types data obtained from TCGA.

Results

Our results show that TMB, computed through the identified signature, strongly correlates with TMB obtained from whole-exome sequencing (WES).

Conclusion

We have rigorously analyzed the effectiveness of our methodology on top of several independent datasets. In particular we conducted a comprehensive testing on: (i) 126 samples sourced from the TCGA database; few independent whole-exome sequencing (WES) datasets linked to colon, breast, and liver cancers, all acquired from the EGA and the ICGC Data Portal. This rigorous evaluation clearly highlights the robustness and practicality of our approach, positioning it as a promising avenue for driving substantial progress within the realm of clinical practice.

The genetic basis of early-onset hereditary ataxia in Iran: results of a national registry of a heterogeneous population

BACKGROUND

To investigate the genetics of early-onset progressive cerebellar ataxia in Iran, we conducted a study at the Children's Medical Center (CMC), the primary referral center for pediatric disorders in the country, over a three-year period from 2019 to 2022. In this report, we provide the initial findings from the national registry.

METHODS

We selected all early-onset patients with an autosomal recessive mode of inheritance to assess their phenotype, paraclinical tests, and genotypes. The clinical data encompassed clinical features, the Scale for the Assessment and Rating of Ataxia (SARA) scores, Magnetic Resonance Imaging (MRI) results, Electrodiagnostic exams (EDX), and biomarker features. Our genetic investigations included single-gene testing, Whole Exome Sequencing (WES), and Whole Genome Sequencing (WGS).

RESULTS

Our study enrolled 162 patients from various geographic regions of our country. Among our subpopulations, we identified known and novel pathogenic variants in 42 genes in 97 families. The overall genetic diagnostic rate was 59.9%. Notably, we observed PLA2G6, ATM, SACS, and SCA variants in 19, 14, 12, and 10 families, respectively. Remarkably, more than 59% of the cases were attributed to pathogenic variants in these genes.

CONCLUSIONS

Iran, being at the crossroad of the Middle East, exhibits a highly diverse genetic etiology for autosomal recessive hereditary ataxia. In light of this heterogeneity, the development of preventive strategies and targeted molecular therapeutics becomes crucial. A national guideline for the diagnosis and management of patients with these conditions could significantly aid in advancing healthcare approaches and improving patient outcomes.

Exome sequencing identified mutations in the WNT1 and COL1A2 genes in osteogenesis imperfecta cases

BACKGROUND

Osteogenesis imperfecta (OI) is a heritable connective tissue disorder characterized by bone deformities, fractures and reduced bone mass. OI can be inherited as a dominant, recessive, or X-linked disorder. The mutational spectrum has shown that autosomal dominant mutations in the type I collagen-encoding genes are responsible for OI in 85% of the cases. Apart from collagen genes, mutations in more than 20 other genes, such as CRTAP, CREB3L1, MBTPS2, P4HB, SEC24D, SPARC, FKBP10, LEPRE1, PLOD2, PPIB, SERPINF1, SERPINH1, SP7, WNT1, BMP1, TMEM38B, and IFITM5 have been reported in OI.

METHODS AND RESULTS

To understand the genetic cause of OI in four cases, we conducted whole exome sequencing, followed by Sanger sequencing. In case #1, we identified a novel c.506delG homozygous mutation in the WNT1 gene, resulting in a frameshift and early truncation of the protein at the 197th amino acid. In cases #2, 3 and 4, we identified a heterozygous c.838G > A mutation in the COL1A2 gene, resulting in a p.Gly280Ser substitution. The clinvar frequency of this mutation is 0.000008 (GnomAD-exomes). This mutation has been identified by other studies as well and appears to be a mutational hot spot. These pathogenic mutations were found to be absent in 96 control samples analyzed for these sites. The presence of these mutations in the cases, their absence in controls, their absence or very low frequency in general population, and their evaluation using various in silico prediction tools suggested their pathogenic nature.

CONCLUSIONS

Mutations in the WNT1 and COL1A2 genes explain these cases of osteogenesis imperfecta.

Statistical methods for assessing the effects of de novo variants on birth defects

With the development of next-generation sequencing technology, de novo variants (DNVs) with deleterious effects can be identified and investigated for their effects on birth defects such as congenital heart disease (CHD). However, statistical power is still limited for such studies because of the small sample size due to the high cost of recruiting and sequencing samples and the low occurrence of DNVs. DNV analysis is further complicated by genetic heterogeneity across diseased individuals. Therefore, it is critical to jointly analyze DNVs with other types of genomic/biological information to improve statistical power to identify genes associated with birth defects. In this review, we discuss the general workflow, recent developments in statistical methods, and future directions for DNV analysis.

Biallelic NAA60 variants with impaired n-terminal acetylation capacity cause autosomal recessive primary familial brain calcifications

Primary familial brain calcification (PFBC) is characterized by calcium deposition in the brain, causing progressive movement disorders, psychiatric symptoms, and cognitive decline. PFBC is a heterogeneous disorder currently linked to variants in six different genes, but most patients remain genetically undiagnosed. Here, we identify biallelic NAA60 variants in ten individuals from seven families with autosomal recessive PFBC. The NAA60 variants lead to loss-of-function with lack of protein N-terminal (Nt)-acetylation activity. We show that the phosphate importer SLC20A2 is a substrate of NAA60 in vitro. In cells, loss of NAA60 caused reduced surface levels of SLC20A2 and a reduction in extracellular phosphate uptake. This study establishes NAA60 as a causal gene for PFBC, provides a possible biochemical explanation of its disease-causing mechanisms and underscores NAA60-mediated Nt-acetylation of transmembrane proteins as a fundamental process for healthy neurobiological functioning.

Whole-Exome Sequencing (WES) Reveals Novel Sex-Specific Gene Variants in Non-Alcoholic Steatohepatitis (MASH)

Non-alcoholic steatohepatitis (NASH, also known as MASH) is a severe form of non-alcoholic fatty liver disease (NAFLD, also known as MASLD). Emerging data indicate that the progression of the disease to MASH is higher in postmenopausal women and that genetic susceptibility increases the risk of MASH-related cirrhosis. This study aimed to investigate the association between genetic polymorphisms in MASH and sexual dimorphism. We applied whole-exome sequencing (WES) to identify gene variants in 8 age-adjusted matched pairs of livers from both male and female patients. Sequencing alignment, variant calling, and annotation were performed using standard methods. Polymerase chain reaction (PCR) coupled with Sanger sequencing and immunoblot analysis were used to validate specific gene variants. cBioPortal and Gene Set Enrichment Analysis (GSEA) were used for actionable target analysis. We identified 148,881 gene variants, representing 57,121 and 50,150 variants in the female and male cohorts, respectively, of which 251 were highly significant and MASH sex-specific (p < 0.0286). Polymorphisms in CAPN14, SLC37A3, BAZ1A, SRP54, MYH11, ABCC1, and RNFT1 were highly expressed in male liver samples. In female samples, Polymorphisms in RGSL1, SLC17A2, HFE, NLRC5, ACTN4, SBF1, and ALPK2 were identified. A heterozygous variant 1151G>T located on 18q21.32 for ALPK2 (rs3809983) was validated by Sanger sequencing and expressed only in female samples. Immunoblot analysis confirmed that the protein level of β-catenin in female samples was 2-fold higher than normal, whereas ALPK2 expression was 0.5-fold lower than normal. No changes in the protein levels of either ALPK2 or β-catenin were observed in male samples. Our study suggests that the perturbation of canonical Wnt/β-catenin signaling observed in postmenopausal women with MASH could be the result of polymorphisms in ALPK2.

Transcriptomic, Proteomic, and Genomic Mutational Fraction Differences Based on HPV Status Observed in Patient-Derived Xenograft Models of Penile Squamous Cell Carcinoma

Metastatic penile squamous cell carcinoma (PSCC) has only a 50% response rate to first-line combination chemotherapies and there are currently no targeted-therapy approaches. Therefore, we have an urgent need in advanced-PSCC treatment to find novel therapies. Approximately half of all PSCC cases are positive for high-risk human papillomavirus (HR-HPV). Our objective was to generate HPV-positive (HPV+) and HPV-negative (HPV-) patient-derived xenograft (PDX) models and to determine the biological differences between HPV+ and HPV- disease. We generated four HPV+ and three HPV- PSCC PDX animal models by directly implanting resected patient tumor tissue into immunocompromised mice. PDX tumor tissue was found to be similar to patient tumor tissue (donor tissue) by histology and short tandem repeat fingerprinting. DNA mutations were mostly preserved in PDX tissues and similar APOBEC (apolipoprotein B mRNA editing catalytic polypeptide) mutational fractions in donor tissue and PDX tissues were noted. A higher APOBEC mutational fraction was found in HPV+ versus HPV- PDX tissues (p = 0.044), and significant transcriptomic and proteomic expression differences based on HPV status included p16 (CDKN2A), RRM2, and CDC25C. These models will allow for the direct testing of targeted therapies in PSCC and determine their response in correlation to HPV status.

Developmental outcome of electroencephalographic findings in SYNGAP1 encephalopathy

SYNGAP1 haploinsufficiency results in a developmental and epileptic encephalopathy (DEE) causing generalized epilepsies accompanied by a spectrum of neurodevelopmental symptoms. Concerning interictal epileptiform discharges (IEDs) in electroencephalograms (EEG), potential biomarkers have been postulated, including changes in background activity, fixation-off sensitivity (FOS) or eye closure sensitivity (ECS). In this study we clinically evaluate a new cohort of 36 SYNGAP1-DEE individuals. Standardized questionnaires were employed to collect clinical, electroencephalographic and genetic data. We investigated electroencephalographic findings, focusing on the cortical distribution of interictal abnormalities and their changes with age. Among the 36 SYNGAP1-DEE cases 18 presented variants in the SYNGAP1 gene that had never been previously reported. The mean age of diagnosis was 8 years and 8 months, ranging from 2 to 17 years, with 55.9% being male. All subjects had global neurodevelopmental/language delay and behavioral abnormalities; 83.3% had moderate to profound intellectual disability (ID), 91.7% displayed autistic traits, 73% experienced sleep disorders and 86.1% suffered from epileptic seizures, mainly eyelid myoclonia with absences (55.3%). A total of 63 VEEGs were revised, observing a worsening of certain EEG findings with increasing age. A disorganized background was observed in all age ranges, yet this was more common among older cases. The main IEDs were bilateral synchronous and asynchronous posterior discharges, accounting for ≥50% in all age ranges. Generalized alterations with maximum amplitude in the anterior region showed as the second most frequent IED (≥15% in all age ranges) and were also more common with increasing age. Finally, diffuse fast activity was much more prevalent in cases with 6 years or older. To the best of our knowledge, this is the first study to analyze EEG features across different age groups, revealing an increase in interictal abnormalities over infancy and adolescence. Our findings suggest that SYNGAP1 haploinsufficiency has complex effects in human brain development, some of which might unravel at different developmental stages. Furthermore, they highlight the potential of baseline EEG to identify candidate biomarkers and the importance of natural history studies to develop specialized therapies and clinical trials.

A Genome-Wide Association Study of Respiratory Syncytial Virus Infection Severity in Infants

BACKGROUND

Respiratory syncytial virus (RSV) is a significant cause of infant morbidity and mortality worldwide. Most children experience at least one 1 RSV infection by the age of two 2 years, but not all develop severe disease. However, the understanding of genetic risk factors for severe RSV is incomplete. Consequently, we conducted a genome-wide association study of RSV severity.

METHODS

Disease severity was assessed by the ReSVinet scale, in a cohort of 251 infants aged 1 week to 1 year. Genotyping data were collected from multiple European study sites as part of the RESCEU Consortium. Linear regression models were used to assess the impact of genotype on RSV severity and gene expression as measured by microarray.

RESULTS

While no SNPs reached the genome-wide statistical significance threshold (P < 5 × 10-8), we identified 816 candidate SNPs with a P-value of <1 × 10-4. Functional annotation of candidate SNPs highlighted genes relevant to neutrophil trafficking and cytoskeletal functions, including LSP1 and RAB27A. Moreover, SNPs within the RAB27A locus significantly altered gene expression (false discovery rate, FDR P < .05).

CONCLUSIONS

These findings may provide insights into genetic mechanisms driving severe RSV infection, offering biologically relevant information for future investigations.

Identification and bioinformatics analysis of a novel variant in the HERC2 gene in a patient with intellectual developmental disorder

HERC2-associated neurodevelopmental-disorders(NDD) encompass a cluster of medical conditions that arise from genetic mutations occurring within the HERC2 gene. These disorders can manifest a spectrum of symptoms that impact the brain and nervous system, including delayed psychomotor development, severe mental retardation, seizures and autistic features. Whole-Exome-Sequencing(WES) was performed on a ten-year-old male patient referred to the genetic center for genetic analysis. Blood samples were collected from the proband, his parents, and his sister to extract DNA. PCR-Sanger-sequencing was utilized to validate the findings obtained from WES. In order to obtain a more thorough understanding of the impact of the mutation, an extensive analysis was conducted using bioinformatics tools. WES data analysis identified a homozygous single nucleotide change(C > T) at position c14215 located in exon ninety-two of the HERC2 gene (NC_000015.10(NM_004667.6):c.14215C > T). The absence of this mutation among our cohort composed of four hundred normal healthy adults from the same ethnic group, and its absence in any other population database, confirms the pathogenicity of the mutation. This study revealed that the substitution of arginine with a stop codon within the Hect domain caused a premature stop codon at position 4739(p.Arg4739Ter). This mutation significantly results in the production of a truncated HERC2 protein with an incomplete HECT domain. In the final stage of ubiquitin attachment, HECT E3 ubiquitin ligases play a catalytic role by creating a thiolester intermediate using their conserved catalytic cysteine (Cys4762). This intermediate is formed before ubiquitin is transferred to a substrate protein. The truncation of the HERC2 protein is expected to disrupt its ability to perform this function, which could potentially hinder important regulatory processes related to the development and maintenance of synapses. The identification of a novel pathogenic variant, NC_000015.10(NM_004667.6):c.14215C > T, located within the ninety-two exon of the HERC2 gene, is notable for its association with an autosomal recessive inheritance pattern in cases of Intellectual Developmental Disorder(IDD). In the end, this variant could potentially play a part in the underlying mechanisms leading to the onset of intellectual developmental disorder.

CRUX, a platform for visualising, exploring and analysing cancer genome cohort data

To better understand how tumours develop, identify prognostic biomarkers and find new treatments, researchers have generated vast catalogues of cancer genome data. However, these datasets are complex, so interpreting their important features requires specialized computational skills and analytical tools, which presents a significant technical challenge. To address this, we developed CRUX, a platform for exploring genomic data from cancer cohorts. CRUX enables researchers to perform common analyses including cohort comparisons, biomarker discovery, survival analysis, and to create visualisations including oncoplots and lollipop charts. CRUX simplifies cancer genome analysis in several ways: (i) it has an easy-to-use graphical interface; (ii) it enables users to create custom cohorts, as well as analyse precompiled public and private user-created datasets; (iii) it allows analyses to be run locally to address data privacy concerns (though an online version is also available) and (iv) it makes it easy to use additional specialized tools by exporting data in the correct formats. We showcase CRUX's capabilities with case studies employing different types of cancer genome analysis, demonstrating how it can be used flexibly to generate valuable insights into cancer biology. CRUX is freely available at https://github.com/CCICB/CRUX and https://ccicb.shinyapps.io/crux (DOI: 10.5281/zenodo.8015714).

Retinoblastoma and polydactyly in a child with 46, XY, 15pstk+ karyotype-A case report and literature review

BACKGROUND

Retinoblastoma (Rb) is the most common intraocular malignancy in childhood, originating from primitive retinal stem cells or cone precursor cells. It can be triggered by mutations of the RB1 gene or amplification of the MYCN gene. Rb may rarely present with polydactyly.

METHODS

We conducted karyotype analysis, copy number variation sequencing, and whole-genome sequencing on the infant proband and his family. The clinical course and laboratory results of the proband's infant were documented and collected. We also reviewed the relevant literature.

RESULTS

A 68-day-old boy presented with preaxial polydactyly and corneal edema. His intraocular pressure (IOP) was 40/19 mmHg, and color Doppler imaging revealed vitreous solid mass-occupying lesions with calcification in the right eye. Ocular CT showed flaky high-density and calcification in the right eye. This was classified as an International Retinoblastoma Staging System group E retinoblastoma with an indication for enucleation. Enucleation and orbital implantation were performed on the child's right eye. Karyotype analysis revealed an abnormal 46, XY, 15pstk+ karyotype, and the mother exhibited diploidy of the short arm of chromosome 15. The Alx-4 development factor, 13q deletion syndrome, and the PAPA2 gene have been reported as potential mechanisms for Rb combined with polydactyly.

CONCLUSION

We report the case of a baby boy with Rb and polydactyly exhibiting a 46, XY, 15pstk+ Karyotype. We discuss potential genetic factors related to both Rb and polydactyly. Furthermore, there is a need for further exploration into the impact of chromosomal polymorphisms in Rb with polydactyly.

Whole genome sequencing followed by functional analysis of genomic deletion encompassing ERCC8 and NDUFAF2 genes in a non-consanguineous Indian family reveals dysfunctional mitochondrial bioenergetics leading to infant mortality

Genomic investigations on an infant who presented with a putative mitochondrial disorder led to identification of compound heterozygous deletion with an overlapping region of ∼142 kb encompassing two nuclear encoded genes namely ERCC8 and NDUFAF2. Investigations on fetal-derived fibroblast culture demonstrated impaired bioenergetics and mitochondrial dysfunction, which explains the phenotype and observed infant mortality in the present study. The genetic findings from this study extended the utility of whole-genome sequencing as it led to development of a MLPA-based assay for carrier screening in the extended family and the prenatal testing aiding in the birth of two healthy children.