Extensive remodeling of DC function by rapid maturation-induced transcriptional silencing.

The activation, or maturation, of dendritic cells (DCs) is crucial for the initiation of adaptive T-cell mediated immune responses. Research on the molecular mechanisms implicated in DC maturation has focused primarily on inducible gene-expression events promoting the acquisition of new functions, such as cytokine production and enhanced T-cell-stimulatory capacity. In contrast, mechanisms that modulate DC function by inducing widespread gene-silencing remain poorly understood. Yet the termination of key functions is known to be critical for the function of activated DCs. Genome-wide analysis of activation-induced histone deacetylation, combined with genome-wide quantification of activation-induced silencing of nascent transcription, led us to identify a novel inducible transcriptional-repression pathway that makes major contributions to the DC-maturation process. This silencing response is a rapid primary event distinct from repression mechanisms known to operate at later stages of DC maturation. The repressed genes function in pivotal processes--including antigen-presentation, extracellular signal detection, intracellular signal transduction and lipid-mediator biosynthesis--underscoring the central contribution of the silencing mechanism to rapid reshaping of DC function. Interestingly, promoters of the repressed genes exhibit a surprisingly high frequency of PU.1-occupied sites, suggesting a novel role for this lineage-specific transcription factor in marking genes poised for inducible repression.

Evolutionary histories of expanded peptidase families in Schistosoma mansoni

Schistosoma mansoni is one of the three main causative agents of human schistosomiasis, a major health problem with a vast socio-economic impact. Recent advances in the proteomic analysis of schistosomes have revealed that peptidases are the main virulence factors involved in the pathogenesis of this disease. In this context, evolutionary studies can be applied to identify peptidase families that have been expanded in genomes over time in response to different selection pressures. Using a phylogenomic approach, we searched for expanded endopeptidase families in the S. mansoni predicted proteome with the aim of contributing to the knowledge of such enzymes as potential therapeutic targets. We found three endopeptidase families that comprise leishmanolysins (metallopeptidase M8 family), cercarial elastases (serine peptidase S1 family) and cathepsin D proteins (aspartic peptidase A1 family). Our results suggest that the Schistosoma members of these families originated from successive gene duplication events in the parasite lineage after its diversification from other metazoans. Overall, critical residues are conserved among the duplicated genes/proteins. Furthermore, each protein family displays a distinct evolutionary history. Altogether, this work provides an evolutionary view of three S. mansoni peptidase families, which allows for a deeper understanding of the genomic complexity and lineage-specific adaptations potentially related to the parasitic lifestyle.

Developmentally regulated expression of P-glycoprotein (multidrug resistance) activity in mouse thymocytes.

P-glycoprotein (P-gly) is the transmembrane efflux pump responsible for multidrug resistance in tumor cells. The activity of P-gly in mature peripheral lymphocytes is lineage specific, with CD8+ T cells and natural killer (NK) cells expressing high levels as compared to CD4+ T cells and B cells. We have now investigated P-gly activity in immature and mature subsets of mouse thymocytes. Our data indicate that P-gly activity is undetectable in immature CD4-8- and CD4+8+ thymocyte subsets. Among mature thymocytes, P-gly activity is absent in the CD4+ subset but present in the more mature (HSAlow) fraction of CD8+ cells. Furthermore, while thymic CD4-8- T cell receptor (TCR) gamma delta cells have little P-gly activity, a minor subset of CD4-8- or CD4+ TCR alpha beta + thymocytes bearing the NK1.1 surface marker expresses high levels of P-gly activity. Collectively, our results indicate that P-gly activity arises late during thymus development and is expressed in a lineage-specific fashion.

Molecular changes underlying mammalian phenotypic innovation

Mammals are characterized by specific phenotypic traits that include lactation, hair, and relatively large brains with unique structures. Individual mammalian lineages have, in turn, evolved characteristic traits that distinguish them from others. These include obvious anatom¬ical differences but also differences related to reproduction, life span, cognitive abilities, be¬havior. and disease susceptibility. However, the molecular basis of the diverse mammalian phenotypes and the selective pressures that shaped their evolution remain largely unknown. In the first part of my thesis, I analyzed the genetic factors associated with the origin of a unique mammalian phenotype lactation and I studied the selective pressures that forged the transition from oviparity to viviparity. Using a comparative genomics approach and evolutionary simulations, I showed that the emergence of lactation, as well as the appear¬ance of the casein gene family, significantly reduced selective pressure on the major egg-yolk proteins (the vitellogenin family). This led to a progressive loss of vitellogenins, which - in oviparous species - act as storage proteins for lipids, amino acids, phosphorous and calcium in the isolated egg. The passage to internal fertilization and placentation in therian mam¬mals rendered vitellogenins completely dispensable, which ended in the loss of the whole gene family in this lineage. As illustrated by the vitellogenin study, changes in gene content are one possible underlying factor for the evolution of mammalian-specific phenotypes. However, more subtle genomic changes, such as mutations in protein-coding sequences, can also greatly affect the phenotypes. In particular, it was proposed that changes at the level of gene reg¬ulation could underlie many (or even most) phenotypic differences between species. In the second part of my thesis, I participated in a major comparative study of mammalian tissue transcriptomes, with the goal of understanding how evolutionary forces affected expression patterns in the past 200 million years of mammalian evolution. I showed that, while com¬parisons of gene expressions are in agreement with the known species phylogeny, the rate of expression evolution varies greatly among lineages. Species with low effective population size, such as monotremes and hominoids, showed significantly accelerated rates of gene expression evolution. The most likely explanation for the high rate of gene expression evolution in these lineages is the accumulation of mildly deleterious mutations in regulatory regions, due to the low efficiency of purifying selection. Thus, our observations are in agreement with the nearly neutral theory of molecular evolution. I also describe substantial differences in evolutionary rates between tissues, with brain being the most constrained (especially in primates) and testis significantly accelerated. The rate of gene expression evolution also varies significantly between chromosomes. In particular, I observed an acceleration of gene expression changes on the X chromosome, probably as a result of adaptive processes associated with the origin of therian sex chromosomes. Lastly, I identified several individual genes as well as co-regulated expression modules that have undergone lineage specific expression changes and likely under¬lie various phenotypic innovations in mammals. The methods developed during my thesis, as well as the comprehensive gene content analyses and transcriptomics datasets made available by our group, will likely prove to be useful for further exploratory analyses of the diverse mammalian phenotypes.

Patterns of evolutionary constraints on genes in humans

Background: Different regions in a genome evolve at different rates depending on structural and functional constraints. Some genomic regions are highly conserved during metazoan evolution, while other regions may evolve rapidly, either in all species or in a lineage-specific manner. A strong or even moderate change in constraints in functional regions, for example in coding regions, can have significant evolutionary consequences. Results: Here we discuss a novel framework, 'BaseDiver', to classify groups of genes in humans based on the patterns of evolutionary constraints on polymorphic positions in their coding regions. Comparing the nucleotide-level divergence among mammals with the extent of deviation from the ancestral base in the human lineage, we identify patterns of evolutionary pressure on nonsynonymous base-positions in groups of genes belonging to the same functional category. Focussing on groups of genes in functional categories, we find that transcription factors contain a significant excess of nonsynonymous base-positions that are conserved in other mammals but changed in human, while immunity related genes harbour mutations at base-positions that evolve rapidly in all mammals including humans due to strong preference for advantageous alleles. Genes involved in olfaction also evolve rapidly in all mammals, and in humans this appears to be due to weak negative selection. Conclusion: While recent studies have identified genes under positive selection in humans, our approach identifies evolutionary constraints on Gene Ontology groups identifying changes in humans relative to some of the other mammals.

ExPASy: SIB bioinformatics resource portal.

