Transcription factor bHLH121 regulates root cortical aerenchyma formation in maize.

Root anatomical phenotypes present a promising yet underexploited avenue to deliver major improvements in yield and climate resilience of crops by improving water and nutrient uptake. For instance, the formation of root cortical aerenchyma (RCA) significantly increases soil exploration and resource capture by reducing the metabolic costs of root tissue. A key bottleneck in studying such phenotypes has been the lack of robust high-throughput anatomical phenotyping platforms. We exploited a phenotyping approach based on laser ablation tomography, termed Anatomics, to quantify variation in RCA formation of 436 diverse maize lines in the field. Results revealed a significant and heritable variation for RCA formation. Genome-wide association studies identified a single-nucleotide polymorphism mapping to a root cortex-expressed gene-encoding transcription factor bHLH121. Functional studies identified that the bHLH121 Mu transposon mutant line and CRISPR/Cas9 loss-of-function mutant line showed reduced RCA formation, whereas an overexpression line exhibited significantly greater RCA formation when compared to the wild-type line. Characterization of these lines under suboptimal water and nitrogen availability in multiple soil environments revealed that bHLH121 is required for RCA formation developmentally as well as under studied abiotic stress. Overall functional validation of the bHLH121 gene's importance in RCA formation provides a functional marker to select varieties with improved soil exploration and thus yield under suboptimal conditions.

Genetic architecture of source-sink-regulated senescence in maize.

Source and sink interactions play a critical but mechanistically poorly understood role in the regulation of senescence. To disentangle the genetic and molecular mechanisms underlying source-sink-regulated senescence (SSRS), we performed a phenotypic, transcriptomic, and systems genetics analysis of senescence induced by the lack of a strong sink in maize (Zea mays). Comparative analysis of genotypes with contrasting SSRS phenotypes revealed that feedback inhibition of photosynthesis, a surge in reactive oxygen species, and the resulting endoplasmic reticulum (ER) stress were the earliest outcomes of weakened sink demand. Multienvironmental evaluation of a biparental population and a diversity panel identified 12 quantitative trait loci and 24 candidate genes, respectively, underlying SSRS. Combining the natural diversity and coexpression networks analyses identified 7 high-confidence candidate genes involved in proteolysis, photosynthesis, stress response, and protein folding. The role of a cathepsin B like protease 4 (ccp4), a candidate gene supported by systems genetic analysis, was validated by analysis of natural alleles in maize and heterologous analyses in Arabidopsis (Arabidopsis thaliana). Analysis of natural alleles suggested that a 700-bp polymorphic promoter region harboring multiple ABA-responsive elements is responsible for differential transcriptional regulation of ccp4 by ABA and the resulting variation in SSRS phenotype. We propose a model for SSRS wherein feedback inhibition of photosynthesis, ABA signaling, and oxidative stress converge to induce ER stress manifested as programed cell death and senescence. These findings provide a deeper understanding of signals emerging from loss of sink strength and offer opportunities to modify these signals to alter senescence program and enhance crop productivity.

A key to totipotency: Wuschel-like homeobox 2a unlocks embryogenic culture response in maize (Zea mays L.).

The ability of plant somatic cells to dedifferentiate, form somatic embryos and regenerate whole plants in vitro has been harnessed for both clonal propagation and as a key component of plant genetic engineering systems. Embryogenic culture response is significantly limited, however, by plant genotype in most species. This impedes advancements in both plant transformation-based functional genomics research and crop improvement efforts. We utilized natural variation among maize inbred lines to genetically map somatic embryo generation potential in tissue culture and identify candidate genes underlying totipotency. Using a series of maize lines derived from crosses involving the culturable parent A188 and the non-responsive parent B73, we identified a region on chromosome 3 associated with embryogenic culture response and focused on three candidate genes within the region based on genetic position and expression pattern. Two candidate genes showed no effect when ectopically expressed in B73, but the gene Wox2a was found to induce somatic embryogenesis and embryogenic callus proliferation. Transgenic B73 cells with strong constitutive expression of the B73 and A188 coding sequences of Wox2a were found to produce somatic embryos at similar frequencies, demonstrating that sufficient expression of either allele could rescue the embryogenic culture phenotype. Transgenic B73 plants were regenerated from the somatic embryos without chemical selection and no pleiotropic effects were observed in the Wox2a overexpression lines in the regenerated T0 plants or in the two independent events which produced T1 progeny. In addition to linking natural variation in tissue culture response to Wox2a, our data support the utility of Wox2a in enabling transformation of recalcitrant genotypes.

