Full-length Transcriptome Sequencing Research Articles

Umbilical Cord Mesenchymal Stem Cell-Derived Small Extracellular Vesicles Deliver miR-21 to Promote Corneal Epithelial Wound Healing through PTEN/PI3K/Akt Pathway.

Objective Rapid restoration of corneal epithelium integrity after injury is particularly important for preserving corneal transparency and vision. Mesenchymal stem cells (MSCs) can be taken into account as the promising regenerative therapeutics for improvement of wound healing processes based on the variety of the effective components. The extracellular vesicles form MSCs, especially exosomes, have been considered as important paracrine mediators though transferring microRNAs into recipient cell. This study investigated the mechanism of human umbilical cord MSC-derived small extracellular vesicles (HUMSC-sEVs) on corneal epithelial wound healing. Methods HUMSC-sEVs were identified by transmission electron microscopy, nanoparticle tracking analysis, and Western blot. Corneal fluorescein staining and histological staining were evaluated in a corneal mechanical wound model. Changes in HCEC proliferation after HUMSC-sEVs or miR-21 mimic treatment were evaluated by CCK-8 and EdU assays, while migration was assessed by in vitro scratch wound assay. Full-length transcriptome sequencing was performed to identify the differentially expressed genes associated with HUMSC-sEVs treatment, followed by validation via real-time PCR and Western blot. Results The sEVs derived from HUMSCs can significantly promote corneal epithelial cell proliferation, migration in vitro, and corneal epithelial wound healing in vivo. Similar effects were obtained after miR-21 transfection, while the beneficial effects of HUMSC-sEVs were partially negated by miR-21 knockdown. Results also show that the benefits are associated with decreased PTEN level and activated the PI3K/Akt signaling pathway in HCECs. Conclusion HUMSC-sEVs could enhance the recovery of corneal epithelial wounds though restraining PTEN by transferring miR-21 and may represent a promising novel therapeutic agent for corneal wound repair.

Full-Length Transcriptome Sequencing Reveals Alternative Splicing and lncRNA Regulation during Nodule Development in Glycine max.

Alternative splicing (AS) is a ubiquitous phenomenon among eukaryotic intron-containing genes, which greatly contributes to transcriptome and proteome diversity. Here we performed the isoform sequencing (Iso-Seq) of soybean underground tissues inoculated and uninoculated with Rhizobium and obtained 200,681 full-length transcripts covering 26,183 gene loci. It was found that 80.78% of the multi-exon loci produced more than one splicing variant. Comprehensive analysis of these identified 7874 differentially splicing events with highly diverse splicing patterns during nodule development, especially in defense and transport-related processes. We further profiled genes with differential isoform usage and revealed that 2008 multi-isoform loci underwent stage-specific or simultaneous major isoform switches after Rhizobium inoculation, indicating that AS is a vital way to regulate nodule development. Moreover, we took the lead in identifying 1563 high-confidence long non-coding RNAs (lncRNAs) in soybean, and 157 of them are differentially expressed during nodule development. Therefore, our study uncovers the landscape of AS during the soybean-Rhizobium interaction and provides systematic transcriptomic data for future study of multiple novel directions in soybean.

Evolutionary Analysis and Functional Identification of Ancient Brassinosteroid Receptors in Ceratopteris richardii

Phytohormones play an important role in the adaptive evolution of terrestrial plants. Brassinosteroids (BRs) are essential hormones that regulate multiple aspects of plant growth and development in angiosperms, but the presence of BR signaling in non-seed plants such as ferns remains unknown. Here, we found that BR promotes the growth of Ceratopteris richardii, while the synthetic inhibitor PCZ inhibits the growth. Using full-length transcriptome sequencing, we identified four BRI1-like receptors. By constructing chimeric receptors, we found that the kinase domains of these four receptors could trigger BR downstream signaling. Further, the extracellular domains of two receptors were functionally interchangeable with that of BRI1. In addition, we identified a co-receptor, CtSERK1, that could phosphorylate with CtBRL2s in vitro. Together, these proved the presence of a receptor complex in Ceratopteris richardii that might perceive BR and activate downstream hormone signaling. Our results shed light on the biological and molecular mechanisms of BR signaling in ferns and the role of BR hormone signaling in the adaptive evolution of terrestrial plants.

