Outlier loci highlight the direction of introgression in oaks

Loci considered to be under selection are generally avoided in attempts to infer past demographic processes as they do not fit neutral model assumptions. However, opportunities to better reconstruct some aspects of past demography might thus be missed. Here we examined genetic differentiation between two sympatric European oak species with contrasting ecological dynamics (Quercus robur and Quercus petraea) with both outlier (i.e. loci possibly affected by divergent selection between species or by hitchhiking effects with genomic regions under selection) and nonoutlier loci. We sampled 855 individuals in six mixed forests in France and genotyped them with a set of 262 SNPs enriched with markers showing high interspecific differentiation, resulting in accurate species delimitation. We identified between 13 and 74 interspecific outlier loci, depending on the coalescent simulation models and parameters used. Greater genetic diversity was predicted in Q. petraea (a late‐successional species) than in Q. robur (an early successional species) as introgression should theoretically occur predominantly from the resident species to the invading species. Remarkably, this prediction was verified with outlier loci but not with nonoutlier loci. We suggest that the lower effective interspecific gene flow at loci showing high interspecific divergence has better preserved the signal of past asymmetric introgression towards Q. petraea caused by the species' contrasting dynamics. Using markers under selection to reconstruct past demographic processes could therefore have broader potential than generally recognized.     

Development and Characterization of Genomic and Gene-Based Microsatellite Markers in North American Red Oak Species

Oaks (Quercus: Fagaceae) are ecological and economic keystones of many forested ecosystems but effective genetic management strategies are hindered by high levels of phenotypic plasticity within species and frequent hybridization among them. These same features, however, make oak communities suited for the study of speciation, hybridization, and genetic adaptation. Efforts to develop new and to adapt existing genomic resources to less-studied members of this genus should not only improve oak conservation and management but also aid the study of fundamental evolutionary processes. Here, we present a suite of 27 highly polymorphic simple sequence repeat (SSR) markers tested in four North American red oak (Quercus section Lobatae) species: Q. rubra, Q. ellipsoidalis, Q. coccinea, and Q. velutina. Five markers are genomic SSRs (gSSRs) — four novel and one previously transferred from Q. petraea — and 22 are gene-based SSRs derived from Q. robur and Q. petraea expressed sequence tags (EST-SSRs). Overall, levels of polymorphism detected with these primer pairs were high, with gene diversity (H _e) averaging 0.66 across all loci in natural populations. In addition, we show that EST-SSR markers may have the potential to detect divergent selection at stress-resistance candidate genes among closely related oak species.     

Heterogeneity and concordance in locus‐specific differentiation and introgression between species of towhees

The maintenance or breakdown of reproductive isolation is an observable outcome of secondary contact between species. In cases where hybrids beyond the F1 are formed, the representation of each species' ancestry can vary dramatically among genomic regions. This genomic heterogeneity in ancestry and introgression can offer insight into evolutionary processes, particularly if introgression is compared in multiple hybrid zones. Similarly, considerable heterogeneity exists across the genome in the extent to which populations and species have diverged, reflecting the combined effects of different evolutionary processes on genetic variation. We studied hybridization across two hybrid zones of two phenotypically well‐differentiated bird species in Mexico (Pipilo maculatus and P. ocai), to investigate genomic heterogeneity in differentiation and introgression. Using genotyping‐by‐sequencing (GBS) and hierarchical Bayesian models, we genotyped 460 birds at over 41 000 single nucleotide polymorphism (SNP) loci. We identified loci exhibiting extreme introgression relative to the genome‐wide expectation using a Bayesian genomic cline model. We also estimated locus‐specific F_ST and identified loci with exceptionally high genetic divergence between the parental species. We found some concordance of locus‐specific introgression in the two independent hybrid zones (6–20% of extreme loci shared across zones), reflecting areas of the genome that experience similar gene flow when the species interact. Additionally, heterogeneity in introgression and divergence across the genome revealed another subset of loci under the influence of locally specific factors. These results are consistent with a history in which reproductive isolation has been influenced by a common set of loci in both hybrid zones, but where local environmental and stochastic factors also lead to genomic differentiation.     

