Stability of Single-Parent Gene Expression Complementation in Maize Hybrids upon Water Deficit Stress

Molecular concepts to explain heterosis in crops

Heterosis describes the superior performance of hybrid plants compared with their genetically distinct parents and is a pillar of global food security. Here we review the current status of the molecular dissection of heterosis. We discuss how extensive intraspecific structural genomic variation between parental genotypes leads to heterosis by genetic complementation in hybrids. Moreover, we survey how global gene expression complementation contributes to heterosis by hundreds of additionally active genes in hybrids and how overdominant single genes mediate heterosis in several species. Furthermore, we highlight the prominent role of the microbiome in improving the performance of hybrids. Taken together, the molecular understanding of heterosis will pave the way to accelerate hybrid productivity and a more sustainable agriculture.

Dynamic patterns of gene expressional and regulatory variations in cotton heterosis

Purpose

Although the application of heterosis has significantly increased crop yield over the past century, the mechanisms underlying this phenomenon still remain obscure. Here, we applied transcriptome sequencing to unravel the impacts of parental expression differences and transcriptomic reprogramming in cotton heterosis.

Methods

A high-quality transcriptomic atlas covering 15 developmental stages and tissues was constructed for XZM2, an elite hybrid of upland cotton (Gossypium hirsutum L.), and its parental lines, CRI12 and J8891. This atlas allowed us to identify gene expression differences between the parents and to characterize the transcriptomic reprogramming that occurs in the hybrid.

Results

Our analysis revealed abundant gene expression differences between the parents, with pronounced tissue specificity; a total of 1,112 genes exhibited single-parent expression in at least one tissue. It also illuminated transcriptomic reprogramming in the hybrid XZM2, which included both additive and non-additive expression patterns. Coexpression networks between parents and hybrid constructed via weighted gene coexpression network analysis identified modules closely associated with fiber development. In particular, key regulatory hub genes involved in fiber development showed high-parent dominant or over dominant patterns in the hybrid, potentially driving the emergence of heterosis. Finally, high-depth resequencing data was generated and allele-specific expression patterns examined in the hybrid, enabling the dissection of cis- and trans-regulation contributions to the observed expression differences.

Conclusion

Parental transcriptional differences and transcriptomic reprogramming in the hybrid, especially the non-additive upregulation of key genes, play an important role in shaping heterosis. Collectively, these findings provide new insights into the molecular basis of heterosis in cotton.

Dynamic relationships among pathways producing hydrocarbons and fatty acids of maize silk cuticular waxes

The hydrophobic cuticle is the first line of defense between aerial portions of plants and the external environment. On maize (Zea mays L.) silks, the cuticular cutin matrix is infused with cuticular waxes, consisting of a homologous series of very long-chain fatty acids (VLCFAs), aldehydes, and hydrocarbons. Together with VLC fatty-acyl-CoAs (VLCFA-CoAs), these metabolites serve as precursors, intermediates, and end-products of the cuticular wax biosynthetic pathway. To deconvolute the potentially confounding impacts of the change in silk microenvironment and silk development on this pathway, we profiled cuticular waxes on the silks of the inbreds B73 and Mo17, and their reciprocal hybrids. Multivariate interrogation of these metabolite abundance data demonstrates that VLCFA-CoAs and total free VLCFAs are positively correlated with the cuticular wax metabolome, and this metabolome is primarily affected by changes in the silk microenvironment and plant genotype. Moreover, the genotype effect on the pathway explains the increased accumulation of cuticular hydrocarbons with a concomitant reduction in cuticular VLCFA accumulation on B73 silks, suggesting that the conversion of VLCFA-CoAs to hydrocarbons is more effective in B73 than Mo17. Statistical modeling of the ratios between cuticular hydrocarbons and cuticular VLCFAs reveals a significant role of precursor chain length in determining this ratio. This study establishes the complexity of the product-precursor relationships within the silk cuticular wax-producing network by dissecting both the impact of genotype and the allocation of VLCFA-CoA precursors to different biological processes and demonstrates that longer chain VLCFA-CoAs are preferentially utilized for hydrocarbon biosynthesis.

Transcriptome analysis reveals the key role of overdominant expression of photosynthetic and respiration-related genes in the formation of tobacco(Nicotiana tabacum L.) biomass heterosis

BACKGROUND

