Control-FREEC: a tool for assessing copy number and allelic content using next-generation sequencing data

Integrated proteomic and glycoproteomic analysis reveals heterogeneity and molecular signatures of brain metastases from lung adenocarcinomas

Brain metastasis is a major cause of poor prognosis and death in lung adenocarcinoma (LUAD); however, the understanding of therapeutic strategies and mechanisms for brain metastases from LUAD (BM-LUAD) remains notably limited, especially at the proteomics levels. To address this issue, we conducted integrated proteomic and glycoproteomic analyses on 49 BM-LUAD tumors, revealing two distinct subtypes of the disease: BM-S1 and BM-S2. Whole exome sequencing analysis revealed that somatic mutations in STK11 and KEAP1, as well as copy number deletions on chr19p13.3, such as STK11, UQCR11, and SLC25A23, were more frequently detected in BM-S2. In BM-S1 tumors, we observed significant infiltration of GFAP + astrocytes, as evidenced by elevated levels of GFAP, GABRA2, GABRG1 and GAP43 proteins and an enrichment of astrocytic signatures in both our proteomic data and external spatial transcriptomic data. Conversely, BM-S2 tumors demonstrated higher levels of PD-1 immune cell infiltration, supported by the upregulation of PD-1 and LAG-3 genes. These findings suggest distinct microenvironmental adaptations required by the different BM-LUAD subtypes. Additionally, we observed unique glycosylation patterns between the subtypes, with increased fucosylation in BM-S1 and enhanced sialylation in BM-S2, primarily affected by glycosylation enzymes such as FUT9, B4GALT1, and ST6GAL1. Specifically, in BM-S2, these sialylation modifications are predominantly localized to the lysosomes, underscoring the critical role of N-glycosylation in the tumor progression of BM-LUAD. Overall, our study not only provides a comprehensive multi-omic data resource but also offers valuable biological insights into BM-LUAD, highlighting potential mechanisms and therapeutic targets for further investigation.

Chromosome aberrations cause tumorigenesis through chromosomal rearrangements in a hepatocarcinogenesis rat model

Chromosome aberrations (CAs), a genotoxic potential of carcinogens, are believed to contribute to tumorigenesis by chromosomal rearrangements through micronucleus formation. However, there is no direct evidence that proves the involvement of CAs in tumorigenesis in vivo. In the current study, we sought to clarify the involvement of CAs in chemical carcinogenesis using a rat model with a pure CA-inducer hepatocarcinogen, acetamide. Whole-genome analysis indicated that hepatic tumors induced by acetamide treatment for 26-30 weeks showed a broad range of copy number alterations in various chromosomes. In contrast, hepatic tumors induced by a typical mutagen (diethylnitrosamine) followed by a nonmutagen (phenobarbital) did not show such mutational patterns. Additionally, structural alterations such as translocations were observed more frequently in the acetamide-induced tumors. Moreover, most of the acetamide-induced tumors expressed c-Myc and/or MDM2 protein due to the copy number gain of each oncogene. These results suggest the occurrence of chromosomal rearrangements and subsequent oncogene amplification in the acetamide-induced tumors. Taken together, the results indicate that CAs are directly involved in tumorigenesis through chromosomal rearrangements in an acetamide-induced hepatocarcinogenesis rat model.

Diagnostic challenges in complicated case of glioblastoma

Glioblastoma is the commonest primary malignant brain tumor, with a very poor prognosis and short overall survival. It is characterized by its high intra- and intertumoral heterogeneity, in terms of both the level of single-nucleotide variants, copy number alterations, and aneuploidy. Therefore, routine diagnosis can be challenging in some cases. We present a complicated case of glioblastoma, which was characterized with five cytogenomic methods: interphase fluorescence in situ hybridization, multiplex ligation-dependent probe amplification, comparative genomic hybridization array and single-nucleotide polymorphism, targeted gene panel, and whole-genome sequencing. These cytogenomic methods revealed classical findings associated with glioblastoma, such as a lack of IDH and TERT mutations, gain of chromosome 7, and loss of chromosome 10. At least three pathological clones were identified, including one with whole-genome duplication, and one with loss of 1p and suspected loss of 19q. Deletion and mutation of the TP53 gene were detected with numerous breakends on 17p and 20q. Based on these findings, we recommend a combined approach to the diagnosis of glioblastoma involving the detection of copy number alterations, mutations, and aneuploidy. The choice of the best combination of methods is based on cost, time required, staff expertise, and laboratory equipment. This integrated strategy could contribute directly to tangible improvements in the diagnosis, prognosis, and prediction of the therapeutic responses of patients with brain tumors.

Comparative analysis of methodologies for detecting extrachromosomal circular DNA

Extrachromosomal circular DNA (eccDNA) is crucial in oncogene amplification, gene transcription regulation, and intratumor heterogeneity. While various analysis pipelines and experimental methods have been developed for eccDNA identification, their detection efficiencies have not been systematically assessed. To address this, we evaluate the performance of 7 analysis pipelines using seven simulated datasets, in terms of accuracy, identity, duplication rate, and computational resource consumption. We also compare the eccDNA detection efficiency of 7 experimental methods through twenty-one real sequencing datasets. Here, we show that Circle-Map and Circle_finder (bwa-mem-samblaster) outperform the other short-read pipelines. However, Circle_finder (bwa-mem-samblaster) exhibits notable redundancy in its outcomes. CReSIL is the most effective pipeline for eccDNA detection in long-read sequencing data at depths higher than 10X. Moreover, long-read sequencing-based Circle-Seq shows superior efficiency in detecting copy number-amplified eccDNA over 10 kb in length. These results offer valuable insights for researchers in choosing the suitable methods for eccDNA research.

Real-World Performance of Integrative Clinical Genomics in Pediatric Precision Oncology

Despite significant improvement in the survival of pediatric patients with cancer, treatment outcomes for high-risk, relapsed, and refractory cancers remain unsatisfactory. Moreover, prolonged survival is frequently associated with long-term adverse effects due to intensive multimodal treatments. Accelerating the progress of pediatric oncology requires both therapeutic advances and strategies to mitigate the long-term cytotoxic side effects, potentially through targeting specific molecular drivers of pediatric malignancies. In this report, we present the results of integrative genomic and transcriptomic profiling of 230 patients with malignant solid tumors (the "primary cohort") and 18 patients with recurrent or otherwise difficult-to-treat nonmalignant conditions (the "secondary cohort"). The integrative workflow for the primary cohort enabled the identification of clinically significant single nucleotide variants, small insertions/deletions, and fusion genes, which were found in 55% and 28% of patients, respectively. For 38% of patients, molecularly informed treatment recommendations were made. In the secondary cohort, known or potentially driving alteration was detected in 89% of cases, including a suspected novel causal gene for patients with inclusion body infantile digital fibromatosis. Furthermore, 47% of findings also brought therapeutic implications for subsequent management. Across both cohorts, changes or refinements to the original histopathological diagnoses were achieved in 4% of cases. Our study demonstrates the efficacy of integrating advanced genomic and transcriptomic analyses to identify therapeutic targets, refine diagnoses, and optimize treatment strategies for challenging pediatric and young adult malignancies and underscores the need for broad implementation of precision oncology in clinical settings.

Machine Learning-Based Detection of Bladder Cancer by Urine cfDNA Fragmentation Hotspots that Capture Cancer-Associated Molecular Features

BACKGROUND

cfDNA fragmentomics-based liquid biopsy is a potential option for noninvasive bladder cancer (BLCA) detection that remains an unmet clinical need.

METHODS

We assessed the diagnostic performance of cfDNA hotspot-driven machine-learning models in a cohort of 55 BLCA patients, 51 subjects with benign conditions, and 11 healthy volunteers. We further performed functional bioinformatics analysis for biological understanding and interpretation of the tool's diagnostic capability.

RESULTS

Urinary cfDNA hotspots-based machine-learning model enabled effective BLCA detection, achieving high performance (area under curve 0.96) and an 87% sensitivity at 100% specificity. It outperformed models using other cfDNA-derived features. In stage-stratified analysis, the sensitivity at 100% specificity of the urine hotspots-based model was 71% and 92% for early (low-grade Ta and T1) and advanced (high-grade T1 and muscle-invasive) disease, respectively. Biologically, cfDNA hotspots effectively retrieved regulatory elements and were correlated with the cell of origin. Urine cfDNA hotspots specifically captured BLCA-related molecular features, including key functional pathways, chromosome loci associated with BLCA risk as identified in genome-wide association studies, or presenting frequent somatic alterations in BLCA tumors, and the transcription factor regulatory landscape.

CONCLUSIONS

Our findings support the applicability of urine cfDNA fragmentation hotspots for noninvasive BLCA diagnosis, as well as for future translational study regarding its molecular pathology and heterogeneity.

Identification of candidate genes associating with soybean cyst nematode in soybean (Glycine max L.) using BSA-seq

