Polyadenylation ligation-mediated sequencing (PALM-Seq) characterizes cell-free coding and non-coding RNAs in human biofluids.

BACKGROUND: Cell-free messenger RNA (cf-mRNA) and long non-coding RNA (cf-lncRNA) are becoming increasingly important in liquid biopsy by providing biomarkers for disease prediction, diagnosis and prognosis, but the simultaneous characterization of coding and non-coding RNAs in human biofluids remains challenging. METHODS: Here, we developed polyadenylation ligation-mediated sequencing (PALM-Seq), an RNA sequencing strategy employing treatment of RNA with T4 polynucleotide kinase to generate cell-free RNA (cfRNA) fragments with 5' phosphate and 3' hydroxyl and RNase H to deplete abundant RNAs, achieving simultaneous quantification and characterization of cfRNAs. RESULTS: Using PALM-Seq, we successfully identified well-known differentially abundant mRNA, lncRNA and microRNA in the blood plasma of pregnant women. We further characterized cfRNAs in blood plasma, saliva, urine, seminal plasma and amniotic fluid and found that the detected numbers of different RNA biotypes varied with body fluids. The profiles of cf-mRNA reflected the function of originated tissues, and immune cells significantly contributed RNA to blood plasma and saliva. Short fragments (<50 nt) of mRNA and lncRNA were major in biofluids, whereas seminal plasma and amniotic fluid tended to retain long RNA. Body fluids showed distinct preferences of pyrimidine at the 3' end and adenine at the 5' end of cf-mRNA and cf-lncRNA, which were correlated with the proportions of short fragments. CONCLUSION: Together, PALM-Seq enables a simultaneous characterization of cf-mRNA and cf-lncRNA, contributing to elucidating the biology and promoting the application of cfRNAs.

Genetic characteristics of human papillomavirus type 16, 18, 52 and 58 in southern China.

Persistent infections of high-risk human papillomaviruses (HPVs) are the leading cause of cervical cancers. We collected cervical exfoliated cell samples from females in Changsha city, Hunan Province and obtained 338 viral genomes of four major HPV types, including HPV 16 (n = 82), 18 (n = 35), 52 (n = 121) and 58 (n = 100). The lineage/sublineage distribution of the four HPVs confirmed previous epidemiological reports, with the predominant prevailing sublineage as A4 (50%), A1 (37%) and A3 (13%) for HPV16, A1 (83%) for HPV18, B2 (86%) for HPV52 and A1 (65%), A3 (19%) and A2 (12%) for HPV58. We also identified two potentially novel HPV18 sublineages, i.e. A6 and A7. Virus mutation analysis further revealed the presence of HPV16 and HPV58 sublineages associated with potentially high oncogenicity. These findings expanded our knowledge of the HPV genetic diversity in China, providing valuable evidence to facilitate HPV DNA screening, vaccine effectiveness evaluation and control strategy development.

Transcriptional Start Site Coverage Analysis in Plasma Cell-Free DNA Reveals Disease Severity and Tissue Specificity of COVID-19 Patients.

Symptoms of coronavirus disease 2019 (COVID-19) range from asymptomatic to severe pneumonia and death. A deep understanding of the variation of biological characteristics in severe COVID-19 patients is crucial for the detection of individuals at high risk of critical condition for the clinical management of the disease. Herein, by profiling the gene expression spectrum deduced from DNA coverage in regions surrounding transcriptional start site in plasma cell-free DNA (cfDNA) of COVID-19 patients, we deciphered the altered biological processes in the severe cases and demonstrated the feasibility of cfDNA in measuring the COVID-19 progression. The up- and downregulated genes in the plasma of severe patient were found to be closely related to the biological processes and functions affected by COVID-19 progression. More importantly, with the analysis of transcriptome data of blood cells and lung cells from control group and cases with severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2) infection, we revealed that the upregulated genes were predominantly involved in the viral and antiviral activity in blood cells, reflecting the intense viral replication and the active reaction of immune system in the severe patients. Pathway analysis of downregulated genes in plasma DNA and lung cells also demonstrated the diminished adenosine triphosphate synthesis function in lung cells, which was evidenced to correlate with the severe COVID-19 symptoms, such as a cytokine storm and acute respiratory distress. Overall, this study revealed tissue involvement, provided insights into the mechanism of COVID-19 progression, and highlighted the utility of cfDNA as a noninvasive biomarker for disease severity inspections.

A Comprehensive RNA Study to Identify circRNA and miRNA Biomarkers for Docetaxel Resistance in Breast Cancer.

