Transcriptomic analyses reveal molecular mechanisms underlying growth heterosis and weakness of rubber tree seedlings

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.

Revisiting the amalgaviral landscapes in plant transcriptomes expands the host range of plant amalgaviruses

Plant amalgaviruses are monopartite, double-stranded RNA viruses, capable of vertical transmission through seeds. An attempt to revisit plant transcriptome-assembled contigs for amalgaviral sequences identified 40 putative novel amalgaviruses in 35 plant species, nearly doubling the number of plant amalgaviruses. Of the 35 plant species, 33 are reported to host amalgaviruses for the first time, including a pteridophytic and two gymnospermic species. Coding-complete genomes of all identified viruses were recovered and the putative +1 programmed ribosomal frameshift (PRF) sites were determined. Genomes of 35 identified amalgaviruses contained the conserved +1 PRF motif 'UUU_CGN', while variant versions were predicted in five genomes. Phylogenetic analysis grouped pteridophyte- and gymnosperm-infecting amalgaviruses together in divergent sub-clades while few of the related angiosperm-infecting amalgaviruses infect members of the same plant family, reiterating the co-evolution of plant amalgaviruses and their hosts. The current study paves way for further studies on understanding biological properties of identified viruses.

Dynamic analysis of Tapping Panel Dryness in Hevea brasiliensis reveals new insights on this physiological syndrome affecting latex production

Tapping Panel Dryness (TPD) is a physiological disorder affecting natural rubber production in Hevea brasiliensis. TPD is associated with clonal susceptibility and overexploitation of rubber trees. Most studies are based on a binary point view of the absence or presence of TPD. This study sets out to characterize the dynamic of the TPD onset through the monthly monitoring of the dry cut length. This reveals the presence of dry spots on the tapped panel of any trees. The frequency of these dry spots increases dramatically in trees developing high level of TPD. Brown bast is an irreversible form of TPD. Brown bast is correlated to a high level of dry cut length. Application of an intensive harvesting system induces early TPD occurrence, which facilitates the study of TPD. Among latex diagnosis parameters, only sucrose content is significantly associated with TPD. Other parameters are more prone to environmental effects and are not reliable as physiological markers. These findings explain the contradictory conclusions of some papers. This study suggests to use intensive harvesting system and monitor the dry cut length for genetic analysis of TPD.

Gene expression variation explains maize seed germination heterosis

BACKGROUND

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

RESULTS

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

Identification and Functional Evaluation of Three Polyubiquitin Promoters from Hevea brasiliensis

Hevea brasiliensis is an economically important tree species that provides the only commercial source of natural rubber. The replacement of the CaMV35S promoter by endogenous polyubiquitin promoters may be a viable way to improve the genetic transformation of this species. However, no endogenous polyubiquitin promoters in Hevea have been reported yet. Here, we identified three Hevea polyubiquitin genes HbUBI10.1, HbUBI10.2 and HbUBI10.3, which encode ubiquitin monomers having nearly identical amino acid sequences to that of AtUBQ10. The genomic fragments upstream of these HbUBI genes, including the signature leading introns, were amplified as putative HbUBI promoters. In silico analysis showed that a number of cis-acting elements which are conserved within strong constitutive polyubiquitin promoters were presented in these HbUBI promoters. Transcriptomic data revealed that HbUBI10.1 and HbUBI10.2 had a constitutive expression in Hevea plants. Semi-quantitative RT-PCR showed that these three HbUBI genes were expressed higher than the GUS gene driven by CaMV35S in transgenic Hevea leaves. All three HbUBI promoters exhibited the capability to direct GFP expression in both transient and stable transformation assays, although they produced lower protoplast transformation efficiencies than the CaMV35S promoter. These HbUBI promoters will expand the availability of promoters for driving the transgene expression in Hevea genetic engineering.

Transcriptomic Characterization of Miscanthus sacchariflorus × M. lutarioriparius and Its Implications for Energy Crop Development in the Semiarid Mine Area

Miscanthus interspecific hybrids have been proved to have better adaptability in marginal lands than their parents. Miscanthus sacchariflorus and Miscanthus lutarioriparius were used as the parents to develop hybrids. We performed the transcriptome for 110 F1 hybrids of Miscanthus sacchariflorus × Miscanthus lutarioriparius and their parents that had been established on the Loess Plateau mine area, to estimate the population’s genetic expression variation, and illuminate the adaptive mechanism of the F1 population. The result speculated that the F1 population has mainly inherited the stress response metabolic pathway of its female parent (M. sacchariflorus), which may be responsible for its higher environmental adaptability and biomass yield compared with male parents. Based on PopART, we assembled a leaf reference transcriptome for M. sacchariflorus (LRTMS) and obtained 8116 high-quality transcripts. When we analyze the differential expression of genes between F1 population and parent, 39 and 56 differentially expressed genes were screened out in the female parent and male parent, respectively. The enrichment analysis showed that pathways of carbohydrate metabolism, lipid metabolism, biosynthesis of secondary metabolites and circadian rhythm-plant played a key role in resisting the harsh environment. The carbohydrate metabolism and lipid metabolism were also significantly enriched, and the synthesis of these substances facilitated the yield. The results provided an insight into breeding Miscanthus hybrids more suited to the harsh environment of the Loess Plateau.

