QTL mapping for tuberous stem formation of kohlrabi (<Emphasis Type="Italic">Brassica oleracea</Emphasis> var. <Emphasis Type="Italic">gongylodes</Emphasis> L.)

The tuberous stem of kohlrabi is an important quantitative trait, which affects its yield and quality. Genetic control of this trait has not yet been unveiled. To identify the QTLs controlling stem swelling of kohlrabi, a BC₁ population of 92 plants was developed from a cross of broccoli DH line GCP04 and kohlrabi var. Seine. A wide range of variation in tuberous stem diameter was observed among the mapping populations. We constructed a genetic map of nine linkage groups (LGs) with different types of markers, spanning a total length of 913.5 cM with an average marker distance of 7.55 cM. Four significant QTLs for radial enlargement of kohlrabi stem, namely, REnBo1, REnBo2, REnBo3, and REnBo4 were detected on C02, C03, C05, and C09, respectively, and accounted for the phenotypic variation of 59% for the stem diameter and 55% for the qualitative grading of tuberous stem in classes. Then, we confirmed the stability of identified QTLs using BC₁S₁ populations derived from the BC₁ plants having heterozygous alleles at the target QTL and homozygous kohlrabi alleles at the remaining QTLs. REnBo1and REnBo2 using 128 plants of BC₁68S₁ and 94 plants of BC₁43S₁, respectively, and REnBo3 and REnBo4 using 152 plants of BC₁57S₁ were detected at the same positions as the respective QTLs of the BC₁ population. Confirmation of QTLs in two successive generations indicates that the QTLs are persistent. The QTLs obtained in this study could be useful in marker-assisted selection of kohlrabi breeding, and to understand the genetic mechanisms of stem swelling and storage organ development in kohlrabi and other Brassica species.     

Genetic marking of glume color in <Emphasis Type="Italic">Triticum spelta</Emphasis> L. var. <Emphasis Type="Italic">caeruleum</Emphasis> using gliadins

Using gliadins as genetic markers, Triticum spelta L. var. caeruleum accessions were analyzed to identify genetic control of the dark color of glumes. The research material was F₂ and BC₁ plants from crosses between spelt accessions and white-glumed common wheat varieties. The segregation for glume color fitted the monogenic control of the trait. The electrophoretic analysis of gliadins in grains from the hybrid plants has shown that the Gli-Alj* allele in the T. spelta var. caeruleum accessions is linked to the allele for the dark (black) color of glumes at the Rg-A1 locus.     

Genetics and fine mapping of a yellow-green leaf gene (<Emphasis Type="Italic">ygl</Emphasis>-<Emphasis Type="Italic">1</Emphasis>) in cabbage (<Emphasis Type="Italic">Brassica oleracea</Emphasis> var. <Emphasis Type="Italic">capitata</Emphasis> L.)

Cabbage (Brassica oleracea var. capitata L.) is one of the most popular cultivated vegetables worldwide. Cabbage has rich phenotypic diversity, including plant height, head shape, head color, leaf shape and leaf color. Leaf color plays an important role in cabbage growth and development. At present, there are few reports on fine mapping of leaf color mutants in B. oleracea. In this study, a naturally occurring yellow-green leaf cabbage mutant (YL-1), derived from the self-pollinated progenies of the hybrid ‘Hosom’, was used for inheritance analysis and gene mapping. Segregation populations including F₂ and BC₁ were generated from the cross of two inbred lines, YL-1 and 01–20. Genetic analysis with the F₂ and BC₁ populations demonstrated that the yellow-green leaf color was controlled by a single recessive nuclear gene, ygl-1. Insertion–deletion (InDel) markers, designed based on the parental re-sequencing data, were used for the preliminary mapping with BSA (bulked segregant analysis) method. A genetic map constructed with 15 InDels indicated that ygl-1 was located on chromosome C01. The ygl-1 gene is flanked by InDel markers ID2 and M8, with genetic distances of 0.4 cM and 0.35 cM, respectively. The interval distance between two markers is 167 kb. Thus, it enables us to locate the ygl-1 gene for the first time in B. oleracea. This study lays the foundation for candidate gene prediction and ygl-1gene cloning.     

QTL analysis of cotton fiber length in advanced backcross populations derived from a cross between <Emphasis Type="Italic">Gossypium hirsutum</Emphasis> and <Emphasis Type="Italic">G. mustelinum</Emphasis>

Key message

