Underestimated Effect Sizes in GWAS: Fundamental Limitations of Single SNP Analysis for Dichotomous Phenotypes

Complex diseases are often highly heritable. However, for many complex traits only a small proportion of the heritability can be explained by observed genetic variants in traditional genome-wide association (GWA) studies. Moreover, for some of those traits few significant SNPs have been identified. Single SNP association methods test for association at a single SNP, ignoring the effect of other SNPs. We show using a simple multi-locus odds model of complex disease that moderate to large effect sizes of causal variants may be estimated as relatively small effect sizes in single SNP association testing. This underestimation effect is most severe for diseases influenced by numerous risk variants. We relate the underestimation effect to the concept of non-collapsibility found in the statistics literature. As described, continuous phenotypes generated with linear genetic models are not affected by this underestimation effect. Since many GWA studies apply single SNP analysis to dichotomous phenotypes, previously reported results potentially underestimate true effect sizes, thereby impeding identification of true effect SNPs. Therefore, when a multi-locus model of disease risk is assumed, a multi SNP analysis may be more appropriate

Echinoderms have bilateral tendencies

Echinoderms take many forms of symmetry. Pentameral symmetry is the major form and the other forms are derived from it. However, the ancestors of echinoderms, which originated from Cambrian period, were believed to be bilaterians. Echinoderm larvae are bilateral during their early development. During embryonic development of starfish and sea urchins, the position and the developmental sequence of each arm are fixed, implying an auxological anterior/posterior axis. Starfish also possess the Hox gene cluster, which controls symmetrical development. Overall, echinoderms are thought to have a bilateral developmental mechanism and process. In this article, we focused on adult starfish behaviors to corroborate its bilateral tendency. We weighed their central disk and each arm to measure the position of the center of gravity. We then studied their turning-over behavior, crawling behavior and fleeing behavior statistically to obtain the center of frequency of each behavior. By joining the center of gravity and each center of frequency, we obtained three behavioral symmetric planes. These behavioral bilateral tendencies might be related to the A/P axis during the embryonic development of the starfish. It is very likely that the adult starfish is, to some extent, bilaterian because it displays some bilateral propensity and has a definite behavioral symmetric plane. The remainder of bilateral symmetry may have benefited echinoderms during their evolution from the Cambrian period to the present

Recurrent Modification of a Conserved Cis-Regulatory Element Underlies Fruit Fly Pigmentation Diversity

The development of morphological traits occurs through the collective action of networks of genes connected at the level of gene expression. As any node in a network may be a target of evolutionary change, the recurrent targeting of the same node would indicate that the path of evolution is biased for the relevant trait and network. Although examples of parallel evolution have implicated recurrent modification of the same gene and cis-regulatory element (CRE), little is known about the mutational and molecular paths of parallel CRE evolution. In Drosophila melanogaster fruit flies, the Bric-à-brac (Bab) transcription factors control the development of a suite of sexually dimorphic traits on the posterior abdomen. Female-specific Bab expression is regulated by the dimorphic element, a CRE that possesses direct inputs from body plan (ABD-B) and sex-determination (DSX) transcription factors. Here, we find that the recurrent evolutionary modification of this CRE underlies both intraspecific and interspecific variation in female pigmentation in the melanogaster species group. By reconstructing the sequence and regulatory activity of the ancestral Drosophila melanogaster dimorphic element, we demonstrate that a handful of mutations were sufficient to create independent CRE alleles with differing activities. Moreover, intraspecific and interspecific dimorphic element evolution proceeded with little to no alterations to the known body plan and sex-determination regulatory linkages. Collectively, our findings represent an example where the paths of evolution appear biased to a specific CRE, and drastic changes in function were accompanied by deep conservation of key regulatory linkages. © 2013 Rogers et al

Expression of Distal-less, dachshund, and optomotor blind in Neanthes arenaceodentata (Annelida, Nereididae) does not support homology of appendage-forming mechanisms across the Bilateria

