Comparative phenotype and microRNAome in developing anthers of wild-type and male-sterile Lycium barbarum L

Integrative transcriptome and metabolome analysis reveals the discrepancy in the accumulation of active ingredients between Lycium barbarum cultivars

MAIN CONCLUSION

The combined analysis of transcriptome and metabolome provided molecular insight into the dynamics of multiple active ingredients biosynthesis and accumulation across different cultivars of Lycium barbarum. Lycium barbarum L. has a high concentration of active ingredients and is well known in traditional Chinese herbal medicine for its therapeutic properties. However, there are many Lycium barbarum cultivars, and the content of active components varies, resulting in inconsistent quality between Lycium barbarum cultivars. At present, few research has been conducted to reveal the difference in active ingredient content among different cultivars of Lycium barbarum at the molecular level. Therefore, the transcriptome of 'Ningqi No.1' and 'Qixin No.1' during the three development stages (G, T, and M) was constructed in this study. A total of 797,570,278 clean reads were obtained. Between the two types of wolfberries, a total of 469, 2394, and 1531 differentially expressed genes (DEGs) were obtained in the 'G1 vs. G10,' 'T1 vs. T10,' and 'M1 vs. M10,' respectively, and were annotated with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) orthology identifiers. Using these transcriptome data, most DEGs related to the metabolism of the active ingredients in 'Ningqi No.1' and 'Qixin No.1' were identified. Moreover, a widely targeted metabolome analysis of the metabolites of 'Ningqi 1' and 'Qixin 1' fruits at the maturity stage revealed 1,135 differentially expressed metabolites (DEMs) in 'M1 vs. M10,' and many DEMs were associated with active ingredients such as flavonoids, alkaloids, terpenoids, and so on. We further quantified the flavonoid, lignin, and carotenoid contents of the two Lycium barbarum cultivars during the three developmental stages. The present outcome provided molecular insight into the dynamics of multiple active ingredients biosynthesis and accumulation across different cultivars of Lycium barbarum, which would provide the basic data for the formation of Lycium barbarum fruit quality and the breeding of outstanding strains.

Transcriptional deciphering of the metabolic pathways associated with the bioactive ingredients of wolfberry species with different quality characteristics

BACKGROUND

Wolfberry is rich in carotenoids, flavonoids, vitamins, alkaloids, betaines and other bioactive ingredients. For over 2,000 years, wolfberry has been used in China as a medicinal and edible plant resource. Nevertheless, the content of bioactive ingredients varies by cultivars, resulting in uneven quality across wolfberry cultivars and species. To date, research has revealed little about the underlying molecular mechanism of the metabolism of flavonoids, carotenoids, and other bioactive ingredients in wolfberry.

RESULTS

In this context, the transcriptomes of the Lycium barbarum L. cultivar 'Ningqi No. 1' and Lycium chinense Miller were compared during the fruit maturity stage using the Illumina NovaSeq 6000 sequencing platform, and subsequently, the changes of the gene expression profiles in two types of wolfberries were analysed. In total, 256,228,924 clean reads were obtained, and 8817 differentially expressed genes (DEGs) were identified, then assembled by Basic Local Alignment Search Tool (BLAST) similarity searches and annotated using Gene Ontology (GO), Clusters of Orthologous Groups of proteins (KOG), and the Kyoto Encyclopedia of Genes and Genomes (KEGG). By combining these transcriptome data with data from the PubMed database, 36 DEGs related to the metabolism of bioactive ingredients and implicated in the metabolic pathway of carotenoids, flavonoids, terpenoids, alkaloids, vitamins, etc., were identified. In addition, among the 9 differentially expressed transcription factors, LbAPL, LbPHL11 and LbKAN4 have raised concerns. The protein physicochemical properties, structure prediction and phylogenetic analysis indicated that LbAPL and LbPHL11 may be good candidate genes involved in regulating the flavonoid metabolism pathway in wolfberry.

CONCLUSIONS

This study provides preliminary evidence for the differences in bioactive ingredient content at the transcription level among different wolfberry species, as well as a research and theoretical basis for the screening, cloning and functional analysis of key genes involved in the metabolism of bioactive ingredients in wolfberry.

Changes in m6A RNA methylation are associated with male sterility in wolfberry

BACKGROUND

N6-methyladenosine (m6A) modification is the most abundant type of RNA modification in eukaryotic cells, playing pivotal roles in multiple plant growth and development processes. Yet the potential role of m6A in conferring the trait of male sterility in plants remains unknown.

RESULTS

