Joint Analysis for Genome-Wide Association Studies in Family-Based Designs

In family-based data, association information can be partitioned into the between-family information and the within-family information. Based on this observation, Steen et al. (Nature Genetics. 2005, 683–691) proposed an interesting two-stage test for genome-wide association (GWA) studies under family-based designs which performs genomic screening and replication using the same data set. In the first stage, a screening test based on the between-family information is used to select markers. In the second stage, an association test based on the within-family information is used to test association at the selected markers. However, we learn from the results of case-control studies (Skol et al. Nature Genetics. 2006, 209–213) that this two-stage approach may be not optimal. In this article, we propose a novel two-stage joint analysis for GWA studies under family-based designs. For this joint analysis, we first propose a new screening test that is based on the between-family information and is robust to population stratification. This new screening test is used in the first stage to select markers. Then, a joint test that combines the between-family information and within-family information is used in the second stage to test association at the selected markers. By extensive simulation studies, we demonstrate that the joint analysis always results in increased power to detect genetic association and is robust to population stratification

Ribbon Crystals

A repetitive crystal-like pattern is spontaneously formed upon the twisting of straight ribbons. The pattern is akin to a tessellation with isosceles triangles, and it can easily be demonstrated with ribbons cut from an overhead transparency. We give a general description of developable ribbons using a ruled procedure where ribbons are uniquely described by two generating functions. This construction defines a differentiable frame, the ribbon frame, which does not have singular points, whereby we avoid the shortcomings of the Frenet-Serret frame. The observed spontaneous pattern is modeled using planar triangles and cylindrical arcs, and the ribbon structure is shown to arise from a maximization of the end-to-end length of the ribbon, i.e. from an optimal use of ribbon length. The phenomenon is discussed in the perspectives of incompatible intrinsic geometries and of the emergence of long-range order

The Oxytocin Receptor (OXTR) Contributes to Prosocial Fund Allocations in the Dictator Game and the Social Value Orientations Task

Background: Economic games observe social decision making in the laboratory that involves real money payoffs. Previously we have shown that allocation of funds in the Dictator Game (DG), a paradigm that illustrates costly altruistic behavior, is partially determined by promoter-region repeat region variants in the arginine vasopressin 1a receptor gene (AVPR1a). In the current investigation, the gene encoding the related oxytocin receptor (OXTR) was tested for association with the DG and a related paradigm, the Social Values Orientation (SVO) task. Methodology/Principal Findings: Association (101 male and 102 female students) using a robust-family based test between 15 single tagging SNPs (htSNPs) across the OXTR was demonstrated with both the DG and SVO. Three htSNPs across the gene region showed significant association with both of the two games. The most significant association was observed with rs1042778 (p = 0.001). Haplotype analysis also showed significant associations for both DG and SVO. Following permutation test adjustment, significance was observed for 2–5 locus haplotypes (p,0.05). A second sample of 98 female subjects was subsequently and independently recruited to play the dictator game and was genotyped for the three significant SNPs found in the first sample. The rs1042778 SNP was shown to be significant for the second sample as well (p = 0.004, Fisher’s exact test). Conclusions: The demonstration that genetic polymorphisms for the OXTR are associated with human prosocial decisio

The Quantitative Methods Boot Camp:Teaching Quantitative Thinking and Computing Skills to Graduate Students in the Life Sciences

<div><p>The past decade has seen a rapid increase in the ability of biologists to collect large amounts of data. It is therefore vital that research biologists acquire the necessary skills during their training to visualize, analyze, and interpret such data. To begin to meet this need, we have developed a “boot camp” in quantitative methods for biology graduate students at Harvard Medical School. The goal of this short, intensive course is to enable students to use computational tools to visualize and analyze data, to strengthen their computational thinking skills, and to simulate and thus extend their intuition about the behavior of complex biological systems. The boot camp teaches basic programming using biological examples from statistics, image processing, and data analysis. This integrative approach to teaching programming and quantitative reasoning motivates students’ engagement by demonstrating the relevance of these skills to their work in life science laboratories. Students also have the opportunity to analyze their own data or explore a topic of interest in more detail. The class is taught with a mixture of short lectures, Socratic discussion, and in-class exercises. Students spend approximately 40% of their class time working through both short and long problems. A high instructor-to-student ratio allows students to get assistance or additional challenges when needed, thus enhancing the experience for students at all levels of mastery. Data collected from end-of-course surveys from the last five offerings of the course (between 2012 and 2014) show that students report high learning gains and feel that the course prepares them for solving quantitative and computational problems they will encounter in their research. We outline our course here which, together with the course materials freely available online under a Creative Commons License, should help to facilitate similar efforts by others.</p></div

Synaptic Dysbindin-1 Reductions in Schizophrenia Occur in an Isoform-Specific Manner Indicating Their Subsynaptic Location

