Search

RNA-DNA hybrids are widespread epigenetic features of genomes that provide a growing range of activities in transcription, chromatin and DNA replication and repair. Understanding of these diverse functions has been advanced by characterising the proteins that interact with the hybrids, with all such studies revealing hundreds of potential interactors. However, all interaction analyses to date have focused on mammalian cells, and so it is unclear if a similar spectrum of RNA-DNA hybrid interactors is found in other eukaryotes, thus limiting our understanding of the conserved and lineage-specific activities linked to these genetic structures. The African trypanosome is a compelling organism in which to address these questions. As a divergent single-cell eukaryotic parasite of the Discoba grouping,Trypanosoma bruceidisplays substantial divergence in several aspects of core biology from its mammalian host and, unusually for a protist, has well-developed tools for molecular genetic analysis. For these reasons, we used DNA-RNA hybrid immunoprecipitation coupled with mass spectrometry to reveal 602 putative interactors inT. bruceimammal- or insect vector-infective stage cells. We show that the approach selects for a subset of the parasite proteome and reveals a range of predicted RNA-DNA hybrid associated activities, some overlapping with similar studies in mammals. We demonstrate that loss of three factors, two putative helicases and a RAD51 paralogue, impact onT. bruceinuclear RNA-DNA hybrid and DNA damage levels. Moreover, loss of each affects the operation of the crucial parasite immune survival mechanism of antigenic variation. Thus, our work reveals the broad range of activities contributed by RNA-DNA hybrids toT. bruceibiology, including new functions in host immune evasion as well as many conserved with mammals, and so likely fundamental to eukaryotic genome function.

RNA-DNA hybrids are epigenetic features of all genomes that intersect with many processes, including transcription, telomere homeostasis and centromere function. Increasing evidence suggests RNA-DNA hybrids can provide two conflicting roles in the maintenance and transmission of genomes: they can be the triggers of DNA damage, leading to genome change, or can aid the DNA repair processes needed to respond to DNA lesions. Evasion of host immunity by African trypanosomes, such asTrypanosoma brucei, relies on targeted recombination of silent Variant Surface Glycoprotein (VSG) genes into a specialised telomeric locus that directs transcription of just oneVSGfrom thousands. How suchVSGrecombination is targeted and initiated is unclear. Here, we show that a key enzyme ofT. bruceihomologous recombination, RAD51, interacts with RNA-DNA hybrids. In addition, we show that RNA-DNA hybrids display a genome- wide co-localisation with DNA breaks, and that this relationship is impaired by mutation of RAD51. Finally, we show that RAD51 acts to repair highly abundant, localised DNA breaks at the single transcribedVSG, and that mutation of RAD51 alters RNA-DNA hybrid abundance both around the transcribedVSGand across the silentVSGarchive. This work reveals a widespread, generalised role for RNA-DNA hybrids in directing RAD51 activity during recombination and uncovers a specialised application of this interplay during targeted DNA break repair needed for the criticalT. bruceiimmune evasion reaction of antigenic variation.

Malaria remains a global health problem causing more than 400,000 deaths annually. Plasmodium parasites, the causative agents of malaria, replicate asexually in red blood cells (RBCs) of their vertebrate host, while a subset differentiates into sexual stages (gametocytes) for mosquito transmission. Parasite replication and gametocyte maturation in the erythropoietic niches of the bone marrow and spleen contribute to pathogenesis and drive transmission, but the mechanisms underlying this organ enrichment remain unknown. Here, we performed a comprehensive analysis of rodent P. berghei infection by flow cytometry and single-cell RNA sequencing. We identified CD71 as a host receptor for reticulocyte invasion and found that parasites metabolically adapt to the host cell environment. Transcriptional analysis and functional assays further revealed a nutrient-dependent tropism for gametocyte formation in reticulocytes. Together, we provide a thorough characterization of host-parasite interactions in erythropoietic niches and define host cell maturation state as the key driver of parasite adaptation.

Plasmodium vivax is the major cause of human malaria in the Americas. How P. vivax infection can lead to poor clinical outcomes, despite low peripheral parasitaemia remains a matter of intense debate. Estimation of total P. vivax biomass based on circulating markers indicates existence of a predominant parasite population outside of circulation. In this study we investigate associations between both peripheral and total parasite biomass and host response in vivax malaria. We analysed parasite and host signatures in a cohort of uncomplicated vivax malaria patients from Manaus, Brazil, combining clinical and parasite parameters, multiplexed analysis of host responses and ex vivo assays. Patterns of clinical features, parasite burden and host signatures measured in plasma across the patient cohort were highly heterogenous. Further data deconvolution revealed two patient clusters, here termed Vivaxlow and Vivaxhigh. These patient subgroups were defined based on differences in total parasite biomass but not peripheral parasitaemia. Overall Vivaxlow patients clustered with healthy donors and Vivaxhigh patients showed more profound alterations in haematological parameters, endothelial cell (EC) activation and glycocalyx breakdown and levels of cytokines regulating different haematopoiesis pathways compared to Vivaxlow. Vivaxhigh patients presented more severe thrombocytopenia and lymphopenia, along with enrichment of neutrophils in the peripheral blood and increased neutrophil-to-lymphocyte ratio (NLCR). When patients’ signatures were combined, high association of total parasite biomass with a subset of markers of EC activation, thrombocytopenia and lymphopenia severity was observed. Finally, machine learning models defined a combination of host parameters measured in the circulation that could predict the extent of parasite infection outside of circulation. Altogether, our data show that total parasite biomass is a better predictor of perturbations in host homeostasis in P. vivax patients than peripheral parasitaemia. This supports the emerging paradigm of a P. vivax tissue reservoir, in particular in the hematopoietic niche of bone marrow and spleen.

