Using genetic variants to evaluate the causal effect of cholesterol lowering on head and neck cancer risk: A Mendelian randomization study

Head and neck squamous cell carcinoma (HNSCC), which includes cancers of the oral cavity and oropharynx, is a cause of substantial global morbidity and mortality. Strategies to reduce disease burden include discovery of novel therapies and repurposing of existing drugs. Statins are commonly prescribed for lowering circulating cholesterol by inhibiting HMG-CoA reductase (HMGCR). Results from some observational studies suggest that statin use may reduce HNSCC risk. We appraised the relationship of genetically-proxied cholesterol-lowering drug targets and other circulating lipid traits with oral (OC) and oropharyngeal (OPC) cancer risk using two-sample Mendelian randomization (MR). For the primary analysis, germline genetic variants in HMGCR, NPC1L1, CETP, PCSK9 and LDLR were used to proxy the effect of low-density lipoprotein cholesterol (LDL-C) lowering therapies. In secondary analyses, variants were used to proxy circulating levels of other lipid traits in a genome-wide association study (GWAS) meta-analysis of 188,578 individuals. Both primary and secondary analyses aimed to estimate the downstream causal effect of cholesterol lowering therapies on OC and OPC risk. The second sample for MR was taken from a GWAS of 6,034 OC and OPC cases and 6,585 controls (GAME-ON). Analyses were replicated in UK Biobank, using 839 OC and OPC cases and 372,016 controls and the results of the GAME-ON and UK Biobank analyses combined in a fixed-effects meta-analysis. We found limited evidence of a causal effect of genetically-proxied LDL-C lowering using HMGCR, NPC1L1, CETP or other circulating lipid traits on either OC or OPC risk. Genetically-proxied PCSK9 inhibition equivalent to a 1 mmol/L (38.7 mg/dL) reduction in LDL-C was associated with an increased risk of OC and OPC combined (OR 1.8 95%CI 1.2, 2.8, p = 9.31 x10^-05), with good concordance between GAME-ON and UK Biobank (I² = 22%). Effects for PCSK9 appeared stronger in relation to OPC (OR 2.6 95%CI 1.4, 4.9) than OC (OR 1.4 95%CI 0.8, 2.4). LDLR variants, resulting in genetically-proxied reduction in LDL-C equivalent to a 1 mmol/L (38.7 mg/dL), reduced the risk of OC and OPC combined (OR 0.7, 95%CI 0.5, 1.0, p = 0.006). A series of pleiotropy-robust and outlier detection methods showed that pleiotropy did not bias our findings. We found limited evidence for a role of cholesterol-lowering in OC and OPC risk, suggesting previous observational results may have been confounded. There was some evidence that genetically-proxied inhibition of PCSK9 increased risk, while lipid-lowering variants in LDLR, reduced risk of combined OC and OPC. This result suggests that the mechanisms of action of PCSK9 on OC and OPC risk may be independent of its cholesterol lowering effects; however, this was not supported uniformly across all sensitivity analyses and further replication of this finding is required.     

Exploiting collider bias to apply two-sample summary data Mendelian randomization methods to one-sample individual level data

Over the last decade the availability of SNP-trait associations from genome-wide association studies has led to an array of methods for performing Mendelian randomization studies using only summary statistics. A common feature of these methods, besides their intuitive simplicity, is the ability to combine data from several sources, incorporate multiple variants and account for biases due to weak instruments and pleiotropy. With the advent of large and accessible fully-genotyped cohorts such as UK Biobank, there is now increasing interest in understanding how best to apply these well developed summary data methods to individual level data, and to explore the use of more sophisticated causal methods allowing for non-linearity and effect modification.In this paper we describe a general procedure for optimally applying any two sample summary data method using one sample data. Our procedure first performs a meta-analysis of summary data estimates that are intentionally contaminated by collider bias between the genetic instruments and unmeasured confounders, due to conditioning on the observed exposure. These estimates are then used to correct the standard observational association between an exposure and outcome. Simulations are conducted to demonstrate the method’s performance against naive applications of two sample summary data MR. We apply the approach to the UK Biobank cohort to investigate the causal role of sleep disturbance on HbA1c levels, an important determinant of diabetes.Our approach can be viewed as a generalization of Dudbridge et al. (Nat. Comm. 10: 1561), who developed a technique to adjust for index event bias when uncovering genetic predictors of disease progression based on case-only data. Our work serves to clarify that in any one sample MR analysis, it can be advantageous to estimate causal relationships by artificially inducing and then correcting for collider bias.     

