Single-cell and single-nucleus RNA-sequencing from paired normal-adenocarcinoma lung samples provide both common and discordant biological insights.

Whether single-cell RNA-sequencing (scRNA-seq) captures the same biological information as single-nucleus RNA-sequencing (snRNA-seq) remains uncertain and likely to be context-dependent. Herein, a head-to-head comparison was performed in matched normal-adenocarcinoma human lung samples to assess biological insights derived from scRNA-seq versus snRNA-seq and better understand the cellular transition that occurs from normal to tumoral tissue. Here, the transcriptome of 160,621 cells/nuclei was obtained. In non-tumor lung, cell type proportions varied widely between scRNA-seq and snRNA-seq with a predominance of immune cells in the former (81.5%) and epithelial cells (69.9%) in the later. Similar results were observed in adenocarcinomas, in addition to an overall increase in cell type heterogeneity and a greater prevalence of copy number variants in cells of epithelial origin, which suggests malignant assignment. The cell type transition that occurs from normal lung tissue to adenocarcinoma was not always concordant whether cells or nuclei were examined. As expected, large differential expression of the whole-cell and nuclear transcriptome was observed, but cell-type specific changes of paired normal and tumor lung samples revealed a set of common genes in the cells and nuclei involved in cancer-related pathways. In addition, we showed that the ligand-receptor interactome landscape of lung adenocarcinoma was largely different whether cells or nuclei were evaluated. Immune cell depletion in fresh specimens partly mitigated the difference in cell type composition observed between cells and nuclei. However, the extra manipulations affected cell viability and amplified the transcriptional signatures associated with stress responses. In conclusion, research applications focussing on mapping the immune landscape of lung adenocarcinoma benefit from scRNA-seq in fresh samples, whereas snRNA-seq of frozen samples provide a low-cost alternative to profile more epithelial and cancer cells, and yield cell type proportions that more closely match tissue content.

A Test to Comprehensively Capture the Known Genetic Component of Familial Pulmonary Fibrosis.

The recent European Respiratory Society statement on familial pulmonary fibrosis supports the need for genetic testing in the care of patients and their relatives. However, no solution (i.e., a concrete test) was provided to implement genetic testing in daily practice. Herein, we tabulated and standardized the nomenclature of 128 genetic variants in 20 genes implicated in adult-onset pulmonary fibrosis. The objective was to develop a laboratory-developed test (LDT) on the basis of standard Sanger sequencing to capture all known familial pulmonary fibrosis-associated variants. Targeted DNA fragments were amplified using harmonized PCR conditions to perform the LDT in a single 96-well plate. The new genetic test was evaluated in 62 sporadic cases of idiopathic pulmonary fibrosis. As expected in this population, we observed a low yield of disease-causing mutations. More important, 100% of targeted variants by the LDT were successfully evaluated. Furthermore, four variants of uncertain significance with in silico-predicted deleterious scores were identified in three patients, suggesting novel pathogenic variants in genes known to cause idiopathic pulmonary fibrosis. Finally, the MUC5B promoter variant rs35705950 was strongly enriched in these patients with a minor allele frequency of 41.1% compared with 10.6% in a matched population-based cohort (n = 29,060), leading to an estimation that this variant may explain up to 35% of the population-attributable risk. This LDT provides a solution for rapid clinical translation. Technical laboratory details are provided so that specialized pulmonary centers can implement the LDT in house to expedite the clinical recommendations of expert panels.

Inoculation with Rhizophagus Irregularis Does Not Alter Arbuscular Mycorrhizal Fungal Community Structure within the Roots of Corn, Wheat, and Soybean Crops.

Little is known about establishment success of the arbuscular mycorrhizal fungal (AMF) inocula and their effects on a soil-indigenous community of AMF. In this study, we assessed the effect of introducing Rhizophagus irregularis DAOM-197198 in soil under field condition on the community composition of indigenous AMF in the roots of corn (Zea mays), soybean (Glycine max), and wheat (Triticum aestivum). Three field trials were conducted with inoculated and non-inoculated plots. Four to ten roots and their rhizosphere soil samples of two growth stages for corn and wheat, and one growing stage of soybean, were collected, totalling 122 root and soil samples. Root colonization was measured microscopically, and the fungal communities were determined by paired-end Illumina MiSeq amplicon sequencing using 18S rDNA marker. After quality trimming and merging of paired ends, 6.7 million sequences could be assigned to 414 different operational taxonomic units. These could be assigned to 68 virtual taxa (VT) using the AMF reference sequence database MaarjAM. The most abundant VT corresponded to R. irregularis. The inoculation treatment did not influence the presence of R. irregularis, or AMF community diversity in roots. This seems to indicate that inoculation with R. irregularis DAOM-197198 does not change the indigenous AMF community composition, probably because it is already present in high abundance naturally.

