Strong Selective Constraints Research Articles

Testing the Neutral Fixation of Hetero-Oligomerism in the Archaeal Chaperonin CCT

The evolutionary transition from homo-oligomerism to hetero-oligomerism in multimeric proteins and its contribution to function innovation and organism complexity remain to be investigated. Here, we undertake the challenge of contributing to this theoretical ground by investigating the hetero-oligomerism in the molecular chaperonin cytosolic chaperonin containing tailless complex polypeptide 1 (CCT) from archaea. CCT is amenable to this study because, in contrast to eukaryotic CCTs where sub-functionalization after gene duplication has been taken to completion, archaeal CCTs present no evidence for subunit functional specialization. Our analyses yield additional information to previous reports on archaeal CCT paralogy by identifying new duplication events. Analyses of selective constraints show that amino acid sites from 1 subunit have fixed slightly deleterious mutations at inter-subunit interfaces after gene duplication. These mutations have been followed by compensatory mutations in nearby regions of the same subunit and in the interface contact regions of its paralogous subunit. The strong selective constraints in these regions after speciation support the evolutionary entrapment of CCTs as hetero-oligomers. In addition, our results unveil different evolutionary dynamics depending on the degree of CCT hetero-oligomerism. Archaeal CCT protein complexes comprising 3 distinct classes of subunits present 2 evolutionary processes. First, slightly deleterious and compensatory mutations were fixed neutrally at inter-subunit regions. Second, sub-functionalization may have occurred at substrate-binding and adenosine triphosphate-binding regions after the 2nd gene duplication event took place. CCTs with 2 distinct types of subunits did not present evidence of sub-functionalization. Our results provide the 1st in silico evidence for the neutral fixation of hetero-oligomerism in archaeal CCTs and provide information on the evolution of hetero-oligomerism toward sub-functionalization in archaeal CCTs.

Complete plastid genome sequences suggest strong selection for retention of photosynthetic genes in the parasitic plant genus Cuscuta

BackgroundPlastid genome content and protein sequence are highly conserved across land plants and their closest algal relatives. Parasitic plants, which obtain some or all of their nutrition through an attachment to a host plant, are often a striking exception. Heterotrophy can lead to relaxed constraint on some plastid genes or even total gene loss. We sequenced plastid genomes of two species in the parasitic genus Cuscuta along with a non-parasitic relative, Ipomoea purpurea, to investigate changes in the plastid genome that may result from transition to the parasitic lifestyle.ResultsAside from loss of all ndh genes, Cuscuta exaltata retains photosynthetic and photorespiratory genes that evolve under strong selective constraint. Cuscuta obtusiflora has incurred substantially more change to its plastid genome, including loss of all genes for the plastid-encoded RNA polymerase. Despite extensive change in gene content and greatly increased rate of overall nucleotide substitution, C. obtusiflora also retains all photosynthetic and photorespiratory genes with only one minor exception.ConclusionAlthough Epifagus virginiana, the only other parasitic plant with its plastid genome sequenced to date, has lost a largely overlapping set of transfer-RNA and ribosomal genes as Cuscuta, it has lost all genes related to photosynthesis and maintains a set of genes which are among the most divergent in Cuscuta. Analyses demonstrate photosynthetic genes are under the highest constraint of any genes within the plastid genomes of Cuscuta, indicating a function involving RuBisCo and electron transport through photosystems is still the primary reason for retention of the plastid genome in these species.

Highly constrained proteins contain an unexpectedly large number of amino acid tandem repeats

Single-amino-acid tandem repeats are very common in mammalian proteins but their function and evolution are still poorly understood. Here we investigate how the variability and prevalence of amino acid repeats are related to the evolutionary constraints operating on the proteins. We find a significant positive correlation between repeat size difference and protein nonsynonymous substitution rate in human and mouse orthologous genes. This association is observed for all the common amino acid repeat types and indicates that rapid diversification of repeat structures, involving both trinucleotide slippage and nucleotide substitutions, preferentially occurs in proteins subject to low selective constraints. However, strikingly, we also observe a significant negative correlation between the number of repeats in a protein and the gene nonsynonymous substitution rate, particularly for glutamine, glycine, and alanine repeats. This implies that proteins subject to strong selective constraints tend to contain an unexpectedly high number of repeats, which tend to be well conserved between the two species. This is consistent with a role for selection in the maintenance of a significant number of repeats. Analysis of the codon structure of the sequences encoding the repeats shows that codon purity is associated with high repeat size interspecific variability. Interestingly, polyalanine and polyglutamine repeats associated with disease show very distinctive features regarding the degree of repeat conservation and the protein sequence selective constraints.

