SMRT sequencing of the full-length transcriptome of the Sunda pangolin (Manis javanica)

Exploring the innate immune system of Urechis unicinctus: Insights from full-length transcriptome analysis

The Echiura worm Urechis unicinctus refers to a common benthic invertebrate found in the intertidal zone of Huanghai as well as Bohai Bay. U. unicinctus is known to contain various physiologically active substances, making it highly valuable in terms of its edibility, medicinal properties, and economic potential. Nonetheless, the limited study on the immune system of U. unicinctus poses difficulties for its aquaculture and artificial reproduction. Marine invertebrates, including shellfish and U. unicinctus, are thought to primarily depend on their innate immune system for disease protection, owing to the severalinnate immune molecules they possess. Herein, we employed PacBio single-molecule real-time (SMRT) sequencing technology to perform the full-length transcriptome analysis of U. unicinctus individuals under five different conditions (room temperature (RT), low temperature (LT), high temperature (HT), without water (DRY), ultraviolet irradiation (UV)). Concequently, we identified 59,371 unigenes that had a 2,779 bp average length, 2,613 long non-coding RNAs (lncRNAs), 59,190 coding sequences (CDSs), 35,166 simple sequence repeats (SSRs), and 1,733 transcription factors (TFs), successfully annotating 90.58 % (53,778) of the unigenes. Subsequently, key factors associated with immune-related processes, such as non-self-recognition, cellular immune defenses, and humoral immune defenses, were searched. Our study also identified pattern recognition receptors (PRRs) that included 17 peptidoglycan recognition proteins (PGRPs), 13 Gram-negative binding proteins (GNBPs), 18 scavenger receptors (SRs), 74 toll-like receptors (TLRs), and 89 C-type lectins (CLTs). Altogether, the high-quality transcriptome obtained data will offer valuable insights for further investigations into U. unicinctus innate immune response, laying the foundation for subsequent molecular biology studies and aquaculture.

Evolution and conservation genetics of pangolins

Pangolins (Pholidota, Manidae) are classified as an evolutionarily distinct and globally endangered mammal due to their unique morphology (nail-like scales and a myrmecophagous diet) and being the victim of heavy poaching and worldwide trafficking. As such, pangolins serve as a textbook example for studying the special phenotypic evolutionary adaptations and conservation genetics of an endangered species. Recent years have demonstrated significant advancements in the fields of molecular genetics and genomics, which have translated to a series of important research achievements and breakthroughs concerning the evolution and conservation genetics of pangolins. This review comprehensively presents the hitherto advances in phylogeny, adaptive evolution, conservation genetics, and conservation genomics that are related to pangolins, which will provide an ample understanding of their diversity, molecular adaptation mechanisms, and evolutionary potentials. In addition, we highlight the priority of investigating species/population diversity among pangolins and suggest several avenues of research that are highly relevant for future pangolin conservation.

Metagenomics uncovers dietary adaptations for chitin digestion in the gut microbiota of convergent myrmecophagous mammals

IMPORTANCE

Myrmecophagous mammals are specialized in the consumption of ants and/or termites. They do not share a direct common ancestor and evolved convergently in five distinct placental orders raising questions about the underlying adaptive mechanisms involved and the relative contribution of natural selection and phylogenetic constraints. Understanding how these species digest their prey can help answer these questions. More specifically, the role of their gut microbial symbionts in the digestion of the insect chitinous exoskeleton has not been investigated in all myrmecophagous orders. We generated 29 new gut metagenomes from nine myrmecophagous species to reconstruct more than 300 bacterial genomes in which we identified chitin-degrading enzymes. Studying the distribution of these chitinolytic bacteria among hosts revealed both shared and specific bacteria between ant-eating species. Overall, our results highlight the potential role of gut symbionts in the convergent dietary adaptation of myrmecophagous mammals and the evolutionary mechanisms shaping their gut microbiota.

SMRT Sequencing Technology Was Used to Construct the Batocera horsfieldi (Hope) Transcriptome and Reveal Its Features

