Whole-genome sequencing and comparative genomics reveal candidate genes associated with quality traits in Dioscorea alata.

BACKGROUND: Quality traits are essential determinants of consumer preferences. Dioscorea alata (Greater Yam), is a starchy tuber crop in tropical regions. However, a comprehensive understanding of the genetic basis underlying yam tuber quality remains elusive. To address this knowledge gap, we employed population genomics and candidate gene association approaches to unravel the genetic factors influencing the quality attributes of boiled yam. METHODS AND RESULTS: Comparative genomics analysis of 45 plant species revealed numerous novel genes absent in the existing D. alata gene annotation. This approach, adding 48% more genes, significantly enhanced the functional annotation of three crucial metabolic pathways associated with boiled yam quality traits: pentose and glucuronate interconversions, starch and sucrose metabolism, and flavonoid biosynthesis. In addition, the whole-genome sequencing of 127 genotypes identified 27 genes under selection and 22 genes linked to texture, starch content, and color through a candidate gene association analysis. Notably, five genes involved in starch content and cell wall composition, including 1,3-beta Glucan synthase, ß-amylase, and Pectin methyl esterase, were common to both approaches and their expression levels were assessed by transcriptomic data. CONCLUSIONS: The analysis of the whole-genome of 127 genotypes of D. alata and the study of three specific pathways allowed the identification of important genes for tuber quality. Our findings provide insights into the genetic basis of yam quality traits and will help the enhancement of yam tuber quality through breeding programs.

Unzipped genome assemblies of polyploid root-knot nematodes reveal unusual and clade-specific telomeric repeats.

Using long-read sequencing, we assembled and unzipped the polyploid genomes of Meloidogyne incognita, M. javanica and M. arenaria, three of the most devastating plant-parasitic nematodes. We found the canonical nematode telomeric repeat to be missing in these and other Meloidogyne genomes. In addition, we find no evidence for the enzyme telomerase or for orthologs of C. elegans telomere-associated proteins, suggesting alternative lengthening of telomeres. Instead, analyzing our assembled genomes, we identify species-specific composite repeats enriched mostly at one extremity of contigs. These repeats are G-rich, oriented, and transcribed, similarly to canonical telomeric repeats. We confirm them as telomeric using fluorescent in situ hybridization. These repeats are mostly found at one single end of chromosomes in these species. The discovery of unusual and specific complex telomeric repeats opens a plethora of perspectives and highlights the evolutionary diversity of telomeres despite their central roles in senescence, aging, and chromosome integrity.

The Stilbene Synthase Family in Arachis: A Genome-Wide Study and Functional Characterization in Response to Stress.

Peanut (Arachis hypogaea) and its wild relatives are among the few species that naturally synthesize resveratrol, a well-known stilbenoid phytoalexin that plays a crucial role in plant defense against biotic and abiotic stresses. Resveratrol has received considerable attention due to its health benefits, such as preventing and treating various human diseases and disorders. Chalcone (CHS) and Stilbene (STS) Synthases are plant-specific type III Polyketide Synthases (PKSs) that share the same substrates and are key branch enzymes in the biosynthesis of flavonoids and stilbenoids, respectively. Although resveratrol accumulation in response to external stimulus has been described in peanut, there are no comprehensive studies of the CHS and STS gene families in the genus Arachis. In the present study, we identified and characterized 6 CHS and 46 STS genes in the tetraploid peanut and an average of 4 CHS and 22 STS genes in three diploid wild species (Arachis duranensis, Arachis ipaënsis and Arachis stenosperma). The CHS and STS gene and protein structures, chromosomal distributions, phylogenetic relationships, conserved amino acid domains, and cis-acting elements in the promoter regions were described for all Arachis species studied. Based on gene expression patterns of wild A. stenosperma STS genes in response to different biotic and abiotic stresses, we selected the candidate AsSTS4 gene, which is strongly induced by ultraviolet (UV) light exposure, for further functional investigation. The AsSTS4 overexpression in peanut hairy roots significantly reduced (47%) root-knot nematode infection, confirming that stilbene synthesis activation in transgenic plants can increase resistance to pathogens. These findings contribute to understanding the role of resveratrol in stress responses in Arachis species and provide the basis for genetic engineering for improved production of valuable secondary metabolites in plants.

Engineering Resistance against Sclerotinia sclerotiorum Using a Truncated NLR (TNx) and a Defense-Priming Gene.

