iTRAQ-based Comparative Proteomic Analysis of Differences in the Protein Profiles of Stems and Leaves From Two Alfalfa Genotypes

doi:10.21203/rs.3.rs-36055/v1

Download PDF

Research article

iTRAQ-based Comparative Proteomic Analysis of Differences in the Protein Profiles of Stems and Leaves From Two Alfalfa Genotypes

https://doi.org/10.21203/rs.3.rs-36055/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 29 Sep, 2020

Read the published version in BMC Plant Biology →

You are reading this older preprint version

Read the latest preprint version →

Background:

To explore the molecular regulatory mechanisms of early stem and leaf development, proteomic analysis was performed on leaves and stems of F genotype alfalfa, with thin stems and small leaves, and M genotype alfalfa, with thick stems and large leaves.

Results:

Based on fold-change thresholds of >1.20 or <0.83 (p<0.05), a large number of proteins were identified as being differentially enriched between the M and F genotypes: 249 downregulated and 139 upregulated in stems and 164 downregulated and 134 upregulated in leaves. The differentially expressed proteins in stems were mainly involved in amino acid biosynthesis, phenylpropanoid biosynthesis, carbon fixation, and phenylalanine metabolism. The differentially expressed proteins in leaves were mainly involved in porphyrin and chlorophyll metabolism, phenylpropanoid biosynthesis, starch and sucrose metabolism, and carbon fixation in photosynthetic organisms. Six differentially enriched proteins were mapped onto the porphyrin and chlorophyll metabolism pathway in leaves of the M genotype, including five upregulated proteins involved in chlorophyll biosynthesis and one downregulated protein involved in chlorophyll degradation. Eleven differentially enriched proteins were mapped onto the phenylpropanoid pathway in stems of the M genotype, including two upregulated proteins and nine downregulated proteins.

Conclusion:

Enhanced chlorophyll synthesis and decreased lignin synthesis provided a reasonable explanation for the larger leaves and lower levels of stem lignification in M genotype alfalfa. This proteomic study aimed to classify the functions of differentially enriched proteins and to provide information on the molecular regulatory networks involved in stem and leaf development.

Plant Molecular Biology and Genetics

Plant Physiology and Morphology

Alfalfa

Proteomics

Stem

Leaf

Lignin synthesis

Photosynthesis

In contrast to animals, land plants constantly produce new tissues and organs to complete their life cycle, and the formation of stems and leaves is a crucial aspect of plant growth and development. The plant shoot, which includes leaves, stems, and flowers, differentiates from the shoot apical meristem (SAM) ^[1]. Leaves, which originate from the SAM^[2], are the central organs of photosynthesis and photoperception. In addition, they are the first barriers for plant defense, protecting plants against photochemical damage from UV rays, maintaining plant mechanical properties, and preventing damage from biotic and abiotic factors ^[3]. As for the stem, plant stem cells are dynamic structures composed of polysaccharides, phenolic compounds, and proteins; stems connect plant organs and transport nutrition through the xylem and phloem ^[4]. The plant stem is not only essential for maintaining plant height and rigidity but also responds to environmental stress through the release of signaling molecules ^[5].

Alfalfa (Medicago sativa L.) is a widely cultivated perennial legume forage with high nutritional value ^[6]. Compared with the gradual decline in nutritional value of the alfalfa stem during plant maturation, the nutritional value of alfalfa leaves remains high throughout the growing season ^[7]. During its growth and development, alfalfa’s nutritional value changes as the ratio of stem to leaf increases^[8]. The ratio of stem to leaf is a relevant index that reflects the quality of alfalfa. Tender plants have a high leaf yield, high crude protein content, and high nutritional value. As the plants grow and age, the leaf ratio decreases, dry matter content increases, cellulose content increases, cell walls gradually thicken, and organic sugars increase. The decrease in relative food value of the stem is due to changes in the deposition of cell wall components ^[9]. To elucidate the differences in cell wall concentration and composition of stems of two alfalfa genotypes, several candidate genes were identified in transcript profiling ^[10]. With declining quality, the leaf is an essential source of nutrition for alfalfa and directly affects the quality of alfalfa hay ^[6,11,12].

In recent years, an extensive collection of genes involved in alfalfa stem development and leaf morphogenesis have been identified. However, the molecular regulation of stem development and leaf morphogenesis is an intricate biological process that has yet to be fully revealed. Classical genetic approaches are insufficient to elucidate the molecular regulation mechanisms of leaf morphogenesis and stem development, and proteomics has been used as a complementary approach for protein functional identification ^[13]. The impressive set of proteomics tools and approaches for studying differentially enriched proteins can provide insight into the mechanisms responsible for differences in leaf and stem architecture. Proteome analysis of different stem regions has shown that a subset of proteins involved in cellulose biosynthesis is significantly enriched in intermediate and mature basal regions, resulting in higher cellulose and lignin contents ^[14].

Here we report on a proteomic analysis of leaves and stems from F genotype alfalfa, with thin stems and small leaves (denoted FS and FL), and M genotype alfalfa, with thick stems and large leaves (denoted MS and ML). In contrast to a study of molecular mechanisms of leaf and stem development in the same genetic background published previously, here leaf and stem samples were collected from two alfalfa genotypes with different genetic backgrounds selected from a panel of low and high forage yields. This study was designed to identify differentially enriched proteins related to leaf morphogenesis and stem development, thereby further exploring the regulatory mechanisms of early stem and leaf development.

Phenotypic differences between alfalfa genotypes

There were significant phenotypic differences between the F and M genotypes (Fig. 1a), especially in biomass. Individual plant biomass in M genotype alfalfa was significantly higher than that in F genotype alfalfa, consistent with differences in stem diameter and leaf area (Fig. 1b). The biomass of the M genotype was approximately 50% higher than that of the F genotype, and the leaf area of the M genotype was 2.6 times higher than that of the F genotype. However, although the height of the M genotype alfalfa was slightly greater than that of the F genotype alfalfa, there was no significant difference in height between the two genotypes. Near-infrared reflectance spectroscopy (NIRS) was used to measure the nutritional values of the two genotypes. The M genotype had greater nutritional value than the F genotype, including higher crude protein concentrations and lower cellulose and lignin contents (Supplementary Fig. S1).

Quantitative identification of proteins in alfalfa stems and leaves using iTRAQ

We obtained 14,315 unique peptides, including 3,784 identified proteins (Supplementary Tables 1 and 2). A fold-change threshold of >1.20 or <0.83 (p<0.05) was used to identify differentially enriched proteins (DEPs) between FS and MS and between FL and ML. Using these criteria, 388 DEPs were identified in stems, including 249 downregulated and 139 upregulated DEPs (Fig. 2a, Supplementary Table 3). Likewise, 298 DEPs were identified in leaves: 164 downregulated and 134 upregulated (Fig. 2b, Supplementary Table 4). Compared with F genotype alfalfa, eighty-six DEPs (14.3%) were identified in stems and leaves, including twenty-seven upregulated and forty-six downregulated (Fig. 3a–c). Whereas, the remaining thirteen DEPs showed the opposite trend in stems and leaves.

