Tests for the statistical significance of protein sequence similarities in data-bank searches

A suite of tests to evaluate the statistical significance of protein sequence similarities is developed for use in data bank searches. The tests are based on the Wilbur-Lipman word-search algorithm, and take into account the sequence lengths and compositions, and optionally the weighting of amino acid matches. The method is extended to allow for the existence of a sequence insertion/deletion within the region of similarity. The accuracy of statistical distributions underlying the tests is validated using randomly generated sequences and real sequences selected at random from the data banks. A computer program to perform the tests is briefly described.     

Generalisation of New Sequence Knowledge Depends on Response Modality

New visuomotor skills can guide behaviour in novel situations. Prior studies indicate that learning a visuospatial sequence via responses based on manual key presses leads to effector- and response-independent knowledge. Little is known, however, about the extent to which new sequence knowledge can generalise, and, thereby guide behaviour, outside of the manual response modality. Here, we examined whether learning a visuospatial sequence either via manual (key presses, without eye movements), oculomotor (obligatory eye movements), or perceptual (covert reorienting of visuospatial attention) responses supported generalisation to direct and indirect tests administered either in the same (baseline conditions) or a novel response modality (transfer conditions) with respect to initial study. Direct tests measured the use of conscious knowledge about the studied sequence, whereas the indirect tests did not ostensibly draw on the study phase and measured response priming. Oculomotor learning supported the use of conscious knowledge on the manual direct tests, whereas manual learning supported generalisation to the oculomotor direct tests but did not support the conscious use of knowledge. Sequence knowledge acquired via perceptual responses did not generalise onto any of the manual tests. Manual, oculomotor, and perceptual sequence learning all supported generalisation in the baseline conditions. Notably, the manual baseline condition and the manual to oculomotor transfer condition differed in the magnitude of general skill acquired during the study phase; however, general skill did not predict performance on the post-study tests. The results demonstrated that generalisation was only affected by the responses used to initially code the visuospatial sequence when new knowledge was applied to a novel response modality. We interpret these results in terms of response-effect distinctiveness, the availability of integrated effector- and motor-plan based information, and discuss their implications for neurocognitive accounts of sequence learning.     

Lateglacial and Holocene vegetation history in the mountain range of central Pyrenees

We present the pollen analysis of a new sedimentary sequence taken at La Pouretère ( 1720 m), in the mountain vegetation zone of the Marcadau valley (central Pyrenees). The Lateglacial and Holocene chronology is supported by six 14C-dating results. The complementary analysis of some vegetal macroremains, stomata, pollen-clusters and the use of pollen influx allows us to elucidate the dynamic of mountain species such as Pinus and specially Abies but also to infer the unusual part played by Betula at the beginning of the Postglacial period.     

Birds perching on bushes: Networks to visualize conflicting phylogenetic signals during early avian radiation

Hybridization is increasingly seen as an important source of adaptive genetic variation and biotic diversity. Recent phylogenetic studies on the early evolution of birds suggest that the early diversification of neoavian orders perhaps involved a period of extensive hybridization or incomplete lineage sorting. Phylogenetic error, saturation, long-branch attraction, and convergence make it difficult to detect ancient hybridization events and differentiate them from incomplete lineage sorting using sequence data. We used recently published retroposon marker data to visualize the early radiation of Neoaves within a phylogenetic network approach, and found that the most basal neoavian taxa indeed show a complex pattern of reticulated relationships. Moreover, the reticulation levels of different parts of the network are consistent with the insertion pattern of the retroposon elements. The use of network-based analyses on homoplasy-free data shows true conflicting signals and the taxa involved that are not represented in trees.     

The amino acid sequence of thiogalactoside transacetylase of Escherichia coli

The amino acid sequence of thiogalactoside transacetylase, a dimer, has been determined. The monomer contains 202 amino acid residues in a single polypeptide chain and has a molecular weight of 22,671. The analysis was carried out by treatment of the carboxymethylated protein with cyanogen bromide and with trypsin. All seven cyanogen bromide peptides were isolated in pure form and were ordered by peptides isolated from tryptic digests. The sequence analysis was aided by determination of the DNA sequence of the lacA gene. The amino terminus of the protein is heterogenous because the initiator methionine is only partially cleaved. Another rather unusual feature of this cytoplasmic protein is a very hydrophobic segment in the center portion of the chain. Comparison of the amino acid sequence of thiogalactoside transacetylase to those of the lac repressor, beta-galactosidase, and lactose permease did not reveal any marked similarities. Therefore, there is no obvious evolutionary relatedness among proteins of the Lactose Operon.     

