Established and emerging strategies to crack the genetic code of obesity

Tremendous progress has been made in the genetic elucidation of obesity over the past two decades, driven largely by technological, methodological and organizational innovations. Current strategies for identifying obesity-predisposing loci/genes, including cytogenetics, linkage analysis, homozygosity mapping, admixture mapping, candidate gene studies, genome-wide association studies, custom genotyping arrays, whole-exome sequencing and targeted exome sequencing, have achieved differing levels of success, and the identified loci in aggregate explain only a modest fraction of the estimated heritability of obesity. This review outlines the successes and limitations of these approaches and proposes novel strategies, including the use of exceptionally large sample sizes, the study of diverse ethnic groups and deep phenotypes and the application of innovative methods and study designs, to identify the remaining obesity-predisposing genes. The use of both established and emerging strategies has the potential to crack the genetic code of obesity in the not-too-distant future. The resulting knowledge is likely to yield improvements in obesity prediction, prevention and care.

Computational modeling of early T-cell precursor acute lymphoblastic leukemia (ETP-ALL) to identify personalized therapy using genomics

Early T-cell precursor acute lymphoblastic leukemia (ETP-ALL) is an aggressive hematological malignancy for which optimal therapeutic approaches are poorly characterized. Using computational biology modeling (CBM) in conjunction with genomic data from cell lines and individual patients, we generated disease-specific protein network maps that were used to identify unique characteristics associated with the mutational profiles of ETP-ALL compared to non-ETP-ALL (T-ALL) cases and simulated cellular responses to a digital library of FDA-approved and investigational agents. Genomics-based classification of ETP-ALL patients using CBM had a prediction sensitivity and specificity of 93% and 87%, respectively. This analysis identified key genomic and pathway characteristics that are distinct in ETP-ALL including deletion of nucleophosmin-1 (NPM1), mutations of which are used to direct therapeutic decisions in acute myeloid leukemia. Computational simulations based on mutational profiles of 62 ETP-ALL patient models identified 87 unique targeted combination therapies in 56 of the 62 patients despite actionable mutations being present in only 37% of ETP-ALL patients. Shortlisted two-drug combinations were predicted to be synergistic in 11 profiles and were validated by in vitro chemosensitivity assays. In conclusion, computational modeling was able to identify unique biomarkers and pathways for ETP-ALL, and identify new drug combinations for potential clinical testing.

Revisiting the evolutionary origins of obesity: lazy versus peppy-thrifty genotype hypothesis

The recent global obesity epidemic is attributed to major societal and environmental changes, such as excessive energy intake and sedentary lifestyle. However, exposure to 'obesogenic' environments does not necessarily result in obesity at the individual level, as 40-75% of body mass index variation in population is attributed to genetic differences. The thrifty genotype theory posits that genetic variants promoting efficient food sequestering and optimal deposition of fat during periods of food abundance were evolutionarily advantageous for the early hunter-gatherer and were positively selected. However, the thrifty genotype is likely too simplistic and fails to provide a justification for the complex distribution of obesity predisposing gene variants and for the broad range of body mass index observed in diverse ethnic groups. This review proposes that gene pleiotropy may better account for the variability in the distribution of obesity susceptibility alleles across modern populations. We outline the lazy-thrifty versus peppy-thrifty genotype hypothesis and detail the body of evidence in the literature in support of this novel concept. Future population genetics and mathematical modelling studies that account for pleiotropy may further improve our understanding of the evolutionary origins of the current obesity epidemic.

Abstract 3253: Droplet digital PCR is a sensitive method for detecting refractory acute myeloid leukemia (AML) clones in peripheral blood and saliva

