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.

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.

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.

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.

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.

Advances in Ecological Speciation: an integrative approach.

The role of natural selection in promoting reproductive isolation has received substantial renewed interest within the last two decades. As a consequence, the study of ecological speciation has become an extremely productive research area in modern evolutionary biology. Recent innovations in sequencing technologies offer an unprecedented opportunity to study the mechanisms involved in ecological speciation. Genome scans provide significant insights but have some important limitations; efforts are needed to integrate them with other approaches to make full use of the sequencing data deluge. An international conference 'Advances in Ecological Speciation' organized by the University of Porto (Portugal) aimed to review current progress in ecological speciation. Using some of the examples presented at the conference, we highlight the benefits of integrating ecological and genomic data and discuss different mechanisms of parallel evolution. Finally, future avenues of research are suggested to advance our knowledge concerning the role of natural selection in the establishment of reproductive isolation during ecological speciation.

Mandated data archiving greatly improves access to research data.

The data underlying scientific papers should be accessible to researchers both now and in the future, but how best can we ensure that these data are available? Here we examine the effectiveness of four approaches to data archiving: no stated archiving policy, recommending (but not requiring) archiving, and two versions of mandating data deposition at acceptance. We control for differences between data types by trying to obtain data from papers that use a single, widespread population genetic analysis, structure. At one extreme, we found that mandated data archiving policies that require the inclusion of a data availability statement in the manuscript improve the odds of finding the data online almost 1000-fold compared to having no policy. However, archiving rates at journals with less stringent policies were only very slightly higher than those with no policy at all. We also assessed the effectiveness of asking for data directly from authors and obtained over half of the requested datasets, albeit with ∼8 d delay and some disagreement with authors. Given the long-term benefits of data accessibility to the academic community, we believe that journal-based mandatory data archiving policies and mandatory data availability statements should be more widely adopted.

Polygenic inheritance and its interplay with smoking history in predicting lung cancer diagnosis: a French-Canadian case-control cohort.

BACKGROUND: The most near-term clinical application of genome-wide association studies in lung cancer is a polygenic risk score (PRS). METHODS: A case-control dataset was generated consisting of 4002 lung cancer cases from the LORD project and 20,010 ethnically matched controls from CARTaGENE. A genome-wide PRS including >1.1 million genetic variants was derived and validated in UK Biobank (n = 5419 lung cancer cases). The predictive ability and diagnostic discrimination performance of the PRS was tested in LORD/CARTaGENE and benchmarked against previous PRSs from the literature. Stratified analyses were performed by smoking status and genetic risk groups defined as low (<20th percentile), intermediate (20-80th percentile) and high (>80th percentile) PRS. FINDINGS: The phenotypic variance explained and the effect size of the genome-wide PRS numerically outperformed previous PRSs. Individuals with high genetic risk had a 2-fold odds of lung cancer compared to low genetic risk. The PRS was an independent predictor of lung cancer beyond conventional clinical risk factors, but its diagnostic discrimination performance was incremental in an integrated risk model. Smoking increased the odds of lung cancer by 7.7-fold in low genetic risk and by 11.3-fold in high genetic risk. Smoking with high genetic risk was associated with a 17-fold increase in the odds of lung cancer compared to individuals who never smoked and with low genetic risk. INTERPRETATION: Individuals at low genetic risk are not protected against the smoking-related risk of lung cancer. The joint multiplicative effect of PRS and smoking increases the odds of lung cancer by nearly 20-fold. FUNDING: This work was supported by the CQDM and the IUCPQ Foundation owing to a generous donation from Mr. Normand Lord.

A new approach to quantify the adaptive potential of gene expression variation in gymnosperms.

