Evolution of Copy Number at the MHC Varies across the Avian Tree of Life.
Ontology highlight
ABSTRACT:
SUBMITTER: Minias P
PROVIDER: S-EPMC6319602 | biostudies-literature | 2019 Jan
REPOSITORIES: biostudies-literature
Publications
Evolution of Copy Number at the MHC Varies across the Avian Tree of Life.
Minias Piotr P Pikus Ewa E Whittingham Linda A LA Dunn Peter O PO
Genome biology and evolution 20190101 1
The evolution of the major histocompatibility complex (MHC) is shaped by frequent gene duplications and deletions, which generate extensive variation in the number of loci (gene copies) between different taxa. Here, we collected estimates of copy number at the MHC for over 250 bird species from 68 families. We found contrasting patterns of copy number evolution between MHC class I and class IIB, which encode receptors for intra- and extracellular pathogens, respectively. Across the avian evoluti ...[more]
Similar Datasets
Project description:All organisms encounter pathogens, and birds are especially susceptible to infection by malaria parasites and other haemosporidians. It is important to understand how immune genes, primarily innate immune genes which are the first line of host defense, have evolved across birds, a highly diverse group of tetrapods. Here, we find that innate immune genes are highly conserved across the avian tree of life and that although most show evidence of positive or diversifying selection within specific lineages or clades, the number of sites is often proportionally low in this broader context of putative constraint. Rather, the evidence shows a much higher level of negative or purifying selection in these innate immune genes - rather than adaptive immune genes - which is consistent with birds' long coevolutionary history with pathogens and the need to maintain a rapid response to infection. We further explored avian responses to haemosporidians by comparing differential gene expression in wild birds (1) uninfected with haemosporidians, (2) infected with Plasmodium, and (3) infected with Haemoproteus (Parahaemoproteus). We found patterns of significant differential expression with some genes unique to infection with each genus and a few shared between "treatment" groups, but none that overlapped with the genes included in the phylogenetic study.
| S-EPMC6890419 | biostudies-literature
Project description: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.
| S-EPMC2847753 | biostudies-literature
Project description:tRNAs are universal decoders that bridge the gap between transcriptome and proteome. They can also be processed into small RNA fragments with regulatory functions. In this work, we show that tRNA copy number is largely controlled by genome size in all cellular organisms, in contrast to what is observed for protein-coding genes that stop expanding between ~20,000 and ~35,000 loci per haploid genome in eukaryotes, regardless of genome size. Our analyses indicate that after the bacteria/archaea split, the tRNA gene pool experienced the evolution of increased anticodon diversity in the archaeal lineage, along with a tRNA gene size increase and mature tRNA size decrease. The evolution and diversification of eukaryotes from archaeal ancestors involved further expansion of the tRNA anticodon repertoire, additional increase in tRNA gene size and decrease in mature tRNA length, along with an explosion of the tRNA gene copy number that emerged coupled with accelerated genome size expansion. Our findings support the notion that macroscopic eukaryotes with a high diversity of cell types, such as land plants and vertebrates, independently evolved a high diversity of tRNA anticodons along with high gene redundancy caused by the expansion of the tRNA copy number. The results presented here suggest that the evolution of tRNA genes played important roles in the early split between bacteria and archaea, and in eukaryogenesis and the later emergence of complex eukaryotes, with potential implications in protein translation and gene regulation through tRNA-derived RNA fragments.
| S-EPMC9858662 | biostudies-literature
Project description:Premise of the studyPolyploidy has profound evolutionary consequences for land plants. Despite the availability of large phylogenetic and chromosomal data sets, estimating the rates of polyploidy and chromosomal evolution across the tree of life remains a challenging, computationally complex problem. We introduce the R package chromploid, which allows scientists to perform inference of chromosomal evolution rates across large phylogenetic trees.Methods and resultschromploid is an open-source package in the R environment that calculates the likelihood function of models of chromosome evolution. Models of discrete character evolution can be customized using chromploid. We demonstrate the performance of the BiChroM model, testing for associations between rates of chromosome doubling (as a proxy for polyploidy) and a binary phenotypic character, within chromploid using simulations and empirical data from Solanum. In simulations, estimated chromosome-doubling rates were unbiased and the variance decreased with larger trees, but distinguishing small differences in rates of chromosome doubling, even from large data sets, remains challenging. In the Solanum data set, a custom model of chromosome number evolution demonstrated higher rates of chromosome doubling in herbaceous species compared to woody.Conclusionschromploid enables researchers to perform robust likelihood-based inferences using complex models of chromosome number evolution across large phylogenies.
| S-EPMC5895187 | biostudies-literature
Project description:The human UGT2B17 gene varies in copy number from zero to two per individual and also differs in mean number between populations from Africa, Europe, and East Asia. We show that such a high degree of geographical variation is unusual and investigate its evolutionary history. This required first reinterpreting the reference sequence in this region of the genome, which is misassembled from the two different alleles separated by an artifactual gap. A corrected assembly identifies the polymorphism as a 117 kb deletion arising by nonallelic homologous recombination between approximately 4.9 kb segmental duplications and allows the deletion breakpoint to be identified. We resequenced approximately 12 kb of DNA spanning the breakpoint in 91 humans from three HapMap and one extended HapMap populations and one chimpanzee. Diversity was unusually high and the time to the most recent common ancestor was estimated at approximately 2.4 or approximately 3.0 million years by two different methods, with evidence of balancing selection in Europe. In contrast, diversity was low in East Asia where a single haplotype predominated, suggesting positive selection for the deletion in this part of the world.
| S-EPMC2556428 | biostudies-literature
Project description:Copy number polymorphism (CNP) is ubiquitous in eukaryotic genomes, but the degree to which this reflects the action of positive selection is poorly understood. The first gene in the Plasmodium folate biosynthesis pathway, GTP-cyclohydrolase I (gch1), shows extensive CNP. We provide compelling evidence that gch1 CNP is an adaptive consequence of selection by antifolate drugs, which target enzymes downstream in this pathway. (1) We compared gch1 CNP in parasites from Thailand (strong historical antifolate selection) with those from neighboring Laos (weak antifolate selection). Two percent of chromosomes had amplified copy number in Laos, while 72% carried multiple (2-11) copies in Thailand, and differentiation exceeded that observed at 73 synonymous SNPs. (2) We found five amplicon types containing one to greater than six genes and spanning 1 to >11 kb, consistent with parallel evolution and strong selection for this gene amplification. gch1 was the only gene occurring in all amplicons suggesting that this locus is the target of selection. (3) We observed reduced microsatellite variation and increased linkage disequilibrium (LD) in a 900-kb region flanking gch1 in parasites from Thailand, consistent with rapid recent spread of chromosomes carrying multiple copies of gch1. (4) We found that parasites bearing dhfr-164L, which causes high-level resistance to antifolate drugs, carry significantly (p = 0.00003) higher copy numbers of gch1 than parasites bearing 164I, indicating functional association between genes located on different chromosomes but linked in the same biochemical pathway. These results demonstrate that CNP at gch1 is adaptive and the associations with dhfr-164L strongly suggest a compensatory function. More generally, these data demonstrate how selection affects multiple enzymes in a single biochemical pathway, and suggest that investigation of structural variation may provide a fast-track to locating genes underlying adaptation.
| S-EPMC2570623 | biostudies-literature
Project description:BackgroundMicroinversions are cytologically undetectable inversions of DNA sequences that accumulate slowly in genomes. Like many other rare genomic changes (RGCs), microinversions are thought to be virtually homoplasy-free evolutionary characters, suggesting that they may be very useful for difficult phylogenetic problems such as the avian tree of life. However, few detailed surveys of these genomic rearrangements have been conducted, making it difficult to assess this hypothesis or understand the impact of microinversions upon genome evolution.ResultsWe surveyed non-coding sequence data from a recent avian phylogenetic study and found substantially more microinversions than expected based upon prior information about vertebrate inversion rates, although this is likely due to underestimation of these rates in previous studies. Most microinversions were lineage-specific or united well-accepted groups. However, some homoplastic microinversions were evident among the informative characters. Hemiplasy, which reflects differences between gene trees and the species tree, did not explain the observed homoplasy. Two specific loci were microinversion hotspots, with high numbers of inversions that included both the homoplastic as well as some overlapping microinversions. Neither stem-loop structures nor detectable sequence motifs were associated with microinversions in the hotspots.ConclusionsMicroinversions can provide valuable phylogenetic information, although power analysis indicates that large amounts of sequence data will be necessary to identify enough inversions (and similar RGCs) to resolve short branches in the tree of life. Moreover, microinversions are not perfect characters and should be interpreted with caution, just as with any other character type. Independent of their use for phylogenetic analyses, microinversions are important because they have the potential to complicate alignment of non-coding sequences. Despite their low rate of accumulation, they have clearly contributed to genome evolution, suggesting that active identification of microinversions will prove useful in future phylogenomic studies.
| S-EPMC3123225 | biostudies-literature
Project description:BACKGROUND:Structural rearrangements of the genome resulting in genic imbalance due to copy number change are often deleterious at the organismal level, but are common in immortalized cell lines and tumors, where they may be an advantage to cells. In order to explore the biological consequences of copy number changes in the Drosophila genome, we resequenced the genomes of 19 tissue-culture cell lines and generated RNA-Seq profiles. RESULTS:Our work revealed dramatic duplications and deletions in all cell lines. We found three lines of evidence indicating that copy number changes were due to selection during tissue culture. First, we found that copy numbers correlated to maintain stoichiometric balance in protein complexes and biochemical pathways, consistent with the gene balance hypothesis. Second, while most copy number changes were cell line-specific, we identified some copy number changes shared by many of the independent cell lines. These included dramatic recurrence of increased copy number of the PDGF/VEGF receptor, which is also over-expressed in many cancer cells, and of bantam, an anti-apoptosis miRNA. Third, even when copy number changes seemed distinct between lines, there was strong evidence that they supported a common phenotypic outcome. For example, we found that proto-oncogenes were over-represented in one cell line (S2-DRSC), whereas tumor suppressor genes were under-represented in another (Kc167). CONCLUSION:Our study illustrates how genome structure changes may contribute to selection of cell lines in vitro. This has implications for other cell-level natural selection progressions, including tumorigenesis.
| S-EPMC4289277 | biostudies-literature
Project description:MHC genes, which code for proteins responsible for presenting pathogen-derived antigens to the host immune system, show remarkable copy-number variation both between and within species. However, the evolutionary forces driving this variation are poorly understood. Here, we use computer simulations to investigate whether evolution of the number of MHC variants in the genome can be shaped by the number of pathogen species the host population encounters (pathogen richness). Our model assumed that while increasing a range of pathogens recognised, expressing additional MHC variants also incurs costs such as an increased risk of autoimmunity. We found that pathogen richness selected for high MHC copy number only when the costs were low. Furthermore, the shape of the association was modified by the rate of pathogen evolution, with faster pathogen mutation rates selecting for increased host MHC copy number, but only when pathogen richness was low to moderate. Thus, taking into account factors other than pathogen richness may help explain wide variation between vertebrate species in the number of MHC genes. Within population, variation in the number of unique MHC variants carried by individuals (INV) was observed under most parameter combinations, except at low pathogen richness. This variance gave rise to positive correlations between INV and host immunocompetence (proportion of pathogens recognised). However, within-population variation in host immunocompetence declined with pathogen richness. Thus, counterintuitively, pathogens can contribute more to genetic variance for host fitness in species exposed to fewer pathogen species, with consequences to predictions from "Hamilton-Zuk" theory of sexual selection.
| S-EPMC6541284 | biostudies-literature
Project description:MotivationCopy number aberrations (CNAs), which delete or amplify large contiguous segments of the genome, are a common type of somatic mutation in cancer. Copy number profiles, representing the number of copies of each region of a genome, are readily obtained from whole-genome sequencing or microarrays. However, modeling copy number evolution is a substantial challenge, because different CNAs may overlap with one another on the genome. A recent popular model for copy number evolution is the copy number distance (CND), defined as the length of a shortest sequence of deletions and amplifications of contiguous segments that transforms one profile into the other. In the CND, all events contribute equally; however, it is well known that rates of CNAs vary by length, genomic position and type (amplification versus deletion).ResultsWe introduce a weighted CND that allows events to have varying weights, or probabilities, based on their length, position and type. We derive an efficient algorithm to compute the weighted CND as well as the associated transformation. This algorithm is based on the observation that the constraint matrix of the underlying optimization problem is totally unimodular. We show that the weighted CND improves phylogenetic reconstruction on simulated data where CNAs occur with varying probabilities, aids in the derivation of phylogenies from ultra-low-coverage single-cell DNA sequencing data and helps estimate CNA rates in a large pan-cancer dataset.Availability and implementationCode is available at https://github.com/raphael-group/WCND.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC7355291 | biostudies-literature