CRISPR-Cas9 Gene Editing of the Sal1 Gene Family in Wheat.

The highly conserved Sal1 encodes a bifunctional enzyme with inositol polyphosphate-1-phosphatase and 3' (2'), 5'-bisphosphate nucleotidase activity and has been shown to alter abiotic stress tolerance in plants when disrupted. Precise gene editing techniques were used to generate Sal1 mutants in hexaploid bread wheat. The CRISPR (Clustered Regulatory Interspaced Short Palindromic Repeats) Cas9 system with three guide RNAs (gRNAs) was used to inactivate six Sal1 homologous genes within the Bobwhite wheat genome. The resulting mutant wheat plants with all their Sal1 genes disabled had slimmer stems, had a modest reduction in biomass and senesced more slowly in water limiting conditions, but did not exhibit improved yield under drought conditions. Our results show that multiplexed gRNAs enabled effective targeted gene editing of the Sal1 gene family in hexaploid wheat. These Sal1 mutant wheat plants will be a resource for further research studying the function of this gene family in wheat.

Exploiting the reference genome sequence of hexaploid wheat: a proteomic study of flour proteins from the cultivar Chinese Spring.

Although the economic value of wheat flour is determined by the complement of gluten proteins, these proteins have been challenging to study because of the complexity of the major protein groups and the tremendous sequence diversity among wheat cultivars. The completion of a high-quality wheat genome sequence from the reference wheat Chinese Spring recently facilitated the assembly and annotation of a complete set of gluten protein genes from a single cultivar, making it possible to link individual proteins in the flour to specific gene sequences. In a proteomic analysis of total wheat flour protein from Chinese Spring using quantitative two-dimensional gel electrophoresis combined with tandem mass spectrometry, gliadins or low-molecular-weight glutenin subunits were identified as the predominant proteins in 72 protein spots. Individual spots were associated with 40 of 56 Chinese Spring gene sequences, including 16 of 26 alpha gliadins, 10 of 11 gamma gliadins, six of seven omega gliadins, one of two delta gliadins, and nine of ten LMW-GS. Most genes that were not associated with protein spots were either expressed at low levels in endosperm or encoded proteins with high similarity to other proteins. A wide range of protein accumulation levels were observed and discrepancies between transcript levels and protein levels were noted. This work together with similar studies using other commercial cultivars should provide new insight into the molecular basis of wheat flour quality and allergenic potential.

Rapid evolution of α-gliadin gene family revealed by analyzing Gli-2 locus regions of wild emmer wheat.

α-Gliadins are a major group of gluten proteins in wheat flour that contribute to the end-use properties for food processing and contain major immunogenic epitopes that can cause serious health-related issues including celiac disease (CD). α-Gliadins are also the youngest group of gluten proteins and are encoded by a large gene family. The majority of the gene family members evolved independently in the A, B, and D genomes of different wheat species after their separation from a common ancestral species. To gain insights into the origin and evolution of these complex genes, the genomic regions of the Gli-2 loci encoding α-gliadins were characterized from the tetraploid wild emmer, a progenitor of hexaploid bread wheat that contributed the AABB genomes. Genomic sequences of Gli-2 locus regions for the wild emmer A and B genomes were first reconstructed using the genome sequence scaffolds along with optical genome maps. A total of 24 and 16 α-gliadin genes were identified for the A and B genome regions, respectively. α-Gliadin pseudogene frequencies of 86% for the A genome and 69% for the B genome were primarily caused by C to T substitutions in the highly abundant glutamine codons, resulting in the generation of premature stop codons. Comparison with the homologous regions from the hexaploid wheat cv. Chinese Spring indicated considerable sequence divergence of the two A genomes at the genomic level. In comparison, conserved regions between the two B genomes were identified that included α-gliadin pseudogenes containing shared nested TE insertions. Analyses of the genomic organization and phylogenetic tree reconstruction indicate that although orthologous gene pairs derived from speciation were present, large portions of α-gliadin genes were likely derived from differential gene duplications or deletions after the separation of the homologous wheat genomes ~ 0.5 MYA. The higher number of full-length intact α-gliadin genes in hexaploid wheat than that in wild emmer suggests that human selection through domestication might have an impact on α-gliadin evolution. Our study provides insights into the rapid and dynamic evolution of genomic regions harboring the α-gliadin genes in wheat.

