Development of a universal double-digest RAD sequencing approach for a group of nonmodel, ecologically and economically important insect and fish taxa

Genetic Variants Underlying Plasticity in Natural Populations of Spadefoot Toads: Environmental Assessment versus Phenotypic Response

Many organisms facultatively produce different phenotypes depending on their environment, yet relatively little is known about the genetic bases of such plasticity in natural populations. In this study, we describe the genetic variation underlying an extreme form of plasticity--resource polyphenism--in Mexican spadefoot toad tadpoles, Spea multiplicata. Depending on their environment, these tadpoles develop into one of two drastically different forms: a carnivore morph or an omnivore morph. We collected both morphs from two ponds that differed in which morph had an adaptive advantage and performed genome-wide association studies of phenotype (carnivore vs. omnivore) and adaptive plasticity (adaptive vs. maladaptive environmental assessment). We identified four quantitative trait loci associated with phenotype and nine with adaptive plasticity, two of which exhibited signatures of minor allele dominance and two of which (one phenotype locus and one adaptive plasticity locus) did not occur as minor allele homozygotes. Investigations into the genetics of plastic traits in natural populations promise to provide novel insights into how such complex, adaptive traits arise and evolve.

Life-history stage and the population genetics of the tiger mosquito Aedes albopictus at a fine spatial scale

As a widespread vector of disease with an expanding range, the mosquito Aedes albopictus Skuse (Diptera: Culicidae) is a high priority for research and management. A. albopictus has a complex life history with aquatic egg, larval and pupal stages, and a terrestrial adult stage. This requires targeted management strategies for each life stage, coordinated across time and space. Population genetics can aid in A. albopictus control by evaluating patterns of genetic diversity and dispersal. However, how life stage impacts population genetic characteristics is unknown. We examined whether patterns of A. albopictus genetic diversity and differentiation changed with life stage at a spatial scale relevant to management efforts. We first conducted a literature review of field-caught A. albopictus population genetic papers and identified 101 peer-reviewed publications, none of which compared results between life stages. Our study uniquely examines population genomic patterns of egg and adult A. albopictus at five sites in Wake County, North Carolina, USA, using 8425 single nucleotide polymorphisms. We found that the level of genetic diversity and connectivity between sites varied between adults and eggs. This warrants further study and is critical for research aimed at informing local management.

Morphometric and Genetic Description of Trophic Adaptations in Cichlid Fishes

Since Darwin, biologists have sought to understand the evolution and origins of phenotypic adaptations. The skull is particularly diverse due to intense natural selection on feeding biomechanics. We investigated the genetic and molecular origins of trophic adaptation using Lake Malawi cichlids, which have undergone an exemplary evolutionary radiation. We analyzed morphological differences in the lateral and ventral head shape among an insectivore that eats by suction feeding, an obligate biting herbivore, and their F2 hybrids. We identified variation in a series of morphological traits—including mandible width, mandible length, and buccal length—that directly affect feeding kinematics and function. Using quantitative trait loci (QTL) mapping, we found that many genes of small effects influence these craniofacial adaptations. Intervals for some traits were enriched in genes related to potassium transport and sensory systems, the latter suggesting co-evolution of feeding structures and sensory adaptations for foraging. Despite these indications of co-evolution of structures, morphological traits did not show covariation. Furthermore, phenotypes largely mapped to distinct genetic intervals, suggesting that a common genetic basis does not generate coordinated changes in shape. Together, these suggest that craniofacial traits are mostly inherited as separate modules, which confers a high potential for the evolution of morphological diversity. Though these traits are not restricted by genetic pleiotropy, functional demands of feeding and sensory structures likely introduce constraints on variation. In all, we provide insights into the quantitative genetic basis of trophic adaptation, identify mechanisms that influence the direction of morphological evolution, and provide molecular inroads to craniofacial variation.

Faster-haplodiploid evolution under divergence-with-gene-flow: simulations and empirical data from pine-feeding hymenopterans

