Disease Resistance Gene Analogs (RGAs) in Plants

Identification of anthracnose race 1 resistance loci in lentil by integrating linkage mapping and genome-wide association study

Anthracnose, caused byColletotrichum lentis, is a devastating disease of lentil (Lens culinaris Medik.) in western Canada. Growing resistant lentil cultivars is the most cost-effective and environmentally friendly approach to prevent seed yield losses that can exceed 70%. To identify loci conferring resistance to anthracnose race 1 in lentil, biparental quantitative trait loci (QTL) mapping of two recombinant inbred line (RIL) populations was integrated with a genome-wide association study (GWAS) using 200 diverse lentil accessions from a lentil diversity panel. A major-effect QTL (qAnt1.Lc-3) conferring resistance to race 1 was mapped to lentil chromosome 3 and colocated on the lentil physical map for both RIL populations. Clusters of candidate nucleotide-binding leucine-rich repeat (NB-LRR) and other defense-related genes were uncovered within the QTL region. A GWAS detected 14 significant single nucleotide polymorphism (SNP) markers associated with race 1 resistance on chromosomes 3, 4, 5, and 6. The most significant GWAS SNPs on chromosome 3 supported qAnt1.Lc-3 and delineated a region of 1.6 Mb containing candidate resistance genes. The identified SNP markers can be directly applied in marker-assisted selection (MAS) to accelerate the introgression of race 1 resistance in lentil breeding.

Transcriptome Analysis of a Cotton Cultivar Provides Insights into the Differentially Expressed Genes Underlying Heightened Resistance to the Devastating Verticillium Wilt

Cotton is an important economic crop worldwide. Verticillium wilt (VW) caused by Verticillium dahliae (V. dahliae) is a serious disease in cotton, resulting in massive yield losses and decline of fiber quality. Breeding resistant cotton cultivars is an efficient but elaborate method to improve the resistance of cotton against V. dahliae infection. However, the functional mechanism of several excellent VW resistant cotton cultivars is poorly understood at present. In our current study, we carried out RNA-seq to discover the differentially expressed genes (DEGs) in the roots of susceptible cotton Gossypium hirsutum cultivar Junmian 1 (J1) and resistant cotton G.hirsutum cultivar Liaomian 38 (L38) upon Vd991 inoculation at two time points compared with the mock inoculated control plants. The potential function of DEGs uniquely expressed in J1 and L38 was also analyzed by GO enrichment and KEGG pathway associations. Most DEGs were assigned to resistance-related functions. In addition, resistance gene analogues (RGAs) were identified and analyzed for their role in the heightened resistance of the L38 cultivar against the devastating Vd991. In summary, we analyzed the regulatory network of genes in the resistant cotton cultivar L38 during V. dahliae infection, providing a novel and comprehensive insight into VW resistance in cotton.

A new comprehensive annotation of leucine-rich repeat-containing receptors in rice

Oryza sativa (rice) plays an essential food security role for more than half of the world's population. Obtaining crops with high levels of disease resistance is a major challenge for breeders, especially today, given the urgent need for agriculture to be more sustainable. Plant resistance genes are mainly encoded by three large leucine-rich repeat (LRR)-containing receptor (LRR-CR) families: the LRR-receptor-like kinase (LRR-RLK), LRR-receptor-like protein (LRR-RLP) and nucleotide-binding LRR receptor (NLR). Using lrrprofiler, a pipeline that we developed to annotate and classify these proteins, we compared three publicly available annotations of the rice Nipponbare reference genome. The extended discrepancies that we observed for LRR-CR gene models led us to perform an in-depth manual curation of their annotations while paying special attention to nonsense mutations. We then transferred this manually curated annotation to Kitaake, a cultivar that is closely related to Nipponbare, using an optimized strategy. Here, we discuss the breakthrough achieved by manual curation when comparing genomes and, in addition to 'functional' and 'structural' annotations, we propose that the community adopts this approach, which we call 'comprehensive' annotation. The resulting data are crucial for further studies on the natural variability and evolution of LRR-CR genes in order to promote their use in breeding future resilient varieties.

Integrative Pre-Breeding for Biotic Resistance in Forest Trees

Climate change is unleashing novel biotic antagonistic interactions for forest trees that may jeopardize populations' persistence. Therefore, this review article envisions highlighting major opportunities from ecological evolutionary genomics to assist the identification, conservation, and breeding of biotic resistance in forest tree species. Specifically, we first discuss how assessing the genomic architecture of biotic stress resistance enables us to recognize a more polygenic nature for a trait typically regarded Mendelian, an expectation from the Fisherian runaway pathogen-host concerted arms-race evolutionary model. Secondly, we outline innovative pipelines to capture and harness natural tree pre-adaptations to biotic stresses by merging tools from the ecology, phylo-geography, and omnigenetics fields within a predictive breeding platform. Promoting integrative ecological genomic studies promises a better understanding of antagonistic co-evolutionary interactions, as well as more efficient breeding utilization of resistant phenotypes.

The whole-genome and expression profile analysis of WRKY and RGAs in Dactylis glomerata showed that DG6C02319.1 and DgWRKYs may cooperate in the immunity against rust

Orchardgrass (Dactylis glomerata) is one of the top four perennial forages worldwide and, despite its large economic advantages, often threatened by various environmental stresses. WRKY transcription factors (TFs) can regulate a variety of plant processes, widely participate in plant responses to biotic and abiotic stresses, and are one of the largest gene families in plants. WRKYs can usually bind W-box elements specifically. In this study, we identified a total of 93 DgWRKY genes and 281 RGAs, including 65, 169 and 47 nucleotide-binding site-leucine-rich repeats (NBS-LRRs), leucine-rich repeats receptor-like protein kinases (LRR-RLKs), and leucine-rich repeats receptor-like proteins (LRR-RLPs), respectively. Through analyzing the expression of DgWRKY genes in orchardgrass under different environmental stresses, it was found that many DgWRKY genes were differentially expressed under heat, drought, submergence, and rust stress. In particular, it was found that the greatest number of genes were differentially expressed under rust infection. Consistently, GO and KEGG enrichment analysis of all genes showed that 78 DgWRKY TFs were identified in the plant–pathogen interaction pathway, with 59 of them differentially expressed. Through cis-acting element prediction, 154 RGAs were found to contain W-box elements. Among them, DG6C02319.1 (a member of the LRR-RLK family) was identified as likely to interact with 14 DGWRKYs. Moreover, their expression levels in susceptible plants after rust inoculation were first up-regulated and then down-regulated, while those in the resistant plants were always up-regulated. In general, DgWRKYs responded to both biotic stress and abiotic stress. DgWRKYs and RGAs may synergistically respond to the response of orchardgrass to rust. This study provides meaningful insight into the molecular mechanisms of WRKY proteins in orchardgrass.

Transcriptome profiling in susceptible and tolerant rubber tree clones in response to cassiicolin Cas1, a necrotrophic effector from Corynespora cassiicola

Corynespora cassiicola, a fungal plant pathogen with a large host range, causes important damages in rubber tree (Hevea brasiliensis), in Asia and Africa. A small secreted protein named cassiicolin was previously identified as a necrotrophic effector required for the virulence of C. cassiicola in specific rubber tree clones. The objective of this study was to decipher the cassiicolin-mediated molecular mechanisms involved in this compatible interaction. We comparatively analyzed the RNA-Seq transcriptomic profiles of leaves treated or not with the purified cassiicolin Cas1, in two rubber clones: PB260 (susceptible) and RRIM600 (tolerant). The reads were mapped against a synthetic transcriptome composed of all available transcriptomic references from the two clones. Genes differentially expressed in response to cassiicolin Cas1 were identified, in each clone, at two different time-points. After de novo annotation of the synthetic transcriptome, we analyzed GO enrichment of the differentially expressed genes in order to elucidate the main functional pathways impacted by cassiicolin. Cassiicolin induced qualitatively similar transcriptional modifications in both the susceptible and the tolerant clones, with a strong negative impact on photosynthesis, and the activation of defense responses via redox signaling, production of pathogenesis-related protein, or activation of the secondary metabolism. In the tolerant clone, transcriptional reprogramming occurred earlier but remained moderate. By contrast, the susceptible clone displayed a late but huge transcriptional burst, characterized by massive induction of phosphorylation events and all the features of a hypersensitive response. These results confirm that cassiicolin Cas1 is a necrotrophic effector triggering a hypersensitive response in susceptible rubber clones, in agreement with the necrotrophic-effector-triggered susceptibility model.

Genes Associated with Foliar Resistance to Septoria Nodorum Blotch of Hexaploid Wheat (Triticum aestivum L.)

