PacBio Single-molecule Real-time Sequencing Research Articles

Integrated PacBio SMRT and Illumina sequencing uncovers transcriptional and physiological responses to drought stress in whole-plant Nitraria tangutorum.

Nitraria tangutorum Bobr., a prominent xerophytic shrub, exhibits remarkable adaptability to harsh environment and plays a significant part in preventing desertification in northwest China owing to its exceptional drought and salinity tolerance. To investigate the drought-resistant mechanism underlying N. tangutorum, we treated 8-week-old seedlings with polyethylene glycol (PEG)-6000 (20%, m/m) to induce drought stress. 27 samples from different tissues (leaves, roots and stems) of N. tangutorum at 0, 6 and 24h after drought stress treatment were sequenced using PacBio single-molecule real-time (SMRT) sequencing and Illumina RNA sequencing to obtain a comprehensive transcriptome. The PacBio SMRT sequencing generated 44,829 non-redundant transcripts and provided valuable reference gene information. In leaves, roots and stems, we identified 1162, 2024 and 232 differentially expressed genes (DEGs), respectively. The Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that plant hormone signaling and mitogen-activated protein kinase (MAPK) cascade played a pivotal role in transmitting stress signals throughout the whole N. tangutorum plant following drought stress. The interconversion of starch and sucrose, as well as the biosynthesis of amino acid and lignin, may represent adaptive strategies employed by N. tangutorum to effectively cope with drought. Transcription factor analysis showed that AP2/ERF-ERF, WRKY, bHLH, NAC and MYB families were mainly involved in the regulation of drought response genes. Furthermore, eight physiological indexes, including content of proline, hydrogen peroxide (H2O2), malondialdehyde (MDA), total amino acid and soluble sugar, and activities of three antioxidant enzymes were all investigate after PEG treatment, elucidating the drought tolerance mechanism from physiological perspective. The weighted gene co-expression network analysis (WGCNA) identified several hub genes serve as key regulator in response to drought through hormone participation, ROS cleavage, glycolysis, TF regulation in N. tangutorum. These findings enlarge genomic resources and facilitate research in the discovery of novel genes research in N. tangutorum, thereby establishing a foundation for investigating the drought resistance mechanism of xerophyte.

Full-Length Transcriptome Construction and Systematic Characterization of Virulence Factor-Associated Isoforms in Vairimorpha (Nosema) Ceranae.

Vairimorpha (Nosema) ceranae is a single-cellular fungus that obligately infects the midgut epithelial cells of adult honeybees, causing bee microsporidiosis and jeopardizing bee health and production. This work aims to construct the full-length transcriptome of V. ceranae and conduct a relevant investigation using PacBio single-molecule real-time (SMRT) sequencing technology. Following PacBio SMRT sequencing, 41,950 circular consensus (CCS) were generated, and 25,068 full-length non-chimeric (FLNC) reads were then detected. After polishing, 4387 high-quality, full-length transcripts were gained. There are 778, 2083, 1202, 1559, 1457, 1232, 1702, and 3896 full-length transcripts that could be annotated to COG, GO, KEGG, KOG, Pfam, Swiss-Prot, eggNOG, and Nr databases, respectively. Additionally, 11 alternative splicing (AS) events occurred in 6 genes were identified, including 1 alternative 5' splice-site and 10 intron retention. The structures of 225 annotated genes in the V. ceranae reference genome were optimized, of which 29 genes were extended at both 5' UTR and 3' UTR, while 90 and 106 genes were, respectively, extended at the 5' UTR as well as 3' UTR. Furthermore, a total of 29 high-confidence lncRNAs were obtained, including 12 sense-lncRNAs, 10 lincRNAs, and 7 antisense-lncRNAs. Taken together, the high-quality, full-length transcriptome of V. ceranae was constructed and annotated, the structures of annotated genes in the V. ceranae reference genome were improved, and abundant new genes, transcripts, and lncRNAs were discovered. Findings from this current work offer a valuable resource and a crucial foundation for molecular and omics research on V. ceranae.

Identification of 43 novel HLA alleles by PacBio single molecule real-time sequencing in haematopoietic cell donors.

