Major impacts of widespread structural variation on sorghum

Genetic diversity is critical to crop breeding and improvement, and dissection of the genomic variation underlying agronomic traits can both assist breeding and give insight into basic biological mechanisms. Although recent genome analyses in plants reveal many structural variants (SVs), most current studies of crop genetic variation are dominated by single-nucleotide polymorphisms (SNPs). The extent of the impact of SVs on global trait variation, as well as their utility in genome-wide selection, is not yet understood. In this study, we built an SV data set based on whole-genome resequencing of diverse sorghum lines (n = 363), validated the correlation of photoperiod sensitivity and variety type, and identified SV hotspots underlying the divergent evolution of cellulosic and sweet sorghum. In addition, we showed the complementary contribution of SVs for heritability of traits related to sorghum adaptation. Importantly, inclusion of SV polymorphisms in association studies revealed genotype-phenotype associations not observed with SNPs alone. Three-way genome-wide association studies (GWAS) based on whole-genome SNP, SV, and integrated SNP + SV data sets showed substantial associations between SVs and sorghum traits. The addition of SVs to GWAS substantially increased heritability estimates for some traits, indicating their important contribution to functional allelic variation at the genome level. Our discovery of the widespread impacts of SVs on heritable gene expression variation could render a plausible mechanism for their disproportionate impact on phenotypic variation. This study expands our knowledge of SVs and emphasizes the extensive impacts of SVs on sorghum.

Topological data analysis expands the genotype to phenotype map for 3D maize root system architecture

A central goal of biology is to understand how genetic variation produces phenotypic variation, which has been described as a genotype to phenotype (G to P) map. The plant form is continuously shaped by intrinsic developmental and extrinsic environmental inputs, and therefore plant phenomes are highly multivariate and require comprehensive approaches to fully quantify. Yet a common assumption in plant phenotyping efforts is that a few pre-selected measurements can adequately describe the relevant phenome space. Our poor understanding of the genetic basis of root system architecture is at least partially a result of this incongruence. Root systems are complex 3D structures that are most often studied as 2D representations measured with relatively simple univariate traits. In prior work, we showed that persistent homology, a topological data analysis method that does not pre-suppose the salient features of the data, could expand the phenotypic trait space and identify new G to P relations from a commonly used 2D root phenotyping platform. Here we extend the work to entire 3D root system architectures of maize seedlings from a mapping population that was designed to understand the genetic basis of maize-nitrogen relations. Using a panel of 84 univariate traits, persistent homology methods developed for 3D branching, and multivariate vectors of the collective trait space, we found that each method captures distinct information about root system variation as evidenced by the majority of non-overlapping QTL, and hence that root phenotypic trait space is not easily exhausted. The work offers a data-driven method for assessing 3D root structure and highlights the importance of non-canonical phenotypes for more accurate representations of the G to P map.

Transformation and gene editing in the bioenergy grass Miscanthus

BACKGROUND

Miscanthus, a C4 member of Poaceae, is a promising perennial crop for bioenergy, renewable bioproducts, and carbon sequestration. Species of interest include nothospecies M. x giganteus and its parental species M. sacchariflorus and M. sinensis. Use of biotechnology-based procedures to genetically improve Miscanthus, to date, have only included plant transformation procedures for introduction of exogenous genes into the host genome at random, non-targeted sites.

RESULTS

We developed gene editing procedures for Miscanthus using CRISPR/Cas9 that enabled the mutation of a specific (targeted) endogenous gene to knock out its function. Classified as paleo-allopolyploids (duplicated ancient Sorghum-like DNA plus chromosome fusion event), design of guide RNAs (gRNAs) for Miscanthus needed to target both homeologs and their alleles to account for functional redundancy. Prior research in Zea mays demonstrated that editing the lemon white1 (lw1) gene, involved in chlorophyll and carotenoid biosynthesis, via CRISPR/Cas9 yielded pale green/yellow, striped or white leaf phenotypes making lw1 a promising target for visual confirmation of editing in other species. Using sequence information from both Miscanthus and sorghum, orthologs of maize lw1 were identified; a multi-step screening approach was used to select three gRNAs that could target homeologs of lw1. Embryogenic calli of M. sacchariflorus, M. sinensis and M. x giganteus were transformed via particle bombardment (biolistics) or Agrobacterium tumefaciens introducing the Cas9 gene and three gRNAs to edit lw1. Leaves on edited Miscanthus plants displayed the same phenotypes noted in maize. Sanger sequencing confirmed editing; deletions in lw1 ranged from 1 to 26 bp in length, and one deletion (433 bp) encompassed two target sites. Confocal microscopy verified lack of autofluorescence (chlorophyll) in edited leaves/sectors.

CONCLUSIONS

We developed procedures for gene editing via CRISPR/Cas9 in Miscanthus and, to the best of our knowledge, are the first to do so. This included five genotypes representing three Miscanthus species. Designed gRNAs targeted all copies of lw1 (homeologous copies and their alleles); results also confirmed lw1 made a good editing target in species other than Z. mays. The ability to target specific loci to enable endogenous gene editing presents a new avenue for genetic improvement of this important biomass crop.

Small RNAs mediate transgenerational inheritance of genome-wide trans-acting epialleles in maize

Background

Hybridization and backcrossing are commonly used in animal and plant breeding to induce heritable variation including epigenetic changes such as paramutation. However, the molecular basis for hybrid-induced epigenetic memory remains elusive.

Results

Here, we report that hybridization between the inbred parents B73 and Mo17 induces trans-acting hypermethylation and hypomethylation at thousands of loci; several hundreds (~ 3%) are transmitted through six backcrossing and three selfing generations. Notably, many transgenerational methylation patterns resemble epialleles of the nonrecurrent parent, despite > 99% of overall genomic loci are converted to the recurrent parent. These epialleles depend on 24-nt siRNAs, which are eliminated in the isogenic hybrid Mo17xB73:mop1-1 that is defective in siRNA biogenesis. This phenomenon resembles paramutation-like events and occurs in both intraspecific (Mo17xB73) and interspecific (W22xTeosinte) hybrid maize populations. Moreover, siRNA abundance and methylation levels of these epialleles can affect expression of their associated epigenes, many of which are related to stress responses.

Conclusion

