The degree of postpartum metabolic challenge in dairy cows is associated with peripheral blood mononuclear cell transcriptome changes of the innate immune system

Dairy cows undergo a nutrient deficit immediately postpartum when lactational demands exceed nutrient intake. This occurs concurrently to an increased challenge due to bacterial and viral infections, yet ability for pathogen clearance is reduced despite a heightened and often host-damaging inflammatory response. We hypothesized that nutrient stress is associated with differences in the immune cell transcriptome. Our objective was therefore to investigate differentially expressed pathways (DEP) by RNA-seq in peripheral blood mononuclear cells harvested 3 weeks before and 1 week after calving from Holstein cows in low (L, n = 3) or high (H, n = 3) postpartum metabolic stress situations. Metabolic stress was defined by differences in circulating concentrations of glucose, fatty acids, and ketones postpartum. Cows in group H showed several upregulated DEP in relation to myeloid cell function and inflammatory response, as well as downregulation of the Th2 pathway. Principal components analysis showed that the transcriptome of group H postpartum samples was most different from all other samples. Differences in DE genes were noted even prepartum albeit fewer DE genes were identified and myeloid cell pathways in group H were generally downregulated at this time compared with group L. Samples within group L showed little difference between the two time points. We conclude that the metabolic phenotype of cows allowed us to identify differences in immune-regulatory pathways and that myeloid immune cells could play a dominant role in identifying these metabolically-associated differences that were demonstrated among a mixed mononuclear cell population.

Dynamic transcriptome profiles within spermatogonial and spermatocyte populations during postnatal testis maturation revealed by single-cell sequencing

Spermatogenesis is the process by which male gametes are formed from a self-renewing population of spermatogonial stem cells (SSCs) residing in the testis. SSCs represent less than 1% of the total testicular cell population in adults, but must achieve a stable balance between self-renewal and differentiation. Once differentiation has occurred, the newly formed and highly proliferative spermatogonia must then enter the meiotic program in which DNA content is doubled, then halved twice to create haploid gametes. While much is known about the critical cellular processes that take place during the specialized cell division that is meiosis, much less is known about how the spermatocytes in the "first-wave" in juveniles compare to those that contribute to long-term, "steady-state" spermatogenesis in adults. Given the strictly-defined developmental process of spermatogenesis, this study explored the transcriptional profiles of developmental cell stages during testis maturation. Using a combination of comprehensive germ cell sampling with high-resolution, single-cell-mRNA-sequencing, we have generated a reference dataset of germ cell gene expression. We show that discrete developmental stages of spermatogenesis possess significant differences in the transcriptional profiles from neonates compared to juveniles and adults. Importantly, these gene expression dynamics are also reflected at the protein level in their respective cell types. We also show differential utilization of many biological pathways with age in both spermatogonia and spermatocytes, demonstrating significantly different underlying gene regulatory programs in these cell types over the course of testis development and spermatogenic waves. This dataset represents the first unbiased sampling of spermatogonia and spermatocytes during testis maturation, at high-resolution, single-cell depth. Not only does this analysis reveal previously unknown transcriptional dynamics of a highly transitional cell population, it has also begun to reveal critical differences in biological pathway utilization in developing spermatogonia and spermatocytes, including response to DNA damage and double-strand breaks.

Transcriptome profiling of alphaherpesvirus-infected cells treated with the HIV-integrase inhibitor raltegravir reveals profound and specific alterations in host transcription

Anti-microbial compounds typically exert their action by directly interfering with one or more stages of the pathogen's life cycle. However, some compounds also have secondary effects on the host that aid in pathogen clearance. Raltegravir is a human immunodeficiency virus (HIV)-integrase inhibitor that has been shown to alter the host immune response to HIV in addition to its direct antiviral effect. Interestingly, raltegravir can also directly inhibit the replication of various herpesviruses. However, the host-targeted effects of this drug in the context of a herpesvirus infection have not been explored. Here, we used felid alphaherpesvirus 1 (FHV-1), a close relative of human alphaherpesvirus 1 (HHV-1) that similarly causes ocular herpes, to characterize the host-targeted effects of raltegravir on corneal epithelial cells during an alphaherpesvirus infection. Using RNA deep sequencing, we found that raltegravir specifically boosts the expression of anti-angiogenic factors and promotes metabolic homeostasis in FHV-1-infected cells. In contrast, few changes in host gene transcription were found in uninfected cells. Importantly, we were able to demonstrate that these effects were specific to raltegravir and independent of the direct-acting antiviral effect of the drug, since treatment with the DNA polymerase inhibitor phosphonoacetic acid did not induce these host-targeted effects. Taken together, these results indicate that raltegravir has profound and specific effects on the host transcription profile of herpesvirus-infected cells that may contribute to the overall antiviral activity of the drug and could provide therapeutic benefits in vivo. Furthermore, this study provides a framework for future efforts evaluating the host-targeted effects of anti-microbial compounds.

