Epithelial zonation along the mouse and human small intestine defines five discrete metabolic domains.

A key aspect of nutrient absorption is the exquisite division of labour across the length of the small intestine, with individual nutrients taken up at different proximal:distal positions. For millennia, the small intestine was thought to comprise three segments with indefinite borders: the duodenum, jejunum and ileum. By examining the fine-scale longitudinal transcriptional patterns that span the mouse and human small intestine, we instead identified five domains of nutrient absorption that mount distinct responses to dietary changes, and three regional stem cell populations. Molecular domain identity can be detected with machine learning, which provides a systematic method to computationally identify intestinal domains in mice. We generated a predictive model of transcriptional control of domain identity and validated the roles of Ppar-δ and Cdx1 in patterning lipid metabolism-associated genes. These findings represent a foundational framework for the zonation of absorption across the mammalian small intestine.

Epithelial zonation along the mouse and human small intestine defines five discrete metabolic domains.

A key aspect of nutrient absorption is the exquisite division of labor across the length of the small intestine, with individual classes of micronutrients taken up at different positions. For millennia, the small intestine was thought to comprise three segments with indefinite borders: the duodenum, jejunum, and ileum. By examining fine-scale longitudinal segmentation of the mouse and human small intestines, we identified transcriptional signatures and upstream regulatory factors that define five domains of nutrient absorption, distinct from the three traditional sections. Spatially restricted expression programs were most prominent in nutrient-absorbing enterocytes but initially arose in intestinal stem cells residing in three regional populations. While a core signature was maintained across mice and humans with different diets and environments, domain properties were influenced by dietary changes. We established the functions of Ppar-ẟ and Cdx1 in patterning lipid metabolism in distal domains and generated a predictive model of additional transcription factors that direct domain identity. Molecular domain identity can be detected with machine learning, representing the first systematic method to computationally identify specific intestinal regions in mice. These findings provide a foundational framework for the identity and control of longitudinal zonation of absorption along the proximal:distal small intestinal axis.

Patterning and folding of intestinal villi by active mesenchymal dewetting.

Tissue folding generates structural motifs critical to organ function. In the intestine, bending of a flat epithelium into a periodic pattern of folds gives rise to villi, the numerous finger-like protrusions that are essential for nutrient absorption. However, the molecular and mechanical mechanisms driving the initiation and morphogenesis of villi remain a matter of debate. Here, we identify an active mechanical mechanism that simultaneously patterns and folds intestinal villi. We find that PDGFRA+ subepithelial mesenchymal cells generate myosin II-dependent forces sufficient to produce patterned curvature in neighboring tissue interfaces. At the cell-level, this occurs through a process dependent upon matrix metalloproteinase-mediated tissue fluidization and altered cell-ECM adhesion. By combining computational models with in vivo experiments, we reveal these cellular features manifest at the tissue-level as differences in interfacial tensions that promote mesenchymal aggregation and interface bending through a process analogous to the active de-wetting of a thin liquid film.

Lymphangiocrine signals are required for proper intestinal repair after cytotoxic injury.

The intestinal epithelium undergoes continuous renewal and has an exceptional capacity to regenerate after injury. Maintenance and proliferation of intestinal stem cells (ISCs) are regulated by their surrounding niche, largely through Wnt signaling. However, it remains unclear which niche cells produce signals during different injury states, and the role of endothelial cells (ECs) as a component of the ISC niche during homeostasis and after injury has been underappreciated. Here, we show that lymphatic endothelial cells (LECs) reside in proximity to crypt epithelial cells and secrete molecules that support epithelial renewal and repair. LECs are an essential source of Wnt signaling in the small intestine, as loss of LEC-derived Rspo3 leads to a lower number of stem and progenitor cells and hinders recovery after cytotoxic injury. Together, our findings identify LECs as an essential niche component for optimal intestinal recovery after cytotoxic injury.

Whole-exome sequencing reveals damaging gene variants associated with hypoalphalipoproteinemia.

