PacBio Single-molecule Real-time Sequencing Research Articles

Genome sequencing, assembly, and characterization of Pichia fermentans Z9Y-3 as a non-Saccharomyces yeast with aroma enhancing potential

Pichia fermentans is an important non-Saccharomyces yeast applied in mixed fermentation to enhance fruity and floral traits of wine aroma. In this study, genome sequencing of P. fermentans Z9Y-3 was performed using the PacBio single-molecule real-time sequencing platform, and was further assembled using MECAT. The sequenced genes are annotated by non-redundant protein (Nr), Gene Ontology (GO), Cluster of Orthologous Groups of proteins (COG), Kyoto Encyclopedia of Genes and Genomes (KEGG) and carbohydrate-active enzymes (CAZymes) database. Results indicated that the genome size was about 13.5 Mb with 42.09% GC content and 4730 protein coding genes. The 95.90% genes were annotated to the non-redundant protein (Nr) database, and then annotated to GO database. The 62.64% genes annotated to the COG database revealed 25 functional groups. The 58.22% genes annotated to the KEGG database revealed 5 groups. The 586 genes were annotated to CAZymes, in which two genes bglH and bglC encoding beta-glucosidase were found and further analyzed its motifs. A terpene and a nonribosomal peptides like secondary metabolic gene clusters were predicted using antiSMASH 7.0.0. This study benefits us to understand physiological and biochemical behavior of P. fermentans Z9Y-3 at a molecular level to better use P. fermentans Z9Y-3 in aroma enhancing winemaking.

Full-Length Transcriptome of the Great Himalayan Leaf-Nosed Bats (Hipposideros armiger) Optimized Genome Annotation and Revealed the Expression of Novel Genes.

The Great Himalayan Leaf-nosed bat (Hipposideros armiger) is one of the most representative species of all echolocating bats and is an ideal model for studying the echolocation system of bats. An incomplete reference genome and limited availability of full-length cDNAs have hindered the identification of alternatively spliced transcripts, which slowed down related basic studies on bats' echolocation and evolution. In this study, we analyzed five organs from H. armiger for the first time using PacBio single-molecule real-time sequencing (SMRT). There were 120 GB of subreads generated, including 1,472,058 full-length non-chimeric (FLNC) sequences. A total of 34,611 alternative splicing (AS) events and 66,010 Alternative Polyadenylation (APA) sites were detected by transcriptome structural analysis. Moreover, a total of 110,611 isoforms were identified, consisting of 52% new isoforms of known genes and 5% of novel gene loci, as well as 2112 novel genes that have not been annotated before in the current reference genome of H. armiger. Furthermore, several key novel genes, including Pol, RAS, NFKB1, and CAMK4, were identified as being associated with nervous, signal transduction, and immune system processes, which may be involved in regulating the auditory nervous perception and immune system that helps bats to regulate in echolocation. In conclusion, the full-length transcriptome results optimized and replenished existing H. armiger genome annotation in multiple ways and offer advantages for newly discovered or previously unrecognized protein-coding genes and isoforms, which can be used as a reference resource.

Long-read sequencing improves diagnostic rate in neuromuscular disorders.

Massive parallel sequencing methods, such as exome, genome, and targeted DNA sequencing, have aided molecular diagnosis of genetic diseases in the last 20 years. However, short-read sequencing methods still have several limitations, such inaccurate genome assembly, the inability to detect large structural variants, and variants located in hard-to-sequence regions like highly repetitive areas. The recently emerged PacBio single-molecule real-time (SMRT) and Oxford nanopore technology (ONT) long-read sequencing (LRS) methods have been shown to overcome most of these technical issues, leading to an increase in diagnostic rate. LRS methods are contributing to the detection of repeat expansions in novel disease-causing genes (e.g., ABCD3, NOTCH2NLC and RILPL1 causing an Oculopharyngodistal myopathy or PLIN4 causing a Myopathy with rimmed ubiquitin-positive autophagic vacuolation), of structural variants (e.g., in DMD), and of single nucleotide variants in repetitive regions (TTN and NEB). Moreover, these methods have simplified the characterization of the D4Z4 repeats in DUX4, facilitating the diagnosis of Facioscapulohumeral muscular dystrophy (FSHD). We review recent studies that have used either ONT or PacBio SMRT sequencing methods and discuss different types of variants that have been detected using these approaches in individuals with neuromuscular disorders.