Differential introgression from a sister species explains high FST outlier loci within a mussel species

Scanning genomes for loci with high levels of population differentiation has become a standard of population genetics. F_ST outlier loci are most often interpreted as signatures of local selection, but outliers might arise for many other reasons too often left unexplored. Here, we tried to identify further the history and genetic basis underlying strong differentiation at F_ST outlier loci in a marine mussel. A genome scan of genetic differentiation has been conducted between Atlantic and Mediterranean populations of Mytilus galloprovincialis. The differentiation was low overall (F_ST = 0.03), but seven loci (2%) were strong F_ST outliers. We then analysed DNA sequence polymorphism at two outlier loci. The genetic structure proved to be the consequence of differential introgression of alleles from the sister‐hybridizing species Mytilus edulis. Surprisingly, the Mediterranean population was the most introgressed at these two loci, although the contact zone between the two species is nowadays localized along the Atlantic coasts of France and the British Isles. A historical contact between M. edulis and Mediterranean M. galloprovincialis should have happened during glacial periods. It proved difficult to disentangle two hypotheses: (i) introgression was adaptive, implying edulis alleles have been favoured in Mediterranean populations, or (ii) the genetic architecture of the barrier to edulis gene flow is different between the two M. galloprovincialis backgrounds. Five of the seven outliers between M. galloprovincialis populations were also outliers between M. edulis and Atlantic M. galloprovincialis, which would support the latter hypothesis. Differential introgression across semi‐permeable barriers to gene flow is a neglected scenario to interpret outlying loci that may prove more widespread than anticipated.     

Gene flow,divergent selection and resistance to introgression in two species of morning glories (Ipomoea)

Gene flow is thought to impede genetic divergence and speciation by homogenizing genomes. Recent theory and research suggest that sufficiently strong divergent selection can overpower gene flow, leading to loci that are highly differentiated compared to others. However, there are also alternative explanations for this pattern. Independent evidence that loci in highly differentiated regions are under divergent selection would allow these explanations to be distinguished, but such evidence is scarce. Here, we present multiple lines of evidence that many of the highly divergent SNPs in a pair of sister morning glory species, Ipomoea cordatotriloba and I. lacunosa, are the result of divergent selection in the face of gene flow. We analysed a SNP data set across the genome to assess the amount of gene flow, resistance to introgression and patterns of selection on loci resistant to introgression. We show that differentiation between the two species is much lower in sympatry than in allopatry, consistent with interspecific gene flow in sympatry. Gene flow appears to be substantially greater from I. lacunosa to I. cordatotriloba than in the reverse direction, resulting in sympatric and allopatric I. cordatotriloba being substantially more different than sympatric and allopatric I. lacunosa. Many SNPs highly differentiated in allopatry have experienced divergent selection, and, despite gene flow in sympatry, resist homogenization in sympatry. Finally, five out of eight floral and inflorescence characteristics measured exhibit asymmetric convergence in sympatry. Consistent with the pattern of gene flow, I. cordatotriloba traits become much more like those of I. lacunosa than the reverse. Our investigation reveals the complex interplay between selection and gene flow that can occur during the early stages of speciation.     

High degree of genetic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene‐based genome‐scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection—as determined with several outlier detection methods—was low (F_ST = 0.033 or 0.011, respectively), whereas average F_ST for directionally selected markers was significantly higher (F_ST = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.     

Lack of gene flow: Narrow and dispersed differentiation islands in a triplet of Leptidea butterfly species

Genome scans in recently separated species can inform on molecular mechanisms and evolutionary processes driving divergence. Large‐scale polymorphism data from multiple species pairs are also key to investigate the repeatability of divergence—whether radiations tend to show parallel responses to similar selection pressures and/or underlying molecular forces. Here, we used whole‐genome resequencing data from six wood white (Leptidea sp.) butterfly populations, representing three closely related species with karyomorph variation, to infer the species' demographic history and characterize patterns of genomic diversity and differentiation. The analyses supported previously established species relationships, and there was no evidence for postdivergence gene flow. We identified significant intraspecific genetic structure, in particular between karyomorph extremes in the wood white (L. sinapis)—a species with a remarkable chromosome number cline across the distribution range. The genomic landscapes of differentiation were erratic, and outlier regions were narrow and dispersed. Highly differentiated (F_ST) regions generally had low genetic diversity (θ_π), but increased absolute divergence (D_XY) and excess of rare frequency variants (low Tajima's D). A minority of differentiation peaks were shared across species and population comparisons. However, highly differentiated regions contained genes with overrepresented functions related to metabolism, response to stimulus and cellular processes, indicating recurrent directional selection on a specific set of traits in all comparisons. In contrast to the majority of genome scans in recently diverged lineages, our data suggest that divergence landscapes in Leptidea have been shaped by directional selection and genetic drift rather than stable recombination landscapes and/or introgression.     