Low levels of high density lipoprotein-cholesterol (HDL-C) are associated with an elevated risk of arteriosclerotic coronary heart disease. Heritability of HDL-C levels is high. In this research discovery study, we used whole-exome sequencing to identify damaging gene variants that may play significant roles in determining HDL-C levels. We studied 204 individuals with a mean HDL-C level of 27.8 ± 6.4 mg/dl (range: 4-36 mg/dl). Data were analyzed by statistical gene burden testing and by filtering against candidate gene lists. We found 120 occurrences of probably damaging variants (116 heterozygous; four homozygous) among 45 of 104 recognized HDL candidate genes. Those with the highest prevalence of damaging variants were ABCA1 (n = 20), STAB1 (n = 9), OSBPL1A (n = 8), CPS1 (n = 8), CD36 (n = 7), LRP1 (n = 6), ABCA8 (n = 6), GOT2 (n = 5), AMPD3 (n = 5), WWOX (n = 4), and IRS1 (n = 4). Binomial analysis for damaging missense or loss-of-function variants identified the ABCA1 and LDLR genes at genome-wide significance. In conclusion, whole-exome sequencing of individuals with low HDL-C showed the burden of damaging rare variants in the ABCA1 and LDLR genes is particularly high and revealed numerous occurrences in HDL candidate genes, including many genes identified in genome-wide association study reports. Many of these genes are involved in cancer biology, which accords with epidemiologic findings of the association of HDL deficiency with increased risk of cancer, thus presenting a new area of interest in HDL genomics.

The role of exome sequencing in newborn screening for inborn errors of metabolism.

Public health newborn screening (NBS) programs provide population-scale ascertainment of rare, treatable conditions that require urgent intervention. Tandem mass spectrometry (MS/MS) is currently used to screen newborns for a panel of rare inborn errors of metabolism (IEMs)1-4. The NBSeq project evaluated whole-exome sequencing (WES) as an innovative methodology for NBS. We obtained archived residual dried blood spots and data for nearly all IEM cases from the 4.5 million infants born in California between mid-2005 and 2013 and from some infants who screened positive by MS/MS, but were unaffected upon follow-up testing. WES had an overall sensitivity of 88% and specificity of 98.4%, compared to 99.0% and 99.8%, respectively for MS/MS, although effectiveness varied among individual IEMs. Thus, WES alone was insufficiently sensitive or specific to be a primary screen for most NBS IEMs. However, as a secondary test for infants with abnormal MS/MS screens, WES could reduce false-positive results, facilitate timely case resolution and in some instances even suggest more appropriate or specific diagnosis than that initially obtained. This study represents the largest, to date, sequencing effort of an entire population of IEM-affected cases, allowing unbiased assessment of current capabilities of WES as a tool for population screening.

Three patients with homozygous familial hypercholesterolemia: Genomic sequencing and kindred analysis.

BACKGROUND: Homozygous Familial Hypercholesterolemia (HoFH) is an inherited recessive condition associated with extremely high levels of low-density lipoprotein (LDL) cholesterol in affected individuals. It is usually caused by homozygous or compound heterozygous functional mutations in the LDL receptor (LDLR). A number of mutations causing FH have been reported in literature and such genetic heterogeneity presents great challenges for disease diagnosis. OBJECTIVE: We aim to determine the likely genetic defects responsible for three cases of pediatric HoFH in two kindreds. METHODS: We applied whole exome sequencing (WES) on the two probands to determine the likely functional variants among candidate FH genes. We additionally applied 10x Genomics (10xG) Linked-Reads whole genome sequencing (WGS) on one of the kindreds to identify potentially deleterious structural variants (SVs) underlying HoFH. A PCR-based screening assay was also established to detect the LDLR structural variant in a cohort of 641 patients with elevated LDL. RESULTS: In the Caucasian kindred, the FH homozygosity can be attributed to two compound heterozygous LDLR damaging variants, an exon 12 p.G592E missense mutation and a novel 3kb exon 1 deletion. By analyzing the 10xG phased data, we ascertained that this deletion allele was most likely to have originated from a Russian ancestor. In the Mexican kindred, the strikingly elevated LDL cholesterol level can be attributed to a homozygous frameshift LDLR variant p.E113fs. CONCLUSIONS: While the application of WES can provide a cost-effective way of identifying the genetic causes of FH, it often lacks sensitivity for detecting structural variants. Our finding of the LDLR exon 1 deletion highlights the broader utility of Linked-Read WGS in detecting SVs in the clinical setting, especially when HoFH patients remain undiagnosed after WES.

