Including Gene Conversion Research Articles

Characterization of the mutation spectrum in a Pakistani cohort of type 3 von Willebrand disease.

Type 3 von Willebrand disease (VWD), a severe autosomal recessive hereditary bleeding disorder, is described by the virtual absence of von Willebrand factor (VWF). In consanguineous populations, for example Pakistan, the disease is reported with a higher incidence rate than the worldwide prevalence. This study aims to characterize molecular pathology and clinical profile of type 3 VWD cohort of Pakistani origin. In total, 48 patients were enrolled in the current study. Initially, the index patients (IPs) were evaluated by a standardized questionnaire for recording bleeding manifestations and by performing conventional coagulation tests. The diagnosis of VWD type 3 was confirmed by VWF antigens less than 5IU/dL. Direct sequencing of VWF gene (VWF) was carried out to identify causative gene variations. We evaluated the potential consequence of novel splice site and missense variations by predictive computational programs and in silico structural analysis. VWF mutations were detected in 46 out of 48 IPs (95.8%), predominantly as homozygous variants. In total, twenty-nine different gene defects were characterized in this cohort from which 10 (34.5%) are novel. The majority of the mutations were null alleles (66%; including gene conversions, nonsense, splice site variations, small deletions and insertions), and 34% of them were missense substitutions. Herein, we reported for the first time, the pattern of gene defects in Pakistani type 3 VWD cohort. We identified a wide heterogeneous mutation spectrum along with variability in the type of bleeding episodes.

Gene tree species tree reconciliation with gene conversion.

Gene tree/species tree reconciliation is a recent decisive progress in phylogenetic methods, accounting for the possible differences between gene histories and species histories. Reconciliation consists in explaining these differences by gene-scale events such as duplication, loss, transfer, which translates mathematically into a mapping between gene tree nodes and species tree nodes or branches. Gene conversion is a frequent and important evolutionary event, which results in the replacement of a gene by a copy of another from the same species and in the same gene tree. Including this event in reconciliation models has never been attempted because it introduces a dependency between lineages, and standard algorithms based on dynamic programming become ineffective. We propose here a novel mathematical framework including gene conversion as an evolutionary event in gene tree/species tree reconciliation. We describe a randomized algorithm that finds, in polynomial running time, a reconciliation minimizing the number of duplications, losses and conversions in the case when their weights are equal. We show that the space of optimal reconciliations includes an analog of the last common ancestor reconciliation, but is not limited to it. Our algorithm outputs any optimal reconciliation with a non-null probability. We argue that this study opens a research avenue on including gene conversion in reconciliation, and discuss its possible importance in biology.

Whole genome sequencing reveals within-host genetic changes in paired meningococcal carriage isolates from Ethiopia

BackgroundMeningococcal colonization is a prerequisite for transmission and disease, but the bacterium only very infrequently causes disease while asymptomatic carriage is common. Carriage is highly dynamic, showing a great variety across time and space within and across populations, but also within individuals. The understanding of genetic changes in the meningococcus during carriage, when the bacteria resides in its natural niche, is important for understanding not only the carriage state, but the dynamics of the entire meningococcal population.ResultsPaired meningococcal isolates, obtained from 50 asymptomatic carriers about 2 months apart were analyzed with whole genome sequencing (WGS). Phylogenetic analysis revealed that most paired isolates from the same individual were closely related, and the average and median number of allelic differences between paired isolates defined as the same strain was 35. About twice as many differences were seen between isolates from different individuals within the same sequence type (ST). In 8%, different strains were detected at different time points. A difference in ST was observed in 6%, including an individual who was found to carry three different STs over the course of 9 weeks. One individual carried different strains from the same ST.In total, 566 of 1605 cgMLST genes had undergone within-host genetic changes in one or more pairs. The most frequently changed cgMLST gene was relA that was changed in 47% of pairs. Across the whole genome, pilE, differed mostly, in 85% of the pairs. The most frequent mechanisms of genetic difference between paired isolates were phase variation and recombination, including gene conversion. Different STs showed variation with regard to which genes that were most frequently changed, mostly due to absence/presence of phase variation.ConclusionsThis study revealed within-host genetic differences in meningococcal isolates during short-term asymptomatic carriage. The most frequently changed genes were genes belonging to the pilin family, the restriction/modification system, opacity proteins and genes involved in glycosylation. Higher resolution genome-wide sequence typing is necessary to resolve the diversity of isolates and reveals genetic differences not discovered by traditional typing schemes, and would be the preferred choice of technology.

