Dual roles for LUBAC signaling in thymic epithelial cell development and survival.

Thymic epithelial cells (TECs) form a unique microenvironment that orchestrates T cell differentiation and immunological tolerance. Despite the importance of TECs for adaptive immunity, there is an incomplete understanding of the signalling networks that support their differentiation and survival. We report that the linear ubiquitin chain assembly complex (LUBAC) is essential for medullary TEC (mTEC) differentiation, cortical TEC survival and prevention of premature thymic atrophy. TEC-specific loss of LUBAC proteins, HOIL-1 or HOIP, severely impaired expansion of the thymic medulla and AIRE-expressing cells. Furthermore, HOIL-1-deficiency caused early thymic atrophy due to Caspase-8/MLKL-dependent apoptosis/necroptosis of cortical TECs. By contrast, deficiency in the LUBAC component, SHARPIN, caused relatively mild defects only in mTECs. These distinct roles for LUBAC components in TECs correlate with their function in linear ubiquitination, NFκB activation and cell survival. Thus, our findings reveal dual roles for LUBAC signaling in TEC differentiation and survival.

Nfkb2 variants reveal a p100-degradation threshold that defines autoimmune susceptibility.

NF-κB2/p100 (p100) is an inhibitor of κB (IκB) protein that is partially degraded to produce the NF-κB2/p52 (p52) transcription factor. Heterozygous NFKB2 mutations cause a human syndrome of immunodeficiency and autoimmunity, but whether autoimmunity arises from insufficiency of p52 or IκB function of mutated p100 is unclear. Here, we studied mice bearing mutations in the p100 degron, a domain that harbors most of the clinically recognized mutations and is required for signal-dependent p100 degradation. Distinct mutations caused graded increases in p100-degradation resistance. Severe p100-degradation resistance, due to inheritance of one highly degradation-resistant allele or two subclinical alleles, caused thymic medullary hypoplasia and autoimmune disease, whereas the absence of p100 and p52 did not. We inferred a similar mechanism occurs in humans, as the T cell receptor repertoires of affected humans and mice contained a hydrophobic signature of increased self-reactivity. Autoimmunity in autosomal dominant NFKB2 syndrome arises largely from defects in nonhematopoietic cells caused by the IκB function of degradation-resistant p100.

A missense mutation in the MLKL brace region promotes lethal neonatal inflammation and hematopoietic dysfunction.

MLKL is the essential effector of necroptosis, a form of programmed lytic cell death. We have isolated a mouse strain with a single missense mutation, MlklD139V, that alters the two-helix 'brace' that connects the killer four-helix bundle and regulatory pseudokinase domains. This confers constitutive, RIPK3 independent killing activity to MLKL. Homozygous mutant mice develop lethal postnatal inflammation of the salivary glands and mediastinum. The normal embryonic development of MlklD139V homozygotes until birth, and the absence of any overt phenotype in heterozygotes provides important in vivo precedent for the capacity of cells to clear activated MLKL. These observations offer an important insight into the potential disease-modulating roles of three common human MLKL polymorphisms that encode amino acid substitutions within or adjacent to the brace region. Compound heterozygosity of these variants is found at up to 12-fold the expected frequency in patients that suffer from a pediatric autoinflammatory disease, chronic recurrent multifocal osteomyelitis (CRMO).

Isolation and Propagation of Mammary Epithelial Stem and Progenitor Cells.

Several methods of mammary gland dissociation have been described that utilize a combined strategy of mechanical and enzymatic dissociation to isolate mammary epithelial cells (MECs) from intact tissue (Smalley et al., J Mammary Gland Biol Neoplasia 17:91-97, 2012). Here we detail a robust method that enables the isolation of all major stem and progenitor MEC populations, which has been successfully used to study stem cell behavior when coupled with transplantation and in vitro assays (Shackleton et al., Nature 439:84-88, 2006; Bouras et al., Cell Stem Cell 3:429-441, 2008; Sheridan et al., BMC Cancer 15:221, 2015; Jamieson et al., Development 144:1065-1071, 2017). Furthermore, we outline two prominent methods for culturing MECs for the purposes of ex vivo manipulation or study: 2D feeder layer cultures and 3D Matrigel colony assays. Importantly, all outlined methods retain stem and progenitor cell behaviors and can be used in combination with downstream in vivo, in vitro, or in silico analyses.

