Reciprocal adaptation of rice and Xanthomonas oryzae pv. oryzae: cross-species 2D GWAS reveals the underlying genetics

Domestication and the evolution of crops: variable syndromes, complex genetic architectures, and ecological entanglements

Domestication can be considered a specialized mutualism in which a domesticator exerts control over the reproduction or propagation (fitness) of a domesticated species to gain resources or services. The evolution of crops by human-associated selection provides a powerful set of models to study recent evolutionary adaptations and their genetic bases. Moreover, the domestication and dispersal of crops such as rice, maize, and wheat during the Holocene transformed human social and political organization by serving as the key mechanism by which human societies fed themselves. Here we review major themes and identify emerging questions in three fundamental areas of crop domestication research: domestication phenotypes and syndromes, genetic architecture underlying crop evolution, and the ecology of domestication. Current insights on the domestication syndrome in crops largely come from research on cereal crops such as rice and maize, and recent work indicates distinct domestication phenotypes can arise from different domestication histories. While early studies on the genetics of domestication often identified single large-effect loci underlying major domestication traits, emerging evidence supports polygenic bases for many canonical traits such as shattering and plant architecture. Adaptation in human-constructed environments also influenced ecological traits in domesticates such as resource acquisition rates and interactions with other organisms such as root mycorrhizal fungi and pollinators. Understanding the ecological context of domestication will be key to developing resource-efficient crops and implementing more sustainable land management and cultivation practices.

Alternative polyadenylation profiles of susceptible and resistant rice (Oryza sativa L.) in response to bacterial leaf blight using RNA-seq

BACKGROUND

Alternative polyadenylation (APA) is an important pattern of post-transcriptional regulation of genes widely existing in eukaryotes, involving plant physiological and pathological processes. However, there is a dearth of studies investigating the role of APA profile in rice leaf blight.

RESULTS

In this study, we compared the APA profile of leaf blight-susceptible varieties (CT 9737-613P-M) and resistant varieties (NSIC RC154) following bacterial blight infection. Through gene enrichment analysis, we found that the genes of two varieties typically exhibited distal poly(A) (PA) sites that play different roles in two kinds of rice, indicating differential APA regulatory mechanisms. In this process, many disease-resistance genes displayed multiple transcripts via APA. Moreover, we also found five polyadenylation factors of similar expression patterns of rice, highlighting the critical roles of these five factors in rice response to leaf blight about PA locus diversity.

CONCLUSION

Notably, the present study provides the first dynamic changes of APA in rice in early response to biotic stresses and proposes a possible functional conjecture of APA in plant immune response, which lays the theoretical foundation for in-depth determination of the role of APA events in plant stress response and other life processes.

Uncovering the Genomic Regions Associated with Yield Maintenance in Rice Under Drought Stress Using an Integrated Meta-Analysis Approach

The complex trait of yield is controlled by several quantitative trait loci (QTLs). Given the global water deficit issue, the development of rice varieties suitable for non-flooded cultivation holds significant importance in breeding programs. The powerful approach of Meta-QTL (MQTL) analysis can be used for the genetic dissection of complicated quantitative traits. In the current study, a comprehensive MQTL analysis was conducted to identify consistent QTL regions associated with drought tolerance and yield-related traits under water deficit conditions in rice. In total, 1087 QTLs from 134 rice populations, published between 2000 to 2021, were utilized in the analysis. Distinct MQTL analysis of the relevant traits resulted in the identification of 213 stable MQTLs. The confidence interval (CI) for the detected MQTLs was between 0.12 and 19.7 cM. The average CI of the identified MQTLs (4.68 cM) was 2.74 times narrower compared to the average CI of the initial QTLs. Interestingly, 63 MQTLs coincided with SNP peak positions detected by genome-wide association studies for yield and drought tolerance-associated traits under water deficit conditions in rice. Considering the genes located both in the QTL-overview peaks and the SNP peak positions, 19 novel candidate genes were introduced, which are associated with drought response index, plant height, panicle number, biomass, and grain yield. Moreover, an inclusive MQTL analysis was performed on all the traits to obtain "Breeding MQTLs". This analysis resulted in the identification of 96 MQTLs with a CI ranging from 0.01 to 9.0 cM. The mean CI of the obtained MQTLs (2.33 cM) was 4.66 times less than the mean CI of the original QTLs. Thirteen MQTLs fulfilling the criteria of having more than 10 initial QTLs, CI < 1 cM, and an average phenotypic variance explained greater than 10%, were designated as "Breeding MQTLs". These findings hold promise for assisting breeders in enhancing rice yield under drought stress conditions.