Adaptive divergence despite strong genetic drift: genomic analysis of the evolutionary mechanisms causing genetic differentiation in the island fox (Urocyon littoralis)

The evolutionary mechanisms generating the tremendous biodiversity of islands have long fascinated evolutionary biologists. Genetic drift and divergent selection are predicted to be strong on islands and both could drive population divergence and speciation. Alternatively, strong genetic drift may preclude adaptation. We conducted a genomic analysis to test the roles of genetic drift and divergent selection in causing genetic differentiation among populations of the island fox (Urocyon littoralis). This species consists of six subspecies, each of which occupies a different California Channel Island. Analysis of 5293 SNP loci generated using Restriction‐site Associated DNA (RAD) sequencing found support for genetic drift as the dominant evolutionary mechanism driving population divergence among island fox populations. In particular, populations had exceptionally low genetic variation, small N_e (range = 2.1–89.7; median = 19.4), and significant genetic signatures of bottlenecks. Moreover, islands with the lowest genetic variation (and, by inference, the strongest historical genetic drift) were most genetically differentiated from mainland grey foxes, and vice versa, indicating genetic drift drives genome‐wide divergence. Nonetheless, outlier tests identified 3.6–6.6% of loci as high F_ST outliers, suggesting that despite strong genetic drift, divergent selection contributes to population divergence. Patterns of similarity among populations based on high F_ST outliers mirrored patterns based on morphology, providing additional evidence that outliers reflect adaptive divergence. Extremely low genetic variation and small N_e in some island fox populations, particularly on San Nicolas Island, suggest that they may be vulnerable to fixation of deleterious alleles, decreased fitness and reduced adaptive potential.     

Bayesian estimation of genomic clines

We developed a Bayesian genomic cline model to study the genetic architecture of adaptive divergence and reproductive isolation between hybridizing lineages. This model quantifies locus‐specific patterns of introgression with two cline parameters that describe the probability of locus‐specific ancestry as a function of genome‐wide admixture. ‘Outlier’ loci with extreme patterns of introgression relative to most of the genome can be identified. These loci are potentially associated with adaptive divergence or reproductive isolation. We simulated genetic data for admixed populations that included neutral introgression, as well as loci that were subject to directional, epistatic or underdominant selection, and analysed these data using the Bayesian genomic cline model. Under many demographic conditions, underdominance or directional selection had detectable and predictable effects on cline parameters, and ‘outlier’ loci were greatly enriched for genetic regions affected by selection. We also analysed previously published genetic data from two transects through a hybrid zone between Mus domesticus and M. musculus. We found considerable variation in rates of introgression across the genome and particularly low rates of introgression for two X‐linked markers. There were similarities and differences in patterns of introgression between the two transects, which likely reflects a combination of stochastic variability because of genetic drift and geographic variation in the genetic architecture of reproductive isolation. By providing a robust framework to quantify and compare patterns of introgression among genetic regions and populations, the Bayesian genomic cline model will advance our understanding of the genetics of reproductive isolation and the speciation process.     

Introgression and selection shaped the evolutionary history of sympatric sister‐species of coral reef fishes (genus: Haemulon)

