Identification of Long Noncoding RNAs That Exert Transcriptional Regulation by Forming RNA-DNA Triplexes in Prostate Cancer.

Long noncoding RNAs (lncRNAs) are involved in transcriptional regulation, and their deregulation is associated with the development of various human cancers, including prostate cancer (PCa). However, their underlying mechanisms remain unclear. In this study, lncRNAs that interact with DNA and regulate mRNA transcription in PCa were screened and identified to promote PCa development. First, 4195 protein-coding genes (PCGs, mRNAs) were obtained from the The Cancer Genome Atlas (TCGA) database, in which 1148 lncRNAs were differentially expressed in PCa. Then, 44,270 pairs of co-expression relationships were calculated between 612 lncRNAs and 2742 mRNAs, of which 42,596 (96%) were positively correlated. Among the 612 lncRNAs, 392 had the potential to interact with the promoter region to form DNA:DNA:RNA triplexes, from which lncRNA AD000684.2(AC002128.1) was selected for further validation. AC002128.1 was highly expressed in PCa. Furthermore, AD000684.2 positively regulated the expression of the correlated genes. In addition, AD000684.2 formed RNA-DNA triplexes with the promoter region of the regulated genes. Functional assays also demonstrated that lncRNA AD000684.2 promotes cell proliferation and motility, as well as inhibits apoptosis, in PCa cell lines. The results suggest that AD000684.2 could positively regulate the transcription of target genes via triplex structures and serve as a candidate prognostic biomarker and target for new therapies in human PCa.

Androgen-Responsive Oncogenic lncRNA RP11-1023L17.1 Enhances c-Myc Protein Stability in Prostate Cancer.

Long noncoding RNAs (lncRNAs) have been found as novel participants in the pathophysiology of prostate cancer (PCa), which is predominantly regulated by androgen and its receptor. The biological function of androgen-responsive lncRNAs remains poorly understood. Here, we identified that lncRNA RP11-1023L17.1, which is highly expressed in PCa. RP11-1023L17.1 expression, can be directly repressed by the androgen receptor in PCa cells. RP11-1023L17.1 depletion inhibited the proliferation, migration, and cell cycle progression, and promoted the apoptosis of PCa cells, indicating that RP11-1023L17.1 acts as an oncogene in PCa cells. Microarray results revealed that RP11-1023L17.1 depletion downregulated the c-Myc transcription signature in PCa cells. RP11-1023L17.1 depletion-induced cellular phenotypes can be overcome by ectopically overexpressed c-Myc. Mechanistically, RP11-1023L17.1 represses FBXO32 mRNA expression, thereby enhancing c-Myc protein stability by blocking FBXO32-mediated c-Myc degradation. Our findings reveal the previously unrecognized roles of RP11-1023L17.1 in c-Myc-dependent PCa tumorigenesis.

Genome-wide identification of lncRNAs and mRNAs differentially expressed in non-functioning pituitary adenoma and construction of an lncRNA-mRNA co-expression network.

The involvement of long non-coding RNAs (lncRNAs) during tumorigenesis is a recent emerging theme. Yet no systematic evaluation of lncRNAs has been previously reported for non-functioning pituitary adenoma (NFPA), a fairly common type of intracranial tumor. Here, we report the first genome-wide expression profile for lncRNAs and mRNAs in NFPA, using formalin-fixed and paraffin-embedded tissue specimens. Using microarray analyses, we identified 113 lncRNAs and 80 mRNAs differentially expressed in NFPA; this list includes lncRNAs previously implicated in a variety of cancers. Using real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR) we further confirmed differential expression in NFPA for ten of the 113 lncRNAs. Using these ten doubly confirmed lncRNAs, we constructed an lncRNA-mRNA co-expression network comprising of 130 specific lncRNA-mRNA co-expression relationships. In addition, we conducted GO and KEGG analyses for the 80 mRNAs differentially expressed in NFPA. Our microarray and qRT-PCR analyses provided a working list of lncRNAs that may be functionally relevant to NFPA tumorigenesis. Our co-expression network in turn connected these largely uncharacterized lncRNAs to specific mRNAs, whose roles we further elucidated via GO and KEGG analyses, thus providing specific, testable hypotheses for the functions of these lncRNAs. Together, our study laid the foundation for future investigation of the specific function and mechanism by which lncRNAs are involved in NFPA tumorigenesis.

