Copy Number Variants Burden Research Articles

Genome-wide copy number variation association study in anorexia nervosa.

This study represents the first large-scale investigation of rare (<1% population frequency) copy number variants (CNVs) in anorexia nervosa (AN). Large, rare CNVs are reported to be causally associated with anthropometric traits, neurodevelopmental disorders, and schizophrenia, yet their role in the genetic basis of AN is unclear. Using genome-wide association study (GWAS) array data from the Anorexia Nervosa Genetics Initiative (ANGI), which included 7414 AN case and 5044 controls, we investigated the association of 67 well-established syndromic CNVs and 178 pleiotropic disease-risk dosage-sensitive CNVs with AN. To identify novel CNV regions (CNVRs) that increase the risk of AN, we conducted genome-wide association studies with a focus on rare CNV-breakpoints (CNV-GWAS). We found no net enrichment of rare CNVs, either deletions or duplications, in AN, and none of the well-established syndromic or pleiotropic CNVs had a significant association with AN status. However, the CNV-GWAS found 21 nominally associated CNVRs that contribute to AN risk, covering protein-coding genes implicated in synaptic function, metabolic/mitochondrial factors, and lipid characteristics, like the CD36 (7q21.11) gene, which transports long-chain fatty acids into cells. CNVRs intersecting genes previously related to neurodevelopmental traits include deletions of NRXN1 intron 5 (2p16.3), IMMP2L (7q31.1), and PTPRD (9p23). Overall, given that our study is well powered to detect the CNV burden level reported for schizophrenia, we can conclude that rare CNVs have a limited role in the etiology of AN, as reported for bipolar disorder. Our nominal associations for the 21 discovered CNVRs are consistent with AN being a metabo-psychiatric trait, as demonstrated by the common genetic architecture of AN, and we provide association results to allow for replication in future research.

A burden of rare copy number variants in obsessive-compulsive disorder.

Current genetic research on obsessive-compulsive disorder (OCD) supports contributions to risk specifically from common single nucleotide variants (SNVs), along with rare coding SNVs and small insertion-deletions (indels). The contribution to OCD risk from rare copy number variants (CNVs), however, has not been formally assessed at a similar scale. Here we describe an analysis of rare CNVs called from genotype array data in 2248 deeply phenotyped OCD cases and 3608 unaffected controls from Sweden and Norway. Cases carry an elevated burden of CNVs ≥30 kb in size (OR = 1.12, P = 1.77 × 10-3). The excess rate of these CNVs in cases versus controls was around 0.07 (95% CI 0.02-0.11, P = 2.58 × 10-3). This signal was largely driven by CNVs overlapping protein-coding regions (OR = 1.19, P = 3.08 × 10-4), particularly deletions impacting loss-of-function intolerant genes (pLI >0.995, OR = 4.12, P = 2.54 × 10-5). We did not identify any specific locus where CNV burden was associated with OCD case status at genome-wide significance, but we noted non-random recurrence of CNV deletions in cases (permutation P = 2.60 × 10-3). In cases where sufficient clinical data were available (n = 1612) we found that carriers of neurodevelopmental duplications were more likely to have comorbid autism (P < 0.001), and that carriers of deletions overlapping neurodevelopmental genes had lower treatment response (P = 0.02). The results demonstrate a contribution of rare CNVs to OCD risk, and suggest that studies of rare coding variation in OCD would have increased power to identify risk genes if this class of variation were incorporated into formal tests.

Genome-wide association study of copy number variations in Parkinson's disease.

Our study investigates the impact of copy number variations (CNVs) on Parkinson's disease (PD) pathogenesis using genome-wide data, aiming to uncover novel genetic mechanisms and improve the understanding of the role of CNVs in sporadic PD. We applied a sliding window approach to perform CNV-GWAS and conducted genome-wide burden analyses on CNV data from 11,035 PD patients (including 2,731 early-onset PD (EOPD)) and 8,901 controls from the COURAGE-PD consortium. We identified 14 genome-wide significant CNV loci associated with PD, including one deletion and 13 duplications. Among these, duplications in 7q22.1, 11q12.3 and 7q33 displayed the highest effect. Two significant duplications overlapped with PD-related genes SNCA and VPS13C, but none overlapped with recent significant SNP-based GWAS findings. Five duplications included genes associated with neurological disease, and four overlapping genes were dosage-sensitive and intolerant to loss-of-function variants. Enriched pathways included neurodegeneration, steroid hormone biosynthesis, and lipid metabolism. In early-onset cases, four loci were significantly associated with EOPD, including three known duplications and one novel deletion in PRKN. CNV burden analysis showed a higher prevalence of CNVs in PD-related genes in patients compared to controls (OR=1.56 [1.18-2.09], p=0.0013), with PRKN showing the highest burden (OR=1.47 [1.10-1.98], p=0.026). Patients with CNVs in PRKN had an earlier disease onset. Burden analysis with controls and EOPD patients showed similar results. This is the largest CNV-based GWAS in PD identifying novel CNV regions and confirming the significant CNV burden in EOPD, primarily driven by the PRKN gene, warranting further investigation.

Contribution of copy number variants on antipsychotic treatment response in Han Chinese patients with schizophrenia

