Language change as an evolutionary process

In the thesis it is discussed in what ways concepts and methodology developed in evolutionary biology can be applied to the explanation and research of language change. The parallel nature of the mechanisms of biological evolution and language change is explored along with the history of the exchange of ideas between these two disciplines. Against this background computational methods developed in evolutionary biology are taken into consideration in terms of their applicability to the study of historical relationships between languages. Different phylogenetic methods are explained in common terminology, avoiding the technical language of statistics. The thesis is on one hand a synthesis of earlier scientific discussion, and on the other an attempt to map out the problems of earlier approaches in addition to finding new guidelines in the study of language change on their basis. Primarily literature about the connections between evolutionary biology and language change, along with research articles describing applications of phylogenetic methods into language change have been used as source material. The thesis starts out by describing the initial development of the disciplines of evolutionary biology and historical linguistics, a process which right from the beginning can be seen to have involved an exchange of ideas concerning the mechanisms of language change and biological evolution. The historical discussion lays the foundation for the handling of the generalised account of selection developed during the recent few decades. This account is aimed for creating a theoretical framework capable of explaining both biological evolution and cultural change as selection processes acting on self-replicating entities. This thesis focusses on the capacity of the generalised account of selection to describe language change as a process of this kind. In biology, the mechanisms of evolution are seen to form populations of genetically related organisms through time. One of the central questions explored in this thesis is whether selection theory makes it possible to picture languages are forming populations of a similar kind, and what a perspective like this can offer to the understanding of language in general. In historical linguistics, the comparative method and other, complementing methods have been traditionally used to study the development of languages from a common ancestral language. Computational, quantitative methods have not become widely used as part of the central methodology of historical linguistics. After the fading of a limited popularity enjoyed by the lexicostatistical method since the 1950s, only in the recent years have also the computational methods of phylogenetic inference used in evolutionary biology been applied to the study of early language history. In this thesis the possibilities offered by the traditional methodology of historical linguistics and the new phylogenetic methods are compared. The methods are approached through the ways in which they have been applied to the Indo-European languages, which is the most thoroughly investigated language family using both the traditional and the phylogenetic methods. The problems of these applications along with the optimal form of the linguistic data used in these methods are explored in the thesis. The mechanisms of biological evolution are seen in the thesis as parallel in a limited sense to the mechanisms of language change, however sufficiently so that the development of a generalised account of selection is deemed as possibly fruiful for understanding language change. These similarities are also seen to support the validity of using phylogenetic methods in the study of language history, although the use of linguistic data and the models of language change employed by these models are seen to await further development.

Genetic studies on recurrent miscarriage

Recurrent miscarriage (RM) is defined as three consecutive pregnancy failures and is estimated to affect ~1% of couples trying to conceive. The cause of RM remains unknown in approximately 50% of cases. In this study, it was hypothesized that some of the underlying factors yet to be discovered are genetic. The aim was to search for mutations in genes AMN, EPCR, TM, and p53 known to cause miscarriage in mouse models and thereby find new genetic causes for unexplained miscarriages in humans. In addition, the mitochondrial genome was studied because mitochondria are involved in processes important in early development. Furthermore, sex chromosome characteristics suggested to underlie miscarriage were also studied. A total of 40 couples and 8 women with unexplained RM were collected for this study and screened for mutations in the candidate genes. Six interesting exonic or potential splice site disrupting variations were detected. However, their phenotypic effects cannot be determined without further investigations. Additionally, an association between the C11992A polymorphism of the p53 gene and RM was detected. The results indicate that women carrying the C/A or A/A genotype have a two-fold higher risk for RM than women with a C/C genotype. This strengthens the results of previous studies reporting that p53 sequence variations may cause miscarriage. The role of variation C11992A in embryonic development is, however, difficult to predict without further studies When screening the mitochondrial genome a heteroplasmic mtDNA variation was found in an unexpected high number of women, as heteroplasmic variations are reported to be rare. One novel variation and 18 previously reported polymorphisms were detected in the mitochondrial genome. Although the detected variations are likely to be neutral polymorphisms, a role in the aetiology of miscarriage cannot be excluded as some mtDNA variations may be pathogenic only when a threshold is reached. Recent publications have reported skewed X chromosome inactivation and Y chromosome microdeletions to be associated with RM. Therefore, these sex chromosome abnormalities in the context of RM were investigated. No associations between skewed X chromosome inactivation or Y chromosome microdeletions and RM in the Finnish patients were detected. Data on ancestral birthplaces of the patients were collected to study any possible geographic clustering, which would indicate a common predisposing factor. The results showed clustering of the birthplaces in eastern Finland in a subset of patients. This suggests a possibility of an enriched susceptibility gene which may contribute to RM.