Soybean cyst nematode disease represents the major soil-borne disease of soybean. Identifying disease-resistant genes in soybean has a substantial impact on breeding of disease-resistant crops and genetic improvement. The present work created the F2 population with the disease-resistant line H-10 and disease-susceptible line Chidou4. 30 respective F2disease-resistant and disease-susceptible individuals for forming two DNA pools for whole-genome re-sequencing were selected. As a result, a total of 11,522,230 single nucleotide polymorphism (SNPs) markers from these two parental lines and two mixed pools were obtained. Accordng to SNP-index based association analysis, there were altogether 741 genes out of 99% confidence interval, which were mainly enriched into regions of 38,524,128∼39,849,988 bp with a total length of 1.33 Mb contain 111 genes on chromosome 2, 27,821,012∼29,612,574 bp with a total length of 1.79 Mb contain 92 genes on chromosome 3, 308∼348,214 bp with a total of length 0.35 Mb contain 34 genes on chromosome 10, and 53,867,581∼58,017, 852 bp with a total length of 4.15 Mb contain 504 genes on chromosome 18. Bulk segregant analysis in F2 generations (BSA-seq) was correlated with a disease resistance interval containing 15 genes. Then, using bioinformatics analysis and differential expression analysis, five candidate genes were identified: Glyma.02G211400, Glyma.18G252800, Glyma.18G285800, Glyma.18G287400 and Glyma.18G298200. Our results provides a key basis for analyzing the soybean resistance mechanism against soybean cyst nematode and cloning soybean resistance genes.

Novel molecular, structural and clinical findings in an Italian cohort of congenital cataract

The current genetic diagnostic workup of congenital cataract (CC) is mainly based on NGS panels, whereas exome sequencing (ES) has occasionally been employed. In this multicentre study, we investigated by ES the detection yield, mutational spectrum and genotype-phenotype correlations in a CC cohort recruited between 2020 and mid-2022. The cohort consisted of 67 affected individuals from 51 unrelated families and included both non-syndromic (75%) and syndromic (25%) phenotypes, with extra-CC ocular/visual features present in both groups (48% and 76%, respectively). The functional effect of variants was predicted by 3D modelling and hydropathy properties changes. Variant clustering was used for the in-depth assessment of genotype-phenotype correlations. A diagnostic (pathogenic or likely pathogenic) variant was identified in 19 out of 51 probands/families (~37%). In a further 14 probands/families a candidate variant was identified: in 12 families a VUS was detected, of which 9 were considered plausibly pathogenic (i.e., 4 or 5 points according to ACMG criteria), while in 2 probands ES identified a single variant in an autosomal recessive gene associated with CC. Eighteen probands/families, manifesting primarily non-syndromic CC (15/18, 83%), remained unsolved. The identified variants (8 P, 12 LP, 10 VUS-PP, and 5 VUS), half of which were unreported in the literature, affected five functional categories of genes involved in transcription/splicing, lens formation/homeostasis (i.e., crystallin genes), membrane signalling, cell-cell interaction, and immune response. A phenotype-specific variant clustering was observed in four genes (KIF1A, MAF, PAX6, SPTAN1), whereas variable expressivity and potential phenotypic expansion in two (BCOR, NHS) and five genes (CWC27, KIF1A, IFIH1, PAX6, SPTAN1), respectively. Finally, ES allowed to detect variants in six genes not commonly included in commercial CC panels. These findings broaden the genotype-phenotype correlations in one of the largest CC cohorts tested by ES, providing novel insights into the underlying pathogenetic mechanisms and emphasising the power of ES as first-tier test.

Mutational Landscape of Gastric Adenocarcinoma of the Fundic Gland Type Revealed by Whole Genome Sequencing

BACKGROUND

Gastric adenocarcinoma of the fundic gland type (GA-FG) is a newly described variant of gastric adenocarcinoma with lack of knowledges regarding its genetic features.

METHODS

We performed whole-genome sequencing (WGS) in formalin-fixed paraffin-embedded (FFPE) tumor tissues and matched adjacent noncancerous specimens from 21 patients with GA-FG, and integrated published datasets from 1105 patients with traditional gastric adenocarcinoma with the purpose of dissecting genetic determinants both common to conventional gastric adenocarcinoma and unique to GA-FG disease.

RESULTS

We characterized the genomic architecture of GA-FG disease, revealing the predominant proportion of C > T substitution among the six types of SNVs. GNAS was the most significantly mutated driver gene (14.29%). 42.8% of samples harbored "Kataegis." Distinct genomic alterations between GA-FG and conventional gastric cancer were identified. Specifically, low mutational burden and relatively moderate mutational frequencies of significantly mutated driver genes, coupled with the absence of non-silent alterations of formerly well-known drivers such as TP53, PIK3CA and KRAS were identified in GA-FG patients. Oncogenic signaling pathway analysis revealed mutational processes associated with focal adhesions and proteoglycans in cancer, highlighting both common and specific procedures during the development of GA-FG and conventional gastric cancer.

CONCLUSION

Our study is the first to comprehensively depict the genomic landscape highlighting the multidimensional perturbations in GA-FG patients. These discoveries offered mechanistic insights for novel diagnostic and therapeutic strategies for patients with such disease.

Gelada genomes highlight events of gene flow, hybridisation and local adaptation that track past climatic changes

Theropithecus gelada, the last surviving species of this genus, occupy a unique and highly specialised ecological niche in the Ethiopian highlands. A subdivision into three geographically defined populations (Northern, Central and Southern) has been tentatively proposed for this species on the basis of genetic analyses, but genomic data have been investigated only for two of these groups (Northern and Central). Here we combined newly generated whole genome sequences of individuals sampled from the population living south of the East Africa Great Rift Valley with available data from the other two gelada populations to reconstruct the evolutionary history of the species. Integrating genomic and paleoclimatic data we found that gene-flow across populations and with Papio species tracked past climate changes. The isolation and climatic conditions experienced by Southern geladas during the Holocene shaped local diversity and generated diet-related genomic signatures.

Whole-Exome Sequencing of Vulvar Squamous Cell Carcinomas Reveals an Impaired Prognosis in Patients With TP53 Mutations and Concurrent CCND1 Gains

Very little information is available on the mutational landscape of vulvar squamous cell carcinoma (VSCC), a disease that mainly affects older women. Studies focusing on the mutational patterns of the currently recognized etiopathogenic types of this tumor (human papillomavirus [HPV]-associated [HPV-A], HPV-independent [HPV-I] with TP53 mutation [HPV-I/TP53mut], and HPV-I with wild-type TP53 [HPV-I/TP53wt]) are particularly rare, and there is almost no information on the prognostic implications of these abnormalities.Whole-exome DNA sequencing of 60 VSCC and matched normal tissues from each patient was performed. HPV detection, immunohistochemistry (IHC) for p16, p53, and mismatch repair proteins were also performed. Ten tumors (16.7%) were classified as HPV-A, 37 (61.7%) as HPV-I/TP53mut, and 13 (21.6%) as HPV-I/TP53wt. TP53 was the most frequently mutated gene (66.7%), followed by FAT1 (28.3%), CDKN2A (25.0%), RNF213 (23.3%), NFE2L2 (20%) and PIK3CA (20%). All the 60 tumors (100%) were DNA mismatch repair proficient. Seventeen tumors (28.3%) showed CCND1 gain. Bivariate analysis, adjusted for International Federation of Gynecology and Obstetrics stage, revealed that TP53 mutation, CCND1 gain, and the combination of the 2 alterations were strongly associated with impaired recurrence-free survival (hazard ratio, 4.4; P < .001) and disease-specific survival (hazard ratio, 6.1; P = .002). Similar results were obtained when p53 IHC status was used instead of TP53 status and when considering only HPV-I VSCC. However, in the latter category, p53 IHC maintained its prognostic impact only in combination with CCND1 gains. All tumors carried at least one potentially actionable genomic alteration. In conclusion, VSCCs with CCND1 gain represent a prognostically adverse category among HPV-I/TP53mut tumors. All patients with VSCCs are potential candidates for targeted therapy.

Single-cell RNA sequencing of pediatric Hodgkin lymphoma to study the inhibition of T cell subtypes

