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.

A single point mutation in the pore region of the epithelial Na+ channel changes ion selectivity by modifying molecular sieving.

The epithelial Na+ channel (ENaC) belongs to a new class of channel proteins called the ENaC/DEG superfamily involved in epithelial Na+ transport, mechanotransduction, and neurotransmission. The role of ENaC in Na+ homeostasis and in the control of blood pressure has been demonstrated recently by the identification of mutations in ENaC beta and gamma subunits causing hypertension. The function of ENaC in Na+ reabsorption depends critically on its ability to discriminate between Na+ and other ions like K+ or Ca2+. ENaC is virtually impermeant to K+ ions, and the molecular basis for its high ionic selectivity is largely unknown. We have identified a conserved Ser residue in the second transmembrane domain of the ENaC alpha subunit (alphaS589), which when mutated allows larger ions such as K+, Rb+, Cs+, and divalent cations to pass through the channel. The relative ion permeability of each of the alphaS589 mutants is related inversely to the ionic radius of the permeant ion, indicating that alphaS589 mutations increase the molecular cutoff of the channel by modifying the pore geometry at the selectivity filter. Proper geometry of the pore is required to tightly accommodate Na+ and Li+ ions and to exclude larger cations. We provide evidence that ENaC discriminates between cations mainly on the basis of their size and the energy of dehydration.

Molecular cophylogenetic relationships between European bats and their ectoparasitic mites (Acari, Spinturnicidae).

Cospeciation between host-parasite species is generally thought to result in mirror-image congruent phylogenies. Incongruence can be explained by mechanisms such as host switching, duplication, failure to speciate and sorting events. To investigate the level of association in the host-parasite relationship between Spinturnicid mites and their bat hosts, we constructed the phylogenetic tree of the genus Spinturnix (Acari, Mesostigmata) and compared it to the host phylogeny. We sequenced 938bp of the mitochondrial 16S rDNA and Cytochrome Oxydase subunit I (COI) genes among eleven morphospecies of Spinturnix collected on 20 European Vespertilionid and Rhinolophid bat species. Phylogenetic reconstruction of hosts and parasites showed statistical evidence for cospeciation and suggested that their evolutionary history involved also failure to speciate events and host switches. The latter seem to be mainly promoted by similar roosting habits of the host. As currently understood, host associations of Spinturnicid mites likely results from a complex interaction between the phylogenetic history of the host and the behaviour and the ecology of both parasite and host.

Genetic isolation of insular populations of the Maghrebian bat, Myotis punicus, in the Mediterranean Basin

We investigate the population genetic structure of the Maghrebian bat, Myotis punicus, between the mainland and islands to assess the island colonization pattern and current gene flow between nearby islands and within the mainland. Location North Africa and the Mediterranean islands of Corsica and Sardinia. Methods We sequenced part of the control region (HVII) of 79 bats across 11 colonies. The phylogeographical pattern was assessed by analysing molecular diversity indices, examining differentiation among populations and estimating divergence time. In addition, we genotyped 182 bats across 10 colonies at seven microsatellite loci. We used analysis of molecular variance and a Bayesian approach to infer nuclear population structure. Finally, we estimated sex-specific dispersal between Corsica and Sardinia. Results Mitochondrial analyses indicated that colonies between Corsica, Sardinia and North Africa are highly differentiated. Within islands there was no difference between colonies, while at the continental level Moroccan and Tunisian populations were highly differentiated. Analyses with seven microsatellite loci showed a similar pattern. The sole difference was the lack of nuclear differentiation between populations in North Africa, suggesting a male-biased dispersal over the continental area. The divergence time of Sardinian and Corsican populations was estimated to date back to the early and mid-Pleistocene. Main conclusions Island colonization by the Maghrebian bats seems to have occurred in a stepping-stone manner and certainly pre-dated human colonization. Currently, open water seems to prevent exchange of bats between the two islands, despite their ability to fly and the narrowness of the strait of Bonifacio. Corsican and Sardinian populations are thus currently isolated from any continental gene pool and must therefore be considered as different evolutionarily significant units (ESU).

RIP4 (DIK/PKK), a novel member of the RIP kinase family, activates NF-kappa B and is processed during apoptosis.

