Long read and single molecule DNA sequencing simplifies genome assembly and TAL effector gene analysis of Xanthomonas translucens

High genomic plasticity and unique features of Xanthomonas translucens pv. graminis revealed through comparative analysis of complete genome sequences

BACKGROUND

Xanthomonas translucens pv. graminis (Xtg) is a major bacterial pathogen of economically important forage grasses, causing severe yield losses. So far, genomic resources for this pathovar consisted mostly of draft genome sequences, and only one complete genome sequence was available, preventing comprehensive comparative genomic analyses. Such comparative analyses are essential in understanding the mechanisms involved in the virulence of pathogens and to identify virulence factors involved in pathogenicity.

RESULTS

In this study, we produced high-quality, complete genome sequences of four strains of Xtg, complementing the recently obtained complete genome sequence of the Xtg pathotype strain. These genomic resources allowed for a comprehensive comparative analysis, which revealed a high genomic plasticity with many chromosomal rearrangements, although the strains were highly related. A high number of transposases were exclusively found in Xtg and corresponded to 413 to 457 insertion/excision transposable elements per strain. These mobile genetic elements are likely to be involved in the observed genomic plasticity and may play an important role in the adaptation of Xtg. The pathovar was found to lack a type IV secretion system, and it possessed the smallest set of type III effectors in the species. However, three XopE and XopX family effectors were found, while in the other pathovars of the species two or less were present. Additional genes that were specific to the pathovar were identified, including a unique set of minor pilins of the type IV pilus, 17 TonB-dependent receptors (TBDRs), and 11 plant cell wall degradative enzymes.

CONCLUSION

These results suggest a high adaptability of Xtg, conferred by the abundance of mobile genetic elements, which could play a crucial role in pathogen adaptation. The large amount of such elements in Xtg compared to other pathovars of the species could, at least partially, explain its high virulence and broad host range. Conserved features that were specific to Xtg were identified, and further investigation will help to determine genes that are essential to pathogenicity and host adaptation of Xtg.

Comparative Genomics of Xanthomonas translucens pv. undulosa Strains Isolated from Weedy Grasses and Cultivated Wild Rice

Bacterial leaf streak (BLS) of wheat (Triticum aestivum), caused by Xanthomonas translucens pv. undulosa, is a disease of major concern in the Northern Great Plains. The host range for X. translucens pv. undulosa is relatively broad, including several small grains and perennial grasses. In Minnesota, X. translucens pv. undulosa was isolated from weedy grasses in and around wheat fields that exhibited BLS symptoms and from cultivated wild rice (Zizania palustris) with symptomatic leaf tissue. Currently, no genomic resources are available for X. translucens pv. undulosa strains isolated from non-wheat hosts. In this study, we sequenced and assembled the complete genomes of five strains isolated from weedy grass hosts, foxtail barley (Hordeum jubatum), green foxtail (Setaria viridis), and wild oat (Avena fatua), and from cultivated wild rice and wheat. These five genomes were compared with the publicly available genomes of seven X. translucens pv. undulosa strains originating from wheat and one genome of an X. translucens pv. secalis strain originating from rye (Secale cereale). Global alignments of the genomes revealed little variation in genomic structures. Average nucleotide identity-based phylogeny and life identification numbers revealed that the strains share ≥99.25% identity. We noted differences in the presence of Type III secreted effectors, including transcription activator-like effectors. Despite differences between strains, we did not identify unique features distinguishing strains isolated from wheat and non-wheat hosts. This study contributes to the availability of genomic data for X. translucens pv. undulosa from non-wheat hosts, thus increasing our understanding of the diversity within the pathogen population.

Epidemiology, Host Resistance, and Genomics of the Small Grain Cereals Pathogen Xanthomonas translucens: New Advances and Future Prospects

Bacterial leaf streak (BLS) primarily affects barley and wheat and is mainly caused by the pathogens Xanthomonas translucens pv. translucens and X. translucens pv. undulosa, respectively. BLS is distributed globally and poses a risk to food security and the supply of malting barley. X. translucens pv. cerealis can infect both wheat and barley but is rarely isolated from these hosts in natural infections. These pathogens have undergone a confusing taxonomic history, and the biology has been poorly understood, making it difficult to develop effective control measures. Recent advancements in the ability and accessibility to sequence bacterial genomes have shed light on phylogenetic relationships between strains and identified genes that may play a role in virulence, such as those that encode Type III effectors. In addition, sources of resistance to BLS have been identified in barley and wheat lines, and ongoing efforts are being made to map these genes and evaluate germplasm. Although there are still gaps in BLS research, progress has been made in recent years to further understand epidemiology, diagnostics, pathogen virulence, and host resistance.

Detection and Diagnosis of Bacterial Leaf Streak on Small Grain Cereals: From Laboratory to Field

Bacterial leaf streak of small-grain cereals is an economically important disease of wheat and barley crops. The disease occurs in many countries across the globe, with particular importance in regions characterized by high precipitation or areas in which sprinkler irrigation is used. Three genetically distinct lineages of the Gram-negative bacterium Xanthomonas translucens (X. translucens pv. undulosa, X. translucens pv. translucens, and X. translucens pv. cerealis) are responsible for most of the bacterial leaf streak infections on wheat and barley crops. Considering the seedborne nature of the pathogens, they are included in the A2 (high-risk) list of quarantine organisms for some European countries; hence, they are under strict quarantine control and zero tolerance. Due to the taxonomic complexities within X. translucens, the exact geographic distribution of each pathovar has not yet been determined. In this mini review, we provide an updated overview of the detection and diagnosis of the bacterial leaf streak pathogens. First, a short history of the leaf streak pathogens is provided, followed by the symptomology and host range of the causal agents. Then, the utility of conventional methods and high-throughput molecular approaches in the precise detection and identification of the pathogens is explained. Finally, we highlight the role of quarantine inspections and early detection of the pathogen in combating the risk of bacterial leaf streak in the 21st century's small-grains cereals' industry.

Comparative Genomics of Barley-Infecting Xanthomonas translucens Shows Overall Genetic Similarity but Globally Distributed Virulence Factor Diversity

Xanthomonas translucens pv. translucens (Xtt) is a global barley patho-gen and a concern for resistance breeding and regulation. Long-read whole genome sequences allow in-depth understanding of pathogen diversity. We have completed long-read PacBio sequencing of two Minnesotan Xtt strains and an in-depth analysis of available Xtt genomes. We found that average nucleotide identity (ANI)-based approaches organize Xtt strains different from the previous standard multilocus sequencing analysis approach. According to ANI, Xtt forms a separate clade from X. translucens pv. undulosa and consists of three main groups which are represented on multiple continents. Some virulence factors, such as 17 Type III-secreted effectors, are highly conserved and offer potential targets for the elicitation of broad resistance. However, there is a high degree of variation in virulence factors, meaning that germplasm should be screened for resistance with a diverse panel of Xtt.

DNA Markers, Pathogenicity Test, and Multilocus Sequence Analysis to Differentiate and Characterize Cereal-Specific Xanthomonas translucens Strains

Xanthomonas translucens contains a group of bacterial pathogens that are closely related and have been divided into several pathovars based on their host range. X. translucens pv. undulosa (Xtu) and X. translucens pv. translucens (Xtt) are two important pathovars that cause bacterial leaf streak disease on wheat and barley, respectively. In this study, DNA markers were developed to differentiate Xtu and Xtt and were then used to characterize a collection of X. translucens strains with diverse origins, followed by confirmation and characterization with pathogenicity tests and multilocus sequence analysis/typing (MLSA/MLST). We first developed cleaved amplified polymorphic sequence markers based on the single-nucleotide polymorphisms within a cereal pathovar-specific DNA sequence. In addition, two Xtt-specific markers, designated Xtt-XopM and Xtt-SP1, were developed from comparative genomics among the sequenced Xtt/Xtu genomes. Using the developed markers, a collection of X. translucens strains were successfully identified as Xtu or Xtt. Pathogenicity tests on wheat and barley plants and MLSA of four housekeeping genes validated the pathovar assignation of those strains. Furthermore, MLSA revealed distinct subclades within both Xtu and Xtt groups. Seven and three sequence types were identified from MLST for Xtu and Xtt strains, respectively. The establishment of efficient Xtt/Xtu differentiation methods and characterization of those strains will be useful in studying disease epidemiology and host-pathogen interactions and breeding programs when screening for sources of resistance for these two important bacterial pathogens.

Comparative Transcriptomic Analysis of Wheat Cultivars in Response to Xanthomonas translucens pv. cerealis and Its T2SS, T3SS, and TALEs Deficient Strains

