Skip to main content
Advertisement

Main menu

  • Home
  • Articles
    • Newest Articles
    • Current Issue
    • Archive
    • Subjects
  • Submit
    • Submit a Manuscript
    • Author Guidelines
    • License, Copyright, Fee
    • FAQ
    • Why Submit
  • About
    • About Us
    • Editors & Staff
    • Board Members
    • Licensing and Reuse
    • Reviewer Guidelines
    • Privacy Policy
    • Advertise
    • Contact Us
    • LSA LLC
  • Alerts
  • Other Publications
    • EMBO Press
    • The EMBO Journal
    • EMBO reports
    • EMBO Molecular Medicine
    • Molecular Systems Biology
    • Rockefeller University Press
    • Journal of Cell Biology
    • Journal of Experimental Medicine
    • Journal of General Physiology
    • Cold Spring Harbor Laboratory Press
    • Genes & Development
    • Genome Research

User menu

  • My alerts

Search

  • Advanced search
Life Science Alliance
  • Other Publications
    • EMBO Press
    • The EMBO Journal
    • EMBO reports
    • EMBO Molecular Medicine
    • Molecular Systems Biology
    • Rockefeller University Press
    • Journal of Cell Biology
    • Journal of Experimental Medicine
    • Journal of General Physiology
    • Cold Spring Harbor Laboratory Press
    • Genes & Development
    • Genome Research
  • My alerts
Life Science Alliance

Advanced Search

  • Home
  • Articles
    • Newest Articles
    • Current Issue
    • Archive
    • Subjects
  • Submit
    • Submit a Manuscript
    • Author Guidelines
    • License, Copyright, Fee
    • FAQ
    • Why Submit
  • About
    • About Us
    • Editors & Staff
    • Board Members
    • Licensing and Reuse
    • Reviewer Guidelines
    • Privacy Policy
    • Advertise
    • Contact Us
    • LSA LLC
  • Alerts
  • Follow lsa Template on Twitter
Resource
Source Data
Transparent Process
Open Access

A carnivorous plant genetic map: pitcher/insect-capture QTL on a genetic linkage map of Sarracenia

View ORCID ProfileRussell L Malmberg  Correspondence email, Willie L Rogers, View ORCID ProfileMagdy S Alabady
Russell L Malmberg
1Department of Plant Biology, Miller Plant Sciences Building, University of Georgia, Athens, GA, USA
2Institute of Bioinformatics, Davison Life Sciences Building, University of Georgia, Athens, GA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Russell L Malmberg
  • For correspondence: malmberg@uga.edu
Willie L Rogers
1Department of Plant Biology, Miller Plant Sciences Building, University of Georgia, Athens, GA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Magdy S Alabady
1Department of Plant Biology, Miller Plant Sciences Building, University of Georgia, Athens, GA, USA
3Georgia Genomics and Bioinformatics Core, University of Georgia, Athens, GA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Magdy S Alabady
Published 29 November 2018. DOI: 10.26508/lsa.201800146
  • Article
  • Figures & Data
  • Info
  • Metrics
  • Reviewer Comments
  • PDF
Loading

Abstract

The study of carnivorous plants can afford insight into their unique evolutionary adaptations and their interactions with prokaryotic and eukaryotic species. For Sarracenia (pitcher plants), we identified 64 quantitative trait loci (QTL) for insect-capture traits of the pitchers, providing the genetic basis for differences between the pitfall and lobster-trap strategies of insect capture. The linkage map developed here is based upon the F2 of a cross between Sarracenia rosea and Sarracenia psittacina; we mapped 437 single nucleotide polymorphism and simple sequence repeat markers. We measured pitcher traits which differ between S. rosea and S. psittacina, mapping 64 QTL for 17 pitcher traits; there are hot-spot locations where multiple QTL map near each other. There are epistatic interactions in many cases where there are multiple loci for a trait. The QTL map uncovered the genetic basis for the differences between pitfall- and lobster-traps, and the changes that occurred during the divergence of these species. The longevity and clonability of Sarracenia plants make the F2 mapping population a resource for mapping more traits and for phenotype-to-genotype studies.

Introduction

Insect-eating plants fascinate scientists and the general public; Darwin was so taken with their unique adaptations that he used the word “wonderful” 15 times in his descriptions of them (Darwin, 1888). The carnivorous behavior of plants evolved independently, possibly as many as nine times (Givnish, 2015; Wheeler & Carstens, 2018). These organisms live in nutrient poor conditions, such as wetland bogs, and are able to obtain minerals, primarily nitrogen and phosphorous, from capturing and digesting insect prey (Darwin, 1888; Ellison & Adamec, 2018; Adlassnig et al, 2012; Bradshaw & Creelman, 1984). Their leaves are specialized to perform multiple functions; secrete attractive scents (Jurgens et al, 2009), capture insects, secrete extracellular digestive enzymes, absorb nutrients, photosynthesize, and develop microbial symbioses. Comparative genomic approaches are beginning to be applied to identify sequences associated with the evolution of carnivory (Wheeler & Carstens, 2018). However, until this report, a classic genetic linkage mapping approach has not been implemented with these systems; we have used pitcher plants of the genus Sarracenia to develop such a genetic map.

Sarracenia pitcher plants are potentially useful for addressing a number of developmental, physiological, ecological, and evolutionary questions. Some of these are: How the pitchers actually digest the insects and obtain nutrition is of interest; for example, the flow of nitrogen and phosphorous (Gallie & Chang, 1997; Ellison, 2006; Butler & Ellison, 2007; Karagatzides et al, 2009, 2012) has been studied. Pitcher plants provide an opportunity to study their unique interactions with microbial communities; they contain a microbiome associated with insect digestion within their pitchers (Buckley et al, 2003; Gotelli & Ellison, 2006; Koopman et al, 2010; Koopman & Carstens, 2011; Baiser et al, 2012), a partially enclosed container offering an experimentally manipulable system suitable for studying eukaryotic-host/microbiome relationships. The plants and their pitchers can be used to study the developmental, evolutionary, and ecological tradeoffs associated with phenotypic plasticity as under high nitrogen conditions, the leaves will grow as flattened phyllodia (a leafstalk flattened into a blade shape which is photosynthetically more efficient than a pitcher [Ellison & Gotelli, 2002]) instead of as a pitcher. There are 14–50 named taxa within the genus (Schnell & McPherson, 2011; Weakley, 2015), and these can all hybridize with each other, yielding abundant material for investigating the role of hybridization in plant evolution. Some of these taxa are endangered, posing conservation genetics issues, such as how much genetic diversity is present in given populations, and what is the best way to preserve it? All of these example areas of inquiry would benefit from the availability of a genetic linkage map.

Pitcher plants are also charismatic attractive plants which can captivate broad audiences and help raise awareness of scientific issues in the general public. For an overview—a lively recent book covers many aspects of Sarracenia and carnivorous plant biology in detail (Ellison & Adamec, 2018).

Some of the difficulties in using Sarracenia as a genetic mapping system are: (a) it takes 3–4 y for a plant to produce seeds, (b) 1 cM on a genetic map would correspond on average to more than 106 basepairs because of the large genome size (Rogers et al, 2010), and (c) the presence of a recent partial genome duplication (Srivastava et al, 2011) will complicate identifying orthologs. On the other hand: (a) the species/taxa can all be crossed with each other, and natural hybrids have been found for most of the combinations; (b) pitcher plants are long-lived perennials and can be readily vegetatively propagated, so that after genotyping individuals, a mapping population can be reused multiple times and the same genotype grown under multiple growth conditions for phenotyping.

We report developing an F2 generation of more than 200 individuals from a cross between Sarracenia rosea (Sarracenia purpurea venosa burkii) and Sarracenia psittacina, and use this F2 to create a genetic map with markers and pitcher quantitative trait loci (QTL). S. rosea was chosen as a parent as it is the Southeastern U.S. variant of S. purpurea, the most common, most familiar, and widely distributed of the Sarracenias in North America. S. psittacina was chosen as it is the most morphologically diverse and distinctive of the Sarracenias, with a decumbent, horizontal, pitcher growth habit; its pitchers have been described as lobster-trap–like with respect to catching insects, meaning that the insects can enter readily but have difficulty exiting. The other Sarracenias can be described as having pitfall traps into ponds. S. rosea and S. psittacina are in different subgroups of the genus, based upon a recently developed phylogeny of the taxa within the genus (Stephens et al, 2015) generated from the sequences of 199 nuclear genes. S. rosea groups with the other S. purpurea-related taxa, whereas S. psittacina groups with Sarracenia minor–related and Sarracenia flava–related taxa; hence, it is reasonable to expect some genetic differentiation between S. rosea and S. psittacina suitable for creating a genetic map. A third sub-grouping within the genus contains Sarracenia oreophila, Sarracenia rubra–related taxa, Sarracenia alata, and Sarracenia leucophylla.

