Multiplexed Shotgun Genotyping Research Articles

First Insights into Population Structure and Genetic Diversity Versus Host Specificity in Trypanorhynch Tapeworms Using Multiplexed Shotgun Genotyping.

Theory predicts relaxed host specificity and high host vagility should contribute to reduced genetic structure in parasites while strict host specificity and low host vagility should increase genetic structure. Though these predictions are intuitive, they have never been explicitly tested in a population genomic framework. Trypanorhynch tapeworms, which parasitize sharks and rays (elasmobranchs) as definitive hosts, are the only order of elasmobranch tapeworms that exhibit considerable variability in their definitive host specificity. This allows for unique combinations of host use and geographic range, making trypanorhynchs ideal candidates for studying how these traits influence population-level structure and genetic diversity. Multiplexed shotgun genotyping (MSG) data sets were generated to characterize component population structure and infrapopulation diversity for a representative of each trypanorhynch suborder: the ray-hosted Rhinoptericola megacantha (Trypanobatoida) and the shark-hosted Callitetrarhynchus gracilis (Trypanoselachoida). Adults of R. megacantha are more host-specific and less broadly distributed than adults of C. gracilis, allowing correlation between these factors and genetic structure. Replicate tapeworm specimens were sequenced from the same host individual, from multiple conspecific hosts within and across geographic regions, and from multiple definitive host species. For R. megacantha, population structure coincided with geography rather than host species. For C. gracilis, limited population structure was found, suggesting a potential link between degree of host specificity and structure. Conspecific trypanorhynchs from the same host individual were found to be as, or more, genetically divergent from one another as from conspecifics from different host individuals. For both species, high levels of homozygosity and positive FIS values were documented.

Multiplexed shotgun genotyping (MSG) data resolve phylogenetic relationships within and among archipelagos in Macaronesian Tolpis.

Plants endemic to oceanic archipelagos are suitable for studying evolution, being isolated on substrates of different ages. Evolution has been recent, rendering traditionally employed sequences insufficiently variable for resolving relationships. This study includes sampling in the genus Tolpis (Asteraceae) from the Azores, Madeira, and Cape Verde, and expands upon an earlier study demonstrating the efficacy of multiplexed shotgun genotyping (MSG) for resolving relationships in Canarian Tolpis. Genomic libraries for 90 accessions of Tolpis and two from the outgroup were generated for genotyping individuals using MSG. Loci were de novo assembled with iPyrad, which clusters MSG loci within and between samples. A maximum likelihood phylogeny was generated with RAxML. Ancestral area reconstruction was inferred using R package BioGeoBEARS. MSG data recovered a highly resolved phylogeny from population to inter-archipelago levels. Ancestral area reconstruction provided biogeographic hypotheses for the radiation of Macaronesian Tolpis. Four major clades were resolved. The Madeiran endemic T. macrorhiza is sister to other Tolpis. Species from the Canaries, Cape Verdes, and the continent are sister to T. succulenta from Madeira, which has a sister subclade of Azorean populations composed of T. succulenta and T. azorica. Population-level resolution suggests unrecognized taxa on several archipelagos. Ancestral reconstruction suggests initial dispersal from the continent to Madeira, with dispersal to the Azores, then dispersal from Madeira to the Canary Islands, with both subsequent dispersal to the Cape Verdes and back-dispersal to the continent. Single-island radiations and inter-island dispersal are implicated in divergence in Macaronesian Tolpis.

Co-evolving wing spots and mating displays are genetically separable traits in Drosophila.

