CONTRA: copy number analysis for targeted resequencing

Genomic landscape defines peritoneal metastatic pattern and related target of peritoneal metastasis in colorectal cancer

The primary objective of this study is to develop a prediction model for peritoneal metastasis (PM) in colorectal cancer by integrating the genomic features of primary colorectal cancer, along with clinicopathological features. Concurrently, we aim to identify potential target implicated in the peritoneal dissemination of colorectal cancer through bioinformatics exploration and experimental validation. By analyzing the genomic landscape of primary colorectal cancer and clinicopathological features from 363 metastatic colorectal cancer patients, we identified 22 differently distributed variables, which were used for subsequent LASSO regression to construct a PM prediction model. The integrated model established by LASSO regression, which incorporated two clinicopathological variables and seven genomic variables, precisely discriminated PM cases (AUC 0.899; 95% CI 0.860-0.937) with good calibration (Hosmer-Lemeshow test p = .147). Model validation yielded AUCs of 0.898 (95% CI 0.896-0.899) and 0.704 (95% CI 0.622-0.787) internally and externally, respectively. Additionally, the peritoneal metastasis-related genomic signature (PGS), which was composed of the seven genes in the integrated model, has prognostic stratification capability for colorectal cancer. The divergent genomic landscape drives the driver genes of PM. Bioinformatic analysis concerning these driver genes indicated SERINC1 may be associated with PM. Subsequent experiments indicate that knocking down of SERINC1 functionally suppresses peritoneal dissemination, emphasizing its importance in CRCPM. In summary, the genomic landscape of primary cancer in colorectal cancer defines peritoneal metastatic pattern and reveals the potential target of SERINC1 for PM in colorectal cancer.

Molecular landscape of ERBB2 alterations in 3000 advanced NSCLC patients

ERBB2 (HER2) represents a newly recognized actionable oncogenic driver in non-small cell lung cancer (NSCLC), with approved targeted therapy available. Understanding the landscape of ERBB2 alterations and co-occurring mutations is essential for guiding treatment decisions. We conducted an analysis involving 3000 NSCLC patients with all types of ERBB2 alterations, drawn from two extensive retrospective cohorts: 1281 from Geneplus (Chinese) and 1719 from Guardant360 (the United States, US). The incidence of all types of ERBB2 alterations was found to be 5.6% in the Chinese group and 5.2% in the US group. In both cohorts, among oncogenic alterations of ERBB2, exon 20 insertion Y772_A775dupYVMA was the most frequent alteration (58% vs 41.6% in the Chinese vs the US), followed by G776delinsVC/LC/VV/IC (10.7% vs 9.7%), and S310X (10.5% vs 15.4%). EGFR ex20 insertions were identified in the A767-V774 region, whereas ERBB2 ex20 insertions were observed in the Y772-P780 region. Notably, EGFR ex20 insertions exhibited greater insertion diversity. Clinical characteristics of EGFR and ERBB2 ex20 NSCLC were similar, characterized by low tumor mutation burden (TMB), a predominant never-smoker population, and a majority of lung adenocarcinoma cases.

Importance of IFT140 in Patients with Polycystic Kidney Disease Without a Family History

Introduction

Recently, the monoallelic loss-of-function IFT140 variant was identified as a causative gene for autosomal dominant polycystic kidney disease (ADPKD). In patients with polycystic kidneys who have a positive family history, >90% have pathogenic variants in PKD1 or PKD2, whereas only 1% have IFT140. However, approximately 40% of patients with polycystic kidneys without a family history do not have any pathogenic variants in PKD1 and PKD2.

Methods

We conducted a comprehensive genetic analysis of 157 adult patients with polycystic kidneys whose parents did not have evident polycystic kidneys. We sequenced up to 92 genes associated with inherited cystic kidney disease, including IFT140.

Results

Of the 157 patients, 7 (4.5%) presented with monoallelic loss-of-function variants in the IFT140 gene, 51 (32.5%) with pathogenic variants in the PKD1 or PKD2 gene, and 7 (4.5%) with pathogenic variants in other genes related to inherited kidney cystic disease. The proportion of monoallelic loss-of-function IFT140 variants in this cohort was higher than that in previously reported cohorts with polycystic kidneys who had a positive family history. None of the patients with monoallelic loss-of-function IFT140 variants had polycystic liver disease (PLD). Furthermore, patients with IFT140 pathogenic variants had a significantly smaller kidney volume and a remarkably higher estimated glomerular filtration rate (eGFR) than those with PKD1 pathogenic variants (P = 0.01 and 0.03, respectively).

Conclusion

Because the phenotype of polycystic kidneys caused by the IFT140 gene is mild, parental kidney disease may be overlooked. Therefore, patients without a positive family history are more likely to carry pathogenic variants in IFT140.

Evaluation of molecular residual disease in operable non-small cell lung cancer with gene fusions, MET exon skipping or de novo MET amplification

Gene fusions and MET alterations are rare and difficult to detect in plasma samples. The clinical detection efficacy of molecular residual disease (MRD) based on circulating tumor DNA (ctDNA) in patients with non-small cell lung cancer (NSCLC) with these mutations remains unknown. This prospective, non-intervention study recruited 49 patients with operable NSCLC with actionable gene fusions (ALK, ROS1, RET, and FGFR1), MET exon 14 skipping or de novo MET amplification. We analyzed 43 tumor tissues and 111 serial perioperative plasma samples using 1021- and 338-gene panels, respectively. Detectable MRD correlated with a significantly higher recurrence rate (P < 0.001), yielding positive predictive values of 100% and 90.9%, and negative predictive values of 82.4% and 86.4% at landmark and longitudinal time points, respectively. Patients with detectable MRD showed reduced disease-free survival (DFS) compared to those with undetectable MRD (P < 0.001). Patients who harbored tissue-derived fusion/MET alterations in their MRD had reduced DFS compared to those who did not (P = 0.05). To our knowledge, this is the first comprehensive study on ctDNA-MRD clinical detection efficacy in operable NSCLC patients with gene fusions and MET alterations. Patients with detectable tissue-derived fusion/MET alterations in postoperative MRD had worse clinical outcomes.

Clinical, genomic and immune microenvironmental determinants of nivolumab response in head and neck squamous cell carcinoma

Background

In view of improving biomarkers predicting the efficacy of immunotherapy for head and neck squamous cell carcinoma (R/M HNSCC), this multicenter retrospective study aimed to identify clinical, tumor microenvironmental, and genomic factors that are related to therapeutic response to the anti- Programmed cell death protein 1 (PD-1) antibody, nivolumab, in patients with R/M HNSCC.

Methods

The study compared 53 responders and 47 non-responders, analyzing formalin-fixed paraffin-embedded samples using 14-marker multiplex immunohistochemistry and targeted gene sequencing.

Results

Of 100 patients included, responders had significantly lower smoking and alcohol index, higher incidence of immune related adverse events, and higher PD-1 ligand (PD-L1) expression in immune cells as well as PD-L1 combined positive score (CPS) than non-responders. The frequency of natural killer cells was associated with nivolumab response in patients with prior cetuximab use, but not in cetuximab-naïve status. Age-stratified analysis showed nivolumab response was linked to high CPS and lymphoid-inflamed profiles in patients aged ≥ 65. In contrast, lower NLR in peripheral blood counts was associated with response in patients aged < 65. Notably, TP53 mutation-positive group had lower CPS and T cell densities, suggesting an immune-excluded microenvironment. Patients with altered tumor suppressor gene pathways, including TP53, CDKN2A, and SMAD4 mutations, had lower CPS, higher smoking index, and were associated with poor responses.

Conclusion

Nivolumab treatment efficacy in HNSCC is influenced by a combination of clinical factors, age, prior treatment, immune environmental characteristics, and gene mutation profiles.

Prognostic value of structural variants in early breast cancer patients

Genomic analysis of structural variants(SVs) in breast cancer (BC) patients has been conducted, but the relationship between genomic alterations and BC prognosis remains unclear. We performed RNA sequencing of 297 early BC fresh-frozen tissues. We identified SVs using three tools (STAR.Arriba, STAR.fusion, and STAR.SEQR) with the COSMIC and Mitelman databases as guide references. We found a median of five to eight fusions per sample. In BC intrinsic subtypes, normal subtype had the fewest fusions (median: 1, interquartile range [IQR]: 0, 3) followed by luminal A (median: 5.5, IQR: 2.75, 10.25), luminal B (median: 9, IQR: 6, 16.5), HER2-enriched (median: 9, IQR: 6, 16.5) and basal (median 10, IQR: 6, 15.5) subtypes (p < 0.05). Intrachromosomal fusion was more frequent observed rather than interchromosomal fusion. In location, chromosome 17 had the most fusions followed by chromosome 1 and 11. When samples were divided into high and low fusion groups based on a cut-off value of 11 fusions, five-year event-free survival (5Y-EFS) was 68.1% in the high fusion group (n = 72) and 80.1% in the low fusion group (n = 125) (p = 0.024) while 75.6% among all patients (95% confidence interval: 0.699, 0.819). Among BC subtype, TNBCs with more fusions had shorter EFS compared to those with fewer fusions (5Y-EFS rate: 65.1% vs. 85.7%; p = 0.013) but no EFS differences were observed in other BC subtypes. ESTIMATE ImmuneScore was also associated with the number of fusions in TNBC (p < 0.005) and TNBCs with high ImmuneScore had better 5Y-EFS compared to those with low ImmuneScore (p = 0.041). In conclusion, diverse fusions were observed by BC subtype, and the number of fusions was associated with BC survival outcome and immune status in TNBC.

