BackgroundMitogen-activated protein kinase kinase kinases (MAPKKKs) are the important components of MAPK cascades, which play the crucial role in plant growth and development as well as in response to diverse stresses. Although this family has been systematically studied in many plant species, little is known about MAPKKK genes in wheat (Triticum aestivum L.), especially those involved in the regulatory network of stress processes.ResultsIn this study, we identified 155 wheat MAPKKK genes through a genome-wide search method based on the latest available wheat genome information, of which 29 belonged to MEKK, 11 to ZIK and 115 to Raf subfamily, respectively. Then, chromosome localization, gene structure and conserved protein motifs and phylogenetic relationship as well as regulatory network of these TaMAPKKKs were systematically investigated and results supported the prediction. Furthermore, a total of 11 homologous groups between A, B and D sub-genome and 24 duplication pairs among them were detected, which contributed to the expansion of wheat MAPKKK gene family. Finally, the expression profiles of these MAPKKKs during development and under different abiotic stresses were investigated using the RNA-seq data. Additionally, 10 tissue-specific and 4 salt-responsive TaMAPKKK genes were selected to validate their expression level through qRT-PCR analysis.ConclusionsThis study for the first time reported the genome organization, evolutionary features and expression profiles of the wheat MAPKKK gene family, which laid the foundation for further functional analysis of wheat MAPKKK genes, and contributed to better understanding the roles and regulatory mechanism of MAPKKKs in wheat.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-2993-7) contains supplementary material, which is available to authorized users.

Background Klebsiella pneumoniae, as a global priority pathogen, is well known for its capability of acquiring mobile genetic elements that carry resistance and/or virulence genes. Its virulence plasmid, previously deemed nonconjugative and restricted within hypervirulent K. pneumoniae (hvKP), has disseminated into classic K. pneumoniae (cKP), particularly carbapenem-resistant K. pneumoniae (CRKP), which poses alarming challenges to public health. However, the mechanism underlying its transfer from hvKP to CRKP is unclear. Methods A total of 28 sequence type (ST) 11 bloodstream infection-causing CRKP strains were collected from Ruijin Hospital in Shanghai, China, and used as recipients in conjugation assays. Transconjugants obtained from conjugation assays were confirmed by XbaI and S1 nuclease pulsed-field gel electrophoresis, PCR detection and/or whole-genome sequencing. The plasmid stability of the transconjugants was evaluated by serial culture. Genetically modified strains and constructed mimic virulence plasmids were employed to investigate the mechanisms underlying mobilization. The level of extracellular polysaccharides was measured by mucoviscosity assays and uronic acid quantification. An in silico analysis of 2608 plasmids derived from 814 completely sequenced K. pneumoniae strains available in GenBank was performed to investigate the distribution of putative helper plasmids and mobilizable virulence plasmids. Results A nonconjugative virulence plasmid was mobilized by the conjugative plasmid belonging to incompatibility group F (IncF) from the hvKP strain into ST11 CRKP strains under low extracellular polysaccharide-producing conditions or by employing intermediate E. coli strains. The virulence plasmid was mobilized via four modes: transfer alone, cotransfer with the conjugative IncF plasmid, hybrid plasmid formation due to two rounds of single-strand exchanges at specific 28-bp fusion sites or homologous recombination. According to the in silico analysis, 31.8% (242) of the putative helper plasmids and 98.8% (84/85) of the virulence plasmids carry the 28-bp fusion site. All virulence plasmids carry the origin of the transfer site. Conclusions The nonconjugative virulence plasmid in ST11 CRKP strains is putatively mobilized from hvKP or E. coli intermediates with the help of conjugative IncF plasmids. Our findings emphasize the importance of raising public awareness of the rapid dissemination of virulence plasmids and the consistent emergence of hypervirulent carbapenem-resistant K. pneumoniae (hv-CRKP) strains.

BackgroundWRKY genes, as the most pivotal transcription factors in plants, play the indispensable roles in regulating various physiological processes, including plant growth and development as well as in response to stresses. Broomcorn millet is one of the most important crops in drought areas worldwide. However, the WRKY gene family in broomcorn millet remains unknown.ResultsA total of 32 PmWRKY genes were identified in this study using computational prediction method. Structural analysis found that PmWRKY proteins contained a highly conserved motif WRKYGQK and two common variant motifs, namely WRKYGKK and WRKYGEK. Phylogenetic analysis of PmWRKYs together with the homologous genes from the representative species could classify them into three groups, with the number of 1, 15, and 16, respectively. Finally, the transcriptional profiles of these 32 PmWRKY genes in various tissues or under different abiotic stresses were systematically investigated using qRT-PCR analysis. Results showed that the expression level of 22 PmWRKY genes varied significantly under one or more abiotic stress treatments, which could be defined as abiotic stress-responsive genes.ConclusionsThis was the first study to identify the organization and transcriptional profiles of PmWRKY genes, which not only facilitates the functional analysis of the PmWRKY genes, and also lays the foundation to reveal the molecular mechanism of stress tolerance in this important crop.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-2677-3) contains supplementary material, which is available to authorized users.