The evolution of sexual traits often involves correlated changes in morphology and behavior. For example, in Drosophila, divergent mating displays are often accompanied by divergent pigment patterns. To better understand how such traits co-evolve, we investigated the genetic basis of correlated divergence in wing pigmentation and mating display between the sibling species Drosophila elegans and Drosophila gunungcola. Drosophila elegans males have an area of black pigment on their wings known as a wing spot and appear to display this spot to females by extending their wings laterally during courtship. By contrast, D. gunungcola lost both of these traits. Using Multiplexed Shotgun Genotyping (MSG), we identified a ∼440 kb region on the X chromosome that behaves like a genetic switch controlling the presence or absence of male-specific wing spots. This region includes the candidate gene optomotor-blind (omb), which plays a critical role in patterning the Drosophila wing. The genetic basis of divergent wing display is more complex, with at least two loci on the X chromosome and two loci on autosomes contributing to its evolution. Introgressing the X-linked region affecting wing spot development from D. gunungcola into D. elegans reduced pigmentation in the wing spots but did not affect the wing display, indicating that these are genetically separable traits. Consistent with this observation, broader sampling of wild D. gunungcola populations confirmed that the wing spot and wing display are evolving independently: some D. gunungcola males performed wing displays similar to D. elegans despite lacking wing spots. These data suggest that correlated selection pressures rather than physical linkage or pleiotropy are responsible for the coevolution of these morphological and behavioral traits. They also suggest that the change in morphology evolved prior to the change in behavior.

Rapid Characterization of the Genetic Loci Controlling Commodity Traits of Chinese Hami Melon (Cucumis melo var. Saccharinensis Naud.) through Multiplexed Shotgun Genotyping

The genetic architecture and the genetic loci controlling commodity traits in this Hami melon have not been characterized. Multiplexed shotgun genotyping (MSG) was used to genotype an F2 population of 370 Chinese Hami melon progeny. A total of 47,609 single nucleotide polymorphism (SNP) markers were obtained after strict filtering. Thebins were used to construct a genetic linkage map with a total length of 1572.954 cM. Quantitative trait locus (QTL) analysis revealed that fruit color was controlled by one major gene about 2 Mb region on chr09, while exocarp color (EC) was controlled by one major gene about 1.9 Mb on chr04, and skin spotting was controlled by two dominant genes, one in the same region of chr04as the EC QTL and the other in the 1031.05 kb region on chr02. Two major QTLs on chr03 and chr05 were related pleiotropically to several quantitative fruit traits, namely, edge sugar content (ES), center sugar content (CS), fruit weight (FW), and fruit length (FL). A further QTL on chr09 also influenced ES, while five other QTLs affected FL. This study was the first to conduct genetic architecture analysis and QTL mapping in Chinese Hami melon with high-density markers and a large target population.

The genetic breakdown of sporophytic self-incompatibility in Tolpis coronopifolia (Asteraceae).

Angiosperm diversity has been shaped by mating system evolution, with the most common transition from outcrossing to self-fertilizing. To investigate the genetic basis of this transition, we performed crosses between two species endemic to the Canary Islands, the self-compatible (SC) species Tolpis coronopifolia and its self-incompatible (SI) relative Tolpis santosii. We scored self-compatibility as self-seed set of recombinant plants within two F2 populations. To map and genetically characterize the breakdown of SI, we built a draft genome sequence of T.coronopifolia, genotyped F2 plants using multiplexed shotgun genotyping (MSG), and located MSG markers to the genome sequence. We identified a single quantitative trait locus (QTL) that explains nearly all variation in self-seed set in both F2 populations. To identify putative causal genetic variants within the QTL, we performed transcriptome sequencing on mature floral tissue from both SI and SC species, constructed a transcriptome for each species, and then located each predicted transcript to the T.coronopifolia genome sequence. We annotated each predicted gene within the QTL and found two strong candidates for SI breakdown. Each gene has a coding sequence insertion/deletion mutation within the SC species that produces a truncated protein. Homologs of each gene have been implicated in pollen development, pollen germination, and pollen tube growth in other species.

Multiplexed shotgun genotyping resolves species relationships within the North American genus Penstemon.

