Assessing Significance in High-Throughput Experiments by Sequential Goodness of Fit and q-Value Estimation

We developed a new multiple hypothesis testing adjustment called SGoF+ implemented as a sequential goodness of fit metatest which is a modification of a previous algorithm, SGoF, taking advantage of the information of the distribution of p-values in order to fix the rejection region. The new method uses a discriminant rule based on the maximum distance between the uniform distribution of p-values and the observed one, to set the null for a binomial test. This new approach shows a better power/pFDR ratio than SGoF. In fact SGoF+ automatically sets the threshold leading to the maximum power and the minimum false non-discovery rate inside the SGoF' family of algorithms. Additionally, we suggest combining the information provided by SGoF+ with the estimate of the FDR that has been committed when rejecting a given set of nulls. We study different positive false discovery rate, pFDR, estimation methods to combine q-value estimates jointly with the information provided by the SGoF+ method. Simulations suggest that the combination of SGoF+ metatest with the q-value information is an interesting strategy to deal with multiple testing issues. These techniques are provided in the latest version of the SGoF+ software freely available at http://webs.uvigo.es/acraaj/SGoF.htm

Identity-by-descent filtering as a tool for the identification of disease alleles in exome sequence data from distant relatives

Large-scale, deep resequencing may be the next logical step in the genetic investigation of common complex diseases. Because each individual is likely to carry many thousands of variants, the identification of causal alleles requires an efficient strategy to reduce the number of candidate variants. Under many genetic models, causal alleles can be expected to reside within identity-by-descent (IBD) regions shared by affected relatives. In distant relatives, IBD regions constitute a small portion of the genome and can thus greatly reduce the search space for causal alleles. However, the effectiveness of this strategy is unknown. We test the simulated mini-exome data set in extended pedigrees provided by Genetic Analysis Workshop 17. At the fourth- and fifth-degree level of relatedness, case-case pairs shared between 1% and 9% of the genome identical by descent. As expected, no genes were shared identical by descent by all case subjects, but 43 genes were shared by many case subjects across at least 50 replicates. We filtered variants in these genes based on population frequency, function, informativeness, and evidence of association using the family-based association test. This analysis highlighted five genes previously implicated in triglyceride, lipid, and cholesterol metabolism. Comparison with the list of true risk alleles revealed that strict IBD filtering followed by association testing of the rarest alleles was the most sensitive strategy. IBD filtering may be a useful strategy for narrowing down the list of candidate variants in exome data, but the optimal degree of relatedness of affected pairs will depend on the genetic architecture of the disease under study

Targeting vaccinations for the licensed dengue vaccine: considerations for serosurvey design

Background The CYD-TDV vaccine was unusual in that the recommended target population for vaccination was originally defined not only by age, but also by transmission setting as defined by seroprevalence. WHO originally recommended countries consider vaccination against dengue with CYD-TDV vaccine in geographic settings only where prior infection with any dengue serotype, as measured by seroprevalence, was >170% in the target age group. Vaccine was not recommended in settings where seroprevalence was <50%. Test-and-vaccinate strategies suggested following new analysis by Sanofi will still require age-stratified seroprevalence surveys to optimise age-group targeting. Here we address considerations for serosurvey design in the context of vaccination program planning. Methods To explore how the design of seroprevalence surveys affects estimates of transmission intensity, 100 age-specific seroprevalence surveys were simulated using a beta-binomial distribution and a simple catalytic model for different combinations of age-range, survey size, transmission setting, and test sensitivity/specificity. We then used a Metropolis-Hastings Markov Chain Monte-Carlo algorithm to estimate the force of infection from each simulated dataset. Results Sampling from a wide age-range led to more accurate estimates than merely increasing sample size in a narrow age-range. This finding was consistent across all transmission settings. The optimum test sensitivity and specificity given an imperfect test differed by setting with high sensitivity being important in high transmission settings and high specificity important in low transmission settings. Conclusions When assessing vaccination suitability by seroprevalence surveys, countries should ensure an appropriate age-range is sampled, considering epidemiological evidence about the local burden of disease

Relative roles of endothelin-1 and angiotensin II in experimental post-ischaemic acute renal failure

