Probabilistic alignments with quality scores: an application to short-read mapping toward accurate SNP/indel detection

MOTIVATION: Recent studies have revealed the importance of considering quality scores of reads generated by next-generation sequence (NGS) platforms in various downstream analyses. It is also known that probabilistic alignments based on marginal probabilities (e.g. aligned-column and/or gap probabilities) provide more accurate alignment than conventional maximum score-based alignment. There exists, however, no study about probabilistic alignment that considers quality scores explicitly, although the method is expected to be useful in SNP/indel callers and bisulfite mapping, because accurate estimation of aligned columns or gaps is important in those analyses. RESULTS: In this study, we propose methods of probabilistic alignment that consider quality scores of (one of) the sequences as well as a usual score matrix. The method is based on posterior decoding techniques in which various marginal probabilities are computed from a probabilistic model of alignments with quality scores, and can arbitrarily trade-off sensitivity and positive predictive value (PPV) of prediction (aligned columns and gaps). The method is directly applicable to read mapping (alignment) toward accurate detection of SNPs and indels. Several computational experiments indicated that probabilistic alignments can estimate aligned columns and gaps accurately, compared with other mapping algorithms e.g. SHRiMP2, Stampy, BWA and Novoalign. The study also suggested that our approach yields favorable precision for SNP/indel calling.     

Receptor for advanced glycation end products binds to phosphatidylserine and assists in the clearance of apoptotic cells

Clearance of apoptotic cells is necessary for tissue development, homeostasis and resolution of inflammation. The uptake of apoptotic cells is initiated by an 'eat-me' signal, such as phosphatidylserine, on the cell surface and phagocytes recognize the signal by using specific receptors. In this study, we show that the soluble form of the receptor for advanced glycation end products (RAGE) binds to phosphatidylserine as well as to the apoptotic thymocytes. RAGE-deficient (Rage(-/-)) alveolar macrophages showed impaired phagocytosis of apoptotic thymocytes and defective clearance of apoptotic neutrophils in Rage(-/-) mice. Our results indicate that RAGE functions as a phosphatidylserine receptor and assists in the clearance of apoptotic cells.     

Pdx1 and Ngn3 overexpression enhances pancreatic differentiation of mouse ES cell-derived endoderm population

In order to define the molecular mechanisms regulating the specification and differentiation of pancreatic β-islet cells, we investigated the effect of upregulating Pdx1 and Ngn3 during the differentiation of the β-islet-like cells from murine embryonic stem (ES) cell-derived activin induced-endoderm. Induced overexpression of Pdx1 resulted in a significant upregulation of insulin (Ins1 and Ins2), and other pancreas-related genes. To enhance the developmental progression from the pancreatic bud to the formation of the endocrine lineages, we induced the overexpression express of Ngn3 together with Pdx1. This combination dramatically increased the level and timing of maximal Ins1 mRNA expression to approximately 100% of that found in the βTC6 insulinoma cell line. Insulin protein and C-peptide expression was confirmed by immunohistochemistry staining. These inductive effects were restricted to c-kit(+) endoderm enriched EB-derived populations suggesting that Pdx1/Ngn3 functions after the specification of pancreatic endoderm. Although insulin secretion was stimulated by various insulin secretagogues, these cells had only limited glucose response. Microarray analysis was used to evaluate the expression of a broad spectrum of pancreatic endocrine cell-related genes as well as genes associated with glucose responses. Taken together, these findings demonstrate the utility of manipulating Pdx1 and Ngn3 expression in a stage-specific manner as an important new strategy for the efficient generation of functionally immature insulin-producing β-islet cells from ES cells.     

Frequent retention of heterozygosity for point mutations in p53 and Ikaros in N-ethyl-N-nitrosourea-induced mouse thymic lymphomas

In agreement with Knudson's two-hit theory, recent findings indicate that the inactivation of tumor suppressor genes is not only mediated by the loss of function but also by the dominant-negative or gain-of-function activity. The former generally accompanies loss of a wild-type allele whereas in the latter a wild-type allele is retained. N-Ethyl-N-nitrosourea (ENU), which efficiently induces point mutations, reportedly leads to the development of tumors by activating ras oncogenes. Little is known about how ENU affects tumor suppressor genes and, therefore, we examined ENU-induced mutations of p53 and Ikaros in thymic lymphomas and compared these with mutations of Kras. In addition, loss of heterozygosity was examined for chromosome 11 to which both p53 and Ikaros were mapped. The frequency of point mutations in p53 and Ikaros was 30% (8/27) and 19% (5/27), respectively, comparable to that observed in Kras (33%: 9/27). In total, 14 of the 27 thymic lymphomas examined (52%) harbored mutations in at least one of these genes. One Ikaros mutation was located at the splice donor site, generating a novel splice isoform lacking zinc finger 3, Ik (F3del). Interestingly, 90% (10/11) of the tumors with point mutations retained wild-type alleles of p53 and Ikaros. Sequence analysis revealed that the most common nucleic acid substitutions were T>A (4/8) in p53, T>C (4/5) in Ikaros and G>A/T (8/9) in Kras, suggesting that the spectrum of mutations was gene dependent. These results suggest that point mutations in tumor suppressor genes without loss of the wild-type allele play an important role in ENU-induced lymphomagenesis.     

Plus and minus sexual agglutinins from Chlamydomonas reinhardtii

Gametes of the unicellular green alga Chlamydomonas reinhardtii undergo sexual adhesion via enormous chimeric Hyp-rich glycoproteins (HRGPs), the plus and minus sexual agglutinins, that are displayed on their flagellar membrane surfaces. We have previously purified the agglutinins and analyzed their structural organization using electron microscopy. We report here the cloning and sequencing of the Sag1 and Sad1 genes that encode the two agglutinins and relate their derived amino acid sequences and predicted secondary structure to the morphology of the purified proteins. Both agglutinin proteins are organized into three distinct domains: a head, a shaft in a polyproline II configuration, and an N-terminal domain. The plus and minus heads are related in overall organization but poorly conserved in sequence except for two regions of predicted hydrophobic alpha-helix. The shafts contain numerous repeats of the PPSPX motif previously identified in Gp1, a cell wall HRGP. We propose that the head domains engage in autolectin associations with the distal termini of their own shafts and suggest ways that adhesion may involve head-head interactions, exolectin interactions between the heads and shafts of opposite type, and antiparallel shaft-shaft interactions mediated by carbohydrates displayed in polyproline II configurations.     

Towards a view of functioning dimeric metabotropic receptors

X-ray crystallography was used to solve the atomic structure of the ligand binding domain of the metabotropic glutamate receptor type1 homo-dimer, making it possible to show the conformational change of this domain upon glutamate binding. Studies of dimeric metabotropic receptors thereafter have focused on the respective roles and interaction of the two subunits, on the activation mechanisms following the structural rearrangements of the ligand-binding domain, and on the functional significance of polyvalent cations, the binding of which was identified in the crystal. The direct interaction between the GABA(B) receptor and the metabotropic glutamate receptor (mGluR1) has also attracted attention. Recently, attention has focused on incorporating these structural features into a functional view of the receptors.     

Roles of isoamylase and ADP-glucose pyrophosphorylase in starch granule synthesis in rice endosperm

Amyloplast-targeted green fluorescent protein (GFP) was used to monitor amyloplast division and starch granule synthesis in the developing endosperm of transgenic rice. Two classical starch mutants, sugary and shrunken, contain reduced activities of isoamylase1 (ISA1) and cytosolic ADP-glucose pyrophosphorylase, respectively. Dividing amyloplasts in the wild-type and shrunken endosperms contained starch granules, whereas those in sugary endosperm did not contain detectable granules, suggesting that ISA1 plays a role in granule synthesis at the initiation step. The transition from phytoglycogen to sugary-amylopectin was gradual in the boundary region between the inner and outer endosperms of sugary. These results suggest that the synthesis of sugary-amylopectin and phytoglycogen involved a stochastic process and that ISA1 activity plays a critical role in the stochastic process in starch synthesis in rice endosperm. The reduction of cytosolic ADP-glucose pyrophosphorylase activity in shrunken endosperm did not inhibit granule initiation but severely restrained the subsequent enlargement of granules. The shrunken endosperm often developed pleomorphic amyloplasts containing a large number of underdeveloped granules or a large cluster of small grains of amyloplasts, each containing a simple-type starch granule. Although constriction-type divisions of amyloplasts were much more frequent, budding-type divisions were also found in the shrunken endosperm. We show that monitoring GFP in developing amyloplasts was an effective means of evaluating the roles of enzymes involved in starch granule synthesis in the rice endosperm.