Promises and challenges of populational proteomics in health and disease

Advances in proteomic assay technologies has significantly increased coverage and throughput, enabling recent increases in the number of large-scale population-based proteomic studies of human plasma and serum. Improvement in multiplexed protein assays have facilitated quantification of thousands of proteins over a large dynamic range, a key requirement for detecting the lowest ranging, and potentially the most disease relevant, blood circulating proteins. In this perspective, we examine how populational proteomic datasets in conjunction with other concurrent omic measures can be leveraged to better understand the genomic and non-genomic correlates of the soluble proteome, constructing biomarker panels for disease prediction, among others. Mass spectrometry workflows are discussed as they are becoming increasingly competitive with the affinity-based array platforms in terms of speed, cost and proteome coverage due to advances in both instrumentation and workflows. Despite much success, there remains considerable challenges such as orthogonal validation and absolute quantification. We also highlight emergent challenges associated with study design, analytical considerations and data integration as population scale studies are run in batches and may involve longitudinal samples collated over many years. Lastly, we take a look at the future of what the nascent next-generation proteomic technologies might provide to the analysis of large sets of blood samples, as well as the difficulties in designing large-scale studies that will likely require participation from multiple and complex funding sources and where data sharing, study designs and financing must be solved.

Rare variant associations with plasma protein levels in the UK Biobank

Integrating human genomics and proteomics can help elucidate disease mechanisms, identify clinical biomarkers and discover drug targets1-4. Because previous proteogenomic studies have focused on common variation via genome-wide association studies, the contribution of rare variants to the plasma proteome remains largely unknown. Here we identify associations between rare protein-coding variants and 2,923 plasma protein abundances measured in 49,736 UK Biobank individuals. Our variant-level exome-wide association study identified 5,433 rare genotype-protein associations, of which 81% were undetected in a previous genome-wide association study of the same cohort5. We then looked at aggregate signals using gene-level collapsing analysis, which revealed 1,962 gene-protein associations. Of the 691 gene-level signals from protein-truncating variants, 99.4% were associated with decreased protein levels. STAB1 and STAB2, encoding scavenger receptors involved in plasma protein clearance, emerged as pleiotropic loci, with 77 and 41 protein associations, respectively. We demonstrate the utility of our publicly accessible resource through several applications. These include detailing an allelic series in NLRC4, identifying potential biomarkers for a fatty liver disease-associated variant in HSD17B13 and bolstering phenome-wide association studies by integrating protein quantitative trait loci with protein-truncating variants in collapsing analyses. Finally, we uncover distinct proteomic consequences of clonal haematopoiesis (CH), including an association between TET2-CH and increased FLT3 levels. Our results highlight a considerable role for rare variation in plasma protein abundance and the value of proteogenomics in therapeutic discovery.

Plasma proteomic associations with genetics and health in the UK Biobank

The Pharma Proteomics Project is a precompetitive biopharmaceutical consortium characterizing the plasma proteomic profiles of 54,219 UK Biobank participants. Here we provide a detailed summary of this initiative, including technical and biological validations, insights into proteomic disease signatures, and prediction modelling for various demographic and health indicators. We present comprehensive protein quantitative trait locus (pQTL) mapping of 2,923 proteins that identifies 14,287 primary genetic associations, of which 81% are previously undescribed, alongside ancestry-specific pQTL mapping in non-European individuals. The study provides an updated characterization of the genetic architecture of the plasma proteome, contextualized with projected pQTL discovery rates as sample sizes and proteomic assay coverages increase over time. We offer extensive insights into trans pQTLs across multiple biological domains, highlight genetic influences on ligand-receptor interactions and pathway perturbations across a diverse collection of cytokines and complement networks, and illustrate long-range epistatic effects of ABO blood group and FUT2 secretor status on proteins with gastrointestinal tissue-enriched expression. We demonstrate the utility of these data for drug discovery by extending the genetic proxied effects of protein targets, such as PCSK9, on additional endpoints, and disentangle specific genes and proteins perturbed at loci associated with COVID-19 susceptibility. This public-private partnership provides the scientific community with an open-access proteomics resource of considerable breadth and depth to help to elucidate the biological mechanisms underlying proteo-genomic discoveries and accelerate the development of biomarkers, predictive models and therapeutics1.

Associations of plasma proteomics with type 2 diabetes and related traits: results from the longitudinal KORA S4/F4/FF4 Study

AIMS/HYPOTHESIS

This study aimed to elucidate the aetiological role of plasma proteins in glucose metabolism and type 2 diabetes development.

METHODS