A total of 43 novel HLA alleles detected in haematopoietic cell donors using single molecule real-time DNA sequencing.

Effects of sorghum varieties on microbial communities and volatile compounds in the fermentation of light-flavor Baijiu.

Light-flavor Baijiu (LFB) fermentation is a representative spontaneous mixed-culture solid-state fermentation process in which sorghum is used as the raw material. Raw materials and microorganisms are crucial to the flavor formation and quality of LFB. However, the microbial and physicochemical dynamics of different sorghum varieties during LFB fermentation, as well as their impact on flavor compounds are still largely unknown. Herein, PacBio single-molecule real-time (SMRT) sequencing and headspace solid-phase microextraction coupled with gas chromatography-mass spectrometry (HS-SPME-GC-MS) were applied to investigate microbial community succession and volatile flavor formation in glutinous/non-glutinous sorghum-based fermented grains during LFB fermentation. Fermented grains made of glutinous sorghum Liangnuo No. 1 (GLN) had higher bacterial α-diversity and lower fungal α-diversity than those with fermented grains prepared with non-glutinous red sorghum (NRS) (p < 0.05). The dominant microbial species were Saccharomyces cerevisiae, Acetobacter pasteurinus, and Lactobacillus helveticus, the latter two of which were the predominant bacteria observed at the end of fermentation in GLN and NRS, respectively. Moisture content and reducing sugar had a more significant impact on the microorganisms in GLN, while amino acid nitrogen, total free amino acids, and residual starch were the main driving factors driving the microbial community in NRS. The correlation network and discriminant analysis indicated that a relatively high content of 4-vinylguaiacol showed a significant positive association with significant differential microbial species in GLN. These results provided valuable insights for improving the quality of LFB.

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.

A Combinatorial Single-Molecule Real-Time and Illumina Sequencing Analysis of Postembryonic Gene Expression in the Asian Citrus Psyllid Diaphorina citri.

Huanglongbing (HLB) is a systemic plant disease caused by 'Candidatus Liberibacter asiaticus (CLas)' and transmitted by Diaphorina citri. D. citri acquires the CLas bacteria in the nymph stage and transmits it in the adult stage, indicating that molting from the nymph to adult stages is crucial for HLB transmission. However, the available D. citri reference genomes are incomplete, and gene function studies have been limited to date. In the current research, PacBio single-molecule real-time (SMRT) and Illumina sequencing were performed to investigate the transcriptome of D. citri nymphs and adults. In total, 10,641 full-length, non-redundant transcripts (FLNRTs), 594 alternative splicing (AS) events, 4522 simple sequence repeats (SSRs), 1086 long-coding RNAs (lncRNAs), 281 transcription factors (TFs), and 4459 APA sites were identified. Furthermore, 3746 differentially expressed genes (DEGs) between nymphs and adults were identified, among which 30 DEGs involved in the Hippo signaling pathway were found. Reverse transcription-quantitative PCR (RT-qPCR) further validated the expression levels of 12 DEGs and showed a positive correlation with transcriptome data. Finally, the spatiotemporal expression pattern of genes involved in the Hippo signaling pathway exhibited high expression in the D. citri testis, ovary, and egg. Silencing of the D. citri transcriptional co-activator (DcYki) gene significantly increased D. citri mortality and decreased the cumulative molting. Our results provide useful information and a reliable data resource for gene function research of D. citri.

HPV genotyping in clinical samples using long-read single-molecule real-time sequencing.

Human papillomavirus (HPV) genotyping is widely used, particularly in combination with high-risk (HR) HPV tests for cervical cancer screening. We developed a genotyping method using sequences of approximately 800 bp in the E6/E7 region obtained by PacBio single molecule real-time sequencing (SMRT) and evaluated its performance against MY09-11 L1 sequencing and after the APTIMA HPV genotyping assay. The levels of concordance of PacBio E6/E7 SMRT sequencing with MY09-11 L1 sequencing and APTIMA HPV genotyping were 100% and 90.8%, respectively. The sensitivity of PacBio E6/EA7 SMRT was slightly greater than that of L1 sequencing and, as expected, lower than that of HR-HPV tests. In the context of cervical cancer screening, PacBio E6/E7 SMRT is then best used after a positive HPV test. PacBio E6/E7 SMRT genotyping is an attractive alternative for HR and LR-HPV genotyping of clinical samples. PacBio SMRT sequencing provides unbiased genotyping and can detect multiple HPV infections and haplotypes within a genotype.