Broad-spectrum antagonistic potential of Bacillus spp. volatiles against Rhizoctonia solani and Xanthomonas oryzae pv. oryzae

Rhizoctonia solani and Xanthomonas oryzae pv. oryzae (Xoo) are the two major diseases affecting the quality and quantity of rice production. In the current study, volatile organic compounds (VOCs) of Bacillus spp. were used as green biocontrol agents for plant diseases. In in vitro experiments, Bacillus spp. FZB42, NMTD17, and LLTC93-VOCs displayed strong antimicrobial volatile activity with inhibition rates of 76, 66, and 78% for R. solani and 78, 81, and 76% for Xoo, respectively, compared to control. The synthetic VOCs, namely Pentadecane (PDC), Benzaldehyde (BDH), 1,2-Benz isothiazol-3(2H)-one (1,2-BIT), and mixture (MIX) of VOCs showed high volatile activity with inhibition rates of 86, 86, 89, and 92% against R. solani and 81, 81, 82, and 86%, respectively, against Xoo as compared to control. In addition, the scanning and transmission electron microscopes (SEM and TEM) analyses were performed to examine the effect of Bacillus and synthetic VOC treatments on R. solani and Xoo morphology. The analysis revealed the deformed and irregularized morphology of R. solani mycelia and Xoo cells after VOC treatments. The microscopic analysis showed that the rapid inhibition was due to severe oxidative productions inside the R. solani mycelia and Xoo cells. By using molecular docking, it was determined that the synthetic VOCs entered the active binding site of trehalase and NADH dehydrogenase proteins, causing R. solani and Xoo cells to die prematurely and an accumulation of ROS. In the greenhouse experiment, FZB42, NMTD17, and LLTC93-VOCs significantly reduced the lesions of R. solani 8, 7, and 6 cm, and Xoo 7, 6, and 6 cm, respectively, then control. The synthetic VOCs demonstrated that the PDC, BDH, 1,2-BIT, and MIX-VOCs significantly reduced R. solani lesions on leaves 6, 6, 6, and 5 cm and Xoo 6, 5, 5, and 4 cm, respectively, as compared to control. Furthermore, plant defence-related genes and antioxidant enzymes were upregulated in rice plants. These findings provide novel mechanisms by which Bacillus antimicrobial VOCs control plant diseases.

Genome-Wide Association Study Identifies Resistance Loci for Bacterial Blight in a Collection of Asian Temperate Japonica Rice Germplasm

Growing resistant rice cultivars is the most effective strategy to control bacterial blight (BB), a devastating disease caused by Xanthomonas oryzae pv. oryzae (Xoo). Screening resistant germplasm and identifying resistance (R) genes are prerequisites for breeding resistant rice cultivars. We conducted a genome-wide association study (GWAS) to detect quantitative trait loci (QTL) associated with BB resistance using 359 East Asian temperate Japonica accessions inoculated with two Chinese Xoo strains (KS6-6 and GV) and one Philippine Xoo strain (PXO99^A). Based on the 55K SNPs Array dataset of the 359 Japonica accessions, eight QTL were identified on rice chromosomes 1, 2, 4, 10, and 11. Four of the QTL coincided with previously reported QTL, and four were novel loci. Six R genes were localized in the qBBV-11.1, qBBV-11.2, and qBBV-11.3 loci on chromosome 11 in this Japonica collection. Haplotype analysis revealed candidate genes associated with BB resistance in each QTL. Notably, LOC_Os11g47290 in qBBV-11.3, encoding a leucine-rich repeat receptor-like kinase, was a candidate gene associated with resistance to the virulent strain GV. Knockout mutants of Nipponbare with the susceptible haplotype of LOC_Os11g47290 exhibited significantly improved BB resistance. These results will be useful for cloning BB resistance genes and breeding resistant rice cultivars.

