Identification of genomic regions, genetic variants and gene networks regulating candidate genes for lipid metabolism in pig muscle.

The intramuscular fat content and fatty acid composition of porcine meat have a significant impact on its quality and nutritional value. This research aimed to investigate the expression of 45 genes involved in lipid metabolism in the longissimus dorsi muscle of three experimental pig backcrosses, with a 25% of Iberian background. To achieve this objective, we conducted an expression Genome-Wide Association Study (eGWAS) using gene expression levels in muscle measured by high-throughput real-time qPCR for 45 target genes and genotypes from the PorcineSNP60 BeadChip or Axiom Porcine Genotyping Array and 65 single nucleotide polymorphisms (SNPs) located in 20 genes genotyped by a custom-designed Taqman OpenArray in a cohort of 354 animals. The eGWAS analysis identified 301 eSNPs associated with 18 candidate genes (ANK2, APOE, ARNT, CIITA, CPT1A, EGF, ELOVL6, ELOVL7, FADS3, FASN, GPAT3, NR1D2, NR1H2, PLIN1, PPAP2A, RORA, RXRA and UCP3). Three cis-eQTL (expression quantitative trait loci) were identified for GPAT3, RXRA, and UCP3 genes, which indicates that a genetic polymorphism proximal to the same gene is affecting its expression. Furthermore, 24 trans-eQTLs were detected, and eight candidate regulatory genes were located in these genomic regions. Additionally, two trans-regulatory hotspots in Sus scrofa chromosomes 13 and 15 were identified. Moreover, a co-expression analysis performed on 89 candidate genes and the fatty acid composition revealed the regulatory role of four genes (FABP5, PPARG, SCD, and SREBF1). These genes modulate the levels of α-linolenic, arachidonic, and oleic acids, as well as regulating the expression of other candidate genes associated with lipid metabolism. The findings of this study offer novel insights into the functional regulatory mechanism of genes involved in lipid metabolism, thereby enhancing our understanding of this complex biological process.

Understanding host-microbiota interactions in the commercial piglet around weaning.

Weaning is a critical period in the life of pigs with repercussions on their health and welfare and on the economy of the swine industry. This study aimed to assess the effect of the commercial early weaning on gut microbiota, intestinal gene expression and serum metabolomic response via an integrated-omic approach combining 16S rRNA gene sequencing, the OpenArray gene expression technology and 1H-NMR spectroscopy. Fourteen piglets from different litters were sampled for blood, jejunum tissue and caecal content two days before (- 2d), and three days after (+ 3d) weaning. A clearly differential ordination of caecal microbiota was observed. Higher abundances of Roseburia, Ruminococcus, Coprococcus, Dorea and Lachnospira genera in weaned piglets compared to prior to weaning showed the quick microbial changes of the piglets' gut microbiota. Downregulation of OCLN, CLDN4, MUC2, MUC13, SLC15A1 and SLC13A1 genes, also evidenced the negative impact of weaning on gut barrier and digestive functions. Metabolomic approach pinpointed significant decreases in choline, LDL, triglycerides, fatty acids, alanine and isoleucine and increases in 3-hydroxybutyrate after weaning. Moreover, the correlation between microbiota and metabolome datasets revealed the existence of metabolic clusters interrelated to different bacterial clusters. Our results demonstrate the impact of weaning stress on the piglet and give insights regarding the associations between gut microbiota and the animal gene activity and metabolic response.

Early socialization and environmental enrichment of lactating piglets affects the caecal microbiota and metabolomic response after weaning.

