Multiomics analysis identifies novel facilitators of human dopaminergic neuron differentiation.

Midbrain dopaminergic neurons (mDANs) control voluntary movement, cognition, and reward behavior under physiological conditions and are implicated in human diseases such as Parkinson's disease (PD). Many transcription factors (TFs) controlling human mDAN differentiation during development have been described, but much of the regulatory landscape remains undefined. Using a tyrosine hydroxylase (TH) human iPSC reporter line, we here generate time series transcriptomic and epigenomic profiles of purified mDANs during differentiation. Integrative analysis predicts novel regulators of mDAN differentiation and super-enhancers are used to identify key TFs. We find LBX1, NHLH1 and NR2F1/2 to promote mDAN differentiation and show that overexpression of either LBX1 or NHLH1 can also improve mDAN specification. A more detailed investigation of TF targets reveals that NHLH1 promotes the induction of neuronal miR-124, LBX1 regulates cholesterol biosynthesis, and NR2F1/2 controls neuronal activity.

Early Pleistocene enamel proteome from Dmanisi resolves Stephanorhinus phylogeny.

The sequencing of ancient DNA has enabled the reconstruction of speciation, migration and admixture events for extinct taxa1. However, the irreversible post-mortem degradation2 of ancient DNA has so far limited its recovery-outside permafrost areas-to specimens that are not older than approximately 0.5 million years (Myr)3. By contrast, tandem mass spectrometry has enabled the sequencing of approximately 1.5-Myr-old collagen type I4, and suggested the presence of protein residues in fossils of the Cretaceous period5-although with limited phylogenetic use6. In the absence of molecular evidence, the speciation of several extinct species of the Early and Middle Pleistocene epoch remains contentious. Here we address the phylogenetic relationships of the Eurasian Rhinocerotidae of the Pleistocene epoch7-9, using the proteome of dental enamel from a Stephanorhinus tooth that is approximately 1.77-Myr old, recovered from the archaeological site of Dmanisi (South Caucasus, Georgia)10. Molecular phylogenetic analyses place this Stephanorhinus as a sister group to the clade formed by the woolly rhinoceros (Coelodonta antiquitatis) and Merck's rhinoceros (Stephanorhinus kirchbergensis). We show that Coelodonta evolved from an early Stephanorhinus lineage, and that this latter genus includes at least two distinct evolutionary lines. The genus Stephanorhinus is therefore currently paraphyletic, and its systematic revision is needed. We demonstrate that sequencing the proteome of Early Pleistocene dental enamel overcomes the limitations of phylogenetic inference based on ancient collagen or DNA. Our approach also provides additional information about the sex and taxonomic assignment of other specimens from Dmanisi. Our findings reveal that proteomic investigation of ancient dental enamel-which is the hardest tissue in vertebrates11, and is highly abundant in the fossil record-can push the reconstruction of molecular evolution further back into the Early Pleistocene epoch, beyond the currently known limits of ancient DNA preservation.

PARK7/DJ-1 deficiency impairs microglial activation in response to LPS-induced inflammation.

BACKGROUND: Specific microglia responses are thought to contribute to the development and progression of neurodegenerative diseases, including Parkinson's disease (PD). However, the phenotypic acquisition of microglial cells and their role during the underlying neuroinflammatory processes remain largely elusive. Here, according to the multiple-hit hypothesis, which stipulates that PD etiology is determined by a combination of genetics and various environmental risk factors, we investigate microglial transcriptional programs and morphological adaptations under PARK7/DJ-1 deficiency, a genetic cause of PD, during lipopolysaccharide (LPS)-induced inflammation. METHODS: Using a combination of single-cell RNA-sequencing, bulk RNA-sequencing, multicolor flow cytometry and immunofluorescence analyses, we comprehensively compared microglial cell phenotypic characteristics in PARK7/DJ-1 knock-out (KO) with wildtype littermate mice following 6- or 24-h intraperitoneal injection with LPS. For translational perspectives, we conducted corresponding analyses in human PARK7/DJ-1 mutant induced pluripotent stem cell (iPSC)-derived microglia and murine bone marrow-derived macrophages (BMDMs). RESULTS: By excluding the contribution of other immune brain resident and peripheral cells, we show that microglia acutely isolated from PARK7/DJ-1 KO mice display a distinct phenotype, specially related to type II interferon and DNA damage response signaling, when compared with wildtype microglia, in response to LPS. We also detected discrete signatures in human PARK7/DJ-1 mutant iPSC-derived microglia and BMDMs from PARK7/DJ-1 KO mice. These specific transcriptional signatures were reflected at the morphological level, with microglia in LPS-treated PARK7/DJ-1 KO mice showing a less amoeboid cell shape compared to wildtype mice, both at 6 and 24 h after acute inflammation, as also observed in BMDMs. CONCLUSIONS: Taken together, our results show that, under inflammatory conditions, PARK7/DJ-1 deficiency skews microglia towards a distinct phenotype characterized by downregulation of genes involved in type II interferon signaling and a less prominent amoeboid morphology compared to wildtype microglia. These findings suggest that the underlying oxidative stress associated with the lack of PARK7/DJ-1 affects microglia neuroinflammatory responses, which may play a causative role in PD onset and progression.

