Search

Abstract Understanding transcriptional responses to chemical perturbations is central to drug discovery, but exhaustive experimental screening of disease-compound combinations is unfeasible. To overcome this limitation, here we introduce PRnet, a perturbation-conditioned deep generative model that predicts transcriptional responses to novel chemical perturbations that have never experimentally perturbed at bulk and single-cell levels. Evaluations indicate that PRnet outperforms alternative methods in predicting responses across novel compounds, pathways, and cell lines. PRnet enables gene-level response interpretation and in-silico drug screening for diseases based on gene signatures. PRnet further identifies and experimentally validates novel compound candidates against small cell lung cancer and colorectal cancer. Lastly, PRnet generates a large-scale integration atlas of perturbation profiles, covering 88 cell lines, 52 tissues, and various compound libraries. PRnet provides a robust and scalable candidate recommendation workflow and successfully recommends drug candidates for 233 diseases. Overall, PRnet is an effective and valuable tool for gene-based therapeutics screening.

Abstract Glioblastoma (GBM) ranks among the most lethal of human cancers, containing glioma stem cells (GSCs) that display therapeutic resistance. Here, we report that the lncRNA INHEG is highly expressed in GSCs compared to differentiated glioma cells (DGCs) and promotes GSC self-renewal and tumorigenicity through control of rRNA 2’-O-methylation. INHEG induces the interaction between SUMO2 E3 ligase TAF15 and NOP58, a core component of snoRNP that guides rRNA methylation, to regulate NOP58 sumoylation and accelerate the C/D box snoRNP assembly. INHEG activation enhances rRNA 2’-O-methylation, thereby increasing the expression of oncogenic proteins including EGFR, IGF1R, CDK6 and PDGFRB in glioma cells. Taken together, this study identifies a lncRNA that connects snoRNP-guided rRNA 2’-O-methylation to upregulated protein translation in GSCs, supporting an axis for potential therapeutic targeting of gliomas.

Abstract AlphaFold2 (AF2) is an artificial intelligence (AI) system developed by DeepMind that can predict three-dimensional (3D) structures of proteins from amino acid sequences with atomic-level accuracy. Protein structure prediction is one of the most challenging problems in computational biology and chemistry, and has puzzled scientists for 50 years. The advent of AF2 presents an unprecedented progress in protein structure prediction and has attracted much attention. Subsequent release of structures of more than 200 million proteins predicted by AF2 further aroused great enthusiasm in the science community, especially in the fields of biology and medicine. AF2 is thought to have a significant impact on structural biology and research areas that need protein structure information, such as drug discovery, protein design, prediction of protein function, et al. Though the time is not long since AF2 was developed, there are already quite a few application studies of AF2 in the fields of biology and medicine, with many of them having preliminarily proved the potential of AF2. To better understand AF2 and promote its applications, we will in this article summarize the principle and system architecture of AF2 as well as the recipe of its success, and particularly focus on reviewing its applications in the fields of biology and medicine. Limitations of current AF2 prediction will also be discussed.

Recent advances in RNA engineering have enabled the development of RNA-based therapeutics for a broad spectrum of applications. Developing RNA therapeutics start with targeted RNA screening and move to the drug design and optimization. However, existing target screening tools ignore noncoding RNAs and their disease-relevant regulatory relationships. And designing therapeutic RNAs encounters high computational complexity of multi-objective optimization to overcome the immunogenicity, instability and inefficient translational production. To unlock the therapeutic potential of noncoding RNAs and enable one-stop screening and design of therapeutic RNAs, we have built the platform TREAT. It incorporates 43,087,953 regulatory relationships between coding and noncoding genes from 81 biological networks under different physiological conditions. TREAT introduces graph representation learning with Random Walk Diffusions to perform disease-relevant target screening, in addition to the commonly utilized Topological Degree and PageRank algorithms. Design and optimization of large RNAs or interfering RNAs are both available. To reduce the computational complexity of multi-objective optimization for large RNA, we stratified the features into local and global features. The local features are evaluated on the fixed-length or dynamic-length local bins, whereas the latter are inspired by AI language models of protein sequence. Then the global assessment is performed on refined candidates, thus reducing the enormous search space. Overall, TREAT is a one-stop platform for the screening and designing of therapeutic RNAs, with particular attention to noncoding RNAs and cutting-edge AI technology embedded, leading the progress of innovative therapeutics for challenging diseases. TREAT is freely accessible at https://rna.org.cn/treat.

