Selection and Biased Gene Conversion in a Multigene Family: Consequences of Interallelic Bias and Threshold Selection

In a previous paper, I investigated the interactions in a gene family of additive selection and biased gene conversion in a finite population when conversion events are rare. Here I extend my "weak-conversion limit" model by allowing biased interallelic conversion (conversion between alleles at the same locus) of arbitrary frequency and various threshold selection schemes for rare interlocus conversion events. I suggest that it is not unreasonable for gene families to experience threshold fitness functions, and show that certain types of thresholds can greatly constrain the rate at which advantageous alleles are fixed as compared to other fitness schemes, such as additive selection. It is also shown that the double sampling process operating on a gene family in a finite population (sampling over the number of genes in the gene family and over the number of individuals in the population) can have interesting consequences. For selectively neutral alleles that experience interallelic bias, the probability of fixation at each single locus may be essentially neutral, but the cumulative effects on the entire gene family of small departures from neutrality can be significant, especially if the gene family is large. Thus, in some situations, gene families can respond to directional forces that are weak in comparison to drift at single loci.     

Biased Gene Conversion, Copy Number, and Apparent Mutation Rate Differences within Chloroplast and Bacterial Genomes

We investigate the possibility that differences between synonymous substitution rates of organelle and bacterial genes differing only in copy number may be due to conversion bias. We find that the rather large observed difference in the synonymous rates between genes in the single copy and inverted-repeat regions of chloroplasts can be accounted for by a very small bias against new mutants. More generally, differences in the within-organelle fixation probability result in different apparent mutation rates as measured by the expected rate of appearance of cells homoplasmic for new mutants. Thus, differences in intracellular population parameters rather than molecular mechanisms can account for some variation in the apparent mutation rates of organelle genes, and possibly in other systems with variable numbers of gene copies. On the other hand, our analysis suggests that conversion bias is not a likely explanation for relatively low mutation rates observed near the replication origin of bacterial chromosomes.     

Duplication, Selection and Gene Conversion in a Drosophila mojavensis Female Reproductive Protein Family

Protein components of the Drosophila male ejaculate, several of which evolve rapidly, are critical modulators of reproductive success. Recent studies of female reproductive tract proteins indicate they also are extremely divergent between species, suggesting that reproductive molecules may coevolve between the sexes. Our current understanding of intersexual coevolution, however, is severely limited by the paucity of genetic and evolutionary studies on the female molecules involved. Physiological evidence of ejaculate–female coadaptation, paired with a promiscuous mating system, makes Drosophila mojavensis an exciting model system in which to study the evolution of reproductive proteins. Here we explore the evolutionary dynamics of a five-paralog gene family of female reproductive proteases within populations of D. mojavensis and throughout the repleta species group. We show that the proteins have experienced ongoing gene duplication and adaptive evolution and further exhibit dynamic patterns of pseudogenation, copy number variation, gene conversion, and selection within geographically isolated populations of D. mojavensis. The integration of these patterns in a single gene family has never before been documented in a reproductive protein.     

Gene Conversion, Linkage, and the Evolution of Multigene Families

The evolution of the probabilities of genetic identity within and between the loci of a multigene family is investigated. Unbiased gene conversion, equal crossing over, random genetic drift, and mutation to new alleles are incorporated. Generations are discrete and nonoverlapping; the diploid, monoecious population mates at random. The linkage map is arbitrary, and the location dependence of the probabilities of identity is formulated exactly. The greatest of the rates of gene conversion, random drift, and mutation is epsilon much less than 1. For interchromosomal conversion, the equilibrium probabilities of identity are within order epsilon [i.e., O(epsilon)] of those in a simple model that has no location dependence and, at equilibrium, no linkage disequilibrium. At equilibrium, the linkage disequilibria are of O(epsilon); they are evaluated explicitly with an error of O(epsilon 2); they may be negative if symmetric heteroduplexes occur. The ultimate rate and pattern of convergence to equilibrium are within O(epsilon 2) and O(epsilon), respectively, of that of the same simple model. If linkage is loose (i.e., all the crossover rates greatly exceed epsilon, though they may still be much less than 1/2), the linkage disequilibria are reduced to O(epsilon) in a time of O(-ln epsilon). If intrachromosomal conversion is incorporated, the same results hold for loose linkage, except that, if the crossover rates are much less than 1/2, then the linkage disequilibria generally exceed those for pure interchromosomal conversion.     

Sequence Identity in an Early Chorion Multigene Family Is the Result of Localized Gene Conversion