Phagocytosis-related NADPH oxidase 2 subunit gp91phox contributes to neurodegeneration after repeated systemic challenge with lipopolysaccharides.

Repeated systemic challenge with lipopolysaccharides (LPS) can induce microglia activation and inflammatory neurodegeneration in the substantia nigra pars compacta region of mice. We now explored the role of mononuclear phagocytes associated nicotinamide adenine dinucleotide phosphate oxidase 2 (NOX-2) in inflammatory neurodegeneration. Cybb-deficient NOX-2 knock-out (KO) and control wild type (WT) mice were treated intraperitoneally daily over four consecutive days with 1 µg/gbw/day LPS. Transcriptome analysis by RNA-seq of total brain tissue indicated increased LPS-induced upregulation of genes belonging to the reactive oxygen species and reactive nitrogen species production, complement and lysosome activation as well as apoptosis and necroptosis in WT compared to NOX-2 KO mice. Validation of up-regulated gene transcripts via qRT-PCR confirmed that LPS-challenged NOX-2 KO mice expressed lower levels of the microglial phagocytosis-related genes Nos2, Cd68, Aif1/Iba1, Cyba, Itgam, and Fcer1g compared to WT mice at Day 5 after systemic inflammatory challenge, but no significant differences in the pro-inflammatory genes Tnfα and Il1b as well as microglial IBA1 and CD68 intensities were observed between both genotypes. Furthermore, loss of tyrosine hydroxylase positive (TH+) and NeuN positive neurons in the substantia nigra pars compacta upon repeated systemic LPS application were attenuated in NOX-2 KO mice. Thus, our data demonstrate that loss of dopaminergic neurons in the substantia nigra pars compacta after repeated systemic challenge with LPS is associated with a microglial phagocytosis-related gene activation profile involving the NADPH oxidase subunit Cybb/gp91phox.

Temporal enhancer profiling of parallel lineages identifies AHR and GLIS1 as regulators of mesenchymal multipotency.

Temporal data on gene expression and context-specific open chromatin states can improve identification of key transcription factors (TFs) and the gene regulatory networks (GRNs) controlling cellular differentiation. However, their integration remains challenging. Here, we delineate a general approach for data-driven and unbiased identification of key TFs and dynamic GRNs, called EPIC-DREM. We generated time-series transcriptomic and epigenomic profiles during differentiation of mouse multipotent bone marrow stromal cell line (ST2) toward adipocytes and osteoblasts. Using our novel approach we constructed time-resolved GRNs for both lineages and identifed the shared TFs involved in both differentiation processes. To take an alternative approach to prioritize the identified shared regulators, we mapped dynamic super-enhancers in both lineages and associated them to target genes with correlated expression profiles. The combination of the two approaches identified aryl hydrocarbon receptor (AHR) and Glis family zinc finger 1 (GLIS1) as mesenchymal key TFs controlled by dynamic cell type-specific super-enhancers that become repressed in both lineages. AHR and GLIS1 control differentiation-induced genes and their overexpression can inhibit the lineage commitment of the multipotent bone marrow-derived ST2 cells.

TREM2 triggers microglial density and age-related neuronal loss.

