A plausible mechanism for auxin patterning along the developing root

Background

Recent experimental work has uncovered some of the genetic components required to maintain the Arabidopsis thaliana root stem cell niche (SCN) and its structure. Two main pathways are involved. One pathway depends on the genes SHORTROOT and SCARECROW and the other depends on the PLETHORA genes, which have been proposed to constitute the auxin readouts. Recent evidence suggests that a regulatory circuit, composed of WOX5 and CLE40, also contributes to the SCN maintenance. Yet, we still do not understand how the niche is dynamically maintained and patterned or if the uncovered molecular components are sufficient to recover the observed gene expression configurations that characterize the cell types within the root SCN. Mathematical and computational tools have proven useful in understanding the dynamics of cell differentiation. Hence, to further explore root SCN patterning, we integrated available experimental data into dynamic Gene Regulatory Network (GRN) models and addressed if these are sufficient to attain observed gene expression configurations in the root SCN in a robust and autonomous manner.

Results

We found that an SCN GRN model based only on experimental data did not reproduce the configurations observed within the root SCN. We developed several alternative GRN models that recover these expected stable gene configurations. Such models incorporate a few additional components and interactions in addition to those that have been uncovered. The recovered configurations are stable to perturbations, and the models are able to recover the observed gene expression profiles of almost all the mutants described so far. However, the robustness of the postulated GRNs is not as high as that of other previously studied networks.

Conclusions

These models are the first published approximations for a dynamic mechanism of the A. thaliana root SCN cellular pattering. Our model is useful to formally show that the data now available are not sufficient to fully reproduce root SCN organization and genetic profiles. We then highlight some experimental holes that remain to be studied and postulate some novel gene interactions. Finally, we suggest the existence of a generic dynamical motif that can be involved in both plant and animal SCN maintenance.     

The Generating Function of CID,ETD, and CID/ETD Pairs of Tandem Mass Spectra: Applications to Database Search