Divergent siRNAs between the hybridizing parents can induce trans-acting epialleles in the hybrids, while the induced epigenetic status is maintained for transgenerational inheritance during backcross and hybrid breeding, which alters epigene expression to enhance growth and adaptation. These genetic and epigenetic principles may apply broadly from plants to animals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02614-0.

Evolutionarily informed machine learning enhances the power of predictive gene-to-phenotype relationships

Inferring phenotypic outcomes from genomic features is both a promise and challenge for systems biology. Using gene expression data to predict phenotypic outcomes, and functionally validating the genes with predictive powers are two challenges we address in this study. We applied an evolutionarily informed machine learning approach to predict phenotypes based on transcriptome responses shared both within and across species. Specifically, we exploited the phenotypic diversity in nitrogen use efficiency and evolutionarily conserved transcriptome responses to nitrogen treatments across Arabidopsis accessions and maize varieties. We demonstrate that using evolutionarily conserved nitrogen responsive genes is a biologically principled approach to reduce the feature dimensionality in machine learning that ultimately improved the predictive power of our gene-to-trait models. Further, we functionally validated seven candidate transcription factors with predictive power for NUE outcomes in Arabidopsis and one in maize. Moreover, application of our evolutionarily informed pipeline to other species including rice and mice models underscores its potential to uncover genes affecting any physiological or clinical traits of interest across biology, agriculture, or medicine.

The importance of dominance and genotype-by-environment interactions on grain yield variation in a large-scale public cooperative maize experiment

High-dimensional and high-throughput genomic, field performance, and environmental data are becoming increasingly available to crop breeding programs, and their integration can facilitate genomic prediction within and across environments and provide insights into the genetic architecture of complex traits and the nature of genotype-by-environment interactions. To partition trait variation into additive and dominance (main effect) genetic and corresponding genetic-by-environment variances, and to identify specific environmental factors that influence genotype-by-environment interactions, we curated and analyzed genotypic and phenotypic data on 1918 maize (Zea mays L.) hybrids and environmental data from 65 testing environments. For grain yield, dominance variance was similar in magnitude to additive variance, and genetic-by-environment variances were more important than genetic main effect variances. Models involving both additive and dominance relationships best fit the data and modeling unique genetic covariances among all environments provided the best characterization of the genotype-by-environment interaction patterns. Similarity of relative hybrid performance among environments was modeled as a function of underlying weather variables, permitting identification of weather covariates driving correlations of genetic effects across environments. The resulting models can be used for genomic prediction of mean hybrid performance across populations of environments tested or for environment-specific predictions. These results can also guide efforts to incorporate high-throughput environmental data into genomic prediction models and predict values in new environments characterized with the same environmental characteristics.

Utility of Climatic Information via Combining Ability Models to Improve Genomic Prediction for Yield Within the Genomes to Fields Maize Project

Genomic prediction provides an efficient alternative to conventional phenotypic selection for developing improved cultivars with desirable characteristics. New and improved methods to genomic prediction are continually being developed that attempt to deal with the integration of data types beyond genomic information. Modern automated weather systems offer the opportunity to capture continuous data on a range of environmental parameters at specific field locations. In principle, this information could characterize training and target environments and enhance predictive ability by incorporating weather characteristics as part of the genotype-by-environment (G×E) interaction component in prediction models. We assessed the usefulness of including weather data variables in genomic prediction models using a naïve environmental kinship model across 30 environments comprising the Genomes to Fields (G2F) initiative in 2014 and 2015. Specifically four different prediction scenarios were evaluated (i) tested genotypes in observed environments; (ii) untested genotypes in observed environments; (iii) tested genotypes in unobserved environments; and (iv) untested genotypes in unobserved environments. A set of 1,481 unique hybrids were evaluated for grain yield. Evaluations were conducted using five different models including main effect of environments; general combining ability (GCA) effects of the maternal and paternal parents modeled using the genomic relationship matrix; specific combining ability (SCA) effects between maternal and paternal parents; interactions between genetic (GCA and SCA) effects and environmental effects; and finally interactions between the genetics effects and environmental covariates. Incorporation of the genotype-by-environment interaction term improved predictive ability across all scenarios. However, predictive ability was not improved through inclusion of naive environmental covariates in G×E models. More research should be conducted to link the observed weather conditions with important physiological aspects in plant development to improve predictive ability through the inclusion of weather data.

Genome biology of the paleotetraploid perennial biomass crop Miscanthus

Miscanthus is a perennial wild grass that is of global importance for paper production, roofing, horticultural plantings, and an emerging highly productive temperate biomass crop. We report a chromosome-scale assembly of the paleotetraploid M. sinensis genome, providing a resource for Miscanthus that links its chromosomes to the related diploid Sorghum and complex polyploid sugarcanes. The asymmetric distribution of transposons across the two homoeologous subgenomes proves Miscanthus paleo-allotetraploidy and identifies several balanced reciprocal homoeologous exchanges. Analysis of M. sinensis and M. sacchariflorus populations demonstrates extensive interspecific admixture and hybridization, and documents the origin of the highly productive triploid bioenergy crop M. × giganteus. Transcriptional profiling of leaves, stem, and rhizomes over growing seasons provides insight into rhizome development and nutrient recycling, processes critical for sustainable biomass accumulation in a perennial temperate grass. The Miscanthus genome expands the power of comparative genomics to understand traits of importance to Andropogoneae grasses.

The perennial grass Miscanthus is a promising biomass crop. Here, via genomics and transcriptomics, the authors reveal its allotetraploid origin, characterize gene expression associated with rhizome development and nutrient recycling, and describe the hybrid origin of the triploid M. x giganteus.

Maize genomes to fields (G2F): 2014-2017 field seasons: genotype, phenotype, climatic, soil, and inbred ear image datasets

OBJECTIVES

Advanced tools and resources are needed to efficiently and sustainably produce food for an increasing world population in the context of variable environmental conditions. The maize genomes to fields (G2F) initiative is a multi-institutional initiative effort that seeks to approach this challenge by developing a flexible and distributed infrastructure addressing emerging problems. G2F has generated large-scale phenotypic, genotypic, and environmental datasets using publicly available inbred lines and hybrids evaluated through a network of collaborators that are part of the G2F's genotype-by-environment (G × E) project. This report covers the public release of datasets for 2014-2017.

DATA DESCRIPTION

Datasets include inbred genotypic information; phenotypic, climatic, and soil measurements and metadata information for each testing location across years. For a subset of inbreds in 2014 and 2015, yield component phenotypes were quantified by image analysis. Data released are accompanied by README descriptions. For genotypic and phenotypic data, both raw data and a version without outliers are reported. For climatic data, a version calibrated to the nearest airport weather station and a version without outliers are reported. The 2014 and 2015 datasets are updated versions from the previously released files [1] while 2016 and 2017 datasets are newly available to the public.

Identification and characterization of a novel stay-green QTL that increases yield in maize

Functional stay-green is a valuable trait that extends the photosynthetic period, increases source capacity and biomass and ultimately translates to higher grain yield. Selection for higher yields has increased stay-green in modern maize hybrids. Here, we report a novel QTL controlling functional stay-green that was discovered in a mapping population derived from the Illinois High Protein 1 (IHP1) and Illinois Low Protein 1 (ILP1) lines, which show very different rates of leaf senescence. This QTL was mapped to a single gene containing a NAC-domain transcription factor that we named nac7. Transgenic maize lines where nac7 was down-regulated by RNAi showed delayed senescence and increased both biomass and nitrogen accumulation in vegetative tissues, demonstrating NAC7 functions as a negative regulator of the stay-green trait. More importantly, crosses between nac7 RNAi parents and two different elite inbred testers produced hybrids with prolonged stay-green and increased grain yield by an average 0.29 megagram/hectare (4.6 bushel/acre), in 2 years of multi-environment field trials. Subsequent RNAseq experiments, one employing nac7 RNAi leaves and the other using leaf protoplasts overexpressing Nac7, revealed an important role for NAC7 in regulating genes in photosynthesis, chlorophyll degradation and protein turnover pathways that each contribute to the functional stay-green phenotype. We further determined the putative target of NAC7 and provided a logical extension for the role of NAC7 in regulating resource allocation from vegetative source to reproductive sink tissues. Collectively, our findings make a compelling case for NAC7 as a target for improving functional stay-green and yields in maize and other crops.

Nutritional evaluation of 3 types of novel ethanol coproducts

Five 48-h precision-fed rooster assays were conducted with the objective to determine true metabolizable energy (TMEn) using conventional roosters and/or standardized amino acid digestibility using cecectomized roosters for distillers dried grains with solubles (DDGS) produced from human food waste at high solids content (FWDDGS), DDGS produced from 4 corn hybrids with increases in grain protein concentration or the concentrations of several dietary indispensable amino acids, and a coproduct that is produced by a process which separates a high protein and yeast fraction from ethanol stillage (Still Pro, Fluid Quip Process Technologies, Cedar Rapids, IA). These results from the first cecectomized rooster assay indicated that the standardized digestibility values for Lys, Met, Cys, Thr, and Val were 61, 75, 70, 70, and 72%, respectively. Using conventional roosters, the TMEn for the FWDDGS was 3,890 kcal/kg DM. The DDGS produced from the high protein mutant corn hybrid had a higher protein content of 34% compared with 28% protein for DDGS from the control corn hybrid. Using Lys as an example, there was a large difference between the 2 samples; the high protein mutant DDGS contained 1.60% Lys vs 1.05% for the control DDGS. Standardized digestibility of amino acids was generally not different for the mutant DDGS and the control DDGS. Similar results were observed for Lys, Arg, and Trp for the DDGS produced from the second mutant corn hybrid. The Still Pro sample was analyzed to contain 53% protein (DM basis) with 2.22% Lys, 1.05% Met, 0.90% Cys, 2.06% Thr, and 3.08% Val. The standardized digestibility values for these amino acids were 84, 92, 87, 86, and 87%, respectively. The TMEn of the Still Pro sample was determined to be 3,372 kcal/kg DM. In conclusion, the results of this study indicate that the nutritional value of the 3 types of novel ethanol coproducts is equal to or superior to the nutritional values generally reported for conventional DDGS.

Transcriptome Profiling and Genome-Wide Association Studies Reveal GSTs and Other Defense Genes Involved in Multiple Signaling Pathways Induced by Herbicide Safener in Grain Sorghum

Herbicide safeners protect cereal crops from herbicide injury by inducing genes and proteins involved in detoxification reactions, such as glutathione S-transferases (GSTs) and cytochrome P450s (P450s). Only a few studies have characterized gene or protein expression profiles for investigating plant responses to safener treatment in cereal crops, and most transcriptome analyses in response to safener treatments have been conducted in dicot model species that are not protected by safener from herbicide injury. In this study, three different approaches were utilized in grain sorghum (Sorghum bicolor (L.) Moench) to investigate mechanisms involved in safener-regulated signaling pathways. An initial transcriptome analysis was performed to examine global gene expression in etiolated shoot tissues of hybrid grain sorghum following treatment with the sorghum safener, fluxofenim. Most upregulated transcripts encoded detoxification enzymes, including P450s, GSTs, and UDP-dependent glucosyltransferases (UGTs). Interestingly, several of these upregulated transcripts are similar to genes involved with the biosynthesis and recycling/catabolism of dhurrin, an important chemical defense compound, in these seedling tissues. Secondly, 761 diverse sorghum inbred lines were evaluated in a genome-wide association study (GWAS) to determine key molecular-genetic factors governing safener-mediated signaling mechanisms and/or herbicide detoxification. GWAS revealed a significant single nucleotide polymorphism (SNP) associated with safener-induced response on chromosome 9, located within a phi-class SbGST gene and about 15-kb from a different phi-class SbGST. Lastly, the expression of these two candidate SbGSTs was quantified in etiolated shoot tissues of sorghum inbred BTx623 in response to fluxofenim treatment. SbGSTF1 and SbGSTF2 transcripts increased within 12-hr after fluxofenim treatment but the level of safener-induced expression differed between the two genes. In addition to identifying specific GSTs potentially involved in the safener-mediated detoxification pathway, this research elucidates a new direction for studying both constitutive and inducible mechanisms for chemical defense in cereal crop seedlings.

Maize Genomes to Fields: 2014 and 2015 field season genotype, phenotype, environment, and inbred ear image datasets

OBJECTIVES