QTLs for fiber length mapped in three generations of advanced backcross populations derived from crossing Gossypium hirsutum and Gossypium mustelinum showed opportunities to improve elite cottons by introgression from wild relatives.

Abstract

The molecular basis of cotton fiber length in crosses between Gossypium hirsutum and Gossypium mustelinum was dissected using 21 BC₃F₂ and 12 corresponding BC₃F_2:3 and BC₃F_2:4 families. Sixty-five quantitative trait loci (QTLs) were detected by one-way analysis of variance. The QTL numbers detected for upper-half mean length (UHM), fiber uniformity index (UI), and short fiber content (SFC) were 19, 20, and 26 respectively. Twenty-three of the 65 QTLs could be detected at least twice near adjacent markers in the same family or near the same markers across different families/generations, and 32 QTLs were detected in both one-way variance analyses and mixed model-based composite interval mapping. G. mustelinum alleles increased UHM and UI and decreased SFC for five, one, and one QTLs, respectively. In addition to the main-effect QTLs, 17 epistatic QTLs were detected which helped to elucidate the genetic basis of cotton fiber length. Significant among-family genotypic effects were detected at 18, 16, and 16 loci for UHM, UI, and SFC, respectively. Six, two, and two loci showed genotype?×?family interaction for UHM, UI and SFC, respectively, illustrating complexities that might be faced in introgression of exotic germplasm into cultivated cotton. Co-location of many QTLs for UHM, UI, and SFC accounted for correlations among these traits, and selection of these QTLs may improve the three traits simultaneously. The simple sequence repeat (SSR) markers associated with G. mustelinum QTLs will assist breeders in transferring and maintaining valuable traits from this exotic source during cultivar development.

     

Multi-environment QTL mapping reveals genetic architecture of fruit cracking in a tomato RIL <Emphasis Type="Italic">Solanum lycopersicum</Emphasis> × <Emphasis Type="Italic">S. pimpinellifolium</Emphasis> population

Key message

QTL and codominant genetic markers for fruit cracking have been identified in a tomato genetic map derived from a RIL population, providing molecular tools for marker-assisted breeding of this trait.

Abstract

In tomato, as well as in other fleshy fruits, one of the main disorders that widely limit quality and production is fruit cracking or splitting of the epidermis that is observed on the fruit skin and flesh at any stage of fruit growth and maturation. To elucidate the genetic basis of fruit cracking, a quantitative trait loci (QTL) analysis was conducted in a recombinant inbred line (RIL) population derived from a cross between tomato (Solanum lycopersicum) and the wild-relative species S. pimpinellifolium. The RIL population was evaluated for fruit cracking during three consecutive growing seasons. Construction of a high-density linkage map based on codominant markers, covering more than 1000 cM of the whole genome, led to the identification of both main and epistatic QTL controlling fruit cracking on the basis of a single-environment as well as multiple-environment analysis. This information will enhance molecular breeding for novel cracking resistant varieties and simultaneously assist the identification of genes underlying these QTL, helping to reveal the genetic basis of fruit cracking in tomato.

     

QTL mapping for haploid male fertility by a segregation distortion method and fine mapping of a key QTL <Emphasis Type="Italic">qhmf4</Emphasis> in maize

Key message

Four QTL related to haploid male fertility were detected by a segregation distortion method and the key QTL qhmf4 was fine mapped to an interval of ~800 kb.

Abstract

Doubled haploid (DH) technology enables rapid development of homozygous lines in maize breeding programs. However, haploid genome doubling is a bottleneck for the commercialization of DH technology and is limited by haploid male fertility (HMF). This is the first study reporting the quantitative trait locus (QTL) analysis of HMF in maize. Four QTL, qhmf1, qhmf2, qhmf3, and qhmf4, controlling HMF have been identified by segregation distortion (SD) loci detection in the selected haploid population derived from ‘Yu87-1/Zheng58’. Three loci, qhmf1, qhmf2, and qhmf4, were also detected in the selected haploid population derived from ‘4F1/Zheng58’. The QTL qhmf4 showed the strongest SD in both haploid populations. Based on the sequence information of ‘Yu87-1’ and ‘Zheng58’, thirteen markers being polymorphic between the two lines were developed to saturate the qhmf4 region. A total of 8168 H₁BC₂ (haploid backcross generation) plants produced from ‘Yu87-1’ and ‘Zheng58’ were screened for recombinants. All the 48 recombinants were backcrossed to ‘Zheng58’ to develop H₁BC₃ progeny. The heterozygous H₁BC₃ individuals were crossed with CAU5 to induce haploids. In each H₁BC₃ progeny, haploids were genotyped and evaluated for anther emergence score (AES). Significant (or no significant) difference (P?<?0.05) between haploids with or without ‘Yu87-1’ donor segment indicated presence or absence of qhmf4 in the donor segment. The analysis of the 48 recombinants narrowed the qhmf4 locus down to an ~800 kb interval flanked by markers IND166 and IND1668.

     

