Gapped BLAST and PSI-BLAST: a new generation of protein database search programs

The BLAST programs are widely used tools for searching protein and DNA databases for sequence similarities. For protein comparisons, a variety of definitional, algorithmic and statistical refinements described here permits the execution time of the BLAST programs to be decreased substantially while enhancing their sensitivity to weak similarities. A new criterion for triggering the extension of word hits, combined with a new heuristic for generating gapped alignments, yields a gapped BLAST program that runs at approximately three times the speed of the original. In addition, a method is introduced for automatically combining statistically significant alignments produced by BLAST into a position-specific score matrix, and searching the database using this matrix. The resulting Position-Specific Iterated BLAST (PSI-BLAST) program runs at approximately the same speed per iteration as gapped BLAST, but in many cases is much more sensitive to weak but biologically relevant sequence similarities. PSI-BLAST is used to uncover several new and interesting members of the BRCT superfamily.

Rapid colorimetric assay for cellular growth and survival: application to proliferation and cytotoxicity assays

A tetrazolium salt has been used to develop a quantitative colorimetric assay for mammalian cell survival and proliferation. The assay detects living, but not dead cells and the signal generated is dependent on the degree of activation of the cells. This method can therefore be used to measure cytotoxicity, proliferation or activation. The results can be read on a multiwell scanning spectrophotometer (ELISA reader) and show a high degree of precision. No washing steps are used in the assay. The main advantages of the colorimetric assay are its rapidity and precision, and the lack of any radioisotope. We have used the assay to measure proliferative lymphokines, mitogen stimulations and complement-mediated lysis.

Electrophoretic transfer of proteins from polyacrylamide gels to nitrocellulose sheets: procedure and some applications

A method has been devised for the electrophoretic transfer of proteins from polyacrylamide gels to nitrocellulose sheets. The method results in quantitative transfer of ribosomal proteins from gels containing urea. For sodium dodecyl sulfate gels, the original band pattern was obtained with no loss of resolution, but the transfer was not quantitative. The method allows detection of proteins by autoradiography and is simpler than conventional procedures. The immobilized proteins were detectable by immunological procedures. All additional binding capacity on the nitrocellulose was blocked with excess protein; then a specific antibody was bound and, finally, a second antibody directed against the first antibody. The second antibody was either radioactively labeled or conjugated to fluorescein or to peroxidase. The specific protein was then detected by either autoradiography, under UV light, or by the peroxidase reaction product, respectively. In the latter case, as little as 100 pg of protein was clearly detectable. It is anticipated that the procedure will be applicable to analysis of a wide variety of proteins with specific reactions or ligands.

Gene ontology: tool for the unification of biology. The Gene Ontology Consortium

