Investigation of the host transcriptional response to intracellular bacterial infection using Dictyostelium discoideum as a host model.

BACKGROUND: During infection by intracellular pathogens, a highly complex interplay occurs between the infected cell trying to degrade the invader and the pathogen which actively manipulates the host cell to enable survival and proliferation. Many intracellular pathogens pose important threats to human health and major efforts have been undertaken to better understand the host-pathogen interactions that eventually determine the outcome of the infection. Over the last decades, the unicellular eukaryote Dictyostelium discoideum has become an established infection model, serving as a surrogate macrophage that can be infected with a wide range of intracellular pathogens. In this study, we use high-throughput RNA-sequencing to analyze the transcriptional response of D. discoideum when infected with Mycobacterium marinum and Legionella pneumophila. The results were compared to available data from human macrophages. RESULTS: The majority of the transcriptional regulation triggered by the two pathogens was found to be unique for each bacterial challenge. Hallmark transcriptional signatures were identified for each infection, e.g. induction of endosomal sorting complexes required for transport (ESCRT) and autophagy genes in response to M. marinum and inhibition of genes associated with the translation machinery and energy metabolism in response to L. pneumophila. However, a common response to the pathogenic bacteria was also identified, which was not induced by non-pathogenic food bacteria. Finally, comparison with available data sets of regulation in human monocyte derived macrophages shows that the elicited response in D. discoideum is in many aspects similar to what has been observed in human immune cells in response to Mycobacterium tuberculosis and L. pneumophila. CONCLUSIONS: Our study presents high-throughput characterization of D. discoideum transcriptional response to intracellular pathogens using RNA-seq. We demonstrate that the transcriptional response is in essence distinct to each pathogen and that in many cases, the corresponding regulation is recapitulated in human macrophages after infection by mycobacteria and L. pneumophila. This indicates that host-pathogen interactions are evolutionary conserved, derived from the early interactions between free-living phagocytic cells and bacteria. Taken together, our results strengthen the use of D. discoideum as a general infection model.

Comparative omics and feeding manipulations in chicken indicate a shift of the endocrine role of visceral fat towards reproduction.

BACKGROUND: The mammalian adipose tissue plays a central role in energy-balance control, whereas the avian visceral fat hardly expresses leptin, the key adipokine in mammals. Therefore, to assess the endocrine role of adipose tissue in birds, we compared the transcriptome and proteome between two metabolically different types of chickens, broilers and layers, bred towards efficient meat and egg production, respectively. RESULTS: Broilers and layer hens, grown up to sexual maturation under free-feeding conditions, differed 4.0-fold in weight and 1.6-fold in ovarian-follicle counts, yet the relative accumulation of visceral fat was comparable. RNA-seq and mass-spectrometry (MS) analyses of visceral fat revealed differentially expressed genes between broilers and layers, 1106 at the mRNA level (FDR ≤ 0.05), and 203 at the protein level (P ≤ 0.05). In broilers, Ingenuity Pathway Analysis revealed activation of the PTEN-pathway, and in layers increased response to external signals. The expression pattern of genes encoding fat-secreted proteins in broilers and layers was characterized in the RNA-seq and MS data, as well as by qPCR on visceral fat under free feeding and 24 h-feed deprivation. This characterization was expanded using available RNA-seq data of tissues from red junglefowl, and of visceral fat from broilers of different types. These comparisons revealed expression of new adipokines and secreted proteins (LCAT, LECT2, SERPINE2, SFTP1, ZP1, ZP3, APOV1, VTG1 and VTG2) at the mRNA and/or protein levels, with dynamic gene expression patterns in the selected chicken lines (except for ZP1; FDR/P ≤ 0.05) and feed deprivation (NAMPT, SFTPA1 and ZP3) (P ≤ 0.05). In contrast, some of the most prominent adipokines in mammals, leptin, TNF, IFNG, and IL6 were expressed at a low level (FPKM/RPKM< 1) and did not show differential mRNA expression neither between broiler and layer lines nor between fed vs. feed-deprived chickens. CONCLUSIONS: Our study revealed that RNA and protein expression in visceral fat changes with selective breeding, suggesting endocrine roles of visceral fat in the selected phenotypes. In comparison to gene expression in visceral fat of mammals, our findings points to a more direct cross talk of the chicken visceral fat with the reproductive system and lower involvement in the regulation of appetite, inflammation and insulin resistance.

