Assessing the clinical utility of cancer genomic and proteomic data across tumor types

Molecular profiling of tumors promises to advance the clinical management of cancer, but the benefits of integrating molecular data with traditional clinical variables have not been systematically studied. Here we retrospectively predict patient survival using diverse molecular data (somatic copy-number alteration, DNA methylation and mRNA, miRNA and protein expression) from 953 samples of four cancer types from The Cancer Genome Atlas project. We found that incorporating molecular data with clinical variables yielded statistically significantly improved predictions (FDR < 0.05) for three cancers but those quantitative gains were limited (2.2–23.9%). Additional analyses revealed little predictive power across tumor types except for one case. In clinically relevant genes, we identified 10,281 somatic alterations across 12 cancer types in 2,928 of 3,277 patients (89.4%), many of which would not be revealed in single-tumor analyses. Our study provides a starting point and resources, including an open-access model evaluation platform, for building reliable prognostic and therapeutic strategies that incorporate molecular data

Restriction associated DNA-genotyping at multiple spatial scales in Arabidopsis lyrata reveals signatures of pathogen-mediated selection

Background: Genome scans based on outlier analyses have revolutionized detection of genes involved in adaptive processes, but reports of some forms of selection, such as balancing selection, are still limited. It is unclear whether high throughput genotyping approaches for identification of single nucleotide polymorphisms have sufficient power to detect modes of selection expected to result in reduced genetic differentiation among populations. In this study, we used Arabidopsis lyrata to investigate whether signatures of balancing selection can be detected based on genomic smoothing of Restriction Associated DNA sequencing (RAD-seq) data. We compared how different sampling approaches (both within and between subspecies) and different background levels of polymorphism (inbreeding or outcrossing populations) affected the ability to detect genomic regions showing key signatures of balancing selection, specifically elevated polymorphism, reduced differentiation and shifts towards intermediate allele frequencies. We then tested whether candidate genes associated with disease resistance (R-gene analogs) were detected more frequently in these regions compared to other regions of the genome. Results: We found that genomic regions showing elevated polymorphism contained a significantly higher density of R-gene analogs predicted to be under pathogen-mediated selection than regions of non-elevated polymorphism, and that many of these also showed evidence for an intermediate site-frequency spectrum based on Tajima’s D. However, we found few genomic regions that showed both elevated polymorphism and reduced FST among populations, despite strong background levels of genetic differentiation among populations. This suggests either insufficient power to detect the reduced population structure predicted for genes under balancing selection using sparsely distributed RAD markers, or that other forms of diversifying selection are more common for the R-gene analogs tested. Conclusions: Genome scans based on a small number of individuals sampled from a wide range of populations were sufficient to confirm the relative scarcity of signatures of balancing selection across the genome, but also identified new potential disease resistance candidates within genomic regions showing signatures of balancing selection that would be strong candidates for further sequencing efforts

LOCAS – A Low Coverage Assembly Tool for Resequencing Projects

Motivation: Next Generation Sequencing (NGS) is a frequently applied approach to detect sequence variations between highly related genomes. Recent large-scale re-sequencing studies as the Human 1000 Genomes Project utilize NGS data of low coverage to afford sequencing of hundreds of individuals. Here, SNPs and micro-indels can be detected by applying an alignment-consensus approach. However, computational methods capable of discovering other variations such as novel insertions or highly diverged sequence from low coverage NGS data are still lacking. Results: We present LOCAS, a new NGS assembler particularly designed for low coverage assembly of eukaryotic genomes using a mismatch sensitive overlap-layout-consensus approach. LOCAS assembles homologous regions in a homologyguided manner while it performs de novo assemblies of insertions and highly polymorphic target regions subsequently to an alignment-consensus approach. LOCAS has been evaluated in homology-guided assembly scenarios with low sequence coverage of Arabidopsis thaliana strains sequenced as part of the Arabidopsis 1001 Genomes Project. While assembling the same amount of long insertions as state-of-the-art NGS assemblers, LOCAS showed best results regarding contig size, error rate and runtime. Conclusion: LOCAS produces excellent results for homology-guided assembly of eukaryotic genomes with short reads and low sequencing depth, and therefore appears to be the assembly tool of choice for the detection of novel sequenc

Toxoplasma seroprevalence in a rural population in France: detection of a household effect