The microglial triggering receptor expressed on myeloid cells 2 (TREM2) signals via the activatory membrane adaptor molecule TYROBP. Genetic variants or mutations of TREM2 or TYROBP have been linked to inflammatory neurodegenerative diseases associated with aging. The typical aging process goes along with microglial changes and mild neuronal loss, but the exact contribution of TREM2 is still unclear. Aged TREM2 knock-out mice showed decreased age-related neuronal loss in the substantia nigra and the hippocampus. Transcriptomic analysis of the brains of 24 months old TREM2 knock-out mice revealed 211 differentially expressed genes mostly downregulated and associated with complement activation and oxidative stress response pathways. Consistently, 24 months old TREM2 knock-out mice showed lower transcription of microglial (Aif1 and Tmem119), oxidative stress markers (Inos, Cyba, and Cybb) and complement components (C1qa, C1qb, C1qc, C3, C4b, Itgam, and Itgb2), decreased microglial numbers and expression of the microglial activation marker Cd68, as well as accumulation of oxidized lipids. Cultured microglia of TREM2 knock-out mice showed reduced phagocytosis and oxidative burst. Thus, microglial TREM2 contributes to age-related microglial changes, phagocytic oxidative burst, and loss of neurons with possible detrimental effects during physiological aging.

Recalibrating Equus evolution using the genome sequence of an early Middle Pleistocene horse.

The rich fossil record of equids has made them a model for evolutionary processes. Here we present a 1.12-times coverage draft genome from a horse bone recovered from permafrost dated to approximately 560-780 thousand years before present (kyr BP). Our data represent the oldest full genome sequence determined so far by almost an order of magnitude. For comparison, we sequenced the genome of a Late Pleistocene horse (43 kyr BP), and modern genomes of five domestic horse breeds (Equus ferus caballus), a Przewalski's horse (E. f. przewalskii) and a donkey (E. asinus). Our analyses suggest that the Equus lineage giving rise to all contemporary horses, zebras and donkeys originated 4.0-4.5 million years before present (Myr BP), twice the conventionally accepted time to the most recent common ancestor of the genus Equus. We also find that horse population size fluctuated multiple times over the past 2 Myr, particularly during periods of severe climatic changes. We estimate that the Przewalski's and domestic horse populations diverged 38-72 kyr BP, and find no evidence of recent admixture between the domestic horse breeds and the Przewalski's horse investigated. This supports the contention that Przewalski's horses represent the last surviving wild horse population. We find similar levels of genetic variation among Przewalski's and domestic populations, indicating that the former are genetically viable and worthy of conservation efforts. We also find evidence for continuous selection on the immune system and olfaction throughout horse evolution. Finally, we identify 29 genomic regions among horse breeds that deviate from neutrality and show low levels of genetic variation compared to the Przewalski's horse. Such regions could correspond to loci selected early during domestication.

Crosstalk between different family members: IL27 recapitulates IFNγ responses in HCC cells, but is inhibited by IL6-type cytokines.

Interleukin-27 (IL27) is a type-I-cytokine of the IL6/IL12 family predominantly secreted by activated macrophages and dendritic cells. In the liver, IL27 expression was observed to be upregulated in patients with hepatitis B, and sera of hepatocellular carcinoma (HCC) patients contain significantly elevated levels of IL27 compared to healthy controls or patients with hepatitis and/or liver cirrhosis. In this study, we show that IL27 induces STAT1 and STAT3 phosphorylation in 5 HCC lines and 3 different types of non-transformed liver cells. We were especially interested in the relevance of the IL27-induced STAT3 activation in liver cells. Thus, we compared the IL27 responses with those induced by IFNγ (STAT1-dominated response) or IL6-type cytokines (IL6, hyper-IL6 (hy-IL6) or OSM) (STAT3-dominated response) by microarray analysis and find that in HCC cells, IL27 induces an IFNγ-like, STAT1-dependent transcriptional response, but we do not find an effective STAT3-dependent response. Validation experiments corroborate the finding from the microarray evaluation. Interestingly, the availability of STAT1 seems critical in the shaping of the IL27 response, as the siRNA knock-down of STAT1 revealed the ability of IL27 to induce the acute-phase protein γ-fibrinogen, a typical IL6 family characteristic. Moreover, we describe a crosstalk between the signaling of IL6-type cytokines and IL27: responses to the gp130-engaging cytokine IL27 (but not those to IFNs) can be inhibited by IL6-type cytokine pre-stimulation, likely by a SOCS3-mediated mechanism. Thus, IL27 recapitulates IFNγ responses in liver cells, but differs from IFNγ by its sensitivity to SOCS3 inhibition.

Loss of Myosin Vb in colorectal cancer is a strong prognostic factor for disease recurrence.