Reading and editing the Pleurodeles waltl genome reveals novel features of tetrapod regeneration.

Salamanders exhibit an extraordinary ability among vertebrates to regenerate complex body parts. However, scarce genomic resources have limited our understanding of regeneration in adult salamanders. Here, we present the ~20 Gb genome and transcriptome of the Iberian ribbed newt Pleurodeles waltl, a tractable species suitable for laboratory research. We find that embryonic stem cell-specific miRNAs mir-93b and mir-427/430/302, as well as Harbinger DNA transposons carrying the Myb-like proto-oncogene have expanded dramatically in the Pleurodeles waltl genome and are co-expressed during limb regeneration. Moreover, we find that a family of salamander methyltransferases is expressed specifically in adult appendages. Using CRISPR/Cas9 technology to perturb transcription factors, we demonstrate that, unlike the axolotl, Pax3 is present and necessary for development and that contrary to mammals, muscle regeneration is normal without functional Pax7 gene. Our data provide a foundation for comparative genomic studies that generate models for the uneven distribution of regenerative capacities among vertebrates.

Evolution of the p53-MDM2 pathway.

BACKGROUND: The p53 signalling pathway, which controls cell fate, has been extensively studied due to its prominent role in tumor development. The pathway includes the tumor supressor protein p53, its vertebrate paralogs p63 and p73, and their negative regulators MDM2 and MDM4. The p53/p63/p73-MDM system is ancient and can be traced in all extant animal phyla. Despite this, correct phylogenetic trees including both vertebrate and invertebrate species of the p53/p63/p73 and MDM families have not been published. RESULTS: Here, we have examined the evolution of the p53/p63/p73 protein family with particular focus on the p53/p63/p73 transactivation domain (TAD) and its co-evolution with the p53/p63/p73-binding domain (p53/p63/p73BD) of MDM2. We found that the TAD and p53/p63/p73BD share a strong evolutionary connection. If one of the domains of the protein is lost in a phylum, then it seems very likely to be followed by loss of function by the other domain as well, and due to the loss of function it is likely to eventually disappear. By focusing our phylogenetic analysis to p53/p63/p73 and MDM proteins from phyla that retain the interaction domains TAD and p53/p63/p73BD, we built phylogenetic trees of p53/p63/p73 and MDM based on both vertebrate and invertebrate species. The trees follow species evolution and contain a total number of 183 and 98 species for p53/p63/p73 and MDM, respectively. We also demonstrate that the p53/p63/p73 and MDM families result from whole genome duplications. CONCLUSIONS: The signaling pathway of the TAD and p53/p63/p73BD in p53/p63/p73 and MDM, respectively, dates back to early metazoan time and has since then tightly co-evolved, or disappeared in distinct lineages.

PiiL: visualization of DNA methylation and gene expression data in gene pathways.

BACKGROUND: DNA methylation is a major mechanism involved in the epigenetic state of a cell. It has been observed that the methylation status of certain CpG sites close to or within a gene can directly affect its expression, either by silencing or, in some cases, up-regulating transcription. However, a vertebrate genome contains millions of CpG sites, all of which are potential targets for methylation, and the specific effects of most sites have not been characterized to date. To study the complex interplay between methylation status, cellular programs, and the resulting phenotypes, we present PiiL, an interactive gene expression pathway browser, facilitating analyses through an integrated view of methylation and expression on multiple levels. RESULTS: PiiL allows for specific hypothesis testing by quickly assessing pathways or gene networks, where the data is projected onto pathways that can be downloaded directly from the online KEGG database. PiiL provides a comprehensive set of analysis features that allow for quick and specific pattern searches. Individual CpG sites and their impact on host gene expression, as well as the impact on other genes present in the regulatory network, can be examined. To exemplify the power of this approach, we analyzed two types of brain tumors, Glioblastoma multiform and lower grade gliomas. CONCLUSION: At a glance, we could confirm earlier findings that the predominant methylation and expression patterns separate perfectly by mutations in the IDH genes, rather than by histology. We could also infer the IDH mutation status for samples for which the genotype was not known. By applying different filtering methods, we show that a subset of CpG sites exhibits consistent methylation patterns, and that the status of sites affect the expression of key regulator genes, as well as other genes located downstream in the same pathways. PiiL is implemented in Java with focus on a user-friendly graphical interface. The source code is available under the GPL license from https://github.com/behroozt/PiiL.git .