We measured 233 proteins at baseline in 1653 participants from the Cooperative Health Research in the Region of Augsburg (KORA) S4 cohort study (median follow-up time: 13.5 years). We used logistic regression in the cross-sectional analysis (n=1300), and Cox regression accounting for interval-censored data in the longitudinal analysis (n=1143). We further applied two-level growth models to investigate associations with repeatedly measured traits (fasting glucose, 2 h glucose, fasting insulin, HOMA-B, HOMA-IR, HbA1c), and two-sample Mendelian randomisation analysis to investigate causal associations. Moreover, we built prediction models using priority-Lasso on top of Framingham-Offspring Risk Score components and evaluated the prediction accuracy through AUC.

RESULTS

We identified 14, 24 and four proteins associated with prevalent prediabetes (i.e. impaired glucose tolerance and/or impaired fasting glucose), prevalent newly diagnosed type 2 diabetes and incident type 2 diabetes, respectively (28 overlapping proteins). Of these, IL-17D, IL-18 receptor 1, carbonic anhydrase-5A, IL-1 receptor type 2 (IL-1RT2) and matrix extracellular phosphoglycoprotein were novel candidates. IGF binding protein 2 (IGFBP2), lipoprotein lipase (LPL) and paraoxonase 3 (PON3) were inversely associated while fibroblast growth factor 21 was positively associated with incident type 2 diabetes. LPL was longitudinally linked with change in glucose-related traits, while IGFBP2 and PON3 were linked with changes in both insulin- and glucose-related traits. Mendelian randomisation analysis suggested causal effects of LPL on type 2 diabetes and fasting insulin. The simultaneous addition of 12 priority-Lasso-selected biomarkers (IGFBP2, IL-18, IL-17D, complement component C1q receptor, V-set and immunoglobulin domain-containing protein 2, IL-1RT2, LPL, CUB domain-containing protein 1, vascular endothelial growth factor D, PON3, C-C motif chemokine 4 and tartrate-resistant acid phosphatase type 5) significantly improved the predictive performance (ΔAUC 0.0219; 95% CI 0.0052, 0.0624).

CONCLUSIONS/INTERPRETATION

We identified new candidates involved in the development of derangements in glucose metabolism and type 2 diabetes and confirmed previously reported proteins. Our findings underscore the importance of proteins in the pathogenesis of type 2 diabetes and the identified putative proteins can function as potential pharmacological targets for diabetes treatment and prevention.

Transcriptomic and Functional Approaches Unveil the Role of tmRNA in Zinc Acetate Mediated Levofloxacin Sensitivity in Helicobacter pylori

Rapid increase in resistance of Helicobacter pylori (H. pylori) has hindered antibiotics-based eradication efforts worldwide and raises the need for additional approaches. Here, we investigate the role of zinc-based compounds in inhibiting H. pylori growth and modulating antibiotic sensitivities, interrogate their downstream transcriptomic changes, and highlight the potential mechanism driving the observed effects. We showed that zinc acetate inhibited H. pylori growth and increased H. pylori sensitivity to levofloxacin. Transcriptomic profiling showed distinct gene expression patterns between zinc acetate treated groups versus controls. In particular, we independently replicated the association between zinc acetate treatment and increased ssrA expression. Knockdown of ssrA restored levofloxacin resistance to levels of the control group. In this study, we first demonstrated the role of zinc acetate in H. pylori growth and antibiotic sensitivities. Additionally, we explored the transcriptomic perturbations of zinc acetate followed by functional knockdown follow-up of differentially expressed ssrA, highlighting the role of tmRNA and trans-translation in H. pylori levofloxacin resistance. Our results provide alternative and complementary strategies for H. pylori treatment and shed light on the underlying mechanisms driving these effects. IMPORTANCE Helicobacter pylori (H. pylori) eradication plays an important role in gastric cancer prevention, but the antimicrobial resistance of H. pylori is fast becoming a growing concern. In this study, we investigated the role of zinc acetate in inhibiting H. pylori growth and modulating antibiotic sensitivities in vitro. Additionally, we explored the transcriptomic perturbations of zinc acetate followed by functional knockdown follow-up of differentially expressed ssrA, highlighting the role of tmRNA and trans-translation in H. pylori levofloxacin resistance. Our results open up a new horizon for the treatment of antibiotic-resistant H. pylori.

Genetic map of regional sulcal morphology in the human brain from UK biobank data

Genetic associations with macroscopic brain structure can provide insights into brain function and disease. However, specific associations with measures of local brain folding are largely under-explored. Here, we conducted large-scale genome- and exome-wide associations of regional cortical sulcal measures derived from magnetic resonance imaging scans of 40,169 individuals in UK Biobank. We discovered 388 regional brain folding associations across 77 genetic loci, with genes in associated loci enriched for expression in the cerebral cortex, neuronal development processes, and differential regulation during early brain development. We integrated brain eQTLs to refine genes for various loci, implicated several genes involved in neurodevelopmental disorders, and highlighted global genetic correlations with neuropsychiatric phenotypes. We provide an interactive 3D visualisation of our summary associations, emphasising added resolution of regional analyses. Our results offer new insights into the genetic architecture of brain folding and provide a resource for future studies of sulcal morphology in health and disease.

