Genomic evolution and insights into agronomic trait innovations of Sesamum species.

Sesame is an ancient oilseed crop with high oil content and quality. However, the evolutionary history and genetic mechanisms of its valuable agronomic traits remain unclear. Here, we report chromosome-scale genomes of cultivated sesame (Sesamum indicum L.) and six wild Sesamum species, representing all three karyotypes within this genus. Karyotyping and genome-based phylogenic analysis revealed the evolutionary route of Sesamum species from n = 13 to n = 16 and revealed that allotetraploidization occurred in the wild species Sesamum radiatum. Early divergence of the Sesamum genus (48.5-19.7 million years ago) during the Tertiary period and its ancient phylogenic position within eudicots were observed. Pan-genome analysis revealed 9164 core gene families in the 7 Sesamum species. These families are significantly enriched in various metabolic pathways, including fatty acid (FA) metabolism and FA biosynthesis. Structural variations in SiPT1 and SiDT1 within the phosphatidyl ethanolamine-binding protein gene family lead to the genomic evolution of plant-architecture and inflorescence-development phenotypes in Sesamum. A genome-wide association study (GWAS) of an interspecific population and genome comparisons revealed a long terminal repeat insertion and a sequence deletion in DIR genes of wild Sesamum angustifolium and cultivated sesame, respectively; both variations independently cause high susceptibility to Fusarium wilt disease. A GWAS of 560 sesame accessions combined with an overexpression study confirmed that the NAC1 and PPO genes play an important role in upregulating oil content of sesame. Our study provides high-quality genomic resources for cultivated and wild Sesamum species and insights that can improve molecular breeding strategies for sesame and other oilseed crops.

Genome-wide characterization of the wall-associated kinase-like (WAKL) family in sesame (Sesamum indicum) identifies a SiWAKL6 gene involved in resistance to Macrophomina Phaseolina.

BACKGROUND: Sesame charcoal rot caused by Macrophomina phaseolina is one of the most serious fungal diseases in sesame production, and threatens the yield and quality of sesame. WAKL genes are important in the plant response to biotic stresses by sensing and transmitting external signals to the intracellular receptor. However, there is still a lack about the WAKL gene family and its function in sesame resistance to M. phaseolina. The aim of this study was to interpret the roles of WAKL genes in sesame resistance to M. phaseolina. RESULTS: In this study, a comprehensive study of the WAKL gene family was conducted and 31 WAKL genes were identified in the sesame genome. Tandem duplication events were the main factor in expansion of the SiWAKL gene family. Phylogenetic analysis showed that the sesame SiWAKL gene family was divided into 4 groups. SiWAKL genes exhibited different expression patterns in diverse tissues. Under M. phaseolina stress, most SiWAKL genes were significantly induced. Notably, SiWAKL6 was strongly induced in the resistant variety "Zhengzhi 13". Functional analysis showed that SiWAKL6 was induced by salicylic acid but not methyl jasmonate in sesame. Overexpression of SiWAKL6 in transgenic Arabidopsis thaliana plants enhanced their resistance to M. phaseolina by inducing the expression of genes involved in the salicylic acid signaling pathway and reconstructing reactive oxygen species homeostasis. CONCLUSIONS: Taken together, the results provide a better understanding of functions about SiWAKL gene family and suggest that manipulation of these SiWAKL genes can improve plant resistance to M. phaseolina. The findings contributed to further understanding of functions of SiWAKL genes in plant immunity.

A 4.43-Kb deletion of chromosomal segment containing an ovate family protein confers long capsule in sesame (Sesamum indicum L.).

