Regulating Mechanisms of Gonadal and Pituitary Tumorigenesis in Mice Producing Human Chorionic Gonadotropin

Human chorionic gonadotropin (hCG) and luteinizing hormone (LH) are structurally and functionally similar glycoprotein hormones acting through the same luteinizing hormone chorionic gonadotropin receptor (LHCGR). The functions of LH in reproduction and hCG in pregnancy are well known. Recently, the expression of LHCGR has been found in many nongonadal tissues and cancers, and this has raised the question of whether LH/hCG could affect the function or tumorigenesis of these nongonadal tissues. We have also previously generated an hCG expressing mouse model presenting nongonadal phenotypes. Using this model it is possible to improve our understanding of nongonadal action of highly elevated LH/hCG. In the current study, we analyzed the effect of moderately and highly elevated hCG levels on male reproductive development and function. The main finding was the appearance of fetal Leydig cell (FLC) adenomas in prepubertal males. However, the development and differentiation of FLCs were not significantly affected. We also show that the function of hCG is different in FLCs and in adult Leydig cells (ALC), because in the latter cells hCG was not able to induce tumorigenesis. In FLCs, LHCGR is not desensitized or downregulated upon ligand binding. In this study, we found that the testicular expression of two G protein-coupled receptor kinases responsible for receptor desensitization or downregulation is increased in adult testis. Results suggest that the lack of LHCGR desensitization or downregulation in FLCs protect testosterone (Te) synthesis, but also predispose FLCs for LH/hCG induced adenomas. However, all the hCG induced nongonadal changes observed in male mice were possible to explain by the elevated Te level found in these males. Our findings indicate that the direct nongonadal effects of elevated LH/hCG in males are not pathophysiologically significant. In female mice, we showed that an elevated hCG level was able to induce gonadal tumorigenesis. hCG also induced the formation of pituitary adenomas (PA), but the mechanism was indirect. Furthermore, we found two new potential risk factors and a novel hormonally induced mechanism for PAs. Increased progesterone (P) levels in the presence of physiological estradiol (E2) levels induced the formation of PAs in female mice. E2 and P induced the expression and nuclear localization of a known cell-cycle regulator, cyclin D1. A calorie restricted diet was also able to prevent the formation of PAs, suggesting that obesity is able to promote the formation of PAs. Hormone replacement therapy after gonadectomy and hormone antagonist therapy showed that the nongonadal phenotypes observed in hCG expressing female mice were due to ovarian hyperstimulation. A slight adrenal phenotype was evident even after gonadectomy in hCG expressing females, but E2 and P replacement was able to induce a similar phenotype in WT females without elevated LH/hCG action. In conclusion, we showed that the direct effects of elevated hCG/LH action are limited only to the gonads of both sexes. The nongonadal phenotypes observed in hCG expressing mice were due to the indirect, gonadal hormone mediated effects of elevated hCG. Therefore, the gonads are the only physiologically significant direct targets of LHCGR signalling.

Salmonella enterica– Mechanisms of fluoroquinolone and macrolide resistance

