Whole-genome sequence of the planarian Dugesia japonica combining Illumina and PacBio data.

Advances in stem cell biology have posed the challenges in revealing the mechanistic themes underlying whole tissues and organs formation during regeneration. The planarian Dugesia japonica is an ideal model organism for the study of regeneration and stem cell biology. However, the genome resources for this species are still limited. Here, we combined single-molecule real-time DNA sequencing platform Pacific Biosciences (PacBio) sequencing, Illumina paired-end sequencing and 10× Genomics linked reads data to obtain the whole-genome sequence of the planarian D. japonica. The final assembled D. japonica genome is 1.13G with contig N50 of 248.44 kb, and scaffold N50 of 652.52 kb. Repeat elements account for 64.92% of the genome, and 12,031 protein coding genes were annotated, of which 10,114 genes had at least one functional annotation, representing 84.07% of the total genes. We present a highly contiguous genome assembly of D. japonica. The D. japonica genome assembly, together with gene annotation and transcriptome data provide a valuable resource to investigate molecular mechanism of planarian and stem cell research.

Author Correction: Genome-wide identification of Arabidopsis long noncoding RNAs in response to the blue light.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Genome-wide identification of Arabidopsis long noncoding RNAs in response to the blue light.

Long non-coding RNAs (lncRNAs) have been shown in animals to play roles in a wide range of biological processes. In plant, light modulates the growth and development as a key external signal. However, little is known about the role of plant lncRNA in response to light. In this study, we sequenced the messenger RNAs (mRNAs), lncRNAs and microRNAs (miRNAs) in Arabidopsis seedlings under blue light for 2 h and 8 h. Compared to dark, we identified 4197 mRNAs, 375 miRNAs and 481 lncRNAs, or 5207 mRNAs, 286 miRNAs and 545 lncRNAs of differential expressions under blue light treatments for 2 h or 8 h respectively. Subsequently, a total of 407 competing endogenous RNA (ceRNA) pairs (lncRNA-mRNA-miRNA) were constructed. We identified a blue light-induced lncRNA which plays roles in blue light-directed plant photomorphogenesis and response to mannitol stress by serving as a ceRNA to sequester miR167 in a type of target mimicry. These results revealed previously unknown roles of the lncRNA in blue light signaling and mannitol stress, and provided useful resources of lncRNAs associated with miRNAs in response to blue light.

Exploring transcription factors reveals crucial members and regulatory networks involved in different abiotic stresses in Brassica napus L.

BACKGROUND: Brassica napus (B. napus) encompasses diverse transcription factors (TFs), but thorough identification and characterization of TF families, as well as their transcriptional responsiveness to multifarious stresses are still not clear. RESULTS: Totally 2167 TFs belonging to five families were genome-widely identified in B. napus, including 518 BnAP2/EREBPs, 252 BnbZIPs, 721 BnMYBs, 398 BnNACs and 278 BnWRKYs, which contained some novel members in comparison with existing results. Sub-genome distributions of BnAP2/EREBPs and BnMYBs indicated that the two families might have suffered from duplication and divergence during evolution. Synteny analysis revealed strong co-linearity between B. napus and its two ancestors, although chromosomal rearrangements have occurred and 85 TFs were lost. About 7.6% and 9.4% TFs of the five families in B. napus were novel genes and conserved genes, which both showed preference on the C sub-genome. RNA-Seq revealed that more than 80% TFs were abiotic stress inducible and 315 crucial differentially expressed genes (DEGs) were screened out. Network analysis revealed that the 315 DEGs are highly co-expressed. The homologous gene network in A. thaliana revealed that a considerable amount of TFs could trigger the differential expression of targeted genes, resulting in a complex clustered network with clusters of genes responsible for targeted stress responsiveness. CONCLUSIONS: We identified and characterized five TF families in B. napus. Some crucial members and regulatory networks involved in different abiotic stresses have been explored. The investigations deepen our understanding of TFs for stress tolerance in B. napus.

Extensive ceRNA-ceRNA interaction networks mediated by miRNAs regulate development in multiple rhesus tissues.

