Phylogeography and Adaptive Divergence of Three-spined Stickleback Populations

Genetic studies on phylogeography and adaptive divergence in Northern Hemisphere fish species such as three-spined stickleback (Gasterosteus aculeatus) provide an excellent opportunity to investigate genetic mechanisms underlying population differentiation. According to the theory, the process of population differentiation results from a complex interplay between random and deterministic processes as well historical factors. The main scope in this thesis was to study how historical factors like the Pleistocene ice ages have shaped the patterns molecular diversity in three-spined stickleback populations in Europe and how this information could be utilized in the conservation genetic context. Furthermore, identifying footprints of natural selection at the DNA level might be used in identifying genes involved in evolutionary change. Overall, the results from phylogeographic studies indicate that the three-spined stickleback has colonized the Atlantic basin relatively recently but constitutes three major evolutionary lineages in Europe. In addition, the colonization of freshwater appears to result from multiple and independent invasions by the marine conspecifics. Molecular data together with morphology suggest that the most divergent freshwater populations are located in the Balkan Peninsula and these populations deserve a special conservation genetic status without warranting further taxonomical classification. In order to investigate the adaptive divergence in Fennoscandian three-spined stickleback populations several approaches were used. First, sequence variability in the Eda-gene, coding for the number of lateral plates, was concordant with the previously observed global pattern. Full plated allele is in high frequencies among marine populations whereas low plated allele dominates in the freshwater populations. Second, a microsatellite based genome scan identified both indications of balancing and directional selection in the three-spined stickleback genome, i.e. loci with unusually similar or unusually different allele frequencies over populations. The directionally selected loci were mainly associated with the adaptation to freshwater. A follow up study conducting a more detailed analysis in a chromosome region containing a putatively selected gene locus identified a fairly large genomic region affected by natural selection. However, this region contained several gene predictions, all of which might be the actual target of natural selection. All in all, the phylogeographic and adaptive divergence studies indicate that most of the genetic divergence has occurred in the freshwater populations whereas the marine populations have remained relatively uniform.

Sweet potato chlorotic stunt virus: Studies on viral synergism and suppression of RNA silencing

The studies presented in this thesis aimed to a better understanding of the molecular biology of Sweet potato chlorotic stunt virus (SPCSV, Crinivirus, Closteroviridae) and its role in the development of synergistic viral diseases. The emphasis was on the severe sweet potato virus disease (SPVD) that results from a synergistic interaction of SPCSV and Sweet potato feathery mottle virus (SPFMV, Potyvirus, Potyviridae). SPVD is the most important disease affecting sweetpotato. It is manifested as a significant increase in symptom severity and SPFMV titres. This is accompanied by a dramatic sweetpotato yield reduction. SPCSV titres remain little affected in the diseased plants. Viral synergistic interactions have been associated with the suppression of an adaptive general defence mechanism discovered in plants and known as RNA silencing. In the studies of this thesis two novel proteins (RNase3 and p22) identified in the genome of a Ugandan SPCSV isolate were shown to be involved in suppression of RNA silencing. RNase3 displayed a dsRNA-specific endonuclease activity that enhanced the RNA-silencing suppression activity of p22. Comparative analyses of criniviral genomes revealed variability in the gene content at the 3´end of the genomic RNA1. Molecular analyses of different isolates of SPCSV indicated a marked intraspecific heterogeneity in this region where the p22 and RNase3 genes are located. Isolates of the East African strain of SPCSV from Tanzania and Peru and an isolate from Israel were missing a 767-nt fragment that included the p22 gene. However, regardless of the absence of p22, all SPCSV isolates acted synergistically with SPFMV in co-infected sweetpotato, enhanced SPFMV titres and caused SPVD. These results showed that p22 is dispensable for development of SPVD. The role of RNase3 in SPVD was then studied by generating transgenic plants expressing the RNase3 protein. These plants had increased titres of SPFMV (ca. 600-fold higher in comparison with nontransgenic plants) 2-3 weeks after graft inoculation and displayed the characteristic SPVD symptoms. RNA silencing suppression (RSS) activity of RNase3 was detected in agroinfiltrated leaves of Nicotiana bethamiana. In vitro studies showed that RNase3 was able to cleave small interferring RNAs (siRNA) to products of ~14-nt. The data thus identified RNase3 as a suppressor of RNA silencing able to cleave siRNAs. RNase3 expression alone was sufficient for breaking down resistance to SPFMV in sweetpotato and for the development of SPVD. Similar RNase III-like genes exist in animal viruses which points out a novel and possibly more general mechanism of RSS by viruses. A reproducible method of sweetpotato transformation was used to target RNA silencing against the SPCSV polymerase region (RdRp) with an intron-spliced hairpin construct. Hence, engineered resistance to SPCSV was obtained. Ten out of 20 transgenic events challenged with SPCSV alone showed significantly reduced virus titres. This was however not sufficient to prevent SPVD upon coinfection with SPFMV. Immunity to SPCSV seems to be required to control SPVD and targeting of different SPCSV regions need to be assessed in further studies. Based on the identified key role of RNase3 in SPVD the possibility to design constructs that target this gene might prove more efficient in future studies.

