Chances and limits of method restriction: a detailed analysis of suicide methods in Switzerland

The objective of this study was to estimate the potential of method restriction as a public health strategy in suicide prevention. Data from the Swiss Federal Statistical Office and the Swiss Institutes of Forensic Medicine from 2004 were gathered and categorized into suicide submethods according to accessibility to restriction of means. Of suicides in Switzerland, 39.2% are accessible to method restriction. The highest proportions were found in private weapons (13.2%), army weapons (10.4%), and jumps from hot-spots (4.6%). The presented method permits the estimation of the suicide prevention potential of a country by method restriction and the comparison of restriction potentials between suicide methods. In Switzerland, reduction of firearm suicides has the highest potential to reduce the total number of suicides.

16S rRNA terminal restriction fragment length polymorphism for the characterization of the nasopharyngeal microbiota

A novel non-culture based 16S rRNA Terminal Restriction Fragment Length Polymorphism (T-RFLP) method using the restriction enzymes Tsp509I and Hpy166II was developed for the characterization of the nasopharyngeal microbiota and validated using recently published 454 pyrosequencing data. 16S rRNA gene T-RFLP for 153 clinical nasopharyngeal samples from infants with acute otitis media (AOM) revealed 5 Tsp509I and 6 Hpy166II terminal fragments (TFs) with a prevalence of >10%. Cloning and sequencing identified all TFs with a prevalence >6% allowing a sufficient description of bacterial community changes for the most important bacterial taxa. The conjugated 7-valent pneumococcal polysaccharide vaccine (PCV-7) and prior antibiotic exposure had significant effects on the bacterial composition in an additive main effects and multiplicative interaction model (AMMI) in concordance with the 16S rRNA 454 pyrosequencing data. In addition, the presented T-RFLP method is able to discriminate S. pneumoniae from other members of the Mitis group of streptococci, which therefore allows the identification of one of the most important human respiratory tract pathogens. This is usually not achieved by current high throughput sequencing protocols. In conclusion, the presented 16S rRNA gene T-RFLP method is a highly robust, easy to handle and a cheap alternative to the computationally demanding next-generation sequencing analysis. In case a lot of nasopharyngeal samples have to be characterized, it is suggested to first perform 16S rRNA T-RFLP and only use next generation sequencing if the T-RFLP nasopharyngeal patterns differ or show unknown TFs.

Locally Efficient Estimation with Current Status Data and High Dimensional Covariates

In biostatistical applications, interest often focuses on the estimation of the distribution of time T between two consecutive events. If the initial event time is observed and the subsequent event time is only known to be larger or smaller than an observed monitoring time, then the data is described by the well known singly-censored current status model, also known as interval censored data, case I. We extend this current status model by allowing the presence of a time-dependent process, which is partly observed and allowing C to depend on T through the observed part of this time-dependent process. Because of the high dimension of the covariate process, no globally efficient estimators exist with a good practical performance at moderate sample sizes. We follow the approach of Robins and Rotnitzky (1992) by modeling the censoring variable, given the time-variable and the covariate-process, i.e., the missingness process, under the restriction that it satisfied coarsening at random. We propose a generalization of the simple current status estimator of the distribution of T and of smooth functionals of the distribution of T, which is based on an estimate of the missingness. In this estimator the covariates enter only through the estimate of the missingness process. Due to the coarsening at random assumption, the estimator has the interesting property that if we estimate the missingness process more nonparametrically, then we improve its efficiency. We show that by local estimation of an optimal model or optimal function of the covariates for the missingness process, the generalized current status estimator for smooth functionals become locally efficient; meaning it is efficient if the right model or covariate is consistently estimated and it is consistent and asymptotically normal in general. Estimation of the optimal model requires estimation of the conditional distribution of T, given the covariates. Any (prior) knowledge of this conditional distribution can be used at this stage without any risk of losing root-n consistency. We also propose locally efficient one step estimators. Finally, we show some simulation results.

A Combinatorial Approach to the Variable Selection in Multiple Linear Regression: Analysis of Selwood et al Data set -A Case Study.

A combinatorial protocol (CP) is introduced here to interface it with the multiple linear regression (MLR) for variable selection. The efficiency of CP-MLR is primarily based on the restriction of entry of correlated variables to the model development stage. It has been used for the analysis of Selwood et al data set [16], and the obtained models are compared with those reported from GFA [8] and MUSEUM [9] approaches. For this data set CP-MLR could identify three highly independent models (27, 28 and 31) with Q2 value in the range of 0.632-0.518. Also, these models are divergent and unique. Even though, the present study does not share any models with GFA [8], and MUSEUM [9] results, there are several descriptors common to all these studies, including the present one. Also a simulation is carried out on the same data set to explain the model formation in CP-MLR. The results demonstrate that the proposed method should be able to offer solutions to data sets with 50 to 60 descriptors in reasonable time frame. By carefully selecting the inter-parameter correlation cutoff values in CP-MLR one can identify divergent models and handle data sets larger than the present one without involving excessive computer time.