Leaves are the nutritional and economic organs of tobacco, and their biomass directly affects tobacco yield and the economic benefits of farmers. In the early stage, our research found that tobacco hybrids have more leaves and larger leaf areas, but the performance and formation reasons of biomass heterosis are not yet clear.

RESULTS

This study selected 5 parents with significant differences in tobacco biomass and paired them with hybrid varieties. It was found that tobacco hybrid varieties have a common biomass heterosis, and 45 days after transplantation is the key period for the formation of tobacco biomass heterosis; By analyzing the biomass heterosis of hybrids, Va116×GDH94 and its parents were selected for transcriptome analysis. 76.69% of the differentially expressed genes between Va116×GDH94 and its parents showed overdominant expression pattern, and these overdominant expression genes were significantly enriched in the biological processes of photosynthesis and TCA cycle; During the process of photosynthesis, the overdominant up-regulation of genes such as Lhc, Psa, and rbcl promotes the progress of photosynthesis, thereby increasing the accumulation of tobacco biomass; During the respiratory process, genes such as MDH, ACO, and OGDH are overedominantly down-regulated, inhibiting the TCA cycle and reducing substrate consumption in hybrid offspring; The photosynthetic characteristics of the hybrid and its parents were measured, and the net photosynthetic capacity of the hybrid was significantly higher than that of the parents.

CONCLUSION

These results indicate that the overdominant expression effect of differentially expressed genes in Va116×GDH94 and its parents plays a crucial role in the formation of tobacco biomass heterosis. The overdominant expression of genes related to photosynthesis and respiration enhances the photosynthetic ability of Va116×GDH94, reduces respiratory consumption, promotes the increase of biomass, and exhibits obvious heterosis.

Multiple trait comparison and global intestine transcriptional provide new insights into bases of heterosis in hybrid tilapia (Oreochromis niloticus × Oreochromis aureus)

Heterosis has been utilized in aquaculture for many years, yet its molecular basis remains elusive. Therefore, a comprehensive analysis of heterosis was conducted by comparing growth, digestion and biochemistry indices, as well as the intestinal gene expression profiles of Nile tilapia, blue tilapia and their hybrids. The results revealed that hybrid tilapia demonstrated an enhanced growth traits and elevated digestive enzyme activity compared to Nile and blue tilapia. Additionally, the hybrid tilapia displayed superior antioxidants and non-specific immune levels, with increased levels of catalase (CAT), alkaline phosphatase (AKP), acid phosphatase (ACP), glutathione (GSH), superoxide dismutase (SOD), total antioxidant capacity (TAOC), lysozyme, and immunoglobulin M (IgM) relative to Nile and blue tilapia. Moreover, 3392, 2470 and 1261 differentially expressed genes (DEGs) were identified in the intestinal tissues when comparing Nile tilapia to blue tilapia, hybrid tilapia to blue tilapia, and hybrid tilapia to Nile tilapia. Upon classifying the differentially expressed genes (DEGs), non-additively expressed DEGs accounted for 68.1 % of the total DEGs, with dominant and over-dominant expressed DEGs comprising 63.7 % and 4.4 % in the intestines, respectively. These non-additively expressed DEGs were primarily associated with metabolic, digestive, growth, and developmental pathways. This enrichment enhances our comprehension of the molecular underpinnings of growth heterosis in aquatic species.

Toward understanding and utilizing crop heterosis in the age of biotechnology

Heterosis, a universal phenomenon in nature, mainly reflected in the superior productivity, quality, and fitness of F1 hybrids compared with their inbred parents, has been exploited in agriculture and greatly benefited human society in terms of food security. However, the flexible and efficient utilization of heterosis has remained a challenge in hybrid breeding systems because of the limitations of "three-line" and "two-line" methods. In the past two decades, rapidly developed biotechnologies have provided unprecedented conveniences for both understanding and utilizing heterosis. Notably, "third-generation" (3G) hybrid breeding technology together with high-throughput sequencing and gene editing greatly promoted the efficiency of hybrid breeding. Here, we review emerging ideas about the genetic or molecular mechanisms of heterosis and the development of 3G hybrid breeding system in the age of biotechnology. In addition, we summarized opportunities and challenges for optimal heterosis utilization in the future.

Molecular dissection of heterosis in cereal roots and their rhizosphere

