EndoTime: non-categorical timing estimates for luteal endometrium.

STUDY QUESTION: Can the accuracy of timing of luteal phase endometrial biopsies based on urinary ovulation testing be improved by measuring the expression of a small number of genes and a continuous, non-categorical modelling approach? SUMMARY ANSWER: Measuring the expression levels of six genes (IL2RB, IGFBP1, CXCL14, DPP4, GPX3 and SLC15A2) is sufficient to obtain substantially more accurate timing estimates and to assess the reliability of timing estimates for each sample. WHAT IS KNOWN ALREADY: Commercially available endometrial timing approaches based on gene expression require large gene sets and use a categorical approach that classifies samples as pre-receptive, receptive or post-receptive. STUDY DESIGN, SIZE, DURATION: Gene expression was measured by RTq-PCR in different sample sets, comprising a total of 664 endometrial biopsies obtained 4-12 days after a self-reported positive home ovulation test. A further 36 endometrial samples were profiled by RTq-PCR as well as RNA-sequencing. PARTICIPANTS/MATERIALS, SETTING, METHODS: A computational procedure, named 'EndoTime', was established that models the temporal profile of each gene and estimates the timing of each sample. Iterating these steps, temporal profiles are gradually refined as sample timings are being updated, and confidence in timing estimates is increased. After convergence, the method reports updated timing estimates for each sample while preserving the overall distribution of time points. MAIN RESULTS AND THE ROLE OF CHANCE: The Wilcoxon rank-sum test was used to confirm that ordering samples by EndoTime estimates yields sharper temporal expression profiles for held-out genes (not used when determining sample timings) than ordering the same expression values by patient-reported times (GPX3: P < 0.005; CXCL14: P < 2.7e-6; DPP4: P < 3.7e-13). Pearson correlation between EndoTime estimates for the same sample set but based on RTq-PCR or RNA-sequencing data showed a high degree of congruency between the two (P = 8.6e-10, R2 = 0.687). Estimated timings did not differ significantly between control subjects and patients with recurrent pregnancy loss or recurrent implantation failure (P > 0.05). LARGE SCALE DATA: The RTq-PCR data files are available via the GitHub repository for the EndoTime software at https://github.com/AE-Mitchell/EndoTime, as is the code used for pre-processing of RTq-PCR data. The RNA-sequencing data are available on GEO (accession GSE180485). LIMITATIONS, REASONS FOR CAUTION: Timing estimates are informed by glandular gene expression and will only represent the temporal state of other endometrial cell types if in synchrony with the epithelium. Methods that estimate the day of ovulation are still required as these data are essential inputs in our method. Our approach, in its current iteration, performs batch correction such that larger sample batches impart greater accuracy to timing estimations. In theory, our method requires endometrial samples obtained at different days in the luteal phase. In practice, however, this is not a concern as timings based on urinary ovulation testing are associated with a sufficient level of noise to ensure that a variety of time points will be sampled. WIDER IMPLICATIONS OF THE FINDINGS: Our method is the first to assay the temporal state of luteal-phase endometrial samples on a continuous domain. It is freely available with fully shared data and open-source software. EndoTime enables accurate temporal profiling of any gene in luteal endometrial samples for a wide range of research applications and, potentially, clinical use. STUDY FUNDING/COMPETING INTEREST(S): This study was supported by a Wellcome Trust Investigator Award (Grant/Award Number: 212233/Z/18/Z) and the Tommy's National Miscarriage Research Centre. None of the authors have any competing interests. J.L. was funded by the Biotechnology and Biological Sciences Research Council (UK) through the Midlands Integrative Biology Training Partnership (MIBTP, BB/M01116X/1).

Characterization of highly proliferative decidual precursor cells during the window of implantation in human endometrium.

