Speeding disease gene discovery by sequence based candidate prioritization

BACKGROUND: Regions of interest identified through genetic linkage studies regularly exceed 30 centimorgans in size and can contain hundreds of genes. Traditionally this number is reduced by matching functional annotation to knowledge of the disease or phenotype in question. However, here we show that disease genes share patterns of sequence-based features that can provide a good basis for automatic prioritization of candidates by machine learning. RESULTS: We examined a variety of sequence-based features and found that for many of them there are significant differences between the sets of genes known to be involved in human hereditary disease and those not known to be involved in disease. We have created an automatic classifier called PROSPECTR based on those features using the alternating decision tree algorithm which ranks genes in the order of likelihood of involvement in disease. On average, PROSPECTR enriches lists for disease genes two-fold 77% of the time, five-fold 37% of the time and twenty-fold 11% of the time. CONCLUSION: PROSPECTR is a simple and effective way to identify genes involved in Mendelian and oligogenic disorders. It performs markedly better than the single existing sequence-based classifier on novel data. PROSPECTR could save investigators looking at large regions of interest time and effort by prioritizing positional candidate genes for mutation detection and case-control association studies

Parallelism in eco-morphology and gene expression despite variable evolutionary and genomic backgrounds in a Holarctic fish

Understanding the extent to which ecological divergence is repeatable is essential for predicting responses of biodiversity to environmental change. Here we test the predictability of evolution, from genotype to phenotype, by studying parallel evolution in a salmonid fish, Arctic charr (Salvelinus alpinus), across eleven replicate sympatric ecotype pairs (benthivorous-planktivorous and planktivorous-piscivorous) and two evolutionary lineages. We found considerable variability in eco-morphological divergence, with several traits related to foraging (eye diameter, pectoral fin length) being highly parallel even across lineages. This suggests repeated and predictable adaptation to environment. Consistent with ancestral genetic variation, hundreds of loci were associated with ecotype divergence within lineages of which eight were shared across lineages. This shared genetic variation was maintained despite variation in evolutionary histories, ranging from postglacial divergence in sympatry (ca. 10-15kya) to pre-glacial divergence (ca. 20-40kya) with postglacial secondary contact. Transcriptome-wide gene expression (44,102 genes) was highly parallel across replicates, involved biological processes characteristic of ecotype morphology and physiology, and revealed parallelism at the level of regulatory networks. This expression divergence was not only plastic but in part genetically controlled by parallel cis-eQTL. Lastly, we found that the magnitude of phenotypic divergence was largely correlated with the genetic differentiation and gene expression divergence. In contrast, the direction of phenotypic change was mostly determined by the interplay of adaptive genetic variation, gene expression, and ecosystem size. Ecosystem size further explained variation in putatively adaptive, ecotype-associated genomic patterns within and across lineages, highlighting the role of environmental variation and stochasticity in parallel evolution. Together, our findings demonstrate the parallel evolution of eco-morphology and gene expression within and across evolutionary lineages, which is controlled by the interplay of environmental stochasticity and evolutionary contingencies, largely overcoming variable evolutionary histories and genomic backgrounds

Defining Electron Bifurcation in the Electron-Transferring Flavoprotein Family

Electron bifurcation is the coupling of exergonic and endergonic redox reactions to simultaneously generate (or utilize) low- and high-potential electrons. It is the third recognized form of energy conservation in biology and was recently described for select electron-transferring flavoproteins (Etfs). Etfs are flavin-containing heterodimers best known for donating electrons derived from fatty acid and amino acid oxidation to an electron transfer respiratory chain via Etf-quinone oxidoreductase. Canonical examples contain a flavin adenine dinucleotide (FAD) that is involved in electron transfer, as well as a non-redox-active AMP. However, Etfs demonstrated to bifurcate electrons contain a second FAD in place of the AMP. To expand our understanding of the functional variety and metabolic significance of Etfs and to identify amino acid sequence motifs that potentially enable electron bifurcation, we compiled 1,314 Etf protein sequences from genome sequence databases and subjected them to informatic and structural analyses. Etfs were identified in diverse archaea and bacteria, and they clustered into five distinct well-supported groups, based on their amino acid sequences. Gene neighborhood analyses indicated that these Etf group designations largely correspond to putative differences in functionality. Etfs with the demonstrated ability to bifurcate were found to form one group, suggesting that distinct conserved amino acid sequence motifs enable this capability. Indeed, structural modeling and sequence alignments revealed that identifying residues occur in the NADH- and FAD-binding regions of bifurcating Etfs. Collectively, a new classification scheme for Etf proteins that delineates putative bifurcating versus nonbifurcating members is presented and suggests that Etf-mediated bifurcation is associated with surprisingly diverse enzymes

