Application of CRISPR/Cas12i.3 for targeted mutagenesis in broomcorn millet (Panicum miliaceum L.).

In the present study, we have successfully established a gene editing platform in broomcorn millet, one of the oldest crops originating from China, by using our CRISPR/Cas12i.3, and we also created new elite germplasm for this crop.

Genomic variation in weedy and cultivated broomcorn millet accessions uncovers the genetic architecture of agronomic traits.

Large-scale genomic variations are fundamental resources for crop genetics and breeding. Here we sequenced 1,904 genomes of broomcorn millet to an average of 40× sequencing depth and constructed a comprehensive variation map of weedy and cultivated accessions. Being one of the oldest cultivated crops, broomcorn millet has extremely low nucleotide diversity and remarkably rapid decay of linkage disequilibrium. Genome-wide association studies identified 186 loci for 12 agronomic traits. Many causative candidate genes, such as PmGW8 for grain size and PmLG1 for panicle shape, showed strong selection signatures during domestication. Weedy accessions contained many beneficial variations for the grain traits that are largely lost in cultivated accessions. Weedy and cultivated broomcorn millet have adopted different loci controlling flowering time for regional adaptation in parallel. Our study uncovers the unique population genomic features of broomcorn millet and provides an agronomically important resource for cereal crops.

ZmELP1, an Elongator complex subunit, is required for the maintenance of histone acetylation and RNA Pol II phosphorylation in maize kernels.

The Elongator complex was originally identified as an interactor of hyperphosphorylated RNA polymerase II (RNAPII) in yeast and has histone acetyltransferase (HAT) activity. However, the genome-wide regulatory roles of Elongator on transcriptional elongation and histone acetylation remain unclear. We characterized a maize miniature seed mutant, mn7 and map-based cloning revealed that Mn7 encodes one of the subunits of the Elongator complex, ZmELP1. ZmELP1 deficiency causes marked reductions in the kernel size and weight. Molecular analyses showed that ZmELP1 interacts with ZmELP3, which is required for H3K14 acetylation (H3K14ac), and Elongator complex subunits interact with RNA polymerase II (RNAPII) C-terminal domain (CTD). Genome-wide analyses indicated that loss of ZmELP1 leads to a significant decrease in the deposition of H3K14ac and the CTD of phosphorylated RNAPII on Ser2 (Ser2P). These chromatin changes positively correlate with global transcriptomic changes. ZmELP1 mutation alters the expression of genes involved in transcriptional regulation and kernel development. We also showed that the decrease of Ser2P depends on the deposition of Elongator complex-mediated H3K14ac. Taken together, our results reveal an important role of ZmELP1 in the H3K14ac-dependent transcriptional elongation, which is critical for kernel development.

Enhancing maize's nitrogen-fixing potential through ZmSBT3, a gene suppressing mucilage secretion.

Maize (Zea mays) requires substantial amounts of nitrogen, posing a challenge for its cultivation. Recent work discovered that some ancient Mexican maize landraces harbored diazotrophic bacteria in mucilage secreted by their aerial roots. To see if this trait is retained in modern maize, we conducted a field study of aerial root mucilage (ARM) in 258 inbred lines. We observed that ARM secretion is common in modern maize, but the amount significantly varies, and only a few lines have retained the nitrogen-fixing traits found in ancient landraces. The mucilage of the high-ARM inbred line HN5-724 had high nitrogen-fixing enzyme activity and abundant diazotrophic bacteria. Our genome-wide association study identified 17 candidate genes associated with ARM across three environments. Knockouts of one candidate gene, the subtilase family gene ZmSBT3, confirmed that it negatively regulates ARM secretion. Notably, the ZmSBT3 knockout lines had increased biomass and total nitrogen accumulation under nitrogen-free culture conditions. High ARM was associated with three ZmSBT3 haplotypes that were gradually lost during maize domestication, being retained in only a few modern inbred lines such as HN5-724. In summary, our results identify ZmSBT3 as a potential tool for enhancing ARM, and thus nitrogen fixation, in maize.

