Sectoral impacts on the diversity of carabid beetles (Coleoptera: Carabidae) and their associated ecosystem service of biocontrol

This thesis aimed to provide an understanding of how human-induced changes in the economic sectors of agriculture and transport affect carabid diversity, potential carabidmediated biocontrol and predator-pest interactions. The research involved both observational and manipulative laboratory and field-based studies. Observational research consisted of two large-scale investigations of (1) the impact of Miscanthus and oilseed rape production (n=45) and (2) the impact of horticultural and ecological based landscaping of roadside verges (n=64). This research is the first record of carabid diversity, potential biocontrol and community assemblage with respect to bioenergy crop production and roadside landscaping in an Irish context and it is also an important addition to the limited knowledge of carabid populations in these ecosystems internationally. Manipulative work involved the examination of the role predator identity, diversity and biomass play in the suppression of pollen beetle larvae (an economically damaging insect pest of oilseed rape in Europe), using a novel experimental design called ‘simplex’. To complement this research, an additional field study on the impact of low and high oilseed rape pesticide management on carabid species richness and abundance, and crop yield, was also conducted. This research is a great contribution to the existing understanding of what constitutes the important components of predator biodiversity and expands the knowledge of the usefulness of carabid predators in the context of pollen beetle larvae control. In particular, the work shows that the abundance or biomass of beetles has an effect that is far larger than the effect of diversity on the capacity of beetles to consume prey. In turn, the field study showed that pesticide applications had little impact on yield, or carabid richness, but that carabid abundance/biomass declined drastically. The work provides compelling evidence that management practices erode the useful components of biodiversity that are essential for the delivery of biocontrol services.

Quantifying heritable variation in fitness-related traits of wild, farmed and hybrid Atlantic salmon families in a wild river environment

Farmed fish are typically genetically different from wild conspecifics. Escapees from fish farms may contribute one-way gene flow from farm to wild gene pools, which can depress population productivity, dilute local adaptations and disrupt coadapted gene complexes. Here, we reanalyse data from two experiments (McGinnity et al., 1997, 2003) where performance of Atlantic salmon (Salmo salar) progeny originating from experimental crosses between farm and wild parents (in three different cohorts) were measured in a natural stream under common garden conditions. Previous published analyses focussed on group-level differences but did not account for pedigree structure, as we do here using modern mixed-effect models. Offspring with one or two farm parents exhibited poorer survival in their first and second year of life compared with those with two wild parents and these group-level inferences were robust to excluding outlier families. Variation in performance among farm, hybrid and wild families was generally similar in magnitude. Farm offspring were generally larger at all life stages examined than wild offspring, but the differences were moderate (5–20%) and similar in magnitude in the wild versus hatchery environments. Quantitative genetic analyses conducted using a Bayesian framework revealed moderate heritability in juvenile fork length and mass and positive genetic correlations (>0.85) between these morphological traits. Our study confirms (using more rigorous statistical techniques) previous studies showing that offspring of wild fish invariably have higher fitness and contributes fresh insights into family-level variation in performance of farm, wild and hybrid Atlantic salmon families in the wild. It also adds to a small, but growing, number of studies that estimate key evolutionary parameters in wild salmonid populations. Such information is vital in modelling the impacts of introgression by escaped farm salmon.

Integrating personality research and animal contest theory:Aggressiveness in the green swordtail xiphophorus helleri

Aggression occurs when individuals compete over limiting resources. While theoretical studies have long placed a strong emphasis on context-specificity of aggression, there is increasing recognition that consistent behavioural differences exist among individuals, and that aggressiveness may be an important component of individual personality. Though empirical studies tend to focus on one aspect or the other, we suggest there is merit in modelling both within- and among-individual variation in agonistic behaviour simultaneously. Here, we demonstrate how this can be achieved using multivariate linear mixed effect models. Using data from repeated mirror trials and dyadic interactions of male green swordtails, Xiphophorus helleri, we show repeatable components of (co)variation in a suite of agonistic behaviour that is broadly consistent with a major axis of variation in aggressiveness. We also show that observed focal behaviour is dependent on opponent effects, which can themselves be repeatable but were more generally found to be context specific. In particular, our models show that within-individual variation in agonistic behaviour is explained, at least in part, by the relative size of a live opponent as predicted by contest theory. Finally, we suggest several additional applications of the multivariate models demonstrated here. These include testing the recently queried functional equivalence of alternative experimental approaches, (e.g., mirror trials, dyadic interaction tests) for assaying individual aggressiveness. © 2011 Wilson et al.