Although haplodiploidy is widespread in nature, the evolutionary consequences of this mode of reproduction are not well characterized. Here, we examine how genome-wide hemizygosity and a lack of recombination in haploid males affects genomic differentiation in populations that diverge via natural selection while experiencing gene flow. First, we simulated diploid and haplodiploid "genomes" (500-kb loci) evolving under an isolation-with-migration model with mutation, drift, selection, migration, and recombination; and examined differentiation at neutral sites both tightly and loosely linked to a divergently selected site. So long as there is divergent selection and migration, sex-limited hemizygosity and recombination cause elevated differentiation (i.e., produce a "faster-haplodiploid effect") in haplodiploid populations relative to otherwise equivalent diploid populations, for both recessive and codominant mutations. Second, we used genome-wide SNP data to model divergence history and describe patterns of genomic differentiation between sympatric populations of Neodiprion lecontei and N. pinetum, a pair of pine sawfly species (order: Hymenoptera; family: Diprionidae) that are specialized on different pine hosts. These analyses support a history of continuous gene exchange throughout divergence and reveal a pattern of heterogeneous genomic differentiation that is consistent with divergent selection on many unlinked loci. Third, using simulations of haplodiploid and diploid populations evolving according to the estimated divergence history of N. lecontei and N. pinetum, we found that divergent selection would lead to higher differentiation in haplodiploids. Based on these results, we hypothesize that haplodiploids undergo divergence-with-gene-flow and sympatric speciation more readily than diploids.

The impact of life stage and pigment source on the evolution of novel warning signal traits

Our understanding of how novel warning color traits evolve in natural populations is largely based on studies of reproductive stages and organisms with endogenously produced pigmentation. In these systems, genetic drift is often required for novel alleles to overcome strong purifying selection stemming from frequency‐dependent predation and positive assortative mating. Here, we integrate data from field surveys, predation experiments, population genomics, and phenotypic correlations to explain the origin and maintenance of geographic variation in a diet‐based larval pigmentation trait in the redheaded pine sawfly (Neodiprion lecontei), a pine‐feeding hymenopteran. Although our experiments confirm that N. lecontei larvae are indeed aposematic—and therefore likely to experience frequency‐dependent predation—our genomic data do not support a historical demographic scenario that would have facilitated the spread of an initially deleterious allele via drift. Additionally, significantly elevated differentiation at a known color locus suggests that geographic variation in larval color is currently maintained by selection. Together, these data suggest that the novel white morph likely spread via selection. However, white body color does not enhance aposematic displays, nor is it correlated with enhanced chemical defense or immune function. Instead, the derived white‐bodied morph is disproportionately abundant on a pine species with a reduced carotenoid content relative to other pine hosts, suggesting that bottom‐up selection via host plants may have driven divergence among populations. Overall, our results suggest that life stage and pigment source can have a substantial impact on the evolution of novel warning signals, highlighting the need to investigate diverse aposematic taxa to develop a comprehensive understanding of color variation in nature.

Relevance of ddRADseq method for species and population delimitation of closely related and widely distributed wolf spiders (Araneae, Lycosidae)

Although species delimitation is often controversial, emerging DNA-based and classical morphology-based methods are rarely compared using large-scale samplings, even less in the case of widely distributed species that have distant, allopatric populations. In the current study, we examined species boundaries within two wolf spider species of the genus Pardosa (Araneae, Lycosidae), P. riparia and P. palustris. Wolf spiders constitute an excellent model for testing the relevance of traditional vs. modern methods in species and population delimitation because several closely related species are distributed over cross-continental geographic ranges. Allopatric populations of the two Pardosa species were sampled across Europe to Far East Russia (latitudinal range > 150°) and several dozen individuals were studied using morphological characters (morphometry of three measures for both sexes, plus five in males only and two in females only), DNA barcoding (COI sequencing) and double-digest restriction site associated DNA sequencing (ddRADseq). The results obtained allow for changing the taxonomic status of two Far East Russian populations to subspecies and ddRADseq proved to be a powerful tool for taxonomic research despite scarce sampling and inherent subjectivity of species delimitation in allopatry. Overall, this study pleads for both multi-criteria and more population-based studies in taxonomy.

