Identification of a genomic region controlling thermotolerance at flowering in maize using a combination of whole genomic re-sequencing and bulked segregant analysis

Loss-of-function of LIGULELESS1 activates the jasmonate pathway and promotes maize resistance to corn leaf aphids

Corn leaf aphids (Rhopalosiphum maidis) are highly destructive pests of maize (Zea mays) that threaten growth and seed yield, but resources for aphid resistance are scarce. Here, we identified an aphid-resistant maize mutant, resistance to aphids 1 (rta1), which is allelic to LIGULELESS1 (LG1). We confirmed LG1's role in aphid resistance using the independent allele lg1-2, allelism tests and LG1 overexpression lines. LG1 interacts with, and increases the stability of ZINC-FINGER PROTEIN EXPRESSED IN INFLORESCENCE MERISTEM (ZIM1), a central component of the jasmonic acid (JA) signalling pathway, by disturbing its interaction with the F-box protein CORONATINE INSENSITIVE 1a (COI1a). Natural variation in the LG1 promoter was associated with aphid resistance among inbred lines. Moreover, a loss-of-function mutant in the LG1-related gene SPL8 in the dicot Arabidopsis thaliana conferred aphid resistance. This study revealed the aphid resistance mechanism of lg1, providing a theoretical basis and germplasm for breeding aphid-resistant crops.

A SLAF-based high-density genetic map construction and genetic architecture of thermotolerant traits in maize (Zea mays L.)

The leaf scorching trait at flowering is a crucial thermosensitive phenotype in maize under high temperature stress (HS), yet the genetic basis of this trait remains poorly understood. In this study, we genotyped a 254 RIL-F2:8 population, derived from the leaf scorch-free parental inbred line Abe2 and the leaf scorching maternal inbred line B73, using the specific-locus amplified fragment sequencing (SLAF-seq) method. A total of 10,112 polymorphic SLAF markers were developed, and a high-density genetic map with a total length of 1,475.88 cM was constructed. The average sequencing depth of the parents was 55.23X, and that of the progeny was 12.53X. Then, we identified a total of 16 QTLs associated with thermotolerant traits at flowering, of which four QTLs of leaf scorching damage (LS) were distributed on chromosomes 1 (qLS1), 2 (qLS2.1, qLS2.2) and 3 (qLS3), which could explain 19.73% of phenotypic variation. Combining one qLS1 locus with QTL-seq results led to the identification of 6 candidate genes. Expression experiments and sequence variation indicated that Zm00001d033328, encoding N-acetyl-gamma-glutamyl-phosphate reductase, was the most likely candidate gene controlling thermotolerant traits at flowering. In summary, the high-density genetic map and genetic basis of thermotolerant traits lay a critical foundation for mapping other complex traits and identifying the genes associated with thermotolerant traits in maize.

Genes and pathways correlated with heat stress responses and heat tolerance in maize kernels

Global warming leads to frequent extreme weather, especially the extreme heat events, which threating the safety of maize production. Here we selected a pair of maize inbred lines, PF5411-1 and LH150, with significant differences in heat tolerance at kernel development stage. The two maize inbred lines were treated with heat stress at kernel development stage. Compared with the control groups, transcriptomic analysis identified 770 common up- and down-regulated genes between PF5411-1 and LH150 under heat stress conditions, and 41 putative TFs were predicted. Based on the interaction term of the two-factorial design, we also identified 6,744 differentially regulated genes between LH150 and PF5411-1, 111 common up-regulated and 141 common down-regulated genes were overlapped with the differentially regulated genes, respectively. Combined with proteins and metabolites data, several key pathways including seven differentially regulated genes were highly correlated with the heat tolerance of maize kernels. The first is the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway ko04141: protein processing in endoplasmic reticulum, four small heat shock protein (sHSP) genes were enriched in this pathway, participating with the process of ER-associated degradation (ERAD). The second one is the myricetin biosynthesis pathway, a differentially regulated protein, flavonoid 3',5'-hydroxylase [EC:1.14.14.81], catalyzed the synthesis of myricetin. The third one is the raffinose metabolic pathway, one differentially regulated gene encoded the raffinose synthase controlled the synthesis of raffinose, high level of raffinose enhances the heat tolerance of maize kernels. And the last one is the ethylene signaling pathway. Taken together, our work identifies many genes responded to heat stress in maize kernels, and finds out seven genes and four pathways highly correlated with heat tolerance of maize kernels.

Functional characterization and allelic mining of OsGLR genes for potential uses in rice improvement

