DNA structural variability as a factor in gene expression and evolution

Redundant DNA can buffer sequence dependent structural deviations from an ideal double helix. Buffering serves a mechanistic function by reducing extraneous conformational effects which could interfere with readout or which would impose energetic constraints on evolution. It also serves an evolutionary function by allowing for gradual variations in conformation-dependent regulation of gene expression. Such gradualism is critical for the rate of evolution. The buffer structure concept provides a new interpretation for repetitive DNA and for exons and introns.

Molecular dynamics simulations of sliding in an Fe-Cu tribopair system

Experimental evidence suggests that high strain rates, stresses, strains and temperatures are experienced near sliding interfaces. The associated microstructural changes are due to several dynamic an interacting phenomena. 3D non-equilibrium molecular dynamics (MD) simulations of sliding were conducted with the aim of understanding the dynamic processes taking place in crystalline tribopairs, with a focus on plastic deformation and microstructural evolution. Embedded atom potentials were employed for simulating sliding of an Fe-Cu tribopair. Sliding velocity, crystal orientation and presence of lattice defects were some of the variables in these simulations. Extensive plastic deformation involving dislocation and twin activity, dynamic recrystallization, amorphization and/or nanocrystallization, mechanical mixing and material transfer were observed. Mechanical mixing in the vicinity of the sliding interface was observed even in the Fe-Cu system, which would cluster under equilibrium conditions, hinting at the ballistic nature of the process. Flow localization was observed at high velocities implying the possible role of adiabatic heating. The presence of preexisting defects (such as dislocations and interfaces) played a pivotal role in determining friction and microstructural evolution. The study also shed light on the relationship between adhesion and plastic deformation, and friction. Comparisons with experiments suggest that such simulations can indeed provide valuable insights that are difficult to obtain from experiments.

Molecular characterisation of dicot-infecting mastreviruses from Australia

Monocotyledonous and dicotyledonous plant infecting mastreviruses threaten various agricultural systems throughout Africa, Eurasia and Australasia. In Australia three distinct mastrevirus species are known to infect dicotyledonous hosts such as chickpea, bean and tobacco. Amongst 34 new "dicot-infecting" mastrevirus full genome sequences obtained from these hosts we discovered one new species, four new strains, and various variants of previously described mastrevirus species. Besides providing additional support for the hypothesis that evolutionary processes operating during dicot-infecting mastrevirus evolution (such as patterns of pervasive homologous and non-homologous recombination, and strong purifying selection acting on all genes) have mostly mirrored those found in their monocot-infecting counterparts, we find that the Australian dicot-infecting viruses display patterns of phylogeographic clustering reminiscent of those displayed by monocot infecting mastrevirus species such as Panicum streak virus and Maize streak virus.

Evolution and detection of cyanobacterial hepatotoxin synthetase genes