A Chromosomal-Level Genome of Dermatophagoides farinae, a Common Allergenic Mite Species.

Genome data have been used to find novel allergen from house dust mites. Here, we aim to construct a chromosome-level genome assembly of Dermatophagoides farinae, a common allergenic mite species. We achieved a chromosome-level assembly of D. farinae's genome by integrating PacBio single-molecule real-time sequencing, Illumina paired-end sequencing, and Hi-C technology, followed by annotating allergens and mapping them to specific chromosomes. A 62.43 Mb genome was assembled with a 0.52% heterozygosity rate and a 36.11 Merqury-estimated quality value. The assembled genome represents 92.1% completeness benchmarking universal single-copy orthologs with a scaffold N50 value of 7.11 Mb. Hi-C scaffolding of the genome resulted in construction of 10 pseudochromosomes. The genome comprises 13.01% (7.66 Mb) repetitive sequences and predicts 10,709 protein-coding genes, 96.57% of which are functionally annotated. Moreover, we identified and located 36 allergen groups on specific chromosomes, including allergens Der f 1, Der f 2, Der f 23, Der f 4, Der f 5, Der f 7, and Der f 21 located on chromosomes 2, 1, 7, 3, 4, 6, and 4, respectively. This comprehensive genomic data provides valuable insights into mite biology and evolutionary adaptations, potentially advancing D. farinae allergy research and treatment strategies.

De Novo Transcriptome Analysis by PacBio SMRT-Seq and Illumina RNA-Seq Provides New Insights into Polyphenol Biosynthesis in Chinese Olive Fruit

Polyphenols play a crucial role in fruit flavor. To elucidate the mechanism of fruit polyphenol metabolism, we constructed a transcriptome atlas through PacBio single-molecule real-time (SMRT) sequencing and Illumina next-generation sequencing (NGS) using Canarium album (Lour.) Raeusch., which is a fantastic fruit rich in polyphenolic compounds. In this work, PacBio full-length transcriptome assembly generated 135,439 isoforms with an average length of all isoforms of 2687.94 bp and an N50 length of 3224 bp. To gain deeper insights into the molecular mechanisms of polyphenol biosynthesis in C. album, we constructed twelve RNA-Seq libraries from four developmental stages of the fruits. We identified a total of 28,658 differentially expressed genes (DEGs). We found that many DEGs were involved in metabolic pathways, biosynthesis of secondary metabolites, biosynthesis of antibiotics, starch and sucrose metabolism, and plant hormone signal transduction. Here, we report the expression profiles of 215 DEGs encoding 27 enzymes involved in the polyphenol biosynthesis pathway in C. album. In addition, 285 differentially expressed transcription factors (TFs) continuously down-regulated in four developmental periods of C. album fruit, which may indicate their potential role in the response to polyphenol metabolism and phenylpropanoid biosynthesis pathways. This report will help us understand polyphenol biosynthesis’s functions and metabolic mechanism in C. album. The transcriptome data provide a valuable resource for genetic and genomics research. They will facilitate future work exploiting C. album and other fruits used as medicine and food.

Physiological responses to drought stress of three pine species and comparative transcriptome analysis of Pinus yunnanensis var. pygmaea