Expression profile reveals novel prognostic biomarkers in hepatocellular carcinoma.

The purpose of this study was to identify and validate novel prognostic biomarkers in human hepatocellular carcinoma (HCC). We analyzed gene expression profiles not only between 33 HCCs and their corresponding noncancerous liver tissues, but also between 25 HCCs and pooled normal liver tissues using cDNA microarrays containing 12800 genes. Functional analysis of differentially expressed genes involved in HCC carcinogenesis and tumor progression revealed that up-regulated and down-regulated genes are mainly associated with cell cycle and immune response, respectively. We detected two regions of cytogenetic changes only in poorly-differentiated HCCs using the expression data. We identified a 9-gene expression signature, which was able to predict differentiation degree and survival of HCC samples. Among the 9 most discriminatory genes, minichromosome maintenance protein 2 (MCM2), a significantly up-regulated gene involved in cell cycle pathway, was selected for further analysis. Overexpression of MCM2 protein related to poor-differentiation in HCC was validated using tissue microarray-based immunohistochemistry containing 96 HCCs. Our studies show that the 9-gene expression signature may serve as promising prognostic biomarkers involved in hepatocarcinogenesis and tumor progression.

System biology analysis of cell cycle pathway involved in hepatocellular carcinoma.

To investigate genetic mechanisms of hepatocarcinogenesis and identify potential anticancer targets in hepatocellular carcinoma (HCC), we analyzed microarray gene expression profiles between 33 HCCs and their corresponding noncancerous liver tissues. Functional analysis of differentially-expressed genes in HCC indicated that cell cycle dysregulation plays an important role in hepatocarcinogenesis. Based on 14 differentially-expressed genes involved in cell cycle in HCC, we applied Structural Equation Modeling (SEM) to establish a potential genetic network which could assist understanding of HCC molecular mechanisms. siRNA-mediated knock-down of two significantly up-regulated genes, minichromosome maintenance protein 2 (MCM2) and cyclin B1 (CCNB1), in HCC cells (SMMC-7721 and QGY-7703) induced G2/M-phase arrest, apoptosis and antiproliferation in HCC. Some up-regulated cell cycle-related genes in HCC were down-regulated following specific depletion of MCM2 or/and CCNB1 in HCC cells, which might well validate and complement the reconstructed cell cycle network. This study may contribute to further disclose hepatocarcinogenesis mechanism through systematically analyzed the HCC-related-cell-cycle pathway. This study also shows that MCM2 and CCNB1 could be promising prognostic and therapeutic targets for HCC.

Event-specific detection of seven genetically modified soybean and maizes using multiplex-PCR coupled with oligonucleotide microarray.

With the increasing development of genetically modified organism (GMO) detection techniques, the polymerase chain reaction (PCR) technique has been the mainstay for GMO detection. An oligonucleotide microarray is a glass chip to the surface of which an array of oligonucleotides was fixed as spots, each containing numerous copies of a sequence-specific probe that is complementary to a gene of interest. So it is used to detect ten or more targets synchronously. In this research, an event-specific detection strategy based on the unique and specific integration junction sequences between the host plant genome DNA and the integrated gene is being developed for its high specificity using multiplex-PCR together with oligonucleotide microarray. A commercial GM soybean (GTS 40-3-2) and six GM maize events (MON810, MON863, Bt176, Bt11, GA21, and T25) were detected by this method. The results indicate that it is a suitable method for the identification of these GM soybean and maizes.

Screening genetically modified organisms using multiplex-PCR coupled with oligonucleotide microarray.

