Jianping Zhang
Peng Zhang
Peter Dodds
Evans Lagudah- 1CSIRO Agriculture & Food, Canberra, ACT, Australia
- 2Plant Breeding Institute Cobbitty, The University of Sydney, Sydney, NSW, Australia
The wheat-rust pathosystem has been well-studied among host–pathogen interactions since last century due to its economic importance. Intensified efforts toward cloning of wheat rust resistance genes commenced in the late 1990s with the first successful isolation published in 2003. Currently, a total of 24 genes have been cloned from wheat that provides resistance to stem rust, leaf rust, and stripe rust. Among them, more than half (15) were cloned over the last 4 years. This rapid cloning of resistance genes from wheat can be largely credited to the development of approaches for reducing the genome complexity as 10 out of the 15 genes cloned recently were achieved by approaches that are summarized as TEnSeq (Target-sequence Enrichment and Sequencing) pipelines in this review. The growing repertoire of cloned rust resistance genes provides new tools to support deployment strategies aimed at achieving durable resistance. This will be supported by the identification of genetic variation in corresponding Avr genes from rust pathogens, which has recently begun. Although developed with wheat resistance genes as the primary targets, TEnSeq approaches are also applicable to other classes of genes as well as for other crops with complex genomes.
Introduction
Wheat crops are afflicted by three major rust diseases, namely stem/black rust, stripe/yellow rust, and leaf rust/brown rust, each caused by a different fungal species in the genus Puccinia. Significant yield losses from each of the three diseases have been reported from almost all major wheat-growing regions worldwide. Losses to leaf rust were approximately $350 million between 2000 and 2004 in the United States (Huerta-Espino et al., 2011). In China, leaf rust causes yield losses estimated at 3 million tons annually (Huerta-Espino et al., 2011). Leaf rust was also reported as a severe threat to wheat crops in Mexico and South Asia in the past, but with the utilization of slow rusting resistance genes in some areas like the United States, damage in recent decades has been reduced substantially (Huerta-Espino et al., 2011). Stripe rust epidemics were previously restricted mainly to cooler and humid regions, such as those in Asia and Europe. However, the appearance and spread of aggressive races that have adapted to warmer climates have expanded the geographic footprint of this disease since 2000, resulting in severe losses in many countries (Wellings et al., 2012). Annual global losses due to stripe rust were recently estimated at USD $979 million (Beddow et al., 2015). Stripe rust is also considered the most damaging wheat rust disease in Australia, with annual economic losses valued at around AUD$127 million (Murray and Brennan, 2009). Stem rust is the most destructive wheat attacking rusts and has historically been especially damaging in Africa, the Americas, Europe, and Australia. Numerous severe stem rust epidemics occurred in the United States during the first half of the 20th century, causing average yield losses of 19.3 to 28.4% in some states (Roelfs, 1978). The wide utilization of resistant cultivars adopted during the Green Revolution and the eradication of the alternate host barberry since 1954 in the United States resulted in much improved global control of stem rust. However, this situation has changed following the emergence of the highly virulent stem rust pathotype Ug99, first detected in Uganda in 1998, and now widespread in parts of Africa and the Middle East (Singh et al., 2015).
Genetic control is considered as the most effective and environmentally friendly strategy to control rust disease and involves breeding effective disease resistance genes into wheat cultivars. Many rust resistance genes have been identified genetically, and introgression into wheat lines is increasingly being facilitated by the development of robust molecular markers. However, the massive and complex genome of wheat presents major challenges for the isolation of individual genes. In the past 17 years, 24 rust resistance genes have been cloned using various strategies, with more than half of these (15) identified only in the last 4 years. This recent accelerated progress was made possible by (i) the public availability of the first high quality reference genome for wheat (Chinese Spring RefSeq v1.0) and (ii) the development of various approaches for reducing the genome complexity to allow targeted resequencing analyses. In particular, 10 out of the 15 genes cloned since 2016 were identified through pipelines involving Target-sequence Enrichment and Sequencing (TEnSeq). In this review, we briefly covered some general features of the wheat-rust pathosystem and the most recent progress in the area of cloning wheat rust resistance and rust fungal effector genes.
Rust Disease of Wheat
Rust fungi are one of the most diverse groups of plant pathogens, consisting of more than 120 genera and 6000 species (Duplessis et al., 2011). Studies of cereal rust diseases go back to Felice Fontana in 1767, who is considered to be the first person to provide a detailed description of cereal rusts and to recognize that rusts are caused by fungi (Fontana and Pirone, 1932; Chester, 1946). The primary causal agents of the wheat stem rust, leaf rust, and stripe rust diseases are Puccinia graminis Pers.: Pers. f. sp. tritici Erikss. & E. Henn (Pgt), P. triticina (syn. P. recondita Rob. ex Desm. f. sp. tritici) (Pt), and P. striiformis Westend. f. sp. tritici Erikss. & E. Henn. (Pst), respectively (Roelfs, 1985; Samborski, 1985; Stubbs, 1985; Knott, 1989; McIntosh et al., 1995). All belong to the genus Puccinia, family Pucciniaceae, order Pucciniales, class Teliomycetes, and phylum Basidiomycota, within the kingdom Fungi. Because of the economic importance of rust diseases, the causal agents are the most intensively studied plant pathogenic fungal species (McIntosh et al., 1995; Figueroa et al., 2020). Moreover, some important principles derived from the pathogenetic rust studies, for example, the gene-for-gene model (Flor, 1971), have found wide applications in other host–pathogen systems.
Rust pathogens are well-known to have great pathogenic variability, and the frequent emergence of new virulent strains that overcome resistance genes present in cultivated wheat varieties has hindered efforts to achieve durable resistance to these pathogens (Chen et al., 2014; Zhao et al., 2016). Studies of pathogenic variability in rust populations have indicated that new virulences can arise through the introduction of exotic genotypes, mutation in clonal lineages, sexual recombination and asexual hybridization (Park, 2007; Li et al., 2019). As an example, the Pgt population of Australia was postulated to be the result of four exotic incursions into the country since 1925 (Zwer et al., 1992; Zhang J. et al., 2017). Recent genetic studies have shown that three of these were derived from southern Africa and represent a single clonal lineage that was first described in the 1920s (Li et al., 2019; Visser et al., 2019). Subsequent to these introductions, stepwise mutations to overcome individual resistance genes have led to the divergence of numerous races with different pathotypes (Park, 2007).
Genetic Resistance of Wheat Rust Disease
Concept of Plant Innate Immunity
Long-term co-evolution between plants and their pathogens has equipped plants with a sophisticated multi-layered immune system to guard themselves against pest and pathogens (Andersen et al., 2018). The development of our understanding of the plant immune system is summarized in chronological order in Figure 1.
Figure 1. The milestones of understanding plant immune system in chronological order. Adapted from Andersen et al. (2018).
The plant immune system is often described in terms of two components, i.e., pattern-triggered immunity (PTI), activated by the recognition of microbial or pathogen-associated molecular patterns (MAMPs or PAMPs), and effector-triggered immunity (ETI), which encompasses “gene-for-gene” type of resistance (Jones and Dangl, 2006; Dodds and Rathjen, 2010). Bacterial flagellin and chitin are classic examples of MAMPs or PAMPs and are recognized by pattern recognition receptors (PRRs) such as Receptor-Like Kinase (RLK) and Receptor-Like Protein (RLP) type transmembrane receptor proteins. ETI is often based on the recognition of cytosolic effectors by immune receptors with a conserved nucleotide-binding domain (NB-ARC) and a leucine-rich repeat domain (LRR), hereafter referred to as NLRs. This type of resistance is normally associated with a hypersensitive response (HR) localized to infection sites. However, the distinction between PTI and ETI is often blurred, and recent models of plant immunity highlight that there are a range of cell-surface receptors recognizing extracellular ligands (which may be PAMPS or extracellular effectors or host derived-patterns), and intracellular receptors that recognize intracellular ligands (Cook et al., 2015; Kanyuka and Rudd, 2019).
Plant Resistance (R) Genes
More than 300 R genes have been cloned from plants to date and Kourelis and van der Hoorn (2018) reviewed the defense mechanisms functions and summarized nine main defense mechanisms based on all cloned R proteins. Most plant R genes are dominant in action and encode immune receptors that recognize pathogen avirulence proteins as described above. However, some genes that confer resistance phenotypes operate via different mechanisms. For instance, Hm1 from maize (Zea mays) was the first cloned R gene and encodes an enzyme that detoxifies a toxin from the fungal pathogen Cochliobolus carbonum (Johal and Briggs, 1992). Two wheat adult plant resistance (APR) genes Lr34/Yr18/Sr57/Pm38 and Lr67/Yr46/Sr55/Pm46 are also examples of non-immunity-mediated resistance genes (Krattinger et al., 2009; Moore et al., 2015). Both encode transporter proteins and confer resistance against multiple pathogens in wheat and can also function in other crops (Risk et al., 2013; Chauhan et al., 2015; Krattinger et al., 2016; Rinaldo et al., 2017; Schnippenkoetter et al., 2017; Sucher et al., 2018). These two genes have been important components of wheat breeding for rust resistance and often show additive or synergistic interactions with the other more “typical” R genes (immune receptors). Thus approaches for identifying effective resistance genes must consider both classical immune receptor class genes as well as other novel classes that may operate via different mechanisms.
Plant NLRs
Most of the cloned R genes from the wheat-rust pathosystems encode immune receptors of the NLR class (19 out of 24, Table 1). Despite the huge evolutionary distance between the plant and animal kingdoms, members of both use the intracellular proteins of the NB-ARC-LRR superfamily to perceive pathogens. However, recent studies suggest that this is the result of convergent evolution and that the domain architecture of the NLRs evolved at least twice (Urbach and Ausubel, 2017). Jones et al. (2016) reviewed the processes associated with intracellular innate immunity in both plants and animals and built an NLR tree to illustrate the proposed evolution of NLR genes following independent pathways for plant and animal species. The authors also proposed that plant and animal NLRs evolved from two distinct derivatives of a common ancestral prokaryotic adenosine triphosphatase (ATPase) represented by the NB-ARC domain class (nucleotide-binding domain shared by APAF-1, plant R proteins, and CED-4) and the NACHT domain class (shared by NAIP, CIITA, HET-E, and TP1). While animal and fungal genomes can contain both NB-ARC and NACHT domains, no NACHT domains have been found in plants. NACHT domains are also absent from some animal taxa such as the nematodes and Drosophila.
Table 1. Cloned rust resistance genes in wheat from 2003 to 2020.
The number of NLRs in a given plant genome can be as high as several 1000 (Appels et al., 2018). As revealed by the increasing number of newly available whole genome sequences and the more precise bioinformatic pipelines developed for identifying NLR genes, the number of NLR genes varies greatly between species. While the number of NLRs is normally proportional to the size of the genome, apple (Malus domestica) is an exception in possessing nearly 1,000 NLRs despite having a relatively small genome (740 Mb). In contrast, the number of NLRs in orchids species (Apostasia shenzhenica) with a genome size of 349 Mb was reported to be normally less than a 100 (Zhang G. Q. et al., 2017; Xue et al., 2020). This appears to indicate that the number of NLRs within a certain genome may also be a result of the selection pressures posed by pathogens during the evolutionary history of that plant lineage and degree of exposure to pathogens (Borrelli et al., 2018).
The Wheat Genome and Its NLRs
The bread wheat (Triticum aestivum) genome is one of the most challenging plant genomes to study. It is highly repetitive (∼85%) and approximately 15.4–15.8 Gbp in size, which is five times larger than the human genome (Appels et al., 2018). T. aestivum is a hexaploid species that arose through natural hybridization of three closely related wild grass species which contributed the A, B, and D genomes of wheat (Salamini et al., 2002). A low coverage survey sequence of the wheat genome became available in 2012 (Brenchley et al., 2012), after which assemblies have been greatly improved as sequencing technologies and bioinformatic analysis pipelines became more powerful. The first version of the wheat genome in the form of chromosome-sized scaffolds (IWGSC RefSeq v1.0) was made publicly available in 2018 (Appels et al., 2018) and an improved version v2.0 in 2019. High-quality reference genome sequences were also recently published for the wild diploid progenitor of the wheat D genome (Aegilops tauschii) and the wild tetraploid progenitor T. diccocoides and cultivated tetraploid wheat T. turgidum cv. durum (Avni et al., 2017; Luo et al., 2017; Maccaferri et al., 2019). Based on an analysis of the IWGSC RefSeq v1.0 assembly Steuernagel et al. (2018) reported a total of 3,400 full-length NLR loci.
TEnSeq Pipelines as a Catalyst in Isolating R Genes in Wheat
The massive and complex genome of wheat has made the isolation of individual genes a challenging task. Among ∼200 rust resistance genes cataloged in wheat, only a small number have been cloned and had their molecular functions studied. A complete list of the rust resistance genes cloned so far from wheat is shown in Table 1. In 17 years since 2003, there are in total 24 wheat rust R genes that have been cloned and published. More than half of the genes were cloned in the last 4 years after the MutRenSeq pipeline was published in 2016 (Steuernagel et al., 2016), and 10 of these were cloned directly or partially using approaches based on some version of genome complexity reduction which we refer to here collectively as TEnSeq (Target-sequence Enrichment and Sequencing).
Various gene cloning strategies have been used to position a gene to its exact location within the genome and identify its nucleotide sequence and the protein that it encodes. The traditional map-based or positional cloning strategy narrows down the gene location by using genetic recombination in biparental populations that segregate for the gene of interest (Keller et al., 2018) (Figure 2A). This approach depends on the availability of a large segregating population that allows mapping of the gene of interest to a small genetic interval. Markers from this genetic interval are then used to screen a physical library, often prepared in bacterial artificial chromosomes (BACs). The large genome of wheat and relatively small size of the BAC clones (100–200 kb) requires that the target gene is mapped to a very small genetic interval that corresponds to only a few overlapping BAC clones. Thus this strategy is not viable for target genes derived from wild relatives of wheat and which are located in introgressed genome segments that do not recombine with wheat chromatin. Applying this strategy on genes that are located in centromeric regions is also extremely challenging, as recombination rates in this region are low. It is not uncommon to spend 5–10 years or even longer on cloning one gene by map-based cloning (Mago et al., 2015; Klymiuk et al., 2018).
Figure 2. Illustrations of the five strategies applied for cloning resistance genes in wheat-rust pathosystem. (A) Map-based cloning; (B) MutRenSeq; (C) AgRenSeq; (D) MutChromSeq; (E) TACCA.
To overcome the limitations of the map-based cloning strategy in the large genome of wheat, alternative approaches were developed and validated by the rapid cloning of several genes using TEnSeq pipelines. These include MutRenSeq (Mutagenesis and the Resistance gene Enrichment and Sequencing) (Figure 2B), AgRenSeq (Association genetics with R gene enrichment Sequencing) (Figure 2C), MutChromSeq (Mutgenesis Chromosome flow sorting and short-read Sequencing) (Figure 2D), and TACCA (Targeted Chromosome-based Cloning via long-range Assembly) (Figure 2E). The common component of these approaches is the intent to reduce the genome complexity prior to the use of next-generation sequencing (NGS). MutRenSeq and AgRenSeq are based on NLR-targeted DNA capture by hybridization, while MutChromSeq and TACCA rely on the purification of individual chromosomes from wheat lines. Comparisons of these newly developed approaches and the classical map-based cloning strategy are outlined in Table 2.
Table 2. Different strategies for cloning rust resistance genes in wheat.
MutRenSeq was the earliest developed TEnSeq pipeline applied to clone rust resistance genes in the wheat genome. It is based on the Resistance gene Enrichment Sequencing (RenSeq) approach, which is an NLR gene-targeted, resistance gene enrichment and sequencing method. Jupe et al. (2013) described the RenSeq approach to improve the annotation of the NB-LRR gene repertoire from sequenced plant genomes, which are often poorly assembled due to the presence of complex gene families. As a proof of concept, they used it to reannotate and map NLRs in potato (Solanum tuberosum). The workflow was to first construct a customized target enrichment library (bait library) comprised of a series of 120-mer biotinylated RNA oligonucleotide. These were designed based on previously annotated NLR-like sequences from potato, tomato, tobacco, and pepper genomes and aimed to fully cover each NLR-gene sequence with bait probes of at least 80% sequence identity. This bait library was then used to capture and enrich the NB-LRR genes from a genomic DNA library of potato which was then sequenced using Illumina technology. Through this approach, the number of annotated NLRs in the potato genome was increased from 438 to 755. Furthermore, they applied RenSeq successfully in identifying SNP markers that co-segregate with the resistance against late blight pathogen Phytophthora infestans in two independent segregating populations of wild Solanum species.
Following the successful application of RenSeq on genetic mapping on the Solanum species, Steuernagel et al. (2016) proposed a three-step method, MutRenSeq, that combines mutagenesis with NLR gene capture and sequencing for rapid identification of resistance genes in wheat. Similarly, they designed a bait library containing 60,000 120-mer RNA probes with ≥ 95% similarity to predicted NLR genes present in the genome and transcriptome sequence data from Triticeae species including barley (Hordeum vulgare), hexaploid wheat (T. aestivum), tetraploid wheat (T. durum), red wild einkorn (T. urartu), domesticated einkorn (T. monococcum), and three goatgrass species (Ae. tauschii, Ae. sharonensis, and Ae. speltoides). The first successful application of MutRenSeq pipeline in cloning wheat rust resistance genes was the rapid cloning of wheat stem rust resistance genes Sr22 and Sr45. Later, its high efficiency was again demonstrated through the successful and rapid identification of Yr7, Yr5, and YrSP, which are the first three cloned major R genes against wheat stripe rust (Marchal et al., 2018). The noticeable advantage of the MutRenSeq compared with the classical map-based cloning method is that it obviates the need for a high resolution segregating mapping family and for building a physical library contig that covers the genetic interval. This not only reduces the time involved in these processes but also makes it an ideal approach for targeting genes that are located in low recombination regions such as those derived from alien species. This application is demonstrated by the recent identification of the stem rust resistance genes Sr26 and Sr61, which are located in a non-recombining introgressed segment from Thinopyrum ponticum (Zhang et al., under review).
The AgRenSeq pipeline also utilizes the NLR-gene capture method but integrates with Genome-Wide Association Studies (GWAS) to permit the cloning of R genes from a host diversity panel. Arora et al. (2019) validated this approach by identifying the wheat stem rust resistance genes Sr46 and SrTA1662 from a panel Ae. tauschii accessions. They also identified the previously cloned Sr33 and Sr45 genes within this diversity panel. In this approach, the GWAS analysis was based on the use of unique K-mer (sub-sequences) markers rather than single nucleotide polymorphisms. K-mers associated with resistance phenotypes were then used to identify the candidate NLRs.
MutChromSeq follows a similar basic principle to MutRenSeq of using multiple mutation events to identify candidate genes. However, the NLR gene-capture is replaced by whole chromosome isolation to reduce genome complexity. This approach requires that the chromosome location of the target gene is known so that the chromosome can be isolated from wild-type and mutant lines to allow sequence comparison. This approach circumvents one of the limitations of MutRenSeq, which is the requirement that a similar gene is represented in the bait library designed from existing pan-genome sequences. This approach was described first for cloning of the wheat powdery mildew resistance gene Pm2 by Sanchez-Martin et al. (2016) and was later used on cloning the leaf rust resistance gene Rph1 from barley by Dracatos et al. (2019). MutChromSeq also has the advantage that it does not rely on an underlying assumption that the resistance gene belongs to the NLR class, and therefore would be appropriate for identification of non-immune mediated resistance genes. Its most recent application is the cloning of Med15 encoded by SuSr-D1, a suppressor gene of stem rust resistance from the wheat cultivar ‘Canthatch’ (Hiebert et al., 2020).
The TACCA pipeline was first described by Thind et al. (2017) in the cloning of the leaf rust resistance gene Lr22a. It is essentially a map-based cloning strategy coupled with a cultivar specific chromosome assembly, which effectively increases the size of physical contigs onto which the genetic interval can be mapped. In this case, a sequence assembly was generated for chromosome 2D of a line carrying Lr22a after chromosome flow sorting and using a combination of Illumina short-read sequencing and proximity ligation of in vitro-reconstituted chromatin (Chicago long-range linkage). This allowed the identification of candidate genes within the physical interval delineated by the position of markers closely linked to Lr22a.
Identifying Wheat Rust Effectors
Fungal effectors are proteins secreted by pathogens that facilitate infection, often by suppressing plant immunity to help the invasion of the host (Uhse and Djamei, 2018). Fungal effectors may act in the host cytoplasm or apoplast and are mostly represented by small secreted proteins. Although 300 amino acids in size is commonly adopted as a size cut-off for effector prediction, some exceptions to this limit exist, notably AvrSr35 from Pgt, which is about 600 aa. The prediction of fungal effectors was facilitated in recent years by the increasing availability of fungal genome sequence data, and especially the development of approaches to assemble the two haploid genomes of these dikaryotic organisms separately in the case of rust fungi (Miller et al., 2018; Schwessinger et al., 2018; Li et al., 2019). It also benefit significantly from the rapid development of machine learning-based effector prediction tools (Schwessinger et al., 2018). For rust fungi, most known avirulence proteins are cytosolic effectors that are delivered into host cells during infection from specialized haustoria structures and recognized by intracellular NLR-type resistance proteins (Garnica et al., 2014).
In the wheat-rust pathosystem two rust effector/Avr genes have been identified: namely AvrSr50 and AvrSr35 from Pgt (Chen et al., 2017; Salcedo et al., 2017), which are recognized by the corresponding Sr50 and Sr35 resistance genes in wheat. These were identified based on whole-genome sequencing of wild-type (avirulent) Pgt isolates and virulent mutants. These two Pgt effectors are quite distinct from each other in sequence, but both are Haustorial Secreted Proteins (HSPs). Coincidentally, the two genes are located adjacent to each other in the Pgt genome (Li et al., 2019). As is often the case for fungal effectors, these two Avr genes are unique to P. graminis, with no homologs in related rust species. Allele mining of AvrSr50 and AvrSr35 has identified numerous genetic variants associated with either virulence or avirulence phenotypes which can serve as predictors of Pgt pathogenicity on wheat lines carrying these resistance genes. An expanded repertoire of identified rust Avr genes will ultimately lead to the prediction of isolate virulence profiles from sequence data with applications in field-based molecular diagnostics.
The Pursuit of Durable Resistance
Resistance Gene Stewardship and Deployment
The increasing number of cloned wheat rust resistance genes in recent years has led to a reconsideration of how to deploy these newly cloned resistance genes in order to escape from the “boom and bust cycle” (Figure 3). Effective resistance gene stewardship refers to the careful and responsible management of resistance genes with the aim of prolonging the resistance effect (Pretorius et al., 2017). Gene pyramids are commonly considered as the best approach to gene stewardship as this minimizes the chance of the pathogen acquiring virulence through mutation. Backcrosses and transgenic gene cassettes are two practical methods for combining multiple resistance genes into the same background.
Figure 3. The illustration of the resistance “Boom and Bust Cycle.”
Effective resistance gene stewardship requires reference to the mechanism of resistance conferred by the available genes and their ability to work in combination. Rust resistance genes may be either race-specific or non-race-specific (Periyannan et al., 2017). Race-specific resistance refers to the resistance that is effective against some but not all races within a pathogen formae specialis, and generally follows the gene-for-gene model, e.g., the resistance occurs only when a specific immune receptor (R gene) encounters its corresponding effector (Avr gene). Most NLRs confers race-specific resistance and are effective in all growth stages of the host, therefore, these are often termed as all stage resistance (ASR). Race non-specific resistance describes resistance that is effective against all races of a pathogen species and sometimes may also be effective against multiple pathogens. These are normally quantitative traits conferring partial resistance that is able to slow down disease development. Wheat stem rust resistance gene Sr2, stripe rust resistance gene Yr36, and three multi-pathogen resistance loci Lr34, Lr46, and Lr67, all fall into this race non-specific category (Ellis et al., 2014). Most of the race non-specific resistance genes are effective only at the adult plant stage of the host and are therefore often described as adult plant resistance (APR) genes. The successful cloning of Lr34, Yr36, and Lr67 since 2009 revealed these APRs encode an ABC transporter, a kinase-START protein, and a hexose transporter, respectively. They appear to each have their own resistance mechanism, function constitutively and often increase the basal level of resistance of the host, which is different from the recognition based NLRs. An important consideration is that some race-specific genes are also only effective at the adult plant stage, such as the NLR-encoding Lr22a (Thind et al., 2017).
Ellis et al. (2014) proposed that the most effective and durable means for genetic control of wheat rusts is the use of combinations of multiple broadly effective ASR and APR genes. Mundt (2018) in his latest review of durable resistance also suggested that resistance was likely to be more durable by pyramiding ASR genes into APR gene backgrounds. The APR genes Sr2, Lr34, and Yr36 all have been reported to have some additive effects when in combination with certain seedling resistance genes. Given the distinct resistance mechanisms between APRs and ASRs based on current knowledge, these synergistic or additive resistance phenotypes between the APR and ASR are likely the result of combining two different modes of resistance.
In terms of gene stewardship, Pretorius et al. (2017) suggested that resistance genes could be classified into three groups with different management strategies required. R genes belonging to Group one do not need stewardship. This includes all the slow-rusting APR genes and the genes have already been extensively distributed in breeding programs and current wheat cultivars. Group two contains R genes that are publicly accessible but are yet to be deployed widely. Stewardship of these genes was strongly encouraged, especially by being deployed together with at least one other ASR gene that has broad resistance spectrum, or several APR genes. R genes in Group three refer to newly identified R genes that have never been deployed but potentially have high economic value. For this group, a minimum of three effective genes should be incorporated in order to withstand at least a double mutation to virulence in a pathogen. Patents or material transfer agreements could be employed to facilitate gene stewardship to act effectively, especially for R genes belonging to Group three. Pretorius et al. (2017) also indicated that stewardship is actually a result of the whole agricultural system, as the stewardship chain is only as strong as the weakest link. For example, an unintended release of a cultivar with a single resistance gene is enough to put this gene at risk of being overcome by the pathogen under strong selection pressure, no matter how much effort has been invested by the breeding community in generating gene pyramids containing this gene.
Conclusion
There are roughly over 200 rust resistance genes that have been officially designated in wheat. Only a handful of these has been cloned to date, while many additional effective genes are expected to exist in wheat relatives. The advances in wheat rust resistance gene cloning reviewed in this paper, in particular, the TEnSeq pipeline which includes the MutRenSeq, MutChromSeq, TACCA, and AgRenSeq, will facilitate the identification of a much broader repertoire of wheat resistance genes. This will provide many more tools for marker-assisted selection in wheat breeding as well as the raw gene sequences to pursue gene stacking via transgenic gene cassettes. Together with advances in identifying genetic variation in rust Avr genes, these new tools should lead to more rational deployment strategies to maximize resistance durability. Although the TEnSeq strategies were initially developed specifically for resistance genes, these approaches can also be applied to genes with other functions and also adapted to other crops with large and complex genomes similar to wheat (Sanchez-Martin et al., 2016; Dracatos et al., 2019; Hiebert et al., 2020).
Author Contributions
JZ wrote the manuscript. All authors reviewed and provided edits.
Funding
JZ was supported by the CSIRO Early Research Career (CERC) postdoctoral fellowship program.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank the Grains Research and Development Corporation (GRDC), Australia, for their financial support.
References
Andersen, E. J., Ali, S., Byamukama, E., Yen, Y., and Nepal, M. P. (2018). Disease resistance mechanisms in plants. Genes 9:339. doi: 10.3390/genes9070339
Appels, R., Eversole, K., Feuillet, C., Keller, B., Rogers, J., Stein, N., et al. (2018). Shifting the limits in wheat research and breeding using a fully annotated reference genome. Science 361:eaar7191. doi: 10.1126/science.aar7191
Arora, S., Steuernagel, B., Gaurav, K., Chandramohan, S., Long, Y. M., Matny, O., et al. (2019). Resistance gene cloning from a wild crop relative by sequence capture and association genetics. Nat. Biotechnol. 37, 139–143. doi: 10.1038/s41587-018-0007-9
Avni, R., Nave, M., Barad, O., Baruch, K., Twardziok, S. O., Gundlach, H., et al. (2017). Wild emmer genome architecture and diversity elucidate wheat evolution and domestication. Science 357, 93–97. doi: 10.1126/science.aan0032
Beddow, J. M., Pardey, P. G., Chai, Y., Hurley, T. M., Kriticos, D. J., Braun, H. J., et al. (2015). Research investment implications of shifts in the global geography of wheat stripe rust. Nat. Plants 1:15132. doi: 10.1038/nplants.2015.132
Borrelli, G. M., Mazzucotelli, E., Marone, D., Crosatti, C., Michelotti, V., Vale, G., et al. (2018). Regulation and evolution of NLR genes: a close interconnection for plant immunity. Int. J. Mol. Sci. 19:1662. doi: 10.3390/ijms19061662
Brenchley, R., Spannagl, M., Pfeifer, M., Barker, G. L. A., D’Amore, R., Allen, A. M., et al. (2012). Analysis of the bread wheat genome using whole-genome shotgun sequencing. Nature. 491, 705–710. doi: 10.1038/nature11650
Chauhan, H., Boni, R., Bucher, R., Kuhn, B., Buchmann, G., Sucher, J., et al. (2015). The wheat resistance gene Lr34 results in the constitutive induction of multiple defense pathways in transgenic barley. Plant J. 84, 202–215. doi: 10.1111/tpj.13001
Chen, J., Upadhyaya, N. M., Ortiz, D., Sperschneider, J., Li, F., Bouton, C., et al. (2017). Loss of AvrSr50 by somatic exchange in stem rust leads to virulence for Sr50 resistance in wheat. Science 358, 1607–1610. doi: 10.1126/science.aao4810
Chen, S., Rouse, M. N., Zhang, W., Zhang, X., Guo, Y., Briggs, J., et al. (2019). Wheat gene Sr60 encodes a protein with two putative kinase domains that confers resistance to stem rust. New Phytol. 225, 948–959. doi: 10.1111/nph.16169
Chen, S., Zhang, W., Bolus, S., Rouse, M. N., and Dubcovsky, J. (2018). Identification and characterization of wheat stem rust resistance gene Sr21 effective against the Ug99 race group at high temperature. PLoS Genet. 14:e1007287. doi: 10.1371/journal.pgen.1007287
Chen, W., Wellings, C., Chen, X., Kang, Z., and Liu, T. (2014). Wheat stripe (yellow) rust caused by Puccinia striiformis f. sp. tritici. Mol. Plant Pathol. 15, 433–446. doi: 10.1111/mpp.12116
Chester, K. S. (1946). The Nature and Prevention of the Cereal Rusts as Exemplified in the Leaf Rust of Wheat. Waltham, MA: Chronica Botanica Co.
Cloutier, S., McCallum, B. D., Loutre, C., Banks, T. W., Wicker, T., Feuillet, C., et al. (2007). Leaf rust resistance gene Lr1, isolated from bread wheat (Triticum aestivum L.) is a member of the large psr567 gene family. Plant Mol. Biol. 65, 93–106. doi: 10.1007/s11103-007-9201-8
Cook, D. E., Mesarich, C. H., and Thomma, B. P. H. J. (2015). Understanding plant immunity as a surveillance system to detect invasion. Annu. Rev. Phytopathol. 53, 541–563. doi: 10.1146/annurev-phyto-080614-120114
Dodds, P. N., and Rathjen, J. P. (2010). Plant immunity: towards an integrated view of plant-pathogen interactions. Nat. Rev. Genet. 11, 539–548. doi: 10.1038/nrg2812
Dracatos, P. M., Bartos, J., Elmansour, H., Singh, D., Karafiatova, M., Zhang, P., et al. (2019). The Coiled-Coil NLR Rph1, confers leaf rust resistance in barley cultivar Sudan. Plant Physiol. 179, 1362–1372. doi: 10.1104/pp.18.01052
Duplessis, S., Cuomo, C. A., Lin, Y. C., Aerts, A., Tisserant, E., Veneault-Fourrey, C., et al. (2011). Obligate biotrophy features unraveled by the genomic analysis of rust fungi. Proc. Natl. Acad. Sci. U.S.A. 108, 9166–9171. doi: 10.1073/pnas.1019315108
Ellis, J. G., Lagudah, E. S., Spielmeyer, W., and Dodds, P. N. (2014). The past, present and future of breeding rust resistant wheat. Front. Plant Sci. 5:641. doi: 10.3389/fpls.2014.00641
Feuillet, C., Travella, S., Stein, N., Albar, L., Nublat, A., and Keller, B. (2003). Map-based isolation of the leaf rust disease resistance gene Lr10 from the hexaploid wheat (Triticum aestivum L.) genome. Proc. Natl. Acad. Sci. U.S.A. 100, 15253–15258. doi: 10.1073/pnas.2435133100
Figueroa, M., Dodds, P. N., and Henningsen, E. C. (2020). Evolution of virulence in rust fungi - multiple solutions to one problem. Curr. Opin. Plant Biol. 56, 20–27. doi: 10.1016/j.pbi.2020.02.007
Flor, H. H. (1971). Current status of the gene-for-gene concept. Annu. Rev. Phytopathol. 9, 275–296. doi: 10.1146/annurev-phyto-072910-095339
Fontana, F., and Pirone, P. P. (1932). Observations on the Rust of Grain. Washington, DC: Hayworth Printing Company.
Fu, D., Uauy, C., Distelfeld, A., Blechl, A., Epstein, L., Chen, X., et al. (2009). A kinase-START gene confers temperature-dependent resistance to wheat stripe rust. Science 323, 1357–1360. doi: 10.1126/science.1166289
Garnica, D. P., Nemri, A., Upadhyaya, N. M., Rathjen, J. P., and Dodds, P. N. (2014). The ins and outs of rust haustoria. PLoS Pathog. 10:e1004329. doi: 10.1371/journal.ppat.1004329
Hiebert, C. W., Moscou, M. J., Hewitt, T., Steuernagel, B., Hernández-Pinzón, I., Green, P., et al. (2020). Stem rust resistance in wheat is suppressed by a subunit of the mediator complex. Nat. Commun. 11:1123. doi: 10.1038/s41467-020-14937-2
Huang, L., Brooks, S. A., Li, W., Fellers, J. P., Trick, H. N., and Gill, B. S. (2003). Map-based cloning of leaf rust resistance gene Lr21 from the large and polyploid genome of bread wheat. Genetics 164, 655–664.
Huerta-Espino, J., Singh, R. P., German, S., McCallum, B. D., Park, R. F., Chen, W. Q., et al. (2011). Global status of wheat leaf rust caused by Puccinia triticina. Euphytica 179, 143–160. doi: 10.1007/s10681-011-0361-x
Johal, G. S., and Briggs, S. P. (1992). Reductase-activity encoded by the Hm1 disease resistance gene in maize. Science 258, 985–987. doi: 10.1126/science.1359642
Jones, J. D., and Dangl, J. L. (2006). The plant immune system. Nature 444, 323–329. doi: 10.1038/nature05286
Jones, J. D., Vance, R. E., and Dangl, J. L. (2016). Intracellular innate immune surveillance devices in plants and animals. Science 354:aaf6395. doi: 10.1126/science.aaf6395
Jupe, F., Witek, K., Verweij, W., Sliwka, J., Pritchard, L., Etherington, G. J., et al. (2013). Resistance gene enrichment sequencing (RenSeq) enables reannotation of the NB-LRR gene family from sequenced plant genomes and rapid mapping of resistance loci in segregating populations. Plant J. 76, 530–544. doi: 10.1111/tpj.12307
Kanyuka, K., and Rudd, J. J. (2019). Cell surface immune receptors: the guardians of the plant’s extracellular spaces. Curr. Opin. Plant Biol. 50, 1–8. doi: 10.1016/j.pbi.2019.02.005
Keller, B., Feuillet, C., and Yahiaoui, N. (2005). Map-based isolation of disease resistance genes from bread wheat: cloning in a supersize genome. Genet. Res. 85, 93–100. doi: 10.1017/S0016672305007391
Keller, B., Wicker, T., and Krattinger, S. G. (2018). Advances in wheat and pathogen genomics: implications for disease control. Annu. Rev. Phytopathol. 56, 67–87. doi: 10.1146/annurev-phyto-080516-035419
Klymiuk, V., Yaniv, E., Huang, L., Raats, D., Fatiukha, A., Chen, S. S., et al. (2018). Cloning of the wheat Yr15 resistance gene sheds light on the plant tandem kinase-pseudokinase family. Nat. Commun. 9:3735. doi: 10.1038/s41467-018-06138-9
Kourelis, J., and van der Hoorn, R. A. L. (2018). Defended to the nines: 25 years of resistance gene cloning identifies nine mechanisms for R protein function. Plant Cell 30, 285–299. doi: 10.1105/tpc.17.00579
Krattinger, S. G., Lagudah, E. S., Spielmeyer, W., Singh, R. P., Huerta-Espino, J., McFadden, H., et al. (2009). A putative ABC transporter confers durable resistance to multiple fungal pathogens in wheat. Science 323, 1360–1363. doi: 10.1126/science.1166453
Krattinger, S. G., Sucher, J., Selter, L. L., Chauhan, H., Zhou, B., Tang, M., et al. (2016). The wheat durable, multipathogen resistance gene Lr34 confers partial blast resistance in rice. Plant Biotechnol. J. 14, 1261–1268. doi: 10.1111/pbi.12491
Li, F., Upadhyaya, N. M., Sperschneider, J., Matny, O., Hoa, N. P., Mago, R., et al. (2019). Emergence of the Ug99 lineage of the wheat stem rust pathogen through somatic hybridisation. Nat. Commun. 10:5068. doi: 10.1038/s41467-019-12927-7
Luo, M. C., Gu, Y. Q., Puiu, D., Wang, H., Twardziok, S. O., Deal, K. R., et al. (2017). Genome sequence of the progenitor of the wheat D genome Aegilops tauschii. Nature 551, 498–502. doi: 10.1038/nature24486
Maccaferri, M., Harris, N. S., Twardziok, S. O., Pasam, R. K., Gundlach, H., Spannagl, M., et al. (2019). Durum wheat genome highlights past domestication signatures and future improvement targets. Nat. Genet. 51, 885–895. doi: 10.1038/s41588-019-0381-3
Mago, R., Zhang, P., Vautrin, S., Simkova, H., Bansal, U., Luo, M. C., et al. (2015). The wheat Sr50 gene reveals rich diversity at a cereal disease resistance locus. Nat. Plants 1:15186. doi: 10.1038/nplants.2015.186
Marchal, C., Zhang, J., Zhang, P., Fenwick, P., Steuernagel, B., Adamski, N. M., et al. (2018). BED-domain-containing immune receptors confer diverse resistance spectra to yellow rust. Nat. Plants 4, 662–668. doi: 10.1038/s41477-018-0236-4
McIntosh, R. A., Wellings, C. R., and Park, R. F. (1995). Wheat Rusts: An Atlas of Resistance Genes. Melbourne: CSIRO Publishing.
Miller, M. E., Zhang, Y., Omidvar, V., Sperschneider, J., Schwessinger, B., Raley, C., et al. (2018). De novo assembly and phasing of dikaryotic genomes from two isolates of Puccinia coronata f. sp avenae, the causal agent of oat crown rust. mBio 9:e01650-17. doi: 10.1128/mBio.01650-17
Moore, J. W., Herrera-Foessel, S., Lan, C., Schnippenkoetter, W., Ayliffe, M., Huerta-Espino, J., et al. (2015). A recently evolved hexose transporter variant confers resistance to multiple pathogens in wheat. Nat. Genet. 47, 1494–1498. doi: 10.1038/ng.3439
Mundt, C. C. (2018). Pyramiding for resistance durability: theory and practice. Phytopathology 108, 792–802. doi: 10.1094/PHYTO-12-17-0426-RVW
Murray, G. M., and Brennan, J. P. (2009). Estimating disease losses to the Australian wheat industry. Australas. Plant Pathol. 38, 558–570. doi: 10.1071/AP09053
Park, R. F. (2007). Stem rust of wheat in Australia. Aust. J. Agric. Res. 58, 558–566. doi: 10.1071/Ar07117
Periyannan, S., Milne, R. J., Figueroa, M., Lagudah, E. S., and Dodds, P. N. (2017). An overview of genetic rust resistance: from broad to specific mechanisms. PLoS Pathog. 13:e1006380. doi: 10.1371/journal.ppat.1006380
Periyannan, S., Moore, J., Ayliffe, M., Bansal, U., Wang, X., Huang, L., et al. (2013). The gene Sr33, an ortholog of barley Mla genes, encodes resistance to wheat stem rust race Ug99. Science 341, 786–788. doi: 10.1126/science.1239028
Pretorius, Z. A., Ayliffe, M., Bowden, R. L., Boyd, L. A., DePauw, R. M., Jin, Y., et al. (2017). “Advances in control of wheat rust,” in Achieving Sustainable Cultivation of Wheat: Ambridge: Breeding, Quality Traits, Pests and Diseases, Vol. 1, ed. P. Langridge (Cambridge: Burleigh Dodds Science Publishing Limited).
Rinaldo, A., Gilbert, B., Boni, R., Krattinger, S. G., Singh, D., Park, R. F., et al. (2017). The Lr34 adult plant rust resistance gene provides seedling resistance in durum wheat without senescence. Plant Biotechnol. J. 15, 894–905. doi: 10.1111/pbi.12684
Risk, J. M., Selter, L. L., Chauhan, H., Krattinger, S. G., Kumlehn, J., Hensel, G., et al. (2013). The wheat Lr34 gene provides resistance against multiple fungal pathogens in barley. Plant Biotechnol. J. 11, 847–854. doi: 10.1111/pbi.12077
Roelfs, A. P. (1978). Estimated Losses Caused by Rust in Small Grain Cereals in the United States-1918-1976: Miscellaneous Publication 1363. Washington, DC: United States Department of Agriculture.
Roelfs, A. P. (1985). “Wheat and rye stem rust,” in The Cereal Rusts, eds A. P. Roelfs and W. R. Bushnell (Orlando, FL: Academic press, Inc).
Saintenac, C., Zhang, W., Salcedo, A., Rouse, M. N., Trick, H. N., Akhunov, E., et al. (2013). Identification of wheat gene Sr35 that confers resistance to Ug99 stem rust race group. Science 341, 783–786. doi: 10.1126/science.1239022
Salamini, F., Ozkan, H., Brandolini, A., Schafer-Pregl, R., and Martin, W. (2002). Genetics and geography of wild cereal domestication in the Near East. Nat. Rev. Genet. 3, 429–441. doi: 10.1038/nrg817
Salcedo, A., Rutter, W., Wang, S., Akhunova, A., Bolus, S., Chao, S., et al. (2017). Variation in the AvrSr35 gene determines Sr35 resistance against wheat stem rust race Ug99. Science 358, 1604–1606. doi: 10.1126/science.aao7294
Samborski, D. J. (1985). “Wheat leaf rust,” in The Cereal Rusts, eds A. P. Roelfs and W. R. Bushnell (Orlando, FL: Academic press, Inc), 469–500.
Sanchez-Martin, J., Steuernagel, B., Ghosh, S., Herren, G., Hurni, S., Adamski, N., et al. (2016). Rapid gene isolation in barley and wheat by mutant chromosome sequencing. Genome Biol. 17:221.
Schnippenkoetter, W., Lo, C., Liu, G., Dibley, K., Chan, W. L., White, J., et al. (2017). The wheat Lr34 multipathogen resistance gene confers resistance to anthracnose and rust in sorghum. Plant Biotechnol. J. 15, 1387–1396. doi: 10.1111/pbi.12723
Schwessinger, B., Sperschneider, J., Cuddy, W. S., Garnica, D. P., Miller, M. E., Taylor, J. M., et al. (2018). A near-complete haplotype-phased genome of the dikaryotic wheat stripe rust fungus Puccinia striiformis f. sp tritici reveals high interhaplotype diversity. mBio 9:e02275-17. doi: 10.1128/mBio.02275-17
Singh, R. P., Hodson, D. P., Jin, Y., Lagudah, E. S., Ayliffe, M. A., Bhavani, S., et al. (2015). Emergence and spread of new races of wheat stem rust fungus: continued threat to food security and prospects of genetic control. Phytopathology 105, 872–884. doi: 10.1094/PHYTO-01-15-0030-FI
Steuernagel, B., Periyannan, S. K., Hernandez-Pinzon, I., Witek, K., Rouse, M. N., Yu, G., et al. (2016). Rapid cloning of disease-resistance genes in plants using mutagenesis and sequence capture. Nat. Biotechnol. 34, 652–655. doi: 10.1038/nbt.3543
Steuernagel, B., Witek, K., Krattinger, S. G., Ramirez-Gonzalez, R. H., Schoonbeek, H.-J., Yu, G., et al. (2018). Physical and transcriptional organisation of the bread wheat intracellular immune receptor repertoire. bioRxiv [Preprint]. doi: 10.1101/339424
Stubbs, R. W. (1985). “Stripe rust,” in The Cereal Rusts, eds A. P. Roelfs and W. R. Bushnel (Orlando, FL: Academic press, Inc).
Sucher, J., Menardo, F., Praz, C. R., Boni, R., Krattinger, S. G., and Keller, B. (2018). Transcriptional profiling reveals no response of fungal pathogens to the durable, quantitative Lr34 disease resistance gene of wheat. Plant Pathol. 67, 792–798. doi: 10.1111/ppa.12797
Thind, A. K., Wicker, T., Simkova, H., Fossati, D., Moullet, O., Brabant, C., et al. (2017). Rapid cloning of genes in hexaploid wheat using cultivar-specific long-range chromosome assembly. Nat. Biotechnol. 35, 793–796. doi: 10.1038/nbt.3877
Uhse, S., and Djamei, A. (2018). Effectors of plant-colonizing fungi and beyond. PLoS Pathog. 14:e1006992. doi: 10.1371/journal.ppat.1006992
Urbach, J. M., and Ausubel, F. M. (2017). The NBS-LRR architectures of plant R-proteins and metazoan NLRs evolved in independent events. Proc. Natl. Acad. Sci. U.S.A. 114, 1063–1068. doi: 10.1073/pnas.1619730114
Visser, B., Meyer, M., Park, R. F., Gilligan, C. A., Burgin, L. E., Hort, M. C., et al. (2019). Microsatellite analysis and urediniospore dispersal simulations support the movement of Puccinia graminis f. sp. tritici from Southern Africa to Australia. Phytopathology 109, 133–144. doi: 10.1094/Phyto-04-18-0110-R
Wellings, C. R., Boyd, L. A., and Chen, X. M. M. (2012). “Resistance to stripe rust in wheat: pathogen biology driving resistance breeding,” in Disease Resistance in Wheat, eds I. Sharma (London: CAB International), 63–83. doi: 10.1079/9781845938185.0000
Xue, J. Y., Zhao, T., Liu, Y., Zhang, Y. X., Zhang, G. Q., Chen, H. F., et al. (2020). Genome- wide analysis of the nucleotide binding site leucine-rich repeat genes of four orchids revealed extremely low numbers of disease resistance genes. Front. Genet. 10:1286. doi: 10.3389/fgene.2019.01286
Zhang, C., Huang, L., Zhang, H., Hao, Q., Lyu, B., Wang, M., et al. (2019). An ancestral NB-LRR with duplicated 3’UTRs confers stripe rust resistance in wheat and barley. Nat. Commun. 10:4023. doi: 10.1038/s41467-019-11872-9
Zhang, G. Q., Liu, K. W., Li, Z., Lohaus, R., Hsiao, Y. Y., Niu, S. C., et al. (2017). The Apostasia genome and the evolution of orchids. Nature 549, 379–383. doi: 10.1038/nature23897
Zhang, J., Zhang, P., Karaoglu, H., and Park, R. F. (2017). Molecular characterization of Australian isolates of Puccinia graminis f. sp. tritici supports long-term clonality but also reveals cryptic genetic variation. Phytopathology 107, 1032–1038. doi: 10.1094/PHYTO-09-16-0334-R
Zhang, W., Chen, S., Abate, Z., Nirmala, J., Rouse, M. N., and Dubcovsky, J. (2017). Identification and characterization of Sr13, a tetraploid wheat gene that confers resistance to the Ug99 stem rust race group. Proc. Natl. Acad. Sci. U.S.A. 114, E9483–E9492. doi: 10.1073/pnas.1706277114
Zhao, J., Wang, M., Chen, X., and Kang, Z. (2016). Role of alternate hosts in epidemiology and pathogen variation of cereal rusts. Annu. Rev. Phytopathol. 54, 207–228. doi: 10.1146/annurev-phyto-080615-095851
Keywords: wheat rust, TEnSeq, plant immunity, gene cloning, durable resistance
Citation: Zhang J, Zhang P, Dodds P and Lagudah E (2020) How Target-Sequence Enrichment and Sequencing (TEnSeq) Pipelines Have Catalyzed Resistance Gene Cloning in the Wheat-Rust Pathosystem. Front. Plant Sci. 11:678. doi: 10.3389/fpls.2020.00678
Received: 04 March 2020; Accepted: 30 April 2020;
Published: 26 May 2020.
Edited by:
Paul Christiaan Struik, Wageningen University & Research, NetherlandsReviewed by:
Xianming Chen, United States Department of Agriculture (USDA), United StatesSylvie Cloutier, Agriculture and Agri-Food Canada (AAFC), Canada
Copyright © 2020 Zhang, Zhang, Dodds and Lagudah. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Peng Zhang, cGVuZy56aGFuZ0BzeWRuZXkuZWR1LmF1; Peter Dodds, cGV0ZXIuZG9kZHNAY3Npcm8uYXU=; Evans Lagudah, ZXZhbnMubGFndWRhaEBjc2lyby5hdQ==