Microsatellites for next-generation ecologists: a post-sequencing bioinformatics pipeline

Genetic insights into Cyphocharax magdalenae (Characiformes: Curimatidae): Microsatellite loci development and population analysis in the Cauca River, Colombia

Cyphocharax magdalenae, a Colombian freshwater fish species, plays a vital role in nutrients distribution and serves as a significant food source for other fish species and local fishing communities. Considered a short-distance migratory species, C. magdalenae populations face substantial extinction risk due to human activities impacting their habitats. To address the lack of knowledge on genetic diversity and population structure, this study used next-generation sequencing technology to develop species-specific microsatellite loci and conducted a population genetics analysis of C. magdalenae in the middle and lower sections of the Cauca River, Colombia. Out of 30 pairs of microsatellite primers evaluated in 324 individuals, 14 loci were found to be polymorphic, at linkage equilibrium and, in at least one population, their genotypic frequencies were in Hardy-Weinberg equilibrium. Results showed high genetic diversity levels compared to other neotropical Characiformes, with inbreeding coefficients similar to those reported for phylogenetically related species. Moreover, C. magdalenae exhibits seasonal population structure (rainy-dry) consisting of two genetic stocks showing bottleneck signals and high effective population sizes. This information is essential for understanding the current species genetics and developing future management programs for this fishery resource.

Forest fragmentation causes an isolated population of the golden takin (Budorcas taxicolor bedfordi Thomas, 1911) (Artiodactyla: Bovidae) in the Qinling Mountains (China)

Budorcas taxicolor bedfordi is a rare animal uniquely distributed in the Qinling Mountains (China). Human disturbance and habitat fragmentation have directly affected the survival of B. t. bedfordi. It is urgent to clarify the genetic diversity and genetic structure of the B. t. bedfordi population and implement effective conservation measures. In this study, 20 new polymorphic microsatellite loci were isolated by Illumina sequencing. The genetic diversity and population structure of 124 B. t. bedfordi individuals from three populations (Niubeliang population, Zhouzhi population, and Foping population) were analysed according to these 20 microsatellite loci. Our results indicated that B. t. bedfordi had a low level of genetic variability and that there was inbreeding in the three populations. The population genetic structure analyses showed that the Niubeliang population had a trend of differentiation from other populations. National roads can affect population dispersal, while ecological corridors can promote population gene exchange. None of the three B. t. bedfordi populations experienced bottleneck effects. For conservation management plans, the Zhouzhi population and Foping population should be considered one management unit, and the Niubeliang population should be considered another management unit. We suggest building an ecological corridor to keep the habitat connected and formulating tourism management measures to reduce the influence of human disturbance on B. t. bedfordi.

Transcriptome of the pygmy grasshopper Formosatettix qinlingensis (Orthoptera: Tetrigidae)

Formosatettix qinlingensis (Zheng, 1982) is a tiny grasshopper endemic to Qinling in China. For further study of its transcriptomic features, we obtained RNA-Seq data by Illumina HiSeq X Ten sequencing platform. Firstly, transcriptomic analysis showed that transcriptome read numbers of two female and one male samples were 25,043,314, 24,429,905, and 25,034,457, respectively. We assembled 65,977 unigenes, their average length was 1,072.09 bp, and the length of N50 was 2,031 bp. The average lengths of F. qinlingensis female and male unigenes were 911.30 bp, and 941.82 bp, and the N50 lengths were 1,745 bp and 1,735 bp, respectively. Eight databases were used to annotate the functions of unigenes, and 23,268 functional unigenes were obtained. Besides, we also studied the body color, immunity and insecticide resistance of F. qinlingensis. Thirty-nine pigment-related genes were annotated. Some immunity genes and signaling pathways were found, such as JAK-STAT and Toll-LIKE receptor signaling pathways. There are also some insecticide resistance genes and signal pathways, like nAChR, GST and DDT. Further, some of these genes were differentially expressed in female and male samples, including pigment, immunity and insecticide resistance. The transcriptomic study of F. qinlingensis will provide data reference for gene prediction and molecular expression study of other Tetrigidae species in the future. Differential genetic screening of males and females provides a basis for studying sex and immune balance in insects.

Genome survey sequencing and mining of genome-wide microsatellite markers in yellow-billed babbler (Turdoides affinis)

Turdoides affinis is a species of group dwelling old world passerine of family Leiothrichidae. Unavailability of genome-wide sequence and species-specific molecular markers have hindered comprehensive understanding of cooperative breeding behaviour in T. affinis. Therefore, we generated genome-wide microsatellite markers through whole genome short read sequencing of T. affinis. A total of 68.8 gigabytes of paired-end raw data were sequenced containing 195,067,054 reads. Total sequenced reads spanned a coverage of 17X with genome size of 1.18 Gb. A large number of microsatellite markers (265,297) were mined in the T. affinis genome using Krait, and 50 most informative markers were identified and validated further. In-silico PCR results validated 47 markers. Of these 47 markers, five were randomly selected and validated in-vitro in twelve individuals of T. affinis. Genotyping data on these five loci estimated observed heterozygosity (H₀) and expected heterozygosity (H_e) ratios between 0.333 - 0.833 and 0.851-0.906, respectively. Effective allele size ranged from 6.698 to 10.667, inbreeding coefficient of the population ranged from 0.080 to 0.631 and null allele frequency was calculated at 0.055 to 0.303. Polymorphic information content of all the five loci varied between 0.850 and 0.906. Probabilities of exclusion and identity across 5 loci was estimated to be 0.95 and 0.0036, respectively. All the loci showed significant adherence to Hardy-Weinberg equilibrium. The microsatellite markers reported in this study will facilitate future population genetics studies on T. affinis and other congeneric species.

Whole-Genome Survey and Microsatellite Marker Detection of Antarctic Crocodile Icefish, Chionobathyscus dewitti

Simple Summary