Extraction of symbolic determinants common to a family of biological sequences

A set of sequences can be defined by their common subsequences, and the length of these is a measure of the overall resemblance of the set. Each subsequence corresponds to a succession of symbols embedded in every sequence, following the same order but not necessarily contiguous. Determining the longest common subsequence (LCS) requires the exhaustive testing of all possible common subsequences, which sum up to about 2L, if L is the length of the shortest sequence. We present a polynomial algorithm (O(n X L4), where n is the number of sequences) for generating strings related to the LCS and constructed with the sequence alphabet and an indetermination symbol. Such strings are iteratively improved by deleting indetermination symbols and concomitantly introducing the greatest number of alphabet symbols. Processed accordingly, nucleic acid and protein sequences lead to key-words encompassing the salient positions of homologous chains, which can be used for aligning or classifying them, as well as for finding related sequences in data banks.     

A permutation procedure to correct for confounders in case-control studies, including tests of rare variation

Many case-control tests of rare variation are implemented in statistical frameworks that make correction for confounders like population stratification difficult. Simple permutation of disease status is unacceptable for resolving this issue because the replicate data sets do not have the same confounding as the original data set. These limitations make it difficult to apply rare-variant tests to samples in which confounding most likely exists, e.g., samples collected from admixed populations. To enable the use of such rare-variant methods in structured samples, as well as to facilitate permutation tests for any situation in which case-control tests require adjustment for confounding covariates, we propose to establish the significance of a rare-variant test via a modified permutation procedure. Our procedure uses Fisher's noncentral hypergeometric distribution to generate permuted data sets with the same structure present in the actual data set such that inference is valid in the presence of confounding factors. We use simulated sequence data based on coalescent models to show that our permutation strategy corrects for confounding due to population stratification that, if ignored, would otherwise inflate the size of a rare-variant test. We further illustrate the approach by using sequence data from the Dallas Heart Study of energy metabolism traits. Researchers can implement our permutation approach by using the R package BiasedUrn.     

Graphical coding of nucleic acid sequences

When, in a nucleic acid sequence, the four letters C, G, A, T (or U) are replaced by suitable graphical symbols, some patterns become immediately apparent. Two sets of symbols, constructed for the analysis of either purine/pyrimidine alternations, or of regions of complementarity within a sequence are shown. In addition, another mode of coding is presented, in which the four letters are represented by vectors. The sequence is thus transformed into a planar trajectory. We show, in the case of the gene for human beta hemoglobin, that such a coding enables an easy discrimination between introns and exons.     

Nonparametric Evaluation of Repeated Measurement Designs by Point-Symmetry Testing

It is suggested to use simultaneous point-symmetry tests in evaluating trivariate repeated measurement designs nonparametrically. With k =3 treatments being applied to each of N individuals in randomized sequence, pairs of treatments are compared as to 3 variables or incremental sign patterns. It is shown that unifold trivariate tests preserve more information than do multifold univariate tests. The method is illustrated by an example from psychophysiological research.     

Applications of binary segmentation to the estimation of quantal response curves and spatial intensity

This paper explores the use of binary segmentation procedures in two applications. The first application is concerned with the estimation of nonparametric quantal response curves. With Bernoulli data and an assumed monotone increasing curve, this gives rise a change-point model where the change points are determined using a sequence of nested hypothesis tests of whether a change point exists. The second application concerns cluster identification and inference for spatial data where the shape of the clusters and the number of clusters is unknown. The procedure involves a sequence of nested hypothesis tests of a single cluster versus a pair of distinct clusters. Examples of both applications are provided.     

Two-sample tests for comparing intra-individual genetic sequence diversity between populations

