Evaluation of a 7-gene genetic profile for athletic endurance phenotype in Ironman championship triathletes

Polygenic profiling has been proposed for elite endurance performance, using an additive model determining the proportion of optimal alleles in endurance athletes. To investigate this model’s utility for elite triathletes, we genotyped seven polymorphisms previously associated with an endurance polygenic profile (ACE Ins/Del, ACTN3 Arg577Ter, AMPD1 Gln12Ter, CKMM 1170bp/985+185bp, HFE His63Asp, GDF8 Lys153Arg and PPARGC1A Gly482Ser) in a cohort of 196 elite athletes who participated in the 2008 Kona Ironman championship triathlon. Mean performance time (PT) was not significantly different in individual marker analysis. Age, sex, and continent of origin had a significant influence on PT and were adjusted for. Only the AMPD1 endurance-optimal Gln allele was found to be significantly associated with an improvement in PT (model p=5.79 x 10-17, AMPD1 genotype p=0.01). Individual genotypes were combined into a total genotype score (TGS); TGS distribution ranged from 28.6 to 92.9, concordant with prior studies in endurance athletes (mean±SD: 60.75±12.95). TGS distribution was shifted toward higher TGS in the top 10% of athletes, though the mean TGS was not significantly different (p=0.164) and not significantly associated with PT even when adjusted for age, sex, and origin. Receiver operating characteristic curve analysis determined that TGS alone could not significantly predict athlete finishing time with discriminating sensitivity and specificity for three outcomes (less than median PT, less than mean PT, or in the top 10%), though models with the age, sex, continent of origin, and either TGS or AMPD1 genotype could. These results suggest three things: that more sophisticated genetic models may be necessary to accurately predict athlete finishing time in endurance events; that non-genetic factors such as training are hugely influential and should be included in genetic analyses to prevent confounding; and that large collaborations may be necessary to obtain sufficient sample sizes for powerful and complex analyses of endurance performance.

Ligand Specificity of Group I Biotin Protein Ligase of Mycobacterium tuberculosis

