Systematic identification of novel cancer genes through analysis of deep shRNA perturbation screens

Systematic perturbation screens provide comprehensive resources for the elucidation of cancer driver genes. The perturbation of many genes in relatively few cell lines in such functional screens necessitates the development of specialized computational tools with sufficient statistical power. Here we developed APSiC (Analysis of Perturbation Screens for identifying novel Cancer genes) to identify genetic drivers and effectors in perturbation screens even with few samples. Applying APSiC to the shRNA screen Project DRIVE, APSiC identified well-known and novel putative mutational and amplified cancer genes across all cancer types and in specific cancer types. Additionally, APSiC discovered tumor-promoting and tumor-suppressive effectors, respectively, for individual cancer types, including genes involved in cell cycle control, Wnt/β-catenin and hippo signalling pathways. We functionally demonstrated that LRRC4B, a putative novel tumor-suppressive effector, suppresses proliferation by delaying cell cycle and modulates apoptosis in breast cancer. We demonstrate APSiC is a robust statistical framework for discovery of novel cancer genes through analysis of large-scale perturbation screens. The analysis of DRIVE using APSiC is provided as a web portal and represents a valuable resource for the discovery of novel cancer genes.     

Preferential Utilization of d-Lactate by Shewanella oneidensis

Shewanella oneidensis couples oxidation of lactate to respiration of many substrates. Here we report that llpR (l-lactate-positive regulator, SO_3460) encodes a positive regulator of l-lactate utilization distinct from previously studied regulators. We also demonstrate d-lactate inhibition of l-lactate utilization in S. oneidensis, resulting in preferential utilization of the d isomer.     

RPM-1 Uses Both Ubiquitin Ligase and Phosphatase-Based Mechanisms to Regulate DLK-1 during Neuronal Development

The Pam/Highwire/RPM-1 (PHR) proteins are key regulators of neuronal development that function in axon extension and guidance, termination of axon outgrowth, and synapse formation. Outside of development, the PHR proteins also regulate axon regeneration and Wallerian degeneration. The PHR proteins function in part by acting as ubiquitin ligases that degrade the Dual Leucine zipper-bearing Kinase (DLK). Here, we show that the Caenorhabditis elegans PHR protein, Regulator of Presynaptic Morphology 1 (RPM-1), also utilizes a phosphatase-based mechanism to regulate DLK-1. Using mass spectrometry, we identified Protein Phosphatase Magnesium/Manganese dependent 2 (PPM-2) as a novel RPM-1 binding protein. Genetic, transgenic, and biochemical studies indicated that PPM-2 functions coordinately with the ubiquitin ligase activity of RPM-1 and the F-box protein FSN-1 to negatively regulate DLK-1. PPM-2 acts on S874 of DLK-1, a residue implicated in regulation of DLK-1 binding to a short, inhibitory isoform of DLK-1 (DLK-1S). Our study demonstrates that PHR proteins function through both phosphatase and ubiquitin ligase mechanisms to inhibit DLK. Thus, PHR proteins are potentially more accurate and sensitive regulators of DLK than originally thought. Our results also highlight an important and expanding role for the PP2C phosphatase family in neuronal development.     

Uncovering Genomic Features and Maternal Origin of Korean Native Chicken by Whole Genome Sequencing

The Korean Native Chicken (KNC) is an important endemic biological resource in Korea. While numerous studies have been conducted exploring this breed, none have used next-generation sequencing to identify its specific genomic features. We sequenced five strains of KNC and identified 10.9 million SNVs and 1.3 million InDels. Through the analysis, we found that the highly variable region common to all 5 strains had genes like PCHD15, CISD1, PIK3C2A, and NUCB2 that might be related to the phenotypic traits of the chicken such as auditory sense, growth rate and egg traits. In addition, we assembled unaligned reads that could not be mapped to the reference genome. By assembling the unaligned reads, we were able to present genomic sequences characteristic to the KNC. Based on this, we also identified genes related to the olfactory receptors and antigen that are common to all 5 strains. Finally, through the reconstructed mitochondrial genome sequences, we performed phylogenomic analysis and elucidated the maternal origin of the artificially restored KNC. Our results revealed that the KNC has multiple maternal origins which are in agreement with Korea''s history of chicken breed imports. The results presented here provide a valuable basis for future research on genomic features of KNC and further understanding of KNC''s origin.     

Evolution of the B3 DNA Binding Superfamily: New Insights into REM Family Gene Diversification