Closely related marine species with large overlapping ranges provide opportunities to study mechanisms of speciation, particularly when there is evidence of gene flow between such lineages. Here, we focus on a case of hybridization between the sympatric sister‐species Haemulon maculicauda and H. flaviguttatum, using Sanger sequencing of mitochondrial and nuclear loci, as well as 2422 single nucleotide polymorphisms (SNPs) obtained via restriction site‐associated DNA sequencing (RADSeq). Mitochondrial markers revealed a shared haplotype for COI and low divergence for CytB and CR between the sister‐species. On the other hand, complete lineage sorting was observed at the nuclear loci and most of the SNPs. Under neutral expectations, the smaller effective population size of mtDNA should lead to fixation of mutations faster than nDNA. Thus, these results suggest that hybridization in the recent past (0.174–0.263 Ma) led to introgression of the mtDNA, with little effect on the nuclear genome. Analyses of the SNP data revealed 28 loci potentially under divergent selection between the two species. The combination of mtDNA introgression and limited nuclear DNA introgression provides a mechanism for the evolution of independent lineages despite recurrent hybridization events. This study adds to the growing body of research that exemplifies how genetic divergence can be maintained in the presence of gene flow between closely related species.     

Phylogeography and lineage‐specific patterns of genetic diversity and molecular evolution in a group of North American skinks

Geography influences the evolutionary trajectory of species by mediating opportunities for hybridization, gene flow, demographic shifts and adaptation. We sought to understand how geography and introgression can generate species‐specific patterns of genetic diversity by examining phylogeographical relationships in the North American skink species Plestiodon multivirgatus and P. tetragrammus (Squamata: Scincidae). Using a multilocus dataset (three mitochondrial genes, four nuclear genes; a total of 3455 bp) we discovered mito‐nuclear discordance, consistent with mtDNA introgression. We further tested for evidence of species‐wide mtDNA introgression by using comparisons of genetic diversity, selection tests and extended Bayesian skyline analyses. Our findings suggest that P. multivirgatus acquired its mitochondrial genome from P. tetragrammus after their initial divergence. This putative species‐wide mitochondrial capture was further evidenced by statistically indistinguishable substitution rates between mtDNA and nDNA in P. multivirgatus. This rate discrepancy was observed in P. multivirgatus but not P. tetragrammus, which has important implications for studies that combine mtDNA and nDNA sequences when inferring time since divergence between taxa. Our findings suggest that by facilitating opportunities for interspecific introgression, geography can alter the course of molecular evolution between recently diverged lineages.     

Speciation genomics and a role for the Z chromosome in the early stages of divergence between Mexican ducks and mallards

Speciation is a continuous and dynamic process, and studying organisms during the early stages of this process can aid in identifying speciation mechanisms. The mallard (Anas platyrhynchos) and Mexican duck (A. [p.] diazi) are two recently diverged taxa with a history of hybridization and controversial taxonomy. To understand their evolutionary history, we conducted genomic scans to characterize patterns of genetic diversity and divergence across the mitochondrial DNA (mtDNA) control region, 3523 autosomal loci and 172 Z‐linked sex chromosome loci. Between the two taxa, Z‐linked loci (Φ_ST = 0.088) were 5.2 times more differentiated than autosomal DNA (Φ_ST = 0.017) but comparable to mtDNA (Φ_ST = 0.092). This elevated Z differentiation deviated from neutral expectations inferred from simulated data that incorporated demographic history and differences in effective population sizes between marker types. Furthermore, 3% of Z‐linked loci, compared to <0.1% of autosomal loci, were detected as outlier loci under divergent selection with elevated relative (Φ_ST) and absolute (d_XY) estimates of divergence. In contrast, the ratio of Z‐linked and autosomal differentiation among the seven Mexican duck sampling locations was close to 1:1 (Φ_ST = 0.018 for both markers). We conclude that between mallards and Mexican ducks, divergence at autosomal markers is largely neutral, whereas greater divergence on the Z chromosome (or some portions thereof) is likely the product of selection that has been important in speciation. Our results contribute to a growing body of literature indicating elevated divergence on the Z chromosome and its likely importance in avian speciation.     

Genetic structure of Quercus rubra L. and Quercus ellipsoidalis E. J. Hill populations at gene-based EST-SSR and nuclear SSR markers