Response to antipsychotic drugs (APD) varies greatly among individuals and is affected by genetic factors. This study aims to demonstrate genome-wide associations between copy number variants (CNVs) and response to APD in patients with schizophrenia. A total of 3030 patients of Han Chinese ethnicity randomly received APD (aripiprazole, olanzapine, quetiapine, risperidone, ziprasidone, haloperidol and perphenazine) treatment for six weeks. This study is a secondary data analysis. Percentage change on the Positive and Negative Syndrome Scale (PANSS) reduction was used to assess APD efficacy, and more than 50% change was considered as APD response. Associations between CNV burden, gene set, CNV loci and CNV break-point and APD efficacy were analysed. Higher CNV losses burden decreased the odds of 6-week APD response (OR=0.66 [0.44, 0.98]). CNV losses in synaptic pathway involved in neurotransmitters were associated with 2-week PANSS reduction rate. CNV involved in sialylation (1p31.1 losses) and cellular metabolism (19q13.32 gains) associated with 6-week PANSS reduction rate at genome-wide significant level. Additional 36 CNVs associated with PANSS factors improvement. The OR of protective CNVs for 6-week APD response was 3.10 (95% CI: 1.33-7.19) and risk CNVs was 8.47 (95% CI: 1.92-37.43). CNV interacted with genetic risk score on APD efficacy (Beta=-1.53, SE=0.66, P=0.021). The area under curve to differ 6-week APD response attained 80.45% (95% CI: 78.07%-82.82%). Copy number variants contributed to poor APD efficacy and synaptic pathway involved in neurotransmitter was highlighted. National Natural Science Foundation of China, National Key R&D Program of China, China Postdoctoral Science Foundation.

Abstract 5211: Somatic copy number variation as prognostic marker for recurrence in never-smokers with early-stage lung adenocarcinoma

Abstract Objective: Lung cancer in never-smokers has a distinct genomic profile with low tumor mutation burden, more tumor promotor mutations and fewer somatic copy number variations (CNV) than smokers. CNVs have shown to be associated with poor outcome in advanced stage disease. However, the relationship between CNV and clinicopathologic features on prognosis in early-stage lung cancer is undefined. To address this knowledge gap, we assessed the impact of CNV on recurrence in never-smokers who underwent curative resection for pathological stage I-II lung adenocarcinoma using an integrated genomic and clinicopathological analysis. Methods: We analyzed a cohort of 210 never-smokers with stage I-II lung adenocarcinoma treated from 2014 to 2022 who met the inclusion criteria (no neoadjuvant therapy, pathological stage I-II, complete resection, next-generation sequencing available). The patient cohort was stratified using unsupervised hierarchical clustering of arm-level CNVs. We assessed cumulative incidence of recurrence (CIR) between CNV groups using competing risk regression analysis, adjusting for SUVmax, pathologic stage, lymphovascular invasion, IASLC grade, and tumor mutation burden. Results: Of the 210 patients, 159 patients were female (76%), and patients were primarily from Caucasian (n=140, 67%) or Asian decent (n=50, 24%). Based on arm-level CNV clustering, the cohort was stratified into three groups: CNV low (n=86, 41%), CNV moderate (n=75, 36%) and CNV high (n=49, 23%). The CNV low tumors were characterized by a low CNV burden. Conversely, the CNV moderate tumors primarily exhibited gains in chromosomes 1q, 5p and 7p, along with loss of heterozygosity in chromosome 9p, 9q and 13. CNV high tumors were characterized by whole-genome doubling. Whole-genome doubling was significantly more common in CNV high tumors (n=33, 66%) compared to CNV low (n=2, 2.3%) and CNV moderate tumors (n=2, 2.7%) (p&amp;lt;.001). EGFR (n=145, 69%) and TP53 (n=61, 29%) were the most frequently altered genes, with EGFR alterations being significantly more prevalent in CNV high tumors (n=41, 84%) compared to CNV low tumors (n=46, 53%) (p&amp;lt;.001). In total, 42 patients (20%) developed recurrences. The CIR was lower in the CNV low group (5y-CIR, 14%, 95%CI: 6.5%-25%) compared to the CNV moderate (5y-CIR, 31%, 95%CI: 18%-44%) and CNV high group (5y-CIR, 48%, 95%CI: 27%-66%) (p=.004), with adjusted hazard ratios of 2.3 (95%CI: 0.93-5.71; p=0.074) and 2.65 (95%CI: 1.07-6.60; p=.037), respectively. The CIR for locoregional (n=11, 5%) and distant recurrence (n=31, 15%) remained lower in the CNV low group compared to the combined CNV moderate &amp; high groups (Gray test p=.020 and p=.049, respectively). Conclusion: In never-smokers who underwent curative treatment of pathological stage I-II lung adenocarcinoma, copy number status seems a predictor of recurrence, with CNV low tumors associated with the best prognosis. Citation Format: Stijn Vanstraelen, Allison Reiner, Brooke H. Mastrogiacomo, Kay See Tan, Joseph Dycoco, Bernard J. Park, Prasad S. Adusumilli, Matthew J. Bott, David R. Jones, Gaetano Rocco. Somatic copy number variation as prognostic marker for recurrence in never-smokers with early-stage lung adenocarcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 5211.

