Copy number variants, diseases and gene expression

Copy number variations within fibroblast growth factor 13 gene influence growth traits and alternative splicing in cattle

Previous researches revealed a copy number variation (CNV) region in the bovine fibroblast growth factor 13 (FGF13) gene. However, its effects remain unknown. This study detected the various copy number types in seven Chinese cattle breeds and analysed their population genetic characteristics and effects on growth traits and transcription levels. Copy number Loss was more frequent in Caoyuan Red cattle and Xianan cattle than in the other breeds. Association analysis between CNV and growth traits of Qinchuan indicated that the CNV was significantly related to chest depth, hip width and hucklebone width (P < 0.05). Additionally, the growth traits of individuals with copy number Loss were significantly inferior to those with copy number Gain or Median (P < 0.05). Besides, we found two splicing isoforms, AS1 and AS2, in FGF13 gene, which resulted from alternative 5' splicing sites of intron 1. These isoforms showed varied expression levels in various tissues. Moreover, CNV was significantly and negatively associated with the mRNA expression of AS1 (r = -0.525, P < 0.05). The CNVs in bovine FGF13 gene negatively regulated growth traits and gene transcription. These observations provide new insights into bovine FGF13 gene, delivering potentially useful information for future Chinese cattle breeding programs.

Genome integrity as a potential index of longevity in Ashkenazi Centenarian's families

The aging process, or senescence, is characterized by age-specific decline in physical and physiological function, and increased frailty and genomic changes, including mutation accumulation. However, the mechanisms through which changes in genomic architecture influence human longevity have remained obscure. Copy number variants (CNVs), an abundant class of genomic variants, offer unique opportunities for understanding age-related genomic changes. Here we report the spectrum of CNVs in a cohort of 670 Ashkenazi Jewish centenarians, their progeny, and unrelated controls. The average ages of these groups were 97.4 ± 2.8, 69.2 ± 9.2, and 66.5 ± 7.0 respectively. For the first time, we compared different size classes of CNVs, from 1 kB to 100 MB in size. Using a high-resolution custom Affymetrix array, targeting 44,639 genomic regions, we identified a total of 12,166, 22,188, and 10,285 CNVs in centenarians, their progeny, and control groups, respectively. Interestingly, the offspring group showed the highest number of unique CNVs, followed by control and centenarians. While both gains and losses were found in all three groups, centenarians showed a significantly higher average number of both total gains and losses relative to their controls (p < 0.0327, 0.0182, respectively). Moreover, centenarians showed a lower total length of genomic material lost, suggesting that they may maintain superior genomic integrity over time. We also observe a significance fold increase of CNVs among the offspring, implying greater genomic integrity and a putative mechanism for longevity preservation. Genomic regions that experienced loss or gains appear to be distributed across many sites in the genome and contain genes involved in DNA transcription, cellular transport, developmental pathways, and metabolic functions. Our findings suggest that the exceptional longevity observed in centenarians may be attributed to the prolonged maintenance of functionally important genes. These genes are intrinsic to specific genomic regions as well as to the overall integrity of the genomic architecture. Additionally, a strong association between longer CNVs and differential gene expression observed in this study supports the notion that genomic integrity could positively influence longevity.

Multiomics Analysis of Interleukin-21 as a Potential Immunologic and Biomarker in Hepatocellular Carcinoma

Interleukin-21 (IL-21) is an important antitumor cytokine that contributes to the proliferation and differentiation of CD8 + T cells. It has been proven to enhance the response to immune checkpoint inhibitors (ICIs) in various solid tumors. However, its role in hepatocellular carcinoma (HCC) has not yet been clarified. In this research, we aimed to investigate the antitumor effect of IL-21 in HCC and its effect on ICI treatment. Through transcriptome sequencing analysis and immunohistochemistry validation, we found that patients with high IL-21 expression had a better prognosis. HCCs with high expression of IL-21 had higher infiltration of CD8 + T cells, increased expression of immune checkpoints, and an improved response to ICI treatment. In conclusion, IL-21 can enhance the efficacy of ICI treatment and improve the prognosis of patients by promoting the infiltration of CD8 + T cells and the expression of immune checkpoint-related genes.

Genome-wide association study reveals 2 copy number variations associated with the variation of plumage color in the white duck hybrid population

Plumage color is an intuitive external poultry characteristic with rich manifestations and complex genetic mechanisms. In our previous study, we observed that there were more dark variations in plumage color in the F2 population derived from the hybridization of 2 white duck varieties. Therefore, based on the statistics of plumage color of 308 F2 populations, we further used the resequencing data of these individuals to detect copy number variations (CNVs) in the whole genome and conducted genome-wide association studies (GWAS) to determine the genetic basis related to plumage color traits. The CNV detection revealed 9,337 CNVs, with an average length of 15,950 bp and a total length of 142.02 MB, accounting for approximately 12.91% of the reference genome. The CNV distribution on the chromosomes was relatively uniform, and the number of CNVs on each chromosome positively correlated with the length of the chromosome. In the pure black plumage group, 2,101 CNVs were only identified, and 1,714 were specifically identified in the pure white plumage group. Ten CNVs were randomly selected for validation using quantitative real-time PCR, and 9 CNVs had the same CNV types as predicted, with an accuracy of 90%. Based on GWAS, we identified 2 CNVs potentially associated with plumage color variations, with the associated CNV regions covering 9 genes. Enrichment analysis of these 9 candidate genes showed significant enrichment of 3 pathways (ribosome biogenesis in eukaryotes, RNA transport, and protein export) and 17 gene ontology terms. Among these, VWA5A can downregulate MITF by binding to the regulatory factors SOX10. The occurrence of CNV may indirectly contribute to duck plumage color variation by affecting the regulatory factors of the switch gene MITF in the melanogenesis pathway. These findings have improved the understanding of the genetic basis of duck plumage color variation and have been beneficial for developing and using plumage color traits in subsequent poultry breeding.

Identification of a rare copy number polymorphic gain at 3q12.2 with candidate genes for familial endometriosis

Endometriosis is a complex disease that affects 10-15% of women of reproductive age. Familial studies show that relatives of affected patients have a higher risk of developing the disease, implicating a genetic role for this disorder. Little is known about the impact of germline genomic copy number variant (CNV) polymorphisms on the heredity of the disease. In this study, we describe a rare CNV identified in two sisters with familial endometriosis, which contain genes that may increase the susceptibility and progression of this disease. We investigated the presence of CNVs from the endometrium and blood of the sisters with endometriosis and normal endometrium of five women as controls without the disease using array-CGH through the Agilent 2x400K platform. We excluded common CNVs that were present in the database of genomic variation. We identified, in both sisters, a rare CNV gain affecting 113kb at band 3q12.2 involving two candidate genes: ADGRG7 and TFG. The CNV gain was validated by qPCR. ADGRG7 is located at 3q12.2 and encodes a G protein-coupled receptor influencing the NF-kappaβ pathway. TFG participates in chromosomal translocations associated with hematologic tumor and soft tissue sarcomas, and is also involved in the NF-kappa B pathway. The CNV gain in this family provides a new candidate genetic marker for future familial endometriosis studies. Additional longitudinal studies of affected families must confirm any associations between this rare CNV gain and genes involved in the NF-kappaβ pathway in predisposition to endometriosis.

Distinct immune escape and microenvironment between RG-like and pri-OPC-like glioma revealed by single-cell RNA-seq analysis

The association of neurogenesis and gliogenesis with glioma remains unclear. By conducting single-cell RNA-seq analyses on 26 gliomas, we reported their classification into primitive oligodendrocyte precursor cell (pri-OPC)-like and radial glia (RG)-like tumors and validated it in a public cohort and TCGA glioma. The RG-like tumors exhibited wild-type isocitrate dehydrogenase and tended to carry EGFR mutations, and the pri-OPC-like ones were prone to carrying TP53 mutations. Tumor subclones only in pri-OPC-like tumors showed substantially down-regulated MHC-I genes, suggesting their distinct immune evasion programs. Furthermore, the two subgroups appeared to extensively modulate glioma-infiltrating lymphocytes in distinct manners. Some specific genes not expressed in normal immune cells were found in glioma-infiltrating lymphocytes. For example, glial/glioma stem cell markers OLIG1/PTPRZ1 and B cell-specific receptors IGLC2/IGKC were expressed in pri-OPC-like and RG-like glioma-infiltrating lymphocytes, respectively. Their expression was positively correlated with those of immune checkpoint genes (e.g., LGALS33) and poor survivals as validated by the increased expression of LGALS3 upon IGKC overexpression in Jurkat cells. This finding indicated a potential inhibitory role in tumor-infiltrating lymphocytes and could provide a new way of cancer immune evasion.

Exploring quantitative traits-associated copy number deletions through reanalysis of UK10K consortium whole genome sequencing cohorts

OBJECTIVES

We performed comprehensive association analyses of common high-confidence gnomAD-reported copy number deletions (CNDs) with 60 quantitative traits from UK10K consortium WGS data.

METHODS

