Defining the Design Principles of Skin Epidermis Postnatal Growth.

During embryonic and postnatal development, organs and tissues grow steadily to achieve their final size at the end of puberty. However, little is known about the cellular dynamics that mediate postnatal growth. By combining in vivo clonal lineage tracing, proliferation kinetics, single-cell transcriptomics, and in vitro micro-pattern experiments, we resolved the cellular dynamics taking place during postnatal skin epidermis expansion. Our data revealed that harmonious growth is engineered by a single population of developmental progenitors presenting a fixed fate imbalance of self-renewing divisions with an ever-decreasing proliferation rate. Single-cell RNA sequencing revealed that epidermal developmental progenitors form a more uniform population compared with adult stem and progenitor cells. Finally, we found that the spatial pattern of cell division orientation is dictated locally by the underlying collagen fiber orientation. Our results uncover a simple design principle of organ growth where progenitors and differentiated cells expand in harmony with their surrounding tissues.

RHOJ controls EMT-associated resistance to chemotherapy.

The resistance of cancer cells to therapy is responsible for the death of most patients with cancer1. Epithelial-to-mesenchymal transition (EMT) has been associated with resistance to therapy in different cancer cells2,3. However, the mechanisms by which EMT mediates resistance to therapy remain poorly understood. Here, using a mouse model of skin squamous cell carcinoma undergoing spontaneous EMT during tumorigenesis, we found that EMT tumour cells are highly resistant to a wide range of anti-cancer therapies both in vivo and in vitro. Using gain and loss of function studies in vitro and in vivo, we found that RHOJ-a small GTPase that is preferentially expressed in EMT cancer cells-controls resistance to therapy. Using genome-wide transcriptomic and proteomic profiling, we found that RHOJ regulates EMT-associated resistance to chemotherapy by enhancing the response to replicative stress and activating the DNA-damage response, enabling tumour cells to rapidly repair DNA lesions induced by chemotherapy. RHOJ interacts with proteins that regulate nuclear actin, and inhibition of actin polymerization sensitizes EMT tumour cells to chemotherapy-induced cell death in a RHOJ-dependent manner. Together, our study uncovers the role and the mechanisms through which RHOJ acts as a key regulator of EMT-associated resistance to chemotherapy.

Pharmacological targeting of netrin-1 inhibits EMT in cancer.

Epithelial-to-mesenchymal transition (EMT) regulates tumour initiation, progression, metastasis and resistance to anti-cancer therapy1-7. Although great progress has been made in understanding the role of EMT and its regulatory mechanisms in cancer, no therapeutic strategy to pharmacologically target EMT has been identified. Here we found that netrin-1 is upregulated in a primary mouse model of skin squamous cell carcinoma (SCC) exhibiting spontaneous EMT. Pharmacological inhibition of netrin-1 by administration of NP137, a netrin-1-blocking monoclonal antibody currently used in clinical trials in human cancer (ClinicalTrials.gov identifier NCT02977195 ), decreased the proportion of EMT tumour cells in skin SCC, decreased the number of metastases and increased the sensitivity of tumour cells to chemotherapy. Single-cell RNA sequencing revealed the presence of different EMT states, including epithelial, early and late hybrid EMT, and full EMT states, in control SCC. By contrast, administration of NP137 prevented the progression of cancer cells towards a late EMT state and sustained tumour epithelial states. Short hairpin RNA knockdown of netrin-1 and its receptor UNC5B in EPCAM+ tumour cells inhibited EMT in vitro in the absence of stromal cells and regulated a common gene signature that promotes tumour epithelial state and restricts EMT. To assess the relevance of these findings to human cancers, we treated mice transplanted with the A549 human cancer cell line-which undergoes EMT following TGFß1 administration8,9-with NP137. Netrin-1 inhibition decreased EMT in these transplanted A549 cells. Together, our results identify a pharmacological strategy for targeting EMT in cancer, opening up novel therapeutic interventions for anti-cancer therapy.

Fat1 deletion promotes hybrid EMT state, tumour stemness and metastasis.

