Temporally resolved growth patterns reveal novel information about the polygenic nature of complex quantitative traits.

Plant height can be an indicator of plant health across environments and used to identify superior genotypes. Typically plant height is measured at a single timepoint when plants reach terminal height. Evaluating plant height using unoccupied aerial vehicles allows for measurements throughout the growing season, facilitating a better understanding of plant-environment interactions and the genetic basis of this complex trait. To assess variation throughout development, plant height data was collected from planting until terminal height at anthesis (14 flights 2018, 27 in 2019, 12 in 2020, and 11 in 2021) for a panel of ~500 diverse maize inbred lines. The percent variance explained in plant height throughout the season was significantly explained by genotype (9-48%), year (4-52%), and genotype-by-year interactions (14-36%) to varying extents throughout development. Genome-wide association studies revealed 717 significant single nucleotide polymorphisms associated with plant height and growth rate at different parts of the growing season specific to certain phases of vegetative growth. When plant height growth curves were compared to growth curves estimated from canopy cover, greater Fréchet distance stability was observed in plant height growth curves than for canopy cover. This indicated canopy cover may be more useful for understanding environmental modulation of overall plant growth and plant height better for understanding genotypic modulation of overall plant growth. This study demonstrated that substantial information can be gained from high temporal resolution data to understand how plants differentially interact with the environment and can enhance our understanding of the genetic basis of complex polygenic traits.

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

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

Description and functional analysis of the transcriptome from malting barley.

The present study aimed to establish an early model of the malting barley transcriptome, which describes the expression of genes and their ontologies, identify the period during malting with the largest dynamic shift in gene expression for future investigation, and to determine the expression patterns of all starch degrading enzyme genes relevant to the malting and brewing industry. Large dynamic increases in gene expression occurred early in malting with differential expressed genes enriched for cell wall and starch hydrolases amongst many malting related categories. Twenty-five of forty starch degrading enzyme genes were differentially expressed in the malting barley transcriptome including eleven α-amylase genes, six ß-amylase genes, three α-glucosidase genes, and all five starch debranching enzyme genes. Four new or novel α-amylase genes, one ß-amylase gene (Bmy3), three α-glucosidase genes, and two isoamylase genes had appreciable expression that requires further exploration into their potential relevance to the malting and brewing industry.

Transposable elements contribute to dynamic genome content in maize.

Transposable elements (TEs) are ubiquitous components of eukaryotic genomes and can create variation in genome organization and content. Most maize genomes are composed of TEs. We developed an approach to define shared and variable TE insertions across genome assemblies and applied this method to four maize genomes (B73, W22, Mo17 and PH207) with uniform structural annotations of TEs. Among these genomes we identified approximately 400 000 TEs that are polymorphic, encompassing 1.6 Gb of variable TE sequence. These polymorphic TEs include a combination of recent transposition events as well as deletions of older TEs. There are examples of polymorphic TEs within each of the superfamilies of TEs and they are found distributed across the genome, including in regions of recent shared ancestry among individuals. There are many examples of polymorphic TEs within or near maize genes. In addition, there are 2380 gene annotations in the B73 genome that are located within variable TEs, providing evidence for the role of TEs in contributing to the substantial differences in annotated gene content among these genotypes. TEs are highly variable in our survey of four temperate maize genomes, highlighting the major contribution of TEs in driving variation in genome organization and gene content. OPEN RESEARCH BADGES: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://github.com/SNAnderson/maizeTE_variation; https://mcstitzer.github.io/maize_TEs.

Localized Genetic and Phenotypic Diversity of Xanthomonas translucens Associated With Bacterial Leaf Streak on Wheat and Barley in Minnesota.