Causal association between heart failure and bone mineral density: Insights from a two-sample bidirectional Mendelian randomization study

In recent times, the association between HF and BMD has attracted enormous interest in the scientific community. However, published epidemiological observational studies on the relationship between HF and BMD remain inconclusive. Herein, we evaluated from the analytical perspective a two-sample bidirectional MR study to analyze the causal association between HF and BMD using a summary-level GWAS Catalog. To select instrumental SNPs strongly associated with exposure, we took a series of rigorous quality control steps at the time of analysis. The causal MR assessment of HF on the risk of BMD was performed first and then in the opposite direction. To make the conclusions more reliable and robust, the fixed-effects IVW, weighted median-based method, MR–Egger, simple mode and weighted mode were utilized. A maximum likelihood model was also used if necessary. MR–Egger regression, IVW “leave-one-out” sensitivity analysis, MR-PRESSO, MR–Egger intercept test and Cochran's Q statistic methods were used to assess heterogeneity and pleiotropy. Our MR studies supported a meaningful causal association between HF and TB-BMD (IVW: OR = 0.78, 95% CI: 0.68–0.87, p = 0.00588). At the same time, we did not find a significant causal relationship between HF and FA-BMD, FN-BMD or LS-BMD. No significant causal relationships between BMD and HF were observed. This bidirectional MR analysis suggested a causal association of HF with only low TB-BMD, while the reverse causality hypothesis was not found. Studies of the prevention and treatment of total bone mineral density decline in patients with heart failure need to be performed.     

Vitamin D levels and risk of type 1 diabetes: A Mendelian randomization study

BackgroundVitamin D deficiency has been associated with type 1 diabetes in observational studies, but evidence from randomized controlled trials (RCTs) is lacking. The aim of this study was to test whether genetically decreased vitamin D levels are causally associated with type 1 diabetes using Mendelian randomization (MR).Methods and findingsFor our two-sample MR study, we selected as instruments single nucleotide polymorphisms (SNPs) that are strongly associated with 25-hydroxyvitamin D (25OHD) levels in a large vitamin D genome-wide association study (GWAS) on 443,734 Europeans and obtained their corresponding effect estimates on type 1 diabetes risk from a large meta-analysis of 12 type 1 diabetes GWAS studies (Ntot = 24,063, 9,358 cases, and 15,705 controls). In addition to the main analysis using inverse variance weighted MR, we applied 3 additional methods to control for pleiotropy (MR-Egger, weighted median, and mode-based estimate) and compared the respective MR estimates. We also undertook sensitivity analyses excluding SNPs with potential pleiotropic effects. We identified 69 lead independent common SNPs to be genome-wide significant for 25OHD, explaining 3.1% of the variance in 25OHD levels. MR analyses suggested that a 1 standard deviation (SD) decrease in standardized natural log-transformed 25OHD (corresponding to a 29-nmol/l change in 25OHD levels in vitamin D–insufficient individuals) was not associated with an increase in type 1 diabetes risk (inverse-variance weighted (IVW) MR odds ratio (OR) = 1.09, 95% CI: 0.86 to 1.40, p = 0.48). We obtained similar results using the 3 pleiotropy robust MR methods and in sensitivity analyses excluding SNPs associated with serum lipid levels, body composition, blood traits, and type 2 diabetes. Our findings indicate that decreased vitamin D levels did not have a substantial impact on risk of type 1 diabetes in the populations studied. Study limitations include an inability to exclude the existence of smaller associations and a lack of evidence from non-European populations.ConclusionsOur findings suggest that 25OHD levels are unlikely to have a large effect on risk of type 1 diabetes, but larger MR studies or RCTs are needed to investigate small effects.

Despoina Manousaki and co-workers investigate vitamin D levels and risk of type I diabetes.     