Pediatric classic Hodgkin lymphoma (cHL) patients have a high survival rate but suffer from severe long-term side effects induced by chemo- and radiotherapy. cHL tumors are characterized by the low fraction (0.1%-10%) of malignant Hodgkin and Reed-Sternberg (HRS) cells in the tumor. The HRS cells depend on the surrounding immune cells for survival and growth. This dependence is leveraged by current treatments that target the PD-1/PD-L1 axis in cHL tumors. The development of more targeted therapies that are specific for the tumor and are therefore less toxic for healthy tissue compared with conventional chemotherapy could improve the quality of life of pediatric cHL survivors. Here, we applied single-cell RNA sequencing (scRNA-seq) on isolated HRS cells and the immune cells from the same cHL tumors. Besides TNFRSF8 (CD30), we identified other genes of cell surface proteins that are consistently overexpressed in HRS cells, such as NRXN3 and LRP8, which can potentially be used as alternative targets for antibody-drug conjugates or CAR T cells. Finally, we identified potential interactions by which HRS cells inhibit T cells, among which are the galectin-1/CD69 and HLA-II/LAG3 interactions. RNAscope was used to validate the enrichment of CD69 and LAG3 expression on T cells near HRS cells and indicated large variability of the interaction strength with the corresponding ligands between patients and between tumor tissue regions. In conclusion, this study identifies new potential therapeutic targets for cHL and highlights the importance of studying heterogeneity when identifying therapy targets, specifically those that target tumor-immune cell interactions.

Implementing next-generation sequencing for diagnosis and management of hereditary hearing impairment: a comprehensive review

INTRODUCTION

Sensorineural hearing impairment (SNHI), a common childhood disorder with heterogeneous genetic causes, can lead to delayed language development and psychosocial problems. Next-generation sequencing (NGS) offers high-throughput screening and high-sensitivity detection of genetic etiologies of SNHI, enabling clinicians to make informed medical decisions, provide tailored treatments, and improve prognostic outcomes.

AREAS COVERED

This review covers the diverse etiologies of HHI and the utility of different NGS modalities (targeted sequencing and whole exome/genome sequencing), and includes HHI-related studies on newborn screening, genetic counseling, prognostic prediction, and personalized treatment. Challenges such as the trade-off between cost and diagnostic yield, detection of structural variants, and exploration of the non-coding genome are also highlighted.

EXPERT OPINION

In the current landscape of NGS-based diagnostics for HHI, there are both challenges (e.g. detection of structural variants and non-coding genome variants) and opportunities (e.g. the emergence of medical artificial intelligence tools). The authors advocate the use of technological advances such as long-read sequencing for structural variant detection, multi-omics analysis for non-coding variant exploration, and medical artificial intelligence for pathogenicity assessment and outcome prediction. By integrating these innovations into clinical practice, precision medicine in the diagnosis and management of HHI can be further improved.

Primary unilateral macronodular adrenal hyperplasia with concomitant glucocorticoid and androgen excess and KDM1A inactivation

BACKGROUND

Primary bilateral macronodular adrenal hyperplasia (PBMAH) is a rare cause of Cushing's syndrome. Individuals with PBMAH and glucose-dependent insulinotropic polypeptide (GIP)-dependent Cushing's syndrome due to ectopic expression of the GIP receptor (GIPR) typically harbor inactivating KDM1A sequence variants. Primary unilateral macronodular adrenal hyperplasia (PUMAH) with concomitant glucocorticoid and androgen excess has never been encountered or studied.

METHODS

We investigated a woman with a large, heterogeneous adrenal mass and severe adrenocorticotropic hormone-independent glucocorticoid and androgen excess, a biochemical presentation typically suggestive of adrenocortical carcinoma. The patient presented during pregnancy (22nd week of gestation) and reported an 18-month history of oligomenorrhea, hirsutism, and weight gain. We undertook an exploratory study with detailed histopathological and genetic analysis of the resected adrenal mass and leukocyte DNA collected from the patient and her parents.

RESULTS

Histopathology revealed benign macronodular adrenal hyperplasia. Imaging showed a persistently normal contralateral adrenal gland. Whole-exome sequencing of 4 representative nodules detected KDM1A germline variants, benign NM_001009999.3:c.136G > A:p.G46S, and likely pathogenic NM_001009999.3:exon6:c.865_866del:p.R289Dfs*7. Copy number variation analysis demonstrated an additional somatic loss of the KDM1A wild-type allele on chromosome 1p36.12 in all nodules. RNA sequencing of a representative nodule showed low/absent KDM1A expression and increased GIPR expression compared with 52 unilateral sporadic adenomas and 4 normal adrenal glands. Luteinizing hormone/chorionic gonadotropin receptor expression was normal. Sanger sequencing confirmed heterozygous KDM1A variants in both parents (father: p.R289Dfs*7 and mother: p.G46S) who showed no clinical features suggestive of glucocorticoid or androgen excess.

CONCLUSIONS

We investigated the first PUMAH associated with severe Cushing's syndrome and concomitant androgen excess, suggesting pathogenic mechanisms involving KDM1A.

A human fetal cerebellar map of the late second trimester reveals developmental molecular characteristics and abnormality in trisomy 21

Our understanding of human fetal cerebellum development during the late second trimester, a critical period for the generation of astrocytes, oligodendrocytes, and unipolar brush cells (UBCs), remains limited. Here, we performed single-cell RNA sequencing (scRNA-seq) in human fetal cerebellum samples from gestational weeks (GWs) 18-25. We find that proliferating UBC progenitors distribute in the subventricular zone of the rhombic lip (RLSVZ) near white matter (WM), forming a layer structure. We also delineate two trajectories from astrogenic radial glia (ARGs) to Bergmann glial progenitors (BGPs) and recognize oligodendrogenic radial glia (ORGs) as one source of primitive oligodendrocyte progenitor cells (PriOPCs). Additionally, our scRNA-seq analysis of the trisomy 21 fetal cerebellum at this stage reveals abnormal upregulated genes in pathways such as the cell adhesion pathway and focal adhesion pathway, which potentially promote neuronal differentiation. Overall, our research provides valuable insights into normal and abnormal development of the human fetal cerebellum.

STAG2 mutations reshape the cohesin-structured spatial chromatin architecture to drive gene regulation in acute myeloid leukemia

Cohesin shapes the chromatin architecture, including enhancer-promoter interactions. Its components, especially STAG2, but not its paralog STAG1, are frequently mutated in myeloid malignancies. To elucidate the underlying mechanisms of leukemogenesis, we comprehensively characterized genetic, transcriptional, and chromatin conformational changes in acute myeloid leukemia (AML) patient samples. Specific loci displayed altered cohesin occupancy, gene expression, and local chromatin activation, which were not compensated by the remaining STAG1-cohesin. These changes could be linked to disrupted spatial chromatin looping in cohesin-mutated AMLs. Complementary depletion of STAG2 or STAG1 in primary human hematopoietic progenitors (HSPCs) revealed effects resembling STAG2-mutant AML-specific changes following STAG2 knockdown, not invoked by the depletion of STAG1. STAG2-deficient HSPCs displayed impaired differentiation capacity and maintained HSPC-like gene expression. This work establishes STAG2 as a key regulator of chromatin contacts, gene expression, and differentiation in the hematopoietic system and identifies candidate target genes that may be implicated in human leukemogenesis.

Autosomal dominant stromal corneal dystrophy associated with a SPARCL1 missense variant

Corneal dystrophies are phenotypically and genetically heterogeneous, often resulting in visual impairment caused by corneal opacification. We investigated the genetic cause of an autosomal dominant corneal stromal dystrophy in a pedigree with eight affected individuals in three generations. Affected individuals had diffuse central stromal opacity, with reduced visual acuity in older family members. Histopathology of affected cornea tissue removed during surgery revealed mild stromal textural alterations with alcianophilic deposits. Whole genome sequence data were generated for four affected individuals. No rare variants (MAF < 0.001) were identified in established corneal dystrophy genes. However, a novel heterozygous missense variant in exon 4 of SPARCL1, NM_004684: c.334G > A; p.(Glu112Lys), which is predicted to be damaging, segregated with disease. SPARC-like protein 1 (SPARCL1) is a secreted matricellular protein involved in cell migration, cell adhesion, tissue repair, and remodelling. Interestingly, SPARCL1 has been shown to regulate decorin. Heterozygous variants in DCN, encoding decorin, cause autosomal dominant congenital stromal corneal dystrophy, suggesting a common pathogenic pathway. Therefore, we performed immunohistochemistry to compare SPARCL1 and decorin localisation in corneal tissue from an affected family member and an unaffected control. Strikingly, the level of decorin was significantly decreased in the corneal stroma of the affected tissue, and SPARCL1 appeared to be retained in the epithelium. In summary, we describe a novel autosomal dominant corneal stromal dystrophy associated with a missense variant in SPARCL1, extending the phenotypic and genetic heterogeneity of inherited corneal disease.

Centrosomal Protein 55 Regulates Chromosomal Instability in Cancer Cells by Controlling Microtubule Dynamics