Computational Methods for Reconstruction and Analysis of Genome-Scale Metabolic Networks

Metabolism is the cellular subsystem responsible for generation of energy from nutrients and production of building blocks for larger macromolecules. Computational and statistical modeling of metabolism is vital to many disciplines including bioengineering, the study of diseases, drug target identification, and understanding the evolution of metabolism. In this thesis, we propose efficient computational methods for metabolic modeling. The techniques presented are targeted particularly at the analysis of large metabolic models encompassing the whole metabolism of one or several organisms. We concentrate on three major themes of metabolic modeling: metabolic pathway analysis, metabolic reconstruction and the study of evolution of metabolism. In the first part of this thesis, we study metabolic pathway analysis. We propose a novel modeling framework called gapless modeling to study biochemically viable metabolic networks and pathways. In addition, we investigate the utilization of atom-level information on metabolism to improve the quality of pathway analyses. We describe efficient algorithms for discovering both gapless and atom-level metabolic pathways, and conduct experiments with large-scale metabolic networks. The presented gapless approach offers a compromise in terms of complexity and feasibility between the previous graph-theoretic and stoichiometric approaches to metabolic modeling. Gapless pathway analysis shows that microbial metabolic networks are not as robust to random damage as suggested by previous studies. Furthermore the amino acid biosynthesis pathways of the fungal species Trichoderma reesei discovered from atom-level data are shown to closely correspond to those of Saccharomyces cerevisiae. In the second part, we propose computational methods for metabolic reconstruction in the gapless modeling framework. We study the task of reconstructing a metabolic network that does not suffer from connectivity problems. Such problems often limit the usability of reconstructed models, and typically require a significant amount of manual postprocessing. We formulate gapless metabolic reconstruction as an optimization problem and propose an efficient divide-and-conquer strategy to solve it with real-world instances. We also describe computational techniques for solving problems stemming from ambiguities in metabolite naming. These techniques have been implemented in a web-based sofware ReMatch intended for reconstruction of models for 13C metabolic flux analysis. In the third part, we extend our scope from single to multiple metabolic networks and propose an algorithm for inferring gapless metabolic networks of ancestral species from phylogenetic data. Experimenting with 16 fungal species, we show that the method is able to generate results that are easily interpretable and that provide hypotheses about the evolution of metabolism.

Mechanisms for alphavirus nonstructural polyprotein processing

