Frequent Gene Conversion Research Articles

Gene conversions are frequent but not under positive selection in the Siglec gene families of primates.

Siglecs are cell surface proteins that belong to the immunoglobulin superfamily and which bind sialic acids. They are composed of two groups, the conserved Siglecs and the CD33-related Siglecs. Previous studies have reported the occurrence of gene conversions between human CD33-related Siglecs and suggested that these conversions are adaptive because they increase the diversity of these immunoglobulin-related genes. Here, we analyze the Siglec genes of five primate species and show that gene conversions are not observed between conserved Siglec genes but that they are frequent between primate CD33-related Siglecs. The gene conversions between CD33-related Siglec genes only occur between similar genes and equally frequently in sialic acid binding and nonbinding domains. Furthermore, dN/dS ratio tests show that most of the Ig-like V-type 1 and the Ig-like C2-type 1 domains of Siglec genes evolve either neutrally or under purifying selection and that gene conversions were not responsible for the positively selected regions detected in the Ig-like V-type1 domain of the human SIGLEC7 and SIGLEC9 genes. Our results suggest that the frequent gene conversions between CD33-related Siglec genes are simply a consequence of the high degree of sequence similarity of these genes and that they are not adaptive.

Allele-mining of rice blast resistance genes at AC134922 locus

The AC134922 locus is one of the most rapidly evolving nucleotide binding site-leucine-rich repeat (NBS-LRR) gene family in rice genome. Six rice blast resistance (R) genes have been cloned from this locus and other two resistance candidate genes, Pi34 and Pi47, are also mapped to this complex locus. Therefore, it seems that more functional R genes could be identified from this locus. In this study, we cloned 22 genes from 12 cultivars based on allele-mining strategy at this locus and identified 6 rice blast R genes with 4 of them recognizing more than one isolates. Our result suggests that gene stacking might be the evolutionary strategy for complex gene locus to interact with rapidly evolving pathogens, which might provide a potential way for the cloning of durable resistance genes. Moreover, the mosaic structure and ambiguous ortholog/paralog relationships of these homologous genes, caused by frequent recombination and gene conversion, indicate that multiple alleles of this complex locus may serve as a reservoir for the evolutionary novelty of these R genes.

Purifying selection against gene conversions in the folate receptor genes of primates

We characterized the gene conversions between the human folate receptor (FOLR) genes and those of five other primate species. We found 26 gene conversions having an average length of 534 nucleotides. The length of these conversions is correlated with sequence similarity, converted regions have a higher GC-content and the average size of converted regions from a functional donor to another functional donor is significantly smaller than the average size from a functional donor to a pseudogene. Furthermore, the few conversions observed in the FOLR1 and FOLR2 genes did not change any amino acids in their coding regions and did not affect their promoter regions. In contrast, the promoter and coding regions of the FOLR3 gene are frequently converted and these conversions changed many amino acids in marmoset. These results suggest that purifying selection is limiting the functional impact that frequent gene conversions have on functional folate receptor genes.

Gene conversions are under purifying selection in the carcinoembryonic antigen immunoglobulin gene families of primates

The carcinoembryonic antigen (CEA) family contains a large number of glycoproteins belonging to the immunoglobulin superfamily. Here, we investigate whether the gene conversions occurring between primate CEA-related genes are adaptive. Our results show that primate CEA-related genes are subject to frequent and repeated gene conversion events. Furthermore, gene conversions occur most frequently between nearby genes sharing similar sequences, are not more frequent in Ig-like V-type 1 domains than in the Ig-like C2-type 1 domains and dN/dS ratio tests shown that both these domains evolve either neutrally or under purifying selection. Our results therefore suggest that CEA-related genes evolve under purifying selection and the frequent gene conversion events we observed likely represent selectively neutral events between genes having similar sequences and functions.

Structural variation and dynamics of the nuclear ribosomal intergenic spacer region in key members of theGibberella fujikuroispecies complex

