Recruitment of an ancient branching program to suppress carpel development in maize flowers.

Carpels in maize undergo programmed cell death in half of the flowers initiated in ears and in all flowers in tassels. The HD-ZIP I transcription factor gene GRASSY TILLERS1 (GT1) is one of only a few genes known to regulate this process. To identify additional regulators of carpel suppression, we performed a gt1 enhancer screen and found a genetic interaction between gt1 and ramosa3 (ra3). RA3 is a classic inflorescence meristem determinacy gene that encodes a trehalose-6-phosphate (T6P) phosphatase (TPP). Dissection of floral development revealed that ra3 single mutants have partially derepressed carpels, whereas gt1;ra3 double mutants have completely derepressed carpels. Surprisingly, gt1 suppresses ra3 inflorescence branching, revealing a role for gt1 in meristem determinacy. Supporting these genetic interactions, GT1 and RA3 proteins colocalize to carpel nuclei in developing flowers. Global expression profiling revealed common genes misregulated in single and double mutant flowers, as well as in derepressed gt1 axillary meristems. Indeed, we found that ra3 enhances gt1 vegetative branching, similar to the roles for the trehalose pathway and GT1 homologs in the eudicots. This functional conservation over ∼160 million years of evolution reveals ancient roles for GT1-like genes and the trehalose pathway in regulating axillary meristem suppression, later recruited to mediate carpel suppression. Our findings expose hidden pleiotropy of classic maize genes and show how an ancient developmental program was redeployed to sculpt floral form.

Evolutionary Variation in MADS Box Dimerization Affects Floral Development and Protein Abundance in Maize.

Interactions between MADS box transcription factors are critical in the regulation of floral development, and shifting MADS box protein-protein interactions are predicted to have influenced floral evolution. However, precisely how evolutionary variation in protein-protein interactions affects MADS box protein function remains unknown. To assess the impact of changing MADS box protein-protein interactions on transcription factor function, we turned to the grasses, where interactions between B-class MADS box proteins vary. We tested the functional consequences of this evolutionary variability using maize (Zea mays) as an experimental system. We found that differential B-class dimerization was associated with subtle, quantitative differences in stamen shape. In contrast, differential dimerization resulted in large-scale changes to downstream gene expression. Differential dimerization also affected B-class complex composition and abundance, independent of transcript levels. This indicates that differential B-class dimerization affects protein degradation, revealing an important consequence for evolutionary variability in MADS box interactions. Our results highlight complexity in the evolution of developmental gene networks: changing protein-protein interactions could affect not only the composition of transcription factor complexes but also their degradation and persistence in developing flowers. Our results also show how coding change in a pleiotropic master regulator could have small, quantitative effects on development.

Ontogeny and anatomy of Bouteloua (Poaceae: Chloridoideae) species display a basipetal branch formation and a novel modified leaf structure in grasses.

BACKGROUND AND AIMS: Shoot ontogenesis in grasses follows a transition from a vegetative phase into a reproductive phase. Current studies provide insight into how branch and spikelet formation occur during the reproductive phase. However, these studies do not explain all the complex diversity of grass inflorescence forms and are mostly focused on model grasses. Moreover, truncated inflorescences of the non-model grass genus Urochloa (Panicoideae) with formation of primary branches have basipetal initiation of branches. Bouteloua species (Chloridoideae) are non-model grasses that form truncated inflorescences of primary branches with apical vestiges of uncertain homology at the tips of branching events and sterile florets above the lowermost fertile floret. Sterile florets are reduced to rudimentary lemmas composed of three large awns diverging from an awn column. Conflict about the awn column identity of this rudimentary lemma is often addressed in species descriptions of this genus. We test if Bouteloua species can display basipetal initiation of branches and explore the identity of vestiges and the awn column of rudimentary lemmas. METHODS: We surveyed the inflorescence ontogeny and branch/awn anatomy of Bouteloua species and compared results with recent ontogenetic studies of chloridoids. KEY RESULTS: Bouteloua arizonica has florets with basipetal maturation. Branches display basipetal branch initiation and maturation. Branch vestiges are formed laterally by meristems during early branching events. The spikelet meristem forms the awn column of rudimentary lemmas. Vestiges and sterile floret awns have anatomical similarities to C4 leaves. CONCLUSIONS: Basipetal initiation of branches is a novel feature for Chloridoideae grasses. Branch vestiges are novel vegetative grass structures. Sterile floret awn columns are likely to be extensions of the rachilla.

