Full-genome identification and characterization of NBS-encoding disease resistance genes in wheat

Selection of Candidate Genes Conferring Blast Resistance and Heat Tolerance in Rice through Integration of Meta-QTLs and RNA-Seq

Due to global warming, high temperature is a significant environmental stress for rice production. Rice (Oryza sativa L.), one of the most crucial cereal crops, is also seriously devastated by Magnaporthe oryzae. Therefore, it is essential to breed new rice cultivars with blast and heat tolerance. Although progress had been made in QTL mapping and RNA-seq analysis in rice in response to blast and heat stresses, there are few reports on simultaneously mining blast-resistant and heat-tolerant genes. In this study, we separately conducted meta-analysis of 839 blast-resistant and 308 heat-tolerant QTLs in rice. Consequently, 7054 genes were identified in 67 blast-resistant meta-QTLs with an average interval of 1.00 Mb. Likewise, 6425 genes were obtained in 40 heat-tolerant meta-QTLs with an average interval of 1.49 Mb. Additionally, using differentially expressed genes (DEGs) in the previous research and GO enrichment analysis, 55 DEGs were co-located on the common regions of 16 blast-resistant and 14 heat-tolerant meta-QTLs. Among, OsChib3H-c, OsJAMyb, Pi-k, OsWAK1, OsMT2b, OsTPS3, OsHI-LOX, OsACLA-2 and OsGS2 were the significant candidate genes to be further investigated. These results could provide the gene resources for rice breeding with excellent resistance to these 2 stresses, and help to understand how plants response to the combination stresses of blast fungus and high temperature.

NBS-LRR Gene TaRPS2 is Positively Associated with the High-Temperature Seedling Plant Resistance of Wheat Against Puccinia striiformis f. sp. tritici

Xiaoyan6 (XY6) is a wheat (Triticum aestivum) cultivar possessing nonrace-specific high-temperature seedling plant (HTSP) resistance against stripe rust, caused by Puccinia striiformis f. sp. tritici. Previously, we identified one particular gene, TaRPS2, for its involvement in the HTSP resistance. To elucidate the role of TaRPS2 in the HTSP resistance, we cloned the full length of TaRPS2 from XY6. The transcriptional expression of TaRPS2 was rapidly upregulated (19.11-fold) under the normal-high-normal temperature treatment that induces the HTSP resistance. The expression level of TaRPS2 in leaves was higher than that in the stems and roots. Quantification of the endogenous hormones in wheat leaves after P. striiformis f. sp. tritici inoculation showed that 1-aminocyclopropane-1-carboxylic acid, salicylic acid (SA), and jasmonic acid were involved in the HTSP resistance. In addition, detection of hydrogen peroxide (H₂O₂) accumulation indicated that reactive oxygen species burst was also associated with the HTSP resistance. Two hours after exogenous H₂O₂ treatment or 0.5 h after SA treatment, the expression level of TaRPS2 was increased by 2.66 and 2.35 times, respectively. The subcellular localization of enhanced green fluorescent protein-TaRPS2 fusion protein was in the nuclei and plasma membranes. Virus-induced gene silencing of TaRPS2 reduced the level of HTSP resistance in XY6. Compared with the nonsilenced leaves, the TaRPS2-silenced leaves had the reduction of necrotic cells but a greater number of uredinia. These results indicated that TaRPS2 positively regulates the HTSP resistance of XY6 against P. striiformis f. sp. tritici and is related to the SA and H₂O₂ signaling pathways.

De Novo Genome Sequence Assembly of Dwarf Coconut (Cocos nucifera L. 'Catigan Green Dwarf') Provides Insights into Genomic Variation Between Coconut Types and Related Palm Species

We report the first whole genome sequence (WGS) assembly and annotation of a dwarf coconut variety, 'Catigan Green Dwarf' (CATD). The genome sequence was generated using the PacBio SMRT sequencing platform at 15X coverage of the expected genome size of 2.15 Gbp, which was corrected with assembled 50X Illumina paired-end MiSeq reads of the same genome. The draft genome was improved through Chicago sequencing to generate a scaffold assembly that results in a total genome size of 2.1 Gbp consisting of 7,998 scaffolds with N50 of 570,487 bp. The final assembly covers around 97.6% of the estimated genome size of coconut 'CATD' based on homozygous k-mer peak analysis. A total of 34,958 high-confidence gene models were predicted and functionally associated to various economically important traits, such as pest/disease resistance, drought tolerance, coconut oil biosynthesis, and putative transcription factors. The assembled genome was used to infer the evolutionary relationship within the palm family based on genomic variations and synteny of coding gene sequences. Data show that at least three (3) rounds of whole genome duplication occurred and are commonly shared by these members of the Arecaceae family. A total of 7,139 unique SSR markers were designed to be used as a resource in marker-based breeding. In addition, we discovered 58,503 variants in coconut by aligning the Hainan Tall (HAT) WGS reads to the non-repetitive regions of the assembled CATD genome. The gene markers and genome-wide SSR markers established here will facilitate the development of varieties with resilience to climate change, resistance to pests and diseases, and improved oil yield and quality.

