Ecological drivers of a vector borne pathogen: distribution and abundance of Borrelia burgdorferi sensu lato and its vector Ixodes ricinus in Scotland

Vector-borne disease emergence in recent decades has been associated with different environmental drivers including changes in habitat, hosts and climate. Lyme borreliosis is among the most important vector-borne diseases in the Northern hemisphere and is an emerging disease in Scotland. Transmitted by Ixodid tick vectors between large numbers of wild vertebrate host species, Lyme borreliosis is caused by bacteria from the Borrelia burgdorferi sensu lato species group. Ecological studies can inform how environmental factors such as host abundance and community composition, habitat and landscape heterogeneity contribute to spatial and temporal variation in risk from B. burgdorferi s.l. In this thesis a range of approaches were used to investigate the effects of vertebrate host communities and individual host species as drivers of B. burgdorferi s.l. dynamics and its tick vector Ixodes ricinus. Host species differ in reservoir competence for B. burgdorferi s.l. and as hosts for ticks. Deer are incompetent transmission hosts for B. burgdorferi s.l. but are significant hosts of all life-stages of I. ricinus. Rodents and birds are important transmission hosts of B. burgdorferi s.l. and common hosts of immature life-stages of I. ricinus. In this thesis, surveys of woodland sites revealed variable effects of deer density on B. burgdorferi prevalence, from no effect (Chapter 2) to a possible ‘dilution’ effect resulting in lower prevalence at higher deer densities (Chapter 3). An invasive species in Scotland, the grey squirrel (Sciurus carolinensis), was found to host diverse genotypes of B. burgdorferi s.l. and may act as a spill-over host for strains maintained by native host species (Chapter 4). Habitat fragmentation may alter the dynamics of B. burgdorferi s.l. via effects on the host community and host movements. In this thesis, there was lack of persistence of the rodent associated genospecies of B. burgdorferi s.l. within a naturally fragmented landscape (Chapter 3). Rodent host biology, particularly population cycles and dispersal ability are likely to affect pathogen persistence and recolonization in fragmented habitats. Heterogeneity in disease dynamics can occur spatially and temporally due to differences in the host community, habitat and climatic factors. Higher numbers of I. ricinus nymphs, and a higher probability of detecting a nymph infected with B. burgdorferi s.l., were found in areas with warmer climates estimated by growing degree days (Chapter 2). The ground vegetation type associated with the highest number of I. ricinus nymphs varied between studies in this thesis (Chapter 2 & 3) and does not appear to be a reliable predictor across large areas. B. burgdorferi s.l. prevalence and genospecies composition was highly variable for the same sites sampled in subsequent years (Chapter 2). This suggests that dynamic variables such as reservoir host densities and deer should be measured as well as more static habitat and climatic factors to understand the drivers of B. burgdorferi s.l. infection in ticks. Heterogeneity in parasite loads amongst hosts is a common finding which has implications for disease ecology and management. Using a 17-year data set for tick infestations in a wild bird community in Scotland, different effects of age and sex on tick burdens were found among four species of passerine bird (Chapter 5). There were also different rates of decline in tick burdens among bird species in response to a long term decrease in questing tick pressure over the study. Species specific patterns may be driven by differences in behaviour and immunity and highlight the importance of comparative approaches. Combining whole genome sequencing (WGS) and population genetics approaches offers a novel approach to identify ecological drivers of pathogen populations. An initial analysis of WGS from B. burgdorferi s.s. isolates sampled 16 years apart suggests that there is a signal of measurable evolution (Chapter 6). This suggests demographic analyses may be applied to understand ecological and evolutionary processes of these bacteria. This work shows how host communities, habitat and climatic factors can affect the local transmission dynamics of B. burgdorferi s.l. and the potential risk of infection to humans. Spatial and temporal heterogeneity in pathogen dynamics poses challenges for the prediction of risk. New tools such as WGS of the pathogen (Chapter 6) and blood meal analysis techniques will add power to future studies on the ecology and evolution of B. burgdorferi s.l.

Phytoplankton versus macrophyte contribution to primary production and biogeochemical cycles of a coastal mesotidal system. A modelling approach