Background: Fatty acids are indispensable constituents of mycolic acids that impart toughness & permeability barrier to the cell envelope of M. tuberculosis. Biotin is an essential co-factor for acetyl-CoA carboxylase (ACC) the enzyme involved in the synthesis of malonyl-CoA, a committed precursor, needed for fatty acid synthesis. Biotin carboxyl carrier protein (BCCP) provides the co-factor for catalytic activity of ACC. Methodology/Principal Findings: BPL/BirA (Biotin Protein Ligase), and its substrate, biotin carboxyl carrier protein (BCCP) of Mycobacterium tuberculosis (Mt) were cloned and expressed in E. coli BL21. In contrast to EcBirA and PhBPL, the similar to 29.5 kDa MtBPL exists as a monomer in native, biotin and bio-5'AMP liganded forms. This was confirmed by molecular weigt profiling by gel filtration on Superdex S-200 and Dynamic Light Scattering (DLS). Computational docking of biotin and bio-5'AMP to MtBPL show that adenylation alters the contact residues for biotin. MtBPL forms 11 H-bonds with biotin, relative to 35 with bio-5'AMP. Docking simulations also suggest that bio-5'AMP hydrogen bonds to the conserved `GRGRRG' sequence but not biotin. The enzyme catalyzed transfer of biotin to BCCP was confirmed by incorporation of radioactive biotin and by Avidin blot. The K-m for BCCP was similar to 5.2 mu M and similar to 420 nM for biotin. MtBPL has low affinity (K-b = 1.06 x 10(-6) M) for biotin relative to EcBirA but their K-m are almost comparable suggesting that while the major function of MtBPL is biotinylation of BCCP, tight binding of biotin/bio-5'AMP by EcBirA is channeled for its repressor activity. Conclusions/Significance: These studies thus open up avenues for understanding the unique features of MtBPL and the role it plays in biotin utilization in M. tuberculosis.

Recognition of Interaction Interface Residues in Low-Resolution Structures of Protein Assemblies Solely from the Positions of C alpha Atoms

Background: The number of available structures of large multi-protein assemblies is quite small. Such structures provide phenomenal insights on the organization, mechanism of formation and functional properties of the assembly. Hence detailed analysis of such structures is highly rewarding. However, the common problem in such analyses is the low resolution of these structures. In the recent times a number of attempts that combine low resolution cryo-EM data with higher resolution structures determined using X-ray analysis or NMR or generated using comparative modeling have been reported. Even in such attempts the best result one arrives at is the very course idea about the assembly structure in terms of trace of the C alpha atoms which are modeled with modest accuracy. Methodology/Principal Findings: In this paper first we present an objective approach to identify potentially solvent exposed and buried residues solely from the position of C alpha atoms and amino acid sequence using residue type-dependent thresholds for accessible surface areas of C alpha. We extend the method further to recognize potential protein-protein interface residues. Conclusion/Significance: Our approach to identify buried and exposed residues solely from the positions of C alpha atoms resulted in an accuracy of 84%, sensitivity of 83-89% and specificity of 67-94% while recognition of interfacial residues corresponded to an accuracy of 94%, sensitivity of 70-96% and specificity of 58-94%. Interestingly, detailed analysis of cases of mismatch between recognition of interface residues from C alpha positions and all-atom models suggested that, recognition of interfacial residues using C alpha atoms only correspond better with intuitive notion of what is an interfacial residue. Our method should be useful in the objective analysis of structures of protein assemblies when positions of only C alpha positions are available as, for example, in the cases of integration of cryo-EM data and high resolution structures of the components of the assembly.

PIM2 Induced COX-2 and MMP-9 Expression in Macrophages Requires PI3K and Notch1 Signaling

Activation of inflammatory immune responses during granuloma formation by the host upon infection of mycobacteria is one of the crucial steps that is often associated with tissue remodeling and breakdown of the extracellular matrix. In these complex processes, cyclooxygenase-2 (COX-2) plays a major role in chronic inflammation and matrix metalloproteinase-9 (MMP-9) significantly in tissue remodeling. In this study, we investigated the molecular mechanisms underlying Phosphatidyl-myo-inositol dimannosides (PIM2), an integral component of the mycobacterial envelope, triggered COX-2 and MMP-9 expression in macrophages. PIM2 triggers the activation of Phosphoinositide-3 Kinase (PI3K) and Notch1 signaling leading to COX-2 and MMP-9 expression in a Toll-like receptor 2 (TLR2)-MyD88 dependent manner. Notch1 signaling perturbations data demonstrate the involvement of the cross-talk with members of PI3K and Mitogen activated protein kinase pathway. Enforced expression of the cleaved Notch1 in macrophages induces the expression of COX-2 and MMP-9. PIM2 triggered significant p65 nuclear factor-kappa B (NF-kappa B) nuclear translocation that was dependent on activation of PI3K or Notch1 signaling. Furthermore, COX-2 and MMP-9 expression requires Notch1 mediated recruitment of uppressor of Hairless (CSL) and NF-kappa B to respective promoters. Inhibition of PIM2 induced COX-2 resulted in marked reduction in MMP-9 expression clearly implicating the role of COX-2 dependent signaling events in driving the MMP-9 expression. Taken together, these data implicate PI3K and Notch1 signaling as obligatory early proximal signaling events during PIM2 induced COX-2 and MMP-9 expression in macrophages.

Spiral-Wave Turbulence and Its Control in the Presence of Inhomogeneities in Four Mathematical Models of Cardiac Tissue

Regular electrical activation waves in cardiac tissue lead to the rhythmic contraction and expansion of the heart that ensures blood supply to the whole body. Irregularities in the propagation of these activation waves can result in cardiac arrhythmias, like ventricular tachycardia (VT) and ventricular fibrillation (VF), which are major causes of death in the industrialised world. Indeed there is growing consensus that spiral or scroll waves of electrical activation in cardiac tissue are associated with VT, whereas, when these waves break to yield spiral- or scroll-wave turbulence, VT develops into life-threatening VF: in the absence of medical intervention, this makes the heart incapable of pumping blood and a patient dies in roughly two-and-a-half minutes after the initiation of VF. Thus studies of spiral- and scroll-wave dynamics in cardiac tissue pose important challenges for in vivo and in vitro experimental studies and for in silico numerical studies of mathematical models for cardiac tissue. A major goal here is to develop low-amplitude defibrillation schemes for the elimination of VT and VF, especially in the presence of inhomogeneities that occur commonly in cardiac tissue. We present a detailed and systematic study of spiral- and scroll-wave turbulence and spatiotemporal chaos in four mathematical models for cardiac tissue, namely, the Panfilov, Luo-Rudy phase 1 (LRI), reduced Priebe-Beuckelmann (RPB) models, and the model of ten Tusscher, Noble, Noble, and Panfilov (TNNP). In particular, we use extensive numerical simulations to elucidate the interaction of spiral and scroll waves in these models with conduction and ionic inhomogeneities; we also examine the suppression of spiral- and scroll-wave turbulence by low-amplitude control pulses. Our central qualitative result is that, in all these models, the dynamics of such spiral waves depends very sensitively on such inhomogeneities. We also study two types of control chemes that have been suggested for the control of spiral turbulence, via low amplitude current pulses, in such mathematical models for cardiac tissue; our investigations here are designed to examine the efficacy of such control schemes in the presence of inhomogeneities. We find that a local pulsing scheme does not suppress spiral turbulence in the presence of inhomogeneities; but a scheme that uses control pulses on a spatially extended mesh is more successful in the elimination of spiral turbulence. We discuss the theoretical and experimental implications of our study that have a direct bearing on defibrillation, the control of life-threatening cardiac arrhythmias such as ventricular fibrillation.

Phenotypic and Functional Characterization of Human Mammary Stem/Progenitor Cells in Long Term Culture

Background: Cancer stem cells exhibit close resemblance to normal stem cells in phenotype as well as function. Hence, studying normal stem cell behavior is important in understanding cancer pathogenesis. It has recently been shown that human breast stem cells can be enriched in suspension cultures as mammospheres. However, little is known about the behavior of these cells in long-term cultures. Since extensive self-renewal potential is the hallmark of stem cells, we undertook a detailed phenotypic and functional characterization of human mammospheres over long-term passages. Methodology: Single cell suspensions derived from human breast `organoids' were seeded in ultra low attachment plates in serum free media. Resulting primary mammospheres after a week (termed T1 mammospheres) were subjected to passaging every 7th day leading to the generation of T2, T3, and T4 mammospheres. Principal Findings: We show that primary mammospheres contain a distinct side-population (SP) that displays a CD24(low)/CD44(low) phenotype, but fails to generate mammospheres. Instead, the mammosphere-initiating potential rests within the CD44(high)/CD24(low) cells, in keeping with the phenotype of breast cancer-initiating cells. In serial sphere formation assays we find that even though primary (T1) mammospheres show telomerase activity and fourth passage T4 spheres contain label-retaining cells, they fail to initiate new mammospheres beyond T5. With increasing passages, mammospheres showed an increase in smaller sized spheres, reduction in proliferation potential and sphere forming efficiency, and increased differentiation towards the myoepithelial lineage. Significantly, staining for senescence-associated beta-galactosidase activity revealed a dramatic increase in the number of senescent cells with passage, which might in part explain the inability to continuously generate mammospheres in culture. Conclusions: Thus, the self-renewal potential of human breast stem cells is exhausted within five in vitro passages of mammospheres, suggesting the need for further improvisation in culture conditions for their long-term maintenance.