RIP1 and its homologs, RIP2 and RIP3, form part of a family of Ser/Thr kinases that regulate signal transduction processes leading to NF-kappa B activation. Here, we identify RIP4 (DIK/PKK) as a novel member of the RIP kinase family. RIP4 contains an N-terminal RIP-like kinase domain and a C-terminal region characterized by the presence of 11 ankyrin repeats. Overexpression of RIP4 leads to activation of NF-kappa B and JNK. Kinase inactive RIP4 or a truncated version containing the ankyrin repeats have a dominant negative (DN) effect on NF-kappa B induction by multiple stimuli. RIP4 binds to several members of the TRAF protein family, and DN versions of TRAF1, TRAF3 and TRAF6 inhibit RIP4-induced NF-kappa B activation. Moreover, RIP4 is cleaved after Asp340 and Asp378 during Fas-induced apoptosis. These data suggest that RIP4 is involved in NF-kappa B and JNK signaling and that caspase-dependent processing of RIP4 may negatively regulate NF-kappa B-dependent pro-survival or pro-inflammatory signals.

CD3 delta establishes a functional link between the T cell receptor and CD8.

T cells expressing T cell receptor (TCR) complexes that lack CD3 delta, either due to deletion of the CD3 delta gene, or by replacement of the connecting peptide of the TCR alpha chain, exhibit severely impaired positive selection and TCR-mediated activation of CD8 single-positive T cells. Because the same defects have been observed in mice expressing no CD8 beta or tailless CD8 beta, we examined whether CD3 delta serves to couple TCR.CD3 with CD8. To this end we used T cell hybridomas and transgenic mice expressing the T1 TCR, which recognizes a photoreactive derivative of the PbCS 252-260 peptide in the context of H-2K(d). We report that, in thymocytes and hybridomas expressing the T1 TCR.CD3 complex, CD8 alpha beta associates with the TCR. This association was not observed on T1 hybridomas expressing only CD8 alpha alpha or a CD3 delta(-) variant of the T1 TCR. CD3 delta was selectively co-immunoprecipitated with anti-CD8 antibodies, indicating an avid association of CD8 with CD3 delta. Because CD8 alpha beta is a raft constituent, due to this association a fraction of TCR.CD3 is raft-associated. Cross-linking of these TCR-CD8 adducts results in extensive TCR aggregate formation and intracellular calcium mobilization. Thus, CD3 delta couples TCR.CD3 with raft-associated CD8, which is required for effective activation and positive selection of CD8(+) T cells.

Molecular identification and characterization of the Arabidopsis delta(3,5),delta(2,4)-dienoyl-coenzyme A isomerase, a peroxisomal enzyme participating in the beta-oxidation cycle of unsaturated fatty acids.

Degradation of unsaturated fatty acids through the peroxisomal beta-oxidation pathway requires the participation of auxiliary enzymes in addition to the enzymes of the core beta-oxidation cycle. The auxiliary enzyme delta(3,5),delta(2,4)-dienoyl-coenzyme A (CoA) isomerase has been well studied in yeast (Saccharomyces cerevisiae) and mammals, but no plant homolog had been identified and characterized at the biochemical or molecular level. A candidate gene (At5g43280) was identified in Arabidopsis (Arabidopsis thaliana) encoding a protein showing homology to the rat (Rattus norvegicus) delta(3,5),delta(2,4)-dienoyl-CoA isomerase, and possessing an enoyl-CoA hydratase/isomerase fingerprint as well as aspartic and glutamic residues shown to be important for catalytic activity of the mammalian enzyme. The protein, named AtDCI1, contains a peroxisome targeting sequence at the C terminus, and fusion of a fluorescent protein to AtDCI1 directed the chimeric protein to the peroxisome in onion (Allium cepa) cells. AtDCI1 expressed in Escherichia coli was shown to have delta(3,5),delta(2,4)-dienoyl-CoA isomerase activity in vitro. Furthermore, using the synthesis of polyhydroxyalkanoate in yeast peroxisomes as an analytical tool to study the beta-oxidation cycle, expression of AtDCI1 was shown to complement the yeast mutant deficient in the delta(3,5),delta(2,4)-dienoyl-CoA isomerase, thus showing that AtDCI1 is also appropriately targeted to the peroxisome in yeast and has delta(3,5),delta(2,4)-dienoyl-CoA isomerase activity in vivo. The AtDCI1 gene is expressed constitutively in several tissues, but expression is particularly induced during seed germination. Proteins showing high homology with AtDCI1 are found in gymnosperms as well as angiosperms belonging to the Monocotyledon or Dicotyledon classes.

