Molecular Characterization of Viruses Causing the Cassava Brown Streak Disease Epidemic in Eastern Africa

Cassava brown streak disease (CBSD) was described for the first time in Tanganyika (now Tanzania) about seven decades ago. Tanganyika (now Tanzania) about seven decades ago. It was endemic in the lowland areas of East Africa and inland parts of Malawi and caused by Cassava brown streak virus (CBSV; genus Ipomovirus; Potyviridae). However, in 1990s CBSD was observed at high altitude areas in Uganda. The causes for spread to new locations were not known.The present work was thus initiated to generate information on genetic variability, clarify the taxonomy of the virus or viruses associated with CBSD in Eastern Africa as well as to understand the evolutionary forces acting on their genes. It also sought to develop a molecular based diagnostic tool for detection of CBSD-associated virus isolates. Comparison of the CP-encoding sequences of CBSD-associated virus isolates collected from Uganda and north-western Tanzania in 2007 and the partial sequences available in Genbank revealed occurrence of two genetically distinct groups of isolates. Two isolates were selected to represent the two groups. The complete genomes of isolates MLB3 (TZ:Mlb3:07) and Kor6 (TZ:Kor6:08) obtained from North-Western (Kagera) and North-Eastern (Tanga) Tanzania, respectively, were sequenced. The genomes were 9069 and 8995 nucleotides (nt), respectively. They translated into polyproteins that were predicted to yield ten mature proteins after cleavage. Nine proteins were typical in the family Potyviridae, namely P1, P3, 6K1, CI, 6K2, VPg, NIa-Pro, NIb and CP, but the viruses did not contain HC-Pro. Interestingly, genomes of both isolates contained a Maf/HAM1-like sequence (HAM1h; 678 nucleotides, 25 kDa) recombined between the NIb and CP domains in the 3’-proximal part of the genomes. HAM1h was also identified in Euphorbia ringspot virus (EuRSV) whose sequence was in GenBank. The HAM1 gene is widely spread in both prokaryotes and eukaryotes. In yeast (Saccharomyces cerevisiae) it is known to be a nucleoside triphosphate (NTP) pyrophosphatase. Novel information was obtained on the structural variation at the N-termini of polyproteins of viruses in the genus Ipomovirus. Cucumber vein yellowing virus (CVYV) and Squash vein yellowing virus (SqVYV) contain a duplicated P1 (P1a and P1b) but lack the HC-Pro. On the other hand, Sweet potato mild mottle virus (SPMMV), has a single but large P1 and has HC-Pro. Both virus isolates (TZ:Mlb3:07 & TZ:Kor6:08) characterized in this study contained a single P1 and lacked the HC-Pro which indicates unique evolution in the family Potyviridae. Comparison of 12 complete genomes of CBSD-associated viruses which included two genomes characterized in this study, revealed genetic identity of 69.0–70.3% (nt) and amino acid (aa) identities of 73.6–74.4% at polyprotein level. Comparison was also made among 68 complete CP sequences, which indicated 69.0-70.3 and 73.6-74.4 % identity at nt and aa levels, respectively. The genetic variation was large enough for dermacation of CBSD-associated virus isolates into two distinct species. The name CBSV was retained for isolates that were related to CBSV isolates available in database whereas the new virus described for the first time in this study was named Ugandan cassava brown streak virus (UCBSV) by the International Committee on Virus Taxonomy (ICTV). The isolates TZ:Mlb3:07 and TZ:Kor6:08 belong to UCBSV and CBSV, respectively. The isolates of CBSV and UCBSV were 79.3-95.5% and 86.3-99.3 % identitical at nt level, respectively, suggesting more variation amongst CBSV isolates. The main sources of variation in plant viruses are mutations and recombination. Signals for recombination events were detected in 50% of isolates of each virus. Recombination events were detected in coding and non-coding (3’-UTR) sequences except in the 5’UTR and P3. There was no evidence for recombination between isolates of CBSV and UCBSV. The non-synonomous (dN) to synonomous (dS) nucleotide substitution ratio (ω) for the HAM1h and CP domains of both viruses were ≤ 0.184 suggesting that most sites of these proteins were evolving under strong purifying selection. However, there were individual amino acid sites that were submitted to adaptive evolution. For instance, adaptive evolution was detected in the HAM1h of UCBSV (n=15) where 12 aa sites were under positive selection (P< 0.05) but not in CBSV (n=12). The CP of CBSV (n=23) contained 12 aa sites (p<0.01) while only 5 aa sites in the CP gene of UCBSV were predicted to be submitted to positive selection pressure (p<0.01). The advantages offered by the aa sites under positive selection could not be established but occurrence of such sites in the terminal ends of UCBSV-HAMIh, for example, was interpreted as a requirement for proteolysis during polyprotein processing. Two different primer pairs that simultaneously detect UCBSV and CBSV isolates were developed in this study. They were used successfully to study distribution of CBSV, UCBSV and their mixed infections in Tanzania and Uganda. It was established that the two viruses co-infect cassava and that incidences of co-infection could be as high as 50% around Lake Victoria on the Tanzanian side. Furthermore, it was revealed for the first time that both UCBSV and CBSV were widely distributed in Eastern Africa. The primer pair was also used to confirm infection in a close relative of cassava, Manihot glaziovii (Müller Arg.) with CBSV. DNA barcoding of M. glaziovii was done by sequencing the matK gene. Two out of seven M. glaziovii from the coastal areas of Korogwe and Kibaha in north eastern Tanzania were shown to be infected by CBSV but not UCBSV isolates. Detection in M. glaziovii has an implication in control and management of CBSD as it is likely to serve as virus reservoir. This study has contributed to the understanding of evolution of CBSV and UCBSV, which cause CBSD epidemic in Eastern Africa. The detection tools developed in this work will be useful in plant breeding, verification of the phytosanitary status of materials in regional and international movement of germplasm, and in all diagnostic activities related to management of CBSD. Whereas there are still many issues to be resolved such as the function and biological significance of HAM1h and its origin, this work has laid a foundation upon which the studies on these aspects can be based.