Centrosomal Protein 55 (CEP55) exhibits various oncogenic activities; it regulates the PI3K-Akt-pathway, midbody abscission, and chromosomal instability (CIN) in cancer cells. Here, we analyzed the mechanism of how CEP55 controls CIN in ovarian and breast cancer (OvCa) cells. Down-regulation of CEP55 reduced CIN in all cell lines analyzed, and CEP55 depletion decreased spindle microtubule (MT)-stability in OvCa cells. Moreover, recombinant CEP55 accelerated MT-polymerization and attenuated cold-induced MT-depolymerization. To analyze a potential relationship between CEP55-controlled CIN and its impact on MT-stability, we identified the CEP55 MT-binding peptides inside the CEP55 protein. Thereafter, a mutant with deficient MT-binding activity was re-expressed in CEP55-depleted OvCa cells and we could show that this mutant did not restore reduced CIN in CEP55-depleted cells. This finding strongly indicates that CEP55 regulates CIN by controlling MT dynamics.

Towards Digital Quantification of Ploidy from Pan-Cancer Digital Pathology Slides using Deep Learning

Abnormal DNA ploidy, found in numerous cancers, is increasingly being recognized as a contributor in driving chromosomal instability, genome evolution, and the heterogeneity that fuels cancer cell progression. Furthermore, it has been linked with poor prognosis of cancer patients. While next-generation sequencing can be used to approximate tumor ploidy, it has a high error rate for near-euploid states, a high cost and is time consuming, motivating alternative rapid quantification methods. We introduce PloiViT, a transformer-based model for tumor ploidy quantification that outperforms traditional machine learning models, enabling rapid and cost-effective quantification directly from pathology slides. We trained PloiViT on a dataset of fifteen cancer types from The Cancer Genome Atlas and validated its performance in multiple independent cohorts. Additionally, we explored the impact of self-supervised feature extraction on performance. PloiViT, using self-supervised features, achieved the lowest prediction error in multiple independent cohorts, exhibiting better generalization capabilities. Our findings demonstrate that PloiViT predicts higher ploidy values in aggressive cancer groups and patients with specific mutations, validating PloiViT potential as complementary for ploidy assessment to next-generation sequencing data. To further promote its use, we release our models as a user-friendly inference application and a Python package for easy adoption and use.

Identification of potent and reversible piperidine carboxamides that are species-selective orally active proteasome inhibitors to treat malaria

Malaria remains a global health concern as drug resistance threatens treatment programs. We identified a piperidine carboxamide (SW042) with anti-malarial activity by phenotypic screening. Selection of SW042-resistant Plasmodium falciparum (Pf) parasites revealed point mutations in the Pf_proteasome β5 active-site (Pfβ5). A potent analog (SW584) showed efficacy in a mouse model of human malaria after oral dosing. SW584 had a low propensity to generate resistance (minimum inoculum for resistance [MIR] >109) and was synergistic with dihydroartemisinin. Pf_proteasome purification was facilitated by His8-tag introduction onto β7. Inhibition of Pfβ5 correlated with parasite killing, without inhibiting human proteasome isoforms or showing cytotoxicity. The Pf_proteasome_SW584 cryoelectron microscopy (cryo-EM) structure showed that SW584 bound non-covalently distal from the catalytic threonine, in an unexplored pocket at the β5/β6/β3 subunit interface that has species differences between Pf and human proteasomes. Identification of a reversible, species selective, orally active series with low resistance propensity provides a path for drugging this essential target.

Drug resistance and genetic characteristics of one Eimeria tenella isolate from Xiantao, Hubei Province, China

Chicken coccidiosis causes retarded growth and low production performance in poultry, resulting in huge economic losses to the poultry industry. In order to prevent and control chicken coccidiosis, great efforts have been made to develop new drugs and vaccines, which require pure isolates of Eimeria spp. In this study, we obtained the Eimeira tenella Xiantao isolate by single oocyst isolation technology and compared its genome with the reference genome GCF_000499545.2_ETH001 of the Houghton strain. The results of the comparative genomic analysis indicated that the genome of this isolate contained 46,888 single nucleotide polymorphisms (SNPs). There were 15,107 small insertion and deletion variations (indels), 1693 structural variations (SV), and 3578 copy number variations (CNV). In addition, 64 broilers were used to determine the resistance profile of Xiantao strain. Drug susceptibility testing revealed that this isolate was completely resistant to monensin, diclazuril, halofuginone, sulfachlorpyrazine sodium, and toltrazuril, but sensitive to decoquinate. These data improve our understanding of drug resistance in avian coccidia.

Selection for robust metabolism in domesticated yeasts is driven by adaptation to Hsp90 stress

Protein folding both promotes and constrains adaptive evolution. We uncover this surprising duality in the role of the protein-folding chaperone heat shock protein 90 (Hsp90) in maintaining the integrity of yeast metabolism amid proteotoxic stressors within industrial domestication niches. Ethanol disrupts critical Hsp90-dependent metabolic pathways and exerts strong selective pressure for redundant duplications of key genes within these pathways, yielding the classical genomic signatures of beer and bread domestication. This work demonstrates a mechanism of adaptive canalization in an ecology of major economic importance and highlights Hsp90-dependent variation as an important source of phantom heritability in complex traits.

The Open Pediatric Cancer Project

Background

In 2019, the Open Pediatric Brain Tumor Atlas (OpenPBTA) was created as a global, collaborative open-science initiative to genomically characterize 1,074 pediatric brain tumors and 22 patient-derived cell lines. Here, we extend the OpenPBTA to create the Open Pediatric Cancer (OpenPedCan) Project, a harmonized open-source multi-omic dataset from 6,112 pediatric cancer patients with 7,096 tumor events across more than 100 histologies. Combined with RNA-Seq from the Genotype-Tissue Expression (GTEx) and The Cancer Genome Atlas (TCGA), OpenPedCan contains nearly 48,000 total biospecimens (24,002 tumor and 23,893 normal specimens).

Findings

We utilized Gabriella Miller Kids First (GMKF) workflows to harmonize WGS, WXS, RNA-seq, and Targeted Sequencing datasets to include somatic SNVs, InDels, CNVs, SVs, RNA expression, fusions, and splice variants. We integrated summarized CPTAC whole cell proteomics and phospho-proteomics data, miRNA-Seq data, and have developed a methylation array harmonization workflow to include m-values, beta-vales, and copy number calls. OpenPedCan contains reproducible, dockerized workflows in GitHub, CAVATICA, and Amazon Web Services (AWS) to deliver harmonized and processed data from over 60 scalable modules which can be leveraged both locally and on AWS. The processed data are released in a versioned manner and accessible through CAVATICA or AWS S3 download (from GitHub), and queryable through PedcBioPortal and the NCI's pediatric Molecular Targets Platform. Notably, we have expanded PBTA molecular subtyping to include methylation information to align with the WHO 2021 Central Nervous System Tumor classifications, allowing us to create research- grade integrated diagnoses for these tumors.

Conclusions

OpenPedCan data and its reproducible analysis module framework are openly available and can be utilized and/or adapted by researchers to accelerate discovery, validation, and clinical translation.

Epigenetic alterations affecting hematopoietic regulatory networks as drivers of mixed myeloid/lymphoid leukemia

Leukemias with ambiguous lineage comprise several loosely defined entities, often without a clear mechanistic basis. Here, we extensively profile the epigenome and transcriptome of a subgroup of such leukemias with CpG Island Methylator Phenotype. These leukemias exhibit comparable hybrid myeloid/lymphoid epigenetic landscapes, yet heterogeneous genetic alterations, suggesting they are defined by their shared epigenetic profile rather than common genetic lesions. Gene expression enrichment reveals similarity with early T-cell precursor acute lymphoblastic leukemia and a lymphoid progenitor cell of origin. In line with this, integration of differential DNA methylation and gene expression shows widespread silencing of myeloid transcription factors. Moreover, binding sites for hematopoietic transcription factors, including CEBPA, SPI1 and LEF1, are uniquely inaccessible in these leukemias. Hypermethylation also results in loss of CTCF binding, accompanied by changes in chromatin interactions involving key transcription factors. In conclusion, epigenetic dysregulation, and not genetic lesions, explains the mixed phenotype of this group of leukemias with ambiguous lineage. The data collected here constitute a useful and comprehensive epigenomic reference for subsequent studies of acute myeloid leukemias, T-cell acute lymphoblastic leukemias and mixed-phenotype leukemias.

Non-homologous end joining shapes the genomic rearrangement landscape of chromothripsis from mitotic errors

Mitotic errors generate micronuclei entrapping mis-segregated chromosomes, which are susceptible to catastrophic fragmentation through chromothripsis. The reassembly of fragmented chromosomes by error-prone DNA double-strand break (DSB) repair generates diverse genomic rearrangements associated with human diseases. How specific repair pathways recognize and process these lesions remains poorly understood. Here we use CRISPR/Cas9 to systematically inactivate distinct DSB repair pathways and interrogate the rearrangement landscape of fragmented chromosomes. Deletion of canonical non-homologous end joining (NHEJ) components substantially reduces complex rearrangements and shifts the rearrangement landscape toward simple alterations without the characteristic patterns of chromothripsis. Following reincorporation into the nucleus, fragmented chromosomes localize within sub-nuclear micronuclei bodies (MN bodies) and undergo ligation by NHEJ within a single cell cycle. In the absence of NHEJ, chromosome fragments are rarely engaged by alternative end-joining or recombination-based mechanisms, resulting in delayed repair kinetics, persistent 53BP1-labeled MN bodies, and cell cycle arrest. Thus, we provide evidence supporting NHEJ as the exclusive DSB repair pathway generating complex rearrangements from mitotic errors.

