CollapsABEL: An R library for detecting compound heterozygote alleles in genome-wide association studies

Background: Compound Heterozygosity (CH) in classical genetics is the presence of two different recessive mutations at a particular gene locus. A relaxed form of CH alleles may account for an essential proportion of the missing heritability, i.e. heritability of phenotypes so far not accounted for by single genetic variants. Methods to detect CH-like effects in genome-wide association studies (GWAS) may facilitate explaining the missing heritability, but to our knowledge no viable software tools for this purpose are currently available. Results: In this work we present the Generalized Compound Double Heterozygosity (GCDH) test and its implementation in the R package CollapsABEL. Time-consuming procedures are optimized for computational efficiency using Java or C++. Intermediate results are stored either in an SQL database or in a so-called big.matrix file to achieve reasonable memory footprint. Our large scale simulation studies show that GCDH is capable of discovering genetic associations due to CH-like interactions with much higher power than a conventional single-SNP approach under various settings, whether the causal genetic variations are available or not. CollapsABEL provides a user-friendly pipeline for genotype collapsing, statistical testing, power estimation, type I error control and graphics generation in the R language. Conclusions: CollapsABEL provides a computationally efficient solution for screening general forms of CH alleles in densely imputed microarray or whole genome sequencing datasets. The GCDH test provides an improved power over single-SNP based methods in detecting the prevalence of CH in human complex phenotypes, offering an opportunity for tackling the missing heritability problem. Binary and source packages of CollapsABEL are available on CRAN (https://cran.r-project.org/web/packages/CollapsABEL) and the website of the GenABEL project (http://www.genabel.org/packages)

Impurity Energy Level Within The Haldane Gap

An impurity bond $J{'}$ in a periodic 1D antiferromagnetic, spin 1 chain with exchange $J$ is considered. Using the numerical density matrix renormalization group method, we find an impurity energy level in the Haldane gap, corresponding to a bound state near the impurity bond. When $J{'}<J$ the level changes gradually from the edge of the Haldane gap to the ground state energy as the deviation $dev=(J-J{'})/J$ changes from 0 to 1. It seems that there is no threshold. Yet, there is a threshold when $J{'}>J$. The impurity level appears only when the deviation $dev=(J{'}-J)/J{'}$ is greater than $B_{c}$, which is near 0.3 in our calculation.Comment: Latex file,9 pages uuencoded compressed postscript including 4 figure

Genetics of skin color variation in Europeans: genome-wide association studies with functional follow-up

In the International Visible Trait Genetics (VisiGen) Consortium, we investigated the genetics of human skin color by combining a series of genome-wide association studies (GWAS) in a total of 17,262 Europeans with functional follow-up of discovered loci. Our GWAS provide the first genome-wide significant evidence for chromosome 20q11.22 harboring the ASIP gene being explicitly associated with skin color in Europeans. In addition, genomic loci at 5p13.2 (SLC45A2), 6p25.3 (IRF4), 15q13.1 (HERC2/OCA2), and 16q24.3 (MC1R) were confirmed to be involved in skin coloration in Europeans. In follow-up gene expression and regulation studies of 22 genes in 20q11.22, we highlighted two novel genes EIF2S2 and GSS, serving as competing functional candidates in this region and providing future research lines. A genetically inferred skin color score obtained from the 9 top-associated SNPs from 9 genes in 940 worldwide samples (HGDP-CEPH) showed a clear gradual pattern in Western Eurasians similar to the distribution of physical skin color, suggesting the used 9 SNPs as suitable markers for DNA prediction of skin color in Europeans and neighboring populations, relevant in future forensic and anthropological investigations

Single Spin Asymmetry ANA_N in Polarized Proton-Proton Elastic Scattering at s=200\sqrt{s}=200 GeV

We report a high precision measurement of the transverse single spin asymmetry $A_N$ at the center of mass energy $\sqrt{s}=200$ GeV in elastic proton-proton scattering by the STAR experiment at RHIC. The $A_N$ was measured in the four-momentum transfer squared $t$ range $0.003 \leqslant |t| \leqslant 0.035$ \GeVcSq, the region of a significant interference between the electromagnetic and hadronic scattering amplitudes. The measured values of $A_N$ and its $t$-dependence are consistent with a vanishing hadronic spin-flip amplitude, thus providing strong constraints on the ratio of the single spin-flip to the non-flip amplitudes. Since the hadronic amplitude is dominated by the Pomeron amplitude at this $\sqrt{s}$, we conclude that this measurement addresses the question about the presence of a hadronic spin flip due to the Pomeron exchange in polarized proton-proton elastic scattering.Comment: 12 pages, 6 figure

Longitudinal double-spin asymmetry and cross section for inclusive neutral pion production at midrapidity in polarized proton collisions at sqrt(s) = 200 GeV

We report a measurement of the longitudinal double-spin asymmetry A_LL and the differential cross section for inclusive Pi0 production at midrapidity in polarized proton collisions at sqrt(s) = 200 GeV. The cross section was measured over a transverse momentum range of 1 < p_T < 17 GeV/c and found to be in good agreement with a next-to-leading order perturbative QCD calculation. The longitudinal double-spin asymmetry was measured in the range of 3.7 < p_T < 11 GeV/c and excludes a maximal positive gluon polarization in the proton. The mean transverse momentum fraction of Pi0's in their parent jets was found to be around 0.7 for electromagnetically triggered events.Comment: 6 pages, 3 figures, submitted to Phys. Rev. D (RC

High pTp_{T} non-photonic electron production in pp+pp collisions at s\sqrt{s} = 200 GeV

We present the measurement of non-photonic electron production at high transverse momentum ($p_T >$ 2.5 GeV/$c$) in $p$ + $p$ collisions at $\sqrt{s}$ = 200 GeV using data recorded during 2005 and 2008 by the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). The measured cross-sections from the two runs are consistent with each other despite a large difference in photonic background levels due to different detector configurations. We compare the measured non-photonic electron cross-sections with previously published RHIC data and pQCD calculations. Using the relative contributions of B and D mesons to non-photonic electrons, we determine the integrated cross sections of electrons ($\frac{e^++e^-}{2}$) at 3 GeV/$c < p_T <~$10 GeV/$c$ from bottom and charm meson decays to be ${d\sigma_{(B\to e)+(B\to D \to e)} \over dy_e}|_{y_e=0}$ = 4.0$\pm0.5$({\rm stat.})$\pm1.1$({\rm syst.}) nb and ${d\sigma_{D\to e} \over dy_e}|_{y_e=0}$ = 6.2$\pm0.7$({\rm stat.})$\pm1.5$({\rm syst.}) nb, respectively.Comment: 17 pages, 17 figure

Evolution of the differential transverse momentum correlation function with centrality in Au+Au collisions at sNN=200\sqrt{s_{NN}} = 200 GeV

We present first measurements of the evolution of the differential transverse momentum correlation function, {\it C}, with collision centrality in Au+Au interactions at $\sqrt{s_{NN}} = 200$ GeV. {\it C} exhibits a strong dependence on collision centrality that is qualitatively similar to that of number correlations previously reported. We use the observed longitudinal broadening of the near-side peak of {\it C} with increasing centrality to estimate the ratio of the shear viscosity to entropy density, $\eta/s$, of the matter formed in central Au+Au interactions. We obtain an upper limit estimate of $\eta/s$ that suggests that the produced medium has a small viscosity per unit entropy.Comment: 7 pages, 4 figures, STAR paper published in Phys. Lett.

Longitudinal scaling property of the charge balance function in Au + Au collisions at 200 GeV

We present measurements of the charge balance function, from the charged particles, for diverse pseudorapidity and transverse momentum ranges in Au + Au collisions at 200 GeV using the STAR detector at RHIC. We observe that the balance function is boost-invariant within the pseudorapidity coverage [-1.3, 1.3]. The balance function properly scaled by the width of the observed pseudorapidity window does not depend on the position or size of the pseudorapidity window. This scaling property also holds for particles in different transverse momentum ranges. In addition, we find that the width of the balance function decreases monotonically with increasing transverse momentum for all centrality classes.Comment: 6 pages, 3 figure

Measurement of the Bottom contribution to non-photonic electron production in p+pp+p collisions at s\sqrt{s} =200 GeV

The contribution of $B$ meson decays to non-photonic electrons, which are mainly produced by the semi-leptonic decays of heavy flavor mesons, in $p+p$ collisions at $\sqrt{s} =$ 200 GeV has been measured using azimuthal correlations between non-photonic electrons and hadrons. The extracted $B$ decay contribution is approximately 50% at a transverse momentum of $p_{T} \geq 5$ GeV/$c$. These measurements constrain the nuclear modification factor for electrons from $B$ and $D$ meson decays. The result indicates that $B$ meson production in heavy ion collisions is also suppressed at high $p_{T}$.Comment: 6 pages, 4 figures, accepted by PR