Baseline circulating tumor DNA predicts long-term survival outcomes for patients with early breast cancer

Background

Circulating tumor DNA (ctDNA) is a potential biomarker not only capable of monitoring the treatment response during neoadjuvant therapy (NAT) or rescue therapy, but also identifying minimal residual disease (MRD) and detecting early relapses after primary treatment. However, it remains uncertain whether the detection of ctDNA at diagnosis, before any treatment, can predict the prognosis for patients with early breast cancer. The objective of our study was to evaluate the predictive value of baseline ctDNA for prognosis in patients with early breast cancer.

Methods

A total of 90 patients with early breast cancer and 24 healthy women were recruited between August 2016 and October 2016. Peripheral blood samples were collected from patients at diagnosis, before any treatment. Blood samples were processed and subjected to targeted deep sequencing with a next-generation sequencing (NGS) panel of 1,021 cancer-related genes. The recurrence-free survival (RFS) and invasive disease-free survival (iDFS) were reported.

Results

The 90 patients with breast cancer included 6 patients with ductal carcinoma in situ (DCIS) and 84 patients with invasive breast cancer. Within the cohort of patients with invasive breast cancer, ctDNA were detected in 57 patients, with a ctDNA detection rate of 67.9%. Meanwhile, no ctDNA was detected in DCIS patients. Among 84 patients with invasive breast cancer, patients with high-level ctDNA had a significantly lower RFS compared to patients with low-level ctDNA (log-rank P=0.0036).

Conclusions

Our study suggested that ctDNA at diagnosis, before any treatment, could potentially serve as a biomarker to predict the prognosis for patients with early breast cancer. However, further follow-up and more studies with large sample sizes are required to confirm these findings.

Intrahepatic cholangiocarcinoma with FGFR alterations: A series of Chinese cases with an emphasis on their clinicopathologic and genetic features

Intrahepatic Cholangiocarcinoma (iCCA) with FGFR alterations is relatively rare, and its identification is important in the era of targeted therapy. We collected a large series of FGFR-altered cases in the Chinese population and characterized their clinicopathological and genetic features. Among the 18 FGFR-altered cases out of 260 iCCAs, 10 were males and 8 were females, ranging in age from 35 to 74 years (mean, 57.3 years; median, 58 years). Pathologically, they include 9 cases of large duct (LD, 50 %) and small duct (SD, 50 %) types each. All of them (100 %, 18/18) showed microsatellite stable (MSS) and low tumor mutation burden (TMB). Genetically, FGFR alterations involved FGFR1 (20 %), FGFR2 (70 %), and FGFR3 (10 %), with FGFR2 rearrangement accounting for the most (11/18). The most frequently altered genes/biological processes were development/proliferation-related pathways (44 %), chromatin organization (20 %), and tumor suppressors (32 %). Our study further revealed the clinicopathological and genetic features of FGFR-altered iCCA and demonstrated that its occurrence may show regional or ethnic variability and is less common in the Chinese population. A significant number of LD-type iCCA cases also have FGFR alterations rather than the SD type.

Molecular patterns of resistance to immune checkpoint blockade in melanoma

Immune checkpoint blockade (ICB) has improved outcome for patients with metastatic melanoma but not all benefit from treatment. Several immune- and tumor intrinsic features are associated with clinical response at baseline. However, we need to further understand the molecular changes occurring during development of ICB resistance. Here, we collect biopsies from a cohort of 44 patients with melanoma after progression on anti-CTLA4 or anti-PD1 monotherapy. Genetic alterations of antigen presentation and interferon gamma signaling pathways are observed in approximately 25% of ICB resistant cases. Anti-CTLA4 resistant lesions have a sustained immune response, including immune-regulatory features, as suggested by multiplex spatial and T cell receptor (TCR) clonality analyses. One anti-PD1 resistant lesion harbors a distinct immune cell niche, however, anti-PD1 resistant tumors are generally immune poor with non-expanded TCR clones. Such immune poor microenvironments are associated with melanoma cells having a de-differentiated phenotype lacking expression of MHC-I molecules. In addition, anti-PD1 resistant tumors have reduced fractions of PD1+ CD8+ T cells as compared to ICB naïve metastases. Collectively, these data show the complexity of ICB resistance and highlight differences between anti-CTLA4 and anti-PD1 resistance that may underlie differential clinical outcomes of therapy sequence and combination.

CFAP47 is a novel causative gene implicated in X-linked polycystic kidney disease

Autosomal dominant polycystic kidney disease (ADPKD) is a well-described condition in which ~80% of cases have a genetic explanation, while the genetic basis of sporadic cystic kidney disease in adults remains unclear in ~30% of cases. This study aimed to identify novel genes associated with polycystic kidney disease (PKD) in patients with sporadic cystic kidney disease in which a clear genetic change was not identified in established genes. A next-generation sequencing panel analyzed known genes related to renal cysts in 118 sporadic cases, followed by whole-genome sequencing on 47 unrelated individuals without identified candidate variants. Three male patients were found to have rare missense variants in the X-linked gene Cilia And Flagella Associated Protein 47 (CFAP47). CFAP47 was expressed in primary cilia of human renal tubules, and knockout mice exhibited vacuolation of tubular cells and tubular dilation, providing evidence that CFAP47 is a causative gene involved in cyst formation. This discovery of CFAP47 as a newly identified gene associated with PKD, displaying X-linked inheritance, emphasizes the need for further cases to understand the role of CFAP47 in PKD.

Mobocertinib in Patients with EGFR Exon 20 Insertion-Positive Non-Small Cell Lung Cancer (MOON): An International Real-World Safety and Efficacy Analysis

EGFR exon 20 (EGFR Ex20) insertion mutations in non-small cell lung cancer (NSCLC) are insensitive to traditional EGFR tyrosine kinase inhibitors (TKIs). Mobocertinib is the only approved TKI specifically designed to target EGFR Ex20. We performed an international, real-world safety and efficacy analysis on patients with EGFR Ex20-positive NSCLC enrolled in a mobocertinib early access program. We explored the mechanisms of resistance by analyzing postprogression biopsies, as well as cross-resistance to amivantamab. Data from 86 patients with a median age of 67 years and a median of two prior lines of treatment were analyzed. Treatment-related adverse events (TRAEs) occurred in 95% of patients. Grade ≥3 TRAEs were reported in 38% of patients and included diarrhea (22%) and rash (8%). In 17% of patients, therapy was permanently discontinued, and two patients died due to TRAEs. Women were seven times more likely to discontinue treatment than men. In the overall cohort, the objective response rate to mobocertinib was 34% (95% CI, 24-45). The response rate in treatment-naïve patients was 27% (95% CI, 8-58). The median progression-free and overall survival was 5 months (95% CI, 3.5-6.5) and 12 months (95% CI, 6.8-17.2), respectively. The intracranial response rate was limited (13%), and one-third of disease progression cases involved the brain. Mobocertinib also showed antitumor activity following EGFR Ex20-specific therapy and vice versa. Potential mechanisms of resistance to mobocertinib included amplifications in MET, PIK3CA, and NRAS. Mobocertinib demonstrated meaningful efficacy in a real-world setting but was associated with considerable gastrointestinal and cutaneous toxicity.

Genetic Diagnosis of Adult Hemodialysis Patients With Unknown Etiology

Introduction

Kidney disease of unknown etiology accounts for 1 in 10 adult end-stage renal disease (ESRD) cases worldwide. The aim of this study is to clarify the genetic background of patients with chronic kidney disease (CKD) of unknown etiology who initiated renal replacement therapy (RRT) in adulthood.

Methods

This is a multicenter cross-sectional cohort study. Of the 1164 patients who attended 4 dialysis clinics in Japan, we first selected patients who started RRT between the ages of 20 and 49 years. After excluding patients with apparent causes of CKD (e.g., diabetic nephropathy, polycystic kidney disease (PKD) with family history, patients who underwent renal biopsy), 90 patients with CKD of unknown cause were included. The 298 genes associated with CKD were analyzed using capture-based targeted next-generation sequencing.

Results