The genetic control of host response to the fungal necrotrophic disease Septoria nodorum blotch (SNB) in bread wheat is complex, involving many minor genes. Quantitative trait loci (QTL) controlling SNB response were previously identified on chromosomes 1BS and 5BL. The aim of this study, therefore, was to align and compare the genetic map representing QTL interval on 1BS and 5BS with the reference sequence of wheat and identify resistance genes (R-genes) associated with SNB response. Alignment of QTL intervals identified significant genome rearrangements on 1BS between parents of the DH population EGA Blanco, Millewa and the reference sequence of Chinese Spring with subtle rearrangements on 5BL. Nevertheless, annotation of genomic intervals in the reference sequence were able to identify and map 13 and 12 R-genes on 1BS and 5BL, respectively. R-genes discriminated co-located QTL on 1BS into two distinct but linked loci. NRC1a and TFIID mapped in one QTL on 1BS whereas RGA and Snn1 mapped in the linked locus and all were associated with SNB resistance but in one environment only. Similarly, Tsn1 and WK35 were mapped in one QTL on 5BL with NETWORKED 1A and RGA genes mapped in the linked QTL interval. This study provided new insights on possible biochemical, cellular and molecular mechanisms responding to SNB infection in different environments and also addressed limitations of using the reference sequence to identify the full complement of functional R-genes in modern varieties.

Whole-genome assembly of Corylus avellana cv'Tonda Gentile delle Langhe' using linked-reads (10X Genomics)

The European hazelnut (Corylus avellana L.; 2n = 2x = 22) is a worldwide economically important tree nut that is cross-pollinated due to sporophytic incompatibility. Therefore, any individual plant is highly heterozygous. Cultivars are clonally propagated using mound layering, rooted suckers, and micropropagation. In recent years, the interest in this crop has increased, due to a growing demand related to the recognized health benefits of nut consumption. C. avellana cv “Tonda Gentile delle Langhe” (“TGdL”) is well-known for its high kernel quality, and the premium price paid for this cultivar is an economic benefit for producers in northern Italy. Assembly of a high-quality genome is a difficult task in many plant species because of the high level of heterozygosity. We assembled a chromosome-level genome sequence of “TGdL” with a two-step approach. First, 10X Genomics Chromium Technology was used to create a high-quality sequence, which was then assembled into scaffolds with cv “Tombul” genome as the reference. Eleven pseudomolecules were obtained, corresponding to 11 chromosomes. A total of 11,046 scaffolds remained unplaced, representing 11% of the genome (46,504,161 bp). Gene prediction, performed with Maker-P software, identified 27,791 genes (AED ≤0.4 and 92% of BUSCO completeness), whose function was analyzed with BlastP and InterProScan software. To characterize “TGdL” specific genetic mechanisms, Orthofinder was used to detect orthologs between hazelnut and closely related species. The “TGdL” genome sequence is expected to be a powerful tool to understand hazelnut genetics and allow detection of markers/genes for important traits to be used in targeted breeding programs.

Genome-Wide Identification and Analysis of CC-NBS-LRR Family in Response to Downy Mildew and Black Rot in Chinese Cabbage

The nucleotide-binding site-leucine-rich repeat (NBS-LRR) gene family is the largest group of plant disease resistance (R) genes widespread in response to viruses, bacteria, and fungi usually involved in effector triggered immunity (ETI). Forty members of the Chinese cabbage CC type NBS-LRR family were investigated in this study. Gene and protein characteristics, such as distributed locations on chromosomes and gene structures, were explored through comprehensive analysis. CC-NBS-LRR proteins were classified according to their conserved domains, and the phylogenetic relationships of CC-NBS-LRR proteins in Brassica rapa, Arabidopsis thaliana, and Oryza sativa were compared. Moreover, the roles of BrCC-NBS-LRR genes involved in pathogenesis-related defense were studied and analyzed. First, the expression profiles of BrCC-NBS-LRR genes were detected by inoculating with downy mildew and black rot pathogens. Second, sensitive and resistant Chinese cabbage inbred lines were screened by downy mildew and black rot. Finally, the differential expression levels of BrCC-NBS-LRR genes were monitored at 0, 1, 3, 6, 12 and 24 h for short and 0, 3, 5, 7, 10 and 14 days for long inoculation time. Our study provides information on BrCC-NBS-LRR genes for the investigation of the functions and mechanisms of CC-NBS-LRR genes in Chinese cabbage.

Isolation of the ZmERS4 Gene From Maize and Its Functional Analysis in Transgenic Plants

A gene encoding a protein similar to ethylene receptor was isolated from maize (Zea mays L.), which was named as ZmERS4.The gene was 1,905 bp in length with an open reading frame that encoded a protein consisting of 634 amino acids. The homologous analysis showed that ZmERS4 shared high similarity with the ethylene receptor protein, OsERS1, from rice (Oryza sativa L.). ZmERS4 grouped into the ETR1 subfamily of ethylene receptors based on its conserved domain and phylogenetic status. Tissue-specific and induced expression analyses indicated that ZmERS4 was differentially expressed in maize tissues, predominantly in the stems and leaves, and was induced by salicylic acid (SA). Overexpression of ZmERS4 in Arabidopsis improved resistance against the bacterial pathogen, PstDC3000, by inducing the expression of SA signaling-related genes. Moreover, treatment with flg22 induced the expression of the defense-related gene, PR1, in maize protoplasts that transiently expressed ZmERS4. Furthermore, the ultra-high-performance liquid chromatography (UPLC) analysis showed that the SA contents in ZmERS4-overexpressing Arabidopsis lines were significantly higher than the control lines. Additionally, the improved resistance of ZmERS4-overexpressing Arabidopsis against PstDC3000 was blocked after pretreatment with the SA biosynthetic inhibitor, ABT. Based on the collective findings, we hypothesize that ZmERS4 plays an important role in disease resistance through SA-mediated signaling pathways.

Transcriptional Regulatory Networks Associate with Early Stages of Potato Virus X Infection of Solanum tuberosum

Potato virus X (PVX) belongs to genus Potexvirus. This study characterizes the cellular transcriptome responses to PVX infection in Russet potato at 2 and 3 days post infection (dpi). Among the 1242 differentially expressed genes (DEGs), 268 genes were upregulated, and 37 genes were downregulated at 2 dpi while 677 genes were upregulated, and 265 genes were downregulated at 3 dpi. DEGs related to signal transduction, stress response, and redox processes. Key stress related transcription factors were identified. Twenty-five pathogen resistance gene analogs linked to effector triggered immunity or pathogen-associated molecular pattern (PAMP)-triggered immunity were identified. Comparative analysis with Arabidopsis unfolded protein response (UPR) induced DEGs revealed genes associated with UPR and plasmodesmata transport that are likely needed to establish infection. In conclusion, this study provides an insight on major transcriptional regulatory networked involved in early response to PVX infection and establishment.

Comparison and evolutionary analysis of Brassica nucleotide binding site leucine rich repeat (NLR) genes and importance for disease resistance breeding

The Brassica genus contains many agriculturally significant oilseed and vegetable crops, however the crop yield is threatened by a range of fungal and bacterial pathogens. Nucleotide Binding Site Leucine Rich Repeat (NLR) genes play important roles in plant innate immunity. The evolution of NLR genes is influenced by genomic processes and pathogen selection. At the whole genome level, whole genome duplications (WGDs) generate abundant gene copies, most of which are lost during genome fractionation. At sub-genomic levels, some retained copies undergo duplication forming clusters which facilitate rapid evolution through recombination. The number, distribution and genetic variations of the NLR genes vary among Brassica species and within populations suggesting differential selection pressure exerted by pathogen populations throughout the evolutionary history. A study of the evolution of disease resistance genes in agriculturally important plants such as Brassicas helps gain insights into their function and inform the identification of resistance genes for breeding of resistant lines.

Transcriptomic Analysis of Wheat Seedling Responses to the Systemic Acquired Resistance Inducer N-Hydroxypipecolic Acid

The fungal pathogen Fusarium graminearum can cause destructive diseases on wheat, such as Fusarium head blight and Fusarium crown rot. However, a solution is still unavailable. Recently, N-hydroxypipecolic acid (NHP) was identified as a potent signaling molecule that is capable of inducing systemic acquired resistance to bacterial, oomycete, and fungal infection in several plant species. However, it is not clear whether NHP works in wheat to resist F. graminearum infection or how NHP affects wheat gene expression. In this report, we showed that pretreatment with NHP moderately increased wheat seedling resistance to F. graminearum. Using RNA sequencing, we found that 17% of wheat-expressed genes were significantly affected by NHP treatment. The genes encoding nucleotide-binding leucine-rich repeat immune receptors were significantly overrepresented in the group of genes upregulated by NHP treatment, while the genes encoding receptor-like kinases were not. Our results suggested that NHP treatment sensitizes a subset of the immune surveillance system in wheat seedlings, thereby facilitating wheat defense against F. graminearum infection.

