Analyses of six homologous proteins of Protochlamydia amoebophila UWE25 encoded by large GC-rich genes (lgr): a model of evolution and concatenation of leucine-rich repeats.

BACKGROUND: Along the chromosome of the obligate intracellular bacteria Protochlamydia amoebophila UWE25, we recently described a genomic island Pam100G. It contains a tra unit likely involved in conjugative DNA transfer and lgrE, a 5.6-kb gene similar to five others of P. amoebophila: lgrA to lgrD, lgrF. We describe here the structure, regulation and evolution of these proteins termed LGRs since encoded by "Large G+C-Rich" genes. RESULTS: No homologs to the whole protein sequence of LGRs were found in other organisms. Phylogenetic analyses suggest that serial duplications producing the six LGRs occurred relatively recently and nucleotide usage analyses show that lgrB, lgrE and lgrF were relocated on the chromosome. The C-terminal part of LGRs is homologous to Leucine-Rich Repeats domains (LRRs). Defined by a cumulative alignment score, the 5 to 18 concatenated octacosapeptidic (28-meric) LRRs of LGRs present all a predicted alpha-helix conformation. Their closest homologs are the 28-residue RI-like LRRs of mammalian NODs and the 24-meres of some Ralstonia and Legionella proteins. Interestingly, lgrE, which is present on Pam100G like the tra operon, exhibits Pfam domains related to DNA metabolism. CONCLUSION: Comparison of the LRRs, enable us to propose a parsimonious evolutionary scenario of these domains driven by adjacent concatenations of LRRs. Our model established on bacterial LRRs can be challenged in eucaryotic proteins carrying less conserved LRRs, such as NOD proteins and Toll-like receptors.

The murine interleukin-2 receptor gamma chain gene: organization, chromosomal localization and expression in the adult thymus.

Defects in the interleukin-2 receptor gamma (IL-2R gamma) chain in the man result in an X-linked severe combined immunodeficiency, SCIDX1, characterized by an absence of T-cell differentiation. This phenotype may result from pertubations in IL-2, IL-4-, IL-7- or IL-15-mediated signaling, as the IL-2R gamma chain forms an integral component of these receptor systems. We have isolated and characterized cDNA and genomic clones for the murine IL-2R gamma. The gene (Il2rg) is well conserved between mouse and man with respect to overall structure and size, and contains regions of high conservation in the promoter region as well. Il2rg maps to mouse X chromosome region 40, in a region of synteny with human Xq12-13.1. We have also explored the expression of the IL-2R gamma during thymocyte development. IL-2R gamma transcripts are detected in the earliest thymocyte precursor cells and persist throughout intrathymic development into the mature peripheral compartment. Genomic clones for the murine IL-2R gamma will allow for further studies on the regulation and function of this gene in vivo.

Copper acquisition by the SenC protein regulates aerobic respiration in Pseudomonas aeruginosa PAO1.

Aerobic respiration of Pseudomonas aeruginosa involves four terminal oxidases belonging to the heme-copper family (that is, three cytochrome c oxidases and one quinol oxidase) plus one copper-independent, cyanide-insensitive quinol oxidase (CIO). The PA0114 gene encoding an SCO1/SenC-type protein, which is known to be important for copper delivery to cytochrome c in yeast, Rhodobacter spp. and Agrobacterium tumefaciens, was found to be important for copper acquisition and aerobic respiration in P. aeruginosa. A PA0114 (senC) mutant grew poorly in low-copper media and had low cytochrome cbb(3)-type oxidase activity, but expressed CIO at increased levels, by comparison with the wild-type PAO1. Addition of copper reversed these phenotypes, suggesting that periplasmic copper capture by the SenC protein helps P. aeruginosa to adapt to copper deprivation.

Postglacial recolonization at a snail's pace (Trochulus villosus): confronting competing refugia hypotheses using model selection.

