External signal-activated liposomal drug delivery systems

Modern drug discovery gives rise to a great number of potential new therapeutic agents, but in some cases the efficient treatment of patient may not be achieved because the delivery of active compounds to the target site is insufficient. Thus, drug delivery is one of the major challenges in current pharmaceutical research. Numerous nanoparticle-based drug carriers, e.g. liposomes, have been developed for enhanced drug delivery and targeting. Drug targeting may enhance the efficiency of the treatment and, importantly, reduce unwanted side effects by decreasing drug distribution to non-target tissues. Liposomes are biocompatible lipid-based carriers that have been studied for drug delivery during the last 40 years. They can be functionalized with targeting ligands and sensing materials for triggered activation. In this study, various external signal-assisted liposomal delivery systems were developed. Signals can be used to modulate drug permeation or release from the liposome formulation, and they provide accurate control of time, place and rate of activation. The study involved three types of signals that were used to trigger drug permeation and release: electricity, heat and light. Electrical stimulus was utilized to enhance the permeation of liposomal DNA across the skin. Liposome/DNA complex-mediated transfections were performed in tight rat epidermal cell model. Various transfection media and current intensities were tested, and transfection efficiency was evaluated non-invasively by monitoring the concentration of secreted reporter protein in cell culture medium. Liposome/DNA complexes produced gene expression, but electrical stimulus did not enhance the transfection efficiency significantly. Heat-sensitive liposomal drug delivery system was developed by coating liposomes with biodegradable and thermosensitive poly(N-(2-hydroxypropyl) methacrylamide-mono/dilactate polymer. Temperature-triggered liposome aggregation and contents release from liposomes were evaluated. The cloud point temperature (CP) of the polymer was set to 42 °C. Polymer-coated liposome aggregation and contents release were observed above CP of the polymer, while non-coated liposomes remained intact. Polymer precipitates above its CP and interacts with liposomal bilayers. It is likely that this induces permeabilization of the liposomal membrane and contents release. Light-sensitivity was introduced to liposomes by incorporation of small (< 5 nm) gold nanoparticles. Hydrophobic and hydrophilic gold nanoparticles were embedded in thermosensitive liposomes, and contents release was investigated upon UV light exposure. UV light-induced lipid phase transitions were examined with small angle X-ray scattering, and light-triggered contents release was shown also in human retinal pigment epithelial cell line. Gold nanoparticles absorb light energy and transfer it into heat, which induces phase transitions in liposomes and triggers the contents release. In conclusion, external signal-activated liposomes offer an advanced platform for numerous applications in drug delivery, particularly in the localized drug delivery. Drug release may be localized to the target site with triggering stimulus that results in better therapeutic response and less adverse effects. Triggering signal and mechanism of activation can be selected according to a specific application.

The Mechanisms of Hedgehog Signal Transduction

Studies in both vertebrates and invertebrates have identified proteins of the Hedgehog (Hh) family of secreted signaling molecules as key organizers of tissue patterning. Initially discovered in Drosophila in 1992, Hh family members have been discovered in animals with body plans as diverse as those of mammals, insects and echinoderms. In humans three related Hh genes have been identified: Sonic, Indian and Desert hedgehog (Shh, Ihh and Dhh). Transduction of the Hh signal to the cytoplasm utilizes an unusual mechanism involving consecutive repressive interactions between Hh and its receptor components, Patched (Ptc) and Smoothened (Smo). Several cytoplasmic proteins involved in Hh signal transduction are known in Drosophila, but mammalian homologs are known only for the Cubitus interruptus (Ci) transcription factor (GLI(1-3)) and for the Ci/GLI-associated protein, Suppressor of Fused (Su(fu)). In this study I analyzed the mechanisms of how the Hh receptor Ptc regulates the signal transducer Smo, and how Smo relays the Shh signal from the cell surface to the cytoplasm ultimately leading to the activation of GLI transcription factors. In Drosophila, the kinesin-like protein Costal2 (Cos2) is required for suppression of Hh target gene expression in the absence of ligand, and loss of Cos2 causes embryonic lethality. Cos2 acts by bridging Smo to the Ci. Another protein, Su(Fu) exerts a weak suppressive influence on Ci activity and loss of Su(Fu) causes subtle changes in Drosophila wing pattern. This study revealed that domains in Smo that are critical for Cos2 binding in Drosophila are dispensable for mammalian Smo function. Furthermore, by analyzing the function of Su(Fu) and the closest mouse homologs of Cos2 by protein overexpression and RNA interference I found that inhibition of the Hh response pathway in the absence of ligand does not require Cos2 activity, but instead critically depends on the activity of Su(Fu). These results indicate that a major change in the mechanism of action of a conserved signaling pathway occurred during evolution, probably through phenotypic drift made possible by the existence in some species of two parallel pathways acting between the Hh receptor and the Ci/GLI transcription factors. In a second approach to unravel Hh signaling we cloned > 90% of all human full-length protein kinase cDNAs and constructed the corresponding kinase-activity deficient mutants. Using this kinome resource as a screening tool, two kinases, MAP3K10 and DYRK2 were found to regulate Shh signaling. DYRK2 directly phosphorylated and induced the proteasome dependent degradation of the key Hh-pathway regulated transcription factor, GLI2. MAP3K10, in turn, affected GLI2 indirectly by modulating the activity of DYRK2.

Methods to improve gene signal: Application to cDNA microarrays

Microarrays are high throughput biological assays that allow the screening of thousands of genes for their expression. The main idea behind microarrays is to compute for each gene a unique signal that is directly proportional to the quantity of mRNA that was hybridized on the chip. A large number of steps and errors associated with each step make the generated expression signal noisy. As a result, microarray data need to be carefully pre-processed before their analysis can be assumed to lead to reliable and biologically relevant conclusions. This thesis focuses on developing methods for improving gene signal and further utilizing this improved signal for higher level analysis. To achieve this, first, approaches for designing microarray experiments using various optimality criteria, considering both biological and technical replicates, are described. A carefully designed experiment leads to signal with low noise, as the effect of unwanted variations is minimized and the precision of the estimates of the parameters of interest are maximized. Second, a system for improving the gene signal by using three scans at varying scanner sensitivities is developed. A novel Bayesian latent intensity model is then applied on these three sets of expression values, corresponding to the three scans, to estimate the suitably calibrated true signal of genes. Third, a novel image segmentation approach that segregates the fluorescent signal from the undesired noise is developed using an additional dye, SYBR green RNA II. This technique helped in identifying signal only with respect to the hybridized DNA, and signal corresponding to dust, scratch, spilling of dye, and other noises, are avoided. Fourth, an integrated statistical model is developed, where signal correction, systematic array effects, dye effects, and differential expression, are modelled jointly as opposed to a sequential application of several methods of analysis. The methods described in here have been tested only for cDNA microarrays, but can also, with some modifications, be applied to other high-throughput technologies. Keywords: High-throughput technology, microarray, cDNA, multiple scans, Bayesian hierarchical models, image analysis, experimental design, MCMC, WinBUGS.

Food intake, growth and social interactions of signal crayfish, Pacifastacus leniusculus (Dana)

The knowledge about the optimal rearing conditions, such as water temperature and quality, photoperiod and density, with the understanding of animal nutritional requirements forms the basis of economically stable aquaculture for freshwater crayfish. However, the shift from a natural environment to effective culture conditions induces several changes, not only at the population level, but also at the individual level. The social contacts between conspecifics increase with increasing animal density. The competition for limited resources (e.g. food, shelter, mates) is more severe with the presence of agonistic behaviour and may lead to unequal distribution of these. The objectives of this study were to: 1) study the distribution of a common food resource between communally reared signal crayfish (Pacifastacus leniusculus) and to assign potential feeding hierarchy on the basis of individual food intake measurements, 2) explore the possibilities of size distribution manipulations to affect population dynamics and food intake to improve growth and survival in culture and 3) study the effect of food ration and spatial distribution on food intake and to explore the effect of temperature and food ration on growth and body composition of freshwater crayfish. The feeding ranks between animals were assigned with a new method for individual food intake measurement of communally reared crayfish. This technique has a high feasibility and a great potential to be applied in crayfish aquaculture studies. In this study, signal crayfish showed high size-related variability in food consumption both among individuals within a group (inter-individual) and within individual day-to-day variation (intra-individual). Increased competition for food led to an unequal distribution of this resource and this may be a reason for large growth differences between animals. The consumption was significantly higher when reared individually in comparison with communal housing. These results suggest that communally housed crayfish form a feeding hierarchy and that the animal size is the major factor controlling the position in this hierarchy. The optimisation of the social environment ( social conditions ) was evaluated in this study as a new approach to crayfish aquaculture. The results showed that the absence of conspecifics (individual rearing vs. communal housing) affects growth rate, food intake and the proportion of injured animals, whereas size variation between animals influences the number and duration of agonistic encounters. In addition, animal size had a strong influence on the fighting success of signal crayfish reared in a social milieu with a wide size variation of conspecifics. Larger individuals initiated and won most of the competitions, which suggests size-based social hierarchy of P. leniusculus. This is further supported by the fact that the length and weight gain of smaller animals increased after size grading, maybe because of a better access to the food resource due to diminished social pressure. However, the high dominance index was not based on size under conditions of limited size variation, e.g. those characteristic of restocked natural populations and aquaculture, indicating the important role of behaviour on social hierarchy.