Xanthomonas translucens pv. cerealis causes bacterial leaf streak disease on small grain cereals. Type II and III secretion systems (T2SS and T3SS) play a pivotal role in the pathogenicity of the bacterium, while no data are available on the transcriptomic profile of wheat cultivars infected with either wild type (WT) or mutants of the pathogen. In this study, WT, TAL-effector mutants, and T2SS/T3SS mutants of X. translucens pv. cerealis strain NXtc01 were evaluated for their effect on the transcriptomic profile of two wheat cultivars, 'Chinese Spring' and 'Yangmai-158', using Illumina RNA-sequencing technology. RNA-Seq data showed that the number of differentially expressed genes (DEGs) was higher in Yangmai-158 than in Chinese Spring, suggesting higher susceptibility of Yangmai-158 to the pathogen. In T2SS, most suppressed DEGs were related to transferase, synthase, oxidase, WRKY, and bHLH transcription factors. The gspD mutants showed significantly decreased disease development in wheat, suggesting an active contribution of T2SS in virulence. Moreover, the gspD mutant restored full virulence and its multiplication in planta by addition of gspD in trans. In the T3SS-deficient strain, downregulated DEGs were associated with cytochrome, peroxidases, kinases, phosphatases, WRKY, and ethylene-responsive transcription factors. In contrast, upregulated DEGs were trypsin inhibitors, cell number regulators, and calcium transporter. Transcriptomic analyses coupled with quantitative real-time-PCR indicated that some genes are upregulated in Δtal1/Δtal2 compared with the tal-free strain, but no direct interaction was observed. These results provide novel insight into wheat transcriptomes in response to X. translucens infection and pave the way for understanding host-pathogen interactions.

Bacterium-enabled transient gene activation by artificial transcription factors for resolving gene regulation in maize

Understanding gene regulatory networks is essential to elucidate developmental processes and environmental responses. Here, we studied regulation of a maize (Zea mays) transcription factor gene using designer transcription activator-like effectors (dTALes), which are synthetic Type III TALes of the bacterial genus Xanthomonas and serve as inducers of disease susceptibility gene transcription in host cells. The maize pathogen Xanthomonas vasicola pv. vasculorum was used to introduce 2 independent dTALes into maize cells to induced expression of the gene glossy3 (gl3), which encodes a MYB transcription factor involved in biosynthesis of cuticular wax. RNA-seq analysis of leaf samples identified, in addition to gl3, 146 genes altered in expression by the 2 dTALes. Nine of the 10 genes known to be involved in cuticular wax biosynthesis were upregulated by at least 1 of the 2 dTALes. A gene previously unknown to be associated with gl3, Zm00001d017418, which encodes aldehyde dehydrogenase, was also expressed in a dTALe-dependent manner. A chemically induced mutant and a CRISPR-Cas9 mutant of Zm00001d017418 both exhibited glossy leaf phenotypes, indicating that Zm00001d017418 is involved in biosynthesis of cuticular waxes. Bacterial protein delivery of dTALes proved to be a straightforward and practical approach for the analysis and discovery of pathway-specific genes in maize.

Pathovar-Specific PCR Method for Detection and Identification of Xanthomonas translucens pv. undulosa

Bacterial leaf streak disease caused by Xanthomonas translucens pv. undulosa is an economically important disease threatening wheat and barley crops around the globe. Thus far, specific PCR-based detection and identification tests for X. translucens pathovars are not available. In this study, we used comparative genomics approach to design a pathovar-specific primer pair for detection of X. translucens pv. undulosa in naturally infected seeds and its differentiation from other pathovars of the species. For this aim, complete genome sequences of strains of different X. translucens pathovars were compared and the specific PCR primer pair XtuF/XtuR was designed. These primers were strictly specific to X. translucens pv. undulosa because the expected 229-bp DNA fragment was not amplified in the closely related pathovars or in other xanthomonads, wheat-pathogenic bacteria, and other plant-pathogenic bacteria. High sensitivity of the primer pair XtuF/XtuR allowed detection of pure DNA of the pathogen in a concentration as low as 4.5 pg/μl. The pathogen was also detected in water suspension at a concentration of 8.6 × 102 CFU/ml. The PCR test was capable of detecting the pathogen in extracts of naturally infected wheat seeds at a concentration of 3.5 × 104 CFU/g while a culture-plate method was able to detect the pathogen at a concentration of 50 × 105 CFU/g of the same seeds. The PCR test developed in this study is a step forward for precise detection and identification of X. translucens pv. undulosa to prevent outbreaks of the bacterial leaf streak disease.

Pangenome insights into the diversification and disease specificity of worldwide Xanthomonas outbreaks

The bacterial genus Xanthomonas is responsible for disease outbreaks in several hundred plant species, many of them economically important crops. In the era of next-generation sequencing, thousands of strains from this genus have now been sequenced as part of isolated studies that focus on outbreak characterization, host range, diversity, and virulence factor identification. However, these data have not been synthesized and we lack a comprehensive phylogeny for the genus, with some species designations in public databases still relying on phenotypic similarities and representative sequence typing. The extent of genetic cohesiveness among Xanthomonas strains, the distribution of virulence factors across strains, and the impact of evolutionary history on host range across the genus are also poorly understood. In this study, we present a pangenome analysis of 1,910 diverse Xanthomonas genomes, highlighting their evolutionary relationships, the distribution of virulence-associated genes across strains, and rates of horizontal gene transfer. We find a number of broadly conserved classes of virulence factors and considerable diversity in the Type 3 Secretion Systems (T3SSs) and Type 3 Secreted Effector (T3SE) repertoires of different Xanthomonas species. We also use these data to re-assign incorrectly classified strains to phylogenetically informed species designations and find evidence of both monophyletic host specificity and convergent evolution of phylogenetically distant strains to the same host. Finally, we explore the role of recombination in maintaining genetic cohesion within the Xanthomonas genus as a result of both ancestral and recent recombination events. Understanding the evolutionary history of Xanthomonas species and the relationship of key virulence factors with host-specificity provides valuable insight into the mechanisms through which Xanthomonas species shift between hosts and will enable us to develop more robust resistance strategies against these highly virulent pathogens.

How do emerging long-read sequencing technologies function in transforming the plant pathology research landscape?

Long-read sequencing technologies are revolutionizing the sequencing and analysis of plant and pathogen genomes and transcriptomes, as well as contributing to emerging areas of interest in plant-pathogen interactions, disease management techniques, and the introduction of new plant varieties or cultivars. Long-read sequencing (LRS) technologies are progressively being implemented to study plants and pathogens of agricultural importance, which have substantial economic effects. The variability and complexity of the genome and transcriptome affect plant growth, development and pathogen responses. Overcoming the limitations of second-generation sequencing, LRS technology has significantly increased the length of a single contiguous read from a few hundred to millions of base pairs. Because of the longer read lengths, new analysis methods and tools have been developed for plant and pathogen genomics and transcriptomics. LRS technologies enable faster, more efficient, and high-throughput ultralong reads, allowing direct sequencing of genomes that would be impossible or difficult to investigate using short-read sequencing approaches. These benefits include genome assembly in repetitive areas, creating more comprehensive and exact genome determinations, assembling full-length transcripts, and detecting DNA and RNA alterations. Furthermore, these technologies allow for the identification of transcriptome diversity, significant structural variation analysis, and direct epigenetic mark detection in plant and pathogen genomic regions. LRS in plant pathology is found efficient for identifying and characterization of effectors in plants as well as known and unknown plant pathogens. In this review, we investigate how these technologies are transforming the landscape of determination and characterization of plant and pathogen genomes and transcriptomes efficiently and accurately. Moreover, we highlight potential areas of interest offered by LRS technologies for future study into plant-pathogen interactions, disease control strategies, and the development of new plant varieties or cultivars.

Characterization of defense responses against bacterial pathogens in duckweeds lacking EDS1

ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1) mediates the induction of defense responses against pathogens in most angiosperms. However, it has recently been shown that a few species have lost EDS1. It is unknown how defense against disease unfolds and evolves in the absence of EDS1. We utilize duckweeds; a collection of aquatic species that lack EDS1, to investigate this question. We established duckweed-Pseudomonas pathosystems and used growth curves and microscopy to characterize pathogen-induced responses. Through comparative genomics and transcriptomics, we show that the copy number of infection-associated genes and the infection-induced transcriptional responses of duckweeds differ from other model species. Pathogen defense in duckweeds has evolved along different trajectories than in other plants, including genomic and transcriptional reprogramming. Specifically, the miAMP1 domain-containing proteins, which are absent in Arabidopsis, showed pathogen responsive upregulation in duckweeds. Despite such divergence between Arabidopsis and duckweed species, we found conservation of upregulation of certain genes and the role of hormones in response to disease. Our work highlights the importance of expanding the pool of model species to study defense responses that have evolved in the plant kingdom independent of EDS1.