Batocera horsfieldi (Hope) (Coleoptera: Cerambycidae) is an important forest pest in China that mainly infests timber and economic forests. This pest primarily causes plant tissue to necrotize, rot, and eventually die by feeding on the woody parts of tree trunks. To gain a deeper understanding of the genetic mechanism of B. horsfieldi, this study employed single-molecule real-time sequencing (SMRT) and Illumina RNA-seq technologies to conduct full-length transcriptome sequencing of the insect. Total RNA extracted from male and female adults was mixed and subjected to SMRT sequencing, generating a complete transcriptome. Transcriptome analysis, prediction of long non-coding RNA (lncRNA), coding sequences (CDs), analysis of simple sequence repeats (SSR), prediction of transcription factors, and functional annotation of transcripts were performed in this study. The collective 20,356,793 subreads (38.26 G, clean reads) were generated, including 432,091 circular consensus sequences and 395,851 full-length non-chimera reads. The full-length non-chimera reads (FLNC) were clustered and redundancies were removed, resulting in 39,912 consensus reads. SSR and ANGEL software v3.0 were used for predicting SSR and CDs. In addition, four tools were used for annotating 6058 lncRNAs, identifying 636 transcription factors. Furthermore, a total of 84,650 transcripts were functionally annotated in seven different databases. This is the first time that the full-length transcriptome of B. horsfieldi has been obtained using SMRT sequencing. This provides an important foundation for investigating the gene regulation underlying the interaction between B. horsfieldi and its host plants through gene editing in the future and provides a scientific basis for the prevention and control of B. horsfieldi.

Comprehensive transcriptome characterization of Grus japonensis using PacBio SMRT and Illumina sequencing

The red-crowned crane (Grus japonensis) is an endangered species distributed across southeast Russia, northeast China, Korea, and Japan. Here, we sequenced for the first time the full-length unreferenced transcriptome of red-crowned crane mixed samples using a PacBio Sequel platform. A total of 359,136 circular consensus sequences (CCS) were obtained via clustering to remove redundancy. A total of 303,544 full-length non-chimeric sequences were identified by judging whether CCS contained 5' and 3' adapters, and the poly(A) tail. Eight samples were sequenced using Illumina, and PacBio sequencing data were corrected according to the collected Illumina data to obtain more accurate full-length transcripts. A total of 4,100 long non-coding RNAs, 13,115 simple sequences repeat loci and 29 transcription factor families were identified. The expression of lncRNAs and TFs in pancreas was lowest comparing with other tissues. Many enriched immune-related transmission pathways (MHC and IL receptors) were identified in the spleen. This study will contribute to a better understanding of the gene structure and post-transcriptional regulatory network, and provide references for future studies on red-crowned cranes.

SMRT Sequencing of the Full-Length Transcriptome of the Coelomactra antiquata

Coelomactra antiquata is an important aquatic economic shellfish with high medicinal value. However, because C. antiquata has no reference genome, a lot of molecular biology research cannot be carried out, so the analysis of its transcripts is an important step to study the regulatory genes of various substances in C. antiquata. In the present study, we conducted the first full-length transcriptome analysis of C. antiquata by using PacBio single-molecule real-time (SMRT) sequencing technology. The results identified a total of 39,209 unigenes with an average length of 2,732 bp, 23,338 CDSs, 251 AS events, 9,881 lncRNAs, 20,106 SSRs, and 2,316 TFs. Subsequently, 59.22% (23,220) of the unigenes were successfully annotated, of which 23,164, 18,711, 15,840, 13,534, and 13,474 unigenes could be annotated using NR, Swiss-prot, KOG, GO, and KEGG databases, respectively. This study lays the foundation for the follow-up research of molecular biology and provides a reference for studying the more medicinal value of C. antiquata.

Full-length transcriptome analysis of pecan (Carya illinoinensis) kernels

Pecan is rich in bioactive components such as fatty acids (FAs) and flavonoids and is an important nut type worldwide. Therefore, the molecular mechanisms of phytochemical biosynthesis in pecan are a focus of research. Recently, a draft genome and several transcriptomes have been published. However, the full-length mRNA transcripts remain unclear, and the regulatory mechanisms behind the quality components biosynthesis and accumulation have not been fully investigated. In this study, single-molecule long-read sequencing technology was used to obtain full-length transcripts of pecan kernels. In total, 37,504 isoforms of 16,702 genes were mapped to the reference genome. The numbers of known isoforms, new isoforms, and novel isoforms were 9013 (24.03%), 26,080 (69.54%), and 2411 (6.51%), respectively. Over 80% of the transcripts (30,751, 81.99%) had functional annotations. A total of 15,465 alternative splicing (AS) events and 65,761 alternative polyadenylation events were detected; wherein, the retained intron was the predominant type (5652, 36.55%) of AS. Furthermore, 1894 long noncoding RNAs and 1643 transcription factors were predicted using bioinformatics methods. Finally, the structural genes associated with FA and flavonoid biosynthesis were characterized. A high frequency of AS accuracy (70.31%) was observed in FA synthesis-associated genes. This study provides a full-length transcriptome data set of pecan kernels, which will significantly enhance the understanding of the regulatory basis of phytochemical biosynthesis during pecan kernel maturation.