Correspondence on Lovell et al.: identification of chicken genes previously assumed to be evolutionarily lost.

Through RNA-Seq analyses, we identified 137 genes that are missing in chicken, including the long-sought-after nephrin and tumor necrosis factor genes. These genes tended to cluster in GC-rich regions that have poor coverage in genome sequence databases. Hence, the occurrence of syntenic groups of vertebrate genes that have not been observed in Aves does not prove the evolutionary loss of such genes.Please see related Research article: http://dx.doi.org/10.1186/s13059-014-0565-1 and Please see response from Lovell et al: https://www.dx.doi.org/10.1186/s13059-017-1234-y.

Combinatorial identification of DNA methylation patterns over age in the human brain.

BACKGROUND: DNA methylation plays a key role in developmental processes, which is reflected in changing methylation patterns at specific CpG sites over the lifetime of an individual. The underlying mechanisms are complex and possibly affect multiple genes or entire pathways. RESULTS: We applied a multivariate approach to identify combinations of CpG sites that undergo modifications when transitioning between developmental stages. Monte Carlo feature selection produced a list of ranked and statistically significant CpG sites, while rule-based models allowed for identifying particular methylation changes in these sites. Our rule-based classifier reports combinations of CpG sites, together with changes in their methylation status in the form of easy-to-read IF-THEN rules, which allows for identification of the genes associated with the underlying sites. CONCLUSION: We utilized machine learning and statistical methods to discretize decision class (age) values to get a general pattern of methylation changes over the lifespan. The CpG sites present in the significant rules were annotated to genes involved in brain formation, general development, as well as genes linked to cancer and Alzheimer's disease.

Coordinated Changes in Gene Expression Throughout Encystation of Giardia intestinalis.

Differentiation into infectious cysts through the process of encystation is crucial for transmission and survival of the intestinal protozoan parasite Giardia intestinalis. Hitherto the majority of studies have focused on the early events, leaving late encystation poorly defined. In order to further study encystation, focusing on the later events, we developed a new encystation protocol that generates a higher yield of mature cysts compared to standard methods. Transcriptome changes during the entire differentiation from trophozoites to cysts were thereafter studied using RNA sequencing (RNA-seq). A high level of periodicity was observed for up- and down-regulated genes, both at the level of the entire transcriptome and putative regulators. This suggests the trajectory of differentiation to be coordinated through developmentally linked gene regulatory activities. Our study identifies a core of 13 genes that are consistently up-regulated during initial encystation. Of these, two constitute previously uncharacterized proteins that we were able to localize to a new type of encystation-specific vesicles. Interestingly, the largest transcriptional changes were seen in the late phase of encystation with the majority of the highly up-regulated genes encoding hypothetical proteins. Several of these were epitope-tagged and localized to further characterize these previously unknown genetic components of encystation and possibly excystation. Finally, we also detected a switch of variant specific surface proteins (VSPs) in the late phase of encystation. This occurred at the same time as nuclear division and DNA replication, suggesting a potential link between the processes.

Broad-scale phylogenomics provides insights into retrovirus-host evolution.

Genomic data provide an excellent resource to improve understanding of retrovirus evolution and the complex relationships among viruses and their hosts. In conjunction with broad-scale in silico screening of vertebrate genomes, this resource offers an opportunity to complement data on the evolution and frequency of past retroviral spread and so evaluate future risks and limitations for horizontal transmission between different host species. Here, we develop a methodology for extracting phylogenetic signal from large endogenous retrovirus (ERV) datasets by collapsing information to facilitate broad-scale phylogenomics across a wide sample of hosts. Starting with nearly 90,000 ERVs from 60 vertebrate host genomes, we construct phylogenetic hypotheses and draw inferences regarding the designation, host distribution, origin, and transmission of the Gammaretrovirus genus and associated class I ERVs. Our results uncover remarkable depths in retroviral sequence diversity, supported within a phylogenetic context. This finding suggests that current infectious exogenous retrovirus diversity may be underestimated, adding credence to the possibility that many additional exogenous retroviruses may remain to be discovered in vertebrate taxa. We demonstrate a history of frequent horizontal interorder transmissions from a rodent reservoir and suggest that rats may have acted as important overlooked facilitators of gammaretrovirus spread across diverse mammalian hosts. Together, these results demonstrate the promise of the methodology used here to analyze large ERV datasets and improve understanding of retroviral evolution and diversity for utilization in wider applications.