IntroductionCentral nervous system (CNS) diseases, such as neurodegenerative disorders and brain diseases caused by acute injuries, are important, yet challenging to study due to disease lesion locations and other complexities.MethodsUtilizing the powerful method of spatial transcriptome analysis together with novel algorithms we developed for the study, we report here for the first time a 3D trajectory map of gene expression changes in the brain following acute neural injury using a mouse model of intraventricular hemorrhage (IVH). IVH is a common and representative complication after various acute brain injuries with severe mortality and mobility implications.ResultsOur data identified three main 3D global pseudospace-time trajectory bundles that represent the main neural circuits from the lateral ventricle to the hippocampus and primary cortex affected by experimental IVH stimulation. Further analysis indicated a rapid response in the primary cortex, as well as a direct and integrated effect on the hippocampus after IVH stimulation.DiscussionThese results are informative for understanding the pathophysiological changes, including the spatial and temporal patterns of gene expression changes, in IVH patients after acute brain injury, strategizing more effective clinical management regimens, and developing novel bioinformatics strategies for the study of other CNS diseases. The algorithm strategies used in this study are searchable via a web service (www.combio-lezhang.online/3dstivh/home).

Currently, oil-coated PVA fibers are the most commonly used material in ECC research. However, the high price limits the application of PVA-ECC in practical engineering. In order to reduce the cost, one of the methods is to partially replace the PVA fibers in ECC. In order to demonstrate the feasibility of PVA/BF-ECC and PVA/PP-ECC, polyvinyl alcohol fibers (PVA), basalt fibers (BFs) and polypropylene fibers (PP) were added at 0.5%, 1.0% and 1.5% by volume of PVA in addition to 1% by volume of PVA. Subsequently, tensile, compression and drop-weight impact tests were conducted on single or hybrid fiber concrete. The results showed that the post-peak compression toughness, tensile strength, and initial cracking impact strength of PVA/BF-ECC and PVA/PP-ECC increased significantly with the increase in the volume ratio of BF and PP fibers, while the performance of PVA-ECC materials with the same fiber volume ratio decreased slightly. Therefore, the cost can be reduced by designing hybrid PVA/BF-ECC materials that meet the performance requirements. The experimental evidence presented in this study demonstrates the feasibility and reasonable prospect of the new hybrid PVA/BF-ECC.

Abstract Background A specific 3-dimensional intrachromosomal architecture of core stem cell factor genes is required to reprogram a somatic cell into pluripotency. As little is known about the epigenetic readers that orchestrate this architectural remodeling, we used a novel chromatin RNA in situ reverse transcription sequencing (CRIST-seq) approach to profile long noncoding RNAs (lncRNAs) in the Oct4 promoter. Results We identify Platr10 as an Oct4 - Sox2 binding lncRNA that is activated in somatic cell reprogramming. Platr10 is essential for the maintenance of pluripotency, and lack of this lncRNA causes stem cells to exit from pluripotency. In fibroblasts, ectopically expressed Platr10 functions in trans to activate core stem cell factor genes and enhance pluripotent reprogramming. Using RNA reverse transcription-associated trap sequencing (RAT-seq), we show that Platr10 interacts with multiple pluripotency-associated genes, including Oct4, Sox2, Klf4, and c-Myc, which have been extensively used to reprogram somatic cells. Mechanistically, we demonstrate that Platr10 helps orchestrate intrachromosomal promoter-enhancer looping and recruits TET1, the enzyme that actively induces DNA demethylation for the initiation of pluripotency. We further show that Platr10 contains an Oct4 binding element that interacts with the Oct4 promoter and a TET1-binding element that recruits TET1. Mutation of either of these two elements abolishes Platr10 activity. Conclusion These data suggest that Platr10 functions as a novel chromatin RNA molecule to control pluripotency in trans by modulating chromatin architecture and regulating DNA methylation in the core stem cell factor network.