Thymospheres Are Formed by Mesenchymal Cells with the Potential to Generate Adipocytes, but Not Epithelial Cells.

Evidence suggests that a stem-cell-driven differentiation hierarchy maintains the dynamic thymic epithelial cell (TEC) network that governs T lymphocyte development. The identification of TEC stem/progenitor cells has been a major focus in the field, and several candidates with contrasting phenotypes have been described. We sought to determine the provenance and function of the only population reported to exhibit TEC stem cell properties in the adult, a Foxn1- EpCAM- cell that generates so-called thymospheres. We provide evidence that the thymosphere-forming cell (TSFC) is not a TEC stem cell but can incorporate bystander TECs into thymospheres, providing an explanation for the epithelial activity ascribed to these structures. TSFCs were found to share a phenotype, transcriptional profile, and developmental origin with thymic fibroblasts and can generate adipocytes. In summary, this study redefines the nature of bipotent TEC stem/progenitor cells in the adult thymus and highlights a potentially important mesenchymal progenitor population.

A critical epithelial survival axis regulated by MCL-1 maintains thymic function in mice.

T-cell differentiation is governed by interactions with thymic epithelial cells (TECs) and defects in this process undermine immune function and tolerance. To uncover new strategies to restore thymic function and adaptive immunity in immunodeficiency, we sought to determine the molecular mechanisms that control life and death decisions in TECs. Guided by gene expression profiling, we created mouse models that specifically deleted prosurvival genes in TECs. We found that although BCL-2 and BCL-XL were dispensable for TEC homeostasis, MCL-1 deficiency impacted on TECs as early as embryonic day 15.5, resulting in early thymic atrophy and T-cell lymphopenia, with near complete loss of thymic tissue by 2 months of age. MCL-1 was not necessary for TEC differentiation but was continually required for the survival of mature cortical and medullary TECs and the maintenance of thymic architecture. A screen of TEC trophic factors in organ cultures showed that epidermal growth factor upregulated MCL-1 via MAPK/ERK kinase activity, providing a molecular mechanism for the support of TEC survival. This signaling axis governing TEC survival and thymic function represents a new target for strategies for thymic protection and regeneration.

Easy and efficient ensemble gene set testing with EGSEA.

Gene set enrichment analysis is a popular approach for prioritising the biological processes perturbed in genomic datasets. The Bioconductor project hosts over 80 software packages capable of gene set analysis. Most of these packages search for enriched signatures amongst differentially regulated genes to reveal higher level biological themes that may be missed when focusing only on evidence from individual genes. With so many different methods on offer, choosing the best algorithm and visualization approach can be challenging. The EGSEA package solves this problem by combining results from up to 12 prominent gene set testing algorithms to obtain a consensus ranking of biologically relevant results.This workflow demonstrates how EGSEA can extend limma-based differential expression analyses for RNA-seq and microarray data using experiments that profile 3 distinct cell populations important for studying the origins of breast cancer. Following data normalization and set-up of an appropriate linear model for differential expression analysis, EGSEA builds gene signature specific indexes that link a wide range of mouse or human gene set collections obtained from MSigDB, GeneSetDB and KEGG to the gene expression data being investigated. EGSEA is then configured and the ensemble enrichment analysis run, returning an object that can be queried using several S4 methods for ranking gene sets and visualizing results via heatmaps, KEGG pathway views, GO graphs, scatter plots and bar plots. Finally, an HTML report that combines these displays can fast-track the sharing of results with collaborators, and thus expedite downstream biological validation. EGSEA is simple to use and can be easily integrated with existing gene expression analysis pipelines for both human and mouse data.

Combining multiple tools outperforms individual methods in gene set enrichment analyses.

Motivation: Gene set enrichment (GSE) analysis allows researchers to efficiently extract biological insight from long lists of differentially expressed genes by interrogating them at a systems level. In recent years, there has been a proliferation of GSE analysis methods and hence it has become increasingly difficult for researchers to select an optimal GSE tool based on their particular dataset. Moreover, the majority of GSE analysis methods do not allow researchers to simultaneously compare gene set level results between multiple experimental conditions. Results: The ensemble of genes set enrichment analyses (EGSEA) is a method developed for RNA-sequencing data that combines results from twelve algorithms and calculates collective gene set scores to improve the biological relevance of the highest ranked gene sets. EGSEA's gene set database contains around 25 000 gene sets from sixteen collections. It has multiple visualization capabilities that allow researchers to view gene sets at various levels of granularity. EGSEA has been tested on simulated data and on a number of human and mouse datasets and, based on biologists' feedback, consistently outperforms the individual tools that have been combined. Our evaluation demonstrates the superiority of the ensemble approach for GSE analysis, and its utility to effectively and efficiently extrapolate biological functions and potential involvement in disease processes from lists of differentially regulated genes. Availability and Implementation: EGSEA is available as an R package at http://www.bioconductor.org/packages/EGSEA/ . The gene sets collections are available in the R package EGSEAdata from http://www.bioconductor.org/packages/EGSEAdata/ . Contacts: monther.alhamdoosh@csl.com.au mritchie@wehi.edu.au. Supplementary information: Supplementary data are available at Bioinformatics online.

Linear ubiquitin chain assembly complex coordinates late thymic T-cell differentiation and regulatory T-cell homeostasis.

The linear ubiquitin chain assembly complex (LUBAC) is essential for innate immunity in mice and humans, yet its role in adaptive immunity is unclear. Here we show that the LUBAC components HOIP, HOIL-1 and SHARPIN have essential roles in late thymocyte differentiation, FOXP3+ regulatory T (Treg)-cell development and Treg cell homeostasis. LUBAC activity is not required to prevent TNF-induced apoptosis or necroptosis but is necessary for the transcriptional programme of the penultimate stage of thymocyte differentiation. Treg cell-specific ablation of HOIP causes severe Treg cell deficiency and lethal immune pathology, revealing an ongoing requirement of LUBAC activity for Treg cell homeostasis. These data reveal stage-specific requirements for LUBAC in coordinating the signals required for T-cell differentiation.

Setdb1-mediated H3K9 methylation is enriched on the inactive X and plays a role in its epigenetic silencing.

BACKGROUND: The presence of histone 3 lysine 9 (H3K9) methylation on the mouse inactive X chromosome has been controversial over the last 15 years, and the functional role of H3K9 methylation in X chromosome inactivation in any species has remained largely unexplored. RESULTS: Here we report the first genomic analysis of H3K9 di- and tri-methylation on the inactive X: we find they are enriched at the intergenic, gene poor regions of the inactive X, interspersed between H3K27 tri-methylation domains found in the gene dense regions. Although H3K9 methylation is predominantly non-genic, we find that depletion of H3K9 methylation via depletion of H3K9 methyltransferase Set domain bifurcated 1 (Setdb1) during the establishment of X inactivation, results in failure of silencing for around 150 genes on the inactive X. By contrast, we find a very minor role for Setdb1-mediated H3K9 methylation once X inactivation is fully established. In addition to failed gene silencing, we observed a specific failure to silence X-linked long-terminal repeat class repetitive elements. CONCLUSIONS: Here we have shown that H3K9 methylation clearly marks the murine inactive X chromosome. The role of this mark is most apparent during the establishment phase of gene silencing, with a more muted effect on maintenance of the silent state. Based on our data, we hypothesise that Setdb1-mediated H3K9 methylation plays a role in epigenetic silencing of the inactive X via silencing of the repeats, which itself facilitates gene silencing through alterations to the conformation of the whole inactive X chromosome.

Integration of microRNA signatures of distinct mammary epithelial cell types with their gene expression and epigenetic portraits.

