Gene expression profile analysis of mouse colon embryonic development

Differential pre-malignant programs and microenvironment chart distinct paths to malignancy in human colorectal polyps

Colorectal cancers (CRCs) arise from precursor polyps whose cellular origins, molecular heterogeneity, and immunogenic potential may reveal diagnostic and therapeutic insights when analyzed at high resolution. We present a single-cell transcriptomic and imaging atlas of the two most common human colorectal polyps, conventional adenomas and serrated polyps, and their resulting CRC counterparts. Integrative analysis of 128 datasets from 62 participants reveals adenomas arise from WNT-driven expansion of stem cells, while serrated polyps derive from differentiated cells through gastric metaplasia. Metaplasia-associated damage is coupled to a cytotoxic immune microenvironment preceding hypermutation, driven partly by antigen-presentation differences associated with tumor cell-differentiation status. Microsatellite unstable CRCs contain distinct non-metaplastic regions where tumor cells acquire stem cell properties and cytotoxic immune cells are depleted. Our multi-omic atlas provides insights into malignant progression of colorectal polyps and their microenvironment, serving as a framework for precision surveillance and prevention of CRC.

Deciphering genomic alterations in colorectal cancer through transcriptional subtype-based network analysis

Both transcriptional subtype and signaling network analyses have proved useful in cancer genomics research. However, these two approaches are usually applied in isolation in existing studies. We reason that deciphering genomic alterations based on cancer transcriptional subtypes may help reveal subtype-specific driver networks and provide insights for the development of personalized therapeutic strategies. In this study, we defined transcriptional subtypes for colorectal cancer (CRC) and identified driver networks/pathways for each subtype. Applying consensus clustering to a patient cohort with 1173 samples identified three transcriptional subtypes, which were validated in an independent cohort with 485 samples. The three subtypes were characterized by different transcriptional programs related to normal adult colon, early colon embryonic development, and epithelial mesenchymal transition, respectively. They also showed statistically different clinical outcomes. For each subtype, we mapped somatic mutation and copy number variation data onto an integrated signaling network and identified subtype-specific driver networks using a random walk-based strategy. We found that genomic alterations in the Wnt signaling pathway were common among all three subtypes; however, unique combinations of pathway alterations including Wnt, VEGF and Notch drove distinct molecular and clinical phenotypes in different CRC subtypes. Our results provide a coherent and integrated picture of human CRC that links genomic alterations to molecular and clinical consequences, and which provides insights for the development of personalized therapeutic strategies for different CRC subtypes.

Developmentally programmed 3' CpG island methylation confers tissue- and cell-type-specific transcriptional activation

During development, a small but significant number of CpG islands (CGIs) become methylated. The timing of developmentally programmed CGI methylation and associated mechanisms of transcriptional regulation during cellular differentiation, however, remain poorly characterized. Here, we used genome-wide DNA methylation microarrays to identify epigenetic changes during human embryonic stem cell (hESC) differentiation. We discovered a group of CGIs associated with developmental genes that gain methylation after hESCs differentiate. Conversely, erasure of methylation was observed at the identified CGIs during subsequent reprogramming to induced pluripotent stem cells (iPSCs), further supporting a functional role for the CGI methylation. Both global gene expression profiling and quantitative reverse transcription-PCR (RT-PCR) validation indicated opposing effects of CGI methylation in transcriptional regulation during differentiation, with promoter CGI methylation repressing and 3' CGI methylation activating transcription. By studying diverse human tissues and mouse models, we further confirmed that developmentally programmed 3' CGI methylation confers tissue- and cell-type-specific gene activation in vivo. Importantly, luciferase reporter assays provided evidence that 3' CGI methylation regulates transcriptional activation via a CTCF-dependent enhancer-blocking mechanism. These findings expand the classic view of mammalian CGI methylation as a mechanism for transcriptional silencing and indicate a functional role for 3' CGI methylation in developmental gene regulation.

Abnormal Wnt and PI3Kinase signaling in the malformed intestine of lama5 deficient mice

Laminins are major constituents of basement membranes and are essential for tissue homeostasis. Laminin-511 is highly expressed in the intestine and its absence causes severe malformation of the intestine and embryonic lethality. To understand the mechanistic role of laminin-511 in tissue homeostasis, we used RNA profiling of embryonic intestinal tissue of lama5 knockout mice and identified a lama5 specific gene expression signature. By combining cell culture experiments with mediated knockdown approaches, we provide a mechanistic link between laminin α5 gene deficiency and the physiological phenotype. We show that laminin α5 plays a crucial role in both epithelial and mesenchymal cell behavior by inhibiting Wnt and activating PI3K signaling. We conclude that conflicting signals are elicited in the absence of lama5, which alter cell adhesion, migration as well as epithelial and muscle differentiation. Conversely, adhesion to laminin-511 may serve as a potent regulator of known interconnected PI3K/Akt and Wnt signaling pathways. Thus deregulated adhesion to laminin-511 may be instrumental in diseases such as human pathologies of the gut where laminin-511 is abnormally expressed as it is shown here.