Understanding the contrasting spatial haplotype patterns of malaria-protective β-globin polymorphisms

The malaria-protective β-globin polymorphisms, sickle-cell (βS) and β0-thalassaemia, are canonical examples of human adaptation to infectious disease. Occurring on distinct genetic backgrounds, they vary markedly in their patterns of linked genetic variation at the population level, suggesting different evolutionary histories. βS is associated with five classical restriction fragment length polymorphism haplotypes that exhibit remarkable specificity in their geographical distributions; by contrast, β0-thalassaemia mutations are found on haplotypes whose distributions overlap considerably. Here, we explore why these two polymorphisms display contrasting spatial haplotypic distributions, despite having malaria as a common selective pressure. We present a meta-population genetic model, incorporating individual-based processes, which tracks the evolution of β-globin polymorphisms on different haplotypic backgrounds. Our simulations reveal that, depending on the rate of mutation, a large population size and/or high population growth rate are required for both the βS- and the β0-thalassaemia-like patterns. However, whilst the βS-like pattern is more likely when population subdivision is high, migration low and long-distance migration absent, the opposite is true for β0-thalassaemia. Including gene conversion has little effect on the overall probability of each pattern; however, when inter-haplotype fitness variation exists, gene conversion is more likely to have contributed to the diversity of haplotypes actually present in the population. Our findings highlight how the contrasting spatial haplotype patterns exhibited by βS and β0-thalassaemia may provide important indications as to the evolution of these adaptive alleles and the demographic history of the populations in which they have evolved.

Revisiting the diffusion approximation to estimate evolutionary rates of gene family diversification

Genetic diversity in multigene families is shaped by multiple processes, including gene conversion and point mutation. Because multi-gene families are involved in crucial traits of organisms, quantifying the rates of their genetic diversification is important. With increasing availability of genomic data, there is a growing need for quantitative approaches that integrate the molecular evolution of gene families with their higher-scale function. In this study, we integrate a stochastic simulation framework with population genetics theory, namely the diffusion approximation, to investigate the dynamics of genetic diversification in a gene family. Duplicated genes can diverge and encode new functions as a result of point mutation, and become more similar through gene conversion. To model the evolution of pairwise identity in a multigene family, we first consider all conversion and mutation events in a discrete manner, keeping track of their details and times of occurrence; second we consider only the infinitesimal effect of these processes on pairwise identity accounting for random sampling of genes and positions. The purely stochastic approach is closer to biological reality and is based on many explicit parameters, such as conversion tract length and family size, but is more challenging analytically. The population genetics approach is an approximation accounting implicitly for point mutation and gene conversion, only in terms of per-site average probabilities. Comparison of these two approaches across a range of parameter combinations reveals that they are not entirely equivalent, but that for certain relevant regimes they do match. As an application of this modelling framework, we consider the distribution of nucleotide identity among VSG genes of African trypanosomes, representing the most prominent example of a multi-gene family mediating parasite antigenic variation and within-host immune evasion.

Rapid Evolution of Enormous, Multichromosomal Genomes in Flowering Plant Mitochondria with Exceptionally High Mutation Rates

Author SummaryA fundamental challenge in evolutionary biology is to explain why organisms exhibit dramatic variation in genome size and complexity. One hypothesis predicts that high rates of mutation in DNA sequence create selection against large and complex genomes, which are more susceptible to mutational disruption. Species of flowering plants in the genus Silene vary by approximately 100-fold in the rates of mutation in their mitochondrial DNA, providing an excellent opportunity to test the predicted effects of high mutation rates on genome evolution. Contrary to expectation, Silene species with elevated mutation rates have experienced dramatic expansions in mitochondrial genome size compared to their slowly evolving relatives, resulting in the largest known mitochondrial genomes. In addition to the increases in size and mutation rate, these genomes also reveal a history of rapid change in genome structure. They have been fragmented into dozens of chromosomes and appear to have experienced major reductions in recombination activity. All of these changes have occurred in just the past few million years. This mitochondrial genome diversity within the genus Silene provides a striking example of rapid genomic change and raises new hypotheses regarding the relationship between mutation rate and genome evolution.

