An official website of the United States government
Abstract The development of high‐throughput sequencing technologies is dramatically increasing the use of single nucleotide polymorphisms (SNPs) across the field of genetics, but most parentage studies of wild populations still rely on microsatellites. We developed a bioinformatic pipeline for identifyingSNPpanels that are informative for parentage analysis from restriction site‐associatedDNAsequencing (RADseq) data. This pipeline includes options for analysis with or without a reference genome, and provides methods to maximize genotyping accuracy and select sets of unlinked loci that have high statistical power. We test this pipeline on small populations of Mexican gray wolf and bighorn sheep, for which parentage analyses are expected to be challenging due to low genetic diversity and the presence of many closely related individuals. We compare the results of parentage analysis acrossSNPpanels generated with or without the use of a reference genome, and betweenSNPs and microsatellites. For Mexican gray wolf, we conducted parentage analyses for 30 pups from a single cohort where samples were available from 64% of possible mothers and 53% of possible fathers, and the accuracy of parentage assignments could be estimated because true identities of parents were known a priori based on field data. For bighorn sheep, we conducted maternity analyses for 39 lambs from five cohorts where 77% of possible mothers were sampled, but true identities of parents were unknown. Analyses with and without a reference genome producedSNPpanels with ≥95% parentage assignment accuracy for Mexican gray wolf, outperforming microsatellites at 78% accuracy. Maternity assignments were completely consistent across allSNPpanels for the bighorn sheep, and were 74.4% consistent with assignments from microsatellites. Accuracy and consistency of parentage analysis were not reduced when using as few as 284SNPs for Mexican gray wolf and 142SNPs for bighorn sheep, indicating our pipeline can be used to developSNPgenotyping assays for parentage analysis with relatively small numbers of loci.
more »
« less
High Diversity of Testate Amoebae (Amoebozoa, Arcellinida) Detected by HTS Analyses in a New England Fen using Newly Designed Taxon‐specific Primers
Abstract Testate (shell‐building) amoebae, such as the Arcellinida (Amoebozoa), are useful bioindicators for climate change. Though past work has relied on morphological analyses to characterize Arcellinida diversity, genetic analyses revealed the presence of multiple cryptic species underlying morphospecies. Here, we design and deploy Arcellinida‐specific primers for theSSU‐rDNAgene to assess the community composition on the molecular level in a pilot study of two samplings from a New England fen: (1) 36‐cm horizontal transects and vertical cores; and (2) 26‐m horizontal transects fractioned into four size classes (2–10, 10–35, 35–100, and 100–300 μm). Analyses of these data show the following: (1) a considerable genetic diversity within Arcellinida, much of which comes from morphospecies lacking sequences on GenBank; (2) communities characterized byDNA(i.e. active + quiescent) are distinct from those characterized byRNA(i.e. active, indicator of biomass); (3) active communities on the surface tend to be more similar to one another than to core communities, despite considerable heterogeneity; and (4) analyses of communities fractioned by size find some lineages (OTUs) that are abundant in disjunct size categories, suggesting the possibility of life‐history stages. Together, these data demonstrate the potential of these primers to elucidate the diversity of Arcellinida communities in diverse habitats.
more »
« less
- Award ID(s):
- 1651908
- PAR ID:
- 10457474
- Publisher / Repository:
- Wiley-Blackwell
- Date Published:
- Journal Name:
- Journal of Eukaryotic Microbiology
- Volume:
- 67
- Issue:
- 4
- ISSN:
- 1066-5234
- Page Range / eLocation ID:
- p. 450-462
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Copy number variation (CNV) is a major part of the genetic diversity segregating within populations, but remains poorly understood relative to single nucleotide variation. Here, we report on atRNAligase gene (Migut.N02091;RLG1a) exhibiting unprecedented, and fitness‐relevant,CNVwithin an annual population of the yellow monkeyflowerMimulus guttatus.RLG1a variation was associated with multiple traits in pooled population sequencing (PoolSeq) scans of phenotypic and phenological cohorts. Resequencing of inbred lines revealed intermediate‐frequency three‐copy variants ofRLG1a (trip+;5/35 = 14%), andtrip+lines exhibited elevatedRLG1a expression under multiple conditions.trip+carriers, in addition to being over‐represented in late‐flowering and large‐flowered PoolSeq populations, flowered later under stressful conditions in a greenhouse experiment (p < 0.05). In wild population samples, we discovered an additional rareRLG1a variant (high+) that carries 250–300 copies ofRLG1a totalling ~5.7 Mb (20–40% of a chromosome). In the progeny of ahigh+carrier, Mendelian segregation of diagnostic alleles andqPCR‐based copy counts indicate thathigh+is a single tandem array unlinked to the single‐copyRLG1a locus. In the wild,high+carriers had highest fitness in two particularly dry and/or hot years (2015 and 2017; bothp < 0.01), while single‐copy individuals were twice as fecund as eitherCNVtype in a lush year (2016:p < 0.005). Our results demonstrate fluctuating selection onCNVs affecting phenological traits in a wild population, suggest that planttRNAligases mediate stress‐responsive life‐history traits, and introduce a novel system for investigating the molecular mechanisms of gene amplification.more » « less
-
Abstract Single‐nucleotide polymorphisms (SNPs) are preferred over microsatellite markers in many evolutionary studies, but have only recently been applied to studies of parentage. Evaluations ofSNPs and microsatellites for assigning parentage have mostly focused on special cases that require a relatively large number of heterozygous loci, such as species with low genetic diversity or with complex social structures. We developed 120SNPmarkers from a transcriptome assembled usingRNA‐sequencing of a songbird with the most common avian mating system—social monogamy. We compared the effectiveness of 97 novelSNPs and six previously described microsatellites for assigning paternity in the black‐throated blue warbler,Setophaga caerulescens. We show that the full panel of 97SNPs (meanHo = 0.19) was as powerful for assigning paternity as the panel of multiallelic microsatellites (meanHo = 0.86). Paternity assignments using the two marker types were in agreement for 92% of the offspring. Filtering individual samples by a 50% call rate andSNPs by a 75% call rate maximized the number of offspring assigned with 95% confidence usingSNPs. We also found that the 40 most heterozygousSNPs (meanHo = 0.37) had similar power to assign paternity as the full panel of 97SNPs. These findings demonstrate that a relatively small number of variableSNPs can be effective for parentage analyses in a socially monogamous species. We suggest that the development ofSNPmarkers is advantageous for studies that require high‐throughput genotyping or that plan to address a range of ecological and evolutionary questions.more » « less
-
Dimethylsulfoniopropionate (DMSP) is produced by many species of marine phytoplankton and has been reported to provide a variety of beneficial functions including osmoregulation. Dinoflagellates are recognized as majorDMSPproducers; however, accumulation has been shown to be highly variable in this group. We explored the effect of hyposaline transfer inGambierdiscus belizeanusbetween ecologically relevant salinities (36 and 31) onDMSPaccumulation, Chla, cell growth, and cell volume, over 12 d. Our results showed thatG. belizeanusmaintained an intracellularDMSPcontent of 16.3 pmol cell−1and concentration of 139 mMin both salinities. Although this intracellular concentration was near the median reported for other dinoflagellates, the cellular content achieved byG. belizeanuswas the highest reported of any dinoflagellate thus far, owing mainly to its large size.DMSPlevels were not significantly affected by salinity treatment but did change over time during the experiment. Salinity, however, did have a significant effect on the ratio ofDMSP:Chla, suggesting that salinity transfer ofG. belizeanusinduced a physiological response other thanDMSPadjustment. A survey ofDMSPcontent in a variety ofGambierdiscusspecies and strains revealed relatively highDMSPconcentrations (1.0–16.4 pmol cell−1) as well as high intrageneric and intraspecific variation. We conclude that, althoughDMSPmay not be involved in long‐term (3–12 d) osmoregulation in this species,G. belizeanusand otherGambierdiscusspecies may be important contributors toDMSPproduction in tropical benthic microalgal communities due to their large size and high cellular content.more » « less
-
Abstract Next‐generation sequencing technologies now allow researchers of non‐model systems to perform genome‐based studies without the requirement of a (often unavailable) closely related genomic reference. We evaluated the role of restriction endonuclease (RE) selection in double‐digest restriction‐site‐associatedDNAsequencing (ddRADseq) by generating reduced representation genome‐wide data using four differentREcombinations. Our expectation was thatREselections targeting longer, more complex restriction sites would recover fewer loci thanREwith shorter, less complex sites. We sequenced a diverse sample of non‐model arachnids, including five congeneric pairs of harvestmen (Opiliones) and four pairs of spiders (Araneae). Sample pairs consisted of either conspecifics or closely related congeneric taxa, and in total 26 sample pair analyses were tested. Sequence demultiplexing, read clustering and variant calling were performed in thepyRADprogram. The 6‐base pair cutterEcoRIcombined with methylated site‐specific 4‐base pair cutterMspIproduced, on average, the greatest numbers of intra‐individual loci and shared loci per sample pair. As expected, the number of shared loci recovered for a sample pair covaried with the degree of genetic divergence, estimated with cytochrome oxidase I sequences, although this relationship was non‐linear. Our comparative results will prove useful in guiding protocol selection for ddRADseq experiments on many arachnid taxa where reference genomes, even from closely related species, are unavailable.more » « less