Rac1 as a therapeutic target in ovarian cancer

Rac1 is a high value therapeutic target for cancer based on its tumor promoting activities, yet clinical applications targeting Rac1 are in their infancy. High expression and hyperactivation of Rac1 in ovarian cancer, along with our identification of R-ketorolac as a novel Rac1 and Cdc42 selective inhibitor with translational potential, prompt us to test the hypothesis that targeting Rac1 has therapeutic utility for ovarian cancer. Ascites tumor cell samples from ovarian cancer patients in a prospective study receiving racemic ketorolac for clinically indicated use in pain relief were previously reported to show time dependent reduction of Rac1 and Cdc42 activities post-treatment. New RNA seq data of these patient samples reveals significant changes of genes involved in cell adhesion, cytokine-mediated signaling and cytokine production pathways. Conversely, the identified downregulated genes were overexpressed and associated with worse survival in ovarian cancer patients analyzed through The Cancer Genome Atlas (TCGA). Among the downregulated genes in the NOD pathway are chemokines and pro-inflammatory cytokines. Follow-up cytokine panels from patients confirm that racemic ketorolac treatment reduces the levels of immunosuppressive cytokines IL-6, IL-10 and RANTES in ascites fluids. Together, these data indicate there may be a benefit to the anti-inflammatory activity of the S- enantiomer, as well as the GTPase inhibitory activity of the R- enantiomer of ketorolac for ovarian cancer treatment

Trauma history and depression predict incomplete adherence to antiretroviral therapies in a low income country.

As antiretroviral therapy (ART) for HIV becomes increasingly available in low and middle income countries (LMICs), understanding reasons for lack of adherence is critical to stemming the tide of infections and improving health. Understanding the effect of psychosocial experiences and mental health symptomatology on ART adherence can help maximize the benefit of expanded ART programs by indicating types of services, which could be offered in combination with HIV care. The Coping with HIV/AIDS in Tanzania (CHAT) study is a longitudinal cohort study in the Kilimanjaro Region that included randomly selected HIV-infected (HIV+) participants from two local hospital-based HIV clinics and four free-standing voluntary HIV counselling and testing sites. Baseline data were collected in 2008 and 2009; this paper used data from 36 month follow-up interviews (N = 468). Regression analyses were used to predict factors associated with incomplete self-reported adherence to ART. INCOMPLETE ART ADHERENCE WAS SIGNIFICANTLY MORE LIKELY TO BE REPORTED AMONGST PARTICIPANTS WHO EXPERIENCED A GREATER NUMBER OF CHILDHOOD TRAUMATIC EVENTS: sexual abuse prior to puberty and the death in childhood of an immediate family member not from suicide or homicide were significantly more likely in the non-adherent group and other negative childhood events trended toward being more likely. Those with incomplete adherence had higher depressive symptom severity and post-traumatic stress disorder (PTSD). In multivariable analyses, childhood trauma, depression, and financial sacrifice remained associated with incomplete adherence.\ud This is the first study to examine the effect of childhood trauma, depression and PTSD on HIV medication adherence in a low income country facing a significant burden of HIV. Allocating spending on HIV/AIDS toward integrating mental health services with HIV care is essential to the creation of systems that enhance medication adherence and maximize the potential of expanded antiretroviral access to improve health and reduce new infections