This study presents an assessment of the contributions of various primary producers to the global annual production and N/P cycles of a coastal system, namely the Arcachon Bay, by means of a numerical model. This 3D model fully couples hydrodynamic with ecological processes and simulates nitrogen, silicon and phosphorus cycles as well as phytoplankton, macroalgae and seagrasses. Total annual production rates for the different components were calculated for different years (2005, 2007 and 2009) during a time period of drastic reduction in seagrass beds since 2005. The total demand of nitrogen and phosphorus was also calculated and discussed with regards to the riverine inputs. Moreover, this study presents the first estimation of particulate organic carbon export to the adjacent open ocean. The calculated annual net production for the Arcachon Bay (except microphytobenthos, not included in the model) ranges between 22,850 and 35,300 tons of carbon. The main producers are seagrasses in all the years considered with a contribution ranging from 56% to 81% of global production. According to our model, the -30% reduction in seagrass bed surface between 2005 and 2007, led to an approximate 55% reduction in seagrass production, while during the same period of time, macroalgae and phytoplankton enhanced their productions by about +83% and +46% respectively. Nonetheless, the phytoplankton production remains about eightfold higher than the macroalgae production. Our results also highlight the importance of remineralisation inside the Bay, since riverine inputs only fulfill at maximum 73% nitrogen and 13% phosphorus demands during the years 2005, 2007 and 2009. Calculated advection allowed a rough estimate of the organic matter export: about 10% of the total production in the bay was exported, originating mainly from the seagrass compartment, since most of the labile organic matter was remineralised inside the bay.

Probabilistic approach of water residence time and connectivity using Markov Chains with application to tidal embayments

Markov Chain analysis was recently proposed to assess the time scales and preferential pathways into biological or physical networks by computing residence time, first passage time, rates of transfer between nodes and number of passages in a node. We propose to adapt an algorithm already published for simple systems to physical systems described with a high resolution hydrodynamic model. The method is applied to bays and estuaries on the Eastern Coast of Canada for their interest in shellfish aquaculture. Current velocities have been computed by using a 2 dimensional grid of elements and circulation patterns were summarized by averaging Eulerian flows between adjacent elements. Flows and volumes allow computing probabilities of transition between elements and to assess the average time needed by virtual particles to move from one element to another, the rate of transfer between two elements, and the average residence time of each system. We also combined transfer rates and times to assess the main pathways of virtual particles released in farmed areas and the potential influence of farmed areas on other areas. We suggest that Markov chain is complementary to other sets of ecological indicators proposed to analyse the interactions between farmed areas - e.g. depletion index, carrying capacity assessment. Markov Chain has several advantages with respect to the estimation of connectivity between pair of sites. It makes possible to estimate transfer rates and times at once in a very quick and efficient way, without the need to perform long term simulations of particle or tracer concentration.

Applied evolution: An integrated approach to studying life history traits in response to drug selection

