Increasing quality, throughput and speed of sample preparation for strand-specific messenger RNA sequencing

Three genome assemblies of Opsariichthys bidens from Yangzte River, Pearl River and Qiantang River basins

Opsariichthys bidens is an endemic minnow, mainly distributed in China as an emerging aquaculture species. In this study, O. bidens collected from Yangzte River (YR), Pearl River (PR) and Qiantang River (QR) basins were used for karyotypic study and genome sequencing. Using PacBio long reads and Hi-C data, we assembled high-quality O. bidens genomes of 852.41 Mb (YR), 843.11 Mb (PR) and 840.94 Mb (QR) with scaffold N50 lengths of 21.01 Mb, 23.62 Mb and 24.75 Mb, respectively, of which 90.39%, 95.67% and 99.01% were anchored to 37, 38 and 39 chromosomes, respectively. 26,556 (YR), 25,036 (PR) and 26,283 (QR) protein-coding genes were predicted, respectively. The karyotype of O. bidens from YR, PR and QR were 2 N = 74 (6 m + 6sm + 4st + 58t), 2 N = 76 (4 m + 6sm + 4st + 62t) and 2 N = 78 (4 m + 4sm + 4st + 66t), respectively. Collinearity analysis and telomere predictions indicated that the observed chromosomal evolution was driven by Robertsonian translocation. These genome assemblies facilitate cryptic species determination, evolutionary study and genetic breeding of genus Opsariichthys.

Transposon insertions regulate genome-wide allele-specific expression and underpin flower colour variations in apple (Malus spp.)

Allele-specific expression (ASE) can lead to phenotypic diversity and evolution. However, the mechanisms regulating ASE are not well understood, particularly in woody perennial plants. In this study, we investigated ASE genes in the apple cultivar 'Royal Gala' (RG). A high quality chromosome-level genome was assembled using a homozygous tetra-haploid RG plant, derived from anther cultures. Using RNA-sequencing (RNA-seq) data from RG flower and fruit tissues, we identified 2091 ASE genes. Compared with the haploid genome of 'Golden Delicious' (GD), a parent of RG, we distinguished the genomic sequences between the two alleles of 817 ASE genes, and further identified allele-specific presence of a transposable element (TE) in the upstream region of 354 ASE genes. These included MYB110a that encodes a transcription factor regulating anthocyanin biosynthesis. Interestingly, another ASE gene, MYB10 also showed an allele-specific TE insertion and was identified using genome data of other apple cultivars. The presence of the TE insertion in both MYB genes was positively associated with ASE and anthocyanin accumulation in apple petals through analysis of 231 apple accessions, and thus underpins apple flower colour evolution. Our study demonstrated the importance of TEs in regulating ASE on a genome-wide scale and presents a novel method for rapid identification of ASE genes and their regulatory elements in plants.

microRNA172 targets APETALA2 to regulate flavonoid biosynthesis in apple (Malus domestica)

MicroRNA172 (miR172) plays a role in regulating a diverse range of plant developmental processes, including flowering, fruit development and nodulation. However, its role in regulating flavonoid biosynthesis is unclear. In this study, we show that transgenic apple plants over-expressing miR172 show a reduction in red coloration and anthocyanin accumulation in various tissue types. This reduction was consistent with decreased expression of APETALA2 homolog MdAP2_1a (a miR172 target gene), MdMYB10, and targets of MdMYB10, as demonstrated by both RNA-seq and qRT-PCR analyses. The positive role of MdAP2_1a in regulating anthocyanin biosynthesis was supported by the enhanced petal anthocyanin accumulation in transgenic tobacco plants overexpressing MdAP2_1a, and by the reduction in anthocyanin accumulation in apple and cherry fruits transfected with an MdAP2_1a virus-induced-gene-silencing construct. We demonstrated that MdAP2_1a could bind directly to the promoter and protein sequences of MdMYB10 in yeast and tobacco, and enhance MdMYB10 promotor activity. In Arabidopsis, over-expression of miR172 reduced flavonoid (including anthocyanins and flavonols) concentration and RNA transcript abundance of flavonoid genes in plantlets cultured on medium containing 7% sucrose. The anthocyanin content and RNA abundance of anthocyanin genes could be partially restored by using a synonymous mutant of MdAP2_1a, which had lost the miR172 target sequences at mRNA level, but not restored by using a WT MdAP2_1a. These results indicate that miR172 inhibits flavonoid biosynthesis through suppressing the expression of an AP2 transcription factor that positively regulates MdMYB10.

A Scalable Strand-Specific Protocol Enabling Full-Length Total RNA Sequencing From Single Cells

RNA sequencing (RNAseq) has been widely used to generate bulk gene expression measurements collected from pools of cells. Only relatively recently have single-cell RNAseq (scRNAseq) methods provided opportunities for gene expression analyses at the single-cell level, allowing researchers to study heterogeneous mixtures of cells at unprecedented resolution. Tumors tend to be composed of heterogeneous cellular mixtures and are frequently the subjects of such analyses. Extensive method developments have led to several protocols for scRNAseq but, owing to the small amounts of RNA in single cells, technical constraints have required compromises. For example, the majority of scRNAseq methods are limited to sequencing only the 3' or 5' termini of transcripts. Other protocols that facilitate full-length transcript profiling tend to capture only polyadenylated mRNAs and are generally limited to processing only 96 cells at a time. Here, we address these limitations and present a novel protocol that allows for the high-throughput sequencing of full-length, total RNA at single-cell resolution. We demonstrate that our method produced strand-specific sequencing data for both polyadenylated and non-polyadenylated transcripts, enabled the profiling of transcript regions beyond only transcript termini, and yielded data rich enough to allow identification of cell types from heterogeneous biological samples.