Spatial and temporal variability in midge (Nematocera) assemblages in shallow Finnish lakes (60−70 °N) : community-based modelling of past environmental change

Multi- and intralake datasets of fossil midge assemblages in surface sediments of small shallow lakes in Finland were studied to determine the most important environmental factors explaining trends in midge distribution and abundance. The aim was to develop palaeoenvironmental calibration models for the most important environmental variables for the purpose of reconstructing past environmental conditions. The developed models were applied to three high-resolution fossil midge stratigraphies from southern and eastern Finland to interpret environmental variability over the past 2000 years, with special focus on the Medieval Climate Anomaly (MCA), the Little Ice Age (LIA) and recent anthropogenic changes. The midge-based results were compared with physical properties of the sediment, historical evidence and environmental reconstructions based on diatoms (Bacillariophyta), cladocerans (Crustacea: Cladocera) and tree rings. The results showed that the most important environmental factor controlling midge distribution and abundance along a latitudinal gradient in Finland was the mean July air temperature (TJul). However, when the dataset was environmentally screened to include only pristine lakes, water depth at the sampling site became more important. Furthermore, when the dataset was geographically scaled to southern Finland, hypolimnetic oxygen conditions became the dominant environmental factor. The results from an intralake dataset from eastern Finland showed that the most important environmental factors controlling midge distribution within a lake basin were river contribution, water depth and submerged vegetation patterns. In addition, the results of the intralake dataset showed that the fossil midge assemblages represent fauna that lived in close proximity to the sampling sites, thus enabling the exploration of within-lake gradients in midge assemblages. Importantly, this within-lake heterogeneity in midge assemblages may have effects on midge-based temperature estimations, because samples taken from the deepest point of a lake basin may infer considerably colder temperatures than expected, as shown by the present test results. Therefore, it is suggested here that the samples in fossil midge studies involving shallow boreal lakes should be taken from the sublittoral, where the assemblages are most representative of the whole lake fauna. Transfer functions between midge assemblages and the environmental forcing factors that were significantly related with the assemblages, including mean air TJul, water depth, hypolimnetic oxygen, stream flow and distance to littoral vegetation, were developed using weighted averaging (WA) and weighted averaging-partial least squares (WA-PLS) techniques, which outperformed all the other tested numerical approaches. Application of the models in downcore studies showed mostly consistent trends. Based on the present results, which agreed with previous studies and historical evidence, the Medieval Climate Anomaly between ca. 800 and 1300 AD in eastern Finland was characterized by warm temperature conditions and dry summers, but probably humid winters. The Little Ice Age (LIA) prevailed in southern Finland from ca. 1550 to 1850 AD, with the coldest conditions occurring at ca. 1700 AD, whereas in eastern Finland the cold conditions prevailed over a longer time period, from ca. 1300 until 1900 AD. The recent climatic warming was clearly represented in all of the temperature reconstructions. In the terms of long-term climatology, the present results provide support for the concept that the North Atlantic Oscillation (NAO) index has a positive correlation with winter precipitation and annual temperature and a negative correlation with summer precipitation in eastern Finland. In general, the results indicate a relatively warm climate with dry summers but snowy winters during the MCA and a cool climate with rainy summers and dry winters during the LIA. The results of the present reconstructions and the forthcoming applications of the models can be used in assessments of long-term environmental dynamics to refine the understanding of past environmental reference conditions and natural variability required by environmental scientists, ecologists and policy makers to make decisions concerning the presently occurring global, regional and local changes. The developed midge-based models for temperature, hypolimnetic oxygen, water depth, littoral vegetation shift and stream flow, presented in this thesis, are open for scientific use on request.