Mouse conjunctival forniceal gene expression during postnatal development and its regulation by Kruppel-like factor 4

PURPOSE

To identify the changes in postnatal mouse conjunctival forniceal gene expression and their regulation by Klf4 during the eye-opening stage when the goblet cells first appear.

METHODS

Laser microdissection (LMD) was used to collect conjunctival forniceal cells from postnatal (PN) day 9, PN14 and PN20 wild-type (WT), and PN14 Klf4-conditional null (Klf4CN) mice, in which goblet cells are absent, developing, present, and missing, respectively. Microarrays were used to compare gene expression among these groups. Expression of selected genes was validated by quantitative RT-PCR, and spatiotemporal expression was assessed by in situ hybridization.

RESULTS

This study identified 668, 251, 1160, and 139 transcripts that were increased and 492, 377, 1419, and 57 transcripts that were decreased between PN9 and PN14, PN14 and PN20, PN9 and PN20, and PN14 WT and Klf4CN conjunctiva, respectively. Transcripts encoding transcription factors Spdef, FoxA1, and FoxA3 that regulate goblet cell development in other mucosal epithelia, and epithelium-specific Ets (ESE) transcription factor family members were increased during conjunctival development. Components of pathways related to the mesenchymal-epithelial transition, glycoprotein biosynthesis, mucosal immunity, signaling, and endocytic and neural regulation were increased during conjunctival development. Conjunctival Klf4 target genes differed significantly from the previously identified corneal Klf4 target genes, implying tissue-dependent regulatory targets for Klf4.

CONCLUSIONS

The changes in gene expression accompanying mouse conjunctival development were identified, and the role of Klf4 in this process was determined. This study provides new probes for examining conjunctival development and function and reveals that the gene regulatory network necessary for goblet cell development is conserved across different mucosal epithelia.

Colitis locus on chromosome 2 impacting the severity of early-onset disease in mice deficient in GPX1 and GPX2

BACKGROUND

Genetic background has a profound effect on inflammatory bowel disease. The Gpx1 and Gpx2 double knockout (GPX1/2-DKO) mice on a mixed C57BL/6 (B6) and 129S1/SvimJ (129) background exhibit spontaneous ileocolitis. The DKO mice on a B6 background have mild ileocolitis. We characterized the 129 DKO mice to identify a genetic locus affecting disease severity.

METHODS

We backcrossed B6;129 DKO mice to 129 and analyzed for ileocolitis penetrance and severity at N5, N7, and N10. By correlating disease severity with single-nucleotide polymorphism (SNP) markers, we identified a colitis locus.

RESULTS

As early as 9 days of age, 129 DKO N5 and N10 mice showed disease signs and morbidity. The N10 DKO mice had the severest colitis with nearly complete penetrance and high morbidity compared with other generations or backgrounds. 129 DKO mice had elevated colonic KC and SAA3 expression, shorter colon length, and cecal E. coli overgrowth compared to B6 DKO mice. Analysis of the B6 loci in 129 N5, N7, and N10 cohorts pointed to a region of chromosome 2: 119 Mbp contributing to mild symptoms.

CONCLUSIONS

GPX1/2-DKO mice on 129 genetic background have the most aggressive colitis compared to B6;129 and B6 colonies. A B6 locus significantly contributing the resistance resides on chromosome 2: 119 Mbp. This region coincides with cytokine-deficiency-induced colitis susceptibility, Cdcs3, identified in the resistant B6 and sensitive C3H/HeJBir (C3Bir) with IL-10 deficiency. A three-way SNP analysis between 129, B6, and C3Bir locus points the major candidate genes to B2m, Dnajc17, Duox2, Pla2g4b, Pla2g4e, Pla2g4f and Slc30a4.

TMPRSS4 induces invasion and epithelial-mesenchymal transition through upregulation of integrin alpha5 and its signaling pathways

