Exploring pleiotropy using principal components.

A standard multivariate principal components (PCs) method was utilized to identify clusters of variables that may be controlled by a common gene or genes (pleiotropy). Heritability estimates were obtained and linkage analyses performed on six individual traits (total cholesterol (Chol), high and low density lipoproteins, triglycerides (TG), body mass index (BMI), and systolic blood pressure (SBP)) and on each PC to compare our ability to identify major gene effects. Using the simulated data from Genetic Analysis Workshop 13 (Cohort 1 and 2 data for year 11), the quantitative traits were first adjusted for age, sex, and smoking (cigarettes per day). Adjusted variables were standardized and PCs calculated followed by orthogonal transformation (varimax rotation). Rotated PCs were then subjected to heritability and quantitative multipoint linkage analysis. The first three PCs explained 73% of the total phenotypic variance. Heritability estimates were above 0.60 for all three PCs. We performed linkage analyses on the PCs as well as the individual traits. The majority of pleiotropic and trait-specific genes were not identified. Standard PCs analysis methods did not facilitate the identification of pleiotropic genes affecting the six traits examined in the simulated data set. In addition, genes contributing 20% of the variance in traits with over 0.60 heritability estimates could not be identified in this simulated data set using traditional quantitative trait linkage analyses. Lack of identification of pleiotropic and trait-specific genes in some cases may reflect their low contribution to the traits/PCs examined or more importantly, characteristics of the sample group analyzed, and not simply a failure of the PC approach itself

Genomic differentiation during speciation-with-gene-flow: Comparing geographic and host-related variation in divergent life history adaptation in rhagoletis pomonella

A major goal of evolutionary biology is to understand how variation within populations gets partitioned into differences between reproductively isolated species. Here, we examine the degree to which diapause life history timing, a critical adaptation promoting population divergence, explains geographic and host-related genetic variation in ancestral hawthorn and recently derived apple-infesting races of Rhagoletis pomonella. Our strategy involved combining experiments on two different aspects of diapause (initial diapause intensity and adult eclosion time) with a geographic survey of genomic variation across four sites where apple and hawthorn flies co-occur from north to south in the Midwestern USA. The results demonstrated that the majority of the genome showing significant geographic and host-related variation can be accounted for by initial diapause intensity and eclosion time. Local genomic differences between sympatric apple and hawthorn flies were subsumed within broader geographic clines; allele frequency differences within the races across the Midwest were two to three-fold greater than those between the races in sympatry. As a result, sympatric apple and hawthorn populations displayed more limited genomic clustering compared to geographic populations within the races. The findings suggest that with reduced gene flow and increased selection on diapause equivalent to that seen between geographic sites, the host races may be recognized as different genotypic entities in sympatry, and perhaps species, a hypothesis requiring future genomic analysis of related sibling species to R. pomonella to test. Our findings concerning the way selection and geography interplay could be of broad significance for many cases of earlier stages of divergence-with-gene flow, including (1) where only modest increases in geographic isolation and the strength of selection may greatly impact genetic coupling and (2) the dynamics of how spatial and temporal standing variation is extracted by selection to generate differences between new and discrete units of biodiversity

Standing geographic variation in eclosion time and the genomics of host race formation in Rhagoletis pomonella fruit flies.

Taxa harboring high levels of standing variation may be more likely to adapt to rapid environmental shifts and experience ecological speciation. Here, we characterize geographic and host-related differentiation for 10,241 single nucleotide polymorphisms in Rhagoletis pomonella fruit flies to infer whether standing genetic variation in adult eclosion time in the ancestral hawthorn (Crataegus spp.)-infesting host race, as opposed to new mutations, contributed substantially to its recent shift to earlier fruiting apple (Malus domestica). Allele frequency differences associated with early vs. late eclosion time within each host race were significantly related to geographic genetic variation and host race differentiation across four sites, arrayed from north to south along a 430-km transect, where the host races co-occur in sympatry in the Midwest United States. Host fruiting phenology is clinal, with both apple and hawthorn trees fruiting earlier in the North and later in the South. Thus, we expected alleles associated with earlier eclosion to be at higher frequencies in northern populations. This pattern was observed in the hawthorn race across all four populations; however, allele frequency patterns in the apple race were more complex. Despite the generally earlier eclosion timing of apple flies and corresponding apple fruiting phenology, alleles on chromosomes 2 and 3 associated with earlier emergence were paradoxically at lower frequency in the apple than hawthorn host race across all four sympatric sites. However, loci on chromosome 1 did show higher frequencies of early eclosion-associated alleles in the apple than hawthorn host race at the two southern sites, potentially accounting for their earlier eclosion phenotype. Thus, although extensive clinal genetic variation in the ancestral hawthorn race exists and contributed to the host shift to apple, further study is needed to resolve details of how this standing variation was selected to generate earlier eclosing apple fly populations in the North

