Population‐genomic approach reveals adaptive floral divergence in discrete populations of a hawk moth‐pollinated violet

Local adaptation to contrasting biotic or abiotic environments is an important evolutionary step that presumably precedes floral diversification at the species level, yet few studies have demonstrated the adaptive nature of intraspecific floral divergence in wild plant populations. We combine a population‐genomic approach with phenotypic information on floral traits to examine whether the differentiation in metric floral traits exhibited by 14 populations of the southern Spanish hawk moth‐pollinated violet Viola cazorlensis reflects adaptive divergence. Screening of many amplified fragment length polymorphism (AFLP) loci using a multiple‐marker‐based neutrality test identified nine outlier loci (2.6% of the total) that departed from neutral expectations and were potentially under selection. Generalized analysis of molecular variance revealed significant relationships between genetic distance and population divergence in three floral traits when genetic distance was based on outlier loci, but not when it was based on neutral ones. Population means of floral traits were closely correlated with population scores on the first principal coordinate axis of the genetic distance matrix using outlier loci, and with the allelic frequencies of four of the outlier loci. Results strongly support the adaptive nature of intraspecific floral divergence exhibited by V. cazorlensis and illustrate the potential of genome scans to identify instances of adaptive divergence when used in combination with phenotypic information.     

Reduced introgression of the Y chromosome between subspecies of the European rabbit (Oryctolagus cuniculus) in the Iberian Peninsula

The role of the Y chromosome in speciation is unclear. Hybrid zones provide natural arenas for studying speciation, as differential introgression of markers may reveal selection acting against incompatibilities. Two subspecies of the European rabbit (Oryctolagus cuniculus) form a hybrid zone in the Iberian Peninsula. Previous work on mitochondrial DNA (mtDNA), Y- and X-linked loci revealed the existence of two divergent lineages in the rabbit genome and that these lineages are largely subspecies-specific for mtDNA and two X-linked loci. Here we investigated the geographic distribution of the two Y chromosome lineages by genotyping two diagnostic single nucleotide polymorphisms in a sample of 353 male rabbits representing both subspecies, and found that Y chromosome lineages are also largely subspecies-specific. We then sequenced three autosomal loci and discovered considerable variation in levels of differentiation at these loci. Finally, we compared estimates of population differentiation between rabbit subspecies at 26 markers and found a surprising bimodal distribution of F(ST)values. The vast majority of loci showed little or no differentiation between rabbit subspecies while a few loci, including the SRY gene, showed little or no introgression across the hybrid zone. Estimates of population differentiation for the Y chromosome were surprisingly high given that there is male-biased dispersal in rabbits. Taken together, these data indicate that there is a clear dichotomy in the rabbit genome and that some loci remain highly differentiated despite extensive gene flow following secondary contact.     

Comparative estimation of effective population sizes and temporal gene flow in two contrasting population systems

Estimation of effective population sizes (N(e)) and temporal gene flow (N(e)m, m) has many implications for understanding population structure in evolutionary and conservation biology. However, comparative studies that gauge the relative performance of N(e), N(e)m or m methods are few. Using temporal genetic data from two salmonid fish population systems with disparate population structure, we (i) evaluated the congruence in estimates and precision of long- and short-term N(e), N(e)m and m from six methods; (ii) explored the effects of metapopulation structure on N(e) estimation in one system with spatiotemporally linked subpopulations, using three approaches; and (iii) determined to what degree interpopulation gene flow was asymmetric over time. We found that long-term N(e) estimates exceeded short-term N(e) within populations by 2-10 times; the two were correlated in the system with temporally stable structure (Atlantic salmon, Salmo salar) but not in the highly dynamic system (brown trout, Salmo trutta). Four temporal methods yielded short-term N(e) estimates within populations that were strongly correlated, and these were higher but more variable within salmon populations than within trout populations. In trout populations, however, these short-term N(e) estimates were always lower when assuming gene flow than when assuming no gene flow. Linkage disequilibrium data generally yielded short-term N(e) estimates of the same magnitude as temporal methods in both systems, but the two were uncorrelated. Correlations between long- and short-term geneflow estimates were inconsistent between methods, and their relative size varied up to eightfold within systems. While asymmetries in gene flow were common in both systems (58-63% of population-pair comparisons), they were only temporally stable in direction within certain salmon population pairs, suggesting that gene flow between particular populations is often intermittent and/or variable. Exploratory metapopulation N(e) analyses in trout demonstrated both the importance of spatial scale in estimating N(e) and the role of gene flow in maintaining genetic variability within subpopulations. Collectively, our results illustrate the utility of comparatively applying N(e), N(e)m and m to (i) tease apart processes implicated in population structure, (ii) assess the degree of continuity in patterns of connectivity between population pairs and (iii) gauge the relative performance of different approaches, such as the influence of population subdivision and gene flow on N(e) estimation. They further reiterate the importance of temporal sampling replication in population genetics, the value of interpreting N(e)or m in light of species biology, and the need to address long-standing assumptions of current N(e), N(e)m or m models more explicitly in future research.     

A 12 000-rad whole-genome radiation hybrid panel in sheep: application to the study of the ovine chromosome 18 region containing a QTL for scrapie susceptibility

Whole-genome radiation hybrid (RH) panels have been constructed for several species, including cattle. RH panels have proven to be an extremely powerful tool to construct high-density maps, which is an essential step in the identification of genes controlling important traits, and they can be used to establish high-resolution comparative maps. Although bovine RH panels can be used with ovine markers to construct sheep RH maps based on bovine genome organization, only some (c. 50%) of the markers available in sheep can be successfully mapped in the bovine genome. So, with the development of genomics and genome sequencing projects, there is a need for a high-resolution RH panel in sheep to map ovine markers. Consequently, we have constructed a 12 000-rad ovine whole-genome RH panel. Two hundred and eight hybrid clones were produced, of which 90 were selected based on their retention frequency. The final panel had an average marker retention frequency of 31.8%. The resolution of this 12 000-rad panel (SheepRH) was estimated by constructing an RH framework map for a 23-Mb region of sheep chromosome 18 (OAR18) that contains a QTL for scrapie susceptibility.     

Reproductive isolation between chromosomal races of the house mouse Mus musculus domesticus in a parapatric contact area revealed by an analysis of multiple unlinked loci

The house mouse, Mus musculus domesticus, exhibits a high level of chromosomal polymorphism because of the occurrence and fast fixation of Robertsonian fusions between telocentric chromosomes. For this reason, it has been considered a classical speciation model to analyse the role of the chromosomal changes in reproductive isolation. In this study, we analysed a parapatric contact area between two metacentric races in central Italy, the Cittaducale race (CD: 2n = 22) and the Ancarano race (ACR: 2n = 24), to estimate gene flow at the boundary. Hybrids between these two races show high levels of structural heterozygosity and are expected to be highly infertile. A sample of 88 mice from 14 sites was used. The mice were genotyped by means of eight microsatellite loci mapped in four different autosomal arms. The results show clear genetic differentiation between the CD and ACR races, as revealed by differences in allele frequencies, factorial correspondence analysis and indexes of genetic population (e.g. F(ST) and R(ST)) along the contact zone. The genetic differentiation between the races was further highlighted by assignation and clustering analyses, in which all the individuals were correctly assigned by their genotypes to the source chromosomal race. This result is particularly interesting in view of the absence of any geographical or ecological barrier in the parapatric contact zone, which occurs within a village. In these conditions, the observed genetic separation suggests an absence of gene flow between the races. The CD-ACR contact area is a rare example of a final stage of speciation between chromosomal races of rodents because of their chromosomal incompatibility.