A cross-sectional study of reproductive performance of smallholder dairy cows in coastal Tanzania

A retrospective cross-sectional study was conducted on 200 randomly selected smallholder farms from a mixed dairy farming system in Tanga, Tanzania, between January and April 1999. We estimated the frequency and determinants of long calving interval (LCI), retention of fetal membrane (RFM), dystocia, and abortion in smallholder crossbred cattle and explored birth trends. The mean calving interval was 500 days and birth rate was 65 per 100 cow-years. Dystocia was reported to affect 58% of calvings, and 17.2% of animals suffered RFM. Using mixed effect models, the variables associated with LCI, RFM and dystocia were breed, level of exotic blood and condition score. Zebu breeding was associated with LCI (odds ratio (OR) = 2.3, p = 0.041) and Friesian breeding with lower odds for RF (OR = 0.26, p = 0.020). Animals with higher levels of exotic blood had lower odds for evidence of dystocia (OR = 0.45, p = 0.021). Evidence of dystocia was significantly associated with poor condition score (beta = -1.10, p = 0.001). Our observations suggest that LCIs are common in smallholder dairy farms in this region and a likely source of economic loss. Dystocia, RFM, poor condition score and mineral deficiency were common problems and were possibly linked to LCI.

Physical activity attenuates the body mass index-increasing influence of genetic variation in the FTO gene 1-3.

BACKGROUND: Intronic variation in the FTO (fat mass and obesity-associated) gene has been unequivocally associated with increased body mass index (BMI; in kg/m(2)) and the risk of obesity in populations of different ethnicity. OBJECTIVE: We examined whether this robust genetic predisposition to obesity can be attenuated by being more physically active. DESIGN: The FTO variant rs1121980 was genotyped in 20,374 participants (39-79 y of age) from the European Prospective Investigation into Cancer and Nutrition-Norfolk Study, an ethnically homogeneous population-based cohort. Physical activity (PA) was assessed with a validated self-reported questionnaire. The interaction between rs1121980 and PA on BMI and waist circumference (WC) was examined by including the interaction term in mixed-effect models. RESULTS: We confirmed that the risk (T) allele of rs1121980 was significantly associated with BMI (0.31-unit increase per allele; P < 0.001) and WC (0.77-cm increase per allele; P < 0.001). The PA level attenuated the effect of rs1121980 on BMI and WC; ie, whereas in active individuals the risk allele increased BMI by 0.25 per allele, the increase in BMI was significantly (P for interaction = 0.004) more pronounced (76%) in inactive individuals (0.44 per risk allele). We observed similar effects for WC (P for interaction = 0.02): the risk allele increased WC by 1.04 cm per allele in inactive individuals but by only 0.64 cm in active individuals. CONCLUSIONS: Our results showed that PA attenuates the effect of the FTO rs1121980 genotype on BMI and WC. This observation has important public health implications because we showed that a genetic susceptibility to obesity induced by FTO variation can be overcome, at least in part, by adopting a physically active lifestyle.

HPA axis related genes and response to psychological therapies: genetics and epigenetics