Photobiological studies of Baltic Sea phytoplankton

Phytoplankton ecology and productivity is one of the main branches of contemporary oceanographic research. Research groups in this branch have increasingly started to utilise bio-optical applications. My main research objective was to critically investigate the advantages and deficiencies of the fast repetition rate (FRR) fluorometry for studies of productivity of phytoplankton, and the responses of phytoplankton towards varying environmental stress. Second, I aimed to clarify the applicability of the FRR system to the optical environment of the Baltic Sea. The FRR system offers a highly dynamic tool for studies of phytoplankton photophysiology and productivity both in the field and in a controlled environment. The FRR metrics obtain high-frequency in situ determinations of the light-acclimative and photosynthetic parameters of intact phytoplankton communities. The measurement protocol is relatively easy to use without phases requiring analytical determinations. The most notable application of the FRR system lies in its potential for making primary productivity (PP) estimations. However, the realisation of this scheme is not straightforward. The FRR-PP, based on the photosynthetic electron flow (PEF) rate, are linearly related to the photosynthetic gas exchange (fixation of 14C) PP only in environments where the photosynthesis is light-limited. If the light limitation is not present, as is usually the case in the near-surface layers of the water column, the two PP approaches will deviate. The prompt response of the PEF rate to the short-term variability in the natural light field makes the field comparisons between the PEF-PP and the 14C-PP difficult to interpret, because this variability is averaged out in the 14C-incubations. Furthermore, the FRR based PP models are tuned to closely follow the vertical pattern of the underwater irradiance. Due to the photoacclimational plasticity of phytoplankton, this easily leads to overestimates of water column PP, if precautionary measures are not taken. Natural phytoplankton is subject to broad-waveband light. Active non-spectral bio-optical instruments, like the FRR fluorometer, emit light in a relatively narrow waveband, which by its nature does not represent the in situ light field. Thus, the spectrally-dependent parameters provided by the FRR system need to be spectrally scaled to the natural light field of the Baltic Sea. In general, the requirement of spectral scaling in the water bodies under terrestrial impact concerns all light-adaptive parameters provided by any active non-spectral bio-optical technique. The FRR system can be adopted to studies of all phytoplankton that possess efficient light harvesting in the waveband matching the bluish FRR excitation. Although these taxa cover the large bulk of all the phytoplankton taxa, one exception with a pronounced ecological significance is found in the Baltic Sea. The FRR system cannot be used to monitor the photophysiology of the cyanobacterial taxa harvesting light in the yellow-red waveband. These taxa include the ecologically-significant bloom-forming cyanobacterial taxa in the Baltic Sea.

Exhaled nitric oxide : Variability and association with bronchial hyperresponsiveness and atopy

