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

Study on Allele Specific Expression of Long-Term Residents in High Altitude Areas

In diploid organisms, half of the chromosomes in each cell come from the father and half from the mother. Through previous studies, it was found that the paternal chromosome and the maternal chromosome can be regulated and expressed independently, leading to the emergence of allele specific expression (ASE). In this study, we analyzed the differential expression of alleles in the high-altitude population and the normal population based on the RNA sequencing data. Through gene cluster analysis and protein interaction network analysis, we found some changes occurred at the gene level, and some negative effects. During the study, we realized that the calmodulin homology domain may have a certain correlation with long-term survival at high altitude. The plateau environment is characterized by hypoxia, low air pressure, strong ultraviolet radiation, and low temperature. Accordingly, the genetic changes in the process of adaptation are mainly reflected in these characteristics. High altitude generation living is also highly related to cancer, immune disease, cardiovascular disease, neurological disease, endocrine disease, and other diseases. Therefore, the medical system in high altitude areas should pay more attention to these diseases.

High-quality genome assembly enables prediction of allele-specific gene expression in hybrid poplar

Poplar (Populus) is a well-established model system for tree genomics and molecular breeding, and hybrid poplar is widely used in forest plantations. However, distinguishing its diploid homologous chromosomes is difficult, complicating advanced functional studies on specific alleles. In this study, we applied a trio-binning design and PacBio high-fidelity long-read sequencing to obtain haplotype-phased telomere-to-telomere genome assemblies for the 2 parents of the well-studied F1 hybrid "84K" (Populus alba × Populus tremula var. glandulosa). Almost all chromosomes, including the telomeres and centromeres, were completely assembled for each haplotype subgenome apart from 2 small gaps on one chromosome. By incorporating information from these haplotype assemblies and extensive RNA-seq data, we analyzed gene expression patterns between the 2 subgenomes and alleles. Transcription bias at the subgenome level was not uncovered, but extensive-expression differences were detected between alleles. We developed machine-learning (ML) models to predict allele-specific expression (ASE) with high accuracy and identified underlying genome features most highly influencing ASE. One of our models with 15 predictor variables achieved 77% accuracy on the training set and 74% accuracy on the testing set. ML models identified gene body CHG methylation, sequence divergence, and transposon occupancy both upstream and downstream of alleles as important factors for ASE. Our haplotype-phased genome assemblies and ML strategy highlight an avenue for functional studies in Populus and provide additional tools for studying ASE and heterosis in hybrids.

Strand-Specific RNA Sequencing Reveals Gene Expression Patterns in F1 Chick Breast Muscle and Liver after Hatching

Heterosis refers to the phenomenon where hybrids exhibit superior performance compared to the parental phenotypes and has been widely utilized in crossbreeding programs for animals and crops, yet the molecular mechanisms underlying this phenomenon remain enigmatic. A better understanding of the gene expression patterns in post-hatch chickens is very important for exploring the genetic basis underlying economically important traits in the crossbreeding of chickens. In this study, breast muscle and liver tissues (n = 36) from full-sib F1 birds and their parental pure lines were selected to identify gene expression patterns and differentially expressed genes (DEGs) at 28 days of age by strand-specific RNA sequencing (ssRNA-seq). This study indicates that additivity is the predominant gene expression pattern in the F1 chicken post-hatch breast muscle (80.6% genes with additivity) and liver (94.2% genes with additivity). In breast muscle, Gene Ontology (GO) enrichment analysis revealed that a total of 11 biological process (BP) terms closely associated with growth and development were annotated in the identified DEG sets and non-additive gene sets, including STAT5A, TGFB2, FGF1, IGF2, DMA, FGF16, FGF12, STAC3, GSK3A, and GRB2. Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation presented that a total of six growth- and development-related pathways were identified, involving key genes such as SLC27A4, GLUL, TGFB2, COX17, and GSK3A, including the PPAR signaling pathway, TGF-beta signaling pathway, and mTOR signaling pathway. Our results may provide a theoretical basis for crossbreeding in domestic animals.

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.

Increased long-distance and homo-trans interactions related to H3K27me3 in Arabidopsis hybrids

In plants, the genome structure of hybrids changes compared with their parents, but the effects of these changes in hybrids remain elusive. Comparing reciprocal crosses between Col × C24 and C24 × Col in Arabidopsis using high-throughput chromosome conformation capture assay (Hi-C) analysis, we found that hybrid three-dimensional (3D) chromatin organization had more long-distance interactions relative to parents, and this was mainly located in promoter regions and enriched in genes with heterosis-related pathways. The interactions between euchromatin and heterochromatin were increased, and the compartment strength decreased in hybrids. In compartment domain (CD) boundaries, the distal interactions were more in hybrids than their parents. In the hybrids of CURLY LEAF (clf) mutants clfCol  × clfC24 and clfC24  × clfCol , the heterosis phenotype was damaged, and the long-distance interactions in hybrids were fewer than in their parents with lower H3K27me3. ChIP-seq data revealed higher levels of H3K27me3 in the region adjacent to the CD boundary and the same interactional homo-trans sites in the wild-type (WT) hybrids, which may have led to more long-distance interactions. In addition, the differentially expressed genes (DEGs) located in the boundaries of CDs and loop regions changed obviously in WT, and the functional enrichment for DEGs was different between WT and clf in the long-distance interactions and loop regions. Our findings may therefore propose a new epigenetic explanation of heterosis in the Arabidopsis hybrids and provide new insights into crop breeding and yield increase.

Pervasive under-dominance in gene expression underlying emergent growth trajectories in Arabidopsis thaliana hybrids

BACKGROUND

Complex traits, such as growth and fitness, are typically controlled by a very large number of variants, which can interact in both additive and non-additive fashion. In an attempt to gauge the relative importance of both types of genetic interactions, we turn to hybrids, which provide a facile means for creating many novel allele combinations.

