Neonatal CD8+ T Cells Resist Exhaustion during Chronic Infection

Chronic viral infections, such as HIV and hepatitis C virus, represent a major public health problem. Although it is well understood that neonates and adults respond differently to chronic viral infections, the underlying mechanisms remain unknown. In this study, we transferred neonatal and adult CD8+ T cells into a mouse model of chronic infection (lymphocytic choriomeningitis virus clone 13) and dissected out the key cell-intrinsic differences that alter their ability to protect the host. Interestingly, we found that neonatal CD8+ T cells preferentially became effector cells early in chronic infection compared with adult CD8+ T cells and expressed higher levels of genes associated with cell migration and effector cell differentiation. During the chronic phase of infection, the neonatal cells retained more immune functionality and expressed lower levels of surface markers and genes related to exhaustion. Because the neonatal cells protect from viral replication early in chronic infection, the altered differentiation trajectories of neonatal and adult CD8+ T cells is functionally significant. Together, our work demonstrates how cell-intrinsic differences between neonatal and adult CD8+ T cells influence key cell fate decisions during chronic infection.

The gene regulatory basis of bystander activation in CD8+ T cells

CD8+ T cells are classically recognized as adaptive lymphocytes based on their ability to recognize specific foreign antigens and mount memory responses. However, recent studies indicate that some antigen-inexperienced CD8+ T cells can respond to innate cytokines alone in the absence of cognate T cell receptor stimulation, a phenomenon referred to as bystander activation. Here, we demonstrate that neonatal CD8+ T cells undergo a robust and diverse program of bystander activation, which corresponds to enhanced innate-like protection against unrelated pathogens. Using a multi-omics approach, we found that the ability of neonatal CD8+ T cells to respond to innate cytokines derives from their capacity to undergo rapid chromatin remodeling, resulting in the usage of a distinct set of enhancers and transcription factors typically found in innate-like T cells. We observed that the switch between innate and adaptive functions in the CD8+ T cell compartment is mediated by changes in the abundance of distinct subsets of cells. The innate CD8+ T cell subset that predominates in early life was also present in adult mice and humans. Our findings provide support for the layered immune hypothesis and indicate that the CD8+ T cell compartment is more functionally diverse than previously thought.

Single-cell transcriptomics of the immune system in ME/CFS at baseline and following symptom provocation

Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) is a serious and poorly understood disease. To understand immune dysregulation in ME/CFS, we use single-cell RNA sequencing (scRNA-seq) to examine immune cells in patient and control cohorts. Postexertional malaise (PEM), an exacerbation of symptoms following strenuous exercise, is a characteristic symptom of ME/CFS. To detect changes coincident with PEM, we applied scRNA-seq on the same cohorts following exercise. At baseline, ME/CFS patients display classical monocyte dysregulation suggestive of inappropriate differentiation and migration to tissue. We identify both diseased and more normal monocytes within patients, and the fraction of diseased cells correlates with disease severity. Comparing the transcriptome at baseline and postexercise challenge, we discover patterns indicative of improper platelet activation in patients, with minimal changes elsewhere in the immune system. Taken together, these data identify immunological defects present at baseline in patients and an additional layer of dysregulation in platelets.

HIV-1 active and latent infections induce disparate chromatin reorganization and transcriptional regulation of mRNAs and lncRNAs in SupT1 cells

IMPORTANCE

HIV-1 infection of T-lymphocytes depends on co-opting cellular transcriptional and translational machineries for viral replication. This requires significant changes in the cellular microenvironment. We have characterized and compared the changes in cellular chromatin structures as well as gene expression landscapes in T cells that are either actively or latently infected with HIV-1. Our results reveal that chromatin accessibility and expression of both protein-coding mRNAs and non-coding lncRNAs are uniquely regulated in HIV-1-infected T cells, depending on whether the virus is actively transcribing or remains in a transcriptionally silent, latent state. HIV-1 latent infection elicits more robust changes in the cellular chromatin organization than active viral infection. Our analysis also identifies the effects of such epigenomic changes on the cellular gene expression and subsequent biological pathways. This study comprehensively characterizes the cellular epigenomic and transcriptomic states that support active and latent HIV-1 infection in an in vitro model of SupT1 cells.

Detection and characterization of novel luchacoviruses, genus Alphacoronavirus, in saliva and feces of meso-carnivores in the northeastern United States

IMPORTANCE

Several coronaviruses (CoVs) have been detected in domesticated, farmed, and wild meso-carnivores, causing a wide range of diseases and infecting diverse species, highlighting their important but understudied role in the epidemiology of these viruses. Assessing the viral diversity hosted in wildlife species is essential to understand their significance in the cross-species transmission of CoVs. Our focus here was on CoV discovery in meso-carnivores in the Northeast United States as a potential "hotspot" area with high density of humans and urban wildlife. This study identifies novel alphacoronaviruses circulating in multiple free-ranging wild and domestic species in this area and explores their potential epidemiological importance based on regions of the Spike gene, which are relevant for virus-host interactions.

Polyomavirus-associated Disseminated T-cell Lymphoma in a Colony of Zebra Finches (Taeniopygia guttata)

Four zebra finches in a closed research colony presented with variable clinical signs, including masses, skin lesions, shivering, and/or ruffled feathers. These birds were not responsive to treatment efforts; 3 died and one was euthanized. All 4 were submitted for necropsy to determine the cause of the clinical signs. Gross necropsy and histopathologic findings from all birds resulted in a diagnosis of round cell neoplasia in multiple organs, including the skin, liver, kidney, and reproductive tract, with intranuclear inclusion bodies in the neoplastic cells. In all 4 cases, immunohistochemical staining showed strong immunoreactivity for CD3 in 70% to 80% of the neoplastic round cells, with a relatively small subset that were immunopositive for Pax5. These findings supported a diagnosis of T-cell lymphoma. Frozen liver tissue from one case was submitted for next-generation sequencing (NGS), which revealed viral RNA with 100% sequence homology to canary polyomavirus strain 34639 that had originally been identified in a European goldfinch. Formalin-fixed paraffin-embedded scrolls from another case were also submitted for NGS, which revealed viral RNA with 97.2% sequence homology to canary polyomavirus strain 37273 that had originally been identified in a canary. To localize the virus in situ, RNAscope hybridization was performed using a probe designed to target the VP1 gene of the sequenced virus in frozen liver tissue. In all 4 cases, disseminated and robust hybridization signals were detected in neoplastic cells. These findings indicate that polyomaviruses have the potential to be oncogenic in zebra finches.

Detection and characterization of novel luchacoviruses, genus Alphacoronavirus, shed in saliva and feces of meso-carnivores in the northeastern United States

Small to mid-sized carnivores, or meso-carnivores, comprise a group of diverse mammals, many of which can adapt to anthropogenically disturbed environments. Wild meso-carnivores living in urban areas may get exposed to or spread pathogens to other species, including stray/feral domestic animals. Several coronaviruses (CoVs) have been detected in domesticated and farmed meso-carnivores, but knowledge of CoVs circulating in free-ranging wild meso-carnivores remains limited. In this study, we analyzed 321 samples collected between 2016 and 2022 from 9 species of free-ranging wild meso-carnivores and stray/feral domestic cats in the northeastern United States. Using a pan-CoV PCR, we screened tissues, feces, and saliva, nasal, and rectal swabs. We detected CoV RNA in fecal and saliva samples of animals in four species: fisher (Pekania pennanti), bobcat (Lynx rufus), red fox (Vulpes vulpes), and domestic cat (Felis catus). Next-generation sequencing revealed that all these viruses belonged to the Luchacovirus subgenus (Alphacoronavirus genus), previously reported only in rodents and lagomorphs (i.e., rabbits). Genetic comparison of the 3'-end of the genome (~12,000bp) revealed that although the viruses detected group with, and have a genetic organization similar to other luchacoviruses, they are genetically distinct from those from rodents and lagomorphs. Genetic characterization of the spike protein revealed that the meso-carnivore luchacoviruses do not have an S1/S2 cleavage motif but do have highly variable structural loops containing cleavage motifs similar to those identified in certain pathogenic CoVs. This study highlights the importance of characterizing the spike protein of CoVs in wild species for further targeted epidemiologic monitoring.

