A visualization of 3D proteome universe: Mapping of a proteome ensemble into 3D space based on the protein-structure composition

To visualize a bird’s-eye view of an ensemble of proteomes for various species, we recently developed a novel method of mapping a proteome ensemble into Three-Dimensional (3D) vector space. In this study, the “proteome” is defined as the entire set of all proteins encoded in a genome sequence, and these proteins were dealt with at the level of the SCOP Fold. First, we represented the proteome of a species s by a 1053-dimensional vector x(s), where its length ∣x(s)∣ represents the overall amount of all the SCOP Folds in the proteome, and its unit vector x(s)/∣x(s)∣ represents the relative composition of the SCOP Folds in the proteome and the size of the dimension, 1053, is the number of all possible Folds in the proteome ensemble given. Second, we mapped the vector x(s) to the 3D vector y(s), based on the two simple principles: (1) ∣y(s)∣ = ∣x(s)∣, and (2) the angle between y(s) and y(t) maximally correlates with the angle between x(s) and x(t). We applied to the mapping of a proteome ensemble for 456 species, which were retrieved from the Genomes TO Protein structures and functions (GTOP) database. As a result, we succeeded in the mapping in that the properties of the 1053-dimensional vectors were quantitatively conserved in the 3D vectors. Particularly, the angles between vectors before and after the mapping highly correlated with each other (correlation coefficients were 0.95–0.96). This new mapping method will allow researchers to intuitively interpret the visual information presented in the maps in a highly effective manner.     

Phylogenetic relationships among eukaryotic kingdoms inferred from ribosomal RNA sequences

Summary Phylogenetic trees among eukaryotic kingdoms were inferred for large- and small-subunit rRNAs by using a maximum-likelihood method developed by Felsenstein. Although Felsenstein's method assumes equal evolutionary rates for transitions and transversions, this is apparently not the case for these data. Therefore, only transversiontype substitutions were taken into account. The molecules used were large-subunit rRNAs fromXenopus laevis (Animalia), rice (Plantae),Saccharomyces cerevisiae (Fungi),Dictyostelium discoideum (Protista), andPhysarum polycephalum (Protista); and small-subunit rRNAs from maize (Plantae),S. cerevisiae, X. laevis, rat (Animalia), andD. discoideum. Only conservative regions of the nucleotide sequences were considered for this study. In the maximum-likelihood trees for both large- and small-subunit rRNAs, Animalia and Fungi were the most closely related eukaryotic kingdoms, and Plantae is the next most closely related kingdom, although other branching orders among Plantae, Animalia, and Fungi were not excluded by this work. These three eukaryotic kingdoms apparently shared a common ancestor after the divergence of the two species of Protista,D. discoideum andP. polycephalum. These two species of Protista do not form a clade, andP. polycephalum diverged first andD. discoideum second from the line leading to the common ancestor of Plantae, Animalia, and Fungi. The sequence data indicate that a drastic change occurred in the nucleotide sequences of rRNAs during the evolutionary separation between prokaryote and eukaryote.     

Characterization of fibrillatory rhythms by ensemble vector directional analysis

