Nucleotide‐binding resistance gene signatures in sugar beet,insights from a new reference genome

Nucleotide‐binding (NB‐ARC), leucine‐rich‐repeat genes (NLRs) account for 60.8% of resistance (R) genes molecularly characterized from plants. NLRs exist as large gene families prone to tandem duplication and transposition, with high sequence diversity among crops and their wild relatives. This diversity can be a source of new disease resistance, but difficulty in distinguishing specific sequences from homologous gene family members hinders characterization of resistance for improving crop varieties. Current genome sequencing and assembly technologies, especially those using long‐read sequencing, are improving resolution of repeat‐rich genomic regions and clarifying locations of duplicated genes, such as NLRs. Using the conserved NB‐ARC domain as a model, 231 tentative NB‐ARC loci were identified in a highly contiguous genome assembly of sugar beet, revealing diverged and truncated NB‐ARC signatures as well as full‐length sequences. The NB‐ARC‐associated proteins contained NLR resistance gene domains, including TIR, CC and LRR, as well as other integrated domains. Phylogenetic relationships of partial and complete domains were determined, and patterns of physical clustering in the genome were evaluated. Comparison of sugar beet NB‐ARC domains to validated R‐genes from monocots and eudicots suggested extensive Beta vulgaris‐specific subfamily expansions. The NLR landscape in the rhizomania resistance conferring Rz region of Chromosome 3 was characterized, identifying 26 NLR‐like sequences spanning 20 MB. This work presents the first detailed view of NLR family composition in a member of the Caryophyllales, builds a foundation for additional disease resistance work in B. vulgaris, and demonstrates an additional nucleic‐acid‐based method for NLR prediction in non‐model plant species.     

The tomato I gene for Fusarium wilt resistance encodes an atypical leucine‐rich repeat receptor‐like protein whose function is nevertheless dependent on SOBIR1 and SERK3/BAK1

We have identified the tomato I gene for resistance to the Fusarium wilt fungus Fusarium oxysporum f. sp. lycopersici (Fol) and show that it encodes a membrane‐anchored leucine‐rich repeat receptor‐like protein (LRR‐RLP). Unlike most other LRR‐RLP genes involved in plant defence, the I gene is not a member of a gene cluster and contains introns in its coding sequence. The I gene encodes a loopout domain larger than those in most other LRR‐RLPs, with a distinct composition rich in serine and threonine residues. The I protein also lacks a basic cytosolic domain. Instead, this domain is rich in aromatic residues that could form a second transmembrane domain. The I protein recognises the Fol Avr1 effector protein, but, unlike many other LRR‐RLPs, recognition specificity is determined in the C‐terminal half of the protein by polymorphic amino acid residues in the LRRs just preceding the loopout domain and in the loopout domain itself. Despite these differences, we show that I/Avr1‐dependent necrosis in Nicotiana benthamiana depends on the LRR receptor‐like kinases (RLKs) SERK3/BAK1 and SOBIR1. Sequence comparisons revealed that the I protein and other LRR‐RLPs involved in plant defence all carry residues in their last LRR and C‐terminal LRR capping domain that are conserved with SERK3/BAK1‐interacting residues in the same relative positions in the LRR‐RLKs BRI1 and PSKR1. Tyrosine mutations of two of these conserved residues, Q922 and T925, abolished I/Avr1‐dependent necrosis in N. benthamiana, consistent with similar mutations in BRI1 and PSKR1 preventing their interaction with SERK3/BAK1.     

Extreme resistance to Potato virus Y in potato carrying the Rysto gene is mediated by a TIR‐NLR immune receptor