Effects of endosperm type and storage length of whole-plant corn silage on nitrogen fraction, fermentation products, zein profile, and starch digestibility.

Zeins are commercially important proteins found in corn endosperms. The objective of this study was to evaluate the effect of altering zein levels in corn inbred lines carrying endosperm mutations with differential allelic dosage and analyze the effects on the composition, nutritive value, and starch digestibility of whole-plant corn silage (WPCS) at 5 storage lengths. Three inbred lines carrying 3 different endosperm modifiers (opaque-2 [o2], floury-2 [fl2], and soft endosperm-1 [h1]) were pollinated with 2 pollen sources to form pairs of near-isogenic lines with either 2 or 3 doses of the mutant allele for each endosperm modifier. The experiment was designed as a split-plot design with 3 replications. Pollinated genotype was the main plot factor, and storage length was the subplot-level factor. Agronomic precautions were taken to mimic hybrid WPCS to the extent possible. Samples were collected at approximately 30% dry matter (DM) using a forage harvester and ensiled in heat-sealed plastic bags for 0, 30, 60, 120, and 240 d. Thus, the experiment consisted of 30 treatments (6 genotypes × 5 storage lengths) and 90 ensiling units (3 replications per treatment). Measurements included nutrient analysis, including crude protein, soluble crude protein, amylase-treated neutral detergent fiber, acid detergent fiber, lignin, starch, fermentation end products, zein concentration, and in vitro starch digestibility (ivSD). The nutritional profile of the inbred-based silage samples was similar to hybrid values reported in literature. Significant differences were found in fresh (unfermented) sample kernels for endosperm vitreousness and zein profiles between and within isogenic pairs. The o2 homozygous (3 doses of mutant allele) had the highest reduction in vitreousness level (74.5 to 38%) and zein concentration (6.2 to 4.7% of DM) compared with the heterozygous counterpart (2 doses of mutant allele). All genotypes showed significant reduction of total zeins and α-zeins during progressive storage length. In vitro starch digestibility increased with storage length and had significant effects of genotype and storage length but not for genotype by storage length interaction, which suggests that the storage period did not attenuate the difference in ivSD between near-isogenic pairs caused by zeins in WPCS. Both total zeins and α-zeins showed a strong negative correlation with ivSD, which agrees with the general hypothesis that the degradation of zeins increases ruminal starch degradability. Homozygous o2 was the only mutant with significantly higher ivSD compared with the heterozygous version, which suggests that, if all other conditions remain constant in a WPCS systems, substantial reductions in endosperm α-zeins are required to significantly improve ivSD in the silo.

Single-parent expression drives dynamic gene expression complementation in maize hybrids.

Single-parent expression (SPE) is defined as gene expression in only one of the two parents. SPE can arise from differential expression between parental alleles, termed non-presence/absence (non-PAV) SPE, or from the physical absence of a gene in one parent, termed PAV SPE. We used transcriptome data of diverse Zea mays (maize) inbreds and hybrids, including 401 samples from five different tissues, to test for differences between these types of SPE genes. Although commonly observed, SPE is highly genotype and tissue specific. A positive correlation was observed between the genetic distance of the two inbred parents and the number of SPE genes identified. Regulatory analysis showed that PAV SPE and non-PAV SPE genes are mainly regulated by cis effects, with a small fraction under trans regulation. Polymorphic transposable element insertions in promoter sequences contributed to the high level of cis regulation for PAV SPE and non-PAV SPE genes. PAV SPE genes were more frequently expressed in hybrids than non-PAV SPE genes. The expression of parentally silent alleles in hybrids of non-PAV SPE genes was relatively rare but occurred in most hybrids. Non-PAV SPE genes with expression of the silent allele in hybrids are more likely to exhibit above high parent expression level than hybrids that do not express the silent allele, leading to non-additive expression. This study provides a comprehensive understanding of the nature of non-PAV SPE and PAV SPE genes and their roles in gene expression complementation in maize hybrids.