Identification of candidate loci for adaptive phenotypic plasticity in natural populations of spadefoot toads

Phenotypic plasticity allows organisms to alter their phenotype in direct response to changes in the environment. Despite growing recognition of plasticity's role in ecology and evolution, few studies have probed plasticity's molecular bases-especially using natural populations. We investigated the genetic basis of phenotypic plasticity in natural populations of spadefoot toads (Spea multiplicata). Spea tadpoles normally develop into an "omnivore" morph that is favored in long-lasting, low-density ponds. However, if tadpoles consume freshwater shrimp or other tadpoles, they can alternatively develop (via plasticity) into a "carnivore" morph that is favored in ephemeral, high-density ponds. By combining natural variation in pond ecology and morph production with population genetic approaches, we identified candidate loci associated with each morph (carnivores vs. omnivores) and loci associated with adaptive phenotypic plasticity (adaptive vs. maladaptive morph choice). Our candidate morph loci mapped to two genes, whereas our candidate plasticity loci mapped to 14 genes. In both cases, the identified genes tended to have functions related to their putative role in spadefoot tadpole biology. Our results thereby form the basis for future studies into the molecular mechanisms that mediate plasticity in spadefoots. More generally, these results illustrate how diverse loci might mediate adaptive plasticity.

Sea Turtle Population Genomic Discovery: Global and Locus-Specific Signatures of Polymorphism, Selection, and Adaptive Potential

In the era of genomics, single-nucleotide polymorphisms (SNPs) have become a preferred molecular marker to study signatures of selection and population structure and to enable improved population monitoring and conservation of vulnerable populations. We apply a SNP calling pipeline to assess population differentiation, visualize linkage disequilibrium, and identify loci with sex-specific genotypes of 45 loggerhead sea turtles (Caretta caretta) sampled from the southeastern coast of the United States, including 42 individuals experimentally confirmed for gonadal sex. By performing reference-based SNP calling in independent runs of Stacks, 3,901–6,998 SNPs and up to 30 potentially sex-specific genotypes were identified. Up to 68 pairs of loci were found to be in complete linkage disequilibrium, potentially indicating regions of natural selection and adaptive evolution. This study provides a valuable SNP diagnostic workflow and a large body of new biomarkers for guiding targeted studies of sea turtle genome evolution and for managing legally protected nonmodel iconic species that have high economic and ecological importance but limited genomic resources.

ddRAD sequencing-based identification of inter-genepool SNPs and association analysis in Brassica juncea

BACKGROUND

Narrow genetic base, complex allo-tetraploid genome and presence of repetitive elements have led the discovery of single nucleotide polymorphisms (SNPs) in Brassica juncea (AABB; 2n = 4x = 36) at a slower pace. Double digest RAD (ddRAD) - a genome complexity reduction technique followed by NGS was used to generate a total of 23 million paired-end reads from three genotypes each of Indian (Pusa Tarak, RSPR-01 and Urvashi) and Exotic (Donskaja IV, Zem 1 and EC287711) genepools.

RESULTS

Sequence data analysis led to the identification of 10,399 SNPs in six genotypes at a read depth of 10x coverage among the genotypes of two genepools. A total of 44 hyper-variable regions (nucleotide variation hotspots) were also found in the genome, of which 93% were found to be a part of coding genes/regions. The functionality of the identified SNPs was estimated by genotyping a subset of SNPs on MassARRAY® platform among a diverse set of B. juncea genotypes. SNP genotyping-based genetic diversity and population studies placed the genotypes into two distinct clusters based mostly on the place of origin. The genotypes were also characterized for six morphological traits, analysis of which revealed a significant difference in the mean values between Indian and Exotic genepools for six traits. The association analysis for six traits identified a total of 45 significant marker-trait associations on 11 chromosomes of A- and B- group of progenitor genomes.

CONCLUSIONS