The intergenic spacer (IGS) region of the ribosomal DNA was cloned and sequenced in eight species within the Gibberella fujikuroi species complex with anamorphs in the genus Fusarium, a group that includes the most relevant toxigenic species. DNA sequence analyses revealed two categories of repeated elements: long repeats and short repeats of 125 and 8 bp, respectively. Long repeats were present in two copies and were conserved in all the species analyzed, whereas different numbers of short repeat elements were observed, leading to species-specific IGS sequences with different length. In Fusarium subglutinans and Fusarium nygamai, these differences seemed to be the result of duplication and deletion events. Here, we propose a model based on unequal crossing over that can explain these processes. The partial IGS sequence of 22 Fusarium proliferatum isolates was also obtained to study variation at the intraspecific level. The results revealed no differences in terms of number or pattern of repeated elements and detected frequent gene conversion events. These results suggest that the homogenization observed at the intraspecific level might not be achieved primarily by unequal crossing-over events but rather by processes associated with recombination such as gene conversion events.

Molecular evolution of the neurohypophysial hormone precursors in mammals: Comparative genomics reveals novel mammalian oxytocin and vasopressin analogues

Among vertebrates the neurohypophysial hormones show considerable variation. However, in eutherian mammals they have been considered rather conserved, with arginine vasopressin (AVP) and oxytocin (OT) in all species except pig and some relatives, where lysine vasopressin replaces AVP. The availability of genomic data for a wide range of mammals makes it possible to assess whether these peptides and their precursors may be more variable in Eutheria than previously suspected. A survey of these data confirms that AVP and OT occur in most eutherians, but with exceptions. In a New-World monkey (marmoset, Callithrix jacchus) and in tree shrew (Tupaia belangeri), Pro8OT replaces OT, confirming a recent report for these species. In armadillo (Dasypus novemcinctus) Leu3OT replaces OT, while in tenrec (Echinops telfairi) Thr4AVP replaces AVP. In these two species there is also evidence for additional genes/pseudogenes, encoding much-modified forms of AVP, but in most other eutherian species there is no evidence for additional neurohypophysial hormone genes. Evolutionary analysis shows that sequences of eutherian neurohypophysial hormone precursors are generally strongly conserved, particularly those regions encoding active peptide and neurophysin. The close association between OT and VP genes has led to frequent gene conversion of sequences encoding neurophysins. A monotreme, platypus (Ornithorhynchus anatinus) has genes for OT and AVP, organized tail-to-tail as in eutherians, but in marsupials 3–4 genes are present for neurohypophysial hormones, organized tail-to-head as in lower vertebrates.

Balancing Selection at the Tomato RCR3 Guardee Gene Family Maintains Variation in Strength of Pathogen Defense

Coevolution between hosts and pathogens is thought to occur between interacting molecules of both species. This results in the maintenance of genetic diversity at pathogen antigens (or so-called effectors) and host resistance genes such as the major histocompatibility complex (MHC) in mammals or resistance (R) genes in plants. In plant–pathogen interactions, the current paradigm posits that a specific defense response is activated upon recognition of pathogen effectors via interaction with their corresponding R proteins. According to the “Guard-Hypothesis,” R proteins (the “guards”) can sense modification of target molecules in the host (the “guardees”) by pathogen effectors and subsequently trigger the defense response. Multiple studies have reported high genetic diversity at R genes maintained by balancing selection. In contrast, little is known about the evolutionary mechanisms shaping the guardee, which may be subject to contrasting evolutionary forces. Here we show that the evolution of the guardee RCR3 is characterized by gene duplication, frequent gene conversion, and balancing selection in the wild tomato species Solanum peruvianum. Investigating the functional characteristics of 54 natural variants through in vitro and in planta assays, we detected differences in recognition of the pathogen effector through interaction with the guardee, as well as substantial variation in the strength of the defense response. This variation is maintained by balancing selection at each copy of the RCR3 gene. Our analyses pinpoint three amino acid polymorphisms with key functional consequences for the coevolution between the guardee (RCR3) and its guard (Cf-2). We conclude that, in addition to coevolution at the “guardee-effector” interface for pathogen recognition, natural selection acts on the “guard-guardee” interface. Guardee evolution may be governed by a counterbalance between improved activation in the presence and prevention of auto-immune responses in the absence of the corresponding pathogen.

Genomic structure and evolution of multigene families: "flowers" on the human genome.