Root angle in maize influences nitrogen capture and is regulated by calcineurin B-like protein (CBL)-interacting serine/threonine-protein kinase 15 (ZmCIPK15).

Crops with reduced nutrient and water requirements are urgently needed in global agriculture. Root growth angle plays an important role in nutrient and water acquisition. A maize diversity panel of 481 genotypes was screened for variation in root angle employing a high-throughput field phenotyping platform. Genome-wide association mapping identified several single nucleotide polymorphisms (SNPs) associated with root angle, including one located in the root expressed CBL-interacting serine/threonine-protein kinase 15 (ZmCIPK15) gene (LOC100285495). Reverse genetic studies validated the functional importance of ZmCIPK15, causing a approximately 10° change in root angle in specific nodal positions. A steeper root growth angle improved nitrogen capture in silico and in the field. OpenSimRoot simulations predicted at 40 days of growth that this change in angle would improve nitrogen uptake by 11% and plant biomass by 4% in low nitrogen conditions. In field studies under suboptimal N availability, the cipk15 mutant with steeper growth angles had 18% greater shoot biomass and 29% greater shoot nitrogen accumulation compared to the wild type after 70 days of growth. We propose that a steeper root growth angle modulated by ZmCIPK15 will facilitate efforts to develop new crop varieties with optimal root architecture for improved performance under edaphic stress.

Integrated Genome-Scale Analysis Identifies Novel Genes and Networks Underlying Senescence in Maize.

Premature senescence in annual crops reduces yield, while delayed senescence, termed stay-green, imposes positive and negative impacts on yield and nutrition quality. Despite its importance, scant information is available on the genetic architecture of senescence in maize (Zea mays) and other cereals. We combined a systematic characterization of natural diversity for senescence in maize and coexpression networks derived from transcriptome analysis of normally senescing and stay-green lines. Sixty-four candidate genes were identified by genome-wide association study (GWAS), and 14 of these genes are supported by additional evidence for involvement in senescence-related processes including proteolysis, sugar transport and signaling, and sink activity. Eight of the GWAS candidates, independently supported by a coexpression network underlying stay-green, include a trehalose-6-phosphate synthase, a NAC transcription factor, and two xylan biosynthetic enzymes. Source-sink communication and the activity of cell walls as a secondary sink emerge as key determinants of stay-green. Mutant analysis supports the role of a candidate encoding Cys protease in stay-green in Arabidopsis (Arabidopsis thaliana), and analysis of natural alleles suggests a similar role in maize. This study provides a foundation for enhanced understanding and manipulation of senescence for increasing carbon yield, nutritional quality, and stress tolerance of maize and other cereals.

Maize sugary enhancer1 (se1) is a gene affecting endosperm starch metabolism.

sugary enhancer1 (se1) is a naturally occurring mutant allele involved in starch metabolism in maize endosperm. It is a recessive modifier of sugary1 (su1) and commercially important in modern sweet corn breeding, but its molecular identity and mode of action remain unknown. Here, we developed a pair of near-isogenic lines, W822Gse (su1-ref/su1-ref se1/se1) and W822GSe (su1-ref/su1-ref Se1/Se1), that Mendelize the se1 phenotype in an su1-ref background. W822Gse kernels have lower starch and higher water soluble polysaccharide and sugars than W822GSe kernels. Using high-resolution genetic mapping, we found that wild-type Se1 is a gene Zm00001d007657 on chromosome 2 and a deletion of this gene causes the se1 phenotype. Comparative metabolic profiling of seed tissue between these 2 isolines revealed the remarkable difference in carbohydrate metabolism, with sucrose and maltose highly accumulated in the mutant. Se1 is predominantly expressed in the endosperm, with low expression in leaf and root tissues. Differential expression analysis identified genes enriched in both starch biosynthesis and degradation processes, indicating a pleiotropic regulatory effect of se1 Repressed expression of Se1 and Su1 in RNA interference-mediated transgenic maize validates that deletion of the gene identified as Se1 is a true causal gene responsible for the se1 phenotype. The findings contribute to our understanding of starch metabolism in cereal crops.