Temporal changes in macrophage phenotype after peripheral nerve injury

BACKGROUND

Macrophages play a key role in peripheral nerve repair and demonstrate complex phenotypes that are highly dependent on microenvironmental cues.

METHODS

We determined temporal changes in macrophage gene expression over time using RNA sequencing after fluorescence-activated cell sorting (FACS) macrophage populations from injured peripheral nerve. We identified key upstream regulators and dominant pathways using ingenuity pathway analysis and confirmed these changes with NanoString technology. We then investigate the effects of extreme polarizers of macrophage phenotype (IL4 and IFNγ) on nerve regeneration. We determined macrophage gene expression in vivo at the site of peripheral nerve injury with NanoString technology, and assessed recovery from sciatic nerve injury by cranial tibial muscle weights and retrograde labeling motor neurons in mice with deletion of IL4 or IFNγ receptors.

RESULTS

We demonstrate that IL4R and IFNγR deletions provide complementary responses to polarization, and alter expression of genes associated with angiogenesis and axonal extension, but do not influence recovery from peripheral nerve transection at 8 weeks after repair.

CONCLUSIONS

Overall, this study provides a framework to evaluate the phenotype of macrophages over time, and provides a broader and more precise assessment of gene expression changes than has previously been commonly used. This data suggests ways in which polarization may be modulated to improve repair.

Developmental Origin Governs CD8+ T Cell Fate Decisions during Infection

Heterogeneity is a hallmark feature of the adaptive immune system in vertebrates. Following infection, naive T cells differentiate into various subsets of effector and memory T cells, which help to eliminate pathogens and maintain long-term immunity. The current model suggests there is a single lineage of naive T cells that give rise to different populations of effector and memory T cells depending on the type and amounts of stimulation they encounter during infection. Here, we have discovered that multiple sub-populations of cells exist in the naive CD8⁺ T cell pool that are distinguished by their developmental origin, unique transcriptional profiles, distinct chromatin landscapes, and different kinetics and phenotypes after microbial challenge. These data demonstrate that the naive CD8⁺ T cell pool is not as homogeneous as previously thought and offers a new framework for explaining the remarkable heterogeneity in the effector and memory T cell subsets that arise after infection.

Maternal Choline Supplementation during Normal Murine Pregnancy Alters the Placental Epigenome: Results of an Exploratory Study

The placental epigenome regulates processes that affect placental and fetal development, and could be mediating some of the reported effects of maternal choline supplementation (MCS) on placental vascular development and nutrient delivery. As an extension of work previously conducted in pregnant mice, the current study sought to explore the effects of MCS on various epigenetic markers in the placenta. RNA and DNA were extracted from placentas collected on embryonic day 15.5 from pregnant mice fed a 1X or 4X choline diet, and were subjected to genome-wide sequencing procedures or mass-spectrometry-based assays to examine placental imprinted gene expression, DNA methylation patterns, and microRNA (miRNA) abundance. MCS yielded a higher (fold change = 1.63-2.25) expression of four imprinted genes (Ampd3, Tfpi2, Gatm and Aqp1) in the female placentas and a lower (fold change = 0.46-0.62) expression of three imprinted genes (Dcn, Qpct and Tnfrsf23) in the male placentas (false discovery rate (FDR) ≤ 0.05 for both sexes). Methylation in the promoter regions of these genes and global placental DNA methylation were also affected (p ≤ 0.05). Additionally, a lower (fold change = 0.3; P_unadjusted = 2.05 × 10^-4; FDR = 0.13) abundance of miR-2137 and a higher (fold change = 1.25-3.92; p < 0.05) expression of its target genes were detected in the 4X choline placentas. These data demonstrate that the placental epigenome is responsive to maternal choline intake during murine pregnancy and likely mediates some of the previously described choline-induced effects on placental and fetal outcomes.