Functional categorization of differentially enriched proteins

GO enrichment was performed to identify the biological functions of the DEPs, categorizing them according to their biological process, cellular component, and molecular function (Supplementary Tables 5 and 6). The primary biological process categories in stems and leaves were metabolic process, cellular process, response to stimulus, and biosynthetic regulation. The most abundant cellular component categories were cell, cell part, and organelle, and the most abundant molecular function categories in stems (Fig. 4a) and leaves (Fig. 4b) were catalytic activity, binding, structural molecule activity, and transport activity.

KEGG analysis was performed to gain insight into the biochemical pathways of the identified DEPs. The 388 DEPs in stems were mapped to 82 pathways, and the 298 DEPs in leaves were mapped to 76 pathways (Supplementary Tables 7 and 8). Functional enrichment analysis demonstrated significant enrichment in pathways associated with metabolism, biosynthesis of secondary metabolites, carbon metabolism, phenylpropanoid biosynthesis, and carbon fixation in photosynthetic organisms in stems (Fig. 5a). By contrast, pathways associated with porphyrin and chlorophyll II metabolism were significantly enriched in leaves (Fig. 5b). According to the KEGG pathway bubble map, identified proteins (p<0.01) primarily mapped to 28 sub-types in stems and 24 sub-types in leaves (Supplementary Tables 9 and 10).

qRT-PCR analysis of gene expression

To confirm the proteomic results, qRT-PCR was used to examine the mRNA expression of randomly selected proteins using specific primers (Supplementary Table 11). Compared with FS, thirteen DEPs (including auxin-binding protein ABP19a (ABP19a), NAD(P)H:quinone oxidoreductase, phosphoenolpyruvate carboxylase (PPPC), fasciclin domain protein, horseradish peroxidase-like protein (HRP), phenylalanine ammonia-lyase-like protein (PAL), cinnamyl alcohol dehydrogenase (CAD), shikimate/quinate hydroxycinnamoyl-transferase (HCT), riboflavin synthase alpha chain, 4-coumarate:CoA ligase-like protein (4CL-like), and caffeic acid O-methyltransferase (CCOAMT) were upregulated and caffeoyl-CoA 3-O-methyltransferase (CCoAOMT) and HAD-family hydrolase IIA) were identified and downregulated in MS (Fig. 6a). Furthermore, although the mRNA levels of CCOAMT and HAD-family hydrolase IIA were inconsistent with their protein levels, the mRNA expression profiles of eleven genes were consistent with their protein levels in stems (Fig. 6a, Supplementary Tables 3 and 4). In addition, compared with FL, ten DEPs (including light-harvesting complex I chlorophyll A/B-binding protein (LHCI), auxin-repressed/dormancy-associated protein (ARP), white-brown-complex ABC transporter family protein (WBC-ABC), red chlorophyll catabolite reductase (RCCR), auxin-induced in root cultures protein (AIR), cytochrome b5-like heme/steroid-binding domain protein (CYB5), and NAD(P)H dehydrogenase B2 (NDB2) were downregulated in FL. Magnesium-protoporphyrin IX monomethyl ester cyclase (MEPC), bark storage-like protein (BSL), and magnesium chelatase subunit ChlI (CHLI)) were significantly enriched and upregulated in ML. What’s more, the mRNA abundances of nine genes were consistent with their protein levels, while the mRNA abundances of two genes (LHCI and CHLI) were opposite to their protein levels in leaves (Fig. 6b, Supplementary Table 5).

Involvement of lignin biosynthesis in the regulation of stem development and secondary growth

As internodes elongate, stem lignification gradually deepens. Increasing stem maturity accelerates stem lignification and subsequently increases the stem/leaf ratio, resulting in decreased digestibility. Previous studies have shown that the maturity of the stem affects the digestibility of forage, and improving digestibility is one of the primary goals of alfalfa breeding ^[15]. Therefore, in order to monitor changes in the digestibility of alfalfa, recent research has focused on the expression of key genes involved in the phenylpropanoid biosynthesis pathway ^[16–17]. Cell wall synthesis, including the synthesis of cellulose, hemicellulose, and lignin, is a complex biological process. Cellulose biosynthesis is performed by a membrane-bound rosette structure of which sucrose synthase is an integral component ^[18]. The presence of lignin, a complex phenolic polymer formed from three alcohol monomers (coumarin, coniferol, and myrosinol), is the main reason for reduced alfalfa digestibility ^[19]. Lignin biosynthesis requires the participation of a variety of enzymes, and lignin is one of the essential products of the phenylpropanoid metabolic pathway. The specific biosynthesis of monolignols begins with the production of phenylalanine in the shikimate pathway ^[20–21] and continues with the formation of 3-p-hydroxyphenyl, guaiacyl, and syringyl units and finally the synthesis of lignin ^[22].

An interaction network of 21 proteins was involved in phenylpropanoid biosynthesis (Fig. 7), and most of the key enzymes of phenylpropanoid biosynthesis were identified in the current study. In total, 11 DEPs were mapped onto the phenylpropanoid pathway (Fig. 8a), including phenylalanine ammonia-lyase-like protein (PAL), cytochrome P450 family cinnamate 4-hydroxylase (C4H), caffeic acid O-methyltransferase (CCoAMT), 4-coumarate:CoA ligase-like protein (4CL), monoglyceride lipase-like protein (MGL), HXXXD-type acyl-transferase family protein (HCT), caffeoyl-CoA 3-O-methyltransferase (CCoAOMT), cinnamyl alcohol dehydrogenase-like protein (CAD), horseradish peroxidase-like protein (HRP), glycoside hydrolase family 1 protein (GH1), and cytochrome P450 family monooxygenase (F5H). Compared with MS, nine of these proteins were significantly upregulated in FS, whereas HXXXD-type acyl-transferase family protein (HCT) and cytochrome P450 family monooxygenase (F5H) were downregulated. These expression changes provide a reasonable explanation for the higher levels of stem lignification in the F genotype (Supplementary Fig. S1).