Airway inflammation is a key feature of bronchial asthma. In asthma management, according to international guidelines, the gold standard is anti-inflammatory treatment. Currently, only conventional procedures (i.e., symptoms, use of rescue medication, PEF-variability, and lung function tests) were used to both diagnose and evaluate the results of treatment with anti-inflammatory drugs. New methods for evaluation of degree of airway inflammation are required. Nitric oxide (NO) is a gas which is produced in the airways of healthy subjects and especially produced in asthmatic airways. Measurement of NO from the airways is possible, and NO can be measured from exhaled air. Fractional exhaled NO (FENO) is increased in asthma, and the highest concentrations are measured in asthmatic patients not treated with inhaled corticosteroids (ICS). Steroid-treated patients with asthma had levels of FENO similar to those of healthy controls. Atopic asthmatics had higher levels of FENO than did nonatopic asthmatics, indicating that level of atopy affected FENO level. Associations between FENO and bronchial hyperresponsiveness (BHR) occur in asthma. The present study demonstrated that measurement of FENO had good reproducibility, and the FENO variability was reasonable both short- and long-term in both healthy subjects and patients with respiratory symptoms or asthma. We demonstrated the upper normal limit for healthy subjects, which was 12 ppb calculated from two different healthy study populations. We showed that patients with respiratory symptoms who did not fulfil the diagnostic criteria of asthma had FENO values significantly higher than in healthy subjects, but significantly lower than in asthma patients. These findings suggest that BHR to histamine is a sensitive indicator of the effect of ICS and a valuable tool for adjustment of corticosteroid treatment in mild asthma. The findings further suggest that intermittent treatment periods of a few weeks’ duration are insufficient to provide long-term control of BHR in patients with mild persistent asthma. Moreover, during the treatment with ICS changes in BHR and changes in FENO were associated. FENO level was associated with BHR measured by a direct (histamine challenge) or indirect method (exercise challenge) in steroid-naïve symptomatic, non-smoking asthmatics. Although these associations could be found only in atopics, FENO level in nonatopic asthma was also increased. It can thus be concluded that assessment of airway inflammation by measuring FENO can be useful for clinical purposes. The methodology of FENO measurements is now validated. Especially in those patients with respiratory symptoms who did not fulfil the diagnostic criteria of asthma, FENO measurement can aid in treatment decisions. Serial measurement of FENO during treatment with ICS can be a complementary or an alternative method for evaluation in patients with asthma.

Adaptive tree multigrids and simplified spherical harmonics approximation in deterministic neutral and charged particle transport

A new deterministic three-dimensional neutral and charged particle transport code, MultiTrans, has been developed. In the novel approach, the adaptive tree multigrid technique is used in conjunction with simplified spherical harmonics approximation of the Boltzmann transport equation. The development of the new radiation transport code started in the framework of the Finnish boron neutron capture therapy (BNCT) project. Since the application of the MultiTrans code to BNCT dose planning problems, the testing and development of the MultiTrans code has continued in conventional radiotherapy and reactor physics applications. In this thesis, an overview of different numerical radiation transport methods is first given. Special features of the simplified spherical harmonics method and the adaptive tree multigrid technique are then reviewed. The usefulness of the new MultiTrans code has been indicated by verifying and validating the code performance for different types of neutral and charged particle transport problems, reported in separate publications.

Evaluation and modelling of the spatial and temporal variability of particulate matter in urban areas

This thesis contains three subject areas concerning particulate matter in urban area air quality: 1) Analysis of the measured concentrations of particulate matter mass concentrations in the Helsinki Metropolitan Area (HMA) in different locations in relation to traffic sources, and at different times of year and day. 2) The evolution of traffic exhaust originated particulate matter number concentrations and sizes in local street scale are studied by a combination of a dispersion model and an aerosol process model. 3) Some situations of high particulate matter concentrations are analysed with regard to their meteorological origins, especially temperature inversion situations, in the HMA and three other European cities. The prediction of the occurrence of meteorological conditions conducive to elevated particulate matter concentrations in the studied cities is examined. The performance of current numerical weather forecasting models in the case of air pollution episode situations is considered. The study of the ambient measurements revealed clear diurnal variation of the PM10 concentrations in the HMA measurement sites, irrespective of the year and the season of the year. The diurnal variation of local vehicular traffic flows seemed to have no substantial correlation with the PM2.5 concentrations, indicating that the PM10 concentrations were originated mainly from local vehicular traffic (direct emissions and suspension), while the PM2.5 concentrations were mostly of regionally and long-range transported origin. The modelling study of traffic exhaust dispersion and transformation showed that the number concentrations of particles originating from street traffic exhaust undergo a substantial change during the first tens of seconds after being emitted from the vehicle tailpipe. The dilution process was shown to dominate total number concentrations. Minimal effect of both condensation and coagulation was seen in the Aitken mode number concentrations. The included air pollution episodes were chosen on the basis of occurrence in either winter or spring, and having at least partly local origin. In the HMA, air pollution episodes were shown to be linked to predominantly stable atmospheric conditions with high atmospheric pressure and low wind speeds in conjunction with relatively low ambient temperatures. For the other European cities studied, the best meteorological predictors for the elevated concentrations of PM10 were shown to be temporal (hourly) evolutions of temperature inversions, stable atmospheric stability and in some cases, wind speed. Concerning the weather prediction during particulate matter related air pollution episodes, the use of the studied models were found to overpredict pollutant dispersion, leading to underprediction of pollutant concentration levels.