The aim of this study was to determine the possible impact of early socialization and an enriched neonatal environment to improve adaptation of piglets to weaning. We hypothesized that changes in the microbiota colonization process and in their metabolic response and intestinal functionality could help the animals face weaning stress. A total of 48 sows and their litters were allotted into a control (CTR) or an enriched treatment (ENR), in which piglets from two adjacent pens were combined and enriched with toys. The pattern of caecal microbial colonization, the jejunal gene expression, the serum metabolome and the intestinal physiology of the piglets were assessed before (-2 d) and after weaning (+ 3d). A differential ordination of caecal microbiota was observed after weaning. Serum metabolome suggested a reduced energetic metabolism in ENR animals, as evidenced by shifts in triglycerides and fatty acids, VLDL/LDL and creatine regions. The TLR2 gene showed to be downregulated in the jejunum of ENR pigs after weaning. The integration of gene expression, metabolome and microbiota datasets confirmed that differences between barren and enriched neonatal environments were evident only after weaning. Our results suggest that improvements in adaptation to weaning could be mediated by a better response to the post-weaning stress.

Transversal gene expression panel to evaluate intestinal health in broiler chickens in different challenging conditions.

There is a high interest on gut health in poultry with special focus on consequences of the intestinal diseases, such as coccidiosis and C. perfringens-induced necrotic enteritis (NE). We developed a custom gene expression panel, which could provide a snapshot of gene expression variation under challenging conditions. Ileum gene expression studies were performed through high throughput reverse transcription quantitative real-time polymerase chain reaction. A deep review on the bibliography was done and genes related to intestinal health were selected for barrier function, immune response, oxidation, digestive hormones, nutrient transport, and metabolism. The panel was firstly tested by using a nutritional/Clostridium perfringens model of intestinal barrier failure (induced using commercial reused litter and wheat-based diets without exogenous supplementation of enzymes) and the consistency of results was evaluated by another experiment under a coccidiosis challenge (orally gavaged with a commercial coccidiosis vaccine, 90× vaccine dose). Growth traits and intestinal morphological analysis were performed to check the gut barrier failure occurrence. Results of ileum gene expression showed a higher expression in genes involved in barrier function and nutrient transport in chickens raised in healthy conditions, while genes involved in immune response presented higher expression in C.perfringens-challenged birds. On the other hand, the Eimeria challenge also altered the expression of genes related to barrier function and metabolism, and increased the expression of genes related to immune response and oxidative stress. The panel developed in the current study gives us an overview of genes and pathways involved in broiler response to pathogen challenge. It also allows us to deep into the study of differences in gene expression pattern and magnitude of responses under either a coccidial vaccine or a NE.

Association of genetic variants and expression levels of porcine FABP4 and FABP5 genes.

The FABP4 and FABP5 genes, coding for fatty acid transport proteins, have long been studied as positional candidate genes for SSC4 QTL affecting fat deposition and composition traits in pigs. Polymorphisms in these genes, FABP4:g.2634_2635insC and FABP5:g.3000T>G, have previously been associated with fatness traits in an Iberian by Landrace cross (IBMAP). The aim of the present work was to evaluate the functional implication of these genetic variants. For this purpose, FABP4 and FABP5 mRNA expression levels in 114 BC1_LD animals (25% Iberian × 75% Landrace) were analyzed using real-time quantitative PCR in backfat and muscle. FABP4 gene expression in backfat, but not in muscle, was associated with FABP4:g.2634_2635insC. In contrast, FABP5:g.3000T>G was not associated with gene expression levels. An expression-based genome-wide association study highlighted the FABP4:g.2634_2635insC polymorphism as the polymorphism most associated with FABP4 gene expression in backfat. Furthermore, other genomic regions associated in trans with the mRNA expression of FABP4 in backfat and FABP5 in muscle were also identified. Finally, two putative transcription binding sites for PPARG and NR4A2 may be affected by the FABP4:g.2634_2635insC polymorphism, modifying FABP4 gene expression. Our results reinforce FABP4 as a candidate gene for fatness traits on SSC4.

Using RNA-Seq SNP data to reveal potential causal mutations related to pig production traits and RNA editing.

