RNA-SEQ ANALYSIS

1. Library preparation and Sequencing

From the RNA samples to the final data, each step, including sample test, library preparation, and sequencing, influences the quality of the data, and data quality directly impacts the analysis results. To guarantee the reliability of the data, quality control (QC) is performed at each step of the procedure. The workflow is as follows:

pipeline


Firstly, ribosomal RNA was removed from total RNA, followed by ethanol precipitation. After fragmentation, the first strand cDNA was synthesized using random hexamer primers. During the second strand cDNA synthesis, dUTPs were replaced with dTTPs in the reaction buffer. The directional library was ready after end repair, A-tailing, adapter ligation, size selection, USER enzyme digestion, amplification, and purification. Following is the workflow of library construction:


library2


The library was checked with Qubit and real-time PCR for quantification and bioanalyzer for size distribution detection. Quantified libraries will be pooled and sequenced on Illumina platforms, according to effective library concentration and data amount required.

2. Results

Data Quality Control: The “e” represents the sequence error rate and Qphred represents the base quality value, Qphred=-10log10(e). The relationship between sequencing error rate (e) and sequencing base quality value (Qphred) is shown in the table below:

Phred scoreError baseRight baseQ-score
101/1090%Q10
201/10099%Q20
301/100099.9%Q30
401/1000099.99%Q40

The distribution of quality score is shown in Fig.1:


VN_0005.QM
VN_0007.QM

Summary of Sequencing Data Information

SampleRaw readsRaw dataEffective(%)Error(%)Q20(%)Q30(%)GC(%)
VN_09207427800411.197.310.0298.0194.4247.88
VN_08226896122210.397.570.0298.0194.3347.51
VN_04136990336410.597.740.0397.7193.6146.46
VN_00126748684210.197.840.0397.8394.0145.30
VN_09176860536610.397.520.0397.8393.9547.16
VN_10827549502611.396.020.0397.0692.1747.00
VN_09257038546610.697.300.0396.8991.8048.07
VN_101310130967415.275.180.0493.2086.6164.00
VN_009010626980215.980.950.0493.4186.6565.99
VN_08126791778410.297.000.0397.7693.6345.62
VN_09247011206810.597.660.0397.8093.7646.47
VN_00856975203810.596.300.0397.6493.6649.85
VN_00056938958210.497.310.0397.8994.1346.18
VN_01886720428210.197.850.0397.8193.8846.42
VN_00076698863410.096.130.0397.7193.7246.74
VN_10836830331210.298.230.0397.6693.3346.82