Search for the decay K+→π+ννˉK^+\to \pi^+ \nu \bar\nu in the momentum region Pπ<195 MeV/cP_\pi < 195 {\rm ~MeV/c}

We have searched for the decay $K^+ \to \pi^+ \nu \bar\nu$ in the kinematic region with pion momentum below the $K^+ \to \pi^+ \pi^0$ peak. One event was observed, consistent with the background estimate of $0.73\pm 0.18$. This implies an upper limit on $B(K^+ \to \pi^+ \nu \bar\nu)< 4.2\times 10^{-9}$ (90% C.L.), consistent with the recently measured branching ratio of $(1.57^{+1.75}_{-0.82}) \times 10^{-10}$, obtained using the standard model spectrum and the kinematic region above the $K^+ \to \pi^+ \pi^0$ peak. The same data were used to search for $K^+ \to \pi^+ X^0$, where $X^0$ is a weakly interacting neutral particle or system of particles with $150 < M_{X^0} < 250 {\rm ~MeV/c^2}$.Comment: 4 pages, 2 figure

First Observation of Coherent π0\pi^0 Production in Neutrino Nucleus Interactions with Eν<E_{\nu}< 2 GeV

The MiniBooNE experiment at Fermilab has amassed the largest sample to date of $\pi^0$s produced in neutral current (NC) neutrino-nucleus interactions at low energy. This paper reports a measurement of the momentum distribution of $\pi^0$s produced in mineral oil (CH$_2$) and the first observation of coherent $\pi^0$ production below 2 GeV. In the forward direction, the yield of events observed above the expectation for resonant production is attributed primarily to coherent production off carbon, but may also include a small contribution from diffractive production on hydrogen. Integrated over the MiniBooNE neutrino flux, the sum of the NC coherent and diffractive modes is found to be (19.5 $\pm$1.1 (stat) $\pm$2.5 (sys))% of all exclusive NC $\pi^0$ production at MiniBooNE. These measurements are of immediate utility because they quantify an important background to MiniBooNE's search for $\nu_{\mu} \to \nu_e$ oscillations.Comment: Submitted to Phys. Lett.

Further search for the decay K+→π+ννˉK^+ \to \pi^+ \nu \bar \nu in the momentum region P < 195 MeV/c

We report the results of a search for the decay $K^+ \to \pi^+ \nu \bar \nu$ in the kinematic region with $\pi^+$ momentum $140 < P < 195$ MeV/c using the data collected by the E787 experiment at BNL. No events were observed. When combined with our previous search in this region, one candidate event with an expected background of $1.22 \pm 0.24$ events results in a 90% C.L. upper limit of $2.2 \times 10^{-9}$ on the branching ratio of $K^+ \to \pi^+ \nu \bar \nu$. We also report improved limits on the rates of $K^+ \to \pi^+ X^0$ and $K^+ \to \pi^+ X^1 X^2$ where $X^0, X^1, X^2$ are hypothetical, massless, long-lived neutral particles.Comment: 5 pages, 3 figures, Accepted for publication in Phys. Rev.

The effects of rare SERPINA1 variants on lung function and emphysema in SPIROMICS

Rationale: The role of PI (protease inhibitor) type Z heterozygotes and additional rare variant genotypes in the gene encoding alpha-1 antitrypsin, SERPINA1 (serpin peptidase inhibitor, clade A,member 1), in determining chronic obstructive pulmonary disease risk and severity is controversial. Objectives: To comprehensively evaluate the effects of rare SERPINA1 variants on lung function and emphysema phenotypes in subjects with significant tobacco smoke exposure using deep gene resequencing and alpha-1 antitrypsin concentrations. Methods: DNA samples from 1,693 non-Hispanic white individuals, 385 African Americans, and 90 Hispanics with >20 pack-years smoking were resequenced for the identification of rare variants (allele frequency,0.05) in 16.9 kB of SERPINA1. Measurements and Main Results: White PI Z heterozygotes confirmed by sequencing (MZ; n = 74) had lower postbronchodilator FEV1 (P = 0.007), FEV1/FVC (P = 0.003), and greater computed tomography-based emphysema (P = 0.02) compared with 1,411 white individuals without PI Z, S, or additional rare variants denoted as VR. PI Z-containing compound heterozygotes (ZS/ZVR; n = 7) had lower FEV1/FVC (P = 0.02) and forced expiratory flow, midexpiratory phase (P = 0.009). Nineteen white heterozygotes for five non-S/Z coding variants associated with lower alpha-1 antitrypsin had greater computed tomography-based emphysema compared with those without rare variants. In African Americans, a 59 untranslated region insertion (rs568223361) was associated with lower alpha-1 antitrypsin and functional small airway disease (P = 0.007). Conclusions: In this integrative deep sequencing study of SERPINA1 with alpha-1 antitrypsin concentrations in a heavy smoker and chronic obstructive pulmonary disease cohort, we confirmed the effects of PI Z heterozygote and compound heterozygote genotypes. We demonstrate the cumulative effects of multiple SERPINA1 variants on alpha-1 antitrypsin deficiency, lung function, and emphysema, thus significantly increasing the frequency of SERPINA1 variation associated with respiratory disease in at-risk smokers

Overlap of Genetic Risk between Interstitial Lung Abnormalities and Idiopathic Pulmonary Fibrosis

Rationale: Interstitial lung abnormalities (ILAs) are associated with the highest genetic risk locus for idiopathic pulmonary fibrosis (IPF); however, the extent to which there are unique associations among individuals with ILAs or additional overlap with IPF is not known.Objectives: To perform a genome-wide association study (GWAS) of ILAs.Methods: ILAs and a subpleural-predominant subtype were assessed on chest computed tomography (CT) scans in the AGES (Age Gene/Environment Susceptibility), COPDGene (Genetic Epidemiology of Chronic Obstructive Pulmonary Disease [COPD]), Framingham Heart, ECLIPSE (Evaluation of COPD Longitudinally to Identify Predictive Surrogate End-points), MESA (Multi-Ethnic Study of Atherosclerosis), and SPIROMICS (Subpopulations and Intermediate Outcome Measures in COPD Study) studies. We performed a GWAS of ILAs in each cohort and combined the results using a meta-analysis. We assessed for overlapping associations in independent GWASs of IPF.Measurements and Main Results: Genome-wide genotyping data were available for 1,699 individuals with ILAs and 10,274 control subjects. The MUC5B (mucin 5B) promoter variant rs35705950 was significantly associated with both ILAs (P = 2.6 × 10-27) and subpleural ILAs (P = 1.6 × 10-29). We discovered novel genome-wide associations near IPO11 (rs6886640, P = 3.8 × 10-8) and FCF1P3 (rs73199442, P = 4.8 × 10-8) with ILAs, and near HTRE1 (rs7744971, P = 4.2 × 10-8) with subpleural-predominant ILAs. These novel associations were not associated with IPF. Among 12 previously reported IPF GWAS loci, five (DPP9, DSP, FAM13A, IVD, and MUC5B) were significantly associated (P < 0.05/12) with ILAs.Conclusions: In a GWAS of ILAs in six studies, we confirmed the association with a MUC5B promoter variant and found strong evidence for an effect of previously described IPF loci; however, novel ILA associations were not associated with IPF. These findings highlight common genetically driven biologic pathways between ILAs and IPF, and also suggest distinct ones