Change in suicide rates in Switzerland before and after firearm restriction resulting from the 2003 “Army XXI” reform

Objective: Firearms are the most common method of suicide among young men in Switzerland. From March 2003 through February 2004, the number of Swiss soldiers was halved as a result of an army reform (Army XXI), leading to a decrease in the availability of guns nationwide. The authors investigated the patterns of the overall suicide rate and the firearm suicide rate before and after the reform. Method: Using a naturalistic study design, the authors compared suicide rates before (1995–2003) and after the intervention (2004–2008) in the affected population (men ages 18–43) and in two comparison groups (women ages 18–44 and men ages 44–53). Data were received from the Swiss Federal Statistical Office. Interrupted time series analysis was used to control for preexisting temporal trends. Alternative methods (Poisson regression, autocorrelation analysis, and surrogate data tests) were used to check validity. Results: The authors found a reduction in both the overall suicide rate and the firearm suicide rate after the Army XXI reform. No significant increases were found for other suicide methods overall. An increase in railway suicides was observed. It was estimated that 22% of the reduction in firearm suicides was substituted by other suicide methods. The attenuation of the suicide rate was not compensated for during the follow-up years. Neither of the comparison groups showed statistically significant changes in firearm suicide rate and overall suicide rate. Conclusions: The restriction of firearm availability in Switzerland resulting from the Army XXI reform was followed by an enduring decrease in the general suicide rate.

A FUNCTIONAL VARIANT OF HLA-DR1 DELINEATED BY T-LYMPHOCYTE CLONE REACTIVITY, PROTEIN CONFORMATION, AND GENOMIC RESTRICTION SITES

This research characterized a serologically indistinguishable form of HLA-DR1 that: (1) cannot stimulate some DR1-restricted or specific T-lymphocyte clones; (2) displays an unusual electrophoretic pattern on two dimensional gels; and (3) is marked by a polymorphic restriction site of the alpha gene. Inefficient stimulation of some DR1-restricted clones was a property of DR1$\sp{+}$ cells that shared HLA-B14 on the same haplotype and/or were carriers of 21-hydroxylase (21-OH) deficiency. Nonclassical 21-OH deficiency frequently demonstrates genetic linkage with HLA-B14;DR1 haplotypes and associates with duplications of C4B and one 21-OH gene. Cells having both stimulatory (DR1$\sb{\rm n}$) and nonstimulatory (DR1$\sb{\rm x}$) parental haplotypes did not mediate proliferation of these clones. However, heterozygous DR1$\sb{\rm x}$, 2 and DR1$\sb{\rm x}$, 7 cells were efficient stimulators of DR2 and DR7 specific clones, respectively, suggesting that a trans acting factor may modify DR1 alleles or products to yield a dominant DR1$\sb{\rm x}$ phenotype. Incompetent stimulator populations did not secrete an intercellular soluble or contact dependent suppressor factor nor did they express interleukin-2 receptors competing for T-cell growth factors. Two dimensional gel analysis of anti-DR immunoprecipitates revealed, in addition to normal DR$\alpha$ and DR$\beta$ chains, a 50kD species from DR1$\sb{\rm x}$ but not from the majority of DR1$\sb{\rm n}$ or non-DR1 cells. The 50kD structure was stable under reducing conditions in SDS and urea, had antigenic homology with DR, and dissociated after boiling into 34kD and 28kD peptide chains apparently identical with DR$\alpha$ and DR$\beta$ as shown by limited digest peptide maps. N-linked glycosylation and sialation of DRgp50 appeared to be unchanged from normal DR$\alpha$ and DR$\beta$. Bg1II digestion and $DR\alpha$ probing of DR1$\sb{\rm x}$ genomic DNA revealed a 4.5kb fragment while DR1$\sb{\rm n}$ DNA yielded 3.8 and 0.76kb fragments; all restriction sites mapped to the 3$\sp\prime$ untranslated region of $DR\alpha$. Collectively, these data suggest that DRgp50 represents a novel combinatorial association between constitutive chains of DR that may interfere with or compete for normal T cell receptor recognition of DR1 as both an alloantigen and restricting element. Furthermore, extensive chromosomal abnormalities previously mapped to the class III region of B14;DR1 haplotypes may extend into the adjacent class II region with consequent intrusion on immune function. ^

Bilateral uterine vessel ligation as a model of intrauterine growth restriction in mice

