Resolving Human Lactation Heterogeneity Using Single Milk-Derived Cells, a Resource at the Ready.

Single cell RNA sequencing (scRNAseq) of human milk-derived cells (HMDCs) makes highly detailed analyses of the biology of human lactation possible. We explore this powerful application as an exciting tool to inspect the cellular composition of human milk. We point out some important challenges unique to this approach and highlight the importance of collaborations between biologists and well-trained bioinformaticians to utilize these data to their maximum potential. We extend this focus by discussing the first two such studies that describe HMDCs via scRNAseq and a variety of important questions in the field that warrant attention through further research. The stage is set to apply scRNAseq in human lactation biology, potentially leading to new insights regarding the molecular and cellular diversity of human secretory mammary epithelial cells.

Single Cell RNA Sequencing of Human Milk-Derived Cells Reveals Sub-Populations of Mammary Epithelial Cells with Molecular Signatures of Progenitor and Mature States: a Novel, Non-invasive Framework for Investigating Human Lactation Physiology.

Cells in human milk are an untapped source, as potential "liquid breast biopsies", of material for investigating lactation physiology in a non-invasive manner. We used single cell RNA sequencing (scRNA-seq) to identify milk-derived mammary epithelial cells (MECs) and their transcriptional signatures in women with diet-controlled gestational diabetes (GDM) with normal lactation. Methodology is described for coordinating milk collections with single cell capture and library preparation via cryopreservation, in addition to scRNA-seq data processing and analyses of MEC transcriptional signatures. We comprehensively characterized 3740 cells from milk samples from two mothers at two weeks postpartum. Most cells (>90%) were luminal MECs (luMECs) expressing lactalbumin alpha and casein beta and positive for keratin 8 and keratin 18. Few cells were keratin 14+ basal MECs and a small immune cell population was present (<10%). Analysis of differential gene expression among clusters identified six potentially distinct luMEC subpopulation signatures, suggesting the potential for subtle functional differences among luMECs, and included one cluster that was positive for both progenitor markers and mature milk transcripts. No expression of pluripotency markers POU class 5 homeobox 1 (POU5F1, encoding OCT4) SRY-box transcription factor 2 (SOX2) or nanog homeobox (NANOG), was observed. These observations were supported by flow cytometric analysis of MECs from mature milk samples from three women with diet-controlled GDM (2-8 mo postpartum), indicating a negligible basal/stem cell population (epithelial cell adhesion molecule (EPCAM)-/integrin subunit alpha 6 (CD49f)+, 0.07%) and a small progenitor population (EPCAM+/CD49f+, 1.1%). We provide a computational framework for others and future studies, as well as report the first milk-derived cells to be analyzed by scRNA-seq. We discuss the clinical potential and current limitations of using milk-derived cells as material for characterizing human mammary physiology.

Compromised global embryonic transcriptome associated with advanced maternal age.

PURPOSE: To investigate the global transcriptome and associated embryonic molecular networks impacted with advanced maternal age (AMA). METHODS: Blastocysts derived from donor oocyte IVF cycles with no male factor infertility (< 30 years of age) and AMA blastocysts (≥ 42 years) with no other significant female factor infertility or male factor infertility were collected with informed patient consent. RNA sequencing libraries were prepared using the SMARTer® Ultra® Low Kit (Clontech Laboratories) and sequenced on the Illumina HiSEQ 4000. Bioinformatics included Ingenuity® Pathway Analysis (Qiagen) with ViiA™ 7 qPCR utilized for gene expression validation (Applied Biosystems). RESULTS: A total of 2688 significant differentially expressed transcripts were identified to distinguish the AMA blastocysts from young, donor controls. 2551 (95%) of these displayed decreased transcription in the blastocysts from older women. Pathway analysis revealed three altered molecular signaling networks known to be critical for embryo and fetal development: CREBBP, ESR1, and SP1. Validation of genes within these networks confirmed the global decreased transcription observed in AMA blastocysts (P < 0.05). CONCLUSIONS: A significant, overall decreased global transcriptome was observed in blastocysts from AMA women. The ESR1/SP1/CREBBP pathway, in particular, was found to be a highly significant upstream regulator impacting biological processes that are vital during embryonic patterning and pre-implantation development. These results provide evidence that AMA embryos are compromised on a cell signaling level which can repress the embryo's ability to proliferate and implant, contributing to a deterioration of reproductive outcomes.

Neonatal intake of Omega-3 fatty acids enhances lipid oxidation in adipocyte precursors.

Establishing metabolic programming begins during fetal and postnatal development, and early-life lipid exposures play a critical role during neonatal adipogenesis. We define how neonatal consumption of a low omega-6 to -3 fatty acid ratio (n6/n3 FA ratio) establishes FA oxidation in adipocyte precursor cells (APCs) before they become adipocytes. In vivo, APCs isolated from mouse pups exposed to the low n6/n3 FA ratio had superior FA oxidation capacity, elevated beige adipocyte mRNAs Ppargc1α, Ucp2, and Runx1, and increased nuclear receptor NR2F2 protein. In vitro, APC treatment with NR2F2 ligand-induced beige adipocyte mRNAs and increased mitochondrial potential but not mass. Single-cell RNA-sequencing analysis revealed low n6/n3 FA ratio yielded more mitochondrial-high APCs and linked APC NR2F2 levels with beige adipocyte signatures and FA oxidation. Establishing beige adipogenesis is of clinical relevance, because fat depots with energetically active, smaller, and more numerous adipocytes improve metabolism and delay metabolic dysfunction.

