The intron has been a big biological mystery since it was first discovered in several aspects. First, all of the completely sequenced eukaryotes harbor introns in the genomic structure, whereas no prokaryotes identified so far carry introns. Second, the amount of total introns varies in different species. Third, the length and number of introns vary in different genes, even within the same species genome. Fourth, all introns are copied into RNAs by transcription and DNAs by replication processes, but intron sequences do not participate in protein-coding sequences. The existence of introns in the genome should be a burden to some cells, because cells have to consume a great deal of energy to copy and excise them exactly at the correct positions with the help of complicated spliceosomal machineries. The existence throughout the long evolutionary history is explained, only if selective advantages of carrying introns are assumed to be given to cells to overcome the negative effect of introns. In that regard, we summarize previous research about the functional roles or benefits of introns. Additionally, several other studies strongly suggesting that introns should not be junk will be introduced.

Genomes of higher eukaryotes have surprisingly long first introns and in some cases, the first introns have been shown to have higher conservation relative to other introns. However, the functional relevance of conserved regions in the first introns is poorly understood. Leveraging the recent ENCODE data, here we assess potential regulatory roles of conserved regions in the first intron of human genes. We first show that relative to other downstream introns, the first introns are enriched for blocks of highly conserved sequences. We also found that the first introns are enriched for several chromatin marks indicative of active regulatory regions and this enrichment of regulatory marks is correlated with enrichment of conserved blocks in the first intron; the enrichments of conservation and regulatory marks in first intron are not entirely explained by a general, albeit variable, bias for certain marks toward the 5’ end of introns. Interestingly, conservation as well as proportions of active regulatory chromatin marks in the first intron of a gene correlates positively with the numbers of exons in the gene but the correlation is significantly weakened in second introns and negligible beyond the second intron. The first intron conservation is also positively correlated with the gene’s expression level in several human tissues. Finally, a gene-wise analysis shows significant enrichments of active chromatin marks in conserved regions of first introns, relative to the conserved regions in other introns of the same gene. Taken together, our analyses strongly suggest that first introns are enriched for active transcriptional regulatory signals under purifying selection.

Single-cell sequencing, which is used to detect clinically important tumor subpopulations, is necessary for understanding tumor heterogeneity. Here, we analyzed transcriptomic data obtained from 34 single cells from human lung adenocarcinoma (LADC) patient-derived xenografts (PDXs). To focus on the intrinsic transcriptomic signatures of these tumors, we filtered out genes that displayed extensive expression changes following xenografting and cell culture. Then, we performed clustering analysis using co-regulated gene modules rather than individual genes to minimize read drop-out errors associated with single-cell sequencing. This combined approach revealed two distinct intra-tumoral subgroups that were primarily distinguished by the gene module G64. The G64 module was predominantly composed of cell-cycle genes. E2F1 was found to be the transcription factor that most likely mediates the expression of the G64 module in single LADC cells. Interestingly, the G64 module also indicated inter-tumoral heterogeneity based on its association with patient survival and other clinical variables such as smoking status and tumor stage. Taken together, these results demonstrate the feasibility of single-cell RNA sequencing and the strength of our analytical pipeline for the identification of tumor subpopulations.

One of the main objectives of the molecular evolution and evolutionary systems biology field is to reveal the underlying principles that dictate protein evolutionary rates. Several studies argue that expression abundance is the most critical component in determining the rate of evolution, especially in unicellular organisms. However, the expression breadth also needs to be considered for multicellular organisms.

Dilated cardiomyopathy (DCM) is one of the main causes of heart failure (called cardiomyopathies) in adults. Alterations in epigenetic regulation (i.e., DNA methylation) have been implicated in the development of DCM. Here, we identified a total of 1828 differentially methylated probes (DMPs) using the Infinium 450K HumanMethylation Bead chip by comparing the methylomes between 18 left ventricles and 9 right ventricles. Alterations in DNA methylation levels were observed mainly in lowly methylated regions corresponding to promoter-proximal regions, which become hypermethylated in severely affected left ventricles. Subsequent mRNA microarray analysis showed that the effect of DNA methylation on gene expression regulation is not unidirectional but is controlled by the functional sub-network context. DMPs were significantly enriched in the transcription factor binding sites (TFBSs) we tested. Alterations in DNA methylation were specifically enriched in the cis-regulatory regions of cardiac development genes, the majority of which are involved in ventricular development (e.g., TBX5 and HAND1).