High overall copy number variation burden by genome-wide methylation profiling holds negative prognostic value in surgically treated pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive type of cancer with an overall 5-year survival of around 10 %. New prognostic tools to stratify patients are needed. Our main aim was to evaluate the prognostic value of overall copy number variation (CNV) burden in surgically treated PDAC. DNA extracted from 108 surgical PDAC specimens was examined to collect data on the genome-wide DNA methylation status of >850,000 CpG sites in promoter, gene body, and enhancer regions (Illumina Infinium Methylation EPIC BeadChip Kit). CNV profiles were obtained and all PDACs were stratified into one of three groups: Low, moderate, or high overall CNV burden. Tumors histologically showing a dominant conventional and/or tubulopapillary pattern in 60 %–100 % and 0–59 % were categorized as Group A and Group B as per Kalimuthu. We also performed targeted next-generation sequencing (NGS) and immunohistochemistry. High overall CNV burden held independent negative prognostic value with poor survival (HR 4.01 (95%CI 1.96–8.19), p = 0.00014) and was more frequent in Group B (p = 0.0003). Most frequent chromosomal arm-level aberrations were gains of 8q (29 %) and 1q (19 %) and losses of 17p (55 %), 18q (43 %), 6q (37 %), 9p (36 %), 6p (26 %), 19p (26 %), and 8p (25 %). Most frequent mutations found were in KRAS (95 %), TP53 (62 %), CDKN2A (24 %), SMAD4 (23 %), ATM (9 %), ARID1A (7 %), RNF43 (7 %), GNAS (6 %), and KDM6A (6 %). Group A PDACs showed more frequently KRAS variants other than Gly12Val and Gly12Asp (p = 0.012). Our data indicate that overall CNV burden using genome-wide methylation profiling may be a useful prognostic tool in surgically treated PDAC. Importantly, our approach, using data from genome-wide methylation profiling for analysis of overall CNV burden, can be performed on formalin-fixed and paraffin embedded PDAC tissues. Future studies should examine the prognostic value of overall CNV burden in unresectable PDAC.

OR20-04 Genetic Profiling Using A Gene Panel Identifies Pathogenic Variants And Copy Number Variations In Pituitary Adenomas

Abstract Disclosure: T. Silva: None. D. Mohan: None. J. Mebarak: None. D. Meredith: None. W. Bi: Speaker; Self; Stryker. R. Carroll: None. U.B. Kaiser: None. A. Marcondes Lerario: None. A. Abreu: None. Introduction: Pituitary adenomas (PA) are the second most common central nervous system neoplasm and arise from different cell lineages. PA are benign tumors, although a subset of them may present with aggressive behavior, resulting in significant morbidity. Several genetic alterations have been implicated in PA tumorigenesis, such as USP8 mutations in Cushing disease (CD) and GNAS mutations in somatotroph adenomas (SA). However, the genetic profiles of PA are not completely understood. The aim of this study was to identify genetic alterations associated with PA using a gene panel. Methods: Targeted next generation exome sequencing (OncoPanel) was performed to identify single nucleotide variants (SNV) and copy number variations (CNV) associated with PA. OncoPanel version 1 included 275 genes, version 2 had 300 genes, and the third, current version 447 genes. USP8 was included only in the most recent version. Patients undergoing surgery at our institution between 2013-2020, who agreed to participate in the OncoPanel initiative and who had samples with good quality DNA were enrolled. Classification of PA into subtypes was performed by combining clinical and histopathological data. A strict algorithm was used to assess SNV pathogenicity. Chi-square, Mann Whitney and Kruskal-Wallis tests were used for statistical analyses. Results: Our cohort included 171 patients with following PA subtypes: TPIT lineage - 18 CD, 14 silent ACTH; PIT-1 lineage - 29 SA, 21 lactotroph and 4 silent poorly differentiated PIT-1 adenomas; SF-1 lineage - 77 gonadotroph and 8 null cell adenomas (negative immunostaining for hormones and TF). GNAS mutations were found in 28% of SA and USP8 mutations in 50% of CD screened for USP8. Pathogenic variants in MECOM were associated with silent ACTH, in FGFR4 and NF2 with CD, and in NOTCH2 with lactotroph adenomas, whereas pathogenic variants in TP53, EGFR, MEN1, BRCA1, FH, ATR, BLM, and XPC were identified in multiple PA subtypes. The number of CNV varied significantly between different subtypes of adenomas and were higher in lactotroph (median, 46), poorly differentiated PIT-1 (41), SA (26) and silent ACTH (17) compared to CD (3), gonadotroph (1) and null cell (0.5) (p&amp;lt; 0.001). MIB-1 proliferation index &amp;gt;3% was associated with higher CNV burden (p=0.022). Interestingly, copy number loss of chromosome 16 and 1q combined were exclusively identified in SA without GNAS mutations. Conclusion: OncoPanel analyses identified pathogenic variants in PA, some well recognized variants such as USP8 and GNAS, and some rare events such as NF2 in CD. In addition, our findings highlight subtype-specific patterns of genomic instability. The association between high CNV burden and the MIB-1 proliferation index may help to predict adenoma behavior. Genomic profile using a gene panel has the potential to better classify PA and contribute to clinical management. Presentation: Saturday, June 17, 2023