Mass occurrences (blooms) of cyanobacteria are common in aquatic environments worldwide. These blooms are often toxic, due to the presence of hepatotoxins or neurotoxins. The most common cyanobacterial toxins are hepatotoxins: microcystins and nodularins. In freshwaters, the main producers of microcystins are Microcystis, Anabaena, and Planktothrix. Nodularins are produced by strains of Nodularia spumigena in brackish waters. Toxic and nontoxic strains of cyanobacteria co-occur and cannot be differentiated by conventional microscopy. Molecular biological methods based on microcystin and nodularin synthetase genes enable detection of potentially hepatotoxic cyanobacteria. In the present study, molecular detection methods for hepatotoxin-producing cyanobacteria were developed, based on microcystin synthetase gene E (mcyE) and the orthologous nodularin synthetase gene F (ndaF) sequences. General primers were designed to amplify the mcyE/ndaF gene region from microcystin-producing Anabaena, Microcystis, Planktothrix, and Nostoc, and nodularin-producing Nodularia strains. The sequences were used for phylogenetic analyses to study how cyanobacterial mcy genes have evolved. The results showed that mcy genes and microcystin are very old and were already present in the ancestor of many modern cyanobacterial genera. The results also suggested that the sporadic distribution of biosynthetic genes in modern cyanobacteria is caused by repeated gene losses in the more derived lineages of cyanobacteria and not by horizontal gene transfer. Phylogenetic analysis also proposed that nda genes evolved from mcy genes. The frequency and composition of the microcystin producers in 70 lakes in Finland were studied by conventional polymerase chain reaction (PCR). Potential microcystin producers were detected in 84% of the lakes, using general mcyE primers, and in 91% of the lakes with the three genus-specific mcyE primers. Potential microcystin-producing Microcystis were detected in 70%, Planktothrix in 63%, and Anabaena in 37% of the lakes. The presence and co-occurrence of potential microcystin producers were more frequent in eutrophic lakes, where the total phosphorus concentration was high. The PCR results could also be associated with various environmental factors by correlation and regression analyses. In these analyses, the total nitrogen concentration and pH were both associated with the presence of multiple microcystin-producing genera and partly explained the probability of occurrence of mcyE genes. In general, the results showed that higher nutrient concentrations increased the occurrence of potential microcystin producers and the risk for toxic bloom formation. Genus-specific probe pairs for microcystin-producing Anabaena, Microcystis, Planktothrix, and Nostoc, and nodularin-producing Nodularia were designed to be used in a DNA-chip assay. The DNA-chip can be used to simultaneously detect all these potential microcystin/nodularin producers in environmental water samples. The probe pairs detected the mcyE/ndaF genes specifically and sensitively when tested with cyanobacterial strains. In addition, potential microcystin/nodularin producers were identified in lake and Baltic Sea samples by the DNA-chip almost as sensitively as by quantitative real-time PCR (qPCR), which was used to validate the DNA-chip results. Further improvement of the DNA-chip assay was achieved by optimization of the PCR, the first step in the assay. Analysis of the mcy and nda gene clusters from various hepatotoxin-producing cyanobacteria was rewarding; it revealed that the genes were ancient. In addition, new methods detecting all the main producers of hepatotoxins could be developed. Interestingly, potential microcystin-producing cyanobacterial strains of Microcystis, Planktothrix, and Anabaena, co-occurred especially in eutrophic and hypertrophic lakes. Protecting waters from eutrophication and restoration of lakes may thus decrease the prevalence of toxic cyanobacteria and the frequency of toxic blooms.

The rph1 Gene Is a Common Contributor to the Evolution of Phosphine Resistance in Independent Field Isolates of Rhyzopertha Dominica

Phosphine is the only economically viable fumigant for routine control of insect pests of stored food products, but its continued use is now threatened by the world-wide emergence of high-level resistance in key pest species. Phosphine has a unique mode of action relative to well-characterised contact pesticides. Similarly, the selective pressures that lead to resistance against field sprays differ dramatically from those encountered during fumigation. The consequences of these differences have not been investigated adequately. We determine the genetic basis of phosphine resistance in Rhyzopertha dominica strains collected from New South Wales and South Australia and compare this with resistance in a previously characterised strain from Queensland. The resistance levels range from 225 and 100 times the baseline response of a sensitive reference strain. Moreover, molecular and phenotypic data indicate that high-level resistance was derived independently in each of the three widely separated geographical regions. Despite the independent origins, resistance was due to two interacting genes in each instance. Furthermore, complementation analysis reveals that all three strains contain an incompletely recessive resistance allele of the autosomal rph1 resistance gene. This is particularly noteworthy as a resistance allele at rph1 was previously proposed to be a necessary first step in the evolution of high-level resistance. Despite the capacity of phosphine to disrupt a wide range of enzymes and biological processes, it is remarkable that the initial step in the selection of resistance is so similar in isolated outbreaks.

The systematics of Endoraecium in Australia

Endoraecium (Raveneliaceae, Pucciniales) is a genus of rust that infects several species of Acacia (Fabaceae) in Australia, south-east Asia and Hawaii. Thirteen species of Endoraecium have been described, including seven species that are endemic to Australia, one species to south-east Asia and five to Hawaii. This study investigated the systematics of Endoraecium from 50 specimens in Australia and south-east Asia with a combined morphological and molecular approach. Phylogenetic analyses were conducted on combined datasets of the SSU, ITS and LSU regions of rDNA. The recovered phylogeny (i) supported a recent division of Endoraecium digitatum into five separate species based on morphology and host specificity and (ii) found lineages that did not correspond with known species.