Of the 90 patients, 10 (11.1%) had pathogenic variants in CKD-causing genes and 17 (18.9%) had variant of unknown significance (VUS). Three patients had PKD1 pathogenic variants, and 1 patient had PKD1 and COL4A4 pathogenic variants. In addition, 2 patients were diagnosed with atypical hemolytic uremic syndrome (aHUS) due to C3 or CFHR5. One patient each was diagnosed with Alport syndrome due to COL4A4 and COL4A3 variants, nephronophthisis due to NPHP1 variants, Fabry disease due to GLA variants, and autosomal-dominant tubulointerstitial kidney disease due to UMOD variants. Genetic diagnoses were not concordant with clinical diagnoses, except for patients with PKD1 variant.

Conclusion

This largest study on genetic analysis in hemodialysis-dependent adults revealed the presence of undiagnosed inherited kidney diseases.

DNA damage response alterations in clear cell renal cell carcinoma: clinical, molecular, and prognostic implications

BACKGROUND

DNA damage repair (DDR) pathways modulate cancer risk, progression, and therapeutic responses. Nonetheless, the characteristics and significance of DDR alterations in clear cell renal cell carcinoma (ccRCC) remain undefined. This study aimed to explore the predictive role, molecular mechanism, and tumor immune profile of DDR genes in ccRCC.

METHODS

We prospectively sequenced 757 tumors and matched blood DNA samples from Chinese patients with ccRCC using next-generation sequencing (NGS) and analyzed data from 537 patients from The Cancer Genome Atlas (TCGA). A comprehensive analysis was performed.

RESULTS

Fifty-two percent of Chinese patients with ccRCC harbored DDR gene mutations and 57% of TCGA patients. The immunotherapy treatment prognosis of patients with DDR gene mutations was superior to that of patients without DDR gene mutations (p = 0.047). DDR gene mutations were associated with more gene mutations and a higher tumor mutation load (TMB, p < 0.001). Moreover, patients with DDR gene mutations have a distinct mutational signature compared with those with wild-type DDR. Furthermore, the DDR-mut group had elevated neoantigen load (including single-nucleotide variants (SNV) and indel neoantigen load, p = 0.037 and p = 0.002, respectively), TCR Shannon (p = 0.025), and neutrophils (p = 0.010). DDR gene mutations exhibited a distinct immune profile with significantly higher expression levels of TNFSF9, CD70, ICAM1, and indoleamine-2,3-dioxygenase (IDO) and lower expression levels of VTCN1 and IL12A.

CONCLUSIONS

Our data suggest that the detection of somatic mutations in DDR genes can predict the efficacy of immunotherapy in patients with ccRCC. Furthermore, we revealed the unique molecular and immune mechanisms underlying ccRCC with DDR gene mutations.

The mutational pattern of homologous recombination repair genes in urothelial carcinoma and its correlation with immunotherapeutic response

BACKGROUND

The mutational pattern of homologous recombination repair (HRR)-associated gene alterations in Chinese urothelial carcinoma (UC) necessitates comprehensive sequencing efforts, and the clinical implications of HRR gene mutations in UC remain to be elucidated.

MATERIALS AND METHODS

We delineated the mutational landscape of 343 Chinese UC patients from West China Hospital and 822 patients from The Cancer Genome Atlas (TCGA) using next-generation sequencing (NGS). Data from 182 metastatic UC patients from MSK-IMPACT cohort were used to assess the association between HRR mutations and immunotherapy efficacy. Comprehensive transcriptomic analysis was performed to explore the impact of HRR mutations on tumor immune microenvironment.

RESULTS

Among Chinese UC patients, 34% harbored HRR gene mutations, with BRCA2, ATM, BRCA1, CDK12, and RAD51C being the most prevalently mutated genes. Mutational signatures contributing to UC differed between patients with and without HRR mutations. Signature 22 for exposure to aristolochic acid was only observed in Chinese UC patients. The presence of HRR mutations was correlated with higher tumor mutational burden, neoantigen burden, and PD-L1 expression. Importantly, patients with HRR mutations exhibited significantly improved prognosis following immunotherapy compared to those without HRR mutations.

CONCLUSIONS

Our findings provide valuable insights into the genomic landscape of Chinese UC patients and underscore the molecular rationale for utilizing immunotherapy in UC patients with HRR mutations.

Dynamic analysis of circulating tumor DNA to predict the prognosis and monitor the treatment response of patients with metastatic triple-negative breast cancer: A prospective study

Background

Limited data are available on applying circulating tumor DNA (ctDNA) in metastatic triple-negative breast cancer (mTNBC) patients. Here, we investigated the value of ctDNA for predicting the prognosis and monitoring the treatment response in mTNBC patients.

Methods

We prospectively enrolled 70 Chinese patients with mTNBC who had progressed after ≤2 lines of chemotherapy and collected blood samples to extract ctDNA for 457-gene targeted panel sequencing.

Results

Patients with ctDNA+, defined by 12 prognosis-relevant mutated genes, had a shorter progression-free survival (PFS) than ctDNA- patients (5.16 months vs. 9.05 months, p=0.001), and ctDNA +was independently associated with a shorter PFS (HR, 95% CI: 2.67, 1.2-5.96; p=0.016) by multivariable analyses. Patients with a higher mutant-allele tumor heterogeneity (MATH) score (≥6.316) or a higher ctDNA fraction (ctDNA%≥0.05) had a significantly shorter PFS than patients with a lower MATH score (5.67 months vs.11.27 months, p=0.007) and patients with a lower ctDNA% (5.45 months vs. 12.17 months, p<0.001), respectively. Positive correlations with treatment response were observed for MATH score (R=0.24, p=0.014) and ctDNA% (R=0.3, p=0.002), but not the CEA, CA125, or CA153. Moreover, patients who remained ctDNA +during dynamic monitoring tended to have a shorter PFS than those who did not (3.90 months vs. 6.10 months, p=0.135).

Conclusions

ctDNA profiling provides insight into the mutational landscape of mTNBC and may reliably predict the prognosis and treatment response of mTNBC patients.

Funding

This work was supported by the National Natural Science Foundation of China (Grant No. 81902713), Natural Science Foundation of Shandong Province (Grant No. ZR2019LZL018), Breast Disease Research Fund of Shandong Provincial Medical Association (Grant No. YXH2020ZX066), the Start-up Fund of Shandong Cancer Hospital (Grant No. 2020-PYB10), Beijing Science and Technology Innovation Fund (Grant No. KC2021-ZZ-0010-1).

Dynamic circulating tumor DNA during chemoradiotherapy predicts clinical outcomes for locally advanced non-small cell lung cancer patients

The value of circulating tumor DNA (ctDNA) during chemoradiotherapy (CRT) remains unclear but is critical for detecting molecular residual disease (MRD). In this prospective study, we sequenced 761 blood samples from 139 patients with locally advanced non-small cell lung cancer treated with definitive radiation therapy (RT). ctDNA concentrations showed a significantly declining trend as CRT progressed at on-RT and after-RT time points versus baseline. Thirty-eight (27.3%) patients with early undetectable ctDNA at both on-RT (RT reached 40 Gy) and after-RT time points, indicating early response to CRT, had better survival outcomes for both with or without consolidation immune checkpoint inhibitors. Longitudinal undetectable MRD was found in 20.1% patients. The 2-year cancer-specific progression-free survival of these patients was 88.4%, corresponding to a potentially cured population. Further analysis revealed that pretreatment ctDNA variants serve as an essential MRD informed source. These data provide clinical insights for ctDNA-MRD detection.

Genetic analysis of fundic gland‑type gastric adenocarcinoma

This study aimed to analyze the molecular characteristics of gastric adenocarcinoma of the fundic-gland type (GAFG) and explore the possible mechanism of tumor development. Samples from 10 Chinese patients with GAFG were collected at the Peking University International Hospital and Liaocheng People's Hospital between January 2015 and March 2022. The nucleic acid sequence of Epstein Barr virus-encoded RNA (EBV-EBER) was detected by in situ hybridization. Genetic mutation information for GNAS, KRAS, NRAS, BRAF, PIK3CA, TP53, APC, CTNNB1, HER2, MLH1, MSH2, MSH6, and PMS2 was obtained by Next-Generation Sequencing, and the relevant literature was reviewed. A total of eight instances of missense mutations were detected, consisting of seven cases with GNAS mutations, two cases with KRAS mutations, and one case with a TP53 mutation. Additionally, two patients had simultaneous missense mutations in GNAS and KRAS. Nonsynonymous mutations in APC, CTNNB1, NRAS, BRAF, PIK3CA, HER2, MLH1, MSH2, MSH6, or PMS2 were not observed in any cases. In addition, all tumors were EBER-negative. GAFG exhibits diversity at the molecular level, and GNAS mutations are more common than KRAS mutations, TP53 mutations, and microsatellite instability. To date, no association between EBV/HER2 and GAFG has been found.

Immune-active tumor-adjacent tissues are associated with favorable prognosis in stage I lung squamous cell carcinoma