Sintilimab plus bevacizumab, oxaliplatin and capecitabine as first-line therapy in RAS-mutant, microsatellite stable, unresectable metastatic colorectal cancer: an open-label, single-arm, phase II trial

Microsatellite stable (MSS) and RAS-mutant metastatic colorectal cancer (mCRC) patients are characterized by an immunosuppressive microenvironment and a low response rate to immunotherapy. Chemotherapy and anti-angiogenesis therapy have been reported to potentially promote immunotherapy response. This study aims to assess the preliminary anti-tumor activity and safety of sintilimab plus bevacizumab, oxaliplatin and capecitabine as a treatment option for patients with RAS-mutant MSS mCRC. This study was an open-label, single-arm, phase II trial in China. Patients with unresectable, RAS-mutant and MSS metastatic colorectal adenocarcinoma received treatment by intravenous sintilimab (200mg, day 1) plus bevacizumab (7.5mg/kg, day 1), oxaliplatin (135mg/m2, day 1) and oral capecitabine (1g/m2, day 1-14) in each 21-day cycle. The primary endpoints included objective response rate (ORR) and adverse events. Biomarker analysis was performed to identify potential predictors of good response to treatment. This study is registered with ClinicalTrials.gov, number NCT04194359. Between April 2021 and December 2021, 25 patients were enrolled. Two (8%) patients showed complete response (CR), 19 (76%) had partial response (PR) and 4 (16%) presented with stable disease. ORR reached 84% (95% CI, 63.9-95.5) and the disease control rate was 100% (95% CI, 86.3-100). The median progression-free survival (PFS) was 18.2 months for the full analysis set. The most common treatment-related adverse events (TRAEs) in all grades were anemia (21/25, 84%), neutropenia (20/25, 80%), and hand-foot syndrome (14/25, 56%). The most frequent grade 3 or 4 TRAEs were neutropenia (3/25, 12%) and increased alanine transaminase (2/25, 8%). No grade 5 adverse events occurred. In the exploration of biomarkers, 5 patients could be characterized as TTN/OBSCN "double-hit" after treatment, and the copy number variants burden was significantly decreased in tumor tissues after treatment compared with the baseline. Nanostring panel RNA sequencing analysis indicated a better tumor immune microenvironment cell infiltration in CR/PR patients compared with non-CR/PR patients as well as the PFS-long (≥12.5 months) group compared with the PFS-short group. Combination treatment with sintilimab plus bevacizumab, oxaliplatin and capecitabine as first-line treatment demonstrated a promising antitumor activity and a manageable safety profile in RAS-mutant, MSS and unresectable mCRC. Exploratory biomarker assessment analysis showed that some RAS-mutant and MSS patients changed into "immune-hot" subtype after the treatment. This study was supported by the Key R&D Program of Zhejiang Province (2021C03125 to Ying Yuan), the National Natural Science Foundation of China (81872481 to Ying Yuan, 82072624 to Kefeng Ding), the Fundamental Research Funds for the Central Universities (No. 226-2022-00009 to Kefeng Ding), and the Zhejiang Provincial Natural Science Foundation of China (No. LY22H160024 to Hanguang Hu).

The implication of chromosomal abnormalities in the surgical outcomes of Chinese pediatric patients with congenital heart disease.

Copy number variations (CNVs) have been shown to be overrepresented in children with congenital heart disease (CHD). Genetic evaluation of CHD is currently underperformed in China. We sought to determine the occurrence of CNVs in CNV regions with disease-causing potential among a large cohort of Chinese pediatric CHD patients and investigate whether these CNVs could be the important critical modifiers of surgical intervention. CNVs screenings were performed in 1,762 Chinese children who underwent at least one cardiac surgery. CNV status at over 200 CNV locus with disease-causing potential was analyzed with a high-throughput ligation-dependent probe amplification (HLPA) assay. We found 378 out of 1,762 samples (21.45%) to have at least one CNV and 2.38% of them were carrying multiple CNVs. The detection rates of ppCNVs (pathogenic and likely pathogenic CNVs) were 9.19% (162/1,762), significantly higher than that of the healthy Han Chinese individuals from The Database of Genomic Variants archive (9.19% vs. 3.63%; P = 0.0012). CHD cases with ppCNVs had a significantly higher proportion of complex surgeries compared to CHD patients with no ppCNVs (62.35% vs. 37.63%, P < 0.001). Duration of cardiopulmonary bypass and aortic cross clamp procedures were significantly longer in CHD cases with ppCNVs (all P < 0.05), while no group differences were identified for complications of surgery and one-month mortality after surgery. The detection rate of ppCNVs in the atrioventricular septal defect (AVSD) subgroup was significantly higher than that in other subgroups (23.10% vs. 9.70%, P = 0.002). CNV burden is an important contributor to Chinese children with CHD. Our study demonstrated the robustness and diagnostic efficiency of HLPA method in the genetic screening of CNVs in CHD patients.

Molecular genetic characteristics of thymic epithelial tumors with distinct histological subtypes