PchR-box recognition by the AraC-type regulator PchR of Pseudomonas aeruginosa requires the siderophore pyochelin as an effector.

Under iron limitation, the opportunistic human pathogen Pseudomonas aeruginosa produces the siderophore pyochelin. When secreted into the extracellular environment, pyochelin complexes ferric ions and delivers them, via the outer membrane receptor FptA, to the bacterial cytoplasm. Extracellular pyochelin also acts as a signalling molecule, inducing the expression of pyochelin biosynthesis and uptake genes by a mechanism involving the AraC-type regulator PchR. We have identified a 32 bp conserved sequence element (PchR-box) in promoter regions of pyochelin-controlled genes and we show that the PchR-box in the pchR-pchDCBA intergenic region is essential for the induction of the pyochelin biosynthetic operon pchDCBA and the repression of the divergently transcribed pchR gene. PchR was purified as a fusion with maltose-binding protein (MBP). Mobility shift assays demonstrated specific binding of MBP-PchR to the PchR-box in the presence, but not in the absence of pyochelin and iron. PchR-box mutations that interfered with pyochelin-dependent regulation in vivo, also affected pyochelin-dependent PchR-box recognition in vitro. We conclude that pyochelin, probably in its iron-loaded state, is the intracellular effector required for PchR-mediated regulation. The fact that extracellular pyochelin triggers this regulation suggests that the siderophore can enter the cytoplasm.

Origin and evolution of homologous repeated sequences in the mitochondrial DNA control region of shrews.

The complete mitochondrial DNA (mtDNA) control region was amplified and directly sequenced in two species of shrew, Crocidura russula and Sorex araneus (Insectivora, Mammalia). The general organization is similar to that found in other mammals: a central conserved region surrounded by two more variable domains. However, we have found in shrews the simultaneous presence of arrays of tandem repeats in potential locations where repeats tend to occur separately in other mammalian species. These locations correspond to regions which are associated with a possible interruption of the replication processes, either at the end of the three-stranded D-loop structure or toward the end of the heavy-strand replication. In the left domain the repeated sequences (R1 repeats) are 78 bp long, whereas in the right domain the repeats are 12 bp long in C. russula and 14 bp long in S. araneus (R2 repeats). Variation in the copy number of these repeated sequences results in mtDNA control region length differences. Southern blot analysis indicates that level of heteroplasmy (more than one mtDNA form within an individual) differs between species. A comparative study of the R2 repeats in 12 additional species representing three shrew subfamilies provides useful indications for the understanding of the origin and the evolution of these homologous tandemly repeated sequences. An asymmetry in the distribution of variants within the arrays, as well as the constant occurrence of shorter repeated sequences flanking only one side of the R2 arrays, could be related to asymmetry in the replication of each strand of the mtDNA molecule. The pattern of sequence and length variation within and between species, together with the capability of the arrays to form stable secondary structures, suggests that the dominant mechanism involved in the evolution of these arrays in unidirectional replication slippage.

Swiss EMBnet node web server.