Variation and Inheritance of Small RNAs in Maize Inbreds and F1 Hybrids.

Small RNAs (sRNAs) regulate gene expression, play important roles in epigenetic pathways, and are hypothesized to contribute to hybrid vigor in plants. Prior investigations have provided valuable insights into associations between sRNAs and heterosis, often using a single hybrid genotype or tissue, but our understanding of the role of sRNAs and their potential value to plant breeding are limited by an incomplete picture of sRNA variation between diverse genotypes and development stages. Here, we provide a deep exploration of sRNA variation and inheritance among a panel of 108 maize (Zea mays) samples spanning five tissues from eight inbred parents and 12 hybrid genotypes, covering a spectrum of heterotic groups, genetic variation, and levels of heterosis for various traits. We document substantial developmental and genotypic influences on sRNA expression, with varying patterns for 21-nucleotide (nt), 22-nt, and 24-nt sRNAs. We provide a detailed view of the distribution of sRNAs in the maize genome, revealing a complex makeup that also shows developmental plasticity, particularly for 22-nt sRNAs. sRNAs exhibited substantially more variation between inbreds as compared with observed variation for gene expression. In hybrids, we identify locus-specific examples of nonadditive inheritance, mostly characterized as partial or complete dominance, but rarely outside the parental range. However, the global abundance of 21-nt, 22-nt, and 24-nt sRNAs varies very little between inbreds and hybrids, suggesting that hybridization affects sRNA expression principally at specific loci rather than on a global scale. This study provides a valuable resource for understanding the potential role of sRNAs in hybrid vigor.

Characterizing introgression-by-environment interactions using maize near isogenic lines.

KEY MESSAGE: Significant introgression-by-environment interactions are observed for traits throughout development from small introgressed segments of the genome. Relatively small genomic introgressions containing quantitative trait loci can have significant impacts on the phenotype of an individual plant. However, the magnitude of phenotypic effects for the same introgression can vary quite substantially in different environments due to introgression-by-environment interactions. To study potential patterns of introgression-by-environment interactions, fifteen near-isogenic lines (NILs) with > 90% B73 genetic background and multiple Mo17 introgressions were grown in 16 different environments. These environments included five geographical locations with multiple planting dates and multiple planting densities. The phenotypic impact of the introgressions was evaluated for up to 26 traits that span different growth stages in each environment to assess introgression-by-environment interactions. Results from this study showed that small portions of the genome can drive significant genotype-by-environment interaction across a wide range of vegetative and reproductive traits, and the magnitude of the introgression-by-environment interaction varies across traits. Some introgressed segments were more prone to introgression-by-environment interaction than others when evaluating the interaction on a whole plant basis throughout developmental time, indicating variation in phenotypic plasticity throughout the genome. Understanding the profile of introgression-by-environment interaction in NILs is useful in consideration of how small introgressions of QTL or transgene containing regions might be expected to impact traits in diverse environments.

Genome-wide association analysis of stalk biomass and anatomical traits in maize.