Full-Length Transcriptomics Reveals Complex Molecular Mechanism of Salt Tolerance in Bromus inermis L.

Bromus inermis L. (commonly known as smooth bromegrass) is a grass species with high nutritional value, great palatability, cold tolerance, and grazing resistance, which has been widely cultivated for pasture and sand fixation in northern and northwestern China. Salt stress is a main environmental factor limiting growth and production of smooth bromegrass. In this study, we performed PacBio Iso-Seq to construct the first full-length transcriptome database for smooth bromegrass under 300 mM NaCl treatment at different time points. Third-generation full-length transcriptome sequencing yielded 19.67 G polymerase read bases, which were assembled into 355,836 full-length transcripts with an average length of 2,542 bp. A total of 116,578 differentially expressed genes were obtained by comparing the results of third-generation sequencing and second-generation sequencing. GO and KEGG enrichment analyses revealed that multiple pathways were differently activated in leaves and roots. In particular, a number of genes participating in the molecular network of plant signal perception, signal transduction, transcription regulation, antioxidant defense, and ion regulation were affected by NaCl treatment. In particular, the CBL-CIPK, MAPK, ABA signaling network, and SOS core regulatory pathways of Ca2+ signal transduction were activated to regulate salt stress response. In addition, the expression patterns of 10 salt-responsive genes were validated by quantitative real-time PCR, which were consistent with those detected by RNA-Seq. Our results reveal the molecular regulation of smooth bromegrass in response to salt stress, which are important for further investigation of critical salt responsive genes and molecular breeding of salt-tolerant smooth bromegrass.

Integrated SMRT and Illumina Sequencing Provide New Insights into Crocin Biosynthesis of Gardenia jasminoides.

Crocins are valuable bioactive components of gardenia fruit, and their biosynthesis and accumulation have attracted widespread interest. Studies have investigated the biosynthesis and accumulation of crocin based on Illumina sequencing, but there is a lack of reports based on full-length transcriptome sequencing. Utilising SMRT sequencing and high-performance liquid chromatography (HPLC), we explored crocin biosynthesis and accumulation in the fruit of Gardenia jasminoides. HPLC analysis showed that crocins specifically exist in fruit and that the content of crocins increases gradually during fruit development. SMRT sequencing generated 46,715 high-quality full-length isoforms, including 5230 novel isoforms that are not present in the G. jasminoides genome. Furthermore, a total of 46 genes and 91 lncRNAs were involved in the biosynthesis and accumulation of crocin. The qRT-PCR indicated that genes involved in crocin biosynthesis reached a peak in the NOV stage. These findings contributed to our understanding of crocin biosynthesis and accumulation.

Full-Length Transcriptome Sequencing Analysis of Differentially Expressed Genes and Pathways After Treatment of Psoriasis With Oxymatrine.

