ParAlign: a parallel sequence alignment algorithm for rapid and sensitive database searches

There is a need for faster and more sensitive algorithms for sequence similarity searching in view of the rapidly increasing amounts of genomic sequence data available. Parallel processing capabilities in the form of the single instruction, multiple data (SIMD) technology are now available in common microprocessors and enable a single microprocessor to perform many operations in parallel. The ParAlign algorithm has been specifically designed to take advantage of this technology. The new algorithm initially exploits parallelism to perform a very rapid computation of the exact optimal ungapped alignment score for all diagonals in the alignment matrix. Then, a novel heuristic is employed to compute an approximate score of a gapped alignment by combining the scores of several diagonals. This approximate score is used to select the most interesting database sequences for a subsequent Smith–Waterman alignment, which is also parallelised. The resulting method represents a substantial improvement compared to existing heuristics. The sensitivity and specificity of ParAlign was found to be as good as Smith–Waterman implementations when the same method for computing the statistical significance of the matches was used. In terms of speed, only the significantly less sensitive NCBI BLAST 2 program was found to outperform the new approach. Online searches are available at http://dna.uio.no/search/

Extending the Microtubule/Microfibril Paradigm1 : Cellulose Synthesis Is Required for Normal Cortical Microtubule Alignment in Elongating Cells

The cortical microtubule array provides spatial information to the cellulose-synthesizing machinery within the plasma membrane of elongating cells. Until now data indicated that information is transferred from organized cortical microtubules to the cellulose-synthesizing complex, which results in the deposition of ordered cellulosic walls. How cortical microtubules become aligned is unclear. The literature indicates that biophysical forces, transmitted by the organized cellulose component of the cell wall, provide a spatial cue to orient cortical microtubules. This hypothesis was tested on tobacco (Nicotiana tabacum L.) protoplasts and suspension-cultured cells treated with the cellulose synthesis inhibitor isoxaben. Isoxaben (0.25–2.5 μm) inhibited the synthesis of cellulose microfibrils (detected by staining with 1 μg mL−1 fluorescent dye and polarized birefringence), the cells failed to elongate, and the cortical microtubules failed to become organized. The affects of isoxaben were reversible, and after its removal microtubules reorganized and cells elongated. Isoxaben did not depolymerize microtubules in vivo or inhibit the polymerization of tubulin in vitro. These data are consistent with the hypothesis that cellulose microfibrils, and hence cell elongation, are involved in providing spatial cues for cortical microtubule organization. These results compel us to extend the microtubule/microfibril paradigm to include the bidirectional flow of information.

Antioxidant Defenses in the Peripheral Cell Layers of Legume Root Nodules1

Ascorbate peroxidase (AP) is a key enzyme that scavenges potentially harmful H2O2 and thus prevents oxidative damage in plants, especially in N2-fixing legume root nodules. The present study demonstrates that the nodule endodermis of alfalfa (Medicago sativa) root nodules contains elevated levels of AP protein, as well as the corresponding mRNA transcript and substrate (ascorbate). Enhanced AP protein levels were also found in cells immediately peripheral to the infected region of soybean (Glycine max), pea (Pisum sativum), clover (Trifolium pratense), and common bean (Phaseolus vulgaris) nodules. Regeneration of ascorbate was achieved by (homo)glutathione and associated enzymes of the ascorbate-glutathione pathway, which were present at high levels. The presence of high levels of antioxidants suggests that respiratory consumption of O2 in the endodermis or nodule parenchyma may be an essential component of the O2-diffusion barrier that regulates the entry of O2 into the central region of nodules and ensures optimal functioning of nitrogenase.

Multiple DNA and protein sequence alignment based on segment-to-segment comparison.

In this paper, a new way to think about, and to construct, pairwise as well as multiple alignments of DNA and protein sequences is proposed. Rather than forcing alignments to either align single residues or to introduce gaps by defining an alignment as a path running right from the source up to the sink in the associated dot-matrix diagram, we propose to consider alignments as consistent equivalence relations defined on the set of all positions occurring in all sequences under consideration. We also propose constructing alignments from whole segments exhibiting highly significant overall similarity rather than by aligning individual residues. Consequently, we present an alignment algorithm that (i) is based on segment-to-segment comparison instead of the commonly used residue-to-residue comparison and which (ii) avoids the well-known difficulties concerning the choice of appropriate gap penalties: gaps are not treated explicity, but remain as those parts of the sequences that do not belong to any of the aligned segments. Finally, we discuss the application of our algorithm to two test examples and compare it with commonly used alignment methods. As a first example, we aligned a set of 11 DNA sequences coding for functional helix-loop-helix proteins. Though the sequences show only low overall similarity, our program correctly aligned all of the 11 functional sites, which was a unique result among the methods tested. As a by-product, the reading frames of the sequences were identified. Next, we aligned a set of ribonuclease H proteins and compared our results with alignments produced by other programs as reported by McClure et al. [McClure, M. A., Vasi, T. K. & Fitch, W. M. (1994) Mol. Biol. Evol. 11, 571-592]. Our program was one of the best scoring programs. However, in contrast to other methods, our protein alignments are independent of user-defined parameters.

Gene recognition via spliced sequence alignment.

Gene recognition is one of the most important problems in computational molecular biology. Previous attempts to solve this problem were based on statistics, and applications of combinatorial methods for gene recognition were almost unexplored. Recent advances in large-scale cDNA sequencing open a way toward a new approach to gene recognition that uses previously sequenced genes as a clue for recognition of newly sequenced genes. This paper describes a spliced alignment algorithm and software tool that explores all possible exon assemblies in polynomial time and finds the multiexon structure with the best fit to a related protein. Unlike other existing methods, the algorithm successfully recognizes genes even in the case of short exons or exons with unusual codon usage; we also report correct assemblies for genes with more than 10 exons. On a test sample of human genes with known mammalian relatives, the average correlation between the predicted and actual proteins was 99%. The algorithm correctly reconstructed 87% of genes and the rare discrepancies between the predicted and real exon-intron structures were caused either by short (less than 5 amino acids) initial/terminal exons or by alternative splicing. Moreover, the algorithm predicts human genes reasonably well when the homologous protein is nonvertebrate or even prokaryotic. The surprisingly good performance of the method was confirmed by extensive simulations: in particular, with target proteins at 160 accepted point mutations (PAM) (25% similarity), the correlation between the predicted and actual genes was still as high as 95%.

Retinal engineering: engrafted neural cell lines locate in appropriate layers.

A major question in central nervous system development, including the neuroretina, is whether migrating cells express cues to find their way and settle at specific locations. We have transplanted quail neuroretinal cell lines QNR/D, a putative amacrine or ganglion cell, and QNR/K2, a putative Müller cell into chicken embryo eyes. Implanted QNR/D cells migrate only to the retinal ganglion and amacrine cell layers and project neurites in the plane of retina; in contrast, QNR/K2 cells migrate through the ganglion and amacrine layers, locate in the inner nuclear layer, and project processes across the retina. These data show that QNR/D and QNR/K2 cell lines represent distinct neural cell types, suggesting that migrating neural cells express distinct address cues. Furthermore, our results raise the possibility that immortalized cell lines can be used for replacement of specific cell types and for the transport of genes to given locations in neuroretina.

Calicheamicin-DNA complexes: warhead alignment and saccharide recognition of the minor groove.

The solution structures of calicheamicin gamma 1I, its cycloaromatized analog (calicheamicin epsilon), and its aryl tetrasaccharide complexed to a common DNA hairpin duplex have been determined by NMR and distance-refined molecular dynamics computations. Sequence specificity is associated with carbohydrate-DNA recognition that places the aryl tetrasaccharide component of all three ligands in similar orientations in the minor groove at the d(T-C-C-T).d(A-G-G-A) segment. The complementary fit of the ligands and the DNA minor groove binding site creates numerous van der Waals contacts as well as hydrogen bonding interactions. Notable are the iodine and sulfur atoms of calicheamicin that hydrogen bond with the exposed amino proton of the 5'- and 3'-guanines, respectively, of the d(A-G-G-A) segment. The sequence-specific carbohydrate binding orients the enediyne aglycone of calicheamicin gamma 1I such that its C3 and C6 proradical centers are adjacent to the cleavage sites. While the enediyne aglycone of calicheamicin gamma 1I is tilted relative to the helix axis and spans the minor groove, the cycloaromatized aglycone is aligned approximately parallel to the helix axis in the respective complexes. Specific localized conformational perturbations in the DNA have been identified from imino proton complexation shifts and changes in specific sugar pucker patterns on complex formation. The helical parameters for the carbohydrate binding site are comparable with corresponding values in B-DNA fibers while a widening of the groove is observed at the adjacent aglycone binding site.