Effect of Temperature on Metronidazole Resistance in Helicobacter pylori

Efficacy of Helicobacter pylori (H. pylori) eradication therapy has declined due to rapid rises in antibiotic resistance. We investigated how increased temperature affected H. pylori (NCTC 11637) growth and its sensitivity to metronidazole in vitro. We performed transcriptomic profiling using RNA-sequencing to identify differentially expressed genes (DEGs) associated with increased temperature. Transcriptional pathways involved in temperature-driven metronidazole resistance changes were analyzed through bioinformatic and literature curation approaches. We showed that H. pylori growth was inhibited at 41°C and inhibition was more apparent with prolonged incubation. Resistance to metronidazole was also reduced—minimum inhibitory concentration for metronidazole decreased from > 256 μg/ml at 37°C to 8 μg/ml at 41°C after culturing for 3 days. RNA-sequencing results, which were highly concordant within treatment conditions, revealed more than one third of genes (583/1,552) to be differentially expressed at increased temperatures with similar proportions up and down-regulated. Quantitative real-time PCR validation for 8 out of 10 DEGs tested gave consistent direction in gene expression changes. We found enrichment for redox and oxygen radical pathways, highlighting a mechanistic pathway driving temperature-related metronidazole resistance. Independent literature review of published genes associated with metronidazole resistance revealed 46 gene candidates, 21 of which showed differential expression and 7 out of 9 DEGs associated with “redox” resistance pathways. Sanger sequencing did not detect any changes in genetic sequences for known resistance genes rdxA, frxA nor fdxB. Our findings suggest that temperature increase can inhibit the growth and reduce H. pylori resistance to metronidazole. Redox pathways are possible potential drivers in metronidazole resistance change induced by temperature. Our study provides insight into potential novel approaches in treating antibiotic resistant H. pylori.

A fast and efficient colocalization algorithm for identifying shared genetic risk factors across multiple traits

Genome-wide association studies (GWAS) have identified thousands of genomic regions affecting complex diseases. The next challenge is to elucidate the causal genes and mechanisms involved. One approach is to use statistical colocalization to assess shared genetic aetiology across multiple related traits (e.g. molecular traits, metabolic pathways and complex diseases) to identify causal pathways, prioritize causal variants and evaluate pleiotropy. We propose HyPrColoc (Hypothesis Prioritisation for multi-trait Colocalization), an efficient deterministic Bayesian algorithm using GWAS summary statistics that can detect colocalization across vast numbers of traits simultaneously (e.g. 100 traits can be jointly analysed in around 1 s). We perform a genome-wide multi-trait colocalization analysis of coronary heart disease (CHD) and fourteen related traits, identifying 43 regions in which CHD colocalized with ≥1 trait, including 5 previously unknown CHD loci. Across the 43 loci, we further integrate gene and protein expression quantitative trait loci to identify candidate causal genes.

[Post-discharge growth of extremely premature infants within corrected age of 24 months]

Objective: To demonstrate the post-discharge catch-up growth of extremely premature infants (EPI) within 24 months of corrected age. Methods: This study retrospectively collected the anthropomorphic measurements of 311 EPI who visited Shenzhen Maternity and Child Healthcare Hospital from August 2013 to April 2020. These infants were stratified according to gestational age at birth (GA): 23-24⁺⁶weeks, 25-26⁺⁶weeks, 27-27⁺⁶weeks; and birth weight:<750 g, 750-999 g, ≥1 000 g. The anthropomorphic measurements, including weight, length, and head circumference for age, were recorded timely from discharge to 24 months of corrected age. And the growth curve stratified by GA and birth weight were fitted in both chronological age and corrected age, which were then compared with the World Health Organization Child Growth Standards for term infant (2006 version), to investigate the catch-up growth pattern of EPI. And appropriate catch-up was defined as the measurements reached the 25^th percentile of WHO growth curve. Results: In these 311 EPI, 184 were males and 127 females, with gestational age of 23-27⁺⁶ weeks and birth weight of 480-1 430 g. Regardless of the GA and birth weight, the growth curves fitted in corrected age failed to overlap with that in chronological age by 24 months of corrected age. The growth velocity of weight, length and head circumference in both corrected and chronological age were all positively correlated with GA and birth weight: the 27-27⁺⁶weeks group showed a preferable growth pattern than the 25-26⁺⁶weeks group, and the curve of the 23-24⁺⁶weeks group was most unfavorable; and the same pattern was observed between the subgroups of different birth weight. Furthermore, the GA had more significant impact on the catch-up growth pattern than birth weight did. When assessed with corrected age curve, the weight and length of both male and female EPIs achieved appropriate catch-up by 24 months, as well as the head circumference of girls; whereas, boys' head circumference reached appropriate catch-up at the corrected age of 9 months, but fell behind the 25^th percentile after that. However, when assessed with chronological age curve, both boys and girls failed to achieve appropriate catch-up in weight, length and head circumference by age 24 months. And no matter in corrected or chronological age, all physical measurements of girls were lower than those of boys. Conclusions: The rapid catch-up growth of EPI happens within 6 months of corrected age. The lower the birth weight and gestational age, the lower the physical measurements at each corresponding month of age, and the longer it takes to achieve appropriate catch-up. Gestational age has a greater impact on the longitudinal catch-up growth than birth weight does. And girls generally grow slower than boys in either correct or actual age. Before 24 months of corrected age, the growth should be assessed with corrected age rather than chronological age.

Phenome-wide Mendelian randomization mapping the influence of the plasma proteome on complex diseases

The human proteome is a major source of therapeutic targets. Recent genetic association analyses of the plasma proteome enable systematic evaluation of the causal consequences of variation in plasma protein levels. Here we estimated the effects of 1,002 proteins on 225 phenotypes using two-sample Mendelian randomization (MR) and colocalization. Of 413 associations supported by evidence from MR, 130 (31.5%) were not supported by results of colocalization analyses, suggesting that genetic confounding due to linkage disequilibrium is widespread in naïve phenome-wide association studies of proteins. Combining MR and colocalization evidence in cis-only analyses, we identified 111 putatively causal effects between 65 proteins and 52 disease-related phenotypes ( https://www.epigraphdb.org/pqtl/ ). Evaluation of data from historic drug development programs showed that target-indication pairs with MR and colocalization support were more likely to be approved, evidencing the value of this approach in identifying and prioritizing potential therapeutic targets.

Interleukin-6 Receptor Signaling and Abdominal Aortic Aneurysm Growth Rates

Supplemental Digital Content is available in the text.

Background:

The Asp358Ala variant (rs2228145; A>C) in the IL (interleukin)-6 receptor (IL6R) gene has been implicated in the development of abdominal aortic aneurysms (AAAs), but its effect on AAA growth over time is not known. We aimed to investigate the clinical association between the IL6R-Asp358Ala variant and AAA growth and to assess the effect of blocking the IL-6 signaling pathway in mouse models of aortic aneurysm rupture or dissection.

Methods:

Using data from 2863 participants with AAA from 9 prospective cohorts, age- and sex-adjusted mixed-effects linear regression models were used to estimate the association between the IL6R-Asp358Ala variant and annual change in AAA diameter (mm/y). In a series of complementary randomized trials in mice, the effect of blocking the IL-6 signaling pathways was assessed on plasma biomarkers, systolic blood pressure, aneurysm diameter, and time to aortic rupture and death.

Results:

After adjusting for age and sex, baseline aneurysm size was 0.55 mm (95% CI, 0.13–0.98 mm) smaller per copy of the minor allele [C] of the Asp358Ala variant. Change in AAA growth was −0.06 mm per year (−0.18 to 0.06) per copy of the minor allele; a result that was not statistically significant. Although all available worldwide data were used, the genetic analyses were not powered for an effect size as small as that observed. In 2 mouse models of AAA, selective blockage of the IL-6 trans-signaling pathway, but not combined blockage of both, the classical and trans-signaling pathways, was associated with improved survival (P<0.05).

Conclusions:

Our proof-of-principle data are compatible with the concept that IL-6 trans-signaling is relevant to AAA growth, encouraging larger-scale evaluation of this hypothesis.

Genome-wide association meta-analyses and fine-mapping elucidate pathways influencing albuminuria

Increased levels of the urinary albumin-to-creatinine ratio (UACR) are associated with higher risk of kidney disease progression and cardiovascular events, but underlying mechanisms are incompletely understood. Here, we conduct trans-ethnic (n = 564,257) and European-ancestry specific meta-analyses of genome-wide association studies of UACR, including ancestry- and diabetes-specific analyses, and identify 68 UACR-associated loci. Genetic correlation analyses and risk score associations in an independent electronic medical records database (n = 192,868) reveal connections with proteinuria, hyperlipidemia, gout, and hypertension. Fine-mapping and trans-Omics analyses with gene expression in 47 tissues and plasma protein levels implicate genes potentially operating through differential expression in kidney (including TGFB1, MUC1, PRKCI, and OAF), and allow coupling of UACR associations to altered plasma OAF concentrations. Knockdown of OAF and PRKCI orthologs in Drosophila nephrocytes reduces albumin endocytosis. Silencing fly PRKCI further impairs slit diaphragm formation. These results generate a priority list of genes and pathways for translational research to reduce albuminuria.

Author Correction: New genetic signals for lung function highlight pathways and chronic obstructive pulmonary disease associations across multiple ancestries

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

New genetic signals for lung function highlight pathways and chronic obstructive pulmonary disease associations across multiple ancestries

Reduced lung function predicts mortality and is key to the diagnosis of chronic obstructive pulmonary disease (COPD). In a genome-wide association study in 400,102 individuals of European ancestry, we define 279 lung function signals, 139 of which are new. In combination, these variants strongly predict COPD in independent populations. Furthermore, the combined effect of these variants showed generalizability across smokers and never smokers, and across ancestral groups. We highlight biological pathways, known and potential drug targets for COPD and, in phenome-wide association studies, autoimmune-related and other pleiotropic effects of lung function-associated variants. This new genetic evidence has potential to improve future preventive and therapeutic strategies for COPD.

ProGeM: a framework for the prioritization of candidate causal genes at molecular quantitative trait loci

Quantitative trait locus (QTL) mapping of molecular phenotypes such as metabolites, lipids and proteins through genome-wide association studies represents a powerful means of highlighting molecular mechanisms relevant to human diseases. However, a major challenge of this approach is to identify the causal gene(s) at the observed QTLs. Here, we present a framework for the 'Prioritization of candidate causal Genes at Molecular QTLs' (ProGeM), which incorporates biological domain-specific annotation data alongside genome annotation data from multiple repositories. We assessed the performance of ProGeM using a reference set of 227 previously reported and extensively curated metabolite QTLs. For 98% of these loci, the expert-curated gene was one of the candidate causal genes prioritized by ProGeM. Benchmarking analyses revealed that 69% of the causal candidates were nearest to the sentinel variant at the investigated molecular QTLs, indicating that genomic proximity is the most reliable indicator of 'true positive' causal genes. In contrast, cis-gene expression QTL data led to three false positive candidate causal gene assignments for every one true positive assignment. We provide evidence that these conclusions also apply to other molecular phenotypes, suggesting that ProGeM is a powerful and versatile tool for annotating molecular QTLs. ProGeM is freely available via GitHub.

Genome-wide mapping of plasma protein QTLs identifies putatively causal genes and pathways for cardiovascular disease

Identifying genetic variants associated with circulating protein concentrations (protein quantitative trait loci; pQTLs) and integrating them with variants from genome-wide association studies (GWAS) may illuminate the proteome's causal role in disease and bridge a knowledge gap regarding SNP-disease associations. We provide the results of GWAS of 71 high-value cardiovascular disease proteins in 6861 Framingham Heart Study participants and independent external replication. We report the mapping of over 16,000 pQTL variants and their functional relevance. We provide an integrated plasma protein-QTL database. Thirteen proteins harbor pQTL variants that match coronary disease-risk variants from GWAS or test causal for coronary disease by Mendelian randomization. Eight of these proteins predict new-onset cardiovascular disease events in Framingham participants. We demonstrate that identifying pQTLs, integrating them with GWAS results, employing Mendelian randomization, and prospectively testing protein-trait associations holds potential for elucidating causal genes, proteins, and pathways for cardiovascular disease and may identify targets for its prevention and treatment.

Genomic atlas of the human plasma proteome

Although plasma proteins have important roles in biological processes and are the direct targets of many drugs, the genetic factors that control inter-individual variation in plasma protein levels are not well understood. Here we characterize the genetic architecture of the human plasma proteome in healthy blood donors from the INTERVAL study. We identify 1,927 genetic associations with 1,478 proteins, a fourfold increase on existing knowledge, including trans associations for 1,104 proteins. To understand the consequences of perturbations in plasma protein levels, we apply an integrated approach that links genetic variation with biological pathway, disease, and drug databases. We show that protein quantitative trait loci overlap with gene expression quantitative trait loci, as well as with disease-associated loci, and find evidence that protein biomarkers have causal roles in disease using Mendelian randomization analysis. By linking genetic factors to diseases via specific proteins, our analyses highlight potential therapeutic targets, opportunities for matching existing drugs with new disease indications, and potential safety concerns for drugs under development.

[Association of sleep-disorders with systemic inflammatory markers and acute exacerbation in chronic obstructive pulmonary disease patients]

Objective: To investigate the association of sleep quality with the levels of systemic inflammatory markers in patients with chronic obstructive pulmonary disease(COPD) and the correlations between the frequency of acute exacerbation of COPD (AECOPD) and Pittsburgh sleep quality index (PSQI). Methods: A total of 198 COPD patients admitted in our hospital from October, 2016 to June, 2017 were screened, and 124 patients were eligible for the study. On the first day of hospitalization, the serum samples and clinical data were collected, including white blood cells, lymphocytes, platelet count, CRP and PSQI. Poor sleep quality was defined as PSQI score >5. Results: The percentage of COPD patients with poor sleep quality was about 68%. Poor sleep quality was positively correlated with the frequency of acute exacerbation in COPD patients. The ratio of neutrophil to lymphocyte (NLR), ratio of platelet to lymphocyte (PLR) and levels of CRP were higher in patients with poor sleep quality than those in the control group. NLR, PLR and CRP in peripheral blood of the patients with poor sleep quality were positively correlated with PSQI score. The CRP levels and PSQI score in COPD patients with poor sleep quality group were positively correlated with the frequency of exacerbations in the past year (r=0.437, r=0.430). Conclusion: A high percentage of COPD patients had poor sleep quality, which was positively correlated with the levels of systemic inflammation as well as the frequency of AECOPD.

Mendelian randomization with fine-mapped genetic data: Choosing from large numbers of correlated instrumental variables

Mendelian randomization uses genetic variants to make causal inferences about the effect of a risk factor on an outcome. With fine-mapped genetic data, there may be hundreds of genetic variants in a single gene region any of which could be used to assess this causal relationship. However, using too many genetic variants in the analysis can lead to spurious estimates and inflated Type 1 error rates. But if only a few genetic variants are used, then the majority of the data is ignored and estimates are highly sensitive to the particular choice of variants. We propose an approach based on summarized data only (genetic association and correlation estimates) that uses principal components analysis to form instruments. This approach has desirable theoretical properties: it takes the totality of data into account and does not suffer from numerical instabilities. It also has good properties in simulation studies: it is not particularly sensitive to varying the genetic variants included in the analysis or the genetic correlation matrix, and it does not have greatly inflated Type 1 error rates. Overall, the method gives estimates that are less precise than those from variable selection approaches (such as using a conditional analysis or pruning approach to select variants), but are more robust to seemingly arbitrary choices in the variable selection step. Methods are illustrated by an example using genetic associations with testosterone for 320 genetic variants to assess the effect of sex hormone related pathways on coronary artery disease risk, in which variable selection approaches give inconsistent inferences.

PhenoScanner: a database of human genotype-phenotype associations

Summary: PhenoScanner is a curated database of publicly available results from large-scale genetic association studies. This tool aims to facilitate ‘phenome scans’, the cross-referencing of genetic variants with many phenotypes, to help aid understanding of disease pathways and biology. The database currently contains over 350 million association results and over 10 million unique genetic variants, mostly single nucleotide polymorphisms. It is accompanied by a web-based tool that queries the database for associations with user-specified variants, providing results according to the same effect and non-effect alleles for each input variant. The tool provides the option of searching for trait associations with proxies of the input variants, calculated using the European samples from 1000 Genomes and Hapmap.

Availability and Implementation: PhenoScanner is available at www.phenoscanner.medschl.cam.ac.uk.

Contact:jrs95@medschl.cam.ac.uk

Supplementary information:Supplementary data are available at Bioinformatics online.

Identification of natural recombinants derived from PCV2a and PCV2b

Porcine circovirus type 2 (PCV2) is considered to be the main pathogen in PC-associated diseases, and significantly affects the global pig-producing industry. PCV2 continuously evolves by point mutations and genome recombinations. In the present study, we aimed to further identify recombinant PCV2 strains. We used polymerase chain reaction to detect PCV2 in the carcasses of pigs with suspected infections from different regions of Guangdong Province in China. DNA was extracted from samples with confirmed infection and full- genome amplification, sequencing, phylogenetic tree construction, gene recombination detection, and sequence alignment were performed in gene recombination analysis. Our results show that recombination occurred between the strains SHC (DQ104421) and ZhuJi2003 (AY579893). The recombination resulted in three recombinants: GD003 (KM503044), GD005 (KM487708), and GD008 (KM487709). Further analyses revealed that these novel recombinants appeared to result from recombination between the PCV2a and PCV2b strains, with crossover regions located in ORF2. This study was a comprehensive analysis that used several different methods, which demonstrated that a cluster of PCV2 strains resulted from the same type of inter-genotypic recombination pattern, with a breakpoint in the structural protein coding region. The results of our study provide both information on the recombination mechanism and disease pathogenesis and useful data for the prevention of PCV2 in the swine industry.

Conventional protein kinase C isoforms differentially regulate ADP- and thrombin-evoked Ca²⁺ signalling in human platelets

Rises in cytosolic Ca(2+) concentration ([Ca(2+)]cyt) are central in platelet activation, yet many aspects of the underlying mechanisms are poorly understood. Most studies examine how experimental manipulations affect agonist-evoked rises in [Ca(2+)]cyt, but these only monitor the net effect of manipulations on the processes controlling [Ca(2+)]cyt (Ca(2+) buffering, sequestration, release, entry and removal), and cannot resolve the source of the Ca(2+) or the transporters or channels affected. To investigate the effects of protein kinase C (PKC) on platelet Ca(2+) signalling, we here monitor Ca(2+) flux around the platelet by measuring net Ca(2+) fluxes to or from the extracellular space and the intracellular Ca(2+) stores, which act as the major sources and sinks for Ca(2+) influx into and efflux from the cytosol, as well as monitoring the cytosolic Na(+) concentration ([Na(+)]cyt), which influences platelet Ca(2+) fluxes via Na(+)/Ca(2+) exchange. The intracellular store Ca(2+) concentration ([Ca(2+)]st) was monitored using Fluo-5N, the extracellular Ca(2+) concentration ([Ca(2+)]ext) was monitored using Fluo-4 whilst [Ca(2+)]cyt and [Na(+)]cyt were monitored using Fura-2 and SFBI, respectively. PKC inhibition using Ro-31-8220 or bisindolylmaleimide I potentiated ADP- and thrombin-evoked rises in [Ca(2+)]cyt in the absence of extracellular Ca(2+). PKC inhibition potentiated ADP-evoked but reduced thrombin-evoked intracellular Ca(2+) release and Ca(2+) removal into the extracellular medium. SERCA inhibition using thapsigargin and 2,5-di(tert-butyl) l,4-benzohydroquinone abolished the effect of PKC inhibitors on ADP-evoked changes in [Ca(2+)]cyt but only reduced the effect on thrombin-evoked responses. Thrombin evokes substantial rises in [Na(+)]cyt which would be expected to reduce Ca(2+) removal via the Na(+)/Ca(2+) exchanger (NCX). Thrombin-evoked rises in [Na(+)]cyt were potentiated by PKC inhibition, an effect which was not due to altered changes in non-selective cation permeability of the plasma membrane as assessed by Mn(2+) quench of Fura-2 fluorescence. PKC inhibition was without effect on thrombin-evoked rises in [Ca(2+)]cyt following SERCA inhibition and either removal of extracellular Na(+) or inhibition of Na(+)/K(+)-ATPase activity by removal of extracellular K(+) or treatment with digoxin. These data suggest that PKC limits ADP-evoked rises in [Ca(2+)]cyt by acceleration of SERCA activity, whilst rises in [Ca(2+)]cyt evoked by the stronger platelet activator thrombin are limited by PKC through acceleration of both SERCA and Na(+)/K(+)-ATPase activity, with the latter limiting the effect of thrombin on rises in [Na(+)]cyt and so forward mode NCX activity. The use of selective PKC inhibitors indicated that conventional and not novel PKC isoforms are responsible for the inhibition of agonist-evoked Ca(2+) signalling.

First Report of Leaf Blotch Caused by Pestalotiopsis clavispora on Rosa chinensis in China

Rosa chinensis Jacq, a traditional Chinese ornamental flower, is an important landscaping plant in northern China. Since July 2013, leaf blotch symptoms were observed in the Tianjin flower nursery (117°09' E 39°17' N). The garden exhibited 40% disease incidence with observable symptoms on basal leaves that were yellowed from the edge to the inside area on infected leaves, in the shape of a V. Yellow lesions covered one third to one half of the leaf. Yellow halos were observed at the junction of the healthy and diseased tissues. Small tissue pieces from the edges of lesions were disinfected in 70% ethyl alcohol for 30 s and 1% hypochlorite for 1 min, rinsed thrice in sterile water, plated on potato dextrose agar (PDA), and incubated at 25°C in lighted incubator for 4 days. Fungal colonies that developed on PDA were white and cottony with concentric rings. Black and globular acervuli appeared after 10 days at 25°C. Conidia (n = 20), which were fusiform, were 9.20 to 31.31 (avg. 26.5) × 4.83 to 9.11 (avg. 6.9) μm. Conidia of all isolates were five celled. Apical and basal cells were colorless, while the three median cells were dark brown. The single basal appendage of conidia was 2.85 to 16.05 μm in length and the two to three apical appendages were 5.93 to 36.23 μm in length. According to colony and conidia morphology (number of cells, number of appendages), the isolates were initially identified as Pestalotiopsis spp (2). A 525-bp band was produced in a conventional PCR assay. Primers ITS1 (5'TCCGTAGGTGAACCTGCGC3') and ITS4 (5'TCCTCCGCTTATTGATATGC3') were used to amplify and sequence the internal transcribed spacer region of rDNA. A BLAST search of the NCBI databases showed that isolate YJYK-1 had 99% homology with Pestalotiopsis clavispora isolate hz-067 (Accession No. FJ517545.1). Pathogenicity tests of the novel isolate YJYK-1 were conducted by placing agar plugs (5 mm in diameter) from an actively growing colony on PDA on surface-disinfected (70% ethyl alcohol, 30 s) leaves (1). Control leaves were inoculated with sterile PDA plugs. Plants with inoculated leaves (five per treatment) were placed in lighted growth chambers at 25°C for 10 days and watered as needed. Symptoms on inoculated leaves were similar to those previously described in the nursery. Black acervuli were easily found on the necrotic tissues. Control plants did not show any symptoms. Cultures isolated from the lesions were similar to those isolated previously from leaves in the nursery. Koch's postulates were confirmed after re-isolation. Although the diversity of endophytic Pestalotiopsis species and its distribution was investigated and the host plants were also listed in China (3), to our knowledge, this is the first report of P. clavispora causing leaf blotch on rose R. chinensis in China. References: (1) M. I. Ahmed et al. Eur. J. Plant Pathol. 135:619, 2013. (2) J. Y. Lu. Diagnosis of plant diseases. Page 194 in: Pestalotiopsis. J. Y. Lu et al., eds. China Agriculture Press, Beijing, 1995. (3) J. G. Wei et al. Mycosystema 24:481, 2005.

Identification of 18 genes encoding necrosis-inducing proteins from the plant pathogen Phytophthora capsici (Pythiaceae: Oomycetes)

Phytophthora capsici is an aggressive plant pathogen that affects solanaceous and cucurbitaceous hosts. Necrosis-inducing Phytophthora proteins (NPPs) are a group of secreted toxins found particularly in oomycetes. Several NPPs from Phytophthora species trigger plant cell death and activate host defense gene expression. We isolated 18 P. capsici NPP genes, of which 12 were active during hypha growth from a Phytophthora stain isolated from pepper (Capsicum annuum) plants in China. The 18 predicted proteins had a sequence homology of 46.26%. The 18 Pcnpp sequences had a conserved GHRHDWE motif and fell into two groups. Eleven sequences in group 1 had two conserved cysteine residues, whereas the other seven sequences in group 2 lacked these two cysteine residues. A phylogenetic tree was constructed on the basis of the alignment of the predicted protein sequences of 52 selected NPP genes from oomycetes, fungi and bacteria from Genbank. The tree did not rigorously follow the taxonomic classification of the species; all the NPPs from oomycetes formed their own clusters, while fungal sequences were grouped into two separate clades, indicating that based on NPPs, we can separate oomycetes from fungi and bacteria, and that expansion of the NPP family was a feature of Phytophthora evolution.

Ductile titanium alloy with low Poisson's ratio

We report a ductile beta-type titanium alloy with body-centered cubic (bcc) crystal structure having a low Poisson's ratio of 0.14. The almost identical ultralow bulk and shear moduli of approximately 24 GPa combined with an ultrahigh strength of approximately 0.9 GPa contribute to easy crystal distortion due to much-weakened chemical bonding of atoms in the crystal, leading to significant elastic softening in tension and elastic hardening in compression. The peculiar elastic and plastic deformation behaviors of the alloy are interpreted as a result of approaching the elastic limit of the bcc crystal under applied stress.

N-type calcium channels and their regulation by GABAB receptors in axons of neonatal rat optic nerve

Axons of neonatal rat optic nerves exhibit fast calcium transients in response to brief action potential stimulation. In response to one to four closely spaced action potentials, evoked calcium transients showed a fast-rising phase followed by a decay with a time constant of approximately 2-3 sec. By selective staining of axons or glial cells with calcium dyes, it was shown that the evoked calcium transient originated from axons. The calcium transient was caused by influx because it was eliminated when bath calcium was removed. Pharmacological profile studies with calcium channel subtype-specific peptides suggested that 58% of the evoked calcium influx was accounted for by N-type calcium channels, whereas L- and P/Q-type calcium channels had little, if any, contribution. The identity of the residual calcium influx remains unclear. GABA application caused a dramatic reduction of the amplitude of the action potential and the associated calcium influx. When GABAA receptors were blocked by bicuculline, the inhibitory effect of GABA on the action potential was eliminated, whereas that on the calcium influx was not, indicating involvement of GABAB receptors. Indeed, the calcium influx was inhibited by the GABAB receptor agonist baclofen. This baclofen effect was occluded by a previous block of N-type calcium channels and was unaffected by the broad-spectrum K+ channel blocker 4-AP. We conclude that neonatal rat optic nerve axons express N-type calcium channels, which are subjected to regulation by G-protein-coupled GABAB receptors. We suggest that receptor-mediated inhibition of axonal calcium channels plays a protective role in neonatal anoxic and/or ischemic injury.