Pregnancy depends on the wholesale transformation of the endometrium, a process driven by differentiation of endometrial stromal cells (EnSC) into specialist decidual cells. Upon embryo implantation, decidual cells impart the tissue plasticity needed to accommodate a rapidly growing conceptus and invading placenta, although the underlying mechanisms are unclear. Here we characterize a discrete population of highly proliferative mesenchymal cells (hPMC) in midluteal human endometrium, coinciding with the window of embryo implantation. Single-cell transcriptomics demonstrated that hPMC express genes involved in chemotaxis and vascular transmigration. Although distinct from resident EnSC, hPMC also express genes encoding pivotal decidual transcription factors and markers, most prominently prolactin. We further show that hPMC are enriched around spiral arterioles, scattered throughout the stroma, and occasionally present in glandular and luminal epithelium. The abundance of hPMC correlated with the in vitro colony-forming unit activity of midluteal endometrium and, conversely, clonogenic cells in culture express a gene signature partially conserved in hPMC. Cross-referencing of single-cell RNA-sequencing data sets indicated that hPMC differentiate into a recently discovered decidual subpopulation in early pregnancy. Finally, we demonstrate that recurrent pregnancy loss is associated with hPMC depletion. Collectively, our findings characterize midluteal hPMC as novel decidual precursors that are likely derived from circulating bone marrow-derived mesenchymal stem/stromal cells and integral to decidual plasticity in pregnancy.

Analysis of chromatin accessibility in decidualizing human endometrial stromal cells.

Spontaneous decidualization of the endometrium in response to progesterone signaling is confined to menstruating species, including humans and other higher primates. During this process, endometrial stromal cells (EnSCs) differentiate into specialized decidual cells that control embryo implantation. We subjected undifferentiated and decidualizing human EnSCs to an assay for transposase accessible chromatin with sequencing (ATAC-seq) to map the underlying chromatin changes. A total of 185,084 open DNA loci were mapped accurately in EnSCs. Altered chromatin accessibility upon decidualization was strongly associated with differential gene expression. Analysis of 1533 opening and closing chromatin regions revealed over-representation of DNA binding motifs for known decidual transcription factors (TFs) and identified putative new regulators. ATAC-seq footprint analysis provided evidence of TF binding at specific motifs. One of the largest footprints involved the most enriched motif-basic leucine zipper-as part of a triple motif that also comprised the estrogen receptor and Pax domain binding sites. Without exception, triple motifs were located within Alu elements, which suggests a role for this primate-specific transposable element (TE) in the evolution of decidual genes. Although other TEs were generally under-represented in open chromatin of undifferentiated EnSCs, several classes contributed to the regulatory DNA landscape that underpins decidual gene expression.-Vrljicak, P., Lucas, E. S., Lansdowne, L., Lucciola, R., Muter, J., Dyer, N. P., Brosens, J. J., Ott, S. Analysis of chromatin accessibility in decidualizing human endometrial stromal cells.

Enhanced Differentiation Potential of Primary Human Endometrial Cells Cultured on 3D Scaffolds.

Novel approaches for culturing primary human cells in vitro are increasingly needed to study cell and tissue physiology and to grow replacement tissue for regenerative medicine. Conventional 2D monolayer cultures of endometrial epithelial and stromal cells fail to replicate the complex 3D architecture of tissue. A fully synthetic scaffold that mimics the microenvironment of the human endometrium can ultimately provide a robust platform for investigating tissue physiology and, hence, take significant steps toward tackling female infertility and IVF failure. In this work, emulsion-templated porous polymers (known as polyHIPEs) were investigated as scaffolds for the culture of primary human endometrial epithelial and stromal cells (HEECs and HESCs). Infiltration of HEECs and HESCs into cell-seeded polyHIPE scaffolds was assessed by histological studies, and phenotype was confirmed by immunostaining. Confocal microscopy revealed that the morphology of HEECs and HESCs is representative of that found in vivo. RNA sequencing was used to investigate transcriptome differences between cells grown on polyHIPE scaffolds and in monolayer cultures. The differentiation status of HEECs and HESCs grown in polyHIPE scaffolds and in monolayer cultures was further evaluated by monitoring the expression of endometrial marker genes. Our observations suggest that a 3D cell culture model that could approximate native human endometrial architecture and function can be developed using tailored polyHIPE scaffolds.

A Notch-dependent transcriptional hierarchy promotes mesenchymal transdifferentiation in the cardiac cushion.

BACKGROUND: Valvuloseptal defects are the most common congenital heart defects. Notch signaling-induced endothelial-to-mesenchymal transition (EMT) in the atrioventricular canal (AVC) cushions at murine embryonic day (E)9.5 is a required step during early valve development. Insights to the transcriptional network that is activated in endocardial cells (EC) during EMT and how these pathways direct valve maturation are lacking. RESULTS: We show that at E11.5, AVC-EC retain the ability to undergo Notch-dependent EMT when explanted on collagen. EC-Notch inhibition at E10.5 blocks expression of known mesenchymal genes in E11.5 AVC-EC. To understand the genetic network and AVC development downstream of Notch signaling beyond E9.5, we constructed Tag-Seq libraries corresponding to different cell types of the E11.5 AVC and atrium in wild-type mice and in EC-Notch inhibited mice. We identified 1,400 potential Notch targets in the AVC-EC, of which 124 are transcription factors (TF). From the 124 TFs, we constructed a transcriptional hierarchy and identify 10 upstream TFs within the network. CONCLUSIONS: We validated 4 of the upstream TFs as Notch targets that are enriched in AVC-EC. Functionally, we show these 4 TFs regulate EMT in AVC explant assays. These novel signaling pathways downstream of Notch are potentially relevant to valve development.

Dynamic chromatin remodeling in cycling human endometrium at single-cell level.

Estrogen-dependent proliferation followed by progesterone-dependent differentiation of the endometrium culminates in a short implantation window. We performed single-cell assay for transposase-accessible chromatin with sequencing on endometrial samples obtained across the menstrual cycle to investigate the regulation of temporal gene networks that control embryo implantation. We identify uniquely accessible chromatin regions in all major cellular constituents of the endometrium, delineate temporal patterns of coordinated chromatin remodeling in epithelial and stromal cells, and gain mechanistic insights into the emergence of a receptive state through integrated analysis of enriched transcription factor (TF) binding sites in dynamic chromatin regions, chromatin immunoprecipitation sequencing analyses, and gene expression data. We demonstrate that the implantation window coincides with pervasive cooption of transposable elements (TEs) into the regulatory chromatin landscape of decidualizing cells and expression of TE-derived transcripts in a spatially defined manner. Our data constitute a comprehensive map of the chromatin changes that control TF activities in a cycling endometrium at cellular resolution.

Lineage skewing and genome instability underlie marrow failure in a zebrafish model of GATA2 deficiency.

Inherited bone marrow failure associated with heterozygous mutations in GATA2 predisposes toward hematological malignancies, but the mechanisms remain poorly understood. Here, we investigate the mechanistic basis of marrow failure in a zebrafish model of GATA2 deficiency. Single-cell transcriptomics and chromatin accessibility assays reveal that loss of gata2a leads to skewing toward the erythroid lineage at the expense of myeloid cells, associated with loss of cebpa expression and decreased PU.1 and CEBPA transcription factor accessibility in hematopoietic stem and progenitor cells (HSPCs). Furthermore, gata2a mutants show impaired expression of npm1a, the zebrafish NPM1 ortholog. Progressive loss of npm1a in HSPCs is associated with elevated levels of DNA damage in gata2a mutants. Thus, Gata2a maintains myeloid lineage priming through cebpa and protects against genome instability and marrow failure by maintaining expression of npm1a. Our results establish a potential mechanism underlying bone marrow failure in GATA2 deficiency.

MicroRNAs emerging coordinate with placental mammals alter pathways in endometrial epithelia important for endometrial function.

We tested the hypothesis that conserved placental mammal-specific microRNAs and their targets facilitate endometrial receptivity to implantation. Expression of miR-340-5p, -542-3p, and -671-5p was regulated by exposure of endometrial epithelial cells to progesterone (10 µg/ml) for 24 h coordinate with 1,713 of their predicted targets. Proteomic analysis of cells transfected with miRNA mimic/inhibitor (48 h: n = 3) revealed 1,745 proteins altered by miR-340-5p (mimic; 1,369, inhibitor; 376) of which 171 were predicted targets and P4-regulated. MiR-542-3p altered 2,353 (mimic; 1,378, inhibitor; 975) 100 which were mirDB predicted, including 46 P4-regulated. MiR-671-5p altered 1,744 proteins (mimic; 1,252, inhibitor; 492) 95 of which were predicted targets and 46 P4-regulated. All miRNAs were detected in luteal phase endometrial biopsies, irrespective of pregnancy outcomes. miR-340-5p expression increased in biopsies from individuals suffering previous and subsequent miscarriage compared to those with subsequent live birth. Dysfunction of these miRNAs and their targets contribute to endometrial-derived recurrent pregnancy loss.

