Silencing LRH-1 in colon cancer cell lines impairs proliferation and alters gene expression programs

Colorectal cancers (CRCs) account for nearly 10% of all cancer deaths in industrialized countries. Recent evidence points to a central role for the nuclear receptor liver receptor homolog-1 (LRH-1) in intestinal tumorigenesis. Interaction of LRH-1 with the Wnt/β-catenin pathway, highly active in a critical subpopulation of CRC cells, underscores the importance of elucidating LRH-1's role in this disease. Reduction of LRH-1 diminishes tumor burden in murine models of CRC; however, it is not known whether LRH-1 is required for tumorigenesis, for proliferation, or for both. In this work, we address this question through shRNA-mediated silencing of LRH-1 in established CRC cell lines. LRH-1 mRNA knockdown results in significantly impaired proliferation in a cell line highly expressing the receptor and more modest impairment in a cell line with moderate LRH-1 expression. Cell-cycle analysis shows prolongation of G0/G1 with LRH-1 silencing, consistent with LRH-1 cell-cycle influences in other tissues. Cluster analysis of microarray gene expression demonstrates significant genome wide alterations with major effects in cell-cycle regulation, signal transduction, bile acid and cholesterol metabolism, and control of apoptosis. This study demonstrates a critical proproliferative role for LRH-1 in established colon cancer cell lines. LRH-1 exerts its effects via multiple signaling networks. Our results suggest that selected CRC patients could benefit from LRH-1 inhibitors

Silencing LRH-1 in colon cancer cell lines impairs proliferation and alters gene expression programs

Colorectal cancers (CRCs) account for nearly 10% of all cancer deaths in industrialized countries. Recent evidence points to a central role for the nuclear receptor liver receptor homolog-1 (LRH-1) in intestinal tumorigenesis. Interaction of LRH-1 with the Wnt/β-catenin pathway, highly active in a critical subpopulation of CRC cells, underscores the importance of elucidating LRH-1’s role in this disease. Reduction of LRH-1 diminishes tumor burden in murine models of CRC; however, it is not known whether LRH-1 is required for tumorigenesis, for proliferation, or for both. In this work, we address this question through shRNA-mediated silencing of LRH-1 in established CRC cell lines. LRH-1 mRNA knockdown results in significantly impaired proliferation in a cell line highly expressing the receptor and more modest impairment in a cell line with moderate LRH-1 expression. Cell-cycle analysis shows prolongation of G0/G1 with LRH-1 silencing, consistent with LRH-1 cell-cycle influences in other tissues. Cluster analysis of microarray gene expression demonstrates significant genome wide alterations with major effects in cell-cycle regulation, signal transduction, bile acid and cholesterol metabolism, and control of apoptosis. This study demonstrates a critical proproliferative role for LRH-1 in established colon cancer cell lines. LRH-1 exerts its effects via multiple signaling networks. Our results suggest that selected CRC patients could benefit from LRH-1 inhibitors

The <em>C. elegans</em> Rab Family: Identification, Classification and Toolkit Construction

<div><p>Rab monomeric GTPases regulate specific aspects of vesicle transport in eukaryotes including coat recruitment, uncoating, fission, motility, target selection and fusion. Moreover, individual Rab proteins function at specific sites within the cell, for example the ER, golgi and early endosome. Importantly, the localization and function of individual Rab subfamily members are often conserved underscoring the significant contributions that model organisms such as <em>Caenorhabditis elegans</em> can make towards a better understanding of human disease caused by Rab and vesicle trafficking malfunction. With this in mind, a bioinformatics approach was first taken to identify and classify the complete <em>C. elegans</em> Rab family placing individual Rabs into specific subfamilies based on molecular phylogenetics. For genes that were difficult to classify by sequence similarity alone, we did a comparative analysis of intron position among specific subfamilies from yeast to humans. This two-pronged approach allowed the classification of 30 out of 31 <em>C. elegans</em> Rab proteins identified here including <em>Rab31/Rab50</em>, a likely member of the last eukaryotic common ancestor (LECA). Second, a molecular toolset was created to facilitate research on biological processes that involve Rab proteins. Specifically, we used Gateway-compatible <em>C. elegans</em> ORFeome clones as starting material to create 44 full-length, sequence-verified, dominant-negative (DN) and constitutive active (CA) <em>rab</em> open reading frames (ORFs). Development of this toolset provided independent research projects for students enrolled in a research-based molecular techniques course at California State University, East Bay (CSUEB).</p> </div

Bootstrap scores for specific terminal clades from 13 additional phylogenetic reconstructions.

<p>Thirteen additional phylogenetic analyses were performed using a combination of statistical methods. Phylogenetic reconstruction methods include Neighbor Joining (NJ), maximum likelihood (ML) or minimum evolution (ME). Amino acid substitution methods include Poisson (Po), JTT, Dayhoff (D), Equal Input (EI) and WAG. Gap deletion treatments include partial (Par) or pairwise (Pai) and rates and patterns of evolution include gamma distributed (+G), invariant sites (+I) or uniform (all others). All phylogenetic reconstructions were performed in MEGA5. Specific orthologous and/or paralogous clades include both worm and human Rab proteins. New orthologous clusters include bootstrap scores supporting the Rab44/4R79.2, Rab45/C33D12.6 and Rab23/ZK669.5 pairs. The new paralogous cluster includes the bootstrap scores that support the Rab44,4R79.2,Rab45 and C33D12.6 terminal clade.</p

A chladogram of Rab family members from <i>C. elegans</i> and <i>H. sapien</i>s.

