Dynamic and epistatic QTL mapping reveals the complex genetic architecture of waterlogging tolerance in chrysanthemum

Comparative transcriptome analysis provides molecular insights into heterosis of waterlogging tolerance in Chrysanthemum indicum

BACKGROUND

Heterosis breeding is one of the most important breeding methods for chrysanthemum. To date, the genetic mechanisms of heterosis for waterlogging tolerance in chrysanthemum are still unclear. This study aims to analyze the expression profiles and potential heterosis-related genes of two hybrid lines and their parents with extreme differences in waterlogging tolerance under control and waterlogging stress conditions by RNA-seq.

RESULTS

A population of 140 F1 progeny derived from Chrysanthemum indicum (Nanchang) (waterlogging-tolerant) and Chrysanthemum indicum (Nanjing) (waterlogging-sensitive) was used to characterize the extent of genetic variation in terms of seven waterlogging tolerance-related traits across two years. Lines 98 and 95, respectively displaying positive and negative overdominance heterosis for the waterlogging tolerance traits together with their parents under control and waterlogging stress conditions, were used for RNA-seq. In consequence, the maximal number of differentially expressed genes (DEGs) occurred in line 98. Gene ontology (GO) enrichment analysis revealed multiple stress-related biological processes for the common up-regulated genes. Line 98 had a significant increase in non-additive genes under waterlogging stress, with transgressive up-regulation and paternal-expression dominant patterns being the major gene expression profiles. Further, GO analysis identified 55 and 95 transgressive up-regulation genes that overlapped with the up-regulated genes shared by two parents in terms of responses to stress and stimulus, respectively. 6,640 genes in total displaying maternal-expression dominance patterns were observed in line 95. In addition, 16 key candidate genes, including SAP12, DOX1, and ERF017 which might be of significant importance for the formation of waterlogging tolerance heterosis in line 98, were highlighted.

CONCLUSION

The current study provides a comprehensive overview of the root transcriptomes among F1 hybrids and their parents under waterlogging stress. These findings lay the foundation for further studies on molecular mechanisms underlying chrysanthemum heterosis on waterlogging tolerance.

Multi-locus genome-wide association studies reveal the dynamic genetic architecture of flowering time in chrysanthemum

KEY MESSAGE

The dynamic genetic architecture of flowering time in chrysanthemum was elucidated by GWAS. Thirty-six known genes and 14 candidate genes were identified around the stable QTNs and QEIs, among which ERF-1 was highlighted. Flowering time (FT) adaptation is one of the major breeding goals in chrysanthemum, a multipurpose ornamental plant. In order to reveal the dynamic genetic architecture of FT in chrysanthemum, phenotype investigation of ten FT-related traits was conducted on 169 entries in 2 environments. The broad-sense heritability of five non-conditional FT traits, i.e., budding (FBD), visible coloring (VC), early opening (EO), full-bloom (OF) and decay period (DP), ranged from 56.93 to 84.26%, which were higher than that of the five derived conditional FT traits (38.51-75.13%). The phenotypic variation coefficients of OF_EO and DP_OF were relatively large ranging from 30.59 to 36.17%. Based on 375,865 SNPs, the compressed variance component mixed linear model 3VmrMLM was applied for a multi-locus genome-wide association study (GWAS). As a result, 313 quantitative trait nucleotides (QTNs) were identified for the non-conditional FT traits in single-environment analysis, while 119 QTNs and 67 QTN-by-environment interactions (QEIs) were identified in multi-environment analysis. As for the conditional traits, 343 QTNs were detected in single-environment analysis, and 119 QTNs and 83 QEIs were identified in multi- environment analysis. Among the genes around stable QTNs and QEIs, 36 were orthologs of known FT genes in Arabidopsis and other plants; 14 candidates were mined by combining the transcriptomics data and functional annotation, including ERF-1, ACA10, and FOP1. Furthermore, the haplotype analysis of ERF-1 revealed six elite accessions with extreme FBD. Our findings contribute to the understanding of dynamic genetic architecture of FT and provide valuable resources for future chrysanthemum molecular breeding programs.