Potato virus Y (PVY) is a major potato (Solanum tuberosum L.) pathogen that causes severe annual crop losses worth billions of dollars worldwide. PVY is transmitted by aphids, and successful control of virus transmission requires the extensive use of environmentally damaging insecticides to reduce vector populations. Ry_sto, from the wild relative S. stoloniferum, confers extreme resistance (ER) to PVY and related viruses and is a valuable trait that is widely employed in potato resistance breeding programmes. Ry_sto was previously mapped to a region of potato chromosome XII, but the specific gene has not been identified to date. In this study, we isolated Ry_sto using resistance gene enrichment sequencing (RenSeq) and PacBio SMRT (Pacific Biosciences single‐molecule real‐time sequencing). Ry_sto was found to encode a nucleotide‐binding leucine‐rich repeat (NLR) protein with an N‐terminal TIR domain and was sufficient for PVY perception and ER in transgenic potato plants. Ry_sto‐dependent extreme resistance was temperature‐independent and requires EDS1 and NRG1 proteins. Ry_sto may prove valuable for creating PVY‐resistant cultivars of potato and other Solanaceae crops.     

The NB‐LRR proteins RGA4 and RGA5 interact functionally and physically to confer disease resistance

Plant resistance proteins of the class of nucleotide‐binding and leucine‐rich repeat domain proteins (NB‐LRRs) are immune sensors which recognize pathogen‐derived molecules termed avirulence (AVR) proteins. We show that RGA4 and RGA5, two NB‐LRRs from rice, interact functionally and physically to mediate resistance to the fungal pathogen Magnaporthe oryzae and accomplish different functions in AVR recognition. RGA4 triggers an AVR‐independent cell death that is repressed in the presence of RGA5 in both rice protoplasts and Nicotiana benthamiana. Upon recognition of the pathogen effector AVR‐Pia by direct binding to RGA5, repression is relieved and cell death occurs. RGA4 and RGA5 form homo‐ and hetero‐complexes and interact through their coiled‐coil domains. Localization studies in rice protoplast suggest that RGA4 and RGA5 localize to the cytosol. Upon recognition of AVR‐Pia, neither RGA4 nor RGA5 is re‐localized to the nucleus. These results establish a model for the interaction of hetero‐pairs of NB‐LRRs in plants: RGA4 mediates cell death activation, while RGA5 acts as a repressor of RGA4 and as an AVR receptor.     

Assembly and annotation of a draft genome sequence for Glycine latifolia,a perennial wild relative of soybean

Glycine latifolia (Benth.) Newell & Hymowitz (2n = 40), one of the 27 wild perennial relatives of soybean, possesses genetic diversity and agronomically favorable traits that are lacking in soybean. Here, we report the 939‐Mb draft genome assembly of G. latifolia (PI 559298) using exclusively linked‐reads sequenced from a single Chromium library. We organized scaffolds into 20 chromosome‐scale pseudomolecules utilizing two genetic maps and the Glycine max (L.) Merr. genome sequence. High copy numbers of putative 91‐bp centromere‐specific tandem repeats were observed in consecutive blocks within predicted pericentromeric regions on several pseudomolecules. No 92‐bp putative centromeric repeats, which are abundant in G. max, were detected in G. latifolia or Glycine tomentella. Annotation of the assembled genome and subsequent filtering yielded a high confidence gene set of 54 475 protein‐coding loci. In comparative analysis with five legume species, genes related to defense responses were significantly overrepresented in Glycine‐specific orthologous gene families. A total of 304 putative nucleotide‐binding site (NBS)‐leucine‐rich‐repeat (LRR) genes were identified in this genome assembly. Different from other legume species, we observed a scarcity of TIR‐NBS‐LRR genes in G. latifolia. The G. latifolia genome was also predicted to contain genes encoding 367 LRR‐receptor‐like kinases, a family of proteins involved in basal defense responses and responses to abiotic stress. The genome sequence and annotation of G. latifolia provides a valuable source of alternative alleles and novel genes to facilitate soybean improvement. This study also highlights the efficacy and cost‐effectiveness of the application of Chromium linked‐reads in diploid plant genome de novo assembly.     

Intragenic allele pyramiding combines different specificities of wheat Pm3 resistance alleles

