According to Chinese tradition Fu Xi (Fu Hsi) was China's first historical ruler (c. 2952-2836 BCE). A man of great diversity, he is often depicted as a serpent with a human's head. Fu Xi was born in the Tianshui region of the Weishui River Valley in Gansu province.
Ancient Chinese legend holds Fu Xi as the first person to discover farming, hunting, fishing with nets, cooking meals in pots over fire, the domestication of animals, the breeding of silk worms, the regulation of marriage, use of the mortar and pestle, and teaching of the eight trigrams (bagua) via divination with yarrow stalks. It is in this last area where the "physical" system of yin/yang dialectics originated in the form of stacked, solid and broken line symbols that were arranged in equal but opposite positions in a circular format.
arrows with bone tips. The Yang-Shao produced red-coiled pottery that were painted with black geometric designs. They also domesticated pigs, practiced silkworm cultivation, made basketry and fiber cloth from hemp and silk.
A typical example of a Yang-Shao settlement was Banpo, which stood on seven acres of land near a small tributary of the Weishui River. Excavated more than 40 years ago, it has been estimated that nearly 400 inhabitants lived here, in round houses constructed of mud and grass, supported by wooden pillars, and in rectangular houses with pitched roofs made of thatch. Meeting houses and cemeteries were also central in these villages as well. Found in many of the grave sites were simple ornaments but there are no traces of images or talismans to accommodate the dead, although some remains indicate shamanistic beliefs.
Additional archaeological evidence taken at a grave site by the Gansu Archaeology Research Foundation* from the Maiji Mountain area of Tianshui also reveals paintings on pottery symbolizing the eight trigrams, fishing nets, crafts and stoneware, a rectangular house foundation, and human remains. Because the contents from this one ancient grave site differed dramatically from other grave sites excavated in the area, it appears that the particular person buried here was unique from the rest of community. Some scholars believe it is the grave of Fu Xi.
*(See Li Shi Min's "Discussion of Fu Hsi Culture and the Maiji Mountains" in The Culture of Fu Xi, Chen Ju-Yuan, ed., 1995, Beijing) This thinking stands to reason, since no other artifacts in this nor any other archeological site in China, have produced trigram paintings. Fishing nets, too, although rare, have been found in other areas. This fact, however does not mean that Fu Hsi didn't invent the fishing net, but rather suggests that he could have and other members of his culture used them as well. Fishing nets do not appear prior to the culture of the Yang-Shao.
A later prehistoric Chinese culture called the Lung-Shan began around the time that the Yang-Shao culture ended and was situated more toward the eastern coast of China. Although the Lung-Shan appear to have been more highly evolved than the Yang-Shao, there are many shared similarities between the two cultures. Many scholars believe that the Lung-Shan may be derived from the Yang-Shao or a synthesis of Yang-Shao and Dawenkou (5000 – 3000 BCE) cultures. During the Dawenkou period, however, we find the first indications of Chinese writing in connection with the use of heated animal bones for fortune-telling. This practice was later elevated to a more ritualized procedure during the Shang dynasty (1850 – 1100 BCE) that followed later. Nowhere in China's prehistoric record, except during Fu Xi's legendary reign in the Tianshui region, have divination practices that utilize trigrams (or yarrow stalks) been found to exist.
"In ancient times, Master Fu Xi ruled over all. He raised his head and observed the heavens, lowered his head and observed the earth, observed the natural patterning of birds and animals, compared them to the disposition of the earth.
** translation by Stephen Selby.
1. Yin and yang energy regulate all physical phenomena.
3. At the extremes, yin changes to yang and yang changes to yin.
Note that there are eight phases to each hexagram's existence before the bottom trigram changes into a different trigram. Each fully changed trigram correlates with one of the "children" originally formed in the first sequence of trigrams. Each of the primary trigrams (or "children") also have physiological attributes based on their level of maturity. Together they represent the full spectrum of the growth and decay cycles of life.
To find what particular level of thinking an individual is at, and hence to seek guidance for particular problems associated with that level of thought, Fu Xi devised a method that could trip time by using a divination method that employed 50 yarrow stalks. This, in a nutshell, is what the original Book of Change (I Ching or Yijing ) is all about. How far these ideas extended into the various regions of China during this time is not known, nor is it known by whom this practice was relegated to. We do know, however, that by the time of the Shang Dynasty we begin to see references to yin and yang in the form of pictographs. These primitive characters have been found engraved into ox bones, tortoise shells, and later, sea shells. These substrates all were used for divination rituals which allowed an official soothsayer to read cracks cause by heating the engraved materials over fire. Divination, during this time was reserved for kings, or the acquisition material control and self protection.
It is also during the time of the Shang dynasty when we begin to see a change from the earth worshiping Neolithic cultures to the worship of a heavenly deity, named Shang Ti ("upper ruler". ) The paternal figure did not have a clearly defined character but his influence over human affairs was so great that no important matter could be attended to unless the advice of Shang Ti was first sought. The psychological advantage of having the advice of Shang Ti communicating only to Shang kings via divination was also a decisive political move. Whatever advice Shang Ti gave the rulers, that advice could not be challenged. As Shang Ti grew in power, the spirit of the individual was reduced — oppression and human sacrifices soon became widespread.
It is important to also note that this change of orientation from earth-centered thinking to sky-centered thinking was a pivotal point in the evolution of Chinese philosophy and religion.
Notice in the King Wen arrangement that the trigrams which were greatest in opposition, the Yielding and the Firm, have been replaced by the Abysmal and the Clear — each containing both solid and broken lines. The bottom, foundation lines of each trigram in this new order also show no conformity in regards to yin/yang placement: there are both "male" and "female" trigrams on both left and right sides of the newly formed central axis — their progression does not change from yin to yang in gradual succession, nor is polar reversal indicated. Hence, the idea of balance within a cycle of growth and decay is completely missing.
In the new hexagram arrangement, we continue to find additional inconsistencies. Remember in King Wen's trigram arrangement the Abysmal was placed in the extreme north. In King Wen's hexagram arrangement the most extreme, "male" hexagram (the Firm) is at the north (yin) instead. This time around, however, the Firm (six solid lines) has its name changed to Heaven. To the right (not opposite) of Heaven is the Yielding (six broken lines,) only this time around her name has been changed to Earth. It is at this point in China's history when we begin to see developing the idea of "heaven" having full male attributes and the "earth" with female attributes, in a position that is secondary to the male.
To differentiate between the Fu Xi and King Wen concepts of yin and yang, George Ohsawa, decided to call the Fu Xi concept physical and King Wen's concept metaphysical. For example; in the macrobiotic description of the physical body's internal organs, all solid organs are listed as yang (heart, lungs, kidneys, liver, spleen/pancreas) and all hollow organs as yin (small intestine, stomach, large intestine, bladder, gall bladder.) In the paradigm of Traditional Chinese Medicine, the yin/yang description of these same organs are in reverse — solid organs are designated as yin and hollow organs as yang.
Along with the new graphs, King Wen attempted to write a text that would explain their main attributes — later to be known as the "judgments." The language of these so-called "judgments" is sparsely phrased — more in line with the cryptic messages used in earlier divination queries. Yet, it is this confusing body of text, called Guazi, which forms the central framework of the I Ching.
After King Wen's death, his son Duke Zhou, appended the Guazi by attempting to explain the meaning of the hexagrams from a line-by-line breakdown using local, historical references. Additional commentaries by other authors — mainly Confucian scholars, were added to the I Ching over the following centuries.
Despite its obvious transformation and break from the Fu Xi order, the revised face of the I Ching managed to maintain its revered position as the official court oracle, using the hexagram system instead of oracle bones, which eventually became obsolete.
Toward the end of the Zhou era political disruption and social instability had become paramount in many regions throughout China. The class barrier system began to disintegrate and the yin/yang paradigm began to spread among the peasantry. Eventually the practice of divination fell from its coveted position in the royalty directly into the hands of commoners. These "commoners", however, were far from common when it came resurrecting ancient hilosophical principles. The revised face of the trigram and hexagrams made little sense from the yin/yang standpoint, yet the idea of the ancient ridge pole that divided equally all phenomena that directly opposed and complimented one another began to grow in popularity.
And the Three gave birth to the ten thousand things.
and wrap their arms around Yang.
Lao Tzu stayed clear of the I Ching, believing that it was a degenerate practice which only confused the underlying message that it originally conveyed.
with the Fu Xi/Lao Tzu physical concept of yin and yang or with the King Wen/Confucian metaphysical concept. The list is long but is clearly highlighted along the long span of time that ends with today.
I'm on day 14 of my juice cleanse and my body has been feeling really, really good for the last few days!
How Much is Enough Plant Based Protein?


Lovely Banks is a future residential township development for Riverlee, 55km Southwest of Melbourne and 10km North of Geelong.
Lovely Banks will create a new township connected to both Geelong and Greater Melbourne, with 10,000 new dwellings, 30,000 residents and 10,000 jobs.
With an additional 37,936 people forecasted to reside in the North Geelong region by 2026 (MACROPLAN DIMASI 2012), the site is an integral part of the region that will bring together sustainable environment, transport, education, and employment practices so that the community and its residents thrive long into the future.


On the afternoon of Friday 22nd September 2017, children, parents and carers gathered in the Junior hall for the "Mystery Bag Exhibition". For the last few weeks children from every class (including Daycare) received a series of eight bags containing recycled material. The children then worked in teams to try to create a work of art using only the contents of each bag.
The children have been working so very hard as you can see in the pictures below. May we take this opportunity to thank all parents for coming to our exhibition. We do hope that you enjoyed the amazing works of art created by our gifted and talented children.


"Dream Journey with Miro" is a digital picture book in which "seeing" is the keyword. Let's peek into a mysterious dream world with Miro. When Miro fell asleep holding her precious book, she awoke to find herself in an unfamiliar place. It was a world inside her dream and she met a butterfly. The story unfolds as she touches the mysterious butterfly. She finds herself in a dessert, at sea, and on the palm of a giant. What is the conclusion to the dream after her mysterious encounters?
Unknown worlds, and mysterious worlds of color. An app for finding worlds that you never knew existed.
This children's app were produced in a class titled "Project & Collaboration Seminar", held for students of the Department of Cross-Disciplinary Art and Design (Field of Media Art and Design) at Joshibi University of Art and Design in Tokyo, Japan. The students were entirely responsible for the planning, storywriting, illustration, animation and authoring of their apps over a six-week period. Children were invited to try out the finished programs, after which successive modifications were made.


text:
→ NAiM Playboy Arch.
Stedelijk Museum Open <
> Interview / Graphic
NAiM Playboy Arch.
Playboy Architecture 1953–1979
Since 2007, we've been responsible for the graphic identity of a Maastricht-based architecture institute which originally went under the name NAi Maastricht, but in 2009 changed its moniker to NAiM / Bureau Europa. (In 2013, the institute even dropped the NAiM-part, and is currently simply known as Bureau Europa).
'Playboy Architecture 1953–1979' took place between September 29, 2012, and February 10, 2013, and was the last NAiM exhibition that took place under the directorship of Guus Beumer. The exhibition was curated by the M+M (Media and Modernity) Program at the Princeton University School of Architecture, a post-graduate research project headed by Beatriz Colomina. The aim of the project was to explore the radical modernist dimensions of Playboy Magazine, and to investigate the magazine's influence on architecture and other forms of design (furniture, interior, etc.).
Exhibition view (click for more images)
As with most of the exhibitions taking place at NAiM / Bureau Europa, the spatial constructions were designed by EventArchitectuur. Our part is then to come up with a (typo-)graphic intervention, specifically responding to these constructions. (Added to that, we are also responsible for all the other printed matter: posters, invitations, window typography, adverts, etc.).
'Playboy Architecture 1953–1979' was a huge exhibition, taking up both halls of the NAiM building. EventArchitectuur constructed a group of 12 separate installations, each space referring to a specific theme within the subject. In the meantime, we came up with a way to include, in these 12 constructions, all the information that was provided to us by the researchers of the M+M Program.
While working on the assignment, our studio was actually filled with a large amount of cardboard boxes, containing all issues of Playboy published between 1953 and 1979. Based on all the individual research projects, our task was to select and scan all the images needed to illustrate the given themes.
Our initial goal was to introduce a 'sub-theme': the theme of the 'blow-up'. We hoped that, by gradually 'blowing up' the images, we would reveal the material base of printed material, in an attempt to undermine the images' illusionary dimension (or at least, to strip away some of the material's overpowering 'sexiness'). Obviously, the 'blow-up' is a theme we worked with before, for example in this print we designed in 2007.
However, since we had to work within a limited amount of space, and had to use a certain number of specific images, the whole 'blow-up' theme couldn't be executed as well as we originally intended. But traces of it can still be found in the typographic layer of the show.
The way we placed the images and text in the wooden constructions refers to the material reality of printed matter as well. By printing all the graphic material on long paper strips of paper, we hoped to invoke the tactility of magazines. It was also a nod to the notion of the printing press, and the way in which, on these presses, long rolls of paper are weaved between the rolls of the machine.
Finally, the long strips of paper also refer to a particular sequence in the movie 'Blow-Up' (Antonioni, 1966): the 'photo-shoot' scene, in which the main character photographs a group of models situated between long, vertical panels:
We just pinned these long strips of paper to the wooden constructions using casual thumbtacks (our intention has always been to come up with 'non-monumental' ways of displaying material, so that's why we often make use of pushpins, thumbtacks, clips, staples or pieces of tape). The long strips of paper in the exhibition were actually printed by QS Graphics in Maastricht.
A central part of the exhibition was a set of huge wooden shelves, displaying all the issues of Playboy published between 1953 and 1979. These issues were shown in chronological order – but we also wanted to introduce a way to explore these issues thematically. That's why we designed a system of colour-coded, bunny-ear-shaped bookmarks, each colour referring to a specific theme within the subject:
As part of the exhibition, one of the researchers (Federica Vannucchi) asked us to make a publication containing material from her research project. Her idea was to make a 'Playboy's Buying Guide', featuring a selection of 'shopping pages' reproduced from issues of Playboy. We designed this publication as a low-key booklet – basically a stack of reproduced pages, printed in black ink on pink paper, overprinted with blue ink for the dates, folded in two, and then stapled. A5-sized, this booklet was printed by Dijkman Offset:
Playboy Buying guide (click for more images)
Another piece of printed matter we designed for the exhibition was the usual set of posters. In our design for these posters, we tried to deconstruct the iconic bunny logo in such a way that the result resembled an actual piece of space-age architecture, the eyes of the bunny resembling windows. Thus, the composition becomes a 'multi-eyed' entity, referring to some of the sub-themes within the exhibition (voyeurism, exhibitionism, panopticonism).
The colour pink refers to the notion of flesh, but is also an attempt to somewhat undermine the 'machismo' of the subject of the exhibition (after all, pink is often seen as an effeminate colour).
The A0-sized poster was screenprinted by Wyber Zeefdruk:
The A2-sized variation was offset-printed by Lenoirschuring:
In the same manner, we designed an A3-sized invitation (to be folded to an A5-sized mailer). This item was again offset-printed by Lenoirschuring:
The window typography (here seen from the back) consisted again of the reconstructed 'bunny-pattern', the eyes of the bunnies forming little round windows, for passers-by to look through:
All the above photographs were either made by ourselves, or by Johannes Schwartz. More photos can be found on the websites of Bureau Europa and EventArchitectuur. The exhibition also included a soundtrack compiled by Femke Dekker (of My Little Underground).
next: Interview / Graphic
prev: Stedelijk Museum Open

redpajama_set_name: RedPajamaCommonCrawl


text:
BMC Ecology and Evolution
Growth genes are implicated in the evolutionary divergence of sympatric piscivorous and insectivorous rainbow trout (Oncorhynchus mykiss)
Jared A. Grummer1,
Michael C. Whitlock1,
Patricia M. Schulte1 &
Eric B. Taylor1,2
BMC Ecology and Evolution volume 21, Article number: 63 (2021) Cite this article
Identifying ecologically significant phenotypic traits and the genomic mechanisms that underly them are crucial steps in understanding traits associated with population divergence. We used genome-wide data to identify genomic regions associated with key traits that distinguish two ecomorphs of rainbow trout (Oncorhynchus mykiss)—insectivores and piscivores—that coexist for the non-breeding portion of the year in Kootenay Lake, southeastern British Columbia. "Gerrards" are large-bodied, rapidly growing piscivores with high metabolic rates that spawn north of Kootenay Lake in the Lardeau River, in contrast to the insectivorous populations that are on average smaller in body size, with lower growth and metabolic rates, mainly forage on aquatic insects, and spawn in tributaries immediately surrounding Kootenay Lake. We used pool-seq data representing ~ 60% of the genome and 80 fish per population to assess the level of genomic divergence between ecomorphs and to identify and interrogate loci that may play functional or selective roles in their divergence.
Genomic divergence was high between sympatric insectivores and piscivores (\(F_{\text{ST}}\) = 0.188), and in fact higher than between insectivorous populations from Kootenay Lake and the Blackwater River (\(F_{\text{ST}}\) = 0.159) that are > 500 km apart. A window-based \(F_{\text{ST}}\) analysis did not reveal "islands" of genomic differentiation; however, the window with highest \(F_{\text{ST}}\) estimate did include a gene associated with insulin secretion. Although we explored the use of the "Local score" approach to identify genomic outlier regions, this method was ultimately not used because simulations revealed a high false discovery rate (~ 20%). Gene ontology (GO) analysis identified several growth processes as enriched in genes occurring in the ~ 200 most divergent genomic windows, indicating many loci of small effect involved in growth and growth-related metabolic processes are associated with the divergence of these ecomorphs.
Our results reveal a high degree of genomic differentiation between piscivorous and insectivorous populations and indicate that the large body piscivorous phenotype is likely not due to one or a few loci of large effect. Rather, the piscivore phenotype may be controlled by several loci of small effect, thus highlighting the power of whole-genome resequencing in identifying genomic regions underlying population-level phenotypic divergences.
Understanding the divergence of populations and the mechanisms responsible for it are important for the fields of speciation, macroevolution, ecology, and conservation [1]. In many cases, population-level divergences are associated with character displacement and phenotypic divergence of sympatric forms, such as in threespine sticklebacks (e.g., [2]) or Nicaraguan crater lake cichlids (e.g., [3]). Recent advancements in sequencing technologies have enabled investigations into the genomics underlying traits associated with phenotypic divergences in non-model systems. In some cases, few highly differentiated loci of large effect can be responsible for population-level phenotypic or behavioural divergence in spite of the majority of genomic differentiation being very low between populations as the result of either recent divergence and/or ongoing gene flow (e.g., Heliconius butterflies; [4, 5]).
Determining the genetic architecture of a phenotypic trait—whether the trait is controlled by few loci of large effect, many loci of small effect, or some combination thereof—is important for understanding how traits respond to evolutionary processes including selection and migration [6]. In fishes, there are a few examples of a small number of loci having an effect on phenotypic divergence between sympatric forms. For instance, Jacobs et al. [7] identified ~ 10 genomic outlier regions in brown trout (Salmo trutta) containing several genes involved in development, growth, and gene expression and regulation that differ between piscivorous (fish-eating) and insectivorous populations. In Atlantic cod (Gadus morhua), a ~ 17 MB (megabase) region containing two chromosomal inversions is associated with population divergence based on migration phenotype [8]. Additionally, in Pacific salmon (genus Oncorhynchus), one small genomic region has been shown to control a vital phenotype for the natural history of two species—adult migration timing [9].
Theory predicts that the distribution of allelic effects on quantitative traits should be nearly exponential, where few loci have large effects and explain most of the variance in traits, but many loci with smaller effects explain the remaining variation in the trait [10, 11]. This pattern was seen for "age at maturity" in Atlantic salmon, where two genes of large effect and > 100 genes of small effect underly this trait [12]. Across taxa, several phenotypes are explained by many loci of small effect. For instance, hundreds of genetic variants explain only ~ 10% of the phenotypic variation in human height [13]. And a genome-wide analysis of 50,000 SNPs (single nucleotide polymorphism) in the collared flycatcher (Ficedula albicollis) revealed that no SNPs were significantly associated with forehead patch size, a sexually selected trait [14]. Furthermore, Kardos et al. [14] report that thousands of individuals and near-complete genome sequencing are necessary to reliably detect large-effect loci.
Identifying the genomic regions underlying divergent ecological forms can provide insights into the genomic architecture of population divergence. Here, we address this question in two ecomorphs of rainbow trout (Oncorhynchus mykiss) in Kootenay Lake of southeastern British Columbia. Rainbow trout in this lake provide an example of divergence into co-existing insectivorous and piscivorous forms, which has repeatedly been observed in trout and char. For example, Arctic charr (Salvelinus alpinus) have repeatedly evolved sympatric forms of zooplanktivores and piscivores on multiple continents [15, 16], including Iceland's Lake Thingvallavatn (e.g., [17]). Additionally, a large and piscivorous form of brown trout ("Ferox" trout; Salmo trutta) occurs throughout Ireland and Scotland [18, 19]. A variety of divergent life-history patterns can be found in populations across the native range of O. mykiss, including migration patterns (resident vs. anadromous), habitat type (lake vs. stream), and diet (piscivory vs. insectivory) [20]. Whereas genomic regions have been found that are associated with migration and habitat type [9, 21], little is known about the evolutionary mechanisms promoting divergence of piscivorous and insectivorous O. mykiss populations.
The rainbow trout ecomorphs in Kootenay Lake differ in a wide variety of ecological and life-history traits. Individuals of the large piscivorous form (often > 60 cm and 5 kg, but up to 12 kg) are locally known as "Gerrards" and have an adult dietary specialization on Kokanee salmon—small, land-locked sockeye salmon (O. nerka)—and spawn near the abandoned town site Gerrard (thus the name "Gerrards") on the Lardeau River that drains into the north arm of Kootenay Lake [22, 23]. In contrast, the smaller form (typically \(\le\) 50 cm and \(\le\) 2 kg) is primarily insectivorous and made up of several populations that spawn in smaller tributaries that drain into Kootenay Lake [23] (Fig. 1). The piscivorous and insectivorous forms are sympatric in the lake for a majority of the year (e.g., non-breeding season) after spending their first 1–2 years separated as juveniles in their respective streams. Piscivorous individuals have significantly larger head and upper jaw sizes, likely due to their specialization on larger prey (Neufeld, BC Ministry of Forests, Lands, Natural Resource Operations and Rural Development, unpublished data). In addition to these morphological differences between ecomorphs, the piscivores have significantly faster growth rates and higher standard metabolic rates relative to those of insectivores, and are behaviourally more bold as juveniles (e.g., took less time to explore novel habitat) [24]. This pattern of differentiation is consistent with the "pace-of-life syndrome" that has often been invoked to explain the divergence of closely related populations and species [25,26,27]. Thus, ecology can have a direct impact on the evolution of particular life-history strategies [28].
Study region of Kootenay Lake in southeastern British Columbia showing sampling locations, including the Blackwater River (shown in inset). See Additional file 1: Table S1 for more detailed sampling location information. Note that all individuals are sympatric in the main portion of Kootenay Lake during the non-breeding portion of the year
Previous research has shown that populations of these ecomorphs in Kootenay Lake are genetically distinct at microsatellite DNA loci (\(F_{\text{ST}}\) = 0.14; [76]), that no interbreeding is occurring, and that phenotypic differences between piscivorous and insectivorous forms show a high degree of heritability [20]. Together, these studies provide evidence against phenotypic plasticity, versus genetic mechanisms, as the source of this sympatric morphological polymorphism [29]. We therefore sought to answer three questions regarding the ecological and evolutionary divergence of the two O. mykiss forms in Kootenay Lake. First, what is the level of genomic divergence between piscivorous and insectivorous populations? Second, can we identify genomic regions responsible for the phenotypic divergence of these two forms? And third, is the phenotypic divergence between forms caused by few loci of large effect or many loci of small effect?
We addressed these questions using a whole-genome resequencing approach ("pool-seq"), which is a cost-effective technique for estimating population-level allelic frequencies and identifying loci putatively involved in functional genetic divergence of closely related populations [30]. This resulted in a dataset composed of millions of single nucleotide polymorphisms (SNPs) sampled from across the genome, with approximately 20× coverage per site. Based on the natural history of these populations, we hypothesized that elevated regions of genetic differentiation would harbour genes that fall into two classes: (1) those underlying metabolic and/or physiological processes due to body size and growth rate differences between the populations, and (2) genes associated with determining feeding behaviour and/or faciocranial morphology in relation to the divergent diet between populations (e.g., [31, 32]).
DNA sequence quality and genomic alignments
To address our hypotheses, we performed pool-seq on three populations—piscivores from Kootenay Lake, insectivores from Kootenay Lake, and insectivores from the Blackwater River (see Fig. 1 for map). Each of the three pools (n = 80 individuals per population) produced ~ 290 million raw read-pairs, with the majority (83–87%) passing quality filters (Table 1). The percentage of the genome covered following quality control and coverage thresholds was 62–63% with mean depth of coverages at 23–24× (Additional file 1: Fig. S1). Sequencing coverage varied across chromosomes, with the chromosomal ends often having higher sequencing coverage (Additional file 1: Fig. S2). Levels of genome-wide nucleotide diversity (\(\pi\)) also showed a similar pattern with elevated levels on ends of chromosomes (Additional file 1: Fig. S3), matching well to regions of known partial tetrasomy retained from the salmonid-specific whole-genome duplication event (Ss4rR; [33, 34]). In the dataset including all three populations, 8,140,802 SNPs met all filters and were included in analyses, whereas the Kootenay Lake piscivores-insectivores dataset had 6,896,554 SNPs. The density of SNPs varied across the genome, with areas of lower density occurring in centromeric regions (Additional file 1: Fig. S4).
Table 1 Sequence read data for the three groups of rainbow trout in this study
Population differentiation
We estimated the fixation index \(F_{\text{ST}}\) both on 100,000 bp non-overlapping windows and single SNPs to identify genomic regions putatively under selection and estimate the overall level of population differentiation between piscivorous and insectivorous ecomorphs. Two methods to estimate \(F_{\text{ST}}\) were used, and the Weir and Cockerham \(F_{\text{ST}}\) estimates were highly congruent with the PoPoolation2 ("Karlsson") results; 19 of the 20 most divergent windows shared and overall pairwise estimates very similar (results not shown). For consistency, only Weir and Cockerham \(F_{\text{ST}}\) estimates are presented. Our non-overlapping window \(F_{\text{ST}}\) analysis revealed high differentiation between all three populations (Fig. 2). The genome-wide divergence was higher between sympatric piscivores and insectivores (\(F_{\text{ST}}\) = 0.188) than between Kootenay Lake insectivores and Blackwater River insectivores (\(F_{\text{ST}}\) = 0.159), a population that is > 500 km (straight line) from Kootenay Lake. Although the modes of the Blackwater–Kootenay insectivore and Kootenay piscivore–Kootenay insectivore \(F_{\text{ST}}\) distributions are approximately the same (Fig. 2), the latter distribution is right-skewed, indicating a larger fraction of the genome has higher divergence in the sympatric population-pair. The Kootenay Lake piscivores and Blackwater River trout were the most divergent populations with an overall \(F_{\text{ST}}\) = 0.288.
Distributions of 100 kbp non-overlapping window Weir and Cockerham \(F_{\text{ST}}\) estimates, 19,512 in total, between all pairs of rainbow trout groups. Orange represents Blackwater insectivores–Kootenay Lake piscivores, blue is Blackwater insectivores–Kootenay Lake insectivores, and red is Kootenay Lake insectivores–Kootenay Lake piscivores. The table shows mean values of each comparison, which are drawn with dashed lines in the graph
Manhattan plot of 100 kbp non-overlapping window-based \(F_{\text{ST}}\) estimates between piscivorous and insectivorous rainbow trout from Kootenay Lake. The "outlier" cutoff line is shown in red, along with the windows that were above this value (0.609)
Highly differentiated genomic regions
Two methods were used to identify highly differentiated ("outlier") genomic regions: non-overlapping \(F_{\text{ST}}\) windows and the "local score" approach described by Fariello et al. [35]. Although single-SNP-based estimates of \(F_{\text{ST}}\) reached 1.0 in the empirical data (Additional file 1: Fig. S5), the maximal window-based \(F_{\text{ST}}\) estimate was 0.746 on chromosome 2 (Fig. 3). In general, the high overall level of genomic differentiation between piscivores and insectivores did not lead to obvious "islands of differentiation" surrounded by regions of low differentiation. The 20 most divergent windows (representing the most divergent 0.10% of all windows) all had \(F_{\text{ST}}\) estimates > 0.609, and were spread across 13 chromosomes (Additional file 1: Table S2).
To examine whether selection or recombination was driving signals of high genetic differentiation, we compared patterns of within-population nucleotide diversity (\(\pi\)) to absolute differentiation (d\(_{xy}\)) across four bins containing windows of increasing genetic divergence (\(F_{\text{ST}}\) > 0.2, 0.3, 0.4, and 0.5). In both piscivores and insectivores, \(\pi\) estimates were lower in the bins containing regions of higher genetic differentiation (Table 2). The decrease between the genome-wide and \(F_{\text{ST}}\) bin > 0.5 \(\pi\) estimates was ~ 73% in the piscivores, whereas the same decrease was only ~ 23% in the insectivores. No pattern existed, however, for d\(_{xy}\) across the bins of increasing genetic divergence (Table 2).
Table 2 Estimates of within population genetic diversity (\(\pi\)) and absolute genetic differentiation between populations (d\(_{xy}\)) of piscivores and insectivores by genomic bin
Table 3 Results from GOrilla gene ontology (GO) analysis of genes in the windows with \(F_{\text{ST}}\) estimates \(\ge\) 0.5, arranged by p-value (smallest to largest)
Given that the local score approach has not been assessed for accuracy and rates of false-positive detection (only power; see [35] for details), we performed simulations with parameters mimicking our empirical dataset, and chose the tuning parameter (\(\xi\)) to equal the 85% quantile of the − log\(_{10}\) p-value distribution (further details below in "Methods"; see Additional file 1: Figs. S6–S7 for p-value distributions of FLK and LK test statistics on simulated data). Assuming a chromosome-wide error correction rate of 0.05 (\(\alpha\)) and 30 chromosomes in a genome (the rainbow trout "version 1" reference genome has 29), the local score test would be predicted to incorrectly identify significant regions (e.g., false positives) on 1–2 chromosomes per genome. In our simulations of only neutral loci, ~ 6 regions per genome on average were identified as significant outliers (5.58 for FLK and 5.50 for LK; Additional file 1: Table S3), leading to a type I error rate of 19% for the FLK test and 18.3% for the LK test. Although FLK and LK are closely related statistics, approximately 7% of the outlier regions were significant based on the results of one test (either FLK or LK), but not the other.
According to a binomial test and false-positive rates calculated from the simulated data, we did not observe more significant local score results than expected by chance in either the FLK (p = 0.0533) or LK tests (p = 0.8094) in the empirical data, meaning that these numbers of significant genomic regions could be all false-positives. Due to the high false-positive rate from this test and weak power to differentiate neutral from "outlier" regions, we did not use the local score to identify outlier regions in our empirical data. We instead focus on divergent regions identified from the \(F_{\text{ST}}\) analysis.
Gene functions in highly differentiated genomic regions
After identifying highly differentiated genomic regions, we sought to determine the biological functions of genes showing strong divergence and determine if they support or refute our hypotheses. We determined biological processes of genes through gene ontology (GO) analyses, where we partitioned our data by level of genomic divergence (genes in windows with \(F_{\text{ST}}\) > 0.2, 0.3, 0.4, 0.5; and the 20 most divergent 100 kb windows). According to our hypotheses, we should see genes related to growth, metabolism, and/or faciocranial development become overrepresented in the more divergent genomic bins.
The rainbow trout genome has annotations for 42,884 protein-coding genes and 71,223 mRNAs. After filtering, our dataset had sequence data from 90% of all annotated genes (38,532) and 93% of all mRNA transcripts. The gene ontology analysis of genes (n = 47) in the 20 windows with the highest \(F_{\text{ST}}\) estimates did not produce any statistically significant results, likely due to too few genes included in the analysis. Gene ontology analyses of the four bins containing \(F_{\text{ST}}\) estimates > 0.2, 0.3, 0.4, and 0.5 resulted in 95, 72, 22, and 19 GO terms/processes, respectively (Fig. 4; Additional file 1: Figs. S8–S10). Whether the target set of genes was included in the background set made no difference on which GO terms were identified (results not shown). Although GO terms for behaviour and growth are present in the bin of least divergent loci (\(F_{\text{ST}}\) > 0.2; Additional file 1: Fig. S8), they are among 90+ other biological processes. However, we observed the process of "behaviour" in the \(F_{\text{ST}}\) > 0.4 bin (among 21 other processes; Additional file 1: Fig. S10), and more noticeably, growth and developmental growth in the bin of most divergent loci as two of 19 processes (Fig. 4; Table 3). Furthermore, the GO term for retinal metabolism, which includes the regulation of genes involved in development and growth, was also enriched in genes in the most divergent genomic regions. Other significantly enriched processes in the most divergent genomic regions (\(F_{\text{ST}}\) > 0.5) included blood vessel formation ("platelet-derived growth factor receptor signaling pathway"), cellular-level assembly of the cytoskeleton ("cytoskeleton organization"), and neuron cell-cell adhesion.
Results from a gene ontology (GO) enrichment analysis of loci from the most divergent 100 kbp windows (\(F_{\text{ST}}\) > 0.5, n = 187) between piscivorous and insectivorous rainbow trout from Kootenay Lake visualized as a "treemap" when the target gene set is not included in the background set. Processes are grouped by functional class (colour), and box size of each process is scaled by the significance (– log\(_{10}\) p-value; see Table 2) of enrichment of that process in relation to all processes in the genome. E.g., a larger tile size indicates a more significant result
We also searched annotations of the five most divergent \(F_{\text{ST}}\) windows, and serine/threonine-protein kinase BRSK2-like was found on chromosome two in the window with the highest \(F_{\text{ST}}\) estimate (Additional file 1: Table S4). This gene functions to regulate insulin secretion in fish, potentially playing a role in differential growth and development between ecomorphs. Lastly, based on gene annotations, we observed no evidence of the functional importance of genes associated with morphological feeding adaptations between these two populations.
Population divergence resulting from divergent selection on distinct phenotypes has been documented in many cases [36,37,38]. The genetic mechanisms underlying ecologically relevant traits, however, are generally not well understood. Here, through a dataset consisting of millions of genome-wide SNPs, we identified divergent loci involved with cell growth and differentiation that are associated with the divergence of two ecomorphs of rainbow trout—piscivores and insectivores. The genetic control of these ecomorphs does not seem to be controlled by a few loci of large effect, but rather by many genomic regions that are highly diverged between the populations. Our results further indicate that these two populations are moderately to highly divergent from one another with a genome-wide \(F_{\text{ST}}\) estimate of 0.18, suggesting little to no gene flow between these populations. These results are in agreement with previous research with microsatellite DNA loci that has shown that the two phenotypes are quite divergent from one another (\(F_{\text{ST}}\) = 0.14; [39]), and that these ecomorphs are largely genetically controlled [20].
Differential growth implicated in divergence of sympatric ecomorphs
Rainbow trout, as well as salmon (both Atlantic and Pacific species), are economically important species in the aquaculture industry that produces meat for human consumption. The Food and Agriculture Organization of the United Nations estimated an average annual production of ~ 800,000 tonnes and > $3 billion from world aquaculture production of O. mykiss alone [40]. Accordingly, much research has been done to determine the genetic mechanisms determining growth and muscle mass (e.g., "fillet weight") in these species (e.g., [41, 42]). Genetic association studies indicate that body weight in O. mykiss is a polygenic trait, with the most important genomic regions typically explaining < 5% of the genetic variance (e.g., [43, 44]). Two previous studies on O. mykiss have identified loci on chromosomes 21 [44] and 27 [43] that explain 2.5 and 1.7%, respectively, of the genetic variance in body weight. One of the 20 most divergent \(F_{\text{ST}}\) windows we identified on chromosome 21 is less 150 kb from the region identified by Neto et al. [44].
We hypothesized that genes either involved with growth and related metabolic processes or faciocranial morphology and feeding behaviours would underly the divergence of sympatric rainbow trout ecomorphs in Kootenay Lake. Identifying single SNPs or genes in the most divergent genomic regions as causal variants related to these traits is difficult given the high genome-wide divergence we estimated between insectivorous and piscivorous populations. Nonetheless, a gene that plays a role in the regulation of insulin secretion—serine/threonine-protein kinase BRSK2-like (100% similar to 987 bases of mRNA sequence)—was found in the most divergent \(F_{\text{ST}}\) window [45]. In teleost fishes, insulin is involved in somatic growth, reproduction, and development [46, 47], indicating that it could play a key role in differential growth between these two populations. This was the only gene of potentially large effect that we found in support of our hypothesis.
In addition to inspecting individual genes for their potential roles in population-level divergence, we also looked at sets of genes and their roles in emergent biological functions through gene ontology (GO) analyses. We hypothesized that biological functions and processes likely playing a role in the divergence of these ecomorphs would make up a greater proportion of the more strongly differentiated loci. In the least divergent bin of genomic regions that we analyzed (\(F_{\text{ST}}\) \(\ge\) 0.2), 238 GO terms were found spanning a variety of biological functions and processes (Additional file 1: Fig. S8). In contrast, only 19 significantly enriched GO terms occurred in the most divergent genomic regions (\(F_{\text{ST}}\) \(\ge\) 0.5), and growth functions accounted for more than half (10/19) of these GO terms. Interestingly, under the same diet, piscivores' growth and standing metabolic rates during ontogenetic development significantly increase relative to those of insectivores [24], corroborating our genetic results here. The gene ontology results indicate that the piscivorous phenotype is controlled by many loci of smaller effect and provide support to our hypothesis that growth plays a key role in leading to ecological divergence between populations.
Beyond the "growth" and "developmental growth" GO terms, gene regions associated with "retinal metabolism" were also significantly enriched in the most divergent \(F_{\text{ST}}\) bin. Retinoids (including retinal and other vitamin A derivatives) have two physiological functions—visual pigments in the eye and hormonal retinoids that regulate the expression of target genes involved in embryonic development, differentiation, and postnatal growth [48, 49]. Regulated amounts of retinoic acid and retinol are needed, particularly within mammalian embryonic tissues, for normal cell differentiation, proliferation, and morphogenesis [50]. Though less is known about the roles of these molecules in fishes, it is conceivable that they play analogous physiological roles and therefore could help explain some of the morphological differences between insectivores and piscivores.
Although not part of our search strategy, we identified via BLAST a combination of four transforming growth factor-\(\beta\) family genes that were present in multiple, un-linked chromosomal regions in the top 20 \(F_{\text{ST}}\) windows. We were not able to confirm, however, that these sequences indeed are from these genes (versus an unidentified repeating element in the rainbow trout genome) and therefore code for these growth proteins. Interestingly, a recent study examining the evolutionary origin of piscivory in brown trout (Salmo trutta) also identified transforming growth factor-\(\beta\) genes associated with the evolution of piscivory [7]. It is possible that the unidentified fragments in our BLAST searches represent paralogs or isoforms of proteins in the TGF-\(\beta\) superfamily. The common ancestor of all salmonids underwent a whole genome duplication event approximately 125 million years ago ("salmonid-specific fourth vertebrate whole-genome duplication, Ss4R"; [51]), in addition to a teleost-specific whole genome duplication event approximately 300–350 million years ago [52, 53]. Therefore, the current rainbow trout genome assembly may be lacking annotations of growth factor paralog isoforms that we have potentially identified in outlier regions. Genes controlling morphological traits are not well annotated as well, potentially explaining why we did not find any such genes significantly diverged between populations. Furthermore, this could be because our approach for detecting outliers is not optimal for detecting loci of small effect.
Evolutionary divergence of sympatric O. mykiss forms
Populations of piscivores and insectivores have been sympatric in Kootenay Lake for at most ~ 15,000 years, since the end of the last glacial maximum. It is not known whether these populations diverged in sympatry or allopatry and subsequently established sympatry, but it is assumed that they evolved in Kootenay Lake and did not originate there via human introduction (some human-mediated introductions of Gerrard trout have occurred to other areas for angling purposes). The high divergence we estimated between these populations (\(F_{\text{ST}}\) = 0.18) implies a long period of isolation, potentially pre-dating postglacial colonization. However, population divergence estimates between trout populations are often high, even for geographically proximate populations (e.g., strong IBD or IBE [isolation-by-environment]; [54]). Conversely, high genome-wide \(F_{\text{ST}}\) estimates can result from population bottlenecks/founder effects that cause reductions in genetic diversity in either or both populations. Although Kootenay Lake is large, the piscivorous population is small with an estimated census size of ~ 300–1000 spawning individuals recorded annually since 1961 [55]. Indeed, we observed on average lower levels of nucleotide diversity in piscivores than in insectivores across most of the genome (Additional file 1: Fig. S3). However, the higher levels of nucleotide diversity we observed in insectivores may be because this ecomorph was sampled from multiple populations around Kootenay Lake, which likely increased estimates of \(\pi\).
With a high overall genome-wide divergence and strong background differentiation between ecomorphs, such as we see in this system, finding clear candidate loci of large effect or "genomic islands" harbouring genes potentially involved in causing population divergence is difficult. In some cases, a suite of genes involved in local adaptation or divergence of phenotypes may be found linked together in a non-recombining block, for instance in an inverted chromosomal region (e.g., [8, 56]). However, we did not find any such candidate genomic regions. While trying to identify genomic outliers, using only \(F_{\text{ST}}\) estimates could be misleading given correlations between divergence and other features such as genetic diversity [57] and recombination [58]. Without a recombination map for wild rainbow trout, we are unable to directly account for this correlation in our analyses. We did, however, examine correlations between within-population genetic diversity (\(\pi\)) and absolute genetic differentiation between populations (d\(_{xy}\)), which help isolate the effects of selection from recombination. The patterns between \(\pi\) and d\(_{xy}\) (Table 2) provided further evidence for evolutionary recent selection driving the genomic divergence of piscivores from insectivores and act as further evidence against recombination being the main factor driving the overall pattern of genomic divergence between these ecomorphs.
Growth was one of the functions enriched in the most divergent loci (\(F_{\text{ST}}\)> 0.5) in our gene ontology (GO) analyses. This is consistent with the hypothesis stemming from the pace-of-life syndrome that closely related populations or species should differ in physiological traits that have co-evolved with the particular life histories of each population/species. This differentiation in growth-related genes could potentially underly the observation that piscivores have higher growth and standard metabolic rates than insectivores Monnet et al. [24]. Results from both our study and Monnet et al. highlight the importance of growth traits in differentiating these ecomorphs, thus producing a hypothesis regarding the functions of these genes that can be tested in an integrated physiological-genetic study.
Through the analysis of pool-seq data, we demonstrated that two feeding ecomorphs of rainbow trout that are sympatric for much of their life spans in Kootenay Lake are highly genetically differentiated from each other. The ecological and evolutionary distinctiveness of these ecomorphs—including contrasting diets, non-overlapping breeding sites, distinct morphologies, physiologies, and divergent genomes—highlights the critical role of managing them separately to promote their persistence into the future. Our results demonstrate an association between genes controlling growth and highly differentiated genomic regions, indicating a potential functional role of these loci in the divergence of these two rainbow trout feeding ecomorphs. Our results inform future studies that should perform transcriptomic analyses of these growth genes to evaluate gene expression during piscivorous O. mykiss development and the ontogenetic shift that accompanies a changing diet from primarily invertebrates to fish [23]. Our study further indicates the potential role of growth genes in causing phenotypic and ecological divergence of natural populations.
Sampling and DNA extraction
Fish were sampled from three separate populations: 80 piscivores from the Lardeau River (spawning site ~ 50 km north of Kootenay Lake, BC), eight insectivores each from 10 tributaries immediately surrounding Kootenay Lake (all 80 individuals combined together as a single "insectivore" population), and 80 insectivorous rainbow trout from the Blackwater River of central British Columbia (Fig. 1; Additional file 1: Table S1). The \(F_{\text{ST}}\) values estimated from microsatellite data between insectivorous populations range from 0.04 to 0.08, and all differ from the piscivorous population by \(F_{\text{ST}}\) = 0.11–0.28 [39]. The Blackwater River population was used as an outgroup for identifying significantly differentiated genomic regions (see below). Spawning adult piscivores were collected by angling in the spring, with some piscivore juveniles collected in the late summer/fall in the Lardeau River. Kootenay Lake insectivores were collected via electrofishing in the late summer/early fall. Fish were collected during the spawning season to facilitate population identification. Whole genomic DNA was extracted from fin tissues, or in a few cases, dried fish scales, with a Qiagen DNeasy extraction kit (Qiagen, Valencia, CA) and quantified with a Qubit fluorometer (Life Technologies, Carlsbad, CA). DNA from 80 individuals per population was combined into a single pool (for a total of three pools) in equimolar ratios for preparation of pool-seq libraries [30], then libraries were sequenced with 125 bp paired-end reads on an Illumina HiSeq 2500 at the Génome Québec Innovation Centre of McGill University.
Bioinformatic dataset assemblies
Raw sequence data (fastq format) were processed using the PoolParty pipeline [59] as follows. We first used FastQC [60] to assess raw read sequence quality. Raw paired-end reads were trimmed using the trim-fastq.pl script (part of PoPoolation v1; [61]) based on a Phred quality cutoff score of 20 and trimmed sequence length > 50 bp. We then used BWA [62] to align reads to the using the "version 1.0" O. mykiss reference genome (NCBI link), then Samblaster [63] and SamTools [64] were used to remove PCR duplicates and ambiguously aligned reads (e.g., reads with a low mapping quality), respectively. Reads were then sorted by coordinate and removed if unpaired in Picard Tools (v2.4.1; http://broadinstitute.github.io/picard/) and Samtools. Lastly, indel regions were then identified using the PoPoolation2 [65] identify-genomic-indel-regions.pl script, and SNPs were discarded if within 5 bp of the indel regions using the filter-sync-by-gtf.pl script in PoPoolation2. SNPs were retained that had (1) two alleles, (2) a minor allele frequency (MAF) cutoff > 0.05 globally (e.g., across all populations), (3) at least two copies of the minor allele at a locus (to mitigate the effect of sequencing errors), and (4) a sequencing depth of 15 \(\le\) x \(\le\) 50 within each population (data filtering script available in the Dryad repository). The minimum bound of 15× increases our confidence in SNP calls, while the upper bound of 50× was established to eliminate paralogs while retaining true homologs.
We used the fixation index (\(F_{\text{ST}}\)) to estimate population differentiation with a window-based approach. We calculated \(F_{\text{ST}}\) using two different approaches: the method of Weir and Cockerham [66] (R script available in the Dryad repository), and the "–karlsson" option within PoPoolation2. Both \(F_{\text{ST}}\) methods were used on two datasets, one containing only Kootenay Lake piscivores and insectivores for identifying outlier regions, and a second dataset containing the Kootenay Lake samples along with Blackwater River samples, to estimate pairwise genome-wide differentiation between the three populations. \(F_{\text{ST}}\) was calculated for 100,000 bp windows (along the reference genome) with a step-size of 100,000 bp (e.g., non-overlapping windows), giving ~ 19,500 \(F_{\text{ST}}\) estimates across the genome. Results by genomic position were visualized as a Manhattan plot in the R package qqman [67].
Identifying genomic regions of high differentiation
Two approaches were used to identify genomic regions that were highly differentiated between piscivores and insectivores and therefore potentially responsible in underlying their phenotypic divergence. The first method was with a window-based \(F_{\text{ST}}\) approach as described above (Weir and Cockerham [66] calculation, 100,000 bp non-overlapping windows). This is a reliable indicator of locus-specific estimates of population divergence, however, determining genomic outliers based solely on \(F_{\text{ST}}\) estimates can result in false-positives that may arise in neutral loci due to other processes such as recombination rate variation [58]. We focused on windows with the highest 0.1% of \(F_{\text{ST}}\) estimates (n = 20) as a first step in identifying the genomic regions that may be associated with the phenotypic differentiation of these two ecomorphs. Non-overlapping windows were used because a sliding-window approach leads to non-independent tests from overlapping windows with no clear statistical adjustment for multiple comparisons.
Outlier detection via \(F_{\text{ST}}\)-based approaches may lead to false-positives when identifying putative regions involved in local adaptation because \(F_{\text{ST}}\) estimates can reflect patterns of recombination rather than selection [58]. In lieu of a recombination map, patterns of nucleotide diversity within populations (\(\pi\)) and absolute differentiation between two populations (d\(_{xy}\), which is independent of levels of genetic diversity within each population) can provide further evidence for or against selection in highly differentiated genomic regions. For instance, when a beneficial mutation arises in a population, \(\pi\) will decrease in this region as the mutation sweeps through the population. If the beneficial mutation has recently evolved, d\(_{xy}\) will not be elevated in selected regions. However, d\(_{xy}\) will be elevated in regions of low \(\pi\) if the beneficial mutation is "ancient" (e.g., more than 10\(^{5}\) generations, given a mutation rate of ~ 10\(^{-8}\)). We therefore compared patterns of \(\pi\) and d\(_{xy}\) in bins of genomic windows ranked by genetic differentiation, with windows divided into four \(F_{\text{ST}}\) bins (\(F_{\text{ST}}\) > 0.2, 0.3, 0.4, and 0.5; same as GO analyses described below), to help distinguish between the evolutionary processes (e.g., recombination or selection) affecting highly differentiated genomic regions between insectivores and piscivores. The median \(\pi\) value was calculated within each population in 100 kb windows in NPStat [68, 69], and the median d\(_{xy}\) was calculated in a custom R script included in the Dryad data repository.
As a complementary approach to outlier identification via \(F_{\text{ST}}\) windows, we also used the "local score" approach developed by Fariello et al. [35]. The local score method identifies genomic regions of high differentiation between populations by using a Lindley process to identify stretches of adjacent/linked SNPs that show significant differentiation as determined by single-locus analyses. This method has the benefit over a fixed window approach in that the outlier region bounds are not user-specified and are instead determined during analysis. As in the \(F_{\text{ST}}\) window-based analysis, this approach takes advantage of linkage by using divergence information from adjacent sites and is therefore superior to single-locus outlier detection methods. The local score method implements a chromosome-wide error rate correction, meaning \(\alpha\) fraction (e.g., 0.05) of the chromosomes in the genome will produce false positives. We used the significance of LK (Lewontin–Krakauer) and FLK (an extension of the LK test using a population kinship matrix "F"; [70]) single-locus test results expressed as p-values that were calculated in R scripts provided by Bonhomme et al. [70], as the input for the local score analysis. We considered loci as outliers at the 0.05 significance level.
The value of the "tuning parameter" (\(\xi\)) is important in local score analyses, which determines the threshold for identifying regions as outliers. The local score approach assumes a uniform distribution of p-values (from either the FLK or LK test statistic). However, in the case of a non-uniform distribution of p-values, \(\xi\) must be adjusted accordingly. Fariello et al. [35] state that \(\xi\) must be between the mean and maximum values of the − log\(_{10}\) p-value distribution, and we explored various values in this range. \(\xi\) values close to the max(− log\(_{10}\) p-value distribution) resulted in no outliers identified, and \(\xi\) values close to mean(− log\(_{10}\) p-value distribution) produced hundreds of significant outlier regions (results not shown). We therefore set \(\xi\) in all analyses to equal the 85% quantile of the − log\(_{10}\) p-value distribution.
One potential downside to the local score approach is that the true type I error rate (e.g., false-positive rate) has not been previously reported. We therefore estimated it through simulations mimicking our empirical dataset as follows. A single ancestral population of 1000 individuals evolved for 1000 generations and then split into three populations, each consisting of 1000 individuals. These three populations evolved independently (e.g., no migration) for 2500 more generations (generating an \(F_{\text{ST}}\) estimate of ~ 0.18 between all three populations) at which point an average of 23 chromosomes (from a truncated distribution of 15–50 and \(\sigma\) = 5) were sampled from each population, based on the average sequencing depth per population from the empirical data (see "Results"). Each chromosome was 75 MB in length and simulated in SLiM v3 [71, 72] with a recombination rate of 1.5 × 10\(^{-8}\) (chromosome was ~ 1 Morgan in length) and mutation rate of 2 × 10\(^{-7}\) to generate ~ 15 million SNPs per genome after implementing a MAF filter (> 0.05); only neutral mutations were introduced. Thirty chromosomes were then combined into a "pseudo-genome" replicate to mimic the empirical dataset (29 chromosomes), from which allelic frequencies were calculated that were then used in FLK and LK analyses using the R scripts provided in Bonhomme et al. [70]. Local score tests were then performed on the FLK and LK results of each genomic replicate using the R scripts provided by Fariello et al. [35] and results were examined at the 0.05 significance level, per genome. One hundred pseudo-genome replicates were simulated and analyzed (e.g., 3000 chromosomal replicates). To determine if significantly more genomic regions were identified in our empirical data than expected given the false discovery rate estimated in our simulations (see "Results" below), we ran an exact binomial test in R with significance determined at the p \(\le\) 0.05 level. The null hypothesis is that we expect as many false positives in the empirical data as are found in the simulations.
Functions of divergent loci
We performed gene ontology (GO) enrichment analyses to identify putative functions of genes in the most divergent genomic regions between piscivores and insectivores. The GO database is designed to be species-neutral, and therefore annotations were transferred from Homo sapiens to O. mykiss [73]. We looked at gene function in two sets of divergent loci. The first set contained the genes in the 20 most divergent \(F_{\text{ST}}\) windows (n genes = 45). And secondly, by separating our data into four bins based on estimates from our window-based \(F_{\text{ST}}\) analyses: windows with \(F_{\text{ST}}\) estimates > 0.2, 0.3, 0.4, and 0.5. These bins represent the 37% (n windows = 7219), 14% (2,737), 4.4% (854), and 0.9% (187) most divergent regions in the genome, respectively. We organized our data into bins based on divergence estimates because we expected the loci associated with biological processes most likely involved in the divergence of these populations to emerge in the more highly differentiated genomic regions (e.g., the genomic bins with higher \(F_{\text{ST}}\) estimates) as compared to the less differentiated genomic regions (genomic bins with lower \(F_{\text{ST}}\) estimates). According to our hypothesis, the most divergent genomic regions would contain genes enriched for processes related to growth and metabolism.
We used bedtools [74] to extract gene and protein annotations from the available rainbow trout reference genome. We used the GOrilla online platform [75, 76] to identify significantly enriched biological processes in differentiated genomic regions with two slightly different methods: (1) including the target gene set (all loci in an \(F_{\text{ST}}\) bin) in the background set (all genes in the rainbow trout genome), vs. (2) excluding the target gene set from the background set. GOrilla uses a minimum hypergeometric (mHG) score to assign the significance of a term occurring in the target set vs. background set of genes. The mHG score reflects the surprise of seeing a particular GO term in the target set compared to its probability of occurrence in the background set, under the null assumption that all GO terms in the background set occur with equal probability. An exact p-value of this score corrected for multiple testing is then calculated. All GO analyses were performed with a significance threshold of p \(\le\) 0.001 (as recommended by the authors). Gene ontology analysis results were visualized with Revigo [77] as "treemaps", where the size of each functional category was scaled by its – log\(_{10}\) p-value (e.g., a smaller p-value equals a larger tile size). Beyond the GO analyses, we also used annotations of genes in the five most divergent \(F_{\text{ST}}\) windows to identify plausible candidates underlying the phenotypic divergence of insectivorous and piscivorous forms.
Data availibility statement
Raw, un-edited DNA Sequences are available in NCBI's Short Read Archive, BioProject ID PRJNA695019. Sampling locality information is available in Additional file 1: Table S1. Scripts for analysis are available in the Dryad repository (Dryad link; https://datadryad.org/stash/share/jYk11SUp6azAZKMNbbdwt3JRZLe-fVdtmLcTby0vbTQ).
Coyne J, Orr H. Speciation. Sunderland: MA p; 2004. p. 281.
Kitano J, Mori S, Peichel CL. Phenotypic divergence and reproductive isolation between sympatric forms of Japanese threespine sticklebacks. Biol J Linnean Soc. 2007;91:671–85.
Barluenga M, Stölting KN, Salzburger W, Muschick M, Meyer A. Sympatric speciation in nicaraguan crater lake cichlid fish. Nature. 2006;439:719–23.
Morris J, Navarro N, Rastas P, Rawlins LD, Sammy J, Mallet J, Dasmahapatra KK. The genetic architecture of adaptation: convergence and pleiotropy in heliconius wing pattern evolution. Heredity. 2019;123:138–52.
Supple MA, Hines HM, Dasmahapatra KK, Lewis JJ, Nielsen DM, Lavoie C, Ray DA, Salazar C, McMillan WO, Counterman BA. Genomic architecture of adaptive color pattern divergence and convergence in heliconius butterflies. Genome Res. 2013;23:1248–57.
Boyle EA, Li YI, Pritchard JK. An expanded view of complex traits: from polygenic to omnigenic. Cell. 2017;169:1177–86.
Jacobs A, Hughes MR, Robinson PC, Adams CE, Elmer KR. The genetic architecture underlying the evolution of a rare piscivorous life history form in brown trout after secondary contact and strong introgression. Genes. 2018;9:280.
Article PubMed Central CAS Google Scholar
Kirubakaran TG, Grove H, Kent MP, Sandve SR, Baranski M, Nome T, De Rosa MC, Righino B, Johansen T, Otterå H, et al. Two adjacent inversions maintain genomic differentiation between migratory and stationary ecotypes of Atlantic cod. Mol Ecol. 2016;25:2130–43.
Prince DJ, O'Rourke SM, Thompson TQ, Ali OA, Lyman HS, Saglam IK, Hotaling TJ, Spidle AP, Miller MR. The evolutionary basis of premature migration in pacific salmon highlights the utility of genomics for informing conservation. Sci Adv. 2017;3:e1603198.
Mackay TF. The genetic architecture of quantitative traits. Ann Rev Genet. 2001;35:303–39.
Robertson A. The nature of quantitative genetic variation. In: Heritage from Mendel, 1967. pp. 265–80.
Sinclair-Waters M, Ødegård J, Korsvoll SA, Moen T, Lien S, Primmer CR, Barson NJ. Beyond large-effect loci: large-scale GWAS reveals a mixed large-effect and polygenic architecture for age at maturity of Atlantic salmon. Genet Sel Evol. 2020;52:1–11.
Allen HL, Estrada K, Lettre G, Berndt SI, Weedon MN, Rivadeneira F, Willer CJ, Jackson AU, Vedantam S, Raychaudhuri S, et al. Hundreds of variants clustered in genomic loci and biological pathways affect human height. Nature. 2010;467:832–8.
Kardos M, Husby A, McFarlane SE, Qvarnström A, Ellegren H. Whole-genome resequencing of extreme phenotypes in collared flycatchers highlights the difficulty of detecting quantitative trait loci in natural populations. Mol Ecol Resour. 2016;16:727–41.
Jonsson B, Jonsson N. Polymorphism and speciation in arctic charr. J Fish Biol. 2001;58:605–38.
Skúlason S. Sympatric morphs, populations and speciation in freshwater fish with emphasis on arctic charr. Evol Biol Divers. 1999:71–92.
Skúlason S, Snorrason SS, Ota D, Noakes DL. Genetically based differences in foraging behaviour among sympatric morphs of arctic charr (pisces: Salmonidae). Anim Behav. 1993;45:1179–92.
Hughes MR, Hooker OE, Van Leeuwen TE, Kettle-White A, Thorne A, Prodöhl P, Adams CE. Alternative routes to piscivory: contrasting growth trajectories in brown trout (Salmo trutta) ecotypes exhibiting contrasting life history strategies. Ecol Freshw Fish. 2019;28:4–10.
Thorne A, MacDonald AI, Thorley JL. The abundance of large, piscivorous ferox trout (Salmo trutta) in loch rannoch, scotland. PeerJ. 2016;4:e2646.
Keeley E, Parkinson E, Taylor E. The origins of ecotypic variation of rainbow trout: a test of environmental vs. genetically based differences in morphology. J Evol Biol. 2007;20:725–36.
Hess JE, Zendt JS, Matala AR, Narum SR. Genetic basis of adult migration timing in anadromous steelhead discovered through multivariate association testing. Proc R Soc B Biol Sci. 2016;283:20153064.
Irvine JR. The Gerrard rainbow trout of Kootenay lake, British Columbia: a discussion of their life history with management, research and enhancement recommendations. BC Hydro Power Authority Fish Management Report. 1978;72:1–62.
Northcote T. Some impacts of man on kootenay lake and its salmonoids. Great Lakes Fish Comm. 1973;25:1–46.
Monnet G, Rosenfeld JS, Richards JG. Adaptive differentiation of growth, energetics and behaviour between piscivore and insectivore juvenile rainbow trout along the pace-of-life continuum. J Anim Ecol. 2020;89:2717–32.
Réale D, Garant D, Humphries MM, Bergeron P, Careau V, Montiglio P-O. Personality and the emergence of the pace-of-life syndrome concept at the population level. Philos Trans R Soc B Biol Sci. 2010;365:4051–63.
Ricklefs RE, Wikelski M. The physiology/life-history nexus. Trends Ecol Evol. 2002;17:462–8.
Wikelski ML, Spinney L, Schelsky W, Scheuerlein A, Gwinner E. Slow pace of life in tropical sedentary birds: a common-garden experiment on four stonechat populations from different latitudes. Proc R Soc Lond Ser B Biol Sci. 2003;270:2383–8.
Roff D. Evolution of life histories: theory and analysis. Berlin: Springer Science & Business Media; 1993.
Skulason S, Smith TB. Resource polymorphisms in vertebrates. Trends Ecol Evol. 1995;10:366–70.
Futschik A, Schlötterer C. The next generation of molecular markers from massively parallel sequencing of pooled DNA samples. Genetics. 2010;186:207–18.
McGirr JA, Martin CH. Novel candidate genes underlying extreme trophic specialization in caribbean pupfishes. Mol Biol Evol. 2017;34:873–88.
Roberts RB, Hu Y, Albertson RC, Kocher TD. Craniofacial divergence and ongoing adaptation via the hedgehog pathway. Proc Natl Acad Sci. 2011;108:13194–9.
Berthelot C, Brunet F, Chalopin D, Juanchich A, Bernard M, Noël B, Bento P, Da Silva C, Labadie K, Alberti A, et al. The rainbow trout genome provides novel insights into evolution after whole-genome duplication in vertebrates. Nat Commun. 2014;5:1–10.
Campbell MA, Hale MC, McKinney GJ, Nichols KM, Pearse DE. Long-term conservation of ohnologs through partial tetrasomy following whole-genome duplication in salmonidae. G3: Genes Genomes Genet. 2019;9:2017–28.
Fariello MI, Boitard S, Mercier S, Robelin D, Faraut T, Arnould C, Recoquillay J, Bouchez O, Salin G, Dehais P, et al. Accounting for linkage disequilibrium in genome scans for selection without individual genotypes: the local score approach. Mol Ecol. 2017;26:3700–14.
McKinnon JS, Mori S, Blackman BK, David L, Kingsley DM, Jamieson L, Chou J, Schluter D. Evidence for ecology's role in speciation. Nature. 2004;429:294–8.
Schluter D, McPhail JD. Ecological character displacement and speciation in sticklebacks. Am Nat. 1992;140:85–108.
Skúlason S, Parsons KJ, Svanbäck R, Räsänen K, Ferguson MM, Adams CE, Amundsen P-A, Bartels P, Bean CW, Boughman JW, et al. A way forward with eco evo devo: an extended theory of resource polymorphism with postglacial fishes as model systems. Biol Rev. 2019;94:1786–808.
Taylor EB, Foley C, Neufeld M. Genetic mixture analyses in support of restoration of a high value recreational fishery for rainbow trout (Oncorhynchus mykiss) from a large lake in interior British Columbia. Conserv Genet. 2019;20:891–902.
Fisheries F, Fishery and aquaculture statistics 2017/fao annuaire. statistiques des pêches et de l'aquaculture 2017/fao anuario. estadísticas de pesca y acuicultura 2017. Global aquaculture production statistics for the year. 2017.
Salem M, Vallejo RL, Leeds TD, Palti Y, Liu S, Sabbagh A, Rexroad CE III, Yao J. RNA-seq identifies SNP markers for growth traits in rainbow trout. PLoS ONE. 2012;7:e36264.
Yoshida GM, Lhorente JP, Carvalheiro R, Yáñez J. Bayesian genome-wide association analysis for body weight in farmed Atlantic salmon (Salmo salar l.). Anim Genet. 2017;48:698–703.
Gonzalez-Pena D, Gao G, Baranski M, Moen T, Cleveland BM, Kenney PB, Vallejo RL, Palti Y, Leeds TD. Genome-wide association study for identifying loci that affect fillet yield, carcass, and body weight traits in rainbow trout (oncorhynchus mykiss). Front Genet. 2016;7:203.
Neto RVR, Yoshida GM, Lhorente JP, Yáñez JM. Genome-wide association analysis for body weight identifies candidate genes related to development and metabolism in rainbow trout (oncorhynchus mykiss). Mol Genet Genom. 2019;294:563–71.
Chen X-Y, Gu X-T, Saiyin H, Wan B, Zhang Y-J, Li J, Wang Y-L, Gao R, Wang Y-F, Dong W-P, et al. Brain-selective kinase 2 (brsk2) phosphorylation on pctaire1 negatively regulates glucose-stimulated insulin secretion in pancreatic β-cells. J Biol Chem. 2012;287:30368–75.
Garcia de la Serrana D, Macqueen DJ. Insulin-like growth factor-binding proteins of teleost fishes. Front Endocrinol. 2018;9:80.
Wood AW, Duan G, Bern HA. Insulin-like growth factor signaling in fish. Int Rev Cytol. 2005;243:215–85.
Mangelsdorf DJ. Vitamin a receptors. Nutr Rev. 1994;52:S32.
Wolf G. Multiple functions of vitamin a. Physiol Rev. 1984;64:873–937.
Morriss-Kay GM. Retinoids in mammalian embryonic development. In: Advances in organ biology, vol. 3. Elsevier. 1997. pp. 79–92.
Macqueen DJ, Johnston IA. A well-constrained estimate for the timing of the salmonid whole genome duplication reveals major decoupling from species diversification. Proc R Soc B Biol Sci. 2014;281:20132881.
Near TJ, Eytan RI, Dornburg A, Kuhn KL, Moore JA, Davis MP, Wainwright PC, Friedman M, Smith WL. Resolution of ray-finned fish phylogeny and timing of diversification. Proc Natl Acad Sci. 2012;109:13698–703.
Vandepoele K, De Vos W, Taylor JS, Meyer A, Van de Peer Y. Major events in the genome evolution of vertebrates: paranome age and size differ considerably between ray-finned fishes and land vertebrates. Proc Natl Acad Sci. 2004;101:1638–43.
Tamkee P, Parkinson E, Taylor E. The influence of wisconsinan glaciation and contemporary stream hydrology on microsatellite DNA variation in rainbow trout (Oncorhynchus mykiss). Can J Fish Aquat Sci. 2010;67:919–35.
Ltd RC. Kootenay lake action plan-2016. Nelson: Lands and Natural Resource Operations, The Ministry of Forests; 2016.
Pearse DE, Barson NJ, Nome T, Gao G, Campbell MA, Abadía-Cardoso A, Anderson EC, Rundio DE, Williams TH, Naish KA, et al. Sex-dependent dominance maintains migration supergene in rainbow trout. Nat Ecol Evol. 2019;3:1731–42.
Cruickshank TE, Hahn MW. Reanalysis suggests that genomic islands of speciation are due to reduced diversity, not reduced gene flow. Mol Ecol. 2014;23:3133–57.
Booker TR, Yeaman S, Whitlock M. Variation in recombination rate affects detection of outliers in genome scans under neutrality. Mol Ecol. 2020;29:4274–9.
Micheletti SJ, Narum SR. Utility of pooled sequencing for association mapping in nonmodel organisms. Mol Ecol Resour. 2018;18:825–37.
Andrews S. Fastqc: a quality control tool for high throughput sequence data. 2010.
Kofler R, Orozco-terWengel P, De Maio N, Pandey RV, Nolte V, Futschik A, Kosiol C, Schlötterer C. Popoolation: a toolbox for population genetic analysis of next generation sequencing data from pooled individuals. PloS ONE. 2011a;6:e15925.
Li H, Durbin R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics. 2009;25:1754–60.
Faust GG, Hall IM. Samblaster: fast duplicate marking and structural variant read extraction. Bioinformatics. 2014;30:2503–5. https://doi.org/10.1093/bioinformatics/btu314.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, et al. The sequence alignment/map format and samtools. Bioinformatics. 2009;25:2078–9.
Article PubMed PubMed Central CAS Google Scholar
Kofler R, Pandey RV, Schlötterer C. Popoolation2: identifying differentiation between populations using sequencing of pooled DNA samples (pool-seq). Bioinformatics. 2011b;27:3435–6.
Weir BS, Cockerham CC. Estimating f-statistics for the analysis of population structure. Evolution. 1984:1358–70.
Turner S. qqman: Q-Q and Manhattan Plots for GWAS Data. R package version 0.1.4. 2017.https://CRAN.R-project.org/package=qqman.
Ferretti L, Ramos-Onsins SE, Pérez-Enciso M. Population genomics from pool sequencing. Mol Ecol. 2013;22:5561–76.
Fracassetti M, Griffin PC, Willi Y. Validation of pooled whole-genome re-sequencing in Arabidopsis lyrata. PLoS ONE. 2015;10:e0140462.
Bonhomme M, Chevalet C, Servin B, Boitard S, Abdallah J, Blott S, SanCristobal M. Detecting selection in population trees: the Lewontin and Krakauer test extended. Genetics. 2010;186:241–62.
Haller BC, Galloway J, Kelleher J, Messer PW, Ralph PL. Tree-sequence recording in slim opens new horizons for forward-time simulation of whole genomes. Mol Ecol Resour. 2019;19:552–66.
Haller BC, Messer PW. Slim 3: Forward genetic simulations beyond the wright-fisher model. Mol Biol Evol. 2019;36:632–7.
Primmer C, Papakostas S, Leder E, Davis M, Ragan M. Annotated genes and nonannotated genomes: cross-species use of gene ontology in ecology and evolution research. Mol Ecol. 2013;22:3216–41.
Quinlan AR, Hall IM. Bedtools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841–2.
Eden E, Lipson D, Yogev S, Yakhini Z. Discovering motifs in ranked lists of DNA sequences. PLoS Comput Biol. 2007;3:e39.
Eden E, Navon R, Steinfeld I, Lipson D, Yakhini Z. Gorilla: a tool for discovery and visualization of enriched go terms in ranked gene lists. BMC Bioinform. 2009;10:48.
Supek F, Bošnjak M, Škunca N, Šmuc T. Revigo summarizes and visualizes long lists of gene ontology terms. PloS ONE. 2011;6:e21800.
We thank S. Narum and the Columbia River Inter-Tribal Fish Commission for DNA library preparations, S. Northrup and the Freshwater Fisheries Society of BS for access to tissues, and S. Micheletti for his troubleshooting assistance with PoolParty. We would also like to give a big thanks T. Booker for his assistance with the Local Score simulations in SLiM.
The study was funded by a Genome Canada Large-Scale Applied Research Project in Natural Resources and the Environment (Project code 242RTE), NSERC Grant no. 81696 to EBT, and NSERC Grant no. RGPIN-2016-03779 to MCW.
Department of Zoology and Biodiversity Research Centre, University of British Columbia, 6270 University Blvd., Vancouver, BC, V6T 1Z4, Canada
Jared A. Grummer, Michael C. Whitlock, Patricia M. Schulte & Eric B. Taylor
Beaty Biodiversity Museum, University of British Columbia, 6270 University Blvd., Vancouver, BC, V6T 1Z4, Canada
Eric B. Taylor
Jared A. Grummer
Michael C. Whitlock
Patricia M. Schulte
JAG designed the study, developed analysis scripts, performed analyses, and wrote the manuscript. MCW, PMS, and EBT designed the study and edited the manuscript. All authors read and approved the final manuscript.
Correspondence to Jared A. Grummer.
All sampling was conducted under scientific permit and was approved by the University of British Columbia Animal Care and Use committee (UBC-AC) under Protocol number A16-0329. Permits issued under the UBC-AC are reviewed by expert committees to ensure ethical treatment of animals and were performed under the Canadian Committee on Animal Care (CCAC) Good Animal Practice certificate. This certificate demonstrates full adherence to CCAC's internationally recognized standards of animal ethics and care in research. This included collection of non-lethal samples and precautions to reduce stress and ensure adequate recovery after release.
Additional tables and figures.
Grummer, J.A., Whitlock, M.C., Schulte, P.M. et al. Growth genes are implicated in the evolutionary divergence of sympatric piscivorous and insectivorous rainbow trout (Oncorhynchus mykiss). BMC Ecol Evo 21, 63 (2021). https://doi.org/10.1186/s12862-021-01795-9
Genome scan
Small-effect loci
\(F_{\text{ST}}\)
GO analysis
Ecomorph
Submission enquiries: bmcecologyandevolution@biomedcentral.com


text:
Literature Meta
Literature Stack Exchange
Answer a question more than 60 days later with score of 5 or more. This badge can be awarded multiple times.
Awarded 242 times.
What evidence is there for the "recession of accent" theory?
Awarded 2 days ago
Gareth Rees
What is the origin of this quote about God, Flies, and Spiders?
Awarded Jan 16 at 14:43
What does "fail us" mean in the book "Pachinko"?
What is known about the content of the Ur-Hamlet?
Awarded Dec 20, 2022 at 18:41
Tsundoku♦
Ray Bradbury says Fahrenheit 451 isn't about censorship. Is he right?
spozun
When did detective fiction become primarily about murder?
Lewis 010
In Great Expectations, who is the man at the pub in Chapter Ten?
Spagirl
16.6k4949 silver badges103103 bronze badges
In Macbeth, why is Fleance 'scaped?
Awarded Dec 5, 2022 at 0:16
Who introduced the term close reading in the context of literature?
Awarded Nov 28, 2022 at 15:17
Do we know what prompted Orwell to write 1984?
Awarded Nov 28, 2022 at 6:09
Faheem Mitha
Why is Georg Büchner considered such an important figure in German literature?
Awarded Nov 6, 2022 at 17:30
What is the name of a poem about relativity?
Awarded Oct 24, 2022 at 11:25
Clara Diaz Sanchez
What is the history of the universe/canon/word-of-god approach to literature?
Book about a man who entered a relationship with a shady woman and then had to kill people
Awarded Oct 20, 2022 at 8:27
What does Borges' reference to the Greek mean in "The Lottery in Babylon"?
Why is black evil and white good?
Chris Sunami
Where did the name Nimitseahpah come from?
Awarded Oct 9, 2022 at 15:18
Why does the Phantom call both Christina and himself "Angel of Music?"
Awarded Sep 30, 2022 at 10:11
DVK
Whose fingerprint is on the cover of "Sapiens: A Brief History of Humankind"?
Mithical♦
21k1616 gold badges8080 silver badges157157 bronze badges
What was the ethnicity of Heathcliff in Wuthering Heights?
Short story about a condemned organ
Quassnoi
If Tolkien disliked allegory, why are there so many allegorical readings of The Lord of the Rings?
Connor Jones
What does "atom" mean in Don Quixote?
Awarded Sep 11, 2022 at 6:40
Do the *books* show that Hogwarts is in Scotland?
Awarded Aug 20, 2022 at 13:02
Rand al'Thor♦
Medieval England story - "Step in, Step Inn", "Good Morning, Good Knight"
Awarded Aug 8, 2022 at 16:47
Who said "Poetry is the art of giving different names to the same thing"?
Meaning of Lorca's "Romance sonámbulo"?
Awarded Jul 30, 2022 at 11:19
b4rtr
What is the meaning of "Not even Jonah's wife would swallow that story" in Marquez's "Love in the Time of Cholera"?
What did Gray mean by "where ignorance is bliss, / 'tis folly to be wise"?
Awarded Jun 21, 2022 at 12:48
DukeZhou
How old is Romeo?
How much time did one canto in Dante's Divine Comedy represent (if any)?
Awarded Jun 5, 2022 at 22:25
WalterVi
Earliest work of fiction in which characters using telepathy can't lie
Awarded May 18, 2022 at 7:08
What's with the beyond-casual introduction of "stone-giants" in The Hobbit?
Auden Young
Significance of the Pardoner's hair style
What is the origin of epigraphs?
Meaning of the title "They Shoot Horses, Don't They?"
Chenmunka
Awarded Apr 1, 2022 at 0:37
Ralph Crown
What was the first book that explicitly used "A Hero's Journey" as a checklist?
Awarded Mar 28, 2022 at 8:03
osterzone
What is the Neverending Story's moral about power?
Awarded Mar 26, 2022 at 20:25
In Huxley's Island, what is a "bass thurible?"
Valorum
Meaning of "I shut my tongue against my fly"
What does it mean to "floor peculiarities parabolous" in conics, from The Pirates of Penzance?
Why does Tom Sawyer forget what "ransom" means?
Awarded Mar 3, 2022 at 3:10
Use of 'Genius' in Nicholas Nickleby?
Awarded Feb 23, 2022 at 3:45
Which traits of Milton's Lucifer from "Paradise Lost" did Neil Gaiman carry to "The Sandman"?
Awarded Feb 20, 2022 at 21:40
What did Scout mean when she said "It takes a woman to do that kind of work"?
Awarded Feb 8, 2022 at 17:05
Reference for author of "The earth has music for those who listen"
Awarded Jan 30, 2022 at 6:56
Is John Milton's Lucifer a tragic hero?
Awarded Jan 27, 2022 at 14:57
Kiteration
What classic mystery novels and stories led to "the butler did it" becoming a cliché?
Term for intentionally impressing different meanings on different audiences
Has a parody of a work of literature ever become more successful than the original work?
Awarded Jan 4, 2022 at 13:08
How can Gilgamesh be one-third man and two-thirds god?
Why does the Iliad start "in the middle"?
Awarded Jan 3, 2022 at 0:11
How did Ibsen's writing in A Doll's House influence the James Joyce character Molly Bloom?
fundagain
Which flood story was first: Genesis or The Epic of Gilgamesh?
Did the recipient receive twelve partridges or just one?
What Confederate counterfactual story features Florida as an independent tourist republic?
Awarded Dec 2, 2021 at 18:50
What's the evidence for "Dr Jekyll and Mr Hyde" being secretly set in Edinburgh?
John Strachan
What is the historical basis for Roark's claim about the Parthenon?
Cahir Mawr Dyffryn æp Ceallach
Did Thoreau respond to claims that he was an anarchist?
HoldOffHunger


text:
Q: React axios api call with user input I'm very new to react and I am trying to make an api call based on user-input using axios. I'm able to get the results with the hard-coded query('cars') on form submit, but I'd like to be able to allow the user to search different items. I'm stuck on how to connect the user input to the api call.
import axios from "axios";

import { useState, useEffect } from "react";

function App() {
  const [data, setData] = useState([]);
  const [userInput, setUserInput] = useState([]);

  const handleSubmit = (event) => {
    event.preventDefault();
    axios({
      url: "https://api.unsplash.com/search/photos",
      method: "GET",
      dataResponse: "json",
      params: {
        client_id: "d54vE8fu6bjJ_JgkBqugaZOt4bwFHGkmiKDGHunnXxc",
        query: "cars",
        per_page: 10,
      },
    }).then((response) => {
      console.log(response.data.results);
      setData(response.data.results);
    });
  };

  return (

      <form action="" onSubmit={handleSubmit}>
        Search
        <input
          id="search"
          type="text"
          value={userInput}
          onChange={(event) => setUserInput(event.target.value)}
        />
      </form>

        {data.map((item) => {
          return (

          );
        })}

  );
}

export default App;

A: You can put query = userInput
const handleSubmit = (event) => {
    event.preventDefault();
    axios({
      url: "https://api.unsplash.com/search/photos",
      method: "GET",
      dataResponse: "json",
      params: {
        client_id: "d54vE8fu6bjJ_JgkBqugaZOt4bwFHGkmiKDGHunnXxc",
        query: userInput,
        per_page: 10,
      },
    }).then((response) => {
      console.log(response.data.results);
      setData(response.data.results);
    });
  };

redpajama_set_name: RedPajamaStackExchange


text:
Home > News > Union Minister Inaugurates Laboratory Complex & NVBDCP Building In Delhi
Union Minister Inaugurates Laboratory Complex & NVBDCP Building In Delhi
2019-08-01 Announcements Other Services Northern India
"National Centre for Disease Control (NCDC) has contributed phenomenally in elimination and eradication of several diseases of public health importance such as smallpox, polio, guinea-worm and yaws." This was stated by Dr Harsh Vardhan, Union Minister of Health and Family Welfare at the 110th Annual Day celebrations of NCDC.
At the function, the Union Health Minister inaugurated the laboratory complex (Lab 3) at NCDC and the new building of National Vector Borne Disease Control Programme (NVBDCP) in presence of Shri Ashwini Kumar Choubey, MoS (HFW). The L3 laboratory complex inaugurated today has five floors and houses 22 Bio safety level(BSL) II laboratories. Setting up of this state-of-art laboratory will give impetus to laboratory diagnosis of various emerging diseases of public health importance. Thanking the Prime Minister Shri Narendra Modi for his visionary leadership and constant guidance, Dr Harsh Vardhan stated that under the charismatic Prime Minister, India is going to witness the elimination of diseases like Tuberculosis, Leprosy, Measles, Malaria and Kala Azar where NCDC is poised to play a very pivotal role. Dr Harsh Vardhan further said that under the leadership of the Prime Minister Shri Narendra Modiji, the Ministry is committed to end TB by 2025, five years ahead of the global target of 2030.
Congratulating NCDC for its contribution in containing Nipah, Dr Harsh Vardhan said that the recent multi-sectoral response co-ordinated by the government for combatting Nipah outbreak in Kerala is another feather in the cap. "It is quite impressive that during the recent reporting of Nipah in Kerala, NCDC has been pivotal in coordination of the response and communication among various sectors and State Govt. It has strengthened the outbreak detection and response capacities and reduced turnaround time for quick decision making. We will take measures to further make it a state-of-art Emergency Operations Centre", the Union Health Minister added. He also commended NCDC's role in supporting the state of Bihar in managing AES.
The NCDC's upgradation plan is also an ambitious plan to keep pace with the changing public health needs of the country and there is a blueprint to build 50 labs in this campus alone, Dr Harsh Vardhan mentioned. He also stated that that six centres of NCDC will soon be set up across the country in Lucknow, Ranchi, Patna, Imphal, Dimapur and Thiruvananthapuram.
Shri Ashwini Kumar Choubey, MoS (HFW) stated that NCDC provides a gamut of services to the people through its laboratories, and strengths in epidemiology, public health capacity building, entomology, etc. NCDC's role in steering the Antimicrobial Resistance (AMR) ContainmentProgramme in the country is commendable. "NCDC can act as a focal point with greater authority and resources for disease surveillance, monitoring of health status, educating the public, providing evidence for public health action and enforcing public health regulations," Shri Choubey mentioned. He also stressed on the importance of public awareness and people's participation in keeping today's lifestyle diseases at bay.
As part of the Annual Day celebrations, a scientific session on "Climate Change: Resilient Health Sector" was also organized by NCDC which was attended by academicians, NCDC and NVBDCP faculty and other eminent experts from the field of public health.
Sh Sanjeeva Kumar, Special Secretary (HFW), Dr AK Saxena, Director General Health Services,Dr Henk Bekedam, WHO Representative to India,Sh Lav Aggarwal, JS(HFW) and Smt Padmaja Singh, JS(VBD) were also present along with other senior officers from NCDC, NVBDCP and the Ministry during the inaugural function.


text:
GDQ is back for COVID-19 relief – Halo, Zelda, and Dark Souls kick things off today
While Summer Games Done Quick has been delayed due to the ongoing pandemic, speedrunners are still getting together under the GDQ banner – this time, to support relief efforts for COVID-19. Corona Relief Done Quick is now underway, kicking off a weekend of speedruns benefiting medical charity Direct Relief.
The marathon is already underway as this story goes live, and is scheduled to last through Sunday, April 19. The day one highlights include 101% Donkey Kong Country kick-off and a co-op Zelda: Ocarina of Time randomizer, but there's plenty of stuff for fans of old PC games today, too – including a Halo 1 run on Legendary, an all bosses Dark Souls effort, and an any% run of Star Wars: Knights of the Old Republic.
You can watch the marathon on Twitch, and donate via the official site. 100% of the donations go straight to Direct Relief. The charity works with health providers worldwide to help distribute resources and provide help as needed – and, of course, the COVID-19 pandemic has strained those resources even further.
Check out the full schedule on the official site, and watch the stream below.
Stick around to enjoy some games and help a good cause.
I'm saving £400 a month but my projected pension at 65 is still puny
I am a bank employee. I work in electrical engineering so don't have any banking or pension knowledge. I have been contributing to the company pension
That amazing Minecraft Harry Potter RPG map has launched
Do you remember that enormous Harry Potter Minecraft map mod that got a playable alpha release earlier this year? The one that recreates the Hogwarts school of Witchcraft and Wizardry as
Google announces good Stadia news, nobody notices because communication has been awful
Stadia is now free. No, it wasn't already free – er, mostly, aside from a free trial, and… well, honestly, I'm not here to tell you
Triumph Studios announces Age of Wonders: Planetfall | PC Invasion
More Age of Wonders is coming from Triumph Studios but Age of Wonders: Planetfall is not quite what you might have expected. Planetfall is a sci-fi themed
Terraria's fourth and final free expansion launches in May
Terraria has moved far beyond its initial appearance as a '2D Minecraft clone', and has remained one of the biggest sandbox games out there for years
The best co-op games on PC in 2020
Co-op games can reward you like no other genre. Games offer up magic in lots of different guises, but working with a friend to achieve a
How to play Magic the Gathering Arena: getting started in MTG
If you're starting out in Magic the Gathering, taking in all the rules of how to play MTG can be a little overwhelming. Even if you've
Monster Hunter World: Iceborne roadmap suffers an indefinite delay due to COVID-19
Monster Hunter World updates are finally in sync across PC and consoles – but that bit of good news comes just as the game's future additions
Free Steam Keys – plus loads more
Everybody loves getting things for free, and us PCGamesN lot are the last people to fly in the face of that tradition. Which is why we
Google has launched its latest Stadia Connect presentation, breaking down what's coming to the platform in the near future and reminding you that, yes, you can