Gene Duplication and Evolution Dynamics in the Homeologous Regions Harboring Multiple Prolamin and Resistance Gene Families in Hexaploid Wheat.

Improving end-use quality and disease resistance are important goals in wheat breeding. The genetic loci controlling these traits are highly complex, consisting of large families of prolamin and resistance genes with members present in all three homeologous A, B, and D genomes in hexaploid bread wheat. Here, orthologous regions harboring both prolamin and resistance gene loci were reconstructed and compared to understand gene duplication and evolution in different wheat genomes. Comparison of the two orthologous D regions from the hexaploid wheat Chinese Spring and the diploid progenitor Aegilops tauschii revealed their considerable difference due to the presence of five large structural variations with sizes ranging from 100 kb to 2 Mb. As a result, 44% of the Ae. tauschii and 71% of the Chinese Spring sequences in the analyzed regions, including 79 genes, are not shared. Gene rearrangement events, including differential gene duplication and deletion in the A, B, and D regions, have resulted in considerable erosion of gene collinearity in the analyzed regions, suggesting rapid evolution of prolamin and resistance gene families after the separation of the three wheat genomes. We hypothesize that this fast evolution is attributed to the co-evolution of the two gene families dispersed within a high recombination region. The identification of a full set of prolamin genes facilitated transcriptome profiling and revealed that the A genome contributes the least to prolamin expression because of its smaller number of expressed intact genes and their low expression levels, while the B and D genomes contribute similarly.

Dynamic Evolution of α-Gliadin Prolamin Gene Family in Homeologous Genomes of Hexaploid Wheat.

Wheat Gli-2 loci encode complex groups of α-gliadin prolamins that are important for breadmaking, but also major triggers of celiac disease (CD). Elucidation of α-gliadin evolution provides knowledge to produce wheat with better end-use properties and reduced immunogenic potential. The Gli-2 loci contain a large number of tandemly duplicated genes and highly repetitive DNA, making sequence assembly of their genomic regions challenging. Here, we constructed high-quality sequences spanning the three wheat homeologous α-gliadin loci by aligning PacBio-based sequence contigs with BioNano genome maps. A total of 47 α-gliadin genes were identified with only 26 encoding intact full-length protein products. Analyses of α-gliadin loci and phylogenetic tree reconstruction indicate significant duplications of α-gliadin genes in the last ~2.5 million years after the divergence of the A, B and D genomes, supporting its rapid lineage-independent expansion in different Triticeae genomes. We showed that dramatic divergence in expression of α-gliadin genes could not be attributed to sequence variations in the promoter regions. The study also provided insights into the evolution of CD epitopes and identified a single indel event in the hexaploid wheat D genome that likely resulted in the generation of the highly toxic 33-mer CD epitope.

New insights into structural organization and gene duplication in a 1.75-Mb genomic region harboring the α-gliadin gene family in Aegilops tauschii, the source of wheat D genome.

Among the wheat prolamins important for its end-use traits, α-gliadins are the most abundant, and are also a major cause of food-related allergies and intolerances. Previous studies of various wheat species estimated that between 25 and 150 α-gliadin genes reside in the Gli-2 locus regions. To better understand the evolution of this complex gene family, the DNA sequence of a 1.75-Mb genomic region spanning the Gli-2 locus was analyzed in the diploid grass, Aegilops tauschii, the ancestral source of D genome in hexaploid bread wheat. Comparison with orthologous regions from rice, sorghum, and Brachypodium revealed rapid and dynamic changes only occurring to the Ae. tauschii Gli-2 region, including insertions of high numbers of non-syntenic genes and a high rate of tandem gene duplications, the latter of which have given rise to 12 copies of α-gliadin genes clustered within a 550-kb region. Among them, five copies have undergone pseudogenization by various mutation events. Insights into the evolutionary relationship of the duplicated α-gliadin genes were obtained from their genomic organization, transcription patterns, transposable element insertions and phylogenetic analyses. An ancestral glutamate-like receptor (GLR) gene encoding putative amino acid sensor in all four grass species has duplicated only in Ae. tauschii and generated three more copies that are interspersed with the α-gliadin genes. Phylogenetic inference and different gene expression patterns support functional divergence of the Ae. tauschii GLR copies after duplication. Our results suggest that the duplicates of α-gliadin and GLR genes have likely taken different evolutionary paths; conservation for the former and neofunctionalization for the latter.