Droplet digital PCR (ddPCR) is a highly sensitive and rapid method for detecting mutant allele frequency (MAF). In preliminary work, our lower limit of detection for common myeloid gene mutations was 0.001% in peripheral blood and bone marrow compared to 0.1% with flow cytometry and 0.01% with real-time quantitative PCR, and turnaround time is 1 day. Furthermore, we detected leukemic mutant alleles in peripheral blood (PB), introducing the possibility of sparing painful bone marrow biopsy procedures to determine treatment response. Thus, we hypothesized that ddPCR is a feasible and accurate method for monitoring leukemic disease burden in PB for the prospective care of patients (pts) with AML. Eighteen patients (pts) with de novo, relapse/refractory, and secondary AML were recruited to an IRB-approved study (NCT02435550) and bone marrow (BM), peripheral blood (PB), and saliva were collected at standard clinical visits. Gene mutations were identified by whole-exome sequencing (WES) of BM specimens at study entry. For ddPCR interrogation, genomic DNA was isolated (Qiagen), and select primers and probes (Bio-Rad/IDT) were developed based on variants identified in WES data. Case-specific primers and probes were validated on archived specimens obtained at study entry. 12/18 pt mutanomes met criteria for primer/probe design. 8 pts are in the primer/probe design and validation stage and 4 have completed validation and serial analyses. WES identified, and ddPCR confirmed, at least 1 mutation per patient at the study entry timepoint. The mutations included NRAS G13R, NRAS G12A, CSF3R T618I, and IDH2 R172K. In 2 cases, we observed a reduction in both PB and saliva MAF that were consistent with the reduction in both BM and PB blasts after treatment, resulting in complete remissions. Although PB blasts were reduced in a third pt receiving ruxolitinib, the persistence of their CSF3R MAF in PB indicated a resistant AML clone. WES revealed the presence of NRAS G13R variant in a secondary AML pt; however, WES did not detect this NRAS G13R variant in a cryopreserved BM specimen obtained at the pts MDS diagnosis. Interestingly, ddPCR was able to detect NRAS G13R variant at 0.1% MAF in a PB sample obtained at MDS diagnosis, demonstrating the ultrasensitive detection of rare variants within a sample, and highlighting the subclonal evolution of this pt's malignancy. Rapid detection of myeloid-related somatic mutations in a variety of tissue sources (i.e., saliva, PB) will allow for noninvasive monitoring of AML tumor burden. ddPCR may be used to observe molecular response to treatment and to detect molecular residual disease and relapse prior to clinically indicated BM biopsies.