De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis.

De novo assembly of RNA-seq data enables researchers to study transcriptomes without the need for a genome sequence; this approach can be usefully applied, for instance, in research on 'non-model organisms' of ecological and evolutionary importance, cancer samples or the microbiome. In this protocol we describe the use of the Trinity platform for de novo transcriptome assembly from RNA-seq data in non-model organisms. We also present Trinity-supported companion utilities for downstream applications, including RSEM for transcript abundance estimation, R/Bioconductor packages for identifying differentially expressed transcripts across samples and approaches to identify protein-coding genes. In the procedure, we provide a workflow for genome-independent transcriptome analysis leveraging the Trinity platform. The software, documentation and demonstrations are freely available from http://trinityrnaseq.sourceforge.net. The run time of this protocol is highly dependent on the size and complexity of data to be analyzed. The example data set analyzed in the procedure detailed herein can be processed in less than 5 h.

Artificially designed promoters: understanding the role of spatial features and canonical binding sites in transcription.

The promoter is a key element in gene transcription and regulation. We previously reported that artificial sequences rich in the dinucleotide CpG are sufficient to drive expression in vitro in mammalian cell lines, without requiring canonical binding sites for transcription factor proteins. Here, we report that introducing a promoter organization that alternates in CpGs and regions rich in A and T further increases expression strength, as well as how insertion of specific binding sites makes such sequences respond to induced levels of the transcription factor NFκB. Our findings further contribute to the mechanistic understanding of promoters, as well as how these sequences might be shaped by evolutionary pressure in living organisms.

Comparative genomics reveals mobile pathogenicity chromosomes in Fusarium.

Fusarium species are among the most important phytopathogenic and toxigenic fungi. To understand the molecular underpinnings of pathogenicity in the genus Fusarium, we compared the genomes of three phenotypically diverse species: Fusarium graminearum, Fusarium verticillioides and Fusarium oxysporum f. sp. lycopersici. Our analysis revealed lineage-specific (LS) genomic regions in F. oxysporum that include four entire chromosomes and account for more than one-quarter of the genome. LS regions are rich in transposons and genes with distinct evolutionary profiles but related to pathogenicity, indicative of horizontal acquisition. Experimentally, we demonstrate the transfer of two LS chromosomes between strains of F. oxysporum, converting a non-pathogenic strain into a pathogen. Transfer of LS chromosomes between otherwise genetically isolated strains explains the polyphyletic origin of host specificity and the emergence of new pathogenic lineages in F. oxysporum. These findings put the evolution of fungal pathogenicity into a new perspective.

High-throughput sequencing of microdissected chromosomal regions.

The linkage of disease gene mapping with DNA sequencing is an essential strategy for defining the genetic basis of a disease. New massively parallel sequencing procedures will greatly facilitate this process, although enrichment for the target region before sequencing remains necessary. For this step, various DNA capture approaches have been described that rely on sequence-defined probe sets. To avoid making assumptions on the sequences present in the targeted region, we accessed specific cytogenetic regions in preparation for next-generation sequencing. We directly microdissected the target region in metaphase chromosomes, amplified it by degenerate oligonucleotide-primed PCR, and obtained sufficient material of high quality for high-throughput sequencing. Sequence reads could be obtained from as few as six chromosomal fragments. The power of cytogenetic enrichment followed by next-generation sequencing is that it does not depend on earlier knowledge of sequences in the region being studied. Accordingly, this method is uniquely suited for situations in which the sequence of a reference region of the genome is not available, including population-specific or tumor rearrangements, as well as previously unsequenced genomic regions such as centromeres.

A cis-acting regulatory mutation causes premature hair graying and susceptibility to melanoma in the horse.

In horses, graying with age is an autosomal dominant trait associated with a high incidence of melanoma and vitiligo-like depigmentation. Here we show that the Gray phenotype is caused by a 4.6-kb duplication in intron 6 of STX17 (syntaxin-17) that constitutes a cis-acting regulatory mutation. Both STX17 and the neighboring NR4A3 gene are overexpressed in melanomas from Gray horses. Gray horses carrying a loss-of-function mutation in ASIP (agouti signaling protein) had a higher incidence of melanoma, implying that increased melanocortin-1 receptor signaling promotes melanoma development in Gray horses. The Gray horse provides a notable example of how humans have cherry-picked mutations with favorable phenotypic effects in domestic animals.