The redox homeostasis system regulates many biological processes, intracellular antioxidant production and redox signaling. However, long noncoding RNAs (lncRNAs) involved in redox regulation have rarely been reported. Herein, we reported that downregulation of MAGI2-AS3 decreased the superoxide level in Human fibroblasts (Fbs), a replicative aging model, as detected by the fluorescent probes dihydroethidium (DHE) and MitoSOX™ Red. RNA pulldown combined with mass spectrometry showed that HSPA8 is a novel interacting protein of MAGI2-AS3, which was further confirmed by photoactivatable ribonucleoside–enhanced crosslinking and immunoprecipitation (PAR-CLIP). Downregulation of MAGI2-AS3 decreased the hydrogen peroxide (H2O2) content by stabilizing the HSPA8 protein level via inhibiting the protesome degradation of HSPA8. Further evidence showed that MAGI2-AS3 interacted with the C-terminal domain (CTD) of HSPA8. Downregulation of MAGI2-AS3 delayed cell senescence, while this antiaging effect was abolished by HSPA8 knockdown. The underlying molecular mechanism by which MAGI2-AS3 knockdown inhibited cell senescence was mediated via suppression of the ROS/MAP2K6/p38 signaling pathway. Taken together, these findings revealed that downregulation of lncRNA MAGI2-AS3 decreased the H2O2 content and delayed cell senescence by stabilizing the HSPA8 protein level, identifying a potential antiaging application.

BCL9 is an activator of oncogenic Wnt/β-catenin. Here, the authors show a β-catenin independent function of BCL9 in a subtype of colorectal cancers, where it interacts with paraspeckle proteins to enhance the mRNA stability of calcium signalling and neural associated genes to promote communication with tumour cells and its microenvironment.

Abstract Advancements in omics‐based technologies over the past few years have led to the discovery of numerous biologically relevant peptides encoded by small open reading frames (smORFs) embedded in long noncoding RNA (lncRNA) transcripts (referred to as microproteins here) in a variety of species. However, the mechanisms and modes of action that underlie the roles of microproteins have yet to be fully characterized. Herein, we provide the first experimental evidence of abundant microproteins in extracellular vesicles (EVs) derived from glioma cancer cells, indicating that the EV‐mediated transfer of microproteins may represent a novel mechanism for intercellular communication. Intriguingly, when examining human plasma, 48, 11 and 3 microproteins were identified from purified EVs, whole plasma and EV‐free plasma, respectively, suggesting that circulating microproteins are primarily enriched in EVs. Most importantly, the preliminary data showed that the expression profile of EV microproteins in glioma patient diverged from the health donors, suggesting that the circulating microproteins in EVs might have potential diagnostic application in identifying patients with glioma.

Non alcoholic fatty liver disease (NAFLD) associates with an elevated risk of developing hepatocellular carcinoma (HCC). Here, the authors find that Nogo-B, an endoplasmic reticulum resident protein, is upregulated by lipid uptake and acts as an oncogene in NAFLD-associated HCC by promoting lipid droplet breakdown by lipophagy and triggering Hippo pathway dysregulation

Abstract Background Increasing evidence has demonstrated that mesenchymal stem cells (MSCs) play a role in the construction of tumor microenvironments. Co-culture between tumor cells and MSCs provides an easy and useful platform for mimicking tumor microenvironments and identifying the important members involved in tumor progress. The long non-coding RNAs (lncRNAs) have been shown to regulate different tumorigenic processes. In this study, we aimed to examine functional lncRNA deregulations associated with breast cancer malignancy instigated by MSC-MCF-7 co-culture. Methods The microarrays were used to profile the expression changes of lncRNAs in MCF-7 cells during epithelial-mesenchymal transition (EMT) induced by co-culture with MSCs. We found that an intergenic lncRNA KB-1732A1.1 (termed LincK, partly overlapped with GASL1) was significantly elevated. To investigate the biological function of LincK, the expression of EMT markers, cell migration, invasion, proliferation, and colony formation were evaluated in vitro and xenograft assay in nude mice were performed in vivo. Furthermore, we detected LincK expression in clinical samples using RNAscope® technology and verified aberrant expression of LincK in breast cancer data sets from The Cancer Genome Atlas (TCGA) by bioinformatic analysis. The underlying mechanisms of LincK were investigated using mRNA microarray analyses, Western blot, RNA pull down, and RNA immunoprecipitation. Results LincK induced an EMT progress in breast cancer cells (BCC) MCF-7, MDA-MB-453, and MDA-MB-231. The depletion of LincK decreased the growth, migration, and invasion in BCC, whereas the overexpression of LincK exerted the opposite effects. Moreover, knockdown of LincK repressed tumorigenesis, and ectopic expression of LincK promoted tumor growth in MCF-7 xenograft model. LincK ablation in MDA-MB-231 cells dramatically impaired lung metastasis when incubated intravenously into nude mice. Further, LincK was frequently elevated in breast cancer compared with normal breast tissue in clinical samples. Mechanistically, LincK may share common miRNA response elements with PBK and ZEB1 and regulate the effects of miR-200 s. Conclusion LincK plays a significant role in regulating EMT and tumor growth and could be a potential therapeutic target in breast cancer.