Recent studies have demonstrated that fibrillatory rhythms are not random phenomena but have definable patterns. However, standard mapping techniques may have limitations in their ability to identify the organization of fibrillation. The purpose of this study was to develop and apply a method, "ensemble vector mapping," for characterizing the spatiotemporal organization of fibrillation. Ventricular fibrillation was induced by burst pacing in normal mongrel dogs. In a separate protocol, atrial fibrillation was induced by epicardial aconitine application. Epicardial electrograms were recorded from a 112-electrode plaque array using a computerized mapping system. Vectors were created by summing orthogonal bipolar electrograms. The magnitude of the vectors was transformed using a logarithmic function, integrated over time, and normalized for local electrogram amplitude to produce an "ensemble vector" index whose magnitude is high when beat-to-beat activation direction is consistent and low when activation direction is variable. The mean index was 137 +/- 36 mV/s during ventricular pacing at a cycle length of 300 ms but only 39 +/- 23 mV/s during ventricular fibrillation (P < 0.001). The ensemble vector index was also lower during atrial fibrillation (60 +/- 54 mV/s) than during atrial pacing (115 +/- 27 mV/s, P < 0.01 vs. atrial fibrillation) but not as low as during ventricular fibrillation (P < 0.05, atrial vs. ventricular fibrillation). The index was also capable of distinguishing atrial tachycardia from atrial fibrillation. Ensemble vector mapping produces an objective assessment of the consistency of myocardial activation during fibrillation. The consistency of activation direction differs in different models of fibrillation and is higher during atrial than ventricular fibrillation. This technique has the potential to rapidly characterize repetitive activation patterns in fibrillatory rhythms and may help distinguish among different characteristics of fibrillatory rhythms.     

Dynamic light scattering from dilute suspensions of thin discs and thin rods as limiting forms of cylinder, ellipsoid and ellipsoidal shell of revolution

A previous formulation of the field correlation function G1(tau) of light quasielastically scattered from suspensions of rigid rods undergoing anisotropic translational as well as rotational diffusion (T. Maeda and S. Fujime, Macromolecules 17 (1984) 1157) was extended to the cases of suspensions of cylinders (length L and radius R), ellipsoids and ellipsoidal shells of revolution (x2/b2 + y2/b2 + z2/a2 = 1). The present formulation includes that for suspensions of rigid rods in the limit of KR 1 or in the limit of b/a 1 and Kb 1 (an extremely prolate ellipsoid), and also that for suspensions of discs in the limit of KL 1 or in the limit of b/a 1 and Ka 1 (an extremely oblate ellipsoid), where K is the length of the scattering vector. Explicit forms of G1(tau), of the first cumulant Gamma of G1(tau) and of the dynamic form factors will be given, and numerical methods suitable for computation of dynamic form factors will be discussed. The present results can be applied to the analysis of experimental data for dilute suspensions of thin rods and thin discs. When the situation is favorable, our method can provide transport coefficients D1, D3, and Theta from dynamic light-scattering data only, where D(1) and D(3) are, respectively, the translational diffusion coefficients parallel with the x (y) and z axes, and Theta the rotational diffusion coefficient around the x (y) axis.     

A genetic linkage map for azuki bean [Vigna angularis (Willd.) Ohwi & Ohashi]

To make progress in genome analysis of azuki bean (Vigna angularis) a genetic linkage map was constructed from a backcross population of (V. nepalensis x V. angularis) x V.angularis consisting of 187 individuals. A total of 486 markers—205 simple sequence repeats (SSRs), 187 amplified fragment length polymorphisms (AFLPs) and 94 restriction fragment length polymorphisms (RFLPs) —were mapped onto 11 linkage groups corresponding to the haploid chromosome number of azuki bean. This map spans a total length of 832.1 cM with an average marker distance of 1.85 cM and is the most saturated map for a Vigna species to date. In addition, RFLP markers from other legumes facilitated finding several orthologous linkage groups based on previously published RFLP linkage maps. Most SSR primers that have been developed from SSR-enriched libraries detected a single locus. The SSR loci identified are distributed throughout the azuki bean genome. This moderately dense linkage map equipped with many SSR markers will be useful for mapping a range of useful traits such as those related to domestication and stress resistance. The mapping population will be used to develop advanced backcross lines for high resolution QTL mapping of these traits. O.K. Han, A. Kaga, T. Isemura have contributed equally to this paper.     

New Molecular Markers for the Distal End of the T-Complex and Their Relationships to Mutations Affecting Mouse Development

