Genome-wide identification of NBS resistance genes in Populus trichocarpa

Phylogenetic, Structural, and Evolutionary Insights into Pepper NBS-LRR Resistance Genes

The comprehensive analysis of NBS-LRR resistance genes in the pepper (Capsicum annuum L.) genome reveals their structural diversity, evolutionary history, and functional importance in plant immunity. A total of 252 NBS-LRR genes were identified, distributed unevenly across all chromosomes, with 54% forming 47 gene clusters. These clusters, driven by tandem duplications and genomic rearrangements, underscore the dynamic evolution of resistance genes. Phylogenetic analysis demonstrated the dominance of the nTNL subfamily over the TNL subfamily, reflecting lineage-specific adaptations and evolutionary pressures. Structural analyses identified six conserved motifs (P-loop, RNBS-A, kinase-2, RNBS-B, RNBS-C, and GLPL) essential for ATP/GTP binding and resistance signaling. Subfamily-specific differences in motif composition and sequence similarity highlight their functional divergence and specialization. Comparative analyses across species further revealed a greater prevalence of nTNL genes in angiosperms, with significant losses of TNL genes in monocots. This study enhances our understanding of the evolution and diversification of plant-resistance genes and provides a foundation for developing disease-resistant crops through targeted breeding strategies.

Genome-Wide Identification and Stress Responses of Cowpea Thaumatin-like Proteins: A Comprehensive Analysis

Cowpea (Vigna unguiculata (L.) Walp.) is an important legume cultivated mainly in regions with limited water availability across the African and American continents. Its productivity is significantly affected by environmental stresses. Thaumatin-like proteins (TLPs), which belong to the PR-5 (pathogenesis-related 5) protein family, are known to be responsive to both biotic and abiotic stresses. However, their role remains controversial, with some TLPs associated with plant defense (particularly against fungal infections) and others associated with abiotic stresses response. In this study, we evaluated the structural diversity and gene expression of TLPs in cowpea (VuTLPs) under different stress conditions, including biotic [mechanical injury followed by inoculation with Cowpea Aphid-borne Mosaic Virus (CABMV) or Cowpea Severe Mosaic Virus (CPSMV)] and abiotic (root dehydration). Genomic anchoring of VuTLPs revealed 34 loci encoding these proteins. Neighbor- joining analysis clustered the VuTLPs into three distinct groups. We identified 15 segmental duplication and 6 tandem duplication gene pairs, with the majority of VuTLP genes found to be under purifying selection. Promoter analysis associated VuTLPs with bHLH, Dof-type, and MYB- related transcription factors, supporting their diverse roles. Diversity in VuTLP function was also observed in their expression profiles under the studied stress conditions. Gene expression data showed that most VuTLPs are recruited within the first minutes after biotic stress imposition. For the root dehydration assay, the most transcripts were up-regulated 150 min post-stress. Moreover, the gene expression data suggested that VuTLPs exhibit functional specialization depending on the stress condition, highlighting their diverse roles and biotechnological potential.

Pangenome Data Analysis Reveals Characteristics of Resistance Gene Analogs Associated with Sclerotinia sclerotiorum Resistance in Sunflower

The sunflower, an important oilseed crop and food source across the world, is susceptible to several pathogens, which cause severe losses in sunflower production. The utilization of genetic resistance is the most economical, effective measure to prevent infectious diseases. Based on the sunflower pangenome, in this study, we explored the variability of resistance gene analogs (RGAs) within the species. According to a comparative analysis of RGA candidates in the sunflower pangenome using the RGAugury pipeline, a total of 1344 RGAs were identified, comprising 1107 conserved, 199 varied, and 38 rare RGAs. We also identified RGAs associated with resistance against Sclerotinia sclerotiorum (S. sclerotiorum) in sunflower at the quantitative trait locus (QTL). A total of 61 RGAs were found to be located at four quantitative trait loci (QTLs). Through a detailed expression analysis of RGAs in one susceptible and two tolerant sunflower inbred lines (ILs) across various time points post inoculation, we discovered that 348 RGAs exhibited differential expression in response to Sclerotinia head rot (SHR), with 17 of these differentially expressed RGAs being situated within the QTL regions. In addition, 15 RGA candidates had gene introgression. Our data provide a better understanding of RGAs, which facilitate genomics-based improvements in disease resistance in sunflower.

Genome-wide characterization of nitric oxide-induced NBS-LRR genes from Arabidopsis thaliana and their association in monocots and dicots

BACKGROUND

Nitric oxide (NO) is pivotal in regulating the activity of NBS-LRR specific R genes, crucial components of the plant's immune system. It is noteworthy that previous research has not included a genome-wide analysis of NO-responsive NBS-LRR genes in plants.

RESULTS

The current study examined 29 NO-induced NBS-LRR genes from Arabidopsis thaliana, along with two monocots (rice and maize) and two dicots (soybean and tomato) using genome-wide analysis tools. These NBS-LRR genes were subjected to comprehensive characterization, including analysis of their physio-chemical properties, phylogenetic relationships, domain and motif identification, exon/intron structures, cis-elements, protein-protein interactions, prediction of S-Nitrosylation sites, and comparison of transcriptomic and qRT-PCR data. Results showed the diverse distribution of NBS-LRR genes across chromosomes, and variations in amino acid number, exons/introns, molecular weight, and theoretical isoelectric point, and they were found in various cellular locations like the plasma membrane, cytoplasm, and nucleus. These genes predominantly harbor the NB-ARC superfamily, LRR, LRR_8, and TIR domains, as also confirmed by motif analysis. Additionally, they feature species-specific PLN00113 superfamily and RX-CC_like domain in dicots and monocots, respectively, both responsive to defense against pathogen attacks. The NO-induced NBS-LRR genes of Arabidopsis reveal the presence of cis-elements responsive to phytohormones, light, stress, and growth, suggesting a wide range of responses mediated by NO. Protein-protein interactions, coupled with the prediction of S-Nitrosylation sites, offer valuable insights into the regulatory role of NO at the protein level within each respective species.

CONCLUSION

These above findings aimed to provide a thorough understanding of the impact of NO on NBS-LRR genes and their relationships with key plant species.

Genome and transcriptome exploration reveals receptor-like kinases as potential resistance gene analogs against bacterial blight in pomegranate

BACKGROUND

Pomegranate (Punica granatum L.) is a tropical fruit crop of pharma-nutritional importance. However, it faces farming challenges due to pests and diseases, particularly bacterial blight and wilt. Developing resistant cultivars is crucial for sustainable pomegranate cultivation, and understanding resistance's genetic basis is essential.

METHODS AND RESULTS

We used an extensive resistance gene analogues (RGA) prediction tool to identify 958 RGAs, classified into Nucleotide Binding Site-leucine-rich repeat (NBS-LRR) proteins, receptor-like kinases (RLKs), receptor-like proteins (RLPs), Transmembrane coiled-coil (TM-CC), and nine non-canonical RGAs. RGAs were distributed across all eight chromosomes, with chromosome 02 containing the most RGAs (161), and chromosome 08 having the highest density (4.42 RGA/Mb). NBS-LRR genes were predominantly present on chromosomes 08 and 02, whereas RLKs and RLPs were primarily located on chromosomes 04 and 07. Gene ontology analysis revealed that 475 RGAs were associated with defence against various biotic stresses. Using RNAseq, we identified 120 differentially expressed RGAs, with RLKs (74) being prominent among the differentially expressed genes.

CONCLUSION

The discovery of these RGAs is a significant step towards breeding pomegranates for pest and disease resistance. The differentially expressed RLKs hold promise for developing resistant cultivars against bacterial blight, thereby contributing to the sustainability of pomegranate cultivation.

Genome-wide identification and expression analysis of the MADS gene family in sweet orange (Citrus sinensis) infested with pathogenic bacteria