Potential Source of Francisella tularensis Live Vaccine Strain Attenuation Determined by Genome Comparison

Francisella tularensis is a bacterial pathogen that causes the zoonotic disease tularemia and is important to biodefense. Currently, the only vaccine known to confer protection against tularemia is a specific live vaccine strain (designated LVS) derived from a virulent isolate of Francisella tularensis subsp. holarctica. The origin and source of attenuation of this strain are not known. To assist with the design of a defined live vaccine strain, we sought to determine the genetic basis of the attenuation of LVS. This analysis relied primarily on the comparison between the genome of LVS and Francisella tularensis holarctica strain FSC200, which differ by only 0.08% of their nucleotide sequences. Under the assumption that the attenuation was due to a loss of function(s), only coding regions were examined in this comparison. To complement this analysis, the coding regions of two slightly more distantly related Francisella tularensis strains were also compared against the LVS coding regions. Thirty-five genes show unique sequence variations predicted to alter the protein sequence in LVS compared to the other Francisella tularensis strains. Due to these polymorphisms, the functions of 15 of these genes are very likely lost or impaired. Seven of these genes were demonstrated to be under stronger selective constraints, suggesting that they are the most probable to be the source of LVS attenuation and useful for a newly defined vaccine.

Towards a crucial role for DC-SIGN in tuberculosis and beyond

The C-type lectin DC-SIGN has recently received considerable attention in the context of tuberculosis and other infectious diseases. Clinical investigations, together with genetic association studies, strongly support the notion that variation in the level of expression of DC-SIGN - but not changes in the structure of the protein - might have a strong impact on the susceptibility to and pathogenesis of several infectious diseases, including tuberculosis. In addition, efforts to decipher the evolutionary history of the gene encoding DC-SIGN clearly demonstrated that this gene is under strong selective constraints that have prevented the accumulation of amino acid changes over time. Altogether, these findings suggest that DC-SIGN might play a crucial part in host immunity to pathogens and possibly beyond, at an early stage of human development.

Functional noncoding sequences derived from SINEs in the mammalian genome

Recent comparative analyses of mammalian sequences have revealed that a large number of nonprotein-coding genomic regions are under strong selective constraint. Here, we report that some of these loci have been derived from a newly defined family of ancient SINEs (short interspersed repetitive elements). This is a surprising result, as SINEs and other transposable elements are commonly thought to be genomic parasites. We named the ancient SINE family AmnSINE1, for Amniota SINE1, because we found it to be present in mammals as well as in birds, and some copies predate the mammalian-bird split 310 million years ago (Mya). AmnSINE1 has a chimeric structure of a 5S rRNA and a tRNA-derived SINE, and is related to five tRNA-derived SINE families that we characterized here in the coelacanth, dogfish shark, hagfish, and amphioxus genomes. All of the newly described SINE families have a common central domain that is also shared by zebrafish SINE3, and we collectively name them the DeuSINE (Deuterostomia SINE) superfamily. Notably, of the approximately 1000 still identifiable copies of AmnSINE1 in the human genome, 105 correspond to loci phylogenetically highly conserved among mammalian orthologs. The conservation is strongest over the central domain. Thus, AmnSINE1 appears to be the best example of a transposable element of which a significant fraction of the copies have acquired genomic functionality.

Structure and age of genetic diversity of dengue virus type 2 in Thailand