Pathogen-Induced Defense Signaling and Signal Crosstalk in Arabidopsis

Erwinia carotovora subsp. carotovora is a bacterial phytopathogen that causes soft rot in various agronomically important crop plants. A genetically specified resistance to E. carotovora has not been defined, and plant resistance to this pathogen is established through nonspecific activation of basal defense responses. This, together with the broad host range, makes this pathogen a good model for studying the activation of plant defenses. Production and secretion of plant cell wall-degrading enzymes (PCWDE) are central to the virulence of E. carotovora. It also possesses the type III secretion system (TTSS) utilized by many Gram-negative bacteria to secrete virulence- promoting effector proteins to plant cells. This study elucidated the role of E. carotovora HrpN (HrpNEcc), an effector protein secreted through TTSS, and the contribution of this protein in the virulence of E. carotovora. Treatment of plants with HrpNEcc was demonstrated to induce a hypersensitive response (HR) as well as resistance to E. carotovora. Resistance induced by HrpNEcc required both salicylic acid (SA)- and jasmonate/ethylene (JA/ET)-dependent defense signaling in Arabidopsis. Simultaneous treatment of Arabidopsis with HrpNEcc and PCWDE polygalacturonase PehA elicited accelerated and enhanced induction of defense genes but also increased production of superoxide and lesion formation. This demonstrates mutual amplification of defense signaling by these two virulence factors of E. carotovora. Identification of genes that are rapidly induced in response to a pathogen can provide novel information about the early events occurring in the plant defense response. CHLOROPHYLLASE 1 (AtCLH1) and EARLY RESPONSIVE TO DEHYDRATION 15 (ERD15) are both rapidly triggered by E. carotovora in Arabidopsis. Characterization of AtCLH1 encoding chlorophyll-degrading enzyme chlorophyllase indicated that it might have a role in chlorophyll degradation during plant tissue damage. Silencing of this gene resulted in increased accumulation of reactive oxygen species (ROS) in response to pathogen infection in a light-dependent manner. This led to enhanced SA-dependent defenses and resistance to E. carotovora. Moreover, crosstalk between different defense signaling pathways was observed; JA-dependent defenses and resistance to fungal pathogen Alternaria brassicicola were impaired, indicating antagonism between SA- and JA-dependent signaling. Characterization of ERD15 suggested that it is a novel, negative regulator of abscisic acid (ABA) signaling in Arabidopsis. Overexpression of ERD15 resulted in insensitivity to ABA and reduced tolerance of the plants to dehydration stress. However, simultaneously, the resistance of the plants to E. carotovora was enhanced. Silencing of ERD15 improved freezing and drought tolerance of transgenic plants. This, together with the reducing effect of ABA on seed germination, indicated hypersensitivity to this phytohormone. ERD15 was hypothesized to act as a capacitor that controls the appropriate activation of ABA responses in Arabidopsis.