Instability of subtelomeric regions during meiosis in Magnaporthe oryzae

Dynamics of chromosomal ends during meiosis was examined using a tetrad F1 population derived from a cross between Setaria and Triticum isolates of Magnaporthe oryzae. Telomeric fragments were liberated by six restriction enzymes and detected with a telomere probe. Each fragment was assigned to one of the 14 chromosomal ends using chromosome-specific markers. Size shifts and non-Mendelian segregation (4:0, 3:1, 1:3, and 0:4) were frequently observed in these telomeric fragments and were considered to be caused by deletion, insertion, point mutation, and gene conversion. Similar results were obtained in another tetrad F1 population derived from a cross between Oryza and Triticum isolates. These results suggest that subtelomeric regions are unstable during meiosis and are prone to various rearrangements including gene conversion.

Evolution of the nuclear ribosomal DNA intergenic spacer in four species of the Daphnia pulex complex

BackgroundConcerted evolution refers to the pattern in which copies of multigene families show high intraspecific sequence homogeneity but high interspecific sequence diversity. Sequence homogeneity of these copies depends on relative rates of mutation and recombination, including gene conversion and unequal crossing over, between misaligned copies. The internally repetitive intergenic spacer (IGS) is located between the genes for the 28S and 18S ribosomal RNAs. To identify patterns of recombination and/or homogenization within IGS repeat arrays, and to identify regions of the IGS that are under functional constraint, we analyzed 13 complete IGS sequences from 10 individuals representing four species in the Daphnia pulex complex.ResultsGene conversion and unequal crossing over between misaligned IGS repeats generates variation in copy number between arrays, as has been observed in previous studies. Moreover, terminal repeats are rarely involved in these events. Despite the occurrence of recombination, orthologous repeats in different species are more similar to one another than are paralogous repeats within species that diverged less than 4 million years ago. Patterns consistent with concerted evolution of these repeats were observed between species that diverged 8-10 million years ago. Sequence homogeneity varies along the IGS; the most homogeneous regions are downstream of the 28S rRNA gene and in the region containing the core promoter. The inadvertent inclusion of interspecific hybrids in our analysis uncovered evidence of both inter- and intrachromosomal recombination in the nonrepetitive regions of the IGS.ConclusionsOur analysis of variation in ribosomal IGS from Daphnia shows that levels of homogeneity within and between species result from the interaction between rates of recombination and selective constraint. Consequently, different regions of the IGS are on substantially different evolutionary trajectories.

Homoeologous nonreciprocal recombination in polyploid cotton

Polyploid formation and processes that create partial genomic duplication generate redundant genomic information, whose fate is of particular interest to evolutionary biologists. Different processes can lead to diversification among duplicate genes, which may be counterbalanced by mechanisms that retard divergence, including gene conversion via nonreciprocal homoeologous exchange. Here, we used genomic resources in diploid and allopolyploid cotton (Gossypium) to detect homoeologous single nucleotide polymorphisms provided by expressed sequence tags from G. arboreum (A genome), G. raimondii (D genome) and G. hirsutum (AD genome), allowing us to identify homoeo-single nucleotide polymorphism patterns indicative of potential homoeologous exchanges. We estimated the proportion of contigs in G. hirsutum that have experienced nonreciprocal homoeologous exchanges since the origin of polyploid cotton 1-2 million years ago (Mya) to be between 1.8% and 1.9%. To address the question of when the intergenomic exchange occurred, we assayed six of the genes affected by homoeo-recombination in all five Gossypium allopolyploids using a phylogenetic approach. This analysis revealed that nonreciprocal homoeologous exchanges have occurred throughout polyploid divergence and speciation, as opposed to saltationally with polyploid formation. In addition, some genomic regions show multiple patterns of homoeologous recombination among species.

Rad51 suppresses gross chromosomal rearrangement at centromere in Schizosaccharomyces pombe

