Search

Abstract Zoologists have adduced morphological convergence among embryonic stages of closely related taxa, which has been called the phylotypic stage of embryogenesis. Transcriptomic analyzes reveal an hourglass pattern of gene expression during plant and animal embryogenesis, characterized by the accumulation of evolutionarily older and conserved transcripts during mid-embryogenesis, whereas younger less-conserved transcripts predominate at earlier and later embryonic stages. In contrast, comparisons of embryonic gene expression among different animal phyla describe an inverse hourglass pattern, where expression is correlated during early and late stages but not during mid-embryo development. Here, multiplexed spatial-transcriptomic analyzes is used to investigate embryogenesis and homology in maize, which has grass-specific morphology. A set of shared, co-expressed genes is identified during initiation of maize embryonic organs, replete for ancient/conserved genes manifesting an hourglass pattern during mid-embryogenesis. Transcriptomic comparisons of maize and Arabidopsis embryogenesis with that of the moss Physcomitrium patens identify an inverse hourglass pattern across plant phyla, as in animals. The data suggest that the phylotypic stages in plants and animals are characterized by expression of ancient and conserved genes during histogenesis, organization of embryonic axes, and initial morphogenesis. We propose a mechanism for gene evolution during the innovation of morphological novelty.

Social Impact Statement Although tremendously successful at feeding humanity, row crop agriculture based on annuals contributes to numerous ecosystem dis‐services, ranging from soil degradation and aquatic eutrophication to greenhouse gas production. In contrast, perennial grain crops (which produce harvests for multiple seasons from single plantings) have the potential to provide valuable regulating and supporting ecosystem services in addition to food production. In particular, losses of ecological capital that threaten permanent food insecurity such as ~1% of global soil per year are expected to be mitigated or even reversed by crops that combine the high yield realized by scientific breeding via multiple cropping cycles from single plantings. Summary Perennial herbaceous may provide food and biomass while preserving ecological capital and reducing energy inputs. Sorghum has two perennial relatives and rich morphological diversity being used to breed for perenniality. We elucidate genetic determinants of rhizomatousness and survival, in two BC1F2 populations totaling 246 genotypes derived from backcrossing different annual Sorghum bicolor X perennial S. halepense F1 plants to a tetraploidized S. bicolor. RNA‐seq assisted in identifying candidate genes for rhizomatousness. Correspondence of rhizomatousness quantitative trait loci (QTLs) with those from two populations derived from crosses between S. halepense progenitors S. bicolor X S. propinquum suggests either the preservation of interspecific polymorphism or the formation of novel alleles following polyploid S. halepense formation. Correspondence of tillering and branching QTLs further supports their developmental. Identification of genes from RNA‐seq study within QTL intervals provides insight toward discovery of causal rhizomatous genes. An unexpected finding from both S. halepense‐ and S. propinquum‐derived populations is that alleles contributing to late flowering are related to reduced rhizomatousness. Twelve of 16 QTL regions conferring rhizomatousness fall in paleo‐duplicated regions tracing to single ancestral regions 96 million years ago, indicating that corresponding genes in these regions have retained similar functions since the duplication event.

The cuticle, a hydrophobic layer of cutin and waxes synthesized by plant epidermal cells, is the major barrier to water loss when stomata are closed at night and under water-limited conditions. Elucidating the genetic architecture of natural variation for leaf cuticular conductance (gc) is important for identifying genes relevant to improving crop productivity in drought-prone environments. To this end, we conducted a genome-wide association study of gc of adult leaves in a maize inbred association panel that was evaluated in four environments (Maricopa, AZ, and San Diego, CA, in 2016 and 2017). Five genomic regions significantly associated with gc were resolved to seven plausible candidate genes (ISTL1, two SEC14 homologs, cyclase-associated protein, a CER7 homolog, GDSL lipase, and β-D-XYLOSIDASE 4). These candidates are potentially involved in cuticle biosynthesis, trafficking and deposition of cuticle lipids, cutin polymerization, and cell wall modification. Laser microdissection RNA sequencing revealed that all these candidate genes, with the exception of the CER7 homolog, were expressed in the zone of the expanding adult maize leaf where cuticle maturation occurs. With direct application to genetic improvement, moderately high average predictive abilities were observed for whole-genome prediction of gc in locations (0.46 and 0.45) and across all environments (0.52). The findings of this study provide novel insights into the genetic control of gc and have the potential to help breeders more effectively develop drought-tolerant maize for target environments.

Bulliform cells comprise specialized cell types that develop on the adaxial (upper) surface of grass leaves, and are patterned to form linear rows along the proximodistal axis of the adult leaf blade. Bulliform cell patterning affects leaf angle and is presumed to function during leaf rolling, thereby reducing water loss during temperature extremes and drought. In this study, epidermal leaf impressions were collected from a genetically and anatomically diverse population of maize inbred lines. Subsequently, convolutional neural networks were employed to measure microscopic, bulliform cell-patterning phenotypes in high-throughput. A genome-wide association study, combined with RNAseq analyses of the bulliform cell ontogenic zone, identified candidate regulatory genes affecting bulliform cell column number and cell width. This study is the first to combine machine learning approaches, transcriptomics, and genomics to study bulliform cell patterning, and the first to utilize natural variation to investigate the genetic architecture of this microscopic trait. In addition, this study provides insight toward the improvement of macroscopic traits such as drought resistance and plant architecture in an agronomically important crop plant.

Abstract Background Plant domestication provides a unique model to study genome evolution. Many studies have been conducted to examine genes, genetic diversity, genome structure, and epigenome changes associated with domestication. Interestingly, domesticated accessions have significantly higher [A] and [T] values across genome-wide polymorphic sites than accessions sampled from the corresponding progenitor species. However, the relative contributions of different genomic regions to this genome divergence pattern and underlying mechanisms have not been well characterized. Results Here, we investigate the genome-wide base-composition patterns by analyzing millions of SNPs segregating among 100 accessions from a teosinte-maize comparison set and among 302 accessions from a wild-domesticated soybean comparison set. We show that non-genic part of the genome has a greater contribution than genic SNPs to the [AT]-increase observed between wild and domesticated accessions in maize and soybean. The separation between wild and domesticated accessions in [AT] values is significantly enlarged in non-genic and pericentromeric regions. Motif frequency and sequence context analyses show the motifs (PyCG) related to solar-UV signature are enriched in these regions, particularly when they are methylated. Additional analysis using population-private SNPs also implicates the role of these motifs in relatively recent mutations. With base-composition across polymorphic sites as a genome phenotype, genome scans identify a set of putative candidate genes involved in UV damage repair pathways. Conclusions The [AT]-increase is more pronounced in genomic regions that are non-genic, pericentromeric, transposable elements; methylated; and with low recombination. Our findings establish important links among UV radiation, mutation, DNA repair, methylation, and genome evolution.

From noble beginnings as a prospective forage, polyploid Sorghum halepense (‘Johnsongrass’) is both an invasive species and one of the world’s worst agricultural weeds. Formed by S. bicolor x S. propinquum hybridization, we show S. halepense to have S. bicolor-enriched allele composition and striking mutations in 5,957 genes that differentiate it from representatives of its progenitor species and an outgroup. The spread of S. halepense may have been facilitated by introgression from closely-related cultivated sorghum near genetic loci affecting rhizome development, seed size, and levels of lutein, a photochemical protectant and abscisic acid precursor. Rhizomes, subterranean stems that store carbohydrates and spawn clonal propagules, have growth correlated with reproductive rather than other vegetative tissues, and increase survival of both temperate cold seasons and tropical dry seasons. Rhizomes of S. halepense are more extensive than those of its rhizomatous progenitor S. propinquum, with gene expression including many alleles from its non-rhizomatous S. bicolor progenitor. The first surviving polyploid in its lineage in ∼96 million years, its post-Columbian spread across six continents carried rich genetic diversity that in the United States has facilitated transition from agricultural to non-agricultural niches. Projected to spread another 200–600 km northward in the coming century, despite its drawbacks S. halepense may offer novel alleles and traits of value to improvement of sorghum.

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.

Abstract Background There are significant limitations in existing methods for the genome-wide identification of genes whose expression patterns affect traits. Results The transcriptomes of five tissues from 27 genetically diverse maize inbred lines were deeply sequenced to identify genes exhibiting high and low levels of expression variation across tissues or genotypes. Transcription factors are enriched among genes with the most variation in expression across tissues, as well as among genes with higher-than-median levels of variation in expression across genotypes. In contrast, transcription factors are depleted among genes whose expression is either highly stable or highly variable across genotypes. We developed a Bayesian-based method for genome-wide association studies (GWAS) in which RNA-seq-based measures of transcript accumulation are used as explanatory variables (eRD-GWAS). The ability of eRD-GWAS to identify true associations between gene expression variation and phenotypic diversity is supported by analyses of RNA co-expression networks, protein–protein interaction networks, and gene regulatory networks. Genes associated with 13 traits were identified using eRD-GWAS on a panel of 369 maize inbred lines. Predicted functions of many of the resulting trait-associated genes are consistent with the analyzed traits. Importantly, transcription factors are significantly enriched among trait-associated genes identified with eRD-GWAS. Conclusions eRD-GWAS is a powerful tool for associating genes with traits and is complementary to SNP-based GWAS. Our eRD-GWAS results are consistent with the hypothesis that genetic variation in transcription factor expression contributes substantially to phenotypic diversity.