Orthology and synteny analysis of receptor-like kinases "RLK" and receptor-like proteins "RLP" in legumes

BACKGROUND

Legume species are an important plant model because of their protein-rich physiology. The adaptability and productivity of legumes are limited by major biotic and abiotic stresses. Responses to these stresses directly involve plasma membrane receptor proteins known as receptor-like kinases and receptor-like proteins. Evaluating the homology relations among RLK and RLP for seven legume species, and exploring their presence among synteny blocks allow an increased understanding of evolutionary relations, physical position, and chromosomal distribution in related species and their shared roles in stress responses.

RESULTS

Typically, a high proportion of RLK and RLP legume proteins belong to orthologous clusters, which is confirmed in this study, where between 66 to 90% of the RLKs and RLPs per legume species were classified in orthologous clusters. One-third of the evaluated syntenic blocks had shared RLK/RLP genes among both legumes and non-legumes. Among the legumes, between 75 and 98% of the RLK/RLP were present in syntenic blocks. The distribution of chromosomal segments between Phaseolus vulgaris and Vigna unguiculata, two species that diverged ~ 8 mya, were highly similar. Among the RLK/RLP synteny clusters, seven experimentally validated resistance RLK/RLP genes were identified in syntenic blocks. The RLK resistant genes FLS2, BIR2, ERECTA, IOS1, and AtSERK1 from Arabidopsis and SLSERK1 from Solanum lycopersicum were present in different pairwise syntenic blocks among the legume species. Meanwhile, only the LYM1- RLP resistant gene from Arabidopsis shared a syntenic blocks with Glycine max.

CONCLUSIONS

The orthology analysis of the RLK and RLP suggests a dynamic evolution in the legume family, with between 66 to 85% of RLK and 83 to 88% of RLP belonging to orthologous clusters among the species evaluated. In fact, for the 10-species comparison, a lower number of singleton proteins were reported among RLP compared to RLK, suggesting that RLP positions are more physically conserved compared to RLK. The identification of RLK and RLP genes among the synteny blocks in legumes revealed multiple highly conserved syntenic blocks on multiple chromosomes. Additionally, the analysis suggests that P. vulgaris is an appropriate anchor species for comparative genomics among legumes.

Transcriptome Analysis of the Fruit of Two Strawberry Cultivars "Sunnyberry" and "Kingsberry" That Show Different Susceptibility to Botrytis cinerea after Harvest

Gray mold (Botrytis cinerea) is a fungal plant pathogen causing postharvest decay in strawberry fruit. Here, we conducted a comparative transcriptome analysis to identify differences in gene expression between the immature-green (IG) and mature-red (MR) stages of the “Sunnyberry” (gray mold-resistant) and “Kingsberry” (gray mold susceptible) strawberry cultivars. Most of the genes involved in lignin and alkane-type wax biosynthesis were relatively upregulated in “Sunnyberry”. However, pathogenesis-related proteins encoding R- and antioxidant-related genes were comparatively upregulated in “Kingsberry”. Analysis of gene expression and physiological traits in the presence and absence of B. cinerea inoculation revealed that the defense response patterns significantly differed between IG and MR rather than the cultivars. “Kingsberry” showed higher antioxidant induction at IG and upregulated hemicellulose-strengthening and R genes at MR. Hence, “Sunnyberry” and “Kingsberry” differed mainly in terms of the expression levels of the genes forming cuticle, wax, and lignin and controlling the defense responses. These discrepancies might explain the relative difference between these strawberry cultivars in terms of their postharvest responses to B. cinerea.

Association of differentially expressed R-gene candidates with leaf spot resistance in peanut (Arachis hypogaea L.)

Early leaf spot (ELS) and late leaf spot (LLS) are major fungal diseases of peanut that can severely reduce yield and quality. Development of acceptable genetic resistance has been difficult due to a strong environmental component and many major and minor QTLs. Resistance genes (R-genes) are an important component of plant immune system and have been identified in peanut. Association of specific R-genes to leaf spot resistance will provide molecular targets for marker-assisted breeding strategies. In this study, advanced breeding lines from different pedigrees were evaluated for leaf spot resistance and 76 candidate R-genes expression study was applied to susceptible and resistant lines. Thirty-six R-genes were differentially expressed and significantly correlated with resistant lines, of which a majority are receptor like kinases (RLKs) and receptor like proteins (RLPs) that sense the presence of pathogen at the cell surface and initiate protection response. The largest group was receptor-like cytoplasmic kinases (RLCKs) VII that are involved in pattern-triggered kinase signaling resulting in the production reactive oxygen species (ROS). Four R-genes were homologous to TMV resistant protein N which has shown to confer resistance against tobacco mosaic virus (TMV). When mapped to peanut genomes, 36 R-genes were represented in most chromosomes except for A09 and B09. Low levels of gene-expression in resistant lines suggest expression is tightly controlled to balance the cost of R-gene expression to plant productively. Identification and association of R-genes involved in leaf spot resistance will facilitate genetic selection of leaf spot resistant lines with good agronomic traits.

Recent Findings Unravel Genes and Genetic Factors Underlying Leptosphaeria maculans Resistance in Brassica napus and Its Relatives

Among the Brassica oilseeds, canola (Brassica napus) is the most economically significant globally. However, its production can be limited by blackleg disease, caused by the fungal pathogen Lepstosphaeria maculans. The deployment of resistance genes has been implemented as one of the key strategies to manage the disease. Genetic resistance against blackleg comes in two forms: qualitative resistance, controlled by a single, major resistance gene (R gene), and quantitative resistance (QR), controlled by numerous, small effect loci. R-gene-mediated blackleg resistance has been extensively studied, wherein several genomic regions harbouring R genes against L. maculans have been identified and three of these genes were cloned. These studies advance our understanding of the mechanism of R gene and pathogen avirulence (Avr) gene interaction. Notably, these studies revealed a more complex interaction than originally thought. Advances in genomics help unravel these complexities, providing insights into the genes and genetic factors towards improving blackleg resistance. Here, we aim to discuss the existing R-gene-mediated resistance, make a summary of candidate R genes against the disease, and emphasise the role of players involved in the pathogenicity and resistance. The comprehensive result will allow breeders to improve resistance to L. maculans, thereby increasing yield.

Transcriptome Analysis Identified Coordinated Control of Key Pathways Regulating Cellular Physiology and Metabolism upon Aspergillus flavus Infection Resulting in Reduced Aflatoxin Production in Groundnut

Aflatoxin-affected groundnut or peanut presents a major global health issue to both commercial and subsistence farming. Therefore, understanding the genetic and molecular mechanisms associated with resistance to aflatoxin production during host–pathogen interactions is crucial for breeding groundnut cultivars with minimal level of aflatoxin contamination. Here, we performed gene expression profiling to better understand the mechanisms involved in reduction and prevention of aflatoxin contamination resulting from Aspergillus flavus infection in groundnut seeds. RNA sequencing (RNA-Seq) of 16 samples from different time points during infection (24 h, 48 h, 72 h and the 7th day after inoculation) in U 4-7-5 (resistant) and JL 24 (susceptible) genotypes yielded 840.5 million raw reads with an average of 52.5 million reads per sample. A total of 1779 unique differentially expressed genes (DEGs) were identified. Furthermore, comprehensive analysis revealed several pathways, such as disease resistance, hormone biosynthetic signaling, flavonoid biosynthesis, reactive oxygen species (ROS) detoxifying, cell wall metabolism and catabolizing and seed germination. We also detected several highly upregulated transcription factors, such as ARF, DBB, MYB, NAC and C2H2 in the resistant genotype in comparison to the susceptible genotype after inoculation. Moreover, RNA-Seq analysis suggested the occurrence of coordinated control of key pathways controlling cellular physiology and metabolism upon A. flavus infection, resulting in reduced aflatoxin production.

Estimating individuals' genetic and non-genetic effects underlying infectious disease transmission from temporal epidemic data