In this research, we developed a multiplex polymerase chain reaction (multiplex-PCR) coupled with a DNA microarray system simultaneously aiming at many targets in a consecutive reaction to detect a genetically modified organism (GMO). There are a total of 20 probes for detecting a GMO in a DNA microarray which can be classified into three categories according to their purpose: the first for screening GMO from un-transgenic plants based on the common elements such as promoter, reporter and terminator genes; the second for specific gene confirmation based on the target gene sequences such as herbicide-resistance or insect-resistance genes; the third for species-specific genes which the sequences are unique for different plant species. To ensure the reliability of this method, different kinds of positive and negative controls were used in DNA microarray. Commercial GM soybean, maize, rapeseed and cotton were identified by means of this method and further confirmed by PCR analysis and sequencing. The results indicate that this method discriminates between the GMOs very quickly and in a cost-saving and more time efficient way. It can detect more than 95% of currently commercial GMO plants and the limits of detection are 0.5% for soybean and 1% for maize. This method is proved to be a new method for routine analysis of GMOs.

Analysis of the factors affecting the accuracy of detection for single base alterations by oligonucleotide microarray.

The oligonucleotide microarray, a high-throughput polymorphism detection technology, holds great promise for the characterization of complex genetic variance. To achieve greater sensitivity and specificity for it to be an effective platform technology we present results and discuss some of the factors influencing signal intensities and single-mismatch discrimination in array-based mutation/SNP detection. Probes with a series of concentrations were spotted onto the slide in order to find the optimal concentration with the identifiable satisfying signals and the stable ratios between matched and mismatched probes. It was found that under our experimental conditions, when the initial probe concentration is higher than the maximum immobilization capability of the slide (7.5 microM), the hybridization signal will be saturated and the ratio between matched and mismatched probes will be more stable than at a lower probe concentration. Considering the cost of probes and the systematic stability, a constant spotting concentration of 10 microM was selected. The stability of different types of mismatched oligo-DNA duplexes on the glass surface was also confirmed. The results show that the order of stability of mismatched oligo-DNA duplexes on a glass surface is in general agreement with previous reports conducted using liquid and polyacrylamide gel pads. This suggests that the influence of the mismatched base pair on the stability of the duplex in a solid hybridization system is similar to that in the solution hybridization environment.

[Microarray DNA chip in analyzing the association between HLA-DRB and advanced hepatosplenic schistosomiasis].

OBJECTIVE: To explore possible associations between host polymorphism of HLA class II genotypes and advanced hepatosplenic schistosomiasis japonica. METHODS: 45 advanced schistosomiasis patients (experimental group) and 44 age- and sex-matched patients with chronic schistosomiasis (control group) from the same area were investigated for their HLA class II gene DRB genotypes by genotyping the alleles using microarray DNA chip. The correlation of allele frequencies to advanced hepatosplenic schistosomiasis was compared for the two groups. RESULTS: HLA-DRB1*04x exhibited markedly higher frequency in advanced patients than that in control group (P < 0.01, RR = 3.928). In contrast, the frequency of HLA-DRB1*15x in advanced patients was much lower when compared with that in control group (P < 0.01, RR = 0.050). Besides, the significant allele HLA-DRB1*15x displayed concurrence with allele DRB5*010x/020x. The linked gene haplotype DRB1*15x-DRB5*010x/020x showed significantly higher incidence in control group than in experimental group (P < 0.01). CONCLUSION: Allele HLA-DRB1*04x is positively, while HLA-DRB1*15x is negatively, correlated with advanced hepatosplenic schistosomiasis.

[Detecting genetic polymorphisms of CYP1 A1 and GSTM1 simultaneously with oligonucleotide microarray].