Some plant resistance genes occur as allelic series, with each member conferring specific resistance against a subset of pathogen races. In wheat, there are 17 alleles of the Pm3 gene. They encode nucleotide‐binding (NB‐ARC) and leucine‐rich‐repeat (LRR) domain proteins, which mediate resistance to distinct race spectra of powdery mildew. It is not known if specificities from different alleles can be combined to create resistance genes with broader specificity. Here, we used an approach based on avirulence analysis of pathogen populations to characterize the molecular basis of Pm3 recognition spectra. A large survey of mildew races for avirulence on the Pm3 alleles revealed that Pm3a has a resistance spectrum that completely contains that of Pm3f, but also extends towards additional races. The same is true for the Pm3b and Pm3c gene pair. The molecular analysis of these allelic pairs revealed a role of the NB‐ARC protein domain in the efficiency of effector‐dependent resistance. Analysis of the wild‐type and chimeric Pm3 alleles identified single residues in the C‐terminal LRR motifs as the main determinant of allele specificity. Variable residues of the N‐terminal LRRs are necessary, but not sufficient, to confer resistance specificity. Based on these data, we constructed a chimeric Pm3 gene by intragenic allele pyramiding of Pm3d and Pm3e that showed the combined resistance specificity and, thus, a broader recognition spectrum compared with the parental alleles. Our findings support a model of stepwise evolution of Pm3 recognition specificities.     

Dissecting virulence function from recognition: cell death suppression in Nicotiana benthamiana by XopQ/HopQ1‐family effectors relies on EDS1‐dependent immunity

Many Gram‐negative plant pathogenic bacteria express effector proteins of the XopQ/HopQ1 family which are translocated into plant cells via the type III secretion system during infection. In Nicotiana benthamiana, recognition of XopQ/HopQ1 proteins induces an effector‐triggered immunity (ETI) reaction which is not associated with strong cell death but renders plants immune against Pseudomonas syringae and Xanthomonas campestris pv. vesicatoria strains. Additionally, XopQ suppresses cell death in N. benthamiana when transiently co‐expressed with cell death inducers. Here, we show that representative XopQ/HopQ1 proteins are recognized similarly, likely by a single resistance protein of the TIR‐NB‐LRR class. Extensive analysis of XopQ derivatives indicates the recognition of structural features. We performed Agrobacterium‐mediated protein expression experiments in wild‐type and EDS1‐deficient (eds1) N. benthamiana leaves, not recognizing XopQ/HopQ1. XopQ recognition limits multiplication of Agrobacterium and attenuates levels of transiently expressed proteins. Remarkably, XopQ fails to suppress cell death reactions induced by different effectors in eds1 plants. We conclude that XopQ‐mediated cell death suppression in N. benthamiana is due to the attenuation of Agrobacterium‐mediated protein expression rather than the cause of the genuine XopQ virulence activity. Thus, our study expands our understanding of XopQ recognition and function, and also challenges the commonly used co‐expression assays for elucidation of in planta effector activities, at least under conditions of ETI induction.     

The Brassica napus receptor‐like protein RLM2 is encoded by a second allele of the LepR3/Rlm2 blackleg resistance locus

Leucine‐rich repeat receptor‐like proteins (LRR‐RLPs) are highly adaptable parts of the signalling apparatus for extracellular detection of plant pathogens. Resistance to blackleg disease of Brassica spp. caused by Leptosphaeria maculans is largely governed by host race‐specific R‐genes, including the LRR‐RLP gene LepR3. The blackleg resistance gene Rlm2 was previously mapped to the same genetic interval as LepR3. In this study, the LepR3 locus of the Rlm2 Brassica napus line ‘Glacier DH24287’ was cloned, and B. napus transformants were analysed for recovery of the Rlm2 phenotype. Multiple B. napus, B. rapa and B. juncea lines were assessed for sequence variation at the locus. Rlm2 was found to be an allelic variant of the LepR3 LRR‐RLP locus, conveying race‐specific resistance to L. maculans isolates harbouring AvrLm2. Several defence‐related LRR‐RLPs have previously been shown to associate with the RLK SOBIR1 to facilitate defence signalling. Bimolecular fluorescence complementation (BiFC) and co‐immunoprecipitation of RLM2‐SOBIR1 studies revealed that RLM2 interacts with SOBIR1 of Arabidopsis thaliana when co‐expressed in Nicotiana benthamiana. The interaction of RLM2 with AtSOBIR1 is suggestive of a conserved defence signalling pathway between B. napus and its close relative A. thaliana.     

Engineering fire blight resistance into the apple cultivar ‘Gala’ using the FB_MR5 CC‐NBS‐LRR resistance gene of Malus × robusta 5

The fire blight susceptible apple cultivar Malus × domestica Borkh. cv. ‘Gala’ was transformed with the candidate fire blight resistance gene FB_MR5 originating from the crab apple accession Malus × robusta 5 (Mr5). A total of five different transgenic lines were obtained. All transgenic lines were shown to be stably transformed and originate from different transgenic events. The transgenic lines express the FB_MR5 either driven by the constitutive CaMV 35S promoter and the ocs terminator or by its native promoter and terminator sequences. Phenotyping experiments were performed with Mr5‐virulent and Mr5‐avirulent strains of Erwinia amylovora, the causal agent of fire blight. Significantly less disease symptoms were detected on transgenic lines after inoculation with two different Mr5‐avirulent E. amylovora strains, while significantly more shoot necrosis was observed after inoculation with the Mr5‐virulent mutant strain ZYRKD3_1. The results of these experiments demonstrated the ability of a single gene isolated from the native gene pool of apple to protect a susceptible cultivar from fire blight. Furthermore, this gene is confirmed to be the resistance determinant of Mr5 as the transformed lines undergo the same gene‐for‐gene interaction in the host–pathogen relationship Mr5–E. amylovora.     

A gene for resistance to the Varroa mite (Acari) in honey bee (Apis mellifera) pupae

Social insect colonies possess a range of defences which protect them against highly virulent parasites and colony collapse. The host–parasite interaction between honey bees (Apis mellifera) and the mite Varroa destructor is unusual, as honey bee colonies are relatively poorly defended against this parasite. The interaction has existed since the mid‐20th Century, when Varroa switched host to parasitize A. mellifera. The combination of a virulent parasite and relatively naïve host means that, without acaricides, honey bee colonies typically die within 3 years of Varroa infestation. A consequence of acaricide use has been a reduced selective pressure for the evolution of Varroa resistance in honey bee colonies. However, in the past 20 years, several natural‐selection‐based breeding programmes have resulted in the evolution of Varroa‐resistant populations. In these populations, the inhibition of Varroa's reproduction is a common trait. Using a high‐density genome‐wide association analysis in a Varroa‐resistant honey bee population, we identify an ecdysone‐induced gene significantly linked to resistance. Ecdysone both initiates metamorphosis in insects and reproduction in Varroa. Previously, using a less dense genetic map and a quantitative trait loci analysis, we have identified Ecdysone‐related genes at resistance loci in an independently evolved resistant population. Varroa cannot biosynthesize ecdysone but can acquire it from its diet. Using qPCR, we are able to link the expression of ecdysone‐linked resistance genes to Varroa's meals and reproduction. If Varroa co‐opts pupal compounds to initiate and time its own reproduction, mutations in the host's ecdysone pathway may represent a key selection tool for honey bee resistance and breeding.     

Genetic dissection of a TIR‐NB‐LRR locus from the wild North American grapevine species Muscadinia rotundifolia identifies paralogous genes conferring resistance to major fungal and oomycete pathogens in cultivated grapevine

The most economically important diseases of grapevine cultivation worldwide are caused by the fungal pathogen powdery mildew (Erysiphe necator syn. Uncinula necator) and the oomycete pathogen downy mildew (Plasmopara viticola). Currently, grapegrowers rely heavily on the use of agrochemicals to minimize the potentially devastating impact of these pathogens on grape yield and quality. The wild North American grapevine species Muscadinia rotundifolia was recognized as early as 1889 to be resistant to both powdery and downy mildew. We have now mapped resistance to these two mildew pathogens in M. rotundifolia to a single locus on chromosome 12 that contains a family of seven TIR‐NB‐LRR genes. We further demonstrate that two highly homologous (86% amino acid identity) members of this gene family confer strong resistance to these unrelated pathogens following genetic transformation into susceptible Vitis vinifera winegrape cultivars. These two genes, designated r esistance to P lasmopara v iticola (MrRPV1) are the first resistance genes to be cloned from a grapevine species. Both MrRUN1 and MrRPV1 were found to confer resistance to multiple powdery and downy mildew isolates from France, North America and Australia; however, a single powdery mildew isolate collected from the south‐eastern region of North America, to which M. rotundifolia is native, was capable of breaking MrRUN1‐mediated resistance. Comparisons of gene organization and coding sequences between M. rotundifolia and the cultivated grapevine V. vinifera at the MrRUN1/MrRPV1 locus revealed a high level of synteny, suggesting that the TIR‐NB‐LRR genes at this locus share a common ancestor.     