Comprehensive analysis of full-length transcriptomes of Schizothorax prenanti by single-molecule long-read sequencing

Schizothorax prenanti (S. prenanti) is one of the most important aquaculture species in the southwest of China. However, information of the full-length transcripts in S. prenanti remains unknown. In this study, single-molecule real-time (SMRT) sequencing was performed to generate full-length transcriptomes of S.prenanti. In total, 23.26 Gb of clean reads were generated. A total of 312,587 circular consensus sequences (CCS) were obtained with average lengths of 2634 bp and 84.16% (270,662) of CCS were full-length non-chimeric reads. After being corrected with Illumina library sequencing, 18,005 contigs were obtained, with 17,797 (98.81%) successfully annotated in eight public databases, including 15,839 complete open reading frames (ORFs) with an average length of 1330 bp. Furthermore, a total of 4152 alternative splicing (AS) events and 250 long non-coding RNA (lncRNA) transcripts were detected. Additionally, a total of 1129 putative transcription factors (TFs) members from 56 TF families and 11,660 simple sequence repeats (SSRs) were identified. This study provided a valuable resource of full-length transcripts for further research on S. prenanti.

BSR-Seq analysis provides insights into the cold stress response of Actinidia arguta F1 populations

BACKGROUND

Freezing injury, which is an important abiotic stress in horticultural crops, influences the growth and development and the production area of kiwifruit (Actinidia Lind1). Among Actinidia species, Actinidia arguta has excellent cold resistance, but knowledge relevant to molecular mechanisms is still limited. Understanding the mechanism underlying cold resistance in kiwifruit is important for breeding cold resistance.

RESULTS

In our study, a population resulting from the cross of A. arguta 'Ruby-3' × 'Kuilv' male was generated for kiwifruit hardiness study, and 20 cold-tolerant and 20 cold-sensitive populations were selected from 492 populations according to their LT50. Then, we performed bulked segregant RNA-seq combined with single-molecule real-time sequencing to identify differentially expressed genes that provide cold hardiness. We found that the content of soluble sucrose and the activity of β-amylase were higher in the cold-tolerant population than in the cold-sensitive population. Upon - 30 °C low-temperature treatment, 126 differentially expressed genes were identify; the expression of 59 genes was up-regulated and that of 67 genes was down-regulated between the tolerant and sensitive pools, respectively. KEGG pathway analysis showed that the DEGs were primarily related to starch and sucrose metabolism, amino sugar and nucleotide sugar metabolism. Ten major key enzyme-encoding genes and two regulatory genes were up-regulated in the tolerant pool, and regulatory genes of the CBF pathway were found to be differentially expressed. In particular, a 14-3-3 gene was down-regulated and an EBF gene was up-regulated. To validate the BSR-Seq results, 24 DEGs were assessed via qRT-PCR, and the results were consistent with those obtained by BSR-Seq.

CONCLUSION

Our research provides valuable insights into the mechanism related to cold resistance in Actinidia and identified potential genes that are important for cold resistance in kiwifruit.

Atlas of ACE2 gene expression reveals novel insights into transmission of SARS-CoV-2

The recent pandemic, COVID-19, is caused by a novel coronavirus, SARS-CoV-2, with elusive origin. SARS-CoV-2 infects mammalian cells via ACE2, a transmembrane protein. Therefore, the conservation and expression patterns of ACE2 may provide valuable insights into tracing the carriers of SARS-CoV-2. In this work, we analyzed the conservation of ACE2 and its expression pattern among various mammalian species that are close to human beings. We show that mammalian ACE2 gene is deeply conserved at both DNA and peptide levels, suggesting that a broad range of mammals can potentially host SARS-CoV-2. We further report that ACE2 expression in certain human tissues are consistent with clinical symptoms of COVID-19 patients. Furthermore, we have built the first atlas of ACE2 expression in various common mammals, which shows that ACE2 expresses in mammalian tissues in a species-specific manner. Most notably, we observe exceptionally high expression of ACE2 in external body parts of cats and dogs, suggesting that these household pet animals could be vulnerable to viral infections and/or may serve as intermediate hosts, thus yielding novel insights into the transmission of SARS-CoV-2.

Full-length transcriptome analysis provides new insights into the innate immune system of Marsupenaeus japonicus

