A comparison of genomic copy number calls by Partek Genomics Suite, Genotyping Console and Birdsuite algorithms to quantitative PCR

CINmetrics: an R package for analyzing copy number aberrations as a measure of chromosomal instability

Genomic instability is an important hallmark of cancer and more recently has been identified in others like neurodegenrative diseases. Chromosomal instability, as a measure of genomic instability, has been used to characterize clinical and biological phenotypes associated with these diseases by measuring structural and numerical chromosomal alterations. There have been multiple chromosomal instability scores developed across many studies in the literature; however, these scores have not been compared because of the lack of a single tool available to calculate and facilitate these various metrics. Here, we provide an R package CINmetrics, that calculates six different chromosomal instability scores and allows direct comparison between them. We also demonstrate how these scores differ by applying CINmetrics to breast cancer data from The Cancer Genome Atlas (TCGA). The package is available on CRAN at https://cran.rproject.org/package=CINmetrics and on GitHub at https://github.com/lasseignelab/CINmetrics.

GSTM3 variant is a novel genetic modifier in Brugada syndrome, a disease with risk of sudden cardiac death

BACKGROUND

Brugada syndrome (BrS) is a rare inherited disease causing sudden cardiac death (SCD). Copy number variants (CNVs) can contribute to disease susceptibility, but their role in Brugada syndrome (BrS) is unknown. We aimed to identify a CNV associated with BrS and elucidated its clinical implications.

METHODS

We enrolled 335 unrelated BrS patients from 2000 to 2018 in the Taiwanese population. Microarray and exome sequencing were used for discovery phase whereas Sanger sequencing was used for the validation phase. HEK cells and zebrafish were used to characterize the function of the CNV variant.

FINDINGS

A copy number deletion of GSTM3 (chr1:109737011-109737301, hg38) containing the eighth exon and the transcription stop codon was observed in 23.9% of BrS patients versus 0.8% of 15,829 controls in Taiwan Biobank (P < 0.001), and 0% in gnomAD. Co-segregation analysis showed that the co-segregation rate was 20%. Patch clamp experiments showed that in an oxidative stress environment, GSTM3 down-regulation leads to a significant decrease of cardiac sodium channel current amplitude. Ventricular arrhythmia incidence was significantly greater in gstm3 knockout zebrafish at baseline and after flecainide, but was reduced after quinidine, consistent with clinical observations. BrS patients carrying the GSTM3 deletion had higher rates of sudden cardiac arrest and syncope compared to those without (OR: 3.18 (1.77-5.74), P<0.001; OR: 1.76 (1.02-3.05), P = 0.04, respectively).

INTERPRETATION

This GSTM3 deletion is frequently observed in BrS patients and is associated with reduced I_Na, pointing to this as a novel potential genetic modifier/risk predictor for the development of the electrocardiographic and arrhythmic manifestations of BrS.

FUNDING

This work was supported by the Ministry of Science and Technology (107-2314-B-002-261-MY3 to J.M.J. Juang), and by grants HL47678, HL138103 and HL152201 from the National Institutes of Health to CA.

Segregating patterns of copy number variations in extended autism spectrum disorder (ASD) pedigrees

Autism spectrum disorder (ASD) is a relatively common childhood onset neurodevelopmental disorder with a complex genetic etiology. While progress has been made in identifying the de novo mutational landscape of ASD, the genetic factors that underpin the ASD's tendency to run in families are not well understood. In this study, nine extended pedigrees each with three or more individuals with ASD, and others with a lesser autism phenotype, were phenotyped and genotyped in an attempt to identify heritable copy number variants (CNVs). Although these families have previously generated linkage signals, no rare CNV segregated with these signals in any family. A small number of clinically relevant CNVs were identified. Only one CNV was identified that segregated with ASD phenotype; namely, a duplication overlapping DLGAP2 in three male offspring each with an ASD diagnosis. This gene encodes a synaptic scaffolding protein, part of a group of proteins known to be pathologically implicated in ASD. On the whole, however, the heritable nature of ASD in the families studied remains poorly understood.

A Genome-Wide Association Study to Identify Potential Germline Copy Number Variants for Sporadic Breast Cancer Susceptibility

Breast cancer (BC) predisposition in populations arises from both genetic and nongenetic risk factors. Structural variations such as copy number variations (CNVs) are heritable determinants for disease susceptibility. The primary objectives of this study are (1) to identify CNVs associated with sporadic BC using a genome-wide association study (GWAS) design; (2) to utilize 2 distinct CNV calling algorithms to identify concordant CNVs as a strategy to reduce false positive associations in the hypothesis-generating GWAS discovery phase, and (3) to identify potential candidate CNVs for follow-up replication studies. We used Affymetrix SNP Array 6.0 data profiled on Caucasian subjects (422 cases/348 controls) to call CNVs using algorithms implemented in Nexus Copy Number and Partek Genomics Suite software. Nexus algorithm identified CNVs associated with BC (731 autosomal CNVs with >5% frequency in the total sample and Q < 0.05). Thirteen CNVs were identified when Partek algorithm-called CNVs were overlapped with Nexus-identified CNVs; these CNVs showed concordances for frequency, effect size, and direction. Coding genes present within BC-associated CNVs were known to play a role in disease etiology and prognosis. Long noncoding RNAs identified within CNVs showed tissue-specific expression, indicating potential functional relevance of the findings. The identified candidate CNVs warrant independent replication.