FAT1, which encodes a protocadherin, is one of the most frequently mutated genes in human cancers1-5. However, the role and the molecular mechanisms by which FAT1 mutations control tumour initiation and progression are poorly understood. Here, using mouse models of skin squamous cell carcinoma and lung tumours, we found that deletion of Fat1 accelerates tumour initiation and malignant progression and promotes a hybrid epithelial-to-mesenchymal transition (EMT) phenotype. We also found this hybrid EMT state in FAT1-mutated human squamous cell carcinomas. Skin squamous cell carcinomas in which Fat1 was deleted presented increased tumour stemness and spontaneous metastasis. We performed transcriptional and chromatin profiling combined with proteomic analyses and mechanistic studies, which revealed that loss of function of FAT1 activates a CAMK2-CD44-SRC axis that promotes YAP1 nuclear translocation and ZEB1 expression that stimulates the mesenchymal state. This loss of function also inactivates EZH2, promoting SOX2 expression, which sustains the epithelial state. Our comprehensive analysis identified drug resistance and vulnerabilities in FAT1-deficient tumours, which have important implications for cancer therapy. Our studies reveal that, in mouse and human squamous cell carcinoma, loss of function of FAT1 promotes tumour initiation, progression, invasiveness, stemness and metastasis through the induction of a hybrid EMT state.

Mechanisms of stretch-mediated skin expansion at single-cell resolution.

The ability of the skin to grow in response to stretching has been exploited in reconstructive surgery1. Although the response of epidermal cells to stretching has been studied in vitro2,3, it remains unclear how mechanical forces affect their behaviour in vivo. Here we develop a mouse model in which the consequences of stretching on skin epidermis can be studied at single-cell resolution. Using a multidisciplinary approach that combines clonal analysis with quantitative modelling and single-cell RNA sequencing, we show that stretching induces skin expansion by creating a transient bias in the renewal activity of epidermal stem cells, while a second subpopulation of basal progenitors remains committed to differentiation. Transcriptional and chromatin profiling identifies how cell states and gene-regulatory networks are modulated by stretching. Using pharmacological inhibitors and mouse mutants, we define the step-by-step mechanisms that control stretch-mediated tissue expansion at single-cell resolution in vivo.

Heterotypic cell-cell communication regulates glandular stem cell multipotency.

Glandular epithelia, including the mammary and prostate glands, are composed of basal cells (BCs) and luminal cells (LCs)1,2. Many glandular epithelia develop from multipotent basal stem cells (BSCs) that are replaced in adult life by distinct pools of unipotent stem cells1,3-8. However, adult unipotent BSCs can reactivate multipotency under regenerative conditions and upon oncogene expression3,9-13. This suggests that an active mechanism restricts BSC multipotency under normal physiological conditions, although the nature of this mechanism is unknown. Here we show that the ablation of LCs reactivates the multipotency of BSCs from multiple epithelia both in vivo in mice and in vitro in organoids. Bulk and single-cell RNA sequencing revealed that, after LC ablation, BSCs activate a hybrid basal and luminal cell differentiation program before giving rise to LCs-reminiscent of the genetic program that regulates multipotency during embryonic development7. By predicting ligand-receptor pairs from single-cell data14, we find that TNF-which is secreted by LCs-restricts BC multipotency under normal physiological conditions. By contrast, the Notch, Wnt and EGFR pathways were activated in BSCs and their progeny after LC ablation; blocking these pathways, or stimulating the TNF pathway, inhibited regeneration-induced BC multipotency. Our study demonstrates that heterotypic communication between LCs and BCs is essential to maintain lineage fidelity in glandular epithelial stem cells.

A slow-cycling LGR5 tumour population mediates basal cell carcinoma relapse after therapy.

Basal cell carcinoma (BCC) is the most frequent cancer in humans and results from constitutive activation of the Hedgehog pathway1. Several Smoothened inhibitors are used to treat Hedgehog-mediated malignancies, including BCC and medulloblastoma2. Vismodegib, a Smoothened inhibitor, leads to BCC shrinkage in the majority of patients with BCC3, but the mechanism by which it mediates BCC regression is unknown. Here we used two genetically engineered mouse models of BCC4 to investigate the mechanisms by which inhibition of Smoothened mediates tumour regression. We found that vismodegib mediates BCC regression by inhibiting a hair follicle-like fate and promoting the differentiation of tumour cells. However, a small population of tumour cells persists and is responsible for tumour relapse following treatment discontinuation, mimicking the situation found in humans5. In both mouse and human BCC, this persisting, slow-cycling tumour population expresses LGR5 and is characterized by active Wnt signalling. Combining Lgr5 lineage ablation or inhibition of Wnt signalling with vismodegib treatment leads to eradication of BCC. Our results show that vismodegib induces tumour regression by promoting tumour differentiation, and demonstrates that the synergy between Wnt and Smoothened inhibitors is a clinically relevant strategy for overcoming tumour relapse in BCC.