Malaria remains a global health problem with over 400,000 deaths annually1. Plasmodium parasites, the causative agents of malaria, replicate asexually in red blood cells (RBCs) of their vertebrate host, while a subset differentiates into sexual stages (gametocytes) for mosquito transmission. Parasite replication and gametocyte maturation in the erythropoietic niches of the bone marrow and spleen contribute to pathogenesis and drive transmission2, but the mechanisms underlying this organ enrichment remain unknown. We performed a comprehensive single cell analysis of rodent P. berghei in spleen, bone marrow and blood to define parasite phenotypes specific to those niches. Single cell RNA-seq analysis of host and parasite cells reveals an interferon-driven host response to infection as well as transcriptional adaptations of Plasmodium to RBC maturation status. We show that P. berghei exhibits a bimodal invasion pattern into either normocytes or early reticulocytes and, using functional assays, identify CD71 as a host receptor for reticulocyte invasion. Importantly, we observe an increased rate of gametocyte formation in reticulocytes that is nutrient-dependent and triggered post invasion (i.e., same cycle sexual commitment). Our data provides a thorough characterisation of host-parasite interactions in erythropoietic niches and defines host cell maturation state as the key driver of parasite adaptation.

Background & Aims Continuing recalcitrance to therapy cements pancreatic cancer (PC) as the most lethal malignancy, which is set to become the second leading cause of cancer death in our society. The study aim was to investigate the association between DNA damage response (DDR), replication stress, and novel therapeutic response in PC to develop a biomarker-driven therapeutic strategy targeting DDR and replication stress in PC. Methods We interrogated the transcriptome, genome, proteome, and functional characteristics of 61 novel PC patient–derived cell lines to define novel therapeutic strategies targeting DDR and replication stress. Validation was done in patient-derived xenografts and human PC organoids. Results Patient-derived cell lines faithfully recapitulate the epithelial component of pancreatic tumors, including previously described molecular subtypes. Biomarkers of DDR deficiency, including a novel signature of homologous recombination deficiency, cosegregates with response to platinum (P < .001) and PARP inhibitor therapy (P < .001) in vitro and in vivo. We generated a novel signature of replication stress that predicts response to ATR (P < .018) and WEE1 inhibitor (P < .029) treatment in both cell lines and human PC organoids. Replication stress was enriched in the squamous subtype of PC (P < .001) but was not associated with DDR deficiency. Conclusions Replication stress and DDR deficiency are independent of each other, creating opportunities for therapy in DDR-proficient PC and after platinum therapy., Graphical abstract

Homologous recombination (HR) has an intimate relationship with genome replication, both during repair of DNA lesions that might prevent DNA synthesis and in tackling stalls to the replication fork. Recent studies led us to ask if HR might have a more central role in replicating the genome of Leishmania, a eukaryotic parasite. Conflicting evidence has emerged regarding whether or not HR genes are essential, and genome-wide mapping has provided evidence for an unorthodox organisation of DNA replication initiation sites, termed origins. To answer this question, we have employed a combined CRISPR/Cas9 and DiCre approach to rapidly generate and assess the effect of conditional ablation of RAD51 and three RAD51-related proteins in Leishmania major. Using this approach, we demonstrate that loss of any of these HR factors is not immediately lethal but in each case growth slows with time and leads to DNA damage and accumulation of cells with aberrant DNA content. Despite these similarities, we show that only loss of RAD51 or RAD51-3 impairs DNA synthesis and causes elevated levels of genome-wide mutation. Furthermore, we show that these two HR factors act in distinct ways, since ablation of RAD51, but not RAD51-3, has a profound effect on DNA replication, causing loss of initiation at the major origins and increased DNA synthesis at subtelomeres. Our work clarifies questions regarding the importance of HR to survival of Leishmania and reveals an unanticipated, central role for RAD51 in the programme of genome replication in a microbial eukaryote., Author summary Homologous recombination plays a key role in genome maintenance during cell division, but loss of factors directing the reaction has not been described as being lethal in any microbe. Here, we have used a genetic strategy to selectively induce loss, singly and doubly, of five genes in Leishmania that act in homologous recombination, revealing two things. First, loss of any gene related to RAD51, which catalyses homologous recombination, is not immediately lethal, but leads to increasing growth impairment and genome damage accumulation. Second, loss of RAD51 causes a pronounced change in the programme of Leishmania genome replication. Thus, we show that homologous recombination in Leishmania can be essential, in part due to an unanticipated role in genome transmission.

Control of DNA methylation level is critical for gene regulation, and the factors that govern hypomethylation at CpG islands (CGIs) are still being uncovered. Here, we provide evidence that G-quadruplex (G4) DNA secondary structures are genomic features that influence methylation at CGIs. We show that the presence of G4 structure is tightly associated with CGI hypomethylation in the human genome. Surprisingly, we find that these G4 sites are enriched for DNA methyltransferase 1 (DNMT1) occupancy, which is consistent with our biophysical observations that DNMT1 exhibits higher binding affinity for G4s as compared to duplex, hemi-methylated, or single-stranded DNA. The biochemical assays also show that the G4 structure itself, rather than sequence, inhibits DNMT1 enzymatic activity. Based on these data, we propose that G4 formation sequesters DNMT1 thereby protecting certain CGIs from methylation and inhibiting local methylation.

