Antimicrobial resistance in bacteria from Swiss veal calves at slaughter

Bacteria with antimicrobial resistance can be transferred from animals to humans and may compromise antimicrobial treatment in case of infection. To determine the antimicrobial resistance situation in bacteria from Swiss veal calves, faecal samples from 500 randomly selected calves originating from 129 farms were collected at four big slaughterhouses. Samples were cultured for Escherichia coli, Enterococcus sp. and Campylobacter sp. and isolated strains were tested for antimicrobial susceptibility to selected antimicrobial agents by the minimal inhibitory concentration technique using the broth microdilution method. From 100 farms, data on farm management, animal husbandry and antimicrobial treatments of the calves were collected by questionnaire. Risk factors associated with antimicrobial resistance were identified by logistic regression. In total, 467 E. coli, 413 Enterococcus sp. and 202 Campylobacter sp. were isolated. Of those, 68.7%, 98.7% and 67.8%, respectively, were resistant to at least one of the tested antimicrobial agents. Resistance was mainly observed to antimicrobials frequently used in farm animals. Prevalence of resistance to antimicrobials important for human treatment was generally low. However, a rather high number of quinupristin/dalfopristin-resistant Enterococcus faecium and ciprofloxacin-resistant Campylobacter sp. were detected. External calf purchase, large finishing groups, feeding of milk by-products and administration of antimicrobials through feed upon arrival of the animals on the farm significantly increased the risk of antimicrobial resistance at farm level. Participation in a quality assurance programme and injection of a macrolide upon arrival of the animals on the farm had a protective effect. The present study showed that veal calves may serve as a reservoir for resistant bacteria. To ensure food safety, veal calves should be included in the national monitoring programme for antimicrobial resistance in farm animals. By improving farm management and calf husbandry the prevalence of resistance may be reduced.

Intrasexual competition among females and the stabilization of a conspicuous colour polymorphism in a Lake Victoria cichlid fish

The maintenance of colour polymorphisms within populations has been a long-standing interest in evolutionary ecology. African cichlid fish contain some of the most striking known cases of this phenomenon. Intrasexual selection can be negative frequency dependent when males bias aggression towards phenotypically similar rivals, stabilizing male colour polymorphisms. We propose that where females are territorial and competitive, aggression biases in females may also promote coexistence of female morphs. We studied a polymorphic population of the cichlid fish Neochromis omnicaeruleus from Lake Victoria, in which three distinct female colour morphs coexist: one plain brown and two blotched morphs. Using simulated intruder choice tests in the laboratory, we show that wild-caught females of each morph bias aggression towards females of their own morph, suggesting that females of all three morphs may have an advantage when their morph is locally the least abundant. This mechanism may contribute to the establishment and stabilization of colour polymorphisms. Next, by crossing the morphs, we generated sisters belonging to different colour morphs. We find no sign of aggression bias in these sisters, making pleiotropy unlikely to explain the association between colour and aggression bias in wild fish, which is maintained in the face of gene flow. We conclude that female-female aggression may be one important force for stabilizing colour polymorphism in cichlid fish.

Relationship between antibiotic concentration in bone and efficacy of treatment of staphylococcal osteomyelitis in rats: azithromycin compared with clindamycin and rifampin

We examined the effect of azithromycin (CP-62,993), a new oral macrolide-like antibiotic, alone and in combination with rifampin, as treatment for experimental staphylococcal osteomyelitis. Clindamycin was used as a comparison drug. Rats (n = 10 to 15 per group) were infected by direct instillation of Staphylococcus aureus into the tibial medullary cavity. After 10 days, 21-day treatments with azithromycin (50 mg/kg of body weight, once daily, by the oral route), rifampin (20 mg/kg, once daily, subcutaneously), or clindamycin (90 mg/kg, three times daily, by the oral route) were started. The drugs were used singly or in combination (azithromycin plus rifampin or clindamycin plus rifampin). Peak azithromycin concentrations in bone were > 30 times higher than levels in serum, but the drug had little effect on final bacterial titers (5.13 +/- 0.46 log10 CFU/g of bone; for controls, 6.54 +/- 0.28 log10 CFU/g). Clindamycin was more active than azithromycin (3.26 +/- 2.14 log10 CFU/g of bone; 20% of sterilized bones), but rifampin was the most active single drug (1.5 +/- 1.92 log10 CFU/g; 53% of sterilized bones). Therapy with rifampin or clindamycin alone was associated with the emergence of resistance. Rifampin plus azithromycin (0.51 +/- 1.08 log10 CFU/g of bone; 80% of sterilized bones) and rifampin plus clindamycin (0.87 +/- 1.34 log10 CFU/g of bone; 66% of sterilized bones) were the most active regimens. Thus, azithromycin is ineffective as a single drug for the treatment of experimental staphylococcal osteomyelitis, despite high levels in bone that markedly exceeded the MIC, but it may be an attractive partner drug for rifampin.

