Next Article in Journal
Added Fructose in Non-Alcoholic Fatty Liver Disease and in Metabolic Syndrome: A Narrative Review
Previous Article in Journal
Selenium Deficiency during Pregnancy in Mice Impairs Exercise Performance and Metabolic Function in Adult Offspring
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Nutrients
Volume 14
Issue 5
10.3390/nu14051126
nutrients-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Soy Formula Is Not Estrogenic and Does Not Result in Reproductive Toxicity in Male Piglets: Results from a Controlled Feeding Study

by
Martin J. J. Ronis
1,*,
Horacio Gomez-Acevedo
2,
Kartik Shankar
3,
Leah Hennings
4,
Neha Sharma
5,
Michael L. Blackburn
5,
Isabelle Miousse
6,
Harry Dawson
7,
Celine Chen
7,
Kelly E. Mercer
5,8 and
Thomas M. Badger
5,8
1
Department of Pharmacology & Experimental Therapeutics, Louisiana State University Health Sciences Center, 1901 Perdido Str., New Orleans, LA 70112, USA
2
Department of Biomedical Informatics, University of Arkansas for Medical Sciences, Little Rock, AR 72205, USA
3
Department of Pediatrics-Nutrition, University of Colorado Anschutz Medical Center, Aurora, CO 80045, USA
4
Department of Pathology, University of Arkansas for Medical Sciences, Little Rock, AR 72205, USA
5
Arkansas Children’s Nutrition Center, Little Rock, AR 72202, USA
6
Department of Biochemistry, University of Arkansas for Medical Sciences, Little Rock, AR 72205, USA
7
USDA ARS Nutrition Center, Diet Genomics and Immunology Laboratory, Beltsville, MD 20705, USA
8
Department of Pediatrics, University of Arkansas for Medical Sciences, Little Rock, AR 72205, USA
*
Author to whom correspondence should be addressed.
Nutrients 2022, 14(5), 1126; https://0-doi-org.brum.beds.ac.uk/10.3390/nu14051126
Submission received: 19 January 2022 / Revised: 17 February 2022 / Accepted: 28 February 2022 / Published: 7 March 2022
(This article belongs to the Section Phytochemicals and Human Health)
Browse Figures
Versions Notes

Abstract