Glutamate-like receptor (GLR) genes are a group of regulatory genes involved in many physiological processes of plants. With 26 members in the rice genome, the functionalities of most rice GLR genes remain unknown. To facilitate their potential uses in rice improvement, an integrated strategy involving CRISPR-Cas9 mediated knockouts, deep mining and analyses of transcriptomic responses to different abiotic stresses/hormone treatments and gene CDS haplotype (gcHap) diversity in 3,010 rice genomes was taken to understand the functionalities of the 26 rice GLR genes, which led us to two conclusions. First, the expansion of rice GLR genes into a large gene family during evolution had gone through repeated gene duplication events occurred primarily in two large GLR gene clusters on rice chromosomes 9 and 6, which was accompanied with considerable functional differentiation. Secondly, except for two extremely conserved ones (OsGLR6.2 and OsGLR6.3), rich gcHap diversity exists at the remaining GLR genes which played important roles in rice population differentiation and rice improvement, evidenced by their very strong sub-specific and population differentiation, by their differentiated responses to day-length and different abiotic stresses, by the large phenotypic effects of five GLR gene knockout mutants on rice yield traits, by the significant association of major gcHaps at most GLR loci with yield traits, and by the strong genetic bottleneck effects and artificial selection on the gcHap diversity in populations Xian (indica) and Geng (japonica) during modern breeding. Our results suggest the potential values of the natural variation at most rice GLR loci for improving the productivity and tolerances to abiotic stresses. Additional efforts are needed to determine the phenotypic effects of major gcHaps at these GLR loci in order to identify 'favorable' alleles at specific GLR loci specific target traits in specific environments to facilitate their application to rice improvement in future.

Comprehensive Analysis of Glutamate Receptor-like Genes in Rice (Oryza sativa L.): Genome-Wide Identification, Characteristics, Evolution, Chromatin Accessibility, gcHap Diversity, Population Variation and Expression Analysis

Glutamate receptors (GLR) are widely present in animals and plants, playing essential roles in regulating plant growth, development and stress response. At present, most studies of GLRs in plants are focused on Arabidopsis thaliana, while there have been few studies on rice. In this study, we identified 26 OsGLR genes in rice (Oryza sativa L.). Then, we analyzed the chromosomal location, physical and chemical properties, subcellular location, transmembrane (TM) helices, signal peptides, three-dimensional (3D) structure, cis-acting elements, evolution, chromatin accessibility, population variation, gene-coding sequence haplotype (gcHap) and gene expression under multiple abiotic stress and hormone treatments. The results showed that out of the 26 OsGLR genes, ten genes had the TM domain, signal peptides and similar 3D structures. Most OsGLRs exhibited high tissue specificity in expression under drought stress. In addition, several OsGLR genes were specifically responsive to certain hormones. The favorable gcHap of many OsGLR genes in modern varieties showed obvious differentiation between Xian/indica and Geng/japonica subspecies. This study, for the first time, comprehensively analyzes the OsGLR genes in rice, and provides an important reference for further research on their molecular function.

Genetic and Molecular Characterization of a Self-Compatible Brassica rapa Line Possessing a New Class II S Haplotype

Most flowering plants have evolved a self-incompatibility (SI) system to maintain genetic diversity by preventing self-pollination. The Brassica species possesses sporophytic self-incompatibility (SSI), which is controlled by the pollen- and stigma-determinant factors SP11/SCR and SRK. However, the mysterious molecular mechanism of SI remains largely unknown. Here, a new class II S haplotype, named BrS-325, was identified in a pak choi line '325', which was responsible for the completely self-compatible phenotype. To obtain the entire S locus sequences, a complete pak choi genome was gained through Nanopore sequencing and de novo assembly, which provided a good reference genome for breeding and molecular research in B. rapa. S locus comparative analysis showed that the closest relatives to BrS-325 was BrS-60, and high sequence polymorphism existed in the S locus. Meanwhile, two duplicated SRKs (BrSRK-325a and BrSRK-325b) were distributed in the BrS-325 locus with opposite transcription directions. BrSRK-325b and BrSCR-325 were expressed normally at the transcriptional level. The multiple sequence alignment of SCRs and SRKs in class II S haplotypes showed that a number of amino acid variations were present in the contact regions (CR II and CR III) of BrSCR-325 and the hypervariable regions (HV I and HV II) of BrSRK-325s, which may influence the binding and interaction between the ligand and the receptor. Thus, these results suggested that amino acid variations in contact sites may lead to the SI destruction of a new class II S haplotype BrS-325 in B. rapa. The complete SC phenotype of '325' showed the potential for practical breeding application value in B. rapa.