Angiogenic factor imbalance precedes complement deposition in placentae of the BPH/5 model of preeclampsia

Preeclampsia (PE), a hypertensive disorder of pregnancy, is a leading cause of maternal and fetal morbidity and mortality. Although the etiology is unknown, PE is thought to be caused by defective implantation and decidualization in pregnancy. Pregnant blood pressure high (BPH)/5 mice spontaneously develop placentopathies and maternal features of human PE. We hypothesized that BPH/5 implantation sites have transcriptomic alterations. Next-generation RNA sequencing of implantation sites at peak decidualization, embryonic day (E)7.5, revealed complement gene up-regulation in BPH/5 vs. controls. In BPH/5, expression of complement factor 3 was increased around the decidual vasculature of E7.5 implantation sites and in the trophoblast giant cell layer of E10.5 placentae. Altered expression of VEGF pathway genes in E5.5 BPH/5 implantation sites preceded complement dysregulation, which correlated with abnormal vasculature and increased placental growth factor mRNA and VEGF₁₆₄ expression at E7.5. By E10.5, proangiogenic genes were down-regulated, whereas antiangiogenic sFlt-1 was up-regulated in BPH/5 placentae. We found that early local misexpression of VEGF genes and abnormal decidual vasculature preceded sFlt-1 overexpression and increased complement deposition in BPH/5 placentae. Our findings suggest that abnormal decidual angiogenesis precedes complement activation, which in turn contributes to the aberrant trophoblast invasion and poor placentation that underlie PE.-Sones, J. L., Merriam, A. A., Seffens, A., Brown-Grant, D.-A., Butler, S. D., Zhao, A. M., Xu, X., Shawber, C. J., Grenier, J. K., Douglas, N. C. Angiogenic factor imbalance precedes complement deposition in placentae of the BPH/5 model of preeclampsia.

Melanocyte Stem Cell Activation and Translocation Initiate Cutaneous Melanoma in Response to UV Exposure

Melanoma is one of the deadliest cancers, yet the cells of origin and mechanisms of tumor initiation remain unclear. The majority of melanomas emerge from clear skin without a precursor lesion, but it is unknown whether these melanomas can arise from melanocyte stem cells (MCSCs). Here we employ mouse models to define the role of MCSCs as melanoma cells of origin, demonstrate that MCSC quiescence acts as a tumor suppressor, and identify the extrinsic environmental and molecular factors required for the critical early steps of melanoma initiation. Specifically, melanomas originate from melanoma-competent MCSCs upon stimulation by UVB, which induces MCSC activation and translocation via an inflammation-dependent process. Moreover, the chromatin-remodeling factor Hmga2 in the skin plays a critical role in UVB-mediated melanomagenesis. These findings delineate melanoma formation from melanoma-competent MCSCs following extrinsic stimuli, and they suggest that abrogation of Hmga2 function in the microenvironment can suppress MCSC-originating cutaneous melanomas.

RNA-seq transcriptome analysis of formalin fixed, paraffin-embedded canine meningioma

Meningiomas are the most commonly reported primary intracranial tumor in dogs and humans and between the two species there are similarities in histology and biologic behavior. Due to these similarities, dogs have been proposed as models for meningioma pathobiology. However, little is known about specific pathways and individual genes that are involved in the development and progression of canine meningioma. In addition, studies are lacking that utilize RNAseq to characterize gene expression in clinical cases of canine meningioma. The primary objective of this study was to develop a technique for which high quality RNA can be extracted from formalin-fixed, paraffin embedded tissue and then used for transcriptome analysis to determine patterns of gene expression. RNA was extracted from thirteen canine meningiomas-eleven from formalin fixed and two flash-frozen. These represented six grade I and seven grade II meningiomas based on the World Health Organization classification system for human meningioma. RNA was also extracted from fresh frozen leptomeninges from three control dogs for comparison. RNAseq libraries made from formalin fixed tissue were of sufficient quality to successfully identify 125 significantly differentially expressed genes, the majority of which were related to oncogenic processes. Twelve genes (AQP1, BMPER, FBLN2, FRZB, MEDAG, MYC, PAMR1, PDGFRL, PDPN, PECAM1, PERP, ZC2HC1C) were validated using qPCR. Among the differentially expressed genes were oncogenes, tumor suppressors, transcription factors, VEGF-related genes, and members of the WNT pathway. Our work demonstrates that RNA of sufficient quality can be extracted from FFPE canine meningioma samples to provide biologically relevant transcriptome analyses using a next-generation sequencing technique, such as RNA-seq.