The similarity in the genetic regulation of arthropod and vertebrate appendage formation has been interpreted as the product of a plesiomorphic gene network that was primitively involved in bilaterian appendage development and co-opted to build appendages (in modern phyla) that are not historically related as structures. Data from lophotrochozoans are needed to clarify the pervasiveness of plesiomorphic appendage forming mechanisms. We assayed the expression of three arthropod and vertebrate limb gene orthologs, Distal-less (Dll), dachshund (dac), and optomotor blind (omb), in direct-developing juveniles of the polychaete Neanthes arenaceodentata. Parapodial Dll expression marks premorphogenetic notopodia and neuropodia, becoming restricted to the bases of notopodial cirri and to ventral portions of neuropodia. In outgrowing cephalic appendages, Dll activity is primarily restricted to proximal domains. Dll expression is also prominent in the brain. dac expression occurs in the brain, nerve cord ganglia, a pair of pharyngeal ganglia, presumed interneurons linking a pair of segmental nerves, and in newly differentiating mesoderm. Domains of omb expression include the brain, nerve cord ganglia, one pair of anterior cirri, presumed precursors of dorsal musculature, and the same pharyngeal ganglia and presumed interneurons that express dac. Contrary to their roles in outgrowing arthropod and vertebrate appendages, Dll, dac, and omb lack comparable expression in Neanthes appendages, implying independent evolution of annelid appendage development. We infer that parapodia and arthropodia are not structurally or mechanistically homologous (but their primordia might be), that Dll’s ancestral bilaterian function was in sensory and central nervous system differentiation, and that locomotory appendages possibly evolved from sensory outgrowths

Simultaneous Analysis of All SNPs in Genome-Wide and Re-Sequencing Association Studies

Testing one SNP at a time does not fully realise the potential of genome-wide association studies to identify multiple causal variants, which is a plausible scenario for many complex diseases. We show that simultaneous analysis of the entire set of SNPs from a genome-wide study to identify the subset that best predicts disease outcome is now feasible, thanks to developments in stochastic search methods. We used a Bayesian-inspired penalised maximum likelihood approach in which every SNP can be considered for additive, dominant, and recessive contributions to disease risk. Posterior mode estimates were obtained for regression coefficients that were each assigned a prior with a sharp mode at zero. A non-zero coefficient estimate was interpreted as corresponding to a significant SNP. We investigated two prior distributions and show that the normal-exponential-gamma prior leads to improved SNP selection in comparison with single-SNP tests. We also derived an explicit approximation for type-I error that avoids the need to use permutation procedures. As well as genome-wide analyses, our method is well-suited to fine mapping with very dense SNP sets obtained from re-sequencing and/or imputation. It can accommodate quantitative as well as case-control phenotypes, covariate adjustment, and can be extended to search for interactions. Here, we demonstrate the power and empirical type-I error of our approach using simulated case-control data sets of up to 500 K SNPs, a real genome-wide data set of 300 K SNPs, and a sequence-based dataset, each of which can be analysed in a few hours on a desktop workstation

Predicting Unobserved Phenotypes for Complex Traits from Whole-Genome SNP Data

Genome-wide association studies (GWAS) for quantitative traits and disease in humans and other species have shown that there are many loci that contribute to the observed resemblance between relatives. GWAS to date have mostly focussed on discovery of genes or regulatory regions habouring causative polymorphisms, using single SNP analyses and setting stringent type-I error rates. Genome-wide marker data can also be used to predict genetic values and therefore predict phenotypes. Here, we propose a Bayesian method that utilises all marker data simultaneously to predict phenotypes. We apply the method to three traits: coat colour, %CD8 cells, and mean cell haemoglobin, measured in a heterogeneous stock mouse population. We find that a model that contains both additive and dominance effects, estimated from genome-wide marker data, is successful in predicting unobserved phenotypes and is significantly better than a prediction based upon the phenotypes of close relatives. Correlations between predicted and actual phenotypes were in the range of 0.4 to 0.9 when half of the number of families was used to estimate effects and the other half for prediction. Posterior probabilities of SNPs being associated with coat colour were high for regions that are known to contain loci for this trait. The prediction of phenotypes using large samples, high-density SNP data, and appropriate statistical methodology is feasible and can be applied in human medicine, forensics, or artificial selection programs