Dynamic transcriptome landscape of developing maize ear.

Seed number and harvesting ability in maize (Zea mays L.) are primarily determined by the architecture of female inflorescence, namely the ear. Therefore, ear morphogenesis contributes to grain yield and as such is one of the key target traits during maize breeding. However, the molecular networks of this highly dynamic and complex grain-bearing inflorescence remain largely unclear. As a first step toward characterizing these networks, we performed a high-spatio-temporal-resolution investigation of transcriptomes using 130 ear samples collected from developing ears with length from 0.1 mm to 19.0 cm. Comparisons of these mRNA populations indicated that these spatio-temporal transcriptomes were clearly separated into four distinct stages stages I, II, III, and IV. A total of 23 793 genes including 1513 transcription factors (TFs) were identified in the investigated developing ears. During the stage I of ear morphogenesis, 425 genes were predicted to be involved in a co-expression network established by eight hub TFs. Moreover, 9714 ear-specific genes were identified in the seven kinds of meristems. Additionally, 527 genes including 59 TFs were identified as especially expressed in ear and displayed high temporal specificity. These results provide a high-resolution atlas of gene activity during ear development and help to unravel the regulatory modules associated with the differentiation of the ear in maize.

A transcription factor ZmGLK36 confers broad resistance to maize rough dwarf disease in cereal crops.

Maize rough dwarf disease (MRDD), caused by maize rough dwarf virus (MRDV) or rice black-streaked dwarf virus (RBSDV), seriously threatens worldwide production of all major cereal crops, including maize, rice, wheat and barley. Here we report fine mapping and cloning of a previously reported major quantitative trait locus (QTL) (qMrdd2) for RBSDV resistance in maize. Subsequently, we show that qMrdd2 encodes a G2-like transcription factor named ZmGLK36 that promotes resistance to RBSDV by enhancing jasmonic acid (JA) biosynthesis and JA-mediated defence response. We identify a 26-bp indel located in the 5' UTR of ZmGLK36 that contributes to differential expression and resistance to RBSDV in maize inbred lines. Moreover, we show that ZmDBF2, an AP2/EREBP family transcription factor, directly binds to the 26-bp indel and represses ZmGLK36 expression. We further demonstrate that ZmGLK36 plays a conserved role in conferring resistance to RBSDV in rice and wheat using transgenic or marker-assisted breeding approaches. Our results provide insights into the molecular mechanisms of RBSDV resistance and effective strategies to breed RBSDV-resistant cereal crops.

Targeted large fragment deletion in plants using paired crRNAs with type I CRISPR system.

The CRISPR-Cas systems have been widely used as genome editing tools, with type II and V systems typically introducing small indels, and type I system mediating long-range deletions. However, the precision of type I systems for large fragment deletion is still remained to be optimized. Here, we developed a compact Cascade-Cas3 Dvu I-C system with Cas11c for plant genome editing. The Dvu I-C system was efficient to introduce controllable large fragment deletion up to at least 20 kb using paired crRNAs. The paired-crRNAs design also improved the controllability of deletions for the type I-E system. Dvu I-C system was sensitive to spacer length and mismatch, which was benefit for target specificity. In addition, we showed that the Dvu I-C system was efficient for generating stable transgenic lines in maize and rice with the editing efficiency up to 86.67%. Overall, Dvu I-C system we developed here is powerful for achieving controllable large fragment deletions.

A complete telomere-to-telomere assembly of the maize genome.

A complete telomere-to-telomere (T2T) finished genome has been the long pursuit of genomic research. Through generating deep coverage ultralong Oxford Nanopore Technology (ONT) and PacBio HiFi reads, we report here a complete genome assembly of maize with each chromosome entirely traversed in a single contig. The 2,178.6 Mb T2T Mo17 genome with a base accuracy of over 99.99% unveiled the structural features of all repetitive regions of the genome. There were several super-long simple-sequence-repeat arrays having consecutive thymine-adenine-guanine (TAG) tri-nucleotide repeats up to 235 kb. The assembly of the entire nucleolar organizer region of the 26.8 Mb array with 2,974 45S rDNA copies revealed the enormously complex patterns of rDNA duplications and transposon insertions. Additionally, complete assemblies of all ten centromeres enabled us to precisely dissect the repeat compositions of both CentC-rich and CentC-poor centromeres. The complete Mo17 genome represents a major step forward in understanding the complexity of the highly recalcitrant repetitive regions of higher plant genomes.

