Yang Cui
Baolian Fan
Xu Xu
Shasha Sheng
Yuhui Xu
Xiaoyun Wang- 1Research Center for Traditional Chinese Medicine Resources and Ethnic Minority Medicine, Jiangxi University of Chinese Medicine, Nanchang, China
- 2Adsen Biotechnology Co., Ltd., Urumchi, China
The gardenia is a traditional medicinal horticultural plant in China, but its molecular genetic research has been largely hysteretic. Here, we constructed an F1 population with 200 true hybrid individuals. Using the genotyping-by-sequencing method, a high-density sex-average genetic map was generated that contained 4,249 SNPs with a total length of 1956.28 cM and an average genetic distance of 0.46 cM. We developed 17 SNP-based Kompetitive Allele-Specific PCR markers and found that 15 SNPs were successfully genotyped, of which 13 single-nucleotide polymorphism genotypings of 96 F1 individuals showed genotypes consistent with GBS-mined genotypes. A genomic collinearity analysis between gardenia and the Rubiaceae species Coffea arabica, Coffea canephora and Ophiorrhiza pumila showed the relativity strong conservation of LG11 with NC_039,919.1, HG974438.1 and Bliw01000011.1, respectively. Lastly, a quantitative trait loci analysis at three phenotyping time points (2019, 2020, and 2021) yielded 18 QTLs for growth-related traits and 31 QTLs for leaf-related traits, of which qBSBN7-1, qCD8 and qLNP2-1 could be repeatably detected. Five QTL regions (qCD8 and qSBD8, qBSBN7 and qSI7, qCD4-1 and qLLLS4, qLNP10 and qSLWS10-2, qSBD10 and qLLLS10) with potential pleiotropic effects were also observed. This study provides novel insight into molecular genetic research and could be helpful for further gene cloning and marker-assisted selection for early growth and development traits in the gardenia.
Introduction
Gardenia (Gardenia jasminoides Ellis, 2n = 22) originated in central China, and it is a perennial shrub in the Rubiaceae family with edible flowers and medicinal fruits. Its dried ripe fruit has high quantities of crocin, geniposide, and genipin compounds (Chen Q. et al., 2020) and therefore possesses anti-inflammatory, antidepressant, anti-diabetes, antioxidant and antihypertensive activities (Qin et al., 2013; Higashino et al., 2014; Khajeh et al., 2020). The fruits are used in many traditional Chinese medicine preparations and formulas to treat different diseases (Chen L. et al., 2020). In addition to applications in traditional Chinese medicine, extracts of gardenia fruit are used as a natural colorant in the food and textile industries (Chen L. et al., 2020). Gardenia has beautiful fragrant flowers and evergreen leaves, so it is widely used for garden decoration. Fresh flowers are also used in China as edible vegetables or used to extract essential oils (Wang et al., 2017). Gardenia has a cultivation history of more than 1,000 years in China and was gradually introduced to Africa, Asia, Australia, Europe, North and South America, and the Pacific islands because of its medicinal, ornamental and industrial value (Xu et al., 2020).
Using traditional phenotypic selection-based breeding methods for genetic improvement is a labor- and time-consuming process because of the long lifecycle and highly heterozygous nature of the gardenia. By contrast, marker-assisted selection (MAS) using tightly linked or functional molecular markers with elite traits is an ideal approach to improving breeding efficiency (Mathew et al., 2014; Pootakham et al., 2015; Dong et al., 2019). However, the current molecular biology research for gardenia falls further behind model species, primarily focusing on phenotype, genetic evaluation or accession discrimination (Tsanakas et al., 2013; Hu et al., 2019; Wei et al., 2019; Li et al., 2021). Very limited studies on molecular marker identification in gardenia have been reported, such as dozens of SSR developments (Xu et al., 2014; Deng et al., 2015). Recently, a chromosomal-level genome assembly for the gardenia was released to dissect the pathway of crocin biosynthesis (Xu et al., 2020). Furthermore, helix-loop-helix (bHLH) transcription factors responsible for crocin biosynthesis were identified based on the gardenia genome (Tian et al., 2020). Genome assembly will undoubtedly accelerate functional genomics studies in gardenia. Nevertheless, the notably shortage of genome-wide molecular marker and the large gap between the phenotyping and genotyping are still bottlenecks for gardenia genetic improvement by molecular breeding and thus restrict the gardenia related industry.
Genetic maps based on the F1 segregating population are a robust tool for identifying the linkage between traits and molecular markers, which have long been applied widely in highly heterozygous species of trees, flowering plants and aquatics (Jorge et al., 2005; Wang et al., 2006; Lambert et al., 2007; Oyant et al., 2007; Sánchez-Pérez et al., 2012; Pacheco et al., 2014). In the next-generation sequencing (NGS) era, sequencing-based technologies can provide novel strategies for genome-wide SNP (single-nucleotide polymorphism) development and help to construct a high-density genetic linkage map for high-resolution QTL (quantitative trait loci) identification (Rehman et al., 2020). SNP markers can be named in many ways, including reduced-representation sequencing, resequencing and transcriptome sequencing. Reduced-representation sequencing has been differentiated into different technologies, including genotyping-by-sequencing (GBS), restriction site-associated DNA sequencing (RAD-Seq), double-digest RAD (ddRAD), specific-locus amplified fragment sequencing (SLAF-seq), ezRAD (Toonen et al., 2013) and 2b-restriction site-associated DNA sequencing (2b-RAD) (Baird et al., 2008; Elshire et al., 2011; Peterson et al., 2012; Wang et al., 2012; Sun et al., 2013; Toonen et al., 2013). Notably, GBS is a feasible SNP discovery method for highly diverse and large genome species, even without reference genome, and it has been widely adopted in genotyping for genetic map construction (İpek et al., 2017; Gabay et al., 2018; Paudel et al., 2018; Robledo et al., 2018; Lewter et al., 2019; Rubio et al., 2020). For instance, a high-density linkage map of coffee, a tree belongs to the same Rubiaceae family with gardenia, was constructed using 848 SSR and SNP markers, of which the SNP markers were developed by GBS (Moncada et al., 2016). Additional high-density genetic maps with 3,000–6,000 SNP markers have been reported in many perennial plants (Pootakham et al., 2015; Zhang et al., 2016; İpek et al., 2016; Tello et al., 2019; Zhang et al., 2019). GBS was also used for genetic diversity analysis in coffee (Anagbogu et al., 2019).
Genetic map can provide chromosome-level variation information across species. In fact, the microsynteny and macrosynteny relationship have long been verified in plants (Paterson et al., 2004; Yan et al., 2004). Comparative mapping can illustrate the co-located molecular marker distribution patterns between different genome of organisms, and further reveal structural variations and collinearity among chromosomes. Using this method, a high degree of colinearity and chromosome recombination and inversion has been found in Salicaceae species (Hanley et al., 2006; Berlin et al., 2010). Similarly, chromosomal translocations and inversions were confirmed by comparing an eggplant genetic map with the genome sequence of both tomato and pepper (Rinaldi et al., 2016). Lately, the genomic evolutionary of Coffea canephora and Ophiorrhiza pumila were investigated (Zhao et al., 2021), and some high collinearity pairs and potential karyotype rearrangement were observed, indicating their chromosomal evolution in genomic differentiation (Kai et al., 2011; Kodama et al., 2014).
QTL mapping is a traditional method to build an association bridge between genotypes and phenotypes. The tightly linked markers in QTL regions can potentially be used for MAS (Chang et al., 2018; Kim et al., 2018; Yamakawa et al., 2021). The phenotypes for typical QTL mapping always focus on specific developmental stages, and the identified QTLs represent the accumulation effect of related gene expression at the phenotyping stages. However, plant growth and development are dynamic, ever-changing processes. Dynamic QTL analysis enables QTL detection for target traits over the entire developmental process, especially for tree species, which require a relatively long time for morphogenesis. Dynamic QTL mapping studies have been published primarily for crops such as rice (Sun et al., 2015), maize (Wang et al., 2019), wheat (Mohler and Stadlmeier, 2019), cotton (Shang et al., 2015) and oilseed rape (Wang et al., 2015). In tree species, however, limited dynamic QTL maps were conducted. Desnoues et al. (2016) reported the dynamic QTL mapping of fresh weight, sugar, acid and enzyme activity at different developmental stages of peach fruit, and observed the effect of allele changes during fruit ripening. Recently, the leaf traits and plant height of Catalpa bungei at five successive time points were investigated, and a total of 33 QTLs were mapped using a high-density genetic map (Lu et al., 2019). In Populus, a total of 311 QTLs for three growth traits at 12 time points were mapped, and many QTLs specific to one time point were identified (Du et al., 2019). These results illustrated the importance of dynamic QTL mapping for the genetic dissection of developmental traits.
The genetic map of the gardenia has not been released to date. In the present study, we used a paternity test-passed F1 population of gardenia, and then employed GBS technology to construct a high-density genetic map for collinearity analysis between Rubiaceae species. Moreover, a high-resolution dynamic QTL mapping analysis was performed on growth and leaf related traits during the vegetative growth stage for three continuous years. This study was the first high-density genetic map-based QTL study in gardenia, laying a foundation for further gene cloning and MAS breeding.
Materials and Methods
Mapping Population Construction and Phenotyping
We previously screened two Gardenia jasminoides Ellis. germplasms that exhibited distinct phenotypes, namely, GD1 with high branches, large fruit, medium leaf widths and a broad crown type and AX5 with dwarf branches, small fruit, narrow leaf widths and a thin crown type. In May 2017, following emasculation at the early stage of flower development, GD1 (♀) and AX5 (♂) were crossed by artificial pollination. The dark red fruits were harvested during the first frost. The hybrid seeds were isolated and then placed on moist germination paper in Petri dishes in November 2017. At the time of radicle protrusion, the seeds were transferred into pots in the greenhouse. During the following year on March 27, the seedlings were transplanted within the Botanic Garden at Jiangxi University of Chinese Medicine (N28N°40′, E115°45′). The two parents and a total of 207 F1 individuals were randomly planted. In April 2019, young leaves from the two parents and all the F1 individuals were harvested and stored in a silica-gel drier for further DNA extraction.
Phenotyping and Data Processing
We measured 12 traits over three continuous years in October 2019, July 2020 and April 2021, and all the traits were measured three times. The detailed measurements are shown in Table 1 and Figure 1. Protractors and Vernier calipers were used to measure the stem inclinations and stem base diameters, respectively. The flexible rule was used to measure the remaining traits. SPSS V17.0 software (SPSS Inc., Chicago, IL, United States) was used for variance analysis. TBtools was used to display the variation and Pearson pairwise correlations graphically among different traits (Chen C. et al., 2020).
TABLE 1. Detailed measurement methods for the 12 agronomic traits.
FIGURE 1. Part of the trait measurement schematic diagram.
Paternity Test
To ensure an expectant hybrid seed panel, a paternity test was conducted using simple sequence repeat (SSR) markers. The total DNA from the two parents and 207 F1 individuals was isolated separately using DNA Rapid Extraction Kit DN1403 (Aidlab Biotechnologies Co., Ltd., Beijing, China). Using the total DNA from the two parents, a total of 25 SSRs from Deng et al. (2015) were used for polymorphic screening. Homozygous and polymorphic SSR markers were selected to genotype the 207 F1 individuals. For example, if the genotypes of the two parents were encoded with “aa” and “bb”, then the genotype of the true F1 offspring was “ab”. The polymerase chain reaction (PCR) system for SSR genotyping was performed in a 10.0 μl volume, with 5 μl 2×Taq MasterMix, 0.2 μl forward primer (F) and reverse primer (R), respectively 1 μl sample genomic DNA and 3.6 μl ddH2O. The system was pre-degenerated at 94°C for 3 min, and then PCR amplification began for 34 cycles of 94°C for 30 s, 55°C for 30 s, 72°C for 30 s, and a final extension at 72°C for 5 min. An 8% denaturing polyacrylamide gel was used to separate the PCR products for further silver staining.
Population GBS Sequencing and Genotyping
Similar to the paternity test, genomic DNA was isolated from the young leaves of GD1, AX5 and 200 true hybrid F1 individuals using the DNA Rapid Extraction Kit DN1403 (Aidlab Biotechnologies Co., Ltd., Beijing, China). The DNA concentration and quality were monitored using a NanoDrop spectrophotometer (ND 2000, Thermo Fisher Scientific, United States) and electrophoresis on 0.85% agarose gels, respectively. Then, GBS libraries were constructed. In brief, the genomic DNA was placed into a combination solution of RsaI and HaeIII for digestion. Products between 429 and 459 bp in length were enriched in 3% agarose gels, and end repair was performed with End Prep Enzyme Mix, followed by 3′A extension and adaptor addition. The dual index for further sample identification was introduced by PCR with eight cycles. Library quantification was performed using an Agilent 2,100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, United States), and all the libraries were mixed into one lane for paired-end sequencing (PE150) at Adsen Biotechnology Co., Ltd (Urumchi, China) using an Illumina NovaSeq 6,000 (Illumina, San Diego, CA, United States). The raw data were filtered to generate high-quality clean data according to Zhao’s criteria (Zhao et al., 2021). The genotyping was processed according to the following steps. First, a Burrows-Wheeler aligner (Li and Durbin, 2009) was used to map the clean reads to the reference genome of gardenia (Xu et al., 2020), followed by duplicate removal (Picard: http://sourceforge.net/projects/picard/). Second, SNPs were called by combining the HaplotypeCaller module of GATK (McKenna et al., 2010) and SAMtools (Li et al., 2009) to guarantee a high-quality SNP dataset. Lastly, dual-detected SNPs with sequencing depths ≥8 in the two parents, segregation distortion p > 0.01 (Chi-square) and integrity ≥60% in the offspring were maintained and encoded into eight genotyping patterns suitable for diploid species (aa × bb, ab × cd, ef ×, e.g., hk × hk, lm × ll, nn × np, ab × cc, and cc × ab). All the genotypes except aa × bb were selected and SMOOTH algorithms (van Os et al., 2005) were used to correct genotypes and imputation for further genetic map construction.
Genetic Linkage Map Construction, QTL Mapping and Gene Annotation Analyses
All the retained SNP markers were assigned into linkage groups (LGs) based on the mapping location on the reference genome of the gardenia (Xu et al., 2020). JoinMap software (V4.1) was applied for linear arrangement within LGs using the mapping function of the cross pollination (CP) model (Van Ooijen 2006). Map distances were estimated using the Kosambi mapping function (Kosambi, 1943). Genetic map visualization was performed using a CheckMatrix heat plot (http://cgpdb.ucdavis.edu/XLinkage/). The Spearman correlation coefficient between the final LGs and the reference genome was calculated and visualized using R (www.r-project.org/). MapQTL V6.0 was used for QTL analyses using the interval mapping (IM) algorithm (Van Ooijen 2009). QTLs were cut off when the LOD (logarithm of odds) values of three continuous SNPs were ≥2.5. Genes underlying stable expressed QTLs were annotated by ANNOVAR (Wang et al., 2010), and functional enrichment analyses were conducted by UniProtKB/Swiss-Prot database (Schneider et al., 2004), Pfam (Bateman et al., 2004), Gene Ontology (Ashburner et al., 2000) and KEGG (kyoto encyclopedia of genes and genomes) (Kanehisa and Goto, 2000).
Genome Synteny Analyses
To explore the evolutionary relationship between gardenia and other Rubiaceae species with chromosome-level genomes, SNP-based high-density genetic maps were aligned to the genomes of Coffea canephora (https://www.ncbi.nlm.nih.gov/genome/12248), Coffea arabica (https://www.ncbi.nlm.nih.gov/genome/?term=Coffea+arabica) and Ophiorrhiza pumila (https://www.ncbi.nlm.nih.gov/genome/97777?genome_assembly_id=1538555) using BLAST (Kent, 2002), and the physical positions of the homologous sequences were used to generate a collinearity diagram in R (www.r-project.org/).
SNP Confirmation by Kompetitive Allele-Specific PCR (KASP)
To confirm the SNPs developed by GBS, we randomly genotyped 96 F1 individuals by KASP using 17 SNPs from five randomly selected QTL regions (Supplementary Table S1). Primer 5.0 was used to design the primers, and BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi?PROGRAM=blastn&PAGE_TYPE=BlastSearch&LINK_LOC=blasthome) was used to check the primer specificity. The primer information is shown in Supplementary Table S1. The KASP genotyping processes were conducted in the GeneMatrix system (HC Scientific, Chengdu, China) according to the following three parts: Matrix Arrayer reaction plate preparation apparatus, Matrix Cycler high-throughput water bath thermal cycler, and Matrix Scanner high-speed fluorescence scanner. The PCR system contained 1 μl 2×KASP Master mix (standard ROX) (LGC Biosearch Technologies, United KIngdom), 0.028 μl KASP primer mix and 1 μl sample DNA (∼50 ng/μl). The detailed KASP thermal cycling program was 94°C for 15 min, followed by 10 cycles of 94°C for 20 s, 61–55°C for 20 s (dropping 0.6°C per cycle), 72°C for 45 s, 30 cycles of 94°C for 20 s, and 55°C for 1 min.
Results
Hybridization Test
In the present study, the parents GD1 and AX5 were used as materials, and 47 published SSRs were used for polymorphism tests. A total of 19 pairs of primers were observed to be polymorphic. We further selected markers that were homozygous and polymorphic in the parents. That is, aa×bb-type polymorphic SSRs eGJ026 and eGJ118 were used for genotyping the progeny. A total of 200 progeny out of 207 individuals in the F1 population were true hybrids with the segregation type “ab” (Supplementary Table S1–S2), indicating that the hybridization experiment was strictly controlled.
Genetic Variations in 12 Phenotypes
A set of relatively wide ranges of variations were observed in the crown diameter (CD), stem inclination (SI), plant height (PH), main stem height (MSH), leaf number per plant (LNP), the longest leaf length on stem (LLLS) and the shortest leaf length on stem (SLLS), while mild variations were present in the remaining five phenotypes (Supplementary Table S2). The coefficient of variation (CV) of the leaf-related trait LLLS and the longest leaf width on stem (LLWS) were primarily stationary at the three time points, which was similar to all six growth-related traits, suggesting minor differences among these three time points. Only the leaf-related traits LNS, LNP, SLLS and SLWS demonstrated acute CV fluctuation (Supplementary Table S2). Among every year for the phenotypes at the three time point, CD, PH and LNP exhibited strong correlations with other phenotypes. However, SI had no highly significant correlations with the other ten traits (except SBD). PH had the strongest positive correlations with MSH, with correlation coefficients of 0.95, 0.93 and 0.88 in 2019, 2020 and 2021, respectively. SLLS and SLWS also exhibited a correlation greater than 0.80 in the 3 years (Figure 2). The correlation analysis implied that there was an independent and interdependent relationship between growth-related traits and leaf-related traits.
FIGURE 2. Variation and Pearson pairwise correlation analyses of growth-related and leaf-related traits of the F1 population. (A), (B) and (C) represent the variation and Pearson pairwise correlations in 2019, 2020 and 2021, respectively. The correlations were calculated with Spearman correlation coefficients, and the p values are indicated as follows: *, p < 0.05; **, p < 0.01; and ***, p < 0.001. The abbreviations given in the histograms are as follows: CD: crown diameter; BSBN: basal stem branch number; SI: stem inclination; PH: plant height; MSH: main stem height; SBD: stem base diameter; LNS: leaf number on stem; LNP: leaf number per plant; LLLS: longest leaf length on stem; LLWS: longest leaf width on stem; SLLS: shortest leaf length on stem; and SLWS: shortest leaf width on stem.
Variation Calling and Genotyping
In total, GBS sequencing generated 29, 630, 679 clean reads after quality control, with 1.77 and 0.83 Gb for the parent AX5 and GD1, respectively. For the offspring, 12, 186, 237 reads (1.81 Gb) were obtained per individual. The statistics showed that the average Q30 was higher than 85%, and the GC content (%) was distributed between 40.29 and 48.18 (Supplementary Table S3). Upon using BWA software to align the sequencing data to the reference genome of gardenia, the mapping rates were 94.38, 93.88 and 96.66% for AX5, GD1 and all the progeny, respectively. These pre-processing procedures indicated a high quality of sequencing data for further analysis.
A total of 154,909 SNPs were detected by combining the GATK and SAMTools, and the genotypes of these SNPs were encoded into eight segregation patterns. Among them, lm×ll, np×nn, hk×hk and aa×bb occupied 47,244, 63,583, 21,048 and 21,834 SNPs, respectively, accounting for 99.23% of the total SNPs (Figure 3). After depth and integrity filtering, the remaining SNPs were used for genetic map construction.
FIGURE 3. The distributions of SNP marker segregation patterns.
High Density Genetic Map
JoinMap was used to construct a female genetic map containing 2,585 markers spanning 2,348.48 cM and a male genetic map containing 1,963 markers spanning 1,348.38 cM, and both consisted of 11 LGs (Table 2). Integrating the female and male maps formed a sex-average genetic map, which included 4,249 SNPs with a total length of 1956.28 cM and an average genetic distance of 0.46 cM (Table 2; Figure 4). Among them, LG2 was the longest (228.19 cM), including 648 SNPs, and the average genetic distance was 0.35 cM. Conversely, LG3 was the shortest group (110.56 cM) with 468 SNP tags but a higher resolution of 0.24 cM between adjacent markers on average. There were 120 SNPs in LG9, which was the least in all 11 LGs, with 1.19 cM, the largest average distance between adjacent markers. Of the 11 LGs, the proportions of genetic gaps (≤ 5 cM) ranged from 92.44 to 99.79%, with 96.86% on average. No genetic gaps larger than 10 cM were observed in LG2, LG4 or LG8. The largest gap presented on LG7, which was 23.86 cM (Table 2). The detailed information of all the SNPs and the corresponding genetic and physical positions were displayed in Supplementary Table S4.
TABLE 2. The basic characteristics of the female genetic map, male genetic map and sex-average genetic map.
FIGURE 4. High-density sex-average genetic map of gardenia.
The Spearman correlation coefficient between the genetic map and the reference genome was approximately 0.722 and up to 0.901, indicating that the marker ordering of the genetic map was basically accurate (Supplementary Table S5).
QTLs for Growth and Leaf-Related Traits
We divided the 12 traits into two categories, namely, phenotypes associated with gardenia growth (CD, BSBN, SI, PH, MSH, and SBD) and leaves (LNS, LNP, LLLS, LLWS, SLLS, and SLWS). Using the high-density genetic map and continuous phenotypic data, 23, 22, and 9 QTLs were mapped in 2019, 2020 and 2021, respectively, of which 18 QTLs were associated with the growth traits, while 31 QTLs were related to the leaf-related traits (Table 3). These QTLs were distributed in all the LGs of gardenia except LG6, with phenotypic variance explained (PVE) values ranging from 5.6 to 11.7%.
TABLE 3. QTL mapping results.
Eighteen QTLs were detected for the six growth traits except PH, including nine for CD, two for BSBN, one for MSH, three for SI, and three for SBD (Table 3). Two major-effect QTLs, qBSBN7-1 and qBSBN7-2 (PVE = 10–10.1%), were identified with LOD>4, and qBSBN7-1 could be detected in both 2020 and 2021 (Table 3). The LOD of qCD4-1 in 2019 was slightly over 2.5 but was supported by 31 SNPs. Notably, qCD8, which was located in the 182.843–184.973 cM interval, could be repeatably mapped over 3 years, indicating that this QTL could be expressed continuously, contributing to the establishment of the crown (Figure 5). Especially in 2020, the LOD value of qCD8 exceeded four, and PVE exceeded 9%. For SI, there were two QTLs (qSI4-1 and qSI4-2) with a spacing of only 1.322 cM, which were supported by 16 and 12 SNPs, respectively.
FIGURE 5. Repeatable QTLs for three dynamic phenotyping time points.
In the QTLs of the leaf-related phenotype, four QTLs for LLLS were mapped to LG4, LG9, LG10 and LG11, of which the highest LOD value of qLLLS9 in 2020 was up to 5.3, and the corresponding PVE was equal to 11.7%. The largest number of identified QTLs belonged to trait LNP, up to 12, and was primarily distributed on LG10, of which qLNP2-1, located at 19.951–28.282 cM, was detected for three consecutive years (Table 3; Figure 5). qLNP2-2 and qLNP2-10 were 0.502 cM apart, implying that they might be the same QTL. A total of five QTLs were responsible for SLLS, of which qSLLS11-1, qSLLS11-2 and qSLLS11-3 gathered between 211.278 and 223.331 cM, with 41 supported SNPs. There were three QTLs for LNS on LG7, LG8 and LG9. With respect to SLWS and LLWS, three and two QTLs were mapped, respectively.
We continued to explore the QTLs among different phenotypes and found that there were five pairs of QTLs with shared regions, including qCD8 and qSBD8, qBSBN7 and qSI7, qCD4-1 and qLLLS4, qLNP10 and qSLWS10-2, qSBD10 and qLLLS10 (Table 3; Figure 6), suggesting that each pair underlying a single QTL and pleiotropism might play a significant role in gardenia morphogenesis and vegetative development. The structural and functional gene annotations of the above stable and potential pleiotropism QTLs were isolated, resulting in 2,514 nonredundant genes and the corresponding annotation information (Supplementary Table S6–S7).
FIGURE 6. Potential pleiotropism QTLs.
KASP-Based SNP Confirmation
We selected 17 SNP-based KASP markers on chromosomes 2, 7, 9, and 11 and a total of 96 samples for SNP accuracy verification. There were 15 out of 17 markers with successful fluorescence signals in the HC KASP platform, accounting for 88.24%. Among these 15 markers, 13 SNP markers showed genotypes consistent with the GBS results of each individual (Supplementary Table S8), indicating the accuracy of the sequencing analysis.
Synteny Analyses
We used this high-density genetic map to investigate the evolutionary relationship of Rubiaceae species, as shown in Figure 7. Different levels of synteny were observed between the LGs of gardenia between C. arabica (A), C. canephora (B) and O. pumila (C). Specifically, relatively strong synteny was consistently noted between LG11 of gardenia and NC_039,919.1 of C. arabica, HG974438.1 of C. canephora, and Bliw01000011.1 of O. pumila, indicating that LG11 was more conserved than other LGs. In addition, this type of stronger collinearity was also found between these pairs: the pair LG10 and C. arabica’s chromosome NC_039,919.1, the pair LG4 and C. canephora’s chromosome HG974438.1, and the pair LG10 and O. pumila’s chromosome BLIW01000011.1. Moreover, chromosomes NC_039,917.1, NC_039,918.1 and NC_039,919.1 demonstrated higher collinearity than other chromosomes in C. arabica. Similarly, chromosomes HG974438.1, HG974437.1 and HG974436.1 in C. canephora and chromosomes BLIW01000010.1 and BLIW01000011.1 in O. pumila displayed stronger synteny than other corresponding chromosomes.
FIGURE 7. Synteny analyses between the genetic map of gardenia and the genomes of Coffea arabica (A), Coffea canephora (B), and Ophiorrhiza pumila (C).
Discussion
Gardenia is a type of gardening species that has medicinal and industrial value. At present, there is less genetic research on this species. In this study, after constructing an F1 segregating population, the first high-density genetic map of gardenia was accomplished using a high-throughput sequencing method.
Three-year dynamic QTL positioning identified a panel of vegetative growth-related QTLs. We believe that this research will open a new avenue for gardenia molecular genetic research.
NGS has greatly accelerated the process of QTL mapping according to forward genetics, as primarily performed by bulk segregation analysis (BSA), high-density genetic map-based QTL mapping and genome-wide association studies (GWASs) in horticultural plants (Ban and Xu, 2020; Ferrão et al., 2020; Song et al., 2020). GBS, which is rooted in NGS, has opened up new possibilities for genome-wide SNP mining without high investments (Chung et al., 2017). At the cost of approximately 1.8 Gb sequencing data per sample, 154,909 dual-calling SNPs were de novo developed in the present study. Hereafter, the sequencing depth, segregation distortion and integrity control were processed to guarantee a high-quality panel of SNPs for high-density genetic map construction. Using an F1 population of 200 plants of gardenia and GBS-based genotyping, a high-density genetic map harboring 4,249 SNPs was constructed, which showed high resolution (0.46 cM per adjacent SNPs) and satisfied marker orders with an approximately 0.8 collinearity compared with the reference genome, a similar standard as in other species (Ji et al., 2018; Lu et al., 2019; Wang et al., 2020). To evaluate the SNP accuracy, we selected 17 SNPs and genotyped 96 F1 individuals, and high consistency between the KASP and GBS genotypes was observed according to Song’s report (Song et al., 2020). These results indicated that a high-quality and high-density genetic map was generated. Furthermore, the genetic map was used for the QTL mapping responsible for early growth and development traits for three dynamic years, and a total of 49 QTLs for 12 traits (CD, BSBN, SI, PH, MSH, SBD, LNS, LNP, LLLS, LLWS, SLLS, and SLWS) were identified. Because the most useful organ was the gardenia fruit, further QTLs associated with fruit-related traits, such as the weight, shape, size or functional substance content, could be expected in 2022 or later, when all F1 gardenia individuals will transfer to reproductive growth. This high-density genetic map might provide a new lesson for molecular genetic research in gardenia.
As mentioned above, 49 QTLs for 12 traits were mapped, including three stably expressed QTLs, qBSBN7-1, qCD8 and qLNP2-1 (Table 3). These QTLs played persistent roles in the corresponding morphogenesis during the vegetative period and were valuable and useful for MAS-based breeding programs (Jamshed et al., 2016; Galiano-Carneiro et al., 2021). In addition, we also found five regions within two QTLs (qCD8 and qSBD8, qBSBN7 and qSI7, qCD4-1 and qLLLS4, qLNP10 and qSLWS10-2, qSBD10 and qLLLS10), suggesting that each region underlies a single QTL with pleiotropism (Sun et al., 2017; Wang et al., 2018). Notably, there was a genetic basis for the phenotypic correlation between CD and SBD, and CD and LLLS, which was consistent with the strong correlations between CD and SBD (0.40***), and CD and LLLS (0.53***) (Figure 2). The SNPs underlying these stable and pleiotropic QTLs could be further converted into KASP markers and potentially used as MAS markers. Further gene cloning may also benefit from the gene structural and functional annotations underlying stable and potential pleiotropism QTLs (Supplementary Table S6, S7).
Generally, gardenia plants are focused on vegetative development for the first 3 years after seed germination, which is called the juvenile period. A gradual declining trend in QTL numbers was observed from 2019, 2020 to 2021, which might be associated with the fact that the channels of vegetative development gradually slowed, while reproductive growth gradually opened. Furthermore, some year-specific QTLs were found except stable and pleiotropism potential QTLs, which could be explained by specific functions in plant growth phases. This phenotyping-derived dynamic QTL mapping based on continuous development time points has recently been performed in peaches (Desnoues et al., 2016), Populus (Du et al., 2019), Catalpa bungei (Lu et al., 2019), and chrysanthemum (Ao et al., 2019). Compared with these studies, one shortcoming of this study is that phenotyping was conducted at only three time points, which might limit the dissection resolution of the developmental trait inheritance. Currently, high-throughput phenomics combining spectral imaging and machine learning methods provides particular insight into the deciphering of dynamic phenotypes in a way that is plant damage-free (Li and Sillanpää 2015; Adams et al., 2020; Streich et al., 2020; Zhu et al., 2021). Novel dynamic QTL perspectives might be enabled by employing phenomic methods in the near future.
The complex synteny of gardenia with C. arabica, C. canephora and O. pumila (Figure 7) implied that widely chromosomal fission and fusion have happened after their divergence from the common ancestor, similar to other species (Luo et al., 2020; Yang et al., 2020). The synteny levels of all LGs in gardenia indicated that different chromosomes underwent different evolutionary process. Strong synteny in LG10, LG11 and LG4 of gardenia demonstrated that these LGs were more conserved than other LGs. The sequences underlying these conserved regions could be potentially used to speculate the corresponding genetic information of species in Salicaceae, and the QTLs on LG4, LG10 and LG11 (qCD4-1, qSBD10, qSLLS11-3, etc.), potentially orthologous QTLs (Rinaldi et al., 2016; Webb et al., 2016), could be applied for comparative mapping in other species of Rubiaceae. Ta°Cken together, the synteny analyses of this paper may lay a foundation for subsequent comparative genomic research.
One ultimate goal of QTL mapping is to perform QTL fine mapping, screening and gene cloning of candidate genes. This goal was commonly achieved in therophyte plants, such as through map-based cloning (Jia et al., 2020; Yu et al., 2020; Sierra-Orozco et al., 2021). However, for many perennial species, it takes more than 6 years to perform hybridization and further backcrossing, and the population size is always restricted to several hundred, which leads to limited recombination. This limitation makes the process of QTL cloning in perennial species slow. To expedite this process, transcriptomics analysis can be used to call RNA variants and differentially expressed genes (DEGs) within QTL regions (Park et al., 2019; Wen et al., 2019). Other omics analyses, such as metabolomics and proteomics, can also provide useful information on the metabolic chemicals or proteins related to phenotypic variations (Szymański et al., 2020; Mou et al., 2021). In future experimental designs, these analyses will be considered to understand the basis of phenotypic variation comprehensively.
5 Conclusion
In this study, we developed a panel of genome-wide high-quality SNPs using the GBS method and providing the first high-density genetic map in the gardenia. SNPs and genetic maps could be useful for further genetic study and evolutionary genomics. Based on this high-density genetic map, 18 and 31 QTLs were identified for growth traits and leaf-related traits at three dynamic phenotyping time points, respectively. Stably expressed QTLs and potential pleiotropism QTLs could be targets for MAS breeding and for further gene cloning.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
XW conceived and designed the experiments; YC performed paternity test, phenotying and first draft writing; BF and SS performed phenotying; YX performed the data analysis; YC wrote the manuscript and XW, YX revised the manuscript. All authors have read and approved the final manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (grant number 31660218) and the Scientific Foundation of Double First-class Discipline Development of TCM (grant number JXSYLXK-ZHYA0028). National Natural Science Foundation of China is a national public funder. Scientific Foundation of Double First-class Discipline Development of TCM is a funder of Jiangxi University of Chinese Medicine.
Conflict of Interest
YX was employed by Adsen Biotechnology Co., Ltd.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.802738/full#supplementary-material
Abbreviations
SNPs, single-nucleotide polymorphisms; QTLs, quantitative trait loci; LOD, logarithm of odds; MAS, marker-assisted selection; NGS, next-generation sequencing; GBS, genotyping-by-sequencing; CD, crown diameter; BSBN, basal stem branch number; SI, stem inclination; PH, plant height; MSH, main stem height; SBD, stem base diameter; LNS, leaf number on stem; LNP, leaf number per plant; LLLS, longest leaf length on stem; LLWS, longest leaf width on stem; SLLS, shortest leaf length on stem; SLWS, shortest leaf width on stem; CV, coefficient of variation; LGs, linkage groups; PVE, phenotypic variance explained; KASP, kompetitive allele-specific PCR; CP, cross pollination; DEGs, differentially expressed genes.
References
Adams, J., Qiu, Y., Xu, Y., and Schnable, J. C. (2020). Plant Segmentation by Supervised Machine Learning Methods. Plant phenome j. 3, e20001. doi:10.1002/ppj2.20001
Anagbogu, C. F., Bhattacharjee, R., Ilori, C., Tongyoo, P., Dada, K. E., Muyiwa, A. A., et al. (2019). Genetic Diversity and Re-classification of Coffee (Coffea Canephora Pierre Ex A. Froehner) from South Western Nigeria through Genotyping-By-Sequencing-Single Nucleotide Polymorphism Analysis. Genet. Resour. Crop Evol. 66, 685–696. doi:10.1007/s10722-019-00744-2
Ao, N., Ma, J., Xu, T., Su, J., Yang, X., Guan, Z., et al. (2019). Genetic Variation and QTL Mapping for Cold Tolerance in a chrysanthemum F1 Population at Different Growth Stages. Euphytica 215, 88. doi:10.1007/s10681-019-2412-7
Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene Ontology: Tool for the Unification of Biology. Nat. Genet. 25, 25–29. doi:10.1038/75556
Baird, N. A., Etter, P. D., Atwood, T. S., Currey, M. C., Shiver, A. L., Lewis, Z. A., et al. (2008). Rapid SNP Discovery and Genetic Mapping Using Sequenced RAD Markers. PLoS One 3, e3376. doi:10.1371/journal.pone.0003376
Ban, S., and Xu, K. (2020). Identification of Two QTLs Associated with High Fruit Acidity in Apple Using Pooled Genome Sequencing Analysis. Hortic. Res. 7, 171. doi:10.1038/s41438-020-00393-y
Bateman, A., Coin, L., Durbin, R., Finn, R. D., Hollich, V., Griffiths-Jones, S., et al. (2004). The Pfam Protein Families Database. Nucleic Acids Res. 32, 138D–141D. doi:10.1093/nar/gkh121
Berlin, S., Lagercrantz, U., von Arnold, S., Öst, T., and Rönnberg-Wästljung, A. (2010). High-density Linkage Mapping and Evolution of Paralogs and Orthologs in Salix and Populus. BMC Genomics 11, 129. doi:10.1186/1471-2164-11-129
Chang, Y., Ding, J., Xu, Y., Li, D., Zhang, W., Li, L., et al. (2018). SLAF-based High-Density Genetic Map Construction and QTL Mapping for Major Economic Traits in Sea Urchin Strongylocentrotus Intermedius. Sci. Rep. 8, 820. doi:10.1038/s41598-017-18768-y
Chen, C., Chen, H., Zhang, Y., Thomas, H. R., Frank, M. H., He, Y., et al. (2020c). TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol. Plant 13, 1194–1202. doi:10.1016/j.molp.2020.06.009
Chen, L., Li, M., Yang, Z., Tao, W., Wang, P., Tian, X., et al. (2020b). Gardenia Jasminoides Ellis: Ethnopharmacology, Phytochemistry, and Pharmacological and Industrial Applications of an Important Traditional Chinese Medicine. J. Ethnopharmacology 257, 112829. doi:10.1016/j.jep.2020.112829
Chen, Q., Xue, G., Ni, Q., Wang, Y., Gao, Q., Zhang, Y., et al. (2020a). Physicochemical and Rheological Characterization of Pectin‐rich Polysaccharides from Gardenia Jasminoides J. Ellis Flower. Food Sci. Nutr. 8, 3335–3345. doi:10.1002/fsn3.1612
Chung, Y. S., Choi, S. C., Jun, T.-H., and Kim, C. (2017). Genotyping-by-Sequencing: a Promising Tool for Plant Genetics Research and Breeding. Hortic. Environ. Biotechnol. 58, 425–431. doi:10.1007/s13580-017-0297-8
Deng, S.-Y., Wang, X.-r., Zhu, P.-l., Wen, Q., and Yang, C.-x. (2015). Development of Polymorphic Microsatellite Markers in the Medicinal Plant Gardenia Jasminoides (Rubiaceae). Biochem. Syst. Ecol. 58, 149–155. doi:10.1016/j.bse.2014.11.009
Desnoues, E., Baldazzi, V., Génard, M., Mauroux, J.-B., Lambert, P., Confolent, C., et al. (2016). Dynamic QTLs for Sugars and Enzyme Activities Provide an Overview of Genetic Control of Sugar Metabolism during Peach Fruit Development. Exbotj 67, 3419–3431. doi:10.1093/jxb/erw169
Dong, M., He, Q., Zhao, J., Zhang, Y., Yuan, D., and Zhang, A. J. (2019). Genetic Mapping of Prince Rupprecht's Larch (Larix Principis-rupprechtii Mayr) by Specific-Locus Amplified Fragment Sequencing. Genes 10, 583. doi:10.3390/genes10080583
Du, Q., Yang, X., Xie, J., Quan, M., Xiao, L., Lu, W., et al. (2019). Time-specific and Pleiotropic Quantitative Trait Loci Coordinately Modulate Stem Growth inPopulus. Plant Biotechnol. J. 17, 608–624. doi:10.1111/pbi.13002
Elshire, R. J., Glaubitz, J. C., Sun, Q., Poland, J. A., Kawamoto, K., Buckler, E. S., et al. (2011). A Robust, Simple Genotyping-By-Sequencing (GBS) Approach for High Diversity Species. PLoS One 6, e19379. doi:10.1371/journal.pone.0019379
Ferrão, L. F. V., Johnson, T. S., Benevenuto, J., Edger, P. P., Colquhoun, T. A., and Munoz, P. R. (2020). Genome‐wide Association of Volatiles Reveals Candidate Loci for Blueberry Flavor. New Phytol. 226, 1725–1737. doi:10.1111/nph.16459
Gabay, G., Dahan, Y., Izhaki, Y., Faigenboim, A., Ben-Ari, G., Elkind, Y., et al. (2018). High-resolution Genetic Linkage Map of European Pear (Pyrus Communis) and QTL fine-mapping of Vegetative Budbreak Time. BMC Plant Biol. 18, 175. doi:10.1186/s12870-018-1386-2
Galiano-Carneiro, A. L., Kessel, B., Presterl, T., and Miedaner, T. (2021). Intercontinental Trials Reveal Stable QTL for Northern Corn Leaf Blight Resistance in Europe and in Brazil. Theor. Appl. Genet. 134, 63–79. doi:10.1007/s00122-020-03682-1
Hanley, S. J., Mallott, M. D., and Karp, A. (2006). Alignment of a Salix Linkage Map to the Populus Genomic Sequence Reveals Macrosynteny between Willow and poplar Genomes. Tree Genet. Genomes 3, 35–48. doi:10.1007/s11295-006-0049-x
Higashino, S., Sasaki, Y., Giddings, J. C., Hyodo, K., Fujimoto Sakata, S., Matsuda, K., et al. (2014). Crocetin, a Carotenoid fromGardenia jasminoidesEllis, Protects against Hypertension and Cerebral Thrombogenesis in Stroke-Prone Spontaneously Hypertensive Rats. Phytother. Res. 28, 1315–1319. doi:10.1002/ptr.5130
Hu, Y., Liu, X., Xia, Q., Yin, T., Bai, C., Wang, Z., et al. (2019). Comparative Anti-arthritic Investigation of Iridoid Glycosides and Crocetin Derivatives from Gardenia Jasminoides Ellis in Freund's Complete Adjuvant-Induced Arthritis in Rats. Phytomedicine 53, 223–233. doi:10.1016/j.phymed.2018.07.005
İpek, A., İpek, M., Ercişli, S., and Tangu, N. A. (2017). Transcriptome-based SNP Discovery by GBS and the Construction of a Genetic Map for Olive. Funct. Integr. Genomics 17, 493–501. doi:10.1007/s10142-017-0552-1
İpek, A., Yılmaz, K., Sıkıcı, P., Tangu, N. A., Öz, A. T., Bayraktar, M., et al. (2016). SNP Discovery by GBS in Olive and the Construction of a High-Density Genetic Linkage Map. Biochem. Genet. 54, 313–325. doi:10.1007/s10528-016-9721-5
Jamshed, M., Jia, F., Gong, J., Palanga, K. K., Shi, Y., Li, J., et al. (2016). Identification of Stable Quantitative Trait Loci (QTLs) for Fiber Quality Traits across Multiple Environments in Gossypium Hirsutum Recombinant Inbred Line Population. BMC Genomics 17, 197. doi:10.1186/s12864-016-2560-2
Ji, F., Wei, W., Liu, Y., Wang, G., Zhang, Q., Xing, Y., et al. (2018). Construction of a SNP-Based High-Density Genetic Map Using Genotyping by Sequencing (GBS) and QTL Analysis of Nut Traits in Chinese Chestnut (Castanea Mollissima Blume). Front. Plant Sci. 9, 816. doi:10.3389/fpls.2018.00816
Jia, H., Li, M., Li, W., Liu, L., Jian, Y., Yang, Z., et al. (2020). A Serine/threonine Protein Kinase Encoding Gene KERNEL NUMBER PER ROW6 Regulates maize Grain Yield. Nat. Commun. 11, 988. doi:10.1038/s41467-020-14746-7
Jorge, V., Dowkiw, A., Faivre-Rampant, P., and Bastien, C. (2005). Genetic Architecture of Qualitative and Quantitative Melampsora Larici-populina Leaf Rust Resistance in Hybrid poplar: Genetic Mapping and QTL Detection. New Phytol. 167, 113–127. doi:10.1111/j.1469-8137.2005.01424.x
Kai, W., Kikuchi, K., Tohari, S., Chew, A. K., Tay, A., Fujiwara, A., et al. (2011). Integration of the Genetic Map and Genome Assembly of Fugu Facilitates Insights into Distinct Features of Genome Evolution in Teleosts and Mammals. Genome Biol. Evol. 3, 424–442. doi:10.1093/gbe/evr041
Kanehisa, M., and Goto, S. (2000). KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 28, 27–30. doi:10.1093/nar/28.1.27
Kent, W. J. (2002). BLAT-the BLAST-like Alignment Tool. Genome Res. 12, 656–664. doi:10.1101/gr.229202
Khajeh, E., Rasmi, Y., Kheradmand, F., Malekinejad, H., Aramwit, P., Saboory, E., et al. (2020). Crocetin Suppresses the Growth and Migration in HCT-116 Human Colorectal Cancer Cells by Activating the P-38 MAPK Signaling Pathway. Res. Pharma Sci. 15, 592–601. doi:10.4103/1735-5362.301344
Kim, B., Hwang, I. S., Lee, H. J., Lee, J. M., Seo, E., Choi, D., et al. (2018). Identification of a Molecular Marker Tightly Linked to Bacterial Wilt Resistance in Tomato by Genome-wide SNP Analysis. Theor. Appl. Genet. 131, 1017–1030. doi:10.1007/s00122-018-3054-1
Kodama, M., Brieuc, M. S. O., Devlin, R. H., Hard, J. J., and Naish, K. A. (2014). Comparative Mapping between Coho salmon (Oncorhynchus kisutch) and Three Other Salmonids Suggests a Role for Chromosomal Rearrangements in the Retention of Duplicated Regions Following a Whole Genome Duplication Event. G3-genes Genom. Genet. 4, 1717–1730. doi:10.1534/g3.114.012294
Kosambi, D. D. (1943). The Estimation of Map Distances from Recombination Values. Ann. Eugen. 12, 172–175. doi:10.1111/j.1469-1809.1943.tb02321.x
Lambert, P., Dicenta, F., Rubio, M., and Audergon, J. M. (2007). QTL Analysis of Resistance to Sharka Disease in the Apricot (Prunus Armeniaca L.) 'Polonais' × 'Stark Early Orange' F1 Progeny. Tree Genet. Genomes 3, 299–309. doi:10.1007/s11295-006-0069-6
Lewter, J., Worthington, M. L., Clark, J. R., Varanasi, A. V., Nelson, L., Owens, C. L., et al. (2019). High-density Linkage Maps and Loci for berry Color and Flower Sex in Muscadine Grape (Vitis Rotundifolia). Theor. Appl. Genet. 132, 1571–1585. doi:10.1007/s00122-019-03302-7
Li, H., and Durbin, R. (2009). Fast and Accurate Short Read Alignment with Burrows-Wheeler Transform. Bioinformatics 25, 1754–1760. doi:10.1093/bioinformatics/btp324
Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., et al. (2009). The Sequence Alignment/Map Format and SAMtools. Bioinformatics 25, 2078–2079. doi:10.1093/bioinformatics/btp352
Li, W., Ren, C., Fei, C., Wang, Y., Xue, Q., Li, L., et al. (2021). Analysis of the Chemical Composition Changes of Gardeniae Fructus before and after Processing Based on Ultra‐high‐performance Liquid Chromatography Quadrupole Time‐of‐flight Mass Spectrometry. J. Sep. Sci. 44, 981–991. doi:10.1002/jssc.202000957
Li, Z., and Sillanpää, M. J. (2015). Dynamic Quantitative Trait Locus Analysis of Plant Phenomic Data. Trends Plant Sci. 20, 822–833. doi:10.1016/j.tplants.2015.08.012
Lu, N., Zhang, M., Xiao, Y., Han, D., Liu, Y., Zhang, Y., et al. (2019). Construction of a High-Density Genetic Map and QTL Mapping of Leaf Traits and Plant Growth in an Interspecific F1 Population of Catalpa Bungei × Catalpa Duclouxii Dode. BMC Plant Biol. 19, 596. doi:10.1186/s12870-019-2207-y
Luo, X., Xu, L., Wang, Y., Dong, J., Chen, Y., Tang, M., et al. (2020). An Ultra‐high‐density Genetic Map Provides Insights into Genome Synteny, Recombination Landscape and Taproot Skin Colour in Radish ( Raphanus Sativus L). Plant Biotechnol. J. 18, 274–286. doi:10.1111/pbi.13195
Mathew, L. S., Spannagl, M., Al-Malki, A., George, B., Torres, M. F., Al-Dous, E. K., et al. (2014). A First Genetic Map of Date palm (Phoenix Dactylifera) Reveals Long-Range Genome Structure Conservation in the Palms. BMC Genomics 15, 285. doi:10.1186/1471-2164-15-285
McKenna, A., Hanna, M., Banks, E., Sivachenko, A., Cibulskis, K., Kernytsky, A., et al. (2010). The Genome Analysis Toolkit: A MapReduce Framework for Analyzing Next-Generation DNA Sequencing Data. Genome Res. 20, 1297–1303. doi:10.1101/gr.107524.110
Mohler, V., and Stadlmeier, M. (2019). Dynamic QTL for Adult Plant Resistance to Powdery Mildew in Common Wheat (Triticum aestivum L). J. Appl. Genet. 60, 291–300. doi:10.1007/s13353-019-00518-7
Moncada, M. D. P., Tovar, E., Montoya, J. C., González, A., Spindel, J., and McCouch, S. (2016). A Genetic Linkage Map of Coffee (Coffea Arabica L.) and QTL for Yield, Plant Height, and Bean Size. Tree Genet. Genomes 12, 5. doi:10.1007/s11295-015-0927-1
Mou, J., Zhang, Z., Qiu, H., Lu, Y., Zhu, X., Fan, Z., et al. (2021). Multiomics-based Dissection of Citrus Flavonoid Metabolism Using a Citrus Reticulata × Poncirus Trifoliata Population. Hortic. Res. 8, 56. doi:10.1038/s41438-021-00472-8
Oyant, L. H.-S., Crespel, L., Rajapakse, S., Zhang, L., and Foucher, F. (2007). Genetic Linkage Maps of Rose Constructed with New Microsatellite Markers and Locating QTL Controlling Flowering Traits. Tree Genet. Genomes 4, 11–23. doi:10.1007/s11295-007-0084-2
Pacheco, I., Bassi, D., Eduardo, I., Ciacciulli, A., Pirona, R., Rossini, L., et al. (2014). QTL Mapping for Brown Rot (Monilinia Fructigena) Resistance in an Intraspecific Peach (Prunus Persica L. Batsch) F1 Progeny. Tree Genet. Genomes 10, 1223–1242. doi:10.1007/s11295-014-0756-7
Park, M., Lee, J.-H., Han, K., Jang, S., Han, J., Lim, J.-H., et al. (2019). A Major QTL and Candidate Genes for Capsaicinoid Biosynthesis in the Pericarp of Capsicum Chinense Revealed Using QTL-Seq and RNA-Seq. Theor. Appl. Genet. 132, 515–529. doi:10.1007/s00122-018-3238-8
Paterson, A. H., Bowers, J. E., and Chapman, B. A. (2004). Ancient Polyploidization Predating Divergence of the Cereals, and its Consequences for Comparative Genomics. Proc. Natl. Acad. Sci. 101, 9903–9908. doi:10.1073/pnas.0307901101
Paudel, D., Kannan, B., Yang, X., Harris-Shultz, K., Thudi, M., Varshney, R. K., et al. (2018). Surveying the Genome and Constructing a High-Density Genetic Map of Napiergrass (Cenchrus Purpureus Schumach). Sci. Rep. 8, 14419. doi:10.1038/s41598-018-32674-x
Peterson, B. K., Weber, J. N., Kay, E. H., Fisher, H. S., and Hoekstra, H. E. (2012). Double Digest RADseq: An Inexpensive Method for De Novo SNP Discovery and Genotyping in Model and Non-model Species. PLoS One 7, e37135. doi:10.1371/journal.pone.0037135
Pootakham, W., Ruang-Areerate, P., Jomchai, N., Sonthirod, C., Sangsrakru, D., Yoocha, T., et al. (2015). Construction of a High-Density Integrated Genetic Linkage Map of Rubber Tree (Hevea Brasiliensis) Using Genotyping-By-Sequencing (GBS). Front. Plant Sci. 6, 367. doi:10.3389/fpls.2015.00367
Qin, F.-m., Meng, L.-j., Zou, H.-l., and Zhou, G.-x. (2013). Three New Iridoid Glycosides from the Fruit of Gardenia Jasminoides Var. Radicans. Chem. Pharm. Bull. 61, 1071–1074. doi:10.1248/cpb.c13-00262
Rehman, F., Gong, H., Li, Z., Zeng, S., Yang, T., Ai, P., et al. (2020). Identification of Fruit Size Associated Quantitative Trait Loci Featuring SLAF Based High-Density Linkage Map of Goji berry (Lycium spp). BMC Plant Biol. 20, 474. doi:10.1186/s12870-020-02567-1
Rinaldi, R., Van Deynze, A., Portis, E., Rotino, G. L., Toppino, L., Hill, T., et al. (2016). New Insights on Eggplant/tomato/pepper Synteny and Identification of Eggplant and Pepper Orthologous QTL. Front. Plant Sci. 7, 1031. doi:10.3389/fpls.2016.01031
Robledo, D., Palaiokostas, C., Bargelloni, L., Martínez, P., and Houston, R. (2018). Applications of Genotyping by Sequencing in Aquaculture Breeding and Genetics. Rev. Aquacult. 10, 670–682. doi:10.1111/raq.12193
Rubio, B., Lalanne-Tisné, G., Voisin, R., Tandonnet, J.-P., Portier, U., Van Ghelder, C., et al. (2020). Characterization of Genetic Determinants of the Resistance to Phylloxera, Daktulosphaira Vitifoliae, and the Dagger Nematode Xiphinema index from Muscadine Background. BMC Plant Biol. 20, 213. doi:10.1186/s12870-020-2310-0
Sánchez-Pérez, R., Dicenta, F., and Martínez-Gómez, P. (2012). Inheritance of Chilling and Heat Requirements for Flowering in almond and QTL Analysis. Tree Genet. Genomes 8, 379–389. doi:10.1007/s11295-011-0448-5
Schneider, M., Tognolli, M., and Bairoch, A. (2004). The Swiss-Prot Protein Knowledgebase and ExPASy: Providing the Plant Community with High Quality Proteomic Data and Tools. Plant Physiol. Biochem. 42, 1013–1021. doi:10.1016/j.plaphy.2004.10.009
Shang, L., Liu, F., Wang, Y., Abduweli, A., Cai, S., Wang, K., et al. (2015). Dynamic QTL Mapping for Plant Height in Upland Cotton (Gossypium Hirsutum). Plant Breed 134, 703–712. doi:10.1111/pbr.12316
Sierra-Orozco, E., Shekasteband, R., Illa-Berenguer, E., Snouffer, A., van der Knaap, E., Lee, T. G., et al. (2021). Identification and Characterization of GLOBE, a Major Gene Controlling Fruit Shape and Impacting Fruit Size and Marketability in Tomato. Hortic. Res. 8, 138. doi:10.1038/s41438-021-00574-3
Song, X., Xu, Y., Gao, K., Fan, G., Zhang, F., Deng, C., et al. (2020). High-density Genetic Map Construction and Identification of Loci Controlling Flower-type Traits in Chrysanthemum (Chrysanthemum × Morifolium Ramat). Hortic. Res. 7, 108. doi:10.1038/s41438-020-0333-1
Streich, J., Romero, J., Gazolla, J. G. F. M., Kainer, D., Cliff, A., Prates, E. T., et al. (2020). Can Exascale Computing and Explainable Artificial Intelligence Applied to Plant Biology Deliver on the United Nations Sustainable Development Goals? Curr. Opin. Biotechnol. 61, 217–225. doi:10.1016/j.copbio.2020.01.010
Sun, B., Zhan, X.-D., Lin, Z.-C., Wu, W.-X., Yu, P., Zhang, Y.-X., et al. (2017). Fine Mapping and Candidate Gene Analysis of qHD5, a Novel Major QTL with Pleiotropism for Yield-Related Traits in rice (Oryza Sativa L). Theor. Appl. Genet. 130, 247–258. doi:10.1007/s00122-016-2787-y
Sun, J., Wang, J., Liu, H., Xie, D., Zheng, H., Zhao, H., et al. (2015). Dynamic QTL Analysis of rice Seedling Height and Tiller Number under Salt Stress. J. Nucl. Agric. Sci. 29, 235–243. (In Chinese). doi:10.11869/j.issn.100-8551.2015.02.0235
Sun, X., Liu, D., Zhang, X., Li, W., Liu, H., Hong, W., et al. (2013). SLAF-seq: An Efficient Method of Large-Scale De Novo SNP Discovery and Genotyping Using High-Throughput Sequencing. PLoS One 8, e58700. doi:10.1371/journal.pone.0058700
Szymański, J., Bocobza, S., Panda, S., Sonawane, P., Cárdenas, P. D., Lashbrooke, J., et al. (2020). Analysis of Wild Tomato Introgression Lines Elucidates the Genetic Basis of Transcriptome and Metabolome Variation Underlying Fruit Traits and Pathogen Response. Nat. Genet. 52, 1111–1121. doi:10.1038/s41588-020-0690-6
Tello, J., Roux, C., Chouiki, H., Laucou, V., Sarah, G., Weber, A., et al. (2019). A Novel High-Density grapevine (Vitis vinifera L.) Integrated Linkage Map Using GBS in a Half-Diallel Population. Theor. Appl. Genet. 132, 2237–2252. doi:10.1007/s00122-019-03351-y
Tian, Y., Pu, X., Yu, H., Ji, A., Gao, R., Hu, Y., et al. (2020). Genome-wide Characterization and Analysis of bHLH Transcription Factors Related to Crocin Biosynthesis in Gardenia Jasminoides Ellis (Rubiaceae). Biomed. Res. Int. 2020, 1–11. doi:10.1155/2020/2903861
Toonen, R. J., Puritz, J. B., Forsman, Z. H., Whitney, J. L., Fernandez-Silva, I., Andrews, K. R., et al. (2013). ezRAD: a Simplified Method for Genomic Genotyping in Non-model Organisms. PeerJ 1, e203. doi:10.7717/peerj.203
Tsanakas, G. F., Polidoros, A. N., and Economou, A. S. (2013). Genetic Variation in Gardenia Grown as Pot Plant in Greece. Scientia Horticulturae 162, 213–217. doi:10.1016/j.scienta.2013.08.020
Van Ooijen, J. W. (2006). JoinMap ® 4, Software for the Calculation of Genetic Linkage Maps in Experimental Populations. Wageningen, Netherlands: Kyazma BV.
Van Ooijen, J. W. (2009). MapQTL ® 6, Software for the Mapping of Quantitative Trait Loci in Experimental Populations of Diploid Species. Wageningen, Netherlands: Kyazma BV.
van Os, H., Stam, P., Visser, R. G. F., and van Eck, H. J. (2005). SMOOTH: a Statistical Method for Successful Removal of Genotyping Errors from High-Density Genetic Linkage Data. Theor. Appl. Genet. 112, 187–194. doi:10.1007/s00122-005-0124-y
Wang, C. M., Lo, L. C., Zhu, Z. Y., and Yue, G. H. (2006). A Genome Scan for Quantitative Trait Loci Affecting Growth-Related Traits in an F1 Family of Asian Seabass (Lates calcarifer). BMC Genomics 7, 274. doi:10.1186/1471-2164-7-274
Wang, F., Miao, M., Xia, H., Yang, L.-G., Wang, S.-K., and Sun, G.-J. (2017). Antioxidant Activities of Aqueous Extracts from 12 Chinese Edible Flowers In Vitro and In Vivo. Food Nutr. Res. 61, 1265324. doi:10.1080/16546628.2017.1265324
Wang, J., Li, Q., Zhong, X., Song, J., Kong, L., and Yu, H. (2018). An Integrated Genetic Map Based on EST-SNPs and QTL Analysis of Shell Color Traits in Pacific Oyster Crassostrea gigas. Aquaculture 492, 226–236. doi:10.1016/j.aquaculture.2018.04.018
Wang, K., Li, M., and Hakonarson, H. (2010). ANNOVAR: Functional Annotation of Genetic Variants from High-Throughput Sequencing Data. Nucleic Acids Res. 38, e164. doi:10.1093/nar/gkq603
Wang, S., Meyer, E., McKay, J. K., and Matz, M. V. (2012). 2b-RAD: a Simple and Flexible Method for Genome-wide Genotyping. Nat. Methods 9, 808–810. doi:10.1038/nmeth.2023
Wang, X., Wang, H., Long, Y., Liu, L., Zhao, Y., Tian, J., et al. (2015). Dynamic and Comparative QTL Analysis for Plant Height in Different Developmental Stages of Brassica Napus L. Theor. Appl. Genet. 128, 1175–1192. doi:10.1007/s00122-015-2498-9
Wang, X., Zhang, R., Song, W., Han, L., Liu, X., Sun, X., et al. (2019). Dynamic Plant Height QTL Revealed in maize through Remote Sensing Phenotyping Using a High-Throughput Unmanned Aerial Vehicle (UAV). Sci. Rep. 9, 3458. doi:10.1038/s41598-019-39448-z
Wang, Y., Wang, C., Han, H., Luo, Y., Wang, Z., Yan, C., et al. (2020). Construction of a High-Density Genetic Map and Analysis of Seed-Related Traits Using Specific Length Amplified Fragment Sequencing for Cucurbita Maxima. Front. Plant Sci. 10, 1782. doi:10.3389/fpls.2019.01782
Webb, A., Cottage, A., Wood, T., Khamassi, K., Hobbs, D., Gostkiewicz, K., et al. (2016). A SNP-Based Consensus Genetic Map for Synteny-Based Trait Targeting in Faba Bean (Vicia faba L). Plant Biotechnol. J. 14, 177–185. doi:10.1111/pbi.12371
Wei, J., Man, Q., Ding, C., Hu, Y., Liu, M., Li, H., et al. (2019). Proteomic Investigations of Transcription Factors Critical in Geniposide-Mediated Suppression of Alcoholic Steatosis and in Overdose-Induced Hepatotoxicity on Liver in Rats. J. Proteome Res. 18, 3821–3830. doi:10.1021/acs.jproteome.9b00140
Wen, J., Jiang, F., Weng, Y., Sun, M., Shi, X., Zhou, Y., et al. (2019). Identification of Heat-Tolerance QTLs and High-Temperature Stress-Responsive Genes through Conventional QTL Mapping, QTL-Seq and RNA-Seq in Tomato. BMC Plant Biol. 19, 398. doi:10.1186/s12870-019-2008-3
Xu, Y. Q., Wei, G. Y., Zhou, Y., Ge, F., and Luo, G. M. (2014). Isolation and Characterization of Twenty-Two Polymorphic Microsatellite Markers from Gardenia Jasminoides (Rubiaceae). J. Genet. 94, e22–e24. doi:10.1007/s12041-014-0348-1
Xu, Z., Pu, X., Gao, R., Demurtas, O. C., Fleck, S. J., Richter, M., et al. (2020). Tandem Gene Duplications Drive Divergent Evolution of Caffeine and Crocin Biosynthetic Pathways in Plants. BMC Biol. 18, 63. doi:10.1186/s12915-020-00795-3
Yamakawa, H., Haque, E., Tanaka, M., Takagi, H., Asano, K., Shimosaka, E., et al. (2021). Polyploid QTL‐seq towards Rapid Development of Tightly Linked DNA Markers for Potato and Sweetpotato Breeding through Whole‐genome Resequencing. Plant Biotechnol. J. 19, 2040–2051. doi:10.1111/pbi.13633
Yan, H. H., Mudge, J., Kim, D.-J., Shoemaker, R. C., Cook, D. R., and Young, N. D. (2004). Comparative Physical Mapping Reveals Features of Microsynteny betweenGlycine max,Medicago Truncatula, andArabidopsis Thaliana. Genome 47, 141–155. doi:10.1139/g03-106
Yang, W., Wang, Y., Jiang, D., Tian, C., Zhu, C., Li, G., et al. (2020). ddRADseq-Assisted Construction of a High-Density SNP Genetic Map and QTL fine Mapping for Growth-Related Traits in the Spotted Scat (Scatophagus argus). BMC genomics 21, 278. doi:10.1186/s12864-020-6658-1
Yu, C., Yan, C., Liu, Y., Liu, Y., Jia, Y., Lavelle, D., et al. (2020). Upregulation of a KN1 Homolog by Transposon Insertion Promotes Leafy Head Development in Lettuce. Proc. Natl. Acad. Sci. USA 117, 33668–33678. doi:10.1073/pnas.2019698117
Zhang, L., Guo, D., Guo, L., Guo, Q., Wang, H., and Hou, X. (2019). Construction of a High-Density Genetic Map and QTLs Mapping with GBS from the Interspecific F1 Population of P. Ostii 'Fengdan Bai' and P. Suffruticosa 'Xin Riyuejin'. Scientia Horticulturae 246, 190–200. doi:10.1016/j.scienta.2018.10.039
Zhang, Z., Wei, T., Zhong, Y., Li, X., and Huang, J. (2016). Construction of a High-Density Genetic Map of Ziziphus Jujuba Mill. Using Genotyping by Sequencing Technology. Tree Genet. Genomes 12, 76. doi:10.1007/s11295-016-1032-9
Zhao, J., Li, H., Xu, Y., Yin, Y., Huang, T., Zhang, B., et al. (2021). A Consensus and Saturated Genetic Map Provides Insight into Genome Anchoring, Synteny of Solanaceae and Leaf- and Fruit-Related QTLs in Wolfberry (Lycium Linn). BMC Plant Biol. 21, 350. doi:10.1186/s12870-021-03115-1
Keywords: genetic map, genotyping-by-sequencing, growth-and leaf-related traits, QTL, synteny, gardenia
Citation: Cui Y, Fan B, Xu X, Sheng S, Xu Y and Wang X (2022) A High-Density Genetic Map Enables Genome Synteny and QTL Mapping of Vegetative Growth and Leaf Traits in Gardenia. Front. Genet. 12:802738. doi: 10.3389/fgene.2021.802738
Received: 27 October 2021; Accepted: 13 December 2021;
Published: 04 January 2022.
Edited by:
Jinyoung Y. Barnaby, Agricultural Research Service (USDA), United StatesReviewed by:
Chengsong Zhu, University of Texas Southwestern Medical Center, United StatesXingbo Wu, Agricultural Research Service (USDA), United States
Copyright © 2022 Cui, Fan, Xu, Sheng, Xu and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yuhui Xu, genetics_2010@163.com; Xiaoyun Wang, wxy20052002@aliyun.com