The transcriptome response of Enterobacter sp. S-33 is modulated by low pH-stress.

BACKGROUND: Acidic environments naturally occur worldwide and uncontrolled use of agricultural practices may also cause acidification of soils. The development of acidic conditions disturbs the establishment of efficient microbial populations in their natural niches. The survival of Enterobacter species under acidic stress remains poorly understood. OBJECTIVE: This study aimed to investigate the survival of an environmental isolate Enterobacter sp. S-33 under acidic stress and to identify the various genes involved in stress protection at the global gene transcription level. The obtained results provide new targets that will allow understanding the in-depth mechanisms involved in the adaptation of bacteria to environmental pH changes. METHODS: We used the next-generation sequencing (NGS) method to analyze the expression (up-regulation & down-regulation) of genes under varying pH conditions. RESULTS: A total of 4214 genes were differentially expressed under acidic conditions (pH 5.0), with 294 up-regulated and 167 down-regulated. At pH 6.0, 50 genes were significantly expressed, of which 34 and 16 were identified as up-regulated and down-regulated, respectively. Many of the up-regulated genes were involved in carbohydrate metabolism, amino acid transport & metabolism, and the most down-regulated genes were related to post-translational modification, lipid transport & metabolism, etc. The observed transcriptomic regulation of genes and pathways identified that Enterobacter reduced its post-translational modification, lipid transport & metabolism, and increased carbohydrate metabolism, amino acid metabolism & transport, energy production & conversion to adapt and grow in acidic stress. CONCLUSIONS: The present work provides in-depth information on the characterization of genes associated with tolerance or adaptation to acidic stress of Enterobacter bacterium.

Protein profiling and immunoinformatic analysis of the secretome of a metal-resistant environmental isolate Pseudomonas aeruginosa S-8.

The bacterial secretome represents a comprehensive catalog of proteins released extracellularly that have multiple important roles in virulence and intercellular communication. This study aimed to characterize the secretome of an environmental isolate Pseudomonas aeruginosa S-8 by analyzing trypsin-digested culture supernatant proteins using nano-LC-MS/MS tool. Using a combined approach of bioinformatics and mass spectrometry, 1088 proteins in the secretome were analyzed by PREDLIPO, SecretomeP 2.0, SignalP 4.1, and PSORTb tool for their subcellular localization and further categorization of secretome proteins according to signal peptides. Using the gene ontology tool, secretome proteins were categorized into different functional categories. KEGG pathway analysis identified the secreted proteins into different metabolic functional pathways. Moreover, our LC-MS/MS data revealed the secretion of various CAZymes into the extracellular milieu, which suggests its strong biotechnological applications to breakdown complex carbohydrate polymers. The identified immunodominant epitopes from the secretome of P. aeruginosa showed the characteristic of being non-allergenic, highly antigenic, nontoxic, and having a low risk of triggering autoimmune responses, which highlights their potential as successful vaccine targets. Overall, the identification of secreted proteins of P. aeruginosa could be important for both diagnostic purposes and the development of an effective candidate vaccine.

Unravelling the transcriptome response of Enterobacter sp. S-33 under varying temperature.

Enterobacter genus includes the bacteria occupying every aspect of environment, however, their adaptability at varying temperature is not clear. In the present study, we analyzed the transcriptome response of Enterobacter sp. S-33 and their functional genes under various temperatures (30-45 ℃) that were expressed and accumulated in cells under temperature-stress. During a temperature shift from 37 to 45 ℃, 165 genes showed differential expression including 112 up-regulated and 53 down-regulated. In particular, heat-shock genes such as CspA, 16 kDa heat shock protein A/B and transcriptional regulators such as LysR, TetR, and LuxR were differentially expressed, indicating the role of complex molecular mechanism of Enterobacter adapting to temperature stress. Similarly, genes associated to signal transduction, ABC transporters, iron homeostasis, and quorum sensing were also induced. The Gene ontology enrichment analysis of differentially expressed genes (DEGs) were categorized into "transmembrane transport", "tRNA binding", "hydrogen sulfide biosynthetic process" and "sulfate assimilation" terms. In addition, Kyoto Encyclopedia of Genes and Genomes pathways showed that ABC transporter as well as quorum sensing pathways were significantly enriched. Overall, current study has contributed to explore the adaptive molecular mechanisms of Enterobacter spp. upon temperature change, which further opens new avenues for future in-depth functional studies.