Crocodile icefish inhabit the deep sea around the Southern Ocean and belong to the family Channichthyidae. The species lacks hemoglobin and has evolved an antifreeze protein, unlike other teleosts. In this study, the whole-genome survey and microsatellite motifs were analyzed, which provide relevant information on genetic diversity, population genetics, and the genomic study of crocodile icefish.

Abstract

The crocodile icefish, Chionobathyscus dewitti, belonging to the family Channichthyidae, is an endemic species of the Southern Ocean. The study of its biological features and genetics is challenging as the fish inhabits the deep sea around Antarctic waters. The icefish, the sole cryopelagic species, shows unique physiological and genetic features, unlike other teleosts. It lacks hemoglobin and has evolved antifreeze proteins. Here, we report the genome sequencing data of crocodile icefish produced using the Illumina Novaseq 6000 platform. The estimated genome size was 0.88 Gb with a K-value of 19, and the unique sequence, heterozygosity, error, and duplication rates were 57.4%, 0.421%, 0.317%, and 0.738%, respectively. A genome assembly of 880.69 Mb, with an N50 scaffold length of 2401 bp, was conducted. We identified 2,252,265 microsatellite motifs from the genome assembly data, and dinucleotide repeats (1,920,127; 85.25%) had the highest rate. We selected 84 primer pairs from the genome survey assembly and randomly selected 30 primer pairs for validation. As a result, 15 primer pairs were validated as microsatellite markers.

Evaluation of a new antimicrobial agent production (RSMM C3) by using metagenomics approaches from Egyptian marine biota

•

Metagenomics technique has the ability for production of novel antimicrobial agents.

•

Marine sediment samples from Alexandria used as a source for production of novel antimicrobial agents.

•

Activity of the RSMM C3 antimicrobial agent was a wide spectrum towards different microorganisms.

•

Molecular analysis and characterization of RSMM C3 antimicrobial agent ensure novelty.

Diseases and epidemics in the current days need new types of antibiotics in order to be able to eliminate them. The goal of this research is to use metagenomics to identify isolated utilitarian gene (s) as antimicrobial specialists. Collection of diverse locations from sea sediment samples from Alexandria and extraction of total DNA, restriction enzyme fragmentation, cloning into pUC19 vector, and expression of the isolated gene(s) in E. coli DH5α were all part of the process. Characterization of Antimicrobial agent was done for the best clone for antimicrobial agent's production to detect efficiency, optimum pH, thermal stability, pH stability, effect of different compounds on antimicrobial activity, and residual activity of product after preservation in room temperature. Amino acid sequence of RSMM C3 gene (1250 bp) was 72% identity with Herbaspirillum sp. The ideal temperature level of the RSMM C3 antimicrobial agent production was 36 °C. The antimicrobial agent RSMM C3′s stability was stable at -20 °Celsius for up to two months without thawing. The antibacterial agent RSMM C3 was stable at 4 °C for 14 days without loss in activity. The ideal pH level of the RSMM C3 antimicrobial agent was 6. Remain activity was gradually decreased at pH 5, 6, 6.5 and 7 (86.1, 96.9, 97.2 and 94.9%, respectively). On the other hand, residual activity was (92 and 84%) at (pH 7.5 and 8) for 8 days. The tested antimicrobial RSMM C3 was stable against 1 mM of different compounds (DMSO, Glycerol, NaCl, CaCl₂, MgCl₂, ZnCl₂, FeSO₄, MnSO₄ and CuSO₄). The research provides for the Metagenomics technique that has the ability for the production of novel antimicrobial agents produced by clone RSMM C3 which has a wide spectrum activity towards different microorganisms comparing to other antibiotics as Ampicillin and Tetracycline.

A phylogeny of the genus Limia (Teleostei: Poeciliidae) suggests a single-lake radiation nested in a Caribbean-wide allopatric speciation scenario

OBJECTIVE

The Caribbean is an important global biodiversity hotspot. Adaptive radiations there lead to many speciation events within a limited period and hence are particularly prominent biodiversity generators. A prime example are freshwater fish of the genus Limia, endemic to the Greater Antilles. Within Hispaniola, nine species have been described from a single isolated site, Lake Miragoâne, pointing towards extraordinary sympatric speciation. This study examines the evolutionary history of the Limia species in Lake Miragoâne, relative to their congeners throughout the Caribbean.

RESULTS

For 12 Limia species, we obtained almost complete sequences of the mitochondrial cytochrome b gene, a well-established marker for lower-level taxonomic relationships. We included sequences of six further Limia species from GenBank (total N  = 18 species). Our phylogenies are in concordance with other published phylogenies of Limia. There is strong support that the species found in Lake Miragoâne in Haiti are monophyletic, confirming a recent local radiation. Within Lake Miragoâne, speciation is likely extremely recent, leading to incomplete lineage sorting in the mtDNA. Future studies using multiple unlinked genetic markers are needed to disentangle the relationships within the Lake Miragoâne clade.

Population genetics of the endangered catfish Pseudoplatystoma magdaleniatum (Siluriformes: Pimelodidae) based on species-specific microsatellite loci

ABSTRACT The Neotropical catfish genus Pseudoplatystoma comprises eight species of large size, widely distributed in South American basins. The endangered species P. magdaleniatum is endemic to Magdalena basin (Colombia), experiences high fishing pressure and its population genetics is relatively unknown. To study the genetic status and structure of P. magdaleniatum, 25 species-specific polymorphic microsatellite loci were developed using next-generation sequencing and then tested in samples collected in the Magdalena-Cauca basin. Based on 15 of these loci, P. magdaleniatum showed a high number of alleles per locus (9-10), high values of observed (0.762-0.798) and expected (0.770-0.791) heterozygosities, recent reduction of population size and gene flow. These findings constitute a baseline to measure potential changes in genetic diversity and structure of this commercially important species in a basin undergoing high anthropogenic activities.