Onset of taste bud cell renewal starts at birth and coincides with a shift in SHH function.

Embryonic taste bud primordia are specified as taste placodes on the tongue surface and differentiate into the first taste receptor cells (TRCs) at birth. Throughout adult life, TRCs are continually regenerated from epithelial progenitors. Sonic hedgehog (SHH) signaling regulates TRC development and renewal, repressing taste fate embryonically, but promoting TRC differentiation in adults. Here, using mouse models, we show TRC renewal initiates at birth and coincides with onset of SHHs pro-taste function. Using transcriptional profiling to explore molecular regulators of renewal, we identified Foxa1 and Foxa2 as potential SHH target genes in lingual progenitors at birth and show that SHH overexpression in vivo alters FoxA1 and FoxA2 expression relevant to taste buds. We further bioinformatically identify genes relevant to cell adhesion and cell locomotion likely regulated by FOXA1;FOXA2 and show that expression of these candidates is also altered by forced SHH expression. We present a new model where SHH promotes TRC differentiation by regulating changes in epithelial cell adhesion and migration.

Whole Exome Sequencing of 23 Multigeneration Idiopathic Scoliosis Families Reveals Enrichments in Cytoskeletal Variants, Suggests Highly Polygenic Disease.

Adolescent idiopathic scoliosis (AIS) is a lateral spinal curvature >10° with rotation that affects 2-3% of healthy children across populations. AIS is known to have a significant genetic component, and despite a handful of risk loci identified in unrelated individuals by GWAS and next-generation sequencing methods, the underlying etiology of the condition remains largely unknown. In this study, we performed exome sequencing of affected individuals within 23 multigenerational families, with the hypothesis that the occurrence of rare, low frequency, disease-causing variants will co-occur in distantly related, affected individuals. Bioinformatic filtering of uncommon, potentially damaging variants shared by all sequenced family members revealed 1448 variants in 1160 genes across the 23 families, with 132 genes shared by two or more families. Ten genes were shared by >4 families, and no genes were shared by all. Gene enrichment analysis showed an enrichment of variants in cytoskeletal and extracellular matrix related processes. These data support a model that AIS is a highly polygenic disease, with few variant-containing genes shared between affected individuals across different family lineages. This work presents a novel resource for further exploration in familial AIS genetic research.

Activity of Combined Androgen Receptor Antagonism and Cell Cycle Inhibition in Androgen Receptor Positive Triple Negative Breast Cancer.

Triple-negative breast cancer (TNBC) is an aggressive subtype, with a peak recurrence rate within the first few years after diagnosis. Few targeted therapies are available to treat this breast cancer subtype, defined by the lack of estrogen receptor (ER) and progesterone receptor and without amplification of human epidermal growth factor receptor 2 (HER2). Although cell cycle cyclin-dependent kinase (CDK) 4/6 inhibitors are approved for treatment of ER-positive (ER+) breast cancer, they have not proven effective as monotherapy in patients with TNBC. The androgen receptor (AR) has emerged as a therapeutic target in a subset of TNBCs and with significant clinical benefit observed in multiple trials. The purpose of this study was to investigate the preclinical activity of the CDK4/6 inhibitor, abemaciclib, in combination with an agent that targets both androgen biosynthesis and AR activity, seviteronel, using TNBC cell lines expressing high AR, cell line xenografts, and an AR-positive (AR+), androgen-responsive TNBC patient-derived xenograft (PDX). Single-cell RNA sequencing demonstrated heterogeneity in AR levels, even in a highly AR+ cell line, and identified cell cycle pathway activation in ARHigh- versus ARLow-expressing cells. Combination treatment with the cell cycle CDK4/6 inhibitor, abemaciclib, and seviteronel showed synergy in an AR+ TNBC model compared with each drug alone. Although cell cycle inhibitors are FDA approved for use in ER+ breast cancer, our studies suggest that they may also be effective in AR+ TNBC, perhaps combined with AR-targeted agents.

Enhancing Hematopoiesis from Murine Embryonic Stem Cells through MLL1-Induced Activation of a Rac/Rho/Integrin Signaling Axis.

The Mixed Lineage Leukemia (MLL1, KMT2A) gene is critical for development and maintenance of hematopoietic stem cells (HSCs), however, whether this protein is limiting for HSC development is unknown due to lack of physiologic model systems. Here, we develop an MLL1-inducible embryonic stem cell (ESC) system and show that induction of wild-type MLL1 during ESC differentiation selectively increases hematopoietic potential from a transitional c-Kit+/Cd41+ population in the embryoid body and also at sites of hematopoiesis in embryos. Single-cell sequencing analysis illustrates inherent heterogeneity of the c-Kit+/Cd41+ population and demonstrates that MLL1 induction shifts its composition toward multilineage hematopoietic identities. Surprisingly, this does not occur through increasing Hox or other canonical MLL1 targets but through an enhanced Rac/Rho/integrin signaling state, which increases responsiveness to Vla4 ligands and enhances hematopoietic commitment. Together, our data implicate a Rac/Rho/integrin signaling axis in the endothelial to hematopoietic transition and demonstrate that MLL1 actives this axis.