The mutational landscape of normal human endometrial epithelium.

All normal somatic cells are thought to acquire mutations, but understanding of the rates, patterns, causes and consequences of somatic mutations in normal cells is limited. The uterine endometrium adopts multiple physiological states over a lifetime and is lined by a gland-forming epithelium1,2. Here, using whole-genome sequencing, we show that normal human endometrial glands are clonal cell populations with total mutation burdens that increase at about 29 base substitutions per year and that are many-fold lower than those of endometrial cancers. Normal endometrial glands frequently carry 'driver' mutations in cancer genes, the burden of which increases with age and decreases with parity. Cell clones with drivers often originate during the first decades of life and subsequently progressively colonize the epithelial lining of the endometrium. Our results show that mutational landscapes differ markedly between normal tissues-perhaps shaped by differences in their structure and physiology-and indicate that the procession of neoplastic change that leads to endometrial cancer is initiated early in life.

Genetic evidence suggests a sense of family, parity and conquest in the Xiongnu Iron Age nomads of Mongolia.

In an effort to characterize the people who composed the groups known as the Xiongnu, nuclear and whole mitochondrial DNA data were generated from the skeletal remains of 52 individuals excavated from the Tamir Ulaan Khoshuu (TUK) cemetery in Central Mongolia. This burial site, attributed to the Xiongnu period, was used from the first century BC to the first century AD. Kinship analyses were conducted using autosomal and Y-chromosomal DNA markers along with complete sequences of the mitochondrial genome. These analyses suggested close kin relationships between many individuals. Nineteen such individuals composed a large family spanning five generations. Within this family, we determined that a woman was of especially high status; this is a novel insight into the structure and hierarchy of societies from the Xiongnu period. Moreover, our findings confirmed that the Xiongnu had a strongly admixed mitochondrial and Y-chromosome gene pools and revealed a significant western component in the Xiongnu group studied. Using a fine-scale approach (haplotype instead of haplogroup-level information), we propose Scytho-Siberians as ancestors of the Xiongnu and Huns as their descendants.

Review of epidemiological factors (other than maternal age) that determine the prevalence of common autosomal trisomies.

The birth prevalence of each common autosomal trisomy (21, 18 and 13) increases with advancing maternal age and this is the most important epidemiological risk factor. Prevalence during pregnancy is also dependent on gestational age. Other factors claimed to influence prevalence include paternal age, ethnicity, family history, premature reproductive aging, parity, twinning, smoking, environmental exposures, maternal medical conditions, and predispositions. We review the evidence for these associations since they may provide insights into causal mechanisms. When investigating potential co-factors it is important to adequately allow for maternal age and minimize its confounding contribution. This is well illustrated by reports of an inverse paternal age effect where there is strong correlation between parental ages. Gestational age at diagnosis, availability of prenatal screening, diagnostic testing, and elective termination of affected pregnancies and healthcare disparities also confound the studies on ethnicity, medical conditions, and predispositions or environmental factors. Data from twin zygosity studies demonstrate the importance of differences in fetal viability for affected pregnancies. We conclude that existing epidemiological evidence for most of the co-factors discussed should currently be considered tenuous; history of Down syndrome, albeit biased, may be an exception. The co-factors may yet provide clues to hitherto poorly understood causal pathways.

Hyperketonemia GWAS and parity-dependent SNP associations in Holstein dairy cows intensively sampled for blood ß-hydroxybutyrate concentration.

