Search

Accurate knowledge of geographic ranges and genetic relationships among populations is important when managing a species or population of conservation concern. Along the western coast of Canada, a subspecies of the northern goshawk (Accipiter gentilis laingi) is legally designated as Threatened. The range and distinctness of this form, in comparison with the broadly distributed North American subspecies (Accipiter gentilis atricapillus), is unclear. Given this morphological uncertainty, we analyzed genomic relationships in thousands of single nucleotide polymorphisms identified using genotyping‐by‐sequencing of high‐quality genetic samples. Results revealed a genetically distinct population of northern goshawks on the archipelago of Haida Gwaii and subtle structuring among other North American sampling regions. We then developed genotyping assays for ten loci that are highly differentiated between the two main genetic clusters, allowing inclusion of hundreds of low‐quality samples and confirming that the distinct genetic cluster is restricted to Haida Gwaii. As the laingi form was originally described as being based on Haida Gwaii (where the type specimen is from), further morphological analysis may result in this name being restricted to the Haida Gwaii genetic cluster. Regardless of taxonomic treatment, the distinct Haida Gwaii genetic cluster along with the small and declining population size of the Haida Gwaii population suggests a high risk of extinction of an ecologically and genetically distinct form of northern goshawk. Outside of Haida Gwaii, sampling regions along the coast of BC and southeast Alaska (often considered regions inhabited by laingi) show some subtle differentiation from other North American regions. These results will increase the effectiveness of conservation management of northern goshawks in northwestern North America. More broadly, other conservation‐related studies of genetic variation may benefit from the two‐step approach we employed that first surveys genomic variation using high‐quality samples and then genotypes low‐quality samples at particularly informative loci.

Biological invasions are accelerating, and invasive species can have large economic impacts as well as severe consequences for biodiversity. During invasions, species can interact, potentially resulting in hybridization. Here, we examined two Cakile species, C. edentula and C. maritima (Brassicaceae), that co-occur and may hybridize during range expansion in separate regions of the globe. Cakile edentula invaded each location first, while C. maritima established later, apparently replacing the former. We assessed the evidence for hybridization in western North America and Australia, where both species have been introduced, and identified source populations with 4561 SNPs using Genotype-by-Sequencing. Our results indicate that C. edentula in Australia originated from one region of eastern North America while in western North America it is probably from multiple sources. Cakile maritima in Australia is derived from at least two different parts of Europe while the introduction in western North America is from one. Although morphological evidence of hybridization is generally limited to mixed species populations in Australia and virtually absent elsewhere, our genetic analysis revealed relatively high levels of hybridization in Australia (58% hybrids using Admixture) and supported the presence of hybrids in western North America (16% hybrids using Admixture) and New Zealand. Hybrids might be commonly overlooked in invaders, as identification based solely on morphological traits may represent only the tip of the iceberg. Our study reveals a repeated pattern of invasion, hybridization and apparent replacement of one species by another, which offers an opportunity to investigate the role of hybridization and introgression during invasion.

Accurate knowledge of geographic ranges and genetic relationships among populations is important when managing a species or population of conservation concern. Along the western coast of Canada, a subspecies of the northern goshawk (Accipiter gentilis laingi) is legally designated as Threatened. The range and distinctness of this form, in comparison with the broadly distributed North American subspecies (Accipiter gentilis atricapillus), is unclear. Given this morphological uncertainty, we analyzed genomic relationships in thousands of single nucleotide polymorphisms identified using genotyping‐by‐sequencing of high‐quality genetic samples. Results revealed a genetically distinct population of northern goshawks on the archipelago of Haida Gwaii and subtle structuring among other North American sampling regions. We then developed genotyping assays for ten loci that are highly differentiated between the two main genetic clusters, allowing inclusion of hundreds of low‐quality samples and confirming that the distinct genetic cluster is restricted to Haida Gwaii. As the laingi form was originally described as being based on Haida Gwaii (where the type specimen is from), further morphological analysis may result in this name being restricted to the Haida Gwaii genetic cluster. Regardless of taxonomic treatment, the distinct Haida Gwaii genetic cluster along with the small and declining population size of the Haida Gwaii population suggests a high risk of extinction of an ecologically and genetically distinct form of northern goshawk. Outside of Haida Gwaii, sampling regions along the coast of BC and southeast Alaska (often considered regions inhabited by laingi) show some subtle differentiation from other North American regions. These results will increase the effectiveness of conservation management of northern goshawks in northwestern North America. More broadly, other conservation‐related studies of genetic variation may benefit from the two‐step approach we employed that first surveys genomic variation using high‐quality samples and then genotypes low‐quality samples at particularly informative loci.