Identification of the tumour transition states occurring during EMT.

In cancer, the epithelial-to-mesenchymal transition (EMT) is associated with tumour stemness, metastasis and resistance to therapy. It has recently been proposed that, rather than being a binary process, EMT occurs through distinct intermediate states. However, there is no direct in vivo evidence for this idea. Here we screen a large panel of cell surface markers in skin and mammary primary tumours, and identify the existence of multiple tumour subpopulations associated with different EMT stages: from epithelial to completely mesenchymal states, passing through intermediate hybrid states. Although all EMT subpopulations presented similar tumour-propagating cell capacity, they displayed differences in cellular plasticity, invasiveness and metastatic potential. Their transcriptional and epigenetic landscapes identify the underlying gene regulatory networks, transcription factors and signalling pathways that control these different EMT transition states. Finally, these tumour subpopulations are localized in different niches that differentially regulate EMT transition states.

Dealing with Macrophage Plasticity to Address Therapeutic Challenges in Head and Neck Cancers.

The head and neck tumor microenvironment (TME) is highly infiltrated with macrophages. More specifically, tumor-associated macrophages (TAM/M2-like) are one of the most critical components associated with poor overall survival in head and neck cancers (HNC). Two extreme states of macrophage phenotypes are described as conducting pro-inflammatory/anti-tumoral (M1) or anti-inflammatory/pro-tumoral (M2) activities. Moreover, specific metabolic pathways as well as oxidative stress responses are tightly associated with their phenotypes and functions. Hence, due to their plasticity, targeting M2 macrophages to repolarize in the M1 phenotype would be a promising cancer treatment. In this context, we evaluated macrophage infiltration in 60 HNC patients and demonstrated the high infiltration of CD68+ cells that were mainly related to CD163+ M2 macrophages. We then optimized a polarization protocol from THP1 monocytes, validated by specific gene and protein expression levels. In addition, specific actors of glutamine pathway and oxidative stress were quantified to indicate the use of glutaminolysis by M2 and the production of reactive oxygen species by M1. Finally, we evaluated and confirmed the plasticity of our model using M1 activators to repolarize M2 in M1. Overall, our study provides a complete reversible polarization protocol allowing us to further evaluate various reprogramming effectors targeting glutaminolysis and/or oxidative stress in macrophages.

Reactivation of multipotency by oncogenic PIK3CA induces breast tumour heterogeneity.

Breast cancer is the most frequent cancer in women and consists of heterogeneous types of tumours that are classified into different histological and molecular subtypes. PIK3CA and P53 (also known as TP53) are the two most frequently mutated genes and are associated with different types of human breast cancers. The cellular origin and the mechanisms leading to PIK3CA-induced tumour heterogeneity remain unknown. Here we used a genetic approach in mice to define the cellular origin of Pik3ca-derived tumours and the impact of mutations in this gene on tumour heterogeneity. Surprisingly, oncogenic Pik3ca(H1047R) mutant expression at physiological levels in basal cells using keratin (K)5-CreER(T2) mice induced the formation of luminal oestrogen receptor (ER)-positive/progesterone receptor (PR)-positive tumours, while its expression in luminal cells using K8-CReER(T2) mice gave rise to luminal ER(+)PR(+) tumours or basal-like ER(-)PR(-) tumours. Concomitant deletion of p53 and expression of Pik3ca(H1047R) accelerated tumour development and induced more aggressive mammary tumours. Interestingly, expression of Pik3ca(H1047R) in unipotent basal cells gave rise to luminal-like cells, while its expression in unipotent luminal cells gave rise to basal-like cells before progressing into invasive tumours. Transcriptional profiling of cells that underwent cell fate transition upon Pik3ca(H1047R) expression in unipotent progenitors demonstrated a profound oncogene-induced reprogramming of these newly formed cells and identified gene signatures characteristic of the different cell fate switches that occur upon Pik3ca(H1047R) expression in basal and luminal cells, which correlated with the cell of origin, tumour type and different clinical outcomes. Altogether our study identifies the cellular origin of Pik3ca-induced tumours and reveals that oncogenic Pik3ca(H1047R) activates a multipotent genetic program in normally lineage-restricted populations at the early stage of tumour initiation, setting the stage for future intratumoural heterogeneity. These results have important implications for our understanding of the mechanisms controlling tumour heterogeneity and the development of new strategies to block PIK3CA breast cancer initiation.