The localization of Last Glacial Maximum (LGM) refugia is crucial information to understand a species' history and predict its reaction to future climate changes. However, many phylogeographical studies often lack sampling designs intensive enough to precisely localize these refugia. The hairy land snail Trochulus villosus has a small range centred on Switzerland, which could be intensively covered by sampling 455 individuals from 52 populations. Based on mitochondrial DNA sequences (COI and 16S), we identified two divergent lineages with distinct geographical distributions. Bayesian skyline plots suggested that both lineages expanded at the end of the LGM. To find where the origin populations were located, we applied the principles of ancestral character reconstruction and identified a candidate refugium for each mtDNA lineage: the French Jura and Central Switzerland, both ice-free during the LGM. Additional refugia, however, could not be excluded, as suggested by the microsatellite analysis of a population subset. Modelling the LGM niche of T. villosus, we showed that suitable climatic conditions were expected in the inferred refugia, but potentially also in the nunataks of the alpine ice shield. In a model selection approach, we compared several alternative recolonization scenarios by estimating the Akaike information criterion for their respective maximum-likelihood migration rates. The 'two refugia' scenario received by far the best support given the distribution of genetic diversity in T. villosus populations. Provided that fine-scale sampling designs and various analytical approaches are combined, it is possible to refine our necessary understanding of species responses to environmental changes.

OpenFluDB, a database for human and animal influenza virus.

Although research on influenza lasted for more than 100 years, it is still one of the most prominent diseases causing half a million human deaths every year. With the recent observation of new highly pathogenic H5N1 and H7N7 strains, and the appearance of the influenza pandemic caused by the H1N1 swine-like lineage, a collaborative effort to share observations on the evolution of this virus in both animals and humans has been established. The OpenFlu database (OpenFluDB) is a part of this collaborative effort. It contains genomic and protein sequences, as well as epidemiological data from more than 27,000 isolates. The isolate annotations include virus type, host, geographical location and experimentally tested antiviral resistance. Putative enhanced pathogenicity as well as human adaptation propensity are computed from protein sequences. Each virus isolate can be associated with the laboratories that collected, sequenced and submitted it. Several analysis tools including multiple sequence alignment, phylogenetic analysis and sequence similarity maps enable rapid and efficient mining. The contents of OpenFluDB are supplied by direct user submission, as well as by a daily automatic procedure importing data from public repositories. Additionally, a simple mechanism facilitates the export of OpenFluDB records to GenBank. This resource has been successfully used to rapidly and widely distribute the sequences collected during the recent human swine flu outbreak and also as an exchange platform during the vaccine selection procedure. Database URL: http://openflu.vital-it.ch.

Control of jasmonate biosynthesis and senescence by miR319 targets.

Considerable progress has been made in identifying the targets of plant microRNAs, many of which regulate the stability or translation of mRNAs that encode transcription factors involved in development. In most cases, it is unknown, however, which immediate transcriptional targets mediate downstream effects of the microRNA-regulated transcription factors. We identified a new process controlled by the miR319-regulated clade of TCP (TEOSINTE BRANCHED/CYCLOIDEA/PCF) transcription factor genes. In contrast to other miRNA targets, several of which modulate hormone responses, TCPs control biosynthesis of the hormone jasmonic acid. Furthermore, we demonstrate a previously unrecognized effect of TCPs on leaf senescence, a process in which jasmonic acid has been proposed to be a critical regulator. We propose that miR319-controlled TCP transcription factors coordinate two sequential processes in leaf development: leaf growth, which they negatively regulate, and leaf senescence, which they positively regulate.

Evolutionary trajectories of primate genes involved in HIV pathogenesis.

The current availability of five complete genomes of different primate species allows the analysis of genetic divergence over the last 40 million years of evolution. We hypothesized that the interspecies differences observed in susceptibility to HIV-1 would be influenced by the long-range selective pressures on host genes associated with HIV-1 pathogenesis. We established a list of human genes (n = 140) proposed to be involved in HIV-1 biology and pathogenesis and a control set of 100 random genes. We retrieved the orthologous genes from the genome of humans and of four nonhuman primates (Pan troglodytes, Pongo pygmaeus abeli, Macaca mulatta, and Callithrix jacchus) and analyzed the nucleotide substitution patterns of this data set using codon-based maximum likelihood procedures. In addition, we evaluated whether the candidate genes have been targets of recent positive selection in humans by analyzing HapMap Phase 2 single-nucleotide polymorphisms genotyped in a region centered on each candidate gene. A total of 1,064 sequences were used for the analyses. Similar median K(A)/K(S) values were estimated for the set of genes involved in HIV-1 pathogenesis and for control genes, 0.19 and 0.15, respectively. However, genes of the innate immunity had median values of 0.37 (P value = 0.0001, compared with control genes), and genes of intrinsic cellular defense had K(A)/K(S) values around or greater than 1.0 (P value = 0.0002). Detailed assessment allowed the identification of residues under positive selection in 13 proteins: AKT1, APOBEC3G, APOBEC3H, CD4, DEFB1, GML, IL4, IL8RA, L-SIGN/CLEC4M, PTPRC/CD45, Tetherin/BST2, TLR7, and TRIM5alpha. A number of those residues are relevant for HIV-1 biology. The set of 140 genes involved in HIV-1 pathogenesis did not show a significant enrichment in signals of recent positive selection in humans (intraspecies selection). However, we identified within or near these genes 24 polymorphisms showing strong signatures of recent positive selection. Interestingly, the DEFB1 gene presented signatures of both interspecies positive selection in primates and intraspecies recent positive selection in humans. The systematic assessment of long-acting selective pressures on primate genomes is a useful tool to extend our understanding of genetic variation influencing contemporary susceptibility to HIV-1.