Crop improvement relies on analysis of phenotypic, genotypic, and environmental data. Given large, well-integrated, multi-year datasets, diverse queries can be made: Which lines perform best in hot, dry environments? Which alleles of specific genes are required for optimal performance in each environment? Such datasets also can be leveraged to predict cultivar performance, even in uncharacterized environments. The maize Genomes to Fields (G2F) Initiative is a multi-institutional organization of scientists working to generate and analyze such datasets from existing, publicly available inbred lines and hybrids. G2F's genotype by environment project has released 2014 and 2015 datasets to the public, with 2016 and 2017 collected and soon to be made available.

DATA DESCRIPTION

Datasets include DNA sequences; traditional phenotype descriptions, as well as detailed ear, cob, and kernel phenotypes quantified by image analysis; weather station measurements; and soil characterizations by site. Data are released as comma separated value spreadsheets accompanied by extensive README text descriptions. For genotypic and phenotypic data, both raw data and a version with outliers removed are reported. For weather data, two versions are reported: a full dataset calibrated against nearby National Weather Service sites and a second calibrated set with outliers and apparent artifacts removed.

Transposable elements, mRNA expression level and strand-specificity of small RNAs are associated with non-additive inheritance of gene expression in hybrid plants

BACKGROUND

Gene expression inheritance patterns in Arabidopsis hybrid plants were investigated for correlation with the presence of transposable elements (TEs) and small RNA profile.

RESULTS

The presence of TEs in a gene and the expression of small RNA matching a gene were both found to be associated with non-additive mRNA inheritance patterns in hybrids. Expression levels below mid-parent values in the hybrids were associated with low mRNA expression in parents, with the presence of small RNA from both strands, and with the presence of TEs. High-parent dominance of mRNA levels was found to be associated with high parental mRNA expression levels, the absence of TEs, and for some genes, with small RNA fragments that are predominantly from the sense strand. These small RNAs exhibit a broader size distribution than siRNA and reduced nucleotide end bias, which are consistent with an origin from degraded mRNA. Thus, increased as well as decreased gene expression in hybrids relative to the parental mean is associated with gene expression levels, TE presence and small RNA fragments with differing characteristics.

CONCLUSIONS

The data presented here is consistent with a role for differential mRNA decay kinetics as one mechanism contributing to high-parent dominance in gene expression. Our evidence is also consistent with trans repression by siRNA and TEs as the cause of low-parent dominance.

Transcriptional responses indicate maintenance of photosynthetic proteins as key to the exceptional chilling tolerance of C4 photosynthesis in Miscanthus × giganteus

Miscanthus × giganteus is exceptional among C4 plants in its ability to acclimate to chilling (≤14 °C) and maintain a high photosynthetic capacity, in sharp contrast to maize, leading to very high productivity even in cool temperate climates. To identify the mechanisms that underlie this acclimation, RNA was isolated from M × giganteus leaves in chilling and nonchilling conditions and hybridized to microarrays developed for its close relative Zea mays. Among 21 000 array probes that yielded robust signals, 723 showed significant expression change under chilling. Approximately half of these were for annotated genes. Thirty genes associated with chloroplast membrane function were all upregulated. Increases in transcripts for the lhcb5 (chlorophyll a/b-binding protein CP26), ndhF (NADH dehydrogenase F, chloroplast), atpA (ATP synthase alpha subunit), psbA (D1), petA (cytochrome f), and lhcb4 (chlorophyll a/b-binding protein CP29), relative to housekeeping genes in M. × giganteus, were confirmed by quantitative reverse-transcription PCR. In contrast, psbo1, lhcb5, psbA, and lhcb4 were all significantly decreased in Z. mays after 14 days of chilling. Western blot analysis of the D1 protein and LHCII type II chlorophyll a/b-binding protein also showed significant increases in M. × giganteus during chilling and significant decreases in Z. mays. Compared to other C4 species, M. × giganteus grown in chilling conditions appears to counteract the loss of photosynthetic proteins and proteins protecting photosystem II typically observed in other species by increasing mRNA levels for their synthesis.

A detailed gene expression study of the Miscanthus genus reveals changes in the transcriptome associated with the rejuvenation of spring rhizomes

BACKGROUND

The Miscanthus genus of perennial C4 grasses contains promising biofuel crops for temperate climates. However, few genomic resources exist for Miscanthus, which limits understanding of its interesting biology and future genetic improvement. A comprehensive catalog of expressed sequences were generated from a variety of Miscanthus species and tissue types, with an emphasis on characterizing gene expression changes in spring compared to fall rhizomes.

RESULTS

Illumina short read sequencing technology was used to produce transcriptome sequences from different tissues and organs during distinct developmental stages for multiple Miscanthus species, including Miscanthus sinensis, Miscanthus sacchariflorus, and their interspecific hybrid Miscanthus × giganteus. More than fifty billion base-pairs of Miscanthus transcript sequence were produced. Overall, 26,230 Sorghum gene models (i.e., ~ 96% of predicted Sorghum genes) had at least five Miscanthus reads mapped to them, suggesting that a large portion of the Miscanthus transcriptome is represented in this dataset. The Miscanthus × giganteus data was used to identify genes preferentially expressed in a single tissue, such as the spring rhizome, using Sorghum bicolor as a reference. Quantitative real-time PCR was used to verify examples of preferential expression predicted via RNA-Seq. Contiguous consensus transcript sequences were assembled for each species and annotated using InterProScan. Sequences from the assembled transcriptome were used to amplify genomic segments from a doubled haploid Miscanthus sinensis and from Miscanthus × giganteus to further disentangle the allelic and paralogous variations in genes.

CONCLUSIONS

This large expressed sequence tag collection creates a valuable resource for the study of Miscanthus biology by providing detailed gene sequence information and tissue preferred expression patterns. We have successfully generated a database of transcriptome assemblies and demonstrated its use in the study of genes of interest. Analysis of gene expression profiles revealed biological pathways that exhibit altered regulation in spring compared to fall rhizomes, which are consistent with their different physiological functions. The expression profiles of the subterranean rhizome provides a better understanding of the biological activities of the underground stem structures that are essentials for perenniality and the storage or remobilization of carbon and nutrient resources.

A framework genetic map for Miscanthus sinensis from RNAseq-based markers shows recent tetraploidy

Background