Citation Format: Kimberly E. Hawkins, Cesia Salan, Madeleine Turcotte, Lauren T. Vaughn, Mei Zhang, Yanping Zhang, Barry Sawicki, Glenda G. Anderson, Nosha Farhadfar, Hemant S. Murthy, Biljana N. Horn, Helen L. Leather, Paul Castillo, Maxim Norkin, John W. Hiemenz, Randy A. Brown, William Slayton, Jack W. Hsu, John R. Wingard, Christopher R. Cogle, Leylah M. Drusbosky. Droplet digital PCR is a sensitive method for detecting refractory acute myeloid leukemia (AML) clones in peripheral blood and saliva [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 3253.

Obesity genetics: insights from the Pakistani population

The Pakistani population is extensively diverse, indicating a genetic admixture of European and Central/West Asian migrants with indigenous South Asian gene pools. Pakistanis are organized in different ethnicities/castes based on cultural, linguistic and geographical origin. While Pakistan is facing a rapid nutritional transition, the rising prevalence of obesity is driving a growing burden of health complications and mortality. This represents a unique opportunity for the research community to study the interplay between obesogenic environmental changes and obesity predisposing genes in the time frame of one generation. This review recapitulates the ancestral origins of Pakistani population, the societal determinants of the rise in obesity and its governmental management. We describe the contribution of syndromic, monogenic non-syndromic and polygenic obesity genes identified in the Pakistani population. We then discuss the utility of gene identification approaches based on large consanguineous families and original gene × environment interaction study designs in discovering new obesity genes and causal pathways. Elucidation of the genetic basis of obesity in the Pakistani population may result in improved methods of obesity prevention and treatment globally.

On the origin of obesity: identifying the biological, environmental and cultural drivers of genetic risk among human populations

Genetic predisposition to obesity presents a paradox: how do genetic variants with a detrimental impact on human health persist through evolutionary time? Numerous hypotheses, such as the thrifty genotype hypothesis, attempt to explain this phenomenon yet fail to provide a justification for the modern obesity epidemic. In this critical review, we appraise existing theories explaining the evolutionary origins of obesity and explore novel biological and sociocultural agents of evolutionary change to help explain the modern-day distribution of obesity-predisposing variants. Genetic drift, acting as a form of 'blind justice,' may randomly affect allele frequencies across generations while gene pleiotropy and adaptations to diverse environments may explain the rise and subsequent selection of obesity risk alleles. As an adaptive response, epigenetic regulation of gene expression may impact the manifestation of genetic predisposition to obesity. Finally, exposure to malnutrition and disease epidemics in the wake of oppressive social systems, culturally mediated notions of attractiveness and desirability, and diverse mating systems may play a role in shaping the human genome. As an important first step towards the identification of important drivers of obesity gene evolution, this review may inform empirical research focused on testing evolutionary theories by way of population genetics and mathematical modelling.

Measurement of B --> K*gamma branching fractions and charge asymmetries

The branching fractions of the exclusive decays B0-->K(*0)gamma and B+-->K(*+)gamma are measured from a sample of (22.74+/-0.36)x10(6) BB decays collected with the BABAR detector at the PEP-II asymmetric e(+)e(-) collider. We find B (B0-->K(*0)gamma) = [4.23+/-0.40(stat)+/-0.22(syst)]x10(-5), B(B+-->K(*+)gamma) = [3.83+/-0.62(stat)+/-0.22(syst)]x10(-5) and constrain the CP-violating charge asymmetry to be -0.170<A(CP)(B-->K(*)gamma)<0.082 at 90% C.L.

Search for the decay B0-->gammagamma

We present a limit on the branching fraction for the decay B0-->gammagamma using data collected at the Upsilon(4S) resonance with the BABAR detector at the PEP-II asymmetric energy e+e- collider. Based on the observation of one event in the signal region, out of a sample of 21.3x10(6) e+e--->Upsilon(4S)-->BB decays, we establish an upper limit on the branching fraction of B(B0-->gammagamma)<1.7x10(-6) at the 90% confidence level. This result substantially improves upon existing limits.

Measurement of the B--> J/psiK*(892) decay amplitudes

We present a measurement of the decay amplitudes in B-->J/psiK*(892) channels using 20.7 fb(-1) of data collected at the Upsilon(4S) resonance with the BABAR detector at PEP-II. We measure a P-wave fraction R(perpendicular) = (16.0 +/- 3.2 +/- 1.4)% and a longitudinal polarization fraction (59.7 +/- 2.8 +/- 2.4)%. The measurement of a relative phase that is neither 0 nor pi, phi = 2.50 +/- 0.20 +/-0.08 radians, favors a departure from the factorization hypothesis. Although the decay B-->/psiK(pi) proceeds mainly via K*(892), there is also evidence for K2*(1430) and K(pi) S-wave contributions.

Generating protein three-dimensional fold signatures using inductive logic programming

Inductive logic programming (ILP) has been applied to automatically discover protein fold signatures. This paper investigates the use of topological information to circumvent problems encountered during previous experiments, namely (1) matching of non-structurally related secondary structures and (2) scaling problems. Cross-validation tests were carried out for 20 folds. The overall estimated accuracy is 73.37+/-0.35%. The new representation allows us to process the complete set of examples, while previously it was necessary to sample the negative examples. Topological information is used in approximately 90% of the rules presented here. Information about the topology of a sheet is present in 63% of the rules. This set of rules presents characteristics of the overall architecture of the fold. In contrast, 26% of the rules contain topological information which is limited to the packing of a restricted number of secondary structures, as such, the later set resembles those found in our previous studies.

Measurements of the branching fractions of exclusive charmless B meson decays with eta(') or omega mesons

We present the results of searches for B decays to charmless two-body final states containing eta(') or omega mesons, based on 20.7 fb(-1) of data collected with the BABAR detector. We find the branching fractions Beta(B(+)-->eta(')K(+)) = (70+/-8+/-5) x 10(-6), Beta(B(0)-->eta(')K(0)) = (42(+13)(-11) +/- 4) x 10(-6), and Beta(B(+)-->omega pi(+)) = (6.6(+2.1)(-1.8) +/- 0.7) x 10(-6), where the first error quoted is statistical and the second is systematic. We give measurements of four additional modes for which the 90% confidence level upper limits are Beta(B(+)-->eta(')pi(+)) < 12 x 10(-6), Beta(B(+)-->omega K(+)) < 4 x 10(-6), Beta(B(0)-->omega K(0)) < 13 x 10(-6), and Beta(B(0)-->omega pi(0)) < 3 x 10(-6).

Measurement of the B(0) and B(+) meson lifetimes with fully reconstructed hadronic final states

The B(0) and B(+) meson lifetimes have been measured in e(+)e(-) annihilation data collected in 1999 and 2000 with the BABAR detector at center-of-mass energies near the Upsilon(4S) resonance. Events are selected in which one B meson is fully reconstructed in a hadronic final state while the second B meson is reconstructed inclusively. A combined fit to the B(0) and the B(+) decay time difference distributions yields tau(B(0)) = 1.546+/-0.032(stat)+/-0.022(syst) ps, tau(B(+)) = 1.673+/-0.032(stat)+/-0.023(syst) ps, and tau(B(+))/tau(B(0)) = 1.082+/-0.026(stat)+/-0.012(syst).

Measurement of J/psi production in continuum e(+)e(-) annihilations near square root of s = 10.6 GeV

The production of J/psi mesons in continuum e(+)e(-) annihilations has been studied with the BABAR detector at energies near the Upsilon(4S) resonance. The mesons are distinguished from J/psi production in B decays through their center-of-mass momentum and energy. We measure the cross section e(+)e(-)-->J/psi X to be 2.52+/-0.21+/-0.21 pb. We set a 90% C.L. upper limit on the branching fraction for direct Upsilon(4S)-->J/psi X decays at 4.7 x 10(-4).

Measurement of the decays B--> phiK and B--> phiK*

We have observed the decays B--> phiK and phiK(*) in a sample of over 45 million B mesons collected with the BABAR detector at the PEP-II collider. The measured branching fractions are B(B+--> phiK+) = (7.7(+1.6)(-1.4)+/-0.8)x10(-6), B(B0--> phiK0) = (8.1(+3.1)(-2.5)+/-0.8)x10(-6), B(B+--> phiK(*+)) = (9.7(+4.2)(-3.4)+/-1.7)x10(-6), and B(B0--> phiK(*0)) = (8.7(+2.5)(-2.1)+/-1.1)x10(-6). We also report the upper limit B(B+--> phipi(+))<1.4x10(-6) ( 90% C.L.).

Measurement of branching fractions and search for CP-violating charge asymmetries in charmless two-body B decays into pions and kaons

We present measurements, based on a sample of approximately 23x10(6) BB pairs, of the branching fractions and a search for CP-violating charge asymmetries in charmless hadronic decays of B mesons into two-body final states of kaons and pions. We find the branching fractions B(B0-->pi(+)pi(-)) = (4.1+/-1.0+/-0.7)x10(-6), B(B0-->K+pi(-)) = (16.7+/-1.6+/-1.3)x10(-6), B(B+-->K+pi(0)) = (10.8(+2.1)(-1.9)+/-1.0)x10(-6), B(B+-->K0pi(+)) = (18.2(+3.3)(-3.0)+/-2.0)x10(-6), B(B0-->K0pi(0)) = (8.2(+3.1)(-2.7)+/-1.2)x10(-6). We also report 90% confidence level upper limits for B meson decays to the pi(+)pi(0), K+K-, and K0K+ final states. In addition, charge asymmetries have been found to be consistent with zero, where the statistical precision is in the range of +/-0.10 to +/-0.18, depending on the decay mode.

Observation of CP violation in the B(0) meson system

We present an updated measurement of time-dependent CP-violating asymmetries in neutral B decays with the BABAR detector at the PEP-II asymmetric B Factory at SLAC. This result uses an additional sample of Upsilon(4S) decays collected in 2001, bringing the data available to 32 x 10(6) BB macro pairs. We select events in which one neutral B meson is fully reconstructed in a final state containing charmonium and the flavor of the other neutral B meson is determined from its decay products. The amplitude of the CP-violating asymmetry, which in the standard model is proportional to sin2 beta, is derived from the decay time distributions in such events. The result sin2 beta = 0.59+/-0.14(stat)+/-0.05(syst) establishes CP violation in the B(0) meson system. We also determine absolute value of lambda = 0.93+/-0.09(stat)+/-0.03(syst), consistent with no direct CP violation.

Measurement of CP-violating asymmetries in B0 decays to CP eigenstates

We present measurements of time-dependent CP-violating asymmetries in neutral B decays to several CP eigenstates. The measurement uses a data sample of 23x10(6) Upsilon(4S)-->BbarB decays collected by the BABAR detector at the PEP-II asymmetric B Factory at SLAC. In this sample, we find events in which one neutral B meson is fully reconstructed in a CP eigenstate containing charmonium and the flavor of the other neutral B meson is determined from its decay products. The amplitude of the CP-violating asymmetry, which in the standard model is proportional to sin2beta, is derived from the decay time distributions in such events. The result is sin2beta = 0.34+/-0.20 (stat)+/-0.05 (syst).

Automated discovery of structural signatures of protein fold and function

There are constraints on a protein sequence/structure for it to adopt a particular fold. These constraints could be either a local signature involving particular sequences or arrangements of secondary structure or a global signature involving features along the entire chain. To search systematically for protein fold signatures, we have explored the use of Inductive Logic Programming (ILP). ILP is a machine learning technique which derives rules from observation and encoded principles. The derived rules are readily interpreted in terms of concepts used by experts. For 20 populated folds in SCOP, 59 rules were found automatically. The accuracy of these rules, which is defined as the number of true positive plus true negative over the total number of examples, is 74% (cross-validated value). Further analysis was carried out for 23 signatures covering 30% or more positive examples of a particular fold. The work showed that signatures of protein folds exist, about half of rules discovered automatically coincide with the level of fold in the SCOP classification. Other signatures correspond to homologous family and may be the consequence of a functional requirement. Examination of the rules shows that many correspond to established principles published in specific literature. However, in general, the list of signatures is not part of standard biological databases of protein patterns. We find that the length of the loops makes an important contribution to the signatures, suggesting that this is an important determinant of the identity of protein folds. With the expansion in the number of determined protein structures, stimulated by structural genomics initiatives, there will be an increased need for automated methods to extract principles of protein folding from coordinates.

A predicted consensus structure for the N-terminal fragment of the heat shock protein HSP90 family

A secondary structure has been predicted for the heat shock protein HSP90 family from an aligned set of homologous protein sequences by using a transparent method in both manual and automated implementation that extracts conformational information from patterns of variation and conservation within the family. No statistically significant sequence similarity relates this family to any protein with known crystal structure. However, the secondary structure prediction, together with the assignment of active site positions and possible biochemical properties, suggest that the fold is similar to that seen in N-terminal domain of DNA gyrase B (the ATPase fragment).

The combination of symbolic and numerical computation for three-dimensional modeling of RNA

Three-dimensional (3-D) structural models of RNA are essential for understanding of the cellular roles played by RNA. Such models have been obtained by a technique based on a constraint satisfaction algorithm that allows for the facile incorporation of secondary and other structural information. The program generates 3-D structures of RNA with atomic-level resolution that can be refined by numerical techniques such as energy minimization. The precision of this technique was evaluated by comparing predicted transfer RNA loop and RNA pseudoknot structures with known or consensus structures. The root-mean-square deviation (2.0 to 3.0 angstroms before minimization) between predicted and control structures reveal this system to be an effective method in modeling RNA.

Similarities of the cataleptic state induced by beta-endorphin and morphine

A variety of behavioral tests were used to characterize the cataleptic state induced by various treatments. Besides catalepsy, posture, locomotion, rigidity and the presence of reflexive responses were assessed. Measures of analgesia and body temperature were taken. The behavioral profiles of beta-endorphin, morphine, etonitazene, haloperidol, arecoline and GABA were compared at the time maximal catalepsy scores were obtained. Results indicated that, for an equivalent degree of catalepsy, the profile of beta-endorphin was similar to that of opiates, except for changes in body temperature; beta-endorphin's profile differed markedly from that of haloperidol, arecoline and GABA. Catalepsy was less pronounced with the latter two drugs. There were similarities in the behavioral profile of haloperidol and arecoline.