SOX2 controls tumour initiation and cancer stem-cell functions in squamous-cell carcinoma.

Cancer stem cells (CSCs) have been reported in various cancers, including in skin squamous-cell carcinoma (SCC). The molecular mechanisms regulating tumour initiation and stemness are still poorly characterized. Here we find that Sox2, a transcription factor expressed in various types of embryonic and adult stem cells, was the most upregulated transcription factor in the CSCs of squamous skin tumours in mice. SOX2 is absent in normal epidermis but begins to be expressed in the vast majority of mouse and human pre-neoplastic skin tumours, and continues to be expressed in a heterogeneous manner in invasive mouse and human SCCs. In contrast to other SCCs, in which SOX2 is frequently genetically amplified, the expression of SOX2 in mouse and human skin SCCs is transcriptionally regulated. Conditional deletion of Sox2 in the mouse epidermis markedly decreases skin tumour formation after chemical-induced carcinogenesis. Using green fluorescent protein (GFP) as a reporter of Sox2 transcriptional expression (SOX2-GFP knock-in mice), we showed that SOX2-expressing cells in invasive SCC are greatly enriched in tumour-propagating cells, which further increase upon serial transplantations. Lineage ablation of SOX2-expressing cells within primary benign and malignant SCCs leads to tumour regression, consistent with the critical role of SOX2-expressing cells in tumour maintenance. Conditional Sox2 deletion in pre-existing skin papilloma and SCC leads to tumour regression and decreases the ability of cancer cells to be propagated upon transplantation into immunodeficient mice, supporting the essential role of SOX2 in regulating CSC functions. Transcriptional profiling of SOX2-GFP-expressing CSCs and of tumour epithelial cells upon Sox2 deletion uncovered a gene network regulated by SOX2 in primary tumour cells in vivo. Chromatin immunoprecipitation identified several direct SOX2 target genes controlling tumour stemness, survival, proliferation, adhesion, invasion and paraneoplastic syndrome. We demonstrate that SOX2, by marking and regulating the functions of skin tumour-initiating cells and CSCs, establishes a continuum between tumour initiation and progression in primary skin tumours.

Skin squamous cell carcinoma propagating cells increase with tumour progression and invasiveness.

Cancer stem cells have been described in various cancers including squamous tumours of the skin by their ability to reform secondary tumours upon transplantation into immunodeficient mice. Here, we used transplantation of limiting dilution of different populations of FACS-isolated tumour cells from four distinct mouse models of squamous skin tumours to investigate the frequency of tumour propagating cells (TPCs) at different stages of tumour progression. We found that benign papillomas, despite growing rapidly in vivo and being clonogenic in vitro, reformed secondary tumours upon transplantation at very low frequency and only when tumour cells were co-transplanted together with tumour-associated fibroblasts or endothelial cells. In two models of skin squamous cell carcinoma (SCC), TPCs increased with tumour invasiveness. Interestingly, the frequency of TPCs increased in CD34(HI) but not in CD34(LO) SCC cells with serial transplantations, while the two populations initially gave rise to secondary tumours with the same frequency. Our results illustrate the progressive increase of squamous skin TPCs with tumour progression and invasiveness and reveal that serial transplantation may be required to define the long-term renewal potential of TPCs.

Cell type and stage specific transcriptional, chromatin and cell-cell communication landscapes in the mammary gland.

The mammary gland (MG) is composed of three main epithelial lineages, the basal cells (BC), the estrogen receptor (ER) positive luminal cells (ER+ LC), and the ER negative LC (ER- LC). Defining the cell identity of each lineage and how it is modulated throughout the different stages of life is important to understand how these cells function and communicate throughout life. Here, we used transgenic mice specifically labelling ER+ LC combined to cell surface markers to isolate with high purity the 3 distinct cell lineages of the mammary gland and defined their expression profiles and chromatin landscapes by performing bulk RNAseq and ATACseq of these isolated populations in puberty, adulthood and mid-pregnancy. Our analysis identified conserved genes, ligands and transcription factor (TF) associated with a specific lineage throughout life as well as genes, ligands and TFs specific for a particular stage of the MG. In summary, our study identified genes and TF network associated with the identity, function and cell-cell communication of the different epithelial lineages of the MG at different stages of life.