Heterosis is already manifested early in root development. Consistent with the dominance model of heterosis, gene expression complementation is a general mechanism that contributes to phenotypic heterosis in maize hybrids. Highly heterozygous F1-hybrids outperform their parental inbred lines, a phenomenon known as heterosis. Utilization of heterosis is of paramount agricultural importance and has been widely applied to increase yield in many crop cultivars. Plant roots display heterosis for many traits and are an important target for further crop improvement. To explain the molecular basis of heterosis, several genetic hypotheses have been proposed. In recent years, high-throughput gene expression profiling techniques have been applied to investigate hybrid vigor. Consistent with the classical genetic dominance model, gene expression complementation has been demonstrated to be a general mechanism to contribute to phenotypic heterosis in diverse maize hybrids. Functional classification of these genes supported the notion that gene expression complementation can dynamically promote hybrid vigor under fluctuating environmental conditions. Hybrids tend to respond differently to available nutrients in the soil. It was hypothesized that hybrid vigor is promoted through a higher nutrient use efficiency which is linked to an improved root system performance of hybrids in comparison to their inbred parents. Recently, the interaction between soil microbes and their plant host was added as further dimension to disentangle heterosis in the belowground part of plants. Soil microbes influenced the performance of maize hybrids as illustrated in comparisons of sterile soil and soil inhabited by beneficial microorganisms.

SiMYBS3, Encoding a Setaria italica Heterosis-Related MYB Transcription Factor, Confers Drought Tolerance in Arabidopsis

Drought is a major limiting factor affecting grain production. Drought-tolerant crop varieties are required to ensure future grain production. Here, 5597 DEGs were identified using transcriptome data before and after drought stress in foxtail millet (Setaria italica) hybrid Zhangza 19 and its parents. A total of 607 drought-tolerant genes were screened through WGCNA, and 286 heterotic genes were screened according to the expression level. Among them, 18 genes overlapped. One gene, Seita.9G321800, encoded MYBS3 transcription factor and showed upregulated expression after drought stress. It is highly homologous with MYBS3 in maize, rice, and sorghum and was named SiMYBS3. Subcellular localization analysis showed that the SiMYBS3 protein was located in the nucleus and cytoplasm, and transactivation assay showed SiMYBS3 had transcriptional activation activity in yeast cells. Overexpression of SiMYBS3 in Arabidopsis thaliana conferred drought tolerance, insensitivity to ABA, and earlier flowering. Our results demonstrate that SiMYBS3 is a drought-related heterotic gene and it can be used for enhancing drought resistance in agricultural crop breeding.

Comparative transcriptomic reveals the molecular mechanism of maize hybrid Zhengdan538 in response to water deficit

Heterosis, known as one of the most successful strategies for increasing grain yield and abiotic/biotic stress tolerance, has been widely exploited in maize breeding. However, the underlying molecular processes are still to be elucidated. The maize hybrid "Zhengdan538" shows high tolerance to drought stress. The transcriptomes of the seedling leaves of its parents, "ZhengA88" and "ZhengT22" and their reciprocal F₁ hybrid under well-watered and water deficit conditions, were analyzed by RNA sequencing (RNA-Seq). Transcriptome profiling of the reciprocal hybrid revealed 2994-4692 differentially expressed genes (DEGs) under well-watered and water-deficit conditions, which were identified by comparing with their parents. The reciprocal hybrid was more closely related to the parental line "ZhengT22" than to the parental line "ZhengA88" in terms of gene expression patterns under water-deficit condition. Furthermore, genes showed expression level dominance (ELD), especially the high-parental ELD (Class 3 and 5), accounted for the largest proportion of DEGs between the reciprocal F₁ hybrid and their parental lines under water deficit. These ELD genes mainly participated in photosynthesis, energy biosynthesis, and metabolism processes. The results indicated that ELD genes played important roles in hybrid tolerance to water deficit. Moreover, a set of important drought-responsive transcription factors were found to be encoded by the identified ELD genes and are thought to function in improving drought tolerance in maize hybrid plants. Our results provide a better understanding of the molecular mechanism of drought tolerance in hybrid maize.

Single-parent expression complementation contributes to phenotypic heterosis in maize hybrids