Polyomavirus-associated Disseminated T-cell Lymphoma in a Colony of Zebra Finches (Taeniopygia guttata)

Four zebra finches in a closed research colony presented with variable clinical signs, including masses, skin lesions,shivering, and/or ruffled feathers. These birds were not responsive to treatment efforts; 3 died and one was euthanized. All4 were submitted for necropsy to determine the cause of the clinical signs. Gross necropsy and histopathologic findings fromall birds resulted in a diagnosis of round cell neoplasia in multiple organs, including the skin, liver, kidney, and reproductivetract, with intranuclear inclusion bodies in the neoplastic cells. In all 4 cases, immunohistochemical staining showed strongimmunoreactivity for CD3 in 70% to 80% of the neoplastic round cells, with a relatively small subset that were immunopositivefor Pax5. These findings supported a diagnosis of T-cell lymphoma. Frozen liver tissue from one case was submittedfor next-generation sequencing (NGS), which revealed viral RNA with 100% sequence homology to canary polyomavirusstrain 34639 that had originally been identified in a European goldfinch. Formalin-fixed paraffin-embedded scrolls fromanother case were also submitted for NGS, which revealed viral RNA with 97.2% sequence homology to canary polyomavirusstrain 37273 that had originally been identified in a canary. To localize the virus in situ, RNAscope hybridizationwas performed using a probe designed to target the VP1 gene of the sequenced virus in frozen liver tissue. In all 4 cases,disseminated and robust hybridization signals were detected in neoplastic cells. These findings indicate that polyomaviruseshave the potential to be oncogenic in zebra finches.

Transcriptomic signatures of feline chronic gingivostomatitis are influenced by upregulated IL6

Feline chronic gingivostomatitis (FCGS) is a relatively common and debilitating disease characterized by bilateral inflammation and ulceration of the caudal oral mucosa, alveolar and buccal mucosa, and varying degrees of periodontal disease. The etiopathogenesis of FCGS remains unresolved. In this study, we performed bulk RNA-seq molecular profiling of affected tissues derived from a cohort of client-owned cats with FCGS compared to tissues from unaffected animals, to identify candidate genes and pathways that can help guide future exploration of novel clinical solutions. We complemented transcriptomic findings with immunohistochemistry and in situ hybridization assays to better understand the biological significance of the results and performed RNA-seq validation of biologically relevant differentially expressed genes using qPCR assays to demonstrate technical reproducibility. Transcriptomic profiles of oral mucosal tissues in cats with FCGS are enriched with immune- and inflammation-related genes and pathways that appear to be largely influenced by IL6, and include NFKB, JAK/STAT, IL-17 and IFN type I and II signaling, offering new opportunities to develop novel clinical applications based on a more rational understanding of the disease.

Genomic landscape and gene expression profiles of feline oral squamous cell carcinoma

Feline oral squamous cell carcinoma (FOSCC) is a cancer of the squamous cell lining in the oral cavity and represents up to 80% of all oral cancers in cats, with a poor prognosis. We have used whole exome sequencing (WES) and RNA sequencing of the tumor to discover somatic mutations and gene expression changes that may be associated with FOSCC occurrence. FOSCC offers a potential comparative model to study human head and neck squamous cell carcinoma (HNSCC) due to its similar spontaneous formation, and morphological and histological features. In this first study using WES to identify somatic mutations in feline cancer, we have identified tumor-associated gene mutations in six cats with FOSCC and found some overlap with identified recurrently mutated genes observed in HNSCC. Four samples each had mutations in TP53, a common mutation in all cancers, but each was unique. Mutations in other cellular growth control genes were also found such as KAT2B and ARID1A. Enrichment analysis of FOSCC gene expression profiles suggests a molecular similarity to human OSCC as well, including alterations in epithelial to mesenchymal transition and IL6/JAK/STAT pathways. In this preliminary study, we present exome and transcriptome results that further our understanding of FOSCC.