We report the results of an extensive investigation of genomic structures in the human genome, with a particular focus on relatively large repeats (>50 kb) in adjacent chromosomal regions. We named such structures “Flowers” because the pattern observed on dot plots resembles a flower. We detected a total of 291 Flowers in the human genome. They were predominantly located in euchromatic regions. Flowers are gene-rich compared to the average gene density of the genome. Genes involved in systems receiving environmental information, such as immunity and detoxification, were overrepresented in Flowers. Within a Flower, the mean number of duplication units was approximately four. The maximum and minimum identities between homologs in a Flower showed different distributions; the maximum identity was often concentrated to 100% identity, while the minimum identity was evenly distributed in the range of 78% to 100%. Using a gene conversion detection test, we found frequent and/or recent gene conversion events within the tested Flowers. Interestingly, many of those converted regions contained protein-coding genes. Computer simulation studies suggest that one role of such frequent gene conversions is the elongation of the life span of gene families in a Flower by the resurrection of pseudogenes.

Independent evolution of functional MHC class II DRB genes in New World bat species

Genes of the major histocompatibility complex (MHC) play a pivotal role in the vertebrate immune system and are attractive markers for functional, fitness-related, genetic variation. Although bats (Chiroptera) represent the second largest mammalian order and are prone to various emerging infectious diseases, little is known about MHC evolution in bats. In the present study, we examined expressed MHC class II DRB sequences (exons 1 to 4) of New World bat species, Saccopteryx bilineata, Carollia perspicillata, Noctilio albiventris and Noctilio leporinus (only exon 2). We found a wide range of copy number variation of DRB loci with one locus detected in the genus Noctilio and up to ten functional loci observed in S. bilineata. Sequence variation between alleles of the same taxa was high with evidence for positive selection. We found statistical support for recombination or gene conversion events among sequences within the same but not between bat species. Phylogenetic relationships among DRB alleles provided strong evidence for independent evolution of the functional MHC class II DRB genes in the three investigated species, either by recent gene duplication, or homogenization of duplicated loci by frequent gene conversion events. Phylogenetic analysis of all available chiropteran DRB exon 2 sequences confirmed their monophyletic origin within families, but revealed a possible trans-species mode of evolution pattern in congeneric bat species, e.g. within the genera Noctilio and Myotis. This is the first study investigating phylogenetic relationships of MHC genes within bats and therefore contributes to a better understanding of MHC evolution in one of the most dominant mammalian order.

Gene conversions in the growth hormone gene family of primates: Stronger homogenizing effects in the Hominidae lineage

In humans, the growth hormone/chorionic somatomammotropin gene family is composed of five highly similar genes. We characterized the gene conversions that occurred between the growth hormone genes of 11 primate species. We detected 48 conversions using GENECONV and others were only detected using phylogenetic analyses. Gene conversions were detected in all species analyzed, their average size (± standard deviation) is 197.8 ± 230.4 nucleotides, the size of the conversions is correlated with sequence similarity and converted regions are significantly more GC-rich than non-converted regions. Gene conversions have a stronger homogenizing effect in Hominidae genes than in other primate species. They are also less frequent in conserved gene regions and towards functionally important genes. This suggests that the high degree of sequence similarity observed between the growth hormone genes of primate species is a consequence of frequent gene conversions in gene regions which are under little selective constraints.

Extensive Gene Conversion Drives the Concerted Evolution of Paralogous Copies of the SRY Gene in European Rabbits

The human Y chromosome consists of ampliconic genes, which are located in palindromes and undergo frequent gene conversion, and single-copy genes including the primary sex-determining locus, SRY. Here, we demonstrate that SRY is duplicated in a large palindrome in the European rabbit (Oryctolagus cuniculus). Furthermore, we show through comparative sequencing that orthologous palindrome arms have diverged 0.40% between rabbit subspecies over at least 2 My, but paralogous palindrome arms have remained nearly identical. This provides clear evidence of gene conversion on the rabbit Y chromosome. Together with previous observations in humans, these results suggest that gene conversion is a general feature of the evolution of the mammalian Y chromosome.

Recombination-Associated Sequence Homogenization of Neighboring Alu Elements: Signature of Nonallelic Gene Conversion

