Genome-wide identification of the NLR gene family in Haynaldia villosa by SMRT-RenSeq

Background

Nucleotide-binding and leucine-rich repeat (NLR) genes have attracted wide attention due to their crucial role in protecting plants from pathogens. SMRT-RenSeq, combining PacBio sequencing after resistance gene enrichment sequencing (RenSeq), is a powerful method for selectively capturing and sequencing full-length NLRs. Haynaldia villosa, a wild grass species with a proven potential for wheat improvement, confers resistance to multiple diseases. So, genome-wide identification of the NLR gene family in Haynaldia villosa by SMRT-RenSeq can facilitate disease resistance genes exploration.

Results

In this study, SMRT-RenSeq was performed to identify the genome-wide NLR complement of H. villosa. In total, 1320 NLRs were annotated in 1169 contigs, including 772 complete NLRs. All the complete NLRs were phylogenetically analyzed and 11 main clades with special characteristics were derived. NLRs could be captured with high efficiency when aligned with cloned R genes, and cluster expansion in some specific gene loci was observed. The physical location of NLRs to individual chromosomes in H. villosa showed a perfect homoeologous relationship with group 1, 2, 3, 5 and 6 of other Triticeae species, however, NLRs physically located on 4VL were largely in silico predicted to be located on the homoeologous group 7. Fifteen types of integrated domains (IDs) were integrated in 52 NLRs, and Kelch and B3 NLR-IDs were found to have expanded in H. villosa, while DUF948, NAM-associated and PRT_C were detected as unique integrated domains implying the new emergence of NLR-IDs after H. villosa diverged from other species.

Conclusion

SMRT-RenSeq is a powerful tool to identify NLR genes from wild species using the baits of the evolutionary related species with reference sequences. The availability of the NLRs from H. villosa provide a valuable library for R gene mining and transfer of disease resistance into wheat.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08334-w.

The Exocyst Complex Subunit EXO70E1-V From Haynaldia villosa Interacts With Wheat Powdery Mildew Resistance Gene CMPG1-V

EXO70 belongs to the exocyst complex subunit that plays a critical role in regulating plant cell polarity establishment and defense response. A previous study proved that the E3 ligase CMPG1-V from Haynaldia villosa, a diploid wheat relative, positively regulates the resistance to wheat powdery mildew (Pm), caused by fungus Blumeria graminis f.sp tritici (Bgt). In this study, a member of EXO70 superfamily named EXO70E1-V was isolated from H. villosa, and EXO70E1-V interacted with CMPG1-V were shown by yeast two-hybrid (Y2H), pull-down assay, bimolecular fluorescence complementation (BiFC) assay, and luciferase complementation imaging (LCI) assay. It is localized in various subcellular organs, i.e., plasma membrane (PM) and endoplasmic reticulum. Co-expression of EXO70E1-V and CMPG1-V showed dot-like structure fluorescence signals that were mainly in PM and nucleus. Expression of EXO70E1-V was relatively higher in leaf and was significantly induced by Bgt infection and exogenous application of hormones such as salicylic acid. Transient or stable overexpression of EXO70E1-V could not enhance/decrease the Pm resistance level, suggesting overexpression of EXO70E1-V alone has no impact on Pm resistance in wheat.

The Regulatory Network of CMPG1-V in Wheat-Blumeria graminis f. sp. tritici Interaction Revealed by Temporal Profiling Using RNA-Seq

Wheat powdery mildew (Pm), caused by Blumeria graminis f. sp. tritici (Bgt), is a prevalent fungal disease. The diploid wheat relative Haynaldia villosa (H. villosa) showed broad-spectrum resistance (BSR) to Pm. A previous study reported an E3 ligase gene, CMPG1-V from H. villosa, showing BSR to Pm. To elucidate the regulatory network mediated by CMPG1-V, in this study, gene expression profiling of CMPG1-V transgenic plant (CMPG1-V_OE) and its receptor Yangmai 158 was analyzed and compared after Bgt inoculation at four infection stages. GO and KEGG analysis revealed obvious reprogramming of SA and ABA signaling, starch/sucrose metabolism, and photosynthesis in CMPG1-V_OE, compared with those in Yangmai 158. Transcripts of SA synthesis genes SARD1 and UGT, signaling factors TGA and PRs, and SnRKs in ABA signaling were specifically upregulated in CMPG1-V_OE rather than Yangmai 158. Transcripts of LHCII in photosynthesis, GLUC and TPP in starch/sucrose metabolism were also induced distinctly in CMPG1-V_OE. WGCNA analysis showed crucial regulatory candidates of CMPG1-V, involving serine/threonine-protein kinase in phosphorylation, glucosyltransferase in flavonoid biosynthesis, defense factor WRKYs, and peroxidase in oxidative stress. Our results facilitate the deciphering of the resistant regulatory network of CMPG1-V and the identification of key candidates which might be employed in breeding programs.