Bacterial leaf streak (BLS) of wheat and barley has been a disease of increasing concern in the Upper Midwest over the past decade. In this study, intra- and interfield genetic and pathogenic diversity of bacteria causing BLS in Minnesota was evaluated. In 2015, 89 strains were isolated from 100 leaf samples collected from two wheat and two barley fields naturally infected with BLS. Virulence assays and multilocus sequence alignments of four housekeeping genes supported pathovar identifications. All wheat strains were pathogenic on wheat and barley and belonged to the same lineage as the Xanthomonas translucens pv. undulosa-type strain. All barley strains were pathogenic on barley but not on wheat. Three lineages of barley strains were detected. The frequency and number of sequence types of each pathovar varied within and between fields. A significant population variance was detected between populations of X. translucens pv. undulosa collected from different wheat fields. Population stratification of X. translucens pv. translucens was not detected. Significant differences in virulence were detected among three dominant sequence types of X. translucens pv. undulosa but not those of X. translucens pv. translucens. Field trials with wheat and barley plants inoculated with strains of known sequence type and virulence did not detect significant race structures within either pathovar. Knowledge of virulence, sequence types, and population structures of X. translucens on wheat and barley can support studies on plant-bacterial interactions and breeding for BLS disease resistance.

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

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

Natural variation for gene expression responses to abiotic stress in maize.

Plants respond to abiotic stress through a variety of physiological, biochemical, and transcriptional mechanisms. Many genes exhibit altered levels of expression in response to abiotic stress, which requires concerted action of both cis- and trans-regulatory features. In order to study the variability in transcriptome response to abiotic stress, RNA sequencing was performed using 14-day-old maize seedlings of inbreds B73, Mo17, Oh43, PH207 and B37 under control, cold and heat conditions. Large numbers of genes that responded differentially to stress between parental inbred lines were identified. RNA sequencing was also performed on similar tissues of the F1 hybrids produced by crossing B73 and each of the three other inbred lines. By evaluating allele-specific transcript abundance in the F1 hybrids, we were able to measure the abundance of cis- and trans-regulatory variation between genotypes for both steady-state and stress-responsive expression differences. Although examples of trans-regulatory variation were observed, cis-regulatory variation was more common for both steady-state and stress-responsive expression differences. The genes with cis-allelic variation for response to cold or heat stress provided an opportunity to study the basis for regulatory diversity.

RNA-directed DNA methylation enforces boundaries between heterochromatin and euchromatin in the maize genome.

The maize genome is relatively large (∼ 2.3 Gb) and has a complex organization of interspersed genes and transposable elements, which necessitates frequent boundaries between different types of chromatin. The examination of maize genes and conserved noncoding sequences revealed that many of these are flanked by regions of elevated asymmetric CHH (where H is A, C, or T) methylation (termed mCHH islands). These mCHH islands are quite short (∼ 100 bp), are enriched near active genes, and often occur at the edge of the transposon that is located nearest to genes. The analysis of DNA methylation in other sequence contexts and several chromatin modifications revealed that mCHH islands mark the transition from heterochromatin-associated modifications to euchromatin-associated modifications. The presence of an mCHH island is fairly consistent in several distinct tissues that were surveyed but shows some variation among different haplotypes. The presence of insertion/deletions in promoters often influences the presence and position of an mCHH island. The mCHH islands are dependent upon RNA-directed DNA methylation activities and are lost in mop1 and mop3 mutants, but the nearby genes rarely exhibit altered expression levels. Instead, loss of an mCHH island is often accompanied by additional loss of DNA methylation in CG and CHG contexts associated with heterochromatin in nearby transposons. This suggests that mCHH islands and RNA-directed DNA methylation near maize genes may act to preserve the silencing of transposons from activity of nearby genes.

Genomic limitations to RNA sequencing expression profiling.