RNA-Seq technology is widely used in quantitative gene expression studies and identification of non-annotated transcripts. However this technology also can be used for polymorphism detection and RNA editing in transcribed regions in an efficient and cost-effective way. This study used SNP data from an RNA-Seq assay to identify genes and mutations underlying production trait variations in an experimental pig population. The hypothalamic and hepatic transcriptomes of nine extreme animals for growth and fatness from an (Iberian × Landrace) × Landrace backcross were analyzed by RNA-Seq methodology, and SNP calling was conducted. More than 125 000 single nucleotide variants (SNVs) were identified in each tissue, and 78% were considered to be potential SNPs, those SNVs segregating in the context of this study. Potential informative SNPs were detected by considering those showing a homozygous or heterozygous genotype in one extreme group and the alternative genotype in the other group. In this way, 4396 and 1862 informative SNPs were detected in hypothalamus and liver respectively. Out of the 32 SNPs selected for validation, 25 (80%) were confirmed as actual SNPs. Association analyses for growth, fatness and premium cut yields with 19 selected SNPs were carried out, and four potential causal genes (RETSAT, COPA, RNMT and PALMD) were identified. Interestingly, new RNA editing modifications were detected and validated for the NR3C1:g.102797 (ss1985401074) and ACSM2B:g.13374 (ss1985401075) positions and for the COG3:g3.4525 (ss1985401087) modification previously identified across vertebrates, which could lead to phenotypic variation and should be further investigated.

Analysis of the porcine APOA2 gene expression in liver, polymorphism identification and association with fatty acid composition traits.

APOA2 is a protein implicated in triglyceride, fatty acid and glucose metabolism. In pigs, the APOA2 gene is located on pig chromosome 4 (SSC4) in a QTL region affecting fatty acid composition, fatness and growth traits. In this study, we evaluated APOA2 as a candidate gene for meat quality traits in an Iberian × Landrace backcross population. The APOA2:c.131T>A polymorphism, located in exon 3 of APOA2 and determining a missense mutation, was associated with the percentage of hexadecenoic acid [C16:1(n-9)], linoleic acid [C18:2(n-6)], α-linolenic acid [C18:3(n-3)], dihomo-gamma-linolenic acid [C20:3(n-6)] and polyunsaturated fatty acids (PUFAs) in backfat. Furthermore, this SNP was associated with the global mRNA expression levels of APOA2 in liver and was used as a marker to determine allelic expression imbalance by pyrosequencing. We determined an overexpression of the T allele in heterozygous samples with a mean ratio of 2.8 (T/A), observing a high variability in the allelic expression among individuals. This result suggests that complex regulatory mechanisms, beyond a single polymorphism (e.g. epigenetic effects or multiple cis-acting polymorphisms), may be regulating APOA2 gene expression.

Transcriptional analysis of intramuscular fatty acid composition in the longissimus thoracis muscle of Iberian × Landrace back-crossed pigs.

This study aimed at identifying differential gene expression conditional on the fatty acid profile of the longissimus thoracis (Lt) muscle, a prime cut of economic relevance for fresh and cured pork production. A population of 110 Iberian (25%) × Landrace (75%) back-crossed pigs was used, because these two breeds exhibit extreme profiles of intramuscular saturated fatty acid, monounsaturated fatty acid (MUFA) and polyunsaturated fatty acid (PUFA) contents. Total RNA from Lt muscle was individually hybridized to GeneChip Porcine Genome arrays (Affymetrix). A principal component analysis was performed with data from the 110 animals to select 40 extreme animals based on the total fatty acid profile and the MUFA composition (MAP). Comparison of global transcription levels between extreme fatty acid profile pigs (n = 40) resulted in 219 differentially expressed probes (false discovery rate <0.10). Gene ontology, pathway and network analysis indicated that animals with higher percentages of PUFA exhibit a shift toward a more oxidative muscular metabolism state, with a raise in mitochondria function (PPARGC1A, ATF2), fatty acid uptake and oxidation (FABP5, MGLL). On the other hand, 87 probes were differentially expressed between MUFA composition groups (n = 40; false discovery rate <0.10). In particular, muscles rich in n-7 MUFA expressed higher levels of genes involved in lipid metabolism (GLUL, CRAT, PLA2G15) and lower levels of fatty acid elongation genes (ELOVL5). Moreover, the chromosomal position of FABP5, PAQR3, MGLL, PPARGC1A, GLUL and ELOVL5 co-localized with very relevant QTL for fat deposition and composition described in the same resource population. This study represents a complementary approach to identifying genes underlying these QTL effects.