Genomic characteristics and evolution of Multicentric Esophageal and gastric Cardiac Cancer

BACKGROUND

Esophageal carcinoma (EC) and gastric cardiac adenocarcinoma (GCA) have high incidence rates in the Chaoshan region of South China. Multifocal esophageal and cardiac cancer (MECC) is commonly observed in this region in clinical practice. However, the genomic characteristics of MECC remains unclear.

MATERIALS AND METHODS

In this study, a total of 2123 clinical samples of EC and GCA were analyzed to determine the frequency of multifocal tumors, as well as their occurrence sites and pathological types. Cox proportional hazards regression was used to model the relationship between age, sex, and tumor state concerning survival in our analysis of the cohort of 541 patients with available follow-up data. We performed whole-genome sequencing on 20 tumor foci and 10 normal samples from 10 MECC patients to infer clonal structure on 6 MECC patients to explore genome characteristics.

RESULT

The MECC rate of EC and GCA was 5.65% (121 of 2123). Age and sex were potential factors that may influence the risk of MECC (p < 0.001). Furthermore, MECC patients showed worse survival compared with single tumor patients. We found that 12 foci from 6 patients were multicentric origin model (MC), which exhibited significant heterogeneity of variations in paired foci and had an increased number of germline mutations in immune genes compared to metastatic model. In MC cases, different lesions in the same patient were driven by distinct mutation and copy number variation (CNV) events. Although TP53 and other driver mutation genes have a high frequency in the samples, their mutation sites show significant heterogeneity in paired tumor specimens. On the other hand, CNV genes exhibited higher concordance in paired samples, especially in the amplification of oncogenes and the deletion of tumor suppressor genes.

CONCLUSIONS

The extent of inter-tumor heterogeneity suggests both monoclonal and polyclonal origins of MECC, which could provide insight into the genome diversity of MECC and guide clinical implementation.

Single-cell transcriptome dissecting the microenvironment remodeled by PD1 blockade combined with photodynamic therapy in a mouse model of oral carcinogenesis

Oral squamous cell carcinoma (OSCC) stands as a predominant and perilous malignant neoplasm globally, with the majority of cases originating from oral potential malignant disorders (OPMDs). Despite this, effective strategies to impede the progression of OPMDs to OSCC remain elusive. In this study, we established mouse models of oral carcinogenesis via 4-nitroquinoline 1-oxide induction, mirroring the sequential transformation from normal oral mucosa to OPMDs, culminating in OSCC development. By intervening during the OPMDs stage, we observed that combining PD1 blockade with photodynamic therapy (PDT) significantly mitigated oral carcinogenesis progression. Single-cell transcriptomic sequencing unveiled microenvironmental dysregulation occurring predominantly from OPMDs to OSCC stages, fostering a tumor-promoting milieu characterized by increased Treg proportion, heightened S100A8 expression, and decreased Fib_Igfbp5 (a specific fibroblast subtype) proportion, among others. Notably, intervening with PD1 blockade and PDT during the OPMDs stage hindered the formation of the tumor-promoting microenvironment, resulting in decreased Treg proportion, reduced S100A8 expression, and increased Fib_Igfbp5 proportion. Moreover, combination therapy elicited a more robust treatment-associated immune response compared with monotherapy. In essence, our findings present a novel strategy for curtailing the progression of oral carcinogenesis.

Cellular adaptation to cancer therapy along a resistance continuum

Advancements in precision oncology over the past decades have led to new therapeutic interventions, but the efficacy of such treatments is generally limited by an adaptive process that fosters drug resistance1. In addition to genetic mutations2, recent research has identified a role for non-genetic plasticity in transient drug tolerance3 and the acquisition of stable resistance4,5. However, the dynamics of cell-state transitions that occur in the adaptation to cancer therapies remain unknown and require a systems-level longitudinal framework. Here we demonstrate that resistance develops through trajectories of cell-state transitions accompanied by a progressive increase in cell fitness, which we denote as the 'resistance continuum'. This cellular adaptation involves a stepwise assembly of gene expression programmes and epigenetically reinforced cell states underpinned by phenotypic plasticity, adaptation to stress and metabolic reprogramming. Our results support the notion that epithelial-to-mesenchymal transition or stemness programmes-often considered a proxy for phenotypic plasticity-enable adaptation, rather than a full resistance mechanism. Through systematic genetic perturbations, we identify the acquisition of metabolic dependencies, exposing vulnerabilities that can potentially be exploited therapeutically. The concept of the resistance continuum highlights the dynamic nature of cellular adaptation and calls for complementary therapies directed at the mechanisms underlying adaptive cell-state transitions.

Genetic mosaic of the Mediterranean fig: comprehensive genomic insights from a gene bank collection

High-depth whole-genome resequencing of 53 diverse fig tree genotypes yielded a rich dataset of genetic variants. We successfully identified 5,501,460 single-nucleotide polymorphisms (SNPs) and 1,228,537 insertions and deletions (InDels), providing a high-density and excellent-quality genetic map of the fig tree. We also performed a detailed population structure analysis, dividing the 53 genotypes into three geographical groups and assessing their genetic diversity and divergence. Analysis of structural variants (SVs) and copy number variations (CNVs) revealed their potential functional impact, particularly in plant-pathogen interaction and secondary metabolism. Metabolomic fingerprinting of fig genotypes uncovered extensive variation in primary metabolites and polyphenolic compounds, highlighting the influence of genotype on fruit quality traits such as nutritional content and bioactive compound composition. The genome-wide association study (GWAS) identified critical SNPs associated with fruit quality and morphological features. The discovery of significant candidate genes, such as AGL62, GDSL, and COBRA-like protein 4 genes, offers promising targets for marker-assisted selection and genome editing approaches to improve fig fruit morphological and quality traits. This extensive genomic analysis of fig trees enhances our understanding of the genetic basis of important agronomic traits and provides a rich resource for future research in this economically and nutritionally significant fruit.

Establishment of a patient-specific avatar organoid model derived from EUS-guided fine-needle biopsy for timely clinical application in pancreatic ductal adenocarcinoma (with video)

BACKGROUND AND AIMS

Pancreatic ductal adenocarcinoma (PDAC) has the worst survival rate among tumors. At the time of diagnosis, more than 80% of PDACs are considered to be surgically unresectable, and there is an unmet need for treatment options in these inoperable PDACs. This study aimed to establish a patient-derived organoid (PDO) platform from EUS-guided fine-needle biopsy (EUS-FNB) collected at diagnosis and to determine its clinical applicability for the timely treatment of unresectable PDAC.

METHODS

Patients with suspected PDAC were prospectively enrolled at the Samsung Medical Center from 2015 to 2019. PDAC tissues were acquired by means of EUS-FNB to establish PDAC PDOs, which were comprehensively analyzed for histology, genomic sequencing, and high-throughput screening (HTS) drug sensitivity test.

RESULTS

PDAC PDOs were established with a success rate of 83.2% (94/113). It took approximately 3 weeks from acquiring minimal EUS-FNB specimens to generating sufficient PDAC PDOs for the simultaneous HTS drug sensitivity test and genomic sequencing. The high concordance between PDAC tissues and matched PDOs was confirmed, and whole-exome sequencing revealed the increased detection of genetic alterations in PDOs compared with EUS-FNB tissues. The HTS drug sensitivity test showed clinical correlation between the ex vivo PDO response and the actual chemotherapeutic response of the study patients in the real world (13 out of 15 cases). In addition, whole-transcriptome sequencing identified candidate genes associated with nab-paclitaxel resistance, such as ITGB7, ANPEP, and ST3GAL1.

CONCLUSIONS

This PDAC PDO platform allows several therapeutic drugs to be tested within a short time window and opens the possibility for timely personalized medicine as a "patient avatar model" in clinical practice.

Tumor mutation burden and FAT3 mutation influence long-term survival in surgically resected small cell lung cancer

Background

Small cell lung cancer (SCLC) is highly malignant and has a higher risk of recurrence even in patients who undergo early surgery. However, a subgroup of patients survived for many years. So far, the factors that determine the long-term survivorship remain largely unknown. To determine the genetic characteristics of long-term survival (LTS) after surgery in SCLC, we performed comprehensive comparative genomic profiling and tumor mutation burden (TMB) analysis of resected tumor tissues from patients with LTS and short-term survival (STS) after surgery.

Methods