Abstract Background The majority of eukaryote genomes can be actively transcribed into non-coding RNAs (ncRNAs), which are functionally important in development and evolution. In the study of maize, an important crop for both humans and animals, aside from microRNAs and long non-coding RNAs, few studies have been conducted on intermediate-size ncRNAs. Results We constructed a homogenized cDNA library of 50–500 nt RNAs in the maize inbred line Chang 7–2. Sequencing revealed 169 ncRNAs, which contained 58 known and 111 novel ncRNAs (including 70 snoRNAs, 27 snRNAs, 13 unclassified ncRNAs and one tRNA). Forty of the novel ncRNAs were specific to the Panicoideae, and 24% of them are located on sense-strand of the 5′ or 3′ terminus of protein coding genes on chromosome. Target site analysis found that 22 snoRNAs can guide to 38 2’-O-methylation and pseudouridylation modification sites of ribosomal RNAs and small nuclear RNAs. Expression analysis showed that 43 ncRNAs exhibited significantly altered expression in different tissues or developmental stages of maize seedlings, eight ncRNAs had tissue-specific expression and five ncRNAs were strictly accumulated in the early stage of leaf development. Further analysis showed that 3 of the 5 stage-specific ncRNAs (Zm-3, Zm-18, and Zm-73) can be highly induced under drought and salt stress, while one snoRNA Zm-8 can be repressed under PEG-simulated drought condition. Conclusions We provided a genome-wide identification and functional analysis of ncRNAs with a size range of 50–500 nt in maize. 111 novel ncRNAs were cloned and 40 ncRNAs were determined to be specific to Panicoideae. 43 ncRNAs changed significantly during maize development, three ncRNAs can be strongly induced under drought and salt stress, suggesting their roles in maize stress response. This work set a foundation for further study of intermediate-size ncRNAs in maize.

Abstract Aberrant regulation of miRNA genes contributes to pathogenesis of a wide range of human diseases, including cancer. The TAR DNA binding protein 43 (TDP-43), a RNA/DNA binding protein associated with neurodegeneration, is involved in miRNA biogenesis. Here, we systematically examined miRNAs regulated by TDP-43 using RNA-Seq coupled with an siRNA-mediated knockdown approach. TDP-43 knockdown affected the expression of a number of miRNAs. In addition, TDP-43 down-regulation led to alterations in the patterns of different isoforms of miRNAs (isomiRs) and miRNA arm selection, suggesting a previously unknown role of TDP-43 in miRNA processing. A number of TDP-43 associated miRNAs, and their candidate target genes, are associated with human cancers. Our data reveal highly complex roles of TDP-43 in regulating different miRNAs and their target genes. Our results suggest that TDP-43 may promote migration of lung cancer cells by regulating miR-423-3p. In contrast, TDP-43 increases miR-500a-3p expression and binds to the mature miR-500a-3p sequence. Reduced expression of miR-500a-3p is associated with poor survival of lung cancer patients, suggesting that TDP-43 may have a suppressive role in cancer by regulating miR-500a-3p. Cancer-associated genes LIF and PAPPA are possible targets of miR-500a-3p. Our work suggests that TDP-43-regulated miRNAs may play multifaceted roles in the pathogenesis of cancer.