Replication and transcription of the RNA genome of alphaviruses relies on a set of virus-encoded nonstructural proteins. They are synthesized as a long polyprotein precursor, P1234, which is cleaved at three processing sites to yield nonstructural proteins nsP1, nsP2, nsP3 and nsP4. All the four proteins function as constitutive components of the membrane-associated viral replicase. Proteolytic processing of P1234 polyprotein is precisely orchestrated and coordinates the replicase assembly and maturation. The specificity of the replicase is also controlled by proteolytic cleavages. The early replicase is composed of P123 polyprotein intermediate and nsP4. It copies the positive sense RNA genome to complementary minus-strand. Production of new plus-strands requires complete processing of the replicase. The papain-like protease residing in nsP2 is responsible for all three cleavages in P1234. This study addressed the mechanisms of proteolytic processing of the replicase polyprotein in two alphaviruses Semliki Forest virus (SFV) and Sindbis virus (SIN) representing different branches of the genus. The survey highlighted the functional relation of the alphavirus nsP2 protease to the papain-like enzymes. A new structural motif the Cys-His catalytic dyad accompanied with an aromatic residue following the catalytic His was described for nsP2 and a subset of other thiol proteases. Such an architecture of the catalytic center was named the glycine specificity motif since it was implicated in recognition of a specific Gly residue in the substrate. In particular, the presence of the motif in nsP2 makes the appearance of this amino acid at the second position upstream of the scissile bond a necessary condition for the cleavage. On top of that, there were four distinct mechanisms identified, which provide affinity for the protease and specifically direct the enzyme to different sites in the P1234 polyprotein. Three factors RNA, the central domain of nsP3 and the N-terminus of nsP2 were demonstrated to be external modulators of the nsP2 protease. Here I suggest that the basal nsP2 protease specificity is inherited from the ancestral papain-like enzyme and employs the recognition of the upstream amino acid signature in the immediate vicinity of the scissile bond. This mechanism is responsible for the efficient processing of the SFV nsP3/nsP4 junction. I propose that the same mechanism is involved in the cleavage of the nsP1/nsP2 junction of both viruses as well. However, in this case it rather serves to position the substrate, whereas the efficiency of the processing is ensured by the capability of nsP2 to cut its own N-terminus in cis. Both types of cleavages are demonstrated here to be inhibited by RNA, which is interpreted as impairing the basal papain-like recognition of the substrate. In contrast, processing of the SIN nsP3/nsP4 junction was found to be activated by RNA and additionally potentiated by the presence of the central region of nsP3 in the protease. The processing of the nsP2/nsP3 junction in both viruses occurred via another mechanism, requiring the exactly processed N-terminus of nsP2 in the protease and insensitive to RNA addition. Therefore, the three processing events in the replicase polyprotein maturation are performed via three distinct mechanisms in each of two studied alphaviruses. Distinct sets of conditions required for each cleavage ensure sequential maturation of P1234 polyprotein: nsP4 is released first, then the nsP1/nsP2 site is cut in cis, and liberation of the nsP2 N-terminus activates the cleavage of the nsP2/nsP3 junction at last. The first processing event occurs differently in SFV and SIN, whereas the subsequent cleavages are found to be similar in the two viruses and therefore, their mechanisms are suggested to be conserved in the genus. The RNA modulation of the alphavirus nonstructural protease activity, discovered here, implies bidirectional functional interplay between the alphavirus RNA metabolism and protease regulation. The nsP2 protease emerges as a signal transmitting moiety, which senses the replication stage and responds with proteolytic cleavages. A detailed hypothetical model of the alphavirus replicase core was inferred from the data obtained in the study. Similar principles in replicase organization and protease functioning are expected to be employed by other RNA viruses.

Human genetic variation in the Baltic Sea region: Features of population history and natural selection

In this thesis, the genetic variation of human populations from the Baltic Sea region was studied in order to elucidate population history as well as evolutionary adaptation in this region. The study provided novel understanding of how the complex population level processes of migration, genetic drift, and natural selection have shaped genetic variation in North European populations. Results from genome-wide, mitochondrial DNA and Y-chromosomal analyses suggested that the genetic background of the populations of the Baltic Sea region lies predominantly in Continental Europe, which is consistent with earlier studies and archaeological evidence. The late settlement of Fennoscandia after the Ice Age and the subsequent small population size have led to pronounced genetic drift, especially in Finland and Karelia but also in Sweden, evident especially in genome-wide and Y-chromosomal analyses. Consequently, these populations show striking genetic differentiation, as opposed to much more homogeneous pattern of variation in Central European populations. Additionally, the eastern side of the Baltic Sea was observed to have experienced eastern influence in the genome-wide data as well as in mitochondrial DNA and Y-chromosomal variation – consistent with linguistic connections. However, Slavic influence in the Baltic Sea populations appears minor on genetic level. While the genetic diversity of the Finnish population overall was low, genome-wide and Y-chromosomal results showed pronounced regional differences. The genetic distance between Western and Eastern Finland was larger than for many geographically distant population pairs, and provinces also showed genetic differences. This is probably mainly due to the late settlement of Eastern Finland and local isolation, although differences in ancestral migration waves may contribute to this, too. In contrast, mitochondrial DNA and Y-chromosomal analyses of the contemporary Swedish population revealed a much less pronounced population structure and a fusion of the traces of ancient admixture, genetic drift, and recent immigration. Genome-wide datasets also provide a resource for studying the adaptive evolution of human populations. This study revealed tens of loci with strong signs of recent positive selection in Northern Europe. These results provide interesting targets for future research on evolutionary adaptation, and may be important for understanding the background of disease-causing variants in human populations.