As invertebrates, shrimp are generally thought to solely rely on their innate immune system to combat invading pathogens. Recently, an increasing number of studies have revealed that the innate immune response of invertebrates exhibits diversity and specificity based on their diverse immune molecules. Herein, a full-length transcriptome analysis of several immune-related tissues (hepatopancreas, gill, hemocytes, stomach and intestine) in the kuruma shrimp (Marsupenaeus japonicus) was conducted to identify immune-related molecules with a focus on transcript variations. In total, 11,222 nonredundant full-length transcripts with an N50 length of 5174 were obtained, and most of these transcripts (94.84%) were successfully annotated. In addition, a total of 147 long noncoding RNAs (lncRNAs) were also predicted. Importantly, transcript variants of several vital immune-related genes were observed, including twenty-five alpha-2-macroglobulins (α2-Ms), ten Toll-like receptors (TLRs), six C-type lectins (CTLs), five M-type lectins (MTLs) and three Down syndrome cell adhesion molecules (Dscams). Furthermore, 509 nonredundant full-length transcripts were predicted to be generated from alternative splicing (AS) events, which contribute to the diversity of immune molecules. Overall, our study provides valuable data on the full-length transcripts of M. japonicus, which will facilitate the exploration of immune molecules in this species. Moreover, numerous transcript variants of immune molecules detected in this study provide clues for further investigating the diversity and specificity of the innate immune response in shrimp.

Comparative Transcriptome Analysis Combining SMRT- and Illumina-Based RNA-Seq Identifies Potential Candidate Genes Involved in Betalain Biosynthesis in Pitaya Fruit

To gain more valuable genomic information about betalain biosynthesis, the full-length transcriptome of pitaya pulp from ‘Zihonglong’ (red pulp) and ‘Jinghonglong’ (white pulp) in four fruit developmental stages was analyzed using Single-Molecule Real-Time (SMRT) sequencing corrected by Illumina RNA-sequence (Illumina RNA-Seq). A total of 65,317 and 91,638 genes were identified in ‘Zihonglong’ and ‘Jinghonglong’, respectively. A total of 11,377 and 15,551 genes with more than two isoforms were investigated from ‘Zihonglong’ and ‘Jinghonglong’, respectively. In total, 156,955 genes were acquired after elimination of redundancy, of which, 120,604 genes (79.63%) were annotated, and 30,875 (20.37%) sequences without hits to reference database were probably novel genes in pitaya. A total of 31,169 and 53,024 simple sequence repeats (SSRs) were uncovered from the genes of ‘Zihonglong’ and ‘Jinghonglong’, and 11,650 long non-coding RNAs (lncRNAs) in ‘Zihonglong’ and 11,113 lncRNAs in ‘Jinghonglong’ were obtained herein. qRT-PCR was conducted on ten candidate genes, the expression level of six novel genes were consistent with the Fragments Per Kilobase of transcript per Million mapped reads (FPKM) values. In conclusion, we firstly undertook SMRT sequencing of the full-length transcriptome of pitaya, and the valuable resource that was acquired through this sequencing facilitated the identification of additional betalain-related genes. Notably, a list of novel putative genes related to the synthesis of betalain in pitaya fruits was assembled. This may provide new insights into betalain synthesis in pitaya.

Full-Length Transcriptome Analysis of the ABCB, PIN/PIN-LIKES, and AUX/LAX Families Involved in Somatic Embryogenesis of Lilium pumilum DC. Fisch

Plant cell totipotency is one of the 25 major topics in current scientific research, and somatic embryos are good experimental material for studying cell totipotency. Polar auxin transport plays an important regulatory role in somatic embryogenesis (SE). However, little is known about the auxin transport genes and their regulatory mechanisms in Lilium SE. In this study, we applied single-molecule real-time (SMRT) sequencing to Lilium pumilum DC. Fisch. for the first time and obtained a total of 119,649 transcripts, of which 14 encoded auxin transport genes. Correlation analyses between somatic embryo induction and gene expression under different treatments revealed that auxin transport genes, especially ATP-binding cassette (ABC) transporter B family member 21 (ABCB21) and PIN-FORMED (PIN) LIKES 7 (PILS7), may be key players in SE, and the necessary duration of picloram (PIC) treatment to induce SE is as short as 3 days. Our research provides valuable genetic information on Lilium pumilum, elucidating the candidate auxin transport genes involved in SE and their influencing factors. This study lays a foundation for elucidating the regulatory mechanism of auxin transport in SE.