To investigate the relationship between non-coding RNAs [especially circular RNAs (circRNAs)] and docetaxel resistance in breast cancer, and to find potential predictive biomarkers for taxane-containing therapies, we have performed transcriptome and microRNA (miRNA) sequencing for two established docetaxel-resistant breast cancer (DRBC) cell lines and their docetaxel-sensitive parental cell lines. Our analyses revealed differences between circRNA signatures in the docetaxel-resistant and -sensitive breast cancer cells, and discovered circRNAs generated by multidrug-resistance genes in taxane-resistant cancer cells. In DRBC cells, circABCB1 was identified and validated as a circRNA that is strongly up-regulated, whereas circEPHA3.1 and circEPHA3.2 are strongly down-regulated. Furthermore, we investigated the potential functions of these circRNAs by bioinformatics analysis, and miRNA analysis was performed to uncover potential interactions between circRNAs and miRNAs. Our data showed that circABCB1, circEPHA3.1 and circEPHA3.2 may sponge up eight significantly differentially expressed miRNAs that are associated with chemotherapy and contribute to docetaxel resistance via the PI3K-Akt and AGE-RAGE signaling pathways. We also integrated differential expression data of mRNA, long non-coding RNA, circRNA, and miRNA to gain a global profile of multi-level RNA changes in DRBC cells, and compared them with changes in DNA copy numbers in the same cell lines. We found that Chromosome 7 q21.12-q21.2 was a common region dominated by multi-level RNA overexpression and DNA amplification, indicating that overexpression of the RNA molecules transcribed from this region may result from DNA amplification during stepwise exposure to docetaxel. These findings may help to further our understanding of the mechanisms underlying docetaxel resistance in breast cancer.

Genetic signatures for lineage/sublineage classification of HPV16, 18, 52 and 58 variants.

Increasing evidences indicate that high-risk HPV variants are heterogeneous in carcinogenicity and ethnic dispersion. In this work, we identified genetic signatures for convenient determination of lineage/sublineage of HPV16, 18, 52 and 58 variants. Using publicly available genomes, we found that E2 of HPV16, L2 of HPV18, L1 and LCR of HPV52, and L2, LCR and E1 of HPV58 contain the proper genetic signature for lineage/sublineage classification. Sets of hierarchical signature nucleotide positions were further confirmed for high accuracy (>95%) by classifying HPV genomes obtained from Chinese females, which included 117 HPV16 variants, 48 HPV18 variants, 117 HPV52 variants and 89 HPV58 variants. The circulation of HPV variants posing higher cancer risk in Eastern China, such as HPV16 A4 and HPV58 A3, calls for continuous surveillance in this region. The marker genes and signature nucleotide positions may facilitate cost-effective diagnostic detections of HPV variants in clinical settings.

Performance of copy number variants detection based on whole-genome sequencing by DNBSEQ platforms.

BACKGROUND: DNBSEQ™ platforms are new massively parallel sequencing (MPS) platforms that use DNA nanoball technology. Use of data generated from DNBSEQ™ platforms to detect single nucleotide variants (SNVs) and small insertions and deletions (indels) has proven to be quite effective, while the feasibility of copy number variants (CNVs) detection is unclear. RESULTS: Here, we first benchmarked different CNV detection tools based on Illumina whole-genome sequencing (WGS) data of NA12878 and then assessed these tools in CNV detection based on DNBSEQ™ sequencing data from the same sample. When the same tool was used, the CNVs detected based on DNBSEQ™ and Illumina data were similar in quantity, length and distribution, while great differences existed within results from different tools and even based on data from a single platform. We further estimated the CNV detection power based on available CNV benchmarks of NA12878 and found similar precision and sensitivity between the DNBSEQ™ and Illumina platforms. We also found higher precision of CNVs shorter than 1 kbp based on DNBSEQ™ platforms than those based on Illumina platforms by using Pindel, DELLY and LUMPY. We carefully compared these two available benchmarks and found a large proportion of specific CNVs between them. Thus, we constructed a more complete CNV benchmark of NA12878 containing 3512 CNV regions. CONCLUSIONS: We assessed and benchmarked CNV detections based on WGS with DNBSEQ™ platforms and provide guidelines for future studies.

Impact of sequencing depth and technology on de novo RNA-Seq assembly.