Drought stress can significantly affect plant growth, development, and yield. Fewer comparative studies have been conducted between different species of pines, particularly involving Pinus yunnanensis var. pygmaea (P. pygmaea). In this study, the physiological indices, photosynthetic pigment and related antioxidant enzyme changes in needles from P. pygmaea, P. elliottii and P. massoniana under drought at 0, 7, 14, 21, 28 and 35 d, as well as 7 days after rehydration, were measured. The PacBio single-molecule real-time (SMRT) and Illumina RNA sequencing were used to uncover the gene expression differences in P. pygmaea under drought and rehydration conditions. The results showed that the total antioxidant capacity (TAOC) of P. pygmaea was significantly higher than P. massoniana and P. elliottii. TAOC showed a continuous increase trend across all species. Soluble sugar (SS), starch content and non-structural carbohydrate (NSC) of all three pines displayed a "W" pattern, declining initially, increasing, and then decreasing again. P. pygmaea exhibits stronger drought tolerance and greater recovery ability under prolonged drought conditions. Through the PacBio SMRT-seq, a total of 50,979 high-quality transcripts were generated, and 6,521 SSR and 5,561 long non-coding RNAs (LncRNAs) were identified. A total of 2310, 1849, 5271, 5947, 7710, and 6854 differentially expressed genes (DEGs) were identified compared to the control (Pp0D) in six pair-wise comparisons of treatment versus control. bHLH, NAC, ERF, MYB_related, C3H transcription factors (TFs) play an important role in drought tolerance of P. pygmaea. KEGG enrichment analysis and Gene set enrichment analysis (GSEA) analysis showed that P. pygmaea may respond to drought by enhancing metabolic processes such as ABA signaling pathway, alpha-linolenic acid. Weighted gene co-expression network analysis (WGCNA) revealed GST, CAT, LEC14B, SEC23 were associated with antioxidant enzyme activity and TAOC. This study provides a basis for further research on drought tolerance differences among coniferous species.

Bioturbation effect of artificial inoculation on the flavor metabolites and bacterial communities in the Chinese Mao-tofu fermentation

A comparison between artificially inoculated Mao-tofu (CC) and naturally fermented Mao-tofu (MM) indicated that artificially adding Mucor plasmaticus to Mao-tofu dramatically enhanced the essential amino acid (EAA) content, as well as umami and sweet amino acids. Gas chromatography-tandem mass spectrometry (GC–MS/MS) analysis revealed that phenol (3.226 μg/g), 1-octen-3-ol (5.031 μg/g), ethyl heptanoate (1.646 μg/g), and indole (3.422 μg/g) were the key flavor components in Mao-tofu. Unlike MM, CC displayed a substantial increase in esters and a considerable decrease in foul odor substances, including sulfur-containing compounds and indole. Lactococcus raffinolactis, Enterobacter sp. 638, and Streptococcus parauberis KCTC 11537 represented the key bacterial species altering the amino acids and flavor of Mao-tofu according to PacBio single-molecule real-time (SMRT) sequencing and correlation analysis. This study presents the technical feasibility of artificially inoculating Mao-tofu to regulate the core bacterial communities and control the quality of fermented soybean products.

Postmortem Muscle Proteome Characteristics of Silver Carp (Hypophthalmichthys molitrix): Insights from Full-Length Transcriptome and Deep 4D Label-Free Proteomic.

The coverage of the protein database directly determines the results of shotgun proteomics. In this study, PacBio single-molecule real-time sequencing technology was performed on postmortem silver carp muscle transcripts. A total of 42.43 Gb clean data, 35,834 nonredundant transcripts, and 15,413 unigenes were obtained. In total, 99.32% of the unigenes were successfully annotated and assigned specific functions. PacBio long-read isoform sequencing (Iso-Seq) analysis can provide more accurate protein information with a higher proportion of complete coding sequences and longer lengths. Subsequently, 2671 proteins were identified in deep 4D proteomics informed by a full-length transcriptomics technique, which has been shown to improve the identification of low-abundance muscle proteins and potential protein isoforms. The feature of the sarcomeric protein profile and information on more than 30 major proteins in the white dorsal muscle of silver carp were reported here for the first time. Overall, this study provides valuable transcriptome data resources and the comprehensive muscle protein information detected to date for further study into the processing characteristic of early postmortem fish muscle, as well as a spectral library for data-independent acquisition and data processing. This batch of muscle-specific dependent acquisition data is available via PRIDE with identifier PXD043702.

Full-length transcriptome characterization and comparative analysis of Gleditsia sinensis

