Comparative genomics of ParaHox clusters of teleost fishes: gene cluster breakup and the retention of gene sets following whole genome duplications.

BACKGROUND: The evolutionary lineage leading to the teleost fish underwent a whole genome duplication termed FSGD or 3R in addition to two prior genome duplications that took place earlier during vertebrate evolution (termed 1R and 2R). Resulting from the FSGD, additional copies of genes are present in fish, compared to tetrapods whose lineage did not experience the 3R genome duplication. Interestingly, we find that ParaHox genes do not differ in number in extant teleost fishes despite their additional genome duplication from the genomic situation in mammals, but they are distributed over twice as many paralogous regions in fish genomes. RESULTS: We determined the DNA sequence of the entire ParaHox C1 paralogon in the East African cichlid fish Astatotilapia burtoni, and compared it to orthologous regions in other vertebrate genomes as well as to the paralogous vertebrate ParaHox D paralogons. Evolutionary relationships among genes from these four chromosomal regions were studied with several phylogenetic algorithms. We provide evidence that the genes of the ParaHox C paralogous cluster are duplicated in teleosts, just as it had been shown previously for the D paralogon genes. Overall, however, synteny and cluster integrity seems to be less conserved in ParaHox gene clusters than in Hox gene clusters. Comparative analyses of non-coding sequences uncovered conserved, possibly co-regulatory elements, which are likely to contain promoter motives of the genes belonging to the ParaHox paralogons. CONCLUSION: There seems to be strong stabilizing selection for gene order as well as gene orientation in the ParaHox C paralogon, since with a few exceptions, only the lengths of the introns and intergenic regions differ between the distantly related species examined. The high degree of evolutionary conservation of this gene cluster's architecture in particular - but possibly clusters of genes more generally - might be linked to the presence of promoter, enhancer or inhibitor motifs that serve to regulate more than just one gene. Therefore, deletions, inversions or relocations of individual genes could destroy the regulation of the clustered genes in this region. The existence of such a regulation network might explain the evolutionary conservation of gene order and orientation over the course of hundreds of millions of years of vertebrate evolution. Another possible explanation for the highly conserved gene order might be the existence of a regulator not located immediately next to its corresponding gene but further away since a relocation or inversion would possibly interrupt this interaction. Different ParaHox clusters were found to have experienced differential gene loss in teleosts. Yet the complete set of these homeobox genes was maintained, albeit distributed over almost twice the number of chromosomes. Selection due to dosage effects and/or stoichiometric disturbance might act more strongly to maintain a modal number of homeobox genes (and possibly transcription factors more generally) per genome, yet permit the accumulation of other (non regulatory) genes associated with these homeobox gene clusters.

Evolution of gene expression in fire ants: the effects of developmental stage, caste, and species.

Ants provide remarkable examples of equivalent genotypes developing into divergent and discrete phenotypes. Diploid eggs can develop either into queens, which specialize in reproduction, or workers, which participate in cooperative tasks such as building the nest, collecting food, and rearing the young. In contrast, the differentiation between males and females generally depends upon whether eggs are fertilized, with fertilized (diploid) eggs giving rise to females and unfertilized (haploid) eggs giving rise to males. To obtain a comprehensive picture of the relative contributions of gender (sex), caste, developmental stage, and species divergence to gene expression evolution, we investigated gene expression patterns in pupal and adult queens, workers, and males of two species of fire ants, Solenopsis invicta and S. richteri. Microarray hybridizations revealed that variation in gene expression profiles is influenced more by developmental stage than by caste membership, sex, or species identity. The second major contributor to variation in gene expression was the combination of sex and caste. Although workers and queens share equivalent diploid nuclear genomes, they have highly distinctive patterns of gene expression in both the pupal and the adult stages, as might be expected given their extraordinary level of phenotypic differentiation. Overall, the difference in the proportion of differentially expressed genes was greater between workers and males than between workers and queens or queens and males, consistent with the fact that workers and males share neither gender nor reproductive capability. Moreover, between-species comparisons revealed that the greatest difference in gene expression patterns occurred in adult workers, a finding consistent with the fact that adult workers most directly experience the distinct external environments characterizing the different habitats occupied by the two species. Thus, much of the evolution of gene expression in ants may occur in the worker caste, despite the fact that these individuals are largely or completely sterile. Analyses of gene expression evolution revealed a combination of positive selection and relaxation of stabilizing selection as important factors driving the evolution of such genes.