The dominance model of heterosis explains the superior performance of F1-hybrids via the complementation of deleterious alleles by beneficial alleles in many genes. Genes active in one parent but inactive in the second lead to single-parent expression (SPE) complementation in maize (Zea mays L.) hybrids. In this study, SPE complementation resulted in approximately 700 additionally active genes in different tissues of genetically diverse maize hybrids on average. We established that the number of SPE genes is significantly associated with mid-parent heterosis (MPH) for all surveyed phenotypic traits. In addition, we highlighted that maternally (SPE_B) and paternally (SPE_X) active SPE genes enriched in gene co-expression modules are highly correlated within each SPE type but separated between these two SPE types. While SPE_B-enriched co-expression modules are positively correlated with phenotypic traits, SPE_X-enriched modules displayed a negative correlation. Gene ontology term enrichment analyses indicated that SPE_B patterns are associated with growth and development, whereas SPE_X patterns are enriched in defense and stress response. In summary, these results link the degree of phenotypic MPH to the prevalence of gene expression complementation observed by SPE, supporting the notion that hybrids benefit from SPE complementation via its role in coordinating maize development in fluctuating environments.

Expression complementation of gene presence/absence polymorphisms in hybrids contributes importantly to heterosis in sunflower

INTRODUCTION

Numerous crops have transitioned to hybrid seed production to increase yields and yield stability through heterosis. However, the molecular mechanisms underlying heterosis and its stability across environments are not yet fully understood.

OBJECTIVES

This study aimed to (1) elucidate the genetic and molecular mechanisms underlying heterosis in sunflower, and (2) determine how heterosis is maintained under different environments.

METHODS

Genome-wide association (GWA) analyses were employed to assess the effects of presence/absence variants (PAVs) and stop codons on 16 traits phenotyped in the sunflower association mapping population at three locations. To link the GWA results to transcriptomic variation, we sequenced the transcriptomes of two sunflower cultivars and their F₁ hybrid (INEDI) under both control and drought conditions and analyzed patterns of gene expression and alternative splicing.

RESULTS

Thousands of PAVs were found to affect phenotypic variation using a relaxed significance threshold, and at most such loci the "absence" allele reduced values of heterotic traits, but not those of non-heterotic traits. This pattern was strengthened for PAVs that showed expression complementation in INEDI. Stop codons were much rarer than PAVs and less likely to reduce heterotic trait values. Hybrid expression patterns were enriched for the GO category, sensitivity to stimulus, but all genotypes responded to drought similarily - by up-regulating water stress response pathways and down-regulating metabolic pathways. Changes in alternative splicing were strongly negatively correlated with expression variation, implying that alternative splicing in this system largely acts to reinforce expression responses.

CONCLUSION

Our results imply that complementation of expression of PAVs in hybrids is a major contributor to heterosis in sunflower, consistent with the dominance model of heterosis. This mechanism can account for yield stability across different environments. Moreover, given the much larger numbers of PAVs in plant vs. animal genomes, it also offers an explanation for the stronger heterotic responses seen in the former.

Transcriptome shock in interspecific F1 allotriploid hybrids between Brassica species

Interspecific hybridization drives the evolution of angiosperms and can be used to introduce novel alleles for important traits or to activate heterosis in crop breeding. Hybridization brings together gene expression networks from two different species, potentially causing global alterations of gene expression in the F1 plants which is called 'transcriptome shock'. Here, we explored such a transcriptome shock in allotriploid Brassica hybrids. We generated interspecific F1 allotriploid hybrids between the allotetraploid species Brassica napus and three accessions of the diploid species Brassica rapa. RNA-seq of the F1 hybrids and the parental plants revealed that 26.34-30.89% of genes were differentially expressed between the parents. We also analyzed expression level dominance and homoeolog expression bias between the parents and the F1 hybrids. The expression-level dominance biases of the Ar, An, and Cn subgenomes was genotype and stage dependent, whereas significant homoeolog expression bias was observed among three subgenomes from different parents. Furthermore, more genes were involved in trans regulation than in cis regulation in allotriploid F1 hybrids. Our findings provide new insights into the transcriptomic responses of cross-species hybrids and hybrids showing heterosis, as well as a new method for promoting the breeding of desirable traits in polyploid Brassica species.

From hybrid genomes to heterotic trait output: Challenges and opportunities

Heterosis (or hybrid vigor) has been widely used in crop seed breeding to improve many key economic traits. Nevertheless, the genetic and molecular basis of this important phenomenon has long remained elusive, constraining its flexible and effective exploitation. Advanced genomic approaches are efficient in characterizing the mechanism of heterosis. Here, we review how the omics approaches, including genomic, transcriptomic, and population genetics methods such as genome-wide association studies, can reveal how hybrid genomes outperform parental genomes in plants. This information opens up opportunities for genomic exploration and manipulation of heterosis in crop breeding.