Temporal variation in spatial genetic structure during population outbreaks: Distinguishing among different potential drivers of spatial synchrony.

Spatial synchrony is a common characteristic of spatio-temporal population dynamics across many taxa. While it is known that both dispersal and spatially autocorrelated environmental variation (i.e., the Moran effect) can synchronize populations, the relative contributions of each, and how they interact, are generally unknown. Distinguishing these mechanisms and their effects on synchrony can help us to better understand spatial population dynamics, design conservation and management strategies, and predict climate change impacts. Population genetic data can be used to tease apart these two processes as the spatio-temporal genetic patterns they create are expected to be different. A challenge, however, is that genetic data are often collected at a single point in time, which may introduce context-specific bias. Spatio-temporal sampling strategies can be used to reduce bias and to improve our characterization of the drivers of spatial synchrony. Using spatio-temporal analyses of genotypic data, our objective was to identify the relative support for these two mechanisms to the spatial synchrony in population dynamics of the irruptive forest insect pest, the spruce budworm (Choristoneura fumiferana), in Quebec (Canada). AMOVA, cluster analysis, isolation by distance, and sPCA were used to characterize spatio-temporal genomic variation using 1,370 SBW larvae sampled over four years (2012-2015) and genotyped at 3,562 SNP loci. We found evidence of overall weak spatial genetic structure that decreased from 2012 to 2015 and a genetic diversity homogenization among the sites. We also found genetic evidence of a long-distance dispersal event over >140 km. These results indicate that dispersal is the key mechanism involved in driving population synchrony of the outbreak. Early intervention management strategies that aim to control source populations have the potential to be effective through limiting dispersal. However, the timing of such interventions relative to outbreak progression is likely to influence their probability of success.

Putative Mitochondrial Sex Determination in the Bivalvia: Insights From a Hybrid Transcriptome Assembly in Freshwater Mussels.

Bivalves exhibit an astonishing diversity of sexual systems, with genetic and environmental determinants of sex, and possibly the only example of mitochondrial genes influencing sex determination pathways in animals. In contrast to all other animal species in which strict maternal inheritance (SMI) of mitochondria is the rule, bivalves possess a system known as doubly uniparental inheritance (DUI) of mitochondria in which maternal and paternal mitochondria (and their corresponding female-transmitted or F mtDNA and male-transmitted or M mtDNA genomes) are transmitted within a species. Species with DUI also possess sex-associated mtDNA-encoded proteins (in addition to the typical set of 13), which have been hypothesized to play a role in sex determination. In this study, we analyzed the sex-biased transcriptome in gonads of two closely-related freshwater mussel species with different reproductive and mitochondrial transmission modes: the gonochoric, DUI species, Utterbackia peninsularis, and the hermaphroditic, SMI species, Utterbackia imbecillis. Through comparative analysis with other DUI and non-DUI bivalve transcriptomes already available, we identify common male and female-specific genes, as well as SMI and DUI-related genes, that are probably involved in sex determination and mitochondrial inheritance in this animal group. Our results contribute to the understanding of what could be the first animal sex determination system involving the mitochondrial genome.

A commercial seaweed extract structured microbial communities associated with tomato and pepper roots and significantly increased crop yield.

Seaweeds have been used as a source of natural fertilizer and biostimulant in agriculture for centuries. However, their effects on soil and crop root microbiota remain unclear. Here, we used a commercially available Ascophyllum nodosum extract (ANE) to test its effect on bacterial and fungal communities of rhizospheric soils and roots of pepper and tomato plants in greenhouse trials. Two independent trials were conducted in a split-block design. We used amplicon sequencing targeting fungal ITS and bacterial 16S rRNA gene to determine microbial community structure changes. We find that productivity parameters of root, shoot and fruit biomass were positively and significantly influenced by the ANE amendment. In addition, a-diversity differed significantly between amended and control plants, but only in some of the experimental conditions. Species composition among sites (b-diversity) differed according to the amendment treatment in all four communities (fungal-root, fungal-soil, bacterial-root and bacterial-soil). Finally, we identified a number of candidate taxa most strongly correlated with crop yield increases. Further studies on isolation and characterization of these microbial taxa linked to the application of liquid seaweed extract may help to enhance crop yield in sustainable agro-ecosystems.

Genome Survey of the Freshwater Mussel Venustaconcha ellipsiformis (Bivalvia: Unionida) Using a Hybrid De Novo Assembly Approach.

Freshwater mussels (Bivalvia: Unionida) serve an important role as aquatic ecosystem engineers but are one of the most critically imperilled groups of animals. Here, we used a combination of sequencing strategies to assemble and annotate a draft genome of Venustaconcha ellipsiformis, which will serve as a valuable genomic resource given the ecological value and unique "doubly uniparental inheritance" mode of mitochondrial DNA transmission of freshwater mussels. The genome described here was obtained by combining high-coverage short reads (65× genome coverage of Illumina paired-end and 11× genome coverage of mate-pairs sequences) with low-coverage Pacific Biosciences long reads (0.3× genome coverage). Briefly, the final scaffold assembly accounted for a total size of 1.54 Gb (366,926 scaffolds, N50 = 6.5 kb, with 2.3% of "N" nucleotides), representing 86% of the predicted genome size of 1.80 Gb, while over one third of the genome (37.5%) consisted of repeated elements and >85% of the core eukaryotic genes were recovered. Given the repeated genetic bottlenecks of V. ellipsiformis populations as a result of glaciations events, heterozygosity was also found to be remarkably low (0.6%), in contrast to most other sequenced bivalve species. Finally, we reassembled the full mitochondrial genome and found six polymorphic sites with respect to the previously published reference. This resource opens the way to comparative genomics studies to identify genes related to the unique adaptations of freshwater mussels and their distinctive mitochondrial inheritance mechanism.

Deciphering the Link between Doubly Uniparental Inheritance of mtDNA and Sex Determination in Bivalves: Clues from Comparative Transcriptomics.

Bivalves exhibit an astonishing diversity of sexual systems and sex-determining mechanisms. They can be gonochoric, hermaphroditic or androgenetic, with both genetic and environmental factors known to determine or influence sex. One unique sex-determining system involving the mitochondrial genome has also been hypothesized to exist in bivalves with doubly uniparental inheritance (DUI) of mtDNA. However, the link between DUI and sex determination remains obscure. In this study, we performed a comparative gonad transcriptomics analysis for two DUI-possessing freshwater mussel species to better understand the mechanisms underlying sex determination and DUI in these bivalves. We used a BLAST reciprocal analysis to identify orthologs between Venustaconcha ellipsiformis and Utterbackia peninsularis and compared our results with previously published sex-specific bivalve transcriptomes to identify conserved sex-determining genes. We also compared our data with other DUI species to identify candidate genes possibly involved in the regulation of DUI. A total of ∼12,000 orthologous relationships were found, with 2,583 genes differentially expressed in both species. Among these genes, key sex-determining factors previously reported in vertebrates and in bivalves (e.g., Sry, Dmrt1, Foxl2) were identified, suggesting that some steps of the sex-determination pathway may be deeply conserved in metazoans. Our results also support the hypothesis that a modified ubiquitination mechanism could be responsible for the retention of the paternal mtDNA in male bivalves, and revealed that DNA methylation could also be involved in the regulation of DUI. Globally, our results suggest that sets of genes associated with sex determination and DUI are similar in distantly-related DUI species.

A multi-disciplinary perspective on emergent and future innovations in peer review.

Peer review of research articles is a core part of our scholarly communication system. In spite of its importance, the status and purpose of peer review is often contested. What is its role in our modern digital research and communications infrastructure? Does it perform to the high standards with which it is generally regarded? Studies of peer review have shown that it is prone to bias and abuse in numerous dimensions, frequently unreliable, and can fail to detect even fraudulent research. With the advent of web technologies, we are now witnessing a phase of innovation and experimentation in our approaches to peer review. These developments prompted us to examine emerging models of peer review from a range of disciplines and venues, and to ask how they might address some of the issues with our current systems of peer review. We examine the functionality of a range of social Web platforms, and compare these with the traits underlying a viable peer review system: quality control, quantified performance metrics as engagement incentives, and certification and reputation. Ideally, any new systems will demonstrate that they out-perform and reduce the biases of existing models as much as possible. We conclude that there is considerable scope for new peer review initiatives to be developed, each with their own potential issues and advantages. We also propose a novel hybrid platform model that could, at least partially, resolve many of the socio-technical issues associated with peer review, and potentially disrupt the entire scholarly communication system. Success for any such development relies on reaching a critical threshold of research community engagement with both the process and the platform, and therefore cannot be achieved without a significant change of incentives in research environments.

