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.

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.

Genetic Architecture of Maize Rind Strength Revealed by the Analysis of Divergently Selected Populations.

The strength of the stalk rind, measured as rind penetrometer resistance (RPR), is an important contributor to stalk lodging resistance. To enhance the genetic architecture of RPR, we combined selection mapping on populations developed by 15 cycles of divergent selection for high and low RPR with time-course transcriptomic and metabolic analyses of the stalks. Divergent selection significantly altered allele frequencies of 3,656 and 3,412 single- nucleotide polymorphisms (SNPs) in the high and low RPR populations, respectively. Surprisingly, only 110 (1.56%) SNPs under selection were common in both populations, while the majority (98.4%) were unique to each population. This result indicated that high and low RPR phenotypes are produced by biologically distinct mechanisms. Remarkably, regions harboring lignin and polysaccharide genes were preferentially selected in high and low RPR populations, respectively. The preferential selection was manifested as higher lignification and increased saccharification of the high and low RPR stalks, respectively. The evolution of distinct gene classes according to the direction of selection was unexpected in the context of parallel evolution and demonstrated that selection for a trait, albeit in different directions, does not necessarily act on the same genes. Tricin, a grass-specific monolignol that initiates the incorporation of lignin in the cell walls, emerged as a key determinant of RPR. Integration of selection mapping and transcriptomic analyses with published genetic studies of RPR identified several candidate genes including ZmMYB31, ZmNAC25, ZmMADS1, ZmEXPA2, ZmIAA41 and hk5. These findings provide a foundation for an enhanced understanding of RPR and the improvement of stalk lodging resistance.

Genome-wide identification, expression profiling, and network analysis of AT-hook gene family in maize.

AT-hook motif nuclear localized (AHL) genes have diverse but poorly understood biological functions. We identified and analyzed 37 AHL genes in maize. We also discovered four and one additional AHLs in rice and sorghum, respectively, besides those reported earlier. The maize AHLs were classified into two clades (A and B) and three distinct types (I, II, and III) as also reported in Arabidopsis. Phylogenetic and ortholog analyses showed that, while the evolutionary classification was conserved in plants, expansion of the AHL gene family in maize was accompanied with new biological functions. Gene structure analysis showed that, while all but one Type-I AHLs lacked an intron, origin of Type-II and Type-III AHLs was associated with the gain of introns suggesting evolutionarily distinct temporal and spatial expression patterns and, likely, neofunctionalization. Gene duplication analysis revealed that AHLs in maize expanded via dispersive duplication further supporting their functional diversity. To discern these functions, we analyzed 71 transcriptomes from diverse tissues and developmental stages of maize and classified AHLs into eight groups with distinct temporal/spatial expression profiles. Coexpression analysis implicated 5 AHLs and 33 novel genes in networks specific to endosperm, seed, root, leaf, and reproductive tissues indicating their role in the development of these organs. Major processes coregulated by AHLs include pollen development, drought response, senescence, and wound response. We also identified interactions of AHL proteins in coregulating important processes including stress response. These novel insights into the role of AHLs in plant development provide a platform for functional analyses in maize and related grasses.

Characterization of natural genetic variation identifies multiple genes involved in salt tolerance in maize.

Progressive decline in irrigation water is forcing farmers to use brackish water which increases soil salinity and exposes the crop plants to salinity. Maize, one of the most important crops, is sensitive to salinity. Salt tolerance is a complex trait controlled by a number of physiological and biochemical processes. Scant information is available on the genetic architecture of salt tolerance in maize. We evaluated 399 inbred lines for six early vigor shoot and root traits upon exposure of 18-day seedlings to salinity (ECiw = 16 dS m-1) stress. Contrasting response of shoot and root growth to salinity indicated a meticulous reprogramming of resource partitioning by the plants to cope with the stress. The genomic analysis identified 57 single nucleotide polymorphisms (SNP) associated with early vigor traits. Candidate genes systematically associated with each SNP include both previously known and novel genes. Important candidates include a late embryogenesis protein, a divalent ion symporter, a proton extrusion protein, an RNA-binding protein, a casein kinase 1, and an AP2/EREBP transcription factor. The late embryogenesis protein is associated with both shoot and root length, indicating a coordinated change in resource allocation upon salt stress. Identification of a casein kinase 1 indicates an important role for Ser/Thr kinases in salt tolerance. Validation of eight candidates based on expression in a salt-tolerant and a salt-sensitive inbred line supported their role in salt tolerance. The candidate genes identified in this investigation provide a foundation for dissecting genetic and molecular regulation of salt tolerance in maize and related grasses.

Sugar partitioning and source-sink interaction are key determinants of leaf senescence in maize.

Source-sink communication is one of the key regulators of senescence; however, the mechanisms underlying such regulation are largely unknown. We analysed senescence induced by the lack of grain sink in maize, termed source-sink regulated senescence (SSRS), and compared the associated physiological and metabolic changes with those accompanying natural senescence. Phenotypic characterization of 31 diverse field-grown inbreds revealed substantial variation for both SSRS and natural senescence. Partitioning of excess carbohydrates to alternative sinks, mainly internodes and husks, emerged as a critical mechanism underlying both SSRS and stay-green. Time-course analyses of SSRS sensitive (B73) and resistant (PHG35) inbreds confirmed the role of sugar partitioning in SSRS and stay-green. Elevated hemicellulose content in PHG35 internodes highlighted the role of the cell wall as a significant alternative sink. Sugar signalling emerged as an important regulator of SSRS as evident from an increased accumulation of trehalose-6-phosphate and decreased transcript levels of snf1-related protein kinase1, two signalling components associated with senescence, in B73. These findings demonstrate a crucial role of sugar partitioning, signalling, and utilization in SSRS. Available genetic variation for SSRS and a better understanding of the underlying mechanisms would help modify sugar partitioning and senescence to enhance the productivity of maize and related grasses.

Molecular characterization and expression analysis of the Na+/H+ exchanger gene family in Medicago truncatula.

One important mechanism plants use to cope with salinity is keeping the cytosolic Na+ concentration low by sequestering Na+ in vacuoles, a process facilitated by Na+/H+ exchangers (NHX). There are eight NHX genes (NHX1 through NHX8) identified and characterized in Arabidopsis thaliana. Bioinformatics analyses of the known Arabidopsis genes enabled us to identify six Medicago truncatula NHX genes (MtNHX1, MtNHX2, MtNHX3, MtNHX4, MtNHX6, and MtNHX7). Twelve transmembrane domains and an amiloride binding site were conserved in five out of six MtNHX proteins. Phylogenetic analysis involving A. thaliana, Glycine max, Phaseolus vulgaris, and M. truncatula revealed that each individual MtNHX class (class I: MtNHX1 through 4; class II: MtNHX6; class III: MtNHX7) falls under a separate clade. In a salinity-stress experiment, M. truncatula exhibited ~ 20% reduction in biomass. In the salinity treatment, sodium contents increased by 178 and 75% in leaves and roots, respectively, and Cl- contents increased by 152 and 162%, respectively. Na+ exclusion may be responsible for the relatively smaller increase in Na+ concentration in roots under salt stress as compared to Cl-. Decline in tissue K+ concentration under salinity was not surprising as some antiporters play an important role in transporting both Na+ and K + . MtNHX1, MtNHX6, and MtNHX7 display high expression in roots and leaves. MtNHX3, MtNHX6, and MtNHX7 were induced in roots under salinity stress. Expression analysis results indicate that sequestering Na+ into vacuoles may not be the principal component trait of the salt tolerance mechanism in M. truncatula and other component traits may be pivotal.

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.

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.

Phenotypic and Transcriptional Analysis of Divergently Selected Maize Populations Reveals the Role of Developmental Timing in Seed Size Determination.

Seed size is a component of grain yield and an important trait in crop domestication. To understand the mechanisms governing seed size in maize (Zea mays), we examined transcriptional and developmental changes during seed development in populations divergently selected for large and small seed size from Krug, a yellow dent maize cultivar. After 30 cycles of selection, seeds of the large seed population (KLS30) have a 4.7-fold greater weight and a 2.6-fold larger size compared with the small seed population (KSS30). Patterns of seed weight accumulation from the time of pollination through 30 d of grain filling showed an earlier onset, slower rate, and earlier termination of grain filling in KSS30 relative to KLS30. This was further supported by transcriptome patterns in seeds from the populations and derived inbreds. Although the onset of key genes was earlier in small seeds, similar maximum transcription levels were observed in large seeds at later stages, suggesting that functionally weaker alleles, rather than transcript abundance, may be the basis of the slow rate of seed filling in KSS30. Gene coexpression networks identified several known genes controlling cellularization and proliferation as well as novel genes that will be useful candidates for biotechnological approaches aimed at altering seed size in maize and other cereals.

Maize Unstable factor for orange1 is required for maintaining silencing associated with paramutation at the pericarp color1 and booster1 loci.

To understand the molecular mechanisms underlying paramutation, we examined the role of Unstable factor for orange1 (Ufo1) in maintaining paramutation at the maize pericarp color1 (p1) and booster1 (b1) loci. Genetic tests revealed that the Ufo1-1 mutation disrupted silencing associated with paramutation at both p1 and b1. The level of up regulation achieved at b1 was lower than that at p1, suggesting differences in the role Ufo1-1 plays at these loci. We characterized the interaction of Ufo1-1 with two silenced p1 epialleles, P1-rr' and P1-pr(TP), that were derived from a common P1-rr ancestor. Both alleles are phenotypically indistinguishable, but differ in their paramutagenic activity; P1-rr' is paramutagenic to P1-rr, while P1-pr(TP) is non-paramutagenic. Analysis of cytosine methylation revealed striking differences within an enhancer fragment that is required for paramutation; P1-rr' exhibited increased methylation at symmetric (CG and CHG) and asymmetric (CHH) sites, while P1-pr(TP) was methylated only at symmetric sites. Both silenced alleles had higher levels of dimethylation of lysine 9 on histone 3 (H3K9me2), an epigenetic mark of silent chromatin, in the enhancer region. Both epialleles were reactivated in the Ufo1-1 background; however, reactivation of P1-rr' was associated with dramatic loss of symmetric and asymmetric cytosine methylation in the enhancer, while methylation of up-regulated P1-pr(TP) was not affected. Interestingly, Ufo1-1-mediated reactivation of both alleles was accompanied with loss of H3K9me2 mark from the enhancer region. Therefore, while earlier studies have shown correlation between H3K9me2 and DNA methylation, our study shows that these two epigenetic marks are uncoupled in the Ufo1-1-reactivated p1 alleles. Furthermore, while CHH methylation at the enhancer region appears to be the major distinguishing mark between paramutagenic and non-paramutagenic p1 alleles, H3K9me2 mark appears to be important for maintaining epigenetic silencing.

Biomechanical phenotyping pipeline for stalk lodging resistance in maize.

Stalk lodging (structural failure crops prior to harvest) significantly reduces annual yields of vital grain crops. The lack of standardized, high throughput phenotyping methods capable of quantifying biomechanical plant traits prevents comprehensive understanding of the genetic architecture of stalk lodging resistance. A phenotyping pipeline developed to enable higher throughput biomechanical measurements of plant traits related to stalk lodging is presented. The methods were developed using principles from the fields of engineering mechanics and metrology and they enable retention of plant-specific data instead of averaging data across plots as is typical in most phenotyping studies. This pipeline was specifically designed to be implemented in large experimental studies and has been used to phenotype over 40,000 maize stalks. The pipeline includes both lab- and field-based phenotyping methodologies and enables the collection of metadata. Best practices learned by implementing this pipeline over the past three years are presented. The specific instruments (including model numbers and manufacturers) that work well for these methods are presented, however comparable instruments may be used in conjunction with these methods as seen fit.•Efficient methods to measure biomechanical traits and record metadata related to stalk lodging.•Can be used in studies with large sample sizes (i.e., > 1,000).

Transcriptional and metabolic analysis of senescence induced by preventing pollination in maize.

Transcriptional and metabolic changes were evaluated during senescence induced by preventing pollination in the B73 genotype of maize (Zea mays). Accumulation of free glucose and starch and loss of chlorophyll in leaf was manifested early at 12 d after anthesis (DAA), while global transcriptional and phenotypic changes were evident only at 24 DAA. Internodes exhibited major transcriptomic changes only at 30 DAA. Overlaying expression data onto metabolic pathways revealed involvement of many novel pathways, including those involved in cell wall biosynthesis. To investigate the overlap between induced and natural senescence, transcriptional data from induced senescence in maize was compared with that reported for Arabidopsis (Arabidopsis thaliana) undergoing natural and sugar-induced senescence. Notable similarities with natural senescence in Arabidopsis included up-regulation of senescence-associated genes (SAGs), ethylene and jasmonic acid biosynthetic genes, APETALA2, ethylene-responsive element binding protein, and no apical meristem transcription factors. However, differences from natural senescence were highlighted by unaltered expression of a subset of the SAGs, and cytokinin, abscisic acid, and salicylic acid biosynthesis genes. Key genes up-regulated during sugar-induced senescence in Arabidopsis, including a cysteine protease (SAG12) and three flavonoid biosynthesis genes (PRODUCTION OF ANTHOCYANIN PIGMENT1 (PAP1), PAP2, and LEUCOANTHOCYANIDIN DIOXYGENASE), were also induced, suggesting similarities in senescence induced by pollination prevention and sugar application. Coexpression analysis revealed networks involving known senescence-related genes and novel candidates; 82 of these were shared between leaf and internode networks, highlighting similarities in induced senescence in these tissues. Insights from this study will be valuable in systems biology of senescence in maize and other grasses.

Salinity stress tolerance prediction for biomass-related traits in maize (Zea mays L.) using genome-wide markers.

Maize (Zea mays L.) is the third most important cereal crop after rice (Oryza sativa) and wheat (Triticum aestivum). Salinity stress significantly affects vegetative biomass and grain yield and, therefore, reduces the food and silage productivity of maize. Selecting salt-tolerant genotypes is a cumbersome and time-consuming process that requires meticulous phenotyping. To predict salt tolerance in maize, we estimated breeding values for four biomass-related traits, including shoot length, shoot weight, root length, and root weight under salt-stressed and controlled conditions. A five-fold cross-validation method was used to select the best model among genomic best linear unbiased prediction (GBLUP), ridge-regression BLUP (rrBLUP), extended GBLUP, Bayesian Lasso, Bayesian ridge regression, BayesA, BayesB, and BayesC. Examination of the effect of different marker densities on prediction accuracy revealed that a set of low-density single nucleotide polymorphisms obtained through filtering based on a combination of analysis of variance and linkage disequilibrium provided the best prediction accuracy for all the traits. The average prediction accuracy in cross-validations ranged from 0.46 to 0.77 across the four derived traits. The GBLUP, rrBLUP, and all Bayesian models except BayesB demonstrated comparable levels of prediction accuracy that were superior to the other modeling approaches. These findings provide a roadmap for the deployment and optimization of genomic selection in breeding for salt tolerance in maize.

Genomes to Fields 2022 Maize genotype by Environment Prediction Competition.

OBJECTIVES: The Genomes to Fields (G2F) 2022 Maize Genotype by Environment (GxE) Prediction Competition aimed to develop models for predicting grain yield for the 2022 Maize GxE project field trials, leveraging the datasets previously generated by this project and other publicly available data. DATA DESCRIPTION: This resource used data from the Maize GxE project within the G2F Initiative [1]. The dataset included phenotypic and genotypic data of the hybrids evaluated in 45 locations from 2014 to 2022. Also, soil, weather, environmental covariates data and metadata information for all environments (combination of year and location). Competitors also had access to ReadMe files which described all the files provided. The Maize GxE is a collaborative project and all the data generated becomes publicly available [2]. The dataset used in the 2022 Prediction Competition was curated and lightly filtered for quality and to ensure naming uniformity across years.

2018-2019 field seasons of the Maize Genomes to Fields (G2F) G x E project.

OBJECTIVES: This report provides information about the public release of the 2018-2019 Maize G X E project of the Genomes to Fields (G2F) Initiative datasets. G2F is an umbrella initiative that evaluates maize hybrids and inbred lines across multiple environments and makes available phenotypic, genotypic, environmental, and metadata information. The initiative understands the necessity to characterize and deploy public sources of genetic diversity to face the challenges for more sustainable agriculture in the context of variable environmental conditions. DATA DESCRIPTION: Datasets include phenotypic, climatic, and soil measurements, metadata information, and inbred genotypic information for each combination of location and year. Collaborators in the G2F initiative collected data for each location and year; members of the group responsible for coordination and data processing combined all the collected information and removed obvious erroneous data. The collaborators received the data before the DOI release to verify and declare that the data generated in their own locations was accurate. ReadMe and description files are available for each dataset. Previous years of evaluation are already publicly available, with common hybrids present to connect across all locations and years evaluated since this project's inception.

2020-2021 field seasons of Maize GxE project within the Genomes to Fields Initiative.

OBJECTIVES: This release note describes the Maize GxE project datasets within the Genomes to Fields (G2F) Initiative. The Maize GxE project aims to understand genotype by environment (GxE) interactions and use the information collected to improve resource allocation efficiency and increase genotype predictability and stability, particularly in scenarios of variable environmental patterns. Hybrids and inbreds are evaluated across multiple environments and phenotypic, genotypic, environmental, and metadata information are made publicly available. DATA DESCRIPTION: The datasets include phenotypic data of the hybrids and inbreds evaluated in 30 locations across the US and one location in Germany in 2020 and 2021, soil and climatic measurements and metadata information for all environments (combination of year and location), ReadMe, and description files for each data type. A set of common hybrids is present in each environment to connect with previous evaluations. Each environment had a collaborator responsible for collecting and submitting the data, the GxE coordination team combined all the collected information and removed obvious erroneous data. Collaborators received the combined data to use, verify and declare that the data generated in their own environments was accurate. Combined data is released to the public with minimal filtering to maintain fidelity to the original data.

Genome-wide atlas of transcription during maize development.

Maize is an important model species and a major constituent of human and animal diets. It has also emerged as a potential feedstock and model system for bioenergy research due to recent worldwide interest in developing plant biomass-based, carbon-neutral liquid fuels. To understand how the underlying genome sequence results in specific plant phenotypes, information on the temporal and spatial transcription patterns of genes is crucial. Here we present a comprehensive atlas of global transcription profiles across developmental stages and plant organs. We used a NimbleGen microarray containing 80,301 probe sets to profile transcription patterns in 60 distinct tissues representing 11 major organ systems of inbred line B73. Of the 30,892 probe sets representing the filtered B73 gene models, 91.4% were expressed in at least one tissue. Interestingly, 44.5% of the probe sets were expressed in all tissues, indicating a substantial overlap of gene expression among plant organs. Clustering of maize tissues based on global gene expression profiles resulted in formation of groups of biologically related tissues. We utilized this dataset to examine the expression of genes that encode enzymes in the lignin biosynthetic pathway, and found that expansion of distinct gene families was accompanied by divergent, tissue-specific transcription patterns of the paralogs. This comprehensive expression atlas represents a valuable resource for gene discovery and functional characterization in maize.

Cross-sectional geometry predicts failure location in maize stalks.

BACKGROUND: Stalk lodging (breaking of agricultural plant stalks prior to harvest) is a multi-billion dollar a year problem. Stalk lodging occurs when high winds induce bending moments in the stalk which exceed the bending strength of the plant. Previous biomechanical models of plant stalks have investigated the effect of cross-sectional morphology on stalk lodging resistance (e.g., diameter and rind thickness). However, it is unclear if the location of stalk failure along the length of stem is determined by morphological or compositional factors. It is also unclear if the crops are structurally optimized, i.e., if the plants allocate structural biomass to create uniform and minimal bending stresses in the plant tissues. The purpose of this paper is twofold: (1) to investigate the relationship between bending stress and failure location of maize stalks, and (2) to investigate the potential of phenotyping for internode-level bending stresses to assess lodging resistance. RESULTS: 868 maize specimens representing 16 maize hybrids were successfully tested in bending to failure. Internode morphology was measured, and bending stresses were calculated. It was found that bending stress is highly and positively associated with failure location. A user-friendly computational tool is presented to help plant breeders in phenotyping for internode-level bending stress. Phenotyping for internode-level bending stresses could potentially be used to breed for more biomechanically optimal stalks that are resistant to stalk lodging. CONCLUSIONS: Internode-level bending stress plays a potentially critical role in the structural integrity of plant stems. Equations and tools provided herein enable researchers to account for this phenotype, which has the potential to increase the bending strength of plants without increasing overall structural biomass.

Progressive loss of DNA methylation releases epigenetic gene silencing from a tandemly repeated maize Myb gene.

Maize pericarp color1 (p1) gene, which regulates phlobaphene biosynthesis in kernel pericarp and cob glumes, offers an excellent genetic system to study tissue-specific gene regulation. A multicopy p1 allele, P1-wr (white pericarp/red cob) is epigenetically regulated. Hypomethylation of P1-wr in the presence of Unstable factor for orange1 (Ufo1), leads to ectopic pigmentation of pericarp and other organs. The Ufo1-induced phenotypes show incomplete penetrance and poor expressivity: gain of pigmentation is observed only in a subset of plants carrying Ufo1 mutation, and the extent of pigmentation is highly variable. We show that Ufo1 induces progressive hypomethylation of P1-wr repeats over generations. After five generations of exposure to Ufo1, a 30-40% decrease in CG and CNG methylation was observed in an upstream enhancer and an intron region of P1-wr. Interestingly, such hypomethylation correlated with an increase in penetrance of the Ufo1-induced pigmentation phenotype from approximately 27 to 61%. Expressivity of the Ufo1-induced phenotype also improved markedly as indicated by increased uniformity of pericarp pigmentation in the later generations. Furthermore, the poor expressivity of the Uo1 is associated with mosaic methylation patterns of the P1-wr upstream enhancer in individual cells and distinct P1-wr gene copies. Finally, comparison of methylation among different tissues indicated that Ufo1 induces rapid CG and CNG hypomethylation of P1-wr repeats during plant development. Together, these data indicate that the poor penetrance and expressivity of Ufo1-induced phenotypes is caused by mosaicism of methylation, and progressive mitotic hypomethylation leads to improved meiotic heritability of the mutant phenotype. In duplicated genomes like maize, loss of an epigenetic regulator may produce mosaic patterns due to redundancy of epigenetic regulators and their target sequences. We show here that multiple developmental cycles may be required for complete disruption of suppressed epigenetic states and appearance of heritable phenotypes.