PacBio Single-molecule Long-read Sequencing Research Articles

Comparative Analysis of Alternative Splicing in Moso Bamboo and Its Dwarf Mutant, Phyllostachys edulis ‘Tubaeformis’

Internode length is a crucial phenotypic trait of bamboo, significantly impacting its processing and utilization. Phyllostachys edulis ‘Tubaeformis’ (Shengyin Bamboo), a variety of Moso bamboo, exhibits drastically shortened internodes, making it a valuable ornamental bamboo species. We used PacBio single-molecule long-read sequencing and second-generation sequencing to identify genome-wide alternative splicing (AS) events in Moso bamboo and its dwarf mutant, Shengyin bamboo, and compared the differences between the two. Our sequencing data unveiled 139,539 AS events, with retained introns as the most prevalent events. A large number of genes were differentially alternatively spliced (DAS) between Moso bamboo and Shengyin bamboo, and genes related to RNA splicing were most significantly enriched. The high expression of SR isoforms in the 24th internode of Moso bamboo is likely the main factor leading to its greater number of alternative splicing events. Alternative splicing affects the functional domains of partial GRF, E2F, and NAM isoforms, leading to the loss of domains in some isoforms and enabling some isoforms to acquire new functional domains, and this phenomenon is more common in Shengyin bamboo. AS modifies the functional domains of certain GRF isoforms, frequently resulting in domain losses or endowing isoforms with novel domains, and this phenomenon is more common in Shengyin bamboo. We used PacBio single-molecule long-read sequencing and second-generation sequencing to identify genome-wide alternative splicing (AS) events in Moso bamboo and its dwarf mutant, Shengyin Bamboo and compared the differences between the two.

DNA-m6A calling and integrated long-read epigenetic and genetic analysis with fibertools.

Long-read DNA sequencing has recently emerged as a powerful tool for studying both genetic and epigenetic architectures at single-molecule and single-nucleotide resolution. Long-read epigenetic studies encompass both the direct identification of native cytosine methylation as well as the identification of exogenously placed DNA N6-methyladenine (DNA-m6A). However, detecting DNA-m6A modifications using single-molecule sequencing, as well as coprocessing single-molecule genetic and epigenetic architectures, is limited by computational demands and a lack of supporting tools. Here, we introduce fibertools, a state-of-the-art toolkit that features a semisupervised convolutional neural network for fast and accurate identification of m6A-marked bases using PacBio single-molecule long-read sequencing, as well as the coprocessing of long-read genetic and epigenetic data produced using either PacBio or Oxford Nanopore sequencing platforms. We demonstrate accurate DNA-m6A identification (>90% precision and recall) along >20 kilobase long DNA molecules with a ~1,000-fold improvement in speed. In addition, we demonstrate that fibertools can readily integrate genetic and epigenetic data at single-molecule resolution, including the seamless conversion between molecular and reference coordinate systems, allowing for accurate genetic and epigenetic analyses of long-read data within structurally and somatically variable genomic regions.

Comprehensive Profiling of Alternative Splicing and Alternative Polyadenylation during Fruit Ripening in Watermelon (Citrullus lanatus).

Fruit ripening is a highly complicated process that is accompanied by the formation of fruit quality. In recent years, a series of studies have demonstrated post-transcriptional control play important roles in fruit ripening and fruit quality formation. Till now, the post-transcriptional mechanisms for watermelon fruit ripening have not been comprehensively studied. In this study, we conducted PacBio single-molecule long-read sequencing to identify genome-wide alternative splicing (AS), alternative polyadenylation (APA) and long non-coding RNAs (lncRNAs) in watermelon fruit. In total, 6,921,295 error-corrected and mapped full-length non-chimeric (FLNC) reads were obtained. Notably, more than 42,285 distinct splicing isoforms were derived from 5,891,183 intron-containing full-length FLNC reads, including a large number of AS events associated with fruit ripening. In addition, we characterized 21,506 polyadenylation sites from 11,611 genes, 8703 of which have APA sites. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that fructose and mannose metabolism, starch and sucrose metabolism and carotenoid biosynthesis were both enriched in genes undergoing AS and APA. These results suggest that post-transcriptional regulation might potentially have a key role in regulation of fruit ripening in watermelon. Taken together, our comprehensive PacBio long-read sequencing results offer a valuable resource for watermelon research, and provide new insights into the molecular mechanisms underlying the complex regulatory networks of watermelon fruit ripening.

