Mutation, methylation, and gene expression profiles in dup(1q)-positive pediatric B-cell precursor acute lymphoblastic leukemia.

High-throughput sequencing was applied to investigate the mutation/methylation patterns on 1q and gene expression profiles in pediatric B-cell precursor acute lymphoblastic leukemia (BCP ALL) with/without (w/wo) dup(1q). Sequencing of the breakpoint regions and all exons on 1q in seven dup(1q)-positive cases revealed non-synonymous somatic single nucleotide variants (SNVs) in BLZF1, FMN2, KCNT2, LCE1C, NES, and PARP1. Deep sequencing of these in a validation cohort w (n = 17)/wo (n = 94) dup(1q) revealed similar SNV frequencies in the two groups (47% vs. 35%; P = 0.42). Only 0.6% of the 36,259 CpGs on 1q were differentially methylated between cases w (n = 14)/wo (n = 13) dup(1q). RNA sequencing of high hyperdiploid (HeH) and t(1;19)(q23;p13)-positive cases w (n = 14)/wo (n = 52) dup(1q) identified 252 and 424 differentially expressed genes, respectively; only seven overlapped. Of the overexpressed genes in the HeH and t(1;19) groups, 23 and 31%, respectively, mapped to 1q; 60-80% of these encode nucleic acid/protein binding factors or proteins with catalytic activity. We conclude that the pathogenetically important consequence of dup(1q) in BCP ALL is a gene-dosage effect, with the deregulated genes differing between genetic subtypes, but involving similar molecular functions, biological processes, and protein classes.

Inferring the evolution of the major histocompatibility complex of wild pigs and peccaries using hybridisation DNA capture-based sequencing.

The major histocompatibility complex (MHC) is a key genomic model region for understanding the evolution of gene families and the co-evolution between host and pathogen. To date, MHC studies have mostly focused on species from major vertebrate lineages. The evolution of MHC classical (Ia) and non-classical (Ib) genes in pigs has attracted interest because of their antigen presentation roles as part of the adaptive immune system. The pig family Suidae comprises over 18 extant species (mostly wild), but only the domestic pig has been extensively sequenced and annotated. To address this, we used a DNA-capture approach, with probes designed from the domestic pig genome, to generate MHC data for 11 wild species of pigs and their closest living family, Tayassuidae. The approach showed good efficiency for wild pigs (~80% reads mapped, ~87× coverage), compared to tayassuids (~12% reads mapped, ~4× coverage). We retrieved 145 MHC loci across both families. Phylogenetic analyses show that the class Ia and Ib genes underwent multiple duplications and diversifications before suids and tayassuids diverged from their common ancestor. The histocompatibility genes mostly form orthologous groups and there is genetic differentiation for most of these genes between Eurasian and sub-Saharan African wild pigs. Tests of selection showed that the peptide-binding region of class Ib genes was under positive selection. These findings contribute to better understanding of the evolutionary history of the MHC, specifically, the class I genes, and provide useful data for investigating the immune response of wild populations against pathogens.

SPlinted Ligation Adapter Tagging (SPLAT), a novel library preparation method for whole genome bisulphite sequencing.

Sodium bisulphite treatment of DNA combined with next generation sequencing (NGS) is a powerful combination for the interrogation of genome-wide DNA methylation profiles. Library preparation for whole genome bisulphite sequencing (WGBS) is challenging due to side effects of the bisulphite treatment, which leads to extensive DNA damage. Recently, a new generation of methods for bisulphite sequencing library preparation have been devised. They are based on initial bisulphite treatment of the DNA, followed by adaptor tagging of single stranded DNA fragments, and enable WGBS using low quantities of input DNA. In this study, we present a novel approach for quick and cost effective WGBS library preparation that is based on splinted adaptor tagging (SPLAT) of bisulphite-converted single-stranded DNA. Moreover, we validate SPLAT against three commercially available WGBS library preparation techniques, two of which are based on bisulphite treatment prior to adaptor tagging and one is a conventional WGBS method.

