Phylogenetic alpha and beta diversities of butterfly communities correlate with climate in the western Swiss Alps

The observation of non-random phylogenetic distribution of traits in communities provides evidence for niche-based community assembly. Environment may influence the phylogenetic structure of communities because traits determining how species respond to prevailing conditions can be phylogenetically conserved. In this study, we investigate the variation of butterfly species richness and of phylogenetic - and -diversities along temperature and plant species richness gradients. Our study indicates that butterfly richness is independently positively correlated to temperature and plant species richness in the study area. However, the variation of phylogenetic - and -diversities is only correlated to temperature. The significant phylogenetic clustering at high elevation suggests that cold temperature filters butterfly lineages, leading to communities mostly composed of closely related species adapted to those climatic conditions. These results suggest that in colder and more severe conditions at high elevations deterministic processes and not purely stochastic events drive the assemblage of butterfly communities.

Birth and expression evolution of mammalian microRNA genes.

MicroRNAs (miRNAs) are major post-transcriptional regulators of gene expression, yet their origins and functional evolution in mammals remain little understood due to the lack of appropriate comparative data. Using RNA sequencing, we have generated extensive and comparable miRNA data for five organs in six species that represent all main mammalian lineages and birds (the evolutionary outgroup) with the aim to unravel the evolution of mammalian miRNAs. Our analyses reveal an overall expansion of miRNA repertoires in mammals, with threefold accelerated birth rates of miRNA families in placentals and marsupials, facilitated by the de novo emergence of miRNAs in host gene introns. Generally, our analyses suggest a high rate of miRNA family turnover in mammals with many newly emerged miRNA families being lost soon after their formation. Selectively preserved mammalian miRNA families gradually evolved higher expression levels, as well as altered mature sequences and target gene repertoires, and were apparently mainly recruited to exert regulatory functions in nervous tissues. However, miRNAs that originated on the X chromosome evolved high expression levels and potentially diverse functions during spermatogenesis, including meiosis, through selectively driven duplication-divergence processes. Overall, our study thus provides detailed insights into the birth and evolution of mammalian miRNA genes and the associated selective forces.

Interaction between primary and secondary sporocysts of Schistosoma mansoni and the internal defence system of Biomphalaria resistant and susceptible to the parasite

The outcome of the interaction between Biomphalaria and Schistosoma mansoni depends on the response of the host internal defence system (IDS) and the escape mechanisms of the parasite. The aim of this study was to evaluate the responsiveness of the IDS (haemocytes and soluble haemolymph factors) of resistant and susceptible Biomphalaria tenagophila lineages and Biomphalaria glabrata lineages in the presence of in vitro-transformed primary sporocysts and secondary sporocysts obtained from infected B. glabrata. To do this, we assayed the cellular adhesion index (CAI), analysed viability/mortality, used fluorescent markers to evaluate the tegumental damage and transplanted secondary sporocysts. B. tenagophila Taim was more effective against primary and secondary sporocystes than the susceptible lineage and B. glabrata. Compared with secondary sporocysts exposed to B. tenagophila, primary sporocysts showed a higher CAI, a greater percentage of dead sporocysts and were labelled by lectin from Glycine max and Alexa-Fluor 488 fluorescent probes at a higher rate than the secondary sporocysts. However, the two B. tenagophila lineages showed no cercarial shedding after inoculation with secondary sporocysts. Our hypothesis that secondary sporocysts can escape the B. tenagophila IDS cannot be confirmed by the transplantation experiments. These data suggest that there are additional mechanisms involved in the lower susceptibilty of B. tenagophila to S. mansoni infection.

Comparison of spoligotyping, mycobacterial interspersed repetitive units typing and IS6110-RFLP in a study of genotypic diversity of Mycobacterium tuberculosis in Delhi, North India

