Comparative expression profiling in meristems of inbred-hybrid triplets of maize based on morphological investigations of heterosis for plant height

Genome-wide chromatin accessibility reveals transcriptional regulation of heterosis in inter-subspecific hybrid rice

The utilization of rice heterosis is essential for ensuring global food security; however, its molecular mechanism remains unclear. In this study, comprehensive analyses of accessible chromatin regions (ACRs), DNA methylation, and gene expression in inter-subspecific hybrid and its parents were performed to determine the potential role of chromatin accessibility in rice heterosis. The hybrid exhibited abundant ACRs, in which the gene ACRs and proximal ACRs were directly related to transcriptional activation rather than the distal ACRs. Regarding the dynamic accessibility contribution of the parents, paternal ZHF1015 transmitted a greater number of ACRs to the hybrid. Accessible genotype-specific target genes were enriched with overrepresented transcription factors, indicating a unique regulatory network of genes in the hybrid. Compared with its parents, the differentially accessible chromatin regions with upregulated chromatin accessibility were much greater than those with downregulated chromatin accessibility, reflecting a stronger regulation in the hybrid. Furthermore, DNA methylation levels were negatively correlated with ACR intensity, and genes were strongly affected by CHH methylation in the hybrid. Chromatin accessibility positively regulated the overall expression level of each genotype. ACR-related genes with maternal Z04A-bias allele-specific expression tended to be enriched during carotenoid biosynthesis, whereas paternal ZHF1015-bias genes were more active in carbohydrate metabolism. Our findings provide a new perspective on the mechanism of heterosis based on chromatin accessibility in inter-subspecific hybrid rice.

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.

Comparative transcriptomic analysis of maize ear heterosis during the inflorescence meristem differentiation stage

BACKGROUND

Heterosis is widely used in many crops and is important for global food safety, and maize is one of the most successful crops to take advantage of heterosis. Gene expression patterns control the development of the maize ear, but the mechanisms by which heterosis affects transcriptional-level control are not fully understood.

RESULTS

In this study, we sampled ear inflorescence meristems (IMs) from the single-segment substitution maize (Zea mays) line lx9801^hlEW2b, which contains the heterotic locus hlEW2b associated with ear width, as well as the receptor parent lx9801, the test parent Zheng58, and their corresponding hybrids Zheng58 × lx9801^hlEW2b (HY) and Zheng58 × lx9801 (CK). After RNA sequencing and transcriptomic analysis, 2531 unique differentially expressed genes (DEGs) were identified between the two hybrids (HY vs. CK). Our results showed that approximately 64% and 48% of DEGs exhibited additive expression in HY and CK, whereas the other genes displayed a non-additive expression pattern. The DEGs were significantly enriched in GO functional categories of multiple metabolic processes, plant organ morphogenesis, and hormone regulation. These essential processes are potentially associated with heterosis performance during the maize ear developmental stage. In particular, 125 and 100 DEGs from hybrids with allele-specific expression (ASE) were specifically identified in HY and CK, respectively. Comparison between the two hybrids suggested that ASE genes were involved in different development-related processes that may lead to the hybrid vigor phenotype during maize ear development. In addition, several critical genes involved in auxin metabolism and IM development were differentially expressed between the hybrids and showed various expression patterns (additive, non-additive, and ASE). Changes in the expression levels of these genes may lead to differences in auxin homeostasis in the IM, affecting the transcription of core genes such as WUS that control IM development.

CONCLUSIONS

Our research suggests that additive, non-additive, and allele-specific expression patterns may fine-tune the expression of crucial DEGs that modulate carbohydrate and protein metabolic processes, nitrogen assimilation, and auxin metabolism to optimal levels, and these transcriptional changes may play important roles in maize ear heterosis. The results provide new information that increases our understanding of the relationship between transcriptional variation and heterosis during maize ear development, which may be helpful for clarifying the genetic and molecular mechanisms of heterosis.

Heterosis of fitness and phenotypic variance in the evolution of a diploid gene regulatory network

Heterosis describes the phenomenon, whereby a hybrid population has higher fitness than an inbred population, which has previously been explained by either Mendelian dominance or overdominance under the general assumption of a simple genotype-phenotype relationship. However, recent studies have demonstrated that genes interact through a complex gene regulatory network (GRN). Furthermore, phenotypic variance is reportedly lower for heterozygotes, and the origin of such variance-related heterosis remains elusive. Therefore, a theoretical analysis linking heterosis to GRN evolution and stochastic gene expression dynamics is required. Here, we investigated heterosis related to fitness and phenotypic variance in a system with interacting genes by numerically evolving diploid GRNs. According to the results, the heterozygote population exhibited higher fitness than the homozygote population, indicating fitness-related heterosis resulting from evolution. In addition, the heterozygote population exhibited lower noise-related phenotypic variance in expression levels than the homozygous population, implying that the heterozygote population is more robust to noise. Furthermore, the distribution of the ratio of heterozygote phenotypic variance to homozygote phenotypic variance exhibited quantitative similarity with previous experimental results. By applying dominance and differential gene expression rather than only a single gene expression model, we confirmed the correlation between heterosis and differential gene expression. We explain our results by proposing that the convex high-fitness region is evolutionarily shaped in the genetic space to gain noise robustness under genetic mixing through sexual reproduction. These results provide new insights into the effects of GRNs on variance-related heterosis and differential gene expression.

Transcriptome and DNA methylome analyses provide insight into the heterosis in flag leaf of inter-subspecific hybrid rice

Flag leaf heterosis of inter-subspecific hybrid rice is suggested to be related to leaf area, gene expression pattern and allele-specific expression, putatively related to DNA methylation differences between the hybrid and its parents. Inter-subspecific hybrid rice combinations of indica × japonica have great potential to broaden genetic diversity and enhance the heterosis. The genetic and epigenetic molecular mechanism of its heterosis is not completely understood. Here, the dissection of gene expression and epigenetic regulation of an elite inter-subspecific hybrid rice were reported. In the hybrid, plant height, flag leaf area and Pn showed significant heterosis at the heading stage. Chloroplast-related differentially expressed genes (DEGs) and 530 allele-specific expression genes in hybrid were identified. Analysis of the genome-wide distribution of DNA methylation (5-methylcytosine, 5mC) and its association with transcription showed that there were variant DNA methylation maps and that the regulation of gene expression levels was negatively regulated by DNA methylation in the inter-subspecific hybrid rice. Differentially methylated DEGs were significantly enriched in photosynthetic functions. Moreover, distinct 5mC sequence contexts and distinct functional elements (promoter/gene body) may have different influences on heterosis related genes. The data identified heterosis related molecular mechanisms in inter-subspecific hybrid rice and suggested that epigenetic changes could extensively influence the flag leaf gene expression and heterosis.

Transcriptome analysis reveals the molecular mechanisms of heterosis on thermal resistance in hybrid abalone

BACKGROUND

Heterosis has been exploited for decades in different animals and crops due to it resulting in dramatic increases in yield and adaptability. Hybridization is a classical breeding method that can effectively improve the genetic characteristics of organisms through heterosis. Abalone has become an increasingly economically important aquaculture resource with high commercial value. However, due to changing climate, abalone is now facing serious threats of high temperature in summer. Interspecific hybrid abalone (Haliotis gigantea ♀ × H. discus hannai ♂, SD) has been cultured at large scale in southern China and has been shown high survival rates under heat stress in summer. Therefore, SD has become a good model material for heterosis research, but the molecular basis of heterosis remains elusive.

RESULTS