Grass inflorescence architecture and meristem determinacy.

The grass inflorescence is striking not only for its beauty and diversity, but also for its developmental complexity. While models of inflorescence architecture have been proposed in both eudicots and grasses, these are inadequate to fully explain the complex branching events that occur during the development of the grass inflorescence. Key to understanding grass inflorescence architecture is the meristem determinacy/indeterminacy decision, which regulates the number of branching events that occur. Here we review what has been learned about meristem determinacy from grass mutants with defects in inflorescence development. A picture is emerging of a complex network of signaling molecules and meristem identity factors that interact to regulate inflorescence meristem activity, many of which have been modified during crop domestication directly affecting yield traits.

The Maize PI/GLO Ortholog Zmm16/sterile tassel silky ear1 Interacts with the Zygomorphy and Sex Determination Pathways in Flower Development.

In monocots and eudicots, B class function specifies second and third whorl floral organ identity as described in the classic ABCE model. Grass B class APETALA3/DEFICIENS orthologs have been functionally characterized; here, we describe the positional cloning and characterization of a maize (Zea mays) PISTILLATA/GLOBOSA ortholog Zea mays mads16 (Zmm16)/sterile tassel silky ear1 (sts1). We show that, similar to many eudicots, all the maize B class proteins bind DNA as obligate heterodimers and positively regulate their own expression. However, sts1 mutants have novel phenotypes that provide insight into two derived aspects of maize flower development: carpel abortion and floral asymmetry. Specifically, we show that carpel abortion acts downstream of organ identity and requires the growth-promoting factor grassy tillers1 and that the maize B class genes are expressed asymmetrically, likely in response to zygomorphy of grass floral primordia. Further investigation reveals that floral phyllotactic patterning is also zygomorphic, suggesting significant mechanistic differences with the well-characterized models of floral polarity. These unexpected results show that despite extensive study of B class gene functions in diverse flowering plants, novel insights can be gained from careful investigation of homeotic mutants outside the core eudicot model species.

Evolutionary Dynamics of Floral Homeotic Transcription Factor Protein-Protein Interactions.

Protein-protein interactions (PPIs) have widely acknowledged roles in the regulation of development, but few studies have addressed the timing and mechanism of shifting PPIs over evolutionary history. The B-class MADS-box transcription factors, PISTILLATA (PI) and APETALA3 (AP3) are key regulators of floral development. PI-like (PI(L)) and AP3-like (AP3(L)) proteins from a number of plants, including Arabidopsis thaliana (Arabidopsis) and the grass Zea mays (maize), bind DNA as obligate heterodimers. However, a PI(L) protein from the grass relative Joinvillea can bind DNA as a homodimer. To ascertain whether Joinvillea PI(L) homodimerization is an anomaly or indicative of broader trends, we characterized PI(L) dimerization across the Poales and uncovered unexpected evolutionary lability. Both obligate B-class heterodimerization and PI(L) homodimerization have evolved multiple times in the order, by distinct molecular mechanisms. For example, obligate B-class heterodimerization in maize evolved very recently from PI(L) homodimerization. A single amino acid change, fixed during domestication, is sufficient to toggle one maize PI(L) protein between homodimerization and obligate heterodimerization. We detected a signature of positive selection acting on residues preferentially clustered in predicted sites of contact between MADS-box monomers and dimers, and in motifs that mediate MADS PPI specificity in Arabidopsis. Changing one positively selected residue can alter PI(L) dimerization activity. Furthermore, ectopic expression of a Joinvillea PI(L) homodimer in Arabidopsis can homeotically transform sepals into petals. Our results provide a window into the evolutionary remodeling of PPIs, and show that novel interactions have the potential to alter plant form in a context-dependent manner.

Grass inflorescence architecture and evolution: the origin of novel signaling centers.

Contents 367 I. 367 II. 368 III. 370 IV. 371 371 References 371 SUMMARY: A central goal of evo-devo is to understand how morphological diversity arises from existing developmental mechanisms, requiring a clear, predictive explanatory framework of the underlying developmental mechanisms. Despite an ever-increasing literature on genes regulating grass inflorescence development, an effective model of inflorescence patterning is lacking. I argue that the existing framework for grass inflorescence development, which invokes homeotic shifts in multiple distinct meristem identities, obscures a recurring theme emerging from developmental genetic studies in grass models, that is that inflorescence branching is regulated by novel localized signaling centers. Understanding the origin and function of these novel signaling centers will be key to future evo-devo work on the grass inflorescence.

From many, one: genetic control of prolificacy during maize domestication.

A reduction in number and an increase in size of inflorescences is a common aspect of plant domestication. When maize was domesticated from teosinte, the number and arrangement of ears changed dramatically. Teosinte has long lateral branches that bear multiple small ears at their nodes and tassels at their tips. Maize has much shorter lateral branches that are tipped by a single large ear with no additional ears at the branch nodes. To investigate the genetic basis of this difference in prolificacy (the number of ears on a plant), we performed a genome-wide QTL scan. A large effect QTL for prolificacy (prol1.1) was detected on the short arm of chromosome 1 in a location that has previously been shown to influence multiple domestication traits. We fine-mapped prol1.1 to a 2.7 kb "causative region" upstream of the grassy tillers1 (gt1) gene, which encodes a homeodomain leucine zipper transcription factor. Tissue in situ hybridizations reveal that the maize allele of prol1.1 is associated with up-regulation of gt1 expression in the nodal plexus. Given that maize does not initiate secondary ear buds, the expression of gt1 in the nodal plexus in maize may suppress their initiation. Population genetic analyses indicate positive selection on the maize allele of prol1.1, causing a partial sweep that fixed the maize allele throughout most of domesticated maize. This work shows how a subtle cis-regulatory change in tissue specific gene expression altered plant architecture in a way that improved the harvestability of maize.

BARREN STALK FASTIGIATE1 is an AT-hook protein required for the formation of maize ears.

Ears are the seed-bearing inflorescences of maize (Zea mays) plants and represent a crucial component of maize yield. The first step in the formation of ears is the initiation of axillary meristems in the axils of developing leaves. In the classic maize mutant barren stalk fastigiate1 (baf1), first discovered in the 1950s, ears either do not form or, if they do, are partially fused to the main stalk. We positionally cloned Baf1 and found that it encodes a transcriptional regulator containing an AT-hook DNA binding motif. Single coorthologs of Baf1 are found in syntenic regions of brachypodium (Brachypodium distachyon), rice (Oryza sativa), and sorghum (Sorghum bicolor), suggesting that the gene is likely present in all cereal species. Protein-protein interaction assays suggest that BAF1 is capable of forming homodimers and heterodimers with other members of the AT-hook family. Another transcriptional regulator required for ear initiation is the basic helix-loop-helix protein BARREN STALK1 (BA1). Genetic and expression analyses suggest that Baf1 is required to reach a threshold level of Ba1 expression for the initiation of maize ears. We propose that Baf1 functions in the demarcation of a boundary region essential for the specification of a stem cell niche.

grassy tillers1 promotes apical dominance in maize and responds to shade signals in the grasses.

The shape of a plant is largely determined by regulation of lateral branching. Branching architecture can vary widely in response to both genotype and environment, suggesting regulation by a complex interaction of autonomous genetic factors and external signals. Tillers, branches initiated at the base of grass plants, are suppressed in response to shade conditions. This suppression of tiller and lateral branch growth is an important trait selected by early agriculturalists during maize domestication and crop improvement. To understand how plants integrate external environmental cues with endogenous signals to control their architecture, we have begun a functional characterization of the maize mutant grassy tillers1 (gt1). We isolated the gt1 gene using positional cloning and found that it encodes a class I homeodomain leucine zipper gene that promotes lateral bud dormancy and suppresses elongation of lateral ear branches. The gt1 expression is induced by shading and is dependent on the activity of teosinte branched1 (tb1), a major domestication locus controlling tillering and lateral branching. Interestingly, like tb1, gt1 maps to a quantitative trait locus that regulates tillering and lateral branching in maize and shows evidence of selection during maize domestication. Branching and shade avoidance are both of critical agronomic importance, but little is known about how these processes are integrated. Our results indicate that gt1 mediates the reduced branching associated with the shade avoidance response in the grasses. Furthermore, selection at the gt1 locus suggests that it was involved in improving plant architecture during the domestication of maize.

QTL analysis of divergent floral morphology traits between Gilia yorkii and G. capitata.

Speciation is a complex process typically accompanied by significant genetic and morphological differences between sister populations. In plants, divergent floral morphologies and pollinator differences can result in reproductive isolation between populations. Here, we explore floral trait differences between two recently diverged species, Gilia yorkii and G. capitata. The distributions of floral traits in parental, F1, and F2 populations are compared, and groups of correlated traits are identified. We describe the genetic architecture of floral traits through a quantitative trait locus (QTL) analysis using an F2 population of 187 individuals. While all identified QTLs were of moderate (10-25%) effect, interestingly, most QTL intervals were non-overlapping, suggesting that, in general, traits do not share a common genetic basis. Our results provide a framework for future identification of genes involved in the evolution of floral morphology.

Transcription factor binding site divergence across maize inbred lines drives transcriptional and phenotypic variation.

Regulatory elements are important constituents of plant genomes that have shaped ancient and modern crops. Their identification, function, and diversity in crop genomes however are poorly characterized, thus limiting our ability to harness their power for further agricultural advances using induced or natural variation. Here, we use DNA affinity purification-sequencing (DAP-seq) to map transcription factor (TF) binding events for 200 maize TFs belonging to 30 distinct families and heterodimer pairs in two distinct inbred lines historically used for maize hybrid plant production, providing empirical binding site annotation for 5.3% of the maize genome. TF binding site comparison in B73 and Mo17 inbreds reveals widespread differences, driven largely by structural variation, that correlate with gene expression changes. TF binding site presence-absence variation helps clarify complex QTL such as vgt1, an important determinant of maize flowering time, and DICE, a distal enhancer involved in herbivore resistance. Modification of TF binding regions via CRISPR-Cas9 mediated editing alters target gene expression and phenotype. Our functional catalog of maize TF binding events enables collective and comparative TF binding analysis, and highlights its value for agricultural improvement.

The maize SBP-box transcription factor encoded by tasselsheath4 regulates bract development and the establishment of meristem boundaries.

Plant architecture consists of repeating units called phytomers, each containing an internode, leaf and axillary meristem. The formation of boundaries within the phytomer is necessary to differentiate and separate these three components, otherwise some will grow at the expense of others. The microRNA-targeted SBP-box transcription factor tasselsheath4 (tsh4) plays an essential role in establishing these boundaries within the inflorescence. tsh4 mutants display altered phyllotaxy, fewer lateral meristems and ectopic leaves that grow at the expense of the meristem. Double-mutant analyses of tsh4 and several highly branched mutants, such as ramosa1-3 and branched silkless1, demonstrated a requirement for tsh4 in branch meristem initiation and maintenance. TSH4 protein, however, was localized throughout the inflorescence stem and at the base of lateral meristems, but not within the meristem itself. Double labeling of TSH4 with the ramosa2, branched silkless1 and knotted1 meristem markers confirmed that TSH4 forms a boundary adjacent to all lateral meristems. Indeed, double labeling of miR156 showed a meristem-specific pattern complementary to that of TSH4, consistent with tsh4 being negatively regulated by this microRNA. Thus, downregulation of TSH4 by a combination of microRNAs and branching pathway genes allows the establishment of lateral meristems and the repression of leaf initiation, thereby playing a major role in defining meristem versus leaf boundaries.

