Classes Alléliques d’Haplotypes et Sélection Positive dans le Génome Humain

L'identification de régions génomiques cibles de la sélection naturelle positive permet de mieux comprendre notre passé évolutif et de trouver des variants génétiques fonctionnels importants. Puisque la fréquence des allèles sélectionnés augmente dans la population, la sélection laisse des traces sur les séquences d'ADN et ces empreintes sont détectées lorsque la variabilité génétique d'une région est différente de celle attendue sous neutralité sélective. On propose une nouvelle approche pour analyser les données de polymorphismes : le calcul des classes alléliques d’haplotypes (HAC), permettant d'évaluer la diversité globale des haplotypes en étudiant leur composition allélique. L'idée de l'approche est de déterminer si un site est sous sélection positive récente en comparant les distributions des HAC obtenues pour les deux allèles de ce site. Grâce à l'utilisation de données simulées, nous avons étudié ces distributions sous neutralité et sous sélection en testant l'effet de différents paramètres populationnels. Pour tester notre approche empiriquement, nous avons analysé la variation génétique au niveau du gène de lactase dans les trois populations inclues dans le projet HapMap.

"Genetic Diversity and Selection in Three Plasmodium vivax Merozoite Surface Protein 7 (Pvmsp-7) Genes in a Colombian Population"

A completely effective vaccine for malaria (one of the major infectious diseases worldwide) is not yet available; different membrane proteins involved in parasite-host interactions have been proposed as candidates for designing it. It has been found that proteins encoded by the merozoite surface protein (msp)-7 multigene family are antibody targets in natural infection; the nucleotide diversity of three Pvmsp-7 genes was thus analyzed in a Colombian parasite population. By contrast with P. falciparum msp-7 loci and ancestral P. vivax msp-7 genes, specie-specific duplicates of the latter specie display high genetic variability, generated by single nucleotide polymorphisms, repeat regions, and recombination. At least three major allele types are present in Pvmsp-7C, Pvmsp-7H and Pvmsp-7I and positive selection seems to be operating on the central region of these msp-7 genes. Although this region has high genetic polymorphism, the C-terminus (Pfam domain ID: PF12948) is conserved and could be an important candidate when designing a subunit-based antimalarial vaccine.

Genetic diversity and selection in three plasmodium vivax merozoite surface protein 7 (pvmsp-7) genes in a colombian population

A completely effective vaccine for malaria (one of the major infectious diseases worldwide) is not yet available; different membrane proteins involved in parasite-host interactions have been proposed as candidates for designing it. It has been found that proteins encoded by the merozoite surface protein (msp)-7 multigene family are antibody targets in natural infection; the nucleotide diversity of three Pvmsp-7 genes was thus analyzed in a Colombian parasite population. By contrast with P. falciparum msp-7 loci and ancestral P. vivax msp-7 genes, specie-specific duplicates of the latter specie display high genetic variability, generated by single nucleotide polymorphisms, repeat regions, and recombination. At least three major allele types are present in Pvmsp-7C, Pvmsp-7H and Pvmsp-7I and positive selection seems to be operating on the central region of these msp-7 genes. Although this region has high genetic polymorphism, the C-terminus (Pfam domain ID: PF12948) is conserved and could be an important candidate when designing a subunit-based antimalarial vaccine.

Genomic data as the hitchhiker's guide to cattle adaptation: tracking the milestones of past selection in the bovine genome

The bovine species have witnessed and played a major role in the drastic socio-economical changes that shaped our culture over the last 10,000 years. During this journey, cattle hitchhiked on human development and colonized the world, facing strong selective pressures such as dramatic environmental changes and disease challenge. Consequently, hundreds of specialized cattle breeds emerged and spread around the globe, making up a rich spectrum of genomic resources. Their DNA still carry the scars left from adapting to this wide range of conditions, and we are now empowered with data and analytical tools to track the milestones of past selection in their genomes. In this review paper, we provide a summary of the reconstructed demographic events that shaped cattle diversity, offer a critical synthesis of popular methodologies applied to the search for signatures of selection (SS) in genomic data, and give examples of recent SS studies in cattle. Then, we outline the potential and challenges of the application of SS analysis in cattle, and discuss the future directions in this field.

Selection and recombination drive the evolution of MHC class II DRB diversity in ungulates

Major histocompatibility complex (MHC) antigen-presenting genes are the most variable loci in vertebrate genomes. Host-parasite co-evolution is assumed to maintain the excessive polymorphism in the MHC loci. However, the molecular mechanisms underlying the striking diversity in the MHC remain contentious. The extent to which recombination contributes to the diversity at MHC loci in natural populations is still controversial, and there have been only few comparative studies that make quantitative estimates of recombination rates. In this study, we performed a comparative analysis for 15 different ungulates species to estimate the population recombination rate, and to quantify levels of selection. As expected for all species, we observed signatures of strong positive selection, and identified individual residues experiencing selection that were congruent with those constituting the peptide-binding region of the human DRB gene. However, in addition for each species, we also observed recombination rates that were significantly different from zero on the basis of likelihood-permutation tests, and in other non-quantitative analyses. Patterns of synonymous and non-synonymous sequence diversity were consistent with differing demographic histories between species, but recent simulation studies by other authors suggest inference of selection and recombination is likely to be robust to such deviations from standard models. If high rates of recombination are common in MHC genes of other taxa, re-evaluation of many inference-based phylogenetic analyses of MHC loci, such as estimates of the divergence time of alleles and trans-specific polymorphism, may be required.

Expression and diversity of Echinococcus multilocularis AgB genes in secondarily infected mice: evaluating the influence of T-cell immune selection on antigenic variation

The T-cell-mediated immune response exhibits a crucial function in the control of the intrahepatic proliferation of Echinococcus multilocularis larvae in mice and humans, both being natural intermediate hosts of the parasite. Antigen B (AgB), a metabolized Echinococcus spp. lipoprotein, contributes to the modulation of the T-cell immune response, and distinct sites of the corresponding AgB1, AgB3 and AgB4 genes were shown to be under positive selection pressure. Since several AgB gene variants are present in a single Echinococcus metacestode, we used secondary E. multilocularis infections in BALB/c and in athymic nude mice (devoid of T-cell responses) to analyze the effect of the cellular immune response on the expression and diversity of EmAgB1-EmAgB4 genes. We demonstrated hereby that EmAgB transcripts were less abundant in nude mice during the early phase of infection (at one month post-infection), and that EmAgB2 is simultaneously down-regulated when compared to the other three genes. A negative relationship exists between the level of transcription and diversity of EmAgB genes. Moreover, no excess of non-synonymous substitutions was found among the distinct EmAgB alleles from a single host. Together, these results pointed to the effect of purifying selection, which seemed to eliminate the detrimental AgB variants generated during the development of the metacestode within the peritoneal cavity of its intermediate host.

Three essays on prices and selection

Competitive Market Segmentation Abstract In a two-firm model where each firm sells a high-quality and a low-quality version of a product, customers differ with respect to their brand preferences and their attitudes towards quality. We show that the standard result of quality-independent markups crucially depends on the assumption that the customers' valuation of quality is identical across firms. Once we relax this assumption, competition across qualities leads to second-degree price discrimination. We find that markups on low-quality products are higher if consuming a low-quality product involves a firm-specific disutility. Likewise, markups on high-quality products are higher if consuming a high-quality product creates a firm-specific surplus. Selection upon Wage Posting Abstract We discuss a model of a job market where firms announce salaries. Thereupon, they decide through the evaluation of a productivity test whether to hire applicants. Candidates for a job are locked in once they have applied at a given employer. Hence, such a market exhibits a specific form of the bargain-then-ripoff principle. With a single firm, the outcome is efficient. Under competition, what might be called "positive selection" leads to market failure. Thus our model provides a rationale for very small employment probabilities in some sectors. Exclusivity Clauses: Enhancing Competition, Raising Prices Abstract In a setting where retailers and suppliers compete for each other by offering binding contracts, exclusivity clauses serve as a competitive device. As a result of these clauses, firms addressed by contracts only accept the most favorable deal. Thus the contract-issuing parties have to squeeze their final customers and transfer the surplus within the vertical supply chain. We elaborate to what extent the resulting allocation depends on the sequence of play and discuss the implications of a ban on exclusivity clauses.

INCREASED EXPRESSION OF ONE OF TWO ADENOSINE DEAMINASE ALLELES IN A HUMAN CHORIOCARCINOMA CELL LINE FOLLOWING SELECTION WITH ADENINE NUCLEOSIDES

The human choriocarcinoma cell line JEG-3 is heterozygous at the adenosine deaminase (ADA) gene locus. Both allelic genes are under strong but incomplete repression causing a very low level expression of the gene locus. Because cytotoxic adenosine analogues such as 9-(beta)-D arabinofuranosyladenine (ara-A) and 9-(beta)-D xylofuranosyladenine (xyl-A) can be specifically detoxified by the action of ADA, these analogues were used to select for JEG-3 derived cells which had increased ADA expression. When JEG-3 cells were subjected to a multi-step, successively increasing dosage of either ara-A or xyl-A, resistant cells with increased ADA expression were generated. This increased ADA expression in the resistant cells was unstable, so that when the selective pressure was removed, cellular ADA expression would decrease. Subclone analysis of xyl-A resistant cells revealed that compared to parental JEG-3 cells, individual resistant cells had either elevated ADA levels or decreased adenosine kinase (ADK) levels or both. This altered ADA and ADK expression in the resistant cells were found to be independent events. Because of high endogenous tissue conversion factor (TCF) expression in the JEG-3 cells, the allelic nature of the increased ADA expression in most of the resistant cells could not be determined. However, several resistant subcloned cells were found to have lost TCF expression. These TCF('-) cells expressed only the ADA*2 allelic gene product. Cell fusion experiments demonstrated that the ADA*1 allelic gene was intact and functional in the A3-1A7 cell line. Chromosomal analysis of the A3-1A7 cells showed that they had no double-minutes or homogeneously staining chromosomal regions, although a pair of new chromosomes were found in these cells. Segregation analysis of the hybrid cells indicated that an ADA*2 allelic gene was probably located on this new chromosome. The analysis of the A3-1A7 cell line suggested that the expression of only ADA 2 in these cells was the result of possibly a cis-deregulation of the ADA gene locus or more probably an amplification of the ADA*2 allelic gene. Two effective positive selection systems for ADA('+) cells were also developed and tested. These selection systems should eventually lead to the isolation of the ADA gene.^

Cross-lineage expression of Ig-β (B29) in thymocytes: Positive and negative gene regulation to establish T cell identity

Developmental commitment involves activation of lineage-specific genes, stabilization of a lineage-specific gene expression program, and permanent inhibition of inappropriate characteristics. To determine how these processes are coordinated in early T cell development, the expression of T and B lineage-specific genes was assessed in staged subsets of immature thymocytes. T lineage characteristics are acquired sequentially, with germ-line T cell antigen receptor-β transcripts detected very early, followed by CD3ɛ and terminal deoxynucleotidyl transferase, then pTα, and finally RAG1. Only RAG1 expression coincides with commitment. Thus, much T lineage gene expression precedes commitment and does not depend on it. Early in the course of commitment to the T lineage, thymocytes lose the ability to develop into B cells. To understand how this occurs, we also examined expression of well defined B lineage-specific genes. Although λ5 and Ig-α are not expressed, the μ0 and Iμ transcripts from the unrearranged IgH locus are expressed early, in distinct patterns, then repressed just before RAG1 expression. By contrast, RNA encoding the B cell receptor component Ig-β was found to be transcribed in all immature thymocyte subpopulations and throughout most thymocyte differentiation. Ig-β expression is down-regulated only during positive selection of CD4+CD8– cells. Thus several key participants in the B cell developmental program are expressed in non-B lineage-committed cells, and one is maintained even through commitment to an alternative lineage, and repressed only after extensive T lineage differentiation. The results show that transcriptional activation of “lymphocyte-specific” genes can occur in uncommitted precursors, and that T lineage commitment is a composite of distinct positive and negative regulatory events.

Isolation and characterization of iron-independent positive dominant mutants of the diphtheria toxin repressor DtxR

It is well known that the functional activity of the diphtheria toxin repressor DtxR is controlled by iron, which serves as an essential cofactor necessary for activation of target DNA binding by this regulatory element. In this communication, we describe the isolation and characterization of a unique series of DtxR mutants that are constitutively active and repress the expression of β-galactosidase from a diphtheria tox promoter/operator–lacZ transcriptional fusion, even in the absence of iron. These self-activating mutants of DtxR (SAD) were isolated through the use of a positive selection system for the cloning of functional dtxR alleles and target DNA operator sites. Of the four independently isolated SAD mutants that were characterized, two (SAD2 and SAD11) were found to carry a single missense mutation (E175K) in their respective C-terminal SH3-like domains. In contrast, the mutant allele encoding SAD3 was found to carry a total of six missense mutations distributed throughout the N- and C-terminal domains of the repressor. Partial diploid analysis of strains carrying both native dtxR and alleles encoding either SAD2 or SAD3 demonstrate that these iron-independent mutants possess a positive dominant phenotype in the regulation of β-galactosidase expression from a diphtheria tox promoter/operator–lacZ transcriptional fusion.

Rapid evolution in plant chitinases: Molecular targets of selection in plant-pathogen coevolution

Many pathogen recognition genes, such as plant R-genes, undergo rapid adaptive evolution, providing evidence that these genes play a critical role in plant-pathogen coevolution. Surprisingly, whether rapid adaptive evolution also occurs in genes encoding other kinds of plant defense proteins is unknown. Unlike recognition proteins, plant chitinases attack pathogens directly, conferring disease resistance by degrading chitin, a component of fungal cell walls. Here, we show that nonsynonymous substitution rates in plant class I chitinase often exceed synonymous rates in the plant genus Arabis (Cruciferae) and in other dicots, indicating a succession of adaptively driven amino acid replacements. We identify individual residues that are likely subject to positive selection by using codon substitution models and determine the location of these residues on the three-dimensional structure of class I chitinase. In contrast to primate lysozymes and plant class III chitinases, structural and functional relatives of class I chitinase, the adaptive replacements of class I chitinase occur disproportionately in the active site cleft. This highly unusual pattern of replacements suggests that fungi directly defend against chitinolytic activity through enzymatic inhibition or other forms of chemical resistance and identifies target residues for manipulating chitinolytic activity. These data also provide empirical evidence that plant defense proteins not involved in pathogen recognition also evolve in a manner consistent with rapid coevolutionary interactions.

T-cell receptor (TCR) usage in Lewis rat experimental autoimmune encephalomyelitis: TCR beta-chain-variable-region V beta 8.2-positive T cells are not essential for induction and course of disease.

Predominant usage of V beta 8.2 gene segments, encoding a T-cell receptor (TCR) beta chain variable region, has been reported for pathogenic Lewis rat T cells reactive to myelin basic protein (MBP). However, up to 75% of the alpha/beta T cells in a panel of MBP-specific T-cell lines did not display TCR V beta 8.2, V beta 8.5, V beta 10, or V beta 16 elements. To further investigate TCR usage, we sorted the T-cell lines for V beta 8.2- and V beta 10-positive T cells or depleted the lines of cells with these TCRs. V beta 8.2-positive T cells and one of the depleted T-cell lines strongly reacted against the MBP peptide MBP-(68-88). The depleted T-cell line caused marked experimental autoimmune encephalomyelitis (EAE) even in Lewis rats in which endogenous V beta 8.2-positive T cells had been eliminated by neonatal treatment with anti-V beta 8.2 monoclonal antibodies. T-cell hybridomas generated from this line predominantly used V beta 3 TCR genes coexpressed with TCR V alpha 2 transcripts, which were also used by V beta 8.2-positive T cells. Furthermore, V beta 10-positive T cells reactive to MBP-(44-67) were encephalitogenic when injected immediately after positive selection. After induction of EAE by sorted V beta 8.2- or V beta 10-positive T-cell lines, immunocytochemical analysis of the spinal cord tissue showed a predominance of the injected TCR or of nontypable alpha/beta T cells after injection of the depleted line. Our results demonstrate heterogeneity of TCR beta-chain usage even for a single autoantigen in an inbred strain. Moreover, V beta 8.2-positive T cells are not essential for the induction and progression of adoptive-transfer EAE.

Selection pressure-driven evolution of the Epstein-Barr virus-encoded oncogene LMP1 in virus isolates from southeast Asia

The geographically constrained distribution of Epstein-Barr virus (EBV)-associated nasopharyngeal carcinoma (NPC) in southeast Asian populations suggests that both viral and host genetics may influence disease risk. Although susceptibility loci have been mapped within the human genome, the role of viral genetics in the focal distribution of NPC remains an enigma. Here we report a molecular phylogenetic analysis of an NPC-associated viral oncogene, LMP1, in a large panel of EBV isolates from southeast Asia and from Papua New Guinea, Africa, and Australia, regions of the world where NPC is and is not endemic, respectively. This analysis revealed that LMP1 sequences show a distinct geographic structure, indicating that the southeast Asian isolates have evolved as a lineage distinct from those of Papua New Guinea, African, and Australian isolates. Furthermore, a likelihood ratio test revealed that the C termini of the LMP1 sequences of the southeast Asian lineage are under significant positive selection pressure, particularly at some sites within the C-terminal activator regions. We also present evidence that although the N terminus and transmembrane region of LMP1 have undergone recombination, the C-terminal region of the gene has evolved without any history of recombination. Based on these observations, we speculate that selection pressure may be driving the LMP1 sequences in virus isolates from southeast Asia towards a more malignant phenotype, thereby influencing the endemic distribution of NPC in this region.

Forward-in-Time, Spatially Explicit Modeling Software to Simulate Genetic Lineages Under Selection

SELECTOR is a software package for studying the evolution of multiallelic genes under balancing or positive selection while simulating complex evolutionary scenarios that integrate demographic growth and migration in a spatially explicit population framework. Parameters can be varied both in space and time to account for geographical, environmental, and cultural heterogeneity. SELECTOR can be used within an approximate Bayesian computation estimation framework. We first describe the principles of SELECTOR and validate the algorithms by comparing its outputs for simple models with theoretical expectations. Then, we show how it can be used to investigate genetic differentiation of loci under balancing selection in interconnected demes with spatially heterogeneous gene flow. We identify situations in which balancing selection reduces genetic differentiation between population groups compared with neutrality and explain conflicting outcomes observed for human leukocyte antigen loci. These results and three previously published applications demonstrate that SELECTOR is efficient and robust for building insight into human settlement history and evolution.

Genetic Variation among Major Human Geographic Groups Supports a Peculiar Evolutionary Trend in PAX9

A total of 172 persons from nine South Amerindian, three African and one Eskimo populations were studied in relation to the Paired box gene 9 (PAX9) exon 3 (138 base pairs) as well as its 5' and 3' flanking intronic segments (232 bp and 220 bp, respectively) and integrated with the information available for the same genetic region from individuals of different geographical origins. Nine mutations were scored in exon 3 and six in its flanking regions; four of them are new South American tribe-specific singletons. Exon3 nucleotide diversity is several orders of magnitude higher than its intronic regions. Additionally, a set of variants in the PAX9 and 101 other genes related with dentition can define at least some dental morphological differences between Sub-Saharan Africans and non-Africans, probably associated with adaptations after the modern human exodus from Africa. Exon 3 of PAX9 could be a good molecular example of how evolvability works.