A search for the decay modes B+/- to h+/- tau l

We present a search for the lepton flavor violating decay modes B+/- to h+/- tau l (h= K,pi; l= e,mu) using the BaBar data sample, which corresponds to 472 million BBbar pairs. The search uses events where one B meson is fully reconstructed in one of several hadronic final states. Using the momenta of the reconstructed B, h, and l candidates, we are able to fully determine the tau four-momentum. The resulting tau candidate mass is our main discriminant against combinatorial background. We see no evidence for B+/- to h+/- tau l decays and set a 90% confidence level upper limit on each branching fraction at the level of a few times 10^-5.Comment: 15 pages, 7 figures, submitted to Phys. Rev.

Reactive oxygen species regulate context-dependent inhibition of NFAT5 target genes

The activation of nuclear factor of activated T cells 5(NFAT5), a well-known osmoprotective factor, can be induced by isotonic stimuli, such as activated Toll-like receptors (TLRs). It is unclear, however, how NFAT5 discriminates between isotonic and hypertonic stimuli. In this study we identified a novel context-dependent suppression of NFAT5 target gene expression in RAW 264.7 macrophages stimulated with lipopolysaccharide (LPS) or a high salt (NaCl) concentration. Although LPS and NaCl both used NFAT5 as a core transcription factor, these stimuli mutually inhibited distinct sets of NFAT5 targets within the cells. Although reactive oxygen species (ROS) are essential for this inhibition, the source of ROS differed depending on the context: mitochondria for high salt and xanthine oxidase for TLRs. Specifically, the high salt-induced suppression of interleukin-6 (IL-6) production was mediated through the ROS-induced inhibition of NFAT5 binding to the IL-6 promoter. The context-dependent inhibition of NFAT5 target gene expression was also confirmed in mouse spleen and kidney tissues that were cotreated with LPS and high salt. Taken together, our data suggest that ROS function as molecular sensors to discriminate between TLR ligation and osmotic stimuli in RAW 264.7 macrophages, directing NFAT5 activity toward proinflammatory or hypertonic responses in a context-dependent manner.open3

Observation and study of baryonic B decays: B -> D(*) p pbar, D(*) p pbar pi, and D(*) p pbar pi pi

We present a study of ten B-meson decays to a D(*), a proton-antiproton pair, and a system of up to two pions using BaBar's data set of 455x10^6 BBbar pairs. Four of the modes (B0bar -> D0 p anti-p, B0bar -> D*0 p anti-p, B0bar -> D+ p anti-p pi-, B0bar -> D*+ p anti-p pi-) are studied with improved statistics compared to previous measurements; six of the modes (B- -> D0 p anti-p pi-, B- -> D*0 p anti-p pi-, B0bar -> D0 p anti-p pi- pi+, B0bar -> D*0 p anti-p pi- pi+, B- -> D+ p anti-p pi- pi-, B- -> D*+ p anti-p pi- pi-) are first observations. The branching fractions for 3- and 5-body decays are suppressed compared to 4-body decays. Kinematic distributions for 3-body decays show non-overlapping threshold enhancements in m(p anti-p) and m(D(*)0 p) in the Dalitz plots. For 4-body decays, m(p pi-) mass projections show a narrow peak with mass and full width of (1497.4 +- 3.0 +- 0.9) MeV/c2, and (47 +- 12 +- 4) MeV/c2, respectively, where the first (second) errors are statistical (systematic). For 5-body decays, mass projections are similar to phase space expectations. All results are preliminary.Comment: 28 pages, 90 postscript figures, submitted to LP0