Comparative transcriptomics and co-expression networks reveal cultivar-specific molecular signatures associated with reproductive-stage cold stress in rice

The resistance of Huaidao5 results from the high constitutive expression of tolerance genes, while that of Huaidao9 is due to the cold-induced resistance in flag leaves and panicles. The regulation mechanism of rice seedlings' cold tolerance is relatively clear, and knowledge of its underlying mechanisms at the reproductive stage is limited. We performed differential expression and co-expression network analyses to transcriptomes from panicle and flag leaf tissues of a cold-tolerant cultivar (Huaidao5), and a sensitive cultivar (Huaidao9), under reproductive-stage cold stress. The results revealed that the expression levels of genes in stress-related pathways such as MAPK signaling pathway, diterpenoid biosynthesis, glutathione metabolism, plant-pathogen interaction and plant hormone signal transduction were constitutively highly expressed in Huaidao5, especially in panicles. Moreover, the Hudaidao5's panicle sample-specific (under cold) module contained some genes related to rice yield, such as GW5L, GGC2, SG1 and CTPS1. However, the resistance of Huaidao9 was derived from the induced resistance to cold in flag leaves and panicles. In the flag leaves, the responses included a series of stress response and signal transduction, while in the panicles nitrogen metabolism was severely affected, especially 66 endosperm-specific genes. Through integrating differential expression with co-expression networks, we predicted 161 candidate genes (79 cold-responsive genes common to both cultivars and 82 cold-tolerance genes associated with differences in cold tolerance between cultivars) potentially affecting cold response/tolerance, among which 85 (52.80%) were known to be cold-related genes. Moreover, 52 (65.82%) cold-responsive genes (e.g., TIFY11C, LSK1 and LPA) could be confirmed by previous transcriptome studies and 72 (87.80%) cold-tolerance genes (e.g., APX5, OsFbox17 and OsSTA109) were located within QTLs associated with cold tolerance. This study provides an efficient strategy for further discovery of mechanisms of cold tolerance in rice.

Genome-wide association studies in plant pathosystems: success or failure?

Harnessing natural genetic variation is an established alternative to artificial genetic variation for investigating the molecular dialog between partners in plant pathosystems. Herein, we review the successes of genome-wide association studies (GWAS) in both plants and pathogens. While GWAS in plants confirmed that the genetic architecture of disease resistance is polygenic, dynamic during the infection kinetics, and dependent on the environment, GWAS shortened the time of identification of quantitative trait loci (QTLs) and revealed both complex epistatic networks and a genetic architecture dependent upon the geographical scale. A similar picture emerges from the few GWAS in pathogens. In addition, the ever-increasing number of functionally validated QTLs has revealed new molecular plant defense mechanisms and pathogenicity determinants. Finally, we propose recommendations to better decode the disease triangle.

Protocatechuic acid, ferulic acid and relevant defense enzymes correlate closely with walnut resistance to Xanthomonas arboricola pv. juglandis

BACKGROUND

Juglans regia L. is an important nut tree that has a wide range of distribution in temperate regions of the world. In some walnut orchards, walnut blight can become a problematic disease that affects the growth of walnut trees. To explore the correlation between biochemical response and walnut resistance, we inoculated four walnut cultivars with Xanthomonas arboricola pv. juglandis (Xaj). The walnut cultivars were, namely, 'Xiangling', 'Xiluo 2', 'Yuanfeng' and 'Xifu 2'. Total phenol content (TPC) and total flavonoid content (TFC) were measured, whereby nine major phenolic compounds and several relevant enzymes were identified.