BACKGROUND: Maize stover is an important source of crop residues and a promising sustainable energy source in the United States. Stalk is the main component of stover, representing about half of stover dry weight. Characterization of genetic determinants of stalk traits provide a foundation to optimize maize stover as a biofuel feedstock. We investigated maize natural genetic variation in genome-wide association studies (GWAS) to detect candidate genes associated with traits related to stalk biomass (stalk diameter and plant height) and stalk anatomy (rind thickness, vascular bundle density and area). RESULTS: Using a panel of 942 diverse inbred lines, 899,784 RNA-Seq derived single nucleotide polymorphism (SNP) markers were identified. Stalk traits were measured on 800 members of the panel in replicated field trials across years. GWAS revealed 16 candidate genes associated with four stalk traits. Most of the detected candidate genes were involved in fundamental cellular functions, such as regulation of gene expression and cell cycle progression. Two of the regulatory genes (Zmm22 and an ortholog of Fpa) that were associated with plant height were previously shown to be involved in regulating the vegetative to floral transition. The association of Zmm22 with plant height was confirmed using a transgenic approach. Transgenic lines with increased expression of Zmm22 showed a significant decrease in plant height as well as tassel branch number, indicating a pleiotropic effect of Zmm22. CONCLUSION: Substantial heritable variation was observed in the association panel for stalk traits, indicating a large potential for improving useful stalk traits in breeding programs. Genome-wide association analyses detected several candidate genes associated with multiple traits, suggesting common regulatory elements underlie various stalk traits. Results of this study provide insights into the genetic control of maize stalk anatomy and biomass.

Draft Assembly of Elite Inbred Line PH207 Provides Insights into Genomic and Transcriptome Diversity in Maize.

Intense artificial selection over the last 100 years has produced elite maize (Zea mays) inbred lines that combine to produce high-yielding hybrids. To further our understanding of how genome and transcriptome variation contribute to the production of high-yielding hybrids, we generated a draft genome assembly of the inbred line PH207 to complement and compare with the existing B73 reference sequence. B73 is a founder of the Stiff Stalk germplasm pool, while PH207 is a founder of Iodent germplasm, both of which have contributed substantially to the production of temperate commercial maize and are combined to make heterotic hybrids. Comparison of these two assemblies revealed over 2500 genes present in only one of the two genotypes and 136 gene families that have undergone extensive expansion or contraction. Transcriptome profiling revealed extensive expression variation, with as many as 10,564 differentially expressed transcripts and 7128 transcripts expressed in only one of the two genotypes in a single tissue. Genotype-specific genes were more likely to have tissue/condition-specific expression and lower transcript abundance. The availability of a high-quality genome assembly for the elite maize inbred PH207 expands our knowledge of the breadth of natural genome and transcriptome variation in elite maize inbred lines across heterotic pools.

A robust, high-throughput method for computing maize ear, cob, and kernel attributes automatically from images.

Grain yield of the maize plant depends on the sizes, shapes, and numbers of ears and the kernels they bear. An automated pipeline that can measure these components of yield from easily-obtained digital images is needed to advance our understanding of this globally important crop. Here we present three custom algorithms designed to compute such yield components automatically from digital images acquired by a low-cost platform. One algorithm determines the average space each kernel occupies along the cob axis using a sliding-window Fourier transform analysis of image intensity features. A second counts individual kernels removed from ears, including those in clusters. A third measures each kernel's major and minor axis after a Bayesian analysis of contour points identifies the kernel tip. Dimensionless ear and kernel shape traits that may interrelate yield components are measured by principal components analysis of contour point sets. Increased objectivity and speed compared to typical manual methods are achieved without loss of accuracy as evidenced by high correlations with ground truth measurements and simulated data. Millimeter-scale differences among ear, cob, and kernel traits that ranged more than 2.5-fold across a diverse group of inbred maize lines were resolved. This system for measuring maize ear, cob, and kernel attributes is being used by multiple research groups as an automated Web service running on community high-throughput computing and distributed data storage infrastructure. Users may create their own workflow using the source code that is staged for download on a public repository.

Autophagic recycling plays a central role in maize nitrogen remobilization.

