RESUMO
PREMISE: A complicating factor in analyzing allopolyploid genomes is the possibility of physical interactions between homoeologous chromosomes during meiosis, resulting in either crossover (homoeologous exchanges) or non-crossover products (homoeologous gene conversion). Homoeologous gene conversion was first described in cotton by comparing SNP patterns in sequences from two diploid progenitors with those from the allopolyploid subgenomes. These analyses, however, did not explicitly consider other evolutionary scenarios that may give rise to similar SNP patterns as homoeologous gene conversion, creating uncertainties about the reality of the inferred gene conversion events. METHODS: Here, we use an expanded phylogenetic sampling of high-quality genome assemblies from seven allopolyploid Gossypium species (all derived from the same polyploidy event), four diploid species (two closely related to each subgenome), and a diploid outgroup to derive a robust method for identifying potential genomic regions of gene conversion and homoeologous exchange. RESULTS: We found little evidence for homoeologous gene conversion in allopolyploid cottons, and that only two of the 40 best-supported events were shared by more than one species. We did, however, reveal a single, shared homoeologous exchange event at one end of chromosome 1, which occurred shortly after allopolyploidization but prior to divergence of the descendant species. CONCLUSIONS: Overall, our analyses demonstrated that homoeologous gene conversion and homoeologous exchanges are uncommon in Gossypium, affecting between zero and 24 genes per subgenome (0.0-0.065%) across the seven species. More generally, we highlighted the potential problems of using simple four-taxon tests to investigate patterns of homoeologous gene conversion in established allopolyploids.
Assuntos
Conversão Gênica , Gossypium , Filogenia , Poliploidia , Gossypium/genética , Genoma de Planta , Polimorfismo de Nucleotídeo Único , Diploide , Genes de Plantas , Cromossomos de Plantas/genéticaRESUMO
CRISPR-Cas9 homing gene drives are designed to induce a targeted double-stranded DNA break at a wild type allele ('recipient'), which, when repaired by the host cell, is converted to the drive allele from the homologous ('donor') chromosome. Germline localisation of this process leads to super-Mendelian inheritance of the drive and the rapid spread of linked traits, offering a novel strategy for population control through the deliberate release of drive individuals. During the homology-based DNA repair, additional segments of the recipient chromosome may convert to match the donor, potentially impacting carrier fitness and strategy success. Using Anopheles gambiae strains with variations around the drive target site, here we assess the extent and nature of chromosomal conversion. We show both homing and meiotic drive contribute as mechanisms of inheritance bias. Additionally, over 80% of homing events resolve within 50 bp of the chromosomal break, enabling rapid gene drive transfer into locally-adapted genetic backgrounds.
Assuntos
Anopheles , Sistemas CRISPR-Cas , Tecnologia de Impulso Genético , Anopheles/genética , Animais , Tecnologia de Impulso Genético/métodos , Feminino , Alelos , Conversão Gênica , Meiose/genética , Masculino , Quebras de DNA de Cadeia Dupla , Cromossomos de Insetos/genéticaRESUMO
The major histocompatibility complex (MHC) multigene family encodes key pathogen-recognition molecules of the vertebrate adaptive immune system. Hyper-polymorphism of MHC genes is de novo generated by point mutations, but new haplotypes may also arise by re-shuffling of existing variation through intra- and inter-locus gene conversion. Although the occurrence of gene conversion at the MHC has been known for decades, we still have limited understanding of its functional importance. Here, I took advantage of extensive genetic resources (~9000 sequences) to investigate broad scale macroevolutionary patterns in gene conversion processes at the MHC across nearly 200 avian species. Gene conversion was found to constitute a universal mechanism in birds, as 83% of species showed footprints of gene conversion at either MHC class and 25% of all allelic variants were attributed to gene conversion. Gene conversion processes were stronger at MHC-II than MHC-I, but inter-specific variation at both MHC classes was explained by similar evolutionary scenarios, reflecting fluctuating selection towards different optima and drift. Gene conversion showed uneven phylogenetic distribution across birds and was driven by gene copy number variation, supporting significant role of inter-locus gene conversion processes in the evolution of the avian MHC. Finally, MHC gene conversion was stronger in species with fast life histories (high fecundity) and in long-distance migrants, likely reflecting variation in population sizes and host-pathogen coevolutionary dynamics. The results provide a robust comparative framework for understanding macroevolutionary variation in gene conversion at the avian MHC and reinforce important contribution of this mechanism to functional MHC diversity.
Assuntos
Aves , Evolução Molecular , Conversão Gênica , Complexo Principal de Histocompatibilidade , Filogenia , Seleção Genética , Animais , Aves/genética , Complexo Principal de Histocompatibilidade/genética , Seleção Genética/genética , Dosagem de Genes , Haplótipos/genética , Variação GenéticaRESUMO
The role of balancing selection is a long-standing evolutionary puzzle. Balancing selection is a crucial evolutionary process that maintains genetic variation (polymorphism) over extended periods of time; however, detecting it poses a significant challenge. Building upon the Polymorphism-aware phylogenetic Models (PoMos) framework rooted in the Moran model, we introduce a PoMoBalance model. This novel approach is designed to disentangle the interplay of mutation, genetic drift, and directional selection (GC-biased gene conversion), along with the previously unexplored balancing selection pressures on ultra-long timescales comparable with species divergence times by analyzing multi-individual genomic and phylogenetic divergence data. Implemented in the open-source RevBayes Bayesian framework, PoMoBalance offers a versatile tool for inferring phylogenetic trees as well as quantifying various selective pressures. The novel aspect of our approach in studying balancing selection lies in polymorphism-aware phylogenetic models' ability to account for ancestral polymorphisms and incorporate parameters that measure frequency-dependent selection, allowing us to determine the strength of the effect and exact frequencies under selection. We implemented validation tests and assessed the model on the data simulated with SLiM and a custom Moran model simulator. Real sequence analysis of Drosophila populations reveals insights into the evolutionary dynamics of regions subject to frequency-dependent balancing selection, particularly in the context of sex-limited color dimorphism in Drosophila erecta.
Assuntos
Conversão Gênica , Modelos Genéticos , Filogenia , Polimorfismo Genético , Seleção Genética , Animais , Teorema de Bayes , Evolução Molecular , Masculino , FemininoRESUMO
The model haloarchaeon Haloferax volcanii is polyploid with about 20 copies of its major chromosome. Recently it has been described that highly efficient intermolecular gene conversion operates in H. volcanii to equalize the chromosomal copies. In the current study, 24 genes were selected that encode proteins with orthologs involved in gene conversion or homologous recombination in archaea, bacteria, or eukaryotes. Single gene deletion strains of 22 genes and a control gene were constructed in two parent strains for a gene conversion assay; only radA and radB were shown to be essential. Protoplast fusions were used to generate strains that were heterozygous for the gene HVO_2528, encoding an enzyme for carotinoid biosynthesis. It was revealed that a lack of six of the proteins did not influence the efficiency of gene conversion, while sixteen mutants had severe gene conversion defects. Notably, lack of paralogous proteins of gene families had very different effects, e.g., mutant Δrad25b had no phenotype, while mutants Δrad25a, Δrad25c, and Δrad25d were highly compromised. Generation of a quadruple rad25 and a triple sph deletion strain also indicated that the paralogs have different functions, in contrast to sph2 and sph4, which cannot be deleted simultaneously. There was no correlation between the severity of the phenotypes and the respective transcript levels under non-stressed conditions, indicating that gene expression has to be induced at the onset of gene conversion. Phylogenetic trees of the protein families Rad3/25, MutL/S, and Sph/SMC/Rad50 were generated to unravel the history of the paralogous proteins of H. volcanii. Taken together, unselected intermolecular gene conversion in H. volcanii involves at least 16 different proteins, the molecular roles of which can be studied in detail in future projects.
Assuntos
Proteínas Arqueais , Conversão Gênica , Haloferax volcanii , Haloferax volcanii/genética , Proteínas Arqueais/genética , Proteínas Arqueais/metabolismo , Poliploidia , Genoma Arqueal/genética , Deleção de Genes , Dosagem de GenesRESUMO
Long-read sequencing, exemplified by PacBio, revolutionizes genomics, overcoming challenges like repetitive sequences. However, the high DNA requirement ( > 1 µg) is prohibitive for small organisms. We develop a low-input (100 ng), low-cost, and amplification-free library-generation method for PacBio sequencing (LILAP) using Tn5-based tagmentation and DNA circularization within one tube. We test LILAP with two Drosophila melanogaster individuals, and generate near-complete genomes, surpassing preexisting single-fly genomes. By analyzing variations in these two genomes, we characterize mutational processes: complex transpositions (transposon insertions together with extra duplications and/or deletions) prefer regions characterized by non-B DNA structures, and gene conversion of transposons occurs on both DNA and RNA levels. Concurrently, we generate two complete assemblies for the endosymbiotic bacterium Wolbachia in these flies and similarly detect transposon conversion. Thus, LILAP promises a broad PacBio sequencing adoption for not only mutational studies of flies and their symbionts but also explorations of other small organisms or precious samples.
Assuntos
Elementos de DNA Transponíveis , Drosophila melanogaster , Genoma de Inseto , Mutação , Wolbachia , Animais , Drosophila melanogaster/genética , Elementos de DNA Transponíveis/genética , Wolbachia/genética , Genoma de Inseto/genética , Sequenciamento de Nucleotídeos em Larga Escala/métodos , Análise de Sequência de DNA/métodos , Genômica/métodos , Conversão GênicaRESUMO
It is commonly thought that the long-term advantage of meiotic recombination is to dissipate genetic linkage, allowing natural selection to act independently on different loci. It is thus theoretically expected that genes with higher recombination rates evolve under more effective selection. On the other hand, recombination is often associated with GC-biased gene conversion (gBGC), which theoretically interferes with selection by promoting the fixation of deleterious GC alleles. To test these predictions, several studies assessed whether selection was more effective in highly recombining genes (due to dissipation of genetic linkage) or less effective (due to gBGC), assuming a fixed distribution of fitness effects (DFE) for all genes. In this study, I directly derive the DFE from a gene's evolutionary history (shaped by mutation, selection, drift, and gBGC) under empirical fitness landscapes. I show that genes that have experienced high levels of gBGC are less fit and thus have more opportunities for beneficial mutations. Only a small decrease in the genome-wide intensity of gBGC leads to the fixation of these beneficial mutations, particularly in highly recombining genes. This results in increased positive selection in highly recombining genes that is not caused by more effective selection. Additionally, I show that the death of a recombination hotspot can lead to a higher dN/dS than its birth, but with substitution patterns biased towards AT, and only at selected positions. This shows that controlling for a substitution bias towards GC is therefore not sufficient to rule out the contribution of gBGC to signatures of accelerated evolution. Finally, although gBGC does not affect the fixation probability of GC-conservative mutations, I show that by altering the DFE, gBGC can also significantly affect nonsynonymous GC-conservative substitution patterns.
Assuntos
Evolução Molecular , Conversão Gênica , Modelos Genéticos , Recombinação Genética , Seleção Genética , Aptidão Genética , Mutação , Composição de Bases , Ligação GenéticaRESUMO
Although gene conversion (GC) in Saccharomyces cerevisiae is the most error-free way to repair double-strand breaks (DSBs), the mutation rate during homologous recombination is 1,000 times greater than during replication. Many mutations involve dissociating a partially copied strand from its repair template and re-aligning with the same or another template, leading to -1 frameshifts in homonucleotide runs, quasipalindrome (QP)-associated mutations and microhomology-mediated interchromosomal template switches. We studied GC induced by HO endonuclease cleavage at MATα, repaired by an HMR::KI-URA3 donor. We inserted into HMR::KI-URA3 an 18-bp inverted repeat where one arm had a 4-bp insertion. Most GCs yield MAT::KI-ura3::QP + 4 (Ura-) outcomes, but template-switching produces Ura+ colonies, losing the 4-bp insertion. If the QP arm without the insertion is first encountered by repair DNA polymerase and is then (mis)used as a template, the palindrome is perfected. When the QP + 4 arm is encountered first, Ura+ derivatives only occur after second-end capture and second-strand synthesis. QP + 4 mutations are suppressed by mismatch repair (MMR) proteins Msh2, Msh3, and Mlh1, but not Msh6. Deleting Rdh54 significantly reduces QP mutations only when events creating Ura+ occur in the context of a D-loop but not during second-strand synthesis. A similar bias is found with a proofreading-defective DNA polymerase mutation (poI3-01). DSB-induced mutations differed in several genetic requirements from spontaneous events. We also created a + 1 frameshift in the donor, expanding a run of 4 Cs to 5 Cs. Again, Ura+ recombinants markedly increased by disabling MMR, suggesting that MMR acts during GC but favors the unbroken, template strand.
Assuntos
Quebras de DNA de Cadeia Dupla , Reparo de Erro de Pareamento de DNA , Mutação da Fase de Leitura , Mutagênese , Proteínas de Saccharomyces cerevisiae , Saccharomyces cerevisiae , Proteínas de Saccharomyces cerevisiae/genética , Proteínas de Saccharomyces cerevisiae/metabolismo , Saccharomyces cerevisiae/genética , Saccharomyces cerevisiae/metabolismo , Conversão Gênica , Proteínas de Ligação a DNA/genética , Proteínas de Ligação a DNA/metabolismo , Proteína 2 Homóloga a MutS/genética , Proteína 2 Homóloga a MutS/metabolismo , Proteína 3 Homóloga a MutS/genética , Proteína 3 Homóloga a MutS/metabolismo , Proteína 1 Homóloga a MutLRESUMO
Molecular dissection of meiotic recombination in mammals, combined with population-genetic and comparative studies, have revealed a complex evolutionary dynamic characterized by short-lived recombination hotspots. Hotspots are chromosome positions containing DNA sequences where the protein PRDM9 can bind and cause crossing-over. To explain these fast evolutionary dynamic, a so-called intra-genomic Red Queen model has been proposed, based on the interplay between two antagonistic forces: biased gene conversion, mediated by double-strand breaks, resulting in hotspot extinction (the hotspot conversion paradox), followed by positive selection favoring mutant PRDM9 alleles recognizing new sequence motifs. Although this model predicts many empirical observations, the exact causes of the positive selection acting on new PRDM9 alleles is still not well understood. In this direction, experiment on mouse hybrids have suggested that, in addition to targeting double strand breaks, PRDM9 has another role during meiosis. Specifically, PRDM9 symmetric binding (simultaneous binding at the same site on both homologues) would facilitate homology search and, as a result, the pairing of the homologues. Although discovered in hybrids, this second function of PRDM9 could also be involved in the evolutionary dynamic observed within populations. To address this point, here, we present a theoretical model of the evolutionary dynamic of meiotic recombination integrating current knowledge about the molecular function of PRDM9. Our modeling work gives important insights into the selective forces driving the turnover of recombination hotspots. Specifically, the reduced symmetrical binding of PRDM9 caused by the loss of high affinity binding sites induces a net positive selection eliciting new PRDM9 alleles recognizing new targets. The model also offers new insights about the influence of the gene dosage of PRDM9, which can paradoxically result in negative selection on new PRDM9 alleles entering the population, driving their eviction and thus reducing standing variation at this locus.
Assuntos
Evolução Molecular , Histona-Lisina N-Metiltransferase , Meiose , Histona-Lisina N-Metiltransferase/genética , Histona-Lisina N-Metiltransferase/metabolismo , Meiose/genética , Animais , Camundongos , Conversão Gênica , Quebras de DNA de Cadeia Dupla , Alelos , Modelos Genéticos , Humanos , Recombinação GenéticaRESUMO
In many mammals, recombination events are concentrated in hotspots directed by a sequence-specific DNA-binding protein named PRDM9. Intriguingly, PRDM9 has been lost several times in vertebrates, and notably among mammals, it has been pseudogenized in the ancestor of canids. In the absence of PRDM9, recombination hotspots tend to occur in promoter-like features such as CpG islands. It has thus been proposed that one role of PRDM9 could be to direct recombination away from PRDM9-independent hotspots. However, the ability of PRDM9 to direct recombination hotspots has been assessed in only a handful of species, and a clear picture of how much recombination occurs outside of PRDM9-directed hotspots in mammals is still lacking. In this study, we derived an estimator of past recombination activity based on signatures of GC-biased gene conversion in substitution patterns. We quantified recombination activity in PRDM9-independent hotspots in 52 species of boreoeutherian mammals. We observe a wide range of recombination rates at these loci: several species (such as mice, humans, some felids, or cetaceans) show a deficit of recombination, while a majority of mammals display a clear peak of recombination. Our results demonstrate that PRDM9-directed and PRDM9-independent hotspots can coexist in mammals and that their coexistence appears to be the rule rather than the exception. Additionally, we show that the location of PRDM9-independent hotspots is relatively more stable than that of PRDM9-directed hotspots, but that PRDM9-independent hotspots nevertheless evolve slowly in concert with DNA hypomethylation.
Assuntos
Histona-Lisina N-Metiltransferase , Recombinação Genética , Animais , Histona-Lisina N-Metiltransferase/genética , Histona-Lisina N-Metiltransferase/metabolismo , Humanos , Mamíferos/genética , Ilhas de CpG/genética , Eutérios/genética , Camundongos , Feminino , Conversão Gênica , Evolução MolecularRESUMO
The rate at which recombination events occur in a population is an indicator of its effective population size and the organism's reproduction mode. It determines the extent of linkage disequilibrium along the genome and, thereby, the efficacy of both purifying and positive selection. The population recombination rate can be inferred using models of genome evolution in populations. Classic methods based on the patterns of linkage disequilibrium provide the most accurate estimates, providing large sample sizes are used and the demography of the population is properly accounted for. Here, the capacity of approaches based on the sequentially Markov coalescent (SMC) to infer the genome-average recombination rate from as little as a single diploid genome is examined. SMC approaches provide highly accurate estimates even in the presence of changing population sizes, providing that (1) within genome heterogeneity is accounted for and (2) classic maximum-likelihood optimization algorithms are employed to fit the model. SMC-based estimates proved sensitive to gene conversion, leading to an overestimation of the recombination rate if conversion events are frequent. Conversely, methods based on the correlation of heterozygosity succeed in disentangling the rate of crossing over from that of gene conversion events, but only when the population size is constant and the recombination landscape homogeneous. These results call for a convergence of these two methods to obtain accurate and comparable estimates of recombination rates between populations.
Assuntos
Desequilíbrio de Ligação , Cadeias de Markov , Modelos Genéticos , Recombinação Genética , Genoma , Algoritmos , Genética Populacional/métodos , Conversão Gênica , Animais , Humanos , Densidade DemográficaRESUMO
Different frequencies amongst codons that encode the same amino acid (i.e. synonymous codons) have been observed in multiple species. Studies focused on uncovering the forces that drive such codon usage showed that a combined effect of mutational biases and translational selection works to produce different frequencies of synonymous codons. However, only few have been able to measure and distinguish between these forces that may leave similar traces on the coding regions. Here, we have developed a codon model that allows the disentangling of mutation, selection on amino acids and synonymous codons, and GC-biased gene conversion (gBGC) which we employed on an extensive dataset of 415 chordates and 191 arthropods. We found that chordates need 15 more synonymous codon categories than arthropods to explain the empirical codon frequencies, which suggests that the extent of codon usage can vary greatly between animal phyla. Moreover, methylation at CpG sites seems to partially explain these patterns of codon usage in chordates but not in arthropods. Despite the differences between the two phyla, our findings demonstrate that in both, GC-rich codons are disfavored when mutations are GC-biased, and the opposite is true when mutations are AT-biased. This indicates that selection on the genomic coding regions might act primarily to stabilize its GC/AT content on a genome-wide level. Our study shows that the degree of synonymous codon usage varies considerably among animals, but is likely governed by a common underlying dynamic.
Assuntos
Artrópodes , Uso do Códon , Seleção Genética , Animais , Artrópodes/genética , Cordados/genética , Mutação , Evolução Molecular , Códon , Modelos Genéticos , Composição de Bases , Conversão GênicaRESUMO
We present a method for efficiently identifying clusters of identical-by-descent haplotypes in biobank-scale sequence data. Our multi-individual approach enables much more computationally efficient inference of identity by descent (IBD) than approaches that infer pairwise IBD segments and provides locus-specific IBD clusters rather than IBD segments. Our method's computation time, memory requirements, and output size scale linearly with the number of individuals in the dataset. We also present a method for using multi-individual IBD to detect alleles changed by gene conversion. Application of our methods to the autosomal sequence data for 125,361 White British individuals in the UK Biobank detects more than 9 million converted alleles. This is 2,900 times more alleles changed by gene conversion than were detected in a previous analysis of familial data. We estimate that more than 250,000 sequenced probands and a much larger number of additional genomes from multi-generational family members would be required to find a similar number of alleles changed by gene conversion using a family-based approach. Our IBD clustering method is implemented in the open-source ibd-cluster software package.
Assuntos
Bancos de Espécimes Biológicos , Conversão Gênica , Humanos , Software , Haplótipos/genética , Cromossomos , Polimorfismo de Nucleotídeo ÚnicoRESUMO
While the malaria parasite Plasmodium falciparum has low average genome-wide diversity levels, likely due to its recent introduction from a gorilla-infecting ancestor (approximately 10,000 to 50,000 years ago), some genes display extremely high diversity levels. In particular, certain proteins expressed on the surface of human red blood cell-infecting merozoites (merozoite surface proteins (MSPs)) possess exactly 2 deeply diverged lineages that have seemingly not recombined. While of considerable interest, the evolutionary origin of this phenomenon remains unknown. In this study, we analysed the genetic diversity of 2 of the most variable MSPs, DBLMSP and DBLMSP2, which are paralogs (descended from an ancestral duplication). Despite thousands of available Illumina WGS datasets from malaria-endemic countries, diversity in these genes has been hard to characterise as reads containing highly diverged alleles completely fail to align to the reference genome. To solve this, we developed a pipeline leveraging genome graphs, enabling us to genotype them at high accuracy and completeness. Using our newly- resolved sequences, we found that both genes exhibit 2 deeply diverged lineages in a specific protein domain (DBL) and that one of the 2 lineages is shared across the genes. We identified clear evidence of nonallelic gene conversion between the 2 genes as the likely mechanism behind sharing, leading us to propose that gene conversion between diverged paralogs, and not recombination suppression, can generate this surprising genealogy; a model that is furthermore consistent with high diversity levels in these 2 genes despite the strong historical P. falciparum transmission bottleneck.
Assuntos
Hominidae , Malária Falciparum , Malária , Parasitos , Animais , Humanos , Plasmodium falciparum/metabolismo , Parasitos/metabolismo , Conversão Gênica , Antígenos de Superfície , Malária/parasitologia , Proteínas de Protozoários/genética , Proteínas de Protozoários/metabolismo , Variação GenéticaRESUMO
Artificial biomolecular condensates are emerging as a versatile approach to organize molecular targets and reactions without the need for lipid membranes. Here we ask whether the temporal response of artificial condensates can be controlled via designed chemical reactions. We address this general question by considering a model problem in which a phase separating component participates in reactions that dynamically activate or deactivate its ability to self-attract. Through a theoretical model we illustrate the transient and equilibrium effects of reactions, linking condensate response and reaction parameters. We experimentally realize our model problem using star-shaped DNA motifs known as nanostars to generate condensates, and we take advantage of strand invasion and displacement reactions to kinetically control the capacity of nanostars to interact. We demonstrate reversible dissolution and growth of DNA condensates in the presence of specific DNA inputs, and we characterize the role of toehold domains, nanostar size, and nanostar valency. Our results will support the development of artificial biomolecular condensates that can adapt to environmental changes with prescribed temporal dynamics.
Assuntos
Condensados Biomoleculares , Empacotamento do DNA , Replicação do DNA , Conversão Gênica , Motivos de NucleotídeosRESUMO
Identification of novel HLA-A*23:128 allele generated by interlocus gene conversion in Brazilian bone marrow donor.
Assuntos
Medula Óssea , Conversão Gênica , Humanos , Brasil , Alelos , Doadores de Tecidos , Antígenos HLA-A/genéticaRESUMO
Population genetic inference of selection on the nucleotide sequence level often proceeds by comparison to a reference sequence evolving only under mutation and population demography. Among the few candidates for such a reference sequence is the 5' part of short introns (5SI) in Drosophila. In addition to mutation and population demography, however, there is evidence for a weak force favouring GC bases, likely due to GC-biased gene conversion (gBGC), and for the effect of linked selection. Here, we use polymorphism and divergence data of Drosophila melanogaster to detect and describe the forces affecting the evolution of the 5SI. We separately analyse mutation classes, compare them between chromosomes, and relate them to recombination rate frequencies. GC-conservative mutations seem to be mainly influenced by mutation and drift, with linked selection mostly causing differences between the central and the peripheral (i.e., telomeric and centromeric) regions of the chromosome arms. Comparing GC-conservative mutation patterns between autosomes and the X chromosome showed differences in mutation rates, rather than linked selection, in the central chromosomal regions after accounting for differences in effective population sizes. On the other hand, GC-changing mutations show asymmetric site frequency spectra, indicating the presence of gBGC, varying among mutation classes and in intensity along chromosomes, but approximately equal in strength in autosomes and the X chromosome.
Assuntos
Drosophila melanogaster , Conversão Gênica , Animais , Drosophila melanogaster/genética , Íntrons , Evolução Molecular , Mutação , Drosophila/genética , Cromossomo X/genética , Seleção GenéticaRESUMO
BACKGROUND: Despite GC variation constitutes a fundamental element of genome and species diversity, the precise mechanisms driving it remain unclear. The abundant sequence data available for the ITS2, a commonly employed phylogenetic marker in plants, offers an exceptional resource for exploring the GC variation across angiosperms. RESULTS: A comprehensive selection of 8666 species, comprising 165 genera, 63 families, and 30 orders were used for the analyses. The alignment of ITS2 sequence-structures and partitioning of secondary structures into paired and unpaired regions were performed using 4SALE. Substitution rates and frequencies among GC base-pairs in the paired regions of ITS2 were calculated using RNA-specific models in the PHASE package. The results showed that the distribution of ITS2 GC contents on the angiosperm phylogeny was heterogeneous, but their increase was generally associated with ITS2 sequence homogenization, thereby supporting the occurrence of GC-biased gene conversion (gBGC) during the concerted evolution of ITS2. Additionally, the GC content in the paired regions of the ITS2 secondary structure was significantly higher than that of the unpaired regions, indicating the selection of GC for thermodynamic stability. Furthermore, the RNA substitution models demonstrated that base-pair transformations favored both the elevation and fixation of GC in the paired regions, providing further support for gBGC. CONCLUSIONS: Our findings highlight the significance of secondary structure in GC investigation, which demonstrate that both gBGC and structure-based selection are influential factors driving angiosperm ITS2 GC content.
Assuntos
Magnoliopsida , Humanos , Magnoliopsida/genética , Filogenia , Conversão Gênica , Composição de Bases , RNA , Evolução MolecularRESUMO
Ultraconserved elements (UCEs) are the most conserved regions among the genomes of evolutionarily distant species and are thought to play critical biological functions. However, some UCEs rapidly evolved in specific lineages, and whether they contributed to adaptive evolution is still controversial. Here, using an increased number of sequenced genomes with high taxonomic coverage, we identified 2191 mammalian UCEs and 5938 avian UCEs from 95 mammal and 94 bird genomes, respectively. Our results show that these UCEs are functionally constrained and that their adjacent genes are prone to widespread expression with low expression diversity across tissues. Functional enrichment of mammalian and avian UCEs shows different trends indicating that UCEs may contribute to adaptive evolution of taxa. Focusing on lineage-specific accelerated evolution, we discover that the proportion of fast-evolving UCEs in nine mammalian and 10 avian test lineages range from 0.19% to 13.2%. Notably, up to 62.1% of fast-evolving UCEs in test lineages are much more likely to result from GC-biased gene conversion (gBGC). A single cervid-specific gBGC region embracing the uc.359 allele significantly alters the expression of Nova1 and other neural-related genes in the rat brain. Combined with the altered regulatory activity of ancient gBGC-induced fast-evolving UCEs in eutherians, our results provide evidence that synergy between gBGC and selection shaped lineage-specific substitution patterns, even in the most constrained regulatory elements. In summary, our results show that gBGC played an important role in facilitating lineage-specific accelerated evolution of UCEs, and further support the idea that a combination of multiple evolutionary forces shapes adaptive evolution.
Assuntos
Conversão Gênica , Mamíferos , Animais , Ratos , Mamíferos/genética , Alelos , Aves/genética , Evolução Molecular , Antígeno Neuro-Oncológico VentralRESUMO
In the absence of recombination, the number of transposable elements (TEs) increases due to less efficient selection, but the dynamics of such TE accumulations are not well characterized. Leveraging a dataset of 21 independent events of recombination cessation of different ages in mating-type chromosomes of Microbotryum fungi, we show that TEs rapidly accumulated in regions lacking recombination, but that TE content reached a plateau at ca. 50% of occupied base pairs by 1.5 million years following recombination suppression. The same TE superfamilies have expanded in independently evolved non-recombining regions, in particular rolling-circle replication elements (Helitrons). Long-terminal repeat (LTR) retrotransposons of the Copia and Ty3 superfamilies also expanded, through transposition bursts (distinguished from gene conversion based on LTR divergence), with both non-recombining regions and autosomes affected, suggesting that non-recombining regions constitute TE reservoirs. This study improves our knowledge of genome evolution by showing that TEs can accumulate through bursts, following non-linear decelerating dynamics.