Location of Pathogenic Bacteria during Persistent Infections: Insights from an Analysis Using Game Theory

Bacterial persistent infections are responsible for a significant amount of the human morbidity and mortality. Unlike acute bacterial infections, it is very difficult to treat persistent bacterial infections (e.g. tuberculosis). Knowledge about the location of pathogenic bacteria during persistent infection will help to treat such conditions by designing novel drugs which can reach such locations. In this study, events of bacterial persistent infections were analyzed using game theory. A game was defined where the pathogen and the host are the two players with a conflict of interest. Criteria for the establishment of Nash equilibrium were calculated for this game. This theoretical model, which is very simple and heuristic, predicts that during persistent infections pathogenic bacteria stay in both intracellular and extracellular compartments of the host. The result of this study implies that a bacterium should be able to survive in both intracellular and extracellular compartments of the host in order to cause persistent infections. This explains why persistent infections are more often caused by intracellular pathogens like Mycobacterium and Salmonella. Moreover, this prediction is in consistence with the results of previous experimental studies.

Structure-Based Phylogeny as a Diagnostic for Functional Characterization of Proteins with a Cupin Fold

Background: The members of cupin superfamily exhibit large variations in their sequences, functions, organization of domains, quaternary associations and the nature of bound metal ion, despite having a conserved beta-barrel structural scaffold. Here, an attempt has been made to understand structure-function relationships among the members of this diverse superfamily and identify the principles governing functional diversity. The cupin superfamily also contains proteins for which the structures are available through world-wide structural genomics initiatives but characterized as ``hypothetical''. We have explored the feasibility of obtaining clues to functions of such proteins by means of comparative analysis with cupins of known structure and function. Methodology/Principal Findings: A 3-D structure-based phylogenetic approach was undertaken. Interestingly, a dendrogram generated solely on the basis of structural dissimilarity measure at the level of domain folds was found to cluster functionally similar members. This clustering also reflects an independent evolution of the two domains in bicupins. Close examination of structural superposition of members across various functional clusters reveals structural variations in regions that not only form the active site pocket but are also involved in interaction with another domain in the same polypeptide or in the oligomer. Conclusions/Significance: Structure-based phylogeny of cupins can influence identification of functions of proteins of yet unknown function with cupin fold. This approach can be extended to other proteins with a common fold that show high evolutionary divergence. This approach is expected to have an influence on the function annotation in structural genomics initiatives.