Evidence for adaptive radiation from a phylogenetic study of plant defenses.

One signature of adaptive radiation is a high level of trait change early during the diversification process and a plateau toward the end of the radiation. Although the study of the tempo of evolution has historically been the domain of paleontologists, recently developed phylogenetic tools allow for the rigorous examination of trait evolution in a tremendous diversity of organisms. Enemy-driven adaptive radiation was a key prediction of Ehrlich and Raven's coevolutionary hypothesis [Ehrlich PR, Raven PH (1964) Evolution 18:586-608], yet has remained largely untested. Here we examine patterns of trait evolution in 51 North American milkweed species (Asclepias), using maximum likelihood methods. We study 7 traits of the milkweeds, ranging from seed size and foliar physiological traits to defense traits (cardenolides, latex, and trichomes) previously shown to impact herbivores, including the monarch butterfly. We compare the fit of simple random-walk models of trait evolution to models that incorporate stabilizing selection (Ornstein-Ulenbeck process), as well as time-varying rates of trait evolution. Early bursts of trait evolution were implicated for 2 traits, while stabilizing selection was implicated for several others. We further modeled the relationship between trait change and species diversification while allowing rates of trait evolution to vary during the radiation. Species-rich lineages underwent a proportionately greater decline in latex and cardenolides relative to species-poor lineages, and the rate of trait change was most rapid early in the radiation. An interpretation of this result is that reduced investment in defensive traits accelerated diversification, and disproportionately so, early in the adaptive radiation of milkweeds.

Mutations in penicillin-binding protein (PBP) genes and in non-PBP genes during selection of penicillin-resistant Streptococcus gordonii.

Penicillin resistance in Streptococcus spp. involves multiple mutations in both penicillin-binding proteins (PBPs) and non-PBP genes. Here, we studied the development of penicillin resistance in the oral commensal Streptococcus gordonii. Cyclic exposure of bacteria to twofold-increasing penicillin concentrations selected for a progressive 250- to 500-fold MIC increase (from 0.008 to between 2 and 4 microg/ml). The major MIC increase (> or = 35-fold) was related to non-PBP mutations, whereas PBP mutations accounted only for a 4- to 8-fold additional increase. PBP mutations occurred in class B PBPs 2X and 2B, which carry a transpeptidase domain, but not in class A PBP 1A, 1B, or 2A, which carry an additional transglycosylase domain. Therefore, we tested whether inactivation of class A PBPs affected resistance development in spite of the absence of mutations. Deletion of PBP 1A or 2A profoundly slowed down resistance development but only moderately affected resistance in already highly resistant mutants (MIC = 2 to 4 microg/ml). Thus, class A PBPs might facilitate early development of resistance by stabilizing penicillin-altered peptidoglycan via transglycosylation, whereas they might be less indispensable in highly resistant mutants which have reestablished a penicillin-insensitive cell wall-building machinery. The contribution of PBP and non-PBP mutations alone could be individualized in DNA transformation. Both PBP and non-PBP mutations conferred some level of intrinsic resistance, but combining the mutations synergized them to ensure high-level resistance (> or = 2 microg/ml). The results underline the complexity of penicillin resistance development and suggest that inhibition of transglycosylase might be an as yet underestimated way to interfere with early resistance development.

Effects of selection and drift on G matrix evolution in a heterogeneous environment: a multivariate Qst-Fst Test with the freshwater snail Galba truncatula.