Miscanthus (subtribe Saccharinae, tribe Andropogoneae, family Poaceae) is a genus of temperate perennial C4 grasses whose high biomass production makes it, along with its close relatives sugarcane and sorghum, attractive as a biofuel feedstock. The base chromosome number of Miscanthus (x = 19) is different from that of other Saccharinae and approximately twice that of the related Sorghum bicolor (x = 10), suggesting large-scale duplications may have occurred in recent ancestors of Miscanthus. Owing to the complexity of the Miscanthus genome and the complications of self-incompatibility, a complete genetic map with a high density of markers has not yet been developed.

Results

We used deep transcriptome sequencing (RNAseq) from two M. sinensis accessions to define 1536 single nucleotide variants (SNVs) for a GoldenGate™ genotyping array, and found that simple sequence repeat (SSR) markers defined in sugarcane are often informative in M. sinensis. A total of 658 SNP and 210 SSR markers were validated via segregation in a full sibling F1 mapping population. Using 221 progeny from this mapping population, we constructed a genetic map for M. sinensis that resolves into 19 linkage groups, the haploid chromosome number expected from cytological evidence. Comparative genomic analysis documents a genome-wide duplication in Miscanthus relative to Sorghum bicolor, with subsequent insertional fusion of a pair of chromosomes. The utility of the map is confirmed by the identification of two paralogous C4-pyruvate, phosphate dikinase (C4-PPDK) loci in Miscanthus, at positions syntenic to the single orthologous gene in Sorghum.

Conclusions

The genus Miscanthus experienced an ancestral tetraploidy and chromosome fusion prior to its diversification, but after its divergence from the closely related sugarcane clade. The recent timing of this tetraploidy complicates discovery and mapping of genetic markers for Miscanthus species, since alleles and fixed differences between paralogs are comparable. These difficulties can be overcome by careful analysis of segregation patterns in a mapping population and genotyping of doubled haploids. The genetic map for Miscanthus will be useful in biological discovery and breeding efforts to improve this emerging biofuel crop, and also provide a valuable resource for understanding genomic responses to tetraploidy and chromosome fusion.

The inheritance pattern of 24 nt siRNA clusters in arabidopsis hybrids is influenced by proximity to transposable elements

Hybrids often display increased size and growth, and thus are widely cultivated in agriculture and horticulture. Recent discoveries demonstrating the important regulatory roles of small RNAs have greatly improved our understanding of many basic biological questions, and could illuminate the molecular basis for the enhanced growth and size of hybrid plants. We profiled small RNAs by deep sequencing to characterize the inheritance patterns of small RNA levels in reciprocal hybrids of two Arabidopsis thaliana accessions, Columbia and Landsberg erecta. We find 24-nt siRNAs predominate among those small RNAs that are differentially expressed between the parents. Following hybridization, the transposable element (TE)-derived siRNAs are often inherited in an additive manner, whereas siRNAs associated with protein-coding genes are often down-regulated in hybrids to the levels observed for the parent with lower relative siRNA levels. Among the protein-coding genes that exhibit this pattern, genes that function in pathogen defense, abiotic stress tolerance, and secondary metabolism are significantly enriched. Small RNA clusters from protein-coding genes where a TE is present within one kilobase show a different predominant inheritance pattern (additive) from those that do not (low-parent dominance). Thus, down-regulation in the form of low-parent dominance is likely the default pattern of inheritance for genic siRNA, and a different inheritance mechanism for TE siRNA is suggested.

Microcollinearity between autopolyploid sugarcane and diploid sorghum genomes

Background

Sugarcane (Saccharum spp.) has become an increasingly important crop for its leading role in biofuel production. The high sugar content species S. officinarum is an octoploid without known diploid or tetraploid progenitors. Commercial sugarcane cultivars are hybrids between S. officinarum and wild species S. spontaneum with ploidy at ~12×. The complex autopolyploid sugarcane genome has not been characterized at the DNA sequence level.

Results

The microsynteny between sugarcane and sorghum was assessed by comparing 454 pyrosequences of 20 sugarcane bacterial artificial chromosomes (BACs) with sorghum sequences. These 20 BACs were selected by hybridization of 1961 single copy sorghum overgo probes to the sugarcane BAC library with one sugarcane BAC corresponding to each of the 20 sorghum chromosome arms. The genic regions of the sugarcane BACs shared an average of 95.2% sequence identity with sorghum, and the sorghum genome was used as a template to order sequence contigs covering 78.2% of the 20 BAC sequences. About 53.1% of the sugarcane BAC sequences are aligned with sorghum sequence. The unaligned regions contain non-coding and repetitive sequences. Within the aligned sequences, 209 genes were annotated in sugarcane and 202 in sorghum. Seventeen genes appeared to be sugarcane-specific and all validated by sugarcane ESTs, while 12 appeared sorghum-specific but only one validated by sorghum ESTs. Twelve of the 17 sugarcane-specific genes have no match in the non-redundant protein database in GenBank, perhaps encoding proteins for sugarcane-specific processes. The sorghum orthologous regions appeared to have expanded relative to sugarcane, mostly by the increase of retrotransposons.

Conclusions

The sugarcane and sorghum genomes are mostly collinear in the genic regions, and the sorghum genome can be used as a template for assembling much of the genic DNA of the autopolyploid sugarcane genome. The comparable gene density between sugarcane BACs and corresponding sorghum sequences defied the notion that polyploidy species might have faster pace of gene loss due to the redundancy of multiple alleles at each locus.

Genomic and small RNA sequencing of Miscanthus x giganteus shows the utility of sorghum as a reference genome sequence for Andropogoneae grasses

Genomic data together with sequencing of tissue specific small RNA libraries reveals insights into the genome content, small RNA repertoire and evolutionary origins of the grass Miscanthus × giganteus.

Background

Miscanthus × giganteus (Mxg) is a perennial grass that produces superior biomass yields in temperate environments. The essentially uncharacterized triploid genome (3n = 57, x = 19) of Mxg is likely critical for the rapid growth of this vegetatively propagated interspecific hybrid.

Results

A survey of the complex Mxg genome was conducted using 454 pyrosequencing of genomic DNA and Illumina sequencing-by-synthesis of small RNA. We found that the coding fraction of the Mxg genome has a high level of sequence identity to that of other grasses. Highly repetitive sequences representing the great majority of the Mxg genome were predicted using non-cognate assembly for de novo repeat detection. Twelve abundant families of repeat were observed, with those related to either transposons or centromeric repeats likely to comprise over 95% of the genome. Comparisons of abundant repeat sequences to a small RNA survey of three Mxg organs (leaf, rhizome, inflorescence) revealed that the majority of observed 24-nucleotide small RNAs are derived from these repetitive sequences. We show that high-copy-number repeats match more of the small RNA, even when the amount of the repeat sequence in the genome is accounted for.