MRLocus: Identifying causal genes mediating a trait through Bayesian estimation of allelic heterogeneity

Expression quantitative trait loci (eQTL) studies are used to understand the regulatory function of non-coding genome-wide association study (GWAS) risk loci, but colocalization alone does not demonstrate a causal relationship of gene expression affecting a trait. Evidence for mediation, that perturbation of gene expression in a given tissue or developmental context will induce a change in the downstream GWAS trait, can be provided by two-sample Mendelian Randomization (MR). Here, we introduce a new statistical method, MRLocus, for Bayesian estimation of the gene-to-trait effect from eQTL and GWAS summary data for loci with evidence of allelic heterogeneity, that is, containing multiple causal variants. MRLocus makes use of a colocalization step applied to each nearly-LD-independent eQTL, followed by an MR analysis step across eQTLs. Additionally, our method involves estimation of the extent of allelic heterogeneity through a dispersion parameter, indicating variable mediation effects from each individual eQTL on the downstream trait. Our method is evaluated against other state-of-the-art methods for estimation of the gene-to-trait mediation effect, using an existing simulation framework. In simulation, MRLocus often has the highest accuracy among competing methods, and in each case provides more accurate estimation of uncertainty as assessed through interval coverage. MRLocus is then applied to five candidate causal genes for mediation of particular GWAS traits, where gene-to-trait effects are concordant with those previously reported. We find that MRLocus’s estimation of the causal effect across eQTLs within a locus provides useful information for determining how perturbation of gene expression or individual regulatory elements will affect downstream traits. The MRLocus method is implemented as an R package available at https://mikelove.github.io/mrlocus.     

Cu2+ Inhibits Photosystem II Activities but Enhances Photosystem I Quantum Yield of Microcystis aeruginosa

Responses of photosystem I and II activities of Microcystis aeruginosa to various concentrations of Cu²⁺ were simultaneously examined using a Dual-PAM-100 fluorometer. Cell growth and contents of chlorophyll a were significantly inhibited by Cu²⁺. Photosystem II activity [Y(II)] and electron transport [rETR_max(II)] were significantly altered by Cu²⁺. The quantum yield of photosystem II [Y(II)] decreased by 29 % at 100 μg L^?1 Cu²⁺ compared to control. On the contrary, photosystem I was stable under Cu²⁺ stress and showed an obvious increase of quantum yield [Y(I)] and electron transport [rETR_max(I)] due to activation of cyclic electron flow (CEF). Yield of cyclic electron flow [Y(CEF)] was enhanced by 17 % at 100 μg L^?1 Cu²⁺ compared to control. The contribution of linear electron flow to photosystem I [Y(II)/Y(I)] decreased with increasing Cu²⁺ concentration. Yield of cyclic electron flow [Y(CEF)] was negatively correlated with the maximal photosystem II photochemical efficiency (F _v/F _m). In summary, photosystem II was the major target sites of toxicity of Cu²⁺, while photosystem I activity was enhanced under Cu²⁺ stress.     

A Balanced Tissue Composition Reveals New Metabolic and Gene Expression Markers in Prostate Cancer

Molecular analysis of patient tissue samples is essential to characterize the in vivo variability in human cancers which are not accessible in cell-lines or animal models. This applies particularly to studies of tumor metabolism. The challenge is, however, the complex mixture of various tissue types within each sample, such as benign epithelium, stroma and cancer tissue, which can introduce systematic biases when cancers are compared to normal samples. In this study we apply a simple strategy to remove such biases using sample selections where the average content of stroma tissue is balanced between the sample groups. The strategy is applied to a prostate cancer patient cohort where data from MR spectroscopy and gene expression have been collected from and integrated on the exact same tissue samples. We reveal in vivo changes in cancer-relevant metabolic pathways which are otherwise hidden in the data due to tissue confounding. In particular, lowered levels of putrescine are connected to increased expression of SRM, reduced levels of citrate are attributed to upregulation of genes promoting fatty acid synthesis, and increased succinate levels coincide with reduced expression of SUCLA2 and SDHD. In addition, the strategy also highlights important metabolic differences between the stroma, epithelium and prostate cancer. These results show that important in vivo metabolic features of cancer can be revealed from patient data only if the heterogeneous tissue composition is properly accounted for in the analysis.     

