Fast and accurate genotype imputation in genome-wide association studies through pre-phasing

The 1000 Genomes Project and disease-specific sequencing efforts are producing large collections of haplotypes that can be used as reference panels for genotype imputation in genome-wide association studies (GWAS). However, imputing from large reference panels with existing methods imposes a high computational burden. We introduce a strategy called 'pre-phasing' that maintains the accuracy of leading methods while reducing computational costs. We first statistically estimate the haplotypes for each individual within the GWAS sample (pre-phasing) and then impute missing genotypes into these estimated haplotypes. This reduces the computational cost because (i) the GWAS samples must be phased only once, whereas standard methods would implicitly repeat phasing with each reference panel update, and (ii) it is much faster to match a phased GWAS haplotype to one reference haplotype than to match two unphased GWAS genotypes to a pair of reference haplotypes. We implemented our approach in the MaCH and IMPUTE2 frameworks, and we tested it on data sets from the Wellcome Trust Case Control Consortium 2 (WTCCC2), the Genetic Association Information Network (GAIN), the Women's Health Initiative (WHI) and the 1000 Genomes Project. This strategy will be particularly valuable for repeated imputation as reference panels evolve.     

Modulation of γδ T cell responses by TLR ligands

Toll-like receptors (TLR) are pattern-recognition receptors that recognize a broad variety of structurally conserved molecules derived from microbes. The recognition of TLR ligands functions as a primary sensor of the innate immune system, leading to subsequent indirect activation of the adaptive immunity as well as none-immune cells. However, TLR are also expressed by several T cell subsets, and the respective ligands can directly modulate their effector functions. The present review summarizes the recent findings of γδ T cell modulation by TLR ligands. TLR1/2/6, 3, and 5 ligands can act directly in combination with T cell receptor (TCR) stimulation to enhance cytokine/chemokine production of freshly isolated human γδ T cells. In contrast to human γδ T cells, murine and bovine γδ T cells can directly respond to TLR2 ligands with increased proliferation and cytokine production in a TCR-independent manner. Indirect stimulatory effects on IFN-γ production of human and murine γδ T cells via TLR-ligand activated dendritic cells have been described for TLR2, 3, 4, 7, and 9 ligands. In addition, TLR3 and 7 ligands indirectly increase tumor cell lysis by human γδ T cells, whereas ligation of TLR8 abolishes the suppressive activity of human tumor-infiltrating Vδ1 γδ T cells on αβ T cells and dendritic cells. Taken together, these data suggest that TLR-mediated signals received by γδ T cells enhance the initiation of adaptive immune responses during bacterial and viral infection directly or indirectly. Moreover, TLR ligands enhance cytotoxic tumor responses of γδ T cells and regulate the suppressive capacity of γδ T cells.     

Optimal temperatures for growth and upper thermal tolerance of juvenile northern leatherside chub

We examined optimal temperatures for growth and the upper thermal tolerance of juvenile northern leatherside chub ( Lepidomeda copei ). We conducted 2 experiments using the acclimated chronic-exposure method to estimate optimal temperature for growth of age-0 northern leatherside chub (range 12.8–28.3 °C). Upper thermal tolerance was estimated using the critical thermal maximum (CTM) and upper incipient lethal temperature (UILT) methods for fish acclimated at 15, 18, 23, and 28 °C. We also measured stream temperatures in Yellow Creek, Summit County, Utah, during July–August 2006 to compare our results to actual summer stream temperatures. Survival in growth tests was not significantly different between treatment temperatures in either experiment (P > 0.098). The optimal temperature for growth in the 1st trial estimated from the 2nd-order polynomial regression was 23.0 °C, falling outside the range of experimental temperatures (12.8–22.2 °C). The estimated optimal temperature for growth in the 2nd trial was 23.2 °C. In the upper thermal tolerance tests, juvenile northern leatherside chub had CTM values between 29.6 and 35.0 °C; CTM values increased as acclimation temperature increased. Upper incipient lethal temperatures (LT50) ranged from 26.5 to 30.2 °C, increasing with acclimation temperature. Summer stream temperatures in Yellow Creek had a lower mean (14.0–18.1 °C) than did the optimal temperature for growth determined in these studies, but these temperatures exhibited diel fluctuations as large as 15.7 °C.     

Arthropod and plant communities as indicators of land rehabilitation effectiveness in a semiarid shrubsteppe

We describe a case study evaluating the ecological impact of Bromus tectorum L. (cheatgrass) invasion following fire disturbance and the effectiveness of revegetation as a means of rehabilitation in a degraded semiarid shrubsteppe system. The effectiveness of rehabilitation efforts was assessed relative to arthropod richness, vegetation and arthropod community composition, and ground-cover characteristics in 3 habitats: undisturbed, burned and weed infested ( B. tectorum ), and burned then rehabilitated with native and nonnative vegetation. Arthropods were collected in each habitat using pitfall traps. Differences in arthropod richness were compared using rarefaction curves. Nonmetric multidimensional scaling and nonparametric multivariate statistical procedures, including analysis of similarity and similarity percentage routines, were used to compare arthropod and vegetation community composition and ground-cover characteristics between habitats. Arthropod communities in the rehabilitated habitat were distinct from those observed in the undisturbed and weed-infested habitats. Rehabilitation in this study resulted in a shift toward conditions observed in an undisturbed habitat and perhaps is an intermediate step to complete restoration. Arthropod richness, arthropod and vegetation community composition, and ground-cover characteristics were all useful indicators but returned slightly different results. Assessing multiple variables yielded the most complete understanding of the habitats studied.     

Plovers breeding in the highlands of Jalisco, Aguascalientes, Zacatecas, and San Luis Potosí, central Mexico