Establishment of tongue microbiota by 18 months of age and determinants of its microbial profile.

Understanding the development of oral microbiota early in life and the factors that influence it is important for preventing the establishment of dysbiotic oral microbiota later in life. This study demonstrates that the tongue microbiota undergoes early development from 4 to 18 months of age and converges into two types of microbiota showing indications of adult characteristics, with either S. salivarius or Neisseria-dominance. Interestingly, their divergence was strongly determined by their weaning status and the dietary frequencies of sweetened beverages, snacks, and fruits, suggesting that dietary habits during this period might influence the establishment of the oral microbiota. These findings may contribute to the development of novel preventive strategies against oral microbiota-related diseases.

PacBio single-molecule long-read sequencing provides new insights into the complexity of full-length transcripts in oriental river prawn, macrobrachium nipponense

BackgroundOriental river prawn (Macrobrachium nipponense) is one of the most dominant species in shrimp farming in China, which is a rich source of protein and contributes to a significant impact on the quality of human life. Thus, more complete and accurate annotation of gene models are important for the breeding research of oriental river prawn.ResultsA full-length transcriptome of oriental river prawn muscle was obtained using the PacBio Sequel platform. Then, 37.99 Gb of subreads were sequenced, including 584,498 circular consensus sequences, among which 512,216 were full length non-chimeric sequences. After Illumina-based correction of long PacBio reads, 6,599 error-corrected isoforms were identified. Transcriptome structural analysis revealed 2,263 and 2,555 alternative splicing (AS) events and alternative polyadenylation (APA) sites, respectively. In total, 620 novel genes (NGs), 197 putative transcription factors (TFs), and 291 novel long non-coding RNAs (lncRNAs) were identified.ConclusionsIn summary, this study offers novel insights into the transcriptome complexity and diversity of this prawn species, and provides valuable information for understanding the genomic structure and improving the draft genome annotation of oriental river prawn.

High-Resolution Detection of Translocation of Oral Bacteria to the Gut

Ectopic enrichment of oral microbes in the gut is a notable alteration in gut microbial balance. These microbes are likely delivered from the oral cavity with saliva and food; however, evidence of oral–gut microbial transmission is insufficient and needs further investigation. In this observational study, we examined 144 pairs of saliva and stool samples collected from community-dwelling adults to verify the oral–gut microbial link and identify the relevant influencing factors on the increased abundance of oral microbes within the gut. The bacterial composition of each sample was determined using PacBio single-molecule long-read sequencing of the full-length 16S ribosomal RNA gene and amplicon sequence variant (ASV) analysis. Although the bacterial compositions of salivary and gut microbiota were distinctly different, at least 1 ASV was shared between salivary and gut microbiota in 72.9% of subjects. Shared ASVs accounted for 0.0% to 63.1% (median 0.14%) of the gut microbiota in each subject and frequently included abundant Streptococcus salivarius and Streptococcus parasanguinis. Their total relative abundance in the gut was significantly higher in older subjects or those with dental plaque accumulation. The gut microbiota with ≥5% of shared ASVs displayed a higher abundance of Streptococcus, Lactobacillus, and Klebsiella and a lower abundance of Faecalibacterium, Blautia, Megamonas, and Parabacteroides. Our study presents evidence for the translocation of oral bacteria to the gut in community-dwelling adults and suggests that aging and dental plaque accumulation contribute to an increased abundance of oral microbes in the gut, which might be relevant to the compositional shift in the gut commensals.

Full-Length Transcriptome Sequencing and Comparative Transcriptomic Analysis Provide Insights Into the Ovarian Maturation of Exopalaemon carinicauda

