Identification of a Differentially Expressed TIR-NBS-LRR Gene in a Major QTL Associated to Leaf Rust Resistance in Salix

Partially Resistant Avocado Rootstock Dusa® Shows Prolonged Upregulation of Nucleotide Binding-Leucine Rich Repeat Genes in Response to Phytophthora cinnamomi Infection

Avocado is an important agricultural food crop in many countries worldwide. Phytophthora cinnamomi, a hemibiotrophic oomycete, remains one of the most devastating pathogens within the avocado industry, as it is near impossible to eradicate from areas where the pathogen is present. A key aspect to Phytophthora root rot disease management is the use of avocado rootstocks partially resistant to P. cinnamomi, which demonstrates an increased immune response following infection. In plant species, Nucleotide binding-Leucine rich repeat (NLR) proteins form an integral part of pathogen recognition and Effector triggered immune responses (ETI). To date, a comprehensive set of Persea americana NLR genes have yet to be identified, though their discovery is crucial to understanding the molecular mechanisms underlying P. americana-P. cinnamomi interactions. In this study, a total of 161 PaNLR genes were identified in the P. americana West-Indian pure accession genome. These putative resistance genes were characterized using bioinformatic approaches and grouped into 13 distinct PaNLR gene clusters, with phylogenetic analysis revealing high sequence similarity within these clusters. Additionally, PaNLR expression levels were analyzed in both a partially resistant (Dusa®) and a susceptible (R0.12) avocado rootstock infected with P. cinnamomi using an RNA-sequencing approach. The results showed that the partially resistant rootstock has increased expression levels of 84 PaNLRs observed up to 24 h post-inoculation, while the susceptible rootstock only showed increased PaNLR expression during the first 6 h post-inoculation. Results of this study may indicate that the partially resistant avocado rootstock has a stronger, more prolonged ETI response which enables it to suppress P. cinnamomi growth and combat disease caused by this pathogen. Furthermore, the identification of PaNLRs may be used to develop resistant rootstock selection tools, which can be employed in the avocado industry to accelerate rootstock screening programs.

Comparative transcriptomics and eQTL mapping of response to Melampsora americana in selected Salix purpurea F2 progeny

Background

Melampsora spp. rusts are the greatest pathogen threat to shrub willow (Salix spp.) bioenergy crops. Genetic resistance is key to limit the effects of these foliar diseases on host response and biomass yield, however, the genetic basis of host resistance has not been characterized. The addition of new genomic resources for Salix provides greater power to investigate the interaction between S. purpurea and M. americana, species commonly found in the Northeast US. Here, we utilize 3′ RNA-seq to investigate host-pathogen interactions following controlled inoculations of M. americana on resistant and susceptible F₂S. purpurea genotypes identified in a recent QTL mapping study. Differential gene expression, network analysis, and eQTL mapping were used to contrast the response to inoculation and to identify associated candidate genes.

Results

Controlled inoculation in a replicated greenhouse study identified 19 and 105 differentially expressed genes between resistant and susceptible genotypes at 42 and 66 HPI, respectively. Defense response gene networks were activated in both resistant and susceptible genotypes and enriched for many of the same defense response genes, yet the hub genes of these common response modules showed greater mean expression among the resistant plants. Further, eight and six eQTL hotspots were identified at 42 and 66 HPI, respectively. The combined results of three analyses highlight 124 candidate genes in the host for further analysis while analysis of pathogen RNA showed differential expression of 22 genes, two of which are candidate pathogen effectors.

Conclusions

We identified two differentially expressed M. americana transcripts and 124 S. purpurea genes that are good candidates for future studies to confirm their role in conferring resistance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08254-1.

Differential Susceptibility of Diverse Salix spp. to Melampsora americana and Melampsora paradoxa

Melampsora spp. willow rust is the most serious disease of shrub willow bioenergy production in the northeastern United States. Recent phylogenetic studies have identified several Melampsora spp. present on willow in the Northeast; however, in-depth understanding of Melampsora spp. host susceptibility remain unresolved. In this study, a panel of 82 rust isolates collected from the northeastern United States were genotyped via ribosomal DNA sequencing and a subset of these isolates were assayed for host susceptibility. This work revealed that Melampsora americana is the most prevalent species in the sampled geographic region and that there is potential for rust resistance breeding using the Salix spp. taxa assayed. Additionally, leaf morphology traits of these Salix spp. hosts were quantified for correlation analysis, revealing that trichome density and stomata density are possible contributors to resistance. This work provides foundational rust pathology information, which is crucial for M. americana resistance breeding.

Identifying Plasmopara viticola resistance Loci in grapevine (Vitis amurensis) via genotyping-by-sequencing-based QTL mapping