BACKGROUND Intrauterine growth restriction (IUGR) occurs in up to 10% of pregnancies and is considered as a major risk to develop various diseases in adulthood, such as cardiovascular diseases, insulin resistance, hypertension or end stage kidney disease. Several IUGR models have been developed in order to understand the biological processes linked to fetal growth retardation, most of them being rat or mouse models and nutritional models. In order to reproduce altered placental flow, surgical models have also been developed, and among them bilateral uterine ligation has been frequently used. Nevertheless, this model has never been developed in the mouse, although murine tools display multiple advantages for biological research. The aim of this work was therefore to develop a mouse model of bilateral uterine ligation as a surgical model of IUGR. RESULTS In this report, we describe the set up and experimental data obtained from three different protocols (P1, P2, P3) of bilateral uterine vessel ligation in the mouse. Ligation was either performed at the cervical end of each uterine horn (P1) or at the central part of each uterine horn (P2 and P3). Time of surgery was E16 (P1), E17 (P2) or E16.5 (P3). Mortality, maternal weight and abortion parameters were recorded, as well as placentas weights, fetal resorption, viability, fetal weight and size. Results showed that P1 in test animals led to IUGR but was also accompanied with high mortality rate of mothers (50%), low viability of fetuses (8%) and high resorption rate (25%). P2 and P3 improved most of these parameters (decreased mortality and improved pregnancy outcomes; improved fetal viability to 90% and 27%, respectively) nevertheless P2 was not associated to IUGR contrary to P3. Thus P3 experimental conditions enable IUGR with better pregnancy and fetuses outcomes parameters that allow its use in experimental studies. CONCLUSIONS Our results show that bilateral uterine artery ligation according to the protocol we have developed and validated can be used as a surgical mouse model of IUGR.

Moment Based Inference with Stratified Data

Many datasets used by economists and other social scientists are collected by stratified sampling. The sampling scheme used to collect the data induces a probability distribution on the observed sample that differs from the target or underlying distribution for which inference is to be made. If this effect is not taken into account, subsequent statistical inference can be seriously biased. This paper shows how to do efficient semiparametric inference in moment restriction models when data from the target population is collected by three widely used sampling schemes: variable probability sampling, multinomial sampling, and standard stratified sampling.

Geochemical data from the USGS #1 Portland core, Canon City, Colorado

During the Cretaceous, widespread black shale deposition occurred during a series of Oceanic Anoxic Events (OAEs). Multiple processes are known to control the deposition of marine black shales, including changes in primary productivity, organic matter preservation, and dilution. OAEs offer an opportunity to evaluate the relative roles of these forcing factors. The youngest of these events-the Coniacian to Santonian OAE 3-resulted in a prolonged organic carbon burial event in shallow and restricted marine environments including the Western Interior Seaway. New high-resolution isotope, organic, and trace metal records from the latest Turonian to early Santonian Niobrara Formation are used to characterize the amount and composition of organic matter preserved, as well as the geochemical conditions under which it accumulated. Redox sensitive metals (Mo, Mn, and Re) indicate a gradual drawdown of oxygen leading into the abrupt onset of organic carbon-rich (up to 8%) deposition. High Hydrogen Indices (HI) and organic carbon to total nitrogen ratios (C:N) demonstrate that the elemental composition of preserved marine organic matter is distinct under different redox conditions. Local changes in d13C indicate that redox-controlled early diagenesis can also significantly alter d13Corg records. These results demonstrate that the development of anoxia is of primary importance in triggering the prolonged carbon burial in the Niobrara Formation. Sea level reconstructions, d18O results, and Mo/total organic carbon ratios suggest that stratification and enhanced bottom water restriction caused the drawdown of bottom water oxygen. Increased nutrients from benthic regeneration and/or continental runoff may have sustained primary productivity.

Model simulation output data (CCSM3) related to uplift of African topography and the intensification of upwelling in the Benguela region

The Benguela Current, located off the west coast of southern Africa, is tied to a highly productive upwelling system**1. Over the past 12 million years, the current has cooled, and upwelling has intensified**2, 3, 4. These changes have been variously linked to atmospheric and oceanic changes associated with the glaciation of Antarctica and global cooling**5, the closure of the Central American Seaway**1, 6 or the further restriction of the Indonesian Seaway**3. The upwelling intensification also occurred during a period of substantial uplift of the African continent**7, 8. Here we use a coupled ocean-atmosphere general circulation model to test the effect of African uplift on Benguela upwelling. In our simulations, uplift in the East African Rift system and in southern and southwestern Africa induces an intensification of coastal low-level winds, which leads to increased oceanic upwelling of cool subsurface waters. We compare the effect of African uplift with the simulated impact of the Central American Seaway closure9, Indonesian Throughflow restriction10 and Antarctic glaciation**11, and find that African uplift has at least an equally strong influence as each of the three other factors. We therefore conclude that African uplift was an important factor in driving the cooling and strengthening of the Benguela Current and coastal upwelling during the late Miocene and Pliocene epochs.