Identification and Characterization of the EXO70 Gene Family in Polyploid Wheat and Related Species

The EXO70 gene family is involved in different biological processes in plants, ranging from plant polar growth to plant immunity. To date, analysis of the EXO70 gene family has been limited in Triticeae species, e.g., hexaploidy Triticum aestivum and its ancestral/related species. By in silico analysis of multiple Triticeae sequence databases, a total of 200 EXO70 members were identified. By homologue cloning approaches, 15 full-length cDNA of EXO70s were cloned from diploid Haynaldia villosa. Phylogenetic relationship analysis of 215 EXO70 members classified them into three groups (EXO70.1, EXO70.2, and EXO70.3) and nine subgroups (EXO70A to EXO70I). The distribution of most EXO70 genes among different species/sub-genomes were collinear, implying their orthologous relationship. The EXO70A subgroup has the most introns (at least five introns), while the remaining seven subgroups have only one intron on average. The expression profiling of EXO70 genes from wheat revealed that 40 wheat EXO70 genes were expressed in at least one tissue (leaf, stem, or root), of which 25 wheat EXO70 genes were in response to at least one biotic stress (stripe rust or powdery mildew) or abiotic stress (drought or heat). Subcellular localization analysis showed that ten EXO70-V proteins had distinct plasma membrane localization, EXO70I1-V showed a distinctive spotted pattern on the membrane. The 15 EXO70-V genes were differentially expressed in three tissue. Apart from EXO70D2-V, the remaining EXO70-V genes were in response to at least one stress (flg22, chitin, powdery mildew, drought, NaCl, heat, or cold) or phytohormones (salicylic acid, methyl jasmonate, ethephon, or abscisic acid) and hydrogen peroxide treatments. This research provides a genome-wide glimpse of the Triticeae EXO70 gene family and those up- or downregulated genes require further validation of their biological roles in response to biotic/abiotic stresses.

Sequencing flow-sorted short arm of Haynaldia villosa chromosome 4V provides insights into its molecular structure and virtual gene order

BACKGROUND

Haynaldia villosa (H. villosa) has been recognized as a species potentially useful for wheat improvement. The availability of its genomic sequences will boost its research and application.

RESULTS

In this work, the short arm of H. villosa chromosome 4V (4VS) was sorted by flow cytometry and sequenced using Illumina platform. About 170.6 Mb assembled sequences were obtained. Further analysis showed that repetitive elements accounted for about 64.6% of 4VS, while the coding fraction, which is corresponding to 1977 annotated genes, represented 1.5% of the arm. The syntenic regions of the 4VS were searched and identified on wheat group 4 chromosomes 4AL, 4BS, 4DS, Brachypodium chromosomes 1 and 4, rice chromosomes 3 and 11, and sorghum chromosomes 1, 5 and 8. Based on genome-zipper analysis, a virtual gene order comprising 735 gene loci on 4VS genome was built by referring to the Brachypodium genome, which was relatively consistent with the scaffold order determined for Ae. tauschii chromosome 4D. The homologous alleles of several cloned genes on wheat group 4 chromosomes including Rht-1 gene were identified.

CONCLUSIONS

The sequences provided valuable information for mapping and positional-cloning genes located on 4VS, such as the wheat yellow mosaic virus resistance gene Wss1. The work on 4VS provided detailed insights into the genome of H. villosa, and may also serve as a model for sequencing the remaining parts of H. villosa genome.

