We present primary results from the Sequencing Quality Control (SEQC) project, coordinated by the United States Food and Drug Administration. Examining Illumina HiSeq, Life Technologies SOLiD and Roche 454 platforms at multiple laboratory sites using reference RNA samples with built-in controls, we assess RNA sequencing (RNA-seq) performance for junction discovery and differential expression profiling and compare it to microarray and quantitative PCR (qPCR) data using complementary metrics. At all sequencing depths, we discover unannotated exon-exon junctions, with >80% validated by qPCR. We find that measurements of relative expression are accurate and reproducible across sites and platforms if specific filters are used. In contrast, RNA-seq and microarrays do not provide accurate absolute measurements, and gene-specific biases are observed, for these and qPCR. Measurement performance depends on the platform and data analysis pipeline, and variation is large for transcript-level profiling. The complete SEQC data sets, comprising >100 billion reads (10Tb), provide unique resources for evaluating RNA-seq analyses for clinical and regulatory settings.

Hepatocellular carcinoma (HCC) is one of the most deadly cancers worldwide and has no effective treatment, yet the molecular basis of hepatocarcinogenesis remains largely unknown. Here we report findings from a whole-genome sequencing (WGS) study of 88 matched HCC tumor/normal pairs, 81 of which are Hepatitis B virus (HBV) positive, seeking to identify genetically altered genes and pathways implicated in HBV-associated HCC. We find beta-catenin to be the most frequently mutated oncogene (15.9%) and TP53 the most frequently mutated tumor suppressor (35.2%). The Wnt/beta-catenin and JAK/STAT pathways, altered in 62.5% and 45.5% of cases, respectively, are likely to act as two major oncogenic drivers in HCC. This study also identifies several prevalent and potentially actionable mutations, including activating mutations of Janus kinase 1 (JAK1), in 9.1% of patients and provides a path toward therapeutic intervention of the disease.

Elucidating the molecular basis of hepatocellular carcinoma (HCC) is crucial to developing targeted diagnostics and therapies for this deadly disease. The landscape of somatic genomic rearrangements (GRs), which can lead to oncogenic gene fusions, remains poorly characterized in HCC. We have predicted 4314 GRs including large-scale insertions, deletions, inversions and translocations based on the whole-genome sequencing data for 88 primary HCC tumor/non-tumor tissues. We identified chromothripsis in 5 HCC genomes (5.7%) recurrently affecting chromosomal arms 1q and 8q. Albumin (ALB) was found to harbor GRs, deactivating mutations and deletions in 10% of cohort. Integrative analysis identified a pattern of paired intra-chromosomal translocations flanking focal amplifications and asymmetrical patterns of copy number variation flanking breakpoints of translocations. Furthermore, we predicted 260 gene fusions which frequently result in aberrant over-expression of the 3' genes in tumors and validated 18 gene fusions, including recurrent fusion (2/88) of ABCB11 and LRP2.

Upland and lowland rice (Oryza sativa L.) represent two of the most important rice ecotypes adapted to ago-ecosystems with contrasting soil-water conditions. Upland rice, domesticated in the water-limited environment, contains valuable drought-resistant characters that can be used in water-saving breeding. Knowledge about the divergence between upland and lowland rice will provide valuable cues for the evolution of drought-resistance in rice. Genetic differentiation between upland and lowland rice was explored by 47 Simple Sequence Repeats (SSRs) located in drought responding expressed sequence tags (ESTs) among 377 rice landraces. The morphological traits of drought-resistance were evaluated in the field experiments. Different outlier loci were detected in the japonica and indica subspecies, respectively. Considerable genetic differentiation between upland and lowland rice on these outlier loci was estimated in japonica (Fst = 0.258) and indica (Fst = 0.127). Furthermore, populations of the upland and lowland ecotypes were clustered separately on these outlier loci. A significant correlation between genetic distance matrices and the dissimilarity matrices of drought-resistant traits was determined, indicating a certain relationship between the upland-lowland rice differentiation and the drought-resistance. Divergent selections occur between upland and lowland rice on the drought-resistance as the Qsts of some drought-resistant traits are significantly higher than the neutral Fst. In addition, the upland- and lowland-preferable alleles responded differently among ecotypes or allelic types under osmotic stress. This shows the evolutionary signature of drought resistance at the gene expression level. The findings of this study can strengthen our understanding of the evolution of drought-resistance in rice with significant implications in the improvement of rice drought-resistance.

DNA methylation changes within the genome can be used to predict human age. However, the existing biological age prediction models based on DNA methylation are predominantly adult-oriented. We established a methylation-based age prediction model for children (9-212 months old) using data from 716 blood samples in 11 DNA methylation datasets. Our elastic net model includes 111 CpG sites, mostly in genes associated with development and aging. The model performed well and exhibited high precision, yielding a 98% correlation between the DNA methylation age and the chronological age, with an error of only 6.7 months. When we used the model to assess age acceleration in children based on their methylation data, we observed the following: first, the aging rate appears to be fastest in mid-childhood, and this acceleration is more pronounced in autistic children; second, lead exposure early in life increases the aging rate in boys, but not in girls; third, short-term recombinant human growth hormone treatment has little effect on the aging rate of children. Our child-specific methylation-based age prediction model can effectively detect epigenetic changes and health imbalances early in life. This may thus be a useful model for future studies of epigenetic interventions for age-related diseases.

BackgroundSmoking leads to the aging of organs. However, no studies have been conducted to quantify the effect of smoking on the aging of respiratory organs and the aging-reversing ability of smoking cessation.ResultsWe collected genome-wide methylation datasets of buccal cells, airway cells, esophagus tissue, and lung tissue from non-smokers, smokers, and ex-smokers. We used the “epigenetic clock” method to quantify the epigenetic age acceleration in the four organs. The statistical analyses showed the following: (1) Smoking increased the epigenetic age of airway cells by an average of 4.9 years and lung tissue by 4.3 years. (2) After smoking ceased, the epigenetic age acceleration in airway cells (but not in lung tissue) slowed to a level that non-smokers had. (3) The epigenetic age acceleration in airway cells and lung tissue showed no gender difference.ConclusionsSmoking can accelerate the epigenetic age of human respiratory organs, but the effect varies among organs and can be reversed by smoking cessation. Our study provides a powerful incentive to reduce tobacco consumption autonomously.

The aim ofthis study is to examine the mental health status ofyoung migrant workers in Shenzhen.Using the Syniptonis Check List-90 (SCL-90), Eysenck Personality Questionnaire, Social Support Scale and Mental Health Questionnaire for Laborers,371 migrant workers who came from inland areas ofChina and 100 local workers Xvere investigated.The SCL-90 profile of migrant workers was also compared to the SCL-90 n o r m provided by general people in China. The SCL-90 results showed that the total scores, the average scores of the positive symptoms, the three factor scores of obscssionality, interpersonal sensitivity and phobia in migrant workers were significantly higher than those in the local workers.According to the multivariate analysis, the aniount ofcontribution to mental health, in descending order, was neuroticisni, psychological pressure, income, home sickness, marital or love problems, extroversion and introversion, living conditions and social status.The mental health status of young migrant workers in Shenzhen was poorer than that of their local counterparts, as \vcU as people in China on the SCL-90. It is recommended that niental health workers should help migrant workers adjust to the new urban environment by providing psychological counseling and other relevant treatment facilities.,

BackgroundPanduratin A extracted from Boesenbergia rotunda is a flavonoid reported to possess a range of medicinal indications which include anti-dengue, anti-HIV, anti-cancer, antioxidant and anti-inflammatory properties. Boesenbergia rotunda is a plant from the Zingiberaceae family commonly used as a food ingredient and traditional medicine in Southeast Asia and China. Reports on the health benefits of secondary metabolites extracted from Boesenbergia rotunda over the last few years has resulted in rising demands for panduratin A. However large scale extraction has been hindered by the naturally low abundance of the compound and limited knowledge of its biosynthetic pathway.ResultsTranscriptome sequencing and digital gene expression (DGE) analysis of native and phenylalanine treated Boesenbergia rotunda cell suspension cultures were carried out to elucidate the key genes differentially expressed in the panduratin A biosynthetic pathway. Based on experiments that show increase in panduratin A production after 14 days post treatment with exogenous phenylalanine, an aromatic amino acid derived from the shikimic acid pathway, total RNA of untreated and 14 days post-phenylalanine treated cell suspension cultures were extracted and sequenced using next generation sequencing technology employing an Illumina-Solexa platform. The transcriptome data generated 101, 043 unigenes with 50, 932 (50.41%) successfully annotated in the public protein databases; including 49.93% (50, 447) in the non-redundant (NR) database, 34.63% (34, 989) in Swiss-Prot, 24,07% (24, 316) in Kyoto Encyclopedia of Genes and Genomes (KEGG) and 16.26% (16, 426) in Clusters of Orthologous Groups (COG). Through DGE analysis, we found that 14, 644 unigenes were up-regulated and 14, 379 unigenes down-regulated in response to exogenous phenylalanine treatment. In the phenylpropanoid pathway leading to the proposed panduratin A production, 2 up-regulated phenylalanine ammonia-lyase (PAL), 3 up-regulated 4-coumaroyl:coenzyme A ligase (4CL) and 1 up-regulated chalcone synthase (CHS) were found.ConclusionsThis is the first report of Boesenbergia rotunda de novo transcriptome data that could serve as a reference for gene or enzyme functional studies in the Zingiberaceae family. Although enzymes that are directly involved in the panduratin A biosynthetic pathway were not completely elucidated, the data provides an overall picture of gene regulation patterns leading to panduratin A production.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-984) contains supplementary material, which is available to authorized users.