Dengue virus type 2 (DENV-2) is a common viral infection and an important health concern in South-East Asia. To determine the molecular evolution of DENV-2 in Thailand, 105 isolates of the E (envelope) gene and 10 complete genomes sampled over a 27 year period were sequenced. Phylogenetic analysis of these data revealed that three genotypes of DENV-2 have circulated in Thailand, although, since 1991, only viruses assigned to Asian genotype I have been sampled from the population. A broader analysis of 35 complete genomes of DENV-2 revealed that most amino acids are subject to strong selective constraints, indicative of widespread purifying selection against deleterious mutations. This was further supported by an analysis of genome-wide substitution rates, which indicated that DENV-2 fixes approximately 10 mutations per genome per year, far lower than expected from its mutational dynamics. Finally, estimates of the age of DENV-2 were remarkably consistent among genes, indicating that the current genetic diversity in this virus probably arose within the last 120 years, concordant with the first determination of the aetiology of dengue disease.

Understanding the Degradation of Hominid Gene Control

Peter D. Keightley, Martin J. Lercher, Adam Eyre-Walker Recently, two groups have examined the level of sequence constraint in noncoding DNA flanking mammalian genes, and appear to have found conflicting results. By comparing 500-bp blocks in mice and rats, we found that mean nucleotide divergence within 2 kb of the start and stop codons of protein-coding genes is substantially lower than that of introns, and decreases when approaching the coding sequence [1]. If nucleotide changes within introns are largely free from selection, this implies that noncoding blocks close to genes evolved under selective constraints, presumably because they contain gene expression control regions. In contrast, we find that upstream sequences in hominids do not evolve slower than introns, while downstream regions are under about half of the constraint seen in murids [1]. By analysing a similar set of noncoding DNA sequences, Bush and Lahn also found that the mean level of selective constraints in upstream regions between humans and chimpanzees is very low. However, their slightly more complex main analysis was to search for 16-bp sequences within upstream regions that are strongly conserved between humans, mice, and either dogs or chickens. They then examined the divergence between humans and chimpanzees at the flanking nucleotides, finding substantially reduced divergence compared with the genomic mean. This demonstrated selective constraints at certain upstream sequences in hominids. An analogous analysis of mouse–rat sequences showed that the selective constraints are about twice as strong in murids as in hominids [2]. These two findings—on one hand, a near absence of selective constraints in blocks upstream of hominid genes [1], and on the other, evidence for strong selective constraints in these regions [2]—appear to contradict each other. How can we square the two sets of results? The answer is rather simple—windows with high conservation scores are relatively rare, and they contribute little to the mean calculated over 500-bp windows (unfortunately, Bush and Lahn do not tell us the fraction of 59 alignments within high conservation scores). Bush and Lahn also suggest that the apparent discrepancy ‘‘likely results from the fact that in large 500-bp blocks, functional elements that are under constraint are mixed with large sections of nonfunctional DNA, which are not under constraint’’ [2]. We believe that this interpretation, while formally correct, obscures important and interesting information that can be gained from combining the two studies. Some sequences outside the conserved 16-mers identified by Bush and Lahn are also likely to be functional, since the same 500-bp regions (largely ‘‘nonfunctional’’ according to Bush and Lahn) show strong evidence of evolutionary constraints between mice and rats [1]. Bush and Lahn also note that constraint, on either side of conserved windows, is greater in murids than in hominids. However, they observe a much smaller difference than that seen in our analysis. This is deceptive because by concentrating attention on regions that are conserved between humans, mice, and dogs, they ignore the fact that there might be many more highly conserved regions in murids than there are in hominids. In summary, there is no conflict between our results and those of Bush and Lahn; they concentrate their attention on a preselected subset of the sites we considered and so have a different perspective on the problem. What is clear from both studies is that there is a qualitative difference in the level of conservation in the 59 flanking sequences between murids and hominids. We have argued that this is likely to be due to the fixation of slightly deleterious mutations in hominids that are otherwise selectively eliminated in rodents. Differences in constraints between hominids and murids demonstrate that the overwhelming majority of changes at upstream regulatory sites have only small effects on fitness. This has counterintuitive consequences: to obtain a comprehensive list of human regulatory sites, it might be better to examine conservation in murid rather than hominid genomes. “

The Heritage of Pathogen Pressures and Ancient Demography in the Human Innate-Immunity CD209/CD209L Region