Electron cryo-microscopy studies of bacteriophage phi8 and archaeal virus SH1

Symmetry is a key principle in viral structures, especially the protein capsid shells. However, symmetry mismatches are very common, and often correlate with dynamic functionality of biological significance. The three-dimensional structures of two isometric viruses, bacteriophage phi8 and the archaeal virus SH1 were reconstructed using electron cryo-microscopy. Two image reconstruction methods were used: the classical icosahedral method yielded high resolution models for the symmetrical parts of the structures, and a novel asymmetric in-situ reconstruction method allowed us to resolve the symmetry mismatches at the vertices of the viruses. Evidence was found that the hexameric packaging enzyme at the vertices of phi8 does not rotate relative to the capsid. The large two-fold symmetric spikes of SH1 were found not to be responsible for infectivity. Both virus structures provided insight into the evolution of viruses. Comparison of the phi8 polymerase complex capsid with those of phi6 and other dsRNA viruses suggests that the quaternary structure in dsRNA bacteriophages differs from other dsRNA viruses. SH1 is unusual because there are two major types of capsomers building up the capsid, both of which seem to be composed mainly of single beta-barrels perpendicular to the capsid surface. This indicates that the beta-barrel may be ancestral to the double beta-barrel fold.

Dung beetle radiations in Madagascar

Speciation on islands is affected by island size and the range of habitats and resources available and often also by limited interactions with other taxa. An ancestral population may evolve into a large number of species via an adaptive radiation. In Madagascar, most groups of animals and plants have radiated on the island, having arrived via oceanic dispersal during the long isolation of Madagascar. Characteristic features of Malagasy biota are exceptionally high level of endemism, high species richness as well as lack of many higher taxa that are dominant on the African mainland. Malagasy dung beetles are dominated by two tribes, Canthonini and Helictopleurina, with more than 250 endemic species. In this thesis I have reconstructed molecular phylogenies for the two tribes using several gene regions and different phylogenetic methods. Evolution of closely related species and among populations of the same species was examined with haplotype networks. The Malagasy Canthonini consists of three large lineages, while Helictopleurina forms a monophyletic group. The ancestors of each of the four clades colonised Madagascar at different times during Cenozoic. The subsequent radiations differ in terms of the number of extant species (from 37 to more than 100) and the level of ecological differentiation. In addition, Onthophagini (6 species) and Scarabaeini (3) have colonised Madagascar several times, but they have not radiated and the few species have not entered forests where Canthonini and Helictopleurina mostly occur. Among the three Canthonini radiations, speciation appears to have been mostly allopatric in the oldest and the youngest clades, while in the Epactoides clade sister species have diverged in their ecologies but have similar geographical distributions, indicating that speciation may have occurred in regional sympatry. The most likely isolating mechanisms have been rivers and forest refugia during dry and cool geological periods. Most species are generalists feeding on both carrion and dung, and competition among ecologically similar species may prevent their coexistence in the same communities. Some species have evolved to forage in the canopy and a few species have shifted to use cattle dung, a new resource in the open habitats following the introduction of cattle 1500 years ago. The latter shift has allowed species to expand their geographical ranges.

Spectral Tuning and Adaptation to Different Light Environments of Mysid Visual Pigments

