Polymorphisms in candidate genes and their association with carcass traits and meat quality in Nellore cattle

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Single nucleotide polymorphisms in candidate genes associated with gastrointestinal nematode infection in goats

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Seleção fenotípica e identificação de snps em regiões genômicas ligados a resistência à ferrugem em Eucalyptus grandis

Pós-graduação em Agronomia (Proteção de Plantas) - FCA

Association between ACTA1 candidate gene and performance, organs and carcass traits in broilers

This study investigates the genetic association of the SNP present in the ACTA1 gene with performance traits, organs and carcass of broilers to help marker-assisted selection of a paternal broiler line (TT) from EMBRAPA Swine and Poultry, Brazil. Genetic and phenotypic data of 1,400 broilers for 68 traits related to body performance, organ weights, weight of carcass parts, and yields as a percentage of organs and carcass parts were used. The maximum likelihood method, considering 4 analytical models, was used to analyze the genetic association between the SNP and these important economic traits. The association analysis was performed using a mixed animal model including the random effect of the animal (polygenic), and the fixed effects of sex (2 levels), hatch (5 levels) and SNP (3 levels), besides the random error. The traits significantly associated (P < 0.05) with the SNP were analyzed, along with body weight at 42 days of age (BW42), by the restricted maximum likelihood method using the multi-trait animal model to estimate genetic parameters. The analysis included the residual and additive genetic random effects and the sex-hatch fixed effect. The additive effects of the SNP were associated with breast meat (BMY), liver yield (LIVY), body weight at 35 days of age (BW35); drumstick skin (DSW), drumstick (DW) and breast (BW) weights. The heritability estimates for these traits, in addition to BW42, ranged from 0.24 ± 0.06 to 0.45 ± 0.08 for LIVY and BW35, respectively. The genetic correlation ranged from 0.02 ± 0.18 for LIVY and BMY to 0.97 ± 0.01 for BW35 and BW42. Based on the results of this study, it can be concluded that ACTA1 gene is associated with performance traits BW35, LIV and BMY, DW, BW and DW adjusted for body weight at 42 days of age. Therefore, the ACTA1 gene is an important molecular marker that could be used together with others already described to increase the economically important traits in broilers.

QUANTITATIVE TRAIT LOCI FOR AGRONOMIC AND END-USE QUALITY PERFORMANCE AND THE EFFECT OF SOILBORNE WHEAT MOSAIC VIRUS IN A HARD WINTER WHEAT POPULATION IN NEBRASKA

To better understand agronomic and end-use quality in wheat (Triticum aestivum L.) we developed a population containing 154 F6:8 recombinant inbred lines (RILs) from the cross TAM107-R7/Arlin. The parental lines and RILs were phenotyped at six environments in Nebraska and differed for resistance to Wheat soilborne mosaic virus (WSBMV), morphological, agronomic, and end-use quality traits. Additionally, a 2300 cM genome-wide linkage map was created for quantitative trait loci (QTL) analysis. Based on our results across multiple environments, the best RILs could be used for cultivar improvement. The population and marker data are publicly available for interested researchers for future research. The population was used to determine the effect of WSBMV on agronomic and end-use quality and for the mapping of a resistance locus. Results from two infected environments showed that all but two agronomic traits were significantly affected by the disease. Specifically, the disease reduced grain yield by 30% of susceptible RILs and they flowered 5 d later and were 11 cm shorter. End-use quality traits were not negatively affected but flour protein content was increased in susceptible RILs. The resistance locus SbmTmr1 mapped to 27.1 cM near marker wPt-5870 on chromosome 5DL using ELISA data. Finally, we investigated how WSBMV affected QTL detection in the population. QTLs were mapped at two WSBMV infected environments, four uninfected environments, and in the resistant and susceptible RIL subpopulations in the infected environments. Fifty-two significant (LOD≥3) QTLs were mapped in RILs at uninfected environments. Many of the QTLs were pleiotropic or closely linked at 6 chromosomal regions. Forty-seven QTLs were mapped in RILs at WSBMV infected environments. Comparisons between uninfected and infected environments identified 20 common QTLs and 21 environmentally specific QTLs. Finally, 24 QTLs were determined to be affected by WSBMV by comparing the subpopulations in QTL analyses within the same environment. The comparisons were statistically validated using marker by disease interactions. These results showed that QTLs can be affected by WSBMV and careful interpretation of QTL results is needed where biotic stresses are present. Finally, beneficial QTLs not affected by WSBMV or the environment are candidates for marker-assisted selection.