Tissue signals imprint ILC2 identity with anticipatory function.

Group 2 innate lymphoid cells (ILC2s) are distributed systemically and produce type 2 cytokines in response to a variety of stimuli, including the epithelial cytokines interleukin (IL)-25, IL-33, and thymic stromal lymphopoietin (TSLP). Transcriptional profiling of ILC2s from different tissues, however, grouped ILC2s according to their tissue of origin, even in the setting of combined IL-25-, IL-33-receptor-, and TSLP-receptor-deficiency. Single-cell profiling confirmed a tissue-organizing transcriptome and identified ILC2 subsets expressing distinct activating receptors, including the major subset of skin ILC2s, which were activated preferentially by IL-18. Tissue ILC2 subsets were unaltered in number and expression in germ-free mice, suggesting that endogenous, tissue-derived signals drive the maturation of ILC2 subsets by controlling expression of distinct patterns of activating receptors, thus anticipating tissue-specific perturbations occurring later in life.

Whole exome and whole genome sequencing with dried blood spot DNA without whole genome amplification.

Newborn screening (NBS) for rare conditions is performed in all 50 states in the USA. We have partnered with the California Department of Public Health Genetic Disease Laboratory to determine whether sufficient DNA can be extracted from archived dried blood spots (DBS) for next-generation sequencing in the hopes that next-generation sequencing can play a role in NBS. We optimized the DNA extraction and sequencing library preparation protocols for residual infant DBS archived over 20 years ago and successfully obtained acceptable whole exome and whole genome sequencing data. This sequencing study using DBS DNA without whole genome amplification prior to sequencing library preparation provides evidence that properly stored residual newborn DBS are a satisfactory source of DNA for genetic studies.

BLIMP1 Induces Transient Metastatic Heterogeneity in Pancreatic Cancer.

Pancreatic ductal adenocarcinoma (PDAC) is one of the most metastatic and deadly cancers. Despite the clinical significance of metastatic spread, our understanding of molecular mechanisms that drive PDAC metastatic ability remains limited. By generating a genetically engineered mouse model of human PDAC, we uncover a transient subpopulation of cancer cells with exceptionally high metastatic ability. Global gene expression profiling and functional analyses uncovered the transcription factor BLIMP1 as a driver of PDAC metastasis. The highly metastatic PDAC subpopulation is enriched for hypoxia-induced genes, and hypoxia-mediated induction of BLIMP1 contributes to the regulation of a subset of hypoxia-associated gene expression programs. These findings support a model in which upregulation of BLIMP1 links microenvironmental cues to a metastatic stem cell character.Significance: PDAC is an almost uniformly lethal cancer, largely due to its tendency for metastasis. We define a highly metastatic subpopulation of cancer cells, uncover a key transcriptional regulator of metastatic ability, and define hypoxia as an important factor within the tumor microenvironment that increases metastatic proclivity. Cancer Discov; 7(10); 1184-99. ©2017 AACR.See related commentary by Vakoc and Tuveson, p. 1067This article is highlighted in the In This Issue feature, p. 1047.

Identification and Targeting of Long-Term Tumor-Propagating Cells in Small Cell Lung Cancer.