Plant Community and Nitrogen Deposition as Drivers of Alpha and Beta Diversities of Prokaryotes in Reconstructed Oil Sand Soils and Natural Boreal Forest Soils.

The Athabasca oil sand deposit is one of the largest single oil deposits in the world. Following surface mining, companies are required to restore soil-like profiles that can support the previous land capabilities. The objective of this study was to assess whether the soil prokaryotic alpha diversity (α-diversity) and ß-diversity in oil sand soils reconstructed 20 to 30 years previously and planted to one of three vegetation types (coniferous or deciduous trees and grassland) were similar to those found in natural boreal forest soils subject to wildfire disturbance. Prokaryotic α-diversity and ß-diversity were assessed using massively parallel sequencing of 16S rRNA genes. The ß-diversity, but not the α-diversity, differed between reconstructed and natural soils. Bacteria associated with an oligotrophic lifestyle were more abundant in natural forest soils, whereas bacteria associated with a copiotrophic lifestyle were more abundant in reconstructed soils. Ammonia-oxidizing archaea were most abundant in reconstructed soils planted with grasses. Plant species were the main factor influencing α-diversity in natural and in reconstructed soils. Nitrogen deposition, pH, and plant species were the main factors influencing the ß-diversity of the prokaryotic communities in natural and reconstructed soils. The results highlight the importance of nitrogen deposition and aboveground-belowground relationships in shaping soil microbial communities in natural and reconstructed soils.IMPORTANCE Covering over 800 km2, land disturbed by the exploitation of the oil sands in Canada has to be restored. Here, we take advantage of the proximity between these reconstructed ecosystems and the boreal forest surrounding the oil sand mining area to study soil microbial community structure and processes in both natural and nonnatural environments. By identifying key characteristics shaping the structure of soil microbial communities, this study improved our understanding of how vegetation, soil characteristics and microbial communities interact and drive soil functions.

Phytogeographic and genetic variation in Sorbus, a traditional antidiabetic medicine-adaptation in action in both a plant and a discipline.

Mountain ash (Sorbus decora and S. americana) is used by the Cree Nation of the James Bay region of Quebec (Eeyou Istchee) as traditional medicine. Its potential as an antidiabetic medicine is thought to vary across its geographical range, yet little is known about the factors that affect its antioxidant capacity. Here, we examined metabolite gene expression in relation to antioxidant activity, linking phytochemistry and medicinal potential. Samples of leaf and bark from S. decora and S. americana were collected from 20 populations at four different latitudes. Two genes known to produce antidiabetic substances, flavonol synthase and squalene synthase, were analyzed using quantitative real time PCR. Gene expression was significantly higher for flavonol synthase compared to squalene synthase and increased in the most Northern latitude. Corresponding differences observed in the antioxidant capacity of ethanolic extracts from the bark of Sorbus spp. confirm that plants at higher latitudes increase production of stress-induced secondary metabolites and support Aboriginal perceptions of their higher medicinal potential. Modern genetic techniques such as quantitative real time PCR offer unprecedented resolution to substantiate and scrutinise Aboriginal medicinal plant perception. Furthermore, it offers valuable insights into how environmental stress can trigger an adaptive response resulting in the accumulation of secondary metabolites with human medicinal properties.

History repeats itself: genomic divergence in copepods.

Press stop, erase everything from now till some arbitrary time in the past and start recording life as it evolves once again. Would you see the same tape of life playing itself over and over, or would a different story unfold every time? The late Steven Jay Gould called this experiment replaying the tape of life and argued that any replay of the tape would lead evolution down a pathway radically different from the road actually taken (Gould 1989). This thought experiment has puzzled evolutionary biologists for a long time: how repeatable are evolutionary events? And if history does indeed repeat itself, what are the factors that may help us predict the path taken? A powerful means to address these questions at a small evolutionary scale is to study closely related populations that have evolved independently, under similar environmental conditions. This is precisely what Pereira et al. (2016) set out to do using marine copepods Tigriopus californicus, and present their results in this issue of Molecular Ecology. They show that evolution can be repeatable and even partly predictable, at least at the molecular level. As expected from theory, patterns of divergence were shaped by natural selection. At the same time, strong genetic drift due to small population sizes also constrained evolution down a similar evolutionary road, and probably contributed to repeatable patterns of genomic divergence.