INTRODUCTION: MicroRNAs (miRNAs) have been implicated in governing lineage specification and differentiation in multiple organs; however, little is known about their specific roles in mammopoiesis. We have determined the global miRNA expression profiles of functionally distinct epithelial subpopulations in mouse and human mammary tissue, and compared these to their cognate transcriptomes and epigenomes. Finally, the human miRNA signatures were used to interrogate the different subtypes of breast cancer, with a view to determining miRNA networks deregulated during oncogenesis. METHODS: RNA from sorted mouse and human mammary cell subpopulations was subjected to miRNA expression analysis using the TaqMan MicroRNA Array. Differentially expressed (DE) miRNAs were correlated with gene expression and histone methylation profiles. Analysis of miRNA signatures of the intrinsic subtypes of breast cancer in The Cancer Genome Atlas (TCGA) database versus those of normal human epithelial subpopulations was performed. RESULTS: Unique miRNA signatures characterized each subset (mammary stem cell (MaSC)/basal, luminal progenitor, mature luminal, stromal), with a high degree of conservation across species. Comparison of miRNA and transcriptome profiles for the epithelial subtypes revealed an inverse relationship and pinpointed key developmental genes. Interestingly, expression of the primate-specific miRNA cluster (19q13.4) was found to be restricted to the MaSC/basal subset. Comparative analysis of miRNA signatures with H3 lysine modification maps of the different epithelial subsets revealed a tight correlation between active or repressive marks for the top DE miRNAs, including derepression of miRNAs in Ezh2-deficient cellular subsets. Interrogation of TCGA-identified miRNA profiles with the miRNA signatures of different human subsets revealed specific relationships. CONCLUSIONS: The derivation of global miRNA expression profiles for the different mammary subpopulations provides a comprehensive resource for understanding the interplay between miRNA networks and target gene expression. These data have highlighted lineage-specific miRNAs and potential miRNA-mRNA networks, some of which are disrupted in neoplasia. Furthermore, our findings suggest that key developmental miRNAs are regulated by global changes in histone modification, thus linking the mammary epigenome with genome-wide changes in the expression of genes and miRNAs. Comparative miRNA signature analyses between normal breast epithelial cells and breast tumors confirmed an important linkage between luminal progenitor cells and basal-like tumors.

A pooled shRNA screen for regulators of primary mammary stem and progenitor cells identifies roles for Asap1 and Prox1.

BACKGROUND: The molecular regulators that orchestrate stem cell renewal, proliferation and differentiation along the mammary epithelial hierarchy remain poorly understood. Here we have performed a large-scale pooled RNAi screen in primary mouse mammary stem cell (MaSC)-enriched basal cells using 1295 shRNAs against genes principally involved in transcriptional regulation. METHODS: MaSC-enriched basal cells transduced with lentivirus pools carrying shRNAs were maintained as non-adherent mammospheres, a system known to support stem and progenitor cells. Integrated shRNAs that altered culture kinetics were identified by next generation sequencing as relative frequency changes over time. RNA-seq-based expression profiling coupled with in vitro progenitor and in vivo transplantation assays was used to confirm a role for candidate genes in mammary stem and/or progenitor cells. RESULTS: Utilizing a mammosphere-based assay, the screen identified several candidate regulators. Although some genes had been previously implicated in mammary gland development, the vast majority of genes uncovered have no known function within the mammary gland. RNA-seq analysis of freshly purified primary mammary epithelial populations and short-term cultured mammospheres was used to confirm the expression of candidate regulators. Two genes, Asap1 and Prox1, respectively implicated in breast cancer metastasis and progenitor cell function in other systems, were selected for further analysis as their roles in the normal mammary gland were unknown. Both Prox1 and Asap1 were shown to act as negative regulators of progenitor activity in vitro, and Asap1 knock-down led to a marked increase in repopulating activity in vivo, implying a role in stem cell activity. CONCLUSIONS: This study has revealed a number of novel genes that influence the activity or survival of mammary stem and/or progenitor cells. Amongst these, we demonstrate that Prox1 and Asap1 behave as negative regulators of mammary stem/progenitor function. Both of these genes have also been implicated in oncogenesis. Our findings provide proof of principle for the use of short-term cultured primary MaSC/basal cells in functional RNAi screens.

