The rate of pubertal development and weaning to estrus interval are correlated and affect reproductive efficiency of swine. Quantitative trait loci (QTL) for age of puberty, nipple number and ovulation rate have been identified in Meishan crosses on pig chromosome 10q (SSC10) near the telomere, which is homologous to human chromosome 10p15 and contains an aldo-keto reductase (AKR) gene cluster with at least six family members. AKRs are tissue-specific hydroxysteroid dehydrogenases that interconvert weak steroid hormones to their more potent counterparts and regulate processes involved in development, homeostasis and reproduction. Because of their location in the swine genome and their implication in reproductive physiology, this gene cluster was characterized and evaluated for effects on reproductive traits in swine.
Screening the porcine CHORI-242 BAC library with a full-length AKR1C4 cDNA identified 7 positive clones and sample sequencing of 5 BAC clones revealed 5 distinct AKR1C genes (AKR1CL2 and AKR1C1 through 4), which mapped to 126–128 cM on SSC10. Using the IMpRH7000rad and IMNpRH212000rad radiation hybrid panels, these 5 genes mapped between microsatellite markers SWR67 and SW2067. Comparison of sequence data with the porcine BAC fingerprint map show that the cluster of genes resides in a 300 kb region. Twelve SNPs were genotyped in gilts observed for age at first estrus and ovulation rate from the F8 and F10 generations of one-quarter Meishan descendants of the USMARC resource population. Age at puberty, nipple number and ovulation rate data were analyzed for association with genotypes by MTDFREML using an animal model. One SNP, a phenylalanine to isoleucine substitution in AKR1C2, was associated with age of puberty (p = 0.07) and possibly ovulation rate (p = 0.102). Two SNP in AKR1C4 were significantly associated with nipple number (p ≤ 0.03) and another possibly associated with age at puberty (p = 0.09).
AKR1C genotypes were associated with nipple number as well as possible effects on age at puberty and ovulation rate. The estimated effects of AKR1C genotypes on these traits suggest that the SNPs are in incomplete linkage disequilibrium with the causal mutations that affect reproductive traits in swine. Further investigations are necessary to identify these mutations and understand how these AKR1C genes affect these important reproductive traits.
The nucleotide sequence data reported have been submitted to GenBank and assigned accession numbers [GenBank:DQ474064–DQ474068, GenBank:DQ494488–DQ494490 and GenBank:DQ487182–DQ487184].
In swine, the rate of pubertal development and successful pregnancy in gilts affects the efficient management of breeding females. Selection for growth rate and leanness in modern commercial pigs has resulted in a delay in the onset of puberty . Age at puberty and weaning to estrus interval (WEI) are positively correlated  and the primary reason for culling sows is failure to return to estrus after weaning. Quantitative trait loci (QTL) have been identified for age of puberty in the pig on different chromosomes [3, 4]. One of the QTL regions is located on the long arm of pig chromosome 10 (SSC10q) near the telomere, which is homologous to human chromosome 10p15 [5, 6]. The q-arm of pig chromosome 10 also has QTL for ovulation rate  and number of nipples [7–10]. In the human, an aldo-keto reductase (AKR) gene family (AKR1C) has been identified near the telomere on chromosome 10p15 that contains at least six aldo-keto reductase family 1, member C genes . The homologous region on mouse chromosome 13 contains a cluster of eight or nine AKR1C genes . These duplicated genes maintain a high degree of sequence similarity, but differ greatly in their substrate specificity and tissue expression patterns.
The AKR superfamily are monomeric oxidoreductases that catalyze the NADP(H)-dependent reduction of a wide variety of substrates, ranging from steroids, prostaglandins, bile acids, carbohydrates and xenobiotics . AKRs are also thought to deactivate damaging reactive oxygen species like carbonyl compounds from lipids and proteins leading to their elimination [14, 15]. Aldo-keto reductases interconvert weak androgens, estrogens, progestins, mineralocorticoids and glucocorticoids to their more potent counterparts by catalyzing the reduction and oxidation of keto- and hydroxysteroids, respectively, thereby regulating a wide range of physiological processes involved in development, homeostasis and reproduction . In this manner, AKRs regulate the occupancy and transactivation of several steroid receptors in target tissues leading to transcription of hormone-responsive genes . These steroid substrates can also act directly through "non-genomic" effects, such as, formation of neuroactive steroids and activation of ion channels, G-protein-coupled receptors and several kinase signalling pathways. The products of AKR activity have been implicated in prostate disease, breast cancer, obesity, polycystic ovary disease and delay in the onset of puberty in humans [17–22]. Because of their location in the swine genome and their implication in directing reproductive physiology, this gene cluster was characterized and evaluated for affecting age at puberty in the pig.
Identification of AKR1C genes
A genomic amplicon spanning exons 4 and 5 of AKR1C4 (Genbank accession number AF473815; ) probed against one third of the porcine RPCI-44 BAC library identified 3 clones (62L11, 69L21 and 125A17). A full-length cDNA for AKR1C4 (TC200328, The Institute for Genomic Research (TIGR)  probed against one third of the porcine CHORI-242 BAC library identified seven clones (203C8, 204L24, 226I21, 264H20, 275P11, 315D4 and 319P22). Three of these clones (203C8, 226I21 and 264H20) are represented in contig 10007 of pig BAC fingerprint map . Southern analysis of a BamHI digest of 8 clones from both libraries using the same cDNA probe showed 5–9 shared bands and dissimilar banding patterns. Five of these BACs from the CHORI-242 library (203C8, 226I21, 264H20, 315D4 and 319P22) were digested with BamHI and HindIII and subcloned for sample sequencing. One BAC clone (CHORI-242-203C8) that contained at least 4 AKR1C genes was nebulized and subcloned into pBluescript to obtain more coverage. The complete AKR1C4 gene was PCR cloned by amplification using exon primers and the RPCI44-125A17 BAC clone as template. Exon sequences for AKR1CL2, AKR1C1, AKR1C2 and AKR1C3 were found in 4 of the 5 CHORI-242 BAC subclone libraries (203C8, 226I21, 315D4 and 319P22) and AKR1C4 sequence was found in CHORI-242-264H20 subclones exclusively. AKR1C2 was found in both sets of BACs in the overlapping region. Sequences corresponding to AKR1CL1 were not found in BAC subclones. The promoter region of AKR1C4 was cloned by probing BamHI/BglII digests of the RPCI BACs with a fragment containing exon 1 and part of intron 1. The promoter contained LSF (late SV40 factor), ERE (estrogen response element) and multiple SP1 and MYC sites [Genbank: DQ494489]. A 3900 bp contig from BAC clone CHORI-242-203C8 upstream of AKR1CL2 contained promoter elements CCAAT, SP1, ETS, GATA, NF-1 and ERE [Genbank: DQ494488] .
Identification of cDNAs for AKR1C
TIGR contigs were identified for four of the five AKR1C genes identified from BAC subclone sequence; only AKR1C2 was not represented in EST libraries sequenced or in contigs assembled by TIGR. These mRNAs were confirmed by RT-PCR of overlapping fragments or by completely sequencing individual clones. Amino acid and nucleotide homologies ranged from about 50–90% and 73–93%, respectively, among the pig AKR1C transcripts and the presence of active site residues common to aldo-keto reductases was conserved (Figure 1). No additional ESTs that would represent AKR1CL1 or other AKR1C genes were identified by sequence similarity analyses. Because of the high homology of these genes to all of the human AKR1C genes they were named by their relative position and conservation of amino acid sequence with human genes. The identity of pig AKR1CL2 could be confidently assigned; pig AKR1CL2 was most similar to human AKR1CL2 and mouse Akr1e1 (Figure 2) and had little similarity to other human AKR1C genes. The pig AKR1C1 followed by AKR1C2 was the most similar to human AKR1C genes (83–85%) and to bovine AKR1C genes; pig AKR1C4 was more closely related to human AKR1C4 and AKR1CL1. Except for AKR1CL2, paralogues were more similar within species than were homologues among species.
Mapping and gene organization
SNPs identified in three of the five genes (AKR1CL2, AKR1C2 and AKR1C4, [Genbank: BV102614, BV680543, AF473815, respectively]) and microsatellites SB88-91 [GenBank: DQ487182-84] found in BAC subclone sequences were linkage mapped to SSC10, positions 126–128 cM. The resolution of the map was not great enough to determine order of genes or markers. Four AKR1C genes (AKR1CL2, AKR1C2, AKR1C3 and AKR1C4) and microsatellite markers SB89, SB90 and SB91 (Table 1) were mapped using the IMpRH7000rad and IMNpRH212000rad radiation hybrid panels and anchored with flanking genes (PRKCQ, GDI2 and IDI1) and microsatellite markers SWR67 and SW2067 (Figure 3). The resolution of the IMNpRH212000rad panel was nearly the same as that of the IMpRH7000rad panel. Comparison of sample sequence data with mRNA sequences and the porcine BAC fingerprint map showed that the complete cluster of genes resides in two overlapping BAC clones that were subcloned, CHORI-242-203C8 and 264H20; clones 226I21, 315D4 and 319P22 completely overlap 203C8 (Figure 4). AKR1CL2, AKR1C1, AKR1C2 and AKR1C3 were all contained within CHORI-242 203C8 and microsatellite markers SB89, SB90 and SB91 were identified from this clone, as well. The gene order was determined to be AKR1C4, AKR1C2, AKR1C1, AKR1C3 and AKR1CL2 from centromere to telomere. BAC clone CHORI-242-264H20 contained AKR1C2, AKR1C4, the urocortin 3 (UCN3), tubulin alpha-like 3 (TUBAL3) and neuroepithelial cell transforming gene 1 (NET1) genes (Figure 4). The orientation of the AKR1C genes was determined by aligning cDNAs or BAC subclone contigs with BES on the BAC fingerprint map . This region corresponds to about 630 kb of human sequence (4.8–5.5 Mb on HSA10) but is contained in only two porcine BAC clones, average insert size of 173 kb, suggesting that this region is about half the size in the pig. A portion of this reduction is possibly due to the lack of a pig homologue to human AKR1CL1, although it is possible that AKR1CL1 was missed in the BAC sequence survey.
Primers used for mapping and tissue-specific expression of pig AKR1C genes.
Forward primer (5'→3')
Reverse primer (5'→3')
1Accession number is for genomic DNA for RH/linkage markers and mRNA sequence for RT-PCR.
Expression of pig AKR1C genes in different adult tissues
Gene-specific cDNA for the five different AKR1C genes was amplified from 16 tissues (Table 1). Six tissues expressed all five genes tested (spleen, lung, ovary, adrenal, kidney, and endometrium) and AKR1CL2 and AKR1C4 were the most widely expressed genes (Figure 5). AKR1CL2 was expressed in all tissues except pancreas and brain. Unlike human AKR1C4, pig AKR1C4 expression was not specific to the liver but was expressed in all tissues and AKR1C2 was the only other AKR1C gene expressed in brain. These expression patterns were also reflected by the number of clones and tissue source of the EST libraries from which these cDNAs were identified .
Association of SNPs with phenotypic traits
Polymorphisms were identified from RT-PCR sequences of Meishan and White composite endometrium cDNA. SNPs identified in the coding region were genotyped across generations F8 and F10 of the resource population. Additional SNPs were found by sequencing genomic DNA from animals of similar breed types and were chosen based on their frequency and potential for being a non-synonymous mutation (Table 2 and Additional File 1). Twelve SNPs were genotyped for association with reproductive traits. Three of these were in the AKR1C2 coding region, one was in AKR1CL2, and eight were the AKR1C4 gene, including one in the promoter region (Table 2). One SNP (49422_42), an isoleucine to phenylalanine substitution in AKR1C2, was associated with age of puberty (p = 0.07) and possibly ovulation rate (p = 0.102). Another SNP in AKR1C4 (49431_198) was possibly associated with age at puberty (p = 0.093; Table 3). Two other SNP were significantly associated with nipple number (p ≤ 0.03; Table 3).
SNPs identified in AKR1C genes genotyped in 1/4 Meishan F8 and F10 population.
F8 & F10
T = 0.63
T = 0.21
T = 0.81
A = 0.77
A = 0.50
A = 1.00
A = 0.62
A = 0.17
A = 0.89
G = 0.85
T = 0.55
T = 0.10
T = 0.94
G = 0.81
G = 1.00
G = 0.67
T = 0.67
T = 0.36
T = 0.86
C = 0.57
C = 0.17
C = 0.92
C = 0.78
C = 0.50
C = 1.00
T = 0.56
C = 0.56
C = 0.25
C = 0.94
T = 0.58
T = 0.50
T = 0.81
1Accession numbers beginning with BV are sequences submitted to dbSTS.
2Location refers to nucleotide position in cDNA.
3Allele frequencies were calculated from F8 and F10 animals genotyped (ranging from 238–420 animals), and founder Meishan (MS) and founder White composite (WC) animals (7–10 animals each) in the original resource population; ND = not determined.
4Polymorphisms are identified by IUB code with flanking genomic sequence.
Tests of significance and estimated effects of AKR1C genotypes on reproductive traits in F8 and F10 animals.
Age at puberty1
Number of nipples3
p > 0.154
p > 0.15
0.148 (TT = 0; CT = 0.170; CC = -0.077)
0.070 (TT = 0; AT = 13.46; AA = 6.21)
0.102 (TT = 0; AT = 1.54; AA = 1.15)
p > 0.15
p > 0.15
p > 0.15
0.064 (TT = 0; GT = 0.280; GG = 0.206)
p > 0.15
p > 0.15
0.024 (CC = 0; CT = 0.476; TT = 0.276)
p > 0.15
p > 0.15
0.030 (CC = 0; CT = 0.304; TT = 0.119)
p > 0.154
p > 0.15
0.058 (TT = 0; GT = 0.284; GG = 0.180)
p > 0.15
p > 0.15
0.106 (CC = 0; CT = 0.263; TT = 0.199)
0.093 (CC = 0; CT = 6.53; TT = -2.17)
p > 0.15
p > 0.15
1The analysis for age at puberty contained 191 observations, mean = 195 ± 27 days (mean ± std dev), range = 135 – 266.
2233 observations for ovulation rate, mean = 13.7 ± 2.4 corpora lutea, range = 9 – 21.
31144 observations for nipple number, mean = 14.6 ± 1.06 nipples, range = 11 – 19. Genotype probabilities were calculated for all F8 and F10 animals in the population using GenoProb .
4Values indicate p-value and estimated genotype effects in parentheses. The allele frequencies for F8 and F10 animals are given in Table 2.
As in other species, this gene cluster is conserved in the pig, although individual family members have undergone some sequence divergence and specialization of tissue expression, possibly due to duplication of function. Because of high sequence similarity of the genes within species and divergence among species, it is difficult to definitively assign homologues for all members of the gene family and not all genes are represented in the pig (i.e., hAKR1CL1) that are found in human or mouse. Gene duplication usually results in tandem duplication of genes or segments along the chromosome  and gene conversion can result in a species paralogues being more closely related than homologues among species . Because these genes are expressed in a multitude of tissues and the expression patterns in pig tissues do not differentiate these genes with expression patterns described in human or mouse, assignment of homologues is even more complicated. In addition, the orientation of genes in this cluster is not identical to human or mouse gene order, while gene order of flanking genes (UCN3, TUBAL3, and NET1) is conserved . As more species are fully sequenced, a clearer picture of the evolutionary process of this gene family can be drawn.
Because this pluripotent family of enzymes regulates steroid hormone action in a tissue-specific manner, they are compelling positional candidates for regulating reproductive functions [16, 28]. Steroid metabolites of AKR1C enzymes rise at the onset of puberty [22, 29] presumably due to increased substrate and enzyme activity. The onset of puberty is marked by hormonal changes directed by neuronal signals that result in activation of the hypothalamic-pituitary-gonadal axis and reproductive maturity . Central to behavioral and gonadal maturity is the release of gonadotropin releasing hormone (GnRH) followed by synthesis and secretion of luteinizing hormone (LH) and follicle stimulating hormone (FSH). An LH surge in turn is essential for stimulating the cascade of events leading to ovulation . The gonadal steroid 3alpha-hydroxy-4-pregnen-20-one (3 alpha HP) produced from progesterone by AKR1C inhibits GnRH activity on gonadotropes and suppresses FSH release from pituitary cells . Modulation of the GnRH pulse frequency could therefore cause variation in the timing of puberty. Furthermore, as GnRH regulates FSH and LH release, it is possible that AKR1C activity may affect the number of ova shed during an estrus. FSH secretion stimulates the development of antral follicles and FSH levels are greater in some lines of gilts with higher ovulation rate [33, 34].
QTL for nipple number have been identified on SSC10q and two SNPs in AKR1C4 were significantly associated with nipple number. Some of these QTL are more proximally located on SSC10q [8–10], but one identified in a Meishan/Pietrain cross maps to the same location as the AKR1C gene cluster . Because androgen and antiandrogen treatment in rodents alters nipple development and retention , and treatment of rats with an inhibitor of 5α-reductase during gestation inhibits male nipple regression , a role for dihydrotestosterone (DHT) is implicated in normal nipple development. AKRs convert DHT, a preferred substrate, to the less active androgen 3α-androstanediol, thereby regulating steroid responsiveness in target tissues  such that variation in AKR1C activity could affect nipple development in the pig. AKR1C4 is the most catalytically active isoform for DHT reduction in human  and has high expression in mammary tissue (Figure 5). Because it is ubiquitously expressed and ESTs have been identified in porcine embryonic libraries , it is probably expressed during embryonic development of mammary tissue.
Considerable support for an association of the AKR1C genotypes with nipple number was detected as well as some indication of an effect on age at puberty and possibly ovulation rate. Age of puberty and ovulation rate at a specific age are negatively correlated traits in Meishan pigs because ovulation rate increases from puberty to later estrus cycles  and animals that reach puberty earlier will have had more cycles at the time of measurement and greater number of ova shed; however, this increase in ovulation rate is less dramatic in occidental pigs. Because there was no selection performed on these animals and this is an area of increased recombination, it is assumed that recombination has greatly reduced linkage disequilibrium in this region, facilitating fine-mapping of reproductive traits in advanced generations (F8 and F10) of this population. While genetic variation in this region appears to affect reproductive traits in swine there is little evidence that the SNP markers tested are causative. Rather, the estimated effects of AKR1C genotypes show an overdominant effect on reproductive traits suggesting that these SNPs may be in incomplete linkage disequilibrium with the causal mutations or possibly there are multiple causative SNPs within this region acting in repulsion. Further investigations are necessary to identify the causal mutations and understand the role AKR1C genes have on these important reproductive traits.
Variation in the aldo-keto reductase gene cluster on pig chromosome 10 may be associated with age of puberty, nipple number and ovulation rate in swine. Future studies will determine if this variation will be useful for selection of breeding females with greater reproductive efficiency in industry populations.
cDNA synthesis, amplification and sequencing
Porcine cDNAs for aldo-keto reductases (AKRs) were identified from EST sequences deposited in GenBank and assembled at The Institute for Genomic Research (TIGR)  or by homology of porcine genomic BAC sequence to human AKRs. Tissues from adult purebred and composite breed animals were collected in RNAlater (Ambion, Austin, TX) and homogenized in Trizol (Invitrogen, Carlsbad, CA) as the source for total RNA synthesis of RT-PCR template for identification of sequence variation. Tissues from a purebred Meishan sow collected at day 25 of gestation and testis from a mature boar were used for differential expression of AKR1C transcripts. cDNA was synthesized with M-MLV reverse transcriptase (Promega, Madison, WI) using 2 μg of total RNA from endometrium, placenta, testis, ovary, liver, lung, adrenal, kidney, spleen, pituitary, hypothalamus, brain, pancreas, small intestine, skeletal muscle and mammary tissue. These reactions were run for 35 cycles with 20 ng of template as described below. Sequences were extended using 3'-RACE or RT-PCR from exon sequences identified in BAC subclones. Full-length cDNA clones were obtained by iterative screening and self-ligation of inverse PCR (SLIP)  of the MARC 1PIG and 2PIG primary libraries before normalization .
BAC screens and subclone libraries
Filters from the RPCI-44 and CHORI-242 porcine BAC libraries were screened using a random-primed nearly full-length cDNA of porcine AKR1C4 as probe. The probe was prepared by PCR of the MARC 2PIG library using primers in exons 1 and 9 and radioactive random-primed labeled (Megaprime DNA Labeling System, Amersham, Piscataway, NJ). Positive clones were grown in 100 ml cultures and processed for BAC DNA using a midi-prep column (Marligen Biosciences, Ijamsville, MD). The BACs were digested with BamHI and HindIII separately and subcloned into pBluescript. One 96-well plate of clones was processed for each BAC and enzyme combination and sequenced with T3 and T7 primer. One CHORI-242 BAC (203C8) was sheared with a nebulizer (Invitrogen, Carlsbad, CA), end-repaired, then cloned into pBluescript and four 384-well plates of clones were sequenced.
PCR and sequencing
Primer pairs for amplification of genomic DNA were designed from porcine AKR1C cDNA sequences and genomic sequence obtained during this study using Primer 3; code available at the Primer 3 Software website . PCR was performed in a PTC-225 DNA engine (MJ Research Inc, Watertown, Mass) using 0.25 U Hot Star®Taq polymerase (Qiagen, Valencia, CA, USA), 1× of supplied buffer, 1.5 mM MgCl2, 200 μM dNTPs, 0.8 μM each primer, and 100 ng of genomic DNA in 25 μl reactions. Five μl of the PCR reaction was electrophoresed in 1.5% agarose gels to determine quality of amplification and the remainder was prepared for sequencing. Chromatograms were imported into the MARC database, bases called with Phred, assembled into contigs with Phrap, polymorphisms identified using Polyphred, and assessed using Consed .
SNPs were mapped using a primer extension assay on the Sequenom MassArray™ system (San Diego, CA, USA). Ten μl PCR reactions contained 10 ng of genomic DNA, 0.25 U HotStar Taq, 1× of supplied buffer, 1.5 mM MgCl2, 200 μM dNTPs, and 0.4 μM forward and reverse tailed primers. The primer extension reaction used 0.6 μM of probe primer and was performed according to the manufacturer's recommendations for hME chemistry (Sequenom, San Diego, CA, USA).
Radiation hybrid and linkage mapping
Genes were mapped using the 118-clone INRA-University of Minnesota porcine Radiation Hybrid IMpRH7000rad panel and the 90-clone IMNpRH212000rad panel [42, 43]. Primers used were described above, or designed from sequences obtained from subclone sequences. Amplifications were performed in 15 μl PCR in duplicate using 12.5 ng panel DNA, 1.5 mM MgCl2, 200 μM dNTPs, 1 μM each primer, 0.25 U Hot Star®Taq and 1× of supplied buffer. The PCR mixture was held at 94°C for 15 min, and cycled 40 times at 94°C for 20 sec, held at the indicated annealing temperature for 30 sec and extension at 72°C for 45–60 sec, followed by a final extension at 72°C for 5 min. One half of the reaction was loaded on 2% agarose gels and manually genotyped. Data were analyzed for two-point and multipoint linkage with the IMpRH mapping tool  and submitted to the IMpRH database . Carthagene [46, 47] was used to estimate multipoint marker distance and order using all public markers on chromosome 10 in the IMpRH database  for the IMpRH7000rad panel and those developed in this study for the IMNpRH212000rad panel to approximate position of mapped markers. Markers run on the IMpRH7000rad panel have been submitted to the IMpRH public database.
Linkage analyses were performed as described  where TWOPOINT analyses were used to indicate the chromosome linkage group and the ALL, FLIPS and FIXED options were used to determine the multipoint position of the marker (CRIMAP v2.4). Multipoint locations for all mapped markers are based on the latest published swine genetic map .
Animals and resource population
Genomic DNA from parents of the MARC reference family (a White composite boar and seven crossbred sows) was used to identify SNPs. Phenotypes and genotypes were collected from animals of generations 8 and 10 (F8 and F10) produced from the original resource population  used to identify QTL. The original animals were from reciprocal backcrosses of 10 purebred Meishan and 10 White composite (composed of Chester White, Landrace, Large White and Yorkshire). The half-Meishan F4 animals were crossed with a Landrace-Large White composite reducing the Meishan influence to one quarter. This population has been inter se mated since the F5 generation and maintained at 30 litters/year. Procedures for the handling of animals complied with those specified in the Guide for the Care and Use of Agricultural Animals in Agricultural Research and Teaching (1999) [1st rev. ed. Savoy, IL: Federation of Animal Science Societies; 1999].
Phenotypic data and statistical methods
The F8 gilts were observed for first estrus beginning when the oldest gilts reached 120 days of age as described by Rohrer et al. , and ovulation rates were determined by counting corpora lutea on the ovaries at slaughter after the third estrus in gilts from the F8 and F10 generations. Genotype probabilities were calculated for all animals using an extended version of GenoProb . Association analyses between phenotypes and genotypes were conducted using MTDFREML . The model fitted included fixed effects for contemporary group and regressions on SNP marker genotype probabilities (i.e., probability of an animal being aa, aA, or AA). Random effects included the animal's polygenic breeding value and residual error. Each marker was analyzed separately. No adjustments to reported p-values were made for multiple comparisons. Twelve SNP markers were analyzed for three traits; age of puberty (191 observations), ovulation rate (233 observations) and nipple number (1144 observations).
AKR1C peptide sequences were aligned with ClustalX, the distance matrix constructed with the program PRODIST using a Dayhoff PAM matrix model, and a neighbor-joining tree constructed in PHYLIP (v 3.65). The tree was constructed with the program NEIGHBOR and rooted at a mid-point with the program RETREE. Bootstrap values were derived from 1000 pseudo-datasets generated in SEQBOOT. The tree was viewed in TREEVIEW (v 1.6.6).
The authors thank Sue Hauver, Bree Quigley and Kris Simmerman for excellent technical assistance, Jim Wray for sequence submissions, Michael Clawson for phylogenetic analysis and Sherry Kluver for manuscript preparation. Mention of trade names or commercial products is solely for the purpose of providing information and does not imply recommendation, endorsement or exclusion of other suitable products by the U.S. Department of Agriculture.
USDA-ARS, U.S. Meat Animal Research Center
Rauw WM, Kanis E, Noordhuizen-Stassen EN, Grommers FJ: Undesirable side effects of selection for high production efficiency in farm animals: a review.Livestock Production Science 1998, 56:15–33.View Article
Sterning M, Rydhmer L, Eliasson-Selling L: Relationships between age at puberty and interval from weaning to estrus and between estrus signs at puberty and after the first weaning in pigs.Journal of Animal Science 1998, 76:353–359.PubMed
Cassady JP, Johnson RK, Pomp D, Rohrer GA, Van Vleck LD, Spiegel EK, Gilson KM: Identification of quantitative trait loci affecting reproduction in pigs.Journal of Animal Science 2001, 79:623–633.PubMed
Rohrer GA, Ford JJ, Wise TH, Vallet JL, Christenson RK: Identification of quantitative trait loci affecting female reproductive traits in a multigeneration Meishan-White composite swine population.Journal of Animal Science 1999, 77:1385–1391.PubMed
Nonneman DJ, Rohrer GA: Comparative mapping of a region on chromosome 10 containing QTL for reproduction in swine.Animal Genetics 2003, 34:42–46.View ArticlePubMed
Nonneman D, Rohrer GA: Comparative mapping of human chromosome 10 to pig chromosomes 10 and 14.Animal Genetics 2004, 35:338–343.View ArticlePubMed
Dragos-Wendrich M, Moser G, Bartenschlager H, Reiner G, Geldermann H: Linkage and QTL mapping for Sus scrofa chromosome 10.Journal of Animal Breeding and Genetics 2003, 120:82–88.View Article
Hirooka H, de Koning DJ, Harlizius B, van Arendonk JA, Rattink AP, Groenen MA, Brascamp EW, Bovenhuis H: A whole-genome scan for quantitative trait loci affecting teat number in pigs.Journal of Animal Science 2001, 79:2320–2326.PubMed
Rodriguez C, Tomas A, Alves E, Ramirez O, Arque M, Munoz G, Barragan C, Varona L, Silio L, Amills M, Noguera JL: QTL mapping for teat number in an Iberian-by-Meishan pig intercross.Animal Genetics 2005, 36:490–496.PubMed
Rohrer GA: Identification of quantitative trait loci affecting birth characters and accumulation of backfat and weight in a Meishan-White Composite resource population.Journal of Animal Science 2000, 78:2547–2553.PubMed
Vergnes L, Phan J, Stolz A, Reue K: A cluster of eight hydroxysteroid dehydrogenase genes belonging to the aldo-keto reductase supergene family on mouse chromosome 13.Journal of Lipid Research 2003, 44:503–511.View ArticlePubMed
Penning TM: Hydroxysteroid dehydrogenases and pre-receptor regulation of steroid hormone action.Human Reproduction Update 2003, 9:193–205.View ArticlePubMed
Fujii J, Iuchi Y, Okada F: Fundamental roles of reactive oxygen species and protective mechanisms in the female reproductive system.Reproductive Biology and Endocrinology 2005, 3:43.View ArticlePubMed
Zhang D, Tai LK, Wong LL, Chiu LL, Sethi SK, Koay ES: Proteomic study reveals that proteins involved in metabolic and detoxification pathways are highly expressed in HER-2/neu-positive breast cancer.Molecular & Cellular Proteomics 2005, 4:1686–1696.View Article
Bauman DR, Steckelbroeck S, Penning TM: The roles of aldo-keto reductases in steroid hormone action.Drug News Perspective 2004, 17:563–578.View Article
Stanbrough M, Bubley GJ, Ross K, Golub TR, Rubin MA, Penning TM, Febbo PG, Balk SP: Increased expression of genes converting adrenal androgens to testosterone in androgen-independent prostate cancer.Cancer Research 2006, 66:2815–2825.View ArticlePubMed
Lewis MJ, Wiebe JP, Heathcote JG: Expression of progesterone metabolizing enzyme genes (AKR1C1, AKR1C2, AKR1C3, SRD5A1, SRD5A2) is altered in human breast carcinoma.BMC Cancer 2004, 4:27.View ArticlePubMed
Lord SJ, Mack WJ, Van Den Berg D, Pike MC, Ingles SA, Haiman CA, Wang W, Parisky YR, Hodis HN, Ursin G: Polymorphisms in genes involved in estrogen and progesterone metabolism and mammographic density changes in women randomized to postmenopausal hormone therapy: results from a pilot study.Breast Cancer Research 2005, 7:R336–344.View ArticlePubMed
Blouin K, Blanchette S, Richard C, Dupont P, Luu-The V, Tchernof A: Expression and activity of steroid aldoketoreductases 1C in omental adipose tissue are positive correlates of adiposity in women.American Journal of Physiology Endocrinology & Metabolism 2005, 288:E398–404.View Article
Qin K, Ehrmann DA, Cox N, Refetoff S, Rosenfield RL: Identification of a functional polymorphism of the human type 5 17beta-hydroxysteroid dehydrogenase gene associated with polycystic ovary syndrome.Journal of Clinical Endocrinology and Metabolism 2006, 91:270–276.View ArticlePubMed
Rittner HL, Lee PD, Blum WF, Doerr HG, Steiss J, Kreuder J, Rascher W, Kiess W: Developmental patterns of serum 3 alpha-androstanediol glucuronide.Journal of Endocrinloogical Investigation 1997, 20:138–143.
Frith MC, Hansen U, Weng Z: Detection of cis-element clusters in higher eukaryotic DNA.Bioinformatics 2001, 17:878–889.View ArticlePubMed
Copley RR, Goodstadt L, Ponting C: Eukaryotic domain evolution inferred from genome comparisons.Current Opinion in Genetics & Development 2003, 13:623–628.View Article
Winter EE, Ponting CP: Mammalian BEX, WEX and GASP genes: coding and non-coding chimaerism sustained by gene conversion events.BMC Evolutionary Biology 2005, 5:54.View ArticlePubMed
Penning TM, Burczynski ME, Jez JM, Hung H-K, Ma H, Moore M, Palackal N, Ratnam K: Human 3α-hydroxysteroid dehydrogenase isoforms (AKR1C1-AKR1C4) of the aldo-keto-reductase superfamily: functional plasticity and tissue distribution reveals roles in the inactivation and formation of male and female sex hormones.Biochemical Journal 2000, 351:67–77.View ArticlePubMed
Riddick LM, Garibaldi LR, Wang ME, Senne AR, Klimah PE, Clark AT, Levine LS, Oberfield SE, Pang SY: 3 alpha-androstanediol glucuronide in premature and normal pubarche.Journal of Clinical Endocrinology & Metabolism 1991, 72:46–50.View Article
Sisk CL, Foster DL: The neural basis of puberty and adolescence.Nature Neuroscience 2004, 7:1040–1047.View ArticlePubMed
Knox RV: Recruitment and selection of ovarian follicles for determination of ovulation rate in the pig.Domestic Animal Endocrinology 2005, 29:385–397.View ArticlePubMed
Griffin LD, Mellon SH: Biosynthesis of the neurosteroid 3 alpha-hydroxy-4-pregnen-20-one (3 alpha hp), a specific inhibitor of FSH release.Endocrinology 2001, 142:4617–4622.View ArticlePubMed
Ford JJ, Zimmerman DR, Wise TH, Leymaster KA, Christenson RK: Increased plasma follicle-stimulating hormone concentrations in prepubertal gilts from lines selected for increased number of corpora lutea.Journal of Animal Science 2001, 79:1877–1882.PubMed
Knox RV, Vatzias G, Naber CH, Zimmerman DR: Plasma gonadotropins and ovarian hormones during the estrous cycle in high compared to low ovulation rate gilts.Journal of Animal Science 2003, 81:249–260.PubMed
Foster PM, McIntyre BS: Endocrine active agents: implications of adverse and non-adverse changes.Toxicol Pathol 2002, 30:59–65.View ArticlePubMed
Imperato-McGinley J, Binienda Z, Gedney J, Vaughan ED Jr: Nipple differentiation in fetal male rats treated with an inhibitor of the enzyme 5 alpha-reductase: definition of a selective role for dihydrotestosterone.Endocrinology 1986, 118:132–137.View ArticlePubMed
Christenson RK: Ovulation rate and embryonic survival in Chinese Meishan and white crossbred pigs.Journal of Animal Science 1993, 71:3060–3066.PubMed
Hoskins RA, Stapleton M, George RA, Yu C, Wan KH, Carlson JW, Celniker SE: Rapid and efficient cDNA library screening by self-ligation of inverse PCR products (SLIP).Nucleic Acids Research 2005, 33:e185.View ArticlePubMed
Fahrenkrug SC, Smith TP, Freking BA, Cho J, White J, Vallet J, Wise T, Rohrer G, Pertea G, Sultana R, Quackenbush J, Keele JW: Porcine gene discovery by normalized cDNA-library sequencing and EST cluster assembly.Mammalian Genome 2002, 13:475–478.View ArticlePubMed
Fahrenkrug SC, Freking BA, Smith TP, Rohrer GA, Keele JW: Single nucleotide polymorphism (SNP) discovery in porcine expressed genes.Animal Genetics 2002, 33:186–195.View ArticlePubMed
Yerle M, Pinton P, Robic A, Alfonso A, Palvadeau Y, Decros C, Hawken R, Alexander L, Beattie C, Schook L, Milan D, Gellin J: Construction of a whole-genome radiation hybrid panel for high-resolution gene mapping in pigs.Cytogenetics and Cell Genetics 1998, 82:182–188.View ArticlePubMed
Yerle M, Pinton P, Delcros C, Aranal N, Milan D, Robic A: Generation and characterization of a 12,000-rad radiation hybrid panel for fine mapping in pig.Cytogenetic and Genome Research 2002, 97:219–228.View ArticlePubMed
Milan D, Hawken R, Cabau C, Leroux S, Genet C, Lahbib Y, Tosser G, Robic A, Hatey F, Alexander L, Beattie C, Schook L, Yerle M, Gellin J: IMpRH server: an RH mapping server available on the Web.Bioinformatics 2000, 16:558–559.View ArticlePubMed
Thallman RM: User's Manual for GenoProb Version 2.000.USDA-ARS, U.S. Meat Animal Research Center 2002.
Boldman KG, Kriese LA, Van Vleck LD, Van Tassel CP, Kachman SD: A manual for the use of MTDFREML. A set of programs to obtain estimates of variance and covariances [DRAFT].USDA-ARS. Washington, DC: U.S. Gov. Print. Office 1995.
Heredia VV, Kruger RG, Penning TM: Steroid-binding site residues dictate optimal substrate positioning in rat 3alpha-hydroxysteroid dehydrogenase (3alpha-HSD or AKR1C9).Chem Biol Interact 2003, 143–144:393–400.View ArticlePubMed
Jez JM, Bennett MJ, Schlegel BP, Lewis M, Penning TM: Comparative anatomy of the aldo-keto reductase superfamily.Biochem J 1997, 326:625–636.PubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.