Transcriptome Analysis of Interspecific Hybrid between Brassica napus and B. rapa Reveals Heterosis for Oil Rape Improvement

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.

Transcriptome Analysis Reveals a Major Gene Expression Pattern and Important Metabolic Pathways in the Control of Heterosis in Chinese Cabbage

Although heterosis is commonly used in Chinese cabbage, its molecular basis is poorly understood. In this study, 16Chinese cabbage hybrids were utilized as test subjects to explore the potential molecular mechanism of heterosis. RNA sequencing revealed 5815-10,252 differentially expressed genes (DEGs) (female parent vs. male parent), 1796-5990 DEGs (female parent-vs-hybrid), and 2244-7063 DEGs (male parent vs. hybrid) in 16 cross combinations at the middle stage of heading. Among of them, 72.83-84.20% DEGs conformed to the dominant expression pattern, which is the predominant expression pattern in hybrids. There were 13 pathways in which DEGs were significantly enriched in most cross combinations. Among them, the plant-pathogen interaction (ko04626) and circadian rhythm-plant (ko04712)were significantly enriched by DEGs in strong heterosis hybrids. WGCNA also proved that the two pathways were significantly related to heterosis in Chinese cabbage.

Identification of the Exercise and Time Effects on Human Skeletal Muscle through Bioinformatics Methods

Background

The human body has more than 600 kinds of skeletal muscles, which accounts for about 40% of the whole weight. Most skeletal muscles can make bones move, and their strength and endurance directly affect their performance during exercise.

Methods

To determine the effects of exercise and time on human skeletal muscle, we downloaded the microarray expression profile of GSE1832 and analyzed it to select differentially expressed genes (DEGs). Then, a protein-protein interaction (PPI) network was established, and the hub genes were identified. Afterwards, DEGs were applied to perform Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Finally, with the help of Gene Set Enrichment Analysis (GSEA), the gene sets in the 7 samples were enriched in the KEGG pathway.

Results

Through a series of bioinformatics analyses, we obtained a total of 271 DEGs. After that, four hub genes were determined through the PPI network, namely, EP300, STAT1, CDKN1A, and RAC2. In addition, we got that these DEGs were enriched in GO, such as regulation of cell population proliferation, cellular water homeostasis, and so on, and in KEGG, namely, hepatitis B, Epstein-Barr virus infection, small cell lung cancer, pathways in cancer, and others. Finally, the gene set in the samples obtained by GSEA was enriched in the cell cycle, chemokine signaling pathway, DNA replication, cytokine receptor interaction, ECM receptor interaction, and focal adhesion in KEGG.

Conclusion

The findings obtained in this study will provide new clues for elucidating the mechanism of exercise and time on human skeletal muscles.

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.

Cross-kingdom transcriptomic trends in the evolution of hybrid gene expression

Hybridization is a route to speciation that occurs widely across the eukaryote tree of life. The success of allopolyploids (hybrid species with increased ploidy) and homoploid hybrids (with unchanged ploidy) is well documented. However, their formation and establishment is not straightforward, with a suite of near‐instantaneous and longer term biological repercussions faced by the new species. Central to these challenges is the rewiring of gene regulatory networks following the merger of distinct genomes inherited from both parental species. Research on the evolution of hybrid gene expression has largely involved studies on a single hybrid species or a few gene families. Here, we present the first standardized transcriptome‐wide study exploring the fates of genes following hybridization across three kingdoms: animals, plants and fungi. Within each kingdom, we pair an allopolyploid system with a closely related homoploid hybrid to decouple the influence of increased ploidy from genome merger. Genome merger, not changes in ploidy, has the greatest effect on posthybridization expression patterns across all study systems. Strikingly, we find that differentially expressed genes in parent species preferentially switch to more similar expression in hybrids across all kingdoms, likely as a consequence of regulatory trans‐acting cross‐talk within the hybrid nucleus. We also highlight the prevalence of gene loss or silencing among extremely differentially expressed genes in hybrid species across all kingdoms. These shared patterns suggest that the evolutionary process of hybridization leads to common high‐level expression outcomes, regardless of the particular species or kingdom.

Genome-Wide Identification and Analysis of the Ascorbate Peroxidase (APX) Gene Family of Winter Rapeseed (Brassica rapa L.) Under Abiotic Stress

Winter Brassica rapa (B. rapa) is an important oilseed crop in northern China, but the mechanism of its cold resistance remains unclear. Ascorbate peroxidase (APX) plays important roles in the response of this plant to abiotic stress and in scavenging free radicals. In this study, the roles of APX proteins in the cold response and superoxide metabolism pathways in rapeseed species were investigated, and a comprehensive analysis of phylogeny, chromosome distribution, motif identification, sequence structure, gene duplication, and RNA-seq expression profiles in the APX gene family was conducted. Most BrAPX genes were specifically expressed under cold stress and behaved significantly differently in cold-tolerant and weakly cold-resistant varieties. Quantitative real-time-PCR (qRT-PCR) was also used to verify the differences in expression between these two varieties under cold, freezing, drought and heat stress. The expression of five BrAPX genes was significantly upregulated in growth cones at 3 h of cold stress, while their expression was significantly lower at 24 h than at 3 h. The expression of Bra015403 and Bra003918 was significantly higher in “Longyou-7” growth cones than in other treatments. Five BrAPXs (Bra035235, Bra003918, Bra033040, Bra017120, and Bra031934) were closely associated with abiotic stress responses in B. rapa. These candidate genes may play important roles in the response of B. rapa to low temperature stress and provide new information for the elucidation of the cold resistance mechanism in B. rapa.

Transcriptome-Wide Analyses Identify Dominant as the Predominantly Non-Conservative Alternative Splicing Inheritance Patterns in F1 Chickens

Transcriptome analysis has been used to investigate many economically traits in chickens; however, alternative splicing still lacks a systematic method of study that is able to promote proteome diversity, and fine-tune expression dynamics. Hybridization has been widely utilized in chicken breeding due to the resulting heterosis, but the dynamic changes in alternative splicing during this process are significant yet unclear. In this study, we performed a reciprocal crossing experiment involving the White Leghorn and Cornish Game chicken breeds which exhibit major differences in body size and reproductive traits, and conducted RNA sequencing of the brain, muscle, and liver tissues to identify the inheritance patterns. A total of 40 515 and 42 612 events were respectively detected in the brain and muscle tissues, with 39 843 observed in the liver; 2807, 4242, and 4538 events significantly different between two breeds were identified in the brain, muscle, and liver tissues, respectively. The hierarchical cluster of tissues from different tissues from all crosses, based on the alternative splicing profiles, suggests high tissue and strain specificity. Furthermore, a comparison between parental strains and hybrid crosses indicated that over one third of alternative splicing genes showed conserved patterns in all three tissues, while the second prevalent pattern was non-additive, which included both dominant and transgressive patterns; this meant that the dominant pattern plays a more important role than suppression. Our study provides an overview of the inheritance patterns of alternative splicing in layer and broiler chickens, to better understand post-transcriptional regulation during hybridization.

Exploring the gene pool of Brassica napus by genomics-based approaches

De novo allopolyploidization in Brassica provides a very successful model for reconstructing polyploid genomes using progenitor species and relatives to broaden crop gene pools and understand genome evolution after polyploidy, interspecific hybridization and exotic introgression. B. napus (AACC), the major cultivated rapeseed species and the third largest oilseed crop in the world, is a young Brassica species with a limited genetic base resulting from its short history of domestication, cultivation, and intensive selection during breeding for target economic traits. However, the gene pool of B. napus has been significantly enriched in recent decades that has been benefit from worldwide effects by the successful introduction of abundant subgenomic variation and novel genomic variation via intraspecific, interspecific and intergeneric crosses. An important question in this respect is how to utilize such variation to breed crops adapted to the changing global climate. Here, we review the genetic diversity, genome structure, and population-level differentiation of the B. napus gene pool in relation to known exotic introgressions from various species of the Brassicaceae, especially those elucidated by recent genome-sequencing projects. We also summarize progress in gene cloning, trait-marker associations, gene editing, molecular marker-assisted selection and genome-wide prediction, and describe the challenges and opportunities of these techniques as molecular platforms to exploit novel genomic variation and their value in the rapeseed gene pool. Future progress will accelerate the creation and manipulation of genetic diversity with genomic-based improvement, as well as provide novel insights into the neo-domestication of polyploid crops with novel genetic diversity from reconstructed genomes.

Comparative Transcriptomic Analysis of Gene Expression Inheritance Patterns Associated with Cabbage Head Heterosis

The molecular mechanism of heterosis or hybrid vigor, where F1 hybrids of genetically diverse parents show superior traits compared to their parents, is not well understood. Here, we studied the molecular regulation of heterosis in four F1 cabbage hybrids that showed heterosis for several horticultural traits, including head size and weight. To examine the molecular mechanisms, we performed a global transcriptome profiling in the hybrids and their parents by RNA sequencing. The proportion of genetic variations detected as single nucleotide polymorphisms and small insertion–deletions as well as the numbers of differentially expressed genes indicated a larger role of the female parent than the male parent in the genetic divergence of the hybrids. More than 86% of hybrid gene expressions were non-additive. More than 81% of the genes showing divergent expressions showed dominant inheritance, and more than 56% of these exhibited maternal expression dominance. Gene expression regulation by cis-regulatory mechanisms appears to mediate most of the gene expression divergence in the hybrids; however, trans-regulatory factors appear to have a higher effect compared to cis-regulatory factors on parental expression divergence. These observations bring new insights into the molecular mechanisms of heterosis during the cabbage head development.

Potential of the C Genome of the Different Variants of Brassica oleracea for Heterosis in Spring B. napus Canola

The genetic base of Brassica napus canola need to be broadened for exploitation of heterosis at a greater level in the breeding of F1 hybrid canola cultivars. In this study, we evaluated 228 inbred B. napus canola lines derived from six B. napus × B. oleracea interspecific crosses and following two breeding methods (F2- and BC1-derived lines) to understand the effect of the B. oleracea alleles on heterosis for different agronomic and seed quality traits. Test hybrids of the inbreds derived from crosses involving vars. botrytis (cauliflower), alboglabra (Chinese kale) and capitata (cabbage) cv. Badger Shipper, on an average, gave about 10% mid-parent heterosis (MPH), and about 67% of the test hybrids gave higher seed yield than the common B. napus parent indicating that B. oleracea alleles can contribute to heterosis for seed yield in spring B. napus canola hybrids. This was also evident from a positive correlation of the genetic distance of the inbred lines from the common B. napus parent with MPH for seed yield (r = 0.31) as well as with hybrid yield (r = 0.26). Almost no correlation was found between genetic distance and MPH for seed oil and protein content as well as with the performance of the test hybrids for these two traits. The occurrence of positive correlation between seed yield of the inbred lines and test hybrids suggested the importance of the genes exerting additive effect for high seed yield in the hybrids. Very little or almost no heterosis was found for the other agronomic traits as well as for seed oil and protein content. While comparing the two breeding methods, no significant difference was found for seed yield of the test hybrids or the level of MPH; however, the BC1-derived inbred and test hybrid populations flowered and matured earlier and had longer grain-filling period than the F2-derived population. Thus, the results suggested that the B. oleracea gene pool can be used in the breeding of spring B. napus canola to improve seed yield in hybrid cultivars.

Effects of parental genetic divergence on gene expression patterns in interspecific hybrids of Camellia

BACKGROUND

The merging of two divergent genomes during hybridization can result in the remodeling of parental gene expression in hybrids. A molecular basis underling expression change in hybrid is regulatory divergence, which may change with the parental genetic divergence. However, there still no unanimous conclusion for this hypothesis.

RESULTS

Three species of Camellia with a range of genetic divergence and their F₁ hybrids were used to study the effect of parental genetic divergence on gene expression and regulatory patterns in hybrids by RNA-sequencing and allelic expression analysis. We found that though the proportion of differentially expressed genes (DEGs) between the hybrids and their parents did not increase, a greater proportion of DEGs would be non-additively (especially transgressively) expressed in the hybrids as genomes between the parents become more divergent. In addition, the proportion of genes with significant evidence of cis-regulatory divergence increased, whereas with trans-regulatory divergence decreased with parental genetic divergence.

