The role of gene conversion in determining sequence variation and divergence in the Est-5 gene family in Drosophila pseudoobscura.

Nucleotide sequences of eight Est-5A and Est-5C genes corresponding to previously sequenced Est-5B genes in Drosophila pseudoobscura were determined to compare patterns of polymorphism and divergence among members of this small gene family. The three esterase genes were also sequenced from D. persimilis and D. miranda for interspecific comparisons. The data provide evidence that gene conversion between loci contributes to polymorphism and to the homogenization of the Est5 genes. For Est-5B, which encodes one of the most highly polymorphic proteins in Drosophila, 12% of the segregating amino acid variants appear to have been introduced via gene conversion from other members of the gene family. Interlocus gene conversion can also explain high sequence similarity, especially at synonymous sites, between Est-5B and Est-5A. Tests of neutrality using interspecific comparisons show that levels of polymorphism conform to neutral expectations at each Est-5 locus. However, McDonald-Kreitman tests based on intraspecific gene comparisons indicate that positive selection on amino acids has accompanied Est-5 gene duplication and divergence in D. pseudoobscura.     

The sorghum photoperiod sensitivity gene, Ma3, encodes a phytochrome B.

The Ma3 gene is one of six genes that regulate the photoperiodic sensitivity of flowering in sorghum (Sorghum bicolor [L.] Moench). The ma3R mutation of this gene causes a phenotype that is similar to plants that are known to lack phytochrome B, and ma3 sorghum lacks a 123-KD phytochrome that predominates in light-grown plants and that is present in non-ma3 plants. A population segregating for Ma3 and ma3 was created and used to identify two randomly amplified polymorphic DNA markers linked to Ma3. These two markers were cloned and mapped in a recombinant inbred population as restriction fragment length polymorphisms. cDNA clones of PHYA and PHYC were cloned and sequenced from a cDNA library prepared from green sorghum leaves. Using a genome-walking technique, a 7941-bp partial sequence of PHYB, was determined from genomic DNA from ma3 sorghum. PHYA, PHYB, and PHYC all mapped to the same linkage group. The Ma3-linked markers mapped with PHYB more than 121 centimorgans from PHYA and PHYC. A frameshift mutation resulting in a premature stop codon was found in the PHYB sequence from ma3 sorghum. Therefore, we conclude that the Ma3 locus in sorghum is a PHYB gene that encodes a 123-kD phytochrome.     

Phytochrome gene expression and phylogenetic analysis in the short-day plant Pharbitis nil (Convolvulaceae): Differential regulation by light and an endogenous clock

To investigate the role of distinct phytochrome pools in photoperiodic timekeeping, we characterized four phytochrome genes in the short-day plant Pharbitis nil. Each PHY gene had different photosensory properties and sensitivity to night break that inhibits flowering. During extended dark periods, PHYE, PHYB, and PHYC mRNA accumulation exhibited a circadian rhythmicity indicative of control by an endogenous clock. Phylogenetic analysis recovered four clades of angiosperm phytochrome genes, phyA, phyB, phyC, and phyE. All except the phyE clade included sequences from both monocots and eudicots. In addition, phyA is sister to phyC and phyE sister to phyB, with gymnosperm sequences sister to either the phyA-phyC clade or to the phyB-phyE clade. These results suggest that a single duplication occurred in an ancestral seed plant before the divergence of extant gymnosperms from angiosperms and that two subsequent duplications occurred in an ancestral angiosperm before the divergence of monocots from eudicots. Thus in P. nil, a multigene family with different patterns of mRNA abundance in light and darkness contributes to the total phytochrome pool: one pool is light labile (phyA), whereas the other is light stable (phyB and phyE). In addition, PHYC mRNA represents a third phytochrome pool with intermediate photosensory properties.     

Molecular evolution of the toll-like receptor multigene family in birds