EMBnet is a consortium of collaborating bioinformatics groups located mainly within Europe (http://www.embnet.org). Each member country is represented by a 'node', a group responsible for the maintenance of local services for their users (e.g. education, training, software, database distribution, technical support, helpdesk). Among these services a web portal with links and access to locally developed and maintained software is essential and different for each node. Our web portal targets biomedical scientists in Switzerland and elsewhere, offering them access to a collection of important sequence analysis tools mirrored from other sites or developed locally. We describe here the Swiss EMBnet node web site (http://www.ch.embnet.org), which presents a number of original services not available anywhere else.

In vitro emergence of rifampicin resistance in Propionibacterium acnes and molecular characterization of mutations in the rpoB gene.

OBJECTIVES: Activity of rifampicin against Propionibacterium acnes biofilms was recently demonstrated, but rifampicin resistance has not yet been described in this organism. We investigated the in vitro emergence of rifampicin resistance in P. acnes and characterized its molecular background. METHODS: P. acnes ATCC 11827 was used (MIC 0.007 mg/L). The mutation rate was determined by inoculation of 10(9) cfu of P. acnes on rifampicin-containing agar plates incubated anaerobically for 7 days. Progressive emergence of resistance was studied by serial exposure to increasing concentrations of rifampicin in 72 h cycles using a low (10(6) cfu/mL) and high (10(8) cfu/mL) inoculum. The stability of resistance was determined after three subcultures of rifampicin-resistant isolates on rifampicin-free agar. For resistant mutants, the whole rpoB gene was amplified, sequenced and compared with a P. acnes reference sequence (NC006085). RESULTS: P. acnes growth was observed on rifampicin-containing plates with mutation rates of 2 ± 1 cfu  ×  10(-9) (4096× MIC) and 12 ± 5 cfu  ×  10(-9) (4 × MIC). High-level rifampicin resistance emerged progressively after 4 (high inoculum) and 13 (low inoculum) cycles. In rifampicin-resistant isolates, the MIC remained >32 mg/L after three subcultures. Mutations were detected in clusters I (amino acids 418-444) and II (amino acids 471-486) of the rpoB gene after sequence alignment with a Staphylococcus aureus reference sequence (CAA45512). The five following substitutions were found: His-437 → Tyr, Ser-442 → Leu, Leu-444 → Ser, Ile-483 → Val and Ser-485 → Leu. CONCLUSION: The rifampicin MIC increased from highly susceptible to highly resistant values. The resistance remained stable and was associated with mutations in the rpoB gene. To our knowledge, this is the first report of the emergence of rifampicin resistance in P. acnes.

Unexpected Role for Connexin37 in Circulating FVIII and Coagulation

The gap junction protein connexin37 (Cx37) plays an important role in cell-cell communication in the vasculature. Cx37 is expressed in endothelial cells, platelets and megakaryocytes. We have recently shown that Cx37 limits thrombus propensity by permitting intercellular signaling between aggregating platelets. Here, we have performed high throughput phage display to identify potential binding partners for the regulatory intracellular C-terminus of Cx37 (Cx37CT). We retrieved 2 consensus binding motifs for Cx37CT: WHK...[K,R]XP... and FH-K...[K,R]XXP.... Sequence alignment against the NCBI protein database indicated 66% homology of one the selected peptides with FVIII B-domain. We performed cross-linking reactions using BS3 and confirmed that an 11-mer peptide of the FVIII B-domain sequence linked to recombinant Cx37CT. In vitro binding of this peptide to Cx37CT was also confirmed by surface plasmon resonance. The dissociation constant of FVIII B-domain peptides to Cx37CT was ~20 uM. Other peptide sequences, designed upstream or downstream of the FVIII B-domain sequence, showed very low or no affinity for Cx37CT. Finally, in vivo studies revealed that thrombin generation in platelet-poor plasma from Cx37-/- mice (endogenous thrombin potential: 634±11 nM min, mean±SEM) was increased compared to Cx37+/+ mice (427±12, P<0.001). Moreover, partial activated thromboplastin time (aPTT) was shorter in Cx37-/- (39.7±1.5 s) than in Cx37+/+ mice (45.9±1.8, P=0.03), whereas prothrombin time was comparable. The shorter aPTT in Cx37-/- mice correlated with higher circulating FVIII activity (46.0±0.7 vs. 53.5±2.7 s for Cx37+/+, P=0.03). Overall, our data show for the first time a functional interaction between FVIII and Cx37. This interaction may be relevant for the control of FVIII secretion and, thereby, in the regulation of levels of FVIII circulating in blood. In addition, these results may open new perspectives to improve the efficiency of recombinant FVIII manufacturing.

Numerous potentially functional but non-genic conserved sequences on human chromosome 21.

The use of comparative genomics to infer genome function relies on the understanding of how different components of the genome change over evolutionary time. The aim of such comparative analysis is to identify conserved, functionally transcribed sequences such as protein-coding genes and non-coding RNA genes, and other functional sequences such as regulatory regions, as well as other genomic features. Here, we have compared the entire human chromosome 21 with syntenic regions of the mouse genome, and have identified a large number of conserved blocks of unknown function. Although previous studies have made similar observations, it is unknown whether these conserved sequences are genes or not. Here we present an extensive experimental and computational analysis of human chromosome 21 in an effort to assign function to sequences conserved between human chromosome 21 (ref. 8) and the syntenic mouse regions. Our data support the presence of a large number of potentially functional non-genic sequences, probably regulatory and structural. The integration of the properties of the conserved components of human chromosome 21 to the rapidly accumulating functional data for this chromosome will improve considerably our understanding of the role of sequence conservation in mammalian genomes.