Characterization and Analysis of the Full-Length Transcriptome Provide Insights into Fruit Quality Formation in Kiwifruit Cultivar Actinidia arguta cv. Qinziyu

Kiwifruit is an economically important horticultural crop with extremely high values in nutrition and health care. However, the molecular mechanisms underlying fruit quality formation remain largely limited for most kiwifruit varieties. Recently, a new kiwifruit cultivar with a high level of soluble solids, Actinidia arguta cv. Qinziyu (full-red flesh) was discovered through the introduction and propagation test. To provide new insights into fruit quality formation in a typical kiwifruit cultivar, we integrated full-length transcriptome surveys based on PacBio single-molecule real-time (SMRT) sequencing, key enzyme genes expression involved in carbohydrate and amino acids metabolism pathways, and bHLH gene family analysis to enhance the understanding of soluble sugar, organic acid, and anthocyanin biosynthesis in A. arguta cv. Qinziyu. A total of 175,913 CCSs were generated, of which 124,789 were identified as FLNC transcripts. In total, 45,923 (86.99%) transcripts were successfully annotated, and more than 76.05% of the transcripts were longer than 1 Kb. KEGG pathway analysis showed that 630 candidate genes encoding 55 enzymes were mainly involved in carbohydrate and amino acid biosynthesis pathways. Further analysis verified the expression of 12 key enzyme genes (e.g., pyruvate kinase (PK), enolase (ENO), hexokinase (HK), and phosphoglycerate kinase (PGK)) in flowers using quantitative real-time PCR. Furthermore, we also screened 10 AabHLH proteins’ function in anthocyanin biosynthesis and characterized the AabHLH gene family in A. arguta cv. Qinziyu. Overall, our research data generated by SMRT technology provide the first set of gene isoforms from a full-length transcriptome in A. arguta cv. Qinziyu and more comprehensive insights into the molecular mechanism of fruit quality formation.

Genome Resource for Acinetobacter schindleri H4-3-C1: An Endophyte of Pseudostellaria heterophylla with Degradation Activity to Toxins Produced by Fungal Pathogens

Acinetobacter schindleri is an endophyte of Pseudostellaria heterophylla, a traditional Chinese herbal plant. It has high degradation activity to toxins produced by fungal pathogen Fusarium graminearum. Here, we deployed PacBio single-molecule real-time long-read sequencing technology to generate a complete genome assembly for the Acinetobacter schindleri H4-3-C1 strain and obtained 1.59 Gb of clean reads. These reads were assembled to a single circular DNA chromosome with a length of 3,265,024 bp, and no plasmid was found in the genome. Totals of 3,193 coding sequences, 91 transfer RNA, 21 ribosomal RNA, and 75 small RNAs were identified in the genome. This high-quality genome assembly and gene annotation resource will facilitate the excavation of the zearalenone degradation gene and provide valuable resources for preventing and controlling toxigenic fungal diseases of P. heterophylla. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

A full-length transcriptome and gene expression analysis of three detoxification gene families in a predatory stink bug, Picromerus lewisi.

