Fishery-induced selection on an Alpine whitefish: quantifying genetic and environmental effects on individual growth rate

Size-selective fishing, environmental changes and reproductive strategies are expected to affect life-history traits such as the individual growth rate. The relative contribution of these factors is not clear, particularly whether size-selective fishing can have a substantial impact on the genetics and hence on the evolution of individual growth rates in wild populations. We analysed a 25-year monitoring survey of an isolated population of the Alpine whitefish Coregonus palaea. We determined the selection differentials on growth rate, the actual change of growth rate over time and indicators of reproductive strategies that may potentially change over time. The selection differential can be reliably estimated in our study population because almost all the fish are harvested within their first years of life, i.e. few fish escape fishing mortality. We found a marked decline in average adult growth rate over the 25 years and a significant selection differential for adult growth, but no evidence for any linear change in reproductive strategies over time. Assuming that the heritability of growth in this whitefish corresponds to what was found in other salmonids, about a third of the observed decline in growth rate would be linked to fishery-induced evolution. Size-selective fishing seems to affect substantially the genetics of individual growth in our study population.

An integrated genomic and expression analysis of 7q deletion in splenic marginal zone lymphoma.

Splenic marginal zone lymphoma (SMZL) is an indolent B-cell lymphoproliferative disorder characterised by 7q32 deletion, but the target genes of this deletion remain unknown. In order to elucidate the genetic target of this deletion, we performed an integrative analysis of the genetic, epigenetic, transcriptomic and miRNomic data. High resolution array comparative genomic hybridization of 56 cases of SMZL delineated a minimally deleted region (2.8 Mb) at 7q32, but showed no evidence of any cryptic homozygous deletion or recurrent breakpoint in this region. Integrated transcriptomic analysis confirmed significant under-expression of a number of genes in this region in cases of SMZL with deletion, several of which showed hypermethylation. In addition, a cluster of 8 miRNA in this region showed under-expression in cases with the deletion, and three (miR-182/96/183) were also significantly under-expressed (P<0.05) in SMZL relative to other lymphomas. Genomic sequencing of these miRNA and IRF5, a strong candidate gene, did not show any evidence of somatic mutation in SMZL. These observations provide valuable guidance for further characterisation of 7q deletion.

Impaired glutathione synthesis in schizophrenia: convergent genetic and functional evidence

Schizophrenia is a complex multifactorial brain disorder with a genetic component. Convergent evidence has implicated oxidative stress and glutathione (GSH) deficits in the pathogenesis of this disease. The aim of the present study was to test whether schizophrenia is associated with a deficit of GSH synthesis. Cultured skin fibroblasts from schizophrenia patients and control subjects were challenged with oxidative stress, and parameters of the rate-limiting enzyme for the GSH synthesis, the glutamate cysteine ligase (GCL), were measured. Stressed cells of patients had a 26% (P = 0.002) decreased GCL activity as compared with controls. This reduction correlated with a 29% (P < 0.001) decreased protein expression of the catalytic GCL subunit (GCLC). Genetic analysis of a trinucleotide repeat (TNR) polymorphism in the GCLC gene showed a significant association with schizophrenia in two independent case-control studies. The most common TNR genotype 7/7 was more frequent in controls [odds ratio (OR) = 0.6, P = 0.003], whereas the rarest TNR genotype 8/8 was three times more frequent in patients (OR = 3.0, P = 0.007). Moreover, subjects with disease-associated genotypes had lower GCLC protein expression (P = 0.017), GCL activity (P = 0.037), and GSH contents (P = 0.004) than subjects with genotypes that were more frequent in controls. Taken together, the study provides genetic and functional evidence that an impaired capacity to synthesize GSH under conditions of oxidative stress is a vulnerability factor for schizophrenia.

Molecular and cellular analysis of genodermatoses