In this study, we performed RNA-sequencing (RNA-Seq) and m6A-sequencing (m6A-Seq) of RNAs obtained from the anther tissue of two wolfberry lines: 'Ningqi No.1' (LB1) and its natural male sterile mutant 'Ningqi No.5' (LB5). Based on the newly assembled transcriptome, we established transcriptome-wide m6A maps for LB1 and LB5 at the single nucleus pollen stage. We found that the gene XLOC_021201, a homolog of m6A eraser-related gene ALKBH10 in Arabidopsis thaliana, was significantly differentially expressed between LB1 and LB5. We also identified 1642 and 563 m6A-modified genes with hypermethylated and hypomethylated patterns, respectively, in LB1 compared with LB5. We found the hypermethylated genes significantly enriched in biological processes related to energy metabolism and lipid metabolism, while hypomethylation genes were mainly linked to cell cycle process, gametophyte development, and reproductive process. Among these 2205 differentially m6A methylated genes, 13.74% (303 of 2205) were differentially expressed in LB1 vis-à-vis LB5.

CONCLUSIONS

This study constructs the first m6A transcriptome map of wolfberry and establishes an association between m6A and the trait of male sterility in wolfberry.

Different MicroRNA Families Involved in Regulating High Temperature Stress Response during Cotton (Gossypium hirsutum L.) Anther Development

MicroRNAs (miRNAs) are small molecule RNAs widely involved in responses to plant abiotic stresses. We performed small RNA sequencing of cotton anthers at four developmental stages under normal and high temperature (NT and HT, respectively) conditions to investigate the stress response characteristics of miRNA to HT. A total of 77 miRNAs, including 33 known miRNAs and 44 novel miRNAs, were identified, and 41 and 28 miRNAs were differentially expressed under NT and HT stress conditions, respectively. The sporogenous cell proliferation (SCP), meiotic phase (MP), microspore release period (MRP), and pollen maturity (PM) stages had 10 (including 12 miRNAs), four (including six miRNAs), four (including five miRNAs), and seven (including 11 miRNAs) HT stress-responsive miRNA families, respectively, which were identified after removing the changes in genotype-specific miRNAs under NT condition. Seven miRNA families (miR2949, miR167, and miR160 at the SCP stage; miR156 and miR172 at the MP stage; miR156 at the MRP stage; and miR393 and miR3476 at the PM stage), which had expression abundance of more than 10% of the total expression abundance, served as the main regulators responding to HT stress with positive or negative regulation patterns. These miRNAs orchestrated the expression of the corresponding target genes and led to different responses in the HT-tolerant and the HT-sensitive lines. The results revealed that the HT stress response of miRNAs in cotton anthers were stage-specific and differed with the development of anthers. Our study may enhance the understanding of the response of miRNAs to HT stress in cotton anthers and may clarify the mechanism of plant tolerance to HT stress.

A SNP-Based High-Density Genetic Map of Leaf and Fruit Related Quantitative Trait Loci in Wolfberry (Lycium Linn.)

Wolfberry (Lycium Linn. 2n = 24) fruit, Gouqizi, is a perennial shrub, traditional food and medicinal plant resource in China. Leaf and fruit related characteristics are economically important traits that are the focus for genetic improvement, but few studies into the molecular genetics of this crop have been reported to date. Here, an F₁ population (302 individuals) derived from a cross between "NO.1 Ningqi" (Lycium barbarum L.) and "Chinese gouqi" (Lycium chinese Mill.) was constructed. We recorded fruit weight, longitude, diameter and index along with leaf length, width and index for three consecutive years from 2015 to 2017. Based on this population and these phenotypic data, we constructed the first high-density genetic map of Lycium using specific length amplified fragment sequencing (SLAF-seq) and analyzed quantitative trait loci (QTLs). The map contains 6733 single nucleotide polymorphisms and 12 linkage groups (LG) with a total map distance of 1702.45 cM and an average map distance of 0.253 cM. A total of 55 QTLs were mapped for more than 2 years, of which 18 stable QTLs for fruit index on LG 11, spanning an interval of 73.492-90.945 cM, were detected. qFI11-15 for fruit index was an impressive QTL with logarithm of odds (LOD) and proportion of variance explained (PEV) values reaching 11.07 and 19.7%, respectively. The QTLs on LG 11 were gathered tightly, having an average interval of less than 1 cM per QTL, suggesting that there might be a cluster region controlling fruit index. Remarkably, qLI10-2 and qLI11-2 for leaf index were detectable for 3 years. These results give novel insight into the genetic control of leaf and fruit related traits in Lycium and provide robust support for undertaking further positional cloning studies and implementing marker-assisted selection in seedlings.