The immunogenomic features of tumor-adjacent lungs (TALs) in stage I lung squamous cell carcinoma (LUSC) are not clear. Multiomics analyses of tumor tissues and paired TALs from 59 stage I LUSC patients were performed. Compared to tumors, TALs exhibited a better-preserved immune contexture indicated by upregulation of immune pathways, increased immune infiltration, and higher expression of immune effector molecules. Notably, TALs had no mutations in PTEN and KEAP1, a lower incidence of human leukocyte antigen (HLA) loss and higher expression of HLA class I genes, major histocompatibility complex (MHC) I chaperones, and interferon (IFN)-γ-associated genes. Digital spatial profiling validated the generally higher immune infiltration in TALs and revealed a higher level of immune heterogeneity in LUSC tumors. Importantly, patients with higher immune infiltration in TALs had significantly longer survival, while high immune heterogeneity was associated with inferior patient survival. Our work can be considered in the selection of patients for adjuvant therapy, especially immunotherapy.

Identification of a second genetic alteration in patients with SHOX deficiency individuals: a potential explanation for phenotype variability

OBJECTIVE

Our study aimed to assess the impact of genetic modifiers on the significant variation in phenotype that is observed in individuals with SHOX deficiency, which is the most prevalent monogenic cause of short stature.

DESIGN AND METHODS

We performed a genetic analysis in 98 individuals from 48 families with SHOX deficiency with a target panel designed to capture the entire SHOX genomic region and 114 other genes that modulate growth and/or SHOX action. We prioritized rare potentially deleterious variants.

RESULTS

We did not identify potential deleterious variants in the promoter or intronic regions of the SHOX genomic locus. In contrast, we found eight heterozygous variants in 11 individuals from nine families in genes with a potential role as genetic modifiers. In addition to a previously described likely pathogenic (LP) variant in CYP26C1 observed in two families, we identified LP variants in PTHLH and ACAN, and variants of uncertain significance in NPR2, RUNX2, and TP53 in more affected individuals from families with SHOX deficiency. Families with a SHOX alteration restricted to the regulatory region had a higher prevalence of a second likely pathogenic variant (27%) than families with an alteration compromising the SHOX coding region (2.9%, P = .04).

CONCLUSION

In conclusion, variants in genes related to the growth plate have a potential role as genetic modifiers of the phenotype in individuals with SHOX deficiency. In individuals with SHOX alterations restricted to the regulatory region, a second alteration could be critical to determine the penetrance and expression of the phenotype.

Mutational profiling of Chinese patients with thyroid cancer

Background

The incidence of thyroid cancer in China has rapidly increased in recent decades. As the genetic profiles of thyroid cancer vary dramatically between different geographical regions, a comprehensive genetic landscape of thyroid cancer in the Chinese population is urgently needed.

Methods

We retrospectively included thyroid cancer patients from three Chinese medical centers between February 2015 and August 2020. To dissect the genomic profiling of these patients, we performed targeted next-generation sequencing on their tumor tissues using a 1,021-gene panel.

Results

A total of 458 Chinese patients with thyroid cancer were enrolled, including four malignant histological subtypes arising from follicular epithelial thyroid cells. BRAF driver mutations were identified in 76.0% of patients, followed by RET rearrangements (7.6%) and RAS driver mutations (4.1%). Tumors with more somatic mutations correlated with worse clinical characteristics, including older age at diagnosis, less differentiation of tumor, larger tumor size, lymph node metastasis and distal metastasis. Subclonal BRAF mutations occurred in 20% (6/30) of patients and were frequent in poorly differentiated or anaplastic tumors (33.3% [2/6] vs. 4.2% [1/24], P = 0.09) and those with distal metastasis (50.0% [2/4] vs. 8.7% [2/23], P = 0.09). Tumors with TERT promoter mutations had significantly more somatic mutations (average: 6.5 vs. 1.8, P < 0.001). Moreover, TERT promoter mutations were not associated with lymph node metastasis but significantly associated with older age at diagnosis and poorly differentiated or anaplastic tumors, regardless of their clonal architecture.

Conclusion

Our results shed light on the molecular pathogenesis and clinical characteristics of thyroid cancer in the Chinese population. The number of somatic mutations, TERT promoter mutations, and the clonal architecture of BRAF mutations should be considered in the risk stratification of thyroid cancer.

Mutational pattern off homologous recombination repair (HRR)-related genes in upper tract urothelial carcinoma

BACKGROUND

Homologous recombination (HR) repair (HRR) has been indicated to be a biomarker for immunotherapy, chemotherapy, and poly-ADP ribose polymerase inhibitors inhibitors (PARPis). Nonetheless, their molecular correlates in upper tract urothelial carcinoma (UTUC) have not been well studied. This study aimed to explore the molecular mechanism and tumor immune profile of HRR genes and the relevance of their prognostic value in patients with UTUC.

MATERIALS AND METHODS

One hundred and ninety-seven tumors and matched blood samples from Chinese UTUC were subjected to next-generation sequencing. A total of 186 patients from The Cancer Genome Atlas were included. Comprehensive analysis was performed.

RESULTS

In Chinese patients with UTUC, 5.01% harbored germline HRR gene mutations, and 1.01% had Lynch syndrome-related genes. A total of 37.6% (74/197) of patients carried somatic or germline HRR gene mutations. There was marked discrepancy in the mutation landscapes, genetic interactions, and driver genes between the HRR-mut cohorts and HRR-wt cohorts. Aristolochic acid signatures and defective DNA mismatch repair signatures only existed in individuals in the HRR-mut cohorts. Inversely, the unknown signature (signature A) and signature SBS55 only existed in patients in the HRR-wt cohorts. HRR gene mutations regulated immune activities by NKT cells, plasmacytoid dendritic cells, hematopoietic stem cell, and M1 macrophages. In patients with local recurrence, patients with HRR gene mutations had poorer DFS rates than patients with wild-type HRR genes.

CONCLUSIONS

Our results imply that the detection of HRR gene mutations can predict recurrence in patients with UC. In addition, this study provides a path to explore the role of HRR-directed therapies, including PARPis, chemotherapy, and immunotherapy.

Genomic landscape and efficacy of HER2-targeted therapy in patients with HER2-mutant non-small cell lung cancer

Background

HER2-targeted therapy provides survival benefits to HER2-mutant non-small cell lung cancer (NSCLC). A better understanding of the clinical and genomic characterization of treatment-naïve HER2-positive NSCLC, as well as the efficacy of and resistance to HER2-targeted therapy in HER2-altered NSCLC, could promote further improvement of HER2 targeted therapy.

Methods

HER2-altered NSCLC patients was retrospectively included and their genomic profiles were performed by next-generation sequencing. The clinical outcomes included overall response rate, disease control rate and progression-free survival.

Results

Among 176 treatment-naïve patients with HER2 alterations, 64.8% harbored HER2 mutations with/without HER2 amplification, and 35.2% carried HER2 amplification only. Molecular characterization was correlated with tumor stage that late-stage NSCLC with HER2 oncogenic mutations showed a higher prevalence of TP53 mutations and a higher tumor mutation burden. However, this correlation was not found in patients with HER2 amplification only. Twenty-one patients with HER2 alterations treated with pyrotinib or afatinib were retrospectively enrolled. Pyrotinib yielded a longer median progression-free survival than afatinib (5.9 [95% CI, 3.8-13.0] vs. 4.0 months [95% CI, 1.9-6.3], P = 0.06) in these patients. Analysis of the genomic profiles before and after anti-HER2 targeted therapies identified de novo HER2 copy number gain and G518W mutation, as well as mutations involving DNA damage repair signaling, SWI–SNF complex, and epigenetic regulations as potential resistance mechanisms.

Conclusion

HER2-mutant NSCLC had different molecular features from HER2-amplified NSCLC, and its genomic profile was dependent of tumor stage. Pyrotinib had superior therapeutic effects than afatinib in HER2-altered NSCLC, although larger cohorts are warranted to validate it. HER2-dependent and -independent resistance mechanisms to afatinib and pyrotinib were unveiled.

Clinical outcome and biomarker assessments of a multi-centre phase II trial assessing niraparib with or without dostarlimab in recurrent endometrial carcinoma

This multi-centre, non-randomized, open-label, phase II trial (NCT03016338), assessed niraparib monotherapy (cohort 1, C1), or niraparib and dostarlimab (cohort 2, C2) in patients with recurrent serous or endometrioid endometrial carcinoma. The primary endpoint was clinical benefit rate (CBR), with ≥5/22 overall considered of interest. Secondary outcomes were safety, objective response rate (ORR), duration of response, progression free survival and overall survival. Translational research was an exploratory outcome. Potential biomarkers were evaluated in archival tissue by immunohistochemistry and next generation sequencing panel. In C1, 25 patients were enrolled, and CBR was 20% (95% CI: 9-39) with median clinical benefit duration of 5.3 months. The ORR was 4% (95% CI: 0-20). In C2, 22 patients were enrolled, and the CBR was 31.8% (95% CI: 16-53) with median clinical benefit duration of 6.8 months. The ORR was 14% (95% CI: 3-35). No new safety signals were detected. No significant association was detected between clinical benefit and IHC markers (PTEN, p53, MMR, PD-L1), or molecular profiling (PTEN, TP53, homologous recombination repair genes). In conclusion, niraparib monotherapy did not meet the efficacy threshold. Niraparib in combination with dostarlimab showed modest activity.