The predatory stink bug P. Lewisi shows potential for Integrated Pest Management programs for controlling Lepidoptera pest insects in crops and forests. The importance of this insect for biological control has stimulated several studies into its biology and ecology. However, P. lewisi has little genetic information available. In the present study, PacBio single-molecule real-time (SMRT) sequencing and Illumina RNA-seq sequencing technologies were used to reveal the full-length transcriptome profiling and tissue-specific expression patterns of P. lewisi. A total of 12,997 high-quality transcripts with an average length of 2,292 bp were obtained from different stages of P. lewisi using SMRT sequencing. Among these, 12,101 were successfully annotated in seven public databases. A total of 67 genes of cytochrome P450 monooxygenases, 43 carboxylesterase genes, and 18 glutathione S-transferase genes were identified, most of which were obtained with full-length ORFs. Then, tissue-specific expression patterns of 5th instar nymphs were analyzed using Illumina sequencing. Several candidate genes related to detoxification of insecticides and other xenobiotics as well as the degradation of odors, were identified in the guts and antennae of P. lewisi. The current study offered in-depth knowledge to understand the biology and ecology of this beneficial predator and related species.

A high-quality genome assembly of Lactarius hatsudake strain JH5.

Lactarius hatsudake is a species of Lactarius commonly found in pine forests, is edible with a delicious and nutritious fruiting body, and exhibits medicinal properties. It is an ideal natural multifunctional food with bioactive components including fungal polysaccharides, crude fiber, unsaturated fatty acids, nucleic acid derivatives, various amino acids, and vitamins. However, biological and genomic analyses of this mycorrhizal mushroom are sparse, thereby hindering large-scale cultivation. Previously, we isolated and screened L. hatsudake JH5 strains and have applied our garnered knowledge to the large-scale cultivation of mycorrhizal seedlings. In this study, we produced a high-quality genome assembly of L. hatsudake JH5 by combining Illumina paired-end and PacBio single molecule real-time sequencing, resulting in PacBio single molecule real-time reads of 7.67 Gb and Illumina Pair-End reads of 1,560 Mb. Based on the distribution of k-mer frequencies, the genome size of this strain was estimated to be 63.84 Mb (1.14% heterozygosity). Based on de novo genome assembly, the final genome size was determined to be 76.7 Mb, with scaffold N50 of 223.2 kb and N90 of 54.5 kb, and a GC content of 54.38%. BUSCO assessment showed that genome completeness was 89.0%. The N50 length of the JH5 genome was 43.6% longer than that of the previously published L. hatsudake MG20 genome. This high-quality L. hatsudake genome assembly will facilitate research on the functional genome, molecular breeding, yield enhancement, and sustainability of L. hatsudake cultivation.

PacBio single molecule real-time sequencing of a full-length transcriptome of the greenfin horse-faced filefish Thamnaconus modestus

PacBio single molecule real-time sequencing of a full-length transcriptome of the greenfin horse-faced filefish Thamnaconus modestus

Strain-level resolution and pneumococcal carriage dynamics by single-molecule real-time (SMRT) sequencing of the plyNCR marker: a longitudinal study in Swiss infants

BackgroundPneumococcal carriage has often been studied from a serotype perspective; however, little is known about the strain-specific carriage and inter-strain interactions. Here, we examined the strain-level carriage and co-colonization dynamics of Streptococcus pneumoniae in a Swiss birth cohort by PacBio single-molecule real-time (SMRT) sequencing of the plyNCR marker.MethodsA total of 872 nasal swab (NS) samples were included from 47 healthy infants during the first year of life. Pneumococcal carriage was determined based on the quantitative real-time polymerase chain reaction (qPCR) targeting the lytA gene. The plyNCR marker was amplified from 214 samples having lytA-based carriage for pneumococcal strain resolution. Amplicons were sequenced using SMRT technology, and sequences were analyzed with the DADA2 pipeline. In addition, pneumococcal serotypes were determined using conventional, multiplex PCR (cPCR).ResultsPCR-based plyNCR amplification demonstrated a 94.2% sensitivity and 100% specificity for Streptococcus pneumoniae if compared to lytA qPCR. The overall carriage prevalence was 63.8%, and pneumococcal co-colonization (≥ 2 plyNCR amplicon sequence variants (ASVs)) was detected in 38/213 (17.8%) sequenced samples with the relative proportion of the least abundant strain(s) ranging from 1.1 to 48.8% (median, 17.2%; IQR, 5.8–33.4%). The median age to first acquisition was 147 days, and having ≥ 2 siblings increased the risk of acquisition.ConclusionThe plyNCR amplicon sequencing is species-specific and enables pneumococcal strain resolution. We therefore recommend its application for longitudinal strain-level carriage studies of Streptococcus pneumoniae.7kNn5E5Dxp8AY6Abk1AWrtVideo