Genomic sequencing has made it clear that a large fraction of the genes specifying the core biological functions are shared by all eukaryotes. Knowledge of the biological role of such shared proteins in one organism can often be transferred to other organisms. The goal of the Gene Ontology Consortium is to produce a dynamic, controlled vocabulary that can be applied to all eukaryotes even as knowledge of gene and protein roles in cells is accumulating and changing. To this end, three independent ontologies accessible on the World-Wide Web (http://www.geneontology.org ) are being constructed: biological process, molecular function and cellular component.

MicroRNAs: genomics, biogenesis, mechanism, and function

MicroRNAs (miRNAs) are endogenous approximately 22 nt RNAs that can play important regulatory roles in animals and plants by targeting mRNAs for cleavage or translational repression. Although they escaped notice until relatively recently, miRNAs comprise one of the more abundant classes of gene regulatory molecules in multicellular organisms and likely influence the output of many protein-coding genes.

A new mathematical model for relative quantification in real-time RT-PCR

Use of the real-time polymerase chain reaction (PCR) to amplify cDNA products reverse transcribed from mRNA is on the way to becoming a routine tool in molecular biology to study low abundance gene expression. Real-time PCR is easy to perform, provides the necessary accuracy and produces reliable as well as rapid quantification results. But accurate quantification of nucleic acids requires a reproducible methodology and an adequate mathematical model for data analysis. This study enters into the particular topics of the relative quantification in real-time RT-PCR of a target gene transcript in comparison to a reference gene transcript. Therefore, a new mathematical model is presented. The relative expression ratio is calculated only from the real-time PCR efficiencies and the crossing point deviation of an unknown sample versus a control. This model needs no calibration curve. Control levels were included in the model to standardise each reaction run with respect to RNA integrity, sample loading and inter-PCR variations. High accuracy and reproducibility (<2.5% variation) were reached in LightCycler PCR using the established mathematical model.

KEGG: kyoto encyclopedia of genes and genomes

KEGG (Kyoto Encyclopedia of Genes and Genomes) is a knowledge base for systematic analysis of gene functions, linking genomic information with higher order functional information. The genomic information is stored in the GENES database, which is a collection of gene catalogs for all the completely sequenced genomes and some partial genomes with up-to-date annotation of gene functions. The higher order functional information is stored in the PATHWAY database, which contains graphical representations of cellular processes, such as metabolism, membrane transport, signal transduction and cell cycle. The PATHWAY database is supplemented by a set of ortholog group tables for the information about conserved subpathways (pathway motifs), which are often encoded by positionally coupled genes on the chromosome and which are especially useful in predicting gene functions. A third database in KEGG is LIGAND for the information about chemical compounds, enzyme molecules and enzymatic reactions. KEGG provides Java graphics tools for browsing genome maps, comparing two genome maps and manipulating expression maps, as well as computational tools for sequence comparison, graph comparison and path computation. The KEGG databases are daily updated and made freely available (http://www. genome.ad.jp/kegg/).

A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences

Some simple formulae were obtained which enable us to estimate evolutionary distances in terms of the number of nucleotide substitutions (and, also, the evolutionary rates when the divergence times are known). In comparing a pair of nucleotide sequences, we distinguish two types of differences; if homologous sites are occupied by different nucleotide bases but both are purines or both pyrimidines, the difference is called type I (or "transition" type), while, if one of the two is a purine and the other is a pyrimidine, the difference is called type II (or "transversion" type). Letting P and Q be respectively the fractions of nucleotide sites showing type I and type II differences between two sequences compared, then the evolutionary distance per site is K = -(1/2) ln [(1-2P-Q) square root of 1-2Q]. The evolutionary rate per year is then given by k = K/(2T), where T is the time since the divergence of the two sequences. If only the third codon positions are compared, the synonymous component of the evolutionary base substitutions per site is estimated by K'S = -(1/2) ln (1-2P-Q). Also, formulae for standard errors were obtained. Some examples were worked out using reported globin sequences to show that synonymous substitutions occur at much higher rates than amino acid-altering substitutions in evolution.

Induction of pluripotent stem cells from mouse embryonic and adult fibroblast cultures by defined factors

Differentiated cells can be reprogrammed to an embryonic-like state by transfer of nuclear contents into oocytes or by fusion with embryonic stem (ES) cells. Little is known about factors that induce this reprogramming. Here, we demonstrate induction of pluripotent stem cells from mouse embryonic or adult fibroblasts by introducing four factors, Oct3/4, Sox2, c-Myc, and Klf4, under ES cell culture conditions. Unexpectedly, Nanog was dispensable. These cells, which we designated iPS (induced pluripotent stem) cells, exhibit the morphology and growth properties of ES cells and express ES cell marker genes. Subcutaneous transplantation of iPS cells into nude mice resulted in tumors containing a variety of tissues from all three germ layers. Following injection into blastocysts, iPS cells contributed to mouse embryonic development. These data demonstrate that pluripotent stem cells can be directly generated from fibroblast cultures by the addition of only a few defined factors.

Isolation of biologically active ribonucleic acid from sources enriched in ribonuclease

Intact ribonucleic acid (RNA) has been prepared from tissues rich in ribonuclease such as the rat pancreas by efficient homogenization in a 4 M solution of the potent protein denaturant guanidinium thiocyanate plus 0.1 M 2-mercaptoethanol to break protein disulfide bonds. The RNA was isolated free of protein by ethanol precipitation or by sedimentation through cesium chloride. Rat pancreas RNA obtained by these means has been used as a source for the purification of alpha-amylase messenger ribonucleic acid.

DADA2: High-resolution sample inference from Illumina amplicon data

We present the open-source software package DADA2 for modeling and correcting Illumina-sequenced amplicon errors (https://github.com/benjjneb/dada2). DADA2 infers sample sequences exactly and resolves differences of as little as 1 nucleotide. In several mock communities, DADA2 identified more real variants and output fewer spurious sequences than other methods. We applied DADA2 to vaginal samples from a cohort of pregnant women, revealing a diversity of previously undetected Lactobacillus crispatus variants.

A technique for radiolabeling DNA restriction endonuclease fragments to high specific activity

A technique for conveniently radiolabeling DNA restriction endonuclease fragments to high specific activity is described. DNA fragments are purified from agarose gels directly by ethanol precipitation and are then denatured and labeled with the large fragment of DNA polymerase I, using random oligonucleotides as primers. Over 70% of the precursor triphosphate is routinely incorporated into complementary DNA, and specific activities of over 10(9) dpm/microgram of DNA can be obtained using relatively small amounts of precursor. These "oligolabeled" DNA fragments serve as efficient probes in filter hybridization experiments.

Initial sequencing and analysis of the human genome

The human genome holds an extraordinary trove of information about human development, physiology, medicine and evolution. Here we report the results of an international collaboration to produce and make freely available a draft sequence of the human genome. We also present an initial analysis of the data, describing some of the insights that can be gleaned from the sequence.

Improved patch-clamp techniques for high-resolution current recording from cells and cell-free membrane patches

1. The extracellular patch clamp method, which first allowed the detection of single channel currents in biological membranes, has been further refined to enable higher current resolution, direct membrane patch potential control, and physical isolation of membrane patches. 2. A description of a convenient method for the fabrication of patch recording pipettes is given together with procedures followed to achieve giga-seals i.e. pipette-membrane seals with resistances of 10(9) - 10(11) omega. 3. The basic patch clamp recording circuit, and designs for improved frequency response are described along with the present limitations in recording the currents from single channels. 4. Procedures for preparation and recording from three representative cell types are given. Some properties of single acetylcholine-activated channels in muscle membrane are described to illustrate the improved current and time resolution achieved with giga-seals. 5. A description is given of the various ways that patches of membrane can be physically isolated from cells. This isolation enables the recording of single channel currents with well-defined solutions on both sides of the membrane. Two types of isolated cell-free patch configurations can be formed: an inside-out patch with its cytoplasmic membrane face exposed to the bath solution, and an outside-out patch with its extracellular membrane face exposed to the bath solution. 6. The application of the method for the recording of ionic currents and internal dialysis of small cells is considered. Single channel resolution can be achieved when recording from whole cells, if the cell diameter is small (less than 20 micrometer). 7. The wide range of cell types amenable to giga-seal formation is discussed.

A pneumonia outbreak associated with a new coronavirus of probable bat origin

Since the outbreak of severe acute respiratory syndrome (SARS) 18 years ago, a large number of SARS-related coronaviruses (SARSr-CoVs) have been discovered in their natural reservoir host, bats1-4. Previous studies have shown that some bat SARSr-CoVs have the potential to infect humans5-7. Here we report the identification and characterization of a new coronavirus (2019-nCoV), which caused an epidemic of acute respiratory syndrome in humans in Wuhan, China. The epidemic, which started on 12 December 2019, had caused 2,794 laboratory-confirmed infections including 80 deaths by 26 January 2020. Full-length genome sequences were obtained from five patients at an early stage of the outbreak. The sequences are almost identical and share 79.6% sequence identity to SARS-CoV. Furthermore, we show that 2019-nCoV is 96% identical at the whole-genome level to a bat coronavirus. Pairwise protein sequence analysis of seven conserved non-structural proteins domains show that this virus belongs to the species of SARSr-CoV. In addition, 2019-nCoV virus isolated from the bronchoalveolar lavage fluid of a critically ill patient could be neutralized by sera from several patients. Notably, we confirmed that 2019-nCoV uses the same cell entry receptor-angiotensin converting enzyme II (ACE2)-as SARS-CoV.

Methods for dietary fiber, neutral detergent fiber, and nonstarch polysaccharides in relation to animal nutrition

There is a need to standardize the NDF procedure. Procedures have varied because of the use of different amylases in attempts to remove starch interference. The original Bacillus subtilis enzyme Type IIIA (XIA) no longer is available and has been replaced by a less effective enzyme. For fiber work, a new enzyme has received AOAC approval and is rapidly displacing other amylases in analytical work. This enzyme is available from Sigma (Number A3306; Sigma Chemical Co., St. Louis, MO). The original publications for NDF and ADF (43, 53) and the Agricultural Handbook 379 (14) are obsolete and of historical interest only. Up to date procedures should be followed. Triethylene glycol has replaced 2-ethoxyethanol because of reported toxicity. Considerable development in regard to fiber methods has occurred over the past 5 yr because of a redefinition of dietary fiber for man and monogastric animals that includes lignin and all polysaccharides resistant to mammalian digestive enzymes. In addition to NDF, new improved methods for total dietary fiber and nonstarch polysaccharides including pectin and beta-glucans now are available. The latter are also of interest in rumen fermentation. Unlike starch, their fermentations are like that of cellulose but faster and yield no lactic acid. Physical and biological properties of carbohydrate fractions are more important than their intrinsic composition.

The neighbor-joining method: a new method for reconstructing phylogenetic trees

A new method called the neighbor-joining method is proposed for reconstructing phylogenetic trees from evolutionary distance data. The principle of this method is to find pairs of operational taxonomic units (OTUs [= neighbors]) that minimize the total branch length at each stage of clustering of OTUs starting with a starlike tree. The branch lengths as well as the topology of a parsimonious tree can quickly be obtained by using this method. Using computer simulation, we studied the efficiency of this method in obtaining the correct unrooted tree in comparison with that of five other tree-making methods: the unweighted pair group method of analysis, Farris's method, Sattath and Tversky's method, Li's method, and Tateno et al.'s modified Farris method. The new, neighbor-joining method and Sattath and Tversky's method are shown to be generally better than the other methods.

Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction

A new method of total RNA isolation by a single extraction with an acid guanidinium thiocyanate-phenol-chloroform mixture is described. The method provides a pure preparation of undegraded RNA in high yield and can be completed within 4 h. It is particularly useful for processing large numbers of samples and for isolation of RNA from minute quantities of cells or tissue samples.

An integrated encyclopedia of DNA elements in the human genome

The human genome encodes the blueprint of life, but the function of the vast majority of its nearly three billion bases is unknown. The Encyclopedia of DNA Elements (ENCODE) project has systematically mapped regions of transcription, transcription factor association, chromatin structure and histone modification. These data enabled us to assign biochemical functions for 80% of the genome, in particular outside of the well-studied protein-coding regions. Many discovered candidate regulatory elements are physically associated with one another and with expressed genes, providing new insights into the mechanisms of gene regulation. The newly identified elements also show a statistical correspondence to sequence variants linked to human disease, and can thereby guide interpretation of this variation. Overall, the project provides new insights into the organization and regulation of our genes and genome, and is an expansive resource of functional annotations for biomedical research.

The ferric reducing ability of plasma (FRAP) as a measure of "antioxidant power": the FRAP assay

A simple, automated test measuring the ferric reducing ability of plasma, the FRAP assay, is presented as a novel method for assessing "antioxidant power." Ferric to ferrous ion reduction at low pH causes a colored ferrous-tripyridyltriazine complex to form. FRAP values are obtained by comparing the absorbance change at 593 nm in test reaction mixtures with those containing ferrous ions in known concentration. Absorbance changes are linear over a wide concentration range with antioxidant mixtures, including plasma, and with solutions containing one antioxidant in purified form. There is no apparent interaction between antioxidants. Measured stoichiometric factors of Trolox, alpha-tocopherol, ascorbic acid, and uric acid are all 2.0; that of bilirubin is 4.0. Activity of albumin is very low. Within- and between-run CVs are <1.0 and <3.0%, respectively, at 100-1000 micromol/liter. FRAP values of fresh plasma of healthy Chinese adults: 612-1634 micromol/liter (mean, 1017; SD, 206; n = 141). The FRAP assay is inexpensive, reagents are simple to prepare, results are highly reproducible, and the procedure is straightforward and speedy. The FRAP assay offers a putative index of antioxidant, or reducing, potential of biological fluids within the technological reach of every laboratory and researcher interested in oxidative stress and its effects.