Grid2: A Program for Rapid Estimation of the Jovian Radiation Environment

Grid2 is a program that utilizes the Galileo Interim Radiation Electron model 2 (GIRE2) Jovian radiation model to compute fluences and doses for Jupiter missions. (Note: The iterations of these two softwares have been GIRE and GIRE2; likewise Grid and Grid2.) While GIRE2 is an important improvement over the original GIRE radiation model, the GIRE2 model can take as long as a day or more to compute these quantities for a complete mission. Grid2 fits the results of the detailed GIRE2 code with a set of grids in local time and position thereby greatly speeding up the execution of the model-minutes as opposed to days. The Grid2 model covers the time period from 1971 to 2050 and distances of 1.03 to 30 Jovian diameters (Rj). It is available as a direct-access database through a FORTRAN interface program. The new database is only slightly larger than the original grid version: 1.5 gigabytes (GB) versus 1.2 GB

Plasma-Based Detector of Outer-Space Dust Particles

A report presents a concept for an instrument to be flown in outer space, where it would detect dust particles - especially those associated with comets. The instrument would include a flat plate that would intercept the dust particles. The anticipated spacecraft/dust-particle relative speeds are so high that the impingement of a dust particle on the plate would generate a plasma cloud. Simple electric dipole sensors located equidistantly along the circumference of the plate would detect the dust particle indirectly by detecting the plasma cloud. The location of the dust hit could be estimated from the timing of the detection pulses of the different dipoles. The mass and composition of the dust particle could be estimated from the shapes and durations of the pulses from the dipoles. In comparison with other instruments for detecting hypervelocity dust particles, the proposed instrument offers advantages of robustness, large collection area, and simplicity

Next-generation sequencing of circulating tumor DNA to predict recurrence in triple-negative breast cancer patients with residual disease after neoadjuvant chemotherapy

Next-generation sequencing to detect circulating tumor DNA is a minimally invasive method for tumor genotyping and monitoring therapeutic response. The majority of studies have focused on detecting circulating tumor DNA from patients with metastatic disease. Herein, we tested whether circulating tumor DNA could be used as a biomarker to predict relapse in triple-negative breast cancer patients with residual disease after neoadjuvant chemotherapy. In this study, we analyzed samples from 38 early-stage triple-negative breast cancer patients with matched tumor, blood, and plasma. Extracted DNA underwent library preparation and amplification using the Oncomine Research Panel consisting of 134 cancer genes, followed by high-coverage sequencing and bioinformatics. We detected high-quality somatic mutations from primary tumors in 33 of 38 patients. TP53 mutations were the most prevalent (82%) followed by PIK3CA (16%). Of the 33 patients who had a mutation identified in their primary tumor, we were able to detect circulating tumor DNA mutations in the plasma of four patients (three TP53 mutations, one AKT1 mutation, one CDKN2A mutation). All four patients had recurrence of their disease (100% specificity), but sensitivity was limited to detecting only 4 of 13 patients who clinically relapsed (31% sensitivity). Notably, all four patients had a rapid recurrence (0.3, 4.0, 5.3, and 8.9 months). Patients with detectable circulating tumor DNA had an inferior disease free survival (p < 0.0001; median disease-free survival: 4.6 mos. vs. not reached; hazard ratio = 12.6, 95% confidence interval: 3.06-52.2). Our study shows that next-generation circulating tumor DNA sequencing of triple-negative breast cancer patients with residual disease after neoadjuvant chemotherapy can predict recurrence with high specificity, but moderate sensitivity. For those patients where circulating tumor DNA is detected, recurrence is rapid

Psoriasis prediction from genome-wide SNP profiles

<p>Abstract</p> <p>Background</p> <p>With the availability of large-scale genome-wide association study (GWAS) data, choosing an optimal set of SNPs for disease susceptibility prediction is a challenging task. This study aimed to use single nucleotide polymorphisms (SNPs) to predict psoriasis from searching GWAS data.</p> <p>Methods</p> <p>Totally we had 2,798 samples and 451,724 SNPs. Process for searching a set of SNPs to predict susceptibility for psoriasis consisted of two steps. The first one was to search top 1,000 SNPs with high accuracy for prediction of psoriasis from GWAS dataset. The second one was to search for an optimal SNP subset for predicting psoriasis. The sequential information bottleneck (sIB) method was compared with classical linear discriminant analysis(LDA) for classification performance.</p> <p>Results</p> <p>The best test harmonic mean of sensitivity and specificity for predicting psoriasis by sIB was 0.674(95% CI: 0.650-0.698), while only 0.520(95% CI: 0.472-0.524) was reported for predicting disease by LDA. Our results indicate that the new classifier sIB performs better than LDA in the study.</p> <p>Conclusions</p> <p>The fact that a small set of SNPs can predict disease status with average accuracy of 68% makes it possible to use SNP data for psoriasis prediction.</p