Background

The B3 DNA binding domain includes five families: auxin response factor (ARF), abscisic acid-insensitive3 (ABI3), high level expression of sugar inducible (HSI), related to ABI3/VP1 (RAV) and reproductive meristem (REM). The release of the complete genomes of the angiosperm eudicots Arabidopsis thaliana and Populus trichocarpa, the monocot Orysa sativa, the bryophyte Physcomitrella patens,the green algae Chlamydomonas reinhardtii and Volvox carteri and the red algae Cyanidioschyzon melorae provided an exceptional opportunity to study the evolution of this superfamily.

Methodology

In order to better understand the origin and the diversification of B3 domains in plants, we combined comparative phylogenetic analysis with exon/intron structure and duplication events. In addition, we investigated the conservation and divergence of the B3 domain during the origin and evolution of each family.

Conclusions

Our data indicate that showed that the B3 containing genes have undergone extensive duplication events, and that the REM family B3 domain has a highly diverged DNA binding. Our results also indicate that the founding member of the B3 gene family is likely to be similar to the ABI3/HSI genes found in C. reinhardtii and V. carteri. Among the B3 families, ABI3, HSI, RAV and ARF are most structurally conserved, whereas the REM family has experienced a rapid divergence. These results are discussed in light of their functional and evolutionary roles in plant development.     

The Linker for Activation of B Cells (LAB)/Non-T Cell Activation Linker (NTAL) Regulates Triggering Receptor Expressed on Myeloid Cells (TREM)-2 Signaling and Macrophage Inflammatory Responses Independently of the Linker for Activation of T Cells

Triggering receptor expressed on myeloid cells-2 (TREM-2) is rapidly emerging as a key regulator of the innate immune response via its regulation of macrophage inflammatory responses. Here we demonstrate that proximal TREM-2 signaling parallels other DAP12-based receptor systems in its use of Syk and Src-family kinases. However, we find that the linker for activation of T cells (LAT) is severely reduced as monocytes differentiate into macrophages and that TREM-2 exclusively uses the linker for activation of B cells (LAB encoded by the gene Lat2^−/−) to mediate downstream signaling. LAB is required for TREM-2-mediated activation of Erk1/2 and dampens proximal TREM-2 signals through a novel LAT-independent mechanism resulting in macrophages with proinflammatory properties. Thus, Lat2^−/− macrophages have increased TREM-2-induced proximal phosphorylation, and lipopolysaccharide stimulation of these cells leads to increased interleukin-10 (IL-10) and decreased IL-12p40 production relative to wild type cells. Together these data identify LAB as a critical, LAT-independent regulator of TREM-2 signaling and macrophage development capable of controlling subsequent inflammatory responses.     

Functional Screening of Metagenome and Genome Libraries for Detection of Novel Flavonoid-Modifying Enzymes

The functional detection of novel enzymes other than hydrolases from metagenomes is limited since only a very few reliable screening procedures are available that allow the rapid screening of large clone libraries. For the discovery of flavonoid-modifying enzymes in genome and metagenome clone libraries, we have developed a new screening system based on high-performance thin-layer chromatography (HPTLC). This metagenome extract thin-layer chromatography analysis (META) allows the rapid detection of glycosyltransferase (GT) and also other flavonoid-modifying activities. The developed screening method is highly sensitive, and an amount of 4 ng of modified flavonoid molecules can be detected. This novel technology was validated against a control library of 1,920 fosmid clones generated from a single Bacillus cereus isolate and then used to analyze more than 38,000 clones derived from two different metagenomic preparations. Thereby we identified two novel UDP glycosyltransferase (UGT) genes. The metagenome-derived gtfC gene encoded a 52-kDa protein, and the deduced amino acid sequence was weakly similar to sequences of putative UGTs from Fibrisoma and Dyadobacter. GtfC mediated the transfer of different hexose moieties and exhibited high activities on flavones, flavonols, flavanones, and stilbenes and also accepted isoflavones and chalcones. From the control library we identified a novel macroside glycosyltransferase (MGT) with a calculated molecular mass of 46 kDa. The deduced amino acid sequence was highly similar to sequences of MGTs from Bacillus thuringiensis. Recombinant MgtB transferred the sugar residue from UDP-glucose effectively to flavones, flavonols, isoflavones, and flavanones. Moreover, MgtB exhibited high activity on larger flavonoid molecules such as tiliroside.     

Ankyrin Repeat Domain Protein 2 and Inhibitor of DNA Binding 3 Cooperatively Inhibit Myoblast Differentiation by Physical Interaction

Ankyrin repeat domain protein 2 (ANKRD2) translocates from the nucleus to the cytoplasm upon myogenic induction. Overexpression of ANKRD2 inhibits C2C12 myoblast differentiation. However, the mechanism by which ANKRD2 inhibits myoblast differentiation is unknown. We demonstrate that the primary myoblasts of mdm (muscular dystrophy with myositis) mice (pMB^mdm) overexpress ANKRD2 and ID3 (inhibitor of DNA binding 3) proteins and are unable to differentiate into myotubes upon myogenic induction. Although suppression of either ANKRD2 or ID3 induces myoblast differentiation in mdm mice, overexpression of ANKRD2 and inhibition of ID3 or vice versa is insufficient to inhibit myoblast differentiation in WT mice. We identified that ANKRD2 and ID3 cooperatively inhibit myoblast differentiation by physical interaction. Interestingly, although MyoD activates the Ankrd2 promoter in the skeletal muscles of wild-type mice, SREBP-1 (sterol regulatory element binding protein-1) activates the same promoter in the skeletal muscles of mdm mice, suggesting the differential regulation of Ankrd2. Overall, we uncovered a novel pathway in which SREBP-1/ANKRD2/ID3 activation inhibits myoblast differentiation, and we propose that this pathway acts as a critical determinant of the skeletal muscle developmental program.     

Identification of Interphase Functions for the NIMA Kinase Involving Microtubules and the ESCRT Pathway

The Never in Mitosis A (NIMA) kinase (the founding member of the Nek family of kinases) has been considered a mitotic specific kinase with nuclear restricted roles in the model fungus Aspergillus nidulans. By extending to A. nidulans the results of a synthetic lethal screen performed in Saccharomyces cerevisiae using the NIMA ortholog KIN3, we identified a conserved genetic interaction between nimA and genes encoding proteins of the Endosomal Sorting Complex Required for Transport (ESCRT) pathway. Absence of ESCRT pathway functions in combination with partial NIMA function causes enhanced cell growth defects, including an inability to maintain a single polarized dominant cell tip. These genetic insights suggest NIMA potentially has interphase functions in addition to its established mitotic functions at nuclei. We therefore generated endogenously GFP-tagged NIMA (NIMA-GFP) which was fully functional to follow its interphase locations using live cell spinning disc 4D confocal microscopy. During interphase some NIMA-GFP locates to the tips of rapidly growing cells and, when expressed ectopically, also locates to the tips of cytoplasmic microtubules, suggestive of non-nuclear interphase functions. In support of this, perturbation of NIMA function either by ectopic overexpression or through partial inactivation results in marked cell tip growth defects with excess NIMA-GFP promoting multiple growing cell tips. Ectopic NIMA-GFP was found to locate to the plus ends of microtubules in an EB1 dependent manner, while impairing NIMA function altered the dynamic localization of EB1 and the cytoplasmic microtubule network. Together, our genetic and cell biological analyses reveal novel non-nuclear interphase functions for NIMA involving microtubules and the ESCRT pathway for normal polarized fungal cell tip growth. These insights extend the roles of NIMA both spatially and temporally and indicate that this conserved protein kinase could help integrate cell cycle progression with polarized cell growth.     

Fast and Rigorous Computation of Gene and Pathway Scores from SNP-Based Summary Statistics

Integrating single nucleotide polymorphism (SNP) p-values from genome-wide association studies (GWAS) across genes and pathways is a strategy to improve statistical power and gain biological insight. Here, we present Pascal (Pathway scoring algorithm), a powerful tool for computing gene and pathway scores from SNP-phenotype association summary statistics. For gene score computation, we implemented analytic and efficient numerical solutions to calculate test statistics. We examined in particular the sum and the maximum of chi-squared statistics, which measure the strongest and the average association signals per gene, respectively. For pathway scoring, we use a modified Fisher method, which offers not only significant power improvement over more traditional enrichment strategies, but also eliminates the problem of arbitrary threshold selection inherent in any binary membership based pathway enrichment approach. We demonstrate the marked increase in power by analyzing summary statistics from dozens of large meta-studies for various traits. Our extensive testing indicates that our method not only excels in rigorous type I error control, but also results in more biologically meaningful discoveries.