Patterns of genome-wide allele-specific expression in hybrid rice and the implications on the genetic basis of heterosis

Seed development-related genes contribute to high yield heterosis in integrated utilization of elite autotetraploid and neo-tetraploid rice

Introduction

Autotetraploid rice holds high resistance to abiotic stress and substantial promise for yield increase, but it could not be commercially used because of low fertility. Thus, our team developed neo-tetraploid rice with high fertility and hybrid vigor when crossed with indica autotetraploid rice. Despite these advances, the molecular mechanisms underlying this heterosis remain poorly understood.

Methods

An elite indica autotetraploid rice line (HD11) was used to cross with neo-tetraploid rice, and 34 hybrids were obtained to evaluate agronomic traits related to yield. WE-CLSM, RNA-seq, and CRISPR/Cas9 were employed to observe endosperm structure and identify candidate genes from two represent hybrids.

Results and discussion

These hybrids showed high seed setting and an approximately 55% increase in 1000-grain weight, some of which achieved grain yields comparable to those of the diploid rice variety. The endosperm observations indicated that the starch grains in the hybrids were more compact than those in paternal lines. A total of 119 seed heterosis related genes (SHRGs) with different expressions were identified, which might contribute to high 1000-grain weight heterosis in neo-tetraploid hybrids. Among them, 12 genes had been found to regulate grain weight formation, including OsFl3, ONAC023, OsNAC024, ONAC025, ONAC026, RAG2, FLO4, FLO11, OsISA1, OsNF-YB1, NF-YC12, and OsYUC9. Haplotype analyses of these 12 genes revealed the various effects on grain weight among different haplotypes. The hybrids could polymerize more dominant haplotypes of above grain weight regulators than any homozygous cultivar. Moreover, two SHRGs (OsFl3 and SHRG2) mutants displayed a significant reduction in 1000-grain weight and an increase in grain chalkiness, indicating that OsFl3 and SHRG2 positively regulate grain weight. Our research has identified a valuable indica autotetraploid germplasm for generating strong yield heterosis in combination with neo-tetraploid lines and gaining molecular insights into the regulatory processes of heterosis in tetraploid rice.

A haplotype-resolved genome assembly of Malus domestica 'Red Fuji'

The 'Red Fuji' apple (Malus domestica), is one of the most important and popular economic crops worldwide in the fruit industry. Using PacBio HiFi long reads and Hi-C reads, we assembled a high-quality haplotype-resolved genome of 'Red Fuji', with sizes of 668.7 and 668.8 Mb, and N50 sizes of 34.1 and 31.4 Mb. About 97.2% of sequences were anchored in 34 chromosomes. We annotated both haploid genomes, identifying a total of 95,439 protein-coding genes in the two haplotype genomes, with 98% functional annotation. The haplotype-resolved genome of 'Red Fuji' apple stands as a precise benchmark for an array of analyses, such as comparative genomics, transcriptomics, and allelic expression studies. This comprehensive resource is paramount in unraveling variations in allelic expression, advancing quality improvements, and refining breeding efforts.

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.

Natural variations of heterosis-related allele-specific expression genes in promoter regions lead to allele-specific expression in maize

BACKGROUND

Heterosis has successfully enhanced maize productivity and quality. Although significant progress has been made in delineating the genetic basis of heterosis, the molecular mechanisms underlying its genetic components remain less explored. Allele-specific expression (ASE), the imbalanced expression between two parental alleles in hybrids, is increasingly being recognized as a factor contributing to heterosis. ASE is a complex process regulated by both epigenetic and genetic variations in response to developmental and environmental conditions.

RESULTS

In this study, we explored the differential characteristics of ASE by analyzing the transcriptome data of two maize hybrids and their parents under four light conditions. On the basis of allele expression patterns in different hybrids under various conditions, ASE genes were divided into three categories: bias-consistent genes involved in basal metabolic processes in a functionally complementary manner, bias-reversal genes adapting to the light environment, and bias-specific genes maintaining cell homeostasis. We observed that 758 ASE genes (ASEGs) were significantly overlapped with heterosis quantitative trait loci (QTLs), and high-frequency variations in the promoter regions of heterosis-related ASEGs were identified between parents. In addition, 10 heterosis-related ASEGs participating in yield heterosis were selected during domestication.

CONCLUSIONS

The comprehensive analysis of ASEGs offers a distinctive perspective on how light quality influences gene expression patterns and gene-environment interactions, with implications for the identification of heterosis-related ASEGs to enhance maize yield.

Evolutionary divergence of subgenomes in common carp provides insights into speciation and allopolyploid success

Hybridization and polyploidization have made great contributions to speciation, heterosis, and agricultural production within plants, but there is still limited understanding and utilization in animals. Subgenome structure and expression reorganization and cooperation post hybridization and polyploidization are essential for speciation and allopolyploid success. However, the mechanisms have not yet been comprehensively assessed in animals. Here, we produced a high-fidelity reference genome sequence for common carp, a typical allotetraploid fish species cultured worldwide. This genome enabled in-depth analysis of the evolution of subgenome architecture and expression responses. Most genes were expressed with subgenome biases, with a trend of transition from the expression of subgenome A during the early stages to that of subgenome B during the late stages of embryonic development. While subgenome A evolved more rapidly, subgenome B contributed to a greater level of expression during development and under stressful conditions. Stable dominant patterns for homoeologous gene pairs both during development and under thermal stress suggest a potential fixed heterosis in the allotetraploid genome. Preferentially expressing either copy of a homoeologous gene at higher levels to confer development and response to stress indicates the dominant effect of heterosis. The plasticity of subgenomes and their shifting of dominant expression during early development, and in response to stressful conditions, provide novel insights into the molecular basis of the successful speciation, evolution, and heterosis of the allotetraploid common carp.

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.

The OsWRKY72-OsAAT30/OsGSTU26 module mediates reactive oxygen species scavenging to drive heterosis for salt tolerance in hybrid rice

Hybrid rice (Oryza sativa) generally outperforms its inbred parents in yield and stress tolerance, a phenomenon termed heterosis, but the underlying mechanism is not completely understood. Here, we combined transcriptome, proteome, physiological, and heterosis analyses to examine the salt response of super hybrid rice Chaoyou1000 (CY1000). In addition to surpassing the mean values for its two parents (mid-parent heterosis), CY1000 exhibited a higher reactive oxygen species scavenging ability than both its parents (over-parent heterosis or heterobeltiosis). Nonadditive expression and allele-specific gene expression assays showed that the glutathione S-transferase gene OsGSTU26 and the amino acid transporter gene OsAAT30 may have major roles in heterosis for salt tolerance, acting in an overdominant fashion in CY1000. Furthermore, we identified OsWRKY72 as a common transcription factor that binds and regulates OsGSTU26 and OsAAT30. The salt-sensitive phenotypes were associated with the OsWRKY72paternal genotype or the OsAAT30maternal genotype in core rice germplasm varieties. OsWRKY72paternal specifically repressed the expression of OsGSTU26 under salt stress, leading to salinity sensitivity, while OsWRKY72maternal specifically repressed OsAAT30, resulting in salinity tolerance. These results suggest that the OsWRKY72-OsAAT30/OsGSTU26 module may play an important role in heterosis for salt tolerance in an overdominant fashion in CY1000 hybrid rice, providing valuable clues to elucidate the mechanism of heterosis for salinity tolerance in hybrid rice.

Near-gapless genome and transcriptome analyses provide insights into fruiting body development in Lentinula edodes

Fruiting body development in macrofungi is an intensive research subject. In this study, high-quality genomes were assembled for two sexually compatible monokaryons from a heterokaryotic Lentinula edodes strain WX1, and variations in L. edodes genomes were analyzed. Specifically, differential gene expression and allele-specific expression (ASE) were analyzed using the two monokaryotic genomes and transcriptome data from four different stages of fruiting body development in WX1. Results revealed that after aeration, mycelia sensed cell wall stress, pheromones, and a decrease in CO2 concentration, leading to up-regulated expression in genes related to cell adhesion, cell wall remodeling, proteolysis, and lipid metabolism, which may promote primordium differentiation. Aquaporin genes and those related to proteolysis, mitosis, lipid, and carbohydrate metabolism may play important roles in primordium development, while genes related to tissue differentiation and sexual reproduction were active in fruiting body. Several essential genes for fruiting body development were allele-specifically expressed and the two nuclear types could synergistically regulate fruiting body development by dominantly expressing genes with different functions. ASE was probably induced by long terminal repeat-retrotransposons. Findings here contribute to the further understanding of the mechanism of fruiting body development in macrofungi.

Haplotype-resolved genome assembly of Phanera championii reveals molecular mechanisms of flavonoid synthesis and adaptive evolution

Phanera championii is a medicinal liana plant that has successfully adapted to hostile karst habitats. Despite extensive research on its medicinal components and pharmacological effects, the molecular mechanisms underlying the biosynthesis of critical flavonoids and its adaptation to karst habitats remain elusive. In this study, we performed high-coverage PacBio and Hi-C sequencing of P. championii, which revealed its high heterozygosity and phased the genome into two haplotypes: Hap1 (384.60 Mb) and Hap2 (383.70 Mb), encompassing a total of 58 612 annotated genes. Comparative genomes analysis revealed that P. championii experienced two whole-genome duplications (WGDs), with approximately 59.59% of genes originating from WGD events, thereby providing a valuable genetic resource for P. championii. Moreover, we identified a total of 112 genes that were strongly positively selected. Additionally, about 81.60 Mb of structural variations between the two haplotypes. The allele-specific expression patterns suggested that the dominant effect of P. championii was the elimination of deleterious mutations and the promotion of beneficial mutations to enhance fitness. Moreover, our transcriptome and metabolome analysis revealed alleles in different tissues or different haplotypes collectively regulate the synthesis of flavonoid metabolites. In summary, our comprehensive study highlights the significance of genomic and morphological adaptation in the successful adaptation of P. championii to karst habitats. The high-quality phased genomes obtained in this study serve as invaluable genomic resources for various applications, including germplasm conservation, breeding, evolutionary studies, and elucidation of pathways governing key biological traits of P. championii.

A haplotype-resolved genome provides insight into allele-specific expression in wild walnut (Juglans regia L.)

Wild germplasm resources are crucial for gene mining and molecular breeding because of their special trait performance. Haplotype-resolved genome is an ideal solution for fully understanding the biology of subgenomes in highly heterozygous species. Here, we surveyed the genome of a wild walnut tree from Gongliu County, Xinjiang, China, and generated a haplotype-resolved reference genome of 562.99 Mb (contig N50 = 34.10 Mb) for one haplotype (hap1) and 561.07 Mb (contig N50 = 33.91 Mb) for another haplotype (hap2) using PacBio high-fidelity (HiFi) reads and Hi-C technology. Approximately 527.20 Mb (93.64%) of hap1 and 526.40 Mb (93.82%) of hap2 were assigned to 16 pseudochromosomes. A total of 41039 and 39744 protein-coding gene models were predicted for hap1 and hap2, respectively. Moreover, 123 structural variations (SVs) were identified between the two haplotype genomes. Allele-specific expression genes (ASEGs) that respond to cold stress were ultimately identified. These datasets can be used to study subgenome evolution, for functional elite gene mining and to discover the transcriptional basis of specific traits related to environmental adaptation in wild walnut.

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.

Haplotype-phased genome unveils the butylphthalide biosynthesis and homoploid hybrid origin of Ligusticum chuanxiong

Butylphthalide is one of the first-line drugs for ischemic stroke therapy, while no biosynthetic enzyme for butylphthalide has been reported. Here, we present a haplotype-resolved genome of Ligusticum chuanxiong, a long-cultivated and phthalide-rich medicinal plant in Apiaceae. On the basis of comprehensive screening, four Fe(II)- and 2-oxoglutarate-dependent dioxygenases and two CYPs were mined and further biochemically verified as phthalide C-4/C-5 desaturases (P4,5Ds) that effectively promoted the forming of (S)-3-n-butylphthalide and butylidenephthalide. The substrate promiscuity and functional redundancy featured for P4,5Ds may contribute to the high phthalide diversity in L. chuanxiong. Notably, comparative genomic evidence supported L. chuanxiong as a homoploid hybrid with Ligusticum sinense as a potential parent. The two haplotypes demonstrated exceptional structure variance and diverged around 3.42 million years ago. Our study is an icebreaker for the dissection of phthalide biosynthetic pathway and reveals the hybrid origin of L. chuanxiong, which will facilitate the metabolic engineering for (S)-3-n-butylphthalide production and breeding for L. chuanxiong.

Allele-dependent expression and functionality of lipid enzyme phospholipid:diacylglycerol acyltransferase affect diatom carbon storage and growth