<p>Abstract</p> <p>Background</p> <p><it>Toxoplasma gondii</it>, the agent of toxoplasmosis, has a complex life cycle. In humans, the parasite may be acquired either through ingestion of contaminated meat or through oocysts present in the environment. The importance of each source of contamination varies locally according to the environment characteristics and to differences concerning human eating habits and the presence of cats; thus, the risk factors may be determined through fine-scale studies. Here, we searched for factors associated with seropositivity in the population of two adjacent villages in Lorraine region, France.</p> <p>Methods</p> <p>All voluntary inhabitants filled out a questionnaire and gave a blood sample. The seroprevalence was estimated globally and according to the inhabitants' ages using a cubic spline regression. A mixed logistic regression model was used to quantify the effect of individual and household factors on the probability of seropositivity.</p> <p>Results</p> <p>Based on serological results from 273 persons, we estimated seroprevalence to be 47% (95% confidence interval: 41 to 53%). That seroprevalence increased with age: the slope was the steepest up to the age of 40 years (OR = 2.48 per 10-year increment, 95% credibility interval: [1.29 to 5.09]), but that increase was not significant afterwards. The probability of seropositivity tended to be higher in men than in women (OR = 2.01, 95% credibility interval: [0.92 to 4.72]) and in subjects eating raw vegetables at least once a week than in the others (OR = 8.4, 95% credibility interval: [0.93 to 72.1]). These effects were close to statistical significance. The multivariable analysis highlighted a significant seroprevalence heterogeneity among households. That seroprevalence varied between 6 and 91% (5^thand 95^thpercentile of the household seropositivity distribution).</p> <p>Conclusion</p> <p>The major finding is the household effect, with a strong heterogeneity of seroprevalence among households. This effect may be explained by common exposures of household members to local risk factors. Future work will quantify the link between the presence of oocysts in the soil and the seroprevalence of exposed households using a spatial analysis.</p

Restriction Site Extension PCR: A Novel Method for High-Throughput Characterization of Tagged DNA Fragments and Genome Walking

BACKGROUND: Insertion mutant isolation and characterization are extremely valuable for linking genes to physiological function. Once an insertion mutant phenotype is identified, the challenge is to isolate the responsible gene. Multiple strategies have been employed to isolate unknown genomic DNA that flanks mutagenic insertions, however, all these methods suffer from limitations due to inefficient ligation steps, inclusion of restriction sites within the target DNA, and non-specific product generation. These limitations become close to insurmountable when the goal is to identify insertion sites in a high throughput manner. METHODOLOGY/PRINCIPAL FINDINGS: We designed a novel strategy called Restriction Site Extension PCR (RSE-PCR) to efficiently conduct large-scale isolation of unknown genomic DNA fragments linked to DNA insertions. The strategy is a modified adaptor-mediated PCR without ligation. An adapter, with complementarity to the 3' overhang of the endonuclease (KpnI, NsiI, PstI, or SacI) restricted DNA fragments, extends the 3' end of the DNA fragments in the first cycle of the primary RSE-PCR. During subsequent PCR cycles and a second semi-nested PCR (secondary RSE-PCR), touchdown and two-step PCR are combined to increase the amplification specificity of target fragments. The efficiency and specificity was demonstrated in our characterization of 37 tex mutants of Arabidopsis. All the steps of RSE-PCR can be executed in a 96 well PCR plate. Finally, RSE-PCR serves as a successful alternative to Genome Walker as demonstrated by gene isolation from maize, a plant with a more complex genome than Arabidopsis. CONCLUSIONS/SIGNIFICANCE: RSE-PCR has high potential application in identifying tagged (T-DNA or transposon) sequence or walking from known DNA toward unknown regions in large-genome plants, with likely application in other organisms as well

Generation of ESTs for Flowering Gene Discovery and SSR Marker Development in Upland Cotton