ExPASy (http://www.expasy.org) has worldwide reputation as one of the main bioinformatics resources for proteomics. It has now evolved, becoming an extensible and integrative portal accessing many scientific resources, databases and software tools in different areas of life sciences. Scientists can henceforth access seamlessly a wide range of resources in many different domains, such as proteomics, genomics, phylogeny/evolution, systems biology, population genetics, transcriptomics, etc. The individual resources (databases, web-based and downloadable software tools) are hosted in a 'decentralized' way by different groups of the SIB Swiss Institute of Bioinformatics and partner institutions. Specifically, a single web portal provides a common entry point to a wide range of resources developed and operated by different SIB groups and external institutions. The portal features a search function across 'selected' resources. Additionally, the availability and usage of resources are monitored. The portal is aimed for both expert users and people who are not familiar with a specific domain in life sciences. The new web interface provides, in particular, visual guidance for newcomers to ExPASy.

Trophic specialization influences the rate of environmental niche evolution in damselfishes (Pomacentridae).

The rate of environmental niche evolution describes the capability of species to explore the available environmental space and is known to vary among species owing to lineage-specific factors. Trophic specialization is a main force driving species evolution and is responsible for classical examples of adaptive radiations in fishes. We investigate the effect of trophic specialization on the rate of environmental niche evolution in the damselfish, Pomacentridae, which is an important family of tropical reef fishes. First, phylogenetic niche conservatism is not detected in the family using a standard test of phylogenetic signal, and we demonstrate that the environmental niches of damselfishes that differ in trophic specialization are not equivalent while they still overlap at their mean values. Second, we estimate the relative rates of niche evolution on the phylogenetic tree and show the heterogeneity among rates of environmental niche evolution of the three trophic groups. We suggest that behavioural characteristics related to trophic specialization can constrain the evolution of the environmental niche and lead to conserved niches in specialist lineages. Our results show the extent of influence of several traits on the evolution of the environmental niche and shed new light on the evolution of damselfishes, which is a key lineage in current efforts to conserve biodiversity in coral reefs.

A burst of segmental duplications in the genome of the African great ape ancestor

It is generally accepted that the extent of phenotypic change between human and great apes is dissonant with the rate of molecular change. Between these two groups, proteins are virtually identical, cytogenetically there are few rearrangements that distinguish ape-human chromosomes, and rates of single-base-pair change and retrotransposon activity have slowed particularly within hominid lineages when compared to rodents or monkeys. Studies of gene family evolution indicate that gene loss and gain are enriched within the primate lineage. Here, we perform a systematic analysis of duplication content of four primate genomes (macaque, orang-utan, chimpanzee and human) in an effort to understand the pattern and rates of genomic duplication during hominid evolution. We find that the ancestral branch leading to human and African great apes shows the most significant increase in duplication activity both in terms of base pairs and in terms of events. This duplication acceleration within the ancestral species is significant when compared to lineage-specific rate estimates even after accounting for copy-number polymorphism and homoplasy. We discover striking examples of recurrent and independent gene-containing duplications within the gorilla and chimpanzee that are absent in the human lineage. Our results suggest that the evolutionary properties of copy-number mutation differ significantly from other forms of genetic mutation and, in contrast to the hominid slowdown of single-base-pair mutations, there has been a genomic burst of duplication activity at this period during human evolution.

Global Genomic Diversity of Human Papillomavirus 6 Based on 724 Isolates and 190 Complete Genome Sequences.

Human papillomavirus type 6 (HPV6) is the major etiological agent of anogenital warts and laryngeal papillomas and has been included in both the quadrivalent and nonavalent prophylactic HPV vaccines. This study investigated the global genomic diversity of HPV6, using 724 isolates and 190 complete genomes from six continents, and the association of HPV6 genomic variants with geographical location, anatomical site of infection/disease, and gender. Initially, a 2,800-bp E5a-E5b-L1-LCR fragment was sequenced from 492/530 (92.8%) HPV6-positive samples collected for this study. Among them, 130 exhibited at least one single nucleotide polymorphism (SNP), indel, or amino acid change in the E5a-E5b-L1-LCR fragment and were sequenced in full. A global alignment and maximum likelihood tree of 190 complete HPV6 genomes (130 fully sequenced in this study and 60 obtained from sequence repositories) revealed two variant lineages, A and B, and five B sublineages: B1, B2, B3, B4, and B5. HPV6 (sub)lineage-specific SNPs and a 960-bp representative region for whole-genome-based phylogenetic clustering within the L2 open reading frame were identified. Multivariate logistic regression analysis revealed that lineage B predominated globally. Sublineage B3 was more common in Africa and North and South America, and lineage A was more common in Asia. Sublineages B1 and B3 were associated with anogenital infections, indicating a potential lesion-specific predilection of some HPV6 sublineages. Females had higher odds for infection with sublineage B3 than males. In conclusion, a global HPV6 phylogenetic analysis revealed the existence of two variant lineages and five sublineages, showing some degree of ethnogeographic, gender, and/or disease predilection in their distribution. IMPORTANCE: This study established the largest database of globally circulating HPV6 genomic variants and contributed a total of 130 new, complete HPV6 genome sequences to available sequence repositories. Two HPV6 variant lineages and five sublineages were identified and showed some degree of association with geographical location, anatomical site of infection/disease, and/or gender. We additionally identified several HPV6 lineage- and sublineage-specific SNPs to facilitate the identification of HPV6 variants and determined a representative region within the L2 gene that is suitable for HPV6 whole-genome-based phylogenetic analysis. This study complements and significantly expands the current knowledge of HPV6 genetic diversity and forms a comprehensive basis for future epidemiological, evolutionary, functional, pathogenicity, vaccination, and molecular assay development studies.

Characterising avian haemosporidian communities : ecological factors, parasite persistence and co-infection

Les parasites jouent un rôle clef dans l'évolution des comportements et des traits d'histoire de vie de leurs hôtes. Le parasitisme s'avère parfois dévastateur à l'échelle de population d'hôtes, et peut également altérer certains traits associés à la valeur sélective d'un individu infecté, tels que son succès reproducteur ou encore son taux de mortalité. La coévolution hôte/parasite, qui représente l'une des forces sélectives les plus puissantes dans l'évolution des organismes, peut également conduire les partenaires de l'association parasitaire à s'adapter localement à des environnements hétérogènes. Cette thèse porte sur l'étude de parasites aviaires, du genre Plasmodium, Haemopro- teus et Leucocytozoon (Haemosporidae), naturellement associés à différentes populations de mésanges charbonnières (Parus major) et d'hirondelles des fenêtres (Delichon ur- bicum). Dans un premier temps, nous avons cherché à déterminer comment se distribuent ces parasites au sein de différentes populations hôtes et si ces communautés de parasites sont structurées. Par la suite, la principale question à laquelle nous voulions répondre était de savoir comment ces parasites, et notamment après coexistence de plusieurs lignées génétiques d'Haemosporidae au sein dun même-individu (i.e. co-infection), affectent la physiologie et le succès de reproducteur des hôtes. Nos résultats suggèrent que la distribution des Haemosporidae est principalement gouvernée par la présence d'insectes vecteurs et que la persistance de l'infection chez les hôtes varie en fonction du genre d'Haemosporidae (Chapitre 1-2). Par ailleurs, nous avons trouvé que des lignées de parasite génétiquement distinctes peuvent avoir des effets contrastés sur leurs hôtes. Par exemple, les hôtes exhibent des différences de parasitémie marquées en fonction des lignées de parasites responsable de l'infection. De plus, le succès reproducteur ainsi que la charge parasitaire des mésanges infectées par Plasmodium ou Haemoproteus n'étaient pas affecté par l'infection simultanée avec Leucocytozoon (Chapitre 2-3). Dans le Chapitre 4, j'ai examiné la capacité immunitaire de mésanges charbonnières infectées par des hémosporidies. Les résultats n'ont pas été concluant, et je suggère fortement une réévaluation de ceux-ci dans de futures études. Les mésanges charbonnières ne semblent pas signaler leur statut infectieux par la coloration de leur plumage (Chapitre 5); toutefois, la coloration noire des plumes reflète l'état de stress oxydatif des mésanges, qui dépend lui-même de l'infection parasitaire. La coloration verte pourrait également indiquer la qualité des soins paxentaux délivrés par les mésanges adultes femelles à leurs petits, comme le suggère la corrélation que nous avons observée entre la masse des jeunes d'une nichée et la coloration de leur mère. Les hirondelles capturées en Algérie souffrent plus de l'infection que celles échantillon¬nées en Europe (Chapitre 6). Les similitudes observées entre les communautés de par¬asites affectant les populations européennes et celles des populations nord-africaines suggèrent que la transmission des parasites a lieu lors de la migration vers le sud. A l'instar de nos observations sur les mésanges dans les chapitres 2 et 3, les hirondelles co-infectées ne montrent pas d'altérations de leur condition physique. Cette thèse démontre qu'il existe, au sein des populations de mésanges charbonnières, des interactions antagonistes entre, d'une part, les parasites et leurs hôtes et d'autre part, entre différent parasites. Le résultat de ces interactions antagonistes varie en fonction des espèces et de la zone géographique considérée. Nous avons démontré que les interactions ne suivent pas toujours la théorie, puisque la coevolution qui, en suivant le concept de la virulence, devrait augmenter la charge parasitaire et diminuer la condition physique des hôtes, ne montre pourtant pas d'impact négatif sur les populations de mésanges. Nous pouvons maintenant concentrer nos efforts à la caractérisation des interactions antagonistes. De plus, grâce aux avancées des méthodes moléculaires, nous pouvons suivre et étudier en détails comment ces interactions se manifestent et quels sont leurs effets sur la condition physique des hôtes. - Parasites are key in shaping various behavioural and life-history traits of their hosts. The influence of parasitism on host populations varies from slight to devastating and might influence such parameters as mortality rates or reproductive success. Host-parasite coevolution is one of the most powerful selective forces in evolution and can lead to local adaptation of parasites and hosts in spatially structured environments. In this thesis, I studied haemosporidian parasites in different populations of great tits (Parus major) and house martins (Delichon urbicum). Firstly, I wanted to determine how parasites are distributed and if parasite communities are structured. The main question I wanted to address hereafter was how parasites, and specifically infection with multiple genera of parasites (i.e. co-infection) influenced host physiology and reproductive success. I found that parasite distribution is environmentally driven and could therefore be closely linked to vector prevalence; and that the stability of parasite infection over time is genus-dependent (Chapter 1 - 2). I further found that different haemosporidian lineages might interact differently with their hosts as parasitaemia was strongly lineage-specific and that the presence of Leucocytozoon parasites showed no correlation to Plasmodium or Haemoproteus parasitaemia, nor to great tit reproductive success (Chapter 2-3). In Chapter 4 I examined immune capacity of haemosporidian-infected great tits. The results proved inconclusive, and I strongly suggest re-evaluation hereof in future work. Great tits do not appear to signal parasite infection through plumage colouration (Chapter 5); however, infection did have a link to oxidative stress resistance which is strongly signalled through the black breast stripe, with darker males being more resistant and darker females less resistant. Females might incur different costs associated with darker stripes. This would allow reversal of signaling function. Green colouration could also serve as a cue for female provisioning quality as indicated by the strong correlation between colouration and chick body mass. Breeding house martins caught in Algeria suffer greater haemosporidian infection than European populations (Chapter 6). Similar parasite communities in European and North-African populations suggest transmission of parasites may occur during southward migration. Similarly to what was observed in great tits in Chapter 2 and 3, no relationship was found between parasite co-infection and Swiss house martin body condition. This thesis demonstrates that host-parasite and inter-parasite antagonistic interac¬tions exist in great tit populations. How these interactions play out is species dependent and varies geographically. I have demonstrated that interactions do not always follow the theory, as co-infection - which under the concept of virulence should increase parasitaemia and decrease body condition - showed no negative impact on great tit populations. We can now concentrate our efforts on characterising these antagonistic interactions, and with the advance in molecular methods, track and investigate how these interactions play out and what the effect on host fitness is.

Staphylococcus aureus host range and human-bovine host shift.

Staphylococcus aureus is a major agent of bovine mastitis. The concomitant emergence of pig-associated methicillin-resistant S. aureus (MRSA) in human carriage and infection requires a reexamination of the host range and specificity of human- and cow-associated S. aureus strains, something which has not been systematically studied previously. The genetic relatedness of 500 S. aureus isolates from bovine mastitis cases, 57 isolates from nasal carriage of farmers, and 133 isolates from nonfarmers was determined by amplified fragment length polymorphism (AFLP) analysis and spa typing. Multilocus sequence typing (MLST) was conducted on a subset of isolates to match AFLP clusters with MLST clonal complexes (CCs). This data set allowed us to study host range and host specificity and to estimate the extent of bovine-to-human transmission. The genotype compositions of S. aureus isolates from farmers and nonfarmers were very similar, while the mastitis isolates were quite distinct. Overall, transmission was low, but specific genotypes did show increased cow-to-human transmission. Unexpectedly, more than one-third of mastitis isolates belonged to CC8, a lineage which has not been considered to be bovine mastitis associated, but it is well known from human carriage and infection (i.e., USA300). Despite the fact that we did detect some transmission of other genotypes from cows to farmers, no transmission of CC8 isolates to farmers was detected, except for one tentative case. This was despite the close genetic relatedness of mastitis CC8 strains to nonfarmer carriage strains. These results suggest that the emergence of the new bovine-adapted genotype was due to a recent host shift from humans to cows concurrent with a loss of the ability to colonize humans. More broadly, our results indicate that host specificity is a lineage-specific trait that can rapidly evolve.

Avian haemosporidian persistence and co-infection in great tits at the individual level.

BACKGROUND: Many studies have tracked the distribution and persistence of avian haemosporidian communities across space and time at the population level, but few studies have investigated these aspects of infection at the individual level over time. Important aspects of parasite infection at the individual level can be missed if only trends at the population level are studied. This study aimed to determine how persistent Haemosporida are in great tit individuals recaptured over several years, whether parasitaemia differed by parasite lineage (mitochondrial cytochrome b haplotype) and how co-infection (i.e. concurrent infection with multiple genera of parasites) affects parasitaemia and body mass. METHODS: Parasite prevalence was determined by polymerase chain reaction (PCR), quantitative PCR were used to assess parasitaemia and sequencing was employed to determine the identity of the lineages using the MalAvi database. RESULTS: Haemosporidian prevalence was high over sampled years with 98% of 55 recaptured individuals showing infection in at least one year of capture. Eighty-two percent of all positive individuals suffered co-infection, with an overall haemosporidian lineage diversity of seventeen. Plasmodium and Haemoproteus parasites were found to be highly persistent, with lineages from these genera consistently found in individuals across years and with no differences in individual parasitaemia being recorded at subsequent captures. Conversely, Leucocytozoon parasites showed higher turnover with regard to lineage changes or transitions in infection status (infected vs non-infected) across years. Parasitaemia was found to be lineage specific and there was no relationship between Plasmodium parasitaemia or host body condition and the presence of Leucocytozoon parasites. CONCLUSIONS: The findings of this study suggest that different genera of haemosporidian parasites interact differently with their host and other co-infecting parasites, influencing parasite persistence most likely through inter-parasite competition or host-parasite immune interactions. Even-though co-infections do not seem to result in increased virulence (higher parasitaemia or poorer host body condition), further investigation into infection potential of these parasites, both individually and as co-infections, is necessary.

BadiRate: Estimating Family Turnover Rates by Likelihood-Based Methods

Motivation: The comparative analysis of gene gain and loss rates is critical for understanding the role of natural selection and adaptation in shaping gene family sizes. Studying complete genome data from closely related species allows accurate estimation of gene family turnover rates. Current methods and software tools, however, are not well designed for dealing with certain kinds of functional elements, such as microRNAs or transcription factor binding sites. Results: Here, we describe BadiRate, a new software tool to estimate family turnover rates, as well as the number of elements in internal phylogenetic nodes, by likelihood-based methods and parsimony. It implements two stochastic population models, which provide the appropriate statistical framework for testing hypothesis, such as lineage-specific gene family expansions or contractions. We have assessed the accuracy of BadiRate by computer simulations, and have also illustrated its functionality by analyzing a representative empirical dataset.

Odorant binding proteins of the red imported fire ant, Solenopsis invicta: an example of the problems facing the analysis of widely divergent proteins.

We describe the odorant binding proteins (OBPs) of the red imported fire ant, Solenopsis invicta, obtained from analyses of an EST library and separate 454 sequencing runs of two normalized cDNA libraries. We identified a total of 18 putative functional OBPs in this ant. A third of the fire ant OBPs are orthologs to honey bee OBPs. Another third of the OBPs belong to a lineage-specific expansion, which is a common feature of insect OBP evolution. Like other OBPs, the different fire ant OBPs share little sequence similarity (∼ 20%), rendering evolutionary analyses difficult. We discuss the resulting problems with sequence alignment, phylogenetic analysis, and tests of selection. As previously suggested, our results underscore the importance for careful exploration of the sensitivity to the effects of alignment methods for data comprising widely divergent sequences.

BadiRate: estimating family turnover rates by likelihood-based methods

Motivation: The comparative analysis of gene gain and loss rates is critical for understanding the role of natural selection and adaptation in shaping gene family sizes. Studying complete genome data from closely related species allows accurate estimation of gene family turnover rates. Current methods and software tools, however, are not well designed for dealing with certain kinds of functional elements, such as microRNAs or transcription factor binding sites. Results: Here, we describe BadiRate, a new software tool to estimate family turnover rates, as well as the number of elements in internal phylogenetic nodes, by likelihood-based methods and parsimony. It implements two stochastic population models, which provide the appropriate statistical framework for testing hypothesis, such as lineage-specific gene family expansions or contractions. We have assessed the accuracy of BadiRate by computer simulations, and have also illustrated its functionality by analyzing a representative empirical dataset.