We did whole-transcriptome sequencing and whole-genome sequencing on nine pairs of Hepatocellular carcinoma (HCC) tumors and matched adjacent tissues to identify RNA editing events. We identified mean 26,982 editing sites with mean 89.5% canonical A→G edits in each sample using an improved bioinformatics pipeline. The editing rate was significantly higher in tumors than adjacent normal tissues. Comparing the difference between tumor and normal tissues of each patient, we found 7 non-synonymous tissue specific editing events including 4 tumor-specific edits and 3 normal-specific edits in the coding region, as well as 292 edits varying in editing degree. The significant expression changes of 150 genes associated with RNA editing were found in tumors, with 3 of the 4 most significant genes being cancer related. Our results show that editing might be related to higher gene expression. These findings indicate that RNA editing modification may play an important role in the development of HCC.

Gene expression profiling is being widely applied in cancer research to identify biomarkers for clinical endpoint prediction. Since RNA-seq provides a powerful tool for transcriptome-based applications beyond the limitations of microarrays, we sought to systematically evaluate the performance of RNA-seq-based and microarray-based classifiers in this MAQC-III/SEQC study for clinical endpoint prediction using neuroblastoma as a model.

Eusocial insects have evolved the capacity to generate adults with distinct morphological, reproductive and behavioural phenotypes from the same genome. Recent studies suggest that RNA editing might enhance the diversity of gene products at the post-transcriptional level, particularly to induce functional changes in the nervous system. Using head samples from the leaf-cutting ant Acromyrmex echinatior, we compare RNA editomes across eusocial castes, identifying ca. 11,000 RNA editing sites in gynes, large workers and small workers. Those editing sites map to 800 genes functionally enriched for neurotransmission, circadian rhythm, temperature response, RNA splicing and carboxylic acid biosynthesis. Most A. echinatior editing sites are species specific, but 8-23% are conserved across ant subfamilies and likely to have been important for the evolution of eusociality in ants. The level of editing varies for the same site between castes, suggesting that RNA editing might be a general mechanism that shapes caste behaviour in ants.

Schizophrenia (SCZ) and bipolar disorder (BPD) are severe mental disorders with high heritability. Clinicians have long noticed the similarities of clinic symptoms between these disorders. In recent years, accumulating evidence indicates some shared genetic liabilities. However, what is shared remains elusive. In this study, we conducted whole transcriptome analysis of post-mortem brain tissues (cingulate cortex) from SCZ, BPD and control subjects, and identified differentially expressed genes in these disorders. We found 105 and 153 genes differentially expressed in SCZ and BPD, respectively. By comparing the t-test scores, we found that many of the genes differentially expressed in SCZ and BPD are concordant in their expression level (q⩽0.01, 53 genes; q⩽0.05, 213 genes; q⩽0.1, 885 genes). Using genome-wide association data from the Psychiatric Genomics Consortium, we found that these differentially and concordantly expressed genes were enriched in association signals for both SCZ (P<10(-7)) and BPD (P=0.029). To our knowledge, this is the first time that a substantially large number of genes show concordant expression and association for both SCZ and BPD. Pathway analyses of these genes indicated that they are involved in the lysosome, Fc gamma receptor-mediated phagocytosis, regulation of actin cytoskeleton pathways, along with several cancer pathways. Functional analyses of these genes revealed an interconnected pathway network centered on lysosomal function and the regulation of actin cytoskeleton. These pathways and their interacting network were principally confirmed by an independent transcriptome sequencing data set of the hippocampus. Dysregulation of lysosomal function and cytoskeleton remodeling has direct impacts on endocytosis, phagocytosis, exocytosis, vesicle trafficking, neuronal maturation and migration, neurite outgrowth and synaptic density and plasticity, and different aspects of these processes have been implicated in SCZ and BPD.

Congenital hearing loss affects approximately 1-2 infants out of every 1000, with 50% of the cases resulting from genetic factors. Targeted gene panels have been widely used for genetic diagnosis of hearing loss. This study aims to reveal new diagnoses via reanalyzing historical data of a multigene panel, and exam the reasons for new diagnoses.

With advances in massive parallel sequencing (MPS) technology, whole-genome sequencing (WGS) has gradually evolved into the first-tier diagnostic test for genetic disorders. However, deployment practice and pipeline testing for clinical WGS are lacking.

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.

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.

Preimplantation genetic testing for monogenic defects (PGT-M) has been available in clinical practice. This study aimed to validate the applicability of targeted capture sequencing in developing personalized PGT-M assay.