The association of both cell-surface PRRs (Pattern Recognition Receptors) and intracellular receptor NLRs (Nucleotide-Binding Leucine-Rich Repeat) in engineered plants have the potential to activate strong defenses against a broad range of pathogens. Here, we describe the identification, characterization, and in planta functional analysis of a novel truncated NLR (TNx) gene from the wild species Arachis stenosperma (AsTIR19), with a protein structure lacking the C-terminal LRR (Leucine Rich Repeat) domain involved in pathogen perception. Overexpression of AsTIR19 in tobacco plants led to a significant reduction in infection caused by Sclerotinia sclerotiorum, with a further reduction in pyramid lines containing an expansin-like B gene (AdEXLB8) potentially involved in defense priming. Transcription analysis of tobacco transgenic lines revealed induction of hormone defense pathways (SA; JA-ET) and PRs (Pathogenesis-Related proteins) production. The strong upregulation of the respiratory burst oxidase homolog D (RbohD) gene in the pyramid lines suggests its central role in mediating immune responses in plants co-expressing the two transgenes, with reactive oxygen species (ROS) production enhanced by AdEXLB8 cues leading to stronger defense response. Here, we demonstrate that the association of potential priming elicitors and truncated NLRs can produce a synergistic effect on fungal resistance, constituting a promising strategy for improved, non-specific resistance to plant pathogens.

NBS-LRR-WRKY genes and protease inhibitors (PIs) seem essential for cowpea resistance to root-knot nematode.

Cowpea (Vigna unguiculata L. Walp) is a legume of great economic importance, however it is highly affected by nematodes. The present work aimed to identify proteins and genes involved in nematode resistance by proteomic and transcriptomic analysis. Plants of a genotype resistant (CE31) to root-knot nematode (Meloidogyne spp.) were collected 12 days after inoculation with Meloidogyne incognita and the total proteins and RNA were extracted from the root samples. Shotgun proteomic analysis was performed using an Orbitrap Elite mass spectrometer and the construction and sequencing of cDNA libraries were carried out in a Hi-Seq 2000 sequencing system. The proteomic and transcriptomic analyses revealed key processes involved in cowpea defense and some interesting candidates were further analyzed by RT-qPCR. Proteins and genes involved in essential biological processes were differentially accumulated such as, regulation of transcription, cell wall stiffening and microtubule-based process. However, the main defense strategies of Vigna unguiculata seem to be focused on the interaction of NBS-LRR and WRKY genes for the activation of R genes, production of protease inhibitors and maintenance of actin cytoskeleton. These are key processes that can culminate in the suppression of giant cell formation and consequently in the development of Meloidogyne incognita. SIGNIFICANCE: In this study, we identified proteins and transcripts regulated in cowpea resistant to the nematode Meloidogyne spp. upon inoculation. The results revealed key candidate genes involved in the activation of R genes, the production of protease inhibitors and maintenance of the actin cytoskeleton. These processes might be essential for cowpea resistance, as they can impede nematode nutrition, giant cell formation and consequently the development of Meloidogyne incognita.

Transcriptome Responses of Wild Arachis to UV-C Exposure Reveal Genes Involved in General Plant Defense and Priming.

Stress priming is an important strategy for enhancing plant defense capacity to deal with environmental challenges and involves reprogrammed transcriptional responses. Although ultraviolet (UV) light exposure is a widely adopted approach to elicit stress memory and tolerance in plants, the molecular mechanisms underlying UV-mediated plant priming tolerance are not fully understood. Here, we investigated the changes in the global transcriptome profile of wild Arachis stenosperma leaves in response to UV-C exposure. A total of 5751 differentially expressed genes (DEGs) were identified, with the majority associated with cell signaling, protein dynamics, hormonal and transcriptional regulation, and secondary metabolic pathways. The expression profiles of DEGs known as indicators of priming state, such as transcription factors, transcriptional regulators and protein kinases, were further characterized. A meta-analysis, followed by qRT-PCR validation, identified 18 metaDEGs as being commonly regulated in response to UV and other primary stresses. These genes are involved in secondary metabolism, basal immunity, cell wall structure and integrity, and may constitute important players in the general defense processes and establishment of a priming state in A. stenosperma. Our findings contribute to a better understanding of transcriptional dynamics involved in wild Arachis adaptation to stressful conditions of their natural habitats.

Defining the combined stress response in wild Arachis.

Nematodes and drought are major constraints in tropical agriculture and often occur simultaneously. Plant responses to these stresses are complex and require crosstalk between biotic and abiotic signaling pathways. In this study, we explored the transcriptome data of wild Arachis species subjected to drought (A-metaDEG) and the root-knot nematode Meloidogyne arenaria (B-metaDEG) via meta-analysis, to identify core-stress responsive genes to each individual and concurrent stresses in these species. Transcriptome analysis of a nematode/drought bioassay (cross-stress) showed that the set of stress responsive DEGs to concurrent stress is distinct from those resulting from overlapping A- and B-metaDEGs, indicating a specialized and unique response to combined stresses in wild Arachis. Whilst individual biotic and abiotic stresses elicit hormone-responsive genes, most notably in the jasmonic and abscisic acid pathways, combined stresses seem to trigger mainly the ethylene hormone pathway. The overexpression of a cross-stress tolerance candidate gene identified here, an endochitinase-encoding gene (AsECHI) from Arachis stenosperma, reduced up to 30% of M. incognita infection and increased post-drought recovery in Arabidopsis plants submitted to both stresses. The elucidation of the network of cross-stress responsive genes in Arachis contributes to better understanding the complex regulation of biotic and abiotic responses in plants facilitating more adequate crop breeding for combined stress tolerance.