Identification of differentially expressed genes in the oviduct of two rabbit lines divergently selected for uterine capacity using suppression subtractive hybridization.

Suppressive subtractive hybridization libraries from oviduct at 62 h post-mating of two lines of rabbits divergently selected for uterine capacity were generated to identify differentially expressed genes. A total of 438 singletons and 126 contigs were obtained by cluster assembly and sequence alignment of 704 expressed sequence tags (ESTs), of which 54% showed homology to known proteins of the non-redundant NCBI databases. Differential screening by dot blot validated 71 ESTs, of which 47 showed similarity to known genes. Transcripts of genes were functionally annotated in the molecular function and the biological process gene ontology categories using the BLAST2GO software and were assigned to reproductive developmental process, immune response, amino acid metabolism and degradation, response to stress and apoptosis terms. Finally, three interesting genes, PGR, HSD17B4 and ERO1L, were identified as overexpressed in the low line using RT-qPCR. Our study provides a list of candidate genes that can be useful to understanding the molecular mechanisms underlying the phenotypic differences observed in early embryo survival and development traits.

Genome-wide association study for intramuscular fatty acid composition in an Iberian × Landrace cross.

The lipid content and fatty acid (FA) profile have an important impact in human health as well as in the technological transformation and nutritional and organoleptic quality of meat. A genome-wide association study (GWAS) on 144 backcross pigs (25% Iberian × 75% Landrace) was performed for 32 traits associated with intramuscular FA composition and indices of FA metabolism. The GWAS was carried out using Qxpak 5.0 and the genotyping information obtained from the Porcine SNP60K BeadChip (Illumina Inc., San Diego, CA). Signals of significant association considering a false- discovery rate (q-value < 0.05) were observed in 15 of the 32 analyzed traits, and a total of 813 trait-associated SNP (TAS), distributed in 43 chromosomal intervals on almost all autosomes, were annotated. According to the clustering analysis based on functional classification, several of the annotated genes are related to FA composition and lipid metabolism. Some interesting positional concordances among TAS and previously reported QTL for FA compositions and/or other lipid traits were also found. These common genomic regions for different traits suggest pleiotropic effects for FA composition and were found primarily on SSC4, SSC8, and SSC16. These results contribute to our understanding of the complex genetic basis of FA composition and FA metabolism.

Recombination rates across porcine autosomes inferred from high-density linkage maps.

Studies of the variation in recombination rate across the genome provide a better understanding of evolutionary genomics and are also an important step towards mapping and dissecting complex traits in domestic animals. With the recent completion of the porcine genome sequence and the availability of a high-density porcine single nucleotide polymorphism (SNP) array, it is now possible to construct a high-density porcine linkage map and estimate recombination rate across the genome. A total of 416 animals were genotyped with the Porcine SNP60BeadChip, and high-density chromosome linkage maps were constructed using CRI-MAP, assuming the physical order of the Sscrofa10 assembly. The total linkage map length was 2018.79 cM, using 658 meioses and 14,503 SNPs. The estimated average recombination rate across the porcine autosomes was 0.86 cM/Mb. However, a large variation in recombination rate was observed among chromosomes. The estimated average recombination rates (cM/Mb) per chromosome ranged from 0.48 in SSC1 to 1.48 in SSC10, displaying a significant negative correlation with the chromosome sizes. In addition, the analysis of the variation in the recombination rates taking 1-Mb sliding windows has allowed us to demonstrate the variation in recombination rates within chromosomes. In general, a larger recombination rate was observed in the extremes than in the centre of the chromosome. Finally, the ratio between female and male recombination rates was also inferred, obtaining a value of 1.38, with the heterogametic sex having the least recombination.