In the present thesis, questions of spectral tuning, the relation of spectral and thermal properties of visual pigments, and evolutionary adaptation to different light environments were addressed using a group of small crustaceans of the genus Mysis as a model. The study was based on microspectrophotometric measurements of visual pigment absorbance spectra, electrophysiological measurements of spectral sensitivities of dark-adapted eyes, and sequencing of the opsin gene retrieved through PCR. The spectral properties were related to the spectral transmission of the respective light environments, as well as to the phylogentic histories of the species. The photoactivation energy (Ea) was estimated from temperature effects on spectral sensitivity in the long-wavelength range, and calculations were made for optimal quantum catch and optimal signal-to-noise ratio in the different light environments. The opsin amino acid sequences of spectrally characterized individuals were compared to find candidate residues for spectral tuning. The general purpose was to clarify to what extent and on what time scale adaptive evolution has driven the functional properties of (mysid) visual pigments towards optimal performance in different light environments. An ultimate goal was to find the molecular mechanisms underlying the spectral tuning and to understand the balance between evolutionary adaptation and molecular constraints. The totally consistent segregation of absorption maxima (λmax) into (shorter-wavelength) marine and (longer-wavelength) freshwater populations suggests that truly adaptive evolution is involved in tuning the visual pigment for optimal performance, driven by selection for high absolute visual sensitivity. On the other hand, the similarity in λmax and opsin sequence between several populations of freshwater M. relicta in spectrally different lakes highlights the limits to adaptation set by evolutionary history and time. A strong inverse correlation between Ea and λmax was found among all visual pigments studied in these respects, including those of M. relicta and 10 species of vertebrate pigments, and this was used to infer thermal noise. The conceptual signal-to-noise ratios thus calculated for pigments with different λmax in the Baltic Sea and Lake Pääjärvi light environments supported the notion that spectral adaptation works towards maximizing the signal-to-noise ratio rather than quantum catch as such. Judged by the shape of absorbance spectra, the visual pigments of all populations of M. relicta and M. salemaai used exclusively the A2 chromophore (3, 4-dehydroretinal). A comparison of amino acid substitutions between M. relicta and M. salemaai indicated that mysid shrimps have a small number of readily available tuning sites to shift between a shorter - and a longer -wavelength opsin. However, phylogenetic history seems to have prevented marine M. relicta from converting back to the (presumably) ancestral opsin form, and thus the more recent reinvention of marine spectral sensitivity has been accomplished by some other novel mechanism, yet to be found

Major histocompatibility complex and coronary artery disease: Special emphasis on Chlamydia pneumoniae, and periodontitis

Most of the genes in the MHC region are involveed in adaptive and innate immunity, with essential function in inflammatory reactions and in protection against infections. These genes might serve as a candidate region for infection and inflammation associated diseases. CAD is an inflammatory disease. The present set of studies was performed to assess whether the MHC region harbors genetic markers for CAD, and whether these genetic markers explain the CAD risk factors: e.g., C. pneumoniae, periodontitis, and periodontal pathogens. Study I was performed using two separate patient materials and age- and sex-matched healthy controls, categorizing them into two independent studies: the HTx and ACS studies. Both studies consistently showed the HLA-A3– B35– DR1 (35 ancestral haplotype) haplotype as a susceptible MHC genetic marker for CAD. HLA-DR1 alone was associated not only with CAD, but also with CAD risk factor diseases, e.g., diabetes mellitus, and hyperlipidemia. The ACS study further showed the HLA-B*07 and -DRB1*15 -related haplotype as a protective MHC haplotype for CAD. Study II showed that patients with CAD showed signs of chronic C. pneumoniae infection when compared to age- and sex-matched healthy controls. HLA-B*35 or -related haplotypes associated with the C. pneumoniae infection markers. Among these haplotype carriers, males and smokers associated with elevated C. pneumoniae infection markers. Study III showed that CAD patients with periodontitis had elevated serum markers of P. gingivalis and occurrence of the pathogen in saliva. LTA+496C strongly associated with periodontitis, while HLA-DRB1*01 with periodontitis and with the elevated serum antibodies of P. gingivalis. Study IV showed that the increased level of C3/C4 ratio was a new risk factor and was associated with recurrent cardiovascular end-points. The increased C3 and decreased C4 concentrations in serum explained the increased level of the C3/C4 ratio. Both the higher than cut-off value (4.53) and the highest quartile of the C3/C4 ratio were also associated with worst survival, increased end-points, and C4 null alleles. The presence of C4 null alleles associated with decreased serum C4 concentration, and increased C3/C4 ratio. In conclusion, the present studies show that the CAD susceptibility haplotype (HLA-A3− B35− DR1 -related haplotypes, Study I) partially explains the development of CAD in patients possessing several recognized and novel risk factors: diabetes mellitus, increased LDL, smoking, C4B*Q0, C. pneumnoiae, periodontitis, P. gingivalis, and complement C3/C4 ratio (Study II, III, and IV).

Evolution of the Neckeraceae (Bryophyta): resolving the backbone phylogeny