Once every cell cycle, DNA replication takes place to allow cells to duplicate their genome and segregate the two resulting copies into offspring cells. In eukaryotes, the number of DNA replication initiation loci, termed origins, is proportional to chromosome size. However, previous studies have suggested that inLeishmania, a group of single-celled eukaryotic parasites, DNA replication starts from just a single origin per chromosome, which is predicted to be insufficient to secure complete genome duplication within S phase. Here, we show that the paucity of origins activated in early S phase is balanced by DNA synthesis activity outside S phase. Simultaneous recruitment of acetylated histone H3 (AcH3), modified base J and the kinetochore factor KKT1 is exclusively found at the origins used in early S phase, while subtelomeric DNA replication can only be linked to AcH3 and displays persistent activity through the cell cycle, including in G2/M and G1 phases. We also show that subtelomeric DNA replication, unlike replication from the previously mapped origins, is sensitive to hydroxyurea and dependent on subunits of the 9-1-1 complex. Our work indicates thatLeishmaniagenome transmission relies on an unconventional DNA replication programme, which may have implications for genome stability in this important parasite.

Next-generation sequencing is enabling molecularly guided therapy for many cancer types, yet failure rates remain relatively high in pancreatic cancer (PC). The aim of this study is to investigate the feasibility of genomic profiling using endoscopic ultrasound (EUS) biopsy samples to facilitate personalised therapy for PC. Ninty-five patients underwent additional research biopsies at the time of diagnostic EUS. Diagnostic formalin-fixed (FFPE) and fresh frozen EUS samples underwent DNA extraction, quantification and targeted gene sequencing. Whole genome (WGS) and RNA sequencing was performed as proof of concept. Only 2 patients (2%) with a diagnosis of PC had insufficient material for targeted sequencing in both FFPE and frozen specimens. Targeted panel sequencing (n=54) revealed mutations in PC genes (KRAS, GNAS, TP53, CDKN2A, SMAD4) in patients with histological evidence of PC, including potentially actionable mutations (BRCA1, BRCA2, ATM, BRAF). WGS (n=5) of EUS samples revealed mutational signatures that are potential biomarkers of therapeutic responsiveness. RNA sequencing (n=35) segregated patients into clinically relevant molecular subtypes based on transcriptome. Integrated multi-omic analysis of PC using standard EUS guided biopsies offers clinical utility to guide personalized therapy and study the molecular pathology in all patients with PC.

Nucleosomes are the basic unit of chromatin that help the packaging of genetic material while controlling access to the genetic information. The underlying DNA sequence, together with transcription-associated proteins and chromatin remodelling complexes, are important factors that influence the organization of nucleosomes. Here, we show that the naturally occurring DNA modification, 5-formylcytosine (5fC) is linked to tissue-specific nucleosome organization. Our study reveals that 5fC is associated with increased nucleosome occupancy in vitro and in vivo. We demonstrate that 5fC-associated nucleosomes at enhancers in the mammalian hindbrain and heart are linked to elevated gene expression. Our study also reveals the formation of a reversible-covalent Schiff base linkage between lysines of histone proteins and 5fC within nucleosomes in a cellular environment. We define their specific genomic loci in mouse embryonic stem cells and look into the biological consequences of these DNA-histone Schiff base sites. Collectively, our findings show that 5fC is a determinant of nucleosome organization and plays a role in establishing distinct regulatory regions that control transcription.

Nucleosomes are the basic unit of chromatin that ensure genome integrity and control access to the genetic information. The organization of nucleosomes is influenced by the underlying DNA sequence itself, transcription factors or other transcriptional machinery associated proteins and chromatin remodeling complexes (1–4). Herein, we show that the naturally occurring DNA modification, 5-formylcytosine (5fC) contributes to the positioning of nucleosomes. We show that the ability of 5fC to position nucleosomesin vitrois associated with the formation of covalent interactions between histone residues and 5fC in the form of Schiff bases. We demonstrate that similar interactions can occur in a cellular environment and define their specific genomic loci in mouse embryonic stem cells. Collectively, our findings identify 5fC as a determinant of nucleosomal organization in which 5fC plays a role in establishing distinct regulatory regions that are linked to gene expression Our study provides a previously unknown molecular mechanism, involving the formation of reversible-covalent bonds between chromatin and DNA that supports a molecular linkage between DNA sequence, DNA base modification and chromatin structure.

The cytosine modifications 5-hydroxymethylcytosine (5hmC) and 5-formylcytosine (5fC) were recently found to exist in the genomic DNA of a wide range of mammalian cell types. It is important to now understand their role in normal biological function and disease. Here we introduce reduced bisulfite sequencing (redBS-Seq), a method to quantitatively decode 5fC in DNA at singlebase resolution, based on a selective chemical reduction of 5fC to 5hmC followed by bisulfite treatment. Following extensive validation on synthetic and genomic DNA, we combined redBS-Seq and oxidative bisulfite sequencing (oxBS-Seq) to generate the first combined genomic map of 5-methylcytosine, 5hmC and 5fC in mouse embryonic stem cells. Our experiments revealed that 5fC is present at relatively high levels in certain genomic locations, in comparison to 5hmC and 5mC. The combination of these chemical methods can quantify and precisely map these three cytosine derivatives in the genome and will help provide insights into their function.