Psoriasis is a recurrent chronic inflammatory skin disease. Unlike many of the latest psoriasis treatments that only confer limited curative effects and have certain side effects, oxymatrine effectively improves severe plaque psoriasis with mild adverse reactions. Here, we explored the genes and pathways underlying the effects of oxymatrine on psoriasis. Briefly, patients with severe plaque psoriasis were treated with oxymatrine and their lesioned skin samples were sequenced by full-length transcriptomics. Next, the differentially expressed genes (DEGs) in psoriatic lesions were identified and compared in oxymatrine-treated patients and healthy controls, their genes were functionally annotated, and protein–protein interaction network analysis and immunohistochemistry were performed. Both Psoriasis Area and Severity Index (PASI) and Body Surface Area (BSA) scores were recovered significantly from all 16 patients (all p < 0.001). The number of DEGs in patients before and after oxymatrine treatment was 4232, and 4105 DEGs were found between the psoriasis group (before oxymatrine treatment) and the normal control group [p < 0.01, |log2 fold change, (FC)| >1.5]. While most of the DEGs recovered significantly after oxymatrine treatment, only 650 DEGs were observed between the psoriasis group (after oxymatrine treatment) and the normal control group (p < 0.01, |log2FC|> 1.5). Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed that 64 pathways were significantly activated after oxymatrine treatment (p < 0.05). Only 12 pathways were statistically significant between after oxymatrine treatment and the normal control group (p < 0 .05). Among all the restored pathways, the improvement of the IL-17 signaling pathway was the most significant (p = 1.18E-06). Gene loci of oxymatrine action was assessed by protein interaction analysis on 205 DEGs that were co-expressed in 5 patients before and after oxymatrine treatment (p < 0.05, FC > 1.5). After oxymatrine treatment, the expression of two mitosis-related genes namely, cyclin dependent kinase 1 (CDK1) and cyclin B1 (CCNB1), that affect cell proliferation recovered significantly. In light of these results, we conclude that oxymatrine likely alters the abnormal expression of some genes and pathways in psoriasis patients. Multipathway and multitarget therapy can greatly ameliorate abnormalities in genes and pathways and effectively treat psoriasis. Importantly, among the DEGs, the proliferation-related genes, such as CDK1 and CCNB1, are likely important targets for treating psoriasis by oxymatrine. We believe that these findings may lead to a new treatment strategy for psoriasis.

The Isolation and Full-Length Transcriptome Sequencing of a Novel Nidovirus and Response of Its Infection in Japanese Flounder (Paralichthys olivaceus).

A novel nidovirus, CSBV Bces-Po19, was isolated from the marine fish, Japanese flounder (Paralichthys olivaceus). The viral genome was 26,597 nucleotides long and shared 98.62% nucleotide identity with CSBV WHQSR4345. PacBio Sequel and Illumina sequencing were used to perform full-length transcriptome sequencing on CSBV Bces-Po19-sensitive (S) and -resistant (R) Japanese flounder. The results of negative staining revealed bacilliform and spherical virions. There were in total 1444 different genes between CSBV Bces-Po19 S and R groups, with 935 being up-regulated and 513 being down-regulated. Metabolism-, immune-, and RNA-related pathways were significantly enriched. Furthermore, CSBV Bces-Po19 infection induced alternative splicing (AS) events in Japanese flounder; the S group had a higher numbers of AS events (12,352) than the R group (11,452). The number of long non-coding RNA (lncRNA) in the S group, on the other hand, was significantly lower than in the R group. In addition to providing valuable information that sheds more light on CSBV Bces-Po19 infection, these research findings provide further clues for CSBV Bces-Po19 prevention and treatment.

Iso-Seq assembly and functional annotation of full-length transcriptome of turbot (Scophthalmus maximus) during bacterial infection.

Turbot (Scophthalmus maximus) is a flatfish, which is not only important in mariculture worldwide with the unique characteristic of body asymmetry, but also as an economically important species in aquaculture. Herein, we performed the first full-length transcriptome sequencing of turbot during the bacterial infection. A total of 307.1 Gb raw reads were obtained and processed with Iso-Seq, generating 187,509 high-quality redundant transcripts with an average length of 3005 base pairs. Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis identified 81.5% complete BUSCOs and only 1.7% fragmented BUSCOs, suggesting a transcript structural completeness and functional diversity of this transcriptome. Moreover, the redundant transcripts were collapsed and compared to ENSEMBL reference with Cupcake and SQANTI3. Among 60,476 collapsed transcripts, we identified 12,059 annotated and 1684 novel genes. 42,956 (71.1%) transcripts provided new evidence for splice junctions identification. Furthermore, the untranslated region (UTR) identification was also benefited from the transcriptome. The open read frames prediction was conducted with PASApipeline. 42,118 transcripts were assigned with known function by aligning against Swiss-Prot or functional domain prediction. Taken together, the full-length transcriptome built in this study could provide important resources for immunologic research on turbot.