The multigene families that encode the chorion (eggshell) of the silk moth, Bombyx mori, are closely linked on one chromosome. We report here the isolation and characterization of two segments, totaling 102 kb of genomic DNA, containing the genes expressed during the early period of choriogenesis. Most of these early genes can be divided into two multigene families, ErA and ErB, organized into five divergently transcribed ErA/ErB gene pairs. Nucleotide sequence identity in the major coding regions of the ErA genes was 96%, while nucleotide sequence identity for the ErB major coding regions was only 63%. Selection pressure on the encoded proteins cannot explain this difference in the level of sequence conservation between the ErA and ErB gene families, since when only fourfold redundant codon positions are considered, the divergence within the ErA genes is 8%, while the divergence within the ErB genes (corrected for multiple substitutions at the same site) is 110%. The high sequence identity of the ErA major exons can be explained by sequence exchange events similar to gene conversion localized to the major exon of the ErA genes. These gene conversions are correlated with the presence of clustered copies of the nucleotide sequence GGXGGX, encoding paired glycine residues. This sequence has previously been correlated with gradients of gene conversion that extend throughout the coding and noncoding regions of the High-cysteine (Hc) chorion genes of B. mori. We suggest that the difference in the extent of the conversion tracts in these gene families reflects a tendency for these recombination events to become localized over time to the protein encoding regions of the major exons.     

The Evolution of Multigene Families under Intrachromosomal Gene Conversion

A model for the evolution of the probabilities of genetic identity within and between loci of a multigene family in a finite population is formulated and investigated. Unbiased intrachromosomal gene conversion, equal crossing over between tandemly repeated genes, random genetic drift and mutation to new alleles are incorporated. Generations are discrete and nonoverlapping; the diploid, monoecious population mates at random. Formulas for the equilibrium values of the probabilities of identity and a cubic equation for the rate of convergence are deduced. Numerical examples indicate the following. The amount of homology at equilibrium generally decreases as the mutation rate, the population size and the number of repeats increase; it may increase or decrease with increasing crossover rate. The intralocus homology has an intermediate minimum, whereas the interlocus homology increases, as the rate of gene conversion increases. The intralocus homology decreases, whereas the interlocus homology increases, as the proportion of symmetric heteroduplexes increases. The characteristic convergence time can be sufficiently short to imply that intrachromosomal gene conversion may be an important mechanism for maintaining sequence homogeneity among repeated genes. The convergence time decreases as the conversion rate and the proportion of symmetric heteroduplexes increase; although exceptions occur, it generally increases as the population size and the number of repeats increase; it may increase or decrease with increasing crossover rate.     

Segregating Sites in a Gene Conversion Model with Mutation

Formulae for the expectation and variance of the number of segregating and homogeneous sites in a sample of two chromosomes are found. The model includes gene conversion and infinitely-many-alleles mutation in a coalescent framework. The corresponding infinitely-many-sites model limits are also found. The formulae for the expectation are extended to any sample size. Comparisons are drawn between the pure mutation model and the model where gene conversion has been added.     

The Roles of Gene Duplication,Gene Conversion and Positive Selection in Rodent Esp and Mup Pheromone Gene Families with Comparison to the Abp Family

Three proteinaceous pheromone families, the androgen-binding proteins (ABPs), the exocrine-gland secreting peptides (ESPs) and the major urinary proteins (MUPs) are encoded by large gene families in the genomes of Mus musculus and Rattus norvegicus. We studied the evolutionary histories of the Mup and Esp genes and compared them with what is known about the Abp genes. Apparently gene conversion has played little if any role in the expansion of the mouse Class A and Class B Mup genes and pseudogenes, and the rat Mups. By contrast, we found evidence of extensive gene conversion in many Esp genes although not in all of them. Our studies of selection identified at least two amino acid sites in β-sheets as having evolved under positive selection in the mouse Class A and Class B MUPs and in rat MUPs. We show that selection may have acted on the ESPs by determining K_a/K_s for Exon 3 sequences with and without the converted sequence segment. While it appears that purifying selection acted on the ESP signal peptides, the secreted portions of the ESPs probably have undergone much more rapid evolution. When the inner gene converted fragment sequences were removed, eleven Esp paralogs were present in two or more pairs with K_a/K_s >1.0 and thus we propose that positive selection is detectable by this means in at least some mouse Esp paralogs. We compare and contrast the evolutionary histories of all three mouse pheromone gene families in light of their proposed functions in mouse communication.     

CBFA2T1,a Gene Rearranged in Human Leukemia, Is a Member of a Multigene Family