The study made use of data generated by the UK10K Consortium. UK10K consortium WGS data consist of TwinsUK (n = 1754, middle-aged females) and ALSPAC (n = 1867, birth to adolescence) cohorts. UK10K consortium called 18,739 CNDs (hg19) with GenomeSTRiP software. After filtering out variants with minor allele frequency < 0.05 or HWE P < 1.0 × 10- 6, 1222 (TwinsUK) and 1211 (ALSPAC) CNDs remained for association analyses with 60 normalized quantitative traits.

RESULTS

We identified 23 genome-wide significant associations at 13 loci, among which 2 associations reached experiment-wide significance. We found that two common deletions in chromosome 4, located between WDR1 and ZNF518B (23.3 kb, dbVar ID:nssv15888957, 4:10211262-10,234,569 and 9.8 kb, dbVar ID:nssv15888975, 4:10392422-10,402,191), were associated with uric acid levels (P = 5.23 × 10- 11 and 2.29 × 10- 8, respectively). We also discovered a novel deletion spanning chromosome 18 (823 bp, dbVar ID: nssv15841628, 8:74347187-74,348,010) associated with low HDL cholesterol levels (P = 4.15 × 10- 7). Additionally, we observed two red blood cell traits-associated loci with genome-wide significance, a 13.2 kb deletion in 7q22.1 (nssv15922542) and a 3.7 kb deletion in 12q24.12 (nssv15813226), both of which were located in regions previously reported to be associated with red blood cell traits. Two deletions in 11q11 (nssv15803200 and nssv15802240), where clusters of multiple olfactory receptor genes exist, and a deletion (nssv15929560) upstream to DOCK5 were associated with childhood obesity. Finally, when defining Trait-Associated copy number Deletions (TADs) as CNDs with phenotype associations at sub-threshold significance (P < 10- 3), we identified 157 (97.5%) out of 161 TADs in non-coding regions, with a mean size of 4 kb (range: 209 - 47,942 bp).

CONCLUSION

We conducted a reanalysis of the UK10K Whole Genome Sequencing cohort, which led to the identification of multiple high confidence copy number deletions associated with quantitative traits. These deletions have standard dbVar IDs and replicate previous findings, as well as reveal novel loci that require further replication studies.

Palmitoylation landscapes across human cancers reveal a role of palmitoylation in tumorigenesis

BACKGROUND

Protein palmitoylation, which is catalyzed by palmitoyl-transferase and de-palmitoyl-transferase, plays a crucial role in various biological processes. However, the landscape and dynamics of protein palmitoylation in human cancers are not well understood.

METHODS

We utilized 23 palmitoyl-acyltransferases and seven de-palmitoyl-acyltransferases as palmitoylation-related genes for protein palmitoylation analysis. Multiple publicly available datasets were employed to conduct pan-cancer analysis, examining the transcriptome, genomic alterations, clinical outcomes, and correlation with c-Myc (Myc) for palmitoylation-related genes. Real-time quantitative PCR and immunoblotting were performed to assess the expression of palmitoylation-related genes and global protein palmitoylation levels in cancer cells treated with Myc depletion or small molecule inhibitors. Protein docking and drug sensitivity analyses were employed to predict small molecules that target palmitoylation-related genes.

RESULTS

We identified associations between palmitoylation and cancer subtype, stage, and patient survival. We discovered that abnormal DNA methylation and oncogenic Myc-driven transcriptional regulation synergistically contribute to the dysregulation of palmitoylation-related genes. This dysregulation of palmitoylation was closely correlated with immune infiltration in the tumor microenvironment and the response to immunotherapy. Importantly, dysregulated palmitoylation was found to modulate canonical cancer-related pathways, thus influencing tumorigenesis. To support our findings, we performed a proof-of-concept experiment showing that depletion of Myc led to reduced expression of most palmitoylation-related genes, resulting in decreased global protein palmitoylation levels. Through mass spectrometry and enrichment analyses, we also identified palmitoyl-acyltransferases ZDHHC7 and ZDHHC23 as significant contributors to mTOR signaling, DNA repair, and immune pathways, highlighting their potential roles in tumorigenesis. Additionally, our study explored the potential of three small molecular (BI-2531, etoposide, and piperlongumine) to modulate palmitoylation by targeting the expression or activity of palmitoylation-related genes or enzymes.

CONCLUSIONS

Overall, our findings underscore the critical role of dysregulated palmitoylation in tumorigenesis and the response to immunotherapy, mediated through classical cancer-related pathways and immune cell infiltration. Additionally, we propose that the aforementioned three small molecule hold promise as potential therapeutics for modulating palmitoylation, thereby offering novel avenues for cancer therapy.

Copy Number Variation Regions Differing in Segregation Patterns Span Different Sets of Genes

Copy number variations regions (CNVRs) can be classified either as segregating, when found in both parents, and offspring, or non-segregating. A total of 65 segregating and 31 non-segregating CNVRs identified in at least 10 individuals within a dense pedigree of the Gochu Asturcelta pig breed was subjected to enrichment and functional annotation analyses to ascertain their functional independence and importance. Enrichment analyses allowed us to annotate 1018 and 351 candidate genes within the bounds of the segregating and non-segregating CNVRs, respectively. The information retrieved suggested that the candidate genes spanned by segregating and non-segregating CNVRs were functionally independent. Functional annotation analyses allowed us to identify nine different significantly enriched functional annotation clusters (ACs) in segregating CNVR candidate genes mainly involved in immunity and regulation of the cell cycle. Up to five significantly enriched ACs, mainly involved in reproduction and meat quality, were identified in non-segregating CNVRs. The current analysis fits with previous reports suggesting that segregating CNVRs would explain performance at the population level, whereas non-segregating CNVRs could explain between-individuals differences in performance.

A comparison of algorithms for identifying copy number variants in family-based whole-exome sequencing data and its implications in inheritance pattern analysis

There remain challenges in accurately identifying constitutional or germline copy number variants (gCNVs) based on whole-exome sequencing data that have implications for genetic diagnosis for 'rare undiagnosed disease' in the clinical setting. Although multiple algorithms have been proposed, a systematic comparison of these algorithms for calling gCNVs and analyzing inherited pattern have yet to be fully conducted. Therefore, we empirically compared seven exome-based algorithms, including XHMM, CLAMMS, CODEX2, ExomeDepth, DECoN, CN.MOPS, and GATK gCNV, for calling gCNVs in 151 individuals from 44 pedigrees, together with the gold standard of genotyping-derived gCNVs in the same cohort for the performance assessment. These algorithms demonstrated varied powers in identifying gCNVs, although the distribution of gCNVs size was similar. The number of shared gCNVs across these algorithms was limited (e.g., only four gCNVs shared among seven algorithms); however, several algorithms showed varying degrees of consistency (e.g., 1,843 gCNVs shared between DECoN and ExomeDepth). CLAMMS and CODEX2 outperformed the remaining algorithms according to a relatively higher F-score (i.e., 0.145 and 0.152, respectively). In addition, these algorithms exhibited different Mendelian inconsistencies of gCNVs and significant challenges remained in inheritance pattern analysis. In conclusion, selecting good algorithms may have important implications in gCNVs-based inheritance pattern analysis for family-based studies.

Whole genome sequencing of simmental cattle for SNP and CNV discovery

BACKGROUD

The single nucleotide polymorphisms (SNPs) and copy number variations (CNVs) are two major genomic variants, which play crucial roles in evolutionary and phenotypic diversity.

RESULTS

In this study, we performed a comprehensive analysis to explore the genetic variations (SNPs and CNVs) of high sperm motility (HSM) and poor sperm motility (PSM) Simmental bulls using the high-coverage (25×) short-read next generation sequencing and single-molecule long reads sequencing data. A total of ~ 15 million SNPs and 2,944 CNV regions (CNVRs) were detected in Simmental bulls, and a set of positive selected genes (PSGs) and CNVRs were found to be overlapped with quantitative trait loci (QTLs) involving immunity, muscle development, reproduction, etc. In addition, we detected two new variants in LEPR, which may be related to the artificial breeding to improve important economic traits. Moreover, a set of genes and pathways functionally related to male fertility were identified. Remarkably, a CNV on SPAG16 (chr2:101,427,468 - 101,429,883) was completely deleted in all poor sperm motility (PSM) bulls and half of the bulls in high sperm motility (HSM), which may play a crucial role in the bull-fertility.

CONCLUSIONS

In conclusion, this study provides a valuable genetic variation resource for the cattle breeding and selection programs.

Large-scale integration of the non-coding RNAs with DNA methylation in human cancers