Toll-like receptors (TLR) are membrane-bound sensors of the innate immune system that recognize invariant and distinctive molecular features of invading microbes and are also essential for initiating adaptive immunity in vertebrates. The genetic variation at TLR genes has been directly related to differential pathogen outcomes in humans and livestock. Nonetheless, new insights about the impact of TLRs polymorphism on the evolutionary ecology of infectious diseases can be gained through the investigation of additional vertebrate groups not yet investigated in detail. In this study, we have conducted the first characterization of the entire TLR multigene family (N = 10 genes) in non-model avian species. Using primers targeting conserved coding regions, we aimed at amplifying large segments of the extracellular domains (275-435 aa) involved in pathogen recognition across seven phylogenetically diverse bird species. Our analyses suggest avian TLRs are dominated by stabilizing selection, suggesting that slow rates of nonsynonymous substitution help preserve biological function. Overall, mean values of ω (= d(n)/d(s)) at each TLR locus ranged from 0.196 to 0.517. However, we also found patterns of positive selection acting on specific amino acid sites that could be linked to species-specific differences in pathogen-associated molecular pattern recognition. Only 39 of 2,875 (～1.35%) of the codons analyzed exhibited significant patterns of positive selection. At least one half of these positively selected codons can be mapped to putative ligand-binding regions, as suggested by crystallographic structures of TLRs and their ligands and mutagenic analyses. We also surveyed TLR polymorphism in wild populations of two bird species, the Lesser Kestrel Falco naumanni and the House Finch Carpodacus mexicanus. In general, avian TLRs displayed low to moderate single nucleotide polymorphism levels and an excess of silent nucleotide substitutions, but also conspicuous instances of positive directional selection. In particular, TLR5 and TLR15 exhibited the highest degree of genetic polymorphism and the highest occurrence of nonconservative amino acid substitutions. This study provides critical primers and a first look at the evolutionary patterns and implications of TLR polymorphism in non-model avian species and extends the list of candidate loci for avian eco-immunogenetics beyond the widely employed genes of the Major Histocompatibility Complex (MHC).     

The relation between the neutrality index for mitochondrial genes and the distribution of mutational effects on fitness

We explore factors affecting patterns of polymorphism and divergence (as captured by the neutrality index) at mammalian mitochondrial loci. To do this, we develop a population genetic model that incorporates a fraction of neutral amino acid sites, mutational bias, and a probability distribution of selection coefficients against new nonsynonymous mutations. We confirm, by reanalyzing publicly available datasets, that the mitochondrial cyt-b gene shows a broad range of neutrality indices across mammalian taxa, and explore the biological factors that can explain this observation. We find that observed patterns of differences in the neutrality index, polymorphism, and divergence are not caused by differences in mutational bias. They can, however, be explained by a combination of a small fraction of neutral amino acid sites, weak selection acting on most amino acid mutations, and differences in effective population size among taxa.     

Reproductive protein evolution within and between species: maintenance of divergent ZP3 alleles in Peromyscus

In a variety of animal taxa, proteins involved in reproduction evolve more rapidly than nonreproductive proteins. Most studies of reproductive protein evolution, however, focus on divergence between species, and little is known about differentiation among populations within a species. Here we investigate the molecular population genetics of the protein ZP3 within two Peromyscus species. ZP3 is an egg coat protein involved in primary binding of egg and sperm and is essential for fertilization. We find that amino acid polymorphism in the sperm-combining region of ZP3 is high relative to silent polymorphism in both species of Peromyscus . In addition, while there is geographical structure at a mitochondrial gene ( Cytb ), a nuclear gene ( Lcat ) and eight microsatellite loci, we find no evidence for geographical structure at Zp3 in Peromyscus truei . These patterns are consistent with the maintenance of ZP3 alleles by balancing selection, possibly due to sexual conflict or pathogen resistance. However, we do not find evidence that reinforcement promotes ZP3 diversification; allelic variation in P. truei is similar among populations, including populations allopatric and sympatric with sibling species. In fact, most alleles are present in all populations sampled across P. truei's range. While additional data are needed to identify the precise evolutionary forces responsible for sequence variation in ZP3, our results suggest that in Peromyscus , selection to maintain divergent alleles within species contributes to the pattern of rapid amino acid substitution observed among species.     

Contrasting modes of natural selection acting on pigmentation genes in the Drosophila dunni subgroup

Genes that encode for divergent adaptive traits may have genealogies that contrast with those from loci that are not functionally involved in differentiation. Here, we examine DNA sequence variation among the species of the eastern Caribbean Drosophila dunni subgroup at two loci, yellow and dopa decaboxylase (Ddc), which both play integral roles in pigmentation patterning of adult Drosophila. Phylogenetic analyses of these loci produce gene genealogies with topologies that mirror those described for other nuclear genes: the six morphologically distinct species within the subgroup are divided into only three lineages, with one lineage containing four species that share extensive ancestral polymorphism. At the Ddc locus these major lineages are delineated only by silent site variation. We observe a significantly higher rate of synonymous site divergence than non-synonymous divergence, consistent with strong purifying selection acting on the locus. In contrast, the yellow locus exhibits patterns of amino acid divergence and nucleotide diversity that are consistent with recent diversifying selection acting in two different lineages. This selection appears to be targeting amino acid variants in the signal sequence of the Yellow protein, a region which is tightly constrained among members of the larger D. cardini radiation. This result highlights not only the potential importance of yellow in the evolution of divergent pigmentation patterns among members of the D. dunni subgroup, but also hints that variation in signal peptide sequences may play a role in phenotypic diversification.     