Whole-genome sequencing of SARS-CoV-2: Comparison of target capture and amplicon single molecule real-time sequencing protocols.

Fast, accurate sequencing methods are needed to identify new variants and genetic mutations of the severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) genome. Single‐molecule real‐time (SMRT) Pacific Biosciences (PacBio) provides long, highly accurate sequences by circular consensus reads. This study compares the performance of a target capture SMRT PacBio protocol for whole‐genome sequencing (WGS) of SARS‐CoV‐2 to that of an amplicon PacBio SMRT sequencing protocol. The median genome coverage was higher (p < 0.05) with the target capture protocol (99.3% [interquartile range, IQR: 96.3–99.5]) than with the amplicon protocol (99.3% [IQR: 69.9–99.3]). The clades of 65 samples determined with both protocols were 100% concordant. After adjusting for C t values, S gene coverage was higher with the target capture protocol than with the amplicon protocol. After stratification on C t values, higher S gene coverage with the target capture protocol was observed only for samples with C t > 17 (p < 0.01). PacBio SMRT sequencing protocols appear to be suitable for WGS, genotyping, and detecting mutations of SARS‐CoV‐2.

Insights into adaptive divergence of Japanese mantis shrimp Oratosquilla oratoria inferred from comparative analysis of full-length transcriptomes

The heterogeneous seascapes in the northwestern Pacific (NWP) can be important selective forces driving adaptive divergence of marine coastal species distributed along the gradients. Here, we tested this hypothesis in Japanese mantis shrimp (Oratosquilla oratoria) with a wide distribution in the NWP and a significant north-south population structure. To this end, the full-length (FL) transcriptomes of northern and southern O. oratoria were firstly sequenced using PacBio single molecule real-time sequencing technology. Based on the FL transcriptome data, we captured large-scale FL transcripts of O. oratoria and predicted the FL transcriptome structure, including coding region, transcription factor and long noncoding RNA. To reveal the divergence between northern and southern O. oratoria, we identified 2,182 pairs of orthologous genes and inferred their sequence divergences. The average differences in coding, 5’ untranslated and 3’ untranslated region were 1.44%, 2.79% and 1.46%, respectively, providing additional support to previous proposition that northern and southern O. oratoria are two species. We provided further evolutionary context to our analysis by identifying positive selected genes (PSGs) between northern and southern O. oratoria. In total, 98 orthologs were found evolving under positive selection and involved several environmentally responsive genes associated with stress response, immunity and cytoskeletal organization, etc. Furthermore, we found PSGs also diverged in gene expression response of northern and southern O. oratoria to heat stress. These findings not only highlight the importance of genetic variation in these genes in adapting to environmental changes in O. oratoria, but also suggest that natural selection may act on the plasticity of gene expression to facilitate O. oratoria adaptation to environmental gradients. Overall, our work contributes to understanding how marine coastal species has evolved to adapt to heterogeneous seascapes in the NWP.

Characterization of the gene expression profile response to drought stress in Haloxylon using PacBio single-molecule real-time and Illumina sequencing.