The aim of the present study was to compare polymerase chain reaction (PCR)-based methods - spoligotyping and mycobacterial interspersed repetitive units (MIRU) typing - with the gold-standard IS6110 restriction fragment length polymorphism (RFLP) analysis in 101 isolates of Mycobacterium tuberculosis to determine the genetic diversity of M. tuberculosis clinical isolates from Delhi, North India. Spoligotyping resulted in 49 patterns (14 clusters); the largest cluster was composed of Spoligotype International Types (SITs)26 [Central-Asian (CAS)1-Delhi lineage], followed by SIT11 [East-African-Indian (EAI) 3-Indian lineage]. A large number of isolates (75%) belonged to genotypic lineages, such as CAS, EAI and Manu, with a high specificity for the Indian subcontinent, emphasising the complex diversity of the phylogenetically coherent M. tuberculosis in North India. MIRU typing, using 11 discriminatory loci, was able to distinguish between all but two strains based on individual patterns. IS6110-RFLP analysis (n = 80 strains) resulted in 67 unique isolates and four clusters containing 13 strains. MIRUs discriminated all 13 strains, whereas spoligotyping discriminated 11 strains. Our results validate the use of PCR-based molecular typing of M. tuberculosis using repetitive elements in Indian isolates and demonstrate the usefulness of MIRUs for discriminating low-IS6110-copy isolates, which accounted for more than one-fifth of the strains in the present study.

The presence of Helicobacter pylori in the liver depends on the Th1, Th17 and Treg cytokine profile of the patient

The hypothesis that Helicobactermight be a risk factor for human liver diseases has arisen after the detection of Helicobacter DNA in hepatic tissue of patients with hepatobiliary diseases. Nevertheless, no explanation that justifies the presence of the bacterium in the human liver has been proposed. We evaluated the presence of Helicobacterin the liver of patients with hepatic diseases of different aetiologies. We prospectively evaluated 147 patients (106 with primary hepatic diseases and 41 with hepatic metastatic tumours) and 20 liver donors as controls. Helicobacter species were investigated in the liver by culture and specific 16S rDNA nested-polymerase chain reaction followed by sequencing. Serum and hepatic levels of representative cytokines of T regulatory cell, T helper (Th)1 and Th17 cell lineages were determined using enzyme linked immunosorbent assay. The data were evaluated using logistic models. Detection of Helicobacter pylori DNA in the liver was independently associated with hepatitis B virus/hepatitis C virus, pancreatic carcinoma and a cytokine pattern characterised by high interleukin (IL)-10, low/absent interferon-γ and decreased IL-17A concentrations (p < 10-3). The bacterial DNA was never detected in the liver of patients with alcoholic cirrhosis and autoimmune hepatitis that are associated with Th1/Th17 polarisation. H. pylori may be observed in the liver of patients with certain hepatic and pancreatic diseases, but this might depend on the patient cytokine profile.

Biological, biochemical and molecular features of Trypanosoma cruzi strains isolated from patients infected through oral transmission during a 2005 outbreak in the state of Santa Catarina, Brazil: its correspondence with the new T. cruzi Taxonomy Consensus (2009)

We examined strains of Trypanosoma cruzi isolated from patients with acute Chagas disease that had been acquired by oral transmission in the state of Santa Catarina, Brazil (2005) and two isolates that had been obtained from a marsupial (Didelphis aurita) and a vector (Triatoma tibiamaculata). These strains were characterised through their biological behaviour and isoenzymic profiles and genotyped according to the new Taxonomy Consensus (2009) based on the discrete typing unities, that is, T. cruzi genotypes I-VI. All strains exhibited the biological behaviour of biodeme type II. In six isolates, late peaks of parasitaemia, beyond the 20th day, suggested a double infection with biodemes II + III. Isoenzymes revealed Z2 or mixed Z1 and Z2 profiles. Genotyping was performed using three polymorphic genes (cytochrome oxidase II, spliced leader intergenic region and 24Sα rRNA) and the restriction fragment length polymorphism of the kDNA minicircles. Based on these markers, all but four isolates were characterised as T. cruzi II genotypes. Four mixed populations were identified: SC90, SC93 and SC97 (T. cruzi I + T. cruzi II) and SC95 (T. cruzi I + T. cruzi VI). Comparison of the results obtained by different methods was essential for the correct identification of the mixed populations and major lineages involved indicating that characterisation by different methods can provide new insights into the relationship between phenotypic and genotypic aspects of parasite behaviour.