Unraveling the effect of selection vs. drift on the evolution of quantitative traits is commonly achieved by one of two methods. Either one contrasts population differentiation estimates for genetic markers and quantitative traits (the Q(st)-F(st) contrast) or multivariate methods are used to study the covariance between sets of traits. In particular, many studies have focused on the genetic variance-covariance matrix (the G matrix). However, both drift and selection can cause changes in G. To understand their joint effects, we recently combined the two methods into a single test (accompanying article by Martin et al.), which we apply here to a network of 16 natural populations of the freshwater snail Galba truncatula. Using this new neutrality test, extended to hierarchical population structures, we studied the multivariate equivalent of the Q(st)-F(st) contrast for several life-history traits of G. truncatula. We found strong evidence of selection acting on multivariate phenotypes. Selection was homogeneous among populations within each habitat and heterogeneous between habitats. We found that the G matrices were relatively stable within each habitat, with proportionality between the among-populations (D) and the within-populations (G) covariance matrices. The effect of habitat heterogeneity is to break this proportionality because of selection for habitat-dependent optima. Individual-based simulations mimicking our empirical system confirmed that these patterns are expected under the selective regime inferred. We show that homogenizing selection can mimic some effect of drift on the G matrix (G and D almost proportional), but that incorporating information from molecular markers (multivariate Q(st)-F(st)) allows disentangling the two effects.

Wound infection after excision and primary midline closure for pilonidal disease: risk factor analysis to improve patient selection.

BACKGROUND: Excision and primary midline closure for pilonidal disease (PD) is a simple procedure; however, it is frequently complicated by infection and prolonged healing. The aim of this study was to analyze risk factors for surgical site infection (SSI) in this context. METHODS: All consecutive patients undergoing excision and primary closure for PD from January 2002 through October 2008 were retrospectively assessed. The end points were SSI, as defined by the Center for Disease Control, and time to healing. Univariable and multivariable risk factor analyses were performed. RESULTS: One hundred thirty-one patients were included [97 men (74%), median age = 24 (range 15-66) years]. SSI occurred in 41 (31%) patients. Median time to healing was 20 days (range 12-76) in patients without SSI and 62 days (range 20-176) in patients with SSI (P < 0.0001). In univariable and multivariable analyses, smoking [OR = 2.6 (95% CI 1.02, 6.8), P = 0.046] and lack of antibiotic prophylaxis [OR = 5.6 (95% CI 2.5, 14.3), P = 0.001] were significant predictors for SSI. Adjusted for SSI, age over 25 was a significant predictor of prolonged healing. CONCLUSION: This study suggests that the rate of SSI after excision and primary closure of PD is higher in smokers and could be reduced by antibiotic prophylaxis. SSI significantly prolongs healing time, particularly in patients over 25 years.

Selection of patients with germline MLH1 mutated Lynch syndrome by determination of MLH1 methylation and BRAF mutation.

Lynch syndrome is one of the most common hereditary colorectal cancer (CRC) syndrome and is caused by germline mutations of MLH1, MSH2 and more rarely MSH6, PMS2, MLH3 genes. Whereas the absence of MSH2 protein is predictive of Lynch syndrome, it is not the case for the absence of MLH1 protein. The purpose of this study was to develop a sensitive and cost effective algorithm to select Lynch syndrome cases among patients with MLH1 immunohistochemical silencing. Eleven sporadic CRC and 16 Lynch syndrome cases with MLH1 protein abnormalities were selected. The BRAF c.1799T> A mutation (p.Val600Glu) was analyzed by direct sequencing after PCR amplification of exon 15. Methylation of MLH1 promoter was determined by Methylation-Sensitive Single-Strand Conformation Analysis. In patients with Lynch syndrome, there was no BRAF mutation and only one case showed MLH1 methylation (6%). In sporadic CRC, all cases were MLH1 methylated (100%) and 8 out of 11 cases carried the above BRAF mutation (73%) whereas only 3 cases were BRAF wild type (27%). We propose the following algorithm: (1) no further molecular analysis should be performed for CRC exhibiting MLH1 methylation and BRAF mutation, and these cases should be considered as sporadic CRC; (2) CRC with unmethylated MLH1 and negative for BRAF mutation should be considered as Lynch syndrome; and (3) only a small fraction of CRC with MLH1 promoter methylation but negative for BRAF mutation should be true Lynch syndrome patients. These potentially Lynch syndrome patients should be offered genetic counselling before searching for MLH1 gene mutations.

