Systematic biases in DNA copy number originate from isolation procedures

Background

The ability to accurately detect DNA copy number variation in both a sensitive and quantitative manner is important in many research areas. However, genome-wide DNA copy number analyses are complicated by variations in detection signal.

Results

While GC content has been used to correct for this, here we show that coverage biases are tissue-specific and independent of the detection method as demonstrated by next-generation sequencing and array CGH. Moreover, we show that DNA isolation stringency affects the degree of equimolar coverage and that the observed biases coincide with chromatin characteristics like gene expression, genomic isochores, and replication timing.

Conclusion

These results indicate that chromatin organization is a main determinant for differential DNA retrieval. These findings are highly relevant for germline and somatic DNA copy number variation analyses.     

Global Analysis of DNA Methylation Variation in Adipose Tissue from Twins Reveals Links to Disease-Associated Variants in Distal Regulatory Elements

Epigenetic modifications such as DNA methylation play a key role in gene regulation and disease susceptibility. However, little is known about the genome-wide frequency, localization, and function of methylation variation and how it is regulated by genetic and environmental factors. We utilized the Multiple Tissue Human Expression Resource (MuTHER) and generated Illumina 450K adipose methylome data from 648 twins. We found that individual CpGs had low variance and that variability was suppressed in promoters. We noted that DNA methylation variation was highly heritable (h²_median = 0.34) and that shared environmental effects correlated with metabolic phenotype-associated CpGs. Analysis of methylation quantitative-trait loci (metQTL) revealed that 28% of CpGs were associated with nearby SNPs, and when overlapping them with adipose expression quantitative-trait loci (eQTL) from the same individuals, we found that 6% of the loci played a role in regulating both gene expression and DNA methylation. These associations were bidirectional, but there were pronounced negative associations for promoter CpGs. Integration of metQTL with adipose reference epigenomes and disease associations revealed significant enrichment of metQTL overlapping metabolic-trait or disease loci in enhancers (the strongest effects were for high-density lipoprotein cholesterol and body mass index [BMI]). We followed up with the BMI SNP rs713586, a cg01884057 metQTL that overlaps an enhancer upstream of ADCY3, and used bisulphite sequencing to refine this region. Our results showed widespread population invariability yet sequence dependence on adipose DNA methylation but that incorporating maps of regulatory elements aid in linking CpG variation to gene regulation and disease risk in a tissue-dependent manner.     

Detecting and Characterizing Genomic Signatures of Positive Selection in Global Populations

Natural selection is a significant force that shapes the architecture of the human genome and introduces diversity across global populations. The question of whether advantageous mutations have arisen in the human genome as a result of single or multiple mutation events remains unanswered except for the fact that there exist a handful of genes such as those that confer lactase persistence, affect skin pigmentation, or cause sickle cell anemia. We have developed a long-range-haplotype method for identifying genomic signatures of positive selection to complement existing methods, such as the integrated haplotype score (iHS) or cross-population extended haplotype homozygosity (XP-EHH), for locating signals across the entire allele frequency spectrum. Our method also locates the founder haplotypes that carry the advantageous variants and infers their corresponding population frequencies. This presents an opportunity to systematically interrogate the whole human genome whether a selection signal shared across different populations is the consequence of a single mutation process followed subsequently by gene flow between populations or of convergent evolution due to the occurrence of multiple independent mutation events either at the same variant or within the same gene. The application of our method to data from 14 populations across the world revealed that positive-selection events tend to cluster in populations of the same ancestry. Comparing the founder haplotypes for events that are present across different populations revealed that convergent evolution is a rare occurrence and that the majority of shared signals stem from the same evolutionary event.     

FLT3 Mutations in Early T-Cell Precursor ALL Characterize a Stem Cell Like Leukemia and Imply the Clinical Use of Tyrosine Kinase Inhibitors

Early T-cell precursor acute lymphoblastic leukemia (ETP-ALL) has been identified as high-risk subgroup of acute T-lymphoblastic leukemia (T-ALL) with a high rate of FLT3-mutations in adults. To unravel the underlying pathomechanisms and the clinical course we assessed molecular alterations and clinical characteristics in a large cohort of ETP-ALL (n = 68) in comparison to non-ETP T-ALL adult patients. Interestingly, we found a high rate of FLT3-mutations in ETP-ALL samples (n = 24, 35%). Furthermore, FLT3 mutated ETP-ALL was characterized by a specific immunophenotype (CD2+/CD5-/CD13+/CD33-), a distinct gene expression pattern (aberrant expression of IGFBP7, WT1, GATA3) and mutational status (absence of NOTCH1 mutations and a low frequency, 21%, of clonal TCR rearrangements). The observed low GATA3 expression and high WT1 expression in combination with lack of NOTCH1 mutations and a low rate of TCR rearrangements point to a leukemic transformation at the pluripotent prothymocyte stage in FLT3 mutated ETP-ALL. The clinical outcome in ETP-ALL patients was poor, but encouraging in those patients with allogeneic stem cell transplantation (3-year OS: 74%). To further explore the efficacy of targeted therapies, we demonstrate that T-ALL cell lines transfected with FLT3 expression constructs were particularly sensitive to tyrosine kinase inhibitors. In conclusion, FLT3 mutated ETP-ALL defines a molecular distinct stem cell like leukemic subtype. These data warrant clinical studies with the implementation of FLT3 inhibitors in addition to early allogeneic stem cell transplantation for this high risk subgroup.     

Scanning electron microscopy observation of host entry by two brown algae endophytic in Laminaria saccharina (Laminariales, Phaeophyceae)

The germination of spores of the endophytic brown algae Laminariocolax aecidioides and Laminarionema elsbetiae (Ectocarpales sensu lato) on their host Laminaria saccharina was studied using scanning electron microscopy. Zoids were the infective units in both taxa. Despite a large difference in size between spores of the species studied, the initial steps of infection were similar. Zoids settled on the host surface, and during secretion of adhesive material from their anterior end, the spores became elongated to rod-shaped and stood erect on the host. A single germ tube developed at the proximal end and penetrated the intact surface of the host. Sharp edges around the entrance hole and the absence of inward deformation of the host surface around the settled zoids suggest an enzymatic rather than a mechanical penetration mechanism.     

A genotype calling algorithm for the Illumina BeadArray platform

MOTIVATION: Large-scale genotyping relies on the use of unsupervised automated calling algorithms to assign genotypes to hybridization data. A number of such calling algorithms have been recently established for the Affymetrix GeneChip genotyping technology. Here, we present a fast and accurate genotype calling algorithm for the Illumina BeadArray genotyping platforms. As the technology moves towards assaying millions of genetic polymorphisms simultaneously, there is a need for an integrated and easy-to-use software for calling genotypes. RESULTS: We have introduced a model-based genotype calling algorithm which does not rely on having prior training data or require computationally intensive procedures. The algorithm can assign genotypes to hybridization data from thousands of individuals simultaneously and pools information across multiple individuals to improve the calling. The method can accommodate variations in hybridization intensities which result in dramatic shifts of the position of the genotype clouds by identifying the optimal coordinates to initialize the algorithm. By incorporating the process of perturbation analysis, we can obtain a quality metric measuring the stability of the assigned genotype calls. We show that this quality metric can be used to identify SNPs with low call rates and accuracy. AVAILABILITY: The C++ executable for the algorithm described here is available by request from the authors.     

The architecture of gene regulatory variation across multiple human tissues: the MuTHER study