Consider a study of two groups of individuals infected with a population of a genetically related heterogeneous mixture of viruses, and multiple viral sequences are sampled from each person. Based on estimates of genetic distances between pairs of aligned viral sequences within individuals, we develop four new tests to compare intra-individual genetic sequence diversity between the two groups. This problem is complicated by two levels of dependency in the data structure: (i) Within an individual, any pairwise distances that share a common sequence are positively correlated; and (ii) for any two pairings of individuals which share a person, the two differences in intra-individual distances between the paired individuals are positively correlated. The first proposed test is based on the difference in mean intra-individual pairwise distances pooled over all individuals in each group, standardized by a variance estimate that corrects for the correlation structure using U-statistic theory. The second procedure is a nonparametric rank-based analog of the first test, and the third test contrasts the set of subject-specific average intra-individual pairwise distances between the groups. These tests are very easy to use and solve correlation problem (i). The fourth procedure is based on a linear combination of all possible U-statistics calculated on independent, identically distributed sequence subdatasets, over the two levels (i) and (ii) of dependencies in the data, and is more complicated than the other tests but can be more powerful. Although the proposed methods are empirical and do not fully utilize knowledge from population genetics, the tests reflect biology through the evolutionary models used to derive the pairwise sequence distances. The new tests are evaluated theoretically and in a simulation study, and are applied to a dataset of 200 HIV sequences sampled from 21 children.     

Pear scab resistance QTLs via a European pear (Pyrus communis) linkage map

Pear scab caused by Venturia pyrina is an economically important disease throughout the world and can cause severe crop loss in susceptible cultivars. The varying range of susceptibility to pear scab in F1 populations has made it possible to identify quantitative trait loci (QTLs). Ninety-five seedlings derived from the cross ‘Abbè Fétel’ (AF) × ‘Max Red Bartlett’ (MRB) were evaluated for scab resistance in greenhouse tests, with 39% being classified as resistant, 33 as moderately susceptible and 28 as highly susceptible. Amplified fragment length polymorphisms (157) and simple sequence repeats (41) were used to construct two maps, one of 908.1 cM (AF) and the other of 879.8 cM (MRB). The analysis of the resistance data collected made it possible to identify two major QTLs on linkage groups 3 and 7 associated with resistance to V. pyrina. Both QTLs explained 88% of the phenotypic variance and the log of odds values were higher than 10, suggesting the involvement of two major genes in pear scab resistance. L. Pierantoni and L. Dondini have contributed equally to this work.     

DNA tests of neutral theory: applications in marine genetics

The principal methods of using DNA sequence information to test the neutral theory of evolution and polymorphism are described. These include the use of synonymous and nonsynonymous substitutions for detecting purifying and positive selection, the analysis of nucleotide diversity, mismatch analysis and the HKA, McDonald-Kreitman, Tajima and Ewens-Watterson tests. Analysis of the covariation of different kinds of molecular markers and the relationship between genetic variation and fitness is also considered. Examples of the use of these approaches in a wide variety of marine organisms are described. It is emphasised that tests of neutral theory, in addition to providing important fundamental knowledge about the action of evolutionary forces, provide valuable information about the influence of environmental and demographic factors.     

ISSR标记技术及其在遗传多样性研究中的应用

ISSR(Inter-Simple Sequence Repeat)技术是在PCR中直接使用微卫星序列进行DNA扩增的一种DNA分子标记。章主要介绍了ISSR标记的原理、方法、特点及其在遗传多样性研究中的应用。ISSR标记方法具有无需知道任何靶标序列的微卫星背景信息、遗传多态性高、检测快速等特点,在遗传多样性研究中具有广泛的应用前景。     

Recent aspects of genetic manipulation in Bacillus thuringiensis

The conjugative plasmid pAM beta 1 was transferred from Streptococcus faecalis to several strains of Bacillus thuringiensis by a filter-mating process. From a transconjugant clone of B. thuringiensis a hybrid plasmid resulting from an in vivo insertion into pAM beta 1 of a 3 Md DNA sequence was isolated. This 3 Md DNA molecule (Th sequence) is related to several host plasmids found in different serotypes of B. thuringiensis. A reciprocal conjugation-like process involving the transfer of pAM beta 1 from B. thuringiensis to S. faecalis was also demonstrated. The comparison of the restriction maps of the crystal genes from plasmid and chromosomal origins of different serotypes, six of which having been cloned in E. coli, revealed the existence of two classes of genes which are very similar in the map corresponding to the N-terminal part of the protein, and which differ essentially in the 3' region. The presence of the transposon-like Th sequence was found in several cases associated with the crystal gene in the same host plasmid, and a model for their structural organization is proposed.     

Scoring functions in protein folding and design