XX Disorder of Sex Development is associated with an insertion on chromosome 9 and downregulation of RSPO1 in dogs (Canis lupus familiaris)

Remarkable progress has been achieved in understanding the mechanisms controlling sex determination, yet the cause for many Disorders of Sex Development (DSD) remains unknown. Of particular interest is a rare XX DSD subtype in which individuals are negative for SRY, the testis determining factor on the Y chromosome, yet develop testes or ovotestes, and both of these phenotypes occur in the same family. This is a naturally occurring disorder in humans (Homo sapiens) and dogs (C. familiaris). Phenotypes in the canine XX DSD model are strikingly similar to those of the human XX DSD subtype. The purposes of this study were to identify 1) a variant associated with XX DSD in the canine model and 2) gene expression alterations in canine embryonic gonads that could be informative to causation. Using a genome wide association study (GWAS) and whole genome sequencing (WGS), we identified a variant on C. familiaris autosome 9 (CFA9) that is associated with XX DSD in the canine model and in affected purebred dogs. This is the first marker identified for inherited canine XX DSD. It lies upstream of SOX9 within the canine ortholog for the human disorder, which resides on 17q24. Inheritance of this variant indicates that XX DSD is a complex trait in which breed genetic background affects penetrance. Furthermore, the homozygous variant genotype is associated with embryonic lethality in at least one breed. Our analysis of gene expression studies (RNA-seq and PRO-seq) in embryonic gonads at risk of XX DSD from the canine model identified significant RSPO1 downregulation in comparison to XX controls, without significant upregulation of SOX9 or other known testis pathway genes. Based on these data, a novel mechanism is proposed in which molecular lesions acting upstream of RSPO1 induce epigenomic gonadal mosaicism.

Neonatal CD8+ T cells provide innate immune protection against unrelated pathogens

We previously found that neonatal and adult CD8+ T cells adopt different fates after infection with specific pathogens because they are derived from a distinct progenitor cell that is distinguished by expression of Lin28b. In this report, we investigated if CD8+ T cells respond differently to unrelated pathogens in early life. To start, we compared gene expression profiles of neonatal and adult CD8+ T cells after stimulation with innate cytokines (IL-12 + IL-18). Interestingly, we found that neonatal CD8+ T cells are more responsive and express a distinct cytokine profile compared to their adult counterparts. To examine the biological significance of these age-related differences, we next adoptively transferred monoclonal populations of neonatal or adult CD8+ T cells into T cell-deficient recipients and compared their ability to mediate immune protection against an irrelevant strain of Listeria monocytogenes. Surprisingly, we found that neonatal donor cell displayed increased proliferation and effector functions, and the pathogen burden was significantly reduced when compared to recipients containing adult donor cells. We also found that Lin28b may be responsible for the enhanced responsiveness to innate cytokines in early life, as ectopic expression of Lin28b in adult cells enabled them to initiate a program of bystander activation that was analogous to neonates. Collectively, these findings indicate that neonatal CD8+ T cells possess innate-like functions during infections at the cost of generating memory, and these unique behaviors may be driven by high levels of Lin28b. We currently are employing ATAC-seq to determine whether Lin28b changes the program of bystander activation by altering the chromatin landscape.

SIRT7 Is an RNA-Activated Protein Lysine Deacylase

Mammalian SIRT7 is a member of the sirtuin family that regulates multiple biological processes including genome stability, metabolic pathways, stress responses, and tumorigenesis. SIRT7 has been shown to be important for ribosome biogenesis and transcriptional regulation. SIRT7 knockout mice exhibit complications associated with fatty liver and increased aging in hematopoietic stem cells. However, the molecular basis for its biological function remains unclear, in part due to the lack of efficient enzymatic activity in vitro. Previously, we have demonstrated that double-stranded DNA could activate SIRT7's deacetylase activity in vitro, allowing it to deacetylate H3K18 in the context of chromatin. Here, we show that RNA can increase the catalytic efficiency of SIRT7 even better and that SIRT7 can remove long chain fatty acyl groups more efficiently than removing acetyl groups. Truncation and mutagenesis studies revealed residues at both the amino and carboxyl termini of SIRT7 that are involved in RNA-binding and important for activity. RNA immunoprecipitation-sequencing (RIP-seq) identified ribosomal RNA (rRNA) as the predominant RNA binding partner of SIRT7. The associated RNA was able to effectively activate the deacetylase and defatty-acylase activities of SIRT7. Knockdown of SIRT7 increased the lysine fatty acylation of several nuclear proteins based on metabolic labeling with an alkyne-tagged fatty acid analog, supporting that the defatty-acylase activity of SIRT7 is physiologically relevant. These findings provide important insights into the biological functions of SIRT7, as well as an improved platform to develop SIRT7 modulators.