While there have been studies exploring regulatory variation in one or more tissues, the complexity of tissue-specificity in multiple primary tissues is not yet well understood. We explore in depth the role of cis-regulatory variation in three human tissues: lymphoblastoid cell lines (LCL), skin, and fat. The samples (156 LCL, 160 skin, 166 fat) were derived simultaneously from a subset of well-phenotyped healthy female twins of the MuTHER resource. We discover an abundance of cis-eQTLs in each tissue similar to previous estimates (858 or 4.7% of genes). In addition, we apply factor analysis (FA) to remove effects of latent variables, thus more than doubling the number of our discoveries (1,822 eQTL genes). The unique study design (Matched Co-Twin Analysis--MCTA) permits immediate replication of eQTLs using co-twins (93%-98%) and validation of the considerable gain in eQTL discovery after FA correction. We highlight the challenges of comparing eQTLs between tissues. After verifying previous significance threshold-based estimates of tissue-specificity, we show their limitations given their dependency on statistical power. We propose that continuous estimates of the proportion of tissue-shared signals and direct comparison of the magnitude of effect on the fold change in expression are essential properties that jointly provide a biologically realistic view of tissue-specificity. Under this framework we demonstrate that 30% of eQTLs are shared among the three tissues studied, while another 29% appear exclusively tissue-specific. However, even among the shared eQTLs, a substantial proportion (10%-20%) have significant differences in the magnitude of fold change between genotypic classes across tissues. Our results underline the need to account for the complexity of eQTL tissue-specificity in an effort to assess consequences of such variants for complex traits.     

Novel loci for adiponectin levels and their influence on type 2 diabetes and metabolic traits: a multi-ethnic meta-analysis of 45,891 individuals

Circulating levels of adiponectin, a hormone produced predominantly by adipocytes, are highly heritable and are inversely associated with type 2 diabetes mellitus (T2D) and other metabolic traits. We conducted a meta-analysis of genome-wide association studies in 39,883 individuals of European ancestry to identify genes associated with metabolic disease. We identified 8 novel loci associated with adiponectin levels and confirmed 2 previously reported loci (P = 4.5×10⁻⁸–1.2×10⁻⁴³). Using a novel method to combine data across ethnicities (N = 4,232 African Americans, N = 1,776 Asians, and N = 29,347 Europeans), we identified two additional novel loci. Expression analyses of 436 human adipocyte samples revealed that mRNA levels of 18 genes at candidate regions were associated with adiponectin concentrations after accounting for multiple testing (p<3×10⁻⁴). We next developed a multi-SNP genotypic risk score to test the association of adiponectin decreasing risk alleles on metabolic traits and diseases using consortia-level meta-analytic data. This risk score was associated with increased risk of T2D (p = 4.3×10⁻³, n = 22,044), increased triglycerides (p = 2.6×10⁻¹⁴, n = 93,440), increased waist-to-hip ratio (p = 1.8×10⁻⁵, n = 77,167), increased glucose two hours post oral glucose tolerance testing (p = 4.4×10⁻³, n = 15,234), increased fasting insulin (p = 0.015, n = 48,238), but with lower in HDL-cholesterol concentrations (p = 4.5×10⁻¹³, n = 96,748) and decreased BMI (p = 1.4×10⁻⁴, n = 121,335). These findings identify novel genetic determinants of adiponectin levels, which, taken together, influence risk of T2D and markers of insulin resistance.     

The Function of the Chemokine Receptor CXCR6 in the T Cell Response of Mice against Listeria monocytogenes

The chemokine receptor CXCR6 is expressed on different T cell subsets and up-regulated following T cell activation. CXCR6 has been implicated in the localization of cells to the liver due to the constitutive expression of its ligand CXCL16 on liver sinusoidal endothelial cells. Here, we analyzed the role of CXCR6 in CD8⁺ T cell responses to infection of mice with Listeria monocytogenes. CD8⁺ T cells responding to listerial antigens acquired high expression levels of CXCR6. However, deficiency of mice in CXCR6 did not impair control of the L. monocytogenes infection. CXCR6-deficient mice were able to generate listeria-specific CD4⁺ and CD8⁺ T cell responses and showed accumulation of T cells in the infected liver. In transfer assays, we detected reduced accumulation of listeria-specific CXCR6-deficient CD8⁺ T cells in the liver at early time points post infection. Though, CXCR6 was dispensable at later time points of the CD8⁺ T cell response. When transferred CD8⁺ T cells were followed for extended time periods, we observed a decline in CXCR6-deficient CD8⁺ T cells. The manifestation of this cell loss depended on the tissue analyzed. In conclusion, our results demonstrate that CXCR6 is not required for the formation of a T cell response to L. monocytogenes and for the accumulation of T cells in the infected liver but CXCR6 appears to influence long-term survival and tissue distribution of activated cells.