An unmatched case controlled study of clinicopathologic abnormalities in dogs with Bartonella infection

We compared clinicopathologic findings in dogs with Bartonella infection to Bartonella spp. negative dogs suspected of a vector-borne disease. Cases (n=47) and controls (n=93) were selected on the basis of positive or negative enrichment culture PCR results, respectively. Signalment, clinicopathologic findings and treatments were extracted from medical records. DNA sequencing identified Bartonella henselae (n=28, 59.6%), Bartonella vinsonii subsp. berkhoffii (n=20, 42.6%), Bartonella koehlerae (n=3, 6.4%), Bartonella volans-like (n=3, 6.4%) and Bartonella bovis (n=1, 2.1%). There were no significant differences in age, breed, size, sex or neuter status between cases and controls. Dogs infected with Bartonella sp. often had a history of weight loss [OR=2.82; 95% CI: 1.08-7.56] and were hypoglobulinemic [OR=4.26; 95% CI: 1.31-14.41]. With the exception of weight loss and hypoglobulinemia, clinicopathologic abnormalities in Bartonella-infected dogs in this study were similar to dogs suspected of other vector-borne infections.

Spectrum and prevalence of mutations involving BrS1- through BrS12-susceptibility genes in a cohort of unrelated patients referred for Brugada syndrome genetic testing: implications for genetic testing

OBJECTIVES The aim of this study was to provide the spectrum and prevalence of mutations in the 12 Brugada syndrome (BrS)-susceptibility genes discovered to date in a single large cohort of unrelated BrS patients. BACKGROUND BrS is a potentially lethal heritable arrhythmia syndrome diagnosed electrocardiographically by coved-type ST-segment elevation in the right precordial leads (V1 to V3; type 1 Brugada electrocardiographic [ECG] pattern) and the presence of a personal/family history of cardiac events. METHODS Using polymerase chain reaction, denaturing high-performance liquid chromatography, and DNA sequencing, comprehensive mutational analysis of BrS1- through BrS12-susceptibility genes was performed in 129 unrelated patients with possible/probable BrS (46 with clinically diagnosed BrS [ECG pattern plus personal/family history of a cardiac event] and 83 with a type 1 BrS ECG pattern only). RESULTS Overall, 27 patients (21%) had a putative pathogenic mutation, absent in 1,400 Caucasian reference alleles, including 21 patients with an SCN5A mutation, 2 with a CACNB2B mutation, and 1 each with a KCNJ8 mutation, a KCND3 mutation, an SCN1Bb mutation, and an HCN4 mutation. The overall mutation yield was 23% in the type 1 BrS ECG pattern-only patients versus 17% in the clinically diagnosed BrS patients and was significantly greater among young men<20 years of age with clinically diagnosed BrS and among patients who had a prolonged PQ interval. CONCLUSIONS We identified putative pathogenic mutations in ∼20% of our BrS cohort, with BrS genes 2 through 12 accounting for <5%. Importantly, the yield was similar between patients with only a type 1 BrS ECG pattern and those with clinically established BrS. The yield approaches 40% for SCN5A-mediated BrS (BrS1) when the PQ interval exceeds 200 ms. Calcium channel-mediated BrS is extremely unlikely in the absence of a short QT interval.

Cardiac channel molecular autopsy: insights from 173 consecutive cases of autopsy-negative sudden unexplained death referred for postmortem genetic testing

OBJECTIVE To perform long QT syndrome and catecholaminergic polymorphic ventricular tachycardia cardiac channel postmortem genetic testing (molecular autopsy) for a large cohort of cases of autopsy-negative sudden unexplained death (SUD). METHODS From September 1, 1998, through October 31, 2010, 173 cases of SUD (106 males; mean ± SD age, 18.4 ± 12.9 years; age range, 1-69 years; 89% white) were referred by medical examiners or coroners for a cardiac channel molecular autopsy. Using polymerase chain reaction, denaturing high-performance liquid chromatography, and DNA sequencing, a comprehensive mutational analysis of the long QT syndrome susceptibility genes (KCNQ1, KCNH2, SCN5A, KCNE1, and KCNE2) and a targeted analysis of the catecholaminergic polymorphic ventricular tachycardia type 1-associated gene (RYR2) were conducted. RESULTS Overall, 45 putative pathogenic mutations absent in 400 to 700 controls were identified in 45 autopsy-negative SUD cases (26.0%). Females had a higher yield (26/67 [38.8%]) than males (19/106 [17.9%]; P<.005). Among SUD cases with exercise-induced death, the yield trended higher among the 1- to 10-year-olds (8/12 [66.7%]) compared with the 11- to 20-year-olds (4/27 [14.8%]; P=.002). In contrast, for those who died during a period of sleep, the 11- to 20-year-olds had a higher yield (9/25 [36.0%]) than the 1- to 10-year-olds (1/24 [4.2%]; P=.01). CONCLUSION Cardiac channel molecular autopsy should be considered in the evaluation of autopsy-negative SUD. Several interesting genotype-phenotype observations may provide insight into the expected yields of postmortem genetic testing for SUD and assist in selecting cases with the greatest potential for mutation discovery and directing genetic testing efforts.

