Patterns of genetic variation and local adaptation of a native herbivore to a lethal invasive plant

Understanding the evolutionary processes that influence fitness is critical to predicting species' responses to selection. Interactions among evolutionary processes including gene flow, drift and the strength of selection can lead to either local adaptation or maladaptation, especially in heterogenous landscapes. Populations experiencing novel environments or resources are ideal for understanding the mechanisms underlying adaptation or maladaptation, specifically in locally co-evolved interactions. We used the interaction between a native herbivore that oviposits on a patchily distributed introduced plant that in turn causes significant mortality to the larvae to test for signatures of local adaptation in areas where the two co-occurred. We used whole-genome sequencing to explore population structure, patterns of gene flow and signatures of local adaptation. We found signatures of local adaptation in response to the introduced plant in the absence of strong population structure with no genetic differentiation and low genetic variation. Additionally, we found localized allele frequency differences within a single population between habitats with and without the lethal plant, highlighting the effects of strong selection. Finally, we identified that selection was acting on larval ability to feed on the plant rather than on females' ability to avoid oviposition, thus uncovering the specific ontogenetic target of selection. Our work highlights the potential for adaptation to occur in a fine-grained landscape in the presence of gene flow and low genetic variation.

Understanding the joint evolution of dispersal and host specialisation using phytophagous arthropods as a model group

Theory generally predicts that host specialisation and dispersal should evolve jointly. Indeed, many models predict that specialists should be poor dispersers to avoid landing on unsuitable hosts while generalists will have high dispersal abilities. Phytophagous arthropods are an excellent group to test this prediction, given extensive variation in their host range and dispersal abilities. Here, we explore the degree to which the empirical literature on this group is in accordance with theoretical predictions. We first briefly outline the theoretical reasons to expect such a correlation. We then report empirical studies that measured both dispersal and the degree of specialisation in phytophagous arthropods. We find a correlation between dispersal and levels of specialisation in some studies, but with wide variation in this result. We then review theoretical attributes of species and environment that may blur this correlation, namely environmental grain, temporal heterogeneity, habitat selection, genetic architecture, and coevolution between plants and herbivores. We argue that theoretical models fail to account for important aspects, such as phenotypic plasticity and the impact of selective forces stemming from other biotic interactions, on both dispersal and specialisation. Next, we review empirical caveats in the study of this interplay. We find that studies use different measures of both dispersal and specialisation, hampering comparisons. Moreover, several studies do not provide independent measures of these two traits. Finally, variation in these traits may occur at scales that are not being considered. We conclude that this correlation is likely not to be expected from large-scale comparative analyses as it is highly context dependent and should not be considered in isolation from the factors that modulate it, such as environmental scale and heterogeneity, intrinsic traits or biotic interactions. A stronger crosstalk between theoretical and empirical studies is needed to understand better the prevalence and basis of the correlation between dispersal and specialisation.

Evolution of chemosensory and detoxification gene families across herbivorous Drosophilidae

Herbivorous insects are exceptionally diverse, accounting for a quarter of all known eukaryotic species, but the genomic basis of adaptations that enabled this dietary transition remains poorly understood. Many studies have suggested that expansions and contractions of chemosensory and detoxification gene families-genes directly mediating interactions with plant chemical defenses-underlie successful plant colonization. However, this hypothesis has been challenging to test because the origins of herbivory in many insect lineages are ancient (>150 million years ago (mya)), obscuring genomic evolutionary patterns. Here, we characterized chemosensory and detoxification gene family evolution across Scaptomyza, a genus nested within Drosophila that includes a recently derived (<15 mya) herbivore lineage of mustard (Brassicales) specialists and carnation (Caryophyllaceae) specialists, and several nonherbivorous species. Comparative genomic analyses revealed that herbivorous Scaptomyza has among the smallest chemosensory and detoxification gene repertoires across 12 drosophilid species surveyed. Rates of gene turnover averaged across the herbivore clade were significantly higher than background rates in over half of the surveyed gene families. However, gene turnover was more limited along the ancestral herbivore branch, with only gustatory receptors and odorant-binding proteins experiencing strong losses. The genes most significantly impacted by gene loss, duplication, or changes in selective constraint were those involved in detecting compounds associated with feeding on living plants (bitter or electrophilic phytotoxins) or their ancestral diet (fermenting plant volatiles). These results provide insight into the molecular and evolutionary mechanisms of plant-feeding adaptations and highlight gene candidates that have also been linked to other dietary transitions in Drosophila.