Transcriptomic signatures of feline chronic gingivostomatitis are influenced by upregulated IL6

Feline chronic gingivostomatitis (FCGS) is a relatively common and debilitating disease characterized by bilateral inflammation and ulceration of the caudal oral mucosa, alveolar and buccal mucosa, and varying degrees of periodontal disease. The etiopathogenesis of FCGS remains unresolved. In this study, we performed bulk RNA-seq molecular profiling of affected tissues derived from a cohort of client-owned cats with FCGS compared to tissues from unaffected animals, to identify candidate genes and pathways that can help guide future exploration of novel clinical solutions. We complemented transcriptomic findings with immunohistochemistry and in situ hybridization assays to better understand the biological significance of the results and performed RNA-seq validation of selected differentially expressed genes using qPCR assays to demonstrate technical reproducibility. Transcriptomic profiles of oral mucosal tissues in cats with FCGS are enriched with immune- and inflammation-related genes and pathways that appear to be largely influenced by IL6 , and include NFKB, JAK/STAT, IL-17 and IFN type I and II signaling, offering new opportunities to develop novel clinical applications based on a more rational understanding of the disease.

Highly recurrent BRAF p.V595E mutation in canine papillary oral squamous cell carcinoma

Oral squamous cell carcinoma (OSCC) is the most common oral epithelial malignancy in dogs. It exhibits locally aggressive biological behaviour with the potential to metastasize, and a reported 1-year survival rate of 0% when left untreated. Expression studies suggest that aberrant MAPK signalling plays a key role in canine OSCC tumorigenesis, which is consistent with BRAF and HRAS MAPK-activating mutations reported in some tumours. Several morphological subtypes of canine OSCC have been described, with papillary, conventional, and basaloid as the most common patterns. We hypothesized that mutational differences may underlie these phenotypic variations. In this study, targeted Sanger sequencing and restriction fragment length polymorphism assays demonstrate that up to 85.7% of canine papillary OSCC (n = 14) harbour a BRAF p.V595E mutation. Assessment of neoplastic epithelial cell proliferation using Ki67 immunolabelling (n = 10) confirmed a relatively high proliferation activity, consistent with their known aggressive clinical behaviour. These findings underscore a consistent genetic feature of canine papillary OSCC and provide a basis for the development of novel diagnostic and targeted therapeutic approaches that can improve the quality of veterinary care.

Delay modulates the immune response to nerve repair

Effective regeneration after peripheral nerve injury requires macrophage recruitment. We investigated the activation of remodeling pathways within the macrophage population when repair is delayed and identified alteration of key upstream regulators of the inflammatory response. We then targeted one of these regulators, using exogenous IL10 to manipulate the response to injury at the repair site. We demonstrate that this approach alters macrophage polarization, promotes macrophage recruitment, axon extension, neuromuscular junction formation, and increases the number of regenerating motor units reaching their target. We also demonstrate that this approach can rescue the effects of delayed nerve graft.

Postnatal oogenesis leads to an exceptionally large ovarian reserve in naked mole-rats

In the long-lived naked mole-rat (NMR), the entire process of oogenesis occurs postnatally. Germ cell numbers increase significantly in NMRs between postnatal days 5 (P5) and P8, and germs cells positive for proliferation markers (Ki-67, pHH3) are present at least until P90. Using pluripotency markers (SOX2 and OCT4) and the primordial germ cell (PGC) marker BLIMP1, we show that PGCs persist up to P90 alongside germ cells in all stages of female differentiation and undergo mitosis both in vivo and in vitro. We identified VASA+ SOX2+ cells at 6 months and at 3-years in subordinate and reproductively activated females. Reproductive activation was associated with proliferation of VASA+ SOX2+ cells. Collectively, our results suggest that highly desynchronized germ cell development and the maintenance of a small population of PGCs that can expand upon reproductive activation are unique strategies that could help to maintain the NMR's ovarian reserve for its 30-year reproductive lifespan.