An online database for einkorn wheat to aid in gene discovery and functional genomics studies.

Diploid A-genome wheat (einkorn wheat) presents a nutrition-rich option as an ancient grain crop and a resource for the improvement of bread wheat against abiotic and biotic stresses. Realizing the importance of this wheat species, reference-level assemblies of two einkorn wheat accessions were generated (wild and domesticated). This work reports an einkorn genome database that provides an interface to the cereals research community to perform comparative genomics, applied genetics and breeding research. It features queries for annotated genes, the use of a recent genome browser release, and the ability to search for sequence alignments using a modern BLAST interface. Other features include a comparison of reference einkorn assemblies with other wheat cultivars through genomic synteny visualization and an alignment visualization tool for BLAST results. Altogether, this resource will help wheat research and breeding. Database URL  https://wheat.pw.usda.gov/GG3/pangenome.

In-depth genome and pan-genome analysis of a metal-resistant bacterium Pseudomonas parafulva OS-1.

A metal-resistant bacterium Pseudomonas parafulva OS-1 was isolated from waste-contaminated soil in Ranchi City, India. The isolated strain OS-1 showed its growth at 25-45°C, pH 5.0-9.0, and in the presence of ZnSO4 (upto 5 mM). Phylogenetic analysis based on 16S rRNA gene sequences revealed that strain OS-1 belonged to the genus Pseudomonas and was most closely related to parafulva species. To unravel the genomic features, we sequenced the complete genome of P. parafulva OS-1 using Illumina HiSeq 4,000 sequencing platform. The results of average nucleotide identity (ANI) analysis indicated the closest similarity of OS-1 to P. parafulva PRS09-11288 and P. parafulva DTSP2. The metabolic potential of P. parafulva OS-1 based on Clusters of Othologous Genes (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) indicated a high number of genes related to stress protection, metal resistance, and multiple drug-efflux, etc., which is relatively rare in P. parafulva strains. Compared with other parafulva strains, P. parafulva OS-1 was found to have the unique ß-lactam resistance and type VI secretion system (T6SS) gene. Additionally, its genomes encode various CAZymes such as glycoside hydrolases and other genes associated with lignocellulose breakdown, suggesting that strain OS-1 have strong biomass degradation potential. The presence of genomic complexity in the OS-1 genome indicates that horizontal gene transfer (HGT) might happen during evolution. Therefore, genomic and comparative genome analysis of parafulva strains is valuable for further understanding the mechanism of resistance to metal stress and opens a perspective to exploit a newly isolated bacterium for biotechnological applications.

Molecular characterization and in-depth genome analysis of Enterobacter sp. S-16.

Enterobacter species are considered to be an opportunistic human pathogen owing to the existence of antibiotic-resistant strains and drug resides; however, the detailed analysis of the antibiotic resistance and virulence features in environmental isolates is poorly characterized. Here, in the study, we characterized the biochemical characteristics, and genome, pan-genome, and comparative genome analyses of an environmental isolate Enterobacter sp. S-16. The strain was identified as Enterobacter spp. by using 16S rRNA gene sequencing. To unravel genomic features, whole genome of Enterobacter sp. S-16 was sequenced using a hybrid assembly approach and genome assembly was performed using the Unicycler tool. The assembled genome contained the single conting size 5.3 Mbp, GC content 55.43%, and 4500 protein-coding genes. The genome analysis revealed the various gene clusters associated with virulence, antibiotic resistance, type VI secretion system (T6SS), and many stress tolerant genes, which may provide important insight for adapting to changing environment conditions. Moreover, different metabolic pathways were identified that potentially contribute to environmental survival. Various hydrolytic enzymes and motility functions equipped the strain S-16 as an active colonizer. The genome analysis confirms the presence of carbohydrate-active enzymes (CAZymes), and non-enzymatic carbohydrate-binding modules (CBMs) involved in the hydrolysis of complex carbohydrate polymers. Moreover, the pan-genome analysis provides detailed information about the core genes and shared genes with the closest related Enterobacter species. The present study is the first report showing the presence of YdhE/NorM in Enterobacter spp. Thus, the elucidation of genome sequencing will increase our understanding of the pathogenic nature of environmental isolate, supporting the One Health Concept.