Joint effects of porcine leptin and leptin receptor polymorphisms on productivity and quality traits.

Leptin signalling plays a fundamental role in growth, fatness and body composition. The aim of this study was to investigate the porcine LEP gene sequence in an Iberian × Landrace experimental cross to identify polymorphisms associated with productivity and quality traits. Because of the documented effects on these traits of the LEPR c.1987C>T polymorphism, the LEP and LEPR c.1987C>T polymorphisms and their interactions have been jointly investigated. The LEP gene sequencing has allowed the identification of 39 polymorphisms, eight of which are novel. Three intronic SNPs, LEP g.1382C>T, LEP g.1387C>T and LEP g.1723A>G, have been genotyped, and association analyses have been carried out. Analyses of LEP g.1387C>T, fully linked to LEP g.1382C>T, have revealed additive effects on live and carcass weights and dominant effects on several backfat thickness measurements. Novel effects of both LEP and LEPR polymorphisms on fatty acid composition in subcutaneous fat have been detected, probably mediated through the effects on fatness. The results reported here suggest that the T alleles of both LEP g.1387C>T and LEPR c.1987C>T, which are fixed in the Iberian pigs, would lead to an increase in growth, fatness and saturated fatty acid content in fat, which could be explained by an increased feed intake.

Evaluation of the porcine ACSL4 gene as a candidate gene for meat quality traits in pigs.

Long-chain acyl-CoA synthetase (ACSL) family members catalyse the formation of long-chain acyl-CoA from fatty acid, ATP and CoA, thus playing an important role in both de novo lipid synthesis and fatty acid catabolism. Previous studies in our group evaluated ACSL4 as a positional candidate gene for quantitative trait loci located on chromosome X in an Iberian × Landrace cross. A DQ144454:c.2645G>A SNP located in the 3' untranslated region of the ACSL4 gene was associated with the percentages of oleic and monounsaturated fatty acids. The aim of the present work was to evaluate the functional implication of this genetic variant. An expression analysis was performed for 120 individuals with different genotypes for the DQ144454:c.2645G>A polymorphism using real-time quantitative PCR. Differences between genotypes were identified in liver, with the ACSL4 mRNA expression levels higher in animals with the G allele than in animals with the A allele. A SNP genome-wide association study with ACSL4 relative expression levels showed significant positions on chromosomes 6 and 12. Description of positional candidate genes for ACSL4 regulation on chromosomes 6 and 12 is provided.

Building phenotype networks to improve QTL detection: a comparative analysis of fatty acid and fat traits in pigs.

Models in QTL mapping can be improved by considering all potential variables, i.e. we can use remaining traits other than the trait under study as potential predictors. QTL mapping is often conducted by correcting for a few fixed effects or covariates (e.g. sex, age), although many traits with potential causal relationships between them are recorded. In this work, we evaluate by simulation several procedures to identify optimum models in QTL scans: forward selection, undirected dependency graph and QTL-directed dependency graph (QDG). The latter, QDG, performed better in terms of power and false discovery rate and was applied to fatty acid (FA) composition and fat deposition traits in two pig F2 crosses from China and Spain. Compared with the typical QTL mapping, QDG approach revealed several new QTL. To the contrary, several FA QTL on chromosome 4 (e.g. Palmitic, C16:0; Stearic, C18:0) detected by typical mapping vanished after adjusting for phenotypic covariates in QDG mapping. This suggests that the QTL detected in typical mapping could be indirect. When a QTL is supported by both approaches, there is an increased confidence that the QTL have a primary effect on the corresponding trait. An example is a QTL for C16:1 on chromosome 8. In conclusion, mapping QTL based on causal phenotypic networks can increase power and help to make more biologically sound hypothesis on the genetic architecture of complex traits.

Partial short-read sequencing of a highly inbred Iberian pig and genomics inference thereof.