Expression dynamics of metabolites in diploid and triploid watermelon in response to flooding

Watermelon (Citrullus lanatus) is an economically important cucurbitaceous crop worldwide. The productivity of watermelon is affected by both biotic and abiotic stresses. Flooding has significant impacts on the growth of watermelons by causing oxygen deficiency and a loss of agricultural productivity. Currently, we used the triploid and diploid watermelon Zhengzhou No.3 to study the dynamics of metabolites in response to flooding stress. Quantification of metabolites was performed by UPLC-ESI-MS/MS at different time intervals i.e., 0, 3, 5 and 7 days under flooding stress. We observed that the activities of oxidants were higher in the diploid watermelon, whereas the higher antioxidant activities in the triploid watermelon makes them more resistant to the flooding stress. We also observed that the root activity and the chlorophyll in the triploid watermelon plants were higher as compared to the diploid watermelon plants. Co-expression network analysis leads to the identification of twenty-four hub metabolites that might be the key metabolites linked to flooding tolerance. Resolving the underlying mechanisms for flooding tolerance and identification of key molecules serving as indicators for breeding criteria are necessary for developing flooding-resistant varieties.

Transcriptome Analysis Reveals Genes Respond to Chlorophyll Deficiency in Green and Yellow Leaves of Chrysanthemum morifolium Ramat

Chlorophyll is vital for photosynthesis to produce sugars and other useful biochemical products in green plants. However, the molecular effects of chlorophyll deficiency in Chrysanthemum are largely unknown. In this study, we identified a bud sport mutant chrysanthemum belonging to the variety ‘Nannong Binyun’, which has yellow branches. Plant physiological studies have shown that the yellow color is revealed due to chlorophyll loss. RNA extracts of yellow and green tissues were analyzed using high-throughput RNA-sequencing, and a total of 11,649 tissue enriched unigenes that respond to chlorophyll deficiency were identified, including 4803 unigenes upregulated in yellow tissues and 6846 unigenes in green tissues. GO analysis revealed that these tissue-enriched genes may involve in the physiological processes of chlorophyll accumulation and photosynthesis. In addition, many DEGs from the families of AP2-EREBP, bHLH, MYB, and FAR1 that are associated with plant development and stress response were detected. Our study found that most of the genes from the GRAS family were downregulated in yellow leaves, indicating their putative roles in stem cell maintenance and possible contribution to leaf size determination.

Genome-wide association study and transcriptome analysis reveal key genes affecting root growth dynamics in rapeseed

BACKGROUND

In terms of global demand, rapeseed is the third-largest oilseed crop after soybeans and palm, which produces vegetable oil for human consumption and biofuel for industrial production. Roots are vital organs for plant to absorb water and attain mineral nutrients, thus they are of great importance to plant productivity. However, the genetic mechanisms regulating root development in rapeseed remain unclear. In the present study, seven root-related traits and shoot biomass traits in 280 Brassica napus accessions at five continuous vegetative stages were measured to establish the genetic basis of root growth in rapeseed.

RESULTS

The persistent and stage-specific genetic mechanisms were revealed by root dynamic analysis. Sixteen persistent and 32 stage-specific quantitative trait loci (QTL) clusters were identified through genome-wide association study (GWAS). Root samples with contrasting (slow and fast) growth rates throughout the investigated stages and those with obvious stage-specific changes in growth rates were subjected to transcriptome analysis. A total of 367 differentially expressed genes (DEGs) with persistent differential expressions throughout root development were identified, and these DEGs were significantly enriched in GO terms, such as energy metabolism and response to biotic or abiotic stress. Totally, 485 stage-specific DEGs with different expressions at specific stage were identified, and these DEGs were enriched in GO terms, such as nitrogen metabolism. Four candidate genes were identified as key persistent genetic factors and eight as stage-specific ones by integrating GWAS, weighted gene co-expression network analysis (WGCNA), and differential expression analysis. These candidate genes were speculated to regulate root system development, and they were less than 100 kb away from peak SNPs of QTL clusters. The homologs of three genes (BnaA03g52990D, BnaA06g37280D, and BnaA09g07580D) out of 12 candidate genes have been reported to regulate root development in previous studies.