Stacking three late blight resistance genes from wild species directly into African highland potato varieties confers complete field resistance to local blight races

Considered responsible for one million deaths in Ireland and widespread famine in the European continent during the 1840s, late blight, caused by Phytophthora infestans, remains the most devastating disease of potato (Solanum tuberosum L.) with about 15%–30% annual yield loss in sub‐Saharan Africa, affecting mainly smallholder farmers. We show here that the transfer of three resistance (R) genes from wild relatives [RB, Rpi‐blb2 from Solanum bulbocastanum and Rpi‐vnt1.1 from S. venturii] into potato provided complete resistance in the field over several seasons. We observed that the stacking of the three R genes produced a high frequency of transgenic events with resistance to late blight. In the field, 13 resistant transgenic events with the 3R‐gene stack from the potato varieties ‘Desiree’ and ‘Victoria’ grew normally without showing pathogen damage and without any fungicide spray, whereas their non‐transgenic equivalent varieties were rapidly killed. Characteristics of the local pathogen population suggest that the resistance to late blight may be long‐lasting because it has low diversity, and essentially consists of the single lineage, 2_A1, which expresses the cognate avirulence effector genes. Yields of two transgenic events from ‘Desiree’ and ‘Victoria’ grown without fungicide to reflect small‐scale farm holders were estimated to be 29 and 45 t/ha respectively. This represents a three to four‐fold increase over the national average. Thus, these late blight resistant potato varieties, which are the farmers’ preferred varieties, could be rapidly adopted and bring significant income to smallholder farmers in sub‐Saharan Africa.     

The Lr34 adult plant rust resistance gene provides seedling resistance in durum wheat without senescence

The hexaploid wheat (Triticum aestivum) adult plant resistance gene, Lr34/Yr18/Sr57/Pm38/Ltn1, provides broad‐spectrum resistance to wheat leaf rust (Lr34), stripe rust (Yr18), stem rust (Sr57) and powdery mildew (Pm38) pathogens, and has remained effective in wheat crops for many decades. The partial resistance provided by this gene is only apparent in adult plants and not effective in field‐grown seedlings. Lr34 also causes leaf tip necrosis (Ltn1) in mature adult plant leaves when grown under field conditions. This D genome‐encoded bread wheat gene was transferred to tetraploid durum wheat (T. turgidum) cultivar Stewart by transformation. Transgenic durum lines were produced with elevated gene expression levels when compared with the endogenous hexaploid gene. Unlike nontransgenic hexaploid and durum control lines, these transgenic plants showed robust seedling resistance to pathogens causing wheat leaf rust, stripe rust and powdery mildew disease. The effectiveness of seedling resistance against each pathogen correlated with the level of transgene expression. No evidence of accelerated leaf necrosis or up‐regulation of senescence gene markers was apparent in these seedlings, suggesting senescence is not required for Lr34 resistance, although leaf tip necrosis occurred in mature plant flag leaves. Several abiotic stress‐response genes were up‐regulated in these seedlings in the absence of rust infection as previously observed in adult plant flag leaves of hexaploid wheat. Increasing day length significantly increased Lr34 seedling resistance. These data demonstrate that expression of a highly durable, broad‐spectrum adult plant resistance gene can be modified to provide seedling resistance in durum wheat.     

Comparison of spleen transcriptomes of two wild rodent species reveals differences in the immune response against Borrelia afzelii