Sympatric hybridizing oak species provide a model system for studying local adaptation. Disjunct populations of Quercus rubra L. and Quercus ellipsoidalis E. J. Hill at the northern edge of their distribution may harbor important reservoirs of adaptive genetic variation. Genic (expressed sequence tag- simple sequence repeat?=?EST-SSR) and non-genic nuclear microsatellite (nuclear SSR?=?nSSR) markers were used to estimate neutral and potentially adaptive genetic variation in these two supposedly interfertile oak species showing different adaptations to drought. Eleven populations of putative Q. rubra and Q. ellipsoidalis located in the Western Upper Peninsula of Michigan were characterized using seven EST-SSRs and eight nSSRs. Bayesian cluster analysis revealed two distinct groups corresponding to each species with evidence of low levels of potential introgression. A comparison of the genetic structure of adult trees and seedlings revealed no evidence for selection against hybrids. Overall, similar levels of genetic variation and differentiation between populations and species were found at both EST-SSRs and nSSRs indicating that most EST-SSRs chosen reflect neutral variation. Two loci, 3A05 (nSSR) and GOT021 (EST-SSR, putative histidine kinase 4-like), were identified as putative outlier loci between species showing largely reduced variation in Q. ellipsoidalis. Future analyses of an increased number of EST-SSRs located in functional genes will allow the identification of genes involved in the reproductive isolation between both species.     

Single nucleotide polymorphisms unravel hierarchical divergence and signatures of selection among Alaskan sockeye salmon (Oncorhynchus nerka) populations

Background  

Disentangling the roles of geography and ecology driving population divergence and distinguishing adaptive from neutral evolution at the molecular level have been common goals among evolutionary and conservation biologists. Using single nucleotide polymorphism (SNP) multilocus genotypes for 31 sockeye salmon (Oncorhynchus nerka) populations from the Kvichak River, Alaska, we assessed the relative roles of geography (discrete boundaries or continuous distance) and ecology (spawning habitat and timing) driving genetic divergence in this species at varying spatial scales within the drainage. We also evaluated two outlier detection methods to characterize candidate SNPs responding to environmental selection, emphasizing which mechanism(s) may maintain the genetic variation of outlier loci.     

Demographic history has shaped the strongly differentiated corkwing wrasse populations in Northern Europe

Understanding the biological processes involved in genetic differentiation and divergence between populations within species is a pivotal aim in evolutionary biology. One particular phenomenon that requires clarification is the maintenance of genetic barriers despite the high potential for gene flow in the marine environment. Such patterns have been attributed to limited dispersal or local adaptation, and to a lesser extent to the demographic history of the species. The corkwing wrasse (Symphodus melops) is an example of a marine fish species where regions of particular strong divergence are observed. One such genetic break occurred at a surprisingly small spatial scale (F_ST ~0.1), over a short coastline (<60 km) in the North Sea‐Skagerrak transition area in southwestern Norway. Here, we investigate the observed divergence and purported reproductive isolation using genome resequencing. Our results suggest that historical events during the post‐glacial recolonization route can explain the present population structure of the corkwing wrasse in the northeast Atlantic. While the divergence across the break is strong, we detected ongoing gene flow between populations over the break suggesting recent contact or negative selection against hybrids. Moreover, we found few outlier loci and no clear genomic regions potentially being under selection. We concluded that neutral processes and random genetic drift e.g., due to founder events during colonization have shaped the population structure in this species in Northern Europe. Our findings underline the need to take into account the demographic process in studies of divergence processes.     

The role of local ecology during hybridization at the initial stages of ecological speciation in a marine snail

