Candidate Luminal B Breast Cancer Genes Identified by Genome,Gene Expression and DNA Methylation Profiling

Breast cancers (BCs) of the luminal B subtype are estrogen receptor-positive (ER+), highly proliferative, resistant to standard therapies and have a poor prognosis. To better understand this subtype we compared DNA copy number aberrations (CNAs), DNA promoter methylation, gene expression profiles, and somatic mutations in nine selected genes, in 32 luminal B tumors with those observed in 156 BCs of the other molecular subtypes. Frequent CNAs included 8p11-p12 and 11q13.1-q13.2 amplifications, 7q11.22-q34, 8q21.12-q24.23, 12p12.3-p13.1, 12q13.11-q24.11, 14q21.1-q23.1, 17q11.1-q25.1, 20q11.23-q13.33 gains and 6q14.1-q24.2, 9p21.3-p24,3, 9q21.2, 18p11.31-p11.32 losses. A total of 237 and 101 luminal B-specific candidate oncogenes and tumor suppressor genes (TSGs) presented a deregulated expression in relation with their CNAs, including 11 genes previously reported associated with endocrine resistance. Interestingly, 88% of the potential TSGs are located within chromosome arm 6q, and seven candidate oncogenes are potential therapeutic targets. A total of 100 candidate oncogenes were validated in a public series of 5,765 BCs and the overexpression of 67 of these was associated with poor survival in luminal tumors. Twenty-four genes presented a deregulated expression in relation with a high DNA methylation level. FOXO3, PIK3CA and TP53 were the most frequent mutated genes among the nine tested. In a meta-analysis of next-generation sequencing data in 875 BCs, KCNB2 mutations were associated with luminal B cases while candidate TSGs MDN1 (6q15) and UTRN (6q24), were mutated in this subtype. In conclusion, we have reported luminal B candidate genes that may play a role in the development and/or hormone resistance of this aggressive subtype.     

De novo genome assembly of the stress tolerant forest species Casuarina equisetifolia provides insight into secondary growth

Casuarina equisetifolia (C. equisetifolia), a conifer‐like angiosperm with resistance to typhoon and stress tolerance, is mainly cultivated in the coastal areas of Australasia. C. equisetifolia, making it a valuable model to study secondary growth associated genes and stress‐tolerance traits. However, the genome sequence is unavailable and therefore wood‐associated growth rate and stress resistance at the molecular level is largely unexplored. We therefore constructed a high‐quality draft genome sequence of C. equisetifolia by a combination of Illumina second‐generation sequencing reads and Pacific Biosciences single‐molecule real‐time (SMRT) long reads to advance the investigation of this species. Here, we report the genome assembly, which contains approximately 300 megabases (Mb) and scaffold size of N50 is 1.06 Mb. Additionally, gene annotation, assisted by a combination of prediction and RNA‐seq data, generated 29 827 annotated protein‐coding genes and 1983 non‐coding genes, respectively. Furthermore, we found that the total number of repetitive sequences account for one‐third of the genome assembly. Here we also construct the genome‐wide map of DNA modification, such as two novel forms N⁶‐adenine (6mA) and N4‐methylcytosine (4mC) at the level of single‐nucleotide resolution using single‐molecule real‐time (SMRT) sequencing. Interestingly, we found that 17% of 6mA modification genes and 15% of 4mC modification genes also included alternative splicing events. Finally, we investigated cellulose, hemicellulose, and lignin‐related genes, which were associated with secondary growth and contained different DNA modifications. The high‐quality genome sequence and annotation of C. equisetifolia in this study provide a valuable resource to strengthen our understanding of the diverse traits of trees.     

A global view of the oncogenic landscape in nasopharyngeal carcinoma: an integrated analysis at the genetic and expression levels