COI barcode versus morphological identification of Culex ( Culex ) (Diptera: Culicidae) species: a case study using samples from Argentina and Brazil

Sequences of the cytochrome c oxidase subunit I (COI) mitochondrial gene from adults of 22 Culex ( Culex ) species from Argentina and Brazil were employed to assess species identification and to test the usefulness of COI for barcoding using the best close match (BCM) algorithm. A pairwise Kimura two-parameter distance matrix including the mean intra and interspecific distances for 71 COI barcode sequences was constructed. Of the 12 COI lineages recovered in the Neighbour-joining topology, five confirmed recognised morphological species ( Cx. acharistus , Cx. chidesteri , Cx. dolosus , Cx. lygrus and Cx. saltanensis ) with intraspecific divergences lower than 1.75%. Cx. bilineatus is formally resurrected from the synonymy of Cx. dolosus . Cx. maxi , Cx. surinamensis and the Coronator group species included were clustered into an unresolved lineage. The intraspecific distance of Cx. pipiens (3%) was almost twice the interspecific between it and Cx. quinquefasciatus (1.6%). Regarding the BCM criteria, the COI barcode successfully identified 69% of all species. The rest of the sequences, approximately 10%, 18% and 3%, remained as ambiguously, mis and unidentified, respectively. The COI barcode does not contain enough information to distinguish Culex ( Cux. ) species.

Genotyping Mycobacterium bovis from cattle in the Central Pampas of Argentina: temporal and regional trends

Mycobacterium bovis is the causative agent of bovine tuberculosis (TB), a disease that affects approximately 5% of Argentinean cattle. Among the molecular methods for genotyping, the most convenient are spoligotyping and variable number of tandem repeats (VNTR). A total of 378 samples from bovines with visible lesions consistent with TB were collected at slaughterhouses in three provinces, yielding 265 M. bovis spoligotyped isolates, which were distributed into 35 spoligotypes. In addition, 197 isolates were also typed by the VNTR method and 54 combined VNTR types were detected. There were 24 clusters and 27 orphan types. When both typing methods were combined, 98 spoligotypes and VNTR types were observed with 27 clusters and 71 orphan types. By performing a meta-analysis with previous spoligotyping results, we identified regional and temporal trends in the population structure of M. bovis. For SB0140, the most predominant spoligotype in Argentina, the prevalence percentage remained high during different periods, varying from 25.5-57.8% (1994-2011). By contrast, the second and third most prevalent spoligotypes exhibited important fluctuations. This study shows that there has been an expansion in ancestral lineages as demonstrated by spoligotyping. However, exact tandem repeat typing suggests dynamic changes in the clonal population of this microorganism.

Antibody recognition of Plasmodium falciparum infected red blood cells by symptomatic and asymptomatic individuals in the Brazilian Amazon

In the Amazon Region, there is a virtual absence of severe malaria and few fatal cases of naturally occurring Plasmodium falciparum infections; this presents an intriguing and underexplored area of research. In addition to the rapid access of infected persons to effective treatment, one cause of this phenomenon might be the recognition of cytoadherent variant proteins on the infected red blood cell (IRBC) surface, including the var gene encoded P. falciparum erythrocyte membrane protein 1. In order to establish a link between cytoadherence, IRBC surface antibody recognition and the presence or absence of malaria symptoms, we phenotype-selected four Amazonian P. falciparum isolates and the laboratory strain 3D7 for their cytoadherence to CD36 and ICAM1 expressed on CHO cells. We then mapped the dominantly expressed var transcripts and tested whether antibodies from symptomatic or asymptomatic infections showed a differential recognition of the IRBC surface. As controls, the 3D7 lineages expressing severe disease-associated phenotypes were used. We showed that there was no profound difference between the frequency and intensity of antibody recognition of the IRBC-exposed P. falciparum proteins in symptomatic vs. asymptomatic infections. The 3D7 lineages, which expressed severe malaria-associated phenotypes, were strongly recognised by most, but not all plasmas, meaning that the recognition of these phenotypes is frequent in asymptomatic carriers, but is not necessarily a prerequisite to staying free of symptoms.

