Weighted gene co-expression network analysis identifies genes related to anthocyanin biosynthesis and functional verification of hub gene SmWRKY44

Transcriptome analysis reveals the synergistic involvement of MGLL and LPIN1 in fatty acid synthesis in broiler pectoral muscles

Fatty acids (FAs) are one of the most important bioactive compounds affecting the quality of meat. In this study, we compared the expression profiles of genes involved in FA production in the breast muscle of Jingxing Yellow chickens at different days of age determined by transcriptomic analysis to identify key genes and pathways regulating the FA composition of the breast muscle. Through clustering analysis of gene expression data, the growth process of broiler chickens can be divided into two stages, namely the growth and development stage at the 35th and 63rd days of age (D35, D63), and the mature stage at the 119th day of age (D119). The content of some important unsaturated fatty acids (UFAs), such as C18:2n6c, C20:4n6, and C22:6n3, in the pectoral muscles, differed significantly between these two stages (p < 0.05). Therefore, we compared the gene expression profiles at D35 and D63 with those at D119, and identified differentially expressed genes (DEGs). The gene modules related to the five UFAs with significant changes were identified by weighted gene co-expression network analysis (WGCNA), and then 150 crossover genes were identified by crossover analysis of the detected DEGs and WGCNA. The results of the pathway enrichment analysis revealed the glycerolipid metabolism pathway related to lipid metabolism, in which the MGLL and LPIN1 genes were particularly enriched. In this study, the expression levels of MGLL and LPIN1 showed an increasing trend during the growth process of broilers, with a negative regulatory effect on the significantly reduced content of C18:2n6c in the pectoral muscle, and a positive regulatory effect on the significantly increased content of C20:4n6. These findings indicated that MGLL and LPIN1 synergistically promote the deposition of FAs, which may further promote the conversion of linoleic acid (C18:2n6c) to arachidonic acid (C20:4n6). Therefore, screening and identifying FA production-related functional genes are key to elucidate the regulatory molecular mechanism of production of FAs in chicken muscle, aiming to provide a theoretical basis for improving chicken meat quality.

A Functional InDel in the WRKY10 Promoter Controls the Degree of Flesh Red Pigmentation in Apple

MYB transcription factors have been linked to anthocyanin synthesis and various color phenotypes in plants. In apple, MYB10 confers a red-flesh phenotype due to a minisatellite insertion in its R6 promoter, but R6:MYB10 genotypes exhibit various degrees of red pigmentation in the flesh, suggesting the involvement of other genetic factors. Here, it is shown that MdWRKY10, a transcription factor identified via DNA pull-down trapping, binds to the promoter of MdMYB10 and activates its transcription. MdWRKY10 specifically interacts with the WDR protein MdTTG1 to join the apple MYB-bHLH-WDR (MBW) complex, which significantly enhances its transcriptional activation activity. A 163-bp InDel detected in the promoter region of the alleles of MdWRKY10 in a hybrid population of identical heterozygous genotypes regarding R6 by structural variation analysis, contains a typical W-box element that MdWRKY10 binds to for transactivation. This leads to increased transcript levels of MdWRKY10 and MdMYB10 and enhanced anthocyanin synthesis in the flesh, largely accounting for the various degrees of flesh red pigmentation in the R6 background. These findings reveal a novel regulatory role of the WRKY-containing protein complex in the formation of red flesh apple phenotypes and provide broader insights into the molecular mechanism governing anthocyanin synthesis in plants.

Transcriptome and weighted gene co-expression network analyses reveal key genes and pathways involved in early fruit ripening in Citrus sinensis

BACKGROUND

The fruit ripening period is an important target trait in fruit tree crop breeding programs. Thus, citrus tree breeders seek to develop extreme early ripening cultivars that allow optimization of citrus maturation periods. In this study, we explored the regulatory network involved in fruit ripening in Citrus sinensis using the 'Newhall' navel orange variety and its early-ripening mutant, 'Gannanzao'. This research will provide a basis for further research on important signaling pathways, gene functions and variety breeding of Citrus sinensis related to fruit ripening period.

RESULTS

Physiological analyses suggested that early fruit ripening in 'Gannanzao' is regulated by early accumulation of abscisic acid (ABA), persistently high levels of jasmonic acid (JA), and higher sucrose content in the pericarp. Pericarp samples from 'Gannanzao' and 'Newhall' navel oranges were sampled for RNA sequencing analysis at 180, 200, and 220 days after flowering; 1430 differentially expressed genes (DEGs) were identified. Functional enrichment analysis indicated that these DEGs were mainly enriched in the plant hormone signal transduction and sugar metabolism pathways, as well as other pathways related to fruit ripening. Important DEGs associated with fruit ripening in 'Gannanzao' included genes involved in ABA and JA metabolism and signal transduction, as well as sugar metabolism. Weighted gene co-expression network analysis showed that the deep pink module had the strongest correlations with ABA content, JA content, and early ripening. Based on gene functionality and gene expression analyses of 37 genes in this module, two candidate hub genes and two ethylene response factor 13 (ERF13) genes (Cs_ont_5g000690 and Cs_ont_5g000700) were identified as key genes regulated by ABA and JA signaling. These findings will help to clarify the mechanisms that underlie early citrus fruit ripening and will lead to the development of excellent genetic resources for further breeding of extreme early-ripening varieties.

