| 1. |
- Burri, Reto, et al.
(author)
-
Linked selection and recombination rate variation drive the evolution of the genomic landscape of differentiation across the speciation continuum of Ficedula flycatchers
- 2015
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 25:11, s. 1656-1665
-
Journal article (peer-reviewed)abstract
- Speciation is a continuous process during which genetic changes gradually accumulate in the genomes of diverging species. Recent studies have documented highly heterogeneous differentiation landscapes, with distinct regions of elevated differentiation ("differentiation islands") widespread across genomes. However, it remains unclear which processes drive the evolution of differentiation islands; how the differentiation landscape evolves as speciation advances; and ultimately, how differentiation islands are related to speciation. Here, we addressed these questions based on population genetic analyses of 200 resequenced genomes from 10 populations of four Ficedula flycatcher sister species. We show that a heterogeneous differentiation landscape starts emerging among populations within species, and differentiation islands evolve recurrently in the very same genomic regions among independent lineages. Contrary to expectations from models that interpret differentiation islands as genomic regions involved in reproductive isolation that are shielded from gene flow, patterns of sequence divergence (d(XY) relative node depth) do not support a major role of gene flow in the evolution of the differentiation landscape in these species. Instead, as predicted by models of linked selection, genome-wide variation in diversity and differentiation can be explained by variation in recombination rate and the density of targets for selection. We thus conclude that the heterogeneous landscape of differentiation in Ficedula flycatchers evolves mainly as the result of background selection and selective sweeps in genomic regions of low recombination. Our results emphasize the necessity of incorporating linked selection as a null model to identify genome regions involved in adaptation and speciation.
|
|
| 2. |
- Böhme, Ulrike, et al.
(author)
-
Complete avian malaria parasite genomes reveal features associated with lineage-specific evolution in birds and mammals
- 2018
-
In: Genome Research. - : Cold Spring Harbor Laboratory Press (CSHL). - 1088-9051 .- 1549-5469. ; 28:4, s. 547-560
-
Journal article (peer-reviewed)abstract
- Avian malaria parasites are prevalent around the world and infect a wide diversity of bird species. Here, we report the sequencing and analysis of high-quality draft genome sequences for two avian malaria species, Plasmodium relictum and Plasmodium gallinaceum. We identify 50 genes that are specific to avian malaria, located in an otherwise conserved core of the genome that shares gene synteny with all other sequenced malaria genomes. Phylogenetic analysis suggests that the avian malaria species form an outgroup to the mammalian Plasmodium species, and using amino acid divergence between species, we estimate the avian- and mammalian-infective lineages diverged in the order of 10 million years ago. Consistent with their phylogenetic position, we identify orthologs of genes that had previously appeared to be restricted to the clades of parasites containing Plasmodium falciparum and Plasmodium vivax, the species with the greatest impact on human health. From these orthologs, we explore differential diversifying selection across the genus and show that the avian lineage is remarkable in the extent to which invasion-related genes are evolving. The subtelomeres of the P. relictum and P. gallinaceum genomes contain several novel gene families, including an expanded surf multigene family. We also identify an expansion of reticulocyte binding protein homologs in P. relictum, and within these proteins, we detect distinct regions that are specific to nonhuman primate, humans, rodent, and avian hosts. For the first time in the Plasmodium lineage, we find evidence of transposable elements, including several hundred fragments of LTR-retrotransposons in both species and an apparently complete LTR-retrotransposon in the genome of P. gallinaceum.
|
|
| 3. |
- Chen, G, et al.
(author)
-
Single-cell analyses of X Chromosome inactivation dynamics and pluripotency during differentiation
- 2016
-
In: Genome research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051. ; 26:10, s. 1342-1354
-
Journal article (peer-reviewed)abstract
- Pluripotency, differentiation, and X Chromosome inactivation (XCI) are key aspects of embryonic development. However, the underlying relationship and mechanisms among these processes remain unclear. Here, we systematically dissected these features along developmental progression using mouse embryonic stem cells (mESCs) and single-cell RNA sequencing with allelic resolution. We found that mESCs grown in a ground state 2i condition displayed transcriptomic profiles diffused from preimplantation mouse embryonic cells, whereas EpiStem cells closely resembled the post-implantation epiblast. Sex-related gene expression varied greatly across distinct developmental states. We also identified novel markers that were highly enriched in each developmental state. Moreover, we revealed that several novel pathways, including PluriNetWork and Focal Adhesion, were responsible for the delayed progression of female EpiStem cells. Importantly, we “digitalized” XCI progression using allelic expression of active and inactive X Chromosomes and surprisingly found that XCI states exhibited profound variability in each developmental state, including the 2i condition. XCI progression was not tightly synchronized with loss of pluripotency and increase of differentiation at the single-cell level, although these processes were globally correlated. In addition, highly expressed genes, including core pluripotency factors, were in general biallelically expressed. Taken together, our study sheds light on the dynamics of XCI progression and the asynchronicity between pluripotency, differentiation, and XCI.
|
|
| 4. |
- Chen, Jenny, et al.
(author)
-
A quantitative framework for characterizing the evolutionary history of mammalian gene expression
- 2019
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 29:1, s. 53-63
-
Journal article (peer-reviewed)abstract
- The evolutionary history of a gene helps predict its function and relationship to phenotypic traits. Although sequence conservation is commonly used to decipher gene function and assess medical relevance, methods for functional inference from comparative expression data are lacking. Here, we use RNA-seq across seven tissues from 17 mammalian species to show that expression evolution across mammals is accurately modeled by the Ornstein-Uhlenbeck process, a commonly proposed model of continuous trait evolution. We apply this model to identify expression pathways under neutral, stabilizing, and directional selection. We further demonstrate novel applications of this model to quantify the extent of stabilizing selection on a gene's expression, parameterize the distribution of each gene's optimal expression level, and detect deleterious expression levels in expression data from individual patients. Our work provides a statistical framework for interpreting expression data across species and in disease.
|
|
| 5. |
- Corcoran, Padraic, et al.
(author)
-
Introgression maintains the genetic integrity of the mating-type determining chromosome of the fungus Neurospora tetrasperma.
- 2016
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 26:4, s. 486-498
-
Journal article (peer-reviewed)abstract
- Genome evolution is driven by a complex interplay of factors, including selection, recombination, and introgression. The regions determining sexual identity are particularly dynamic parts of eukaryotic genomes that are prone to molecular degeneration associated with suppressed recombination. In the fungus Neurospora tetrasperma, it has been proposed that this molecular degeneration is counteracted by the introgression of nondegenerated DNA from closely related species. In this study, we used comparative and population genomic analyses of 92 genomes from eight phylogenetically and reproductively isolated lineages of N. tetrasperma, and its three closest relatives, to investigate the factors shaping the evolutionary history of the genomes. We found that suppressed recombination extends across at least 6 Mbp (similar to 63%) of the mating-type (mat) chromosome in N. tetrasperma and is associated with decreased genetic diversity, which is likely the result primarily of selection at linked sites. Furthermore, analyses of molecular evolution revealed an increased mutational load in this region, relative to recombining regions. However, comparative genomic and phylogenetic analyses indicate that the mat chromosomes are temporarily regenerated via introgression from sister species; six of eight lineages show introgression into one of their mat chromosomes, with multiple Neurospora species acting as donors. The introgressed tracts have been fixed within lineages, suggesting that they confer an adaptive advantage in natural populations, and our analyses support the presence of selective sweeps in at least one lineage. Thus, these data strongly support the previously hypothesized role of introgression as a mechanism for the maintenance of mating-type determining chromosomal regions.
|
|
| 6. |
- Edmonson, MN, et al.
(author)
-
Pediatric Cancer Variant Pathogenicity Information Exchange (PeCanPIE): a cloud-based platform for curating and classifying germline variants
- 2019
-
In: Genome research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051. ; 29:9, s. 1555-1565
-
Journal article (peer-reviewed)abstract
- Variant interpretation in the era of massively parallel sequencing is challenging. Although many resources and guidelines are available to assist with this task, few integrated end-to-end tools exist. Here, we present the Pediatric Cancer Variant Pathogenicity Information Exchange (PeCanPIE), a web- and cloud-based platform for annotation, identification, and classification of variations in known or putative disease genes. Starting from a set of variants in variant call format (VCF), variants are annotated, ranked by putative pathogenicity, and presented for formal classification using a decision-support interface based on published guidelines from the American College of Medical Genetics and Genomics (ACMG). The system can accept files containing millions of variants and handle single-nucleotide variants (SNVs), simple insertions/deletions (indels), multiple-nucleotide variants (MNVs), and complex substitutions. PeCanPIE has been applied to classify variant pathogenicity in cancer predisposition genes in two large-scale investigations involving >4000 pediatric cancer patients and serves as a repository for the expert-reviewed results. PeCanPIE was originally developed for pediatric cancer but can be easily extended for use for nonpediatric cancers and noncancer genetic diseases. Although PeCanPIE's web-based interface was designed to be accessible to non-bioinformaticians, its back-end pipelines may also be run independently on the cloud, facilitating direct integration and broader adoption. PeCanPIE is publicly available and free for research use.
|
|
| 7. |
- Elvers, Ingegerd, et al.
(author)
-
Exome sequencing of lymphomas from three dog breeds reveals somatic mutation patterns reflecting genetic background
- 2015
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 25:11, s. 1634-1645
-
Journal article (peer-reviewed)abstract
- Lymphoma is the most common hematological malignancy in developed countries. Outcome is strongly determined by molecular subtype, reflecting a need for new and improved treatment options. Dogs spontaneously develop lymphoma, and the predisposition of certain breeds indicates genetic risk factors. Using the dog breed structure, we selected three lymphoma predisposed breeds developing primarily T-cell (boxer), primarily B-cell (cocker spaniel), and with equal distribution of B- and T-cell lymphoma (golden retriever), respectively. We investigated the somatic mutations in B- and T-cell lymphomas from these breeds by exome sequencing of tumor and normal pairs. Strong similarities were evident between B-cell lymphomas from golden retrievers and cocker spaniels, with recurrent mutations in TRAF3-MAP3K14 (28% of all cases), FBXW7 (25%), and POT1 (17%). The FBXW7 mutations recurrently occur in a specific codon; the corresponding codon is recurrently mutated in human cancer. In contrast, T-cell lymphomas from the predisposed breeds, boxers and golden retrievers, show little overlap in their mutation pattern, sharing only one of their 15 most recurrently mutated genes. Boxers, which develop aggressive T-cell lymphomas, are typically mutated in the PTEN-mTOR pathway. T-cell lymphomas in golden retrievers are often less aggressive, and their tumors typically showed mutations in genes involved in cellular metabolism. We identify genes with known involvement in human lymphoma and leukemia, genes implicated in other human cancers, as well as novel genes that could allow new therapeutic options.
|
|
| 8. |
- Fabbri, M, et al.
(author)
-
Decrypting noncoding RNA interactions, structures, and functional networks
- 2019
-
In: Genome research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051. ; 29:9, s. 1377-1388
-
Journal article (peer-reviewed)abstract
- The world of noncoding RNAs (ncRNAs) is composed of an enormous and growing number of transcripts, ranging in length from tens of bases to tens of kilobases, involved in all biological processes and altered in expression and/or function in many types of human disorders. The premise of this review is the concept that ncRNAs, like many large proteins, have a multidomain architecture that organizes them spatially and functionally. As ncRNAs are beginning to be imprecisely classified into functional families, we review here how their structural properties might inform their functions with focus on structural architecture–function relationships. We will describe the properties of “interactor elements” (IEs) involved in direct physical interaction with nucleic acids, proteins, or lipids and of “structural elements” (SEs) directing their wiring within the “ncRNA interactor networks” through the emergence of secondary and/or tertiary structures. We suggest that spectrums of “letters” (ncRNA elements) are assembled into “words” (ncRNA domains) that are further organized into “phrases” (complete ncRNA structures) with functional meaning (signaling output) through complex “sentences” (the ncRNA interactor networks). This semiotic analogy can guide the exploitation of ncRNAs as new therapeutic targets through the development of IE-blockers and/or SE-lockers that will change the interactor partners’ spectrum of proteins, RNAs, DNAs, or lipids and consequently influence disease phenotypes.
|
|
| 9. |
- Forsberg, Lars A., et al.
(author)
-
Signatures of post-zygotic structural genetic aberrations in the cells of histologically normal breast tissue that can predispose to sporadic breast cancer
- 2015
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 25:10, s. 1521-1535
-
Journal article (peer-reviewed)abstract
- Sporadic breast cancer (SBC) is a common disease without robust means of early risk prediction in the population. We studied 282 females with SBC, focusing on copy number aberrations in cancer-free breast tissue (uninvolved margin, UM) outside the primary tumor (PT). In total, 1162 UMs (1-14 per breast) were studied. Comparative analysis between UM(s), PT(s), and blood/skin from the same patient as a control is the core of the study design. We identified 108 patients with at least one aberrant UM, representing 38.3% of cases. Gains in gene copy number were the principal type of mutations in microscopically normal breast cells, suggesting that oncogenic activation of genes via increased gene copy number is a predominant mechanism for initiation of SBC pathogenesis. The gain of ERBB2, with overexpression of HER2 protein, was the most common aberration in normal cells. Five additional growth factor receptor genes (EGFR, FGFR1, IGF1R, LIFR, and NGFR) also showed recurrent gains, and these were occasionally present in combination with the gain of ERBB2. All the aberrations found in the normal breast cells were previously described in cancer literature, suggesting their causative, driving role in pathogenesis of SBC. We demonstrate that analysis of normal cells from cancer patients leads to identification of signatures that may increase risk of SBC and our results could influence the choice of surgical intervention to remove all predisposing cells. Early detection of copy number gains suggesting a predisposition toward cancer development, long before detectable tumors are formed, is a key to the anticipated shift into a preventive paradigm of personalized medicine for breast cancer.
|
|
| 10. |
- Gandin, V, et al.
(author)
-
nanoCAGE reveals 5' UTR features that define specific modes of translation of functionally related MTOR-sensitive mRNAs
- 2016
-
In: Genome research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051. ; 26:5, s. 636-648
-
Journal article (peer-reviewed)abstract
- The diversity of MTOR-regulated mRNA translation remains unresolved. Whereas ribosome-profiling suggested that MTOR almost exclusively stimulates translation of the TOP (terminal oligopyrimidine motif) and TOP-like mRNAs, polysome-profiling indicated that MTOR also modulates translation of mRNAs without the 5′ TOP motif (non-TOP mRNAs). We demonstrate that in ribosome-profiling studies, detection of MTOR-dependent changes in non-TOP mRNA translation was obscured by low sensitivity and methodology biases. Transcription start site profiling using nano-cap analysis of gene expression (nanoCAGE) revealed that not only do many MTOR-sensitive mRNAs lack the 5′ TOP motif but that 5′ UTR features distinguish two functionally and translationally distinct subsets of MTOR-sensitive mRNAs: (1) mRNAs with short 5′ UTRs enriched for mitochondrial functions, which require EIF4E but are less EIF4A1-sensitive; and (2) long 5′ UTR mRNAs encoding proliferation- and survival-promoting proteins, which are both EIF4E- and EIF4A1-sensitive. Selective inhibition of translation of mRNAs harboring long 5′ UTRs via EIF4A1 suppression leads to sustained expression of proteins involved in respiration but concomitant loss of those protecting mitochondrial structural integrity, resulting in apoptosis. Conversely, simultaneous suppression of translation of both long and short 5′ UTR mRNAs by MTOR inhibitors results in metabolic dormancy and a predominantly cytostatic effect. Thus, 5′ UTR features define different modes of MTOR-sensitive translation of functionally distinct subsets of mRNAs, which may explain the diverse impact of MTOR and EIF4A inhibitors on neoplastic cells.
|
|
| 11. |
- Glemin, Sylvain, et al.
(author)
-
Quantification of GC-biased gene conversion in the human genome
- 2015
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 25:8, s. 1215-1228
-
Journal article (peer-reviewed)abstract
- Much evidence indicates that GC-biased gene conversion (gBGC) has a major impact on the evolution of mammalian genomes. However, a detailed quantification of the process is still lacking. The strength of gBGC can be measured from the analysis of derived allele frequency spectra (DAF), but this approach is sensitive to a number of confounding factors. In particular, we show by simulations that the inference is pervasively affected by polymorphism polarization errors and by spatial heterogeneity in gBGC strength. We propose a new general method to quantify gBGC from DAF spectra, incorporating polarization errors, taking spatial heterogeneity into account, and jointly estimating mutation bias. Applying it to human polymorphism data from the 1000 Genomes Project, we show that the strength of gBGC does not differ between hypermutable CpG sites and non-CpG sites, suggesting that in humans gBGC is not caused by the base-excision repair machinery. Genome-wide, the intensity of gBGC is in the nearly neutral area. However, given that recombination occurs primarily within recombination hotspots, 1%-2% of the human genome is subject to strong gBGC. On average, gBGC is stronger in African than in non-African populations, reflecting differences in effective population sizes. However, due to more heterogeneous recombination landscapes, the fraction of the genome affected by strong gBGC is larger in non-African than in African populations. Given that the location of recombination hotspots evolves very rapidly, our analysis predicts that, in the long term, a large fraction of the genome is affected by short episodes of strong gBGC.
|
|
| 12. |
- Guschanski, Katerina, et al.
(author)
-
The evolution of duplicate gene expression in mammalian organs
- 2017
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 27:9, s. 1461-1474
-
Journal article (peer-reviewed)abstract
- Gene duplications generate genomic raw material that allows the emergence of novel functions, likely facilitating adaptive evolutionary innovations. However, global assessments of the functional and evolutionary relevance of duplicate genes in mammals were until recently limited by the lack of appropriate comparative data. Here, we report a large-scale study of the expression evolution of DNA-based functional gene duplicates in three major mammalian lineages (placental mammals, marsupials, egg-laying monotremes) and birds, on the basis of RNA sequencing (RNA-seq) data from nine species and eight organs. We observe dynamic changes in tissue expression preference of paralogs with different duplication ages, suggesting differential contribution of paralogs to specific organ functions during vertebrate evolution. Specifically, we show that paralogs that emerged in the common ancestor of bony vertebrates are enriched for genes with brain-specific expression and provide evidence for differential forces underlying the preferential emergence of young testis-and liver-specific expressed genes. Further analyses uncovered that the overall spatial expression profiles of gene families tend to be conserved, with several exceptions of pronounced tissue specificity shifts among lineage-specific gene family expansions. Finally, we trace new lineage-specific genes that may have contributed to the specific biology of mammalian organs, including the little-studied placenta. Overall, our study provides novel and taxonomically broad evidence for the differential contribution of duplicate genes to tissue-specific transcriptomes and for their importance for the phenotypic evolution of vertebrates.
|
|
| 13. |
- Hadfield, James, et al.
(author)
-
Comprehensive global genome dynamics of Chlamydia trachomatis show ancient diversification followed by contemporary mixing and recent lineage expansion
- 2017
-
In: Genome Research. - : Cold Spring Harbor Laboratory Press. - 1088-9051 .- 1549-5469. ; 27:7, s. 1220-1229
-
Journal article (peer-reviewed)abstract
- Chlamydia trachomatis is the world's most prevalent bacterial sexually transmitted infection and leading infectious cause of blindness, yet it is one of the least understood human pathogens, in part due to the difficulties of in vitro culturing and the lack of available tools for genetic manipulation. Genome sequencing has reinvigorated this field, shedding light on the contemporary history of this pathogen. Here, we analyze 563 full genomes, 455 of which are novel, to show that the history of the species comprises two phases, and conclude that the currently circulating lineages are the result of evolution in different genomic ecotypes. Temporal analysis indicates these lineages have recently expanded in the space of thousands of years, rather than the millions of years as previously thought, a finding that dramatically changes our understanding of this pathogen's history. Finally, at a time when almost every pathogen is becoming increasingly resistant to antimicrobials, we show that there is no evidence of circulating genomic resistance in C. trachomatis.
|
|
| 14. |
- Hagey, Daniel W., et al.
(author)
-
Distinct transcription factor complexes act on a permissive chromatin landscape to establish regionalized gene expression in CNS stem cells
- 2016
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 26:7, s. 908-917
-
Journal article (peer-reviewed)abstract
- Spatially distinct gene expression profiles in neural stem cells (NSCs) are a prerequisite to the formation of neuronal diversity, but how these arise from the regulatory interactions between chromatin accessibility and transcription factor activity has remained unclear. Here, we demonstrate that, despite their distinct gene expression profiles, NSCs of the mouse cortex and spinal cord share the majority of their DNase I hypersensitive sites (DHSs). Regardless of this similarity, domain-specific gene expression is highly correlated with the relative accessibility of associated DHSs, as determined by sequence read density. Notably, the binding pattern of the general NSC transcription factor SOX2 is also largely cell type specific and coincides with an enrichment of LHX2 motifs in the cortex and HOXA9 motifs in the spinal cord. Interestingly, in a zebrafish reporter gene system, these motifs were critical determinants of patterned gene expression along the rostral-caudal axis. Our findings establish a predictive model for patterned NSC gene expression, whereby domain-specific expression of LHX2 and HOX proteins act on their target motifs within commonly accessible cis-regulatory regions to specify SOX2 binding. In turn, this binding correlates strongly with these DHSs relative accessibility-a robust predictor of neighboring gene expression.
|
|
| 15. |
- Han, Fan, et al.
(author)
-
Gene flow, ancient polymorphism, and ecological adaptation shape the genomic landscape of divergence among Darwin's finches
- 2017
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 27:6, s. 1004-1015
-
Journal article (peer-reviewed)abstract
- Genomic comparisons of closely related species have identified "islands" of locally elevated sequence divergence. Genomic islands may contain functional variants involved in local adaptation or reproductive isolation and may therefore play an important role in the speciation process. However, genomic islands can also arise through evolutionary processes unrelated to speciation, and examination of their properties can illuminate how new species evolve. Here, we performed scans for regions of high relative divergence (FST) in 12 species pairs of Darwin's finches at different genetic distances. In each pair, we identify genomic islands that are, on average, elevated in both relative divergence (FST) and absolute divergence (dXY). This signal indicates that haplotypes within these genomic regions became isolated from each other earlier than the rest of the genome. Interestingly, similar numbers of genomic islands of elevated dXY are observed in sympatric and allopatric species pairs, suggesting that recent gene flow is not a major factor in their formation. We find that two of the most pronounced genomic islands contain the ALX1 and HMGA2 loci, which are associated with variation in beak shape and size, respectively, suggesting that they are involved in ecological adaptation. A subset of genomic island regions, including these loci, appears to represent anciently diverged haplotypes that evolved early during the radiation of Darwin's finches. Comparative genomics data indicate that these loci, and genomic islands in general, have exceptionally low recombination rates, which may play a role in their establishment.
|
|
| 16. |
- Hashimoto, K, et al.
(author)
-
CAGE profiling of ncRNAs in hepatocellular carcinoma reveals widespread activation of retroviral LTR promoters in virus-induced tumors
- 2015
-
In: Genome research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051. ; 25:12, s. 1812-1824
-
Journal article (peer-reviewed)abstract
- An increasing number of noncoding RNAs (ncRNAs) have been implicated in various human diseases including cancer; however, the ncRNA transcriptome of hepatocellular carcinoma (HCC) is largely unexplored. We used CAGE to map transcription start sites across various types of human and mouse HCCs with emphasis on ncRNAs distant from protein-coding genes. Here, we report that retroviral LTR promoters, expressed in healthy tissues such as testis and placenta but not liver, are widely activated in liver tumors. Despite HCC heterogeneity, a subset of LTR-derived ncRNAs were more than 10-fold up-regulated in the vast majority of samples. HCCs with a high LTR activity mostly had a viral etiology, were less differentiated, and showed higher risk of recurrence. ChIP-seq data show that MYC and MAX are associated with ncRNA deregulation. Globally, CAGE enabled us to build a mammalian promoter map for HCC, which uncovers a new layer of complexity in HCC genomics.
|
|
| 17. |
- Ju, Young Seok, et al.
(author)
-
Frequent somatic transfer of mitochondrial DNA into the nuclear genome of human cancer cells.
- 2015
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051. ; 25:6, s. 814-824
-
Journal article (peer-reviewed)abstract
- Mitochondrial genomes are separated from the nuclear genome for most of the cell cycle by the nuclear double membrane, intervening cytoplasm, and the mitochondrial double membrane. Despite these physical barriers, we show that somatically acquired mitochondrial-nuclear genome fusion sequences are present in cancer cells. Most occur in conjunction with intranuclear genomic rearrangements, and the features of the fusion fragments indicate that nonhomologous end joining and/or replication-dependent DNA double-strand break repair are the dominant mechanisms involved. Remarkably, mitochondrial-nuclear genome fusions occur at a similar rate per base pair of DNA as interchromosomal nuclear rearrangements, indicating the presence of a high frequency of contact between mitochondrial and nuclear DNA in some somatic cells. Transmission of mitochondrial DNA to the nuclear genome occurs in neoplastically transformed cells, but we do not exclude the possibility that some mitochondrial-nuclear DNA fusions observed in cancer occurred years earlier in normal somatic cells.
|
|
| 18. |
- Karmin, Monika, et al.
(author)
-
A recent bottleneck of Y chromosome diversity coincides with a global change in culture.
- 2015
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 25:4
-
Journal article (peer-reviewed)abstract
- It is commonly thought that human genetic diversity in non-African populations was shaped primarily by an out-of-Africa dispersal 50-100 thousand yr ago (kya). Here, we present a study of 456 geographically diverse high-coverage Y chromosome sequences, including 299 newly reported samples. Applying ancient DNA calibration, we date the Y-chromosomal most recent common ancestor (MRCA) in Africa at 254 (95% CI 192-307) kya and detect a cluster of major non-African founder haplogroups in a narrow time interval at 47-52 kya, consistent with a rapid initial colonization model of Eurasia and Oceania after the out-of-Africa bottleneck. In contrast to demographic reconstructions based on mtDNA, we infer a second strong bottleneck in Y-chromosome lineages dating to the last 10 ky. We hypothesize that this bottleneck is caused by cultural changes affecting variance of reproductive success among males.
|
|
| 19. |
- Lee, Hun-Goo, et al.
(author)
-
Genome-wide activities of Polycomb complexes control pervasive transcription
- 2015
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 25:8, s. 1170-1181
-
Journal article (peer-reviewed)abstract
- Polycomb group (PcG) complexes PRC1 and PRC2 are well known for silencing specific developmental genes. PRC2 is a methyltransferase targeting histone H3K27 and producing H3K27me3, essential for stable silencing. Less well known but quantitatively much more important is the genome-wide role of PRC2 that dimethylates similar to 70% of total H3K27. We show that H3K27me2 occurs in inverse proportion to transcriptional activity in most non-PcG target genes and intergenic regions and is governed by opposing roaming activities of PRC2 and complexes containing the H3K27 demethylase UTX. Surprisingly, loss of H3K27me2 results in global transcriptional derepression proportionally greatest in silent or weakly transcribed intergenic and genic regions and accompanied by an increase of H3K27ac and H3K4me1. H3K27me2 therefore sets a threshold that prevents random, unscheduled transcription all over the genome and even limits the activity of highly transcribed genes. PRC1-type complexes also have global roles. Unexpectedly, we find a pervasive distribution of histone H2A ubiquitylated at lysine 118 (H2AK118ub) outside of canonical PcG target regions, dependent on the RING/Sce subunit of PRC1-type complexes. We show, however, that H2AK118ub does not mediate the global PRC2 activity or the global repression and is predominantly produced by a new complex involving L(3) 73Ah, a homolog of mammalian PCGF3.
|
|
| 20. |
- Noorani, A, et al.
(author)
-
A comparative analysis of whole genome sequencing of esophageal adenocarcinoma pre- and post-chemotherapy
- 2017
-
In: Genome research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051. ; 27:6, s. 902-912
-
Journal article (peer-reviewed)abstract
- The scientific community has avoided using tissue samples from patients that have been exposed to systemic chemotherapy to infer the genomic landscape of a given cancer. Esophageal adenocarcinoma is a heterogeneous, chemoresistant tumor for which the availability and size of pretreatment endoscopic samples are limiting. This study compares whole-genome sequencing data obtained from chemo-naive and chemo-treated samples. The quality of whole-genomic sequencing data is comparable across all samples regardless of chemotherapy status. Inclusion of samples collected post-chemotherapy increased the proportion of late-stage tumors. When comparing matched pre- and post-chemotherapy samples from 10 cases, the mutational signatures, copy number, and SNV mutational profiles reflect the expected heterogeneity in this disease. Analysis of SNVs in relation to allele-specific copy-number changes pinpoints the common ancestor to a point prior to chemotherapy. For cases in which pre- and post-chemotherapy samples do show substantial differences, the timing of the divergence is near-synchronous with endoreduplication. Comparison across a large prospective cohort (62 treatment-naive, 58 chemotherapy-treated samples) reveals no significant differences in the overall mutation rate, mutation signatures, specific recurrent point mutations, or copy-number events in respect to chemotherapy status. In conclusion, whole-genome sequencing of samples obtained following neoadjuvant chemotherapy is representative of the genomic landscape of esophageal adenocarcinoma. Excluding these samples reduces the material available for cataloging and introduces a bias toward the earlier stages of cancer.
|
|
| 21. |
- Pettersson, Mats, et al.
(author)
-
A chromosome-level assembly of the Atlantic herring : detection of a supergene and other signals of selection
- 2019
-
In: Genome Research. - : Cold Spring Harbor Laboratory Press (CSHL). - 1088-9051 .- 1549-5469. ; 29:11, s. 1919-1928
-
Journal article (peer-reviewed)abstract
- The Atlantic herring is a model species for exploring the genetic basis for ecological adaptation, due to its huge population size and extremely low genetic differentiation at selectively neutral loci. However, such studies have so far been hampered because of a highly fragmented genome assembly. Here, we deliver a chromosome-level genome assembly based on a hybrid approach combining a de novo Pacific Biosciences (PacBio) assembly with Hi-C-supported scaffolding. The assembly comprises 26 autosomes with sizes ranging from 12.4 to 33.1 Mb and a total size, in chromosomes, of 726 Mb, which has been corroborated by a high-resolution linkage map. A comparison between the herring genome assembly with other high-quality assemblies from bony fishes revealed few inter-chromosomal but frequent intra-chromosomal rearrangements. The improved assembly facilitates analysis of previously intractable large-scale structural variation, allowing, for example, the detection of a 7.8-Mb inversion on Chromosome 12 underlying ecological adaptation. This supergene shows strong genetic differentiation between populations. The chromosome-based assembly also markedly improves the interpretation of previously detected signals of selection, allowing us to reveal hundreds of independent loci associated with ecological adaptation.
|
|
| 22. |
- Schmitges, FW, et al.
(author)
-
Multiparameter functional diversity of human C2H2 zinc finger proteins
- 2016
-
In: Genome research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051. ; 26:12, s. 1742-1752
-
Journal article (peer-reviewed)abstract
- C2H2 zinc finger proteins represent the largest and most enigmatic class of human transcription factors. Their C2H2-ZF arrays are highly variable, indicating that most will have unique DNA binding motifs. However, most of the binding motifs have not been directly determined. In addition, little is known about whether or how these proteins regulate transcription. Most of the ∼700 human C2H2-ZF proteins also contain at least one KRAB, SCAN, BTB, or SET domain, suggesting that they may have common interacting partners and/or effector functions. Here, we report a multifaceted functional analysis of 131 human C2H2-ZF proteins, encompassing DNA binding sites, interacting proteins, and transcriptional response to genetic perturbation. We confirm the expected diversity in DNA binding motifs and genomic binding sites, and provide motif models for 78 previously uncharacterized C2H2-ZF proteins, most of which are unique. Surprisingly, the diversity in protein–protein interactions is nearly as high as diversity in DNA binding motifs: Most C2H2-ZF proteins interact with a unique spectrum of co-activators and co-repressors. Thus, multiparameter diversification likely underlies the evolutionary success of this large class of human proteins.
|
|
| 23. |
- Shah, MY, et al.
(author)
-
Cancer-associated rs6983267 SNP and its accompanying long noncoding RNA CCAT2 induce myeloid malignancies via unique SNP-specific RNA mutations
- 2018
-
In: Genome research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051. ; 28:4, s. 432-447
-
Journal article (peer-reviewed)abstract
- The cancer-risk-associated rs6983267 single nucleotide polymorphism (SNP) and the accompanying long noncoding RNA CCAT2 in the highly amplified 8q24.21 region have been implicated in cancer predisposition, although causality has not been established. Here, using allele-specific CCAT2 transgenic mice, we demonstrate that CCAT2 overexpression leads to spontaneous myeloid malignancies. We further identified that CCAT2 is overexpressed in bone marrow and peripheral blood of myelodysplastic/myeloproliferative neoplasms (MDS/MPN) patients. CCAT2 induces global deregulation of gene expression by down-regulating EZH2 in vitro and in vivo in an allele-specific manner. We also identified a novel non-APOBEC, non-ADAR, RNA editing at the SNP locus in MDS/MPN patients and CCAT2-transgenic mice. The RNA transcribed from the SNP locus in malignant hematopoietic cells have different allelic composition from the corresponding genomic DNA, a phenomenon rarely observed in normal cells. Our findings provide fundamental insights into the functional role of rs6983267 SNP and CCAT2 in myeloid malignancies.
|
|
| 24. |
- Smeds, Linnea, et al.
(author)
-
Direct estimate of the rate of germline mutation in a bird
- 2016
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 26:9, s. 1211-1218
-
Journal article (peer-reviewed)abstract
- The fidelity of DNA replication together with repair mechanisms ensure that the genetic material is properly copied from one generation to another. However, on extremely rare occasions when damages to DNA or replication errors are not repaired, germline mutations can be transmitted to the next generation. Because of the rarity of these events, studying the rate at which new mutations arise across organisms has been a great challenge, especially in multicellular nonmodel organisms with large genomes. We sequenced the genomes of 11 birds from a three-generation pedigree of the collared flycatcher (Ficedula albicollis) and used highly stringent bioinformatic criteria for mutation detection and used several procedures to validate mutations, including following the stable inheritance of new mutations to subsequent generations. We identified 55 de novo mutations with a 10-fold enrichment of mutations at CpG sites and with only a modest male mutation bias. The estimated rate of mutation per site per generation was 4.6 x 10(-9), which corresponds to 2.3 x 10(-9) mutations per site per year. Compared to mammals, this is similar to mouse but about half of that reported for humans, which may be due to the higher frequency of male mutations in humans. We confirm that mutation rate scales positively with genome size and that there is a strong negative relationship between mutation rate and effective population size, in line with the drift-barrier hypothesis. Our study illustrates that it should be feasible to obtain direct estimates of the rate of mutation in essentially any organism from which family material can be obtained.
|
|
| 25. |
- Smid, Marcel, et al.
(author)
-
The circular RNome of primary breast cancer
- 2019
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 29:3, s. 356-366
-
Journal article (peer-reviewed)abstract
- Circular RNAs (circRNAs) are a class of RNAs that is under increasing scrutiny, although their functional roles are debated. We analyzed RNA-seq data of 348 primary breast cancers and developed a method to identify circRNAs that does not rely on unmapped reads or known splice junctions. We identified 95,843 circRNAs, of which 20,441 were found recurrently. Of the circRNAs that match exon boundaries of the same gene, 668 showed a poor or even negative (R <0.2) correlation with the expression level of the linear gene. In silico analysis showed only a minority (8.5%) of circRNAs could be explained by known splicing events. Both these observations suggest that specific regulatory processes for circRNAs exist. We confirmed the presence of circRNAs of CNOT2, CREBBP, and RERE in an independent pool of primary breast cancers. We identified circRNA profiles associated with subgroups of breast cancers and with biological and clinical features, such as amount of tumor lymphocytic infiltrate and proliferation index. siRNA-mediated knockdown of circCNOT2 was shown to significantly reduce viability of the breast cancer cell lines MCF-7 and BT-474, further underlining the biological relevance of circRNAs. Furthermore, we found that circular, and not linear, CNOT2 levels are predictive for progression-free survival time to aromatase inhibitor (AI) therapy in advanced breast cancer patients, and found that circCNOT2 is detectable in cell-free RNA from plasma. We showed that circRNAs are abundantly present, show characteristics of being specifically regulated, are associated with clinical and biological properties, and thus are relevant in breast cancer.
|
|
| 26. |
- Sun, Song, et al.
(author)
-
An extended set of yeast-based functional assays accurately identifies human disease mutations
- 2016
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 26:5, s. 670-680
-
Journal article (peer-reviewed)abstract
- We can now routinely identify coding variants within individual human genomes. A pressing challenge is to determine which variants disrupt the function of disease-associated genes. Both experimental and computational methods exist to predict pathogenicity of human genetic variation. However, a systematic performance comparison between them has been lacking. Therefore, we developed and exploited a panel of 26 yeast-based functional complementation assays to measure the impact of 179 variants (101 disease-and 78 non-disease-associated variants) from 22 human disease genes. Using the resulting reference standard, we show that experimental functional assays in a 1-billion-year diverged model organism can identify pathogenic alleles with significantly higher precision and specificity than current computational methods.
|
|
| 27. |
- Sun, Yu, et al.
(author)
-
Intron Evolution in Neurospora : the role of mutational bias and selection
- 2015
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 25:1, s. 100-110
-
Journal article (peer-reviewed)abstract
- We used comparative and population genomics to study intron evolutionary dynamics in the fungal model genus Neurospora. For our investigation, we used well-annotated genomes of N. crassa, N. discreta, and N. tetrasperma, and 92 resequenced genomes of N. tetrasperma from natural populations. By analyzing the four well-annotated genomes, we identified 9495 intron sites in 7619 orthologous genes. Our data supports nonhomologous end joining (NHEJ) and tandem duplication as mechanisms for intron gains in the genus and the RT-mRNA process as a mechanism for intron loss. We found a moderate intron gain rate (5.78-6.89x10(-13) intron gains per nucleotide site per year) and a high intron loss rate (7.53-13.76x10(-10) intron losses per intron sites per year) as compared to other eukaryotes. The derived intron gains and losses are skewed to high frequencies, relative to neutral SNPs, in natural populations of N. tetrasperma, suggesting that selection is involved in maintaining a high intron turnover. Furthermore, our analyses of the association between intron population-level frequency and genomic features suggest that selection is involved in shaping a 5' intron position bias and a low intron GC content. However, intron sequence analyses suggest that the gained introns were not exposed to recent selective sweeps. Taken together, this work contributes to our understanding of the importance of mutational bias and selection in shaping the intron distribution in eukaryotic genomes.
|
|
| 28. |
- Svensson, JP, et al.
(author)
-
A nucleosome turnover map reveals that the stability of histone H4 Lys20 methylation depends on histone recycling in transcribed chromatin
- 2015
-
In: Genome research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051. ; 25:6, s. 872-883
-
Journal article (peer-reviewed)abstract
- Nucleosome composition actively contributes to chromatin structure and accessibility. Cells have developed mechanisms to remove or recycle histones, generating a landscape of differentially aged nucleosomes. This study aimed to create a high-resolution, genome-wide map of nucleosome turnover in Schizosaccharomyces pombe. The recombination-induced tag exchange (RITE) method was used to study replication-independent nucleosome turnover through the appearance of new histone H3 and the disappearance or preservation of old histone H3. The genome-wide location of histones was determined by chromatin immunoprecipitation–exonuclease methodology (ChIP-exo). The findings were compared with diverse chromatin marks, including histone variant H2A.Z, post-translational histone modifications, and Pol II binding. Finally, genome-wide mapping of the methylation states of H4K20 was performed to determine the relationship between methylation (mono, di, and tri) of this residue and nucleosome turnover. Our analysis showed that histone recycling resulted in low nucleosome turnover in the coding regions of active genes, stably expressed at intermediate levels. High levels of transcription resulted in the incorporation of new histones primarily at the end of transcribed units. H4K20 was methylated in low-turnover nucleosomes in euchromatic regions, notably in the coding regions of long genes that were expressed at low levels. This transcription-dependent accumulation of histone methylation was dependent on the histone chaperone complex FACT. Our data showed that nucleosome turnover is highly dynamic in the genome and that several mechanisms are at play to either maintain or suppress stability. In particular, we found that FACT-associated transcription conserves histones by recycling them and is required for progressive H4K20 methylation.
|
|
| 29. |
- Teppo, Susanna, et al.
(author)
-
Genome-wide repression of eRNA and target gene loci by the ETV6-RUNX1 fusion in acute leukemia
- 2016
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 26:11, s. 1468-1477
-
Journal article (peer-reviewed)abstract
- Approximately 20%-25% of childhood acute lymphoblastic leukemias carry the ETV6-RUNX1 (E/R) fusion gene, a fusion of two central hematopoietic transcription factors, ETV6 (TEL) and RUNX1 (AML1). Despite its prevalence, the exact genomic targets of E/R have remained elusive. We evaluated gene loci and enhancers targeted by E/R genome-wide in precursor B acute leukemia cells using global run-on sequencing (GRO-seq). We show that expression of the E/R fusion leads to widespread repression of RUNX1 motif-containing enhancers at its target gene loci. Moreover, multiple super-enhancers from the CD19(+)/CD20(+)-lineage were repressed, implicating a role in impediment of lineage commitment. In effect, the expression of several genes involved in B cell signaling and adhesion was down-regulated, and the repression depended on the wild-type DNA-binding Runt domain of RUNX1. We also identified a number of E/R-regulated annotated and de novo noncoding genes. The results provide a comprehensive genome-wide mapping between E/R-regulated key regulatory elements and genes in precursor B cell leukemia that disrupt normal B lymphopoiesis.
|
|
| 30. |
- Ungerbäck, Jonas, et al.
(author)
-
Pioneering, chromatin remodeling, and epigenetic constraint in early T-cell gene regulation by SPIT (PU.1)
- 2018
-
In: Genome Research. - : COLD SPRING HARBOR LAB PRESS, PUBLICATIONS DEPT. - 1088-9051 .- 1549-5469. ; 28:10, s. 1508-1519
-
Journal article (peer-reviewed)abstract
- SPI1 (also known as PU.1) is a dominant but transient regulator in early T-cell precursors and a potent transcriptional controller of developmentally important pro-T-cell genes. Before T-lineage commitment, open chromatin is frequently occupied by PU.1, and many PU.1 sites lose accessibility when PU.1 is later down-regulated. Pioneering activity of PU.1 was tested in this developmentally dynamic context by quantitating the relationships between PU.1 occupancy and site quality and accessibility as PU.1 levels naturally declined in pro-T-cell development and by using stage-specific gain- and loss-offunction perturbations to relate binding to effects on target genes. PU.1 could bind closed genomic sites, but rapidly opened many of them, despite the absence of its frequent collaborator, CEBPA. RUNX motifs and RUNX1 binding were often linked to PU.1 at open sites, but highly expressed PU.1 could bind its sites without RUNX1 The dynamic properties of PU.1 engagements implied that PU.1 binding affinity and concentration determine its occupancy choices, but with quantitative trade-offs for occupancy between site sequence quality and stage-dependent site accessibility in chromatin. At nonpromoter sites, PU.1 binding criteria were more stringent than at promoters, and PU.1 was also much more effective as a transcriptional regulator at non promoter sites where local chromatin accessibility depended on the presence of PU.1. Notably, closed chromatin presented a qualitative barrier to occupancy by the PU.1 DNA-binding domain alone. Thus, effective pioneering at closed chromatin sites also depends on requirements beyond site recognition, served by non-DNA-binding domains of PU.1.
|
|
| 31. |
- Varshney, Gaurav K., et al.
(author)
-
High-throughput gene targeting and phenotyping in zebrafish using CRISPR/Cas9
- 2015
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 25:7, s. 1030-1042
-
Journal article (peer-reviewed)abstract
- The use of CRISPR/Cas9 as a genome-editing tool in various model organisms has radically changed targeted mutagenesis. Here, we present a high-throughput targeted mutagenesis pipeline using CRISPR/Cas9 technology in zebrafish that will make possible both saturation mutagenesis of the genome and large-scale phenotyping efforts. We describe a cloning-free single-guide RNA (sgRNA) synthesis, coupled with streamlined mutant identification methods utilizing fluorescent PCR and multiplexed, high-throughput sequencing. We report germline transmission data from 162 loci targeting 83 genes in the zebrafish genome, in which we obtained a 99% success rate for generating mutations and an average germline transmission rate of 28%. We verified 678 unique alleles from 58 genes by high-throughput sequencing. We demonstrate that our method can be used for efficient multiplexed gene targeting. We also demonstrate that phenotyping can be done in the F-1 generation by inbreeding two injected founder fish, significantly reducing animal husbandry and time. This study compares germline transmission data from CRISPR/Cas9 with those of TALENs and ZFNs and shows that efficiency of CRISPR/Cas9 is sixfold more efficient than other techniques. We show that the majority of published "rules" for efficient sgRNA design do not effectively predict germline transmission rates in zebrafish, with the exception of a GG or GA dinucleotide genomic match at the 5' end of the sgRNA. Finally, we show that predicted off-target mutagenesis is of low concern for in vivo genetic studies.
|
|
| 32. |
- Wallerman, Ola, et al.
(author)
-
A chromosome-level assembly of the Atlantic herring genome-detection of a supergene and other signals of selection
- 2019
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 29, s. 1919-1928
-
Journal article (peer-reviewed)abstract
- The Atlantic herring is a model species for exploring the genetic basis for ecological adaptation, due to its huge population size and extremely low genetic differentiation at selectively neutral loci. However, such studies have so far been hampered because of a highly fragmented genome assembly. Here, we deliver a chromosome-level genome assembly based on a hybrid approach combining a de novo Pacific Biosciences (PacBio) assembly with Hi-C-supported scaffolding. The assembly comprises 26 autosomes with sizes ranging from 12.4 to 33.1 Mb and a total size, in chromosomes, of 726 Mb, which has been corroborated by a high-resolution linkage map. A comparison between the herring genome assembly with other high-quality assemblies from bony fishes revealed few inter-chromosomal but frequent intra-chromosomal rearrangements. The improved assembly facilitates analysis of previously intractable large-scale structural variation, allowing, for example, the detection of a 7.8-Mb inversion on Chromosome 12 underlying ecological adaptation. This supergene shows strong genetic differentiation between populations. The chromosome-based assembly also markedly improves the interpretation of previously detected signals of selection, allowing us to reveal hundreds of independent loci associated with ecological adaptation.
|
|
| 33. |
- Warnefors, Maria, et al.
(author)
-
Sex-biased microRNA expression in mammals and birds reveals underlying regulatory mechanisms and a role in dosage compensation
- 2017
-
In: Genome Research. - : COLD SPRING HARBOR LAB PRESS, PUBLICATIONS DEPT. - 1088-9051 .- 1549-5469. ; 27:12, s. 1961-1973
-
Journal article (peer-reviewed)abstract
- Sexual dimorphism depends on sex-biased gene expression, but the contributions of microRNAs (miRNAs) have not been globally assessed. We therefore produced an extensive small RNA sequencing data set to analyze male and female miRNA expression profiles in mouse, opossum, and chicken. Our analyses uncovered numerous cases of somatic sex-biased miRNA expression, with the largest proportion found in the mouse heart and liver. Sex-biased expression is explained by miRNA-specific regulation, including sex-biased chromatin accessibility at promoters, rather than piggybacking of intronic miRNAs on sex-biased protein-coding genes. In mouse, but not opossum and chicken, sex bias is coordinated across tissues such that autosomal testis-biased miRNAs tend to be somatically male-biased, whereas autosomal ovary-biased miRNAs are female-biased, possibly due to broad hormonal control. In chicken, which has a Z/W sex chromosome system, expression output of genes on the Z Chromosome is expected to be male-biased, since there is no global dosage compensation mechanism that restores expression in ZW females after almost all genes on the W Chromosome decayed. Nevertheless, we found that the dominant liver miRNA, miR-122-5p, is Z-linked but expressed in an unbiased manner, due to the unusual retention of a W-linked copy. Another Z-linked miRNA, the male-biased miR-2954-3p, shows conserved preference for dosage-sensitive genes on the Z Chromosome, based on computational and experimental data from chicken and zebra finch, and acts to equalize male-to-female expression ratios of its targets. Unexpectedly, our findings thus establish miRNA regulation as a novel gene-specific dosage compensation mechanism.
|
|
| 34. |
- Wei, W, et al.
(author)
-
Chromatin-sensitive cryptic promoters putatively drive expression of alternative protein isoforms in yeast
- 2019
-
In: Genome research. - : Cold Spring Harbor Laboratory. - 1549-5469 .- 1088-9051. ; 29:12, s. 1974-1984
-
Journal article (peer-reviewed)abstract
- Cryptic transcription is widespread and generates a heterogeneous group of RNA molecules of unknown function. To improve our understanding of cryptic transcription, we investigated their transcription start site (TSS) usage, chromatin organization, and posttranscriptional consequences in Saccharomyces cerevisiae. We show that TSSs of chromatin-sensitive internal cryptic transcripts retain comparable features of canonical TSSs in terms of DNA sequence, directionality, and chromatin accessibility. We define the 5′ and 3′ boundaries of cryptic transcripts and show that, contrary to RNA degradation–sensitive ones, they often overlap with the end of the gene, thereby using the canonical polyadenylation site, and associate to polyribosomes. We show that chromatin-sensitive cryptic transcripts can be recognized by ribosomes and may produce truncated polypeptides from downstream, in-frame start codons. Finally, we confirm the presence of the predicted polypeptides by reanalyzing N-terminal proteomic data sets. Our work suggests that a fraction of chromatin-sensitive internal cryptic promoters initiates the transcription of alternative truncated mRNA isoforms. The expression of these chromatin-sensitive isoforms is conserved from yeast to human, expanding the functional consequences of cryptic transcription and proteome complexity.
|
|
| 35. |
- Weissensteiner, Matthias H., et al.
(author)
-
Combination of short-read, long-read, and optical mapping assemblies reveals large-scale tandem repeat arrays with population genetic implications
- 2017
-
In: Genome Research. - : COLD SPRING HARBOR LAB PRESS, PUBLICATIONS DEPT. - 1088-9051 .- 1549-5469. ; 27:5, s. 697-708
-
Journal article (peer-reviewed)abstract
- Accurate and contiguous genome assembly is key to a comprehensive understanding of the processes shaping genomic diversity and evolution. Yet, it is frequently constrained by constitutive heterochromatin, usually characterized by highly repetitive DNA. As a key feature of genome architecture associated with centromeric and subtelomeric regions, it locally influences meiotic recombination. In this study, we assess the impact of large tandem repeat arrays on the recombination rate landscape in an avian speciation model, the Eurasian crow. We assembled two high-quality genome references using single-molecule real-time sequencing (long-read assembly [LR]) and single-molecule optical maps (optical map assembly [ OM]). A three-way comparison including the published short-read assembly (SR) constructed for the same individual allowed assessing assembly properties and pinpointing misassemblies. By combining information from all three assemblies, we characterized 36 previously unidentified large repetitive regions in the proximity of sequence assembly breakpoints, the majority of which contained complex arrays of a 14-kb satellite repeat or its 1.2-kb subunit. Using whole-genome population resequencing data, we estimated the population-scaled recombination rate (rho) and found it to be significantly reduced in these regions. These findings are consistent with an effect of low recombination in regions adjacent to centromeric or subtelomeric heterochromatin and add to our understanding of the processes generating widespread heterogeneity in genetic diversity and differentiation along the genome. By combining three different technologies, our results highlight the importance of adding a layer of information on genome structure that is inaccessible to each approach independently.
|
|
| 36. |
- Zhang, Tongwu, et al.
(author)
-
Cell-type-specific eQTL of primary melanocytes facilitates identification of melanoma susceptibility genes
- 2018
-
In: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 28:11, s. 1621-1635
-
Journal article (peer-reviewed)abstract
- Most expression quantitative trait locus (eQTL) studies to date have been performed in heterogeneous tissues as opposed to specific cell types. To better understand the cell-type-specific regulatory landscape of human melanocytes, which give rise to melanoma but account for <5% of typical human skin biopsies, we performed an eQTL analysis in primary melanocyte cultures from 106 newborn males. We identified 597,335 cis-eQTL SNPs prior to linkage disequilibrium (LD) pruning and 4997 eGenes (FDR < 0.05). Melanocyte eQTLs differed considerably from those identified in the 44 GTEx tissue types, including skin. Over a third of melanocyte eGenes, including key genes in melanin synthesis pathways, were unique to melanocytes compared to those of GTEx skin tissues or TCGA melanomas. The melanocyte data set also identified trans-eQTLs, including those connecting a pigmentation-associated functional SNP with four genes, likely through cis-regulation of IRF4. Melanocyte eQTLs are enriched in cis-regulatory signatures found in melanocytes as well as in melanoma-associated variants identified through genome-wide association studies. Melanocyte eQTLs also colocalized with melanoma GWAS variants in five known loci. Finally, a transcriptome-wide association study using melanocyte eQTLs uncovered four novel susceptibility loci, where imputed expression levels of five genes (ZFP90, HEBP1, MSC, CBWD1, and RP11-383H13.1) were associated with melanoma at genome-wide significant P-values. Our data highlight the utility of lineage-specific eQTL resources for annotating GWAS findings, and present a robust database for genomic research of melanoma risk and melanocyte biology.
|
|