Integrated Analysis of mRNA and Non-coding RNA Transcriptome in Pepper (Capsicum chinense) Hybrid at Seedling and Flowering Stages

Pepper is an important vegetable in the world. In this work, mRNA and ncRNA transcriptome profiles were applied to understand the heterosis effect on the alteration in the gene expression at the seedling and flowering stages between the hybrid and its parents in Capsicum chinense. Our phenotypic data indicated that the hybrid has dominance in leaf area, plant scope, plant height, and fruit-related traits. Kyoto Encyclopedia of Genes and Genomes analysis showed that nine members of the plant hormone signal transduction pathway were upregulated in the seedling and flowering stages of the hybrid, which was supported by weighted gene coexpression network analysis and that BC332_23046 (auxin response factor 8), BC332_18317 (auxin-responsive protein IAA20), BC332_13398 (ethylene-responsive transcription factor), and BC332_27606 (ethylene-responsive transcription factor WIN1) were candidate hub genes, suggesting the important potential role of the plant hormone signal transduction in pepper heterosis. Furthermore, some transcription factor families, including bHLH, MYB, and HSF were greatly over-dominant. We also identified 2,525 long ncRNAs (lncRNAs), 47 micro RNAs (miRNAs), and 71 circle RNAs (circRNAs) in the hybrid. In particular, downregulation of miR156, miR169, and miR369 in the hybrid suggested their relationship with pepper growth vigor. Moreover, we constructed some lncRNA–miRNA–mRNA regulatory networks that showed a multi-dimension to understand the ncRNA relationship with heterosis. These results will provide guidance for a better understanding of the molecular mechanism involved in pepper heterosis.

Comparative analysis of latex transcriptomes reveals the potential mechanisms underlying rubber molecular weight variations between the Hevea brasiliensis clones RRIM600 and Reyan7-33-97

BACKGROUND

The processabilities and mechanical properties of natural rubber depend greatly on its molecular weight (MW) and molecular weight distribution (MWD). However, the mechanisms underlying the regulation of molecular weight during rubber biosynthesis remain unclear.

RESULTS

In the present study, we determined the MW and particle size of latex from 1-year-old virgin trees and 30-year-old regularly tapped trees of the Hevea clones Reyan7-33-97 and RRIM600. The results showed that both the MW and the particle size of latex varied between these two clones and increased with tree age. Latex from RRIM600 trees had a smaller average particle size than that from Reyan7-33-97 trees of the same age. In 1-year-old trees, the Reyan7-33-97 latex displayed a slightly higher MW than that of RRIM600, whereas in 30-year-old trees, the RRIM600 latex had a significantly higher MW than the Reyan7-33-97 latex. Comparative analysis of the transcriptome profiles indicated that the average rubber particle size is negatively correlated with the expression levels of rubber particle associated proteins, and that the high-MW traits of latex are closely correlated with the enhanced expression of isopentenyl pyrophosphate (IPP) monomer-generating pathway genes and downstream allylic diphosphate (APP) initiator-consuming non-rubber pathways. By bioinformatics analysis, we further identified a group of transcription factors that potentially regulate the biosynthesis of IPP.

CONCLUSIONS

Altogether, our results revealed the potential regulatory mechanisms involving gene expression variations in IPP-generating pathways and the non-rubber isoprenoid pathways, which affect the ratios and contents of IPP and APP initiators, resulting in significant rubber MW variations among same-aged trees of the Hevea clones Reyan7-33-97 and RRIM600. Our findings provide a better understanding of rubber biosynthesis and lay the foundation for genetic improvement of rubber quality in H. brasiliensis.

Comparative analyses of hypothalamus transcriptomes reveal fertility-, growth-, and immune-related genes and signal pathways in different ploidy cyprinid fish