Summary Skin is the essential interface between our body and its environment; not only does it prevent water loss and protect us from external insults it also plays an essential role in the central nervous system acting as a major sense organ primarily for touch and pain. The main cell type present in skin, keratinocyte, undergoes a differentiation process leading to the formation of this protecting barrier. This work is intended to contribute to the understanding of how keratinocyte differentiates and skin functions. To do this, we studied two genetic skin diseases: Erythrokeratodermia variabilis and Mal de Meleda. Our approach was to examine the expression and localization of proteins implicated in these two pathologies in normal and diseased tissues and to determine the influence of mutant proteins at the molecular and cellular levels. Connexins are major components of gap junctions, channels allowing direct communication between cells. Our laboratory has identified mutations in both connexin 30.3 (Cx30.3) and 31 (Cx31) to be causally involved in erythrokeratodermia variabilis (EKV), an autosomal dominant disorder of keratinization. In the first chapter, we show a new mutation of Cx31, L209P-Cx31, in 3 EKV patients, extending the field of EKV-causing mutations although the mechanism by which connexin mutations lead to the disease is unclear. In the second chapter, we studied the effect of F137L-Cx30.3 on expression, trafficking and localization of cotransfected Cx31 and Cx30.3 in connexin-deficient HeLa cells. The F137 amino acid, highly conserved in connexin family, is oriented towards the channel pore and F137L mutation in either Cx30.3 or Cx31 lead to EKV. As two genes can lead to EKV when mutated, our hypothesis was that Cx31 and Cx30.3 might cooperate at a molecular level. We were able to demonstrate a physical interaction between Cx31 and Cx30.3. The presence of F137L-Cx30.3 disturbed the trafficking of both connexins, less connexins were integrated into gap junctions and thus, the coupling between cell was diminished. Connexins formed in the presence of F137L-Cx30.3 are degraded at their exit from the endoplasmic reticulum. In conclusion, our results indicate that the genetic heterogeneity of EKV is due to mutations in two interacting proteins. F137L-Cx30.3 has a dominant negative effect and affects Cx31, disturbing cellular communication in epidermal cells. Mal de Meleda is an autosomal recessive inflammatory and a keratotic palmoplantar skin disorder due to mutations in SLURP1 (secreted LY6/PLAUR-related protein 1). SLURP1 belongs to the LY6/PLAUR family of proteins and has the particularity of being secreted instead of being GPI-anchored. The high degree of structural similarity between SLURP1 and the three fingers motif of snake neurotoxins and LYNX 1-C suggests that this protein could interact with the neuronal acetylcholine receptors. In the third chapter, we show that SLURP1 potentiates responses of the a7 nicotinic acetylcholine receptor (nAchR) to acetylcholine. These results identify SLURP1 as a secreted epidermal neuromodulator that is likely to be essential for palmoplantar skin. In the fourth chapter, we show that SLURP1 is expressed in the granular layer of the epidermis but is absent from skin biopsies of Mal de Meleda patients. SLURP1 is also present in secretions such as sweat, tears or saliva. An in vitro analysis on two mutant of SLURP-I demonstrates that W15R-SLURP1 is absent in cells while G86R-SLURP1 is expressed and secreted, suggesting that SLURP1 can lead to the disease by either an absent or an abnormal protein. Finally, in the fifth chapter, we analyse the expression and biological properties of other LY6/PLAUR members, clustered around SLURP] on chromosome 8. Their GPI-anchored or secreted status were analysed in vitro. SLURP1, LYNX1-A and -B are secreted while LYPDC2 and LYNX 1-C are GPI anchored. Three of these proteins are expressed in the epidermis and in cultured keratinocytes. These results suggest that these LY6/PLAUR members may have an important role in skin homeostasis. Résumé Résumé La peau est la barrière essentielle entre notre corps et l'environnement, nous protégeant des agressions extérieures, de la déshydratation et assurant aussi un rôle dans le système nerveux central en tant qu'organe du toucher et de la douleur. Le principal type de cellules présent dans la peau est le kératinocyte qui suit un processus de différenciation aboutissant à la formation de cette barrière protectrice. Ce travail est destiné à comprendre la différenciation des kératinocytes et le fonctionnement de la peau. Pour cela, nous avons étudié deux maladies génodermatoses : l'Erthrokeratodermia Variabilis (EKV) et le Mal de Meleda. Nous avons examiné l'expression et la localisation des protéines impliquées dans ces deux pathologies dans des tissus normaux et malades puis déterminé l'influence des protéines mutantes aux niveaux moléculaires et cellulaires. Les connexines (Cx) sont les composants majeurs des jonctions communicantes, canaux permettant la communication directe entre les cellules. Notre laboratoire a identifié des mutations dans les Cx30.3 et Cx31 comme responsables de l'EKV, génodermatose de transmission autosomique dominante. Dans le ler chapitre, nous décrivons une nouvelle mutation de Cx31, L209-Cx31, et contribuons à l'établissement du catalogue des mutations de Cx31 entraînant cette maladie. Cependant, le mécanisme par lequel les mutations de Cx31 et C3x0.3 provoquent l'EKV est inconnu. Dans le 2ème chapitre, nous étudions les effets de la mutation F137L-Cx30.3 sur l'expression, le trafic et la localisation des Cx31 et Cx30.3 transfectées dans des cellules HeLa, déficientes en connexines. Comme deux gènes peuvent causer une EKV quand ils sont mutés, notre hypothèse était que Cx31 et Cx30.3 pourraient coopérer au niveau moléculaire. Nous avons montré l'existence d'une interaction physique entre ces deux connexines. La présence de la mutation F137L-Cx30.3 perturbe le trafic des deux connexines, moins de connexines sont intégrées dans les jonctions communicantes et donc le couplage entre les cellules est diminué. Les connexons formés en présence de cette mutation sont dégradés à leur sortie du réticulum endoplasmique. En conclusion, nos résultats indiquent que l'hétérogénéité génétique de EKV est due à des mutations dans deux protéines qui interagissent. F137L-Cx30.3 a un effet dominant négatif et affecte Cx31, perturbant la communication entre les cellules épidermiques. Le Mal de Meleda est une maladie récessive de la peau palmoplantaire due à des mutations dans SLURP1. SLURP1 appartient à la famille des protéines contenant un domaine LY6/PLAUR et a la particularité d'être sécrétée. La grande homologie de structure existant entre SLURP1, les neurotoxines de serpent et LYNX1-C suggère que la protéine pourrait interagir avec des récepteurs à acétylcholine (Ach). Dans le 3ème chapitre, nous montrons que SLURP1 module la réponse à l'Ach du récepteur nicotinique α7. Ces résultats identifient SLURP1 comme un neuromodulateur épidermique sécrété, probablement essentiel pour la peau palmoplantaire. Dans le 4ème chapitre, nous montrons que SLURP1 est exprimé dans la couche granuleuse de l'épiderme et qu'il est absent des biopsies des patients. SLURP1 a aussi été détecté dans des sécrétions telles que la sueur, les lamies et la salive. Une analyse in vitro de deux mutants de SLURP1 a montré que W15R-SLURP1 est absent des cellules tandis que G86R-SLURP1 est exprimé et sécrété, suggérant qu'une absence ou une anomalie de SLURP1 peuvent causer la maladie. Finalement, dans le 5ème chapitre, nous analysons l'expression et les propriétés biologiques d'autres membres de la famille LY6/PLAUR localisés autour de SLURP1 sur le chromosome 8. Leur statut de protéines sécrétées ou liées à la membrane par une ancre GPI est analysé in vitro. SLURP1, LYNXI-A et -B sont sécrétées alors que LYPDC2 et LYNX1-C sont liés à la membrane. Trois de ces protéines sont exprimées dans l'épiderme et dans des kératinocytes cultivés. Ces résultats suggèrent que la famille LY6/PLAUR pourrait avoir un rôle important dans l'homéostasie de la peau.