Summary: Current genome browsers are designed to work via graphical user interfaces (GUIs), which, however intuitive, are not amenable to operate within console terminals and therefore are difficult to streamline or integrate in scripts. To circumvent these limitations, ASCIIGenome runs exclusively via command line interface to display genomic data directly in a terminal window. By following the same philosophy of UNIX tools, ASCIIGenome aims to be easily integrated with the command line, including batch processing of data, and therefore enables an effective exploration of the data. Implementation: ASCIIGenome is written in Java. Consequently, it is a cross-platform tool and requires minimal or no installation. Some of the common genomic data types are supported and data access on remote ftp servers is possible. Speed and memory footprint are comparable to or better than those of common genome browsers. Availability and Implementation: Software and source code (MIT License) are available at https://github.com/dariober/ASCIIGenome with detailed documentation at http://asciigenome.readthedocs.io. Contact: Dario.beraldi@cruk.cam.ac.uk Supplementary information: Supplementary data are available at Bioinformatics online.

Background Genome-wide methylation of cytosine can be modulated in the presence of TET and thymine DNA glycosylase (TDG) enzymes. TET is able to oxidise 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC). TDG can excise the oxidative products 5fC and 5caC, initiating base excision repair. These modified bases are stable and detectable in the genome, suggesting that they could have epigenetic functions in their own right. However, functional investigation of the genome-wide distribution of 5fC has been restricted to cell culture-based systems, while its in vivo profile remains unknown. Results Here, we describe the first analysis of the in vivo genome-wide profile of 5fC across a range of tissues from both wild-type and Tdg-deficient E11.5 mouse embryos. Changes in the formylation profile of cytosine upon depletion of TDG suggest TET/TDG-mediated active demethylation occurs preferentially at intron-exon boundaries and reveals a major role for TDG in shaping 5fC distribution at CpG islands. Moreover, we find that active enhancer regions specifically exhibit high levels of 5fC, resulting in characteristic tissue-diagnostic patterns, which suggest a role in embryonic development. Conclusions The tissue-specific distribution of 5fC can be regulated by the collective contribution of TET-mediated oxidation and excision by TDG. The in vivo profile of 5fC during embryonic development resembles that of embryonic stem cells, sharing key features including enrichment of 5fC in enhancer and intragenic regions. Additionally, by investigating mouse embryo 5fC profiles in a tissue-specific manner, we identify targeted enrichment at active enhancers involved in tissue development. Electronic supplementary material The online version of this article (doi:10.1186/s13059-016-1001-5) contains supplementary material, which is available to authorized users.

Identifying the genes underlying phenotypic variation in natural populations can provide novel insight into the evolutionary process. The candidate gene approach has been applied to studies of a number of traits in various species, in an attempt to elucidate their genetic basis. Here, we test the application of the candidate gene approach to identify the loci involved in variation in gastrointestinal parasite burden, a complex trait likely to be controlled by many loci, in a wild population of Soay sheep. A comprehensive literature review, Gene Ontology databases, and comparative genomics resources between cattle and sheep were used to generate a list of candidate genes. In a pilot study, these candidates, along with 50 random genes, were then sequenced in two pools of Soay sheep; one with low gastrointestinal nematode burden and the other high, using a NimbleGen sequence capture experiment. Further candidates were identified from single nucleotide polymorphisms (SNPs) that were highly differentiated between high- and low-resistance sheep breeds. A panel of 192 candidate and control SNPs were then typed in 960 individual Soay sheep to examine whether they individually explained variation in parasite burden, as measured as faecal egg count, as well as two immune measures (Teladorsagia circumcincta-specific antibodies and antinuclear antibodies). The cumulative effect of the candidate and control SNPs were estimated by fitting genetic relationship matrices (GRMs) as random effects in animal models of the three traits. No more significant SNPs were identified in the pilot sequencing experiment and association study than expected by chance. Furthermore, no significant difference was found between the proportions of candidate or control SNPs that were found to be significantly associated with parasite burden/immune measures. No significant effect of the candidate or control gene GRMs was found. There is thus little support for the candidate gene approach to the identification of loci explaining variation in parasitological and immunological traits in this population. However, a number of SNPs explained significant variation in multiple traits and significant correlations were found between the proportions of variance explained by individual SNPs across multiple traits. The significant SNPs identified in this study may still, therefore, merit further investigation.

Understanding the genetic architecture of phenotypic variation in natural populations is a fundamental goal of evolutionary genetics. Wild Soay sheep (Ovis aries) have an inherited polymorphism for horn morphology in both sexes, controlled by a single autosomal locus, Horns. The majority of males have large normal horns, but a small number have vestigial, deformed horns, known as scurs; females have either normal horns, scurs or no horns (polled). Given that scurred males and polled females have reduced fitness within each sex, it is counterintuitive that the polymorphism persists within the population. Therefore, identifying the genetic basis of horn type will provide a vital foundation for understanding why the different morphs are maintained in the face of natural selection. We conducted a genome-wide association study using ∼36000 single nucleotide polymorphisms (SNPs) and determined the main candidate for Horns as RXFP2, an autosomal gene with a known involvement in determining primary sex characters in humans and mice. Evidence from additional SNPs in and around RXFP2 supports a new model of horn-type inheritance in Soay sheep, and for the first time, sheep with the same horn phenotype but different underlying genotypes can be identified. In addition, RXFP2 was shown to be an additive quantitative trait locus (QTL) for horn size in normal-horned males, accounting for up to 76% of additive genetic variation in this trait. This finding contrasts markedly from genome-wide association studies of quantitative traits in humans and some model species, where it is often observed that mapped loci only explain a modest proportion of the overall genetic variation.