Previous studies have reported that the tumour cells of nasopharyngeal carcinoma (NPC) exhibit recurrent chromosome abnormalities. These genetic changes are broadly assumed to lead to changes in gene expression which are important for the pathogenesis of this tumour. However, this assumption has yet to be formally tested at a global level. Therefore a genome wide analysis of chromosome copy number and gene expression was performed in tumour cells micro-dissected from the same NPC biopsies. Cellular tumour suppressor and tumour-promoting genes (TSG, TPG) and Epstein-Barr Virus (EBV)-encoded oncogenes were examined. The EBV-encoded genome maintenance protein EBNA1, along with the putative oncogenes LMP1, LMP2 and BARF1 were expressed in the majority of NPCs that were analysed. Significant downregulation of expression in an average of 76 cellular TSGs per tumour was found, whilst a per-tumour average of 88 significantly upregulated, TPGs occurred. The expression of around 60% of putative TPGs and TSGs was both up-and down-regulated in different types of cancer, suggesting that the simplistic classification of genes as TSGs or TPGs may not be entirely appropriate and that the concept of context-dependent onco-suppressors may be more extensive than previously recognised. No significant enrichment of TPGs within regions of frequent genomic gain was seen but TSGs were significantly enriched within regions of frequent genomic loss. It is suggested that loss of the FHIT gene may be a driver of NPC tumourigenesis. Notwithstanding the association of TSGs with regions of genomic loss, on a gene by gene basis and excepting homozygous deletions and high-level amplification, there is very little correlation between chromosomal copy number aberrations and expression levels of TSGs and TPGs in NPC.     

Exploring genome-wide DNA methylation profiles altered in hepatocellular carcinoma using Infinium HumanMethylation 450 BeadChips

Hepatocellular carcinoma (HCC) incidence has increased in the US and also has one of the fastest growing death rates of any cancer. The purpose of the current study was to discover novel genome-wide aberrant DNA methylation patterns in HCC tumors that are predominantly HCV-related. Infinium HumanMethylation 450K BeadChip arrays were used to examine genome-wide DNA methylation profiles in 66 pairs of HCC tumor and adjacent non-tumor tissues. After Bonferroni adjustment, a total of 130,512 CpG sites significantly differed in methylation level in tumor compared with non-tumor tissues, with 28,017 CpG sites hypermethylated and 102,495 hypomethylated in tumor tissues. Absolute tumor/non-tumor methylation differences ≥ 20% were found in 24.9% of the hypermethylated and 43.1% of the hypomethylated CpG sites; almost 10,000 CpG sites have ≥ 30% DNA methylation differences. Most (60.1%) significantly hypermethylated CpG sites are located in CpG islands, with 21.6% in CpG shores and 3.6% in shelves. In contrast, only a small proportion (8.2%) of significantly hypomethylated CpG sites are situated in islands, while most are found in open sea (60.2%), shore (17.3%) or shelf (14.3%) regions. A total of 2,568 significant CpG sites (2,441 hypermethylated and 127 hypomethylated) covering 589 genes are located within 684 differentially methylated regions defined as regions with at least two significant CpG sites displaying > 20% methylation differences in the same direction within 250-bp. The top 500 significant CpG sites can significantly distinguish HCC tumor from adjacent tissues with one misclassification. Within adjacent non-tumor tissues, we also identified 75 CpG sites significantly associated with gender, 228 with HCV infection, 17,207 with cirrhosis, and 56 with both HCV infection and cirrhosis after multiple comparisons adjustment. Aberrant DNA methylation profiles across the genome were identified in tumor tissues from US HCC cases that are predominantly related to HCV infection. These results demonstrate the significance of aberrant DNA methylation in HCC tumorigenesis.     

Identification of genomic functional hotspots with copy number alteration in liver cancer