MTG8(HGMW-approved symbolCBFA2T1) was originally identified as one of the loci involved in the t(8;21)(q22;q22) of acute myeloid leukemia. We characterize two humanMTG8-related genes,MTGR1andMTGR2(HGMW-approved symbolsCBFA2T2andCBFA2T3). The former is duplicated in mouse, one locus possibly being a retroposon. MultipleMTG8-related sequences are found in several vertebrate species, from fish to mammals, albeit not in a urodele.MTGR2maps to 16q24 and, likeMTG8andMTGR1,is close to one of three loci encoding a syntrophin (dystrophin-associated proteins). Moreover, an alternativeMTGR1promoter/5′ exon is contained within the α1-syntrophin locus. Thus, the two classes of genes may define novel paralogous groups.MTGR1is expressed mainly in brain, whileMTGR2is expressed in the thymus and possibly in monocytes. LikeMTG8, MTGR1is transcribed into a number of isoforms due to alternative splicing of different 5′ exons onto a common splice acceptor site. Comparison of the three predicted human MTG8-related polypeptides to theirDrosophilacounterpart (nervy) highlights four separate regions of sequence conservation that may correspond to distinct domains. The most NH₂-terminal of these is proportionately more conserved among the human polypeptides, presumably due to specific structural/functional constraints.     

GC-Content Evolution in Bacterial Genomes: The Biased Gene Conversion Hypothesis Expands

The characterization of functional elements in genomes relies on the identification of the footprints of natural selection. In this quest, taking into account neutral evolutionary processes such as mutation and genetic drift is crucial because these forces can generate patterns that may obscure or mimic signatures of selection. In mammals, and probably in many eukaryotes, another such confounding factor called GC-Biased Gene Conversion (gBGC) has been documented. This mechanism generates patterns identical to what is expected under selection for higher GC-content, specifically in highly recombining genomic regions. Recent results have suggested that a mysterious selective force favouring higher GC-content exists in Bacteria but the possibility that it could be gBGC has been excluded. Here, we show that gBGC is probably at work in most if not all bacterial species. First we find a consistent positive relationship between the GC-content of a gene and evidence of intra-genic recombination throughout a broad spectrum of bacterial clades. Second, we show that the evolutionary force responsible for this pattern is acting independently from selection on codon usage, and could potentially interfere with selection in favor of optimal AU-ending codons. A comparison with data from human populations shows that the intensity of gBGC in Bacteria is comparable to what has been reported in mammals. We propose that gBGC is not restricted to sexual Eukaryotes but also widespread among Bacteria and could therefore be an ancestral feature of cellular organisms. We argue that if gBGC occurs in bacteria, it can account for previously unexplained observations, such as the apparent non-equilibrium of base substitution patterns and the heterogeneity of gene composition within bacterial genomes. Because gBGC produces patterns similar to positive selection, it is essential to take this process into account when studying the evolutionary forces at work in bacterial genomes.     

Positive Darwinian Selection Promotes Heterogeneity Among Members of the Antifreeze Protein Multigene Family

A variety of organisms have independently evolved proteins exhibiting antifreeze activity that allows survival at subfreezing temperatures. The antifreeze proteins (AFPs) bind ice nuclei and depress the freezing point by a noncolligative absorption–inhibition mechanism. Many organisms have a heterogeneous suite of AFPs with variation in primary sequence between paralogous loci. Here, we demonstrate that the diversification of the AFP paralogues is promoted by positive Darwinian selection in two independently evolved AFPs from fish and beetle. First, we demonstrate an elevated rate of nonsynonymous substitutions compared to synonymous substitutions in the mature protein coding region. Second, we perform phylogeny-based tests of selection to demonstrate a subset of codons is subjected to positive selection. When mapped onto the three-dimensional structure of the fish antifreeze type III antifreeze structure, these codons correspond to amino acid positions that surround but do not interrupt the putative ice-binding surface. The selective agent may be related to efficient binding to diverse ice surfaces or some other aspect of AFP function. Received: 27 February 2001 / Accepted: 12 September 2001     

Biased Gene Conversion Skews Allele Frequencies in Human Populations,Increasing the Disease Burden of Recessive Alleles

Gene conversion results in the nonreciprocal transfer of genetic information between two recombining sequences, and there is evidence that this process is biased toward G and C alleles. However, the strength of GC-biased gene conversion (gBGC) in human populations and its effects on hereditary disease have yet to be assessed on a genomic scale. Using high-coverage whole-genome sequences of African hunter-gatherers, agricultural populations, and primate outgroups, we quantified the effects of GC-biased gene conversion on population genomic data sets. We find that genetic distances (F_ST and population branch statistics) are modified by gBGC. In addition, the site frequency spectrum is left-shifted when ancestral alleles are favored by gBGC and right-shifted when derived alleles are favored by gBGC. Allele frequency shifts due to gBGC mimic the effects of natural selection. As expected, these effects are strongest in high-recombination regions of the human genome. By comparing the relative rates of fixation of unbiased and biased sites, the strength of gene conversion was estimated to be on the order of Nb ≈ 0.05 to 0.09. We also find that derived alleles favored by gBGC are much more likely to be homozygous than derived alleles at unbiased SNPs (+42.2% to 62.8%). This results in a curse of the converted, whereby gBGC causes substantial increases in hereditary disease risks. Taken together, our findings reveal that GC-biased gene conversion has important population genetic and public health implications.     