Integration of mRNA and miRNA Profiling Reveals Heterosis in Oreochromis niloticus × O. aureus Hybrid Tilapia

Heterosis is a widespread biological phenomenon in fishes, in which hybrids have superior traits to parents. However, the underlying molecular basis for heterosis remains uncertain. Heterosis in growth and survival rates is apparent in hybrid tilapia (Oreochromis niloticus ♀ × O. aureus ♂). Comparisons of growth and hematological biochemical characteristics and mRNA and miRNA transcriptional analyses were performed in hybrid and parents tilapia stocks to investigate the underlying molecular basis for heterosis. Growth characteristics and hematological glucose and cholesterol parameters were significantly improved in hybrids. Of 3097 differentially expressed genes (DEGs) and 120 differentially expressed miRNAs (DEMs) identified among three stocks (O. niloticus, O. aureus, and hybrids), 1598 DEGs and 62 DEMs were non-additively expressed in hybrids. Both expression level dominance and overdominance patterns occurred for DEGs and DEMs, indicating that dominance and overdominance models are widespread in the transcriptional and post-transcriptional regulation of genes involved in growth, metabolism, immunity, and antioxidant capacity in hybrid tilapia. Moreover, potential negative regulation networks between DEMs and predicted target DEGs revealed that most DEGs from miRNA-mRNA pairs are up-regulated. Dominance and overdominance models in levels of transcriptome and miRNAome facilitate the integration of advantageous parental alleles into hybrids, contributing to heterosis of growth and improved survival. The present study provides new insights into molecular heterosis in hybrid tilapia, advancing our understanding of the complex mechanisms involved in this phenomenon in aquatic animals.

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.

Single-Parent Expression of Anti-sense RNA Contributes to Transcriptome Complementation in Maize Hybrid

Anti-sense transcription is increasingly being recognized as an important regulator of gene expression. But the transcriptome complementation of anti-sense RNA in hybrid relative to their inbred parents was largely unknown. In this study, we profiled strand-specific RNA sequencing (RNA-seq) in a maize hybrid and its inbred parents (B73 and Mo17) in two tissues. More anti-sense transcripts were present in the hybrid compared with the parental lines. We detected 293 and 242 single-parent expression of anti-sense (SPEA) transcripts in maize immature ear and leaf tissues, respectively. There was little overlap of the SPEA transcripts between the two maize tissues. These results suggested that SPEA is a general mechanism that drives extensive complementation in maize hybrids. More importantly, extremely high-level expression of anti-sense transcripts was associated with low-level expression of the cognate sense transcript by reducing the level of histone H3 lysine 36 methylation (H3K36me3). In summary, these SPEA transcripts increased our knowledge about the transcriptomic complementation in hybrid.

Meta Gene Regulatory Networks in Maize Highlight Functionally Relevant Regulatory Interactions

The regulation of gene expression is central to many biological processes. Gene regulatory networks (GRNs) link transcription factors (TFs) to their target genes and represent maps of potential transcriptional regulation. Here, we analyzed a large number of publically available maize (Zea mays) transcriptome data sets including >6000 RNA sequencing samples to generate 45 coexpression-based GRNs that represent potential regulatory relationships between TFs and other genes in different populations of samples (cross-tissue, cross-genotype, and tissue-and-genotype samples). While these networks are all enriched for biologically relevant interactions, different networks capture distinct TF-target associations and biological processes. By examining the power of our coexpression-based GRNs to accurately predict covarying TF-target relationships in natural variation data sets, we found that presence/absence changes rather than quantitative changes in TF gene expression are more likely associated with changes in target gene expression. Integrating information from our TF-target predictions and previous expression quantitative trait loci (eQTL) mapping results provided support for 68 TFs underlying 74 previously identified trans-eQTL hotspots spanning a variety of metabolic pathways. This study highlights the utility of developing multiple GRNs within a species to detect putative regulators of important plant pathways and provides potential targets for breeding or biotechnological applications.

Dominance Effects and Functional Enrichments Improve Prediction of Agronomic Traits in Hybrid Maize