Impact of a cell cycle and an extracellular matrix remodeling transcriptional signature on tumor progression and correlation with EZH2 expression in meningioma

OBJECTIVE

The authors searched for genetic and transcriptional signatures associated with tumor progression and recurrence in their cohort of patients with meningiomas, combining the analysis of targeted exome, NF2-LOH, transcriptome, and protein expressions.

METHODS

The authors included 91 patients who underwent resection of intracranial meningioma at their institution between June 2000 and November 2007. The search of somatic mutations was performed by Next Generation Sequencing through a customized panel and multiplex ligation-dependent probe amplification for NF2 loss of heterozygosity. The transcriptomic profile was analyzed by QuantSeq 3' mRNA-Seq. The differentially expressed genes of interest were validated at the protein level analysis by immunohistochemistry.

RESULTS

The transcriptomic analysis identified an upregulated set of genes related to metabolism and cell cycle and downregulated genes related to immune response and extracellular matrix remodeling in grade 2 (atypical) meningiomas, with a significant difference in recurrent compared with nonrecurrent cases. EZH2 nuclear positivity associated with grade 2, particularly with recurrent tumors and EZH2 gene expression level, correlated positively with the expression of genes related to cell cycle and negatively to genes related to immune response and regulation of cell motility.

CONCLUSIONS

The authors identified modules of dysregulated genes in grade 2 meningiomas related to the activation of oxidative metabolism, cell division, cell motility due to extracellular remodeling, and immune evasion that were predictive of survival and exhibited significant correlations with EZH2 expression.

Next-generation sequencing shows the genomic features of ovarian clear cell cancer and compares the genetic architectures of high-grade serous ovarian cancer and clear cell carcinoma in ovarian and endometrial tissues

Ovarian clear cell carcinoma (OCCC) is a special histological type of epithelial ovarian cancer (EOC) that is not derived from epithelial cells of the ovarian or fallopian tube as the most common type of ovarian cancer, high-grade serous ovarian carcinoma (HGSOC), but is closely related to endometriosis and similar to endometrial clear cell carcinoma (ECCC) at morphologic and phenotypic features. However, limited data was shown in OCCC genomic features and compared with that in OCCC, HGSOC and ECCC. Herein, we utilized next-generation sequencing analysis of a panel of 1,021 genes to profile the mutational alterations in 34 OCCC and compared them to those from HGSOC (402 cases) and ECCC (30 cases). In result, the ARID1A and PIK3CA are high-frequency mutations of OCCC. Clonal architectures showed that all the mutations of genes occur in the later stage in the OCCC progress, whereas KRAS mutation is the earlier event compared with mutation of ARID1A or PIK3CA, which usually occurs in a group of ARID1A or PIK3CA mutations. The mutation frequency of main driver genes is similar between OCCC and ECCC, while TP53 is the main mutation in HGSOC and ECCC. Shared mutational signatures between OCCC and ECCC tissues with commonly observed a C>T change indicated a common carcinogens-exposed between these two carcinomas, but HGSOC and ECCC have common and distinct mutational signatures across cohorts respectively. In addition, we identified some novel CNV gains in NF1, ASXL1, TCF7L2, CREBBP and LRP1B and loss in ATM, FANCM, RB1 and FLT in OCCC. Our study offered a new perspective for OCCC tumorigenesis from two organs, the ovary and uterus, at genomic architectures and revealed novel CNV events for helping to provide theoretical support for OCCC treatment.

Genetic epidemiology of inherited retinal diseases in a large patient cohort followed at a single center in Italy

Inherited retinal diseases (IRDs) are the leading cause of vision loss in the working-age population. We performed a retrospective epidemiological study to determine the genetic basis of IRDs in a large Italian cohort (n = 2790) followed at a single referral center. We provided, mainly by next generation sequencing, potentially conclusive molecular diagnosis for 2036 patients (from 1683 unrelated families). We identified a total of 1319 causative sequence variations in 132 genes, including 353 novel variants, and 866 possibly actionable genotypes for therapeutic approaches. ABCA4 was the most frequently mutated gene (n = 535; 26.3% of solved cases), followed by USH2A (n = 228; 11.2%) and RPGR (n = 102; 5.01%). The other 129 genes had a lower contribution to IRD pathogenesis (e.g. CHM 3.5%, RHO 3.5%; MYO7A 3.4%; CRB1 2.7%; RPE65 2%, RP1 1.8%; GUCY2D 1.7%). Seventy-eight genes were mutated in five patients or less. Mitochondrial DNA variants were responsible for 2.1% of cases. Our analysis confirms the complex genetic etiology of IRDs and reveals the high prevalence of ABCA4 and USH2A mutations. This study also uncovers genetic associations with a spectrum of clinical subgroups and highlights a valuable number of cases potentially eligible for clinical trials and, ultimately, for molecular therapies.

Diagnostic yield of a multigene sequencing approach in children classified as idiopathic short stature

Objective

Most children with short stature remain without an etiologic diagnosis after extensive clinical and laboratory evaluation and are classified as idiopathic short stature (ISS). This study aimed to determine the diagnostic yield of a multigene analysis in children classified as ISS.

Design and methods

We selected 102 children with ISS and performed the genetic analysis as part of the initial investigation. We developed customized targeted panel sequencing, including all genes already implicated in the isolated short-stature phenotype. Rare and deleterious single nucleotide or copy number variants were assessed by bioinformatic tools.

Results

We identified 20 heterozygous pathogenic (P) or likely pathogenic (LP) genetic variants in 17 of 102 patients (diagnostic yield = 16.7%). Three patients had more than one P/LP genetic alteration. Most of the findings were in genes associated with the growth plate differentiation: IHH (n = 4), SHOX (n = 3), FGFR3 (n = 2), NPR2 (n = 2), ACAN (n = 2), and COL2A1 (n = 1) or involved in the RAS/MAPK pathway: NF1 (n = 2), PTPN11 (n = 1), CBL (n = 1), and BRAF (n = 1). None of these patients had clinical findings to guide a candidate gene approach. The diagnostic yield was higher among children with severe short stature (35% vs 12.2% for height SDS ≤ or > -3; P = 0.034). The genetic diagnosis had an impact on clinical management for four children.

Conclusion

A multigene sequencing approach can determine the genetic etiology of short stature in up to one in six children with ISS, removing the term idiopathic from their clinical classification.

Utility of cell-free DNA from bronchial washing fluid in diagnosis and genomic determination for radiology-suspected pulmonary nodules

BACKGROUND

Bronchial washing fluid (BWF) is a less-invasive specimen. Due to the limited sensitivity of BWF cellular component diagnosis, the aim of this study was to explore the potential role of BWF supernatant as a source of liquid biopsy of lung cancer.

METHODS

This prospective study enrolled 76 suspected and 5 progressed lung cancer patients. Transbronchial biopsy tissues, BWF supernatant (BWF_Sup) and BWF precipitant (BWF_Pre) were tested by a targeted panel of 1021 genes.

RESULTS

BWF_Sup cell-free DNA (cfDNA) was superior to tissue biopsy and BWF_Pre in determining mutational allele frequency, tumour mutational burden, and chromosomal instability. Moreover, BWF_Sup and BWF_Pre achieved comparable efficacy to tissue samples in differentiating malignant and benign patients, but only BWF_Sup persisted differentiated performance after excluding 55 malignancies pathologically diagnosed by bronchoscopic biopsy. Among 67 malignant patients, 82.1% and 71.6% of tumour-derived mutations (TDMs) were detected in BWF_Sup and BWF_Pre, respectively, and the detectability of TDMs in BWF_Sup was independent of the cytological examination of BWF. BWF_Sup outperformed BWF_Pre in providing more subclonal information and thus might yield advantage in tracking drug-resistant markers.

CONCLUSIONS

BWF_Sup cfDNA is a reliable medium for lung cancer diagnosis and genomic profiles and may provide important information for subsequent therapeutic regimens.

varAmpliCNV: analyzing variance of amplicons to detect CNVs in targeted NGS data

MOTIVATION

Computational identification of copy number variants (CNVs) in sequencing data is a challenging task. Existing CNV-detection methods account for various sources of variation and perform different normalization strategies. However, their applicability and predictions are restricted to specific enrichment protocols. Here, we introduce a novel tool named varAmpliCNV, specifically designed for CNV-detection in amplicon-based targeted resequencing data (Haloplex™ enrichment protocol) in the absence of matched controls. VarAmpliCNV utilizes principal component analysis (PCA) and/or metric dimensional scaling (MDS) to control variances of amplicon associated read counts enabling effective detection of CNV signals.