TMPRSS4 is a novel type II transmembrane serine protease that is highly expressed on the cell surface in pancreatic, thyroid and other cancer tissues, although its oncogenic significance and molecular mechanisms are unknown. Previously, we have shown that TMPRSS4 promotes invasion, migration and metastasis of human tumor cells by facilitating an epithelial-mesenchymal transition (EMT). In this study, we explored the molecular basis underlying TMPRSS4-mediated effects. We show that multiple downstream signaling pathways, including focal adhesion kinase (FAK), extracellular signal-regulated kinase (ERK), Akt, Src and Rac1, are activated by TMPRSS4 expression and that FAK signaling and ERK activation are required for TMPRSS4-induced invasiveness and EMT, including cadherin switch. Inhibition of PI3K or Src reduced invasiveness and actin rearrangement mediated by TMPRSS4 without restoring E-cadherin expression. Downregulation of E-cadherin was required for TMPRSS4-mediated effects but was not sufficient to induce EMT and invasion. TMPRSS4 induced integrin alpha5 expression and its signal transduction, leading to invasiveness and EMT accompanied by downregulation of E-cadherin. Functional blocking confirmed that integrin alpha5beta1 is a critical signaling molecule that is sufficient to induce TMPRSS4-mediated effects. Immunohistochemical analysis showed that TMPRSS4 expression was significantly higher in human colorectal cancer tissues from advanced stages than in that of early stage. Furthermore, upregulation of TMPRSS4 was correlated with enhanced integrin alpha5 expression. These observations implicate integrin alpha5 upregulation as a molecular mechanism by which TMPRSS4 induces invasion and contributes to cancer progression.

A unified mixed effects model for gene set analysis of time course microarray experiments

Methods for gene set analysis test for coordinated changes of a group of genes involved in the same biological process or molecular pathway. Higher statistical power is gained for gene set analysis by combining weak signals from a number of individual genes in each group. Although many gene set analysis methods have been proposed for microarray experiments with two groups, few can be applied to time course experiments. We propose a unified statistical model for analyzing time course experiments at the gene set level using random coefficient models, which fall into the more general class of mixed effects models. These models include a systematic component that models the mean trajectory for the group of genes, and a random component (the random coefficients) that models how each gene's trajectory varies about the mean trajectory. We show that the proposed model (1) outperforms currently available methods at discriminating gene sets differentially changed over time from null gene sets; (2) provides more stable results that are less affected by sampling variations; (3) models dependency among genes adequately and preserves type I error rate; and (4) allows for gene ranking based on predicted values of the random effects. We describe simulation studies using gene expression data with "real life" correlations and we demonstrate the proposed random coefficient model using a mouse colon development time course dataset. The agreement between results of the proposed random coefficient model and the previous reports for this proof-of-concept trial further validates this methodology, which provides a unified statistical model for systems analysis of microarray experiments with complex experimental designs when re-sampling based methods are difficult to apply.

Expression patterns of the aquaporin gene family during renal development: influence of genetic variability

High-throughput analyses have shown that aquaporins (AQPs) belong to a cluster of genes that are differentially expressed during kidney organogenesis. However, the spatiotemporal expression patterns of the AQP gene family during tubular maturation and the potential influence of genetic variation on these patterns and on water handling remain unknown. We investigated the expression patterns of all AQP isoforms in fetal (E13.5 to E18.5), postnatal (P1 to P28), and adult (9 weeks) kidneys of inbred (C57BL/6J) and outbred (CD-1) mice. Using quantitative polymerase chain reaction (PCR), we evidenced two mRNA patterns during tubular maturation in C57 mice. The AQPs 1-7-11 showed an early (from E14.5) and progressive increase to adult levels, similar to the mRNA pattern observed for proximal tubule markers (Megalin, NaPi-IIa, OAT1) and reflecting the continuous increase in renal cortical structures during development. By contrast, AQPs 2-3-4 showed a later (E15.5) and more abrupt increase, with transient postnatal overexpression. Most AQP genes were expressed earlier and/or stronger in maturing CD-1 kidneys. Furthermore, adult CD-1 kidneys expressed more AQP2 in the collecting ducts, which was reflected by a significant delay in excreting a water load. The expression patterns of proximal vs. distal AQPs and the earlier expression in the CD-1 strain were confirmed by immunoblotting and immunostaining. These data (1) substantiate the clustering of important genes during tubular maturation and (2) demonstrate that genetic variability influences the regulation of the AQP gene family during tubular maturation and water handling by the mature kidney.

Transcriptional recapitulation and subversion of embryonic colon development by mouse colon tumor models and human colon cancer

Colon tumors from four independent mouse models and 100 human colorectal cancers all exhibited striking recapitulation of embryonic colon gene expression from embryonic days 13.5-18.5.

Background

The expression of carcino-embryonic antigen by colorectal cancer is an example of oncogenic activation of embryonic gene expression. Hypothesizing that oncogenesis-recapitulating-ontogenesis may represent a broad programmatic commitment, we compared gene expression patterns of human colorectal cancers (CRCs) and mouse colon tumor models to those of mouse colon development embryonic days 13.5-18.5.