Due to the low incidence and histological heterogeneity, the molecular features and underlying carcinogenic mechanisms of thymic epithelial tumors (TETs) are yet to be fully elucidated, especially for different subtypes of TETs. Tumor tissue samples of 43 TETs with distinct histological subtypes were collected. We analyzed the molecular characteristics in different subtypes based on whole exome sequencing data. The mutational profiles of the different subtypes of TETs varied. Compared with thymomas, thymic carcinomas (TCs) had a higher mutation frequency of MYO16 (33% vs. 3%, p= 0.024) and a lower frequency of ZNF729 mutations (0% vs. 35%, p=0.044). No significant difference was observed in the median tumor mutation burden across different subtypes. The value of copy number variation burden, weighted genome instability index, and the number of amplified segments were all higher in TCs than thymomas, and they also tended to be higher in B3 thymoma than in non-B3 thymomas, while they had no significant differences between B3 thymoma and TCs. Clustering analyses revealed that Wnt, MAPK, Hedgehog, AMPK, and cell junction assembly signaling pathways were exclusively enriched in non-B3 thymomas, lysine degradation pathway in B3 thymoma, and extracellular matrix-receptor (ECM-receptor) interaction, positive regulation of cell cycle process, and activation of innate immune response pathways in TCs. This study revealed distinct molecular landscapes of different subtypes of TETs, suggesting diverse pathogenesis of non-B3 thymomas, B3 thymomas, and TCs. Our findings warrant further validation in future large-scale studies and may provide a theoretical basis for potential personalized therapeutic strategies.

Diagnosis of prostate cancer using cell-free DNA methylation profiles from expressed prostatic secretions.

389 Background: Urinary and seminal cell-free DNA (cfDNA) have been recognized as promising biomarkers of prostate cancer. DNA methylation signal in tumor is increasingly used as tumor diagnosis and longitudinal monitoring indicator. However, the clinical utility of cfDNA from expressed prostatic secretions (EPS) remains unknown. Methods: The prospective study includes 50 prostate cancer (PC) patients and equivalent benign prostatic hyperplasia (BPH) patients, where EPS samples were collected after the prostatic massage of each patient. Cell-free DNA from EPS was extracted and treated by PredicineEPIC, a liquid biopsy comprehensive DNA methylation assay. The study developed integrated bioinformatic algorithms to profile the genome-wide epigenomic characteristics and investigate tumor-specific methylation patterns. Results: The initial exploratory cohort included 5 PC and 11 BPH cases. PredicineEPIC whole-genome DNA methylation of EPS identified over 300 differentially methylated regions (DMRs). Prostate cancers were distinguished from the BPH group by an unsupervised clustering method, suggesting that the identified DMRs embody the specific profiles of malignant tumors. Additionally, this study detected genome-wide copy number variation burden (CNB) in parallel. A risk model was built for cancer risk assessment based on methylation and CNB characteristics. Conclusions: This study demonstrated the feasibility of methylation profiling of cfDNA in EPS in prostate cancer. This non-invasive liquid biopsy approach could sensitively navigate prostate cancer from benign prostatic hyperplasia, suggesting future direction of methylation-based liquid biopsy in detecting prostate cancer.

Rare copy number variation analysis identifies disease-related variants in atrioventricular septal defect patients.

Atrioventricular septal defect (AVSD) is a deleterious subtype of congenital heart diseases (CHD) characterized by atrioventricular canal defect. The pathogenic genetic changes of AVSD remain elusive, particularly for copy number variation (CNV), a large segment variation of the genome, which is one of the major forms of genetic variants resulting in congenital heart diseases. In the present study, we recruited 150 AVSD cases and 100 healthy subjects as controls for whole exome sequencing (WES). We identified total 4255 rare CNVs using exon Hidden Markov model (XHMM) and screened rare CNVs by eliminating common CNVs based on controls and Database of Genomic Variants (DGV). Each patient contained at least 9 CNVs, and the CNV burden was prominently presented in chromosomes 19,22,21&16. Small CNVs (<500kb) were frequently observed. By leveraging gene-based burden test, we further identified 20 candidate AVSD-risk genes. Among them, DYRK1A, OBSCN and TTN were presented in the core disease network of CHD and highly and dynamically expressed in the heart during the development, which indicated they possessed the high potency to be AVSD-susceptible genes. These findings not only provided a roadmap for finally unveiling the genetic cause of AVSD, but also provided more resources and proofs for clinical genetics.

Esophageal Squamous Cancer from 4NQO-Induced Mice Model: CNV Alterations.

Squamous esophageal carcinoma is a common pathological type of esophageal carcinoma around the world. The prognosis of esophageal carcinoma is usually poor and diagnosed at late stages. Recently, research suggested that genomic instability occurred in esophageal cells during the development of esophageal squamous cell carcinoma (ESCC). Identifying prognostic and specific genomic characteristics, especially at the early hyperplasia stage, is critical. Mice were given 4-nitroquinoline 1-oxide (4NQO) with drinking water to induce esophageal cancer. The immortalized human esophageal epithelial cell line (NE2) was also treated with 4NQO. We performed histologic analyses, immunofluorescence, and immunohistochemical staining to detect DNA damage at different time points. Whole-exome sequencing was accomplished on the esophagus tissues at different pathological stages to detect single-nucleotide variants and copy number variation (CNV) in the genome. Our findings indicate that all mice were tumor-forming, and a series of changes from simple hyperplasia (ESSH) to intraepithelial neoplasia (IEN) to esophageal squamous cell carcinoma (ESCC) was seen at different times. The expression of γ-H2AX increased from ESSH to ESCC. In addition, mutations of the Muc4 gene were detected throughout the pathological stages. Furthermore, CNV burden appeared in the esophageal tissues from the beginning of ESSH and accumulated more in cancer with the deepening of the lesions. This study demonstrates that mutations caused by the early appearance of DNA damage may appear in the early stage of malignant tissue before the emergence of atypia. The detection of CNV and mutations of the Muc4 gene may be used as an ultra-early screening indicator for esophageal cancer.