DNA methylome analysis of acute lymphoblastic leukemia cells reveals stochastic de novo DNA methylation in CpG islands.

AIM: To identify regions of aberrant DNA methylation in acute lymphoblastic leukemia (ALL) cells of different subtypes on a genome-wide scale. MATERIALS & METHODS: Whole-genome bisulfite sequencing (WGBS) was used to determine the DNA methylation levels in cells from four pediatric ALL patients of different subtypes. The findings were confirmed by 450k DNA methylation arrays in a large patient set. RESULTS: Compared with mature B or T cells WGBS detected on average 82,000 differentially methylated regions per patient. Differentially methylated regions are enriched to CpG poor regions, active enhancers and transcriptional start sites. We also identified approximately 8000 CpG islands with variable intermediate DNA methylation that seems to occur as a result of stochastic de novo methylation. CONCLUSION: WGBS provides an unbiased view and novel insights into the DNA methylome of ALL cells.

Mitogenomic analysis of a 50-generation chicken pedigree reveals a rapid rate of mitochondrial evolution and evidence for paternal mtDNA inheritance.

Mitochondrial genomes represent a valuable source of data for evolutionary research, but studies of their short-term evolution have typically been limited to invertebrates, humans and laboratory organisms. Here we present a detailed study of 12 mitochondrial genomes that span a total of 385 transmissions in a well-documented 50-generation pedigree in which two lineages of chickens were selected for low and high juvenile body weight. These data allowed us to test the hypothesis of time-dependent evolutionary rates and the assumption of strict maternal mitochondrial transmission, and to investigate the role of mitochondrial mutations in determining phenotype. The identification of a non-synonymous mutation in ND4L and a synonymous mutation in CYTB, both novel mutations in Gallus, allowed us to estimate a molecular rate of 3.13 × 10(-7) mutations/site/year (95% confidence interval 3.75 × 10(-8)-1.12 × 10(-6)). This is substantially higher than avian rate estimates based upon fossil calibrations. Ascertaining which of the two novel mutations was present in an additional 49 individuals also revealed an instance of paternal inheritance of mtDNA. Lastly, an association analysis demonstrated that neither of the point mutations was strongly associated with the phenotypic differences between the two selection lines. Together, these observations reveal the highly dynamic nature of mitochondrial evolution over short time periods.

The mutational landscape in pediatric acute lymphoblastic leukemia deciphered by whole genome sequencing.

Genomic characterization of pediatric acute lymphoblastic leukemia (ALL) has identified distinct patterns of genes and pathways altered in patients with well-defined genetic aberrations. To extend the spectrum of known somatic variants in ALL, we performed whole genome and transcriptome sequencing of three B-cell precursor patients, of which one carried the t(12;21)ETV6-RUNX1 translocation and two lacked a known primary genetic aberration, and one T-ALL patient. We found that each patient had a unique genome, with a combination of well-known and previously undetected genomic aberrations. By targeted sequencing in 168 patients, we identified KMT2D and KIF1B as novel putative driver genes. We also identified a putative regulatory non-coding variant that coincided with overexpression of the growth factor MDK. Our results contribute to an increased understanding of the biological mechanisms that lead to ALL and suggest that regulatory variants may be more important for cancer development than recognized to date. The heterogeneity of the genetic aberrations in ALL renders whole genome sequencing particularly well suited for analysis of somatic variants in both research and diagnostic applications.

Accurate detection of subclonal single nucleotide variants in whole genome amplified and pooled cancer samples using HaloPlex target enrichment.

BACKGROUND: Target enrichment and resequencing is a widely used approach for identification of cancer genes and genetic variants associated with diseases. Although cost effective compared to whole genome sequencing, analysis of many samples constitutes a significant cost, which could be reduced by pooling samples before capture. Another limitation to the number of cancer samples that can be analyzed is often the amount of available tumor DNA. We evaluated the performance of whole genome amplified DNA and the power to detect subclonal somatic single nucleotide variants in non-indexed pools of cancer samples using the HaloPlex technology for target enrichment and next generation sequencing. RESULTS: We captured a set of 1528 putative somatic single nucleotide variants and germline SNPs, which were identified by whole genome sequencing, with the HaloPlex technology and sequenced to a depth of 792-1752. We found that the allele fractions of the analyzed variants are well preserved during whole genome amplification and that capture specificity or variant calling is not affected. We detected a large majority of the known single nucleotide variants present uniquely in one sample with allele fractions as low as 0.1 in non-indexed pools of up to ten samples. We also identified and experimentally validated six novel variants in the samples included in the pools. CONCLUSION: Our work demonstrates that whole genome amplified DNA can be used for target enrichment equally well as genomic DNA and that accurate variant detection is possible in non-indexed pools of cancer samples. These findings show that analysis of a large number of samples is feasible at low cost, even when only small amounts of DNA is available, and thereby significantly increases the chances of indentifying recurrent mutations in cancer samples.

Genome-wide signatures of differential DNA methylation in pediatric acute lymphoblastic leukemia.

BACKGROUND: Although aberrant DNA methylation has been observed previously in acute lymphoblastic leukemia (ALL), the patterns of differential methylation have not been comprehensively determined in all subtypes of ALL on a genome-wide scale. The relationship between DNA methylation, cytogenetic background, drug resistance and relapse in ALL is poorly understood. RESULTS: We surveyed the DNA methylation levels of 435,941 CpG sites in samples from 764 children at diagnosis of ALL and from 27 children at relapse. This survey uncovered four characteristic methylation signatures. First, compared with control blood cells, the methylomes of ALL cells shared 9,406 predominantly hypermethylated CpG sites, independent of cytogenetic background. Second, each cytogenetic subtype of ALL displayed a unique set of hyper- and hypomethylated CpG sites. The CpG sites that constituted these two signatures differed in their functional genomic enrichment to regions with marks of active or repressed chromatin. Third, we identified subtype-specific differential methylation in promoter and enhancer regions that were strongly correlated with gene expression. Fourth, a set of 6,612 CpG sites was predominantly hypermethylated in ALL cells at relapse, compared with matched samples at diagnosis. Analysis of relapse-free survival identified CpG sites with subtype-specific differential methylation that divided the patients into different risk groups, depending on their methylation status. CONCLUSIONS: Our results suggest an important biological role for DNA methylation in the differences between ALL subtypes and in their clinical outcome after treatment.

Novel insights into the bovine polled phenotype and horn ontogenesis in Bovidae.

Despite massive research efforts, the molecular etiology of bovine polledness and the developmental pathways involved in horn ontogenesis are still poorly understood. In a recent article, we provided evidence for the existence of at least two different alleles at the Polled locus and identified candidate mutations for each of them. None of these mutations was located in known coding or regulatory regions, thus adding to the complexity of understanding the molecular basis of polledness. We confirm previous results here and exhaustively identify the causative mutation for the Celtic allele (PC) and four candidate mutations for the Friesian allele (PF). We describe a previously unreported eyelash-and-eyelid phenotype associated with regular polledness, and present unique histological and gene expression data on bovine horn bud differentiation in fetuses affected by three different horn defect syndromes, as well as in wild-type controls. We propose the ectopic expression of a lincRNA in PC/p horn buds as a probable cause of horn bud agenesis. In addition, we provide evidence for an involvement of OLIG2, FOXL2 and RXFP2 in horn bud differentiation, and draw a first link between bovine, ovine and caprine Polled loci. Our results represent a first and important step in understanding the genetic pathways and key process involved in horn bud differentiation in Bovidae.

Expression of carnitine palmitoyl-CoA transferase-1B is influenced by a cis-acting eQTL in two chicken lines selected for high and low body weight.

Carnitine palmitoyl-CoA transferase-1B is a mitochondrial enzyme in the fatty acid oxidation pathway. In a previous study, CPT1B was identified as differentially expressed in the hypothalamus of two lines of chickens established by long-term selection for high (HWS) or low (LWS) body weight. Mammals have three paralogs (CPT1a, b and c) while nonmammalian vertebrates only have two (CPT1A, B). CPT1A is expressed in liver and CPT1B in muscle. CPT1c is expressed in hypothalamus, where it regulates feeding and energy expenditure. We identified an intronic length polymorphism, fixed for different alleles in the two populations, and mapped the hitherto missing CPT1B locus in the chicken genome assembly, to the distal tip of chromosome 1p. Based on molecular phylogeny and gene synteny we suggest that chicken CPT1B is pro-orthologous of the mammalian CPT1c. Chicken CPT1B was differentially expressed in both muscle and hypothalamus but in opposite directions: higher levels in hypothalamus but lower levels in muscle in the HWS than in the LWS line. Using an advanced intercross population of the lines, we found CPT1B expression to be influenced by a cis-acting expression quantitative trait locus in muscle. The increased expression in hypothalamus and reduced expression in muscle is consistent with an increased food intake in the HWS line and at the same time reduced fatty acid oxidation in muscle yielding a net accumulation of energy intake and storage. The altered expression of CPT1B in hypothalamus and peripheral tissue is likely to be a mechanism contributing to the remarkable difference between lines.

Analysis of porcine MHC using microarrays.

The major histocompatibility complex (MHC) in Mammals is one of the most gene dense regions of the genome and contains the polymorphic histocompatibility gene families known to be involved in pathogen response and control of auto-immunity. The MHC is a complex genetic system that provides an interesting model system to study genome expression regulation and genetic diversity at the megabase scale. The pig MHC or SLA (Swine Leucocyte Antigen) complex spans 2.4 megabases and 151 loci have been annotated. We will review key results from previous RNA expression studies using microarrays containing probes specific to annotated loci within SLA and in addition present novel data obtained using high-density tiling arrays encompassing the whole SLA complex. We have focused on transcriptome modifications of porcine peripheral blood mononuclear cells stimulated with a mixture of phorbol myristate acetate and ionomycin known to activate B and T cell proliferation. Our results show that numerous loci mapping to the SLA complex are affected by the treatment. A general decreased level of expression for class I and II genes and an up-regulation of genes involved in peptide processing and transport were observed. Tiling array-based experiments contributed to refined gene annotations as presented for one SLA class I gene referred to as SLA-11. In conclusion, high-density tiling arrays can serve as an excellent tool to draw comprehensive transcription maps, and improve genome annotations for the SLA complex. We are currently studying their relevance to characterize SLA genetic diversity in combination with high throughput next generation sequencing.

Fine mapping and replication of QTL in outbred chicken advanced intercross lines.

BACKGROUND: Linkage mapping is used to identify genomic regions affecting the expression of complex traits. However, when experimental crosses such as F(2) populations or backcrosses are used to map regions containing a Quantitative Trait Locus (QTL), the size of the regions identified remains quite large, i.e. 10 or more Mb. Thus, other experimental strategies are needed to refine the QTL locations. Advanced Intercross Lines (AIL) are produced by repeated intercrossing of F(2) animals and successive generations, which decrease linkage disequilibrium in a controlled manner. Although this approach is seen as promising, both to replicate QTL analyses and fine-map QTL, only a few AIL datasets, all originating from inbred founders, have been reported in the literature. METHODS: We have produced a nine-generation AIL pedigree (n = 1529) from two outbred chicken lines divergently selected for body weight at eight weeks of age. All animals were weighed at eight weeks of age and genotyped for SNP located in nine genomic regions where significant or suggestive QTL had previously been detected in the F(2) population. In parallel, we have developed a novel strategy to analyse the data that uses both genotype and pedigree information of all AIL individuals to replicate the detection of and fine-map QTL affecting juvenile body weight. RESULTS: Five of the nine QTL detected with the original F(2) population were confirmed and fine-mapped with the AIL, while for the remaining four, only suggestive evidence of their existence was obtained. All original QTL were confirmed as a single locus, except for one, which split into two linked QTL. CONCLUSIONS: Our results indicate that many of the QTL, which are genome-wide significant or suggestive in the analyses of large intercross populations, are true effects that can be replicated and fine-mapped using AIL. Key factors for success are the use of large populations and powerful statistical tools. Moreover, we believe that the statistical methods we have developed to efficiently study outbred AIL populations will increase the number of organisms for which in-depth complex traits can be analyzed.

Genetic analysis of metabolic traits in an intercross between body weight-selected chicken lines.

A network of four interacting loci has been reported previously to influence growth in two lines of chickens divergently selected for body weight at 56 days of age. Located on chromosomes 3 (Growth4), 4 (Growth6), 7 (Growth9), and 20 (Growth12), they explained nearly half of the difference in body weight at selection age between the two lines. The original study reported effects on body weight and fat deposition, but no attempts were made to explore the effects of the network on other phenotypes measured in the F(2) population. In this study we conducted further analyses to evaluate the specific effects of the four-locus network on other metabolic traits as well as refining results from the original study by including a larger number of genetic markers in the quantitative trait locus (QTL) regions. We confirm the previously described effect of the epistatic network on body weight and show that the network increases the total amount of muscle and fat as well as the weight of the internal organs. The network as a whole did not change the relative content of any studied organs or tissues in the body. There was, however, a significant interaction between the loci on chromosomes 3 and 7 that changed the relative proportion of abdominal fat and breast muscle in the chicken by increasing abdominal fat weight without a corresponding increase in muscle mass.

Genetic analysis of an F(2) intercross between two chicken lines divergently selected for body-weight.

BACKGROUND: We have performed Quantitative Trait Loci (QTL) analysis of an F(2) intercross between two chicken lines divergently selected for juvenile body-weight. In a previous study 13 identified loci with effects on body-weight, only explained a small proportion of the large variation in the F(2) population. Epistatic interaction analysis however, indicated that a network of interacting loci with large effect contributed to the difference in body-weight of the parental lines. This previous analysis was, however, based on a sparse microsatellite linkage map and the limited coverage could have affected the main conclusions. Here we present a revised QTL analysis based on a high-density linkage map that provided a more complete coverage of the chicken genome. Furthermore, we utilized genotype data from ~13,000 SNPs to search the genome for potential selective sweeps that have occurred in the selected lines. RESULTS: We constructed a linkage map comprising 434 genetic markers, covering 31 chromosomes but leaving seven microchromosomes uncovered. The analysis showed that seven regions harbor QTL that influence growth. The pair-wise interaction analysis identified 15 unique QTL pairs and notable is that nine of those involved interactions with a locus on chromosome 7, forming a network of interacting loci. The analysis of ~13,000 SNPs showed that a substantial proportion of the genetic variation present in the founder population has been lost in either of the two selected lines since ~60% of the SNPs polymorphic among lines showed fixation in one of the lines. With the current marker coverage and QTL map resolution we did not observe clear signs of selective sweeps within QTL intervals. CONCLUSION: The results from the QTL analysis using the new improved linkage map are to a large extent in concordance with our previous analysis of this pedigree. The difference in body-weight between the parental chicken lines is caused by many QTL each with a small individual effect. Although the increased chromosomal marker coverage did not lead to the identification of additional QTL, we were able to refine the localization of QTL. The importance of epistatic interaction as a mechanism contributing significantly to the remarkable selection response was further strengthened because additional pairs of interacting loci were detected with the improved map.

A high-density SNP-based linkage map of the chicken genome reveals sequence features correlated with recombination rate.