Unexplained drownings and the cardiac channelopathies: a molecular autopsy series.

OBJECTIVE To determine the prevalence and spectrum of mutations associated with long QT syndrome (LQTS) and catecholaminergic polymorphic ventricular tachycardia (CPVT) in a seemingly unexplained drowning cohort. PATIENTS AND METHODS From September 1, 1998, through October 31, 2010, 35 unexplained drowning victims (23 male and 12 female; mean ± SD age, 17±12 years [range, 4-69 years]) were referred for a cardiac channel molecular autopsy. Of these, 28 (20 male and 8 female) drowned while swimming, and 7 (3 male and 4 female) were bathtub submersions. Polymerase chain reaction, denaturing high-performance liquid chromatography, and DNA sequencing were used for a comprehensive mutational analysis of the 3 major LQTS-susceptibility genes (KCNQ1, KCNH2, and SCN5A), and a targeted analysis of the CPVT1-associated, RYR2-encoded cardiac ryanodine receptor was conducted. RESULTS Of the 28 victims of swimming-related drowning, 8 (28.6%) were mutation positive, including 2 with KCNQ1 mutations (L273F, AAPdel71-73 plus V524G) and 6 with RYR2 mutations (R414C, I419F, R1013Q, V2321A, R2401H, and V2475F). None of the bathtub victims were mutation positive. Of the 28 victims who drowned while swimming, women were more likely to be mutation positive than men (5/8 [62.5%] vs 3/20 [15%]; P=.02). Although none of the mutation-positive, swimming-related drowning victims had a premortem diagnosis of LQTS or CPVT, a family history of cardiac arrest, family history of prior drowning, or QT prolongation was present in 50%. CONCLUSION Nearly 30% of the victims of swimming-related drowning hosted a cardiac channel mutation. Genetic testing should be considered in the postmortem evaluation of an unexplained drowning, especially if a positive personal or family history is elicited.

Loss-of-function mutations in the KCNJ8-encoded Kir6.1 K(ATP) channel and sudden infant death syndrome.

BACKGROUND Approximately 10% of sudden infant death syndrome (SIDS) may stem from cardiac channelopathies. The KCNJ8-encoded Kir6.1 (K(ATP)) channel critically regulates vascular tone and cardiac adaptive response to systemic metabolic stressors, including sepsis. KCNJ8-deficient mice are prone to premature sudden death, particularly with infection. We determined the spectrum, prevalence, and function of KCNJ8 mutations in a large SIDS cohort. METHODS AND RESULTS Using polymerase chain reaction, denaturing high-performance liquid chromatography, and DNA sequencing, comprehensive open reading frame/splice-site mutational analysis of KCNJ8 was performed on genomic DNA isolated from necropsy tissue on 292 unrelated SIDS cases (178 males, 204 white; age, 2.9±1.9 months). KCNJ8 mutations were coexpressed heterologously with SUR2A in COS-1 cells and characterized using whole-cell patch-clamp. Two novel KCNJ8 mutations were identified. A 5-month-old white male had an in-frame deletion (E332del) and a 2-month-old black female had a missense mutation (V346I). Both mutations localized to Kir6.1's C-terminus, involved conserved residues and were absent in 400 and 200 ethnic-matched reference alleles respectively. Both cases were negative for mutations in established channelopathic genes. Compared with WT, the pinacidil-activated K(ATP) current was decreased 45% to 68% for Kir6.1-E332del and 40% to 57% for V346I between -20 mV and 40 mV. CONCLUSIONS Molecular and functional evidence implicated loss-of-function KCNJ8 mutations as a novel pathogenic mechanism in SIDS, possibly by predisposition of a maladaptive cardiac response to systemic metabolic stressors akin to the mouse models of KCNJ8 deficiency.