RESULTS

We focus on the interaction between alleles of the same locus, i.e., dominance, and perform a transcriptomic study involving 141 random crosses between different accessions of the plant model species Arabidopsis thaliana. Additivity is rare, consistently observed for only about 300 genes enriched for roles in stress response and cell death. Regulatory rare-allele burden affects the expression level of these genes but does not correlate with F1 rosette size. Non-additive, dominant gene expression in F1 hybrids is much more common, with the vast majority of genes (over 90%) being expressed below the parental average. Unlike in the additive genes, regulatory rare-allele burden in the dominant gene set is strongly correlated with F1 rosette size, even though it only mildly covaries with the expression level of these genes.

CONCLUSIONS

Our study underscores under-dominance as the predominant gene action associated with emergence of rosette growth trajectories in the A. thaliana hybrid model. Our work lays the foundation for understanding molecular mechanisms and evolutionary forces that lead to dominance complementation of rare regulatory alleles.

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.

A complete telomere-to-telomere assembly of the maize genome

A complete telomere-to-telomere (T2T) finished genome has been the long pursuit of genomic research. Through generating deep coverage ultralong Oxford Nanopore Technology (ONT) and PacBio HiFi reads, we report here a complete genome assembly of maize with each chromosome entirely traversed in a single contig. The 2,178.6 Mb T2T Mo17 genome with a base accuracy of over 99.99% unveiled the structural features of all repetitive regions of the genome. There were several super-long simple-sequence-repeat arrays having consecutive thymine-adenine-guanine (TAG) tri-nucleotide repeats up to 235 kb. The assembly of the entire nucleolar organizer region of the 26.8 Mb array with 2,974 45S rDNA copies revealed the enormously complex patterns of rDNA duplications and transposon insertions. Additionally, complete assemblies of all ten centromeres enabled us to precisely dissect the repeat compositions of both CentC-rich and CentC-poor centromeres. The complete Mo17 genome represents a major step forward in understanding the complexity of the highly recalcitrant repetitive regions of higher plant genomes.

Gene expression variations and allele-specific expression of two rice and their hybrid in caryopses at single-nucleus resolution

Seeds are an indispensable part of the flowering plant life cycle and a critical determinant of agricultural production. Distinct differences in the anatomy and morphology of seeds separate monocots and dicots. Although some progress has been made with respect to understanding seed development in Arabidopsis, the transcriptomic features of monocotyledon seeds at the cellular level are much less understood. Since most important cereal crops, such as rice, maize, and wheat, are monocots, it is essential to study transcriptional differentiation and heterogeneity during seed development at a finer scale. Here, we present single-nucleus RNA sequencing (snRNA-seq) results of over three thousand nuclei from caryopses of the rice cultivars Nipponbare and 9311 and their intersubspecies F₁ hybrid. A transcriptomics atlas that covers most of the cell types present during the early developmental stage of rice caryopses was successfully constructed. Additionally, novel specific marker genes were identified for each nuclear cluster in the rice caryopsis. Moreover, with a focus on rice endosperm, the differentiation trajectory of endosperm subclusters was reconstructed to reveal the developmental process. Allele-specific expression (ASE) profiling in endosperm revealed 345 genes with ASE (ASEGs). Further pairwise comparisons of the differentially expressed genes (DEGs) in each endosperm cluster among the three rice samples demonstrated transcriptional divergence. Our research reveals differentiation in rice caryopsis from the single-nucleus perspective and provides valuable resources to facilitate clarification of the molecular mechanism underlying caryopsis development in rice and other monocots.

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.

Characterization of Genes That Exhibit Genotype-Dependent Allele-Specific Expression and Its Implications for the Development of Maize Kernel

Heterosis or hybrid vigor refers to the superior phenotypic traits of hybrids relative to their parental inbred lines. An imbalance between the expression levels of two parental alleles in the F1 hybrid has been suggested as a mechanism of heterosis. Here, based on genome-wide allele-specific expression analysis using RNA sequencing technology, 1689 genes exhibiting genotype-dependent allele-specific expression (genotype-dependent ASEGs) were identified in the embryos, and 1390 genotype-dependent ASEGs in the endosperm, of three maize F1 hybrids. Of these ASEGs, most were consistent in different tissues from one hybrid cross, but nearly 50% showed allele-specific expression from some genotypes but not others. These genotype-dependent ASEGs were mostly enriched in metabolic pathways of substances and energy, including the tricarboxylic acid cycle, aerobic respiration, and energy derivation by oxidation of organic compounds and ADP binding. Mutation and overexpression of one ASEG affected kernel size, which indicates that these genotype-dependent ASEGs may make important contributions to kernel development. Finally, the allele-specific methylation pattern on genotype-dependent ASEGs indicated that DNA methylation plays a potential role in the regulation of allelic expression for some ASEGs. In this study, a detailed analysis of genotype-dependent ASEGs in the embryo and endosperm of three different maize F1 hybrids will provide an index of genes for future research on the genetic and molecular mechanism of heterosis.

Overdominant expression of genes plays a key role in root growth of tobacco hybrids

Heterosis has greatly improved the yield and quality of crops. However, previous studies often focused on improving the yield and quality of the shoot system, while research on the root system was neglected. We determined the root numbers of 12 F₁ hybrids, all of which showed strong heterosis, indicating that tobacco F₁ hybrids have general heterosis. To understand its molecular mechanism, we selected two hybrids with strong heterosis, GJ (G70 × Jiucaiping No.2) and KJ (K326 × Jiucaiping No.2), and their parents for transcriptome analysis. There were 84.22% and 90.25% of the differentially expressed genes were overdominantly expressed. The enrichment analysis of these overdominantly expressed genes showed that "Plant hormone signal transduction", "Phenylpropanoid biosynthesis", "MAPK signaling pathway - plant", and "Starch and sucrose metabolism" pathways were associated with root development. We focused on the analysis of the biosynthetic pathways of auxin(AUX), cytokinins(CTK), abscisic acid(ABA), ethylene(ET), and salicylic acid(SA), suggesting that overdominant expression of these hormone signaling pathway genes may enhance root development in hybrids. In addition, Nitab4.5_0011528g0020、Nitab4.5_0003282g0020、Nitab4.5_0004384g0070 may be the genes involved in root growth. Genome-wide comparative transcriptome analysis enhanced our understanding of the regulatory network of tobacco root development and provided new ideas for studying the molecular mechanisms of tobacco root development.