BACKGROUND: Selecting the most beneficial treatment regimens for colorectal cancer (CRC) patients remains challenging due to a lack of prognostic markers. Members of the Myosin family, proteins recognised to have a major role in trafficking and polarisation of cells, have recently been reported to be closely associated with several types of cancer and might thus serve as potential prognostic markers in the context of CRC. METHODS: We used a previously established meta-analysis of publicly available gene expression data to analyse the expression of different members of the Myosin V family, namely MYO5A, 5B, and 5C, in CRC. Using laser-microdissected material as well as tissue microarrays from paired human CRC samples, we validated both RNA and protein expression of Myosin Vb (MYO5B) and its known adapter proteins (RAB8A and RAB25) in an independent patient cohort. Finally, we assessed the prognostic value of both MYO5B and its adapter-coupled combinatorial gene expression signatures. RESULTS: The meta-analysis as well as an independent patient cohort study revealed a methylation-independent loss of MYO5B expression in CRC that matched disease progression. Although MYO5B mutations were identified in a small number of patients, these cannot be solely responsible for the common downregulation observed in CRC patients. Significantly, CRC patients with low MYO5B expression displayed shorter overall, disease-, and metastasis-free survival, a trend that was further reinforced when RAB8A expression was also taken into account. CONCLUSIONS: Our data identify MYO5B as a powerful prognostic biomarker in CRC, especially in early stages (stages I and II), which might help stratifying patients with stage II for adjuvant chemotherapy.

Speciation with gene flow in equids despite extensive chromosomal plasticity.

Horses, asses, and zebras belong to a single genus, Equus, which emerged 4.0-4.5 Mya. Although the equine fossil record represents a textbook example of evolution, the succession of events that gave rise to the diversity of species existing today remains unclear. Here we present six genomes from each living species of asses and zebras. This completes the set of genomes available for all extant species in the genus, which was hitherto represented only by the horse and the domestic donkey. In addition, we used a museum specimen to characterize the genome of the quagga zebra, which was driven to extinction in the early 1900s. We scan the genomes for lineage-specific adaptations and identify 48 genes that have evolved under positive selection and are involved in olfaction, immune response, development, locomotion, and behavior. Our extensive genome dataset reveals a highly dynamic demographic history with synchronous expansions and collapses on different continents during the last 400 ky after major climatic events. We show that the earliest speciation occurred with gene flow in Northern America, and that the ancestor of present-day asses and zebras dispersed into the Old World 2.1-3.4 Mya. Strikingly, we also find evidence for gene flow involving three contemporary equine species despite chromosomal numbers varying from 16 pairs to 31 pairs. These findings challenge the claim that the accumulation of chromosomal rearrangements drive complete reproductive isolation, and promote equids as a fundamental model for understanding the interplay between chromosomal structure, gene flow, and, ultimately, speciation.

Prehistoric genomes reveal the genetic foundation and cost of horse domestication.

The domestication of the horse ∼ 5.5 kya and the emergence of mounted riding, chariotry, and cavalry dramatically transformed human civilization. However, the genetics underlying horse domestication are difficult to reconstruct, given the near extinction of wild horses. We therefore sequenced two ancient horse genomes from Taymyr, Russia (at 7.4- and 24.3-fold coverage), both predating the earliest archeological evidence of domestication. We compared these genomes with genomes of domesticated horses and the wild Przewalski's horse and found genetic structure within Eurasia in the Late Pleistocene, with the ancient population contributing significantly to the genetic variation of domesticated breeds. We furthermore identified a conservative set of 125 potential domestication targets using four complementary scans for genes that have undergone positive selection. One group of genes is involved in muscular and limb development, articular junctions, and the cardiac system, and may represent physiological adaptations to human utilization. A second group consists of genes with cognitive functions, including social behavior, learning capabilities, fear response, and agreeableness, which may have been key for taming horses. We also found that domestication is associated with inbreeding and an excess of deleterious mutations. This genetic load is in line with the "cost of domestication" hypothesis also reported for rice, tomatoes, and dogs, and it is generally attributed to the relaxation of purifying selection resulting from the strong demographic bottlenecks accompanying domestication. Our work demonstrates the power of ancient genomes to reconstruct the complex genetic changes that transformed wild animals into their domesticated forms, and the population context in which this process took place.