Autophagy is a primary route for nutrient recycling in plants by which superfluous or damaged cytoplasmic material and organelles are encapsulated and delivered to the vacuole for breakdown. Central to autophagy is a conjugation pathway that attaches AUTOPHAGY-RELATED8 (ATG8) to phosphatidylethanolamine, which then coats emerging autophagic membranes and helps with cargo recruitment, vesicle enclosure, and subsequent vesicle docking with the tonoplast. A key component in ATG8 function is ATG12, which promotes lipidation upon its attachment to ATG5. Here, we fully defined the maize (Zea mays) ATG system transcriptionally and characterized it genetically through atg12 mutants that block ATG8 modification. atg12 plants have compromised autophagic transport as determined by localization of a YFP-ATG8 reporter and its vacuolar cleavage during nitrogen or fixed-carbon starvation. Phenotypic analyses showed that atg12 plants are phenotypically normal and fertile when grown under nutrient-rich conditions. However, when nitrogen-starved, seedling growth is severely arrested, and as the plants mature, they show enhanced leaf senescence and stunted ear development. Nitrogen partitioning studies revealed that remobilization is impaired in atg12 plants, which significantly decreases seed yield and nitrogen-harvest index. Together, our studies demonstrate that autophagy, while nonessential, becomes critical during nitrogen stress and severely impacts maize productivity under suboptimal field conditions.

Genome-wide associations with flowering time in switchgrass using exome-capture sequencing data.

Flowering time is a major determinant of biomass yield in switchgrass (Panicum virgatum), a perennial bioenergy crop, because later flowering allows for an extended period of vegetative growth and increased biomass production. A better understanding of the genetic regulation of flowering time in switchgrass will aid the development of switchgrass varieties with increased biomass yields, particularly at northern latitudes, where late-flowering but southern-adapted varieties have high winter mortality. We use genotypes derived from recently published exome-capture sequencing, which mitigates challenges related to the large, highly repetitive and polyploid switchgrass genome, to perform genome-wide association studies (GWAS) using flowering time data from a switchgrass association panel in an effort to characterize the genetic architecture and genes underlying flowering time regulation in switchgrass. We identify associations with flowering time at multiple loci, including in a homolog of FLOWERING LOCUS T and in a locus containing TIMELESS, a homolog of a key circadian regulator in animals. Our results suggest that flowering time variation in switchgrass is due to variation at many positions across the genome. The relationship of flowering time and geographic origin indicates likely roles for genes in the photoperiod and autonomous pathways in generating switchgrass flowering time variation.

Insights into the maize pan-genome and pan-transcriptome.

Genomes at the species level are dynamic, with genes present in every individual (core) and genes in a subset of individuals (dispensable) that collectively constitute the pan-genome. Using transcriptome sequencing of seedling RNA from 503 maize (Zea mays) inbred lines to characterize the maize pan-genome, we identified 8681 representative transcript assemblies (RTAs) with 16.4% expressed in all lines and 82.7% expressed in subsets of the lines. Interestingly, with linkage disequilibrium mapping, 76.7% of the RTAs with at least one single nucleotide polymorphism (SNP) could be mapped to a single genetic position, distributed primarily throughout the nonpericentromeric portion of the genome. Stepwise iterative clustering of RTAs suggests, within the context of the genotypes used in this study, that the maize genome is restricted and further sampling of seedling RNA within this germplasm base will result in minimal discovery. Genome-wide association studies based on SNPs and transcript abundance in the pan-genome revealed loci associated with the timing of the juvenile-to-adult vegetative and vegetative-to-reproductive developmental transitions, two traits important for fitness and adaptation. This study revealed the dynamic nature of the maize pan-genome and demonstrated that a substantial portion of variation may lie outside the single reference genome for a species.

Genetic perturbation of the maize methylome.

DNA methylation can play important roles in the regulation of transposable elements and genes. A collection of mutant alleles for 11 maize (Zea mays) genes predicted to play roles in controlling DNA methylation were isolated through forward- or reverse-genetic approaches. Low-coverage whole-genome bisulfite sequencing and high-coverage sequence-capture bisulfite sequencing were applied to mutant lines to determine context- and locus-specific effects of these mutations on DNA methylation profiles. Plants containing mutant alleles for components of the RNA-directed DNA methylation pathway exhibit loss of CHH methylation at many loci as well as CG and CHG methylation at a small number of loci. Plants containing loss-of-function alleles for chromomethylase (CMT) genes exhibit strong genome-wide reductions in CHG methylation and some locus-specific loss of CHH methylation. In an attempt to identify stocks with stronger reductions in DNA methylation levels than provided by single gene mutations, we performed crosses to create double mutants for the maize CMT3 orthologs, Zmet2 and Zmet5, and for the maize DDM1 orthologs, Chr101 and Chr106. While loss-of-function alleles are viable as single gene mutants, the double mutants were not recovered, suggesting that severe perturbations of the maize methylome may have stronger deleterious phenotypic effects than in Arabidopsis thaliana.

Diversity and population structure of northern switchgrass as revealed through exome capture sequencing.

Panicum virgatum L. (switchgrass) is a polyploid, perennial grass species that is native to North America, and is being developed as a future biofuel feedstock crop. Switchgrass is present primarily in two ecotypes: a northern upland ecotype, composed of tetraploid and octoploid accessions, and a southern lowland ecotype, composed of primarily tetraploid accessions. We employed high-coverage exome capture sequencing (~2.4 Tb) to genotype 537 individuals from 45 upland and 21 lowland populations. From these data, we identified ~27 million single-nucleotide polymorphisms (SNPs), of which 1 590 653 high-confidence SNPs were used in downstream analyses of diversity within and between the populations. From the 66 populations, we identified five primary population groups within the upland and lowland ecotypes, a result that was further supported through genetic distance analysis. We identified conserved, ecotype-restricted, non-synonymous SNPs that are predicted to affect the protein function of CONSTANS (CO) and EARLY HEADING DATE 1 (EHD1), key genes involved in flowering, which may contribute to the phenotypic differences between the two ecotypes. We also identified, relative to the near-reference Kanlow population, 17 228 genes present in more copies than in the reference genome (up-CNVs), 112 630 genes present in fewer copies than in the reference genome (down-CNVs) and 14 430 presence/absence variants (PAVs), affecting a total of 9979 genes, including two upland-specific CNV clusters. In total, 45 719 genes were affected by an SNP, CNV, or PAV across the panel, providing a firm foundation to identify functional variation associated with phenotypic traits of interest for biofuel feedstock production.

Evidence for maternal control of seed size in maize from phenotypic and transcriptional analysis.

Seed size is an important component of grain yield and a key determinant trait for crop domestication. The Krug Yellow Dent long-term selection experiment for large and small seed provides a valuable resource to dissect genetic and phenotypic changes affecting seed size within a common genetic background. In this study, inbred lines derived from Krug Large Seed (KLS) and Krug Small Seed (KSS) populations and reciprocal F1 crosses were used to investigate developmental and molecular mechanisms governing seed size. Seed morphological characteristics showed striking differences between KLS and KSS inbred lines, and the reciprocal cross experiment revealed a strong maternal influence on both seed weight and seed size. Quantification of endosperm area, starchy endosperm cell size, and kernel dry mass accumulation indicated a positive correlation between seed size, endosperm cell number, and grain filling rate, and patterns of grain filling in reciprocal crosses mirrored that of the maternal parent. Consistent with the maternal contribution to seed weight, transcriptome profiling of reciprocal F1 hybrids showed substantial similarities to the maternal parent. A set of differentially expressed genes between KLS and KSS inbreds were found, which fell into a broad number of functional categories including DNA methylation, nucleosome assembly, and heat stress response. In addition, gene co-expression network analysis of parental inbreds and reciprocal F1 hybrids identified co-expression modules enriched in ovule development and DNA methylation, implicating these two processes in seed size determination. These results expand our understanding of seed size regulation and help to uncover the developmental and molecular basis underlying maternal control of seed size in maize.