A mixed model QTL analysis for sugarcane multiple-harvest-location trial data

Sugarcane-breeding programs take at least 12 years to develop new commercial cultivars. Molecular markers offer a possibility to study the genetic architecture of quantitative traits in sugarcane, and they may be used in marker-assisted selection to speed up artificial selection. Although the performance of sugarcane progenies in breeding programs are commonly evaluated across a range of locations and harvest years, many of the QTL detection methods ignore two- and three-way interactions between QTL, harvest, and location. In this work, a strategy for QTL detection in multi-harvest-location trial data, based on interval mapping and mixed models, is proposed and applied to map QTL effects on a segregating progeny from a biparental cross of pre-commercial Brazilian cultivars, evaluated at two locations and three consecutive harvest years for cane yield (tonnes per hectare), sugar yield (tonnes per hectare), fiber percent, and sucrose content. In the mixed model, we have included appropriate (co)variance structures for modeling heterogeneity and correlation of genetic effects and non-genetic residual effects. Forty-six QTLs were found: 13 QTLs for cane yield, 14 for sugar yield, 11 for fiber percent, and 8 for sucrose content. In addition, QTL by harvest, QTL by location, and QTL by harvest by location interaction effects were significant for all evaluated traits (30 QTLs showed some interaction, and 16 none). Our results contribute to a better understanding of the genetic architecture of complex traits related to biomass production and sucrose content in sugarcane.

Association of single nucleotide polymorphisms in the bovine leptin and leptin receptor genes with growth and ultrasound carcass traits in Nellore cattle