Crosstalk between RNAs mediated by shared microRNAs (miRNAs) represents a novel layer of gene regulation, which plays important roles in development. In this study, we analyzed time series expression data for coding genes and long non-coding RNAs (lncRNAs) to identify thousands of interactions among competitive endogenous RNAs (ceRNAs) in four rhesus tissues. The ceRNAs exhibited dynamic expression and regulatory patterns during each tissue development process, which suggests that ceRNAs might work synergistically during different developmental stages or tissues to control specific functions. In addition, lncRNAs exhibit higher specificity as ceRNAs than coding-genes and their functions were predicted based on their competitive coding-gene partners to discover their important developmental roles. In addition to the specificity of tissue development, functional analyses demonstrated that the combined effects of multiple ceRNAs can have major impacts on general developmental and metabolic processes in multiple tissues, especially transcription-related functions where competitive interactions. Moreover, ceRNA interactions could sequentially and/or synergistically mediate the crosstalk among different signaling pathways during brain development. Analyzing ceRNA interactions during the development of multiple tissues will provideinsights in the regulation of normal development and the dysregulation of key mechanisms during pathogenesis.

Crosstalk of dynamic functional modules in lung development of rhesus macaques.

Lung development follows a complex series of dynamic histogenic events that depend on the fluctuation of gene expression, and the disruption of gene regulation could lead to devastating consequences, such as diseases in adulthood. In order to investigate the mechanism of lung development, we performed RNA sequencing by Illumina HiSeq™ 2000 to measure mRNA expression in lung tissues of nine rhesus macaques spanning from foetuses at gestation of 45 days to postnatal at 7 days. This development period was divided into three developmental stages, including the early stage (45-100 gestational days), the middle stage (137-163 gestational days) and the late stage (after birth at 4-7 days). Firstly, we identified stage-specific genes, based on which we found that the principle biological processes of the early stage were mainly associated with internal growth signalling, while the middle and late stage-specific genes controlled the external stress signalling. Then, we constructed a stage-specific protein-protein interaction (PPI) subnetwork, extracted dynamic modules, and identified crosstalk between modules. Moreover, we found four active pathways that could mediate the crosstalk, including the Notch signalling pathway, cell cycle, NOD-like receptor signalling pathway, and Toll-like receptor signalling pathway. These pathways not only played crucial roles in lung development, but also were implicated in lung diseases. Finally, some important bridgers, such as PSEN2, HSP90AA1 and CASP8, were discovered to explain the potential mechanism of crosstalk. Therefore, our study presents the landscape of gene expression of lung development of rhesus macaques, and provides an extended insight into the lung development mechanism.

Quantitative proteomics reveals the novel co-expression signatures in early brain development for prognosis of glioblastoma multiforme.

Although several researches have explored the similarity across development and tumorigenesis in cellular behavior and underlying molecular mechanisms, not many have investigated the developmental characteristics at proteomic level and further extended to cancer clinical outcome. In this study, we used iTRAQ to quantify the protein expression changes during macaque rhesus brain development from fetuses at gestation 70 days to after born 5 years. Then, we performed weighted gene co-expression network analysis (WGCNA) on protein expression data of brain development to identify co-expressed modules that highly expressed on distinct development stages, including early stage, middle stage and late stage. Moreover, we used the univariate cox regression model to evaluate the prognostic potentials of these genes in two independent glioblastoma multiforme (GBM) datasets. The results showed that the modules highly expressed on early stage contained more reproducible prognostic genes, including ILF2, CCT7, CCT4, RPL10A, MSN, PRPS1, TFRC and APEX1. These genes were not only associated with clinical outcome, but also tended to influence chemoresponse. These signatures identified from embryonic brain development might contribute to precise prediction of GBM prognosis and identification of novel drug targets in GBM therapies. Thus, the development could become a viable reference model for researching cancers, including identifying novel prognostic markers and promoting new therapies.

Gekko japonicus genome reveals evolution of adhesive toe pads and tail regeneration.