Evolutionarily conserved plant genes responsive to root-knot nematodes identified by comparative genomics.

Root-knot nematodes (RKNs, genus Meloidogyne) affect a large number of crops causing severe yield losses worldwide, more specifically in tropical and sub-tropical regions. Several plant species display high resistance levels to Meloidogyne, but a general view of the plant immune molecular responses underlying resistance to RKNs is still lacking. Combining comparative genomics with differential gene expression analysis may allow the identification of widely conserved plant genes involved in RKN resistance. To identify genes that are evolutionary conserved across plant species, we used OrthoFinder to compared the predicted proteome of 22 plant species, including important crops, spanning 214 Myr of plant evolution. Overall, we identified 35,238 protein orthogroups, of which 6,132 were evolutionarily conserved and universal to all the 22 plant species (PLAnts Common Orthogroups-PLACO). To identify host genes responsive to RKN infection, we analyzed the RNA-seq transcriptome data from RKN-resistant genotypes of a peanut wild relative (Arachis stenosperma), coffee (Coffea arabica L.), soybean (Glycine max L.), and African rice (Oryza glaberrima Steud.) challenged by Meloidogyne spp. using EdgeR and DESeq tools, and we found 2,597 (O. glaberrima), 743 (C. arabica), 665 (A. stenosperma), and 653 (G. max) differentially expressed genes (DEGs) during the resistance response to the nematode. DEGs' classification into the previously characterized 35,238 protein orthogroups allowed identifying 17 orthogroups containing at least one DEG of each resistant Arachis, coffee, soybean, and rice genotype analyzed. Orthogroups contain 364 DEGs related to signaling, secondary metabolite production, cell wall-related functions, peptide transport, transcription regulation, and plant defense, thus revealing evolutionarily conserved RKN-responsive genes. Interestingly, the 17 DEGs-containing orthogroups (belonging to the PLACO) were also universal to the 22 plant species studied, suggesting that these core genes may be involved in ancestrally conserved immune responses triggered by RKN infection. The comparative genomic approach that we used here represents a promising predictive tool for the identification of other core plant defense-related genes of broad interest that are involved in different plant-pathogen interactions.

Contrasting Effects of Wild Arachis Dehydrin Under Abiotic and Biotic Stresses.

Plant dehydrins (DNHs) belong to the LEA (Late Embryogenesis Abundant) protein family and are involved in responses to multiple abiotic stresses. DHNs are classified into five subclasses according to the organization of three conserved motifs (K-; Y-; and S-segments). In the present study, the DHN protein family was characterized by molecular phylogeny, exon/intron organization, protein structure, and tissue-specificity expression in eight Fabaceae species. We identified 20 DHN genes, encompassing three (YnSKn, SKn, and Kn) subclasses sharing similar gene organization and protein structure. Two additional low conserved DHN Φ-segments specific to the legume SKn-type of proteins were also found. The in silico expression patterns of DHN genes in four legume species (Arachis duranensis, A. ipaënsis, Glycine max, and Medicago truncatula) revealed that their tissue-specific regulation is associated with the presence or absence of the Y-segment. Indeed, DHN genes containing a Y-segment are mainly expressed in seeds, whereas those without the Y-segment are ubiquitously expressed. Further qRT-PCR analysis revealed that, amongst stress responsive dehydrins, a SKn-type DHN gene from A. duranensis (AdDHN1) showed opposite response to biotic and abiotic stress with a positive regulation under water deficit and negative regulation upon nematode infection. Furthermore, transgenic Arabidopsis lines overexpressing (OE) AdDHN1 displayed improved tolerance to multiple abiotic stresses (freezing and drought) but increased susceptibility to the biotrophic root-knot nematode (RKN) Meloidogyne incognita. This contradictory role of AdDHN1 in responses to abiotic and biotic stresses was further investigated by qRT-PCR analysis of transgenic plants using a set of stress-responsive genes involved in the abscisic acid (ABA) and jasmonic acid (JA) signaling pathways and suggested an involvement of DHN overexpression in these stress-signaling pathways.

Comparative root transcriptome of wild Arachis reveals NBS-LRR genes related to nematode resistance.

BACKGROUND: The Root-Knot Nematode (RKN), Meloidogyne arenaria, significantly reduces peanut grain quality and yield worldwide. Whilst the cultivated species has low levels of resistance to RKN and other pests and diseases, peanut wild relatives (Arachis spp.) show rich genetic diversity and harbor high levels of resistance to many pathogens and environmental constraints. Comparative transcriptome analysis can be applied to identify candidate resistance genes. RESULTS: Transcriptome analysis during the early stages of RKN infection of two peanut wild relatives, the highly RKN resistant Arachis stenosperma and the moderately susceptible A. duranensis, revealed genes related to plant immunity with contrasting expression profiles. These included genes involved in hormone signaling and secondary metabolites production and also members of the NBS-LRR class of plant disease resistance (R) genes. From 345 NBS-LRRs identified in A.duranensis reference genome, 52 were differentially expressed between inoculated and control samples, with the majority occurring in physical clusters unevenly distributed on eight chromosomes with preferential tandem duplication. The majority of these NBS-LRR genes showed contrasting expression behaviour between A. duranensis and A. stenosperma, particularly at 6 days after nematode inoculation, coinciding with the onset of the Hypersensitive Response in the resistant species. The physical clustering of some of these NBS-LRR genes correlated with their expression patterns in the contrasting genotypes. Four NBS-LRR genes exclusively expressed in A. stenosperma are located within clusters on chromosome Aradu. A09, which harbors a QTL for RKN resistance, suggesting a functional role for their physical arrangement and their potential involvement in this defense response. CONCLUSION: The identification of functional novel R genes in wild Arachis species responsible for triggering effective defense cascades can contribute to the crop genetic improvement and enhance peanut resilience to RKN.

Early responses to dehydration in contrasting wild Arachis species.

Wild peanut relatives (Arachis spp.) are genetically diverse and were selected throughout evolution to a range of environments constituting, therefore, an important source of allelic diversity for abiotic stress tolerance. In particular, A. duranensis and A. stenosperma, the parents of the reference Arachis A-genome genetic map, show contrasting transpiration behavior under limited water conditions. This study aimed to build a comprehensive gene expression profile of these two wild species under dehydration stress caused by the withdrawal of hydroponic nutrient solution. For this purpose, roots of both genotypes were collected at seven time-points during the early stages of dehydration and used to construct cDNA paired-end libraries. Physiological analyses indicated initial differences in gas exchange parameters between the drought-tolerant genotype of A. duranensis and the drought-sensitive genotype of A. stenosperma. High-quality Illumina reads were mapped against the A. duranensis reference genome and resulted in the identification of 1,235 and 799 Differentially Expressed Genes (DEGs) that responded to the stress treatment in roots of A. duranensis and A. stenosperma, respectively. Further analysis, including functional annotation and identification of biological pathways represented by these DEGs confirmed the distinct gene expression behavior of the two contrasting Arachis species genotypes under dehydration stress. Some species-exclusive and common DEGs were then selected for qRT-PCR analysis, which corroborated the in silico expression profiling. These included genes coding for regulators and effectors involved in drought tolerance responses, such as activation of osmosensing molecular cascades, control of hormone and osmolyte content, and protection of macromolecules. This dataset of transcripts induced during the dehydration process in two wild Arachis genotypes constitute new tools for the understanding of the distinct gene regulation processes in these closely related species but with contrasting drought responsiveness. In addition, our findings provide insights into the nature of drought tolerance in wild germoplasm, which might be explored as novel sources of diversity and useful wild alleles to develop climate-resilient crop varieties.

Genome-wide analysis of expansin superfamily in wild Arachis discloses a stress-responsive expansin-like B gene.

Expansins are plant cell wall-loosening proteins involved in adaptive responses to environmental stimuli and various developmental processes. The first genome-wide analysis of the expansin superfamily in the Arachis genus identified 40 members in A. duranensis and 44 in A. ipaënsis, the wild progenitors of cultivated peanut (A. hypogaea). These expansins were further characterized regarding their subfamily classification, distribution along the genomes, duplication events, molecular structure, and phylogeny. A RNA-seq expression analysis in different Arachis species showed that the majority of these expansins are modulated in response to diverse stresses such as water deficit, root-knot nematode (RKN) infection, and UV exposure, with an expansin-like B gene (AraEXLB8) displaying a highly distinct stress-responsive expression profile. Further analysis of the AraEXLB8 coding sequences showed high conservation across the Arachis genotypes, with eight haplotypes identified. The modulation of AraEXLB8 expression in response to the aforementioned stresses was confirmed by qRT-PCR analysis in distinct Arachis genotypes, whilst in situ hybridization revealed transcripts in different root tissues according to the stress imposed. The overexpression of AraEXLB8 in soybean (Glycine max) composite plants remarkably decreased the number of galls in transformed hairy roots inoculated with RKN. This study improves the current understanding of the molecular evolution, divergence, and gene expression of expansins in Arachis, and provides molecular and functional insights into the role of expansin-like B, the less-studied plant expansin subfamily.