Despite dramatic reduction in sequencing costs with the advent of next generation sequencing technologies, obtaining a complete mammalian genome sequence at sufficient depth is still costly. An alternative is partial sequencing. Here, we have sequenced a reduced representation library of an Iberian sow from the Guadyerbas strain, a highly inbred strain that has been used in numerous QTL studies because of its extreme phenotypic characteristics. Using the Illumina Genome Analyzer II (San Diego, CA, USA), we resequenced ∼ 1% of the genome with average 4 × depth, identifying 68,778 polymorphisms. Of these, 55,457 were putative fixed differences with respect to the assembly, based on the genome of a Duroc pig, and 13,321 were heterozygous positions within Guadyerbas. Despite being highly inbred, the estimate of heterozygosity within Guadyerbas was ∼ 0.78 kb(-1) in autosomes, after correcting for low depth. Nucleotide variability was consistently higher at the telomeric regions than on the rest of the chromosome, likely a result of increased recombination rates. Further, variability was 50% lower in the X-chromosome than in autosomes, which may be explained by a recent bottleneck or by selection. We divided the whole genome in 500 kb windows and we analyzed overrepresented gene ontology terms in regions of low and high variability. Multi organism process, pigmentation and cell killing were overrepresented in high variability regions and metabolic process ontology, within low variability regions. Further, a genome wide Hudson-Kreitman-Aguadé test was carried out per window; overall, variability was in agreement with neutral expectations.

Evolutionary study of a potential selection target region in the pig.

Domestication, modern breeding and artificial selection have shaped dramatically the genomic variability of domestic animals. In livestock, the so-called FAT1 quantitative trait locus (QTL) in porcine chromosome 4 was the first QTL uncovered although, to date, its precise molecular nature has remained elusive. Here, we characterize the nucleotide variability of 13 fragments of ∼500 bp equally spaced in a 2 Mb region in the vicinity of the FAT1 region in a wide-diversity panel of 32 pigs. Asian and European animals, including local Mediterranean and international pig breeds, were sequenced. Patterns of genetic variability were very complex and varied largely across loci and populations; they did not reveal overall a clear signal of a selective sweep in any breed, although FABP4 fragment showed a significantly higher diversity. We used an approximate Bayesian computation approach to infer the evolutionary history of this SSC4 region. Notably, we found that European pig populations have a much lower effective size than their Asian counterparts: in the order of hundreds vs hundreds of thousands. We show also an important part of extant European variability is actually due to introgression of Asian germplasm into Europe. This study shows how a potential loss in diversity caused by bottlenecks and possible selective sweeps associated with domestication and artificial selection can be counterbalanced by migration, making it much more difficult the identification of selection footprints based on naive demographic assumptions. Given the small fragment analyzed here, it remains to be studied how these conclusions apply to the rest of the genome.

Nucleotide variability of the porcine SERPINA6 gene and the origin of a putative causal mutation associated with meat quality.

The serpin peptidase inhibitor, clade A, member 6 gene (SERPINA6), also known as corticosteroid-binding globulin or CBG, is involved in obesity and stress sensitivity. Previous studies have reported putative causal mutations within that gene in the porcine species. To characterize a hypothetical selective footprint, we have resequenced approximately 6 kb of coding and non-coding fragments in 20 pigs comprising domestic breeds and wild boars from Asia and Europe. Nucleotide variability was found to be far greater within Asian pig breeds than European breeds (π = 1% vs. 0.05%, respectively), which is consistent with pig evolutionary history. The putative causal amino acid substitution p.Gly307Arg (SNP c.919G>A) associated with meat quality (drip loss) was only detected in European domestic pig breeds, suggesting a very recent mutation that appeared after domestication in Europe. No support for positive selection was detected, as no reduction in levels of diversity surrounding the mutation was found in lean breeds with respect to wild boar.

Investigation of the oviductal glycoprotein 1 (OVGP1) gene associated with embryo survival and development in the rabbit.