Pervasive positive selection on duplicated and nonduplicated vertebrate protein coding genes.

A stringent branch-site codon model was used to detect positive selection in vertebrate evolution. We show that the test is robust to the large evolutionary distances involved. Positive selection was detected in 77% of 884 genes studied. Most positive selection concerns a few sites on a single branch of the phylogenetic tree: Between 0.9% and 4.7% of sites are affected by positive selection depending on the branches. No functional category was overrepresented among genes under positive selection. Surprisingly, whole genome duplication had no effect on the prevalence of positive selection, whether the fish-specific genome duplication or the two rounds at the origin of vertebrates. Thus positive selection has not been limited to a few gene classes, or to specific evolutionary events such as duplication, but has been pervasive during vertebrate evolution.

Aging and motor programming: differential effects according to the selection of action

There are controversial reports about the effect of aging on movement preparation, and it is unclear to which extent cognitive and/or motor related cerebral processes may be affected. This study examines the age effects on electro-cortical oscillatory patterns during various motor programming tasks, in order to assess potential differences according to the mode of action selection. Twenty elderly (EP, 60-84 years) and 20 young (YP, 20-29 years) participants with normal cognition underwent 3 pre-cued response tasks (S1-S2 paradigm). S1 carried either complete information on response side (Full; stimulus-driven motor preparation), no information (None; general motor alertness), or required free response side selection (Free; internally-driven motor preparation). Electroencephalogram (EEG) was recorded using 64 surface electrodes. Alpha (8-12 Hz) desynchronization (ERD)/synchronization (ERS) and motor-related amplitude asymmetries (MRAA) were analyzed during the S1-S2 interval. Reaction times (RTs) to S2 were slower in EP than YP, and in None than in the other 2 tasks. There was an Age x Task interaction due to increased RTs in Free compared to Full in EP only. Central bilateral and midline activation (alpha ERD) was smaller in EP than YP in None. In Full just before S2, readiness to move was reflected by posterior midline inhibition (alpha ERS) in both groups. In Free, such inhibition was present only in YP. Moreover, MRAA showed motor activity lateralization in both groups in Full, but only in YP in Free. The results indicate reduced recruitment of motor regions for motor alertness in the elderly. They further show less efficient cerebral processes subtending free selection of movement in elders, suggesting reduced capacity for internally-driven action with age.

Inconsistency between different measures of sexual selection.