Motor development in infancy and spine shape in early old age: findings from a British birth cohort study

Spine shape changes dramatically in early life, influenced by attainment of developmental milestones such as independent walking. Whether these associations persist across life is unknown. Therefore, we investigated associations between developmental milestones and spine shape, as determined using statistical shape models (SSMs) of lumbar spine from DXA scans in 1327 individuals (688 female) at 60‐64y in the MRC National Survey of Health and Development. Lumbar lordosis angle (L4 inferior endplate to T12 superior endplate) was measured using the two‐line Cobb method. In analyses adjusted for sex, height, lean and fat mass, socioeconomic position and birthweight, later walking age was associated with greater lordosis described by SSM1 (regression coefficient 0.023, 95%CI 0.000‐0.047, p=0.05) and direct angle measurement. Modest associations between walking age and less variation in anterior‐posterior vertebral size caudally (SSM6) were also observed (0.021, 95%CI ‐0.002‐0.044, p=0.07). Sex interactions showed that later walking was associated with larger relative vertebral anterior‐posterior dimensions in men (SSM3; ‐0.043, 95%CI ‐0.075‐0.01, p=0.01) but not women (0.018, 95%CI ‐0.0007‐0.043, p=0.17). Similar associations were observed between age at independent standing and SSMs but there was little evidence of association between sitting age and spine shape. Unadjusted associations between walking age and SSMs 1 and 6 remained similar after adjustment for potential confounders and mediators. This suggests that these associations may be explained by altered mechanical loading of the spine during childhood growth, although other factors could contribute. Early life motor development, particularly walking, may have a lasting effect on features of spine morphology with clinical significance

Age at onset of walking in infancy is associated with hip shape in early old age

Bones' shapes and structures adapt to the muscle and reaction forces they experience during everyday movements. Onset of independent walking, at approximately 12 months, represents the first postnatal exposure of the lower limbs to the large forces associated with bipedal movements, and, accordingly, earlier walking is associated with greater bone strength. However, associations between early life loading and joint shape have not been explored. We therefore examined associations between walking age and hip shape at age 60-64y in 1423 individuals (740 women) from the MRC National Survey of Health and Development, a nationally-representative British birth cohort. Walking age in months was obtained from maternal interview at age 2y. Ten modes of variation in hip shape (HM1-HM10), described by statistical shape models, were ascertained from dual-energy X-ray absorptiometry (DXA) images. In sex-adjusted analyses, earlier walking age was associated with higher HM1 and HM7 scores; these associations were maintained after further adjustment for height, body composition and socioeconomic position. Earlier walking was also associated with lower HM2 scores in women only, and lower HM4 scores in men only. Taken together, this suggests that earlier walkers have proportionately larger (HM4) and flatter (HM1,4) femoral heads, wider (HM1,4,7) and flatter (HM1, 7) femoral necks, smaller neck-shaft angle (HM1,4), anteversion (HM2,7) and development of osteophytes (HM1). These results suggest that age at onset of walking in infancy is associated with variations in hip shape in older age. Early walkers have a larger femoral head and neck and smaller neck-shaft angle; these features are associated with reduced hip fracture risk, but also represent an osteoarthritic-like phenotype. Unlike results of previous studies of walking age and bone mass, associations in this study were not affected by adjustment for lean mass suggesting that associations may relate directly to skeletal loading in early life when joint shape changes rapidly. This article is protected by copyright. All rights reserved.</p

International Lessons in New Methods for Grading and Integrating Cost Effectiveness Evidence into Clinical Practice Guidelines