Role of Calprotectin as a Modulator of the IL27-Mediated Proinflammatory Effect on Endothelial Cells.

An underlying endothelial dysfunction plays a fundamental role in the pathogenesis of cardiovascular events and is the central feature of atherosclerosis. The protein-based communication between leukocytes and inflamed endothelial cells leading to diapedesis has been largely investigated and several key players such as IL6, TNFα, or the damage associated molecular pattern molecule (DAMP) calprotectin are now well identified. However, regarding cytokine IL27, the controversial current knowledge about its inflammatory role and the involved regulatory elements requires clarification. Therefore, we examined the inflammatory impact of IL27 on primary endothelial cells and the potentially modulatory effect of calprotectin on both transcriptome and proteome levels. A qPCR-based screening demonstrated high IL27-mediated gene expression of IL7, IL15, CXCL10, and CXCL11. Calprotectin time-dependent downregulatory effects were observed on IL27-induced IL15 and CXCL10 gene expression. A mass spectrometry-based approach of IL27 ± calprotectin cell stimulation enlightened a calprotectin modulatory role in the expression of 28 proteins, mostly involved in the mechanism of leukocyte transmigration. Furthermore, we showed evidence for STAT1 involvement in this process. Our findings provide new evidence about the IL27-dependent proinflammatory signaling which may be under the control of calprotectin and highlight the need for further investigations on molecules which might have antiatherosclerotic functions.

True single-molecule DNA sequencing of a pleistocene horse bone.

Second-generation sequencing platforms have revolutionized the field of ancient DNA, opening access to complete genomes of past individuals and extinct species. However, these platforms are dependent on library construction and amplification steps that may result in sequences that do not reflect the original DNA template composition. This is particularly true for ancient DNA, where templates have undergone extensive damage post-mortem. Here, we report the results of the first "true single molecule sequencing" of ancient DNA. We generated 115.9 Mb and 76.9 Mb of DNA sequences from a permafrost-preserved Pleistocene horse bone using the Helicos HeliScope and Illumina GAIIx platforms, respectively. We find that the percentage of endogenous DNA sequences derived from the horse is higher among the Helicos data than Illumina data. This result indicates that the molecular biology tools used to generate sequencing libraries of ancient DNA molecules, as required for second-generation sequencing, introduce biases into the data that reduce the efficiency of the sequencing process and limit our ability to fully explore the molecular complexity of ancient DNA extracts. We demonstrate that simple modifications to the standard Helicos DNA template preparation protocol further increase the proportion of horse DNA for this sample by threefold. Comparison of Helicos-specific biases and sequence errors in modern DNA with those in ancient DNA also reveals extensive cytosine deamination damage at the 3' ends of ancient templates, indicating the presence of 3'-sequence overhangs. Our results suggest that paleogenomes could be sequenced in an unprecedented manner by combining current second- and third-generation sequencing approaches.

mapDamage2.0: fast approximate Bayesian estimates of ancient DNA damage parameters.

MOTIVATION: Ancient DNA (aDNA) molecules in fossilized bones and teeth, coprolites, sediments, mummified specimens and museum collections represent fantastic sources of information for evolutionary biologists, revealing the agents of past epidemics and the dynamics of past populations. However, the analysis of aDNA generally faces two major issues. Firstly, sequences consist of a mixture of endogenous and various exogenous backgrounds, mostly microbial. Secondly, high nucleotide misincorporation rates can be observed as a result of severe post-mortem DNA damage. Such misincorporation patterns are instrumental to authenticate ancient sequences versus modern contaminants. We recently developed the user-friendly mapDamage package that identifies such patterns from next-generation sequencing (NGS) sequence datasets. The absence of formal statistical modeling of the DNA damage process, however, precluded rigorous quantitative comparisons across samples. RESULTS: Here, we describe mapDamage 2.0 that extends the original features of mapDamage by incorporating a statistical model of DNA damage. Assuming that damage events depend only on sequencing position and post-mortem deamination, our Bayesian statistical framework provides estimates of four key features of aDNA molecules: the average length of overhangs (λ), nick frequency (ν) and cytosine deamination rates in both double-stranded regions ( ) and overhangs ( ). Our model enables rescaling base quality scores according to their probability of being damaged. mapDamage 2.0 handles NGS datasets with ease and is compatible with a wide range of DNA library protocols. AVAILABILITY: mapDamage 2.0 is available at ginolhac.github.io/mapDamage/ as a Python package and documentation is maintained at the Centre for GeoGenetics Web site (geogenetics.ku.dk/publications/mapdamage2.0/). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