Some of the pitcher traits differing between S. rosea and S. psittacina have previously been investigated to identify their roles in insect attraction and capture (Horner et al, 2018). For example, the striking patterns of red stripes on pitchers, and other carnivorous plants, have been suggested to be involved in insect attraction (Schaefer & Ruxton, 2008). Many Sarracenia also have white windows or fenestrations in patterns on the pitchers which have been hypothesized to confuse trapped insects (Moran et al, 2012). There are morphological differences between S. psittacina pitchers and the pitchers of other Sarracenia, including the size and orientation of the orifice, the shapes of the hoods above and around the orifice, the presence of tissue and hairs apparently meant to trap the insects within the pitchers, and the overall small size of the S. psittacina pitcher (Naczi, 2018); Naczi made the intriguing observation and suggestion that S. psittacina pitchers are more similar to juvenile pitchers of other Sarracenia, and hence that their developmental evolution may be characterized as heterochrony. These pitcher differences between S. rosea and S. psittacina are thus some of the possible candidate traits for QTL genetic mapping.

We used the methods of genotyping by sequencing and restriction site–associated RNA-sequencing (RARseq [Alabady et al, 2015]) to generate linkage groups and a genetic map with more than 430 mapped markers covering 2017 cM. We identified and mapped 64 QTL for 17 pitcher traits which differ between the two parental species and which give insights into the genetic bases of the pitfall versus lobster-trap insect-capture strategies of Sarracenia.

Results

Linkage map

Hecht (1949) has described the chromosomes in root tip metaphases of six Sarracenia species, including S. purpurea and S. psittacina, as 2N = 26 or 1N = 13; a finished genetic map might thus have 13 linkage groups. Figs 1 and 2 show the genetic linkage map we developed, with linkage groups shown that are greater than 10 cM, or that are shorter than 10 cM but which also have a QTL mapping to them. Most of the markers we used for genetic mapping were single nucleotide polymorphisms (SNPs) identified by sequencing from RNA (RARseq) or DNA (RADseq), plus there were a few simple sequence repeat markers (SSRs). The final total of 437 is made up of 343 RNA-based SNPs, 78 DNA-based SNPs, and 16 SSRs. These are distributed into 42 linkage groups with a total length of 2,017 cM; 12 of the linkage groups contain a small number of markers and are less than 10 cM in length. The overall average density is 1 marker every 4.6 cM. The 13 largest linkage groups contain 273 markers with a total length of 1,446 cM, and with the shortest of the 13 having a length of 76.7 cM. The DNA-based markers were generally interspersed among the RNA-based markers, with 39 of them either a single isolated DNA marker or a doublet of two DNA markers. There were a few clusters of DNA-based markers: linkage group 24 was comprised entirely of 6 RADseq markers covering 20 cM whereas linkage group 27 was entirely 3 RADseq markers over 13 cM; the end of linkage group 1 contained 5 RADseq markers over 37 cM; near one end of linkage group 8 there was a group of 5 RADseq markers and 1 SSR over 28 cM.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1. Linkage groups with QTLs mapped (continued in Fig 2).

Linkage groups are shown and numbered in decreasing order of length in cM. Linkage groups 1–30 are shown as they are greater than 10 cM in length; in addition groups 31, 34, and 39 are shown as QTLs mapped to them. QTLs are indicated by their name and by vertical lines which indicate a Bayesian position confidence interval, as calculated by R/qtl (Broman et al, 2003; Broman & Sen, 2009; Arends et al, 2010). When there are multiple loci for the same quantitative trait, these are indicated by a numbered trait name, which are numbered in decreasing order of %PVE such that number one indicates the locus with the highest %PVE. The genetic map was rendered by the R/LinkageMapView package (Ouellette et al, 2018).

Source data are available for this figure.

Source Data for Figure 1[LSA-2018-00146_SdataF1.zip]

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2. Linkage groups with QTLs mapped (continued from Fig 1).

Linkage groups are shown and numbered in decreasing order of length in cM. Linkage groups 1–30 are shown as they are greater than 10 cM in length; in addition groups 31, 34, and 39 are shown as QTLs mapped to them. QTLs are indicated by their name and by vertical lines which indicate a Bayesian position confidence interval, as calculated by R/qtl (Broman et al, 2003; Broman & Sen, 2009; Arends et al, 2010). When there are multiple loci for the same quantitative trait, these are indicated by a numbered trait name, which are numbered in decreasing order of %PVE such that number one indicates the locus with the highest %PVE. The genetic map was rendered by the R/LinkageMapView package (Ouellette et al, 2018).

Source data are available for this figure.

Source Data for Figure 2[LSA-2018-00146_SdataF2.zip]

The STACKS program (Catchen et al, 2013) generated sequence tags of 143 bases for each of the SNPs reported, as listed in Supplemental Information 5. We performed Blast2Go (Conesa et al, 2005; Gotz et al, 2008) analyses of these sequence tags for the SNPs (Supplemental Information 7). We could not detect a particular pattern in the functional annotations that would suggest anything about the nature of the sequences containing SNPs either overall or just the SNPs located near the mapped QTLs.

We did notice five instances of SNP tag sequence similarity where the SNPs were located in different regions of the genetic map, suggesting gene duplications. These included:

  • R7834 (lg3-pos32) and RI736 (lg3-pos59) with 100% identity over 143/143 bases;

  • R6009 (lg3-pos139) and R6000 (lg16-pos27) with 90% identity over 72/143 bases;

  • R8080 (lg4-pos40) and R8078 (lg4-pos70) with 100% identity over 120/143 bases;

  • R7818 (lg7-pos99) and R7820 (lg22-pos20) with 90% identity over 128/143 bases;

  • R1768 (lg8-pos0) and R1769 (lg15-pos0) with 100% identity over 97/143 bases.

Pitcher traits mapped in F2

We measured a variety of pitcher traits that are different between S. rosea and S. psittacina to attempt to map their genetic bases. Fig 3 shows the parental plants and F1, giving an overall impression of the plants. For each quantitative trait mapped, the R/qtl fitqtl() function gives an estimate of the total phenotypic variance which is explained (%PVE, the heritability of the trait) in the F2 population by the postulated loci and the pairwise interactions, epistasis, between them. We describe here the traits mapped and the genetic loci and interactions involved, giving them in descending order of percent phenotypic variance explained; the eight letter codes shown are those used in the data files and in Table 1. When more than one locus is involved, the loci are numbered with locus-1 having the highest individual %PVE, and then descending in order; thus for traitx, traitx-1 has the highest individual %PVE, traitx-2 has the second highest, and so on. Figs 4, 5, 6, 7, 8, and 9 illustrate each of the individual pitcher traits measured. Table 2 summarizes the results for each trait and locus, whereas Table 3 lists the epistatic interactions detected.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3. Pitchers of the parents and the F1.

(A) S. purpurea. (B) F1. (C) S. psittacina.

Source data are available for this figure.

Source Data for Figure 3[LSA-2018-00146_SdataF3.zip]

View this table:
  • View inline
  • View popup
Table 1.

The 17 quantitative traits where a locus or loci mapped, listed in the order discussed in the text.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4. Examples of traits scored.

Listed in order of decreasing %PVE, the same order as discussed in the text. (A–C) Dots indicate pitchers in rosette to be counted. (D, E) The solid arrow points to the glossy zone free of hairs. (F) The dashed arrow points to hairs throughout. (G) The solid arrow points to the smooth zone in pitcher otherwise more similar to S. psittacina. (H) The arrows point to absent or weak windowing effect. (I, J) The arrow points to white window effect. (K) Well-defined white windows due to distinct color borders.

Source data are available for this figure.

Source Data for Figure 4[LSA-2018-00146_SdataF4.zip]

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5. Examples of traits scored.

Listed in order of decreasing %PVE, the same order as discussed in the text. (A) The lip is curled outward. (B) The lip extends slightly upward. (C) The lip extends upward and inward. (D) There is no visible tissue above lip. (E) There is a small amount of tissue upward from lip. (F) There is tissue extending upward and inward from lip.

Source data are available for this figure.

Source Data for Figure 5[LSA-2018-00146_SdataF5.zip]

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6. Examples of traits scored.

