The analysis of introgression of genomic regions between divergent populations provides an excellent opportunity to determine the genetic basis of reproductive isolation during the early stages of speciation. However, hybridization and subsequent gene flow must be relatively common in order to localize individual loci that resist introgression. In this study, we used next-generation sequencing to study genome-wide patterns of genetic differentiation between two hybridizing subspecies of rabbits (Oryctolagus cuniculus algirus and O. c. cuniculus) that are known to undergo high rates of gene exchange. Our primary objective was to identify specific genes or genomic regions that have resisted introgression and are likely to confer reproductive barriers in natural conditions. On the basis of 326,000 polymorphisms, we found low to moderate overall levels of differentiation between subspecies, and fewer than 200 genomic regions dispersed throughout the genome showing high differentiation consistent with a signature of reduced gene flow. Most differentiated regions were smaller than 200 Kb and contained very few genes. Remarkably, 30 regions were each found to contain a single gene, facilitating the identification of candidate genes underlying reproductive isolation. This gene-level resolution yielded several insights into the genetic basis and architecture of reproductive isolation in rabbits. Regions of high differentiation were enriched on the X-chromosome and near centromeres. Genes lying within differentiated regions were often associated with transcription and epigenetic activities, including chromatin organization, regulation of transcription, and DNA binding. Overall, our results from a naturally hybridizing system share important commonalities with hybrid incompatibility genes identified using laboratory crosses in mice and flies, highlighting general mechanisms underlying the maintenance of reproductive barriers.
Most research into the genetic underpinnings of speciation is based on laboratory studies of specific reproductively isolating phenotypes [1], [2]. A small number of these studies have identified individual genes that confer lower fitness when in a foreign or admixed genetic background [3]–[6]. However, it is not clear whether these genes and corresponding phenotypes - often identified in crosses between highly divergent species pairs - were ever relevant to restricting gene flow in nature, or whether they simply represent an inevitable byproduct of functional divergence between species [1], [7]. Moreover, it is not clear that hybrid incompatibilities accumulated long after the cessation of gene flow are strongly representative of the early stages of speciation. The direct analysis of naturally hybridizing populations provides an alternative and arguably more direct way to study the genetic basis of reproductive isolation [8]. This approach assures relevance to the early stages of speciation and allows the fitness of hybrid genotypes to be evaluated under natural conditions [9], [10]. Natural hybridization has been extensively studied in the last few years revealing the mosaic nature of gene flow across species boundaries [11]. It is now well established that genomic regions harboring hybrid incompatibilities are expected to display low levels of introgression, whereas regions not harboring incompatibilities should introgress more freely [12]–[14]. In spite of these general expectations, the identification of individual genes causing reproductive isolation in nature has remained elusive for at least three reasons. First, many hybrid zone studies have lacked mapping resolution due to limited sampling of the genome. Second, numerous historical and biological factors can inherently limit the fine-scale resolution of individual incompatibilities within many hybrid zones. For example, time since secondary contact, rates of dispersal and hybridization, and genomic architecture will all influence the effective rate of recombination and thus mapping resolution within a hybrid zone. Third, in many species, widespread differentiation across much of the genome limits the ability to detect loci contributing to isolation. In principle, some of these limitations can be overcome by applying high-throughput sequencing to the study of introgression between partially-isolated populations that are characterized by a long history of extensive gene flow. Genome-wide data from such populations promises to provide insight into the genetic basis of reproductive isolation in nature [15]–[19] and thereby provide a powerful counterpart to laboratory studies of speciation. Two naturally hybridizing subspecies of the European rabbit (Oryctolagus cuniculus algirus and O. c. cuniculus) are a good model for the early stages of reproductive isolation. These two subspecies currently occur in parapatry, hybridize over a large area in the central part of the Iberian Peninsula (Figure S1), and have diverged approximately 1. 8 million years ago [20], [21]. The current secondary contact is thought to result from a range expansion associated with the amelioration of climatic conditions after the last glacial maximum (∼18000 ya) [20], although previous glacial/interglacial cycles during the Pleistocene may have provided multiple opportunities for vicariant and past hybridization events. Both O. c. algirus and O. c. cuniculus span the three major climatic zones in the Iberian Peninsula (Mediterranean, Oceanic, and semi-arid) with little evidence for morphological differentiation along these ecological transitions [22], . The two subspecies can hybridize in the lab and produce viable hybrids, but hybrid males suffer from reduced fertility (Ferrand, unpublished results). Genetic studies using samples along the distribution range of both subspecies and based on a few dozen nuclear markers have revealed loci of relatively high divergence (0. 3–1. 2%) between subspecies embedded in a genome otherwise characterized by low levels of differentiation and high levels of bidirectional gene flow [21], [24]–[26], likely facilitated by high effective population sizes, high dispersal, and a relatively short generation time [21], [24]–[26]. Differentiation appears most pronounced on the sex chromosomes and in centromeric regions, albeit based on limited and biased genomic sampling. A more recent study using individuals sampled across the hybrid zone and 22 markers showed that some loci are characterized by concordant and stepped clines that do not co-localize with any obvious physical or ecological barriers [27], suggesting an overall barrier to gene flow likely to be primarily maintained by selection against intrinsically unfit hybrids. This same study revealed large inter-locus variation in cline width and cline center among some loci, consistent with substantial variation in patterns and levels of introgression as indicated by previous work [21], [24]–[26]. Long tails of introgression were often detected such that populations from the extreme ends of the Iberian Peninsula, which are well removed from the primary hybrid zone, were not genetically pure [27] (Figure S2). Thus, introgression seems not only to affect a large portion of the genome but frequently occurs through the entire range of parental populations. These previous studies provided important information regarding the dynamics of contact and hybridization. However, they were not designed to delineate regions of the genome resisting introgression nor investigate the frequency and types of genes involved in genetic incompatibilities, because they were based on a relatively small number of loci (less than 50) and the sampling of the genome was non-random (e. g. overrepresentation of the X-chromosome or regions near centromeres). The combination of extensive introgression between rabbit subspecies together with the availability of a rabbit reference genome sequence provides ideal conditions to resolve the genetic basis and architecture of persistent reproductive barriers between hybridizing taxa in the early stages of divergence. In this study, we used Illumina sequencing to provide a high-resolution map of genetic differentiation between rabbit subspecies. Our main goal was to identify specific genes or genomic regions that have resisted introgression and are likely to confer reproductive barriers in natural conditions. We further investigated the proportion of the genome that is resisting introgression and the number and size of such regions. We also tested whether previously hypothesized positional effects on genetic differentiation between rabbit subspecies (i. e. centromeres and X-chromosome) are robust to a more dense and unbiased sampling of the genome. We used custom microarrays to capture and Illumina sequence 6 Mb of targeted intronic DNA, and combined this dataset with published protein coding sequences derived from transcriptome sequencing [28]–[29]. The combination of genome-wide data and population genetic inference allowed the fine-scale identification of genes likely to underlie reproductive barriers in nature.
We analyzed a large dataset consisting of ∼326,000 SNPs (∼177,000 from the targeted DNA capture of intronic sequence and ∼149,000 from the transcriptome resequencing) distributed on all rabbit chromosomes, providing an unbiased picture of genetic differentiation between rabbit subspecies at the genome-wide level. This dataset was obtained in rabbits sampled at some distance from the primary contact zone (Figure S1), but in geographic regions that nonetheless show high levels of introgression [21], [24]–[27] (see Figure S2). We found that overall levels of genetic differentiation between the rabbit subspecies were low to moderate (Figure 1). Fixed differences in our dataset accounted for 0. 3% of the total number of mutations, and shared polymorphisms were ∼100 fold more numerous (31. 2%). The percentage of fixed differences in nature may be even lower given the relatively small number of individuals used in our study. We estimated an average fixation index (FST) value of 0. 084 across all polymorphisms and of 0. 169 across polymorphisms with a minor allele frequency (MAF) threshold of 10% (132,199 SNPs). Thus, the overall baseline of differentiation is moderate even after excluding low frequency variants, which are more likely than high-frequency variants to represent sequencing and genotyping errors [30]. The mean levels of differentiation reported here (FST autosomes = 0. 082; FST X-chromosome = 0. 185) were lower than previously reported values using Sanger sequencing of small intronic fragments (FST autosomes = 0. 146; FST X-chromosome = 0. 447) [26] and remained lower for the X-chromosome after the removal of low frequency variants (FST autosomes MAF>10% = 0. 164; FST X-chromosome MAF>10% = 0. 368). This difference likely reflects the biased sampling of the genome in previous studies towards regions of higher differentiation (e. g. centromeric regions). For comparison, our genome-wide estimates in rabbits are similar to that between human population groups, namely Africans, Caucasians, and Asians (FST≈0. 120) [31]. There are very few examples in the literature describing such low levels of genetic differentiation in the context of hybridizing taxa showing partial reproductive barriers. Anopheles mosquitoes provide perhaps one of the most extreme examples. Recent results describe numerous regions of differentiation between forms but large portions of the genome remain lowly differentiated [16], [32], [33]. However, they are thought to have diverged very recently - within the past 10,000 years - associated with the development of agriculture. In contrast, genetic divergence at multiple loci suggests that rabbits started diverging approximately 1. 8 million years ago [20], [21]. A more comparable example is found in sunflowers. A recent study found that genetic differentiation was much lower throughout a large fraction of the genome between the now sympatric species pair Helianthus annuus and H. petiolaris, which have diverged for ∼1. 0 million years, than between more closely related species with non-overlapping geographical ranges [34], [35]. The shared variation between rabbit subspecies is likely due to both gene flow and retention of ancestral variation. Evidence for high levels of gene flow comes from several sources. First, previous estimates of gene flow between rabbit subspecies using an Isolation-with-Migration model [36], [37] and a similar sampling scheme to the one used here suggest moderate to high levels of gene exchange (2Nm≈1. 2 for the autosomes averaged in both directions) [26]. These estimates are at the higher end of the range of values reported in a recent review of similar studies of hybridizing taxa [11]. Second, the spatial distribution of allelic frequencies across the hybrid zone revealed a subset of coincident, narrow and stepped clines consistent with reproductive barriers acting to prevent introgression in a fraction of the genome, but also that cline width varied by a factor of 50 with introgressed alleles often reaching the distribution ends of both subspecies [27] (Figure S2). Evidence for unsorted ancestral variation comes from a consideration of the divergence time and the effective population size (Ne). Effective population sizes (Ne) in rabbits are on the order of 106 [21]. Since the average coalescence time for alleles within a population is 4Ne generations for the autosomes and 3Ne for the X-chromosome, ancestral variation is still expected to be segregating between subspecies, even in the absence of gene flow (Ne≈106, divergence time≈1. 8 Mya, and a generation time of one year) [21], [38]. It is exceptionally difficult to resolve the relative contribution of these two forces in structuring genetic variation between the subspecies. Nonetheless, given high levels of gene flow (2Nm>1) and persistence of ancestral variation, fixed differences between the subspecies are expected to be very rare across the genome. To quantify this expectation given the number of individuals sampled, we conducted simulations to estimate the expected number of fixed differences between the subspecies assuming previous estimates of the degree of gene flow [21], . We found that both the autosomes and X-chromosome were significantly enriched for fixed differences relative to the simulated neutral expectation (P<0. 0001; Figure S3). Thus, genome-wide patterns of variation are characterized by extremely low levels of differentiation punctuated by an unusually high number of fixed differences. Given this pattern, genomic regions that are enriched for fixed differences are likely to be maintained by selection against introgression and thus contain genes involved in hybrid incompatibilities (or loci experiencing selection at linked sites within subspecies; see below). To investigate patterns of differentiation across the genome we first used a 100 kb sliding window approach (Figure 2; values are detailed in Dataset S1). Levels of differentiation varied greatly across the genome. Similar to our results for individual SNPs, the great majority of windows showed low to moderate differentiation (median FST MAF>10% = 0. 139), and extended portions of the genome were completely devoid of fixed differences. Several windows, however, were strongly differentiated relative to the genome-wide average. For the following analyses, we used SNPs with MAF>10%, although similar results are obtained when all SNPs are used (see Materials and Methods). For each window, we calculated the Z-score or the number of standard deviations departing from the median FST over all windows, the proportion of fixed differences versus shared polymorphisms, and also the absolute number of fixed differences. 3. 2% of windows showed a Z-score greater than three, which corresponded to an FST value of 0. 590 or higher. Intervals defined by a high Z-score also tended to show an enrichment of the absolute number of fixed differences and the proportion of fixed differences relative to shared polymorphisms (Figure 2). We found a 100% overlap between windows falling in the top 1% of the distribution of the absolute number of fixed differences and windows showing a significantly elevated proportion of fixed differences to shared polymorphisms when compared to the genome-wide average. We then used a randomized permutation test across SNPs within the genome (1000 replicates) to test if the clustering in levels of differentiation was unusual, given the genome-wide distribution of FST MAF>10% values for individual SNPs and fixed differences. None of the random permutations exceeded the observed differentiation (P<0. 001), indicating that these regions of differentiation are unlikely to reflect simple stochastic variation in differentiation across the genome due to the reduced number of individuals sampled. The spatial distribution of differentiation across the genome is potentially informative about the genetic architecture of reproductive isolation. To further demarcate contiguous genomic regions of high differentiation (i. e. , genomic islands of differentiation), we identified all single or consecutive windows with significantly elevated ratios of fixed differences to shared polymorphisms. Using these criteria, we identified 140 independent regions of differentiation that together encompass ∼46. 8 Mb of genomic sequence (∼1. 8% of the genome; Table S1). Alternatively, we identified 173 regions (60. 6 Mb and ∼2. 3% of the genome; Table S2) when we used a Z-score of three or higher to define highly differentiated regions. Importantly, we found a high degree of overlap between these two metrics of differentiation (102 overlapping regions) (Table S3) and most differentiated regions were well demarcated (i. e. appearing as single peaks surrounded by background levels of much lower differentiation; Figure 2). Our approach to identifying islands of divergence relies, in part, in levels of variation within populations and therefore is potentially sensitive to any process that reduces diversity within either subspecies, and thus high divergence may not be directly associated with signatures of reduced gene flow due to incompatibilities. For example, regions of high differentiation may be identified in our study due to selection at linked sites acting independently in each subspecies reducing within population variation relative to divergence between populations [39], [40]. We used two separate approaches to address this issue. First, we used a measure of differentiation, relative node depth (RND) [41], which does not depend on a within-population component of variation and therefore should be useful for distinguishing between both scenarios [42], [43]. RND takes into account possible mutation rate differences among genomic regions by incorporating divergence to a third more distantly related species, and is expected to be inversely proportional to the amount of gene flow. To generate appropriate outgroup data, we used a custom microarray to capture and sequence the same 5,000 introns used for rabbits in a single Hare (Lepus timidus) and then calculated RND for each region. We found much higher variance in RND values for intronic fragments consisting of a reduced number of sites, likely due to stochastic variation in the mutation process. To reduce sampling noise, we restricted our analyses to the 1,665 fragments with more than 1. 0 kb of aligned Lepus data. Consistent with variable levels of introgression playing a substantial role in determining patterns of genomic differentiation, we found a significant increase in mean RND values for fragments residing within regions of differentiation defined both using the Z-score or the proportion of fixed differences versus shared polymorphisms (P<0. 05 for all comparisons; Figure 3). Second, we tested for evidence of more frequent selection within islands of differentiation. If regions of high differentiation reflect the independent action of positive selection within either subspecies (i. e. , subspecific selective sweeps), then we would expect an excess of high frequency derived alleles within regions of differentiation [44]. We did not find significant differences or deviations in a consistent direction in the frequency of derived alleles in islands of divergence within O. c. cuniculus or in either subspecies when defining differentiation based on the proportion of fixed differences versus shared polymorphisms (Figure 3); however, intervals of high FST in O. c. algirus did show a significant skew towards high frequency derived alleles. Our results do not rule out that a simple hitchhiking model or a combination of both factors (i. e. hitchhiking plus reduced introgression) may explain some of the regions of differentiation, particular within O. c. algirus. Indeed, several admixture scenarios [44], [45] and selection at linked sites are expected to generate an excess of high frequency derived alleles, but the observation that most islands of differentiation are not enriched for high frequency derived alleles suggests that selection at linked sites is likely to have a relatively small contribution overall to the occurrence of areas of high differentiation. Moreover, the fact that the inferred skew towards high frequency derived alleles was only detected in one subspecies is mostly consistent with asymmetric migration. Taken together, our results indicate that both FST (i. e. Z-score) and the ratio of shared polymorphisms to fixed differences are capturing signatures of reduced gene flow and that regions of high differentiation indentified in our study are likely to be enriched for regions involved in reproductive barriers. Given that combining different summary statistics is likely to reduce the number of false positives, for the remainder of the analysis we defined highly differentiated regions as the 102 genomic intervals that show significant differentiation for both metrics. Our results further show that signatures of reduced introgression are spread across the genome, suggesting that reproductive barriers in rabbits appear to have a complex genetic basis controlled by many unlinked genes. Genome-wide patterns of isolation do not provide direct information regarding the mechanisms or phenotypes that contribute to reproductive barriers, underscoring both the power of this population genetic approach because it is not constrained by choice of a specific phenotype, but also its limitations. Information regarding hybrid phenotypes in rabbits is scarce because wild animals are exceptionally difficult to maintain in captivity. Nonetheless preliminary data from crosses between subspecies indicate that hybrid males are less fertile than hybrid females or parental individuals. Hybrid male sterility often evolves very quickly and may show a complex genetic basis even at the earliest stages of divergence [e. g. 46], which could account for reduced introgression at multiple regions of the genome as inferred here. Yet, we cannot discard that other hybrid phenotypes, both resulting from intrinsic or extrinsic causes, also act as reproductive barriers between rabbit subspecies in nature. The physical size of differentiated regions was highly heterogeneous (Figure 4). The largest region of differentiation spanned more than ∼1. 3 Mb but the majority of differentiation occurred on a physical scale smaller than 200 kb. The large number of small and independent regions of elevated differentiation in a background of low differentiation indicates that the time since secondary contact and subsequent hybridization has been sufficient for recombination to decouple most of the genome from adjacent isolation factors. Since the average gene density in a mammalian genome is ∼1 gene per 100 Kb, this scale of differentiation approaches gene-level resolution (see below). Differentiated regions and their size were not randomly distributed across the genome. First, the X-chromosome exhibited particularly strong genetic differentiation (Figure 2). Mean FST values were significantly higher on the X-chromosome (FST MAF>10% = 0. 368) when compared to the autosomes (FST MAF>10% = 0. 164; Mann-Whitney U test, P<0. 001), as was the ratio of fixed to shared polymorphisms (0. 283 vs 0. 007; Fisher' s Exact Test, P<0. 001). This elevated differentiation was reflected in a significantly larger size of differentiated regions on the X than that on the autosomes (MeanAutosomes = 202 Kb; MeanX-chromosome = 379 Kb; Mann-Whitney U test P<0. 001; Figure 4). Moreover, while the X-chromosome represents ∼5% of the rabbit genome, it contained a significant enrichment of differentiated regions (27. 5%, P<0. 001). In contrast, the mean RND value was significantly lower for the X-chromosome (RND = 0. 176) when compared to the autosomes (RND = 0. 209; P<0. 001). A similar pattern between chimpanzees and humans has been interpreted as evidence for complex speciation scenarios involving differential gene flow [47], but several authors have pointed out that several other neutral and simpler explanations were not considered [e. g. 48]. RND values are corrected for mutation rate differences among loci using divergence to an outgroup; however, the X-chromosomes spends only one third of its evolutionary history in males. Thus, variation in the strength of male-biased mutation across time due to changes in reproductive system or generation time in the rabbits versus Lepus comparison could contribute in part to the autosome/X-chromosome difference. Additionally, the ancestral Ne is expected to be smaller for the X chromosome than for the autosomes and this is expected to lead to lower RND for the former than for the latter. Second, we observed strong differentiation close to several centromeric regions, which are known to have reduced rates of recombination in several species [e. g. 49]. Mean FST and RND were significantly higher within 5 Mb of centromeres (FST MAF>10% = 0. 246; RND = 0. 228) when compared to the remainder of the genome (FST MAF>10% = 0. 164; RND = 0. 213; Mann-Whitney U test, P<0. 05 for both tests). Two of the most notable such regions are found on chromosomes 4 and X (Figure 2). Regions of differentiation found within 5 Mb of centromeres were larger when compared to the rest of the genome (Meancentromeres = 339 kb; Meangenome = 225 kb; Mann-Whitney U test P<0. 07; Figure 5), and significantly more numerous since they represent ∼9% of the genome and contain 22. 5% of the total number of differentiated regions (P<0. 01). We note that differentiation was not universally elevated near centromeres and many differentiated regions were found along chromosomal arms (Figure 2), including some larger than 500 Kb (Figure 5). These genome-wide data confirm and extend previously described positional effects on genetic differentiation in rabbits [21], [24]–[26]. These findings underscore the central role that the X chromosome plays in the evolution of reproductive isolation [2], [50], and support the theoretical prediction that low recombination might facilitate species divergence in the face of gene flow [reviewed in 51]. According to these hypotheses, these regions should be enriched for genes affecting reproductive isolation. However, we cannot exclude the possibility that some of the regions of higher differentiation on the X or in centromeric regions are driven by selection at linked sites or reduced Ne [39], [40]. The availability of an annotated rabbit genome provides a good opportunity to access candidate genes within differentiated regions. Overall, these regions contained 410 known or predicted genes (Table S3), and of these genes, 337 were annotated as protein-coding genes, 2 as retrotransposed genes, 57 as noncoding RNAs, and 14 as pseudogenes. Consistent with the variation in size we also found a heterogeneous number of genes contained within each region (ranging from 1–31). Nine regions (8. 8%) contained more than 10 genes each. In such regions it will be difficult to identify candidate variants underlying reproductive isolation. However, most differentiated regions contained very few annotated genes: 59 segments (57. 8%) contained three or fewer genes and 30 (29. 4%) localized to a single gene. Thus, while many regions throughout the genome appear to contribute to reproductive isolation, the substrate of selection within most regions is likely to be one or a few genes. These findings further underscore the utility of the rabbit hybrid zone for fine-scale mapping of candidate genes likely to be involved in the early stages of reproductive isolation. Due to the reduced number of genes within most islands of isolation, our results allow us to explore the functional underpinnings of putative reproductive barriers at a higher resolution than most comparable genomic studies. We found that genes within islands of differentiation were enriched for several functional classes (Table S4), and most terms were related to transcription and epigenetic pathways (e. g. chromatin organization and modification, regulation of transcription and translation, DNA binding). These results are consistent with the emerging paradigm that epigenetic and transcription regulation frequently underlies the evolution of postzygotic hybrid dysfunctions [4], [5], [52], [53]. Contrary to several case studies on hybrid incompatibility genes [54]–[57], we found no statistical evidence for more rapid protein evolution within islands of differentiation (Table 1). Single-gene islands provide the strongest candidates for hybrid incompatibilities. Given this, we performed additional examination of their molecular functions (Table 2). We found a wide variety of molecular functions among these candidates, and interactions between distinct elements (e. g. DNA binding, RNA binding, and protein binding), and several transcription activities (e. g. transcription corepressor activity, transcription factor binding, and histone binding), were particularly common. Moreover, several candidate genes were associated with male sterility, which is the only hybrid phenotype detected so far in crosses between rabbit subspecies. For example, EIF4G3 encodes a eukaryotic translation initiation factor, and mutations in this gene have been shown to cause male limited infertility in mice by inducing meiosis and spermatogenesis arrest in meiotic prophase [58]. CDKN2C is a member of the INK4 family of cyclin-dependent kinase inhibitors and its expression is largely confined to spermatocytes undergoing meiosis in seminiferous tubules. Male mice bearing a null CDKN2C mutant suffer from increased germ cell apoptosis, resulting in abnormal spermatogenesis and oligozoospermia [59]. Finally, when we consider our findings in the context of hybrid male sterility genes previously uncovered in laboratory crosses, some intriguing commonalities are apparent. Prdm9 was the first hybrid male sterility gene discovered in vertebrates and was identified from crosses between mouse subspecies [5]. It contains a KRAB motif, a histone H3 Lysine-4-methyltransferase domain, and several zinc fingers that likely mediate sequence-specific binding to DNA. OdsH, a hybrid male sterility gene identified from crosses between Drosophila mauritiana and D. simulans, is thought to encode a heterochromatin-binding protein, and differential DNA binding properties between alleles result in altered heterochromatic localization and chromosome decondensation [4]. Ovd, another hybrid male sterility gene in Drosophila, is also thought to function as a DNA-binding protein [6]. It is noteworthy that several genes in our candidate gene list have nucleic-acid binding properties, and at least five genes contain zinc-fingers and two interact with histones similarly to Prdm9 (Table 2). Crossing experiments and follow-up functional studies will ultimately be necessary to confirm the involvement of these genes in reproductive isolation and to determine the specific phenotypes which they are controlling. The findings presented here provide a manageable list of candidate genes with explicit phenotypic predictions.
We have used genome-wide data and a population genetic approach to dissect the genetic basis of reproductive isolation in the European rabbit by taking advantage of natural hybridization between subspecies in the Iberian Peninsula. While most of the genome is lowly differentiated between subspecies, we identified numerous regions of strong differentiation, suggesting that the genetic basis of reproductive isolation may be highly polygenic. In addition, the architecture of differentiation is such that these regions were of small size and contained very few genes. This observation and the availability of a reference rabbit genome allowed us to identify several candidate genes for reproductive isolation in this system that suggest an important role for epigenetic and transcription regulation in the maintenance of reproductive barriers. Most interestingly, the molecular functions of some genes identified here using patterns of genetic differentiation in nature strongly parallel those of the hybrid male sterility genes identified using laboratory crosses. Genome-wide analyses of differentiation, such as has been pursued here, are now feasible in a large number of organisms, and are poised to provide new genetic insights into the early stages of speciation.
A number of laboratory studies of speciation have uncovered individual genes that confer lower fitness when in a divergent genetic background. It is, however, unclear whether these genes were ever relevant to restricting gene flow in nature, or whether they are indirect consequences of functional divergence between species that in most cases can no longer hybridize in natural conditions. Analysis of introgression across the genome of divergent populations provides an alternative approach to determine the genetic basis of reproductive isolation. Using two subspecies of rabbits as a model for the early stages of speciation, we provide a genome-wide map of genetic differentiation. Our study revealed important aspects of the genetic architecture of differentiation in the early stages of divergence and allowed the identification of genomic regions that resist introgression and are likely to confer reproductive barriers in natural conditions. The gene content of these regions, which in several cases reached gene-level resolution, suggests an important role for epigenetic and transcription regulation in the maintenance of reproductive barriers. Some of the genes identified here in a natural system are similar in function to the hybrid male sterility genes identified in laboratory studies of speciation.
genome sequencing
genome analysis tools
gene ontologies
genomics
genome scans
speciation
introgression
population genetics
hybridization
biology
evolutionary biology
gene flow
evolutionary processes
transcriptomes
The Genomic Architecture of Population Divergence between Subspecies of the European Rabbit