Complete Genome Assemblies of All Xanthomonas translucens Pathotype Strains Reveal Three Genetically Distinct Clades

The Xanthomonas translucens species comprises phytopathogenic bacteria that can cause serious damage to cereals and to forage grasses. So far, the genomic resources for X. translucens were limited, which hindered further understanding of the host–pathogen interactions at the molecular level and the development of disease-resistant cultivars. To this end, we complemented the available complete genome sequence of the X. translucens pv. translucens pathotype strain DSM 18974 by sequencing the genomes of all the other 10 X. translucens pathotype strains using PacBio long-read technology and assembled complete genome sequences. Phylogeny based on average nucleotide identity (ANI) revealed three distinct clades within the species, which we propose to classify as clades Xt-I, Xt-II, and Xt-III. In addition to 2,181 core X. translucens genes, a total of 190, 588, and 168 genes were found to be exclusive to each clade, respectively. Moreover, 29 non-transcription activator-like effector (TALE) and 21 TALE type III effector classes were found, and clade- or strain-specific effectors were identified. Further investigation of these genes could help to identify genes that are critically involved in pathogenicity and/or host adaptation, setting the grounds for the development of new resistant cultivars.

Cassava diseases caused by Xanthomonas phaseoli pv. manihotis and Xanthomonas cassavae

Xanthomonas phaseoli pv. manihotis (Xpm) and X. cassavae (Xc) are two bacterial pathogens attacking cassava. Cassava bacterial blight (CBB) is a systemic disease caused by Xpm, which might have dramatic effects on plant growth and crop production. Cassava bacterial necrosis is a nonvascular disease caused by Xc with foliar symptoms similar to CBB, but its impacts on the plant vigour and the crop are limited. In this review, we describe the epidemiology and ecology of the two pathogens, the impacts and management of the diseases, and the main research achievements for each pathosystem. Because Xc data are sparse, our main focus is on Xpm and CBB.

Xanthomonas translucens pv. undulosa Identified on Common Weedy Grasses in Naturally Infected Wheat Fields in Minnesota

Bacterial leaf streak (BLS) of wheat, caused by Xanthomonas translucens pv. undulosa, has been a notable disease in Minnesota wheat fields over the past decade. Potential sources of the pathogen include infested seed and crop debris. Perennial weeds are also considered a possible inoculum source, but no surveys have been conducted to evaluate which X. translucens pathovars are present on weedy grasses that are common in Minnesota wheat fields. Multilocus sequence analysis (MLSA) of four housekeeping genes (rpoD, dnaK, fyuA, and gyrB) was used to identify 77 strains isolated from six weedy grass species, wheat, and barley in and around naturally infected wheat fields in Minnesota. The MLSA phylogeny identified all strains originating from weedy grass species, except smooth brome, as X. translucens pv. undulosa, whereas strains isolated from smooth brome were determined to be X. translucens pv. cerealis. In planta character states corroborated these identifications on a subset of 41 strains, as all strains from weedy grasses caused water-soaking on wheat and barley in greenhouse assays. Multilocus sequence typing was used to evaluate genetic diversity and revealed that sequence types of X. translucens pv. undulosa originating from weedy grass hosts are similar to those found on wheat. This study identifies both annual and perennial poaceous weeds common in Minnesota that harbor X. translucens pv. undulosa and expands our understanding of the diversity of the pathogen population.

Genomics-Enabled Novel Insight Into the Pathovar-Specific Population Structure of the Bacterial Leaf Streak Pathogen Xanthomonas translucens in Small Grain Cereals

The Gram-negative bacterium Xanthomonas translucens infects a wide range of gramineous plants with a notable impact on small grain cereals. However, genomics-informed intra-species population structure and virulence repertories of the pathogen have rarely been investigated. In this study, the complete genome sequences of seven X. translucens strains representing an entire set of genetic diversity of two pathovars X. translucens pv. undulosa and X. translucens pv. translucens is provided and compared with those of seven publicly available complete genomes of the pathogen. Organization of the 25 type III secretion system genes in all the 14 X. translucens strains was exactly the same, while TAL effector genes localized singly or in clusters across four loci in X. translucens pv. translucens and five to six loci in X. translucens pv. undulosa. Beside two previously unreported endogenous plasmids in X. translucens pv. undulosa, and variations in repeat variable diresidue (RVD) of the 14 strains, tal1a of X. translucens pv. translucens strain XtKm8 encode the new RVDs HE and YI which have not previously been reported in xanthomonads. Further, a number of truncated tal genes were predicted among the 14 genomes lacking conserved BamHI site at N-terminus and SphI site at C-terminus. Our data have doubled the number of complete genomes of X. translucens clarifying the population structure and genomics of the pathogen to pave the way in the small grain cereals industry for disease resistance breeding in the 21st century's agriculture.

Functional Annotation of a Full-Length Transcriptome and Identification of Genes Associated with Flower Development in Rhododendronsimsii (Ericaceae)

Rhododendronsimsii is one of the top ten famous flowers in China. Due to its historical value and high aesthetic, it is widely popular among Chinese people. Various colors are important breeding objectives in Rhododendron L. The understanding of the molecular mechanism of flower color formation can provide a theoretical basis for the improvement of flower color in Rhododendron L. To generate the R.simsii transcriptome, PacBio sequencing technology has been used. A total of 833,137 full-length non-chimeric reads were obtained and 726,846 high-quality full-length transcripts were found. Moreover, 40,556 total open reading frames were obtained; of which 36,018 were complete. In gene annotation analyses, 39,411, 18,565, 16,102 and 17,450 transcriptions were allocated to GO, Nr, KEGG and COG databases, correspondingly. To identify long non-coding RNAs (lncRNAs), we utilized four computational methods associated with Protein families (Pfam), Cooperative Data Classification (CPC), Coding Assessing Potential Tool (CPAT) and Coding Non Coding Index (CNCI) databases and observed 6170, 2265, 4084 and 1240 lncRNAs, respectively. Based on the results, most genes were enriched in the flavonoid biosynthetic pathway. The eight key genes on the anthocyanin biosynthetic pathway were further selected and analyzed by qRT-PCR. The F3'H and ANS showed an upward trend in the developmental stages of R. simsii. The highest expression of F3'5'H and FLS in the petal color formation of R. simsii was observed. This research provided a huge number of full-length transcripts, which will help to proceed genetic analyses of R.simsii. native, which is a semi-deciduous shrub.

Novel Xanthomonas Species From the Perennial Ryegrass Seed Microbiome - Assessing the Bioprotection Activity of Non-pathogenic Relatives of Pathogens

The productivity of the Australian dairy industry is underpinned by pasture grasses, and importantly perennial ryegrass. The performance of these pasture grasses is supported by the fungal endophyte Epichloë spp. that has bioprotection activities, however, the broader microbiome is not well characterized. In this study, we characterized a novel bioprotectant Xanthomonas species isolated from perennial ryegrass (Lolium perenne L. cv. Alto). In vitro and in planta bioassays against key fungal pathogens of grasses (Sclerotium rolfsii, Drechslera brizae and Microdochium nivale) indicated strong bioprotection activities. A complete circular chromosome of ∼5.2 Mb was generated for three strains of the novel Xanthomonas sp. Based on the 16S ribosomal RNA gene, the strains were closely related to the plant pathogen Xanthomonas translucens, however, comparative genomics of 22 closely related xanthomonad strains indicated that these strains were a novel species. The comparative genomics analysis also identified two unique gene clusters associated with the production of bioprotectant secondary metabolites including one associated with a novel nonribosomal peptide synthetase and another with a siderophore. The analysis also identified genes associated with an endophytic lifestyle (e.g., Type VI secretion system), while no genes associated with pathogenicity were identified (e.g., Type III secretion system and effectors). Overall, these results indicate that these strains represent a novel, bioactive, non-pathogenic species of the genus Xanthomonas. Strain GW was the designated type strain of this novel Xanthomonas sp.

Mechanistic insights into host adaptation, virulence and epidemiology of the phytopathogen Xanthomonas

Xanthomonas is a well-studied genus of bacterial plant pathogens whose members cause a variety of diseases in economically important crops worldwide. Genomic and functional studies of these phytopathogens have provided significant understanding of microbial-host interactions, bacterial virulence and host adaptation mechanisms including microbial ecology and epidemiology. In addition, several strains of Xanthomonas are important as producers of the extracellular polysaccharide, xanthan, used in the food and pharmaceutical industries. This polymer has also been implicated in several phases of the bacterial disease cycle. In this review, we summarise the current knowledge on the infection strategies and regulatory networks controlling virulence and adaptation mechanisms from Xanthomonas species and discuss the novel opportunities that this body of work has provided for disease control and plant health.

