Hepatocellular carcinoma (HCC) is the most common type of liver cancer worldwide and one of the deadliest cancers in Asia. But at present, effective targets for HCC clinical therapy are still limited. The “guilt by association” rule suggests that interacting proteins share the same or similar functions and hence may be involved in the same pathway. This assumption can be used to identify disease related genes from protein association networks constructed from existing PPI data. Given the close association between Hepatitis B virus and Hepatitis B which may lead to HCC, here we develop a computational method to identify hepatocellular carcinoma related genes based on k-th shortest paths in the protein–protein interaction (PPI) network (we set k = 1, 2 in this study). Finally, we found 33 genes whose p-values were less than 0.05, and most of them have been reported to be involved in HCC tumorigenesis and development. The results also provide a new reference for research into HCC oncogenesis and for development of new strategies for HCC clinical therapies.

A single mammalian transcript normally encodes one protein, but the transcript of GNAS (G-protein α-subunit) contains two reading frames and produces two structurally unrelated proteins, XLαs and ALEX. No other confirmed GNAS-like dual-coding transcripts have been reported to date, even though many such candidate genes have been predicted by bioinformatics analysis. In this study, we constructed a series of vectors to test how two protein products were translated from a single transcript in vitro. The length of the ORF (open reading frame), position of the first AUG and the Kozak motif were found to be important factors. These factors, as well as 55-bp NMD (nonsense-mediated mRNA decay) rule, were used in a bioinformatics search for candidate dual-coding transcripts. A total of 1307, 750 and 474 two-ORF-containing transcripts were found in human, mouse and rat, respectively, of which 170, 89 and 70, respectively, were found to be potential dual-coding transcripts. Most transcripts showed low conservation among species. Interestingly, dual-coding transcripts were significantly enriched for transcripts from the zinc-finger protein family, which are usually DNA-binding proteins involved in regulation of the transcription process.

There is increasing realization that regulation of genes is done partly at the RNA level by sense-antisense binding. Studies typically concentrate on the role of non-coding RNAs in regulating coding RNA. But the majority of transcripts in a cell are likely to be coding. Is it possible that coding RNA might regulate other coding RNA by short perfect sense-antisense binding? Here we compare all well-described human protein coding mRNAs against all others to identify sites 15-25 bp long that could potentially perfectly match sense-antisense.From 24,968 protein coding mRNA RefSeq sequences, none failed to find at least one match in the transcriptome. By randomizations generating artificial transcripts matched for G+C content and length, we found that there are more such trans short sense-antisense pairs than expected. Several further features are consistent with functionality of some of the putative matches. First, transcripts with more potential partners have lower expression levels, and the pair density of tissue specific genes is significantly higher than that of housekeeping genes. Further, the single nucleotide polymorphism density is lower in short pairing regions than it is in flanking regions. We found no evidence that the sense-antisense pairing regions are associated with small RNAs derived from the protein coding genes.Our results are consistent with the possibility of common short perfect sense-antisense pairing between transcripts of protein coding genes.

Human isolated gingival fibromatosis is an oral disorder characterized by a slowly progressive benign enlargement of gingival tissues. The most common genetic form, hereditary gingival fibromatosis (HGF), is usually transmitted as an autosomal dominant trait. We report here for the first time a newly identified maternally inherited gingival fibromatosis in two unrelated Chinese families and mapped this disease locus to human chromosome 11p15 with a maximum two point LOD score of 8.70 at D11S4046 (θ = 0) for family 1 and of 6.02 at D11S1318 for family 2. Haplotype analysis placed the critical region in the interval defined by D11S1984 and D11S1338. A cluster of maternally expressed genes is within this critical region. We screened individuals in these two families for mutations for all known maternally expressed genes within this region. None was found either within the coding sequence or at the intron–exon boundary of these genes. Neither did we detect any loss of imprinting in three informative imprinted genes including H19, KCNQ1 downstream neighbor (KCNQ1DN) and cyclin-dependent kinase inhibitor 1C (CDKN1C). However, gene expression profile analysis revealed reduced expression of hemoglobin beta (HBB), hemoglobin delta (HBD), hemoglobin gamma A (HBG1) and hemoglobin gamma G (HBG2) genes at disease locus in HGF patients. This study suggests that genome imprinting might affect the development of HGF.

Ichthyosis vulgaris (IV) is an inherited scaling skin disorder with a prevalence estimated at 2.29% in China. The gene responsible for this disorder has not been elucidated. To find the disease gene, we ascertained two Chinese IV families. Linkage analysis identified an IV locus on chromosome 1q22 with a maximum two-point Lod score of 2.47 at D1S1653 (=0.00). Haplotype analysis placed the critical region in a 7-cM interval defined by D1S1653 and D1S2675. These results provide the basis for further identifying the gene responsible for IV disorder.

MicroRNA (miRNA)-mediated RNA silencing and nonsense-mediated decay (NMD) are two conserved RNA-level regulatory pathways. Although they are mechanically different, both can regulate target genes by RNA degradation and translational repression. Moreover, studies of individual target genes indicated that these two pathways can be involved in the same processes (e.g., development and stress responses). These facts raise an important question that whether these two systems are cooperative, interchangeable or optimal for regulation of different sorts of genes. We addressed this by comparing miRNA and NMD targets in Arabidopsis thaliana at the genome-wide scale. We find no more overlap in the genes targeted by both systems than expected by chance. Moreover, the sorts of genes or pathways regulated by these systems are categorically different on several cross-correlating fronts. While miRNA targets show enrichment in the process of development, metabolism and transcription, NMD targets are associated with stress responses but otherwise poorly annotated. Validated miRNA targets are more highly expressed, less variably expressed and slower evolving. These differences suggest that the modes of regulation need not be interchangeable. Instead, we suggest that miRNA genes are commonly dose-sensitive and require fine control of levels through weak pull-down by miRNAs. This is consistent with miRNA-regulated genes being more likely to be involved in protein–protein interactions. Many NMD-regulated genes, by contrast, have properties consistent with them being rapid emergency response “fire-fighter” genes. If true, the lack of annotation of NMD targets suggests that we poorly understand the emergencies plants face in the wild.

Although there is an accumulating appreciation of the key roles that long intergenic non-coding RNAs (lincRNAs) play in diverse cellular processes, our knowledge of how lincRNAs function in cancer remains sparse. Here, we present a comprehensive landscape of RNA-seq transcriptome profiles of lung adenocarcinomas and their paired normal counterparts to unravel gene regulation rules of lincRNAs. Consistent with previous findings of co-expression between neighboring protein-coding genes, lincRNAs were typically co-expressed with their neighboring genes, which was found in both cancerous and normal tissues. By building a mathematical model based on correlated gene expression, we distinguished an additional subset of lincRNAs termed “regulatory lincRNAs”, representing their dominant roles in gene regulation. The number of regulatory lincRNAs was significantly higher in cancerous compared to normal tissues, and most of them positively regulated protein-coding genes in trans. Functional validation, using knockdown, determined that regulatory lincRNA, GAS5, affected its predicted protein-coding targets. Moreover, we discovered hundreds of differentially expressed regulatory lincRNAs with inclusion of some cancer-associated lincRNAs. Our integrated analysis reveals enhanced regulatory effects of lincRNAs and provides a resource for the study of regulatory lincRNAs that play critical roles in lung adenocarcinoma.