BACKGROUND: RNA-Seq data is inherently nonuniform for different transcripts because of differences in gene expression. This makes it challenging to decide how much data should be generated from each sample. How much should one spend to recover the less expressed transcripts? The sequencing technology used is another consideration, as there are inevitably always biases against certain sequences. To investigate these effects, we first looked at high-depth libraries from a set of well-annotated organisms to ascertain the impact of sequencing depth on de novo assembly. We then looked at libraries sequenced from the Universal Human Reference RNA (UHRR) to compare the performance of Illumina HiSeq and MGI DNBseq™ technologies. RESULTS: On the issue of sequencing depth, the amount of exomic sequence assembled plateaued using data sets of approximately 2 to 8 Gbp. However, the amount of genomic sequence assembled did not plateau for many of the analyzed organisms. Most of the unannotated genomic sequences are single-exon transcripts whose biological significance will be questionable for some users. On the issue of sequencing technology, both of the analyzed platforms recovered a similar number of full-length transcripts. The missing "gap" regions in the HiSeq assemblies were often attributed to higher GC contents, but this may be an artefact of library preparation and not of sequencing technology. CONCLUSIONS: Increasing sequencing depth beyond modest data sets of less than 10 Gbp recovers a plethora of single-exon transcripts undocumented in genome annotations. DNBseq™ is a viable alternative to HiSeq for de novo RNA-Seq assembly.

Intra-host Ebola viral adaption during human infection.

The onsite next generation sequencing (NGS) of Ebola virus (EBOV) genomes during the 2013-2016 Ebola epidemic in Western Africa provides an opportunity to trace the origin, transmission, and evolution of this virus. Herein, we have diagnosed a cohort of EBOV patients in Sierra Leone in 2015, during the late phase of the outbreak. The surviving EBOV patients had a recovery process characterized by decreasing viremia, fever, and biochemical parameters. EBOV genomes sequenced through the longitudinal blood samples of these patients showed dynamic intra-host substitutions of the virus during acute infection, including the previously described short stretches of 13 serial T>C mutations. Remarkably, within individual patients, samples collected during the early phase of infection possessed Ts at these nucleotide sites, whereas they were replaced by Cs in samples collected in the later phase, suggesting that these short stretches of T>C mutations could emerge independently. In addition, up to a total of 35 nucleotide sites spanning the EBOV genome were mutated coincidently. Our study showed the dynamic intra-host adaptation of EBOV during patient recovery and gave more insight into the complex EBOV-host interactions.

Cytokine-induced killer cell/dendritic cell-cytokine-induced killer cell immunotherapy for the postoperative treatment of gastric cancer: A systematic review and meta-analysis.

BACKGROUND: Immunotherapy is emerging as a new treatment strategy for gastric cancer(GC). However, the efficacy and safety of this technique remain unclear. This meta-analysis aimed to assess the effect of cytokine-induced killer cell (CIK)/dendritic cell-cytokine-induced killer cell (DC-CIK) treatment for GC after surgery. METHODS: Hazard ratio (HR), overall survival (OS) rates, and disease-free survival (DFS) rates were calculated using a Mantel-Haenszel (M-H) fixed-effects model (FEM), and results were displayed using forest plots. Publication bias was assessed by Begg test, and data were presented using funnel plots. Date robustness was assessed by the trim and fill method. Descriptive analysis was performed on T lymphocytes and adverse effects. RESULTS: In total, 9 trials, including 1216 patients, were eligible for inclusion in this meta-analysis. Compared with the control group, the HR for OS was 0.712 (95% confidence interval [CI] 0.594-0.854) and 0.66 (95% CI 0.546-0.797) for overall (DFS). The risk ratio (RR) of the 3 and 5-year OS rate was 1.29 (95% CI 1.15-1.46) and 1.73 (95% CI 1.36-2.19), respectively. The RR for the 3 and 5-year DFS rate 1.40 (95% CI 1.19-1.65) and 2.10 (95% CI1.53-2.87), respectively. The proportion of patients who were CD3+, CD4+, and CD4+/CD8+ increased in the cellular therapy groups. No fatal adverse reactions were noted. CONCLUSION: Chemotherapy combined with CIK/DC-CIK therapy after surgery resulted in low HR, and significantly increasing OS rates, DFS rates, and T-lymphocyte responses in patients with GC.

Germline and somatic variant identification using BGISEQ-500 and HiSeq X Ten whole genome sequencing.