Results

We report here that 39 colon tumors from four independent mouse models and 100 human CRCs encompassing all clinical stages shared a striking recapitulation of embryonic colon gene expression. Compared to normal adult colon, all mouse and human tumors over-expressed a large cluster of genes highly enriched for functional association to the control of cell cycle progression, proliferation, and migration, including those encoding MYC, AKT2, PLK1 and SPARC. Mouse tumors positive for nuclear β-catenin shifted the shared embryonic pattern to that of early development. Human and mouse tumors differed from normal embryonic colon by their loss of expression modules enriched for tumor suppressors (EDNRB, HSPE, KIT and LSP1). Human CRC adenocarcinomas lost an additional suppressor module (IGFBP4, MAP4K1, PDGFRA, STAB1 and WNT4). Many human tumor samples also gained expression of a coordinately regulated module associated with advanced malignancy (ABCC1, FOXO3A, LIF, PIK3R1, PRNP, TNC, TIMP3 and VEGF).

Conclusion

Cross-species, developmental, and multi-model gene expression patterning comparisons provide an integrated and versatile framework for definition of transcriptional programs associated with oncogenesis. This approach also provides a general method for identifying pattern-specific biomarkers and therapeutic targets. This delineation and categorization of developmental and non-developmental activator and suppressor gene modules can thus facilitate the formulation of sophisticated hypotheses to evaluate potential synergistic effects of targeting within- and between-modules for next-generation combinatorial therapeutics and improved mouse models.

Small proline-rich proteins (SPRR) function as SH3 domain ligands, increase resistance to injury and are associated with epithelial-mesenchymal transition (EMT) in cholangiocytes

BACKGROUND/AIMS

Deficient biliary epithelial cell (BEC) expression of small proline-rich protein (SPRR) 2A in IL-6(-/-) mice is associated with defective biliary barrier function after bile duct ligation. And numerous gene array expression studies show SPRR2A to commonly be among the most highly up-regulated genes in many non-squamous, stressed and remodeling barrier epithelia. Since the function of SPRR in these circumstances is unknown, we tested the exploratory hypothesis that BEC SPRR2A expression contributes to BEC barrier function and wound repair.

METHODS

The effect of SPRR2A expression was studied in primary mouse BEC cultures; in a BEC cell line after forced overexpression of SPRR2A; and in human livers removed at the time of liver transplantation.

RESULTS

Forced SPRR2A overexpression showed that it functions as a SH3 domain ligand that increases resistance to oxidative injury and promotes wound restitution by enhancing migration and acquisition of mesenchymal characteristics. Low confluency non-neoplastic mouse BEC cultures show a phenotype similar to the stable transfectants, as did spindle-shaped BEC participating in atypical ductular reactions in primary biliary cirrhosis.

CONCLUSIONS

These observations suggest that SPRR2A-related BEC barrier modifications represent a novel, but widely utilized and evolutionarily conserved, response to stress that is worthy of further study.

Map of differential transcript expression in the normal human large intestine

While there is considerable research related to using differential gene expression to predict disease phenotype classification, e.g., neoplastic tissue from nonneoplastic controls, there is little understanding of the range of expression in normal tissues. Understanding patterns of gene expression in nonneoplastic tissue, including regional anatomic expression changes within an organ, is vital to understanding gene expression changes in diseased tissue. To explore the gene expression change along the proximal-distal axis of the large intestine, we analyzed microarray data in 184 normal human specimens using univariate and multivariate techniques. We found 219 probe sets that were differentially expressed between the proximal and distal colorectal regions and 115 probe sets that were differentially expressed between the terminal segments, i.e., the cecum and rectum. We did not observe any probe sets that were statistically different between any two contiguous colorectal segments. The dominant expression pattern (65 probe sets) follows a dichotomous expression pattern consistent with the midgut-hindgut embryonic origins of the gut while a second pattern (50 probe sets) depicts a gradual change in transcript levels from the cecum to the rectum. While the dichotomous pattern includes roughly equal numbers of probe sets that are elevated proximally and distally, nearly all probe sets that show a gradual change demonstrate increasing expression levels moving from proximal to distal segments. These patterns describe an expression map of individual transcript variation as well as multigene expression patterns along the large intestine. This is the first gene expression map of an entire human organ.

Hedgehog signaling in development and homeostasis of the gastrointestinal tract

