DnaSP, DNA polymorphism analyses by the coalescent and other methods

DnaSP is a software package for the analysis of DNA polymorphism data. Present version introduces several new modules and features which, among other options allow: (1) handling big data sets (~5 Mb per sequence); (2) conducting a large number of coalescent-based tests by Monte Carlo computer simulations; (3) extensive analyses of the genetic differentiation and gene flow among populations; (4) analysing the evolutionary pattern of preferred and unpreferred codons; (5) generating graphical outputs for an easy visualization of results. Availability: The software package, including complete documentation and examples, is freely available to academic users from: http://www.ub.es/dnasp

Odorant Receptor (Or) genes: polymorphism and divergence in the D. melanogaster and D. pseudoobscura Lineages

Background: In insects, like in most invertebrates, olfaction is the principal sensory modality, which provides animals with essential information for survival and reproduction. Odorant receptors are involved in this response, mediating interactions between an individual and its environment, as well as between individuals of the same or different species. The adaptive importance of odorant receptors renders them good candidates for having their variation shaped by natural selection. Methodology/Principal Findings: We analyzed nucleotide variation in a subset of eight Or genes located on the 3L chromosomal arm of Drosophila melanogaster in a derived population of this species and also in a population of Drosophila pseudoobscura. Some heterogeneity in the silent polymorphism to divergence ratio was detected in the D. melanogaster/D. simulans comparison, with a single gene (Or67b) contributing ~37% to the test statistic. However, no other signals of a very recent selective event were detected at this gene. In contrast, at the speciation timescale, the MK test uncovered the footprint of positive selection driving the evolution of two of the encoded proteins in both D. melanogaster ¿OR65c and OR67a ¿and D. pseudoobscura ¿OR65b1 and OR67c. Conclusions: The powerful polymorphism/divergence approach provided evidence for adaptive evolution at a rather high proportion of the Or genes studied after relatively recent speciation events. It did not provide, however, clear evidence for very recent selective events in either D. melanogaster or D. pseudoobscura.

Changes in chromosomal polymorphism and global warming: The case of Drosophila subobscura from Apatin (Serbia)

Abstract In this study, chromosomal inversion polymorphism data for a natural population of Drosophila subobscura from a swampy region near the town of Apatin (Serbia) were compared with data for the same population collected approximately 15 years earlier. The pattern of chromosomal inversion polymorphism changed over time. There were significant increases in the frequency of characteristic southern latitude ("warm" adapted) chromosomal arrangements and significant decreases in the frequency of characteristic northern latitude ("cold" adapted) chromosomal arrangements in the O and U chromosomes. The chromosomal arrangements O3+4 and O3+4+22 (derived from the O3+4 arrangement)showed significant increases in 2008 and 2009 with regard to the 1994 sample. There was also a significant increase (~50%) in the U1+2 arrangement, while U1+8+2 (a typical southern arrangement) was detected for the first time. Since the Apatin swampy population ofD. subobscura has existed for a long time in a stable habitat with high humidity that has not been changed by man our results indicate that natural selection has produced chromosomal changes in response to the increase in temperature that has occurred in the Balkan Peninsula of central southeastern European. Key words: chromosomal inversions, Drosophila subobscura, global warming, karyotypes.

Adaptive colour polymorphism of Acrida ungarica H. (Orthoptera: Acrididae) in a spatially heterogeneous environment

Intra-specific colour polymorphism provides a cryptic camouflage from predators in heterogeneous habitats. The orthoptera species, Acrida ungarica (Herbst, 1786) possess two well-distinguished colour morphs: brown and green and displays several disruptive colouration patterns within each morph to improve the crypsis. This study focused on how the features of the background environment relate to the proportion of the two morphs and to the intensity of disruptive colouration patterns in A. ungarica. As the two sexes are very distinct with respect to mass and length, we also distinctively tested the relationship for each sex. In accordance with the background matching hypothesis, we found that, for both sexes, the brown morph was in higher proportion at sites with a brown-dominant environment, and green morphs were in higher proportion in green-dominant environments. Globally, individuals in drier sites and in the drier year also had more intense disruptive colouration patterns, and brown morphs and females were also more striped. Colour patterns differed largely between populations and were significantly correlated with relevant environmental features. Even if A. ungarica is a polymorphic specialist, disruptive colouration still appears to provide strong benefits, particularly in some habitats. Moreover, because females are larger, they are less able to flee, which might explain the difference between sexes

Genome-wide DNA polymorphism analyses using VariScan