Background: An increasing number of studies report associations between variation in DTNBP1, a top candidate gene in schizophrenia, and both the clinical symptoms of the disorder and its cognitive deficits. DTNBP1 encodes dysbindin-1, reduced levels of which have been found in synaptic fields of schizophrenia cases. This study determined whether such synaptic reductions are isoform-specific. Methodology/Principal Findings: Using Western blotting of tissue fractions, we first determined the synaptic localization of the three major dysbindin-1 isoforms (A, B, and C). All three were concentrated in synaptosomes of multiple brain areas, including auditory association cortices in the posterior half of the superior temporal gyrus (pSTG) and the hippocampal formation (HF). Tests on the subsynaptic tissue fractions revealed that each isoform is predominantly, if not exclusively, associated with synaptic vesicles (dysbindin-1B) or with postsynaptic densities (dysbindin-1A and -1C). Using Western blotting on pSTG (n = 15) and HF (n = 15) synaptosomal fractions from schizophrenia cases and their matched controls, we discovered that synaptic dysbindin-1 is reduced in an isoform-specific manner in schizophrenia without changes in levels of synaptophysin or PSD-95. In pSTG, about 92% of the schizophrenia cases displayed synaptic dysbindin-1A reductions averaging 48% (p = 0.0007) without alterations in other dysbindin-1 isoforms. In the HF, by contrast, schizophrenia cases displayed normal levels of synaptic dysbindin-1A, but 67% showed synaptic reductions in dysbindin-1B averaging 33% (p = 0.0256), while 80% showed synaptic reductions in dysbindin-1C averaging 35% (p = 0.0171). Conclusions/Significance: Given the distinctive subsynaptic localization of dysbindin-1A, -1B, and -1C across brain regions, the observed pSTG reductions in dysbindin-1A are postsynaptic and may promote dendritic spine loss with consequent disruption of auditory information processing, while the noted HF reductions in dysbindin-1B and -1C are both presynaptic and postsynaptic and could promote deficits in spatial working memory

Search for strongly interacting massive particles generating trackless jets in proton-proton collisions at s = 13 TeV

A search for dark matter in the form of strongly interacting massive particles (SIMPs) using the CMS detector at the LHC is presented. The SIMPs would be produced in pairs that manifest themselves as pairs of jets without tracks. The energy fraction of jets carried by charged particles is used as a key discriminator to suppress efficiently the large multijet background, and the remaining background is estimated directly from data. The search is performed using proton-proton collision data corresponding to an integrated luminosity of 16.1 fb - 1 , collected with the CMS detector in 2016. No significant excess of events is observed above the expected background. For the simplified dark matter model under consideration, SIMPs with masses up to 100 GeV are excluded and further sensitivity is explored towards higher masses

Search for dark matter produced in association with a Higgs boson decaying to a pair of bottom quarks in proton-proton collisions at root s=13TeV

A search for dark matter produced in association with a Higgs boson decaying to a pair of bottom quarks is performed in proton-proton collisions at a center-of-mass energy of 13 TeV collected with the CMS detector at the LHC. The analyzed data sample corresponds to an integrated luminosity of 35.9 fb(-1). The signal is characterized by a large missing transverse momentum recoiling against a bottom quark-antiquark system that has a large Lorentz boost. The number of events observed in the data is consistent with the standard model background prediction. Results are interpreted in terms of limits both on parameters of the type-2 two-Higgs doublet model extended by an additional light pseudoscalar boson a (2HDM+a) and on parameters of a baryonic Z simplified model. The 2HDM+a model is tested experimentally for the first time. For the baryonic Z model, the presented results constitute the most stringent constraints to date.Peer reviewe

Evidence for Top Quark Production in Nucleus-Nucleus Collisions

Peer reviewe

Observation of the B-s(0) -> X(3872)phi Decay

Using a data sample of proton-proton collisions at root s = 13 TeV, corresponding to an integrated luminosity of 140 fb(-1) collected by the CMS experiment in 2016-2018, the B-s(0) -> X(3872)phi decay is observed. Decays into J/psi pi(+)pi(-) and K+K- are used to reconstruct, respectively, the X(3872) and phi. The ratio of the product of branching fractions B[B-s(0) -> X(3872)phi]B[X(3872) -> J/psi pi(+)pi(-)] to the product B[B-s(0) ->psi(2S)phi]B[psi(2S) -> J/psi pi(+)pi(-)] is measured to be [2.21 +/- 0.29(stat) +/- 0.17(syst)]%. The ratio B[B-s(0) -> X(3872)phi]/B[B-0 -> X(3872)K-0] is found to be consistent with one, while the ratio B[B-s(0) -> X(3872)phi]/B[B+-> X(3872)K+] is two times smaller. This suggests a difference in the production dynamics of the X(3872) in B-0 and B(0)s meson decays compared to B+. The reported observation may shed new light on the nature of the X(3872) particle.Peer reviewe