Factorial estimating assembly base errors using k-mer abundance difference (KAD) between short reads and genome assembled sequences

Genome sequences provide genomic maps with a single-base resolution for exploring genetic contents. Sequencing technologies, particularly long reads, have revolutionized genome assemblies for producing highly continuous genome sequences. However, current long-read sequencing technologies generate inaccurate reads that contain many errors. Some errors are retained in assembled sequences, which are typically not completely corrected by using either long reads or more accurate short reads. The issue commonly exists, but few tools are dedicated for computing error rates or determining error locations. In this study, we developed a novel approach, referred to as k-mer abundance difference (KAD), to compare the inferred copy number of each k-mer indicated by short reads and the observed copy number in the assembly. Simple KAD metrics enable to classify k-mers into categories that reflect the quality of the assembly. Specifically, the KAD method can be used to identify base errors and estimate the overall error rate. In addition, sequence insertion and deletion as well as sequence redundancy can also be detected. Collectively, KAD is valuable for quality evaluation of genome assemblies and, potentially, provides a diagnostic tool to aid in precise error correction. KAD software has been developed to facilitate public uses.

Single-molecule long-read sequencing of the full-length transcriptome of Rhododendron lapponicum L

Rhododendron lapponicum L. is a familiar ornamental plant worldwide with important ornamental and economic value. However, a full-length R. lapponicum transcriptome is still lacking. In the present study, we used the Pacific Biosciences single-molecule real-time sequencing technology to generate the R. lapponicum transcriptome. A total of 346,270 full-length non-chimeric reads were generated, from which we obtained 75,002 high-quality full-length transcripts. We identified 55,255 complete open reading frames, 7,140 alternative splicing events and 2,011 long non-coding RNAs. In gene annotation analyses, 71,155, 33,653, 30,359 and 31,749 transcripts were assigned to the Nr, GO, COG and KEGG databases, respectively. Additionally, 3,150 transcription factors were detected. KEGG pathway analysis showed that 96 transcripts were identified coding for the enzymes associated with anthocyanin synthesis. Furthermore, we identified 64,327 simple sequence repeats from 45,319 sequences, and 150 pairs of primers were randomly selected to develop SSR markers. This study provides a large number of full-length transcripts, which will facilitate the further study of the genetics of R. lapponicum.

Identification of a virulence tal gene in the cotton pathogen, Xanthomonas citri pv. malvacearum strain Xss-V2-18

Background

Bacterial blight of cotton (BBC), which is caused by the bacterium Xanthomonas citri pv. malvacearum (Xcm), is a destructive disease in cotton. Transcription activator-like effectors (TALEs), encoded by tal-genes, play critical roles in the pathogenesis of xanthomonads. Characterized strains of cotton pathogenic Xcm harbor 8–12 different tal genes and only one of them is functionally decoded. Further identification of novel tal genes in Xcm strains with virulence contributions are prerequisite to decipher the Xcm-cotton interactions.

Results

In this study, we identified six tal genes in Xss-V₂–18, a highly-virulent strain of Xcm from China, and assessed their role in BBC. RFLP-based Southern hybridization assays indicated that Xss-V₂–18 harbors the six tal genes on a plasmid. The plasmid-encoded tal genes were isolated by cloning BamHI fragments and screening clones by colony hybridization. The tal genes were sequenced by inserting a Tn5 transposon in the DNA encoding the central repeat region (CRR) of each tal gene. Xcm TALome evolutionary relationship based on TALEs CRR revealed relatedness of Xss-V₂–18 to MSCT1 and MS14003 from the United States. However, Tal2 of Xss-V₂–18 differs at two repeat variable diresidues (RVDs) from Tal6 and Tal26 in MSCT1 and MS14003, respectively, inferred functional dissimilarity. The suicide vector pKMS1 was then used to construct tal deletion mutants in Xcm Xss-V₂–18. The mutants were evaluated for pathogenicity in cotton based on symptomology and growth in planta. Four mutants showed attenuated virulence and all contained mutations in tal2. One tal2 mutant designated M2 was further investigated in complementation assays. When tal2 was introduced into Xcm M2 and expressed in trans, the mutant was complemented for both symptoms and growth in planta, thus indicating that tal2 functions as a virulence factor in Xcm Xss-V₂–18.

Conclusions

Overall, the results demonstrated that Tal2 is a major pathogenicity factor in Xcm strain Xss-V₂–18 that contributes significantly in BBC. This study provides a foundation for future efforts aimed at identifying susceptibility genes in cotton that are targeted by Tal2.

High-Quality Genome Resource of Xanthomonas hyacinthi Generated via Long-Read Sequencing

The bacterial plant pathogen Xanthomonas hyacinthi is the causal agent of yellow disease of Hyacinthus and other ornamental plant genera. There is no available complete genome for X. hyacinthi, limiting basic research for this pathogen. Here, we release a high-quality complete genome sequence for the X. hyacinthi type strain, CFBP 1156. Single-molecule real-time (SMRT) sequencing with a mean coverage of 306× revealed two contigs of 4,918,645 and 44,381 bp in size. This was the first characterized plant-disease-causing species of Xanthomonas and this genome provides a resource to better understand the biology of yellow disease of hyacinth.

The translucens group of Xanthomonas translucens: Complicated and important pathogens causing bacterial leaf streak on cereals

Xanthomonas translucens is a group of gram-negative bacteria that can cause important diseases in cereal crops and forage grasses. Different pathovars have been defined according to their host ranges, and molecular and biochemical characteristics. Pathovars have been placed into two major groups: translucens and graminis. The translucens group contains the pathovars causing bacterial leaf streak (BLS) on cereal crops such as wheat, barley, triticale, rye, and oat. In recent years, BLS has re-emerged as a major problem for many wheat- and barley-producing areas worldwide. The biology of the pathogens and the host-pathogen interactions in cereal BLS diseases were poorly understood. However, recent genome sequence data have provided an insight into the bacterial phylogeny and identification and pathogenicity/virulence. Furthermore, identification of sources of resistance to BLS and mapping of the resistance genes have been initiated.

TAXONOMY

Kingdom Bacteria; Phylum Proteobacteria; Class Gammaproteobacteria; Order Xanthomonadales; Family Xanthomonadaceae; Genus Xanthomonas; Species X. translucens; translucens group pathovars: undulosa, translucens, cerealis, hordei, and secalis; graminis group pathovars: arrhenatheri, graminis, poae, phlei; newly established pathovar: pistaciae.

HOST RANGE

X. translucens mainly infects plant species in the Poaceae with the translucens group on cereal crop species and the graminis group on forage grass species. However, some strains have been isolated from, and are able to infect, ornamental asparagus and pistachio trees. Most pathovars have a narrow host range, while a few can infect a broad range of hosts.

GENOME

The complete genome sequence is available for two X. translucens pv. undulosa strains and one pv. translucens strain. A draft genome sequence is also available for at least one strain from each pathovar. The X. translucens pv. undulosa strain Xt4699 was the first to have its complete genome sequenced, which consists of 4,561,137 bp with total GC content approximately at 68% and 3,528 predicted genes.

VIRULENCE MECHANISMS

Like most xanthomonads, X. translucens utilizes a type III secretion system (T3SS) to deliver a suite of T3SS effectors (T3Es) inside plant cells. Transcription activator-like effectors, a special group of T3Es, have been identified in most of the X. translucens genomes, some of which have been implicated in virulence. Genetic factors determining host range virulence have also been identified.

Analysis of Alternative Splicing and Alternative Polyadenylation in Populus alba var. pyramidalis by Single-Molecular Long-Read Sequencing

Poplars are worldwidely cultivated with ecologically and economically important value. Populus alba var. pyramidalis (= P. bolleana) is a main tree of the farmland shelter-belt system in the arid region of Northwest China due to its rapid growth, erect stems, and high biomass production. However, the full-length messenger RNA (mRNA) sequences and complete structure of P. alba var. pyramidalis remain unclear. In this study, using single-molecular real-time (SMRT) and next-generation high-throughput sequencing (NGS) platform, we sequenced transcripts from leaf, root, xylem, and phloem of P. alba var. pyramidalis, to obtain the full-length mRNA transcripts and annotate the complete structure. In total, 86,327 mapped full-length non-chimeric (FLNC) reads were identified, with 705 previously unannotated loci and 3,410 long noncoding RNAs (lncRNAs) and 174 fusion genes found. Alternative spicing (AS) events were detected in 7,536 genes, of which 4,652 genes had multiple AS events. A total of 10,213 alternative polyadenylation (APA) sites were identified, with two or more APA sites observed in 2,212 genes. Our transcriptome data provided the full-length sequences and gene isoforms of transcripts for P. alba var. pyramidalis, which will be helpful in improving our understanding for the genome annotation and gene structures of P. alba var. pyramidalis.