Evolution of chemosensory and detoxification gene families across herbivorous Drosophilidae

Herbivorous insects are exceptionally diverse, accounting for a quarter of all known eukaryotic species, but the genetic basis of adaptations that enabled this dietary transition remains poorly understood. Many studies have suggested that expansions and contractions of chemosensory and detoxification gene families - genes directly mediating interactions with plant chemical defenses - underlie successful plant colonization. However, this hypothesis has been challenging to test because the origins of herbivory in many lineages are ancient (>150 million years ago [mya]), obscuring genomic evolutionary patterns. Here, we characterized chemosensory and detoxification gene family evolution across Scaptomyza, a genus nested within Drosophila that includes a recently derived (<15 mya) herbivore lineage of mustard (Brassicales) specialists and carnation (Caryophyllaceae) specialists, and several non-herbivorous species. Comparative genomic analyses revealed that herbivorous Scaptomyza have among the smallest chemosensory and detoxification gene repertoires across 12 drosophilid species surveyed. Rates of gene turnover averaged across the herbivore clade were significantly higher than background rates in over half of the surveyed gene families. However, gene turnover was more limited along the ancestral herbivore branch, with only gustatory receptors and odorant binding proteins experiencing strong losses. The genes most significantly impacted by gene loss, duplication, or changes in selective constraint were those involved in detecting compounds associated with feeding on plants (bitter or electrophilic phytotoxins) or their ancestral diet (yeast and fruit volatiles). These results provide insight into the molecular and evolutionary mechanisms of plant-feeding adaptations and highlight strong gene candidates that have also been linked to other dietary transitions in Drosophila .

Niche Specificity, Polygeny, and Pleiotropy in Herbivorous Insects

AbstractWhat causes host use specificity in herbivorous insects? Population genetic models predict specialization when habitat preference can evolve and there is antagonistic pleiotropy at a performance-affecting locus. But empirically for herbivorous insects, host use performance is governed by many genetic loci, and antagonistic pleiotropy seems to be rare. Here, we use individual-based quantitative genetic simulation models to investigate the role of pleiotropy in the evolution of sympatric host use specialization when performance and preference are quantitative traits. We look first at pleiotropies affecting only host use performance. We find that when the host environment changes slowly, the evolution of host use specialization requires levels of antagonistic pleiotropy much higher than what has been observed in nature. On the other hand, with rapid environmental change or pronounced asymmetries in productivity across host species, the evolution of host use specialization readily occurs without pleiotropy. When pleiotropies affect preference as well as performance, even with slow environmental change and host species of equal productivity, we observe fluctuations in host use breadth, with mean specificity increasing with the pervasiveness of antagonistic pleiotropy. Thus, our simulations show that pleiotropy is not necessary for specialization, although it can be sufficient, provided it is extensive or multifarious.

Taste adaptations associated with host specialization in the specialist Drosophila sechellia

Chemosensory-driven host plant specialization is a major force mediating insect ecological adaptation and speciation. Drosophila sechellia, a species endemic to the Seychelles islands, feeds and oviposits on Morinda citrifolia almost exclusively. This fruit is harmless to D. sechellia but toxic to other Drosophilidae, including the closely related generalists D. simulans and D. melanogaster, because of its high content of fatty acids. While several olfactory adaptations mediating D. sechellia's preference for its host have been uncovered, the role of taste has been much less examined. We found that D. sechellia has reduced taste and feeding aversion to bitter compounds and host fatty acids that are aversive to D. melanogaster and D. simulans. The loss of aversion to canavanine, coumarin and fatty acids arose in the D. sechellia lineage, as its sister species D. simulans showed responses akin to those of D. melanogaster. Drosophila sechellia has increased taste and feeding responses towards M. citrifolia. These results are in line with D. sechellia's loss of genes that encode bitter gustatory receptors (GRs) in D. melanogaster. We found that two GR genes which are lost in D. sechellia, GR39a.a and GR28b.a, influence the reduction of aversive responses to some bitter compounds. Also, D. sechellia has increased appetite for a prominent host fatty acid compound that is toxic to its relatives. Our results support the hypothesis that changes in the taste system, specifically a reduction of sensitivity to bitter compounds that deter generalist ancestors, contribute to the specialization of D. sechellia for its host.