RESULTS

The results showed that the most resistant and susceptible walnut varieties were 'Xiluo 2' and 'Xifu 2' respectively. The reaction of walnut to Xaj was characterized by the early accumulation of phenolic compounds in the infected site. After inoculation with Xaj, we found that the resistant variety 'Xiluo 2' show the significant differences with other varieties at different time points through the determination of related antioxidant enzymes such as catalase (CAT) and peroxidase (POD). Meanwhile, the phenylalanine ammonia lyase (PAL) of 'Xiluo 2' increased significantly at 8 day post infection (dpi) and made differences from the control samples, while other varieties changed little. And the polyphenol oxidase (PPO) was significantly higher than in the control at 16 dpi, maintaining the highest and the lowest activity in 'Xiluo 2' and 'Xifu 2' respectively. It was also found that the content of protocatechuic acid in all cultivars increased significantly at 4 dpi, and 'Xiluo 2' was significantly higher than that of the control. In the early stage of the disease, ferulic acid content increased significantly in 'Xiluo 2'.

CONCLUSION

Our findings confirmed that the metabolism of phenolic compounds and related defense enzymes are of great significance in the response of walnut to Xaj.

Integrating plant and fungal quantitative genetics to improve the ecological and agricultural applications of mycorrhizal symbioses

Finding and targeting genes that quantitatively contribute to agricultural and ecological processes progresses food production and conservation efforts. Typically, quantitative genetic approaches link variants in a single organism's genome with a trait of interest. Recently, genome-to-genome mapping has found genome variants interacting between species to produce the result of a multiorganism (including multikingdom) interaction. These were plant and bacterial pathogen genome interactions; plant-fungal coquantitative genetics have not yet been applied. Plant-mycorrhizae symbioses exist across most biomes, for a majority of land plants, including crop plants, and manipulate many traits from single organisms to ecosystems for which knowing the genetic basis would be useful. The availability of Rhizophagus irregularis mycorrhizal isolates, with genomic information, makes dual-genome methods with beneficial mutualists accessible and imminent.

CRISPR/Cas-mediated editing of cis-regulatory elements for crop improvement

To improve future agricultural production, major technological advances are required to increase crop production and yield. Targeting the coding region of genes via the Clustered Regularly Interspaced Short Palindromic Repeats/CRISPR-associated Protein (CRISPR/Cas) system has been well established and has enabled the rapid generation of transgene-free plants, which can lead to crop improvement. The emergence of the CRISPR/Cas system has also enabled scientists to achieve cis-regulatory element (CRE) editing and, consequently, engineering endogenous critical CREs to modulate the expression of target genes. Recent genome-wide association studies have identified the domestication of natural CRE variants to regulate complex agronomic quantitative traits and have allowed for their engineering via the CRISPR/Cas system. Although engineering plant CREs can be advantageous to drive gene expression, there are still many limitations to its practical application. Here, we review the current progress in CRE editing and propose future strategies to effectively target CREs for transcriptional regulation for crop improvement.

A Genome-Wide Association study in Arabidopsis thaliana to decipher the adaptive genetics of quantitative disease resistance in a native heterogeneous environment

Pathogens are often the main selective agents acting in plant communities, thereby influencing the distribution of polymorphism at loci affecting resistance within and among natural plant populations. In addition, the outcome of plant-pathogen interactions can be drastically affected by abiotic and biotic factors at different spatial and temporal grains. The characterization of the adaptive genetic architecture of disease resistance in native heterogeneous environments is however still missing. In this study, we conducted an in situ Genome-Wide Association study in the spatially heterogeneous native habitat of a highly genetically polymorphic local mapping population of Arabidopsis thaliana, to unravel the adaptive genetic architecture of quantitative disease resistance. Disease resistance largely differed among three native soils and was affected by the presence of the grass Poa annua. The observation of strong crossing reactions norms among the 195 A. thaliana genotypes for disease resistance among micro-habitats, combined with a negative fecundity-disease resistance relationship in each micro-habitat, suggest that alternative local genotypes of A. thaliana are favored under contrasting environmental conditions at the scale of few meters. A complex genetic architecture was detected for disease resistance and fecundity. However, only few QTLs were common between these two traits. Heterogeneous selection in this local population should therefore promote the maintenance of polymorphism at only few candidate resistance genes.

