Tree shape‐based approaches for the comparative study of cophylogeny

Cophylogeny is the congruence of phylogenetic relationships between two different groups of organisms due to their long‐term interaction. We investigated the use of tree shape distance measures to quantify the degree of cophylogeny. We implemented a reverse‐time simulation model of pathogen phylogenies within a fixed host tree, given cospeciation probability, host switching, and pathogen speciation rates. We used this model to evaluate 18 distance measures between host and pathogen trees including two kernel distances that we developed for labeled and unlabeled trees, which use branch lengths and accommodate different size trees. Finally, we used these measures to revisit published cophylogenetic studies, where authors described the observed associations as representing a high or low degree of cophylogeny. Our simulations demonstrated that some measures are more informative than others with respect to specific coevolution parameters especially when these did not assume extreme values. For real datasets, trees’ associations projection revealed clustering of high concordance studies suggesting that investigators are describing it in a consistent way. Our results support the hypothesis that measures can be useful for quantifying cophylogeny. This motivates their usage in the field of coevolution and supports the development of simulation‐based methods, i.e., approximate Bayesian computation, to estimate the underlying coevolutionary parameters.     

Deciphering the global phylogeography of a coastal shrub (Scaevola taccada) reveals the influence of multiple forces on contemporary population structure

The phylogeography of coastal plant species is heavily influenced by past sealevel fluctuations, dispersal barriers, and life-history traits, such as long-distance dispersal ability of the propagules. Unlike the widely studied mangroves, phylogeographic patterns have remained mostly obscure for other coastal plant species. In this study, we sampled 42 populations of Scaevola taccada (Gaertn.) Roxb., a coastal shrub of the family Goodeniaceae, from 17 countries across its distribution range. We used five chloroplast DNA (cpDNA) and 14 nuclear microsatellite (simple sequence repeat [SSR]) markers to assess the influence of abiotic factors and population genetic processes on the phylogeographic pattern of the species. Geographical distribution of cpDNA haplotypes suggests that the species originated in Australia, followed by historical dispersal and expansion of its geographic range. Multiple abiotic factors, including the sealevel changes during the Pleistocene, the presence of landmasses like the Malay Peninsula, and contemporary oceanic circulation patterns, restricted gene flow between geographically distinct populations, thereby creating low haplotype diversity and a strong population structure. Population genetic processes acted on these isolated populations, leading to high nuclear genetic diversity and population differentiation, as revealed from analyzing the polymorphic SSR loci. Although genetic divergence was mostly concordant between cpDNA and SSR data, asymmetrical gene flow and ancestral polymorphism could explain the discordance in the detailed genetic structure. Overall, our findings indicate that abiotic factors and population genetic processes interactively influenced the evolutionary history and current phylogeographic pattern of S. taccada across its distribution range.     

CpG/CpNpG motifs in the coding region are preferred sites for mutagenesis in the breast cancer susceptibility genes

The range of BRCA1/BRCA2 gene mutations is diverse and the mechanism accounting for this heterogeneity is obscure. To gain insight into the endogenous mutational mechanisms involved, we evaluated the association of specific sequences (i.e. CpG/CpNpG motifs, homonucleotides, short repeats) and mutations within the genes. We classified 1337 published mutations in BRCA1 (1765 BRCA2 mutations) for each specific sequence, and employed computer simulation combined with mathematical calculations to estimate the true underlying tendency of mutation occurrence. Interestingly, we found no mutational bias to homonucleotides and repeats in deletions/insertions and substitutions but striking bias to CpG/CpNpG in substitutions in both genes. This suggests that methylation-dependent DNA alterations would be a major mechanism for mutagenesis.     

Human TRBP and PACT directly interact with each other and associate with dicer to facilitate the production of small interfering RNA

Mammalian Dicer interacts with double-stranded RNA-binding protein TRBP or PACT to mediate RNA interference and micro-RNA processing. TRBP and PACT are structurally related but exert opposite regulatory activities on PKR. It is not understood whether TRBP and PACT are simultaneously required for Dicer. Here we show that TRBP directly interacts with PACT in vitro and in mammalian cells. TRBP and PACT form a triple complex with Dicer and facilitate the production of small interfering RNA (siRNA) by Dicer. Knockdown of both TRBP and PACT in cultured cells leads to significant inhibition of gene silencing mediated by short hairpin RNA but not by siRNA, suggesting that TRBP and PACT function primarily at the step of siRNA production. Taken together, these findings indicate that human TRBP and PACT directly interact with each other and associate with Dicer to stimulate the cleavage of double-stranded or short hairpin RNA to siRNA. Our work significantly alters the current model for the assembly and function of the Dicer-containing complex that generates siRNA and micro-RNA in human.     

Functional promiscuity correlates with conformational heterogeneity in A-class glutathione S-transferases

