Graduate Theses & Dissertations

Pages

Investigating the sources and fate of monomethylmercury and dimethylmercury in the Arctic marine boundary layer and waters
Monomethylmercury (MMHg), the most bioavailable form of mercury (Hg) and a potent neurotoxin, is present at elevated concentrations in Arctic marine mammals posing serious health threats to the local populations relying on marine food for their subsistence living. The sources of MMHg in the Arctic Ocean surface water and the role of dimethylmercury (DMHg) as a source of MMHg remain unclear. The objective of this research was to determine the sources and fate of methylated Hg species (MMHg and DMHg) in the marine ecosystem by investigating processes controlling the presence of methylated Hg species in the Arctic Ocean marine boundary layer (MBL) and surface waters. A method based on solid phase adsorption on Bond Elut ENV was developed and successfully used for unprecedented measurement of methylated Hg species in the MBL in Hudson Bay (HB) and the Canadian Arctic Archipelago (CAA). MMHg and DMHg concentrations averaged 2.9 ± 3.6 (mean ± SD) and 3.8 ± 3.1 pg m-3, respectively, and varied significantly among sampling sites. MMHg in the MBL is suspected to be the product of marine DMHg degradation in the atmosphere. MMHg summer (June to September) atmospheric wet deposition rates were estimated to be 188 ± 117.5 ng m-2 and 37 ± 21.7 ng m-2 for HB and CAA, respectively, sustaining MMHg concentrations available for bio-magnification in the pelagic food web. The production and loss of methylated Hg species in surface waters was assessed using enriched stable isotope tracers. MMHg production in surface water was observed from methylation of inorganic Hg (Hg(II)) and, for the first time, from DMHg demethylation with experimentally derived rate constants of 0.92 ± 0.82 x 10-3 d-1 and 0.04 ± 0.02 d-1 respectively. DMHg demethyation rate constant (0.98 ± 0.51 d-1) was higher than that of MMHg (0.35 ± 0.25 d-1). Furthermore, relationships with environmental parameters suggest that methylated Hg species transformations in surface water are mainly biologically driven. We propose that in addition to Hg(II) methylation, the main processes controlling MMHg production in the Arctic Ocean surface waters are DMHg demethylation and deposition of atmospheric MMHg. These results are valuable for a better understanding of the cycle of methylated Hg in the Arctic marine environment. Author Keywords: Arctic Ocean, Atmosphere, Demethylation, Dimethylmercury, Methylation, Monomethylmercury
Effect Assessment of Binary Metal Mixtures of Ni, Cu, Zn, and Cd to Daphnia magna
Mixtures of metals occur in surface waters, toxicity of which has drawn world-wide attention due to their crucial role in both ecotoxicology and regulations. The present research was undertaken to study the acute toxicity of binary mixtures of Ni, Cu, Zn, and Cd to the freshwater organism, Daphnia magna. The experimental approach included single and binary metal toxicity tests based on the 48h acute toxicity bioassay of Environment Canada. The acute toxicity of single metals followed the order of Cd > Cu > Zn > Ni. Based on the calculated 48h EC50 value of single metals, a toxic unit (TU) approach was used to combine two metals in a binary mixture, in which 1TU was equal to the 48h EC50 value of a metal in single exposure. The toxicity of binary metal mixtures to D. magna followed the order of Cu-Cd > Cu-Zn > Zn-Cd > Cu-Ni > Zn-Ni > Cd-Ni, which demonstrated three types of toxicity (i.e., less than additive, additive, and greater additive). Predictions from additivity models (including concentration addition (CA) and independent action (IA) models), a generalized linear model (GLM), and a biotic-ligand-like model (BLM-like) were compared to the bioassay results. The CA and the RA models also predicted three types of toxicity of the binary metal mixtures (i.e., less than additive, additive, and greater than additive). However, the CA model mostly overestimated the toxicity of binary mixtures. Predictions from the GLM supported the inclusion of the interaction between two metals in a mixture to predict the toxicity of binary metal mixtures. The binary metal toxicity was also predicted using a BLM-like model based on the calculated concentrations of free ionic forms of the metals, affinity constants, and toxic potency of each metal. In this model, it was hypothesized that the toxicity of metal mixture is the result of competition of metals with Ca2+ at biotic ligands, which can lead to whole-body deficiency of Ca2+ in D. magna. The BLM-like model provided the toxic potency of single metals with the following order, Cu > Cd > Zn > Ni. Although the prediction of the BLM-like model was not in good agreement with the observed toxicity of binary metal mixtures, an overestimation of risk of mixture toxicity was obtained using this model, which could be promising for use in environmental risk assessment. Author Keywords: biotic ligand model, concentration addition, Daphnia magna, independent action, metal toxicity, modeling
collaborative ecotoxicological risk assessment of in-place pollutants in Owen Sound Bay, Lake Huron within the Saugeen Ojibway Nation Territory
Owen Sound Bay, which is located within the traditional territory of the Saugeen Ojibway Nation (SON), is contaminated as a result of historical industrial and shipping activity. Gross contamination of the sediments in the inner part of the Bay (i.e., Owen Sound Harbour) includes high concentrations of polycyclic aromatic hydrocarbons (PAHs) and other organic compounds, as well as metals that may pose a risk to the SON fishery for lake whitefish (Coregonus clupeaformis). However, evaluating the environmental risks posed by contaminated sediments is a challenge, as these risks are dependent upon several factors and require multiple lines of evidence. Including Indigenous communities in environmental risk assessment and the management of those risks is vital for sustaining ecosystem integrity, as well as respecting Treaty Rights. In this study, a risk assessment framework was developed that included several risk assessment tools used in Western science and also encompassed the concerns and values of the SON, including the application of SON-ecological knowledge. Methods for risk evaluation included gathering lines of evidence though community workshops, as well as field sampling in the Bay to determine the concentrations of PAHs and other organic contaminants in sediments and in the water column. Laboratory studies of toxicity to early life stages of lake whitefish and Japanese medaka (Oryzias latipes) and sediment disturbance simulations to evaluate biological responses in juvenile lake whitefish were also completed as lines of evidence. The results indicate that leaving the harbour “as is” without a thorough analysis of remediation options fails to address the concerns of the people within the SON communities. Overall, this research demonstrated a successful process for developing a collaborative risk assessment framework that recognizes the sovereignty of Indigenous peoples and promotes Nation-to-Nation decision making. Author Keywords: biomarkers, Coregonus clupeaformis, Indigenous knowledge, polycyclic aromatic hydrocarbons, risk assessment, source tracking
Hydroclimatic and spatial controls on stream nutrient export from forested catchments
Winter nutrient export from forested catchments is extremely variable from year-to-year and across the landscape of south-central Ontario. Understanding the controls on this variability is critical, as what happens during the winter sets up the timing and nature of the spring snowmelt, the major period of export for water and nutrients from seasonally snow-covered forests. Furthermore, winter processes are especially vulnerable to changes in climate, particularly to shifts in precipitation from snow to rain as air temperatures rise. The objective of this thesis was to assess climatic and topographic controls on variability in stream nutrient export from a series of forested catchments in south-central Ontario. The impacts of climate on the timing and magnitude of winter stream nutrient export, with particular focus on the impact of winter rain-on-snow (ROS) events was investigated through a) analysis of long-term hydrological, chemical and meteorological records and b) high frequency chemical and isotopic measurements of stream and snow samples over two winters. The relationship between topography and variability in stream chemistry among catchments was investigated through a) a series of field and laboratory incubations to measure rates and discern controls on nitrogen mineralization and nitrification and b) analysis of high resolution spatial data to assess relationships between topographic metrics and seasonal stream chemistry. Warmer winters with more ROS events were shown to shift the bulk of nitrate (NO3-N) export earlier in the winter at the expense of spring export; this pattern was not observed in other nutrients [i.e. dissolved organic carbon (DOC), total phosphorus (TP), sulphate (SO4), calcium (Ca)]. Hydrograph separation revealed the majority of ROS flow came from baseflow, but the NO3-N concentrations in rainfall and melting snow were so high that the majority of NO3-N export was due to these two sources. Other nutrient concentrations did not show such a great separation between sources, and thus event export of these nutrients was not as great. Proportionally, catchments with varying topography responded similarly to ROS events, but the absolute magnitude of export varied substantially, due to differences in baseflow NO3-N concentrations. Field and laboratory incubations revealed differences in rates of net NO3-N production between wetland soils and upland soils, suggesting that topographic differences amongst catchments may be responsible for differences in baseflow NO3-N. Spatial analysis of digital elevation models revealed strong relationships between wetland coverage and DOC and dissolved organic nitrogen (DON) concentrations in all seasons, but relationships between topography and NO3-N were often improved by considering only the area within 50 or 100m of the stream channel. This suggests nutrient cycling processes occurring near the stream channel may exert a stronger control over NO3-N stream outflow chemistry. Overall, topography and climate exert strong controls over spatial and temporal variability in stream chemistry at forested catchments; it is important to consider the interaction of these two factors when predicting the effects of future changes in climate or deposition. Author Keywords: biogeochemistry, forest, nitrate, south-central Ontario, stream chemistry, winter
effects of Dissolved Organic Matter (DOM) sources on Pb2+, Zn2+ and Cd2+ binding
Metal binding to dissolved organic matter (DOM) determines metal speciation and strongly influences potential toxicity. The understanding of this process, however, is challenged by DOM source variation, which is not always considered by most existing metal speciation models. Source determines the molecular structure of DOM, including metal binding functional groups. This study has experimentally showed that the allochthonous-dominant DOM (i.e. more aromatic and humic) consistently has higher level of Pb binding than the autochthonous-dominant DOM (i.e. more aliphatic and proteinaceous) by more than two orders of magnitude. This source-discrimination, however, is less noticeable for Zn and Cd, although variation still exceeds a factor of four for both metals. The results indicate that metal binding is source-dependent, but the dependency is metal-specific. Accordingly, metal speciation models, such as the Windermere Humic Aqueous Model (WHAM), needs to consider DOM source variations. The WHAM input of active fraction of DOM participating in metal binding (f) is sensitive to DOM source. The commonly-used f = 0.65 substantially overestimated the Pb and Zn binding to autochthonous-dominant DOM, indicating f needs to be adjusted specifically. The optimal f value (fopt) linearly correlates with optical indexes, showing a potential to estimate fopt using simple absorbance and/or fluorescence measurements. Other DOM properties not optically-characterized may be also important to determine fopt, such as thiol, which shows strong affinity to most toxic metals and whose concentrations are appreciably high in natural waters (< 0.1 to 400 nmol L-1). Other analytical techniques rather than Cathodic Stripping Voltammetry (CSV) are required to accurately quantify thiol concentration for DOM with concentration > 1 mg L-1. To better explain the DOM-source effects, the conditional affinity spectrum (CAS) was calculated using a Fully Optimized ContinUous Spectrum (FOCUS) method. This method not only provides satisfactory goodness-of-fit, but also unique CAS solution. The allochthonous-dominant DOM consistently shows higher Pb affinity than autochthonous-dominant DOM. This source-discrimination is not clearly observed for Zn and Cd. Neither the variability of affinity nor capacity can be fully explained by the variability of individual DOM properties, indicating multiple properties may involve simultaneously. Together, the results help improve WHAM prediction of metal speciation, and consequently, benefit geochemical modelling of metal speciation, such as Biotic Ligand Model for predicting metal toxicity. Author Keywords: Dissolved organic matter, Metal binding, Source, Windermere Humic Aqueous Model
Characterization of a Zn(II)2Cys6 transcription factor in Ustilago maydis and its role in pathogenesis
Ustilago maydis (D.C.) Corda is a biotrophic pathogen that secretes effectors to establish and maintain a relationship with its host, Zea mays. In this pathosystem, the molecular function of effectors is well-studied, but the regulation of effector gene expression remains largely unknown. This study characterized Zfp1, a putative U. maydis Zn(II)2Cys6 transcription factor, as a modulator of effector gene expression. The amino acid sequence of Zfp1 indicated the presence of a GAL4-like zinc binuclear cluster as well as a fungal specific transcription factor domain. Nuclear localization was confirmed by tagging Zfp1 with enhanced green fluorescent protein. Deletion of zfp1 resulted in attenuated hyphal growth, reduced infection frequency, an arrest in pathogenic development, and decreased anthocyanin production. This phenotype can be attributed to the altered transcript levels of genes encoding predicted and confirmed U. maydis effectors in the zfp1 deletion strain during pathogenic growth. Complementation of zfp1 deletion strain with tin2, an effector involved in anthocyanin induction, suggested this effector is downstream of Zfp1 and its expression is influenced by this transcription factor during in planta growth. When wild-type zfp1 was ectopically inserted in the zfp1 deletion strain, pathogenesis and virulence were partially restored. This, coupled with zfp1 over-expression strains having a similar phenotype as the deletion strains, suggested Zfp1 may interact with other proteins for full function. These findings show that Zfp1, in conjunction with one or more binding partners, contributes to U. maydis pathogenesis, virulence, and anthocyanin production through the regulation of effector gene expression. Author Keywords: effector, pathogenesis, transcription factor, Ustilago maydis, Zea mays, zinc finger
Comparative phylogeography in conservation biology
Phylogeographic histories of taxa around the Great Lakes region in North America are relevant to a range of ongoing issues including conservation management and biological invasions. In this thesis I investigated the comparative phylogeographic histories of plant species with disjunct distributions and plant species with continuous distributions around the Great Lakes region; this is a very dynamic geographic area with relatively recent colonisation histories that have been influenced by a range of factors including postglacial landscape modifications, and more recently, human-mediated dispersion. I first characterized four species that have disjunct populations in the Great Lakes region: (Bartonia paniculata subsp. paniculata, Empetrum nigrum, Sporobolus heterolepis, and Carex richardsonii). Through comparisons of core and disjunct populations, I found that a range of historical processes have resulted in two broad scenarios: in the first scenario, genetically distinct disjunct and core populations diverged prior to the last glacial cycle, and in the second scenario more recent vicariant events have resulted in genetically similar core and disjunct populations. The former scenario has important implications for conservation management. I then characterized the Typha species complex (T. latifolia, T. angustifolia, T. x glauca), which collectively represent species with continuous distributions. Recent microevolutionary processes, including hybridization, introgression, and intercontinental dispersal, obscure the phylogeographic patterns and complicate the evolutionary history of Typha spp. around the Great Lakes region, and have resulted in the growing dominance of non-native lineages. A broader geographical comparison of Typha spp. lineages from around the world identified repeated cryptic dispersal and long-distant movement as important phylogeographic influences. This research has demonstrated that comparisons of regional and global evolutionary histories can provide insight into historical and contemporary processes useful for management decisions in conservation biology and invasive species. Author Keywords: chloroplast DNA, conservation genetics, disjunct populations, invasive species, phylogeography, postglacial recolonisation
Size and fluorescence properties of allochthonous dissolved organic matter
Dissolved organic matter (DOM) is a mixture of molecules with dynamic structure and composition that are ubiquitous in aquatic systems. DOM has several important functions in both natural and engineered systems, such as supporting microorganisms, governing the toxicity of metals and other pollutants, and controlling the fate of dissolved carbon. The structure and composition of DOM determine its reactivity, and hence its effectiveness in these ecosystem functions. While the structure, composition, and reactivity of riverine and marine DOM have been previously investigated, those of allochthonous DOM collected prior to exposure to microbes and sunlight have received scant attention. The following dissertation constitutes the first in-depth study of the structure, composition, and reactivity of allochthonous DOM at its point of origin (i.e. leaf leachates, LLDOM), as detected by measuring its size and optical properties. Concomitantly, novel chemometric methods were developed to interpret size-resolved data obtained using asymmetrical flow field-flow fractionation, including spectral deconvolution and the application of machine learning algorithms such as self-organizing maps to fluorescence data using a dataset of more than 1000 fluorescence excitation-emission matrices. The size and fluorescence properties of LLDOM are highly distinct. Indeed, LLDOM was correctly classified as one of 13 species/sources with 92.5% accuracy based on its fluorescence composition, and LLDOM was distinguished from riverine DOM sampled from eight different rivers with 98.3% accuracy. Additionally, both fluorescence and size properties were effective conservative tracers of DOC contribution in pH-controlled mixtures of leaf leachates and riverine DOM over two weeks. However, the structure of LLDOM responded differently to pH changes for leaves/needles from different tree species, and for older needles. Structural changes were non-reversible. Copper-binding strength (log K) differed for the different fluorescent components of DOM in a single allochthonous source by more than an order of magnitude (4.73 compared to 6.11). Biotransformation preferentially removed protein/polyphenol-like fluorescence and altered copper-binding parameters: log K increased from 4.7 to 5.5 for one fluorescent component measured by fluorescence quenching, but decreased from 7.2 to 5.8 for the overall DOM, as measured using voltammetry. The complexing capacity of DOM increased in response to biotransformation for both fluorescent and total DOM. The relationship between fluorescence and size properties was consistent for fresh allochthonous DOM, but differed in aged material. Since the size and fluorescence properties of LLDOM are strikingly different from those of riverine DOM, deeper investigation into transformative pathways and mixing processes is required to elucidate the contribution of riparian plant species to DOM signatures in rivers. Author Keywords: Analytical chemistry, Chemometrics, Dissolved organic matter (DOM), Field-flow fractionation, Fluorescence spectroscopy, Parallel factor analysis (PARAFAC)
Interactome Study of Giardia Intestinalis Cytochromes B5
Giardia intestinalis is an anaerobic protozoan that lacks common eukaryotic heme-dependent respiratory complexes and does not encode any proteins involved in heme biosynthesis. Nevertheless, the parasite encodes several hemeproteins, including three members of the Type II cytochrome b5 sub-group of electron transport proteins found in anaerobic protist and amitochondriate organisms. Unlike the more well-characterized cytochrome b5s of animals, no function has been ascribed to any of the Type II proteins. To explore the functions of these Giardia cytochromes (gCYTB5s), I used bioinformatics, immunofluorescence microscopy (IFM) and co-immunoprecipitation assays. The protein-protein interaction in silico prediction tool, STRING, failed to identify relevant interacting partners for any of the Type II cytochromes b5 from Giardia or other organisms. Differential cellular localization of the gCYTB5s was detected by IFM: gCYTB5-I in the perinuclear space; gCYTB5-II in the cytoplasm with a staining pattern similar to peripheral vacuole-associated protein; and gCYTB5-III in the nucleus. Co-immunoprecipitation with the gCYTB5s as bait identified potential interacting proteins for each isotype. The most promising candidate is the uncharacterized protein GL50803_9861, which was identified in the immunoprecipitate of both gCYTB5-I and II, and which co-localizes with both. Structural analysis of GL50803_9861 using Swiss Model, Phyre2, I-TASSER and RaptorX predicts the presence of a nucleotide-binding domain, which is consistent with a potential redox role involving nicotinamide or flavin-containing cofactors. Finally, the protein GL50803_7204 which contains a RNA/DNA binding domain was identified a potential partner of gCYTB5-III. These findings represent the first steps in the discovery of the roles played by these proteins in Giardia. Author Keywords: Cytochrome b5, Giardia intestinalis, Heme, Interactome, Protein structure prediction
Genome annotation, gene characterization, and the functional analysis of natural antisense transcripts in the fungal plant pathogen Ustilago maydis
Ustilago maydis (DC) Corda is the causal agent of 'common smut of corn'. Completion of the U. maydis lifecycle is dependent on development inside its host, Zea mays. Symptoms of U. maydis infection include chlorosis and the formation of tumours on all aerial corn tissues. Within the tumours, thick-walled diploid teliospores form; these are the reproductive and dispersal agent for the fungus. U. maydis is the model to study basidiomycete biotrophic plant-pathogen interactions. It holds this status in part because of the completely sequenced 20.5 Mb genome; however, thorough genome annotation is required to fully realize the value of this resource. The research presented here improved U. maydis genome annotation through the analysis of cDNA library sequences and comparative genomics. These analyses identified and characterized pathogenesis-related genes, and identified putative meiosis genes. This enabled the use of U. maydis as a model for investigating 'host-induced' meiosis. Further, the cDNA library analyses identified non-coding RNAs (ncRNAs) and natural antisense transcripts (NATs). NATs are endogenous RNA molecules with regions complementary to a protein-coding transcript. Although NATs have been identified in a wide variety of mammals, plants, and fungi, very few have been functionally characterized. Over 200 U. maydis NATs were annotated by analyzing full-length cDNA sequences. NAT structural features were characterized. Strand-specific RT-PCR was used to detect NATs in U. maydis and in a related smut fungus, U. hordei. The data supported a common role for NATs in smut teliospore development, independent of the RNA interference pathway. Analysis of the expression of one U. maydis NAT, as-um02151, in haploid cells, led to a model for NAT function in U. maydis during teliospore dormancy. This model proposed NATs facilitate the maintenance of stored mRNAs through the formation of double-stranded RNA. In testing this model, it was determined that the deletion of two separate upstream regulatory regions, one of which contained a ncRNA (ncRNA1), altered NAT levels and decreased pathogenesis. These studies strengthened U. maydis as a model organism, and began the functional investigation of NATs in U. maydis, which identified a new class of fungal pathogenesis genes. Author Keywords: cDNA library analysis, genome annotation, mRNA stability, natural antisense transcripts, pathogenesis, Ustilago maydis
Characterization of Synthetic and Natural Se8 and Related SenSm Compounds by Gas Chromatography-Mass Spectrometry
Elemental selenium has been extensively quantitatively measured in sediments; however, its physical composition is largely unknown, despite it being the dominant selenium species in some reducing environments. Here, for the first time, it is shown that small, cyclic selenium compounds can account for a quantitatively-relevant fraction of the total elemental selenium present. A new method was developed to analyze for cyclooctaselenium (Se8) in both synthetic samples and selenium-impacted sediments. Despite some analytical limitations, this gas chromatography-mass spectrometry (GC-MS) method is the first GC-MS method developed to identify and quantify Se8 in sediments. Once this method was established, it was then applied to more complex systems: first, the identification of compounds in mixed selenium-sulfur melt solutions, and then the determination of SenSm in selenium-impacted sediments. Despite complications arising from pronounced fragmentation in the ion source, assignment of definitive molecular formulae to chromatographically-resolved peaks was possible for five compounds. Developing a fully quantitative method to obtain elemental ratio information can aid in the assignment of molecular formulae to chromatographically-resolved SeS-containing chromatographic peaks. Coupling the existing gas chromatography method to an inductively coupled plasma-mass spectrometer (ICP-MS) system should accomplish this. However, due to a number of complications, this was not completed successfully during the duration of this thesis project. High detection limits for sulfur, retention time discrepancies, and inconsistent injection results between the GC-MS and GC-ICP-MS system led to difficulties in comparing results between both analytical methods. Despite these limitations, GC-ICP-MS remains the most promising method for the identification and quantification of SenSm compounds in synthetic melt mixtures and selenium impacted sediments. Author Keywords: gas chromatography-mass spectrometry, sediments, selenium
Carbon Exchange along a Natural Gradient of Deciduous Shrub Coverage in the Low-Arctic
Arctic terrestrial ecosystems have experienced substantial structural and compositional changes in response to warming climate in recent decades, especially the expansion of shrub species in Arctic tundra. Climatic and vegetation changes could feedback to the global climate by changing the carbon balance of Arctic tundra. The objective of this thesis was to investigate the influence of increased shrub coverage on carbon exchange processes between atmosphere and the Arctic tundra ecosystem. In this study a space-for-time substitution was used, referred to as a shrub expansion “chronosequence”, with three sites along a natural gradient of deciduous shrub coverage in the Canadian low Arctic. Leaf-level photosynthetic capacity (Amax) of dominating birch shrub Betula glandulosa (Michx.) was significantly higher (P<0.05) at the site where shrubs were more abundant and taller than at the other sites. For all sites, mean Amax in 2014 was significantly lower than in 2013, in part potentially due to differences in precipitation distribution. Bulk soil respiration (RS) rate was significantly higher (P<0.05) at the site with more shrubs compared with the other sites. The differences in RS across sites appeared to be driven by differences in soil physiochemical properties, such as soil nitrogen and soil bulk density rather than soil microclimate factors (e.g. soil temperature, moisture). The three sites were either annual CO2 sources (NEP<0) to the atmosphere or CO2 neutral, with strongest annual CO2 sources (-44.1±7.0 g C m-2) at the site with most shrubs. Overall this study suggests that shrubs tundra carbon balance will change with shrub expansion and that shrub ecosystems in the Arctic currently act as annual carbon sources or neutral to the atmospheric CO2 and further shrub expansion might strengthen the CO2 emissions, causing a positive feedback to the warming climate. Author Keywords: arctic tundra, carbon exchange, climate change, photosynthetic capacity, shrub expansion, soil respiration

Pages

Search Our Digital Collections

Query

Enabled Filters

  • (-) ≠ Nocera
  • (-) = Environmental and Life Sciences
  • (-) ≠ Ecology
  • (-) ≠ Water resources management
  • (-) = Doctor of Philosophy

Filter Results

Date

2004 - 2024
(decades)
Specify date range: Show
Format: 2024/05/22