Genome-wide analysis of transcriptome and histone modifications in Brassica napus hybrid

Although utilization of heterosis has largely improved the yield of many crops worldwide, the underlying molecular mechanism of heterosis, particularly for allopolyploids, remains unclear. Here, we compared epigenome and transcriptome data of an elite hybrid and its parental lines in three assessed tissues (seedling, flower bud, and silique) to explore their contribution to heterosis in allopolyploid B. napus. Transcriptome analysis illustrated that a small proportion of non-additive genes in the hybrid compared with its parents, as well as parental expression level dominance, might have a significant effect on heterosis. We identified histone modification (H3K4me3 and H3K27me3) variation between the parents and hybrid, most of which resulted from the differences between parents. H3K4me3 variations were positively correlated with gene expression differences among the hybrid and its parents. Furthermore, H3K4me3 and H3K27me3 were rather stable in hybridization and were mainly inherited additively in the B. napus hybrid. Together, our data revealed that transcriptome reprogramming and histone modification remodeling in the hybrid could serve as valuable resources for better understanding heterosis in allopolyploid crops.

Allele-specific expression and chromatin accessibility contribute to heterosis in tea plants (Camellia sinensis)

Heterosis is extensively used to improve crop productivity, yet its allelic and chromatin regulation remains unclear. Based on our resolved genomes of the maternal TGY and paternal HD, we analyzed the contribution of allele-specific expression (ASE) and chromatin accessibility of JGY and HGY, the artificial hybrids of oolong tea with the largest cultivated area in China. The ASE genes (ASEGs) of tea hybrids with maternal-biased were mainly related to the energy and terpenoid metabolism pathways, whereas the ASEGs with paternal-biased tend to be enriched in glutathione metabolism, and these parental bias of hybrids may coordinate and lead to the acquisition of heterosis in more biological pathways. ATAC-seq results showed that hybrids have significantly higher accessible chromatin regions (ACRs) compared with their parents, which may confer broader and stronger transcriptional activity of genes in hybrids. The number of ACRs with significantly increased accessibility in hybrids was much greater than decreased, and the associated alleles were also affected by differential ACRs across different parents, suggesting enhanced positive chromatin regulation and potential genetic effects in hybrids. Core ASEGs of terpene and purine alkaloid metabolism pathways with significant positive heterosis have greater chromatin accessibility in hybrids, and were potentially regulated by several members of the MYB, DOF and TRB families. The binding motif of CsMYB85 in the promoter ACR of the rate-limiting enzyme CsDXS was verified by DAP-seq. These results suggest that higher numbers and more accessible ACRs in hybrids contribute to the regulation of ASEGs, thereby affecting the formation of heterotic metabolites.

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.

Solving the grand challenge of phenotypic integration: allometry across scales

Phenotypic integration is a concept related to the cascade of trait relationships from the lowest organizational levels, i.e. genes, to the highest, i.e. whole-organism traits. However, the cause-and-effect linkages between traits are notoriously difficult to determine. In particular, we still lack a mathematical framework to model the relationships involved in the integration of phenotypic traits. Here, we argue that allometric models developed in ecology offer testable mathematical equations of trait relationships across scales. We first show that allometric relationships are pervasive in biology at different organizational scales and in different taxa. We then present mechanistic models that explain the origin of allometric relationships. In addition, we emphasized that recent studies showed that natural variation does exist for allometric parameters, suggesting a role for genetic variability, selection and evolution. Consequently, we advocate that it is time to examine the genetic determinism of allometries, as well as to question in more detail the role of genome size in subsequent scaling relationships. More broadly, a possible-but so far neglected-solution to understand phenotypic integration is to examine allometric relationships at different organizational levels (cell, tissue, organ, organism) and in contrasted species.

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.

Full-Length Transcriptome Reconstruction Reveals the Genetic Mechanisms of Eyestalk Displacement and Its Potential Implications on the Interspecific Hybrid Crab (Scylla serrata ♀ × S. paramamosain ♂)

Simple Summary

The eyestalk is a key organ in crustaceans that produces neurohormones and regulates a range of physiological functions. Eyestalk displacement was discovered in some first-generation (F1) offspring of the novel interspecific hybrid crab (Scylla serrata ♀ × S. paramamosain ♂). To uncover the genetic mechanism underlying eyestalk displacement and its potential implications, high-quality transcriptome was reconstructed using single-molecule real-time (SMRT) sequencing. A total of 37 significantly differential alternative splicing (DAS) events (17 up-regulated and 20 down-regulated) and 1475 significantly differential expressed transcripts (DETs) (492 up-regulated and 983 down-regulated) were detected in hybrid crabs with displaced eyestalks (DH). The most significant DAS events and DETs were annotated as being endoplasmic reticulum chaperone BiP and leucine-rich repeat protein lrrA-like isoform X2. In addition, the top ten significant gene ontology (GO) terms were related to the cuticle or chitin. Overall, this study highlights the underlying genetic mechanisms of eyestalk displacement and provide useful knowledge for mud crab (Scylla spp.) crossbreeding.

Abstract