A conserved mechanism of bract suppression in the grass family.

Suppression of inflorescence leaf, or bract, growth has evolved multiple times in diverse angiosperm lineages, including the Poaceae and Brassicaceae. Studies of Arabidopsis thaliana mutants have revealed several genes involved in bract suppression, but it is not known if these genes play a similar role in other plants with suppressed bracts. We identified maize (Zea mays) tassel sheath (tsh) mutants, characterized by the loss of bract suppression, that comprise five loci (tsh1-tsh5). We used map-based cloning to identify Tsh1 and found that it encodes a GATA zinc-finger protein, a close homolog of HANABA TARANU (HAN) of Arabidopsis. The bract suppression function of Tsh1 is conserved throughout the grass family, as we demonstrate that the rice (Oryza sativa) NECK LEAF1 (NL1) and barley (Hordeum vulgare) THIRD OUTER GLUME (TRD) genes are orthologous with Tsh1. Interestingly, NL1/Tsh1/TRD expression and function are not conserved with HAN. The existence of paralogous NL1/Tsh1/TRD-like genes in the grasses indicates that the NL1/Tsh1/TRD lineage was created by recent duplications that may have facilitated its neofunctionalization. A comparison with the Arabidopsis genes regulating bract suppression further supports the hypothesis that the convergent evolution of bract suppression in the Poaceae involved recruitment of a distinct genetic pathway.

Chromosome-Scale Genome Assembly of Gilia yorkii Enables Genetic Mapping of Floral Traits in an Interspecies Cross.

Substantial morphological variation in land plants remains inaccessible to genetic analysis because current models lack variation in important ecological and agronomic traits. The genus Gilia was historically a model for biosystematics studies and includes variation in morphological traits that are poorly understood at the genetic level. We assembled a chromosome-scale reference genome of G. yorkii and used it to investigate genome evolution in the Polemoniaceae. We performed QTL (quantitative trait loci) mapping in a G. yorkii×G. capitata interspecific population for traits related to inflorescence architecture and flower color. The genome assembly spans 2.75 Gb of the estimated 2.80-Gb genome, with 96.7% of the sequence contained in the nine largest chromosome-scale scaffolds matching the haploid chromosome number. Gilia yorkii experienced at least one round of whole-genome duplication shared with other Polemoniaceae after the eudicot paleohexaploidization event. We identified QTL linked to variation in inflorescence architecture and petal color, including a candidate for the major flower color QTL-a tandem duplication of flavanol 3',5'-hydroxylase. Our results demonstrate the utility of Gilia as a forward genetic model for dissecting the evolution of development in plants including the causal loci underlying inflorescence architecture transitions.

Boundary domain genes were recruited to suppress bract growth and promote branching in maize.

Grass inflorescence development is diverse and complex and involves sophisticated but poorly understood interactions of genes regulating branch determinacy and leaf growth. Here, we use a combination of transcript profiling and genetic and phylogenetic analyses to investigate tasselsheath1 (tsh1) and tsh4, two maize genes that simultaneously suppress inflorescence leaf growth and promote branching. We identify a regulatory network of inflorescence leaf suppression that involves the phase change gene tsh4 upstream of tsh1 and the ligule identity gene liguleless2 (lg2). We also find that a series of duplications in the tsh1 gene lineage facilitated its shift from boundary domain in nongrasses to suppressed inflorescence leaves of grasses. Collectively, these results suggest that the boundary domain genes tsh1 and lg2 were recruited to inflorescence leaves where they suppress growth and regulate a nonautonomous signaling center that promotes inflorescence branching, an important component of yield in cereal grasses.

Integration of high-density genetic mapping with transcriptome analysis uncovers numerous agronomic QTL and reveals candidate genes for the control of tillering in sorghum.