The present study screened 11 patients from 52 patients with SCLC who underwent surgery at Zhejiang Cancer Hospital from April 2008 to December 2017. A total of six LTS patients (≥4 years) with stage IIB or IIIA SCLC and five STS patients (<2 years) with stage IA or IB SCLC were included in the study. The STS patients were used as a control. All the patients underwent resection without neoadjuvant therapy. We assessed the genomic profiles of the resected tumor tissues and calculated the TMB using next-generation sequencing. We then analyzed and compared the molecular characteristics between the LTS and STS groups.

Results

Our data indicated that tumor tissues from patients with LTS harbor a high TMB. The median TMB for LTS patients was high (approximately 16.4 mutations/Mb), while that for STS patients was low (approximately 8.5 mutations/Mb). The median TMB of patients with LTS and STS showed a trend of significant difference (P=0.08). Gene alterations characterized the survival differences between the two groups. The FAT3 mutation was only found in the LTS group, and the P value determined by Fisher's exact test was 0.06.

Conclusions

A high non-synonymous TMB and the FAT3 mutation could potentially influence LTS after SCLC resection. This study provides valuable information about the molecular differences between LTS and STS patients. Studies with larger sample sizes need to be conducted to confirm our findings in the future.

Evaluation of somatic copy number variation detection by NGS technologies and bioinformatics tools on a hyper-diploid cancer genome

BACKGROUND

Copy number variation (CNV) is a key genetic characteristic for cancer diagnostics and can be used as a biomarker for the selection of therapeutic treatments. Using data sets established in our previous study, we benchmark the performance of cancer CNV calling by six most recent and commonly used software tools on their detection accuracy, sensitivity, and reproducibility. In comparison to other orthogonal methods, such as microarray and Bionano, we also explore the consistency of CNV calling across different technologies on a challenging genome.

RESULTS

While consistent results are observed for copy gain, loss, and loss of heterozygosity (LOH) calls across sequencing centers, CNV callers, and different technologies, variation of CNV calls are mostly affected by the determination of genome ploidy. Using consensus results from six CNV callers and confirmation from three orthogonal methods, we establish a high confident CNV call set for the reference cancer cell line (HCC1395).

CONCLUSIONS

NGS technologies and current bioinformatics tools can offer reliable results for detection of copy gain, loss, and LOH. However, when working with a hyper-diploid genome, some software tools can call excessive copy gain or loss due to inaccurate assessment of genome ploidy. With performance matrices on various experimental conditions, this study raises awareness within the cancer research community for the selection of sequencing platforms, sample preparation, sequencing coverage, and the choice of CNV detection tools.

Structure and transcription of integrated HPV DNA in vulvar carcinomas

HPV infections are associated with a fraction of vulvar cancers. Through hybridization capture and DNA sequencing, HPV DNA was detected in five of thirteen vulvar cancers. HPV16 DNA was integrated into human DNA in three of the five. The insertions were in introns of human NCKAP1, C5orf67, and LRP1B. Integrations in NCKAP1 and C5orf67 were flanked by short direct repeats in the human DNA, consistent with HPV DNA insertions at sites of abortive, staggered, endonucleolytic incisions. The insertion in C5orf67 was present as a 36 kbp, human-HPV-hetero-catemeric DNA as either an extrachromosomal circle or a tandem repeat within the human genome. The human circularization/repeat junction was defined at single nucleotide resolution. The integrated viral DNA segments all retained an intact upstream regulatory region and the adjacent viral E6 and E7 oncogenes. RNA sequencing revealed that the only HPV genes consistently transcribed from the integrated viral DNAs were E7 and E6*I. The other two HPV DNA+ tumors had coinfections, but no evidence for integration. HPV-positive and HPV-negative vulvar cancers exhibited contrasting human, global gene expression patterns partially overlapping with previously observed differences between HPV-positive and HPV-negative cervical and oropharyngeal cancers. A substantial fraction of the differentially expressed genes involved immune system function. Thus, transcription and HPV DNA integration in vulvar cancers resemble those in other HPV-positive cancers. This study emphasizes the power of hybridization capture coupled with DNA and RNA sequencing to identify a broad spectrum of HPV types, determine human genome integration status of viral DNAs, and elucidate their structures.

In vitro evolution and whole genome analysis to study chemotherapy drug resistance in haploid human cells

In vitro evolution and whole genome analysis has proven to be a powerful method for studying the mechanism of action of small molecules in many haploid microbes but has generally not been applied to human cell lines in part because their diploid state complicates the identification of variants that confer drug resistance. To determine if haploid human cells could be used in MOA studies, we evolved resistance to five different anticancer drugs (doxorubicin, gemcitabine, etoposide, topotecan, and paclitaxel) using a near-haploid cell line (HAP1) and then analyzed the genomes of the drug resistant clones, developing a bioinformatic pipeline that involved filtering for high frequency alleles predicted to change protein sequence, or alleles which appeared in the same gene for multiple independent selections with the same compound. Applying the filter to sequences from 28 drug resistant clones identified a set of 21 genes which was strongly enriched for known resistance genes or known drug targets (TOP1, TOP2A, DCK, WDR33, SLCO3A1). In addition, some lines carried structural variants that encompassed additional known resistance genes (ABCB1, WWOX and RRM1). Gene expression knockdown and knockout experiments of 10 validation targets showed a high degree of specificity and accuracy in our calls and demonstrates that the same drug resistance mechanisms found in diverse clinical samples can be evolved, discovered and studied in an isogenic background.

Whole genome resequencing unveils low-temperature stress tolerance specific genomic variations in jute (Corchorus sp.)

Jute (Corchorus sp.), a commercially important and eco-friendly crop, is widely cultivated in Bangladesh, India, and China. Some varieties of this tropical plant such as the Corchorus olitorius. Variety accession no. 2015 (acc. 2015) has been found to be low-temperature tolerant. The current study was designed to explore the genome-wide variations present in the tolerant plant acc. 2015 in comparison to the sensitive farmer popular variety Corchorus olitorius var. O9897 using the whole genome resequencing technique. Among different variations, intergenic Single Nucleotide Polymorphism (SNPs) and Insertion-Deletion (InDels) were found in the highest percentage whereas approximately 3% SNPs and 2% InDels were found in exonic regions in both plants. Gene enrichment analysis indicated the presence of acc. 2015 specific SNPs in the genes encoding peroxidase, ER lumen protein retaining receptor, and hexosyltransferase involved in stress response (GO:0006950) which were not present in sensitive variety O9897. Besides, distinctive copy number variation regions (CNVRs) comprising 120 gene loci were found in acc. 2015 with a gain of function from multiple copy numbers but absent in O9897. Gene ontology analysis revealed these gene loci to possess different receptors like kinases, helicases, phosphatases, transcription factors especially Myb transcription factors, regulatory proteins containing different binding domains, annexin, laccase, acyl carrier protein, potassium transporter, and vesicular transporter proteins that are responsible for low temperature induced adaptation pathways in plants. This work of identifying genomic variations linked to cold stress tolerance traits will help to develop successful markers that will pave the way to develop genetically modified cold-resistant jute lines for year-round cultivation to meet the demand for a sustainable fiber crop economy.

Unveiling the population genetic structure of Iranian horses breeds by whole-genome resequencing analysis

Preserving genetic diversity is pivotal for enhancing genetic improvement and facilitating adaptive responses to selection. This study focuses on identifying key genetic variants, including single nucleotide polymorphisms (SNPs), insertion/deletion polymorphisms (INDELs), and copy number variants (CNVs), while exploring the genomic evolutionary connectedness among seven Iranian horses representing five indigenous breeds: Caspian, Turkemen, DareShuri, Kurdish, and Asil. Using whole-genome resequencing, we generated 2.7 Gb of sequence data, with raw reads ranging from 1.2 Gb for Caspian horses to 0.38 Gb for Turkoman horses. Post-filtering, approximately 1.9 Gb of reads remained, with ~ 1.5 Gb successfully mapped to the horse reference genome (EquCab3.0), achieving mapping rates between 76.4% (Caspian) and 98.35% (Turkoman). We identified 2,909,816 SNPs in Caspian horses, constituting around 0.1% of the genome. Notably, 71% of these SNPs were situated in intergenic regions, while 8.5 and 6.8% were located upstream and downstream, respectively. A comparative analysis of SNPs between Iranian and non-Iranian horse breeds showed that Caspian horses had the lowest number of shared SNPs with Turkoman horses. Instead, they showed a closer genetic relationship with DareShuri, Quarter, Arabian, Standardbred, and Asil breeds. Hierarchical clustering highlighted Caspian horses as a distinct cluster, underscoring their distinctive genomic signature. Caspian horses exhibit a unique genetic profile marked by an enrichment of private mutations in neurological genes, influencing sensory perception and awareness. This distinct genetic makeup shapes mating preferences and signifies a separate evolutionary trajectory. Additionally, significant non-synonymous single nucleotide polymorphisms (nsSNPs) in reproductive genes offer intervention opportunities for managing Caspian horses. These findings reveal the population genetic structure of Iranian horse breeds, contributing to the advancement of knowledge in areas such as conservation, performance traits, climate adaptation, reproduction, and resistance to diseases in equine science.