Phenylalanine ammonia-lyase-like protein (PAL), an enzyme associated with lignification in primary and secondary tissues, catalyzes the deamination of phenylalanine to initiate phenylpropanoid metabolism ^[23]. The expression of PAL varies coordinately with condensed tannins (CTs) accumulation during the primary to secondary growth transition in stems, and PAL is mainly expressed in non-lignifying cells of stems, leaves, and roots ^[24]. In the current study, PAL was downregulated in MS compared with FS, consistent with the maturation of the stem (Fig. 6a), suggesting that PAL played an essential role in phenylpropanoid metabolism as well as in stem development. A previous study has reported that both coumarate 3-hydroxylase (C3H) and C4H are involved in the early steps of monolignol biosynthesis and exert a negative effect on stem digestibility ^[17]. In the current study, the content of C4H was significantly lower in MS than in FS (Fig. 6a), consistent with lower monolignol biosynthesis and higher stem digestibility in MS.

CAD catalyzes the final step in monolignol synthesis and is therefore a crucial enzyme for the synthesis of S-, G-, and H-lignin ^[25]. 4CL plays a part in the biosynthesis of lignin monomers, particularly guaiacyl (G) lignin ^[26], and there is an overlap in the expression patterns of 4CL family genes ^[27]. In the current study, CAD, 4CL, and C4H were enriched abundantly in the phenylpropanoid biosynthesis pathway. Furthermore, all of these proteins exhibited consistently lower protein and transcript levels in MS than in FS (Fig. 6a), consistent with greater lignin accumulation in FS. Previous studies have reported that the downregulation of CCoAOMT and CCoAMT leads to lower lignin levels and a reduction in G units ^[28]. In the current study, CCoAOMT and CCOAMT were significantly downregulated in MS relative to FS (Fig. 6a), consistent with lignin biosynthesis in MS.

Previous studies have reported that shikimate/quinate HCT plays an essential role in lignin biosynthesis, and there is a strong correlation between HCT accumulation and lignin properties ^[29]. In the current study, protein and transcript levels of shikimate/quinate HCT were lower in MS than in FS, consistent with lower lignin accumulation in MS. Moreover, four lipid-transfer proteins, eight serine proteases, and six fasciclin domain proteins (including 3 fasciclin-like arabinogalactan protein, FLA) were also enriched in stems, and these proteins have been reported to participate in cell wall maturation and secondary growth ^[30–31].

Taking all of the genes mentioned above into consideration, DEPs were involved in lignin biosynthesis and the regulation of stem development and secondary growth. Although the M genotype had a thicker stem (Fig. 1a), the expression of DEPs associated with lignin synthesis was lower in M than in F, suggesting that M accumulated less lignin. This indicates that the M genotype may have higher biological yields and a higher nutritional value than the F genotype (Fig. 1b, Supplementary Fig. S1).

Differences in metabolism-related DEPs between leaves of two alfalfa genotypes

Leaves originate from the SAM and ultimately become flat organs specialized to facilitate light capture and photosynthesis. Leaf morphogenesis, which gives rise to a wide variety of sizes and shapes, is an intricate process that is regulated by many genes from multiple pathways ^[32–35]. Early leaf development is arbitrarily divided into three stages: the initiation of the leaf primordium, the establishment of leaf adaxial-abaxial polarity, and the expansion of the leaf blade ^[36–37]. Previously, KANADI and YABBY transcription factors have been reported to be responsible for the development of abaxial tissues ^[38^–^41]. Initial asymmetric leaf development is regulated by polar YABBY expression ^[42]. Furthermore, gain-of-function alleles of KANADI and YABBY3 (YAB3) result in radial abaxialized organs ^[39–40] and abaxial tissue differentiation, respectively^[38,41].

Leaf development is regulated by internal genetic mechanisms and external environmental cues. Phytohormones, especially auxin, regulate the entire process of leaf development ^[43]. Plant auxin homeostasis is mainly maintained by three processes: de novo IAA biosynthesis, IAA degradation, and IAA conjugation/deconjugation. In addition, IAA carboxyl methyltransferase has been reported to have an essential role in the regulation of auxin homeostasis through the conversion of IAA to methyl-IAA ester (MeIAA) ^[44]. In the current study, indole-3-acetaldehyde oxidase, which is involved in the biosynthesis of IAA, was lower in ML than in FL, suggesting that less auxin was produced in ML ^[45]. In addition, the trend in indole-3-acetaldehyde oxidase levels was consistent with that of ARP and AIR (Fig. 6b), which are involved in auxin homeostasis and play vital roles in leaf development.

Photosynthesis is a process used by plants and other organisms to convert light energy into chemical energy that can later be released to fuel the organisms' activities. Photosystems are functional and structural units of protein complexes involved in photosynthesis that together carry out the primary photochemistry of photosynthesis: the absorption of light and the transfer of energy and electrons. LHCI influences the capture of light energy by photosystem II, which is the key step in the light reactions ^[46]. The chloroplast NAD (P) H dehydrogenase (NDH) complex, as an electron donor, mediates photosystem I (PSI) cyclic and chlororespiratory electron transport ^[47^–^48]. Efficient operation of NDH requires supercomplex formation via minor LHCI in Arabidopsis, and both play essential roles in photosynthesis ^[49]. In the current study, LHCI and NDH were downregulated in ML relative to FL (Fig. 6b). The reduced protein abundance may have resulted in a decreased photosynthetic rate per unit leaf area in ML. However, compared with F genotype alfalfa, larger leaf area in M genotype alfalfa (Fig. 1b) might have compensated for the decrease of photosynthetic rate in ML, which provided a reasonable explanation for higher leaf biomass in ML.

Green plants obtain most of their energy from sunlight via photosynthesis by chloroplasts (including chlorophylls a and b), which gives them their green color. Interestingly, DEPs associated with carbon fixation metabolism and photosynthesis were significantly enriched in leaves, which supply energy for plant development. Chlorophyll, the primary pigment in plant leaves, mainly participates in photosynthesis, and chlorophyll biosynthesis plays an essential role in leaf development. The insertion of magnesium into the chlorophyll molecule is primarily controlled by the activity of magnesium chelatase subunit CHLI, which performs a critical step in the chlorophyll biosynthetic pathway ^[50–51]. In the current study, six DEPs were mapped onto the porphyrin and chlorophyll metabolism pathway, including porphobilinogen deaminase (PBGD), uroporphyrinogen decarboxylase (UPOD), bacteriochlorophyll synthase (BchG), red chlorophyll catabolite reductase (RCCR), CHLI, and magnesium-protoporphyrin IX monomethyl ester cyclase (MPEC), which were significantly enriched in leaves (Fig. 8b). Moreover, both CHLI and MPEC ^[52] were significantly upregulated in ML relative to FL, suggesting that chlorophyll biosynthesis and photosynthetic efficiency were higher in ML (Fig. 6b). In addition, a higher content of granule-bound starch synthase I in ML also suggested the accumulation of more photosynthetic products ^[53], consistent with higher biomass in ML than in FL.