A Global Regulatory Network for Dysregulated Gene Expression and Abnormal Metabolic Signaling in Immune Cells in the Microenvironment of Graves' Disease and Hashimoto's Thyroiditis.

BackgroundAlthough the pathogenetic mechanisms of Hashimoto’s thyroiditis (HT) and Graves’ disease (GD) have been elucidated, the molecular mechanisms by which the abnormal immune function of cellular subpopulations trigger an autoimmune attack on thyroid tissue largely remains unexplained.MethodsThe study included 2 HT patients, 2 GD patients, and 1 control donor. The thyroid samples were extracted for single-cell RNA sequencing, whole transcriptome, full-length transcriptome (Oxford Nanopore Technologies), and metabolome sequencing. Identification of immune cells with dysregulated gene expression and abnormal metabolic signaling was performed in the microenvironment, both at the bulk and single-cell levels. Based on functional enrichment analysis, the biological processes and pathways involved in abnormal immune cells were further explored. Finally, according to cell communication analysis, the global regulatory network of immune cells was constructed.ResultsCD4+ T cells, CD8+ T cells, and macrophages were abnormally increased in patients with HT and GD. The differentially expressed genes of these cells were significantly involved in signaling pathways, including Th1 and Th2 cell differentiation, Th17 cell differentiation, cytokine–cytokine receptor interaction, and NF-kappa B signaling pathway. Moreover, in HT, CD4+ T cells interact with macrophages via the IL16-CCR5/FGF10-FGFR1/CXCL13-CXCR3 axis, and macrophages interact with CD8+ T cells via the CD70-CD27 axis, thereby activating the T-cell receptor signaling pathway and NF-kappa B signaling pathway. In GD, CD4+ T cells interact with macrophages via the CXCR3-CXCL10/PKM-CD44/MHCII-NFKBIE axis, and macrophages interact with CD8+ T cells via the IFNG-IFNGR1/CCR7-CCL21 axis, thereby activating T-cell receptor signaling pathway, Th1 and Th2 cell differentiation, and chemokine signaling pathway.ConclusionIn HT and GD, immune dysregulated cells interact and activate relevant immune pathways and further aggravate the immune response. This may trigger the immune cells to target the thyroid tissue and influence the development of the disease.

Full-Length Transcriptome Sequencing Analysis and Characterization of Gene Isoforms Involved in Flavonoid Biosynthesis in the Seedless Kiwifruit Cultivar ‘Chengxiang’ (Actinidia arguta)

Kiwifruit an important horticultural crop that is widely cultivated and is known as the king of fruits. Recently, a new seedless kiwifruit cultivar, ‘Chengxiang’ (Actinidia arguta), was discovered by field transplantation. It exhibited distinguishable characteristics such as parthenocarpy, and a unique flavor and appearance when compared to other cultivated type. Flavonoids are known to play an important role in fertility and parthenocarpy in plants. However, the genes responsible for flavonoid biosynthesis in seedless kiwifruit remain largely unknown. Especially, chalcone synthase (CHS), as a key enzyme catalyzing the first committed step in the flavonoid pathway, remains a mystery. In this study, we combined a full-length transcriptome survey by PacBio single-molecule real-time (SMRT) sequencing, CHS gene family analysis, and analysis of the gene expression involved in flavonoid pathways to further enhance the understanding of parthenocarpy. Based on SMRT, we obtained 80,615 high-quality full-length consensus transcripts. In total, 52,406 (90.79%) transcripts were functionally annotated, and more than 80% of the transcripts were longer than 1Kb. Among them, 39,117 (74.64%) transcripts were assigned to GO terms, the majority of which were associated with the cell (19,089, 48.80%) and metabolic process (19,859, 50.77%). Furthermore, 25,289 (48.26%) transcripts were mapped into 129 KEGG pathways. We identified the majority of putative genes as being involved in the flavonoid biosynthesis pathway, including 14 key enzyme gene families, such as CHS, chalcone isomerase (CHI), flavonol synthase (FLS), and so on. Moreover, we also identified 13 CHS genes and characterized the CHS gene family in seedless kiwifruit. We further evaluated the expression pattern of 10 flavonoid-related key enzyme genes in flowers using quantitative real-time PCR. This is the first time that the full-length transcriptome have been studied in seedless kiwifruit, and the findings enhance our understanding the molecular mechanisms of parthenocarpy.