Centromere that plays a pivotal role in chromosome segregation is composed of repetitive elements in many eukaryotes. Although chromosomal regions containing repeats are the hotspots of rearrangements, little is known about the stability of centromere repeats. Here, by using a minichromosome that has a complete set of centromere sequences, we have developed a fission yeast system to detect gross chromosomal rearrangements (GCRs) that occur spontaneously. Southern and comprehensive genome hybridization analyses of rearranged chromosomes show two types of GCRs: translocation between homologous chromosomes and formation of isochromosomes in which a chromosome arm is replaced by a copy of the other. Remarkably, all the examined isochromosomes contain the breakpoint in centromere repeats, showing that isochromosomes are produced by centromere rearrangement. Mutations in the Rad3 checkpoint kinase increase both types of GCRs. In contrast, the deletion of Rad51 recombinase preferentially elevates isochromosome formation. Chromatin immunoprecipitation analysis shows that Rad51 localizes at centromere around S phase. These data suggest that Rad51 suppresses rearrangements of centromere repeats that result in isochromosome formation.

Analysis of the human Alu Ye lineage

BackgroundAlu elements are short (~300 bp) interspersed elements that amplify in primate genomes through a process termed retroposition. The expansion of these elements has had a significant impact on the structure and function of primate genomes. Approximately 10 % of the mass of the human genome is comprised of Alu elements, making them the most abundant short interspersed element (SINE) in our genome. The majority of Alu amplification occurred early in primate evolution, and the current rate of Alu retroposition is at least 100 fold slower than the peak of amplification that occurred 30–50 million years ago. Alu elements are therefore a rich source of inter- and intra-species primate genomic variation.ResultsA total of 153 Alu elements from the Ye subfamily were extracted from the draft sequence of the human genome. Analysis of these elements resulted in the discovery of two new Alu subfamilies, Ye4 and Ye6, complementing the previously described Ye5 subfamily. DNA sequence analysis of each of the Alu Ye subfamilies yielded average age estimates of ~14, ~13 and ~9.5 million years old for the Alu Ye4, Ye5 and Ye6 subfamilies, respectively. In addition, 120 Alu Ye4, Ye5 and Ye6 loci were screened using polymerase chain reaction (PCR) assays to determine their phylogenetic origin and levels of human genomic diversity.ConclusionThe Alu Ye lineage appears to have started amplifying relatively early in primate evolution and continued propagating at a low level as many of its members are found in a variety of hominoid (humans, greater and lesser ape) genomes. Detailed sequence analysis of several Alu pre-integration sites indicated that multiple types of events had occurred, including gene conversions, near-parallel independent insertions of different Alu elements and Alu-mediated genomic deletions. A potential hotspot for Alu insertion in the Fer1L3 gene on chromosome 10 was also identified.

Hypervariable and Highly Divergent Intron?Exon Organizations in the Chordate Oikopleura dioica

Oikopleura dioica is a pelagic tunicate with a very small genome and a very short life cycle. In order to investigate the intron-exon organizations in Oikopleura, we have isolated and characterized ribosomal protein EF-1alpha, Hox, and alpha-tubulin genes. Their intron positions have been compared with those of the same genes from various invertebrates and vertebrates, including four species with entirely sequenced genomes. Oikopleura genes, like Caenorhabditis genes, have introns at a large number of nonconserved positions, which must originate from late insertions or intron sliding of ancient insertions. Both species exhibit hypervariable intron-exon organization within their alpha-tubulin gene family. This is due to localization of most nonconserved intron positions in single members of this gene family. The hypervariability and divergence of intron positions in Oikopleura and Caenorhabditis may be related to the predominance of short introns, the processing of which is not very dependent upon the exonic environment compared to large introns. Also, both species have an undermethylated genome, and the control of methylation-induced point mutations imposes a control on exon size, at least in vertebrate genes. That introns placed at such variable positions in Oikopleura or C. elegans may serve a specific purpose is not easy to infer from our current knowledge and hypotheses on intron functions. We propose that new introns are retained in species with very short life cycles, because illegitimate exchanges including gene conversion are repressed. We also speculate that introns placed at gene-specific positions may contribute to suppressing these exchanges and thereby favor their own persistence.

The effect of gene conversion on the divergence between duplicated genes.