QTL Analysis Reveals Conserved and Differential Genetic Regulation of Maize Lateral Angles above the Ear.

Improving the density tolerance and planting density has great importance for increasing maize production. The key to promoting high density planting is breeding maize with a compact canopy architecture, which is mainly influenced by the angles of the leaves and tassel branches above the ear. It is still unclear whether the leaf angles of different stem nodes and tassel branches are controlled by similar genetic regulatory mechanisms, which limits the ability to breed for density-tolerant maize. Here, we developed a population with 571 double haploid lines derived from inbred lines, PHBA6 and Chang7-2, showing significant differences in canopy architecture. Phenotypic and QTL analyses revealed that the genetic regulation mechanism was largely similar for closely adjacent leaves above the ears. In contrast, the regulation mechanisms specifying the angles of distant leaves and the angles of leaves vs. tassel branches are largely different. The liguless1 gene was identified as a candidate gene for QTLs co-regulating the angles of different leaves and the tassel branch, consistent with its known roles in regulating plant architecture. Our findings can be used to develop strategies for the improvement of leaf and tassel architecture through the introduction of trait-specific or pleiotropic genes, thus benefiting the breeding of maize with increased density tolerance in the future.

HEAT SHOCK PROTEIN 90.6 interacts with carbon and nitrogen metabolism components during seed development.

Carbon and nitrogen are the two main nutrients in maize (Zea mays L.) kernels, and kernel filling and metabolism determine seed formation and germination. However, the molecular mechanisms underlying the relationship between kernel filling and corresponding carbon and nitrogen metabolism remain largely unknown. Here, we found that HEAT SHOCK PROTEIN 90.6 (HSP90.6) is involved in both seed filling and the metabolism processes of carbon and nitrogen. A single-amino acid mutation within the HATPase_c domain of HSP90.6 led to small kernels. Transcriptome profiling showed that the expression of amino acid biosynthesis- and carbon metabolism-related genes was significantly downregulated in the hsp90.6 mutant. Further molecular evidence showed strong interactions between HSP90.6 and the 26S proteasome subunits REGULATORY PARTICLE NON-ATPASE6 (RPN6) and PROTEASOME BETA SUBUNITD2 (PBD2). The mutation of hsp90.6 significantly reduced the activity of the 26S proteasome, resulting in the accumulation of ubiquitinated proteins and defects in nitrogen recycling. Moreover, we verified that HSP90.6 is involved in carbon metabolism through interacting with the 14-3-3 protein GENERAL REGULATORY FACTOR14-4 (GF14-4). Collectively, our findings revealed that HSP90.6 is involved in seed filling and development by interacting with the components controlling carbon and nitrogen metabolism.

A teosinte-derived allele of an HKT1 family sodium transporter improves salt tolerance in maize.

The sodium cation (Na+ ) is the predominant cation with deleterious effects on crops in salt-affected agricultural areas. Salt tolerance of crop can be improved by increasing shoot Na+ exclusion. Therefore, it is crucial to identify and use genetic variants of various crops that promote shoot Na+ exclusion. Here, we show that a HKT1 family gene ZmNC3 (Zea mays L. Na+ Content 3; designated ZmHKT1;2) confers natural variability in shoot-Na+ accumulation and salt tolerance in maize. ZmHKT1;2 encodes a Na+ -preferential transporter localized in the plasma membrane, which mediates shoot Na+ exclusion, likely by withdrawing Na+ from the root xylem flow. A naturally occurring nonsynonymous SNP (SNP947-G) increases the Na+ transport activity of ZmHKT1;2, promoting shoot Na+ exclusion and salt tolerance in maize. SNP947-G first occurred in the wild grass teosinte (at a allele frequency of 43%) and has become a minor allele in the maize population (allele frequency 6.1%), suggesting that SNP947-G is derived from teosinte and that the genomic region flanking SNP947 likely has undergone selection during domestication or post-domestication dispersal of maize. Moreover, we demonstrate that introgression of the SNP947-G ZmHKT1;2 allele into elite maize germplasms reduces shoot Na+ content by up to 80% and promotes salt tolerance. Taken together, ZmNC3/ZmHKT1;2 was identified as an important QTL promoting shoot Na+ exclusion, and its favourable allele provides an effective tool for developing salt-tolerant maize varieties.