Behavioural Flexibility: evolutionary past and its role in a changing world

A través de la historia de la vida, gran parte de los organismos han desarrollado estrategias para responder a un mundo en constante cambio. Hoy en día, las actividades humanas producen cambios ambientales a una velocidad sin precedentes, lo cual se traduce en grandes desafíos para la persistencia de biodiversidad. Esta investigación evalúa las respuesta de los animales a los cambios ambientales enfocándose en la flexibilidad del comportamiento como estrategia adaptativa. En una primera aproximación a una escala evolutiva, se otorgan evidencias del vínculo hasta ahora tenue entre la cognición e historias de vida, entregando un claro apoyo a la relación entre longevidad, vida reproductiva y el tamaño del cerebro en mamíferos. La longevidad es el centro de muchas hipótesis respecto a las ventajas de desarrollar un cerebro grande, como por ejemplo en la hipótesis del buffer cognitivo y las respuestas flexibles frente a nuevos ambientes. En un segundo nivel, se abordan factores extrínsecos e intrínsecos que podrían explicar las diferencias individuales en innovación, un componente clave en la flexibilidad del comportamiento. Por medio de una aproximación experimental, se evalúan potenciales escenarios que podrían conducir a consistentes diferencias individuales en uno de los principales factores subyacentes a la innovación (i.e. la motivación), y el potencial control endocrino sobre estos escenarios. Posteriormente, con el objetivo de evaluar la respuesta de los animales frente a los cambios ambientales actuales, se explora la respuesta de los animales frente a una de las actividades humanas mas disruptivas sobre los ecosistemas, la urbanización. Por medio de un analisis filogenetico comparativo a nivel global en aves se abordan los mecanismos implicados en la perdida de biodiversidad observada en ambientes urbanos. Los resultados entregan evidencias sobre la importancia de procesos de dispersión local junto con el papel clave de los rasgos de historia de vida, pero en un sentido diferente al clasicamente pensado. Finalmente por medio de una revisión bibliográfica se entregan evidencias teóricas y empíricas que respaldan el rol clave de la flexibilidad del comportamiento en confrontar los desafíos de una vida urbana. La integración de estos resultados muestra cómo el pasado evolutivo contribuye a hacer frente a los retos ambientales actuales, y pone de relieve posibles consecuencias ante un planeta más cambiante que nunca.

Characterization of microsatellite loci and reliable genotyping in a polyploid plant, Mercurialis perennis (Euphorbiaceae).

For many applications in population genetics, codominant simple sequence repeats (SSRs) may have substantial advantages over dominant anonymous markers such as amplified fragment length polymorphisms (AFLPs). In high polyploids, however, allele dosage of SSRs cannot easily be determined and alleles are not easily attributable to potentially diploidized loci. Here, we argue that SSRs may nonetheless be better than AFLPs for polyploid taxa if they are analyzed as effectively dominant markers because they are more reliable and more precise. We describe the transfer of SSRs developed for diploid Mercurialis huetii to the clonal dioecious M. perennis. Primers were tested on a set of 54 male and female plants from natural decaploid populations. Eight of 65 tested loci produced polymorphic fragments. Binary profiles from 4 different scoring routines were used to define multilocus lineages (MLLs). Allowing for fragment differences within 1 MLL, all analyses revealed the same 14 MLLs without conflicting with merigenet, sex, or plot assignment. For semiautomatic scoring, a combination of as few as 2 of the 4 most polymorphic loci resulted in unambiguous discrimination of clones. Our study demonstrates that microsatellite fingerprinting of polyploid plants is a cost efficient and reliable alternative to AFLPs, not least because fewer loci are required than for diploids.