Characterization and in-depth genome analysis of a halotolerant probiotic bacterium Paenibacillus sp. S-12, a multifarious bacterium isolated from Rauvolfia serpentina.

BACKGROUND: Members of Paenibacillus genus from diverse habitats have attracted great attention due to their multifarious properties. Considering that members of this genus are mostly free-living in soil, we characterized the genome of a halotolerant environmental isolate belonging to the genus Paenibacillus. The genome mining unravelled the presence of CAZymes, probiotic, and stress-protected genes that suggested strain S-12 for industrial and agricultural purposes. RESULTS: Molecular identification by 16 S rRNA gene sequencing showed its closest match to other Paenibacillus species. The complete genome size of S-12 was 5.69 Mb, with a GC-content 46.5%. The genome analysis of S-12 unravelled the presence of an open reading frame (ORF) encoding the functions related to environmental stress tolerance, adhesion processes, multidrug efflux systems, and heavy metal resistance. Genome annotation identified the various genes for chemotaxis, flagellar motility, and biofilm production, illustrating its strong colonization ability. CONCLUSION: The current findings provides the in-depth investigation of a probiotic Paenibacillus bacterium that possessed various genome features that enable the bacterium to survive under diverse conditions. The strain shows the strong ability for probiotic application purposes.

First Report on the Versatile Secretome of an Environmental Isolate Acinetobacter pittii S-30.

Acinetobacter species is currently ranked as high-priority pathogen for their extraordinary ability to become resistant to almost all existing antibiotics. The diverse range of effectors secreted by Acinetobacter spp. constitutes a significant proportion of the virulence arsenal. Therefore, our study aims to characterize the secretome of Acinetobacter pittii S-30. Analysis of extracellular secreted proteins of A. pittii S-30 revealed the presence of transporter proteins, outer membrane proteins, molecular chaperones, porins, and several proteins of unknown function. Additionally, proteins related to metabolism, as well as those involved in gene expression and protein translation, type VI secretion system (T6SS) proteins, and stress response-related proteins were also identified in the secretome. The comprehensive analysis of secretome revealed putative protein antigens which could elicit substantial immune response. The limited availability of effective antibiotics and the worldwide growth of secretome data make this approach appealing in the development of effective vaccines against Acinetobacter and other bacterial pathogens.

Secretome of Paenibacillus sp. S-12 provides an insight about its survival and possible pathogenicity.

Our aim in this study was to characterize and investigate the secretome of Paenibacillus sp. S-12 by nanoLC-MS/MS tool-based analysis of trypsin digested culture supernatant proteins. Using a bioinformatics and combined approach of mass spectrometry, we identified 657 proteins in the secretome. Bioinformatic tools such as PREDLIPO, SecretomeP 2.0, SignalP 4.1, and PSORTb were used for the subcellular localization and categorization of secretome on basis of signal peptides. Among the identified proteins, more than 25% of the secretome proteins were associated with virulence proteins including flagellar, adherence, and immune modulators. Gene ontology analysis using Blast2GO tools categorized 60 proteins of the secretome into biological processes, cellular components, and molecular functions. KEGG pathway analysis identified the enzymes or proteins involved in various biosynthesis and degradation pathways. Functional analysis of secretomes reveals a large number of proteins involved in the uptake and exchange of nutrients, colonization, and chemotaxis. A good number of proteins were involved in survival and defense mechanism against oxidative stress, the production of toxins and antimicrobial compounds. The present study is the first report of the in-depth protein profiling of Paenibacillus bacterium. In summary, the current findings of Paenibacillus sp. S-12 secretome provide basic information to understand its survival and the possible pathogenic mechanism.

