1. |
- Alfoeldi, Jessica, et al.
(author)
-
Comparative genomics as a tool to understand evolution and disease
- 2013
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 23:7, s. 1063-1068
-
Journal article (peer-reviewed)abstract
- When the human genome project started, the major challenge was how to sequence a 3 billion letter code in an organized and cost-effective manner. When completed, the project had laid the foundation for a huge variety of biomedical fields through the production of a complete human genome sequence, but also had driven the development of laboratory and analytical methods that could produce large amounts of sequencing data cheaply. These technological developments made possible the sequencing of many more vertebrate genomes, which have been necessary for the interpretation of the human genome. They have also enabled large-scale studies of vertebrate genome evolution, as well as comparative and human medicine. In this review, we give examples of evolutionary analysis using a wide variety of time frames-from the comparison of populations within a species to the comparison of species separated by at least 300 million years. Furthermore, we anticipate discoveries related to evolutionary mechanisms, adaptation, and disease to quickly accelerate in the coming years.
|
|
2. |
- Andersen, M. R., et al.
(author)
-
Comparative genomics of citric-acid-producing Aspergillus niger ATCC 1015 versus enzyme-producing CBS 513.88
- 2011
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 21:6, s. 885-897
-
Journal article (peer-reviewed)abstract
- The filamentous fungus Aspergillus niger exhibits great diversity in its phenotype. It is found globally, both as marine and terrestrial strains, produces both organic acids and hydrolytic enzymes in high amounts, and some isolates exhibit pathogenicity. Although the genome of an industrial enzyme-producing A. niger strain (CBS 513.88) has already been sequenced, the versatility and diversity of this species compel additional exploration. We therefore undertook wholegenome sequencing of the acidogenic A. niger wild-type strain (ATCC 1015) and produced a genome sequence of very high quality. Only 15 gaps are present in the sequence, and half the telomeric regions have been elucidated. Moreover, sequence information from ATCC 1015 was used to improve the genome sequence of CBS 513.88. Chromosome-level comparisons uncovered several genome rearrangements, deletions, a clear case of strain-specific horizontal gene transfer, and identification of 0.8 Mb of novel sequence. Single nucleotide polymorphisms per kilobase (SNPs/kb) between the two strains were found to be exceptionally high (average: 7.8, maximum: 160 SNPs/kb). High variation within the species was confirmed with exo-metabolite profiling and phylogenetics. Detailed lists of alleles were generated, and genotypic differences were observed to accumulate in metabolic pathways essential to acid production and protein synthesis. A transcriptome analysis supported up-regulation of genes associated with biosynthesis of amino acids that are abundant in glucoamylase A, tRNA-synthases, and protein transporters in the protein producing CBS 513.88 strain. Our results and data sets from this integrative systems biology analysis resulted in a snapshot of fungal evolution and will support further optimization of cell factories based on filamentous fungi
|
|
3. |
- Axelsson, Erik, et al.
(author)
-
Death of PRDM9 coincides with stabilization of the recombination landscape in the dog genome
- 2011
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 22:1, s. 51-63
-
Journal article (peer-reviewed)abstract
- Analysis of diverse eukaryotes has revealed that recombination events cluster in discrete genomic locations known as hotspots. In humans, a zinc-finger protein, PRDM9, is believed to initiate recombination in >40% of hotspots by binding to a specific DNA sequence motif. However, the PRDM9 coding sequence is disrupted in the dog genome assembly, raising questions regarding the nature and control of recombination in dogs. By analyzing the sequences of PRDM9 orthologs in a number of dog breeds and several carnivores, we show here that this gene was inactivated early in canid evolution. We next use patterns of linkage disequilibrium using more than 170,000 SNP markers typed in almost 500 dogs to estimate the recombination rates in the dog genome using a coalescent-based approach. Broad-scale recombination rates show good correspondence with an existing linkage-based map. Significant variation in recombination rate is observed on the fine scale, and we are able to detect over 4000 recombination hotspots with high confidence. In contrast to human hotspots, 40% of canine hotspots are characterized by a distinct peak in GC content. A comparative genomic analysis indicates that these peaks are present also as weaker peaks in the panda, suggesting that the hotspots have been continually reinforced by accelerated and strongly GC biased nucleotide substitutions, consistent with the long-term action of biased gene conversion on the dog lineage. These results are consistent with the loss of PRDM9 in canids, resulting in a greater evolutionary stability of recombination hotspots. The genetic determinants of recombination hotspots in the dog genome may thus reflect a fundamental process of relevance to diverse animal species.
|
|
4. |
- Backström, Niclas, et al.
(author)
-
The recombination landscape of the zebra finch Taeniopygia guttata genome
- 2010
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 20:4, s. 485-495
-
Journal article (peer-reviewed)abstract
- Understanding the causes and consequences of variation in the rate of recombination is essential since this parameter is considered to affect levels of genetic variability, the efficacy of selection, and the design of association and linkage mapping studies. However, there is limited knowledge about the factors governing recombination rate variation. We genotyped 1920 single nucleotide polymorphisms in a multigeneration pedigree of more than 1000 zebra finches (Taeniopygia guttata) to develop a genetic linkage map, and then we used these map data together with the recently available draft genome sequence of the zebra finch to estimate recombination rates in 1 Mb intervals across the genome. The average zebra finch recombination rate (1.5 cM/Mb) is higher than in humans, but significantly lower than in chicken. The local rates of recombination in chicken and zebra finch were only weakly correlated, demonstrating evolutionary turnover of the recombination landscape in birds. The distribution of recombination events was heavily biased toward ends of chromosomes, with a stronger telomere effect than so far seen in any organism. In fact, the recombination rate was as low as 0.1 cM/Mb in intervals up to 100 Mb long in the middle of the larger chromosomes. We found a positive correlation between recombination rate and GC content, as well as GC-rich sequence motifs. Levels of linkage disequilibrium (LD) were significantly higher in regions of low recombination, showing that heterogeneity in recombination rates have left a footprint on the genomic landscape of LD in zebra finch populations.
|
|
5. |
- Beyan, Huriya, et al.
(author)
-
Guthrie card methylomics identifies temporally stable epialleles that are present at birth in humans
- 2012
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051. ; 22:11, s. 2138-2145
-
Journal article (peer-reviewed)abstract
- A major concern in common disease epigenomics is distinguishing causal from consequential epigenetic variation. One means of addressing this issue is to identify the temporal origins of epigenetic variants via longitudinal analyses. However, prospective birth-cohort studies are expensive and time consuming. Here, we report DNA methylomics of archived Guthrie cards for the retrospective longitudinal analyses of in-utero-derived DNA methylation variation. We first validate two methodologies for generating comprehensive DNA methylomes from Guthrie cards. Then, using an integrated epigenomic/genomic analysis of Guthrie cards and follow-up samplings, we identify interindividual DNA methylation variation that is present both at birth and 3 yr later. These findings suggest that disease-relevant epigenetic variation could be detected at birth, i.e., before overt clinical disease. Guthrie card methylomics offers a potentially powerful and cost-effective strategy for studying the dynamics of interindividual epigenomic variation in a range of common human diseases.
|
|
6. |
- Denver, Dee R, et al.
(author)
-
Selective sweeps and parallel mutation in the adaptive recovery from deleterious mutation in Caenorhabditis elegans.
- 2010
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 20:12
-
Journal article (peer-reviewed)abstract
- Deleterious mutation poses a serious threat to human health and the persistence of small populations. Although adaptive recovery from deleterious mutation has been well-characterized in prokaryotes, the evolutionary mechanisms by which multicellular eukaryotes recover from deleterious mutation remain unknown. We applied high-throughput DNA sequencing to characterize genomic divergence patterns associated with the adaptive recovery from deleterious mutation using a Caenorhabditis elegans recovery-line system. The C. elegans recovery lines were initiated from a low-fitness mutation-accumulation (MA) line progenitor and allowed to independently evolve in large populations (N ∼ 1000) for 60 generations. All lines rapidly regained levels of fitness similar to the wild-type (N2) MA line progenitor. Although there was a near-zero probability of a single mutation fixing due to genetic drift during the recovery experiment, we observed 28 fixed mutations. Cross-generational analysis showed that all mutations went from undetectable population-level frequencies to a fixed state in 10-20 generations. Many recovery-line mutations fixed at identical timepoints, suggesting that the mutations, if not beneficial, hitchhiked to fixation during selective sweep events observed in the recovery lines. No MA line mutation reversions were detected. Parallel mutation fixation was observed for two sites in two independent recovery lines. Analysis using a C. elegans interactome map revealed many predicted interactions between genes with recovery line-specific mutations and genes with previously accumulated MA line mutations. Our study suggests that recovery-line mutations identified in both coding and noncoding genomic regions might have beneficial effects associated with compensatory epistatic interactions.
|
|
7. |
- Ekdahl, Ylva, et al.
(author)
-
A-to-I editing of microRNAs in the mammalian brain increases during development
- 2012
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 22:8, s. 1477-1487
-
Journal article (peer-reviewed)abstract
- Adenosine-to-inosine (A-to-I) RNA editing targets double-stranded RNA stem-loop structures in the mammalian brain. It has previously been shown that miRNAs are substrates for A-to-I editing. For the first time, we show that for several definitions of edited miRNA, the level of editing increases with development, thereby indicating a regulatory role for editing during brain maturation. We use high-throughput RNA sequencing to determine editing levels in mature miRNA, from the mouse transcriptome, and compare these with the levels of editing in pri-miRNA. We show that increased editing during development gradually changes the proportions of the two miR-376a isoforms, which previously have been shown to have different targets. Several other miRNAs that also are edited in the seed sequence show an increased level of editing through development. By comparing editing of pri-miRNA with editing and expression of the corresponding mature miRNA, we also show an editing-induced developmental regulation of miRNA expression. Taken together, our results imply that RNA editing influences the miRNA repertoire during brain maturation.
|
|
8. |
- Ellegren, Hans
(author)
-
Emergence of male-biased genes on the chicken Z-chromosome : Sex-chromosome contrasts between male and female heterogametic systems
- 2011
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 21:12, s. 2082-2086
-
Journal article (peer-reviewed)abstract
- There has been extensive traffic of male-biased genes out of the mammalian and Drosophila X-chromosomes, and there are also reports of an under-representation of male-biased genes on the X. This may reflect an adaptive process driven by natural selection where an autosomal location of male-biased genes is favored since male genes are only exposed to selection one-third of the time when X-linked. However, there are several alternative explanations to "out-of-the-X'' gene movement, including mutational bias and a means for X-linked genes to escape meiotic sex chromosome inactivation (MSCI) during spermatogenesis. As a critical test of the hypothesis that genomic relocation of sex-biased genes is an adaptive process, I examined the emergence, and loss, of genes on the chicken Z-chromosome, i.e., a female heterogametic system (males ZZ, females ZW). Here, the analogous prediction would be an emergence of male-biased genes onto, not a loss from, the Z-chromosome because Z is found more often in males than autosomes are. I found that genes expressed in testis but not in ovary are highly over-represented among genes that have emerged on the Z-chromosome during avian evolution. Moreover, genes with male-biased expression are similarly over-represented among new Z-chromosomal genes. Interestingly, genes with female-biased expression have more often moved from than to the Z-chromosome. These observations show that male and female heterogametic organisms display opposing directionalities in the emergence and loss of sex-biased genes on sex chromosomes. This is consistent with theoretical models on the evolution of sexually antagonistic genes in which new mutations are at least partly dominant.
|
|
9. |
- Hansen, KD, et al.
(author)
-
Large-scale hypomethylated blocks associated with Epstein-Barr virus-induced B-cell immortalization
- 2014
-
In: Genome research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051. ; 24:2, s. 177-184
-
Journal article (peer-reviewed)abstract
- Altered DNA methylation occurs ubiquitously in human cancer from the earliest measurable stages. A cogent approach to understanding the mechanism and timing of altered DNA methylation is to analyze it in the context of carcinogenesis by a defined agent. Epstein-Barr virus (EBV) is a human oncogenic herpesvirus associated with lymphoma and nasopharyngeal carcinoma, but also used commonly in the laboratory to immortalize human B-cells in culture. Here we have performed whole-genome bisulfite sequencing of normal B-cells, activated B-cells, and EBV-immortalized B-cells from the same three individuals, in order to identify the impact of transformation on the methylome. Surprisingly, large-scale hypomethylated blocks comprising two-thirds of the genome were induced by EBV immortalization but not by B-cell activation per se. These regions largely corresponded to hypomethylated blocks that we have observed in human cancer, and they were associated with gene-expression hypervariability, similar to human cancer, and consistent with a model of epigenomic change promoting tumor cell heterogeneity. We also describe small-scale changes in DNA methylation near CpG islands. These results suggest that methylation disruption is an early and critical step in malignant transformation.
|
|
10. |
- Hillmer, AM, et al.
(author)
-
Comprehensive long-span paired-end-tag mapping reveals characteristic patterns of structural variations in epithelial cancer genomes
- 2011
-
In: Genome research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051. ; 21:5, s. 665-675
-
Journal article (peer-reviewed)abstract
- Somatic genome rearrangements are thought to play important roles in cancer development. We optimized a long-span paired-end-tag (PET) sequencing approach using 10-Kb genomic DNA inserts to study human genome structural variations (SVs). The use of a 10-Kb insert size allows the identification of breakpoints within repetitive or homology-containing regions of a few kilobases in size and results in a higher physical coverage compared with small insert libraries with the same sequencing effort. We have applied this approach to comprehensively characterize the SVs of 15 cancer and two noncancer genomes and used a filtering approach to strongly enrich for somatic SVs in the cancer genomes. Our analyses revealed that most inversions, deletions, and insertions are germ-line SVs, whereas tandem duplications, unpaired inversions, interchromosomal translocations, and complex rearrangements are over-represented among somatic rearrangements in cancer genomes. We demonstrate that the quantitative and connective nature of DNA–PET data is precise in delineating the genealogy of complex rearrangement events, we observe signatures that are compatible with breakage-fusion-bridge cycles, and we discover that large duplications are among the initial rearrangements that trigger genome instability for extensive amplification in epithelial cancers.
|
|
11. |
- Inaki, K, et al.
(author)
-
Transcriptional consequences of genomic structural aberrations in breast cancer
- 2011
-
In: Genome research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051. ; 21:5, s. 676-687
-
Journal article (peer-reviewed)abstract
- Using a long-span, paired-end deep sequencing strategy, we have comprehensively identified cancer genome rearrangements in eight breast cancer genomes. Herein, we show that 40%–54% of these structural genomic rearrangements result in different forms of fusion transcripts and that 44% are potentially translated. We find that single segmental tandem duplication spanning several genes is a major source of the fusion gene transcripts in both cell lines and primary tumors involving adjacent genes placed in the reverse-order position by the duplication event. Certain other structural mutations, however, tend to attenuate gene expression. From these candidate gene fusions, we have found a fusion transcript (RPS6KB1–VMP1) recurrently expressed in ∼30% of breast cancers associated with potential clinical consequences. This gene fusion is caused by tandem duplication on 17q23 and appears to be an indicator of local genomic instability altering the expression of oncogenic components such as MIR21 and RPS6KB1.
|
|
12. |
- Islam, S, et al.
(author)
-
Characterization of the single-cell transcriptional landscape by highly multiplex RNA-seq
- 2011
-
In: Genome research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051. ; 21:7, s. 1160-1167
-
Journal article (peer-reviewed)abstract
- Our understanding of the development and maintenance of tissues has been greatly aided by large-scale gene expression analysis. However, tissues are invariably complex, and expression analysis of a tissue confounds the true expression patterns of its constituent cell types. Here we describe a novel strategy to access such complex samples. Single-cell RNA-seq expression profiles were generated, and clustered to form a two-dimensional cell map onto which expression data were projected. The resulting cell map integrates three levels of organization: the whole population of cells, the functionally distinct subpopulations it contains, and the single cells themselves—all without need for known markers to classify cell types. The feasibility of the strategy was demonstrated by analyzing the transcriptomes of 85 single cells of two distinct types. We believe this strategy will enable the unbiased discovery and analysis of naturally occurring cell types during development, adult physiology, and disease.
|
|
13. |
- Jolma, A, et al.
(author)
-
Multiplexed massively parallel SELEX for characterization of human transcription factor binding specificities
- 2010
-
In: Genome research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051. ; 20:6, s. 861-873
-
Journal article (peer-reviewed)abstract
- The genetic code—the binding specificity of all transfer-RNAs—defines how protein primary structure is determined by DNA sequence. DNA also dictates when and where proteins are expressed, and this information is encoded in a pattern of specific sequence motifs that are recognized by transcription factors. However, the DNA-binding specificity is only known for a small fraction of the ∼1400 human transcription factors (TFs). We describe here a high-throughput method for analyzing transcription factor binding specificity that is based on systematic evolution of ligands by exponential enrichment (SELEX) and massively parallel sequencing. The method is optimized for analysis of large numbers of TFs in parallel through the use of affinity-tagged proteins, barcoded selection oligonucleotides, and multiplexed sequencing. Data are analyzed by a new bioinformatic platform that uses the hundreds of thousands of sequencing reads obtained to control the quality of the experiments and to generate binding motifs for the TFs. The described technology allows higher throughput and identification of much longer binding profiles than current microarray-based methods. In addition, as our method is based on proteins expressed in mammalian cells, it can also be used to characterize DNA-binding preferences of full-length proteins or proteins requiring post-translational modifications. We validate the method by determining binding specificities of 14 different classes of TFs and by confirming the specificities for NFATC1 and RFX3 using ChIP-seq. Our results reveal unexpected dimeric modes of binding for several factors that were thought to preferentially bind DNA as monomers.
|
|
14. |
- Kristell, Carolina, et al.
(author)
-
Nitrogen depletion in the fission yeast Schizosaccharomyces pombe causes nucleosome loss in both promoters and coding regions of activated genes
- 2010
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 20:3, s. 361-371
-
Journal article (peer-reviewed)abstract
- Gene transcription is associated with local changes in chromatin, both in nucleosome positions and in chemical modifications of the histones. Chromatin dynamics has mostly been studied on a single-gene basis. Those genome-wide studies that have been made primarily investigated steady-state transcription. However, three studies of genome-wide changes in chromatin during the transcriptional response to heat shock in the budding yeast Saccharomyces cerevisiae revealed nucleosome eviction in promoter regions but only minor effects in coding regions. Here, we describe the short-term response to nitrogen starvation in the fission yeast Schizosaccharomyces pombe. Nitrogen depletion leads to a fast induction of a large number of genes in S. pombe and is thus suitable for genome-wide studies of chromatin dynamics during gene regulation. After 20 min of nitrogen removal, 118 transcripts were up-regulated. The distribution of regulated genes throughout the genome was not random; many up-regulated genes were found in clusters, while large parts of the genome were devoid of up-regulated genes. Surprisingly, this up-regulation was associated with nucleosome eviction of equal magnitudes in the promoters and in the coding regions. The nucleosome loss was not limited to induction by nitrogen depletion but also occurred during cadmium treatment. Furthermore, the lower nucleosome density persisted for at least 60 min after induction. Two highly induced genes, urg1(+) and urg2(+), displayed a substantial nucleosome loss, with only 20% of the nucleosomes being left in the coding region. We conclude that nucleosome loss during transcriptional activation is not necessarily limited to promoter regions.
|
|
15. |
- Meunier, Julien, et al.
(author)
-
Birth and expression evolution of mammalian microRNA genes
- 2013
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 23:1, s. 34-45
-
Journal article (peer-reviewed)abstract
- MicroRNAs (miRNAs) are major post-transcriptional regulators of gene expression, yet their origins and functional evolution in mammals remain little understood due to the lack of appropriate comparative data. Using RNA sequencing, we have generated extensive and comparable miRNA data for five organs in six species that represent all main mammalian lineages and birds (the evolutionary outgroup) with the aim to unravel the evolution of mammalian miRNAs. Our analyses reveal an overall expansion of miRNA repertoires in mammals, with threefold accelerated birth rates of miRNA families in placentals and marsupials, facilitated by the de novo emergence of miRNAs in host gene introns. Generally, our analyses suggest a high rate of miRNA family turnover in mammals with many newly emerged miRNA families being lost soon after their formation. Selectively preserved mammalian miRNA families gradually evolved higher expression levels, as well as altered mature sequences and target gene repertoires, and were apparently mainly recruited to exert regulatory functions in nervous tissues. However, miRNAs that originated on the X chromosome evolved high expression levels and potentially diverse functions during spermatogenesis, including meiosis, through selectively driven duplication-divergence processes. Overall, our study thus provides detailed insights into the birth and evolution of mammalian miRNA genes and the associated selective forces.
|
|
16. |
- Mondal, Tanmoy, et al.
(author)
-
Characterization of the RNA content of chromatin
- 2010
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 20:7, s. 899-907
-
Journal article (peer-reviewed)abstract
- Noncoding RNA (ncRNA) constitutes a significant portion of the mammalian transcriptome. Emerging evidence suggests that it regulates gene expression in cis or trans by modulating the chromatin structure. To uncover the functional role of ncRNA in chromatin organization, we deep sequenced chromatin-associated RNAs (CARs) from human fibroblast (HF) cells. This resulted in the identification of 141 intronic regions and 74 intergenic regions harboring CARs. The intronic and intergenic CARs show significant conservation across 44 species of placental mammals. Functional characterization of one of the intergenic CARs, Intergenic10, revealed that it regulates gene expression of neighboring genes through modulating the chromatin structure in cis. Our data suggest that ncRNA is an integral component of chromatin and that it may regulate various biological functions through fine-tuning of the chromatin architecture.
|
|
17. |
- Navin, N, et al.
(author)
-
Inferring tumor progression from genomic heterogeneity
- 2010
-
In: Genome research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051. ; 20:1, s. 68-80
-
Journal article (peer-reviewed)abstract
- Cancer progression in humans is difficult to infer because we do not routinely sample patients at multiple stages of their disease. However, heterogeneous breast tumors provide a unique opportunity to study human tumor progression because they still contain evidence of early and intermediate subpopulations in the form of the phylogenetic relationships. We have developed a method we call Sector-Ploidy-Profiling (SPP) to study the clonal composition of breast tumors. SPP involves macro-dissecting tumors, flow-sorting genomic subpopulations by DNA content, and profiling genomes using comparative genomic hybridization (CGH). Breast carcinomas display two classes of genomic structural variation: (1) monogenomic and (2) polygenomic. Monogenomic tumors appear to contain a single major clonal subpopulation with a highly stable chromosome structure. Polygenomic tumors contain multiple clonal tumor subpopulations, which may occupy the same sectors, or separate anatomic locations. In polygenomic tumors, we show that heterogeneity can be ascribed to a few clonal subpopulations, rather than a series of gradual intermediates. By comparing multiple subpopulations from different anatomic locations, we have inferred pathways of cancer progression and the organization of tumor growth.
|
|
18. |
- Nestor, Colm, 1977-, et al.
(author)
-
Tissue type is a major modifier of the 5-hydroxymethylcytosine content of human genes
- 2012
-
In: Genome Research. - : Cold Spring Harbor Laboratory Press (CSHL). - 1088-9051 .- 1549-5469. ; 22:3, s. 467-477
-
Journal article (peer-reviewed)abstract
- The discovery of substantial amounts of 5-hydroxymethylcytosine (5hmC), formed by the oxidation of 5-methylcytosine (5mC), in various mouse tissues and human embryonic stem (ES) cells has necessitated a reevaluation of our knowledge of 5mC/5hmC patterns and functions in mammalian cells. Here, we investigate the tissue specificity of both the global levels and locus-specific distribution of 5hmC in several human tissues and cell lines. We find that global 5hmC content of normal human tissues is highly variable, does not correlate with global 5mC content, and decreases rapidly as cells from normal tissue adapt to cell culture. Using tiling microarrays to map 5hmC levels in DNA from normal human tissues, we find that 5hmC patterns are tissue specific; unsupervised hierarchical clustering based solely on 5hmC patterns groups independent biological samples by tissue type. Moreover, in agreement with previous studies, we find 5hmC associated primarily, but not exclusively, with the body of transcribed genes, and that within these genes 5hmC levels are positively correlated with transcription levels. However, using quantitative 5hmC-qPCR, we find that the absolute levels of 5hmC for any given gene are primarily determined by tissue type, gene expression having a secondary influence on 5hmC levels. That is, a gene transcribed at a similar level in several different tissues may have vastly different levels of 5hmC (>20-fold) dependent on tissue type. Our findings highlight tissue type as a major modifier of 5hmC levels in expressed genes and emphasize the importance of using quantitative analyses in the study of 5hmC levels.
|
|
19. |
- Nissenbaum, J, et al.
(author)
-
Susceptibility to chronic pain following nerve injury is genetically affected by CACNG2
- 2010
-
In: Genome research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051. ; 20:9, s. 1180-1190
-
Journal article (peer-reviewed)abstract
- Chronic neuropathic pain is affected by specifics of the precipitating neural pathology, psychosocial factors, and by genetic predisposition. Little is known about the identity of predisposing genes. Using an integrative approach, we discovered that CACNG2 significantly affects susceptibility to chronic pain following nerve injury. CACNG2 encodes for stargazin, a protein intimately involved in the trafficking of glutamatergic AMPA receptors. The protein might also be a Ca2+ channel subunit. CACNG2 has previously been implicated in epilepsy. Initially, using two fine-mapping strategies in a mouse model (recombinant progeny testing [RPT] and recombinant inbred segregation test [RIST]), we mapped a pain-related quantitative trait locus (QTL) (Pain1) into a 4.2-Mb interval on chromosome 15. This interval includes 155 genes. Subsequently, bioinformatics and whole-genome microarray expression analysis were used to narrow the list of candidates and ultimately to pinpoint Cacng2 as a likely candidate. Analysis of stargazer mice, a Cacng2 hypomorphic mutant, provided electrophysiological and behavioral evidence for the gene's functional role in pain processing. Finally, we showed that human CACNG2 polymorphisms are associated with chronic pain in a cohort of cancer patients who underwent breast surgery. Our findings provide novel information on the genetic basis of neuropathic pain and new insights into pain physiology that may ultimately enable better treatments.
|
|
20. |
- Parts, Leopold, et al.
(author)
-
Revealing the genetic structure of a trait by sequencing a population under selection.
- 2011
-
In: Genome research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051. ; 21:7, s. 1131-8
-
Journal article (peer-reviewed)abstract
- One approach to understanding the genetic basis of traits is to study their pattern of inheritance among offspring of phenotypically different parents. Previously, such analysis has been limited by low mapping resolution, high labor costs, and large sample size requirements for detecting modest effects. Here, we present a novel approach to map trait loci using artificial selection. First, we generated populations of 10-100 million haploid and diploid segregants by crossing two budding yeast strains of different heat tolerance for up to 12 generations. We then subjected these large segregant pools to heat stress for up to 12 d, enriching for beneficial alleles. Finally, we sequenced total DNA from the pools before and during selection to measure the changes in parental allele frequency. We mapped 21 intervals with significant changes in genetic background in response to selection, which is several times more than found with traditional linkage methods. Nine of these regions contained two or fewer genes, yielding much higher resolution than previous genomic linkage studies. Multiple members of the RAS/cAMP signaling pathway were implicated, along with genes previously not annotated with heat stress response function. Surprisingly, at most selected loci, allele frequencies stopped changing before the end of the selection experiment, but alleles did not become fixed. Furthermore, we were able to detect the same set of trait loci in a population of diploid individuals with similar power and resolution, and observed primarily additive effects, similar to what is seen for complex trait genetics in other diploid organisms such as humans.
|
|
21. |
- Picelli, S, et al.
(author)
-
Tn5 transposase and tagmentation procedures for massively scaled sequencing projects
- 2014
-
In: Genome research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051. ; 24:12, s. 2033-2040
-
Journal article (peer-reviewed)abstract
- Massively parallel DNA sequencing of thousands of samples in a single machine-run is now possible, but the preparation of the individual sequencing libraries is expensive and time-consuming. Tagmentation-based library construction, using the Tn5 transposase, is efficient for generating sequencing libraries but currently relies on undisclosed reagents, which severely limits development of novel applications and the execution of large-scale projects. Here, we present simple and robust procedures for Tn5 transposase production and optimized reaction conditions for tagmentation-based sequencing library construction. We further show how molecular crowding agents both modulate library lengths and enable efficient tagmentation from subpicogram amounts of cDNA. The comparison of single-cell RNA-sequencing libraries generated using produced and commercial Tn5 demonstrated equal performances in terms of gene detection and library characteristics. Finally, because naked Tn5 can be annealed to any oligonucleotide of choice, for example, molecular barcodes in single-cell assays or methylated oligonucleotides for bisulfite sequencing, custom Tn5 production and tagmentation enable innovation in sequencing-based applications.
|
|
22. |
- Riddle, Nicole C, et al.
(author)
-
Plasticity in patterns of histone modifications and chromosomal proteins in Drosophila heterochromatin
- 2011
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 21:2, s. 147-163
-
Journal article (peer-reviewed)abstract
- Eukaryotic genomes are packaged in two basic forms, euchromatin and heterochromatin. We have examined the composition and organization of Drosophila melanogaster heterochromatin in different cell types using ChIP-array analysis of histone modifications and chromosomal proteins. As anticipated, the pericentric heterochromatin and chromosome 4 are on average enriched for the "silencing" marks H3K9me2, H3K9me3, HP1a, and SU(VAR)3-9, and are generally depleted for marks associated with active transcription. The locations of the euchromatin-heterochromatin borders identified by these marks are similar in animal tissues and most cell lines, although the amount of heterochromatin is variable in some cell lines. Combinatorial analysis of chromatin patterns reveals distinct profiles for euchromatin, pericentric heterochromatin, and the 4th chromosome. Both silent and active protein-coding genes in heterochromatin display complex patterns of chromosomal proteins and histone modifications; a majority of the active genes exhibit both "activation" marks (e.g., H3K4me3 and H3K36me3) and "silencing" marks (e.g., H3K9me2 and HP1a). The hallmark of active genes in heterochromatic domains appears to be a loss of H3K9 methylation at the transcription start site. We also observe complex epigenomic profiles of intergenic regions, repeated transposable element (TE) sequences, and genes in the heterochromatic extensions. An unexpectedly large fraction of sequences in the euchromatic chromosome arms exhibits a heterochromatic chromatin signature, which differs in size, position, and impact on gene expression among cell types. We conclude that patterns of heterochromatin/euchromatin packaging show greater complexity and plasticity than anticipated. This comprehensive analysis provides a foundation for future studies of gene activity and chromosomal functions that are influenced by or dependent upon heterochromatin.
|
|
23. |
- Rowe, Helen M., et al.
(author)
-
TRIM28 repression of retrotransposon-based enhancers is necessary to preserve transcriptional dynamics in embryonic stem cells
- 2013
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051. ; 23:3, s. 452-461
-
Journal article (peer-reviewed)abstract
- TRIM28 is critical for the silencing of endogenous retroviruses (ERVs) in embryonic stem (ES) cells. Here, we reveal that an essential impact of this process is the protection of cellular gene expression in early embryos from perturbation by cis-acting activators contained within these retroelements. In TRIM28-depleted ES cells, repressive chromatin marks at ERVs are replaced by histone modifications typical of active enhancers, stimulating transcription of nearby cellular genes, notably those harboring bivalent promoters. Correspondingly, ERV-derived sequences can repress or enhance expression from an adjacent promoter in transgenic embryos depending on their TRIM28 sensitivity in ES cells. TRIM28-mediated control of ERVs is therefore crucial not just to prevent retrotransposition, but more broadly to safeguard the transcriptional dynamics of early embryos.
|
|
24. |
- Schwartz, Yuri B, et al.
(author)
-
Nature and function of insulator protein binding sites in the Drosophila genome
- 2012
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 22, s. 2188-2198
-
Journal article (peer-reviewed)abstract
- Chromatin insulator elements and associated proteins have been proposed to partition eukaryotic genomes into sets of independently regulated domains. Here we test this hypothesis by quantitative genome-wide analysis of insulator protein binding to Drosophila chromatin. We find distinct combinatorial binding of insulator proteins to different classes of sites and uncover a novel type of insulator element that binds CP190 but not any other known insulator proteins. Functional characterization of different classes of binding sites indicates that only a small fraction act as robust insulators in standard enhancer-blocking assays. We show that insulators restrict the spreading of the H3K27me3 mark but only at a small number of Polycomb target regions and only to prevent repressive histone methylation within adjacent genes that are already transcriptionally inactive. RNAi knockdown of insulator proteins in cultured cells does not lead to major alterations in genome expression. Taken together these observations argue against the concept of a genome partitioned by specialized boundary elements and suggest that insulators are reserved for specific regulation of selected genes.
|
|
25. |
- Seth-Smith, Helena M. B., et al.
(author)
-
Whole-genome sequences of Chlamydia trachomatis directly from clinical samples without culture
- 2013
-
In: Genome Research. - : Cold Spring Harbor Laboratory Press (CSHL). - 1088-9051 .- 1549-5469. ; 23:5, s. 855-866
-
Journal article (peer-reviewed)abstract
- The use of whole-genome sequencing as a tool for the study of infectious bacteria is of growing clinical interest. Chlamydia trachomatis is responsible for sexually transmitted infections and the blinding disease trachoma, which affect hundreds of millions of people worldwide. Recombination is widespread within the genome of C. trachomatis, thus whole-genome sequencing is necessary to understand the evolution, diversity, and epidemiology of this pathogen. Culture of C trachomatis has, until now, been a prerequisite to obtain DNA for whole-genome sequencing; however, as C trachomatis is an obligate intracellular pathogen, this procedure is technically demanding and time consuming. Discarded clinical samples represent a large resource for sequencing the genomes of pathogens, yet clinical swabs frequently contain very low levels of C trachomatis DNA and large amounts of contaminating microbial and human DNA. To determine whether it is possible to obtain whole-genome sequences from bacteria without the need for culture, we have devised an approach that combines immunomagnetic separation (IMS) for targeted bacterial enrichment with multiple displacement amplification (MDA) for whole-genome amplification. Using IMS-MDA ill conjunction with high-throughput multiplexed Illumina sequencing, we have produced the first whole bacterial genome sequences direct from clinical samples. We also show that this method can be used to generate genome data from nonviable archived samples. This method will prove a useful tool in answering questions relating to the biology of many difficult-to-culture or fastidious bacteria of clinical concern.
|
|
26. |
- Song, Lingyun, et al.
(author)
-
Open chromatin defined by DNaseI and FAIRE identifies regulatory elements that shape cell-type identity
- 2011
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 21:10, s. 1757-1767
-
Journal article (peer-reviewed)abstract
- The human body contains thousands of unique cell types, each with specialized functions. Cell identity is governed in large part by gene transcription programs, which are determined by regulatory elements encoded in DNA. To identify regulatory elements active in seven cell lines representative of diverse human cell types, we used DNase-seq and FAIRE-seq (Formaldehyde Assisted Isolation of Regulatory Elements) to map open chromatin.'' Over 870,000 DNaseI or FAIRE sites, which correspond tightly to nucleosome-depleted regions, were identified across the seven cell lines, covering nearly 9% of the genome. The combination of DNaseI and FAIRE is more effective than either assay alone in identifying likely regulatory elements, as judged by coincidence with transcription factor binding locations determined in the same cells. Open chromatin common to all seven cell types tended to be at or near transcription start sites and to be coincident with CTCF binding sites, while open chromatin sites found in only one cell type were typically located away from transcription start sites and contained DNA motifs recognized by regulators of cell-type identity. We show that open chromatin regions bound by CTCF are potent insulators. We identified clusters of open regulatory elements (COREs) that were physically near each other and whose appearance was coordinated among one or more cell types. Gene expression and RNA Pol II binding data support the hypothesis that COREs control gene activity required for the maintenance of cell-type identity. This publicly available atlas of regulatory elements may prove valuable in identifying noncoding DNA sequence variants that are causally linked to human disease.
|
|
27. |
- Vihinen, Mauno
(author)
-
Variation Ontology for annotation of variation effects and mechanisms.
- 2013
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051.
-
Journal article (peer-reviewed)abstract
- Ontology organizes and formally conceptualizes information in a knowledge domain with a controlled vocabulary having defined terms and relationships between them. Several ontologies have been used to annotate numerous databases in biology and medicine. Due to their unambiguous nature, ontological annotations facilitate systematic description and data organization, data integration and mining, pattern recognition and statistics, as well as development of analysis and prediction tools. The Variation Ontology was developed to allow the annotation of effects, consequences and mechanisms of DNA, RNA and protein variations. Variation types are systematically organized and a detailed description of effects and mechanisms is possible. VariO is for annotating the variant, not the normal state features or properties, and requires a reference (e.g. reference sequence, reference state property, activity etc) compared to which the changes are indicated. VariO is versatile and can be used for variations ranging from genomic multiplications to single nucleotide or amino acid changes whether of genetic or non-genetic origin. VariO annotations are position specific and can be used for variations in any organism.
|
|
28. |
- Völker, Martin, et al.
(author)
-
Copy number variation, chromosome rearrangement, and their association with recombination during avian evolution
- 2010
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 20:4, s. 503-511
-
Journal article (peer-reviewed)abstract
- Chromosomal rearrangements and copy number variants (CNVs) play key roles in genome evolution and genetic disease; however, the molecular mechanisms underlying these types of structural genomic variation are not fully understood. The availability of complete genome sequences for two bird species, the chicken and the zebra finch, provides, for the first time, an ideal opportunity to analyze the relationship between structural genomic variation (chromosomal and CNV) and recombination on a genome-wide level. The aims of this study were therefore threefold: (1) to combine bioinformatics, physical mapping to produce comprehensive comparative maps of the genomes of chicken and zebra finch. In so doing, this allowed the identification of evolutionary chromosomal rearrangements distinguishing them. The previously reported interchromosomal conservation of synteny was confirmed, but a larger than expected number of intrachromosomal rearrangements were reported; (2) to hybridize zebra finch genomic DNA to a chicken tiling path microarray and identify CNVs in the zebra finch genome relative to chicken; 32 interspecific CNVs were identified; and (3) to test the hypothesis that there is an association between CNV, chromosomal rearrangements, and recombination by correlating data from (1) and (2) with recombination rate data from a high-resolution genetic linkage map of the zebra finch. We found a highly significant association of both chromosomal rearrangements and CNVs with elevated recombination rates. The results thus provide support for the notion of recombination-based processes playing a major role in avian genome evolution.
|
|
29. |
- Laabei, M., et al.
(author)
-
Predicting the virulence of MRSA from its genome sequence
- 2014
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051. ; 24:5, s. 839-849
-
Journal article (peer-reviewed)abstract
- Microbial virulence is a complex and often multifactorial phenotype, intricately linked to a pathogen's evolutionary trajectory. Toxicity, the ability to destroy host cell membranes, and adhesion, the ability to adhere to human tissues, are the major virulence factors of many bacterial pathogens, including Staphylococcus aureus. Here, we assayed the toxicity and adhesiveness of 90 MRSA (methicillin resistant S. aureus) isolates and found that while there was remarkably little variation in adhesion, toxicity varied by over an order of magnitude between isolates, suggesting different evolutionary selection pressures acting on these two traits. We performed a genome-wide association study (GWAS) and identified a large number of loci, as well as a putative network of epistatically interacting loci, that significantly associated with toxicity. Despite this apparent complexity in toxicity regulation, a predictive model based on a set of significant single nucleotide polymorphisms (SNPs) and insertion and deletions events (indels) showed a high degree of accuracy in predicting an isolate's toxicity solely from the genetic signature at these sites. Our results thus highlight the potential of using sequence data to determine clinically relevant parameters and have further implications for understanding the microbial virulence of this opportunistic pathogen.
|
|