Digital gene expression (DGE) analysis offers a route to gene discovery which by-passes the need to develop bespoke arrays for nonmodel species, and is therefore a potentially valuable tool for molecular ecologists. Scottish blackface sheep, which vary in resistance to the common abomasal parasitic nematode Teladorsagia circumcincta, were trickle-infected with L3 larvae over 3 months to mimic the natural progression of infection. DGE was performed on abomasal lymph node tissue after the resolution of infection in resistant animals. Susceptible (low resistance) animals showed a large number of differentially expressed genes associated with inflammation and cell activation, but generally few differentially regulated genes in either the susceptible or the resistant group were directly involved in the adaptive immune function. Our results are consistent with the hypothesis that both resistance and susceptibility are active responses to infection and that susceptibility is associated with dysfunction in T cell differentiation and regulation.

Wing length is a key character for essential behaviours related to bird flight such as migration and foraging. In the present study, we initiate the search for the genes underlying wing length in birds by studying a long-distance migrant, the great reed warbler (Acrocephalus arundinaceus). In this species wing length is an evolutionary interesting trait with pronounced latitudinal gradient and sex-specific selection regimes in local populations. We performed a quantitative trait locus (QTL) scan for wing length in great reed warblers using phenotypic, genotypic, pedigree and linkage map data from our long-term study population in Sweden. We applied the linkage analysis mapping method implemented in GridQTL (a new web-based software) and detected a genome-wide significant QTL for wing length on chromosome 2, to our knowledge, the first detected QTL in wild birds. The QTL extended over 25 cM and accounted for a substantial part (37%) of the phenotypic variance of the trait. A genome scan for tarsus length (a body-size-related trait) did not show any signal, implying that the wing-length QTL on chromosome 2 was not associated with body size. Our results provide a first important step into understanding the genetic architecture of avian wing length, and give opportunities to study the evolutionary dynamics of wing length at the locus level.

The evolution of male weaponry in animals is driven by sexual selection, which is predicted to reduce the genetic variability underlying such traits. Soay sheep have an inherited polymorphism for horn type in both sexes, with males presenting with either large, normal horns or small, deformed horns (scurs). In addition, there is additive genetic variation in horn length among males with normal horns. Given that scurred males cannot win conflicts with normal-horned males, it is unusual that genes conferring scurs should persist in the population. Identifying the genetic basis of these traits should help us in understanding their evolution. We developed microsatellite markers in a targeted region of the Soay sheep genome and refined the location of the Horns locus (Ho) to a approximately 7.4 cM interval on chromosome 10 (LOD=8.78). We then located quantitative trait loci spanning a 34 cM interval with a peak centred close to Ho, which explained the majority of the genetic variation for horn length and base circumference in normal-horned males (LOD=2.51 and LOD=1.04, respectively). Therefore, the genetic variation in both horn type and horn length is attributable to the same chromosomal region. Understanding the maintenance of horn type and length variation will require an investigation of selection on genotypes that (co)determine both traits.

One of the biggest challenges facing evolutionary biologists is to identify and understand loci that explain fitness variation in natural populations. This review describes how genetic (linkage) mapping with single nucleotide polymorphism (SNP) markers can lead to great progress in this area. Strategies for SNP discovery and SNP genotyping are described and an overview of how to model SNP genotype information in mapping studies is presented. Finally, the opportunity afforded by new generation sequencing and typing technologies to map fitness genes by genome-wide association studies is discussed.

We searched for quantitative trait loci (QTL) underlying fitness-related traits in a free-living pedigree of 588 Soay sheep in which a genetic map using 251 markers with an average spacing of 15 cM had been established previously. Traits examined included birth date and weight, considered both as maternal and offspring traits, foreleg length, hindleg length, and body weight measured on animals in August and jaw length and metacarpal length measured on cleaned skeletal material. In some cases the data were split to consider different age classes separately, yielding a total of 15 traits studied. Genetic and environmental components of phenotypic variance were estimated for each trait and, for those traits showing nonzero heritability (N= 12), a QTL search was conducted by comparing a polygenic model with a model including a putative QTL. Support for a QTL at genome-wide significance was found on chromosome 11 for jaw length; suggestive QTL were found on chromosomes 2 and 5 (for birth date as a trait of the lamb), 8 (birth weight as a trait of the lamb), and 15 (adult hindleg length). We discuss the prospects for refining estimates of QTL position and effect size in the study population, and for QTL searches in free-living pedigrees in general.

Background The FANTOM5 consortium used Cap Analysis of Gene Expression (CAGE) tag sequencing to produce a comprehensive atlas of promoters and enhancers within the human and mouse genomes. We reasoned that the mapping of these regulatory elements to the pig genome could provide useful annotation and evidence to support assignment of orthology. Results For human transcription start sites (TSS) associated with annotated human-mouse orthologs, 17% mapped to the pig genome but not to the mouse, 10% mapped only to the mouse, and 55% mapped to both pig and mouse. Around 17% did not map to either species. The mapping percentages were lower where there was not clear orthology relationship, but in every case, mapping to pig was greater than to mouse, and the degree of homology was also greater. Combined mapping of mouse and human CAGE-defined promoters identified at least one putative conserved TSS for >16,000 protein-coding genes. About 54% of the predicted locations of regulatory elements in the pig genome were supported by CAGE and/or RNA-Seq analysis from pig macrophages. Conclusions Comparative mapping of promoters and enhancers from humans and mice can provide useful preliminary annotation of other animal genomes. The data also confirm extensive gain and loss of regulatory elements between species, and the likelihood that pigs provide a better model than mice for human gene regulation and function. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-2111-2) contains supplementary material, which is available to authorized users.

The modified base 5-formylcytosine (5fC) was recently identified in mammalian DNA and might be considered as the “seventh” base of the genome. This nucleotide has been implicated in active demethylation mediated by the base excision repair enzyme thymine DNA glycosylase (TDG). Genomics and proteomics studies have suggested a further role for 5fC in transcription regulation through chromatin remodeling. Herein we propose how 5fC might signal these processes through its effect on DNA conformation. Biophysical and structural analysis revealed that 5fC alters the structure of the DNA double helix leading to a conformation unique amongst known DNA structures including those comprising other cytosine modifications. The 1.4 Å resolution X-ray crystal structure of a DNA dodecamer comprising three 5fCpG sites shown how 5fC changes the geometry of the grooves and base pairs associated with the modified base, which lead to helical under-winding.

Identifying the genes that underlie phenotypic variation in natural populations is a central objective of evolutionary genetics. Here, we report the identification of the gene and causal mutation underlying coat colour variation in a free-living population of Soay sheep (Ovis aries). We targeted tyrosinase-related protein 1 (TYRP1), a positional candidate gene based on previous work that mapped theCoat colourlocus to an approximately 15 cM window on chromosome 2. We identified a non-synonymous substitution in exon IV that was perfectly associated with coat colour. This polymorphism is predicted to cause the loss of a cysteine residue that is highly evolutionarily conserved and likely to be of functional significance. We eliminated the possibility that this association is due to the presence of strong linkage disequilibrium with an unknown regulatory mutation by demonstrating that there is no difference in relativeTYRP1expression between colour morphs. Analysis of this putative causal mutation in a complex pedigree of more than 500 sheep revealed almost perfect co-segregation with coat colour (Χ2-test,pCoat colourandTYRP1(LOD=29.50).

The availability of accurate linkage maps is an important step for the localization of genetic variants of interest. However, most studies in livestock assume the published map is applicable in their population despite the large differences between the breeds of a species. A region of sheep Chromosome 1 was previously identified as providing evidence for a marker order inconsistent with the published linkage map. In this study the identified region was investigated in more detail. Four microsatellite markers covering the central 5 cM of the inconsistent region and two flanking markers were genotyped in three sheep breeds, a commercial population (Charollais), an experimental population (Scottish Blackface), and a feral population (Soay). With the inclusion of the published linkage map, this provided evidence for three different marker orders across four sheep populations. Evidence for selection in this region was investigated using both a single-point allelic competition model and a multipoint allele-sharing statistic. Only the Charollais population provided evidence for selection, with significant transmission bias observed at marker BM7145. The implications of variation in linkage maps on the design and analysis of fine-mapping studies are discussed.

The evolutionary changes that occur over a small number of generations in natural populations often run counter to what is expected on the basis of the heritability of traits and the selective forces acting upon them. In Soay sheep, dark coat color is associated with large size, which is heritable and positively correlated with fitness, yet the frequency of dark sheep has decreased. This unexpected microevolutionary trend is explained by genetic linkage between the causal mutation underlying the color polymorphism and quantitative trait loci with antagonistic effects on size and fitness. As a consequence, homozygous dark sheep are large, but have reduced fitness relative to phenotypically indistinguishable dark heterozygotes and light sheep. This result demonstrates the importance of understanding the genetic basis of fitness variation when making predictions about the microevolutionary consequences of selection.

To uncover the function of and interplay between the mammalian cytosine modifications 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC), new techniques and advances in current technology are needed. To this end, we have developed oxidative bisulfite sequencing (oxBS-seq), which can quantitatively locate 5mC and 5hmC marks at single-base resolution in genomic DNA. In bisulfite sequencing (BS-seq), both 5mC and 5hmC are read as cytosines and thus cannot be discriminated; however, in oxBS-seq, specific oxidation of 5hmC to 5-formylcytosine (5fC) and conversion of the newly formed 5fC to uracil (under bisulfite conditions) means that 5hmC can be discriminated from 5mC. A positive readout of actual 5mC is gained from a single oxBS-seq run, and 5hmC levels are inferred by comparison with a BS-seq run. Here we describe an optimized second-generation protocol that can be completed in 2 d.