Economic evidence is influential in health technology assessment world-wide. Clinical Practice Guidelines (CPG) can enable economists to include economic information on health care provision. Application of economic evidence in CPGs, and its integration into clinical practice and national decision making is hampered by objections from professions, paucity of economic evidence or lack of policy commitment. The use of state-of-art economic methodologies will improve this. Economic evidence can be graded by 'checklists' to establish the best evidence for decision making given methodological rigor. New economic evaluation checklists, Multi-Criteria Decision Analyses (MCDA) and other decision criteria enable health economists to impact on decision making world-wide. We analyse the methodologies for integrating economic evidence into CPG agencies globally, including the Agency of Health Research and Quality (AHRQ) in the USA, National Health and Medical Research Council (NHMRC) and Australian political reforms. The Guidelines and Economists Network International (GENI) Board members from Australia, UK, Canada and Denmark presented the findings at the conference of the International Health Economists Association (IHEA) and we report conclusions and developments since. The Consolidated Guidelines for the Reporting of Economic Evaluations (CHEERS) 24 item check list can be used by AHRQ, NHMRC, other CPG and health organisations, in conjunction with the Drummond ten-point check list and a questionnaire that scores that checklist for grading studies, when assessing economic evidence. Cost-effectiveness Analysis (CEA) thresholds, opportunity cost and willingness-to-pay (WTP) are crucial issues for decision rules in CEA generally, including end-of-life therapies. Limitations of inter-rater reliability in checklists can be addressed by including more than one assessor to reach a consensus, especially when impacting on treatment decisions. We identify priority areas to generate economic evidence for CPGs by NHMRC, AHRQ, and other agencies. The evidence may cover demand for care issues such as involved time, logistics, innovation price, price sensitivity, substitutes and complements, WTP, absenteeism and presentism. Supply issues may include economies of scale, efficiency changes, and return on investment. Involved equity and efficiency measures may include cost-of-illness, disease burden, quality-of-life, budget impact, cost-effective ratios, net benefits and disparities in access and outcomes.. Priority setting remains essential and trade-off decisions between policy criteria can be based on MCDA, both in evidence based clinical medicine and in health planning

Ribosomal oxygenases are structurally conserved from prokaryotes to humans

2-Oxoglutarate (2OG)-dependent oxygenases have important roles in the regulation of gene expression via demethylation of N-methylated chromatin components1,2 and in the hydroxylation of transcription factors3 and splicing factor proteins4. Recently, 2OG-dependent oxygenases that catalyse hydroxylation of transfer RNA5,6,7 and ribosomal proteins8 have been shown to be important in translation relating to cellular growth, TH17-cell differentiation and translational accuracy9,10,11,12. The finding that ribosomal oxygenases (ROXs) occur in organisms ranging from prokaryotes to humans8 raises questions as to their structural and evolutionary relationships. In Escherichia coli, YcfD catalyses arginine hydroxylation in the ribosomal protein L16; in humans, MYC-induced nuclear antigen (MINA53; also known as MINA) and nucleolar protein 66 (NO66) catalyse histidine hydroxylation in the ribosomal proteins RPL27A and RPL8, respectively. The functional assignments of ROXs open therapeutic possibilities via either ROX inhibition or targeting of differentially modified ribosomes. Despite differences in the residue and protein selectivities of prokaryotic and eukaryotic ROXs, comparison of the crystal structures of E. coli YcfD and Rhodothermus marinus YcfD with those of human MINA53 and NO66 reveals highly conserved folds and novel dimerization modes defining a new structural subfamily of 2OG-dependent oxygenases. ROX structures with and without their substrates support their functional assignments as hydroxylases but not demethylases, and reveal how the subfamily has evolved to catalyse the hydroxylation of different residue side chains of ribosomal proteins. Comparison of ROX crystal structures with those of other JmjC-domain-containing hydroxylases, including the hypoxia-inducible factor asparaginyl hydroxylase FIH and histone Nε-methyl lysine demethylases, identifies branch points in 2OG-dependent oxygenase evolution and distinguishes between JmjC-containing hydroxylases and demethylases catalysing modifications of translational and transcriptional machinery. The structures reveal that new protein hydroxylation activities can evolve by changing the coordination position from which the iron-bound substrate-oxidizing species reacts. This coordination flexibility has probably contributed to the evolution of the wide range of reactions catalysed by oxygenases