Scribble modulates the MAPK/Fra1 pathway to disrupt luminal and ductal integrity and suppress tumour formation in the mammary gland.

Polarity coordinates cell movement, differentiation, proliferation and apoptosis to build and maintain complex epithelial tissues such as the mammary gland. Loss of polarity and the deregulation of these processes are critical events in malignant progression but precisely how and at which stage polarity loss impacts on mammary development and tumourigenesis is unclear. Scrib is a core polarity regulator and tumour suppressor gene however to date our understanding of Scrib function in the mammary gland has been limited to cell culture and transplantation studies of cell lines. Utilizing a conditional mouse model of Scrib loss we report for the first time that Scrib is essential for mammary duct morphogenesis, mammary progenitor cell fate and maintenance, and we demonstrate a critical and specific role for Scribble in the control of the early steps of breast cancer progression. In particular, Scrib-deficiency significantly induced Fra1 expression and basal progenitor clonogenicity, which resulted in fully penetrant ductal hyperplasia characterized by high cell turnover, MAPK hyperactivity, frank polarity loss with mixing of apical and basolateral membrane constituents and expansion of atypical luminal cells. We also show for the first time a role for Scribble in mammalian spindle orientation with the onset of mammary hyperplasia being associated with aberrant luminal cell spindle orientation and a failure to apoptose during the final stage of duct tubulogenesis. Restoring MAPK/Fra1 to baseline levels prevented Scrib-hyperplasia, whereas persistent Scrib deficiency induced alveolar hyperplasia and increased the incidence, onset and grade of mammary tumours. These findings, based on a definitive genetic mouse model provide fundamental insights into mammary duct maturation and homeostasis and reveal that Scrib loss activates a MAPK/Fra1 pathway that alters mammary progenitor activity to drive premalignancy and accelerate tumour progression.

edgeR: a versatile tool for the analysis of shRNA-seq and CRISPR-Cas9 genetic screens.

Pooled library sequencing screens that perturb gene function in a high-throughput manner are becoming increasingly popular in functional genomics research. Irrespective of the mechanism by which loss of function is achieved, via either RNA interference using short hairpin RNAs (shRNAs) or genetic mutation using single guide RNAs (sgRNAs) with the CRISPR-Cas9 system, there is a need to establish optimal analysis tools to handle such data. Our open-source processing pipeline in edgeR provides a complete analysis solution for screen data, that begins with the raw sequence reads and ends with a ranked list of candidate genes for downstream biological validation. We first summarize the raw data contained in a fastq file into a matrix of counts (samples in the columns, genes in the rows) with options for allowing mismatches and small shifts in sequence position. Diagnostic plots, normalization and differential representation analysis can then be performed using established methods to prioritize results in a statistically rigorous way, with the choice of either the classic exact testing methodology or generalized linear modeling that can handle complex experimental designs. A detailed users' guide that demonstrates how to analyze screen data in edgeR along with a point-and-click implementation of this workflow in Galaxy are also provided. The edgeR package is freely available from http://www.bioconductor.org.

Generation of tissue organoids by compaction reaggregation.

Cellular reaggregation methods are commonly used to generate tissue organoids for use in biological studies. Using a modified method termed "compaction reaggregation," it is possible to establish reaggregates of reproducible size from defined input cell numbers with ease and without specialist equipment. Importantly, this method is suitable for the study of tissues that have proved refractory to reaggregation by other methods. With the option of juxtaposing cell populations, this method is useful for studies of tissue organization and structure.

Global changes in the mammary epigenome are induced by hormonal cues and coordinated by Ezh2.