CONCLUSIONS

Sixteen QTL clusters and four candidate genes controlling persistently root development, and 32 QTL clusters and eight candidate genes stage-specifically regulating root growth in rapeseed were detected in this study. Our results provide new insights into the temporal genetic mechanisms of root growth by identifying key candidate QTL/genes in rapeseed.

Linkage Map Development by EST-SSR Markers and QTL Analysis for Inflorescence and Leaf Traits in Chrysanthemum (Chrysanthemum morifolium Ramat.)

Chrysanthemums (Chrysanthemum morifolium Ramat.) are famous ornamental crops with high medicinal and industrial values. The inflorescence and leaf traits are key factors that affect the yield and quality of chrysanthemum. However, the genetic improvement of those traits is slow within chrysanthemum because of its hexaploidy, high heterozygosity and enormous genome. To study the genetic control of the important traits and facilitate marker-assisted selection (MAS) in chrysanthemum, it is desirable to populate the genetic maps with an abundance of transferrable markers such as microsatellites (SSRs). A genetic map was constructed with expressed sequence tag–simple sequence repeat (EST-SSR) markers in an F₁ progeny of 192 offspring. A total of 1000 alleles were generated from 223 EST-SSR primer pairs. The preliminary maternal and paternal maps consisted of 265 marker alleles arranged into 49 and 53 linkage groups (LGs), respectively. The recombined parental maps covered 906.3 and 970.1 cM of the genome, respectively. Finally, 264 polymorphic loci were allocated to nine LGs. The integrated map spanned 954.5 cM in length with an average genetic distance of 3.6 cM between two neighbouring loci. Quantitative trait loci (QTLs) analysis was performed using the integrated map for inflorescence diameter (ID), central disc flower diameter (CDFD), number of whorls of ray florets (NWRF), number of ray florets (NRF), number of disc florets (NDF), number of florets (NF), ray floret length (RFL), ray floret width (RFW), ray floret length/width (RFL/W), leaf length (LL), leaf width (LW) and leaf length/width (LL/W). Overall, 36 (21 major) QTLs were identified. The successful mapping of inflorescence and leaf traits QTL demonstrated the utility of the new integrated linkage map. This study is the first report of a genetic map based on EST-SSR markers in chrysanthemum. The EST-SSR markers, genetic map and QTLs reported here could be valuable resources in implementing MAS for chrysanthemums in breeding programs.

High-density genetic map construction and identification of loci controlling flower-type traits in Chrysanthemum (Chrysanthemum × morifolium Ramat.)

Flower type is an important and extremely complicated trait of chrysanthemum. The corolla tube merged degree (CTMD) and the relative number of ray florets (RNRF) are the two key factors affecting chrysanthemum flower type. However, few reports have clarified the inheritance of these two complex traits, which limits directed breeding for flower-type improvement. In this study, 305 F1 hybrids were obtained from two parents with obvious differences in CTMD and RNRF performance. Using specific-locus amplified fragment sequencing (SLAF-seq) technology, we constructed a high-density genetic linkage map with an average map distance of 0.76 cM. Three major QTLs controlling CTMD and four major QTLs underlying RNRF were repeatedly detected in the 2 years. Moreover, the synteny between the genetic map and other Compositae species was investigated, and weak collinearity was observed. In QTL regions with a high degree of genomic collinearity, eight annotated genes were probed in the Helianthus annuus L. and Lactuca sativa L. var. ramosa Hort. genomes. Furthermore, 20 and 11 unigenes were identified via BLAST searches between the SNP markers of the QTL regions and the C. vestitum and C. lavandulifolium transcriptomes, respectively. These results lay a foundation for molecular marker-assisted breeding and candidate gene exploration in chrysanthemum without a reference assembly.

Current achievements and future prospects in the genetic breeding of chrysanthemum: a review