Point mutations in the stem region and the fourth AAA domain of cytoplasmic dynein heavy chain partially suppress the phenotype of NUDF/LIS1 loss in Aspergillus nidulans

Cytoplasmic dynein performs multiple cellular tasks but its regulation remains unclear. The dynein heavy chain has a N-terminal stem that binds to other subunits and a C-terminal motor unit that contains six AAA (ATPase associated with cellular activities) domains and a microtubule-binding site located between AAA4 and AAA5. In Aspergillus nidulans, NUDF (a LIS1 homolog) functions in the dynein pathway, and two nudF6 partial suppressors were mapped to the nudA dynein heavy chain locus. Here we identified these two mutations. The nudAL1098F mutation resides in the stem region, and nudAR3086C is in the end of AAA4. These mutations partially suppress the phenotype of nudF deletion but do not suppress the phenotype exhibited by mutants of dynein intermediate chain and Arp1. Surprisingly, the stronger DeltanudF suppressor, nudAR3086C, causes an obvious decrease in the basal level of dynein's ATPase activity and an increase in dynein's distribution along microtubules. Thus, suppression of the DeltanudF phenotype may result from mechanisms other than simply the enhancement of dynein's ATPase activity. The fact that a mutation in the end of AAA4 negatively regulates dynein's ATPase activity but partially compensates for NUDF loss indicates the importance of the AAA4 domain in dynein regulation in vivo.

A feasibility study on model-based evaluation of kidney perfusion measured by means of FAIR prepared true-FISP arterial spin labeling (ASL) on a 3-T MR scanner

RATIONALE AND OBJECTIVES: A feasibility study on measuring kidney perfusion by a contrast-free magnetic resonance (MR) imaging technique is presented. MATERIALS AND METHODS: A flow-sensitive alternating inversion recovery (FAIR) prepared true fast imaging with steady-state precession (TrueFISP) arterial spin labeling sequence was used on a 3.0-T MR-scanner. The basis for quantification is a two-compartment exchange model proposed by Parkes that corrects for diverse assumptions in single-compartment standard models. RESULTS: Eleven healthy volunteers (mean age, 42.3 years; range 24-55) were examined. The calculated mean renal blood flow values for the exchange model (109 +/- 5 [medulla] and 245 +/- 11 [cortex] ml/min - 100 g) are in good agreement with the literature. Most important, the two-compartment exchange model exhibits a stabilizing effect on the evaluation of perfusion values if the finite permeability of the vessel wall and the venous outflow (fast solution) are considered: the values for the one-compartment standard model were 93 +/- 18 (medulla) and 208 +/- 37 (cortex) ml/min - 100 g. CONCLUSION: This improvement will increase the accuracy of contrast-free imaging of kidney perfusion in treatment renovascular disease.

Outcome and patients' satisfaction after functional treatment of acute lateral ankle injuries at emergency departments versus family doctor offices