Background Hypothalamic–pituitary–adrenal (HPA) axis functioning has been implicated in the development of stress-related psychiatric diagnoses and response to adverse life experiences. This study aimed to investigate the association between genetic and epigenetics in HPA axis and response to cognitive behavior therapy (CBT). Methods Children with anxiety disorders were recruited into the Genes for Treatment project (GxT, N = 1,152). Polymorphisms of FKBP5 and GR were analyzed for association with response to CBT. Percentage DNA methylation at the FKBP5 and GR promoter regions was measured before and after CBT in a subset (n = 98). Linear mixed effect models were used to investigate the relationship between genotype, DNA methylation, and change in primary anxiety disorder severity (treatment response). Results Treatment response was not associated with FKBP5 and GR polymorphisms, or pretreatment percentage DNA methylation. However, change in FKBP5 DNA methylation was nominally significantly associated with treatment response. Participants who demonstrated the greatest reduction in severity decreased in percentage DNA methylation during treatment, whereas those with little/no reduction in severity increased in percentage DNA methylation. This effect was driven by those with one or more FKBP5 risk alleles, with no association seen in those with no FKBP5 risk alleles. No significant association was found between GR methylation and response. Conclusions Allele-specific change in FKBP5 methylation was associated with treatment response. This is the largest study to date investigating the role of HPA axis related genes in response to a psychological therapy. Furthermore, this is the first study to demonstrate that DNA methylation changes may be associated with response to psychological therapies in a genotype-dependent manner.

Seasonal- and sex-specific correlations between dispersal and exploratory behaviour in the great tit

While the importance of personality in explaining individual dispersal strategies is increasingly recognized, limited information is still available on how patterns of personality-dependent dispersal may develop, hampering our understanding of the ecological significance of behavioural dispersal syndromes. Here, we examine the relative importance of personality at different stages of dispersal in the great tit, by analysing the sex-specific relationship between exploratory behaviour (EB; quantified in a novel environment) and dispersal distances in different seasons over the course of the first year of life (summer, autumn, winter, and until the first breeding attempt). In females, we found that EB was an important predictor of dispersal distances in summer and autumn, but only a weak to moderate correlation remained for females captured in winter or for natal dispersal distances based on first breeding records. We obtained a contrasting pattern at the population level, whereby male (but not female) immigrants captured in summer and autumn had higher EB scores than locally born birds, while this was not the case in birds captured in winter and those recruited as breeders into the population. In addition to providing further evidence for the existence of a behaviour dispersal syndrome in birds, our results show that correlations between EB and dispersal appear strongest at the early stages of the dispersal process, rather than being developed gradually. These findings show the importance of analysing the effect of phenotypic attributes on dispersal across different stages of the dispersal phenomenon and in each sex separately.

Influence of intrinsic variation on foraging behaviour of adult female Australian fur seals

Phenotypic variation and individual experience can create behavioural and/or dietary variation within a population. This may reduce intra-specific competition, creating a buffer to environmental change. This study examined how intrinsic variation affects foraging behaviour of Australian fur seals. Foraging movements of 29 female Australian fur seals were recorded using FastLoc GPS and dive behaviour recorders. For each individual, body mass, flipper length and axis length were recorded, a tooth was sampled to determine age, and milk was collected for diet analysis. Clustering of fatty acid dietary analysis revealed 5 distinct groups in the population. Behaviour was described using 19 indices, which were then reduced to 7 principal components (>80% of the behavioural variation). Bayesian mixed effect models were developed to describe the relationship between these components and intrinsic variation. No association was found between diet and age or body shape; however, age had a negative relationship with component 1 (27% of variation). Older females spent less time at-sea and foraged nearer to the colony. Age had an effect on component 5 (7% of variation), which represented haul-outs and dive depth; older females made fewer visits to haul-out sites and dived deeper to the benthos. This suggests that as animals age they are able to utilise prior knowledge to exploit nearby foraging sites that younger animals are either unaware of, or have yet to gain the experience required to efficiently utilise. Mass had a negative effect on components representing the directedness of a foraging trip, suggesting heavier individuals were more likely to travel directly to a foraging site.

Growth response to antiretroviral treatment in HIV-infected children: a cohort study from Lilongwe, Malawi