Assessing the reproducibility of discriminant function analyses.

Data are the foundation of empirical research, yet all too often the datasets underlying published papers are unavailable, incorrect, or poorly curated. This is a serious issue, because future researchers are then unable to validate published results or reuse data to explore new ideas and hypotheses. Even if data files are securely stored and accessible, they must also be accompanied by accurate labels and identifiers. To assess how often problems with metadata or data curation affect the reproducibility of published results, we attempted to reproduce Discriminant Function Analyses (DFAs) from the field of organismal biology. DFA is a commonly used statistical analysis that has changed little since its inception almost eight decades ago, and therefore provides an opportunity to test reproducibility among datasets of varying ages. Out of 100 papers we initially surveyed, fourteen were excluded because they did not present the common types of quantitative result from their DFA or gave insufficient details of their DFA. Of the remaining 86 datasets, there were 15 cases for which we were unable to confidently relate the dataset we received to the one used in the published analysis. The reasons ranged from incomprehensible or absent variable labels, the DFA being performed on an unspecified subset of the data, or the dataset we received being incomplete. We focused on reproducing three common summary statistics from DFAs: the percent variance explained, the percentage correctly assigned and the largest discriminant function coefficient. The reproducibility of the first two was fairly high (20 of 26, and 44 of 60 datasets, respectively), whereas our success rate with the discriminant function coefficients was lower (15 of 26 datasets). When considering all three summary statistics, we were able to completely reproduce 46 (65%) of 71 datasets. While our results show that a majority of studies are reproducible, they highlight the fact that many studies still are not the carefully curated research that the scientific community and public expects.

The Accumulation of Deleterious Mutations as a Consequence of Domestication and Improvement in Sunflowers and Other Compositae Crops.

For populations to maintain optimal fitness, harmful mutations must be efficiently purged from the genome. Yet, under circumstances that diminish the effectiveness of natural selection, such as the process of plant and animal domestication, deleterious mutations are predicted to accumulate. Here, we compared the load of deleterious mutations in 21 accessions from natural populations and 19 domesticated accessions of the common sunflower using whole-transcriptome single nucleotide polymorphism data. Although we find that genetic diversity has been greatly reduced during domestication, the remaining mutations were disproportionally biased toward nonsynonymous substitutions. Bioinformatically predicted deleterious mutations affecting protein function were especially strongly over-represented. We also identify similar patterns in two other domesticated species of the sunflower family (globe artichoke and cardoon), indicating that this phenomenon is not due to idiosyncrasies of sunflower domestication or the sunflower genome. Finally, we provide unequivocal evidence that deleterious mutations accumulate in low recombining regions of the genome, due to the reduced efficacy of purifying selection. These results represent a conundrum for crop improvement efforts. Although the elimination of harmful mutations should be a long-term goal of plant and animal breeding programs, it will be difficult to weed them out because of limited recombination.

Genomics of homoploid hybrid speciation: diversity and transcriptional activity of long terminal repeat retrotransposons in hybrid sunflowers.

Hybridization is thought to play an important role in plant evolution by introducing novel genetic combinations and promoting genome restructuring. However, surprisingly little is known about the impact of hybridization on transposable element (TE) proliferation and the genomic response to TE activity. In this paper, we first review the mechanisms by which homoploid hybrid species may arise in nature. We then present hybrid sunflowers as a case study to examine transcriptional activity of long terminal repeat retrotransposons in the annual sunflowers Helianthus annuus, Helianthus petiolaris and their homoploid hybrid derivatives (H. paradoxus, H. anomalus and H. deserticola) using high-throughput transcriptome sequencing technologies (RNAseq). Sampling homoploid hybrid sunflower taxa revealed abundant variation in TE transcript accumulation. In addition, genetic diversity for several candidate genes hypothesized to regulate TE activity was characterized. Specifically, we highlight one candidate chromatin remodelling factor gene with a direct role in repressing TE activity in a hybrid species. This paper shows that TE amplification in hybrid lineages is more idiosyncratic than previously believed and provides a first step towards identifying the mechanisms responsible for regulating and repressing TE expansions.