Molecular diversity of Chickpea chlorotic dwarf virus in Sudan: High rates of intra-species recombination – a driving force in the emergence of new strains

In Sudan Chickpea chlorotic dwarf virus (CpCDV, genus Mastrevirus, family Geminiviridae) is an important pathogen of pulses that are grown both for local consumption, and for export. Although a few studies have characterised CpCDV genomes from countries in the Middle East, Africa and the Indian subcontinent, little is known about CpCDV diversity in any of the major chickpea production areas in these regions. Here we analyse the diversity of 146 CpCDV isolates characterised from pulses collected across the chickpea growing regions of Sudan. Although we find that seven of the twelve known CpCDV strains are present within the country, strain CpCDV-H alone accounted for ∼73% of the infections analysed. Additionally we identified four new strains (CpCDV-M, -N, -O and -P) and show that recombination has played a significant role in the diversification of CpCDV, at least in this region. Accounting for observed recombination events, we use the large amounts of data generated here to compare patterns of natural selection within protein coding regions of CpCDV and other dicot-infecting mastrevirus species.

Characterization of the molecular components and function of BARE-1, Hin-Mu and Mu transposition machineries

Transposable elements, transposons, are discrete DNA segments that are able to move or copy themselves from one locus to another within or between their host genome(s) without a requirement for DNA homology. They are abundant residents in virtually all the genomes studied, for instance, the genomic portion of TEs is approximately 3% in Saccharomyces cerevisiae, 45% in humans, and apparently more than 70% in some plant genomes such as maize and barley. Transposons plays essential role in genome evolution, in lateral transfer of antibiotic resistance genes among bacteria and in life cycle of certain viruses such as HIV-1 and bacteriophage Mu. Despite the diversity of transposable elements they all use a fundamentally similar mechanism called transpositional DNA recombination (transposition) for the movement within and between the genomes of their host organisms. The DNA breakage and joining reactions that underlie their transposition are chemically similar in virtually all known transposition systems. The similarity of the reactions is also reflected in the structure and function of the catalyzing enzymes, transposases and integrases. The transposition reactions take place within the context of a transposition machinery, which can be particularly complex, as in the case of the VLP (virus like particle) machinery of retroelements, which in vivo contains RNA or cDNA and a number of element encoded structural and catalytic proteins. Yet, the minimal core machinery required for transposition comprises a multimer of transposase or integrase proteins and their binding sites at the element DNA ends only. Although the chemistry of DNA transposition is fairly well characterized, the components and function of the transposition machinery have been investigated in detail for only a small group of elements. This work focuses on the identification, characterization, and functional studies of the molecular components of the transposition machineries of BARE-1, Hin-Mu and Mu. For BARE-1 and Hin-Mu transpositional activity has not been shown previously, whereas bacteriophage Mu is a general model of transposition. For BARE-1, which is a retroelement of barley (Hordeum vulgare), the protein and DNA components of the functional VLP machinery were identified from cell extracts. In the case of Hin-Mu, which is a Mu-like prophage in Haemophilus influenzae Rd genome, the components of the core machinery (transposase and its binding sites) were characterized and their functionality was studied by using an in vitro methodology developed for Mu. The function of Mu core machinery was studied for its ability to use various DNA substrates: Hin-Mu end specific DNA substrates and Mu end specific hairpin substrates. The hairpin processing reaction by MuA was characterized in detail. New information was gained of all three machineries. The components or their activity required for functional BARE-1 VLP machinery and retrotransposon life cycle were present in vivo and VLP-like structures could be detected. The Hin-Mu core machinery components were identified and shown to be functional. The components of the Mu and Hin-Mu core machineries were partially interchangeable, reflecting both evolutionary conservation and flexibility within the core machineries. The Mu core machinery displayed surprising flexibility in substrate usage, as it was able to utilize Hin-Mu end specific DNA substrates and to process Mu end DNA hairpin substrates. This flexibility may be evolutionarily and mechanistically important.

trans gene regulation in adaptive evolution: A genetic algorithm model