2-Hydroxyglutarate modulates histone methylation at specific loci and alters gene expression via Rph1 inhibition.

2-Hydroxyglutarate (2-HG) is an oncometabolite that accumulates in certain cancers. Gain-of-function mutations in isocitrate dehydrogenase lead to 2-HG accumulation at the expense of alpha-ketoglutarate. Elevated 2-HG levels inhibit histone and DNA demethylases, causing chromatin structure and gene regulation changes with tumorigenic consequences. We investigated the effects of elevated 2-HG levels in Saccharomyces cerevisiae, a yeast devoid of DNA methylation and heterochromatin-associated histone methylation. Our results demonstrate genetic background-dependent gene expression changes and altered H3K4 and H3K36 methylation at specific loci. Analysis of histone demethylase deletion strains indicated that 2-HG inhibits Rph1 sufficiently to induce extensive gene expression changes. Rph1 is the yeast homolog of human KDM4 demethylases and, among the yeast histone demethylases, was the most sensitive to the inhibitory effect of 2-HG in vitro. Interestingly, Rph1 deficiency favors gene repression and leads to further down-regulation of already silenced genes marked by low H3K4 and H3K36 trimethylation, but abundant in H3K36 dimethylation. Our results provide novel insights into the genome-wide effects of 2-HG and highlight Rph1 as its preferential demethylase target.

Proteomic analysis of a pleistocene mammoth femur reveals more than one hundred ancient bone proteins.

We used high-sensitivity, high-resolution tandem mass spectrometry to shotgun sequence ancient protein remains extracted from a 43 000 year old woolly mammoth ( Mammuthus primigenius ) bone preserved in the Siberian permafrost. For the first time, 126 unique protein accessions, mostly low-abundance extracellular matrix and plasma proteins, were confidently identified by solid molecular evidence. Among the best characterized was the carrier protein serum albumin, presenting two single amino acid substitutions compared to extant African ( Loxodonta africana ) and Indian ( Elephas maximus ) elephants. Strong evidence was observed of amino acid modifications due to post-mortem hydrolytic and oxidative damage. A consistent subset of this permafrost bone proteome was also identified in more recent Columbian mammoth ( Mammuthus columbi ) samples from temperate latitudes, extending the potential of the approach described beyond subpolar environments. Mass spectrometry-based ancient protein sequencing offers new perspectives for future molecular phylogenetic inference and physiological studies on samples not amenable to ancient DNA investigation. This approach therefore represents a further step into the ongoing integration of different high-throughput technologies for identification of ancient biomolecules, unleashing the field of paleoproteomics.

Improving ancient DNA read mapping against modern reference genomes.

BACKGROUND: Next-Generation Sequencing has revolutionized our approach to ancient DNA (aDNA) research, by providing complete genomic sequences of ancient individuals and extinct species. However, the recovery of genetic material from long-dead organisms is still complicated by a number of issues, including post-mortem DNA damage and high levels of environmental contamination. Together with error profiles specific to the type of sequencing platforms used, these specificities could limit our ability to map sequencing reads against modern reference genomes and therefore limit our ability to identify endogenous ancient reads, reducing the efficiency of shotgun sequencing aDNA. RESULTS: In this study, we compare different computational methods for improving the accuracy and sensitivity of aDNA sequence identification, based on shotgun sequencing reads recovered from Pleistocene horse extracts using Illumina GAIIx and Helicos Heliscope platforms. We show that the performance of the Burrows Wheeler Aligner (BWA), that has been developed for mapping of undamaged sequencing reads using platforms with low rates of indel-types of sequencing errors, can be employed at acceptable run-times by modifying default parameters in a platform-specific manner. We also examine if trimming likely damaged positions at read ends can increase the recovery of genuine aDNA fragments and if accurate identification of human contamination can be achieved using a strategy previously suggested based on best hit filtering. We show that combining our different mapping and filtering approaches can increase the number of high-quality endogenous hits recovered by up to 33%. CONCLUSIONS: We have shown that Illumina and Helicos sequences recovered from aDNA extracts could not be aligned to modern reference genomes with the same efficiency unless mapping parameters are optimized for the specific types of errors generated by these platforms and by post-mortem DNA damage. Our findings have important implications for future aDNA research, as we define mapping guidelines that improve our ability to identify genuine aDNA sequences, which in turn could improve the genotyping accuracy of ancient specimens. Our framework provides a significant improvement to the standard procedures used for characterizing ancient genomes, which is challenged by contamination and often low amounts of DNA material.