Given the important role of leptin in metabolism, we looked for a possible association of leptin and leptin receptor polymorphisms with carcass and growth traits in Nellore cattle. We examined associations of leptin and leptin receptor SNPs with ultrasound carcass (longissimus dorsi muscle area (ribeye area), backfat thickness and rump fat thickness and growth traits (weaning weight adjusted to 210 days of age, yearling weight adjusted to 550 days of age, weight gain of weaning to yearling and scrotal circumference adjusted to 550 days of age) of 2162 Bos primigenius indicus (Nellore) animals. Allele and genotypic frequencies were calculated for each marker. Allele substitution, additive and dominance effects of the polymorphisms were also evaluated. Some alleles of the molecular markers had low frequencies, lower than 1%, in the sample analyzed, although the same polymorphisms described for B. p. taurus cattle were found. Due to very low allelic frequencies, the E2JW, A59V and UASMS2 markers were not included in the analysis, because they were almost fixed. E2FB was found to be significantly associated with weight gain, ribeye area and backfat thickness. The promoter region markers, C963T and UASMS1, were also found to be significantly associated with ribeye area. T945M was significantly associated with weight gain. We conclude that the leptin and receptor gene markers would be useful for marker-assisted selection.

Mapping of angular leaf spot resistance QTL in common bean (Phaseolus vulgaris L.) under different environments

Background: Common bean (Phaseolus vulgaris L.) is the most important grain legume for human diet worldwide and the angular leaf spot (ALS) is one of the most devastating diseases of this crop, leading to yield losses as high as 80%. In an attempt to breed resistant cultivars, it is important to first understand the inheritance mode of resistance and to develop tools that could be used in assisted breeding. Therefore, the aim of this study was to identify quantitative trait loci (QTL) controlling resistance to ALS under natural infection conditions in the field and under inoculated conditions in the greenhouse. Results: QTL analyses were made using phenotypic data from 346 recombinant inbreed lines from the IAC-UNA x CAL 143 cross, gathered in three experiments, two of which were conducted in the field in different seasons and one in the greenhouse. Joint composite interval mapping analysis of QTL x environment interaction was performed. In all, seven QTLs were mapped on five linkage groups. Most of them, with the exception of two, were significant in all experiments. Among these, ALS10.1(DG,UC) presented major effects (R-2 between 16% - 22%). This QTL was found linked to the GATS11b marker of linkage group B10, which was consistently amplified across a set of common bean lines and was associated with the resistance. Four new QTLs were identified. Between them the ALS5.2 showed an important effect (9.4%) under inoculated conditions in the greenhouse. ALS4.2 was another major QTL, under natural infection in the field, explaining 10.8% of the variability for resistance reaction. The other QTLs showed minor effects on resistance. Conclusions: The results indicated a quantitative inheritance pattern of ALS resistance in the common bean line CAL 143. QTL x environment interactions were observed. Moreover, the major QTL identified on linkage group B10 could be important for bean breeding, as it was stable in all the environments. Thereby, the GATS11b marker is a potential tool for marker assisted selection for ALS resistance.

QTL detection of yield-related traits of cashew

The identification of quantitative trait loci (QTL) and marker-assisted selection with a view to breeding programs have aroused great interest, including for cashew improvement. This study identified QTL for yield-related traits: nut weight, male and hermaphrodite flowers. The traits were evaluated in 71 F-1 genotypes of the cross CCP 1001 x CP 96. The methods of interval mapping and multiple QTL mapping were applied to identify QTL. Eleven QTL were detected: three for nut weight, four for male flowers and four for hermaphrodite flowers. The QTL accounted for 3.79 to 12.98 % of the total phenotypic variance and had phenotypic effects of -31.81 to 34.25 %. The potential for marker-assisted selection of the QTL hf-2f and hf-3m is great and the phenotypic effects and percentage of phenotypic variation higher than of the others.

Genetic analysis of kernel oil content in tropical maize with design III and QTL mapping

Oil content and grain yield in maize are negatively correlated, and so far the development of high-oil high-yielding hybrids has not been accomplished. Then a fully understand of the inheritance of the kernel oil content is necessary to implement a breeding program to improve both traits simultaneously. Conventional and molecular marker analyses of the design III were carried out from a reference population developed from two tropical inbred lines divergent for kernel oil content. The results showed that additive variance was quite larger than the dominance variance, and the heritability coefficient was very high. Sixteen QTL were mapped, they were not evenly distributed along the chromosomes, and accounted for 30.91% of the genetic variance. The average level of dominance computed from both conventional and QTL analysis was partial dominance. The overall results indicated that the additive effects were more important than the dominance effects, the latter were not unidirectional and then heterosis could not be exploited in crosses. Most of the favorable alleles of the QTL were in the high-oil parental inbred, which could be transferred to other inbreds via marker-assisted backcross selection. Our results coupled with reported information indicated that the development of high-oil hybrids with acceptable yields could be accomplished by using marker-assisted selection involving oil content, grain yield and its components. Finally, to exploit the xenia effect to increase even more the oil content, these hybrids should be used in the Top Cross((TM)) procedure.

MUC1 gene polymorphism in three Nelore lines selected for growth and its association with growth and carcass traits

The objective of this study was to describe the VNTR polymorphism of the mucin 1 gene (MUC1) in three Nelore lines selected for yearling weight to determine whether allele and genotype frequencies of this polymorphism were affected by selection for growth. In addition, the effects of the polymorphism on growth and carcass traits were evaluated. Birth, weaning and yearling weights, rump height, Longissimus muscle area, backfat thickness, and rump fat thickness, were analyzed. A total of 295 Nelore heifers from the Beef Cattle Research Center, Instituto de Zootecnia de Sertozinho, were used, including 41 of the control line, 102 of the selection line and 152 of the traditional. The selection and traditional lines comprise animals selected for higher yearling weight, whereas control line animals are selected for yearling weight close to the average. Five alleles were identified, with allele 1 being the most frequent in the three lines, especially in the lines selected for higher means for yearling weight. Heterozygosity was significantly higher in the control line. Association analyses showed significant effects of allele 1 on birth weight and weaning weight while the allele 3 exert significant effects on yearling weight and back fat thickness. Despite these findings, application of this marker to marker-assisted selection requires more consistent results based on the genotyping of a larger number of animals in order to increase the accuracy of the statistical analyses.

Genetic diversity and population structure of Musa accessions in ex situ conservation

Abstract Background Banana cultivars are mostly derived from hybridization between wild diploid subspecies of Musa acuminata (A genome) and M. balbisiana (B genome), and they exhibit various levels of ploidy and genomic constitution. The Embrapa ex situ Musa collection contains over 220 accessions, of which only a few have been genetically characterized. Knowledge regarding the genetic relationships and diversity between modern cultivars and wild relatives would assist in conservation and breeding strategies. Our objectives were to determine the genomic constitution based on Internal Transcribed Spacer (ITS) regions polymorphism and the ploidy of all accessions by flow cytometry and to investigate the population structure of the collection using Simple Sequence Repeat (SSR) loci as co-dominant markers based on Structure software, not previously performed in Musa. Results From the 221 accessions analyzed by flow cytometry, the correct ploidy was confirmed or established for 212 (95.9%), whereas digestion of the ITS region confirmed the genomic constitution of 209 (94.6%). Neighbor-joining clustering analysis derived from SSR binary data allowed the detection of two major groups, essentially distinguished by the presence or absence of the B genome, while subgroups were formed according to the genomic composition and commercial classification. The co-dominant nature of SSR was explored to analyze the structure of the population based on a Bayesian approach, detecting 21 subpopulations. Most of the subpopulations were in agreement with the clustering analysis. Conclusions The data generated by flow cytometry, ITS and SSR supported the hypothesis about the occurrence of homeologue recombination between A and B genomes, leading to discrepancies in the number of sets or portions from each parental genome. These phenomenons have been largely disregarded in the evolution of banana, as the “single-step domestication” hypothesis had long predominated. These findings will have an impact in future breeding approaches. Structure analysis enabled the efficient detection of ancestry of recently developed tetraploid hybrids by breeding programs, and for some triploids. However, for the main commercial subgroups, Structure appeared to be less efficient to detect the ancestry in diploid groups, possibly due to sampling restrictions. The possibility of inferring the membership among accessions to correct the effects of genetic structure opens possibilities for its use in marker-assisted selection by association mapping.

Functional markers for gene mapping and genetic diversity studies in sugarcane

Abstract Background The database of sugarcane expressed sequence tags (EST) offers a great opportunity for developing molecular markers that are directly associated with important agronomic traits. The development of new EST-SSR markers represents an important tool for genetic analysis. In sugarcane breeding programs, functional markers can be used to accelerate the process and select important agronomic traits, especially in the mapping of quantitative traits loci (QTL) and plant resistant pathogens or qualitative resistance loci (QRL). The aim of this work was to develop new simple sequence repeat (SSR) markers in sugarcane using the sugarcane expressed sequence tag (SUCEST database). Findings A total of 365 EST-SSR molecular markers with trinucleotide motifs were developed and evaluated in a collection of 18 genotypes of sugarcane (15 varieties and 3 species). In total, 287 of the EST-SSRs markers amplified fragments of the expected size and were polymorphic in the analyzed sugarcane varieties. The number of alleles ranged from 2-18, with an average of 6 alleles per locus, while polymorphism information content values ranged from 0.21-0.92, with an average of 0.69. The discrimination power was high for the majority of the EST-SSRs, with an average value of 0.80. Among the markers characterized in this study some have particular interest, those that are related to bacterial defense responses, generation of precursor metabolites and energy and those involved in carbohydrate metabolic process. Conclusions These EST-SSR markers presented in this work can be efficiently used for genetic mapping studies of segregating sugarcane populations. The high Polymorphism Information Content (PIC) and Discriminant Power (DP) presented facilitate the QTL identification and marker-assisted selection due the association with functional regions of the genome became an important tool for the sugarcane breeding program.

Identificazione di un QTL principale per resistenza a ruggine bruna sul cromosoma 7B di frumento duro

Leaf rust caused by Puccinia triticina is a serious disease of durum wheat (Triticum durum) worldwide. However, genetic and molecular mapping studies aimed at characterizing leaf rust resistance genes in durum wheat have been only recently undertaken. The Italian durum wheat cv. Creso shows a high level of resistance to P. triticina that has been considered durable and that appears to be due to a combination of a single dominant gene and one or more additional factors conferring partial resistance. In this study, the genetic basis of leaf rust resistance carried by Creso was investigated using 176 recombinant inbred lines (RILs) from the cross between the cv. Colosseo (C, leaf rust resistance donor) and Lloyd (L, susceptible parent). Colosseo is a cv. directly related to Creso with the leaf rust resistance phenotype inherited from Creso, and was considered as resistance donor because of its better adaptation to local (Emilia Romagna, Italy) cultivation environment. RILs have been artificially inoculated with a mixture of 16 Italian P. triticina isolates that were characterized for virulence to seedlings of 22 common wheat cv. Thatcher isolines each carrying a different leaf rust resistance gene, and for molecular genotypes at 15 simple sequence repeat (SSR) loci, in order to determine their specialization with regard to the host species. The characterization of the leaf rust isolates was conducted at the Cereal Disease Laboratory of the University of Minnesota (St. Paul, USA) (Chapter 2). A genetic linkage map was constructed using segregation data from the population of 176 RILs from the cross CL. A total of 662 loci, including 162 simple sequence repeats (SSRs) and 500 Diversity Arrays Technology markers (DArTs), were analyzed by means of the package EasyMap 0.1. The integrated SSR-DArT linkage map consisted of 554 loci (162 SSR and 392 DArT markers) grouped into 19 linkage blocks with an average marker density of 5.7 cM/marker. The final map spanned a total of 2022 cM, which correspond to a tetraploid genome (AABB) coverage of ca. 77% (Chapter 3). The RIL population was phenotyped for their resistance to leaf rust under artificial inoculation in 2006; the percentage of infected leaf area (LRS, leaf rust susceptibility) was evaluated at three stages through the disease developmental cycle and the area under disease progress curve (AUDPC) was then calculated. The response at the seedling stage (infection type, IT) was also investigated. QTL analysis was carried out by means of the Composite Interval Mapping method based on a selection of markers from the CL map. A major QTL (QLr.ubo-7B.2) for leaf rust resistance controlling both the seedling and the adult plant response, was mapped on the distal region of chromosome arm 7BL (deletion bin 7BL10-0.78-1.00), in a gene-dense region known to carry several genes/QTLs for resistance to rusts and other major cereal fungal diseases in wheat and barley. QLr.ubo-7B.2 was identified within a supporting interval of ca. 5 cM tightly associated with three SSR markers (Xbarc340.2, Xgwm146 e Xgwm344.2), and showed an R2 and an LOD peak value for the AUDPC equal to 72.9% an 44.5, respectively. Three additional minor QTLs were also detected (QLr.ubo-7B.1 on chr. 7BS; QLr.ubo-2A on chr. 2AL and QLr.ubo-3A on chr. 3AS) (Chapter 4). The presence of the major QTL (QLr.ubo-7B.2) was validated by a linkage disequilibrium (LD)-based test using field data from two different plant materials: i) a set of 62 advanced lines from multiple crosses involving Creso and his directly related resistance derivates Colosseo and Plinio, and ii) a panel of 164 elite durum wheat accessions representative of the major durum breeding program of the Mediterranean basin. Lines and accessions were phenotyped for leaf rust resistance under artificial inoculation in two different field trials carried out at Argelato (BO, Italy) in 2006 and 2007; the durum elite accessions were also evaluated in two additional field experiments in Obregon (Messico; 2007 and 2008) and in a green-house experiment (seedling resistance) at the Cereal Disease Laboratory (St. Paul, USA, 2008). The molecular characterization involved 14 SSR markers mapping on the 7BL chromosome region found to harbour the major QTL. Association analysis was then performed with a mixed-linear-model approach. Results confirmed the presence of a major QTL for leaf rust resistance, both at adult plant and at seedling stage, located between markers Xbarc340.2, Xgwm146 and Xgwm344.2, in an interval that coincides with the supporting interval (LOD-2) of QLr.ubo-7B.2 as resulted from the RIL QTL analysis. (Chapter 5). The identification and mapping of the major QTL associated to the durable leaf rust resistance carried by Creso, together with the identification of the associated SSR markers, will enhance the selection efficiency in durum wheat breeding programs (MAS, Marker Assisted Selection) and will accelerate the release of cvs. with durable resistance through marker-assisted pyramiding of the tagged resistance genes/QTLs most effective against wheat fungal pathogens.