Background The domestic pig is known as an excellent model for human immunology and the two species share many pathogens. Susceptibility to infectious disease is one of the major constraints on swine performance, yet the structure and function of genes comprising the pig immunome are not well-characterized. The completion of the pig genome provides the opportunity to annotate the pig immunome, and compare and contrast pig and human immune systems. Results The Immune Response Annotation Group (IRAG) used computational curation and manual annotation of the swine genome assembly 10.2 (Sscrofa10.2) to refine the currently available automated annotation of 1,369 immunity-related genes through sequence-based comparison to genes in other species. Within these genes, we annotated 3,472 transcripts. Annotation provided evidence for gene expansions in several immune response families, and identified artiodactyl-specific expansions in the cathelicidin and type 1 Interferon families. We found gene duplications for 18 genes, including 13 immune response genes and five non-immune response genes discovered in the annotation process. Manual annotation provided evidence for many new alternative splice variants and 8 gene duplications. Over 1,100 transcripts without porcine sequence evidence were detected using cross-species annotation. We used a functional approach to discover and accurately annotate porcine immune response genes. A co-expression clustering analysis of transcriptomic data from selected experimental infections or immune stimulations of blood, macrophages or lymph nodes identified a large cluster of genes that exhibited a correlated positive response upon infection across multiple pathogens or immune stimuli. Interestingly, this gene cluster (cluster 4) is enriched for known general human immune response genes, yet contains many un-annotated porcine genes. A phylogenetic analysis of the encoded proteins of cluster 4 genes showed that 15% exhibited an accelerated evolution as compared to 4.1% across the entire genome. Conclusions This extensive annotation dramatically extends the genome-based knowledge of the molecular genetics and structure of a major portion of the porcine immunome. Our complementary functional approach using co-expression during immune response has provided new putative immune response annotation for over 500 porcine genes. Our phylogenetic analysis of this core immunome cluster confirms rapid evolutionary change in this set of genes, and that, as in other species, such genes are important components of the pig’s adaptation to pathogen challenge over evolutionary time. These comprehensive and integrated analyses increase the value of the porcine genome sequence and provide important tools for global analyses and data-mining of the porcine immune response.

Human and mouse monocyte can be divided into two different subpopulations based on surface marker expression: CD14/16 and Ly6C/CX3CR1, respectively. Monocyte subpopulations in the pig were identified based on reciprocal expression of CD14 and the scavenger receptor CD163. The two populations, CD14hi-CD163low and CD14low-CD163hi, show approximately equal abundance in the steady-state. Culture of pig PBMCs in CSF1 indicates that the two populations are a maturation series controlled by this growth factor. Gene expression in pig monocyte subpopulations was profiled using the newly developed and annotated pig whole genome snowball microarray. Previous studies have suggested a functional equivalence between human and mouse subsets, but certain genes such as CD36, CLEC4E, or TREM-1 showed human-specific expression. The same genes were expressed selectively in pig monocyte subsets. However, the profiles suggest that the pig CD14low-CD163high cells are actually equivalent to intermediate human monocytes, and there is no CD14− CD16+ “nonclassical” population. The results are discussed in terms of the relevance of the pig as a model for understanding human monocyte function.

The G-quadruplex is an alternative DNA structural motif that is considered to be functionally important in the mammalian genome for transcriptional regulation, DNA replication and genome stability, but the nature and distribution of G-quadruplexes across the genome remains elusive. Here, we address the hypothesis that G-quadruplex structures exist within double-stranded genomic DNA and can be explicitly identified using a G-quadruplex-specific probe. An engineered antibody is employed to enrich for DNA containing G-quadruplex structures, followed by deep sequencing to detect and map G-quadruplexes at high resolution in genomic DNA from human breast adenocarcinoma cells. Our high sensitivity structure-based pull-down strategy enables the isolation of genomic DNA fragments bearing single, as well as multiple G-quadruplex structures. Stable G-quadruplex structures are found in sub-telomeres, gene bodies and gene regulatory regions. For a sample of identified target genes, we show that G-quadruplex-stabilizing ligands can modulate transcription. These results confirm the existence of G-quadruplex structures and their persistence in human genomic DNA.

Background The draft genome of the domestic pig (Sus scrofa) has recently been published permitting refined analysis of the transcriptome. Pig breeds have been reported to differ in their resistance to infectious disease. In this study we examine whether there are corresponding differences in gene expression in innate immune cells Results We demonstrate that macrophages can be harvested from three different compartments of the pig (lungs, blood and bone-marrow), cryopreserved and subsequently recovered and differentiated in CSF-1. We have performed surface marker analysis and gene expression profiling on macrophages from these compartments, comparing twenty-five animals from five different breeds and their response to lipopolysaccharide. The results provide a clear distinction between alveolar macrophages (AM) and monocyte-derived (MDM) and bone-marrow-derived macrophages (BMDM). In particular, the lung macrophages express the growth factor, FLT1 and its ligand, VEGFA at high levels, suggesting a distinct pathway of growth regulation. Relatively few genes showed breed-specific differential expression, notably CXCR2 and CD302 in alveolar macrophages. In contrast, there was substantial inter-individual variation between pigs within breeds, mostly affecting genes annotated as being involved in immune responses. Conclusions Pig macrophages more closely resemble human, than mouse, in their set of macrophage-expressed and LPS-inducible genes. Future research will address whether inter-individual variation in macrophage gene expression is heritable, and might form the basis for selective breeding for disease resistance.