Signatures of selection and sex-specific expression variation of a novel duplicate during the evolution of the Drosophila desaturase gene family

The tempo and mode of evolution of loci with a large effect on adaptation and reproductive isolation will influence the rate of evolutionary divergence and speciation. Desaturase loci are involved in key biochemical changes in long-chain fatty acids. In insects, these have been shown to influence adaptation to starvation or desiccation resistance and in some cases act as important pheromones. The desaturase gene family of Drosophila is known to have evolved by gene duplication and diversification, and at least one locus shows rapid evolution of sex-specific expression variation. Here, we examine the evolution of the gene family in species representing the Drosophila phylogeny. We find that the family includes more loci than have been previously described. Most are represented as single-copy loci, but we also find additional examples of duplications in loci which influence pheromone blends. Most loci show patterns of variation associated with purifying selection, but there are strong signatures of diversifying selection in new duplicates. In the case of a new duplicate of desat1 in the obscura group species, we show that strong selection on the coding sequence is associated with the evolution of sex-specific expression variation. It seems likely that both sexual selection and ecological adaptation have influenced the evolution of this gene family in Drosophila.     

Adaptive evolution in the photosensory domain of phytochrome A in early angiosperms

Flowering plant diversity now far exceeds the combined diversity of all other plant groups. Recently identified extant remnants of the earliest-diverging lines suggest that the first angiosperms may have lived in shady, disturbed, and moist understory habitats, and that the aquatic habit also arose early. This would have required the capacity to begin life in dimly lit environments. If so, evolution in light-sensing mechanisms may have been crucial to their success. The photoreceptor phytochrome A is unique among angiosperm phytochromes in its capacity to serve a transient role under conditions where an extremely high sensitivity is required. We present evidence of altered functional constraints between phytochrome A (PHYA) and its paralog, PHYC. Tests for selection suggest that an elevation in nonsynonymous rates resulted from an episode of selection along the branch leading to all angiosperm PHYA sequences. Most nucleotide sites (95%) are selectively constrained, and the ratio of nonsynonymous to synonymous substitutions on branches within the PHYA clade does not differ from the ratio on the branches in the PHYC clade. Thus, positive selection at a handful of sites, rather than relaxation of selective constraints, apparently has played a major role in the evolution of the photosensory domain of phytochrome A. The episode of selection occurred very early in the history of flowering plants, suggesting that innovation in phyA may have given the first angiosperms some adaptive advantage.     

Phylogenetic relationships of B-related phytochromes in the Brassicaceae: Redundancy and the persistence of phytochrome D

Plants use phytochrome (phy) photoreceptors to detect and respond to changes in the quantities and proportions of red (R) and far-red (FR) light in their environments. The principal mediators of responses to R and FR in Arabidopsis thaliana are phyA and phyB, which are found in all angiosperms surveyed. The present study is concerned with a phytochrome gene pair in Arabidopsis, PHYB and PHYD, which are of relatively recent origin, share high sequence identity, and are partially redundant. Our data suggest that the duplication occurred after the mustard family (Brassicaceae) diverged from its closest relatives but before the radiation of extant Brassicaceae, and that both copies have persisted for up to 40myr. We detected no evidence of positive selection in the divergence of PHYD from PHYB; the evolution of both sequences is constrained by purifying selection. Levels of diversity at both loci are among the lowest observed at nuclear genes in A. thaliana. In common with other loci in A. thaliana, PHYB and PHYD showed elevated levels of intraspecific replacement variation, and each showed an excess of rare nucleotide polymorphisms, consistent with a recent, rapid population expansion. Our results are consistent with the functional importance of amino acid divergence in the central regions of phyB and phyD and suggest specific sites for mutagenesis that may yield insights into the functional differences of phyB and phyD.     

Clines and adaptive evolution in the methuselah gene region in Drosophila melanogaster