Construction of a Z-DNA-specific restriction endonuclease

Novel restriction enzymes can be created by fusing the nuclease domain of FokI endonuclease with defined DNA binding domains. Recently, we have characterized a domain (Zα) from the N-terminal region of human double-stranded RNA adenosine deaminase (hADAR1), which binds the Z-conformation with high specificity. Here we report creation of a conformation-specific endonuclease, Zα nuclease, which is a chimera of Zα and FokI nuclease. Purified Zα nuclease cleaves negatively supercoiled plasmids only when they contain a Z-DNA forming insert, such as (dC-dG)13. The precise location of the cleavage sites was determined by primer extension. Cutting has been mapped to the edge of the B-Z junction, suggesting that Zα nuclease binds within the Z-DNA insert, but cleaves in the nearby B-DNA, by using a mechanism similar to type IIs restriction enzymes. These data show that Zα binds Z-DNA in an environment similar to that in a cell. Zα nuclease, a structure-specific restriction enzyme, may be a useful tool for further study of the biological role of Z-DNA.

Metal ion-mediated substrate-assisted catalysis in type II restriction endonucleases

The 2.15-Å resolution cocrystal structure of EcoRV endonuclease mutant T93A complexed with DNA and Ca2+ ions reveals two divalent metals bound in one of the active sites. One of these metals is ligated through an inner-sphere water molecule to the phosphate group located 3′ to the scissile phosphate. A second inner-sphere water on this metal is positioned approximately in-line for attack on the scissile phosphate. This structure corroborates the observation that the pro-SP phosphoryl oxygen on the adjacent 3′ phosphate cannot be modified without severe loss of catalytic efficiency. The structural equivalence of key groups, conserved in the active sites of EcoRV, EcoRI, PvuII, and BamHI endonucleases, suggests that ligation of a catalytic divalent metal ion to this phosphate may occur in many type II restriction enzymes. Together with previous cocrystal structures, these data allow construction of a detailed model for the pretransition state configuration in EcoRV. This model features three divalent metal ions per active site and invokes assistance in the bond-making step by a conserved lysine, which stabilizes the attacking hydroxide ion nucleophile.

A family of phase-variable restriction enzymes with differing specificities generated by high-frequency gene rearrangements

The hsd genes of Mycoplasma pulmonis encode restriction and modification enzymes exhibiting a high degree of sequence similarity to the type I enzymes of enteric bacteria. The S subunits of type I systems dictate the DNA sequence specificity of the holoenzyme and are required for both the restriction and the modification reactions. The M. pulmonis chromosome has two hsd loci, both of which contain two hsdS genes each and are complex, site-specific DNA inversion systems. Embedded within the coding region of each hsdS gene are a minimum of three sites at which DNA inversions occur to generate extensive amino acid sequence variations in the predicted S subunits. We show that the polymorphic hsdS genes produced by gene rearrangement encode a family of functional S subunits with differing DNA sequence specificities. In addition to creating polymorphisms in hsdS sequences, DNA inversions regulate the phase-variable production of restriction activity because the other genes required for restriction activity (hsdR and hsdM) are expressed only from loci that are oriented appropriately in the chromosome relative to the hsd promoter. These data cast doubt on the prevailing paradigms that restriction systems are either selfish or function to confer protection from invasion by foreign DNA.

GenEST, a powerful bidirectional link between cDNA sequence data and gene expression profiles generated by cDNA-AFLP

The release of vast quantities of DNA sequence data by large-scale genome and expressed sequence tag (EST) projects underlines the necessity for the development of efficient and inexpensive ways to link sequence databases with temporal and spatial expression profiles. Here we demonstrate the power of linking cDNA sequence data (including EST sequences) with transcript profiles revealed by cDNA-AFLP, a highly reproducible differential display method based on restriction enzyme digests and selective amplification under high stringency conditions. We have developed a computer program (GenEST) that predicts the sizes of virtual transcript-derived fragments (TDFs) of in silico-digested cDNA sequences retrieved from databases. The vast majority of the resulting virtual TDFs could be traced back among the thousands of TDFs displayed on cDNA-AFLP gels. Sequencing of the corresponding bands excised from cDNA-AFLP gels revealed no inconsistencies. As a consequence, cDNA sequence databases can be screened very efficiently to identify genes with relevant expression profiles. The other way round, it is possible to switch from cDNA-AFLP gels to sequences in the databases. Using the restriction enzyme recognition sites, the primer extensions and the estimated TDF size as identifiers, the DNA sequence(s) corresponding to a TDF with an interesting expression pattern can be identified. In this paper we show examples in both directions by analyzing the plant parasitic nematode Globodera rostochiensis. Various novel pathogenicity factors were identified by combining ESTs from the infective stage juveniles with expression profiles of ∼4000 genes in five developmental stages produced by cDNA-AFLP.