Background This work describes the first genome-wide analysis of the transcriptional landscape of the pig. A new porcine Affymetrix expression array was designed in order to provide comprehensive coverage of the known pig transcriptome. The new array was used to generate a genome-wide expression atlas of pig tissues derived from 62 tissue/cell types. These data were subjected to network correlation analysis and clustering. Results The analysis presented here provides a detailed functional clustering of the pig transcriptome where transcripts are grouped according to their expression pattern, so one can infer the function of an uncharacterized gene from the company it keeps and the locations in which it is expressed. We describe the overall transcriptional signatures present in the tissue atlas, where possible assigning those signatures to specific cell populations or pathways. In particular, we discuss the expression signatures associated with the gastrointestinal tract, an organ that was sampled at 15 sites along its length and whose biology in the pig is similar to human. We identify sets of genes that define specialized cellular compartments and region-specific digestive functions. Finally, we performed a network analysis of the transcription factors expressed in the gastrointestinal tract and demonstrate how they sub-divide into functional groups that may control cellular gastrointestinal development. Conclusions As an important livestock animal with a physiology that is more similar than mouse to man, we provide a major new resource for understanding gene expression with respect to the known physiology of mammalian tissues and cells. The data and analyses are available on the websites http://biogps.org and http://www.macrophages.com/pig-atlas.

Mouse bone marrow-derived macrophages (BMDM) grown in M-CSF (CSF-1) have been used widely in studies of macrophage biology and the response to TLR agonists. We investigated whether similar cells could be derived from the domestic pig using human rCSF-1 and whether porcine macrophages might represent a better model of human macrophage biology. Cultivation of pig bone marrow cells for 5–7 d in presence of human rCSF-1 generated a pure population of BMDM that expressed the usual macrophage markers (CD14, CD16, and CD172a), were potent phagocytic cells, and produced TNF in response to LPS. Pig BMDM could be generated from bone marrow cells that had been stored frozen and thawed so that multiple experiments can be performed on samples from a single animal. Gene expression in pig BMDM from outbred animals responding to LPS was profiled using Affymetrix microarrays. The temporal cascade of inducible and repressible genes more closely resembled the known responses of human than mouse macrophages, sharing with humans the regulation of genes involved in tryptophan metabolism (IDO, KYN), lymphoattractant chemokines (CCL20, CXCL9, CXCL11, CXCL13), and the vitamin D3-converting enzyme, Cyp27B1. Conversely, in common with published studies of human macrophages, pig BMDM did not strongly induce genes involved in arginine metabolism, nor did they produce NO. These results establish pig BMDM as an alternative tractable model for the study of macrophage transcriptional control.

Understanding the genetic architecture of phenotypic variation in natural populations is a fundamental goal of evolutionary genetics. Wild Soay sheep (Ovis aries) have an inherited polymorphism for horn morphology in both sexes, controlled by a single autosomal locus, Horns. The majority of males have large normal horns, but a small number have vestigial, deformed horns, known as scurs; females have either normal horns, scurs or no horns (polled). Given that scurred males and polled females have reduced fitness within each sex, it is counterintuitive that the polymorphism persists within the population. Therefore, identifying the genetic basis of horn type will provide a vital foundation for understanding why the different morphs are maintained in the face of natural selection. We conducted a genome-wide association study using ∼36000 single nucleotide polymorphisms (SNPs) and determined the main candidate for Horns as RXFP2, an autosomal gene with a known involvement in determining primary sex characters in humans and mice. Evidence from additional SNPs in and around RXFP2 supports a new model of horn-type inheritance in Soay sheep, and for the first time, sheep with the same horn phenotype but different underlying genotypes can be identified. In addition, RXFP2 was shown to be an additive quantitative trait locus (QTL) for horn size in normal-horned males, accounting for up to 76% of additive genetic variation in this trait. This finding contrasts markedly from genome-wide association studies of quantitative traits in humans and some model species, where it is often observed that mapped loci only explain a modest proportion of the overall genetic variation.

[Maloney et al . (2009)][1] propose that the decline in frequency of dark Soay sheep on St Kilda, first reported by [Gratten et al . (2008)][2], is due to climate change. As supporting evidence, they report that the proportion of dark sheep is negatively correlated with average minimum winter

SummaryHistorical information can be used, in addition to pedigree, traits and genotypes, to map quantitative trait locus (QTL) in general populations via maximum likelihood estimation of variance components. This analysis is known as linkage disequilibrium (LD) and linkage mapping, because it exploits both linkage in families and LD at the population level. The search for QTL in the wild population of Soay sheep on St. Kilda is a proof of principle. We analysed the data from a previous study and confirmed some of the QTLs reported. The most striking result was the confirmation of a QTL affecting birth weight that had been reported using association tests but not when using linkage-based analyses.

Bridging the genotype-phenotype gap for traits of ecological and evolutionary importance in natural populations can provide a novel insight into the origin and maintenance of variation. Here, we identify the gene and putative causal mutations underlying a recessive colour pattern phenotype ('self' or uniform colour) in a wild population of primitive Soay sheep. We targeted the agouti signalling protein (ASIP) gene, a positional candidate based on previous study that mapped the Coat pattern locus to a presumptive region on chromosome 13. We found evidence for three recessive mutations, including two functional changes in the coding sequence and a putative third cis-regulatory mutation that inactivates the promoter. These mutations define up to five haplotypes in Soays, which collectively explained the coat pattern in all but one member of a complex multi-generational pedigree containing 621 genotyped individuals. The functional mutations are in strong linkage disequilibrium in the study population, and are identical to those known to underlie the self phenotype in domestic sheep. This is indicative of a recent (and simultaneous) origin in Soay sheep, possibly as a consequence of past interbreeding with modern domestic breeds. This is only the second study in which ASIP has been linked to variation in pigmentation in a natural population. Knowledge of the genetic basis of self-colour pattern in Soay sheep, and the recognition that several mutations are segregating in the population, will aid future studies investigating the role of selection in the maintenance of the polymorphism.