In the acyl-CoA-independent pathway of triacylglycerol (TAG) synthesis unique to plants, fungi, and algae, TAG formation is catalyzed by the enzyme phospholipid:diacylglycerol acyltransferase (PDAT). The unique PDAT gene of the model diatom Phaeodactylum tricornutum strain CCMP2561 boasts 47 single nucleotide variants within protein coding regions of the alleles. To deepen our understanding of TAG synthesis, we observed the allele-specific expression of PDAT by the analysis of 87 published RNA-sequencing (RNA-seq) data and experimental validation. The transcription of one of the two PDAT alleles, Allele 2, could be specifically induced by decreasing nitrogen concentrations. Overexpression of Allele 2 in P. tricornutum substantially enhanced the accumulation of TAG by 44% to 74% under nutrient stress; however, overexpression of Allele 1 resulted in little increase of TAG accumulation. Interestingly, a more serious growth inhibition was observed in the PDAT Allele 1 overexpression strains compared with Allele 2 counterparts. Heterologous expression in yeast (Saccharomyces cerevisiae) showed that enzymes encoded by PDAT Allele 2 but not Allele 1 had TAG biosynthetic activity, and 7 N-terminal and 3 C-terminal amino acid variants between the 2 allele-encoded proteins substantially affected enzymatic activity. P. tricornutum PDAT, localized in the innermost chloroplast membrane, used monogalactosyldiacylglycerol and phosphatidylcholine as acyl donors as demonstrated by the increase of the 2 lipids in PDAT knockout lines, which indicated a common origin in evolution with green algal PDATs. Our study reveals unequal roles among allele-encoded PDATs in mediating carbon storage and growth in response to nitrogen stress and suggests an unsuspected strategy toward lipid and biomass improvement for biotechnological purposes.

Non-Additive and Asymmetric Allelic Expression of p38 mapk in Hybrid Tilapia (Oreochromis niloticus ♀ × O. aureus ♂)

Hybridization is a widely used breeding technique in fish species that enhances desirable traits in cultured species through heterosis. However, the mechanism by which hybrids alter gene expression to form heterosis remains unclear. In this study, a group of hybrid tilapia was used to elucidate heterosis through interspecies crossing. Specifically, p38 was analyzed to describe the regulation of gene expression variation in hybrid tilapia. Transcripts from the Nile tilapia allele were found to be significantly higher than those from the blue tilapia allele in hybrid individuals, indicating that the expression of p38 was dominated by Nile tilapia sub-genomic alleles. The study also found a compensatory interaction of cis- and trans-acting elements of the Nile tilapia and blue tilapia sub-genomes, inducing a non-additive expression of p38 in hybrids. Eight specific SNPs were identified in the p38 promoter regions of Nile tilapia and blue tilapia, and were found to be promoter differences leading to differences in gene expression efficiencies between parental alleles using a dual-luciferase reporter system. This study provides insights into the non-additive expression patterns of key functional genes in fish hybrids related to growth and immunity response.

Heterosis in crop improvement

Heterosis, also known as hybrid vigor, is the phenomenon wherein a progeny exhibits superior traits relative to one or both parents. In terms of crop breeding, this usually refers to the yield advantage of F1 hybrids over both inbred parents. The development of high-yielding hybrid cultivars across a wider range of crops is key to meeting future food demands. However, conventional hybrid breeding strategies are proving to be exceptionally challenging to apply commercially in many self-pollinating crops, particularly wheat and barley. Currently in these crops, the relative performance advantage of hybrids over inbred line cultivars does not outweigh the cost of hybrid seed production. Here, we review the genetic basis of heterosis, discuss the challenges in hybrid breeding, and propose a strategy to recruit multiple heterosis-associated genes to develop lines with improved agronomic characteristics. This strategy leverages modern genetic engineering tools to synthesize supergenes by fusing multiple heterotic alleles across multiple heterosis-associated loci. We outline a plan to assess the feasibility of this approach to improve line performance using barley (Hordeum vulgare) as the model self-pollinating crop species, and a few heterosis-associated genes. The proposed method can be applied to all crops for which heterotic gene combinations can be identified.

Nonadditive and allele-specific expression of ghrelin in hybrid tilapia

Background

Interspecies hybridization is an important breeding method to generate fishes with heterosis in aquaculture. Using this method, hybrid Nile tilapia (Oreochromis niloticus, ♀) × blue tilapia (Oreochromis aureus, ♂) has been produced and widely farmed due to its growth and appetite superiorities. However, the genetic mechanism of these advanced traits is still not well understood. Ghrelin is a crucial gene that regulates growth and appetite in fishes. In the present study, we focused on the expression characteristics and its regulation of ghrelin in the hybrid.

Results

The tissue distribution analysis showed that ghrelin was predominantly expressed in the stomach in the hybrid. Ghrelin was more highly expressed in the stomach in the hybrid and Nile tilapia, compared to blue tilapia, showing a nonadditive pattern. Two single-nucleotide polymorphism (SNP) sites were identified including T/C and C/G from the second exon in the ghrelin gene from Nile tilapia and blue tilapia. By pyrosequencing based on the SNP sites, the allele-specific expression (ASE) of ghrelin in the hybrid was assayed. The result indicated that ghrelin in the hybrid showed higher maternal allelic transcript ratios. Fasting significantly increased ghrelin overall expression at 4, 8, 12, 24, and 48 h. In addition, higher maternal allelic transcript ratios were not changed in the fasting hybrids at 48 h. The cis and trans effects were determined by evaluating the overall expression and ASE values in the hybrid. The expression of ghrelin was mediated by compensating cis and trans effects in hybrid.

Conclusion

In summary, the present lines of evidence showed the nonadditive expression of ghrelin in the hybrid tilapia and its regulation by subgenomes, offering new insight into gene expression characteristics in hybrids.

Transcriptome Shock in Developing Embryos of a Brassica napus and Brassica rapa Hybrid

Interspecific crosses that fuse the genomes of two different species may result in overall gene expression changes in the hybrid progeny, called 'transcriptome shock'. To better understand the expression pattern after genome merging during the early stages of allopolyploid formation, we performed RNA sequencing analysis on developing embryos of Brassica rapa, B. napus, and their synthesized allotriploid hybrids. Here, we show that the transcriptome shock occurs in the developing seeds of the hybrids. Of the homoeologous gene pairs, 17.1% exhibit expression bias, with an overall expression bias toward B. rapa. The expression level dominance also biases toward B. rapa, mainly induced by the expression change in homoeologous genes from B. napus. Functional enrichment analysis revealed significant differences in differentially expressed genes (DEGs) related to photosynthesis, hormone synthesis, and other pathways. Further study showed that significant changes in the expression levels of the key transcription factors (TFs) could regulate the overall interaction network in the developing embryo, which might be an essential cause of phenotype change. In conclusion, the present results have revealed the global changes in gene expression patterns in developing seeds of the hybrid between B. rapa and B. napus, and provided novel insights into the occurrence of transcriptome shock for harnessing heterosis.

Differential subgenome expression underlies biomass accumulation in allotetraploid Pennisetum giganteum

BACKGROUND