Extensive local adaptation within the chemosensory system following Drosophila melanogaster's global expansion

How organisms adapt to new environments is of fundamental biological interest, but poorly understood at the genetic level. Chemosensory systems provide attractive models to address this problem, because they lie between external environmental signals and internal physiological responses. To investigate how selection has shaped the well-characterized chemosensory system of Drosophila melanogaster, we have analysed genome-wide data from five diverse populations. By couching population genomic analyses of chemosensory protein families within parallel analyses of other large families, we demonstrate that chemosensory proteins are not outliers for adaptive divergence between species. However, chemosensory families often display the strongest genome-wide signals of recent selection within D. melanogaster. We show that recent adaptation has operated almost exclusively on standing variation, and that patterns of adaptive mutations predict diverse effects on protein function. Finally, we provide evidence that chemosensory proteins have experienced relaxed constraint, and argue that this has been important for their rapid adaptation over short timescales.

Fruit flies gain valuable information about their environment by sensing chemicals. Here, Arguello et al. show strong signals of recent selection on the chemosensory system of the fruit fly Drosophila melanogaster, consistent with the adaptation of populations to their local chemical environment.

miRNA-embedded shRNAs for Lineage-specific BCL11A Knockdown and Hemoglobin F Induction

RNA interference (RNAi) technology using short hairpin RNAs (shRNAs) expressed via RNA polymerase (pol) III promoters has been widely exploited to modulate gene expression in a variety of mammalian cell types. For certain applications, such as lineage-specific knockdown, embedding targeting sequences into pol II-driven microRNA (miRNA) architecture is required. Here, using the potential therapeutic target BCL11A, we demonstrate that pol III-driven shRNAs lead to significantly increased knockdown but also increased cytotoxcity in comparison to pol II-driven miRNA adapted shRNAs (shRNA(miR)) in multiple hematopoietic cell lines. We show that the two expression systems yield mature guide strand sequences that differ by a 4 bp shift. This results in alternate seed sequences and consequently influences the efficacy of target gene knockdown. Incorporating a corresponding 4 bp shift into the guide strand of shRNA(miR)s resulted in improved knockdown efficiency of BCL11A. This was associated with a significant de-repression of the hemoglobin target of BCL11A, human γ-globin or the murine homolog Hbb-y. Our results suggest the requirement for optimization of shRNA sequences upon incorporation into a miRNA backbone. These findings have important implications in future design of shRNA(miR)s for RNAi-based therapy in hemoglobinopathies and other diseases requiring lineage-specific expression of gene silencing sequences.

Behavioral idiosyncrasy reveals genetic control of phenotypic variability

Quantitative genetics has primarily focused on describing genetic effects on trait means and largely ignored the effect of alternative alleles on trait variability, potentially missing an important axis of genetic variation contributing to phenotypic differences among individuals. To study the genetic effects on individual-to-individual phenotypic variability (or intragenotypic variability), we used Drosophila inbred lines and measured the spontaneous locomotor behavior of flies walking individually in Y-shaped mazes, focusing on variability in locomotor handedness, an assay optimized to measure variability. We discovered that some lines had consistently high levels of intragenotypic variability among individuals, whereas lines with low variability behaved as although they tossed a coin at each left/right turn decision. We demonstrate that the degree of variability is itself heritable. Using a genome-wide association study (GWAS) for the degree of intragenotypic variability as the phenotype across lines, we identified several genes expressed in the brain that affect variability in handedness without affecting the mean. One of these genes, Ten-a, implicates a neuropil in the central complex of the fly brain as influencing the magnitude of behavioral variability, a brain region involved in sensory integration and locomotor coordination. We validated these results using genetic deficiencies, null alleles, and inducible RNAi transgenes. Our study reveals the constellation of phenotypes that can arise from a single genotype and shows that different genetic backgrounds differ dramatically in their propensity for phenotypic variabililty. Because traditional mean-focused GWASs ignore the contribution of variability to overall phenotypic variation, current methods may miss important links between genotype and phenotype.