Gain-of-function mutation S422L in the KCNJ8-encoded cardiac K(ATP) channel Kir6.1 as a pathogenic substrate for J-wave syndromes.

BACKGROUND J-wave syndromes have emerged conceptually to encompass the pleiotropic expression of J-point abnormalities including Brugada syndrome (BrS) and early repolarization syndrome (ERS). KCNJ8, which encodes the cardiac K(ATP) Kir6.1 channel, recently has been implicated in ERS following identification of the functionally uncharacterized missense mutation S422L. OBJECTIVE The purpose of this study was to further explore KCNJ8 as a novel susceptibility gene for J-wave syndromes. METHODS Using polymerase chain reaction, denaturing high-performance liquid chromatography, and direct DNA sequencing, comprehensive open reading frame/splice site mutational analysis of KCNJ8 was performed in 101 unrelated patients with J-wave syndromes, including 87 with BrS and 14 with ERS. Six hundred healthy individuals were examined to assess the allelic frequency for all variants detected. KCNJ8 mutation(s) was engineered by site-directed mutagenesis and coexpressed heterologously with SUR2A in COS-1 cells. Ion currents were recorded using whole-cell configuration of the patch-clamp technique. RESULTS One BrS case and one ERS case hosted the identical missense mutation S422L, which was reported previously. KCNJ8-S422L involves a highly conserved residue and was absent in 1,200 reference alleles. Both cases were negative for mutations in all known BrS and ERS susceptibility genes. K(ATP) current of the Kir6.1-S422L mutation was increased significantly over the voltage range from 0 to 40 mV compared to Kir6.1-WT channels (n = 16-21; P <.05). CONCLUSION These findings further implicate KCNJ8 as a novel J-wave syndrome susceptibility gene and a marked gain of function in the cardiac K(ATP) Kir6.1 channel secondary to KCNJ8-S422L as a novel pathogenic mechanism for the phenotypic expression of both BrS and ERS.

Sudden infant death syndrome-associated mutations in the sodium channel beta subunits.

BACKGROUND Approximately 10% of sudden infant death syndrome (SIDS) cases may stem from potentially lethal cardiac channelopathies, with approximately half of channelopathic SIDS involving the Na(V)1.5 cardiac sodium channel. Recently, Na(V) beta subunits have been implicated in various cardiac arrhythmias. Thus, the 4 genes encoding Na(V) beta subunits represent plausible candidate genes for SIDS. OBJECTIVE This study sought to determine the spectrum, prevalence, and functional consequences of sodium channel beta-subunit mutations in a SIDS cohort. METHODS In this institutional review board-approved study, mutational analysis of the 4 beta-subunit genes, SCN1B to 4B, was performed using polymerase chain reaction, denaturing high-performance liquid chromatography, and direct DNA sequencing of DNA derived from 292 SIDS cases. Engineered mutations were coexpressed with SCN5A in HEK 293 cells and were whole-cell patch clamped. One of the putative SIDS-associated mutations was similarly studied in adenovirally transduced adult rat ventricular myocytes. RESULTS Three rare (absent in 200 to 800 reference alleles) missense mutations (beta3-V36M, beta3-V54G, and beta4-S206L) were identified in 3 of 292 SIDS cases. Compared with SCN5A+beta3-WT, beta3-V36M significantly decreased peak I(Na) and increased late I(Na), whereas beta3-V54G resulted in a marked loss of function. beta4-S206L accentuated late I(Na) and positively shifted the midpoint of inactivation compared with SCN5A+beta4-WT. In native cardiomyocytes, beta4-S206L accentuated late I(Na) and increased the ventricular action potential duration compared with beta4-WT. CONCLUSION This study provides the first molecular and functional evidence to implicate the Na(V) beta subunits in SIDS pathogenesis. Altered Na(V)1.5 sodium channel function due to beta-subunit mutations may account for the molecular pathogenic mechanism underlying approximately 1% of SIDS cases.

Alpha1-syntrophin mutations identified in sudden infant death syndrome cause an increase in late cardiac sodium current.