<p>The evolutionary history was inferred using the Neighbor-Joining phylogenetic reconstruction method. The tree is rooted with the natural outlying clade, Rab28. The optimal tree is shown with the percentage of replicate trees (>40) in which the associated genes cluster together in the bootstrap test (500 replicates) provided next to each branch. The tree is drawn to emphasize topology. The evolutionary distances were computed using the JTT amino acid substitution method and are in the units of the number of amino acid differences per site. Evolutionary analyses were conducted using MEGA5. Clades marked with red, orange or yellow circles indicate their degree of stability under a variety of phylogenetic reconstruction parameters (see text and methods for details). Red = 14/14, orange = 13/14 and yellow = 12/14 trees. Genes highlighted with black circles represent putative orphan <i>C. elegan</i> Rabs (lacking a human ortholog). For simplicity, closely related splice variants and well-supported human-specific clades were deleted (see methods for details).</p

Comparative analysis of intron position among diverse Rab subfamily members.

<p>A) Cladogram indicating evolutionary relationships of 18 species examined here <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0049387#pone.0049387-Dunn1" target="_blank">[128]</a>–<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0049387#pone.0049387-Csuros1" target="_blank">[130]</a>. Ecdy. = Ecdysozoa, Chromal. = Chromalveolata. For species abbreviations see Methods. B) An ML tree of Opisthokonts created from the MSA used to map intron positions. Bootstrap support (100 replicates) is indicated for each subfamily cluster. C) For each subfamily, the number of times a <u>S</u>ubfamily <u>S</u>pecific <u>C</u>onserved <u>I</u>ntron <u>P</u>osition (SSCIP) involving the indicated number of species was observed (gray bars), compared to what is expected by chance (black diamonds). The difference between observed and expected is statistically significant where indicated. *P(Monte Carlo) <0.05. ***P(Monte Carlo)≤0.00001. The Rab<i>31, 6, 5, 22, 34, 21</i> and <i>23</i> subfamilies include 17, 18, 17, 9, 9, 10, 14 and 12 species, respectively. D) and E) Heat map indicating number of introns within <i>Rab31</i> (D) or <i>Rab6</i> (E) that match SSCIPs from <i>Rab31, 5, 22, 21, 6, 34</i> and <i>23</i>. The circled number indicates the number of introns present in the MSA for each gene. % equals the percentage of introns that are shared with the true SSCIP. <i>C56E6.2</i> (D) and <i>Y71H2AM.12</i> (E) are highlighted red. Genbank Descriptions (if any) and RABDB! classifications are included. Classification abbreviations include: HypoRabX1 (H.RabX1), HypoRabX2 (H.RabX2), HypoRabX3 (H.RabX3) and MetazoaRabX3 (M.RabX3). F) A pairwise comparison of intron position conservation between specific genes (Rab31 at left, Rab6 at right) and their corresponding set of SSCIPs. Black diamonds plot the probability that a specific number of intron matches would be expected by chance for each set of conditions. Chart 1 plots a comparison of 5 introns with 7 SSCIPs (5×7). Chart 2∶4×7. Chart 3∶2×7. Chart 4∶4×6. Observed values for a subset of genes are indicated with P values estimated from the Monte Carlo simulation data (See text and methods). Species abbreviations are as in A. C) and F) 72 protosplice sites assumed.</p

Full-length ORF sequence of <i>tag-312</i> (Rab45) isolates reveals a different splice pattern than predicted.

<p>A) A nucleotide alignment of <i>tag-312</i> (FlOCS), genomic DNA and the predicted <i>tag-312</i> ORF at two regions where exon/intron splice junction differences were found. In the top alignment, FlOCS reveals a new intron, splitting predicted exon 7 into exons7a and 7b. In the bottom alignment, FlOCS reveals an alternate 5′ splice donor and 3′ splice acceptor for intron 8. Compare FlOCS-supported 5′ and 3′ splice sites boxed in bold to predicted 5′ and 3′ splices sites (boxed, not bold). B) Two multiple sequence alignments of Rab45 subfamily members spanning the two regions described in 5A above demonstrate the impact of FLOCS-supported gene structure differences. The alignment includes proteins from the nematodes, <i>Caenhorabditis elegans</i>, <i>tag-312</i>(FlOCS) and NP508523.1, <i>Caenhorabditis brenneri, Caenhorabditis remanei</i> and <i>Ascaris suum</i> in addition to <i>Xenopus laevis</i> (frog), <i>Homo sapiens</i> (human), <i>Monodelphis domestica</i> (opossum) and <i>Anolis carolinensis</i> (lizard). The intron that splits exon 7 into two creates a 15 amino acid deletion that is conserved among all species examined (top). The alternate intron 8 creates an Indel in a region of Rab45 that is conserved among nematodes only.</p

A list and description of Rab isolates created for the <i>C. elegans</i> ORFeome-based toolkit.

<p>WT, DN and CA clones included in the Rab Toolkit are given isolate names otherwise an explanation for its absence is provided. Subfamily classifications are based on Diekmann et al. 2011 <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0049387#pone.0049387-Diekmann1" target="_blank">[33]</a> and/or data presented here. The majority of <i>C. elegans rab</i> genes are predicted to have only one splice variant with the following exceptions: WormBase describes two splice variants for <i>rab-3</i> that code for proteins 233 and 219 amino acids (aa) in length. ORFeome project primers were designed to amplify the shorter isoform only. WormBase describes two splice variants for <i>4R79.2 (Rab44)</i> that code for proteins 311 aa and 395 aa in length. ORFeome project primers were designed to amplify the longer isoform only. Names listed under “other” are from WormBase or Pereira-Leal and Seabra (2001). Finally, while <i>rab-37</i> shows 100% identity with the Refseq protein NP_001041293 it contains an additional 5 amino acids at its N-terminus. See text for details.</p