Dgcr8 and Dicer are essential for sex chromosome integrity during meiosis in males

Small RNAs play crucial roles in regulating gene expression during mammalian meiosis. To investigate the function of microRNAs (miRNAs) and small interfering RNAs (siRNAs) during meiosis in males, we generated germ-cell-specific conditional deletions of Dgcr8 and Dicer in mice. Analysis of spermatocytes from both conditional knockout lines revealed that there were frequent chromosomal fusions during meiosis, always involving one or both sex chromosomes. RNA sequencing indicates upregulation of Atm in spermatocytes from miRNA-deficient mice, and immunofluorescence imaging demonstrates an increased abundance of activated ATM kinase and mislocalization of phosphorylated MDC1, an ATM phosphorylation substrate. The Atm 3′UTR contains many potential microRNA target sites, and, notably, target sites for several miRNAs depleted in both conditional knockout mice were highly effective at promoting repression. RNF8, a telomere-associated protein whose localization is controlled by the MDC1–ATM kinase cascade, normally associates with the sex chromosomes during pachytene, but in both conditional knockouts redistributed to the autosomes. Taken together, these results suggest that Atm dysregulation in microRNA-deficient germ lines contributes to the redistribution of proteins involved in chromosomal stability from the sex chromosomes to the autosomes, resulting in sex chromosome fusions during meiotic prophase I.

Highlighted Article: miRNA-deficient spermatocytes display frequent sex chromosome fusions and fail to progress through meiosis in a process that is probably mediated by dysregulation of Atm.

Receptor tyrosine kinases fall into distinct classes based on their inferred signaling networks

Although many anticancer drugs that target receptor tyrosine kinases (RTKs) provide clinical benefit, their long-term use is limited by resistance that is often attributed to increased abundance or activation of another RTK that compensates for the inhibited receptor. To uncover common and unique features in the signaling networks of RTKs, we measured time-dependent signaling in six isogenic cell lines, each expressing a different RTK as downstream proteins were systematically perturbed by RNA interference. Network models inferred from the data revealed a conserved set of signaling pathways and RTK-specific features that grouped the RTKs into three distinct classes: (i) an EGFR/FGFR1/c-Met class constituting epidermal growth factor receptor, fibroblast growth factor receptor 1, and the hepatocyte growth factor receptor c-Met; (ii) an IGF-1R/NTRK2 class constituting insulin-like growth factor 1 receptor and neurotrophic tyrosine receptor kinase 2; and (iii) a PDGFRβ class constituting platelet-derived growth factor receptor β. Analysis of cancer cell line data showed that many RTKs of the same class were coexpressed and that increased abundance of an RTK or its cognate ligand frequently correlated with resistance to a drug targeting another RTK of the same class. In contrast, abundance of an RTK or ligand of one class generally did not affect sensitivity to a drug targeting an RTK of a different class. Thus, classifying RTKs by their inferred networks and then therapeutically targeting multiple receptors within a class may delay or prevent the onset of resistance.

A comprehensive platform for highly multiplexed mammalian functional genetic screens

Background

Genome-wide screening in human and mouse cells using RNA interference and open reading frame over-expression libraries is rapidly becoming a viable experimental approach for many research labs. There are a variety of gene expression modulation libraries commercially available, however, detailed and validated protocols as well as the reagents necessary for deconvolving genome-scale gene screens using these libraries are lacking. As a solution, we designed a comprehensive platform for highly multiplexed functional genetic screens in human, mouse and yeast cells using popular, commercially available gene modulation libraries. The Gene Modulation Array Platform (GMAP) is a single microarray-based detection solution for deconvolution of loss and gain-of-function pooled screens.

Results