Listed in order of decreasing %PVE, the same order as discussed in the text. (A–C) Do the veins and windows have similar or different color patterns? (D–F) The degree of openness of pitcher. (G–I) The length of the region containing hairs from base upward.

Source data are available for this figure.

Source Data for Figure 6[LSA-2018-00146_SdataF6.zip]

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7. Examples of traits scored.

Listed in order of decreasing %PVE, the same order as discussed in the text. (A–C) How well defined and easily traceable are the veins? (D–F) The length of tissue above the lip line. (G–I) The width of the leaf directly below the mouth opening.

Source data are available for this figure.

Source Data for Figure 7[LSA-2018-00146_SdataF7.zip]

Figure 8.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 8. Examples of traits scored.

Listed in order of decreasing %PVE, the same order as discussed in the text. (A–D) The fraction of pitchers open. (E–G) Green versus red/purple leaf color. (H–K) The fraction of pitchers with green versus red/purple leaves.

Source data are available for this figure.

Source Data for Figure 8[LSA-2018-00146_SdataF8.zip]

Figure 9.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 9. Examples of traits scored.

Listed in order of decreasing %PVE, the same order as discussed in the text. (A–C) The width of the leaf at its widest spot. (D–F) The angle between a flat area on the back of the pitcher and the tip of the operculum. (G–J) The degree of radial symmetry of rosette, with asymmetry seen in F2.

Source data are available for this figure.

Source Data for Figure 9[LSA-2018-00146_SdataF9.zip]

View this table:
  • View inline
  • View popup
Table 2.

Quantitative traits mapped and their loci.

View this table:
  • View inline
  • View popup
Table 3.

Interactions between loci for the same trait with >1% phenotypic variance explained.

There were nine traits measured which failed to show any heritability; these are described in the Materials and Methods section for their value as negative data.

  • numpitch: Number of pitchers/leaves per rosette; 100% PVE; 9 loci plus 36 of 36 pairwise interactions (Fig 4A–C).

This is a count of the leaves/pitchers on a single rosette. S. rosea usually totals 4–6 leaves whereas S. psittacina can have up to 30. There are nine loci, and all pairwise interactions are significant. The pairwise interactions were individually small in their %PVE but totaled to 39% for the 36 pairs. The dominance relationships were: no dominance—loci 4, 6, 7, and 9; S. rosea dominance—locus 8; S. psittacina dominance—loci 1, 2, 3, and 5. The allele effect values were mixed across the multiple loci, with some in the direction of S. rosea and some in the direction of S. psittacina.

  • smoothzn: Smooth area lacking hairs (trichomes) within pitcher; 97% PVE; 9 loci plus 36 of 36 pairwise interactions (Fig 4D–G).

In S. rosea there is a glossy smooth zone with a polished look directly above the hairy zone at the base of the pitcher. S. psittacina lacks this smooth zone entirely; sometimes at the lip there can be a very small area of waxy surface but never this smooth zone. The F2s may have sporadic or plentiful hairs in a smooth zone, or pores for the hairs but no actual hairs seen. Nine significant loci were identified, all with similar %PVE. All the pairwise interactions were significant; however, the sizes of the %PVE for each of the interactions were small. For the 36 interactions, the sum of their %PVE was 32%. The dominance patterns were: no dominance—loci 2, 4, and 7; over/under-dominance—loci 1, 6, and 9; S. rosea dominance—locus 8; S. psittacina dominance—loci 3 and 5. The allele effect values were mixed across the multiple loci, with some in the directions of each parent, but more in the direction of S. rosea.

  • fenestra: White windows on the pitcher; 58%–99% PVE; 15 loci plus 65/105 interactions (Fig 4H–K).

S. psittacina pitchers have a windowing-effect pattern of white areas, whereas S. rosea pitchers do not. This was scored on a 6-point scale with 0 being S. psittacina–like, and 6 being S. rosea–like. A total of 13 significant loci were identified with 2 additional loci included as having significant interactions with other loci. The number of possible pairwise interactions among the 15 loci was too large to describe in a single genetic model to be evaluated by fitqtl(); the 55% PVE listed for the full model, and the values shown for the individual loci were obtained from a purely additive model. We explored the pairwise interactions by testing subsets of them in groups together with the additive effects; in some cases the final %PVE was as high as 99%; 65/105 interactions were significant. In the purely additive model, the locus fenestra-1 explained 15% of the phenotypic variance, whereas other loci explained 7% or less. The dominance relationships were: no dominance—loci 1, 2, 3, 4, 7, 10, and 15; over/under-dominance—loci 8, 11, and 12; S. rosea dominance—loci 5 and 14; S. psittacina dominance—loci 6, 9, and 13. The allele effect values were mixed across the loci, with some in the directions of each parent.

  • periotuk: Periostome in or out; 50% PVE; 5 loci plus 1 of 10 pairwise interactions (Fig 5A–F).

Is the periostome tucked inward, S. psittacina–like as part of the lobster-trap morphology, or protruded, S. rosea–like as part of the pitfall morphology? This was scored on a 4 point scale with S. psittacina–like inward tucking as 0, S. rosea–like protrusion as 3, and with intermediates scored as 1 or 2. There were five significant loci, with one of the possible pairs, periotuk-2 and periotuk-3, having a significant interaction. Periotuk-1 and periotuk-5 have dominance by the S. psittacina allele; periotuk-2 and periotuk-3 show no dominance; periotuk-4 has dominance by the S. rosea allele. The allele effect values were generally in the direction of S. rosea.

  • heigwing: Height of wing; 43% PVE; five loci (Fig 5G–I).

This is the height of the wing perpendicular to the pitcher from its highest point directly back to where it meets and fuses with the leaf/body of the leaf. Heigwing-1 and heigwing-3 had no dominance; heigwing-2 and heigwing-5 had over- or underdominance; heigwing-4 had dominance by the S. rosea allele. The allele effect values were mixed in the directions of both parents.

  • ptnwindw: Pattern of Color; 30% PVE; three loci plus two of three pairwise interactions (Fig 6A–C).

This is a measure of the patterning of the color with respect to veins and windows (fenestrations) on the pitcher. The color can be either green or a shade of purple. It was scored on a 0 to 3 point scale, where 0 is a solid color with veins not easily distinguishable (S. rosea), three indicates the veins and windows are highly noticeable (S. psittacina), and one and two are intermediate states. Ptnwindow-1 and ptnwindow-3 have dominant S. psittacina alleles, whereas ptnwindow-2 shows no dominance. Ptnwindow-1 has significant interactions with each of ptnwindow-2 and ptwindow-3. The allele effect values were in the direction of S. psittacina.

  • openness: Degree of openness of pitchers; 26% PVE; three loci (Fig 6D–F).

This measures the degree of openness of the operculum to the pitcher mouth, and is related to, but different, from openfrac which is the fraction of pitchers open. In S. rosea these are always open (pitfall trap), whereas in S. psittacina they are nearly always sealed closed except for a small opening (lobster-trap). The degree of openness was measured on a six point scale with five being S. rosea–like and zero being S. psittacina–like. There were three significant loci. Openness-1 and openness-2 both have dominance by the S. psittacina allele, whereas openness-3 shows no dominance. The allele effect values were in the directions of S. psittacina.

  • lengintr: Internal length hairy region; 23% PVE; two loci plus zero of one pairwise interaction (Fig 6G–I).

This is the distance from the base of the leaf where it attaches to the meristem upward until where the line of hairs ends, measured in pitchers which have been sliced open. In S. rosea this is at the very base of the leaf and does not extend very far upward; there is usually a pronounced clearly visible line at the base signaling the end of the hairy section. In S. psittacina this includes the entire length of the leaf right up to the peristome. The F2s have a wide range of region lengths and can have sporadic hairs appearing out of place. There are two loci, with the lengintr-1 showing no dominance and lengintr-2 showing partial dominance by the S. psittacina direction.

  • veinprom: Well-defined veins; 18% PVE; one locus (Fig 7A–C).

S. psittacina tends to have well-defined easily traceable veins, whereas S. rosea can have veins which are less clearly visible in the context of the leaf. This was scored on a 3 point scale with 0 being most well-defined, and 2 being least defined. There was one significant locus with no dominance.

  • Abovelip: Tissue above the lip line; 17% PVE; two loci plus one of one pairwise interaction (Fig 7D–F).