As an economically important tree, Gleditsia sinensis Lam. is widely planted. A lack of background genetic information on G. sinensis hinders molecular breeding. Based on PacBio single-molecule real-time (SMRT) sequencing and analysis of G. sinensis, a total of 95,183 non-redundant transcript sequences were obtained, of which 93,668 contained complete open reading frames (ORFs), 2,858 were long non-coding RNAs (LncRNAs) and 18,855 alternative splicing (AS) events were identified. Genes orthologous to different Gleditsia species pairs were identified, stress-related genes had been positively selected during the evolution. AGA, AGG, and CCA were identified as the universal optimal codon in the genus of Gleditsia. EIF5A was selected as a suitable fluorescent quantitative reference gene. 315 Cytochrome P450 monooxygenases (CYP450s) and 147 uridine diphosphate (UDP)-glycosyltransferases (UGTs) were recognized through the PacBio SMRT transcriptome. Randomized selection of GsIAA14 for cloning verified the reliability of the PacBio SMRT transcriptome assembly sequence. In conclusion, the research data lay the foundation for further analysis of the evolutionary mechanism and molecular breeding of Gleditsia.

Comparative transcriptome analysis of dioecious floral development in Trachycarpus fortunei using Illumina and PacBio SMRT sequencing

BackgroundTrachycarpus fortunei is a plant with significant economic and ornamental value. Both male and female flowers of T. fortunei originate as bisexual flowers, and selective abortion occurs during floral development. However, the regulatory mechanisms underlying this process remain unclear in T. fortunei. In this study, transcriptome sequencing with Illumina and Pacific BioSciences (PacBio) single-molecule real-time (SMRT) platforms were used to investigate gene expression differences between male and female T. fortunei plants.ResultsA total of 833,137 full-length non-chimeric (FLNC) reads were obtained, and 726,846 high-quality full-length transcripts were identified. A total of 159 genes were differentially expressed between male and female flowers at all development stages. Some of the differentially expressed genes (DEGs) showed male bias, including serine/threonine-protein kinase (STPK), THUMP1 homolog and other genes. Through single-nucleotide polymorphisms(SNPs) identification, 28 genes were considered as potential sex-associated SNPs. Time-Ordered Gene Co-expression Network (TO-GCN) analysis revealed that MADS2 and MADS26 may play important roles in the development of female and male flowers T. fortune plants, respectively.ConclusionsThese findings provide a genetic basis for flower development and differentiation in T. fortunei, and improve our understanding of the mechanisms underlying sexual differentiation in T. fortunei.

Long-Read High-Throughput Sequencing (HTS) Revealed That the Sf-Rhabdovirus X+ Genome Contains a 3.7 kb Internal Duplication.

We previously reported a novel rhabdovirus produced from the Spodoptera frugiperda Sf9 cell line, designated as Sf-rhabdovirus X+ since it contained a unique accessory gene X. The Sf-rhabdovirus X+ genome sequence was generated using Sanger sequencing and short-read high-throughput sequencing (HTS). In this study, we have used long-read HTS technologies, PacBio's single-molecule real-time sequencing and Oxford's Nanopore RNA direct sequencing, to analyze the parent Sf9 cell line transcriptome and the virus RNA produced from an X+ cell clone, respectively. A unique 3.7 kb duplication was identified in the L gene between nucleotide position 8523 and 8524, preceded by a GA dinucleotide insertion. This duplication contained a partial G gene, the complete X gene, and a partial L gene, which extended from nucleotide positions 4767-8523 in the X+ virus. Thus, the X+ genome length is 17,361 nucleotides, and we have re-designated the virus as Sf-rhabdovirus X+3.7. The 3.7 kb duplication was found in all Sf9 cell clones producing the X+ variant virus. Furthermore, the Sf-rhabdovirus X+3.7 genome was stable at passage 30, which was the highest passage tested. These results highlight the importance of combining short-read and long-read technologies for accurately sequencing virus genomes using HTS.

High-quality genome of a pioneer mangrove Laguncularia racemosa explains its advantages for intertidal zone reforestation.