Experiments with specially constructed lentiviral-based plasmid pools containing ~78,000 shRNAs demonstrated that the GMAP is capable of deconvolving genome-wide shRNA "dropout" screens. Further experiments with a larger, ~90,000 shRNA pool demonstrate that equivalent results are obtained from plasmid pools and from genomic DNA derived from lentivirus infected cells. Parallel testing of large shRNA pools using GMAP and next-generation sequencing methods revealed that the two methods provide valid and complementary approaches to deconvolution of genome-wide shRNA screens. Additional experiments demonstrated that GMAP is equivalent to similar microarray-based products when used for deconvolution of open reading frame over-expression screens.

Conclusion

Herein, we demonstrate four major applications for the GMAP resource, including deconvolution of pooled RNAi screens in cells with at least 90,000 distinct shRNAs. We also provide detailed methodologies for pooled shRNA screen readout using GMAP and compare next-generation sequencing to GMAP (i.e. microarray) based deconvolution methods.

Unbiased reconstruction of a mammalian transcriptional network mediating pathogen responses

Models of mammalian regulatory networks controlling gene expression have been inferred from genomic data but have largely not been validated. We present an unbiased strategy to systematically perturb candidate regulators and monitor cellular transcriptional responses. We applied this approach to derive regulatory networks that control the transcriptional response of mouse primary dendritic cells to pathogens. Our approach revealed the regulatory functions of 125 transcription factors, chromatin modifiers, and RNA binding proteins, which enabled the construction of a network model consisting of 24 core regulators and 76 fine-tuners that help to explain how pathogen-sensing pathways achieve specificity. This study establishes a broadly applicable, comprehensive, and unbiased approach to reveal the wiring and functions of a regulatory network controlling a major transcriptional response in primary mammalian cells.

Integrative genomic approaches identify IKBKE as a breast cancer oncogene

The karyotypic chaos exhibited by human epithelial cancers complicates efforts to identify mutations critical for malignant transformation. Here we integrate complementary genomic approaches to identify human oncogenes. We show that activation of the ERK and phosphatidylinositol 3-kinase (PI3K) signaling pathways cooperate to transform human cells. Using a library of activated kinases, we identify several kinases that replace PI3K signaling and render cells tumorigenic. Whole genome structural analyses reveal that one of these kinases, IKBKE (IKKepsilon), is amplified and overexpressed in breast cancer cell lines and patient-derived tumors. Suppression of IKKepsilon expression in breast cancer cell lines that harbor IKBKE amplifications induces cell death. IKKepsilon activates the nuclear factor-kappaB (NF-kappaB) pathway in both cell lines and breast cancers. These observations suggest a mechanism for NF-kappaB activation in breast cancer, implicate the NF-kappaB pathway as a downstream mediator of PI3K, and provide a framework for integrated genomic approaches in oncogene discovery.

A lentiviral RNAi library for human and mouse genes applied to an arrayed viral high-content screen

To enable arrayed or pooled loss-of-function screens in a wide range of mammalian cell types, including primary and nondividing cells, we are developing lentiviral short hairpin RNA (shRNA) libraries targeting the human and murine genomes. The libraries currently contain 104,000 vectors, targeting each of 22,000 human and mouse genes with multiple sequence-verified constructs. To test the utility of the library for arrayed screens, we developed a screen based on high-content imaging to identify genes required for mitotic progression in human cancer cells and applied it to an arrayed set of 5,000 unique shRNA-expressing lentiviruses that target 1,028 human genes. The screen identified several known and approximately 100 candidate regulators of mitotic progression and proliferation; the availability of multiple shRNAs targeting the same gene facilitated functional validation of putative hits. This work provides a widely applicable resource for loss-of-function screens, as well as a roadmap for its application to biological discovery.

Exploiting the enzymatic recognition of an unnatural base pair to develop a universal genetic analysis system

BACKGROUND

With the invention of the DNA chip, genome-wide analysis is now a reality. Unfortunately, solid-phase detection systems such as the DNA chip suffer from a narrow range in quantification and sensitivity. Today the best methodology for sensitive, wide dynamic range quantification and genotyping of nucleic acids is real-time PCR. However, multiplexed real-time PCR technologies require complicated and costly design and manufacturing of separate detection probes for each new target.

METHODS

We developed a novel real-time PCR technology that uses universal energy transfer probes constructed from An Expanded Genetic Information System (AEGIS) for both quantification and genotyping analyses.

RESULTS