Technological innovation and increased affordability have contributed to the widespread adoption of genome sequencing technologies in biomedical research. In particular large cancer research consortia have embraced next generation sequencing, and have used the technology to define the somatic mutation landscape of multiple cancer types. These studies have primarily utilised the Illumina HiSeq platforms. In this study we performed whole genome sequencing of three malignant pleural mesothelioma and matched normal samples using a new platform, the BGISEQ-500, and compared the results obtained with Illumina HiSeq X Ten. Germline and somatic, single nucleotide variants and small insertions or deletions were independently identified from data aligned human genome reference. The BGISEQ-500 and HiSeq X Ten platforms showed high concordance for germline calls with genotypes from SNP arrays (>99%). The germline and somatic single nucleotide variants identified in both sequencing platforms were highly concordant (86% and 72% respectively). These results indicate the potential applicability of the BGISEQ-500 platform for the identification of somatic and germline single nucleotide variants by whole genome sequencing. The BGISEQ-500 datasets described here represent the first publicly-available cancer genome sequencing performed using this platform.

Identification of balanced chromosomal rearrangements previously unknown among participants in the 1000 Genomes Project: implications for interpretation of structural variation in genomes and the future of clinical cytogenetics.

PURPOSE: Recent studies demonstrate that whole-genome sequencing enables detection of cryptic rearrangements in apparently balanced chromosomal rearrangements (also known as balanced chromosomal abnormalities, BCAs) previously identified by conventional cytogenetic methods. We aimed to assess our analytical tool for detecting BCAs in the 1000 Genomes Project without knowing which bands were affected. METHODS: The 1000 Genomes Project provides an unprecedented integrated map of structural variants in phenotypically normal subjects, but there is no information on potential inclusion of subjects with apparent BCAs akin to those traditionally detected in diagnostic cytogenetics laboratories. We applied our analytical tool to 1,166 genomes from the 1000 Genomes Project with sufficient physical coverage (8.25-fold). RESULTS: With this approach, we detected four reciprocal balanced translocations and four inversions, ranging in size from 57.9 kb to 13.3 Mb, all of which were confirmed by cytogenetic methods and polymerase chain reaction studies. One of these DNAs has a subtle translocation that is not readily identified by chromosome analysis because of the similarity of the banding patterns and size of exchanged segments, and another results in disruption of all transcripts of an OMIM gene. CONCLUSION: Our study demonstrates the extension of utilizing low-pass whole-genome sequencing for unbiased detection of BCAs including translocations and inversions previously unknown in the 1000 Genomes Project.

Comparative performance of the BGISEQ-500 vs Illumina HiSeq2500 sequencing platforms for palaeogenomic sequencing.

Ancient DNA research has been revolutionized following development of next-generation sequencing platforms. Although a number of such platforms have been applied to ancient DNA samples, the Illumina series are the dominant choice today, mainly because of high production capacities and short read production. Recently a potentially attractive alternative platform for palaeogenomic data generation has been developed, the BGISEQ-500, whose sequence output are comparable with the Illumina series. In this study, we modified the standard BGISEQ-500 library preparation specifically for use on degraded DNA, then directly compared the sequencing performance and data quality of the BGISEQ-500 to the Illumina HiSeq2500 platform on DNA extracted from 8 historic and ancient dog and wolf samples. The data generated were largely comparable between sequencing platforms, with no statistically significant difference observed for parameters including level (P = 0.371) and average sequence length (P = 0718) of endogenous nuclear DNA, sequence GC content (P = 0.311), double-stranded DNA damage rate (v. 0.309), and sequence clonality (P = 0.093). Small significant differences were found in single-strand DNA damage rate (δS; slightly lower for the BGISEQ-500, P = 0.011) and the background rate of difference from the reference genome (θ; slightly higher for BGISEQ-500, P = 0.012). This may result from the differences in amplification cycles used to polymerase chain reaction-amplify the libraries. A significant difference was also observed in the mitochondrial DNA percentages recovered (P = 0.018), although we believe this is likely a stochastic effect relating to the extremely low levels of mitochondria that were sequenced from 3 of the samples with overall very low levels of endogenous DNA. Although we acknowledge that our analyses were limited to animal material, our observations suggest that the BGISEQ-500 holds the potential to represent a valid and potentially valuable alternative platform for palaeogenomic data generation that is worthy of future exploration by those interested in the sequencing and analysis of degraded DNA.

A reference human genome dataset of the BGISEQ-500 sequencer.