BACKGROUND: In some Western countries, more and more patients seek initial treatment even for minor injuries at emergency units of hospitals. The initial evaluation and treatment as well as aftercare of these patients require large amounts of personnel and logistical resources, which are limited and costly, especially if compared to treatment by a general practitioner. In this study, we investigated whether outsourcing from our level 1 trauma center to a general practitioner has an influence on patient satisfaction and compliance. METHODS: This prospective, randomized study, included n = 100 patients who suffered from a lateral ankle ligament injury grade I-II (16, 17). After radiological exclusion of osseous lesions, the patients received early functional treatment and were shown physical therapy exercises to be done at home, without immobilization or the use of stabilizing ortheses. The patients were randomly assigned into two groups of 50 patients each: Group A (ER): Follow-up and final examination in the hospital's emergency unit. Group B (GP): Follow-up by general practitioner, final examination at hospital's emergency unit. The patients were surveyed regarding their satisfaction with the treatment and outcome of the treatment. RESULTS: Female and male patients were equally represented in both groups. The age of the patients ranged from 16 - 64 years, with a mean age of 34 years (ER) and 35 years (GP). 98% (n = 98) of all patients were satisfied with their treatment, and 93% (n = 93) were satisfied with the outcome. For these parameters no significant difference between the two groups could be noted (p = 0.7406 and 0.7631 respectively). 39% of all patients acquired stabilizing ortheses like ankle braces (Aircast, Malleoloc etc.) on their own initiative. There was a not significant tendency for more self-acquired ortheses in the group treated by general practicioners (p = 0,2669). CONCLUSION: Patients who first present at the ER with a lateral ankle ligament injury grade I-II can be referred to a general practitioner for follow-up treatment without affecting patient satisfaction regarding treatment and treatment outcome.

Posttranslational modification of the AU-rich element binding protein HuR by protein kinase Cdelta elicits angiotensin II-induced stabilization and nuclear export of cyclooxygenase 2 mRNA

The mRNA stabilizing factor HuR is involved in the posttranscriptional regulation of many genes, including that coding for cyclooxygenase 2 (COX-2). Employing RNA interference technology and actinomycin D experiments, we demonstrate that in human mesangial cells (hMC) the amplification of cytokine-induced COX-2 by angiotensin II (AngII) occurs via a HuR-mediated increase of mRNA stability. Using COX-2 promoter constructs with different portions of the 3' untranslated region of COX-2, we found that the increase in COX-2 mRNA stability is attributable to a distal class III type of AU-rich element (ARE). Likewise, the RNA immunoprecipitation assay showed AngII-induced binding of HuR to this ARE. Using the RNA pulldown assay, we demonstrate that the AngII-caused HuR assembly with COX-2 mRNA is found in free and cytoskeleton-bound polysomes indicative of an active RNP complex. Mechanistically, the increased HuR binding to COX-2-ARE by AngII is accompanied by increased nucleocytoplasmic HuR shuttling and depends on protein kinase Cdelta (PKCdelta), which physically interacts with nuclear HuR, thereby promoting its phosphorylation. Mapping of phosphorylation sites identified serines 221 and 318 as critical target sites for PKCdelta-triggered HuR phosphorylation and AngII-induced HuR export to the cytoplasm. Posttranslational modification of HuR by PKCdelta represents an important novel mode of HuR activation implied in renal COX-2 regulation.

Automated design of self-stabilization

Self-stabilization is a property of a distributed system such that, regardless of the legitimacy of its current state, the system behavior shall eventually reach a legitimate state and shall remain legitimate thereafter. The elegance of self-stabilization stems from the fact that it distinguishes distributed systems by a strong fault tolerance property against arbitrary state perturbations. The difficulty of designing and reasoning about self-stabilization has been witnessed by many researchers; most of the existing techniques for the verification and design of self-stabilization are either brute-force, or adopt manual approaches non-amenable to automation. In this dissertation, we first investigate the possibility of automatically designing self-stabilization through global state space exploration. In particular, we develop a set of heuristics for automating the addition of recovery actions to distributed protocols on various network topologies. Our heuristics equally exploit the computational power of a single workstation and the available parallelism on computer clusters. We obtain existing and new stabilizing solutions for classical protocols like maximal matching, ring coloring, mutual exclusion, leader election and agreement. Second, we consider a foundation for local reasoning about self-stabilization; i.e., study the global behavior of the distributed system by exploring the state space of just one of its components. It turns out that local reasoning about deadlocks and livelocks is possible for an interesting class of protocols whose proof of stabilization is otherwise complex. In particular, we provide necessary and sufficient conditions – verifiable in the local state space of every process – for global deadlock- and livelock-freedom of protocols on ring topologies. Local reasoning potentially circumvents two fundamental problems that complicate the automated design and verification of distributed protocols: (1) state explosion and (2) partial state information. Moreover, local proofs of convergence are independent of the number of processes in the network, thereby enabling our assertions about deadlocks and livelocks to apply on rings of arbitrary sizes without worrying about state explosion.

Computational study of internal and external condensing flows and experimental synthesis to investigate their attainability and stability in ground-based and space-based environments