The lack of high-quality juvenile crabs is the greatest impediment to the growth of the mud crab (Scylla paramamosain) industry. To obtain high-quality hybrid offspring, a novel hybrid mud crab (S. serrata ♀ × S. paramamosain ♂) was successfully produced in our previous study. Meanwhile, an interesting phenomenon was discovered, that some first-generation (F1) hybrid offspring’s eyestalks were displaced during the crablet stage I. To uncover the genetic mechanism underlying eyestalk displacement and its potential implications, both single-molecule real-time (SMRT) and Illumina RNA sequencing were implemented. Using a two-step collapsing strategy, three high-quality reconstructed transcriptomes were obtained from purebred mud crabs (S. paramamosain) with normal eyestalks (SPA), hybrid crabs with normal eyestalks (NH), and hybrid crabs with displaced eyestalks (DH). In total, 37 significantly differential alternative splicing (DAS) events (17 up-regulated and 20 down-regulated) and 1475 significantly differential expressed transcripts (DETs) (492 up-regulated and 983 down-regulated) were detected in DH. The most significant DAS events and DETs were annotated as being endoplasmic reticulum chaperone BiP and leucine-rich repeat protein lrrA-like isoform X2. In addition, the top ten significant GO terms were related to the cuticle or chitin. Overall, high-quality reconstructed transcriptomes were obtained for the novel interspecific hybrid crab and provided valuable insights into the genetic mechanisms of eyestalk displacement in mud crab (Scylla spp.) crossbreeding.

Transcriptome profiling of two super hybrid rice provides insights into the genetic basis of heterosis

BACKGROUND

Heterosis is a phenomenon that hybrids show superior performance over their parents. The successful utilization of heterosis has greatly improved rice productivity, but the molecular basis of heterosis remains largely unclear.

RESULTS

Here, the transcriptomes of young panicles and leaves of the two widely grown two-line super hybrid rice varieties (Jing-Liang-You-Hua-Zhan (JLYHZ) and Long-Liang-You-Hua-Zhan (LLYHZ)) and their parents were analyzed by RNA-seq. Transcriptome profiling of the hybrids revealed 1,778 ~ 9,404 differentially expressed genes (DEGs) in two tissues, which were identified by comparing with their parents. GO, and KEGG enrichment analysis showed that the pathways significantly enriched in both tissues of two hybrids were all related to yield and resistance, like circadian rhythm (GO:0,007,623), response to water deprivation (GO:0,009,414), and photosynthetic genes (osa00196). Allele-specific expression genes (ASEGs) were also identified in hybrids. The ASEGs were most significantly enriched in ionotropic glutamate receptor signaling pathway, which was hypothesized to be potential amino acid sensors in plants. Moreover, the ASEGs were also differentially expressed between parents. The number of variations in ASEGs is higher than expected, especially for large effect variations. The DEGs and ASEGs are the potential reasons for the formation of heterosis in the two elite super hybrid rice.

CONCLUSIONS

Our results provide a comprehensive understanding of the heterosis of two-line super hybrid rice and facilitate the exploitation of heterosis in hybrid rice breeding with high yield heterosis.

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.

Gene expression variation explains maize seed germination heterosis

BACKGROUND

Heterosis has been extensively utilized in plant breeding, however, the underlying molecular mechanism remains largely elusive. Maize (Zea mays), which exhibits strong heterosis, is an ideal material for studying heterosis.

RESULTS

In this study, there is faster imbibition and development in reciprocal crossing Zhengdan958 hybrids than in their parent lines during seed germination. To investigate the mechanism of heterosis of maize germination, comparative transcriptomic analyses were conducted. The gene expression patterns showed that 1324 (47.27%) and 1592 (66.44%) of the differential expression genes between hybrids and either parental line display parental dominance up or higher levels in the reciprocal cross of Zhengdan958, respectively. Notably, these genes were mainly enriched in metabolic pathways, including carbon metabolism, glycolysis/gluconeogenesis, protein processing in endoplasmic reticulum, etc. CONCLUSION: Our results provide evidence for the higher expression level genes in hybrid involved in metabolic pathways acting as main contributors to maize seed germinating heterosis. These findings provide new insights into the gene expression variation of maize embryos and improve the understanding of maize seed germination heterosis.

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.

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.

Transcriptomic Variations and Network Hubs Controlling Seed Size and Weight During Maize Seed Development

To elucidate the mechanisms underlying seed development in maize, comprehensive RNA-seq analyses were conducted on Zhengdan1002 (ZD1002), Zhengdan958 (ZD958), and their parental lines during seven seed developmental stages. We found that gene expression levels were largely nonadditive in hybrids and that cis-only or trans × cis pattern played a large role in hybrid gene regulation during seed developmental stage. Weighted gene co-expression network (WGCNA) analysis showed that 36 modules were highly correlated (r = -0.90-0.92, p < 0.05) with kernel weight, length, and width during seed development. Forty-five transcription factors and 38 ribosomal protein genes were identified as major hub genes determining seed size/weight. We also described a network hub, Auxin Response Factor 12 of maize (ZmARF12), a member of a family of transcription factor that mediate gene expression in response to auxin, potentially links auxin signal pathways, cell division, and the size of the seeds. The ZmARF12 mutant exhibited larger seed size and higher grain weight. ZmARF12 transcription was negatively associated with cell division during seed development, which was confirmed by evaluating the yield of protoplasts that isolated from the kernels of the mutant and other inbred lines. Transient knock-down of ZmARF12 in maize plants facilitated cell expansion and division, whereas transient silencing of its potential interactor ZmIAA8 impaired cell division. ZmIAA8 expression was repressed in the ZmARF12 over-expressed protoplasts. The mutant phenotype and the genetics studies presented here illustrated evidence that ZmARF12 is a cell division repressor, and potentially determines the final seed size.

Genome-wide analysis of deletions in maize population reveals abundant genetic diversity and functional impact