BACKGROUND: Upland cotton, Gossypium hirsutum L., is one of the world's most important economic crops. In the absence of the entire genomic sequence, a large number of expressed sequence tag (EST) resources of upland cotton have been generated and used in several studies. However, information about the flower development of this species is rare. METHODOLOGY/PRINCIPAL FINDINGS: To clarify the molecular mechanism of flower development in upland cotton, 22,915 high-quality ESTs were generated and assembled into 14,373 unique sequences consisting of 4,563 contigs and 9,810 singletons from a normalized and full-length cDNA library constructed from pooled RNA isolated from shoot apexes, squares, and flowers. Comparative analysis indicated that 5,352 unique sequences had no high-degree matches to the cotton public database. Functional annotation showed that several upland cotton homologs with flowering-related genes were identified in our library. The majority of these genes were specifically expressed in flowering-related tissues. Three GhSEP (G. hirsutum L. SEPALLATA) genes determining floral organ development were cloned, and quantitative real-time PCR (qRT-PCR) revealed that these genes were expressed preferentially in squares or flowers. Furthermore, 670 new putative microsatellites with flanking sequences sufficient for primer design were identified from the 645 unigenes. Twenty-five EST-simple sequence repeats were randomly selected for validation and transferability testing in 17 Gossypium species. Of these, 23 were identified as true-to-type simple sequence repeat loci and were highly transferable among Gossypium species. CONCLUSIONS/SIGNIFICANCE: A high-quality, normalized, full-length cDNA library with a total of 14,373 unique ESTs was generated to provide sequence information for gene discovery and marker development related to upland cotton flower development. These EST resources form a valuable foundation for gene expression profiling analysis, functional analysis of newly discovered genes, genetic linkage, and quantitative trait loci analysis

Bias adjustment and ensemble recalibration methods for seasonal forecasting: a comprehensive intercomparison using the C3S dataset

This work presents a comprehensive intercomparison of diferent alternatives for the calibration of seasonal forecasts, ranging from simple bias adjustment (BA)-e.g. quantile mapping-to more sophisticated ensemble recalibration (RC) methods- e.g. non-homogeneous Gaussian regression, which build on the temporal correspondence between the climate model and the corresponding observations to generate reliable predictions. To be as critical as possible, we validate the raw model and the calibrated forecasts in terms of a number of metrics which take into account diferent aspects of forecast quality (association, accuracy, discrimination and reliability). We focus on one-month lead forecasts of precipitation and temperature from four state-of-the-art seasonal forecasting systems, three of them included in the Copernicus Climate Change Service dataset (ECMWF-SEAS5, UK Met Ofce-GloSea5 and Météo France-System5) for boreal winter and summer over two illustrative regions with diferent skill characteristics (Europe and Southeast Asia). Our results indicate that both BA and RC methods efectively correct the large raw model biases, which is of paramount importance for users, particularly when directly using the climate model outputs to run impact models, or when computing climate indices depending on absolute values/thresholds. However, except for particular regions and/or seasons (typically with high skill), there is only marginal added value-with respect to the raw model outputs-beyond this bias removal. For those cases, RC methods can outperform BA ones, mostly due to an improvement in reliability. Finally, we also show that whereas an increase in the number of members only modestly afects the results obtained from calibration, longer hindcast periods lead to improved forecast quality, particularly for RC methods.This work has been funded by the C3S activity on Evaluation and Quality Control for seasonal forecasts. JMG was partially supported by the project MULTI-SDM (CGL2015-66583-R, MINECO/FEDER). FJDR was partially funded by the H2020 EUCP project (GA 776613)

NVP-AUY922: a small molecule HSP90 inhibitor with potent antitumor activity in preclinical breast cancer models

INTRODUCTION:Heat shock protein 90 (HSP90) is a key component of a multichaperone complex involved in the post-translational folding of a large number of client proteins, many of which play essential roles in tumorigenesis. HSP90 has emerged in recent years as a promising new target for anticancer therapies.METHODS:The concentrations of the HSP90 inhibitor NVP-AUY922 required to reduce cell numbers by 50% (GI50 values) were established in a panel of breast cancer cell lines and patient-derived human breast tumors. To investigate the properties of the compound in vivo, the pharmacokinetic profile, antitumor effect, and dose regimen were established in a BT-474 breast cancer xenograft model. The effect on HSP90-p23 complexes, client protein degradation, and heat shock response was investigated in cell culture and breast cancer xenografts by immunohistochemistry, Western blot analysis, and immunoprecipitation.RESULTS:We show that the novel small molecule HSP90 inhibitor NVP-AUY922 potently inhibits the proliferation of human breast cancer cell lines with GI50 values in the range of 3 to 126 nM. NVP-AUY922 induced proliferative inhibition concurrent with HSP70 upregulation and client protein depletion � hallmarks of HSP90 inhibition. Intravenous acute administration of NVP-AUY922 to athymic mice (30 mg/kg) bearing subcutaneous BT-474 breast tumors resulted in drug levels in excess of 1,000 times the cellular GI50 value for about 2 days. Significant growth inhibition and good tolerability were observed when the compound was administered once per week. Therapeutic effects were concordant with changes in pharmacodynamic markers, including HSP90-p23 dissociation, decreases in ERBB2 and P-AKT, and increased HSP70 protein levels.CONCLUSION:NVP-AUY922 is a potent small molecule HSP90 inhibitor showing significant activity against breast cancer cells in cellular and in vivo settings. On the basis of its mechanism of action, preclinical activity profile, tolerability, and pharmaceutical properties, the compound recently has entered clinical phase I breast cancer trials