The resolution of the chicken consensus linkage map has been dramatically improved in this study by genotyping 12,945 single nucleotide polymorphisms (SNPs) on three existing mapping populations in chicken: the Wageningen (WU), East Lansing (EL), and Uppsala (UPP) mapping populations. As many as 8599 SNPs could be included, bringing the total number of markers in the current consensus linkage map to 9268. The total length of the sex average map is 3228 cM, considerably smaller than previous estimates using the WU and EL populations, reflecting the higher quality of the new map. The current map consists of 34 linkage groups and covers at least 29 of the 38 autosomes. Sex-specific analysis and comparisons of the maps based on the three individual populations showed prominent heterogeneity in recombination rates between populations, but no significant heterogeneity between sexes. The recombination rates in the F(1) Red Jungle fowl/White Leghorn males and females were significantly lower compared with those in the WU broiler population, consistent with a higher recombination rate in purebred domestic animals under strong artificial selection. The recombination rate varied considerably among chromosomes as well as along individual chromosomes. An analysis of the sequence composition at recombination hot and cold spots revealed a strong positive correlation between GC-rich sequences and high recombination rates. The GC-rich cohesin binding sites in particular stood out from other GC-rich sequences with a 3.4-fold higher density at recombination hot spots versus cold spots, suggesting a functional relationship between recombination frequency and cohesin binding.

QTL analysis of body composition and metabolic traits in an intercross between chicken lines divergently selected for growth.

The high- and low-growth lines of chickens have been developed from a single founder population by divergent selection for body weight at 56 days of age for more than 40 generations. The two lines show a ninefold difference in body weight at selection age and several interesting correlated selection responses such as altered body composition and metabolic differences. We have generated a reciprocal intercross comprising >800 F2 birds. In a previous study, we reported the detection of 13 quantitative trait loci (QTLs) affecting growth. Here we report QTLs for body composition (fat deposition, muscle development), weight of internal organs, and metabolic traits (plasma concentrations of glucose, insulin, cholesterol, glucagon, triglycerides, and IGF-I). Most of the QTLs with convincing statistical support mapped in the vicinity of growth QTLs. One of the most interesting observations was that the type of reciprocal cross had highly significant effects on body weight at hatch and on plasma concentrations of glucose, cholesterol, insulin, and IGF-I, but it had no significant effect on body weight at 56 days of age. The reciprocal cross explained between 15 and 35% of the phenotypic variance for weight at hatch and for plasma concentrations of glucose and insulin. The observed pattern indicated that these effects were caused by maternal effects or by genetic differences in mitochondrial DNA.

Many QTLs with minor additive effects are associated with a large difference in growth between two selection lines in chickens.

Two growth-selected lines in chickens have been developed from a single founder population by divergent selection for body weight at 56 days of age. After more than 40 generations of selection they show a nine-fold difference in body weight at selection age and large differences in growth rate, appetite, fat deposition and metabolic characteristics. We have generated a large intercross between these lines comprising more than 800 F2 birds. QTL mapping revealed 13 loci affecting growth. The most striking observation was that the allele in the high weight line in all cases was associated with enhanced growth, but each locus explained only a small proportion of the phenotypic variance using a standard QTL model (1.3-3.1%). This result is in sharp contrast to our previous study where we reported that the two-fold difference in adult body size between the red junglefowl and White Leghorn domestic chickens is explained by a small number of QTLs with large additive effects. Furthermore, no QTLs for anorexia or antibody response were detected despite large differences for these traits between the founder lines. The result is an excellent example where a large phenotypic difference between populations occurs in the apparent absence of any single locus with large phenotypic effects. The study underscores the need for powerful experimental designs in genetic studies of multifactorial traits. No QTL at all would have reached genome-wide significance using a less powerful design (e.g. approx. 200 F2 individuals) regardless of the nine-fold phenotypic difference between the founder lines for the selected trait.