Pennisetum giganteum (AABB, 2n = 4x = 28) is a C4 plant in the genus Pennisetum with origin in Africa but currently also grown in Asia and America. It is a crucial forage and potential energy grass with significant advantages in yield, stress resistance, and environmental adaptation. However, the mechanisms underlying these advantageous traits remain largely unexplored. Here, we present a high-quality genome assembly of the allotetraploid P. giganteum aiming at providing insights into biomass accumulation.

RESULTS

Our assembly has a genome size 2.03 Gb and contig N50 of 88.47 Mb that was further divided into A and B subgenomes. Genome evolution analysis revealed the evolutionary relationships across the Panicoideae subfamily lineages and identified numerous genome rearrangements that had occurred in P. giganteum. Comparative genomic analysis showed functional differentiation between the subgenomes. Transcriptome analysis found no subgenome dominance at the overall gene expression level; however, differentially expressed homoeologous genes and homoeolog-specific expressed genes between the two subgenomes were identified, suggesting that complementary effects between the A and B subgenomes contributed to biomass accumulation of P. giganteum. Besides, C4 photosynthesis-related genes were significantly expanded in P. giganteum and their sequences and expression patterns were highly conserved between the two subgenomes, implying that both subgenomes contributed greatly and almost equally to the highly efficient C4 photosynthesis in P. giganteum. We also identified key candidate genes in the C4 photosynthesis pathway that showed sustained high expression across all developmental stages of P. giganteum.

CONCLUSIONS

Our study provides important genomic resources for elucidating the genetic basis of advantageous traits in polyploid species, and facilitates further functional genomics research and genetic improvement of P. giganteum.

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.

Haplotype-Resolution Transcriptome Analysis Reveals Important Responsive Gene Modules and Allele-Specific Expression Contributions under Continuous Salt and Drought in Camellia sinensis

The tea plant, Camellia sinensis (L.) O. Kuntze, is one of the most important beverage crops with significant economic and cultural value. Global climate change and population growth have led to increased salt and drought stress, negatively affecting tea yield and quality. The response mechanism of tea plants to these stresses remains poorly understood due to the lack of reference genome-based transcriptional descriptions. This study presents a high-quality genome-based transcriptome dynamic analysis of C. sinensis' response to salt and drought stress. A total of 2244 upregulated and 2164 downregulated genes were identified under salt and drought stress compared to the control sample. Most of the differentially expression genes (DEGs) were found to involve divergent regulation processes at different time points under stress. Some shared up- and downregulated DEGs related to secondary metabolic and photosynthetic processes, respectively. Weighted gene co-expression network analysis (WGCNA) revealed six co-expression modules significantly positively correlated with C. sinensis' response to salt or drought stress. The MEpurple module indicated crosstalk between the two stresses related to ubiquitination and the phenylpropanoid metabolic regulation process. We identified 1969 salt-responsive and 1887 drought-responsive allele-specific expression (ASE) genes in C. sinensis. Further comparison between these ASE genes and tea plant heterosis-related genes suggests that heterosis likely contributes to the adversity and stress resistance of C. sinensis. This work offers new insight into the underlying mechanisms of C. sinensis' response to salt and drought stress and supports the improved breeding of tea plants with enhanced salt and drought tolerance.

Transcriptomic and methylomic analyses provide insights into the molecular mechanism and prediction of heterosis in rice

Heterosis has been widely used in multiple crops. However, the molecular mechanism and prediction of heterosis remains elusive. We generated five F1 hybrids [four showing better-parent heterosis (BPH) and one showing mid-parent heterosis], and performed the transcriptomic and methylomic analyses to identify the candidate genes for BPH and explore the molecular mechanism of heterosis and the potential predictors for heterosis. Transcriptomic results showed that most of the differentially expressed genes shared in the four better-parent hybrids were significantly enriched into the terms of molecular function, and the additive and dominant effects played crucial roles for BPH. DNA methylation level, especially in CG context, significantly and positively correlated with grain yield per plant. The ratios of differentially methylated regions in CG context in exons to transcription start sites between the parents exhibited significantly negative correlation with the heterosis levels of their hybrids, as was further confirmed in 24 pairwise comparisons of other rice lines, implying that this ratio could be a feasible predictor for heterosis level, and this ratio of less than 5 between parents in early growth stages might be a critical index for judging that their F1 hybrids would show BPH. Additionally, we identified some important genes showing differential expression and methylation, such as OsDCL2, Pi5, DTH2, DTH8, Hd1 and GLW7 in the four better-parent hybrids as the candidate genes for BPH. Our findings helped shed more light on the molecular mechanism and heterosis prediction.

Design, execution, and interpretation of plant RNA-seq analyses

Genomics has transformed our understanding of the genetic architecture of traits and the genetic variation present in plants. Here, we present a review of how RNA-seq can be performed to tackle research challenges addressed by plant sciences. We discuss the importance of experimental design in RNA-seq, including considerations for sampling and replication, to avoid pitfalls and wasted resources. Approaches for processing RNA-seq data include quality control and counting features, and we describe common approaches and variations. Though differential gene expression analysis is the most common analysis of RNA-seq data, we review multiple methods for assessing gene expression, including detecting allele-specific gene expression and building co-expression networks. With the production of more RNA-seq data, strategies for integrating these data into genetic mapping pipelines is of increased interest. Finally, special considerations for RNA-seq analysis and interpretation in plants are needed, due to the high genome complexity common across plants. By incorporating informed decisions throughout an RNA-seq experiment, we can increase the knowledge gained.

Genetic Analysis for the Flag Leaf Heterosis of a Super-Hybrid Rice WFYT025 Combination Using RNA-Seq

The photosynthetic capacity of flag leaf plays a key role in grain yield in rice. Nevertheless, there are few studies on the heterosis of the rice flag leaf. Therefore, this study focuses on investigating the genetic basis of heterosis for flag leaf in the indica super hybrid rice combination WFYT025 in China using a high-throughput next-generation RNA-seq strategy. We analyzed the gene expression of flag leaf in different environments and different time periods between WFYT025 and its female parent. After obtaining the gene expression profile of the flag leaf, we further investigated the gene regulatory network. Weighted gene expression network analysis (WGCNA) was used to identify the co-expressed gene sets, and a total of 5000 highly expressed genes were divided into 24 co-expression groups. In CHT025, we found 13 WRKY family transcription factors in SDG_hps under the environment of early rice and 16 WRKY family genes in SDG_hps of under the environment of middle rice. We found nine identical transcription factors in the two stages. Except for five reported TFs, the other four TFs might play an important role in heterosis for grain number and photosynthesis. Transcription factors such as WRKY3, WRKY68, and WRKY77 were found in both environments. To eliminate the influence of the environment, we examined the metabolic pathway with the same SDG_hp (SSDG_hp) in two environments. There were 312 SSDG_hps in total. These SSDG_hps mainly focused on the phosphorus metallic process, phosphorylation, plasma membrane, etc. These results provide resources for studying heterosis during super hybrid rice flag leaf development.

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.