Haloxylon ammodendron and Haloxylon persicum are important drought-tolerant plants in northwest China. The whole-genome sequencing of H. ammodendron and H. persicum grown in their natural environment is incomplete, and their transcriptional regulatory network in response to drought environment remains unclear. To reveal the transcriptional responses of H. ammodendron and H. persicum to an arid environment, we performed single-molecule real-time (SMRT) and Illumina RNA sequencing. In total, 20,246,576 and 908,053 subreads and 435,938 and 210,334 circular consensus sequencing (CCS) reads were identified by SMRT sequencing of H. ammodendron and H. persicum, and 15,238 and 10,135 unigenes, respectively, were successfully obtained. In addition, 9,794 and 7,330 simple sequence repeats (SSRs) and 838 and 71 long non-coding RNAs were identified. In an arid environment, the growth of H. ammodendron was restricted; plant height decreased significantly; basal and branch diameters became thinner and hydrogen peroxide (H2O2) content and peroxidase (POD) activity were increased. Under dry and wet conditions, 11,803 and 15,217 differentially expressed genes (DEGs) were identified in H. ammodendron and H. persicum, respectively. There were 319 and 415 DEGs in the signal transduction pathways related to drought stress signal perception and transmission, including the Ca2+ signal pathway, the ABA signal pathway, and the MAPK signal cascade. In addition, 217 transcription factors (TFs) and 398 TFs of H. ammodendron and H. persicum were differentially expressed, including FAR1, MYB, and AP2/ERF. Bioinformatic analysis showed that under drought stress, the expression patterns of genes related to active oxygen [reactive oxygen species (ROS)] scavenging, functional proteins, lignin biosynthesis, and glucose metabolism pathways were altered. Thisis the first full-length transcriptome report concerning the responses of H. ammodendron and H. persicum to drought stress. The results provide a foundation for further study of the adaptation to drought stress. The full-length transcriptome can be used in genetic engineering research.

Nasal carriage of CTX-M-55-producing Escherichia coli ST8369 in a healthy cohort in the city of Yangzhou, China.

This study aimed to investigate the prevalence and diversity of extended-spectrum β-lactamases (ESBL)-producing Escherichia coli isolates from healthy individuals in a community and to elucidate their dissemination mechanism. Cefotaxime-resistant E. coli were isolated from 95 samples of healthy persons from one community in Yangzhou, China, and were tested for minimal inhibitory concentrations of 14 antimicrobial agents. The isolates were subjected to whole genome sequencing by Illumina Hiseq or PacBio single-molecule real-time sequencing. A total of 30 cefotaxime-resistant E. coli isolates were obtained, carrying bla CTX-M (n=29) or bla DHA (n=1), of which the bla CTX-M-55 (n=19) was the most predominant genotype. One novel bla CTX-M variant bla CTX-M-252 was identified. Thirteen CTX-M-55-producing E. coli isolates belonged to ST8369 from nasal (n=12) or faecal (n=1) samples shared the identical cgMLST type, resistance profiles, resistance genes, plasmid replicons, and a 5,053-bp bla CTX-M-55 structure ΔIS26-ΔISEcp1-bla CTX-M-55-Δorf477-ΔTn2. The bla CTX-M-55 gene was located on IncHI2/ST3 plasmid in E. coli ST8369. The lengths of bla CTX-M/bla DHA-carrying contigs in the remaining 17 E. coli strains ranged from 1,663 to 382,836 bp, located on chromosome (n=4) or plasmids (n=5); the location of the other eight contigs could not be determined due to incomplete assembly. The bla CTX-M was associated with ISEcp1 as previously reported. Nasal colonization of CTX-M-55-producing ST8369 E. coli strains has occurred among healthy individuals in one community. There is a potential risk of antimicrobial resistance dissemination between humans within one community through close contact or environment via aerosols or dust. Therefore, surveillance of nasal carriage of bla CTX-M in communities is warranted to further monitor the spread of the antimicrobial resistance genes in China.

The genome of the rice planthopper egg parasitoid wasps Anagrus nilaparvatae casts light on the chemo- and mechanosensation in parasitism

