Analysis of lncRNA-Associated ceRNA Network Reveals Potential lncRNA Biomarkers in Human Colon Adenocarcinoma

Background/Aims: Long non-coding RNAs (lncRNAs) acting as competing endogenous RNAs (ceRNAs) play significant roles in the development of tumors, but the functions of specific lncRNAs and lncRNA-related ceRNA networks have not been fully elucidated for colon adenocarcinoma (COAD). In this study, we aimed to clarify the lncRNA-microRNA (miRNA)-mRNA ceRNA network and potential lncRNA biomarkers in COAD. Methods: We extracted data from The Cancer Genome Atlas (TCGA) and identified COAD-specific mRNAs, miRNAs, and lncRNAs. The biological processes in Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) were analyzed for COAD-specific mRNAs. We then constructed a ceRNA network of COAD-specific mRNAs, miRNAs and lncRNAs and analyzed the correlation between expression patterns and clinical features of the lncRNAs involved. After identifying potential mRNA targets of 4 lncRNAs related to overall survival (OS), we conducted stepwise analysis of these targets through GO and KEGG. Using tissue samples from our own patients, we also verified certain analytical results using quantitative real-time PCR (qRT-PCR). Results: Data from 521 samples (480 tumor tissue and 41 adjacent non-tumor tissue samples) were extracted from TCGA. A total of 258 specific lncRNAs, 206 specific miRNAs, and 1467 specific mRNAs were identified (absolute log2 [fold change] > 2, false discovery rate < 0.01). Analysis of KEGG revealed that specific mRNAs were enriched in cancer-related pathways. The ceRNA network was constructed with 64 lncRNAs, 18 miRNAs, and 42 mRNAs. Among these lncRNAs involved in the network, 3 lncRNAs (LINC00355, HULC, and IGF2-AS) were confirmed to be associated with certain clinical features and 4 lncRNAs (HOTAIR, LINC00355, KCNQ1OT1, and TSSC1-IT1) were found to be negatively linked to OS (log-rank p < 0.05). KEGG showed that the potential mRNA targets of these 4 lncRNAs may be concentrated in the MAPK pathway. Certain results were validated by qRT-PCR. Conclusion: This study providing novel insights into the lncRNA-miRNA-mRNA ceRNA network and reveals potential lncRNA biomarkers in COAD

An ALS-Linked Mutant SOD1 Produces a Locomotor Defect Associated with Aggregation and Synaptic Dysfunction When Expressed in Neurons of Caenorhabditis elegans

The nature of toxic effects exerted on neurons by misfolded proteins, occurring in a number of neurodegenerative diseases, is poorly understood. One approach to this problem is to measure effects when such proteins are expressed in heterologous neurons. We report on effects of an ALS-associated, misfolding-prone mutant human SOD1, G85R, when expressed in the neurons of Caenorhabditis elegans. Stable mutant transgenic animals, but not wild-type human SOD1 transgenics, exhibited a strong locomotor defect associated with the presence, specifically in mutant animals, of both soluble oligomers and insoluble aggregates of G85R protein. A whole-genome RNAi screen identified chaperones and other components whose deficiency increased aggregation and further diminished locomotion. The nature of the locomotor defect was investigated. Mutant animals were resistant to paralysis by the cholinesterase inhibitor aldicarb, while exhibiting normal sensitivity to the cholinergic agonist levamisole and normal muscle morphology. When fluorescently labeled presynaptic components were examined in the dorsal nerve cord, decreased numbers of puncta corresponding to neuromuscular junctions were observed in mutant animals and brightness was also diminished. At the EM level, mutant animals exhibited a reduced number of synaptic vesicles. Neurotoxicity in this system thus appears to be mediated by misfolded SOD1 and is exerted on synaptic vesicle biogenesis and/or trafficking

Data from: Effect of mutation mechanisms on variant composition and distribution in Caenorhabditis elegans

Genetic diversity is maintained by continuing generation and removal of variants. While examining over 800,000 DNA variants in wild isolates of Caenorhabditis elegans, we made a discovery that the proportions of variant types are not constant across the C. elegans genome. The variant proportion is defined as the fraction of a specific variant type (e.g. single nucleotide polymorphism (SNP) or indel) within a broader set of variants (e.g. all variants or all non-SNPs). The proportions of most variant types show a correlation with the recombination rate. These correlations can be explained as a result of a concerted action of two mutation mechanisms, which we named Morgan and Sanger mechanisms. The two proposed mechanisms act according to the distinct components of recombination rate, specifically the genetic and physical distance. Regression analysis was used to explore the characteristics and contributions of the two mutation mechanisms. According to our model, ~20–40% of all mutations in C. elegans wild populations are derived from programmed meiotic double strand breaks, which precede chromosomal crossovers and thus may be the point of origin for the Morgan mechanism. A substantial part of the known correlation between the recombination rate and variant distribution appears to be caused by the mutations generated by the Morgan mechanism. Mathematically integrating the mutation model with background selection model gives a more complete depiction of how the variant landscape is shaped in C. elegans. Similar analysis should be possible in other species by examining the correlation between the recombination rate and variant landscape within the context of our mutation model

Effect of recombination on genetic diversity of Caenorhabditis elegans

Abstract Greater molecular divergence and genetic diversity are present in regions of high recombination in many species. Studies describing the correlation between variant abundance and recombination rate have long focused on recombination in the context of linked selection models, whereby interference between linked sites under positive or negative selection reduces genetic diversity in regions of low recombination. Here, we show that indels, especially those of intermediate sizes, are enriched relative to single nucleotide polymorphisms in regions of high recombination in C. elegans. To explain this phenomenon, we reintroduce an alternative model that emphasizes the mutagenic effect of recombination. To extend the analysis, we examine the variants with a phylogenetic context and discuss how different models could be examined together. The number of variants generated by recombination in natural populations could be substantial including possibly the majority of some indel subtypes. Our work highlights the potential importance of a mutagenic effect of recombination, which could have a significant role in the shaping of natural genetic diversity

HwangWangPLOScomp.tar

Contains R scripts used for the analysis presented in the associated manuscript along with the source files needed to run the scripts

The landscape of variants by chromosomal position and recombination rate.

<p>(A-D) The proportions of i40-699 out of all variants (A, C) and SNPs out of all variants (B, D) are shown by chromosomal position (A, B) and relative to the recombination rate (C, D). (E) The distribution of all variants relative to the recombination rate. Black lines indicate the line of best fit by linear regression.</p

Fitting the correlation between the recombination rate and variant composition using the mutation model.

<p>The proportion of i40-699 out of all variants (A), SNPs out of all variants (B), i40-699 out of non-SNPs (C), and transitions (Ts) out of SNPs (D) for all segments with >300 variants outside exons are shown relative to the recombination rate (gray circles). Black circles indicate the mean proportions for different intervals. Red line indicates the line of best fit using all segments according to least squares regression using the mutation model.</p

Correlation between variant composition and recombination rate and other genomic features.

<p>(A) Box plot of Pearson's r absolute values for the correlation between variant composition and various genomic features. (B-E) Proportions of SNPs out of all variants (B, C) and of Ts out of all SNPs (D, E) relative to recombination rate for the variants that do not affect exons (B, D) and for the variants that affect exons (C, E).</p