Structural, Functional, and Genomic Diversity of Plant NLR Proteins: An Evolved Resource for Rational Engineering of Plant Immunity

Plants employ a diverse intracellular system of NLR (nucleotide binding-leucine-rich repeat) innate immune receptors to detect pathogens of all types. These receptors represent valuable agronomic traits that plant breeders rely on to maximize yield in the face of devastating pathogens. Despite their importance, the mechanistic underpinnings of NLR-based disease resistance remain obscure. The rapidly increasing numbers of plant genomes are revealing a diverse array of NLR-type immune receptors. In parallel, mechanistic studies are describing diverse functions for NLR immune receptors. In this review, we intend to broadly describe how the structural, functional, and genomic diversity of plant immune receptors can provide a valuable resource for rational engineering of plant immunity.

Genome-Wide Association Study Identifies NBS-LRR-Encoding Genes Related with Anthracnose and Common Bacterial Blight in the Common Bean

Nucleotide-binding site and leucine-rich repeat (NBS-LRR) genes represent the largest and most important disease resistance genes in plants. The genome sequence of the common bean (Phaseolus vulgaris L.) provides valuable data for determining the genomic organization of NBS-LRR genes. However, data on the NBS-LRR genes in the common bean are limited. In total, 178 NBS-LRR-type genes and 145 partial genes (with or without a NBS) located on 11 common bean chromosomes were identified from genome sequences database. Furthermore, 30 NBS-LRR genes were classified into Toll/interleukin-1 receptor (TIR)-NBS-LRR (TNL) types, and 148 NBS-LRR genes were classified into coiled-coil (CC)-NBS-LRR (CNL) types. Moreover, the phylogenetic tree supported the division of these PvNBS genes into two obvious groups, TNL types and CNL types. We also built expression profiles of NBS genes in response to anthracnose and common bacterial blight using qRT-PCR. Finally, we detected nine disease resistance loci for anthracnose (ANT) and seven for common bacterial blight (CBB) using the developed NBS-SSR markers. Among these loci, NSSR24, NSSR73, and NSSR265 may be located at new regions for ANT resistance, while NSSR65 and NSSR260 may be located at new regions for CBB resistance. Furthermore, we validated NSSR24, NSSR65, NSSR73, NSSR260, and NSSR265 using a new natural population. Our results provide useful information regarding the function of the NBS-LRR proteins and will accelerate the functional genomics and evolutionary studies of NBS-LRR genes in food legumes. NBS-SSR markers represent a wide-reaching resource for molecular breeding in the common bean and other food legumes. Collectively, our results should be of broad interest to bean scientists and breeders.

Epistatic determinism of durum wheat resistance to the wheat spindle streak mosaic virus

KEY MESSAGE

The resistance of durum wheat to the Wheat spindle streak mosaic virus (WSSMV) is controlled by two main QTLs on chromosomes 7A and 7B, with a huge epistatic effect. Wheat spindle streak mosaic virus (WSSMV) is a major disease of durum wheat in Europe and North America. Breeding WSSMV-resistant cultivars is currently the only way to control the virus since no treatment is available. This paper reports studies of the inheritance of WSSMV resistance using two related durum wheat populations obtained by crossing two elite cultivars with a WSSMV-resistant emmer cultivar. In 2012 and 2015, 354 recombinant inbred lines (RIL) were phenotyped using visual notations, ELISA and qPCR and genotyped using locus targeted capture and sequencing. This allowed us to build a consensus genetic map of 8568 markers and identify three chromosomal regions involved in WSSMV resistance. Two major regions (located on chromosomes 7A and 7B) jointly explain, on the basis of epistatic interactions, up to 43% of the phenotypic variation. Flanking sequences of our genetic markers are provided to facilitate future marker-assisted selection of WSSMV-resistant cultivars.

Genome-wide identification and characterization of NB-ARC resistant genes in wheat (Triticum aestivum L.) and their expression during leaf rust infection