Association between variants of the leptin receptor gene (LEPR) and overweight: a systematic review and an analysis of the CoLaus study.

BACKGROUND: Three non-synonymous single nucleotide polymorphisms (Q223R, K109R and K656N) of the leptin receptor gene (LEPR) have been tested for association with obesity-related outcomes in multiple studies, showing inconclusive results. We performed a systematic review and meta-analysis on the association of the three LEPR variants with BMI. In addition, we analysed 15 SNPs within the LEPR gene in the CoLaus study, assessing the interaction of the variants with sex. METHODOLOGY/PRINCIPAL FINDINGS: We searched electronic databases, including population-based studies that investigated the association between LEPR variants Q223R, K109R and K656N and obesity- related phenotypes in healthy, unrelated subjects. We furthermore performed meta-analyses of the genotype and allele frequencies in case-control studies. Results were stratified by SNP and by potential effect modifiers. CoLaus data were analysed by logistic and linear regressions and tested for interaction with sex. The meta-analysis of published data did not show an overall association between any of the tested LEPR variants and overweight. However, the choice of a BMI cut-off value to distinguish cases from controls was crucial to explain heterogeneity in Q223R. Differences in allele frequencies across ethnic groups are compatible with natural selection of derived alleles in Q223R and K109R and of the ancient allele in K656N in Asians. In CoLaus, the rs10128072, rs3790438 and rs3790437 variants showed interaction with sex for their association with overweight, waist circumference and fat mass in linear regressions. CONCLUSIONS: Our systematic review and analysis of primary data from the CoLaus study did not show an overall association between LEPR SNPs and overweight. Most studies were underpowered to detect small effect sizes. A potential effect modification by sex, population stratification, as well as the role of natural selection should be addressed in future genetic association studies.