Small cell lung cancer (SCLC) is a neuroendocrine lung cancer characterized by fast growth, early dissemination, and rapid resistance to chemotherapy. We identified a population of long-term tumor-propagating cells (TPCs) in a mouse model of SCLC. This population, marked by high levels of EpCAM and CD24, is also prevalent in human primary SCLC tumors. Murine SCLC TPCs are numerous and highly proliferative but not intrinsically chemoresistant, indicating that not all clinical features of SCLC are linked to TPCs. SCLC TPCs possess a distinct transcriptional profile compared to non-TPCs, including elevated MYC activity. Genetic and pharmacological inhibition of MYC in SCLC cells to non-TPC levels inhibits long-term propagation but not short-term growth. These studies identify a highly tumorigenic population of SCLC cells in mouse models, cell lines, and patient tumors and a means to target them in this most fatal form of lung cancer.

Medulloblastoma-associated DDX3 variant selectively alters the translational response to stress.

DDX3X encodes a DEAD-box family RNA helicase (DDX3) commonly mutated in medulloblastoma, a highly aggressive cerebellar tumor affecting both children and adults. Despite being implicated in several facets of RNA metabolism, the nature and scope of DDX3's interactions with RNA remain unclear. Here, we show DDX3 collaborates extensively with the translation initiation machinery through direct binding to 5'UTRs of nearly all coding RNAs, specific sites on the 18S rRNA, and multiple components of the translation initiation complex. Impairment of translation initiation is also evident in primary medulloblastomas harboring mutations in DDX3X, further highlighting DDX3's role in this process. Arsenite-induced stress shifts DDX3 binding from the 5'UTR into the coding region of mRNAs concomitant with a general reduction of translation, and both the shift of DDX3 on mRNA and decreased translation are blunted by expression of a catalytically-impaired, medulloblastoma-associated DDX3R534H variant. Furthermore, despite the global repression of translation induced by arsenite, translation is preserved on select genes involved in chromatin organization in DDX3R534H-expressing cells. Thus, DDX3 interacts extensively with RNA and ribosomal machinery to help remodel the translation landscape in response to stress, while cancer-related DDX3 variants adapt this response to selectively preserve translation.

Comprehensive genomic profiles of small cell lung cancer.

We have sequenced the genomes of 110 small cell lung cancers (SCLC), one of the deadliest human cancers. In nearly all the tumours analysed we found bi-allelic inactivation of TP53 and RB1, sometimes by complex genomic rearrangements. Two tumours with wild-type RB1 had evidence of chromothripsis leading to overexpression of cyclin D1 (encoded by the CCND1 gene), revealing an alternative mechanism of Rb1 deregulation. Thus, loss of the tumour suppressors TP53 and RB1 is obligatory in SCLC. We discovered somatic genomic rearrangements of TP73 that create an oncogenic version of this gene, TP73Δex2/3. In rare cases, SCLC tumours exhibited kinase gene mutations, providing a possible therapeutic opportunity for individual patients. Finally, we observed inactivating mutations in NOTCH family genes in 25% of human SCLC. Accordingly, activation of Notch signalling in a pre-clinical SCLC mouse model strikingly reduced the number of tumours and extended the survival of the mutant mice. Furthermore, neuroendocrine gene expression was abrogated by Notch activity in SCLC cells. This first comprehensive study of somatic genome alterations in SCLC uncovers several key biological processes and identifies candidate therapeutic targets in this highly lethal form of cancer.

Inhibition of pluripotency networks by the Rb tumor suppressor restricts reprogramming and tumorigenesis.

Mutations in the retinoblastoma tumor suppressor gene Rb are involved in many forms of human cancer. In this study, we investigated the early consequences of inactivating Rb in the context of cellular reprogramming. We found that Rb inactivation promotes the reprogramming of differentiated cells to a pluripotent state. Unexpectedly, this effect is cell cycle independent, and instead reflects direct binding of Rb to pluripotency genes, including Sox2 and Oct4, which leads to a repressed chromatin state. More broadly, this regulation of pluripotency networks and Sox2 in particular is critical for the initiation of tumors upon loss of Rb in mice. These studies therefore identify Rb as a global transcriptional repressor of pluripotency networks, providing a molecular basis for previous reports about its involvement in cell fate pliability, and implicate misregulation of pluripotency factors such as Sox2 in tumorigenesis related to loss of Rb function.