Heterosis in thermal tolerance of SD was verified through Arrhenius break temperatures (ABT) of cardiac performance in this study. Then RNA-Sequencing was conducted to obtain gene expression patterns and alternative splicing events at control temperature (20 °C) and heat stress temperature (30 °C). A total of 356 (317 genes), 476 (435genes), and 876 (726 genes) significantly diverged alternative splicing events were identified in H. discus hannai (DD), H. gigantea (SS), and SD in response to heat stress, respectively. In the heat stress groups, 93.37% (20,512 of 21,969) of the expressed genes showed non-additive expression patterns, and over-dominance expression patterns of genes account for the highest proportion (40.15%). KEGG pathway enrichment analysis showed that the overlapping genes among common DEGs and NAGs were significantly enriched in protein processing in the endoplasmic reticulum, mitophagy, and NF-κB signaling pathway. In addition, we found that among these overlap genes, 39 genes had undergone alternative splicing events in SD. These pathways and genes may play an important role in the thermal resistance of hybrid abalone.

CONCLUSION

More alternative splicing events and non-additive expressed genes were detected in hybrid under heat stress and this may contribute to its thermal heterosis. These results might provide clues as to how hybrid abalone has a better physiological regulation ability than its parents under heat stress, to increase our understanding of heterosis in abalone.

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.

Assessment of hybrid vigour, dominance effect and hybrids regeneration potential in the genus Citrullus

The current research was directed to explore hybrid vigour for hybrids resulting from crosses among six genotypes of the genus Citrullus, including C. mucosospermus. For such purpose, the mid parent heterosis and the best parent heterosis were assessed. Non parametric method related to homogeneity Chi-square at 5 % likelihood was applied to compare the regeneration potential of progenies. Student's parametric test at 5% was used to separe two means. Six parental genotypes and 16 hybrid families were evaluated for heterosis. The results showed a heterosis effect for all the characters studied nevertheless, this one varies according to the crossing. The observed hybrid vigour would be due to a superdominance or partial dominance effect. In addition, progeny from crosses have the same regenerative potential on both study sites.

Transcriptome Analysis by RNA-Seq Reveals Genes Related to Plant Height in Two Sets of Parent-hybrid Combinations in Easter lily (Lilium longiflorum)

In this study, two different hybrids of Easter lily (Lilium longiflorum), obtained from two cross combinations, along with their four parents were sequenced by high–throughput RNA–sequencing (RNA–Seq) to find out differentially expressed gene in parent-hybrid combinations. The leaf mRNA profiles of two hybrids and their four parents were RNA–sequenced with a view to identify the potential candidate genes related to plant height heterosis. In both cross combinations, based to morphological traits mid–parent heterosis (MPH) was higher than high–parent heterosis (HPH) for plant height, leaf length, and number of flowers whereas HPH was higher than MPH for flowering time. A total of 4,327 differentially expressed genes (DEGs) were identified through RNA–Seq between the hybrids and their parents based on fold changes (FC) ≥ 2 for up– and ≤ –2 for down–regulation. Venn diagram analysis revealed that there were 703 common DEGs in two hybrid combinations, those were either up– or down–regulated. Most of the commonly expressed DEGs exhibited higher non–additive effects especially overdominance (75.9%) rather than additive (19.4%) and dominance (4.76%) effects. Among the 384 functionally annotated DEGs identified through Blast2GO tool, 12 DEGs were up–regulated and 16 of them were down–regulated in a similar fashion in both hybrids as revealed by heat map analysis. These 28 universally expressed DEGs were found to encode different types of proteins and enzymes those might regulate heterosis by modulating growth, development and stress–related functions in lily. In addition, gene ontology (GO) analysis of 260 annotated DEGs revealed that biological process might play dominant role in heterotic expression. In this first report of transcriptome sequencing in Easter lily, the notable universally up-regulated DEGs annotated ABC transporter A family member–like, B3 domain–containing, disease resistance RPP13/1, auxin–responsive SAUR68–like, and vicilin–like antimicrobial peptides 2–2 proteins those were perhaps associated with plant height heterosis. The genes expressed universally due to their overdominace function perhaps influenced MPH for greater plant height― largely by modulating biological processes involved therein. The genes identified in this study might be exploited in heterosis breeding for plant height of L. longiflorum.

Comparative transcriptome analysis between inbred and hybrids reveals molecular insights into yield heterosis of upland cotton

BACKGROUND

Utilization of heterosis has greatly improved the productivity of many crops worldwide. Understanding the potential molecular mechanism about how hybridization produces superior yield in upland cotton is critical for efficient breeding programs.

RESULTS

In this study, high, medium, and low hybrids varying in the level of yield heterosis were screened based on field experimentation of different years and locations. Phenotypically, high hybrid produced a mean of 14% more seed cotton yield than its better parent. Whole-genome RNA sequencing of these hybrids and their four inbred parents was performed using different tissues of the squaring stage. Comparative transcriptomic differences in each hybrid parent triad revealed a higher percentage of differentially expressed genes (DEGs) in each tissue. Expression level dominance analysis identified majority of hybrids DEGs were biased towards parent like expressions. An array of DEGs involved in ATP and protein binding, membrane, cell wall, mitochondrion, and protein phosphorylation had more functional annotations in hybrids. Sugar metabolic and plant hormone signal transduction pathways were most enriched in each hybrid. Further, these two pathways had most mapped DEGs on known seed cotton yield QTLs. Integration of transcriptome, QTLs, and gene co-expression network analysis discovered genes Gh_A03G1024, Gh_D08G1440, Gh_A08G2210, Gh_A12G2183, Gh_D07G1312, Gh_D08G1467, Gh_A03G0889, Gh_A08G2199, and Gh_D05G0202 displayed a complex regulatory network of many interconnected genes. qRT-PCR of these DEGs was performed to ensure the accuracy of RNA-Seq data.

CONCLUSIONS

Through genome-wide comparative transcriptome analysis, the current study identified nine key genes and pathways associated with biological process of yield heterosis in upland cotton. Our results and data resources provide novel insights and will be useful for dissecting the molecular mechanism of yield heterosis in cotton.

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.

Gibberellins and Heterosis in Crops and Trees: An Integrative Review and Preliminary Study with Brassica

Heterosis, or hybrid vigor, has contributed substantially to genetic improvements in crops and trees and its physiological basis involves multiple processes. Four associations with the phytohormone gibberellin (GA) indicate its involvement in the regulation of heterosis for shoot growth in maize, sorghum, wheat, rice, tomato and poplar. (1) Inbreds somewhat resemble GA-deficient dwarfs and are often highly responsive to exogenous GA₃. (2) Levels of endogenous GAs, including the bioeffector GA₁, its precursors GA₁₉ and GA₂₀, and/or its metabolite GA₈, are higher in some fast-growing hybrids than parental genotypes. (3) Oxidative metabolism of applied [³H]GAs is more rapid in vigorous hybrids than inbreds, and (4) heterotic hybrids have displayed increased expression of GA biosynthetic genes including GA 20-oxidase and GA 3-oxidase. We further investigated Brassicarapa, an oilseed rape, by comparing two inbreds (AO533 and AO539) and their F₁ hybrid. Seedling emergence was faster in the hybrid and potence ratios indicated dominance for increased leaf number, area and mass, and stem mass. Overdominance (heterosis) was displayed for root mass, leading to slight heterosis for total plant mass. Stem contents of GA_19,20,1 were similar across the Brassica genotypes and increased prior to bolting; elongation was correlated with endogenous GA but heterosis for shoot growth was modest. The collective studies support a physiological role for GAs in the regulation of heterosis for shoot growth in crops and trees, and the Brassica study encourages further investigation of heterosis for root growth.

Transcriptome and Metabolite Profiling Reveal Novel Insights into Volatile Heterosis in the Tea Plant (Camellia Sinensis)