Characterizing influences of DNA methylation (DNAm) on non-coding RNAs (ncRNAs) is important to understand the mechanisms of gene regulation and cancer outcome. In our study, we describe the results of ncRNA quantitative trait methylation sites (ncQTM) analyses on 8,545 samples from The Cancer Genome Atlas (TCGA), 763 samples from the Clinical Proteomic Tumor Analysis Consortium (CPTAC), and 516 samples from Genotype-Tissue Expression (GTEx) to identify the significant associations between DNAm sites and ncRNAs (miRNA, long non-coding RNA [lncRNA], small nuclear RNA [snRNA], small nucleolar RNA [snoRNA], and rRNA) across 32 cancer types. With more than 22 billion tests, we identify 302,764 cis-ncQTMs (6.28% of all tested) and 79,841,728 trans-ncQTMs (1.15% of all tested). Most DNAm sites (70.6% on average) are in trans association, while only 25.2% DNAm sites are in cis association. Further, we develop a subtype named ncm_cluster based on cancer-specific ncRNAs thatis associated with tumor microenvironment, clinical outcome, and biological pathways. To comprehensively describe the ncQTM patterns, we developed a database named Pancan-ncQTM (http://bigdata.njmu.edu.cn/Pancan-ncQTM/).

Saxitoxins-producing potential of the marine dinoflagellate Alexandrium affine and its environmental implications revealed by toxins and transcriptome profiling

The marine dinoflagellate Alexandrium occurs widely in coastal waters, and some of them can produce saxitoxins (STXs) that cause paralytic shellfish poisoning (PSP). Alexandrium affine is a harmful algal bloom (HAB)-forming species off the coast of Asia; however, its ability to produce STXs has been controversial. In the present study, we detected STXs in A. affine Alex02 isolated from the southern coast of Korea. The total STXs equivalent (STXs eq) and profiles of Alex02 varied depending on the tested environmental conditions, including the temperature and nitrate concentrations. STXs toxicity levels of A. affine Alex02 (<0.8 STXs eq fmol cell^-1) were significantly lower than those of toxic A. catenella Alex03 and A. pacificum Alex05. On a genetic basis, we identified all the STX biosynthesis sxt genes, except sxtX in A. affine, via large-scale transcriptome analysis. Interestingly, the two proteins, sxtA4 and sxtG, were similar in sequence and domain structure to those of other toxic dinoflagellates and cyanobacteria; however, their transcript levels were extremely low. Our results suggest that A. affine has the potential to produce STXs, while its toxicity is much lower or negligible, which is unlikely to cause PSP incidents in marine environments.

Upregulation of TUBG1 expression promotes hepatocellular carcinoma development

Tubulin γ-1 (TUBG1) is a highly conserved component of the centrosome and its deregulation is involved in the development of several types of cancer. However, the role of TUBG1 in hepatocellular carcinoma (HCC) remains unclear. In this study, we found that TUBG1 was upregulated in human HCC cells and tissues and that TUBG1 upregulation was associated with promoter hypomethylation in HCC tissues. TUBG1 knockdown suppressed the proliferation, invasion, and migration of HCC cells. While TUBG1 expression was positively correlated with CD4 + memory T lymphocyte infiltration, it was negatively correlated with CD4 + regulatory T-cell infiltration in human HCC tissues. Furthermore, TUBG1 expression was positively correlated with the expression of genes involved in cell division. Noticeably, high expression of TUBG1 was associated with poor prognosis in patients with HCC. Overall, our findings revealed that TUBG1 promotes hepatocarcinogenesis by increasing proliferation, invasion, and migration of HCC cells and may regulate T lymphocyte infiltration. The current findings provide important insights into TUBG1 regulation in HCC, which could provide new therapeutic targets for hepatocarcinoma which has a very high incidence and mortality rate worldwide.

Assessment of copy number variants in three Brazilian locally adapted cattle breeds using whole-genome re-sequencing data

Further characterization of genetic structural variations should strongly focus on small and endangered local breeds given their role in unraveling genes and structural variants underlying selective pressures and phenotype variation. A comprehensive genome-wide assessment of copy number variations (CNVs) based on whole-genome re-sequencing data was performed on three Brazilian locally adapted cattle breeds (Caracu Caldeano, Crioulo Lageano, and Pantaneiro) using the ARS-UCD1.2 genome assembly. Data from 36 individuals with an average coverage depth of 14.07× per individual was used. A total of 24 945 CNVs were identified distributed among the breeds (Caracu Caldeano = 7285, Crioulo Lageano = 7297, and Pantaneiro = 10 363). Deletion events were 1.75-2.07-fold higher than duplications, and the total length of CNVs is composed mostly of a high number of segments between 10 and 30 kb. CNV regions (CNVRs) are not uniformly scattered throughout the genomes (n = 463), and 105 CNVRs were found overlapping among the studied breeds. Functional annotation of the CNVRs revealed variants with high consequence on protein sequence harboring relevant genes, in which we highlighted the BOLA-DQB, BOLA-DQA5, CD1A, β-defensins, PRG3, and ULBP21 genes. Enrichment analysis based on the gene list retrieved from the CNVRs disclosed over-represented terms (p < 0.01) strongly associated with immunity and cattle resilience to harsh environments. Additionally, QTL associated with body conformation and dairy-related traits were also unveiled within the CNVRs. These results provide better understanding of the selective forces shaping the genome of such cattle breeds and identify traces of natural selection pressures by which these populations have been exposed to challenging environmental conditions.

rAAV-CRISPRa therapy corrects Rai1 haploinsufficiency and rescues selective disease features in Smith-Magenis syndrome mice

Haploinsufficiency in retinoic acid induced 1 (RAI1) causes Smith-Magenis syndrome (SMS), a severe neurodevelopmental disorder characterized by neurocognitive deficits and obesity. Currently, curative treatments for SMS do not exist. Here, we take a recombinant adeno-associated virus (rAAV)-clustered regularly interspaced short palindromic repeats activation (CRISPRa) approach to increase expression of the remaining intact Rai1 allele. Building upon our previous work that found the paraventricular nucleus of hypothalamus plays a central role in SMS pathogenesis, we performed paraventricular nucleus of hypothalamus-specific rAAV-CRISPRa therapy by increasing endogenous Rai1 expression in SMS (Rai1^±) mice. We found that rAAV-CRISPRa therapy rescues excessive repetitive behavior, delays the onset of obesity, and partially reduces hyperphagia in SMS mice. Our work provides evidence that rAAV-CRISPRa therapy during early adolescence can boost the expression of healthy Rai1 allele and modify disease progression in a mouse model of Smith-Magenis syndrome.

Identification of N7-methylguanosine related subtypes and construction of prognostic model in gastric cancer

Background

N⁷-methylguanosine (m7G), one of the most common post-transcriptional modifications, can be present in tRNA, mRNA, and miRNA to mediate the progression of various tumors. However, the possible role of m7G in gastric cancer (GC) is still unknown.

Materials and Methods

In this study, SNVs (single nucleotide variations), CNVs (copy number variations), and methylation of m7G-related genes (m7GRGs) were analyzed. The relationship between them and the expression of m7GRGs and prognosis of GC patients was explored. Based on 13 prognostic-related m7GRGs, 567 GC samples were classified into three subtypes using the ConsensusClusterPlus package. we compared survival status, clinical traits, immune cell infiltration, immune checkpoints, tumor microenvironment (TME), tumor immune dysfunction and exclusion (TIDE), and potential biological pathways among the three subtypes. Then, patients were again grouped into different genetic subtypes based on the DEGs among the three subtypes. In addition, a prognostic m7GRG_Score was constructed using five risk genes applicable to patients of any age, gender and stage. We also assessed tumor mutational burden (TMB), microsatellite instability (MSI), cancer stem cell (CSC) index, sensitivity of antineoplastic drugs, efficacy of anti-PD-1 and anti-CTLA4 immunotherapy between high and low m7GRG_Score groups. Finally, we established a nomogram based on m7GRG_Score and tumor stage to enhance the clinical application of the model. miRNAs and lncRNAs that could regulate expression of risk genes were searched.

Results

SNVs, CNVs, and methylation of m7GRGs were associated with m7GRGs expression. However, they did not significantly affect the survival of GC patients. Our results also confirmed that patients in subtypes B and C and low m7GRG_Score groups had longer survival time, better clinical stage, more immune cell infiltration, fewer immune escape and dysfunction compared to subtype A and high m7GRG_Score groups. A low m7GRG_score was featured with increased microsatellite instability-high (MSI-H), TMB, and efficacy of immunotherapy.

Conclusion

The m7GRG_Score model may become a beneficial tool for predicting prognosis and guiding personalized treatment in GC patients. These findings will improve our knowledge of m7G in GC and provide new methods for more effective treatment strategies.

Complement component C4 structural variation and quantitative traits contribute to sex-biased vulnerability in systemic sclerosis

Copy number (CN) polymorphisms of complement C4 play distinct roles in many conditions, including immune-mediated diseases. We investigated the association of C4 CN with systemic sclerosis (SSc) risk. Imputed total C4, C4A, C4B, and HERV-K CN were analyzed in 26,633 individuals and validated in an independent cohort. Our results showed that higher C4 CN confers protection to SSc, and deviations from CN parity of C4A and C4B augmented risk. The protection contributed per copy of C4A and C4B differed by sex. Stronger protection was afforded by C4A in men and by C4B in women. C4 CN correlated well with its gene expression and serum protein levels, and less C4 was detected for both in SSc patients. Conditioned analysis suggests that C4 genetics strongly contributes to the SSc association within the major histocompatibility complex locus and highlights classical alleles and amino acid variants of HLA-DRB1 and HLA-DPB1 as C4-independent signals.

CiberAMP: An R Package to Identify Differential mRNA Expression Linked to Somatic Copy Number Variations in Cancer Datasets

Simple Summary

The ability to establish accurate correlations between the number of copies of genes and the expression levels of their encoded transcripts remains a challenge despite the extensive progress made in the understanding of the genome of cancer cells. Here, we describe a new algorithm that does so by integrating both genomics and transcriptomics data from the Cancer Genome Atlas. In addition to explaining the step-by-step basis of this new method, we provide examples of how this new algorithm can help identify functionally meaningful gene copy alterations that are recurrently detected in cancer patients.

Abstract

Somatic copy number variations (SCNVs) are genetic alterations frequently found in cancer cells. These genetic alterations can lead to concomitant perturbations in the expression of the genes included in them and, as a result, promote a selective advantage to cancer cells. However, this is not always the case. Due to this, it is important to develop in silico tools to facilitate the accurate identification and functional cataloging of gene expression changes associated with SCNVs from pan-cancer data. Here, we present a new R-coded tool, designated as CiberAMP, which utilizes genomic and transcriptomic data contained in the Cancer Genome Atlas (TCGA) to identify such events. It also includes information on the genomic context in which such SCNVs take place. By doing so, CiberAMP provides clues about the potential functional relevance of each of the SCNV-associated gene expression changes found in the interrogated tumor samples. The main features and advantages of this new algorithm are illustrated using glioblastoma data from the TCGA database.

Genome-wide detection of copy number variation in American mink using whole-genome sequencing

BACKGROUND

Copy number variations (CNVs) represent a major source of genetic diversity and contribute to the phenotypic variation of economically important traits in livestock species. In this study, we report the first genome-wide CNV analysis of American mink using whole-genome sequence data from 100 individuals. The analyses were performed by three complementary software programs including CNVpytor, DELLY and Manta.

RESULTS

A total of 164,733 CNVs (144,517 deletions and 20,216 duplications) were identified representing 5378 CNV regions (CNVR) after merging overlapping CNVs, covering 47.3 Mb (1.9%) of the mink autosomal genome. Gene Ontology and KEGG pathway enrichment analyses of 1391 genes that overlapped CNVR revealed potential role of CNVs in a wide range of biological, molecular and cellular functions, e.g., pathways related to growth (regulation of actin cytoskeleton, and cAMP signaling pathways), behavior (axon guidance, circadian entrainment, and glutamatergic synapse), lipid metabolism (phospholipid binding, sphingolipid metabolism and regulation of lipolysis in adipocytes), and immune response (Wnt signaling, Fc receptor signaling, and GTPase regulator activity pathways). Furthermore, several CNVR-harbored genes associated with fur characteristics and development (MYO5A, RAB27B, FGF12, SLC7A11, EXOC2), and immune system processes (SWAP70, FYN, ORAI1, TRPM2, and FOXO3).

CONCLUSIONS

This study presents the first genome-wide CNV map of American mink. We identified 5378 CNVR in the mink genome and investigated genes that overlapped with CNVR. The results suggest potential links with mink behaviour as well as their possible impact on fur quality and immune response. Overall, the results provide new resources for mink genome analysis, serving as a guideline for future investigations in which genomic structural variations are present.

Prognostic significance of copy number variation in B-cell acute lymphoblastic leukemia

Copy number variations (CNVs) are widespread in both pediatric and adult cases of B-cell acute lymphoblastic leukemia (B-ALL); however, their clinical significance remains unclear. This review primarily discusses the most prevalent CNVs in B-ALL to elucidate their clinical value and further personalized management of this population. The discovery of the molecular mechanism of gene deletion and the development of targeted drugs will further enhance the clinical prognosis of B-ALL.

Multiple Recurrent Copy Number Variations (CNVs) in Chromosome 22 Including 22q11.2 Associated with Autism Spectrum Disorder

Introduction

Autism spectrum disorder (ASD) is a developmental disorder that can cause substantial social, communication, and behavioral challenges. Genetic factors play a significant role in ASD, where the risk of ASD has been increased for unclear reasons. Twin studies have shown important evidence of both genetic and environmental contributions in ASD, where the level of contribution of these factors has not been proven yet. It has been suggested that copy number variation (CNV) duplication and the deletion of many genes in chromosome 22 (Ch22) may have a strong association with ASD. This study screened the CNVs in Ch22 in autistic Saudi children and assessed the candidate gene in the CNVs region of Ch22 that is most associated with ASD.

Methods

This study included 15 autistic Saudi children as well as 4 healthy children as controls; DNA was extracted from samples and analyzed using array comparative genomic hybridization (aCGH) and DNA sequencing.

Results

The aCGH detected (in only 6 autistic samples) deletion and duplication in many regions of Ch22, including some critical genes. Moreover, DNA sequencing determined a genetic mutation in the TBX1 gene sequence in autistic samples. This study, carried out using aCGH, found that six autistic patients had CNVs in Ch22, and DNA sequencing revealed mutations in the TBX1 gene in autistic samples but none in the control.

Conclusion

CNV deletion and the duplication of the TBX1 gene could be related to ASD; therefore, this gene needs more analysis in terms of expression levels.

Skewed X-Chromosome Inactivation as a Possible Marker of X-Linked CNV in Women with Pregnancy Loss

Skewed X-chromosome inactivation (sXCI) can be a marker of lethal genetic variants on the X chromosome in a woman since sXCI modifies the pathological phenotype. The aim of this study was to search for CNVs in women with miscarriages and sXCI. XCI was assayed using the classical method based on the amplification of highly polymorphic exon 1 of the androgen receptor (AR) gene. The XCI status was analysed in 313 women with pregnancy loss and in 87 spontaneously aborted embryos with 46,XX karyotype, as well as in control groups of 135 women without pregnancy loss and 64 embryos with 46,XX karyotype from induced abortions in women who terminated a normal pregnancy. The frequency of sXCI differed significantly between women with miscarriages and women without pregnancy losses (6.3% and 2.2%, respectively; p = 0.019). To exclude primary causes of sXCI, sequencing of the XIST and XACT genes was performed. The XIST and XACT gene sequencing revealed no known pathogenic variants that could lead to sXCI. Molecular karyotyping was performed using aCGH, followed by verification of X-linked CNVs by RT-PCR and MLPA. Microdeletions at Xp11.23 and Xq24 as well as gains of Xq28 were detected in women with sXCI and pregnancy loss.

Sequence polymorphisms of PR39 cathelicidins and extensive copy variations in commercial pig breeds

Copy number polymorphisms (CNPs) of antimicrobial peptides (AMPs) in livestock can influence the innate immune response of individuals. We conducted a high-resolution analysis of the genomic variations of porcine cathelicidin PR39 using cloned PR39 amplicons corresponding to the 5' untranslated region (UTR) to 3' UTR from four individuals of three different pig breeds. We identified 15 different sequences corresponding to 9 different coding domain sequences (CDSs), encoding 7 different protein sequences consisting of 3 functional and 4 non-functional forms. Subsequently, we developed a PR39 CNP typing method using real-time polymerase chain reaction (PCR) and analyzed the PR39 copy numbers from 44 pigs of six breeds. Significant variations in PR39 copies ranging from 2 to 10 copies, with a mean copy number of 5, were observed among all commercial breeds, except the wild boar. Among the different breeds, the PR39 copy number was highest (10) in Korean native pigs. Gene expression analysis showed that PR39 expression was correlated with the copy number. Moreover, the comparative analysis of the cathelicidin cluster-containing region among eight mammalian species showed the complete evolutionary conservation of the region, except for differences in the degree of cathelicidin expansion in each species. Therefore, characterization of CNPs in AMP genes could aid in improving the genetic potential of innate immune responses in livestock animals.

The Role of Copy Number Variants in Gene Co-Expression Patterns for Luminal B Breast Tumors

Gene co-expression networks have become a usual approach to integrate the vast amounts of information coming from gene expression studies in cancer cohorts. The reprogramming of the gene regulatory control and the molecular pathways depending on such control are central to the characterization of the disease, aiming to unveil the consequences for cancer prognosis and therapeutics. There is, however, a multitude of factors which have been associated with anomalous control of gene expression in cancer. In the particular case of co-expression patterns, we have previously documented a phenomenon of loss of long distance co-expression in several cancer types, including breast cancer. Of the many potential factors that may contribute to this phenomenology, copy number variants (CNVs) have been often discussed. However, no systematic assessment of the role that CNVs may play in shaping gene co-expression patterns in breast cancer has been performed to date. For this reason we have decided to develop such analysis. In this study, we focus on using probabilistic modeling techniques to evaluate to what extent CNVs affect the phenomenon of long/short range co-expression in Luminal B breast tumors. We analyzed the co-expression patterns in chromosome 8, since it is known to be affected by amplifications/deletions during cancer development. We found that the CNVs pattern in chromosome 8 of Luminal B network does not alter the co-expression patterns significantly, which means that the co-expression program in this cancer phenotype is not determined by CNV structure. Additionally, we found that region 8q24.3 is highly dense in interactions, as well as region p21.3. The most connected genes in this network belong to those cytobands and are associated with several manifestations of cancer in different tissues. Interestingly, among the most connected genes, we found MAF1 and POLR3D, which may constitute an axis of regulation of gene transcription, in particular for non-coding RNA species. We believe that by advancing on our knowledge of the molecular mechanisms behind gene regulation in cancer, we will be better equipped, not only to understand tumor biology, but also to broaden the scope of diagnostic, prognostic and therapeutic interventions to ultimately benefit oncologic patients.

Whole-genome analysis of human embryonic stem cells enables rational line selection based on genetic variation

Despite their widespread use in research, there has not yet been a systematic genomic analysis of human embryonic stem cell (hESC) lines at a single-nucleotide resolution. We therefore performed whole-genome sequencing (WGS) of 143 hESC lines and annotated their single-nucleotide and structural genetic variants. We found that while a substantial fraction of hESC lines contained large deleterious structural variants, finer-scale structural and single-nucleotide variants (SNVs) that are ascertainable only through WGS analyses were present in hESC genomes and human blood-derived genomes at similar frequencies. Moreover, WGS allowed us to identify SNVs associated with cancer and other diseases that could alter cellular phenotypes and compromise the safety of hESC-derived cellular products transplanted into humans. As a resource to enable reproducible hESC research and safer translation, we provide a user-friendly WGS data portal and a data-driven scheme for cell line maintenance and selection.

The Psychoemotional Stress-Induced Changes in the Abundance of SatIII (1q12) and Telomere Repeats, but Not Ribosomal DNA, in Human Leukocytes

INTRODUCTION. As shown earlier, copy number variations (CNV) in the human satellite III (1q12) fragment (f-SatIII) and the telomere repeat (TR) reflects the cell’s response to oxidative stress. The contents of f-SatIII and TR in schizophrenic (SZ) patients were found to be lower than in healthy controls (HC) in previous studies. The major question of this study was: ‘What are the f-SatIII and TR CNV dynamic changes in human leukocytes, depending on psychoemotional stress?’ MATERIALS AND METHODS. We chose a model of psychoemotional stress experienced by second-year medical students during their exams. Blood samples were taken in stressful conditions (exams) and in a control non-stressful period. Biotinylated probes were used for f-SatIII, rDNA, and TR quantitation in leukocyte DNA by non-radioactive quantitative hybridization in SZ patients (n = 97), HC (n = 97), and medical students (n = 17, n = 42). A flow cytometry analysis was used for the oxidative stress marker (NOX4, 8-oxodG, and γH2AX) detection in the lymphocytes of the three groups. RESULTS. Oxidative stress markers increased significantly in the students’ lymphocytes during psychoemotional stress. The TR and f-SatIII, but not the rDNA, contents significantly changed in the DNA isolated from human blood leukocytes. After a restoration period (post-examinational vacations), the f-SatIII content decreased, and the TR content increased. Changes in the blood cells of students during examinational stress were similar to those in SZ patients during an exacerbation of the disease. CONCLUSIONS. Psychoemotional stress in students during exams triggers a universal mechanism of oxidative stress. The oxidative stress causes significant changes in the f-SatIII and TR contents, while the ribosomal repeat content remains stable. A hypothesis is proposed to explain the quantitative polymorphisms of f-SatIII and TR contents under transient (e.g., students’ exams) or chronic (in SZ patients) stress. The changes in the f-SatIII and TR copy numbers are non-specific events, irrespective of the source of stress. Thus, our findings suggest that the psychoemotional stress, common in SZ patients and healthy students during exams, but not in a schizophrenia-specific event, was responsible for the changes in the repeat contents that we observed earlier in SZ patients.

A CTCF-Binding Element and Histone Deacetylation Cooperatively Maintain Chromatin Loops, Linking to Long-Range Gene Regulation in Cancer Genomes

Background

Genes spanning long chromosomal domains are coordinately regulated in human genome, which contribute to global gene dysregulation and carcinogenesis in cancer. It has been noticed that epigenetic modification and chromatin architecture may participate in the regulation process. However, the regulation patterns and functional elements of long-range gene regulation are unclear.

Methods

Based on the clinical transcriptome data from different tumor sets, a novel expressional correlation analysis pipeline was performed to classify the co-regulated regions and subsets of intercorrelated regions. The GLAM2 program was used to predict conserved DNA elements that enriched in regions. Two conserved elements were selected to delete in Ishikawa and HeLa cells by CRISPR-Cas9. SAHA treatment and HDAC knockdown were used to change the histone acetylation status. Using qPCR, MTT, and scratch healing assay, we evaluate the effect on gene expression and cancer cell phenotype. By DNA pull-down and ChIP, the element-binding proteins were testified. 3C and 3D-FISH were performed to depict the alteration in chromatin architecture.

Results

In multiple cancer genomes, we classified subsets of coordinately regulated regions (sub-CRRs) that possibly shared the same regulatory mechanisms and exhibited similar expression patterns. A new conserved DNA element (CRE30) was enriched in sub-CRRs and associated with cancer patient survival. CRE30 could restrict gene regulation in sub-CRRs and affect cancer cell phenotypes. DNA pull-down showed that multiple proteins including CTCF were recruited on the CRE30 locus, and ChIP assay confirmed the CTCF-binding signals. Subsequent results uncovered that as an essential element, CRE30 maintained chromatin loops and mediated a compact chromatin architecture. Moreover, we found that blocking global histone deacetylation induced chromatin loop disruption and CTCF dropping in the region containing CRE30, linked to promoted gene regulation. Additionally, similar effects were observed with CRE30 deletion in another locus of chromosome 8.

Conclusions

Our research clarified a new functional element that recruits CTCF and collaborates with histone deacetylation to maintain high-order chromatin organizations, linking to long-range gene regulation in cancer genomes. The findings highlight a close relationship among conserved DNA element, epigenetic modification, and chromatin architecture in long-range gene regulation process.

Contribution of schizophrenia polygenic burden to longitudinal phenotypic variance in 22q11.2 deletion syndrome

While the recurrent 22q11.2 deletion is one of the strongest genetic risk factors for schizophrenia (SCZ), variability of its associated neuropsychiatric endophenotypes reflects its incomplete penetrance for psychosis development. To assess whether this phenotypic variability is linked to common variants associated with SCZ, we studied the association between SCZ polygenic risk score (PRS) and longitudinally acquired phenotypic information of the Swiss 22q11.2DS cohort (n = 97, 50% females, mean age 17.7 yr, mean visit interval 3.8 yr). The SCZ PRS with the best predictive performance was ascertained in the Estonian Biobank (n = 201,146) with LDpred. The infinitesimal SCZ PRS model showed the strongest capacity in discriminating SCZ cases from controls with one SD difference in SCZ PRS corresponding to an odds ratio (OR) of 1.73 (95% CI 1.57-1.90, P = 1.47 × 10^-29). In 22q11.2 patients, random-effects ordinal regression modelling using longitudinal data showed SCZ PRS to have the strongest effect on social anhedonia (OR = 2.09, P = 0.0002), and occupational functioning (OR = 1.82, P = 0.0003) within the negative symptoms course, and dysphoric mood (OR = 2.00, P = 0.002) and stress intolerance (OR = 1.76, P = 0.0002) within the general symptoms course. Genetic liability for SCZ was additionally associated with full scale cognitive decline (β = -0.25, P = 0.02) and with longitudinal volumetric reduction of the right and left hippocampi (β = -0.28, P = 0.005; β = -0.23, P = 0.02, respectively). Our results indicate that the polygenic contribution to SCZ acts upon the threshold-lowering first hit (i.e., the deletion). It modifies the endophenotypes of 22q11.2DS and augments the derailment of developmental trajectories of negative and general symptoms, cognition, and hippocampal volume.

DCD-chip designed for the digital and ultraprecise quantification of copy number variation

The progress of CNV detection in DCD-chip.

Integrative analysis of prognostic long non-coding RNAs with copy number variation in bladder cancer

Copy number variations (CNVs), which can affect the role of long non-coding RNAs (lncRNAs), are important genetic changes seen in some malignant tumors. We analyzed lncRNAs with CNV to explore the relationship between lncRNAs and prognosis in bladder cancer (BLCA). Messenger RNA (mRNA) expression levels, DNA methylation, and DNA copy number data of 408 BLCA patients were subjected to integrative bioinformatics analysis. Cluster analysis was performed to obtain different subtypes and differently expressed lncRNAs and coding genes. Weighted gene co-expression network analysis (WGCNA) was performed to identify the co-expression gene and lncRNA modules. CNV-associated lncRNA data and their influence on cancer prognosis were assessed with Kaplan-Meier survival curve. Multi-omics integration analysis revealed five prognostic lncRNAs with CNV, namely NR2F1-AS1, LINC01138, THUMPD3-AS1, LOC101928489,and TMEM147-AS1,and a risk-score signature related to overall survival in BLCA was identified. Moreover, validated results in another independent Gene Expression Omnibus (GEO) dataset, GSE31684, were consistent with these results. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that the mitogen-activated protein kinase (MAPK) signaling pathway, focal adhesion pathway, and Janus kinase-signal transducers and activators of transcription (JAK-STAT) signaling pathway were enriched in a high-risk score pattern, suggesting that imbalance in these pathways is closely related to tumor development. We revealed the prognosis-related lncRNAs by analyzing the expression profiles of lncRNAs and CNVs, which can be used as prognostic biomarkers for BLCA.

Copy number variation and expression of exportin-4 associates with severity of fibrosis in metabolic associated fatty liver disease

Background

Liver fibrosis risk is a heritable trait, the outcome of which is the net deposition of extracellular matrix by hepatic stellate cell-derived myofibroblasts. Whereas nucleotide sequence variations have been extensively studied in liver fibrosis, the role of copy number variations (CNV) in which genes exist in abnormal numbers of copies (mostly due to duplication or deletion) has had limited exploration.

Methods

The impact of the XPO4 CNV on histological liver damage was examined in a cohort comprised 646 Caucasian patients with biopsy-proven MAFLD and 170 healthy controls. XPO4 expression was modulated and function was examined in human and animal models.

Findings

Here we demonstrate in a cohort of 816 subjects, 646 with biopsy-proven metabolic associated liver disease (MAFLD) and 170 controls, that duplication in the exportin 4 (XPO4) CNV is associated with the severity of liver fibrosis. Functionally, this occurs via reduced expression of hepatic XPO4 that maintains sustained activation of SMAD3/SMAD4 and promotes TGF-β1-mediated HSC activation and fibrosis. This effect was mediated through termination of nuclear SMAD3 signalling. XPO4 demonstrated preferential binding to SMAD3 compared to other SMADs and led to reduced SMAD3-mediated responses as shown by attenuation of TGFβ1 induced SMAD transcriptional activity, reductions in the recruitment of SMAD3 to target gene promoters following TGF-β1, as well as attenuation of SMAD3 phosphorylation and disturbed SMAD3/SMAD4 complex formation.

Interpretation

We conclude that a CNV in XPO4 is a critical mediator of fibrosis severity and can be exploited as a therapeutic target for liver fibrosis.

Funding

ME and JG are supported by the Robert W. Storr Bequest to the Sydney Medical Foundation, University of Sydney; a National Health and Medical Research Council of Australia (NHMRC) Program Grant (APP1053206) and Project and ideas grants (APP2001692, APP1107178 and APP1108422). AB is supported by an Australian Government Research Training Program (RTP) scholarship. EB is supported by Horizon 2020 under grant 634413 for the project EPoS.

Establishment of a novel CNV-related prognostic signature predicting prognosis in patients with breast cancer

Background

Copy number variation (CNVs) is a key factor in breast cancer development. This study determined prognostic molecular characteristics to predict breast cancer through performing a comprehensive analysis of copy number and gene expression data.

Methods

Breast cancer expression profiles, CNV and complete information from The Cancer Genome Atlas (TCGA) dataset were collected. Gene Expression Omnibus (GEO) chip data sets (GSE20685 and GSE31448) containing breast cancer samples were used as external validation sets. Univariate survival COX analysis, multivariate survival COX analysis, least absolute shrinkage and selection operator (LASSO), Chi square, Kaplan-Meier (KM) survival curve and receiver operating characteristic (ROC) analysis were applied to build a gene signature model and assess its performance.

Results

A total of 649 CNV related-differentially expressed gene obtained from TCGA-breast cancer dataset were related to several cancer pathways and functions. A prognostic gene sets with 9 genes were developed to stratify patients into high-risk and low-risk groups, and its prognostic performance was verified in two independent patient cohorts (n = 327, 246). The result uncovered that 9-gene signature could independently predict breast cancer prognosis. Lower mutation of PIK3CA and higher mutation of TP53 and CDH1 were found in samples with high-risk score compared with samples with low-risk score. Patients in the high-risk group showed higher immune score, malignant clinical features than those in the low-risk group. The 9-gene signature developed in this study achieved a higher AUC.

Conclusion

The current research established a 5-CNV gene signature to evaluate prognosis of breast cancer patients, which may innovate clinical application of prognostic assessment.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13048-021-00823-y.

CNVIntegrate: the first multi-ethnic database for identifying copy number variations associated with cancer

Human copy number variations (CNVs) and copy number alterations (CNAs) are DNA segments (>1000 base pairs) of duplications or deletions with respect to the reference genome, potentially causing genomic imbalance leading to diseases such as cancer. CNVs further cause genetic diversity in healthy populations and are predominant drivers of gene/genome evolution. Initiatives have been taken by the research community to establish large-scale databases to comprehensively characterize CNVs in humans. Exome Aggregation Consortium (ExAC) is one such endeavor that catalogs CNVs, of nearly 60 000 healthy individuals across five demographic clusters. Furthermore, large projects such as the Catalogue of Somatic Mutations in Cancer (COSMIC) and the Cancer Cell Line Encyclopedia (CCLE) combine CNA data from cancer-affected individuals and large panels of human cancer cell lines, respectively. However, we lack a structured and comprehensive CNV/CNA resource including both healthy individuals and cancer patients across large populations. CNVIntegrate is the first web-based system that hosts CNV and CNA data from both healthy populations and cancer patients, respectively, and concomitantly provides statistical comparisons between copy number frequencies of multiple ethnic populations. It further includes, for the first time, well-cataloged CNV and CNA data from Taiwanese healthy individuals and Taiwan Breast Cancer data, respectively, along with imported resources from ExAC, COSMIC and CCLE. CNVIntegrate offers a CNV/CNA-data hub for structured information retrieval for clinicians and scientists towards important drug discoveries and precision treatments. Database URL: http://cnvintegrate.cgm.ntu.edu.tw/.

Genome-wide copy number variations in a large cohort of bantu African children

BACKGROUND

Copy number variations (CNVs) account for a substantial proportion of inter-individual genomic variation. However, a majority of genomic variation studies have focused on single-nucleotide variations (SNVs), with limited genome-wide analysis of CNVs in large cohorts, especially in populations that are under-represented in genetic studies including people of African descent.

METHODS

We carried out a genome-wide copy number analysis in > 3400 healthy Bantu Africans from Tanzania. Signal intensity data from high density (> 2.5 million probes) genotyping arrays were used for CNV calling with three algorithms including PennCNV, DNAcopy and VanillaICE. Stringent quality metrics and filtering criteria were applied to obtain high confidence CNVs.

RESULTS

We identified over 400,000 CNVs larger than 1 kilobase (kb), for an average of 120 CNVs (SE = 2.57) per individual. We detected 866 large CNVs (≥ 300 kb), some of which overlapped genomic regions previously associated with multiple congenital anomaly syndromes, including Prader-Willi/Angelman syndrome (Type1) and 22q11.2 deletion syndrome. Furthermore, several of the common CNVs seen in our cohort (≥ 5%) overlap genes previously associated with developmental disorders.

CONCLUSIONS

These findings may help refine the phenotypic outcomes and penetrance of variations affecting genes and genomic regions previously implicated in diseases. Our study provides one of the largest datasets of CNVs from individuals of African ancestry, enabling improved clinical evaluation and disease association of CNVs observed in research and clinical studies in African populations.

Genetic Variants of Alcohol Metabolizing Enzymes and Alcohol-Related Liver Cirrhosis Risk

Alcohol-related liver disease (ARLD) is a major public health issue caused by excessive alcohol consumption. ARLD encompasses a wide range of chronic liver lesions, alcohol-related liver cirrhosis being the most severe and harmful state. Variations in the genes encoding the enzymes, which play an active role in ethanol metabolism, might influence alcohol exposure and hence be considered as risk factors of developing cirrhosis. We conducted a case-control study in which 164 alcohol-related liver cirrhosis patients and 272 healthy controls were genotyped for the following functional single nucleotide variations (SNVs): ADH1B gene, rs1229984, rs1041969, rs6413413, and rs2066702; ADH1C gene, rs35385902, rs283413, rs34195308, rs1693482, and rs35719513; CYP2E1 gene, rs3813867. Furthermore, copy number variations (CNVs) for ADH1A, ADH1B, ADH1C, and CYP2E1 genes were analyzed. A significant protective association with the risk of developing alcohol-related liver cirrhosis was observed between the mutant alleles of SNVs ADH1B rs1229984 (P_c value = 0.037) and ADH1C rs283413 (P_c value = 0.037). We identified CNVs in all genes studied, ADH1A gene deletions being more common in alcohol-related liver cirrhosis patients than in control subjects, although the association lost statistical significance after multivariate analyses. Our findings support that susceptibility to alcohol-related liver cirrhosis is related to variations in alcohol metabolism genes.

TPQCI: A topology potential-based method to quantify functional influence of copy number variations

Copy number variation (CNV) is a major type of chromosomal structural variation that play important roles in many diseases including cancers. Due to genome instability, a large number of CNV events can be detected in diseases such as cancer. Therefore, it is important to identify the functionally important CNVs in diseases, which currently still poses a challenge in genomics. One of the critical steps to solve the problem is to define the influence of CNV. In this paper, we provide a topology potential based method, TPQCI, to quantify this kind of influence by integrating statistics, gene regulatory associations, and biological function information. We used this metric to detect functionally enriched genes on genomic segments with CNV in breast cancer and multiple myeloma and discovered biological functions influenced by CNV. Our results demonstrate that, by using our proposed TPQCI metric, we can detect disease-specific genes that are influenced by CNVs. Source codes of TPQCI are provided in Github (https://github.com/usos/TPQCI).

The Copy Number Variation of OsMTD1 Regulates Rice Plant Architecture

Copy number variation (CNV) may have phenotypic effects by altering the expression level of the gene(s) or regulatory element(s) contained. It is believed that CNVs play pivotal roles in controlling plant architecture and other traits in plant. However, the effects of CNV contributing to special traits remain largely unknown. Here we report a CNV involved in rice architecture by modulating tiller number and leaf angle. In the genome of Oryza sativa ssp. japonica cv. Nipponbare, we found a locus Loc_Os08g34249 is derived from a 13,002-bp tandem duplication in the nearby region of OsMTD1, a gene regulating tillering in rice. Further survey of 230 rice cultivars showed that the duplication occurred in only 13 japonica rice cultivars. Phenotypic investigation indicated that this CNV region may contribute to tiller number. Moreover, we revealed that OsMTD1 not only influences rice tiller number and leaf angle, but also represses pri-miR156f transcription in the CNV region. Intriguingly, this CNV performs function through both the dosage and position effects on OsMTD1 and pri-miR156f. Thus, our work identified a CNV and revealed a molecular regulatory basis for its effects on plant architecture, implying this CNV may possess importance and application potential in molecular breeding in rice.

Molecular Mechanisms of MYCN Dysregulation in Cancers

MYCN, a member of MYC proto-oncogene family, encodes a basic helix-loop-helix transcription factor N-MYC. Abnormal expression of N-MYC is correlated with high-risk cancers and poor prognosis. Initially identified as an amplified oncogene in neuroblastoma in 1983, the oncogenic effect of N-MYC is expanded to multiple neuronal and nonneuronal tumors. Direct targeting N-MYC remains challenge due to its "undruggable" features. Therefore, alternative therapeutic approaches for targeting MYCN-driven tumors have been focused on the disruption of transcription, translation, protein stability as well as synthetic lethality of MYCN. In this review, we summarize the latest advances in understanding the molecular mechanisms of MYCN dysregulation in cancers.

Copy-number variation of the NPHP1 gene in patients with juvenile Nephronophthisis

Objective: Juvenile nephronophthisis (NPHP) is an autosomal recessive cystic disease of the kidney. It represents the most frequent genetic cause of chronic renal failure in children. Methods: we investigated clinical and molecular features in two children with Juvenile nephronophthisis using firstly Multiplex ligation-dependent probe amplification (MLPA) and secondly multiplex PCR. Results: we report a homozygous NPHP1 deletion in two children. Conclusion: NPHP1 deletion analysis using diagnostic methods (e.g. MLPA, Multiplex PCR) should always be considered in patients with nephronophthisis, especially from consanguineous families. Our results provide insights into genotype-phenotype correlations in juvenile nephronophthisis that can be utilized in genetic counseling.

RNA m6A Methylation Regulators Subclassify Luminal Subtype in Breast Cancer

RNA N⁶-methyladenosine (m⁶A) methylation is the most prevalent epitranscriptomic modification in mammals, with a complex and fine-tuning regulatory system. Recent studies have illuminated the potential of m⁶A regulators in clinical applications including diagnosis, therapeutics, and prognosis. Based on six datasets of breast cancer in The Cancer Genome Atlas (TCGA) database and two additional proteomic datasets, we provide a comprehensive view of all the known m⁶A regulators in their gene expression, copy number variations (CNVs), DNA methylation status, and protein levels in breast tumors and their association with prognosis. Among four breast cancer subtypes, basal-like subtype exhibits distinct expression and genomic alteration in m⁶A regulators from other subtypes. Accordingly, four representative regulators (IGF2BP2, IGF2BP3, YTHDC2, and RBM15) are identified as basal-like subtype-featured genes. Notably, luminal A/B samples are subclassified into two clusters based on the methylation status of those four genes. In line with its similarity to basal-like subtype, cluster1 shows upregulation in immune-related genes and cell adhesion molecules, as well as an increased number of tumor-infiltrating lymphocytes. Besides, cluster1 has worse disease-free and progression-free survival, especially among patients diagnosed with stage II and luminal B subtype. Together, this study highlights the potential functions of m⁶A regulators in the occurrence and malignancy progression of breast cancer. Given the heterogeneity within luminal subtype and high risk of recurrence and metastasis in a portion of patients, the prognostic stratification of luminal A/B subtypes utilizing basal-featured m⁶A regulators may help to improve the accuracy of diagnosis and therapeutics of breast cancer.

Genome-wide association study between copy number variation regions and carcass- and meat-quality traits in Nellore cattle

Context Indicine breeds are the main source of beef products in tropical and subtropical regions. However, genetic improvement for carcass- and meat-quality traits in zebu cattle have been limited and genomics studies concerning structural variations that influence these traits are essential. Aim The aim of this study was to perform a genome-wide association study between copy number variation regions (CNVRs) and carcass- and meat quality-traits in Nellore cattle. Methods In total, 3794 animals, males and females included, were genotyped using a 777962 single-nucleotide polymorphism platform of BovineHD BeadChip (777k; Illumina Inc.). Of these, 1751 Nellore bulls were slaughtered at 24 months of age for further carcass beef analysis. The following traits were studied: beef tenderness, marbling, rib-eye area, backfat thickness and meat colour (lightness, redness and yellowness). The CNV detection was conducted through PennCNV software. The association analyses were performed using CNVRuler software. Key results Several identified genomic regions were linked to quantitative trait loci associated with fat deposition (FABP7) and lipid metabolism (PPARA; PLA2 family; BCHE), extracellular matrix (INS; COL10A1), contraction (SLC34A3; TRDN) and muscle development (CAPZP). The gene-enrichment analyses highlighted biological mechanisms directly related to the metabolism and synthesis of lipids and fatty acids. Conclusions The large number of potential candidate genes identified within the CNVRs, as well as the functions and pathways identified, should help better elucidate the genetic mechanisms involved in the expression of beef and carcass traits in Nellore cattle. Several CNVRs harboured genes that might have a functional impact to improve the beef and carcass traits. Implications The results obtained contribute to upgrade the sensorial and organoleptic attributes of Nellore cattle and make feasible the genetic improvement of carcass- and meat-quality traits.

Genomic Copy Number Variation Study of Nine Macaca Species Provides New Insights into Their Genetic Divergence, Adaptation, and Biomedical Application

Copy number variation (CNV) can promote phenotypic diversification and adaptive evolution. However, the genomic architecture of CNVs among Macaca species remains scarcely reported, and the roles of CNVs in adaptation and evolution of macaques have not been well addressed. Here, we identified and characterized 1,479 genome-wide hetero-specific CNVs across nine Macaca species with bioinformatic methods, along with 26 CNV-dense regions and dozens of lineage-specific CNVs. The genes intersecting CNVs were overrepresented in nutritional metabolism, xenobiotics/drug metabolism, and immune-related pathways. Population-level transcriptome data showed that nearly 46% of CNV genes were differentially expressed across populations and also mainly consisted of metabolic and immune-related genes, which implied the role of CNVs in environmental adaptation of Macaca. Several CNVs overlapping drug metabolism genes were verified with genomic quantitative polymerase chain reaction, suggesting that these macaques may have different drug metabolism features. The CNV-dense regions, including 15 first reported here, represent unstable genomic segments in macaques where biological innovation may evolve. Twelve gains and 40 losses specific to the Barbary macaque contain genes with essential roles in energy homeostasis and immunity defense, inferring the genetic basis of its unique distribution in North Africa. Our study not only elucidated the genetic diversity across Macaca species from the perspective of structural variation but also provided suggestive evidence for the role of CNVs in adaptation and genome evolution. Additionally, our findings provide new insights into the application of diverse macaques to drug study.

A Comprehensive Analysis of Alterations in DNA Damage Repair Pathways Reveals a Potential Way to Enhance the Radio-Sensitivity of Esophageal Squamous Cell Cancer

Esophageal squamous cell cancer (ESCC) is a common malignancy with a poor 5-year overall survival in China. Altered DNA damage repair (DDR) pathways are associated with a predisposition to cancer and contribute to therapeutic response and resistance in cancers. However, alterations of DDR pathway genes in ESCC are still largely unknown. In this study, we employed genome sequencing data of 192 samples, comparative genomic hybridization data of 123 cases, and gene expression microarray data of 119 patients to firstly perform a comprehensive analysis of the gene alterations of 7 DDR pathways in ESCC. Gene mutations and copy number variations (CNVs) were observed in all 7 DDR pathways, and especially, CNVs were the dominant alteration types. Compared with other pathways, two DNA double-strand break (DSB) repair pathways homologous recombination (HR) and non-homologous end joining (NHEJ), carried significant gene mutations and CNVs especially gene amplifications. Most genes including RAD54B, NBS1, RAD51B, and PRKDC were significantly amplified and over-expressed in ESCC. Amplification and high expression of DSB repair pathway genes were associated with poorer overall survival. Gene set variation analysis further showed that DSB repair pathways were up-regulated in ESCC. Besides, we firstly demonstrated that combination of mirin and NU7441, two inhibitors for HR and NHEJ respectively, with ionizing radiation treatment significantly enhanced DSBs, reduced clonogenic cell survival, inhibited cell proliferation, and promoted cell apoptosis in ESCC cells with DSB pathway gene amplification. These findings suggest that DSB repair pathways were significantly altered in ESCC and inhibiting DSB repair pathways might enhance the radio-sensitivity of ESCC with DSB repair up-regulation.

Development of MR-based preoperative nomograms predicting DNA copy number subtype in lower grade gliomas with prognostic implication

OBJECTIVES

We aimed to determine the value of MR-based preoperative nomograms in predicting DNA copy number (CN) subtype in lower grade glioma (LGG) patients.

METHODS

The overall survival (OS) data were analyzed. MRI data of 170 subjects were retrospectively analyzed. The correlation was explored by univariate and multivariate regression analysis.

RESULTS

CN2 subtype was associated with shortest median OS (CN2 subtype vs. others: 46.8 vs. 221.7 months, p < 0.05). The time-dependent receiver operating characteristic for the CN2 subtype was 0.80 (95% CI: 0.74-0.85) for survival at 1 year, 0.80 (95% CI: 0.75-0.85) for survival at 2 years, and 0.77 (95% CI: 0.73-0.83) for survival at 3 years. On multivariate analysis, hemorrhage (OR: 0.118; p < 0.001; 95% CI: 0.037-0.376), poorly defined margin (OR: 4.592; p < 0.001; 95% CI: 1.965-10.730), extranodular growth (OR: 0.247; p = 0.006; 95% CI: 0.091-0.671), and volume ≥ 60 cm³ (OR: 4.734.256; p < 0.001; 95% CI: 2.051-10.924) were associated with CN1 subtype (AUC: 0.781). Proportion CE tumor (OR: 5.905; p = 0.007; 95% CI: 1.622-21.493), extranodular growth (OR: 9.047; p = 0.001; 95% CI: 2.349-34.846), width ≥ median (OR: 0.231; p = 0.049; 95% CI: 0.054-0.998), and depth ≥ median (OR: 0.192; p = 0.023; 95% CI: 0.046-0.799) were associated with CN2 subtype (AUC: 0.854). Necrosis/cystic (OR: 6.128; p = 0.007; 95% CI: 1.635-22.968), hemorrhage (OR: 5.752; p = 0.002; 95% CI: 1.953-16.942), poorly defined margin (OR: 0.164; p < 0.001; 95% CI: 0.063-0.427), and volume ≥ median (OR: 4.422; p < 0.001; 95% CI: 1.925-10.160) were associated with CN3 subtype (AUC: 0.808). All three nomograms showed good discrimination and calibration. Decision curve analysis supported that all nomograms were clinically useful. The average accuracy of the tenfold cross-validation was 0.680 (CN1), 0.794 (CN2), and 0.894 (CN3), respectively.

CONCLUSIONS

The shortest OS was observed in patients with CN2 subtype. This preliminary radiogenomics analysis revealed that the MR-based preoperative nomograms provide individualized prediction of DNA copy number subtype in LGG patients.

KEY POINTS

• This preliminary radiogenomics analysis of LGG revealed that the MR-based preoperative nomograms provide individualized prediction of DNA copy number subtype in LGG patients. • The AUC for the ROC curve was 0.781 for CN1 subtype, 0.854 for CN2 subtype, and 0.808 for CN3 subtype. Decision curve analysis supported that all nomograms were clinically useful. • The sensitivity was 0.779 (CN1), 0.731 (CN2), and 0.851 (CN3), respectively. The specificity was 0.664 (CN1), 0.872 (CN2), and 0.625 (CN3), respectively. And the accuracy was 0.717 (CN1), 0.849 (CN2), and 0.692 (CN3), respectively.

A multiple genomic data fused SF2 prediction model, signature identification, and gene regulatory network inference for personalized radiotherapy

Radiotherapy is one of the most important cancer treatments, but its response varies greatly among individual patients. Therefore, the prediction of radiosensitivity, identification of potential signature genes, and inference of their regulatory networks are important for clinical and oncological reasons. Here, we proposed a novel multiple genomic fused partial least squares deep regression method to simultaneously analyze multi-genomic data. Using 60 National Cancer Institute cell lines as examples, we aimed to identify signature genes by optimizing the radiosensitivity prediction model and uncovering regulatory relationships. A total of 113 signature genes were selected from more than 20,000 genes. The root mean square error of the model was only 0.0025, which was much lower than previously published results, suggesting that our method can predict radiosensitivity with the highest accuracy. Additionally, our regulatory network analysis identified 24 highly important ‘hub’ genes. The data analysis workflow we propose provides a unified and computational framework to harness the full potential of large-scale integrated cancer genomic data for integrative signature discovery. Furthermore, the regression model, signature genes, and their regulatory network should provide a reliable quantitative reference for optimizing personalized treatment options, and may aid our understanding of cancer progress mechanisms.

Copy number variation of IL17RA gene and its association with the ankylosing spondylitis risk in Iranian patients: a case-control study

Background

Ankylosing spondylitis (AS) is considered as a subtype of spondyloarthritis (SpA) that mainly leads to fatigue, stiffness, spinal ankylosis, and impaired physical functions with reduced quality of life. Interleukin (IL)-17A provokes additional inflammatory mediators and recruits immune cells to the inflamed site. IL17 expression increased in various inflammatory disorders including psoriasis, rheumatoid arthritis, multiple sclerosis, crohn’s disease, and ankylosing spondylitis. The current study aimed to evaluate the association of IL17RA copy number changes with the susceptibility to AS and their correlation to IL17RA expression in Iranian population.

Methods

IL17RA copy number genotyping assessments were carried out in 455 AS patients and 450 healthy controls, using custom TaqMan CNV assays. TaqMan primers and probe were located in Chr.22:17109553 based on pre-designed IL17RA Copy Number Assay ID, Hs02339506_cn. mRNA expression of IL17RA was also measured by SYBR Green real-time polymerase chain reaction (PCR).

Results

A IL17RA copy number loss (< 2) was associated with AS compared to 2 copies as reference (OR:2.18, 95% CI: (1.38–3.44), P-value < 0.001) and increased the risk of AS. IL17RA mRNA expression showed a significant increase in peripheral blood mononuclear cells (PBMCs) of all AS individuals than controls. The mRNA expression level of 2 copies was significantly higher in AS patients.

Conclusions

Our findings revealed that a low copy number of IL17RA might confer a susceptibility risk to AS. However, it is probably not directly involved in the regulation of IL17RA mRNA expression. Epigenetic mechanisms like DNA methylation, post-transcriptional, and -translational modifications that regulate the expression of the genes may contribute in upregulation of IL17RA mRNA expression in the loss of gene copy number condition.

Computational identification and characterization of glioma candidate biomarkers through multi-omics integrative profiling

BACKGROUND

Glioma is one of the most common malignant brain tumors and exhibits low resection rate and high recurrence risk. Although a large number of glioma studies powered by high-throughput sequencing technologies have led to massive multi-omics datasets, there lacks of comprehensive integration of glioma datasets for uncovering candidate biomarker genes.

RESULTS

In this study, we collected a large-scale assemble of multi-omics multi-cohort datasets from worldwide public resources, involving a total of 16,939 samples across 19 independent studies. Through comprehensive molecular profiling across different datasets, we revealed that PRKCG (Protein Kinase C Gamma), a brain-specific gene detectable in cerebrospinal fluid, is closely associated with glioma. Specifically, it presents lower expression and higher methylation in glioma samples compared with normal samples. PRKCG expression/methylation change from high to low is indicative of glioma progression from low-grade to high-grade and high RNA expression is suggestive of good survival. Importantly, PRKCG in combination with MGMT is effective to predict survival outcomes in a more precise manner.

CONCLUSIONS

PRKCG bears the great potential for glioma diagnosis, prognosis and therapy, and PRKCG-like genes may represent a set of important genes associated with different molecular mechanisms in glioma tumorigenesis. Our study indicates the importance of computational integrative multi-omics data analysis and represents a data-driven scheme toward precision tumor subtyping and accurate personalized healthcare.