Localized Genetic and Phenotypic Diversity of Xanthomonas translucens Associated With Bacterial Leaf Streak on Wheat and Barley in Minnesota

Bacterial leaf streak (BLS) of wheat and barley has been a disease of increasing concern in the Upper Midwest over the past decade. In this study, intra- and interfield genetic and pathogenic diversity of bacteria causing BLS in Minnesota was evaluated. In 2015, 89 strains were isolated from 100 leaf samples collected from two wheat and two barley fields naturally infected with BLS. Virulence assays and multilocus sequence alignments of four housekeeping genes supported pathovar identifications. All wheat strains were pathogenic on wheat and barley and belonged to the same lineage as the Xanthomonas translucens pv. undulosa-type strain. All barley strains were pathogenic on barley but not on wheat. Three lineages of barley strains were detected. The frequency and number of sequence types of each pathovar varied within and between fields. A significant population variance was detected between populations of X. translucens pv. undulosa collected from different wheat fields. Population stratification of X. translucens pv. translucens was not detected. Significant differences in virulence were detected among three dominant sequence types of X. translucens pv. undulosa but not those of X. translucens pv. translucens. Field trials with wheat and barley plants inoculated with strains of known sequence type and virulence did not detect significant race structures within either pathovar. Knowledge of virulence, sequence types, and population structures of X. translucens on wheat and barley can support studies on plant-bacterial interactions and breeding for BLS disease resistance.

Genome Resource of Barley Bacterial Blight and Leaf Streak Pathogen Xanthomonas translucens pv. translucens strain UPB886

Xanthomonas translucens pv. translucens causes bacterial leaf streak and bacterial blight diseases of barley. This pathogen limits barley production globally but remains understudied, with limited genomic resources. To better understand the biology of this X. translucens subgroup, we sequenced the complete genome of the X. translucens pv. translucens strain UPB886.

TargeTALE: A Web Resource to Identify TALEs in Xanthomonas Genomes and Their Respective Targets

The Xanthomonas genus, comprises more than 30 species of gram-negative bacteria, most of which are pathogens of plants with high economic value, such as rice, common bean, and maize. Transcription activator-like effectors (TALEs), which act by regulating the host gene expression, are some of the major virulence factors of these bacteria. We present a novel tool to identify TALE genes in the genome of Xanthomonas strains and their respective targets. The analysis of the results obtained by TargeTALE in a proof-of-concept validation demonstrate that, at optimum setting, approximately 93% of the predicted target genes with available expression data were confirmed as upregulated during the infection, indicating that the tool might be useful for researchers in the field.

Single-molecule real-time sequencing facilitates the analysis of transcripts and splice isoforms of anthers in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

BACKGROUND

Anther development has been extensively studied at the transcriptional level, but a systematic analysis of full-length transcripts on a genome-wide scale has not yet been published. Here, the Pacific Biosciences (PacBio) Sequel platform and next-generation sequencing (NGS) technology were combined to generate full-length sequences and completed structures of transcripts in anthers of Chinese cabbage.

RESULTS

Using single-molecule real-time sequencing (SMRT), a total of 1,098,119 circular consensus sequences (CCSs) were generated with a mean length of 2664 bp. More than 75% of the CCSs were considered full-length non-chimeric (FLNC) reads. After error correction, 725,731 high-quality FLNC reads were estimated to carry 51,501 isoforms from 19,503 loci, consisting of 38,992 novel isoforms from known genes and 3691 novel isoforms from novel genes. Of the novel isoforms, we identified 407 long non-coding RNAs (lncRNAs) and 37,549 open reading frames (ORFs). Furthermore, a total of 453,270 alternative splicing (AS) events were identified and the majority of AS models in anther were determined to be approximate exon skipping (XSKIP) events. Of the key genes regulated during anther development, AS events were mainly identified in the genes SERK1, CALS5, NEF1, and CESA1/3. Additionally, we identified 104 fusion transcripts and 5806 genes that had alternative polyadenylation (APA).

CONCLUSIONS

Our work demonstrated the transcriptome diversity and complexity of anther development in Chinese cabbage. The findings provide a basis for further genome annotation and transcriptome research in Chinese cabbage.

Broad-spectrum resistance to bacterial blight in rice using genome editing

Bacterial blight of rice is an important disease in Asia and Africa. The pathogen, Xanthomonas oryzae pv. oryzae (Xoo), secretes one or more of six known transcription-activator-like effectors (TALes) that bind specific promoter sequences and induce, at minimum, one of the three host sucrose transporter genes SWEET11, SWEET13 and SWEET14, the expression of which is required for disease susceptibility. We used CRISPR-Cas9-mediated genome editing to introduce mutations in all three SWEET gene promoters. Editing was further informed by sequence analyses of TALe genes in 63 Xoo strains, which revealed multiple TALe variants for SWEET13 alleles. Mutations were also created in SWEET14, which is also targeted by two TALes from an African Xoo lineage. A total of five promoter mutations were simultaneously introduced into the rice line Kitaake and the elite mega varieties IR64 and Ciherang-Sub1. Paddy trials showed that genome-edited SWEET promoters endow rice lines with robust, broad-spectrum resistance.

Molecular Typing Reveals High Genetic Diversity of Xanthomonas translucens Strains Infecting Small-Grain Cereals in Iran

This study provides a phylogeographic insight into the population diversity of Xanthomonas translucens strains causing bacterial leaf streak disease of small-grain cereals in Iran. Among the 65 bacterial strains isolated from wheat, barley, and gramineous weeds in eight Iranian provinces, multilocus sequence analysis and typing (MLSA and MLST) of four housekeeping genes (dnaK, fyuA, gyrB, and rpoD), identified 57 strains as X. translucens pv. undulosa, while eight strains were identified as X. translucens pv. translucens. Although the pathogenicity patterns on oat and ryegrass weed species varied among the strains, all X. translucens pv. undulosa strains were pathogenic on barley, Harding's grass, rye (except for XtKm35) and wheat, and all X. translucens pv. translucens strains were pathogenic on barley and Harding's grass, while none of the latter group was pathogenic on rye or wheat (except for XtKm18). MLST using the 65 strains isolated in Iran, as well as the sequences of the four genes from 112 strains of worldwide origin retrieved from the GenBank database, revealed higher genetic diversity (i.e., haplotype frequency, haplotype diversity, and percentage of polymorphic sites) among the Iranian population of X. translucens than among the North American strains of the pathogen. High genetic diversity of the BLS pathogen in Iran was in congruence with the fact that the Iranian Plateau is considered the center of origin of cultivated wheat. However, further studies using larger collections of strains are warranted to precisely elucidate the global population diversity and center of origin of the pathogen.IMPORTANCE Bacterial leaf streak (BLS) of small-grain cereals (i.e., wheat and barley) is one of the economically important diseases of gramineous crops worldwide. The disease occurs in many countries across the globe, with particular importance in regions characterized by high levels of precipitation. Two genetically distinct xanthomonads-namely, Xanthomonas translucens pv. undulosa and X. translucens pv. translucens-have been reported to cause BLS disease on small-grain cereals. As seed-borne pathogens, the causal agents are included in the A2 list of quarantine pathogens by the European and Mediterranean Plant Protection Organization (EPPO). Despite its global distribution and high economic importance, the population structure, genetic diversity, and phylogeography of X. translucens remain undetermined. This study, using MLSA and MLST, provides a global-scale phylogeography of X. translucens strains infecting small-grain cereals. Based on the diversity parameters, neutrality indices, and population structure, we observe higher genetic diversity of the BLS pathogen in Iran, which is geographically close to the center of origin of common wheat, than has so far been observed in other areas of the world, including North America. The results obtained in this study provide a novel insight into the genetic diversity and population structure of the BLS pathogen of small-grain cereals on a global scale.

An efficient method to clone TAL effector genes from Xanthomonas oryzae using Gibson assembly