RESULTS

Performance of VarAmpliCNV was compared against three existing methods (ConVaDING, ONCOCNV and DECoN) on data of 167 samples run with an aortic aneurysm gene panel (n = 30), including 9 positive control samples. Additionally, we validated the performance on a large deafness gene panel (n = 145) run on 138 samples, containing 4 positive controls. VarAmpliCNV achieved higher sensitivity (100%) and specificity (99.78%) in comparison to competing methods. In addition, unsupervised clustering of CNV segments and visualization plots of amplicons spanning these regions are included as a downstream strategy to filter out false positives.

AVAILABILITY AND IMPLEMENTATION

The tool is freely available through galaxy toolshed and at: https://hub.docker.com/r/cmgantwerpen/varamplicnv. Supplementary Data File S1: https://tinyurl.com/2yzswyhh; Supplementary Data File S2: https://tinyurl.com/ycyf2fb4.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

SWI/SNF complex gene variations are associated with a higher tumor mutational burden and a better response to immune checkpoint inhibitor treatment: a pan-cancer analysis of next-generation sequencing data corresponding to 4591 cases

Background

Genes related to the SWItch/sucrose nonfermentable (SWI/SNF) chromatin remodeling complex are frequently mutated across cancers. SWI/SNF-mutant tumors are vulnerable to synthetic lethal inhibitors. However, the landscape of SWI/SNF mutations and their associations with tumor mutational burden (TMB), microsatellite instability (MSI) status, and response to immune checkpoint inhibitors (ICIs) have not been elucidated in large real-world Chinese patient cohorts.

Methods

The mutational rates and variation types of six SWI/SNF complex genes (ARID1A, ARID1B, ARID2, SMARCA4, SMARCB1, and PBRM1) were analyzed retrospectively by integrating next-generation sequencing data of 4591 cases covering 18 cancer types. Thereafter, characteristics of SWI/SNF mutations were depicted and the TMB and MSI status and therapeutic effects of ICIs in the SWI/SNF-mutant and SWI/SNF-non-mutant groups were compared.

Results

SWI/SNF mutations were observed in 21.8% of tumors. Endometrial (54.1%), gallbladder and biliary tract (43.4%), and gastric (33.9%) cancers exhibited remarkably higher SWI/SNF mutational rates than other malignancies. Further, ARID1A was the most frequently mutated SWI/SNF gene, and ARID1A D1850fs was identified as relatively crucial. The TMB value, TMB-high (TMB-H), and MSI-high (MSI-H) proportions corresponding to SWI/SNF-mutant cancers were significantly higher than those corresponding to SWI/SNF-non-mutant cancers (25.8 vs. 5.6 mutations/Mb, 44.3% vs. 10.3%, and 16.0% vs. 0.9%, respectively; all p  < 0.0001). Furthermore, these indices were even higher for tumors with co-mutations of SWI/SNF genes and MLL2/3. Regarding immunotherapeutic effects, patients with SWI/SNF variations showed significantly longer progression-free survival (PFS) rates than their SWI/SNF-non-mutant counterparts (hazard ratio [HR], 0.56 [95% confidence interval {CI} 0.44–0.72]; p < 0.0001), and PBRM1 mutations were associated with relatively better ICI treatment outcomes than the other SWI/SNF gene mutations (HR, 0.21 [95% CI 0.12–0.37]; p = 0.0007). Additionally, patients in the SWI/SNF-mutant + TMB-H (HR, 0.48 [95% CI 0.37–0.54]; p  < 0.0001) cohorts had longer PFS rates than those in the SWI/SNF-non-mutant + TMB-low cohort.

Conclusions

SWI/SNF complex genes are frequently mutated and are closely associated with TMB-H status, MSI-H status, and superior ICI treatment response in several cancers, such as colorectal cancer, gastric cancer, and non-small cell lung cancer. These findings emphasize the necessity and importance of molecular-level detection and interpretation of SWI/SNF complex mutations.

Genetic and methylation profiles distinguish benign, malignant and spitzoid melanocytic tumors

Accurate classification of melanocytic tumors is important for prognostic evaluation, treatment and follow-up protocols of patients. The majority of melanocytic proliferations can be classified solely based on clinical and pathological criteria, however in select cases a definitive diagnostic assessment remains challenging and additional diagnostic biomarkers would be advantageous. We analyzed melanomas, nevi, Spitz nevi and atypical spitzoid tumors using parallel sequencing (exons of 611 genes and 507 gene translocation analysis) and methylation arrays (850k Illumina EPIC). By combining detailed genetic and epigenetic analysis with reference-based and reference-free DNA methylome deconvolution we compared Spitz nevi to nevi and melanoma and assessed the potential for these methods in classifying challenging spitzoid tumors. Results were correlated with clinical and histologic features. Spitz nevi were found to cluster independently of nevi and melanoma and demonstrated a different mutation profile. Multiple copy number alterations and TERT promoter mutations were identified only in melanomas. Genome-wide methylation in Spitz nevi was comparable to benign nevi while the Leukocytes UnMethylation for Purity (LUMP) algorithm in Spitz nevi was comparable to melanoma. Histologically difficult to classify Spitz tumor cases were assessed which, based on methylation arrays, clustered between Spitz nevi and melanoma and in terms of genetic profile or copy number variations demonstrated worrisome features suggesting a malignant neoplasm. Comprehensive sequencing and methylation analysis verify Spitz nevi as an independent melanocytic entity distinct from both nevi and melanoma. Combined genetic and methylation assays can offer additional insights in diagnosing difficult to classify Spitzoid tumors.

NGS-Based Molecular Karyotyping of Multiple Myeloma: Results from the GEM12 Clinical Trial

Simple Summary

Multiple Myeloma (MM) is considered an incurable chronic disease, which prognosis depends on the presence of different genomic alterations. To accomplish a complete molecular diagnosis in a single essay, we have designed and validated a capture-based NGS approach to reliably identify pathogenic mutations (SNVs and indels), genomic alterations (CNVs and chromosomic translocations), and IGH rearrangements. We have observed a good correlation of the results obtained using our capture panel with data obtained by both FISH and WES techniques. In this study, the molecular classification performed using our approach was significantly associated with the stratification and outcome of MM patients. Additionally, this panel has been proven to detect specific IGH rearrangements that could be used as biomarkers in patient follow-ups through minimal residual disease (MRD) assays. In conclusion, we think that MM patients could benefit from the use of this capture-based NGS approach with a more accurate, single-essay molecular diagnosis.

Abstract

Next-generation sequencing (NGS) has greatly improved our ability to detect the genomic aberrations occurring in multiple myeloma (MM); however, its transfer to routine clinical labs and its validation in clinical trials remains to be established. We designed a capture-based NGS targeted panel to identify, in a single assay, known genetic alterations for the prognostic stratification of MM. The NGS panel was designed for the simultaneous study of single nucleotide and copy number variations, insertions and deletions, chromosomal translocations and V(D)J rearrangements. The panel was validated using a cohort of 149 MM patients enrolled in the GEM2012MENOS65 clinical trial. The results showed great global accuracy, with positive and negative predictive values close to 90% when compared with available data from fluorescence in situ hybridization and whole-exome sequencing. While the treatments used in the clinical trial showed high efficacy, patients defined as high-risk by the panel had shorter progression-free survival (p = 0.0015). As expected, the mutational status of TP53 was significant in predicting patient outcomes (p = 0.021). The NGS panel also efficiently detected clonal IGH rearrangements in 81% of patients. In conclusion, molecular karyotyping using a targeted NGS panel can identify relevant prognostic chromosomal abnormalities and translocations for the clinical management of MM patients.

HRD-MILN: Accurately estimate tumor homologous recombination deficiency status from targeted panel sequencing data

Homologous recombination deficiency (HRD) is a critical feature guiding drug and treatment selection, mainly for ovarian and breast cancers. As it cannot be directly observed, HRD status is estimated on a small set of genomic instability features from sequencing data. The existing methods often perform poorly when handling targeted panel sequencing data; however, the targeted panel is the most popular sequencing strategy in clinical practices. Thus, we proposed HRD-MILN to overcome the computational challenges from targeted panel sequencing. HRD-MILN incorporated a multi-instance learning framework to discover as many loss of heterozygosity (LOH) associated with HRD status to cluster as possible. Then the HRD score is obtained based on the association between the LOHs and the cluster in the sample to be estimated, and finally, the HRD status is estimated based on the score.

In comparison experiments on targeted panel sequencing data, the Precision of HRD-MILN could achieve 87%, significantly improved from 63% reported by the existing methods, where the highest margin of improvement reached 14%. It also presented advantages on whole exome sequencing data. Based on our best knowledge, HRD-MILN is the first practical tool for estimating HRD status from targeted panel sequencing data and could benefit clinical applications.