This measures the length of the tissue above the lip line. When one slices the leaves open (down the back of the leaf and filet the leaf open), S. psittacina has tissue extending upward protruding into the bell that makes the lobster-trap. When one does the same thing to S. rosea there is no tissue extending upward but instead tissue curled down and under to form the lip rim as part of the pitfall trap. There are two significant loci with a significant interaction between them. Considered individually, abovelip-1 has dominance by the S. psittacina allele, whereas abovelip-2 has dominance by the S. purpurea allele; the highest value occurs in individuals heterozygous for both loci.

  • widtleaf: Leaf width, 17% PVE, two loci plus one of one pairwise interaction (Fig 7G–I).

This is the width of the leaf directly below the mouth opening. There are two significant loci, and the interaction between the two is also significant. The locus effects were in opposing directions with the largest widths in pitchers homozygous for the widtleaf-1 locus from S. psittacina and widtleaf-2 locus from S. rosea. There is no significant dominance.

  • openfrac: Fraction of pitchers open (pitfall) versus mostly closed (lobster); 16% PVE; two loci plus zero of one pairwise interaction (Fig 8A–D).

This is the fraction of pitchers which are open, S. rosea pitfall-like, versus mostly closed with just a small periostome opening at the base, S. psittacina or lobster-trap–like. There are two loci, with no interactions; openfrac-2 shows dominance by the S. psittacina allele.

  • leafcolr: Leaf color, green versus purple, ignoring veins; 16% PVE; two loci plus one of one pairwise interaction (Fig 8E–G).

This is the degree of purpleness, with green scored as 0, intense purple scored as 3, and intermediate shades as 1 or 2. Two significant genetic loci were found, with an interaction between them. The most intense color was found when the two loci were both homozygous for the S. psittacina allele, but the least intense was when leafcolr-1 was homozygous for the S. psittacina allele and leafcolr-2 was homozygous for the S. rosea allele. Leafcolr-1 had some dominance by S. psittacina, and leafcolr-2 showed no dominance.

  • deggreen: Degree of green color; 8% PVE; one locus (Fig 8H–K).

This is an overall estimate of how much green versus red there is in an individual. In some F2 plants, individual pitchers could vary significantly from each other. Plants were scored on a 0 to 3 point scale where 0 is all green and 3 is no green. There was one locus found, explaining 8% of the phenotypic variation, with dominance by the S. rosea allele.

  • widtspot: Widest leaf width; 8% PVE; one locus (Fig 9A–C).

This is the width of the leaf at its widest spot, which may be either above or below the mouth, whereas the widtleaf measure is always just below the mouth opening, and maps to different loci. The one locus found shows dominance by the S. psittacina allele.

  • curvatur: Pitcher curvature, 6% PVE; one locus (Fig 9D–F).

This was measured as the angle between a flat area on the back of the pitcher and the tip of the operculum. For S. rosea, this number is usually 90° but for S. psittacina it is always less than 50°. One significant locus was found which explained 6% of the phenotypic variation. There is no evidence for dominance.

  • radalsym: Rosette radial symmetry, 4% PVE, one locus (Fig 9G–J).

This is a visual observation of whether or not the rosettes have radial symmetrical looking downward at them. In some non-symmetric F2 plants, the pitchers bunch together more on one side of the rosette than the other; or the pitchers may project upward and outward at irregular angles. This was scored as a binary trait, 1 or 0. A single significant genetic locus was found which explained 4% of the phenotypic variation. S. rosea has greater symmetry than S. psittacina; and this trait had significant dominance by the S. psittacina direction.

Traits in the F1

The F1 plants were largely intermediate in all traits examined between the two parental species (Fig 3). Many of the traits we measured had multiple loci detected, with a mix of dominance relationships and allele effects; hence, the intermediate state of the F1 is not surprising.

Discussion

We have constructed a linkage map of Sarracenia using an F2 generation of a cross between two species, S. psittacina and a subspecies of S. purpurea, S. rosea. For 17 pitcher traits which differ between these taxa, we identified 64 loci which were placed on the genetic map. This is the first genetic linkage map and QTL mapping for a carnivorous plant.

A genetic linkage map for Sarracenia

Our linkage map contains 437 markers across 42 linkage groups, with a total length of 2,017 cM. Twelve of these linkage groups were 10 cM or less in length. Previous cytogenetic work (Hecht, 1949) indicated that there are 1N = 13 chromosomes. The 13 longest of the linkage groups we have identified range from 189.9 cM to 76.6 cM with a total length of 1,446 cM. For comparison, grape, Vitis, is a species with extensive sequencing, and its transcript sequences are some of the most similar to those of Sarracenia transcripts as detected by Blast comparisons (Srivastava et al, 2011). A recent Vitis linkage map (Wang et al, 2017) had 6,000 markers with a total length of 2,186 cM in 19 linkage groups with lengths ranging from 86 cM to 237 cM. A linkage map of maize (Davis et al, 1999), certainly one of the genetically most well-studied plants, contains 10 linkage groups representing 10 chromosomes, and has a total length of 1,723 cM; the maize genome is about 70% of the size of the Sarracenia genome. These comparisons suggest that our current linkage map covers most of the total genetic map, more than 80% of it, and lead to a rough global estimate of 1.4 × 106 bp per cM for Sarracenia.

There were several difficulties in constructing the genetic map: it took a long time, more than 10 y, to generate the F2 plants. The parents were heterozygous at multiple loci as we did not start with pure lines as parents (Mendel, 1951). Sarracenia has a relatively large genome size estimated to be 3.6 × 109 basepairs (Rogers et al, 2010). There is a recent partial genome duplication which dates to about the time of the divergence of the species within the genus (Srivastava et al, 2011), perhaps 2 million years ago. These difficulties likely lead to the map being incomplete in spite of our multiple rounds of sequencing efforts. Although we were successful in generating some SNP markers which were able to be mapped with the RADseq protocol (Elshire et al, 2011) (78 SNPs), the RARseq (Alabady et al, 2015) method of generating markers from RNA was more productive (343 SNPs). A few of the RNA-based markers we mapped showed evidence of duplications. The existence of high quality reference genome sequences from both the two parents, S. rosea and S. psittacina would likely greatly improve the genetic map.

Pitcher trait inheritance

The focus of our trait measurements has been on pitcher traits as insectivory is a hallmark of the genus, and insectivorous adaptations are the reason for the general fascination with these plants. There are a wide diversity of pitcher morphologies within the Sarracenia genus, which may be associated with their insect-eating lifestyle, and may represent differences in insect-catching and digestion strategies across the species. We picked 26 different pitcher traits to measure; for 17 of these traits we were able to map significant QTLs, whereas nine of them failed to map. Our current study established the genetic system and was an initial exploration of quantitative trait heritability in one particular setting. Any of these traits may behave differently when the plants are grown in a different environment. Some of the traits which showed no inheritance or low %PVE could have higher heritability if studied somewhat differently; conversely a different environment might reduce the %PVE for other traits.

The %PVE varied from 3.7% (radalsym) to 96.5% (smoothzn) and 99.6% (numpitch) for the full models. Radalsym was presence or absence of a lack of radial symmetry in the pitcher rosettes (Fig 9G–J); both parents have greater radial symmetry than some of the F2 plants which showed the most extreme phenotype. Smoothzn is the glossy internal region free of hairs found in S. rosea (Fig 4D–G), whereas numpitch is the number of leaves/pitchers per rosette (Fig 3A–C). The next highest %PVE after numpitch and smoothzn was 55.3% (periotuk) (Fig 5A–F). The high %PVE may partially be due to these two traits being straightforward to score in a quantitative fashion. Other traits with lower %PVE may show more inheritance if scored differently or if the plants were grown differently. For single loci, the %PVE fell between 2 and 21%. The highest value was 20.4% for ptnwindw-1, one locus for a trait (Fig 6A–C) which measured the pattern of color in veins and non-vein areas. Some of the traits with multiple loci identified had a single locus which explained most of the variation together with additional loci of smaller effect; other traits had several loci with similarly sized %PVE. There was not a discernible pattern to the dominance relationships: semidominance, over- or under-dominance, and dominance from both parents were found. For the allele effects, when multiple loci were found for some traits, the pattern of effects sometimes pointed to one parent, but sometimes was mixed with direction pointing to both of the two parents.

Table 3 summarizes pairwise interactions among the loci where these were significant and the total %PVE explained by the interaction was greater than 1%. For the traits smoothzn, numpitch, and fenestra, which had large numbers of loci identified, the %PVE for any one interaction was small, on the order of 1%–2%.