Individuals differ widely in their contribution to the spread of infection within and across populations. Three key epidemiological host traits affect infectious disease spread: susceptibility (propensity to acquire infection), infectivity (propensity to transmit infection to others) and recoverability (propensity to recover quickly). Interventions aiming to reduce disease spread may target improvement in any one of these traits, but the necessary statistical methods for obtaining risk estimates are lacking. In this paper we introduce a novel software tool called SIRE (standing for "Susceptibility, Infectivity and Recoverability Estimation"), which allows for the first time simultaneous estimation of the genetic effect of a single nucleotide polymorphism (SNP), as well as non-genetic influences on these three unobservable host traits. SIRE implements a flexible Bayesian algorithm which accommodates a wide range of disease surveillance data comprising any combination of recorded individual infection and/or recovery times, or disease diagnostic test results. Different genetic and non-genetic regulations and data scenarios (representing realistic recording schemes) were simulated to validate SIRE and to assess their impact on the precision, accuracy and bias of parameter estimates. This analysis revealed that with few exceptions, SIRE provides unbiased, accurate parameter estimates associated with all three host traits. For most scenarios, SNP effects associated with recoverability can be estimated with highest precision, followed by susceptibility. For infectivity, many epidemics with few individuals give substantially more statistical power to identify SNP effects than the reverse. Importantly, precise estimates of SNP and other effects could be obtained even in the case of incomplete, censored and relatively infrequent measurements of individuals' infection or survival status, albeit requiring more individuals to yield equivalent precision. SIRE represents a new tool for analysing a wide range of experimental and field disease data with the aim of discovering and validating SNPs and other factors controlling infectious disease transmission.

Identification of QTLs conferring resistance to scald (Rhynchosporium commune) in the barley nested association mapping population HEB-25

BACKGROUND

Barley scald, caused by the fungus Rhynchosporium commune, is distributed worldwide to all barley growing areas especially in cool and humid climates. Scald is an economically important leaf disease resulting in yield losses of up to 40%. To breed resistant cultivars the identification of quantitative trait loci (QTLs) conferring resistance to scald is necessary. Introgressing promising resistance alleles of wild barley is a way to broaden the genetic basis of scald resistance in cultivated barley. Here, we apply nested association mapping (NAM) to map resistance QTLs in the barley NAM population HEB-25, comprising 1420 lines in BC1S3 generation, derived from crosses of 25 wild barley accessions with cv. Barke.

RESULTS

In scald infection trials in the greenhouse variability of resistance across and within HEB-25 families was found. NAM based on 33,005 informative SNPs resulted in the identification of eight reliable QTLs for resistance against scald with most wild alleles increasing resistance as compared to cv. Barke. Three of them are located in the region of known resistance genes and two in the regions of QTLs, respectively. The most promising wild allele was found at Rrs17 in one specific wild donor. Also, novel QTLs with beneficial wild allele effects on scald resistance were detected.

CONCLUSIONS

To sum up, wild barley represents a rich resource for scald resistance. As the QTLs were linked to the physical map the identified candidate genes will facilitate cloning of the scald resistance genes. The closely linked flanking molecular markers can be used for marker-assisted selection of the respective resistance genes to integrate them in elite cultivars.

Resistance Gene Analogs in the Brassicaceae: Identification, Characterization, Distribution, and Evolution

The Brassicaceae consists of a wide range of species, including important Brassica crop species and the model plant Arabidopsis (Arabidopsis thaliana). Brassica spp. crop diseases impose significant yield losses annually. A major way to reduce susceptibility to disease is the selection in breeding for resistance gene analogs (RGAs). Nucleotide binding site-leucine rich repeats (NLRs), receptor-like kinases (RLKs), and receptor-like proteins (RLPs) are the main types of RGAs; they contain conserved domains and motifs and play specific roles in resistance to pathogens. Here, all classes of RGAs have been identified using annotation and assembly-based pipelines in all available genome annotations from the Brassicaceae, including multiple genome assemblies of the same species where available (total of 32 genomes). The number of RGAs, based on genome annotations, varies within and between species. In total 34,065 RGAs were identified, with the majority being RLKs (21,691), then NLRs (8,588) and RLPs (3,786). Analysis of the RGA protein sequences revealed a high level of sequence identity, whereby 99.43% of RGAs fell into several orthogroups. This study establishes a resource for the identification and characterization of RGAs in the Brassicaceae and provides a framework for further studies of RGAs for an ultimate goal of assisting breeders in improving resistance to plant disease.

Computational identification of receptor-like kinases "RLK" and receptor-like proteins "RLP" in legumes

Background

In plants, the plasma membrane is enclosed by the cell wall and anchors RLK and RLP proteins, which play a fundamental role in perception of developmental and environmental cues and are crucial in plant development and immunity. These plasma membrane receptors belong to large gene/protein families that are not easily classified computationally. This detailed analysis of these plasma membrane proteins brings a new source of information to the legume genetic, physiology and breeding research communities.

Results

A computational approach to identify and classify RLK and RLP proteins is presented. The strategy was evaluated using experimentally-validated RLK and RLP proteins and was determined to have a sensitivity of over 0.85, a specificity of 1.00, and a Matthews correlation coefficient of 0.91. The computational approach can be used to develop a detailed catalog of plasma membrane receptors (by type and domains) in several legume/crop species. The exclusive domains identified in legumes for RLKs are WaaY, APH Pkinase_C, LRR_2, and EGF, and for RLP are L-lectin LPRY and PAN_4. The RLK-nonRD and RLCK subclasses are also discovered by the methodology. In both classes, less than 20% of the total RLK predicted for each species belong to this class. Among the 10-species evaluated ~ 40% of the proteins in the kinome are RLKs. The exclusive legume domain combinations identified are B-Lectin/PR5K domains in G. max, M. truncatula, V. angularis, and V. unguiculata and a three-domain combination B-lectin/S-locus/WAK in C. cajan, M. truncatula, P. vulgaris, V. angularis. and V. unguiculata.

Conclusions

The analysis suggests that about 2% of the proteins of each genome belong to the RLK family and less than 1% belong to RLP family. Domain diversity combinations are greater for RLKs compared with the RLP proteins and LRR domains, and the dual domain combination LRR/Malectin were the most frequent domain for both groups of plasma membrane receptors among legume and non-legume species. Legumes exclusively show Pkinase extracellular domains, and atypical domain combinations in RLK and RLP compared with the non-legumes evaluated. The computational logic approach is statistically well supported and can be used with the proteomes of other plant species.

Whole genome resequencing of four Italian sweet pepper landraces provides insights on sequence variation in genes of agronomic value

Sweet pepper (Capsicum annuum L.) is a high value crop and one of the most widely grown vegetables belonging to the Solanaceae family. In addition to commercial varieties and F₁ hybrids, a multitude of landraces are grown, whose genetic combination is the result of hundreds of years of random, environmental, and farmer selection. High genetic diversity exists in the landrace gene pool which however has scarcely been studied, thus bounding their cultivation. We re-sequenced four pepper inbred lines, within as many Italian landraces, which representative of as many fruit types: big sized blocky with sunken apex ('Quadrato') and protruding apex or heart shaped ('Cuneo'), elongated ('Corno') and smaller sized sub-spherical ('Tumaticot'). Each genomic sequence was obtained through Illumina platform at coverage ranging from 39 to 44×, and reconstructed at a chromosome scale. About 35.5k genes were predicted in each inbred line, of which 22,017 were shared among them and the reference genome (accession 'CM334'). Distinctive variations in miRNAs, resistance gene analogues (RGAs) and susceptibility genes (S-genes) were detected. A detailed survey of the SNP/Indels occurring in genes affecting fruit size, shape and quality identified the highest frequencies of variation in regulatory regions. Many structural variations were identified as presence/absence variations (PAVs), notably in resistance gene analogues (RGAs) and in the capsanthin/capsorubin synthase (CCS) gene. The large allelic diversity observed in the four inbred lines suggests their potential use as a pre-breeding resource and represents a one-stop resource for C. annuum genomics and a key tool for dissecting the path from sequence variation to phenotype.

Breeding Wheat for Durable Leaf Rust Resistance in Southern Africa: Variability, Distribution, Current Control Strategies, Challenges and Future Prospects

Leaf or brown rust of wheat caused by Puccinia triticina (Pt) is one of the most damaging diseases globally. Considerable progress has been made to control leaf rust through crop protection chemicals and host plant resistance breeding in southern Africa. However, frequent changes in the pathogen population still present a major challenge to achieve durable resistance. Disease surveillance and monitoring of the pathogen have revealed the occurrence of similar races across the region, justifying the need for concerted efforts by countries in southern Africa to develop and deploy more efficient and sustainable strategies to manage the disease. Understanding the genetic variability and composition of Pt is a pre-requisite for cultivar release with appropriate resistance gene combinations for sustainable disease management. This review highlights the variability and distribution of the Pt population, and the current control strategies, challenges and future prospects of breeding wheat varieties with durable leaf rust resistance in southern Africa. The importance of regular, collaborative and efficient surveillance of the pathogen and germplasm development across southern Africa is discussed, coupled with the potential of using modern breeding technologies to produce wheat cultivars with durable resistance.