Alignment of a 1.2-Mb chromosomal region from three strains of Rhodobacter capsulatus reveals a significantly mosaic structure.

High-resolution physical maps of the genomes of three Rhodobacter capsulatus strains, derived from ordered cosmid libraries, were aligned. The 1.2-Mb segment of the SB1003 genome studied here is adjacent to a 1-Mb region analyzed previously [Fonstein, M., Nikolskaya, T. & Haselkorn, H. (1995) J. Bacteriol. 177, 2368-2372]. Probes derived from the ordered cosmid set of R. capsulatus SB1003 were used to link cosmids from the St. Louis and 2.3.1 strain libraries. Cosmids selected this way did not merge into a single contig but formed several unlinked groups. EcoRV restriction maps of the ordered cosmids were then constructed using lambda terminase and fused to derive fragments of the chromosomal map. In order to link these fragments, their ends were transcribed to produce secondary probes for hybridization to gridded cosmid libraries of the same strains. This linking reduced the number of subcontigs to three for the St. Louis strain and one for the 2.3.1 strain. Hybridization of the same probes back to the ordered cosmid set of SB1003 positioned the subcontigs on the high-resolution physical map of SB1003. The final alignment of the restriction maps shows numerous large and small translocations in this 1.2-Mb chromosomal region of the three Rhodobacter strains. In addition, the chromosomes of the three strains, whose fine-structure maps can now be compared over 2.2 Mb, are seen to contain regions of 15-80 kb in which restriction sites are highly polymorphic, interspersed among regions in which the positions of restriction sites are highly conserved.

Rat growth hormone gene introns stimulate nucleosome alignment in vitro and in transgenic mice.

Average hepatic expression (mRNA per cell per gene) of a metallothionein-rat growth hormone (rGH) gene with its natural introns was about 15-fold higher than an intronless version when tested in transgenic mice. We examined the idea that intron removal leads to an alteration in chromatin structure that might be responsible for this effect. Using an in vitro chromatin assembly system, we observed that nucleosomes were aligned in a characteristic ordered array over the gene and promoter when all introns were present. Linker histones were necessary for this alignment to occur. In contrast, nucleosome alignment was perturbed in constructs lacking some or all of the introns. A similar disruption of nucleosome alignment was observed when comparing chromatin from livers of transgenic mice carrying rGH transgenes with or without introns. In vitro, sequences at the 3' end of the rGH gene position nucleosomes and facilitate nucleosome alignment upstream; however, nucleosome alignment does not occur on the approximately 3 kb of downstream flanking rat sequence. These observations suggest that signals present in genomic rGH DNA may serve to establish appropriate nucleosome alignment during development and, possibly, to restore nucleosome alignment to the transcribed region after disruption incurred by the passage of an RNA polymerase molecule, thereby facilitating subsequent rounds of transcription.

Oxide Semiconductors For Organic Opto-electronic Devices

Development of transparent oxide semiconductors (TOS) from Earth-abundant materials is of great interest for cost-effective thin film device applications, such as solar cells, light emitting diodes (LEDs), touch-sensitive displays, electronic paper, and transparent thin film transistors. The need of inexpensive or high performance electrode might be even greater for organic photovoltaic (OPV), with the goal to harvest renewable energy with inexpensive, lightweight, and cost competitive materials. The natural abundance of zinc and the wide bandgap ($sim$3.3 eV) of its oxide make it an ideal candidate. In this dissertation, I have introduced various concepts on the modulations of various surface, interface and bulk opto-electronic properties of ZnO based semiconductor for charge transport, charge selectivity and optimal device performance. I have categorized transparent semiconductors into two sub groups depending upon their role in a device. Electrodes, usually 200 to 500 nm thick, optimized for good transparency and transporting the charges to the external circuit. Here, the electrical conductivity in parallel direction to thin film, i.e bulk conductivity is important. And contacts, usually 5 to 50 nm thick, are optimized in case of solar cells for providing charge selectivity and asymmetry to manipulate the built in field inside the device for charge separation and collection. Whereas in Organic LEDs (OLEDs), contacts provide optimum energy level alignment at organic oxide interface for improved charge injections. For an optimal solar cell performance, transparent electrodes are designed with maximum transparency in the region of interest to maximize the light to pass through to the absorber layer for photo-generation, plus they are designed for minimum sheet resistance for efficient charge collection and transport. As such there is need for material with high conductivity and transparency. Doping ZnO with some common elements such as B, Al, Ga, In, Ge, Si, and F result in n-type doping with increase in carriers resulting in high conductivity electrode, with better or comparable opto-electronic properties compared to current industry-standard indium tin oxide (ITO). Furthermore, improvement in mobility due to improvement on crystallographic structure also provide alternative path for high conductivity ZnO TCOs. Implementing these two aspects, various studies were done on gallium doped zinc oxide (GZO) transparent electrode, a very promising indium free electrode. The dynamics of the superimposed RF and DC power sputtering was utilized to improve the microstructure during the thin films growth, resulting in GZO electrode with conductivity greater than 4000 S/cm and transparency greater than 90 %. Similarly, various studies on research and development of Indium Zinc Tin Oxide and Indium Zinc Oxide thin films which can be applied to flexible substrates for next generation solar cells application is presented. In these new TCO systems, understanding the role of crystallographic structure ranging from poly-crystalline to amorphous phase and the influence on the charge transport and optical transparency as well as important surface passivation and surface charge transport properties. Implementation of these electrode based on ZnO on opto-electronics devices such as OLED and OPV is complicated due to chemical interaction over time with the organic layer or with ambient. The problem of inefficient charge collection/injection due to poor understanding of interface and/or bulk property of oxide electrode exists at several oxide-organic interfaces. The surface conductivity, the work function, the formation of dipoles and the band-bending at the interfacial sites can positively or negatively impact the device performance. Detailed characterization of the surface composition both before and after various chemicals treatment of various oxide electrode can therefore provide insight into optimization of device performance. Some of the work related to controlling the interfacial chemistry associated with charge transport of transparent electrodes are discussed. Thus, the role of various pre-treatment on poly-crystalline GZO electrode and amorphous indium zinc oxide (IZO) electrode is compared and contrasted. From the study, we have found that removal of defects and self passivating defects caused by accumulation of hydroxides in the surface of both poly-crystalline GZO and amorphous IZO, are critical for improving the surface conductivity and charge transport. Further insight on how these insulating and self-passivating defects cause charge accumulation and recombination in an device is discussed. With recent rapid development of bulk-heterojunction organic photovoltaics active materials, devices employing ZnO and ZnO based electrode provide air stable and cost-competitive alternatives to traditional inorganic photovoltaics. The organic light emitting diodes (OLEDs) have already been commercialized, thus to follow in the footsteps of this technology, OPV devices need further improvement in power conversion efficiency and stable materials resulting in long device lifetimes. Use of low work function metals such as Ca/Al in standard geometry do provide good electrode for electron collection, but serious problems using low work-function metal electrodes originates from the formation of non-conductive metal oxide due to oxidation resulting in rapid device failure. Hence, using low work-function, air stable, conductive metal oxides such as ZnO as electrons collecting electrode and high work-function, air stable metals such as silver for harvesting holes, has been on the rise. Devices with degenerately doped ZnO functioning as transparent conductive electrode, or as charge selective layer in a polymer/fullerene based heterojunction, present useful device structures for investigating the functional mechanisms within OPV devices and a possible pathway towards improved air-stable high efficiency devices. Furthermore, analysis of the physical properties of the ZnO layers with varying thickness, crystallographic structure, surface chemistry and grain size deposited via various techniques such as atomic layer deposition, sputtering and solution-processed ZnO with their respective OPV device performance is discussed. We find similarity and differences in electrode property for good charge injection in OLEDs and good charge collection in OPV devices very insightful in understanding physics behind device failures and successes. In general, self-passivating surface of amorphous TCOs IZO, ZTO and IZTO forms insulating layer that hinders the charge collection. Similarly, we find modulation of the carrier concentration and the mobility in electron transport layer, namely zinc oxide thin films, very important for optimizing device performance.

Search for signatures of magnetically-induced alignment in the arrival directions measured by the Pierre Auger Observatory