Dissecting the genetic basis of heterosis in elite super-hybrid rice

Y900 is one of the top hybrid rice (Oryza sativa) varieties, with its yield exceeding 15 t·hm-2. To dissect the mechanism of heterosis, we sequenced the male parent line R900 and female parent line Y58S using long-read and Hi-C technology. High-quality reference genomes of 396.41 Mb and 398.24 Mb were obtained for R900 and Y58S, respectively. Genome-wide variations between the parents were systematically identified, including 1,367,758 single-nucleotide polymorphisms, 299,149 insertions/deletions, and 4,757 structural variations. The level of variation between Y58S and R900 was the lowest among the comparisons of Y58S with other rice genomes. More than 75% of genes exhibited variation between the two parents. Compared with other two-line hybrids sharing the same female parent, the portion of Geng/japonica (GJ)-type genetic components from different male parents increased with yield increasing in their corresponding hybrids. Transcriptome analysis revealed that the partial dominance effect was the main genetic effect that constituted the heterosis of Y900. In the hybrid, both alleles from the two parents were expressed, and their expression patterns were dynamically regulated in different tissues. The cis-regulation was dominant for young panicle tissues, while trans-regulation was more common in leaf tissues. Overdominance was surprisingly prevalent in stems and more likely regulated by the trans-regulation mechanism. Additionally, R900 contained many excellent GJ haplotypes, such as NARROW LEAF1, Oryza sativa SQUAMOSA PROMOTER BINDING PROTEIN-LIKE13, and Grain number, plant height, and heading date8, making it a good complement to Y58S. The fine-tuned mechanism of heterosis involves genome-wide variation, GJ introgression, key functional genes, and dynamic gene/allele expression and regulation pattern changes in different tissues and growth stages.

Haplotype mapping of H3K27me3-associated chromatin interactions defines topological regulation of gene silencing in rice

Histone modification H3K27me3 is an important chromatin mark that plays vital roles in repressing expression of developmental genes. Here, we construct high-resolution 3D genome maps using long-read chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) and characterize H3K27me3-associated chromatin interactions in an elite rice hybrid, Shanyou 63. We find that many H3K27me3-marked regions may function as silencer-like regulatory elements. The silencer-like elements can come into proximity with distal target genes via forming chromatin loops in 3D space of the nuclei, regulating gene silencing and plant traits. Natural and induced deletion of silencers upregulate expression of distal connected genes. Furthermore, we identify extensive allele-specific chromatin loops. We find that genetic variations alter allelic chromatin topology, thus modulating allelic gene imprinting in rice hybrids. In conclusion, the characterization of silencer-like regulatory elements and haplotype-resolved chromatin interaction maps provide insights into the understanding of molecular mechanisms underlying allelic gene silencing and plant trait controlling.

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.

An enhanced network of energy metabolism, lysine acetylation, and growth-promoting protein accumulation is associated with heterosis in elite hybrid rice

Heterosis refers to the superior performance of a hybrid compared with its parental lines. Although several genetic and molecular models have been proposed to explain heterosis, it remains unclear how hybrid cells integrate complementary gene expression or activity to drive heterotic growth. In this work, we show that accumulation of growth-promoting and energy metabolism proteins, enhanced energy metabolism activities, and increased protein lysine acetylation were associated with superior growth of the panicle meristem in the elite hybrid rice Shanyou 63 relative to its parental varieties. Metabolism of nuclear/cytosolic acetyl-coenzyme A was also enhanced in the hybrid, which paralleled increases in histone H3 acetylation to selectively target the expression of growth-promoting and metabolic genes. Lysine acetylation of cellular proteins, including TARGET OF RAPAMYCIN complex 1, ribosomal proteins, and energy metabolism enzymes, was also augmented and/or remodeled to modulate their activities. The data indicate that an enhanced network of energy-producing metabolic activity and growth-promoting histone acetylation/gene expression in the hybrid could contribute to its superior growth rate and may constitute a model to explain 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.

Sequencing and Assembly of Polyploid Genomes

Polyploidy has been observed throughout major eukaryotic clades and has played a vital role in the evolution of angiosperms. Recent polyploidizations often result in highly complex genome structures, posing challenges to genome assembly and phasing. Recent advances in sequencing technologies and genome assembly algorithms have enabled high-quality, near-complete chromosome-level assemblies of polyploid genomes. Advances in novel sequencing technologies include highly accurate single-molecule sequencing with HiFi reads, chromosome conformation capture with Hi-C technique, and linked reads sequencing. Additionally, new computational approaches have also significantly improved the precision and reliability of polyploid genome assembly and phasing, such as HiCanu, hifiasm, ALLHiC, and PolyGembler. Herein, we review recently published polyploid genomes and compare the various sequencing, assembly, and phasing approaches that are utilized in these genome studies. Finally, we anticipate that accurate and telomere-to-telomere chromosome-level assembly of polyploid genomes could ultimately become a routine procedure in the near future.

Identification of regulatory variants of carboxylesterase 1 (CES1): A proof-of-concept study for the application of the Allele-Specific Protein Expression (ASPE) assay in identifying cis-acting regulatory genetic polymorphisms

It is challenging to study regulatory genetic variants as gene expression is affected by both genetic polymorphisms and non-genetic regulators. The mRNA allele-specific expression (ASE) assay has been increasingly used for the study of cis-acting regulatory variants because cis-acting variants affect gene expression in an allele-specific manner. However, poor correlations between mRNA and protein expressions were observed for many genes, highlighting the importance of studying gene expression regulation at the protein level. In the present study, we conducted a proof-of-concept study to utilize a recently developed allele-specific protein expression (ASPE) assay to identify the cis-acting regulatory variants of CES1 using a large set of human liver samples. The CES1 gene encodes for carboxylesterase 1 (CES1), the most abundant hepatic hydrolase in humans. Two cis-acting regulatory variants were found to be significantly associated with CES1 ASPE, CES1 protein expression, and its catalytic activity on enalapril hydrolysis in human livers. Compared to conventional gene expression-based approaches, ASPE demonstrated an improved statistical power to detect regulatory variants with small effect sizes since allelic protein expression ratios are less prone to the influence of non-genetic regulators (e.g., diseases and inducers). This study suggests that the ASPE approach is a powerful tool for identifying cis-regulatory variants.

Insertion of a TRIM-like sequence in MdFLS2-1 promoter is associated with its allele-specific expression in response to Alternaria alternata in apple

Alternaria blotch disease, caused by Alternaria alternata apple pathotype (AAAP), is one of the major fungal diseases in apple. Early field observations revealed, the anther-derived homozygote Hanfu line (HFTH1) was highly susceptible to AAAP, whereas Hanfu (HF) exhibited resistance to AAAP. To understand the molecular mechanisms underlying the difference in sensitivity of HF and HFTH1 to AAAP, we performed allele-specific expression (ASE) analysis and comparative transcriptomic analysis before and after AAAP inoculation. We reported an important immune gene, namely, MdFLS2, which displayed strong ASE in HF with much lower expression levels of HFTH1-derived alleles. Transient overexpression of the dominant allele of MdFLS2-1 from HF in GL-3 apple leaves could enhance resistance to AAAP and induce expression of genes related to salicylic acid pathway. In addition, MdFLS2-1 was identified with an insertion of an 85-bp terminal-repeat retrotransposon in miniature (TRIM) element-like sequence in the upstream region of the nonreference allele. In contrast, only one terminal direct repeat (TDR) from TRIM-like sequence was present in the upstream region of the HFTH1-derived allele MdFLS2-2. Furthermore, the results of luciferase and β-glucuronidase reporter assays demonstrated that the intact TRIM-like sequence has enhancer activity. This suggested that insertion of the TRIM-like sequence regulates the expression level of the allele of MdFLS2, in turn, affecting the sensitivity of HF and HFTH1 to AAAP.

Large-scale genomic and transcriptomic profiles of rice hybrids reveal a core mechanism underlying heterosis

BACKGROUND

Heterosis is widely used in agriculture. However, its molecular mechanisms are still unclear in plants. Here, we develop, sequence, and record the phenotypes of 418 hybrids from crosses between two testers and 265 rice varieties from a mini-core collection.

RESULTS

Phenotypic analysis shows that heterosis is dependent on genetic backgrounds and environments. By genome-wide association study of 418 hybrids and their parents, we find that nonadditive QTLs are the main genetic contributors to heterosis. We show that nonadditive QTLs are more sensitive to the genetic background and environment than additive ones. Further simulations and experimental analysis support a novel mechanism, homo-insufficiency under insufficient background (HoIIB), underlying heterosis. We propose heterosis in most cases is not due to heterozygote advantage but homozygote disadvantage under the insufficient genetic background.

CONCLUSION

The HoIIB model elucidates that genetic background insufficiency is the intrinsic mechanism of background dependence, and also the core mechanism of nonadditive effects and heterosis. This model can explain most known hypotheses and phenomena about heterosis, and thus provides a novel theory for hybrid rice breeding in future.

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.

Analysis of the Whole-Genome Sequences from an Equus Parent-Offspring Trio Provides Insight into the Genomic Incompatibilities in the Hybrid Mule

Interspecific hybridization often shows negative effects on hybrids. However, only a few multicellular species, limited to a handful of plants and animals, have shown partial genetic mechanisms by which hybridization leads to low fitness in hybrids. Here, to explore the outcome of combining the two genomes of a horse and donkey, we analyzed the whole-genome sequences from an Equus parent-offspring trio using Illumina platforms. We generated 41.39× and 46.21× coverage sequences for the horse and mule, respectively. For the donkey, a 40.38× coverage sequence was generated and stored in our laboratory. Approximately 24.86 million alleles were discovered that varied from the reference genome. Single nucleotide polymorphisms were used as polymorphic markers for assigning alleles to their parental genomic inheritance. We identified 25,703 Mendelian inheritance error single nucleotide polymorphisms in the mule genome that were not inherited from the parents through Mendelian inheritance. A total of 555 de novo single nucleotide polymorphisms were also identified. The rate of de novo single nucleotide polymorphisms was 2.21 × 10^-7 in the mule from the Equus parent-offspring trio. This rate is obviously higher than the natural mutation rate for Equus, which is also consistent with the previous hypothesis that interracial crosses may have a high mutation rate. The genes associated with these single nucleotide polymorphisms are mainly involved in immune processes, DNA repair, and cancer processes. The results of the analysis of three genomes from an Equus parent-offspring trio improved our knowledge of the consequences of the integration of parental genomes in mules.

Development of the maize 5.5K loci panel for genomic prediction through genotyping by target sequencing

Genotyping platforms are important for genetic research and molecular breeding. In this study, a low-density genotyping platform containing 5.5K SNP markers was successfully developed in maize using genotyping by target sequencing (GBTS) technology with capture-in-solution. Two maize populations (Pop1 and Pop2) were used to validate the GBTS panel for genetic and molecular breeding studies. Pop1 comprised 942 hybrids derived from 250 inbred lines and four testers, and Pop2 contained 540 hybrids which were generated from 123 new-developed inbred lines and eight testers. The genetic analyses showed that the average polymorphic information content and genetic diversity values ranged from 0.27 to 0.38 in both populations using all filtered genotyping data. The mean missing rate was 1.23% across populations. The Structure and UPGMA tree analyses revealed similar genetic divergences (76-89%) in both populations. Genomic prediction analyses showed that the prediction accuracy of reproducing kernel Hilbert space (RKHS) was slightly lower than that of genomic best linear unbiased prediction (GBLUP) and three Bayesian methods for general combining ability of grain yield per plant and three yield-related traits in both populations, whereas RKHS with additive effects showed superior advantages over the other four methods in Pop1. In Pop1, the GBLUP and three Bayesian methods with additive-dominance model improved the prediction accuracies by 4.89-134.52% for the four traits in comparison to the additive model. In Pop2, the inclusion of dominance did not improve the accuracy in most cases. In general, low accuracies (0.33-0.43) were achieved for general combing ability of the four traits in Pop1, whereas moderate-to-high accuracies (0.52-0.65) were observed in Pop2. For hybrid performance prediction, the accuracies were moderate to high (0.51-0.75) for the four traits in both populations using the additive-dominance model. This study suggests a reliable genotyping platform that can be implemented in genomic selection-assisted breeding to accelerate maize new cultivar development and improvement.

Dominant complementation of biological pathways in maize hybrid lines is associated with heterosis

Allele-specific expressed genes (ASEGs) are widespread in maize hybrid lines and play important roles of complementation of biological pathways in heterosis. Heterosis (hybrid vigor) is an important phenomenon with both theoretical and practical value. However, our understanding of the genetic and molecular mechanisms behind heterosis is still limited. Here, we analyzed a comprehensive dataset of maize (Zea mays L.), including RNA-seq data from three hybrid-parent triplets (HPTs) and acetylated protein data from one HPT. The gene expression patterns exhibited extensive variation between the hybrids and their parents, and a substantial number of allele-specific expressed genes (ASEGs) were identified in the hybrids. Notably, ASEGs from different HPTs were significantly enriched in various conserved pathways. The parental alleles of ASEGs with fewer deleterious single-nucleotide polymorphisms were more likely to be expressed in hybrid lines than other parental alleles. ASEGs were mainly enriched in the functional gene ontology terms protein biosynthesis, photosynthesis, and metabolism. In addition, the ASEGs across the three HPTs were involved in key photosynthetic pathways and might enhance the photosynthetic efficiency of the hybrids. These findings suggest that ASEGs involved in complementary biological pathways in maize hybrids contribute to heterosis, shedding new light on the molecular mechanism of heterosis.

Transcriptome analysis provides critical answers to the "variants of uncertain significance" conundrum

While whole-genome and exome sequencing have transformed our collective understanding of genetics' role in disease pathogenesis, there are certain conditions and populations for whom DNA-level data fails to identify the underlying genetic etiology. Specifically, patients of non-White race and non-European ancestry are disproportionately affected by "variants of unknown/uncertain significance" (VUS), limiting the scope of precision medicine for minority patients and perpetuating health disparities. VUS often include deep intronic and splicing variants which are difficult to interpret from DNA data alone. RNA analysis can illuminate the consequences of VUS, thereby allowing for their reclassification as pathogenic versus benign. Here we review the critical role transcriptome analysis plays in clarifying VUS in both neoplastic and non-neoplastic diseases.

Transcriptome and Metabolome Analysis Provides Insights into the Heterosis of Yield and Quality Traits in Two Hybrid Rice Varieties (Oryza sativa L.)

Heterosis is a common biological phenomenon that is useful for breeding superior lines. Using heterosis to increase the yield and quality of crops is one of the main achievements of modern agricultural science. In this study, we analysed the transcriptome and metabolome of two three-line hybrid rice varieties, Taiyou 871 (TY871), and Taiyou 398 (TY398) and the parental grain endosperm using RNA-seq (three biological repeats per variety) and untargeted metabolomic (six biological repeats per variety) methods. TY871 and TY398 showed specific heterosis in yield and quality. Transcriptome analysis of the hybrids revealed 638 to 4059 differentially expressed genes in the grain when compared to the parents. Metabolome analysis of the hybrids revealed 657 to 3714 differential grain metabolites when compared to the parents. The honeydew1 and grey60 module core genes Os04g0350700 and Os05g0154700 are involved in the regulation of awn development, grain size, and grain number, as well as the regulation of grain length and plant height, respectively. Rice grain length may be an important indicator for improving the quality of three-line hybrid rice. In addition, the rice quality-related metabolite NEG_M341T662 was highly connected to the module core genes Os06g0254300 and Os03g0168100. The functions of Os06g0254300 and Os03g0168100 are EF-hand calcium binding protein and late embroideries absolute protein repeat containing protein, respectively. These genes may play a role in the formation of rice quality. We constructed a gene and metabolite coexpression network, which provides a scientific basis for the utilization of heterosis in producing high-yield and high-quality hybrid rice.

A Ghd7-centered regulatory network provides a mechanistic approximation to optimal heterosis in an elite rice hybrid

Heterosis refers to the superior performance of hybrids over their parents, which is a general phenomenon occurring in diverse organisms. Many commercial hybrids produce high yield without delayed flowering, which we refer to as optimal heterosis and is desired in hybrid breeding. Here, we attempted to illustrate the genomic basis of optimal heterosis by reinvestigating the single-locus quantitative trait loci and digenic interactions of two traits, the number of spikelets per panicle (SP) and heading date (HD), using recombinant inbred lines and 'immortalized F₂ s' derived from the elite rice (Oryza sativa) hybrid Shanyou 63. Our analysis revealed a regulatory network that may provide an approximation to the genetic constitution of the optimal heterosis observed in this hybrid. In this network, Ghd7 works as the core element, and three other genes, Ghd7.1, Hd1, and Hd3a/RFT1, also have major roles. The effects of positive dominance by Ghd7 and Ghd7.1 and negative dominance by Hd1 and Hd3a/RFT1 in the hybrid background contribute the major part to the high SP without delaying HD; numerous epistatic interactions, most of which involve Ghd7, also play important roles collectively. The results expand our understanding of the genic interaction networks underlying hybrid rice breeding programs, which may be very useful in future crop genetic improvement.

Non-coding RNA expression analysis revealed the molecular mechanism of flag leaf heterosis in inter-subspecific hybrid rice

Heterosis has been used widespread in agriculture, but its molecular mechanism is inadequately understood. Plants have a large number of non-coding RNAs (ncRNAs), among them, functional ncRNAs that have been studied widely containing long non-coding RNA (lncRNA) and circular RNA (circRNA) that play a role in varied biological processes, as well as microRNA (miRNA), which can not only regulate the post-transcriptional expression of target genes, but also target lncRNA and circRNA then participate the competing endogenous RNA (ceRNA) regulatory network. However, the influence of these three ncRNAs and their regulatory relationships on heterosis is unknown in rice. In this study, the expression profile of ncRNAs and the ncRNA regulatory network related to heterosis were comprehensively analyzed in inter-subspecific hybrid rice. A total of 867 miRNAs, 3,278 lncRNAs and 2,521 circRNAs were identified in the hybrid and its parents. Analysis of the global profiles of these three types of ncRNAs indicated that significant differences existed in the distribution and sequence characteristics of the corresponding genes. The numbers of miRNA and lncRNA in hybrid were higher than those in its parents. A total of 784 ncRNAs (169 miRNAs, 573 lncRNAs and 42 circRNAs) showed differentially expressed in the hybrid, and their target/host genes were vital in stress tolerance, growth and development in rice. These discoveries suggested that the expression plasticity of ncRNA has an important role of inter-subspecific hybrid rice heterosis. It is worth mentioning that miRNAs exhibited substantially more variations between hybrid and parents compared with observed variation for lncRNA and circRNA. Non-additive expression ncRNAs and allele-specific expression genes-related ncRNAs in hybrid were provided in this study, and multiple sets of ncRNA regulatory networks closely related to heterosis were obtained. Meanwhile, heterosis-related regulatory networks of ceRNA (lncRNA and circRNA) and miRNA were also demonstrated.

Genome-wide identification of the AlkB homologs gene family, PagALKBH9B and PagALKBH10B regulated salt stress response in Populus