The use of chemical control measures to reduce the impact of parasite and pest species has frequently resulted in the development of resistance. Thus, resistance management has become a key concern in human and veterinary medicine, and in agricultural production. Although it is known that factors such as gene flow between susceptible and resistant populations, drug type, application methods, and costs of resistance can affect the rate of resistance evolution, less is known about the impacts of density-dependent eco-evolutionary processes that could be altered by drug-induced mortality. The overall aim of this thesis was to take an experimental evolution approach to assess how life history traits respond to drug selection, using a free-living dioecious worm (Caenorhabditis remanei) as a model. In Chapter 2, I defined the relationship between C. remanei survival and Ivermectin dose over a range of concentrations, in order to control the intensity of selection used in the selection experiment described in Chapter 4. The dose-response data were also used to appraise curve-fitting methods, using Akaike Information Criterion (AIC) model selection to compare a series of nonlinear models. The type of model fitted to the dose response data had a significant effect on the estimates of LD50 and LD99, suggesting that failure to fit an appropriate model could give misleading estimates of resistance status. In addition, simulated data were used to establish that a potential cost of resistance could be predicted by comparing survival at the upper asymptote of dose-response curves for resistant and susceptible populations, even when differences were as low as 4%. This approach to dose-response modeling ensures that the maximum amount of useful information relating to resistance is gathered in one study. In Chapter 3, I asked how simulations could be used to inform important design choices used in selection experiments. Specifically, I focused on the effects of both within- and between-line variation on estimated power, when detecting small, medium and large effect sizes. Using mixed-effect models on simulated data, I demonstrated that commonly used designs with realistic levels of variation could be underpowered for substantial effect sizes. Thus, use of simulation-based power analysis provides an effective way to avoid under or overpowering a study designs incorporating variation due to random effects. In Chapter 4, I 3 investigated how Ivermectin dosage and changes in population density affect the rate of resistance evolution. I exposed replicate lines of C. remanei to two doses of Ivermectin (high and low) to assess relative survival of lines selected in drug-treated environments compared to untreated controls over 10 generations. Additionally, I maintained lines where mortality was imposed randomly to control for differences in density between drug treatments and to distinguish between the evolutionary consequences of drug treatment versus ecological processes affected by changes in density-dependent feedback. Intriguingly, both drug-selected and random-mortality lines showed an increase in survivorship when challenged with Ivermectin; the magnitude of this increase varied with the intensity of selection and life-history stage. The results suggest that interactions between density-dependent processes and life history may mediate evolved changes in susceptibility to control measures, which could result in misleading conclusions about the evolution of heritable resistance following drug treatment. In Chapter 5, I investigated whether the apparent changes in drug susceptibility found in Chapter 4 were related to evolved changes in life-history of C. remanei populations after selection in drug-treated and random-mortality environments. Rapid passage of lines in the drug-free environment had no effect on the measured life-history traits. In the drug-free environment, adult size and fecundity of drug-selected lines increased compared to the controls but drug selection did not affect lifespan. In the treated environment, drug-selected lines showed increased lifespan and fecundity relative to controls. Adult size of randomly culled lines responded in a similar way to drug-selected lines in the drug-free environment, but no change in fecundity or lifespan was observed in either environment. The results suggest that life histories of nematodes can respond to selection as a result of the application of control measures. Failure to take these responses into account when applying control measures could result in adverse outcomes, such as larger and more fecund parasites, as well as over-estimation of the development of genetically controlled resistance. In conclusion, my thesis shows that there may be a complex relationship between drug selection, density-dependent regulatory processes and life history of populations challenged with control measures. This relationship could have implications for how resistance is monitored and managed if life histories of parasitic species show such eco-evolutionary responses to drug application.

Integrating Environmental Justice and Social-Ecological Resilience for Successful Adaptation to Climate Change: Lessons from African American Communities on the Eastern Shore of the Chesapeake Bay

This research concerns the conceptual and empirical relationship between environmental justice and social-ecological resilience as it relates to climate change vulnerability and adaptation. Two primary questions guided this work. First, what is the level of resilience and adaptive capacity for social-ecological systems that are characterized by environmental injustice in the face of climate change? And second, what is the role of an environmental justice approach in developing adaptation policies that will promote social-ecological resilience? These questions were investigated in three African American communities that are particularly vulnerable to flooding from sea-level rise on the Eastern Shore of the Chesapeake Bay. Using qualitative and quantitative methods, I found that in all three communities, religious faith and the church, rootedness in the landscape, and race relations were highly salient to community experience. The degree to which these common aspects of the communities have imparted adaptive capacity has changed over time. Importantly, a given social-ecological factor does not have the same effect on vulnerability in all communities; however, in all communities political isolation decreases adaptive capacity and increases vulnerability. This political isolation is at least partly due to procedural injustice, which occurs for a number of interrelated reasons. This research further revealed that while all stakeholders (policymakers, environmentalists, and African American community members) generally agree that justice needs to be increased on the Eastern Shore, stakeholder groups disagree about what a justice approach to adaptation would look like. When brought together at a workshop, however, these stakeholders were able to identify numerous challenges and opportunities for increasing justice. Resilience was assessed by the presence of four resilience factors: living with uncertainty, nurturing diversity, combining different types of knowledge, and creating opportunities for self-organization. Overall, these communities seem to have low resilience; however, there is potential for resilience to increase. Finally, I argue that the use of resilience theory for environmental justice communities is limited by the great breadth and depth of knowledge required to evaluate the state of the social-ecological system, the complexities of simultaneously promoting resilience at both the regional and local scale, and the lack of attention to issues of justice.

Ecological effects of flow and temperature regimes on fish communities in temperate rivers