This doctoral thesis presents the computational work and synthesis with experiments for internal (tube and channel geometries) as well as external (flow of a pure vapor over a horizontal plate) condensing flows. The computational work obtains accurate numerical simulations of the full two dimensional governing equations for steady and unsteady condensing flows in gravity/0g environments. This doctoral work investigates flow features, flow regimes, attainability issues, stability issues, and responses to boundary fluctuations for condensing flows in different flow situations. This research finds new features of unsteady solutions of condensing flows; reveals interesting differences in gravity and shear driven situations; and discovers novel boundary condition sensitivities of shear driven internal condensing flows. Synthesis of computational and experimental results presented here for gravity driven in-tube flows lays framework for the future two-phase component analysis in any thermal system. It is shown for both gravity and shear driven internal condensing flows that steady governing equations have unique solutions for given inlet pressure, given inlet vapor mass flow rate, and fixed cooling method for condensing surface. But unsteady equations of shear driven internal condensing flows can yield different “quasi-steady” solutions based on different specifications of exit pressure (equivalently exit mass flow rate) concurrent to the inlet pressure specification. This thesis presents a novel categorization of internal condensing flows based on their sensitivity to concurrently applied boundary (inlet and exit) conditions. The computational investigations of an external shear driven flow of vapor condensing over a horizontal plate show limits of applicability of the analytical solution. Simulations for this external condensing flow discuss its stability issues and throw light on flow regime transitions because of ever-present bottom wall vibrations. It is identified that laminar to turbulent transition for these flows can get affected by ever present bottom wall vibrations. Detailed investigations of dynamic stability analysis of this shear driven external condensing flow result in the introduction of a new variable, which characterizes the ratio of strength of the underlying stabilizing attractor to that of destabilizing vibrations. Besides development of CFD tools and computational algorithms, direct application of research done for this thesis is in effective prediction and design of two-phase components in thermal systems used in different applications. Some of the important internal condensing flow results about sensitivities to boundary fluctuations are also expected to be applicable to flow boiling phenomenon. Novel flow sensitivities discovered through this research, if employed effectively after system level analysis, will result in the development of better control strategies in ground and space based two-phase thermal systems.

Interferon-beta stabilizes barrier characteristics of the blood-brain barrier in four different species in vitro

BACKGROUND: Blood-brain barrier (BBB) breakdown is an early event in the pathogenesis of multiple sclerosis (MS). In a previous study we have found a direct stabilization of barrier characteristics after treatment of bovine brain capillary endothelial cells (BCECs) with human recombinant interferon-beta-1a (IFN-beta-1a) in an in vitro BBB model. In the present study we examined the effect of human recombinant IFN-beta-1a on the barrier properties of BCECs derived from four different species including humans to predict treatment efficacy of IFN-beta-1a in MS patients. METHODS: We used primary bovine and porcine BCECs, as well as human and murine BCEC cell lines. We investigated the influence of human recombinant IFN-beta-1a on the paracellular permeability for 3H-inulin and 14C-sucrose across monolayers of bovine, human, and murine BCECs. In addition, the transendothelial electrical resistance (TEER) was determined in in vitro systems applying porcine and murine BCECS. RESULTS: We found a stabilizing effect on the barrier characteristics of BCECs after pretreatment with IFN-beta-1a in all applied in vitro models: addition of IFN-beta-1a resulted in a significant decrease of the paracellular permeability across monolayers of human, bovine, and murine BCECs. Furthermore, the TEER was significantly increased after pretreatment of porcine and murine BCECs with IFN-beta-1a. CONCLUSION: Our data suggest that BBB stabilization by IFN-beta-1a may contribute to its beneficial effects in the treatment of MS. A human in vitro BBB model might be useful as bioassay for testing the treatment efficacy of drugs in MS.

SPACECRAFT ATTITUDE CONTROL USING MAGNETIC ACTUATORS

This report presents a study on the problem of spacecraft attitude control using magnetic actuators. Several existing approaches are reviewed and one control strategy is implemented and simulated. A time-varying feedback control law achieving inertial pointing for magnetically actuated spacecraft is implemented. The report explains the modeling of the spacecraft rigid body dynamics, kinematics and attitude control in detail. Besides the fact that control laws have been established for stabilization around local equilibrium, this report presents the results of a control law that yields a generic, global solution for attitude stabilization of a magnetically actuated spacecraft. The report also involves the use MATLAB as a tool for both modeling and simulation of the spacecraft and controller. In conclusion, the simulation outlines the performance of the controller in independently stabilizing the spacecraft in three mutually perpendicular directions.