Phylogenetic studies of cyanobacterial lichens

Phylogenetic studies of cyanobacterial lichens Lichens are symbiotic assemblages between fungi (mycobiont) and green algae (phycobiont) or/and cyanobacteria (cyanobiont). Fossil records show that lichen-like symbioses occurred already 600 million years ago. Lichen symbiosis has since then become an important life strategy for the Fungi, particularly for species in the phylum Ascomycota as approximately 98% of the lichenized fungal species are ascomycetes. The taxonomy of lichen associations is based on the mycobiont. We reconstructed, using DNA sequence data, hypotheses of phylogenetic relationships of lichen-forming fungi that include species associated with cyanobacteria. These hypotheses of phylogeny should form the basis for the taxonomy. They also allowed studies of the origin and the evolution of specific symbioses. Genetic diversity and phylogenetic relationships of symbiotic cyanobionts were also studied in order to examine selectivity of cyanobionts and mycobionts as well as possible co-evolution between partners involved in lichen associations. The suggested circumscription of the family Stereocaulaceae to include Stereocaulon and Lepraria is supported. The recently described crustose Stereocaulon species seem to be correctly placed in the genus, although Stereocaulon traditionally included only fruticose species. The monospecific crustose genus Muhria is also shown to be best placed in Stereocaulon. Family Lobariaceae as currently delimited is monophyletic. Within Lobariaceae genus Sticta including Dendriscocaulon dendroides form a monophyletic group while the genera Lobaria and Pseudocyphellaria are non-monophyletic. A new classification of Lobariaceae is obviously needed. Further studies are however required before a final proposal for a new classification can be made. Our results show that the cyanobacterial symbiotic state has been gained repeatedly in the Ascomycota while losses of symbiotic cyanobacteria appear to be rare. The symbiosis with green algae is confirmed to have been gained repeatedly in Ascomycota but also repeatedly lost. Cyanobacterial symbioses therefore seem to be more stable than green algal associations. Cyanobacteria are perhaps more beneficial for the lichen fungi and therefore maintained. The results indicate a dynamic association of the lichen symbiosis. This evolutionary instability will perhaps be important for the lichen fungi as the utilization of options will perhaps enable lichens to colonize new substrates and survive environmental changes. Some cyanobacterial lichen genera seem to be highly selective towards the cyanobiont while others form symbioses with a broad spectrum of cyanobacteria. No evidence of co-evolution between fungi and cyanobacteria in cyanolichens could be demonstrated.