Les organismes aquatiques sont adaptés à une grande variabilité hydrique et thermique des rivières. Malgré ceci, la régulation des eaux suscite des changements aux débits qui peuvent provoquer des impacts négatifs sur la biodiversité et les processus écologiques en rivière. Celle-ci peut aussi causer des modifications au niveau des régimes thermiques et des caractéristiques de l’habitat du poisson. Des données environnementales et biologiques décrivant l’habitat du poisson existent, mais elles sont incomplètes pour plusieurs rivières au Canada et de faible qualité, limitant les relations quantitatives débit-température-poissons à un petit nombre de rivières ou à une région étudiée. La recherche menée dans le cadre de mon doctorat concerne les impacts de la génération d'hydroélectricité sur les rivières; soit les changements aux régimes hydriques et thermiques reliés à la régulation des eaux sur la variation des communautés ichtyologiques qui habitent les rivières régulées et naturelles au Canada. Suite à une comparaison d’échantillonnage de pêche, une méthode constante pour obtenir des bons estimés de poisson (richesse, densité et biomasse des espèces) a été établie pour évaluer la structure de la communauté de poissons pour l’ensemble des rivières ciblées par l’étude. Afin de mieux comprendre ces changements environnementaux, les principales composantes décrivant ces régimes ont été identifiées et l’altération des régimes hydriques pour certaines rivières régulées a été quantifiée. Ces résultats ont servi à établir la relation significative entre le degré de changement biotique et le degré de changement hydrique pour illustrer les différences entre les régimes de régulation. Pour faire un complément aux indices biotiques déjà calculés pour l’ensemble des communautés de poissons (diversité, densité et biomasse des espèces par rivière), les différences au niveau des guildes de poissons ont été quantifiées pour expliquer les divers effets écologiques dus aux changements de régimes hydriques et thermiques provenant de la gestion des barrages. Ces derniers résultats servent à prédire pour quels traits écologiques ou groupes d’espèces de poissons les composantes hydriques et thermiques sont importantes. De plus, ces derniers résultats ont servi à mettre en valeur les variables décrivant les régimes thermiques qui ne sont pas toujours inclues dans les études hydro-écologiques. L’ensemble des résultats de cette thèse ont des retombées importantes sur la gestion des rivières en évaluant, de façon cohérente, l’impact de la régulation des rivières sur les communautés de poissons et en développant des outils de prévision pour la restauration des écosystèmes riverains.

First metabolomic approach of the epiphytic bacteria-marine diatom Haslea ostrearia relationships

The marine diatom Haslea ostrearia [1] produces a water-soluble blue-pigment named marennine [2] of economic interest. But the lack of knowledge of the ecological conditions, under which this microalga develops in its natural ecosystem, more especially bacteria H. ostrearia interactions, prevents any optimization of its culture in well-controlled conditions. The structure of the bacterial community was analyzed by PCR-TTGE before and after the isolation of H. ostrearia cells recovered from 4 localities, to distinguish the relative part of the biotope and the biocenose and eventually to describe the temporal dynamic of the structure of the bacterial community at two time-scales. The differences in genetic fingerprints, more especially high between two H. ostrearia isolates (HO-R and HO-BM) showed also the highest differences in the bacterial structure [3] as the result of specific metabolomics profiles. The non-targeted metabolomic investigation showed that these profiles were more distinct in case of bacteria-alga associations than for the H. ostrearia monoculture Here we present a Q-TOF LC/MS metabolomic fingerprinting approach [3]: - to investigate differential metabolites of axenic versus non axenic H. ostrearia cultures. - to focus on the specific metabolites of a bacterial surrounding associated with the activation or inhibition of the microalga growing. The Agilent suite of data processing software makes feature finding, statistical analysis, and identification easier. This enables rapid transformation of complex raw data into biologically relevant metabolite information.

Preliminary metabolomic approach on the bacterial structure community of Haslea ostrearia

The marine diatom Haslea ostrearia [1] produces a water-soluble blue-pigment named marennine [2] of economic interest. But the lack of knowledge of the ecological conditions, under which this microalga develops in its natural ecosystem, more especially bacteria H. ostrearia interactions, prevents any optimization of its culture in well-controlled conditions. The structure of the bacterial community was analyzed by PCR-TTGE before and after the isolation of H. ostrearia cells recovered from 4 localities, to distinguish the relative part of the biotope and the biocenose and eventually to describe the temporal dynamic of the structure of the bacterial community at two time-scales. The differences in genetic fingerprints, more especially high between two H. ostrearia isolates (HO-R and HO-BM) showed also the highest differences in the bacterial structure [3] as the result of specific metabolomics profiles. The non-targeted metabolomic investigation showed that these profiles were more distinct in case of bacteria-alga associations than for the H. ostrearia monoculture Here we present a Q-TOF LC/MS metabolomic fingerprinting approach [3]: - to investigate differential metabolites of axenic versus non axenic H. ostrearia cultures. - to focus on the specific metabolites of a bacterial surrounding associated with the activation or inhibition of the microalga growing. The Agilent suite of data processing software makes feature finding, statistical analysis, and identification easier. This enables rapid transformation of complex raw data into biologically relevant metabolite information.