Tea aroma is a key indicator for evaluating tea quality. Although notable success in tea aroma improvement has been achieved with heterosis breeding technology, the molecular basis underlying heterosis remains largely unexplored. Thus, the present report studies the tea plant volatile heterosis using a high-throughput next-generation RNA-seq strategy and gas chromatography-mass spectrometry. Phenotypically, we found higher terpenoid volatile and green leaf volatile contents by gas chromatography-mass spectrometry in the F1 hybrids than in their parental lines. Volatile heterosis was obvious in both F1 hybrids. At the molecular level, the comparative transcriptomics analysis revealed that approximately 41% (9027 of 21,995) of the genes showed non-additive expression, whereas only 7.83% (1723 of 21,995) showed additive expression. Among the non-additive genes, 42.1% showed high parental dominance and 17.6% showed over-dominance. Among different expression genes with high parental dominance and over-dominance expression patterns, KEGG and GO analyses found that plant hormone signal transduction, tea plant physiological process related pathways and most pathways associated with tea tree volatiles were enriched. In addition, we identified multiple genes (CsDXS, CsAATC2, CsSPLA2, etc.) and transcription factors (CsMYB1, CsbHLH79, CsWRKY40, etc.) that played important roles in tea volatile heterosis. Based on transcriptome and metabolite profiling, we conclude that non-additive action plays a major role in tea volatile heterosis. Genes and transcription factors involved in tea volatiles showing over-dominance expression patterns can be considered candidate genes and provide novel clues for breeding high-volatile tea varieties.

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.

Transcriptome analysis of growth heterosis in pearl oyster Pinctada fucata martensii

Heterosis improves growth and survival of shellfish species. Although breeders have widely exploited heterosis, its underlying molecular mechanisms remain unclear. In this study, a 2 × 2 complete diallel cross was facilitated between two full-sib families to produce two inbred families (A and D) and their reciprocal hybrid families (B and C) of pearl oyster Pinctada fucata martensii. Growth traits of the four families were compared at the adult stages. Transcriptome analysis was conducted on the four families using an Illumina sequencing platform. The results revealed that the growth traits of the four families significantly varied (P < 0.05). The mid-parent heterosis values of shell length, shell height, shell width, shell weight, and total weight were 12.9%, 14.9%, 18.2%, 17.2%, and 33.2%, respectively. The B- and C-inbred (A and D) triads had 79 and 68 differentially expressed genes (DEGs), respectively, which were dominantly nonadditive, including overdominance, underdominance, and low-parent dominance. Gene ontology term analysis showed that the DEGs in the B- and C-inbred triads were enriched for metabolic process, cellular process cell part, binding, and catalytic activity. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis indicated that the DEGs in the B- and C-inbred triads were involved in focal adhesion, the P13K-Akt signaling pathway, the mRNA surveillance pathway, and the focal adhesion pathway. The reliability of the sequencing data was confirmed by real-time polymerase chain reaction analysis of six growth-related genes. The findings of this study provide new insights into heterosis for growth traits and the design of genetic breeding programs for this species.

Heterosis-related genes under different planting densities in maize

Heterosis and increasing planting density have contributed to improving maize grain yield (GY) for several decades. As planting densities increase, the GY per plot also increases, whereas the contribution of heterosis to GY decreases. There are trade-offs between heterosis and planting density, and the transcriptional characterization of heterosis may explain the mechanism involved. In this study, 48 transcriptome libraries were sequenced from four inbred Chinese maize lines and their F1 hybrids. They were planted at densities of 45000 and 67500 plants ha-1. Maternal-effect differentially expressed genes (DEGs) played important roles in processes related to photosynthesis and carbohydrate biosynthesis and metabolism. Paternal-effect DEGs participated in abiotic/biotic stress response and plant hormone production under high planting density. Weighted gene co-expression network analysis revealed that high planting density induced heterosis-related genes regulating abiotic/biotic stress response, plant hormone biosynthesis, and ubiquitin-mediated proteolysis, but repressed other genes regulating energy formation. Under high planting density, maternal genes were mainly enriched in the photosynthesis reaction center, while paternal genes were mostly concentrated in the peripheral antenna system. Four important genes were identified in maize heterosis and high planting density, all with functions in photosynthesis, starch biosynthesis, auxin metabolism, gene silencing, and RNAi.

Transcriptome analysis of wheat seedling and spike tissues in the hybrid Jingmai 8 uncovered genes involved in heterosis

Transcriptome analysis was carried out for wheat seedlings and spikes from hybrid Jingmai 8 and both inbred lines to unravel mechanisms underlying heterosis. Heterosis, known as one of the most successful strategies for increasing crop yield, has been widely exploited in plant breeding systems. Despite its great importance, the molecular mechanism underlying heterosis remains elusive. In the present study, RNA sequencing (RNA-seq) was performed on the seedling and spike tissues of the wheat (Triticum aestivum) hybrid Jingmai 8 (JM8) and its homozygous parents to unravel the underlying mechanisms of wheat heterosis. In total, 1686 and 2334 genes were identified as differentially expressed genes (DEGs) between the hybrid and the two inbred lines in seedling and spike tissues, respectively. Gene Ontology analysis revealed that DEGs from seedling tissues were significantly enriched in processes involved in photosynthesis and carbon fixation, and the majority of these DEGs expressed at a higher level in JM8 compared to both inbred lines. In addition, cell wall biogenesis and protein biosynthesis-related pathways were also significantly represented. These results confirmed that a combination of different pathways could contribute to heterosis. The DEGs between the hybrid and the two inbred progenitors from the spike tissues were significantly enriched in biological processes related to transcription, RNA biosynthesis and molecular function categories related to transcription factor activities. Furthermore, transcription factors such as NAC, ERF, and TIF-IIA were highly expressed in the hybrid JM8. These results may provide valuable insights into the molecular mechanisms underlying wheat heterosis.

Genome-wide proteomic profiling reveals the role of dominance protein expression in heterosis in immature maize ears

Heterosis refers to the phenomenon in which hybrid progeny show superior performance relative to their parents. Early maize ear development shows strong heterosis in ear architecture traits and greatly affects grain yield. To explore the underlying molecular mechanisms, genome-wide proteomics of immature ears of maize hybrid ZD909 and its parents were analyzed using tandem mass tag (TMT) technology. A total of 9,713 proteins were identified in all three genotypes. Among them, 3,752 (38.6%) proteins were differentially expressed between ZD909 and its parents. Multiple modes of protein action were discovered in the hybrid, while dominance expression patterns accounted for 63.6% of the total differentially expressed proteins (DEPs). Protein pathway enrichment analysis revealed that high parent dominance proteins mainly participated in carbon metabolism and nitrogen assimilation processes. Our results suggested that the dominant expression of favorable alleles related to C/N metabolism in the hybrid may be essential for ZD909 ear growth and heterosis formation. Integrated analysis of proteomic and quantitative trait locus (QTL) data further support our DEP identification and provide useful information for the discovery of genes associated with ear development. Our study provides comprehensive insight into the molecular mechanisms underlying heterosis in immature maize ears from a proteomic perspective.

Variable Level of Dominance of Candidate Genes Controlling Drought Functional Traits in Maize Hybrids