Genetic and environmental effects on the covariation betweencolour polymorphism and a life history trait

Variation in coloration with a strong underlying genetic basis is frequently found within animal populations but little is known about its function. Covariation between colour polymorphism and life-history traits can arise because morphs perform differently among environments or because they possess alternative alleles coding for key life-history traits. To test these two hypotheses, we studied a population of tawny owls Strix aluco, a bird displaying red, brown and grey morphs. We assessed the colour morph of breeding females, swapped eggs or hatchlings between pairs of nests, and examined how body condition in 3-week-old nestlings covaries with coloration of foster and genetic mothers. Redder foster and genetic mothers produced young in better condition. Because in two other years we observed that greyish females produced offspring in better condition than those of red females, the present study suggests that colour polymorphism signals genetic and phenotypic adaptations to cope with a fluctuating environment.

Early onset collagen VI myopathies: Genetic and clinical correlations.

OBJECTIVE: Mutations in the genes encoding the extracellular matrix protein collagen VI (ColVI) cause a spectrum of disorders with variable inheritance including Ullrich congenital muscular dystrophy, Bethlem myopathy, and intermediate phenotypes. We extensively characterized, at the clinical, cellular, and molecular levels, 49 patients with onset in the first 2 years of life to investigate genotype-phenotype correlations. METHODS: Patients were classified into 3 groups: early-severe (18%), moderate-progressive (53%), and mild (29%). ColVI secretion was analyzed in patient-derived skin fibroblasts. Chain-specific transcript levels were quantified by quantitative reverse transcriptase polymerase chain reaction (qRT-PCR), and mutation identification was performed by sequencing of complementary DNA. RESULTS: ColVI secretion was altered in all fibroblast cultures studied. We identified 56 mutations, mostly novel and private. Dominant de novo mutations were detected in 61% of the cases. Importantly, mutations causing premature termination codons (PTCs) or in-frame insertions strikingly destabilized the corresponding transcripts. Homozygous PTC-causing mutations in the triple helix domains led to the most severe phenotypes (ambulation never achieved), whereas dominant de novo in-frame exon skipping and glycine missense mutations were identified in patients of the moderate-progressive group (loss of ambulation). INTERPRETATION: This work emphasizes that the diagnosis of early onset ColVI myopathies is arduous and time-consuming, and demonstrates that quantitative RT-PCR is a helpful tool for the identification of some mutation-bearing genes. Moreover, the clinical classification proposed allowed genotype-phenotype relationships to be explored, and may be useful in the design of future clinical trials.

RNA-based gene duplication: mechanistic and evolutionary insights.

Gene copies that stem from the mRNAs of parental source genes have long been viewed as evolutionary dead-ends with little biological relevance. Here we review a range of recent studies that have unveiled a significant number of functional retroposed gene copies in both mammalian and some non-mammalian genomes. These studies have not only revealed previously unknown mechanisms for the emergence of new genes and their functions but have also provided fascinating general insights into molecular and evolutionary processes that have shaped genomes. For example, analyses of chromosomal gene movement patterns via RNA-based gene duplication have shed fresh light on the evolutionary origin and biology of our sex chromosomes.

A comparative phenotypic and genomic analysis of C57BL/6J and C57BL/6N mouse strains.

BACKGROUND: The mouse inbred line C57BL/6J is widely used in mouse genetics and its genome has been incorporated into many genetic reference populations. More recently large initiatives such as the International Knockout Mouse Consortium (IKMC) are using the C57BL/6N mouse strain to generate null alleles for all mouse genes. Hence both strains are now widely used in mouse genetics studies. Here we perform a comprehensive genomic and phenotypic analysis of the two strains to identify differences that may influence their underlying genetic mechanisms. RESULTS: We undertake genome sequence comparisons of C57BL/6J and C57BL/6N to identify SNPs, indels and structural variants, with a focus on identifying all coding variants. We annotate 34 SNPs and 2 indels that distinguish C57BL/6J and C57BL/6N coding sequences, as well as 15 structural variants that overlap a gene. In parallel we assess the comparative phenotypes of the two inbred lines utilizing the EMPReSSslim phenotyping pipeline, a broad based assessment encompassing diverse biological systems. We perform additional secondary phenotyping assessments to explore other phenotype domains and to elaborate phenotype differences identified in the primary assessment. We uncover significant phenotypic differences between the two lines, replicated across multiple centers, in a number of physiological, biochemical and behavioral systems. CONCLUSIONS: Comparison of C57BL/6J and C57BL/6N demonstrates a range of phenotypic differences that have the potential to impact upon penetrance and expressivity of mutational effects in these strains. Moreover, the sequence variants we identify provide a set of candidate genes for the phenotypic differences observed between the two strains.

