Combining QTL mapping and RNA-Seq Unravels candidate genes for Alfalfa (Medicago sativa L.) leaf development
BMC Plant Biology volume 22, Article number: 485 (2022)
Leaf size affects crop canopy morphology and photosynthetic efficiency, which can influence forage yield and quality. It is of great significance to mine the key genes controlling leaf development for breeding new alfalfa varieties. In this study, we mapped leaf length (LL), leaf width (LW), and leaf area (LA) in an F1 mapping population derived from a cultivar named ZhongmuNo.1 with larger leaf area and a landrace named Cangzhou with smaller leaf area.
This study showed that the larger LW was more conducive to increasing LA. A total of 24 significant quantitative trait loci (QTL) associated with leaf size were identified on both the paternal and maternal linkage maps. Among them, nine QTL explained about 11.50–22.45% phenotypic variation. RNA-seq analysis identified 2,443 leaf-specific genes and 3,770 differentially expressed genes. Combining QTL mapping, RNA-seq alalysis, and qRT-PCR, we identified seven candidate genes associated with leaf development in five major QTL regions.
Our study will provide a theoretical basis for marker-assisted breeding and lay a foundation for further revealing molecular mechanism of leaf development in alfalfa.
Photosynthesis has become a major goal for breeders aiming to further improve yields in field crops , which is directly affected plant growth and grain yield . Leaves are the main photosynthetic organs in plants. It has been reported that leaf size has a significant effect on photosynthetic rate, transpiration, and carbon fixation, which contribute to plant biomass [3,4,5]. Therefore, it is feasible to produce new cultivars with high yield by optimizing the leaf-related traits of crop plants.
Given its obvious effect on plant growth, the genetic mechanism that control leaf size has been investigated intensively in model plants [6, 7]. The leaf size and growth are controlled by several regulators, including transcriptional factors, phytohormones and microRNAs . In Arabidopsis, phytohormones (auxin, gibberellin, cytokinin, and brassinosteroid) biosynthesis or signaling genes have been cloned and identified in leaf area regulation . For instance, the disturbance to auxin homeostasis and signaling achieved by over‐expressing PID resulted in plants which produce plants exhibiting smaller than normal leaves . Meanwhile, QTL and candidate genes associated with leaf size have been reported in grass plant species, such as maize , rice [10, 11], wheat , and barley . Using a recombinant inbred line population with a single-nucleotide polymorphism (SNP) genetic linkage map, a major QTL for leaf width of maize was further fine-mapped to an interval of 55 kb and two candidate genes were identified in developing leaves . In rice, enhanced expression of a QTL related to leaf length, qLL9, significantly increased the grain size, thousand-grain weight, and grain yield . Also, legumes, such as soybean , common bean , and peanuts  have been investigated for the genetic basis of variation in leaf size. In pea, 37 bulked segregant analysis markers associated with leaf traits were identified by using QTL mapping combined with BSA .
Alfalfa (Medicago sativa L.), a perennial forage legume, has been cultivated for hay, pasture, and silage worldwide [20, 21]. One of the primary goals of alfalfa breeding programs is to produce new cultivars with high biomass, as the primary commercial determinant of alfalfa production is its above-ground biomass. Therefore, it is feasible to improve photosynthetic efficiency and yield by optimizing the leaf-related traits of alfalfa.
Autotetraploid alfalfa, a self-incompatibly cross-pollinated crop, have a very complex genome . This complex genetic nature, eg the complex loci segregate, hinders the construction of genetic linkage maps to improve important agricultural traits in alfalfa. The strategy of using limited biallelic markers segregating in specific patterns are suitable for constructing genetic linkage maps of the tetraploid crop. Single dose allele (SDA) markers unique to each parent with segregation of 1:1 (Aaaa x aaaa) can be used to construct genetic linkage maps for an F1 populations, which was derived from two heterozygous alfalfa genotypes . As such, JoinMap software, which was designed for diploid species, were used to construct linkage maps in tetraploid alfalfa [23,24,25]. Also, there have been reported on alfalfa QTL associated with important agronomic traits, including biomass-related traits [25, 26], spring vigor , as well as stress-related traits including winter hardness , and leaf rust resistance .
The genetic and genomic basis of leaf size has been investigated extensively in cereals crops, whereas there is scarce on herbaceous perennials, such as alfalfa. Although a number of QTL have been identified in alfalfa, the genetic architecture and main-effect QTL/genes related to leaf development remain unclear. Integration of RNA-seq and QTL mapping had been considered to be a reliable method to rapidly identify potential candidate genes. Although it has been applied in rice and maize [30, 31]. The strategy has not been reported in alfalfa leaf development related genes. The objectives of this study were: (1) to identify QTL associated with alfalfa leaf development, and (2) to locate potential candidate genes according to QTL mapping, differentially expressed genes (DEGs) analysis, and functional annotation. The identified QTL and potential candidate genes would be helpful for alfalfa breeding program.
Materials and methods
Mapping population and Field location were as described previously . Briefly, the the Medicago. sativa subsp. sativa landrace (Cangzhou) with smaller leaf area was used as paternal parent (P1), and the Medicago. sativa subsp. sativa cultivar (ZhongmuNo.1) with larger leaves was used as maternal parent (P2). In 2018, 388 of the original progenies (F1 population) generated by crossing survived and under phenotypic investigation for this study. Both the parents and the F1 were grown at two locations; Langfang (LF), Hebei Province, China (39.59°N, 116.59°E) and Changping (CP), Beijing, China (40.18°N, 116.24°E), using a randomized complete block design with three replications.
Phenotypic data collection and analysis
Phenotypic data were collected at the two locations (CP and LF) in the year of 2018 (CP18 and LF18) and 2019 (CP19 and LF19). In order to ensure the consistency of sampling at different locations in different years, leaves of all genotypes were sampled at the initial flowering stage (the first flower appeared). The leaf from the middle leaflet of the third or fourth fully expanded trifoliate leaf from the shoot tip of each plant were chosen for sampling. For each plant, three leaves were sampled. Leaf size related phenotype, leaf length (LL), leaf width (LW), and leaf area (LA) were measured using a LAM-B type handheld leaf area meter (Shijiazhuang Fansheng Technology Co.; Ltd). Best linear unbiased predictions (BLUPs) were analyzed of the phenotypic data using PROC MIXED . Interaction of a given genotype with different environments was estimated using PROC GLM (SAS Institute, 2010). Correlation of leaf-related traits was analyzed using the R package R/corrplot. The broad-sense heritability (H2) was calculated using the AOV (Analysis of variance) function of QTL IciMapping .
Genetic linkage map construction and QTL mapping
The methods used for linkage map construction and QTL mapping were described previously . Briefly, DNA samples of the progenies and both parents were extracted according to the manufacturer’s protocol from ~ 100 mg of fresh leaf tissues using the CWBIO Plant Genomic DNA Kit (CoWin Biosciences, Beijing, China). Single dose allele (SDA) markers with a segregation ratio of less than 2:1 among F1 progenies were used to construct a genetic linkage map using JoinMap 4.0 according to a method published previously .
QTL detection was performed by inclusive composite interval mapping (ICIM) using the QTL IciMapping Software (version 4.1) [34, 35]. For the phenotypic data in single environment and the BLUP value, QTL were analyzed using the ICIM-ADD mapping method with the default parameters (PIN = 0.001, Walking step = 1.0 cM, type I error = 0.05), and the significance threshold values (P ≤ 0.05) were determined through 1000 permutation tests. The result of QTL analysis was obtained and imported to the software MapChart to display results .
RNA-seq data analysis
In previous study, RNA-Seq project of leaves of two different Medicago sativa genotype (Medicago sativa subsp. sativa cv. Beaver and M. sativa subsp. falcata accession PI 641,381) was performed . Significant difference was noted between genotypes in terms of leaf size, while they had no difference in plant height and biomass. The first expanded leaflet was harvested from four biological replicates of ‘sativa’ and ‘falcata’, respectively. Also, a RNA-seq project of alfalfa different tissues (leaves, roots, nodules, flowers, elongating internodes, and post-elongation internodes) was performed in previous . Leaves were harvested 28 days after the second cutback. In this study, to identify leaf-specific expression genes and DEGs between different leaf size genotypes, the two RNA-seq projects were download. RNA-Seq data quality was estimated by FastQC v0.11.9 (http://bioinformatics.bbsrc.ac.uk/projects/fastqc), and low quality sequences were filtered using fastp v0.12.4 . Then, clean reads were mapped to the alfalfa reference genome using hisat2 v4.8.2 , and Samtools was used to translate SAM file to BAM file and then to sort BAM files using default parameters . FeatureCounts v2.0.1 was used to generate read counts for each sample .
Leaf-specific genes retrieval
We used the R package TissueEnrich to identify leaf-specific genes . The genes were divided into six groups based on their gene expression across the tissues in this package. In this research, leaf-specific expression genes include two groups (leaf-enriched and group-enriched). Leaf-enriched genes had an expression level in leaf tissue greater than or equal to 1 (FPKM), and also have at least five-fold higher expression levels in leaf tissue compared to all other tissues. Group-enriched genes with a FPKM value ≥ 1 that also have at least five-fold higher expression levels in leaf and another tissue compared to all other tissues, and that were not considered as leaf-enriched.
Differentially expressed genes identification
DEGs, defined as |log2 fold change|≥ 5.259 and an adjusted P value of < 0.01, were screened using the R package DESeq2 .
Potential candidate genes identification
Local blast search was conducted using TBtools software . Briefly, based on the flanking marker sequences, the physical positions of the detected QTL were located on alfalfa genome (cultivar XingJiangDaYe) . Within the physical region, sequences encoding protein products were extracted. Leaf-specific genes/DEGs within QTL intervals were annotated based on BLSATP search in NCBI (https://www.ncbi.nlm.nih.gov/) and Ensembl (https://ensembl.gramene.org/).
Quantitative real-time PCR (qRT-PCR) were used to determine whether candidate genes were differentially expressed in the two parents. qRT-PCR was conducted using unexpanded trifoliate leaves from three biological replicates of each parent. Six accessions with long leaf length and six accessions with short leaf length were selected to identify the relative expression levels of MS.gene07851 in germplasms with different leaf length. These materials were grown in Langfang and leaf-related phenotypes were collected in 2018, 2019, and 2020. The information of these accessions was in placed in Supplementary Table 1. On a 7500 Real-Time PCR System (Applied Biosystems, CA, USA), qRT-PCR was implemented using the SYBR Premix Ex Taq (Takara, Japan). The relative gene expression level was calculated by the 2-△△Ct method.
In order to determine the relationship between polymorphic sites of the candidate genes and leaf development, variations were detected in each parent by comparing the resequencing data to the XingJiangDaYe reference genome. The process of obtaining BAM files for each genotype was consistent with RNA-seq data. GATK4.0 was used to perform variant calling , and the Heterozygous variants were filtered out. The resequencing data of the two parents have been deposited in the NCBI Sequence Read Archive under BioProject accession number PRJNA861857.
For F1 population and their parental plants, leaf size related traits including leaf length (LL), leaf width (LW), and leaf area (LA) in the tested environments were summarized. The average of leaf size of the paternal parent was significantly (P < 0.01) smaller than that of the maternal parent in the tested environments (Table 1). In the F1 population, three traits exhibited continuous distributions with transgressive segregations (Fig. 1, Table 1), suggesting that QTL analysis was feasible to dissect the leaf size components in alfalfa. The frequency distribution histogram of each trait generally conformed to a normal distribution in all environments (Fig. 1). Pearson’s correlation coefficient indicated highly significant (P < 0.001) degrees of correlation between two of the three traits at both locations (Fig. 1). The correlation coefficient was 0.93 between LA and LW, and 0.89 between LA and LL, suggesting that the LW contributes more to LA than LL does. The effects of genotype (G), environment (E) and their interaction (G x E) were highly significant (P < 0.001) for each trait analyzed. (Table S1). The average broad-sense heritability (H2) for LL, LW, and LA was 60%, 63%, and 64%, respectively, indicating that genetic factors played an important role in determining these traits (Table S1).
QTL analysis of leaf-related traits
Within the 32 homologs of the 8 alfalfa chromosomes, we detected 24 significant (P ≤ 0.05) QTL for leaf-related traits on both male and female linkage maps (Table 2). P1 contributed 4 of the 24 QTL, and 20 QTL were contributed by P2. Thirteen QTL were identified at the location of CP, seven at the location of LF, and four QTL were identified using BLUPs. These QTL contributed between 2.37%-22.45% of the phenotypic variance at different locations.
QTL associated with leaf length. Four significant QTL (qLL-2B, qLL-6D-1, qLL-6D-2, and qLL-6D-3) associated with LL were mapped on the maternal linkage maps. The QTL were detected on linkage group (LG) 2B, and 6D (Figs. S1 and S2). The phenotypic variation explained (PVE) by a single QTL ranged from 3.95% to 6.88% in different environments (Table 2). Among them, the QTL (qLL-2B) was on LG 2B at the position of 40.5–42.5 cM with PVE value of 5.38%. The rest three QTL (qLL-6D-1, qLL-6D-2, and qLL-6D-3), located on LG 6D of P2 from the 38.5 cM to 40.5 cM position, were considered to be the same QTL.
QTL associated with leaf width. Ten QTL for LW were located on homologs 2C (1 QTL), 6B (1 QTL), 6D (5 QTL), 7D (1 QTL), 8A (1 QTL), and 8B (1 QTL) with an individual QTL accounting for 2.37–22.45% of phenotypic variations (Table 2, Figs. S1 and S2). Among them, two QTL were identified in P1 and 8 QTL in P2 with five QTL (qLW-6D-1, qLW-6D-2, qLW-6D-3, qLW-6D-4, and qLW-6D-5) located on chromosome 6D. These five QTL explained 11.50%-22.45% of the phenotypic variation for LW.
QTL associated with leaf area. For LA, 10 significant QTL were identified on LG 4A, 6D, 7D, and 8C. These QTL explained 3.49 to 18.88% of the phenotypic variation in different environments (Table 2, Figs. S1 and S2). The LSI (1-LOD support interval in cM unit) of these QTL was from 1 to 7 cM. Five QTL (qLA-6D-1, qLA-6D-2, qLA-6D-3, qLA-6D-4, and qLA-6D-5) were co-localized on chromosome 6D in P2 with higher phenotypic variation than the most other QTL explained averagely for 7.01% for LA. Among them, the QTL qLA-6D-5 explained the highest (18.88%) phenotypic variance for LA, implying a higher reliability of the loci on 6D for leaf area in alfalfa.
Overlapped QTL intervals. Four overlapped QTL intervals, which were identified at the same position in different environments/phenotype, were detected. All overlapped intervals were identified in P2 parent on chromosome 6D. qLL-6D-1, qLL-6D-2, and qLL-6D-3 located in 38.5–40.5 cM on chromosome 6D. qLW-6D-2 and qLA-6D-1 were identified at CP and LF, respectively. And they were overlapped with qLA-6D-2. qLW-6D-3, qLW-6D-4, and qLA-6D-3 were located at the position of 28 cM with an average PVE value of 14.53%. The last overlapped interval (qLW-6D-5 and qLA-6D-4) was at the position of 34 cM with an average PVE value of 13.95%.
We focused on six major QTL intervals, including four overlapped QTL intervals and two QTLs with high PVE values (≥ 10%). The physical interval of the co-located QTL, including qLW-6D-5 and qLA-6D-4, was really large (25.1 ~ 74.6 Mb). So we selected ± 2 Mb physical intervals of two flanking markers to predict candidate genes. Based on local-BLAST, a total of 1,573 genes were identified in seven physical regions, which were used for further analysis in our study.
Leaf-specific genes within QTL regions
A total of 2,443 leaf-specific genes were identified, including 1,294 leaf-enriched and 1,149 group-enriched (Fig. 2A). Among the significantly KEGG enriched pathways, metabolism and biosynthesis pathway contained more genes than the others, indicating that these metabolism and biosynthesis pathway might have a modulating effect on the function and development of alfalfa leaves (Fig. 2B). Leaf-specific analysis and QTL mapping revealed 29 common genes, and four of them were predicted as candidates based on BLAST-P (Fig. 2C and 2D). These genes are involved in cell proliferation and growth, and their homologous have been reported to play a role in leaf development. The expression levels of these four genes were further investigated in the young leaves of the two parents. The result showed that MS.gene08405 and MS.gene05412 were differently expressed in the two parents (Fig. 3). They were located in two different QTL intervals of qLA-7D-1 (MS.gene08405), qLW-6D-5 and qLA-6D-4 (MS.gene05412) (Table S5).
DEGs within QTL regions
Compared to ‘sativa’, leaf size was 43.3% relative decrease in ‘falcata’. DEG analysis identified 2,061 up-regulated and 1,709 down-regulated DEGs (‘sativa’ Vs ‘falcata’) in leaf tissue of the two genotype (Fig. 4A and 4B; Table S3; Table S4). A total of 41 DEGs located within the QTL physical intervals (Fig. 4C). And nine DEGs, whose homologous genes were reported to have functions in cell proliferation and cell expansion, were annotated as candidates associated with leaf development in plants. Among them, five genes, which had different expression patterns in two parents, were identified as candidate genes (Fig. 3). The five candidates were located in four intervals, including: qLW-6D-1 (MS.gene54171 and MS.gene80652), qLW-6D-5, qLA-6D-4 (MS.gene05433), qLW-6D-3, qLW-6D-4, and qLA-6D-3 (MS.gene44488), qLL-6D-1, qLL-6D-2, qLL-6D-3 (MS.gene07851) (Table S5). MS.gene07851, which were located in three LL-related QTLs, had a significant expression difference (P < 0.0001) between the two parents. In addition, long leaf length accessions displayed high expression levels of MS.gene07851, which was in accordance with their LL values (Fig. 5). Thus, MS.gene07851 might be the causal gene that underlies variation in the leaf length.
Variation of candidate genes between the parental genotypes
We investigated the polymorphic sites among the seven candidate genes by comparing their whole-genome re-sequencing reads with the reference genome. By mapping to the XingJiangDaYe genome, 26 SNPs and 3 InDels were detected among the seven candidate genes between two genotypes (Table S5). Among them, five SNPs caused amino acid changes in three candidate genes of MS.gene05433 (three), MS.gene44488 (one), and MS.gene07851.t1 (one). Notably, the gene MS.gene07851 showed a SNP (C/G) variant in the two parental genotypes, which may correspond to its different function in leaf development. These results revealed multiple variants in the seven candidate genes, suggesting that the variants may contribute to different protein functions and differences in leaf size between the two parents.
Leaf size is an important target to increase yield in plant breeding because it affects crop canopy morphology and photosynthetic efficiency . For forage crop, alfalfa, increasing leaf size is an important approach to improve its yield and quality. However, the genetic mechanisms of alfalfa are intricate due to its homologous polyploidy, cross-pollinated, and highly heterozygous in nature . To date, there have been less reports regarding leaf development genes in alfalfa.
RNA-sequencing (RNA-seq), a high resolution technology, can provide far more massive amounts of information regarding the gene expression levels of various species in different conditions [48, 49]. Based on quantifications of expressed genes, differentially expressed genes are a group of genes which have a statistically significant difference at the expression level. DEG analyses can provide reliable insight into the genetic mechanisms in crops that have different phenotype, including plant at different development stages and in different growing conditions [50, 51]. The tissue-specific genes refer to genes specifically expressed in one or several tissues/cell types, which involve regulate tissues-specific morphological structures or physiological functions. Identification of these genes can help us better understanding of tissue-gene relationship. In alfalfa, RNA-seq analysis uncovered molecular mechanisms underlying abiotic stress, including salt [52, 53], drought [54, 55], and cold tolerance . Meanwhile, pathways involved in early flowering  and differentially expressed genes related to defoliation traits  had also been identified in alfalfa by transcriptome analysis. In the mapping population, the two parents not only had a significant difference in leaf size but also in other traits such as plant height and yield . ‘Beaver’ and ‘PI 641,381’, appeared morphologically similar overall, and no significant differences were noted between genotypes in terms of plant height, aboveground biomass, or root length. ‘PI 641,381’ had a significant reduction in leaf area compared to 'Beaver' .
QTL mapping, a powerful method to identify trait associated QTLs/genes, has been reports on alfalfa important agronomic traits [25, 26, 60]. Here, we mapped 24 QTL for leaf size using the high density genetic maps. The effects of environment were highly significant (P < 0.001) for each trait. No QTL was identified at the same genetic position by using the data from different cultivated condition. Most of the detected QTLs for LA, LL and LW were different. Only three QTL intervals were detected for both LW and LA (Table 2). In this study, the genetic linkage maps were constructed by GBS-SNP markers. Although, using GBS-SNP mapping greatly facilitates genetic dissection of quantitative traits in alfalfa. There are some limitations of this strategy. GBS missed a lot of genotypic information, and many polymorphic markers not classified as SDA were not used for map construction, resulting in reduced power of QTL detection .
The traditional QTL fine mapping and map-based cloning in alfalfa are limited by constructing near-isogenic lines. So far, QTL mapping studies related to alfalfa important agronomic traits were mainly focused on the QTL identification . Few studies have made an effort to identify candidate genes within the QTL ranges. Transcriptome analysis can be used to characterize gene expression levels in various species. In recent years, RNA-seq has been used in QTL mapping studies to rapidly identify potential candidate genes in crops [30, 31, 62]. The strategy of QTL mapping combined with RNA-seq allowed for rapid detection of candidate genes in target regions. A total of seven candidates associated with leaf development were finally detected in our study.
The processes that mainly drive leaf development are cell proliferation and cell expansion, and numerous genes have been identified in Arabidopsis, such as CYP78A5, DA1, and MYBH [7, 63, 64]. CYP78A5, a member of cytochrome P450 protein family, which played an important role in cell proliferation and expansion , was a homologous gene of one candidate (MS.gene08405). It was annotated with the function of stimulating cell proliferation during leaf development . Interestingly, MS.gene07851 was annotated as protein tyrosine and serine/threonine kinase, which contained conserved domain PK_Tyr_Ser-Thr. Protein kinase has conserved function in several cellular processes, including division, proliferation, apoptosis, and differentiation [66,67,68]. The family members have been reported to regulate leaf development in Arabidopsis [68, 69]. SNPs are the major determinant of phenotypic differences . Thus, we can hypothesize that the identified SNP in MS.gene07851, may affect their expression and function. At the same time, the results of qRT-PCR showed that the expression levels of MS.gene07851 in long-leaf germplasms were higher than that in short-leaf germplasms, indicating that it may be involved in the regulation of leaf length development. The candidate genes identified in our study should be confirmed by genetic transformation and further studies.
In summary, QTL mapping and transcriptome analysis were performed to identify the genes for leaf development in an F1 population derived from a larger leaf area cultivar named Zhong mu No.1 and a smaller leaf area accession named Cangzhou. Six important QTL regions were identified, which contained seven candidate genes. This study provides valuable resources for enhancing breeding programs aimed at improvement of leaf development in alfalfa.
Availability of data and materials
The raw data of GBS were submitted to the NCBI Sequence Read Archive with bioproject ID: PRJNA522887 (https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA522887). RNA-seq raw data used in this study was downloaded from the SRA database in NCBI with the Bioproject accession numbers of PRJNA765383 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA765383) and PRJNA276155 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA276155). The Resequencing data have been deposited in the NCBI with the BioProject accession number PRJNA861857 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA861857).
Best linear unbiased predictions
Differentially expressed genes
Genotyping by sequencing
- G x E:
Genotype by environment
Composite interval mapping with additive effect
Logarithm of odds for each QTL
Phenotypic variation explained
Quantitative trait loci
Single dose allele
Single nucleotide polymorphism
Wu A, Hammer GL, Doherty A, von Caemmerer S, Farquhar GD. Quantifying impacts of enhancing photosynthesis on crop yield. Nat Plants. 2019;5(4):380–8.
Zelitch I. The close relationship between net photosynthesis and crop yield. Bioscience. 1982;32(10):796–802.
Walter A, Silk WK, Schurr U. Environmental effects on spatial and temporal patterns of leaf and root growth. Annu Rev Plant Biol. 2009;60:279–304.
Giuliani R, Koteyeva N, Voznesenskaya E, Evans MA, Cousins AB, Edwards GE. Coordination of Leaf Photosynthesis, Transpiration, and Structural Traits in Rice and Wild Relatives (Genus Oryza). Plant Physiol. 2013;162(3):1632–51.
Xu P, Ali A, Han B, Wu X. Current Advances in Molecular Basis and Mechanisms Regulating Leaf Morphology in Rice. Front Plant Sci. 2018;9:1528.
Wang H, Kong F, Zhou C. From genes to networks: The genetic control of leaf development. J Integr J Integr Plant Biol. 2021;63(7):1181–96.
Vercruysse J, Baekelandt A, Gonzalez N, Inzé D. Molecular networks regulating cell division during Arabidopsis leaf growth. J Exp Bot. 2020;71(8):2365–78.
Saini K, Markakis MN, Zdanio M, Balcerowicz DM, Beeckman T, De Veylder L, Prinsen E, Beemster GTS, Vissenberg K. Alteration in Auxin Homeostasis and Signaling by Overexpression Of PINOID Kinase Causes Leaf Growth Defects in Arabidopsis thaliana. Front Plant Sci. 2017;8:1009.
Fu Y, Xu G, Chen H, Wang X, Chen Q, Huang C, Li D, Xu D, Tian J, Wu W. QTL mapping for leaf morphology traits in a large maize-teosinte population. Mol Breed. 2019;39:103.
Wu J, Qi Y, Hu G, Li J, Li Z, Zhang H. Genetic architecture of flag leaf length and width in rice (Oryza sativa L.) revealed by association mapping. Genes & Genomics. 2017;39:341–52.
Tang X, Gong R, Sun W, Zhang C, Yu S. Genetic dissection and validation of candidate genes for flag leaf size in rice (Oryza sativa L.). Theor Appl Genet. 2018;131(4):801–15.
Zhao C, Bao Y, Wang X, Yu H, Ding A, Guan C, Cui J, Wu Y, Sun H, Li X. QTL for flag leaf size and their influence on yield-related traits in wheat (Triticum aestivum L.). Mol Breed. 2015;35:24.
Liu L, Sun G, Ren X, Li C, Sun D. Identification of QTL underlying physiological and morphological traits of flag leaf in barley. BMC Genet. 2015;16:29.
Wang B, Zhu Y, Zhu J, Liu Z, Liu H, Dong X, Guo J, Li W, Chen J, Gao C, Zheng X, E L, Lai J, Zhao H, Song W. Identification and Fine-Mapping of a Major Maize Leaf Width QTL in a Re-sequenced Large Recombinant Inbred Lines Population. Front Plant Sci. 2018;9:101.
Fu X, Xu J, Zhou M, Chen M, Shen L, Li T, Zhu Y, Wang J, Hu J, Zhu L, Gao Z, Dong G, Guo L, Ren D, Chen G, Lin J, Qian Q, Zhang G. Enhanced Expression of QTL qLL9/DEP1 Facilitates the Improvement of Leaf Morphology and Grain Yield in Rice. Int J Mol Sci. 2019;20(4):866.
Wang L, Cheng Y, Ma Q, Mu Y, Huang Z, Xia Q, Zhang G, Nian H. QTL fine-mapping of soybean (Glycine max L.) leaf type associated traits in two RILs populations. BMC Genomics. 2019;20(1):260.
Keller B, Manzanares C, Jara C, Lobaton JD, Studer B, Raatz B. Fine-mapping of a major QTL controlling angular leaf spot resistance in common bean (Phaseolus vulgaris L.). Theor Appl Genet. 2015;128(5):813–26.
Han S, Yuan M, Clevenger JP, Li C, Hagan A, Zhang X, Chen C, He G. A SNP-Based Linkage Map Revealed QTLs for Resistance to Early and Late Leaf Spot Diseases in Peanut (Arachis hypogaea L.). Front Plant Sci. 2018;9:1012.
Zheng Y, Xu F, Li Q, Wang G, Liu N, Gong Y, Li L, Chen ZH, Xu S. QTL Mapping Combined With Bulked Segregant Analysis Identify SNP Markers Linked to Leaf Shape Traits in Pisum sativum Using SLAF Sequencing. Front Genet. 2018;9:615.
Adhikari L, Missaoui AM. Nodulation response to molybdenum supplementation in alfalfa and its correlation with root and shoot growth in low pH soil. J Plant Nutr. 2017;40(16):2290–302.
Li X, Brummer EC. Applied genetics and genomics in alfalfa breeding. Agron. 2012;2(1):40–61.
Chen H, Zeng Y, Yang Y, Huang L, Tang B, Zhang H, Hao F, Liu W, Li Y, Liu Y, Zhang X, Zhang R, Zhang Y, Li Y, Wang K, He H, Wang Z, Fan G, Yang H, Bao A, Shang Z, Chen J, Wang W, Qiu Q. Allele-aware chromosome-level genome assembly and efficient transgene-free genome editing for the autotetraploid cultivated alfalfa. Nat Commun. 2020;11(1):2494.
Li X, Wei Y, Acharya A, Jiang Q, Kang J, Brummer EC. A saturated genetic linkage map of autotetraploid alfalfa (Medicago sativa L.) developed using genotyping-by-sequencing is highly syntenous with the Medicago truncatula genome. G3 (Bethesda). 2014;4(10):1971–9.
Ooijen VAN, JW. Multipoint maximum likelihood mapping in a full-sib family of an outbreeding species. Genet Res (Camb). 2011;93(5):343–9.
Zhang F, Kang J, Long R, Yu LX, Wang Z, Zhao Z, Zhang T, Yang Q. High-density linkage map construction and mapping QTL for yield and yield components in autotetraploid alfalfa using RAD-seq. BMC Plant Biol. 2019;19(1):165.
Adhikari L, Makaju SO, Missaoui AM. QTL mapping of flowering time and biomass yield in tetraploid alfalfa (Medicago sativa L.). BMC Plant Biol. 2019;19(1):359.
Mccord P, Gordon V, Saha G, Hellinga J, Vandemark G, Larsen R, Smith M, Miller D. Detection of QTL for forage yield lodging resistance and spring vigor traits in alfalfa (Medicago sativa L.). Euphytica. 2014;200:269–79.
Adhikari L, Lindstrom OM, Markham J, Missaoui AM. Dissecting Key Adaptation Traits in the Polyploid Perennial Medicago sativa Using GBS-SNP Mapping. Front Plant Sci. 2018;9:934.
Adhikari L, Missaoui AM. Quantitative trait loci mapping of leaf rust resistance in tetraploid alfalfa. Physiol Mol Plant Pathol. 2019;106:238–45.
Han Q, Zhu Q, Shen Y, Lee M, Lübberstedt T, Zhao G. QTL Mapping Low-Temperature Germination Ability in the Maize IBM Syn10 DH Population. Plants (Basel). 2022;11(2):214.
Lei L, Zheng H, Bi Y, Yang L, Liu H, Wang J, Sun J, Zhao H, Li X, Li J, Lai Y, Zou D. Identification of a Major QTL and Candidate Gene Analysis of Salt Tolerance at the Bud Burst Stage in Rice (Oryza sativa L.) Using QTL-Seq and RNA-Seq. Rice (N Y). 2020;13(1):55.
Zhang F, Kang J, Long R, Yu LX, Sun Y, Wang Z, Zhao Z, Zhang T, Yang Q. Construction of high-density genetic linkage map and mapping quantitative trait loci (QTL) for flowering time in autotetraploid alfalfa (Medicago sativa L.) using genotyping by sequencing. Plant Genome. 2020;13(3):e20045.
Wolfinger RD, Federer WT, Cordero-Brana O. Recovering information in augmented designs using SAS PROC GLM and PROC MIXED. Agron J. 1997;89:856–9.
Meng L, Li H, Zhang L, Wang J. QTL IciMapping: Integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental populations. Crop J. 2015;3(3):269–83.
Li S, Wang J, Zhang L. Inclusive Composite Interval Mapping of QTL by Environment Interactions in Biparental Populations. PLoS ONE. 2015;10(7):e0132414.
Voorrips RE. MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered. 2002;93(1):77–8.
Singer SD, Subedi U, Lehmann M, Burton Hughes K, Feyissa BA, Hannoufa A, Shan B, Chen G, Kader K, Ortega Polo R, Schwinghamer T, Kaur Dhariwal G, Acharya S. Identification of Differential Drought Response Mechanisms in Medicago sativa subsp. sativa and falcata through Comparative Assessments at the Physiological, Biochemical, and Transcriptional Levels. Plants (Basel). 2021;10(10):2107.
O’Rourke JA, Fu F, Bucciarelli B, Yang SS, Samac DA, Lamb JF, Monteros MJ, Graham MA, Gronwald JW, Krom N, Li J, Dai X, Zhao PX, Vance CP. The Medicago sativa gene index 1.2: a web-accessible gene expression atlas for investigating expression differences between Medicago sativa subspecies. BMC Genomics. 2015;16(1):502.
Chen S, Zhou Y, Chen Y, Gu J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34(17):i884–90.
Kim D, Paggi JM, Park C, Bennett C, Salzberg SL. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol. 2019;37(8):907–15.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R; 1000 Genome Project Data Processing Subgroup. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30.
Jain A, Tuteja G. TissueEnrich: Tissue-specific gene enrichment analysis. Bioinformatics. 2019;35(11):1966–7.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, Xia R. TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol Plant. 2020;13(8):1194–202.
Van der Auwera GA, Carneiro MO, Hartl C, Poplin R, Del Angel G, Levy-Moonshine A, Jordan T, Shakir K, Roazen D, Thibault J, Banks E, Garimella KV, Altshuler D, Gabriel S, DePristo MA. From FastQ data to high confidence variant calls: the Genome Analysis Toolkit best practices pipeline. Curr Protoc Bioinformatics. 2013;43(1110):11.10.1-11.10.33.
Simkin AJ. Genetic Engineering for Global Food Security: Photosynthesis and Biofortification. Plants (Basel). 2019;8(12):586.
Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10(1):57–63.
Ozsolak F, Milos PM. RNA sequencing: advances, challenges and opportunities. Nat Rev Genet. 2011;12(2):87–98.
Oshlack A, Robinson MD, Young MD. From RNA-seq reads to differential expression results. Genome biol. 2010;11(12):1–10.
Schurch NJ, Schofield P, Gierliński M, Cole C, Sherstnev A, Singh V, Wrobel N, Gharbi K, Simpson GG, Owen-Hughes T. How many biological replicates are needed in an RNA-seq experiment and which differential expression tool should you use? RNA. 2016;22(6):839–51.
Kaundal R, Duhan N, Acharya BR, Pudussery MV, Ferreira JFS, Suarez DL, Sandhu D. Transcriptional profiling of two contrasting genotypes uncovers molecular mechanisms underlying salt tolerance in alfalfa. Sci Rep. 2021;11(1):5210.
Lei Y, Xu Y, Hettenhausen C, Lu C, Shen G, Zhang C, Li J, Song J, Lin H, Wu J. Comparative analysis of alfalfa (Medicago sativa L.) leaf transcriptomes reveals genotype-specific salt tolerance mechanisms. BMC Plant Biol. 2018;18(1):35.
Feyissa BA, Renaud J, Nasrollahi V, Kohalmi SE, Hannoufa A. Transcriptome-IPMS analysis reveals a tissue-dependent miR156/SPL13 regulatory mechanism in alfalfa drought tolerance. BMC Genomics. 2020;21(1):721.
Medina CA, Samac DA, Yu LX. Pan-transcriptome identifying master genes and regulation network in response to drought and salt stresses in Alfalfa (Medicago sativa L.). Sci Rep. 2021;11(1):17203.
Zhou Q, Luo D, Chai X, Wu Y, Wang Y, Nan Z, Yang Q, Liu W, Liu Z. Multiple Regulatory Networks Are Activated during Cold Stress in Medicago sativa L. Int J Mol Sci. 2018;19(10):3169.
Ma D, Liu B, Ge L, Weng Y, Weng Y, Cao X, Liu F, Mao P, Mao X. Identification and characterization of regulatory pathways involved in early flowering in the new leaves of alfalfa (Medicago sativa L.) by transcriptome analysis. BMC Plant Biol. 2021;21(1):8.
Cheng Q, Bai S, Ge G, Li P, Liu L, Zhang C, Jia Y. Study on differentially expressed genes related to defoliation traits in two alfalfa varieties based on RNA-Seq. BMC Genomics. 2018;19(1):807.
He F, Long R-C, Zhang T-J, Zhang F, Wang Z, Yang X-J, Jiang X-Q, Yang C-F, Zhi X-X, Li M-N, Yu L-X, Kang J-M, Yang Q-C. Quantitative trait locus mapping of yield and plant height in autotetraploid alfalfa (Medicago sativa L.).The Crop Journal, 2020;8(5):812–818.
Jiang X-Q, Zhang F, Wang Z, Long R-C, Li M.-N, He F, Yang X-J, Yang C-F, Jiang X, Yang Q-C, Wang Q-Z, Kang J-M. Detection of quantitative trait loci (QTL) associated with spring regrowth in alfalfa (Medicago sativa L.). J Integr Agric. 2022;21(3):812–8.
Hawkins C, Yu L-X. Recent progress in alfalfa (Medicago sativa L.) genomics and genomic selection. Crop J. 2018;6(6):565–75.
Derakhshani B, Jafary H, Maleki Zanjani B, Hasanpur K, Mishina K, Tanaka T, Kawahara Y, Oono Y. Combined QTL mapping and RNA-Seq profiling reveals candidate genes associated with cadmium tolerance in barley. PLoS ONE. 2020;15(4):e0230820.
Vanhaeren H, Nam YJ, De Milde L, Chae E, Storme V, Weigel D, Gonzalez N, Inzé D. Forever Young: The Role of Ubiquitin Receptor DA1 and E3 Ligase BIG BROTHER in Controlling Leaf Growth and Development. Plant Physiol. 2017;173(2):1269–82.
Anastasiou E, Kenz S, Gerstung M, MacLean D, Timmer J, Fleck C, Lenhard M. Control of plant organ size by KLUH/CYP78A5-dependent intercellular signaling. Dev Cell. 2007;13(6):843–56.
Li Y, Wei K. Comparative functional genomics analysis of cytochrome P450 gene superfamily in wheat and maize. BMC Plant Biol. 2020;20(1):93.
Tabata R, Sumida K, Yoshii T, Ohyama K, Shinohara H, Matsubayashi Y. Perception of root-derived peptides by shoot LRR-RKs mediates systemic N-demand signaling. Science. 2014;346(6207):343–6.
Domagalska MA, Schomburg FM, Amasino RM, Vierstra RD, Nagy F, Davis SJ. Attenuation of brassinosteroid signaling enhances FLC expression and delays flowering. Development. 2007;134(15):2841–50.
DeYoung BJ, Bickle KL, Schrage KJ, Muskett P, Patel K, Clark SE. The CLAVATA1-related BAM1, BAM2 and BAM3 receptor kinase-like proteins are required for meristem function in Arabidopsis. Plant J. 2006;45(1):1–16.
Gonzalez N, De Bodt S, Sulpice R, Jikumaru Y, Chae E, Dhondt S, Van Daele T, De Milde L, Weigel D, Kamiya Y, Stitt M, Beemster GT, Inzé D. Increased leaf size: different means to an end. Plant Physiol. 2010;153(3):1261–79.
Rafalski A. Applications of single nucleotide polymorphisms in crop genetics. Curr Opin Plant Biol. 2002;5(2):94–100.
We thank Wenshan Guo and Tiejun Zhang for their technical help in the field work and data collection.
This work was supported by the National Natural Science Foundation of China (32071868) and the earmarked fund for China Agriculture Research System (CARS-35–04) and the key research project of Ningxia province for alfalfa breeding program (2019NYY203) and Agricultural Science and Technology Innovation Program (ASTIP-IAS14).
Ethics approval and consent to participate
Study complied with local and national regulations for using plants.
Consent for publication
The authors declare that they have no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Leaf- related QTLs on 32 linkage groups from a genetic linkage map of paternal parent (P1).
Leaf- related QTLs on 32 linkage groups from a genetic linkage map of maternal parent (P2).
Analysis of variance for LL, LW and LA in an F1 population using mixed model for two years in two locations (CP and LF)
List of 2,443 leaf-specific genes.
List of 2,061 up-regulated genes (‘sativa’ Vs ‘falcata’).
List of 1,709 down-regulated genes (‘sativa’ Vs ‘falcata’).
The information of the seven candidates.
About this article
Cite this article
Jiang, X., Yang, X., Zhang, F. et al. Combining QTL mapping and RNA-Seq Unravels candidate genes for Alfalfa (Medicago sativa L.) leaf development. BMC Plant Biol 22, 485 (2022). https://doi.org/10.1186/s12870-022-03864-7