KEY MESSAGE: A 4.43-Kb structural variation in the sesame genome results in the deletion of the Siofp1 gene and induces the long capsule length trait. Capsule length (CL) has a positive effect on seed weight and yield in various agronomically important species; however, the molecular mechanism underlying long capsule trait regulation in sesame remains unknown. The inheritance analysis showed that long capsule traits (CL > 4.0 cm) were dominant over normal length (average CL = 3.0 cm) and were controlled by a single gene pair. Association mapping with a RIL population and 259 natural sesame germplasm accessions indicated that the target interval was 52,830-730,961 bp of SiChr.10 in sesame. Meanwhile, the structural variation (SV) of the association mapping revealed that only SV_414325 on chromosome 10 was significantly associated with the CL trait, with a P value of 1.1135E-19. SV_414325 represents a 4430-bp deletion from 414,325 to 418,756 bp on SiChr.10, covering Sindi_2155000 (named SiOFP1). In the normal length type, Siofp1 encodes 411 amino acids of the ovate family proteins and is highly expressed in the leaf, stem, bud, and capsule tissues of sesame. In accordance with the transcriptional repressor character, Siofp1 overexpression in transgenic Arabidopsis (T0 and T1 generations) induced a 25-39% greater shortening of silique length than the wild type (P < 0.05), as well as round cauline leaves and short carpels. These results confirm that SiOFP1 plays a key role in regulating CL trait in sesame and other flowering plants. These findings provide a theoretical and material basis for sesame capsule development and high-yield breeding research.

Genome-wide analysis of the U-box E3 ubiquitin ligase family role in drought tolerance in sesame (Sesamum indicum L.).

Plant U-box (PUB) proteins belong to a class of ubiquitin ligases essential in various biological processes. Sesame (Sesamum indicum L.) is an important and worldwide cultivated oilseed crop. However few studies have been conducted to explore the role of PUBs in drought tolerance in sesame. This study identified a total of 56 members of the sesame PUB family (SiPUB) genes distributed unevenly across all 13 chromosomes. Based on phylogenetic analysis, all 56 SiPUB genes were classified into six groups with various structures and motifs. Cis-acting element analysis suggested that the SiPUB genes are involved in response to various stresses including drought. Based on RNA-seq analysis and quantitative real-time PCR, we identified nine SiPUB genes with significantly different expression profiles under drought stress. The expression patterns of six SiPUB genes in root, leaf and stem tissues corroborated the reliability of the RNA-seq datasets. These findings underscore the importance of SiPUB genes in enhancing drought tolerance in sesame plants. Our study provides novel insights into the evolutionary patterns and variations of PUB genes in sesame and lays the foundation for comprehending the functional characteristics of SiPUB genes under drought-induced stress conditions.

Comprehensive and functional analyses reveal the genomic diversity and potential toxicity of Microcystis.

Microcystis is a cyanobacteria that is widely distributed across the world. It has attracted great attention because it produces the hepatotoxin microcystin (MC) that can inhibit eukaryotic protein phosphatases and pose a great risk to animal and human health. Due to the high diversity of morphospecies and genomes, it is still difficult to classify Microcystis species. In this study, we investigated the pangenome of 23 Microcystis strains to detect the genetic diversity and evolutionary dynamics. Microcystis revealed an open pangenome containing 22,009 gene families and exhibited different functional constraints. The core-genome phylogenetic analysis accurately differentiated the toxic and nontoxic strains and could be used as a taxonomic standard at the genetic level. We also investigated the functions of HGT events, of which were mostly conferred from cyanobacteria and closely related species. In order to detect the potential toxicity of Microcystis, we searched and characterized MC biosynthetic gene clusters and other secondary metabolite gene clusters. Our work provides insights into the genetic diversity, evolutionary dynamics, and potential toxicity of Microcystis, which could benefit the species classification and development of new methods for drinking water quality control and management of bloom formation in the future.

Multi-Omics Approach Reveals the Potential Core Vaccine Targets for the Emerging Foodborne Pathogen Campylobacter jejuni.