The sugar transporter ZmSUGCAR1 of the nitrate transporter 1/peptide transporter family is critical for maize grain filling.

Maternal-to-filial nutrition transfer is central to grain development and yield. nitrate transporter 1/peptide transporter (NRT1-PTR)-type transporters typically transport nitrate, peptides, and ions. Here, we report the identification of a maize (Zea mays) NRT1-PTR-type transporter that transports sucrose and glucose. The activity of this sugar transporter, named Sucrose and Glucose Carrier 1 (SUGCAR1), was systematically verified by tracer-labeled sugar uptake and serial electrophysiological studies including two-electrode voltage-clamp, non-invasive microelectrode ion flux estimation assays in Xenopus laevis oocytes and patch clamping in HEK293T cells. ZmSUGCAR1 is specifically expressed in the basal endosperm transfer layer and loss-of-function mutation of ZmSUGCAR1 caused significantly decreased sucrose and glucose contents and subsequent shrinkage of maize kernels. Notably, the ZmSUGCAR1 orthologs SbSUGCAR1 (from Sorghum bicolor) and TaSUGCAR1 (from Triticum aestivum) displayed similar sugar transport activities in oocytes, supporting the functional conservation of SUGCAR1 in closely related cereal species. Thus, the discovery of ZmSUGCAR1 uncovers a type of sugar transporter essential for grain development and opens potential avenues for genetic improvement of seed-filling and yield in maize and other grain crops.

Maize sterility gene DRP1 encodes a desiccation-related protein that is critical for Ubisch bodies and pollen exine development.

Desiccation tolerance is a remarkable feature of pollen, seeds, and resurrection-type plants. Exposure to desiccation stress can cause sporophytic defects, resulting in male sterility. Here, we report the novel maize sterility gene DRP1 (Desiccation-Related Protein 1), which was identified by bulked-segregant analysis sequencing and encodes a desiccation-related protein. Loss of function of DRP1 results in abnormal Ubisch bodies, defective tectum of the pollen exine, and complete male sterility. Our results suggest that DRP1 may facilitate anther dehydration to maintain appropriate water status. DRP1 is a secretory protein that is specifically expressed in the tapetum and microspore from the tetrad to the uninucleate microspore stage. Differentially expressed genes in drp1 are enriched in Gene Ontology terms for pollen exine formation, polysaccharide catabolic process, extracellular region, and response to heat. In addition, DRP1 is a target of selection that appears to have played an important role in the spread of maize from tropical/subtropical to temperate regions. Taken together, our results suggest that DRP1 encodes a desiccation-related protein whose loss of function causes male sterility. Our findings provide a potential genetic resource that may be used to design crops for heterosis utilization.

Identification of two new QTLs of maize (Zea mays L.) underlying kernel row number using the HNAU-NAM1 population.