The proboscis extension reflex (PER) is a common classical conditioned reflex which is widely used in the neurology and ethology. In honeybees, PER experiments can train bees to associate an odor with a reward or punishment. Here we present a variation of the PER experiment in Apis mellifera that trains honeybees using only one antenna. This variation on the PER paradigm could assist research efforts in fields which study lateralization within the nervous system.

C/EBPα is a critical transcriptional regulator of adipogenesis. How C/EBPα transcription is itself regulated is poorly understood, however, and remains a key question that needs to be addressed for a complete understanding of adipogenic development. Here, we identify a lncRNA, ADINR (adipogenic differentiation induced noncoding RNA), transcribed from a position ∼450 bp upstream of the C/EBPα gene, that orchestrates C/EBPα transcription in vivo. Depletion of ADINR leads to a severe adipogenic defect that is rescued by overexpression of C/EBPα. Moreover, we reveal that ADINR RNA specifically binds to PA1 and recruits MLL3/4 histone methyl-transferase complexes so as to increase H3K4me3 and decrease H3K27me3 histone modification in the C/EBPα locus during adipogenesis. These results show that ADINR plays important roles in regulating the differentiation of human mesenchymal stem cells into adipocytes by modulating C/EBPα in cis.

Intermediate-size noncoding RNAs (is-ncRNAs) have been shown to play important regulatory roles in the development of several eukaryotic organisms. However, they have not been thoroughly explored in Plasmodium falciparum, which is the most virulent malaria parasite infecting human being. By using Illumina/Solexa paired-end sequencing of an is-ncRNA-specific library, we performed a systematic identification of novel is-ncRNAs in intraerythrocytic P. falciparum, strain 3D7. A total of 1,198 novel is-ncRNA candidates, including antisense, intergenic, and intronic is-ncRNAs, were identified. Bioinformatics analyses showed that the intergenic is-ncRNAs were the least conserved among different Plasmodium species, and antisense is-ncRNAs were more conserved than their sense counterparts. Twenty-two novel snoRNAs were identified, and eight potential novel classes of P. falciparum is-ncRNAs were revealed by clustering analysis. The expression of randomly selected novel is-ncRNAs was confirmed by RT-PCR and northern blotting assays. An obvious different expressional profile of the novel is-ncRNA between the early and late intraerythrocytic developmental stages of the parasite was observed. The expression levels of the antisense RNAs correlated with those of their cis-encoded sense RNA counterparts, suggesting that these is-ncRNAs are involved in the regulation of gene expression of the parasite. In conclusion, we accomplished a deep profiling analysis of novel is-ncRNAs in P. falciparum, analysed the conservation and structural features of these novel is-ncRNAs, and revealed their differential expression patterns during the development of the parasite. These findings provide important information for further functional characterisation of novel is-ncRNAs during the development of P. falciparum.

Social caste determination in the honey bee is assumed to be determined by the dietary status of the young larvae and translated into physiological and epigenetic changes through nutrient-sensing pathways. We have employed Illumina/Solexa sequencing to examine the small RNA content in the bee larval food, and show that worker jelly is enriched in miRNA complexity and abundance relative to royal jelly. The miRNA levels in worker jelly were 7-215 fold higher than in royal jelly, and both jellies showed dynamic changes in miRNA content during the 4(th) to 6(th) day of larval development. Adding specific miRNAs to royal jelly elicited significant changes in queen larval mRNA expression and morphological characters of the emerging adult queen bee. We propose that miRNAs in the nurse bee secretions constitute an additional element in the regulatory control of caste determination in the honey bee.