Downy mildew, caused by Plasmopara viticola, is a major disease that affects grapevines, and a few resistance (R) genes have been identified thus far. In order to identify R genes, we investigated F₁ progeny from a cross between the downy mildew-resistant Vitis amurensis 'Shuang Hong' and the susceptible Vitis vinifera 'Cabernet Sauvignon'. The P. viticola-resistance of the progeny varied continuously and was segregated as a quantitative trait. Genotyping-by-sequencing was used to construct linkage maps. The integrated map spanned 1898.09 cM and included 5603 single nucleotide polymorphisms on 19 linkage groups (LGs). Linkage analysis identified three quantitative trait loci (QTLs) for P. viticola resistance: 22 (Rpv22) on LG 02, Rpv23 on LG15, and Rpv24 on LG18. The phenotypic variance contributed by these three QTLs ranged from 15.9 to 30.0%. qRT-PCR analysis of R-gene expression in these QTLs revealed a CC-NBS-LRR disease resistance gene RPP8, two LRR receptor-like serine/threonine-protein kinases, a serine/threonine-protein kinase BLUS1, a glutathione peroxidase 8, an ethylene-responsive transcription factor ERF038, and a transcription factor bZIP11 were induced by P. viticola, and these genes may play important role in P. viticola response.

Association genetics identifies a specifically regulated Norway spruce laccase gene, PaLAC5, linked to Heterobasidion parviporum resistance

It is important to improve the understanding of the interactions between the trees and pathogens and integrate this knowledge about disease resistance into tree breeding programs. The conifer Norway spruce (Picea abies) is an important species for the forest industry in Europe. Its major pathogen is Heterobasidion parviporum, causing stem and root rot. In this study, we identified 11 Norway spruce QTLs (Quantitative trait loci) that correlate with variation in resistance to H. parviporum in a population of 466 trees by association genetics. Individual QTLs explained between 2.1 and 5.2% of the phenotypic variance. The expression of candidate genes associated with the QTLs was analysed in silico and in response to H. parviporum hypothesizing that (a) candidate genes linked to control of fungal sapwood growth are more commonly expressed in sapwood, and; (b) candidate genes associated with induced defences are respond to H. parviporum inoculation. The Norway spruce laccase PaLAC5 associated with control of lesion length development is likely to be involved in the induced defences. Expression analyses showed that PaLAC5 responds specifically and strongly in close proximity to the H. parviporum inoculation. Thus, PaLAC5 may be associated with the lignosuberized boundary zone formation in bark adjacent to the inoculation site.

Combining QTL Mapping with Genome Resequencing Identifies an Indel in an R Gene that is Associated with Variation in Leaf Rust Disease Resistance in Poplar

Poplar trees (Populus spp.) are important and are widely grown worldwide. However, the extensive occurrence of leaf rust disease caused by Melampsora spp. seriously inhibits their growth and reduces their biomass. In our previous study, a high-quality genetic map was constructed for the poplar F1 population I-69 × XYY by using next-generation sequencing-based genotyping-by-sequencing. Here, we collected phenotypic data on leaf rust disease resistance on three different dates for all 300 progenies of the F1 population. Combining a high-quality genetic map and phenotypic data, we were able to detect 11 major quantitative trait loci (QTLs) for leaf rust disease resistance. Among these 11 QTLs, two pairs were detected on at least two dates. In the corresponding genomic sequence, we found that resistance (R) gene clusters were located in these two QTL regions. By using genome resequencing, PCR confirmation and statistical analysis, a 611-bp deletion within an R gene in one QTL region was found to be associated with variation in leaf rust disease resistance. A PCR-based examination of this 611-bp deletion was performed. This 611-bp deletion was also found to affect mRNA splicing and form a new protein with the loss of some key protein domains. Based on this study, we were able to determine the genetic architecture of variation in poplar leaf rust disease resistance, and the 611-bp deletion in the R gene could be used as a diagnostic marker for future poplar molecular breeding.

Expression polymorphism at the ARPC4 locus links the actin cytoskeleton with quantitative disease resistance to Sclerotinia sclerotiorum in Arabidopsis thaliana

Quantitative disease resistance (QDR) is a form of plant immunity widespread in nature, and the only one active against broad host range fungal pathogens. The genetic determinants of QDR are complex and largely unknown, and are thought to rely partly on genes controlling plant morphology and development. We used genome-wide association mapping in Arabidopsis thaliana to identify ARPC4 as associated with QDR against the necrotrophic fungal pathogen Sclerotinia sclerotiorum. Mutants impaired in ARPC4 showed enhanced susceptibility to S. sclerotiorum, defects in the structure of the actin filaments and in their responsiveness to S. sclerotiorum. Disruption of ARPC4 also alters callose deposition and the expression of defense-related genes upon S. sclerotiorum infection. Analysis of ARPC4 diversity in A. thaliana identified one haplotype (ARPC4^R ) showing a c. 1 kbp insertion in ARPC4 regulatory region and associated with higher level of QDR. Accessions from the ARPC4^R haplotype showed enhanced ARPC4 expression upon S. sclerotiorum challenge, indicating that polymorphisms in ARPC4 regulatory region are associated with enhanced QDR. This work identifies a novel actor of plant QDR against a fungal pathogen and provides a prime example of genetic mechanisms leading to the recruitment of cell morphology processes in plant immunity.