Two read depth methods were jointly used in next-generation sequencing data to identify deletions in maize population. GWAS by deletions were analyzed for gene expression pattern and classical traits, respectively. Many studies have confirmed that structural variation (SV) is pervasive throughout the maize genome. Deletion is one type of SV that may impact gene expression and cause phenotypic changes in quantitative traits. In this study, two read count approaches were used to analyze the deletions in the whole-genome sequencing data of 270 maize inbred lines. A total of 19,754 deletion windows overlapped 12,751 genes, which were unevenly distributed across the genome. The deletions explained population structure well and correlated with genomic features. The deletion proportion of genes was determined to be negatively correlated with its expression. The detection of gene expression quantitative trait loci (eQTL) indicated that local eQTL were fewer but had larger effects than distant ones. The common associated genes were related to basic metabolic processes, whereas unique associated genes with eQTL played a role in the stress or stimulus responses in multiple tissues. Compared with the eQTL detected by SNPs derived from the same sequencing data, 89.4% of the associated genes could be detected by both markers. The effect of top eQTL detected by SNPs was usually larger than that detected by deletions for the same gene. A genome-wide association study (GWAS) on flowering time and plant height illustrated that only a few loci could be consistently captured by SNPs, suggesting that combining deletion and SNP for GWAS was an excellent strategy to dissect trait architecture. Our findings will provide insights into characteristic and biological function of genome-wide deletions in maize.

Genome-wide differences in gene expression and alternative splicing in developing embryo and endosperm, and between F1 hybrids and their parental pure lines in sorghum

Developing embryo and endosperm of sorghum show substantial and multifaceted differences in gene expression and alternative splicing, which are potentially relevant to heterosis. Differential regulation of gene expression and alternative splicing (AS) are major molecular mechanisms dictating plant growth and development, as well as underpinning heterosis in F1 hybrids. Here, using deep RNA-sequencing we analyzed differences in genome-wide gene expression and AS between developing embryo and endosperm, and between F1 hybrids and their pure-line parents in sorghum. We uncover dramatic differences in both gene expression and AS between embryo and endosperm with respect to gene features and functions, which are consistent with the fundamentally different biological roles of the two tissues. Accordingly, F1 hybrids showed substantial and multifaceted differences in gene expression and AS compared with their pure-line parents, again with clear tissue specificities including extents of difference, genes involved and functional enrichments. Our results provide useful transcriptome resources as well as novel insights for further elucidation of seed yield heterosis in sorghum and related crops.

Advances in Research on the Mechanism of Heterosis in Plants

Heterosis is a common biological phenomenon in nature. It substantially contributes to the biomass yield and grain yield of plants. Moreover, this phenomenon results in high economic returns in agricultural production. However, the utilization of heterosis far exceeds the level of theoretical research on this phenomenon. In this review, the recent progress in research on heterosis in plants was reviewed from the aspects of classical genetics, parental genetic distance, quantitative trait loci, transcriptomes, proteomes, epigenetics (DNA methylation, histone modification, and small RNA), and hormone regulation. A regulatory network of various heterosis-related genes under the action of different regulatory factors was summarized. This review lays a foundation for the in-depth study of the molecular and physiological aspects of this phenomenon to promote its effects on increasing the yield of agricultural production.

Transcriptome Reveals Allele Contribution to Heterosis in Maize

Heterosis, which has greatly increased maize yields, is associated with gene expression patterns during key developmental stages that enhance hybrid phenotypes relative to parental phenotypes. Before heterosis can be more effectively used for crop improvement, hybrid maize developmental gene expression patterns must be better understood. Here, six maize hybrids, including the popular hybrid Zhengdan958 (ZC) from China, were studied. Maize hybrids created in-house were generated using an incomplete diallel cross (NCII)-based strategy from four elite inbred parental lines. Differential gene expression (DEG) profiles corresponding to three developmental stages revealed that hybrid partial expression patterns exhibited complementarity of expression of certain parental genes, with parental allelic expression patterns varying both qualitatively and quantitatively in hybrids. Single-parent expression (SPE) and parent-specific expression (PSE) types of qualitative variation were most prevalent, 43.73 and 41.07% of variation, respectively. Meanwhile, negative super-dominance (NSD) and positive super-dominance (PSD) types of quantitative variation were most prevalent, 31.06 and 24.30% of variation, respectively. During the early reproductive growth stage, the gene expression pattern differed markedly from other developmental stage patterns, with allelic expression patterns during seed development skewed toward low-value parental alleles in hybrid seeds exhibiting significant quantitative variation-associated superiority. Comparisons of qualitative gene expression variation rates between ZC and other hybrids revealed proportions of SPE-DEGs (41.36%) in ZC seed DEGs that significantly exceeded the average proportion of SPE-DEGs found in seeds of other hybrids (28.36%). Importantly, quantitative gene expression variation rate comparisons between ZC and hybrids, except for transgressive expression, revealed that the ZC rate exceeded the average rate for other hybrids, highlighting the importance of partial gene expression in heterosis. Moreover, enriched ZC DEGs exhibiting distinct tissue-specific expression patterns belonged to four biological pathways, including photosynthesis, plant hormone signal transduction, biology metabolism and biosynthesis. These results provide valuable technical insights for creating hybrids exhibiting strong heterosis.

Unraveling regulatory divergence, heterotic malleability, and allelic imbalance switching in rice due to drought stress