Ecological restoration of mangrove ecosystems that became susceptible to recent habitat perturbations is crucial for tropical coast conservation. The white mangrove Laguncularia racemosa, a pioneer species inhabiting intertidal environments of the Atlantic East Pacific (AEP) region, has been used for reforestation in China for decades. However, the molecular mechanisms underlying its fast growth and high adaptive potential remain unknown. Using PacBio single-molecule real-time sequencing, we completed a high-quality L. racemosa genome assembly covering 1105 Mb with scaffold N50 of 3.46 Mb. Genomic phylogeny shows that L. racemosa invaded intertidal zones during a period of global warming. Multi-level genomic convergence analyses between L. racemosa and three native dominant mangrove clades show that they experienced convergent changes in genes involved in nutrient absorption and high salinity tolerance. This may explain successful L. racemosa adaptation to stressful intertidal environments after introduction. Without recent whole-genome duplications or activated transposable elements, L. racemosa has retained many tandem gene duplications. Some of them are involved in auxin biosynthesis, intense light stress and cold stress response pathways, associated with L. racemosa's ability to grow fast under high light or cold conditions when used for reforestation. In summary, our study identifies shared mechanisms of intertidal environmental adaptation and unique genetic changes underlying fast growth in mangrove-unfavourable conditions and sheds light on the molecular mechanisms of the white mangrove utility in ecological restoration.

Comprehensive transcripts analysis based on single-molecule real-time sequencing and Illumina sequencing provides insights into the mining of Toll-like receptor family in Schizothorax lissolabiatus.

Schizothorax lissolabiatus is an economically important cold-water fish species in southwestern China. Because of water pollution and habitat destruction, the number of wild populations has dramatically decreased. In this study, we used PacBio single-molecule real-time (SMRT) sequencing and Illumina sequencing to generate the first full-length transcriptome and transcriptome, respectively. A total of 19310 polished consensus reads (PC) were obtained, with an average length of 1379 bp and an N50 length of 1485 bp. Meanwhile, 12253 transcripts were successfully annotated as known homologous genes. The pathway annotation indicated that the enrichment and expression of most genes were mainly related to membrane, signal transduction and binding, and immune response. Furthermore, we identified 16 Toll-like receptors (TLRs) by mining the data from the transcripts. Phylogeny analysis showed that S. lissolabiatus TLR genes (slTLRs) supported the classification of TLRs into six families as in other vertebrates. Selection pressure analyses showed that 16 slTLRs revealed purification selection at the overall evolutionary selection. Further, positive selection signals were still detected in eight slTLRs, and most of the positive selection sites were located in the leucine-rich repeat region (LRR domain) associated with the recognition of pathogenic microorganisms, indicating that the function of these slTLR genes may be affected. Tissue specific expression analysis showed all slTLRs are present in kidney, spleen and liver but the relative expression varied among tissues. In conclusion, this study not only provided a valuable resource of transcripts for further research on S. lissolabiatus, but also contributed to improve the current understanding of the evolutionary history of immune-related genes and the TLR gene family in S. lissolabiatus.

Cellulase–lactic acid bacteria synergy action regulates silage fermentation of woody plant

BackgroundFeed shortage is an important factor limiting livestock production in the world. To effectively utilize natural woody plant resources, we used wilting and microbial additives to prepare an anaerobic fermentation feed of mulberry, and used PacBio single-molecule real-time (SMRT) sequencing technology to analyse the “enzyme–bacteria synergy” and fermentation mechanism.ResultsThe fresh branches and leaves of mulberry have high levels of moisture and nutrients, and also contain a diverse range of epiphytic microorganisms. After ensiling, the microbial diversity decreased markedly, and the dominant bacteria rapidly shifted from Gram-negative Proteobacteria to Gram-positive Firmicutes. Lactic acid bacteria (LAB) emerged as the dominant microbial population, resulting in increased in the proportion of the carbohydrate metabolism and decreased in the proportion of the amino acid and “global and overview map” (GOM) metabolism categories. The combination of cellulase and LAB exhibited a synergistic effect, through which cellulases such as glycanase, pectinase, and carboxymethyl cellulase decomposed cellulose and hemicellulose into sugars. LAB converted these sugars into lactic acid through the glycolytic pathway, thereby improving the microbial community structure, metabolism and fermentation quality of mulberry silage. The GOM, carbohydrate metabolism, and amino acid metabolism were the main microbial metabolic categories during ensiling. The presence of LAB had an important effect on the microbial community and metabolic pathways during silage fermentation. A “co-occurrence microbial network” formed with LAB, effectively inhibiting the growth of harmful microorganisms, and dominating the anaerobic fermentation process.ConclusionsIn summary, PacBio SMRT was used to accurately analyse the microbial network information and regulatory mechanism of anaerobic fermentation, which provided a scientific basis for the study of woody silage fermentation theory. This study reveals for the first time the main principle of the enzyme–bacteria synergy in a woody silage fermentation system, which provides technical support for the development and utilization of woody feed resources, and achieves sustainable livestock production.