Communication pathways from the anti-codon region to the aminoacylation site in Methionyl tRNA Synthetase: Molecular dynamics simulations and the structure network analysis

The enzymes of the family of tRNA synthetases perform their functions with high precision by synchronously recognizing the anticodon region and the aminoacylation region, which are separated by ?70 in space. This precision in function is brought about by establishing good communication paths between the two regions. We have modeled the structure of the complex consisting of Escherichia coli methionyl-tRNA synthetase (MetRS), tRNA, and the activated methionine. Molecular dynamics simulations have been performed on the modeled structure to obtain the equilibrated structure of the complex and the cross-correlations between the residues in MetRS have been evaluated. Furthermore, the network analysis on these simulated structures has been carried out to elucidate the paths of communication between the activation site and the anticodon recognition site. This study has provided the detailed paths of communication, which are consistent with experimental results. Similar studies also have been carried out on the complexes (MetRS + activated methonine) and (MetRS + tRNA) along with ligand-free native enzyme. A comparison of the paths derived from the four simulations clearly has shown that the communication path is strongly correlated and unique to the enzyme complex, which is bound to both the tRNA and the activated methionine. The details of the method of our investigation and the biological implications of the results are presented in this article. The method developed here also could be used to investigate any protein system where the function takes place through long-distance communication.

Conformational studies of the triphenylphosphazenyl side chain in cyclophosphazenes. I. Crystal and molecular structure of N3P3C15(NPPh3)

The title compound, C t8H~sC15NaP4, crystallizes in the monoclinic space group P2~/n with a = 20.14 (2), b = 8.69 (1), c = 14.92 (2) A, fl = 98.8 (3) ° , Z = 4. The structure was determined from visual data and refined to R = 0-069 for 1450 reflections. The cyclophosphazene ring is non-planar. The exocyclic NPPh 3 group exhibits type I conformation [R. A. Shaw (1975). Pure Appl. Chem. 44, 317-341] in which the N-P bond is perpendicular to the adjacent P-CI bond.

Molecular genetics of Usher syndrome -inherited deafness and blindness