Seropositivity and Risk Factors Associated with Toxoplasma gondii Infection in Wild Birds from Spain

Toxoplasma gondii is a zoonotic intracellular protozoan parasite of worldwide distribution that infects many species of warm-blooded animals, including birds. To date, there is scant information about the seropositivity of T. gondii and the risk factors associated with T. gondii infection in wild bird populations. In the present study, T. gondii infection was evaluated on sera obtained from 1079 wild birds belonging to 56 species (including Falconiformes (n = 610), Strigiformes (n = 260), Ciconiiformes (n = 156), Gruiformes (n = 21), and other orders (n = 32), from different areas of Spain. Antibodies to T. gondii (modified agglutination test, MAT titer ≥1∶25) were found in 282 (26.1%, IC95%:23.5–28.7) of the 1079 birds. This study constitute the first extensive survey in wild birds species in Spain and reports for the first time T. gondii antibodies in the griffon vulture (Gyps fulvus), short-toed snake-eagle (Circaetus gallicus), Bonelli's eagle (Aquila fasciata), golden eagle (Aquila chrysaetos), bearded vulture (Gypaetus barbatus), osprey (Pandion haliaetus), Montagu's harrier (Circus pygargus), Western marsh-harrier (Circus aeruginosus), peregrine falcon (Falco peregrinus), long-eared owl (Asio otus), common scops owl (Otus scops), Eurasian spoonbill (Platalea leucorodia), white stork (Ciconia ciconia), grey heron (Ardea cinerea), common moorhen (Gallinula chloropus); in the International Union for Conservation of Nature (IUCN) “vulnerable” Spanish imperial eagle (Aquila adalberti), lesser kestrel (Falco naumanni) and great bustard (Otis tarda); and in the IUCN “near threatened” red kite (Milvus milvus). The highest seropositivity by species was observed in the Eurasian eagle owl (Bubo bubo) (68.1%, 98 of 144). The main risk factors associated with T. gondii seropositivity in wild birds were age and diet, with the highest exposure in older animals and in carnivorous wild birds. The results showed that T. gondii infection is widespread and can be at a high level in many wild birds in Spain, most likely related to their feeding behaviour

The HSP90 Inhibitor NVP-AUY922 Radiosensitizes by Abrogation of Homologous Recombination Resulting in Mitotic Entry with Unresolved DNA Damage

Heat shock protein 90 (HSP90) is a molecular chaperone responsible for the conformational maintenance of a number of client proteins that play key roles in cell cycle arrest, DNA damage repair and apoptosis following radiation. HSP90 inhibitors exhibit antitumor activity by modulating the stabilisation and activation of HSP90 client proteins. We sought to evaluate NVP-AUY922, the most potent HSP90 inhibitor yet reported, in preclinical radiosensitization studies.NVP-AUY922 potently radiosensitized cells in vitro at low nanomolar concentrations with a concurrent depletion of radioresistance-linked client proteins. Radiosensitization by NVP-AUY922 was verified for the first time in vivo in a human head and neck squamous cell carcinoma xenograft model in athymic mice, as measured by delayed tumor growth and increased surrogate end-point survival (p = <0.0001). NVP-AUY922 was shown to ubiquitously inhibit resolution of dsDNA damage repair correlating to delayed Rad51 foci formation in all cell lines tested. Additionally, NVP-AUY922 induced a stalled mitotic phenotype, in a cell line-dependent manner, in HeLa and HN5 cell lines irrespective of radiation exposure. Cell cycle analysis indicated that NVP-AUY922 induced aberrant mitotic entry in all cell lines tested in the presence of radiation-induced DNA damage due to ubiquitous CHK1 depletion, but resultant downstream cell cycle effects were cell line dependent.These results identify NVP-AUY922 as the most potent HSP90-mediated radiosensitizer yet reported in vitro, and for the first time validate it in a clinically relevant in vivo model. Mechanistic analysis at clinically achievable concentrations demonstrated that radiosensitization is mediated by the combinatorial inhibition of cell growth and survival pathways, ubiquitous delay in Rad51-mediated homologous recombination and CHK1-mediated G(2)/M arrest, but that the contribution of cell cycle perturbation to radiosensitization may be cell line specific