The risk of pathogenic bacterial invasion in plantations has increased dramatically due to high environmental climate change and has seriously affected sweet orange fruit quality. MADS genes allow plants to develop increased resistance, but functional genes for resistance associated with pathogen invasion have rarely been reported. MADS gene expression profiles were analyzed in sweet orange leaves and fruits infested with Lecanicillium psalliotae and Penicillium digitatum, respectively. Eighty-two MADS genes were identified from the sweet orange genome, and they were classified into five prime subfamilies concerning the Arabidopsis MADS gene family, of which the MIKC subfamily could be subdivided into 13 minor subfamilies. Protein structure analysis showed that more than 93% of the MADS protein sequences of the same subfamily between sweet orange and Arabidopsis were very similar in tertiary structure, with only CsMADS8 and AG showing significant differences. The variability of MADS genes protein structures between sweet orange and Arabidopsis subgroups was less than the variabilities of protein structures within species. Chromosomal localization and covariance analysis showed that these genes were unevenly distributed on nine chromosomes, with the most genes on chromosome 9 and the least on chromosome 2, with 36 and two, respectively. Four pairs of tandem and 28 fragmented duplicated genes in the 82 MADS gene sequences were found in sweet oranges. GO (Gene Ontology) functional enrichment and expression pattern analysis showed that the functional gene CsMADS46 was strongly downregulated of sweet orange in response to biotic stress adversity. It is also the first report that plants' MADS genes are involved in the biotic stress responses of sweet oranges. For the first time, L. psalliotae was experimentally confirmed to be the causal agent of sweet orange leaf spot disease, which provides a reference for the research and control of pathogenic L. psalliotae.

Comparative genomics and bioinformatics approaches revealed the role of CC-NBS-LRR genes under multiple stresses in passion fruit

Passion fruit is widely cultivated in tropical, subtropical regions of the world. The attack of bacterial and fungal diseases, and environmental factors heavily affect the yield and productivity of the passion fruit. The CC-NBS-LRR (CNL) gene family being a subclass of R-genes protects the plant against the attack of pathogens and plays a major role in effector-triggered immunity (ETI). However, no information is available regarding this gene family in passion fruit. To address the underlying problem a total of 25 and 21 CNL genes have been identified in the genome of purple (Passiflora edulis Sims.) and yellow (Passiflora edulis f. flavicarpa) passion fruit respectively. Phylogenetic tree was divided into four groups with PeCNLs present in 3 groups only. Gene structure analysis revealed that number of exons ranged from 1 to 9 with 1 being most common. Most of the PeCNL genes were clustered at the chromosome 3 and underwent strong purifying selection, expanded through segmental (17 gene pairs) and tandem duplications (17 gene pairs). PeCNL genes contained cis-elements involved in plant growth, hormones, and stress response. Transcriptome data indicated that PeCNL3, PeCNL13, and PeCNL14 were found to be differentially expressed under Cucumber mosaic virus and cold stress. Three genes were validated to be multi-stress responsive by applying Random Forest model of machine learning. To comprehend the biological functions of PeCNL proteins, their 3D structure and gene ontology (GO) enrichment analysis were done. Our research analyzed the CNL gene family in passion fruit to understand stress regulation and improve resilience. This study lays the groundwork for future investigations aimed at enhancing the genetic composition of passion fruit to ensure robust growth and productivity in challenging environments.

Genome-wide identification and characterization of NBLRR genes in finger millet (Eleusine coracana L.) and their expression in response to Magnaporthe grisea infection

BACKGROUND

The nucleotide binding site leucine rich repeat (NBLRR) genes significantly regulate defences against phytopathogens in plants. The genome-wide identification and analysis of NBLRR genes have been performed in several species. However, the detailed evolution, structure, expression of NBLRRs and functional response to Magnaporthe grisea are unknown in finger millet (Eleusine coracana (L.) Gaertn.).

RESULTS

The genome-wide scanning of the finger millet genome resulted in 116 NBLRR (EcNBLRRs1-116) encompassing 64 CC-NB-LRR, 47 NB-LRR and 5 CCR-NB-LRR types. The evolutionary studies among the NBLRRs of five Gramineae species, viz., purple false brome (Brachypodium distachyon (L.) P.Beauv.), finger millet (E. coracana), rice (Oryza sativa L.), sorghum (Sorghum bicolor L. (Moench)) and foxtail millet (Setaria italica (L.) P.Beauv.) showed the evolution of NBLRRs in the ancestral lineage of the target species and subsequent divergence through gene-loss events. The purifying selection (Ka/Ks < 1) shaped the expansions of NBLRRs paralogs in finger millet and orthologs among the target Gramineae species. The promoter sequence analysis showed various stress- and phytohormone-responsive cis-acting elements besides growth and development, indicating their potential role in disease defence and regulatory mechanisms. The expression analysis of 22 EcNBLRRs in the genotypes showing contrasting responses to Magnaporthe grisea infection revealed four and five EcNBLRRs in early and late infection stages, respectively. The six of these nine candidate EcNBLRRs proteins, viz., EcNBLRR21, EcNBLRR26, EcNBLRR30, EcNBLRR45, EcNBLRR55 and EcNBLRR76 showed CC, NB and LRR domains, whereas the EcNBLRR23, EcNBLRR32 and EcNBLRR83 showed NB and LRR somains.

CONCLUSION

The identification and expression analysis of EcNBLRRs showed the role of EcNBLRR genes in assigning blast resistance in finger millet. These results pave the foundation for in-depth and targeted functional analysis of EcNBLRRs through genome editing and transgenic approaches.

Editing Metabolism, Sex, and Microbiome: How Can We Help Poplar Resist Pathogens?

Poplar (Populus) is a genus of woody plants of great economic value. Due to the growing economic importance of poplar, there is a need to ensure its stable growth by increasing its resistance to pathogens. Genetic engineering can create organisms with improved traits faster than traditional methods, and with the development of CRISPR/Cas-based genome editing systems, scientists have a new highly effective tool for creating valuable genotypes. In this review, we summarize the latest research data on poplar diseases, the biology of their pathogens and how these plants resist pathogens. In the final section, we propose to plant male or mixed poplar populations; consider the genes of the MLO group, transcription factors of the WRKY and MYB families and defensive proteins BbChit1, LJAMP2, MsrA2 and PtDef as the most promising targets for genetic engineering; and also pay attention to the possibility of microbiome engineering.

Comparative genomics profiling revealed multi-stress responsive roles of the CC-NBS-LRR genes in three mango cultivars

The nucleotide-binding site-leucine-rich repeat (NBS-LRR) gene family is the largest group of disease resistance (R) genes in plants and is active in response to viruses, bacteria, and fungi usually involved in effector-triggered immunity (ETI). Pangenome-wide studies allow researchers to analyze the genetic diversity of multiple species or their members simultaneously, providing a comprehensive understanding of the evolutionary relationships and diversity present among them. The draft pan-genome of three Mangifera indica cultivars (Alphonso, Hong Xiang Ya, and Tommy atkins) was constructed and Presence/absence variants (PAVs) were filtered through the ppsPCP pipeline. As a result, 2823 genes and 5907 PAVs from H. Xiang Ya, and 1266 genes and 2098 PAVs from T. atkins were added to the reference genome. For the identification of CC-NBS-LRR (CNL) genes in these mango cultivars, this draft pan-genome study has successfully identified 47, 27, and 36 members in Alphonso, H. Xiang Ya, and T. atkins respectively. The phylogenetic analysis divided MiCNL proteins into four distinct subgroups. All MiCNL genes are unevenly distributed on chromosomes. Both tandem and segmental duplication events played a significant role in the expansion of the CNL gene family. These genes contain cis-elements related to light, stress, hormone, and development. The analysis of protein-protein interactions (PPI) revealed that MiCNL proteins interacted with other defense-responsive proteins. Gene Ontology (GO) analysis indicated that MiCNL genes play a role in defense mechanisms within the organism. The expression level of the identified genes in fruit peel was observed under disease and cold stress which showed that Mi_A_CNL13 and 14 were up-regulated while Mi_A_CNL15, 25, 30, 31, and 40 were down-regulated in disease stress. On the other hand, Mi_A_CNL2, 14, 41, and 45 were up-regulated and Mi_A_CNL47 is down-regulated in cold stress. Subsequently, the Random Forest (RF) classifier was used to assess the multi-stress response of MiCNLs. It was found that Mi_A_CNL14 is a gene that responds to multiple stress conditions. The CNLs have similar protein structures which show that they are involved in the same function. The above findings provide a foundation for a deeper understanding of the functional characteristics of the mango CNL gene family.

In silico structural-functional characterization of three differentially expressed resistance gene analogs identified in Dalbergia sissoo against dieback disease reveals their role in immune response regulation