Chlorophyll degradation occurs during leaf senescence, the final stage of leaf development that is regulated by transcription factors and receptor kinases through signal perception and transduction ^[54–55]. Recent research has shown that chlorophyllase, magnesium-chelating substance, and RCCR participate in chlorophyll breakdown ^[56–57]. In the current study, RCCR, a major inducer of cell death ^[58], had consistently lower protein and transcript levels in ML than in FL (Fig. 6b), probably resulting in RCCR accumulation in FL. Enhanced chlorophyll degradation likely further reduced photosynthetic rate, contributing to the differences in biomass observed between the alfalfa genotypes.

Comparative proteomics analysis of DEPs in leaves and stems

Substance synthesis and energy metabolism provide an essential nutrient supply for plants and are indispensable for the completion of normal growth and development. Previous studies have reported that carbohydrate metabolism directly affects plant growth status ^[59]. Eighty-six DEPs (including 27 upregulated and 46 downregulated, and the remaining 13 DEPs showed the opposite trend between stems and leaves) were shared in leaves and stems; (Fig. 3a–c). KEGG enrichment analysis of these proteins revealed significant overlap in enriched KEGG pathways, including energy metabolism, lignin biosynthesis, and phenylalanine metabolism. Likewise, the largest group of DEPs between the alfalfa genotypes was also related to metabolism in leaves and stems.

DEPs involved in glycolysis/gluconeogenesis, the pentose phosphate pathway, and pyruvate metabolism were identified in leaves (Fig. 5b). Similarly, a large number of DEPs involved in the pentose phosphate pathway, the citrate cycle (TCA cycle), and glyoxylate and dicarboxylate metabolism were identified in stems of the alfalfa genotypes (Fig. 5a). Interestingly, carbon metabolism and energy metabolism showed similar trends in stems and leaves. In addition, most of the DEPs identified and upregulated in MS and ML were involved in carbohydrate metabolism and amino acid biosynthesis, indicating that primary metabolism was enhanced to facilitate leaf and stem development and promote increased biomass in the M genotype.

Stem and leaf tissue provide necessary nutrients and material reserves for flower development at the budding stage. In the current study, linoleate 13S-lipoxygenase 2-1 and seed linoleate 9S-lipoxygenase involved in the inositol phospholipid signaling pathway were identified in stems and leaves. Previous studies have shown that various components of the inositol phospholipid signaling system participate in vacuolar changes during pollen development and in vesicle transport during pollen tube growth^[60]. The decreased protein abundance of linoleate 13S-lipoxygenase 2-1 and seed linoleate 9S-lipoxygenase in the M genotype suggested that this genotype could maintain longer vegetative growth and avoid earlier reproductive growth, which was conducive to the accumulation of more biomass.

Anthocyanoin biosynthesis plays an indispensable role in pollen development, especially for alfalfa that is cross-pollinated. To avoid self-pollination in alfalfa, floral pigmentation is conducive to attracting insects and transmitting pollen. Dihydroflavonol 4-reductase (DFR), an essential enzyme for anthocyanin biosynthesis that catalyzes the reaction of dihydroflanovol to leuco- cyanidin (-delphynidin and -pelargonidin) ^[61], was significantly downregulated in M genotype alfalfa, further demonstrating that M genotype alfalfa could maintain longer vegetative growth. Moreover, as the final storage form of photosynthetic products, dark storage protein was upregulated and enriched in the M genotype, providing a further explanation for the increased biomass of the M genotype.

The biological functions of 302 DEPs identified only in stems and 212 DEPs identified only in leaves were further analyzed. In terms of lignin biosynthesis and phenylalanine metabolism, although the DEPs were identified in leaves and stems, almost all the enzymes associated with phenylalanine pathways were significantly enriched in stems, whereas few accumulated in leaves, suggesting that lignin synthesis and phenylalanine metabolism mainly play a role in stem development. However, the differences in biomass between the two genotypes derived from differences in photosynthetic efficiency. Several DEPs involved in porphyrin and chlorophyll metabolism (including PBGD, UPOD, BchG, ChlI, and MPEC) and granule-bound starch synthase I were only identified in the leaf proteome and were upregulated in ML, providing a reasonable explanation for the higher biomass of the M genotype.

Leaf and stem development, responsible for plant morphology and establishment, is a complicated process that runs throughout the whole life history of plants. Especially for alfalfa, there is always a paradox between yield and quality. To obtain high yields and high-quality alfalfa, it is imperative to study the molecular mechanisms that regulate alfalfa stem and leaf development. In the current study, six DEPs were mapped onto the porphyrin and chlorophyll metabolism pathway in ML, including five upregulated proteins involved in chlorophyll biosynthesis and one downregulated protein involved in chlorophyll degradation. At the same time, eleven DEPs were mapped onto the phenylpropanoid pathway in MS, including two upregulated proteins and nine downregulated proteins. Enhanced chlorophyll synthesis and decreased lignin synthesis provided a reasonable explanation for the larger leaves and lower levels of stem lignification in M genotype alfalfa. This analysis of stem and leaf protein expression profiles at the budding stage contributes to our understanding of the molecular mechanisms that underlie alfalfa stem and leaf development.

Plant material. Two alfalfa genotypes, F genotype alfalfa with thin stems and small leaves derived from the local variety “Cangzhou” (CF000735) (F) and M genotype alfalfa with thick stems and large leaves derived from the cultivated variety Zhongmu No.1 (CF032020) (M), were individually selected from a panel of low and high forage yields, respectively. The plant material was collected and cloned by our laboratory, and the approximate dormancy classification of the two alfalfa genotypes was level 3. During the early branching stage in spring 2014, one plant each was chosen from the F and M genotype alfalfa and cloned by cuttage propagation. The cloned plants were then planted in the field of the CAAS research station in Langfang, Hebei Province (39.59°N, 116.59°E). Three individual plants from the ten cloned plants obtained from each genotype were randomly selected for further analysis. Agronomic traits of three randomly selected individual plants were measured at the budding stage of the first stubble under natural conditions in 2017: aboveground biomass, plant height, leaf area, and stem diameter (Fig. 1a). Aboveground biomass was obtained by measuring the dry weight of three individual plants, and height was obtained by measuring the highest branch of three individual plants. The diameters of five randomly selected branches at the shoot base of each plant were measured after harvesting. Leaf areas were obtained by measuring the terminal leaflet of trifoliate leaves located at the third node of three individual plants. Due to differences in the initial flowering period of the two materials, branches were sampled at the same physiological stage (i.e., when the flower buds had not yet opened and were still in bud). Stem and leaf samples were collected simultaneously at the budding stage. Samples of stems (approximately the first three internodes from the stem apex) and leaves (including petioles and stipules and located within the first three internodes) were collected from three individual cloned plants from each genotype and used for proteome and gene expression analysis. The M genotype stem, M genotype leaves, F genotype stem, and F genotype leaves were labeled MS, ML, FS, and FL, respectively. All samples were immediately frozen in liquid nitrogen and stored at −80 °C prior to protein and RNA extraction. Random samples were collected from the two alfalfa genotypes and dried to constant weight at 105 °C for 48 h. The dry samples were crushed into pieces, sieving for 1 mm, and stored under conditions of low temperature, dark and heat insulation. The quality of the dried samples was analyzed with a near infrared spectrometer system (FOSS NIR System 5000, FOSS, Denmark) over the 1100–2500 nm wavelength range at 2 nm wavelength intervals. Measurements were performed under stable working conditions at room temperature (25 °C), and WinI-SI III calibration software was used. Each sample was loaded and scanned three times, the average value was taken, and the spectral data was recorded in the form of log (1/R). Under the same working parameters and conditions, composition analysis (including crude protein, NDF, ADF, and lignin) of the two alfalfa genotypes was performed with same calibration software. All experiments were conducted with three biological replicates. The experimental design and technical process are presented as a technical roadmap in Supplementary Fig. S2.