Microbial communities and their correlation with flavor compound formation during the mechanized production of light-flavor Baijiu

Light-flavor Baijiu fermentation is a typical spontaneous solid-state fermentation process fueled by a variety of microorganisms. Mechanized processes have been increasingly employed in Baijiu production to replace traditional manual operation processes, however, the microbiological and physicochemical dynamics in mechanized processes remain largely unknown. Here, we investigated the microbial community succession and flavor compound formation during a whole mechanized fermentation process of light-flavor Baijiu using the conventional dilution plating method, PacBio single-molecule real-time (SMRT) sequencing and headspace solid-phase microextraction coupled with gas chromatography-mass spectrometry. The results showed that largely different fungal and bacterial communities were involved in the soaking and fermentation processes. A clear succession from Pantoea agglomerans to Bacillus (B.) smithii and B. coagulans in dominant bacterial species and from Cladosporium exasperatum to Saccharomyces cerevisiae and Lichtheimia ramosa in dominant fungal species occurred in the soaking processes. In the fermentation process, the most dominant bacterial species was shifted from Pantoea agglomerans to Lactobacillus (La.) acetotolerans and the most dominant fungal species were shifted from Lichtheimia ramose and Rhizopus arrhizus to Saccharomyces cerevisiae. The bacterial and fungal species positively associated with acidity and the formation of ethanol and different flavor compounds were specified. The microbial species exhibited strong co-occurrence or co-exclusion relationships were also identified. The results are helpful for the improvement of mechanized fermentation process of light-flavor Baijiu production.

SARS-CoV-2 Co-Infections and Recombinations Identified by Long-Read Single-Molecule Real-Time Sequencing.

Co-infection with at least 2 strains of virus is the prerequisite for recombination, one of the means of genetic diversification. Little is known about the prevalence of these events in SARS-CoV-2, partly because it is difficult to detect them. We used long-read PacBio single-molecule real-time (SMRT) sequencing technology to sequence whole genomes and targeted regions for haplotyping. We identified 17 co-infections with SARS-CoV-2 strains belonging to different clades in 6829 samples sequenced between January and October, 2022 (prevalence 0.25%). There were 3 Delta/Omicron co-infections and 14 Omicron/Omicron co-infections (4 cases of 21K/21L, 1 case of 21L/22A, 2 cases of 21L/22B, 4 cases of 22A/22B, 2 cases of 22B/22C and 1 case of 22B/22E). Four of these patients (24%) also harbored recombinant minor haplotypes, including one with a recombinant virus that was selected in the viral quasispecies over the course of his chronic infection. While co-infections remain rare among SARS-CoV-2-infected individuals, long-read SMRT sequencing is a useful tool for detecting them as well as recombinant events, providing the basis for assessing their clinical impact, and a precise indicator of epidemic evolution. IMPORTANCE SARS-CoV-2 variants have been responsible for the successive waves of infection over the 3 years of pandemic. While co-infection followed by recombination is one driver of virus evolution, there have been few reports of co-infections, mainly between Delta and Omicron variants or between the first 2 Omicron variants 21K_BA.1 and 21L_BA.2. The 17 co-infections we detected during 2022 included cases with the recent clades of Omicron 22A, 22B, 22C, and 22E; 24% harbored recombinant variants. This study shows that long-read SMRT sequencing is well suited to SARS-CoV-2 genomic surveillance.