Objective  Malnutrition is common in HIV-infected children in Africa and an indication for antiretroviral treatment (ART). We examined anthropometric status and response to ART in children treated at a large public-sector clinic in Malawi. Methods  All children aged <15 years who started ART between January 2001 and December 2006 were included and followed until March 2008. Weight and height were measured at regular intervals from 1 year before to 2 years after the start of ART. Sex- and age-standardized z-scores were calculated for weight-for-age (WAZ) and height-for-age (HAZ). Predictors of growth were identified in multivariable mixed-effect models. Results  A total of 497 children started ART and were followed for 972 person-years. Median age (interquartile range; IQR) was 8 years (4–11 years). Most children were underweight (52% of children), stunted (69%), in advanced clinical stages (94% in WHO stages 3 or 4) and had severe immunodeficiency (77%). After starting ART, median (IQR) WAZ and HAZ increased from −2.1 (−2.7 to −1.3) and −2.6 (−3.6 to −1.8) to −1.4 (−2.1 to −0.8) and −1.8 (−2.4 to −1.1) at 24 months, respectively (P < 0.001). In multivariable models, baseline WAZ and HAZ scores were the most important determinants of growth trajectories on ART. Conclusions  Despite a sustained growth response to ART among children remaining on therapy, normal values were not reached. Interventions leading to earlier HIV diagnosis and initiation of treatment could improve growth response.

An in vitro study on spontaneous myometrial contractility in the mare during estrus and diestrus

Uterine smooth muscle specimens were collected from euthanatized mares in estrus and diestrus. Longitudinal and circular specimens were mounted in organ baths and the signals transcribed to a Grass polygraph. After equilibration time and 2 g preload, their physiologic isometric contractility was recorded for a continuous 2.0 h. Area under the curve, frequency and time occupied by contractions were studied. Differences between cycle phases, between muscle layers, and over the recorded time periods were statistically evaluated using linear mixed-effect models. In the mare, physiologic contractility of the uterus decreased significantly over time for all variables evaluated (time as covariate on a continuous scale). For area under the curve, there was a significant effect of muscle layer (longitudinal > circular). For frequency, higher values were recorded in estrus for circular smooth muscle layer, whereas higher values were seen in longitudinal smooth muscle layers during diestrus. In longitudinal layer and in diestrus, more time was occupied by contractions than in circular layer, and in estrus. This study is describing physiologic myometrial motility in the organ bath depending on cycle phase.

Zidovudine impairs immunological recovery on first-line antiretroviral therapy: collaborative analysis of cohort studies in southern Africa

OBJECTIVES Zidovudine (ZDV) is recommended for first-line antiretroviral therapy (ART) in resource-limited settings. ZDV may, however, lead to anemia and impaired immunological response. We compared CD4+ cell counts over 5 years between patients starting ART with and without ZDV in southern Africa. DESIGN Cohort study. METHODS Patients aged at least 16 years who started first-line ART in South Africa, Botswana, Zambia, or Lesotho were included. We used linear mixed-effect models to compare CD4+ cell count trajectories between patients on ZDV-containing regimens and patients on other regimens, censoring follow-up at first treatment change. Impaired immunological recovery, defined as a CD4+ cell count below 100 cells/μl at 1 year, was assessed in logistic regression. Analyses were adjusted for baseline CD4+ cell count and hemoglobin level, age, sex, type of regimen, viral load monitoring, and calendar year. RESULTS A total of 72,597 patients starting ART, including 19,758 (27.2%) on ZDV, were analyzed. Patients on ZDV had higher CD4+ cell counts (150 vs.128 cells/μl) and hemoglobin level (12.0 vs. 11.0 g/dl) at baseline, and were less likely to be women than those on other regimens. Adjusted differences in CD4+ cell counts between regimens containing and not containing ZDV were -16 cells/μl [95% confidence interval (CI) -18 to -14] at 1 year and -56 cells/μl (95% CI -59 to -52) at 5 years. Impaired immunological recovery was more likely with ZDV compared to other regimens (odds ratio 1.40, 95% CI 1.22-1.61). CONCLUSION In southern Africa, ZDV is associated with inferior immunological recovery compared to other backbones. Replacing ZDV with another nucleoside reverse transcriptase inhibitor could avoid unnecessary switches to second-line ART.