The indica ecotypes, IR64, an elite drought-susceptible variety adapted to irrigated ecosystem, and Apo (IR55423-01 or NSIC RC9), a moderate drought-tolerant upland genotype together with their hybrid (IR64 × Apo) were exposed to non- and water-stress conditions. By sequencing (RNA-seq) these genotypes, we were able to map genes diverging in cis and/or trans factors. Under non-stress condition, cis dominantly explains (11.2%) regulatory differences, followed by trans (8.9%). Further analysis showed that water-limiting condition largely affects trans and cis + trans factors. On the molecular level, cis and/or trans regulatory divergence explains their genotypic differences and differential drought response. Between the two parental genotypes, Apo appears to exhibit more photosynthetic efficiency even under water-limiting condition and is ascribed to trans. Statistical analyses showed that regulatory divergence is significantly influenced by environmental conditions. Likewise, the mode of parental expression inheritance which drives heterosis (HET) is significantly affected by environmental conditions indicating the malleability of heterosis to external factors. Further analysis revealed that the HET class, dominance, was significantly enriched under water-stress condition. We also identified allelic imbalance switching in which several genes prefer IR64- (or Apo-) specific allele under non-stress condition but switched to Apo- (or IR64-) specific allele when exposed to water-stress condition.

Parental variation in CHG methylation is associated with allelic-specific expression in elite hybrid rice

Heterosis refers to the superior performance of hybrid lines over inbred parental lines. Besides genetic variation, epigenetic differences between parental lines are suggested to contribute to heterosis. However, the precise nature and extent of differences between the parental epigenomes and the reprograming in hybrids that govern heterotic gene expression remain unclear. In this work, we analyzed DNA methylomes and transcriptomes of the widely cultivated and genetically studied elite hybrid rice (Oryza sativa) SY63, the reciprocal hybrid, and the parental varieties ZS97 and MH63, for which high-quality reference genomic sequences are available. We showed that the parental varieties displayed substantial variation in genic methylation at CG and CHG (H = A, C, or T) sequences. Compared with their parents, the hybrids displayed dynamic methylation variation during development. However, many parental differentially methylated regions (DMRs) at CG and CHG sites were maintained in the hybrid. Only a small fraction of the DMRs displayed non-additive DNA methylation variation, which, however, showed no overall correlation relationship with gene expression variation. In contrast, most of the allelic-specific expression (ASE) genes in the hybrid were associated with DNA methylation, and the ASE negatively associated with allelic-specific methylation (ASM) at CHG. These results revealed a specific DNA methylation reprogramming pattern in the hybrid rice and pointed to a role for parental CHG methylation divergence in ASE, which is associated with phenotype variation and hybrid vigor in several plant species.

Molecular basis of heterosis and related breeding strategies reveal its importance in vegetable breeding

Heterosis has historically been exploited in plants; however, its underlying genetic mechanisms and molecular basis remain elusive. In recent years, due to advances in molecular biotechnology at the genome, transcriptome, proteome, and epigenome levels, the study of heterosis in vegetables has made significant progress. Here, we present an extensive literature review on the genetic and epigenetic regulation of heterosis in vegetables. We summarize six hypotheses to explain the mechanism by which genes regulate heterosis, improve upon a possible model of heterosis that is triggered by epigenetics, and analyze previous studies on quantitative trait locus effects and gene actions related to heterosis based on analyses of differential gene expression in vegetables. We also discuss the contributions of yield-related traits, including flower, fruit, and plant architecture traits, during heterosis development in vegetables (e.g., cabbage, cucumber, and tomato). More importantly, we propose a comprehensive breeding strategy based on heterosis studies in vegetables and crop plants. The description of the strategy details how to obtain F1 hybrids that exhibit heterosis based on heterosis prediction, how to obtain elite lines based on molecular biotechnology, and how to maintain heterosis by diploid seed breeding and the selection of hybrid simulation lines that are suitable for heterosis research and utilization in vegetables. Finally, we briefly provide suggestions and perspectives on the role of heterosis in the future of vegetable breeding.

An enhanced photosynthesis and carbohydrate metabolic capability contributes to heterosis of the cotton (Gossypium hirsutum) hybrid 'Huaza Mian H318', as revealed by genome-wide gene expression analysis

Background

Heterosis has been exploited for decades in different crops due to resulting in dramatic increases in yield, but relatively little molecular evidence on this topic was reported in cotton.

Results

The elite cotton hybrid variety ‘Huaza Mian H318’ (H318) and its parental lines were used to explore the source of its yield heterosis. A four-year investigation of yield-related traits showed that the boll number of H318 showed higher stability than that of its two parents, both in suitable and unsuitable climate years. In addition, the hybrid H318 grew faster and showed higher fresh and dry weights than its parental lines at the seedling stage. Transcriptome analysis of seedlings identified 17,308 differentially expressed genes (DEGs) between H318 and its parental lines, and 3490 extremely changed DEGs were screened out for later analysis. Most DEGs (3472/3490) were gathered between H318 and its paternal line (4–5), and only 64 DEGs were found between H318 and its maternal line (B0011), which implied that H318 displays more similar transcriptional patterns to its maternal parent at the seedling stage. GO and KEGG analyses showed that these DEGs were highly enriched in photosynthesis, lipid metabolic, carbohydrate metabolic and oxidation-reduction processes, and the expression level of these DEGs was significantly higher in H318 relative to its parental lines, which implied that photosynthesis, metabolism and stress resistances were enhanced in H318.

Conclusion

The enhanced photosynthesis, lipid and carbohydrate metabolic capabilities contribute to the heterosis of H318 at the seedling stage, and establishes a material foundation for subsequent higher boll-setting rates in complex field environments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07580-8.

Transcriptome-wide analysis of epitranscriptome and translational efficiency associated with heterosis in maize