Abstract 15501: Rare Genomic Copy Number Variants Implicate New Candidate Genes for Bicuspid Aortic Valve

Introduction: Bicuspid Aortic Valve (BAV), the most common adult congenital heart defect, is a major cause of aortic insufficiency or stenosis requiring valve replacement and thoracic aortic aneurysms predisposing to acute aortic dissections. The spectrum of BAV ranges from early onset valve and aortic complications to sporadic late onset disease. We determined the frequency and gene content of rare genomic copy number variants (CNVs) in isolated BAV cases. Methods: We performed genome-wide SNP microarray analysis of familial BAV cases with early onset complications in the UTHealth BAV Research Registry (EBAV, n=394) and elderly probands from the International BAV Consortium (BAVCON, n=4216). CNVs were detected in Illumina genotypes using the PennCNV, QuantiSNP and cnvPartition algorithms. BAV cases were compared to 16,576 controls without known cardiovascular disease from the Database of Genotypes and Phenotypes. Only CNV calls with &gt;6 consecutive variants, &gt;20 kb in length, identified by at least two algorithms and intersected with known genes were included. For comparison, we assessed a cohort with left ventricular outflow tract obstructive lesions including BAV (LVOTO, n=1561) using identical methods. CNV burden and associations were tested using PLINK. Results: We identified 84 large recurrent CNVs in BAV cases that are absent or rare (&lt;0.1%) in controls. 34 rare CNVs overlap between BAV and LVOTO cases and involve candidate genes that interact with each other during heart development. The largest and most prevalent of the recurrent CNVs are 8p23 duplications involving GATA4 (OR 340, CI: 42-2700). CNVs of genes that cause BAV when mutated were significantly enriched in cases (n=135, P =1x10 -27 ). The overall burden of rare genic CNVs, specifically large deletions, was also increased in BAV cases (empiric P &lt; 1x10 -5 ) and was higher in younger EBAV cases than in older BAVCON cases. Conclusion: We identified likely pathogenic CNVs in more than 10% of BAV cases, implicating alterations of candidate genes at these loci in the pathogenesis of BAV.

Identification and functional interpretation of miRNAs affected by rare CNVs in CAKUT

Rare copy number variants (CNVs) are among the most common genomic disorders underlying CAKUT. miRNAs located in rare CNVs represent well-founded functional variants for human CAKUT research. The study aimed to identify and functionally interpret miRNAs most frequently affected by rare CNVs in CAKUT and to estimate the overall burden of rare CNVs on miRNA genes in CAKUT. The additional aim of this study was to experimentally confirm the effect of a rare CNV in CAKUT on candidate miRNA’s expression and the subsequent change in mRNA levels of selected target genes. A database of CAKUT-associated rare CNV regions, created by literature mining, was used for mapping of the miRNA precursors. miRNAs and miRNA families, most frequently affected by rare CAKUT-associated CNVs, have been subjected to bioinformatic analysis. CNV burden analysis was performed to identify chromosomes with over/underrepresentation of miRNA genes in rare CNVs associated with CAKUT. A functional study was performed on HEK293 MIR484+/- KO and HEK293 WT cell lines, followed by the analysis of relative miRNA and mRNA target gene levels. 80% of CAKUT patients with underlying rare CNV had at least one miRNA gene overlapping the identified CNV. Network analysis of the most frequently affected miRNAs has revealed the dominant regulation of the two miRNAs, hsa-miR-484 and hsa-miR-185-5p. Additionally, miR-548 family members have shown substantial enrichment in rare CNVs in CAKUT. An over/underrepresentation of miRNA genes in rare CNVs associated with CAKUT was observed in multiple chromosomes, such as chr16, chr20, and chr21. A significant 0.37 fold downregulation of hsa-miR-484, followed by a notable upregulation of MDM2 and APAF1 and downregulation of NOTCH3 was detected in HEK293 MIR484+/- KO compared to HEK293 WT cell lines, supporting the study hypothesis. miRNA genes are frequently affected by rare CNVs in CAKUT patients. Understanding the potential of CNV-affected miRNAs to participate in CAKUT as genetic drivers represent a crucial implication for the development of novel therapeutic approaches.

Rare copy number variation in posttraumatic stress disorder