Copy number alterations (CNAs) can be observed in most of cancer patients. Several oncogenes and tumor suppressor genes with CNAs have been identified in different kinds of tumor. However, the systematic survey of CNA-affected functions is still lack. By employing systems biology approaches, instead of examining individual genes, we directly identified the functional hotspots on human genome. A total of 838 hotspots on human genome with 540 enriched Gene Ontology functions were identified. Seventy-six aCGH array data of hepatocellular carcinoma (HCC) tumors were employed in this study. A total of 150 regions which putatively affected by CNAs and the encoded functions were identified. Our results indicate that two immune related hotspots had copy number alterations in most of patients. In addition, our data implied that these immune-related regions might be involved in HCC oncogenesis. Also, we identified 39 hotspots of which copy number status were associated with patient survival. Our data implied that copy number alterations of the regions may contribute in the dysregulation of the encoded functions. These results further demonstrated that our method enables researchers to survey biological functions of CNAs and to construct regulation hypothesis at pathway and functional levels.     

Identification and Validation of Novel Immunogenic Cell Death- and DNA Damage Response-Related Molecular Patterns Correlated with Immune Status and Prognosis in Hepatocellular Carcinoma

Immunogenic cell death (ICD) and DNA damage response (DDR) are involved in cancer progression and prognosis. Currently, chemotherapy is the first-line treatment for intermediate or advanced hepatocellular carcinoma (HCC), which is mostly based on platinum and anthracyclines that induce DNA damage and ICD. With the treatment of HCC with immune checkpoint inhibitors (ICIs), it is important to understand the molecular characteristics and prognostic values of ICD and DDR-related genes (IDRGs). We aimed to explore the characteristics of ICD and DDR-related molecular patterns, immune status, and the association of immunotherapy and prognosis with IDRGs in HCC. We identified IDRGs in HCC and evaluated their differential expression, biological behaviors, molecular characteristics, immune cell infiltration, and prognostic value. Prognostic IDRGs and subtypes were identified and validated. FFAR3, DDX1, POLR3G, FANCL, ADA, PI3KR1, DHX58, TPT1, MGMT, SLAMF6, and EIF2AK4 were determined as risk factors for HCC, and the biological experiments indicated that high FANCL expression is harmful to the treatment and prognosis. HCC was classified into high- and low-risk groups based on the median values of the risk factors to construct a predictive nomogram. These findings provide novel insights into the treatment and prognosis of HCC and provide a new research direction for HCC.     

Localized structural frustration for evaluating the impact of sequence variants

Population-scale sequencing is increasingly uncovering large numbers of rare single-nucleotide variants (SNVs) in coding regions of the genome. The rarity of these variants makes it challenging to evaluate their deleteriousness with conventional phenotype–genotype associations. Protein structures provide a way of addressing this challenge. Previous efforts have focused on globally quantifying the impact of SNVs on protein stability. However, local perturbations may severely impact protein functionality without strongly disrupting global stability (e.g. in relation to catalysis or allostery). Here, we describe a workflow in which localized frustration, quantifying unfavorable local interactions, is employed as a metric to investigate such effects. Using this workflow on the Protein Databank, we find that frustration produces many immediately intuitive results: for instance, disease-related SNVs create stronger changes in localized frustration than non-disease related variants, and rare SNVs tend to disrupt local interactions to a larger extent than common variants. Less obviously, we observe that somatic SNVs associated with oncogenes and tumor suppressor genes (TSGs) induce very different changes in frustration. In particular, those associated with TSGs change the frustration more in the core than the surface (by introducing loss-of-function events), whereas those associated with oncogenes manifest the opposite pattern, creating gain-of-function events.     

Genome-wide profiling of sperm DNA methylation in relation to buffalo (Bubalus bubalis) bull fertility

The DNA methylation pattern in spermatozoa of buffalo bulls of different fertility status was investigated. Spermatozoa isolated DNA from two groups of buffalo bulls (n = 5), selected based on their artificial insemination–generated conception rate data followed by IVF efficiency, were studied for global methylation changes using a custom-designed 180 K buffalo (Bubalus bubalis) CpG island/promoter microarray. A total of 96 individual genes with another 55 genes covered under CpG islands were found differentially methylated in sperm of high-fertile and subfertile buffalo bulls. Important genes associated with biological processes, cellular components, and functions were identified to be differentially methylated in buffalo bulls with differential fertility status. The identified differentially methylated genes were found to be involved in germ cell development, spermatogenesis, capacitation, and embryonic development. The observations hint that methylation defects of sperm DNA may play a crucial role in determining the fertility of breeding bulls. This growing field of sperm epigenetics will be of great benefit in understanding the graded fertility conditions of breeding bulls in commercial livestock production system.     