Heterosis has been extensively utilized to increase productivity in crops, yet the underlying molecular mechanisms remain largely elusive. Here, we generated transcriptome-wide profiles of mRNA abundance, m6A methylation, and translational efficiency from the maize F1 hybrid B73×Mo17 and its two parental lines to ascertain the contribution of each regulatory layer to heterosis at the seedling stage. We documented that although the global abundance and distribution of m6A remained unchanged, a greater number of genes had gained an m6A modification in the hybrid. Superior variations were observed at the m6A modification and translational efficiency levels when compared with mRNA abundance between the hybrid and parents. In the hybrid, the vast majority of genes with m6A modification exhibited a non-additive expression pattern, the percentage of which was much higher than that at levels of mRNA abundance and translational efficiency. Non-additive genes involved in different biological processes were hierarchically coordinated by discrete combinations of three regulatory layers. These findings suggest that transcriptional and post-transcriptional regulation of gene expression make distinct contributions to heterosis in hybrid maize. Overall, this integrated multi-omics analysis provides a valuable portfolio for interpreting transcriptional and post-transcriptional regulation of gene expression in hybrid maize, and paves the way for exploring molecular mechanisms underlying hybrid vigor.

Improved redox homeostasis owing to the up-regulation of one-carbon metabolism and related pathways is crucial for yeast heterosis at high temperature

Heterosis or hybrid vigor is a common phenomenon in plants and animals; however, the molecular mechanisms underlying heterosis remain elusive, despite extensive studies on the phenomenon for more than a century. Here we constructed a large collection of F1 hybrids of Saccharomyces cerevisiae by spore-to-spore mating between homozygous wild strains of the species with different genetic distances and compared growth performance of the F1 hybrids with their parents. We found that heterosis was prevalent in the F1 hybrids at 40°C. A hump-shaped relationship between heterosis and parental genetic distance was observed. We then analyzed transcriptomes of selected heterotic and depressed F1 hybrids and their parents growing at 40°C and found that genes associated with one-carbon metabolism and related pathways were generally up-regulated in the heterotic F1 hybrids, leading to improved cellular redox homeostasis at high temperature. Consistently, genes related with DNA repair, stress responses, and ion homeostasis were generally down-regulated in the heterotic F1 hybrids. Furthermore, genes associated with protein quality control systems were also generally down-regulated in the heterotic F1 hybrids, suggesting a lower level of protein turnover and thus higher energy use efficiency in these strains. In contrast, the depressed F1 hybrids, which were limited in number and mostly shared a common aneuploid parental strain, showed a largely opposite gene expression pattern to the heterotic F1 hybrids. We provide new insights into molecular mechanisms underlying heterosis and thermotolerance of yeast and new clues for a better understanding of the molecular basis of heterosis in plants and 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.

In-Depth Observation on the Microbial and Fungal Community Structure of Four Contrasting Tomato Cultivation Systems in Soil Based and Soilless Culture Systems

As soil and soilless culture systems are highly dynamic environments, the structure of rhizosphere microbial communities is consistently adapting. There is a knowledge gap between the microbial community structure of soil based and soilless culture systems and thus we aimed at surveying their impact on diversity and composition of bacterial communities across a 10-month period in a tomato cultivation system. We compared community metrics between an soil based culture system fertilized with malt sprouts and blood meal, known for its slow and high mineralization rate, respectively and a soilless culture system fertilized with fish effluent or supplemented with an liquid organic fertilizer. Bacterial and fungal community composition was followed over time using two complementary techniques, phospholipid fatty acid analysis and 16S rRNA amplicon sequencing. Nitrogen dynamics and plant performance were assessed to provide insight on how bacterial diversity of soil and soilless microbial communities ultimately impacts productivity. Similar plant performance was observed in soilless culture systems and soil based system and yield was the highest with the aquaponics-derived fertilizer. Soil and soilless cultivating systems supplemented with different nitrogen-rich fertilizers differed on its characteristics throughout the experimental period. Fast-paced fluctuations in pH(H₂O) and nutrient cycling processes were observed in growing medium. Physicochemical characteristics changed over time and interacted with bacterial community metrics. Multivariate analysis showed that plant length, pH, Flavisolibacter, phosphorus, chloride, ammonium, potassium, calcium, magnesium, sodium, electrical conductivity, nitrate, sulfate, and the bacterial genera Desulfotomaculum, Solirubrobacter, Dehalococcoides, Bythopirellula, Steroidobacter, Litorilinea, Nonomuraea were the most significant factors discriminating between natural soils supplemented with animal and plant by-products. Long-term fertilizer regimes significantly changed the PLFA fingerprints in both the soilless culture and soil based culture system. The use of these by-products in the soil was positively associated with arbuscular mycorrhizal fungi (AMF), which may influence rhizosphere communities through root exudates and C translocation. Community structure was distinct and consistently different over time, despite the fertilizer supplementation. The fungal microbial community composition was less affected by pH, while the composition of the bacterial communities (Actinomycetes, Gram-negative bacteria, and Gram-positive bacteria) was closely defined by soil pH, demonstrating the significance of pH as driver of bacterial community composition. Fertilizer application may be responsible for variations over time in the ecosystem. Knowledge about the microbial interactions in tomato cultivating systems opens a window of opportunity for designing targeted fertilizers supporting sustainable crop production.

The comparative gene expression concern to the seed pigmentation in maize (Zea mays L.)

Maize seed pigmentation is one of the important issue to develop maize seed breeding. The differently gene expression was characterized and compared for three inbred lines, such as the pigment accumulated seed (CM22) and non-pigmented seed (CM5 and CM19) at 10 days after pollination. We obtained a total of 63,870, 82,496, and 54,555 contigs by de novo assembly to identify gene expression in the CM22, CM5, and CM19, respectably. In differentially expressed gene analysis, it was revealed that 7,044 genes were differentially expressed by at least two-fold, with 4,067 upregulated in colored maize inbred lines and 2,977 upregulated in colorless maize inbred lines. Of them,18 genes were included to the anthocyanin biosynthesis pathways, while 15 genes were upregulated in both CM22/5 and CM22/19. Additionally, 37 genes were detected in the metabolic pathway concern to the seed pigmentation by BINs analysis using MAPMAN software. Finally, these differently expressed genes may aid in the research on seed pigmentation in maize breeding programs.

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.

The Role of Transcriptional Regulation in Hybrid Vigor