Hyperketonemia (HYK) is a metabolic disorder that affects early postpartum dairy cows; however, there has been limited success in identifying genomic variants contributing to HYK susceptibility. We conducted a genome-wide association study (GWAS) using HYK phenotypes based on an intensive screening protocol, interrogated genotype interactions with parity group (GWIS), and evaluated the enrichment of annotated metabolic pathways. Holstein cows were enrolled into the experiment after parturition, and blood samples were collected at four timepoints between 5 and 18 days postpartum. Concentration of blood ß-hydroxybutyrate (BHB) was quantified cow-side via a handheld BHB meter. Cows were labeled as a HYK case when at least one blood sample had BHB ≥ 1.2 mmol/L, and all other cows were considered non-HYK controls. After quality control procedures, 1,710 cows and 58,699 genotypes were available for further analysis. The GWAS and GWIS were performed using the forward feature select linear mixed model method. There was evidence for an association between ARS-BFGL-NGS-91238 and HYK susceptibility, as well as parity-dependent associations to HYK for BovineHD0600024247 and BovineHD1400023753. Candidate genes annotated to these single nuclear polymorphism associations have been previously associated with obesity, diabetes, insulin resistance, and fatty liver in humans and rodent models. Enrichment analysis revealed focal adhesion and axon guidance as metabolic pathways contributing to HYK etiology, while genetic variation in pathways related to insulin secretion and sensitivity may affect HYK susceptibility in a parity-dependent matter. In conclusion, the present work proposes several novel marker associations and metabolic pathways contributing to genetic risk for HYK susceptibility.

Does parity matter in women's risk of dementia? A COSMIC collaboration cohort study.

BACKGROUND: Dementia shows sex difference in its epidemiology. Childbirth, a distinctive experience of women, is associated with the risk for various diseases. However, its association with the risk of dementia in women has rarely been studied. METHODS: We harmonized and pooled baseline data from 11 population-based cohorts from 11 countries over 3 continents, including 14,792 women aged 60 years or older. We investigated the association between parity and the risk of dementia using logistic regression models that adjusted for age, educational level, hypertension, diabetes mellitus, and cohort, with additional analyses by region and dementia subtype. RESULTS: Across all cohorts, grand multiparous (5 or more childbirths) women had a 47% greater risk of dementia than primiparous (1 childbirth) women (odds ratio [OR] = 1.47, 95% confidence interval [CI] = 1.10-1.94), while nulliparous (no childbirth) women and women with 2 to 4 childbirths showed a comparable dementia risk to primiparous women. However, there were differences associated with region and dementia subtype. Compared to women with 1 to 4 childbirths, grand multiparous women showed a higher risk of dementia in Europe (OR = 2.99, 95% CI = 1.38-6.47) and Latin America (OR = 1.49, 95% CI = 1.04-2.12), while nulliparous women showed a higher dementia risk in Asia (OR = 2.15, 95% CI = 1.33-3.47). Grand multiparity was associated with 6.9-fold higher risk of vascular dementia in Europe (OR = 6.86, 95% CI = 1.81-26.08), whereas nulliparity was associated with a higher risk of Alzheimer disease (OR = 1.91, 95% CI 1.07-3.39) and non-Alzheimer non-vascular dementia (OR = 3.47, 95% CI = 1.44-8.35) in Asia. CONCLUSION: Parity is associated with women's risk of dementia, though this is not uniform across regions and dementia subtypes.

Genetic parameters estimates for ewes' behavioural reactivity towards their litter after lambing.

In livestock, improving maternal reactivity towards the litter is an important issue in breeding strategies to promote production and animal welfare. As of yet, no studies have investigated the within-breed genetic variation of maternal reactivity in sheep. The objective of this study was to estimate the genetic parameters of maternal reactivity traits. A total of 1,095 primiparous and 1,441 multiparous Romane ewes were phenotyped 24 hr postlambing using a behavioural test (arena test, AT) over a 10-year experimental period. The test consisted of three successive phases evaluating the ewe's attraction to her litter, reactivity to separation from her litter, and reactivity to a conflict between attraction to her litter and avoidance of a motionless human. The ewes were reared exclusively on rangelands (South of France) and lambed outdoors in the spring. High-pitched bleating and low-pitched bleating in the AT were mostly highly heritable (0.39-0.46). Heritabilities were moderate for proximity to the litter in the presence of a human (0.27) and low for locomotion and vigilance in the AT (0.09-0.15). The measurements of a given behaviour in the three phases of the AT were highly genetically correlated. Few genetic correlations were found between the different behavioural traits in the AT, the highest correlations being between high-pitched bleating and low-pitched bleating (-0.43 to -0.77). In conclusion, our findings demonstrate moderate-to-high heritability for maternal reactivity traits. These traits could be included in genetic selection schemes to enhance maternal attachment provided there is no unfavourable link with other production traits.

Genome-wide association for milk production and lactation curve parameters in Holstein dairy cows.

The aim of this study was to identify genomic regions associated with 305-day milk yield and lactation curve parameters on primiparous (n = 9,910) and multiparous (n = 11,158) Holstein cows. The SNP solutions were estimated using a weighted single-step genomic BLUP approach and imputed high-density panel (777k) genotypes. The proportion of genetic variance explained by windows of 50 consecutive SNP (with an average of 165 Kb) was calculated, and regions that accounted for more than 0.50% of the variance were used to search for candidate genes. Estimated heritabilities were 0.37, 0.34, 0.17, 0.12, 0.30 and 0.19, respectively, for 305-day milk yield, peak yield, peak time, ramp, scale and decay for primiparous cows. Genetic correlations of 305-day milk yield with peak yield, peak time, ramp, scale and decay in primiparous cows were 0.99, 0.63, 0.20, 0.97 and -0.52, respectively. The results identified three windows on BTA14 associated with 305-day milk yield and the parameters of lactation curve in primi- and multiparous cows. Previously proposed candidate genes for milk yield supported by this work include GRINA, CYHR1, FOXH1, TONSL, PPP1R16A, ARHGAP39, MAF1, OPLAH and MROH1, whereas newly identified candidate genes are MIR2308, ZNF7, ZNF34, SLURP1, MAFA and KIFC2 (BTA14). The protein lipidation biological process term, which plays a key role in controlling protein localization and function, was identified as the most important term enriched by the identified genes.

Estimation of litter size variability phenotypes in Large White sows.

The objective of this study was to obtain new phenotypes of phenotypic variability for the total number born (TNB) in pigs using the residual variance of TNB. The analysis was based on 246,799 Large White litter observations provided by Topigs Norsvin. Three animal models were used to obtain estimates of residual variance for TNB: the basic model (BM) containing fixed effects of farm-year and season and random effects of animal and permanent environmental sow, the basic model with an additional fixed effect of parity (BMP) and a random regression model (RRM). The within-individual variance of the residuals was calculated and log-transformed to obtain three new variability traits: LnVarBM, LnVarBMP and LnVarRRM. Then, (co)variance components, heritability, the genetic coefficient of variation at the standard deviation level (GCVSDe ) and genetic correlations between the three LnVar's and between the LnVar's and mean total number born (mTNB) were estimated with uni-, bi- and trivariate models. Results indicated that genetically LnVar's are the same trait and are positively correlated with the mTNB (~0.60). Thus, both traits should be included in breeding programmes to avoid an increase in TNB variability while selecting for increased TNB. Heritability of the LnVar's was estimated at 0.021. The GCVSDe for LnVar's showed that a change of 8% in residual standard deviation of TNB could be obtained per generation. Those results indicate that phenotypic variability of litter size is under genetic control, thus it may be improved by selection.

Imprint of parity and age at first pregnancy on the genomic landscape of subsequent breast cancer.

BACKGROUND: Although parity and age at first pregnancy are among the most known extrinsic factors that modulate breast cancer risk, their impact on the biology of subsequent breast cancer has never been explored in depth. Recent data suggest that pregnancy-induced tumor protection is different according to breast cancer subtypes, with parity and young age at first pregnancy being associated with a marked reduction in the risk of developing luminal subtype but not triple negative breast cancer. In this study, we investigated the imprint of parity and age at first pregnancy on the pattern of somatic mutations, somatic copy number alterations, transcriptomic profiles, and tumor immune microenvironment by assessing tumor-infiltrating lymphocytes (TILs) levels of subsequent breast cancer. METHODS: A total of 313 patients with primary breast cancer with available whole genome, RNA sequencing, and TILs data were included in this study. We used a multivariate analysis adjusted for age at diagnosis, pathological stage, molecular subtypes, and histological subtypes. We compared nulliparous vs. parous, late parous vs. early parous, and nulliparous vs. pregnancy-associated breast cancer (PABC) patients. Late and early parous patients were grouped by using the median age at first pregnancy. PABC was defined as patients diagnosed up to 10 years postpartum. RESULTS: Genomic alterations of breast cancer were associated with age at first pregnancy but not with parity status alone. Independently of clinicopathological features, early parous patients developed tumors characterized by a higher number of Indels (Padj = 0.002), a lower frequency of CDH1 mutations (1.2% vs. 12.7%; Padj = 0.013), a higher frequency of TP53 mutations (50% vs. 22.5%; Padj = 0.010), and MYC amplification (28% vs. 7%; Padj = 0.008). PABC were associated with increased TILs infiltration (Padj = 0.0495). CONCLUSIONS: These findings highlight an unprecedented link between reproductive history and the genomic landscape of subsequent breast cancer. We further hypothesize that TP53-mutant premalignant lesions could be less susceptible to the protective effect of an early parity, which might explain the difference of parity-induced protection according to breast cancer subtypes. This work also advocates that reproductive history should be routinely collected in future large-scale genomic studies addressing the biology of female cancers.

Germline TP53 mutation spectrum in Sudanese premenopausal breast cancer patients: correlations with reproductive factors.

PURPOSE: The role of non-genetic factors as modifiers of TP53-related hereditary breast cancer (BC) risk is debated. In this regard, little is known about the impact of germline TP53 mutations on BC in sub-Saharan Africa, where the disease often presents in non-contraceptive multiparous premenopausal women with extended history of breastfeeding. Herein, we report the germline TP53 mutations found in a series of 92 Sudanese premenopausal BC patients characterized for reproductive history. METHODS: The entire TP53 coding sequence, including intron-exon boundaries and UTRs, was analyzed via DHPLC and direct sequencing, and the association of TP53 genotypes with BC risk and with individual lifetime exposures to reproductive factors was investigated with statistical tools. RESULTS: The germline TP53 mutation spectrum comprised 20 variants, 15 in the non-coding and 5 in the coding region. The latter included a deleterious missense mutation, c.817C>T (p.Arg273Cys), in a unique patient, and the common and functionally relevant coding polymorphism at amino acid 72 [Pro72Arg (rs1042522)]. The non-coding mutations included c.919+1G>A, a known deleterious splice site mutation, also in a unique patient. Notably, the 2 carriers of deleterious TP53 mutations clustered in the subset of cases with stronger reproductive history relative to childbearing age. When analyzed in comparison to population controls, the codon 72 polymorphism did not reveal associations with BC. CONCLUSIONS: Our study suggests that the codon 72 Arg>Pro polymorphism is not implicated in premenopausal BC susceptibility, whereas multiparity and breastfeeding might be BC risk factors for carriers of deleterious TP53 mutations.

Reproduction, DNA methylation and biological age.

STUDY QUESTION: Are reproductive characteristics associated with genome-wide DNA methylation and epigenetic age? SUMMARY ANSWER: Our data suggest that increasing parity is associated with differences in blood DNA methylation and small increases in epigenetic age. WHAT IS KNOWN ALREADY: A study of 397 young Filipino women (ages 20-22) observed increasing epigenetic age with an increasing number of pregnancies. STUDY DESIGN, SIZE, DURATION: We used data from 2356 non-Hispanic white women (ages 35-74) enrolled in the Sister Study cohort. PARTICIPANTS/MATERIALS, SETTING, METHODS: Data on reproductive history were ascertained via questionnaire. Of the 2356 women, 1897 (81%) reported at least one live birth. Among parous women, 487 (26%) women reported ever experiencing a pregnancy complication. Three epigenetic clocks (i.e. Hannum, Horvath and Levine) and genome-wide methylation were measured in DNA from whole blood using Illumina's HumanMethylation450 BeadChip. We estimated association ß-values and 95% CIs using linear regression. MAIN RESULTS AND THE ROLE OF CHANCE: All three epigenetic clocks showed weak associations between number of births and epigenetic age (per live birth; Hannum: ß = 0.16, 95% CI = 0.02, 0.29, P = 0.03; Horvath: ß = 0.12, 95% CI = -0.04, 0.27, P = 0.14; Levine: ß = 0.27, 95% CI = 0.08, 0.45, P = 0.01); however, additional adjustment for current BMI attenuated the associations. Among parous women, a history of abnormal glucose tolerance during pregnancy was associated with increased epigenetic age by the Hannum clock (ß = 0.96; 95% CI = 0.10, 1.81; P = 0.03) and Levine clocks (ß = 1.69; 95% CI = 0.54, 2.84; P < 0.01). In epigenome-wide analysis, increasing parity was associated with methylation differences at 17 CpG sites (Bonferroni corrected P≤ 1.0 × 10-7). LIMITATIONS, REASONS FOR CAUTION: We relied on retrospective recall to ascertain reproductive history and pregnancy complications. WIDER IMPLICATIONS OF THE FINDINGS: Our findings suggest that parity is associated with small increases in epigenetic age and with DNA methylation at multiple sites in the genome. STUDY FUNDING/COMPETING INTEREST(S): This research was supported by the Intramural Research program of the NIH, National Institute of Environmental Health Sciences (Z01-ES049033, Z01-ES049032 and Z01-ES044055). None of the authors have a conflict of interest. TRIAL REGISTRATION NUMBER: Not applicable.

NRDR inhibits estradiol synthesis and is associated with changes in reproductive traits in pigs.

Cumulus cells secreting steroid hormones have important functions in oocyte development. Several members of the short-chain dehydrogenase/reductase (SDR) family are critical to the biosynthesis of steroid hormones. NADPH-dependent retinol dehydrogenase/reductase ( NRDR), a member of the SDR superfamily, is overexpressed in pig breeds that also show high levels of androstenone. However, the potential functions and regulatory mechanisms of NRDR in pig ovaries have not been reported to date. The present study demonstrated that NRDR is highly expressed in pig ovaries and is specifically located in cumulus granulosa cells. Functional studies showed that NRDR inhibition increased estradiol synthesis. Both pregnant mare serum gonadotropin and human chorionic gonadotropin downregulated the expression of NRDR in pig cumulus granulosa cells. When the relationship between reproductive traits and single-nucleotide polymorphisms (SNPs) of the NRDR gene was examined, we found that two SNPs affected reproductive traits. SNP rs701332503 was significantly associated with a decrease in the total number of piglets born during multiparity, and rs326982309 was significantly associated with an increase in the average birth weight during primiparity. Thus, NRDR has an important role in steroid hormone biosynthesis in cumulus granulosa cells, and NRDR SNPs are associated with changes in porcine reproduction traits.

Reproductive characteristics are associated with gene-specific promoter methylation status in breast cancer.

BACKGROUND: Reproductive characteristics are well-established risk factors for breast cancer, but the underlying mechanisms are not fully resolved. We hypothesized that altered DNA methylation, measured in tumor tissue, could act in concert with reproductive factors to impact breast carcinogenesis. METHODS: Among a population-based sample of women newly diagnosed with first primary breast cancer, reproductive history was assessed using a life-course calendar approach in an interviewer-administered questionnaire. Methylation-specific polymerase chain reaction and Methyl Light assays were used to assess gene promotor methylation status (methylated vs. unmethylated) for 13 breast cancer-related genes in archived breast tumor tissue. We used case-case unconditional logistic regression to estimate adjusted odds ratios (ORs) and 95% confidence intervals (CIs) for associations with age at menarche and parity (among 855 women), and age at first birth and lactation (among a subset of 736 parous women) in association with methylation status. RESULTS: Age at first birth > 27 years, compared with < 23 years, was associated with lower odds of methylation of CDH1 (OR = 0.44, 95% CI = 0.20-0.99) and TWIST1 (OR = 0.48, 95% CI = 0.28-0.82), and higher odds of methylation of BRCA1 (OR = 1.63, 95% CI = 1.14-2.35). Any vs. no lactation was associated with higher odds of methylation of the PGR gene promoter (OR = 1.59, 95% CI = 1.01-2.49). No associations were noted for parity and methylation in any of the genes assayed. CONCLUSIONS: Our findings indicate that age at first birth, lactation and, perhaps age at menarche, are associated with gene promoter methylation in breast cancer, and should be confirmed in larger studies with robust gene coverage.

A genome wide association study for the number of animals born dead in domestic pigs.

BACKGROUND: The number of animals born dead, which includes the number of mummified (NM) and stillborn (NS) animals, is the most important trait to directly quantify the reproductive loss in domestic pigs. In this study, 282 Landrace sows and 250 Large White sows were genotyped by sequencing (GBS). A total of 816 and 1068 litter records for NM and NS were collected from them. A genome-wide association study (GWAS) was conducted to reveal the genetic difference between NM and NS. RESULTS: A total of 248 and 10 genome-wide significant SNPs were detected for NM and NS across numerous parities in Landrace pigs. The corresponding numbers for Large White pigs were 175 and 6, respectively. All of the detected SNPs were parity specific for both NM and NS in two breeds. Based on significant SNPs, in total 242 (146 for Landrace pig, 96 for Large White pig) and 10 significant chromosome regions (8 for Landrace pigs, 2 for Large White pigs) were found for NM and NS, respectively. Among them, 237 (142 for Landrace pig, 95 for Large White pig) and 8 significant chromosome regions (6 for Landrace pigs, 2 for Large White pigs) for NM and NS were not reported in previous studies. A list of candidate genes at the identified loci was proposed, including HMGB1, SOX5, KCNJ8, ABCC9 and YY1 for NM, ASTN1 for NS. CONCLUSION: This is the first time when GBS data was used to identify genetic regions affecting NM and NS in Landrace and Large White pigs. Many identified informative SNPs and candidate genes advance our understanding of the genetic architecture of NM and NS in pigs. However, further studies are needed to validate using larger populations with more breeds.

Histone deacetylase inhibitor treatment induces postpartum-like maternal behavior and immediate early gene expression in the maternal neural pathway in virgin mice.

The peripartum period is associated with the onset of behaviors that shelter, feed and protect young offspring from harm. The neural pathway that regulates caregiving behaviors has been mapped in female rats and is conserved in mice. However, rats rely on late gestational hormones to shift their perception of infant cues from aversive to attractive, whereas laboratory mice are "spontaneously" maternal, but their level of responding depends on experience. For example, pup-naïve virgin female mice readily care for pups in the home cage, but avoid pups in a novel environment. In contrast, pup-experienced virgin mice care for pups in both contexts. Thus, virgin mice rely on experience to shift their perception of infant cues from aversive to attractive in a novel context. We hypothesize that alterations in immediate early gene activation may underlie the experience-driven shift in which neural pathways (fear/avoidance versus maternal/approach) are activated by pups to modulate context-dependent changes in maternal responding. Here we report that the effects of sodium butyrate, a drug that allows for an amplification of experience-induced histone acetylation and gene expression in virgins, are comparable to the natural onset of caregiving behaviors in postpartum mice and induce postpartum-like patterns of immediate early gene expression across brain regions. These data suggest that pups can activate a fear/defensive circuit in mice and experience-driven improvements in caregiving behavior could be regulated in part through decreased activation of this pathway.

Genetic parameters for female fertility in Nordic Holstein and Red Cattle dairy breeds.

Genetic evaluation of female fertility in Danish, Finnish, and Swedish dairy cows was updated in 2015 to multiple-trait animal model evaluation, where heifer and cow fertility up to third parity are considered as separate traits. A model for conception rate was also developed, which required variance component estimation for Nordic Holstein and Nordic Red Dairy Cattle. We used a multiple-trait multiple-lactation sire model to determine variance components for interval from calving to first insemination, length of service period, and conception rate. Monte Carlo Expectation Maximization REML allowed estimation of all 11 traits simultaneously. Study data were sampled from Swedish Holstein (n = 140,040) and Red Dairy Cattle (n = 101,315) heifers and cows. Conception rate observations are binomial observations with various numbers of failures preceding an observation of success. Using a simulation study, we confirmed that including a service number effect into the conception rate model allowed us to model the change in expectation of successful AI with increasing number of services. Heifers outperformed cows in all fertility traits according to the phenotypic means in the records. Heritabilities for the traits varied from 3 to 7% for interval from calving to first insemination, from 1 to 5% for length of service period, and from 1 to 3% for conception rate. Genetic correlations within traits (i.e., between parities) were favorable, ranging from moderate to high; genetic correlations between heifer and cow traits were lower than between cow traits in different parities. Lowest genetic correlations between traits were for interval from calving to first insemination and conception rate, intermediate for interval from calving to first insemination and length of service period, and highest for length of service period and conception rate. The variance components estimated in this study have been used in Nordic fertility breeding value evaluations since 2016.

Identifying SNPs and candidate genes for three litter traits using single-step GWAS across six parities in Landrace and Large White pigs.

Total number born (TNB), number born alive (NBA), and litter weight born alive (LWB) are critically important traits in pig production. The sow's parity is one of the major factors influencing litter traits. Because of monogenic or polygenic contributions and the presence of temporal gene effects in different sows' parities, it is difficult to clarify the biological and genetic background. To systematically explore the genetic mechanism of litter traits, we conducted 18 GWASs using single-step GWAS (ssGWAS) based on two breeds (908 Landrace and 1,130 Large White sow litter records) for each litter trait in different parities. A total of 300 Landrace and 300 Large White sows were genotyped by sequencing (GBS). ssGWAS was performed separately for each breed and each parity due to population stratification and temporal gene effect. In summary, we identified 80 (15 for Landrace and 65 for Large White), 227 (52 for Landrace, 175 for Large White), and 187 (34 for Landrace, 153 for Large White) single nucleotide polymorphisms (SNPs) affecting TNB, NBA, and LWB, respectively. Of them, we suggest that a total of 22 loci (SSC1: 125098202, SSC1: 117560058, SSC14: 147794697, SSC8: 84823302, SSC9: 143554876, and SSC9: 138766097 for Landrace; SSC1: 4023577, SSC1: 3859573, SSC1: 4891063, SSC16: 5197665, SSC10: 32050819, SSC13: 13552924, SSC13: 92819, SSC17: 3579607, SSC13: 196698221, SSC7: 30918403, SSC16: 46221484, SSC16: 46169204, SSC2: 41988642, SSC2: 44475457, SSC2: 42521875, and SSC7: 58411951 for Large White) are shared by TNB, NBA, and LWB. These results indicate the existence of gene temporal effect in each parity. Furthermore, our findings suggest four interesting candidate genes (FBXL7, ALDH1A2, LEPR, and DDX1) associated with litter traits in different parities that have a major effect on embryonic development progression. In conclusion, 22 crucial SNPs and four interesting candidate genes were identified for three litter traits across six parities. These findings advance our understanding of the genetic architecture of litter traits and confirm the presence of temporal gene effects in different parities. Importantly, functional validation studies for findings of particular interest are recommended in litter traits.

Reproductive risk factors in breast cancer and genetic hormonal pathways: a gene-environment interaction in the MCC-Spain project.