Molecular characterization and expression analysis of pitaya (Hylocereus polyrhizus) HpLRR genes in response to Neoscytalidium dimidiatum infection

BACKGROUND

Canker disease caused by Neoscytalidium dimidiatum is a devastating disease resulting in a major loss to the pitaya industry. However, resistance proteins in plants play crucial roles to against pathogen infection. Among resistance proteins, the leucine-rich repeat (LRR) protein is a major family that plays crucial roles in plant growth, development, and biotic and abiotic stress responses, especially in disease defense.

RESULTS

In the present study, a transcriptomics analysis identified a total of 272 LRR genes, 233 of which had coding sequences (CDSs), in the plant pitaya (Hylocereus polyrhizus) in response to fungal Neoscytalidium dimidiatum infection. These genes were divided into various subgroups based on specific domains and phylogenetic analysis. Molecular characterization, functional annotation of proteins, and an expression analysis of the LRR genes were conducted. Additionally, four LRR genes (CL445.Contig4_All, Unigene28_All, CL28.Contig2_All, and Unigene2712_All, which were selected because they had the four longest CDSs were further assessed using quantitative reverse transcription PCR (qRT-PCR) at different fungal infection stages in different pitaya species (Hylocereus polyrhizus and Hylocereus undatus), in different pitaya tissues, and after treatment with salicylic acid (SA), methyl jasmonate (MeJA), and abscisic acid (ABA) hormones. The associated protein functions and roles in signaling pathways were identified.

CONCLUSIONS

This study provides a comprehensive overview of the HpLRR family genes at transcriptional level in pitaya in response to N. dimidiatum infection, it will be helpful to understand the molecular mechanism of pitaya canker disease, and lay a strong foundation for further research.

Gene expression profiling identifies pathways involved in seed maturation of Jatropha curcas

Background

Jatropha curcas, a tropical shrub, is a promising biofuel crop, which produces seeds with high content of oil and protein. To better understand the maturation process of J. curcas seeds and to improve its agronomic performance, a two-step approach was performed in six different maturation stages of seeds: 1) generation of the entire transcriptome of J. curcas seeds using 454-Roche sequencing of a cDNA library, 2) comparison of transcriptional expression levels using a custom Agilent 8x60K oligonucleotide microarray.

Results

A total of 793,875 high-quality reads were assembled into 19,382 unique full-length contigs, of which 13,507 could be annotated with Gene Ontology (GO) terms. Microarray data analysis identified 9111 probes (out of 57,842 probes), which were differentially expressed between the six maturation stages. The expression results were validated for 75 selected transcripts based on expression levels, predicted function, pathway, and length.

Result from cluster analyses showed that transcripts associated with fatty acid, flavonoid, and phenylpropanoid biosynthesis were over-represented in the early stages, while those of lipid storage were over-represented in the late stages. Expression analyses of different maturation stages of J. curcas seed showed that most changes in transcript abundance occurred between the two last stages, suggesting that the timing of metabolic pathways during seed maturation in J. curcas occurs in late stages. The co-expression results showed that the hubs (CB5-D, CDR1, TT8, DFR, HVA22) with the highest number of edges, associated with fatty acid and flavonoid biosynthesis, are showing a decrease in their expression during seed maturation. Furthermore, seed development and hormone pathways are significantly well connected.

Conclusion

The obtained results revealed differentially expressed sequences (DESs) regulating important pathways related to seed maturation, which could contribute to the understanding of the complex regulatory network during seed maturation with the focus on lipid, flavonoid and phenylpropanoid biosynthesis. This study provides detailed information on transcriptional changes during J. curcas seed maturation and provides a starting point for a genomic survey of seed quality traits. The results highlighted specific genes and processes relevant to the molecular mechanisms involved in Jatropha seed maturation. These data can also be utilized regarding other Euphorbiaceae species.

Characterization of disease resistance genes in the Brassica napus pangenome reveals significant structural variation

Methods based on single nucleotide polymorphism (SNP), copy number variation (CNV) and presence/absence variation (PAV) discovery provide a valuable resource to study gene structure and evolution. However, as a result of these structural variations, a single reference genome is unable to cover the entire gene content of a species. Therefore, pangenomics analysis is needed to ensure that the genomic diversity within a species is fully represented. Brassica napus is one of the most important oilseed crops in the world and exhibits variability in its resistance genes across different cultivars. Here, we characterized resistance gene distribution across 50 B. napus lines. We identified a total of 1749 resistance gene analogs (RGAs), of which 996 are core and 753 are variable, 368 of which are not present in the reference genome (cv. Darmor-bzh). In addition, a total of 15 318 SNPs were predicted within 1030 of the RGAs. The results showed that core R-genes harbour more SNPs than variable genes. More nucleotide binding site-leucine-rich repeat (NBS-LRR) genes were located in clusters than as singletons, with variable genes more likely to be found in clusters. We identified 106 RGA candidates linked to blackleg resistance quantitative trait locus (QTL). This study provides a better understanding of resistance genes to target for genomics-based improvement and improved disease resistance.

NLGenomeSweeper: A Tool for Genome-Wide NBS-LRR Resistance Gene Identification

Although there are a number of bioinformatic tools to identify plant nucleotide-binding leucine-rich repeat (NLR) disease resistance genes based on conserved protein sequences, only a few of these tools have attempted to identify disease resistance genes that have not been annotated in the genome. The overall goal of the NLGenomeSweeper pipeline is to annotate NLR disease resistance genes, including RPW8, in the genome assembly with high specificity and a focus on complete functional genes. This is based on the identification of the complete NB-ARC domain, the most conserved domain of NLR genes, using the BLAST suite. In this way, the tool has a high specificity for complete genes and relatively intact pseudogenes. The tool returns all candidate NLR gene locations as well as InterProScan ORF and domain annotations for manual curation of the gene structure.

Analysis of the Rdr1 gene family in different Rosaceae genomes reveals an origin of an R-gene cluster after the split of Rubeae within the Rosoideae subfamily

The Rdr1 gene confers resistance to black spot in roses and belongs to a large TNL gene family, which is organized in two major clusters at the distal end of chromosome 1. We used the recently available chromosome scale assemblies for the R. chinensis 'Old Blush' genome, re-sequencing data for nine rose species and genome data for Fragaria, Rubus, Malus and Prunus to identify Rdr1 homologs from different taxa within Rosaceae. Members of the Rdr1 gene family are organized into two major clusters in R. chinensis and at a syntenic location in the Fragaria genome. Phylogenetic analysis indicates that the two clusters existed prior to the split of Rosa and Fragaria and that one cluster has a more recent origin than the other. Genes belonging to cluster 2, such as the functional Rdr1 gene muRdr1A, were subject to a faster evolution than genes from cluster 1. As no Rdr1 homologs were found in syntenic positions for Prunus persica, Malus x domestica and Rubus occidentalis, a translocation of the Rdr1 clusters to the current positions probably happened after the Rubeae split from other groups within the Rosoideae approximately 70-80 million years ago during the Cretaceous period.

Physical Mapping of Pm57, a Powdery Mildew Resistance Gene Derived from Aegilops searsii

Powdery mildew caused by Blumeria graminis f. sp. tritici (Bgt) is one of many severe diseases that threaten bread wheat (Triticum aestivum L.) yield and quality worldwide. The discovery and deployment of powdery mildew resistance genes (Pm) can prevent this disease epidemic in wheat. In a previous study, we transferred the powdery mildew resistance gene Pm57 from Aegilops searsii into common wheat and cytogenetically mapped the gene in a chromosome region with the fraction length (FL) 0.75-0.87, which represents 12% segment of the long arm of chromosome 2Ss#1. In this study, we performed RNA-seq using RNA extracted from leaf samples of three infected and mock-infected wheat-Ae. searsii 2Ss#1 introgression lines at 0, 12, 24, and 48 h after inoculation with Bgt isolates. Then we designed 79 molecular markers based on transcriptome sequences and physically mapped them to Ae. searsii chromosome 2Ss#1- in seven intervals. We used these markers to identify 46 wheat-Ae. searsii 2Ss#1 recombinants induced by ph1b, a deletion mutant of pairing homologous (Ph) genes. After analyzing the 46 ph1b-induced 2Ss#1L recombinants in the region where Pm57 is located with different Bgt-responses, we physically mapped Pm57 gene on the long arm of 2Ss#1 in a 5.13 Mb genomic region, which was flanked by markers X67593 (773.72 Mb) and X62492 (778.85 Mb). By comparative synteny analysis of the corresponding region on chromosome 2B in Chinese Spring (T. aestivum L.) with other model species, we identified ten genes that are putative plant defense-related (R) genes which includes six coiled-coil nucleotide-binding site-leucine-rich repeat (CNL), three nucleotide-binding site-leucine-rich repeat (NL) and a leucine-rich receptor-like repeat (RLP) encoding proteins. This study will lay a foundation for cloning of Pm57, and benefit the understanding of interactions between resistance genes of wheat and powdery mildew pathogens.