Full-Length Transcriptome Sequencing Reveals the Impact of Cold Stress on Alternative Splicing in Quinoa.

Quinoa is a cold-resistant and nutrient-rich crop. To decipher the cold stress response of quinoa, the full-length transcriptomes of the cold-resistant quinoa variety CRQ64 and the cold-sensitive quinoa variety CSQ5 were compared. We identified 55,389 novel isoforms and 6432 novel genes in these transcriptomes. Under cold stress, CRQ64 had more differentially expressed genes (DEGs) and differentially alternative splicing events compared to non-stress conditions than CSQ5. DEGs that were specifically present only in CRQ64 were significantly enriched in processes which contribute to osmoregulation and ROS homeostasis in plants, such as sucrose metabolism and phenylpropanoid biosynthesis. More genes with differential alternative splicing under cold stress were enriched in peroxidase functions in CRQ64. In total, 5988 transcription factors and 2956 long non-coding RNAs (LncRNAs) were detected in this dataset. Many of these had altered expression patterns under cold stress compared to non-stress conditions. Our transcriptome results demonstrate that CRQ64 undergoes a wider stress response than CSQ5 under cold stress. Our results improved the annotation of the quinoa genome and provide new insight into the mechanisms of cold resistance in quinoa.

Full-Length Transcriptional Analysis of the Same Soybean Genotype With Compatible and Incompatible Reactions to Heterodera glycines Reveals Nematode Infection Activating Plant Defense Response.

Full-length transcriptome sequencing with long reads is a powerful tool to analyze transcriptional and post-transcriptional events; however, it has not been applied on soybean (Glycine max). Here, a comparative full-length transcriptome analysis was performed on soybean genotype 09-138 infected with soybean cyst nematode (SCN, Heterodera glycines) race 4 (SCN4, incompatible reaction) and race 5 (SCN5, compatible reaction) using Oxford Nanopore Technology. Each of 9 full-length samples collected 8 days post inoculation with/without nematodes generated an average of 6.1 GB of clean data and a total of 65,038 transcript sequences. After redundant transcripts were removed, 1,117 novel genes and 41,096 novel transcripts were identified. By analyzing the sequence structure of the novel transcripts, a total of 28,759 complete open reading frame (ORF) sequences, 5,337 transcription factors, 288 long non-coding RNAs, and 40,090 novel transcripts with function annotation were predicted. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of differentially expressed genes (DEGs) revealed that growth hormone, auxin-activated signaling pathway and multidimensional cell growth, and phenylpropanoid biosynthesis pathway were enriched by infection with both nematode races. More DEGs associated with stress response elements, plant-hormone signaling transduction pathway, and plant–pathogen interaction pathway with more upregulation were found in the incompatible reaction with SCN4 infection, and more DEGs with more upregulation involved in cell wall modification and carbohydrate bioprocess were detected in the compatible reaction with SCN5 infection when compared with each other. Among them, overlapping DEGs with a quantitative difference was triggered. The combination of protein–protein interaction with DEGs for the first time indicated that nematode infection activated the interactions between transcription factor WRKY and VQ (valine-glutamine motif) to contribute to soybean defense. The knowledge of the SCN–soybean interaction mechanism as a model will present more understanding of other plant–nematode interactions.

Combination analysis of single-molecule long-read and Illumina sequencing provides insights into the anthocyanin accumulation mechanism in an ornamental grass, Pennisetum setaceum cv. Rubrum