A practical problem with Egger regression in Mendelian randomization

Mendelian randomization (MR) is an instrumental variable (IV) method using genetic variants such as single nucleotide polymorphisms (SNPs) as IVs to disentangle the causal relationship between an exposure and an outcome. Since any causal conclusion critically depends on the three valid IV assumptions, which will likely be violated in practice, MR methods robust to the IV assumptions are greatly needed. As such a method, Egger regression stands out as one of the most widely used due to its easy use and perceived robustness. Although Egger regression is claimed to be robust to directional pleiotropy under the instrument strength independent of direct effect (InSIDE) assumption, it is known to be dependent on the orientations/coding schemes of SNPs (i.e. which allele of an SNP is selected as the reference group). The current practice, as recommended as the default setting in some popular MR software packages, is to orientate the SNPs to be all positively associated with the exposure, which however, to our knowledge, has not been fully studied to assess its robustness and potential impact. We use both numerical examples (with both real data and simulated data) and analytical results to demonstrate the practical problem of Egger regression with respect to its heavy dependence on the SNP orientations. Under the assumption that InSIDE holds for some specific (and unknown) coding scheme of the SNPs, we analytically show that other coding schemes would in general lead to the violation of InSIDE. Other related MR and IV regression methods may suffer from the same problem. Cautions should be taken when applying Egger regression (and related MR and IV regression methods) in practice.     

Identifying drug targets for neurological and psychiatric disease via genetics and the brain transcriptome

Discovering drugs that efficiently treat brain diseases has been challenging. Genetic variants that modulate the expression of potential drug targets can be utilized to assess the efficacy of therapeutic interventions. We therefore employed Mendelian Randomization (MR) on gene expression measured in brain tissue to identify drug targets involved in neurological and psychiatric diseases. We conducted a two-sample MR using cis-acting brain-derived expression quantitative trait loci (eQTLs) from the Accelerating Medicines Partnership for Alzheimer’s Disease consortium (AMP-AD) and the CommonMind Consortium (CMC) meta-analysis study (n = 1,286) as genetic instruments to predict the effects of 7,137 genes on 12 neurological and psychiatric disorders. We conducted Bayesian colocalization analysis on the top MR findings (using P<6x10^-7 as evidence threshold, Bonferroni-corrected for 80,557 MR tests) to confirm sharing of the same causal variants between gene expression and trait in each genomic region. We then intersected the colocalized genes with known monogenic disease genes recorded in Online Mendelian Inheritance in Man (OMIM) and with genes annotated as drug targets in the Open Targets platform to identify promising drug targets. 80 eQTLs showed MR evidence of a causal effect, from which we prioritised 47 genes based on colocalization with the trait. We causally linked the expression of 23 genes with schizophrenia and a single gene each with anorexia, bipolar disorder and major depressive disorder within the psychiatric diseases and 9 genes with Alzheimer’s disease, 6 genes with Parkinson’s disease, 4 genes with multiple sclerosis and two genes with amyotrophic lateral sclerosis within the neurological diseases we tested. From these we identified five genes (ACE, GPNMB, KCNQ5, RERE and SUOX) as attractive drug targets that may warrant follow-up in functional studies and clinical trials, demonstrating the value of this study design for discovering drug targets in neuropsychiatric diseases.     

Effect of yttrium on calcium-dependent processes in vertebrate myocardium

Inoptopic effect of yttrium acetate (Y³⁺) on myocardium of the marsh frog Rana ridibunda and its effect on ion transport across the inner mitochondrial membrane (IMM) of rat heart was studied. Y³⁺ was found to decrease the rate of heart contractions and to stimulate ion transport in the rat heart mitochondria in media with 10 mM glutamate and 2 mM malate. Presence of Y³⁺ induced inhibition of energy-dependent Ca²⁺ transport into mitochondria, which was expressed as a marked decrease of their swelling in the media containing 125 mM NH₄NO₃ and Ca²⁺ or 25 mM potassium acetate, 100 mM sucrose and Ca²⁺. It is suggested that the Y³⁺-induced decrease in rat muscle contractions is determined not only by direct suppressing effect of Y³⁺ on potential-modulated Ca²⁺-channels of pacemaker and contractile cardiomyocytes (CM), but also by its indirect effect on Ca²⁺-carrier in IMM. The data confirming that Y³⁺ activates energy-dependent K⁺ transport catalyzed by mitochondrial uniporter and blocks Ca²⁺-channels in the mitochondrial membrane are important for more complete understanding of mechanisms of the Y³⁺ action on vertebrates and human CM.     