Evolutionary radiations provide opportunities to examine large-scale patterns in diversification and character evolution, yet are often recalcitrant to phylogenetic resolution due to rapid speciation events. The plant genus Penstemon has been difficult to resolve using Sanger sequence-based markers, leading to the hypothesis that it represents a recent North American radiation. The current study demonstrates the utility of multiplexed shotgun genotyping (MSG), a style of restriction site-associated DNA sequencing (RADseq), to infer phylogenetic relationships within a subset of species in this genus and provide insight into evolutionary patterns. We sampled genomic DNA, primarily from herbarium material, and subjected it to MSG library preparation and Illumina sequencing. The resultant sequencing reads were clustered into homologous loci, aligned, and concatenated into data matrices that differed according to clustering similarity and amount of missing data. We performed phylogenetic analyses on these matrices using maximum likelihood (RAxML) and a species tree approach (SVDquartets). MSG data provide a highly resolved estimate of species relationships within Penstemon. While most species relationships were highly supported, the position of certain taxa remains ambiguous, suggesting that increased taxonomic sampling or additional methodologies may be required. The data confirm that evolutionary shifts from hymenopteran- to hummingbird-adapted flowers have occurred independently many times. This study demonstrates that phylogenomic approaches yielding thousands of variable sites can greatly improve species-level resolution of recent and rapid radiations. Similar to other studies, we found that less conservative similarity and missing data thresholds resulted in more highly supported topologies.

A Major Locus Controls a Genital Shape Difference Involved in Reproductive Isolation Between Drosophila yakuba and Drosophila santomea.

Rapid evolution of genitalia shape, a widespread phenomenon in animals with internal fertilization, offers the opportunity to dissect the genetic architecture of morphological evolution linked to sexual selection and speciation. Most quantitative trait loci (QTL) mapping studies of genitalia divergence have focused on Drosophila melanogaster and its three most closely related species, D. simulans, D. mauritiana, and D. sechellia, and have suggested that the genetic basis of genitalia evolution involves many loci. We report the first genetic study of male genitalia evolution between D. yakuba and D. santomea, two species of the D. melanogaster species subgroup. We focus on male ventral branches, which harm females during interspecific copulation. Using landmark-based geometric morphometrics, we characterized shape variation in parental species, F1 hybrids, and backcross progeny and show that the main axis of shape variation within the backcross population matches the interspecific variation between parental species. For genotyping, we developed a new molecular method to perform multiplexed shotgun genotyping (MSG), which allowed us to prepare genomic DNA libraries from 365 backcross individuals in a few days using little DNA. We detected only three QTL, one of which spans 2.7 Mb and exhibits a highly significant effect on shape variation that can be linked to the harmfulness of the ventral branches. We conclude that the genetic architecture of genitalia morphology divergence may not always be as complex as suggested by previous studies.

The Genetics of Resistance to Morinda Fruit Toxin During the Postembryonic Stages in Drosophila sechellia.

Although a great deal has been learned regarding the genetic changes that give rise to adaptation in bacteria and yeast, an understanding of how new complex traits arise in multicellular organisms is far less complete. Many phytophagous insect species are ecological specialists that have adapted to utilize a single host plant. Drosophila sechellia is a specialist that utilizes the ripe fruit of Morinda citrifolia, which is toxic to its sibling species, D. simulans. Here we apply multiplexed shotgun genotyping and QTL analysis to examine the genetic basis of resistance to M. citrifolia fruit toxin in interspecific hybrids. We identify a locus of large effect on the third chromosome (QTL-IIIsima) in the D. simulans backcross that was not detected in previous analyses. We also identify a highly significant QTL of large effect on the X chromosome, QTL-Xsim. Additional smaller-effect loci were also identified in the D. simulans and D. sechellia backcrosses. We did not detect significant epistasis between loci. Instead, our analysis reveals large and smaller-effect loci that contribute to M. citrifolia resistance additively. The additive effect of each locus suggests that partial resistance to lower levels of M. citrifolia toxin could be passed through introgression from D. sechellia to D. simulans in nature. The identification of the major effect loci, QTL-IIIsima and QTL-Xsim, is an important step toward identifying the molecular basis of adaptation in a multicellular organism.

High throughput preparation of fly genomic DNA in 96-well format using a paint-shaker

Sample homogenization is an essential step for genomic DNA extraction, with multiple downstream applications in Molecular Biology. Genotyping hundreds or thousands of samples requires an automation of this homogenization step, and high throughput homogenizer equipment currently costs 7000 euros or more. We present an apparatus for homogenization of individual Drosophila adult flies in 96-well micro-titer dishes, which was built from a small portable paint-shaker (F5 portable paint-shaker, Ushake). Single flies are disrupted in each well that contains extraction buffer and a 4-mm metal ball. Our apparatus can hold up to five 96-well micro-titer plates. Construction of the homogenizer apparatus takes about 3–4 days, and all equipment can be obtained from a home improvement store. The total material cost is approximately 700 euros including the paint-shaker. We tested the performance of our apparatus using the ZR-96 Quick-gDNA™ kit (Zymo Research) homogenization buffer and achieved nearly complete tissue homogenization after 15 minutes of shaking. PCR tests did not detect any cross contamination between samples of neighboring wells. We obtained on average 138 ng of genomic DNA per fly, and DNA quality was adequate for standard PCR applications. In principle, our tissue homogenizer can be used for isolation of DNA suitable for library production and high throughput genotyping by Multiplexed Shotgun Genotyping (MSG), as well as RNA isolation from single flies. The sample adapter can also hold and shake other items, such as centrifuge tubes (15–50 mL) or small bottles.

Multiplexed-shotgun-genotyping data resolve phylogeny within a very recently derived insular lineage.

Endemic plants on oceanic islands have long served as model systems for studying patterns and processes of evolution. However, phylogenetic studies of island plants frequently illustrate a decoupling of molecular divergence and ecological/morphological diversity, resulting in phylogenies lacking the resolution required to interpret patterns of evolution in a phylogenetic context. The current study uses the primarily Macaronesian flowering plant genus Tolpis to illustrate the utility of multiplexed shotgun genotyping (MSG) for resolving relationships at relatively deep (among archipelagos) and very shallow (within archipelagos) nodes in this small, yet diverse insular plant lineage that had not been resolved with other molecular markers. Genomic libraries for 27 accessions of Macaronesian Tolpis were generated for genotyping individuals using MSG, a form of reduced-representation sequencing, similar to restriction-site-associated DNA markers (RADseq). The resulting data files were processed using the program pyRAD, which clusters MSG loci within and between samples. Phylogenetic analyses of the aligned data matrix were conducted using RAxML. Analysis of MSG data recovered a highly resolved phylogeny with generally strong support, including the first robust inference of relationships within the highly diverse Canary Island clade of Tolpis. The current study illustrates the utility of MSG data for resolving relationships in lineages that have undergone recent, rapid diversification resulting in extensive ecological and morphological diversity. We suggest that a similar approach may prove generally useful for other rapid plant radiations where resolving phylogeny has been difficult.

Construction of high-density genetic linkage maps for orange-spotted grouper Epinephelus coioidesusing multiplexed shotgun genotyping

BackgroundOrange-spotted grouper, Epinephelus coioides, is one of the most valuable fish species in China. Commercial production of orange-spotted grouper could be increased by developing higher growth rates and improving commercially important traits. Information on genetic markers associated with quantitative trait loci (QTL) can be used in breeding programs to identify and select individuals carrying desired traits. A high-density genetic linkage map is the basis for QTL study, and multiplexed shotgun genotyping (MSG) facilitates the development of single nucleotide polymorphisms (SNPs) and genotyping. In this study, the first high-density genetic linkage maps for groupers were generated on the basis of the MSG method.ResultsThe sex-averaged map contained a total of 4,608 SNPs, which spanned 1581.7 cM, with a mean distance between SNPs of 0.34 cM. The 4,608 SNPs were located in 2,849 unique locations on the linkage map, with an average inter-location space at 0.56 cM. There were 2,516 SNPs on the female map, and the number of unique locus was 1,902. However, the male map contained more numbers of SNP (2,939) and unique locations (2,005). The total length of the female and male maps was 1,370.9 and 1,335.5 cM, respectively.ConclusionsThe high-resolution genetic linkage maps will be very useful for QTL analyses and marker-assisted selection (MAS) for economically important traits in molecular breeding of the orange-spotted grouper.

Genetic analysis of an elite super-hybrid rice parent using high-density SNP markers.