Integrative analysis of RUNX1 downstream pathways and target genes

BACKGROUND: The RUNX1 transcription factor gene is frequently mutated in sporadic myeloid and lymphoid leukemia through translocation, point mutation or amplification. It is also responsible for a familial platelet disorder with predisposition to acute myeloid leukemia (FPD-AML). The disruption of the largely unknown biological pathways controlled by RUNX1 is likely to be responsible for the development of leukemia. We have used multiple microarray platforms and bioinformatic techniques to help identify these biological pathways to aid in the understanding of why RUNX1 mutations lead to leukemia. RESULTS: Here we report genes regulated either directly or indirectly by RUNX1 based on the study of gene expression profiles generated from 3 different human and mouse platforms. The platforms used were global gene expression profiling of: 1) cell lines with RUNX1 mutations from FPD-AML patients, 2) over-expression of RUNX1 and CBFbeta, and 3) Runx1 knockout mouse embryos using either cDNA or Affymetrix microarrays. We observe that our datasets (lists of differentially expressed genes) significantly correlate with published microarray data from sporadic AML patients with mutations in either RUNX1 or its cofactor, CBFbeta. A number of biological processes were identified among the differentially expressed genes and functional assays suggest that heterozygous RUNX1 point mutations in patients with FPD-AML impair cell proliferation, microtubule dynamics and possibly genetic stability. In addition, analysis of the regulatory regions of the differentially expressed genes has for the first time systematically identified numerous potential novel RUNX1 target genes. CONCLUSION: This work is the first large-scale study attempting to identify the genetic networks regulated by RUNX1, a master regulator in the development of the hematopoietic system and leukemia. The biological pathways and target genes controlled by RUNX1 will have considerable importance in disease progression in both familial and sporadic leukemia as well as therapeutic implications.

LIFE CYCLE ASSESSMENTS (LCAs) OF PYROLYSIS-BASED GASOLINE AND DIESEL FROM DIFFERENT REGIONAL FEEDSTOCKS: CORN STOVER, SWITCHGRASS, SUGAR CANE BAGASSE, WASTE WOOD, GUINEA GRASS, ALGAE, AND ALBIZIA

Renewable hydrocarbon biofuels are being investigated as possible alternatives to conventional liquid transportation fossil fuels like gasoline, kerosene (aviation fuel), and diesel. A diverse range of biomass feedstocks such as corn stover, sugarcane bagasse, switchgrass, waste wood, and algae, are being evaluated as candidates for pyrolysis and catalytic upgrading to produce drop-in hydrocarbon fuels. This research has developed preliminary life cycle assessments (LCA) for each feedstock-specific pathway and compared the greenhouse gas (GHG) emissions of the hydrocarbon biofuels to current fossil fuels. As a comprehensive study, this analysis attempts to account for all of the GHG emissions associated with each feedstock pathway through the entire life cycle. Emissions from all stages including feedstock production, land use change, pyrolysis, stabilizing the pyrolysis oil for transport and storage, and upgrading the stabilized pyrolysis oil to a hydrocarbon fuel are included. In addition to GHG emissions, the energy requirements and water use have been evaluated over the entire life cycle. The goal of this research is to help understand the relative advantages and disadvantages of the feedstocks and the resultant hydrocarbon biofuels based on three environmental indicators; GHG emissions, energy demand, and water utilization. Results indicate that liquid hydrocarbon biofuels produced through this pyrolysis-based pathway can achieve greenhouse gas emission savings of greater than 50% compared to petroleum fuels, thus potentially qualifying these biofuels under the US EPA RFS2 program. GHG emissions from biofuels ranged from 10.7-74.3 g/MJ from biofuels derived from sugarcane bagasse and wild algae at the extremes of this range, respectively. The cumulative energy demand (CED) shows that energy in every biofuel process is primarily from renewable biomass and the remaining energy demand is mostly from fossil fuels. The CED for biofuel range from 1.25-3.25 MJ/MJ from biofuels derived from sugarcane bagasse to wild algae respectively, while the other feedstock-derived biofuels are around 2 MJ/MJ. Water utilization is primarily from cooling water use during the pyrolysis stage if irrigation is not used during the feedstock production stage. Water use ranges from 1.7 - 17.2 gallons of water per kg of biofuel from sugarcane bagasse to open pond algae, respectively.