Though massively parallel sequencing has been widely applied to noninvasive prenatal screen for common trisomy, the clinical use of massively parallel sequencing to noninvasive prenatal diagnose monogenic disorders is limited. This study was to develop a method for directly determining paternal haplotypes for noninvasive prenatal diagnosis of monogenic disorders without requiring proband's samples.

Lynch syndrome is a genetically and clinically heterogeneous disorder; it is caused by a germline mutation in DNA mismatch repair (MMR) genes. Individuals with a heterozygous mutation in MLH1 have an increased risk for developing colorectal cancer. Here we described a 5-generation Chinese Lynch syndrome family with different severity and onset age. A novel heterozygous germline mutation (c.3G>T, p.Met1Ile) in MLH1 gene was discovered by next generation sequencing. Our study also revealed by qPCR that the MLH1 mRNA expression in peripheral blood of patients in this family was remarkably lower than that of the unaffected carriers and non-carriers. The research results indicated that the mRNA expression level may provide predictive suggestions of treatment and management for carriers with the initiation codon mutation of MLH1 in this family. Further studies are undertaken in this family as well as other families with Lynch syndrome to interrogate the exact reasons affecting the MLH1 mRNA expression level and whether mRNA expression in peripheral blood could be a significant factor for early diagnosis and surveillance of Lynch syndrome.

Panduratin 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.

RNA-Seq, a method using next generation sequencing technologies to sequence the transcriptome, facilitates genome-wide analysis of splice junction sites. In this paper, we introduce SOAPsplice, a robust tool to detect splice junctions using RNA-Seq data without using any information of known splice junctions. SOAPsplice uses a novel two-step approach consisting of first identifying as many reasonable splice junction candidates as possible, and then, filtering the false positives with two effective filtering strategies. In both simulated and real datasets, SOAPsplice is able to detect many reliable splice junctions with low false positive rate. The improvement gained by SOAPsplice, when compared to other existing tools, becomes more obvious when the depth of sequencing is low. SOAPsplice is freely available at http://soap.genomics.org.cn/soapsplice.html.

Genetically engineered mouse models have proven to be essential tools for unraveling fundamental aspects of cancer biology and for testing novel therapeutic strategies. To optimally serve these goals, it is essential that the mouse model faithfully recapitulates the human disease. Recently, novel mouse models for neuroblastoma have been developed. Here, we report on the further genomic characterization through exome sequencing and DNA copy number analysis of four of the currently available murine neuroblastoma model systems (ALK, Th-MYCN, Dbh-MYCN and Lin28b). The murine tumors revealed a low number of genomic alterations - in keeping with human neuroblastoma - and a positive correlation of the number of genetic lesions with the time to onset of tumor formation was observed. Gene copy number alterations are the hallmark of both murine and human disease and frequently affect syntenic genomic regions. Despite low mutational load, the genes mutated in murine disease were found to be enriched for genes mutated in human disease. Taken together, our study further supports the validity of the tested mouse models for mechanistic and preclinical studies of human neuroblastoma.

Inherited loss-of-function mutations in the tumor suppressor BRCA2 gene are associated with a high risk of ovarian cancer in the Chinese population. The current case report discusses a novel heterozygous insertion in BRCA2 gene, c.3195_3196insA, in a 54-year-old Chinese female with hereditary ovarian cancer. This frameshift mutation generates a premature stop codon at amino acid 1,076, which leads to a truncated BRCA2 protein instead of a wild-type BRCA2 protein with 3,418 amino acids. According to the Breast Cancer Information Core database, this mutation has not been previously reported. However, germline mutations of BRCA2 are a more prevalent cause of ovarian cancer in Chinese females compared with females in Western populations. The present study expands the mutational spectra of BRCA2 that is associated with ovarian cancer.

Spinal muscular atrophy (SMA) is one of the most common and severe genetic diseases. SMA carrier screening is an effective way to identify couples at risk of having affected children. Next-generation sequencing (NGS)-based expanded carrier screening could detect SMN1 gene copy number without extra experiment and with high cost performance. However, its performance has not been fully evaluated. Here we conducted a systematic comparative study to evaluate the performance of three common methods. 478 samples were analyzed with multiplex ligation probe amplification (MLPA), real-time quantitative polymerase chain reaction (qPCR) and NGS, simultaneously. Taking MLPA-based results as the reference, for 0 copy, 1 copy and ≥ 2 copy SMN1 analysis with NGS, the sensitivity, specificity and precision were all 100%. Using qPCR method, the sensitivity was 100%, 97.52% and 94.30%, respectively; 98.63%, 95.48% and 100% for specificity; and 72.72%, 88.72% and 100% for precision. NGS repeatability was higher than that of qPCR. Moreover, among three methods, NGS had the lowest retest rate. Thus, NGS is a relatively more reliable method for SMN1 gene copy number detection. In expanded carrier screening, compared with the combination of multiple methods, NGS method could reduce the test cost and simplify the screening process.