The homeodomain-leucine zipper (HD-Zip) gene family, as plant-specific transcription factors, plays an important role in plant development and growth as well as in the response to diverse stresses. Although HD-Zip genes have been extensively studied in many plants, they had not yet been studied in wheat, especially those involved in response to abiotic stresses. In this study, 46 wheat HD-Zip genes were identified using a genome-wide search method. Phylogenetic analysis classified these genes into four groups, numbered 4, 5, 17 and 20 respectively. In total, only three genes with A, B and D homoeologous copies were identified. Furthermore, the gene interaction networks found that the TaHDZ genes played a critical role in the regulatory pathway of organ development and osmotic stress. Finally, the expression profiles of the wheat HD-Zips in different tissues and under various abiotic stresses were investigated using the available RNA sequencing (RNA-Seq) data and then validated by quantitative real-time polymerase chain reaction (qRT-PCR) to obtain the tissue-specific and stress-responsive candidates. This study systematically identifies the HD-Zip gene family in wheat at the genome-wide level, providing important candidates for further functional analysis and contributing to the better understanding of the molecular basis of development and stress tolerance in wheat.

VRprofile2 is an updated pipeline that rapidly identifies diverse mobile genetic elements in bacterial genome sequences. Compared with the previous version, three major improvements were made. First, the user-friendly visualization could aid users in investigating the antibiotic resistance gene cassettes in conjunction with various mobile elements in the multiple resistance region with mosaic structure. VRprofile2 could compare the predicted mobile elements to the collected known mobile elements with similar architecture. A new mobilome indicator was proposed to give an overall estimation of the mobilome size in individual bacterial genomes. Second, the relationship between antibiotic resistance genes, mobile elements, and host strains would be efficiently examined with the aid of predicted strain's sequence typing, the incompatibility group and the transferability of plasmids. Finally, the updated back-end database, MobilomeDB2, now collected nearly a thousand active mobile elements retrieved from literature or based on prediction. The pre-computed results of the antibiotic resistance gene-carrying mobile elements of >5500 ESKAPEE genomes were also provided. We expect that VRprofile2 will provide better support for researchers interested in bacterial mobile elements and the dissemination of antibiotic resistance. VRprofile2 is freely available to all users without any login requirement at https://tool2-mml.sjtu.edu.cn/VRprofile.

Virulence plasmids of hypervirulent Klebsiella pneumoniae (hvKp) have the potential to transfer to drug-resistant strains or integrate with other plasmids, facilitating the genome evolution of threatening pathogens. We conducted an in-depth analysis of the publicly available 156 complete genome sequences of hvKp together with a multi-region clinical cohort of 171 hvKp strains from China to provide evidence for the virulence plasmid evolution. Virulence plasmids were frequently detected in the ST23 and ST11 K. pneumoniae strains. Multidrug-resistant hvKp (MDR-hvKp) occupied a large proportion of hvKp, and the coexistence of virulence and resistance plasmids may be the major cause. Virulence plasmids commonly possessed multiple replicons, of which IncFIB K was the most prevalent (84.6%). We identified 49 IncFIB K alleles among 583 IncFIB K plasmids, and they could be divided into Clades I, II, and III. We further observed that conjugative and non-conjugative virulence plasmids could be distinguished by IncFIB K genetic diversity, and IncFIB K subtyping could also indirectly indicate a chimeric preference of conjugative virulence plasmids. On this basis, we developed an open-access web tool called KpVR for IncFIB K subtyping. In conclusion, the genetic diversity of IncFIB K virulence plasmids could be used for tracking the evolution of virulence plasmids, and further preventing the emergence of MDR-hvKp strains.

Microsatellites are an important constituent of plant genome and distributed across entire genome. In this study, genome-wide analysis of microsatellites in 8 Triticeae species and 9 model plants revealed that microsatellite characteristics were similar among the Triticeae species. Furthermore, genome-wide microsatellite markers were designed in wheat and then used to analyze the evolutionary relationship of wheat and other Triticeae species. Results displayed that Aegilops tauschii was found to be the closest species to Triticum aestivum, followed by Triticum urartu, Triticum turgidum and Aegilops speltoides, while Triticum monococcum, Aegilops sharonensis and Hordeum vulgare showed a relatively lower PCR amplification effectivity. Additionally, a significantly higher PCR amplification effectivity was found in chromosomes at the same subgenome than its homoeologous when these markers were subjected to search against different chromosomes in wheat. After a rigorous screening process, a total of 20,666 markers showed high amplification and polymorphic potential in wheat and its relatives, which were integrated with the public available wheat markers and then anchored to the genome of wheat (CS). This study not only provided the useful resource for SSR markers development in Triticeae species, but also shed light on the evolution of polyploid wheat from the perspective of microsatellites.

The mitogen-activated protein kinase (MAPK) cascade is a universal signal transduction module that plays a vital role in regulating growth and development, as well as environmental stress responses in plants. Wheat is one of the most important crops worldwide. Although the MAPK kinase kinase (MAP3K) family in wheat has been investigated, the MAPK and MAPK kinase (MAP2K) gene families remain unknown at present. Here, 54 MAPK and 18 MAPKK genes were identified in wheat using recent genomic information. Phylogenetic analysis of Triticum aestivum L. MAPKs and MAPKKs (TaMAPKs and TaMAPKKs) together with homologous genes from other species classified them into four groups, and the clustering was consistent with the genomic exon/intron structures. Conserved motif analysis found that MAPK proteins contained a typical TXY phosphorylation site and MAPKK proteins contained an S/T-X5-S/T motif. RNA-seq data mapping analysis showed that MAPK and MAPKK genes in group IV had tissue-specific expression profiles, whereas each group I member showed relatively high expression in all organs. Expression patterns of TaMAPK and TaMAPKK genes under stress conditions were also investigated and stress-responsive candidates were identified. Co-expression network analysis identified 11 TaMAPK genes and 6 TaMAPKK genes involved in the interaction network pathway. Overall, this study provided useful information for evolutionary and functional surveys of MAPK and MAPKK gene families in wheat and beyond.