The mammary epithelium is a dynamic, highly hormone-responsive tissue. To explore chromatin modifications underlying its lineage specification and hormone responsiveness, we determined genome-wide histone methylation profiles of mammary epithelial subpopulations in different states. The marked differences in H3K27 trimethylation between subpopulations in the adult gland suggest that epithelial cell-fate decisions are orchestrated by polycomb-complex-mediated repression. Remarkably, the mammary epigenome underwent highly specific changes in different hormonal contexts, with a profound change being observed in the global H3K27me3 map of luminal cells during pregnancy. We therefore examined the role of the key H3K27 methyltransferase Ezh2 in mammary physiology. Its expression and phosphorylation coincided with H3K27me3 modifications and peaked during pregnancy, driven in part by progesterone. Targeted deletion of Ezh2 impaired alveologenesis during pregnancy, preventing lactation, and drastically reduced stem/progenitor cell numbers. Taken together, these findings reveal that Ezh2 couples hormonal stimuli to epigenetic changes that underpin progenitor activity, lineage specificity, and alveolar expansion in the mammary gland.

Isolation of mouse mammary epithelial subpopulations: a comparison of leading methods.

Isolation of mammary epithelial subpopulations, including stem and progenitor cells, has become a standard technique in recent years. However, a number of methods and approaches for this have developed and the relative benefits of the different approaches, and the reason for their development, have not always been clear. Here, three of the leading laboratories working on the separation of mammary cell subpopulations have summarised their methods, highlighted their differences and similarities and also discussed the reasoning behind the approaches they have taken. This article will assist workers establishing mammary cell separation protocols in their laboratories to make informed choices about the methods they should use.

Foxn1 regulates lineage progression in cortical and medullary thymic epithelial cells but is dispensable for medullary sublineage divergence.

The forkhead transcription factor Foxn1 is indispensable for thymus development, but the mechanisms by which it mediates thymic epithelial cell (TEC) development are poorly understood. To examine the cellular and molecular basis of Foxn1 function, we generated a novel and revertible hypomorphic allele of Foxn1. By varying levels of its expression, we identified a number of features of the Foxn1 system. Here we show that Foxn1 is a powerful regulator of TEC differentiation that is required at multiple intermediate stages of TE lineage development in the fetal and adult thymus. We find no evidence for a role for Foxn1 in TEC fate-choice. Rather, we show it is required for stable entry into both the cortical and medullary TEC differentiation programmes and subsequently is needed at increasing dosage for progression through successive differentiation states in both cortical and medullary TEC. We further demonstrate regulation by Foxn1 of a suite of genes with diverse roles in thymus development and/or function, suggesting it acts as a master regulator of the core thymic epithelial programme rather than regulating a particular aspect of TEC biology. Overall, our data establish a genetics-based model of cellular hierarchies in the TE lineage and provide mechanistic insight relating titration of a single transcription factor to control of lineage progression. Our novel revertible hypomorph system may be similarly applied to analyzing other regulators of development.

Control of mammary stem cell function by steroid hormone signalling.

The ovarian hormones oestrogen and progesterone profoundly influence breast cancer risk, underpinning the benefit of endocrine therapies in the treatment of breast cancer. Modulation of their effects through ovarian ablation or chemoprevention strategies also significantly decreases breast cancer incidence. Conversely, there is an increased risk of breast cancer associated with pregnancy in the short term. The cellular mechanisms underlying these observations, however, are poorly defined. Here we demonstrate that mouse mammary stem cells (MaSCs) are highly responsive to steroid hormone signalling, despite lacking the oestrogen and progesterone receptors. Ovariectomy markedly diminished MaSC number and outgrowth potential in vivo, whereas MaSC activity increased in mice treated with oestrogen plus progesterone. Notably, even three weeks of treatment with the aromatase inhibitor letrozole was sufficient to reduce the MaSC pool. In contrast, pregnancy led to a transient 11-fold increase in MaSC numbers, probably mediated through paracrine signalling from RANK ligand. The augmented MaSC pool indicates a cellular basis for the short-term increase in breast cancer incidence that accompanies pregnancy. These findings further indicate that breast cancer chemoprevention may be achieved, in part, through suppression of MaSC function.