SYNTHESIS OF GRAPHENE AND ITS APPLICATIONS FOR DYE-SENSITIZED SOLAR CELLS

Graphene, which is a two-dimensional carbon material, exhibits unique properties that promise its potential applications in photovoltaic devices. Dye-sensitized solar cell (DSSC) is a representative of the third generation photovoltaic devices. Therefore, it is important to synthesize graphene with special structures, which possess excellent properties for dye-sensitized solar cells. This dissertation research was focused on (1) the effect of oxygen content on the structure of graphite oxide, (2) the stability of graphene oxide solution, (3) the application of graphene precipitate from graphene oxide solution as counter electrode for DSSCs, (4) the development of a novel synthesis method for the three-dimensional graphene with honeycomb-like structure, and (5) the exploration of honeycomb structured graphene (HSG) as counter electrodes for DSSCs. Graphite oxide is a crucial precursor to synthesize graphene sheets via chemical exfoliation method. The relationship between the oxygen content and the structures of graphite oxides was still not explored. In this research, the oxygen content of graphite oxide is tuned by changing the oxidation time and the effect of oxygen content on the structure of graphite oxide was evaluated. It has been found that the saturated ratio of oxygen to carbon is 0.47. The types of functional groups in graphite oxides, which are epoxy, hydroxyl, and carboxylgroups, are independent of oxygen content. However, the interplanar space and BET surface area of graphite oxide linearly increases with increasing O/C ratio. Graphene oxide (GO) can easily dissolve in water to form a stable homogeneous solution, which can be used to fabricate graphene films and graphene based composites. This work is the first research to evaluate the stability of graphene oxide solution. It has been found that the introduction of strong electrolytes (HCl, LiOH, LiCl) into GO solution can cause GO precipitation. This indicates that the electrostatic repulsion plays a critical role in stabilizing aqueous GO solution. Furthermore, the HCl-induced GO precipitation is a feasible approach to deposit GO sheets on a substrate as a Pt-free counter electrode for a dye-sensitized solar cell (DSSC), which exhibited 1.65% of power conversion efficiency. To explore broad and practical applications, large-scale synthesis with controllable integration of individual graphene sheets is essential. A novel strategy for the synthesis of graphene sheets with three-dimensional (3D) Honeycomb-like structure has been invented in this project based on a simple and novel chemical reaction (Li2O and CO to graphene and Li2CO3). The simultaneous formation of Li2CO3 with graphene not only can isolate graphene sheets from each other to prevent graphite formation during the process, but also determine the locally curved shape of graphene sheets. After removing Li2CO3, 3D graphene sheets with a honeycomb-like structure were obtained. This would be the first approach to synthesize 3D graphene sheets with a controllable shape. Furthermore, it has been demonstrated that the 3D Honeycomb-Structured Graphene (HSG) possesses excellent electrical conductivity and high catalytic activity. As a result, DSSCs with HSG counter electrodes exhibit energy conversion efficiency as high as 7.8%, which is comparable to that of an expensive noble Pt electrode.

Selection on floral traits through male fertility in a natural plant population

Most studies on selection in plants estimate female fitness components and neglect male mating success, although the latter might also be fundamental to understand adaptive evolution. Information from molecular genetic markers can be used to assess determinants of male mating success through parentage analyses. We estimated paternal selection gradients on floral traits in a large natural population of the herb Mimulus guttatus using a paternity probability model and maximum likelihood methods. This analysis revealed more significant selection gradients than a previous analysis based on regression of estimated male fertilities on floral traits. There were differences between results of univariate and multivariate analyses most likely due to the underlying covariance structure of the traits. Multivariate analysis, which corrects for the covariance structure of the traits, indicated that male mating success declined with distance from and depended on the direction to the mother plants. Moreover, there was directional selection for plants with fewer open flowers which have smaller corollas, a smaller anther-stigma separation, more red dots on the corolla and a larger fluctuating asymmetry therein. For most of these traits, however, there was also stabilizing selection indicating that there are intermediate optima for these traits. The large number of significant selection gradients in this study shows that even in relatively large natural populations where not all males can be sampled, it is possible to detect significant paternal selection gradients, and that such studies can give us valuable information required to better understand adaptive plant evolution.