Spectral Tuning and Adaptation to Different Light Environments of Mysid Visual Pigments

In the present thesis, questions of spectral tuning, the relation of spectral and thermal properties of visual pigments, and evolutionary adaptation to different light environments were addressed using a group of small crustaceans of the genus Mysis as a model. The study was based on microspectrophotometric measurements of visual pigment absorbance spectra, electrophysiological measurements of spectral sensitivities of dark-adapted eyes, and sequencing of the opsin gene retrieved through PCR. The spectral properties were related to the spectral transmission of the respective light environments, as well as to the phylogentic histories of the species. The photoactivation energy (Ea) was estimated from temperature effects on spectral sensitivity in the long-wavelength range, and calculations were made for optimal quantum catch and optimal signal-to-noise ratio in the different light environments. The opsin amino acid sequences of spectrally characterized individuals were compared to find candidate residues for spectral tuning. The general purpose was to clarify to what extent and on what time scale adaptive evolution has driven the functional properties of (mysid) visual pigments towards optimal performance in different light environments. An ultimate goal was to find the molecular mechanisms underlying the spectral tuning and to understand the balance between evolutionary adaptation and molecular constraints. The totally consistent segregation of absorption maxima (λmax) into (shorter-wavelength) marine and (longer-wavelength) freshwater populations suggests that truly adaptive evolution is involved in tuning the visual pigment for optimal performance, driven by selection for high absolute visual sensitivity. On the other hand, the similarity in λmax and opsin sequence between several populations of freshwater M. relicta in spectrally different lakes highlights the limits to adaptation set by evolutionary history and time. A strong inverse correlation between Ea and λmax was found among all visual pigments studied in these respects, including those of M. relicta and 10 species of vertebrate pigments, and this was used to infer thermal noise. The conceptual signal-to-noise ratios thus calculated for pigments with different λmax in the Baltic Sea and Lake Pääjärvi light environments supported the notion that spectral adaptation works towards maximizing the signal-to-noise ratio rather than quantum catch as such. Judged by the shape of absorbance spectra, the visual pigments of all populations of M. relicta and M. salemaai used exclusively the A2 chromophore (3, 4-dehydroretinal). A comparison of amino acid substitutions between M. relicta and M. salemaai indicated that mysid shrimps have a small number of readily available tuning sites to shift between a shorter - and a longer -wavelength opsin. However, phylogenetic history seems to have prevented marine M. relicta from converting back to the (presumably) ancestral opsin form, and thus the more recent reinvention of marine spectral sensitivity has been accomplished by some other novel mechanism, yet to be found

Phylogenetics, taxon delimitation and phylogenetic diversity

Throughout the history of Linnean taxonomy, species have been described with varying degrees of justification. Many descriptions have been based on only a few ambiguous morphological characters. Moreover, species have been considered natural, well-defined units whereas higher taxa have been treated as disparate, non-existent creations. In the present thesis a few such cases were studied in detail. Often the species-level descriptions were based on only a few specimens and the variation previously thought to be interspecific was found to be intraspecific. In some cases morphological characters were sufficient to resolve the evolutionary relationships between the taxa, but generally more resolution was gained by the addition of molecular evidence. However, both morphological and molecular data were found to be deceptive in some cases. The DNA sequences of morphologically similar specimens were found to differ distinctly in some cases, whereas in other closely related species the morphology of specimens with identical DNA sequences differed substantially. This study counsels caution when evolutionary relationships are being studied utilizing only one source of evidence or a very limited number of characters (e.g. barcoding). Moreover, it emphasizes the importance of high quality data as well as the utilization of proper methods when making scientific inferences. Properly conducted analyses produce robust results that can be utilized in numerous interesting ways. The present thesis considered two such extensions of systematics. A novel hypothesis on the origin of bioluminescence in Elateriformia beetles is presented, tying it to the development of the clicking mechanism in the ancestors of these animals. An entirely different type of extension of systematics is the proposed high value of the white sand forests in maintaining the diversity of beetles in the Peruvian Amazon. White sand forests are under growing pressure from human activities that lead to deforestation. They were found to harbor an extremely diverse beetle fauna and many taxa were specialists living only in this unique habitat. In comparison to the predominant clay soil forests, considerably more elateroid beetles belonging to all studied taxonomic levels (species, genus, tribus, and subfamily) were collected in white sand forests. This evolutionary diversity is hypothesized to be due to a combination of factors: (1) the forest structure, which favors the fungus-plant interactions important for the elateroid beetles, (2) the old age of the forest type favoring survival of many evolutionary lineages and (3) the widespread distribution and fragmentation of the forests in the Miocene, favoring speciation.

Atrial Electric Signal During Sinus Rhythm in Lone Paroxysmal Atrial Fibrillation

Atrial fibrillation (AF) is the most common tachyarrhythmia and is associated with substantial morbidity, increased mortality and cost. The treatment modalities of AF have increased, but results are still far from optimal. More individualized therapy may be beneficial. Aiming for this calls improved diagnostics. Aim of this study was to find non-invasive parameters obtained during sinus rhythm reflecting electrophysiological patterns related to propensity to AF and particularly to AF occurring without any associated heart disease, lone AF. Overall 240 subjects were enrolled, 136 patients with paroxysmal lone AF and 104 controls (mean age 45 years, 75% males). Signal measurements were performed by non-invasive magnetocardiography (MCG) and by invasive electroanatomic mapping (EAM). High-pass filtering techniques and a new method based on a surface gradient technique were adapted to analyze atrial MCG signal. The EAM was used to elucidate atrial activation in patients and as a reference for MCG. The results showed that MCG mapping is an accurate method to detect atrial electrophysiologic properties. In lone paroxysmal AF, duration of the atrial depolarization complex was marginally prolonged. The difference was more obvious in women and was also related to interatrial conduction patterns. In the focal type of AF (75%), the root mean square (RMS) amplitudes of the atrial signal were normal, but in AF without demonstrable triggers the late atrial RMS amplitudes were reduced. In addition, the atrial characteristics tended to remain similar even when examined several years after the first AF episodes. The intra-atrial recordings confirmed the occurrence of three distinct sites of electrical connection from right to left atrium (LA): the Bachmann bundle (BB), the margin of the fossa ovalis (FO), and the coronary sinus ostial area (CS). The propagation of atrial signal could also be evaluated non-invasively. Three MCG atrial wave types were identified, each of which represented a distinct interatrial activation pattern. In conclusion, in paroxysmal lone AF, active focal triggers are common, atrial depolarization is slightly prolonged, but with a normal amplitude, and the arrhythmia does not necessarily lead to electrical or mechanical dysfunction of the atria. In women the prolongation of atrial depolarization is more obvious. This may be related to gender differences in presentation of AF. A significant minority of patients with lone AF lack frequent focal triggers, and in them, the late atrial signal amplitude is reduced, possibly signifying a wider degenerative process in the LA. In lone AF, natural impulse propagation to LA during sinus rhythm goes through one or more of the principal pathways described. The BB is the most common route, but in one-third, the earliest LA activation occurs outside the BB. Susceptibility to paroxysmal lone AF is associated with propagation of the atrial signal via the margin of the FO or via multiple pathways. When conduction occurs via the BB, it is related with prolonged atrial activation. Thus, altered and alternative conduction pathways may contribute to pathogenesis of lone AF. There is growing evidence of variability in genesis of AF also within lone paroxysmal AF. Present study suggests that this variation may be reflected in cardiac signal pattern. Recognizing the distinct signal profiles may assist in understanding the pathogenesis of AF and identifying subgroups for patient-tailored therapy.