Transcription Activator-Like effectors (TALes) represent the largest family of type III effectors among pathogenic bacteria and play a critical role in the process of infection. Strains of Xanthomonas oryzae pv. oryzae (Xoo) and some strains of other Xanthomonas pathogens contain large numbers of TALe genes. Previous techniques to clone individual or a complement of TALe genes through conventional strategies are inefficient and time-consuming due to multiple genes (up to 29 copies) in a given genome, and technically challenging due to the repetitive sequences (up to 33 nearly identical 102-nucleotide repeats) of individual TALe genes. Thus, only a limited number of TALe genes have been molecularly cloned and characterized, and the functions of most TALe genes remain unknown. Here, we present an easy and efficient cloning technique to clone TALe genes selectively through in vitro homologous recombination and single-strand annealing, and demonstrate the feasibility of this approach with four different Xoo strains. Based on the Gibson assembly strategy, two complementary vectors with scaffolds that can preferentially capture all TALe genes from a pool of genomic fragments were designed. Both vector systems enabled cloning of a full complement of TALe genes from each of four Xoo strains and functional analysis of individual TALes in rice in approximately 1 month compared to 3 months by previously used methods. The results demonstrate a robust tool to advance TALe biology and a potential for broad usage of this approach to clone multiple copies of highly competitive DNA elements in any genome of interest.

Xanthomonas translucens commandeers the host rate-limiting step in ABA biosynthesis for disease susceptibility

Pathogenic bacteria acquire new virulence strategies for exploiting their hosts. This work reveals that the bacterial wheat pathogen Xanthomonas translucens uses a transcription activation-like (TAL) effector to promote virulence by directly activating the host gene 9-cis-epoxycarotenoid dioxygenase, the rate-limiting enzyme in biosynthesis of abscisic acid that is normally involved in water management within the host plant. Evolutionarily, TAL effectors are a relatively new class of virulence factors limited to a few species of pathogenic bacteria, and this work adds to the diversity of host susceptibility genes that can be exploited by pathogens through TAL effector gene function.

Plants are vulnerable to disease through pathogen manipulation of phytohormone levels, which otherwise regulate development, abiotic, and biotic responses. Here, we show that the wheat pathogen Xanthomonas translucens pv. undulosa elevates expression of the host gene encoding 9-cis-epoxycarotenoid dioxygenase (TaNCED-5BS), which catalyzes the rate-limiting step in the biosynthesis of the phytohormone abscisic acid and a component of a major abiotic stress-response pathway, to promote disease susceptibility. Gene induction is mediated by a type III transcription activator-like effector. The induction of TaNCED-5BS results in elevated abscisic acid levels, reduced host transpiration and water loss, enhanced spread of bacteria in infected leaves, and decreased expression of the central defense gene TaNPR1. The results represent an appropriation of host physiology by a bacterial virulence effector.

Tal1NXtc01 in Xanthomonas translucens pv. cerealis Contributes to Virulence in Bacterial Leaf Streak of Wheat

Xanthomonas translucens pv. cerealis (Xtc) causes bacterial leaf streak (BLS) of important cereal crops, including wheat (Triticum aestivum) and barley (Hordeum vulgare). Transcription activator-like effectors (TALEs) play vital roles in many plant diseases caused by Xanthomonas spp., however, TALEs have not been previously characterized in Xtc. In this study, the whole genome of NXtc01, a virulent strain of Xtc from Xinjiang, China, was sequenced and compared with genomes of other Xanthomonas spp. Xtc NXtc01 consists of a single 4,622,298 bp chromosome that encodes 4,004 genes. Alignment of the NXtc01 sequence with the draft genome of Xtc strain CFBP 2541 (United States) revealed a single giant inversion and differences in the location of two tal genes, which were designated tal1 and tal2. In NXtc01, both tal genes are located on the chromosome, whereas tal2 is plasmid-encoded in CFBP 2541. The repeat variable diresidues (RVDs) at the 12th and 13th sites within Tal2 repeat units were identical in both strains, whereas Tal1 showed differences in the third RVD. Xtc NXtc01 and CFBP 2541 encoded 35 and 33 non-TALE type III effectors (T3Es), respectively. tal1, tal2, and tal-free deletion mutants of Xtc NXtc01 were constructed and evaluated for virulence. The tal1 and tal-free deletion mutants were impaired with respect to symptom development and growth in wheat, suggesting that tal1 is a virulence factor in NXtc01. This was confirmed in gain-of-function experiments that showed the introduction of tal1, but not tal2, restored virulence to the tal-free mutant. Furthermore, we generated a hrcC deletion mutant of NXtc01; the hrcC mutant was non-pathogenic on wheat and unable to elicit a hypersensitive response in the non-host Nicotiana benthamiana. Our data provide a platform for exploring the roles of both TALEs and non-TALEs in promoting BLS on wheat.

Zinc binding proteome of a phytopathogen Xanthomonas translucens pv. undulosa

Xanthomonas translucens pv. undulosa (Xtu) is a proteobacteria which causes bacterial leaf streak (BLS) or bacterial chaff disease in wheat and barley. The constant competition for zinc (Zn) metal nutrients contributes significantly in plant-pathogen interactions. In this study, we have employed a systematic in silico approach to study the Zn-binding proteins of Xtu. From the whole proteome of Xtu, we have identified approximately 7.9% of proteins having Zn-binding sequence and structural motifs. Further, 115 proteins were found homologous to plant-pathogen interaction database. Among these 115 proteins, 11 were predicted as putative secretory proteins. The functional diversity in Zn-binding proteins was revealed by functional domain, gene ontology and subcellular localization analysis. The roles of Zn-binding proteins were found to be varied in the range from metabolism, proteolysis, protein biosynthesis, transport, cell signalling, protein folding, transcription regulation, DNA repair, response to oxidative stress, RNA processing, antimicrobial resistance, DNA replication and DNA integration. This study provides preliminary information on putative Zn-binding proteins of Xtu which may further help in designing new metal-based antimicrobial agents for controlling BLS and bacterial chaff infections on staple crops.

A Pathovar of Xanthomonas oryzae Infecting Wild Grasses Provides Insight Into the Evolution of Pathogenicity in Rice Agroecosystems

Xanthomonas oryzae (Xo) are globally important rice pathogens. Virulent lineages from Africa and Asia and less virulent strains from the United States have been well characterized. Xanthomonas campestris pv. leersiae (Xcl), first described in 1957, causes bacterial streak on the perennial grass, Leersia hexandra, and is a close relative of Xo. L. hexandra, a member of the Poaceae, is highly similar to rice phylogenetically, is globally ubiquitous around rice paddies, and is a reservoir of pathogenic Xo. We used long read, single molecule real time (SMRT) genome sequences of five strains of Xcl from Burkina Faso, China, Mali, and Uganda to determine the genetic relatedness of this organism with Xo. Novel transcription activator-like effectors (TALEs) were discovered in all five strains of Xcl. Predicted TALE target sequences were identified in the Leersia perrieri genome and compared to rice susceptibility gene homologs. Pathogenicity screening on L. hexandra and diverse rice cultivars confirmed that Xcl are able to colonize rice and produce weak but not progressive symptoms. Overall, based on average nucleotide identity (ANI), type III (T3) effector repertoires, and disease phenotype, we propose to rename Xcl to X. oryzae pv. leersiae (Xol) and use this parallel system to improve understanding of the evolution of bacterial pathogenicity in rice agroecosystems.

Analysis of transcripts and splice isoforms in Medicago sativa L. by single-molecule long-read sequencing

The full-length transcriptome of alfalfa was analyzed with PacBio single-molecule long-read sequencing technology. The transcriptome data provided full-length sequences and gene isoforms of transcripts in alfalfa, which will improve genome annotation and enhance our understanding of the gene structure of alfalfa. As an important forage, alfalfa (Medicago sativa L.) is world-wide planted. For its complexity of genome and unfinished whole genome sequencing, the sequences and complete structure of mRNA transcripts remain unclear in alfalfa. In this study, single-molecule long-read sequencing was applied to investigate the alfalfa transcriptome using the Pacific Biosciences platform, and a total of 113,321 transcripts were obtained from young, mature and senescent leaves. We identified 72,606 open reading frames including 46,616 full-length ORFs, 1670 transcription factors from 54 TF families and 44,040 simple sequence repeats from 30,797 sequences. A total of 7568 alternative splicing events was identified and the majority of alternative splicing events in alfalfa was intron retention. In addition, we identified 17,740 long non-coding RNAs. Our results show the feasibility of deep sequencing full-length RNA from alfalfa transcriptome on a single-molecule level.

Genetic Diversity and Virulence of Wheat and Barley Strains of Xanthomonas translucens from the Upper Midwestern United States