Accurate knowledge of geographic ranges and genetic relationships among populations is important when managing a species or population of conservation concern. In the western Canadian province of British Columbia, a subspecies of the northern goshawk (Accipiter gentilis laingi) is designated as Threatened under the Canadian Species at Risk Act. Historically, the range of this bird of prey has been ambiguous and its genetic distinctness from the other North American subspecies (Accipiter gentilis atricapillus) has not been well established. Given the uncertainty in using morphological traits to assign individual goshawks to these two subspecies, we analyzed genomic relationships in tens of thousands of single nucleotide polymorphisms identified using genotyping-bysequencing of high-quality genetic samples. This genome-wide analysis revealed a genetically distinct population of northern goshawks on the archipelago of Haida Gwaii and subtle genetic structuring among the remainder of our sampling sites within North America. Following from this analysis, we developed targeted genotyping assays for ten loci that are highly differentiated between the two main genetic clusters, allowing the addition of hundreds of low-quality samples to our analysis. This additional information confirmed that the distinct genetic cluster on Haida Gwaii is restricted to that archipelago. As the laingi form was originally described as being based in Haida Gwaii, where the type specimen of that form is from, further study (especially of morphological traits) may indicate a need to restrict this name to the Haida Gwaii genetic cluster. Regardless of taxonomic treatment, our finding of a distinct Haida Gwaii genetic cluster along with the small and historically declining population size of the Haida Gwaii population suggests a high risk of extinction of an ecologically and genetically distinct form of northern goshawk. Outside of Haida Gwaii, sampling regions along the coast of BC and southeast Alaska (often considered regions inhabited by laingi) show some subtle differentiation from other North American regions. We anticipate that these results will increase the effectiveness of conservation management of northern goshawks in northwestern North America. More broadly, other conservation-related studies of genetic variation may benefit from the two-step approach we employed that first surveys genomic variation using high-quality samples and then genotypes low-quality samples at particularly informative loci.