The trait fenestra (Fig 4H–K), white windows on the pitchers, had the most loci (15) in the initial analysis. Because of computational limitations, we were not able to test all possible pairwise interactions among these in the same way as we could for the traits with smaller numbers of loci. We therefore tested a number of QTL models which included the additive and random subsets of the interaction possibilities. The values shown in Table 2 for fenestra are from evaluation of an additive model, not allowing any interactions. The bottom of Table 3 lists interactions among the fenestra loci which were potentially significant in the exploratory analyses.

The genetic map, Figs 1 and 2, shows the traits mapped with a Bayesian confidence window, as calculated by R/qtl (Broman et al, 2003; Arends et al, 2010). The position confidence intervals vary from being very large, poorly defined and taking up most of a linkage group, to appearing very small and highly defined. There is an inverse relationship between the %PVE for whole trait, and the smallness of the confidence interval. Thus, radalsym, with the lowest %PVE, has its position poorly defined on a linkage group. In contrast, the numpitch and smoothzn loci appear to almost have point estimates of their position. The same F2 plants were scored for all the traits. A trait with a low %PVE effectively has less of its variation available to indicate the genetic map position than a trait with a high %PVE and hence has a more poorly defined position.

There are some intriguing overlapping localizations of loci from multiple traits appearing in the same map region. Linkage group 3 has 12 loci from smoothzn, fenestra, periotuk, heigwing, ptnwindow, openness, veinprom, abovelip, openfrac, and deggreen mapping in the same general vicinity. Several other such QTL hotspots seem to exist. Assuming that these hotspots are meaningful and not by chance, a variety of mechanisms might potentially be involved in creating and maintaining them: pleiotropy, the influence of a single gene on seemingly unrelated phenotypic traits (Plate, 1910; Stearns, 2010) can create multiple phenotypes from a small number of genes; selective sweeps, when a selectively favored allele increases the frequencies of nearby genetic variants by hitchhiking (Smith & Haigh, 1974), can lead to co-segregation of the linked genes, associating the phenotypes; chromosomal rearrangements, such as inversions, can suppress recombination within a region, also causing association of phenotypes.

Pitfall and lobster traps

The genetic cross we performed for mapping was between S. psittacina and S. rosea. S. rosea is the southeastern variant of S. purpurea, and is a member of the medium-height group of Sarracenia taxa. S. psittacina pitchers have been described as lobster-trap–like with respect to catching insects, whereas the other Sarracenia are described as having pitfall traps.

Some of the loci we have identified are part of the genetic basis for the lobster-trap versus pitfall morphologies. The periotuk trait describes the inward folded tissue of the lobster-trap, and we identified five loci involved in this, whereas abovelip measures the amount of tissue in this folding, and we found one locus for this. These are both lobster-trap traits from the S. psittacina parent. The two open traits openness and openfrac describe how wide the pitcher opening is, with the wider pitcher being an S. rosea pitfall trait, and between which we found a total of four loci. The inside surfaces of pitchers of S. purpurea and S. rosea have a smooth, glossy, hair-free zone above a more hirsute region, whereas S. psittacina has hairs (trichomes) throughout. The presence of hairs is possibly important for S. psittacina to help keep the insects trapped inside, in combination with the narrow opening and inward periostome, whereas S. rosea has its pitcher pond, and thus may not need hairs in that region. Fig 4G is an example of an F2 pitcher which is largely S. psittacina–like in shape, but which has some smoothzn. Fenestration, the pattern of white windows, has been interpreted as being involved in insect prey-attraction (Schaefer & Ruxton, 2014). This is found in S. psittacina, not S. rosea; we found the fenestra quantitative trait has 13–15 loci and potentially a large number of interactions among them. The genetic loci underpinning the periotuk, openness, openfrac, smoothzn, and fenestra traits are thus involved in the pitfall versus lobster-trap differences in insect-capture strategy.

Evolution

The phylogeny of the genus from Stephens et al (2015) divides the taxa into two main groups—tall pitcher plants (S. oreophila, Sarracenia alabamensis, S. alata, S. rosea, and related taxa) and mixed-sized pitcher plants (S. purpurea, S. minor, S. psittacina, S. flava and other taxa); S. psittacina is the shortest, smallest, and most horizontal in its growth habit. The radiation of the taxa within the genus was relatively recent, possibly 2 million years ago at a time coincident with a partial genome duplication (Srivastava et al, 2011). Our analysis detected 64 trait loci which differ between S. rosea and S. psittacina. The analysis of crosses between similar plant species has shown examples of both a small number of major effect loci, and also a larger number of minor effect loci for traits contributing to the evolution or domestication of these species (as examples [Doebley, John, 1992; Wills & Burke, 2007; Brandvain et al, 2014; Fishman et al, 2014; Wang et al, 2015; Garner et al, 2016; Kenney & Sweigart, 2016; Badouin et al, 2017]). The pattern we have seen in the pitcher traits examined so far is more similar to having a large number of minor effect loci. We have not examined floral traits whose inheritance might be expected to provide insight into the divergence of S. rosea and S. psittacina; superficially, the flowers of the two parental species differ modestly in size and color, but their overall morphology is very similar. Particularly where there are multiple loci for the same trait, some of the genetic differences we have found might be related to duplications of portions of the genome, beyond the five pairs of loci listed which map at different locations. The multiple loci may be identifying paralogous genes.

Prospects

We have generated a linkage map which is 80% complete; the successful mapping of a number of pitcher QTLs demonstrates some of the potential of this system. We have studied a variety of pitcher morphological traits, but there are many other questions that could be similarly addressed with these F2 plants by measuring appropriate traits, provided they differ between S. rosea and S. psittacina. As examples: The roles of the various pitcher components in insect capture and digestion might be studied in F2 plants which contain various combinations of the traits, such as hairs (trichomes), fenestrations, and periostome structures, to determine the effect of each; specific combinations of QTL could be used to test results from the presence and absence of particular pitcher features. Our pitcher QTL were measured at only a relatively full grown stage of the pitchers; comparison of pitchers of S. rosea, S. psittacina, and the F2 plants at various developmental stages might help examine the hypothesis (Naczi, 2018) that heterochrony is a part of the S. psittacina developmental evolution. Hybridization can occur between all species or taxa within the genus, raising questions about what is the nature of a species in Sarracenia, and what evolutionary genetic mechanisms exist which might maintain a species in the face of rampant hybridization. The genetic map and the F2 population could be tools which contribute to solving these problems.

One of the attractions of this system is that the F2 plants are long-lived and are easily subdivided for vegetative propagation. Given a set of F2 plants which have been genotyped, it should be possible to readily map multiple additional QTLs by scoring the trait of interest across these plants, then performing a QTL analysis with the same marker and map files. For any traits of interest, these could be scored under multiple environmental conditions by growing clones of the characterized genotypes. Similarly, additional analyses might improve the number of markers in each F2 plant, and then it would be straightforward to reanalyze the data with a more detailed genetic map.

Materials and Methods

Cross

The initial cross was between the southeastern variety of S. purpurea, S. rosea (S. purpurea venosa burkii) (male), and S. psittacina (female). F1 plants were grown and selfed to generate F2 individuals. All plants were grown in the Plant Biology Department greenhouses, University of Georgia, Athens, Georgia, USA (33°55′45.7″N 83°21′49.6″W).

RNA-based markers

We used the restriction site associated RARseq method, as described by Alabady et al (2015), to generate RNA-based SNP markers. RNA was isolated from each of the F2 plants then these were reverse transcribed into cDNA using the methods described (Alabady et al, 2015). The cDNA samples were digested with restriction MseI, then the resulting fragments were size fractionated on an agarose gel to produce fragments of 250–600 bp. These fragments were ligated to oligonucleotide adapters which contained unique barcodes for each individual plant and which prepared the samples as libraries for Illumina sequencing. The sequencing was performed on the NextSeq platform.

We generated reference transcriptome sequences from S. rosea and S. psittacina with deep coverage, using both nextseq (Illumina) and long read Iso-Seq (PacBio) platforms. SNP identification was performed using the program STACKS (Catchen et al, 2013) version 1.48; we performed reference-based SNP detection by comparison of the sequences from the F2 plants to the reference transcriptomes. In the final marker set used for mapping, the SNPs detected from the NextSeq ranscriptome reference begin with the letter R followed by the number assigned by STACKS. The actual SNPs in their sequence contexts are given in the Supplemental Information 5.

DNA-based markers