Earlier phylogenetic studies, including species belonging to the Neckeraceae, have indicated that this pleurocarpous moss family shares a strongly supported sister group relationship with the Lembophyllaceae, but the family delimitation of the former needs adjustment. To test the monophyly of the Neckeraceae, as well as to redefine the family circumscription and to pinpoint its phylogenetic position in a larger context, a phylogenetic study based on molecular data was carried out. Sequence data were compiled, combining data from all three genomes: nuclear ITS1 and 2, plastid trnS-rps4-trnT-trnL-trnF and rpl16, and mitochondrial nad5 intron. The Neckeraceae have sometimes been divided into the two families, Neckeraceae and Thamnobryaceae, a division rejected here. Both parsimony and Bayesian analyses of molecular data revealed that the family concept of the Neckeraceae needs several further adjustments, such as the exclusion of some individual species and smaller genera as well as the inclusion of the Leptodontaceae. Within the family three well-supported clades (A, B and C) can be distinguished. Members of clade A are mainly non-Asiatic and nontropical. Most species have a weak costa and immersed capsules with reduced peristomes (mainly Neckera spp.) and the teeth at the leaf margins are usually unicellular. Clade B members are also mainly non-Asiatic. They are typically fairly robust, distinctly stipilate, having a single, at least relatively strong costa, long setae (capsules exserted), and the peristomes are well developed or only somewhat reduced. Members of clade C are essentially Asiatic and tropical. The species of this clade usually have a strong costa and a long seta, the seta often being mammillose in its upper part. The peristome types in this clade are mixed, since both reduced and unreduced types are found. Several neckeraceous genera that were recognised on a morphological basis are polyphyletic (e.g. Neckera, Homalia, Thamnobryum, Porotrichum). Ancestral state reconstructions revealed that currently used diagnostic traits, such as the leaf asymmetry and costa strength are highly homoplastic. Similarly, the reconstructions revealed that the 'reduced' sporophyte features have evolved independently in each of the three clades.

Phospholipid cytochrome c interactions

The ability of the peripherally associated membrane protein cytochrome c (cyt c) to bind phospholipids in vitro was studied using fluorescence spectroscopy and large unilamellar liposomes. Previous work has shown that cyt c can bind phospholipids using two distinct mecha- nisms and sites, the A-site and the C-site. This binding is mediated by electrostatic or hydrophobic interactions, respectively. Here, we focus on the mechanism underlying these interactions. A chemically modified cyt c mutant Nle91 was used to study the ATP-binding site, which is located near the evolutionarily invariant Arg 91 on the protein surface. This site was also demonstrated to mediate phospholipid binding, possibly by functioning as a phospholipid binding site. Circular dichroism spectroscopy, time resolved fluorescence spectroscopy of zinc- porphyrin modified [Zn2+-heme] cyt c and liposome binding studies of the Nle91 mutant were used to demonstrate that ATP induces a conformational change in membrane- bound cyt c. The ATP-induced conformational changes were mediated by Arg 91 and were most pronounced in cyt c bound to phospholipids via the C-site. It has been previously reported that the hydrophobic interaction between phospho- lipids and cyt c (C-site) includes the binding of a phospholipid acyl chain inside the protein. In this mechanism, which is known as extended phospholipid anchorage, the sn-2 acyl chain of a membrane phospholipid protrudes out of the membrane surface and is able to bind in a hydrophobic cavity in cyt c. Direct evidence for this type of bind- ing mechanism was obtained by studying cyt c/lipid interaction using fluorescent [Zn2+- heme] cyt c and fluorescence quenching of brominated fatty acids and phospholipids. Under certain conditions, cyt c can form fibrillar protein-lipid aggregates with neg- atively charged phospholipids. These aggregates resemble amyloid fibrils, which are involved in the pathogenesis of many diseases. Congo red staining of these fibers con- firmed the presence of amyloid structures. A set of phospholipid-binding proteins was also found to form similar aggregates, suggesting that phospholipid-induced amyloid formation could be a general mechanism of amyloidogenesis.

Hiekkatympösen (Hebeloma cylindrosporum) L-aminohappo-oksidaasi -geenin eristäminen