Conclusions

We show that major repeats are present within the triploid Mxg genome and are actively producing small RNAs. We also confirm the hypothesized origins of Mxg, and suggest that while the repeat content of Mxg differs from sorghum, the sorghum genome is likely to be of utility in the assembly of a gene-space sequence of Mxg.

Cool C4 photosynthesis: pyruvate Pi dikinase expression and activity corresponds to the exceptional cold tolerance of carbon assimilation in Miscanthus x giganteus

The bioenergy feedstock grass Miscanthus x giganteus is exceptional among C(4) species for its high productivity in cold climates. It can maintain photosynthetically active leaves at temperatures 6 degrees C below the minimum for maize (Zea mays), which allows it a longer growing season in cool climates. Understanding the basis for this difference between these two closely related plants may be critical in adapting maize to colder weather. When M. x giganteus and maize grown at 25 degrees C were transferred to 14 degrees C, light-saturated CO(2) assimilation and quantum yield of photosystem II declined by 30% and 40%, respectively, in the first 48 h in these two species. The decline continued in maize but arrested and then recovered partially in M. x giganteus. Within 24 h of the temperature transition, the pyruvate phosphate dikinase (PPDK) protein content per leaf area transiently declined in M. x giganteus but then steadily increased, such that after 7 d the enzyme content was significantly higher than in leaves growing in 25 degrees C. By contrast it declined throughout the chilling period in maize leaves. Rubisco levels remained constant in M. x giganteus but declined in maize. Consistent with increased PPDK protein content, the extractable PPDK activity per unit leaf area (V(max)(,ppdk)) in cold-grown M. x giganteus leaves was higher than in warm-grown leaves, while V(max,ppdk) was lower in cold-grown than in warm-grown maize. The rate of light activation of PPDK was also slower in cold-grown maize than M. x giganteus. The energy of activation (E(a)) of extracted PPDK was lower in cold-grown than warm-grown M. x giganteus but not in maize. The specific activities and E(a) of purified recombinant PPDK from M. x giganteus and maize cloned into Escherichia coli were similar. The increase in PPDK protein in the M. x giganteus leaves corresponded to an increase in PPDK mRNA level. These results indicate that of the two enzymes known to limit C(4) photosynthesis, increase of PPDK, not Rubisco content, corresponds to the recovery and maintenance of photosynthetic capacity. Functionally, increased enzyme concentration is shown to increase stability of M. x giganteus PPDK at low temperature. The results suggest that increases in either PPDK RNA transcription and/or the stability of this RNA are important for the increase in PPDK protein content and activity in M. x giganteus under chilling conditions relative to maize.

Can the cold tolerance of C4 photosynthesis in Miscanthus x giganteus relative to Zea mays be explained by differences in activities and thermal properties of Rubisco?

The previous investigations show that the amount and activity of Rubisco appears the major limitation to effective C(4) photosynthesis at low temperatures. The chilling-tolerant and bioenergy feedstock species Miscanthus x giganteus (M. x giganteus) is exceptionally productive among C(4) grasses in cold climates. It is able to develop photosynthetically active leaves at temperatures 6 degrees C below the minimum for maize, and achieves a productivity even at 52 degrees N that exceeds that of the most productive C(3) crops at this latitude. This study investigates whether this unusual low temperature tolerance can be attributed to differences in the amount or kinetic properties of Rubisco relative to maize. An efficient protocol was developed to purify large amounts of functional Rubisco from C(4) leaves. The maximum carboxylation activities (V(max)), activation states, catalytic rates per active site (K(cat)) and activation energies (E(a)) of purified Rubisco and Rubisco in crude leaf extracts were determined for M. x giganteus grown at 14 degrees C and 25 degrees C, and maize grown at 25 degrees C. The sequences of M. x giganteus Rubisco small subunit mRNA are highly conserved, and 91% identical to those of maize. Although there were a few differences between the species in the translated protein sequences, there were no significant differences in the catalytic properties (V(max), K(cat), and E(a)) for purified Rubisco, nor was there any effect of growth temperature in M. x giganteus on these kinetic properties. Extracted activities were close to the observed rates of CO(2) assimilation by the leaves in vivo. On a leaf area basis the extracted activities and activation state of Rubisco did not differ significantly, either between the two species or between growth temperatures. The activation state of Rubisco in leaf extracts showed no significant difference between warm and cold-grown M. x giganteus. In total, these results suggest that the ability of M. x giganteus to be productive and maintain photosynthetically competent leaves at low temperature does not result from low temperature acclimation or adaptation of the catalytic properties of Rubisco.

microRNA172 down-regulates glossy15 to promote vegetative phase change in maize

Shoot development in many higher plant species is characterized by phase change, where meristems and organs transition from one set of identities to another. The transition from a juvenile to adult leaf identity in maize is regulated by the APETALA2-like gene glossy15 (gl15). We demonstrate here that increasing gl15 activity in transgenic maize not only increases the number of leaves expressing juvenile traits, but also delays the onset of reproductive development, indicating that gl15 plays a primary role in the maintenance of the juvenile phase. We also show that the accumulation of a maize microRNA homologous to miR172 increases during shoot development and mediates gl15 mRNA degradation. These data indicate that vegetative phase change in maize is regulated by the opposing actions of gl15 and miR172, with gl15 maintaining the juvenile phase and miR172 promoting the transition to the adult phase by down-regulation of gl15. Our results also suggest that the balance of activities between APETALA2-like genes and miR172 may be a general mechanism for regulating vegetative phase change in higher plants.

Amino acid metabolism in maize earshoots. Implications for assimilate preconditioning and nitrogen signaling