Genome-Wide Mining of Disease Resistance Gene Analogs Using Conserved Domains

The production of legume crop species is severely affected by disease, imposing a significant yield loss annually worldwide. Plant resistance gene analogs (RGAs) play specific roles in plant resistance responses, and their identification and subsequent application in breeding programs help to reduce this yield loss. RGAs contain conserved domains and motifs, which can be used for their identification and classification. Nucleotide-binding site-leucine-rich repeat (NLR), receptor like kinase (RLK), and receptor like protein (RLP) genes are the main types of RGAs. Computational identification and characterization of RGAs has been performed successfully among different plant species. Here, we explain the computational workflow for genome-wide RGA identification in legumes.

Deployment of Genetic and Genomic Tools Toward Gaining a Better Understanding of Rice-Xanthomonas oryzae pv. oryzae Interactions for Development of Durable Bacterial Blight Resistant Rice

Rice is the most important food crop worldwide and sustainable rice production is important for ensuring global food security. Biotic stresses limit rice production significantly and among them, bacterial blight (BB) disease caused by Xanthomonas oryzae pv. oryzae (Xoo) is very important. BB reduces rice yields severely in the highly productive irrigated and rainfed lowland ecosystems and in recent years; the disease is spreading fast to other rice growing ecosystems as well. Being a vascular pathogen, Xoo interferes with a range of physiological and biochemical exchange processes in rice. The response of rice to Xoo involves specific interactions between resistance (R) genes of rice and avirulence (Avr) genes of Xoo, covering most of the resistance genes except the recessive ones. The genetic basis of resistance to BB in rice has been studied intensively, and at least 44 genes conferring resistance to BB have been identified, and many resistant rice cultivars and hybrids have been developed and released worldwide. However, the existence and emergence of new virulent isolates of Xoo in the realm of a rapidly changing climate necessitates identification of novel broad-spectrum resistance genes and intensification of gene-deployment strategies. This review discusses about the origin and occurrence of BB in rice, interactions between Xoo and rice, the important roles of resistance genes in plant's defense response, the contribution of rice resistance genes toward development of disease resistance varieties, identification and characterization of novel, and broad-spectrum BB resistance genes from wild species of Oryza and also presents a perspective on potential strategies to achieve the goal of sustainable disease management.

Inheritance Pattern and Molecular Markers for Resistance to Blackleg Disease in Cabbage

The inheritance and causal loci for resistance to blackleg, a devastating disease of Brassicaceous crops, are yet to be known in cabbage (Brassica oleracea L.). Here, we report the pattern of inheritance and linked molecular marker for this trait. A segregating BC₁ population consisting of 253 plants was raised from resistant and susceptible parents, L29 (♀) and L16 (♂), respectively. Cotyledon resistance bioassay of BC₁ population, measured based on a scale of 0-9 at 12 days after inoculation with Leptosphaeria maculans isolate 03-02 s, revealed the segregation of resistance and ratio, indicative of dominant monogenic control of the trait. Investigation of potential polymorphism in the previously identified differentially expressed genes within the collinear region of 'B. napus blackleg resistant loci Rlm1' in B. oleracea identified two insertion/deletion (InDel) mutations in the intron and numerous single nucleotide polymorphisms (SNPs) throughout the LRR-RLK gene Bol040029, of which six SNPs in the first exon caused the loss of two LRR domains in the susceptible line. An InDel marker, BLR-C-InDel based on the InDel mutations, and a high resolution melting (HRM) marker, BLR-C-2808 based on the SNP C²⁸⁰⁸T in the second exon were developed, which predicated the resistance status of the BC₁ population with 80.24%, and of 24 commercial inbred lines with 100% detection accuracy. This is the first report of inheritance and molecular markers linked with blackleg resistance in cabbage. This study will enhance our understanding of the trait, and will be helpful in marker assisted breeding aiming at developing resistant cabbage varieties.

Plant virus interaction mechanism and associated pathways in mosaic disease of small cardamom (Elettaria cardamomum Maton) by RNA-Seq approach

Small cardamom (Elettaria cardamomum), grown in limited coastal tropical countries is one of the costliest and widely exported agri-produce having global turnover of >10 billion USD. Mosaic/marble disease is one of the major impediments that requires understanding of disease at molecular level. Neither whole genome sequence nor any genomic resources are available, thus RNA seq approach can be a rapid and economical alternative. De novo transcriptome assembly was done with Illumina Hiseq data. A total of 5317 DEGs, 2267 TFs, 114 pathways and 175,952 genic region putative markers were obtained. Gene regulatory network analysis deciphered molecular events involved in marble disease. This is the first transcriptomic report revealing disease mechanism mediated by perturbation in auxin homeostasis and ethylene signalling leading to senescence. The web-genomic resource (SCMVTDb) catalogues putative molecular markers, candidate genes and transcript information. SCMVTDb can be used in germplasm improvement against mosaic disease in endeavour of small cardamom productivity. Availability of genomic resource, SCMVTDb: http://webtom.cabgrid.res.in/scmvtdb/.

Genome survey of resistance gene analogs in sugarcane: genomic features and differential expression of the innate immune system from a smut-resistant genotype

BACKGROUND

Resistance genes composing the two-layer immune system of plants are thought as important markers for breeding pathogen-resistant crops. Many have been the attempts to establish relationships between the genomic content of Resistance Gene Analogs (RGAs) of modern sugarcane cultivars to its degrees of resistance to diseases such as smut. However, due to the highly polyploid and heterozygous nature of sugarcane genome, large scale RGA predictions is challenging.

RESULTS

We predicted, searched for orthologs, and investigated the genomic features of RGAs within a recently released sugarcane elite cultivar genome, alongside the genomes of sorghum, one sugarcane ancestor (Saccharum spontaneum), and a collection of de novo transcripts generated for six modern cultivars. In addition, transcriptomes from two sugarcane genotypes were obtained to investigate the roles of RGAs differentially expressed (RGADE) in their distinct degrees of resistance to smut. Sugarcane references lack RGAs from the TNL class (Toll-Interleukin receptor (TIR) domain associated to nucleotide-binding site (NBS) and leucine-rich repeat (LRR) domains) and harbor elevated content of membrane-associated RGAs. Up to 39% of RGAs were organized in clusters, and 40% of those clusters shared synteny. Basically, 79% of predicted NBS-encoding genes are located in a few chromosomes. S. spontaneum chromosome 5 harbors most RGADE orthologs responsive to smut in modern sugarcane. Resistant sugarcane had an increased number of RGAs differentially expressed from both classes of RLK (receptor-like kinase) and RLP (receptor-like protein) as compared to the smut-susceptible. Tandem duplications have largely contributed to the expansion of both RGA clusters and the predicted clades of RGADEs.

CONCLUSIONS

Most of smut-responsive RGAs in modern sugarcane were potentially originated in chromosome 5 of the ancestral S. spontaneum genotype. Smut resistant and susceptible genotypes of sugarcane have a distinct pattern of RGADE. TM-LRR (transmembrane domains followed by LRR) family was the most responsive to the early moment of pathogen infection in the resistant genotype, suggesting the relevance of an innate immune system. This work can help to outline strategies for further understanding of allele and paralog expression of RGAs in sugarcane, and the results should help to develop a more applied procedure for the selection of resistant plants in sugarcane.

Whole-Genome and Transposed Duplication Contributes to the Expansion and Diversification of TLC Genes in Maize