NB-ARC domain-containing resistance genes from the wheat genome were identified, characterized and localized on chromosome arms that displayed differential yet positive response during incompatible and compatible leaf rust interactions. Wheat (Triticum aestivum L.) is an important cereal crop; however, its production is affected severely by numerous diseases including rusts. An efficient, cost-effective and ecologically viable approach to control pathogens is through host resistance. In wheat, high numbers of resistance loci are present but only few have been identified and cloned. A comprehensive analysis of the NB-ARC-containing genes in complete wheat genome was accomplished in this study. Complete NB-ARC encoding genes were mined from the Ensembl Plants database to predict 604 NB-ARC containing sequences using the HMM approach. Genome-wide analysis of orthologous clusters in the NB-ARC-containing sequences of wheat and other members of the Poaceae family revealed maximum homology with Oryza sativa indica and Brachypodium distachyon. The identification of overlap between orthologous clusters enabled the elucidation of the function and evolution of resistance proteins. The distributions of the NB-ARC domain-containing sequences were found to be balanced among the three wheat sub-genomes. Wheat chromosome arms 4AL and 7BL had the most NB-ARC domain-containing contigs. The spatio-temporal expression profiling studies exemplified the positive role of these genes in resistant and susceptible wheat plants during incompatible and compatible interaction in response to the leaf rust pathogen Puccinia triticina. Two NB-ARC domain-containing sequences were modelled in silico, cloned and sequenced to analyze their fine structures. The data obtained in this study will augment isolation, characterization and application NB-ARC resistance genes in marker-assisted selection based breeding programs for improving rust resistance in wheat.

Rapid evolutionary dynamics in a 2.8-Mb chromosomal region containing multiple prolamin and resistance gene families in Aegilops tauschii

Prolamin and resistance gene families are important in wheat food use and in defense against pathogen attacks, respectively. To better understand the evolution of these multi-gene families, the DNA sequence of a 2.8-Mb genomic region, representing an 8.8 cM genetic interval and harboring multiple prolamin and resistance-like gene families, was analyzed in the diploid grass Aegilops tauschii, the D-genome donor of bread wheat. Comparison with orthologous regions from rice, Brachypodium, and sorghum showed that the Ae. tauschii region has undergone dramatic changes; it has acquired more than 80 non-syntenic genes and only 13 ancestral genes are shared among these grass species. These non-syntenic genes, including prolamin and resistance-like genes, originated from various genomic regions and likely moved to their present locations via sequence evolution processes involving gene duplication and translocation. Local duplication of non-syntenic genes contributed significantly to the expansion of gene families. Our analysis indicates that the insertion of prolamin-related genes occurred prior to the separation of the Brachypodieae and Triticeae lineages. Unlike in Brachypodium, inserted prolamin genes have rapidly evolved and expanded to encode different classes of major seed storage proteins in Triticeae species. Phylogenetic analyses also showed that the multiple insertions of resistance-like genes and subsequent differential expansion of each R gene family. The high frequency of non-syntenic genes and rapid local gene evolution correlate with the high recombination rate in the 2.8-Mb region with nine-fold higher than the genome-wide average. Our results demonstrate complex evolutionary dynamics in this agronomically important region of Triticeae species.

Large-scale bioinformatic analysis of the regulation of the disease resistance NBS gene family by microRNAs in Poaceae

In the present study, we have screened 71, 713, 525, 119 and 241 mature miRNA variants from Hordeum vulgare, Oryza sativa, Brachypodium distachyon, Triticum aestivum, and Sorghum bicolor, respectively, and classified them with respect to their conservation status and expression levels. These Poaceae non-redundant miRNA species (1,669) were distributed over a total of 625 MIR families, among which only 54 were conserved across two or more plant species, confirming the relatively recent evolutionary differentiation of miRNAs in grasses. On the other hand, we have used 257 H. vulgare, 286T. aestivum, 119 B. distachyon, 269 O. sativa, and 139 S. bicolor NBS domains, which were either mined directly from the annotated proteomes, or predicted from whole genome sequence assemblies. The hybridization potential between miRNAs and their putative NBS genes targets was analyzed, revealing that at least 454 NBS genes from all five Poaceae were potentially regulated by 265 distinct miRNA species, most of them expressed in leaves and predominantly co-expressed in additional tissues. Based on gene ontology, we could assign these probable miRNA target genes to 16 functional groups, among which three conferring resistance to bacteria (Rpm1, Xa1 and Rps2), and 13 groups of resistance to fungi (Rpp8,13, Rp3, Tsn1, Lr10, Rps1-k-1, Pm3, Rpg5, and MLA1,6,10,12,13). The results of the present analysis provide a large-scale platform for a better understanding of biological control strategies of disease resistance genes in Poaceae, and will serve as an important starting point for enhancing crop disease resistance improvement by means of transgenic lines with artificial miRNAs.