Verification and fine mapping of <Emphasis Type="Italic">qGW1.05</Emphasis>, a major QTL for grain weight in maize (<Emphasis Type="Italic">Zea mays</Emphasis> L.)

Grain weight, one of the important factors to determine corn yield, is a typical quantitative inheritance trait. However, the molecular genetic basis of grain weight still remains limited. In our previous researches, a major QTL associated with grain weight, qGW1.05, has been identified between SSR markers umc1601 and umc1754 at bin locus 1.05–1.06 in maize. Here, its genetic and environmental stabiliteis were verified using a BC₃F₂ population to identify the effect of qGW1.05 on grain weight. Further, qGW1.05-NILs were obtained by MAS successfully. Via a large BC₆F₂ segregation population, together with polymorphic microsatellite markers developed between the parents to screen the genotype of the recombinant plants, qGW1.05 was positioned to a 1.11 Mb genome interval. Furthermore, the progenies of 15 recombinants were tested to confirm the effect of qGW1.05 on grain weight. Combining collinearity among cereal crops and genome annotation, the several candidate genes taking part in grain development were identified in the qGW1.05 region. In this study, qGW1.05 was limited to a 1.11 Mb region on chromosome 1, which established the foundation for understanding the molecular basis underlying kernel development and improving grain weight through MAS using the tightly flanking molecular markers in maize.     

Breeding lines of the Indian mega-rice variety,MTU 1010, possessing protein kinase <Emphasis Type="Italic">OsPSTOL</Emphasis> (<Emphasis Type="Italic">Pup1</Emphasis>), show better root system architecture and higher yield in soils with low phosphorus

MTU 1010 is a high-yielding mega-variety of rice grown extensively in India. However, it does not perform well in soils with low phosphorus (P) levels. With an objective to improve MTU 1010 for tolerance to low soil P, we have transferred Pup1, a major quantitative trait locus (QTL) associated with tolerance from another mega-variety, Swarna, through marker-assisted backcross breeding (MABB). Foreground selection of the F₁ and backcross plants was performed with the co-dominant, closely linked CAPS marker, K20-2, while two flanking markers RM28011 and RM28157 were utilized for recombinant selection. At each backcross generation, positive plants were also analyzed with a set of 85 parental polymorphic SSR markers to identify the QTL-positive plants possessing maximum introgression of MTU 1010 genome. At BC₂F₁, the best backcross plant was selfed to generate BC₂F₂s. Among them, the plants homozygous for Pup1 (n?=?22) were reconfirmed using the functional marker for Pup1, viz., K46-1, and they were advanced through pedigree method of selection until BC₂F₆ generation. A total of five elite BC₂F₆ lines, possessing Pup1 and phenotypically similar to MTU 1010, were screened in the low soil P plot and normal plot (with optimum soil P levels) during wet season, 2016. All the selected lines showed better performance under low P soil with more number of productive tillers, better root system architecture, and significantly higher yield (>?390%) as compared to MTU 1010. Further, under normal soil, the lines were observed to be similar to or better than MTU 1010 for most of the agro-morphological traits and yield. This study represents the successful application of marker-assisted selection for improvement of tolerance to low soil P in a high-yielding Indian rice variety.     

A major QTL (<Emphasis Type="Italic">qFT12.1</Emphasis>) allele from wild soybean delays flowering time