Resilience to Climate Change: An Ethnographic Approach

Global projections for climate change impacts produce a startling picture of the future for low-lying coastal communities. The United States’ Chesapeake Bay region and especially marginalized and rural communities will be severely impacted by sea level rise and other changes over the next one hundred years. The concept of resilience has been theorized as a measure of social-ecological system health and as a unifying framework under which people can work together towards climate change adaptation. But it has also been critiqued for the way in which it does not adequately take into account local perspective and experiences, bringing into question the value of this concept as a tool for local communities. We must be sure that the concerns, weaknesses, and strengths of particular local communities are part of the climate change adaptation, decision-making, and planning process in which communities participate. An example of this type of planning process is the Deal Island Marsh and Community Project (DIMCP), a grant funded initiative to build resilience within marsh ecosystems and communities of the Deal Island Peninsula area of Maryland (USA) to environmental and social impacts from climate change. I argue it is important to have well-developed understandings of vulnerabilities and resiliencies identified by local residents and others to accomplish this type of work. This dissertation explores vulnerability and resilience to climate change using an engaged and ethnographic anthropological perspective. Utilizing participant observation, semi-structured and structured interviews, text analysis, and cultural domain analysis I produce an in-depth perspective of what vulnerability and resilience means to the DIMCP stakeholder network. Findings highlight significant vulnerabilities and resiliencies inherent in the local area and how these interface with additional vulnerabilities and resiliencies seen from a nonlocal perspective. I conclude that vulnerability and resilience are highly dynamic and context-specific for the local community. Vulnerabilities relate to climate change and other social and environmental changes. Resilience is a long-standing way of life, not a new concept related specifically to climate change. This ethnographic insight into vulnerability and resilience provides a basis for stronger engagement in collaboration and planning for the future.

Harmonised framework for ecological risk assessment of sediments from ports and estuarine zones of North and South Atlantic

This paper presents a harmonised framework of sediment quality assessment and dredging material characterisation for estuaries and port zones of North and South Atlantic. This framework, based on the weight-of-evidence approach, provides a structure and a process for conducting sediment/dredging material assessment that leads to a decision. The main structure consists of step 1 (examination of available data); step 2 (chemical characterisation and toxicity assessment); decision 1 (any chemical level higher than reference values? are sediments toxic?); step 3 (assessment of benthic community structure); step 4 (integration of the results); decision 2 (are sediments toxic or benthic community impaired?); step 5 (construction of the decision matrix) and decision 3 (is there environmental risk?). The sequence of assessments may be interrupted when the information obtained is judged to be sufficient for a correct characterisation of the risk posed by the sediments/dredging material. This framework brought novel features compared to other sediment/dredging material risk assessment frameworks: data integration through multivariate analysis allows the identification of which samples are toxic and/or related to impaired benthic communities; it also discriminates the chemicals responsible for negative biological effects; and the framework dispenses the use of a reference area. We demonstrated the successful application of this framework in different port and estuarine zones of the North (Gulf of Cadiz) and South Atlantic (Santos and Paranagua Estuarine Systems).

Ecological infrastructures in sustainable olive growing: studies about Prays oleae (Bernard) and its natural enemies

Doutoramento em Engenharia Agronómica - Instituto Superior de Agronomia - UL

An integrated approach to delimiting species borders in the genus Chrysotoxum Meigen, 1803 (Diptera: Syrphidae), with description of two new species

Integrative taxonomy tests the validity of taxa using methods additional to traditional morphology. The existence of two different morphotypes in specimens identified as Chrysotoxum vernale Loew (Diptera: Syrphidae) prompted their taxonomic study using an integrative approach that included morphology, wing and male-surstylus geometric morphometrics, genetic and ecological analyses. As a result, a new species is recognised, Chrysotoxum montanum Nedeljković & Vujić sp. nov., and C. vernale is re-defined. A lectotype and paralectotypes are designated for C. vernale to stabilize this concept. An additional species, Chrysotoxum orthostylum Vujić sp. nov., with distinctive male genitalia is also described. The three species share an antenna with the basoflagellomere shorter than the scape plus pedicel and terga with yellow fasciae not reaching the lateral margins. This study confirms the value of integrative approach for resolving species boundaries.