We used standard restriction site–associated DNA sequencing (RADseq) methods (Elshire et al, 2011) to generate a subset of genomic fragments for sequencing. A double digest with MspI and SbfI succeeded in generating sequences with a reasonable probability of being matched across the F2 individuals. After the size selection, the fragments were ligated to oligonucleotide adapters which contained unique barcodes for each individual plant and which prepared the samples as libraries for Illumina sequencing. The sequencing was carried out on the NextSeq platform. SNPs were detected by the software STACKS as a denovo, non-referenced–based, pipeline. In the final marker set used for mapping, the SNPs detected from RADseq begin with the letter D followed by the number assigned by STACKS. The actual SNPs in their sequence contexts are given in the Supplemental Information 5.

A total of 16 SSRs (microsatellites), isolated as described in Rogers et al (2010) were also used as markers in the map construction. The SSR primer sequences are given in the Supplemental Information 6.

Traits

The F2 plants were greenhouse grown. The plant traits discussed were measured in May 2017, when most of the F2 plants were 5 y old. The traits are illustrated in Fig 4 (parents) and in Fig 5 (F2 examples).

In addition to the traits presented in the results and discussion, there were nine traits which failed to show any genetic component to the phenotypic variance: plant height (plheight); a combination of pitcher width and wing height measured in a profile (winglfpf); the distance from the back of the leaf to the tip of the operculum measured in a profile (pfloperc); the width of the pitcher at its widest point above the mouth opening, periostome, (widtbell); a measure of the location of green in a leaf of mixed color (locgreen); a measure of the color of the veins (veincolr); a measure of the thickness of the periostome lip (periowid); total length from base to crown (lengbscr); and one measure of whether the periostome was more S. rosea or S. psittacina—like (perioinv).

Map construction

R version 3.4.3 2017-11-30 (R Core Team, 2013), and the packages that are associated with this version were used in this analysis. The R package onemap (Margarido et al, 2007) was used to construct the map, with some additional analyses of the linkage map performed by the R package qtl (Broman et al, 2003; Broman & Sen, 2009; Arends et al, 2010). Markers were discarded if they appeared in fewer than 10% of the total genetic lines; F2 genotypes were discarded if they had fewer than 10% of the total number of markers; any marker showing segregation distortion was discarded. The final dataset contained 281 F2 lines (genotypes), and 572 markers. The onemap function group() formed linkage groups, then order_seq() and ripple_seq() were used to form the maps. The onemap package suggested an initial LOD score criterion of 6.1 for the mapping process; we actually used LOD 6.5.

Blast2Go (Conesa et al, 2005; Gotz et al, 2008) analysis of the sequence tags for the mapped SNP markers was performed using version 5.1.12 with a standard annotation pipeline of blast and interpro searches followed by go mapping and annotation. The summary results table (Supplemental Information 7) is organized in genetic map order with the marker and QTL positions indicated on the genetic map.

QTL mapping

The R package qtl (Broman et al, 2003; Broman & Sen, 2009; Arends et al, 2010) was used for mapping the quantitative traits. The qtl package functions scanone() and scantwo() were used to scan for candidate loci, and then makeqtl() and fitqtl() were used to evaluate the qtl models. 2,000 permutations were used to set significance levels. Significant QTLs were identified for 17 of 26 quantitative traits. The ANOVA approach of the fitqtl() function was used to determine the significance of single or multiple loci per trait, and to report the significance of all pairwise interactions among these loci, or pairwise epistasis. Effectplot() was also used to help assess allele effects and dominance; examples are given in Fig S1. The QTL effects were calculated by R/qtl as 0.5 * (rosea_value - psittacina_value).

Figure S1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S1. Examples of dominance scored from effect plots.

AA is an S. psittacina homozygote, AB is a heterozygote, and BB is an S. rosea homozygote.

Source data are available for this figure.

Source Data for Figure S1[LSA-2018-00146_SdataF10.zip]

We tested some of the mqm multiple qtl mapping functions; however, the mqmaugment() function never finished its analysis, which suggests that this approach does not work well with our dataset (Broman in Google Groups R/qtl discussion, 2017). The traits and their scoring are discussed in more detail in the Results section. The package LinkageMapView (Ouellette et al, 2018) works together with the qtl package, and was used to generate the genetic map and QTL diagrams.

Data Availability