The field of genomics has grown rapidly with the advent of massively parallel sequencing technologies, allowing for novel biological insights with regards to genomic, transcriptomic, and epigenomic variation. One widely utilized application of high-throughput sequencing is transcriptional profiling using RNA sequencing (RNAseq). Understanding the limitations of a technology is critical for accurate biological interpretations, and clear interpretation of RNAseq data can be difficult in species with complex genomes. To understand the limitations of accurate profiling of expression levels we simulated RNAseq reads from annotated gene models in several plant species including Arabidopsis, brachypodium, maize, potato, rice, soybean, and tomato. The simulated reads were aligned using various parameters such as unique versus multiple read alignments. This allowed the identification of genes recalcitrant to RNAseq analyses by having over- and/or under-estimated expression levels. In maize, over 25% of genes deviated by more than 20% from the expected count values, suggesting the need for cautious interpretation of RNAseq data for certain genes. The reasons identified for deviation from expected expression varied between species due to differences in genome structure including, but not limited to, genes encoding short transcripts, overlapping gene models, and gene family size. Utilizing existing empirical datasets we demonstrate the potential for biological misinterpretation resulting from inclusion of 'flagged genes' in analyses. While RNAseq is a powerful tool for understanding biology, there are limitations to this technology that need to be understood in order to improve our biological interpretations.

Conservation and purifying selection of transcribed genes located in a rice centromere.

Recombination is strongly suppressed in centromeric regions. In chromosomal regions with suppressed recombination, deleterious mutations can easily accumulate and cause degeneration of genes and genomes. Surprisingly, the centromere of chromosome8 (Cen8) of rice (Oryza sativa) contains several transcribed genes. However, it remains unclear as to what selective forces drive the evolution and existence of transcribed genes in Cen8. Sequencing of orthologous Cen8 regions from two additional Oryza species, Oryza glaberrima and Oryza brachyantha, which diverged from O. sativa 1 and 10 million years ago, respectively, revealed a set of seven transcribed Cen8 genes conserved across all three species. Chromatin immunoprecipitation analysis with the centromere-specific histone CENH3 confirmed that the sequenced orthologous regions are part of the functional centromere. All seven Cen8 genes have undergone purifying selection, representing a striking phenomenon of active gene survival within a recombination-free zone over a long evolutionary time. The coding sequences of the Cen8 genes showed sequence divergence and mutation rates that were significantly reduced from those of genes located on the chromosome arms. This suggests that Oryza has a mechanism to maintain the fidelity and functionality of Cen8 genes, even when embedded in a sea of repetitive sequences and transposable elements.

Optimized Methods for Applying and Assessing Heat, Drought, and Nutrient Stress of Maize Seedlings in Controlled Environment Experiments.

Maize (Zea mays), also known as corn, is an important crop that plays a crucial role in global agriculture. The economic uses of maize are numerous, including for food, feed, fiber, and fuel. It has had a significant historical importance in research as well, with important discoveries made in maize regarding plant domestication, transposons, heterosis, genomics, and epigenetics. Unfortunately, environmental stresses cause substantial yield loss to maize crops each year. Yield losses are predicted to increase in future climate scenarios, posing a threat to food security and other sectors of the global economy. Developing efficient methods to study maize abiotic stress responses is a crucial step toward a more resilient and productive agricultural system. This review describes the importance of and methods for studying the effects of heat, drought, and nutrient deficiency on early developmental stages of maize grown in controlled environments. Studying the early effects of environmental stressors in controlled environments allows researchers to work with a variety of environmental conditions with low environmental variance, which can inform future field-based research. We highlight the current knowledge of physiological responses of maize to heat, drought, and nutrient stress; remaining knowledge gaps and challenges; and information on how standardized protocols can address these issues.

Maize Abiotic Stress Treatments in Controlled Environments.

Maize (Zea mays) is one of the world's most important crops, providing food for humans and livestock and serving as a bioenergy source. Climate change and the resulting abiotic stressors in the field reduce crop yields, threatening food security and the global economy. Water deficit (i.e., drought), heat, and insufficient nutrients (e.g., nitrogen and phosphorus) are major environmental stressors that affect maize yields, and impact growth and development at all stages of the plant life cycle. Understanding the biological processes underlying these responses in maize has the potential to increase yields in the face of abiotic stress. Optimizing individual or combined abiotic stress treatments in controlled environments reduces potential noise in data collection that can be present under less controlled growth conditions. Here, we describe methods and conditions for controlled abiotic stress treatments and associated controls during early vegetative growth of maize, conducted in greenhouses or growth chambers. This includes the environmental conditions, equipment, soil preparation, and intensity and duration of heat, drought, nitrogen deficiency, and phosphorous deficiency. Controlled experiments at early growth stages are informative for future in-field studies that require greater labor and inputs, saving researchers time and growing space, and thus research funds, before testing plants across later stages of development. We suggest that stress treatments be severe enough to result in a measurable phenotype, but not so severe that all plants die prior to sample collection. This protocol is designed to set important standards for replicable research in maize.