Analysis of POFUT1 Gene Mutation in a Chinese Family with Dowling-Degos Disease

Dowling-Degos disease (DDD) is an autosomal dominant genodermatosis characterized by reticular pigmented anomaly mainly affecting flexures. Though KRT5 has been identified to be the causal gene of DDD, the heterogeneity of this disease was displayed: for example, POFUT1 and POGLUT1 were recently identified and confirmed to be additional pathogenic genes of DDD. To identify other DDD causative genes, we performed genome-wide linkage and exome sequencing analyses in a multiplex Chinese DDD family, in which the KRT5 mutation was absent. Only a novel 1-bp deletion (c.246+5delG) in POFUT1 was found. No other novel mutation or this deletion was detected in POFUT1 in a second DDD family and a sporadic DDD case by Sanger Sequencing. The result shows the genetic-heterogeneity and complexity of DDD and will contribute to the further understanding of DDD genotype/phenotype correlations and to the pathogenesis of this disease.     

Genetic Differentiation of Transposable Elements under Mutation and Unbiased Gene Conversion

A model is developed to predict the extent of genetic differentiation in a family of transposable elements under the combined effects of genetic drift, transposition, mutation and unbiased gene conversion. The model is based on simplifying assumptions that are valid when transposition is always to new sites and copy number per site is low. In the absence of gene conversion, the degree of differentiation as measured by the probability of identity of different elements is the same as at a single locus with the same mutation rate but in a population of effective size Nc/2, where N is the population size and c is the number of copies per individual. The inclusion of unbiased gene conversion does not significantly change this result. If, as seems to be the case, families of transposable elements are relatively homogeneous, then the model implies either that mutation rates for transposable elements are much lower than at comparable single-copy loci or that some other force, such as natural selection or biased gene conversion, is at work. Transposition is a very ineffective force for homogenizing a family of transposable elements.     

Phylogenetic and Structural Relationships of the PR5 Gene Family Reveal an Ancient Multigene Family Conserved in Plants and Select Animal Taxa

Pathogenesis-related group 5 (PR5) plant proteins include thaumatin, osmotin, and related proteins, many of which have antimicrobial activity. The recent discovery of PR5-like (PR5-L) sequences in nematodes and insects raises questions about their evolutionary relationships. Using complete plant genome data and discovery of multiple insect PR5-L sequences, phylogenetic comparisons among plants and animals were performed. All PR5/PR5-L protein sequences were mined from genome data of a member of each of two main angiosperm groups—the eudicots (Arabidoposis thaliana) and the monocots (Oryza sativa)—and from the Caenorhabditis nematode (C. elegans and C. briggsase). Insect PR5-L sequences were mined from EST databases and GenBank submissions from four insect orders: Coleoptera (Diaprepes abbreviatus and Biphyllus lunatus), Orthoptera (Schistocerca gregaria), Hymenoptera (Lysiphlebus testaceipes), and Hemiptera (Toxoptera citricida). Parsimony and Bayesian phylogenetic analyses showed that the PR5 family is paraphyletic in plants, likely arising from 10 genes in a common ancestor to monocots and eudicots. After evolutionary divergence of monocots and eudicots, PR5 genes increased asymmetrically among the 10 clades. Insects and nematodes contain multiple sequences (seven PR5-Ls in nematodes and at least three in some insects) all related to the same plant clade, with nematode and insect sequences separating as two clades. Protein structural homology modeling showed strong similarity among animal and plant PR5/PR5-Ls, with divergence only in surface-exposed loops. Sequence and structural conservation among PR5/PR5-Ls suggests an important and conserved role throughout the evolutionary divergence of the diverse organisms from which they reside. [Reviewing Editor: Dr. Rafael Zardoya]     

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.     

Targeted Heritable Mutation and Gene Conversion by Cas9-CRISPR in Caenorhabditis elegans

We have achieved targeted heritable genome modification in Caenorhabditis elegans by injecting mRNA of the nuclease Cas9 and Cas9 guide RNAs. This system rapidly creates precise genomic changes, including knockouts and transgene-instructed gene conversion.