Many mutations affecting mouse development have been mapped to the t-complex of mouse chromosome 17. We have obtained 17 cosmid clones as molecular markers for this region by screening a hamster-mouse chromosome 17 and 18 cell hybrid cosmid library with mouse-specific repetitive elements and mapping positive clones via t-haplotype vs. C3H restriction fragment length polymorphism (RFLP) analysis. Twelve of the clones mapping distal to Leh66B in t-haplotypes are described here. Using standard RFLP analysis or simple sequence length polymorphism between t-haplotypes, exceptional partial t-haplotypes and nested sets of inter-t-haplotype recombinants, five cosmids have been mapped in or around In(17)3 and seven in the most distal inversion In17(4). More precise mapping of four of the cosmids from In(17)4 shows that they will be useful in the molecular identification of some of the recessive lethals mapped to the t-complex: two cosmids map between H-2K and Crya-1, setting a distal limit in t-haplotypes for the position of the tw5 lethal, one is inseparable from the tw12 lethal, and one maps distal to tf near the t0(t6) lethal and cld.     

Utilization of deletion bins to anchor and order sequences along the wheat 7B chromosome

Key message

A total of 3,671 sequence contigs and scaffolds were mapped to deletion bins on wheat chromosome 7B providing a foundation for developing high-resolution integrated physical map for this chromosome.

Abstract

Bread wheat (Triticum aestivum L.) has a large, complex and highly repetitive genome which is challenging to assemble into high quality pseudo-chromosomes. As part of the international effort to sequence the hexaploid bread wheat genome by the international wheat genome sequencing consortium (IWGSC) we are focused on assembling a reference sequence for chromosome 7B. The successful completion of the reference chromosome sequence is highly dependent on the integration of genetic and physical maps. To aid the integration of these two types of maps, we have constructed a high-density deletion bin map of chromosome 7B. Using the 270 K Nimblegen comparative genomic hybridization (CGH) array on a set of cv. Chinese spring deletion lines, a total of 3,671 sequence contigs and scaffolds (~7.8 % of chromosome 7B physical length) were mapped into nine deletion bins. Our method of genotyping deletions on chromosome 7B relied on a model-based clustering algorithm (Mclust) to accurately predict the presence or absence of a given genomic sequence in a deletion line. The bin mapping results were validated using three different approaches, viz. (a) PCR-based amplification of randomly selected bin mapped sequences (b) comparison with previously mapped ESTs and (c) comparison with a 7B genetic map developed in the present study. Validation of the bin mapping results suggested a high accuracy of the assignment of 7B sequence contigs and scaffolds to the 7B deletion bins.     

Rat mutations cvd and hob with cerebellar malformations map to chromosome 2.

In this paper, we executed genome mapping and comparative mapping analyses for cvd and hob, autosomal recessive mutations with cerebellar vermis defect and cerebellar dysplasia in the rat. For the linkage analysis, we produced three sets of backcross progeny, (ACI x CVD)F(1) and (F344 x CVD)F(1) females crossed to a cvd homozygous male rat, and (HOB x WKY)F(1) males crossed to hob homozygous female rats. Analysis of the segregation patterns of simple sequence length polymorphism (SSLP) markers scanning the whole rat genome allowed the mapping of these autosomal recessive mutations to rat Chromosome (Chr) 2. The most likely gene order is D2Mgh12 - D2Rat86 - D2Mit15 - D2Rat185 - cvd - D2Rat66 - D2Mgh13, and D2Mit18 - Fga -D2Mit14 - D2Rat16 - hob - D2Mgh13. Crossing test between a proven cvd heterozygous and a hob heterozygous rats demonstrated their allelism. Furthermore, comparative mapping indicated the cvd locus corresponds to mouse chromosome 3 and a strong candidate gene Unc5h3, a causative gene for the rostral cerebellar malformation mouse, was implicated.     

Phylogenetic position of Dictyostelium inferred from multiple protein data sets