Plant immunity includes enemy recognition, signal transduction, and defensive response against pathogens. We experimented to identify the genes that contribute resistance against dieback disease to Dalbergia sissoo, an economically important timber tree. In this study, we investigated the role of three differentially expressed genes identified in the dieback-induced transcriptome in Dalbergia sissoo. The transcriptome was probed using DOP-rtPCR analysis. The identified RGAs were characterized in silico as the contributors of disease resistance that switch on under dieback stress. Their predicted fingerprints revealed involvement in stress response. Ds-DbRCaG-02-Rga.a, Ds-DbRCaG-04-Rga.b, and Ds-DbRCaG-06-Rga.c showed structural homology with the Transthyretin-52 domain, EAL associated YkuI_C domain, and Src homology-3 domain respectively, which are the attributes of signaling proteins possessing a role in regulating immune responses in plants. Based on in-silico structural and functional characterization, they were predicted to have a role in immune response regulation in D. sissoo.

Identification, characterization, and validation of NBS-encoding genes in grass pea

Grass pea is a promising crop with the potential to provide food and fodder, but its genomics has not been adequately explored. Identifying genes for desirable traits, such as drought tolerance and disease resistance, is critical for improving the plant. Grass pea currently lacks known R-genes, including the nucleotide-binding site-leucine-rich repeat (NBS-LRR) gene family, which plays a key role in protecting the plant from biotic and abiotic stresses. In our study, we used the recently published grass pea genome and available transcriptomic data to identify 274 NBS-LRR genes. The evolutionary relationships between the classified genes on the reported plants and LsNBS revealed that 124 genes have TNL domains, while 150 genes have CNL domains. All genes contained exons, ranging from 1 to 7. Ten conserved motifs with lengths ranging from 16 to 30 amino acids were identified. We found TIR-domain-containing genes in 132 LsNBSs, with 63 TIR-1 and 69 TIR-2, and RX-CCLike in 84 LsNBSs. We also identified several popular motifs, including P-loop, Uup, kinase-GTPase, ABC, ChvD, CDC6, Rnase_H, Smc, CDC48, and SpoVK. According to the gene enrichment analysis, the identified genes undergo several biological processes such as plant defense, innate immunity, hydrolase activity, and DNA binding. In the upstream regions, 103 transcription factors were identified that govern the transcription of nearby genes affecting the plant excretion of salicylic acid, methyl jasmonate, ethylene, and abscisic acid. According to RNA-Seq expression analysis, 85% of the encoded genes have high expression levels. Nine LsNBS genes were selected for qPCR under salt stress conditions. The majority of the genes showed upregulation at 50 and 200 μM NaCl. However, LsNBS-D18, LsNBS-D204, and LsNBS-D180 showed reduced or drastic downregulation compared to their respective expression levels, providing further insights into the potential functions of LsNBSs under salt stress conditions. They provide valuable insights into the potential functions of LsNBSs under salt stress conditions. Our findings also shed light on the evolution and classification of NBS-LRR genes in legumes, highlighting the potential of grass pea. Further research could focus on the functional analysis of these genes, and their potential use in breeding programs to improve the salinity, drought, and disease resistance of this important crop.

Pistachio genomes provide insights into nut tree domestication and ZW sex chromosome evolution

Pistachio is a nut crop domesticated in the Fertile Crescent and a dioecious species with ZW sex chromosomes. We sequenced the genomes of Pistacia vera cultivar (cv.) Siirt, the female parent, and P. vera cv. Bagyolu, the male parent. Two chromosome-level reference genomes of pistachio were generated, and Z and W chromosomes were assembled. The ZW chromosomes originated from an autosome following the first inversion, which occurred approximately 8.18 Mya. Three inversion events in the W chromosome led to the formation of a 12.7-Mb (22.8% of the W chromosome) non-recombining region. These W-specific sequences contain several genes of interest that may have played a pivotal role in sex determination and contributed to the initiation and evolution of a ZW sex chromosome system in pistachio. The W-specific genes, including defA, defA-like, DYT1, two PTEN1, and two tandem duplications of six VPS13A paralogs, are strong candidates for sex determination or differentiation. Demographic history analysis of resequenced genomes suggest that cultivated pistachio underwent severe domestication bottlenecks approximately 7640 years ago, dating the domestication event close to the archeological record of pistachio domestication in Iran. We identified 390, 211, and 290 potential selective sweeps in 3 cultivar subgroups that underlie agronomic traits such as nut development and quality, grafting success, flowering time shift, and drought tolerance. These findings have improved our understanding of the genomic basis of sex determination/differentiation and horticulturally important traits and will accelerate the improvement of pistachio cultivars and rootstocks.

Genome-wide identification, characterization, and expression profile ofNBS-LRRgene family in sweet orange (Citrussinensis)

BACKGROUND

The NBS-LRR (nucleotide-binding site-leucine-rich repeat gene) gene family, known as the plant R (resistance) gene family with the most members, plays a significant role in plant resistance to various external adversity stresses. The NBS-LRR gene family has been researched in many plant species. Citrus is one of the most vital global cash crops, the number one fruit group, and the third most traded agricultural product world wild. However, as one of the largest citrus species, a comprehensive study of the NBS-LRR gene family has not been reported on sweet oranges.

METHODS

In this study, NBS-LRR genes were identified from the Citrus sinensis genome (v3.0), with a comprehensive analysis of this gene family performed, including phylogenetic analysis, gene structure, cis-acting element of a promoter, and chromosomal localization, among others. The expression pattern of NBS-LRR genes was analyzed when sweet orange fruits were infected by Penicillium digitatum, employing experimental data from our research group. It first reported the expression patterns of NBS-LRR genes under abiotic stresses, using three transcript data from NCBI (National Center for Biotechnology Information).

RESULTS

In this study, 111 NBS-LRR genes were identified in the C. sinensis genome (v3.0) and classified into seven subfamilies according to their N-terminal and C-terminal domains. The phylogenetic tree results indicate that genes containing only the NBS structural domain are more ancient in the sweet orange NBS-LRR gene family. The chromosome localization results showed that 111 NBS-LRR genes were distributed unevenly on nine chromosomes, with the most genes distributed on chromosome 1. In addition, we identified a total of 18 tandem duplication gene pairs in the sweet orange NBS-LRR gene family, and based on the Ka/Ks ratio, all of the tandem duplication genes underwent purifying selection. Transcriptome data analysis showed a significant number of NBS-LRR genes expressed under biotic and abiotic stresses, and some reached significantly different levels of expression. It indicates that the NBS-LRR gene family is vital in resistance to biotic and abiotic stresses in sweet oranges.

CONCLUSION

Our study provides the first comprehensive framework on the NBS-LRR family of genes, which provides a basis for further in-depth studies on the biological functions of NBS-LRR in growth, development, and response to abiotic stresses in sweet orange.

The responses of poplars to fungal pathogens: A review of the defensive pathway

Long-lived tree species need to cope with changing environments and pathogens during their lifetime. Fungal diseases cause damage to trees growth and forest nurseries. As model system for woody plants, poplars are also hosts of a large variety of fungus. The defense strategies to fungus are generally associated with the type of fungus, therefore, the defense strategies of poplar against necrotrophic and biotrophic fungus are different. Poplars initiate constitutive defenses and induced defenses based on recognition of the fungus, hormone signaling network cascades, activation of defense-related genes and transcription factors and production of phytochemicals. The means of sensing fungus invasion in poplars are similar with herbs, both of which are mediated by receptor proteins and resistance (R) proteins, leading to pattern-triggered immunity (PTI) and effector-triggered immunity (ETI), but poplars have evolved some unique defense mechanisms compared with Arabidopsis due to their longevity. In this paper, current researches on poplar defensive responses to necrotrophic and biotrophic fungus, which mainly include the physiological and genetic aspects, and the role of noncoding RNA (ncRNA) in fungal resistance are reviewed. This review also provides strategies to enhance poplar disease resistance and some new insights into future research directions.

Comprehensive Analysis and Functional Verification of the Pinus massoniana NBS-LRR Gene Family Involved in the Resistance to Bursaphelenchus xylophilus