CONCLUSIONS

Through analyses of the 'Newhall' navel orange cultivar and its early-ripening mutant 'Gannanzao', we identified genes involved in ABA and JA metabolism, signal transduction, and sugar metabolism that were related to fruit ripening. Among these, two ERF13 genes were inferred to be key genes in the regulation of fruit ripening. These findings provide insights into the genetic architecture related to early fruit ripening in C. sinensis.

Genetic and Biotechnological Approaches to Improve Fruit Bioactive Content: A Focus on Eggplant and Tomato Anthocyanins

Anthocyanins are a large group of water-soluble flavonoid pigments. These specialized metabolites are ubiquitous in the plant kingdom and play an essential role not only in plant reproduction and dispersal but also in responses to biotic and abiotic stresses. Anthocyanins are recognized as important health-promoting and chronic-disease-preventing components in the human diet. Therefore, interest in developing food crops with improved levels and compositions of these important nutraceuticals is growing. This review focuses on work conducted to elucidate the genetic control of the anthocyanin pathway and modulate anthocyanin content in eggplant (Solanum melongena L.) and tomato (Solanum lycopersicum L.), two solanaceous fruit vegetables of worldwide relevance. While anthocyanin levels in eggplant fruit have always been an important quality trait, anthocyanin-based, purple-fruited tomato cultivars are currently a novelty. As detailed in this review, this difference in the anthocyanin content of the cultivated germplasm has largely influenced genetic studies as well as breeding and transgenic approaches to improve the anthocyanin content/profile of these two important solanaceous crops. The information provided should be of help to researchers and breeders in devising strategies to address the increasing consumer demand for nutraceutical foods.

Genome-wide characterization of SmZHD gene family and the role of SmZHD12 in regulating anthocyanin biosynthesis in eggplant (Solanum melongena L.)

KEY MESSAGE

SmZHDs was highly expressed in anthocyanin-rich parts of eggplant. SmZHD12 can activate the expression of SmCHS, SmANS, SmDFR and SmF3H. Overexpression of SmZHD12 promotes anthocyanin biosynthesis in Arabidopsis. The Zinc finger-homeodomain (ZHD) proteins family genes are known to play a significant role in plant development and physiological processes. However, the evolutionary history and function of the ZHD gene family in eggplant remain largely unexplored. This study categorizes a total of 15 SmZHD genes into SmMIF and SmZHD subfamilies based on conserved domains. The phylogeny, gene structure, conserved motifs, promoter elements, and chromosomal locations of the SmZHD genes were comprehensively analyzed. Tissue expression profiles indicate that the majority of SmZHD genes are expressed in anthocyanin-rich areas. qRT-PCR assays revealed distinct expression patterns of SmZHD genes in response to various treatments, indicating their potential involvement in multiple signaling pathways. Analysis of transcriptomic data from light-treated eggplant peel identified SmZHD12 as the most light-responsive gene among the 15 SmZHD genes. Consequently, this study provides further evidence that SmZHD12 facilitates anthocyanin accumulation in Arabidopsis leaves by upregulating the expression of anthocyanin biosynthesis structural genes, as confirmed by dual-luciferase assays and Arabidopsis genetic transformation. Our study will lay a solid foundation for the in-depth study of the involvement of SmZHD genes in the regulation of anthocyanin biosynthesis.

Transcriptome analysis of green and purple fruited pepper provides insight into novel regulatory genes in anthocyanin biosynthesis

Background

Pepper (Capsicum annuum L.) is a valuable horticultural crop with economic significance, and its purple fruit color is attributed to anthocyanin, a phytonutrient known for its health-promoting benefits. However, the mechanisms regulating anthocyanin biosynthesis in pepper have yet to be fully elucidated.

Methods

RNA sequencing (RNA-seq) was utilized to analyze the transcriptome of fruits from three purple-fruited varieties (HN191, HN192, and HN005) and one green-fruited variety (EJT) at various developmental stages. To determine the relationships between samples, Pearson correlation coefficients (PCC) and principal component analysis (PCA) were calculated. Differential expression analysis was performed using the DESeq2 package to identify genes that were expressed differently between two samples. Transcription factors (TF) were predicted using the iTAK program. Heatmaps of selected genes were generated using Tbtools software.

Results