Pohjoisella havumetsävyöhykkeellä typpi on usein kasvien kasvua rajoittava tekijä. Metsämaan typpivarannot koostuvat pääasiassa orgaaniseen ainekseen sitoutuneista typpiyhdisteistä, erityisesti aminohapoista. Ektomykorritsasienet osallistuvat metsämaassa tapahtuvaan typenkiertoon hajottamalla orgaanisia typpiyhdisteitä ja kuljettamalla niitä kasvien käytettäväksi. Sienisolun sisällä tapahtuvasta aminohappojen mineralisaatiosta tiedetään toistaiseksi melko vähän. Aminohappo-oksidaasit katalysoivat aminohappojen mineralisaatiota. Eräissä ektomykorritsaa muodostavien kantasienten suvuissa on osoitettu L-aminohappo-oksidaaseja (LAO). Toistaiseksi LAO-geeniä ei tunneta kantasienistä. Työssä kuvattiin ensimmäistä kertaa LAO-geeni kantasienistä. Hiekkatympösen LAO1- geenin cDNA:n 5´ ja 3´ päiden emäsjärjestykset määritettiin RACE-PCR -menetelmällä, josta saatujen sekvenssien perusteella suunniteltiin alukkeet koko geenin cDNA:n ja genomisen DNA:n monistamiseksi. Genomisen DNA ja cDNA -sekvenssien perusteella määritettiin hiekkatympösen LAO1-geenin rakenne. Hiekkatympösen LAO1-geeni koostuu viidestä eksonista ja neljästä intronista. Hiekkatympösen LAO1-geenin yläpuoliselta alueelta löydettiin typpimetabolian säätelyyn osallistuvan proteiinin sitoutumiskohta. LAO1-geeniä edeltävä geenin osittainen genominen DNA-sekvenssi määritettiin. Kangaslohisienen genomissa LAO1-geeniä edeltävä geeni oli ennustettu pyruvaattidekarboksylaasiksi. Lisäksi työssä määritettiin hiekkatympösen toisen LAOhomologin cDNA:n osittainen emäsjärjestys. Työssä tunnistettiin myös toisen kantasienen, kangaslohisienen, LAO-geeni. LAO-geeniksi tunnistettu kangaslohisienen geenimalli oli aiemmin ennustettu NCBI:n tietokannassa toiminnaltaan tuntemattomaksi proteiiniksi. Proteiinien sukupuun perusteella hiekkatympösen ja kangaslohisienen LAO:n kantamuoto on kahdentunut. Työstä saatu tutkimustulos tuo täysin uutta tietoa molekyylibiologian tasolla ektomykorritsasienten aminohappojen katabolisista reaktioista. Aminohappojen mineralisaation seurauksen muodostuneet ammoniumionit saattavat olla merkittävä typen lähde myös maan muille mikrobeille ja kasveille. On mahdollista, että ektomykorritsasienten LAO-entsyymi on yksi merkittävä tekijä metsämaan typenkierrossa.

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.

Phylogeography and demographic history of Lacerta lepida in the Iberian Peninsula

The Iberian Peninsula is recognized as an important refugial area for species survival and diversification during the climatic cycles of the Quaternary. Recent phylogeographic studies have revealed Iberia as a complex of multiple refugia. However, most of these studies have focused either on species with narrow distributions within the region or species groups that, although widely distributed, generally have a genetic structure that relates to pre-Quaternary cladogenetic events. In this study we undertake a detailed phylogeographic analysis of the lizard species, Lacerta lepida, whose distribution encompasses the entire Iberian Peninsula. We attempt to identify refugial areas, recolonization routes, zones of secondary contact and date demographic events within this species. Results support the existence of 6 evolutionary lineages (phylogroups) with a strong association between genetic variation and geography, suggesting a history of allopatric divergence in different refugia. Diversification within phylogroups is concordant with the onset of the Pleistocene climatic oscillations. The southern regions of several phylogroups show a high incidence of ancestral alleles in contrast with high incidence of recently derived alleles in northern regions. All phylogroups show signs of recent demographic and spatial expansions. We have further identified several zones of secondary contact, with divergent mitochondrial haplotypes occurring in narrow zones of sympatry. The concordant patterns of spatial and demographic expansions detected within phylogroups, together with the high incidence of ancestral haplotypes in southern regions of several phylogroups, suggests a pattern of contraction of populations into southern refugia during adverse climatic conditions from which subsequent northern expansions occurred. This study supports the emergent pattern of multiple refugia within Iberia but adds to it by identifying a pattern of refugia coincident with the southern distribution limits of individual evolutionary lineages. These areas are important in terms of long-term species persistence and therefore important areas for conservation.