Recent emergence of new mass spectrometry techniques (e.g. electron transfer dissociation, ETD) and improved availability of additional proteases (e.g. Lys-N) for protein digestion in high-throughput experiments raised the challenge of designing new algorithms for interpreting the resulting new types of tandem mass (MS/MS) spectra. Traditional MS/MS database search algorithms such as SEQUEST and Mascot were originally designed for collision induced dissociation (CID) of tryptic peptides and are largely based on expert knowledge about fragmentation of tryptic peptides (rather than machine learning techniques) to design CID-specific scoring functions. As a result, the performance of these algorithms is suboptimal for new mass spectrometry technologies or nontryptic peptides. We recently proposed the generating function approach (MS-GF) for CID spectra of tryptic peptides. In this study, we extend MS-GF to automatically derive scoring parameters from a set of annotated MS/MS spectra of any type (e.g. CID, ETD, etc.), and present a new database search tool MS-GFDB based on MS-GF. We show that MS-GFDB outperforms Mascot for ETD spectra or peptides digested with Lys-N. For example, in the case of ETD spectra, the number of tryptic and Lys-N peptides identified by MS-GFDB increased by a factor of 2.7 and 2.6 as compared with Mascot. Moreover, even following a decade of Mascot developments for analyzing CID spectra of tryptic peptides, MS-GFDB (that is not particularly tailored for CID spectra or tryptic peptides) resulted in 28% increase over Mascot in the number of peptide identifications. Finally, we propose a statistical framework for analyzing multiple spectra from the same precursor (e.g. CID/ETD spectral pairs) and assigning p values to peptide-spectrum-spectrum matches.Since the introduction of electron capture dissociation (ECD)1 in 1998 (1), electron-based peptide dissociation technologies have played an important role in analyzing intact proteins and post-translational modifications (2). However, until recently, this research-grade technology was available only to a small number of laboratories because it was commercially unavailable, required experience for operation, and could be implemented only with expensive FT-ICR instruments. The discovery of electron-transfer dissociation (ETD) (3) enabled an ECD-like technology to be implemented in (relatively cheap) ion-trap instruments. Nowadays, many researchers are employing the ETD technology for tandem mass spectra generation (4–9).Although the hardware technologies to generate ETD spectra are maturing rapidly, software technologies to analyze ETD spectra are still in infancy. There are two major approaches to analyzing tandem mass spectra: de novo sequencing and database search. Both approaches find the best-scoring peptide either among all possible peptides (de novo sequencing) or among all peptides in a protein database (database search). Although de novo sequencing is emerging as an alternative to database search, database search remains a more accurate (and thus preferred) method of spectral interpretation, so here we focus on the database search approach.Numerous database search engines are currently available, including SEQUEST (10), Mascot (11), OMSSA (12), X!Tandem (13), and InsPecT (14). However, most of them are inadequate for the analysis of ETD spectra because they are optimized for collision induced dissociation (CID) spectra that show different fragmentation propensities than those of ETD spectra. Additionally, the existing tandem mass spectrometry (MS/MS) tools are biased toward the analysis of tryptic peptides because trypsin is usually used for CID, and thus not suitable for the analysis of nontryptic peptides that are common for ETD. Therefore, even though some database search engines support the analysis of ETD spectra (e.g. SEQUEST, Mascot, and OMSSA), their performance remains suboptimal when it comes to analyzing ETD spectra. Recently, an ETD-specific database search tool (Z-Core) was developed; however it does not significantly improve over OMSSA (15).We present a new database search tool (MS-GFDB) that significantly outperforms existing database search engines in the analysis of ETD spectra, and performs equally well on nontryptic peptides. MS-GFDB employs the generating function approach (MS-GF) that computes rigorous p values of peptide-spectrum matches (PSMs) based on the spectrum-specific score histogram of all peptides (16).² MS-GF p values are dependent only on the PSM (and not on the database), thus can be used as an alternative scoring function for the database search.Computing p values requires a scoring model evaluating qualities of PSMs. MS-GF adopts a probabilistic scoring model (MS-Dictionary scoring model) described in Kim et al., 2009 (17), considering multiple features including product ion types, peak intensities and mass errors. To define the parameters of this scoring model, MS-GF only needs a set of training PSMs.³ This set of PSMs can be obtained in a variety of ways: for example, one can generate CID/ETD pairs and use peptides identified by CID to form PSMs for ETD. Alternatively, one can generate spectra from a purified protein (when PSMs can be inferred from the accurate parent mass alone) or use a previously developed (not necessary optimal) tool to generate training PSMs. From these training PSMs, MS-GF automatically derives scoring parameters without assuming any prior knowledge about the specifics of a particular peptide fragmentation method (e.g. ETD, CID, etc.) and/or proteolytic origin of the peptides. MS-GF was originally designed for the analysis of CID spectra, but now it has been extended to other types of spectra generated by various fragmentation techniques and/or various enzymes. We show that MS-GF can be successfully applied to novel types of spectra (e.g. ETD of Lys-N peptides (18, 19)) by simply retraining scoring parameters without any modification. Note that although the same scoring model is used for different types of spectra, the parameters derived to score different types of spectra are dissimilar.We compared the performance of MS-GFDB with Mascot on a large ETD data set and found that it generated many more peptide identifications for the same false discovery rates (FDR). For example, at 1% peptide level FDR, MS-GFDB identified 9450 unique peptides from 81,864 ETD spectra of Lys-N peptides whereas Mascot only identified 3672 unique peptides, ≈160% increase in the number of peptide identifications (a similar improvement is observed for ETD spectra of tryptic peptides).⁴ MS-GFDB also showed a significant 28% improvement in the number of identified peptides from CID spectra of tryptic peptides (16,203 peptides as compared with 12,658 peptides identified by Mascot).The ETD technology complements rather than replaces CID because both technologies have some advantages: CID for smaller peptides with small charges, ETD for larger and multiply charged peptides (20, 21). An alternative way to utilize ETD is to use it in conjunction with CID because CID and ETD generate complementary sequence information (20, 22, 23). ETD-enabled instruments often support generating both CID and ETD spectra (CID/ETD pairs) for the same peptide. Although the CID/ETD pairs promise a great improvement in peptide identification, the full potential of such pairs has not been fully realized yet. In the case of de novo sequencing, de novo sequencing tools utilizing CID/ETD pairs indeed result in more accurate de novo peptide sequencing than traditional CID-based algorithms (23, 24, 25). However, in the case of database search, the argument that the use of CID/ETD pairs improves peptide identifications remains poorly substantiated. A few tools are developed to use CID/ETD (or CID/ECD) pairs for the database search but they are limited to preprocessing/postprocessing of the spectral data before or following running a traditional database search tool (26, 27). Nielsen et al., 2005 (22) pioneered the combined use of CID and ECD for the database search. Given a CID/ECD pair, they generated a combined spectrum comprised only of complementary pairs of peaks, and searched it with Mascot.⁵ However, this approach is hard to generalize to less accurate CID/ETD pairs generated by ion-trap instruments because there is a higher chance that the identified complementary pairs of peaks are spurious. More importantly, using traditional MS/MS tools (such as Mascot) for the database search of the combined spectrum is inappropriate, because they are not optimized for analyzing such combined spectra; a better approach would be to develop a new database search tool tailored for the combined spectrum. Recently, Molina et al., 2008 (26) studied database search of CID/ETD pairs using Spectrum Mill (Agilent Technologies, Santa Clara, CA) and came to a counterintuitive conclusion that using only CID spectra identifies 12% more unique peptides than using CID/ETD pairs. We believe that it is an acknowledgment of limitations of the traditional MS/MS database search tools for the analysis of multiple spectra generated from a single peptide.In this paper, we modify the generating function approach for interpreting CID/ETD pairs and further apply it to improve the database search with CID/ETD pairs. In contrast to previous approaches, our scoring is specially designed to interpret CID/ETD pairs and can be generalized to analyzing any type of multiple spectra generated from a single peptide. When CID/ETD pairs from trypsin digests are used, MS-GFDB identified 13% and 27% more peptides compared with the case when only CID spectra and only ETD spectra are used, respectively. The difference was even more prominent when CID/ETD pairs from Lys-N digests were used, with 41% and 33% improvement over CID only and ETD only, respectively.Assigning a p value to a PSM greatly helped researchers to evaluate the quality of peptide identifications. We now turn to the problem of assigning a p value to a peptide-spectrum-spectrum match (PS²M) when two spectra in PS²M are generated by different fragmentation technologies (e.g. ETD and CID). We argue that assigning statistical significance to a PS²M (or even PSⁿM) is a prerequisite for rigorous CID/ETD analyses. To our knowledge, MS-GFDB is the first tool to generate statistically rigorous p values of PSⁿMs.The MS-GFDB executable and source code is available at the website of Center for Computational Mass Spectrometry at UCSD (http://proteomics.ucsd.edu). It takes a set of spectra (CID, ETD, or CID/ETD pairs) and a protein database as an input and outputs peptide matches. If the input is a set of CID/ETD pairs, it outputs the best scoring peptide matches and their p values (1) using only CID spectra, (2) using only ETD spectra, and (3) using combined spectra of CID/ETD pairs.     

Innate-Like Control of Human iNKT Cell Autoreactivity via the Hypervariable CDR3β Loop

Invariant Natural Killer T cells (iNKT) are a versatile lymphocyte subset with important roles in both host defense and immunological tolerance. They express a highly conserved TCR which mediates recognition of the non-polymorphic, lipid-binding molecule CD1d. The structure of human iNKT TCRs is unique in that only one of the six complementarity determining region (CDR) loops, CDR3β, is hypervariable. The role of this loop for iNKT biology has been controversial, and it is unresolved whether it contributes to iNKT TCR:CD1d binding or antigen selectivity. On the one hand, the CDR3β loop is dispensable for iNKT TCR binding to CD1d molecules presenting the xenobiotic alpha-galactosylceramide ligand KRN7000, which elicits a strong functional response from mouse and human iNKT cells. However, a role for CDR3β in the recognition of CD1d molecules presenting less potent ligands, such as self-lipids, is suggested by the clonal distribution of iNKT autoreactivity. We demonstrate that the human iNKT repertoire comprises subsets of greatly differing TCR affinity to CD1d, and that these differences relate to their autoreactive functions. These functionally different iNKT subsets segregate in their ability to bind CD1d-tetramers loaded with the partial agonist α-linked glycolipid antigen OCH and structurally different endogenous β-glycosylceramides. Using surface plasmon resonance with recombinant iNKT TCRs and different ligand-CD1d complexes, we demonstrate that the CDR3β sequence strongly impacts on the iNKT TCR affinity to CD1d, independent of the loaded CD1d ligand. Collectively our data reveal a crucial role for CDR3β for the function of human iNKT cells by tuning the overall affinity of the iNKT TCR to CD1d. This mechanism is relatively independent of the bound CD1d ligand and thus forms the basis of an inherent, CDR3β dependent functional hierarchy of human iNKT cells.     

Homing endonuclease mediated gene targeting in Anopheles gambiae cells and embryos

Homing endonuclease genes (HEGs) are ‘selfish’ genetic elements that combine the capability to selectively disrupt specific gene sequences with the ability to rapidly spread from a few individuals to an entire population through homologous recombination repair events. Because of these properties, HEGs are regarded as promising candidates to transfer genetic modifications from engineered laboratory mosquitoes to wild-type populations including Anopheles gambiae the vector of human malaria. Here we show that I-SceI and I-PpoI homing endonucleases cleave their recognition sites with high efficiency in A. gambiae cells and embryos and we demonstrate HEG-induced homologous and non-homologous repair events in a variety of functional assays. We also propose a gene drive system for mosquitoes that is based on our finding that I-PpoI cuts genomic rDNA located on the X chromosome in A. gambiae, which could be used to selectively incapacitate X-carrying spermatozoa thereby imposing a severe male-biased sex ratio.     

Electrophysical characteristics of Azospirillum brasilense Sp245 during interaction with antibodies to various cell surface epitopes

This work was undertaken to examine the electrooptical characteristics of cells of Azospirillum brasilense Sp245 during their interaction with antibodies developed to various cell surface epitopes. We used the dependences of the cell suspension optical density changes induced by electroorientation on the orienting field frequency (740, 1000, 1450, 2000, and 2800kHz). Cell interactions with homologous strain-specific antibodies to the A. brasilense Sp245 O antigen and with homologous antibodies to whole bacterial cells brought about considerable changes in the electrooptical properties of the bacterial suspension. When genus-specific antibodies to the flagellin of the Azospirillum sheathed flagellum and antibodies to the serologically distinct O antigen of A. brasilense Sp7 were included in the A. brasilense Sp245 suspension, the changes caused in the electrooptical signal were slight and had values close to those for the above changes. These findings agree well with the immunochemical characteristics of the Azospirillum O antigens and with the data on the topographical distribution of the Azospirillum major cell surface antigens. The obtained results can serve as a basis for the development of a rapid test for the intraspecies detection of microorganisms.     

Protolytic reactions on reduction of cytochrome c oxidase studied by ATR-FTIR spectroscopy

Reduction of cytochrome c oxidase is coupled to proton uptake, and the reduced-minus-oxidized FTIR spectrum should include signatures of protonation of protolytic centers. The major part of the spectrum shows only small differences between acidic and alkaline conditions, which is consistent with the rather weak pH dependence of the proton uptake stoichiometry. Here we aim at revealing redox state-dependent protonatable sites and present a comprehensive investigation over a wide pH range. The reduced-minus-oxidized transition of cytochrome c oxidase from Paracoccus denitrificans was studied by means of Fourier transform infrared spectroscopy in the pH range 5.2-9.5. Effects of pH were analyzed as the difference between reduced-minus-oxidized FTIR spectra at different pH values. Two pH-dependent processes with apparent pKa values of 6.6 and 8.4 and Hill coefficients 0.9 and 0.1, respectively, were found by this methodology. A sharp OH band appears in the IR "water region" on reduction of the enzyme, independent of pH in the range 6.5-9.0, and downshifted by approximately 940 cm-1 on changing the solvent to D2O and by 10 cm-1 on H216O/H218O isotope exchange. This feature of an asymmetric water molecule may belong to water that is produced in the binuclear center upon reduction or to a structured water molecule that loses a hydrogen bond.     

Distribution of cephalopods in the upper epipelagic northwestern Bering Sea in autumn

In the northwestern Bering Sea in autumn, the epipelagic cephalopod community was represented by the boreal fauna, and was found to be composed of three families and nine species of the order Teuthida: Gonatidae (Berryteuthis magister, Boreoteuthis borealis, Gonatopsis japonicus, Gonatus madokai, Gonatus kamtschaticus, Gonatus onyx, and Gonatus pyros), Chiroteuthidae (Chiroteuthis calyx) and Onychoteuthidae (Onychoteuthis borealijaponica). Two pelagic gonatid species (B. borealis and G. kamtschaticus) dominated the cephalopod community in the upper 50 m. The distribution patterns of B. borealis and G. onyx were associated with diel vertical migrations of these squid. The distribution of two distinct size groups of G. kamtschaticus suggested ontogenetic migration of larger squid to deeper layers, and adds to previous data suggesting that this species may be a heterogeneous assemblage. Demersal B. magister rarely occurred in the surface waters. The occurrence of maturing O. borealijaponica in the southern marine area indicated that these were occasional seasonal migrants from the ocean. The occurrence of juvenile C. calyx suggested that these squid may conduct vertical forage migrations from deep waters to the surface layers.     

Comprehensive proteomic mass spectrometric characterization of human cannabinoid CB2 receptor

The CB1 and CB2 cannabinoid receptors belong to the GPCR superfamily and are associated with a variety of physiological and pathophysiological processes. Both receptors, with several lead compounds at different phases of development, are potentially useful targets for drug discovery. For this reason, fully elucidating the structural features of these membrane-associated proteins would be extremely valuable in designing more selective, novel therapeutic drug molecules. As a first step toward obtaining information on the structural features of the drug-receptor complex, we describe the full mass spectrometric (MS) analysis of the recombinant human cannabinoid CB2 receptor. This first complete proteomic characterization of a GPCR protein beyond rhodopsin was accomplished by a combination of several LC/MS approaches involving nanocapillary liquid chromatography, coupled with either a quadrupole-linear ion trap or linear ion trap-FTICR mass spectrometer. The CB2 receptor, with incorporated N-terminal FLAG and C-terminal HIS6 epitope tags, was functionally expressed in baculovirus cells and purified using a single step of anti-FLAG M2 affinity chromatography. To overcome the difficulties involved with in-gel digestion, due to the highly hydrophobic nature of this membrane-associated protein, we conducted in-solution trypsin and chymotrypsin digestions of purified and desalted samples in the presence of a low concentration of CYMAL5. This was followed by nanoLC peptide separation and analysis using a nanospray ESI source operated in the positive mode. The results can be reported confidently, based on the overlapping sequence data obtained using the highly mass accurate LTQ-FT and the 4000 Q-Trap mass spectrometers. Both instruments gave very similar patterns of identified peptides, with full coverage of all transmembrane helices, resulting in the complete characterization of the cannabinoid CB2 receptor. Mass spectrometric identification of all amino acid residues in the cannabinoid CB2 receptor is a key step toward the "Ligand Based Structural Biology" approach developed in our laboratory for characterizing ligand binding sites in GPCRs using a variety of covalent cannabinergic ligands.     

Phylogenetic diversity in sulphate-reducing bacterial communities from oxidised and reduced bottom sediments of the Barents Sea

In the bottom sediments from a number of the Barents Sea sites, including coastal areas of the Novaya Zemlya, Franz Josef Land, and Svalbard archipelagos, sulphate reduction rates were measured and the phylogenetic composition of sulphate-reducing bacterial (SRB) communities was analysed for the first time. Molecular genetic analysis of the sequences of the 16S rRNA and dsrB genes (the latter encodes the β-subunit of dissimilatory (bi)sulphite reductase) revealed significant differences in the composition of bacterial communities in different sampling stations and sediment horizons of the Barents Sea depending on the physicochemical conditions. The major bacteria involved in reduction of sulphur compounds in Arctic marine bottom sediments belonged to Desulfobulbaceae, Desulfobacteraceae, Desulfovibrionaceae, Desulfuromonadaceae, and Desulfarculaceae families, as well as to uncultured clades SAR324 and Sva0485. Desulfobulbaceae and Desulfuromonadaceae predominated in the oxidised (E_h = 154–226 mV) upper layers of the sediments (up to 9% and 5.9% from all reads of the 16S rRNA gene sequences in the sample, correspondingly), while in deeper, more reduced layers (E_h = ?210 to ?105 mV) the share of Desulfobacteraceae in the SRB community was also significant (up to 5%). The highest relative abundance of members of Desulfarculaceae family (3.1%) was revealed in reduced layers of sandy-clayey sediments from the Barents Sea area affected by currents of transformed (mixed, with changed physicochemical characteristics) Atlantic waters.