lac Repressor Is an Antivirulence Factor of Salmonella enterica: Its Role in the Evolution of Virulence in Salmonella

The genus Salmonella includes many pathogens of great medical and veterinary importance. Bacteria belonging to this genus are very closely related to those belonging to the genus Escherichia. lacZYA operon and lacI are present in Escherichia coli, but not in Salmonella enterica. It has been proposed that Salmonella has lost lacZYA operon and lacI during evolution. In this study, we have investigated the physiological and evolutionary significance of the absence of lacI in Salmonella enterica. Using murine model of typhoid fever, we show that the expression of Lacl causes a remarkable reduction in the virulence of Salmonella enterica. Lacl also suppresses the ability of Salmonella enterica to proliferate inside murine macrophages. Microarray analysis revealed that Lacl interferes with the expression of virulence genes of Salmonella pathogenicity island 2. This effect was confirmed by RT-PCR and Western blot analysis. Interestingly, we found that SBG0326 of Salmonella bongori is homologous to lacI of Escherichia coli. Salmonella bongori is the only other species of the genus Salmonella and it lacks the virulence genes of Salmonella pathogenicity island 2. Overall, our results demonstrate that Lacl is an antivirulence factor of Salmonella enterica and suggest that absence of lacI has facilitated the acquisition of virulence genes of Salmonella pathogenicity island 2 in Salmonella enterica making it a successful systemic pathogen.

A Search for Energy Minimized Sequences of Proteins

In this paper, we present numerical evidence that supports the notion of minimization in the sequence space of proteins for a target conformation. We use the conformations of the real proteins in the Protein Data Bank (PDB) and present computationally efficient methods to identify the sequences with minimum energy. We use edge-weighted connectivity graph for ranking the residue sites with reduced amino acid alphabet and then use continuous optimization to obtain the energy-minimizing sequences. Our methods enable the computation of a lower bound as well as a tight upper bound for the energy of a given conformation. We validate our results by using three different inter-residue energy matrices for five proteins from protein data bank (PDB), and by comparing our energy-minimizing sequences with 80 million diverse sequences that are generated based on different considerations in each case. When we submitted some of our chosen energy-minimizing sequences to Basic Local Alignment Search Tool (BLAST), we obtained some sequences from non-redundant protein sequence database that are similar to ours with an E-value of the order of 10(-7). In summary, we conclude that proteins show a trend towards minimizing energy in the sequence space but do not seem to adopt the global energy-minimizing sequence. The reason for this could be either that the existing energy matrices are not able to accurately represent the inter-residue interactions in the context of the protein environment or that Nature does not push the optimization in the sequence space, once it is able to perform the function.