Long noncoding RNA EWSAT1-mediated gene repression facilitates Ewing sarcoma oncogenesis.

Chromosomal translocation that results in fusion of the genes encoding RNA-binding protein EWS and transcription factor FLI1 (EWS-FLI1) is pathognomonic for Ewing sarcoma. EWS-FLI1 alters gene expression through mechanisms that are not completely understood. We performed RNA sequencing (RNAseq) analysis on primary pediatric human mesenchymal progenitor cells (pMPCs) expressing EWS-FLI1 in order to identify gene targets of this oncoprotein. We determined that long noncoding RNA-277 (Ewing sarcoma-associated transcript 1 [EWSAT1]) is upregulated by EWS-FLI1 in pMPCs. Inhibition of EWSAT1 expression diminished the ability of Ewing sarcoma cell lines to proliferate and form colonies in soft agar, whereas EWSAT1 inhibition had no effect on other cell types tested. Expression of EWS-FLI1 and EWSAT1 repressed gene expression, and a substantial fraction of targets that were repressed by EWS-FLI1 were also repressed by EWSAT1. Analysis of RNAseq data from primary human Ewing sarcoma further supported a role for EWSAT1 in mediating gene repression. We identified heterogeneous nuclear ribonucleoprotein (HNRNPK) as an RNA-binding protein that interacts with EWSAT1 and found a marked overlap in HNRNPK-repressed genes and those repressed by EWS-FLI1 and EWSAT1, suggesting that HNRNPK participates in EWSAT1-mediated gene repression. Together, our data reveal that EWSAT1 is a downstream target of EWS-FLI1 that facilitates the development of Ewing sarcoma via the repression of target genes.

Organ size control is dominant over Rb family inactivation to restrict proliferation in vivo.

In mammals, a cell's decision to divide is thought to be under the control of the Rb/E2F pathway. We previously found that inactivation of the Rb family of cell cycle inhibitors (Rb, p107, and p130) in quiescent liver progenitors leads to uncontrolled division and cancer initiation. Here, we show that, in contrast, deletion of the entire Rb gene family in mature hepatocytes is not sufficient for their long-term proliferation. The cell cycle block in Rb family mutant hepatocytes is independent of the Arf/p53/p21 checkpoint but can be abrogated upon decreasing liver size. At the molecular level, we identify YAP, a transcriptional regulator involved in organ size control, as a factor required for the sustained expression of cell cycle genes in hepatocytes. These experiments identify a higher level of regulation of the cell cycle in vivo in which signals regulating organ size are dominant regulators of the core cell cycle machinery.

A meta-analysis of lung cancer gene expression identifies PTK7 as a survival gene in lung adenocarcinoma.

Lung cancer remains the most common cause of cancer-related death worldwide and it continues to lack effective treatment. The increasingly large and diverse public databases of lung cancer gene expression constitute a rich source of candidate oncogenic drivers and therapeutic targets. To define novel targets for lung adenocarcinoma, we conducted a large-scale meta-analysis of genes specifically overexpressed in adenocarcinoma. We identified an 11-gene signature that was overexpressed consistently in adenocarcinoma specimens relative to normal lung tissue. Six genes in this signature were specifically overexpressed in adenocarcinoma relative to other subtypes of non-small cell lung cancer (NSCLC). Among these genes was the little studied protein tyrosine kinase PTK7. Immunohistochemical analysis confirmed that PTK7 is highly expressed in primary adenocarcinoma patient samples. RNA interference-mediated attenuation of PTK7 decreased cell viability and increased apoptosis in a subset of adenocarcinoma cell lines. Further, loss of PTK7 activated the MKK7-JNK stress response pathway and impaired tumor growth in xenotransplantation assays. Our work defines PTK7 as a highly and specifically expressed gene in adenocarcinoma and a potential therapeutic target in this subset of NSCLC.