The ridgetail white shrimp Exopalaemon carinicauda has the potential to be used as a useful experimental organism in the field of crustacean research due to its transparent body, large egg, and short reproductive cycle. However, little is known about the regulatory mechanisms of ovarian maturation in E. carinicauda. In this study, we applied PacBio single-molecule long-read sequencing technology to unveil the whole transcriptome landscape of Exopalaemon carinicauda to better understand the molecular mechanisms of ovarian development. A total of 49.05 G polymerase read bases were generated, finally, 15, 151 unigenes were successfully identified with a mean length of 4, 257 bp and an N50 of 4, 884 bp from ovary tissues (five different ovarian developmental stages). A total of 14, 007 unigenes were successfully annotated in at least one public database. In addition, 8, 861 coding sequences, 4, 594 lncRNAs and 895 transcription factors were identified. Furthermore, we performed RNA-Seq analysis combined with Iso-Seq results to investigate ovarian maturation regulation mechanism and validated the results by quantitative real-time PCR. We annotated five vitellogenin isoforms in the full-length transcriptome, meanwhile these five genes and vitellogenin receptor suggested an important role during previtellogenesis of E. carinicauda. Genes were involved in PI3K-Akt signaling pathway, retinol metabolism, cell cycle and hedgehog signaling pathway, which showed obvious associations with ovarian maturation of E. carinicauda. We identified 40 DEGs, such as forkhead box protein L2, which may be related to ovarian development in the comparisons. The expression profiles of genes such as ecdysone receptor, ecdysone-induced protein 74EF, ecdysone inducible gene E75 and post-molt protein suggest that they were significantly associated with reproductive molting in the ovarian mature stage. In conclusion, the present study identified important genes and pathways involved in ovarian maturation, which might be useful for studying the reproductive regulation and mechanisms of ovarian maturation in E. carinicauda.

Single-molecule long-read sequencing analysis improves genome annotation and sheds new light on the transcripts and splice isoforms of Zoysia japonica

BackgroundZoysia japonica is an important warm-season turfgrass used worldwide. Although the draft genome sequence and a vast amount of next-generation sequencing data have been published, the current genome annotation and complete mRNA structure remain incomplete. Therefore, to analyze the full-length transcriptome of Z. japonica, we used the PacBio single-molecule long-read sequencing method in this study.ResultsFirst, we generated 37,056 high-confidence non-redundant transcripts from 16,005 gene loci. Next, 32,948 novel transcripts, 913 novel gene loci, 8035 transcription factors, 89 long non-coding RNAs, and 254 fusion transcripts were identified. Furthermore, 15,675 alternative splicing events and 5325 alternative polyadenylation sites were detected. In addition, using bioinformatics analysis, the underlying transcriptional mechanism of senescence was explored based on the revised reference transcriptome.ConclusionThis study provides a full-length reference transcriptome of Z. japonica using PacBio single-molecule long-read sequencing for the first time. These results contribute to our knowledge of the transcriptome and improve the knowledge of the reference genome of Z. japonica. This will also facilitate genetic engineering projects using Z. japonica.

Dynamic Microbial Shifts and Functional Analysis of Saliva Microbial Communities with Caries Children

To observe the dynamic changes in the salivary microbiota of children with dental caries and those who were caries-free and to analyze the functional differences in the oral microecology of the two groups during the course of sugar metabolism and the synthesis and transport of multiple amino acids. Ten children with dental caries and 10 caries-free children were enrolled. We employed Illumina metagenomics technology to analyze the composition and function of salivary microbiome in children with and without caries. Six months later, PacBio single-molecule long-read sequencing technology was used to analyze the changes over time in the oral microbial communities of the two groups. We studied the patterns of change in the oral microbial communities under diseased or healthy conditions and attempted to offer a comprehensive interpretation of children's oral microbiota in terms of its composition and functions. The composition of the oral microbiota of children with or without dental caries changed significantly over time. At the phylum level, changing trends in the salivary microbial communities of children with dental caries were in line with those in caries-free children. In these microbial communities, increased proportions of Firmicutes and decreased proportions of Actinobacteria and Bacteroidetes were found in the two groups. At the genus level, however, the two groups showed significantly different changes of the salivary microbial communities. Upward trends in the abundance of Lactobacillus, Methylobacterium, and Megasphaera were found in the caries group, while the abundance of these genera in the caries-free group showed downward trends. At the species level, L. fermentum, L. gasseri, L. oris, S. downei, and some other species belonging to the genus Lactobacillus showed upward trends in the saliva of children with caries, while they consistently stayed at a relative low level of abundance in caries-free children. The abundance of S. gordonii and S. mutans decreased to a certain extent in children with dental caries, but the abundance of S. gordonii and S. mutans in caries-free children were always at a low level. Species such as S. mutans and C. gracilis were positively correlated to the sum of decayed, missing and filled teeth (dmft), while N. flavescens was negatively correlated to dmft. gltA, icd, and mqo, the key genes related to tricarboxylic acid (TCA) cycle, gudB, a glutamate synthesis-related gene, and argAB/C/J, arginine synthesis-related genes, were significantly increased in caries-free children. In addition, the abundance of the NADH dehydrogenase-related gene nuoB/C/D/E/H/I/J/K/L/M in the electron transport chain increased significantly in caries-free children. Dynamic changes were found in the oral microbiota of children with or without caries. The trends of microbial shifts over time were associated with the oral health status. Oxidative phosphorylation and the synthesis and transport of amino acids such as glutamate and arginine in the oral microecology were more active in caries-free children.