Using IRAP markers for analysis of genetic variability in populations of resource and rare species of plants

Species specific LTR retrotransposons were first cloned in five rare relic species of drug plants located in the Perm’ region. Sequences of LTR retrotransposons were used for PCR analysis based on amplification of repeated sequences from LTR or other sites of retrotransposons (IRAP). Genetic diversity was studied in six populations of rare relic species of plants Adonis vernalis L. by means of the IRAP method; 125 polymorphic IRAP markers were analyzed. Parameters for DNA polymorphism and genetic diversity of A. vernalis populations were determined.

Using IRAP markers for analysis of genetic variability in populations of resource and rare species of plants

Species specific LTR retrotransposons were first cloned in five rare relic species of drug plants located in the Perm’ region. Sequences of LTR retrotransposons were used for PCR analysis based on amplification of repeated sequences from LTR or other sites of retrotransposons (IRAP). Genetic diversity was studied in six populations of rare relic species of plants Adonis vernalis L. by means of the IRAP method; 125 polymorphic IRAP markers were analyzed. Parameters for DNA polymorphism and genetic diversity of A. vernalis populations were determined.

Models of Dispersal and Diversification

Ecology and evolutionary biology is the study of life on this planet. One of the many methods applied to answering the great diversity of questions regarding the lives and characteristics of individual organisms, is the utilization of mathematical models. Such models are used in a wide variety of ways. Some help us to reason, functioning as aids to, or substitutes for, our own fallible logic, thus making argumentation and thinking clearer. Models which help our reasoning can lead to conceptual clarification; by expressing ideas in algebraic terms, the relationship between different concepts become clearer. Other mathematical models are used to better understand yet more complicated models, or to develop mathematical tools for their analysis. Though helping us to reason and being used as tools in the craftmanship of science, many models do not tell us much about the real biological phenomena we are, at least initially, interested in. The main reason for this is that any mathematical model is a simplification of the real world, reducing the complexity and variety of interactions and idiosynchracies of individual organisms. What such models can tell us, however, both is and has been very valuable throughout the history of ecology and evolution. Minimally, a model simplifying the complex world can tell us that in principle, the patterns produced in a model could also be produced in the real world. We can never know how different a simplified mathematical representation is from the real world, but the similarity models do strive for, gives us confidence that their results could apply. This thesis deals with a variety of different models, used for different purposes. One model deals with how one can measure and analyse invasions; the expanding phase of invasive species. Earlier analyses claims to have shown that such invasions can be a regulated phenomena, that higher invasion speeds at a given point in time will lead to a reduction in speed. Two simple mathematical models show that analysis on this particular measure of invasion speed need not be evidence of regulation. In the context of dispersal evolution, two models acting as proof-of-principle are presented. Parent-offspring conflict emerges when there are different evolutionary optima for adaptive behavior for parents and offspring. We show that the evolution of dispersal distances can entail such a conflict, and that under parental control of dispersal (as, for example, in higher plants) wider dispersal kernels are optimal. We also show that dispersal homeostasis can be optimal; in a setting where dispersal decisions (to leave or stay in a natal patch) are made, strategies that divide their seeds or eggs into fractions that disperse or not, as opposed to randomized for each seed, can prevail. We also present a model of the evolution of bet-hedging strategies; evolutionary adaptations that occur despite their fitness, on average, being lower than a competing strategy. Such strategies can win in the long run because they have a reduced variance in fitness coupled with a reduction in mean fitness, and fitness is of a multiplicative nature across generations, and therefore sensitive to variability. This model is used for conceptual clarification; by developing a population genetical model with uncertain fitness and expressing genotypic variance in fitness as a product between individual level variance and correlations between individuals of a genotype. We arrive at expressions that intuitively reflect two of the main categorizations of bet-hedging strategies; conservative vs diversifying and within- vs between-generation bet hedging. In addition, this model shows that these divisions in fact are false dichotomies.

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.