Seasonal differences in migration patterns of a soaring bird in relation to environmental conditions: a multi-scale approach

Many studies suggest that migratory birds are expected to travel more quickly during spring, when they are en route to the breeding grounds, in order to ensure a high-quality territory. Using data recorded by means of Global Positioning System satellite tags, we analysed at three temporal scales (hourly, daily and overall journey) seasonal differences in migratory performance of the booted eagle (Aquila pennata), a soaring raptor migrating between Europe and tropical Africa, taking into account environmental conditions such as wind, thermal uplift and day length. Unexpectedly, booted eagles showed higher travel rates (hourly speed, daily distance, overall migration speed and overall straightness) during autumn, even controlling for abiotic factors, probably thanks to higher hourly speeds, more straight routes and less non-travelling days during autumn. Tailwinds were the main environmental factor affecting daily distance. During spring, booted eagles migrated more quickly when flying over the Sahara desert. Our results raise new questions about which ecological and behavioural reasons promote such unexpected faster speeds in autumn and not during spring and how events occurring in very different regions can affect migratory performance, interacting with landscape characteristics, weather conditions and flight behaviour.

Dairy cattle-slurry treatment technologies: towards an integrating approach of crop production, animal feeding and environmental impacts

Dairy cattle farms have a well-known environmental impact that affects all ecological compartments: air, soil, water and biosphere [1]. Dairy cattle farming are a significant source of anthropogenic gases from enteric fermentation, manure storage and land application, mainly ammonia (NH3), nitric oxide (NO), nitrous oxide (N2O), carbon dioxide (CO2) and methane (CH4). The emission of such gases represents not only an environmental problem but also leads to energy and nitrogen (N) losses in ruminant production systems [2-5]. Several efforts are required on the development of new technologies and strategies that mitigate gaseous emissions, N losses and improve the efficiency of the energy and N cycles [6, 7]. In the Northwest of Portugal, dairy cattle production has a major impact on the economy, with strong repercussions at national scale. Therefore, our Ph.D. thesis project aims to: a) Study natural supplements as additives in the dairy cattle diet towards a decrease in GHG emissions from feeding operations; b) Compare commercial dairy cattle diets with and without additives on gaseous emissions from manure deposited in a simulated concrete floor; c) Assess the concentrations and emissions of NH3 and greenhouse gases from commercial dairy cattle facilities; d) Evaluate the effects of different additives on lowering gaseous emissions from dairy cattle excreta, using a laboratory system simulating a dairy house concrete floor.

Variations between post- and pre-harvest seasons in stunting, wasting, and Infant and Young Child Feeding (IYCF) practices among children 6-23 months of age in lowland and midland agro-ecological zones of rural Ethiopia

Introduction: Food availability and access are strongly affected by seasonality in Ethiopia. However, there are little data on seasonal variation in Infant and Young Child Feeding (IYCF) practices and malnutrition among 6-23 months old children in different agro-ecological zones of rural Ethiopia. Methods: Socio-demographic, anthropometry and IYCF indicators were assessed in post- and pre-harvest seasons among children aged 6–23 months of age randomly selected from rural villages of lowland and midland agro-ecological zones. Results: Child stunting and underweight increased from prevalence of 39.8% and 26.9% in post-harvest to 46.0% and 31.8% in pre-harvest seasons, respectively. The biggest increase in prevalence of stunting and underweight between post- and pre-harvest seasons was noted in the midland zone. Wasting decreased from 11.6% post-harvest to 8.5% pre-harvest, with the biggest decline recorded in the lowland zone. Minimum meal frequency, minimum acceptable diet and poor dietary diversity increased considerably in pre-harvest compared to post-harvest season in the lowland zone. Feeding practices and maternal age were predictors of wasting, while women’s dietary diversity and children age was predictor of child dietary diversity in both seasons. Conclusion: There is seasonal variation in malnutrition and IYCF practices among children 6-23 months of age with more pronounced effect in midland agro-ecological zone. A major contributing factor for child malnutrition may be poor feeding practices. Health information strategies focused on both IYCF practices and dietary diversity of mothers could be a sensible approach to reduce the burden of child malnutrition in rural Ethiopia.