Wheat-Net: An Automatic Dense Wheat Spike Segmentation Method Based on an Optimized Hybrid Task Cascade Model.

Precise segmentation of wheat spikes from a complex background is necessary for obtaining image-based phenotypic information of wheat traits such as yield estimation and spike morphology. A new instance segmentation method based on a Hybrid Task Cascade model was proposed to solve the wheat spike detection problem with improved detection results. In this study, wheat images were collected from fields where the environment varied both spatially and temporally. Res2Net50 was adopted as a backbone network, combined with multi-scale training, deformable convolutional networks, and Generic ROI Extractor for rich feature learning. The proposed methods were trained and validated, and the average precision (AP) obtained for the bounding box and mask was 0.904 and 0.907, respectively, and the accuracy for wheat spike counting was 99.29%. Comprehensive empirical analyses revealed that our method (Wheat-Net) performed well on challenging field-based datasets with mixed qualities, particularly those with various backgrounds and wheat spike adjacence/occlusion. These results provide evidence for dense wheat spike detection capabilities with masking, which is useful for not only wheat yield estimation but also spike morphology assessments.

Hyperspectral imaging and improved feature variable selection for automated determination of deoxynivalenol in various genetic lines of barley kernels for resistance screening.

Fusarium head blight (FHB), a fungus disease of small grain cereal crops, results in reduced yields and diminished value of harvested grain due to the presence of deoxynivalenol (DON), a mycotoxin produced by the causal pathogen Fusarium graminearum. DON and other tricothecene mycotoxins pose serious health risks to both humans and livestock, especially swine. Due to these health concerns, barley used for malting, food or feed is routinely assayed for DON levels. Various methods are available for assaying DON levels in grain samples including enzyme-linked immunosorbent assay (ELISA) and gas chromatography-mass spectrometry (GC-MS). ELISA and GC-MS are very accurate; however, assaying grain samples by these techniques are laborious, expensive and destructive. In this study, we explored the feasibility of using hyperspectral imaging (382-1030 nm) to develop a rapid and non-destructive protocol for assaying DON in barley kernels. Samples of 888 and 116 from various genetic lines were selected for calibration and prediction. Full-wavelength locally weighted partial least squares regression (LWPLSR) achieved high accuracy with the coefficient of determination in prediction (R2P) of 0.728 and root mean square error of prediction (RMSEP) of 3.802. Competitive adaptive reweighted sampling (CARS) was used to choose potential feature wavelengths, and these selected variables were further optimized using the iterative selection of successive projections algorithm (ISSPA). The CARS-ISSPA-LWPLSR model developed using 7 feature variables yielded R2P of 0.680 and RMSEP of 4.213 in DON content prediction. Based on the 7 wavelengths selected by CARS-ISSPA, partial least square discriminant analysis (PLSDA) discriminated barley kernels having lower DON (less than1.25 mg/kg) levels from those with higher levels (including 1.25-3 mg/kg, 3-5 mg/kg, and 5-10 mg/kg), with Matthews correlation coefficient in cross-validation (M-RCV) of as high as 0.931. The results demonstrate that hyperspectral imaging have potential for accelerating non-destructive DON assays of barley samples.

Identification of Candidate Susceptibility Genes to Puccinia graminis f. sp. tritici in Wheat.