CONCLUSIONS

The discordance within hybrid would intensify as the parents become more divergent, manifesting as more DEGs would be non-additively expressed. Trans-regulatory divergence contributed more to the additively inherited genes than cis, however, its contribution to expression difference would be weakened as cis mutations accumulated over time; and this might be an important reason for that the more divergent the parents are, the greater proportion of DEGs would be non-additively expressed in hybrid.

Transcriptome analysis reveals hybridization-induced genome shock in an interspecific F1 hybrid from Camellia

The combination of two divergent genomes during hybridization can result in "genome shock". Although genome shock has been reported in the hybrids of some herbaceous plants, the pattern and the principle it follows are far from understood, especially in woody plants. Here, the gene expression patterns were remodeled in the F₁ hybrid from the crossing of Camellia azalea × Camellia amplexicaulis compared with the parents as revealed by RNA-seq. About 54.5% of all unigenes were differentially expressed between the F₁ hybrid and at least one of the parents, including 6404 unigenes with the highest expression level in the F₁ hybrid. A series of genes, related to flower development, essential for RNA-directed DNA methylation and histone methylation, as well as 223 transposable elements, were enriched; and most of them exhibited a higher level of expression in the F₁ hybrid. These results indicated that the genome shock induced by interspecific hybridization in Camellia could indeed result in changes of gene expression patterns, potentially through regulating DNA methylation and histone methylation which may be helpful for the maintaining of genome stability and even related to the unique phenotype of the F₁ hybrid.

IntroMap: a signal analysis based method for the detection of genomic introgressions

BACKGROUND

Breeding programs often rely on marker-assisted tests or variant calling of next generation sequence (NGS) data to identify regions of genomic introgression arising from the hybridization of two plant species. In this paper we present IntroMap, a bioinformatics pipeline for the screening of candidate plants through the application of signal processing techniques to NGS data, using alignment to a reference genome sequence (annotation is not required) that shares homology with the recurrent parental cultivar, and without the need for de novo assembly of the read data or variant calling.

RESULTS

We show the accurate identification of introgressed genomic regions using both in silico simulated genomes, and a hybridized cultivar data set using our pipeline. Additionally we show, through targeted marker-based assays, validation of the IntroMap predicted regions for the hybrid cultivar.

CONCLUSIONS

This approach can be used to automate the screening of large populations, reducing the time and labor required, and can improve the accuracy of the detection of introgressed regions in comparison to a marker-based approach. In contrast to other approaches that generally rely upon a variant calling step, our method achieves accurate identification of introgressed regions without variant calling, relying solely upon alignment.

Transcriptome and functional analysis reveals hybrid vigor for oil biosynthesis in oil palm

Oil palm is the most productive oil crop in the world and composes 36% of the world production. However, the molecular mechanisms of hybrids vigor (or heterosis) between Dura, Pisifera and their hybrid progeny Tenera has not yet been well understood. Here we compared the temporal and spatial compositions of lipids and transcriptomes for two oil yielding organs mesocarp and endosperm from Dura, Pisifera and Tenera. Multiple lipid biosynthesis pathways are highly enriched in all non-additive expression pattern in endosperm, while cytokinine biosynthesis and cell cycle pathways are highly enriched both in endosperm and mesocarp. Compared with parental palms, the high oil content in Tenera was associated with much higher transcript levels of EgWRI1, homolog of Arabidopsis thaliana WRINKLED1. Among 338 identified genes in lipid synthesis, 207 (61%) has been identified to contain the WRI1 specific binding AW motif. We further functionally identified EgWRI1-1, one of three EgWRI1 orthologs, by genetic complementation of the Arabidopsis wri1 mutant. Ectopic expression of EgWRI1-1 in plant produced dramatically increased seed mass and oil content, with oil profile changed. Our findings provide an explanation for EgWRI1 as an important gene contributing hybrid vigor in lipid biosynthesis in oil palm.