Genomewide discovery of DNA polymorphisms in rice cultivars with contrasting drought and salinity stress response and their functional relevance

Next‐generation sequencing technologies provide opportunities to understand the genetic basis of phenotypic differences, such as abiotic stress response, even in the closely related cultivars via identification of large number of DNA polymorphisms. We performed whole‐genome resequencing of three rice cultivars with contrasting responses to drought and salinity stress (sensitive IR64, drought‐tolerant Nagina 22 and salinity‐tolerant Pokkali). More than 356 million 90‐bp paired‐end reads were generated, which provided about 85% coverage of the rice genome. Applying stringent parameters, we identified a total of 1 784 583 nonredundant single‐nucleotide polymorphisms (SNPs) and 154 275 InDels between reference (Nipponbare) and the three resequenced cultivars. We detected 401 683 and 662 509 SNPs between IR64 and Pokkali, and IR64 and N22 cultivars, respectively. The distribution of DNA polymorphisms was found to be uneven across and within the rice chromosomes. One‐fourth of the SNPs and InDels were detected in genic regions, and about 3.5% of the total SNPs resulted in nonsynonymous changes. Large‐effect SNPs and InDels, which affect the integrity of the encoded protein, were also identified. Further, we identified DNA polymorphisms present in the differentially expressed genes within the known quantitative trait loci. Among these, a total of 548 SNPs in 232 genes, located in the conserved functional domains, were identified. The data presented in this study provide functional markers and promising target genes for salinity and drought tolerance and present a valuable resource for high‐throughput genotyping and molecular breeding for abiotic stress traits in rice.     

QTL identification of flowering time at three different latitudes reveals homeologous genomic regions that control flowering in soybean

Since the genetic control of flowering time is very important in photoperiod-sensitive soybean (Glycine max (L.) Merr.), genes affecting flowering under different environment conditions have been identified and described. The objectives were to identify quantitative trait loci (QTLs) for flowering time in different latitudinal and climatic regions, and to understand how chromosomal rearrangement and genome organization contribute to flowering time in soybean. Recombinant inbred lines from a cross between late-flowering ‘Jinpumkong 2’ and early-flowering ‘SS2-2’ were used to evaluate the phenotypic data for days to flowering (DF) collected from Kamphaeng Saen, Thailand (14°01′N), Suwon, Korea (37°15′N), and Longjing, China (42°46′N). A weakly positive phenotypic correlation (r = 0.36) was found between DF in Korea and Thailand; however, a strong correlation (r = 0.74) was shown between Korea and China. After 178 simple sequence repeat (SSR) markers were placed on a genetic map spanning 2,551.7 cM, four independent DF QTLs were identified on different chromosomes (Chrs). Among them, three QTLs on Chrs 9, 13 and 16 were either Thailand- or Korea-specific. The DF QTL on Chr 6 was identified in both Korea and China, suggesting it is less environment-sensitive. Comparative analysis of four DF QTL regions revealed a syntenic relationship between two QTLs on Chrs 6 and 13. All five duplicated gene pairs clustered in the homeologous genomic regions were found to be involved in the flowering. Identification and comparative analysis of multiple DF QTLs from different environments will facilitate the significant improvement in soybean breeding programs with respect to control of flowering time.     

Integration of DNA copy number alterations and transcriptional expression analysis in human gastric cancer

Background

Genomic instability with frequent DNA copy number alterations is one of the key hallmarks of carcinogenesis. The chromosomal regions with frequent DNA copy number gain and loss in human gastric cancer are still poorly defined. It remains unknown how the DNA copy number variations contributes to the changes of gene expression profiles, especially on the global level.

Principal Findings