The Hedgehog family of secreted morphogenetic proteins acts through a complex evolutionary conserved signaling pathway to regulate patterning events during development and in the adult organism. In this review I discuss the role of Hedgehog signaling in the development, postnatal maintenance, and carcinogenesis of the gastrointestinal tract. Three mammalian hedgehog genes, sonic hedgehog (Shh), indian hedgehog (Ihh), and desert hedgehog (Dhh) have been identified. Shh and Ihh are important endodermal signals in the endodermal-mesodermal cross-talk that patterns the developing gut tube along different axes. Mutations in Shh, Ihh, and downstream signaling molecules lead to a variety of gross malformations of the murine gastrointestinal tract including esophageal atresia, tracheoesophageal fistula, annular pancreas, midgut malrotation, and duodenal and anal atresia. These congenital malformations are also found in varying constellations in humans, suggesting a possible role for defective Hedgehog signaling in these patients. In the adult, Hedgehog signaling regulates homeostasis in several endoderm-derived epithelia, for example, the stomach, intestine, and pancreas. Finally, growth of carcinomas of the proximal gastrointestinal tract such as esophageal, gastric, biliary duct, and pancreatic cancers may depend on Hedgehog signaling offering a potential avenue for novel therapy for these aggressive cancers.

Using genomics to understand intestinal biology

The use of microarrays to evaluate the transcriptome has transformed our view of biology. In addition to the focused, hypothesis-testing studies that we have traditionally conducted in cell biology, we are now able to see global changes within the entire system of the cell in response to a treatment. By examining a biological question under multiple complementary perturbations model systems (e.g. yeast, C. Elegans) have revealed new complexity that would have been impossible to see on a gene-by-gene approach. Unfortunately, beyond the use of transcript profiles to define the molecular signature of diseases (e.g. cancer), transcriptomics has not been extensively used to study intestinal biology. This review will provide a roadmap for effective use of gene expression profiling for biological research and will review some of the microarray work that has been done to better understand the nature of intestinal development and enterocyte differentiation.

Targeted imaging of colonic tumors in smad3-/- mice discriminates cancer and inflammation

The peripheral benzodiazepine receptor (PBR) is a trans-mitochondrial membrane protein that modulates steroid biosynthesis. Recently, up-regulation and nuclear localization of PBR has been shown to be associated with colon, prostate, and breast cancer. PBR has been targeted by the exogenous synthetic ligand, PK11195, for various purposes including imaging. To capitalize on these observations, we developed a high-throughput, noninvasive, in vivo imaging approach to detect spontaneously arising colonic tumors in mice using a novel PBR-targeted molecular imaging agent (NIR-conPK11195). NIR-conPK11195 localized and was retained in colonic adenomas and carcinomas in Smad3(-/-) mice but not in non-neoplastic hamartomas or chronically inflamed colonic tissue. Using a fluorescence signal-to-noise ratio of > or =4-fold 13 h after injection of the agent, we detected colonic tumors with a sensitivity of 67% and a specificity of 86% in a cohort of 37 Smad3(-/-) mice and control littermates. Furthermore, using oral administration of dextran sulfate to induce colonic inflammation, we showed that the clearance profile of NIR-conPK11195 distinguished transient uptake in inflammatory tissue from longer term retention in tumors. Taken together, these results indicate that NIR-conPK11195 is a promising optical molecular imaging tool to rapidly screen for colonic tumors in mice and to discriminate inflammation from cancer.

Differences in the placental membrane cytokine response: a possible explanation for the racial disparity in preterm birth

PROBLEM

The prematurity rate is higher in African-Americans (AA) compared with Caucasians (C). As spontaneous preterm labor has been hypothesized to be a host inflammatory response disease racial differences in human placental membrane inflammatory cytokine and prostaglandin pathway gene expression patterns between AA and C were examined in this report.

METHOD OF STUDY

Placental membranes (amniochorion) collected from AA and C women from cesareans at term were maintained in an organ explant system and stimulated with endotoxin (lipopolysaccharide, LPS). Microarray analysis and enzyme-linked immunosorbent assay was performed on mRNAs and culture media from AA- and C-derived membranes to document any differences in mRNA expression and protein production of IL-1, IL-6, IL-8, IL-10 and expression of cyclooxygenase 1 (COX-1), COX-2 and 15-hydroxyprostaglandin dehydrogenase (PGDH).

RESULTS

Increased mRNA expression of IL-1, IL-8 and COX-2 in AA and IL-6, IL-10, COX-1 and PGDH in C were documented after LPS stimulation. Concentration of IL-1 was significantly higher in media derived from AA whereas IL-6 and IL-10 concentrations were higher in C with no differences observed in IL-8 after LPS stimulation compared with respective unstimulated controls.

CONCLUSION

These data document ethnic diversity in placental membrane immune response.