Introduction. Small noncoding RNAs have important regulatory functions in different cell pathways. It is believed that most of them mainly play role in gene post-transcriptional regulation in the cytoplasm. Recent evidence suggests miRNA and siRNA activity in the nucleus. Here, we show distinct genome-wide sub-cellular localization distribution profiles of small noncoding RNAs in human breast cancer cells. Methods. We separated breast cancer cell nuclei from cytoplasm, and identified small RNA sequences using a high-throughput sequencing platform. To determine the relationship between miRNA sub-cellular distribution and cancer progression, we used microarray analysis to examine the miRNA expression levels in nucleus and cytoplasm of three human cell lines, one normal breast cell line and two breast cancer cell lines. Logistic regression and SVM were used for further analysis. Results. The sub-cellular distribution of small noncoding RNAs shows that numerous miRNAs and their isoforms (isomiR) not only locate to the cytoplasm but also appeare in the nucleus. Subsequent microarray analyses indicated that the miRNA nuclear-cytoplasmic-ratio is a significant characteristic of different cancer cell lines. Conclusions. Our results indicate that the sub-cellular distribution is important for miRNA function, and that the characterization of the small RNAs sub-cellular localizome may contribute to cancer research and diagnosis.

BACKGROUND: A network of DNA damage response (DDR) mechanisms functions coordinately to maintain genome integrity and prevent disease. The Nucleotide Excision Repair (NER) pathway is known to function in the response to UV-induced DNA damage. Although numbers of coding genes and miRNAs have been identified and reported to participate in UV-induced DNA damage response (UV-DDR), the precise role of non-coding RNAs (ncRNAs) in UV-DDR remains largely unknown. METHODOLOGY/PRINCIPAL FINDINGS: We used high-throughput RNA-sequencing (RNA-Seq) to discover intermediate-size (70-500 nt) ncRNAs (is-ncRNAs) in C. elegans, using the strains of L4 larvae of wild-type (N2), UV-irradiated (N2/UV100) and NER-deficient mutant (xpa-1), and 450 novel non-coding transcripts were initially identified. A customized microarray assay was then applied to examine the expression profiles of both novel transcripts and known is-ncRNAs, and 57 UV-DDR-related is-ncRNA candidates showed expression variations at different levels between UV irradiated strains and non- irradiated strains. The top ranked is-ncRNA candidates with expression differences were further validated by qRT-PCR analysis, of them, 8 novel is-ncRNAs were significantly up-regulated after UV irradiation. Knockdown of two novel is-ncRNAs, ncRNA317 and ncRNA415, by RNA interference, resulted in higher UV sensitivity and significantly decreased expression of NER-related genes in C. elegans. CONCLUSIONS/SIGNIFICANCE: The discovery of above two novel is-ncRNAs in this study indicated the functional roles of is-ncRNAs in the regulation of UV-DDR network, and aided our understanding of the significance of ncRNA involvement in the UV-induced DNA damage response.

BACKGROUND: Upwards of 1200 miRNA loci have hitherto been annotated in the human genome. The specific features defining a miRNA precursor and deciding its recognition and subsequent processing are not yet exhaustively described and miRNA loci can thus not be computationally identified with sufficient confidence. RESULTS: We rendered pre-miRNA and non-pre-miRNA hairpins as strings of integrated sequence-structure information, and used the software Teiresias to identify sequence-structure motifs (ss-motifs) of variable length in these data sets. Using only ss-motifs as features in a Support Vector Machine (SVM) algorithm for pre-miRNA identification achieved 99.2% specificity and 97.6% sensitivity on a human test data set, which is comparable to previously published algorithms employing combinations of sequence-structure and additional features. Further analysis of the ss-motif information contents revealed strongly significant deviations from those of the respective training sets, revealing important potential clues as to how the sequence and structural information of RNA hairpins are utilized by the miRNA processing apparatus. CONCLUSION: Integrated sequence-structure motifs of variable length apparently capture nearly all information required to distinguish miRNA precursors from other stem-loop structures.

Somatic cells can be reprogrammed to a pluripotent state by over-expression of defined factors, and pluripotency has been confirmed by the tetraploid complementation assay. However, especially in human cells, estimating the quality of Induced Pluripotent Stem Cell(iPSC) is still difficult. Here, we present a novel supervised method for the assessment of the quality of iPSCs by estimating the gene expression profile using a 2-D "Differentiation-index coordinate", which consists of two "developing lines" that reflects the directions of ES cell differentiation and the changes of cell states during differentiation. By applying a novel liner model to describe the differentiation trajectory, we transformed the ES cell differentiation time-course expression profiles to linear "developing lines"; and use these lines to construct the 2-D "Differentiation-index coordinate" of mouse and human. We compared the published gene expression profiles of iPSCs, ESCs and fibroblasts in mouse and human "Differentiation-index coordinate". Moreover, we defined the Distance index to indicate the qualities of iPS cells, which based on the projection distance of iPSCs-ESCs and iPSCs-fibroblasts. The results indicated that the "Differentiation-index coordinate" can distinguish differentiation states of the different cells types. Furthermore, by applying this method to the analysis of expression profiles in the tetraploid complementation assay, we showed that the Distance index which reflected spatial distributions correlated the pluripotency of iPSCs. We also analyzed the significantly changed gene sets of "developing lines". The results suggest that the method presented here is not only suitable for the estimation of the quality of iPS cells based on expression profiles, but also is a new approach to analyze time-resolved experimental data.