Salmonella enterica – Fluorokinoloni- ja makrolidiresistenssimekanisimit Vakavia salmonellainfektioita on pitkään hoidettu fluorokinoloniantibiooteilla, kuten siprofloksasiinilla. Fluorokinolonien runsas käyttö niin ihmisillä kuin eläimilläkin on kuitenkin johtanut fluorokinoloniresistenttien salmonellakantojen lisääntymiseen. Vuoteen 2002 asti kaikki matalan tason fluorokinoloniresistenssiä ilmentävät salmonellakannat olivat resistenttejä nalidiksiinihapolle, joka on vanha ensimmäisen polven kinoloniantibiootti jota ei enää käytetä infektioiden hoidossa. Vuonna 2003 havaitsimme aivan uudentyyppisen resistenssifenotyypin salmonelloissa. Kaikki uuden fenotyypin kannat osoittivat matalaa fluorokinoloniresistenssiä (MIC ≥0.125 mg/L), mutta useat kannat olivat yllättäen aikaisempaa herkempiä nalidiksiinihapolle (MIC ≤32 mg/L). Ilmiöllä on suuri merkitys salmonellan antibioottiherkkyyksien määrittämisessä, sillä jos kanta on ollut nalidiksiinihapolle herkkä, sitä on pidetty herkkänä myös fluorokinoloneille. Väitöskirjatyössä määritettiin vuosina 2003–2007 Suomessa kerättyjen kotimaisten ja ulkomaalaisten S. enterica -kantojen fluorokinoloniresistenssiä sekä tutkittiin uuden salmonellafenotyypin epidemiologiaa ja resistenssimekanismeja. Lisäksi tutkittiin salmonellan hoidossa mahdollisesti käyttökelpoisen makrolidiantibioottijohdannaisen, atsitromysiinin tehoa salmonelloihin ja erityisesti matalaa fluorokinoloniresistenssiä ilmentäviin kantoihin. Tutkimuksessa havaittiin, että matalaa fluorokinoloniresistenssiä osoittavien salmonellakantojen määrä vähenee. Lasku oli voimakkainta Kaakkois-Aasiasta tuoduissa kannoissa. Uusi resistenssifenotyyppi on plasmidivälitteinen ja qnr-geenit olivat ainoa plasmidivälitteinen kinoloniresistenssimekanismi, joka kannoista löydettiin. Myöskään kromosomaalisten gyrA, gyrB ja parE -geenien QRDR-alueelta ei löydetty fluorokinoloniresistenssiä aiheuttavia mutaatioita. Transformaatiolla osoitettiin qnr-plasmidien olevan siirtyviä ja uusi resistenssifenotyyppi saatiin ilmennettyä myös herkässä vastaanottajakannassa. Nämä tulokset osoittavat, että vaikka S. enterican qnr-fenotyyppi on toistaiseksi levinnyt pääasiassa Kaakkois-Aasiaan, se siirtyy helposti bakteerista toiseen ja tulee todennäköisesti aiheuttamaan hoito-ongelmia myös muualla maailmassa. Uudentyyppinen qnr-fenotyyppi voi olla vaikea havaita perinteisellä herkkyysmäärityksellä. Siksi laboratorioissa tulisi aina määrittää sekä siprofloksasiiniettä nalidiksiinihappoherkkyydet. Atsitromysiinin osoitettiin olevan herkkyysmääritysten mukaan tehokas salmonelloja kohtaan mukaanlukien matala-asteista fluorokinoloniresistenssiä ilmentävät bakteerikannat.

The intrinsic mechanisms of softwood fiber damage in brown stock fiber line unit operations Lappeenranta 2010

The effects of pulp processing on softwood fiber properties strongly influence the properties of wet and dry paper webs. Pulp strength delivery studies have provided observations that much of the strength potential of long fibered pulp is lost during brown stock fiber line operations where the pulp is merely washed and transferred to the subsequent processing stages. The objective of this work was to study the intrinsic mechanisms which maycause fiber damage in the different unit operations of modern softwood brown stock processing. The work was conducted by studying the effects of industrial machinery on pulp properties with some actions of unit operations simulated in laboratory scale devices under controlled conditions. An optical imaging system was created and used to study the orientation of fibers in the internal flows during pulp fluidization in mixers and the passage of fibers through the screen openings during screening. The qualitative changes in fibers were evaluated with existing and standardized techniques. The results showed that each process stage has its characteristic effects on fiber properties: Pulp washing and mat formation in displacement washers introduced fiber deformations especially if the fibers entering the stage were intact, but it did not decrease the pulp strength properties. However, storage chests and pulp transfer after displacement washers contributed to strength deterioration. Pulp screening proved to be quite gentle, having the potential of slightly evening out fiber deformations from very deformed pulps and vice versa inflicting a marginal increase in the deformation indices if the fibers were previously intact. Pulp mixing in fluidizing industrial mixers did not have detrimental effects on pulp strength and had the potential of slightly evening out the deformations, provided that the intensity of fluidization was high enough to allow fiber orientation with the flow and that the time of mixing was short. The chemical and mechanical actions of oxygen delignification had two distinct effects on pulp properties: chemical treatment clearly reduced pulp strength with and without mechanical treatment, and the mechanical actions of process machinery introduced more conformability to pulp fibers, but did not clearly contribute to a further decrease in pulp strength. The chemical composition of fibers entering the oxygen stage was also found to affect the susceptibility of fibers to damage during oxygen delignification. Fibers with the smallest content of xylan were found to be more prone to irreversibledeformations accompanied with a lower tensile strength of the pulp. Fibers poor in glucomannan exhibited a lower fiber strength while wet after oxygen delignification as compared to the reference pulp. Pulps with the smallest lignin content on the other hand exhibited improved strength properties as compared to the references.

Kinetic modeling of mechanisms of industrially important organic reactions in gas and liquid phase

This dissertation is based on 5 articles which deal with reaction mechanisms of the following selected industrially important organic reactions: 1. dehydrocyclization of n-butylbenzene to produce naphthalene 2. dehydrocyclization of 1-(p-tolyl)-2-methylbutane (MB) to produce 2,6-dimethylnaphthalene 3. esterification of neopentyl glycol (NPG) with different carboxylic acids to produce monoesters 4. skeletal isomerization of 1-pentene to produce 2-methyl-1-butene and 2-methyl-2-butene The results of initial- and integral-rate experiments of n-butylbenzene dehydrocyclization over selfmade chromia/alumina catalyst were applied when investigating reaction 2. Reaction 2 was performed using commercial chromia/alumina of different acidity, platina on silica and vanadium/calcium/alumina as catalysts. On all catalysts used for the dehydrocyclization, major reactions were fragmentation of MB and 1-(p-tolyl)-2-methylbutenes (MBes), dehydrogenation of MB, double bond transfer, hydrogenation and 1,6-cyclization of MBes. Minor reactions were 1,5-cyclization of MBes and methyl group fragmentation of 1,6- cyclization products. Esterification reactions of NPG were performed using three different carboxylic acids: propionic, isobutyric and 2-ethylhexanoic acid. Commercial heterogeneous gellular (Dowex 50WX2), macroreticular (Amberlyst 15) type resins and homogeneous para-toluene sulfonic acid were used as catalysts. At first NPG reacted with carboxylic acids to form corresponding monoester and water. Then monoester esterified with carboxylic acid to form corresponding diester. In disproportionation reaction two monoester molecules formed NPG and corresponding diester. All these three reactions can attain equilibrium. Concerning esterification, water was removed from the reactor in order to prevent backward reaction. Skeletal isomerization experiments of 1-pentene were performed over HZSM-22 catalyst. Isomerization reactions of three different kind were detected: double bond, cis-trans and skeletal isomerization. Minor side reaction were dimerization and fragmentation. Monomolecular and bimolecular reaction mechanisms for skeletal isomerization explained experimental results almost equally well. Pseudohomogeneous kinetic parameters of reactions 1 and 2 were estimated by usual least squares fitting. Concerning reactions 3 and 4 kinetic parameters were estimated by the leastsquares method, but also the possible cross-correlation and identifiability of parameters were determined using Markov chain Monte Carlo (MCMC) method. Finally using MCMC method, the estimation of model parameters and predictions were performed according to the Bayesian paradigm. According to the fitting results suggested reaction mechanisms explained experimental results rather well. When the possible cross-correlation and identifiability of parameters (Reactions 3 and 4) were determined using MCMC method, the parameters identified well, and no pathological cross-correlation could be seen between any parameter pair.

Role and regulatory mechanisms of heat shock factor 2

Protein homeostasis is essential for cells to prosper and survive. Various forms of stress, such as elevated temperatures, oxidative stress, heavy metals or bacterial infections cause protein damage, which might lead to improper folding and formation of toxic protein aggregates. Protein aggregation is associated with serious pathological conditions such as Alzheimer’s and Huntington’s disease. The heat shock response is a defense mechanism that protects the cell against protein-damaging stress. Its ancient origin and high conservation among eukaryotes suggest that the response is crucial for survival. The main regulator of the heat shock response is the transcription factor heat shock factor 1 (HSF1), which induces transcription of genes encoding protective molecular chaperones. In vertebrates, a family of four HSFs exists (HSF1-4), with versatile functions not only in coping with acute stress, but also in development, longevity and cancer. Thus, knowledge of the HSFs will aid in our understanding on how cells survive suboptimal circumstances, but will also provide insights into normal physiological processes as well as diseaseassociated conditions. In this study, the function and regulation of HSF2 have been investigated. Earlier gene inactivation experiments in mice have revealed roles for HSF2 in development, particularly in corticogenesis and spermatogenesis. Here, we demonstrate that HSF2 holds a role also in the heat shock response and influences stress-induced expression of heat shock proteins. Intriguingly, DNA-binding activity of HSF2 upon stress was dependent on the presence of intact HSF1, suggesting functional interplay between HSF1 and HSF2. The underlying mechanism for this phenomenon could be configuration of heterotrimers between the two factors, a possibility that was experimentally verified. By changing the levels of HSF2, the expression of HSF1-HSF2 heterotrimer target genes was altered, implementing HSF2 as a modulator of HSF-mediated transcription. The results further indicate that HSF2 activity is dependent on its concentration, which led us to ask the question of how accurate HSF2 levels are achieved. Using mouse spermatogenesis as a model system, HSF2 was found to be under direct control of miR-18, a miRNA belonging to the miR-17~92 cluster/Oncomir-1 and whose physiological function had remained unclear. Investigations on spermatogenesis are severely hampered by the lack of cell systems that would mimic the complex differentiation processes that constitute male germ cell development. Therefore, to verify that HSF2 is regulated by miR-18 in spermatogenesis, a novel method named T-GIST (Transfection of Germ cells in Intact Seminiferous Tubules) was developed. Employing this method, the functional consequences of miR-18-mediated regulation in vivo were demonstrated; inhibition of miR- 18 led to increased expression of HSF2 and altered the expression of HSF2 target genes Ssty2 and Speer4a. Consequently, the results link miR-18 to HSF2-mediated processes such as germ cell maturation and quality control and provide miR-18 with a physiological role in gene expression during spermatogenesis.Taken together, this study presents compelling evidence that HSF2 is a transcriptional regulator in the heat shock response and establishes the concept of physical interplay between HSF2 and HSF1 and functional consequences thereof. This is also the first study describing miRNA-mediated regulation of an HSF.