The unripe fruits of HN191, HN192, and HN005, at the stages of 10, 20, and 30 days after anthesis (DAA), display a purple color, whereas the unripe fruits of variety EJT remain green. To understand the molecular basis of this color difference, five transcriptome comparisons between green and purple fruits were conducted: HN191-10 vs EJT-10, HN191-20 vs EJT-20, HN191-30 vs EJT-30, HN192-30 vs EJT-30, and HN005-30 vs EJT-30. Through this analysis, 503 common differentially expressed genes (DEGs) were identified. Among these DEGs, eight structural genes related to the anthocyanin biosynthesis pathway and 24 transcription factors (TFs) were detected. Notably, one structural gene (MSTRG.12525) and three TFs (T459_25295, T459_06113, T459_26036) exhibited expression patterns that suggest they may be novel candidate genes involved in anthocyanin biosynthesis. These results provide new insights into the regulation of anthocyanin biosynthesis in purple pepper fruit and suggest potential candidate genes for future genetic improvement of pepper germplasm with enhanced anthocyanin accumulation.

Comprehensive dissection of variation and accumulation of free amino acids in tea accessions

Free amino acids (FAAs) positively determine the tea quality, notably theanine (Thea), endowing umami taste of tea infusion, which is the profoundly prevalent research in albino tea genetic resources. Therefore, 339 tea accessions were collected to study FAAs level for deciphering its variation and accumulation mechanism. Interestingly, alanine (Ala) and Thea which had the highest diversity index (H') value among three varieties of Camellia sinensis (L.) O. Kuntze were significantly higher than wild relatives (P < 0.05). The intraspecific arginine (Arg) and glutamine (Gln) contents in C. sinensis var. assamica were significantly lower than sinensis and pubilimba varieties. Moreover, the importance of interdependencies operating across FAAs and chlorophyll levels were highlighted via the cell ultrastructure, metabolomics, and transcriptome analysis. We then determined that the association between phytochrome interacting factor 1 (CsPIF1) identified by weighted gene co-expression network analysis (WGCNA) and Thea content. Intriguingly, transient knock-down CsPIF1 expression increased Thea content in tea plant, and the function verification of CsPIF1 in Arabidopsis also indicated that CsPIF1 acts as a negative regulator of Thea content by mainly effecting the genes expression related to Thea biosynthesis, transport, and hydrolysis, especially glutamate synthase (CsGOGAT), which was validated to be associated with Thea content with a nonsynonymous SNP by Kompetitive Allele-Specific PCR (KASP). We also investigated the interspecific and geographical distribution of this SNP. Taken together, these results help us to understand and clarify the variation and profile of major FAAs in tea germplasms and promote efficient utilization in tea genetic improvement and breeding.

RNA-sequencing analysis reveals novel genes involved in the different peel color formation in eggplant

Eggplant (Solanum melongena L.) is a highly nutritious vegetable. Here, the molecular mechanism of color formation in eggplants was determined using six eggplant cultivars with different peel colors and two SmMYB113-overexpressing transgenic eggplants with a purple peel and pulp. Significant differentially expressed genes (DEGs) were identified by RNA-sequencing analysis using the following criteria: log2(sample1/sample2) ≥ 0.75 and q-value ≤ 0.05. Two analytical strategies were used to identify genes related to the different peel color according to the peel color, flavonoids content, delphinidins/flavonoids ratio, and the content of anthocyanins. Finally, 27 novel genes were identified to be related to the color difference among eggplant peels and 32 novel genes were identified to be related to anthocyanin biosynthesis and regulated by SmMYB113. Venn analysis revealed that SmCytb5, SmGST, SmMATE, SmASAT3, and SmF3'5'M were shared among both sets of novel genes. Transient expression assay in tobacco suggested that these five genes were not sufficient for inducing anthocyanin biosynthesis alone, but they play important roles in anthocyanin accumulation in eggplant peels. Yeast one-hybrid, electrophoretic mobility shift assay and dual-luciferase assays indicated that the expression of the five genes could be directly activated by SmMYB113 protein. Finally, a regulatory model for the mechanism of color formation in eggplant was proposed. Overall, the results of this study provide useful information that enhances our understanding of the molecular mechanism underlying the different color formation in eggplant.

Molecular Regulatory Network of Anthocyanin Accumulation in Black Radish Skin as Revealed by Transcriptome and Metabonome Analysis

To understand the coloring mechanism in black radish, the integrated metabolome and transcriptome analyses of root skin from a black recombinant inbred line (RIL 1901) and a white RIL (RIL 1911) were carried out. A total of 172 flavonoids were detected, and the analysis results revealed that there were 12 flavonoid metabolites in radish root skin, including flavonols, flavones, and anthocyanins. The relative concentrations of most flavonoids in RIL 1901 were higher than those in RIL 1911. Meanwhile, the radish root skin also contained 16 types of anthocyanins, 12 of which were cyanidin and its derivatives, and the concentration of cyanidin 3-o-glucoside was very high at different development stages of black radish. Therefore, the accumulation of cyanidin and its derivatives resulted in the black root skin of radish. In addition, a module positively related to anthocyanin accumulation and candidate genes that regulate anthocyanin synthesis was identified by the weighted gene co-expression network analysis (WGCNA). Among them, structural genes (RsCHS, RsCHI, RsDFR, and RsUGT75C1) and transcription factors (TFs) (RsTT8, RsWRKY44L, RsMYB114, and RsMYB308L) may be crucial for the anthocyanin synthesis in the root skin of black radish. The anthocyanin biosynthesis pathway in the root skin of black radish was constructed based on the expression of genes related to flavonoid and anthocyanin biosynthesis pathways (Ko00941 and Ko00942) and the relative expressions of metabolites. In conclusion, this study not only casts new light on the synthesis and accumulation of anthocyanins in the root skin of black radish but also provides a molecular basis for accelerating the cultivation of new black radish varieties.

Metabolomic and transcriptomic analyses reveal the effects of grafting on blood orange quality

Introduction

Blood orange (Citrus sinensis L.) is a valuable source of nutrition because it is enriched in anthocyanins and has high organoleptic properties. Grafting is commonly used in citriculture and has crucial effects on various phenotypes of the blood orange, including its coloration, phenology, and biotic and abiotic resistance. Still, the underlying genetics and regulatory mechanisms are largely unexplored.

Methods

In this study, we investigated the phenotypic, metabolomic, and transcriptomic profiles at eight developmental stages of the lido blood orange cultivar (Citrus sinensis L. Osbeck cv. Lido) grafted onto two rootstocks.

Results and discussion

The Trifoliate orange rootstock provided the best fruit quality and flesh color for Lido blood orange. Comparative metabolomics suggested significant differences in accumulation patterns of metabolites and we identified 295 differentially accumulated metabolites. The major contributors were flavonoids, phenolic acids, lignans and coumarins, and terpenoids. Moreover, transcriptome profiling resulted in the identification of 4179 differentially expressed genes (DEGs), and 54 DEGs were associated with flavonoids and anthocyanins. Weighted gene co-expression network analysis identified major genes associated to 16 anthocyanins. Furthermore, seven transcription factors (C2H2, GANT, MYB-related, AP2/ERF, NAC, bZIP, and MYB) and five genes associated with anthocyanin synthesis pathway (CHS, F3H, UFGT, and ANS) were identified as key modulators of the anthocyanin content in lido blood orange. Overall, our results revealed the impact of rootstock on the global transcriptome and metabolome in relation to fruit quality in lido blood orange. The identified key genes and metabolites can be further utilized for the quality improvement of blood orange varieties.

Genome-wide identification and bioinformatics analysis of the WRKY transcription factors and screening of candidate genes for anthocyanin biosynthesis in azalea (Rhododendron simsii)

WRKY transcription factors have been demonstrated to influence the anthocyanin biosynthesis in many plant species. However, there is limited knowledge about the structure and function of WRKY genes in the major ornamental plant azalea (Rhododendron simsii). In this study, we identified 57 RsWRKY genes in the R. simsii genome and classified them into three main groups and several subgroups based on their structural and phylogenetic characteristics. Comparative genomic analysis suggested WRKY gene family has significantly expanded during plant evolution from lower to higher species. Gene duplication analysis indicated that the expansion of the RsWRKY gene family was primarily due to whole-genome duplication (WGD). Additionally, selective pressure analysis (Ka/Ks) suggested that all RsWRKY duplication gene pairs underwent purifying selection. Synteny analysis indicated that 63 and 24 pairs of RsWRKY genes were orthologous to Arabidopsis thaliana and Oryza sativa, respectively. Furthermore, RNA-seq data was used to investigate the expression patterns of RsWRKYs, revealing that 17 and 9 candidate genes may be associated with anthocyanin synthesis at the bud and full bloom stages, respectively. These findings provide valuable insights into the molecular mechanisms underlying anthocyanin biosynthesis in Rhododendron species and lay the foundation for future functional studies of WRKY genes.

Anthocyanins in metabolites of purple corn

Purple corn (Zea mays L.) is a special variety of corn, rich in a large amount of anthocyanins and other functional phytochemicals, and has always ranked high in the economic benefits of the corn industry. However, most studies on the stability of agronomic traits and the interaction between genotype and environment in cereal crops focus on yield. In order to further study the accumulation and stability of special anthocyanins in the growth process of purple corn, this review starts with the elucidation of anthocyanins in purple corn, the biosynthesis process and the gene regulation mechanism behind them, points out the influence of anthocyanin metabolism on anthocyanin metabolism, and introduces the influence of environmental factors on anthocyanin accumulation in detail, so as to promote the multi-field production of purple corn, encourage the development of color corn industry and provide new opportunities for corn breeders and growers.

CsMYC2 is involved in the regulation of phenylpropanoid biosynthesis induced by trypsin in cucumber (Cucumis sativus) during storage

Trypsin has a new activity of scavenging superoxide anion and generating hydrogen peroxide. Trypsin can significantly improve the storage quality of C. sativus. To illustrate the mechanism of trypsin-induced resistance in fruits and vegetables, an integrated analysis of widely targeted metabolomics and transcriptomics was carried out. Transcriptomic results showed that 1068 genes highly related to phenylpropanoid biosynthesis gathered in the brown module were obtained by WGCNA. In KEGG analysis, differentially expressed genes (DEGs) were also highly enriched in EIP (Environmental Information Processing) pathways "Plant hormone signal transduction (map04075)" and "MAPK signaling pathway-plant (map04016)". Next, 87 genes were identified as the leading edge by GSEA analysis. So far, CsMYC2 was highlighted as a key transcription factor that regulates phenylpropanoid biosynthesis identified by GSEA and WGCNA. Furthermore, the major route of biosynthesis of phenylpropanoid compounds including coumarins, lignins, chlorogenic acid, flavonoids, and derivatives regulated by trypsin was also illustrated by both transcriptomic and metabolomic data. Results of O2PLS showed that CsMYC2 was positively correlated with Rosmarinic acid-3-O-glucoside, Epigallocatechin, Quercetin-3-O-sophoroside (Baimaside), and so on. Correlation between CsMYC2, phenylpropanoid related genes, and metabolites in C. sativus was illustrated by co-expression networks. Roles of CsMYC2 were further checked in C. sativus by VIGS. The results of this study might give new insight into the exploration of the postharvest resistance mechanism of C. sativus induced by trypsin and provide useful information for the subsequent mining of resistance genes in C. sativus.

A Compendium for Novel Marker-Based Breeding Strategies in Eggplant

The worldwide production of eggplant is estimated at about 58 Mt, with China, India and Egypt being the major producing countries. Breeding efforts in the species have mainly focused on increasing productivity, abiotic and biotic tolerance/resistance, shelf-life, the content of health-promoting metabolites in the fruit rather than decreasing the content of anti-nutritional compounds in the fruit. From the literature, we collected information on mapping quantitative trait loci (QTLs) affecting eggplant's traits following a biparental or multi-parent approach as well as genome-wide association (GWA) studies. The positions of QTLs were lifted according to the eggplant reference line (v4.1) and more than 700 QTLs were identified, here organized into 180 quantitative genomic regions (QGRs). Our findings thus provide a tool to: (i) determine the best donor genotypes for specific traits; (ii) narrow down QTL regions affecting a trait by combining information from different populations; (iii) pinpoint potential candidate genes.

Combined widely targeted metabolomics and transcriptomics analysis reveals differentially accumulated metabolites and the underlying molecular bases in fleshy taproots of distinct radish genotypes

Radish is an important taproot crop with medicinal and edible uses that is cultivated worldwide. However, the differences in metabolites and the underlying molecular bases among different radish types remain largely unknown. In the present study, a combined analysis of liquid chromatography-electrospray ionization-tandem mass spectrometry (LC-ESI-MS/MS) and RNA-Seq data was conducted to uncover important differentially accumulated metabolites (DAMs) among radish accessions with green, white and red taproot flesh colours. A total of 657 metabolites were identified and 138 DAMs were commonly present in red vs. green and red vs. white accessions. Red accessions were rich in phenolic compounds, while green and white accessions had more amino acids. Additionally, 41 metabolites and 98 genes encoding 37 enzymes were enriched in the shikimate and phenolic biosynthesis pathways. CHS is the rate-limiting enzyme determining flavonoid differences among accessions. A total of 119 candidate genes might contribute to red accession-specific accumulated metabolites. Specifically, one gene cluster consisting of 16 genes, including one RsMYB1, RsMYBL2, RsTT8, RsDFR, RsANS, Rs4CL3, RsSCPL10, Rs3AT1 and RsSAP2 gene, two RsTT19 and RsWRKY44 genes and three RsUGT genes, might be involved in anthocyanin accumulation in red radish fleshy taproots. More importantly, an InDel marker was developed based on an RsMYB1 promoter mutation, and the accuracy reached 95.9% when it was used to select red-fleshed radishes. This study provides comprehensive insights into the metabolite differences and underlying molecular mechanisms in fleshy taproots among different radish genotypes and will be beneficial for the genetic improvement of radish nutritional quality.

Integrated Metabolome and Transcriptome Analysis Reveals a Regulatory Network of Fruit Peel Pigmentation in Eggplant (Solanum melongena L.)

The color of fruit peel is an economically important character of eggplant, and black-purple eggplant has received much attention for being rich in anthocyanin. However, the reason why different fruit peel colors form in eggplant is not well understood. In the present study, an integrative analysis of the metabolome and transcriptome profiles was performed in five eggplant varieties with different fruit colors. A total of 260 flavonoids were identified, and most of them showed significantly higher abundance in black-purple varieties than in other varieties. The transcriptome analysis indicated the activation of early phenylpropanoid biosynthesis genes (SmPAL, SmC4H, and Sm4CL) was more responsible for anthocyanin accumulation, while SmF3′5′H was the key factor for the formation of a purple color. Furthermore, two transcription factors, SmGL2 and SmGATA26, were identified as new hub genes associated with anthocyanin accumulation. The silencing of SmGL2 and SmGATA26 reduced anthocyanin accumulation in eggplant fruit peels, suggesting the possible involvement of SmGL2 and SmGATA26 in regulating anthocyanin biosynthesis. In addition, the pathway of plant hormone signal transduction was significantly enriched, indicating that phytohormones may cooperatively interact to modulate flavonoid biosynthesis. This study provides comprehensive information of flavonoid metabolites and new insights into the regulatory network of fruit coloration, which might be useful for the molecular breeding of eggplant.

Transcriptome analysis and mining of genes related to shade tolerance in foxtail millet (Setaria italica (L.) P. Beauv.)

A stereo interplanting system with foxtail millet beneath chestnut trees is an effective planting method to raise the utilization of land in chestnut orchards, increase yields and improve quality of chestnut nuts. Consequently, exploration of genes involved in shade tolerance response in foxtail millet and breeding shade-tolerant varieties have become urgent issues. In this study, RNA-seq of leaf samples from two shade-tolerant varieties and three shade-intolerant varieties of foxtail millet at the booting stage was performed. Comparisons between the varieties revealed that 70 genes were commonly differentially expressed. Moreover, the ratio of net photosynthetic rate under shaded environment to that under light environment could be used as an indicator of shade tolerance. Subsequently, weighted gene co-expression network analysis was employed to construct a co-expression network and modules were correlated with this ratio. A total of 375 genes were identified as potentially relevant to shade tolerance, among which nine genes were also present in the 70 differentially expressed genes, which implied that they were good candidates for genes involved in shade tolerance. Our results provide valuable resources for elucidation of the molecular mechanisms underlying shade tolerance and will contribute to breeding of shade-tolerant foxtail millet that are adapted to the shaded environment under chestnut trees.

Mining Candidate Genes Related to Heavy Metals in Mature Melon (Cucumis melo L.) Peel and Pulp Using WGCNA

The content of metal ions in fruits is inseparable from plant intake of trace elements and health effects in the human body. To understand metal ion content in the fruit and pericarp of melon (Cucumis melo L.) and the candidate genes responsible for controlling this process, we analyzed the metal ion content in distinct parts of melon fruit and pericarp and performed RNA-seq. The results showed that the content of metal ions in melon fruit tissue was significantly higher than that in the pericarp. Based on transcriptome expression profiling, we found that the fruit and pericarp contained elevated levels of DEGs. GO functional annotations included cell surface receptor signaling, signal transduction, organic substance metabolism, carbohydrate derivative binding, and hormone-mediated signaling pathways. KEGG pathways included pectate lyase, pentose and glucuronate interconversions, H⁺-transporting ATPase, oxidative phosphorylation, plant hormone signal transduction, and MAPK signaling pathways. We also analyzed the expression patterns of genes and transcription factors involved in hormone biosynthesis and signal transduction. Using weighted gene co-expression network analysis (WGCNA), a co-expression network was constructed to identify a specific module that was significantly correlated with the content of metal ions in melon, after which the gene expression in the module was measured. Connectivity and qRT–PCR identified five candidate melon genes, LOC103501427, LOC103501539, LOC103503694, LOC103504124, and LOC107990281, associated with metal ion content. This study provides a theoretical basis for further understanding the molecular mechanism of heavy metal ion content in melon fruit and peel and provides new genetic resources for the study of heavy metal ion content in plant tissues.

Response of Anthocyanin Accumulation in Pepper (Capsicum annuum) Fruit to Light Days

Light is the key factor affecting the synthesis of anthocyanins in pepper. In this study, pepper fruit under different light days was used as experimental material to explore the synthesis of anthocyanins in purple pepper. A total of 38 flavonoid metabolites were identified in the purple pepper germplasm HNUCA21 by liquid chromatography–tandem mass spectrometry (LC-MS/MS), of which 30 belong to anthocyanins. The detected anthocyanin with the highest content was Delphinidin-3-O-glucoside (17.13 µg/g), which reached the maximum after 168 h of light treatment. Through weighted gene co-expression network analysis (WGCNA), the brown module was identified to be related to the early synthesis of anthocyanins. This module contains many structural genes related to flavonoid synthesis, including chalcone synthase (CHS 107871256, 107864266), chalcone isomerase (CHI 107871144, 107852750), dihydroflavonol 4-reductase (DFR 107860031), flavonoid 3′ 5′-hydroxylase (F3’5’H 107848667), flavonoid 3′-monooxygenase (F3M 107862334), leucoanthocyanidin dioxygenase (LDOX 107866341), and trans-cinnamate 4-monooxygenase (TCM 107875406, 107875407). The module also contained some genes related to anthocyanin transport function, such as glutathione S-transferase (GST 107861273), anthocyanidin 3-O-glucosyltransferase (UDPGT 107861697, 107843659), and MATE (107863234, 107844661), as well as some transcription factors, such as EGL1 (107865400), basic helix-loop-helix 104 (bHLH104 107864591), and WRKY44 (107843538, 107843524). The co-expression regulatory network indicated the involvement of CHS, DFR, CHI, and EGL1, as well as two MATE and two WRKY44 genes in anthocyanin synthesis. The identified genes involved in early, middle, and late light response provided a reference for the further analysis of the regulatory mechanism of anthocyanin biosynthesis in pepper.

Transcriptomic and Metabolomic Analysis of Quality Changes during Sweet Cherry Fruit Development and Mining of Related Genes

Sweet cherries are economically important fruit trees, and their quality changes during development need to be determined. The mechanism of fruit quality changes in sweet cherries were determined by analyzing sweet cherry fruits at 12 developmental stages. The results showed that the soluble sugar, anthocyanin content, and hormones of sweet cherries all changed drastically during the color transition. Therefore, the fruits at the beginning of color conversion, at the end of color conversion, and at the ripening state were selected for the comprehensive analysis of their metabolome and transcriptome. Different sugars, such as D-glucose, sucrose, and trehalose, were identified in the metabolome. Dihydroquercetin, delphinidin-3-glucoside, cyanidin-3-rutincoside, and other flavonoid species were also identified. D-glucose and cyanidin-3-rutinoside were among the most important components of sweet cherry soluble sugars and anthocyanins, respectively. The transcriptional analysis identified key structural genes and nine transcription factors involved in the ABA, sugar, organic acid, and anthocyanin synthesis pathways, with the following specific regulatory patterns. NAC71, WRKY57, and WRKY3 regulate fruit sugar accumulation mainly by acting on INV, SPS, and SUS. MYC2 is involved in the synthesis of anthocyanin precursors by activating PAL and C4H, whereas TCP7 mainly regulates CHI and F3H. WRKY3, NAC71, and WRKY57 have important positive regulatory significance on anthocyanin accumulation, mainly by activating the expression of DFR, ANS, and 3GT.

Combining QTL mapping and gene co-expression network analysis for prediction of candidate genes and molecular network related to yield in wheat

BACKGROUND

Wheat (Triticum aestivum L.) is an important cereal crop. Increasing grain yield for wheat is always a priority. Due to the complex genome of hexaploid wheat with 21 chromosomes, it is difficult to identify underlying genes by traditional genetic approach. The combination of genetics and omics analysis has displayed the powerful capability to identify candidate genes for major quantitative trait loci (QTLs), but such studies have rarely been carried out in wheat. In this study, candidate genes related to yield were predicted by a combined use of linkage mapping and weighted gene co-expression network analysis (WGCNA) in a recombinant inbred line population.

RESULTS

QTL mapping was performed for plant height (PH), spike length (SL) and seed traits. A total of 68 QTLs were identified for them, among which, 12 QTLs were stably identified across different environments. Using RNA sequencing, we scanned the 99,168 genes expression patterns of the whole spike for the recombinant inbred line population. By the combined use of QTL mapping and WGCNA, 29, 47, 20, 26, 54, 46 and 22 candidate genes were predicted for PH, SL, kernel length (KL), kernel width, thousand kernel weight, seed dormancy, and seed vigor, respectively. Candidate genes for different traits had distinct preferences. The known PH regulation genes Rht-B and Rht-D, and the known seed dormancy regulation genes TaMFT can be selected as candidate gene. Moreover, further experiment revealed that there was a SL regulatory QTL located in an interval of about 7 Mbp on chromosome 7A, named TaSL1, which also involved in the regulation of KL.

CONCLUSIONS

A combination of QTL mapping and WGCNA was applied to predicted wheat candidate genes for PH, SL and seed traits. This strategy will facilitate the identification of candidate genes for related QTLs in wheat. In addition, the QTL TaSL1 that had multi-effect regulation of KL and SL was identified, which can be used for wheat improvement. These results provided valuable molecular marker and gene information for fine mapping and cloning of the yield-related trait loci in the future.

Functional differentiation of two lhx8 paralogs and possible regulatory role of lhx8a in Japanese flounder (Paralichthys olivaceus)

Lhx8, belonging to the LIM-Homebox family, is involved in the tooth, nervous system, and primordial follicles development in mammals. However, little is known about the regulatory roles of lhx8 in teleosts. In this study, two lhx8 duplicates were identified in Paralichthys olivaceus, termed Polhx8a and Polhx8b, respectively. Bioinformatic analysis showed that Polhx8a was more likely to be a teleost-specific paralog. According to expression analysis, Polhx8a transcripts were almost exclusively concentrated in the oocytes, while Polhx8b was weakly expressed in the spleen, gill, and some facial organs, indicating sub-functionalization of this gene pair during evolution. Furthermore, Polhx8a mRNA level elevated from perinucleolar oocyte (PNO) stage to vitellogenic oocyte (VO) stage transition and changed after exogenous hormone stimulation, proving that Polhx8a was involved in the oocyte development and could be regulated by sex hormones. Yeast two-hybrid, bimolecular fluorescence complementation (BiFC) and co-immunoprecipitation (co-IP) experiments captured the positive protein interactions between PoLhx8a and the other two oocyte-specific transcription factors: PoFigla and PoNobox. After knocking down lhx8a in embryos or adult ovaries in vivo, the expression of oocyte-associated genes was significantly down-regulated (P < 0.05). Our findings suggest the evolution and functional differentiation of lhx8 genes, and shed light on the potential role of lhx8a in protein interactions and gene regulation in teleosts.

Identification of key gene networks controlling anthocyanin biosynthesis in peach flower

Flavonoids, particularly anthocyanin is the main pigment that determined the red color of peach flowers, and help the plant to attract pollinators, protect the reproductive organs of flower from photo-oxidative effects of light and various non-communicable diseases. Through weightage gene coexpression network analysis (WGCNA) we identified a network of 15 hub genes that co-expressed throughout peach flower development including 5 genes coded for the key enzymes (CHI, F3'H, DFR, LAR and UFGT) of flavonoid biosynthetic pathway and 1 gene Prupe.1G111700 identified as R2R3 family transcription factor MYB108. Over expression of PpMYB108 significantly increased anthocyanin biosynthesis in Tobacco flowers. Moreover, the expression correlation between PpMYB108 and PpDFR, suggests that PpMYB108 play the role of transcriptional activator for PpDFR. This was further supported by a 6 bp insertion of MYB biding site in the core promoter region of PpDFR in red flower. The positive interaction of PpMYB108 with PpDFR promoter from red flower was confirmed in yeast one hybrid assay. These findings may be helpful in peach breeding programs as well as in identifying anthocyanin related genes in other species.

Unraveling the Mechanism of Purple Leaf Formation in Brassica napus by Integrated Metabolome and Transcriptome Analyses

Brassica napus as both oilseed and vegetable, is widely cultivated in China. The purple leaf of B. napus is rich in anthocyanins and can provide valuable nutrients. Although several high-anthocyanin cultivars have been reported, the molecular mechanism underlying anthocyanin biosynthesis in B. napus remains lesser-known. Therefore, in this study, we conducted integrative metabolome and transcriptome analyses in three B. napus cultivars with different leaf colors. Overall, 39 flavonoids were identified (including 35 anthocyanins), and 22 anthocyanins were differentially accumulated in the leaves, contributing to the different leaf colors. Cyanidin-3,5,3'-O-triglucoside was confirmed as the main contributor of the purple leaf phenotype. Meanwhile, other anthocyanins may play important roles in deepening the color of B. napus leaves. A total of 5,069 differentially expressed genes (DEGs) and 32 overlapping DEGs were identified by RNA-sequencing; hence, the correlation between anthocyanin content and DEG expression levels was explored. Two structural genes (DFR and ANS), three GSTs (homologous to TT19), and 68 differentially expressed transcription factors (TFs), especially MYB-related TFs and WRKY44, were identified in three B. napus varieties characterized by different leaf color, thereby indicating that these genes may contribute to anthocyanin biosynthesis, transport, or accumulation in B. napus leaves. The findings of study provide important insights that may contribute to gaining a better understanding of the transcriptional regulation of anthocyanin metabolism in B. napus.

Comparative transcriptomic analysis reveals key components controlling spathe color in Anthurium andraeanum (Hort.)

Anthurium andraeanum (Hort.) is an important ornamental in the tropical cut-flower industry. However, there is currently insufficient information to establish a clear connection between the genetic model(s) proposed and the putative genes involved in the differentiation between colors. In this study, 18 cDNA libraries related to the spathe color and developmental stages of A. andraeanum were characterized by transcriptome sequencing (RNA-seq). For the de novo transcriptome, a total of 114,334,082 primary sequence reads were obtained from the Illumina sequencer and were assembled into 151,652 unigenes. Approximately 58,476 transcripts were generated and used for comparative transcriptome analysis between three cultivars that differ in spathe color (‘Sasha’ (white), ‘Honduras’ (red), and ‘Rapido’ (purple)). A large number of differentially expressed genes (8,324), potentially involved in multiple biological and metabolic pathways, were identified, including genes in the flavonoid and anthocyanin biosynthetic pathways. Our results showed that the chalcone isomerase (CHI) gene presented the strongest evidence for an association with differences in color and the highest correlation with other key genes (flavanone 3-hydroxylase (F3H), flavonoid 3’5’ hydroxylase (F3’5’H)/ flavonoid 3’-hydroxylase (F3’H), and leucoanthocyanidin dioxygenase (LDOX)) in the anthocyanin pathway. We also identified a differentially expressed cytochrome P450 gene in the late developmental stage of the purple spathe that appeared to determine the difference between the red- and purple-colored spathes. Furthermore, transcription factors related to putative MYB-domain protein that may control anthocyanin pathway were identified through a weighted gene co-expression network analysis (WGCNA). The results provided basic sequence information for future research on spathe color, which have important implications for this ornamental breeding strategies.