Mesp1 controls the chromatin and enhancer landscapes essential for spatiotemporal patterning of early cardiovascular progenitors.

The mammalian heart arises from various populations of Mesp1-expressing cardiovascular progenitors (CPs) that are specified during the early stages of gastrulation. Mesp1 is a transcription factor that acts as a master regulator of CP specification and differentiation. However, how Mesp1 regulates the chromatin landscape of nascent mesodermal cells to define the temporal and spatial patterning of the distinct populations of CPs remains unknown. Here, by combining ChIP-seq, RNA-seq and ATAC-seq during mouse pluripotent stem cell differentiation, we defined the dynamic remodelling of the chromatin landscape mediated by Mesp1. We identified different enhancers that are temporally regulated to erase the pluripotent state and specify the pools of CPs that mediate heart development. We identified Zic2 and Zic3 as essential cofactors that act with Mesp1 to regulate its transcription-factor activity at key mesodermal enhancers, thereby regulating the chromatin remodelling and gene expression associated with the specification of the different populations of CPs in vivo. Our study identifies the dynamics of the chromatin landscape and enhancer remodelling associated with temporal patterning of early mesodermal cells into the distinct populations of CPs that mediate heart development.

NR2F2 controls malignant squamous cell carcinoma state by promoting stemness and invasion and repressing differentiation.

The nongenetic mechanisms required to sustain malignant tumor state are poorly understood. During the transition from benign tumors to malignant carcinoma, tumor cells need to repress differentiation and acquire invasive features. Using transcriptional profiling of cancer stem cells from benign tumors and malignant skin squamous cell carcinoma (SCC), we identified the nuclear receptor NR2F2 as uniquely expressed in malignant SCC. Using genetic gain of function and loss of function in vivo, we show that NR2F2 is essential for promoting the malignant tumor state by controlling tumor stemness and maintenance in mouse and human SCC. We demonstrate that NR2F2 promotes tumor cell proliferation, epithelial-mesenchymal transition and invasive features, while repressing tumor differentiation and immune cell infiltration by regulating a common transcriptional program in mouse and human SCCs. Altogether, we identify NR2F2 as a key regulator of malignant cancer stem cell functions that promotes tumor renewal and restricts differentiation to sustain a malignant tumor state.

Nuclei Isolation from Whole Tissue using a Detergent and Enzyme-Free Method.

High-throughput transcriptome and epigenome profiling requires preparation of a single cell or single nuclei suspension. Preparation of the suspension with intact cell or nuclei involves dissociation and permeabilization, steps that can introduce unwanted noise and undesirable damage. Particularly, certain cell-types such as neurons are challenging to dissociate into individual cells. Additionally, permeabilization of the cellular membrane to release nuclei requires optimization by trial-and-error, which can be time consuming, labor intensive and financially nonviable. To enhance the robustness and reproducibility of sample preparation for high-throughput sequencing, we describe a rapid enzyme and detergent-free column-based nuclei isolation method. The protocol enables efficient isolation of nuclei from the entire zebrafish brain within 20 minutes. The isolated nuclei display intact nuclear morphology and low propensity to aggregate. Further, flow cytometry allows nuclei enrichment and clearance of cellular debris for downstream application. The protocol, which should work on soft tissues and cultured cells, provides a simple and accessible method for sample preparation that can be utilized for high-throughput profiling, simplifying the steps required for successful single-nuclei RNA-seq and ATAC-seq experiments.

Chemerin Added to Endothelin-1 Promotes Rat Pulmonary Artery Smooth Muscle Cell Proliferation and Migration.