A drug repositioning approach identifies tricyclic antidepressants as inhibitors of small cell lung cancer and other neuroendocrine tumors.

UNLABELLED: Small cell lung cancer (SCLC) is an aggressive neuroendocrine subtype of lung cancer with high mortality. We used a systematic drug repositioning bioinformatics approach querying a large compendium of gene expression profiles to identify candidate U.S. Food and Drug Administration (FDA)-approved drugs to treat SCLC. We found that tricyclic antidepressants and related molecules potently induce apoptosis in both chemonaïve and chemoresistant SCLC cells in culture, in mouse and human SCLC tumors transplanted into immunocompromised mice, and in endogenous tumors from a mouse model for human SCLC. The candidate drugs activate stress pathways and induce cell death in SCLC cells, at least in part by disrupting autocrine survival signals involving neurotransmitters and their G protein-coupled receptors. The candidate drugs inhibit the growth of other neuroendocrine tumors, including pancreatic neuroendocrine tumors and Merkel cell carcinoma. These experiments identify novel targeted strategies that can be rapidly evaluated in patients with neuroendocrine tumors through the repurposing of approved drugs. SIGNIFICANCE: Our work shows the power of bioinformatics-based drug approaches to rapidly repurpose FDA-approved drugs and identifies a novel class of molecules to treat patients with SCLC, a cancer for which no effective novel systemic treatments have been identified in several decades. In addition, our experiments highlight the importance of novel autocrine mechanisms in promoting the growth of neuroendocrine tumor cells.

A rare population of CD24(+)ITGB4(+)Notch(hi) cells drives tumor propagation in NSCLC and requires Notch3 for self-renewal.

Sustained tumor progression has been attributed to a distinct population of tumor-propagating cells (TPCs). To identify TPCs relevant to lung cancer pathogenesis, we investigated functional heterogeneity in tumor cells isolated from Kras-driven mouse models of non-small-cell lung cancer (NSCLC). CD24(+)ITGB4(+)Notch(hi) cells are capable of propagating tumor growth in both a clonogenic and an orthotopic serial transplantation assay. While all four Notch receptors mark TPCs, Notch3 plays a nonredundant role in tumor cell propagation in two mouse models and in human NSCLC. The TPC population is enriched after chemotherapy, and the gene signature of mouse TPCs correlates with poor prognosis in human NSCLC. The role of Notch3 in tumor propagation may provide a therapeutic target for NSCLC.

Cross-species functional analysis of cancer-associated fibroblasts identifies a critical role for CLCF1 and IL-6 in non-small cell lung cancer in vivo.

Cancer-associated fibroblasts (CAF) have been reported to support tumor progression by a variety of mechanisms. However, their role in the progression of non-small cell lung cancer (NSCLC) remains poorly defined. In addition, the extent to which specific proteins secreted by CAFs contribute directly to tumor growth is unclear. To study the role of CAFs in NSCLCs, a cross-species functional characterization of mouse and human lung CAFs was conducted. CAFs supported the growth of lung cancer cells in vivo by secretion of soluble factors that directly stimulate the growth of tumor cells. Gene expression analysis comparing normal mouse lung fibroblasts and mouse lung CAFs identified multiple genes that correlate with the CAF phenotype. A gene signature of secreted genes upregulated in CAFs was an independent marker of poor survival in patients with NSCLC. This secreted gene signature was upregulated in normal lung fibroblasts after long-term exposure to tumor cells, showing that lung fibroblasts are "educated" by tumor cells to acquire a CAF-like phenotype. Functional studies identified important roles for CLCF1-CNTFR and interleukin (IL)-6-IL-6R signaling in promoting growth of NSCLCs. This study identifies novel soluble factors contributing to the CAF protumorigenic phenotype in NSCLCs and suggests new avenues for the development of therapeutic strategies.