Bacterial leaf streak (BLS) of wheat and barley, caused by Xanthomonas translucens pv. undulosa and X. translucens pv. translucens, has been of growing concern in small grains production in the Upper Midwestern United States. To optimize disease resistance breeding, a greater awareness is needed of the pathovars and genetic diversity within the pathogens causing BLS in the region. Multilocus sequencing typing (MLST) and analysis (MLSA) of four common housekeeping genes (rpoD, dnaK, fyuA, and gyrB) was used to evaluate the genetic diversity of 82 strains of X. translucens isolated between 2006 and 2013 from wheat, barley, rye, and intermediate wheatgrass. In addition, in planta disease assays were conducted on 75 strains to measure relative virulence in wheat and barley. All strains were determined by MLSA to be related to X. translucens pv. undulosa and X. translucens pv. translucens. Clustering of strains based on Bayesian, network, and minimum spanning trees correlated with relative virulence levels in inoculated wheat and barley. Thus, phylogeny based on rpoD, dnaK, fyuA, and gyrB correlated with host of isolation and was an effective means for predicting virulence of strains belonging to X. translucens pv. translucens and X. translucens pv. undulosa.

Genetic mapping of a major gene in triticale conferring resistance to bacterial leaf streak

A major gene conferring resistance to bacterial leaf streak was mapped to chromosome 5R in triticale. Bacterial leaf streak (BLS), caused by Xanthomonas translucens pv. undulosa (Xtu), is an important disease of wheat and triticale around the world. Although resistance to BLS is limited in wheat, several triticale accessions have high levels of resistance. To characterize the genetic basis of this resistance, we developed triticale mapping populations using a resistant accession (Siskiyou) and two susceptible accessions (UC38 and Villax St. Jose). Bulked segregant analysis in an F₂ population derived from the cross of Siskiyou × UC38 led to the identification of a simple sequence repeat (SSR) marker (XSCM138) on chromosome 5R that co-segregated with the resistance gene. The cross of Siskiyou × Villax St. Jose was advanced into an F_2:5 recombinant inbred line population and evaluated for BLS reaction. Genetic linkage maps on this population were assembled with markers generated using genotyping-by-sequencing as well as several SSR markers previously identified on 5R. Quantitative trait locus (QTL) mapping revealed a single major QTL on chromosome 5R, underlined by the same SSR marker as in the Siskiyou × UC38 population. The F₁ hybrids of the two crosses were highly resistant to BLS, indicating that resistance is largely dominant. This work will facilitate introgression of this rye-derived BLS resistance gene into the wheat genome by molecular marker-mediated chromosome engineering.

Analysis of Extreme Phenotype Bulk Copy Number Variation (XP-CNV) Identified the Association of rp1 with Resistance to Goss's Wilt of Maize

Goss's wilt (GW) of maize is caused by the Gram-positive bacterium Clavibacter michiganensis subsp. nebraskensis (Cmn) and has spread in recent years throughout the Great Plains, posing a threat to production. The genetic basis of plant resistance is unknown. Here, a simple method for quantifying disease symptoms was developed and used to select cohorts of highly resistant and highly susceptible lines known as extreme phenotypes (XP). Copy number variation (CNV) analyses using whole genome sequences of bulked XP revealed 141 genes containing CNV between the two XP groups. The CNV genes include the previously identified common rust resistant locus rp1. Multiple Rp1 accessions with distinct rp1 haplotypes in an otherwise susceptible accession exhibited hypersensitive responses upon inoculation. GW provides an excellent system for the genetic dissection of diseases caused by closely related subspecies of C. michiganesis. Further work will facilitate breeding strategies to control GW and provide needed insight into the resistance mechanism of important related diseases such as bacterial canker of tomato and bacterial ring rot of potato.

Complete genome sequence of the sesame pathogen Ralstonia solanacearum strain SEPPX 05

Ralstonia solanacearum is a soil-borne phytopathogen associated with bacterial wilt disease of sesame. R. solanacearum is the predominant agent causing damping-off from tropical to temperate regions. Because bacterial wilt has decreased the sesame industry yield, we sequenced the SEPPX05 genome using PacBio and Illumina HiSeq 2500 systems and revealed that R. solanacearum strain SEPPX05 carries a bipartite genome consisting of a 3,930,849 bp chromosome and a 2,066,085 bp megaplasmid with 66.84% G+C content that harbors 5,427 coding sequences. Based on the whole genome, phylogenetic analysis showed that strain SEPPX05 is grouped with two phylotype I strains (EP1 and GMI1000). Pan-genomic analysis shows that R. solanacearum is a complex species with high biological diversity and was able to colonize various environments during evolution. Despite deletions, insertions, and inversions, most genes of strain SEPPX05 have relatively high levels of synteny compared with strain GMI1000. We identified 104 genes involved in virulence-related factors in the SEPPX05 genome and eight absent genes encoding T3Es of GMI1000. Comparing SEPPX05 with other species, we found highly conserved secretion systems central to modulating interactions of host bacteria. These data may provide important clues for understanding underlying pathogenic mechanisms of R. solanacearum and help in the control of sesame bacterial wilt.

Pacbio sequencing of copper-tolerant Xanthomonas citri reveals presence of a chimeric plasmid structure and provides insights into reassortment and shuffling of transcription activator-like effectors among X. citri strains

BACKGROUND

Xanthomonas citri, a causal agent of citrus canker, has been a well-studied model system due to recent availability of whole genome sequences of multiple strains from different geographical regions. Major limitations in our understanding of the evolution of pathogenicity factors in X. citri strains sequenced by short-read sequencing methods have been tracking plasmid reshuffling among strains due to inability to accurately assign reads to plasmids, and analyzing repeat regions among strains. X. citri harbors major pathogenicity determinants, including variable DNA-binding repeat region containing Transcription Activator-like Effectors (TALEs) on plasmids. The long-read sequencing method, PacBio, has allowed the ability to obtain complete and accurate sequences of TALEs in xanthomonads. We recently sequenced Xanthomonas citri str. Xc-03-1638-1-1, a copper tolerant A group strain isolated from grapefruit in 2003 from Argentina using PacBio RS II chemistry. We analyzed plasmid profiles, copy number and location of TALEs in complete genome sequences of X. citri strains.

RESULTS

We utilized the power of long reads obtained by PacBio sequencing to enable assembly of a complete genome sequence of strain Xc-03-1638-1-1, including sequences of two plasmids, 249 kb (plasmid harboring copper resistance genes) and 99 kb (pathogenicity plasmid containing TALEs). The pathogenicity plasmid in this strain is a hybrid plasmid containing four TALEs. Due to the intriguing nature of this pathogenicity plasmid with Tn3-like transposon association, repetitive elements and multiple putative sites for origins of replication, we might expect alternative structures of this plasmid in nature, illustrating the strong adaptive potential of X. citri strains. Analysis of the pathogenicity plasmid among completely sequenced X. citri strains, coupled with Southern hybridization of the pathogenicity plasmids, revealed clues to rearrangements of plasmids and resulting reshuffling of TALEs among strains.

CONCLUSIONS

We demonstrate in this study the importance of long-read sequencing for obtaining intact sequences of TALEs and plasmids, as well as for identifying rearrangement events including plasmid reshuffling. Rearrangement events, such as the hybrid plasmid in this case, could be a frequent phenomenon in the evolution of X. citri strains, although so far it is undetected due to the inability to obtain complete plasmid sequences with short-read sequencing methods.

Xanthomonas adaptation to common bean is associated with horizontal transfers of genes encoding TAL effectors

Background

Common bacterial blight is a devastating bacterial disease of common bean (Phaseolus vulgaris) caused by Xanthomonas citri pv. fuscans and Xanthomonas phaseoli pv. phaseoli. These phylogenetically distant strains are able to cause similar symptoms on common bean, suggesting that they have acquired common genetic determinants of adaptation to common bean. Transcription Activator-Like (TAL) effectors are bacterial type III effectors that are able to induce the expression of host genes to promote infection or resistance. Their capacity to bind to a specific host DNA sequence suggests that they are potential candidates for host adaption.

Results

To study the diversity of tal genes from Xanthomonas strains responsible for common bacterial blight of bean, whole genome sequences of 17 strains representing the diversity of X. citri pv. fuscans and X. phaseoli pv. phaseoli were obtained by single molecule real time sequencing. Analysis of these genomes revealed the existence of four tal genes named tal23A, tal20F, tal18G and tal18H, respectively. While tal20F and tal18G were chromosomic, tal23A and tal18H were carried on plasmids and shared between phylogenetically distant strains, therefore suggesting recent horizontal transfers of these genes between X. citri pv. fuscans and X. phaseoli pv. phaseoli strains. Strikingly, tal23A was present in all strains studied, suggesting that it played an important role in adaptation to common bean. In silico predictions of TAL effectors targets in the common bean genome suggested that TAL effectors shared by X. citri pv. fuscans and X. phaseoli pv. phaseoli strains target the promoters of genes of similar functions. This could be a trace of convergent evolution among TAL effectors from different phylogenetic groups, and comforts the hypothesis that TAL effectors have been implied in the adaptation to common bean.

