Genome-wide Copy Number Variation Research Articles

Two Genetic Mechanisms in Two Siblings with Intellectual Disability, Autism Spectrum Disorder, and Psychosis.

Intellectual disability (ID) and autism spectrum disorder (ASD) are complex neurodevelopmental disorders with high heritability. To search for the genetic deficits in two siblings affected with ID and ASD in a family, we first performed a genome-wide copy number variation (CNV) analysis using chromosomal microarray analysis (CMA). We found a 3.7 Mb microdeletion at 22q13.3 in the younger sister. This de novo microdeletion resulted in the haploinsufficiency of SHANK3 and several nearby genes involved in neurodevelopment disorders. Hence, she was diagnosed with Phelan–McDermid syndrome (PMS, OMIM#606232). We further performed whole-genome sequencing (WGS) analysis in this family. We did not detect pathogenic mutations with significant impacts on the phenotypes of the elder brother. Instead, we identified several rare, likely pathogenic variants in seven genes implicated in neurodevelopmental disorders: KLHL17, TDO2, TRRAP, EIF3F, ATP10A, DICER1, and CDH15. These variants were transmitted from his unaffected parents, indicating these variants have only moderate clinical effects. We propose that these variants worked together and led to the clinical phenotypes in the elder brother. We also suggest that the combination of multiple genes with moderate effects is part of the genetic mechanism of neurodevelopmental disorders.

Abstract 1177: Introduction of the GA4GH Variation Representation Specification (VRS) and supporting tools for discovery and exchange of clinical genomic and cytogenomic knowledge in cancers

Abstract Precision oncology is the practice of interpreting the clinical significance of observed molecular changes in patient neoplasms, potentially impacting medical decision making and care. This process is labor-intensive and (among other challenges) involves accurately translating between variation representation conventions from one resource to the next. For example, differences in representations of Copy Number Variation (CNV) from genomic regions, cytogenomic bands, or gene features create challenges in knowledge matching due to lack of standards covering all of these modalities of observed variation.The Global Alliance for Genomics and Health (GA4GH; ga4gh.org) is an international collaborative of genomic data sharing initiatives (Driver Projects) developing genomic data sharing standards within a human rights framework. GA4GH recently published the Variation Representation Specification (VRS; pronounced “verse”), a standard for the computational representation of biomolecular variation. VRS is a terminology, schema, and associated conventions for creating uniquely identifiable and federatable representations of molecular variation. VRS has formal data classes well-suited to differentiating between variation on a single molecule (e.g. tandem duplications) from variation measured at a systemic level (e.g. genome-wide copy number variation). In addition to molecular sequence variation, VRS also supports variation on cytogenetic coordinate systems and genes, making it well-suited to representing variation associated with cancer biomarkers.We demonstrate the use of VRS to model reported gene-associated CNVs from the AACR Project GENIE cohort, to aid in the computational discovery of evidence from clinico-genomic knowledgebases with genomic or cytogenomic CNV representations. We highlight the use case of knowledge matching to the Atlas of Genetics and Cytogenetics in Oncology and Haematology (“the Atlas”; atlasgeneticsoncology.org), a cytogenetics resource historically driven by user website navigation. Using VRS search tools we developed for the Variant Interpretation for Cancer Consortium (VICC; cancervariants.org) GA4GH Driver Project, we found that 64% of GENIE samples with reported CNVs matched clinically relevant knowledge in the Atlas. This work was enabled by programmatic search tools leveraging standard VRS object structures, demonstrating how VRS enables collection of real-world evidence across more resources without manual interpretation or custom normalization methods. We conclude with a survey of open-source tools supporting this analysis as well as search of other clinico-genomic knowledgebases with VRS, including CIViC (civicdb.org), BRCA Exchange (brcaexchange.org), and the Molecular Oncology Almanac (moalmanac.org). Citation Format: Matthew Cannon, Kori Kuzma, James Stevenson, Jiachen Liu, Colin O'Sullivan, Bimal P. Chaudhari, Matthew Brush, Robert R. Freimuth, Tristan Nelson, Michael Baudis, Obi L. Griffith, Malachi Griffith, Lawrence Babb, Melissa S. Cline, Xuelu Liu, Brian Walsh, Alex H. Wagner. Introduction of the GA4GH Variation Representation Specification (VRS) and supporting tools for discovery and exchange of clinical genomic and cytogenomic knowledge in cancers [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 1177.

Abstract 540: Genome-wide circulating tumor DNA for monitoring treatment response and metastatic relapse in bladder cancer

Abstract Background: Neoadjuvant chemotherapy (NAC) followed by radical cystectomy (CX), is gold standard treatment in localized muscle-invasive bladder cancer (MIBC). About 45% of patients with MIBC develop metastatic relapse within 2 years after CX. The response rate to chemotherapy and immune checkpoint inhibitors (ICI) is relatively low, and biomarker tests for monitoring response are needed. Furthermore, biomarkers for early detection of minimal residual disease (MRD) after CX is needed to enable earlier treatment initiation. Tumor-informed detection of mutations in cell-free DNA (cfDNA) from peripheral blood has shown promising results in its ability to monitor MRD. However, the low tumor fraction after surgery and limited input material obtained from a typical plasma sample limits the probability of detecting low metastatic burden scenarios. Here we implemented and applied locally a whole-genome sequencing (WGS) approach to circulating tumor DNA (ctDNA) monitoring for improving ctDNA detection. Methods: A total of 140 MIBC patients undergoing NAC and CX were enrolled, including a test cohort (n=19) and a validation cohort (n=120). cfDNA was extracted from ~1mL plasma (n=1100) and procured from longitudinal plasma sampling during NAC (response measure), pre-cystectomy (response measure), post-surgery (relapse monitoring) and during immunotherapy (ICI treatment). WGS was applied to tumor/germline pairs (coverage >30x/20x) and plasma cfDNA (>20x) facilitating detection of genome wide genomic alterations and quantification of ctDNA using the MRDetect method. Results: We developed a personalized tumor-informed WGS model by integrating genome-wide mutation and copy number variation data coupled with advanced signal processing and AI-based error suppression. Patient-specific somatic variant patterns were then used for detecting and measuring the ctDNA levels in low-input blood samples by WGS. The assay sensitivity allowed for detection of tumor fractions down to 8*10-5. Furthermore, in our test cohort of 19 patients, we detected ctDNA after CX in 7 of 8 patients with clinical relapse (88% sensitivity) and detected no ctDNA in 11 of 11 patients with no clinical relapse (100% specificity). We observed a positive lead-time for MRD-based recurrence detection compared to CT-based reccurence detection (9 months on average). The full dataset is currently being processed and will be presented at the AACR 2022 meeting. Conclusions: For precision oncology, we need to develop quantitative and non-invasive methodologies to help tailor the treatments to individual patients and monitor them for further clinical decision-making. The results indicate the clinical potential of personalized genome-wide mutation integration as an ultra-sensitive, non-invasive method for MRD detection and treatment response monitoring which could aid in clinical management of patients with bladder cancer. Citation Format: Iver Nordentoft, Karin Birkenkamp-Demtröder, Emil Christensen, Sunil Deochand, Dillon Maloney, Danielle Afterman, Tomer Lauterman, Noah Friedman, Imane Bourzgui, Nidhi Ramaraj, Zohar Donenhirsh, Ronel Veksler, Sia Viborg, Mads Agerbæk, Jørgen Bjerggaard Jensen, Jonathan Rosenfeld, Ravi Kandasamy, Iman Tavassoly, Boris Oklander, Asaf Zviran, Lars Dyrskjøt. Genome-wide circulating tumor DNA for monitoring treatment response and metastatic relapse in bladder cancer [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 540.

Abstract 2171: The spatial landscape of clonal somatic mutations in benign and malignant tissue

Abstract Introduction: Defining the transition from benign to malignant tissue is fundamental to improve early diagnosis of cancer. In order to obtain spatial information of clonal genetic events, prior studies have used methods such as laser capture microdissection, which results in assessment of small regions or even single cells. These studies have an inherent bias as a limited number of regions per tissue section can be retrieved and examined. Furthermore, because investigators have selected such regions based on morphology, previous studies have limited their analyses to histologically defined tumour areas while excluding regions populated by benign cells. The possibility to perform unsupervised genome and tissue-wide analysis would therefore provide an important contribution to delineate clonal events. We sought study spatial genome integrity in situ to gain molecular insight into clonal relationships. Materials and Methods: We employed spatially resolved transcriptomics (Visium, 10x Genomics) to infer spatial copy number variations in >120 000 spatial regions across multiple organs, including three whole axial prostates and additional tissues from skin, breast and brain tumors. We used this information to deduce clonal relationships between regions harboring 5-20 cells. Results: We demonstrate that genome-wide copy number variation reveals distinct clonal patterns within tumours and in nearby benign tissue. We perform an in-depth spatial analysis of cancers that includes an unprecedented interrogation of up to 50,000 tissue domains in a single patient, and 120,000 tissue domains across 10 patients. In a prostate section, we observed that many CNVs occurred in histologically benign luminal epithelial cells, most notably in chromosomes 8 and 10. This clone constituted a region of exclusively benign acinar cells branching off a duct lined by largely copy neutral cells. The changes in these cells were shared with the nearby intermediate risk prostate cancer cells in the same tissue section. We observed similar findings in another patient’s cutaneous squamous cell carcinoma (cSCC), wherein benign squamous epithelial had alterations in chromosomes 1 and 12 that were shared with nearby cSCC. Our results suggest a model for how genomic instability arises in histo-pathologically benign tissue that may represent early events in cancer evolution. Furthermore the spatial information allowed us to identify small clonal units not evident from morphology and hence would be overlooked by pathologists. Conclusions: We present the first large-scale, comprehensive atlas of genomic evolution at high spatial resolution in prostate cancer. Our study adds an important new approach to the armamentarium of cancer molecular pathology. We highlight the power of an unsupervised approach to capture the molecular and spatial continuums in a tissue context and challenge the rationale for focal therapy in prostate cancer. Citation Format: Andrew Erickson, Emelie Berglund, Mengxiao He, Maja Marklund, Reza Mirzazadeh, Niklas Schultz, Linda Kvastad, Alma Andersson, Ludvig Bergenstråhle, Joseph Bergenstråhle, Ludvig Larsson, Alia Shamikh, Elisa Basmaci, Teresita Diaz De Ståhl, Timothy Rajakumar, Kim Thrane, Andrew L. Ji, Paul A. Khavari, Firaz Tarish, Anna Tanoglidi, Jonas Maaskola, Richard Colling, Tuomas Mirtti, Freddie Hamdy, Dan J. Woodcock, Thomas Helleday, Ian G. Mills, Alastair D. Lamb, Joakim Lundenberg. The spatial landscape of clonal somatic mutations in benign and malignant tissue [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 2171.

Identification of MMACHC and ZEB2 mutations causing coexistent cobalamin C disease and Mowat-Wilson syndrome in a 2-year-old girl

Cobalamin C (cblC) disease and Mowat-Wilson syndrome (MWS) are rare hereditary diseases. To date, there have been no reports of people suffering from these two genetic diseases, or whether there is any correlation between the two diseases. We reported a 2-year-old girl with both cblC disease and MWS. The patient initially manifested as slow weight gain, hypotonia, broad nasal bridge, high forehead, high palate arch, ear crease, patent ductus arteriosus, atrial and ventricular septal defect and bilateral mild ventriculomegaly in the neonatal period. However, as the baby grew older, the typical facial features became more prominent, and overall developmental delays were noted at the subsequent follow-up, with the motor and cognitive development significantly lagging behind that of other children of the same age. At 26 days old, laboratory tests revealed remarkably elevated levels of serum homocysteine, C3/C2 and urine organic acid. Whole-exome sequencing detected compound heterozygous variants in MMACHC, including one previously reported mutation [c.609G > A (p.W203X) and a novel missense mutation[ c.643 T > C (p.Y215H)]. The computer simulations of the protein structure analysis of the novel missense mutation showed the variant p.Y215H replaced a neutral amino acid with a strongly basic lysine, which broken the local structure by changing the carbon chain skeleton and decreasing the interaction with adjacent amino acid. This is expected to damage the utilization of vitamin B12 and influence the synthesis of AdoCbl and MeCbl, contributing to its pathogenicity. Thus, clinical and genetic examinations confirmed the cblC disease. Another heterozygous variant in ZEB2 [NM_014795; loss1(exon:2–10)(all); 127901 bp] was detected by whole-exome sequencing. The heterozygous 3.04 Mb deletion in EB2 [GRCH37]del(2)(q22.2q22.3) (chr2:142237964–145274917) was also confirmed by genome-wide copy number variations (CNVs) scan, which was pathogenic and led to the diagnosis of Mowat-Wilson syndrome. The biochemical indicators associated with cblC disease in the patient were well controlled after treatment with vitamin B12 and betaine. Here, a patient with coexisting cblC disease and MWS caused by different pathogenic genes was reported, which enriched the clinical research on these two rare genetic diseases.

Analytical evaluation of a commercial in-house homologous recombination deficiency (HRD) assay for patients with epithelial ovarian cancers and its concordance with a reference standard.

e17592 Background: Homologous recombination deficiency (HRD) is a common feature of epithelial ovarian cancers (EOCs) as approximately half of EOCs have molecular aberrations in homologous recombination repair (HRR) genes, including BRCA1/2 gene mutations. Poly ADP-ribose Polymerase inhibitors (PARPi) have reformed the management of BRCA1/2 mutant ovarian cancer patients and demonstrated activity also in tumor with HRD status. The evaluation of HRD together with BRCA1/2 gene alterations is nowadays critical to guide the therapy decision for all patients with EOC. Different methods to assessed HRD have been reported and different assays are being implemented in clinical practice. Here, we sought to evaluate the analytical validity of a commercial in house HRD assay as compared to a reference standard method. Methods: 59 consecutive case of EOCs were included in this study. DNA extracted from representative formalin-fixed paraffin-embedded (FFPE) tissue block of each case were subjected to next-generation sequencing (NGS) analysis of BRCA1/2 genes alterations and HRD status using the AmoyDx HRD Focus Panel (Amoy Diagnostics Ltd, Xiamen, Fujian, China). HRD status was assesses through AmoyDx Genomic Scar Score (GSS) considering the length, type and location of the genome-wide chromosomal copy number variation through SNP probes spread throughout the genome. Cases with a GSS > 50 were considered as HRD positive. For each case the same FFPE tissue block was sent to Myriad Genetic Laboratories for myChoice CDx testing as reference standard. Assay concordance was measured as positive (PPA), negative (NPA) and overall (OPA) percentage of agreement according to the different scores proposed by the manufacturers. Results: HRD status data were available for 52 (88.1%) and 56 (94.9%) samples considering in-house and reference standard assay, respectively, and 49 cases was available for concordance analyses. AmoyDx HRD Focus Panel identified 23/49 samples (46.9%) as HRD positive and 26 (53.1%) as HRD negative. The concordance of HRD status between AmoyDx HRD Focus Panel and myChoice CDx was 100% and 81,2% for positive (PPA) and negative (NPA) cases respectively, with an overall analytical concordance (OPA) of 87.8%. No discordance was seen in BRCA1/2 status between in-house and reference standard assay (OPA: 100%). Conclusions: The HRD status assessed by AmoyDx HRD Focus Panel is highly concordant with that obtained using the current standard reference method. Commercial in-house test may offer a cost-effective, reliable and robust alternative to the centralized assay for the evaluation of HRD status in patient with EOC.

Somatic copy number variant mutations in alpha-synuclein and genome-wide in brains of synucleinopathy cases

Synucleinopathies are mostly sporadic neurodegenerative disorders, and include Parkinson’s disease (PD), and multiple system atrophy (MSA). Inherited copy number variants (CNVs) ofSNCA(α-synuclein) are rare causes of familial disease. There is increasing evidence for somatic mutations in the human brain, and we have hypothesized a role for these in synucleinopathies.We further investigated somatic CNVs in brains from synucleinopathy cases, using Fluorescent in-situ Hybri- disation forSNCAgains, and single-cell whole genome sequencing for the first time in a synucleinopathy. In the cingulate cortex, mosaicism forSNCAgains was higher in MSA and PD than controls in neurons (> 2% in both diseases), and for MSA also in non-neurons. In MSA substantia nigra (SN), we notedSNCAgains in> 3% of cells. Cells with CNVs were more likely to have α-synuclein inclusions, in a pattern corresponding to cell categories mostly relevant to the disease: dopaminergic neurons in Lewy-body cases, and other cells in striatonigral degeneration-dominant MSA (MSA-SND). HigherSNCACNV mosaicism in typical MSA-SND SN non-dopaminergic cells may correlate with younger onset, and in PD cingulate neurons with younger death. Whole genome sequencing of 169 single cells from two MSA cases showed somatic CNVs in ~ 30%. Neurons had gains and losses, and other cells almost exclusively gains.We propose that somaticSNCACNVs contribute to the aetiology and pathogenesis of synucleinopathies, and that genome-wide somatic CNVs merit further studyc.proukakis@ucl.ac.uk

Abstract PR016: The spatial landscape of clonal somatic mutations in benign and malignant tissue

Abstract Introduction: Defining the transition from benign to malignant tissue is fundamental to improve early diagnosis of cancer. In order to obtain spatial information of clonal genetic events, prior studies have used methods such as laser capture microdissection, which results in assessment of small regions or even single cells. These studies have an inherent bias as a limited number of regions per tissue section can be retrieved and examined. Furthermore, because investigators have selected such regions based on morphology, previous studies have limited their analyses to histologically defined tumor areas while excluding regions populated by benign cells. The possibility to perform unsupervised genome and tissue-wide analysis would therefore provide an important contribution to delineate clonal events. We sought study spatial genome integrity in situ to gain molecular insight into clonal relationships. Materials and Methods: We employed spatially resolved transcriptomics (Visium, 10x Genomics) to infer spatial copy number variations in >120 000 spatial regions across multiple organs, including three whole axial prostates and additional tissues from skin, breast and brain tumors. We additionally performed in silico assessment of spatial copy number inference. We used this information to deduce clonal relationships between regions harboring 5-20 cells. Results: We demonstrate that genome-wide copy number variation reveals distinct clonal patterns within tumors and in nearby benign tissue. We perform an in-depth spatial analysis of cancers that includes an unprecedented interrogation of up to 50,000 tissue domains in a single patient, and 120,000 tissue domains across 10 patients. In a prostate section, we observed that many CNVs occurred in histologically benign luminal epithelial cells, most notably in chromosomes 8 and 10. This clone constituted a region of exclusively benign acinar cells branching off a duct lined by largely copy neutral cells. The changes in these cells were shared with the nearby intermediate risk prostate cancer cells in the same tissue section. We observed similar findings in another patient’s cutaneous squamous cell carcinoma (cSCC), wherein benign squamous epithelial had alterations in chromosomes 1 and 12 that were shared with nearby cSCC. Our results suggest a model for how genomic instability arises in histo-pathologically benign tissue that may represent early events in cancer evolution. Furthermore the spatial information allowed us to identify small clonal units not evident from morphology and hence would be overlooked by pathologists. Conclusions: We present the first large-scale, comprehensive atlas of genomic evolution at high spatial resolution in prostate cancer. Our study adds an important new approach to the armamentarium of cancer molecular pathology. We highlight the power of an unsupervised approach to capture the molecular and spatial continuums in a tissue context and challenge the rationale for focal therapy in prostate cancer Citation Format: Andrew Erickson, Emelie Berglund, Mengxiao He, Maja Marklund, Reza Mirzazadeh, Niklas Schultz, Ludvig Bergenstråhle, Linda Kvastad, Alma Andersson, Joseph Bergenstråhle, Ludvig Larsson, Alia Shamikh, Elisa Basmaci, Teresita Diaz De Ståhl, Timothy Rajakumar, Kim Thrane, Andrew L. Ji, Paul A. Khavari, Firaz Tarish, Anna Tanoglidi, Jonas Maaskola, Richard Colling, Tuomas Mirtti, Freddie C. Hamdy, Dan J. Woodcock, Thomas Helleday, Ian G. Mills, Alastair D. Lamb, Joakim Lundeberg. The spatial landscape of clonal somatic mutations in benign and malignant tissue [abstract]. In: Proceedings of the AACR Special Conference on the Evolutionary Dynamics in Carcinogenesis and Response to Therapy; 2022 Mar 14-17. Philadelphia (PA): AACR; Cancer Res 2022;82(10 Suppl):Abstract nr PR016.

MP45-14 THE SPATIAL LANDSCAPE OF CLONAL SOMATIC MUTATIONS IN BENIGN AND MALIGNANT PROSTATE EPITHELIA

MP45-14 THE SPATIAL LANDSCAPE OF CLONAL SOMATIC MUTATIONS IN BENIGN AND MALIGNANT PROSTATE EPITHELIA

Genome-wide elucidation of CNV regions and their association with production and reproduction traits in composite Vrindavani cattle

The present study was aimed to analyze the genome-wide copy number variations (CNVs) in Vrindavani composite cattle and concatenate them into CNV regions (CNVRs), and finally test the association of CNVRs with different production and reproduction traits. Genotypic data, generated on BovineSNP50 Beadchip (v3) array for 96 Vrindavani animals, was used to elucidate the CNVs at the genome level. Intensity data covering over 53,218 SNP genotypes on bovine genome was used. Algorithm based on Hidden Markov Model was employed in PennCNV program to detect, normalize and filter CNVs across the genome. 252 putative CNVs, detected via PennCNV program, in different individuals were concatenated into 71 CNV regions (CNVRs) using CNVRuler program. Association of CNVRs with important (re)production traits in Vrindavani animals was assessed using linear regression. Five CNVRs were found to be significantly associated with ten important (re)production traits. The genes harbored in these regions provided useful insights into the association of CNVRs with genes and ultimately the variation at phenotype level. Important genes that overlapped with CNVRs included WASHC4, HS6ST3, MBNL2, TOLLIP, PIDD1 and TSPAN4. Furthermore, the CNVRs were found to overlap with important QTLs available in AnimalQTL database which affect milk yield and composition along with reproduction and immune function traits. The copy number states of three enes were validated using digital droplet PCR technique. The results from the present study significantly enhance the understanding about CNVs in Vrindavani cattle and should help establish its CNV map. The study will also enable further investigation on association of these variants with important traits of economic interest including disease incidence.

Schizophrenia Risk Mediated by microRNA Target Genes Overlapped by Genome-Wide Rare Copy Number Variation in 22q11.2 Deletion Syndrome.

The 22q11.2 deletion is associated with >20-fold increased risk for schizophrenia. The presence of gene DGCR8 in the 22q11.2 deletion region has suggested microRNA (miRNA) dysregulation as possibly contributing to this risk. We therefore investigated the role of miRNA target genes in the context of previously identified genome-wide risk for schizophrenia conveyed by additional copy number variation (CNV) in 22q11.2 deletion syndrome (22q11.2DS). Using a cohort of individuals with 22q11.2DS and documented additional rare CNVs overlapping protein coding genes, we compared those with schizophrenia (n = 100) to those with no psychotic illness (n = 118), assessing for rare CNVs that overlapped experimentally supported miRNA target genes. We further characterized the contributing miRNA target genes using gene set enrichment analyses and identified the miRNAs most implicated. Consistent with our hypothesis, we found a significantly higher proportion of individuals in the schizophrenia than in the non-psychotic group to have an additional rare CNV that overlapped one or more miRNA target genes (odds ratio = 2.12, p = 0.0138). Gene set analyses identified an enrichment of FMRP targets and genes involved in nervous system development and postsynaptic density amongst these miRNA target genes in the schizophrenia group. The miRNAs most implicated included miR-17-5p, miR-34a-5p and miR-124-3p. These results provide initial correlational evidence in support of a possible role for miRNA perturbation involving genes affected by rare genome-wide CNVs in the elevated risk for schizophrenia in 22q11.2DS, consistent with the multi-hit and multi-layered genetic mechanisms implicated in this and other forms of schizophrenia.

Identification of rare mutations of the vasoactive intestinal peptide receptor 2 gene in schizophrenia.

Studies showed that rare copy number variations (CNVs) encompassing the vasoactive intestinal peptide receptor 2 gene (VIPR2) were associated with schizophrenia, indicating VIPR2 is a risk gene for schizophrenia. We hypothesized that besides CNV, rare pathogenic single-nucleotide variant (SNV) or small insertion/deletion (Indel) of VIPR2 might be present in some patients and contribute to the pathogenesis of schizophrenia. We performed genome-wide CNV analysis to screen CNV at the VIPR2 locus and targeted sequencing of all the exons of VIPR2 to search for SNV and indel in a sample of patients with chronic schizophrenia from Taiwan. We detected a 230-kb microduplication encompassing the VIPR2 in 1 out of 200 patients. Furthermore, we identified six ultrarare SNVs, including one splicing SNV and five missense SNVs, in 516 patients. In-silico analyses showed these SNVs had a damaging effect on the function of VIPR2. Our findings support the idea that besides CNV, rare pathogenic SNVs of VIPR2 might contribute to the pathogenesis of schizophrenia in some patients.

SavvyCNV: Genome-wide CNV calling from off-targetreads.

Identifying copy number variants (CNVs) can provide diagnoses to patients and provide important biological insights into human health and disease. Current exome and targeted sequencing approaches cannot detect clinically and biologically-relevant CNVs outside their target area. We present SavvyCNV, a tool which uses off-target read data from exome and targeted sequencing data to call germline CNVs genome-wide. Up to 70% of sequencing reads from exome and targeted sequencing fall outside the targeted regions. We have developed a new tool, SavvyCNV, to exploit this 'free data' to call CNVs across the genome. We benchmarked SavvyCNV against five state-of-the-art CNV callers using truth sets generated from genome sequencing data and Multiplex Ligation-dependent Probe Amplification assays. SavvyCNV called CNVs with high precision and recall, outperforming the five other tools at calling CNVs genome-wide, using off-target or on-target reads from targeted panel and exome sequencing. We then applied SavvyCNV to clinical samples sequenced using a targeted panel and were able to call previously undetected clinically-relevant CNVs, highlighting the utility of this tool within the diagnostic setting. SavvyCNV outperforms existing tools for calling CNVs from off-target reads. It can call CNVs genome-wide from targeted panel and exome data, increasing the utility and diagnostic yield of these tests. SavvyCNV is freely available at https://github.com/rdemolgen/SavvySuite.

Guiding the global evolution of cytogenetic testing for hematologic malignancies

AbstractCytogenetics has long represented a critical component in the clinical evaluation of hematologic malignancies. Chromosome banding studies provide a simultaneous snapshot of genome-wide copy number and structural variation, which have been shown to drive tumorigenesis, define diseases, and guide treatment. Technological innovations in sequencing have ushered in our present-day clinical genomics era. With recent publications highlighting novel sequencing technologies as alternatives to conventional cytogenetic approaches, we, an international consortium of laboratory geneticists, pathologists, and oncologists, describe herein the advantages and limitations of both conventional chromosome banding and novel sequencing technologies and share our considerations on crucial next steps to implement these novel technologies in the global clinical setting for a more accurate cytogenetic evaluation, which may provide improved diagnosis and treatment management. Considering the clinical, logistic, technical, and financial implications, we provide points to consider for the global evolution of cytogenetic testing.

Abstract CC07-01: Whole exome and whole genome methylation sequencing of low-input cfDNA to implement precision medicine in metastatic castration resistant prostate cancer

Abstract Purpose Liquid biopsy has become increasingly important in cancer diagnosis, personalized medicine, and disease progression monitoring. Conventional liquid biopsy relies on targeted cancer gene panels which often contain fewer than 500 genes. Despite the revolutionary impact of liquid biopsy on cancer research and patient care, targeted gene panels may miss key novel mutations involved in cancer development and drug response, and other novel genomic alternations underlining cancer development, such as whole genome structural or epigenetic changes. Additionally, current commonly used methodologies require relatively large amounts of cell free DNA (cfDNA) input material, thus limiting its application when only small amounts of cfDNA input is available. Experimental Design cfDNA was isolated from whole blood collected from commercially available prostate cancer patients or age-matched healthy donors and evaluated using Predicine’s platform for whole exome sequencing, plasma WES, and compared to the targeted 500 gene PredicineATLAS panel (cfDNA input 10 or 30 ng). Additionally, low-input cfDNA (cfDNA input 10 or 20 ng) was also evaluated for whole genome methylation sequencing using PredicineECM and compared with standard whole genome bisulfite sequencing (WGBS). Results Plasma WES profiling of 30 or 10 ng cfDNA detected variants with mutation allele frequency (MAF) of 1% and above in the entire exome region, which enabled more accurate genome wide copy number variation prediction. All mutations detected using the PredicineATLAS panel that were ≥1% allele frequency were also detected by plasma WES, although low input required greater depth (10 ng, 3,500x) compared with high input (30 ng, 2500x depth). Furthermore, plasma WES identified 1.5 to 15.3x additional mutations compared with the PredicineATLAS panel. PredicineECM methylation analysis was superior to WGBS in reducing DNA damage and GC bias, resulting in increased NGS read mapping rate and quality score. Profiling of 2.5, 5 and 10 ng cfDNA using PredicineECM showed higher mapping rate and quality compared with WGBS profiling of 10 ng cfDNA. Additionally, clustering analysis of CpG methylation using PredicineECM separated tumor samples from normal samples and may serve as additional markers for minimal residue disease (MRD) tracing or early cancer detection. Conclusion In conclusion, we have successfully developed a multi-omics analysis platform for low-input cfDNA that integrates whole exome mutation and whole genome methylation profiling to obtain a comprehensive picture of cancer genomics to enable precision medicine in prostate cancer. This technology can also be applied to other cancers. Citation Format: Raju Kandimalla, Binggang Xiang, Chao Dai, Amy Wang, Pan Du, Shujun Luo, Debbie Liao. Whole exome and whole genome methylation sequencing of low-input cfDNA to implement precision medicine in metastatic castration resistant prostate cancer [abstract]. In: Proceedings of the AACR-NCI-EORTC Virtual International Conference on Molecular Targets and Cancer Therapeutics; 2021 Oct 7-10. Philadelphia (PA): AACR; Mol Cancer Ther 2021;20(12 Suppl):Abstract nr CC07-01.

26P Genome-wide copy number variation of circulating cell-free DNA as a biomarker in head and neck cancer patients treated with immunotherapy

26P Genome-wide copy number variation of circulating cell-free DNA as a biomarker in head and neck cancer patients treated with immunotherapy

Comparative analysis of genome-wide copy number variations between Tibetan sheep and White Suffolk sheep

The DNA copy number variations (CNVs) are widely involved in affecting various kinds of biological functions, such as environmental adaptation. Tibetan sheep and White Suffolk sheep are two representative indigenous and exotic breeds raised in Sichuan, China, and both of them have many contrasting biological characteristics. In this study, we employed high-throughput sequencing approach to investigate genome-wide CNVs between the two sheep breeds. A total of 11,135 CNV regions (CNVRs) consisting of 6,488 deletions and 4,647 duplications were detected, whose length ranged from 1,599 bp to 0.56 Mb with the mean of 4,658 bp. There were 281 CNVRs segregated between Tibetan sheep and White Suffolk sheep, and 18 of them have been fixed within both breeds. Functional analyses of candidate genes within the segregating CNVRs revealed the thyroid hormone signaling pathway and CTNNB1 gene that would be responsible for differential biological characteristics of breeds, such as energy metabolism, seasonal reproduction, and litter size. Furthermore, the segregating CNVRs identified in this study were overlapped with many known quantitative trait loci that are associated with growth, testis weight, and reproductive seasonality. In conclusion, these results help us better understanding differential biological characteristics between Tibetan sheep and White Suffolk sheep.

Whole Transcriptome Analysis Identifies Distinct Gene Expression Profiles between SF3B1mut and SF3B1 wt Myelodysplastic Syndrome with Ring Sideroblasts

Background and aims - The 2016 revised WHO classification incorporated somatic mutation in SF3B1 spliceosome gene within the diagnostic criteria of myelodysplastic syndrome (MDS) with ring sideroblasts (RS). However, SF3B1wt MDS-RS display significantly different clinical features and outcome from those of SF3B1mut MDS-RS. Recently, the recognition of SF3B1-mutant MDS as a distinct nosologic entity has been proposed to overcome this limitation.Methods - To evaluate the biological relevance of this proposal, we studied a consecutive cohort of 132 MDS with RS >5% using a pangenomic approach (targeted-DNA sequencing, genome-wide copy number variation analysis and bulk RNA-sequencing of CD34+ bone marrow mononuclear cells). 16 age-matched healthy individuals and 43 MDS-SLD/MLD negative for both splicing mutation and RS were included in this study as controls.Results - Unsupervised clustering analysis based on mutation profiles identified two major clusters predicted by SF3B1 mutation (87 MDS-RS-SF3B1mut and 45 MDS-RS-SF3B1wt). The most recurrently mutated genes in MDS-RS-SF3B1wt were TP53(40%), SRSF2(38%), TET2 (33%), ASXL1 (21%) and DNMT3A (12%). SRSF2 and TP53 mutations were found to be mutually exclusive with SF3B1 (p-value <0.05), whereas no difference was found in TET2, DNMT3A and ASXL1 frequencies between MDS-RS-SF3B1 mutand MDS-RS-SF3B1wt subgroups. TP53-mutated MDS-RS exhibited shorter overall survival (median 1.3 years, log-rank p-value<0.0001) compared to SF3B1mut (median 7.6 years), SRSF2mut (median 3.4 years) and MDS-RS without the abovementioned aberrations (median 4.4 years). Notably, allelic imbalances analysis of oncogenic variants identified a significant enrichment of TP53 biallelic inactivation in MDS-RS-SF3B1wt (64% vs 0, p<0.01).Differential gene expression analysis results were incorporated into a specific expression signature highly predictive of MDS-RS-SF3B1 mutand MDS-RS-SF3B1wt subgroups (Figure 1). The resulting gene clusters were classified in RS-specific genes (cluster 1 and 2) and SF3B1-specific genes (cluster 3 and 4). RS-specific genes comprising heme and hypoxia genes were enriched (Figure 2AB) and correlated with RS percentage (p<0.01). Pathway analysis revealed a specific downregulation of adhesion molecules and an upregulation of G-protein coupled receptor signaling molecules in MDS-RS-SF3B1mut. Among SF3B1-specific genes (Figure 2CD), we confirmed ABCB7 downregulation and identified new molecular targets that may concur to the pathophysiology of SF3B1-driven myeloid neoplasms.Conclusions - This study contributes to unveil molecular features of SF3B1-mutant MDS and provides further evidence to support recognition of somatic SF3B1 mutation as a disease-defining genetic lesion. DisclosuresPapaemmanuil: Isabl Technologies: Divested equity in a private or publicly-traded company in the past 24 months; Kyowa Hakko Kirin Pharma: Consultancy.

Genome-wide CNV analysis reveals variants associated with high-altitude adaptation and meat traits in Qaidam cattle

Qaidam cattle are local breeds that habitats in northwest China. It has many excellent characteristics, such as high cold and roughage tolerance, low oxygen adaptability, and tender meat quality. Copy number variation (CNV) can induce phenotypic changes in animals by a variety of effects, and thus affects the biological functions of the animals. To explore the molecular mechanism of its adaptation to extreme cold weather and muscle fat development, the CNV variations in the genome of three Qaidam cattle were detected by whole-genome sequencing, in this study. A total of 16,743 CNVs and 9498 copy number variable regions (CNVRs) were obtained after the screening, which accounts for 2.18% of the bovine genome. The CNVR length detected ranged from 0.3 KB to 10.77 KB, with a total length of 58.17 MB and an average length of 6.12 KB/ CNVR. Through functional enrichment of CNVR related genes, LDHB, and ME1 genes were screened as the key genes for Qaidam cattle to adapt to the cold and low oxygen environments, whereas KIT and FGF18 genes might be related to the coat color and growth. In the CNVR overlapped with QTLs, variation in CAPN1 and CAST genes might be closely related to the tender meat quality of Qaidam cattle. Therefore, this study provides new genetic insights on the environmental adaptability and important economic traits of Qaidam cattle. How to cite: Guo S, Wu X, Pei J, et al. Genome-wide CNV analysis reveals variants associated with high-altitude adaptation and meat traits in Qaidam cattle. Electron J Biotechnol 2021;54. https://doi.org/10.1016/j.ejbt.2021.07.006

Identifying a cervical cancer survival signature based on mRNA expression and genome-wide copy number variations.

Cervical cancer mortality is the second highest in gynecological cancers. This study developed a new model based on copy number variation data and mRNA data for overall survival prediction of cervical cancer. Differentially expressed genes from The Cancer Genome Atlas dataset detected by univariate Cox regression analysis were further simplified to six by least absolute shrinkage and selection operator (Lasso) and stepwise Akaike information criterion (stepAIC). The study developed a six-gene signature, which was further verified in independent dataset. Association between immune infiltration and risk score was investigated by immune score. The relation between the signature and functional pathways was examined by gene set enrichment analysis. Ninety-nine differentially expressed genes were detected, and C11orf80, FOXP3, GSN, HCCS, PGAM5, and RIBC2 were identified as key genes to construct a six-gene signature. The prognostic signature showed a significant correlation with overall survival (hazard ratio, HR = 3.45, 95% confidence interval (CI) = 2.08-5.72, p < 0.00001). Immune score showed a negative correlation with the risk score calculated by the signature (p < 0.05). Four immune-related pathways were closely associated with risk score (p < 0.0001). The six-gene prognostic signature was an effective tool to predict overall survival of cervical cancer. In conclusion, the newly identified six genes may be considered as new drug targets for cervical cancer treatment.