We present the results of an analysis of data recorded at the Pierre Auger Observatory in which we search for groups of directionally-aligned events (or 'multiplets') which exhibit a correlation between arrival direction and the inverse of the energy. These signatures are expected from sets of events coming from the same source after having been deflected by intervening coherent magnetic fields. The observation of several events from the same source would open the possibility to accurately reconstruct the position of the source and also measure the integral of the component of the magnetic field orthogonal to the trajectory of the cosmic rays. We describe the largest multiplets found and compute the probability that they appeared by chance from an isotropic distribution. We find no statistically significant evidence for the presence of multiplets arising from magnetic deflections in the present data. (C) 2011 Elsevier B.V. All rights reserved.

Examining the impact of pre-k and kindergarten classroom alignment on school readiness outcomes

This study explored children’s experiences of instructional alignment from prekindergarten to kindergarten and analyzed the impact of those alignment experiences on children’s school readiness outcomes. The study answered the following overarching research question: Does the alignment of children’s learning experiences between prekindergarten and kindergarten impact school readiness outcomes? Three sub-questions drove the research design: (1) How do children’s prekindergarten and kindergarten learning experiences align; (2) To what extent does the alignment of early learning experiences predict children’s school readiness outcomes; and (3) Does the quality of prekindergarten classroom teacher interactions moderate the impact of any PK-K alignment effects? Using cluster analysis and hierarchical linear modeling (HLM) to analyze data from over 1,300 children in the 2009 Head Start Family and Child Experiences Survey (FACES), the study found that children have distinct and definable experiences of PK-K alignment. Results also indicated a disparity in children’s PK-K alignment experiences, with Hispanic/Latino children more likely to attend Head Start programs with poor systems transition practices followed by kindergartens with poor classroom structures. The study found that growth in the use of instructional activity centers from prekindergarten to kindergarten is predictive of better literacy and math outcomes. Findings further suggested that boys, minority students, and children from lower income households are predicted to score lower than girls, white classmates, and higher-income peers across school readiness measures. Findings support the need for equitable transition and alignment practices for children from all racial and ethnic groups. They also argue for an increase in child-directed activity centers in kindergarten. With one exception, the current findings did not support the hypothesis that prekindergarten teacher quality is a moderator of alignment effects on children’s school readiness outcomes. The study presents suggestions for further research.

Spin alignment of extra electrons in K-phenanthrene clusters taken from the crystalline tripotassium-intercalated phenanthrene structure

The appearance of ferromagnetic correlations among π electrons of phenanthrene (C14H10) molecules in the herringbone structure is proven for K doped clusters both by ab initio quantum-chemistry calculations and by the direct solution of the many-body Pariser-Parr-Pople Hamiltonian. Magnetic ground states are predicted for one or three additional electrons per phenanthrene molecule. These results are a consequence of the small overlap between the lowest unoccupied molecular orbitals (and lowest unoccupied molecular orbitals + 1) of neutral neighboring phenanthrene molecules, which makes the gain in energy by delocalization similar to the corresponding increase due to the Coulomb interaction.

Structure, gas-barrier properties and overall migration of poly(lactic acid) films coated with hydrogenated amorphous carbon layers

Hydrogenated amorphous carbon (a-C:H) films were grown on a poly(lactic acid) (PLA) substrate by means of a radiofrequency plasma-enhanced chemical vapour deposition (rf-PECVD) technique with different deposition times (5, 20 and 40 min). The main goal of this treatment was to increase the barrier properties of PLA, maintaining its original transparency and colour as well as controlling interactions with food simulants for packaging applications. Morphological, chemical, and mechanical properties of PLA/a-C:H systems were evaluated while permeability and overall migration tests were performed in order to determine the effect of the plasma treatment on the gas-barrier properties of PLA films and their application in food packaging. Morphological results suggested a good adhesion of the deposited layers onto the polymer surface and the samples treated for 5 and 20 min only slightly darkened the PLA film. X-ray photoelectron spectroscopy revealed that the structural properties of the carbon layer deposited onto the PLA film depend on the exposure time. PLA/a-C:H system treated for 5 min showed the highest barrier properties, while none of the studied samples exceeded the migration limit established by the current legislation, suggesting the suitability of these materials in packaging applications.