Combination analysis of single-molecule long-read and Illumina sequencing provide full-length transcriptome information and shed new light on the anthocyanin accumulation mechanism of Pennisetum setaceum cv. 'Rubrum'. Pennisetum setaceum cv. 'Rubrum' is an ornamental grass with purple leaves widely used in landscaping. However, the current next-generation sequencing (NGS) transcriptome information of this species is not satisfactory due to the difficulties in obtaining full-length transcripts. Furthermore, the molecular mechanisms of anthocyanin accumulation in P. setaceum have not been thoroughly studied. In this study, we used PacBio full-length transcriptome sequencing (SMRT) combined with NGS technology to build and improve the transcriptomic datasets and reveal the molecular mechanism of anthocyanin accumulation in P. setaceum cv. 'Rubrum'. Therefore, 280,413 full-length non-chimeric reads sequences were obtained using the SMRT technology. We obtained 97,450 high-quality non-redundant transcripts and identified 5352 alternative splicing events. In addition, 93,066 open reading frames (ORFs), including 57,457 full ORFs and 2910 long non-coding RNA (lncRNAs) were screened out. Furthermore, 10,795 differentially expressed genes were identified using NGS. We also explored key genes, synthesis pathways, and detected lncRNA involved in anthocyanin accumulation, providing new insights into anthocyanin accumulation in P. setaceum cv. 'Rubrum'. To our best knowledge, we provided the full-length transcriptome information of P. setaceum cv. 'Rubrum' for the first time. The results of this study will provide baseline information for gene function studies and pave the way for future P. setaceum cv. 'Rubrum' breeding projects.

Full-length transcriptome sequencing reveals the molecular mechanism of potato seedlings responding to low-temperature

BackgroundPotato (Solanum tuberosum L.) is one of the world's most important crops, the cultivated potato is frost-sensitive, and low-temperature severely influences potato production. However, the mechanism by which potato responds to low-temperature stress is unclear. In this research, we apply a combination of second-generation sequencing and third-generation sequencing technologies to sequence full-length transcriptomes in low-temperature-sensitive cultivars to identify the important genes and main pathways related to low-temperature resistance.ResultsIn this study, we obtained 41,016 high-quality transcripts, which included 15,189 putative new transcripts. Amongst them, we identified 11,665 open reading frames, 6085 simple sequence repeats out of the potato dataset. We used public available genomic contigs to analyze the gene features, simple sequence repeat, and alternative splicing event of 24,658 non-redundant transcript sequences, predicted the coding sequence and identified the alternative polyadenylation. We performed cluster analysis, GO, and KEGG functional analysis of 4518 genes that were differentially expressed between the different low-temperature treatments. We examined 36 transcription factor families and identified 542 transcription factors in the differentially expressed genes, and 64 transcription factors were found in the AP2 transcription factor family which was the most. We measured the malondialdehyde, soluble sugar, and proline contents and the expression genes changed associated with low temperature resistance in the low-temperature treated leaves. We also tentatively speculate that StLPIN10369.5 and StCDPK16 may play a central coordinating role in the response of potatoes to low temperature stress.ConclusionsOverall, this study provided the first large-scale full-length transcriptome sequencing of potato and will facilitate structure–function genetic and comparative genomics studies of this important crop.

Full-length transcriptome analysis of Adiantum flabellulatum gametophyte.

Ferns are important components of plant communities on earth, but their genomes are generally very large, with many redundant genes, making whole genome sequencing of ferns prohibitively expensive and time-consuming. This means there is a significant lack of fern reference genomes, making molecular biology research difficult. The gametophytes of ferns can survive independently, are responsible for sexual reproduction and the feeding of young sporophytes, and play an important role in the alternation of generations. For this study, we selected Adiantum flabellulatum as it has both ornamental and medicinal value and is also an indicator plant of acidic soil. The full-length transcriptome sequencing of its gametophytes was carried out using PacBio three-generation sequencing technology. A total of 354,228 transcripts were obtained, and 231,705 coding sequences (CDSs) were predicted, including 5,749 transcription factors (TFs), 2,214 transcription regulators (TRs) and 4,950 protein kinases (PKs). The transcripts annotated by non-redundant protein sequence database (NR), Kyoto encyclopedia of genes and genomes (KEGG), eukaryotic ortholog groups (KOG), Swissprot, protein family (Pfma), nucleotide sequence database (NT) and gene ontology (GO) were 251,501, 197,474, 193,630, 194,639, 195,956, 113,069 and 197,883, respectively. In addition, 138,995 simple sequence repeats (SSRs) and 111,793 long non-coding RNAs (lncRNAs) were obtained. We selected nine chlorophyll synthase genes for qRT-PCR, and the results showed that the full-length transcript sequences and the annotation information were reliable. This study can provide a reference gene set for subsequent gene expression quantification.