Delayed puberty, gonadotropin abnormalities and subfertility in male Padi2/Padi4 double knockout mice

BACKGROUND

Peptidylarginine deiminase enzymes (PADs) convert arginine residues to citrulline in a process called citrullination or deimination. Recently, two PADs, PAD2 and PAD4, have been linked to hormone signaling in vitro and the goal of this study was to test for links between PAD2/PAD4 and hormone signaling in vivo.

METHODS

Preliminary analysis of Padi2 and Padi4 single knockout (SKO) mice did not find any overt reproductive defects and we predicted that this was likely due to genetic compensation. To test this hypothesis, we created a Padi2/Padi4 double knockout (DKO) mouse model and tested these mice along with wild-type FVB/NJ (WT) and both strains of SKO mice for a range of reproductive defects.

RESULTS

Controlled breeding trials found that male DKO mice appeared to take longer to have their first litter than WT controls. This tendency was maintained when these mice were mated to either DKO or WT females. Additionally, unsexed 2-day old DKO pups and male DKO weanlings both weighed significantly less than their WT counterparts, took significantly longer than WT males to reach puberty, and had consistently lower serum testosterone levels. Furthermore, 90-day old adult DKO males had smaller testes than WT males with increased rates of germ cell apoptosis.

CONCLUSIONS

The Padi2/Padi4 DKO mouse model provides a new tool for investigating PAD function and outcomes from our studies provide the first in vivo evidence linking PADs with hormone signaling.

Chemical, Molecular, and Single-nucleus Analysis Reveal Chondroitin Sulfate Proteoglycan Aberrancy in Fibrolamellar Carcinoma

Fibrolamellar carcinoma (FLC) is an aggressive liver cancer with no effective therapeutic options. The extracellular environment of FLC tumors is poorly characterized and may contribute to cancer growth and/or metastasis. To bridge this knowledge gap, we assessed pathways relevant to proteoglycans, a major component of the extracellular matrix. We first analyzed gene expression data from FLC and nonmalignant liver tissue (n = 27) to identify changes in glycosaminoglycan (GAG) biosynthesis pathways and found that genes associated with production of chondroitin sulfate, but not other GAGs, are significantly increased by 8-fold. We then implemented a novel LC/MS-MS based method to quantify the abundance of different types of GAGs in patient tumors (n = 16) and found that chondroitin sulfate is significantly more abundant in FLC tumors by 6-fold. Upon further analysis of GAG-associated proteins, we found that versican (VCAN) expression is significantly upregulated at the mRNA and protein levels, the latter of which was validated by IHC. Finally, we performed single-cell assay for transposase-accessible chromatin sequencing on FLC tumors (n = 3), which revealed for the first time the different cell types in FLC tumors and also showed that VCAN is likely produced not only from FLC tumor epithelial cells but also activated stellate cells. Our results reveal a pathologic aberrancy in chondroitin (but not heparan) sulfate proteoglycans in FLC and highlight a potential role for activated stellate cells.

Significance

This study leverages a multi-disciplinary approach, including state-of-the-art chemical analyses and cutting-edge single-cell genomic technologies, to identify for the first time a marked chondroitin sulfate aberrancy in FLC that could open novel therapeutic avenues in the future.

Comparative transcriptional profiling of canine acanthomatous ameloblastoma and homology with human ameloblastoma

Ameloblastomas are odontogenic tumors that are rare in people but have a relatively high prevalence in dogs. Because canine acanthomatous ameloblastomas (CAA) have clinicopathologic and molecular features in common with human ameloblastomas (AM), spontaneous CAA can serve as a useful translational model of disease. However, the molecular basis of CAA and how it compares to AM are incompletely understood. In this study, we compared the global genomic expression profile of CAA with AM and evaluated its dental origin by using a bulk RNA-seq approach. For these studies, healthy gingiva and canine oral squamous cell carcinoma served as controls. We found that aberrant RAS signaling, and activation of the epithelial-to-mesenchymal transition cellular program are involved in the pathogenesis of CAA, and that CAA is enriched with genes known to be upregulated in AM including those expressed during the early stages of tooth development, suggesting a high level of molecular homology. These results support the model that domestic dogs with spontaneous CAA have potential for pre-clinical assessment of targeted therapeutic modalities against AM.