We analyzed DNA copy number alterations in 64 human gastric cancer samples and 8 gastric cancer cell lines using bacterial artificial chromosome (BAC) arrays based comparative genomic hybridization (aCGH). Statistical analysis was applied to correlate previously published gene expression data obtained from cDNA microarrays with corresponding DNA copy number variation data to identify candidate oncogenes and tumor suppressor genes. We found that gastric cancer samples showed recurrent DNA copy number variations, including gains at 5p, 8q, 20p, 20q, and losses at 4q, 9p, 18q, 21q. The most frequent regions of amplification were 20q12 (7/72), 20q12–20q13.1 (12/72), 20q13.1–20q13.2 (11/72) and 20q13.2–20q13.3 (6/72). The most frequent deleted region was 9p21 (8/72). Correlating gene expression array data with aCGH identified 321 candidate oncogenes, which were overexpressed and showed frequent DNA copy number gains; and 12 candidate tumor suppressor genes which were down-regulated and showed frequent DNA copy number losses in human gastric cancers. Three networks of significantly expressed genes in gastric cancer samples were identified by ingenuity pathway analysis.

Conclusions

This study provides insight into DNA copy number variations and their contribution to altered gene expression profiles during human gastric cancer development. It provides novel candidate driver oncogenes or tumor suppressor genes for human gastric cancer, useful pathway maps for the future understanding of the molecular pathogenesis of this malignancy, and the construction of new therapeutic targets.     

Cross-Species Array Comparative Genomic Hybridization Identifies Novel Oncogenic Events in Zebrafish and Human Embryonal Rhabdomyosarcoma

Human cancer genomes are highly complex, making it challenging to identify specific drivers of cancer growth, progression, and tumor maintenance. To bypass this obstacle, we have applied array comparative genomic hybridization (array CGH) to zebrafish embryonal rhabdomyosaroma (ERMS) and utilized cross-species comparison to rapidly identify genomic copy number aberrations and novel candidate oncogenes in human disease. Zebrafish ERMS contain small, focal regions of low-copy amplification. These same regions were commonly amplified in human disease. For example, 16 of 19 chromosomal gains identified in zebrafish ERMS also exhibited focal, low-copy gains in human disease. Genes found in amplified genomic regions were assessed for functional roles in promoting continued tumor growth in human and zebrafish ERMS – identifying critical genes associated with tumor maintenance. Knockdown studies identified important roles for Cyclin D2 (CCND2), Homeobox Protein C6 (HOXC6) and PlexinA1 (PLXNA1) in human ERMS cell proliferation. PLXNA1 knockdown also enhanced differentiation, reduced migration, and altered anchorage-independent growth. By contrast, chemical inhibition of vascular endothelial growth factor (VEGF) signaling reduced angiogenesis and tumor size in ERMS-bearing zebrafish. Importantly, VEGFA expression correlated with poor clinical outcome in patients with ERMS, implicating inhibitors of the VEGF pathway as a promising therapy for improving patient survival. Our results demonstrate the utility of array CGH and cross-species comparisons to identify candidate oncogenes essential for the pathogenesis of human cancer.     

Whole genome DNA methylation profiling of oral cancer in ethnic population of Meghalaya,North East India reveals novel genes

Oral Squamous Cell Carcinoma (OSCC) is a serious and one of the most common and highly aggressive malignancies. Epigenetic factors such as DNA methylation have been known to be implicated in a number of cancer etiologies. The main objective of this study was to investigate physiognomies of Promoter DNA methylation patterns associated with oral cancer epigenome with special reference to the ethnic population of Meghalaya, North East India. The present study identifies 27,205 CpG sites and 3811 regions that are differentially methylated in oral cancer when compared to matched normal. 45 genes were found to be differentially methylated within the promoter region, of which 38 were hypermethylated and 7 hypomethylated. 14 of the hypermethylated genes were found to be similar to that of the TCGA-HNSCC study some of which are TSGs and few novel genes which may serve as candidate methylation biomarkers for OSCC in this poorly characterized ethnic group.