Seed Germination Mechanism of Carex rigescens Under Variable Temperature Determinded Using Integrated Single-Molecule Long-Read and Illumina Sequence Analysis.

The turfgrass species Carex rigescens has broad development and utilization prospects in landscaping construction. However, seed dormancy and a low germination rate have inhibited its application. Furthermore, the molecular mechanisms of seed germination in C. rigescens have not been thoroughly studied. Therefore, in the present study, PacBio full-length transcriptome sequencing combined with Illumina sequencing was employed to elucidate the germination mechanism of C. rigescens seeds under variable temperatures. In general, 156,750 full-length non-chimeric sequences, including those for 62,086 high-quality transcripts, were obtained using single-molecule long read sequencing. In total, 40,810 high-quality non-redundant, 1,675 alternative splicing, 28,393 putative coding sequences, and 1,052 long non-coding RNAs were generated. Based on the newly constructed full-length reference transcriptome, 23,147 differentially expressed genes were identified. We screened four hub genes participating in seed germination using weighted gene co-expression network analysis. Combining these results with the physiological observations, the important roles of sucrose and starch metabolic pathways in germination are further discussed. In conclusion, we report the first full-length transcriptome of C. rigescens, and investigated the physiological and transcriptional mechanisms of seed germination under variable temperatures. Our results provide valuable information for future transcriptional analyses and gene function studies of C. rigescens.

PRR-Mediated Immune Response and Intestinal Flora Profile in Soybean Meal-Induced Enteritis of Pearl Gentian Groupers, Epinephelus fuscoguttatus♀ × Epinephelus lanceolatus♂.

Pattern recognition receptors (PRRs) can recognize microbial-specific pathogen-associated molecular patterns, initiate signal cascade transduction, activate the expressions of host immunity and proinflammatory genes, and, ultimately, trigger an immune response against identified pathogens. The present study focused on two outcomes of feeding pearl gentian groupers with high levels of soybean meal (SBM): (1) growth performance and (2) the intestinal environment, including tissue structure, flora profile, and immune responses. Some 720 groupers were randomly divided into three groups (n = 4): (1) controls, fed a 50% fish meal feed (FM), (2) with 20% of the FM substituted with SBM (SBM20), and (3) 40% of the FM substituted with SBM (SBM40). The fish were fed these iso-nitrogenous and iso-lipidic diets for 10 weeks. They were kept in containers with 1 m3 of water under natural light and temperature levels. The experimental results demonstrate that the SBM diets significantly degraded growth performance and intestinal physiology. Typical enteritis characteristics and immune fluctuations appeared, as reflected by the enzyme activities of total superoxide dismutase and lysozyme, and the contents of immunoglobulin M, complement 3, and complement 4. 16SrDNA high-throughput sequencing showed that the intestinal flora was significantly affected, with the abundance of harmful bacteria, such as Vibrio and Streptococcus, increasing with dietary SBM level. Based on “3 + 2” full-length transcriptome sequencing, three triggered PRRs were found in the intestine: the RIG-like receptor, NOD-like receptor, and Toll-like receptor signaling pathways. The intestinal flora variations were significantly correlated with the activation of the three PRR signaling pathways by canonical correlation analysis. These culminated in the transcriptome activation of NF-κB, IRFs, and costimulatory molecules, ultimately promoting the expressions of proinflammatory cytokines, interferons (IFNs), chemokines, and other molecules vital to the innate and/or adaptive immune responses. This study provides new information for diagnosing and preventing SBMIE in aquaculture fish.