PacBio Single-Molecule Long-Read Sequencing Provides New Light on the Complexity of Full-Length Transcripts in Cattle.

Cattle (Bos taurus) is one of the most widely distributed livestock species in the world, and provides us with high-quality milk and meat which have a huge impact on the quality of human life. Therefore, accurate and complete transcriptome and genome annotation are of great value to the research of cattle breeding. In this study, we used error-corrected PacBio single-molecule real-time (SMRT) data to perform whole-transcriptome profiling in cattle. Then, 22.5 Gb of subreads was generated, including 381,423 circular consensus sequences (CCSs), among which 276,295 full-length non-chimeric (FLNC) sequences were identified. After correction by Illumina short reads, we obtained 22,353 error-corrected isoforms. A total of 305 alternative splicing (AS) events and 3,795 alternative polyadenylation (APA) sites were detected by transcriptome structural analysis. Furthermore, we identified 457 novel genes, 120 putative transcription factors (TFs), and 569 novel long non-coding RNAs (lncRNAs). Taken together, this research improves our understanding and provides new insights into the complexity of full-length transcripts in cattle.

Draft Genome Assembly and Transcriptome Dataset for European Turnip (Brassica rapa L. ssp. rapifera), ECD4 Carrying Clubroot Resistance.

Draft Genome Assembly and Transcriptome Dataset for European Turnip (Brassica rapa L. ssp. rapifera), ECD4 Carrying Clubroot Resistance.

PacBio single molecule long-read sequencing provides insight into the complexity and diversity of the Pinctada fucata martensii transcriptome

BackgroundThe pearl oyster Pinctada fucata martensii is an economically valuable shellfish for seawater pearl production, and production of pearls depends on its growth. To date, the molecular mechanisms of the growth of this species remain poorly understood. The transcriptome sequencing has been considered to understanding of the complexity of mechanisms of the growth of P. f. martensii. The recently released genome sequences of P. f. martensii, as well as emerging Pacific Bioscience (PacBio) single-molecular sequencing technologies, provide an opportunity to thoroughly investigate these molecular mechanisms.ResultsHerein, the full-length transcriptome was analysed by combining PacBio single-molecule long-read sequencing (PacBio sequencing) and Illumina sequencing. A total of 20.65 Gb of clean data were generated, including 574,561 circular consensus reads, among which 443,944 full-length non-chimeric (FLNC) sequences were identified. Through transcript clustering analysis of FLNC reads, 32,755 consensus isoforms were identified, including 32,095 high-quality consensus sequences. After removing redundant reads, 16,388 transcripts were obtained, and 641 fusion transcripts were derived by performing fusion transcript prediction of consensus sequences. Alternative splicing analysis of the 16,388 transcripts was performed after accounting for redundancy, and 9097 gene loci were detected, including 1607 new gene loci and 14,946 newly discovered transcripts. The original boundary of 11,235 genes on the chromosomes was corrected, 12,025 complete open reading frame sequences and 635 long non-coding RNAs (LncRNAs) were predicted, and functional annotation of 13,482 new transcripts was achieved. Two thousand three hundred eighteen alternative splicing events were detected. A total of 228 differentially expressed transcripts (DETs) were identified between the largest (L) and smallest (S) pearl oysters. Compared with the S, the L showed 99 and 129 significantly up-and down-regulated DETs, respectively. Six of these DETs were further confirmed by quantitative real-time RT-PCR (RT-qPCR) in independent experiment.ConclusionsOur results significantly improve existing gene models and genome annotations, optimise the genome structure, and in-depth understanding of the complexity and diversity of the differential growth patterns of P. f. martensii.