Usher syndrome (USH) is an inherited blindness and deafness disorder with variable vestibular dysfunction. The syndrome is divided into three subtypes according to the progression and severity of clinical symptoms. The gene mutated in Usher syndrome type 3 (USH3), clarin 1 (CLRN1), was identified in Finland in 2001 and two mutations were identified in Finnish patients at that time. Prior to this thesis study, the two CLRN1 gene mutations were the only USH mutations identified in Finnish USH patients. To further clarify the Finnish USH mutation spectrum, all nine USH genes were studied. Seven mutations were identified: one was a previously known mutation in CLRN1, four were novel mutations in myosin VIIa (MYO7A) and two were a novel and a previously known mutation in usherin (USH2A). Another aim of this thesis research was to further study the structure and function of the CLRN1 gene, and to clarify the effects of mutations on protein function. The search for new splice variants resulted in the identification of eight novel splice variants in addition to the three splice variants that were already known prior to this study. Studies of the possible promoter regions for these splice variants showed the most active region included the 1000 bases upstream of the translation start site in the first exon of the main three exon splice variant. The 232 aa CLRN1 protein encoded by the main (three-exon) splice variant was transported to the plasma membrane when expressed in cultured cells. Western blot studies suggested that CLRN1 forms dimers and multimers. The CLRN1 mutant proteins studied were retained in the endoplasmic reticulum (ER) and some of the USH3 mutations caused CLRN1 to be unstable. During this study, two novel CLRN1 sequence alterations were identified and their pathogenicity was studied with cell culture protein expression. Previous studies with mice had shown that Clrn1 is expressed in mouse cochlear hair cells and spiral ganglion cells, but the expression profile in mouse retina remained unknown. The Clrn1 knockout mice display cochlear cell disruption/death, but do not have a retinal phenotype. The zebrafish, Danio rerio, clrn1 was found to be expressed in hair cells associated with hearing and balance. Clrn1 expression was also found in the inner nuclear layer (INL), photoreceptor layer and retinal pigment epithelium layer (RPE) of the zebrafish retina. When Clrn1 production was knocked down with injected morpholino oligonucleotides (MO) targeting Clrn1 translation or correct splicing, the zebrafish larvae showed symptoms similar to USH3 patients. These larvae had balance/hearing problems and reduced response to visual stimuli. The knowledge this thesis research has provided about the mutations in USH genes and the Finnish USH mutation spectrum are important in USH patient diagnostics. The extended information about the structure and function of CLRN1 is a step further in exploring USH3 pathogenesis caused by mutated CLRN1 as well as a step in finding a cure for the disease.

Ruthenium-Oxygen Coordination-Driven Self-Assembly of a Ru-8(II) Incomplete Prism: Synthesis, Structure, and Shape-Selective Molecular Recognition Study

Coordination-driven self-assembly of 1,3,5-benzenetricarboxylate (tma; 1) and oxalato-bridged p-cymeneruthenium(II) building block Ru-2(mu-eta(4)-C2O4)(MeOH)(2)(eta(6)-p-cymene)(2)](O3SCF3)(2) (2) affords an unusual octanuclear incomplete prism Ru-8(eta(6)-p-cymene)(8)(tma)(2)(mu-eta(4)-C2O4)(2)(OMe)(4)](O3SCF3)( 2) (3), which exhibits a remarkable shape-selective binding affinity for neutral phenolic compounds via hydrogen-bonding interactions (p-cymene = p-(PrC6H4Me)-Pr-i). Such a binding was confirmed by single-crystal X-ray diffraction analysis using 1,3,5-trihydroxybenzene as an analyte.

Frequency-comb-referenced molecular spectroscopy in the mid-infrared region

A simple method for absolute frequency measurements of molecular transitions in the mid-infrared region is reported. The method is based on a cw singly-resonant optical parametric oscillator (SRO), which is tunable from 3.2 to 3.45 µm. The mid- infrared frequency of the SRO is referenced to an optical frequency comb through its pump and signal beams. Sub-Doppler spectroscopy and absolute frequency measurement of the P(7) transition of the ν3 band of CH4 are demonstrated.

Molecular studies on LIM-domain protein CRP1

Various intrinsic and external factors are constantly attacking the cells causing damage to DNA and to other cellular structures. Cells in turn have evolved with different kinds of mechanisms to protect against the attacks and to repair the damage. Ultraviolet radiation (UVR) is one of the major environmental genotoxic carcinogens that causes inflammation, mutations, immunosuppression, accelerated aging of the skin and skin cancers. Epidermis is the outermost layer of the skin consisting mostly of keratinocytes, whose primary function is to protect the skin against e.g. UV radiation. LIM domain proteins are a group of proteins involved in regulation of cell growth, damage signalling, cell fate determination and signal transduction. Despite their two zinc fingers, LIM domains do not bind to DNA, but rather mediate protein-protein interactions and function as modular protein binding interfaces. We initially identified CSRP1 as UVR-regulated transcript by using expression profiling. Here we have further studied the regulation and function of CRP1, a representative of cysteine rich protein- family consisting of two LIM domains. We find that CRP1 is increased by UVR in primary human keratinocytes and in normal human skin fibroblasts. Ectopic expression of CRP1 protected the cells against UVR and provided a survival advantage, whereas silencing of CRP1 rendered the cells more photosensitive. Actinic keratosis is a premalignant lesion of skin caused by excess exposure to sunlight and sunburn, which may lead to formation of squamous cell carcinoma. The expression of CRP1 was increased in basal keratinocytes of Actinic keratosis patient specimens suggesting that CRP1 may be increased by constant exposure to UVR and may provide survival advantage for the cells also in vivo. In squamous cell carcinoma, CRP1 was only expressed in the fibroblasts surrounding the tumour. Moreover, we found that ectopic expression of CRP1 suppresses cell proliferation. Transforming growth factor beta (TGFbeta) is a multifunctional cytokine that regulates several functions in cell including growth, apoptosis and differentiation, and plays important roles in pathological disorders like cancer and fibrosis. We found that TGFbeta-signalling pathway regulates CRP1 at protein, but not at transcriptional level. The increase was mediated both through Smad and non-Smad signalling pathways involving MAPK/p38. Furthermore, we found that TGFbeta-mediated increase in CRP1 was associated with myofibroblast differentiation, and that CRP1 was significantly more expressed in idiopathic pulmonary fibrosis as compared to normal lung specimens. Since cell contractility is a distinct feature of myofibroblasts, and CRP1 is associated with actin cytoskeleton, we studied the role of CRP1 in cell contractility. CRP1 was found to localize to stress fibres that mediate contractility and to mediate myofibroblast contraction. These studies identify CRP1 as a stress responsive and cytokine regulated cytoskeletal protein that participates in pathological processes involved in fibrotic diseases and cancer.