Deep sequencing and transcriptome analyses to identify genes involved in secoiridoid biosynthesis in the Tibetan medicinal plant Swertia mussotii.

Swertia mussotii Franch. is an important traditional Tibetan medicinal plant with pharmacological properties effective in the treatment of various ailments including hepatitis. Secoiridoids are the major bioactive compounds in S. mussotii. To better understand the secoiridoid biosynthesis pathway, we generated transcriptome sequences from the root, leaf, stem, and flower tissues, and performed de novo sequence assembly, yielding 98,613 unique transcripts with an N50 of 1,085 bp. Putative functions could be assigned to 35,029 transcripts (35.52%) based on BLAST searches against annotation databases including GO and KEGG. The expression profiles of 39 candidate transcripts encoding the key enzymes for secoiridoid biosynthesis were examined in different S. mussotii tissues, validated by qRT-PCR, and compared with the homologous genes from S. japonica, a species in the same family, unveiling the gene expression, regulation, and conservation of the pathway. The examination of the accumulated levels of three bioactive compounds, sweroside, swertiamarin, and gentiopicroside, revealed their considerable variations in different tissues, with no significant correlation with the expression profiles of key genes in the pathway, suggesting complex biological behaviours in the coordination of metabolite biosynthesis and accumulation. The genomic dataset and analyses presented here lay the foundation for further research on this important medicinal plant.

Specialization-generalization trade-off in a Bradyrhizobium symbiosis with wild legume hosts.

BACKGROUND: Specialized interactions help structure communities, but persistence of specialized organisms is puzzling because a generalist can occupy more environments and partake in more beneficial interactions. The "Jack-of-all-trades is a master of none" hypothesis asserts that specialists persist because the fitness of a generalist utilizing a particular habitat is lower than that of a specialist adapted to that habitat. Yet, there are many reasons to expect that mutualists will generalize on partners.Plant-soil feedbacks help to structure plant and microbial communities, but how frequently are soil-based symbiotic mutualistic interactions sufficiently specialized to influence species distributions and community composition? To address this question, we quantified realized partner richness and phylogenetic breadth of four wild-grown native legumes (Lupinus bicolor, L. arboreus, Acmispon strigosus and A. heermannii) and performed inoculation trials to test the ability of two hosts (L. bicolor and A. strigosus) to nodulate (fundamental partner richness), benefit from (response specificity), and provide benefit to (effect specificity) 31 Bradyrhizobium genotypes. RESULTS: In the wild, each Lupinus species hosted a broader genetic range of Bradyrhizobium than did either Acmispon species, suggesting that Acmispon species are more specialized. In the greenhouse, however, L. bicolor and A. strigosus did not differ in fundamental association specificity: all inoculated genotypes nodulated both hosts. Nevertheless, A. strigosus exhibited more specificity, i.e., greater variation in its response to, and effect on, Bradyrhizobium genotypes. Lupinus bicolor benefited from a broader range of genotypes but averaged less benefit from each. Both hosts obtained more fitness benefit from symbionts isolated from conspecific hosts; those symbionts in turn gained greater fitness benefit from hosts of the same species from which they were isolated. CONCLUSIONS: This study affirmed two important tenets of evolutionary theory. First, as predicted by the Jack-of-all-trades is a master of none hypothesis, specialist A. strigosus obtained greater benefit from its beneficial symbionts than did generalist L. bicolor. Second, as predicted by coevolutionary theory, each test species performed better with partner genotypes isolated from conspecifics. Finally, positive fitness feedback between the tested hosts and symbionts suggests that positive plant-soil feedback could contribute to their patchy distributions in this system.

The Arabidopsis downy mildew resistance gene RPP8 is induced by pathogens and salicylic acid and is regulated by W box cis elements.