Protein extraction and purification. The samples were ground to a powder in liquid nitrogen and mixed thoroughly with lysis buffer (7 M urea, 2 M thiourea, 4% CHAPS, 50 ml protease inhibitor cocktail). The mixture was vigorously vortexed and ultrasonicated for 60 s (0.2 s on /2 s off) with an amplitude of 22%. The extract was centrifuged at 15,000 × g for 1 h at 10 °C after extraction for 30 min at room temperature. After centrifugation, the supernatant was collected into a new micro-centrifuge tube that contained 4× volume of precooled 10% trichloroacetic acid (TCA)/acetone, and proteins were precipitated overnight at −20 °C. The samples were centrifuged at 13,000 × g for 10 min at 4 °C, and the supernatant was discarded. The precipitate was washed with acetone and centrifuged. Protein concentration was measured using a Bradford assay kit (Bio-Rad Laboratory), and bovine serum albumin (BSA) was used to generate the standard curve.

Labeling of peptides with iTRAQ reagents. Proteins were digested using the filter-aided sample preparation (FASP) workflow as follows ^[62]. DTT was added to 200 μg of protein sample to a final concentration of 25 mM at 60 °C for 1 h. 50 mM iodoacetamide was added to the mixture at room temperature for 10 min. The reductive alkylated protein solution was centrifuged at 12,000 × g for 20 min, and the filtrate was discarded. 100 μL of dissolution buffer was added, the mixture was centrifuged at 12,000 × g, and the filtrate was finally discarded. The protein samples were digested with trypsin overnight (Promega, Madison, WI, USA) at 37 °C. The peptides of the six samples were labelled with isobaric tags from the iTRAQ Reagent-8plex Multiplex Kit (AB Sciex, USA) according to the manufacturer’s recommended procedure ^[63]. After 1 h of incubation at room temperature, the iTRAQ® reagents were centrifuged to the bottom of the tube. 150 μl of isopropyl alcohol was added to the iTRAQ® reagents, and the mixture was vortexed and centrifuged to the bottom of the tube. 50 μl of sample (100 μg of enzymolytic product) was transferred to a new centrifuge tube. The iTRAQ reagent was added to the sample, the mixture was centrifuged, and the reaction was allowed to proceed at room temperature for 2 h before 100 μl of water was added to stop the reaction. Labeled samples were mixed, vortexed, and centrifuged. Ultimately, equal amounts of protein samples were pooled and concentrated by evaporation using a SpeedVac (Heto, Saskatoon, SK, Canada). After mixing, the labeled samples were dissolved in 100 μl of flowing phase A and centrifuged at 14,000 × g for 20 min.

Liquid chromatography, mass spectrometry and database searching. 100 μl aliquots of labeled samples were loaded at a flow rate of 0.7 ml/min. The components obtained from the reverse phase separation with high pH were redissolved in 20 μl 2% methanol and 0.1% formic acid, followed by centrifugation at 12,000 × g. 10 μl of supernatant was loaded at a flow rate of 0.35 ml/min for 15 min. Separation velocity was set at 300 nl/min.

Database searches were performed as described in a previous publication ^[64]. Tandem mass spectra were extracted by ProteoWizard (version 3.0.8789). All MS/MS spectra were analyzed using Mascot (Matrix Science, London, UK; version 2.6.0). Mascot was set up to search the Medicago truncatula database (version 4.0, 57,693 entries) with a fragment ion mass tolerance of 0.020 Da and a parent ion tolerance of 10.0 PPM, assuming the digestion enzyme trypsin. Carbamidomethyl of cysteine and iTRAQ8plex of lysine and the N-terminus were specified in Mascot as fixed modifications. Oxidation of methionine, acetyl of the n-terminus and iTRAQ8plex of tyrosine were specified in Mascot as variable modifications.

Quantitative data analysis. Label Based Quantitation (iTRAQ, TMT, SILAC, etc.) of identified peptides and proteins was performed with Scaffold Q⁺ (version Scaffold_4.6.2, Proteome Software Inc., Portland, OR) as described in a previous publication ^[63]. Peptide identification was established with a probability of more than 81.0% and an FDR of less than 1.0%. Protein recognition was established with a probability of more than 67.0% and an FDR of less than 10.0%. Protein probabilities were assigned by the protein prophet algorithm ^[65]. Proteins that contained similar peptides and could not be differentiated based on MS/MS analysis alone were grouped to satisfy the principles of parsimony. As described in the previous study, normalization was iterated over intensity ^[66]. Means were used for differential analysis. Spectral data were log-transformed, pruned of those that matched multiple proteins or were missing a reference value, and weighted by an adaptive intensity weighting algorithm.