Posttraumatic stress disorder (PTSD) is a heritable (h2 = 24–71%) psychiatric illness. Copy number variation (CNV) is a form of rare genetic variation that has been implicated in the etiology of psychiatric disorders, but no large-scale investigation of CNV in PTSD has been performed. We present an association study of CNV burden and PTSD symptoms in a sample of 114,383 participants (13,036 cases and 101,347 controls) of European ancestry. CNVs were called using two calling algorithms and intersected to a consensus set. Quality control was performed to remove strong outlier samples. CNVs were examined for association with PTSD within each cohort using linear or logistic regression analysis adjusted for population structure and CNV quality metrics, then inverse variance weighted meta-analyzed across cohorts. We examined the genome-wide total span of CNVs, enrichment of CNVs within specified gene-sets, and CNVs overlapping individual genes and implicated neurodevelopmental regions. The total distance covered by deletions crossing over known neurodevelopmental CNV regions was significant (beta = 0.029, SE = 0.005, P = 6.3 × 10−8). The genome-wide neurodevelopmental CNV burden identified explains 0.034% of the variation in PTSD symptoms. The 15q11.2 BP1-BP2 microdeletion region was significantly associated with PTSD (beta = 0.0206, SE = 0.0056, P = 0.0002). No individual significant genes interrupted by CNV were identified. 22 gene pathways related to the function of the nervous system and brain were significant in pathway analysis (FDR q < 0.05), but these associations were not significant once NDD regions were removed. A larger sample size, better detection methods, and annotated resources of CNV are needed to explore this relationship further.

Molecular evolutionary process of advanced gastric cancer during sequential chemotherapy detected by circulating tumor DNA

BackgroundEfficacy of conventional sequential chemotherapy paradigm for advanced gastric cancer (AGC) patients has largely plateaued. Dynamic molecular changes during and after sequential chemotherapy have not been fully delineated. We aimed to profile the molecular evolutionary process of AGC patients during sequential chemotherapy by next generation sequencing (NGS) of plasma circulating tumor DNA (ctDNA).MethodsA total of 30 chemo-naïve patients who were diagnosed with unresectable advanced or metastatic stomach adenocarcinoma were enrolled. All patients received sequential chemotherapy regimens following the clinical guideline. One hundred and eight serial peripheral blood samples were collected at baseline, radiographical assessment and disease progression. Plasma ctDNA was isolated and a customized NGS panel was used to detect the genomic features of ctDNA including single nucleotide variants (SNVs) and gene-level copy number variations (CNVs). KEGG pathway enrichment analysis was performed.ResultsPlatinum-based combination chemotherapy was administrated as first-line regimen. Objective response rate was 50% (15/30). Patients with higher baseline values of copy number instability (CNI), CNVs and variant allel frequency (VAF) were more sensitive to platinum-based first-line regimens. Tumor mutation burden (TMB), CNI and CNV burden at partial response and stable disease were significantly lower than those at baseline, where at progressive disease they recovered to baseline levels. Dynamic change of TMB (ΔTMB) was correlated with progression-free survival of first-line treatment. Fluctuating changes of SNVs and gene-level CNVs could be observed during sequential chemotherapy. Under the pressure of conventional chemotherapy, the number of novel gene-level CNVs were found to be higher than that of novel SNVs. Such novel molecular alterations could be enriched into multiple common oncologic signaling pathways, including EGFR tyrosine kinase inhibitor resistance and platinum drug resistance pathways, where their distributions were found to be highly heterogenous among patients. The impact of subsequent regimens, including paclitaxel-based and irinotecan-based regimens, on the molecular changes driven by first-line therapy was subtle.ConclusionBaseline and dynamic changes of genomic features of ctDNA could be biomarkers for predicting response of platinum-based first-line chemotherapy in AGC patients. After treatment with standard chemotherapy regimens, convergent oncologic pathway enrichment was identified, which is yet characterized by inter-patient heterogenous gene-level CNVs.

The cross-talk of cancer-associated fibroblasts assist in prognosis and immunotherapy in patients with breast carcinoma.

The association between cancer-associated fibroblasts (CAFs) and tumor microenvironment (TME) is a key factor in promoting tumor progression. However, the correlation between CAFs and TME in breast carcinoma has not been elucidated. Thus, further study about the cross-effect between CAFs and TME can provide novel strategies for breast carcinoma treatment, particularly targeted immunotherapy. First, we systematically analyzed cell communication in a single-cell dataset and identified the interacted genes between CAFs and TME components. Then, a robust fibroblast-related score (FRS) model was developed using the LASSO algorithm. The FRS can be a reliable adverse prognostic factor in three cohorts with breast carcinoma. Functional enrichment analysis and single-sample Gene Set Enrichment Analysis showed that patients with a high FRS had cold tumors with active proliferation and immunosuppression. Patients with a low FRS presented with hot tumors with active immune and cell-killing functions. Genomic variation analysis revealed that patients with a low FRS had a higher somatic mutation load and copy number variation burden. Finally, patients with a low FRS were more sensitive to chemotherapy and immunotherapy, particularly anti-PD-1 therapy. In conclusion, a reliable FRS model was constructed not only reliable for predicting prognosis but also competent to estimate clinical immunotherapy and chemotherapy response for patients with BRCA, which might provide significant clinical implications for guiding clinical decision-making for patients with BRCA.

Exploration of Novel Immunological Terms in Lung Cancer With Large Populations: Implications for Immunotherapy.