This is a continuation of earlier studies on the evolution of infinite populations of haploid genotypes within a genetic algorithm framework. We had previously explored the evolutionary consequences of the existence of indeterminate—“plastic”—loci, where a plastic locus had a finite probability in each generation of functioning (being switched “on”) or not functioning (being switched “off”). The relative probabilities of the two outcomes were assigned on a stochastic basis. The present paper examines what happens when the transition probabilities are biased by the presence of regulatory genes. We find that under certain conditions regulatory genes can improve the adaptation of the population and speed up the rate of evolution (on occasion at the cost of lowering the degree of adaptation). Also, the existence of regulatory loci potentiates selection in favour of plasticity. There is a synergistic effect of regulatory genes on plastic alleles: the frequency of such alleles increases when regulatory loci are present. Thus, phenotypic selection alone can be a potentiating factor in a favour of better adaptation.

Sostdc1 Plays an Essential Role in Mammalian Tooth Patterning : Insight into the Rodent Dental Evolution

This thesis work focuses on the role of TGF-beta family antagonists during the development of mouse dentition. Tooth develops through an interaction between the dental epithelium and underlying neural crest derived mesenchyme. The reciprocal signaling between these tissues is mediated by soluble signaling molecules and the balance between activatory and inhibitory signals appears to be essential for the pattern formation. We showed the importance of Sostdc1 in the regulation of tooth shape and number. The absence of Sostdc1 altered the molar cusp patterning and led to supernumerary tooth formation both in the molar and incisor region. We showed that initially, Sostdc1 expression is in the mesenchyme, suggesting that dental mesenchyme may limit supernumerary tooth induction. We tested this in wild-type incisors by minimizing the amount of mesenchymal tissue surrounding the incisor tooth germs prior to culture in vitro. The cultured teeth phenocopied the extra incisor phenotype of the Sostdc1-deficient mice. Furthermore, we showed that minimizing the amount of dental mesenchyme in cultured Sostdc1-deficient incisors caused the formation of additional de novo incisors that resembled the successional incisor development resulting from activated Wnt signaling. Sostdc1 seemed to be able to inhibit both mesenchymal BMP4 and epithelial canonical Wnt signaling, which thus allows Sostdc1 to restrict the enamel knot size and regulate the tooth shape and number. Our work emphasizes the dual role for the tooth mesenchyme as a suppressor as well as an activator during tooth development. We found that the placode, forming the thick mouse incisor, is prone to disintegration during initiation of tooth development. The balance between two mesenchymal TGF-beta family signals, BMP4 and Activin is essential in this regulation. The inhibition of BMP4 or increase in Activin signaling led to the splitting of the large incisor placode into two smaller placodes resulting in thin incisors. These two signals appeared to have different effects on tooth epithelium and the analysis of the double null mutant mice lacking Sostdc1 and Follistatin indicated that these TGF-beta inhibitors regulate the mutual balance of BMP and Activin in vivo. In addition, this work provides an alternative explanation for the issue of incisor identity published in Science by Tucker et al. in 1998 and proposes that the molar like morphology that can be obtained by inhibiting BMP signaling is due to partial splitting of the incisor placodes and not due to change in tooth identity from the incisor to the molar. This thesis work presents possible molecular mechanisms that may have modified the mouse dental pattern during evolution leading to the typical rodent dentition of modern mouse. The rodent dentition is specialized for gnawing and consists of two large continuously growing incisors and toothless diastema region separating the molars and incisors. The ancestors of rodents had higher number of more slender incisors together with canines and premolars. Additionally, murine rodents, which include the mouse, have lost their ability for tooth replacement. This work has revealed that the inhibitory molecules appear to play a role in the tooth number suppression by delineating the spatial and temporal action of the inductive signals. The results suggest that Sostdc1 plays an essential role in several stages of tooth development through the regulation of both the BMP and Wnt pathway. The work shows a dormant sequential tooth forming potential present in wild type mouse incisor region and gives a new perspective on tooth suppression by dental mesenchyme. It reveals as well a novel mechanism to create a large mouse incisor through the regulation of mesenchymal balance between inductive and inhibitory signals.

Ab initio molecular dynamics simulations of excited hydrogen halides and methyl halides

Using excited-state ab initio molecular dynamics simulations employing the complete-active-space self-consistent-field approach, we study the mechanism of photodissociation in terms of time evolution of structure, kinetic energy, charges and potential energy for the first excited state of hydrogen halides and methyl halides. Although the hydrogen halides and methyl halides are similar in the photodissociation mechanism, their dynamics are slightly different. The presence of the methyl group causes delay in photodissociation as compared to hydrogen halides.

Molecular epidemiology of human rhinoviruses

The first part of this work investigates the molecular epidemiology of a human enterovirus (HEV), echovirus 30 (E-30). This project is part of a series of studies performed in our research team analyzing the molecular epidemiology of HEV-B viruses. A total of 129 virus strains had been isolated in different parts of Europe. The sequence analysis was performed in three different genomic regions: 420 nucleotides (nt) in the VP4/VP2 capsid protein coding region, the entire VP1 capsid protein coding gene of 876 nt, and 150 nt in the VP1/2A junction region. The analysis revealed a succession of dominant sublineages within a major genotype. The temporally earlier genotypes had been replaced by a genetically homogenous lineage that has been circulating in Europe since the late 1970s. The same genotype was found by other research groups in North America and Australia. Globally, other cocirculating genetic lineages also exist. The prevalence of a dominant genotype makes E-30 different from other previously studied HEVs, such as polioviruses and coxsackieviruses B4 and B5, for which several coexisting genetic lineages have been reported. The second part of this work deals with molecular epidemiology of human rhinoviruses (HRVs). A total of 61 field isolates were studied in the 420-nt stretch in the capsid coding region of VP4/VP2. The isolates were collected from children under two years of age in Tampere, Finland. Sequences from the clinical isolates clustered in the two previously known phylogenetic clades. Seasonal clustering was found. Also, several distinct serotype-like clusters were found to co-circulate during the same epidemic season. Reappearance of a cluster after disappearing for a season was observed. The molecular epidemiology of the analyzed strains turned out to be complex, and we decided to continue our studies of HRV. Only five previously published complete genome sequences of HRV prototype strains were available for analysis. Therefore, all designated HRV prototype strains (n=102) were sequenced in the VP4/VP2 region, and the possibility of genetic typing of HRV was evaluated. Seventy-six of the 102 prototype strains clustered in HRV genetic group A (HRV-A) and 25 in group B (HRV-B). Serotype 87 clustered separately from other HRVs with HEV species D. The field strains of HRV represented as many as 19 different genotypes, as judged with an approximate demarcation of a 20% nt difference in the VP4/VP2 region. The interserotypic differences of HRV were generally similar to those reported between different HEV serotypes (i.e. about 20%), but smaller differences, less than 10%, were also observed. Because some HRV serotypes are genetically so closely related, we suggest that the genetic typing be performed using the criterion "the closest prototype strain". This study is the first systematic genetic characterization of all known HRV prototype strains, providing a further taxonomic proposal for classification of HRV. We proposed to divide the genus Human rhinoviruses into HRV-A and HRV-B. The final part of the work comprises a phylogenetic analysis of a subset (48) of HRV prototype strains and field isolates (12) in the nonstructural part of the genome coding for the RNA-dependent RNA polymerase (3D). The proposed division of the HRV strains in the species HRV-A and HRV-B was also supported by 3D region. HRV-B clustered closer to HEV species B, C, and also to polioviruses than to HRV-A. Intraspecies variation within both HRV-A and HRV-B was greater in the 3D coding region than in the VP4/VP2 coding region, in contrast to HEV. Moreover, the diversity of HRV in 3D exceeded that of HEV. One group of HRV-A, designated HRV-A', formed a separate cluster outside other HRV-A in the 3D region. It formed a cluster also in the capsid region, but located within HRV-A. This may reflect a different evolutionary history of distinct genomic regions among HRV-A. Furthermore, the tree topology within HRV-A in the 3D region differed from that in the VP4/VP2, suggesting possible recombination events in the evolution of the strains. No conflicting phylogenies were observed in any of the 12 field isolates. Possible recombination was further studied using the Similarity and Bootscanning analyses of the complete genome sequences of HRV available in public databases. Evidence for recombination among HRV-A was found, as HRV2 and HRV39 showed higher similarity in the nonstructural part of the genome. Whether HRV2 and HRV39 strains - and perhaps also some other HRV-A strains not yet completely sequenced - are recombinants remains to be determined.