Background. The relative roles of endothelin (ET)-1 and angiotensin (ANG) II in post-ischaemic acute renal failure (ARF) have not been fully established so far. With the aim of contributing to this goal, we assessed in this study the effect of ANG II and ET-1 blockade on the course of post-ischaemic-ARF. Methods. Anaesthetized Wistar rats received i.v. either bosentan (a dual ET receptor antagonist; 10 mg/kg body weight) or losartan [ANG II type 1 (AT(1)) receptor antagonist; 5 or 10 mg/kg body weight] or both, 20 min before, during and 20 min after ischaemia. Rats in the control group received the vehicle via the same route. Survival and renal function were monitored up to 8 days after the ischaemic challenge, while haemodynamic parameters were measured 24 h after ARF. Results. Our results demonstrate that bosentan treatment has a more beneficial effect on experimental ARF than losartan. The survival rate was remarkably higher in bosentan-treated rats than in both rat groups treated with losartan. In the ARF group treated with bosentan, renal blood flow (RBF) was increased by 129% in comparison with the untreated ARF group, whereas in the losartan-treated ARF groups, RBF was only similar to35 or 38% higher than in control ARF rats. The glomerular filtration rate was markedly higher in bosentan-treated rats than in all other ARF groups on the first and second day after ischaemia. Tubular cell injury was less severe in bosentan-treated rats than in the control ARF rats, but in losartan-treated groups it was similar to that in the ARF group. Concurrent blockade of both ET and AT(1) receptors did not improve ARF because this treatment induced a marked decrease in blood pressure. Conclusions. These results suggest that ET-1 blockade is more efficient in improving the early course of post-ischaemic renal injury than ANG II inhibition, and that blockade of ET-1 might be effective in prophylaxis of ischaemic ARF

Population density, water supply, and the risk of dengue fever in Vietnam: cohort study and spatial analysis.

BACKGROUND: Aedes aegypti, the major vector of dengue viruses, often breeds in water storage containers used by households without tap water supply, and occurs in high numbers even in dense urban areas. We analysed the interaction between human population density and lack of tap water as a cause of dengue fever outbreaks with the aim of identifying geographic areas at highest risk. METHODS AND FINDINGS: We conducted an individual-level cohort study in a population of 75,000 geo-referenced households in Vietnam over the course of two epidemics, on the basis of dengue hospital admissions (n = 3,013). We applied space-time scan statistics and mathematical models to confirm the findings. We identified a surprisingly narrow range of critical human population densities between around 3,000 to 7,000 people/km² prone to dengue outbreaks. In the study area, this population density was typical of villages and some peri-urban areas. Scan statistics showed that areas with a high population density or adequate water supply did not experience severe outbreaks. The risk of dengue was higher in rural than in urban areas, largely explained by lack of piped water supply, and in human population densities more often falling within the critical range. Mathematical modeling suggests that simple assumptions regarding area-level vector/host ratios may explain the occurrence of outbreaks. CONCLUSIONS: Rural areas may contribute at least as much to the dissemination of dengue fever as cities. Improving water supply and vector control in areas with a human population density critical for dengue transmission could increase the efficiency of control efforts. Please see later in the article for the Editors' Summary

A specific case in the classification of woods by FTIR and chemometric: discrimination of Fagales from Malpighiales

Fourier transform infrared (FTIR) spectroscopic data was used to classify wood samples from nine species within the Fagales and Malpighiales using a range of multivariate statistical methods. Taxonomic classification of the family Fagaceae and Betulaceae from Angiosperm Phylogenetic System Classification (APG II System) was successfully performed using supervised pattern recognition techniques. A methodology for wood sample discrimination was developed using both sapwood and heartwood samples. Ten and eight biomarkers emerged from the dataset to discriminate order and family, respectively. In the species studied FTIR in combination with multivariate analysis highlighted significant chemical differences in hemicelluloses, cellulose and guaiacyl (lignin) and shows promise as a suitable approach for wood sample classification

From Re-Emergence to Hyperendemicity: The Natural History of the Dengue Epidemic in Brazil

The spread of dengue virus is a major public health problem. Though the burden of dengue has historically been concentrated in Southeast Asian countries, Brazil has become the country that reports the largest number of cases in the world. While prior to 2007 the disease affected mostly adults, during the 2007 epidemic the number of dengue hemorrhagic fever cases more than doubled, and over 53% of cases were in children under 15 years of age. In this paper, we propose that the conditions for the shift were being set gradually since the re-introduction of dengue in 1986 and that they represent the transition from re-emergence to hyperendemicity. Using data from an age stratified seroprevalence study conducted in Recife, we estimated the force of infection (a measure of transmission intensity) between 1986–2006 and used these estimates to simulate the accumulation of immunity since the re-emergence. As the length of time that dengue has circulated increases, adults have a lower probability of remaining susceptible to primary or secondary infection and thus, cases become on average younger. If in fact the shift represents the transition from re-emergence to hyperendemicity, similar shifts are likely to be observed in the rest of Brazil, the American continent and other regions where transmission emerges