Transcriptome of Edwardsiella anguillarum in vivo and in vitro revealed two-component system, ABC transporter and flagellar assembly are three pathways pathogenic to European eel (Anguilla anguilla)

Edwardsiella anguillarum is one of the common bacterial pathogens for the cultivated eels in China. The aim of this study was to reveal the cause of E. anguillarum pathogenic to European eel (Anguilla anguilla) from the perspective of the transcriptome. In this study, we first prepared E. anguillarum cultured in vitro and analysed the whole transcriptome after extracting the total RNA. Then, eels were i.p injected with E. anguillarum, and total RNA were extracted from the liver of European eels 48 h after the infection. After sequencing the transcriptome, we obtained average 1.97 × 10⁸ clean reads cultured in vitro and 1.36 × 10⁵ clean reads located in vivo after annotating all reads into the genome of E. anguillarum. The whole transcriptome showed, compared to the E. anguillarum cultured in vitro, 503 significantly up and 657 significantly down-regulated different expressed genes (DEGs) were observed. KEGG analysis showed that 38 DEGs of Two-Component System, 41 DEGs of ABC transporter, and 10 DEGs flagellar assembly pathways were highly upregulated in E. anguillarum located in vivo. Then, we designed primers to analyse the up-regulated DEGs through qRT-PCR and confirmed some up-regulated DEGs. The results of this study provide important reference for the further study of pathogen-host interaction between E. anguillarum and European eel.

Identifying inaccuracies in gene expression estimates from unstranded RNA-seq data

RNA-seq methods are widely utilized for transcriptomic profiling of biological samples. However, there are known caveats of this technology which can skew the gene expression estimates. Specifically, if the library preparation protocol does not retain RNA strand information then some genes can be erroneously quantitated. Although strand-specific protocols have been established, a significant portion of RNA-seq data is generated in non-strand-specific manner. We used a comprehensive stranded RNA-seq dataset of 15 blood cell types to identify genes for which expression would be erroneously estimated if strand information was not available. We found that about 10% of all genes and 2.5% of protein coding genes have a two-fold or higher difference in estimated expression when strand information of the reads was ignored. We used parameters of read alignments of these genes to construct a machine learning model that can identify which genes in an unstranded dataset might have incorrect expression estimates and which ones do not. We also show that differential expression analysis of genes with biased expression estimates in unstranded read data can be recovered by limiting the reads considered to those which span exonic boundaries. The resulting approach is implemented as a package available at https://github.com/mikpom/uslcount.

Evaluation of protocols for rRNA depletion-based RNA sequencing of nanogram inputs of mammalian total RNA

Next generation RNA-sequencing (RNA-seq) is a flexible approach that can be applied to a range of applications including global quantification of transcript expression, the characterization of RNA structure such as splicing patterns and profiling of expressed mutations. Many RNA-seq protocols require up to microgram levels of total RNA input amounts to generate high quality data, and thus remain impractical for the limited starting material amounts typically obtained from rare cell populations, such as those from early developmental stages or from laser micro-dissected clinical samples. Here, we present an assessment of the contemporary ribosomal RNA depletion-based protocols, and identify those that are suitable for inputs as low as 1-10 ng of intact total RNA and 100-500 ng of partially degraded RNA from formalin-fixed paraffin-embedded tissues.

Transcriptome Analysis of Differentially Expressed Genes Related to the Growth and Development of the Jinghai Yellow Chicken

The growth traits are important traits in chickens. Compared to white feather broiler breeds, Chinese local broiler breeds have a slow growth rate. The main genes affecting the growth traits of local chickens in China are still unclear and need to be further explored. This experiment used fast-growth and slow-growth groups of the Jinghai Yellow chicken as the research objects. Three males and three females with similar body weights were selected from the two groups at four weeks old and eight weeks old, respectively, with a total of 24 individuals selected. After slaughter, their chest muscles were taken for transcriptome sequencing. In the differentially expressed genes screening, all of the genes obtained were screened by fold change ≥ 2 and false discovery rate (FDR) < 0.05. For four-week-old chickens, a total of 172 differentially expressed genes were screened in males, where there were 68 upregulated genes and 104 downregulated genes in the fast-growth group when compared with the slow-growth group. A total of 31 differentially expressed genes were screened in females, where there were 11 upregulated genes and 20 downregulated genes in the fast-growth group when compared with the slow-growth group. For eight-week-old chickens, a total of 37 differentially expressed genes were screened in males. The fast-growth group had 28 upregulated genes and 9 downregulated genes when compared with the slow-growth group. A total of 44 differentially expressed genes were screened in females. The fast-growth group had 13 upregulated genes and 31 downregulated genes when compared with the slow-growth group. Through gene ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, many genes were found to be related to cell proliferation and differentiation, muscle growth, and cell division such as SNCG, MCL1, ARNTL, PLPPR4, VAMP1, etc. Real-time PCR results were consistent with the RNA-Seq data and validated the findings. The results of this study will help to understand the regulation mechanism of the growth and development of Jinghai Yellow chicken and provide a theoretical basis for improving the growth rate of Chinese local chicken breeds.