BACKGROUND Sudden infant death syndrome (SIDS) is a leading cause of death during the first 6 months after birth. About 5% to 10% of SIDS may stem from cardiac channelopathies such as long-QT syndrome. We recently implicated mutations in alpha1-syntrophin (SNTA1) as a novel cause of long-QT syndrome, whereby mutant SNTA1 released inhibition of associated neuronal nitric oxide synthase by the plasma membrane Ca-ATPase PMCA4b, causing increased peak and late sodium current (I(Na)) via S-nitrosylation of the cardiac sodium channel. This study determined the prevalence and functional properties of SIDS-associated SNTA1 mutations. METHODS AND RESULTS Using polymerase chain reaction, denaturing high-performance liquid chromatography, and DNA sequencing of SNTA1's open reading frame, 6 rare (absent in 800 reference alleles) missense mutations (G54R, P56S, T262P, S287R, T372M, and G460S) were identified in 8 (approximately 3%) of 292 SIDS cases. These mutations were engineered using polymerase chain reaction-based overlap extension and were coexpressed heterologously with SCN5A, neuronal nitric oxide synthase, and PMCA4b in HEK293 cells. I(Na) was recorded using the whole-cell method. A significant 1.4- to 1.5-fold increase in peak I(Na) and 2.3- to 2.7-fold increase in late I(Na) compared with controls was evident for S287R-, T372M-, and G460S-SNTA1 and was reversed by a neuronal nitric oxide synthase inhibitor. These 3 mutations also caused a significant depolarizing shift in channel inactivation, thereby increasing the overlap of the activation and inactivation curves to increase window current. CONCLUSIONS Abnormal biophysical phenotypes implicate mutations in SNTA1 as a novel pathogenic mechanism for the subset of channelopathic SIDS. Functional studies are essential to distinguish pathogenic perturbations in channel interacting proteins such as alpha1-syntrophin from similarly rare but innocuous ones.

Unique mixed phenotype and unexpected functional effect revealed by novel compound heterozygosity mutations involving SCN5A.

BACKGROUND Functional characterization of mutations involving the SCN5A-encoded cardiac sodium channel has established the pathogenic mechanisms for type 3 long QT syndrome and type 1 Brugada syndrome and has provided key insights into the physiological importance of essential structure-function domains. OBJECTIVE This study sought to present the clinical and biophysical phenotypes discerned from compound heterozygosity mutations in SCN5A on different alleles in a toddler diagnosed with QT prolongation and fever-induced ventricular arrhythmias. METHODS A 22-month-old boy presented emergently with fever and refractory ventricular tachycardia. Despite restoration of sinus rhythm, the infant sustained profound neurological injury and died. Using polymerase chain reaction, denaturing high-performance liquid chromatography, and direct DNA sequencing, comprehensive open-reading frame/splice mutational analysis of the 12 known long QT syndrome susceptibility genes was performed. RESULTS The infant had 2 SCN5A mutations: a maternally inherited N-terminal frame shift/deletion (R34fs/60) and a paternally inherited missense mutation, R1195H. The mutations were engineered by site-directed mutagenesis and heterologously expressed transiently in HEK293 cells. As expected, the frame-shifted and prematurely truncated peptide, SCN5A-R34fs/60, showed no current. SCN5A-R1195H had normal peak and late current but abnormal voltage-dependent gating parameters. Surprisingly, co-expression of SCN5A-R34fs/60 with SCN5A-R1195H elicited a significant increase in late sodium current, whereas co-expression of SCN5A-WT with SCN5A-R34fs/60 did not. CONCLUSIONS A severe clinical phenotype characterized by fever-induced monomorphic ventricular tachycardia and QT interval prolongation emerged in a toddler with compound heterozygosity involving SCN5A: R34fs/60, and R1195H. Unexpectedly, the 94-amino-acid fusion peptide derived from the R34fs/60 mutation accentuated the late sodium current of R1195H-containing Na(V)1.5 channels in vitro.

Extended-spectrum cephalosporin-resistant Escherichia coli in community, specialized outpatient clinic and hospital settings in Switzerland