Nonindependent evolution of duplicated genes is called concerted evolution. In this article, we study the evolutionary process of duplicated regions that involves concerted evolution. The model incorporates mutation and gene conversion: the former increases d, the divergence between two duplicated regions, while the latter decreases d. It is demonstrated that the process consists of three phases. Phase I is the time until d reaches its equilibrium value, d(0). In phase II d fluctuates around d(0), and d increases again in phase III. Our simulation results demonstrate that the length of concerted evolution (i.e., phase II) is highly variable, while the lengths of the other two phases are relatively constant. It is also demonstrated that the length of phase II approximately follows an exponential distribution with mean tau, which is a function of many parameters including gene conversion rate and the length of gene conversion tract. On the basis of these findings, we obtain the probability distribution of the level of divergence between a pair of duplicated regions as a function of time, mutation rate, and tau. Finally, we discuss potential problems in genomic data analysis of duplicated genes when it is based on the molecular clock but concerted evolution is common.

The Rad52–Rad59 complex interacts with Rad51 and replication protein A

The RAD52 gene is essential for homology-dependent repair of double-strand breaks in Saccharomyces cerevisiae. Rad52 forms complexes with Rad51, replication protein A (RPA) or Rad59 and its presence is essential for the formation of Rad51–Rad52–Rad59 and RPA–Rad52–Rad59 complexes. The N-terminal region of Rad52, which is required for self-interaction to form a ring structure, is required for interaction with Rad59. Rad59 also shows self-interaction suggesting the formation of heteromeric and homomeric rings of Rad52 and Rad59. In wild-type cells, we propose the Rad51–Rad52–Rad59 complex is involved in conservative recombination events, including gene conversion and reciprocal recombination, whereas the Rad52–Rad59 complex participates in single-strand annealing.

Genetic requirements for RAD51- and RAD54-independent break-induced replication repair of a chromosomal double-strand break.

Broken chromosomes can be repaired by several homologous recombination mechanisms, including gene conversion and break-induced replication (BIR). In Saccharomyces cerevisiae, an HO endonuclease-induced double-strand break (DSB) is normally repaired by gene conversion. Previously, we have shown that in the absence of RAD52, repair is nearly absent and diploid cells lose the broken chromosome; however, in cells lacking RAD51, gene conversion is absent but cells can repair the DSB by BIR. We now report that gene conversion is also abolished when RAD54, RAD55, and RAD57 are deleted but BIR occurs, as with rad51Delta cells. DSB-induced gene conversion is not significantly affected when RAD50, RAD59, TID1 (RDH54), SRS2, or SGS1 is deleted. Various double mutations largely eliminate both gene conversion and BIR, including rad51Delta rad50Delta, rad51Delta rad59Delta, and rad54Delta tid1Delta. These results demonstrate that there is a RAD51- and RAD54-independent BIR pathway that requires RAD59, TID1, RAD50, and presumably MRE11 and XRS2. The similar genetic requirements for BIR and telomere maintenance in the absence of telomerase also suggest that these two processes proceed by similar mechanisms.

Factors affecting ectopic gene conversion in mice.

Duplicated genes and repetitive sequences are distributed throughout the genomes of complex organisms. The homology between related sequences can promote nonallelic (ectopic) recombination, including gene conversion and reciprocal exchange. Resolution of these events can result in translocations, deletions, or other harmful rearrangements. In yeast, ectopic recombination between sequences on nonhomologous chromosomes occurs at high frequency. Because the mammalian genome is replete with duplicated sequences and repetitive elements, high levels of ectopic exchange would cause aneuploidy and genome instability. To understand the factors regulating ectopic recombination in mice, we evaluated the effects of homology length on gene conversion between unlinked sequences in the male germline. Previously, we found high levels of gene conversion between lacZ transgenes containing 2557 bp of homology. We report here that genetic background can play a major role in ectopic recombination; frequency of gene conversion was reduced by more than an order of magnitude by transferring the transgenes from a CF1 strain background to C57BL/6J. Additionally, conversion rates decreased as the homology length decreased. Sequences sharing 1214 bp of sequence identity underwent ectopic conversion less frequently than a pair sharing 2557 bp of identity, while 624 bp was insufficient to catalyze gene conversion at significant levels. These results suggest that the germline recombination machinery in mammals has evolved in a way that prevents high levels of ectopic recombination between smaller classes of repetitive sequences, such as the Alu family. Additionally, genomic location appeared to influence the availability of sequences for ectopic recombination.