Despite narrow diversity, the ddRAD sequencing was able to identify large number of nucleotide polymorphisms between the two genepools. Association analysis led to the identification of common SNPs/genomic regions associated between flowering and maturity traits, thereby underscoring the possible role of common chromosomal regions-harboring genes controlling flowering and maturity in Brassica juncea.

Patterns of genetic divergence among populations of Aedes aegypti L. (Diptera: Culicidae) in the southeastern USA

Background

The yellow fever mosquito, Aedes aegypti is a public health concern in the USA, especially in the wake of emergent diseases such as Zika and chikungunya. Aedes aegypti populations dwindled after the invasion of Aedes albopictus in the 1980s and many populations were extirpated. However, in some areas Ae. aegypti persisted in small populations and there are reports of recent resurgences of Ae. aegypti in Florida, Louisiana, Nevada and California. We assessed the population genetic structure of Ae. aegypti in Florida and Georgia, which has concomitant consequences related to mosquito dispersal, pesticide resistance and vectorial capacity.

Methods

We collected Ae. aegypti across Florida and in Georgia using ovitraps. We hatched the eggs and reared them to adults, and after sacrifice we extracted their DNA. We then probed each individual for variation in 6 microsatellite markers, which we used to address population genetic characteristics.

Results

We collected Ae. aegypti and genotyped seven Florida populations and one Georgia population using microsatellite markers. We found evidence of isolation by distance model of gene flow supported by driving distance among cities within Florida and two theoretic genetic clusters.

Conclusions

Significant genetic structure between some populations with substantial gene flow between geographically distant cities suggests regional genetic structuring of Ae. aegypti in Florida. This study provides information on the genetic exchange between populations of Ae. aegypti in the southeastern USA and suggests potential routes of spread of this species.

Mitonuclear discordance in wolf spiders: Genomic evidence for species integrity and introgression

Systematists and taxonomists have benefited greatly from the emergence of molecular methods. Species identification has become straightforward through DNA barcoding and the rapid build-up of massive DNA barcode reference libraries. In animals, mitonuclear discordance can significantly complicate the process of species identification and delimitation. The causes of mitonuclear discordance are either biological (e.g., introgression, incomplete lineage sorting, horizontal gene transfer androgenesis) or induced by operational factors (e.g., human error with specimen misidentification or incorrect species delimitation). Moreover, endosymbionts may play an important role in promoting fixation of mitochondrial genomes. Here, we study the mitonuclear discordance of wolf spiders species (Lycosidae) (independent cases from Alopecosa aculeata and Pardosa pullata groups) that share identical COI DNA barcodes. We approached the case utilizing double-digest restriction site-associated DNA sequencing (ddRADseq) to obtain and analyse genomic-scale data. Our results suggest that the observed cases of mitonuclear discordance are not due to operational reasons but result from biological processes. Further analysis indicated introgression and that incomplete lineage sorting is unlikely to have been responsible for the observed discrepancy. Additional survey of endosymbionts provided ideas on further research and their role in shaping mitochondrial DNA distribution patterns. Thus, ddRADseq grants an efficient way to study the taxonomy of problematic groups with insight into underlying evolutionary processes.

Developing genome-wide SNPs and constructing an ultrahigh-density linkage map in oil palm

Oil palm (Elaeis guineensis Jacq.) is the leading oil-producing crops and the most important edible oil resource worldwide. DNA markers and genetic linkage maps are essential resources for marker-assisted selection to accelerate genetic improvement. We conducted RAD-seq on an Illumina NextSeq500 to discover genome-wide SNPs, and used the SNPs to construct a linkage map for an oil palm (Tenera) population derived from a cross between a Deli Dura and an AVROS Pisifera. The RAD-seq produced 1,076 million single-end reads across the breeding population containing 155 trees. Mining this dataset detected 510,251 loci. After filtering out loci with low accuracy and more than 20% missing data, 11,394 SNPs were retained. Using these SNPs, in combination with 188 anchor SNPs and 123 microsatellites, we constructed a linkage map containing 10,023 markers covering 16 chromosomes. The map length is 2,938.2 cM with an average marker space of 0.29 cM. The large number of SNPs will supply ample choices of DNA markers in analysing the genetic diversity, population structure and evolution of oil palm. This high-density linkage map will contribute to mapping quantitative trait loci (QTL) for important traits, thus accelerating oil palm genetic improvement.