The immune landscape of SARS-CoV-2-associated Multisystem Inflammatory Syndrome in Children (MIS-C) from acute disease to recovery.

Multisystem inflammatory syndrome in children (MIS-C) is a life-threatening disease occurring several weeks after severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Deep immune profiling showed acute MIS-C patients had highly activated neutrophils, classical monocytes and memory CD8+ T-cells, with increased frequencies of B-cell plasmablasts and double-negative B-cells. Post treatment samples from the same patients, taken during symptom resolution, identified recovery-associated immune features including increased monocyte CD163 levels, emergence of a new population of immature neutrophils and, in some patients, transiently increased plasma arginase. Plasma profiling identified multiple features shared by MIS-C, Kawasaki Disease and COVID-19 and that therapeutic inhibition of IL-6 may be preferable to IL-1 or TNF-α. We identified several potential mechanisms of action for IVIG, the most commonly used drug to treat MIS-C. Finally, we showed systemic complement activation with high plasma C5b-9 levels is common in MIS-C suggesting complement inhibitors could be used to treat the disease.

Genomic analysis distinguishes phases of early development of the mouse atrio-ventricular canal.

Valve formation during embryonic heart development involves a complex interplay of regional specification, cell transformations, and remodeling events. While many studies have addressed the role of specific genes during this process, a global understanding of the genetic basis for the regional specification and development of the heart valves is incomplete. We have undertaken genome-wide transcriptional profiling of the developing heart valves in the mouse. Four Serial Analysis of Gene Expression libraries were generated and analyzed from the mouse atrio-ventricular canal (AVC) at embryonic days 9.5-12.5, covering the stages from initiation of endothelial to mesenchymal transition (EMT) through to the beginning of endocardial cushion remodeling. We identified 14 distinct temporal patterns of gene expression during AVC development. These were associated with specific functions and signaling pathway members. We defined the temporal distribution of mesenchyme genes during the EMT process and of specific Notch and transforming growth factor-beta targets. This work provides the first comprehensive temporal dataset during the formation of heart valves. These results identify molecular signatures that distinguish different phases of early heart valve formation allowing gene expression and function to be further investigated.

Impact of Sustained Transforming Growth Factor-ß Receptor Inhibition on Chromatin Accessibility and Gene Expression in Cultured Human Endometrial MSC.

Endometrial mesenchymal stem cells (eMSC) drive the extraordinary regenerative capacity of the human endometrium. Clinical application of eMSC for therapeutic purposes is hampered by spontaneous differentiation and cellular senescence upon large-scale expansion in vitro. A83-01, a selective transforming growth factor-ß receptor (TGFß-R) inhibitor, promotes expansion of eMSC in culture by blocking differentiation and senescence, but the underlying mechanisms are incompletely understood. In this study, we combined RNA-seq and ATAC-seq to study the impact of sustained TGFß-R inhibition on gene expression and chromatin architecture of eMSC. Treatment of primary eMSC with A83-01 for 5 weeks resulted in differential expression of 1,463 genes. Gene ontology analysis showed enrichment of genes implicated in cell growth whereas extracellular matrix genes and genes involved in cell fate commitment were downregulated. ATAC-seq analysis demonstrated that sustained TGFß-R inhibition results in opening and closure of 3,555 and 2,412 chromatin loci, respectively. Motif analysis revealed marked enrichment of retinoic acid receptor (RAR) binding sites, which was paralleled by the induction of RARB, encoding retinoic acid receptor beta (RARß). Selective RARß inhibition attenuated proliferation and clonogenicity of A83-01 treated eMSC. Taken together, our study provides new insights into the gene networks and genome-wide chromatin changes that underpin maintenance of an undifferentiated phenotype of eMSC in prolonged culture.

Recurrent pregnancy loss is associated with a pro-senescent decidual response during the peri-implantation window.

During the implantation window, the endometrium becomes poised to transition to a pregnant state, a process driven by differentiation of stromal cells into decidual cells (DC). Perturbations in this process, termed decidualization, leads to breakdown of the feto-maternal interface and miscarriage, but the underlying mechanisms are poorly understood. Here, we reconstructed the decidual pathway at single-cell level in vitro and demonstrate that stromal cells first mount an acute stress response before emerging as DC or senescent DC (snDC). In the absence of immune cell-mediated clearance of snDC, secondary senescence transforms DC into progesterone-resistant cells that abundantly express extracellular matrix remodelling factors. Additional single-cell analysis of midluteal endometrium identified DIO2 and SCARA5 as marker genes of a diverging decidual response in vivo. Finally, we report a conspicuous link between a pro-senescent decidual response in peri-implantation endometrium and recurrent pregnancy loss, suggesting that pre-pregnancy screening and intervention may reduce the burden of miscarriage.

The Glycosyltransferase EOGT Regulates Adropin Expression in Decidualizing Human Endometrium.

In pregnancy, resistance of endometrial decidual cells to stress signals is critical for the integrity of the fetomaternal interface and, by extension, survival of the conceptus. O-GlcNAcylation is an essential posttranslational modification that links glucose sensing to cellular stress resistance. Unexpectedly, decidualization of primary endometrial stromal cells (EnSCs) was associated with a 60% reduction in O-linked ß-N-acetylglucosamine (O-GlcNAc)‒modified proteins, reflecting downregulation of the enzyme that adds O-GlcNAc to substrates (O-GlcNAc transferase; OGT) but not the enzyme that removes the modification (O-GlcNAcase). Notably, epidermal growth factor domain-specific O-linked GlcNAc transferase (EOGT), an endoplasmic reticulum-specific OGT that modifies a limited number of secreted and membrane proteins, was markedly induced in differentiating EnSCs. Knockdown of EOGT perturbed a network of decidual genes involved in multiple cellular functions. The most downregulated gene upon EOGT knockdown in decidualizing cells was the energy homeostasis-associated gene (ENHO), which encodes adropin, a metabolic hormone involved in energy homeostasis and glucose and fatty acid metabolism. Analysis of midluteal endometrial biopsies revealed an inverse correlation between endometrial EOGT and ENHO expression and body mass index. Taken together, our findings revealed that obesity impairs the EOGT-adropin axis in decidual cells, which in turn points toward a mechanistic link between metabolic disorders and adverse pregnancy outcome.

Clearance of senescent decidual cells by uterine natural killer cells in cycling human endometrium.

In cycling human endometrium, menstruation is followed by rapid estrogen-dependent growth. Upon ovulation, progesterone and rising cellular cAMP levels activate the transcription factor Forkhead box O1 (FOXO1) in endometrial stromal cells (EnSCs), leading to cell cycle exit and differentiation into decidual cells that control embryo implantation. Here we show that FOXO1 also causes acute senescence of a subpopulation of decidualizing EnSCs in an IL-8 dependent manner. Selective depletion or enrichment of this subpopulation revealed that decidual senescence drives the transient inflammatory response associated with endometrial receptivity. Further, senescent cells prevent differentiation of endometrial mesenchymal stem cells in decidualizing cultures. As the cycle progresses, IL-15 activated uterine natural killer (uNK) cells selectively target and clear senescent decidual cells through granule exocytosis. Our findings reveal that acute decidual senescence governs endometrial rejuvenation and remodeling at embryo implantation, and suggest a critical role for uNK cells in maintaining homeostasis in cycling endometrium.

Loss of Endometrial Sodium Glucose Cotransporter SGLT1 is Detrimental to Embryo Survival and Fetal Growth in Pregnancy.

Embryo implantation requires a hospitable uterine environment. A key metabolic change that occurs during the peri-implantation period, and throughout early pregnancy, is the rise in endometrial glycogen content. Glycogen accumulation requires prior cellular uptake of glucose. Here we show that both human and murine endometrial epithelial cells express the high affinity Na+-coupled glucose carrier SGLT1. Ussing chamber experiments revealed electrogenic glucose transport across the endometrium in wild type (Slc5a1 +/+) but not in SGLT1 deficient (Slc5a1 -/-) mice. Endometrial glycogen content, litter size and weight of offspring at birth were significantly lower in Slc5a1 -/- mice. In humans, SLC5A1 expression was upregulated upon decidualization of primary endometrial stromal cells. Endometrial SLC5A1 expression during the implantation window was attenuated in patients with recurrent pregnancy loss when compared with control subjects. Our findings reveal a novel mechanism establishing adequate endometrial glycogen stores for pregnancy. Disruption of this histiotrophic pathway leads to adverse pregnancy outcome.

Genome-Wide Analysis of Transposon and Retroviral Insertions Reveals Preferential Integrations in Regions of DNA Flexibility.