OBJECTIVES Resistance to extended-spectrum cephalosporins (ESCs) in Escherichia coli can be due to the production of ESBLs, plasmid-mediated AmpCs (pAmpCs) or chromosomal AmpCs (cAmpCs). Information regarding type and prevalence of β-lactamases, clonal relations and plasmids associated with the bla genes for ESC-R E. coli (ESC-R-Ec) detected in Switzerland is lacking. Moreover, data focusing on patients referred to the specialized outpatient clinics (SOCs) are needed. METHODS We analysed 611 unique E. coli isolated during September-December 2011. ESC-R-Ec were studied with microarrays, PCR/DNA sequencing for blaESBLs, blapAmpCs, promoter region of blacAmpC, IS elements, plasmid incompatibility group, and also implementing transformation, aIEF, rep-PCR and MLST. RESULTS The highest resistance rates were observed in the SOCs, whereas those in the hospital and community were lower (e.g. quinolone resistance of 22.6%, 17.2% and 9.0%, respectively; P = 0.003 for SOCs versus community). The prevalence of ESC-R-Ec in the three settings was 5.3% (n = 11), 7.8% (n = 22) and 5.7% (n = 7), respectively. Thirty isolates produced CTX-M ESBLs (14 were CTX-M-15), 5 produced CMY-2 pAmpC and 5 hyper-expressed cAmpCs due to promoter mutations. Fourteen isolates were of sequence type 131 (ST131; 10 with CTX-M-15). blaCTX-M and blaCMY-2 were associated with an intact or truncated ISEcp1 and were mainly carried by IncF, IncFII and IncI1plasmids. CONCLUSIONS ST131 producing CTX-M-15 is the predominant clone. The prevalence of ESC-R-Ec (overall 6.5%) is low, but an unusual relatively high frequency of AmpC producers (25%) was noted. The presence of ESC-R-Ec in the SOCs and their potential ability to be exchanged between hospital and community should be taken into serious consideration.

Non-phenotypic tests to detect and characterize antibiotic resistance mechanisms in Enterobacteriaceae

In the past 2 decades, we have observed a rapid increase of infections due to multidrug-resistant Enterobacteriaceae. Regrettably, these isolates possess genes encoding for extended-spectrum β-lactamases (e.g., blaCTX-M, blaTEM, blaSHV) or plasmid-mediated AmpCs (e.g., blaCMY) that confer resistance to last-generation cephalosporins. Furthermore, other resistance traits against quinolones (e.g., mutations in gyrA and parC, qnr elements) and aminoglycosides (e.g., aminoglycosides modifying enzymes and 16S rRNA methylases) are also frequently co-associated. Even more concerning is the rapid increase of Enterobacteriaceae carrying genes conferring resistance to carbapenems (e.g., blaKPC, blaNDM). Therefore, the spread of these pathogens puts in peril our antibiotic options. Unfortunately, standard microbiological procedures require several days to isolate the responsible pathogen and to provide correct antimicrobial susceptibility test results. This delay impacts the rapid implementation of adequate antimicrobial treatment and infection control countermeasures. Thus, there is emerging interest in the early and more sensitive detection of resistance mechanisms. Modern non-phenotypic tests are promising in this respect, and hence, can influence both clinical outcome and healthcare costs. In this review, we present a summary of the most advanced methods (e.g., next-generation DNA sequencing, multiplex PCRs, real-time PCRs, microarrays, MALDI-TOF MS, and PCR/ESI MS) presently available for the rapid detection of antibiotic resistance genes in Enterobacteriaceae. Taking into account speed, manageability, accuracy, versatility, and costs, the possible settings of application (research, clinic, and epidemiology) of these methods and their superiority against standard phenotypic methods are discussed.

Occurrence and Genetic Characteristics of Third-Generation Cephalosporin-Resistant Escherichia coli in Swiss Retail Meat

Prevalence and genetic relatedness were determined for third-generation cephalosporin-resistant Escherichia coli (3GC-R-Ec) detected in Swiss beef, veal, pork, and poultry retail meat. Samples from meat-packing plants (MPPs) processing 70% of the slaughtered animals in Switzerland were purchased at different intervals between April and June 2013 and analyzed. Sixty-nine 3GC-R-Ec isolates were obtained and characterized by microarray, PCR/DNA sequencing, Multi Locus Sequence Typing (MLST), and plasmid replicon typing. Plasmids of selected strains were transformed by electroporation into E. coli TOP10 cells and analyzed by plasmid MLST. The prevalence of 3GC-R-Ec was 73.3% in chicken and 2% in beef meat. No 3GC-R-Ec were found in pork and veal. Overall, the blaCTX-M-1 (79.4%), blaCMY-2 (17.6%), blaCMY-4 (1.5%), and blaSHV-12 (1.5%) β-lactamase genes were detected, as well as other genes conferring resistance to chloramphenicol (cmlA1-like), sulfonamides (sul), tetracycline (tet), and trimethoprim (dfrA). The 3GC-R-Ec from chicken meat often harbored virulence genes associated with avian pathogens. Plasmid incompatibility (Inc) groups IncI1, IncFIB, IncFII, and IncB/O were the most frequent. A high rate of clonality (e.g., ST1304, ST38, and ST93) among isolates from the same MPPs suggests that strains persist at the plant and spread to meat at the carcass-processing stage. Additionally, the presence of the blaCTX-M-1 gene on an IncI1 plasmid sequence type 3 (IncI1/pST3) in genetically diverse strains indicates interstrain spread of an epidemic plasmid. The blaCMY-2 and blaCMY-4 genes were located on IncB/O plasmids. This study represents the first comprehensive assessment of 3GC-R-Ec in meat in Switzerland. It demonstrates the need for monitoring contaminants and for the adaptation of the Hazard Analysis and Critical Control Point concept to avoid the spread of multidrug-resistant bacteria through the food chain.