Liquid biopsy using ascitic fluid and pleural effusion supernatants for genomic profiling in gastrointestinal and lung cancers

Background

Precision medicine highlights the importance of incorporating molecular genetic testing into standard clinical care. Next-generation sequencing can detect cancer-specific gene mutations, and molecular-targeted drugs can be designed to be effective for one or more specific gene mutations. For patients with special site metastases, it is particularly important to use appropriate samples for genetic profiling. This study aimed to determine whether genomic profiling using ASC and PE is effective in detecting genetic mutations.

Methods

Tissues, plasma, ascites (ASC) supernatants, and pleural effusion (PE) samples from gastrointestinal cancer patients with peritoneal metastasis and lung cancer patients with pleural metastasis were collected for comprehensive genomic profiling. The samples were subjected to next-generation sequencing using a panel of 59 or 1021 cancer-relevant genes panel.

Results

A total of 156 tissues, 188 plasma samples, 45 ASC supernatants, and 1 PE samples from 304 gastrointestinal cancer patients and 446 PE supernatants, 122 tissues, 389 plasma samples, and 45 PE sediments from 407 lung cancer patients were analyzed. The MSAF was significantly higher in ASC and PE supernatant than that in plasma ctDNA (50.00% vs. 3.00%, p < 0.0001 and 28.5% vs. 1.30%, p < 0.0001, respectively). The ASC supernatant had a higher actionable mutation rate and more actionable alterations than the plasma ctDNA in 26 paired samples. The PE supernatant had a higher total actionable mutation rate than plasma (80.3% vs. 48.4%, p < 0.05). The PE supernatant had a higher frequency of uncommon variations than the plasma regardless of distant organ metastasis.

Conclusion

ASC and PE supernatants could be better alternative samples when tumor tissues are not available, especially in patients with only peritoneal or pleural metastases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09922-5.

Unraveling the Structural Variations of Early-Stage Mycosis Fungoides-CD3 Based Purification and Third Generation Sequencing as Novel Tools for the Genomic Landscape in CTCL

Mycosis fungoides (MF) is the most common cutaneous T-cell lymphoma (CTCL). At present, knowledge of genetic changes in early-stage MF is insufficient. Additionally, low tumor cell fraction renders calling of copy-number variations as the predominant mutations in MF challenging, thereby impeding further investigations. We show that enrichment of T cells from a biopsy of a stage I MF patient greatly increases tumor fraction. This improvement enables accurate calling of recurrent MF copy-number variants such as ARID1A and CDKN2A deletion and STAT5 amplification, undetected in the unprocessed biopsy. Furthermore, we demonstrate that application of long-read nanopore sequencing is especially useful for the structural variant rich CTCL. We detect the structural variants underlying recurrent MF copy-number variants and show phasing of multiple breakpoints into complex structural variant haplotypes. Additionally, we record multiple occurrences of templated insertion structural variants in this sample. Taken together, this study suggests a workflow to make the early stages of MF accessible for genetic analysis, and indicates long-read sequencing as a major tool for genetic analysis for MF.

PEcnv: accurate and efficient detection of copy number variations of various lengths

Copy number variation (CNV) is a class of key biomarkers in many complex traits and diseases. Detecting CNV from sequencing data is a substantial bioinformatics problem and a standard requirement in clinical practice. Although many proposed CNV detection approaches exist, the core statistical model at their foundation is weakened by two critical computational issues: (i) identifying the optimal setting on the sliding window and (ii) correcting for bias and noise. We designed a statistical process model to overcome these limitations by calculating regional read depths via an exponentially weighted moving average strategy. A one-run detection of CNVs of various lengths is then achieved by a dynamic sliding window, whose size is self-adopted according to the weighted averages. We also designed a novel bias/noise reduction model, accompanied by the moving average, which can handle complicated patterns and extend training data. This model, called PEcnv, accurately detects CNVs ranging from kb-scale to chromosome-arm level. The model performance was validated with simulation samples and real samples. Comparative analysis showed that PEcnv outperforms current popular approaches. Notably, PEcnv provided considerable advantages in detecting small CNVs (1 kb–1 Mb) in panel sequencing data. Thus, PEcnv fills the gap left by existing methods focusing on large CNVs. PEcnv may have broad applications in clinical testing where panel sequencing is the dominant strategy. Availability and implementation: Source code is freely available at https://github.com/Sherwin-xjtu/PEcnv

Genetic and Methylation Analysis of CTNNB1 in Benign and Malignant Melanocytic Lesions

Simple Summary

Recurrent CTNNB1 exon 3 mutations have been recognized in the distinct group of melanocytic tumors showing deep penetrating nevus-like morphology and in 1–2% of advanced melanoma. We performed a detailed genetic analysis of difficult-to-classify nevi and melanomas with CTNNB1 mutations and found that benign tumors (nevi) show characteristic morphological, genetic and epigenetic traits, which distinguish them from other nevi and melanoma. Malignant CTNNB1-mutant tumors (melanoma) demonstrated a different genetic profile, grouping clearly with other non-CTNNB1 melanomas in methylation assays. To further evaluate the role of CTNNB1 mutations in melanoma, we assessed a large cohort of clinically sequenced melanomas, identifying 38 tumors with CTNNB1 exon 3 mutations, including recurrent S45 (n = 13, 34%), G34 (n = 5, 13%), and S27 (n = 5, 13%) mutations. Locations and histological subtype of CTNNB1-mutated melanoma varied; none were reported as showing deep penetrating nevus-like morphology. The most frequent concurrent activating mutations were BRAF V600 (55%) and NRAS Q61 (34%).

Abstract

Melanocytic neoplasms have been genetically characterized in detail during the last decade. Recurrent CTNNB1 exon 3 mutations have been recognized in the distinct group of melanocytic tumors showing deep penetrating nevus-like morphology. In addition, they have been identified in 1–2% of advanced melanoma. Performing a detailed genetic analysis of difficult-to-classify nevi and melanomas with CTNNB1 mutations, we found that benign tumors (nevi) show characteristic morphological, genetic and epigenetic traits, which distinguish them from other nevi and melanoma. Malignant CTNNB1-mutant tumors (melanomas) demonstrated a different genetic profile, instead grouping clearly with other non-CTNNB1 melanomas in methylation assays. To further evaluate the role of CTNNB1 mutations in melanoma, we assessed a large cohort of clinically sequenced melanomas, identifying 38 tumors with CTNNB1 exon 3 mutations, including recurrent S45 (n = 13, 34%), G34 (n = 5, 13%), and S27 (n = 5, 13%) mutations. Locations and histological subtype of CTNNB1-mutated melanoma varied; none were reported as showing deep penetrating nevus-like morphology. The most frequent concurrent activating mutations were BRAF V600 (n = 21, 55%) and NRAS Q61 (n = 13, 34%). In our cohort, four of seven (58%) and one of nine (11%) patients treated with targeted therapy (BRAF and MEK Inhibitors) or immune-checkpoint therapy, respectively, showed disease control (partial response or stable disease). In summary, CTNNB1 mutations are associated with a unique melanocytic tumor type in benign tumors (nevi), which can be applied in a diagnostic setting. In advanced disease, no clear characteristics distinguishing CTNNB1-mutant from other melanomas were observed; however, studies of larger, optimally prospective, cohorts are warranted.

Clinical and pathological characteristics of familial melanoma with germline TERT promoter variants

Around 10% of melanoma occur in patients with a suspected familial predisposition. TERT promoter mutations are the most common somatic hotspot mutations in human cancers. However, only two families with germline mutations have been identified to date. We present detailed histological, clinical and molecular pathologic analyses of affected patients and details of newly identified individuals in one of these previously reported families. TERT (NM_198253.3) Chr.5:1,295,161T>C (c.-57 T>C) promoter variants were detected in all melanoma-affected (n=18) and one non-diseased family member. Median age at diagnosis was 30 years (n=18, range 16-46 years, 2 unknown). While most primary melanoma arose on the upper extremities (n=7, 21%) and were superficial spreading melanoma (SSM, n=8, 24%), many primary melanoma also originated from non UV-exposed mucosal (n=2, 6%) and acral (n=4, 12%) locations. One SSM sample harboured a Chr.5:1,295,228C>T TERT promoter region in addition to the germline Chr.5:1,295,161T>C variant, arguing additional pathway activation can support tumor pathogenesis. Patients treated with BRAF inhibitor and/or immune checkpoint inhibition (ICI) showed responses, although of limited duration. One mucosal melanoma harboured both a KIT copy number gain and an activating c.1727 p.Leu576Pro mutation. Following modest response to ICI, subsequent KIT inhibitor (imatinib) therapy demonstrated an ongoing complete pathological response (currently 7 months).

Algorithmic improvements for discovery of germline copy number variants in next-generation sequencing data