Single-cross hybrids have been critical to the improvement of maize (Zea mays L.), but the characterization of their genetic architectures remains challenging. Previous studies of hybrid maize have shown the contribution of within-locus complementation effects (dominance) and their differential importance across functional classes of loci. However, they have generally considered panels of limited genetic diversity, and have shown little benefit from genomic prediction based on dominance or functional enrichments. This study investigates the relevance of dominance and functional classes of variants in genomic models for agronomic traits in diverse populations of hybrid maize. We based our analyses on a diverse panel of inbred lines crossed with two testers representative of the major heterotic groups in the U.S. (1106 hybrids), as well as a collection of 24 biparental populations crossed with a single tester (1640 hybrids). We investigated three agronomic traits: days to silking (DTS), plant height (PH), and grain yield (GY). Our results point to the presence of dominance for all traits, but also among-locus complementation (epistasis) for DTS and genotype-by-environment interactions for GY. Consistently, dominance improved genomic prediction for PH only. In addition, we assessed enrichment of genetic effects in classes defined by genic regions (gene annotation), structural features (recombination rate and chromatin openness), and evolutionary features (minor allele frequency and evolutionary constraint). We found support for enrichment in genic regions and subsequent improvement of genomic prediction for all traits. Our results suggest that dominance and gene annotations improve genomic prediction across diverse populations in hybrid maize.

Shared Genetic Control of Root System Architecture between Zea mays and Sorghum bicolor

Determining the genetic control of root system architecture (RSA) in plants via large-scale genome-wide association study (GWAS) requires high-throughput pipelines for root phenotyping. We developed Core Root Excavation using Compressed-air (CREAMD), a high-throughput pipeline for the cleaning of field-grown roots, and Core Root Feature Extraction (COFE), a semiautomated pipeline for the extraction of RSA traits from images. CREAMD-COFE was applied to diversity panels of maize (Zea mays) and sorghum (Sorghum bicolor), which consisted of 369 and 294 genotypes, respectively. Six RSA-traits were extracted from images collected from >3,300 maize roots and >1,470 sorghum roots. Single nucleotide polymorphism (SNP)-based GWAS identified 87 TAS (trait-associated SNPs) in maize, representing 77 genes and 115 TAS in sorghum. An additional 62 RSA-associated maize genes were identified via expression read depth GWAS. Among the 139 maize RSA-associated genes (or their homologs), 22 (16%) are known to affect RSA in maize or other species. In addition, 26 RSA-associated genes are coregulated with genes previously shown to affect RSA and 51 (37% of RSA-associated genes) are themselves transe-quantitative trait locus for another RSA-associated gene. Finally, the finding that RSA-associated genes from maize and sorghum included seven pairs of syntenic genes demonstrates the conservation of regulation of morphology across taxa.

Dynamic patterns of circular and linear RNAs in maize hybrid and parental lines

Hybrid vigor, also known as heterosis, has been widely utilized in agronomic production of maize (Zea mays L.) and other crops. However, the molecular mechanisms underlying heterosis are still not fully understood. To provide a more complete understanding of the transcriptomic dynamics associated with heterosis, we collected a comprehensive set of sequence data on linear mRNA transcripts and circular RNAs (circRNAs) from seedling leaves of two widely used maize inbred lines and their F₁ hybrid at the V4 growth stage. We detected over 25,000 expressed genes with more than 1200 circRNAs that showed dramatic and distinct variations in expression level across the three genotypes. Although most linear and circular transcripts exhibited additive expression in the hybrid, the expression of circRNAs was more likely to be nonadditive. Interestingly, the levels of linear transcripts and their corresponding circRNAs from the same loci showed a significant relationship and coordinated expression mode across all three genotypes. Notably, in the hybrid, allele-specific expression of linear transcripts was significantly associated with the expression of circRNAs from the same locus, suggesting potential regulatory cross talk between linear and circular transcripts. Our study provides a deeper understanding of dynamic variations for both the linear and circular transcriptome in a classical hybrid triplet of maize.

Robust non-syntenic gene expression patterns in diverse maize hybrids during root development

Distantly related maize (Zea mays L.) inbred lines exhibit an exceptional degree of structural genomic diversity, which is probably unique among plants. This study systematically investigated the developmental and genotype-dependent regulation of the primary root transcriptomes of a genetically diverse panel of maize F1-hybrids and their parental inbred lines. While we observed substantial transcriptomic changes during primary root development, we demonstrated that hybrid-associated gene expression patterns, including differential, non-additive, and allele-specific transcriptome profiles, are particularly robust to these developmental fluctuations. For instance, differentially expressed genes with preferential expression in hybrids were highly conserved during development in comparison to their parental counterparts. Similarly, in hybrids a major proportion of non-additively expressed genes with expression levels between the parental values were particularly conserved during development. Importantly, in these expression patterns non-syntenic genes that evolved after the separation of the maize and sorghum lineages were systemically enriched. Furthermore, non-syntenic genes were substantially linked to the conservation of all surveyed gene expression patterns during primary root development. Among all F1-hybrids, between ~40% of the non-syntenic genes with unexpected allelic expression ratios and ~60% of the non-syntenic differentially and non-additively expressed genes were conserved and therefore robust to developmental changes. Hence, the enrichment of non-syntenic genes during primary root development might be involved in the developmental adaptation of maize roots and thus the superior performance of hybrids.

A systems biology approach uncovers a gene co-expression network associated with cell wall degradability in maize

Understanding the mechanisms triggering variation of cell wall degradability is a prerequisite to improving the energy value of lignocellulosic biomass for animal feed or biorefinery. Here, we implemented a multiscale systems approach to shed light on the genetic basis of cell wall degradability in maize. We demonstrated that allele replacement in two pairs of near-isogenic lines at a region encompassing a major quantitative trait locus (QTL) for cell wall degradability led to phenotypic variation of a similar magnitude and sign to that expected from a QTL analysis of cell wall degradability in the F271 × F288 recombinant inbred line progeny. Using DNA sequences within the QTL interval of both F271 and F288 inbred lines and Illumina RNA sequencing datasets from internodes of the selected near-isogenic lines, we annotated the genes present in the QTL interval and provided evidence that allelic variation at the introgressed QTL region gives rise to coordinated changes in gene expression. The identification of a gene co-expression network associated with cell wall-related trait variation revealed that the favorable F288 alleles exploit biological processes related to oxidation-reduction, regulation of hydrogen peroxide metabolism, protein folding and hormone responses. Nested in modules of co-expressed genes, potential new cell-wall regulators were identified, including two transcription factors of the group VII ethylene response factor family, that could be exploited to fine-tune cell wall degradability. Overall, these findings provide new insights into the regulatory mechanisms by which a major locus influences cell wall degradability, paving the way for its map-based cloning in maize.

Comparative transcriptome analysis reveals important roles of nonadditive genes in maize hybrid An'nong 591 under heat stress

BACKGROUND

Heterosis is the superior performance of F₁ hybrids relative to their parental lines for a wide range of traits. In this study, expression profiling and heterosis associated genes were analyzed by RNA sequencing (RNA-Seq) in seedlings of the maize hybrid An'nong 591 and its parental lines under control and heat stress conditions.

RESULTS

Through performing nine pairwise comparisons, the maximum number of differentially expressed genes (DEGs) was detected between the two parental lines, and the minimum number was identified between the F₁ hybrid and the paternal lines under both conditions, which suggested greater genetic contribution of the paternal line to heat stress tolerance. Gene Ontology (GO) enrichment analysis of the 4518 common DEGs indicated that GO terms associated with diverse stress responses and photosynthesis were highly overrepresented in the 76 significant terms of the biological process category. A total of 3970 and 7653 genes exhibited nonadditive expression under control and heat stress, respectively. Among these genes, 2253 (56.8%) genes overlapped, suggesting that nonadditive genes tend to be conserved in expression. In addition, 5400 nonadditive genes were found to be specific for heat stress condition, and further GO analysis indicated that terms associated with stress responses were significantly overrepresented, and 60 genes were assigned to the GO term response to heat. Pathway enrichment analysis indicated that 113 genes were involved in spliceosome metabolic pathways. Nineteen of the 33 overlapping genes assigned to the GO term response to heat showed significantly higher number of alternative splicing (AS) events under heat stress than under control conditions, suggesting that AS of these genes play an important role in response to heat stress.

CONCLUSIONS

This study reveals the transcriptomic divergence of the maize F₁ hybrid and its parental lines under control and heat stress conditions, and provides insight into the underlying molecular mechanisms of heterosis and the response to heat stress in maize.

Dynamic Patterns of Gene Expression Additivity and Regulatory Variation throughout Maize Development