In an effort to characterize further the patterns of selection and adaptive evolution at the methuselah locus in Drosophila species, we extended an analysis of geographical variation to include single nucleotide polymorphisms (SNPs) in adjacent genes on either side of the mth locus, and examined the molecular variation in a neighbouring methuselah paralogue (mth2). An analysis of 13 SNPs spanning a region of nearly 19 kilobases surrounding the mth locus demonstrated that a clinal pattern associated with the most common mth haplotype does not extend to adjacent gene loci, providing compelling evidence that the clinal pattern results from selection on as yet unidentified sites associated with the functional mth locus. mth2 exhibited a significant pattern of adaptive divergence among D. melanogaster, D. simulans and D. yakuba similar to that seen at mth. However, Ka : Ks ratios indicate a difference in levels of functional constraint at the two methuselah, loci with mth2 exhibiting a five- to six-fold reduction in levels of amino acid divergence relative to mth.     

Comparative population genetics of the panicoid grasses: sequence polymorphism, linkage disequilibrium and selection in a diverse sample of sorghum bicolor

Levels of genetic variation and linkage disequilibrium (LD) are critical factors in association mapping methods as well as in identification of loci that have been targets of selection. Maize, an outcrosser, has a high level of sequence variation and a limited extent of LD. Sorghum, a closely related but largely self-pollinating panicoid grass, is expected to have higher levels of LD. As a first step in estimation of population genetic parameters in sorghum, we surveyed 27 diverse S. bicolor accessions for sequence variation at a total of 29,186 bp in 95 short regions derived from genetically mapped RFLPs located throughout the genome. Consistent with its higher level of inbreeding, the extent of LD is at least severalfold greater in sorghum than in maize. Total sequence variation in sorghum is about fourfold lower than that in maize, while synonymous variation is fivefold lower, suggesting a smaller effective population size in sorghum. Because we surveyed a species-wide sample, the mating system, which primarily affects population-level diversity, may not be primarily responsible for this difference. Comparisons of polymorphism and divergence suggest that both directional and diversifying selection have played important roles in shaping variation in the sorghum genome.     

Evidence that the phytochrome gene family in black cottonwood has one PHYA locus and two PHYB loci but lacks members of the PHYC/F and PHYE subfamilies

The phytochrome photoreceptors play important roles in the photoperiodic control of vegetative bud set, growth cessation, dormancy induction, and cold-hardiness in trees. Interestingly, ecotypic differences in photoperiodic responses are observed in many temperate- zone tree species. Northern and southern ecotypes of black cottonwood (Populus trichocarpa Torr. & Gray), for example, exhibit marked differences in the timing of short-day-induced bud set and growth cessation, and these responses are controlled by phytochrome. Therefore, as a first step toward determining the molecular genetic basis of photoperiodic ecotypes in trees, we characterized the phytochrome gene (PHY) family in black cottonwood. We recovered fragments of one PHYA and two PHYB using PCR-based cloning and by screening a genomic library. Results from Southern analyses confirmed that black cottonwood has one PHYA locus and two PHYB loci, which we arbitrarily designated PHYB1 and PHYB2. Phylogenetic analyses which included PHY from black cottonwood, Arabidopsis thaliana and tomato (Solanum lycopersicum) suggest that the PHYB/D duplications in these species occurred independently. When Southern blots were probed with PHYC, PHYE, and PHYE heterologous probes, the strongest bands that we detected were those of black cottonwood PHYA and/or PHYB. These results suggest that black cottonwood lacks members of the PHYC/F and PHYE subfamilies. Although black cottonwood could contain additional PHY that are distantly related to known angiosperm PHY, our results imply that the PHY family of black cottonwood is less complex than that of other well-characterized dicot species such as Arabidopsis and tomato. Based on Southern analyses of five black cottonwood genotypes representing three photoperiodic ecotypes, substantial polymorphism was detected for at least one of the PHYB loci but not for the PHYA locus. The novel character of the PHY family in black cottonwood, as well as the differences in polymorphism we observed between the PHYA and PHYB subfamilies, indicates that a number of fundamental macro- and microevolutionary questions remain to be answered about the PHY family in dicots.      

Sequence divergence at the putative flowering time locus COL1 in Brassicaceae

An insertion/deletion polymorphism (Ind2) in the Brassica nigra CONSTANS LIKE 1 (Bni COL1) gene was previously found to be associated with variation in flowering time. In the present study we examine the inter-specific divergence of COL1 in the family Brassicaceae. Analysis of codon substitution models did not reveal evidence of positive Darwinian selection, but comparisons of the COL1 gene in different species revealed a surprising number of indels. A total of 24 indels were found in the 650 bp of the middle variable region of the gene. This high number of indels could reflect a lack of constraint on length of this region of the protein, or the effect of positive selection. The number of indels was close to that expected in non-coding DNA, but the indels were longer in COL1 than those observed in non-coding regions. Reconstruction of indel evolution indicated that most indels resulted from deletions rather than insertions. The Ind2 indel that has shown association with flowering time in Brassica nigra exhibited a remarkable distribution in the Brassicaceae family, indicating that the polymorphism may have persisted more than ten million years. Considering presumed historic populations sizes of Brassicaceae species, such a long persistence time seems unlikely for a neutral polymorphism.     