Copy number deletion of PLA2G4A affects the susceptibility and clinical phenotypes of schizophrenia

Phospholipase A2(PLA2) superfamily is recognized as being involved in the pathogenesis of schizophrenia by affecting lipid homeostasis in cell membranes. We hypothesized that PLA2 gene copy number variation (CNV) may affect PLA2 enzyme expression and be associated with schizophrenia risk. This study indicated that in the discovery stage, an increased copy number of PLA2G6 and the deletion of PLA2G3, PLA2G4A, PLA2G4F and PLA2G12F was associated with increased risk of schizophrenia. CNV segments involving six PLA2 genes were detected in publicly available datasets, including two deletion segments specific to the PLA2G4A gene. The relationship between the deletion of PLA2G4A and susceptibility to schizophrenia was then reaffirmed in the validation group of 806 individuals. There was a significant correlation between PLA2G4A deletion and the symptoms of poverty of thought in male patients and erotomanic delusion in females. Furthermore, ELISA results demonstrate a significant decrease in peripheral blood cytosolic PLA2(cPLA2) levels in patients with the PLA2G4A deletion genotype compared to those with normal and copy number duplicate genotypes. These data suggest that the functional copy number deletion in the PLA2G4A gene is associated with the risk of schizophrenia and clinical phenotypes by reducing the expression of cPLA2, which may be an indicator of susceptibility to schizophrenia.

Age-specific induction of mutant p53 drives clonal hematopoiesis and acute myeloid leukemia in adult mice

The investigation of the mechanisms behind p53 mutations in acute myeloid leukemia (AML) has been limited by the lack of suitable mouse models, which historically have resulted in lymphoma rather than leukemia. This study introduces two new AML mouse models. One model induces mutant p53 and Mdm2 haploinsufficiency in early development, showing the role of Mdm2 in myeloid-biased hematopoiesis and AML predisposition, independent of p53. The second model mimics clonal hematopoiesis by inducing mutant p53 in adult hematopoietic stem cells, demonstrating that the timing of p53 mutation determines AML vs. lymphoma development. In this context, age-related changes in hematopoietic stem cells (HSCs) collaborate with mutant p53 to predispose toward myeloid transformation rather than lymphoma development. Our study unveils new insights into the cooperative impact of HSC age, Trp53 mutations, and Mdm2 haploinsufficiency on clonal hematopoiesis and the development of myeloid malignancies.

HATCHet2: clone- and haplotype-specific copy number inference from bulk tumor sequencing data

Bulk DNA sequencing of multiple samples from the same tumor is becoming common, yet most methods to infer copy-number aberrations (CNAs) from this data analyze individual samples independently. We introduce HATCHet2, an algorithm to identify haplotype- and clone-specific CNAs simultaneously from multiple bulk samples. HATCHet2 extends the earlier HATCHet method by improving identification of focal CNAs and introducing a novel statistic, the minor haplotype B-allele frequency (mhBAF), that enables identification of mirrored-subclonal CNAs. We demonstrate HATCHet2's improved accuracy using simulations and a single-cell sequencing dataset. HATCHet2 analysis of 10 prostate cancer patients reveals previously unreported mirrored-subclonal CNAs affecting cancer genes.

Whole genome and transcriptome integrated analyses guide clinical care of pediatric poor prognosis cancers

The role for routine whole genome and transcriptome analysis (WGTA) for poor prognosis pediatric cancers remains undetermined. Here, we characterize somatic mutations, structural rearrangements, copy number variants, gene expression, immuno-profiles and germline cancer predisposition variants in children and adolescents with relapsed, refractory or poor prognosis malignancies who underwent somatic WGTA and matched germline sequencing. Seventy-nine participants with a median age at enrollment of 8.8 y (range 6 months to 21.2 y) are included. Germline pathogenic/likely pathogenic variants are identified in 12% of participants, of which 60% were not known prior. Therapeutically actionable variants are identified by targeted gene report and whole genome in 32% and 62% of participants, respectively, and increase to 96% after integrating transcriptome analyses. Thirty-two molecularly informed therapies are pursued in 28 participants with 54% achieving a clinical benefit rate; objective response or stable disease ≥6 months. Integrated WGTA identifies therapeutically actionable variants in almost all tumors and are directly translatable to clinical care of children with poor prognosis cancers.

Longitudinal molecular profiling elucidates immunometabolism dynamics in breast cancer

Although metabolic reprogramming within tumor cells and tumor microenvironment (TME) is well described in breast cancer, little is known about how the interplay of immune state and cancer metabolism evolves during treatment. Here, we characterize the immunometabolic profiles of tumor tissue samples longitudinally collected from individuals with breast cancer before, during and after neoadjuvant chemotherapy (NAC) using proteomics, genomics and histopathology. We show that the pre-, on-treatment and dynamic changes of the immune state, tumor metabolic proteins and tumor cell gene expression profiling-based metabolic phenotype are associated with treatment response. Single-cell/nucleus RNA sequencing revealed distinct tumor and immune cell states in metabolism between cold and hot tumors. Potential drivers of NAC based on above analyses were validated in vitro. In summary, the study shows that the interaction of tumor-intrinsic metabolic states and TME is associated with treatment outcome, supporting the concept of targeting tumor metabolism for immunoregulation.

Genome-wide somatic mutation analysis of sinonasal adenocarcinoma with and without wood dust exposure

BACKGROUND

Sinonasal adenocarcinoma is a rare cancer, encompassing two different entities, the intestinal-type sinonasal adenocarcinoma (ITAC) and the non-intestinal-type sinonasal adenocarcinoma (non-ITAC). Occurrence of ITAC is strongly associated with exposure to hardwood dusts. In countries with predominant exposure to softwood dust the occurrence of sinonasal adenocarcinomas is lower and the relative amount of non-ITACs to ITACs is higher. The molecular mechanisms behind the tumorigenic effects of wood dust remain largely unknown.

METHODS

We carried out whole-genome sequencing of formalin-fixed paraffin-embedded (FFPE) samples of sinonasal adenocarcinomas from ten wood dust-exposed and six non-exposed individuals, with partial tobacco exposure data. Sequences were analyzed for the presence of mutational signatures matching COSMIC database signatures. Driver mutations and CN variant regions were characterized.

RESULTS

Mutation burden was higher in samples of wood dust-exposed patients (p = 0.016). Reactive oxygen species (ROS) damage-related mutational signatures were almost exclusively identified in ITAC subtype samples (p = 0.00055). Tobacco smoke mutational signatures were observed in samples of patients with tobacco exposure or missing information, but not in samples from non-exposed patients. A tetraploidy copy number (CN) signature was enriched in ITAC subtype (p = 0.042). CN variation included recurrent gains in COSMIC Cancer Gene Census genes TERT, SDHA, RAC1, ETV1, PCM1, and MYC. Pathogenic variants were observed most frequently in TP53, NF1, CHD2, BRAF, APC, and LRP1B. Driver mutations and copy number gains did not segregate by subtype.

CONCLUSIONS

Our analysis identified distinct mutational characteristics in ITAC and non-ITAC. Mutational signature analysis may eventually become useful for documentation of occupation-related cancer, while the exact mechanisms behind wood dust-driven carcinogenesis remain elusive. The presence of homologous recombination deficiency signatures implies a novel opportunity for treatment, but further studies are needed.

Genomic variation in weedy and cultivated broomcorn millet accessions uncovers the genetic architecture of agronomic traits

Large-scale genomic variations are fundamental resources for crop genetics and breeding. Here we sequenced 1,904 genomes of broomcorn millet to an average of 40× sequencing depth and constructed a comprehensive variation map of weedy and cultivated accessions. Being one of the oldest cultivated crops, broomcorn millet has extremely low nucleotide diversity and remarkably rapid decay of linkage disequilibrium. Genome-wide association studies identified 186 loci for 12 agronomic traits. Many causative candidate genes, such as PmGW8 for grain size and PmLG1 for panicle shape, showed strong selection signatures during domestication. Weedy accessions contained many beneficial variations for the grain traits that are largely lost in cultivated accessions. Weedy and cultivated broomcorn millet have adopted different loci controlling flowering time for regional adaptation in parallel. Our study uncovers the unique population genomic features of broomcorn millet and provides an agronomically important resource for cereal crops.

Intergenic sequences harboring potential enhancer elements contribute to Axenfeld-Rieger syndrome by regulating PITX2