An association study was performed in rabbits between early embryo survival and development, and the nonconservative SNP 12944C>G located in exon 11 and the triallellic microsatellite [(GT)(15)T(G)(5), (GT)(14)T(G)(5), and (GT)(11)T(G)(7))] located in the promoter region of the oviductal glycoprotein 1 (OVGP1) gene. We analyzed an F(2) cross of 2 lines of rabbits divergently selected for uterine capacity. A total of 172 and 159 females were slaughtered at 48 and 72 h of gestation, respectively, to determine whether OVGP1 influences ovulation rate, fertilization rate, early embryo survival, and embryonic stage of development. The results of the SNP indicated that all genotypes showed similar early embryo survival and a similar embryonic stage of development at 48 h of gestation. However, at 72 h of gestation, the GG genotype showed greater early embryo survival than the CC genotype (0.56 embryos) and their embryos presented less embryonic development. Analysis of the microsatellite was performed to ascertain the presence or absence of the allele (GT)(14)T(G)(5). At both stages of gestation, the (GT)(14)T(G)(5)/(GT)(14)T(G)(5) genotype showed greater early embryo survival (0.94 and 1.54 embryos at 48 and 72 h of gestation, respectively) and less embryonic development than the homozygous genotypes without the allele (GT)(14)T(G)(5).

Candidate gene analysis for reproductive traits in two lines of rabbits divergently selected for uterine capacity.

The objective of this work was to analyze 3 functional candidate genes for reproduction in 2 lines of rabbits divergently selected by uterine capacity. Both lines were selected for 10 generations. The selection was then relaxed until the 17th generation, when it was compounded by 61 and 63 does of the High and Low lines, respectively. We sequenced the SCGB1A1 gene, which encodes the main protein secreted by the rabbit in the uterus and seems to play an important role in implantation. We found 6 SNP in the promoter region cosegregating in 2 haplotypes in both lines with similar frequency. We also analyzed IGF1 mRNA because of its effects on embryo development, but we did not find any polymorphism between individuals of the 2 lines. The third gene analyzed was the TIMP1, which encodes a protein involved in many biological processes related to reproduction. We determined the sequence of its promoter region and found 1 SNP (g.1423A>G) segregating with different frequencies in both lines (0.60 for allele A in the High line and 0.82 for allele G in the Low line). The association study performed in an F(2) population (n = 598) generated by the cross of the 2 lines of rabbits revealed that the AA genotype had 0.88 embryos more than the GG genotype at 72 h of gestation. The difference increased to 2.23 embryos at implantation, but no difference was found between genotypes at birth. These results suggest that TIMP1 could be a candidate gene for embryo implantation and embryo survival.

Expression of progesterone receptor related to the polymorphism in the PGR gene in the rabbit reproductive tract.

The association of the 2464G > A SNP found in the promoter region of the rabbit progesterone receptor gene with progesterone receptor (PR) expression was evaluated by Western blot analysis. This SNP was associated with 2 lines divergently selected for uterine capacity, the high line selected to increase uterine capacity and the low line selected to decrease uterine capacity. Two progesterone isoforms were obtained using a commercial monoclonal antibody: the PR-B isoform described previously in rabbits, and the PR-A isoform, not described previously in rabbits. The GG genotype, the genotype more frequent in the high line, showed less PR-B and PR-A expression than the AA genotype in the oviduct (GG/AA(PR-B) = 0.81 and GG/AA(PR-A) = 0.73) and uterus (around 0.70 in both isoforms). The GA genotype showed similar PR-A expression in both tissues and also similar PR-B expression in the oviduct to the GG genotype. Conversely, the GG genotype showed less PR-B expression than the GA genotype in the uterus (GG/GA(PR-B) = 0.79). Similar expression of both PR isoforms was found in the uterus at d 2 and 3 of gestation; meanwhile, an increase of both isoforms was observed in the oviduct. Similar PR-A expression was observed in the ampulla and isthmus; meanwhile, the PR-B expression in the isthmus was double that in the ampulla.