The full-length transcriptome sequencing of three sturgeons reveals the occurrence of whole genome duplication event.

Sturgeon is one of the ancient fish species with the fossil records over 200 million years ago. However, the sturgeon genomes are rather complicated, owing to several rounds of genome polyploidizations, which makes it difficult to assemble a high quality genome so that only a few studies analyze their whole genome evolutions (WGDs). This study applied the PacBio sequencing technology to sequence the full-length transcriptomes of blood samples of octoploid Huso dauricus, octoploid Acipenser schrenckii and tetraploid Acipenser stellatus. A total of non-redundant 33,407 transcript assemblies of the three sturgeons were finally generated and fully annotated with different gene function databases, as well as the predictions of simple sequence repeats, transcription factors and long non-coding RNAs. Data valuation exhibited that annotations of the three transcriptomes had the higher quality than those of the previously sequenced sturgeon genomes. A phylogeny of eight typical sturgeons was constructed using the low copy genes, which illustrated that Acipenser and Huso might share a single origin. The sturgeon WGDs were parallelly timed by calculating the divergence times of 2-member gene clusters and heterozygous-SNP-carried genes, which both suggested the occurrence of the WGDs later than 150 mya. The full-length transcriptome sequences were supposed to be the reliable reference data for the study of sturgeon.

Screening and Identification of the First Non-CRISPR/Cas9-Treated Chinese Miniature Pig With Defective Porcine Endogenous Retrovirus pol Genes.

Pig to human xenotransplantation is considered to be a possible approach to alleviate the shortage of human allografts. Porcine endogenous retrovirus (PERV) is the most significant pathogen in xenotransplantation. We screened for pigs that consistently did not transmit human-tropic replication competent PERVs (HTRC PERVs), namely, non-transmitting pigs. Then, we conducted whole-genome resequencing and full-length transcriptome sequencing to further investigate the sequence characteristics of one non-transmitting pig. Using in vitro transmission assays, we found 5 (out of 105) pigs of the Chinese Wuzhishan minipig inbred line that did not transmit PERV to human cells, i.e., non-transmitting pigs. Whole-genome resequencing and full-length transcriptome sequencing of one non-transmitting pig showed that all of the pol genes were defective at both the genome and transcript levels. We speculate that the defective PERV pol genes in this pig might be attributable to the long-term inbreeding process. This discovery is promising for the development of a strain of highly homozygous and genetically stable pigs with defective PERV pol genes as a source animal species for xenotransplantation.

Comparative transcriptomic analysis reveals genes related to the rapid accumulation of oleic acid in Camellia chekiangoleosa, an oil tea plant with early maturity and large fruit

Camellia chekiangoleosa has a higher oleic acid content and a shorter reproductive cycle than typical oil tea plants. It was intensively sampled over six C. chekiangoleosa seed development stages. The content of fatty acids determined by GC showed that the accumulation of fatty acids gradually increased from the S1 to S5 stages, and the maximum concentration was reached in S5. Then, fatty acids declined slightly in S6. The main fatty acid component showed the same accumulation trend as the total fatty acids, except linolenic acid, which remained at a low level throughout seed developmental stages. Changes in the expression of fatty acid accumulation-related genes were monitored using second-generation and SMRT full-length transcriptome sequencing. Finally, 18.92 G accurate and reliable data were obtained. Differential expression analysis and weighted coexpression analysis revealed two “gene modules” significantly associated with oleic acid and linoleic acid contents, and the high expression of ENR, KAS I, and KAS II, which accumulate substrates for oleic acid synthesis, was thought to be responsible for the rapid accumulation of fatty acids in the early stage. The rapid increase in fatty acids in the second stage may be closely related to the synergy between the high expression of SAD and low expression of FAD2. In addition, many transcription factors, such as ERF, GRAS, GRF, MADS, MYB and WRKY, may be involved in the fatty acid synthesis. Our data provide a rich resource for further studies on the regulation of fatty acid synthesis in C. chekiangoleosa.