BACKGROUND: Maize kernel row number (KRN) is one of the most important yield traits and has changed greatly during maize domestication and selection. Elucidating the genetic basis of KRN will be helpful to improve grain yield in maize. RESULTS: Here, we measured KRN in four environments using a nested association mapping (NAM) population named HNAU-NAM1 with 1,617 recombinant inbred lines (RILs) that were derived from 12 maize inbred lines with a common parent, GEMS41. Then, five consensus quantitative trait loci (QTLs) distributing on four chromosomes were identified in at least three environments along with the best linear unbiased prediction (BLUP) values by the joint linkage mapping (JLM) method. These QTLs were further validated by the separate linkage mapping (SLM) and genome-wide association study (GWAS) methods. Three KRN genes cloned through the QTL assay were found in three of the five consensus QTLs, including qKRN1.1, qKRN2.1 and qKRN4.1. Two new QTLs of KRN, qKRN4.2 and qKRN9.1, were also identified. On the basis of public RNA-seq and genome annotation data, five genes highly expressed in ear tissue were considered candidate genes contributing to KRN. CONCLUSIONS: This study carried out a comprehensive analysis of the genetic architecture of KRN by using a new NAM population under multiple environments. The present results provide solid information for understanding the genetic components underlying KRN and candidate genes in qKRN4.2 and qKRN9.1. Single-nucleotide polymorphisms (SNPs) closely linked to qKRN4.2 and qKRN9.1 could be used to improve inbred yield during molecular breeding in maize.

Maize PPR278 Functions in Mitochondrial RNA Splicing and Editing.

Pentatricopeptide repeat (PPR) proteins are a large protein family in higher plants and play important roles during seed development. Most reported PPR proteins function in mitochondria. However, some PPR proteins localize to more than one organelle; functional characterization of these proteins remains limited in maize (Zea mays L.). Here, we cloned and analyzed the function of a P-subfamily PPR protein, PPR278. Loss-function of PPR278 led to a lower germination rate and other defects at the seedling stage, as well as smaller kernels compared to the wild type. PPR278 was expressed in all investigated tissues. Furthermore, we determined that PPR278 is involved in the splicing of two mitochondrial transcripts (nad2 intron 4 and nad5 introns 1 and 4), as well as RNA editing of C-to-U sites in 10 mitochondrial transcripts. PPR278 localized to the nucleus, implying that it may function as a transcriptional regulator during seed development. Our data indicate that PPR278 is involved in maize seed development via intron splicing and RNA editing in mitochondria and has potential regulatory roles in the nucleus.

Geometry of sequence working memory in macaque prefrontal cortex.

How the brain stores a sequence in memory remains largely unknown. We investigated the neural code underlying sequence working memory using two-photon calcium imaging to record thousands of neurons in the prefrontal cortex of macaque monkeys memorizing and then reproducing a sequence of locations after a delay. We discovered a regular geometrical organization: The high-dimensional neural state space during the delay could be decomposed into a sum of low-dimensional subspaces, each storing the spatial location at a given ordinal rank, which could be generalized to novel sequences and explain monkey behavior. The rank subspaces were distributed across large overlapping neural groups, and the integration of ordinal and spatial information occurred at the collective level rather than within single neurons. Thus, a simple representational geometry underlies sequence working memory.

The RppC-AvrRppC NLR-effector interaction mediates the resistance to southern corn rust in maize.

Southern corn rust (SCR), caused by the fungal pathogen Puccinia polysora, is a major threat to maize production worldwide. Efficient breeding and deployment of resistant hybrids are key to achieving durable control of SCR. Here, we report the molecular cloning and characterization of RppC, which encodes an NLR-type immune receptor and is responsible for a major SCR resistance quantitative trait locus. Furthermore, we identified the corresponding avirulence effector, AvrRppC, which is secreted by P. polysora and triggers RppC-mediated resistance. Allelic variation of AvrRppC directly determines the effectiveness of RppC-mediated resistance, indicating that monitoring of AvrRppC variants in the field can guide the rational deployment of RppC-containing hybrids in maize production. Currently, RppC is the most frequently deployed SCR resistance gene in China, and a better understanding of its mode of action is critical for extending its durability.

Genetic dissection of maize plant architecture using a novel nested association mapping population.