Chromosomal evolution and patterns of introgression in helianthus.

Knowledge of the nature and extent of karyotypic differences between species provides insight into the evolutionary history of the genomes in question and, in the case of closely related species, the potential for genetic exchange between taxa. We constructed high-density genetic maps of the silverleaf sunflower (Helianthus argophyllus) and Algodones Dune sunflower (H. niveus ssp. tephrodes) genomes and compared them to a consensus map of cultivated sunflower (H. annuus) to identify chromosomal rearrangements between species. The genetic maps of H. argophyllus and H. niveus ssp. tephrodes included 17 linkage groups each and spanned 1337 and 1478 cM, respectively. Comparative analyses revealed greater divergence between H. annuus and H. niveus ssp. tephrodes (13 inverted segments, 18 translocated segments) than between H. annuus and H. argophyllus (10 inverted segments, 8 translocated segments), consistent with their known phylogenetic relationships. Marker order was conserved across much of the genome, with 83 and 64% of the H. argophyllus and H. niveus ssp. tephrodes genomes, respectively, being syntenic with H. annuus. Population genomic analyses between H. annuus and H. argophyllus, which are sympatric across a portion of the natural range of H. annuus, revealed significantly elevated genetic structure in rearranged portions of the genome, indicating that such rearrangements are associated with restricted gene flow between these two species.

RNA-seq reveals transcriptomic shock involving transposable elements reactivation in hybrids of young lake whitefish species.

Identifying the molecular basis of reproductive isolation among diverging lineages represents an essential step toward understanding speciation in natural populations. Postzygotic barriers can lead to hybrid breakdown, a syndrome that has been documented in several systems, potentially involving the reactivation of transposable elements. In northeastern North America, two lake whitefish lineages have repeatedly colonized postglacial lakes ~12,000 years ago, and a dwarf limnetic species has evolved multiple times from the normal benthic species. Reproductive isolation is incomplete between them; viable hybrids can be generated in the laboratory but significant mortality occurs and is associated with a malformed phenotype in backcross embryos, thus revealing a hybrid breakdown syndrome. By means of RNA-seq analyses, the objective of this study was to determine which genes were misregulated in hybrids and rigorously test the hypothesis of transposable element reactivation. We compared the transcriptomic landscape in pure embryos, F1-hybrids, and healthy and malformed backcrosses at the late embryonic stage. Extensive expression differences consistent with previously documented adaptive divergence between pure normal and dwarf embryos were identified for the first time. Pronounced transcriptome-wide deregulation in malformed backcrosses was observed, with over 15% of transcripts differentially expressed in all comparisons, compared with 1.5% between pure parental forms. Convincing evidence of transposable elements and noncoding transcripts reactivation in malformed backcrosses is presented. We propose that hybrid breakdown likely results from extensive genomic incompatibilities, plausibly encompassing transposable elements. Combined with previous studies, these results reveal synergy among many reproductive barriers, thus maintaining divergence between these two young whitefish species.

Shared selective pressure and local genomic landscape lead to repeatable patterns of genomic divergence in sunflowers.

The repeated evolution of traits in organisms facing similar environmental conditions is considered to be fundamental evidence for the role of natural selection in moulding phenotypes. Yet, aside from case studies of parallel evolution and its genetic basis, the repeatability of evolution at the level of the whole genome remains poorly characterized. Here, through the use of transcriptome sequencing, we examined genomic divergence for three pairs of sister species of sunflowers. Two of the pairs (Helianthus petiolaris - H. debilis and H. annuus - H. argophyllus) have diverged along a similar latitudinal gradient and presumably experienced similar selective pressure. In contrast, a third species pair (H. exilis - H. bolanderi) diverged along a longitudinal gradient. Analyses of divergence, as measured in terms of FST, indicated little repeatability across the three pairs of species for individual genetic markers (SNPs), modest repeatability at the level of individual genes and the highest repeatability when large regions of the genome were compared. As expected, higher repeatability was observed for the two species pairs that have diverged along a similar latitudinal gradient, with genes involved in flowering time among the most divergent genes. Genes showing extreme low or high differentiation were more similar than genes showing medium levels of divergence, implying that both purifying and divergent selection contributed to repeatable patterns of divergence. The location of a gene along the chromosome also predicted divergence levels, presumably because of shared heterogeneity in both recombination and mutation rates. In conclusion, repeated genome evolution appeared to result from both similar selective pressures and shared local genomic landscapes.