Unraveling the Virulence Factors and Secreted Proteins of an Environmental Isolate Enterobacter sp. S-16.

Members of the Enterobacter genus include many pathogenic microbes of humans and plants, secrete proteins that contribute to the interactions of bacteria and their environment. Therefore, understanding of secreted proteins is vital to understand bacterial physiology and behavior. Here, we explored the secretome of an environmental isolate Enterobacter sp. S-16 by nanoLC-MS/MS and identified 572 proteins in the culture supernatant. Gene ontology (GO) analysis indicated that proteins were related to biological processes, molecular as well as cellular functions. The majority of the identified proteins are involved in microbial metabolism, chemotaxis & motility, flagellar hook-associated proteins, biosynthesis of antibiotics, and molecular chaperones to assist the protein folding. Bioinformatics analysis of the secretome revealed the presence of type I and type VI secretion system proteins. Presence of these diverse secretion system proteins in Enterobacter sp. S-16 are likely to be involved in the transport of various proteins including nutrient acquisition, adhesion, colonization, and homeostasis maintenance. Among the secreted bacterial proteins with industrial importance, lignocellulolytic enzymes play a major role, therefore, we analyzed our secretome results for any presence of glycoside hydrolases (GHs) and other hydrolytic enzymes (CAZymes). Overall, the secreted proteins may be considered an attractive reservoir of potential antigens for drug development, diagnostic markers, and other biomedical applications.

Secretome analysis of an environmental isolate Enterobacter sp. S-33 identifies proteins related to pathogenicity.

Enterobacter species are responsible for causing infections of the lower respiratory tract, urinary tract, meninges, etc. Proteins secreted by these species may act as determinants of host-pathogen interaction and play a role in virulence. Among the secreted proteins, the Type VI secretion system (T6SS) acts as a molecular nanomachine to deliver many effector proteins directly into prey cells in a contact-dependent manner. The secreted proteins may provide an idea for the interaction of bacteria to their environment and an understanding of the role of these proteins for their role in bacterial physiology and behaviour. Therefore, aim of this study was to characterize the secreted proteins in the culture supernatant by a T6SS bacterium Enterobacter sp. S-33 using nano-LC-MS/MS tool. Using a combined mass spectrometry and bioinformatics approach, we identified a total of 736 proteins in the secretome. Bioinformatics analysis predicting subcellular localization identified 110 of the secreted proteins possessed signal sequences. By gene ontology analysis, more than 80 proteins of the secretome were classified into biological or molecular functions. More than 20 percent of secretome proteins were virulence proteins including T6SS proteins, proteins involved in adherence and fimbriae formation, molecular chaperones, outer membrane proteins, serine proteases, antimicrobial, biofilm, exotoxins, etc. In summary, the results of the present study of the S-33 secretome provide a basis for understanding the possible pathogenic mechanisms and future investigation by detailed experimental approach will provide a confirmation of secreted virulence proteins in the exact role of virulence using the in vivo model.

Genome wide identification of lncRNAs and circRNAs having regulatory role in fruit shelf life in health crop cucumber (Cucumis sativus L.).