Improving the performance of true single molecule sequencing for ancient DNA.

BACKGROUND: Second-generation sequencing technologies have revolutionized our ability to recover genetic information from the past, allowing the characterization of the first complete genomes from past individuals and extinct species. Recently, third generation Helicos sequencing platforms, which perform true Single-Molecule DNA Sequencing (tSMS), have shown great potential for sequencing DNA molecules from Pleistocene fossils. Here, we aim at improving even further the performance of tSMS for ancient DNA by testing two novel tSMS template preparation methods for Pleistocene bone fossils, namely oligonucleotide spiking and treatment with DNA phosphatase. RESULTS: We found that a significantly larger fraction of the horse genome could be covered following oligonucleotide spiking however not reproducibly and at the cost of extra post-sequencing filtering procedures and skewed %GC content. In contrast, we showed that treating ancient DNA extracts with DNA phosphatase improved the amount of endogenous sequence information recovered per sequencing channel by up to 3.3-fold, while still providing molecular signatures of endogenous ancient DNA damage, including cytosine deamination and fragmentation by depurination. Additionally, we confirmed the existence of molecular preservation niches in large bone crystals from which DNA could be preferentially extracted. CONCLUSIONS: We propose DNA phosphatase treatment as a mechanism to increase sequence coverage of ancient genomes when using Helicos tSMS as a sequencing platform. Together with mild denaturation temperatures that favor access to endogenous ancient templates over modern DNA contaminants, this simple preparation procedure can improve overall Helicos tSMS performance when damaged DNA templates are targeted.

mapDamage: testing for damage patterns in ancient DNA sequences.

SUMMARY: Ancient DNA extracts consist of a mixture of contaminant DNA molecules, most often originating from environmental microbes, and endogenous fragments exhibiting substantial levels of DNA damage. The latter introduce specific nucleotide misincorporations and DNA fragmentation signatures in sequencing reads that could be advantageously used to argue for sequence validity. mapDamage is a Perl script that computes nucleotide misincorporation and fragmentation patterns using next-generation sequencing reads mapped against a reference genome. The Perl script outputs are further automatically processed in embedded R script in order to detect typical patterns of genuine ancient DNA sequences. AVAILABILITY AND IMPLEMENTATION: The Perl script mapDamage is freely available with documentation and example files at http://geogenetics.ku.dk/all_literature/mapdamage/. The script requires prior installation of the SAMtools suite and R environment and has been validated on both GNU/Linux and MacOSX operating systems.

Focal adhesion kinase plays a dual role in TRAIL resistance and metastatic outgrowth of malignant melanoma.

Despite remarkable advances in therapeutic interventions, malignant melanoma (MM) remains a life-threating disease. Following high initial response rates to targeted kinase-inhibition metastases quickly acquire resistance and present with enhanced tumor progression and invasion, demanding alternative treatment options. We show 2nd generation hexameric TRAIL-receptor-agonist IZI1551 (IZI) to effectively induce apoptosis in MM cells irrespective of the intrinsic BRAF/NRAS mutation status. Conditioning to the EC50 dose of IZI converted the phenotype of IZI-sensitive parental MM cells into a fast proliferating and invasive, IZI-resistant metastasis. Mechanistically, we identified focal adhesion kinase (FAK) to play a dual role in phenotype-switching. In the cytosol, activated FAK triggers survival pathways in a PI3K- and MAPK-dependent manner. In the nucleus, the FERM domain of FAK prevents activation of wtp53, as being expressed in the majority of MM, and consequently intrinsic apoptosis. Caspase-8-mediated cleavage of FAK as well as FAK knockdown, and pharmacological inhibition, respectively, reverted the metastatic phenotype-switch and restored IZI responsiveness. FAK inhibition also re-sensitized MM cells isolated from patient metastasis that had relapsed from targeted kinase inhibition to cell death, irrespective of the intrinsic BRAF/NRAS mutation status. Hence, FAK-inhibition alone or in combination with 2nd generation TRAIL-receptor agonists may be recommended for treatment of initially resistant and relapsed MM, respectively.