In order to uncover the genetic basis of phenotypic trait variation, we used 448 unrelated wild accessions of black cottonwood (Populus trichocarpa) from much of its range in western North America. Extensive data from large-scale trait phenotyping (with spatial and temporal replications within a common garden) and genotyping (with a 34 K Populus single nucleotide polymorphism (SNP) array) of all accessions were used for gene discovery in a genome-wide association study (GWAS). We performed GWAS with 40 biomass, ecophysiology and phenology traits and 29,355 filtered SNPs representing 3518 genes. The association analyses were carried out using a Unified Mixed Model accounting for population structure effects among accessions. We uncovered 410 significant SNPs using a Bonferroni-corrected threshold (P

Summary The increasing ecological and economical importance of Populus species and hybrids has stimulated research into the investigation of the natural variation of the species and the estimation of the extent of genetic control over its wood quality traits for traditional forestry activities as well as the emerging bioenergy sector. A realized kinship matrix based on informative, high-density, biallelic single nucleotide polymorphism (SNP) genetic markers was constructed to estimate trait variance components, heritabilities, and genetic and phenotypic correlations. Seventeen traits related to wood chemistry and ultrastructure were examined in 334 9-yr-old Populus trichocarpa grown in a common-garden plot representing populations spanning the latitudinal range 44° to 58.6°. In these individuals, 9342 SNPs that conformed to Hardy–Weinberg expectations were employed to assess the genomic pair-wise kinship to estimate narrow-sense heritabilities and genetic correlations among traits. The range-wide phenotypic variation in all traits was substantial and several trait heritabilities were > 0.6. In total, 61 significant genetic and phenotypic correlations and a network of highly interrelated traits were identified. The high trait variation, the evidence for moderate to high heritabilities and the identification of advantageous trait combinations of industrially important characteristics should aid in providing the foundation for the enhancement of poplar tree breeding strategies for modern industrial use.

In the early stages of reproductive isolation, genomic regions of reduced recombination are expected to show greater levels of differentiation, either because gene flow between species is reduced in these regions or because the effects of selection at linked sites within species is enhanced in these regions. Here we study patterns of DNA sequence variation at 27 autosomal loci among populations of Mus musculus musculus, M. m. domesticus, and M. m. castaneus, three subspecies of house mice with co-linear genomes. We found that some loci exhibit considerable shared variation among subspecies, while others exhibit fixed differences. We used an isolation-with-gene-flow model to estimate divergence times, effective population sizes (Ne) and to disentangle ancestral variation from gene flow. Estimates of divergence time indicate that all three subspecies diverged from one another within a very short period of time roughly 350,000 years ago. Overall, Ne for each subspecies was associated with the degree of genetic differentiation: M. m. musculus had the smallest Ne and the greatest proportion of monophyletic gene genealogies, while M. m. castaneus had the largest Ne and the smallest proportion of monophyletic gene genealogies. M. m. domesticus and M. m. musculus were more differentiated from each other than either was from M. m. castaneus, consistent with greater reproductive isolation between M. m. domesticus and M. m. musculus. FST was significantly greater at loci experiencing low recombination rates compared to loci experiencing high recombination rates in comparisons between M. m. castaneus and M. m. musculus or M. m. domesticus. These results provide evidence that genomic regions with less recombination show greater differentiation, even in the absence of chromosomal rearrangements.

The western black cottonwood (Populus trichocarpa) was the first tree to have its genome fully sequenced and has emerged as the model species for the study of secondary growth and wood formation. It is also a good candidate species for the production of lignocellulosic biofuels. Here, we present and make available to the research community the results of the sequencing of the transcriptome of developing xylem in 20 accessions with high-throughput next generation sequencing technology. We found over 0.5 million putative single nucleotide polymorphisms (SNPs) in 26,595 genes that are expressed in developing secondary xylem. More than two-thirds of all SNPs were found in annotated exons, with 18% and 14% in regions of the genome annotated as introns and intergenic, respectively, where only 3% and 4% of sequence reads mapped. This suggests that the current annotation of the poplar genome is remarkably incomplete and that there are many transcripts and novel genes waiting to be annotated. We hope that this resource will stimulate further research in expression profiling, detection of alternative splicing and adaptive evolution in poplar.

International audience; Understanding the genetic structure of domestic species provides a window into the process of domestication and motivates the design of studies aimed at making links between genotype and phenotype. Rabbits exhibit exceptional phenotypic diversity, are of great commercial value, and serve as important animal models in biomedical research. Here, we provide the first comprehensive survey of nucleotide polymorphism and linkage disequilibrium (LD) within and among rabbit breeds. We resequenced 16 genomic regions in population samples of both wild and domestic rabbits and additional 35 fragments in 150 rabbits representing six commonly used breeds. Patterns of genetic variation suggest a single origin of domestication in wild populations from France, supporting historical records that place rabbit domestication in French monasteries. Levels of nucleotide diversity both within and among breeds were ;0.2%, but only 60% of the diversity present in wild populations from France was captured by domestic rabbits. Despite the recent origin of most breeds, levels of population differentiation were high (FST 5 17.9%), but the majority of polymorphisms were shared and thus transferable among breeds. Coalescent simulations suggest that domestication began with a small founding population of less than 1,200 individuals. Taking into account the complex demographic history of domestication with two successive bottlenecks, two loci showed deviations that were consistent with artificial selection, including GPC4, which is known to be associated with growth rates in humans. Levels of diversity were not significantly different between autosomal and X-linked loci, providing no evidence for differential contributions of males and females to the domesticated gene pool. The structure of LD differed substantially within and among breeds. Within breeds, LD extends over large genomic distances. Markers separated by 400 kb typically showed r2 higher than 0.2, and some LD extended up to 3,200 kb. Much less LD was found among breeds. This advantageous LD structure holds great promise for reducing the interval of association in future mapping studies.

The human Y chromosome consists of ampliconic genes, which are located in palindromes and undergo frequent gene conversion, and single-copy genes including the primary sex-determining locus, SRY. Here, we demonstrate that SRY is duplicated in a large palindrome in the European rabbit (Oryctolagus cuniculus). Furthermore, we show through comparative sequencing that orthologous palindrome arms have diverged 0.40% between rabbit subspecies over at least 2 My, but paralogous palindrome arms have remained nearly identical. This provides clear evidence of gene conversion on the rabbit Y chromosome. Together with previous observations in humans, these results suggest that gene conversion is a general feature of the evolution of the mammalian Y chromosome.

The house mouse is a well-established model organism, particularly for studying the genetics of complex traits. However, most studies of mice use classical inbred strains, whose genomes derive from multiple species. Relatively little is known about the distribution of genetic variation among these species or how variation among strains relates to variation in the wild. We sequenced intronic regions of five X-linked loci in large samples of wild Mus domesticus and M. musculus, and we found low levels of nucleotide diversity in both species. We compared these data to published data from short portions of six X-linked and 18 autosomal loci in wild mice. We estimate that M. domesticus and M. musculus diverged

Hybrid zones provide an excellent opportunity for studying the consequences of genetic changes between closely related taxa. Here we investigate patterns of genetic variability and gene flow at four X-linked loci within and between the two subspecies of European rabbit (Oryctolagus cuniculus cuniculus and O. c. algirus). Two of these genes are located near the centromere and two are located near the telomeres. We observed a deep split in the genealogy of each gene with the root located along the deepest branch in each case, consistent with the evolution of these subspecies in allopatry. The two centromeric loci showed low levels of variability, high levels of linkage disequilibrium, and little introgression between subspecies. In contrast, the two telomeric loci showed high levels of variability, low levels of linkage disequilibrium, and considerable introgression between subspecies. These data are consistent with suppression of recombination near the centromere of the rabbit X chromosome. These observations support a view of speciation where genomic incompatibilities at different loci in the genome create localized differences in levels of gene flow between nascent species.

All species of the genus Populus (poplar, aspen) are dioecious, suggesting an ancient origin of this trait. Despite some empirical counter examples, theory suggests that nonrecombining sex-linked regions should quickly spread, eventually becoming heteromorphic chromosomes. In contrast, we show using whole-genome scans that the sex-associated region in Populus trichocarpa is small and much younger than the age of the genus. This indicates that sex determination is highly labile in poplar, consistent with recent evidence of ‘turnover’ of sex-determination regions in animals. We performed whole-genome resequencing of 52 P. trichocarpa (black cottonwood) and 34 Populus balsamifera (balsam poplar) individuals of known sex. Genomewide association studies in these unstructured populations identified 650 SNPs significantly associated with sex. We estimate the size of the sex-linked region to be ~100 kbp. All SNPs significantly associated with sex were in strong linkage disequilibrium despite the fact that they were mapped to six different chromosomes (plus 3 unmapped scaffolds) in version 2.2 of the reference genome. We show that this is likely due to genome misassembly. The segregation pattern of sex-associated SNPs revealed this to be an XY sex-determining system. Estimated divergence times of X and Y haplotype sequences (6–7 Ma) are much more recent than the divergence of P. trichocarpa (poplar) and Populus tremuloides (aspen). Consistent with this, in P. tremuloides, we found no XY haplotype divergence within the P. trichocarpa sex-determining region. These two species therefore have a different genomic architecture of sex, suggestive of at least one turnover event in the recent past.

A ‘genes-to-ecosystems’ approach has been proposed as a novel avenue for integrating the consequences of intraspecific genetic variation with the underlying genetic architecture of a species to shed light on the relationships among hierarchies of ecological organization (genes → individuals → communities → ecosystems). However, attempts to identify genes with major effect on the structure of communities and/or ecosystem processes have been limited and a comprehensive test of this approach has yet to emerge. Here, we present an interdisciplinary field study that integrated a common garden containing different genotypes of a dominant, riparian tree, Populus trichocarpa, and aquatic mesocosms to determine how intraspecific variation in leaf litter alters both terrestrial and aquatic communities and ecosystem functioning. Moreover, we incorporate data from extensive trait screening and genome-wide association studies estimating the heritability and genes associated with litter characteristics. We found that tree genotypes varied considerably in the quality and production of leaf litter, which contributed to variation in phytoplankton abundances, as well as nutrient dynamics and light availability in aquatic mesocosms. These ‘after-life’ effects of litter from different genotypes were comparable to the responses of terrestrial communities associated with the living foliage. We found that multiple litter traits corresponding with aquatic community and ecosystem responses differed in their heritability. Moreover, the underlying genetic architecture of these traits was complex, and many genes contributed only a small proportion to phenotypic variation. Our results provide further evidence that genetic variation is a key component of aquatic–terrestrial linkages, but challenge the ability to predict community or ecosystem responses based on the actions of one or a few genes. Fil: Crutsinger, Gregory M.. University Of British Columbia; Canadá Fil: Rudman, Seth M.. University Of British Columbia; Canadá Fil: Rodriguez Cabal, Mariano Alberto. University Of British Columbia; Canadá. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Patagonia Norte. Instituto de Investigación en Biodiversidad y Medioambiente; Argentina Fil: Mckown, Athena D.. University Of British Columbia; Canadá Fil: Sato, Takuya. Kobe University; Fil: Macdonald, Andrew M.. University Of British Columbia; Canadá Fil: Heavyside, Julian. University Of British Columbia; Canadá Fil: Geraldes, Arnaldo. University Of British Columbia; Canadá Fil: Hart, Edmund M.. University Of British Columbia; Canadá Fil: Leroy, Carri J.. Evergreen State College; Canadá Fil: El-Sabaawi, Rana W.. University Of Victoria; Canadá

Populus trichocarpa is an ecologically important tree across western North America. We used a large population sample of 498 accessions over a wide geographical area genotyped with a 34K Populus SNP array to quantify geographical patterns of genetic variation in this species (landscape genomics). We present evidence that three processes contribute to the observed patterns: (1) introgression from the sister species P. balsamifera, (2) isolation by distance (IBD), and (3) natural selection. Introgression was detected only at the margins of the species' distribution. IBD was significant across the sampled area as a whole, but no evidence of restricted gene flow was detected in a core of drainages from southern British Columbia (BC). We identified a large number of FST outliers. Gene Ontology analyses revealed that FST outliers are overrepresented in genes involved in circadian rhythm and response to red/far-red light when the entire dataset is considered, whereas in southern BC heat response genes are overrepresented. We also identified strong correlations between geoclimate variables and allele frequencies at FST outlier loci that provide clues regarding the selective pressures acting at these loci.

• Populus trichocarpa is widespread across western North America spanning extensive variation in photoperiod, growing season and climate. We investigated trait variation in P. trichocarpa using over 2000 trees from a common garden at Vancouver, Canada, representing replicate plantings of 461 genotypes originating from 136 provenance localities. • We measured 40 traits encompassing phenological events, biomass accumulation, growth rates, and leaf, isotope and gas exchange-based ecophysiology traits. With replicated plantings and 29,354 single nucleotide polymorphisms (SNPs) from 3518 genes, we estimated both broad-sense trait heritability (H(2)) and overall population genetic structure from principal component analysis. • Populus trichocarpa had high phenotypic variation and moderate/high H(2) for many traits. H(2) ranged from 0.3 to 0.9 in phenology, 0.3 to 0.8 in biomass and 0.1 to 0.8 in ecophysiology traits. Most traits correlated strongly with latitude, maximum daylength and temperature of tree origin, but not necessarily with elevation, precipitation or heat : moisture indices. Trait H(2) values reflected trait correlation strength with geoclimate variables. The population genetic structure had one significant principal component (PC1) which correlated with daylength and showed enrichment for genes relating to circadian rhythm and photoperiod. • Robust relationships between traits, population structure and geoclimate in P. trichocarpa reflect patterns which suggest that range-wide geographical and environment gradients have shaped its genotypic and phenotypic variability.

Summary • Plant population genomics informs evolutionary biology, breeding, conservation and bioenergy feedstock development. For example, the detection of reliable phenotype–genotype associations and molecular signatures of selection requires a detailed knowledge about genome-wide patterns of allele frequency variation, linkage disequilibrium and recombination. • We resequenced 16 genomes of the model tree Populus trichocarpa and genotyped 120 trees from 10 subpopulations using 29 213 single-nucleotide polymorphisms. • Significant geographic differentiation was present at multiple spatial scales, and range-wide latitudinal allele frequency gradients were strikingly common across the genome. The decay of linkage disequilibrium with physical distance was slower than expected from previous studies in Populus, with r 2 dropping below 0.2 within 3–6 kb. Consistent with this, estimates of recent effective population size from linkage disequilibrium (Ne 4000–6000) were remarkably low relative to the large census sizes of P. trichocarpa stands. Fine-scale rates of recombination varied widely across the genome, but were largely predictable on the basis of DNA sequence and methylation features. • Our results suggest that genetic drift has played a significant role in the recent evolutionary history of P. trichocarpa. Most importantly, the extensive linkage disequilibrium detected suggests that genome-wide association studies and genomic selection in undomesticated populations may be more feasible in Populus than previously assumed.

Genetic mapping of quantitative traits requires genotypic data for large numbers of markers in many individuals. For such studies, the use of large single nucleotide polymorphism (SNP) genotyping arrays still offers the most cost-effective solution. Herein we report on the design and performance of a SNP genotyping array for Populus trichocarpa (black cottonwood). This genotyping array was designed with SNPs pre-ascertained in 34 wild accessions covering most of the species latitudinal range. We adopted a candidate gene approach to the array design that resulted in the selection of 34 131 SNPs, the majority of which are located in, or within 2 kb of, 3543 candidate genes. A subset of the SNPs on the array (539) was selected based on patterns of variation among the SNP discovery accessions. We show that more than 95% of the loci produce high quality genotypes and that the genotyping error rate for these is likely below 2%. We demonstrate that even among small numbers of samples (n = 10) from local populations over 84% of loci are polymorphic. We also tested the applicability of the array to other species in the genus and found that the number of polymorphic loci decreases rapidly with genetic distance, with the largest numbers detected in other species in section Tacamahaca. Finally, we provide evidence for the utility of the array to address evolutionary questions such as intraspecific studies of genetic differentiation, species assignment and the detection of natural hybrids.

Flax (Linum usitatissimum) is an ancient crop that is widely cultivated as a source of fiber, oil and medicinally relevant compounds. To accelerate crop improvement, we performed whole-genome shotgun sequencing of the nuclear genome of flax. Seven paired-end libraries ranging in size from 300 bp to 10 kb were sequenced using an Illumina genome analyzer. A de novo assembly, comprised exclusively of deep-coverage (approximately 94× raw, approximately 69× filtered) short-sequence reads (44-100 bp), produced a set of scaffolds with N50 = 694 kb, including contigs with N50 = 20.1 kb. The contig assembly contained 302 Mb of non-redundant sequence representing an estimated 81% genome coverage. Up to 96% of published flax ESTs aligned to the whole-genome shotgun scaffolds. However, comparisons with independently sequenced BACs and fosmids showed some mis-assembly of regions at the genome scale. A total of 43 384 protein-coding genes were predicted in the whole-genome shotgun assembly, and up to 93% of published flax ESTs, and 86% of A. thaliana genes aligned to these predicted genes, indicating excellent coverage and accuracy at the gene level. Analysis of the synonymous substitution rates (Ks) observed within duplicate gene pairs was consistent with a recent (5-9 MYA) whole-genome duplication in flax. Within the predicted proteome, we observed enrichment of many conserved domains (Pfam-A) that may contribute to the unique properties of this crop, including agglutinin proteins. Together these results show that de novo assembly, based solely on whole-genome shotgun short-sequence reads, is an efficient means of obtaining nearly complete genome sequence information for some plant species. © 2012 Blackwell Publishing Ltd.

The western black cottonwood (Populus trichocarpa) was the first tree to have its genome fully sequenced and has emerged as the model species for the study of secondary growth and wood formation. It is also a good candidate species for the production of lignocellulosic biofuels. Here, we present and make available to the research community the results of the sequencing of the transcriptome of developing xylem in 20 accessions with high-throughput next generation sequencing technology. We found over 0.5 million putative single nucleotide polymorphisms (SNPs) in 26,595 genes that are expressed in developing secondary xylem. More than two-thirds of all SNPs were found in annotated exons, with 18% and 14% in regions of the genome annotated as introns and intergenic, respectively, where only 3% and 4% of sequence reads mapped. This suggests that the current annotation of the poplar genome is remarkably incomplete and that there are many transcripts and novel genes waiting to be annotated. We hope that this resource will stimulate further research in expression profiling, detection of alternative splicing and adaptive evolution in poplar.

Reconstructing the history of populations is a longstanding goal of molecular ecologists. In addition to a better understanding of the past, it is hoped that this knowledge would also facilitate predictions regarding species' responses to future events such as climate change. The traditional way of doing this is through the fossil record, but these historical records are often incomplete. Inferring historical demography from patterns of nucleotide variability can help to fill these gaps. In this issue of Molecular Ecology, Holliday et al. (2010) glimpse into the demographic past of Sitka spruce, Picea sitchensis, an economically and ecologically important species native to northwestern United States and Canada, by examining the site frequency spectrum (SFS) of 153 loci in six populations covering the species entire range.

Patterns of genetic differentiation among taxa at early stages of divergence provide an opportunity to make inferences about the history of speciation. Here, we conduct a survey of DNA-sequence polymorphism and divergence at loci on the autosomes, X chromosome, Y chromosome and mitochondrial DNA in samples of Mus domesticus, M. musculus and M. castaneus. We analyzed our data under a divergence with gene flow model and estimate that the effective population size of M. castaneus is 200 000–400 000, of M. domesticus is 100 000–200 000 and of M. musculus is 60 000–120 000. These data also suggest that these species started to diverge approximately 500 000 years ago. Consistent with this recent divergence, we observed considerable variation in the genealogical patterns among loci. For some loci, all alleles within each species formed a monophyletic group, while at other loci, species were intermingled on the phylogeny of alleles. This intermingling probably reflects both incomplete lineage sorting and gene flow after divergence. Likelihood ratio tests rejected a strict allopatric model with no gene flow in comparisons between each pair of species. Gene flow was asymmetric: no gene flow was detected into M. domesticus, while significant gene flow was detected into both M. castaneus and M. musculus. Finally, most of the gene flow occurred at autosomal loci, resulting in a significantly higher ratio of fixed differences to polymorphisms at the X and Y chromosomes relative to autosomes in some comparisons, or just the X chromosome in others, emphasizing the important role of the sex chromosomes in general and the X chromosome in particular in speciation.

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 FST 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., This work was supported by Fundação para a Ciência e a Tecnologia (SFRH/BD/4621/2001 and SFRH/BPD/24743/2005 grants to A.G. and Research Projects POCTI/BSE/40280/2001 and PTDC/BIA-BDE/72304/2006), JCCM PA106-0170 and INIA FAU 2006-00014-C02-02, and by a National Science Foundation grant to M.W.N.

Summary We have sequenced 2388 bp of the European rabbit sex determining region Y (SRY) gene. These data provide a 10-fold increase in the coverage of the Y chromosome in this species, including the entire open reading frame of the SRY, the polyadenylation signal, and two repetitive sequences in the 5¢ -region. A survey of 2021 bp of this gene in eight domestic breeds and four wild individuals revealed a total of nine single nucleotide polymorphisms and one indel, defining two deeply divergent lineages. The resulting estimation of nucleotide diversity (p ¼ 1.34 · 10 )3 ) is very high when compared with other species, but no variability was detected among the domestic breeds. This study represents a first step in the characterization of the European rabbit Y chromosome and its variability. These sequences can be used in additional phylogeographical analyses of the European rabbit and other Leporid species, as well as in evolutionary studies of sex determination and the Y chromosome in wild species.

"Chantier qualité spécifique "Auteurs Externes" département de Génétique animale : uniquement liaison auteur au référentiel HR-Access "; International audience

In this work we report the genetic polymorphism of a 7-bp insertion in the 3' untranslated region of the rabbit SRY gene. The polymorphic GAATTAA motif was found exclusively in one of the two divergent rabbit Y-chromosomal lineages, suggesting that its origin is more recent than the separation of the O. c. algirus and O. c. cuniculus Y-chromosomes. In addition, the remarkable observation of haplotypes exhibiting 0, 1 and 2 7-bp inserts in essentially all algirus populations suggests that the rabbit SRY gene is duplicated and evolving under concerted evolution., This work was partly supported by Fundação para a Ciência e a Tecnologia (SFRH/BD/4621/2001 Ph.D. grant to Armando Geraldes and Research Project POCTI/BSE/40280/2001).