Population genetics of the freshwater fish Prochilodus magdalenae (Characiformes: Prochilodontidae), using species-specific microsatellite loci

Prochilodus magdalenae is a freshwater fish endemic to the Colombian Magdalena-Cauca and Caribbean hydrographic basins. The genetic structure patterns of populations of different members of Prochilodus and the historic restocking of its depleted natural populations suggest that P. magdalenae exhibits genetic stocks that coexist and co-migrate throughout the rivers Magdalena, Cauca, Cesar, Sinú and Atrato. To test this hypothesis and explore the levels of genetic diversity and population demography of 725 samples of P. magdalenae from the studied rivers, we developed a set of 11 species-specific microsatellite loci using next-generation sequencing, bioinformatics, and experimental tests of the levels of diversity of the microsatellite loci. The results evidenced that P. magdalenae exhibits high genetic diversity, significant inbreeding coefficient ranging from 0.162 to 0.202, and signs of erosion of the genetic pool. Additionally, the population genetic structure constitutes a mixture of genetic stocks heterogeneously distributed along the studied rivers, and moreover, a highly divergent genetic stock was detected in Chucurí, Puerto Berrío and Palagua that may result from restocking practices. This study provides molecular tools and a wide framework regarding the genetic diversity and structure of P. magdalenae, which is crucial to complement its baseline information, diagnosis and monitoring of populations, and to support the implementation of adequate regulation, management, and conservation policies.

Microsatellite loci development for three catfish species from northwestern South America

ABSTRACT The Neotropical catfish species Ageneiosus pardalis, Pimelodus grosskopfii, and Sorubim cuspicaudus are important fishery resources in Colombia that show historical declines in their capture. This study used next-generation sequencing with 454 FLX technology (Roche Applied Science) and bioinformatics analysis to develop between 18 and 24 microsatellite loci for these species. The novel microsatellite loci showed high values of polymorphic information content -PIC (A. pardalis: 0.601-0.903, P. grosskopfii: 0.748-0.946 and S. cuspicaudus: 0.383-0.876), and the average number of alleles/locus ranged from 7-15 for A. pardalis, 9-30 for P. grosskopfii and 5-14 for S. cuspicaudus. The average observed and expected heterozygosities were respectively, 0.757 ± 0.035 and 0.834 ± 0.015 for A. pardalis; 0.596 ± 0.040 and 0.881 ± 0.009 for P. grosskopfii; and 0.747 ± 0.031 and 0.757 ± 0.025 for S. cuspicaudus. For future studies, these loci can be useful to estimate the genetic diversity and population structure in these three Neotropical catfishes.

A behavioral and genetic study of multiple paternity in a polygamous marine invertebrate, Octopus oliveri

Octopus oliveri is a widespread and common rocky intertidal cephalopod that mates readily in the laboratory, but for which mating behavior has not been reported previously. Four sets of behavioral experiments were recorded wherein three males, small, medium & large in varying order, were introduced to each of six females, for a total of 24 individual females and 12 individual males utilized in the experiments. Video analysis shows that successful mating occurred in each of the mount, reach and beak-to-beak positions. Mating was observed for all males, regardless of size relative to the female, or order of introduction. Females showed preference for the first male to which they were introduced in experimental pairings rather than any specific male trait, and mating time increased significantly with increasing female size. Five novel microsatellite markers were developed and used to test paternity in the eleven broods resulting from these experimental pairings. We found skewed paternity in each brood, with early male precedence and male size being the best predictors of parentage. Multiple paternity was observed in every experimental cross but was estimated to be comparatively low in the field, suggesting that sperm limitation might be common in this species. We saw no evidence of direct sperm competition in Octopus oliveri, but larger males produced significantly more offspring. This study contributes to the growing research on cephalopod mating systems and indicates that octopus mating dynamics might be more variable and complex than thought previously.

Microsatellite development from genome skimming and transcriptome sequencing: comparison of strategies and lessons from frog species

Background

Even though microsatellite loci frequently have been isolated using recently developed next-generation sequencing (NGS) techniques, this task is still difficult because of the subsequent polymorphism screening requires a substantial amount of time. Selecting appropriate polymorphic microsatellites is a critical issue for ecological and evolutionary studies. However, the extent to which assembly strategy, read length, sequencing depth, and library layout produce a measurable effect on microsatellite marker development remains unclear. Here, we use six frog species for genome skimming and two frog species for transcriptome sequencing to develop microsatellite markers, and investigate the effect of different isolation strategies on the yield of microsatellites.

Results

The results revealed that the number of isolated microsatellites increases with increased data quantity and read length. Assembly strategy could influence the yield and the polymorphism of microsatellite development. Larger k-mer sizes produced fewer total number of microsatellite loci, but these loci had a longer repeat length, suggesting greater polymorphism. However, the proportion of each type of nucleotide repeats was not affected; dinucleotide repeats were always the dominant type. Finally, the transcriptomic microsatellites displayed lower levels of polymorphisms and were less abundant than genomic microsatellites, but more likely to be functionally linked loci.

Conclusions

These observations provide deep insight into the evolution and distribution of microsatellites and how different isolation strategies affect microsatellite development using NGS.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5329-y) contains supplementary material, which is available to authorized users.

Microsatellite loci development and population genetics in Neotropical fish Curimata mivartii (Characiformes: Curimatidae)

The Curimatidae family plays an ecological role in the recycling and distribution of nutrients and constitutes a major food source for several commercially important fishes. Curimata mivartii, a member of this family, is considered a short-distance migratory species (≤100 km), categorized by the International Union for Conservation of Nature as a near threatened species, based on its declining population densities and habitat disturbance and fragmentation. Since population genetics and species-specific molecular tools remain unknown for all members of the Curimatidae family, this study developed a set of microsatellite loci and studied the population genetics of C. mivartii in the lower section of the Colombian Magdalena-Cauca basin. The results showed high levels of genetic diversity and evidence of gene flow even between locations separated over 350 km. This information provides a baseline for designing conservation and management programs for C.mivartii and constitutes the first study of population genetics in Curimatidae.

Genetic Diversity and Population Structure of Pythium irregulare from Soybean and Corn Production Fields in Ohio

High levels of genetic diversity have been described within the Pythium irregulare complex from several host plants; however, little is known about the population structure in fields used for grain production. Therefore, the objective of this study was to evaluate the genetic diversity and population structure of 53 isolates baited from 28 soybean and corn production fields from 25 counties in Ohio. Genetic diversity was characterized based on sequence analysis of the internal transcribed spacer (ITS1-5.8S-ITS2) region and with 21 simple sequence repeat (SSR) markers. In addition, aggressiveness on soybean, optimum growth temperature, and sensitivity to metalaxyl fungicide were determined. ITS sequence analysis indicated that four isolates clustered with P. cryptoirregulare, whereas the remaining isolates clustered with P. irregulare that was subdivided into two groups (1 and 2). Cluster analysis of SSR data revealed a similar subdivision, which was also supported by structure analysis. The isolates from group 2 grew at a slower rate, but both groups of P. irregulare and P. cryptoirregulare recovered in this study had the same optimum growth at 27°C. Variability of aggressiveness and sensitivity toward metalaxyl fungicide was also observed among isolates within each group. The results from this study will help in the selection of isolates to be used in screening for resistance, assessment of fungicide efficacy, and disease management recommendations.

Phoretic dispersal influences parasite population genetic structure

Dispersal is a fundamental component of the life history of most species. Dispersal influences fitness, population dynamics, gene flow, genetic drift and population genetic structure. Even small differences in dispersal can alter ecological interactions and trigger an evolutionary cascade. Linking such ecological processes with evolutionary patterns is difficult, but can be carried out in the proper comparative context. Here, we investigate how differences in phoretic dispersal influence the population genetic structure of two different parasites of the same host species. We focus on two species of host-specific feather lice (Phthiraptera: Ischnocera) that co-occur on feral rock pigeons (Columba livia). Although these lice are ecologically very similar, "wing lice" (Columbicola columbae) disperse phoretically by "hitchhiking" on pigeon flies (Diptera: Hippoboscidae), while "body lice" (Campanulotes compar) do not. Differences in the phoretic dispersal of these species are thought to underlie observed differences in host specificity, as well as the degree of host-parasite cospeciation. These ecological and macroevolutionary patterns suggest that body lice should exhibit more genetic differentiation than wing lice. We tested this prediction among lice on individual birds and among lice on birds from three pigeon flocks. We found higher levels of genetic differentiation in body lice compared to wing lice at two spatial scales. Our results indicate that differences in phoretic dispersal can explain microevolutionary differences in population genetic structure and are consistent with macroevolutionary differences in the degree of host-parasite cospeciation.

A biting commentary: Integrating tooth characters with molecular data doubles known species diversity in a lineage of sea slugs that consume "killer algae"

Predicting biotic resistance to highly invasive strains of "killer algae" (Caulerpa spp.) requires understanding the diversity and feeding preferences of native consumers, including sea slugs in family Oxynoidae. Past studies reported low algal host specificity for Oxynoe (6 spp.) and Lobiger (4 spp.), but these taxonomically challenging slugs may represent species complexes of unrecognized specialists that prefer different Caulerpa spp. Here, we assess global diversity of these genera by integrating gene sequences with morphological data from microscopic teeth and internal shells, the only hard parts in these soft-bodied invertebrates. Four delimitation methods applied to datasets comprising mtDNA and/or nuclear alleles yielded up to 16 species hypotheses for samples comprising five nominal taxa, including five highly divergent species in Lobiger and five in Oxynoe. Depending on the analysis, a further four to six species were recovered in the O. antillarum-viridis complex, a clade in which mitochondrial divergence was low and nuclear alleles were shared among lineages. Bayesian species delimitation using only morphological data supported most candidate species, however, and integrative analyses combining morphological and genetic data fully supported all complex members. Collectively, our findings double the recognized biodiversity in Oxynoidae, and illustrate the value of including data from traits that mediate fast-evolving ecological interactions during species delimitation. Preference for Caulerpa spp. and radular tooth characteristics covaried among newly delimited species, highlighting an unappreciated degree of host specialization and coevolution in these taxa that may help predict their role in containing outbreaks of invasive algae.

Clone wars: asexual reproduction dominates in the invasive range of Tubastraea spp. (Anthozoa: Scleractinia) in the South-Atlantic Ocean

Although the invasive azooxanthellate corals Tubastraea coccinea and T. tagusensis are spreading quickly and outcompeting native species in the Atlantic Ocean, there is little information regarding the genetic structure and path of introduction for these species. Here we present the first data on genetic diversity and clonal structure from these two species using a new set of microsatellite markers. High proportions of clones were observed, indicating that asexual reproduction has a major role in the local population dynamics and, therefore, represents one of the main reasons for the invasion success. Although no significant population structure was found, results suggest the occurrence of multiple invasions for T. coccinea and also that both species are being transported along the coast by vectors such as oil platforms and monobouys, spreading these invasive species. In addition to the description of novel microsatellite markers, this study sheds new light into the invasive process of Tubastraea.

Development and characterization of genic SSR markers from low depth genome sequence of Clarias batrachus (magur)

Indian magur (Clarias batrachus) is an important freshwater catfish, which is listed as endangered under A3cde+4acde ver. 3.1 categories by the IUCN (2015) due to decreasing population trend. Microsatellites or short sequence repeats (SSRs) tagged to genes have been utilized as gene marker. In the present study, 31,814 SSRs of C. batrachus (magur) were identified using microsatellite identification tool programme from the next generation sequencing data generated on Roche 454 and Ion Torrent platforms. A bioinformatics pipeline, with stringent criteria resulted in selection of 1672 microsatellite loci falling in the genic region. Initially, a total of 30 loci were selected for primer development; and of these 14 were successfully amplified and five were found to be polymorphic in 30 individuals of C. batrachus (magur). The observed as well as expected heterozygosity ranged from 0.038 to 0.526 and 0.434 to 0.784, respectively, and the number of observed alleles ranged from three to five. The study reported the application of next generation sequencing technologies for rapid development of microsatellite loci in Indian catfish species, C. batrachus (magur).

Characterization, development and multiplexing of microsatellite markers in three commercially exploited reef fish and their application for stock identification

Thirty-four microsatellite loci were isolated from three reef fish species; golden snapper Lutjanus johnii, blackspotted croaker Protonibea diacanthus and grass emperor Lethrinus laticaudis using a next generation sequencing approach. Both IonTorrent single reads and Illumina MiSeq paired-end reads were used, with the latter demonstrating a higher quality of reads than the IonTorrent. From the 1–1.5 million raw reads per species, we successfully obtained 10–13 polymorphic loci for each species, which satisfied stringent design criteria. We developed multiplex panels for the amplification of the golden snapper and the blackspotted croaker loci, as well as post-amplification pooling panels for the grass emperor loci. The microsatellites characterized in this work were tested across three locations of northern Australia. The microsatellites we developed can detect population differentiation across northern Australia and may be used for genetic structure studies and stock identification.

Development and characterization of 24 polymorphic microsatellite loci for the freshwater fish Ichthyoelephas longirostris (Characiformes: Prochilodontidae)

The Neotropical freshwater fish Ichthyoelephas longirostris (Characiformes: prochilodontidae) is a short-distance migratory species endemic to Colombia. This study developed for the first time a set of 24 polymorphic microsatellite loci by using next-generation sequencing to explore the population genetics of this commercially exploited species. Nineteen of these loci were used to assess the genetic diversity and structure of 193 I. longirostris in three Colombian rivers of the Magdalena basin. Results showed that a single genetic stock circulates in the Cauca River, whereas other single different genetic stock is present in the rivers Samaná Norte and San Bartolomé-Magdalena. Additionally, I. longirostris was genetically different among and across rivers. This first insight about the population genetic structure of I. longirostris is crucial for monitoring the genetic diversity, the management and conservation of its populations, and complement the genetic studies in Prochilodontidae.

De Novo Assembly and Characterization of the Transcriptome of Grasshopper Shirakiacris shirakii

Background: The grasshopper Shirakiacris shirakii is an important agricultural pest and feeds mainly on gramineous plants, thereby causing economic damage to a wide range of crops. However, genomic information on this species is extremely limited thus far, and transcriptome data relevant to insecticide resistance and pest control are also not available. Methods: The transcriptome of S. shirakii was sequenced using the Illumina HiSeq platform, and we de novo assembled the transcriptome. Results: Its sequencing produced a total of 105,408,878 clean reads, and the de novo assembly revealed 74,657 unigenes with an average length of 680 bp and N50 of 1057 bp. A total of 28,173 unigenes were annotated for the NCBI non-redundant protein sequences (Nr), NCBI non-redundant nucleotide sequences (Nt), a manually-annotated and reviewed protein sequence database (Swiss-Prot), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Based on the Nr annotation results, we manually identified 79 unigenes encoding cytochrome P450 monooxygenases (P450s), 36 unigenes encoding carboxylesterases (CarEs) and 36 unigenes encoding glutathione S-transferases (GSTs) in S. shirakii. Core RNAi components relevant to miroRNA, siRNA and piRNA pathways, including Pasha, Loquacious, Argonaute-1, Argonaute-2, Argonaute-3, Zucchini, Aubergine, enhanced RNAi-1 and Piwi, were expressed in S. shirakii. We also identified five unigenes that were homologous to the Sid-1 gene. In addition, the analysis of differential gene expressions revealed that a total of 19,764 unigenes were up-regulated and 4185 unigenes were down-regulated in larvae. In total, we predicted 7504 simple sequence repeats (SSRs) from 74,657 unigenes. Conclusions: The comprehensive de novo transcriptomic data of S. shirakii will offer a series of valuable molecular resources for better studying insecticide resistance, RNAi and molecular marker discovery in the transcriptome.

Characterization and analysis of a de novo transcriptome from the pygmy grasshopper Tetrix japonica

The pygmy grasshopper Tetrix japonica is a common insect distributed throughout the world, and it has the potential for use in studies of body colour polymorphism, genomics and the biology of Tetrigoidea (Insecta: Orthoptera). However, limited biological information is available for this insect. Here, we conducted a de novo transcriptome study of adult and larval T. japonica to provide a better understanding of its gene expression and develop genomic resources for future work. We sequenced and explored the characteristics of the de novo transcriptome of T. japonica using Illumina HiSeq 2000 platform. A total of 107 608 206 paired-end clean reads were assembled into 61 141 unigenes using the trinity software; the mean unigene size was 771 bp, and the N50 length was 1238 bp. A total of 29 225 unigenes were functionally annotated to the NCBI nonredundant protein sequences (Nr), NCBI nonredundant nucleotide sequences (Nt), a manually annotated and reviewed protein sequence database (Swiss-Prot), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. A large number of putative genes that are potentially involved in pigment pathways, juvenile hormone (JH) metabolism and signalling pathways were identified in the T. japonica transcriptome. Additionally, 165 769 and 156 796 putative single nucleotide polymorphisms occurred in the adult and larvae transcriptomes, respectively, and a total of 3162 simple sequence repeats were detected in this assembly. This comprehensive transcriptomic data for T. japonica will provide a usable resource for gene predictions, signalling pathway investigations and molecular marker development for this species and other pygmy grasshoppers.

Microsatellite markers for the notothenioid fish Lepidonotothen nudifrons and two congeneric species

BACKGROUND

Loss of genetic variability due to environmental changes, limitation of gene flow between pools of individuals or putative selective pressure at specific markers, were previously documented for Antarctic notothenioid fish species. However, so far no studies were performed for the Gaudy notothen Lepidonotothen nudifrons. Starting from a species-specific spleen transcriptome library, we aimed at isolating polymorphic microsatellites (Type I; i.e. derived from coding sequences) suitable to quantify the genetic variability in this species, and additionally to assess the population genetic structure and demography in nototheniids.

RESULTS

We selected 43,269 transcripts resulting from a MiSeq sequencer run, out of which we developed 19 primer pairs for sequences containing microsatellite repeats. Sixteen loci were successfully amplified in L. nudifrons. Eleven microsatellites were polymorphic and allele numbers per locus ranged from 2 to 17. In addition, we amplified loci identified from L. nudifrons in two other congeneric species (L. squamifrons and L. larseni). Thirteen loci were highly transferable to the two congeneric species. Differences in polymorphism among species were detected.

CONCLUSIONS

Starting from a transcriptome of a non-model organism, we were able to identify promising polymorphic nuclear markers that are easily transferable to other closely related species. These markers can be a key instrument to monitor the genetic structure of the three Lepidonotothen species if genotyped in larger population samples. When compared with anonymous loci isolated in other notothenioids, i.e. Type II (isolated from genomic libraries), they offer the possibility to test how the effects of occurring environmental change influence the population genetic structure in each species and subsequently the composition of the entire ecosystem.

De Novo Assembly and Characterization of Four Anthozoan (Phylum Cnidaria) Transcriptomes

Many nonmodel species exemplify important biological questions but lack the sequence resources required to study the genes and genomic regions underlying traits of interest. Reef-building corals are famously sensitive to rising seawater temperatures, motivating ongoing research into their stress responses and long-term prospects in a changing climate. A comprehensive understanding of these processes will require extending beyond the sequenced coral genome (Acropora digitifera) to encompass diverse coral species and related anthozoans. Toward that end, we have assembled and annotated reference transcriptomes to develop catalogs of gene sequences for three scleractinian corals (Fungia scutaria, Montastraea cavernosa, Seriatopora hystrix) and a temperate anemone (Anthopleura elegantissima). High-throughput sequencing of cDNA libraries produced ~20-30 million reads per sample, and de novo assembly of these reads produced ~75,000-110,000 transcripts from each sample with size distributions (mean ~1.4 kb, N50 ~2 kb), comparable to the distribution of gene models from the coral genome (mean ~1.7 kb, N50 ~2.2 kb). Each assembly includes matches for more than half the gene models from A. digitifera (54-67%) and many reasonably complete transcripts (~5300-6700) spanning nearly the entire gene (ortholog hit ratios ≥0.75). The catalogs of gene sequences developed in this study made it possible to identify hundreds to thousands of orthologs across diverse scleractinian species and related taxa. We used these sequences for phylogenetic inference, recovering known relationships and demonstrating superior performance over phylogenetic trees constructed using single mitochondrial loci. The resources developed in this study provide gene sequences and genetic markers for several anthozoan species. To enhance the utility of these resources for the research community, we developed searchable databases enabling researchers to rapidly recover sequences for genes of interest. Our analysis of de novo assembly quality highlights metrics that we expect will be useful for evaluating the relative quality of other de novo transcriptome assemblies. The identification of orthologous sequences and phylogenetic reconstruction demonstrates the feasibility of these methods for clarifying the substantial uncertainties in the existing scleractinian phylogeny.

Temporal Stability of Genetic Structure in a Mesopelagic Copepod

Although stochasticity in oceanographic conditions is known to be an important driver of temporal genetic change in many marine species, little is known about whether genetically distinct plankton populations can persist in open ocean habitats. A prior study demonstrated significant population genetic structure among oceanic gyres in the mesopelagic copepod Haloptilus longicornis in both the Atlantic and Pacific Oceans, and we hypothesized that populations within each gyre represent distinct gene pools that persist over time. We tested this expectation through basin-scale sampling across the Atlantic Ocean in 2010 and 2012. Using both mitochondrial (mtCOII) and microsatellite markers (7 loci), we show that the genetic composition of populations was stable across two years in both the northern and southern subtropical gyres. Genetic variation in this species was partitioned among ocean gyres (F_CT = 0.285, P < 0.0001 for mtCOII, F_CT = 0.013, P < 0.0001 for microsatellites), suggesting strong spatial population structure, but no significant partitioning was found among sampling years. This temporal persistence of population structure across a large geographic scale was coupled with chaotic genetic patchiness at smaller spatial scales, but the magnitude of genetic differentiation was an order of magnitude lower at these smaller scales. Our results demonstrate that genetically distinct plankton populations persist over time in highly-dispersive open ocean habitats, and this is the first study to rigorously test for temporal stability of large scale population structure in the plankton.

Multilocus evidence for globally distributed cryptic species and distinct populations across ocean gyres in a mesopelagic copepod

Zooplanktonic taxa have a greater number of distinct populations and species than might be predicted based on their large population sizes and open-ocean habitat, which lacks obvious physical barriers to dispersal and gene flow. To gain insight into the evolutionary mechanisms driving genetic diversification in zooplankton, we developed eight microsatellite markers to examine the population structure of an abundant, globally distributed mesopelagic copepod, Haloptilus longicornis, at 18 sample sites across the Atlantic and Pacific Oceans (n = 761). When comparing our microsatellite results with those of a prior study that used a mtDNA marker (mtCOII, n = 1059, 43 sample sites), we unexpectedly found evidence for the presence of a cryptic species pair. These species were globally distributed and apparently sympatric, and were separated by relatively weak genetic divergence (reciprocally monophyletic mtCOII lineages 1.6% divergent; microsatellite FST ranging from 0.28 to 0.88 across loci, P < 0.00001). Using both mtDNA and microsatellite data for the most common of the two species (n = 669 for microsatellites, n = 572 for mtDNA), we also found evidence for allopatric barriers to gene flow within species, with distinct populations separated by continental landmasses and equatorial waters in both the Atlantic and Pacific Ocean basins. Our study shows that oceanic barriers to gene flow can act as a mechanism promoting allopatric diversification in holoplanktonic taxa, despite the high potential dispersal abilities and pelagic habitat for these species.

Multiplexing with three-primer PCR for rapid and economical microsatellite validation

The next generation sequencing revolution has enabled rapid discovery of genetic markers, however, development of fully functioning new markers still requires a long and costly process of marker validation. This study reports a rapid and economical approach for the validation and deployment of polymorphic microsatellite markers obtained from a 454 pyrosequencing library of Atlantic cod, Gadus morhua, Linnaeus 1758. Primers were designed from raw reads to amplify specific amplicon size ranges, allowing effective PCR multiplexing. Multiplexing was combined with a three-primer PCR approach using four universal tails to label amplicons with separate fluorochromes. A total of 192 primer pairs were tested, resulting in 73 polymorphic markers. Of these, 55 loci were combined in six multiplex panels each containing between six and eleven markers. Variability of the loci was assessed on G. morhua from the Celtic Sea (n = 46) and the Scotian Shelf (n = 46), two locations that have shown genetic differentiation in previous studies. Multilocus F(ST) between the two samples was estimated at 0.067 (P = 0.001). After three loci potentially under selection were excluded, the global F(ST) was estimated at 0.043 (P = 0.001). Our technique combines three-primer and multiplex PCR techniques, allowing simultaneous screening and validation of relatively large numbers of microsatellite loci.

Challenges in analysis and interpretation of microsatellite data for population genetic studies

Advancing technologies have facilitated the ever-widening application of genetic markers such as microsatellites into new systems and research questions in biology. In light of the data and experience accumulated from several years of using microsatellites, we present here a literature review that synthesizes the limitations of microsatellites in population genetic studies. With a focus on population structure, we review the widely used fixation (F ST) statistics and Bayesian clustering algorithms and find that the former can be confusing and problematic for microsatellites and that the latter may be confounded by complex population models and lack power in certain cases. Clustering, multivariate analyses, and diversity-based statistics are increasingly being applied to infer population structure, but in some instances these methods lack formalization with microsatellites. Migration-specific methods perform well only under narrow constraints. We also examine the use of microsatellites for inferring effective population size, changes in population size, and deeper demographic history, and find that these methods are untested and/or highly context-dependent. Overall, each method possesses important weaknesses for use with microsatellites, and there are significant constraints on inferences commonly made using microsatellite markers in the areas of population structure, admixture, and effective population size. To ameliorate and better understand these constraints, researchers are encouraged to analyze simulated datasets both prior to and following data collection and analysis, the latter of which is formalized within the approximate Bayesian computation framework. We also examine trends in the literature and show that microsatellites continue to be widely used, especially in non-human subject areas. This review assists with study design and molecular marker selection, facilitates sound interpretation of microsatellite data while fostering respect for their practical limitations, and identifies lessons that could be applied toward emerging markers and high-throughput technologies in population genetics.

tropiTree: an NGS-based EST-SSR resource for 24 tropical tree species