Structural Relationships in 2,3-Bis-n-decyloxyanthracene and 12-Hydroxystearic Acid Molecular Gels and Aerogels Processed in Supercritical CO2

Supercritical carbon dioxide is used to prepare aerogels of two reference molecular organogelators, 2,3-bis-n-decyloxyanthracene (DDOA) (luminescent molecule) and 12-hydroxystearic acid (HSA). Electron microscopy reveals the fibrillar morphology of the aggregates generated by the protocol. SAXS and SANS measurements show that DDOA aerogels are crystalline materials exhibiting three morphs: (1) arrangements of the crystalline solid (2D p6m), (2) a second hexagonal morph slightly more compact, and (3) a packing specific of the fibers in the gel. Aggregates specific of the aerogel (volume fraction being typically phi approximate to 0.60) are developed over larger distances (similar to 1000 angstrom) and bear fewer defaults and residual strains than aggregates in the crystalline and gel phases. Porod, Scherrer and Debye-Bueche analyses of the scattering data have been performed. The first five diffraction peaks show small variations in position and intensity assigned to the variation of the number of fibers and their degree of vicinity within hexagonal bundles of the related SAFIN according to the Oster model. Conclusions are supported by the guidelines offered by the analysis of the situation in HSA aerogels for which the diffraction pattern can be described by two coexisting lamellar-like arrangements. The porosity of the aerogel, as measured by its specific surface extracted from the scattering invariant analysis, is only 1.8 times less than that of the swollen gel and is characteristic of a very porous material.

A nanosecond molecular dynamics study of antiparallel d(G)(7) quadruplex structures: Effect of the coordinated cations

Nanosecond scale molecular dynamics simulations have been performed on antiparallel Greek key type d(G(7)) quadruplex structures with different coordinated ions, namely Na+ and K+ ion, water and Na+ counter ions, using the AMBER force field and Particle Mesh Ewald technique for electrostatic interactions. Antiparallel structures are stable during the simulation, with root mean square deviation values of similar to1.5 Angstrom from the initial structures. Hydrogen bonding patterns within the G-tetrads depend on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate different cations. However, alternating syn-anti arrangement of bases along a chain as well as in a quartet is maintained through out the MD simulation. Coordinated Na+ ions, within the quadruplex cavity are quite mobile within the central channel and can even enter or exit from the quadruplex core, whereas coordinated K+ ions are quite immobile. MD studies at 400 K indicate that K+ ion cannot come out from the quadruplex core without breaking the terminal G-tetrads. Smaller grooves in antiparallel structures are better binding sites for hydrated counter ions, while a string of hydrogen bonded water molecules are observed within both the small and large grooves. The hydration free energy for the K+ ion coordinated structure is more favourable than that for the Na+ ion coordinated antiparallel quadruplex structure.

Effect of coordinated ions on structure and flexiblity of parallel G-quandruplexes: a molecular dynamics study

Single tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G7) quadruplex structures with Na+ or K+ ions coordinated in the cavity formed by the 06 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. These quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 A from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Na+ counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. In the absence of any coordinated ion, due to strong mutual repulsion, 06 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na+ ions, within the quadruplex cavity, are more mobile than coordinated K+ ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures

A molecular dynamics study of a PbO.SiO2 glass and melt

The structures of a PbO.SiO2 glass and melt have been studied using molecular dynamics simulation employing Born-Mayer-Huggins pair potentials. Various pair distribution functions are presented and discussed. Pb-Pb correlations persist in the melt, in agreement with experimental observations. The calculated and experimental radial distribution functions are compared.