Cucumber is an extremely perishable vegetable; however, under room conditions, the fruits become unfit for consumption 2-3 days after harvesting. One natural variant, DC-48 with an extended shelf-life was identified, fruits of which can be stored up to 10-15 days under room temperature. The genes involved in this economically important trait are regulated by non-coding RNAs. The study aims to identify the long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) by taking two contrasting genotypes, DC-48 and DC-83, at two different fruit developmental stages. The upper epidermis of the fruits was collected at 5 days and 10 days after pollination (DAP) for high throughput RNA sequencing. The differential expression analysis was performed to identify differentially expressed (DE) lncRNAs and circRNAs along with the network analysis of lncRNA, miRNA, circRNA, and mRNA interactions. A total of 97 DElncRNAs were identified where 18 were common under both the developmental stages (8 down regulated and 10 upregulated). Based on the back-spliced reads, 238 circRNAs were found to be distributed uniformly throughout the cucumber genomes with the highest numbers (71) in chromosome 4. The majority of the circRNAs (49%) were exonic in origin followed by inter-genic (47%) and intronic (4%) origin. The genes related to fruit firmness, namely, polygalacturonase, expansin, pectate lyase, and xyloglucan glycosyltransferase were present in the target sites and co-localized networks indicating the role of the lncRNA and circRNAs in their regulation. Genes related to fruit ripening, namely, trehalose-6-phosphate synthase, squamosa promoter binding protein, WRKY domain transcription factors, MADS box proteins, abscisic stress ripening inhibitors, and different classes of heat shock proteins (HSPs) were also found to be regulated by the identified lncRNA and circRNAs. Besides, ethylene biosynthesis and chlorophyll metabolisms were also found to be regulated by DElncRNAs and circRNAs. A total of 17 transcripts were also successfully validated through RT PCR data. These results would help the breeders to identify the complex molecular network and regulatory role of the lncRNAs and circRNAs in determining the shelf-life of cucumbers.

Groundnut Bud Necrosis Virus Modulates the Expression of Innate Immune, Endocytosis, and Cuticle Development-Associated Genes to Circulate and Propagate in Its Vector, Thrips palmi.

Thrips palmi (Thysanoptera: Thripidae) is the predominant tospovirus vector in Asia-Pacific region. It transmits economically damaging groundnut bud necrosis virus (GBNV, family Tospoviridae) in a persistent propagative manner. Thrips serve as the alternate host, and virus reservoirs making tospovirus management very challenging. Insecticides and host plant resistance remain ineffective in managing thrips-tospoviruses. Recent genomic approaches have led to understanding the molecular interactions of thrips-tospoviruses and identifying novel genetic targets. However, most of the studies are limited to Frankliniella species and tomato spotted wilt virus (TSWV). Amidst the limited information available on T. palmi-tospovirus relationships, the present study is the first report of the transcriptome-wide response of T. palmi associated with GBNV infection. The differential expression analyses of the triplicate transcriptome of viruliferous vs. nonviruliferous adult T. palmi identified a total of 2,363 (1,383 upregulated and 980 downregulated) significant transcripts. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses showed the abundance of differentially expressed genes (DEGs) involved in innate immune response, endocytosis, cuticle development, and receptor binding and signaling that mediate the virus invasion and multiplication in the vector system. Also, the gene regulatory network (GRN) of most significant DEGs showed the genes like ABC transporter, cytochrome P450, endocuticle structural glycoprotein, gamma-aminobutyric acid (GABA) receptor, heat shock protein 70, larval and pupal cuticle proteins, nephrin, proline-rich protein, sperm-associated antigen, UHRF1-binding protein, serpin, tyrosine-protein kinase receptor, etc., were enriched with higher degrees of interactions. Further, the expression of the candidate genes in response to GBNV infection was validated in reverse transcriptase-quantitative real-time PCR (RT-qPCR). This study leads to an understanding of molecular interactions between T. palmi and GBNV and suggests potential genetic targets for generic pest control.

Fruit transcriptional profiling of the contrasting genotypes for shelf life reveals the key candidate genes and molecular pathways regulating post-harvest biology in cucumber.

Cucumber fruits are perishable in nature and become unfit for market within 2-3 days of harvesting. A natural variant, DC-48 with exceptionally high shelf life was developed and used to dissect the genetic architecture and molecular mechanism for extended shelf life through RNA-seq for first time. A total of 1364 DEGs were identified and cell wall degradation, chlorophyll and ethylene metabolism related genes played key role. Polygalacturunase (PG), Expansin (EXP) and xyloglucan were down regulated determining fruit firmness and retention of fresh green colour was mainly attributed to the low expression level of the chlorophyll catalytic enzymes (CCEs). Gene regulatory networks revealed the hub genes and cross-talk associated with wide variety of the biological processes. Large number of SSRs (21524), SNPs (545173) and InDels (126252) identified will be instrumental in cucumber improvement. A web genomic resource, CsExSLDb developed will provide a platform for future investigation on cucumber post-harvest biology.