The leaf angle (LA), plant height (PH), and ear height (EH) are key plant architectural traits influencing maize (Zea mays L.) yield. However, their genetic determinants have not yet been well-characterized. Here, we developed a maize advanced backcross-nested association mapping population in Henan Agricultural University (HNAU-NAM1) comprised of 1,625 BC1 F4 /BC2 F4 lines. These were obtained by crossing a diverse set of 12 representative inbred lines with the common GEMS41 line, which were then genotyped using the MaizeSNP9.4K array. Genetic diversity and phenotypic distribution analyses showed considerable levels of genetic variation. We obtained 18-88 quantitative trait loci (QTLs) associated with LA, PH, and EH by using three complementary mapping methods, named as separate linkage mapping, joint linkage mapping, and genome-wide association studies. Our analyses enabled the identification of ten QTL hot-spot regions associated with the three traits, which were distributed on nine different chromosomes. We further selected 13 major QTLs that were simultaneously detected by three methods and deduced the candidate genes, of which eight were not reported before. The newly constructed HNAU-NAM1 population in this study will further broaden our insights into understanding of genetic regulation of plant architecture, thus will help to improve maize yield and provide an invaluable resource for maize functional genomics and breeding research.

Loss-of-function alleles of ZmPLD3 cause haploid induction in maize.

Doubled haploid technology has been widely applied to multiple plant species and is recognized as one of the most important technologies for improving crop breeding efficiency. Although mutations in MATRILINEAL/Zea mays PHOSPHOLIPASE A1/NOT LIKE DAD (MTL/ZmPLA1/NLD) and Zea mays DOMAIN OF UNKNOWN FUNCTION 679 MEMBRANE PROTEIN (ZmDMP) have been shown to generate haploids in maize, knowledge of the genetic basis of haploid induction (HI) remains incomplete. Therefore, cloning of new genes underlying HI is important for further elucidating its genetic architecture. Here, we found that loss-of-function mutations of Zea mays PHOSPHOLIPASE D3 (ZmPLD3), one of the members from the phospholipase D subfamily, could trigger maternal HI in maize. ZmPLD3 was identified through a reverse genetic strategy based on analysis of pollen-specifically expressed phospholipases, followed by validation through the clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR-Cas9) system. Mutations of ZmPLD3 resulted in a haploid induction rate (HIR) similar to that of mtl/zmpla1/nld and showed synergistic effects rather than functional redundancy on tripling the HIR (from 1.19% to 4.13%) in the presence of mtl/zmpla1/nld. RNA-seq profiling of mature pollen indicated that a large number of pollen-specific differentially expressed genes were enriched in processes related to gametogenesis development, such as pollen tube development and cell communication, during the double-fertilization process. In addition, ZmPLD3 is highly conserved among cereals, highlighting the potential application of these in vivo haploid-inducer lines for other important crop plant species. Collectively, our discovery identifies a novel gene underlying in vivo maternal HI and provides possibility of breeding haploid inducers with further improved HIR.

ZmCTLP1 is required for the maintenance of lipid homeostasis and the basal endosperm transfer layer in maize kernels.

Maize kernel weight is influenced by the unloading of nutrients from the maternal placenta and their passage through the transfer tissue of the basal endosperm transfer layer (BETL) and the basal intermediate zone (BIZ) to the upper part of the endosperm. Here, we show that Small kernel 10 (Smk10) encodes a choline transporter-like protein 1 (ZmCTLP1) that facilitates choline uptake and is located in the trans-Golgi network (TGN). Its loss of function results in reduced choline content, leading to smaller kernels with a lower starch content. Mutation of ZmCTLP1 disrupts membrane lipid homeostasis and the normal development of wall in-growths. Expression levels of Mn1 and ZmSWEET4c, two kernel filling-related genes, are downregulated in the smk10, which is likely to be one of the major causes of incompletely differentiated transfer cells. Mutation of ZmCTLP1 also reduces the number of plasmodesmata (PD) in transfer cells, indicating that the smk10 mutant is impaired in PD formation. Intriguingly, we also observed premature cell death in the BETL and BIZ of the smk10 mutant. Together, our results suggest that ZmCTLP1-mediated choline transport affects kernel development, highlighting its important role in lipid homeostasis, wall in-growth formation and PD development in transfer cells.