Recently, researchers have begun to recognize that, in order to establish neutral models for disease association and evolutionary genomics studies, it is crucial to have a clear understanding of the genomic impact of nonallelic gene conversion. Drawing on previous successes in characterizing this phenomenon over protein-coding gene families, we undertook a computational analysis of neighboring Alu sequences in the genome scale. For this purpose, we developed adjusted comutation rate (aCMR), a novel statistical method measuring the excess number of identical point mutations shared by adjacent Alu sequences, vis-à-vis random pairs. Using aCMR, we uncovered a remarkable genome-wide sequence homogenization of neighboring Alus, with the strongest signal observed in the pseudoautosomal regions of the X and Y chromosomes. The magnitude of sequence homogenization between Alu pairs is greater with shorter interlocus distance, higher sequence identity, and parallel orientation. Moreover, shared substitutions show a strong directionality toward GC nucleotides, with multiple substitutions tending to cluster within the Alu sequence. Taken together, these observed recombination-associated sequence homogenization patterns are best explained by frequent ubiquitous gene conversion events between neighboring Alus. We believe that these observations help to illuminate the nature and impact of the enigmatic phenomenon of gene conversion.

Evolution of the sex-related locus and genomic features shared in microsporidia and fungi.

BackgroundMicrosporidia are obligate intracellular, eukaryotic pathogens that infect a wide range of animals from nematodes to humans, and in some cases, protists. The preponderance of evidence as to the origin of the microsporidia reveals a close relationship with the fungi, either within the kingdom or as a sister group to it. Recent phylogenetic studies and gene order analysis suggest that microsporidia share a particularly close evolutionary relationship with the zygomycetes.Methodology/Principal FindingsHere we expanded this analysis and also examined a putative sex-locus for variability between microsporidian populations. Whole genome inspection reveals a unique syntenic gene pair (RPS9-RPL21) present in the vast majority of fungi and the microsporidians but not in other eukaryotic lineages. Two other unique gene fusions (glutamyl-prolyl tRNA synthetase and ubiquitin-ribosomal subunit S30) that are present in metazoans, choanoflagellates, and filasterean opisthokonts are unfused in the fungi and microsporidians. One locus previously found to be conserved in many microsporidian genomes is similar to the sex locus of zygomycetes in gene order and architecture. Both sex-related and sex loci harbor TPT, HMG, and RNA helicase genes forming a syntenic gene cluster. We sequenced and analyzed the sex-related locus in 11 different Encephalitozoon cuniculi isolates and the sibling species E. intestinalis (3 isolates) and E. hellem (1 isolate). There was no evidence for an idiomorphic sex-related locus in this Encephalitozoon species sample. According to sequence-based phylogenetic analyses, the TPT and RNA helicase genes flanking the HMG genes are paralogous rather than orthologous between zygomycetes and microsporidians.Conclusion/SignificanceThe unique genomic hallmarks between microsporidia and fungi are independent of sequence based phylogenetic comparisons and further contribute to define the borders of the fungal kingdom and support the classification of microsporidia as unusual derived fungi. And the sex/sex-related loci appear to have been subject to frequent gene conversion and translocations in microsporidia and zygomycetes.

Frequent gene conversion events between the X and Y homologous chromosomal regions in primates

BackgroundMammalian sex-chromosomes originated from a pair of autosomes. A step-wise cessation of recombination is necessary for the proper maintenance of sex-determination and, consequently, generates a four strata structure on the X chromosome. Each stratum shows a specific per-site nucleotide sequence difference (p-distance) between the X and Y chromosomes, depending on the time of recombination arrest. Stratum 4 covers the distal half of the human X chromosome short arm and the p-distance of the stratum is ~10%, on average. However, a 100-kb region, which includes KALX and VCX, in the middle of stratum 4 shows a significantly lower p-distance (1-5%), suggesting frequent sequence exchanges or gene conversions between the X and Y chromosomes in humans. To examine the evolutionary mechanism for this low p-distance region, sequences of a corresponding region including KALX/Y from seven species of non-human primates were analyzed.ResultsPhylogenetic analysis of this low p-distance region in humans and non-human primate species revealed that gene conversion like events have taken place at least ten times after the divergence of New World monkeys and Catarrhini (i.e., Old World monkeys and hominoids). A KALY-converted KALX allele in white-handed gibbons also suggests a possible recent gene conversion between the X and Y chromosomes. In these primate sequences, the proximal boundary of this low p-distance region is located in a LINE element shared between the X and Y chromosomes, suggesting the involvement of this element in frequent gene conversions. Together with a palindrome on the Y chromosome, a segmental palindrome structure on the X chromosome at the distal boundary near VCX, in humans and chimpanzees, may mediate frequent sequence exchanges between X and Y chromosomes.ConclusionGene conversion events between the X and Y homologous regions have been suggested, mainly in humans. Here, we found frequent gene conversions in the evolutionary course of primates. An insertion of a LINE element at the proximal end of the region may be a cause for these frequent conversions. This gene conversion in humans may also be one of the genetic causes of Kallmann syndrome.