Variation in patterns of gene expression contributes to phenotypic diversity and can ultimately predict adaptive responses. However, in many cases, the consequences of regulatory mutations on patterns of gene expression and ultimately phenotypic differences remain elusive. A standard way to study the genetic architecture of expression variation in model systems has been to map gene expression variation to genetic loci (Fig. 1a). At the same time, in many nonmodel species, especially for long-lived organisms, controlled crosses are not feasible. If we are to expand our understanding of the role of regulatory mutations on phenotypes, we need to develop new methodologies to study species under ecologically relevant conditions. In this issue of Molecular Ecology, Verta et al. (2013) present a new approach to analyse gene expression variation and regulatory networks in gymnosperms (Fig. 1b). They capitalized on the fact that gymnosperm seeds contain an energy storage tissue (the megagametophyte) that is directly derived from a single haploid cell (the megaspore). The authors identified over 800 genes for which expression segregated in this maternally inherited haploid tissue. Based on the observed segregation patterns, these genes (Mendelian Expression Traits) are most probably controlled by biallelic variants at a single locus. Most of these genes also belonged to different regulatory networks, except for one large group of 180 genes under the control of a putative trans-acting factor. In addition, the approach developed here may also help to uncover the effect of rare recessive mutations, which usually remain hidden in a heterozygous state in diploid individuals. The appeal of the work by Verta et al. (2013) to study gene expression variation is in its simplicity, which circumvents several of the hurdles behind traditional expression quantitative trait locus (eQTL) studies, and could potentially be applied to a large number of species.

Targeted sequence capture and resequencing implies a predominant role of regulatory regions in the divergence of a sympatric lake whitefish species pair (Coregonus clupeaformis).

Latest technological developments in evolutionary biology bring new challenges in documenting the intricate genetic architecture of species in the process of divergence. Sympatric populations of lake whitefish represent one of the key systems to investigate this issue. Despite the value of random genotype-by-sequencing methods and decreasing cost of sequencing technologies, it remains challenging to investigate variation in coding regions, especially in the case of recently duplicated genomes as in salmonids, as this greatly complicates whole genome resequencing. We thus designed a sequence capture array targeting 2773 annotated genes to document the nature and the extent of genomic divergence between sympatric dwarf and normal whitefish. Among the 2728 genes successfully captured, a total of 2182 coding and 10,415 noncoding putative single-nucleotide polymorphisms (SNPs) were identified after applying a first set of basic filters. A genome scan with a quality-refined selection of 2203 SNPs identified 267 outlier SNPs in 210 candidate genes located in genomic regions potentially involved in whitefish divergence and reproductive isolation. We found highly heterogeneous FST estimates among SNP loci. There was an overall low level of coding polymorphism, with a predominance of noncoding mutations among outliers. The heterogeneous patterns of divergence among loci confirm the porous nature of genomes during speciation with gene flow. Considering that few protein-coding mutations were identified as highly divergent, our results, along with previous transcriptomic studies, imply that changes in regulatory regions most likely had a greater role in the process of whitefish population divergence than protein-coding mutations. This study is the first to demonstrate the efficiency of large-scale targeted resequencing for a nonmodel species with such a large and unsequenced genome.

Recommendations for utilizing and reporting population genetic analyses: the reproducibility of genetic clustering using the program STRUCTURE.

Reproducibility is the benchmark for results and conclusions drawn from scientific studies, but systematic studies on the reproducibility of scientific results are surprisingly rare. Moreover, many modern statistical methods make use of 'random walk' model fitting procedures, and these are inherently stochastic in their output. Does the combination of these statistical procedures and current standards of data archiving and method reporting permit the reproduction of the authors' results? To test this, we reanalysed data sets gathered from papers using the software package STRUCTURE to identify genetically similar clusters of individuals. We find that reproducing structure results can be difficult despite the straightforward requirements of the program. Our results indicate that 30% of analyses were unable to reproduce the same number of population clusters. To improve this, we make recommendations for future use of the software and for reporting STRUCTURE analyses and results in published works.

Contemporary hybrid speciation in sculpins (Cottus spp.).

