SAMFIRE: multi-locus variant calling for time-resolved sequence data.

UNLABELLED: An increasingly common method for studying evolution is the collection of time-resolved short-read sequence data. Such datasets allow for the direct observation of rapid evolutionary processes, as might occur in natural microbial populations and in evolutionary experiments. In many circumstances, evolutionary pressure acting upon single variants can cause genomic changes at multiple nearby loci. SAMFIRE is an open-access software package for processing and analyzing sequence reads from time-resolved data, calling important single- and multi-locus variants over time, identifying alleles potentially affected by selection, calculating linkage disequilibrium statistics, performing haplotype reconstruction and exploiting time-resolved information to estimate the extent of uncertainty in reported genomic data. AVAILABILITY AND IMPLEMENTATION: C ++ code may be found at https://github.com/cjri/samfire/ CONTACT: chris.illingworth@gen.cam.ac.uk SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

In vitro selection of miltefosine resistance in promastigotes of Leishmania donovani from Nepal: genomic and metabolomic characterization.

In this study, we followed the genomic, lipidomic and metabolomic changes associated with the selection of miltefosine (MIL) resistance in two clinically derived Leishmania donovani strains with different inherent resistance to antimonial drugs (antimony sensitive strain Sb-S; and antimony resistant Sb-R). MIL-R was easily induced in both strains using the promastigote-stage, but a significant increase in MIL-R in the intracellular amastigote compared to the corresponding wild-type did not occur until promastigotes had adapted to 12.2 µM MIL. A variety of common and strain-specific genetic changes were discovered in MIL-adapted parasites, including deletions at the LdMT transporter gene, single-base mutations and changes in somy. The most obvious lipid changes in MIL-R promastigotes occurred to phosphatidylcholines and lysophosphatidylcholines and results indicate that the Kennedy pathway is involved in MIL resistance. The inherent Sb resistance of the parasite had an impact on the changes that occurred in MIL-R parasites, with more genetic changes occurring in Sb-R compared with Sb-S parasites. Initial interpretation of the changes identified in this study does not support synergies with Sb-R in the mechanisms of MIL resistance, though this requires an enhanced understanding of the parasite's biochemical pathways and how they are genetically regulated to be verified fully.

Computational studies of Family A and Family B GPCRs.

A full picture of the similarities between Family A and Family B GPCRs (G-protein coupled receptors) has been frustrated by the lack of clear homology between the respective sequences. Here, we review previous computational studies on GPCR dimerization in which the putative dimerization interfaces have been analysed using entropy, the ET (evolutionary trace) method and related methods. The results derived from multiple sequence alignments of Family A subfamilies have been mapped on to the rhodopsin crystal structure using standard alignments. Similarly, the results for the Family B alignments have been mapped on to the rhodopsin crystal structure using the 'cold-spot' alignment. For both Family A and Family B GPCRs, the sequence analysis indicates that there are functional sites on essentially all transmembrane helices, consistent with the parallel daisy chain model of GPCR oligomerization in which each GPCR makes interactions with a number of neighbouring GPCRs. The results are not too sensitive to the quality of the alignment. Molecular Dynamics simulations of the activation process within a single transmembrane bundle of the rhodopsin and the beta(2)-adrenergic receptor have been reviewed; the key observation, which is consistent with other computational studies, is that there is a translation and bending of helix 6, which contributes to a significant opening out of the intracellular face of the receptor, as shown in the accompanying movies. The simulations required the application of specific experiment-derived harmonic and half-harmonic distance restraints and so the application of such simulations to Family B GPCRs requires considerable care because of the alignment problem. Thus, in order to address the alignment problem, we have exploited the observation that GCR1, a plant GPCR, has homology with Family A, Family B and Family E GPCRs. The resulting alignment for transmembrane helix 3 is presented.