Minisatellite Origins in Yeast and Humans

Variable numbers of tandem repeats are valuable markers in genetic studies. The arrays of interest are simple microsatellites, containing repetitions of 1-5 nucleotides, and minisatellites, with multiple iterations of approximately 10 to 100 bp. Microsatellite origins can be explained by replication errors in regions fortuitously containing two or more adjacent short repeats. Microsatellite variation arises by replication errors in the absence of mismatch correction (R. Parsons et al., 1993, Cell 75: 1227-1236; M. Strand et al., 1993, Nature 365: 274-276). Variation in the size of minisatellites is thought to involve homologous recombination processes, including gene conversions (J. Buard and G. Vergnaud, 1994, EMBO J. 13: 3203-3210; A. J. Jeffreys et al., 1994, Nature Genet. 6: 136-145) and possibly unequal exchanges among repeats. However, the origins of minisatellites are less obvious. The probability of finding a direct tandem repeat of minisatellite size by chance alone is very low [<4-(10 to 100)]. Here we report the finding of short direct repeats of 5 to 10 bp flanking many yeast and human minisatellites that may be involved in their origins through replication slippage or unequal crossings over.

Further characterization of the yeast pso4-1 mutant: interaction with rad51 and rad52 mutants after photoinduced psoralen lesions.

The pso4-1 mutant was characterized as deficient in some types of recombination, including gene conversion, crossing over, and intrachromosomal recombination. The mode of interaction between pso4-1 and rad51 and between pso4-1 and rad52 mutants indicated that the PSO4 gene belongs to the RAD52 epistasis group for strand-break repair. Moreover, the presence of the pso4-1 mutation decreased 8-MOP-photoinduced mutagenesis of the rad51 and rad52 mutants. Complementation tests using heterozygous diploid strains showed that the pso4 protein might interact with the rad52 protein during repair of 8-mop photolesions. The pso4-1 mutant, even though defective in inter- and intra-chromosomal recombination, conserves the ability for plasmid integration of circular and linear plasmid DNA. On the other hand, similar to the rad51 mutant, pso4-1 was able to incise but did not restore high-molecular-weight DNA during the repair of cross links induced by 8-MOP plus UVA. These results, together with those of previous reports, indicate that the PSO4 gene belongs to the RAD52 DNA repair group and its product participates in the DNA rejoining step of the repair of cross-link lesions, which are crucial for induced mutagenesis and recombinogenesis.

Recombination and mammalian genome evolution

Several lines of evidence are presented which suggest that sequence G + C content and recombination frequency are related in mammals: (i) chromosome G + C content is positively correlated to chiasmata density; (ii) the non-pairing region of the Y chromosome has one of the lowest G + C contents of any chromosomal segment; (iii) a reduction in the rate of recombination at several loci is mirrored by a decrease in G + C content; and (iv) when compared with humans, mice have a lower variance in chiasmata density which is reflected in a lower variance in G + C content. The observed relation between recombination frequency and sequence G + C content provides an elegant explanation of why gene density is higher in G + C rich isochores than in other parts of the genome, and why long interspersed elements (LINES) are exclusive to G + C poor isochores. However, the cause of the relation is as yet unknown. Several possibilities are considered, including gene conversion.

Distinguishing the Forces Controlling Variation at the Xdh Locus in Drosophila Pseudoobscura

Fifty-eight isochromosomal lines sampled from two natural populations of Drosophila pseudoobscura in California and one from Bogota, Colombia, were examined using four-cutter restriction mapping. A 4.6-kb region of the xanthine dehydrogenase locus was probed and 66 of 135 restriction sites scored were polymorphic. This predicts that on average every 12th bp would be polymorphic in this region for the genes surveyed if polymorphism occurred randomly along the coding region. In addition, there were 12 insertion/deletion polymorphisms. Forty-nine distinct haplotypes were recognized in the 58 lines examined. The most common haplotype obtained a frequency of only 5%. Measures of base pair heterozygosity (0.0097) and linkage disequilibrium lead to a predicted population size in the range of 1.2-2.4 X 10(6) for the species. High levels of recombination (including gene conversion) can be inferred from the presence of all four gametic types in the data set.