A super pan-genomic landscape of rice

Pan-genomes from large natural populations can capture genetic diversity and reveal genomic complexity. Using de novo long-read assembly, we generated a graph-based super pan-genome of rice consisting of a 251-accession panel comprising both cultivated and wild species of Asian and African rice. Our pan-genome reveals extensive structural variations (SVs) and gene presence/absence variations. Additionally, our pan-genome enables the accurate identification of nucleotide-binding leucine-rich repeat genes and characterization of their inter- and intraspecific diversity. Moreover, we uncovered grain weight-associated SVs which specify traits by affecting the expression of their nearby genes. We characterized genetic variants associated with submergence tolerance, seed shattering and plant architecture and found independent selection for a common set of genes that drove adaptation and domestication in Asian and African rice. This super pan-genome facilitates pinpointing of lineage-specific haplotypes for trait-associated genes and provides insights into the evolutionary events that have shaped the genomic architecture of various rice species.

Current Methods for Recombination Detection in Bacteria

The role of genetic exchanges, i.e., homologous recombination (HR) and horizontal gene transfer (HGT), in bacteria cannot be overestimated for it is a pivotal mechanism leading to their evolution and adaptation, thus, tracking the signs of recombination and HGT events is importance both for fundamental and applied science. To date, dozens of bioinformatics tools for revealing recombination signals are available, however, their pros and cons as well as the spectra of solvable tasks have not yet been systematically reviewed. Moreover, there are two major groups of software. One aims to infer evidence of HR, while the other only deals with horizontal gene transfer (HGT). However, despite seemingly different goals, all the methods use similar algorithmic approaches, and the processes are interconnected in terms of genomic evolution influencing each other. In this review, we propose a classification of novel instruments for both HR and HGT detection based on the genomic consequences of recombination. In this context, we summarize available methodologies paying particular attention to the type of traceable events for which a certain program has been designed.

New Data and New Features of the FunRiceGenes (Functionally Characterized Rice Genes) Database: 2021 Update

As a major food crop and model organism, rice has been mostly studied with the largest number of functionally characterized genes among all crops. We previously built the funRiceGenes database including ~ 2800 functionally characterized rice genes and ~ 5000 members of different gene families. Since being published, the funRiceGenes database has been accessed by more than 54,400 users with over 540,000 pageviews. The funRiceGenes database has been continuously updated with newly cloned rice genes and newly published literature, based on the progress of rice functional genomics studies. Up to Nov 2021, ~ 4100 functionally characterized rice genes and ~ 6000 members of different gene families were collected in funRiceGenes, accounting for 22.3% of the 39,045 annotated protein-coding genes in the rice genome. Here, we summarized the update of the funRiceGenes database with new data and new features in the last 5 years.

Twenty years of rice genomics research: From sequencing and functional genomics to quantitative genomics

Since the completion of the rice genome sequencing project in 2005, we have entered the era of rice genomics, which is still in its ascendancy. Rice genomics studies can be classified into three stages: structural genomics, functional genomics, and quantitative genomics. Structural genomics refers primarily to genome sequencing for the construction of a complete map of rice genome sequence. This is fundamental for rice genetics and molecular biology research. Functional genomics aims to decode the functions of rice genes. Quantitative genomics is large-scale sequence- and statistics-based research to define the quantitative traits and genetic features of rice populations. Rice genomics has been a transformative influence on rice biological research and contributes significantly to rice breeding, making rice a good model plant for studying crop sciences.