A transcriptome atlas of silkworm silk glands revealed by PacBio single-molecule long-read sequencing.

The silk gland of the silkworm Bombyx mori is a specialized organ where silk proteins are efficiently synthesized under precise regulation that largely determines the properties of silk fibers. To understand the genes involved in the regulation of silk protein synthesis, considerable research has focused on the transcripts expressed in silk glands; however, the complete transcriptome profile of this organ has yet to be elucidated. Here, we report a full-length silk gland transcriptome obtained by PacBio single-molecule long-read sequencing technology. In total, 11,697 non-redundant transcripts were identified in mixed samples of silk glands dissected from larvae at five developmental stages. When compared with the published reference, the full-length transcripts optimized the structures of 3002 known genes, and a total of 9061 novel transcripts with an average length of 2171bp were detected. Among these, 1403 (15.5%) novel transcripts were computationally revealed to be lncRNAs, 8135 (89.8%) novel transcripts were annotated to different protein and nucleotide databases, and 5655 (62.4%) novel transcripts were predicted to have complete ORFs. Furthermore, we found 1867 alternative splicing events, 2529 alternative polyadenylation events, 784 fusion events and 6596 SSRs. This study provides a comprehensive set of reference transcripts and greatly revises and expands the available silkworm transcript data. In addition, these data will be very useful for studying the regulatory mechanisms of silk protein synthesis.

PacBio single-molecule long-read sequencing shed new light on the transcripts and splice isoforms of the perennial ryegrass.

Perennial ryegrass (Lolium perenne), one of the most widely used forage and cool-season turfgrass worldwide, has a breeding history of more than 100years. However, the current draft genome annotation and transcriptome characterization are incomplete mainly because of the enormous difficulty in obtaining full-length transcripts. To explore the complete structure of the mRNA and improve the current draft genome, we performed PacBio single-molecule long-read sequencing for full-length transcriptome sequencing in perennial ryegrass. We generated 29,175 high-confidence non-redundant transcripts from 15,893 genetic loci, among which more than 66.88% of transcripts and 24.99% of genetic loci were not previously annotated in the current reference genome. The re-annotated 18,327 transcripts enriched the reference transcriptome. Particularly, 6709 alternative splicing events and 23,789 alternative polyadenylation sites were detected, providing a comprehensive landscape of the post-transcriptional regulation network. Furthermore, we identified 218 long non-coding RNAs and 478 fusion genes. Finally, the transcriptional regulation mechanism of perennial ryegrass in response to drought stress based on the newly updated reference transcriptome sequences was explored, providing new information on the underlying transcriptional regulation network. Taken together, we analyzed the full-length transcriptome of perennial ryegrass by PacBio single-molecule long-read sequencing. These results improve our understanding of the perennial ryegrass transcriptomes and refined the annotation of the reference genome.

Characterization of Full-Length Transcriptome Sequences and Splice Variants of Lateolabrax maculatus by Single-Molecule Long-Read Sequencing and Their Involvement in Salinity Regulation.

Transcriptome complexity plays crucial roles in regulating the biological functions of eukaryotes. Except for functional genes, alternative splicing and fusion transcripts produce a vast expansion of transcriptome diversity. In this study, we applied PacBio single-molecule long-read sequencing technology to unveil the whole transcriptome landscape of Lateolabrax maculatus. We obtained 28,809 high-quality non-redundant transcripts, including 18,280 novel isoforms covering 8,961 annotated gene loci within the current reference genome and 3,172 novel isoforms. A total of 10,249 AS events were detected, and intron retention was the predominant AS event. In addition, 1,359 alternative polyadenylation events, 3,112 lncRNAs, 29,609 SSRs, 365 fusion transcripts, and 1,194 transcription factors were identified in this study. Furthermore, we performed RNA-Seq analysis combined with Iso-Seq results to investigate salinity regulation mechanism at the transcripts level. A total of 518 transcripts were differentially expressed, which were further divided into 8 functional groups. Notably, transcripts from the same genes exhibited similar or opposite expression patterns. Our study provides a comprehensive view of the transcriptome complexity in L. maculatus, which significantly improves current gene models. Moreover, the diversity of the expression patterns of transcripts may enhance the understanding of salinity regulatory mechanism in L. maculatus and other euryhaline teleosts.