The innate immunity system constitutes the first line of host defense against pathogens. Two closely related innate immunity genes, CD209 and CD209L, are particularly interesting because they directly recognize a plethora of pathogens, including bacteria, viruses, and parasites. Both genes, which result from an ancient duplication, possess a neck region, made up of seven repeats of 23 amino acids each, known to play a major role in the pathogen-binding properties of these proteins. To explore the extent to which pathogens have exerted selective pressures on these innate immunity genes, we resequenced them in a group of samples from sub-Saharan Africa, Europe, and East Asia. Moreover, variation in the number of repeats of the neck region was defined in the entire Human Genome Diversity Panel for both genes. Our results, which are based on diversity levels, neutrality tests, population genetic distances, and neck-region length variation, provide genetic evidence that CD209 has been under a strong selective constraint that prevents accumulation of any amino acid changes, whereas CD209L variability has most likely been shaped by the action of balancing selection in non-African populations. In addition, our data point to the neck region as the functional target of such selective pressures: CD209 presents a constant size in the neck region populationwide, whereas CD209L presents an excess of length variation, particularly in non-African populations. An additional interesting observation came from the coalescent-based CD209 gene tree, whose binary topology and time depth (approximately 2.8 million years ago) are compatible with an ancestral population structure in Africa. Altogether, our study has revealed that even a short segment of the human genome can uncover an extraordinarily complex evolutionary history, including different pathogen pressures on host genes as well as traces of admixture among archaic hominid populations.

Aflatoxin biosynthesis gene clusters and flanking regions

To compare the biosynthetic gene cluster sequences of the main aflatoxin (AF)-producing Aspergillus species. Sequencing was on fosmid clones selected by homology to Aspergillus parasiticus sequence. Alignments revealed that gene order is conserved among AF gene clusters of Aspergillus nomius, A. parasiticus, two sclerotial morphotypes of Aspergillus flavus, and an unnamed Aspergillus sp. Phylogenetic relationships were established using the maximum likelihood method implemented in PAUP. Based on the Eurotiomycete/Sordariomycete divergence time, the A. flavus-type cluster has been maintained for at least 25 million years. Such conservation of the genes and gene order reflects strong selective constraints on rearrangement. Phylogenetic comparison of individual genes in the cluster indicated that ver-1, which has homology to a melanin biosynthesis gene, experienced selective forces distinct from the other pathway genes. Sequences upstream of the polyketide synthase-encoding gene vary among the species, but a four-gene sugar utilization cluster at the distal end is conserved, indicating a functional relationship between the two adjacent clusters. The high conservation of cluster components needed for AF production suggests there is an adaptive value for AFs in character-shaping niches important to those taxa. This is the first comparison of the complete nucleotide sequences of gene clusters harbouring the AF biosynthesis genes of the main AF-producing species. Such a comparison will aid in understanding how AF biosynthesis is regulated in experimental and natural environments.

The α-Globin Gene Family of an Australian Marsupial, Macropus eugenii: The Long Evolutionary History of the θ-Globin Gene and Its Functional Status in Mammals