Conclusions

Altogether, our results favour a model where plasmidic TAL effectors are able to contribute to host adaptation by being horizontally transferred between distant lineages.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4087-6) contains supplementary material, which is available to authorized users.

Complete Genome Sequencing and Targeted Mutagenesis Reveal Virulence Contributions of Tal2 and Tal4b of Xanthomonas translucens pv. undulosa ICMP11055 in Bacterial Leaf Streak of Wheat

Bacterial leaf streak caused by Xanthomonas translucens pv. undulosa (Xtu) is an important disease of wheat (Triticum aestivum) and barley (Hordeum vulgare) worldwide. Transcription activator-like effectors (TALEs) play determinative roles in many of the plant diseases caused by the different species and pathovars of Xanthomonas, but their role in this disease has not been characterized. ICMP11055 is a highly virulent Xtu strain from Iran. The aim of this study was to better understand genetic diversity of Xtu and to assess the role of TALEs in bacterial leaf streak of wheat by comparing the genome of this strain to the recently completely sequenced genome of a U.S. Xtu strain, and to several other draft X. translucens genomes, and by carrying out mutational analyses of the TALE (tal) genes the Iranian strain might harbor. The ICMP11055 genome, including its repeat-rich tal genes, was completely sequenced using single molecule, real-time technology (Pacific Biosciences). It consists of a single circular chromosome of 4,561,583 bp, containing 3,953 genes. Whole genome alignment with the genome of the United States Xtu strain XT4699 showed two major re-arrangements, nine genomic regions unique to ICMP11055, and one region unique to XT4699. ICMP110055 harbors 26 non-TALE type III effector genes and seven tal genes, compared to 25 and eight for XT4699. The tal genes occur singly or in pairs across five scattered loci. Four are identical to tal genes in XT4699. In addition to common repeat-variable diresidues (RVDs), the tal genes of ICMP11055, like those of XT4699, encode several RVDs rarely observed in Xanthomonas, including KG, NF, Y^∗, YD, and YK. Insertion and deletion mutagenesis of ICMP11055 tal genes followed by genetic complementation analysis in wheat cv. Chinese Spring revealed that Tal2 and Tal4b of ICMP11055 each contribute individually to the extent of disease caused by this strain. A largely conserved ortholog of tal2 is present in XT4699, but for tal4b, only a gene with partial, fragmented RVD sequence similarity can be found. Our results lay the foundation for identification of important host genes activated by Xtu TALEs as targets for the development of disease resistant varieties.

The genome of the cotton bacterial blight pathogen Xanthomonas citri pv. malvacearum strain MSCT1

Xanthomonas citri pv. malvacearum is a major pathogen of cotton, Gossypium hirsutum L.. In this study we report the complete genome of the X. citri pv. malvacearum strain MSCT1 assembled from long read DNA sequencing technology. The MSCT1 genome is the first X. citri pv. malvacearum genome with complete coding regions for X. citri pv. malvacearum transcriptional activator-like effectors. In addition functional and structural annotations are presented in this study that will provide a foundation for future pathogenesis studies with MSCT1.

Comparative Genomics Identifies a Novel Conserved Protein, HpaT, in Proteobacterial Type III Secretion Systems that Do Not Possess the Putative Translocon Protein HrpF

Xanthomonas translucens is the causal agent of bacterial leaf streak, the most common bacterial disease of wheat and barley. To cause disease, most xanthomonads depend on a highly conserved type III secretion system, which translocates type III effectors into host plant cells. Mutagenesis of the conserved type III secretion gene hrcT confirmed that the X. translucens type III secretion system is required to cause disease on the host plant barley and to trigger a non-host hypersensitive response (HR) in pepper leaves. Type III effectors are delivered to the host cell by a surface appendage, the Hrp pilus, and a translocon protein complex that inserts into the plant cell plasma membrane. Homologs of the Xanthomonas HrpF protein, including PopF from Ralstonia solanacearum and NolX from rhizobia, are thought to act as a translocon protein. Comparative genomics revealed that X. translucens strains harbor a noncanonical hrp gene cluster, which rather shares features with type III secretion systems from Ralstonia solanacearum, Paraburkholderia andropogonis, Collimonas fungivorans, and Uliginosibacterium gangwonense than other Xanthomonas spp. Surprisingly, none of these bacteria, except R. solanacearum, encode a homolog of the HrpF translocon. Here, we aimed at identifying a candidate translocon from X. translucens. Notably, genomes from strains that lacked hrpF/popF/nolX instead encode another gene, called hpaT, adjacent to and co-regulated with the type III secretion system gene cluster. An insertional mutant in the X. translucens hpaT gene, which is the first gene of a two-gene operon, hpaT-hpaH, was non-pathogenic on barley and did not cause the HR or programmed cell death in non-host pepper similar to the hrcT mutant. The hpaT mutant phenotypes were partially complemented by either hpaT or the downstream gene, hpaH, which has been described as a facilitator of translocation in Xanthomonas oryzae. Interestingly, the hpaT mutant was also complemented by the hrpF gene from Xanthomonas euvesicatoria. These findings reveal that both HpaT and HpaH contribute to the injection of type III effectors into plant cells.

Characterization of the Xanthomonas translucens Complex Using Draft Genomes, Comparative Genomics, Phylogenetic Analysis, and Diagnostic LAMP Assays

Prevalence of Xanthomonas translucens, which causes cereal leaf streak (CLS) in cereal crops and bacterial wilt in forage and turfgrass species, has increased in many regions in recent years. Because the pathogen is seedborne in economically important cereals, it is a concern for international and interstate germplasm exchange and, thus, reliable and robust protocols for its detection in seed are needed. However, historical confusion surrounding the taxonomy within the species has complicated the development of accurate and reliable diagnostic tools for X. translucens. Therefore, we sequenced genomes of 15 X. translucens strains representing six different pathovars and compared them with additional publicly available X. translucens genome sequences to obtain a genome-based phylogeny for robust classification of this species. Our results reveal three main clusters: one consisting of pv. cerealis, one consisting of pvs. undulosa and translucens, and a third consisting of pvs. arrhenatheri, graminis, phlei, and poae. Based on genomic differences, diagnostic loop-mediated isothermal amplification (LAMP) primers were developed that clearly distinguish strains that cause disease on cereals, such as pvs. undulosa, translucens, hordei, and secalis, from strains that cause disease on noncereal hosts, such as pvs. arrhenatheri, cerealis, graminis, phlei, and poae. Additional LAMP assays were developed that selectively amplify strains belonging to pvs. cerealis and poae, distinguishing them from other pathovars. These primers will be instrumental in diagnostics when implementing quarantine regulations to limit further geographic spread of X. translucens pathovars.

Comparative genomics of host adaptive traits in Xanthomonas translucens pv. graminis

BACKGROUND

Xanthomonas translucens pathovars differ in their individual host ranges among Poaceae. As the causal agent of bacterial wilt in Italian ryegrass (Lolium multiflorum Lam.), X. translucens pv. graminis (Xtg) is one of the most important bacterial pathogens in temperate grassland regions. The genomes of six Xtg strains from Switzerland, Norway, and New Zealand were sequenced in order to gain insight into conserved genomic traits from organisms covering a wide geographical range. Subsequent comparative analysis with previously published genome data of seven non-graminis X. translucens strains including the pathovars arrhenatheri, poae, phlei, cerealis, undulosa, and translucens was conducted to identify candidate genes linked to the host adaptation of Xtg to Italian ryegrass.

RESULTS

Phylogenetic analysis revealed a tight clustering of Xtg strains, which were found to share a large core genome. Conserved genomic traits included a non-canonical type III secretion system (T3SS) and a type IV pilus (T4P), which both revealed distinct primary structures of the pilins when compared to the non-graminis X. translucens strains. Xtg-specific traits that had no homologues in the other X. translucens strains were further found to comprise several hypothetical proteins, a TonB-dependent receptor, transporters, and effector proteins as well as toxin-antitoxin systems and DNA methyltransferases. While a nearly complete flagellar gene cluster was identified in one of the sequenced Xtg strains, phenotypic analysis pointed to swimming-deficiency as a common trait of the pathovar graminis.

CONCLUSION

Our study suggests that host adaptation of X. translucens pv. graminis may be conferred by a combination of pathovar-specific effector proteins, regulatory mechanisms, and adapted nutrient acquisition. Sequence deviations of pathogen-associated molecular patterns (PAMPs), as observed for the pilins of the T4P and T3SS, are moreover likely to impede perception by the plant defense machinery and thus facilitate successful host colonization of Italian ryegrass.