BackgroundMymaridae is an ancient insect group and is a basal lineage of the superfamily Chalcidoidea. Species of Mymaridae have great potential for biological control. Anagrus nilaparvatae, a representative species of Mymaridae, is ideal for controlling rice planthopper due to its high rate of parasitism and ability to find hosts efficiently in paddy ridges and fields.ResultsUsing both PacBio single-molecule real-time and Illumina sequencing, we sequenced and assembled the whole genome of A. nilaparvatae, a first for the family Mymaridae. The assembly consists of 394 scaffolds, totaling 488.8 Mb. The assembly is of high continuity and completeness, indicated by the N50 value of 25.4 Mb and 98.2% mapping rate of Benchmarking Universal Single-Copy Orthologs. In total, 16,894 protein-coding genes in the genome were annotated. A phylogenomic tree constructed for A. nilaparvatae and other 12 species of Hymenoptera confirmed that the family Mymaridae is sister to all remaining chalcidoids. The divergence time between A. nilaparvatae and the other seven Chalcidoidea species was dated at ~ 126.9 Mya. Chemoreceptor and mechanoreceptor genes are important in explaining parasitic behavior. We identified 17 odorant binding proteins, 11 chemosensory proteins, four Niemann-Pick type C2 proteins, 88 olfactory receptors, 12 gustatory receptors, 22 ionotropic receptors and 13 sensory neuron membrane proteins in the genome of A. nilaparvatae, which are associated with the chemosensory functions. Strikingly, there is only one pickpocket receptors and nine transient receptor potential genes in the genome that have a mechanosensory function.ConclusionsWe obtained a high-quality genome assembly for A. nilaparvatae using PacBio single-molecule real-time sequencing, which provides phylogenomic insights for its evolutionary history. The small numbers of chemo- and mechanosensory genes in A. nilaparvatae indicate the species-specific host detection and oviposition behavior of A. nilaparvatae might be regulated by relatively simple molecular pathways.

Whole-genome long-read TAPS deciphers DNA methylation patterns at base resolution using PacBio SMRT sequencing technology.

Long-read sequencing provides valuable information on difficult-to-map genomic regions, which can complement short-read sequencing to improve genome assembly, yet limited methods are available to accurately detect DNA methylation over long distances at a whole-genome scale. By combining our recently developed TET-assisted pyridine borane sequencing (TAPS) method, which enables direct detection of 5-methylcytosine and 5-hydroxymethylcytosine, with PacBio single-molecule real-time sequencing, we present here whole-genome long-read TAPS (wglrTAPS). To evaluate the performance of wglrTAPS, we applied it to mouse embryonic stem cells as a proof of concept, and an N50 read length of 3.5 kb is achieved. By sequencing wglrTAPS to 8.2× depth, we discovered a significant proportion of CpG sites that were not covered in previous 27.5× short-read TAPS. Our results demonstrate that wglrTAPS facilitates methylation profiling on problematic genomic regions with repetitive elements or structural variations, and also in an allelic manner, all of which are extremely difficult for short-read sequencing methods to resolve. This method therefore enhances applications of third-generation sequencing technologies for DNA epigenetics.

Integration analysis of PacBio SMRT- and Illumina RNA-seq reveals P450 genes involved in thiamethoxam detoxification in Bradysia odoriphaga

The sciarid fly Bradysia odoriphaga is a serious pest of Chinese chive (Liliaceae). Neonicotinoid insecticides including thiamethoxam have been used for B. odoriphaga control. However, thiamethoxam resistance in B. odoriphaga has developed in recent years. To identify potential genes involved in detoxification metabolism of thiamethoxam in B. odoriphaga, a PacBio single-molecule real-time (SMRT) transcriptome sequencing and Illumina RNA-seq analysis on thiamethoxam treated B. odoriphaga were performed to explore differentially expressed genes in B. odoriphaga. After SMRT sequencing, analysis of Illumina RNA-Seq data showed a total of 172 differentially expressed genes (DEGs) after thiamethoxam treatment, among which eight upregulated DEGs were P450 genes that may be related to thiamethoxam metabolism. The qRT-PCR results of the eight up-regulated P450 unigenes after thiamethoxam treatment were consistent with RNA-Seq data. Furthermore, oral delivery mediated RNA interference of the eight upregulated P450 transcripts followed by insecticide bioassay was conducted, and three P450 unigenes were verified to be related to thiamethoxam detoxification in B. odoriphaga. This study provides new information about the P450 genes involved in thiamethoxam detoxification in B. odoriphaga.