DNA transposons and retroviruses are important transgenic tools for genome engineering. An important consideration affecting the choice of transgenic vector is their insertion site preferences. Previous large-scale analyses of Ds transposon integration sites in plants were done on the basis of reporter gene expression or germ-line transmission, making it difficult to discern vertebrate integration preferences. Here, we compare over 1300 Ds transposon integration sites in zebrafish with Tol2 transposon and retroviral integration sites. Genome-wide analysis shows that Ds integration sites in the presence or absence of marker selection are remarkably similar and distributed throughout the genome. No strict motif was found, but a preference for structural features in the target DNA associated with DNA flexibility (Twist, Tilt, Rise, Roll, Shift, and Slide) was observed. Remarkably, this feature is also found in transposon and retroviral integrations in maize and mouse cells. Our findings show that structural features influence the integration of heterologous DNA in genomes, and have implications for targeted genome engineering.

A Multifunctional Mutagenesis System for Analysis of Gene Function in Zebrafish.

Since the sequencing of the human reference genome, many human disease-related genes have been discovered. However, understanding the functions of all the genes in the genome remains a challenge. The biological activities of these genes are usually investigated in model organisms such as mice and zebrafish. Large-scale mutagenesis screens to generate disruptive mutations are useful for identifying and understanding the activities of genes. Here, we report a multifunctional mutagenesis system in zebrafish using the maize Ds transposon. Integration of the Ds transposable element containing an mCherry reporter for protein trap events and an EGFP reporter for enhancer trap events produced a collection of transgenic lines marking distinct cell and tissue types, and mutagenized genes in the zebrafish genome by trapping and prematurely terminating endogenous protein coding sequences. We obtained 642 zebrafish lines with dynamic reporter gene expression. The characterized fish lines with specific expression patterns will be made available through the European Zebrafish Resource Center (EZRC), and a database of reporter expression is available online (http://fishtrap.warwick.ac.uk/). Our approach complements other efforts using zebrafish to facilitate functional genomic studies in this model of human development and disease.

Twist1 transcriptional targets in the developing atrio-ventricular canal of the mouse.

Malformations of the cardiovascular system are the most common type of birth defect in humans, frequently affecting the formation of valves and septa. During heart valve and septa formation, cells from the atrio-ventricular canal (AVC) and outflow tract (OFT) regions of the heart undergo an epithelial-to-mesenchymal transformation (EMT) and invade the underlying extracellular matrix to give rise to endocardial cushions. Subsequent maturation of newly formed mesenchyme cells leads to thin stress-resistant leaflets. TWIST1 is a basic helix-loop-helix transcription factor expressed in newly formed mesenchyme cells of the AVC and OFT that has been shown to play roles in cell survival, cell proliferation and differentiation. However, the downstream targets of TWIST1 during heart valve formation remain unclear. To identify genes important for heart valve development downstream of TWIST1, we performed global gene expression profiling of AVC, OFT, atria and ventricles of the embryonic day 10.5 mouse heart by tag-sequencing (Tag-seq). Using this resource we identified a novel set of 939 genes, including 123 regulators of transcription, enriched in the valve forming regions of the heart. We compared these genes to a Tag-seq library from the Twist1 null developing valves revealing significant gene expression changes. These changes were consistent with a role of TWIST1 in controlling differentiation of mesenchymal cells following their transformation from endothelium in the mouse. To study the role of TWIST1 at the DNA level we performed chromatin immunoprecipitation and identified novel direct targets of TWIST1 in the developing heart valves. Our findings support a role for TWIST1 in the differentiation of AVC mesenchyme post-EMT in the mouse, and suggest that TWIST1 can exert its function by direct DNA binding to activate valve specific gene expression.

A mouse atlas of gene expression: large-scale digital gene-expression profiles from precisely defined developing C57BL/6J mouse tissues and cells.

We analyzed 8.55 million LongSAGE tags generated from 72 libraries. Each LongSAGE library was prepared from a different mouse tissue. Analysis of the data revealed extensive overlap with existing gene data sets and evidence for the existence of approximately 24,000 previously undescribed genomic loci. The visual cortex, pancreas, mammary gland, preimplantation embryo, and placenta contain the largest number of differentially expressed transcripts, 25% of which are previously undescribed loci.