Phenotypes such as branching, photoperiod sensitivity, and height were modified during plant domestication and crop improvement. Here, we perform quantitative trait locus (QTL) mapping of these and other agronomic traits in a recombinant inbred line (RIL) population derived from an interspecific cross between Sorghum propinquum and Sorghum bicolor inbred Tx7000. Using low-coverage Illumina sequencing and a bin-mapping approach, we generated ∼1920 bin markers spanning ∼875 cM. Phenotyping data were collected and analyzed from two field locations and one greenhouse experiment for six agronomic traits, thereby identifying a total of 30 QTL. Many of these QTL were penetrant across environments and co-mapped with major QTL identified in other studies. Other QTL uncovered new genomic regions associated with these traits, and some of these were environment-specific in their action. To further dissect the genetic underpinnings of tillering, we complemented QTL analysis with transcriptomics, identifying 6189 genes that were differentially expressed during tiller bud elongation. We identified genes such as Dormancy Associated Protein 1 (DRM1) in addition to various transcription factors that are differentially expressed in comparisons of dormant to elongating tiller buds and lie within tillering QTL, suggesting that these genes are key regulators of tiller elongation in sorghum. Our study demonstrates the usefulness of this RIL population in detecting domestication and improvement-associated genes in sorghum, thus providing a valuable resource for genetic investigation and improvement to the sorghum community.

Development of an Evolutionary Tree Concept Inventory.

Despite the importance of tree-thinking and evolutionary trees to biology, no appropriately developed concept inventory exists to measure student understanding of these important concepts. To address this need, we developed a multiple-choice concept inventory consisting of 24 pairs of items, and we provide evidence to support its use among undergraduate students. A set of learning outcomes was developed to guide the creation of the concept inventory. The learning outcomes, student interviews, and student responses were used to develop and revise inventory items. Supporting evidence was gathered from traditional item analysis, exploratory factor analysis, confirmatory factor analysis, traditional reliability analyses, and comparisons to alternative assessments. Appropriate implementation and utility of the concept inventory are discussed.

The regulatory landscape of a core maize domestication module controlling bud dormancy and growth repression.

Many domesticated crop plants have been bred for increased apical dominance, displaying greatly reduced axillary branching compared to their wild ancestors. In maize, this was achieved through selection for a gain-of-function allele of the TCP transcription factor teosinte branched1 (tb1). The mechanism for how a dominant Tb1 allele increased apical dominance, is unknown. Through ChIP seq, RNA seq, hormone and sugar measurements on 1 mm axillary bud tissue, we identify the genetic pathways putatively regulated by TB1. These include pathways regulating phytohormones such as gibberellins, abscisic acid and jasmonic acid, but surprisingly, not auxin. In addition, metabolites involved in sugar sensing such as trehalose 6-phosphate were increased. This suggests that TB1 induces bud suppression through the production of inhibitory phytohormones and by reducing sugar levels and energy balance. Interestingly, TB1 also putatively targets several other domestication loci, including teosinte glume architecture1, prol1.1/grassy tillers1, as well as itself. This places tb1 on top of the domestication hierarchy, demonstrating its critical importance during the domestication of maize from teosinte.

Bulked-Segregant Analysis Coupled to Whole Genome Sequencing (BSA-Seq) for Rapid Gene Cloning in Maize.

Forward genetics remains a powerful method for revealing the genes underpinning organismal form and function, and for revealing how these genes are tied together in gene networks. In maize, forward genetics has been tremendously successful, but the size and complexity of the maize genome made identifying mutant genes an often arduous process with traditional methods. The next generation sequencing revolution has allowed for the gene cloning process to be significantly accelerated in many organisms, even when genomes are large and complex. Here, we describe a bulked-segregant analysis sequencing (BSA-Seq) protocol for cloning mutant genes in maize. Our simple strategy can be used to quickly identify a mapping interval and candidate single nucleotide polymorphisms (SNPs) from whole genome sequencing of pooled F2 individuals. We employed this strategy to identify narrow odd dwarf as an enhancer of teosinte branched1, and to identify a new allele of defective kernel1 Our method provides a quick, simple way to clone genes in maize.