Multisite Evaluation of Next-Generation Methods for Small RNA Quantification

Small RNAs (smRNAs) are important regulators of many biologic processes and are now most frequently characterized using Illumina sequencing. However, although standard RNA sequencing library preparation has become routine in most sequencing facilities, smRNA sequencing library preparation has historically been challenging because of high input requirements, laborious protocols involving gel purifications, inability to automate, and a lack of benchmarking standards. Additionally, studies have suggested that many of these methods are nonlinear and do not accurately reflect the amounts of smRNAs in vivo. Recently, a number of new kits have become available that permit lower input amounts and less laborious, gel-free protocol options. Several of these new kits claim to reduce RNA ligase-dependent sequence bias through novel adapter modifications and to lessen adapter-dimer contamination in the resulting libraries. With the increasing number of smRNA kits available, understanding the relative strengths of each method is crucial for appropriate experimental design. In this study, we systematically compared 9 commercially available smRNA library preparation kits as well as NanoString probe hybridization across multiple study sites. Although several of the new methodologies do reduce the amount of artificially over- and underrepresented microRNAs (miRNAs), we observed that none of the methods was able to remove all of the bias in the library preparation. Identical samples prepared with different methods show highly varied levels of different miRNAs. Even so, many methods excelled in ease of use, lower input requirement, fraction of usable reads, and reproducibility across sites. These differences may help users select the most appropriate methods for their specific question of interest.

Dual RNA-Seq of Mtb-Infected Macrophages In Vivo Reveals Ontologically Distinct Host-Pathogen Interactions

Dissecting the in vivo host-pathogen interplay is crucial to understanding the molecular mechanisms governing control or progression of intracellular infections. In this work, we explore the in vivo molecular dynamics of Mtb infection by performing dual RNA-seq on Mycobacterium tuberculosis-infected, ontogenetically distinct macrophage lineages isolated directly from murine lungs. We first define an in vivo signature of 180 genes specifically upregulated by Mtb in mouse lung macrophages, then we uncover a divergent transcriptional response of the bacteria between alveolar macrophages that appear to sustain Mtb growth through increased access to iron and fatty acids and interstitial macrophages that restrict Mtb growth through iron sequestration and higher levels of nitric oxide. We use an enrichment protocol for bacterial transcripts, which enables us to probe Mtb physiology at the host cell level in an in vivo environment, with broader application in understanding the infection dynamics of intracellular pathogens in general.

In this study Pisu et al. performed dual RNA-seq on Mycobacterium tuberculosis-infected, ontogenetically distinct macrophage lineages isolated directly from infected murine lungs. The transcriptional response of host and bacteria diverged between alveolar macrophages that sustain Mtb growth and interstitial macrophages that restrict Mtb growth.

Dual RNA-Sequencing of Mycobacterium tuberculosis-Infected Cells from a Murine Infection Model

Dual RNA-sequencing is a powerful technique to assess both bacterial and host transcriptomes in an unbiased way. We developed a protocol to perform Dual RNA-seq on in vivo-derived macrophage populations infected with Mycobacterium tuberculosis. Here, we provide a practical step-by-step guide to execute the protocol on Mtb-infected cells from a murine infection model. Our protocol can also be easily applied to perform Dual RNA-seq on in vitro-derived cells as well as different Mtb-infected host cell types.

For complete details on the use and execution of this protocol, please refer to Pisu et al. (2020)

•

Dual RNA-seq of M. tuberculosis-infected cells

•

Simultaneous profiling of host and pathogen transcripts from an in vivo infection model