Breeding maize for drought tolerance necessitates the knowledge on tolerant genotypes, molecular basis of drought tolerance mechanism, action, and expression pattern of genes. Studying the expression pattern and gene action of candidate genes during drought stress in the hybrids will help in choosing target genes for drought tolerance breeding. In the present investigation, a set of five hybrids and their seven parents with a variable level of tolerance to drought stress was selected to study the magnitude and the direction of 52 drought-responsive candidate genes distributed across various biological functions, viz., stomatal regulation, root development, detoxification, hormone signaling, photosynthesis, and sugar metabolism. The tolerant parents, HKI1105 and CML425, and their hybrid, ADWLH2, were physiologically active under drought stress, since vital parameters viz., chlorophyll, root length and relative water content, were on par with the respective well-watered control. All the genes were up-regulated in ADWLH2, many were down-regulated in HM8 and HM9, and most were down-regulated in PMH1 and PMH3 in the shoots and roots. The nature of the gene action was controlled by the parental combination rather than the parent per se. The differentially expressed genes in all five hybrids explained a mostly non-additive gene action over additivity, which was skewed toward any of the parental lines. Tissue-specific gene action was also noticed in many of the genes. The non-additive gene action is driven by genetic diversity, allele polymorphism, events during gene regulation, and small RNAs under the stress condition. Differential regulation and cross-talk of genes controlling various biological functions explained the basis of drought tolerance in subtropical maize hybrids. The nature of the gene action and the direction of the expression play crucial roles in designing introgression and hybrid breeding programmes to breed drought tolerant maize hybrids.

Transcriptome profiling and comparison of maize ear heterosis during the spikelet and floret differentiation stages

BACKGROUND

Hybridization is a prominent process in the evolution of crop plants that can give rise to gene expression variation, phenotypic novelty and heterosis. Maize is the most successful crop in utilizing heterosis. The development of hybrid maize ears exhibits strong heterotic vigor and greatly affects maize yield. However, a comprehensive perspective on transcriptional variation and its correlation with heterosis during maize ear development is not available.

RESULTS

Using RNA sequencing technology, we investigated the transcriptome profiles of maize ears in the spikelet and floret differentiation stages of hybrid ZD808 and its parents CL11 and NG5. Our results revealed that 53.9% (21,258) of maize protein-coding genes were transcribed in at least one genotype. In both development stages, significant numbers of genes were differentially expressed between the hybrid and its parents. Gene expression inheritance analysis revealed approximately 80% of genes were expressed additively, which suggested that the complementary effect may play a foundation role in maize ear heterosis. Among non-additively expressed genes, NG5-dominant genes were predominant. Analyses of the allele-specific gene expression in hybrid identified pervasive allelic imbalance and significant preferential expression of NG5 alleles in both developmental stages. The results implied that NG5 may provide beneficial alleles that contribute greatly to heterosis. Further comparison of parental and hybrid allele-specific expression suggested that gene expression variation is largely attributable to cis-regulatory variation in maize. The cis-regulatory variations tend to preserve the allelic expression levels in hybrid and result in additive expression. Comparison between the two development stages revealed that allele-specific expression and cis-/trans-regulatory variations responded differently to developmental cues, which may lead to stage-specific vigor phenotype during maize ear development.

CONCLUSION

Our research suggests that cis-regulated additive expression may fine-tune gene expression level into an optimal status and play a foundation role in maize ear heterosis. Our work provides a comprehensive insight into transcriptional variation and its correlation with heterosis during maize ear development. The knowledge gained from this study presents novel opportunity to improve our maize varieties.

Transcriptome shock invokes disruption of parental expression-conserved genes in tetraploid wheat

Allopolyploidy often triggers phenotypic novelty and gene expression remolding in the resulting polyploids. In this study, we employed multiple phenotypic and genetic approaches to investigate the nature and consequences of allotetraploidization between A- and S-subgenome of tetraploid wheat. Results showed that karyotype of the nascent allopolyploid plants (AT2) is stable but they showed clear novelty in multiple morphological traits which might have positively contributed to the initial establishment of the tetraploids. Further microarray-based transcriptome profiling and gene-specific cDNA-pyrosequencing have documented that transcriptome shock was exceptionally strong in AT2, but a substantial proportion of the induced expression changes was rapidly stabilized in early generations. Meanwhile, both additive and nonadditive expression genes showed extensive homeolog expression remodeling and which have led to the subgenome expression dominance in leaf and young inflorescence of AT2. Through comparing the homeolog-expressing patterns between synthetic and natural tetraploid wheats, it appears that the shock-induced expression changes at both the total expression level and subgenome homeolog partitioning are evolutionarily persistent. Together, our study shed new light on how gene expression changes have rapidly occurred at the initial stage following allotetraploidization, as well as their evolutionary relevance, which may have implications for wheat improvements.

Comparative transcriptome analysis among parental inbred and crosses reveals the role of dominance gene expression in heterosis in Drosophila melanogaster

We observed heteroses for body weight in Drosophila melanogaster after generating hybrids from three inbred lines. To better understand the mechanism for this phenomenon at the mRNA level, we compared the mRNA profiles of the parental and hybrid lines using high-throughput RNA-seq. A total of 5877 differentially expressed genes (DEGs) were found and about 92% of these exhibited parental expression level dominance. Genes in the dominance category were functionally characterized using the Kyoto Encyclopedia of Genes and Genomes (KEGG) and the gene classifications offered by the Gene Ontology (GO) Consortium. The analysis identified genes associated with crucial processes such as development and growth in all three crosses. Functional assignments involving aminoglycan metabolism, starch and sucrose metabolism, and galactose metabolism are significantly overrepresented amongst the 215 common dominance DEGs. We conclude that dominance DEGs are important in heteroses in Drosophila melanogaster and contribute specifically to body weight heterosis.

Ectopic expression of a maize hybrid up-regulated gene, ErbB-3 binding Protein 1 (ZmEBP1), increases organ size by promoting cell proliferation in Arabidopsis

The alteration of gene expression in hybrids may be an important factor promoting phenotypic variation and plasticity. To provide insight into the underlying molecular basis of maize heterosis in terms of the kernel number per ear, we established DGE profiles for the immature ears of maize hybrid Zong3/87-1 and its parental lines at the floral organ differentiation stage. Among 4,337 identified differentially expressed genes, 4,021 (92%) exhibited nonadditive expression patterns in the hybrid. Notably, the maize homolog of Arabidopsis EBP1, designated ZmEBP1, displayed an overdominant expression pattern in the Zong3/87-1 hybrid. Moreover, the results of qRT-PCR revealed that the ZmEBP1 gene was upregulated in the immature ears of the reciprocal hybrids Zong3/87-1 and 87-1/Zong3 at different developmental stages. Additionally, ectopic expression of ZmEBP1 in Arabidopsis increased organ size, which was mainly attributed to an increase in cell numbers, rather than cell size. Considering all of these findings together, we speculate that upregulation of ZmEBP1 in maize hybrids may accelerate cell proliferation and promote ear development.

Dissecting repulsion linkage in the dwarfing gene Dw3 region for sorghum plant height provides insights into heterosis

Heterosis is a main contributor to yield increase in many crop species. Different mechanisms have been proposed for heterosis: dominance, overdominance, epistasis, epigenetics, and protein metabolite changes. However, only limited examples of molecular dissection and validation of these mechanisms are available. Here, we present an example of discovery and validation of heterosis generated by a combination of repulsion linkage and dominance. Using a recombinant inbred line population, a separate quantitative trait locus (QTL) for plant height (qHT7.1) was identified near the genomic region harboring the known auxin transporter Dw3 gene. With two loci having repulsion linkage between two inbreds, heterosis in the hybrid can appear as a single locus with an overdominance mode of inheritance (i.e., pseudo-overdominance). Individually, alleles conferring taller plant height exhibited complete dominance over alleles conferring shorter height. Detailed analyses of different height components demonstrated that qHT7.1 affects both the upper and lower parts of the plant, whereas Dw3 affects only the part below the flag leaf. Computer simulations show that repulsion linkage could influence QTL detection and estimation of effect in segregating populations. Guided by findings in linkage mapping, a genome-wide association study of plant height with a sorghum diversity panel pinpointed genomic regions underlying the trait variation, including Dw1, Dw2, Dw3, Dw4, and qHT7.1. Multilocus mixed model analysis confirmed the advantage of complex trait dissection using an integrated approach. Besides identifying a specific genetic example of heterosis, our research indicated that integrated molecular dissection of complex traits in different population types can enable plant breeders to fine tune the breeding process for crop production.