The phylogenetic position of Dictyostelium inferred from 18S rRNA data contradicts that from protein data. Protein trees always show the close affinity of Dictyostelium with animals, fungi, and plants, whereas in 18S rRNA trees the branching of Dictyostelium is placed at a position before the massive radiation of protist groups including the divergence of the three kingdoms. To settle this controversial issue and to determine the correct position of Dictyostelium, we inferred the phylogenetic relationship among Dictyostelium and the three kingdoms Animalia, Fungi, and Plantae by a maximum-likelihood method using 19 different protein data sets. It was shown at the significance level of 1 SE that the branching of Dictyostelium antedates the divergence of Animalia and Fungi, and Plantae is an outgroup of the Animalia-Fungi-Dictyostelium clade.Correspondence to: T. Miyata     

A comparison of genetic maps constructed from haploid and BC1 mapping populations from the same crossing between Gossypium hirsutum L. and Gossypium barbadense L.

Simple sequence repeat (SSR) genetic maps have been separately constructed based on doubled haploid (DH) and (or) haploid and BC1 populations from the same cross between Gossypium hirsutum L. 'TM-1' and Gossypium barbadense L. 'Hai7124'. The BC1 population was produced by pollinating individual plants of the 'TM-1' x 'Hai7124' F1 with 'TM-1', whereas the DH and (or) haploid population developed from the offspring of Vsg x ('TM-1' x 'Hai7124'). Vsg is a virescently marked semigamy line of Gossypium barbadense L. Pima. The BC1 map included 34 linkage groups with an average distance between markers of 9.80 cM (Kosambi, K) and covered 4331.2 cM (K) or approximately 78.7% of the tetraploid cotton genome constructed using 440 SSR and 2 morphological marker genes. Among them, 26 were assigned to 20 chromosomes, 7 to A or D subgenomes, and 1 was unassigned. The haploid map comprised 444 SSR markers mapped to 40 linkage groups with an average distance of 7.35 cM (K) between markers, covering 3262.9 cM (K) or approximately 60.0% of the tetraploid genome. Twenty-nine linkage groups were assigned to all 19 identified chromosomes, 10 to A or D subgenomes, and 1 was unassigned. Fairly good collinearity of marker order was observed along most of the chromosomes or linkage groups. Significant differences in recombination between maps was observed at the chromosomal and genomic level and possible reasons were discussed. Map comparison and combined data provided an essential basis for further mapping of interested genes and QTLs and for studies of diversity, population structure, and phylogeny in Gossypium species.     

Mapping an initiation region of DNA replication at a single-copy chromosomal locus in Drosophila melanogaster cells by two-dimensional gel methods and PCR-mediated nascent-strand analysis: multiple replication origins in a broad zone.

We have mapped an initiation region of DNA replication at a single-copy chromosomal locus in exponentially proliferating Drosophila tissue culture cells, using two-dimensional (2D) gel replicon mapping methods and PCR-mediated analysis of nascent strands. The initiation region was first localized downstream of the DNA polymerase alpha gene by determining direction of replication forks with the neutral/alkaline 2D gel method. Distribution of replication origins in the initiation region was further analyzed by using two types of 2D gel methods (neutral/neutral and neutral/alkaline) and PCR-mediated nascent-strand analysis. Results obtained by three independent methods were essentially consistent with each other and indicated that multiple replication origins are distributed in a broad zone of approximately 10 kb. The nucleotide sequence of an approximately 20-kb region that encompasses the initiation region was determined and searched for sequence elements potentially related to function of replication origins.     

Dynamic light scattering investigations of human erythrocyte spectrin.

Dynamic light scattering measurements were performed on spectrin from human erythrocytes in 25 mM Tris buffer at pH 7.6 with 100 mM NaCl and 5 mM EDTA. Measurements were made on spectrin solutions prepared as dimers and tetramers over the temperature range from 23 to 41 degrees C, as a function of the square of the scattering vector (K2) over the range of 0.7 x 10(10) cm-2 less than or equal to K1 less than or equal to 20 x 10(10) cm-2. Analysis of the autocorrelation functions collected for these solutions revealed the presence of two predominant motional components over the entire range of K2. Plots of the diffusion coefficients (D20) of these components, with viscosity and temperature corrected to water at 20 degrees C, as a function of K2 indicated three rather distinct regions, flat regions at low and high K2 joined by a sloping intermediate region. At small K2 (less than or equal to 4 x 10(10) cm-2) the D20 values were (7.3 +/- 2.0) x 10(-8) cm2/s for the slow component and (20.3 +/- 2.0) x 10(-8) cm2/s for the fast component. At large K2 (greater than or equal to 10 x 10(10) cm-2) the values increased to (13.0 +/- 2.0) x 10(-8) cm2/s for the slow component and (39.4 +/- 2.0) x 10(-8) cm2/s for the fast component. In the intermediate K2 region, D20 is a linear function of K2 and appears as a transition between the low and high K2 regions.(ABSTRACT TRUNCATED AT 250 WORDS)     

Predictability of sequence homologies among lysine-rich histones by immunological distance

The relationship between immunological distance (I.D.) measured by microcomplement fixation and amino acid sequence difference for lysine-rich histones was tested using antisera to lysine-rich histones of known sequence, chicken H1 and H5, goose H5, and trout H1 as well as to trout H5. The best relationship between I.D. (y) and percent sequence difference (x) for lysine-rich histones, y = 2x, applies as well to other histones of known sequence but it differs from y = 5x, reported for other proteins and often used for histones. Although deviations indicate that I.D. is a poor predictor of primary sequence differences among histones, it suggests that trout H5 is more closely related to H1 than to chicken H5.     

Genetic mapping of 66 new microsatellite (SSR) loci in bread wheat

In hexaploid bread wheat ( Triticum aestivum L. em. Thell), ten members of the IWMMN ( International Wheat Microsatellites Mapping Network) collaborated in extending the microsatellite (SSR = simple sequence repeat) genetic map. Among a much larger number of microsatellite primer pairs developed as a part of the WMC ( Wheat Microsatellite Consortium), 58 out of 176 primer pairs tested were found to be polymorphic between the parents of the ITMI ( International Triticeae Mapping Initiative) mapping population W7984 x Opata 85 (ITMI pop). This population was used earlier for the construction of RFLP ( Restriction Fragment Length Polymorphism) maps in bread wheat (ITMI map). Using the ITMI pop and a framework map (having 266 anchor markers) prepared for this purpose, a total of 66 microsatellite loci were mapped, which were distributed on 20 of the 21 chromosomes (no marker on chromosome 6D). These 66 mapped microsatellite (SSR) loci add to the existing 384 microsatellite loci earlier mapped in bread wheat.     

模块神经网络及其性能

提出一个不同于实现ｙ＝ｆ（ｘ）的ＢＰ网络的神经网络模型,给出了网络的结构及并行动力学方程,证明了其动力学的稳定性。通过学习算法的建立,证明网络能精确实现输入矢量对（ｘ,ｙ）映入成相联系的输出矢量ｚ,最重要的是网络能同时存诸依时变化的时序模式与静态模式。此外并给出动力学学习算法,证明此学习算法的收敛性,计算机仿真证实理论结果,最后讨论了某些可能的应用。     

Mapping Accuracy of Short Reads from Massively Parallel Sequencing and the Implications for Quantitative Expression Profiling

Background

Massively parallel sequencing offers an enormous potential for expression profiling, in particular for interspecific comparisons. Currently, different platforms for massively parallel sequencing are available, which differ in read length and sequencing costs. The 454-technology offers the highest read length. The other sequencing technologies are more cost effective, on the expense of shorter reads. Reliable expression profiling by massively parallel sequencing depends crucially on the accuracy to which the reads could be mapped to the corresponding genes.

Methodology/Principal Findings