The sequencing data from this publication have been deposited to the NCBI SRA database (https://www.ncbi.nlm.nih.gov/bioproject/) and assigned the identifier PRJNA288354 with SRA identifiers SRR8096761 and SRR8096762. The data required for genetic marker and QTL mapping are contained in Supplemental Information 2, 3, 4, 5, 6, 7 accompanying this manuscript.

Supplemental Information 2

The data for genetic mapping of the molecular markers and the quantitative traits in a mapmaker .raw format file suitable for R/onemap (Margarido et al, 2007) or R/qtl (Broman et al, 2003).[Supplemental_Information_2.txt]

Supplemental Information 3

Genetic map data in a format suitable for R/qtl (Broman et al, 2003).[Supplemental_Information_3.txt]

Supplemental Information 4

The data for genetic mapping of the molecular markers and the quantitative traits, including the names of the F2 mapping lines, in .csv format suitable for viewing in a spreadsheet.[Supplemental_Information_4.txt]

Supplemental Information 5

The single nucleotide polymorphisms (SNPs) used in the genetic mapping, including their sequence tags and positions within the sequence tags, 0-index-based, as generated by the STACKs program (Catchen et al, 2013), in .csv format suitable for viewing in a spreadsheet.[Supplemental_Information_5.txt]

Supplemental Information 6

The simple sequence repeats (SSR, microsatellites) used in the genetic mapping including the PCR primers needed, in .csv format suitable for viewing in a spreadsheet.[Supplemental_Information_6.txt]

Supplemental Information 7

Blast2Go (Conesa et al, 2005; Gotz et al, 2008) analysis results for the SNP sequence tags, in .csv format suitable for viewing in a spreadsheet.[Supplemental_Information_7.txt]

Acknowledgements

The project was completed over a number of years and we had help and useful discussions with a number of individuals, including: Jess Stephens, Na Wang, Ed McAssey, Jim Leebens-Mack, Jim Hamrick, Greg Cousins, Mary Jo Godt, Lisa Donovan, Ron Determann, Kelly Dawe, Jenny Cruse-Sanders, Michael Boyd, and John Burke. Two manuscript reviewers provided helpful comments and suggestions. Funding supporting this research was received from a University of Georgia Distinguished University Professorship research fund to RL Malmberg.

Author Contributions

  • RL Malmberg: conceptualization, resources, data curation, software, formal analysis, supervision, funding acquisition, validation, investigation, visualization, methodology, project administration, and writing—original draft, review, and editing.

  • WL Rogers: data curation, investigation, methodology, project administration, and writing—review and editing.

  • MS Alabady: resources, data curation, software, formal analysis, investigation, methodology, and writing—review and editing.

Conflict of Interest Statement

The authors declare that they have no conflicts of interest.

  • Received August 3, 2018.
  • Revision received November 16, 2018.
  • Accepted November 19, 2018.
  • © 2018 Malmberg et al.
Creative Commons logoCreative Commons logohttps://creativecommons.org/licenses/by/4.0/

This article is available under a Creative Commons License (Attribution 4.0 International, as described at https://creativecommons.org/licenses/by/4.0/).

References

  1. ↵
    1. Adlassnig W,
    2. Koller-Peroutka M,
    3. Bauer S,
    4. Koshkin E,
    5. Lendl T,
    6. Lichtscheidl IK
    (2012) Endocytotic uptake of nutrients in carnivorous plants. Plant J 71: 303–313. doi:10.1111/j.1365-313X.2012.04997.x
    OpenUrlCrossRefPubMed
  2. ↵
    1. Alabady MS,
    2. Rogers WL,
    3. Malmberg RL
    (2015) Development of transcriptomic markers for population analysis using restriction site associated RNA sequencing (RARseq). PLoS One 10: e0134855. doi:10.1371/journal.pone.0134855
    OpenUrlCrossRef
  3. ↵
    1. Arends D,
    2. Prins P,
    3. Jansen RC,
    4. Broman KW
    (2010) R/qtl: High-throughput multiple QTL mapping. Bioinformatics 26: 2990–2992. doi:10.1093/bioinformatics/btq565
    OpenUrlCrossRefPubMed
  4. ↵
    1. Badouin H,
    2. Gouzy J,
    3. Grassa CJ,
    4. Murat F,
    5. Staton SE,
    6. Cottret L,
    7. Lelandais-Brière C,
    8. Owens GL,
    9. Carrère S,
    10. Mayjonade B, et al.
    (2017) The sunflower genome provides insights into oil metabolism, flowering and Asterid evolution. Nature 546: 148–152. doi:10.1038/nature22380
    OpenUrlCrossRefPubMed
  5. ↵
    1. Baiser B,
    2. Gotelli NJ,
    3. Buckley HL,
    4. Miller TE,
    5. Ellison AM
    (2012) Geographic variation in network structure of a nearctic aquatic food web. Glob Ecol Biogeogr 21: 579–591. doi:10.1111/j.1466-8238.2011.00705.x
    OpenUrlCrossRef
  6. ↵
    1. Bradshaw WE,
    2. Creelman RA
    (1984) Mutualism between the carnivorous purple pitcher plant and its inhabitants. Am Midl Nat 112: 294–304. doi:10.2307/2425436
    OpenUrlCrossRef
  7. ↵
    1. Brandvain Y,
    2. Kenney AM,
    3. Flagel L,
    4. Coop G,
    5. Sweigart AL
    (2014) Speciation and introgression between Mimulus nasutus and Mimulus guttatus. PLoS Genet 10: e1004410: doi:10.1371/journal.pgen.1004410
    OpenUrlCrossRefPubMed
  8. ↵
    1. Broman in Google Groups R/qtl Discussion
    (2017) R/qtl discussion: mqmaugment running forever. https://groups.google.com/forum/#!searchin/rqtl-disc/mqmaugment$20running$20forever%7Csort:date/rqtl-disc/qtlh7DoM-sk/JjqvZ17QAQAJ
  9. ↵
    1. Broman KW,
    2. Sen S
    (2009) Guide to QTL mapping with R/qtl. 1–396. New York, NY: Springer-Verlag
  10. ↵
    1. Broman KW,
    2. Wu H,
    3. Sen S,
    4. Churchill GA
    (2003) R/qtl: QTL mapping in experimental crosses. Bioinformatics 19: 889–890. doi:10.1093/bioinformatics/btg112
    OpenUrlCrossRefPubMed
  11. ↵
    1. Buckley HL,
    2. Miller TE,
    3. Ellison AM,
    4. Gotelli NJ
    (2003) Reverse latitudinal trends in species richness of pitcher-plant food webs. Ecol Lett 6: 825–829. doi:10.1046/j.1461-0248.2003.00504.x
    OpenUrlCrossRef
  12. ↵
    1. Butler JL,
    2. Ellison AM
    (2007) Nitrogen cycling dynamics in the carnivorous northern pitcher plant, Sarracenia purpurea. Funct Ecol 21: 835–843. doi:10.1111/j.1365-2435.2007.01307.x
    OpenUrlCrossRef
  13. ↵
    1. Catchen J,
    2. Hohenlohe PA,
    3. Bassham S,
    4. Amores A,
    5. Cresko WA
    (2013) Stacks: An analysis tool set for population genomics. Mol Ecol 22: 3124–3140. doi:10.1111/mec.12354
    OpenUrlCrossRefPubMed
  14. ↵
    1. Conesa A,
    2. Gotz S,
    3. Garcia-Gomez JM,
    4. Terol J,
    5. Talon M,
    6. Robles M
    (2005) Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21: 3674–3676. doi:10.1093/bioinformatics/bti610
    OpenUrlCrossRefPubMed
  15. ↵
    1. Darwin F
    1. Darwin CR
    (1888) Insectivorous Plants, Darwin F (ed), 2nd edn, London, United Kingdom: John Murray
  16. ↵
    1. Davis GL,
    2. McMullen MD,
    3. Baysdorfer C,
    4. Musket T,
    5. Grant D,
    6. Staebell M,
    7. Xu G,
    8. Polacco M,
    9. Koster L,
    10. Melia-Hancock S, et al.
    (1999) A maize map standard with sequenced core markers, grass genome reference points and 932 expressed sequence tagged sites (ESTs) in a 1736-locus map. Genetics 152: 1137–1172
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Doebley J
    (1992) Mapping the genes that made maize. Trends Genet 8: 302–307. doi:10.1016/0168-9525(92)90261-2
    OpenUrlCrossRefPubMed
  18. ↵
    1. Ellison AM,
    2. Adamec L
    (eds) (2018) Carnivorous Plants: Physiology, Ecology, and Evolution. Oxford, UK: Oxford University Press
  19. ↵
    1. Ellison AM
    (2006) Nutrient limitation and stoichiometry of carnivorous plants. Plant Biol 8: 740–747. doi:10.1055/s-2006-923956
    OpenUrlCrossRefPubMed
  20. ↵
    1. Ellison AM,
    2. Gotelli NJ
    (2002) Nitrogen availability alters the expression of carnivory in the northern pitcher plant, Sarracenia purpurea. Proc Natl Acad Sci USA 99: 4409–4412. doi:10.1073/pnas.022057199
    OpenUrlAbstract/FREE Full Text
  21. ↵
    1. Elshire RJ,
    2. Glaubitz JC,
    3. Sun Q,
    4. Poland JA,
    5. Kawamoto K,
    6. Buckler ES,
    7. Mitchell SE
    (2011) A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. 8: e62137. doi:10.1371/journal.pone.0062137
  22. ↵
    1. Fishman L,
    2. Sweigart AL,
    3. Kenney AM,
    4. Campbell S
    (2014) Major quantitative trait loci control divergence in critical photoperiod for flowering between selfing and outcrossing species of monkeyflower (Mimulus). New Phytol 201: 1498–1507. doi:10.1111/nph.12618
    OpenUrlCrossRefPubMed
  23. ↵
    1. Gallie DR,
    2. Chang SC
    (1997) Signal transduction in the carnivorous plant Sarracenia purpurea: Regulation of secretory hydrolase expression during development and in response to resources. Plant Physiol 115: 1461–1471. doi:10.1104/pp.115.4.1461
    OpenUrlAbstract
  24. ↵
    1. Garner AG,
    2. Kenney AM,
    3. Fishman L,
    4. Sweigart AL
    (2016) Genetic loci with parent-of-origin effects cause hybrid seed lethality in crosses between Mimulus species. New Phytol 211: 319–331. doi:10.1111/nph.13897
    OpenUrlCrossRefPubMed
  25. ↵
    1. Givnish TJ
    (2015) New evidence on the origin of carnivorous plants. Proc Natl Acad Sci USA 112: 10–11. doi:10.1073/pnas.1422278112
    OpenUrlFREE Full Text
  26. ↵
    1. Gotelli NJ,
    2. Ellison AM
    (2006) Food-web models predict species abundances in response to habitat change. PLos Biol 4: 1869–1873. doi:10.1371/journal.pbio.0040324
    OpenUrlCrossRef
  27. ↵
    1. Gotz S,
    2. Garcia-Gomez JM,
    3. Terol J,
    4. Williams TD,
    5. Nagaraj SH,
    6. Nueda MJ,
    7. Robles M,
    8. Talon M,
    9. Dopazo J,
    10. Conesa A
    (2008) High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res 36: 3420–3435. doi:10.1093/nar/gkn176
    OpenUrlCrossRefPubMed
  28. ↵
    1. Hecht A
    (1949) The Somatic chromosomes of Sarracenia. Bull Torrey Bot Club 76: 7–9. doi:10.2307/2481882
    OpenUrlCrossRef
  29. ↵
    1. Ellison AM,
    2. Adamec L
    1. Horner JD,
    2. Plachno BJ,
    3. Bauer U,
    4. Di Giusto B
    (2018) Chapter 12 attraction of prey. In Carnivorous Plants: Physiology, Ecology, and Evolution, Ellison AM, Adamec L (eds) 157–166. Oxford, United Kingdom: Oxford University Press
  30. ↵
    1. Jurgens A,
    2. El-Sayed AM,
    3. Suckling DM
    (2009) Do carnivorous plants use volatiles for attracting prey insects? Funct Ecol 23: 875–887. doi:10.1111/j.1365-2435.2009.01626.x
    OpenUrlCrossRef
  31. ↵
    1. Karagatzides JD,
    2. Butler JL,
    3. Ellison AM
    (2009) The pitcher plant Sarracenia purpurea can directly acquire organic nitrogen and short-circuit the inorganic nitrogen cycle. PLoS One 4: e6164. doi:10.1371/journal.pone.0006164
    OpenUrlCrossRefPubMed
  32. ↵
    1. Stewart P,
    2. Globig S
    1. Karagatzides JD,
    2. Butler JL,
    3. Ellison AM
    (2012) Pitcher plant Sarracenia purpurea and the inorganic nitrogen cycle. In Plant Physiology, Stewart P, Globig S (eds) 15–28. Palm Bay, FL: Apple Academic Press
  33. ↵
    1. Kenney AM,
    2. Sweigart AL
    (2016) Reproductive isolation and introgression between sympatric Mimulus species. Mol Ecol 25: 2499–2517. doi:10.1111/mec.13630
    OpenUrlCrossRef
  34. ↵
    1. Koopman MM,
    2. Carstens BC
    (2011) The microbial Phyllogeography of the carnivorous plant Sarracenia alata. Microb Ecol 61: 750–758. doi:10.1007/s00248-011-9832-9
    OpenUrlCrossRefPubMed
  35. ↵
    1. Koopman MM,
    2. Fuselier DM,
    3. Hird S,
    4. Carstens BC
    (2010) The carnivorous pale pitcher plant harbors diverse, distinct, and time-dependent bacterial communities. Appl Environ Microbiol 76: 1851–1860. doi:10.1128/aem.02440-09
    OpenUrlAbstract/FREE Full Text
  36. ↵
    1. Margarido GRA,
    2. Souza AP,
    3. Garcia AAF
    (2007) OneMap: Software for genetic mapping in outcrossing species. Hereditas 144: 78–79. doi:10.1111/j.2007.0018-0661.02000.x
    OpenUrlCrossRefPubMed
  37. ↵
    1. Von Mendel G
    (1951) Versuche Über Pflanzen-Hybriden, reprinted in Journal of Heredity, Volume 42, pp 3–4, doi:10.1093/oxfordjournals.jhered.a106151
  38. ↵
    1. Moran JA,
    2. Clarke C,
    3. Gowen BE
    (2012) The use of light in prey capture by the tropical pitcher plant Nepenthes aristolochioides. Plant Signal Behav 7: 957–960. doi:10.4161/psb.20912
    OpenUrlCrossRefPubMed
  39. ↵
    1. Ellison AM,
    2. Adamec L
    1. Naczi RFC
    (2018) Systematics and evolution of Sarraceniaceae. In Carnivorous Plants: Physiology, Ecology, and Evolution, Ellison AM, Adamec L (eds), Chapter 9, pp 105–119. Oxford, UK: Oxford University Press
    OpenUrl
  40. ↵
    1. Ouellette LA,
    2. Reid RW,
    3. Blanchard SG,
    4. Brouwer CR
    (2018) LinkageMapView-rendering high-resolution linkage and QTL maps. Bioinformatics 34: 306–307. doi:10.1093/bioinformatics/btx576
    OpenUrlCrossRef
  41. ↵
    1. Plate L
    (1910) Genetics and evolution. In Festschrift Zum Sechzigsten Geburtstag Richard Hertwigs (in German), pp 536–610. Jena, Germany: Fischer
  42. ↵
    1. R Core Team
    (2017) R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/.
  43. ↵
    1. Rogers WL,
    2. Cruse-Sanders JM,
    3. Determann R,
    4. Malmberg RL
    (2010) Development and characterization of microsatellite markers in Sarracenia L. (pitcher plant) species. Conserv Genet Resour 2: 75–79. doi:10.1007/s12686-009-9165-x
    OpenUrlCrossRefPubMed
  44. ↵
    1. Schaefer HM,
    2. Ruxton GD
    (2008) Fatal attraction: Carnivorous plants roll out the red carpet to lure insects. Biol Lett 4: 153–155. doi:10.1098/rsbl.2007.0607
    OpenUrlAbstract/FREE Full Text
  45. ↵
    1. Schaefer HM,
    2. Ruxton GD
    (2014) Fenestration: A window of opportunity for carnivorous plants. Biol Lett 10: doi:10.1098/rsbl.2014.0134
    OpenUrlCrossRef
  46. ↵
    1. Schnell S,
    2. McPherson DE
    (2011) Sarraceniaceae of North America Poole: Redfern Natural History Productions. Poole, UK
  47. ↵
    1. Smith JM,
    2. Haigh J
    (1974) The hitch-hiking effect of a favourable gene. Genet Res 23: 23. doi:10.1017/S0016672300014634
    OpenUrlCrossRefPubMed
  48. ↵
    1. Srivastava A,
    2. Rogers WL,
    3. Breton CM,
    4. Cai LM,
    5. Malmberg RL
    (2011) Transcriptome analysis of Sarracenia, an insectivorous plant. DNA Res 18: 253–261. doi:10.1093/dnares/dsr014
    OpenUrlCrossRefPubMed
  49. ↵
    1. Stearns FW
    (2010) One hundred years of pleiotropy: A retrospective. Genetics 186: 767–773. doi:10.1534/genetics.110.122549
    OpenUrlAbstract/FREE Full Text
  50. ↵
    1. Stephens JD,
    2. Rogers WL,
    3. Heyduk K,
    4. Cruse-Sanders JM,
    5. Determann RO,
    6. Glenn TC,
    7. Malmberg RL
    (2015) Resolving phylogenetic relationships of the recently radiated carnivorous plant genus Sarracenia using target enrichment. Mol Phylogenet Evol 85: 76–87. doi:10.1016/j.ympev.2015.01.015
    OpenUrlCrossRefPubMed
  51. ↵
    1. Wang H,
    2. Studer AJ,
    3. Zhao Q,
    4. Meeley R,
    5. Doebley JF
    (2015) Evidence that the origin of naked kernels during maize domestication was caused by a single amino acid substitution in tga1. Genetics 200: 965–974. doi:10.1534/genetics.115.175752
    OpenUrlAbstract/FREE Full Text
  52. ↵
    1. Wang J,
    2. Su K,
    3. Guo Y,
    4. Xing H,
    5. Zhao Y,
    6. Liu Z,
    7. Li K,
    8. Guo X
    (2017) Construction of a high-density genetic map for grape using specific length amplified fragment (SLAF) sequencing. PLoS One 12: e0181728. doi:10.1371/journal.pone.0181728
    OpenUrlCrossRef
  53. ↵
    1. Weakley AS
    (2015) Flora of the southern and mid-Atlantic states, working draft of May 2015. University of North Carolina Herbarium, North Carolina Botanical Garden, Chapel Hill, NC. http://www.herbarium.unc.edu/flora.htm.
  54. ↵
    1. Wheeler GL,
    2. Carstens BC
    (2018) Evaluating the adaptive evolutionary convergence of carnivorous plant taxa through functional genomics. PeerJ 6: e4322. doi:10.7717/peerj.4322
    OpenUrlCrossRef
  55. ↵
    1. Wills DM,
    2. Burke JM
    (2007) Quantitative trait locus analysis of the early domestication of sunflower. Genetics 176: 2589–2599. doi:10.1534/genetics.107.075333
    OpenUrlAbstract/FREE Full Text
View Abstract
PreviousNext
Back to top
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for your interest in spreading the word on Life Science Alliance.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
A carnivorous plant genetic map: pitcher/insect-capture QTL on a genetic linkage map of Sarracenia
(Your Name) has sent you a message from Life Science Alliance
(Your Name) thought you would like to see the Life Science Alliance web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Sarracenia genetic linkage map with QTL
Russell L Malmberg, Willie L Rogers, Magdy S Alabady
Life Science Alliance Nov 2018, 1 (6) e201800146; DOI: 10.26508/lsa.201800146

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
Sarracenia genetic linkage map with QTL
Russell L Malmberg, Willie L Rogers, Magdy S Alabady
Life Science Alliance Nov 2018, 1 (6) e201800146; DOI: 10.26508/lsa.201800146
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
Issue Cover

In this Issue

Volume 1, No. 6
December 2018
  • Table of Contents
  • Cover (PDF)
  • About the Cover
  • Masthead (PDF)
Advertisement

Jump to section

  • Article
    • Abstract
    • Introduction
    • Results
    • Discussion
    • Materials and Methods
    • Data Availability
    • Acknowledgements
    • References
  • Figures & Data
  • Info
  • Metrics
  • Reviewer Comments
  • PDF

Subjects

  • Genomics & Functional Genomics
  • Evolution
  • Ecology

EMBO Press LogoRockefeller University Press LogoCold Spring Harbor Logo

Content

  • Home
  • Newest Articles
  • Current Issue
  • Archive
  • Subject Collections

For Authors

  • Submit a Manuscript
  • Author Guidelines
  • License, copyright, Fee

Other Services

  • Alerts
  • Twitter
  • RSS Feeds

More Information

  • Editors & Staff
  • Reviewer Guidelines
  • Feedback
  • Licensing and Reuse
  • Privacy Policy

ISSN: 2575-1077
© 2021 Life Science Alliance LLC

Life Science Alliance is registered as a trademark in the U.S. Patent and Trade Mark Office and in the European Union Intellectual Property Office.