TACI, unlike BAFF-R, is solely activated by oligomeric BAFF and APRIL to support survival of activated B cells and plasmablasts.

The cytokine BAFF binds to the receptors TACI, BCMA, and BAFF-R on B cells, whereas APRIL binds to TACI and BCMA only. The signaling properties of soluble trimeric BAFF (BAFF 3-mer) were compared with those of higher-order BAFF oligomers. All forms of BAFF bound BAFF-R and TACI, and elicited BAFF-R-dependent signals in primary B cells. In contrast, signaling through TACI in mature B cells or plasmablasts was only achieved by higher-order BAFF and APRIL oligomers, all of which were also po-tent activators of a multimerization-dependent reporter signaling pathway. These results indicate that, although BAFF-R and TACI can provide B cells with similar signals, only BAFF-R, but not TACI, can respond to soluble BAFF 3-mer, which is the main form of BAFF found in circulation. BAFF 60-mer, an efficient TACI agonist, was also detected in plasma of BAFF transgenic and nontransgenic mice and was more than 100-fold more active than BAFF 3-mer for the activation of multimerization-dependent signals. TACI supported survival of activated B cells and plasmablasts in vitro, providing a rational basis to explain the immunoglobulin deficiency reported in TACI-deficient persons.

TCRep 3D: an automated in silico approach to study the structural properties of TCR repertoires.

TCRep 3D is an automated systematic approach for TCR-peptide-MHC class I structure prediction, based on homology and ab initio modeling. It has been considerably generalized from former studies to be applicable to large repertoires of TCR. First, the location of the complementary determining regions of the target sequences are automatically identified by a sequence alignment strategy against a database of TCR Vα and Vβ chains. A structure-based alignment ensures automated identification of CDR3 loops. The CDR are then modeled in the environment of the complex, in an ab initio approach based on a simulated annealing protocol. During this step, dihedral restraints are applied to drive the CDR1 and CDR2 loops towards their canonical conformations, described by Al-Lazikani et. al. We developed a new automated algorithm that determines additional restraints to iteratively converge towards TCR conformations making frequent hydrogen bonds with the pMHC. We demonstrated that our approach outperforms popular scoring methods (Anolea, Dope and Modeller) in predicting relevant CDR conformations. Finally, this modeling approach has been successfully applied to experimentally determined sequences of TCR that recognize the NY-ESO-1 cancer testis antigen. This analysis revealed a mechanism of selection of TCR through the presence of a single conserved amino acid in all CDR3β sequences. The important structural modifications predicted in silico and the associated dramatic loss of experimental binding affinity upon mutation of this amino acid show the good correspondence between the predicted structures and their biological activities. To our knowledge, this is the first systematic approach that was developed for large TCR repertoire structural modeling.

Species selectivity of mixed-lineage leukemia/trithorax and HCF proteolytic maturation pathways.

Site-specific proteolytic processing plays important roles in the regulation of cellular activities. The histone modification activity of the human trithorax group mixed-lineage leukemia (MLL) protein and the cell cycle regulatory activity of the cell proliferation factor herpes simplex virus host cell factor 1 (HCF-1) are stimulated by cleavage of precursors that generates stable heterodimeric complexes. MLL is processed by a protease called taspase 1, whereas the precise mechanisms of HCF-1 maturation are unclear, although they are known to depend on a series of sequence repeats called HCF-1(PRO) repeats. We demonstrate here that the Drosophila homologs of MLL and HCF-1, called Trithorax and dHCF, are both cleaved by Drosophila taspase 1. Although highly related, the human and Drosophila taspase 1 proteins display cognate species specificity. Thus, human taspase 1 preferentially cleaves MLL and Drosophila taspase 1 preferentially cleaves Trithorax, consistent with coevolution of taspase 1 and MLL/Trithorax proteins. HCF proteins display even greater species-specific divergence in processing: whereas dHCF is cleaved by the Drosophila taspase 1, human and mouse HCF-1 maturation is taspase 1 independent. Instead, human and Xenopus HCF-1PRO repeats are cleaved in vitro by a human proteolytic activity with novel properties. Thus, from insects to humans, HCF proteins have conserved proteolytic maturation but evolved different mechanisms.

Molecular basis of leukocyte rolling on PSGL-1. Predominant role of core-2 O-glycans and of tyrosine sulfate residue 51.

Interactions between the leukocyte adhesion receptor L-selectin and P-selectin glycoprotein ligand-1 play an important role in regulating the inflammatory response by mediating leukocyte tethering and rolling on adherent leukocytes. In this study, we have examined the effect of post-translational modifications of PSGL-1 including Tyr sulfation and presentation of sialylated and fucosylated O-glycans for L-selectin binding. The functional importance of these modifications was determined by analyzing soluble L-selectin binding and leukocyte rolling on CHO cells expressing various glycoforms of PSGL-1 or mutant PSGL-1 targeted at N-terminal Thr or Tyr residues. Simultaneous expression of core-2 beta1,6-N-acetylglucosaminyltransferase and fucosyltransferase VII was required for optimal L-selectin binding to PSGL-1. Substitution of Thr-57 by Ala but not of Thr-44, strongly decreased L-selectin binding and leukocyte rolling on PSGL-1. Substitution of Tyr by Phe revealed that PSGL-1 Tyr-51 plays a predominant role in mediating L-selectin binding and leukocyte rolling whereas Tyr-48 has a minor role, an observation that contrasts with the pattern seen for the interactions between PSGL-1 and P-selectin where Tyr-48 plays a key role. Molecular modeling analysis of L-selectin and P-selectin interactions with PSGL-1 further supported these observations. Additional experiments showed that core-2 O-glycans attached to Thr-57 were also of critical importance in regulating the velocity and stability of leukocyte rolling. These observations pinpoint the structural characteristics of PSGL-1 that are required for optimal interactions with L-selectin and may be responsible for the specific kinetic and mechanical bond properties of the L-selectin-PSGL-1 adhesion receptor-counterreceptor pair.

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.

O-GlcNAc transferase catalyzes site-specific proteolysis of HCF-1.

The human epigenetic cell-cycle regulator HCF-1 undergoes an unusual proteolytic maturation process resulting in stably associated HCF-1(N) and HCF-1(C) subunits that regulate different aspects of the cell cycle. Proteolysis occurs at six centrally located HCF-1(PRO)-repeat sequences and is important for activation of HCF-1(C)-subunit functions in M phase progression. We show here that the HCF-1(PRO) repeat is recognized by O-linked β-N-acetylglucosamine transferase (OGT), which both O-GlcNAcylates the HCF-1(N) subunit and directly cleaves the HCF-1(PRO) repeat. Replacement of the HCF-1(PRO) repeats by a heterologous proteolytic cleavage signal promotes HCF-1 proteolysis but fails to activate HCF-1(C)-subunit M phase functions. These results reveal an unexpected role of OGT in HCF-1 proteolytic maturation and an unforeseen nexus between OGT-directed O-GlcNAcylation and proteolytic maturation in HCF-1 cell-cycle regulation.

The expansion of amino-acid repeats is not associated to adaptive evolution in mammalian genes.

BACKGROUND: The expansion of amino acid repeats is determined by a high mutation rate and can be increased or limited by selection. It has been suggested that recent expansions could be associated with the potential of adaptation to new environments. In this work, we quantify the strength of this association, as well as the contribution of potential confounding factors. RESULTS: Mammalian positively selected genes have accumulated more recent amino acid repeats than other mammalian genes. However, we found little support for an accelerated evolutionary rate as the main driver for the expansion of amino acid repeats. The most significant predictors of amino acid repeats are gene function and GC content. There is no correlation with expression level. CONCLUSIONS: Our analyses show that amino acid repeat expansions are causally independent from protein adaptive evolution in mammalian genomes. Relaxed purifying selection or positive selection do not associate with more or more recent amino acid repeats. Their occurrence is slightly favoured by the sequence context but mainly determined by the molecular function of the gene.