Gene expression variation is a key component underlying phenotypic variation and heterosis. Transcriptome profiling was performed on 23 different tissues or developmental stages of two maize inbreds, B73 and Mo17, as well as their F₁ hybrid. The obtained large-scale datasets provided opportunities to monitor the developmental dynamics of differential expression, additivity for gene expression, and regulatory variation. The transcriptome can be divided into ∼30 000 genes that are expressed in at least one tissue of one inbred and an additional ∼10 000 ″silent" genes that are not expressed in any tissue of any genotype, 90% of which are non-syntenic relative to other grasses. Many (∼74%) of the expressed genes exhibit differential expression in at least one tissue. However, the majority of genes with differential expression do not exhibit consistent differential expression in different tissues. These genes often exhibit tissue-specific differential expression with equivalent expression in other tissues, and in many cases they switch the directionality of differential expression in different tissues. This suggests widespread variation for tissue-specific regulation of gene expression between the two maize inbreds B73 and Mo17. Nearly 5000 genes are expressed in only one parent in at least one tissue (single parent expression) and 97% of these genes are expressed at mid-parent levels or higher in the hybrid, providing extensive opportunities for hybrid complementation in heterosis. In general, additive expression patterns are much more common than non-additive patterns, and this trend is more pronounced for genes with strong differential expression or single parent expression. There is relatively little evidence for non-additive expression patterns that are maintained in multiple tissues. The analysis of allele-specific expression allowed classification of cis- and trans-regulatory variation. Genes with cis-regulatory variation often exhibit additive expression and tend to have more consistent regulatory variation throughout development. In contrast, genes with trans-regulatory variation are enriched for non-additive patterns and often show tissue-specific differential expression. Taken together, this study provides a deeper understanding of regulatory variation and the degree of additive gene expression throughout maize development. The dynamic nature of differential expression, additivity, and regulatory variation imply abundant variability for tissue-specific regulatory mechanisms and suggest that connections between transcriptome and phenome will require expression data from multiple tissues.

Single-Parent Expression Is a General Mechanism Driving Extensive Complementation of Non-syntenic Genes in Maize Hybrids

Maize (Zea mays L.) displays an exceptional degree of structural genomic diversity [1, 2]. In addition, variation in gene expression further contributes to the extraordinary phenotypic diversity and plasticity of maize. This study provides a systematic investigation on how distantly related homozygous maize inbred lines affect the transcriptomic plasticity of their highly heterozygous F₁ hybrids. The classical dominance model of heterosis explains the superiority of hybrid plants by the complementation of deleterious parental alleles by superior alleles of the second parent at many loci [3]. Genes active in one inbred line but inactive in another represent an extreme instance of allelic diversity defined as single-parent expression [4]. We observed on average ∼1,000 such genes in all inbred line combinations during primary root development. These genes consistently displayed expression complementation (i.e., activity) in their hybrid progeny. Consequently, extreme expression complementation is a general mechanism that results on average in ∼600 additionally active genes and their encoded biological functions in hybrids. The modern maize genome is complemented by a set of non-syntenic genes, which emerged after the separation of the maize and sorghum lineages and lack syntenic orthologs in any other grass species [5]. We demonstrated that non-syntenic genes are the driving force of gene expression complementation in hybrids. Among those, the highly diversified families of bZIP and bHLH transcription factors [6] are systematically overrepresented. In summary, extreme gene expression complementation extensively shapes the transcriptomic plasticity of maize hybrids and might therefore be one factor controlling the developmental plasticity of hybrids.

Comparison of chemical, color stability, and phenolic composition from pericarp of nine colored corn unique varieties in a beverage model

The objective was to compare the chemical stability and color of nine unique anthocyanin-rich colored corn varieties named/coded as V1, V2, V3… V9. Extracts were added to a beverage model and stored at 4 °C, 22°C, or 32°C for 12weeks. After 12 weeks of storage at 32°C, variety V6 [high condensed form (CF), high cyanidin-3-O-glucoside (C3G)] had the longest anthocyanin half-life, based on the quantification by HPLC. V3 [high pelargonidin (Pg), high acylated form (C3-mal)] and V5 (high CF, high C3G, high C3-mal) had the most favorable hue. V5 and V6 had some of the smallest changes in color over time. These findings suggest that an abundance of condensed forms with C3G in corn extracts could contribute to the improved stability. Beverage storage parameters also influenced color parameters; low temperatures and low pH enhanced color and anthocyanin stability. The most promising corn varieties for future experiments are V3, V5, and V6 based on color retention.