High-density genetic map construction and quantitative trait loci identification for growth traits in (Taxodium distichum var. distichum × T. mucronatum) × T. mucronatum

BACKGROUND

'Zhongshanshan' is the general designation for the superior interspecific hybrid clones of Taxodium species, which is widely grown for economic and ecological purposes in southern China. Growth is the priority objective in 'Zhongshanshan' tree improvement. A high-density linkage map is vital to efficiently identify key quantitative trait loci (QTLs) that affect growth.

RESULTS

In total, 403.16 Gb of data, containing 2016,336 paired-end reads, was obtained after preprocessing. The average sequencing depth was 28.49 in T. distichum var. distichum, 25.18 in T. mucronatum, and 11.12 in each progeny. In total, 524,662 high-quality SLAFs were detected, of which 249,619 were polymorphic, and 6166 of the polymorphic markers met the requirements for use in constructing a genetic map. The final map harbored 6156 SLAF markers on 11 linkage groups, and was 1137.86 cM in length, with an average distance of 0.18 cM between adjacent markers. Separate QTL analyses of traits in different years by CIM detected 7 QTLs. While combining multiple-year data, 13 QTLs were detected by ICIM. 5 QTLs were repeatedly detected by the two methods, and among them, 3 significant QTLs (q6-2, q4-2 and q2-1) were detected in at least two traits. Bioinformatic analysis discoveried a gene annotated as a leucine-rich repeat receptor-like kinase gene within q4-2.

CONCLUSIONS

This map is the most saturated one constructed in a Taxodiaceae species to date, and would provide useful information for future comparative mapping, genome assembly, and marker-assisted selection.

Genetic Diversity and Genetic Relationships of Purple Willow (Salix purpurea L.) from Natural Locations

In this study, the genetic diversity and structure of 13 natural locations of Salix purpurea were determined with the use of AFLP (amplified length polymorphism), RAPD (randomly amplified polymorphic DNA) and ISSR (inter-simple sequence repeats). The genetic relationships between 91 examined S. purpurea genotypes were evaluated by analyses of molecular variance (AMOVA), principal coordinates analyses (PCoA) and UPGMA (unweighted pair group method with arithmetic mean) dendrograms for both single marker types and a combination of all marker systems. The locations were assigned to distinct regions and the analysis of AMOVA (analysis of molecular variance) revealed a high genetic diversity within locations. The genetic diversity between both regions and locations was relatively low, but typical for many woody plant species. The results noted for the analyzed marker types were generally comparable with few differences in the genetic relationships among S. purpurea locations. A combination of several marker systems could thus be ideally suited to understand genetic diversity patterns of the species. This study makes the first attempt to broaden our knowledge of the genetic parameters of the purple willow (S. purpurea) from natural location for research and several applications, inter alia breeding purposes.

Genetic Diversity and Genetic Relationships of Purple Willow (Salix purpurea L.) from Natural Locations

In this study, the genetic diversity and structure of 13 natural locations of Salix purpurea were determined with the use of AFLP (amplified length polymorphism), RAPD (randomly amplified polymorphic DNA) and ISSR (inter-simple sequence repeats). The genetic relationships between 91 examined S. purpurea genotypes were evaluated by analyses of molecular variance (AMOVA), principal coordinates analyses (PCoA) and UPGMA (unweighted pair group method with arithmetic mean) dendrograms for both single marker types and a combination of all marker systems. The locations were assigned to distinct regions and the analysis of AMOVA (analysis of molecular variance) revealed a high genetic diversity within locations. The genetic diversity between both regions and locations was relatively low, but typical for many woody plant species. The results noted for the analyzed marker types were generally comparable with few differences in the genetic relationships among S. purpurea locations. A combination of several marker systems could thus be ideally suited to understand genetic diversity patterns of the species. This study makes the first attempt to broaden our knowledge of the genetic parameters of the purple willow (S. purpurea) from natural location for research and several applications, inter alia breeding purposes.

Differential expression of NBS-LRR-encoding genes in the root transcriptomes of two Solanum phureja genotypes with contrasting resistance to Globodera rostochiensis

BACKGROUND

The characterization of major resistance genes (R genes) in the potato remains an important task for molecular breeding. However, R genes are rapidly evolving and frequently occur in genomes as clusters with complex structures, and their precise mapping and identification are complicated and time consuming.

RESULTS

Comparative analysis of root transcriptomes of Solanum phureja genotypes with contrasting resistance to Globodera rostochiensis revealed a number of differentially expressed genes. However, compiling a list of candidate R genes for further segregation analysis was hampered by their scarce annotation. Nevertheless, combination of transcriptomic analysis with data on predicted potato NBS-LRR-encoding genes considerably improved the quality of the results and provided a reasonable number of candidate genes that provide S. phureja with strong resistance to the potato golden cyst nematode.

CONCLUSION

Combination of comparative analyses of tissue-specific transcriptomes in resistant and susceptible genotypes may be used as an approach for the rapid identification of candidate potato R genes for co-segregation analysis and may be used in parallel with more sophisticated studies based on genome resequencing.