Dissecting the molecular mechanisms of breast cancer bone metastasis for therapeutic targeting

Breast cancer that has metastasized to bone is currently an incurable disease, causing significant morbidity and mortality. The aim of this thesis work was to elucidate molecular mechanisms of bone metastasis and thereby gain insights into novel therapeutic approaches. First, we found that L‐serine biosynthesis genes, phosphoglycerate dehydrogenase (PHGDH), phosphoserine aminotransferase 1 (PSAT1) and phosphoserine phosphatase (PSPH), were up‐regulated in highly bone metastatic MDA‐MB‐231(SA) cells as compared with the parental breast cancer cell line. Knockdown of serine biosynthesis inhibited proliferation of MDA‐MB‐231(SA) cells, and L‐serine was essential for the formation of bone resorbing osteoclasts. Clinical data demonstrated that high expression of PHGDH and PSAT1 was associated with decreased relapse‐free and overall survival and with features typical of poor outcome in breast cancer. Second, RNA interference screening pointed out heparan sulfate 6‐O‐sulfotransferase 2 (HS6ST2) as a critical gene for transforming growth factor β (TGF‐β)‐induced interleukin 11 (IL‐11) production in MDA‐MB‐231(SA) cells. Exogenous heparan sulfate glycosaminoglycans heparin and K5‐NSOS also inhibited TGF‐β‐induced IL‐11 production in MDA‐MB‐231(SA) cells. Furthermore, K5‐NSOS decreased osteolytic lesion area and tumor burden in bone in mice. Third, we discovered that the microRNAs miR‐204, ‐211 and ‐379 inhibited IL‐11 expression in MDA‐MB‐231(SA) cells through direct targeting of the IL‐11 mRNA. MiR‐379 also inhibited Smad‐mediated signaling. Gene expression profiling of miR‐204 and ‐379 transfected cells indicated that these microRNAs down‐regulate several bone metastasis‐relevant genes, including prostaglandin‐endoperoxide synthase 2 (PTGS2). Taken together, this study identified three potential treatment strategies for bone metastatic breast cancer: inhibition of serine biosynthesis, heparan sulfate glycosaminoglycans and restoration of miR‐204/‐211/‐379.

Molecular Mechanisms of Sexual Dimorphism in Threespine Stickleback

Sexual dimorphism is commonly understood as differences in external features, such as morphological features or coloration. However, it can more broadly encompass behavior and physiology and at the core of these differences is the genetic mechanism – mRNA and protein expression. How, and which, molecular mechanisms influence sexually dimorphic features is not well understood thus far. DNA, RNA and proteins are the template required to create the phenotype of an individual, and they are connected to each other via processes of transcription and translation. As the genome of males and females are almost identical with the exception of the few genes on the sex chromosome or the sex-determining alleles (in the case of organisms without sex chromosomes), it is likely that many of the downstream processes resulting in sexual dimorphism are produced by changes in gene regulation and result from a regulatory cascade and not from a vastly different gene composition. Thus, in this thesis a systems biology approach is used to understand sexual dimorphism at all molecular levels and how different genomic features, e.g. sex chromosome evolution, can affect the interplay of these molecules. The threespine stickleback, Gasterosteus aculeatus, is used as the model to investigate molecular mechanisms of sexual dimorphism. It has well-characterized ecology and behavior, especially in the breeding season when sexual dimorphism is high. Moreover, threespine stickleback has a recently evolved Y chromosome in the early stages of sex chromosome evolution, characterized by a lack of recombination leading to degeneration (i.e. gene loss). The aim of my thesis is to investigate how the genotype links to the molecular phenotype and relates to differences in molecular expression between males and females. Based on previous research on sex differences in mRNA expression, I investigated sex-biased protein expression in adult fish outside the breeding season to see if differences persisted after translation. As sex-biased expression also prevailed in the proteome and previous transcription expression seemed to be related to the sex chromosomes, I investigated the genome level with a particular focus on the sex-chromosomes. I characterized the status of Y chromosome degeneration in the threespine stickleback and its effects on gene function. Furthermore, since the degeneration process leaves genes in a single copy in males, I examined whether the resulting dosage difference of messenger RNA for hemizygous genes is compensated as it is in other organisms. In addition, threespine sticklebacks have wellcharacterized behavioral differences related to the male’s social status during the breeding season. To understand the connection between the genotype and behavior, I examined gene expression patterns related to breeding behavior using dominant and subordinate males as well as female

Natural born weight gainers: The mechanisms of obesity in transgenic mice overexpressing neuropeptide Y

Neuropeptide Y (NPY) is a neurotransmitter promoting energy storage by activating Y-receptors and thus affecting food intake, thermogenesis and adipose tissue metabolism. NPY is expressed both in the central and sympathetic nervous system. Hypothalamic NPY is known to stimulate feeding, but the effects of noradrenergic neuron NPY are more ambiguous. Chronic stress stimulates fat accumulation via NPY release from noradrenergic neurons. Furthermore, polymorphism in the human Npy gene has been associated with metabolic disturbances and increased NPY secretion after sympathetic stimulation. The main objective of this study was to clarify the mechanisms of noradrenergic neuron NPY in the development of obesity. The metabolic phenotype of a homozygous mouse overexpressing NPY in the brain noradrenergic neurons and sympathetic nervous system (OE-NPYDβH mouse) was characterized. OE-NPYDβH mice had an increased fat mass and body weight, which caused impairments of glucose metabolism and hyperinsulinaemia with age. There were no differences in energy intake or expenditure, but the sympathetic tone was down-regulated and the endocannabinoid system activated. Furthermore, peripheral Y2-receptors in energy-rich conditions played an important role in mediating the fat-accumulating effect of NPY. These results indicate that noradrenergic neuron NPY promotes obesity via direct effects in the periphery and by modulating the sympatho-adrenal and endocannabinoid systems. Additionally, NPY in the central noradrenergic neurons is believed to possess many important roles. The phenotype of the OE-NPYDβH mouse resembles the situations of chronic stress and Npy gene polymorphism and thus these mice may be exploited in testing novel drug candidates for the treatment of obesity.

Two-dimensional drift-diffusion simulation of organic solar cells

The aim of this master's thesis is to develop a two-dimensional drift-di usion model, which describes charge transport in organic solar cells. The main bene t of a two-dimensional model compared to a one-dimensional one is the inclusion of the nanoscale morphology of the active layer of a bulk heterojunction solar cell. The developed model was used to study recombination dynamics at the donor-acceptor interface. In some cases, it was possible to determine e ective parameters, which reproduce the results of the two-dimensional model in the one-dimensional case. A summary of the theory of charge transport in semiconductors was presented and discussed in the context of organic materials. Additionally, the normalization and discretization procedures required to nd a numerical solution to the charge transport problem were outlined. The charge transport problem was solved by implementing an iterative scheme called successive over-relaxation. The obtained solution is given as position-dependent electric potential, free charge carrier concentrations and current densities in the active layer. An interfacial layer, separating the pure phases, was introduced in order to describe charge dynamics occurring at the interface between the donor and acceptor. For simplicity, an e ective generation of free charge carriers in the interfacial layer was implemented. The pure phases simply act as transport layers for the photogenerated charges. Langevin recombination was assumed in the two-dimensional model and an analysis of the apparent recombination rate in the one-dimensional case is presented. The recombination rate in a two-dimensional model is seen to e ectively look like reduced Langevin recombination at open circuit. Replicating the J-U curves obtained in the two-dimensional model is, however, not possible by introducing a constant reduction factor in the Langevin recombination rate. The impact of an acceptor domain in the pure donor phase was investigated. Two cases were considered, one where the acceptor domain is isolated and another where it is connected to the bulk of the acceptor. A comparison to the case where no isolated domains exist was done in order to quantify the observed reduction in the photocurrent. The results show that all charges generated at the isolated domain are lost to recombination, but the domain does not have a major impact on charge transport. Trap-assisted recombination at interfacial trap states was investigated, as well as the surface dipole caused by the trapped charges. A theoretical expression for the ideality factor n_id as a function of generation was derived and shown to agree with simulation data. When the theoretical expression was fitted to simulation data, no interface dipole was observed.

Predicting diffusion of innovations with self-organisation and machine learning

The main subject of this master's thesis was predicting diffusion of innovations. The prediction was done in a special case: product has been available in some countries, and based on its diffusion in those countries the prediction is done for other countries. The prediction was based on finding similar countries with Self-Organizing Map~(SOM), using parameters of countries. Parameters included various economical and social key figures. SOM was optimised for different products using two different methods: (a) by adding diffusion information of products to the country parameters, and (b) by weighting the country parameters based on their importance for the diffusion of different products. A novel method using Differential Evolution (DE) was developed to solve the latter, highly non-linear optimisation problem. Results were fairly good. The prediction method seems to be on a solid theoretical foundation. The results based on country data were good. Instead, optimisation for different products did not generally offer clear benefit, but in some cases the improvement was clearly noticeable. The weights found for the parameters of the countries with the developed SOM optimisation method were interesting, and most of them could be explained by properties of the products.

The molecular mechanisms and evolutionary significance of plumage colour variation in pied flycatchers (Ficedula hypoleuca)

A common feature of natural populations is that individuals differ in morphology, physiologyand behavior (i.e .phenotype). A thorough understanding of the molecular mechanisms and evolutionary forces behind this phenotypic variation is a prerequisite for understanding evolution.This thesis examines the molecular mechanism and the roles of the different evolutionary forces in plumage colour variation in pied flycatchers (Ficedulahypoleuca). Malepied flycatchers exhibit marked variation in both pigmentary and structural plumage colourand the trait has repeatedly been suggested to be of adaptive significance. An examination of plumage colour variation on reproductive output trevealed that structural colouration, and more specifically the degree of ultraviolet (UV) reflectance had an effect on number of young sired. Paternity analyses of breeding males revealed that males that had been cuckolded by their social mate tended to be less UV reflectant than males that had not been cuckolded.Neither pigment-based norstructural colouration was found to affect the probability of siring young in other nests. Phenotypic differentiation was found to be markedly greater than differentiation at neutralgenetic markers across the pied flycatcher breeding range. Furthermore patterns of differentiationin phenotypes and selectively neutral genes were not uniform. Outlier tests searching for genomic footprints of selection revealed elevated levels of genetic divergence in a gene associated with feather development (and thus potentially structural colouration) and ultraviolet vision. Th eobserved differentiation in allelic frequencies was particularly pronounced in the Spanish piedflycatcher populations. Examining gene expression during feather development indicated that the TYRP1 gene (known to be involved in the production of black pigment) may be relevant in generating phenotypic variation in pied flycatcher plumage. Also, energy homeostasis related genesfeatured prominently among the genes found to be expressed in one extreme phenotype but not the other. This is of particular interest in light of what is known about the pleiotropy ofthe melanocortin system which underlies brown-black pigment production. The melanocortinsystem is also associated with energy homeostasis (among a number of other physiological functions) and thus the results could be pointing to the signalling function of brown-blackplumage. Plumage colour variation in pied flycatchers, both structural and pigmentary, can thus beconcluded to be exhibiting signals of non-neutral evolution. Structural colouration was found to play a role in sexual selection and putative signals of selection were further detected in acandidate gene for this trait. Evidence for non-neutral evolution of pigmentary colouration was also detected. These findings, together with the fact that preliminary evidence for an energy balance associated signalling function for plumage was found, present good starting points for further investigations into the meaning and mechanisms of plumage colour variation in piedflycatchers.