qRT–PCR analysis. Total RNA was isolated with the SV total RNA extraction system (Promega, Madison, WI, USA) according to the manufacturer’s instructions. First-strand cDNA was synthesized using the Prime Script^TM RT reagent kit (Takara, Beijing, China). qRT–PCR was performed on an ABI Prism 7300 detection system (Applied Biosystems) using the Takara TB Green™ Premix Ex Taq™ II (Tli RNaseH Plus) kit (Takara, Beijing, China). The mRNA transcription levels of eleven DEPs (including red chlorophyll catabolite reductase (XP_003617849.2), white-brown-complex ABC transporter family protein (XP_003630284.1), pectinesterase/pectinesterase inhibitor (XP_003622699.1), NAD(P)H dehydrogenase B2 (XP_013455573.1 (+1)), magnesium-protoporphyrin IX monomethyl ester cyclase (XP_003624856.1), magnesium-chelatase subunit ChlI (XP_003593716.1), light-harvesting complex I chlorophyll A/B-binding protein (XP_013450965.1), auxin-induced in root cultures protein (XP_003626890.1), auxin-repressed/dormancy-associated protein (XP_003593620.1), bark storage-like protein (XP_013462016.1), and cytochrome b5-like heme/steroid-binding domain protein (XP_003627519.1)) identified in leaves and thirteen DEPs (including 4-coumarate:CoA ligase-like protein (XP_003604127.1), auxin-binding protein ABP19a (XP_003627273.1), NAD(P)H: quinone oxidoreductase (XP_003615608.1), horseradish peroxidase-like protein (XP_003602461.1), phosphoenolpyruvate carboxylase (XP_003597124.1), phenylalanine ammonia-lyase-like protein (XP_003591877.1), fasciclin domain protein (XP_003597210.2), cinnamyl alcohol dehydrogenase (XP_003628601.1), HAD-family hydrolase IIA (XP_003608097.1), caffeoyl-CoA 3-O-methyltransferase (XP_003607927.1), riboflavin synthase alpha chain (XP_013443595.1), shikimate/quinate hydroxycinnamoyl-transferase (XP_013454529.1), and caffeic acid O-methyltransferase (XP_003602396.1)) identified in stems were measured, and the β-actin gene (mRNA, EU664318.1) was used as an internal reference. Relative gene expression was quantified using the 2⁻^△^△^CT method ^[67]. The analysis was performed with three biological replicates, and all primer sequences are provided in Supplementary Table 11, and the raw data was listed in Supplementary Table 12-13.