Phylogenetic analyses reveal high levels of polyphyly among pleurocarpous lineages as well as novel clades

Phylogenetic analyses of the Hypnales usually show the same picture of poorly resolved trees with a large number of polyphyletic taxa and low support for the few reconstructed clades. One odd clade, however, consisting of three genera that are currently treated either within the Leskeaceae (Miyabea) or Neckeraceae (Homaliadelphus and Bissetia), was retrieved in a previously published phylogeny based on chloroplast rbcL. In order to elucidate the reliability of the observed Homaliadelphus - Miyabea - Bissetia - clade (HMB-clade) and to reveal its phylogenetic relationships a molecular study based on a representative set of hypnalean taxa was performed. Sequence data from all three genomes, namely the ITS1 and 2 (nuclear), the trnS-rps4-trnT-trnL-trnF cluster (plastid), the nad5 intron (mitochondrial), were analyzed. Although the phylogenetic reconstruction of the combined data set was not fully resolved regarding the backbone it clearly indicated the polyphyletic nature of various hypnalean families, such as the Leskeaceae, Hypnaceae, Hylocomiaceae, Neckeraceae, Leptodontaceae and Anomodontaceae with respect to the included taxa. In addition the results favor the inclusion of the Leptodontaceae and Thamnobryaceae in the Neckeraceae. The maximally supported HMB-clade consisting of the three genera Homaliadelphus (2-3 species), Miyabea (3 species) and Bissetia (1 species) is resolved sister to a so far unnamed clade comprising Taxiphyllum aomoriense, Glossadelphus ogatae and Leptopterigynandrum. The well-resolved and supported HMB-clade, here formally described as the Miyabeaceae, fam. nov. is additionally supported by morphological characters such as strongly incrassate, porose leaf cells, a relatively weak and diffuse costa and the presence of dwarf males. The latter are absent in the Neckeraceae and the Leskeaceae. It is essentially an East Asian family, with one species occurring in North America.

The many ways to delimit species: hairs, genes and surface chemistry

Species identification forms the basis for understanding the diversity of the living world, but it is also a prerequisite for understanding many evolutionary patterns and processes. The most promising approach for correctly delimiting and identifying species is to integrate many types of information in the same study. Our aim was to test how cuticular hydro- carbons, traditional morphometrics, genetic polymorphisms in nuclear markers (allozymes and DNA microsatellites) and DNA barcoding (partial mitochondrial COI gene) perform in delimiting species. As an example, we used two closely related Formica ants, F. fusca and F. lemani, sampled from a sympatric population in the northern part of their distribu- tion. Morphological characters vary and overlap in different parts of their distribution areas, but cuticular hydrocarbons include a strong taxonomic signal and our aim is to test the degree to which morphological and genetic data correspond to the chemical data. In the morphological analysis, species were best separated by the combined number of hairs on pro- notum and mesonotum, but individual workers overlapped in hair numbers, as previously noted by several authors. Nests of the two species were separated but not clustered according to species in a Principal Component Analysis made on nuclear genetic data. However, model-based Bayesian clustering resulted in perfect separation of the species and gave no indication of hybridization. Furthermore, F. lemani and F. fusca did not share any mitochondrial haplotypes, and the species were perfectly separated in a phylogenetic tree. We conclude that F. fusca and F. lemani are valid species that can be separated in our study area relatively well with all methods employed. However, the unusually small genetic differen- tiation in nuclear markers (FST = 0.12) shows that they are closely related, and occasional hybridization between F. fusca and F. lemani cannot be ruled out.