Chrysanthemum (Chrysanthemum morifolium Ramat.) is a leading flower with applied value worldwide. Developing new chrysanthemum cultivars with novel characteristics such as new flower colors and shapes, plant architectures, flowering times, postharvest quality, and biotic and abiotic stress tolerance in a time- and cost-efficient manner is the ultimate goal for breeders. Various breeding strategies have been employed to improve the aforementioned traits, ranging from conventional techniques, including crossbreeding and mutation breeding, to a series of molecular breeding methods, including transgenic technology, genome editing, and marker-assisted selection (MAS). In addition, the recent extensive advances in high-throughput technologies, especially genomics, transcriptomics, proteomics, metabolomics, and microbiomics, which are collectively referred to as omics platforms, have led to the collection of substantial amounts of data. Integration of these omics data with phenotypic information will enable the identification of genes/pathways responsible for important traits. Several attempts have been made to use emerging molecular and omics methods with the aim of accelerating the breeding of chrysanthemum. However, applying the findings of such studies to practical chrysanthemum breeding remains a considerable challenge, primarily due to the high heterozygosity and polyploidy of the species. This review summarizes the recent achievements in conventional and modern molecular breeding methods and emerging omics technologies and discusses their future applications for improving the agronomic and horticultural characteristics of chrysanthemum.

Temporal genetic patterns of root growth in Brassica napus L. revealed by a low-cost, high-efficiency hydroponic system

Application of a low-cost and high-efficiency hydroponic system in a rapeseed population verified two types of genetic factors ("persistent" and "stage-specific") that control root development. The root system is a vital plant component for nutrient and water acquisition and is targeted to enhance plant productivity. Genetic dissection of the root system generally focuses on a single stage, but roots grow continuously during plant development. To reveal the temporal genetic patterns of root development, we measured nine root-related traits in a rapeseed recombinant inbred line population at six continuous stages during vegetative growth, using a modified hydroponic system with low-cost and high-efficiency features that could synchronize plant growth under controlled conditions. Phenotypic correlation and growth dynamic analysis suggested the existence of two types of genetic factors ("persistent" and "stage-specific") that control root development. Dynamic (unconditional and conditional) quantitative trait loci (QTL) mapping detected 28 stage-specific and 23 persistent QTLs related to root growth. Among them, 13 early stage-specific, 19 persistent and 8 later stage-specific QTLs were detected at 7 DAS (days after sowing), 16 DAS and 5 EL (expanding leaf stage), respectively, providing efficient and adaptable stages for QTL identification. The effective prediction of biomass accumulation using root morphological traits (up to 96.6% or 92.64% at a specific stage or the final stage, respectively) verified that root growth allocation with maximum root uptake area facilitated biomass accumulation. Furthermore, marker-assistant selection, which combined the "persistent" and "stage-specific" QTLs, proved their effectiveness for root improvement with an excellent uptake area. Our results highlight the potential of high-throughput and precise phenotyping to assess the dynamic genetics of root growth and provide new insights into ideotype root system-based biomass breeding.

Genome-wide association study identifies favorable SNP alleles and candidate genes for waterlogging tolerance in chrysanthemums

Chrysanthemums are sensitive to waterlogging stress, and the development of screening methods for tolerant germplasms or genes and the breeding of tolerant new varieties are of great importance in chrysanthemum breeding. To understand the genetic basis of waterlogging tolerance (WT) in chrysanthemums, we performed a genome-wide association study (GWAS) using 92,811 single nucleotide polymorphisms (SNPs) in a panel of 88 chrysanthemum accessions, including 64 spray cut and 24 disbud chrysanthemums. The results showed that the average MFVW (membership function value of waterlogging) of the disbud type (0.65) was significantly higher than that of the spray type (0.55) at P < 0.05, and the MFVW of the Asian accessions (0.65) was significantly higher than that of the European accessions (0.48) at P < 0.01. The GWAS performed using the general linear model (GLM) and mixed linear model (MLM) identified 137 and 14 SNP loci related to WT, respectively, and 11 associations were commonly predicted. By calculating the phenotypic effect values for 11 common SNP loci, six highly favorable SNP alleles that explained 12.85-21.85% of the phenotypic variations were identified. Furthermore, the dosage-pyramiding effects of the favorable alleles and the significant linear correlations between the numbers of highly favorable alleles and phenotypic values were identified (r ² = 0.45; P < 0.01). A major SNP locus (Marker6619-75) was converted into a derived cleaved amplified polymorphic sequence (dCAPS) marker that cosegregated with WT with an average efficiency of 78.9%. Finally, four putative candidate genes in the WT were identified via quantitative real-time PCR (qRT-PCR). The results presented in this study provide insights for further research on WT mechanisms and the application of molecular marker-assisted selection (MAS) in chrysanthemum WT breeding programs.