Shared Genomic Regions Between Derivatives of a Large Segregating Population of Maize Identified Using Bulked Segregant Analysis Sequencing and Traditional Linkage Analysis

Delayed transition from the vegetative stage to the reproductive stage of development and increased plant height have been shown to increase biomass productivity in grasses. The goal of this project was to detect quantitative trait loci using extremes from a large synthetic population, as well as a related recombinant inbred line mapping population for these two traits. Ten thousand individuals from a B73 × Mo17 noninbred population intermated for 14 generations (IBM Syn14) were grown at a density of approximately 16,500 plants ha(-1). Flowering time and plant height were measured within this population. DNA was pooled from the 46 most extreme individuals from each distributional tail for each of the traits measured and used in bulk segregant analysis (BSA) sequencing. Allelic divergence at each of the ∼1.1 million SNP loci was estimated as the difference in allele frequencies between the selected extremes. Additionally, 224 intermated B73 × Mo17 recombinant inbred lines were concomitantly grown at a similar density adjacent to the large synthetic population and were assessed for flowering time and plant height. Using the BSA sequencing method, 14 and 13 genomic regions were identified for flowering time and plant height, respectively. Linkage mapping with the RIL population identified eight and three regions for flowering time and plant height, respectively. Of the regions identified, three colocalized between the two populations for flowering time and two colocalized for plant height. This study demonstrates the utility of using BSA sequencing for the dissection of complex quantitative traits important for production of lignocellulosic ethanol.

Hierarchical additive effects on heterosis in rice (Oryza sativa L.)

Exploitation of heterosis in crops has contributed greatly to improvement in global food and energy production. In spite of the pervasive importance of heterosis, a complete understanding of its mechanisms has remained elusive. In this study, a small test-crossed rice population was constructed to investigate the formation mechanism of heterosis for 13 traits. The results of the relative mid-parent heterosis and modes of inheritance of all investigated traits demonstrated that additive effects were the foundation of heterosis for complex traits in a hierarchical structure, and multiplicative interactions among the component traits were the framework of heterosis in complex traits. Furthermore, new balances between unit traits and related component traits provided hybrids with the opportunity to achieve an optimal degree of heterosis for complex traits. This study dissected heterosis of both reproductive and vegetative traits from the perspective of hierarchical structure for the first time. Additive multiplicative interactions of component traits were proven to be the origin of heterosis in complex traits. Meanwhile, more attention should be paid to component traits, rather than complex traits, in the process of revealing the mechanism of heterosis.

Investigating the molecular genetic basis of heterosis for internode expansion in maize by microRNA transcriptomic deep sequencing

Heterosis has been used widely in the breeding of maize and other crops and plays an important role in increasing yield, improving quality, and enhancing stress resistance, but its molecular mechanism is far from clear. To determine whether microRNA (miRNA)-dependent gene regulation is responsible for heterosis of elongating internodes below the ear and ear height in maize, a deep-sequencing strategy was applied to the elite hybrid Xundan20, which is currently cultivated widely in China, and its two parents. RNA was extracted from the eighth internode because it shows clear internode length heterosis. A total of 99 conserved maize miRNAs were detected in both the hybrid and parental lines. Most of these miRNAs were expressed nonadditively in the hybrid compared with its parental lines. These results indicated that miRNAs might participate in heterosis during internode expansion in maize and exert an influence on ear and plant height via the repression of their target genes. In total, eight novel miRNAs belonging to four miRNA families were predicted in the expanding internode. Global repression of miRNAs in the hybrid, which might result in enhanced gene expression, might be one reason why the hybrid shows longer internodes and taller seedlings compared with its parental lines.

Heterosis in early maize ear inflorescence development: a genome-wide transcription analysis for two maize inbred lines and their hybrid

Heterosis, or hybrid vigor, contributes to superior agronomic performance of hybrids compared to their inbred parents. Despite its importance, little is known about the genetic and molecular basis of heterosis. Early maize ear inflorescences formation affects grain yield, and are thus an excellent model for molecular mechanisms involved in heterosis. To determine the parental contributions and their regulation during maize ear-development-genesis, we analyzed genome-wide digital gene expression profiles in two maize elite inbred lines (B73 and Mo17) and their F₁ hybrid using deep sequencing technology. Our analysis revealed 17,128 genes expressed in these three genotypes and 22,789 genes expressed collectively in the present study. Approximately 38% of the genes were differentially expressed in early maize ear inflorescences from heterotic cross, including many transcription factor genes and some presence/absence variations (PAVs) genes, and exhibited multiple modes of gene action. These different genes showing differential expression patterns were mainly enriched in five cellular component categories (organelle, cell, cell part, organelle part and macromolecular complex), five molecular function categories (structural molecule activity, binding, transporter activity, nucleic acid binding transcription factor activity and catalytic activity), and eight biological process categories (cellular process, metabolic process, biological regulation, regulation of biological process, establishment of localization, cellular component organization or biogenesis, response to stimulus and localization). Additionally, a significant number of genes were expressed in only one inbred line or absent in both inbred lines. Comparison of the differences of modes of gene action between previous studies and the present study revealed only a small number of different genes had the same modes of gene action in both maize seedlings and ear inflorescences. This might be an indication that in different tissues or developmental stages, different global expression patterns prevail, which might nevertheless be related to heterosis. Our results support the hypotheses that multiple molecular mechanisms (dominance and overdominance modes) contribute to heterosis.

What is crop heterosis: new insights into an old topic

Heterosis (or hybrid vigor) refers to a natural phenomenon whereby hybrid offspring of genetically diverse individuals out-perform their parents in multiple traits including yield, adaptability and resistances to biotic and abiotic stressors. Innovations in technology and research continue to clarify the mechanisms underlying crop heterosis, however the intrinsic relationship between the biological basis of heterosis remain unclear. In this review, we aim to provide insight into the molecular genetic basis of heterosis by presenting recent advances in the 'omics' of heterosis and the role of non-coding regions, particularly in relation to energy-use efficiency. We propose that future research should focus on integrating the expanding datasets from different species and hybrid combinations, to mine key heterotic genes and unravel interactive 'omics' networks associated with heterosis. Improved understanding of heterosis and the biological basis for its manipulation in agriculture should help to streamline its use in enhancing crop productivity.

Comparative transcriptional profiling of three super-hybrid rice combinations

Utilization of heterosis has significantly increased rice yields. However, its mechanism remains unclear. In this study, comparative transcriptional profiles of three super-hybrid rice combinations, LY2163, LY2186 and LYP9, at the flowering and filling stages, were created using rice whole-genome oligonucleotide microarray. The LY2163, LY2186 and LYP9 hybrids yielded 1193, 1630 and 1046 differentially expressed genes (DGs), accounting for 3.2%, 4.4% and 2.8% of the total number of genes (36,926), respectively, after using the z-test (p < 0.01). Functional category analysis showed that the DGs in each hybrid combination were mainly classified into the carbohydrate metabolism and energy metabolism categories. Further analysis of the metabolic pathways showed that DGs were significantly enriched in the carbon fixation pathway (p < 0.01) for all three combinations. Over 80% of the DGs were located in rice quantitative trait loci (QTLs) of the Gramene database, of which more than 90% were located in the yield related QTLs in all three combinations, which suggested that there was a correlation between DGs and rice heterosis. Pathway Studio analysis showed the presence of DGs in the circadian regulatory network of all three hybrid combinations, which suggested that the circadian clock had a role in rice heterosis. Our results provide information that can help to elucidate the molecular mechanism underlying rice heterosis.