Host association, environment, and geography underlie genomic differentiation in a major forest pest

Diverse geographic, environmental, and ecological factors affect gene flow and adaptive genomic variation within species. With recent advances in landscape ecological modelling and high‐throughput DNA sequencing, it is now possible to effectively quantify and partition their relative contributions. Here, we use landscape genomics to identify determinants of genomic differentiation in the forest tent caterpillar, Malacosoma disstria, a widespread and irruptive pest of numerous deciduous tree species in North America. We collected larvae from multiple populations across Eastern Canada, where the species experiences a diversity of environmental gradients and feeds on a number of different host tree species, including trembling aspen (Populus tremuloides), sugar maple (Acer saccharum), red oak (Quercus rubra), and white birch (Betula papyrifera). Using a combination of reciprocal causal modelling (RCM) and distance‐based redundancy analyses (dbRDA), we show that differentiation of thousands of genome‐wide single nucleotide polymorphisms (SNPs) among individuals is best explained by a combination of isolation by distance, isolation by environment (spatial variation in summer temperatures and length of the growing season), and differences in host association. Configuration of suitable habitat inferred from ecological niche models was not significantly related to genomic differentiation, suggesting that M. disstria dispersal is agnostic with respect to habitat quality. Although population structure was not discretely related to host association, our modelling framework provides the first molecular evidence of host‐associated differentiation in M. disstria, congruent with previous documentation of reduced growth and survival of larvae moved between natal host species. We conclude that ecologically mediated selection is contributing to variation within M. disstria, and that divergent adaptation related to both environmental conditions and host association should be considered in ongoing research and management of this important forest pest.

Evolutionary Ecology of Plant-Arthropod Interactions in Light of the "Omics" Sciences: A Broad Guide

Aboveground plant-arthropod interactions are typically complex, involving herbivores, predators, pollinators, and various other guilds that can strongly affect plant fitness, directly or indirectly, and individually, synergistically, or antagonistically. However, little is known about how ongoing natural selection by these interacting guilds shapes the evolution of plants, i.e., how they affect the differential survival and reproduction of genotypes due to differences in phenotypes in an environment. Recent technological advances, including next-generation sequencing, metabolomics, and gene-editing technologies along with traditional experimental approaches (e.g., quantitative genetics experiments), have enabled far more comprehensive exploration of the genes and traits involved in complex ecological interactions. Connecting different levels of biological organization (genes to communities) will enhance the understanding of evolutionary interactions in complex communities, but this requires a multidisciplinary approach. Here, we review traditional and modern methods and concepts, then highlight future avenues for studying the evolution of plant-arthropod interactions (e.g., plant-herbivore-pollinator interactions). Besides promoting a fundamental understanding of plant-associated arthropod communities' genetic background and evolution, such knowledge can also help address many current global environmental challenges.

Support for the adaptive decoupling hypothesis from whole-transcriptome profiles of a hypermetamorphic and sexually dimorphic insect, Neodiprion lecontei

Though seemingly bizarre, the dramatic morphological and ecological transformation that occurs when immature life stages metamorphose into reproductive adults is one of the most successful developmental strategies on the planet. The adaptive decoupling hypothesis (ADH) proposes that metamorphosis is an adaptation for breaking developmental links between traits expressed in different life stages, thereby facilitating their independent evolution when exposed to opposing selection pressures. Here, we draw inspiration from the ADH to develop a conceptual framework for understanding changes in gene expression across ontogeny. We hypothesized that patterns of stage-biased and sex-biased gene expression are the product of both decoupling mechanisms and selection history. To test this hypothesis, we characterized transcriptome-wide patterns of gene-expression traits for three ecologically distinct larval stages (all male) and adult males and females of a hypermetamorphic insect (Neodiprion lecontei). We found that stage-biased gene expression was most pronounced between larval and adult males, which is consistent with the ADH. However, even in the absence of a metamorphic transition, considerable stage-biased expression was observed among morphologically and behaviourally distinct larval stages. Stage-biased expression was also observed across ecologically relevant Gene Ontology categories and genes, highlighting the role of ecology in shaping patterns of gene expression. We also found that the magnitude and prevalence of stage-biased expression far exceeded adult sex-biased expression. Overall, our results highlight how the ADH can shed light on transcriptome-wide patterns of gene expression in organisms with complex life cycles. For maximal insight, detailed knowledge of organismal ecology is also essential.