Mapping the risk of soil-transmitted helminthic infections in the Philippines

Background In order to increase the efficient allocation of soil-transmitted helminth (STH) disease control resources in the Philippines, we aimed to describe for the first time the spatial variation in the prevalence of A. lumbricoides, T. trichiura and hookworm across the country, quantify the association between the physical environment and spatial variation of STH infection and develop predictive risk maps for each infection. Methodology/Principal Findings Data on STH infection from 35,573 individuals across the country were geolocated at the barangay level and included in the analysis. The analysis was stratified geographically in two major regions: 1) Luzon and the Visayas and 2) Mindanao. Bayesian geostatistical models of STH prevalence were developed, including age and sex of individuals and environmental variables (rainfall, land surface temperature and distance to inland water bodies) as predictors, and diagnostic uncertainty was incorporated. The role of environmental variables was different between regions of the Philippines. This analysis revealed that while A. lumbricoides and T. trichiura infections were widespread and highly endemic, hookworm infections were more circumscribed to smaller foci in the Visayas and Mindanao. Conclusions/Significance This analysis revealed significant spatial variation in STH infection prevalence within provinces of the Philippines. This suggests that a spatially targeted approach to STH interventions, including mass drug administration, is warranted. When financially possible, additional STH surveys should be prioritized to high-risk areas identified by our study in Luzon.

Characterization of the Entamoeba histolytica Ornithine Decarboxylase-Like Enzyme

Background: The polyamines putrescine, spermidine, and spermine are organic cations that are required for cell growth and differentiation. Ornithine decarboxylase (ODC), the first and rate-limiting enzyme in the polyamine biosynthetic pathway, is a highly regulated enzyme. Methodology and Results: To use this enzyme as a potential drug target, the gene encoding putative ornithine decarboxylase (ODC)-like sequence was cloned from Entamoeba histolytica, a protozoan parasite causing amoebiasis. DNA sequence analysis revealed an open reading frame (ORF) of similar to 1,242 bp encoding a putative protein of 413 amino acids with a calculated molecular mass of 46 kDa and a predicted isoelectric point of 5.61. The E. histolytica putative ODC-like sequence has 33% sequence identity with human ODC and 36% identity with the Datura stramonium ODC. The ORF is a single-copy gene located on a 1.9-Mb chromosome. The recombinant putative ODC protein (48 kDa) from E. histolytica was heterologously expressed in Escherichia coli. Antiserum against recombinant putative ODC protein detected a band of anticipated size similar to 46 kDa in E. histolytica whole-cell lysate. Difluoromethylornithine (DFMO), an enzyme-activated irreversible inhibitor of ODC, had no effect on the recombinant putative ODC from E. histolytica. Comparative modeling of the three-dimensional structure of E. histolytica putative ODC shows that the putative binding site for DFMO is disrupted by the substitution of three amino acids-aspartate-332, aspartate-361, and tyrosine-323-by histidine-296, phenylalanine-305, and asparagine-334, through which this inhibitor interacts with the protein. Amino acid changes in the pocket of the E. histolytica enzyme resulted in low substrate specificity for ornithine. It is possible that the enzyme has evolved a novel substrate specificity. Conclusion: To our knowledge this is the first report on the molecular characterization of putative ODC-like sequence from E. histolytica. Computer modeling revealed that three of the critical residues required for binding of DFMO to the ODC enzyme are substituted in E. histolytica, resulting in the likely loss of interactions between the enzyme and DFMO.

Bayesian estimation of small effects in exercise and sports science