Soybean is highly sensitive to photoperiod. To improve the adaptability and productivity of soybean, it is essential to understand the molecular mechanisms regulating flowering time. To identify new flowering time QTLs, we evaluated a BC₃F₅ population consisting of 120 chromosome segment substitution lines (CSSLs) over 2 years under field conditions. CSSLs were derived from a cross between the cultivated soybean cultivar Jackson and the wild soybean accession JWS156-1, followed by continuous backcrossing using Jackson as the recurrent parent. Four QTLs (qFT07.1, qFT12.1, qFT12.2, and qFT19.1) were detected on three chromosomes. Of these, qFT12.1 showed the highest effect, accounting for 36.37–38.27% of the total phenotypic variation over 2 years. This QTL was further confirmed in the F₇ recombinant inbred line population (n?=?94) derived from the same cross (Jackson × JWS156-1). Analysis of the qFT12.1 BC₃F₅ residual heterozygous line RHL509 validated the allele effect of qFT12.1 and revealed that the recessive allele of qFT12.1 resulted in delayed flowering. Evaluating the qFT12.1 near-isogenic lines (NILs) under different growth conditions showed that NILs with the wild soybean genotype always showed later flowering than those with the cultivated soybean genotype. qFT12.1 was delimited to a 2703-kb interval between the markers BARCSOYSSR_12_0220 and BARCSOYSSR_12_0368 on chromosome 12. qFT12.1 may be a new flowering time gene locus in soybean.     

Fine mapping of QTL <Emphasis Type="Italic">qCTB10</Emphasis>-<Emphasis Type="Italic">2</Emphasis> that confers cold tolerance at the booting stage in rice

Key message

The QTL qCTB10 - 2 controlling cold tolerance at the booting stage in rice was delimited to a 132.5 kb region containing 17 candidate genes and 4 genes were cold-inducible.

Abstract

Low temperature at the booting stage is a major abiotic stress-limiting rice production. Although some QTL for cold tolerance in rice have been reported, fine mapping of those QTL effective at the booting stage is few. Here, the near-isogenic line ZL31-2, selected from a BC₇F₂ population derived from a cross between cold-tolerant variety Kunmingxiaobaigu (KMXBG) and the cold-sensitive variety Towada, was used to map a QTL on chromosome 10 for cold tolerance at the booting stage. Using BC₇F₃ and BC₇F₄ populations, we firstly confirmed qCTB10-2 and gained confidence that it could be fine mapped. QTL qCTB10-2 explained 13.9 and 15.9% of the phenotypic variances in those two generations, respectively. Using homozygous recombinants screened from larger BC₇F₄ and BC₇F₅ populations, qCTB10-2 was delimited to a 132.5 kb region between markers RM25121 and MM0568. 17 putative predicted genes were located in the region and only 5 were predicted to encode expressed proteins. Expression patterns of these five genes demonstrated that, except for constant expression of LOC_Os10g11820, LOC_Os10g11730, LOC_Os10g11770, and LOC_Os10g11810 were highly induced by cold stress in ZL31-2 compared to Towada, while LOC_Os10g11750 showed little difference. Our results provide a basis for identifying the genes underlying qCTB10-2 and indicate that markers linked to the qCTB10-2 locus can be used to improve the cold tolerance of rice at the booting stage by marker-assisted selection.

     

A major locus for resistance to <Emphasis Type="Italic">Botryosphaeria dothidea</Emphasis> in <Emphasis Type="Italic">Prunus</Emphasis>

Species in the fungal family Botryosphaeriaceae are significant pathogens of peach. The climatic conditions in the Southeastern USA are conducive to the development of peach fungal gummosis (PFG) with an estimated yield reduction of up to 40% in severe cases. Genotypes with resistance to this PFG were identified in interspecific crosses and segregating backcross populations generated using Kansu peach (Prunus kansuensis Rehder), almond [Prunus dulcis (Mill.) D.A. Webb], and peach [Prunus persica (L.) Batsch]. Hybrids were evaluated for four consecutive years in field conditions. Data generated was validated in different environments using clonal replicates of the hybrids. The F₁ and BC₁F₁ segregation population data suggest a dominant allele for PFG resistance originating from almond. Segregation and mapping analysis located the PFG resistance locus on a chimeric linkage groups 6–8 near the leaf color locus. The molecular markers identified will facilitate marker-assisted selection (MAS) and introgression of this resistance trait into commercial peach germplasm.     

Inheritance and gene mapping of spotted to non-spotted trait gene <Emphasis Type="Italic">CmSp-1</Emphasis> in melon (<Emphasis Type="Italic">Cucumis melo</Emphasis> L<Emphasis Type="Italic">.</Emphasis> var. <Emphasis Type="Italic">chinensis</Emphasis> Pangalo)

Melon (Cucumis melo L.) is one of the most popular and highly nutritious vegetable species within Cucurbitaceae. Because appearance is used as an important indicator of quality, the spotted to non-spotted trait associated with this product somewhat influences the buying habits of consumers. We tested a six-generation family to determine the inheritance and genetic basis of this trait. Genetic groups F₁, F₂, BC₁P₁, and BC₁P₂ were from a cross between “IM16559” (non-spotted) and “IM16553” (spotted). Our genetic analysis showed that the spotted to non-spotted trait was controlled by a single dominant gene that we named CmSp-1. Whole-genome resequencing-bulked segregant analysis (WG-BSA) demonstrated that this gene was located on the end of chromosome 2, in the intersections of 22,160,000 to 22,180,000 bp and 22,260,000 to 26,180,000 bp, an interval distance of 3.94 Mb. Insertion-deletion (InDel) markers designed based on WG-BSA data were used to map this gene. Using 13 InDel markers, we produced a genetic map indicating that CmSp-1 was tightly linked to markers I734-2 and I757, with genetic distances of 1.8 and 0.4 cM and an interval distance of 280.872 kb. The closest marker was I757. Testing of 107 different melon genotypes presented an accuracy of 84.11% in predicting the phenotype. By being able to locate CmSp-1 in melon, we can now use the findings to identify potential targets for further marker-assisted breeding and cloning projects.     

Improvement of blast resistance of the popular high-yielding,medium slender-grain type,bacterial blight resistant rice variety,Improved Samba Mahsuri by marker-assisted breeding

Improved Samba Mahsuri (ISM) is a popular, high-yielding, bacterial blight resistant rice variety possessing medium-slender grain type. As ISM is highly susceptible to blast disease of rice, through the present study we have transferred two major blast resistance genes, Pi2 and Pi54 into the elite variety by marker-assisted backcross breeding. The two blast resistance genes were transferred to ISM through sets of backcrosses. In every backcross generation, PCR-based markers, specific for the blast resistance genes (Pi2 and Pi54) and bacterial blight resistance genes (Xa21, xa13 and xa5) were utilized for foreground selection, while a set of 144 parental polymorphic SSR markers were used for background selection and backcrossing was carried out until BC₂ generation. A solitary BC₂F₁ plant possessing Pi2 or Pi54 along with Xa21, xa13 and xa5 and >?90% recovery of ISM genome was selected from the two sets of backcrosses were crossed and the intercross F₁s (ICF₁s) thus obtained were selfed to generate ICF₂s. Homozygous ICF₂ plants carrying all the five resistance genes were identified through markers and advanced through selfing till ICF₅ generation by adopting pedigree method of selection. Three best lines at ICF₅, possessing excellent resistance against bacterial blight and blast and closely resembling or superior to ISM in terms of grain quality: yield and agro-morphological traits have been identified and advanced for multi-location trials.     

Discovery and introgression of the wild sunflower-derived novel downy mildew resistance gene <Emphasis Type="Italic">Pl</Emphasis><Subscript><Emphasis Type="Italic">19</Emphasis></Subscript> in confection sunflower (<Emphasis Type="Italic">Helianthus annuus</Emphasis> L.)

Key message

A new downy mildew resistance gene, Pl ₁₉ , was identified from wild Helianthus annuus accession PI 435414, introduced to confection sunflower, and genetically mapped to linkage group 4 of the sunflower genome.

Abstract

Wild Helianthus annuus accession PI 435414 exhibited resistance to downy mildew, which is one of the most destructive diseases to sunflower production globally. Evaluation of the 140 BC₁F_2:3 families derived from the cross of CMS CONFSCLB1 and PI 435414 against Plasmopara halstedii race 734 revealed that a single dominant gene controls downy mildew resistance in the population. Bulked segregant analysis conducted in the BC₁F₂ population with 860 simple sequence repeat (SSR) markers indicated that the resistance derived from wild H. annuus was associated with SSR markers located on linkage group (LG) 4 of the sunflower genome. To map and tag this resistance locus, designated Pl ₁₉ , 140 BC₁F₂ individuals were used to construct a linkage map of the gene region. Two SSR markers, ORS963 and HT298, were linked to Pl ₁₉ within a distance of 4.7 cM. After screening 27 additional single nucleotide polymorphism (SNP) markers previously mapped to this region, two flanking SNP markers, NSA_003564 and NSA_006089, were identified as surrounding the Pl ₁₉ gene at a distance of 0.6 cM from each side. Genetic analysis indicated that Pl ₁₉ is different from Pl ₁₇ , which had previously been mapped to LG4, but is closely linked to Pl ₁₇ . This new gene is highly effective against the most predominant and virulent races of P. halstedii currently identified in North America and is the first downy mildew resistance gene that has been transferred to confection sunflower. The selected resistant germplasm derived from homozygous BC₂F₃ progeny provides a novel gene for use in confection sunflower breeding programs.

     

Cloning,identification, and expression analysis of a Dicer-Like gene family from Solanum lycopersicum

Dicer proteins belong to the RNase III family of proteins, which are key components in small RNA biogenesis. In Solanum lycopersicum, seven Dicer-like (DCL) genes have been identified and have been named SlDCL. In this study, we cloned the full-length sequence of the SlDCL genes including untranslated regions using RNA ligase-mediated rapid amplification of cDNA ends. Our analysis indicates that 7 SlDCLs were located on 5 tomato chromosomes (6, 7, 8, 10, and 11). The gene structure of the SlDCLs covered long genomic regions and contained more than 20 exons. Phylogenetic analysis divided the seven SlDCL members into four subgroups. In general, all seven SlDCLs were expressed in all organs but more in flowers and fruits than in the other parts. Moreover, the expressions of some genes changed slightly after treatment with ethylene or 1-methylcyclopropene suggesting their likely roles in plant responses to ethylene. Our findings provide essential information on SlDCL genes in tomato and will aid in the functional classification of DCL families in plants.     

Alterations in the plasma membrane polypeptide pattern of tomato roots (Lycopersicon esculentum) during the development of arbuscular mycorrhiza

Changes induced by arbuscular mycorrhizal (AM) formation in the plasma membrane polypeptide pattern of tomato roots have been assessed by 2D-PAGE analysis. Plasma membrane fractions were isolated by aqueous two-phase partitioning from control and mycorrhizal tomato root microsomes. Analysis of 2D-PAGE gels revealed that AM colonization induces at the plasma membrane level two major changes in protein synthesis: down-regulation of some constitutive polypeptides and synthesis of new polypeptides or endomycorrhizins. A comparison of changes induced by two different levels of AM colonization showed that 16 polypeptides were differentially displayed at both AM colonization stages, while some others were transiently regulated. Five of the differentially displayed plasma membrane polypeptides at both AM colonization stages were selected for N-terminal amino acid sequencing. Reliable sequences were obtained for two of the selected spots. Sequence alignment search indicated that one of the sequenced polypeptides showed 75% identity to the N-terminal sequence of the 69 kDa catalytic subunit of the vacuolar type H(+)-ATPase of several plants. The possible significance of these findings is discussed in relation to the functioning of the AM symbiosis.     

Comparative proteomic analysis of leaf tissue from tomato plants colonized with <Emphasis Type="Italic">Rhizophagus irregularis</Emphasis>

A comparative proteomic approach was performed to analyze the differential accumulation of leaf proteins in response to the symbiosis between Solanum lycopersicum and the arbuscular mycorrhizal fungus (AMF) Rhizophagus irregularis. Protein profiling was examined in leaves from tomato plants colonized with AMF (M), as well as non-colonized plants fertilized with low phosphate (20 μM P; NM-LP) and non-colonized plants fertilized with regular phosphate Hoagland’s solution (200 μM P; NM-RP). Comparisons were made between these groups, and 2D-SDS-PAGE revealed that 27 spots were differentially accumulated in M vs. NM-LP. Twenty-three out of the 27 spots were successfully identified by mass spectrometry. Two of these proteins, 2-methylene-furan-3-one reductase and auxin-binding protein ABP19a, were up-accumulated in M plants. The down-accumulated proteins in M plants were associated mainly with photosynthesis, redox, and other molecular functions. Superoxide dismutase, harpin binding protein, and thioredoxin peroxidase were down-accumulated in leaves of M tomato plants when compared to NM-LP and NM-RP, indicating that these proteins are responsive to AMF colonization independently of the phosphate regime under which they were grown. 14-3-3 protein was up-accumulated in NM-RP vs. NM-LP plants, whereas it was down-accumulated in M vs. NM-LP and M vs. NM-RP, regardless of their phosphate nutrition. This suggests a possible regulation by P nutrition and AMF colonization. Our results demonstrate AMF-induced systemic changes in the expression of tomato leaf proteins, including the down-accumulation of proteins related to photosynthesis and redox function.