Gene loss and evolutionary rates following whole-genome duplication in teleost fishes.

Teleost fishes provide the first unambiguous support for ancient whole-genome duplication in an animal lineage. Studies in yeast or plants have shown that the effects of such duplications can be mediated by a complex pattern of gene retention and changes in evolutionary pressure. To explore such patterns in fishes, we have determined by phylogenetic analysis the evolutionary origin of 675 Tetraodon duplicated genes assigned to chromosomes, using additional data from other species of actinopterygian fishes. The subset of genes, which was retained in double after the genome duplication, is enriched in development, signaling, behavior, and regulation functional categories. The evolutionary rate of duplicate fish genes appears to be determined by 3 forces: 1) fish proteins evolve faster than mammalian orthologs; 2) the genes kept in double after genome duplication represent the subset under strongest purifying selection; and 3) following duplication, there is an asymmetric acceleration of evolutionary rate in one of the paralogs. These results show that similar mechanisms are at work in fishes as in yeast or plants and provide a framework for future investigation of the consequences of duplication in fishes and other animals.

Pharmacokinetic and pharmacodynamic analysis of efavirenz dose reduction using an in vitro-in vivo extrapolation model.

The pharmacokinetics (PK) of efavirenz (EFV) is characterized by marked interpatient variability that correlates with its pharmacodynamics (PD). In vitro-in vivo extrapolation (IVIVE) is a "bottom-up" approach that combines drug data with system information to predict PK and PD. The aim of this study was to simulate EFV PK and PD after dose reductions. At the standard dose, the simulated probability was 80% for viral suppression and 28% for central nervous system (CNS) toxicity. After a dose reduction to 400 mg, the probabilities of viral suppression were reduced to 69, 75, and 82%, and those of CNS toxicity were 21, 24, and 29% for the 516 GG, 516 GT, and 516 TT genotypes, respectively. With reduction of the dose to 200 mg, the probabilities of viral suppression decreased to 54, 62, and 72% and those of CNS toxicity decreased to 13, 18, and 20% for the 516 GG, 516 GT, and 516 TT genotypes, respectively. These findings indicate how dose reductions might be applied in patients with favorable genetic characteristics.

Functional interaction between the estrogen receptor and CTF1: analysis of the vitellogenin gene B1 promoter in yeast.

Eukaryotic gene expression depends on a complex interplay between the transcriptional apparatus and chromatin structure. We report here a yeast model system for investigating the functional interaction between the human estrogen receptor (hER) and CTF1, a member of the CTF/NFI transcription factor family. We show that a CTF1-fusion protein and the hER transactivate a synthetic promoter in yeast in a synergistic manner. This interaction requires the proline-rich transactivation domain of CTF1. When the natural estrogen-dependent vitellogenin B1 promoter is tested in yeast, CTF1 and CTF1-fusion proteins are unable to activate transcription, and no synergy is observed between hER, which activates the B1 promoter, and these factors. Chromatin structure analysis on this promoter reveals positioned nucleosomes at -430 to -270 (+/-20 bp) and at -270 to - 100 (+/-20 bp) relative to the start site of transcription. The positions of the nucleosomes remain unchanged upon hormone-dependent transcriptional activation of the promoter, and the more proximal nucleosome appears to mask the CTF/NFI site located at - 101 to -114. We conclude that a functional interaction of hER with the estrogen response element located upstream of a basal promoter occurs in yeast despite the nucleosomal organization of this promoter, whereas the interaction of CTF1 with its target site is apparently precluded by a nucleosome.

Genetic and karyotypic structure in the shrews of the Sorex araneus group: are they independent?