Natural hybridization between closely related taxa is frequent in many organismal groups, yet it has long been perceived as a force preventing diversification and speciation, especially so in animals. In recent years, growing evidence in favour of hybridization facilitating adaptive divergence has accumulated (Mallet 2007; Mavárez & Linares 2008; Nolte & Tautz 2010). Homoploid hybrid speciation (the formation of hybrid lineages without changes in chromosome number) occurs when distinct species come into contact, hybridize, and at least in part of their range, produce hybrid swarms. If the hybrid genotypes can then colonize areas of the adaptive landscape inaccessible to ancestral species, they may eventually form new distinct lineages, reproductively isolated from their ancestors. Invasive sculpins (Cottus sp.) are one of a few good examples of homoploid hybrid speciation in animals. In this issue, Stemshorn et al. (2011) identified three distinct hybrid lineages, which have emerged out of a secondary contact situation of Cottus rhenanus and Cottus perifretum. Hybrids have recently invaded large river habitats unsuitable to ancestral species. Through the use of genetic mapping, the authors established that contrary to expectations, chromosomal rearrangements were not apparent in the hybrid lineages. In addition, different population genetic models were tested and the results suggest that contemporary gene flow from ancestral species represents an important component of the system. As such, recent and ongoing hybridization appears to be promoting the appearance of phenotypes adapted to novel environments. The examination of partially isolated lineages such as invasive hybrid sculpins should permit to identify early adaptive genetic changes before they become confounded by differences arising once speciation is complete.

SNP signatures of selection on standing genetic variation and their association with adaptive phenotypes along gradients of ecological speciation in lake whitefish species pairs (Coregonus spp.).

As populations adapt to novel environments, divergent selection will promote heterogeneous genomic differentiation via reductions in gene flow for loci underlying adaptive traits. Using a data set of over 100 SNP markers, genome scans were performed to investigate the effect of natural selection maintaining differentiation in five lakes harbouring sympatric pairs of normal and dwarf lake whitefish (Coregonus clupeaformis). A variable proportion of SNPs (between 0% and 12%) was identified as outliers, which corroborated the predicted intensity of competitive interactions unique to each lake. Moreover, strong reduction in heterozygosity was typically observed for outlier loci in dwarf but not in normal whitefish, indicating that directional selection has been acting on standing genetic variation more intensively in dwarf whitefish. SNP associations in backcross hybrid progeny identified 16 genes exhibiting genotype-phenotype associations for four adaptive traits (growth, swimming activity, gill rakers and condition factor). However, neither simple relationship between elevated levels of genetic differentiation with adaptive phenotype nor conspicuous genetic signatures for parallelism at outlier loci were detected, which underscores the importance of independent evolution among lakes. The integration of phenotypic, transcriptomic and functional genomic information identified two candidate genes (sodium potassium ATPase and triosephosphate isomerase) involved in the recent ecological divergence of lake whitefish. Finally, the identification of several markers under divergent selection suggests that many genes, in an environment-specific manner, are recruited by selection and ultimately contributed to the repeated ecological speciation of a dwarf phenotype.

Mining transcriptome sequences towards identifying adaptive single nucleotide polymorphisms in lake whitefish species pairs (Coregonus spp. Salmonidae).

Next-generation sequencing allows the discovery of large numbers of single nucleotide polymorphisms (SNPs) in species where little genomic information was previously available. Here, we assembled, de novo, over 130 mb of non-normalized cDNA using 454 pyrosequencing data from dwarf and normal lake whitefish and backcross hybrids. Our main goals were to gather a large data set of SNP markers, document their distribution within coding regions, evaluate the effect of species divergence on allele frequencies and combine results with previous genomic studies to identify candidate genes underlying the adaptive divergence of lake whitefish. We identified 6094 putative SNPs in 2674 contigs (mean size: 576 bp, range: 101-6116) and 1540 synonymous and 1734 non-synonymous mutations for a genome-wide non-synonymous to synonymous substitution rate ratio (p(N)/p(S)) of 0.37. As expected based on the young age (<15 000 years) of whitefish species pair, the overall level of divergence between them was relatively weak. Yet, 89 SNPs showed pronounced allele frequency differences between sympatric normal and dwarf whitefish. Among these, SNPs in genes annotated to energy metabolic functions were the most abundant and this, in addition to previous experimental data at the gene expression and phenotypic level, brings compelling evidence that genes involved in energy metabolism are prime candidates explaining the adaptive divergence of lake whitefish species pairs. Finally, we unexpectedly identified 44 contigs annotated to transposable elements and these were predominantly composed of backcross hybrids sequences. This indicates an elevated activity of transposable elements, which could potentially contribute to the reduced fitness of hybrids previously documented.