From diabetes to heart disease : studies on dyslipidemia-induced B-cell dysfunction, immunomodulation in heart transplantation, and cardiac stem cell therapy

Diabetes is a growing epidemic with devastating human, social and economic impact. It is associated with significant changes in plasma concentrations of lipoproteins. We tested the hypothesis that lipoproteins modulate the function and survival of insulin-secreting cells. We first detected the presence of several receptors that participate in the binding and processing of plasma lipoproteins and confirmed the internalization of fluorescent LDL and HDL particles in insulin-secreting β-cells. Purified human VLDL and LDL particles reduced insulin mRNA levels and β-cell proliferation, and induced a dose-dependent increase in the rate of apoptosis. In mice lacking the LDL receptor, islets showed a dramatic decrease in LDL uptake and were partially resistant to apoptosis caused by LDL. VLDL-induced apoptosis of β-cells involved caspase-3 cleavage and reduction in levels of the c-Jun N-terminal (JNK) Interacting Protein-1 (IB1/JIP-1). In contrast, the pro-apoptotic signaling of lipoproteins was antagonized by HDL particles or by a small peptide inhibitor of JNK. The protective effects of HDL were mediated, in part, by inhibition of caspase-3 cleavage and activation of the protein kinase Akt/PKB. Heart disease is a major cause of morbidity and mortality among patients with diabetes. When heart failure is refractory to medical therapy and cannot be improved by electrical resynchronization, percutaneous angioplasty or coronary graft bypass surgery, heart transplantation remains a "last resort" therapy. Nevertheless, it is limited by the side effects of immunosuppressive drugs and chronic rejection. Localized expression of immunomodulatory genes in the donor organ can create a state of immune privilege within the graft, and was performed in rodent hearts by infecting cells with an adenovirus encoding indoleamine 2,3-dioxygenase (IDO), the rate-limiting enzyme in the catabolism of tryptophane. Other strategies are based on genetic manipulation of dendritic cells (DCs) with immunosuppressive genes and in vitro exposure of DCs to agents that prevent their maturation by inflammatory cytokines. Finally, we used 5-bromo-2'-deoxyuridine, which is incorporated into DNA and diluted with cell division, to identify long-term label retaining cells in the adult rodent heart. The majority of these cells were positive for the stem cell antigen-1 (Sca-1) and negative for the endothelial precursor marker CD31. They formed cardiospheres in vitro and showed differentiation potential into mesenchymal cell lineages. When cultured in cardiomyogenic differentiation medium, they expressed cardiac-specific genes. Taken together, these data provide evidence of slow-cycling stem cells in the rodent heart. Chronic shortage of donor organs opens the way to cardiac stem cell therapy in humans, although the long way from animal experimentation to routine therapy in patients may still take several years. - Du diabète de type 2 à la maladie coronarienne : trois études sur les dysfonctions de la cellule sécrétrice d'insuline induites par les dyslipidémies, l'immunomodulation dans la transplantation cardiaque, et la thérapie par des cellules souches myocardiques. Le diabète de type 2 a pris les dimensions d'une épidémie, avec des conséquences sociales et économiques dont nous n'avons pas encore pris toute la mesure. La maladie s'accompagne souvent d'une dyslipidémie caractérisée par une hypertriglycéridémie, des taux abaissés de cholestérol HDL, et des concentrations de cholestérol LDL à la limite supérieure de ce qui est considéré comme acceptable. L'hypothèse à la base de cette étude est qu'une modification des taux plasmatiques de lipoprotéines pourrait avoir une influence directe sur la cellule β sécrétrice d'insuline en modifiant sa fonction, sa durée de vie et son taux de régénération. Dans un premier temps, nous avons mis en évidence, sur la cellule β, la présence de plusieurs récepteurs impliqués dans la captation des lipoprotéines. Nous avons confirmé la fonctionnalité de ces récepteurs en suivant l'internalisation de LDL et de HDL marqués. En présence de VLDL ou de LDL humains, nous avons observé une diminution de la transcription du gène de l'insuline, une prolifération cellulaire réduite, et une augmentation de l'apoptose, toutes fonctions de la dose et du temps d'exposition. L'apoptose induite par les VLDL passe par une activation de la caspase-3 et une réduction du taux de la protéine IB1/JIP-1 (Islet Brain1/JNK Interacting Protein 1), dont une mutation est associée à une forme monogénique de diabète de type 2. Par opposition, les HDL, ainsi que des peptides inhibiteurs de JNK, sont capables de contrer la cascade pro-apoptotique déclenchée, respectivement, par les LDL et les VLDL. Ces effets protecteurs comprennent l'inhibition du clivage de la caspase-3 et l'activation de la protéine kinase Akt/PKB. En conclusion, les lipoprotéines sont des éléments clés de la survie de la cellule β, et pourraient contribuer au dysfonctionnement observé dans le pancréas endocrine au cours du développement du diabète. La maladie cardiaque, et plus particulièrement la maladie coronarienne, est une cause majeure de morbidité et de mortalité chez les patients atteints de diabète. Plusieurs stratégies sont utilisées quotidiennement pour pallier les atteintes cardiaques: traitements médicamenteux, électromécaniques par resynchronisation électrique, ou communément appelés « interventionnels » lorsqu'ils font appel à l'angioplastie percutanée. La revascularisation du myocarde par des pontages coronariens donne également de très bons résultats dans certaines situations. Il existe toutefois des cas où plus aucune de ces approches n'est suffisante. La transplantation cardiaque est alors la thérapie de choix pour un nombre restreint de patients. La thérapie génique, en permettant l'expression locale de gènes immunomodulateurs dans l'organe greffé, permet de diminuer les réactions de rejet inhérentes à toute transplantation (à l'exception de celles réalisées entre deux jumeaux homozygotes). Nous avons appliqué chez des rongeurs cette stratégie en infectant le coeur greffé avec un adénovirus codant pour l'enzyme indoleamine 2,3-dioxygénase (IDO), une enzyme clé dans le catabolisme du tryptophane. Nous avons procédé de manière identique in vitro en surexprimant IDO dans les cellules dendritiques, dont le rôle est de présenter les antigènes aux lymphocytes Τ du receveur. Des expériences similaires ont été réalisées en traitant les cellules dendritiques avec des substances capables de prévenir, en partie du moins, leur maturation par des agents pro-inflammatoires. Finalement, nous avons exploré une stratégie utilisée couramment en hématologie, mais qui n'en est encore qu'à ses débuts au niveau cardiaque : la thérapie par des cellules souches. En traitant des rongeurs avec un marqueur qui s'incorpore dans l'ADN nucléaire, le 5-bromo- 2'-deoxyuridine, nous avons identifié une population cellulaire se divisant rarement, positive en grande partie pour l'antigène embryonnaire Sca-1 et négative pour le marqueur endothélial CD31. En culture, ces cellules forment des cardiosphères et sont capables de se différencier dans les principaux types tissulaires mésenchymateux. Dans un milieu de differentiation adéquat, ces cellules expriment des gènes cardiomyocytaires. En résumé, ces données confirment la présence chez le rongeur d'une population résidente de précurseurs myocardiques. En addenda, on trouvera deux publications relatives à la cellule β productrice d'insuline. Le premier article démontre le rôle essentiel joué par la complexine dans l'insulino-sécrétion, tandis que le second souligne l'importance de la protéine IB1/JIP-1 dans la protection contre l'apoptose de la cellule β induite par certaines cytokines.

Using parthenogenetic lineages to identify advantages of sex

The overwhelming predominance of sexual reproduction in nature is surprising given that sex is expected to confer profound costs in terms of production of males and the breakup of beneficial allele combinations. Recognition of these theoretical costs was the inspiration for a large body of empirical research-typically focused on comparing sexual and asexual organisms, lineages, or genomes-dedicated to identifying the advantages and maintenance of sex in natural populations. Despite these efforts, why sex is so common remains unclear. Here, we argue that we can generate general insights into the advantages of sex by taking advantage of parthenogenetic taxa that differ in such characteristics as meiotic versus mitotic offspring production, ploidy level, and single versus multiple and hybrid versus non-hybrid origin. We begin by evaluating benefits that sex can confer via its effects on genetic linkage, diversity, and heterozygosity and outline how the three classes of benefits make different predictions for which type of parthenogenetic lineage would be favored over others. Next, we describe the type of parthenogenetic model system (if any) suitable for testing whether the hypothesized benefit might contribute to the maintenance of sex in natural populations, and suggest groups of organisms that fit the specifications. We conclude by discussing how empirical estimates of characteristics such as time since derivation and number of independent origins of asexual lineages from sexual ancestors, ploidy levels, and patterns of molecular evolution from representatives of these groups can be used to better understand which mechanisms maintain sex in natural populations.

Vitellogenin Underwent Subfunctionalization to Acquire Caste and Behavioral Specific Expression in the Harvester Ant Pogonomyrmex barbatus.

The reproductive ground plan hypothesis (RGPH) proposes that the physiological pathways regulating reproduction were co-opted to regulate worker division of labor. Support for this hypothesis in honeybees is provided by studies demonstrating that the reproductive potential of workers, assessed by the levels of vitellogenin (Vg), is linked to task performance. Interestingly, contrary to honeybees that have a single Vg ortholog and potentially fertile nurses, the genome of the harvester ant Pogonomyrmex barbatus harbors two Vg genes (Pb_Vg1 and Pb_Vg2) and nurses produce infertile trophic eggs. P. barbatus, thus, provides a unique model to investigate whether Vg duplication in ants was followed by subfunctionalization to acquire reproductive and non-reproductive functions and whether Vg reproductive function was co-opted to regulate behavior in sterile workers. To investigate these questions, we compared the expression patterns of P. barbatus Vg genes and analyzed the phylogenetic relationships and molecular evolution of Vg genes in ants. qRT-PCRs revealed that Pb_Vg1 is more highly expressed in queens compared to workers and in nurses compared to foragers. By contrast, the level of expression of Pb_Vg2 was higher in foragers than in nurses and queens. Phylogenetic analyses show that a first duplication of the ancestral Vg gene occurred after the divergence between the poneroid and formicoid clades and subsequent duplications occurred in the lineages leading to Solenopsis invicta, Linepithema humile and Acromyrmex echinatior. The initial duplication resulted in two Vg gene subfamilies preferentially expressed in queens and nurses (subfamily A) or in foraging workers (subfamily B). Finally, molecular evolution analyses show that the subfamily A experienced positive selection, while the subfamily B showed overall relaxation of purifying selection. Our results suggest that in P. barbatus the Vg gene underwent subfunctionalization after duplication to acquire caste- and behavior- specific expression associated with reproductive and non-reproductive functions, supporting the validity of the RGPH in ants.

Perturbation expansions of multilocus fixation probabilities for frequency-dependent selection with applications to the Hill-Robertson effect and to the joint evolution of helping and punishment.

Natural populations are of finite size and organisms carry multilocus genotypes. There are, nevertheless, few results on multilocus models when both random genetic drift and natural selection affect the evolutionary dynamics. In this paper we describe a formalism to calculate systematic perturbation expansions of moments of allelic states around neutrality in populations of constant size. This allows us to evaluate multilocus fixation probabilities (long-term limits of the moments) under arbitrary strength of selection and gene action. We show that such fixation probabilities can be expressed in terms of selection coefficients weighted by mean first passages times of ancestral gene lineages within a single ancestor. These passage times extend the coalescence times that weight selection coefficients in one-locus perturbation formulas for fixation probabilities. We then apply these results to investigate the Hill-Robertson effect and the coevolution of helping and punishment. Finally, we discuss limitations and strengths of the perturbation approach. In particular, it provides accurate approximations for fixation probabilities for weak selection regimes only (Ns < or = 1), but it provides generally good prediction for the direction of selection under frequency-dependent selection.