Plants disease resistance (R) genes encode specialized receptors that are quantitative, rate-limiting defense regulators. R genes must be expressed at optimum levels to function properly. If expression is too low, downstream defense responses are not activated efficiently. Conversely, overexpression of R genes can trigger autoactivation of defenses with deleterious consequences for the plant. Little is known about R gene regulation, particularly under defense-inducing conditions. We examined regulation of the Arabidopsis thaliana gene RPP8 (resistance to Hyaloperonospora arabidopsidis, isolate Emco5). RPP8 was induced in response to challenge with H. arabidopsidis or application of salicylic acid, as shown with RPP8-Luciferase transgenic plants and quantitative reverse-transcription polymerase chain reaction of endogenous alleles. The RPP1 and RPP4 genes were also induced by H. arabidopsidis and salicylic acid, suggesting that some RPP genes are subject to feedback amplification. The RPP8 promoter contains three W box cis elements. Site-directed mutagenesis of all three W boxes greatly diminished RPP8 basal expression, inducibility, and resistance in transgenic plants. Motif searches indicated that the W box is the only known cis element that is statistically overrepresented in Arabidopsis nucleotide-binding leucine-rich repeat promoters. These results indicate that WRKY transcription factors can regulate expression of surveillance genes at the top of the defense-signaling cascade.

Role of recombination in the evolution of the model plant pathogen Pseudomonas syringae pv. tomato DC3000, a very atypical tomato strain.

Pseudomonas syringae pv. tomato strain DC3000 (PtoDC3000) is one of the most intensively studied bacterial plant pathogens today. Here we report a thorough investigation into PtoDC3000 and close relatives isolated from Antirrhinum majus (snapdragon), Apium graveolens (celery), and Solanaceae and Brassicaceae species. Multilocus sequence typing (MLST) was used to resolve the precise phylogenetic relationship between isolates and to determine the importance of recombination in their evolution. MLST data were correlated with an analysis of the locus coding for the type III secreted (T3S) effector AvrPto1 to investigate the role of recombination in the evolution of effector repertoires. Host range tests were performed to determine if closely related isolates from different plants have different host ranges. It was found that PtoDC3000 is located in the same phylogenetic cluster as isolates from several Brassicaceae and Solanaceae species and that these isolates have a relatively wide host range that includes tomato, Arabidopsis thaliana, and cauliflower. All other analyzed tomato isolates from three different continents form a distinct cluster and are pathogenic only on tomato. Therefore, PtoDC3000 is a very unusual tomato isolate. Several recombination breakpoints were detected within sequenced gene fragments, and population genetic tests indicate that recombination contributed more than mutation to the variation between isolates. Moreover, recombination may play an important role in the reassortment of T3S effectors between strains. The data are finally discussed from a taxonomic standpoint, and P. syringae pv. tomato is proposed to be divided into two pathovars.

Naturally occurring nonpathogenic isolates of the plant pathogen Pseudomonas syringae lack a type III secretion system and effector gene orthologues.

Pseudomonas syringae causes plant diseases, and the main virulence mechanism is a type III secretion system (T3SS) that translocates dozens of effector proteins into plant cells. Here we report the existence of a subgroup of P. syringae isolates that do not cause disease on any plant species tested. This group is monophyletic and most likely evolved from a pathogenic P. syringae ancestor through loss of the T3SS. In the nonpathogenic isolate P. syringae 508 the genomic region that in pathogenic P. syringae strains contains the hrp-hrc cluster coding for the T3SS and flanking effector genes is absent. P. syringae 508 was also surveyed for the presence of effector orthologues from the closely related pathogenic strain P. syringae pv. syringae B728a, but none were detected. The absence of the hrp-hrc cluster and effector orthologues was confirmed for other nonpathogenic isolates. Using the AvrRpt2 effector as reporter revealed the inability of P. syringae 508 to translocate effectors into plant cells. Adding a plasmid-encoded T3SS and the P. syringae pv. syringae 61 effector gene hopA1 increased in planta growth almost 10-fold. This suggests that P. syringae 508 supplemented with a T3SS could be used to determine functions of individual effectors in the context of a plant infection, avoiding the confounding effect of other effectors with similar functions present in effector mutants of pathogenic isolates.