The species of the common shrew (Sorex araneus) group are morphologically very similar but exhibit high levels of karyotypic variation. Here we used genetic variation at 10 microsatellite markers in a data set of 212 individuals mostly sampled in the western Alps and composed of five karyotypic taxa (Sorex coronatus, Sorex antinorii and the S. araneus chromosome races Cordon, Bretolet and Vaud) to investigate the concordance between genetic and karyotypic structure. Bayesian analysis confirmed the taxonomic status of the three sampled species since individuals consistently grouped according to their taxonomical status. However, introgression can still be detected between S. antinorii and the race Cordon of S. araneus. This observation is consistent with the expected low karyotypic complexity of hybrids between these two taxa. Geographically based cryptic substructure was discovered within S. antinorii, a pattern consistent with the different postglaciation recolonization routes of this species. Additionally, we detected two genetic groups within S. araneus notwithstanding the presence of three chromosome races. This pattern can be explained by the probable hybrid status of the Bretolet race but also suggests a relatively low impact of chromosomal differences on genetic structure compared to historical factors. Finally, we propose that the current data set (available at http://www.unil.ch/dee/page7010_en.html#1) could be used as a reference by those wanting to identify Sorex individuals sampled in the western Alps.

In vitro susceptibility of Actinobaculum schaalii to 12 antimicrobial agents and molecular analysis of fluoroquinolone resistance.

OBJECTIVES: To assess the in vitro susceptibility of Actinobaculum schaalii to 12 antimicrobial agents as well as to dissect the genetic basis of fluoroquinolone resistance. METHODS: Forty-eight human clinical isolates of A. schaalii collected in Switzerland and France were studied. Each isolate was identified by 16S rRNA sequencing. MICs of amoxicillin, ceftriaxone, gentamicin, vancomycin, clindamycin, linezolid, ciprofloxacin, levofloxacin, moxifloxacin, co-trimoxazole, nitrofurantoin and metronidazole were determined using the Etest method. Interpretation of results was made according to EUCAST clinical breakpoints. The quinolone-resistance-determining regions (QRDRs) of gyrA and parC genes were also identified and sequence analysis was performed for all 48 strains. RESULTS: All isolates were susceptible to amoxicillin, ceftriaxone, gentamicin, clindamycin (except three), vancomycin, linezolid and nitrofurantoin, whereas 100% and 85% were resistant to ciprofloxacin/metronidazole and co-trimoxazole, respectively. Greater than or equal to 90% of isolates were susceptible to the other tested fluoroquinolones, and only one strain was highly resistant to levofloxacin (MIC ?32 mg/L) and moxifloxacin (MIC 8 mg/L). All isolates that were susceptible or low-level resistant to levofloxacin/moxifloxacin (n?=?47) showed identical GyrA and ParC amino acid QRDR sequences. In contrast, the isolate exhibiting high-level resistance to levofloxacin and moxifloxacin possessed a unique mutation in GyrA, Ala83Val (Escherichia coli numbering), whereas no mutation was present in ParC. CONCLUSIONS: When an infection caused by A. schaalii is suspected, there is a risk of clinical failure by treating with ciprofloxacin or co-trimoxazole, and ?-lactams should be preferred. In addition, acquired resistance to fluoroquinolones more active against Gram-positive bacteria is possible.

Clinal variation in a brown lemur (Eulemur spp.) hybrid zone: Combining morphological, genetic and climatic data to examine stability.

Studies of hybrid zones can inform our understanding of reproductive isolation and speciation. Two species of brown lemur (Eulemur rufifrons and E. cinereiceps) form an apparently stable hybrid zone in the Andringitra region of south-eastern Madagascar. The aim of this study was to identify factors that contribute to this stability. We sampled animals at 11 sites along a 90-km transect through the hybrid zone and examined variation in 26 microsatellites, the D-loop region of mitochondrial DNA, six pelage and nine morphological traits; we also included samples collected in more distant allopatric sites. Clines in these traits were noncoincident, and there was no increase in either inbreeding coefficients or linkage disequilibrium at the centre of the zone. These results could suggest that the hybrid zone is maintained by weak selection against hybrids, conforming to either the tension zone or geographical selection-gradient model. However, a closer examination of clines in pelage and microsatellites indicates that these clines are not sigmoid or stepped in shape but instead plateau at their centre. Sites within the hybrid zone also occur in a distinct habitat, characterized by greater seasonality in precipitation and lower seasonality in temperature. Together, these findings suggest that the hybrid zone may follow the bounded superiority model, with exogenous selection favouring hybrids within the transitional zone. These findings are noteworthy, as examples supporting the bounded superiority model are rare and may indicate a process of ecologically driven speciation without geographical isolation.