Why are there not more herbivorous insect species?

Insect species richness is estimated to exceed three million species, of which roughly half is herbivorous. Despite the vast number of species and varied life histories, the proportion of herbivorous species among plant-consuming organisms is lower than it could be due to constraints that impose limits to their diversification. These include ecological factors, such as vague interspecific competition; anatomical and physiological limits, such as neural limits and inability of handling a wide range of plant allelochemicals; phylogenetic constraints, like niche conservatism; and most importantly, a low level of concerted genetic variation necessary to a phyletic conversion. It is suggested that diversification ultimately depends on what we call the intrinsic trend of diversification of the insect genome. In support of the above, we survey the major types of host-specificity, the mechanisms and constraints of host specialization, possible pathways of speciation, and hypotheses concerning insect diversification.

Soft Selective Sweep on Chemosensory Genes Correlates with Ancestral Preference for Toxic Noni in a Specialist Drosophila Population

Understanding how organisms adapt to environmental changes is a major question in evolution and ecology. In particular, the role of ancestral variation in rapid adaptation remains unclear because its trace on genetic variation, known as soft selective sweep, is often hardly recognizable from genome-wide selection scans. Here, we investigate the evolution of chemosensory genes in Drosophila yakuba mayottensis, a specialist subspecies on toxic noni (Morinda citrifolia) fruits on the island of Mayotte. We combine population genomics analyses and behavioral assays to evaluate the level of divergence in chemosensory genes and perception of noni chemicals between specialist and generalist subspecies of D. yakuba. We identify a signal of soft selective sweep on a handful of genes, with the most diverging ones involving a cluster of gustatory receptors expressed in bitter-sensing neurons. Our results highlight the potential role of ancestral genetic variation in promoting host plant specialization in herbivorous insects and identify a number of candidate genes underlying behavioral adaptation.

Evaluating insect-host interactions as a driver of species divergence in palm flower weevils

Plants and their specialized flower visitors provide valuable insights into the evolutionary consequences of species interactions. In particular, antagonistic interactions between insects and plants have often been invoked as a major driver of diversification. Here we use a tropical community of palms and their specialized insect flower visitors to test whether antagonisms lead to higher population divergence. Interactions between palms and the insects visiting their flowers range from brood pollination to florivory and commensalism, with the latter being species that feed on decaying-and presumably undefended-plant tissues. We test the role of insect-host interactions in the early stages of diversification of nine species of beetles sharing host plants and geographical ranges by first delimiting cryptic species and then using models of genetic isolation by environment. The degree to which insect populations are structured by the genetic divergence of plant populations varies. A hierarchical model reveals that this variation is largely uncorrelated with the kind of interaction, showing that antagonistic interactions are not associated with higher genetic differentiation. Other aspects of host use that affect plant-associated insects regardless of the outcomes of their interactions, such as sensory biases, are likely more general drivers of insect population divergence.

What We Don't Know About Diet-Breadth Evolution in Herbivorous Insects

Half a million species of herbivorous insects have been described. Most of them are diet specialists, using only a few plant species as hosts. Biologists suspect that their specificity is key to their diversity. But why do herbivorous insects tend to be diet specialists? In this review, we catalog a broad range of explanations. We review the evidence for each and suggest lines of research to obtain the evidence we lack. We then draw attention to a second major question, namely how changes in diet breadth affect the rest of a species’ biology. In particular, we know little about how changes in diet breadth feed back on genetic architecture, the population genetic environment, and other aspects of a species’ ecology. Knowing more about how generalists and specialists differ should go a long way toward sorting out potential explanations of specificity, and yield a deeper understanding of herbivorous insect diversity.

The genetic architecture of a host shift: An adaptive walk protected an aphid and its endosymbiont from plant chemical defenses