Conservation and divergence of transcriptomic and epigenomic variation in maize hybrids

Background

Recent genome-wide studies suggested that in addition to genetic variations, epigenetic variations may also be associated with differential gene expression and growth vigor in plant hybrids. Maize is an ideal model system for the study of epigenetic variations in hybrids given the significant heterotic performance, the well-known complexity of the genome, and the rich history in epigenetic studies. However, integrated comparative transcriptomic and epigenomic analyses in different organs of maize hybrids remain largely unexplored.

Results

Here, we generated integrated maps of transcriptomes and epigenomes of shoots and roots of two maize inbred lines and their reciprocal hybrids, and globally surveyed the epigenetic variations and their relationships with transcriptional divergence between different organs and genotypes. We observed that whereas histone modifications vary both between organs and between genotypes, DNA methylation patterns are more distinguishable between genotypes than between organs. Histone modifications were associated with transcriptomic divergence between organs and between hybrids and parents. Further, we show that genes up-regulated in both shoots and roots of hybrids were significantly enriched in the nucleosome assembly pathway. Interestingly, 22- and 24-nt siRNAs were shown to be derived from distinct transposable elements, and for different transposable elements in both shoots and roots, the differences in siRNA activity between hybrids and patents were primarily driven by different siRNA species.

Conclusions

These results suggest that despite variations in specific genes or genomic loci, similar mechanisms may account for the genome-wide epigenetic regulation of gene activity and transposon stability in different organs of maize hybrids.

Mendelian and non-Mendelian regulation of gene expression in maize

Transcriptome variation plays an important role in affecting the phenotype of an organism. However, an understanding of the underlying mechanisms regulating transcriptome variation in segregating populations is still largely unknown. We sought to assess and map variation in transcript abundance in maize shoot apices in the intermated B73×Mo17 recombinant inbred line population. RNA–based sequencing (RNA–seq) allowed for the detection and quantification of the transcript abundance derived from 28,603 genes. For a majority of these genes, the population mean, coefficient of variation, and segregation patterns could be predicted by the parental expression levels. Expression quantitative trait loci (eQTL) mapping identified 30,774 eQTL including 96 trans-eQTL “hotspots,” each of which regulates the expression of a large number of genes. Interestingly, genes regulated by a trans-eQTL hotspot tend to be enriched for a specific function or act in the same genetic pathway. Also, genomic structural variation appeared to contribute to cis-regulation of gene expression. Besides genes showing Mendelian inheritance in the RIL population, we also found genes whose expression level and variation in the progeny could not be predicted based on parental difference, indicating that non-Mendelian factors also contribute to expression variation. Specifically, we found 145 genes that show patterns of expression reminiscent of paramutation such that all the progeny had expression levels similar to one of the two parents. Furthermore, we identified another 210 genes that exhibited unexpected patterns of transcript presence/absence. Many of these genes are likely to be gene fragments resulting from transposition, and the presence/absence of their transcripts could influence expression levels of their ancestral syntenic genes. Overall, our results contribute to the identification of novel expression patterns and broaden the understanding of transcriptional variation in plants.

Phenotypes are determined by the expression of genes, the environment, and the interaction of gene expression and the environment. However, a complete understanding of the inheritance of and genome-wide regulation of gene expression is lacking. One approach, called expression quantitative trait locus (eQTL) mapping provides the opportunity to examine the genome-wide inheritance and regulation of gene expression. In this paper, we conducted high-throughput sequencing of gene transcripts to examine gene expression in the shoot apex of a maize biparental mapping population. We quantified expression levels from 28,603 genes in the population and showed that the vast majority of genes exhibited the expected pattern of Mendelian inheritance. We genetically mapped the expression patterns and identified genomic regions associated with gene expression. Notably, we detected gene expression patterns that exhibited non-Mendelian inheritance. These included 145 genes that exhibited expression patterns in the progeny that were similar to only one of the parents and 210 genes with unexpected presence/absence expression patterns. The findings of non-Mendelian inheritance underscore the complexity of gene expression and provide a framework for understanding these complexities.

Progress toward understanding heterosis in crop plants

Although heterosis, or hybrid vigor, is widely exploited in agriculture, a complete description of its molecular underpinnings has remained elusive despite extensive investigation. It appears that there is not a single, simple explanation for heterosis. Instead, it is likely that heterosis arises in crosses between genetically distinct individuals as a result of a diversity of mechanisms. Heterosis generally results from the action of multiple loci, and different loci affect heterosis for different traits and in different hybrids. Hence, multigene models are likely to prove most informative for understanding heterosis. Complementation of allelic variation, as well as complementation of variation in gene content and gene expression patterns, is likely to be an important contributor to heterosis. Epigenetic variation has the potential to interact in hybrid genotypes via novel mechanisms. Several other intriguing hypotheses are also under investigation. In crops, heterosis must be considered within the context of the genomic impacts of prior selection for agronomic traits.

Heterosis: emerging ideas about hybrid vigour

Perceived by Charles Darwin in many vegetable plants and rediscovered by George H Shull and Edward M East in maize, heterosis or hybrid vigour is one of the most widely utilized phenomena, not only in agriculture but also in animal breeding. Although, numerous studies have been carried out to understand its genetic and/or molecular basis in the past 100 years, our knowledge of the underlying molecular processes that results in hybrid vigour can best be defined as superficial. Even after century long deliberations, there is no consensus on the relative/individual contribution of the genetic/epigenetic factors in the manifestation of heterosis. However, with the recent advancements in functional genomics, transcriptomics, proteomics, and metabolomics-related technologies, the riddle of heterosis is being reinvestigated by adopting systems-level approaches to understand the underlying molecular mechanisms. A number of intriguing hypotheses are converging towards the idea of a cumulative positive effect of the differential expression of a variety of genes, on one or several yield-affecting metabolic pathways or overall energy-use efficiency, as the underlying mechanism for the manifestation of heterosis. Presented here is a brief account of clues gathered from various investigative approaches targeted towards better scientific understanding of this process.

Complementation contributes to transcriptome complexity in maize (Zea mays L.) hybrids relative to their inbred parents

Typically, F₁-hybrids are more vigorous than their homozygous, genetically distinct parents, a phenomenon known as heterosis. In the present study, the transcriptomes of the reciprocal maize (Zea mays L.) hybrids B73×Mo17 and Mo17×B73 and their parental inbred lines B73 and Mo17 were surveyed in primary roots, early in the developmental manifestation of heterotic root traits. The application of statistical methods and a suitable experimental design established that 34,233 (i.e., 86%) of all high-confidence maize genes were expressed in at least one genotype. Nearly 70% of all expressed genes were differentially expressed between the two parents and 42%–55% of expressed genes were differentially expressed between one of the parents and one of the hybrids. In both hybrids, ∼10% of expressed genes exhibited nonadditive gene expression. Consistent with the dominance model (i.e., complementation) for heterosis, 1124 genes that were expressed in the hybrids were expressed in only one of the two parents. For 65 genes, it could be shown that this was a consequence of complementation of genomic presence/absence variation. For dozens of other genes, alleles from the inactive inbred were activated in the hybrid, presumably via interactions with regulatory factors from the active inbred. As a consequence of these types of complementation, both hybrids expressed more genes than did either parental inbred. Finally, in hybrids, ∼14% of expressed genes exhibited allele-specific expression (ASE) levels that differed significantly from the parental-inbred expression ratios, providing further evidence for interactions of regulatory factors from one parental genome with target genes from the other parental genome.

Specific changes in total and mitochondrial proteomes are associated with higher levels of heterosis in maize hybrids