BackgroundWith an increasing world population and a gradual decline in the amount of arable land, food security remains a global challenge. Continued increases in rice yield will be required to break through the barriers to grain output. In order to transition from hybrid rice to super-hybrid rice, breeding demands cannot be addressed through traditional heterosis. Therefore, it is necessary to incorporate high yield loci from other rice genetic groups and to scientifically utilize intersubspecific heterosis in breeding lines. In this study, 781 lines from a segregating F2 population constructed by crossing the indica variety, “Giant Spike Rice” R1128 as trait donor with the japonica cultivar ‘Nipponbare’, were re-sequenced using high-throughout multiplexed shotgun genotyping (MSG) technology. In combination with high-density single nucleotide polymorphisms, quantitative trait locus (QTL) mapping and genetic effect analysis were performed for five yield factors (spikelet number per panicle, primary branches per panicle, secondary branches per panicle, plant height, and panicle length) to explore the genetic mechanisms underlying the formation of the giant panicle of R1128. Also, they were preformed to locate new high-yielding rice genetic intervals, providing data for super-high-yielding rice breeding.ResultsQTL mapping and genetic effect analysis for five yield factors in the population gave the following results: 49 QTLs for the five yield factors were distributed on 11 of 12 chromosomes. The super-hybrid line R1128 carries multiple major genes for good traits, including Sd1 for plant height, Hd1 and Ehd1 for heading date, Gn1a for spikelet number and IPA1 for ideal plant shape. These genes accounted for 44.3%, 21.9%, 6.2%, 12.9% and 10.6% of the phenotypic variation in the individual traits. Six novel QTLs, qph1-2, qph9-1, qpl12-1, qgn3-1, qgn11-1 and qsbn11-1 are reported here for the first time.ConclusionsHigh-throughout sequencing technology makes it convenient to study rice genomics and makes the QTL/gene mapping direct, efficient, and more reliable. The genetic regions discovered in this study will be valuable for breeding in rice varieties because of the diverse genetic backgrounds of the rice.Electronic supplementary materialThe online version of this article (doi:10.1186/1939-8433-6-21) contains supplementary material, which is available to authorized users.

Evolution of Multiple Additive Loci Caused Divergence between Drosophila yakuba and D. santomea in Wing Rowing during Male Courtship

In Drosophila, male flies perform innate, stereotyped courtship behavior. This innate behavior evolves rapidly between fly species, and is likely to have contributed to reproductive isolation and species divergence. We currently understand little about the neurobiological and genetic mechanisms that contributed to the evolution of courtship behavior. Here we describe a novel behavioral difference between the two closely related species D. yakuba and D. santomea: the frequency of wing rowing during courtship. During courtship, D. santomea males repeatedly rotate their wing blades to face forward and then back (rowing), while D. yakuba males rarely row their wings. We found little intraspecific variation in the frequency of wing rowing for both species. We exploited multiplexed shotgun genotyping (MSG) to genotype two backcross populations with a single lane of Illumina sequencing. We performed quantitative trait locus (QTL) mapping using the ancestry information estimated by MSG and found that the species difference in wing rowing mapped to four or five genetically separable regions. We found no evidence that these loci display epistasis. The identified loci all act in the same direction and can account for most of the species difference.

Multiplexed shotgun genotyping for rapid and efficient genetic mapping

We present a new approach to genotyping based on multiplexed shotgun sequencing that can identify recombination breakpoints in a large number of individuals simultaneously at a resolution sufficient for most mapping purposes, such as quantitative trait locus (QTL) mapping and mapping of induced mutations. We first describe a simple library construction protocol that uses just 10 ng of genomic DNA per individual and makes the approach accessible to any laboratory with standard molecular biology equipment. Sequencing this library results in a large number of sequence reads widely distributed across the genomes of multiplexed bar-coded individuals. We develop a Hidden Markov Model to estimate ancestry at all genomic locations in all individuals using these data. We demonstrate the utility of the approach by mapping a dominant marker allele in D. simulans to within 105 kb of its true position using 96 F1-backcross individuals genotyped in a single lane on an Illumina Genome Analyzer. We further demonstrate the utility of our method by genetically mapping more than 400 previously unassembled D. simulans contigs to linkage groups and by evaluating the quality of targeted introgression lines. At this level of multiplexing and divergence between strains, our method allows estimation of recombination breakpoints to a median of 38-kb intervals. Our analysis suggests that higher levels of multiplexing and/or use of strains with lower levels of divergence are practicable.