The AlkB homologs (ALKBH) gene family regulates N⁶-methyladenosine (m⁶A) RNA methylation and is involved in plant growth and the abiotic stress response. Poplar is an important model plant for studying perennial woody plants. Poplars typically have a long juvenile period of 7-10 years, requiring long periods of time for studies of flowering or mature wood properties. Consequently, functional studies of the ALKBH genes in Populus species have been limited. Based on AtALKBHs sequence similarity with Arabidopsis thaliana, 23 PagALKBHs were identified in the genome of the poplar 84K hybrid genotype (P. alba  ×  P. tremula var. glandulosa), and gene structures and conserved domains were confirmed between homologs. The PagALKBH proteins were classified into six groups based on conserved sequence compared with human, Arabidopsis, maize, rice, wheat, tomato, barley, and grape. All homologs of PagALKBHs were tissue-specific; most were highly expressed in leaves. ALKBH9B and ALKBH10B are m⁶A demethylases and overexpression of their homologs PagALKBH9B and PagALKBH10B reduced m⁶A RNA methylation in transgenic lines. The number of adventitious roots and the biomass accumulation of transgenic lines decreased compared with WT. Therefore, PagALKBH9B and PagALKBH10B mediate m⁶A RNA demethylation and play a regulatory role in poplar growth and development. Overexpression of PagALKBH9B and PagALKBH10B can reduce the accumulation of H₂O₂ and oxidative damage by increasing the activities of SOD, POD, and CAT, and enhancing protection for Chl a/b, thereby increasing the salt tolerance of transgenic lines. However, overexpression lines were more sensitive to drought stress due to reduced proline content. This research revealed comprehensive information about the PagALKBH gene family and their roles in growth and development and responsing to salt stress of poplar.

The Genetic Basis of Gene Expression Divergence in Antennae of Two Closely Related Moth Species, Helicoverpa armigera and Helicoverpa assulta

The closely related species Helicoverpa armigera (H. armigera) and Helicoverpa assulta (H. assulta) have different host plant ranges and share two principal components of sex pheromones but with reversed ratios. The antennae are the main olfactory organ of insects and play a crucial role in host plant selection and mate seeking. However, the genetic basis for gene expression divergence in the antennae of the two species is unclear. We performed an allele-specific expression (ASE) analysis in the antennal transcriptomes of the two species and their F₁ hybrids, examining the connection between gene expression divergence and phenotypic differences. The results show that the proportion of genes classified as all cis was higher than that of all trans in males and reversed in females. The contribution of regulatory patterns to gene expression divergence in males was less than that in females, which explained the functional differentiation of male and female antennae. Among the five groups of F₁ hybrids, the fertile males from the cross of H. armigera female and H. assulta male had the lowest proportion of misexpressed genes, and the inferred regulatory patterns were more accurate. By using this group of F₁ hybrids, we discovered that cis-related regulations play a crucial role in gene expression divergence of sex pheromone perception-related proteins. These results are helpful for understanding how specific changes in the gene expression of olfactory-related genes can contribute to rapid evolutionary changes in important olfactory traits in closely related moths.

Transcriptome Analyses Indicate Significant Association of Increased Non-Additive and Allele-Specific Gene Expression with Hybrid Weakness in Rice (Oryza sativa L.)

The heterosis in hybrid rice is highly affected by the environment and hybrid weakness occurs frequently depending on the genotypes of the hybrid and its parents. Hybrid weakness was also observed in our field experiments on nine rice hybrids produced by 3 × 3 incomplete diallel crosses. Among the nine hybrids, five displayed mid-parent heterosis (MPH) for grain yield per plant, while four showed mid-parent hybrid weakness (MPHW). A sequencing analysis of transcriptomes in panicles at the seed-filling stage revealed a significant association between enhanced non-additive gene expression (NAE) and allele-specific gene expression (ASE) with hybrid weakness. High proportions of ASE genes, with most being of mono-allele expression, were detected in the four MPHW hybrids, ranging from 22.65% to 45.97%; whereas only 4.80% to 5.69% of ASE genes were found in the five MPH hybrids. Moreover, an independence test indicated that the enhancements of NAE and ASE in the MPHW hybrids were significantly correlated. Based on the results of our study, we speculated that an unfavorable environment might cause hybrid weakness by enhancing ASE and NAE at the transcriptome level.

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.

Transposon insertions regulate genome-wide allele-specific expression and underpin flower colour variations in apple (Malus spp.)

Allele-specific expression (ASE) can lead to phenotypic diversity and evolution. However, the mechanisms regulating ASE are not well understood, particularly in woody perennial plants. In this study, we investigated ASE genes in the apple cultivar 'Royal Gala' (RG). A high quality chromosome-level genome was assembled using a homozygous tetra-haploid RG plant, derived from anther cultures. Using RNA-sequencing (RNA-seq) data from RG flower and fruit tissues, we identified 2091 ASE genes. Compared with the haploid genome of 'Golden Delicious' (GD), a parent of RG, we distinguished the genomic sequences between the two alleles of 817 ASE genes, and further identified allele-specific presence of a transposable element (TE) in the upstream region of 354 ASE genes. These included MYB110a that encodes a transcription factor regulating anthocyanin biosynthesis. Interestingly, another ASE gene, MYB10 also showed an allele-specific TE insertion and was identified using genome data of other apple cultivars. The presence of the TE insertion in both MYB genes was positively associated with ASE and anthocyanin accumulation in apple petals through analysis of 231 apple accessions, and thus underpins apple flower colour evolution. Our study demonstrated the importance of TEs in regulating ASE on a genome-wide scale and presents a novel method for rapid identification of ASE genes and their regulatory elements in plants.

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.

Favoring Expression of Yak Alleles in Interspecies F1 Hybrids of Cattle and Yak Under High-Altitude Environments

Both cis- and trans-regulation could cause differential expression between the parental alleles in diploid species that might have broad biological implications. Due to the relatively distant genetic divergence between cattle and yak, as well as their differential adaptation to high-altitude environments, we investigated genome-wide allelic differential expression (ADE) in their F1 hybrids using Nanopore long-read RNA-seq technology. From adult F1 hybrids raised in high-altitude, ten lung and liver tissues were individually sequenced for producing 31.6 M full-length transcript sequences. Mapping against autosomal homologous regions between cattle and yak, we detected 17,744 and 14,542 protein-encoding genes expressed in lung and liver tissues, respectively. According to the parental assignments of transcript sequences, a total of 3,381 genes were detected to show ADE in at least one sample. There were 186 genes showing ubiquitous ADE in all the studied animals, and among them 135 and 37 genes had consistent higher expression of yak and cattle alleles, respectively. Functional analyses revealed that the genes with favoring expression of yak alleles have been involved in the biological progresses related with hypoxia adaptation and immune response. In contrast, the genes with favoring expression of cattle alleles have been enriched into different biological progresses, such as secretion of endocrine hormones and lipid metabolism. Our results would support unequal contribution of parental genes to environmental adaptation in the F1 hybrids of cattle and yak.