Comparative evolutionary analyses of gene families among divergent lineages can provide information on the order and timing of major gene duplication events and evolution of gene function. Here we investigate the evolutionary history of the alpha-globin gene family in mammals by isolating and characterizing alpha-like globin genes from an Australian marsupial, the tammar wallaby, Macropus eugenii. Sequence and phylogenetic analyses indicate that the tammar alpha-globin family consists of at least four genes including a single adult-expressed gene (alpha), two embryonic/neonatally expressed genes (zeta and zeta'), and theta-globin, each orthologous to the respective alpha-, zeta-, and theta-globin genes of eutherian mammals. The results suggest that the theta-globin lineage arose by duplication of an ancestral adult alpha-globin gene and had already evolved an unusual promoter region, atypical of all known alpha-globin gene promoters, prior to the divergence of the marsupial and eutherian lineages. Evolutionary analyses, using a maximum likelihood approach, indicate that theta-globin, has evolved under strong selective constraints in both marsupials and the lineage leading to human theta-globin, suggesting a long-term functional status. Overall, our results indicate that at least a four-gene cluster consisting of three alpha-like and one beta-like globin genes linked in the order 5'-zeta-alpha-theta-omega-3' existed in the common ancestor of marsupials and eutherians. However, results are inconclusive as to whether the two tammar zeta-globin genes arose by duplication prior to the radiation of the marsupial and eutherian lineages, with maintenance of exon sequences by gene conversion, or more recently within marsupials.

Intron Size and Exon Evolution in Drosophila

We have found a negative correlation between evolutionary rate at the protein level (as measured by d(N)) and intron size in Drosophila. Although such a relation is expected if introns reduce Hill-Robertson interference within genes, it seems more likely to be explained by the higher abundance of cis-regulatory elements in introns (especially first introns) in genes under strong selective constraints.

Mammalian housekeeping genes evolve more slowly than tissue-specific genes.

Do housekeeping genes, which are turned on most of the time in almost every tissue, evolve more slowly than genes that are turned on only at specific developmental times or tissues? Recent large-scale gene expression studies enable us to have a better definition of housekeeping genes and to address the above question in detail. In this study, we examined 1581 human-mouse orthologous gene pairs for their patterns of sequence evolution, contrasting housekeeping genes with tissue-specific genes. Our results show that, in comparison to tissue-specific genes, housekeeping genes on average evolve more slowly and are under stronger selective constraints as reflected by significantly smaller values of Ka/Ks. Besides stronger purifying selection, we explored several other factors that can possibly slow down nonsynonymous rates in housekeeping genes. Although mutational bias might slightly slow the nonsynonymous rates in housekeeping genes, it is unlikely to be the major cause of the rate difference between the two types of genes. The codon usage pattern of housekeeping genes does not seem to differ from that of tissue-specific genes. Moreover, contrary to the old textbook concept, we found that approximately 74% of the housekeeping genes in our study belong to multigene families, not significantly different from that of the tissue-specific genes ( approximately 70%). Therefore, the stronger selective constraints on housekeeping genes are not due to a lower degree of genetic redundancy.

Reconstructing the diversification of subtilisins in the pathogenic fungus Metarhizium anisopliae.

Fungi secrete subtilisin proteinases to acquire nutrients and breach host barriers. Here we sought a global characterization of the diversity of subtilisins in the insect pathogen Metarhizium anisopliae. Expressed sequence tag (EST) analyses showed that a broad host range strain of M. anisopliae sf. anisopliae (strain 2575) expressed 11 subtilisins during growth on insect cuticle, the largest number of subtilisins reported from any fungus. Polymerase chain reaction amplified 10 of their orthologs from a second strain with multiple hosts (strain 820) and seven from the locust specialist M. anisopliae sf. acridum (strain 324). Analyses based on sequence similarities and exon-intron structure grouped M. anisopliae subtilisins into four clusters-a class I ("bacterial") subtilisin (Pr1C), and three clusters of proteinase K-like class II subtilisins: extracellular subfamily 1 (Pr1A, Pr1B, Pr1G, Pr1I and Pr1K), extracellular subfamily 2 (Pr1D, Pr1E, Pr1F and Pr1J) and an endocellular subtilisin (Pr1H). Phylogenetic analysis of homologous sequences from other genera revealed that this subdivision of proteinase K-like subtilisins into three subfamilies preceded speciation of major fungal lineages. However, diversification has continued during the evolution of Metarhizium subtilisins with evidence of gene duplication events after divergence of M. anisopliae sf. anisopliae and M. anisopliae sf. acridum. Comparing alignments and nonsynonymous/synonymous rates for Pr1 isoenzymes within a lineage and between lineages showed that while overall divergence of subtilisins followed neutral expectations, amino acids involved in catalysis were under strong selective constraint. This suggests that each Pr1 paralog contributes to the pathogens fitness. Furthermore, homology modeling predicted differences between the Pr1's in their secondary substrate specificities, adsorption properties to cuticle and alkaline stability, indicative of functional differences.

A molecular epidemiological study of Australian bat lyssavirus.

The genetic diversity of Australian bat lyssavirus (ABL) was investigated by comparing 24 ABL isolate glycoprotein (G) gene nucleotide sequences with those of 37 lyssaviruses representing Lyssavirus genotypes 1-6. Phylogenetic analyses indicated that ABL forms a monophyletic group separate from other lyssaviruses. This group differentiates into two clades: one associated with Pteropus (flying fox) species, the other with the insectivorous bat Saccolaimus flaviventris. Calculation of percentage nucleotide identities between isolates of the two clades revealed up to 18.7 % nucleotide sequence divergence between the two ABL variants. These observations suggest that ABL is a separate lyssavirus species with a similar epidemiology to chiropteran rabies virus (RV), where two distinct ABL variants co-exist in Australia in bat species with dissimilar ecology. Analyses of selection pressures in ABL G gene sequences provided some evidence of weak positive selection within the endodomain at amino acids 499 and 501, although in general the dominant evolutionary process observed was purifying selection. This intimates that, in nature, isolates of ABL, like those of RV, are subject to relatively strong selective constraints, suggesting a stability of host species, cell tropisms and ecological conditions.

Evolution of the sulfide-binding function within the globin multigenic family of the deep-sea hydrothermal vent tubeworm Riftia pachyptila.

The giant extracellular hexagonal bilayer hemoglobin (HBL-Hb) of the deep-sea hydrothermal vent tube worm Riftia pachyptila is able to transport simultaneously O(2) and H(2)S in the blood from the gills to a specific organ: the trophosome that harbors sulfide-oxidizing endosymbionts. This vascular HBL-Hb is made of 144 globins from which four globin types (A1, A2, B1, and B2) coevolve. The H(2)S is bound at a specific location (not on the heme site) onto two of these globin types. In order to understand how such a function emerged and evolved in vestimentiferans and other related annelids, six partial cDNAs corresponding to the six globins known to compose the multigenic family of R. pachyptila have been identified and sequenced. These partial sequences (ca. 120 amino acids, i.e., 80% of the entire protein) were used to reconstruct molecular phylogenies in order to trace duplication events that have led to the family organization of these globins and to locate the position of the free cysteine residues known to bind H(2)S. From these sequences, only two free cysteine residues have been found to occur, at positions Cys + 1 (i.e., 1 a.a. from the well-conserved distal histidine) and Cys + 11 (i.e., 11 a.a. from the same histidine) in globins B2 and A2, respectively. These two positions are well conserved in annelids, vestimentiferans, and pogonophorans, which live in sulfidic environments. The structural comparison of the hydrophobic environment that surrounds these cysteine residues (the sulfide-binding domain) using hydrophobic cluster analysis plots, together with the cysteine positions in paralogous strains, suggests that the sulfide-binding function might have emerged before the annelid radiation in order to detoxify this toxic compound. Moreover, globin evolutionary rates are highly different between paralogous strains. This suggests that either the two globin subfamilies involved in the sulfide-binding function (A2 and B2) have evolved under strong directional selective constraints (negative selection) and that the two other globins (A1 and B1) have accumulated more substitutions through positive selection or have evolved neutrally after a relaxation of selection pressures. A likely scenario on the evolution of this multigenic family is proposed and discussed from this data set.

Coding sequence divergence between two closely related plant species: Arabidopsis thaliana and Brassica rapa ssp. pekinensis.

To characterize the coding-sequence divergence of closely related genomes, we compared DNA sequence divergence between sequences from a Brassica rapa ssp. pekinensis EST library isolated from flower buds and genomic sequences from Arabidopsis thaliana. The specific objectives were (i) to determine the distribution of and relationship between K(a) and K(s), (ii) to identify genes with the lowest and highest K(a): K(s) values, and (iii) to evaluate how codon usage has diverged between two closely related species. We found that the distribution of K(a): K(s) was unimodal, and that substitution rates were more variable at nonsynonymous than synonymous sites, and detected no evidence that K(a) and K(s) were positively correlated. Several genes had K(a): K(s) values equal to or near zero, as expected for genes that have evolved under strong selective constraint. In contrast, there were no genes with K(a): K(s) >1 and thus we found no strong evidence that any of the 218 sequences we analyzed have evolved in response to positive selection. We detected a stronger codon bias but a lower frequency of GC at synonymous sites in A. thaliana than B. rapa. Moreover, there has been a shift in the profile of most commonly used synonymous codons since these two species diverged from one another. This shift in codon usage may have been caused by stronger selection acting on codon usage or by a shift in the direction of mutational bias in the B. rapa phylogenetic lineage.

Evidence for strong selective constraint acting on the nucleotide composition of 16S ribosomal RNA genes.

Previous studies have shown that the guanine plus cytosine (G+C) content of ribosomal RNAs (rRNAs) is highly correlated with bacterial growth temperatures. This correlation is strongest in the double-stranded stem regions of the rRNA, a fact that can be explained by selection for increased structural stability at high growth temperatures. In this study, we examined the single-stranded regions of 16S rRNAs. We reasoned that, since these regions of the molecule are subject to less structural constraint than the stem regions, their nucleotide content might simply reflect the overall nucleotide content of the genome. Contrary to this expectation, however, we found that all of the single-stranded regions are characterized by very high adenine (A) and relatively low cytosine (C) contents. Moreover, the nucleotide content of these single-stranded regions is surprisingly constant between species, despite dramatic differences in optimal growth temperatures, and despite large differences in the overall genomic G+C content. This provides compelling evidence for strong stabilizing selection acting on 16S rRNA single-stranded regions. We found that selection favors purines (A+G), and especially adenine (A), in the single-stranded regions of these rRNAs.

Phenotypic plasticity is the major determinant of changes in phenotypic integration in Arabidopsis.

• The way in which novel genetic variation affects the patterns of phenotypic integration in natural populations is addressed here. • An experimental study is presented of the variability in integration caused by interpopulation hybridization and consequent genetic reshuffling, as well as by changes in the physical environment in the model system Arabidopsis thaliana (Brassicaceae). • Our results show a basic invariance of sets of covarying traits in A. thaliana, with changes in nutrient availability as the principal factor accounting for major departures from the general pattern and where differences in the genetic background of the recombinant lines are less important. In A. thaliana, the relationships among vegetative and reproductive traits form distinct clusters in multivariate space. A high degree of congruence was found between differences in the multivariate mean phenotype and the pattern of phenotypic integration, as expected on the basis of recent theoretical models. • This relationship might indicate strong selective constraints acting on the specialized life history of these populations, which are spring ephemerals inhabiting ruderal habitats and prone to competition avoidance.

The molecular basis of an avian plumage polymorphism in the wild: a melanocortin-1-receptor point mutation is perfectly associated with the melanic plumage morph of the bananaquit, Coereba flaveola.

Evolution depends on natural selection acting on phenotypic variation, but the genes responsible for phenotypic variation in natural populations of vertebrates are rarely known. The molecular genetic basis for plumage color variation has not been described in any wild bird. Bananaquits (Coereba flaveola) are small passerine birds that occur as two main plumage variants, a widespread yellow morph with dark back and yellow breast and a virtually all black melanic morph. A candidate gene for this color difference is the melanocortin-1 receptor (MC1R), a key regulator of melanin synthesis in feather melanocytes. We sequenced the MC1R gene from four Caribbean populations of the bananaquit; two populations of the yellow morph and two populations containing both the yellow morph and the melanic morph. A point mutation resulting in the replacement of glutamate with lysine was present in at least one allele of the MC1R gene in all melanic birds and was absent in all yellow morph birds. This substitution probably causes the color variation, as the same substitution is responsible for melanism in domestic chickens and mice. The evolutionary relationships among the MC1R haplotypes show that the melanic alleles on Grenada and St. Vincent had a single origin. The low prevalence of nonsynonymous substitutions among yellow haplotypes suggests that they have been under stabilizing selection, whereas strong selective constraint on melanic haplotypes is absent. We conclude that a mutation in the MC1R is responsible for the plumage polymorphism in a wild bird population and that the melanic MC1R alleles in Grenada and St. Vincent bananaquit populations have a single evolutionary origin from a yellow allele.