Wheat stem rust disease caused by Puccinia graminis f. sp. tritici (Pgt) is a global threat to wheat production. Fast evolving populations of Pgt limit the efficacy of plant genetic resistance and constrain disease management strategies. Understanding molecular mechanisms that lead to rust infection and disease susceptibility could deliver novel strategies to deploy crop resistance through genetic loss of disease susceptibility. We used comparative transcriptome-based and orthology-guided approaches to characterize gene expression changes associated with Pgt infection in susceptible and resistant Triticum aestivum genotypes as well as the non-host Brachypodium distachyon. We targeted our analysis to genes with differential expression in T. aestivum and genes suppressed or not affected in B. distachyon and report several processes potentially linked to susceptibility to Pgt, such as cell death suppression and impairment of photosynthesis. We complemented our approach with a gene co-expression network analysis to identify wheat targets to deliver resistance to Pgt through removal or modification of putative susceptibility genes.

Lineage-specific adaptive evolution of the centromeric protein CENH3 in diploid and allotetraploid Oryza species.

Centromeres in eukaryotic species are defined by the presence of a centromere-specific histone H3 variant, CENH3. CENH3 plays a key role in recruiting other centromeric proteins; thus, it is the central component in kinetochore formation and centromere function. The CENH3 proteins in several plant and animal species were found to be under positive selection, which was hypothesized to respond to the rapid changing of the repetitive DNA sequences associated with the centromeres. Here, we report the expression and evolution of the CenH3 genes in two allotetraploid rice species as well as their representative diploid progenitor species. Both copies of the CenH3 genes were transcribed in the two allotetraploid species and showed a nonpreferential expression pattern. Contrasting positive and stabilizing selection of the CenH3 genes was associated with different diploid Oryza species. This lineage-specific adaptive evolution of CENH3 was maintained in the two allotetraploid species. Thus, we demonstrate that the allopolyploidization events did not alter the expression or evolutionary patterns of the CenH3 genes in the Oryza species.

Towards Improved Molecular Identification Tools in Fine Fescue (Festuca L., Poaceae) Turfgrasses: Nuclear Genome Size, Ploidy, and Chloroplast Genome Sequencing.

Fine fescues (Festuca L., Poaceae) are turfgrass species that perform well in low-input environments. Based on morphological characteristics, the most commonly-utilized fine fescues are divided into five taxa: three are subspecies within F. rubra L. and the remaining two are treated as species within the F. ovina L. complex. Morphologically, these five taxa are very similar; both identification and classification of fine fescues remain challenging. In an effort to develop identification methods for fescues, we used flow cytometry to estimate genome size and ploidy level and sequenced the chloroplast genome of all five taxa. Fine fescue chloroplast genome sizes ranged from 133,331 to 133,841 bp and contained 113-114 genes. Phylogenetic relationship reconstruction using whole chloroplast genome sequences agreed with previous work based on morphology. Comparative genomics suggested unique repeat signatures for each fine fescue taxon that could potentially be used for marker development for taxon identification.

Dynamic Patterns of Transcript Abundance of Transposable Element Families in Maize.

Transposable Elements (TEs) are mobile elements that contribute the majority of DNA sequences in the maize genome. Due to their repetitive nature, genomic studies of TEs are complicated by the difficulty of properly attributing multi-mapped short reads to specific genomic loci. Here, we utilize a method to attribute RNA-seq reads to TE families rather than particular loci in order to characterize transcript abundance for TE families in the maize genome. We applied this method to assess per-family expression of transposable elements in >800 published RNA-seq libraries representing a range of maize development, genotypes, and hybrids. While a relatively small proportion of TE families are transcribed, expression is highly dynamic with most families exhibiting tissue-specific expression. A large number of TE families were specifically detected in pollen and endosperm, consistent with reproductive dynamics that maintain silencing of TEs in the germ line. We find that B73 transcript abundance is a poor predictor of TE expression in other genotypes and that transcript levels can differ even for shared TEs. Finally, by assessing recombinant inbred line and hybrid transcriptomes, complex patterns of TE transcript abundance across genotypes emerged. Taken together, this study reveals a dynamic contribution of TEs to maize transcriptomes.