The brown hare (Lepus europaeus) as a novel intermediate host for Echinococcus multilocularis in Europe

A typical multivesiculated metacestode tissue has been found in the liver of a European brown hare (Lepus europaeus) originating from a northern area of Switzerland. In this study, the causative species was identified as Echinococcus multilocularis by appropriate histological and molecular analyses and corresponding DNA sequencing. This is the first confirmation of larval E. multilocularis from hares in central Europe. The metacestode tissue contained protoscolices, suggesting that the hare may contribute to the transmission of E. multilocularis in Switzerland.

Chromosomally and extrachromosomally mediated high-level gentamicin resistance in Streptococcus agalactiae

Streptococcus agalactiae (group B Streptococcus, GBS) is a leading cause of sepsis in neonates. The rate of invasive GBS disease in non-pregnant adults also continues to climb. Aminoglycosides alone have little or no effect on GBS, but synergistic killing with penicillin has been shown in vitro. High-level gentamicin resistance (HLGR) in GBS isolates, however, leads to loss of a synergistic effect. We therefore performed a multicentre study to determine the frequency of HLGR GBS isolates and to elucidate the molecular mechanisms leading to gentamicin resistance. From eight centres in four countries, 1128 invasive and colonizing GBS isolates were pooled and investigated for the presence of HLGR. We identified two strains that displayed HLGR (BSU1203 and BSU452), both of which carried the aacA-aphD gene, typically conferring HLGR. Though, only one strain (BSU1203) also carried the previously described chromosomal gentamicin resistance transposon, designated Tn3706. In the other strain (BSU452), plasmid purification and subsequent DNA sequencing resulted in the detection of plasmid pIP501 carrying a remnant of a Tn3 family transposon. Its ability to confer HLGR was proven by transfer into an Enterococcus faecalis isolate. Conversely, loss of HLGR was documented after curing both GBS BSU452 and the transformed E. faecalis strain from the plasmid. This is the first report showing a plasmid mediated HLGR in GBS. Thus, in our clinical GBS isolates HLGR is mediated both chromosomally and extrachromosomally.

Genomics of Rapid Incipient Speciation in Sympatric Threespine Stickleback

Ecological speciation is the process by which reproductively isolated populations emerge as a consequence of divergent natural or ecologically-mediated sexual selection. Most genomic studies of ecological speciation have investigated allopatric populations, making it difficult to infer reproductive isolation. The few studies on sympatric ecotypes have focused on advanced stages of the speciation process after thousands of generations of divergence. As a consequence, we still do not know what genomic signatures of the early onset of ecological speciation look like. Here, we examined genomic differentiation among migratory lake and resident stream ecotypes of threespine stickleback reproducing in sympatry in one stream, and in parapatry in another stream. Importantly, these ecotypes started diverging less than 150 years ago. We obtained 34,756 SNPs with restriction-site associated DNA sequencing and identified genomic islands of differentiation using a Hidden Markov Model approach. Consistent with incipient ecological speciation, we found significant genomic differentiation between ecotypes both in sympatry and parapatry. Of 19 islands of differentiation resisting gene flow in sympatry, all were also differentiated in parapatry and were thus likely driven by divergent selection among habitats. These islands clustered in quantitative trait loci controlling divergent traits among the ecotypes, many of them concentrated in one region with low to intermediate recombination. Our findings suggest that adaptive genomic differentiation at many genetic loci can arise and persist in sympatry at the very early stage of ecotype divergence, and that the genomic architecture of adaptation may facilitate this.