Bioinformatic analysis of proteins. Functional category analysis was performed with BLAST2GO software (http://www.geneontology.org) ^[68]. The KEGG database (http://www.genome.jp/kegg/pathway.html) was used to obtain current information on biochemical pathways and molecular interactions ^[69]. After expression abundances were standardized, hierarchical clustering analysis was performed in R (version 3.2.2). GO and KEGG enrichment analysis were performed to determine the functional subcategories and metabolic pathways of the differentially enriched proteins. A protein–protein interaction network of the DEPs involved in phenylpropanoid biosynthesis was generated using the STRING database (http://string-db.org) ^[70]. Gene expression data were analyzed for statistical significance using a one-way ANOVA (* indicates p<0.05, ** indicates p<0.01, Duncan’s multiple range test) in SPSS 21.0.

iTRAQ: isobaric tags for relative and absolute quantification; NIRS: near-infrared reflectance spectroscopy; NDB2: NAD(P)H dehydrogenase B2; PAL: phenylalanine ammonia-lyase-like protein; C4H: cytochrome P450 family cinnamate 4-hydroxylase; COMT: caffeic acid O-methyltransferase; 4CL: 4-coumarate:CoA ligase-like protein; MGL: monoglyceride lipase-like protein; HCT: HXXXD-type acyl-transferase family protein; CCoAOMT: caffeoyl-CoA 3-O-methyltransferase; CAD: cinnamyl alcohol dehydrogenase-like protein; HRP: horseradish peroxidase-like protein; GH1: glycoside hydrolase family 1 protein; F5H: cytochrome P450 family monooxygenase; PBGD: Porphobilinogen deaminase; UPOD: uroporphyrinogen decarboxylase; BchG: bacteriochlorophyll synthase; RCCR: red chlorophyll catabolite reductase; MPEC: magnesium-protoporphyrin IX monomethyl ester cyclase; BSL: bark storage-like protein; ARP: auxin-repressed/dormancy-associated protein; AIR: auxin-induced in root cultures protein.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Availability of data and material

All data and materials are included within the manuscript.

Competing interests:

The authors declare no competing interests.

Funding

This work was supported by the National Key Basic Research Program of China (973 Program) (2015CB943500), the earmarked fund for China Agriculture Research System (CARS-34), and Agricultural Science and Technology Innovation Program (ASTIP-IAS14). The funding body played no role in the design of the study, the collection, analysis, and interpretation of data and the writing of the manuscript.

Author Contributions

R.L. and Q.Y. conceived and designed the experiments. H.S., J.Y., R.L., and F.Z. performed the experiments. R.L., H.S., J.Y., F.Z., M.L., W.L., J.K., Z.W., and J.Z. analyzed the data. H.S., R.L., J.Y. and Q.Y. wrote the paper. All authors read and approved the final manuscript.

Acknowledgements

This work was supported by the China Forage and Grass Research System (CARS-34), the National Key Basic Research Program of China (973 Program) (2015CB943500), and Agricultural Science and Technology Innovation Program (ASTIP-IAS14). We thank Wenshan Guo and Changfu Yang for their technical help in the field work and data collection.

Rodriguez RE, Debernardi JM, Palatnik JF. Morphogenesis of simple leaves: regulation of leaf size and shape. Wiley Interdisciplinary Reviews Developmental Biology. 2014; 3: 41-57.
Reinhardt D, Kuhlemeier C. Auxin regulates the initiation and radial position of plant lateral organs. Plant Cell. 2000; 12: 507-518.
Petridis A, Therios I, Samouris G, Koundouras S, Giannakoula A. Effect of water deficit on leaf phenolic composition, gas exchange, oxidative damage and antioxidant activity of four Greek olive (Olea europaea L.) cultivars. Plant Physiology & Biochemistry. 2012; 60: 1-11.
Miyashima S, Sebastian J, Lee JY, Helariutta Y. Stem cell function during plant vascular development. EMBO J. 2009; 20: 1097-1106 ().
Baluška F, Volkmann D, Menzel D. Cytoskeleton-plasma membrane-cell wall continuum in plants. Emerging Links Revisited. Plant Physiol. 2003; 133: 482-491.
Vallet C, Lemaire G, Monties B, Chabbert B. Cell wall fractionation of alfalfa stem in relation to internode development: biochemistry aspect. Journal of Agricultural & Food Chemistry. 1998; 46: 3458-3467.
Marković JP. Changes in lignin structure with maturation of alfalfa leaf and stem in relation to ruminants nutrition. African Journal of Agricultural Research. 2012; 7: 257-264.
Nordkvist E, Åman P. Changes during growth in anatomical and chemical composition and in-vitro degradability of Lucerne. Journal of the Science of Food & Agriculture. 1986; 37: 1-7.
Buxton DR, Hornstein JS. Cell-wall concentration and components in stratified canopies of alfalfa, birdsfoot trefoil, and red clover 1. Crop Sci. 1986; 26: 180-184.
Yang SS, Xu WW, Tesfaye M, Lamb JF, Jung HG, VandenBosch KA, Vance CP, Gronwald JW. Transcript profiling of two alfalfa genotypes with contrasting cell wall composition in stems using a cross-species platform: optimizing analysis by masking biased probes. BMC Genomics. 2010; 11: 323.
Engels FM, Jung HG. Alfalfa stem tissues: cell-wall development and lignification. Ann. Bot.-London. 1998; 82: 561-568.
Kumar S. Biotechnological advancements in alfalfa improvement. J. Appl. Genet. 2011; 52: 111-124.
Boschetti E, Righetti PG. Plant proteomics methods to reach low-abundance proteins. Methods in Molecular Biology. 2014; 1072, 111-129.
Printz B, Guerriero G, Sergeant K, Renaut J, Lutts S, Hausman J. Ups and downs in alfalfa: Proteomic and metabolic changes occurring in the growing stem. Plant Sci. 2015; 238: 13-25.
Casler M. Accomplishments and impact from breeding for increased forage nutritional value. Crop Sci. 1999; 39: 12-20.
Jung HJG, Samac DA, Sarath G. Modifying crops to increase cell wall digestibility. Plant Sci. 2012; 65: 185-186.
Reddy MSS, Chen F, Shadle G, Jackson L, Aljoe H, Dixon RA. Targeted down-regulation of cytochrome P450 enzymes for forage quality improvement in alfalfa (Medicago sativa L.). P. Natl. Acad. Sci. USA. 2005, 102(46):16573-16578.
Fujii S, Hayashi T, Mizuno K. Sucrose synthase is an integral component of the cellulose synthesis machinery. Plant & Cell Physiology. 2010; 51: 294.
Boudet AM, Lapierre C, Grimapettenati J. Biochemistry and molecular biology of lignification. New Phytologist. 2010; 129: 203-236.
Vanholme R, Demedts B, Morreel K, Ralph J, Boerjan W. Lignin biosynthesis and structure. Plant Physiol. 2010; 153: 895-905 ().
Vanholme R, Morreel K, Darrah C, Oyarce P, Grabber JH, Ralph J, Boerjan W. Metabolic engineering of novel lignin in biomass crops. New Phytologist. 2012; 196: 978-1000 ().
Wang Y, Maxime C, Richard S, Simon H. Plant cell wall lignification and monolignol metabolism. Front. Plant Sci. 2013; 4: 220.
Kitamura H. Immunocytochemical localization of phenylalanine ammonia-lyase in tissues of Populus kitakamiensis. Planta. 1996; 200: 13-9.
Kao YY Harding SA, Tsai CJ. Differential expression of two distinct phenylalanine ammonia-lyase genes in condensed tannin-accumulating and lignifying cells of quaking aspen. Plant Physiol. 2002; 130: 796-807.
Wyrambik D, Grisebach H. Purification and properties of isoenzymes of cinnamyl-alcohol dehydrogenase from soybean-cell-suspension cultures. FEBS J. 1975; 59: 9-15.
Lee D, Meyer K, Douglas CCAC. Antisense suppression of 4-coumarate: coenzyme a ligase activity in Arabidopsis leads to altered lignin subunit composition. Plant Cell. 1998; 10: 309-309.
Moinuddin SG. Characterization in vitro and in vivo of the putative multigene 4-coumarate: CoA ligase network in Arabidopsis: syringyl lignin and sinapate/sinapyl alcohol derivative formation. Phytochemistry. 2005; 66: 2072-2091 ().
Guo D, Chen F, Inoue K, Blount JW, Dixon RA. Downregulation of caffeic acid 3-O-methyltransferase and caffeoyl CoA 3-O-methyltransferase in transgenic alfalfa: impacts on lignin structure and implications for the biosynthesis of G and S lignin. Plant Cell. (2001); 13: 73-88.
Hoffmann L, Besseau S, Geoffroy P, Ritzenthaler C, Meyer D, Lapierre C, Pollet B, Legrand M. Silencing of hydroxyl-cinnamoyl-coenzyme a shikimate/quinate hydroxyl-cinnamoyl-transferase affects phenylpropanoid biosynthesis. Plant Cell. 2004; 16: 1446-65.
Johnson K, Jones B, Schultz C, Bacic A. “Non-enzymic cell wall (glyco) proteins.” in The Plant Cell Wall ed. Rose JKC, editor. (Boca Raton: CRC Press LLC;) 2003;111–154.
Nieuwland J. Lipid transfer proteins enhance cell wall extension in tobacco. Plant Cell. 2005; 17: 2009-2019.
Uppalapati SR, Ishiga Y, Doraiswamy V, Bedair M, Mysore KS. Loss of abaxial leaf epicuticular wax in Medicago Truncatula irg1/palm1 mutants results in reduced spore differentiation of anthracnose and nonhost rust pathogens. Plant Cell. 2012; 24, 353-370.
Tadege M, Lin H, Bedair M, Berbel A, Wen J, Rojas CM, Niu L, Tang Y, Sumner L, Ratet P. Stenofolia regulates blade outgrowth and leaf vascular patterning in Medicago Truncatula and Nicotiana Sylvestris. Plant cell. 2011; 23, 2125-2142.
Zhou C, Han L, Fu C, Wen J, Cheng X, Nakashima J, Ma J, Tang Y, Tan Y, Tadege M. The trans-acting short interfering RNA3 pathway and no apical meristem antagonistically regulate leaf margin development and lateral organ separation, as revealed by analysis of an argonaute7/lobed leaflet1 mutant in Medicago Truncatula. Plant Cell. 2013; 25, 4845-4862.
Wang H, Chen J, Wen J, Tadege M, Li G, Liu Y, Mysore KS, Chen PRAR. Control of compound leaf development by floricaula/leafy ortholog single leaflet1 in Medicago Truncatula. Plant Physiol. 2008; 146, 1759-1772.
Smith GL. The initiation and determination of leaves. Plant Cell. 1992; 4, 1017-1027.
Tsukaya H. Leaf shape: genetic controls and environmental factors. Int. J. Dev. Biol. 2005; 49, 547-555.
Eshed Y, Baum SF, Bowman AJL. Distinct mechanisms promote polarity establishment in carpels of Arabidopsis. Cell. 1999; 99: 199-209.
Eshed Y, Baum SF, Perea JV, Bowman JL. Establishment of polarity in lateral organs of plants. Curr. Biol. 2001; 11: 1251-1260.
Kerstetter RA, Bollman K, Taylor RA, Bomblies K, Poethig RS. KANADI regulates organ polarity in Arabidopsis. Nature. 2001; 411: 706-709.
Siegfried, K.R. et al. Members of the YABBY gene family specify abaxial cell fate in Arabidopsis. Development. 1999; 126: 4117-4128.
Duruflé H, Clemente HS, Balliau T, Zivy M, Dunand C, Jamet E. et al. Cell wall proteome analysis of Arabidopsis thaliana mature stems. Proteomics. 2017;17: 1600449.
Hay A, Barkoulas M, Tsiantis M. ASYMMETRIC LEAVES1 and auxin activities converge to repress BREVIPEDICELLUS expression and promote leaf development in Arabidopsis. Development. 2006; 133: 3955-61.
Qin G. Gu H, Zhao Y, Ma Z, Shi G, Yang Y, Pichersky E, Chen H, Liu M, Chen Z. An indole-3-acetic acid carboxyl methyltransferase regulates Arabidopsis leaf development. Plant Cell. 2005; 17: 2693-704.
Ljung K, Hull AK, Kowalczyk M, Marchant A, Celenza J, Cohen JD, Sandberg G. Biosynthesis, conjugation, catabolism and homeostasis of indole-3-acetic acid in Arabidopsis thaliana. Plant Mol. Biol. 2002; 49: 249-272.
Leutwiler LS, Meyerowitz EM, Tobin EM. Structure and expression of three light-harvesting chlorophyll a/b-binding protein genes in Arabidopsis thaliana. Nucleic Acids Res. 1986; 14: 4051-64.
Matsuo M, Endo T, Asada K. Properties of the respiratory NAD (P) H dehydrogenase isolated from the Cyanobacterium Synechocystis sp.. Plant & Cell Physiology. 1998; 39: 263-267.
Mi H, Endo T, Ogawa T, Asada K. Thylakoid membrane-bound, NADPH-specific pyridine nucleotide dehydrogenase complex mediates cyclic electron transport in the Cyanobacterium Synechocystis sp.. Plant & Cell Physiology. 1995; 36: 661-668.
Peng L, Fukao Y, Fujiwara M, Takami T, Shikanai T. Efficient operation of NAD (P) H dehydrogenase requires supercomplex formation with photosystem I via minor LHCI in Arabidopsis. Plant Cell. 2009; 21: 3623.
Gibson LC, Jensen PE, Hunter CN. Magnesium chelatase from Rhodobacter sphaeroides: initial characterization of the enzyme using purified subunits and evidence for a BchI-BchD complex. Biochemical Journal. 1999; 337: 243.
Gibson LC, Willows RD, Kannangara CG, von Wettstein D, Hunter CN. Magnesium-protoporphyrin chelatase of Rhodobacter sphaeroides: reconstitution of activity by combining the products of the bchH, -I, and -D genes expressed in Escherichia coli. P. Natl. Acad. Sci. USA. 1995; 92: 1941 ().
Sheng Z, Lv Y, Li W, Luo R, Wei X, Xie L, Hu P. Yellow-Leaf 1 encodes a magnesium-protoporphyrin IX monomethyl ester cyclase, involved in chlorophyll biosynthesis in rice (Oryza sativa L.). PLoS One. 2017; 12: e0177989.
Wattebled F. Granule-bound starch synthase I. FEBS J. 2002; 269: 3810-3820.
Lim PO, Kim HJ, Nam HG. Leaf senescence. Annu. Rev. Plant Biol. 2007; 58: 115-136.
Thomas H, Ougham HJ, Wagstaff C, Stead AD. Defining senescence and death. J. Exp. Bot. 2003; 54: 1127-1132.
Suzuki T, Shioi Y. Re-examination of Mg-dechelation reaction in the degradation of chlorophylls using chlorophyllin a as a substrate. Photosynth. Res. 2002; 74: 217-223.
Takamiya KI, Tsuchiya T, Ohta H. Degradation pathway (s) of chlorophyll: what has gene cloning revealed? Trends Plant Sci. 2000; 5: 426-431.
Pruzinska A, Anders I, Aubry S, Schenk N, Tapernoux-Luthi E, Muller T, Krautler B, Hortensteiner S. In vivo participation of red chlorophyll catabolite reductase in chlorophyll breakdown. Plant Cell. 2007; 19: 369.
Fernando C, Fernie AR. Metabolic regulation underlying tomato fruit development. Journal of Experimental Botany. 2006; 57:1883-1897.
Munnik T, Vermeer JEM (2010) Osmotic stress-induced phosphoinositide and inositol phosphate signalling in plants. Plant Cell Environ 33:655–669
Johnson ET, Ryu S, Yi H, Shin B, Cheong H, Choi G (2010) Alteration of a single amino acid changes the substrate specificity of dihydroflavonol 4-reductase. Plant J 25:325–333
Wisniewski J, Zougman A, Nagaraj N, Mann M. Universal sample preparation method for proteome analysis. Nature Methods. 2009; 6(5):359-362.
Abdallah C, Sergeant K, Guillier C, Dumas-Gaudot E, Leclercq CC, Renaut J. Optimization of iTRAQ labelling coupled to OFFGEL fractionation as a proteomic workflow to the analysis of microsomal proteins of Medicago truncatula roots. Proteome science. 2012; 10(1):37.
Sun H, Long R, Zhang F, Zhang T, Yang Q. Proteomic analysis of shoot tips from two alfalfa cultivars with different florescence. Plant Mol. Biol. Rep. 2019; e46907.
Nesvizhskii AI, Keller A, Kolker E, Aebersold RA. statistical model for identifying proteins by tandem mass spectrometry. Anal. Chem. 2003; 75: 4646-58.
Oberg AL, Mahoney DW, Eckel-Passow JE, Malone CJ, Wolfinger RD, Hill EG, Cooper LT, Onuma OK, Spiro C, Therneau TM. Statistical analysis of relative labeled mass spectrometry data from complex samples using ANOVA. J. Proteome Res. 2008; 7: 225-33 ().
Bustin SA, Benes V, Garson JA. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clinical Chemistry. 2009; 55: 611-622.
Conesa A, Götz S. Blast2GO: A comprehensive suite for functional analysis in plant genomics. International Journal of Plant Genomics, 2008; 2008: 619832.
Kanehisa M. Post-genome informatics. Oxfoed University press. 2000.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva NT, Morris JH, Bork P. String V11: Protein–Protein Association Networks with Increased Coverage, Supporting Functional Discovery in Genome-Wide Experimental Datasets. Nucleic Acids Res. 2019; 47, D607-D613.

Download PDF

Journal Publication

published 29 Sep, 2020

Read the published version in BMC Plant Biology →

Editorial decision: Major revision
27 Jul, 2020
Review #2 received at journal
26 Jul, 2020
Review #1 received at journal
17 Jul, 2020
Reviewer #1 agreed at journal
14 Jul, 2020
Reviewer #2 agreed at journal
14 Jul, 2020
Reviewers invited by journal
13 Jul, 2020
Editor assigned by journal
05 Jul, 2020
Submission checks completed at journal
03 Jul, 2020
Editor invited by journal
03 Jul, 2020
First submitted to journal
18 Jun, 2020

You are reading this older preprint version

Read the latest preprint version →

iTRAQ-based Comparative Proteomic Analysis of Differences in the Protein Profiles of Stems and Leaves From Two Alfalfa Genotypes

Status:

Journal Publication

Version 1

Abstract

Figures

Background

Results

Discussion

Conclusions

Methods

Abbreviations

Declarations

References

Supplementary Files

Status:

Journal Publication

Version 1