Development of intron targeting (IT) markers specific for chromosome arm 4VS of Haynaldia villosa by chromosome sorting and next-generation sequencing

Background

Haynaldia villosa (L.) Schur (syn. Dasypyrum villosum L. Candargy, 2n = 14, genome VV) is the tertiary gene pool of wheat, and thus a potential resource of genes for wheat improvement. Among other, wheat yellow mosaic (WYM) resistance gene Wss1 and a take-all resistance gene were identified on the short arm of chromosome 4 V (4VS) of H. villosa. We had obtained introgressions on 4VS chromosome arm, with the objective of utilizing the target genes. However, monitoring these introgressions has been a daunting task because of inadequate knowledge as to H.villosa genome, as reflected by the lack of specific markers.

Results

This study aims to develop 4VS-specific markers by combination of chromosome sorting and next-generation sequencing. The short arm of chromosome 4VS of H.villosa was flow-sorted using a FACSVantage SE flow cytometer and sorter, and then sequenced by Illumina sequencing. The sequence of H. villosa 4VS was assembled by the software Hecate, and then was compared with the sequence assemblies of wheat chromosome arms 4AL, 4BS and 4DS and Ae. tauschii 4DS, with the objectives of identifying exon-exon junctions and localizing introns on chromosome 4VS of H. villosa. The intron length polymorphisms suitable for designing H. villosa primers were evaluated with criteria. Consequently, we designed a total of 359 intron targeting (IT) markers, among which 232 (64.62%) markers were specific for tracing the 4VS chromatin in the wheat background.

Conclusion

The combination of chromosome sorting and next-generation sequencing to develop specific IT markers for 4VS of H. villosa has high success rate and specificity, thus being applicable for the development of chromosome-specific markers for alien chromatin in wheat breeding.

Electronic supplementary material

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

Haynaldia villosa NAM-V1 is linked with the powdery mildew resistance gene Pm21 and contributes to increasing grain protein content in wheat

BACKGROUND

The 6AL/6VS translocation lines, carrying the wheat powdery mildew resistance gene Pm21, are planted on more than 3.4 million hectares. The NAM-A1 gene, located on chromosome 6AS of hexaploid wheat, has been implicated with increased wheat grain protein content (GPC). However, the NAM-A1 gene was removed from the 6AL/6VS translocation lines after the original chromosome 6AS was replaced by chromosome 6VS of Haynaldia villosa. The present study aimed to clone the NAM homologous gene from chromosome 6VS, to analyze the changes of GPC in the 6AL/6VS translocation lines, and to develop related molecular markers for wheat molecular breeding.

RESULTS

A new NAM family gene, NAM-V1, was cloned from 6VS of H. villosa (GenBank ACC. no. KR873101). NAM-V1 contained an intact open reading frame (ORF) and putatively encodes a protein of 407 amino acids. Phylogenetic analysis indicated that NAM-V1 was an orthologous gene of NAM-A1, B1, and D1. The determination of GPC in four Pm21 F2 segregation populations demonstrated that the replacement of NAM-A1 by NAM-V1 confers increased GPC in hexaploid wheat. Multiple sequence alignment of NAM-A1, B1, B2, D1, D2, and V1 showed the single nucleotide polymorphism (SNP) sites for each of the NAM genes, allowing us to develop a molecular marker, CauNAM-V1, for the specific detection of NAM-V1 gene. Our results indicate that CauNAM-V1 can be used as a novel DNA marker for NAM-V1, and can also be used for selecting Pm21 in wheat breeding programs. Further, we developed a marker, CauNAM-ABD, for the amplification and simultaneously distinguish among the NAM-A1, NAM-B1, NAM-B2, NAM-D1, and NAM-D2 genes in a single step. CauNAM-ABD enabled us to develop an efficient "one-marker-for-five-genes" procedure for identifying genes and its copy numbers related with grain protein content.

CONCLUSION

Here, we report the isolation of the NAM-V1 gene of H. villosa. This gene contributes to increasing GPC in 6AL/6VS translocation wheat lines. We developed a molecular marker for the specific detection of NAM-V1 and a molecular marker that can be used to simultaneously distinguished among the NAM-A1, NAM-B1, NAM-B2, NAM-D1, and NAM-D2 genes in a single step.