We present an analysis of the assumptions behind some of the most commonly used methods for evaluating the goodness of the fit between a sequence and a structure. Our studies on a lattice model show that methods based on statistical considerations are easy to use and can capture some of the features of protein-like sequences and their corresponding native states, but unfortunately are incapable of recognizing, with certainty, the native-like conformation of a sequence among a set of decoys. Meanwhile, an optimization method, entailing the determination of the parameters of an effective free energy of interaction, is much more reliable in recognizing the native state of a sequence. However, the statistical method is shown to perform quite well in tests of protein design.     

Morphological correlates of the first metacarpal proximal articular surface with manipulative capabilities in apes,humans and South African early hominins

This study quantifies the metacarpal 1 (MC 1) proximal articular surface using three-dimensional morphometrics in extant hominids and fossil hominins (SK 84, cf. Paranthropus robustus/Homo erectus and StW 418, Australopithecus africanus) to understand which characteristics of the proximal metacarpal 1 are potentially correlated with human manipulative abilities and if they can be used in a paleoanthropological setting. A principal components (PC) analysis was used to compare MC 1 proximal articular surface shape and ANOVA and Tukey's HSD post-hoc tests were conducted to determine differences among groups. Homo is significantly different from nonhuman hominids having a less radioulnarly and dorsovolarly curved articular surface. All nonhuman hominids have more curved articular surface with Gorilla showing the most curved joint. Moreover, this study highlights the presence of a radially extended surface in Homo that may be related to the greater thumb abduction in human manipulation activities. Both fossils analyzed show a great ape-like MC 1 proximal articular surface which, associated with recent trabecular and archaeological evidence, may indicate that the ability to make/use stone tools preceded the morphological adaptations associated today with such behavior.     

Detecting strong positive selection in the genome

New statistical tests have been developed in the past decade that enable us to infer evidence of recent strong positive selection from genome-wide data on single-nucleotide polymorphism and to localize the targets of selection in the genome. Based on these tests, past demographic events that led to distortions of the site-frequency spectrum of variation can be distinguished from selection, in particular if linkage disequilibrium is taken into account. These methods have been successfully applied to species from which complete sequence information and polymorphism data are available, including Drosophila melanogaster, humans, and several plant species. To make full use of the available data, however, the tests that were primarily designed for panmictic populations need to be extended to spatially structured populations.     

DNA条形码在中国金合欢属中的应用

根据形态特征难以准确地辨别金合欢属植物,DNA条形码技术提供了一种准确地鉴定物种的方法。本文利用条形码技术对中国金合欢属物种的序列(psbA trnH、matK、rbcL和ITS)及其不同组合进行比较,通过计算种内和种间变异进行barcoding gap分析,运用Wilcoxon秩和检验比较不同序列的变异性,构建系统树。结果表明：4个片段均存在barcoding gap,ITS序列种间变异率较psbA trnH、rbcL和matK序列有明显优势,单片段ITS正确鉴定率最高,ITS+rbcL片段联合条码的正确鉴定率最高,因此我们认为ITS片段或条形码组合ITS+rbcL是金合欢属的快速鉴别最理想的条码。     

Structural analysis of a cDNA clone from Xenopus laevis containing a repetitive sequence element

Total polysomal RNA from Xenopus laevis stage 40 embryos was probed for the presence of repetitive sequences by Northern blot analysis with a genomic DNA fragment which had previously been shown to contain several repetitive sequence elements (Spohr et al., 1981). The analysis revealed that various presumptive mRNAs contain sequences complementary to the repetitive probe. Consequently, a cDNA library was constructed and screened with the same probe. Forty-eight positive recombinants containing eucaryotic inserts of 300–700 base pairs were isolated and one such clone was characterized in detail. Analysis of its nucleotide sequence revealed the presence of an open reading frame for 118 amino acids. Comparison of nucleotide sequences located 3′ to this presumptive protein coding region with the sequence of the genomic DNA fragment used as a probe clearly identifies and allows one to define the exact location of the repetitive element in the cloned cDNA. This analysis shows furthermore that one portion of the repeated sequence is highly conserved in the two members of this repetitive sequence family, whereas the other part is more divergent. In this area blocks of oligonucleotides are scattered between nonhomologous DNA stretches. The occurrence frequency of the presumptive mRNAs which carry repetitive elements homologous to the used repetitive probe is suggested to be close to that of rare mRNAs.