Campylobacter jejuni is a leading cause of bacterial gastroenteritis in humans around the world. The emergence of bacterial resistance is becoming more serious; therefore, development of new vaccines is considered to be an alternative strategy against drug-resistant pathogen. In this study, we investigated the pangenome of 173 C. jejuni strains and analyzed the phylogenesis and the virulence factor genes. In order to acquire a high-quality pangenome, genomic relatedness was firstly performed with average nucleotide identity (ANI) analyses, and an open pangenome of 8,041 gene families was obtained with the correct taxonomy genomes. Subsequently, the virulence property of the core genome was analyzed and 145 core virulence factor (VF) genes were obtained. Upon functional genomics and immunological analyses, five core VF proteins with high antigenicity were selected as potential core vaccine targets for humans. Furthermore, functional annotations indicated that these proteins are involved in important molecular functions and biological processes, such as adhesion, regulation, and secretion. In addition, transcriptome analysis in human cells and pig intestinal loop proved that these vaccine target genes are important in the virulence of C. jejuni in different hosts. Comprehensive pangenome and relevant animal experiments will facilitate discovering the potential core vaccine targets with improved efficiency in reverse vaccinology. Likewise, this study provided some insights into the genetic polymorphism and phylogeny of C. jejuni and discovered potential vaccine candidates for humans. Prospective development of new vaccines using the targets will be an alternative to the use of antibiotics and prevent the development of multidrug-resistant C. jejuni in humans and even other animals.

Structure and genetics of Escherichia coli O antigens.

Escherichia coli includes clonal groups of both commensal and pathogenic strains, with some of the latter causing serious infectious diseases. O antigen variation is current standard in defining strains for taxonomy and epidemiology, providing the basis for many serotyping schemes for Gram-negative bacteria. This review covers the diversity in E. coli O antigen structures and gene clusters, and the genetic basis for the structural diversity. Of the 187 formally defined O antigens, six (O31, O47, O67, O72, O94 and O122) have since been removed and three (O34, O89 and O144) strains do not produce any O antigen. Therefore, structures are presented for 176 of the 181 E. coli O antigens, some of which include subgroups. Most (93%) of these O antigens are synthesized via the Wzx/Wzy pathway, 11 via the ABC transporter pathway, with O20, O57 and O60 still uncharacterized due to failure to find their O antigen gene clusters. Biosynthetic pathways are given for 38 of the 49 sugars found in E. coli O antigens, and several pairs or groups of the E. coli antigens that have related structures show close relationships of the O antigen gene clusters within clades, thereby highlighting the genetic basis of the evolution of diversity.

The Roles of Two Type VI Secretion Systems in Cronobacter sakazakii ATCC 12868.

The type VI secretion system (T6SS), which has been found in 25% of gram-negative bacteria, is a crucial virulence factor in several pathogens. Although T6SS gene loci have been discovered in Cronobacter species, one of the major opportunistic foodborne pathogens, its function has not been elucidated. In this study, the roles of two phylogenetically distinct T6SS gene clusters in Cronobacter sakazakii ATCC12868 were investigated. Analysis of 138 genome sequences of C. sakazakii strains, we found that one T6SS gene cluster (T6SS-1) was ubiquitous in all examined strains, whereas another (T6SS-2) was absent or degenerated in a large proportion of the strains (n = 97). In addition, we confirmed the T6SS-1 antibacterial function through an in-frame deletion in the vasK and hcp genes. Compared with the wild-type strain, the T6SS-2-deficient mutant presented a much stronger colonization of organs when infecting neonatal rats. Thus, we proposed that T6SS-2 plays a role in pathogenic processes. This is the first study to investigate the functions of T6SS in C. sakazakii, and the results will extend our understanding of the pathogenic and phylogenetic characteristics of C. sakazakii.

Establishment of a suspension array for Pseudomonas aeruginosa O-antigen serotyping.

Pseudomonas aeruginosa (P. aeruginosa) is widespread in the environment. It readily infects hospital and immunocompromised patients. Traditional biochemical and immunological diagnoses of P. aeruginosa infection cannot meet clinical demands. The variability of O-antigens is the primary basis for the serotyping schemes of many gram-negative bacteria, which is the most widely used method for pathogenic epidemiological purposes. In this study, we developed a suspension array that can accurately diagnose and identify 19 different P. aeruginosa O-antigen serotypes except O15, whose O-antigen gene cluster has not been characterized. In this assay, wzx/wzy genes were selected as the sero-specific gene for 18 serotypes other than O6, for which the glycosyltransferase gene wbpP was utilized. Meanwhile, the wzyß gene was added to separate O2/16 from O5/18/20, and the insertion sequence (IS) in wzx was used to separate O17 from O11. Eighty-two clinical isolates were screened to test our assay. A total of 65 isolates (79.3%) could be serotyped, and the result were confirmed to be correct by sequencing. Sensitivity analysis indicated that at least 5 ng DNA or 103 CFU cells could be detected using our suspension array. To our knowledge, this is the first report on serotyping P. aeruginosa by suspension array and may be of great value in the clinical diagnostics of P. aeruginosa infection.