Measuring the intensity of sexual selection is of fundamental importance to the study of sexual dimorphism, population dynamics, and speciation. Several indices, pools of individuals, and fitness proxies are used in the literature, yet their relative performances are strongly debated. Using 12 independent common lizard populations, we manipulated the adult sex ratio, a potentially important determinant of the intensity of sexual selection at a particular time and place. We investigated differences in the intensity of sexual selection, as estimated using three standard indices of sexual selection-the standardized selection gradient (β'), the opportunity of selection (I), and the Bateman gradient (βss)--calculated for different pools of individuals and different fitness proxies. We show that results based on estimates of I were the opposite of those derived from the other indices, whereas results based on estimates of β' were consistent with predictions derived from knowledge about the species' mating system. In addition, our estimates of the strength and direction of sexual selection depended on both the fitness proxy used and the pool of individuals included in the analysis. These observations demonstrate inconsistencies in distinct measures of sexual selection and underscore the need for caution when comparing studies and species.

Mitochondrial targeting adaptation of the hominoid-specific glutamate dehydrogenase driven by positive Darwinian selection.

Many new gene copies emerged by gene duplication in hominoids, but little is known with respect to their functional evolution. Glutamate dehydrogenase (GLUD) is an enzyme central to the glutamate and energy metabolism of the cell. In addition to the single, GLUD-encoding gene present in all mammals (GLUD1), humans and apes acquired a second GLUD gene (GLUD2) through retroduplication of GLUD1, which codes for an enzyme with unique, potentially brain-adapted properties. Here we show that whereas the GLUD1 parental protein localizes to mitochondria and the cytoplasm, GLUD2 is specifically targeted to mitochondria. Using evolutionary analysis and resurrected ancestral protein variants, we demonstrate that the enhanced mitochondrial targeting specificity of GLUD2 is due to a single positively selected glutamic acid-to-lysine substitution, which was fixed in the N-terminal mitochondrial targeting sequence (MTS) of GLUD2 soon after the duplication event in the hominoid ancestor approximately 18-25 million years ago. This MTS substitution arose in parallel with two crucial adaptive amino acid changes in the enzyme and likely contributed to the functional adaptation of GLUD2 to the glutamate metabolism of the hominoid brain and other tissues. We suggest that rapid, selectively driven subcellular adaptation, as exemplified by GLUD2, represents a common route underlying the emergence of new gene functions.

Sex-dependent selection on an autosomal melanic female ornament promotes the evolution of sex ratio bias.

Sex-dependent selection often leads to spectacularly different phenotypes in males and females. In species in which sexual dimorphism is not complete, it is unclear which benefits females and males derive from displaying a trait that is typical of the other sex. In barn owls (Tyto alba), females exhibit on average larger black eumelanic spots than males but members of the two sexes display this trait in the same range of possible values. In a 12-year study, we show that selection exerted on spot size directly or on genetically correlated traits strongly favoured females with large spots and weakly favoured males with small spots. Intense directional selection on females caused an increase in spot diameter in the population over the study period. This increase is due to a change in the autosomal genes underlying the expression of eumelanic spots but not of sex-linked genes. Female-like males produced more daughters than sons, while male-like females produced more sons than daughters when mated to a small-spotted male. These sex ratio biases appear adaptive because sons of male-like females and daughters of female-like males had above-average survival. This demonstrates that selection exerted against individuals displaying a trait that is typical of the other sex promoted the evolution of specific life history strategies that enhance their fitness. This may explain why in many organisms sexual dimorphism is often not complete.

Decreasing population selection rates of resistance mutation K65R over time in HIV-1 patients receiving combination therapy including tenofovir.

OBJECTIVES: The use of tenofovir is highly associated with the emergence of mutation K65R, which confers broad resistance to nucleoside/nucleotide analogue reverse transcriptase inhibitors (NRTIs), especially when tenofovir is combined with other NRTIs also selecting for K65R. Although recent HIV-1 treatment guidelines discouraging these combinations resulted in reduced K65R selection with tenofovir, updated information on the impact of currently recommended regimens on the population selection rate of K65R is presently lacking. METHODS: In this study, we evaluated changes over time in the selection rate of resistance mutation K65R in a large population of 2736 HIV-1-infected patients failing combination antiretroviral treatment between 2002 and 2010. RESULTS: The K65R resistance mutation was detected in 144 patients, a prevalence of 5.3%. A large majority of observed K65R cases were explained by the use of tenofovir, reflecting its wide use in clinical practice. However, changing patterns over time in NRTIs accompanying tenofovir resulted in a persistent decreasing probability of K65R selection by tenofovir-based therapy. The currently recommended NRTI combination tenofovir/emtricitabine was associated with a low probability of K65R emergence. For any given dual NRTI combination including tenofovir, higher selection rates of K65R were consistently observed with a non-nucleoside reverse transcriptase inhibitor than with a protease inhibitor as the third agent. DISCUSSION: Our finding of a stable time trend of K65R despite elevated use of tenofovir illustrates increased potency of current HIV-1 therapy including tenofovir.