Nitrogen (N) is an essential requirement for kernel growth in maize (Zea mays); however, little is known about how N assimilates are metabolized in young earshoots during seed development. The objective of this study was to assess amino acid metabolism in cob and spikelet tissues during the critical 2 weeks following silking. Two maize hybrids were grown in the field for 2 years at two levels of supplemental N fertilizer (0 and 168 kg N/ha). The effects of the reproductive sink on cob N metabolism were examined by comparing pollinated to unpollinated earshoots. Earshoots were sampled at 2, 8, 14, and 18 d after silking; dissected into cob, spikelet, and/or pedicel and kernel fractions; then analyzed for amino acid profiles and key enzyme activities associated with amino acid metabolism. Major amino acids in the cob were glutamine (Gln), aspartic acid (Asp), asparagine (Asn), glutamate, and alanine. Gln concentrations dropped dramatically from 2 to 14 d after silking in both pollinated and unpollinated cobs, whereas all other measured amino acids accumulated over time in unpollinated spikelets and cobs, especially Asn. N supply had a variable effect on individual amino acid levels in young cobs and spikelets, with Asn being the most notably enhanced. We found that the cob performs significant enzymatic interconversions among Gln, alanine, Asp, and Asn during early reproductive development, which may precondition the N assimilate supply for sustained kernel growth. The measured amino acid profiles and enzymatic activities suggest that the Asn to Gln ratio in cobs may be part of a signal transduction pathway involving aspartate aminotransferase, Gln synthetase, and Asn synthetase to indicate plant N status for kernel development.

The maize macrohairless1 locus specifically promotes leaf blade macrohair initiation and responds to factors regulating leaf identity

The leaf surfaces of almost all plant species possess specialized epidermal cell types that form hairs or trichomes. Maize leaves produce three distinct types of hairs, the most prominent being the macrohairs that serve as a marker for adult leaf identity and may contribute to insect resistance. This report describes the maize macrohairless1 (mhl1) locus, which promotes macrohair initiation specifically in the leaf blade. Each of seven recessive mhl1 mutant alleles significantly reduces or eliminates macrohairs in the leaf blade. The mhl1 mutations block macrohair initiation rather than interfering with macrohair morphogenesis. Genetic mapping placed mhl1 within bin 4 on chromosome 9. A second independently segregating locus was found to partially suppress the mhl1 mutant phenotype in certain genetic backgrounds. Macrohair density was observed to increase during early adult vegetative development and then progressively decline, suggesting macrohair initiation frequency is affected by factors that act throughout shoot development. Genetic analyses demonstrated that mhl1 acts in the same pathway but downstream of factors that either promote or repress adult leaf identity. Thus, mhl1 plays a key role in integrating developmental programs that regulate leaf identity during shoot development with those that specify macrohair initiation within the leaf blade.

The maize O2 and PBF proteins act additively to promote transcription from storage protein gene promoters in rice endosperm cells

A transient expression assay system was employed to investigate the possible use of the maize Opaque 2 (O2) and prolamin box binding factor (PBF) proteins as transcriptional activators of rice and wheat storage protein gene promoters. When assayed in developing rice endosperm cells, either O2 or PBF alone could increase transcription from the promoter of the rice glutelin gene, Gt1. However, mutant forms of O2 and PBF that are defective in DNA binding could not. Co-transfection with both transcriptional activators resulted in an additive increase in transactivation of the Gt1 promoter. Co-bombardment of a Gt1::GUS construct with plasmids expressing the DNA binding domains of O2 and PBF in antisense orientation resulted in a decrease of GUS expression below background levels. Similar stimulatory and additive effects of O2 and PBF could be observed on the promoters from other storage protein genes including rice globulin (Glb), prolamins (RP6 and PG5a) and a wheat glutenin (Bx7). However, responsiveness of the promoters from non-storage protein genes like rice actin and CaMV 35S to O2 and PBF was insignificant. Our results indicate that the maize O2 and PBF proteins can act singly or additively as effective stimulators of heterologous storage protein promoters in developing rice endosperm cells. These data support the use of well-characterized transcription factors from maize as an effective means of increasing the expression level of recombinant proteins in developing rice seeds.

Maize selection passes the century mark: a unique resource for 21st century genomics

The Illinois Long-Term Selection Experiment for grain protein and oil concentration in maize (Zea mays) is the longest continuous genetics experiment in higher plants. A total of 103 cycles of selection have produced nine related populations that exhibit phenotypic extremes for grain composition and a host of correlated traits. The use of functional genomics tools in this unique genetic resource provides exciting opportunities not only to discover the genes that contribute to phenotypic differences but also to investigate issues such as the response of plant genomes to artificial selection, the genetic architecture of quantitative traits and the source of continued genetic variation within domesticated crop genomes.

Cold tolerance of C4 photosynthesis in Miscanthus x giganteus: adaptation in amounts and sequence of C4 photosynthetic enzymes

Field-grown Miscanthus x giganteus maintains high photosynthetic quantum yields and biomass productivity in cool temperate climates. It is related to maize (Zea mays) and uses the same NADP-malic enzyme C(4) pathway. This study tests the hypothesis that M. x giganteus, in contrast to maize, forms photosynthetically competent leaves at low temperatures with altered amounts of pyruvate orthophosphate dikinase (PPDK) and Rubisco or altered properties of PPDK. Both species were grown at 25 degrees C/20 degrees C or 14 degrees C/11 degrees C (day/night), and leaf photosynthesis was measured from 5 degrees C to 38 degrees C. Protein and steady-state transcript levels for Rubisco, PPDK, and phosphoenolpyruvate carboxylase were assessed and the sequence of C(4)-PPDK from M. x giganteus was compared with other C(4) species. Low temperature growth had no effect on photosynthesis in M. x giganteus, but decreased rates by 80% at all measurement temperatures in maize. Amounts and expression of phosphoenolpyruvate carboxylase were affected little by growth temperature in either species. However, PPDK and Rubisco large subunit decreased >50% and >30%, respectively, in cold-grown maize, whereas these levels remained unaffected by temperature in M. x giganteus. Differences in protein content in maize were not explained by differences in steady-state transcript levels. Several different M. x giganteus C(4)-PPDK cDNA sequences were found, but putative translated protein sequences did not show conservation of amino acids contributing to cold stability in Flaveria brownii C(4)-PPDK. The maintenance of PPDK and Rubisco large subunit amounts in M. x giganteus is consistent with the hypothesis that these proteins are critical to maintaining high rates of C(4) photosynthesis at low temperature.

Conserved noncoding sequences among cultivated cereal genomes identify candidate regulatory sequence elements and patterns of promoter evolution

Surveys for conserved noncoding sequences (CNS) among genes from monocot cereal species were conducted to assess the general properties of CNS in grass genomes and their correlation with known promoter regulatory elements. Initial comparisons of 11 orthologous maize-rice gene pairs found that previously defined regulatory motifs could be identified within short CNS but could not be distinguished reliably from random sequence matches. Among the different phylogenetic footprinting algorithms tested, the VISTA tool yielded the most informative alignments of noncoding sequence. VISTA was used to survey for CNS among all publicly available genomic sequences from maize, rice, wheat, barley, and sorghum, representing >300 gene comparisons. Comparisons of orthologous maize-rice and maize-sorghum gene pairs identified 20 bp as a minimal length criterion for a significant CNS among grass genes, with few such CNS found to be conserved across rice, maize, sorghum, and barley. The frequency and length of cereal CNS as well as nucleotide substitution rates within CNS were consistent with the known phylogenetic distances among the species compared. The implications of these findings for the evolution of cereal gene promoter sequences and the utility of using the nearly completed rice genome sequence to predict candidate regulatory elements in other cereal genes by phylogenetic footprinting are discussed.

Tandemly duplicated Safener-induced glutathione S-transferase genes from Triticum tauschii contribute to genome- and organ-specific expression in hexaploid wheat

Glutathione S-transferase (GST) gene expression was examined in several Triticum species, differing in genome constitution and ploidy level, to determine genome contribution to GST expression in cultivated, hexaploid bread wheat (Triticum aestivum). Two tandemly duplicated tau class GST genes (TtGSTU1 and TtGSTU2) were isolated from a single bacterial artificial chromosome clone in a library constructed from the diploid wheat and D genome progenitor to cultivated wheat, Triticum tauschii. The genes are very similar in genomic structure and their encoded proteins are 95% identical. Gene-specific reverse transcriptase-polymerase chain reaction analysis revealed differential transcript accumulation of TtGSTU1 and TtGSTU2 in roots and shoots. Expression of both genes was induced by herbicide safeners, 2,4-dichlorophenoxyacetic acid and abscisic acid, in the shoots of T. tauschii; however, expression of TtGSTU1 was always higher than TtGSTU2. In untreated seedlings, TtGSTU1 was expressed in both shoots and roots, whereas TtGSTU2 expression was only detected in roots. RNA gel-blot analysis of ditelosomic, aneuploid lines that are deficient for 6AS, 6BS, or 6DS chromosome arms of cultivated, hexaploid bread wheat showed differential genome contribution to safener-induced GST expression in shoots compared with roots. The GST genes from the D genome of hexaploid wheat contribute most to safener-induced expression in the shoots, whereas GSTs from the B and D genomes contribute to safener-induced expression in the roots.

A maize zinc-finger protein binds the prolamin box in zein gene promoters and interacts with the basic leucine zipper transcriptional activator Opaque2

The prolamin box (P-box) is a highly conserved 7-bp sequence element (5'-TGTAAAG-3') found in the promoters of many cereal seed storage protein genes. Nuclear factors from maize endosperm specifically interact with the P-box present in maize prolamin genes (zeins). The presence of the P-box in all zein gene promoters suggests that interactions between endosperm DNA binding proteins and the P-box may play an important role in the coordinate activation of zein gene expression during endosperm development. We have cloned an endosperm-specific maize cDNA, named prolamin-box binding factor (PBF), that encodes a member of the recently described Dof class of plant Cys2-Cys2 zinc-finger DNA binding proteins. When tested in gel shift assays, PBF exhibits the same sequence-specific binding to the P-box as factors present in maize endosperm nuclei. Additionally, PBF interacts in vitro with the basic leucine zipper protein Opaque2, a known transcriptional activator of zein gene expression whose target site lies 20 bp downstream of the P-box in the 22-kDa zein gene promoter. The isolation of the PBF gene provides an essential tool to further investigate the functional role of the highly conserved P-box in regulating cereal storage protein gene expression.

Glossy15, an APETALA2-like gene from maize that regulates leaf epidermal cell identity

Vegetative development in many plants progresses through distinct juvenile and adult phases. In maize, the transition from juvenile to adult shoot development affects a variety of leaf epidermal cell traits. These include epicuticular waxes, leaf hairs, and cell wall characteristics. Previous genetic and phenotypic analyses have shown that the maize Glossy15 (Gl15) gene is required for the expression of juvenile epidermal traits after leaf 2. We report here the molecular cloning of the Gl15 gene using a defective Suppressor-Mutator (dSpm) element insertion as a transposon-tag. Consistent with the gl15 mutant phenotype, the pattern of Gl15 mRNA expression was correlated with a juvenile leaf epidermal cell identity and was regulated by upstream factors such as Corngrass1. The Gl15 gene encodes a putative transcription factor with significant sequence similarity to the Arabidopsis regulatory genes APETALA2 and AINTEGUMENTA, which act primarily to regulate floral organ identity and ovule development. This finding expands the known functions of APETALA2-related genes to include the control of both vegetative and reproductive lateral organ identity and provides molecular support for the hypothesis that leaves and floral organs are related structures derived from a common growth plan.

Characterization of an immunoglobulin binding protein homolog in the maize floury-2 endosperm mutant

The maize b-70 protein is an endoplasmic reticulum protein overproduced in the floury-2 (fl2) endosperm mutant. The increase in b-70 levels in fl2 plants occurs during seed maturation and is endosperm specific. We have used amino acid sequence homology to identify b-70 as a homolog of mammalian immunoglobulin binding protein (BiP). Purified b-70 fractions contain two 75-kilodalton polypeptides with pl values of 5.3 and 5.4. Both 75-kilodalton polypeptides share several properties with BiP, including the ability to bind ATP and localization within the lumen of the endoplasmic reticulum. In addition, both b-70 polypeptides can be induced in maize cell cultures with tunicamycin treatment. Like BiP, the pl 5.3 form of b-70 is post-translationally modified by phosphorylation and ADP-ribosylation. However, modification of the pl 5.4 species was not detected in vitro or in vivo. Although the b-70 gene is unlinked to fl2, b-70 overproduction is positively correlated with the fl2 gene and is regulated at the mRNA level. In contrast, the fl2 allele negatively affects the accumulation of the major endosperm storage proteins. The physical similarity of b-70 to BiP and its association with abnormal protein accumulation in fl2 endoplasmic reticulum may reflect a biological function to mediate protein folding and assembly in maize endosperm.