Dataset Information

Genome-Scale Profiling Reveals Noncoding Loci Carry Higher Proportions of Concordant Data.

ABSTRACT: Many evolutionary relationships remain controversial despite whole-genome sequencing data. These controversies arise, in part, due to challenges associated with accurately modeling the complex phylogenetic signal coming from genomic regions experiencing distinct evolutionary forces. Here, we examine how different regions of the genome support or contradict well-established relationships among three mammal groups using millions of orthologous parsimony-informative biallelic sites (PIBS) distributed across primate, rodent, and Pecora genomes. We compared PIBS concordance percentages among locus types (e.g. coding sequences (CDS), introns, intergenic regions), and contrasted PIBS utility over evolutionary timescales. Sites derived from noncoding sequences provided more data and proportionally more concordant sites compared with those from CDS in all clades. CDS PIBS were also predominant drivers of tree incongruence in two cases of topological conflict. PIBS derived from most locus types provided surprisingly consistent support for splitting events spread across the timescales we examined, although we find evidence that CDS and intronic PIBS may, respectively and to a limited degree, inform disproportionately about older and younger splits. In this era of accessible wholegenome sequence data, these results:1) suggest benefits to more intentionally focusing on noncoding loci as robust data for tree inference and 2) reinforce the importance of accurate modeling, especially when using CDS data.

SUBMITTER: Literman R

PROVIDER: S-EPMC8136493 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:MotivationWe have proposed a mixture model based approach to the concordant integrative analysis of multiple large-scale two-sample expression datasets. Since the mixture model is based on the transformed differential expression test P-values (z-scores), it is generally applicable to the expression data generated by either microarray or RNA-seq platforms. The mixture model is simple with three normal distribution components for each dataset to represent down-regulation, up-regulation and no differential expression. However, when the number of datasets increases, the model parameter space increases exponentially due to the component combination from different datasets.ResultsIn this study, motivated by the well-known generalized estimating equations (GEEs) for longitudinal data analysis, we focus on the concordant components and assume that the proportions of non-concordant components follow a special structure. We discuss the exchangeable, multiset coefficient and autoregressive structures for model reduction, and their related expectation-maximization (EM) algorithms. Then, the parameter space is linear with the number of datasets. In our previous study, we have applied the general mixture model to three microarray datasets for lung cancer studies. We show that more gene sets (or pathways) can be detected by the reduced mixture model with the exchangeable structure. Furthermore, we show that more genes can also be detected by the reduced model. The Cancer Genome Atlas (TCGA) data have been increasingly collected. The advantage of incorporating the concordance feature has also been clearly demonstrated based on TCGA RNA sequencing data for studying two closely related types of cancer.Availability and implementationAdditional results are included in a supplemental file. Computer program R-functions are freely available at http://home.gwu.edu/?ylai/research/Concordance.Contactylai@gwu.edu.Supplementary informationSupplementary data are available at Bioinformatics online.

Project description:It is increasingly being recognized that physical interactions between distant chromosomal loci influence the expressivity of the genome. Thus, multiple enhancers may converge on a single gene promoter 1,2 and a single enhancer can stochastically loop to multiple promoters 3. Regulatory elements from neighboring domains or from other chromosomes may interact to generate chromatin structures poised for transcription4-8 or repression 9-11. However, a broader perspective on the formation and function of such networks remains largely unexplored. Using the high-resolution 4C (Circular Chromosome Conformation Capture) analysis, we identify here a genome-wide physical network between genomically imprinted domains impinging H19 imprinting control region (ICR). Sequences interacting with the H19 ICR constitutively map at or within most of known imprinted domains both in mouse embryonic stem cells and in vitro-derived embryoid bodies, despite distinct differences in chromatin environments. Moreover, all-to-all 3D DNA FISH analyses revealed that imprinted domains from 7 different chromosomes co-localize with each other in pair-wise patterns, highlighting the dynamic nature of these interactions. The observed interaction network is dependent on CTCF binding sites within the maternally inherited H19 ICR allele and delays reprogramming of epigenetic features at imprinted domains during male germline development. As other imprinted regions can interact with each other in an H19 ICR-dependent manner, the H19 ICR emerges as an organizer of chromosomal networks by specifically recognizing and reconfiguring imprinted domains. We propose that an ancient key regulatory region, such as the H19 ICR, gradually impacted epigenetic modifications via long-range interactions in the germline to facilitate the recruitment of imprinted genes and their formation in clusters during mammalian radiation. Furthermore, the perspective of epigenetic lesions de-stabilizing such chromosomal networks may advance our understanding of how quantitative traits influence human diseases, such as cancer.

Dataset Information

Genome-Scale Profiling Reveals Noncoding Loci Carry Higher Proportions of Concordant Data.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets