A probability-conserving cross-section biasing mechanism for variance reduction in Monte Carlo particle transport calculations

In Monte Carlo particle transport codes, it is often important to adjust reaction cross sections to reduce the variance of calculations of relatively rare events, in a technique known as non-analogous Monte Carlo. We present the theory and sample code for a Geant4 process which allows the cross section of a G4VDiscreteProcess to be scaled, while adjusting track weights so as to mitigate the effects of altered primary beam depletion induced by the cross section change. This makes it possible to increase the cross section of nuclear reactions by factors exceeding 10^4 (in appropriate cases), without distorting the results of energy deposition calculations or coincidence rates. The procedure is also valid for bias factors less than unity, which is useful, for example, in problems that involve computation of particle penetration deep into a target, such as occurs in atmospheric showers or in shielding

Connectedness among herds of beef cattle bred under natural service

Background: A procedure to measure connectedness among herds was applied to a beef cattle population bred by natural service. It consists of two steps: (a) computing coefficients of determination (CDs) of comparisons among herds; and (b) building sets of connected herds. Methods: The CDs of comparisons among herds were calculated using a sampling-based method that estimates empirical variances of true and predicted breeding values from a simulated n-sample. Once the CD matrix was estimated, a clustering method that can handle a large number of comparisons was applied to build compact clusters of connected herds of the Bruna dels Pirineus beef cattle. Since in this breed, natural service is predominant and there are almost no links with reference sires, to estimate CDs, an animal model was used taking into consideration all pedigree information and, especially, the connections with dams. A sensitivity analysis was performed to contrast single-trait sire and animal model evaluations with different heritabilities, multiple-trait animal model evaluations with different degrees of genetic correlations and models with maternal effects. Results: Using a sire model, the percentage of connected herds was very low even for highly heritable traits whereas with an animal model, most of the herds of the breed were well connected and high CD values were obtained among them, especially for highly heritable traits (the mean of average CD per herd was 0.535 for a simulated heritability of 0.40). For the lowly heritable traits, the average CD increased from 0.310 in the single-trait evaluation to 0.319 and 0.354 in the multi-trait evaluation with moderate and high genetic correlations, respectively. In models with maternal effects, the average CD per herd for the direct effects was similar to that from single-trait evaluations. For the maternal effects, the average CD per herd increased if the maternal effects had a high genetic correlation with the direct effects, but the percentage of connected herds for maternal effects was very low, less than 12%. Conclusions: The degree of connectedness in a bovine population bred by natural service mating, such as Bruna del Pirineus beef cattle, measured as the CD of comparisons among herds, is high. It is possible to define a pool of animals for which estimated breeding values can be compared after an across-herds genetic evaluation, especially for highly heritable traits

PGA: power calculator for case-control genetic association analyses

<p>Abstract</p> <p>Background</p> <p>Statistical power calculations inform the design and interpretation of genetic association studies, but few programs are tailored to case-control studies of single nucleotide polymorphisms (SNPs) in unrelated subjects.</p> <p>Results</p> <p>We have developed the "Power for Genetic Association analyses" (PGA) package which comprises algorithms and graphical user interfaces for sample size and minimum detectable risk calculations using SNP or haplotype effects under different genetic models and study constrains. The software accounts for linkage disequilibrium and statistical multiple comparisons. The results are presented in graphs or tables and can be printed or exported in standard file formats.</p> <p>Conclusion</p> <p>PGA is user friendly software that can facilitate decision making for association studies of candidate genes, fine-mapping studies, and whole-genome scans. Stand-alone executable files and a Matlab toolbox are available for download at: <url>http://dceg.cancer.gov/bb/tools/pga</url></p

Fast index based algorithms and software for matching position specific scoring matrices

BACKGROUND: In biological sequence analysis, position specific scoring matrices (PSSMs) are widely used to represent sequence motifs in nucleotide as well as amino acid sequences. Searching with PSSMs in complete genomes or large sequence databases is a common, but computationally expensive task. RESULTS: We present a new non-heuristic algorithm, called ESAsearch, to efficiently find matches of PSSMs in large databases. Our approach preprocesses the search space, e.g., a complete genome or a set of protein sequences, and builds an enhanced suffix array that is stored on file. This allows the searching of a database with a PSSM in sublinear expected time. Since ESAsearch benefits from small alphabets, we present a variant operating on sequences recoded according to a reduced alphabet. We also address the problem of non-comparable PSSM-scores by developing a method which allows the efficient computation of a matrix similarity threshold for a PSSM, given an E-value or a p-value. Our method is based on dynamic programming and, in contrast to other methods, it employs lazy evaluation of the dynamic programming matrix. We evaluated algorithm ESAsearch with nucleotide PSSMs and with amino acid PSSMs. Compared to the best previous methods, ESAsearch shows speedups of a factor between 17 and 275 for nucleotide PSSMs, and speedups up to factor 1.8 for amino acid PSSMs. Comparisons with the most widely used programs even show speedups by a factor of at least 3.8. Alphabet reduction yields an additional speedup factor of 2 on amino acid sequences compared to results achieved with the 20 symbol standard alphabet. The lazy evaluation method is also much faster than previous methods, with speedups of a factor between 3 and 330. CONCLUSION: Our analysis of ESAsearch reveals sublinear runtime in the expected case, and linear runtime in the worst case for sequences not shorter than | [Formula: see text] |(m )+ m - 1, where m is the length of the PSSM and [Formula: see text] a finite alphabet. In practice, ESAsearch shows superior performance over the most widely used programs, especially for DNA sequences. The new algorithm for accurate on-the-fly calculations of thresholds has the potential to replace formerly used approximation approaches. Beyond the algorithmic contributions, we provide a robust, well documented, and easy to use software package, implementing the ideas and algorithms presented in this manuscript

Epistatic Module Detection for Case-Control Studies: A Bayesian Model with a Gibbs Sampling Strategy

The detection of epistatic interactive effects of multiple genetic variants on the susceptibility of human complex diseases is a great challenge in genome-wide association studies (GWAS). Although methods have been proposed to identify such interactions, the lack of an explicit definition of epistatic effects, together with computational difficulties, makes the development of new methods indispensable. In this paper, we introduce epistatic modules to describe epistatic interactive effects of multiple loci on diseases. On the basis of this notion, we put forward a Bayesian marker partition model to explain observed case-control data, and we develop a Gibbs sampling strategy to facilitate the detection of epistatic modules. Comparisons of the proposed approach with three existing methods on seven simulated disease models demonstrate the superior performance of our approach. When applied to a genome-wide case-control data set for Age-related Macular Degeneration (AMD), the proposed approach successfully identifies two known susceptible loci and suggests that a combination of two other loci—one in the gene SGCD and the other in SCAPER—is associated with the disease. Further functional analysis supports the speculation that the interaction of these two genetic variants may be responsible for the susceptibility of AMD. When applied to a genome-wide case-control data set for Parkinson's disease, the proposed method identifies seven suspicious loci that may contribute independently to the disease

Enhanced genetic maps from family-based disease studies: population-specific comparisons

Abstract Background Accurate genetic maps are required for successful and efficient linkage mapping of disease genes. However, most available genome-wide genetic maps were built using only small collections of pedigrees, and therefore have large sampling errors. A large set of genetic studies genotyped by the NHLBI Mammalian Genotyping Service (MGS) provide appropriate data for generating more accurate maps. Results We collected a large sample of uncleaned genotype data for 461 markers generated by the MGS using the Weber screening sets 9 and 10. This collection includes genotypes for over 4,400 pedigrees containing over 17,000 genotyped individuals from different populations. We identified and cleaned numerous relationship and genotyping errors, as well as verified the marker orders. We used this dataset to test for population-specific genetic maps, and to re-estimate the genetic map distances with greater precision; standard errors for all intervals are provided. The map-interval sizes from the European (or European descent), Chinese, and Hispanic samples are in quite good agreement with each other. We found one map interval on chromosome 8p with a statistically significant size difference between the European and Chinese samples, and several map intervals with significant size differences between the African American and Chinese samples. When comparing Palauan with European samples, a statistically significant difference was detected at the telomeric region of chromosome 11p. Several significant differences were also identified between populations in chromosomal and genome lengths. Conclusions Our new population-specific screening set maps can be used to improve the accuracy of disease-mapping studies. As a result of the large sample size, the average length of the 95% confidence interval (CI) for a 10 cM map interval is only 2.4 cM, which is considerably smaller than on previously published maps.http://deepblue.lib.umich.edu/bitstream/2027.42/112826/1/12881_2010_Article_748.pd

Parametric design optimisation of proximal humerus plates based on finite element method

Optimal treatment of proximal humerus fractures remains controversial. Locking plates offer theoretical advantages but are associated with complications in the clinic. This study aimed to perform parametric design optimisation of proximal humerus plates to enhance their mechanical performance. A finite element (FE) model was developed that simulated a two-part proximal humerus fracture that had been treated with a Spatial Subchondral Support (S3) plate and subjected to varus bending. The FE model was validated against in vitro biomechanical test results. The predicted load required to apply 5 mm cantilever varus bending was only 0.728% lower. The FE model was then used to conduct a parametric optimisation study to determine the orientations of inferomedial plate screws that would yield minimum fracture gap change (i.e. optimal stability). The feasible design space was automatically identified by imposing clinically relevant constraints, and the creation process of each FE model for the design optimisation was automated. Consequently, 538 FE models were generated, from which the obtained optimal model had 4.686% lower fracture gap change (0.156 mm) than that of the manufacturer’s standard plate. Whereas its screws were oriented towards the inferomedial region and within the range of neck-shaft angle of a healthy subject. The methodology presented in this study promises future applications in patient-specific design optimisation of implants for other regions of the human body

Spiders do not escape reproductive manipulations by Wolbachia

<p>Abstract</p> <p>Background</p> <p>Maternally inherited bacteria that reside obligatorily or facultatively in arthropods can increase their prevalence in the population by altering their hosts' reproduction. Such reproductive manipulations have been reported from the major arthropod groups such as insects (in particular hymenopterans, butterflies, dipterans and beetles), crustaceans (isopods) and mites. Despite the observation that endosymbiont bacteria are frequently encountered in spiders and that the sex ratio of particular spider species is strongly female biased, a direct relationship between bacterial infection and sex ratio variation has not yet been demonstrated for this arthropod order.</p> <p>Results</p> <p>Females of the dwarf spider <it>Oedothorax gibbosus </it>exhibit considerable variation in the sex ratio of their clutches and were infected with at least three different endosymbiont bacteria capable of altering host reproduction i.e. <it>Wolbachia</it>, <it>Rickettsia </it>and <it>Cardinium</it>. Breeding experiments show that sex ratio variation in this species is primarily maternally inherited and that removal of the bacteria by antibiotics restores an unbiased sex ratio. Moreover, clutches of females infected with <it>Wolbachia </it>were significantly female biased while uninfected females showed an even sex ratio. As female biased clutches were of significantly smaller size compared to non-distorted clutches, killing of male embryos appears to be the most likely manipulative effect.</p> <p>Conclusions</p> <p>This represents to our knowledge the first direct evidence that endosymbiont bacteria, and in particular <it>Wolbachia</it>, might induce sex ratio variation in spiders. These findings are pivotal to further understand the diversity of reproductive phenotypes observed in this arthropod order.</p

Mating dynamics in a nematode with three sexes and its evolutionary implications

Nematodes have diverse reproductive strategies, which make them ideal subjects for comparative studies to address how mating systems evolve. Here we present the sex ratios and mating dynamics of the free-living nematode Rhabditis sp. SB347, in which males, females and hermaphrodites co-exist. The three sexes are produced by both selfing and outcrossing, and females tend to appear early in a mother’s progeny. Males prefer mating with females over hermaphrodites, which our results suggest is related to the female-specific production of the sex pheromones ascr#1 and ascr#9. We discuss the parallels between this system and that of parasitic nematodes that exhibit alternation between uniparental and biparental reproduction

Rare Variants in Ischemic Stroke: An Exome Pilot Study

The genetic architecture of ischemic stroke is complex and is likely to include rare or low frequency variants with high penetrance and large effect sizes. Such variants are likely to provide important insights into disease pathogenesis compared to common variants with small effect sizes. Because a significant portion of human functional variation may derive from the protein-coding portion of genes we undertook a pilot study to identify variation across the human exome (i.e., the coding exons across the entire human genome) in 10 ischemic stroke cases. Our efforts focused on evaluating the feasibility and identifying the difficulties in this type of research as it applies to ischemic stroke. The cases included 8 African-Americans and 2 Caucasians selected on the basis of similar stroke subtypes and by implementing a case selection algorithm that emphasized the genetic contribution of stroke risk. Following construction of paired-end sequencing libraries, all predicted human exons in each sample were captured and sequenced. Sequencing generated an average of 25.5 million read pairs (75 bp×2) and 3.8 Gbp per sample. After passing quality filters, screening the exomes against dbSNP demonstrated an average of 2839 novel SNPs among African-Americans and 1105 among Caucasians. In an aggregate analysis, 48 genes were identified to have at least one rare variant across all stroke cases. One gene, CSN3, identified by screening our prior GWAS results in conjunction with our exome results, was found to contain an interesting coding polymorphism as well as containing excess rare variation as compared with the other genes evaluated. In conclusion, while rare coding variants may predispose to the risk of ischemic stroke, this fact has yet to be definitively proven. Our study demonstrates the complexities of such research and highlights that while exome data can be obtained, the optimal analytical methods have yet to be determined