Mosaic gene conversion after a tandem duplication of mtDNA sequence in Diomedeidae (albatrosses)

Although the tandem duplication of mitochondrial (mt) sequences, especially those of the control region (CR), has been detected in metazoan species, few studies have focused on the features of the duplicated sequence itself, such as the gene conversion rate, distribution patterns of the variation, and relative rates of evolution between the copies. To investigate the features of duplicated mt sequences, we partially sequenced the mt genome of 16 Phoebastria albatrosses belonging to three species (P. albatrus, P. nigripes, and P. immutabilis). More than 2,300 base pairs of tandemly-duplicated sequence were shared by all three species. The observed gene arrangement was shared in the three Phoebastria albatrosses and suggests that the duplication event occurred in the common ancestor of the three species. Most of the copies in each individual were identical or nearly identical, and were maintained through frequent gene conversions. By contrast, portions of CR domains I and III had different phylogenetic signals, suggesting that gene conversion had not occurred in those sections after the speciation of the three species. Several lines of data, including the heterogeneity of the rate of molecular evolution, nucleotide differences, and putative secondary structures, suggests that the two sequences in CR domain I are maintained through selection; however, additional studies into the mechanisms of gene conversion and mtDNA synthesis are required to confirm this hypothesis.

Gene Conversion Maintains Nonfunctional Transposable Elements in an Obligate Mutualistic Endosymbiont

Long-term bacterial endosymbionts typically exhibit reduced genomes, lack genes encoding recombination functions and transposable elements, such as insertion sequences (ISs). In sharp contrast, I found that ISs constitute 2.4% of the genome of the obligate mutualistic endosymbiont Wolbachia wBm. Although no IS copy is transpositionally functional, I show that ISs persist in wBm because of frequent recombinational gene conversion (GC) homogenizing homologous IS sequences. These results not only indicate that there exists a functional recombination molecular machinery in wBm, but they also suggest that, by slowing down the rate of IS degradation and loss, GC may represent a major force influencing reductive evolution in wBm.

Gene Conversion and Positive Selection Driving the Evolution of the Caenorhabditis ssp. ZIM/HIM-8 Protein Family

In C. elegans, four C2H2 zinc-finger proteins (ZIM-1, ZIM-2, ZIM-3, and HIM-8), which are arranged in tandem, mediate chromosome-specific pairing and synapsis during meiosis. The zim/him-8 genes from three Caenorhabditis species were contrasted in an effort to investigate the mechanisms driving their evolution. Here it is shown that the preservation of higher degree of sequence similarity in the N-terminal portion, particularly in several regions within the second exon between paralogous zim genes (especially between zim-1 and zim-3), is due to independent interparalogue gene conversions. However, the evolutionary force is not uniformly strong across species. The present data reveal that more frequent gene conversion events have occurred in C. elegans, whereas only gene conversions between zim-1 and zim-3 are detected in C. remanei. Although gene conversions are predicted to be present among zim-1, zim-2, and zim-3 in C. briggsae, the conversion tracts between zim-1/zim-2 and zim-2/zim-3 are very short. Moreover, positive selection analysis was performed on the basis of the significantly discordant phylogenies reconstructed using the N- and C-terminal sequences, respectively. Several codon sites located in the regions that are supposed not to have experienced gene conversions are predicted to be under the influence of positive selection. In comparison, stronger positive selection has acted on the C-terminal region relative to the N-terminal region. Thus, the zim/him-8 genes that evolve concertedly have also been shown to undergo adaptive diversifying selection.