Identification of a major QTL, Parth6.1 associated with parthenocarpic fruit development in slicing cucumber genotype, Pusa Parthenocarpic Cucumber-6.

Parthenocarpy is an extremely important trait that revolutionized the worldwide cultivation of cucumber under protected conditions. Pusa Parthenocarpic Cucumber-6 (PPC-6) is one of the important commercially cultivated varieties under protected conditions in India. Understanding the genetics of parthenocarpy, molecular mapping and the development of molecular markers closely associated with the trait will facilitate the introgression of parthenocarpic traits into non-conventional germplasm and elite varieties. The F1, F2 and back-crosses progenies with a non-parthenocarpic genotype, Pusa Uday indicated a single incomplete dominant gene controlling parthenocarpy in PPC-6. QTL-seq comprising of the early parthenocarpy and non-parthenocarpic bulks along with the parental lines identified two major genomic regions, one each in chromosome 3 and chromosome 6 spanning over a region of 2.7 Mb and 7.8 Mb, respectively. Conventional mapping using F2:3 population also identified two QTLs, Parth6.1 and Parth6.2 in chromosome 6 which indicated the presence of a major effect QTL in chromosome 6 determining parthenocarpy in PPC-6. The flanking markers, SSR01148 and SSR 01012 for Parth6.1 locus and SSR10476 and SSR 19174 for Parth6.2 locus were identified and can be used for introgression of parthenocarpy through the marker-assisted back-crossing programme. Functional annotation of the QTL-region identified two major genes, Csa_6G396640 and Csa_6G405890 designated as probable indole-3-pyruvate monooxygenase YUCCA11 and Auxin response factor 16, respectively associated with auxin biosynthesis as potential candidate genes. Csa_6G396640 showed only one insertion at position 2179 in the non-parthenocarpic parent. In the case of Csa_6G405890, more variations were observed between the two parents in the form of SNPs and InDels. The study provides insight about genomic regions, closely associated markers and possible candidate genes associated with parthenocarpy in PPC-6 which will be instrumental for functional genomics study and better understanding of parthenocarpy in cucumber.

Development of Genome-Wide Functional Markers Using Draft Genome Assembly of Guava (Psidium guajava L.) cv. Allahabad Safeda to Expedite Molecular Breeding.

Guava (Psidium guajava L.), a rich source of nutrients, is an important tropical and subtropical fruit of the Myrtaceae family and exhibits magnificent diversity. Genetic diversity analysis is the first step toward the identification of parents for hybridization, genetic mapping, and molecular breeding in any crop species. A diversity analysis based on whole-genome functional markers increases the chances of identifying genetic associations with agronomically important traits. Therefore, here, we sequenced the genome of guava cv. Allahabad Safeda on an Illumina platform and generated a draft assembly of ~304 MB. The assembly of the Allahabad Safeda genome constituted >37.95% repeat sequences, gene prediction with RNA-seq data as evidence identified 14,115 genes, and BLAST n/r, Interproscan, PfamScan, BLAST2GO, and KEGG annotated 13,957 genes. A comparative protein transcript analysis of tree species revealed the close relatedness of guava with Eucalyptus. Comparative transcriptomics-based SSR/InDel/SNP-PCR ready genome-wide markers in greenish-yellow skinned and white fleshed-Allahabad Safeda to four contrasting cultivars viz apple-color-skinned and white-fleshed-Lalima, greenish-yellow-skinned and pink-fleshed-Punjab Pink, purple-black-skinned and purple-fleshed-Purple Local and widely used rootstock-Lucknow-49 were developed. The molecular markers developed here revealed a high level of individual heterozygosity within genotypes in 22 phenotypically diverse guava cultivars. Principal coordinate, STRUCTURE clustering, and neighbor-joining-based genetic diversity analysis identified distinct clusters associated with fruit skin and flesh color. The genome sequencing of guava, functional annotation, comparative transcriptomics-based genome-wide markers, and genetic diversity analysis will expand the knowledge of genomes of climacteric fruits, facilitating trait-based molecular breeding and diversifying the nutritional basket.