A candidate gene approach to identify DNA markers associated with meat quality and production traits: investigation of several cathepsin genes in Italian heavy pigs.

The cathepsin enzymes represent an important family of lysosomal proteinases with a broad spectrum of functions in many, if not in all, tissues and cell types. In addition to their primary role during the normal protein turnover, they possess highly specific proteolytic activities, including antigen processing in the immune response and a direct role in the development of obesity and tumours. In pigs, the involvement of cathepsin enzymes in proteolytic processes have important effects during the conversion of muscle to meat, due to their influence on meat texture and sensory characteristics, mainly in seasoned products. Their contribution is fundamental in flavour development of dry-curing hams. However, several authors have demonstrated that high cathepsin activity, in particular of cathepsin B, is correlated to defects of these products, such as an excessive meat softness together with abnormal free tyrosine content, astringent or metallic aftertastes and formation of a white film on the cut surface. Thus, investigation of their genetic variability could be useful to identify DNA markers associated with these dry cured hams parameters, but also with meat quality, production and carcass traits in Italian heavy pigs. Unfortunately, no association has been found between cathepsin markers and meat quality traits so far, in particular with cathepsin B activity, suggesting that other genes, besides these, affect meat quality parameters. Nevertheless, significant associations were observed with several carcass and production traits in pigs. A recent study has demonstrated that different single nucleotide polymorphisms (SNPs) localized in cathepsin D (CTSD), F (CTSF), H and Z genes were highly associated with growth, fat deposition and production traits in an Italian Large White pig population. The aim of this thesis was to confirm some of these results in other pig populations and identify new cathepsin markers in order to evaluate their effects on cathepsin activity and other production traits. Furthermore, starting from the data obtained in previous studies on CTSD gene, we also analyzed the known polymorphism located in the insulin-like growth factor 2 gene (IGF2 intron3-g.3072G>A). This marker is considered the causative mutation for the quantitative trait loci (QTL) affecting muscle mass and fat deposition in pigs. Since IGF2 maps very close to CTSD on porcine chromosome (SSC) 2, we wanted to clarify if the effects of the CTSD marker were due to linkage disequilibrium with the IGF2 intron3-g.3072G>A mutation or not. In the first chapter, we reported the results from these two SSC2 gene markers. First of all, we evaluated the effects of the IGF2 intron3-g.3072G>A polymorphism in the Italian Large White breed, for which no previous studies have analysed this marker. Highly significant associations were identified with all estimated breeding values for production and carcass traits (P<0.00001), while no effects were observed for meat quality traits. Instead, the IGF2 intron3-g.3072G>A mutation did not show any associations with the analyzed traits in the Italian Duroc pigs, probably due to the low level of variability at this polymorphic site for this breed. In the same Duroc pig population, significant associations were obtained for the CTSD marker for all production and carcass traits (P < 0.001), after excluding possible confounding effects of the IGF2 mutation. The effects of the CTSD g.70G>A polymorphism were also confirmed in a group of Italian Large White pigs homozygous for the IGF2 intron3-g.3072G allele G (IGF2 intron3-g.3072GG) and by haplotype analysis between the markers of the two considered genes. Taken together, all these data indicated that the IGF2 intron3-g.3072G>A mutation is not the only polymorphism affecting fatness and muscle deposition in pigs. In the second chapter, we reported the analysis of two new SNPs identified in cathepsin L (CTSL) and cathepsin S (CTSS) genes and the association results with meat quality parameters (including cathepsin B activity) and several production traits in an Italian Large White pig population. Allele frequencies of these two markers were evaluated in 7 different pig breeds. Furthermore, we mapped using a radiation hybrid panel the CTSS gene on SSC4. Association studies with several production traits, carried out in 268 Italian Large White pigs, indicated positive effects of the CTSL polymorphism on average daily gain, weight of lean cuts and backfat thickness (P<0.05). The results for these latter traits were also confirmed using a selective genotype approach in other Italian Large White pigs (P<0.01). In the 268 pig group, the CTSS polymorphism was associated with feed:gain ratio and average daily gain (P<0.05). Instead, no association was observed between the analysed markers and meat quality parameters. Finally, we wanted to verify if the positive results obtained for the cathepsin L and S markers and for other previous identified SNPs (cathepsin F, cathepsin Z and their inhibitor cystatin B) were confirmed in the Italian Duroc pig breed (third chapter). We analysed them in two groups of Duroc pigs: the first group was made of 218 performance-tested pigs not selected by any phenotypic criteria, the second group was made of 100 Italian Duroc pigs extreme and divergent for visible intermuscular fat trait. In the first group, the CTSL polymorphism was associated with weight of lean cuts (P<0.05), while suggestive associations were obtained for average daily gain and backfat thickness (P<0.10). Allele frequencies of the CTSL gene marker also differed positively among the visible intermuscular extreme tails. Instead, no positive effects were observed for the other DNA markers on the analysed traits. In conclusion, in agreement with the present data and for the biological role of these enzymes, the porcine CTSD and CTSL markers: a) may have a direct effect in the biological mechanisms involved in determining fat and lean meat content in pigs, or b) these markers could be very close to the putative functional mutation(s) present in other genes. These findings have important practical applications, in particular the CTSD and CTSL mutations could be applied in a marker assisted selection (MAS) both in the Italian Large White and Italian Duroc breeds. Marker assisted selection could also increase in efficiency by adding information from the cathepsin S genotype, but only in the Italian Large White breed.