BACKGROUND: DNA sequence polymorphisms analysis can provide valuable information on the evolutionary forces shaping nucleotide variation, and provides an insight into the functional significance of genomic regions. The recent ongoing genome projects will radically improve our capabilities to detect specific genomic regions shaped by natural selection. Current available methods and software, however, are unsatisfactory for such genome-wide analysis. RESULTS: We have developed methods for the analysis of DNA sequence polymorphisms at the genome-wide scale. These methods, which have been tested on a coalescent-simulated and actual data files from mouse and human, have been implemented in the VariScan software package version 2.0. Additionally, we have also incorporated a graphical-user interface. The main features of this software are: i) exhaustive population-genetic analyses including those based on the coalescent theory; ii) analysis adapted to the shallow data generated by the high-throughput genome projects; iii) use of genome annotations to conduct a comprehensive analyses separately for different functional regions; iv) identification of relevant genomic regions by the sliding-window and wavelet-multiresolution approaches; v) visualization of the results integrated with current genome annotations in commonly available genome browsers. CONCLUSION: VariScan is a powerful and flexible suite of software for the analysis of DNA polymorphisms. The current version implements new algorithms, methods, and capabilities, providing an important tool for an exhaustive exploratory analysis of genome-wide DNA polymorphism data.

Linkage disequilibrium between a melanin-based colour polymorphism and tail length in the barn owl

A linkage disequilibrium between sexually selected and life history traits can be explained by three mutually non-exclusive mechanisms. Genes coding for two traits may be located close on the same chromosome, genes responsible for variation in one of the trait may pleiotropically alter the other, and non-random pairing with respect to two traits may generate a non-physical linkage disequilibrium between their genes. Knowledge of which of these three mechanisms is responsible for a covariation between two traits is of interest to understand why differently ornamented individuals differ in several phenotypic aspects. In Switzerland, barn owls Tyto alba mate randomly with respect to a colour polymorphism generating a large range of variants between reddish-brown and white, males being lighter coloured than females. Several studies have shown that plumage coloration is not neutral with respect to some life history components. To test whether coloration is genetically associated with body size, partial cross-fostering experiments were performed by exchanging some hatchlings between nests. These experiments showed that darker biological fathers produce longer-tailed offspring. This sex-specific pattern is consistent with the hypothesis of non-physical linkage disequilibrium. In line with this hypothesis, darker coloured males were mated with longer-tailed females, whereas female coloration was not associated with tail length of their mate. The finding that dark nestlings had a longer tail than their pale siblings also supports the physical linkage and pleiotropy hypotheses. Therefore, non-random pairing can generate or strengthen a genetic covariation between a secondary sexual character and a morphological trait.

An 18-kDa translocator protein (TSPO) polymorphism explains differences in binding affinity of the PET radioligand PBR28

[(11)C]PBR28 binds the 18-kDa Translocator Protein (TSPO) and is used in positron emission tomography (PET) to detect microglial activation. However, quantitative interpretations of signal are confounded by large interindividual variability in binding affinity, which displays a trimodal distribution compatible with a codominant genetic trait. Here, we tested directly for an underlying genetic mechanism to explain this. Binding affinity of PBR28 was measured in platelets isolated from 41 human subjects and tested for association with polymorphisms in TSPO and genes encoding other proteins in the TSPO complex. Complete agreement was observed between the TSPO Ala147Thr genotype and PBR28 binding affinity phenotype (P value=3.1 x 10(-13)). The TSPO Ala147Thr polymorphism predicts PBR28 binding affinity in human platelets. As all second-generation TSPO PET radioligands tested hitherto display a trimodal distribution in binding affinity analogous to PBR28, testing for this polymorphism may allow quantitative interpretation of TSPO PET studies with these radioligands.

Expression hierarchy of T cell epitopes from melanoma differentiation antigens: unexpected high level presentation of tyrosinase-HLA-A2 Complexes revealed by peptide-specific, MHC-restricted, TCR-like antibodies.

Peptide Ags presented by class I MHC molecules on human melanomas and that are recognized by CD8(+) T cells are the subjects of many studies of antitumor immunity and represent attractive candidates for therapeutic approaches. However, no direct quantitative measurements exist to reveal their expression hierarchy on the cell surface. Using novel recombinant Abs which bind these Ags with a peptide-specific, MHC-restricted manner, we demonstrate a defined pattern of expression hierarchy of peptide-HLA-A2 complexes derived from three major differentiation Ags: gp100, Melan-A/Mart-1, and tyrosinase. Studying melanoma cell lines derived from multiple patients, we reveal a surprisingly high level of presentation of tyrosinase-derived complexes and moderate to very low expression of complexes derived from other Ags. No correlation between Ag presentation and mRNA expression was found; however, protein stability may play a major role. These results provide new insights into the characteristics of Ag presentation and are particularly important when such targets are being considered for immunotherapy. These results may shed new light on relationships between Ag presentation and immune response to cancer Ags.

Interferon-γ Induces Expression of MHC Class II on Intestinal Epithelial Cells and Protects Mice from Colitis

Immune responses against intestinal microbiota contribute to the pathogenesis of inflammatory bowel diseases (IBD) and involve CD4(+) T cells, which are activated by major histocompatibility complex class II (MHCII) molecules on antigen-presenting cells (APCs). However, it is largely unexplored how inflammation-induced MHCII expression by intestinal epithelial cells (IEC) affects CD4(+) T cell-mediated immunity or tolerance induction in vivo. Here, we investigated how epithelial MHCII expression is induced and how a deficiency in inducible epithelial MHCII expression alters susceptibility to colitis and the outcome of colon-specific immune responses. Colitis was induced in mice that lacked inducible expression of MHCII molecules on all nonhematopoietic cells, or specifically on IECs, by continuous infection with Helicobacter hepaticus and administration of interleukin (IL)-10 receptor-blocking antibodies (anti-IL10R mAb). To assess the role of interferon (IFN)-γ in inducing epithelial MHCII expression, the T cell adoptive transfer model of colitis was used. Abrogation of MHCII expression by nonhematopoietic cells or IECs induces colitis associated with increased colonic frequencies of innate immune cells and expression of proinflammatory cytokines. CD4(+) T-helper type (Th)1 cells - but not group 3 innate lymphoid cells (ILCs) or Th17 cells - are elevated, resulting in an unfavourably altered ratio between CD4(+) T cells and forkhead box P3 (FoxP3)(+) regulatory T (Treg) cells. IFN-γ produced mainly by CD4(+) T cells is required to upregulate MHCII expression by IECs. These results suggest that, in addition to its proinflammatory roles, IFN-γ exerts a critical anti-inflammatory function in the intestine which protects against colitis by inducing MHCII expression on IECs. This may explain the failure of anti-IFN-γ treatment to induce remission in IBD patients, despite the association of elevated IFN-γ and IBD.

Fas (CD95/APO-1) limits the expansion of T lymphocytes in an environment of limited T-cell antigen receptor/MHC contacts.

Fas-deficient mice (Fas(lpr/lpr)) and humans have profoundly dysregulated T lymphocyte homeostasis, which manifests as an accumulation of CD4(+) and CD8(+) T cells as well as an unusual population of CD4(-)CD8(-)TCRαβ(+) T cells. To date, no unifying model has explained both the increased T-cell numbers and the origin of the CD4(-)CD8(-)TCRαβ(+) T cells. As Fas(lpr/lpr) mice raised in a germ-free environment still manifest lymphadenopathy, we considered that this process is primarily driven by recurrent low-avidity TCR signaling in response to self-peptide/MHC as occurs during homeostatic proliferation. In these studies, we developed two independent systems to decrease the number of self-peptide/MHC contacts. First, expression of MHC class I was reduced in OT-I TCR transgenic mice. Although OT-I Fas(lpr/lpr) mice did not develop lymphadenopathy characteristic of Fas(lpr/lpr) mice, in the absence of MHC class I, OT-I Fas(lpr/lpr) T cells accumulated as both CD8(+) and CD4(-)CD8(-) T cells. In the second system, re-expression of β(2)m limited to thymic cortical epithelial cells of Fas(lpr/lpr) β(2)m-deficient mice yielded a model in which polyclonal CD8(+) thymocytes entered a peripheral environment devoid of MHC class I. These mice accumulated significantly greater numbers of CD4(-)CD8(-)TCRαβ(+) T cells than conventional Fas(lpr/lpr) mice. Thus, Fas shapes the peripheral T-cell repertoire by regulating the survival of a subset of T cells proliferating in response to limited self-peptide/MHC contacts.

Functional analysis of monocyte MHC class II compartments.

Circulating monocytes, as dendritic cell and macrophage precursors, exhibit several functions usually associated with antigen-presenting cells, such as phagocytosis and presence of endosomal/lysosomal degradative compartments particularly enriched in Lamp-1, MHC class II molecules, and other proteins related to antigen processing and MHC class II loading [MHC class II compartments (MIICs)]. Ultrastructural analysis of these organelles indicates that, differently from the multivesicular bodies present in dendritic cells, in monocytes the MIICs are characterized by a single perimetral membrane surrounding an electron-dense core. Analysis of their content reveals enrichment in myeloperoxidase, an enzyme classically associated with azurophilic granules in granulocytes and mast cell secretory lysosomes. Elevation in intracellular free calcium levels in monocytes induced secretion of beta-hexosaminidase, cathepsins, and myeloperoxidase in the extracellular milieu; surface up-regulation of MHC class II molecules; and appearance of lysosomal resident proteins. The Ca(2+)-regulated surface transport mechanism of MHC class II molecules observed in monocytes is different from the tubulovesicular organization of the multivesicular bodies previously reported in dendritic cells and macrophages. Hence, in monocytes, MHC class II-enriched organelles combine degradative functions typical of lysosomes and regulated secretion typical of secretory lysosomes. More important, Ca(2+)-mediated up-regulation of surface MHC class II molecules is accompanied by extracellular release of lysosomal resident enzymes.

MHC class I-related chain A conjugated to antitumor antibodies can sensitize tumor cells to specific lysis by natural killer cells.

PURPOSE: As a first step for the development of a new cancer immunotherapy strategy, we evaluated whether antibody-mediated coating by MHC class I-related chain A (MICA) could sensitize tumor cells to lysis by natural killer (NK) cells. EXPERIMENTAL DESIGN: Recombinant MICA (rMICA) was chemically conjugated to Fab' fragments from monoclonal antibodies specific for tumor-associated antigens, such as carcinoembryonic antigen, HER2, or CD20. RESULTS: Flow cytometry analysis showed an efficient coating of MICA-negative human cancer cell lines with the Fab-rMICA conjugates. This was strictly dependent on the expression of the appropriate tumor-associated antigens in the target cells. Importantly, preincubation of the tumor cells with the appropriate Fab-rMICA conjugate resulted in NK cell-mediated tumor cell lysis. Antibody blocking of the NKG2D receptor in NK cells prevented conjugate-mediated tumor cell lysis. CONCLUSIONS: These results open the way to the development of immunotherapy strategies based on antibody-mediated targeting of MICA.

MHC class II expression is differentially regulated in plasmacytoid and conventional dendritic cells.

Major histocompatibility complex (MHC) class II-restricted antigen presentation is essential for the function of dendritic cells (DCs). We show here that plasmacytoid DCs (pDCs) differ from all other DC subsets with respect to expression of CIITA, the 'master regulator' of MHC class II genes. The gene encoding CIITA is controlled by three cell type-specific promoters: pI, pIII and pIV. With gene targeting in mice, we demonstrate that pDCs rely strictly on the B cell promoter pIII, whereas macrophages and all other DCs depend on pI. The molecular mechanisms driving MHC class II expression in pDCs are thus akin to those operating in lymphoid rather than myeloid cells.

aIr-1, a newly found locus on mouse chromosome 16 encoding a trans-acting activator factor for MHC class II gene expression.

RJ 2.2.5 is a human B cell line that has lost the capacity to express MHC class II genes. The human class II-positive phenotype is restored in somatic cell hybrids between RJ 2.2.5 and mouse spleen cells. By karyotype and molecular studies of an informative family of hybrids we have now shown that the reexpression of human class II gene products, as well as the maintenance of the mouse class II-positive phenotype, correlates with the presence of mouse chromosome 16. Thus, the existence on this mouse chromosome of a newly found locus, designated by us aIr-1, that determines a trans-acting activator function for class II gene expression, is established. Possible implications of this finding are discussed.

Interleukin-1 receptor antagonist gene (IL-1RN) polymorphism is a predictive factor of clinical pregnancy after IVF.

BACKGROUND: Only 25% of IVF transfer cycles lead to a clinical pregnancy, calling for continued technical progress but also more in depth analysis of patients' individual characteristics. The interleukin-1 (IL-1) system and matrix metalloproteinases (MMPs) are strongly implicated in embryo implantation. The genes coding for IL-1Ra (gene symbol IL-1RN), IL-1beta, MMP2 and MMP9 bear functional polymorphisms. We analysed the maternal genetic profile at these polymorphic sites in IVF patients, to determine possible correlations with IVF outcome. METHODS: One hundred and sixty women undergoing an IVF cycle were enrolled and a buccal smear was obtained. The presence of IL-1RN variable number of tandem repeats and IL-1B + 3953, MMP2-1306 and MMP9-1562 single nucleotide substitutions were determined. Patients were divided into pregnancy failures (119), biochemical pregnancies (8) and clinical pregnancies (33). RESULTS: There was a 40% decrease in IL-1RN*2 allele frequency (P = 0.024) and a 45% decrease in IL-1RN*2 carrier status in the clinical pregnancy group as compared to the pregnancy failure group (P = 0.017). This decrease was still statistically significant after a multivariate logistic regression analysis. The likelihood of a clinical pregnancy was decreased accordingly in IL-1RN*2 carriers: odds ratio = 0.349, 95% confidence interval = 0.2-0.8, P = 0.017. The IL-1B, MMP2 and MMP9 polymorphisms showed no correlation with IVF outcome. CONCLUSIONS: IL-1RN*2 allele carriage is associated with a poor prognosis of achieving a pregnancy after IVF.