BackgroundIdeal biomarkers to predict the response to immunotherapy in lung cancer are still lacking. Therefore, there is a need to explore effective biomarkers in large populations.ObjectiveThe objective of this study is to explore novel immunological classifications that are associated with immunotherapy response through the ssGSEA algorithm.MethodsSix independent lung cancer cohorts were collected for analysis including The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), and the EMBL-EBI database. The ssGSEA algorithm was performed to extract immune terms. Then, TCGA samples were involved as a training group and other cohorts were used as a validation group. After LASSO and Cox regression, prognostic associated immune terms were extracted and an immune-related risk score (IRS) signature was constructed. Furthermore, the association between IRS signature and clinical data, genome features, stemness indices analysis, tumor immune microenvironment, immunotherapy efficiency, and targeted therapy response was also investigated.ResultsA total of 1,997 samples were enrolled in this study including six large lung cancer cohorts. Fifty-four immune terms were calculated through the ssGSEA algorithm in combined cohorts. Then, a nine-immune-term risk score model named IRS signature was established to predict the prognosis in combined cohorts. We classified patients into high-risk and low-risk subgroups according to the cutoff point. Subsequently, analysis of clinical data and genome features indicated that the patients in the high-IRS group tend to have advanced clinical features (clinical stage and T classification), as well as a higher level of copy number variation burden, higher tumor burden mutation, and higher tumor stemness indices. Immune landscape analysis demonstrated that high-IRS groups exhibited lower immune cell infiltration and immune-suppressive state. More importantly, the predicted result of the Tumor Immune Dysfunction and Exclusion analysis showed that high-IRS groups might be more insensitive to immunotherapy. Meanwhile, we have also identified that high-IRS groups were associated with better efficiency of several targeted drugs.ConclusionTo summarize, we identified a novel IRS model based on nine immune terms, which was quantified by the ssGSEA algorithm. This model had good efficacy in predicting overall survival and immunotherapy response in non-small cell lung cancer patients, which might be an underlying biomarker.

Abstract 1452: Analyses of whole exome sequencing (WES) data of hereditary lung cancer families identify germline copy number variations (CNVs) in multiple genes

Abstract Lung cancer (LC) is the most common cause of cancer mortality and the third most common cancer by incidence in the United States. While environmental factors play an important role in its development, LC risk also exhibits a high degree of heritability. Linkage analyses and genome-wide association studies have identified multiple loci associated with increased LC risk; however, much of the heritability has yet to be explained by these loci. Structural mutations, including copy number variations (CNVs), contribute to phenotypic diversity through dosage and/or cis-regulatory effects. Unbiased CNV mapping is only possible with the use of high-depth, short-read sequencing and remains much more complicated and prone to false positives than SNP detection. As the tools necessary for CNV detection have improved, both somatic and germline CNVs have been shown to play an important role in disease, especially cancer, and often they account for a significant proportion of pathogenic variants. Little has been done to understand the role of germline CNVs in the biological pathways of hereditary lung cancer (HLC). The goal of the current project is to utilize the whole exome sequencing (WES) data in identifying the CNVs in the HLC families (≥3 LC/family) recruited by the Genetic Epidemiology of Lung Cancer Consortium (GELCC). This work uses germline WES data from 203 individuals (60 with a LC diagnosis) from 25 HLC families (≥3 LC affected/family). We limited our investigation to CNVs called by two independent tools using read depth to infer copy number changes and genomic breakpoints: GATK 4 (https://gatk.broadinstitute.org) and XHMM (https://atgu.mgh.harvard.edu/xhmm/). We have identified CNVs (deletions) in the following protein coding genes: KIR2DL1 (19q13.42; OR=1.77), DMBT1 (10q26.13; OR=1.56), SIBPB1 (20p13; OR=1.53), LILRA6 (20p13; OR=1.53), and HLA-DRB5 (6p21.32; OR=1.64) in ≥3 families in at least two individuals/family. All of these genes with the exception of HLA-DRB5 were previously reported in LC studies; HLA-DRB5 was reported in other cancer studies. DMBT1 was reported to be a candidate tumor suppressor gene and frequent loss of heterozygosity (LOH) was observed in several human tumors. It is highly expressed in lung and loss of expression was reported in LC cell lines. All CNVs identified in our study are consistent with Mendelian expectations and appear to be uncommon in the general population. CNVs in these genes may contribute to the lung cancer risk in the HLC families. Further validation of the identified variants is ongoing. We will consider overall CNV burden and disease association with identified deletions and duplications. The enrichment of rare variants in oncogenically associated genes that co-segregate with LC provides specific mutations for future work and improves our understanding of the inheritance and pathogenesis of LC. Citation Format: John Waldron, Kirsten W. Termine, Anthony M. Musolf, Mariza de Andrade, Colette Gaba, Ramaswamy Govindan, Ping Yang, Ming You, Ellen Jaeger, Angelle F. Bencaz, Marshall W. Anderson, Ann G. Schwartz, Susan M. Pinney, Christopher I. Amos, Joan E. Bailey-Wilson, Diptasri M. Mandal. Analyses of whole exome sequencing (WES) data of hereditary lung cancer families identify germline copy number variations (CNVs) in multiple genes [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 1452.