Host shifts can lead to ecological speciation and the emergence of new pests and pathogens. However, the mutational events that facilitate the exploitation of novel hosts are poorly understood. Here, we characterize an adaptive walk underpinning the host shift of the aphid Myzus persicae to tobacco, including evolution of mechanisms that overcame tobacco chemical defenses. A series of mutational events added as many as 1.5 million nucleotides to the genome of the tobacco-adapted subspecies, M. p. nicotianae, and yielded profound increases in expression of an enzyme that efficiently detoxifies nicotine, both in aphid gut tissue and in the bacteriocytes housing the obligate aphid symbiont Buchnera aphidicola. This dual evolutionary solution overcame the challenge of preserving fitness of a mutualistic symbiosis during adaptation to a toxic novel host. Our results reveal the intricate processes by which genetic novelty can arise and drive the evolution of key innovations required for ecological adaptation.

Avoidance and aggregation create consistent egg distribution patterns of congeneric caddisflies across spatially variable oviposition landscapes

Amongst oviparous animals, the spatial distribution of individuals is often set initially by where females lay eggs, with potential implications for populations and species coexistence. Do the spatial arrangements of oviposition sites or female behaviours determine spatial patterns of eggs? The consequences of spatial patterns may be context independent if strong behaviours drive patterns; context dependent if the local environment dominates. We tested these ideas using a guild of stream-dwelling caddisflies that oviposit on emergent rocks, focussing on genera with contrasting behaviours. In naturally occurring oviposition landscapes (riffles with emergent rocks), we surveyed the spatial arrangement and environmental characteristics of all emergent rocks, identified and enumerated egg masses on each. Multiple riffles were surveyed to test for spatially invariant patterns and behaviours. In landscapes, we tested for spatial clumping of oviposition sites exploited by each species and for segregation of congeneric species. At oviposition sites, we characterised the frequency distributions of egg masses and tested for species associations. Genus-specific behaviours produced different spatial patterns of egg masses in the same landscapes. Congregative behaviour of Ulmerochorema spp. at landscape scales and an aggregative response at preferred oviposition sites led to clumped patterns, local aggregation and species overlap. In contrast, avoidance behaviours by congeners of Apsilochorema resulted in no or weak clumping, and species segregation in some landscapes. Spatial patterns were consistent across riffles that varied in area and oviposition site density. These results suggest that quite different oviposition behaviours may be context independent, and the consequences of spatial patterns may be spatially invariant also.

Evolutionary genomics of host plant adaptation: insights from Drosophila

Variation in gene expression in response to the use of alternate host plants can reveal genetic and physiological mechanisms explaining why insect-host relationships vary from host specialism to generalism. Interpreting transcriptome variation relies on well-annotated genomes, making drosophilids valuable model systems, particularly those species with tractable ecological associations. Patterns of whole genome expression and alternate gene splicing in response to growth on different hosts have revealed expression of gene networks of known detoxification genes as well as novel functionally enriched genes of diverse metabolic and structural functions. Integrating trancriptomic responses with fitness differences and levels of phenotypic plasticity in response to alternate hosts will help to reveal the general nature of genotype-phenotype relationships.

Gene Expression and Diet Breadth in Plant-Feeding Insects: Summarizing Trends

Transcriptomic studies lend insights into the role of transcriptional plasticity in adaptation and specialization. Recently, there has been growing interest in understanding the relationship between variation in herbivorous insect gene expression and the evolution of diet breadth. We review the studies that have emerged on insect gene expression and host plant use, and outline the questions and approaches in the field. Many candidate genes underlying herbivory and specialization have been identified, and a few key studies demonstrate increased transcriptional plasticity associated with generalist compared with specialist species. Addressing the roles that transcriptional variation plays in insect diet breadth will have important implications for our understanding of the evolution of specialization and the genetic and environmental factors that govern insect-plant interactions.

Host plant-related genomic differentiation in the European cherry fruit fly, Rhagoletis cerasi