The phenomenon of hybrid vigor (heterosis) has long been harnessed by plant breeders to improve world food production. However, the changes that are essential for heterotic responses and the mechanisms responsible for heterosis remain undefined. Large increases in biomass and yield in high-heterosis hybrids suggest that alterations in bioenergetic processes may contribute to heterosis. Progeny from crosses between various inbred lines vary in the extent of vigor observed. Field-grown maize F₁ hybrids that consistently exhibited either low or high heterosis across a variety of environments were examined for changes in proteins that may be correlated with increased plant vigor and yield. Unpollinated ears at the time of flowering (ear shoots) were selected for the studies because they are metabolically active, rich in mitochondria, and the sizes of the ears are diagnostic of yield heterosis. Total protein and mitochondrial proteomes were compared among low- and higher-heterosis hybrids. Two-dimensional difference gel electrophoresis was used to identify allelic and/or isoform differences linked to heterosis. Identification of differentially regulated spots by mass spectrometry revealed proteins involved in stress responses as well as primary carbon and protein metabolism. Many of these proteins were identified in multiple spots, but analysis of their abundances by label-free mass spectrometry suggested that most of the expression differences were due to isoform variation rather than overall protein amount. Thus, our proteomics studies suggest that expression of specific alleles and/or post-translational modification of specific proteins correlate with higher levels of heterosis.

Heterosis manifestation during early Arabidopsis seedling development is characterized by intermediate gene expression and enhanced metabolic activity in the hybrids

Heterosis-associated cellular and molecular processes were analyzed in seeds and seedlings of Arabidopsis thaliana accessions Col-0 and C24 and their heterotic hybrids. Microscopic examination revealed no advantages in terms of hybrid mature embryo organ sizes or cell numbers. Increased cotyledon sizes were detectable 4 days after sowing. Growth heterosis results from elevated cell sizes and numbers, and is well established at 10 days after sowing. The relative growth rates of hybrid seedlings were most enhanced between 3 and 4 days after sowing. Global metabolite profiling and targeted fatty acid analysis revealed maternal inheritance patterns for a large proportion of metabolites in the very early stages. During developmental progression, the distribution shifts to dominant, intermediate and heterotic patterns, with most changes occurring between 4 and 6 days after sowing. The highest incidence of heterotic patterns coincides with establishment of size differences at 4 days after sowing. In contrast, overall transcript patterns at 4, 6 and 10 days after sowing are characterized by intermediate to dominant patterns, with parental transcript levels showing the largest differences. Overall, the results suggest that, during early developmental stages, intermediate gene expression and higher metabolic activity in the hybrids compared to the parents lead to better resource efficiency, and therefore enhanced performance in the hybrids.

Characterization of heterotic quantitative trait loci in maize by evaluation of near-isogenic lines and their crosses at two competition levels

In a previous study on a maize (Zea mays L.) population of recombinant inbreds derived from B73 × H99, we identified several quantitative trait loci (QTL) for agronomic traits with high dominance-additive ratio. Then, for four of these QTL, we developed families of near-isogenic lines (NILs) homozygous either for the QTL allele from B73 (BB) or from H99 (HH); for two of these QTL, the NILs' families were produced in two different genetic backgrounds. The present study was conducted to: (1) characterize these QTL for agronomic traits and (2) verify whether their effects were influenced by the genetic background, inbreeding level and plant density (PD). The six NILs' families were tested across 3 years and in three experiments at different inbreeding levels as NILs per se and their reciprocal crosses (Experiment 1), NILs crossed to related inbreds B73 and H99 (Experiment 2) and NILs crossed to four unrelated inbreds (Experiment 3). Experiment 2 was conducted at two PDs (4.5 and 9.0 plants m(-2)). Results of Experiments 1 and 2 confirmed previous findings as to QTL effects, with dominance-additive ratio superior to 1 for several traits; as a tendency, dominance effects were more pronounced in Experiment 1. The QTL effects were also confirmed in Experiment 3. The interactions involving QTL effects, families and PD were generally negligible, suggesting a certain stability of the QTL. Results emphasize the importance of dominance effects for these QTL, suggesting that they might deserve further studies, using the NILs' families and their crosses as base materials.

A unifying theory for general multigenic heterosis: energy efficiency, protein metabolism, and implications for molecular breeding

Hybrids between genetically diverse varieties display enhanced growth, and increased total biomass, stress resistance and grain yield. Gene expression and metabolic studies in maize, rice and other species suggest that protein metabolism plays a role in the growth differences between hybrids and inbreds. Single trait heterosis can be explained by the existing theories of dominance, overdominance and epistasis. General multigenic heterosis is observed in a wide variety of different species and is likely to share a common underlying biological mechanism. This review presents a model to explain differences in growth and yield caused by general multigenic heterosis. The model describes multigenic heterosis in terms of energy-use efficiency and faster cell cycle progression where hybrids have more efficient growth than inbreds because of differences in protein metabolism. The proposed model is consistent with the observed variation of gene expression in different pairs of inbred lines and hybrid offspring as well as growth differences in polyploids and aneuploids. It also suggests an approach to enhance yield gains in both hybrid and inbred crops via the creation of an appropriate computational analysis pipeline coupled to an efficient molecular breeding program.

Nonadditive protein accumulation patterns in Maize (Zea mays L.) hybrids during embryo development

Heterosis describes the superior performance of heterozygous F(1)-hybrid plants compared to their homozygous parental inbred lines. In the present study, heterosis was detected for length, weight, and the time point of seminal root primordia initiation in maize (Zea mays L.) embryos of the reciprocal F(1)-hybrids UH005xUH250 and UH250xUH005. A two-dimensional gel electrophoresis (2-DE) proteome survey of the most abundant proteins of the reciprocal hybrids and their parental inbred lines 25 and 35 days after pollination revealed that 141 of 597 detected proteins (24%) exhibited nonadditive accumulation in at least one hybrid. Approximately 44% of all nonadditively accumulated proteins displayed an expression pattern that was not distinguishable from the low parent value. Electrospray ionization-tandem mass spectrometry (ESI-MS/MS) analyses and subsequent functional classification of the 141 proteins revealed that development, protein metabolism, redox-regulation, glycolysis, and amino acid metabolism were the most prominent functional classes among nonadditively accumulated proteins. In 35-day-old embryos of the hybrid UH250xUH005, a significant up-regulation of enzymes related to glucose metabolism which often exceeded the best parent values was observed. A comparison of nonadditive protein accumulation between rice and maize embryo data sets revealed a significant overlap of nonadditively accumulated proteins suggesting conserved organ- or tissue-specific regulatory mechanisms in monocots related to heterosis.

Heterosis in early seed development: a comparative study of F1 embryo and endosperm tissues 6 days after fertilization

Heterosis specifies the superior performance of heterozygous individuals and although used in plant breeding the underlying molecular mechanisms still remain largely elusive. In this study, we demonstrate the manifestation of heterosis in hybrid maize embryo and endosperm tissue 6 days after fertilization in crosses of several inbred lines. We provide a comparative analysis of heterosis-associated gene expression in these tissues by a combined approach of suppression subtractive hybridization and microarray hybridizations. Non-additive expression pattern indicated a trans-regulatory mechanism to act early after fertilization in hybrid embryo and endosperm although the majority of genes showed mid-parental expression levels in embryo and dosage dependent expression levels in endosperm. The consistent expression pattern within both tissues and both inbred line genotype combinations of genes coding for chromatin related proteins pointed to heterosis-related epigenetic processes. These and genes involved in other biological processes, identified in this study, might provide entry points for the investigation of regulatory networks associated with the specification of heterosis.

Gene expression analysis at the intersection of ploidy and hybridity in maize