The transcriptomics of sympatric dwarf and normal lake whitefish (Coregonus clupeaformis spp., Salmonidae) divergence as revealed by next-generation sequencing.

Gene expression divergence is one of the mechanisms thought to be involved in the emergence of incipient species. Next-generation sequencing has become an extremely valuable tool for the study of this process by allowing whole transcriptome sequencing, or RNA-Seq. We have conducted a 454 GS-FLX pyrosequencing experiment to refine our understanding of adaptive divergence between dwarf and normal lake whitefish species (Coregonus clupeaformis spp.). The objectives were to: (i) investigate transcriptomic divergence as measured by liver RNA-Seq; (ii) test the correlation between divergence in expression and sequence polymorphism; and (iii) investigate the extent of allelic imbalance. We also compared the results of RNA-seq with those of a previous microarray study performed on the same fish. Following de novo assembly, results showed that normal whitefish overexpressed more contigs associated with protein synthesis while dwarf fish overexpressed more contigs related to energy metabolism, immunity and DNA replication and repair. Moreover, 63 SNPs showed significant allelic imbalance, and this phenomenon prevailed in the recently diverged dwarf whitefish. Results also showed an absence of correlation between gene expression divergence as measured by RNA-Seq and either polymorphism rate or sequence divergence between normal and dwarf whitefish. This study reiterates an important role for gene expression divergence, and provides evidence for allele-specific expression divergence as well as evolutionary decoupling of regulatory and coding sequences in the adaptive divergence of normal and dwarf whitefish. It also demonstrates how next-generation sequencing can lead to a more comprehensive understanding of transcriptomic divergence in a young species pair.

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.

Divergence in gene regulation at young life history stages of whitefish (Coregonus sp.) and the emergence of genomic isolation.

BACKGROUND: The evolution of barriers to reproduction is of key interest to understand speciation. However, there may be a current bias towards studying intrinsic postzygotic isolation in old species pairs as compared to the emergence of barriers to gene flow through adaptive divergence. This study evaluates the relative importance of both processes in the evolution of genomic isolation in incipient species of whitefish (Coregonus clupeaformis) for which preliminary data suggest that postzygotic isolation emerges with intrinsic factors acting at embryo stages but also due to extrinsic factors during adult life. RESULTS: Gene expression data were screened using cDNA microarrays to identify regulatory changes at embryo and juvenile stages that provide evidence for genomic divergence at the underlying genetic factors. A comparison of different life history stages shows that 16-week old juvenile fish have 14 times more genes displaying significant regulatory divergence than embryos. Furthermore, regulatory changes in juvenile fish match patterns in adult fish suggesting that gene expression divergence is established early in juvenile fish and persists throughout the adult phase. Comparative analyses with results from previous studies on dwarf-normal species pairs show that at least 26 genetic factors identified in juvenile fish are candidate traits for adaptive divergence in adult fish. Eight of these show parallel directions of gene expression divergence independent of tissue type or age of the fish. The latter are associated with energy metabolism, a complex trait known to drive adaptive divergence in dwarf and normal whitefish. CONCLUSION: Although experimental evidence suggests the existence of genetic factors that cause intrinsic postzygotic isolation acting in embryos, the analysis presented here provided few candidate genes in embryos, which also corroborate previous studies showing a lack of ecological divergence between sympatric dwarf and normal whitefish at the larval stage. In contrast, gene expression divergence in juveniles can be linked to adaptive traits and seems to be driven by positive selection. The results support the idea that adaptive differentiation may be more important in explaining the emergence of barriers to gene flow in an early phase of speciation by providing a broad genomic basis for extrinsic postzygotic isolation rather than intrinsic barriers.