Identifying genetic diversity of O antigens in Aeromonas hydrophila for molecular serotype detection.

Aeromonas hydrophila is a globally occurring, potentially virulent, gram-negative opportunistic pathogen that is known to cause water and food-borne diseases around the world. In this study, we use whole genome sequencing and in silico analyses to identify 14 putative O antigen gene clusters (OGCs) located downstream of the housekeeping genes acrB and/or oprM. We have also identified 7 novel OGCs by analyzing 15 publicly available genomes of different A. hydrophila strains. From the 14 OGCs identified initially, we have deduced that O antigen processing genes involved in the wzx/wzy pathway and the ABC transporter (wzm/wzt) pathway exhibit high molecular diversity among different A. hydrophila strains. Using these genes, we have developed a multiplexed Luminex-based array system that can identify up to 14 A. hydrophila strains. By combining our other results and including the sequences of processing genes from 13 other OGCs (7 OGCs identified from publicly available genome sequences and 6 OGCs that were previously published), we also have the data to create an array system that can identify 25 different A. hydrophila serotypes. Although clinical detection, epidemiological surveillance, and tracing of pathogenic bacteria are typically done using serotyping methods that rely on identifying bacterial surface O antigens through agglutination reactions with antisera, molecular methods such as the one we have developed may be quicker and more cost effective. Our assay shows high specificity, reproducibility, and sensitivity, being able to classify A. hydrophila strains using just 0.1 ng of genomic DNA. In conclusion, our findings indicate that a molecular serotyping system for A. hydrophila could be developed based on specific genes, providing an important molecular tool for the identification of A. hydrophila serotypes.

Global regulatory function of the low oxygen-induced transcriptional regulator LoiA in Salmonella Typhimurium revealed by RNA sequencing.

Salmonella enterica serovar Typhimurium (S. Typhimurium) is a major intestinal pathogen that can infect both humans and a variety of animals. LoiA, a novel virulence-regulating protein encoded in Salmonella pathogenicity island (SPI)-14, has been shown to be induced under low oxygen conditions and contribute to S. Typhimurium invasion into intestinal epithetical cells by activating the SPI-1 invasion genes. However, the global regulatory network of LoiA remains unknown. Here, we used high-throughput RNA sequencing (RNA-seq) technology to investigate the regulatory function of LoiA in S. Typhimurium under low oxygen conditions. A total of 1250 genes were differentially expressed between the loiA mutant and the wild-type strain; 413 genes were up-regulated and 837 were down-regulated. SPI-1 gene expression was down-regulated in the loiA mutant, consistent with previous results. SPI-2 gene expression was not affected by deletion of loiA; the expression of most genes involved in flagellar basal body and hook biosynthesis was up-regulated in the loiA mutant, while the expression of genes associated with flagellin, motility, and chemotaxis was down-regulated; the expression of lon, encoding an ATP-dependent protease, was up-regulated in the mutant. This study indicates that LoiA regulates a variety of virulence-associated genes in S. Typhimurium. The negative regulation of Lon protease by LoiA indicates that LoiA can regulates several virulence-associated genes in S. Typhimurium via the Lon protease.

Development of rapid and simple experimental and in silico serotyping systems for Citrobacter.

AIM: Members of the genus Citrobacter are important opportunistic pathogens responsible for high mortality rate. Therefore, in this study, we aimed to develop efficient and accurate Citrobacter typing schemes for clinical detection and epidemiological surveillance. MATERIALS & METHODS: Using genomic and experimental analyses, we located the O-antigen biosynthesis gene clusters in Citrobacter genome for the first time, and used comparative genomic analyses to reveal the specific genes in different Citrobacter serotypes. RESULTS: Based on the specific genes in O-antigen biosynthesis gene clusters of Citrobacter, we established experimental and in silico serotyping systems for this bacterium. CONCLUSION: Both serotyping tools are reliable, and our observations are biologically and clinically relevant for understanding and managing Citrobacter infection.