Pinus massoniana Lamb. is a crucial timber and resin conifer in China, but its plantation industry is threatened by outbreaks of pine wilt disease (PWD) caused by Bursaphelenchus xylophilus (pinewood nematode; PWN). However, as of yet, there is no comprehensive analysis of NBS-LRR genes in P. massoniana involved in its defense against PWN. In this study, 507 NBS genes were identified in the transcriptome of resistant and susceptible P. masoniana inoculated with the PWN. The phylogenetic analysis and expression profiles of resistant and susceptible P. massoniana revealed that the up-regulated PmNBS-LRR97 gene was involved in conferring resistance to PWN. The results of real-time quantitative PCR (qRT-PCR) showed that PmNBS-LRR97 was significantly up-regulated after PWN infection, especially in the stems. Subcellular localization indicated that PmNBS-LRR97 located to the cell membrane. PmNBS-LRR97 significantly activated the expression of reactive oxygen species (ROS)-related genes in P. massoniana. In addition, the overexpression of PmNBS-LRR97 was capable of promoting the production of ROS, aiding in plant growth and development. In summary, PmNBS-LRR97 participates in the defense response to PWN and plays an active role in conferring resistance in P. massoniana. This finding provides new insight into the regulatory mechanism of the R gene in P. massoniana.

Bioinformatics analysis and function prediction of NBS-LRR gene family in Broussonetia papyrifera

Most of the currently available disease resistance (R) genes have NBS (nucleotide-binding site) and LRR (leucine-rich-repeat) domain which belongs to the NBS-LRR gene family. The whole genome sequencing of Broussonetia papyrifera provides an important bioinformatics database for the study of the NBS-LRR gene family. In this study, 328 NBS-LRR family genes were identified and classified in B. papyrifera according to different classification schemes, where there are 92 N types, 47 CN type, 54 CNL type, 29 NL types, 55 TN type, and 51 TNL type. Subsequently, we conducted bioinformatics analysis of the NBS-LRR gene family. Classification, motif analysis of protein sequences, and phylogenetic tree studies of the NBS-LRR genes in B. papyrifera provide important basis for the functional study of NBS-LRR family genes. Additionally, we performed structural analysis of the chromosomal location, physicochemical properties, and sequences identified by genetic characterization. In addition, through the analysis of GO enrichment, it was found that NBS-LRR genes were involved in defense responses and were significantly enriched in biological stimulation, immune response, and abiotic stress. In addition, we found that Bp06g0955 was the most sensitive to low temperature and encoded the RPM1 protein by analyzing the low temperature transcriptome data of B. papyrifera. Quantitative results of gene expression after 48 h of Fusarium infection showed that Bp01g3293 increased 14 times after infection, which encodes RPM1 protein. The potential of NBS-LRR gene responsive to biotic and abiotic stresses can be exploited to improve the resistance of B. papyrifera.

Analysis of R Genes Related to Blackcurrant Reversion Virus Resistance in the Comparative Transcriptome of Ribes nigrum cv. Aldoniai

Blackcurrant reversion virus (BRV) is the most destructive mite-transmitted pathogen in blackcurrants. The understanding of the resistance to BRV is limited, hindering and delaying the selection process. To identify the resistance (R) gene for BRV resistance, a gene expression analysis based on de novo blackcurrant cv. Aldoniai comparative transcriptome analysis (mock- and BRV-inoculated samples at 2 and 4 days post-inoculation (dpi)) was performed. In this study, 111 annotated clusters associated with pathogenesis according to conservative R gene domains were identified. In virus-infected samples, only Cluster-12591.33361 showed significant expression at 4 dpi. The expression profiles of this cluster were significantly associated with the presence of BRV particles in plant tissues, making it a putative R gene in the dominant resistance strategy in the BRV-Ribes nigrum interaction. The newly identified gene R.nigrum_R belongs to the CC-NBS-LRR class and has 63.9% identity with RPM1 in Populus spp. This study provides new insights on dominant putative R genes related to resistance to BRV in R. nigrum, which could aid targeted research and genetic improvement in breeding programs of blackcurrants.

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.

Genome-Wide Analysis of the Homeobox Gene Family and Identification of Drought-Responsive Members in Populus trichocarpa

Homeobox (HB) genes play critical roles in the regulation of plant morphogenesis, growth and development. Here, we identified a total of 156 PtrHB genes from the Populus trichocarpa genome. According to the topologies and taxonomy of the phylogenetic tree constructed by Arabidopsis thaliana HB members, all PtrHB proteins were divided into six subgroups, namely HD-ZIP, ZF-HD, HB-PHD, TALE, WOX and HB-OTHERS. Multiple alignments of conserved homeodomains (HDs) revealed the conserved loci of each subgroup, while gene structure analysis showed similar exon–intron gene structures, and motif analysis indicated the similarity of motif number and pattern in the same subgroup. Promoter analysis indicated that the promoters of PtrHB genes contain a series of cis-acting regulatory elements involved in responding to various abiotic stresses, indicating that PtrHBs had potential functions in these processes. Collinearity analysis revealed that there are 96 pairs of 127 PtrHB genes mainly distributing on Chromosomes 1, 2, and 5. We analyzed the spatio-temporal expression patterns of PtrHB genes, and the virus-induced gene silencing (VIGS) of PtrHB3 gene resulted in the compromised tolerance of poplar seedlings to mannitol treatment. The bioinformatics on PtrHB family and preliminary exploration of drought-responsive genes can provide support for further study of the family in woody plants, especially in drought-related biological processes. It also provides a direction for developing new varieties of poplar with drought resistance. Overall, our results provided significant information for further functional analysis of PtrHB genes in poplar and demonstrated that PtrHB3 is a dominant gene regulating tolerance to water stress treatment in poplar seedlings.

NB-LRR genes: characteristics in three Solanum species and transcriptional response to Ralstonia solanacearum in tomato

NB-LRR genes in the three Solanum species showed specific constitution characteristics and evolved multiple clusters and duplicates. Some genes could respond to biotic stresses such as tomato bacterial wilt. Nucleotide-binding and leucine-rich repeat (NB-LRR, NLR) is a largest resistance gene family in plants, which plays a key role in response to biotic stresses. In this study, NB-LRR genes in cultivated tomato Solanum lycopersicum (Sl) and its wild relatives S. pennellii (Spe) and S. pimpinellifolium (Spi) were analyzed using bioinformatics approaches. In total, 238, 202 and 217 NB-LRR genes of 8 different types were found in Sl, Spe and Spi, respectively. The three species showed similar genomic characteristics. The NB-LRR genes were mainly distributed on chromosomes 4, 5 and 11 and located at the distal zones, forming multiple clusters and tandem duplicates. A large number of homologs appeared through gene expansion, with most Ka/Ks values being less than 1, indicating that purifying selection had occurred in evolution. These genes were mainly expressed in root and could respond to different biotic stresses. RT-qPCR analysis revealed that SlNLR genes could respond to tomato bacterial wilt, with SlNLR1 probably involved in the resistance response, whereas others being the opposite. The transcription factors (TFs) and interaction proteins that regulate target genes were mainly Dof, NAC and MYB families and kinases. The results provide a basis for the isolation and application of related genes in plant disease resistance breeding.

Mass spectrometry-based forest tree metabolomics

Research in forest tree species has advanced slowly when compared with other agricultural crops and model organisms, mainly due to the long-life cycles, large genome sizes, and lack of genomic tools. Additionally, trees are complex matrices, and the presence of interferents (e.g., oleoresins and cellulose) challenges the analysis of tree tissues with mass spectrometry (MS)-based analytical platforms. In this review, advances in MS-based forest tree metabolomics are discussed. Given their economic and ecological significance, particular focus is given to Pinus, Quercus, and Eucalyptus forest tree species to better understand their metabolite responses to abiotic and biotic stresses in the current climate change scenario. Furthermore, MS-based metabolomics technologies produce large and complex datasets that require expertize to adequately manage, process, analyze, and store the data in dedicated repositories. To ensure that the full potential of forest tree metabolomics data are translated into new knowledge, these data should comply with the FAIR principles (i.e., Findable, Accessible, Interoperable, and Re-usable). It is essential that adequate standards are implemented to annotate metadata from forest tree metabolomics studies as is already required by many science and governmental agencies and some major scientific publishers. © 2019 John Wiley & Sons Ltd. Mass Spec Rev 40:126-157, 2021.

Label-Free Quantitative Proteomics Analysis in Susceptible and Resistant Brassica napus Cultivars Infected with Xanthomonas campestris pv. campestris

Black rot, caused by Xanthomonas campestris pv. campestris (Xcc), is the main disease of cruciferous vegetables. To characterize the resistance mechanism in the Brassica napus–Xcc pathosystem, Xcc-responsive proteins in susceptible (cv. Mosa) and resistant (cv. Capitol) cultivars were investigated using gel-free quantitative proteomics and analysis of gene expression. This allowed us to identify 158 and 163 differentially expressed proteins following Xcc infection in cv. Mosa and cv. Capitol, respectively, and to classify them into five major categories including antioxidative systems, proteolysis, photosynthesis, redox, and innate immunity. All proteins involved in protein degradation such as the protease complex, proteasome subunits, and ATP-dependent Clp protease proteolytic subunits, were upregulated only in cv. Mosa, in which higher hydrogen peroxide accumulation concurred with upregulated superoxide dismutase. In cv. Capitol, photosystem II (PS II)-related proteins were downregulated (excepting PS II 22 kDa), whereas the PS I proteins, ATP synthase, and ferredoxin-NADP⁺ reductase, were upregulated. For redox-related proteins, upregulation of thioredoxin, 2-cys peroxiredoxin, and glutathione S-transferase occurred in cv. Capitol, consistent with higher NADH-, ascorbate-, and glutathione-based reducing potential, whereas the proteins involved in the C₂ oxidative cycle and glycolysis were highly activated in cv. Mosa. Most innate immunity-related proteins, including zinc finger domain (ZFD)-containing protein, glycine-rich RNA-binding protein (GRP) and mitochondrial outer membrane porin, were highly enhanced in cv. Capitol, concomitant with enhanced expression of ZFD and GRP genes. Distinguishable differences in the protein profile between the two cultivars deserves higher importance for breeding programs and understanding of disease resistance in the B. napus–Xcc pathosystem.

Plant Thaumatin-like Proteins: Function, Evolution and Biotechnological Applications

Thaumatin-like proteins (TLPs) are a highly complex protein family associated with host defense and developmental processes in plants, animals, and fungi. They are highly diverse in angiosperms, for which they are classified as the PR-5 (Pathogenesis-Related-5) protein family. In plants, TLPs have a variety of properties associated with their structural diversity. They are mostly associated with responses to biotic stresses, in addition to some predicted activities under drought and osmotic stresses. The present review covers aspects related to the structure, evolution, gene expression, and biotechnological potential of TLPs. The efficiency of the discovery of new TLPs is below its potential, considering the availability of omics data. Furthermore, we present an exemplary bioinformatics annotation procedure that was applied to cowpea (Vigna unguiculata) transcriptome, including libraries of two tissues (root and leaf), and two stress types (biotic/abiotic) generated using different sequencing approaches. Even without using genomic sequences, the pipeline uncovered 56 TLP candidates in both tissues and stresses. Interestingly, abiotic stress (root dehydration) was associated with a high number of modulated TLP isoforms. The nomenclature used so far for TLPs was also evaluated, considering TLP structure and possible functions identified to date. It is clear that plant TLPs are promising candidates for breeding purposes and for plant transformation aiming a better performance under biotic and abiotic stresses. The development of new therapeutic drugs against human fungal pathogens also deserves attention. Despite that, applications derived from TLP molecules are still below their potential, as it is evident in our review.

Rapid Mining of Candidate Genes for Verticillium Wilt Resistance in Cotton Based on BSA-Seq Analysis

Cotton is a globally important cash crop. Verticillium wilt (VW) is commonly known as "cancer" of cotton and causes serious loss of yield and fiber quality in cotton production around the world. Here, we performed a BSA-seq analysis using an F_2:3 segregation population to identify the candidate loci involved in VW resistance. Two QTLs (qvw-D05-1 and qvw-D05-2) related to VW resistance in cotton were identified using two resistant/susceptible bulks from the F₂ segregation population constructed by crossing the resistant cultivar ZZM2 with the susceptible cultivar J11. A total of 30stop-lost SNPs and 42 stop-gained SNPs, which included 17 genes, were screened in the qvw-D05-2 region by SnpEff analysis. Further analysis of the transcriptome data and qRT-PCR revealed that the expression level of Ghir_D05G037630 (designated as GhDRP) varied significantly at certain time points after infection with V. dahliae. The virus-induced gene silencing of GhDRP resulted in higher susceptibility of the plants to V. dahliae than the control, suggesting that GhDRP is involved in the resistance to V. dahlia infection. This study provides a method for rapid mining of quantitative trait loci and screening of candidate genes, as well as enriches the genomic information and gene resources for the molecular breeding of disease resistance in cotton.

Comprehensive Genome-Wide Analysis of Thaumatin-Like Gene Family in Four Cotton Species and Functional Identification of GhTLP19 Involved in Regulating Tolerance to Verticillium dahlia and Drought

Thaumatin-like proteins (TLPs) present in the form of large multigene families play important roles in biotic stress and abiotic stress. However, there has been no systematic analysis of the TLPs in cotton. In this study, comprehensive identification and evolutionary analysis of TLPs in four species of cotton were conducted. In total, 50, 48, 91, and 90 homologous sequences were identified in Gossypium raimondii, G. arboreum, G. barbadense, and G. hirsutum, respectively. Gene structure, protein motifs, and gene expression were further investigated. Transcriptome and quantitative real-time PCR analysis indicated that GhTLPs participate in abiotic, biotic stress and cotton fiber development. GhTLP19 on chromosome At05 was selected as a candidate gene for further study. When GhTLP19 was silenced by virus-induced gene silencing (VIGS) in cotton, with the increase of malondialdehyde (MDA) content and the decrease of catalase (CAT) content, and as the increase of disease index (DI) and hyphae accumulation, the plants were more sensitive to drought and Verticillium dahliae. Furthermore, the GhTLP19 overexpressing Arabidopsis transgenic lines exhibited higher proline content, thicker and longer trichomes and more tolerance to drought when compared to wild type. This study will provide a basis and reference for future research on their roles in stress tolerance and fiber development.

Transcriptome analysis of contrasting resistance to herbivory by Empoasca fabae in two shrub willow species and their hybrid progeny

Short rotation woody biomass cultivars developed from fast-growing shrub species of willow (Salix spp.) have superior properties as perennial energy crops for the Northeast and Midwest US. However, the insect pest potato leafhopper (PLH) Empoasca fabae (Harris) can cause serious damage and reduce yield of susceptible genotypes. Currently, the willow cultivars in use display varying levels of susceptibility under PLH infestation. However, genes and markers for resistance to PLH are not yet available for marker-assisted selection in breeding. In this study, transcriptome differences between a resistant genotype 94006 (S. purpurea) and a susceptible cultivar ‘Jorr’ (S. viminalis), and their hybrid progeny were determined. Over 600 million RNA-Seq reads were generated and mapped to the Salix purpurea reference transcriptome. Gene expression analyses revealed the unique defense mechanism in resistant genotype 94006 that involves PLH-induced secondary cell wall modification. In the susceptible genotypes, genes involved in programed cell death were highly expressed, explaining the necrosis symptoms after PLH feeding. Overall, the discovery of resistance genes and defense mechanisms provides new resources for shrub willow breeding and research in the future.

Identification and Distribution of NBS-Encoding Resistance Genes of Dactylis glomerata L. and Its Expression Under Abiotic and Biotic Stress

Orchardgrass (Dactylis glomerata L.) is drought resistant and tolerant to barren landscapes, making it one of the most important forages for animal husbandry, as well as ecological restoration of rocky landscapes that are undergoing desertification. However, orchardgrass is susceptible to rust, which can significantly reduce its yield and quality. Therefore, understanding the genes that underlie resistance against rust in orchardgrass is critical. The evolution, cloning of plant disease resistance genes, and the analysis of pathogenic bacteria induced expression patterns are important contents in the study of interaction between microorganisms and plants. Genes with nucleotide binding site (NBS) structure are disease-resistant genes ubiquitous in plants and play an important role in plant attacks against various pathogens. Using sequence analysis and re-annotation, we identified 413 NBS resistance genes in orchardgrass. Similar to previous studies, NBS resistance genes containing TIR (toll/interleukin-1 receptor) domain were not found in orchardgrass. The NBS resistance genes can be divided into four types: NBS (up to 264 homologous genes, accounting for 64% of the total number of NBS genes in orchardgrass), NBS-LRR, CC-NBS, and CC-NBS-LRR (minimum of 26 homologous genes, only 6% of the total number of NBS genes in orchardgrass). These 413 NBS resistance genes were unevenly distributed across seven chromosomes where chromosome 5 had up to 99 NBS resistance genes. There were 224 (54%) NBS resistance genes expressed in different tissues (roots, stems, leaves, flowers, and spikes), and we did not detect expression for 45 genes (11%). The remaining 145 (35%) were expressed in some tissues. And we found that 11 NBS resistance genes were differentially expressed under waterlogging stress, 5 NBS resistance genes were differentially expressed under waterlogging and drought stress, and 1 NBS resistance was is differentially expressed under waterlogging and heat stress. Most importantly, we found that 65 NBS resistance genes were significantly expressed in different control groups. On the 7th day of inoculation, 23 NBS resistance genes were differentially expressed in high resistance materials alone, of which 7 NBS resistance genes regulate the "plant-pathogen interaction" pathway by encoding RPM1. At the same time, 2 NBS resistance genes that were differentially expressed in the high resistance material after inoculation were also differentially expressed in abiotic stress. In summary, the NBS resistance gene plays a crucial role in the resistance of orchardgrass to rust.

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.

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.

The chromosome-scale assembly of the willow genome provides insight into Salicaceae genome evolution

Salix suchowensis is an early-flowering shrub willow that provides a desirable system for studies on the basic biology of woody plants. The current reference genome of S. suchowensis was assembled with 454 sequencing reads. Here, we report a chromosome-scale assembly of S. suchowensis generated by combining PacBio sequencing with Hi-C technologies. The obtained genome assemblies covered a total length of 356 Mb. The contig N50 of these assemblies was 263,908 bp, which was ~65-fold higher than that reported previously. The contiguity and completeness of the genome were significantly improved. By applying Hi-C data, 339.67 Mb (95.29%) of the assembled sequences were allocated to the 19 chromosomes of haploid willow. With the chromosome-scale assembly, we revealed a series of major chromosomal fissions and fusions that explain the genome divergence between the sister genera of Salix and Populus. The more complete and accurate willow reference genome obtained in this study provides a fundamental resource for studying many genetic and genomic characteristics of woody plants.

Identification of Genes Underlying the Resistance to Melampsora larici-populina in an R Gene Supercluster of the Populus deltoides Genome

Identification of the particular genes in an R genes supercluster underlying resistance to the rust fungus Melampsora larici-populina in poplar genome remains challenging. Based on the de novo assembly of the Populus deltoides genome, all of the detected major genetic loci conferring resistance to M. larici-populina were confined to a 3.5-Mb region on chromosome 19. The transcriptomes of the resistant and susceptible genotypes were sequenced for a timespan from 0 to 168 hours postinoculation. By mapping the differentially expressed genes to the target genomic region, we identified two constitutive expression R genes and one inducible expression R gene that might confer resistance to M. larici-populina. Nucleotide variations were predicted based on the reconstructed haplotypes for each allele of the candidate genes. We also confirmed that salicylic acid was the phytohormone mediating signal transduction pathways, and PR-1 was identified as a key gene inhibiting rust reproduction. Finally, quantitative reverse transcription PCR assay revealed consistent expressions with the RNA-sequencing data for the detected key genes. This study presents an efficient approach for the identification of particular genes underlying phenotype of interest by the combination of genetic mapping, transcriptome profiling, and candidate gene sequences dissection. The identified key genes would be useful for host resistance diagnosis and for molecular breeding of elite poplar cultivars exhibiting resistance to M. larici-populina infection. The detected R genes are also valuable for testing whether the combination of individual R genes can induce durable quantitative resistance.

Genomic Organization and Comparative Phylogenic Analysis of NBS-LRR Resistance Gene Family in Solanum pimpinellifolium and Arabidopsis thaliana

NBS-LRR (nucleotide-binding site and leucine-rich repeat) is one of the largest resistance gene families in plants. The completion of the genome sequencing of wild tomato Solanum pimpinellifolium provided an opportunity to conduct a comprehensive analysis of the NBS-LRR gene superfamily at the genome-wide level. In this study, gene identification, chromosome mapping, and phylogenetic analysis of the NBS-LRR gene family were analyzed using the bioinformatics methods. The results revealed 245 NBS-LRRs in total, similar to that in the cultivated tomato. These genes are unevenly distributed on 12 chromosomes, and ~59.6% of them form gene clusters, most of which are tandem duplications. Phylogenetic analysis divided the NBS-LRRs into 2 subfamilies (CNL-coiled-coil NBS-LRR and TNL-TIR NBS-LRR), and the expansion of the CNL subfamily was more extensive than the TNL subfamily. Novel conserved structures were identified through conserved motif analysis between the CNL and TNL subfamilies. Compared with the NBS-LRR sequences from the model plant Arabidopsis thaliana, wide genetic variation occurred after the divergence of S. pimpinellifolium and A thaliana. Species-specific expansion was also found in the CNL subfamily in S. pimpinellifolium. The results of this study provide the basis for the deeper analysis of NBS-LRR resistance genes and contribute to mapping and isolation of candidate resistance genes in S. pimpinellifolium.

Genome-wide identification, characterization and expression analysis of the TLP gene family in melon (Cucumis melo L.)

Thaumatin-like proteins (TLPs), which belong to pathogenesis-related (PR) protein family 5 (PR5), are involved in plant host defense and various developmental processes. The functions of the TLP family have been extensively discussed in multiple organisms, whereas the detailed information of this family in melon has not been reported yet. In this study, we identified 28 TLP genes in the melon genome and a N-terminal signal peptide was found highly conserved within each member of this family. Phylogeny analysis indicated that TLPs from melon and other plant species were clustered into ten groups. Twelve segmental and seven tandem duplication gene pairs that underwent purifying selection were identified. TLP genes expressed differentially in different tissues/organs, and were significantly induced after Podosphaera xanthii infection. TLPs in breeding line MR-1 tend to express early after pathogen infection compared with cultivar Top Mark. Our study provides a comprehensive understanding of the melon TLP family and demonstrates their potential roles in disease resistance, therefore provides more reference for further research.

Genome-wide characterization revealed role of NBS-LRR genes during powdery mildew infection in Vitis vinifera

NBS-LRR comprises a large class of disease resistance (R) proteins that play a widespread role in plant protection against pathogens. In grapevine, powdery mildew cause significant losses in its productivity and efforts are being directed towards finding of resistance loci or genes imparting resistance/tolerance against such fungal diseases. In the present study, we performed genome-wide analysis of NBS-LRR genes during PM infection in grapevine. We identified 18, 23, 12, 16, 10, 10, 9, 20 and 14 differentially expressed NBS-LRR genes in response to PM infection in seven partially PM-resistant (DVIT3351.27, Husseine, Karadzhandal, Khalchili, Late vavilov, O34-16, Sochal) and 2 PM-susceptible (Carignan and Thompson seedless) V. vinifera accessions. Further, the identified sequences were characterized based on chromosomal locations, physicochemical properties, gene structure and motif analysis, and functional annotation by Gene Ontology (GO) mapping. The NBS-LRR genes responsive to powdery mildew could potentially be exploited to improve resistance in grapes.

Colonisation of Oncidium orchid roots by the endophyte Piriformospora indica restricts Erwinia chrysanthemi infection, stimulates accumulation of NBS-LRR resistance gene transcripts and represses their targeting micro-RNAs in leaves

BACKGROUND

Erwinia chrysanthemi (Ec) is a destructive pathogen which causes soft-rot diseases in diverse plant species including orchids. We investigated whether colonization of Oncidium roots by the endophytic fungus Piriformospora indica (Pi) restricts Ec-induced disease development in leaves, and whether this might be related to the regulation of nucleotide binding site-leucine rich repeat (NBS-LRR) Resistance (R) genes.

RESULTS

Root colonization of Oncidium stackings by Pi restricts progression of Ec-induced disease development in the leaves. Since Pi does not inhibit Ec growth on agar plates, we tested whether NBS-LRR R gene transcripts and the levels of their potential target miRNAs in Oncidium leaves might be regulated by Pi. Using bioinformatic tools, we first identified NBS-LRR R gene sequences from Oncidium, which are predicted to be targets of miRNAs. Among them, the expression of two R genes was repressed and the accumulation of several regulatory miRNA stimulated by Ec in the leaves of Oncidium plants. This correlated with the progression of disease development, jasmonic and salicylic acid accumulation, ethylene synthesis and H₂O₂ production after Ec infection of Oncidium leaves. Interestingly, root colonization by Pi restricted disease development in the leaves, and this was accompanied by higher expression levels of several defense-related R genes and lower expression level of their target miRNA.

CONCLUSION

Based on these data we propose that Pi controls the levels of NBS-LRR R mRNAs and their target miRNAs in leaves. This regulatory circuit correlates with the protection of Oncidium plants against Ec infection, and molecular and biochemical investigations will demonstrate in the future whether, and if so, to what extent these two observations are related to each other.

Joint linkage and association mapping of complex traits in shrub willow (Salix purpurea L.)

BACKGROUND AND AIMS

Increasing energy demands and the necessity to reduce greenhouse gas emissions are key motivating factors driving the development of lignocellulosic crops as an alternative to non-renewable energy sources. The effects of global climate change will require a better understanding of the genetic basis of complex adaptive traits to breed more resilient bioenergy feedstocks, like willow (Salix spp.). Shrub willow is a sustainable and dedicated bioenergy crop, bred to be fast-growing and high-yielding on marginal land without competing with food crops. In a rapidly changing climate, genomic advances will be vital for the sustained improvement of willow and other non-model bioenergy crops. Here, joint genetic mapping was used to exploit genetic variation garnered from both recent and historical recombination events in S. purpurea.

METHODS

A panel of North American naturalized S. purpurea accessions and full-sib F2S. purpurea population were genotyped and phenotyped for a suite of morphological, physiological, pest and disease resistance, and wood chemical composition traits, collected from multi-environment and multi-year replicated field trials. Controlling for population stratification and kinship in the association panel and spatial variation in the F2, a comprehensive mixed model analysis was used to dissect the complex genetic architecture and plasticity of these important traits.

KEY RESULTS

Individually, genome-wide association (GWAS) models differed in terms of power, but the combined approach, which corrects for yearly and environmental co-factors across datasets, improved the overall detection and resolution of associated loci. Although there were few significant GWAS hits located within support intervals of QTL for corresponding traits in the F2, many large-effect QTL were identified, as well as QTL hotspots.

CONCLUSIONS

This study provides the first comparison of linkage analysis and linkage disequilibrium mapping approaches in Salix, and highlights the complementarity and limits of these two methods for elucidating the genetic architecture of complex bioenergy-related traits of a woody perennial breeding programme.

The large repertoire of conifer NLR resistance genes includes drought responsive and highly diversified RNLs

The NLRs or NBS-LRRs (nucleotide-binding, leucine-rich-repeat) form the largest resistance gene family in plants, with lineage-specific contingents of TNL, CNL and RNL subfamilies and a central role in resilience to stress. The origin, evolution and distribution of NLR sequences has been unclear owing in part to the variable size and diversity of the RNL subfamily and a lack of data in Gymnosperms. We developed, searched and annotated transcriptomes assemblies of seven conifers and identified a resource of 3816 expressed NLR sequences. Our analyses encompassed sequences data spanning the major groups of land plants and determinations of NLR transcripts levels in response to drought in white spruce. We showed that conifers have among the most diverse and numerous RNLs in tested land plants. We report an evolutionary swap in the formation of RNLs, which emerged from the fusion of an RPW8 domain to a NB-ARC domain of CNL. We uncovered a quantitative relationship between RNLs and TNLs across all land plants investigated, with an average ratio of 1:10. The conifer RNL repertoire harbours four distinct groups, with two that differ from Angiosperms, one of which contained several upregulated sequences in response to drought while the majority of responsive NLRs are downregulated.

Stronger When Together: Clustering of Plant NLR Disease resistance Genes

Gene clustering is rare in eukaryotes. However, nucleotide-binding leucine-rich repeat (NLR)-encoding disease resistance (R) genes show consistent clustering in plant genomes. These arrangements are likely to provide coregulatory benefits, as suggested by growing evidence that the gene products of both paired and larger clusters of NLRs act together in triggering immunity. Head-to-head gene pairs where one of the encoded NLRs includes an integrated decoy domain appear to behave differently than clusters evolved from closely related typical NLRs. These patterns may help to explain the broad resistance that most plants have despite their finite number of R genes. By taking into consideration the relationship between genomic arrangement and function, we can improve our understanding of and ability to predict plant immune detection.

Comparative Genomics Analysis in Grass Species Reveals Two Distinct Evolutionary Strategies Adopted by R Genes

Resistance genes play an important role in the defense of plants against the invasion of pathogens. In Setaria italica and closely related grass species, R genes have been identified through genetic mapping and genome-wide homologous/domain searching. However, there has been to date no systematic analysis of the evolutionary features of R genes across all sequenced grass genomes. Here, we determined and comprehensively compared R genes in all 12 assembled grass genomes and an outgroup species (Arabidopsis thaliana) through synteny and selection analyses of multiple genomes. We found that the two groups of nucleotide binding site (NBS) domains containing R genes—R tandem duplications (TD) and R singletons—adopted different strategies and showed different features in their evolution. Based on K_a/K_s analysis between syntenic R loci pairs of TDs or singletons, we conclude that R singletons are under stronger purifying selection to be conserved among different grass species than R TDs, while R genes located at TD arrays have evolved much faster through diversifying selection. Furthermore, using the variome datasets of S. italica populations, we scanned for selection signals on genes and observed that a part of R singleton genes have been under purifying selection in populations of S. italica, which is consistent with the pattern observed in syntenic R singletons among different grass species. Additionally, we checked the synteny relationships of reported R genes in grass species and found that the functionally mapped R genes for novel resistance traits are prone to appear in TDs and are heavily divergent from their syntenic orthologs in other grass species, such the black streak R gene Rxo1 in Z. mays and the blast R gene Pi37 in O. sativa. These findings indicate that the R genes from TDs adopted tandem duplications to evolve faster and accumulate more mutations to facilitate functional innovation to cope with variable threats from a fluctuating environment, while R singletons provide a way for R genes to maintain sequence stability and retain conservation of function.

Sequencing of a Wild Apple (Malus baccata) Genome Unravels the Differences Between Cultivated and Wild Apple Species Regarding Disease Resistance and Cold Tolerance

Malus baccata is one of four wild apple species that can hybridize with the cultivated apple species (Malus domestica). It is widely used in high-latitude apple-producing areas as a rootstock and breeding resource because of its disease resistance, and cold tolerance. A lack of a reference genome has limited the application of M. baccata for apple breeding. We present a draft reference genome for M. baccata. The assembled sequence consisting of 665 Mb, with a scaffold N50 value of 452 kb, included transposable elements (413 Mb) and 46,114 high-quality protein-coding genes. According to a genetic map derived from 390 sibling lines, 72% of the assembly and 85% of the putative genes were anchored to 17 linkage groups. Many of the M. baccata genes under positive selection pressure were associated with plant–pathogen interaction pathways. We identified 2,345 Transcription factor-encoding genes in 58 families in the M. baccata genome. Genes related to disease defense and cold tolerance were also identified. A total of 462 putative nucleotide-binding site (NBS)-leucine-rich-repeat (LRR) genes, 177 Receptor-like kinase (RLK) and 51 receptor-like proteins (RLP) genes were identified in this genome assembly. The M. baccata genome contained 3978 cold-regulated genes, and 50% of these gene promoter containing DREB motif which can be induced by CBF gene. We herein present the first M. baccata genome assembly, which may be useful for exploring genetic variations in diverse apple germplasm, and for facilitating marker-assisted breeding of new apple cultivars exhibiting resistance to disease and cold stress.

The poplar pangenome provides insights into the evolutionary history of the genus

The genus Populus comprises a complex amalgam of ancient and modern species that has become a prime model for evolutionary and taxonomic studies. Here we sequenced the genomes of 10 species from five sections of the genus Populus, identified 71 million genomic variations, and observed new correlations between the single-nucleotide polymorphism-structural variation (SNP-SV) density and indel-SV density to complement the SNP-indel density correlation reported in mammals. Disease resistance genes (R genes) with heterozygous loss-of-function (LOF) were significantly enriched in the 10 species, which increased the diversity of poplar R genes during evolution. Heterozygous LOF mutations in the self-incompatibility genes were closely related to the self-fertilization of poplar, suggestive of genomic control of self-fertilization in dioecious plants. The phylogenetic genome-wide SNPs tree also showed possible ancient hybridization among species in sections Tacamahaca, Aigeiros, and Leucoides. The pangenome resource also provided information for poplar genetics and breeding.

Identification and expression profiling analysis of NBS-LRR genes involved in Fusarium oxysporum f.sp. conglutinans resistance in cabbage