Triploid crucian carp (TCC) is obtained by hybridization of female diploid red crucian carp (Carassius auratus red var., RCC) and male allotetraploid hybrids. In this study, high-throughput sequencing was used to conduct the transcriptome analysis of the female hypothalamus of diploid RCC, diploid common carp (Cyprinus carpio L., CC) and TCC. The key functional expression genes of the hypothalamus were obtained through functional gene annotation and differential gene expression screening. A total of 71.56 G data and 47,572 genes were obtained through sequencing and genome mapping, respectively. The Fuzzy Analysis Clustering assigned the differentially expressed genes (DEGs) into eight groups, two of which, overdominance expression (6005, 12.62%) and underdominance expression (3849, 8.09%) in TCC were further studied. KEGG enrichment analysis showed that the DEGs in overdominance were mainly enriched in four pathways. The expression of several fertility-related genes was lower levels in TCC, whereas the expression of several growth-related genes and immune-related genes was higher levels in TCC. Besides, 15 DEGs were verified by quantitative real-time PCR (qPCR). The present study can provide a reference for breeding sterility, fast-growth, and disease-resistant varieties by distant hybridization.

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

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

Transcriptomic Profiling Provides Insights into Inbreeding Depression in Yesso Scallop Patinopecten yessoensis

Inbreeding often causes a decline in biological fitness, known as inbreeding depression. In genetics study, inbreeding coefficient f gives the proportion by which the heterozygosity of an individual is reduced by inbreeding. With the development of high-throughput sequencing, researchers were able to perform deep approaches to investigate which genes are affected by inbreeding and reveal some molecular underpinnings of inbreeding depression. As one commercially important species, Yesso scallop Patinopecten yessoensis confront the same dilemma of inbreeding depression. To examine how inbreeding affects gene expression, we compared the transcriptome of two experimentally selfing families with inbreeding coefficient f reached 0.5 as well as one natural population (f ≈ 0) of P. yessoensis. A total of 24 RNA-Seq libraries were constructed using scallop adductor muscle, and eventually 676.56 M (96.85%) HQ reads were acquired. Based on differential gene analysis, we were able to identify nine common differentially expressed genes (DEGs) across the top-ranked 30 DEGs in both selfing families in comparation with the natural population. Remarkable, through weighted gene co-expression network analysis (WGCNA), five common DEGs were found enriched in the most significant inbreeding related functional module M14 (FDR = 1.64E-156), including SREBP1, G3BP2, SBK1, KIAA1161, and AATs-Glupro. These five genes showed significantly higher expression in self-bred progeny. Suggested by the genetic functional analysis, up-regulated SREBP1, G3BP2, and KIAA1161 may suggest a perturbing lipid metabolism, a severe inframammary reaction or immune response, and a stress-responsive behavior. Besides, the significant higher SBK1 and AATs-Glupro may reflect the abnormal cellular physiological situation. Together, these genetic aberrant transcriptomic performances may contribute to inbreeding depression in P. yessoensis, deteriorating the stress tolerance and survival phenotype in self-bred progeny. Our results would lay a foundation for further comprehensive understanding of bivalve inbreeding depression, which may potentially benefit the genetic breeding for scallop aquaculture.

Comparative transcriptome analysis provides key insights into seedling development in switchgrass (Panicum virgatum L.)

BACKGROUND

Switchgrass (Panicum virgatum L.), a warm-season perennial C4 plant, can be used as a forage plant, a soil and water conservation plant, a windbreak plant, and as a good source of biofuels and alternative energy with low planting costs. However, switchgrass exhibits low rates of seedling development compared to other crops, which means it is typically out-competed by weeds. There is a large variation in seedling development rates among different plantlets in switchgrass, which limits its usefulness for large-scale cultivation. Little is currently known about the molecular reasons for slow seedling growth.

RESULTS

Characterization of the seedling development process via growth indices indicated a relatively stagnant growth stage in switchgrass. A total of 678 differentially expressed genes (DEGs) were identified from the comparison of transcriptomes from slowly developed (sd) and rapidly developed (rd) switchgrass seedlings. Gene ontology and pathway enrichment analysis showed that DEGs were enriched in diterpenoid biosynthesis, thiamine metabolism, and circadian rhythm. Transcription factor enrichment and expression analyses showed MYB-related, bHLH and NAC family genes were essential for seedling growth. The transcriptome results were consistent with those of quantitative real-time polymerase chain reaction. Then, the expression profiles of maize and switchgrass were compared during seedling leaf development. A total of 128 DEGs that play key roles in seedling growth were aligned to maize genes. Transcriptional information and physiological indices suggested that several genes involved in the circadian rhythm, thiamine metabolism, energy metabolism, gibberellic acid biosynthesis, and signal transduction played important roles in seedling development.

CONCLUSIONS

The seedling development process of switchgrass was characterized, and the molecular differences between slowly developed and rapidly developed seedlings were discussed. This study provides new insights into the reasons for slow seedling development in switchgrass and will be useful for the genetic improvement of switchgrass and other crops.