Development of a molecular serotyping scheme and a multiplexed luminex-based array for Providencia.

Providencia is an opportunistic human pathogen that belongs to the Enterobacteriaceae family. The bacterial cell surface O-antigen is one of the most structurally variable cell constituents and serves as a basis for serotyping gram-negative bacteria. In this work, the genomes of 12 Providencia strains were sequenced, and genes driving O-antigen biosynthesis were analyzed. The O-antigen-synthesizing genes of Providencia are located in the O-antigen gene cluster (OGC) between the cpxA and yibK genes. The gene functions predicted in silico agreed with the known O-antigen structures. All clusters were found to contain both wzx and wzy and exhibit a high degree of heterogeneity. Based on the sero-specific genes, we developed a molecular serotyping system to detect 23 serotypes (from the present and previous studies) for the first time. Five Proteus strains, five Morganella strains, five uropathogenic Escherichia coli (UPEC) strains and 32 Providencia strains with other serotypes were used to assess the specificity of our multiplexed Luminex-based array. Five serogroups (O3, O8, O19, O38 and O52 strains) were used to determine the sensitivity of the suspension array. The detection sensitivity was 0.1 ng genomic DNA, 103 CFU/ml in pure culture, or 104 CFU/ml in mock urine specimens. Furthermore, 29 publicly available Providencia genomes (which have not been serotyped) were analyzed, and 23 novel putative OGC types were identified. In total, we identified 35 new OGCs and developed a molecular serotyping system based on the sero-specific genes. The established classification system can support promising applications in basic research, clinical diagnosis, and epidemiological surveillance.

RNA-seq-based monitoring of gene expression changes of viable but non-culturable state of Vibrio cholerae induced by cold seawater.

Vibrio cholerae O1 is a natural inhabitant of aquatic environments and causes the acute diarrheal disease cholera. Entry into a viable but non-culturable (VBNC) state is a survival strategy by which V. cholerae withstands natural stresses and is important for the transition between the aquatic and host environments during the V. cholerae life cycle. In this study, the formation of VBNC V. cholerae induced by cold seawater exposure was investigated using RNA sequencing (RNA-seq). The analysis revealed that the expression of 1420 genes was changed on VBNC state formation. In the VBNC cells, genes related to biofilm formation, chitin utilization and stress responses were upregulated, whereas those related to cell division, morphology and ribosomal activity were mainly downregulated. The concurrent acquisition of a carbon source and the arrest of cell division in cells with low metabolic activity help bacteria increase their resistance to unfavourable environments. Moreover, two transcriptional regulators, SlmA and MetJ, were found to play roles in both VBNC formation and intestinal colonization, suggesting that some genes may function in both processes. This acquired knowledge will improve our understanding of the molecular mechanisms of stress tolerance and may help control future cholera infections and outbreaks.

Genetic diversity of the O antigens of Proteus species and the development of a suspension array for molecular serotyping.

Proteus species are well-known opportunistic pathogens frequently associated with skin wound and urinary tract infections in humans and animals. O antigen diversity is important for bacteria to adapt to different hosts and environments, and has been used to identify serotypes of Proteus isolates. At present, 80 Proteus O-serotypes have been reported. Although the O antigen structures of most Proteus serotypes have been identified, the genetic features of these O antigens have not been well characterized. The O antigen gene clusters of Proteus species are located between the cpxA and secB genes. In this study, we identified 55 O antigen gene clusters of different Proteus serotypes. All clusters contain both the wzx and wzy genes and exhibit a high degree of heterogeneity. Potential functions of O antigen-related genes were proposed based on their similarity to genes in available databases. The O antigen gene clusters and structures were compared, and a number of glycosyltransferases were assigned to glycosidic linkages. In addition, an O serotype-specific suspension array was developed for detecting 31 Proteus serotypes frequently isolated from clinical specimens. To our knowledge, this is the first comprehensive report to describe the genetic features of Proteus O antigens and to develop a molecular technique to identify different Proteus serotypes.