•

Gain in-depth understanding of host-pathogen interactions at the molecular level

CDK2 kinase activity is a regulator of male germ cell fate

The ability of men to remain fertile throughout their lives depends upon establishment of a spermatogonial stem cell (SSC) pool from gonocyte progenitors, and thereafter balancing SSC renewal versus terminal differentiation. Here, we report that precise regulation of the cell cycle is crucial for this balance. Whereas cyclin-dependent kinase 2 (Cdk2) is not necessary for mouse viability or gametogenesis stages prior to meiotic prophase I, mice bearing a deregulated allele (Cdk2Y15S ) are severely deficient in spermatogonial differentiation. This allele disrupts an inhibitory phosphorylation site (Tyr15) for the kinase WEE1. Remarkably, Cdk2Y15S/Y15S mice possess abnormal clusters of mitotically active SSC-like cells, but these are eventually removed by apoptosis after failing to differentiate properly. Analyses of lineage markers, germ cell proliferation over time, and single cell RNA-seq data revealed delayed and defective differentiation of gonocytes into SSCs. Biochemical and genetic data demonstrated that Cdk2Y15S is a gain-of-function allele causing elevated kinase activity, which underlies these differentiation defects. Our results demonstrate that precise regulation of CDK2 kinase activity in male germ cell development is crucial for the gonocyte-to-spermatogonia transition and long-term spermatogenic homeostasis.

Ultra-frequent HRAS p.Q61R somatic mutation in canine acanthomatous ameloblastoma reveals pathogenic similarities with human ameloblastoma

Ameloblastoma is a locally aggressive odontogenic tumour that occurs in humans and dogs. Most ameloblastomas (AM) in humans harbour mutually-exclusive driving mutations in BRAF, HRAS, KRAS, NRAS or FGFR2 that activate MAPK signalling, and in SMO that activates Hedgehog signalling. The remarkable clinical and histological similarities between canine acanthomatous ameloblastoma (CAA) and AM suggest they may harbour similar driving mutations. In this study, aimed at characterizing the mutational status of SMO, BRAF, HRAS, KRAS, NRAS and FGFR2 in CAA, we used RNA sequencing, Sanger sequencing and restriction fragment length polymorphism assays to demonstrate that 94% of CAA (n = 16) harbour a somatic HRAS p.Q61R mutation. The similarities in MAPK-activating mutational profiles between CAA and AM implicate conserved molecular mechanisms of tumorigenesis, thus, qualifying the dog as a potentially useful model of disease. Given the relevance of RAS mutations in the pathogenesis of odontogenic tumours and other types of cancer, the results of this study are of comparative, translational, and veterinary value.

The Notch signaling pathway promotes basophil responses during helminth-induced type 2 inflammation

Basophils promote type 2 inflammation that mediates worm clearance during murine infection with the gastrointestinal helminth parasite Trichuris muris. Webb et al. show for the first time that basophil–intrinsic Notch signaling is required for basophil gene expression and a functional program that support helminth expulsion.

Type 2 inflammation drives the clearance of gastrointestinal helminth parasites, which infect over two billion people worldwide. Basophils are innate immune cells that support host-protective type 2 inflammation during murine infection with the helminth Trichuris muris. However, the mechanisms required for basophil function and gene expression regulation in this context remain unclear. We show that during T. muris infection, basophils localized to the intestine and up-regulated Notch receptor expression, rendering them sensitive to Notch signals that rapidly regulate gene expression programs. In vitro, Notch inhibition limited basophil cytokine production in response to cytokine stimulation. Basophil-intrinsic Notch signaling was required for T. muris–elicited changes in genome-wide basophil transcriptional programs. Mice lacking basophil-intrinsic functional Notch signaling had impaired worm clearance, decreased intestinal type 2 inflammation, altered basophil localization in the intestine, and decreased CD4⁺ T helper 2 cell responses following infection. These findings demonstrate that Notch is required for basophil gene expression and effector function associated with helminth expulsion during type 2 inflammation.

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.