Accelerating Haplotype-Based Genome-Wide Association Study Using Perfect Phylogeny and Phase-Known Reference Data

The genome-wide association study (GWAS) has become a routine approach for mapping disease risk loci with the advent of large-scale genotyping technologies. Multi-allelic haplotype markers can provide superior power compared with single-SNP markers in mapping disease loci. However, the application of haplotype-based analysis to GWAS is usually bottlenecked by prohibitive time cost for haplotype inference, also known as phasing. In this study, we developed an efficient approach to haplotype-based analysis in GWAS. By using a reference panel, our method accelerated the phasing process and reduced the potential bias generated by unrealistic assumptions in phasing process. The haplotype-based approach delivers great power and no type I error inflation for association studies. With only a medium-size reference panel, phasing error in our method is comparable to the genotyping error afforded by commercial genotyping solutions

Combining μXANES and μXRD mapping to analyse the heterogeneity in calcium carbonate granules excreted by the earthworm Lumbricus terrestris

The use of fluorescence full spectral micro-X-ray absorption near-edge structure (μXANES) mapping is becoming more widespread in the hard energy regime. This experimental method using the Ca K-edge combined with micro-X-ray diffraction (μXRD) mapping of the same sample has been enabled on beamline I18 at Diamond Light Source. This combined approach has been used to probe both long- and short-range order in calcium carbonate granules produced by the earthworm Lumbricus terrestris. In granules produced by earthworms cultured in a control artificial soil, calcite and vaterite are observed in the granules. However, granules produced by earthworms cultivated in the same artificial soil amended with 500 p.p.m. Mg also contain an aragonite. The two techniques, μXRD and μXANES, probe different sample volumes but there is good agreement in the phase maps produced

Graph-based analysis of the metabolic exchanges between two co-resident intracellular symbionts, baumannia cicadellinicola and sulcia muelleri with their insect host, homalodisca coagulata

International audienceEndosymbiotic bacteria from different species can live inside cells of the same eukaryotic organism. Metabolic exchanges occur between host and bacteria but also between different endocytobionts. Since a complete genome annotation is available for both, we built the metabolic network of two endosymbiotic bacteria, Sulcia muelleri and Baumannia cicadellinicola, that live inside specific cells of the sharpshooter Homalodisca coagulata and studied the metabolic exchanges involving transfers of carbon atoms between the three. We automatically determined the set of metabolites potentially exogenously acquired (seeds) for both metabolic networks. We show that the number of seeds needed by both bacteria in the carbon metabolism is extremely reduced. Moreover, only three seeds are common to both metabolic networks, indicating that the complementarity of the two metabolisms is not only manifested in the metabolic capabilities of each bacterium, but also by their different use of the same environment. Furthermore, our results show that the carbon metabolism of S. muelleri may be completely independent of the metabolic network of B. cicadellinicola. On the contrary, the carbon metabolism of the latter appears dependent on the metabolism of S. muelleri, at least for two essential amino acids, threonine and lysine. Next, in order to define which subsets of seeds (precursor sets) are sufficient to produce the metabolites involved in a symbiotic function, we used a graph-based method, PITUFO, that we recently developed. Our results highly refine our knowledge about the complementarity between the metabolisms of the two bacteria and their host. We thus indicate seeds that appear obligatory in the synthesis of metabolites are involved in the symbiotic function. Our results suggest both B. cicadellinicola and S. muelleri may be completely independent of the metabolites provided by the co-resident endocytobiont to produce the carbon backbone of the metabolites provided to the symbiotic system (., thr and lys are only exploited by B. cicadellinicola to produce its proteins)

Biomineralisation by earthworms: an investigation into the stability and distribution of amorphous calcium carbonate