The aim of this paper is to provide a Bayesian formulation of the so-called magnitude-based inference approach to quantifying and interpreting effects, and in a case study example provide accurate probabilistic statements that correspond to the intended magnitude-based inferences. The model is described in the context of a published small-scale athlete study which employed a magnitude-based inference approach to compare the effect of two altitude training regimens (live high-train low (LHTL), and intermittent hypoxic exposure (IHE)) on running performance and blood measurements of elite triathletes. The posterior distributions, and corresponding point and interval estimates, for the parameters and associated effects and comparisons of interest, were estimated using Markov chain Monte Carlo simulations. The Bayesian analysis was shown to provide more direct probabilistic comparisons of treatments and able to identify small effects of interest. The approach avoided asymptotic assumptions and overcame issues such as multiple testing. Bayesian analysis of unscaled effects showed a probability of 0.96 that LHTL yields a substantially greater increase in hemoglobin mass than IHE, a 0.93 probability of a substantially greater improvement in running economy and a greater than 0.96 probability that both IHE and LHTL yield a substantially greater improvement in maximum blood lactate concentration compared to a Placebo. The conclusions are consistent with those obtained using a ‘magnitude-based inference’ approach that has been promoted in the field. The paper demonstrates that a fully Bayesian analysis is a simple and effective way of analysing small effects, providing a rich set of results that are straightforward to interpret in terms of probabilistic statements.

The effects of community attachment and information seeking on displaced disaster victims’ decision making

This paper uses original survey data of the Great East Japan earthquake disaster victims to examine their decision to apply for the temporary housing as well as the timing of application. We assess the effects of victims’ attachment to their locality as well as variation in victims’ information seeking behavior. We additionally consider various factors such as income, age, employment and family structure that are generally considered to affect the decision to choose temporary housing as victims’ solution for their displacement. Empirical results indicate that, ceteris paribus, as the degree of attachment increases, victims are more likely to apply for the temporary housing but attachment does not affect the timing of application. On the other hand, the victims who actively seek information and are able to collect higher quality information are less likely to apply for the temporary housing and if they do apply then they apply relatively later.

Transcription of Human Resistin Gene Involves an Interaction of Sp1 with Peroxisome Proliferator-Activating Receptor Gamma (PPAR gamma)

Background: Resistin is a cysteine rich protein, mainly expressed and secreted by circulating human mononuclear cells. While several factors responsible for transcription of mouse resistin gene have been identified, not much is known about the factors responsible for the differential expression of human resistin.Methodology/Principal Finding: We show that the minimal promoter of human resistin lies within similar to 80 bp sequence upstream of the transcriptional start site (-240) whereas binding sites for cRel, CCAAT enhancer binding protein alpha (C/EBP-alpha), activating transcription factor 2 (ATF-2) and activator protein 1 (AP-1) transcription factors, important for induced expression, are present within sequences up to -619. Specificity Protein 1(Sp1) binding site (-276 to -295) is also present and an interaction of Sp1 with peroxisome proliferator activating receptor gamma (PPAR gamma) is necessary for constitutive expression in U937 cells. Indeed co-immunoprecipitation assay demonstrated a direct physical interaction of Sp1 with PPAR gamma in whole cell extracts of U937 cells. Phorbol myristate acetate (PMA) upregulated the expression of resistin mRNA in U937 cells by increasing the recruitment of Sp1, ATF-2 and PPAR gamma on the resistin gene promoter. Furthermore, PMA stimulation of U937 cells resulted in the disruption of Sp1 and PPAR gamma interaction. Chromatin immunoprecipitation (ChIP) assay confirmed the recruitment of transcription factors phospho ATF-2, Sp1, Sp3, PPAR gamma, chromatin modifier histone deacetylase 1 (HDAC1) and the acetylated form of histone H3 but not cRel, C/EBP-alpha and phospho c-Jun during resistingene transcription.Conclusion: Our findings suggest a complex interplay of Sp1 and PPAR gamma along with other transcription factors that drives the expression of resistin in human monocytic U937 cells.