Molecular signatures of divergence and selection in closely related pine taxa

Efforts to detect loci under selection in plants have mostly focussed on single species. However, assuming that intraspecific divergence may lead to speciation, comparisons of genetic variation within and among recently diverged taxa can help to locate such genes. In this study, coalescent and outlier detection methods were used to assess nucleotide polymorphism and divergence at 79 nuclear gene fragments (1212 SNPs) in 16 populations (153 individuals) of the closely related, but phenotypically and ecologically distinct, pine taxa Pinus mugo, P. uliginosa and P. uncinata across their European distributions. Simultaneously, mitochondrial DNA markers, which are maternally inherited in pines and distributed by seeds at short geographic distance, were used to assess genetic relationships of the focal populations and taxa. The majority of nuclear loci showed homogenous patterns of variation between the taxa due to a high number of shared SNPs and haplotypes, similar levels of polymorphism, and low net divergence. However, against this common genetic background and an overall low population structure within taxa at mitochondrial markers, we identified several genes showing signatures of selection, accompanied by significant intra- and interspecific divergence. Our results indicate that loci involved in species divergence may be involved in intraspecific local adaptation.     

Evidence for directional selection acting on pheromone-binding proteins in the genus Choristoneura

Patterns of nucleotide variation consistent with the action of natural selection have been discovered at a number of different gene loci. Here, pheromone-binding proteins (PBPs) are examined to determine if selection has acted to fix amino acid changes in PBPs in lineages in which pheromone changes have occurred. PBPs from five different species of moths in the genus Choristoneura were sequenced, along with the PBP of Argyrotaenia velutinana, which serves as an outgroup. Three independent major pheromone changes are represented within this group of five Choristoneura species. Two different lineages show evidence for selection based on polymorphism and divergence comparisons and comparisons of rates of replacement evolution to silent and noncoding evolution. Along one of these lineages, leading to Choristoneura fumiferana, there has been a change to an aldehyde pheromone from an acetate pheromone. The second branch does not appear to be associated with a major pheromone change. Other branches in the tree show a trend toward greater replacement fixation than expected under neutrality. This trend could reflect undetected selective events within this group of PBPs. Selection appears to have acted to fix amino acid changes in the PBP of moths from the genus Choristoneura, but it is not clear that this selection is due to pheromone changes between species.     

Bayesian analysis suggests that most amino acid replacements in Drosophila are driven by positive selection

One of the principal goals of population genetics is to understand the processes by which genetic variation within species (polymorphism) becomes converted into genetic differences between species (divergence). In this transformation, selective neutrality, near neutrality, and positive selection may each play a role, differing from one gene to the next. Synonymous nucleotide sites are often used as a uniform standard of comparison across genes on the grounds that synonymous sites are subject to relatively weak selective constraints and so may, to a first approximation, be regarded as neutral. Synonymous sites are also interdigitated with nonsynonymous sites and so are affected equally by genomic context and demographic factors. Hence a comparison of levels of polymorphism and divergence between synonymous sites and amino acid replacement sites in a gene is potentially informative about the magnitude of selective forces associated with amino acid replacements. We have analyzed 56 genes in which polymorphism data from D. simulans are compared with divergence from a reference strain of D. melanogaster. The framework of the analysis is Bayesian and assumes that the distribution of selective effects (Malthusian fitnesses) is Gaussian with a mean that differs for each gene. In such a model, the average scaled selection intensity (gamma = N(e)s) of amino acid replacements eligible to become polymorphic or fixed is -7.31, and the standard deviation of selective effects within each locus is 6.79 (assuming homoscedasticity across loci). For newly arising mutations of this type that occur in autosomal or X-linked genes, the average proportion of beneficial mutations is 19.7%. Among the amino acid polymorphisms in the sample, the expected average proportion of beneficial mutations is 47.7%, and among amino acid replacements that become fixed the average proportion of beneficial mutations is 94.3%. The average scaled selection intensity of fixed mutations is +5.1. The presence of positive selection is pervasive with the single exception of kl-5, a Y-linked fertility gene. We find no evidence that a significant fraction of fixed amino acid replacements is neutral or nearly neutral or that positive selection drives amino acid replacements at only a subset of the loci. These results are model dependent and we discuss possible modifications of the model that might allow more neutral and nearly neutral amino acid replacements to be fixed.