Cytochrome P4501A1 plays a major role in the bioactivation of a number of tobacco procarcinogens. Glutathione S-transferase( GSTM1), a member of the class of GST gene family, has been shown to be polymorphic because of gene deletion resulting in a failure to express the GSTM1 gene in 50% approximately 60% of individuals. Some CYP1 A1/GSTM1 null genotype combinations seem to predispose the lung, esophagus, and oral cavity of smokers to an even higher risk for cancer or DNA damage, requiring, however, confirmation. An easy and reliable oligonuleotide microarray approach validated through direct sequencing method is developed in order to accurately detect single nucleotide polymorphisms of CYP1 A1 gene and discriminate the presence and absence of GSTM1 gene. The m1 (Msp I) and m2 (Ile462Val) polymorphisms of CYP1 A1 gene and GSTM1 null genotype were also determined in a random population of 84 healthy, unrelated volunteers with developed microarray-based method. Of 84 cases, 47.6% were calssified as GSTM1 null, close to the published data. It's interesting that there lack three genotypes of m1 -m2 locus in the population: TT-AG, TT-GG and TC-GG. However, according to the data of the genotype frequencies independently happened at both m1 and m2 site, the combination frequencies of above three genotypes are 11.4%, 2.6%, and 3.1% respectively. Therefore we assume that the haplotypes of m1 -m2 are only T-A, C-A and C-G, but not T-G, as it were,there is no recombination happened between m1 site and m2 site. The frequencies of three haplotypes of T-A, C-A and C-G, calculated through corresponding genotypes, are 69.6%, 7.7% and 22.6% respectively.

Comparison of hybridization behavior between double and single strands of targets and the application of asymmetric PCR targets in cDNA microarray.

Double stranded targets on the cDNA microarray contain representatives of both the coding and noncoding strands, which will introduce hybridization competition with probes. Here, the effect of double and single strands of targets on the signal intensity and the ratios of Cy5/Cy3 within the same slide were compared. The results show that single stranded targets can increase the hybridization efficiency without changing the Cy5/Cy3 ratio. Based on these results, a new strategy was established by generating cDNA targets with asymmetric PCR, instead of conventional PCR, to increase the sensitivity of the cDNA microarray. Furthermore, the feasibility of this approach was validated. The results indicate that the cDNA microarray system based on asymmetric PCR is more sensitive, with no decrease in the reliability and reproducibility as compared with that based on conventional symmetric PCR.

Systematic comparison of the fidelity of aRNA, mRNA and T-RNA on gene expression profiling using cDNA microarray.

In cDNA microarray technology, there are three main reverse transcription based RNA labeling methods, using total RNA (T-RNA), mRNA, and amplified antisense RNA (aRNA), respectively. However, despite the common use of the three types of RNAs, limited data are available regarding their differences and concordances. In this report, we compared the three methods through two sets of self-comparison experiments using the same RNA sample in all cases. Within each method, duplicate hybridizations are highly reproducible with low biases, which are randomly produced. When combining different RNAs within a single array, correlation coefficients between the two channels are rather low, while the discrepancies are persistent. Furthermore, the fidelity of aRNA and mRNA microarrays in the expression profile study shows no significant difference with standard T-RNA based labeling methods. These results suggest that some RNA abundance are selectively changed during aRNA amplification/mRNA purification processes, but it will not affect the gene expression ratio of the two samples if the same type RNA are used. Therefore all three types of RNAs can be used in expression profiling analysis as long as the test and reference samples are generated by identical method within single study.

[Validation of cDNA microarray technology].

cDNA microarray is a technological approach that has the potential to globally measure changes in mRNA expression levels. Self-comparison experiments with the same kind of tissue and differential expression experiments with the different kinds of tissue have been done to verify the reproducibility and the accuracy of this technique. The parameter of the reliability and the reproducibility of the microarray data were analyzed by correlation coefficient (R), coefficient of variation (CV) and false positive rate (FPR) etc. Meanwhile, the error resource also has been inspected. These results showed that generally the correlation coefficient of data from this cDNA microarray system was more than 0.9, the coefficient of variation was about 15%, and the false positive rate was below 3%. The result proves the accuracy of the cDNA microarray data. Consistence rate (CR) was advanced here as a new parameter to evaluate the reproducibility of two replicate experiments. It has some advantages over correlation coefficient and coefficient of variation. The influence of some important factors in the experiments, such as different concentration of spotted DNA, mRNA and total RNA, different batches of slides and different processes of labeling, have been investigated by comparing the results. It was shown that most of the false position produced by the experiment system could be reduced by replicate experiments.

[Optimization of T7-based RNA amplification system for cDNA microarray].

cDNA microarrays are powerful parallel tools for gene expression profiling analysis, which help us to understand the molecular mechanism of diseases and to identify potential targets for therapeutic intervention. However, their broader application are hampered by the large amount of RNA required: up to 200 microg of total RNA or 5 microg of mRNA for one chip, making analysis of small samples difficult. In this work, combined with a template switching effect, the T7 RNA linear amplification procedure was optimized, providing multiple copies of anti-sense RNA with reduced sample inputs: no more than 3 microg total RNA for one chip. Using the same RNA sample in all cases, the anti-sense RNA labeling method was compared with standard total RNA and mRNA methods by two sets of self-comparison experiments. Furthermore, the three methods in profiling analysis were compared with the same pair RNA samples. All the results indicated that the fidelity, reproducibility, and reliability showed no significant difference with conventional total RNA or mRNA microarrays.

Discovery and analysis of hepatocellular carcinoma genes using cDNA microarrays.

PURPOSE: Microarray analysis on a genomic scale was used to profile changes in gene expression accompanying hepatocellular carcinoma. METHODS: Gene expression profiles of liver tissues from twelve hepatocellular carcinoma samples relative to the gene expression profile of the normal liver tissue were analyzed using 4096 chips and 12800 chips. The results of microarray experiments were verified by the Northern blot technique. RESULTS: A group of 1,820 genes with altered expression were identified in more than 50% of the patients examined. This highly concordant expression profile included human genes encoding proteins involved in the function of peroxisomes, serum control, polycyclic aromatic hydrocarbon carcinogenesis, cell growth and differentiation, metastasis, the function of the immune system, apoptosis, and remodeling of the cytoskeleton. CONCLUSIONS: The newly identified genes afford a quantitative view of the changes that accompany liver cancer at the genomic level, enable deeper insights into the molecular basis of disease, and provide an extensive list of potential early-onset molecular markers for improved diagnosis.

A New Strategy for Optimizing the Translation Initiation Rate of Foreign Gene Expression in E. coli.

The translation initiation rate is greatly affected by the secondary structure of the translation initiation region (TIR) of mRNA. A novel system was established for improving the translation initiation rate of a foreign gene in E. coli. As a model, the 5' 114 bp coding sequence (38 amino acids from the start codon) of human proliferating cell nuclear antigen (PCNA) gene was fused with the lacZ' gene at its 5' end in vector pTZ19R. A Shine/Dalgarno (SD) sequence GAGGT was inserted to the -8 position of AUG by site-directed mutation. Then the flanking sequences of SD, which were the 6 nucleotides upstream the SD and the 7 nucleotides between SD and AUG, were randomly changed by PCR using a synthetic primer with partially random sequences. This random mutation led to potential variations in the secondary structure of the TIR of mRNA through base pairing with the 5' coding sequence. The 5' PCNA-LacZ' mRNA could be efficiently and specifically transcribed by inducible T7 RNA polymerase in E. coli strain JM 109 (DE3). There were 269 clones of 5' PCNA-LacZ' fusion plasmid selected first by the blue color on X-gal plate and then by hybridization with a 5' PCNA probe. Eight clones with different blue colors among the 269 clones were chosen for beta-galactosidase activity assay. The results showed that the difference between their enzyme activities was more than 20 fold, but there seemed no apparent difference at the transcription level as assayed by RNA dot hybridization, which suggested that the difference in the expression of the fusion protein was due to the different rate of translation initiation. Thus, by this strategy, an effective translation initiation region for the high expression of human PCNA gene in E. coli was obtained.