RNA quantification by reverse transcription-PCR was linear over four orders of magnitude for the simultaneous analysis of beta-actin messenger RNA and 18S ribosomal RNA. A single trial validation study of 176 previously genotyped clinical specimens was performed by endpoint analysis for factor V Leiden and prothrombin 20210A mutation detection. There was concordance for 173 samples between the genotyping results from Invader tests and the AEGIS universal energy transfer probe system for both factor V Leiden and prothrombin G20210A. Two prothrombin and one factor V sample gave indeterminate results (no calls).

CONCLUSION

The AEGIS universal probe system allows for rapid development of PCR assays for nucleic acid quantification and genotyping.

Functional evolution of the Ultrabithorax protein

The Hox genes have been implicated as central to the evolution of animal body plan diversity. Regulatory changes both in Hox expression domains and in Hox-regulated gene networks have arisen during the evolution of related taxa, but there is little knowledge of whether functional changes in Hox proteins have also contributed to morphological evolution. For example, the evolution of greater numbers of differentiated segments and body parts in insects, compared with the simpler body plans of arthropod ancestors, may have involved an increase in the spectrum of biochemical interactions of individual Hox proteins. Here, we compare the in vivo functions of orthologous Ultrabithorax (Ubx) proteins from the insect Drosophila melanogaster and from an onychophoran, a member of a sister phylum with a more primitive and homonomous body plan. These Ubx proteins, which have been diverging in sequence for over 540 million years, can generate many of the same gain-of-function tissue transformations and can activate and repress many of the same target genes when expressed during Drosophila development. However, the onychophora Ubx (OUbx) protein does not transform the segmental identity of the embryonic ectoderm or repress the Distal-less target gene. This functional divergence is due to sequence changes outside the conserved homeodomain region. The inability of OUbx to function like Drosophila Ubx (DUbx) in the embryonic ectoderm indicates that the Ubx protein may have acquired new cofactors or activity modifiers since the divergence of the onychophoran and insect lineages.

Hox genes in brachiopods and priapulids and protostome evolution

Understanding the early evolution of animal body plans requires knowledge both of metazoan phylogeny and of the genetic and developmental changes involved in the emergence of particular forms. Recent 18S ribosomal RNA phylogenies suggest a three-branched tree for the Bilateria comprising the deuterostomes and two great protostome clades, the lophotrochozoans and ecdysozoans. Here, we show that the complement of Hox genes in critical protostome phyla reflects these phylogenetic relationships and reveals the early evolution of developmental regulatory potential in bilaterians. We have identified Hox genes that are shared by subsets of protostome phyla. These include a diverged pair of posterior (Abdominal-B-like) genes in both a brachiopod and a polychaete annelid, which supports the lophotrochozoan assemblage, and a distinct posterior Hox gene shared by a priapulid, a nematode and the arthropods, which supports the ecdysozoan clade. The ancestors of each of these two major protostome lineages had a minimum of eight to ten Hox genes. The major period of Hox gene expansion and diversification thus occurred before the radiation of each of the three great bilaterian clades.

Evolution of the entire arthropod Hox gene set predated the origin and radiation of the onychophoran/arthropod clade

BACKGROUND

Dramatic changes in body size and pattern occurred during the radiation of many taxa in the Cambrian, and these changes are best documented for the arthropods. The sudden appearance of such diverse body plans raises the fundamental question of when the genes and the developmental control systems that regulate these designs evolved. As Hox genes regulate arthropod body patterns, the evolution of these genes may have played a role in the origin and diversification of the arthropod body plan from a homonomous ancestor. To trace the origin of arthropod Hox genes, we examined their distribution in a myriapod and in the Onychophora, a sister group to the arthropods.

RESULTS

Despite the limited segmental diversity within myriapods and Onychophora, all insect Hox genes are present in both taxa, including the trunk Hox genes Ultrabithorax and abdominal-A as well as an ortholog of the fushi tarazu gene. Comparative analysis of Hox gene deployment revealed that the anterior boundary of expression of trunk Hox genes has shifted dramatically along the anteroposterior axis between Onychophora and different arthropod classes. Furthermore, we found that repression of expression of the Hox target gene Distal-less is unique to the insect lineage.

CONCLUSIONS

A complete arthropod Hox gene family existed in the ancestor of the onychophoran/arthropod clade. No new Hox genes were therefore required to catalyze the arthropod radiation; instead, arthropod body-plan diversity arose through changes in the regulation of Hox genes and their downstream targets.