LES-based Study of the Roughness Effects on the Wake of a Circular Cylinder from Subcritical to Transcritical Reynolds Numbers

This paper investigates the effects of surface roughness on the flow past a circular cylinder at subcritical to transcritical Reynolds numbers. Large eddy simulations of the flow for sand grain roughness of size k/D = 0.02 are performed (D is the cylinder diameter). Results show that surface roughness triggers the transition to turbulence in the boundary layer at all Reynolds numbers, thus leading to an early separation caused by the increased momentum deficit, especially at transcritical Reynolds numbers. Even at subcritical Reynolds numbers, boundary layer instabilities are triggered in the roughness sublayer and eventually lead to the transition to turbulence. The early separation at transcritical Reynolds numbers leads to a wake topology similar to that of the subcritical regime, resulting in an increased drag coefficient and lower Strouhal number. Turbulent statistics in the wake are also affected by roughness; the Reynolds stresses are larger due to the increased turbulent kinetic energy production in the boundary layer and separated shear layers close to the cylinder shoulders.We acknowledge “Red Española de Surpercomputación” (RES) for awarding us access to the MareNostrum III machine based in Barcelona, Spain (Ref. FI-2015-2-0026 and FI-2015-3-0011). We also acknowledge PRACE for awarding us access to Fermi and Marconi Supercomputers at Cineca, Italy (Ref. 2015133120). Oriol Lehmkuhl acknowledges a PDJ 2014 Grant by AGAUR (Generalitat de Catalunya). Ugo Piomelli acknowledges the support of the Natural Sciences and Engineering Research Council (NSERC) of Canada under the Discovery Grant Programme (Grant No. RGPIN-2016-04391). Ricard Borrell acknowledges a Juan de la Cierva postdoctoral grant (IJCI-2014-21034). Ivette Rodriguez, Oriol Lehmkuhl, Ricard Borrell and Assensi Oliva acknowledge Ministerio de Economía y Competitividad, Secretaría de Estado de Investigación, Desarrollo e Innovación, Spain (ref. ENE2014-60577-R).Peer ReviewedPostprint (author's final draft

The K2K SciBar Detector

A new near detector, SciBar, for the K2K long-baseline neutrino oscillation expe riment was installed to improve the measurement of neutrino energy spectrum and to study neutrino interactions in the energy region around 1 GeV. SciBar is a 'fully active' tracking detector with fine segmentation consisting of plastic scintillator bars. The detector was constructed in summer 2003 and is taking data since October 2003. The basic design and initial performance is presented.Comment: 7 pages, 4figures, Contributed to Proceedings of the 10th Vienna Conference on Instrumentation, Vienna, February 16-21, 200

Targeted 'Next-Generation' sequencing in anophthalmia and microphthalmia patients confirms SOX2, OTX2 and FOXE3 mutations

<p>Abstract</p> <p>Background</p> <p>Anophthalmia/microphthalmia (A/M) is caused by mutations in several different transcription factors, but mutations in each causative gene are relatively rare, emphasizing the need for a testing approach that screens multiple genes simultaneously. We used next-generation sequencing to screen 15 A/M patients for mutations in 9 pathogenic genes to evaluate this technology for screening in A/M.</p> <p>Methods</p> <p>We used a pooled sequencing design, together with custom single nucleotide polymorphism (SNP) calling software. We verified predicted sequence alterations using Sanger sequencing.</p> <p>Results</p> <p>We verified three mutations - c.542delC in S<it>OX2</it>, resulting in p.Pro181Argfs*22, p.Glu105X in <it>OTX2 </it>and p.Cys240X in <it>FOXE3</it>. We found several novel sequence alterations and SNPs that were likely to be non-pathogenic - p.Glu42Lys in <it>CRYBA4</it>, p.Val201Met in <it>FOXE3 </it>and p.Asp291Asn in <it>VSX2</it>. Our analysis methodology gave one false positive result comprising a mutation in <it>PAX6 </it>(c.1268A > T, predicting p.X423LeuextX*15) that was not verified by Sanger sequencing. We also failed to detect one 20 base pair (bp) deletion and one 3 bp duplication in <it>SOX2</it>.</p> <p>Conclusions</p> <p>Our results demonstrated the power of next-generation sequencing with pooled sample groups for the rapid screening of candidate genes for A/M as we were correctly able to identify disease-causing mutations. However, next-generation sequencing was less useful for small, intragenic deletions and duplications. We did not find mutations in 10/15 patients and conclude that there is a need for further gene discovery in A/M.</p