Background Many biominerals form from amorphous calcium carbonate (ACC), but this phase is highly unstable when synthesised in its pure form inorganically. Several species of earthworm secrete calcium carbonate granules which contain highly stable ACC. We analysed the milky fluid from which granules form and solid granules for amino acid (by liquid chromatography) and functional group (by Fourier transform infrared (FTIR) spectroscopy) compositions. Granule elemental composition was determined using inductively coupled plasma-optical emission spectroscopy (ICP-OES) and electron microprobe analysis (EMPA). Mass of ACC present in solid granules was quantified using FTIR and compared to granule elemental and amino acid compositions. Bulk analysis of granules was of powdered bulk material. Spatially resolved analysis was of thin sections of granules using synchrotron-based μ-FTIR and EMPA electron microprobe analysis. Results The milky fluid from which granules form is amino acid-rich (≤ 136 ± 3 nmol mg−1 (n = 3; ± std dev) per individual amino acid); the CaCO3 phase present is ACC. Even four years after production, granules contain ACC. No correlation exists between mass of ACC present and granule elemental composition. Granule amino acid concentrations correlate well with ACC content (r ≥ 0.7, p ≤ 0.05) consistent with a role for amino acids (or the proteins they make up) in ACC stabilisation. Intra-granule variation in ACC (RSD = 16%) and amino acid concentration (RSD = 22–35%) was high for granules produced by the same earthworm. Maps of ACC distribution produced using synchrotron-based μ-FTIR mapping of granule thin sections and the relative intensity of the ν2: ν4 peak ratio, cluster analysis and component regression using ACC and calcite standards showed similar spatial distributions of likely ACC-rich and calcite-rich areas. We could not identify organic peaks in the μ-FTIR spectra and thus could not determine whether ACC-rich domains also had relatively high amino acid concentrations. No correlation exists between ACC distribution and elemental concentrations determined by EMPA. Conclusions ACC present in earthworm CaCO3 granules is highly stable. Our results suggest a role for amino acids (or proteins) in this stability. We see no evidence for stabilisation of ACC by incorporation of inorganic components

Clique-Finding for Heterogeneity and Multidimensionality in Biomarker Epidemiology Research: The CHAMBER Algorithm

Commonly-occurring disease etiology may involve complex combinations of genes and exposures resulting in etiologic heterogeneity. We present a computational algorithm that employs clique-finding for heterogeneity and multidimensionality in biomedical and epidemiological research (the "CHAMBER" algorithm).This algorithm uses graph-building to (1) identify genetic variants that influence disease risk and (2) predict individuals at risk for disease based on inherited genotype. We use a set-covering algorithm to identify optimal cliques and a Boolean function that identifies etiologically heterogeneous groups of individuals. We evaluated this approach using simulated case-control genotype-disease associations involving two- and four-gene patterns. The CHAMBER algorithm correctly identified these simulated etiologies. We also used two population-based case-control studies of breast and endometrial cancer in African American and Caucasian women considering data on genotypes involved in steroid hormone metabolism. We identified novel patterns in both cancer sites that involved genes that sulfate or glucuronidate estrogens or catecholestrogens. These associations were consistent with the hypothesized biological functions of these genes. We also identified cliques representing the joint effect of multiple candidate genes in all groups, suggesting the existence of biologically plausible combinations of hormone metabolism genes in both breast and endometrial cancer in both races.The CHAMBER algorithm may have utility in exploring the multifactorial etiology and etiologic heterogeneity in complex disease

Inferring viral quasispecies spectra from 454 pyrosequencing reads

<p>Abstract</p> <p>Background</p> <p>RNA viruses infecting a host usually exist as a set of closely related sequences, referred to as quasispecies. The genomic diversity of viral quasispecies is a subject of great interest, particularly for chronic infections, since it can lead to resistance to existing therapies. High-throughput sequencing is a promising approach to characterizing viral diversity, but unfortunately standard assembly software was originally designed for single genome assembly and cannot be used to simultaneously assemble and estimate the abundance of multiple closely related quasispecies sequences.</p> <p>Results</p> <p>In this paper, we introduce a new <b>Vi</b>ral <b>Sp</b>ectrum <b>A</b>ssembler (ViSpA) method for quasispecies spectrum reconstruction and compare it with the state-of-the-art ShoRAH tool on both simulated and real 454 pyrosequencing shotgun reads from HCV and HIV quasispecies. Experimental results show that ViSpA outperforms ShoRAH on simulated error-free reads, correctly assembling 10 out of 10 quasispecies and 29 sequences out of 40 quasispecies. While ShoRAH has a significant advantage over ViSpA on reads simulated with sequencing errors due to its advanced error correction algorithm, ViSpA is better at assembling the simulated reads after they have been corrected by ShoRAH. ViSpA also outperforms ShoRAH on real 454 reads. Indeed, 7 most frequent sequences reconstructed by ViSpA from a real HCV dataset are viable (do not contain internal stop codons), and the most frequent sequence was within 1% of the actual open reading frame obtained by cloning and Sanger sequencing. In contrast, only one of the sequences reconstructed by ShoRAH is viable. On a real HIV dataset, ShoRAH correctly inferred only 2 quasispecies sequences with at most 4 mismatches whereas ViSpA correctly reconstructed 5 quasispecies with at most 2 mismatches, and 2 out of 5 sequences were inferred without any mismatches. ViSpA source code is available at <url>http://alla.cs.gsu.edu/~software/VISPA/vispa.html</url>.</p> <p>Conclusions</p> <p>ViSpA enables accurate viral quasispecies spectrum reconstruction from 454 pyrosequencing reads. We are currently exploring extensions applicable to the analysis of high-throughput sequencing data from bacterial metagenomic samples and ecological samples of eukaryote populations.</p