BACKGROUND: Reproductive factors are well known risk factors for breast cancer; however, little is known about how genetic variants in hormonal pathways interact with that relationship. METHODS: One thousand one hundred thirty nine cases of breast cancer in women and 1322 frequency-matched controls were compared. Genetic variants in hormonal pathways (identified in the Kyoto Encyclopedia of Genes and Genomes) were screened according to their relationship with breast cancer using the Cochran-Armitage statistic. Information on reproductive factors was obtained using a face-to-face questionnaire. The interaction among the selected genetic variants and reproductive factors was tested with logistic regression. RESULTS: Concerning C allele in rs2229712, compared to nulliparity in non-carriers the ORs for 1-2 and > 2 deliveries were 0.48 (0.28-0.81) and 0.34 (0.19-0.59), and in C carriers they were 0.92 (0.42-1.98) and 0.71 (0.31-1.61). Similar results were found in women carrying the C allele in rs1269851. Carriers of Allele T in rs35652107 and allele C in rs6018027 had the delivery number effect more pronounced. CONCLUSIONS: The number of deliveries had a dose-response protective effect on breast cancer; women carrying C allele in rs2229712 did not benefit from this protective effect.

FOXA1 hypermethylation: link between parity and ER-negative breast cancer in African American women?

BACKGROUND: Reproductive factors, particularly parity, have differential effects on breast cancer risk according to estrogen receptor (ER) status, especially among African American (AA) women. One mechanism could be through DNA methylation, leading to altered expression levels of genes important in cell fate decisions. METHODS: Using the Illumina 450K BeadChip, we compared DNA methylation levels in paraffin-archived tumor samples from 383 AA and 350 European American (EA) women in the Women's Circle of Health Study (WCHS). We combined 450K profiles with RNA-seq data and prioritized genes based on differential methylation by race, correlation between methylation and gene expression, and biological function. We measured tumor protein expression and assessed its relationship to DNA methylation. We evaluated associations between reproductive characteristics and DNA methylation using linear regression. RESULTS: 410 loci were differentially methylated by race, with the majority unique to ER- tumors. FOXA1 was hypermethylated in tumors from AA versus EA women with ER- cancer, and increased DNA methylation correlated with reduced RNA and protein expression. Importantly, parity was positively associated with FOXA1 methylation among AA women with ER- tumors (P = 0.022), as was number of births (P = 0.026), particularly among those who did not breastfeed (P = 0.008). These same relationships were not observed among EA women, although statistical power was more limited. CONCLUSIONS: Methylation and expression of FOXA1 is likely impacted by parity and breastfeeding. Because FOXA1 regulates a luminal gene expression signature in progenitor cells and represses the basal phenotype, this could be a mechanism that links these reproductive exposures with ER- breast cancer.

Effect of metabolic status on conceptus-maternal interactions on day 19 in dairy cattle: II. Effects on the endometrial transcriptome.

The aim of this study was to test the hypothesis that the metabolic stresses associated with lactation alter the ability of the endometrium to respond appropriately to the conceptus by examining endometrial gene expression on day 19 of pregnancy. Immediately after calving, primiparous Holstein cows with similar production and fertility estimated breeding values were randomly divided into two groups and either dried off (i.e. never milked) immediately or milked twice daily. Approximately 65-75 days postpartum, grade 1 blastocysts recovered from superovulated Holstein heifer donors (n = 5) were transferred (1 per recipient) into lactating (n = 11) and nonlactating (n = 11) recipients. Control nulliparous Holstein heifers (n = 6) were artificially inseminated. RNA-sequencing was performed on intercaruncular endometrial samples recovered at slaughter from confirmed pregnant animals on day 19 (n = 5 lactating and nonlactating cows; n = 4 heifers). Differentially expressed genes (DEGs) were identified between both postpartum groups compared to heifers and between lactating and nonlactating cows. Functional annotation of DEGs between cows and heifers revealed over-representation of categories, including endosome, cytoplasmic vesicle, endocytosis, regulation of exocytosis, and cytokine receptor activity. Functional categories including transcription factor binding sites, cell motility, and cell migration were enriched for DEGs between endometria from lactating and nonlactating cows. In conclusion, while the evidence for a major effect of lactation on the endometrial transcriptome is relatively weak, these data suggest that the metabolic status of the animal (heifer vs cow) modulates the response of the endometrium to the developing conceptus.