Using PacBio SMRT Sequencing Technology and Metabolomics to Explore the Microbiota-Metabolome Interaction Related to Silage Fermentation of Woody Plant

Silage fermentation is a dynamic process involving the succession of microbial communities and changes in metabolites. Fresh branched and leaves of paper mulberry were used to prepared silage. Crop by-products including corn bran, rice bran, and wheat bran were used as exogenous additives. Pacific Biosciences single-molecule real-time (SMRT) sequencing technology and metabolomics are used to explore the interaction mechanism of microbial structure and metabolites during woody silage fermentation and to verify the principle that exogenous additives can modulate silage fermentation. Under the dual stress of anaerobic and acidic environment of silage fermentation, the microbial community changed from Gram-negative to Gram-positive bacteria, and a large amount of lactic acid and volatile fatty acid were produced, which lowered the pH value and caused the rapid death of aerobic bacteria with thin cell walls. The exogenous additives of corn bran and wheat bran accelerated the dynamic succession of lactic acid bacteria as the dominant microbial community in silage fermentation, increased the metabolic pathways of lactic acid, unsaturated fatty acids, citric acid, L-malic acid and other flavoring agents, and inhibited the growth of Clostridium and Enterobacter, thereby improving the flavor and quality of the silage. However, because rice bran contained butyric acid spore bacteria, it can multiply in an anaerobic environment, led to butyric acid fermentation, and promoted protein degradation and ammonia nitrogen production, thereby reduced the fermentation quality of woody silage. The results showed that during the silage fermentation process, the microbial community and the metabolome can interact, and exogenous additives can affect the fermentation quality of silage. SMRT sequencing technology and untargeted metabolomics revealed the microbiota-metabolome interaction during silage fermentation. Changes in the structure of the microbial community can affect the metabolic pathways, and the final metabolites can inhibit the growth of microorganisms that are not conducive to silage fermentation. Exogenous carbohydrate additives can change the fermentation substrate and affect microbial community structure, thus modulate the silage fermentation.

SMRT Sequencing Reveals Candidate Genes and Pathways With Medicinal Value in Cipangopaludina chinensis.

Cipangopaludina chinensis is an economically important aquatic snail with high medicinal value. However, molecular biology research on C. chinensis is limited by the lack of a reference genome, so the analysis of its transcripts is an important step to study the regulatory genes of various substances in C. chinensis. Herein, we conducted the first full-length transcriptome analysis of C. chinensis using PacBio single-molecule real-time (SMRT) sequencing technology. We identified a total of 26,312 unigenes with an average length of 2,572 bp, of which the largest number of zf-c2h2 transcription factor families (120,18.24%) were found, and also observed that the majority of the 8,058 SSRs contained 4-7 repeat units, which provided data for subsequent work on snail genetics Subsequently, 91.86% (24,169) of the genes were successfully annotated to the four major databases, while the highest homology was observed with Pomacea canaliculata. Functional annotation revealed that the majority of transcripts were enriched in metabolism, signal transduction and Immune-related pathways, and several candidate genes involved in drug metabolism and immune response were identified (e.g., CYP1A1, CYP2J, CYP2U1, GST, ,PIK3, PDE3A, PRKAG). This study lays a foundation for future molecular biology research and provides a reference for studying genes associated with the medicinal value of C. chinensis.

Identification of Abietane-Type Diterpenoids and Phenolic Acids Biosynthesis Genes in Salvia apiana Jepson Through Full-Length Transcriptomic and Metabolomic Profiling.

Salvia apiana (S. apiana) Jepson is a medicinal plant that is frequently used by the Chumash Indians in southern California as a diaphoretic, calmative, diuretic, or antimicrobial agent. Abietane-type diterpenoids (ATDs) and phenolic acids (PAs) are the main bioactive ingredients in S. apiana. However, few studies have looked into the biosynthesis of ATDs and PAs in S. apiana. In this study, using metabolic profiling focused on the ATDs and PAs in the roots and leaves of S. apiana, we found a distinctive metabolic feature with all-around accumulation of ATDs, but absence of salvianolic acid B. To identify the candidate genes involved in these biosynthesis pathways, full-length transcriptome was performed by PacBio single-molecule real-time (SMRT) sequencing. A total of 50 and 40 unigenes were predicted to be involved in ATDs and PAs biosynthesis, respectively. Further transcriptional profile using Illumina HiSeq sequencing showed that the transcriptional variations of these pathways were consistent with the accumulation patterns of corresponding metabolites. A plant kingdom-wide phylogenetic analysis of cytochromes (CYPs) identified two CYP76AK and two CYP76AH subfamily genes that might contribute for the specific ATDs biosynthesis in S. apiana. We also noticed that the clade VII laccase gene family was significantly expanded in Salvia miltiorrhiza compared with that of S. apiana, indicating their involvements in the formation of salvianolic acid B. In conclusion, our results will enable the further understanding of ATDs and PAs biosynthesis in S. apiana and Salvia genus.

Complete genome sequences and genomic characterization of five plasmids harbored by environmentally persistent Cronobacter sakazakii strains ST83 H322 and ST64 GK1025B obtained from powdered infant formula manufacturing facilities

BackgroundCronobacter sakazakii is a foodborne pathogen that causes septicemia, meningitis, and necrotizing enterocolitis in neonates and infants. The current research details the full genome sequences of two extremely persistent C. sakazakii strains (H322 and GK1025B) isolated from powdered infant formula (PIF) manufacturing settings. In addition, the genetic attributes associated with five plasmids, pH322_1, pH322_2, pGK1025B_1, pGK1025B_2, and pGK1025B_3 are described.Materials and MethodsUsing PacBio single-molecule real-time (SMRT®) sequencing technology, whole genome sequence (WGS) assemblies of C. sakazakii H322 [Sequence type (ST)83, clonal complex [CC] 83) and GK1025B (ST64, CC64) were generated. Plasmids, also sequenced, were aligned with phylogenetically related episomes to determine, and identify conserved and missing genomic regions.ResultsA truncated ~ 13 Kbp type 6 secretion system (T6SS) gene cluster harbored on virulence plasmids pH322_2 and pGK1025B_2, and a second large deletion (~ 6 Kbp) on pH322_2, which included genes for a tyrosine-type recombinase/integrase, a hypothetical protein, and a phospholipase D was identified. Within the T6SS of pH322_2 and pGK1025B_2, an arsenic resistance operon was identified which is in common with that of plasmids pSP291_1 and pESA3. In addition, PHASTER analysis identified an intact 96.9 Kbp Salmonella SSU5 prophage gene cluster in pH322_1 and pGK1025B_1 and showed that these two plasmids were phylogenetically related to C. sakazakii plasmids: pCS1, pCsa767a, pCsaC757b, pCsaC105731a. Plasmid pGK1025B_3 was identified as a novel conjugative Cronobacter plasmid. Furthermore, WGS analysis identified a ~ 16.4 Kbp type 4 secretion system gene cluster harbored on pGK1025B_3, which contained a phospholipase D gene, a key virulence factor in several host–pathogen diseases.ConclusionThese data provide high resolution information on C. sakazakii genomes and emphasizes the need for furthering surveillance studies to link genotype to phenotype of strains from previous investigations. These results provide baseline data necessary for future in-depth investigations of C. sakazakii that colonize PIF manufacturing facility settings and genomic analyses of these two C. sakazakii strains and five associated plasmids will contribute to a better understanding of this pathogen's survival and persistence within various “built environments” like PIF manufacturing facilities.