RNA sequencing reveals two major classes of gene expression levels in metazoan cells

The expression level of a gene is often used as a proxy for determining whether the protein or RNA product is functional in a cell or tissue. Therefore, it is of fundamental importance to understand the global distribution of gene expression levels, and to be able to interpret it mechanistically and functionally. Here we use RNA sequencing (RNA‐seq) of mouse Th2 cells, coupled with a range of other techniques, to show that all genes can be separated, based on their expression abundance, into two distinct groups: one group comprised of lowly expressed and putatively non‐functional mRNAs, and the other of highly expressed mRNAs with active chromatin marks at their promoters. These observations are confirmed in many other microarray and RNA‐seq data sets of metazoan cell types.     

Comparative transcriptional profiling analysis of the effect of heat waves during embryo incubation on the hatchlings of the Chinese soft‐shelled turtle (Pelodiscus sinensis)

Temperature is one of most the important environmental factors that affect the ontogenesis of organisms. In this study, we incubated Chinese soft‐shelled turtle eggs at 28°C (control temperature, C treatment), a temperature with a 16°C cold shock and a 36°C heat shock twice per week (S treatment) or a ramp‐programmed temperature of 29 ± 9°C (with 12 hr (+) and 12 hr (?) every day) (F treatment). The incubation period, hatching success, hatchling weight, and locomotor performance were significantly different between the controls and the different heat treatment groups. The pathogen challenge results illustrated that hatchlings from the S treatment group were more resistant to bacterial infection, whereas hatchlings from the F treatment group were more vulnerable. We used RNA‐seq quantification analysis to identify differentially expressed genes (DEGs) of hatchlings in the S treatment group. Based on the functional annotation results for the DEGs, 9 genes were chosen to verify the RNA‐seq results. The background expression of DEGs was also analyzed for the three treatments, as was the regulation of the pathogen challenge. The results showed that 8 DEGs were related to the immune response after pathogen challenge and that temperature was an important factor in differential regulation of the immunity pathways.     

De novo assembly of a chromosome‐level reference genome of red‐spotted grouper (Epinephelus akaara) using nanopore sequencing and Hi‐C

The red‐spotted grouper Epinephelus akaara (E. akaara) is one of the most economically important marine fish in China, Japan and South‐East Asia and is a threatened species. The species is also considered a good model for studies of sex inversion, development, genetic diversity and immunity. Despite its importance, molecular resources for E. akaara remain limited and no reference genome has been published to date. In this study, we constructed a chromosome‐level reference genome of E. akaara by taking advantage of long‐read single‐molecule sequencing and de novo assembly by Oxford Nanopore Technology (ONT) and Hi‐C. A red‐spotted grouper genome of 1.135 Gb was assembled from a total of 106.29 Gb polished Nanopore sequence (GridION, ONT), equivalent to 96‐fold genome coverage. The assembled genome represents 96.8% completeness (BUSCO) with a contig N50 length of 5.25 Mb and a longest contig of 25.75 Mb. The contigs were clustered and ordered onto 24 pseudochromosomes covering approximately 95.55% of the genome assembly with Hi‐C data, with a scaffold N50 length of 46.03 Mb. The genome contained 43.02% repeat sequences and 5,480 noncoding RNAs. Furthermore, combined with several RNA‐seq data sets, 23,808 (99.5%) genes were functionally annotated from a total of 23,923 predicted protein‐coding sequences. The high‐quality chromosome‐level reference genome of E. akaara was assembled for the first time and will be a valuable resource for molecular breeding and functional genomics studies of red‐spotted grouper in the future.     

RNA‐sequencing profiles hippocampal gene expression in a validated model of cancer‐induced depression

To investigate the pathophysiology of cancer‐induced depression (CID), we have recently developed a validated CID mouse model. Given that the efficacy of antidepressants in cancer patients is controversial, it remains unclear whether CID is a biologically distinct form of depression. We used RNA‐sequencing (RNA‐seq) to investigate differentially expressed genes (DEGs) in hippocampi of animals from our CID model relative a positive control model of depressive‐like behavior induced with chronic corticosterone (CORT). To validate RNA‐seq results, we performed quantitative real‐time RT‐PCR (qRT‐PCR) on a subset of DEGs. Enrichment analysis using DAVID was performed on DEGs to identify enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and biological process gene ontologies (GO:BP). qRT‐PCR results significantly predicted RNA‐seq results. RNA‐seq revealed that most DEGs identified in the CORT model overlapped with the CID model. Enrichment analyses identified KEGG pathways and GO:BP terms associated with ion homeostasis and neuronal communication for both the CORT and CID model. In addition, CID DEGs were enriched in pathways and terms relating to neuronal development, intracellular signaling, learning and memory. This study is the first to investigate CID at the mRNA level. We have shown that most hippocampal mRNA changes that are associated with a depressive‐like state are also associated with cancer. Several other changes occur at the mRNA level in cancer, suggesting that the CID model may represent a biologically distinct form of a depressive‐like state.     

Single cell RNA‐seq in the sea urchin embryo show marked cell‐type specificity in the Delta/Notch pathway

Sea urchin embryos are excellent for in vivo functional studies because of their transparency and tractability in manipulation. They are also favorites for pharmacological approaches since they develop in an aquatic environment and addition of test substances is straightforward. A concern in many pharmacological tests though is the potential for pleiotropic effects that confound the conclusions drawn from the results. Precise cellular interpretations are often not feasible because the impact of the perturbant is not known. Here we use single‐cell mRNA (messenger RNA) sequencing as a metric of cell types in the embryo and to determine the selectivity of two commonly used inhibitors, one each for the Wnt and the Delta‐Notch pathways, on these nascent cell types. We identified 11 distinct cell types based on mRNA profiling, and that the cell lineages affected by Wnt and Delta/Notch inhibition were distinct from each other. These data support specificity and distinct effects of these signaling pathways in the embryo and illuminate how these conserved pathways selectively regulate cell lineages at a single cell level. Overall, we conclude that single cell RNA‐seq analysis in this embryo is revealing of the cell types present during development, of the changes in the gene regulatory network resulting from inhibition of various signaling pathways, and of the selectivity of these pathways in influencing developmental trajectories.     

Whole‐transcriptome analysis reveals genetic factors underlying flowering time regulation in rapeseed (Brassica napus L.)

Rapeseed (Brassica napus L.), one of the most important sources of vegetable oil and protein‐rich meals worldwide, is adapted to different geographical regions by modification of flowering time. Rapeseed cultivars have different day length and vernalization requirements, which categorize them into winter, spring, and semiwinter ecotypes. To gain a deeper insight into genetic factors controlling floral transition in B. napus, we performed RNA sequencing (RNA‐seq) in the semiwinter doubled haploid line, Ningyou7, at different developmental stages and temperature regimes. The expression profiles of more than 54,000 gene models were compared between different treatments and developmental stages, and the differentially expressed genes were considered as targets for association analysis and genetic mapping to confirm their role in floral transition. Consequently, 36 genes with association to flowering time, seed yield, or both were identified. We found novel indications for neofunctionalization in homologs of known flowering time regulators like VIN3 and FUL. Our study proved the potential of RNA‐seq along with association analysis and genetic mapping to identify candidate genes for floral transition in rapeseed. The candidate genes identified in this study could be subjected to genetic modification or targeted mutagenesis and genotype building to breed rapeseed adapted to certain environments.