BACKGROUND: While chemerin has been shown to increase proliferation and migration of systemic vascular smooth muscle cells (SMCs) contributing therefore to the development of hypertension, this remains to be clarified for the pulmonary circulation. METHODS: Expression of chemerin and its three receptors (CMKRL1, CCRL2, GPR1) was examined by immunohistochemistry and RTq-PCR in lungs, pulmonary artery, and thoracic aorta from Wistar rats. Primary cultured rat pulmonary artery and thoracic aorta SMCs treated with recombinant chemerin (tested from 5.10-9 to 10-7 mol/L) were assessed for proliferation and migration (both with 10-7 mol/L endothelin-1), as well as for staurosporine-induced apoptosis. RESULTS: In pulmonary artery and thoracic aorta, CMKLR1 expression was detected in both endothelial cells and SMCs. In primary cultured pulmonary artery SMCs, chemerin and its three receptors were expressed, and CMKLR1 expression was higher than those of CCRL2 and GPR1. Chemerin added to endothelin-1 increased pulmonary artery SMC proliferation, while chemerin or endothelin-1 alone did not. This effect was less pronounced in thoracic aorta SMCs. Chemerin induced pulmonary artery and thoracic aorta SMC migration, which was exacerbated by endothelin-1 and more pronounced in thoracic aorta SMCs. Chemerin concentration-dependently reduced staurosporine-induced apoptosis in both pulmonary artery and thoracic aorta SMCs. In pulmonary artery SMCs, endothelin-1 treatment increased the expression of CMKLR1, CCRL2, and GPR1, while these expressions were not altered in thoracic aorta SMCs. CONCLUSION: Chemerin/CMKRL1 signaling, in conjunction with a key mediator in the pathogenesis of pulmonary hypertensive diseases, endothelin-1, stimulated proliferation and migration, and increased resistance to apoptosis in rat primary cultured pulmonary artery SMCs. Our results suggest that this signaling could play a role in pulmonary artery remodeling observed in pulmonary hypertension.

Epidermal autonomous VEGFA/Flt1/Nrp1 functions mediate psoriasis-like disease.

Psoriasis is a common chronic skin disorder characterized by keratinocyte hyperproliferation with altered differentiation accompanied by inflammation and increased angiogenesis. It remains unclear whether the first events that initiate psoriasis development occur in keratinocytes or inflammatory cells. Here, using different psoriasis mouse models, we showed that conditional deletion of Flt1 or Nrp1 in epidermal cells inhibited psoriasis mediated by Vegfa overexpression or c-Jun/JunB deletion. Administration of anti-Nrp1 antibody reverted the psoriasis phenotype. Using transcriptional and chromatin profiling of epidermal cells following Vegfa overexpression together with Flt1 or Nrp1 deletion, we identified the gene regulatory network regulated by Vegfa/Nrp1/Flt1 during psoriasis development and uncovered a key role of Fosl1 in regulating the chromatin remodeling mediated by Vegfa overexpression in keratinocytes. In conclusion, our study identifies an epidermal autonomous function of Vegfa/Nrp1/Flt1 that mediates psoriatic-like disease and demonstrates the clinical relevance of blocking Vegfa/Nrp1/Flt1 axis in psoriasis.

Context Dependency of Epithelial-to-Mesenchymal Transition for Metastasis.

Epithelial-to-mesenchymal transition (EMT) has been proposed to be important for metastatic dissemination. However, recent studies have challenged the requirement of EMT for metastasis. Here, we assessed in different models of primary skin squamous cell carcinomas (SCCs) whether EMT is associated with metastasis. The incidence of metastasis was much higher in SCCs presenting EMT compared to SCCs without EMT, supporting the notion that a certain degree of EMT is required to initiate the metastatic cascade in primary skin SCCs. Most circulating tumor cells presented EMT, whereas most lung metastasis did not present EMT, showing that mesenchymal-to-epithelial transition is important for metastatic colonization. In contrast, immunodeficient mice transplanted with SCCs, whether displaying EMT or not, presented metastasis. Altogether, our data demonstrate that the association of EMT and metastasis is model dependent, and metastasis of primary skin SCCs is associated with EMT.

Publisher Correction: Early lineage segregation of multipotent embryonic mammary gland progenitors.

In the version of this Article originally published, ref. 52 was incorrectly only attributed to its corresponding author, Fre, S., and an older title was used. The correct citation should have been: Lilja, A. M. et al. Clonal analysis of Notch1-expressing cells reveals the existence of unipotent stem cells that retain long-term plasticity in the embryonic mammary gland. Nat. Cell Biol. https://doi.org/10.1038/s41556-018-0108-1 (2018)'. This has now been amended in all online versions of the Article.