Heterosis and polyploidy are two important aspects of plant evolution. To examine these issues, we conducted a global gene expression study of a maize ploidy series as well as a set of tetraploid inbred and hybrid lines. This gene expression analysis complements an earlier phenotypic study of these same materials. We find that ploidy change affects a large fraction of the genome, albeit at low levels; gene expression changes rarely exceed 2-fold and are typically not statistically significant. The most common gene expression profile we detected is greater than linear increase from monoploid to diploid, and reductions from diploid to triploid and from triploid to tetraploid, a trend that mirrors plant stature. When examining heterosis in tetraploid maize lines, we found a large fraction of the genome impacted but the majority of changes were not statistically significant at 2-fold or less. Non-additive expression was common in the hybrids, and the extent of non-additivity increased both in number and magnitude from duplex to quadruplex hybrids. Overall, we find that gene expression trends mirror observations from the phenotypic studies; however, obvious mechanistic connections remain unknown.

Correlation between parental transcriptome and field data for the characterization of heterosis in Zea mays L

Heterosis is widely exploited in plant breeding, although its molecular basis is still not fully understood. For the characterization of this phenomenon and the development of transcriptome-based methods to predict hybrid performance (HP), we applied a microarray (46k) analysis of 21 European maize (Zea mays L.), 14 dent and 7 flint parental inbred lines. Expression profiles of the parental inbreds at the seedling stage were correlated with grain yield (GY) and grain dry matter content (GDMC) of 98 flint x dent factorial crosses at six locations. We observed highly significant correlations of the parental expression levels of certain differentially expressed genes with heterosis and HP for GY and also with HP for GDMC. This strong correlation provided first evidence toward a prediction potential of the genes and their expression levels. The identified gene set based on the parental transcriptome data revealed functional characteristics of HP and heterosis. Gene ontology (GO) analyses were performed to compare genes correlated with their expression pattern to HP for GY and GDMC, respectively. Between these gene groups, mostly different functional classes of genes were found to be enriched or underrepresented. The phenomenon of heterosis was characterized by the over- and underrepresentation of specific GO terms among heterosis-correlated genes.

Background correction of two-colour cDNA microarray data using spatial smoothing methods

The analysis of two-colour cDNA microarray data usually involves subtracting background values from foreground values prior to normalization and further analysis. This approach has the advantage of reducing bias and the disadvantage of blowing up the variance of lower abundant spots. Whenever background subtraction is considered, it implicitly assumes locally constant background values. In practice, this assumption is often not met, which casts doubts on the usefulness of simple background subtraction. In order to improve background correction, we propose local background smoothing within the pre-processing pipeline of cDNA microarray data prior to background correction. For this purpose, we employ a geostatistical framework with ordinary kriging using both isotropic and anisotropic models of spatial correlation and 2-D locally weighted regression. We show that application of local background smoothing prior to background correction is beneficial in comparison to using raw background estimates. This is done using data of a self-versus-self experiment in Arabidopsis where subsets of differentially expressed genes were simulated. Using locally smoothed background values in conjunction with existing background correction methods increases the power, increases the accuracy and decreases the number of false positive results.

Intraspecific hybrids of Arabidopsis thaliana revealed no gross alterations in endopolyploidy, DNA methylation, histone modifications and transcript levels

Arabidopsis accessions Col-0 and C24 and their reciprocal hybrids were employed as a model system to investigate the potential relationship between changes in DNA methylation, chromatin structure, endopolyploidization and gene expression in heterotic genotypes. Nucleolus size, endopolyploidization level and distribution of DNA and histone H3 methylation at the microscopic level does not differ between parents and their hybrids. Methylation sensitive amplified polymorphism revealed a largely constant pattern of DNA methylation (97% of signals analyzed) after intraspecific crosses. The parental expression profile of selected genes was maintained in hybrid offspring. No correlation was found between expression pattern and DNA methylation levels at restriction sites within 5' regulatory regions. Thus, the results revealed only minor changes of chromatin properties and other nuclear features in response to intraspecific hybridization in Arabidopsis thaliana.

Molecular dissection of heterosis manifestation during early maize root development

Heterosis is of paramount agronomic importance and has been successfully exploited in maize hybrid breeding for decades. Nevertheless, the molecular basis of heterosis remains elusive. Heterosis is not only observed in adult traits like yield or plant height, but is already detected during embryo and seedling development. Hence, the maize (Zea mays L.) primary root which is the first organ that emerges after germination is a suitable model to study heterosis manifestation. Various seedling root traits including primary root length and lateral root density display heterosis. Microarray studies suggest organ specific patterns of nonadditive gene expression in maize hybrids. Moreover, such experiments support the notion that global expression trends in maize primary roots are conserved between different hybrids. Furthermore, nonadditive expression patterns of specific genes such as a SUPEROXIDE DISMUTASE 2 might contribute to the early manifestation of heterosis. Proteome profiling experiments of maize hybrid primary roots revealed nonadditive accumulation patterns that were distinct from the corresponding RNA profiles underscoring the importance of posttranscriptional processes such as protein modifications that might be related to heterosis. Finally, analysis of selected metabolites imply that a subtle regulation of particular biochemical pathways such as the phenylpropanoid pathway in hybrids might contribute to the manifestation of heterosis in maize primary roots. In the future, recently developed molecular tools will facilitate the analysis of the molecular principles underlying heterosis in maize roots.

The genetics of inbreeding depression

Inbreeding depression - the reduced survival and fertility of offspring of related individuals - occurs in wild animal and plant populations as well as in humans, indicating that genetic variation in fitness traits exists in natural populations. Inbreeding depression is important in the evolution of outcrossing mating systems and, because intercrossing inbred strains improves yield (heterosis), which is important in crop breeding, the genetic basis of these effects has been debated since the early twentieth century. Classical genetic studies and modern molecular evolutionary approaches now suggest that inbreeding depression and heterosis are predominantly caused by the presence of recessive deleterious mutations in populations.

Comparative gene expression profiles between heterotic and non-heterotic hybrids of tetraploid Medicago sativa

BACKGROUND

Heterosis, the superior performance of hybrids relative to parents, has clear agricultural value, but its genetic control is unknown. Our objective was to test the hypotheses that hybrids expressing heterosis for biomass yield would show more gene expression levels that were different from midparental values and outside the range of parental values than hybrids that do not exhibit heterosis.

RESULTS

We tested these hypotheses in three Medicago sativa (alfalfa) genotypes and their three hybrids, two of which expressed heterosis for biomass yield and a third that did not, using Affymetrix M. truncatula GeneChip arrays. Alfalfa hybridized to approximately 47% of the M. truncatula probe sets. Probe set signal intensities were analyzed using MicroArray Suite v.5.0 (MAS) and robust multi-array average (RMA) algorithms. Based on MAS analysis, the two heterotic hybrids performed similarly, with about 27% of genes showing differential expression among the parents and their hybrid compared to 12.5% for the non-heterotic hybrid. At a false discovery rate of 0.15, 4.7% of differentially expressed genes in hybrids (approximately 300 genes) showed nonadditive expression compared to only 0.5% (16 genes) in the non-heterotic hybrid. Of the nonadditively expressed genes, approximately 50% showed expression levels that fell outside the parental range in heterotic hybrids, but only one of 16 showed a similar profile in the non-heterotic hybrid. Genes whose expression differed in the parents were three times more likely to show nonadditive expression than genes whose parental transcript levels were equal.

CONCLUSION

The higher proportions of probe sets with expression level that differed from the parental midparent value and that were more extreme than either parental value in the heterotic hybrids compared to a non-heterotic hybrid were also found using RMA. We conclude that nonadditive expression of transcript levels may contribute to heterosis for biomass yield in alfalfa.