As one of the most important resistance (R) gene families in plants, the NBS-LRR genes, encoding proteins with nucleotide-binding site (NBS) and leucine-rich repeat (LRR) domains, play significant roles in resisting pathogens. The published genomic data for cabbage (Brassica oleracea L.) provide valuable data to identify and characterize the genomic organization of cabbage NBS-LRR genes. Ultimately, we identified 105 TIR (N-terminal Toll/interleukin-1 receptor)-NBS-LRR (TNL) genes and 33 CC (coiled-coil)-NBS-LRR (CNL) genes. Further research indicated that 50.7% of the 138 NBS-LRR genes exist in 27 clusters and there are large differences among the gene structures and protein characteristics. Conserved motif and phylogenetic analysis showed that the structures of TNLs and CNLs were similar, with some differences. These NBS-LRRs are evolved under negative selection and mostly arose from whole-genome duplication events during evolution. Tissue-expression profiling of NBS-LRR genes revealed that 37.1% of the TNL genes are highly or specifically expressed in roots, especially the genes on chromosome 7 (76.5%). Digital gene expression and reverse transcription PCR analyses revealed the expression patterns of the NBS-LRR genes upon challenge by Fusarium oxysporum f.sp. conglutinans: nine genes were upregulated, and five were downregulated. The major resistance gene Foc1 probably works together with the other four genes in the same cluster to resist F. oxysporum infection.

Identification and Functional Analysis of NB-LRR-Type Virus Resistance Genes: Overview and Functional Analysis of Candidate Genes

Coexpression of a plant NB-LRR-type resistance (R) gene and corresponding viral avirulent (Avr) gene introduced in Nicotiana benthamiana using Agrobacterium tumefaciens confers hypersensitive response (HR). Such Agrobacterium-mediated transient gene expression methods have contributed to the identification of new plant R genes and facilitated the analysis of their functions. Here we describe a model method, by which several tobamovirus R genes from Solanaceous plants have been successfully identified and characterized molecularly.

Evolutionary Divergence of TNL Disease-Resistant Proteins in Soybean (Glycine max) and Common Bean (Phaseolus vulgaris)

Disease-resistant genes (R genes) encode proteins that are involved in protecting plants from their pathogens and pests. Availability of complete genome sequences from soybean and common bean allowed us to perform a genome-wide identification and analysis of the Toll interleukin-1 receptor-like nucleotide-binding site leucine-rich repeat (TNL) proteins. Hidden Markov model (HMM) profiling of all protein sequences resulted in the identification of 117 and 77 regular TNL genes in soybean and common bean, respectively. We also identified TNL gene homologs with unique domains, and signal peptides as well as nuclear localization signals. The TNL genes in soybean formed 28 clusters located on 10 of the 20 chromosomes, with the majority found on chromosome 3, 6 and 16. Similarly, the TNL genes in common bean formed 14 clusters located on five of the 11 chromosomes, with the majority found on chromosome 10. Phylogenetic analyses of the TNL genes from Arabidopsis, soybean and common bean revealed less divergence within legumes relative to the divergence between legumes and Arabidopsis. Syntenic blocks were found between chromosomes Pv10 and Gm03, Pv07 and Gm10, as well as Pv01 and Gm14. The gene expression data revealed basal level expression and tissue specificity, while analysis of available microRNA data showed 37 predicted microRNA families involved in targeting the identified TNL genes in soybean and common bean.

Genome-Wide Identification of NBS-Encoding Resistance Genes in Sunflower (Helianthus annuus L.)

Nucleotide Binding Site—Leucine-Rich Repeat (NBS-LRR) genes encode disease resistance proteins involved in plants’ defense against their pathogens. Although sunflower is affected by many diseases, only a few molecular details have been uncovered regarding pathogenesis and resistance mechanisms. Recent availability of sunflower whole genome sequences in publicly accessible databases allowed us to accomplish a genome-wide identification of Toll-interleukin-1 receptor-like Nucleotide-binding site Leucine-rich repeat (TNL), Coiled Coil (CC)-NBS-LRR (CNL), Resistance to powdery mildew8 (RPW8)-NBS-LRR (RNL) and NBS-LRR (NL) protein encoding genes. Hidden Markov Model (HMM) profiling of 52,243 putative protein sequences from sunflower resulted in 352 NBS-encoding genes, among which 100 genes belong to CNL group including 64 genes with RX_CC like domain, 77 to TNL, 13 to RNL, and 162 belong to NL group. We also identified signal peptides and nuclear localization signals present in the identified genes and their homologs. We found that NBS genes were located on all chromosomes and formed 75 gene clusters, one-third of which were located on chromosome 13. Phylogenetic analyses between sunflower and Arabidopsis NBS genes revealed a clade-specific nesting pattern in CNLs, with RNLs nested in the CNL-A clade, and species-specific nesting pattern for TNLs. Surprisingly, we found a moderate bootstrap support (BS = 50%) for CNL-A clade being nested within TNL clade making both the CNL and TNL clades paraphyletic. Arabidopsis and sunflower showed 87 syntenic blocks with 1049 high synteny hits between chromosome 5 of Arabidopsis and chromosome 6 of sunflower. Expression data revealed functional divergence of the NBS genes with basal level tissue-specific expression. This study represents the first genome-wide identification of NBS genes in sunflower paving avenues for functional characterization and potential crop improvement.

Comparative Genomics, Whole-Genome Re-sequencing and Expression Profile Analysis of Nucleobase:Cation Symporter 2 (NCS2) Genes in Maize

Nucleobase:cation symporter 2 (NCS2) proteins are important for the transport of free nucleobases, participating in diverse plant growth and developmental processes, as well as response to abiotic stress. To date, a comprehensive analysis of the NCS2 gene family has not been performed in maize. In this study, we conducted a comparative genomics analysis of NCS2 genes in 28 plant species, ranging from aquatic algae to land plants, concentrating mainly on maize. Gene duplication events contributed to the expansion of NCS2 genes from lower aquatic plants to higher angiosperms, and whole-genome/segmental and single-gene duplication events were responsible for the expansion of the maize NCS2 gene family. Phylogenetic construction showed three NCS2 subfamilies, I, II, and III. According to homology-based relationships, members of subfamily I are NCS2/AzgA-like genes, whereas those in subfamilies II and III are NCS2/NATs. Moreover, subfamily I exhibited ancient origins. A motif compositional analysis showed that one symbolic motif (motif 4) of the NCS2/NAT genes was absent in subfamily I. In maize, three NCS2/AzgA-like and 21 NCS2/NAT genes were identified, and purifying selection influenced the duplication of maize NCS2 genes. Additionally, a population genetic analysis of NCS2 genes revealed that ZmNCS2-21 showed the greatest diversity between the 78 inbred and 22 wild surveyed maize populations. An expression profile analysis using transcriptome data and quantitative real-time PCR revealed that NCS2 genes in maize are involved in diverse developmental processes and responses to abiotic stresses, including abscisic acid, salt (NaCl), polyethylene glycol, and low (4°C) and high (42°C) temperatures. ZmNCS2 genes with relatively close relationships had similar expression patterns, strongly indicating functional redundancy. Finally, ZmNCS2-16 and ZmNCS2-23 localize in the plasma membrane, which confirmed their predicted membrane structures. These results provide a foundation for future studies regarding the functions of ZmNCS2 proteins, particularly those with potentially important roles in plant responses to abiotic stresses.

Genome-Wide Screening and Characterization of the Dof Gene Family in Physic Nut (Jatropha curcas L.)

Physic nut (Jatropha curcas L.) is a species of flowering plant with great potential for biofuel production and as an emerging model organism for functional genomic analysis, particularly in the Euphorbiaceae family. DNA binding with one finger (Dof) transcription factors play critical roles in numerous biological processes in plants. Nevertheless, the knowledge about members, and the evolutionary and functional characteristics of the Dof gene family in physic nut is insufficient. Therefore, we performed a genome-wide screening and characterization of the Dof gene family within the physic nut draft genome. In total, 24 JcDof genes (encoding 33 JcDof proteins) were identified. All the JcDof genes were divided into three major groups based on phylogenetic inference, which was further validated by the subsequent gene structure and motif analysis. Genome comparison revealed that segmental duplication may have played crucial roles in the expansion of the JcDof gene family, and gene expansion was mainly subjected to positive selection. The expression profile demonstrated the broad involvement of JcDof genes in response to various abiotic stresses, hormonal treatments and functional divergence. This study provides valuable information for better understanding the evolution of JcDof genes, and lays a foundation for future functional exploration of JcDof genes.