Polygenic sex determination in the cichlid fish Astatotilapia burtoni

BACKGROUND

The East African riverine cichlid species Astatotilapia burtoni serves as an important laboratory model for sexually dimorphic physiology and behavior, and also serves as an outgroup species for the explosive adaptive radiations of cichlid species in Lake Malawi and Lake Victoria. An astounding diversity of genetic sex determination systems have been revealed within the adaptive radiation of East African cichlids thus far, including polygenic sex determination systems involving the epistatic interaction of multiple, independently segregating sex determination alleles. However, sex determination has remained unmapped in A. burtoni. Here we present mapping results supporting the presence of multiple, novel sex determination alleles, and thus the presence of polygenic sex determination in A. burtoni.

RESULTS

Using mapping in small families in conjunction with restriction-site associated DNA sequencing strategies, we identify associations with sex at loci on linkage group 13 and linkage group 5-14. Inheritance patterns support an XY sex determination system on linkage group 5-14 (a chromosome fusion relative to other cichlids studied), and an XYW system on linkage group 13, and these associations are replicated in multiple families. Additionally, combining our genetic data with comparative genomic analysis identifies another fusion that is unassociated with sex, with linkage group 8-24 and linkage group 16-21 fused in A. burtoni relative to other East African cichlid species.

CONCLUSIONS

We identify genetic signals supporting the presence of three previously unidentified sex determination alleles at two loci in the species A. burtoni, strongly supporting the presence of polygenic sex determination system in the species. These results provide a foundation for future mapping of multiple sex determination genes and their interactions. A better understanding of sex determination in A. burtoni provides important context for their use in behavioral studies, as well as studies of the evolution of genetic sex determination and sexual conflicts in East African cichlids.

Development of a universal and simplified ddRAD library preparation approach for SNP discovery and genotyping in angiosperm plants

BACKGROUND

The double digest restriction-site associated DNA sequencing technology (ddRAD-seq) is a reduced representation sequencing technology by sampling genome-wide enzyme loci developed on the basis of next-generation sequencing. ddRAD-seq has been widely applied to SNP marker development and genotyping on animals, especially on marine animals as the original ddRAD protocol is mainly built and trained based on animal data. However, wide application of ddRAD-seq technology in plant species has not been achieved so far. Here, we aim to develop an optimized ddRAD library preparation protocol be accessible to most angiosperm plant species without much startup pre-experiment and costs.

RESULTS

We first tested several combinations of enzymes by in silico analysis of 23 plant species covering 17 families of angiosperm and 1 family of bryophyta and found AvaII + MspI enzyme pair produced consistently higher number of fragments in a broad range of plant species. Then we removed two purifying and one quantifying steps of the original protocol, replaced expensive consumables and apparatuses by conventional experimental apparatuses. Besides, we shortened P1 adapter from 37 to 25 bp and designed a new barcode-adapter system containing 20 pairs of barcodes of varying length. This is an optimized ddRAD strategy for angiosperm plants that is economical, time-saving and requires little technical expertise or investment in laboratory equipment. We refer to this simplified protocol as MiddRAD and we demonstrated the utility and flexibility of our approach by resolving phylogenetic relationships of two genera of woody bamboos (Dendrocalamus and Phyllostachys). Overall our results provide empirical evidence for using this method on different model and non-model plants to produce consistent data.

CONCLUSIONS

As MiddRAD adopts an enzyme pair that works for a broad range of angiosperm plants, simplifies library constructing procedure and requires less DNA input, it will greatly facilitate designing a ddRAD project. Our optimization of this method may make ddRAD be widely used in fields of plant population genetics, phylogenetics, phylogeography and molecular breeding.