The structurally related glutathione S-transferase isoforms GSTA1-1 and GSTA4-4 differ greatly in their relative catalytic promiscuity. GSTA1-1 is a highly promiscuous detoxification enzyme. In contrast, GSTA4-4 exhibits selectivity for congeners of the lipid peroxidation product 4-hydroxynonenal. The contribution of protein dynamics to promiscuity has not been studied. Therefore, hydrogen/deuterium exchange mass spectrometry (H/DX) and fluorescence lifetime distribution analysis were performed with glutathione S-transferases A1-1 and A4-4. Differences in local dynamics of the C-terminal helix were evident as expected on the basis of previous studies. However, H/DX demonstrated significantly greater solvent accessibility throughout most of the GSTA1-1 sequence compared with GSTA4-4. A Phe-111/Tyr-217 aromatic-aromatic interaction in A4-4, which is not present in A1-1, was hypothesized to increase core packing. "Swap" mutants that eliminate this interaction from A4-4 or incorporate it into A1-1 yield H/DX behavior that is intermediate between the wild type templates. In addition, the single Trp-21 residue of each isoform was exploited to probe the conformational heterogeneity at the intrasubunit domain-domain interface. Excited state fluorescence lifetime distribution analysis indicates that this core residue is more conformationally heterogeneous in GSTA1-1 than in GSTA4-4, and this correlates with greater stability toward urea denaturation for GSTA4-4. The fluorescence distribution and urea sensitivity of the mutant proteins were intermediate between the wild type templates. The results suggest that the differences in protein dynamics of these homologs are global. The results suggest also the possible importance of extensive conformational plasticity to achieve high levels of functional promiscuity, possibly at the cost of stability.     

Methods of using click chemistry in the discovery of enzyme inhibitors

This protocol describes the step-by-step procedures for the efficient assembly of bidentate inhibitor libraries of a target enzyme, using the so-called 'click chemistry' between an alkyne-bearing core group and an azide-modified peripheral group, followed by direct biological screening for the identification of potential 'hits'. The reaction is highlighted by its modularity, high efficiency (approximately 100% yield in most cases) and tolerance toward many functional groups present in the fragments, as well as biocompatibility (typically carried out in aqueous conditions with small amounts of biocompatible catalysts). The approach consists of three steps: (i) chemical synthesis of alkyne-bearing protein tyrosine phosphatase or matrix metalloprotease core groups and diverse azide-modified peripheral groups; (ii) click chemistry to assemble the bidentate inhibitor libraries; and (iii) direct screening of the libraries with target enzymes using 384-well microplate assays. Following the chemical synthesis of the core and peripheral groups and optimization of the click chemistry conditions (approximately 1 week), steps (ii) and (iii) take 3 d to complete (approximately 1-2 d for library assembly and 1 d for inhibitor screening).     

Synthesis of cationic single-isomer cyclodextrins for the chiral separation of amino acids and anionic pharmaceuticals

We describe a protocol for the synthesis of mono-6(A)-(1-butyl-3-imidazolium)-6(A)-deoxy-beta-cyclodextrin chloride (BIMCD), a cationic, water-soluble cyclodextrin used in the chiral separation of amino acids and anionic pharmaceuticals by capillary electrophoresis. Starting from commercially available chemicals, BIMCD is synthesized in five steps. The first step involves a nucleophilic substitution between p-toluenesulfonyl chloride and imidazole to afford 1-(p-toluenesulfonyl)imidazole (A). In the second step, a nucleophilic substitution between beta-cyclodextrin and A affords mono-6(A)-(p-toluenesulfonyl)-6(A)-deoxy-beta-cyclodextrin (B). In the third step, a nucleophilic substitution between 1-bromobutane and imidazole affords 1-butylimidazole (C). In the fourth step, a nucleophilic addition between A and C affords BIMCD tosylate. In the final step, anion exchange using an ion-exchange resin yields BIMCD as a highly water-soluble solid. Each step takes up to 2 d, including the time required for product purification. The overall protocol requires approximately 6 d.     

Power to detect risk alleles using genome-wide tag SNP panels

Advances in high-throughput genotyping and the International HapMap Project have enabled association studies at the whole-genome level. We have constructed whole-genome genotyping panels of over 550,000 (HumanHap550) and 650,000 (HumanHap650Y) SNP loci by choosing tag SNPs from all populations genotyped by the International HapMap Project. These panels also contain additional SNP content in regions that have historically been overrepresented in diseases, such as nonsynonymous sites, the MHC region, copy number variant regions and mitochondrial DNA. We estimate that the tag SNP loci in these panels cover the majority of all common variation in the genome as measured by coverage of both all common HapMap SNPs and an independent set of SNPs derived from complete resequencing of genes obtained from SeattleSNPs. We also estimate that, given a sample size of 1,000 cases and 1,000 controls, these panels have the power to detect single disease loci of moderate risk (λ ~ 1.8–2.0). Relative risks as low as λ ~ 1.1–1.3 can be detected using 10,000 cases and 10,000 controls depending on the sample population and disease model. If multiple loci are involved, the power increases significantly to detect at least one locus such that relative risks 20%–35% lower can be detected with 80% power if between two and four independent loci are involved. Although our SNP selection was based on HapMap data, which is a subset of all common SNPs, these panels effectively capture the majority of all common variation and provide high power to detect risk alleles that are not represented in the HapMap data.