We performed an in silico analysis to evaluate whether incorrect mapping of the sequence reads results in a biased expression pattern. A comparison of six available mapping software tools indicated a considerable heterogeneity in mapping speed and accuracy. Independently of the software used to map the reads, we found that for compact genomes both short (35 bp, 50 bp) and long sequence reads (100 bp) result in an almost unbiased expression pattern. In contrast, for species with a larger genome containing more gene families and repetitive DNA, shorter reads (35–50 bp) produced a considerable bias in gene expression. In humans, about 10% of the genes had fewer than 50% of the sequence reads correctly mapped. Sequence polymorphism up to 9% had almost no effect on the mapping accuracy of 100 bp reads. For 35 bp reads up to 3% sequence divergence did not affect the mapping accuracy strongly. The effect of indels on the mapping efficiency strongly depends on the mapping software.

Conclusions/Significance

In complex genomes, expression profiling by massively parallel sequencing could introduce a considerable bias due to incorrectly mapped sequence reads if the read length is short. Nevertheless, this bias could be accounted for if the genomic sequence is known. Furthermore, sequence polymorphisms and indels also affect the mapping accuracy and may cause a biased gene expression measurement. The choice of the mapping software is highly critical and the reliability depends on the presence/absence of indels and the divergence between reads and the reference genome. Overall, we found SSAHA2 and CLC to produce the most reliable mapping results.     

Cloning and genomic localization of the murine LPS-induced CXC chemokine (LIX) gene,Scyb5

LPS-induced CXC chemokine (LIX) is a murine chemokine similar to two human chemokines, ENA-78 (CXCL5) and GCP-2 (CXCL6). To clarify the relationship of LIX to human ENA-78 and GCP-2, we cloned and mapped the LIX gene. The organization of the LIX gene ( Scyb5) is similar to those of the human ENA-78 ( SCYB5) and GCP-2 ( SCYB6) genes. The intron-exon boundaries of the three genes are exactly conserved, and the introns have similar sizes. The first 100 bp of the 5' flanking regions are highly similar, with conserved NF-kappaB and GATA sites in identical positions in all three genes. Further 5', the Lix flanking region sequence diverges from those of ENA-78 and GCP-2, which remain highly similar for 350 bp preceding the start sites. Using a (C57BL/6 J x Mus spretus) F1 x C57BL/6J backcross panel, Lix was mapped to a locus near D5Ucla5 at 49.0 cM on Chromosome (Chr) 5. Mapping with the T31 radiation hybrid panel placed Lix between D5Mit360 and D5Mit6. Physical maps of the CXC chemokine clusters on murine Chr 5 and human Chr 21 were constructed using the Celera mouse genome database and the public human genome database. The sequence and mapping data suggest that the human ENA78-PBP-PF4 and GCP2- psi PBP-PF4V1 loci arose from an evolutionarily recent duplication of an ancestral locus related to the murine Lix-Pbp-Pf4 locus.     

Are red algae plants?

For 200 years prior to the 1938 publication of H. F. Copeland, all authorities (with one exception) classified red algae (Rhodophyta) within Kingdom Plantae or its equivalent. Copeland's reclassification of red algae within Kingdom Protista or Protoctista drew from an alternative tradition, dating to Cohn in 1867, in which red algae were viewed as the earliest or simplest eukaryotes. Analyses of ribosomal RNA (rRNA) sequence data initially favoured Copeland's reclassification. Many more rRNA gene (rDNA) sequences are now available from the eukaryote lineages most closely related to red algae, and based on these data, the hypothesis that red algae and green plants are sister groups cannot be rejected. An increasing body of sequence, intron-location and functional data from nuclear- and mitochondrially encoded proteins likewise supports a sister-group relationship between red algae and green plants. Submerging Kingdoms Plantae, Animalia and Fungi into Eukarya would provide a more natural framework for the eventual resolution of whether red algae are plants or prorists.