Impact of Baseline Renal Function on Tenofovir-containing Antiretroviral Therapy Outcomes in Zambia

Background. Although tenofovir (TDF) use has increased as part of first-line antiretroviral therapy (ART) across sub-Saharan Africa, renal outcomes among patients receiving TDF remain poorly understood. We assessed changes in renal function and mortality in patients starting TDF- or non-TDF-containing ART in Lusaka, Zambia. Methods. We included patients aged ≥16 years who started ART from 2007 onward, with documented baseline weight and serum creatinine. Renal dysfunction was categorized as mild (eGFR 60-89 mL/min), moderate (30-59 mL/min) or severe (<30 mL/min) using the CKD-EPI formula. Differences in eGFR during ART were analyzed using linear mixed-effect models, the odds of developing moderate or severe eGFR decrease with logistic regression and mortality with competing risk regression. Results. We included 62,230 adults, of which 38,716 (62%) initiated a TDF-based regimen. The proportion with moderate or severe renal dysfunction at baseline was lower in the TDF compared to the non-TDF group (1.9% vs. 4.0%). Among patients with no or mild renal dysfunction, those on TDF were more likely to develop moderate (adjusted OR: 3.11; 95%CI: 2.52-3.87) or severe eGFR decrease (adjusted OR: 2.43; 95%CI: 1.80-3.28), although the incidence of such episodes was low. Among patients with moderate or severe renal dysfunction at baseline, renal function improved independently of ART regimen and mortality was similar in both treatment groups. Conclusions. TDF use did not attenuate renal function recovery or increase mortality in patients with renal dysfunction. Further studies are needed to determine the role of routine renal function monitoring before and during ART use in Africa.

Habitat Choice and Female Preference in a Polymorphic Stickleback Population

Background: A small pond, c. 90 years old, near Bern, Switzerland contains a population of threespine stickleback (Gasterosteus aculeatus) with two distinct male phenotypes. Males of one type are large, and red, and nest in the shallow littoral zone. The males of the other are small and orange, and nest offshore at slightly greater depth. The females in this population are phenotypically highly variable but cannot easily be assigned to either male type. Question: Is the existence of two sympatric male morphs maintained by substrate-associated male nest site choice and facilitated by female mate preferences? Organisms: Male stickleback caught individually at their breeding sites. Females caught with minnow traps. Methods: In experimental tanks, we simulated the slope and substrate of the two nesting habitats. We then placed individual males in a tank and observed in which habitat the male would build his nest. In a simultaneous two-stimulus choice design, we gave females the choice between a large, red male and a small, orange one. We measured female morphology and used linear mixed effect models to determine whether female preference correlated with female morphology. Results: Both red and orange males preferred nesting in the habitat that simulated the slightly deeper offshore condition. This is the habitat occupied by the small, orange males in the pond itself. The proportion of females that chose a small orange male was similar to that which chose a large red male. Several aspects of female phenotype correlated with the male type that a female preferred.

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.

A Poisson mixed model with nonnormal random effect distribution

In this paper, we propose a random intercept Poisson model in which the random effect is assumed to follow a generalized log-gamma (GLG) distribution. This random effect accommodates (or captures) the overdispersion in the counts and induces within-cluster correlation. We derive the first two moments for the marginal distribution as well as the intraclass correlation. Even though numerical integration methods are, in general, required for deriving the marginal models, we obtain the multivariate negative binomial model from a particular parameter setting of the hierarchical model. An iterative process is derived for obtaining the maximum likelihood estimates for the parameters in the multivariate negative binomial model. Residual analysis is proposed and two applications with real data are given for illustration. (C) 2011 Elsevier B.V. All rights reserved.