Housekeeping genes (HKGs) generally have fundamental functions in basic biochemical processes in organisms, and usually have relatively steady expression levels across various tissues. They play an important role in the normalization of microarray technology. Using Fourier analysis we transformed gene expression time-series from a Hela cell cycle gene expression dataset into Fourier spectra, and designed an effective computational method for discriminating between HKGs and non-HKGs using the support vector machine (SVM) supervised learning algorithm which can extract significant features of the spectra, providing a basis for identifying specific gene expression patterns. Using our method we identified 510 human HKGs, and then validated them by comparison with two independent sets of tissue expression profiles. Results showed that our predicted HKG set is more reliable than three previously identified sets of HKGs.

The involvement of noncoding RNAs (ncRNAs) in the development of the human brain remains largely unknown. Applying a cloning strategy for detection of intermediate size (50-500 nt) ncRNAs (is-ncRNAs) we have identified 82 novel transcripts in human fetal brain tissue. Most of the novel is-ncRNAs are not well conserved in vertebrates, and several transcripts were only found in primates. Northern blot and microarray analysis indicated considerable variation in expression across human fetal brain development stages and fetal tissues for both novel and known is-ncRNAs. Expression of several of the novel is-ncRNAs was conspicuously absent in one or two brain cancer cell lines, and transient overexpression of some transcripts in cancer cells significantly inhibited cell proliferation. Overall, our results suggest that is-ncRNAs play important roles in the development and tumorigenesis of human brain.

Genome tiling microarray studies have consistently documented rich transcriptional activity beyond the annotated genes. However, systematic characterization and transcriptional profiling of the putative novel transcripts on the genome scale are still lacking. We report here the identification of 25,352 and 27,744 transcriptionally active regions (TARs) not encoded by annotated exons in the rice (Oryza. sativa) subspecies japonica and indica, respectively. The non-exonic TARs account for approximately two thirds of the total TARs detected by tiling arrays and represent transcripts likely conserved between japonica and indica. Transcription of 21,018 (83%) japonica non-exonic TARs was verified through expression profiling in 10 tissue types using a re-array in which annotated genes and TARs were each represented by five independent probes. Subsequent analyses indicate that about 80% of the japonica TARs that were not assigned to annotated exons can be assigned to various putatively functional or structural elements of the rice genome, including splice variants, uncharacterized portions of incompletely annotated genes, antisense transcripts, duplicated gene fragments, and potential non-coding RNAs. These results provide a systematic characterization of non-exonic transcripts in rice and thus expand the current view of the complexity and dynamics of the rice transcriptome.

Transcription factors with a large number of target genes--transcription hub(s), or THub(s)--are usually crucial components of the regulatory system of a cell, and the different patterns through which they transfer the transcriptional signal to downstream cascades are of great interest. By profiling normalized abundances (A(N)) of basic regulatory patterns of individual THubs in the yeast Saccharomyces cerevisiae transcriptional regulation network under five different cellular states and environmental conditions, we have investigated their preferences for different basic regulatory patterns. Subgraph-normalized abundances downstream of individual THubs often differ significantly from that of the network as a whole, and conversely, certain over-represented subgraphs are not preferred by any THub. The THub preferences changed substantially when the cellular or environmental conditions changed. This switching of regulatory pattern preferences suggests that a change in conditions does not only elicit a change in response by the regulatory network, but also a change in the mechanisms by which the response is mediated. The THub subgraph preference profile thus provides a novel tool for description of the structure and organization between the large-scale exponents and local regulatory patterns.