Comparative Analysis and Identification of miRNAs and Their Target Genes Responsive to Salt Stress in Diploid and Tetraploid Paulownia fortunei Seedlings

Salt stress is a global environmental problem that affects plant growth and development. Paulownia fortunei is an adaptable and fast-growing deciduous tree native to China that is environmentally and economically important. MicroRNAs (miRNAs) play important regulatory roles in growth, development, and stress responses in plants. MiRNAs that respond to biotic stresses have been identified; however, how miRNAs in P. fortunei respond to salt stress has not yet been reported. To identify salt-stress-responsive miRNAs and predict their target genes, four small RNA and four degradome libraries were constructed from NaCl-treated and NaCl-free leaves of P. fortunei seedlings. The results indicated that salt stress had different physiological effects on diploid and tetraploid P. fortunei. We detected 53 conserved miRNAs belonging to 17 miRNA families and 134 novel miRNAs in P. fortunei. Comparing their expression levels in diploid and tetraploid P. fortunei, we found 10 conserved and 10 novel miRNAs that were significantly differentially expressed under salt treatment, among them eight were identified as miRNAs probably associated with higher salt tolerance in tetraploid P. fortunei than in diploid P. fortunei. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses were performed to predict the functions of the target genes of the conserved and novel miRNAs. The expressions of 10 differentially expressed miRNAs were validated by quantitative real-time polymerase chain reaction (qRT-PCR). This is the first report on P. fortunei miRNAs and their target genes under salt stress. The results provided information at the physiological and molecular levels for further research into the response mechanisms of P. fortunei to salt stress.

A genome-wide analysis of the auxin/indole-3-acetic acid gene family in hexaploid bread wheat (Triticum aestivum L.)

The Auxin/indole-3-acetic acid (Aux/IAA) gene family plays key roles in the primary auxin-response process and controls a number of important traits in plants. However, the characteristics of the Aux/IAA gene family in hexaploid bread wheat (Triticum aestivum L.) have long been unknown. In this study, a comprehensive identification of the Aux/IAA gene family was performed using the latest draft genome sequence of the bread wheat "Chinese Spring." Thirty-four Aux/IAA genes were identified, 30 of which have duplicated genes on the A, B or D sub-genome, with a total of 84 Aux/IAA sequences. These predicted Aux/IAA genes were non-randomly distributed in all the wheat chromosomes except for chromosome 2D. The information of wheat Aux/IAA proteins is also described. Based on an analysis of phylogeny, expression and adaptive evolution, we prove that the Aux/IAA family in wheat has been replicated twice in the two allopolyploidization events of bread wheat, when the tandem duplication also occurred. The duplicated genes have undergone an evolutionary process of purifying selection, resulting in the high conservation of copy genes among sub-genomes and functional redundancy among several members of the TaIAA family. However, functional divergence probably existed in most TaIAA members due to the diversity of the functional domain and expression pattern. Our research provides useful information for further research into the function of Aux/IAA genes in wheat.

Strategies for transferring resistance into wheat: from wide crosses to GM cassettes

The domestication of wheat in the Fertile Crescent 10,000 years ago led to a genetic bottleneck. Modern agriculture has further narrowed the genetic base by introducing extreme levels of uniformity on a vast spatial and temporal scale. This reduction in genetic complexity renders the crop vulnerable to new and emerging pests and pathogens. The wild relatives of wheat represent an important source of genetic variation for disease resistance. For nearly a century farmers, breeders, and cytogeneticists have sought to access this variation for crop improvement. Several barriers restricting interspecies hybridization and introgression have been overcome, providing the opportunity to tap an extensive reservoir of genetic diversity. Resistance has been introgressed into wheat from at least 52 species from 13 genera, demonstrating the remarkable plasticity of the wheat genome and the importance of such natural variation in wheat breeding. Two main problems hinder the effective deployment of introgressed resistance genes for crop improvement: (1) the simultaneous introduction of genetically linked deleterious traits and (2) the rapid breakdown of resistance when deployed individually. In this review, we discuss how recent advances in molecular genomics are providing new opportunities to overcome these problems.