Causal inference for heritable phenotypic risk factors using heterogeneous genetic instruments

Over a decade of genome-wide association studies (GWAS) have led to the finding of extreme polygenicity of complex traits. The phenomenon that “all genes affect every complex trait” complicates Mendelian Randomization (MR) studies, where natural genetic variations are used as instruments to infer the causal effect of heritable risk factors. We reexamine the assumptions of existing MR methods and show how they need to be clarified to allow for pervasive horizontal pleiotropy and heterogeneous effect sizes. We propose a comprehensive framework GRAPPLE to analyze the causal effect of target risk factors with heterogeneous genetic instruments and identify possible pleiotropic patterns from data. By using GWAS summary statistics, GRAPPLE can efficiently use both strong and weak genetic instruments, detect the existence of multiple pleiotropic pathways, determine the causal direction and perform multivariable MR to adjust for confounding risk factors. With GRAPPLE, we analyze the effect of blood lipids, body mass index, and systolic blood pressure on 25 disease outcomes, gaining new information on their causal relationships and potential pleiotropic pathways involved.     

Robust inference of bi-directional causal relationships in presence of correlated pleiotropy with GWAS summary data

To infer a causal relationship between two traits, several correlation-based causal direction (CD) methods have been proposed with the use of SNPs as instrumental variables (IVs) based on GWAS summary data for the two traits; however, none of the existing CD methods can deal with SNPs with correlated pleiotropy. Alternatively, reciprocal Mendelian randomization (MR) can be applied, which however may perform poorly in the presence of (unknown) invalid IVs, especially for bi-directional causal relationships. In this paper, first, we propose a CD method that performs better than existing CD methods regardless of the presence of correlated pleiotropy. Second, along with a simple but yet effective IV screening rule, we propose applying a closely related and state-of-the-art MR method in reciprocal MR, showing its almost identical performance to that of the new CD method when their model assumptions hold; however, if the modeling assumptions are violated, the new CD method is expected to better control type I errors. Notably bi-directional causal relationships impose some unique challenges beyond those for uni-directional ones, and thus requiring special treatments. For example, we point out for the first time several scenarios where a bi-directional relationship, but not a uni-directional one, can unexpectedly cause the violation of some weak modeling assumptions commonly required by many robust MR methods. We also offer some numerical support and a modeling justification for the application of our new methods (and more generally MR) to binary traits. Finally we applied the proposed methods to 12 risk factors and 4 common diseases, confirming mostly well-known uni-directional causal relationships, while identifying some novel and plausible bi-directional ones such as between body mass index and type 2 diabetes (T2D), and between diastolic blood pressure and stroke.     

Combining the strengths of inverse-variance weighting and Egger regression in Mendelian randomization using a mixture of regressions model

With the increasing availability of large-scale GWAS summary data on various traits, Mendelian randomization (MR) has become commonly used to infer causality between a pair of traits, an exposure and an outcome. It depends on using genetic variants, typically SNPs, as instrumental variables (IVs). The inverse-variance weighted (IVW) method (with a fixed-effect meta-analysis model) is most powerful when all IVs are valid; however, when horizontal pleiotropy is present, it may lead to biased inference. On the other hand, Egger regression is one of the most widely used methods robust to (uncorrelated) pleiotropy, but it suffers from loss of power. We propose a two-component mixture of regressions to combine and thus take advantage of both IVW and Egger regression; it is often both more efficient (i.e. higher powered) and more robust to pleiotropy (i.e. controlling type I error) than either IVW or Egger regression alone by accounting for both valid and invalid IVs respectively. We propose a model averaging approach and a novel data perturbation scheme to account for uncertainties in model/IV selection, leading to more robust statistical inference for finite samples. Through extensive simulations and applications to the GWAS summary data of 48 risk factor-disease pairs and 63 genetically uncorrelated trait pairs, we showcase that our proposed methods could often control type I error better while achieving much higher power than IVW and Egger regression (and sometimes than several other new/popular MR methods). We expect that our proposed methods will be a useful addition to the toolbox of Mendelian randomization for causal inference.     

The Genetic and Cytological Organization of the Y Chromosome of DROSOPHILA MELANOGASTER

Cytological and genetic analyses of 121 translocations between the Y chromosome and the centric heterochromatin of the X chromosome have been used to define and localize six regions on the Y chromosome of Drosophila melanogaster necessary for male fertility. These regions are associated with nonfluorescent blocks of the Y chromosome, as revealed using Hoechst 33258 or quinacrine staining. Each region appears to contain but one functional unit, as defined by failure of complementation among translocations with breakpoints within the same block. The distribution of translocation breakpoints examined appears to be nonrandom, in that breaks occur preferentially in the nonfluorescent blocks and not in the large fluorescent blocks.     

Genome-Wide Prediction of Traits with Different Genetic Architecture Through Efficient Variable Selection

In genome-based prediction there is considerable uncertainty about the statistical model and method required to maximize prediction accuracy. For traits influenced by a small number of quantitative trait loci (QTL), predictions are expected to benefit from methods performing variable selection [e.g., BayesB or the least absolute shrinkage and selection operator (LASSO)] compared to methods distributing effects across the genome [ridge regression best linear unbiased prediction (RR-BLUP)]. We investigate the assumptions underlying successful variable selection by combining computer simulations with large-scale experimental data sets from rice (Oryza sativa L.), wheat (Triticum aestivum L.), and Arabidopsis thaliana (L.). We demonstrate that variable selection can be successful when the number of phenotyped individuals is much larger than the number of causal mutations contributing to the trait. We show that the sample size required for efficient variable selection increases dramatically with decreasing trait heritabilities and increasing extent of linkage disequilibrium (LD). We contrast and discuss contradictory results from simulation and experimental studies with respect to superiority of variable selection methods over RR-BLUP. Our results demonstrate that due to long-range LD, medium heritabilities, and small sample sizes, superiority of variable selection methods cannot be expected in plant breeding populations even for traits like FRIGIDA gene expression in Arabidopsis and flowering time in rice, assumed to be influenced by a few major QTL. We extend our conclusions to the analysis of whole-genome sequence data and infer upper bounds for the number of causal mutations which can be identified by LASSO. Our results have major impact on the choice of statistical method needed to make credible inferences about genetic architecture and prediction accuracy of complex traits.     

Approximation of bias and mean-squared error in two-sample Mendelian randomization analyses

Mendelian randomization (MR) is a type of instrumental variable (IV) analysis that uses genetic variants as IVs for a risk factor to study its causal effect on an outcome. Extensive investigations on the performance of IV analysis procedures, such as the one based on the two-stage least squares (2SLS) procedure, have been conducted under the one-sample scenario, where measures on IVs, the risk factor, and the outcome are assumed to be available for each study participant. Recent MR analysis usually is performed with data from two independent or partially overlapping genetic association studies (two-sample setting), with one providing information on the association between the IVs and the outcome, and the other on the association between the IVs and the risk factor. We investigate the performance of 2SLS in the two-sample–based MR when the IVs are weakly associated with the risk factor. We derive closed form formulas for the bias and mean squared error of the 2SLS estimate and verify them with numeric simulations under realistic circumstances. Using these analytic formulas, we can study the pros and cons of conducting MR analysis under one-sample and two-sample settings and assess the impact of having overlapping samples. We also propose and validate a bias-corrected estimator for the causal effect.     

Production and Characterization of Bacillus thuringiensis Cry1Ac-Resistant Cotton Bollworm Helicoverpa zea (Boddie)

Laboratory-selected Bacillus thuringiensis-resistant colonies are important tools for elucidating B. thuringiensis resistance mechanisms. However, cotton bollworm, Helicoverpa zea, a target pest of transgenic corn and cotton expressing B. thuringiensis Cry1Ac (Bt corn and cotton), has proven difficult to select for stable resistance. Two populations of H. zea (AR and MR), resistant to the B. thuringiensis protein found in all commercial Bt cotton varieties (Cry1Ac), were established by selection with Cry1Ac activated toxin (AR) or MVP II (MR). Cry1Ac toxin reflects the form ingested by H. zea when feeding on Bt cotton, whereas MVP II is a Cry1Ac formulation used for resistance selection and monitoring. The resistance ratio (RR) for AR exceeded 100-fold after 11 generations and has been maintained at this level for nine generations. This is the first report of stable Cry1Ac resistance in H. zea. MR crashed after 11 generations, reaching only an RR of 12. AR was only partially cross-resistant to MVP II, suggesting that MVP II does not have the same Cry1Ac selection pressure as Cry1Ac toxin against H. zea and that proteases may be involved with resistance. AR was highly cross-resistant to Cry1Ab toxin but only slightly cross-resistant to Cry1Ab expressing corn leaf powder. AR was not cross-resistant to Cry2Aa2, Cry2Ab2-expressing corn leaf powder, Vip3A, and cypermethrin. Toxin-binding assays showed no significant differences, indicating that resistance was not linked to a reduction in binding. These results aid in understanding why this pest has not evolved B. thuringiensis resistance, and highlight the need to choose carefully the form of B. thuringiensis protein used in experiments.     

Doubly robust tests of exposure effects under high-dimensional confounding

After variable selection, standard inferential procedures for regression parameters may not be uniformly valid; there is no finite-sample size at which a standard test is guaranteed to approximately attain its nominal size. This problem is exacerbated in high-dimensional settings, where variable selection becomes unavoidable. This has prompted a flurry of activity in developing uniformly valid hypothesis tests for a low-dimensional regression parameter (eg, the causal effect of an exposure A on an outcome Y) in high-dimensional models. So far there has been limited focus on model misspecification, although this is inevitable in high-dimensional settings. We propose tests of the null that are uniformly valid under sparsity conditions weaker than those typically invoked in the literature, assuming working models for the exposure and outcome are both correctly specified. When one of the models is misspecified, by amending the procedure for estimating the nuisance parameters, our tests continue to be valid; hence, they are doubly robust. Our proposals are straightforward to implement using existing software for penalized maximum likelihood estimation and do not require sample splitting. We illustrate them in simulations and an analysis of data obtained from the Ghent University intensive care unit.     

Pleiotropy Can Be Effectively Estimated Without Counting Phenotypes Through the Rank of a Genotype–Phenotype Map

Although pleiotropy, the capability of a gene to affect multiple phenotypes, has been well known as one of the common gene properties, a quantitative estimation remains a great challenge, simply because of the phenotype complexity. Not surprisingly, it is hard for general readers to understand how, without counting phenotypes, gene pleiotropy can be effectively estimated from the genetics data. In this article we extensively discuss the Gu-2007 method that estimated pleiotropy from the protein sequence analysis. We show that this method is actually to estimate the rank (K) of genotype–phenotype mapping that can be concisely written as K = min(r, P_min), where P_min is the minimum pleiotropy among all legitimate measures including the fitness components, and r is the rank of mutational effects of an amino acid site. Together, the effective gene pleiotropy (K_e) estimated by the Gu-2007 method has the following meanings: (i) K_e is an estimate of K = min(r, P_min), the rank of a genotype–phenotype map; (ii) K_e is an estimate for the minimum pleiotropy P_min only if P_min < r; (iii) the Gu-2007 method attempted to estimate the pleiotropy of amino acid sites, a conserved proxy to the true gene pleiotropy; (iv) with a sufficiently large phylogeny such that the rank of mutational effects at an amino acid site is r → 19, one can estimate P_min between 1 and 19; and (v) K_e is a conserved estimate of K because those slightly affected components in fitness have been effectively removed by the estimation procedure. In addition, we conclude that mutational pleiotropy (number of traits affected by a single mutation) cannot be estimated without knowing the phenotypes.