Background

Copy number variants (CNVs) play a significant role in human heredity and disease. However, sensitive and specific characterization of germline CNVs from NGS data has remained challenging, particularly for hybridization-capture data in which read counts are the primary source of copy number information.

Results

We describe two algorithmic adaptations that improve CNV detection accuracy in a Hidden Markov Model (HMM) context. First, we present a method for computing target- and copy number-specific emission distributions. Second, we demonstrate that the Pointwise Maximum a posteriori (PMAP) HMM decoding procedure yields improved sensitivity for small CNV calls compared to the more common Viterbi HMM decoder. We develop a prototype implementation, called Cobalt, and compare it to other CNV detection tools using sets of simulated and previously detected CNVs with sizes spanning a single exon to a full chromosome.

Conclusions

In both the simulation and previously detected CNV studies Cobalt shows similar sensitivity but significantly fewer false positive detections compared to other callers. Overall sensitivity is 80–90% for deletion CNVs spanning 1–4 targets and 90–100% for larger deletion events, while sensitivity is somewhat lower for small duplication CNVs.

Application of next-generation sequencing in the diagnosis of gastric cancer

Objectives: Gastric cancer (GC) is a disease with high mortality, poor prognosis and numerous risk factors. GC has an asymptomatic nature in early stages of the diseases, making timely diagnosis complicated using common conventional approaches, namely pathological examinations and imaging tests. Recently, molecular profiling of GC using next generation sequencing (NGS) has opened new doors to efficient prognostic, diagnostic, and therapeutic strategies. The current review aims to thoroughly discuss and compare the current NGS techniques and commercial platforms utilized for GC diagnosis and treatment, highlighting the most recent NGS-based GC studies. Furthermore, this review addresses the challenges of clinical implementation of NGS in GC.Materials and methods: This review was conducted according to the eligible studies identified via search of Web of Science, PubMed, Scopus, Embase and the Cochrane Library. In the present study, data on gastric cancer patients and NGS methods used to diagnose the disease were reviewed.Conclusion: Given the ever-rising advancements in NGS technologies, bioinformatics, healthcare guidelines and refined classifications, it is hoped that these technologies can actualize their advantages and optimize GC patients' experience.

Computational analysis of cancer genome sequencing data

Distilling biologically meaningful information from cancer genome sequencing data requires comprehensive identification of somatic alterations using rigorous computational methods. As the amount and complexity of sequencing data have increased, so has the number of tools for analysing them. Here, we describe the main steps involved in the bioinformatic analysis of cancer genomes, review key algorithmic developments and highlight popular tools and emerging technologies. These tools include those that identify point mutations, copy number alterations, structural variations and mutational signatures in cancer genomes. We also discuss issues in experimental design, the strengths and limitations of sequencing modalities and methodological challenges for the future.

Genomic and TCR profiling data reveal the distinct molecular traits in epithelial ovarian cancer histotypes

Epithelial ovarian cancer (EOC) is classified into five major histotypes: high-grade serous (HGSOC), low-grade serous (LGSOC), clear cell (CCOC), endometrioid (ENOC), and mucinous (MOC). However, the landscape of molecular and immunological alterations in these histotypes, especially LGSOC, CCOC, ENOC, and MOC, is largely uncharacterized. We collected 101 treatment-naive EOC patients. The resected tumor tissues and paired preoperative peripheral blood samples were collected and subjected to target sequencing of 1021 cancer-associated genes and T cell repertoire sequencing. Distinct characteristics of mutations were identified among the five histotypes. Furthermore, tumor mutation burden (TMB) was found to be higher in CCOC and ENOC, but lower in LGSOC and HGSOC. Alterations associated with DNA damage repair (DDR) pathways and homologous recombination deficiencies (HRD) were prevalent in five histotypes. CCOC demonstrated increased level of T cell clonality compared with HSGOC. Interestingly, the proportion of the 100 most common T cell clones was associated with TMB and tumor neoantigen burden in CCOC, highlighting more sensitive anti-tumor responses in this histotype, which was also evidenced by the enhanced convergent recombination of T cell clones. These findings shed light on the molecular traits of genomic alteration and T cell repertoire in the five major EOC histotypes and may help optimize clinical management of EOC with different histotypes.

Contribution of Clinical and Genetic Approaches for Diagnosing 209 Index Cases With 46,XY Differences of Sex Development

CONTEXT

Massively parallel sequencing (MPS) technologies have emerged as a first-tier approach for diagnosing several pediatric genetic syndromes. However, MPS has not been systematically integrated into the diagnostic workflow along with clinical/biochemical data for diagnosing 46,XY differences of sex development (DSD).

OBJECTIVE

To analyze the contribution of phenotypic classification either alone or in association with genetic evaluations, mainly MPS, for diagnosing a large cohort of 46,XY DSD patients.

DESIGN/PATIENTS

209 nonsyndromic 46,XY DSD index cases from a Brazilian DSD center were included. Patients were initially classified into 3 subgroups according to clinical and biochemical data: gonadal dysgenesis (GD), disorders of androgen secretion/action, and DSD of unknown etiology. Molecular genetic studies were performed by Sanger sequencing and/or MPS.

RESULTS

Clinical/biochemical classification into either GD or disorders of hormone secretion/action was obtained in 68.4% of the index cases. Among these, a molecular diagnosis was obtained in 36% and 96.5%, respectively. For the remainder 31.6% classified as DSD of clinically unknown etiology, a molecular diagnosis was achieved in 31.8%. Overall, the molecular diagnosis was achieved in 59.3% of the cohort. The combination of clinical/biochemical and molecular approaches diagnosed 78.9% of the patients. Clinical/biochemical classification matched with the genetic diagnosis in all except 1 case. DHX37 and NR5A1 variants were the most frequent genetic causes among patients with GD and DSD of clinical unknown etiology, respectively.

CONCLUSIONS

The combination of clinical/biochemical with genetic approaches significantly improved the diagnosis of 46,XY DSD. MPS potentially decreases the complexity of the diagnostic workup as a first-line approach for diagnosing 46,XY DSD.

CoverageMaster: comprehensive CNV detection and visualization from NGS short reads for genetic medicine applications

CoverageMaster (CoM) is a copy number variation (CNV) calling algorithm based on depth-of-coverage maps designed to detect CNVs of any size in exome [whole exome sequencing (WES)] and genome [whole genome sequencing (WGS)] data. The core of the algorithm is the compression of sequencing coverage data in a multiscale Wavelet space and the analysis through an iterative Hidden Markov Model. CoM processes WES and WGS data at nucleotide scale resolution and accurately detects and visualizes full size range CNVs, including single or partial exon deletions and duplications. The results obtained with this approach support the possibility for coverage-based CNV callers to replace probe-based methods such as array comparative genomic hybridization and multiplex ligation-dependent probe amplification in the near future.

DECoN: A Detection and Visualization Tool for Exonic Copy Number Variants

Detection of copy number variants from targeted sequencing, including whole-exome sequencing, can be particularly difficult since the break points of the CNV are not always captured. Here we describe DECoN, a software tool which uses changes in read depth to identify CNVs that affect whole exons. It is optimized for clinical use and allows for interactive visualization of CNVs identified.

Detection of Small CYP11B1 Deletions and One Founder Chimeric CYP11B2/CYP11B1 Gene in 11β-Hydroxylase Deficiency

OBJECTIVE

11β-Hydroxylase deficiency (11β-OHD) caused by mutations in the CYP11B1 gene is the second most common form of congenital adrenal hyperplasia. Both point mutations and genomic rearrangements of CYP11B1 are important causes of 11β-OHD. However, the high degree of sequence identity between CYP11B1 and its homologous gene CYP11B2, presents unique challenges for molecular diagnosis of suspected 11β-OHD. The aim of this study was to detect the point mutation, indel, small deletion of CYP11B1 and chimeric CYP11B2/CYP11B1 gene in a one-tube test, improving the genetic diagnosis of 11β-OHD.

METHODS

Optimized custom-designed target sequencing strategy was performed in three patients with suspected 11β-OHD, in which both the coverage depth of paired-end reads and the breakpoint information of split reads from sequencing data were analysed in order to detect genomic rearrangements covering CYP11B1. Long-range PCR was peformed to validate the speculated CYP11B1 rearrangements with the breakpoint-specifc primers.

RESULTS

Using the optimized target sequencing approach, we detected two intragenic/intergenic deletions of CYP11B1 and one chimeric CYP11B2/CYP11B1 gene from three suspected patients with 11β-OHD besides three pathogenic heterozygous point mutation/indels. Furthermore, we mapped the precise breakpoint of this chimeric CYP11B2/CYP11B1 gene located on chr8:143994517 (hg19) and confirmed it as a founder rearrangement event in the Chinese population.

CONCLUSIONS

Our optimized target sequencing approach improved the genetic diagnosis of 11β-OHD.