Genome-wide association study identifies favorable SNP alleles and candidate genes for waterlogging tolerance in chrysanthemums

Chrysanthemums are sensitive to waterlogging stress, and the development of screening methods for tolerant germplasms or genes and the breeding of tolerant new varieties are of great importance in chrysanthemum breeding. To understand the genetic basis of waterlogging tolerance (WT) in chrysanthemums, we performed a genome-wide association study (GWAS) using 92,811 single nucleotide polymorphisms (SNPs) in a panel of 88 chrysanthemum accessions, including 64 spray cut and 24 disbud chrysanthemums. The results showed that the average MFVW (membership function value of waterlogging) of the disbud type (0.65) was significantly higher than that of the spray type (0.55) at P < 0.05, and the MFVW of the Asian accessions (0.65) was significantly higher than that of the European accessions (0.48) at P < 0.01. The GWAS performed using the general linear model (GLM) and mixed linear model (MLM) identified 137 and 14 SNP loci related to WT, respectively, and 11 associations were commonly predicted. By calculating the phenotypic effect values for 11 common SNP loci, six highly favorable SNP alleles that explained 12.85-21.85% of the phenotypic variations were identified. Furthermore, the dosage-pyramiding effects of the favorable alleles and the significant linear correlations between the numbers of highly favorable alleles and phenotypic values were identified (r ² = 0.45; P < 0.01). A major SNP locus (Marker6619-75) was converted into a derived cleaved amplified polymorphic sequence (dCAPS) marker that cosegregated with WT with an average efficiency of 78.9%. Finally, four putative candidate genes in the WT were identified via quantitative real-time PCR (qRT-PCR). The results presented in this study provide insights for further research on WT mechanisms and the application of molecular marker-assisted selection (MAS) in chrysanthemum WT breeding programs.

Autophagy Stimulus Promotes Early HuR Protein Activation and p62/SQSTM1 Protein Synthesis in ARPE-19 Cells by Triggering Erk1/2, p38MAPK, and JNK Kinase Pathways

RNA-binding protein dysregulation and altered expression of proteins involved in the autophagy/proteasome pathway play a role in many neurodegenerative disease onset/progression, including age-related macular degeneration (AMD). HuR/ELAVL1 is a master regulator of gene expression in human physiopathology. In ARPE-19 cells exposed to the proteasomal inhibitor MG132, HuR positively affects at posttranscriptional level p62 expression, a stress response gene involved in protein aggregate clearance with a role in AMD. Here, we studied the early effects of the proautophagy AICAR + MG132 cotreatment on the HuR-p62 pathway. We treated ARPE-19 cells with Erk1/2, AMPK, p38^MAPK, PKC, and JNK kinase inhibitors in the presence of AICAR + MG132 and evaluated HuR localization/phosphorylation and p62 expression. Two-hour AICAR + MG132 induces both HuR cytoplasmic translocation and threonine phosphorylation via the Erk1/2 pathway. In these conditions, p62 mRNA is loaded on polysomes and its translation in de novo protein is favored. Additionally, for the first time, we report that JNK can phosphorylate HuR, however, without modulating its localization. Our study supports HuR's role as an upstream regulator of p62 expression in ARPE-19 cells, helps to understand better the early events in response to a proautophagy stimulus, and suggests that modulation of the autophagy-regulating kinases as potential therapeutic targets for AMD may be relevant.