Recent studies have uncovered that noncoding sequence variants may relate to Axenfeld-Rieger syndrome (ARS), a rare developmental anomaly with genetic heterogeneity. However, how these genomic regions are functionally and structurally associated with ARS is still unclear. In this study, we performed genome-wide linkage analysis and whole-genome sequencing in a Chinese family with ARS and identified a heterozygous deletion of about 570 kb (termed LOH-1) in the intergenic sequence between paired-like homeodomain transcription factor 2 (PITX2) and family with sequence similarity 241 member A. Knockout of LOH-1 homologous sequences caused ARS phenotypes in mice. RNA-Seq and real-time quantitative PCR revealed a significant reduction in Pitx2 gene expression in LOH-1-/- mice, while forkhead box C1 expression remained unchanged. ChIP-Seq and bioinformatics analysis identified a potential enhancer region (LOH-E1) within LOH-1. Deletion of LOH-E1 led to a substantial downregulation of the PITX2 gene. Mechanistically, we found a sequence (hg38 chr4:111,399,594-111,399,691) that is on LOH-E1 could regulate PITX2 by binding to RAD21, a critical component of the cohesin complex. Knockdown of RAD21 resulted in reduced PITX2 expression. Collectively, our findings indicate that a potential enhancer sequence that is within LOH-1 may regulate PITX2 expression remotely through cohesin-mediated loop domains, leading to ARS when absent.

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.

Chromophobe Renal Cell Carcinoma With Extensive Retraction Artifact: A Potential Diagnostic Pitfall From Micropapillary Urothelial Carcinoma

In addition to "classic" and eosinophilic subtype, chromophobe renal cell carcinoma (RCC) is well-known to demonstrate various morphological patterns including adenomatoid, microcystic, pigmented, multicystic, papillary, neuroendocrine-like, and small cell-like, all of which are important to appreciate for accurate diagnosis. Herein, we expand on a unique chromophobe RCC morphology not previously described consisting of tumor cells with extensive stromal retraction, mimicking upper urothelial tract micropapillary carcinoma (MPC). Twelve MPC-like chromophobe RCC nephrectomies were reviewed with clinicopathological features recorded; molecular testing was performed on 7 of 12 tumors. Patients were mostly men (n=10) with a mean age of 65 years. Mean tumor size was 6.4 cm with pathological stage distribution as follows: 4 (33%) T1a, 2 (17%) T1b, 1 (8%) T2b, and 3 (25%) T3a. The extent of MPC-like chromophobe RCC foci ranged from 10% to 40% (mean=26%; there was no correlation between the extent of MPC-like chromophobe RCC foci and tumor stage). Other chromophobe RCC morphological patterns were not identified. When performed, all (100%) tumors depicted prototypic chromophobe RCC staining pattern of KIT positivity/KRT7 positivity. Molecular showed 6 of 7 (86%) with multiple chromosomal losses. Clinically significant mutations were identified in NF1, TP53, FLCN (likely somatic), CHEK2, and ZFHX3 genes. Follow up available in 9 patients showed no evidence of disease (mean=23 months). Although the etiology behind the extensive stromal retraction in our tumors is unknown, this may likely be artifactual in nature. Nonetheless, it is important to include MPC-like chromophobe RCC in the spectrum of "variant" morphologies to avoid diagnostic pitfalls from micropapillary carcinoma.

Ancient and recent origins of shared polymorphisms in yeast

Shared genetic polymorphisms between populations and species can be ascribed to ancestral variation or to more recent gene flow. Here, we mapped shared polymorphisms in Saccharomyces cerevisiae and its sister species Saccharomyces paradoxus, which diverged 4-6 million years ago. We used a dense map of single-nucleotide diagnostic markers (mean distance 15.6 base pairs) in 1,673 sequenced S. cerevisiae isolates to catalogue 3,852 sequence blocks (≥5 consecutive markers) introgressed from S. paradoxus, with most being recent and clade-specific. The highly diverged wild Chinese S. cerevisiae lineages were depleted of introgressed blocks but retained an excess of individual ancestral polymorphisms derived from incomplete lineage sorting, perhaps due to less dramatic population bottlenecks. In the non-Chinese S. cerevisiae lineages, we inferred major hybridization events and detected cases of overlapping introgressed blocks across distinct clades due to either shared histories or convergent evolution. We experimentally engineered, in otherwise isogenic backgrounds, the introgressed PAD1-FDC1 gene pair that independently arose in two S. cerevisiae clades and revealed that it increases resistance against diverse antifungal drugs. Overall, our study retraces the histories of divergence and secondary contacts across S. cerevisiae and S. paradoxus populations and unveils a functional outcome.

DELE1 haploinsufficiency causes resistance to mitochondrial stress-induced apoptosis in monosomy 5/del(5q) AML

Monosomy 5 and deletions of the chromosome 5q (-5/del(5q)) are recurrent events in de novo adult acute myeloid leukemia (AML), reaching up to 40% of cases in secondary AML. These chromosome anomalies are associated with TP53 mutations and with very poor prognosis. Using the large Leucegene genomic and transcriptomic dataset composed of 48 -5/del(5q) patient specimens and 367 control AML, we identified DELE1 - located in the common deleted region - as the most consistently downregulated gene in these leukemias. DELE1 encodes a mitochondrial protein recently characterized as the relay of mitochondrial stress to the cytosol through a newly defined OMA1-DELE1-HRI pathway which ultimately leads to the activation of ATF4, the master transcription factor of the integrated stress response. Here, we showed that the partial loss of DELE1 expression observed in -5/del(5q) patients was sufficient to significantly reduce the sensitivity to mitochondrial stress in AML cells. Overall, our results suggest that DELE1 haploinsufficiency could represent a new driver mechanism in -5/del(5q) AML.

Mutational signature and prognosis in adenocarcinoma of the bladder

Adenocarcinoma of the bladder is a rare urinary bladder carcinoma with limited therapy options due to lack of molecular characterization. Here, we aimed to reveal the mutational and transcriptomic landscapes of adenocarcinoma of the bladder and assess any relationship with prognosis. Between February 2015 and June 2021, a total of 23 patients with adenocarcinoma of the bladder were enrolled. These included 16 patients with primary bladder adenocarcinomas and seven patients with urachal adenocarcinoma. Whole exome sequencing (16 patients), whole genome sequencing (16 patients), bulk RNA sequencing (RNA-seq) (19 patients), and single-cell RNA-seq (5 patients) were conducted for the specimens. Correlation analysis, survival analysis, and t-tests were also performed. Prevalent T>A substitutions were observed among somatic mutations, and major trinucleotide contexts included 5'-CTC-3' and 5'-CTG-3'. This pattern was mainly contributed by COSMIC signature 22 related to chemical carcinogen exposure (probably aristolochic acid), which has not been reported in bladder adenocarcinoma. Moreover, genes with copy number changes were also enriched in the KEGG term 'chemical carcinogenesis'. Transcriptomic analysis suggested high immune cell infiltration and luminal-like features in the majority of samples. Interestingly, a small fraction of samples with an APOBEC-derived mutational signature exhibited a higher risk of disease progression compared with samples with only a chemical carcinogen-related signature, confirming the molecular and prognostic heterogeneity of bladder adenocarcinoma. This study presents mutational and transcriptomic landscapes of bladder adenocarcinoma, and indicates that a chemical carcinogen-related mutational signature may be related to a better prognosis compared with an APOBEC signature in adenocarcinoma of the bladder. © 2024 The Pathological Society of Great Britain and Ireland.

Mechanisms of action and resistance in histone methylation-targeted therapy

Epigenomes enable the rectification of disordered cancer gene expression, thereby providing new targets for pharmacological interventions. The clinical utility of targeting histone H3 lysine trimethylation (H3K27me3) as an epigenetic hallmark has been demonstrated1-7. However, in actual therapeutic settings, the mechanism by which H3K27me3-targeting therapies exert their effects and the response of tumour cells remain unclear. Here we show the potency and mechanisms of action and resistance of the EZH1-EZH2 dual inhibitor valemetostat in clinical trials of patients with adult T cell leukaemia/lymphoma. Administration of valemetostat reduced tumour size and demonstrated durable clinical response in aggressive lymphomas with multiple genetic mutations. Integrative single-cell analyses showed that valemetostat abolishes the highly condensed chromatin structure formed by the plastic H3K27me3 and neutralizes multiple gene loci, including tumour suppressor genes. Nevertheless, subsequent long-term treatment encounters the emergence of resistant clones with reconstructed aggregate chromatin that closely resemble the pre-dose state. Acquired mutations at the PRC2-compound interface result in the propagation of clones with increased H3K27me3 expression. In patients free of PRC2 mutations, TET2 mutation or elevated DNMT3A expression causes similar chromatin recondensation through de novo DNA methylation in the H3K27me3-associated regions. We identified subpopulations with distinct metabolic and gene translation characteristics implicated in primary susceptibility until the acquisition of the heritable (epi)mutations. Targeting epigenetic drivers and chromatin homeostasis may provide opportunities for further sustained epigenetic cancer therapies.