The genetic basis of hybrid vigor in plants remains largely unsolved but strong evidence suggests that variation in transcriptional regulation can explain many aspects of this phenomenon. Natural variation in transcriptional regulation is highly abundant in virtually all species and thus a potential source of heterotic variability. Allele Specific Expression (ASE), which is tightly linked to parent of origin effects and modulated by complex interactions in cis and in trans, is generally considered to play a key role in explaining the differences between hybrids and parental lines. Here we discuss the recent developments in elucidating the role of transcriptional variation in a number of aspects of hybrid vigor, thereby bridging old paradigms and hypotheses with contemporary research in various species.

Genome-wide analysis of spatiotemporal allele-specific expression in F1 hybrids of meat- and egg-type chickens

In diploid organisms, each gene locus is composed of two parental alleles, which would interact with each other for determining the phenotypic variation. Better understanding of the allele-specific expression (ASE) in farm animals is much important to explore the genetic basis underlying economically important traits, which have been poorly understood yet. In this study, genome-wide analysis was applied to explore the spatiotemporal pattern of ASE in the F1 hybrids of chicken. First, meat- and egg-type chickens were selected for producing a full-sib F1 hybrid population (n = 57). Then, genome resequencing of two parents and 38 offspring were performed and liver and breast muscle samples (n = 38) were subjected to strand-specific RNA sequencing (ssRNA-seq) for ASE detection at 1, 28, and 56 days of age, respectively. The results accurately identified a total of 465 informative genes that could be distinguished with respect to their parental origins. There were 0.4% - 4.1% of informative genes showing ASE, and 57 of them were found across different tissues and time points. Besides, most ASE genes in chickens were tissue-specific, and no matter what the time-point pattern of one ASE gene, the same parental allele of this gene almost showed consistently higher or lower expression across all time points in the same type tissue. In conclusion, this study indicated that most of ASE genes were tissue-specific and time-dependent.

Variant phasing and haplotypic expression from long-read sequencing in maize

Haplotype phasing maize genetic variants is important for genome interpretation, population genetic analysis and functional analysis of allelic activity. We performed an isoform-level phasing study using two maize inbred lines and their reciprocal crosses, based on single-molecule, full-length cDNA sequencing. To phase and analyze transcripts between hybrids and parents, we developed IsoPhase. Using this tool, we validated the majority of SNPs called against matching short-read data from embryo, endosperm and root tissues, and identified allele-specific, gene-level and isoform-level differential expression between the inbred parental lines and hybrid offspring. After phasing 6907 genes in the reciprocal hybrids, we annotated the SNPs and identified large-effect genes. In addition, we identified parent-of-origin isoforms, distinct novel isoforms in maize parent and hybrid lines, and imprinted genes from different tissues. Finally, we characterized variation in cis- and trans-regulatory effects. Our study provides measures of haplotypic expression that could increase accuracy in studies of allelic expression.

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.

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

Background

The cellular machinery for cell wall synthesis and metabolism is encoded by members of large multi-gene families. Maize is both a genetic model for grass species and a potential source of lignocellulosic biomass from crop residues. Genetic improvement of maize for its utility as a bioenergy feedstock depends on identification of the specific gene family members expressed during secondary wall development in stems.

Results

High-throughput sequencing of transcripts expressed in developing rind tissues of stem internodes provided a comprehensive inventory of cell wall-related genes in maize (Zea mays, cultivar B73). Of 1239 of these genes, 854 were expressed among the internodes at ≥95 reads per 20 M, and 693 of them at ≥500 reads per 20 M. Grasses have cell wall compositions distinct from non-commelinid species; only one-quarter of maize cell wall-related genes expressed in stems were putatively orthologous with those of the eudicot Arabidopsis. Using a slope-metric algorithm, five distinct patterns for sub-sets of co-expressed genes were defined across a time course of stem development. For the subset of genes associated with secondary wall formation, fifteen sequence motifs were found in promoter regions. The same members of gene families were often expressed in two maize inbreds, B73 and Mo17, but levels of gene expression between them varied, with 30% of all genes exhibiting at least a 5-fold difference at any stage. Although presence-absence and copy-number variation might account for much of these differences, fold-changes of expression of a CADa and a FLA11 gene were attributed to polymorphisms in promoter response elements.

Conclusions

Large genetic variation in maize as a species precludes the extrapolation of cell wall-related gene expression networks even from one common inbred line to another. Elucidation of genotype-specific expression patterns and their regulatory controls will be needed for association panels of inbreds and landraces to fully exploit genetic variation in maize and other bioenergy grass species.

Tomato plants rather than fertilizers drive microbial community structure in horticultural growing media

Synthetic fertilizer production is associated with a high environmental footprint, as compounds typically dissolve rapidly leaching emissions to the atmosphere or surface waters. We tested two recovered nutrients with slower release patterns, as promising alternatives for synthetic fertilizers: struvite and a commercially available organic fertilizer. Using these fertilizers as nitrogen source, we conducted a rhizotron experiment to test their effect on plant performance and nutrient recovery in juvenile tomato plants. Plant performance was significantly improved when organic fertilizer was provided, promoting higher shoot biomass. Since the microbial community influences plant nitrogen availability, we characterized the root-associated microbial community structure and functionality. Analyses revealed distinct root microbial community structure when different fertilizers were supplied. However, plant presence significantly increased the similarity of the microbial community over time, regardless of fertilization. Additionally, the presence of the plant significantly reduced the potential ammonia oxidation rates, implying a possible role of the rhizosheath microbiome or nitrification inhibition by the plant. Our results indicate that nitrifying community members are impacted by the type of fertilizer used, while tomato plants influenced the potential ammonia-oxidizing activity of nitrogen-related rhizospheric microbial communities. These novel insights on interactions between recovered fertilizers, plant and associated microbes can contribute to develop sustainable crop production systems.

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.