Reptiles are the most morphologically and physiologically diverse tetrapods, and have undergone 300 million years of adaptive evolution. Within the reptilian tetrapods, geckos possess several interesting features, including the ability to regenerate autotomized tails and to climb on smooth surfaces. Here we sequence the genome of Gekko japonicus (Schlegel's Japanese Gecko) and investigate genetic elements related to its physiology. We obtain a draft G. japonicus genome sequence of 2.55 Gb and annotated 22,487 genes. Comparative genomic analysis reveals specific gene family expansions or reductions that are associated with the formation of adhesive setae, nocturnal vision and tail regeneration, as well as the diversification of olfactory sensation. The obtained genomic data provide robust genetic evidence of adaptive evolution in reptiles.

Spatiotemporal-specific lncRNAs in the brain, colon, liver and lung of macaque during development.

Genome-wide expression profiling during development provides useful information to uncover the potential molecular mechanisms of development in mammals. Recent studies have revealed that a subset of lncRNAs can regulate major biological processes during development. Here, we sequenced four tissues, including brain, colon, liver and lung, using RNA-seq across three developmental stages (early, middle and late stage), and then constructed genome-wide expression profiles during macaque development. In each tissue, we identified developmental time-specific lncRNA and mRNA clusters that displayed diverse expression alteration patterns, including a gradual increase, a gradual decrease, or a reversal of expression. These lncRNAs showed more specific functional associations with their corresponding tissues relative to the developmental time-specific mRNAs. Furthermore, we identified 20 spatiotemporal-specific co-modules including 101 lncRNAs and 609 mRNAs distributed at different developmental stages in different tissues. Our findings suggested that lncRNAs could play critical roles in the development of macaques through close cooperation with mRNAs. Finally, we predicted the functions of the spatiotemporal-specific lncRNAs by their spatiotemporal cooperation with mRNAs and further validated our findings using gene knockdown data of mouse. Our study reveals the spatiotemporal characteristics of lncRNAs and provides a functional map of the spatiotemporal-specific lncRNAs during the development of macaques.

A host plant genome (Zizania latifolia) after a century-long endophyte infection.

Despite the importance of host-microbe interactions in natural ecosystems, agriculture and medicine, the impact of long-term (especially decades or longer) microbial colonization on the dynamics of host genomes is not well understood. The vegetable crop 'Jiaobai' with enlarged edible stems was domesticated from wild Zizania latifolia (Oryzeae) approximately 2000 years ago as a result of persistent infection by a fungal endophyte, Ustilago esculenta. Asexual propagation via infected rhizomes is the only means of Jiaobai production, and the Z. latifolia-endophyte complex has been maintained continuously for two centuries. Here, genomic analysis revealed that cultivated Z. latifolia has a significantly smaller repertoire of immune receptors compared with wild Z. latifolia. There are widespread gene losses/mutations and expression changes in the plant-pathogen interaction pathway in Jiaobai. These results show that continuous long-standing endophyte association can have a major effect on the evolution of the structural and transcriptomic components of the host genome.

Molecular cloning and characterization of PtrLAR3, a gene encoding leucoanthocyanidin reductase from Populus trichocarpa, and its constitutive expression enhances fungal resistance in transgenic plants.

The flavonoid-derived proanthocyanidins (PAs) are one class of the major defence phenolics in poplar leaves. Transcriptional activation of PA biosynthetic genes, resulting in PA accumulation in leaves, was detected following infection by the fungal Marssonina brunnea f.sp. multigermtubi using digital gene expression analysis. In order to study PA biosynthesis and its induction by fungi, a putative leucoanthocyanidin reductase gene, PtrLAR3, was isolated from Populus trichocarpa. Sequence comparison of PtrLAR3 with other known leucoanthocyanidin reductase proteins revealed high amino acid sequence similarity. Semi-quantitative reverse-transcription (RT) PCR and quantitative real-time PCR analysis demonstrated that PtrLAR3 was expressed in various tissues and the highest level of expression was observed in roots. Overexpression of PtrLAR3 in Chinese white poplar (Populus tomentosa Carr.) led to a significant plant-wide increase in PA levels. In vitro assays showed that crude leaf extracts from 35S:PtrLAR3 transformants were able to inhibit significantly the hyphal growth of M. brunnea f.sp. multigermtubi compared to the extracts from control plants. The transgenic 35S:PtrLAR3 poplar plants displayed a significant (P < 0.05) reduction in their disease symptoms compared with the control. RT-PCR analysis showed that PtrLAR3 expression was up-regulated in all transformants. These results suggested that constitutive expression of endogenous PtrLAR3 could be exploited to improve resistance to fungal pathogens in poplar.

A systematic analysis on DNA methylation and the expression of both mRNA and microRNA in bladder cancer.

BACKGROUND: DNA methylation aberration and microRNA (miRNA) deregulation have been observed in many types of cancers. A systematic study of methylome and transcriptome in bladder urothelial carcinoma has never been reported. METHODOLOGY/PRINCIPAL FINDINGS: The DNA methylation was profiled by modified methylation-specific digital karyotyping (MMSDK) and the expression of mRNAs and miRNAs was analyzed by digital gene expression (DGE) sequencing in tumors and matched normal adjacent tissues obtained from 9 bladder urothelial carcinoma patients. We found that a set of significantly enriched pathways disrupted in bladder urothelial carcinoma primarily related to "neurogenesis" and "cell differentiation" by integrated analysis of -omics data. Furthermore, we identified an intriguing collection of cancer-related genes that were deregulated at the levels of DNA methylation and mRNA expression, and we validated several of these genes (HIC1, SLIT2, RASAL1, and KRT17) by Bisulfite Sequencing PCR and Reverse Transcription qPCR in a panel of 33 bladder cancer samples. CONCLUSIONS/SIGNIFICANCE: We characterized the profiles between methylome and transcriptome in bladder urothelial carcinoma, identified a set of significantly enriched key pathways, and screened four aberrantly methylated and expressed genes. Conclusively, our findings shed light on a new avenue for basic bladder cancer research.

Comparative mRNA and microRNA expression profiling of three genitourinary cancers reveals common hallmarks and cancer-specific molecular events.

BACKGROUND: Genome-wide gene expression profile using deep sequencing technologies can drive the discovery of cancer biomarkers and therapeutic targets. Such efforts are often limited to profiling the expression signature of either mRNA or microRNA (miRNA) in a single type of cancer. METHODOLOGY: Here we provided an integrated analysis of the genome-wide mRNA and miRNA expression profiles of three different genitourinary cancers: carcinomas of the bladder, kidney and testis. PRINCIPAL FINDINGS: Our results highlight the general or cancer-specific roles of several genes and miRNAs that may serve as candidate oncogenes or suppressors of tumor development. Further comparative analyses at the systems level revealed that significant aberrations of the cell adhesion process, p53 signaling, calcium signaling, the ECM-receptor and cell cycle pathways, the DNA repair and replication processes and the immune and inflammatory response processes were the common hallmarks of human cancers. Gene sets showing testicular cancer-specific deregulation patterns were mainly implicated in processes related to male reproductive function, and general disruptions of multiple metabolic pathways and processes related to cell migration were the characteristic molecular events for renal and bladder cancer, respectively. Furthermore, we also demonstrated that tumors with the same histological origins and genes with similar functions tended to group together in a clustering analysis. By assessing the correlation between the expression of each miRNA and its targets, we determined that deregulation of 'key' miRNAs may result in the global aberration of one or more pathways or processes as a whole. CONCLUSIONS: This systematic analysis deciphered the molecular phenotypes of three genitourinary cancers and investigated their variations at the miRNA level simultaneously. Our results provided a valuable source for future studies and highlighted some promising genes, miRNAs, pathways and processes that may be useful for diagnostic or therapeutic applications.

Integrated profiling of microRNAs and mRNAs: microRNAs located on Xq27.3 associate with clear cell renal cell carcinoma.

BACKGROUND: With the advent of second-generation sequencing, the expression of gene transcripts can be digitally measured with high accuracy. The purpose of this study was to systematically profile the expression of both mRNA and miRNA genes in clear cell renal cell carcinoma (ccRCC) using massively parallel sequencing technology. METHODOLOGY: The expression of mRNAs and miRNAs were analyzed in tumor tissues and matched normal adjacent tissues obtained from 10 ccRCC patients without distant metastases. In a prevalence screen, some of the most interesting results were validated in a large cohort of ccRCC patients. PRINCIPAL FINDINGS: A total of 404 miRNAs and 9,799 mRNAs were detected to be differentially expressed in the 10 ccRCC patients. We also identified 56 novel miRNA candidates in at least two samples. In addition to confirming that canonical cancer genes and miRNAs (including VEGFA, DUSP9 and ERBB4; miR-210, miR-184 and miR-206) play pivotal roles in ccRCC development, promising novel candidates (such as PNCK and miR-122) without previous annotation in ccRCC carcinogenesis were also discovered in this study. Pathways controlling cell fates (e.g., cell cycle and apoptosis pathways) and cell communication (e.g., focal adhesion and ECM-receptor interaction) were found to be significantly more likely to be disrupted in ccRCC. Additionally, the results of the prevalence screen revealed that the expression of a miRNA gene cluster located on Xq27.3 was consistently downregulated in at least 76.7% of ∼50 ccRCC patients. CONCLUSIONS: Our study provided a two-dimensional map of the mRNA and miRNA expression profiles of ccRCC using deep sequencing technology. Our results indicate that the phenotypic status of ccRCC is characterized by a loss of normal renal function, downregulation of metabolic genes, and upregulation of many signal transduction genes in key pathways. Furthermore, it can be concluded that downregulation of miRNA genes clustered on Xq27.3 is associated with ccRCC.