Expression profiles of cell-wall related genes vary broadly between two common maize inbreds during stem development

The heat shock factor 20-HSF4-cellulose synthase A2 module regulates heat stress tolerance in maize

Temperature shapes the geographical distribution and behavior of plants. Understanding the regulatory mechanisms underlying the plant heat stress response is important for developing climate-resilient crops, including maize (Zea mays). To identify transcription factors (TFs) that may contribute to the maize heat stress response, we generated a dataset of short- and long-term transcriptome changes following a heat treatment time course in the inbred line B73. Co-expression network analysis highlighted several TFs, including the class B2a heat shock factor (HSF) ZmHSF20. Zmhsf20 mutant seedlings exhibited enhanced tolerance to heat stress. Furthermore, DNA affinity purification sequencing and Cleavage Under Targets and Tagmentation assays demonstrated that ZmHSF20 binds to the promoters of Cellulose synthase A2 (ZmCesA2) and three class A Hsf genes, including ZmHsf4, repressing their transcription. We showed that ZmCesA2 and ZmHSF4 promote the heat stress response, with ZmHSF4 directly activating ZmCesA2 transcription. In agreement with the transcriptome analysis, ZmHSF20 inhibited cellulose accumulation and repressed the expression of cell wall-related genes. Importantly, the Zmhsf20 Zmhsf4 double mutant exhibited decreased thermotolerance, placing ZmHsf4 downstream of ZmHsf20. We proposed an expanded model of the heat stress response in maize, whereby ZmHSF20 lowers seedling heat tolerance by repressing ZmHsf4 and ZmCesA2, thus balancing seedling growth and defense.

COMPILE: a GWAS computational pipeline for gene discovery in complex genomes

BACKGROUND

Genome-Wide Association Studies (GWAS) are used to identify genes and alleles that contribute to quantitative traits in large and genetically diverse populations. However, traits with complex genetic architectures create an enormous computational load for discovery of candidate genes with acceptable statistical certainty. We developed a streamlined computational pipeline for GWAS (COMPILE) to accelerate identification and annotation of candidate maize genes associated with a quantitative trait, and then matches maize genes to their closest rice and Arabidopsis homologs by sequence similarity.

RESULTS

COMPILE executed GWAS using a Mixed Linear Model that incorporated, without compression, recent advancements in population structure control, then linked significant Quantitative Trait Loci (QTL) to candidate genes and RNA regulatory elements contained in any genome. COMPILE was validated using published data to identify QTL associated with the traits of α-tocopherol biosynthesis and flowering time, and identified published candidate genes as well as additional genes and non-coding RNAs. We then applied COMPILE to 274 genotypes of the maize Goodman Association Panel to identify candidate loci contributing to resistance of maize stems to penetration by larvae of the European Corn Borer (Ostrinia nubilalis). Candidate genes included those that encode a gene of unknown function, WRKY and MYB-like transcriptional factors, receptor-kinase signaling, riboflavin synthesis, nucleotide-sugar interconversion, and prolyl hydroxylation. Expression of the gene of unknown function has been associated with pathogen stress in maize and in rice homologs closest in sequence identity.

CONCLUSIONS

The relative speed of data analysis using COMPILE allowed comparison of population size and compression. Limitations in population size and diversity are major constraints for a trait and are not overcome by increasing marker density. COMPILE is customizable and is readily adaptable for application to species with robust genomic and proteome databases.

Transcriptome Mining Provides Insights into Cell Wall Metabolism and Fiber Lignification in Agave tequilana Weber

Resilience of growing in arid and semiarid regions and a high capacity of accumulating sugar-rich biomass with low lignin percentages have placed Agave species as an emerging bioenergy crop. Although transcriptome sequencing of fiber-producing agave species has been explored, molecular bases that control wall cell biogenesis and metabolism in agave species are still poorly understood. Here, through RNAseq data mining, we reconstructed the cellulose biosynthesis pathway and the phenylpropanoid route producing lignin monomers in A. tequilana, and evaluated their expression patterns in silico and experimentally. Most of the orthologs retrieved showed differential expression levels when they were analyzed in different tissues with contrasting cellulose and lignin accumulation. Phylogenetic and structural motif analyses of putative CESA and CAD proteins allowed to identify those potentially involved with secondary cell wall formation. RT-qPCR assays revealed enhanced expression levels of AtqCAD5 and AtqCESA7 in parenchyma cells associated with extraxylary fibers, suggesting a mechanism of formation of sclerenchyma fibers in Agave similar to that reported for xylem cells in model eudicots. Overall, our results provide a framework for understanding molecular bases underlying cell wall biogenesis in Agave species studying mechanisms involving in leaf fiber development in monocots.

SLAG: A Program for Seeded Local Assembly of Genes in Complex Genomes

Although finished genomes have become more common, there is still a need for assemblies of individual genes or chromosomal regions when only unassembled reads are available. slag (Seeded Local Assembly of Genes) fulfils this need by performing iterative local assembly based on cycles of matching‐read retrieval with blast and assembly with cap3, phrap, spades, canu or unicycler. The target sequence can be nucleotide or protein. Read fragmentation allows slag to use phrap or cap3 to assemble long reads at lower coverage (e.g., 5×) than is possible with canu or unicycler. In simple, nonrepetitive genomes, a slag assembly can cover a whole chromosome, but in complex genomes the growth of target‐matching contigs is limited as additional reads are consumed by consensus contigs consisting of repetitive elements. Apart from genomic complexity, contig length and correctness depend on read length and accuracy. With pyrosequencing or Illumina reads, slag‐assembled contigs are accurate enough to allow design of PCR primers, while contigs assembled from Oxford Nanopore or pre‐HiFi Pacific Biosciences long reads are generally only accurate enough to design baiting sequences for further targeted sequencing. In an application with real reads, slag successfully extended sequences for four wheat genes, which were verified by cloning and Sanger sequencing of overlapping amplicons. slag is a robust alternative to atram2 for local assemblies, especially for read sets with less than 20× coverage. slag is freely available at https://github.com/cfcrane/SLAG.

Necrotic upper tips1 mimics heat and drought stress and encodes a protoxylem-specific transcription factor in maize

Maintaining sufficient water transport during flowering is essential for proper organ growth, fertilization, and yield. Water deficits that coincide with flowering result in leaf wilting, necrosis, tassel browning, and sterility, a stress condition known as "tassel blasting." We identified a mutant, necrotic upper tips1 (nut1), that mimics tassel blasting and drought stress and reveals the genetic mechanisms underlying these processes. The nut1 phenotype is evident only after the floral transition, and the mutants have difficulty moving water as shown by dye uptake and movement assays. These defects are correlated with reduced protoxylem vessel thickness that indirectly affects metaxylem cell wall integrity and function in the mutant. nut1 is caused by an Ac transposon insertion into the coding region of a unique NAC transcription factor within the VND clade of Arabidopsis NUT1 localizes to the developing protoxylem of root, stem, and leaf sheath, but not metaxylem, and its expression is induced by flowering. NUT1 downstream target genes function in cell wall biosynthesis, apoptosis, and maintenance of xylem cell wall thickness and strength. These results show that maintaining protoxylem vessel integrity during periods of high water movement requires the expression of specialized, dynamically regulated transcription factors within the vasculature.