Hybrid zones of ecologically divergent populations are ideal systems to study the interaction between natural selection and gene flow during the initial stages of speciation. Here, we perform an amplified fragment length polymorphism (AFLP) genome scan in parallel hybrid zones between divergent ecotypes of the marine snail Littorina saxatilis, which is considered a model case for the study of ecological speciation. Ridged‐Banded (RB) and Smooth‐Unbanded (SU) ecotypes are adapted to different shore levels and microhabitats, although they present a sympatric distribution at the mid‐shore where they meet and mate (partially assortatively). We used shell morphology, outlier and nonoutlier AFLP loci from RB, SU and hybrid specimens captured in sympatry to determine the level of phenotypic and genetic introgression. We found different levels of introgression at parallel hybrid zones and nonoutlier loci showed more gene flow with greater phenotypic introgression. These results were independent from the phylogeography of the studied populations, but not from the local ecological conditions. Genetic variation at outlier loci was highly correlated with phenotypic variation. In addition, we used the relationship between genetic and phenotypic variation to estimate the heritability of morphological traits and to identify potential Quantitative Trait Loci to be confirmed in future crosses. These results suggest that ecology (exogenous selection) plays an important role in this hybrid zone. Thus, ecologically based divergent natural selection is responsible, simultaneously, for both ecotype divergence and hybridization. On the other hand, genetic introgression occurs only at neutral loci (nonoutliers). In the future, genome‐wide studies and controlled crosses would give more valuable information about this process of speciation in the face of gene flow.     

Outlier SNP markers reveal fine‐scale genetic structuring across European hake populations (Merluccius merluccius)

Shallow population structure is generally reported for most marine fish and explained as a consequence of high dispersal, connectivity and large population size. Targeted gene analyses and more recently genome‐wide studies have challenged such view, suggesting that adaptive divergence might occur even when neutral markers provide genetic homogeneity across populations. Here, 381 SNPs located in transcribed regions were used to assess large‐ and fine‐scale population structure in the European hake (Merluccius merluccius), a widely distributed demersal species of high priority for the European fishery. Analysis of 850 individuals from 19 locations across the entire distribution range showed evidence for several outlier loci, with significantly higher resolving power. While 299 putatively neutral SNPs confirmed the genetic break between basins (F_CT = 0.016) and weak differentiation within basins, outlier loci revealed a dramatic divergence between Atlantic and Mediterranean populations (F_CT range 0.275–0.705) and fine‐scale significant population structure. Outlier loci separated North Sea and Northern Portugal populations from all other Atlantic samples and revealed a strong differentiation among Western, Central and Eastern Mediterranean geographical samples. Significant correlation of allele frequencies at outlier loci with seawater surface temperature and salinity supported the hypothesis that populations might be adapted to local conditions. Such evidence highlights the importance of integrating information from neutral and adaptive evolutionary patterns towards a better assessment of genetic diversity. Accordingly, the generated outlier SNP data could be used for tackling illegal practices in hake fishing and commercialization as well as to develop explicit spatial models for defining management units and stock boundaries.     

Interspecies introgressive hybridization in spiny frogs Quasipaa (Family Dicroglossidae) revealed by analyses on multiple mitochondrial and nuclear genes

Introgression may lead to discordant patterns of variation among loci and traits. For example, previous phylogeographic studies on the genus Quasipaa detected signs of genetic introgression from genetically and morphologically divergent Quasipaa shini or Quasipaa spinosa. In this study, we used mitochondrial and nuclear DNA sequence data to verify the widespread introgressive hybridization in the closely related species of the genus Quasipaa, evaluate the level of genetic diversity, and reveal the formation mechanism of introgressive hybridization. In Longsheng, Guangxi Province, signs of asymmetrical nuclear introgression were detected between Quasipaa boulengeri and Q. shini. Unidirectional mitochondrial introgression was revealed from Q. spinosa to Q. shini. By contrast, bidirectional mitochondrial gene introgression was detected between Q. spinosa and Q. shini in Lushan, Jiangxi Province. Our study also detected ancient hybridizations between a female Q. spinosa and a male Q. jiulongensis in Zhejiang Province. Analyses on mitochondrial and nuclear genes verified three candidate cryptic species in Q. spinosa, and a cryptic species may also exist in Q. boulengeri. However, no evidence of introgressive hybridization was found between Q. spinosa and Q. boulengeri. Quasipaa exilispinosa from all the sampling localities appeared to be deeply divergent from other communities. Our results suggest widespread introgressive hybridization in closely related species of Quasipaa and provide a fundamental basis for illumination of the forming mechanism of introgressive hybridization, classification of species, and biodiversity assessment in Quasipaa.