Molecular evolution of the keratin associated protein gene family in mammals, role in the evolution of mammalian hair

BackgroundHair is unique to mammals. Keratin associated proteins (KRTAPs), which contain two major groups: high/ultrahigh cysteine and high glycine-tyrosine, are one of the major components of hair and play essential roles in the formation of rigid and resistant hair shafts.ResultsThe KRTAP family was identified as being unique to mammals, and near-complete KRTAP gene repertoires for eight mammalian genomes were characterized in this study. An expanded KRTAP gene repertoire was found in rodents. Surprisingly, humans have a similar number of genes as other primates despite the relative hairlessness of humans. We identified several new subfamilies not previously reported in the high/ultrahigh cysteine KRTAP genes. Genes in many subfamilies of the high/ultrahigh cysteine KRTAP genes have evolved by concerted evolution with frequent gene conversion events, yielding a higher GC base content for these gene sequences. In contrast, the high glycine-tyrosine KRTAP genes have evolved more dynamically, with fewer gene conversion events and thus have a lower GC base content, possibly due to positive selection.ConclusionMost of the subfamilies emerged early in the evolution of mammals, thus we propose that the mammalian ancestor should have a diverse KRTAP gene repertoire. We propose that hair content characteristics have evolved and diverged rapidly among mammals because of rapid divergent evolution of KRTAPs between species. In contrast, subfamilies of KRTAP genes have been homogenized within each species due to concerted evolution.

Evidence that strong positive selection drives neofunctionalization in the tandemly duplicated polyhomeotic genes in Drosophila

The polyhomeotic (ph) locus in Drosophila melanogaster consists of the two tandemly duplicated genes ph-d (distal) and ph-p (proximal). They code for transcriptional repressors belonging to the Polycomb group proteins, which regulate homeotic genes and hundreds of other loci. Although the duplication of ph occurred at least 25 million to 30 million years ago, both copies are very similar to each other at both the DNA and the protein levels, probably because of the action of frequent gene conversion. Despite this homogenizing force, differential regulation of both transcriptional units suggests that the functions of the duplicates have begun to diverge. Here, we provide evidence that this functional divergence is driven by positive selection. Based on resequencing of an approximately 30-kb region around the ph locus in an African sample of D. melanogaster X chromosomes, we identified a selective sweep, estimated its age and the strength of selection, and mapped the target of selection to a narrow interval of the ph-p gene. This noncoding region contains a large intron with several regulatory elements that are absent in the ph-d duplicate. Our results suggest that neofunctionalization has been achieved in the Drosophila ph genes through the action of strong positive selection and the inactivation of gene conversion in part of the gene.

Adaptive Evolution and Frequent Gene Conversion in the Brain Expressed X-Linked Gene Family in Mammals

This work examines the molecular evolution of a brain-expressed X-linked gene family in the mammalian genomes of human, chimp, macaque, mouse, rat, dog, and cow. The gene structures are well conserved across family members and among the mammals in that all five members have three exons with the first two exons untranslated. Furthermore, the five members are arranged tandemly on chromosome X with Bex5, Bex1, Bex2 on the negative strand and Bex4, Bex3 on the positive strand, and this physical arrangement remains conserved among species. Sequence analyses indicate that gene conversion has been frequent and ongoing among Bex1-4, occurring in multiple species independently. All gene conversions in different species between Bex1 and Bex4, and between Bex2 and Bex3, appear to be limited to the upstream regions of the third exon, whereas the gene conversions occurred independently in different species between Bex1 and Bex2 and cover only the third exon. Bex5 appears to have little exchange of genetic information with other members, possibly due to its distance from other members. The GC content decreases from 5'-UTR, intron 1, intron 2, coding region, to 3'-UTR, reflecting faithfully the frequency of gene conversion in different regions of the Bex genes. Sequence analyses also suggest that both relaxed selective constraint and positive selection have acted on the Bex members after duplication. In particular, Bex3 shows strong evidence of positive selection and seems to have evolved a new gene function after gene duplication.