Low-pass genome sequencing (LP GS) has shown distinct advantages over traditional methods for the detection of mosaicism. However, no study has systematically evaluated the accuracy of LP GS in the detection of mosaic aneuploidies and copy number variants (CNVs) in prenatal diagnosis. Moreover, the influence of sequencing depth on mosaicism detection of LP GS has not been fully evaluated.

The landscape of genetic alterations in lung adenocarcinoma derived from Asian patients is largely uncharacterized. Here we present an integrated genomic and transcriptomic analysis of 335 primary lung adenocarcinomas and 35 corresponding lymph node metastases from Chinese patients. Altogether 13 significantly mutated genes are identified, including the most commonly mutated gene TP53 and novel mutation targets such as RHPN2, GLI3 and MRC2. TP53 mutations are furthermore significantly enriched in tumours from patients harbouring metastases. Genes regulating cytoskeleton remodelling processes are also frequently altered, especially in metastatic samples, of which the high expression level of IQGAP3 is identified as a marker for poor prognosis. Our study represents the first large-scale sequencing effort on lung adenocarcinoma in Asian patients and provides a comprehensive mutational landscape for both primary and metastatic tumours. This may thus form a basis for personalized medical care and shed light on the molecular pathogenesis of metastatic lung adenocarcinoma.

Expanded carrier screening (ECS) has been demonstrated to increase the detection rate of carriers compared with traditional tests. The aim of this study was to assess the potential value of ECS for clinical application in Southern China, a region with high prevalence of thalassemia and with diverse ethnic groups, and to provide a reference for future implementations in areas with similar population characteristics. A total of 10,476 prenatal/preconception couples from 34 self-reported ethnic groups were simultaneously tested and analyzed anonymously for 11 Mendelian disorders using targeted next-generation sequencing. Overall, 27.49% of individuals without self-reported family history of disorders were found to be carriers of at least 1 of the 11 conditions, and the carrier frequency varied greatly between ethnic groups, ranging from 4.15% to 81.35%. Furthermore, 255 couples (2.43%) were identified as carrier couples at an elevated risk having an affected baby, sixty-five of which would not have been identified through the existing screening strategy, which only detects thalassemia. The modeled risk of fetuses being affected by any of the selected disorders was 531 per 100,000 (95% CI, 497-567 per 100,000). Our data demonstrate the feasibility of ECS, and provide evidence that ECS is a promising alternative to traditional one-condition screening strategies. The lessons learned from this experience should be applicable for other countries or regions with diverse ethnic groups.

Understanding of the RNA editing process has been broadened considerably by the next generation sequencing technology; however, several issues regarding this regulatory step remain unresolved--the strategies to accurately delineate the editome, the mechanism by which its profile is maintained, and its evolutionary and functional relevance. Here we report an accurate and quantitative profile of the RNA editome for rhesus macaque, a close relative of human. By combining genome and transcriptome sequencing of multiple tissues from the same animal, we identified 31,250 editing sites, of which 99.8% are A-to-G transitions. We verified 96.6% of editing sites in coding regions and 97.5% of randomly selected sites in non-coding regions, as well as the corresponding levels of editing by multiple independent means, demonstrating the feasibility of our experimental paradigm. Several lines of evidence supported the notion that the adenosine deamination is associated with the macaque editome--A-to-G editing sites were flanked by sequences with the attributes of ADAR substrates, and both the sequence context and the expression profile of ADARs are relevant factors in determining the quantitative variance of RNA editing across different sites and tissue types. In support of the functional relevance of some of these editing sites, substitution valley of decreased divergence was detected around the editing site, suggesting the evolutionary constraint in maintaining some of these editing substrates with their double-stranded structure. These findings thus complement the "continuous probing" model that postulates tinkering-based origination of a small proportion of functional editing sites. In conclusion, the macaque editome reported here highlights RNA editing as a widespread functional regulation in primate evolution, and provides an informative framework for further understanding RNA editing in human.