Background: BGISEQ-500 is a new desktop sequencer developed by BGI. Using DNA nanoball and combinational probe anchor synthesis developed from Complete Genomics™ sequencing technologies, it generates short reads at a large scale. Here, we present the first human whole-genome sequencing dataset of BGISEQ-500. The dataset was generated by sequencing the widely used cell line HG001 (NA12878) in two sequencing runs of paired-end 50 bp (PE50) and two sequencing runs of paired-end 100 bp (PE100). We also include examples of the raw images from the sequencer for reference. Finally, we identified variations using this dataset, estimated the accuracy of the variations, and compared to that of the variations identified from similar amounts of publicly available HiSeq2500 data. We found similar single nucleotide polymorphism (SNP) detection accuracy for the BGISEQ-500 PE100 data (false positive rate [FPR] = 0.00020%, sensitivity = 96.20%) compared to the PE150 HiSeq2500 data (FPR = 0.00017%, sensitivity = 96.60%) better SNP detection accuracy than the PE50 data (FPR = 0.0006%, sensitivity = 94.15%). But for insertions and deletions (indels), we found lower accuracy for BGISEQ-500 data (FPR = 0.00069% and 0.00067% for PE100 and PE50 respectively, sensitivity = 88.52% and 70.93%) than the HiSeq2500 data (FPR = 0.00032%, sensitivity = 96.28%). Our dataset can serve as the reference dataset, providing basic information not just for future development, but also for all research and applications based on the new sequencing platform.

cPAS-based sequencing on the BGISEQ-500 to explore small non-coding RNAs.

BACKGROUND: We present the first sequencing data using the combinatorial probe-anchor synthesis (cPAS)-based BGISEQ-500 sequencer. Applying cPAS, we investigated the repertoire of human small non-coding RNAs and compared it to other techniques. RESULTS: Starting with repeated measurements of different specimens including solid tissues (brain and heart) and blood, we generated a median of 30.1 million reads per sample. 24.1 million mapped to the human genome and 23.3 million to the miRBase. Among six technical replicates of brain samples, we observed a median correlation of 0.98. Comparing BGISEQ-500 to HiSeq, we calculated a correlation of 0.75. The comparability to microarrays was similar for both BGISEQ-500 and HiSeq with the first one showing a correlation of 0.58 and the latter one correlation of 0.6. As for a potential bias in the detected expression distribution in blood cells, 98.6% of HiSeq reads versus 93.1% of BGISEQ-500 reads match to the 10 miRNAs with highest read count. After using miRDeep2 and employing stringent selection criteria for predicting new miRNAs, we detected 74 high-likely candidates in the cPAS sequencing reads prevalent in solid tissues and 36 candidates prevalent in blood. CONCLUSIONS: While there is apparently no ideal platform for all challenges of miRNome analyses, cPAS shows high technical reproducibility and supplements the hitherto available platforms.

Benzo[a]pyrene promotes gastric cancer cell proliferation and metastasis likely through the Aryl hydrocarbon receptor and ERK-dependent induction of MMP9 and c-myc.

Gastric cancer (GC) is the fifth most common cancer worldwide and the third leading cause of global cancer-related death. Benzo[a]pyrene (BaP), a Group Ⅰ carcinogen categorized by the IARC, is a cumulative foodborne carcinogen and ubiquitous environmental pollutant with potent carcinogenic properties. However, the function and mechanism of BaP exposure on GC progression remains unclear. We investigated the role of BaP in human GC progression to identify potential mechanism underlining its carcinogenic activity. After exposure to various concentrations of BaP, human GC cells SGC-7901 and MNK-45 showed an increased capability of proliferation, migration and invasion. Further study indicated that BaP promotes the expression of matrix metalloproteinase-9 (MMP9) and c-myc at mRNA and protein level, and activates Aryl hydrocarbon receptor (AhR) and ERK pathway. Moreover, BaP-induced overexpression of MMP9 and c-myc were attenuated by the ERK inhibitor U0126 and AhR inhibitor resveratrol, respectively. These data suggest that BaP promotes proliferation and metastasis of GC cells through upregulation of MMP9 and c-myc expression, and this was likely mediated via the AhR and ERK signaling pathway.

OTG-snpcaller: an optimized pipeline based on TMAP and GATK for SNP calling from ion torrent data.

Because the new Proton platform from Life Technologies produced markedly different data from those of the Illumina platform, the conventional Illumina data analysis pipeline could not be used directly. We developed an optimized SNP calling method using TMAP and GATK (OTG-snpcaller). This method combined our own optimized processes, Remove Duplicates According to AS Tag (RDAST) and Alignment Optimize Structure (AOS), together with TMAP and GATK, to call SNPs from Proton data. We sequenced four sets of exomes captured by Agilent SureSelect and NimbleGen SeqCap EZ Kit, using Life Technology's Ion Proton sequencer. Then we applied OTG-snpcaller and compared our results with the results from Torrent Variants Caller. The results indicated that OTG-snpcaller can reduce both false positive and false negative rates. Moreover, we compared our results with Illumina results generated by GATK best practices, and we found that the results of these two platforms were comparable. The good performance in variant calling using GATK best practices can be primarily attributed to the high quality of the Illumina sequences.