Cutting Edge: Chronic NF-κB Activation in CD4+ T Cells in Rheumatoid Arthritis Is Genetically Determined by HLA Risk Alleles

Of identified genetic variants, HLA polymorphisms confer the greatest risk for developing autoimmune diseases, including rheumatoid arthritis (HLA-DRB1*04). There are strong influences of HLA polymorphisms on cell type-specific gene expression in B cells and monocytes. Their influence on gene expression in CD4(+) T cells is not known. We determined transcript and proteins levels of target genes in lymphocyte/monocyte subsets in healthy controls and rheumatoid arthritis subjects as a function of HLA-DRB1*04 haplotype. We identified gene expression dependent on HLA-DRB1*04 genotype in CD4(+) T cells. NF-κB activity in CD4(+) T cells was also dependent on HLA-DRB1*04 genotype, and blocking HLA-DR inhibited NF-κB activity in CD4(+) T cells and normalized gene expression, as did pharmacologic inhibition of NF-κB. We conclude that interactions between TCR and MHC class II encoded by HLA-DRB1*04 create a proinflammatory "hum" altering CD4(+) T cell phenotype.

Copy number variation distribution in six monozygotic twin pairs discordant for schizophrenia

We have evaluated copy number variants (CNVs) in six monozygotic twin pairs discordant for schizophrenia. The data from Affymetrix® Human SNP 6.0 arrays™ were analyzed using Affymetrix® Genotyping Console™, Partek® Genomics Suite™, PennCNV, and Golden Helix SVS™. This yielded both program-specific and overlapping results. Only CNVs called by Affymetrix Genotyping Console, Partek Genomics Suite, and PennCNV were used in further analysis. This analysis included an assessment of calls in each of the six twin pairs towards identification of unique CNVs in affected and unaffected co-twins. Real time polymerase chain reaction (PCR) experiments confirmed one CNV loss at 7q11.21 that was found in the affected patient but not in the unaffected twin. The results identified CNVs and genes that were previously implicated in mental abnormalities in four of the six twin pairs. It included PYY (twin pairs 1 and 5), EPHA3 (twin pair 3), KIAA1211L (twin pair 4), and GPR139 (twin pair 5). They represent likely candidate genes and CNVs for the discordance of four of the six monozygotic twin pairs for this heterogeneous neurodevelopmental disorder. An explanation for these differences is ontogenetic de novo events that differentiate in the monozygotic twins during development.

Use of autocorrelation scanning in DNA copy number analysis

MOTIVATION

Data quality is a critical issue in the analyses of DNA copy number alterations obtained from microarrays. It is commonly assumed that copy number alteration data can be modeled as piecewise constant and the measurement errors of different probes are independent. However, these assumptions do not always hold in practice. In some published datasets, we find that measurement errors are highly correlated between probes that interrogate nearby genomic loci, and the piecewise-constant model does not fit the data well. The correlated errors cause problems in downstream analysis, leading to a large number of DNA segments falsely identified as having copy number gains and losses.

METHOD

We developed a simple tool, called autocorrelation scanning profile, to assess the dependence of measurement error between neighboring probes.

RESULTS

Autocorrelation scanning profile can be used to check data quality and refine the analysis of DNA copy number data, which we demonstrate in some typical datasets.

CONTACT

lzhangli@mdanderson.org.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Copy number variation analysis of matched ovarian primary tumors and peritoneal metastasis

Ovarian cancer is the most deadly gynecological cancer. The high rate of mortality is due to the large tumor burden with extensive metastatic lesion of the abdominal cavity. Despite initial chemosensitivity and improved surgical procedures, abdominal recurrence remains an issue and results in patients' poor prognosis. Transcriptomic and genetic studies have revealed significant genome pathologies in the primary tumors and yielded important information regarding carcinogenesis. There are, however, few studies on genetic alterations and their consequences in peritoneal metastatic tumors when compared to their matched ovarian primary tumors. We used high-density SNP arrays to investigate copy number variations in matched primary and metastatic ovarian cancer from 9 patients. Here we show that copy number variations acquired by ovarian tumors are significantly different between matched primary and metastatic tumors and these are likely due to different functional requirements. We show that these copy number variations clearly differentially affect specific pathways including the JAK/STAT and cytokine signaling pathways. While many have shown complex involvement of cytokines in the ovarian cancer environment we provide evidence that ovarian tumors have specific copy number variation differences in many of these genes.