Elucidating the mechanisms and conditions facilitating the formation of biodiversity are central topics in evolutionary biology. A growing number of studies imply that divergent ecological selection may often play a critical role in speciation by counteracting the homogenising effects of gene flow. Several examples involve phytophagous insects, where divergent selection pressures associated with host plant shifts may generate reproductive isolation, promoting speciation. Here, we use ddRADseq to assess the population structure and to test for host‐related genomic differentiation in the European cherry fruit fly, Rhagoletis cerasi (L., 1758) (Diptera: Tephritidae). This tephritid is distributed throughout Europe and western Asia, and has adapted to two different genera of host plants, Prunus spp. (cherries) and Lonicera spp. (honeysuckle). Our data imply that geographic distance and geomorphic barriers serve as the primary factors shaping genetic population structure across the species range. Locally, however, flies genetically cluster according to host plant, with consistent allele frequency differences displayed by a subset of loci between Prunus and Lonicera flies across four sites surveyed in Germany and Norway. These 17 loci display significantly higher F_ST values between host plants than others. They also showed high levels of linkage disequilibrium within and between Prunus and Lonicera flies, supporting host‐related selection and reduced gene flow. Our findings support the existence of sympatric host races in R. cerasi embedded within broader patterns of geographic variation in the fly, similar to the related apple maggot, Rhagoletis pomonella, in North America.

Predicting evolutionary predictability

The observation that phenotypic convergence and genetic convergence are widespread in nature implies that evolution is at least somewhat predictable. But to what extent and under what circumstances? In other words, how predictable is evolutionary predictability? Answering this question requires going beyond documenting examples of repeated evolution to actually quantifying predictability at different hierarchical levels. At present, few such studies exist. In this issue of Molecular Ecology, Chaturvedi et al. () quantify the predictability of genomewide changes that accompany shifts to an introduced host plant (alfalfa) in populations of the Melissa blue butterfly (Lycaeides melissa). They evaluate predictability in two contexts: (i) overlap in host-associated loci among populations that have independently colonized alfalfa, and (ii) overlap between host-associated loci in nature and loci associated with host performance in laboratory experiments. Overall, they find that the genomic changes that accompany host shifts in this system are indeed somewhat predictable. However, the degree of predictability depends on the type of comparison (among natural populations vs. between natural and experimental populations), type of convergence (specific genomic locations vs. direction of allele frequency change), geographic scale (rangewide vs. specific population pairs) and location in the genome (autosomes vs. sex chromosomes). Together with a handful of comparable data sets, Chaturvedi et al.'s () work suggests that the relative contribution of stochastic and deterministic processes to genomewide responses to novel selection pressures may be highly variable, but possibly predictably so.

Genomic transitions during host race and species formation

Darwin recognized species as discontinuous, yet considered them to be formed by an incremental process of natural selection. Recent theoretical work on 'genome-wide congealing' is bridging this gap between the gradualism of divergent selection and rapid genome-wide divergence, particularly during ecological speciation-with-gene-flow. Host races and species of phytophagous insects, displaying a spectrum of divergence and gene flow among member taxa, provide model systems for testing predicted non-linear transitions from 'genic' divergence at a few uncoupled loci to 'genomic' divergence with genome-wide coupling of selected loci and strong reproductive isolation. Integrating across natural history, genomics, and evolutionary theory, emerging research suggests a tipping point from 'genic' to 'genomic' divergence between host races and species, during both sympatric speciation and secondary contact.

Positive selection at sites of chemosensory genes is associated with the recent divergence and local ecological adaptation in cactophilic Drosophila

Background

Adaptation to new hosts in phytophagous insects often involves mechanisms of host recognition by genes of sensory pathways. Most often the molecular evolution of sensory genes has been explained in the context of the birth-and-death model. The role of positive selection is less understood, especially associated with host adaptation and specialization. Here we aim to contribute evidence for this latter hypothesis by considering the case of Drosophila mojavensis, a species with an evolutionary history shaped by multiple host shifts in a relatively short time scale, and its generalist sister species, D. arizonae.

Results

We used a phylogenetic and population genetic analysis framework to test for positive selection in a subset of four chemoreceptor genes, one gustatory receptor (Gr) and three odorant receptors (Or), for which their expression has been previously associated with host shifts. We found strong evidence of positive selection at several amino acid sites in all genes investigated, most of which exhibited changes predicted to cause functional effects in these transmembrane proteins. A significant portion of the sites identified as evolving positively were largely found in the cytoplasmic region, although a few were also present in the extracellular domains.

Conclusions

The pattern of substitution observed suggests that some of these changes likely had an effect on signal transduction as well as odorant recognition and protein-protein interactions. These findings support the role of positive selection in shaping the pattern of variation at chemosensory receptors, both during the specialization onto one or a few related hosts, but as well as during the evolution and adaptation of generalist species into utilizing several hosts.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1250-x) contains supplementary material, which is available to authorized users.