Different host species often differ considerably in susceptibility to a given pathogen, but the causes of such differences are rarely known. The natural hosts of the tick‐transmitted bacterium Borrelia afzelii, which is one of causative agents of Lyme borreliosis in humans, include a variety of small mammals like voles and mice. Previous studies have shown that B. afzelii‐infected bank voles (Myodes glareolus) have about ten times higher bacterial load than infected yellow‐necked mice (Apodemus flavicollis), indicating that these two species differ in resistance. In this study, we compared the immune response to B. afzelii infection in these host species by using RNA sequencing to quantify gene expression in spleen. Gene set enrichment analysis (GSEA) showed that several immune pathways were down‐regulated in infected animals in both bank voles and yellow‐necked mice. Moreover, IFNα response was up‐regulated in B. afzelii‐infected yellow‐necked mice, while IL6 signaling and the complement pathway were down‐regulated in infected bank voles; differences in regulation of these three pathways between bank voles and yellow‐necked mice could thus contribute to the difference in resistance to B. afzelii between the species. This study provides knowledge of gene expression induced by a zoonotic pathogen in its natural host, and possible species‐specific regulation of immune responses associated with resistance.     

Non‐GMO potato lines,synthesizing increased amylose and resistant starch,are mainly deficient in isoamylase debranching enzyme

Solanum tuberosum potato lines with high amylose content were generated by crossing with the wild potato species Solanum sandemanii followed by repeated backcrossing to Solanum tuberosum lines. The trait, termed increased amylose (IAm), was recessive and present after three generations of backcrossing into S. tuberosum lines (6.25% S. sandemanii genes). The tubers of these lines were small, elongated and irregular with small and misshaped starch granules and high sugar content. Additional backcrossing resulted in less irregular tuber morphology, increased starch content (4.3%–9.5%) and increased amylose content (29%–37.9%) but indifferent sugar content. The amylose in the IAm starch granules was mainly located in peripheral spots, and large cavities were found in the granules. Starch pasting was suppressed, and the digestion‐resistant starch (RS) content was increased. Comprehensive microarray polymer profiling (CoMPP) analysis revealed specific alterations of major pectic and glycoprotein cell wall components. This complex phenotype led us to search for candidate IAm genes exploiting its recessive trait. Hence, we sequenced genomic DNA of a pool of IAm lines, identified SNPs genome wide against the draft genome sequence of potato and searched for regions of decreased heterozygosity. Three regions, located on chromosomes 3, 7 and 10, respectively, displayed markedly less heterozygosity than average. The only credible starch metabolism‐related gene found in these regions encoded the isoamylase‐type debranching enzyme Stisa1. Decreased expression of mRNA (>500 fold) and reduced enzyme activity (virtually absent from IAm lines) supported Stisa1 as a candidate gene for IAm.     

Genome‐wide association analysis and gene set enrichment analysis with SNP data identify genes associated with 305‐day milk yield in Holstein dairy cows

Milk production traits, such as 305‐day milk yield (305MY), have been under direct selection to improve production in dairy cows. Over the past 50 years, the average milk yield has nearly doubled, and over 56% of the increase is attributable to genetic improvement. As such, additional improvements in milk yield are still possible as new loci are identified. The objectives of this study were to detect SNPs and gene sets associated with 305MY in order to identify new candidate genes contributing to variation in milk production. A population of 781 primiparous Holstein cows from six central Washington dairies with records of 305MY and energy corrected milk were used to perform a genome‐wide association analysis (GWAA) using the Illumina BovineHD BeadChip (777 962 SNPs) to identify QTL associated with 305MY (P < 1.0 × 10^?5). A gene set enrichment analysis with SNP data (GSEA‐SNP) was performed to identify gene sets (normalized enrichment score > 3.0) and leading edge genes (LEGs) influencing 305MY. The GWAA identified three QTL comprising 34 SNPs and 30 positional candidate genes. In the GSEA‐SNP, five gene sets with 58 unique and 24 shared LEGs contributed to 305MY. Identification of QTL and LEGs associated with 305MY can provide additional targets for genomic selection to continue to improve 305MY in dairy cattle.