From April to July 2008, we surveyed for breeding plovers at 32 sites in the semiarid highlands of Jalisco, Aguascalientes, Zacatecas, and San Luis Potosí, in the Central Mexican High Plateau. We documented evidence or presumption of breeding Snowy Plovers ( Charadrius alexandrinus ) at 3 sites, Killdeer ( C. vociferus ) at 15 sites, and Mountain Plovers ( C. montanus ) at 1 site. Our surveys showed that the region is important breeding ground for only the Killdeer. We documented an apparent breeding range extension of the Mountain Plover to slightly more than 200 km south of its previously known breeding range.     

Changes in kit fox–coyote–prey relationships in the Great Basin Desert, Utah

Variation in kit fox ( Vulpes macrotis ) population parameters can be influenced by vegetative cover and the distribution and abundance of other predator and prey species. Dramatic changes to Great Basin Desert habitats, which can potentially impact mammalian species, have occurred in some areas in Utah. We examined kit fox demographics and prey populations from 1999 to 2001 on Dugway Proving Ground (DPG), a U.S. Army facility in Utah, and compared some parameters to historical levels (1956–1958, 1966–1969). Adult survival rates were fairly consistent between 1999 and 2000 and between 1999 and 2001; however, survival was greater in 2001 than in 2000. Reproductive rates ranged from 1.0 to 3.8 pups per female in 1999–2000 and were similar to historical numbers (1.0–4.2 pups per female). We found a decrease in pre-whelping kit fox density from the 1960s (0.12 foxes ? km –2 ) to 1999–2001 (0.04 foxes ? km –2 ); however, densities were similar between the current study and the 1950s (0.08 foxes ? km –2 ). Using 9 years of data, we found density dependence between reproductive rates of the current year and annual fox density from the previous year. Using 7 years of data, we found a slight correlation between kit fox annual density and a 1-year lag in leporid abundance, even though leporid abundance was lower during the present study than it was historically. Compared to historical levels, current small mammal abundance and species composition has changed in several habitats. Kit fox breeding density and annual density were inversely correlated with coyote ( Canis latrans ) density. Changes to the landscape at DPG, especially due to invasion of cheatgrass ( Bromus tectorum ) and addition of artificial water sources, have caused a change in available kit fox habitat and prey species, and have increased the abundance of coyotes, the kit fox's major competitor.     

Genome-wide association study reveals three susceptibility loci for common migraine in the general population

Migraine is a common, heterogeneous and heritable neurological disorder. Its pathophysiology is incompletely understood, and its genetic influences at the population level are unknown. In a population-based genome-wide analysis including 5,122 migraineurs and 18,108 non-migraineurs, rs2651899 (1p36.32, PRDM16), rs10166942 (2q37.1, TRPM8) and rs11172113 (12q13.3, LRP1) were among the top seven associations (P < 5 × 10(-6)) with migraine. These SNPs were significant in a meta-analysis among three replication cohorts and met genome-wide significance in a meta-analysis combining the discovery and replication cohorts (rs2651899, odds ratio (OR) = 1.11, P = 3.8 × 10(-9); rs10166942, OR = 0.85, P = 5.5 × 10(-12); and rs11172113, OR = 0.90, P = 4.3 × 10(-9)). The associations at rs2651899 and rs10166942 were specific for migraine compared with non-migraine headache. None of the three SNP associations was preferential for migraine with aura or without aura, nor were any associations specific for migraine features. TRPM8 has been the focus of neuropathic pain models, whereas LRP1 modulates neuronal glutamate signaling, plausibly linking both genes to migraine pathophysiology.     

Common variants at MS4A4/MS4A6E, CD2AP, CD33 and EPHA1 are associated with late-onset Alzheimer's disease

The Alzheimer Disease Genetics Consortium (ADGC) performed a genome-wide association study of late-onset Alzheimer disease using a three-stage design consisting of a discovery stage (stage 1) and two replication stages (stages 2 and 3). Both joint analysis and meta-analysis approaches were used. We obtained genome-wide significant results at MS4A4A (rs4938933; stages 1 and 2, meta-analysis P (P(M)) = 1.7 × 10(-9), joint analysis P (P(J)) = 1.7 × 10(-9); stages 1, 2 and 3, P(M) = 8.2 × 10(-12)), CD2AP (rs9349407; stages 1, 2 and 3, P(M) = 8.6 × 10(-9)), EPHA1 (rs11767557; stages 1, 2 and 3, P(M) = 6.0 × 10(-10)) and CD33 (rs3865444; stages 1, 2 and 3, P(M) = 1.6 × 10(-9)). We also replicated previous associations at CR1 (rs6701713; P(M) = 4.6 × 10(-10), P(J) = 5.2 × 10(-11)), CLU (rs1532278; P(M) = 8.3 × 10(-8), P(J) = 1.9 × 10(-8)), BIN1 (rs7561528; P(M) = 4.0 × 10(-14), P(J) = 5.2 × 10(-14)) and PICALM (rs561655; P(M) = 7.0 × 10(-11), P(J) = 1.0 × 10(-10)), but not at EXOC3L2, to late-onset Alzheimer's disease susceptibility.     

Truncating mutations in the last exon of NOTCH2 cause a rare skeletal disorder with osteoporosis

Hajdu-Cheney syndrome is a rare autosomal dominant skeletal disorder with facial anomalies, osteoporosis and acro-osteolysis. We sequenced the exomes of six unrelated individuals with this syndrome and identified heterozygous nonsense and frameshift mutations in NOTCH2 in five of them. All mutations cluster to the last coding exon of the gene, suggesting that the mutant mRNA products escape nonsense-mediated decay and that the resulting truncated NOTCH2 proteins act in a gain-of-function manner.