The development of genetic tools for non-model organisms has been hampered by cost, but advances in next-generation sequencing (NGS) have created new opportunities. In ecological research, this raises the prospect for developing molecular markers to simultaneously study important genetic processes such as gene flow in multiple non-model plant species within complex natural and anthropogenic landscapes. Here, we report the use of bar-coded multiplexed paired-end Illumina NGS for the de novo development of expressed sequence tag-derived simple sequence repeat (EST-SSR) markers at low cost for a range of 24 tree species. Each chosen tree species is important in complex tropical agroforestry systems where little is currently known about many genetic processes. An average of more than 5,000 EST-SSRs was identified for each of the 24 sequenced species, whereas prior to analysis 20 of the species had fewer than 100 nucleotide sequence citations. To make results available to potential users in a suitable format, we have developed an open-access, interactive online database, tropiTree (http://bioinf.hutton.ac.uk/tropiTree), which has a range of visualisation and search facilities, and which is a model for the efficient presentation and application of NGS data.

QDD version 3.1: a user-friendly computer program for microsatellite selection and primer design revisited: experimental validation of variables determining genotyping success rate

Microsatellite marker development has been greatly simplified by the use of high-throughput sequencing followed by in silico microsatellite detection and primer design. However, the selection of markers designed by the existing pipelines depends either on arbitrary criteria, or older studies on PCR success. Based on wet laboratory experiments, we have identified the following factors that are most likely to influence genotyping success rate: alignment score between the primers and the amplicon; the distance between primers and microsatellites; the length of the PCR product; target region complexity and the number of reads underlying the sequence. The QDD pipeline has been modified to include these most pertinent factors in the output to help the selection of markers. Furthermore, new features are also included in the present version: (i) not only raw sequencing reads are accepted as input, but also contigs, allowing the analysis of assembled high-coverage data; (ii) input data can be both in fasta and fastq format to facilitate the use of Illumina and IonTorrent reads; (iii) A comparison to known transposable elements allows their detection; (iv) A contamination check can be carried out by BLASTing potential markers against the nucleotide (nt) database of NCBI; (v) QDD3 is now also available imbedded into a virtual machine making installation easier and operating system independent. It can be used both on command-line version as well as integrated into a Galaxy server, providing a user-friendly interface, as well as the possibility to utilize a large variety of NGS tools.

The effects of read length, quality and quantity on microsatellite discovery and primer development: from Illumina to PacBio

The advent of next-generation sequencing (NGS) technologies has transformed the way microsatellites are isolated for ecological and evolutionary investigations. Recent attempts to employ NGS for microsatellite discovery have used the 454, Illumina, and Ion Torrent platforms, but other methods including single-molecule real-time DNA sequencing (Pacific Biosciences or PacBio) remain viable alternatives. We outline a workflow from sequence quality control to microsatellite marker validation in three plant species using PacBio circular consensus sequencing (CCS). We then evaluate the performance of PacBio CCS in comparison with other NGS platforms for microsatellite isolation, through simulations that focus on variations in read length, read quantity and sequencing error rate. Although quality control of CCS reads reduced microsatellite yield by around 50%, hundreds of microsatellite loci that are expected to have improved conversion efficiency to functional markers were retrieved for each species. The simulations quantitatively validate the advantages of long reads and emphasize the detrimental effects of sequencing errors on NGS-enabled microsatellite development. In view of the continuing improvement in read length on NGS platforms, sequence quality and the corresponding strategies of quality control will become the primary factors to consider for effective microsatellite isolation. Among current options, PacBio CCS may be optimal for rapid, small-scale microsatellite development due to its flexibility in scaling sequencing effort, while platforms such as Illumina MiSeq will provide cost-efficient solutions for multispecies microsatellite projects.

DNA fingerprinting in zoology: past, present, future

In 1962, Thomas Kuhn famously argued that the progress of scientific knowledge results from periodic 'paradigm shifts' during a period of crisis in which new ideas dramatically change the status quo. Although this is generally true, Alec Jeffreys' identification of hypervariable repeat motifs in the human beta-globin gene, and the subsequent development of a technology known now as 'DNA fingerprinting', also resulted in a dramatic shift in the life sciences, particularly in ecology, evolutionary biology, and forensics. The variation Jeffreys recognized has been used to identify individuals from tissue samples of not just humans, but also of many animal species. In addition, the technology has been used to determine the sex of individuals, as well as paternity/maternity and close kinship. We review a broad range of such studies involving a wide diversity of animal species. For individual researchers, Jeffreys' invention resulted in many ecologists and evolutionary biologists being given the opportunity to develop skills in molecular biology to augment their whole organism focus. Few developments in science, even among the subsequent genome discoveries of the 21st century, have the same wide-reaching significance. Even the later development of PCR-based genotyping of individuals using microsatellite repeats sequences, and their use in determining multiple paternity, is conceptually rooted in Alec Jeffreys' pioneering work.

De novo assembly, gene annotation, and marker discovery in stored-product pest Liposcelis entomophila (Enderlein) using transcriptome sequences

Background

As a major stored-product pest insect, Liposcelis entomophila has developed high levels of resistance to various insecticides in grain storage systems. However, the molecular mechanisms underlying resistance and environmental stress have not been characterized. To date, there is a lack of genomic information for this species. Therefore, studies aimed at profiling the L. entomophila transcriptome would provide a better understanding of the biological functions at the molecular levels.

Methodology/Principal Findings

We applied Illumina sequencing technology to sequence the transcriptome of L. entomophila. A total of 54,406,328 clean reads were obtained and that de novo assembled into 54,220 unigenes, with an average length of 571 bp. Through a similarity search, 33,404 (61.61%) unigenes were matched to known proteins in the NCBI non-redundant (Nr) protein database. These unigenes were further functionally annotated with gene ontology (GO), cluster of orthologous groups of proteins (COG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. A large number of genes potentially involved in insecticide resistance were manually curated, including 68 putative cytochrome P450 genes, 37 putative glutathione S-transferase (GST) genes, 19 putative carboxyl/cholinesterase (CCE) genes, and other 126 transcripts to contain target site sequences or encoding detoxification genes representing eight types of resistance enzymes. Furthermore, to gain insight into the molecular basis of the L. entomophila toward thermal stresses, 25 heat shock protein (Hsp) genes were identified. In addition, 1,100 SSRs and 57,757 SNPs were detected and 231 pairs of SSR primes were designed for investigating the genetic diversity in future.

Conclusions/Significance

We developed a comprehensive transcriptomic database for L. entomophila. These sequences and putative molecular markers would further promote our understanding of the molecular mechanisms underlying insecticide resistance or environmental stress, and will facilitate studies on population genetics for psocids, as well as providing useful information for functional genomic research in the future.

Improved short-sequence-repeat genotyping of Mycobacterium avium subsp. paratuberculosis by using matrix-assisted laser desorption ionization-time of flight mass spectrometry

Accurate sequence analysis of mononucleotide repeat regions is difficult, complicating the use of short sequence repeats (SSRs) as a tool for bacterial strain discrimination. Although multiple SSR loci in the genome of Mycobacterium avium subsp. paratuberculosis allow genotyping of M. avium subsp. paratuberculosis isolates with high discriminatory power, further characterization of the most discriminatory loci is limited due to inherent difficulties in sequencing mononucleotide repeats. Here, a method was evaluated using matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) as an alternative to Sanger sequencing to further differentiate the dominant mycobacterial interspersed repetitive-unit (MIRU)-variable-number tandem-repeat (VNTR) M. avium subsp. paratuberculosis type (n = 37) in Canadian dairy herds by targeting a highly discriminatory mononucleotide SSR locus. First, PCR-amplified DNA was digested with two restriction enzymes to yield a sufficiently small fragment containing the SSR locus. Second, MALDI-TOF MS was performed to identify the mass, and thus repeat length, of the target. Sufficiently intense, discriminating spectra were obtained to determine repeat lengths up to 15, an improvement over the limit of 11 using traditional sequencing techniques. Comparison to synthetic oligonucleotides and Sanger sequencing results confirmed a valid and reproducible assay that increased discrimination of the dominant M. avium subsp. paratuberculosis MIRU-VNTR type. Thus, MALDI-TOF MS was a reliable, fast, and automatable technique to accurately resolve M. avium subsp. paratuberculosis genotypes based on SSRs.