TLC (TRAM/LAG/CRN8) proteins play important roles in ceramide metabolism and mycotoxin resistance. Herein a comparative genomics analysis of TLCs was performed in 31 plant and 3 species from other kingdoms, with an emphasis mainly on maize. TLCs were conserved across kingdoms and expanded in angiosperms, largely due to whole-genome/segmental duplication (WGD/SD) under purifying selection. Phylogeny reconstruction by maximum-likelihood method uncovered five TLC clades, subsequently named as TRAM/LAG, CLN8, PS-TLC, TM136 and TLCD clades. Each clade of TLCs shared specific transmembrane regions and motif composition. Divisions of conserved motifs to subunits may have occurred in TM136-type TLCs. Focusing on maize, five WGD and two DNA-mediated transposed duplication (TD) pairs were discovered, accounting for 61.11% ZmTLCs. Combined with further expression analysis, significant divergence was found in expression patterns between most maize WGD pairs, indicating subfunctionalization or/and neofunctionalization. Moreover, ZmTLC5, a deduced parental copy in a TD pair, was highly induced under FB₁ and fungus pathogen injection and exhibited potential capacity to respond to environmental stimuli. Additionally, population genetics analysis showed that ZmTLC10 in the CLN8-clade may have experienced significant positive selection and differentiated between wild and inbred maize populations. Overall, our results help to decipher the evolutionary history of TLCs in maize and plants, facilitating further functional analysis of them.

Molecular dissection of resistance gene cluster and candidate gene identification of Pl17 and Pl19 in sunflower by whole-genome resequencing

Sunflower (Helianthus annuus L.) production is challenged by different biotic and abiotic stresses, among which downy mildew (DM) is a severe biotic stress that is detrimental to sunflower yield and quality in many sunflower-growing regions worldwide. Resistance against its infestation in sunflower is commonly regulated by single dominant genes. Pl₁₇ and Pl₁₉ are two broad-spectrum DM resistance genes that have been previously mapped to a gene cluster spanning a 3.2 Mb region at the upper end of sunflower chromosome 4. Using a whole-genome resequencing approach combined with a reference sequence-based chromosome walking strategy and high-density mapping populations, we narrowed down Pl₁₇ to a 15-kb region flanked by SNP markers C4_5711524 and SPB0001. A prospective candidate gene HanXRQChr04g0095641 for Pl₁₇ was identified, encoding a typical TNL resistance gene protein. Pl₁₉ was delimited to a 35-kb region and was approximately 1 Mb away from Pl₁₇, flanked by SNP markers C4_6676629 and C4_6711381. The only gene present within the delineated Pl₁₉ locus in the reference genome, HanXRQChr04g0095951, was predicted to encode an RNA methyltransferase family protein. Six and eight SNP markers diagnostic for Pl₁₇ and Pl₁₉, respectively, were identified upon evaluation of 96 diverse sunflower lines, providing a very useful tool for marker-assisted selection in sunflower breeding programs.

Transcriptome analysis of a nematode resistant and susceptible upland cotton line at two critical stages of Meloidogyne incognita infection and development

Host plant resistance is the most practical approach to control the Southern root-knot nematode (Meloidogyne incognita; RKN), which has emerged as one of the most serious economic pests of Upland cotton (Gossypium hirsutum L.). Previous QTL analyses have identified a resistance locus on chromosome 11 (qMi-C11) affecting galling and another locus on chromosome-14 (qMi-C14) affecting egg production. Although these two QTL regions were fine mapped and candidate genes identified, expression profiling of genes would assist in further narrowing the list of candidate genes in the QTL regions. We applied the comparative transcriptomic approach to compare expression profiles of genes between RKN susceptible and resistance genotypes at an early stage of RKN development that coincides with the establishment of a feeding site and at the late stage of RKN development that coincides with RKN egg production. Sequencing of cDNA libraries produced over 315 million reads of which 240 million reads (76%) were mapped on to the Gossypium hirsutum genome. A total of 3,789 differentially expressed genes (DEGs) were identified which were further grouped into four clusters based on their expression profiles. A large number of DEGs were found to be down regulated in the susceptible genotype at the late stage of RKN development whereas several genes were up regulated in the resistant genotype. Key enriched categories included transcription factor activity, defense response, response to phyto-hormones, cell wall organization, and protein serine/threonine kinase activity. Our results also show that the DEGs in the resistant genotype at qMi-C11 and qMi-C14 loci displayed higher expression of defense response, detoxification and callose deposition genes, than the DEGs in the susceptible genotype.

A Novel QTL for Root-Knot Nematode Resistance is Identified from a South African Sweet Sorghum Line

Southern root-knot nematodes, Meloidogyne incognita, feed on the underground portions of hundreds of plant species and affect nutrient partitioning and water uptake of the host plants. Sorghum (Sorghum bicolor) is often not significantly damaged by southern root-knot nematodes (RKN) but some sorghum genotypes support greater population densities of RKN than other genotypes. These higher nematode populations increase the risk of damage to subsequently planted susceptible crops. A previous study identified a major quantitative trait locus (QTL) for RKN resistance on sorghum chromosome (chr.) 3. To maintain long-term resistance, multiple resistance genes should be pyramided in a cultivar. In this study, we identified a new source of RKN resistance, created a mapping population, and identified single-nucleotide polymorphism markers using genotyping-by-sequencing of the segregating population. Use of single-marker analysis and composite interval mapping identified a single QTL on chr. 5 that was associated with egg number and egg number per gram of root from the resistant sweet sorghum line PI 144134. This region on chr. 5 and the prior QTL on chr. 3 can be potentially moved from PI 144134 and Honey Drip, respectively, into elite sorghum germplasm via marker-assisted selection for more durable resistance.

Determination of the Key Resistance Gene Analogs Involved in Ascochyta rabiei Recognition in Chickpea

Chickpea (Cicer arietinum L.) is an important cool season food legume, however, its production is severely constrained by the foliar disease Ascochyta blight caused by the fungus Ascochyta rabiei (syn. Phoma rabiei). Several disease management options have been developed to control the pathogen, including breeding for host plant resistance. However, the pathogen population is evolving to produce more aggressive isolates. For host resistance to be effective, the plant must quickly recognize the pathogen and instigate initial defense mechanisms, optimally at the point of contact. Given that the most resistant host genotypes display rapid pathogen recognition and response, the approach taken was to assess the type, speed and pattern of recognition via Resistance Gene Analog (RGA) transcription among resistant and susceptible cultivated chickpea varieties. RGAs are key factors in the recognition of plant pathogens and the signaling of inducible defenses. In this study, a suite of RGA loci were chosen for further investigation from both published literature and from newly mined homologous sequences within the National Center for Biotechnology Information (NCBI) database. Following their validation in the chickpea genome, 10 target RGAs were selected for differential expression analysis in response to A. rabiei infection. This was performed in a set of four chickpea varieties including two resistant cultivars (ICC3996 and PBA Seamer), one moderately resistant cultivar (PBA HatTrick) and one susceptible cultivar (Kyabra). Gene expression at each RGA locus was assessed via qPCR at 2, 6, and 24 h after A. rabiei inoculation with a previously characterized highly aggressive isolate. As a result, all loci were differentially transcribed in response to pathogen infection in at least one genotype and at least one time point after inoculation. Among these, the differential expression of four RGAs was significant and consistently increased in the most resistant genotype ICC3996 immediately following inoculation, when spore germination began and ahead of penetration into the plant's epidermal tissues. Further in silico analyses indicated that the differentially transcribed RGAs function through ADP-binding within the pathogen recognition pathway. These represent clear targets for future functional validation and potential for selective resistance breeding for introgression into elite cultivars.

Introgression of resistance to Rhopalosiphum padi L. from wild barley into cultivated barley facilitated by doubled haploid and molecular marker techniques

Long-term pre-breeding using Hordeum vulgare ssp. spontaneum as a donor of bird cherry-oat aphid resistance has resulted in agronomically improved resistance sources of barley along with easy-to-use molecular markers. Bird cherry-oat aphid (Rhopalosiphum padi L.) is a pest and a virus vector in barley to which there are no bred-resistant cultivars. The present study describes how resistance from Hordeum vulgare ssp. spontaneum has been introgressed in cultivated barley via five successive crosses with the same cultivar Lina (BC) and in parallel with other more modern barley cultivars. Most of the selections for resistance are based on measurements of individual aphid growth in the laboratory. This very slow phenotyping method has been complemented by molecular marker evaluation and application in part of the breeding material. Doubled haploid production in each generation has been crucial for more precise selection of lines with the quantitatively expressed resistance. A field trial of selected "BC₃"-generation lines essentially confirmed the laboratory results, so did genotyping of the whole pedigree of parents and selected "BC₂" and "BC₄" offspring lines. The Infinium iSelect 50 K SNP assay confirmed relationships between lines and discerned several new markers for a resistance QTL on chromosome 2H.

GWAS analysis of sorghum association panel lines identifies SNPs associated with disease response to Texas isolates of Colletotrichum sublineola

SNPs identify prospective genes related to response to Colletotrichum sublineola (anthracnose) in the sorghum association panel lines. Sorghum association panel (SAP) lines were scored over several years for response to Colletotrichum sublineola, the causal agent of the disease anthracnose. Known resistant and susceptible lines were included each year to verify successful inoculation. Over 79,000 single-nucleotide polymorphic (SNP) loci from a publicly available genotype by sequencing dataset available for the SAP lines were used with TASSEL association mapping software to identify chromosomal locations associated with differences in disease response. When the top-scoring SNPs were mapped to the published sorghum genome, in each case, the nearest annotated gene has precedence for a role in host defense.

Variation in abundance of predicted resistance genes in the Brassica oleracea pangenome

Brassica oleracea is an important agricultural species encompassing many vegetable crops including cabbage, cauliflower, broccoli and kale; however, it can be susceptible to a variety of fungal diseases such as clubroot, blackleg, leaf spot and downy mildew. Resistance to these diseases is meditated by specific disease resistance genes analogs (RGAs) which are differently distributed across B. oleracea lines. The sequenced reference cultivar does not contain all B. oleracea genes due to gene presence/absence variation between individuals, which makes it necessary to search for RGA candidates in the B. oleracea pangenome. Here we present a comparative analysis of RGA candidates in the pangenome of B. oleracea. We show that the presence of RGA candidates differs between lines and suggests that in B. oleracea, SNPs and presence/absence variation drive RGA diversity using separate mechanisms. We identified 59 RGA candidates linked to Sclerotinia, clubroot, and Fusarium wilt resistance QTL, and these findings have implications for crop breeding in B. oleracea, which may also be applicable in other crops species.

Pooled DNA sequencing to identify SNPs associated with a major QTL for bacterial wilt resistance in Italian ryegrass (Lolium multiflorum Lam.)

SNPs and candidate genes associated with bacterial wilt resistance in Italian ryegrass were identified by sequencing the parental plants and pooled F₁ progeny of a segregating population. Italian ryegrass (Lolium multiflorum Lam.) is one of the most important forage grass species in temperate regions. Its yield, quality and persistency can significantly be reduced by bacterial wilt, a serious disease caused by Xanthomonas translucens pv. graminis. Although a major QTL for bacterial wilt resistance has previously been reported, detailed knowledge on underlying genes and DNA markers to allow for efficient resistance breeding strategies is currently not available. We used pooled DNA sequencing to characterize a major QTL for bacterial wilt resistance of Italian ryegrass and to develop inexpensive sequence-based markers to efficiently target resistance alleles for marker-assisted recurrent selection. From the mapping population segregating for the QTL, DNA of 44 of the most resistant and 44 of the most susceptible F₁ individuals was pooled and sequenced using the Illumina HiSeq 2000 platform. Allele frequencies of 18 × 10⁶ single nucleotide polymorphisms (SNP) were determined in the resistant and susceptible pool. A total of 271 SNPs on 140 scaffold sequences of the reference parental genome showed significantly different allele frequencies in both pools. We converted 44 selected SNPs to KASP™ markers, genetically mapped these proximal to the major QTL and thus validated their association with bacterial wilt resistance. This study highlights the power of pooled DNA sequencing to efficiently target binary traits in biparental mapping populations. It delivers genome sequence data, SNP markers and potential candidate genes which will allow to implement marker-assisted strategies to fix bacterial wilt resistance in outcrossing breeding populations of Italian ryegrass.

Divergence between bread wheat and Triticum militinae in the powdery mildew resistance QPm.tut-4A locus and its implications for cloning of the resistance gene

A segment of Triticum militinae chromosome 7G harbors a gene(s) conferring powdery mildew resistance which is effective at both the seedling and the adult plant stages when transferred into bread wheat (T. aestivum). The introgressed segment replaces a piece of wheat chromosome arm 4AL. An analysis of segregating materials generated to positionally clone the gene highlighted that in a plant heterozygous for the introgression segment, only limited recombination occurs between the introgressed region and bread wheat 4A. Nevertheless, 75 genetic markers were successfully placed within the region, thereby confining the gene to a 0.012 cM window along the 4AL arm. In a background lacking the Ph1 locus, the localized rate of recombination was raised 33-fold, enabling the reduction in the length of the region containing the resistance gene to a 480 kbp stretch harboring 12 predicted genes. The substituted segment in the reference sequence of bread wheat cv. Chinese Spring is longer (640 kbp) and harbors 16 genes. A comparison of the segments' sequences revealed a high degree of divergence with respect to both their gene content and nucleotide sequence. Of the 12 T. militinae genes, only four have a homolog in cv. Chinese Spring. Possible candidate genes for the resistance have been identified based on function predicted from their sequence.

Flowering plant immune repertoires expand under mycorrhizal symbiosis

Immune perception in flowering plants is mediated by a repertoire of cytoplasmic and cell-surface receptors that detect invading microbes and their effects on cells. Here, we show that several large families of immune receptors exhibit size variations related to a plant's competence to host symbiotic root fungi (mycorrhiza). Plants that do not participate in mycorrhizal associations have significantly smaller immune repertoires, while the most promiscuous symbiotic hosts (ectomycorrhizal plant species) have significantly larger immune repertoires. By contrast, we find no significant increase in immune repertoire size among legumes competent to form a symbiosis with nitrogen-fixing bacteria (rhizobia). To explain these observations, we hypothesize that plant immune repertoire size expands with symbiote species diversity.

Genome-Wide Association Studies for Pasmo Resistance in Flax (Linum usitatissimum L.)

Pasmo is one of the most widespread diseases threatening flax production. To identify genetic regions associated with pasmo resistance (PR), a genome-wide association study was performed on 370 accessions from the flax core collection. Evaluation of pasmo severity was performed in the field from 2012 to 2016 in Morden, MB, Canada. Genotyping-by-sequencing has identified 258,873 single nucleotide polymorphisms (SNPs) distributed on all 15 flax chromosomes. Marker-trait associations were identified using ten different statistical models. A total of 692 unique quantitative trait nucleotides (QTNs) associated with 500 putative quantitative trait loci (QTL) were detected from six phenotypic PR datasets (five individual years and average across years). Different QTNs were identified with various statistical models and from individual PR datasets, indicative of the complementation between analytical methods and/or genotype × environment interactions of the QTL effects. The single-locus models tended to identify large-effect QTNs while the multi-loci models were able to detect QTNs with smaller effects. Among the putative QTL, 67 had large effects (3-23%), were stable across all datasets and explained 32-64% of the total variation for PR in the various datasets. Forty-five of these QTL spanned 85 resistance gene analogs including a large toll interleukin receptor, nucleotide-binding site, leucine-rich repeat (TNL) type gene cluster on chromosome 8. The number of QTL with positive-effect or favorite alleles (NPQTL) in accessions was significantly correlated with PR (R 2 = 0.55), suggesting that these QTL effects are mainly additive. NPQTL was also significantly associated with morphotype (R 2 = 0.52) and major QTL with positive effect alleles were present in the fiber type accessions. The 67 large effect QTL are suited for marker-assisted selection and the 500 QTL for effective genomic prediction in PR molecular breeding.

Development of Clustered Resistance Gene Analogs-Based Markers of Resistance to Phytophthora capsici in Chili Pepper

The soil-borne pathogen Phytophthora capsici causes severe destruction of Capsicum spp. Resistance in Capsicum against P. capsici is controlled by numerous minor quantitative trait loci (QTLs) and a consistent major QTL on chromosome 5. Molecular markers on Capsicum chromosome 5 have been developed to identify the predominant genetic contributor to resistance but have achieved little success. In this study, previously reported molecular markers were used to reanalyze the major QTL region on chromosome 5 (6.2 Mbp to 139.2 Mbp). Candidate resistance gene analogs (RGAs) were identified in the extended major QTL region including 14 nucleotide binding site leucine-rich repeats, 3 receptor-like kinases, and 1 receptor-like protein. Sequence comparison of the candidate RGAs was performed between two Capsicum germplasms that are resistant and susceptible, respectively, to P. capsici. 11 novel RGA-based markers were developed through high-resolution melting analysis which were closely linked to the major QTL for P. capsici resistance. Among the markers, CaNB-5480 showed the highest cosegregation rate at 86.9% and can be applied to genotyping of the germplasms that were not amenable by previous markers. With combination of three markers such as CaNB-5480, CaRP-5130 and CaNB-5330 increased genotyping accuracy for 61 Capsicum accessions. These could be useful to facilitate high-throughput germplasm screening and further characterize resistance genes against P. capsici in pepper.