PacBio single-molecule long-read sequencing shed new light on the complexity of the Carex breviculmis transcriptome

BackgroundCarex L., a grass genus commonly known as sedges, is distributed worldwide and contributes constructively to turf management, forage production, and ecological conservation. The development of next-generation sequencing (NGS) technologies has considerably improved our understanding of transcriptome complexity of Carex L. and provided a valuable genetic reference. However, the current transcriptome is not satisfactory mainly because of the enormous difficulty in obtaining full-length transcripts.ResultsIn this study, we employed PacBio single-molecule long-read sequencing (SMRT) technology for whole-transcriptome profiling in Carex breviculmis. We generated 60,353 high-confidence non-redundant transcripts with an average length of 2302-bp. A total of 3588 alternative splicing events, and 1273 long non-coding RNAs were identified. Furthermore, 40,347 complete coding sequences were predicted, providing an informative reference transcriptome. In addition, the transcriptional regulation mechanism of C. breviculmis in response to shade stress was further explored by mapping the NGS data to the reference transcriptome constructed by SMRT sequencing.ConclusionsThis study provided a full-length reference transcriptome of C. breviculmis using the SMRT sequencing method for the first time. The transcriptome atlas obtained will not only facilitate future functional genomics studies but also pave the way for further selective and genic engineering breeding projects for C. breviculmis.

Analysis of transcripts and splice isoforms in Medicago sativa L. by single-molecule long-read sequencing.

The full-length transcriptome of alfalfa was analyzed with PacBio single-molecule long-read sequencing technology. The transcriptome data provided full-length sequences and gene isoforms of transcripts in alfalfa, which will improve genome annotation and enhance our understanding of the gene structure of alfalfa. As an important forage, alfalfa (Medicago sativa L.) is world-wide planted. For its complexity of genome and unfinished whole genome sequencing, the sequences and complete structure of mRNA transcripts remain unclear in alfalfa. In this study, single-molecule long-read sequencing was applied to investigate the alfalfa transcriptome using the Pacific Biosciences platform, and a total of 113,321 transcripts were obtained from young, mature and senescent leaves. We identified 72,606 open reading frames including 46,616 full-length ORFs, 1670 transcription factors from 54 TF families and 44,040 simple sequence repeats from 30,797 sequences. A total of 7568 alternative splicing events was identified and the majority of alternative splicing events in alfalfa was intron retention. In addition, we identified 17,740 long non-coding RNAs. Our results show the feasibility of deep sequencing full-length RNA from alfalfa transcriptome on a single-molecule level.

A transcriptome atlas of rabbit revealed by PacBio single-molecule long-read sequencing

It is widely acknowledged that transcriptional diversity largely contributes to biological regulation in eukaryotes. Since the advent of second-generation sequencing technologies, a large number of RNA sequencing studies have considerably improved our understanding of transcriptome complexity. However, it still remains a huge challenge for obtaining full-length transcripts because of difficulties in the short read-based assembly. In the present study we employ PacBio single-molecule long-read sequencing technology for whole-transcriptome profiling in rabbit (Oryctolagus cuniculus). We totally obtain 36,186 high-confidence transcripts from 14,474 genic loci, among which more than 23% of genic loci and 66% of isoforms have not been annotated yet within the current reference genome. Furthermore, about 17% of transcripts are computationally revealed to be non-coding RNAs. Up to 24,797 alternative splicing (AS) and 11,184 alternative polyadenylation (APA) events are detected within this de novo constructed transcriptome, respectively. The results provide a comprehensive set of reference transcripts and hence contribute to the improved annotation of rabbit genome.