:
Soy infant formula which is fed to over half a million infants per year contains isoflavones such as genistein, which have been shown to be estrogenic at high concentrations. The developing testis is sensitive to estrogens, raising concern that the use of soy formulas may result in male reproductive toxicity. In the current study, male White-Dutch Landrace piglets received either sow milk (Sow), or were provided milk formula (Milk), soy formula (Soy), milk formula supplemented with 17-beta-estradiol (2 mg/kg/d) (M + E2) or supplemented with genistein (84 mg/L of diet; (M + G) from postnatal day 2 until day 21. E2 treatment reduced testis weight (p < 0.05) as percentage of body weight, significantly suppressed serum androgen concentrations, increased tubule area, Germ cell and Sertoli cell numbers (p < 0.05) relative to those of Sow or Milk groups. Soy formula had no such effects relative to Sow or Milk groups. mRNAseq revealed 103 differentially expressed genes in the M + E2 group compared to the Milk group related to endocrine/metabolic disorders. However, little overlap was observed between the other treatment groups. These data suggest soy formula is not estrogenic in the male neonatal piglet and that soy formula does not significantly alter male reproductive development.

1. Introduction

Breastfeeding is the best choice for infant nutrition, and the American Academy of Pediatrics (AAP) suggests exclusive breastfeeding [1,2]. However, many women are unable to breastfeed and utilize infant formulas. AAP recommends feeding infant formulas for the first year. Of those infants fed formula, cow’s milk formula accounts for 69%, followed by soy formula at 12% [3]. Soy foods have been safely consumed in China and other Far Eastern countries for thousands of years. Moreover, soy infant formulas have been widely utilized for over 50 years, and many clinical studies report normal growth and development in infants consuming soy formula [4,5,6,7,8]. Nevertheless, some have suggested that soy formula use can result in adverse reproductive health outcomes out of concern for estrogen-like effects resulting from the presence of isoflavone polyphenols such as genistein bound to the soy protein [9,10,11,12]. As a result, in 2011, the National Toxicology Program raised its level of concern about the potential reproductive toxicity of soy formula from negligible to minimal [13]. Isoflavones are similar structurally to 17β estradiol (E2) [14,15]. Genistein has been shown to be estrogenic in certain species and/or at high concentrations in cell culture [5,14,16,17]. In addition, genistein and the other soy-associated isoflavones, daidzein and equol, have the capacity to bind weakly to estrogen receptors (ERs) α (ESRA) and β (ESRB), with a preference for ESRB [18,19,20]. These compounds are considered endocrine-disrupting compounds, which have the potential to interfere with reproductive development and fertility [21].
Testicular Sertoli, Leydig, and germ cells express ESRA and ESRB [18,22,23], and exposure to estrogens during early male development is well-known to interfere with growth and development of the testis, suppress androgen production, and reduce fertility [24]. With respect to soy isoflavones, Wisniewski et al. [25] reported that dietary genistein fed through gestation and lactation reduces testis size and testosterone concentrations in rats and disrupts reproductive behavior in adulthood. More recently, Meena et al. [26] have reported prenatal, i.p., genistein exposure-reduced fertility in adult male rats associated with a deterioration of testicular architecture, dose-dependent reductions in sperm production, and decreased plasma testosterone. In addition, Sharpe and co-workers [11,12] reported reduced serum testosterone and increased Leydig cell numbers in marmoset monkeys following postnatal feeding of soy infant formula, although no lasting effects were observed at adulthood.
In contrast, we have reported no effects on development of male reproductive organs or effects on fertility in rats fed soy protein isolate (SPI), the sole protein source in soy infant formula, throughout development [5,27]. Moreover, we have observed different gene expression profiles in many estrogen-responsive organs of adult and prepubertal rats and in the mammary glands of female neonatal piglets after feeding SPI and soy infant formula relative to those elicited by E2 [28,29,30,31,32,33,34]. These data suggest that the endocrine effects of soy products and E2 differ significantly.
The neonatal piglet is a good surrogate for many aspects of infant physiology and nutrition [34,35,36]. Piglets can be fed commercially available infant formulas. In the current study, we utilized this model to test the hypothesis that soy formula or supplementation of cow’s milk formula with genistein at levels found in soy formula, would activate ER-signaling pathways in the testis, resulting in feminization, and the results determined whether that hypothesis should be rejected. Therefore, groups of male piglets aged two days were fed either cow’s milk formula, cow’s milk formula supplemented with genistein at levels producing serum concentrations similar to those of soy formula-fed infants, soy formula, or sow milk via bowl feeding until postnatal day (PND) 21. To serve as positive control, an additional group was fed cow’s milk formula supplemented with E2 to activate ER signaling pathways for comparison purposes.

2. Materials and Methods

2.1. Animal Experiments

Piglets were housed at the Arkansas Children’s Hospital Research Institute (ACHRI). Animal maintenance and experimental treatments were conducted in accordance with the ethical guidelines for animal research established and approved by the Institutional Animal Care and Use Committee at the University of Arkansas for Medical Sciences. A detailed description of this protocol has been previously reported, and additional data from this study has been previously published [34]. In brief, litters containing 8–11 piglets from Dutch Landrace Duroc sows fed a soy-free diet were used. Piglets were selected at random from many different litters born on the pig farm on the same day to form the different diet groups and to avoid litter effects. At birth, male piglets were allowed to suckle at the farm for 2 days, then transferred to the animal facility at ACHRI. Piglets were then randomly assigned to one of 5 groups and fed commercially obtained sow milk by bowl feeding (Sow, n = 5), a cow’s milk-based formula (Milk, n = 6) (Similac Advance powder, Ross Products, Abbott Laboratories, Abbott Park, IL, USA), a soy-based formula (Soy, n = 6) (Enfamil Prosobee Lipil powder, Mead Johnson Nutritionals, River Point, Chicago, IL, USA), the same cow’s milk-based formula (Similac Advance, Ross Products, Abbott Laboratories, Abbott Park, IL, USA) supplemented with 2 mg/kg/d E2 (M + E2, n = 6) as a positive control, or the Similac Advance with pure genistein (84 mg/L diet) (M + G, n = 6) to produce serum levels comparable to those found in soy formula-fed piglets and infants for 21 days. Formula diets were modified to meet the needs of the growing pig. A detailed description of the nutrient composition of the sow milk and formula diets has previously been published [35]. Piglets were all housed individually and trained to drink from small bowls on a fixed schedule: the first week was every 2 h, the second week was every 4 h, and the third week was every 6 h, to provide 1.04 MJ/kg/d until day 21. In the M + E2 group, E2 was added to the food bowls at each feeding to achieve the serum E2 concentrations achieved in peripubertal piglets (15.6–30 pmol/L) [34]. In the M + G group, genistein was added to the food bowls to achieve a concentration (84 mg/L of diet) that was comparable to the genistein concentrations in SPI, 54 mg/L, (Dupont Nutrition and Biosciences, St. Louis, MO, USA). Piglets were euthanized after anesthetization, 6–8 h after the final feeding period. Blood was collected, and serum was stored at −20 °C until analysis for isoflavones and hormones. Testis, prostate, and numerous other tissues were weighed, and tissue samples were fixed in formalin or snap-frozen in liquid nitrogen and stored at −70 °C until analysis.

2.2. Serum Soy Isoflavones

Serum isoflavones were extracted and analyzed by LC-MS as previously published [34].

2.3. Serum Hormones

Follicle stimulating hormone (FSH) and luteinizing hormone (LH) and testosterone concentrations were measured using ELISA kits, NB-ES0109, BG-POR11481, and NB-ES104, respectively, from Novateinbio (Cambridge, MA, USA). Androstenedione and progesterone were measured using ELISA kits 07-109202 and 07-170102 from MP Biomedicals (Solon, OH, USA). Dihydrotestosterone (DHT) was measured using kit DSL-9600 from Diagnostic Systems Labs Inc. (Webester, TX, USA), and dehydroepiandrosterone (DHEA) and DHEA-sulfate (DHEA-S) concentrations were also measured by ELIA kits PN0980 and PN0953 from TSZ (Framingham, MA, USA). E2 concentrations were measured using a customized MULTI-SPOT®96-well Custom Steroid Hormone Panel Plate (Lot # Z0055407, Meso Scale Diagnostics, Rockville, MD, USA) using the MESO SQ120 QuickPlate instrument.

2.4. Testis Morphology and Cellularity

Bouin’s fixed testis were sectioned, paraffin-embedded, and stained with H&E. Testicular morphology (tubular volume/total volume) and cellularity (numbers of Leydig, Sertoli, and germ cells) were assessed in a blinded fashion in 3 fields from 4 different 1.5 × 1.5 cm sections per animal by a veterinary pathologist (L.H.) as described by Lanning et al. [36].

3. Gene Expression Analysis by mRNAseq

Total RNA was isolated from frozen testes through RNeasy-mini columns (Qiagen, Valencia, CA, USA), including on-column DNAse digestion. cDNA library construction from the mRNA was carried out using NEB-Next reagents (New England Biolabs, Ipswich, MA, USA) following the manufacturer’s instructions. Sequencing of the cDNA Libraries was performed with the Genome Analyzer IIX using the TruSeq v5 reagents (Illumina Inc., San Diego, CA, USA). For details of cDNA library preparation and sequencing procedures, please refer to Ronis et al. [37]. On average, approximately 10 million 36 bp single-end reads per sample were generated. The FASTQ files with raw data were submitted to the National Center for Biotechnology Information (NCBI) Sequence Read Archive database (SRA) under the accession numbers GSE168627.
Nucleotides below Q25 or reads containing more than two ambiguous nucleotides were removed before sequence alignments performed by the CLC Genomics Workbench version 20.01 (QIAGEN Bioinformatics, Redwood City, CA, USA). For gene expression calculation, reads were first mapped to a custom, manually curated nonredundant (NR) 13,547 gene library [38]. The sequences of these genes are in the Porcine Translational Research Database (PTR, http://tinyurl.com/hxxq.3ur, accessed on 10 January 2022). The remaining unmapped reads were mapped again to the Ensembl pig genome build 11.1 (WG) in search of expressed genes that are not covered by the PTR database.
Transcriptomes built from the mapping results were subjected to differential expression analysis. Group correlations estimated by the Pairwise method (JMP Genomics, version 9, SAS, Cary, NC, USA) showed that sample “Sow 2” had correlation values below 0.9. Therefore, it was excluded from further analysis. The statistical analyses were carried out with the exact tests from the Bioconductor package “edgeR” (version 3.30.0; run on RStudio, version 4.0.3, Boston, MA, USA) [39]. Genes were considered differentially expressed with the thresholds of a false discovery rate (FDR) ≤ 0.05 and an absolute fold change ≥ 2.0. Pathway analysis was carried out with DAVID (https://david.ncifcrf.gov/) (accessed on 10 January 2022) to explore KEGG pathways.

4. mRNAseq Confirmation by Real-Time RT-PCR

Real-time RT-PCR was used to confirm expression changes of some genes regulated by E2, feeding of soy vs. milk formula or regulated by genistein. For each sample, total RNA was reverse-transcribed into cDNA and subsequent real-time RT-PCR analysis was carried out using SYBR green and ABI 7500 sequence detection system (Applied Biosystems, Waltham, MA) as described previously [35]. Gene expression was analyzed for each piglet using the 2−∆∆ct method relative to a housekeeping gene 18S gene amplification and expressed as fold change compared with the Sow group. Primer sequences for: Star F 5′-TGCTCAGCATTGACCTCAAGGGAT-3′, R 5′-TTTCGAAGGTGATTGGCAAACTCC-3′; Cyp17a1 F 5′-GGTGCCCAGACCACAATTTA-3′, R 5′-CTTTACCACAGAGGCAGAAGTC3′, and Cyp19 F 5′GAATTCATGAGGGTCTGGATAGG-3′, R 5′-CCAAATCGGCACGTGTAATG-3′.

5. Statistical Analysis

All experimental outcomes other than the data from mRNAseq analysis are summarized as mean ± SE. Differences in isoflavone concentrations between Soy and M + G groups were determined by 2-tailed Student’s t-test, p < 0.05. With respect to body weight, tissue weight, testis morphology, hormone concentrations and real-time RT-PCR gene expression, comparisons between sow and the formula diets and between MILK and M + E2 or M + G supplemented groups were accomplished by one-way ANOVA followed by a Student-Newman-Keuls post hoc analysis p < 0.05 except in cases of non-normal distribution where statistical differences between groups was verified by one-way ANOVA of Ranks followed by Dunn’s post hoc analysis.

6. Results

6.1. Serum Soy Isoflavone Concentrations Detected in the Soy- and M + G-Fed Piglets

Serum isoflavone concentrations in male piglets from this study were determined by LC-MS as previously described [34]. In the Soy-fed group we detected, genistein, daidzein and small amounts of glycitein, and o-desmethylangolansin (DMA) in the serum comparable to circulating soy isoflavones, reported in infants aged 4 months consuming soy formula [15]. Genistein values were 1712 ± 212 pg/mL and daidzein values were 1673 ± 374 pg/mL (Table 1). As observed in human neonates fed soy formula, levels of serum equol (conjugated or unconjugated) were undetectable in all diet groups nor was there any evidence for the presence of the equol precursor dihydrodaidzein. We did detect genistein in the M + G group (1898 ± 515 pg/mL). Genistein concentrations between Soy and M + G groups did not differ significantly.

6.2. Effects of Formula Feeding, Genistein and E2 on Body and Organ Weights

Differences in body weight and testis and prostate weights in the different groups of male piglets on PND 21 at the end of feeding are shown in Table 2. Direct comparisons were made of the effects of diet between sow-fed and formula-fed piglets and between milk formula and the effects of supplementation with pure genistein or E2. No significant differences were observed in body weight. Consumption of soy formula resulted in a small increase in absolute testis weight relative to those in the Sow and Milk groups (p < 0.05). However, significance was lost when corrected for body weight. E2 addition to cow’s milk formula significantly reduced testis size in both absolute terms and when corrected for body weight compared to all other groups (p < 0.05). In contrast, neither diet nor addition of genistein or E2 affected piglet prostate weight.

6.3. Effects of Formula, Genistein, and E2 on Endocrine Parameters in Neonatal Male Piglets

Differences between treatment groups in serum gonadotropins and androgen concentrations are shown in Table 3. Effects on additional gonadal and adrenal steroids are shown in Table 4. Both feeding soy formula and addition of genistein to milk formula resulted in significant increases in serum FSH concentrations relative to soy milk- or cow’s milk-fed groups (p < 0.05). In contrast, addition of E2 to milk formula had the opposite effect and the Milk + E2 group had lower serum FSH than the Milk + Genistein group (p < 0.05). (Table 4).
No significant effects of formula feeding, genistein, or E2 were observed on serum LH concentrations, which were low and quite variable. The testosterone precursor androstenedione was slightly lower in cow’s milk formula-fed piglets compared to that in sow milk- or soy milk-fed animals (p < 0.05), and this was reduced further by addition of E2. Likewise, serum testosterone concentrations were significantly suppressed by E2 treatment relative to all other treatment groups (p < 0.05). No significant effects of formula feeding, genistein, or E2 treatment were observed on serum concentrations of DHT (Table 3). As far as effects of diet on other steroids are concerned, feeding either formula reduced serum progesterone concentrations relative to that when feeding with sow milk (p < 0.05) (Table 4). No significant effects of diet were observed on serum concentrations of E2 or estrone or the adrenal steroid DHEA (Table 4). However, soy formula did slightly elevate circulating concentrations of the major prepubertal steroid hormone DHEA sulfate (DHEA-S) (p < 0.05) (Table 4). Addition of genistein to cow’s milk formula had no significant effects on any of these steroids. However, as expected, addition of E2 to milk formula resulted in substantial increases in circulating E2 and its metabolite estrone (Table 4).

6.4. Effects of Formula Feeding, Genistein, and E2 on Testicular Morphology and Cellularity in Neonatal Male Piglets

Effects of neonatal diet on testis morphology are shown in Figure 1. Although tubule area/total area was reduced (p < 0.05) in soy formula-fed piglets relative to those of the sow milk-fed group, this parameter did not differ significantly between the two formula-fed groups. In contrast, addition of E2 to milk formula resulted in a significant increase in tubule area/total area relative to those of the cow’s milk formula group (p < 0.05). Addition of genistein resulted in no significant effect relative to feeding cow’s milk by itself (Figure 1). No significant effects of either formula or addition of E2 or genistein to cow’s milk formula were observed.
Effects on testicular Leydig cell number are shown in Figure 2. Formula feeding from either cow’s milk or soy both resulted in reduced Sertoli cell numbers relative to feeding sow milk (p < 0.05) but soy formula and cow’s milk formula groups did not differ significantly (Figure 3).
In contrast, E2 treatment increased Sertoli cell number in the M + E2 group vs. the Milk group (p < 0.05), while addition of genistein had no effect (Figure 3). Consumption of cow’s milk formula reduced germ cell numbers relative to those with soy formula feeding (p < 0.05) and feeding of sow milk, while sow and soy formula groups did not differ significantly. In contrast, germ cell numbers were increased by addition of E2 (p < 0.05) to cow’s milk formula relative to those of the Milk group (Figure 3).

6.5. Effects of Formula Feeding, Genistein, and E2 on Expression of mRNAs Encoding Steroid Biosynthetic Enzymes in Neonatal Male Piglet Testis

Given the suppression of piglet testicular growth and circulating androgenic steroids following E2 supplementation of cow’s milk formula, we examined the expression of testicular steroidogenic enzymes by real-time RT-PCR. These data are shown in Figure 4 and Figure 5. mRNA encoding the steroidogenic acute regulatory protein (StAR), which regulates cholesterol transport within the mitochondria and which is considered a rate-limiting step in steroid synthesis, was suppressed in the M + E2 group relative to that in the Milk and M + G groups (p = 0.099), but its expression was completely unaffected by either formula diet relative to that of piglets fed sow milk (Figure 4).
In contrast, no significant effects of neonatal diet, E2, or genistein supplementation were observed on expression of pregnenolone synthase (CYP11A1) (Figure 4) or steroid 17α-hydroxylase/17–20 lyase (CYP17A1) (Figure 5). In contrast, E2 supplementation of cow’s milk formula resulted in suppression of expression of aromatase (CYP19) mRNA (p < 0.05), which encodes the enzyme converting androgens to estrogens. No effects of neonatal diet or genistein supplementation were observed on CYP19 mRNA expression (Figure 5).

6.6. RNAseq Analysis of Male Piglet Testis

RNAseq analysis was conducted on RNA extracted from the testes of all experimental groups. Principal component analysis on the log scale for RPKM values revealed that all the experimental groups were different. Cow’s milk and soy formula groups did not segregate together but were both very different from the sow milk group. Addition of either E2 or genistein to cow’s milk resulted in different RNA expression relative to that of the Sow, but neither pattern was similar to that of the soy formula group (Figure 6). Primary direct comparisons were made between the Sow group and other treatments. A comparison between groups of genes regulated by more than two-fold showed some overlap between E2 and genistein groups relative to the Sow group, but very little overlap between the E2 and soy formula group (Figure 7).
Only 75 of 2555 regulated genes relative to Sow (2.9%) were in common between the E2-treated and soy-fed groups. There were 896 of 2555 regulated genes (35%) in common between E2- and genistein-supplemented groups. Seventy-five of 94 soy formula-regulated genes relative to Sow (80%) were also regulated by genistein, but these only represented 6% of genistein-regulated genes. The top KEGG pathways for the E2-supplemented group relative to the sow group were oxidative phosphorylation and ribosome, but other pathways such as metabolism of xenobiotics by cytochrome P450 and steroid biosynthesis were also affected (p-value < 0.1) (Table 5). Excel files reporting up- and downregulated testicular genes by E2, GE, and soy vs. sow formula-fed are found in Tables S1–S3, and all the RNAseq data are publicly available via the NCBI SRA database (GSE 1688627). The most upregulated gene in E2 was mysterin (RNF213), which contains an AAA domain associated with ATPase activity and has been recently reported as a regulator of cytoplasmic lipid droplets [40]. Another highly induced mRNA by E2 was the transcription factor FOXL2, which has been shown to be important in female sex determination [41].
E2 downregulated genes included the immune regulators GZMA and LECT2 [42,43,44,45]. The most upregulated gene by soy was SMAD6, which is an inhibitor of BMPs and TGF-b signaling, and in vitro experiments showed that SMAD6 overexpression blocks BMPs-induced osteoblast and chondrocyte differentiation [46]. Of note, men with deletions in the Y-chromosome DAZ cluster experience impaired spermatogenesis [47], and only in the M + G group did we find downregulation of the DAZ homolog (DAZL) [48].
Secondary direct comparisons were made between gene expression profiles in the M + E2, M + G, and soy formula groups vs. the Milk formula-fed group. We identified 103 genes differentially expressed in the M + E2 group vs. Milk. Three of these genes were pig-specific, namely DMBT1L, RPL21Ps, and UNK37, and other novel genes classified as long noncoding RNA ENSSSCG00000003503, ENSSSCG00000046393, and ENSSSCG00000049464. We conducted a pathway analysis with IPA (Figure 8). The networks related to endocrine system disorders and metabolic disorders had the larger number of members related directly or indirectly to our set. Figure 8 depicts the network with colored genes indicating the directionality of the comparison M + E2 vs. Milk. In contrast, there were only three genes highly upregulated in M + G vs. Milk, namely Defensin Beta 115 (DEFB115), Epithelial cell adhesion molecule (EPCAM), and Keratin 19 (KRT19). In rats, defensins are highly expressed in testis [49], suggesting that genistein may potentiate their expression. On the other hand, according to the Human Phenotype Oncology, EPCAM gene has been related to abnormalities of the genital system (HP: 000078). Finally, Keratin 19 gene (also known as CK19) has been shown as a biomarker for a certain type of testis tumors [50]. The comparison of Milk + E2 vs. soy formula groups revealed only four differentially regulated genes C-X-C motif chemokine ligand 13 (CXCL13), the long coding RNA ENSSSCG00000046140, Galactosidase beta 1 like 2 (GLIB1L2), and the mitochondrial gene ENSSSCG00000018062. None of these overlapped with the effects of E2 or genistein supplementation. The former two genes were downregulated by soy formula, whereas the last two were upregulated with respect to milk formula. Interestingly, it has been reported that blocking CXCL13 expression prevents the emergence of aggressive and metastatic prostate cancers treated with hormone ablation therapies [51].

7. Discussion

There is considerable discordance in the toxicological literature with regards to the potential estrogenicity and reproductive toxicity of soy infant formula as the result of the presence of the isoflavone phytoestrogens genistein and daidzein. This is especially true of effects on male reproductive development. Neonatal estrogens are well-known to reduce testis size, suppress steroidogenesis, and decrease testosterone concentrations [21,22,23,24]. Similar effects have been observed in rodents following perinatal treatment with genistein- or soy-containing diets [21,25,26]. Moreover, Sharpe et al. have reported decreased serum testosterone and increased Leydig cell numbers in marmoset monkeys after postnatal feeding of soy formula, suggesting development of compensated Leydig cell failure [11,12]. In contrast, studies in our laboratory have observed no male reproductive toxicity following lifetime feeding of soy protein isolate (SPI), the sole protein source in soy formula, and no pattern of estrogen-regulated gene in prenatal rats fed SPI in direct comparison to the effects of E2 [27,30,34]. Some of these differences may be related to the different effects of soy foods, which are complex mixtures of proteins and over 100 phytochemicals in addition to isoflavones, to differences in developmental window and route of administration, and to species’ differences in isoflavone metabolism [52,53]. Rodents and marmoset monkeys metabolize daidzein to the more potent estrogenic metabolite equol, whereas human infants do not make equol until after weaning [53,54]. The current study utilized the piglet, a more appropriate animal model for male neonatal development, which does not make equol [53], and directly compared the effects of soy and cow’s milk formula feeding with that of sow milk and the effects of cow’s milk formula supplementation with pure genistein levels comparable with those in soy formula or with E2 as a positive control.
E2 supplementation of cow’s milk formula resulted in the expected decrease in testis weight, steroidogenesis, and serum testosterone. In addition, E2 resulted in morphological changes in the neonatal pig testis, including increased tubular area, Sertoli and germ cell numbers, and significantly reduced serum FSH concentrations, and produced expected changes in testis gene expression profiles linked to endocrine/metabolic disorders. In contrast, neither soy formula nor pure genistein supplementation at levels comparable to those seen in soy formula had similar effects. In fact, soy formula significantly decreased tubular area and Sertoli cell numbers relative to those in the Sow group, and soy formula and genistein significantly increased FSH.
It is unclear if the effect on neonatal FSH has physiological significance, but this effect is clearly not associated with the activation of estrogenic signaling. In fact, these effects on testicular cellularity and FSH suggest possible anti-estrogenic actions as a result of competition between low-affinity isoflavones and high-affinity endogenous estrogens at the level of the estrogen receptors or estrogen-independent actions. In this regard, anti-estrogenic effects of SPI have previously been reported in the rat mammary gland [32]. The only similarity between SPI and E2 effects resided in a small increase in germ cell numbers in soy formula-fed vs. cow’s milk formula-fed groups which did not differ significantly from the germ cell numbers observed in sow milk-fed piglets. These results also differ significantly from those described by Sharpe et al. in marmoset monkeys fed soy formula [11,12] with no effects on testis size, serum testosterone, or Leydig cell numbers and no effect on serum LH concentrations. These effects may reflect species’ differences in isoflavone metabolism given the extensive production of equol in monkeys and lack of equol in piglets and human infants prior to weaning [53]. A more detailed characterization of testicular effects of diet, genistein, and E2 using RNAseq revealed clear E2 effects on multiple testicular pathways in the M + E2 group vs. the Milk group, but little overlap between E2 and soy formula-fed groups or between the genistein-supplemented and soy formula-fed groups. These data are consistent with previous studies on genome-wide comparisons of the effects of estradiol and soy protein isolate/soy formula in estrogen-responsive tissues of rat and piglets of different sexes and ages indicating that soy acts as a selective estrogen receptor modulator (SERM) rather than a complete estrogen, only overlapping with a small subset of estrogen-responsive pathways [28,29,30,31,32,33,34]. The strengths of this study are that it compared the effects of physiologically relevant levels of estrogens with that of soy formula in an animal model where commercial infant formula could be fed orally during the developmental window, corresponding to formula feeding in infants and where the pattern of isoflavone metabolites was the same as that found in infants. The limitations are that the number of animals per group was small, that extrapolation between animal models and humans must always be treated with caution, and that it is very difficult to definitively prove a negative.
In this regard, it is important to note that our data are also consistent with those of two recent human longitudinal clinical studies of breast-fed and formula-fed male infants. Clinical studies of soy formula-fed infants utilizing ultrasound analysis of breast bud, testis, and prostate development at ages of one and five years in the “Beginnings Study” at the Arkansas Children’s Nutrition Center found no significant effects of soy formula vs. cow’s milk formula in male infants [6,7]. Likewise, a more recent study at the NIEHS also found no significant effects of soy formula vs. cow’s milk formula on breast bud development or urethral epithelial maturation index or estrogen or FSH hormone trajectories in soy formula-fed boys during the first seven months of life [10]. These clinical studies differed in results relating to formula-fed girls. In the “Beginnings Study”, no significant effects of soy formula feeding were observed on breast bud development and uterine or ovarian development using ultrasound analysis of girls fed soy vs. cow’s milk formula at ages of one and five years [6,7]. In contrast, the NIEHS study observed a significantly higher vaginal cell mitotic index and slower decreases in uterine volume in soy formula-fed vs. those in cow’s formula-fed girls [10]. However, these investigators also observed no effects of soy on breast bud dimeter or hormone concentrations in infant girls. It is unclear why the results of these clinical studies differ. In particular, it is unclear why, if soy indeed acts as an estrogen, the latter study did not also observe significant effects on female breast bud development or on male parameters, both of which have been shown to be highly sensitive to estrogens. An additional epigenetic analysis of the vaginal cells from the formula-fed girls in the NIEHS study has been conducted [9]. This study observed changes in methylation of CpG islands near the transcription start site of one estrogen-regulated gene proline rich 5-like (PRR5L) in cells from soy formula-fed girls vs. that of cow’s milk formula group. Such a result is not inconsistent with the hypothesis that soy acts as a SERM interacting with a small subset of estrogen-regulated genes.

8. Conclusions

These data strongly suggest that neither soy formula nor genistein were estrogenic in the neonatal porcine testis and that soy formula feeding during neonatal development did not significantly affect male reproductive development in this model. Moreover, these data are consistent with the conclusions of three recent comprehensive reviews of the clinical literature on the effects of dietary isoflavones and soy products [17,55,56] that there is little to no evidence that genistein or other phytoestrogens in soy formula consumed during the postnatal period significantly alter male reproductive development. Collectively, these studies are in agreement with the support of experts in the U.S. for the use of soy formula as a safe and cost-effective alternative to cow’s milk formula [55].

Supplementary Materials

The following supporting information can be downloaded at: www.mdpi.com/article/10.3390/nu14051126/s1, Tables S1–S3: Up- and downregulated testicular genes by M + E2, M + G and soy vs. sow formula-fed.

Author Contributions

All authors reviewed the manuscript and performed data analysis. M.J.J.R., K.S., I.M. and T.M.B. participated in research design. I.M., N.S., K.E.M. and M.L.B. conducted experiments and analysis. L.H. conducted blinded pathological analysis. H.G.-A., C.C. and H.D. conducted bioinformatics analysis. M.J.J.R. and H.G.-A. wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This article was funded in part by United States Department of Agriculture Service Project 6026-51000-010-05S, which had no role in study design, data collection and analysis, interpretation, writing of the manuscript, or in the decision to submit it for publication.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Animal Care and Use Committee at the University of Arkansas for Medical Sciences (#3122, approved 8 March 2010).

Informed Consent Statement

Not applicable.

Data Availability Statement

All the RNAseq data are publicly available via the NCBI SRA database (GSE 1688627).

Acknowledgments

We wish to thank Matthew Ferguson, Tammy Dalrari, Rachel Nicolson, Mark Robinette, and Jamie Vantrese for technical assistance in the preparation of this manuscript. The license to Ingenuity Pathways was kindly sponsored by Shuk-Mei Ho, Vice Chancellor for Research and Innovation at UAMS.

Conflicts of Interest

The authors declare no conflict of interest with respect to the funding of this research, authorship, and/or publication of this article.

References

  1. American Academy of Pediatrics. Promotion of breast feeding: Recommendations of the Councils of the Society for Pediatric Research (SPR) and American Pediatric Society (APS), and of the American Academy of Pediatrics (AAP). Pediatr. Res. 1982, 16, 264–265. [Google Scholar] [CrossRef] [Green Version]
  2. Sayres, S.; Visentin, L. Breastfeeding: Uncovering barriers and offering solutions. Curr. Opin. Pediatr. 2018, 30, 591–596. [Google Scholar] [CrossRef]
  3. Martin, C.R.; Ling, P.R.; Blackburn, G.L. Review of Infant Feeding: Key Features of Breast Milk and Infant Formula. Nutrients 2016, 8, 279. [Google Scholar] [CrossRef] [Green Version]
  4. Merritt, R.J.; Jenks, B.H. Safety of soy-based infant formulas containing isoflavones: The clinical evidence. J. Nutr. 2004, 134, 1220S–1224S. [Google Scholar] [CrossRef]
  5. Badger, T.M.; Gilchrist, J.M.; Pivik, R.T.; Andres, A.; Shankar, K.; Chen, J.R.; Ronis, M.J. The health implications of soy infant formula. Am. J. Clin. Nutr. 2009, 89, 1668S–1672S. [Google Scholar] [CrossRef]
  6. Gilchrist, J.M.; Moore, M.B.; Andres, A.; Estroff, J.A.; Badger, T.M. Ultrasonographic patterns of reproductive organs in infants fed soy formula: Comparisons to infants fed breast milk and milk formula. J. Pediatr. 2010, 156, 215–220. [Google Scholar] [CrossRef]
  7. Andres, A.; Moore, M.B.; Linam, L.E.; Casey, P.H.; Cleves, M.A.; Badger, T.M. Compared with feeding infants breast milk or cow-milk formula, soy formula feeding does not affect subsequent reproductive organ size at 5 years of age. J. Nutr. 2015, 145, 871–875. [Google Scholar] [CrossRef] [Green Version]
  8. Vandenplas, Y.; De Greef, E.; Devreker, T.; Hauser, B. Soy infant formula: Is it that bad? Acta Paediatr. 2011, 100, 162–166. [Google Scholar] [CrossRef]
  9. Harlid, S.; Adgent, M.; Jefferson, W.N.; Panduri, V.; Umbach, D.M.; Xu, Z.; Stallings, V.A.; Williams, C.J.; Rogan, W.J.; Taylor, J.A. Soy Formula and Epigenetic Modifications: Analysis of Vaginal Epithelial Cells from Infant Girls in the IFED Study. Environ. Health Perspect. 2017, 125, 447–452. [Google Scholar] [CrossRef] [Green Version]
  10. Adgent, M.A.; Umbach, D.M.; Zemel, B.S.; Kelly, A.; Schall, J.I.; Ford, E.G.; James, K.; Darge, K.; Botelho, J.C.; Vesper, H.W.; et al. A Longitudinal Study of Estrogen-Responsive Tissues and Hormone Concentrations in Infants Fed Soy Formula. J. Clin. Endocrinol. Metab. 2018, 103, 1899–1909. [Google Scholar] [CrossRef]
  11. Sharpe, R.M.; Martin, B.; Morris, K.; Greig, I.; McKinnell, C.; McNeilly, A.S.; Walker, M. Infant feeding with soy formula milk: Effects on the testis and on blood testosterone levels in marmoset monkeys during the period of neonatal testicular activity. Hum. Reprod. 2002, 17, 1692–1703. [Google Scholar] [CrossRef] [Green Version]
  12. Tan, K.A.; Walker, M.; Morris, K.; Greig, I.; Mason, J.I.; Sharpe, R.M. Infant feeding with soy formula milk: Effects on puberty progression, reproductive function and testicular cell numbers in marmoset monkeys in adulthood. Hum. Reprod. 2006, 21, 896–904. [Google Scholar] [CrossRef] [Green Version]
  13. McCarver, G.; Bhatia, J.; Chambers, C.; Clarke, R.; Etzel, R.; Foster, W.; Hoyer, P.; Leeder, J.S.; Peters, J.M.; Rissman, E.; et al. NTP-CERHR expert panel report on the developmental toxicity of soy infant formula. Birth Defects Res. B Dev. Reprod. Toxicol. 2011, 92, 421–468. [Google Scholar] [CrossRef]
  14. Messina, M.J.; Persky, V.; Setchell, K.D.; Barnes, S. Soy intake and cancer risk: A review of the in vitro and in vivo data. Nutr. Cancer 1994, 21, 113–131. [Google Scholar] [CrossRef]
  15. Setchell, K.D.; Zimmer-Nechemias, L.; Cai, J.; Heubi, J.E. Exposure of infants to phyto-oestrogens from soy-based infant formula. Lancet 1997, 350, 23–27. [Google Scholar] [CrossRef]
  16. Vitale, D.C.; Piazza, C.; Melilli, B.; Drago, F.; Salomone, S. Isoflavones: Estrogenic activity, biological effect and bioavailability. Eur. J. Drug Metab. Pharmacokinet. 2013, 38, 15–25. [Google Scholar] [CrossRef]
  17. Messina, M.; Barnes, S.; Cassidy, A.; Duncan, A.; Kurzur, M.; Meija, S.B.; Ronis, M.; Rowland, I.; Sievenpiper, J. Neither Soyfoods nor isoflavones warrant classification as endocrine disruptors: A review of clinical and observational data. Crit. Rev. Food Sci. Nutr. 2021, 27, 1–57. [Google Scholar] [CrossRef]
  18. Kuiper, G.G.; Lemmen, J.G.; Carlsson, B.; Corton, J.C.; Safe, S.H.; van der Saag, P.T.; van der Burg, B.; Gustafsson, J.A. Interaction of estrogenic chemicals and phytoestrogens with estrogen receptor beta. Endocrinology 1998, 139, 4252–4263. [Google Scholar] [CrossRef]
  19. Morito, K.; Aomori, T.; Hirose, T.; Kinjo, J.; Hasegawa, J.; Ogawa, S.; Inoue, S.; Muramatsu, M.; Masamune, Y. Interaction of phytoestrogens with estrogen receptors alpha and beta (II). Biol. Pharm. Bull. 2002, 25, 48–52. [Google Scholar] [CrossRef] [Green Version]
  20. Ye, H.; Shaw, I.C. Food flavonoid ligand structure/estrogen receptor-α affinity relationships—Toxicity or food functionality? Food Chem. Toxicol. 2019, 129, 328–336. [Google Scholar] [CrossRef]
  21. Cederroth, C.R.; Zimmermann, C.; Nef, S. Soy, phytoestrogens and their impact on reproductive health. Mol. Cell. Endocrinol. 2012, 355, 192–200. [Google Scholar] [CrossRef] [PubMed]
  22. Akingbemi, B.T. Estrogen regulation of testicular function. Reprod. Biol. Endocrinol. 2005, 3, 51. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Nielsen, M.; Bjornsdottir, S.; Hoyer, P.E.; Byskov, A.G. Ontogeny of oestrogen receptor alpha in gonads and sex ducts of fetal and newborn mice. J. Reprod. Fertil. 2001, 118, 195–204. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Atanassova, N.; McKinnell, C.; Walker, M.; Turner, K.J.; Fisher, J.S.; Morley, M.; Millar, M.R.; Groome, N.P.; Sharpe, R.M. Permanent effects of neonatal estrogen exposure in rats on reproductive hormone levels, Sertoli cell number, and the efficiency of spermatogenesis in adulthood. Endocrinology 1999, 140, 5364–5373. [Google Scholar] [CrossRef]
  25. Wisniewski, A.B.; Klein, S.L.; Lakshmanan, Y.; Gearhart, J.P. Exposure to genistein during gestation and lactation demasculinizes the reproductive system in rats. J. Urol. 2003, 169, 1582–1586. [Google Scholar] [CrossRef]
  26. Meena, R.; Supriya, C.; Pratap Reddy, K.; Sreenivasula Reddy, P. Altered spermatogenesis, steroidogenesis and suppressed fertility in adult male rats exposed to genistein, a non-steroidal phytoestrogen during embryonic development. Food Chem. Toxicol. 2017, 99, 70–77. [Google Scholar] [CrossRef]
  27. Badger, T.M.; Ronis, M.J.; Hakkak, R. Developmental effects and health aspects of soy protein isolate, casein, and whey in male and female rats. Int. J. Toxicol. 2001, 20, 165–174. [Google Scholar] [CrossRef]
  28. Singhal, R.; Shankar, K.; Badger, T.M.; Ronis, M.J. Hepatic gene expression following consumption of soy protein isolate in female Sprague Dawley rats differs from that produced by 17β-estradiol treatments. J. Endocrinol. 2009, 202, 141–152. [Google Scholar] [CrossRef]
  29. Zhang, J.; Lazarenko, O.P.; Wu, X.; Tong, Y.; Blackburn, M.L.; Gomez-Acevedo, H.; Shankar, K.; Badger, T.M.; Ronis, M.J.; Chen, J.R. Differential effects of short-term feeding of a soy protein isolate diet and estrogen treatment on bone in the pre-pubertal rat. PLoS ONE 2012, 7, e35736. [Google Scholar] [CrossRef]
  30. Miousse, I.R.; Sharma, N.; Blackburn, M.; Vantrease, J.; Gomez-Acevedo, H.; Hennings, L.; Shankar, K.; Cleves, M.A.; Badger, T.M.; Ronis, M.J. Feeding soy protein isolate and treatment with estradiol have different effects on mammary gland morphology and gene expression in weanling male and female rats. Physiol. Genom. 2013, 45, 1072–1083. [Google Scholar] [CrossRef] [Green Version]
  31. Ronis, M.J.; Shankar, K.; Gomez-Acevedo, H.; Hennings, L.; Singhal, R.; Blackburn, M.L.; Badger, T.M. Mammary gland morphology and gene expression differ in female rats treated with 17β-estradiol or fed soy protein isolate. Endocrinology 2012, 153, 6021–6032. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Ronis, M.J.; Gomez-Acevedo, H.; Blackburn, M.L.; Cleves, M.A.; Singhal, R.; Badger, T.M. Uterine responses to feeding soy protein isolate and treatment with 17β-estradiol differ in ovariectomized female rats. Toxicol. Appl. Pharmacol. 2016, 297, 68–80. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Ronis, M.J.; Gomez-Acevedo, H.; Shankar, K.; Sharma, N.; Blackburn, M.; Singhal, R.; Mercer, K.E.; Badger, T.M. Soy protein isolate feeding does not result in reproductive toxicity in the pre-pubertal rat testis. Exp. Biol. Med. (Maywood) 2018, 243, 695–707. [Google Scholar] [CrossRef] [PubMed]
  34. Mercer, K.E.; Bhattacharyya, S.; Sharma, N.; Chaudhury, M.; Lin, H.; Yeruva, L.; Ronis, M.J. Infant Formula Feeding Changes the Proliferative Status in Piglet Neonatal Mammary Glands Independently of Estrogen Signaling. J. Nutr. 2020, 150, 730–738. [Google Scholar] [CrossRef]
  35. Chen, J.R.; Lazarenko, O.P.; Blackburn, M.; Badeaux, J.V.; Badger, T.M.; Ronis, M.J. Infant formula promotes bone growth in neonatal piglets by enhancing osteoblastogenesis through bone morphogenic protein signaling. J. Nutr. 2009, 139, 1839–1847. [Google Scholar] [CrossRef] [Green Version]
  36. Lanning, L.L.; Creasy, D.M.; Chapin, R.E.; Mann, P.C.; Barlow, N.J.; Regan, K.S.; Goodman, D.S. Recommended approaches for the evaluation of testicular and epididymal toxicity. Toxicol. Pathol. 2002, 30, 507–520. [Google Scholar] [CrossRef] [Green Version]
  37. Ronis, M.J.; Chen, Y.; Shankar, K.; Gomez-Acevedo, H.; Cleves, M.A.; Badeaux, J.; Blackburn, M.L.; Badger, T.M. Formula feeding alters hepatic gene expression signature, iron and cholesterol homeostasis in the neonatal pig. Physiol. Genom. 2011, 23, 1281–1293. [Google Scholar] [CrossRef] [Green Version]
  38. Chen, C.; Perry, T.L.l.; Chitko-McKown, C.G.; Smith, A.D.; Cheung, L.; Beshah, E.; Urban, J.F.; Dawson, H.D. The regulatory actions of retinoic acid on M2 polarization of porcine macrophages. Dev. Comp. Immunol. 2019, 98, 20–33. [Google Scholar] [CrossRef]
  39. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [Green Version]
  40. Sugihara, M.; Morito, D.; Ainuki, S.; Hirano, Y.; Ogino, K.; Kitamura, A.; Hirata, H.; Nagata, K. The AAA+ ATPase/ubiquitin ligase mysterin stabilizes cytoplasmic lipid droplets. J. Cell Biol. 2019, 218, 949–960. [Google Scholar] [CrossRef] [Green Version]
  41. Huang, S.; Ye, L.; Chen, H. Sex determination and maintenance: The role of DMRT1 and FOXL2. Asian J. Androl. 2017, 19, 619–624. [Google Scholar] [CrossRef] [PubMed]
  42. Phan, N.N.; Wang, C.Y.; Chen, C.F.; Sun, Z.; Lai, M.D.; Lin, Y.C. Voltage-gated calcium channels: Novel targets for cancer therapy. Oncol. Lett. 2017, 14, 2059–2074. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Kumar, S.; Kang, H.; Park, E.; Park, H.S.; Lee, K. The expression of CKLFSF2B is regulated by GATA1 and CREB in the Leydig cells, which modulates testicular steroidogenesis. Biochim. Biophys. Acta Gene Regul. Mech. 2018, 1861, 1063–1075. [Google Scholar] [CrossRef] [PubMed]
  44. Dudem, S.; Large, R.J.; Kulkarni, S.; McClafferty, H.; Tikhonova, I.G.; Sergeant, G.P.; Thornbury, K.D.; Shipston, M.J.; Perrino, B.A.; Hollywood, M.A. LINGO1 is a regulatory subunit of large conductance, Ca 2+-activated potassium channels. Proc. Natl. Acad. Sci. USA 2020, 117, 2194–2200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Dutta, D.; Park, I.; Guililat, H.; Sang, S.; Talapatra, A.; Singhal, B.; Mills, N.C. Testosterone regulates granzyme K expression in rat testes. Endocr. Regul. 2017, 51, 193–204. [Google Scholar] [CrossRef] [Green Version]
  46. Wu, M.; Chen, G.; Li, Y.-P. TGF-β and BMP signaling in osteoblast, skeletal development, and bone formation, homeostasis and disease. Bone Res. 2016, 4, 16009. [Google Scholar] [CrossRef]
  47. Lu, X.J.; Chen, J.; Yu, C.H.; Shi, Y.H.; He, Y.Q.; Zhang, R.C.; Huang, Z.A.; Lv, J.N.; Zhang, S.; Xu, L. LECT2 protects mice against bacterial sepsis by activating macrophages via the CD209a receptor. J. Exp. Med. 2013, 210, 5–13. [Google Scholar] [CrossRef] [Green Version]
  48. Moore, F.L.; Jaruzelska, J.; Fox, M.S.; Urano, J.; Firpo, M.T.; Ture, P.J.; Dorfman, D.M.; Reijo Pera, R.A. Human Pumilio-2 is expressed in embryonic stem cells and germ cells and interacts with DAZ (Deleted in Azoospermia) and DAZ-Like proteins. Proc. Natl. Acad. Sci. USA 2002, 100, 538–543. [Google Scholar] [CrossRef] [Green Version]
  49. Patil, A.A.; Cai, Y.; Sang, Y.; Blecha, F.; Zhang, G. Cross-species analysis of the mammalian beta-defensin gene family: Presence of syntenic gene clusters and preferential expression in the male reproductive tract. Physiol. Genom. 2005, 23, 5–17. [Google Scholar] [CrossRef]
  50. Bremmer, F.; Strobel, P.; Jarry, H.; Strecker, J.; Gaisa, N.; Strauss, A.; Schweyer, S.; Radzun, H.J.; Behnes, C.L. CK19 is a sensitive marker for yolk sac tumours of the testis. Diagn. Pathol. 2015, 10, 7. [Google Scholar] [CrossRef] [Green Version]
  51. Ammirante, M.; Shalapour, S.; Kang, Y.; Jamieson, C.A.; Karin, M. Tissue injury and hypoxia promote malignant progression of prostate cancer by inducing CXCL13 expression in tumor myofibroblasts. Proc. Natl. Acad. Sci. USA 2014, 111, 14776–14781. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Fang, N.; Yu, S.; Badger, T.M. Comprehensive phytochemical profile of soy protein isolate. J. Agric. Food Chem. 2004, 52, 4012–4020. [Google Scholar] [CrossRef] [PubMed]
  53. Gu, L.; House, S.E.; Prior, R.L.; Fang, N.; Ronis, M.J.; Clarkson, T.B.; Wilson, M.E.; Badger, T.M. Metabolic phenotype of isoflavones differ among female rats, pigs, monkeys, and women. J. Nutr. 2006, 136, 1215–1221. [Google Scholar] [CrossRef] [PubMed]
  54. Brown, N.M.; Galandi, S.L.; Summer, S.S.; Zhao, X.; Heubi, J.E.; King, E.C.; Setchell, K.D.R. S-(-) equol production is developmentally regulated and related to early diet composition. Nutr. Res. 2014, 34, 401–409. [Google Scholar] [CrossRef] [Green Version]
  55. Suen, A.A.; Kenan, A.C.; Williams, C.J. Developmental exposure to phytoestrogens in soy: New findings and clinical implications. Biochem. Pharmacol. 2022, 195, 114848. [Google Scholar] [CrossRef]
  56. Dominguez-Lopez, I.; Yago-Aragon, M.; Salas-Huetos, A.; Tresserra-Rimbau, A.; Hurtado-Barroso, S. Effects of dietary phytoestrogens on hormones throughout a human lifespan: A review. Nutrients 2020, 12, 2456. [Google Scholar] [CrossRef]
Nutrients 14 01126 g001 550
Figure 1. Effects of infant formula or estradiol (E2) and genistein on piglet testicular tubule area/total area. Data are mean ± SEM based on assessment of 4 fields from 3 sections of each animal (n = 6/group). a < b statistically different.
Figure 1. Effects of infant formula or estradiol (E2) and genistein on piglet testicular tubule area/total area. Data are mean ± SEM based on assessment of 4 fields from 3 sections of each animal (n = 6/group). a < b statistically different.
Nutrients 14 01126 g001
Nutrients 14 01126 g002 550
Figure 2. Effects of infant formula or estradiol (E2) and genistein on piglet testicular Leydig cell numbers. Data are mean ± SEM based on assessment of 4 fields from 3 sections of each animal (n = 6/group). There were no significant differences between groups.
Figure 2. Effects of infant formula or estradiol (E2) and genistein on piglet testicular Leydig cell numbers. Data are mean ± SEM based on assessment of 4 fields from 3 sections of each animal (n = 6/group). There were no significant differences between groups.
Nutrients 14 01126 g002
Nutrients 14 01126 g003 550
Figure 3. Effects of infant formula or estradiol (E2) and genistein on piglet testicular Sertoli cell and germ cell numbers. Data are mean ± SEM based on assessment of 4 fields from 3 sections of each animal (n = 6/group). a < b statistically different.
Figure 3. Effects of infant formula or estradiol (E2) and genistein on piglet testicular Sertoli cell and germ cell numbers. Data are mean ± SEM based on assessment of 4 fields from 3 sections of each animal (n = 6/group). a < b statistically different.
Nutrients 14 01126 g003
Nutrients 14 01126 g004 550
Figure 4. Effects of infant formula or estradiol (E2) and genistein on piglet testicular steroidogenic enzyme mRNA expression I: steroid-activated receptor (StAR) and cytochrome P450 CYP11A1. Data are mRNA/18S relative to sow milk (Sow) (formulas) or relative to cow’s milk formula (MILK) (E2/genistein) (n = 6/group). One-way ANOVA followed by Student–Newman–Keuls post hoc analysis only showed a trend (p = 0.99) for reduction of StAR mRNA expression by MILK + E2 relative to those of MILK or MILK + genistein. Means with the same superscript letter a do not differ statistically.
Figure 4. Effects of infant formula or estradiol (E2) and genistein on piglet testicular steroidogenic enzyme mRNA expression I: steroid-activated receptor (StAR) and cytochrome P450 CYP11A1. Data are mRNA/18S relative to sow milk (Sow) (formulas) or relative to cow’s milk formula (MILK) (E2/genistein) (n = 6/group). One-way ANOVA followed by Student–Newman–Keuls post hoc analysis only showed a trend (p = 0.99) for reduction of StAR mRNA expression by MILK + E2 relative to those of MILK or MILK + genistein. Means with the same superscript letter a do not differ statistically.
Nutrients 14 01126 g004
Nutrients 14 01126 g005 550
Figure 5. Effects of infant formula or estradiol (E2) and genistein on piglet testicular steroidogenic enzyme mRNA expression II: cytochrome P450 17A1 (CYP17A1) and cytochrome P450 19 (aromatase, CYP19). Data are mRNA/18S relative to sow milk (Sow) (formulas) or relative to cow’s milk formula (MILK) (E2/genistein) (n = 6/group). One-way ANOVA followed by Student–Newman–Keuls post hoc analysis only showed a trend (p = 0.062) for increased CYP19 mRNA expression by soy formula relative to that of sow milk or cow’s milk but b < a, p < 0.05 for E2 in cow’s milk relative to that in cow’s milk.
Figure 5. Effects of infant formula or estradiol (E2) and genistein on piglet testicular steroidogenic enzyme mRNA expression II: cytochrome P450 17A1 (CYP17A1) and cytochrome P450 19 (aromatase, CYP19). Data are mRNA/18S relative to sow milk (Sow) (formulas) or relative to cow’s milk formula (MILK) (E2/genistein) (n = 6/group). One-way ANOVA followed by Student–Newman–Keuls post hoc analysis only showed a trend (p = 0.062) for increased CYP19 mRNA expression by soy formula relative to that of sow milk or cow’s milk but b < a, p < 0.05 for E2 in cow’s milk relative to that in cow’s milk.
Nutrients 14 01126 g005
Nutrients 14 01126 g006 550
Figure 6. Principal component analysis(PCA) of pig testicular mRNAseq data demonstrating lack of clustering between soy formula-, genistein-, or E2-treated groups.
Figure 6. Principal component analysis(PCA) of pig testicular mRNAseq data demonstrating lack of clustering between soy formula-, genistein-, or E2-treated groups.
Nutrients 14 01126 g006
Nutrients 14 01126 g007 550
Figure 7. Pie charts illustrating overlap in gene expression profiles between E2 and genistein groups; E2 and soy formula groups and genistein and soy formula groups.
Figure 7. Pie charts illustrating overlap in gene expression profiles between E2 and genistein groups; E2 and soy formula groups and genistein and soy formula groups.
Nutrients 14 01126 g007
Nutrients 14 01126 g008 550
Figure 8. Ingenuity Pathways Assist analysis of M + E2 effects on gene expression and biochemical pathways relative to those of the Milk group. Estrogen induced and inhibited testicular gene expression. Gene names abbreviated according to the NCBI database (GeneCards.org).
Figure 8. Ingenuity Pathways Assist analysis of M + E2 effects on gene expression and biochemical pathways relative to those of the Milk group. Estrogen induced and inhibited testicular gene expression. Gene names abbreviated according to the NCBI database (GeneCards.org).
Nutrients 14 01126 g008
Table
Table 1. Serum isoflavone concentrations in neonatal piglets fed sow milk or formula.
Table 1. Serum isoflavone concentrations in neonatal piglets fed sow milk or formula.
DietGenistein (pg/mL)Daidzein (pg/mL)Equal
(pg/mL)		0-DMA (pg/mL)Glycitein (pg/mL)
Sow Milk56 ± 8 a256 ± 7 a5 ± 10 a139 ± 87 a
Milk Formula23 ± 9 a2 ± 1 a4 ± 10 a0 a
Soy Formula1712 ± 212 b1673 ± 374 b3 ± 12995 ± 714 b4567 ± 867 b
Milk + Genistein1898 ± 515 b7 ± 1 a9 ± 30 a379 ± 241 a
Milk + E2142 ± 1 a3 ± 2 a5 ± 12 ± 1 a9 ± 1 a
Data are mean ± SEM (n = 6). Means with different letters are significantly different p < 0.05, a < b by one-way ANOVA followed by Student–Newman–Keuls post hoc analysis.
Table
Table 2. Body and sex organ weights of neonatal pigs fed sow milk or formula.
Table 2. Body and sex organ weights of neonatal pigs fed sow milk or formula.
Body WeightTestis Weight%BWProstate Weight%BW
(kg)(g) (g)
Sow Milk7.30 ± 0.186.64 ± 0.76 a0.0090 ± 0.00090.38 ± 0.020.0005 ± 0.0001
Milk Formula8.33 ± 0.277.93 ± 0.57 a,*0.0096 ± 0.0010 *0.37 ± 0.030.0005 ± 0.0001
Soy Formula8.52 ± 0.3810.54 ± 1.37 b0.0122 ± 0.00130.31 ± 0.040.0004 ± 0.0001
Milk + Genistein7.88 ± 0.429.46 ± 1.50 a,*0.0122 ± 0.0014 *0.36 ± 0.050.0005 ± 0.0001
Milk + E28.27 ± 0.474.31 ± 0.5 a,#0.0052 ± 0.0006 #0.44 ± 0.050.0005 ± 0.0001
Data are mean ± SEM (n = 6). Means with different letters are significantly different from that of the sow milk group p < 0.05, a < b; means with different symbols are significantly different from that of the milk formula group p < 0.05, # < * by one-way ANOVA followed by Student–Newman–Keuls post hoc analysis.
Table
Table 3. Gonadotropin–testis axis hormones in serum of male piglets fed sow milk or formula.
Table 3. Gonadotropin–testis axis hormones in serum of male piglets fed sow milk or formula.
DietFSH (mIU/mL)LH (mIU/mL)Androsterone (ng/mL)Testosterone (ng/mL)DHT (ng/mL)
Sow Milk0.10 ± 0.05 a0.007 ± 0.00776 ± 8 b0.52 ± 0.10 b1.65 ± 0.17
Milk Formula0.13 ± 0.04 a,#0.047 ± 0.0242 ± 6 a,#0.35 ± 0.02 b,*1.33 ± 0.17
Soy Formula0.52 ± 0.04 b0.001 ± 0.00186 ± 5 b0.50 ± 0.10 b1.46 ± 0.14
Milk + Genistein0.34 ± 0.04 b,*0.052 ± 0.00364 ± 8 b,*0.42 ± 0.08 b,*1.24 ± 0.20
Milk + E20.07 ± 0.04 a,#0.040 ± 0.0328 ± 6 a,#0.09 ± 0.02 a,#0.88 ± 0.26
ANOVA results for group comparison and post hoc analysis. Data are mean ± SEM (n = 6), and means with differing letters differ significantly from that of the sow milk group p < 0.05, a < b; means with different symbols within the milk groups differ from those of the milk formula group p < 0.05, # < *.
Table
Table 4. Serum levels of non-androgenic steroid hormones in male piglets fed sow milk or formula.
Table 4. Serum levels of non-androgenic steroid hormones in male piglets fed sow milk or formula.
DietDHEA (ng/mL)DHEA-S (ng/mL)Progesterone (ng/mL)Estradiol (pg/mL)Estrone (pg/mL)
Sow Milk3.9 ± 0.229.7 ± 1.1 a4.20 ± 0.25 b115 ± 22 a354 ± 42 a
Milk Formula4.5 ± 0.230.5 ± 1.6 a0.59 ± 0.05 a78 ± 18 a,#240 ± 23 a,#
Soy Formula4.2 ± 0.236.7 ± 1.2 b1.12 ± 0.23 a113 ± 12 a335 ± 31 a
Milk + Genistein4.2 ± 0.432.1 ± 1.0 a,b0.54 ± 0.26 a108 ± 12 a284 ± 21 a,#
Milk + E24.5 ± 0.432.7 ± 2.2 a,b0.85 ± 0.21 a1376 ± 418 b,*1615 ± 321 b,*
Data are mean ± SEM (n = 6). Means with differing letters are significantly different from the sow milk group p < 0.05, a < b. Means with different symbols are significantly different from the milk formula group within milk-fed groups p < 0.05, # < * by one-way ANOVA followed by Student–Newman–Keuls post hoc analysis.
Table
Table 5. Top KEGG pathways for Milk + E2 vs. Sow groups.
Table 5. Top KEGG pathways for Milk + E2 vs. Sow groups.
PATHWAYp-ValueFDR
ssc05204: Chemical carcinogenesis0.0270657370.481770112
ssc04146: Peroxisome0.0308648230.515056728
ssc00980: Metabolism of xenobiotics by cytochrome P4500.0366793940.576082249
ssc01212: Fatty acid metabolism0.0516892750.766724247
ssc00100: Steroid biosynthesis0.0714731270.970909091
ssc05205: Proteoglycans in cancer0.0923737270.970909091
ssc03018: RNA degradation0.0972393740.970909091
ssc04260: Cardiac muscle contraction0.0972393740.970909091
ssc00650: Butanoate metabolism0.09876392500.970909091
FDR: false discovery rate.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ronis, M.J.J.; Gomez-Acevedo, H.; Shankar, K.; Hennings, L.; Sharma, N.; Blackburn, M.L.; Miousse, I.; Dawson, H.; Chen, C.; Mercer, K.E.; et al. Soy Formula Is Not Estrogenic and Does Not Result in Reproductive Toxicity in Male Piglets: Results from a Controlled Feeding Study. Nutrients 2022, 14, 1126. https://0-doi-org.brum.beds.ac.uk/10.3390/nu14051126

AMA Style

Ronis MJJ, Gomez-Acevedo H, Shankar K, Hennings L, Sharma N, Blackburn ML, Miousse I, Dawson H, Chen C, Mercer KE, et al. Soy Formula Is Not Estrogenic and Does Not Result in Reproductive Toxicity in Male Piglets: Results from a Controlled Feeding Study. Nutrients. 2022; 14(5):1126. https://0-doi-org.brum.beds.ac.uk/10.3390/nu14051126

Chicago/Turabian Style

Ronis, Martin J. J., Horacio Gomez-Acevedo, Kartik Shankar, Leah Hennings, Neha Sharma, Michael L. Blackburn, Isabelle Miousse, Harry Dawson, Celine Chen, Kelly E. Mercer, and et al. 2022. "Soy Formula Is Not Estrogenic and Does Not Result in Reproductive Toxicity in Male Piglets: Results from a Controlled Feeding Study" Nutrients 14, no. 5: 1126. https://0-doi-org.brum.beds.ac.uk/10.3390/nu14051126

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop