Microsatellites, from molecules to populations and back

Assessment of population genetic structure in the arbovirus vector midge, Culicoides brevitarsis (Diptera: Ceratopogonidae), using multi-locus DNA microsatellites

Bluetongue virus (BTV) is a major pathogen of ruminants that is transmitted by biting midges (Culicoides spp.). Australian BTV serotypes have origins in Asia and are distributed across the continent into two distinct episystems, one in the north and another in the east. Culicoides brevitarsis is the major vector of BTV in Australia and is distributed across the entire geographic range of the virus. Here, we describe the isolation and use of DNA microsatellites and gauge their ability to determine population genetic connectivity of C. brevitarsis within Australia and with countries to the north. Eleven DNA microsatellite markers were isolated using a novel genomic enrichment method and identified as useful for genetic analyses of sampled populations in Australia, northern Papua New Guinea (PNG) and Timor-Leste. Significant (P < 0.05) population genetic subdivision was observed between all paired regions, though the highest levels of genetic sub-division involved pair-wise tests with PNG (PNG vs. Australia (FST = 0.120) and PNG vs. Timor-Leste (FST = 0.095)). Analysis of multi-locus allelic distributions using STRUCTURE identified a most probable two-cluster population model, which separated PNG specimens from a cluster containing specimens from Timor-Leste and Australia. The source of incursions of this species in Australia is more likely to be Timor-Leste than PNG. Future incursions of BTV positive C. brevitarsis into Australia may be genetically identified to their source populations using these microsatellite loci. The vector's panmictic genetic structure within Australia cannot explain the differential geographic distribution of BTV serotypes.

Assessment of population genetic structure in the arbovirus vector midge, Culicoides brevitarsis (Diptera: Ceratopogonidae), using multi-locus DNA microsatellites

Bluetongue virus (BTV) is a major pathogen of ruminants that is transmitted by biting midges (Culicoides spp.). Australian BTV serotypes have origins in Asia and are distributed across the continent into two distinct episystems, one in the north and another in the east. Culicoides brevitarsis is the major vector of BTV in Australia and is distributed across the entire geographic range of the virus. Here, we describe the isolation and use of DNA microsatellites and gauge their ability to determine population genetic connectivity of C. brevitarsis within Australia and with countries to the north. Eleven DNA microsatellite markers were isolated using a novel genomic enrichment method and identified as useful for genetic analyses of sampled populations in Australia, northern Papua New Guinea (PNG) and Timor-Leste. Significant (P < 0.05) population genetic subdivision was observed between all paired regions, though the highest levels of genetic sub-division involved pair-wise tests with PNG (PNG vs. Australia (F_ST = 0.120) and PNG vs. Timor-Leste (F_ST = 0.095)). Analysis of multi-locus allelic distributions using STRUCTURE identified a most probable two-cluster population model, which separated PNG specimens from a cluster containing specimens from Timor-Leste and Australia. The source of incursions of this species in Australia is more likely to be Timor-Leste than PNG. Future incursions of BTV positive C. brevitarsis into Australia may be genetically identified to their source populations using these microsatellite loci. The vector’s panmictic genetic structure within Australia cannot explain the differential geographic distribution of BTV serotypes.

Population structure of the Indonesian giant tiger shrimp Penaeus monodon: a window into evolutionary similarities between paralogous mitochondrial DNA sequences and their genomes

Here we used both microsatellites and mtCR (mitochondrial DNA control region) sequences as genetic markers to examine the genetic diversity and population structure of Penaeus monodon shrimp from six Indonesian regions. The microsatellite data showed that shrimp from the Indian and the Pacific Ocean were genetically distinct from each other. It has been reported previously that P. monodon mtCR sequences from the Indo-Pacific group into two major paralogous clades of unclear origin. Here we show that the population structure inferred from mtCR sequences matches the microsatellite-based population structure for one of these clades. This is consistent with the notion that this mtCR clade shares evolutionary history with nuclear DNA and may thus represent nuclear mitochondrial pseudogenes (Numts).

Genome-Wide Survey and Analysis of Microsatellite Sequences in Bovid Species

Microsatellites or simple sequence repeats (SSRs) have become the most popular source of genetic markers, which are ubiquitously distributed in many eukaryotic and prokaryotic genomes. This is the first study examining and comparing SSRs in completely sequenced genomes of the Bovidae. We analyzed and compared the number of SSRs, relative abundance, relative density, guanine-cytosine (GC) content and proportion of SSRs in six taxonomically different bovid species: Bos taurus, Bubalus bubalis, Bos mutus, Ovis aries, Capra hircus, and Pantholops hodgsonii. Our analysis revealed that, based on our search criteria, the total number of perfect SSRs found ranged from 663,079 to 806,907 and covered from 0.44% to 0.48% of the bovid genomes. Relative abundance and density of SSRs in these Bovinae genomes were non-significantly correlated with genome size (Pearson, r < 0.420, p > 0.05). Perfect mononucleotide SSRs were the most abundant, followed by the pattern: perfect di- > tri- > penta- > tetra- > hexanucleotide SSRs. Generally, the number of SSRs, relative abundance, and relative density of SSRs decreased as the motif repeat length increased in each species of Bovidae. The most GC-content was in trinucleotide SSRs and the least was in the mononucleotide SSRs in the six bovid genomes. The GC-contents of tri- and pentanucleotide SSRs showed a great deal of similarity among different chromosomes of B. taurus, O. aries, and C. hircus. SSR number of all chromosomes in the B. taurus, O.aries, and C. hircus is closely positively correlated with chromosome sequence size (Pearson, r > 0.980, p < 0.01) and significantly negatively correlated with GC-content (Pearson, r < -0.638, p < 0.01). Relative abundance and density of SSRs in all chromosomes of the three species were significantly negatively correlated with GC-content (Pearson, r < -0.333, P < 0.05) but not significantly correlated with chromosome sequence size (Pearson, r < -0.185, P > 0.05). Relative abundances of the same nucleotide SSR type showed great similarity among different chromosomes of B. taurus, O. aries, and C. hircus.

De Novo Assembly of Bitter Gourd Transcriptomes: Gene Expression and Sequence Variations in Gynoecious and Monoecious Lines

Bitter gourd (Momordica charantia L.) is a nutritious vegetable crop of Asian origin, used as a medicinal herb in Indian and Chinese traditional medicine. Molecular breeding in bitter gourd is in its infancy, due to limited molecular resources, particularly on functional markers for traits such as gynoecy. We performed de novo transcriptome sequencing of bitter gourd using Illumina next-generation sequencer, from root, flower buds, stem and leaf samples of gynoecious line (Gy323) and a monoecious line (DRAR1). A total of 65,540 transcripts for Gy323 and 61,490 for DRAR1 were obtained. Comparisons revealed SNP and SSR variations between these lines and, identification of gene classes. Based on available transcripts we identified 80 WRKY transcription factors, several reported in responses to biotic and abiotic stresses; 56 ARF genes which play a pivotal role in auxin-regulated gene expression and development. The data presented will be useful in both functions studies and breeding programs in bitter gourd.

Multilocus microsatellite typing of Leishmania and clinical applications: a review

Microsatellite markers have been used for Leishmania genetic studies worldwide, giving useful insight into leishmaniasis epidemiology. Understanding the geographic distribution, dynamics of Leishmania populations, and disease epidemiology improved markedly with this tool. In endemic foci, the origins of antimony-resistant strains and multidrug treatment failures were explored with multilocus microsatellite typing (MLMT). High genetic variability was detected but no association between parasite genotypes and drug resistance was established. An association between MLMT profiles and clinical disease manifestations was highlighted in only three studies and this data needs further confirmation. At the individual level, MLMT provided information on relapse and reinfection when multiple leishmaniasis episodes occurred. This information could improve knowledge of epidemiology and guide therapeutic choices for active chronic visceral leishmaniasis, the disease form in some HIV-positive patients.

Admixture of Eastern and Western European Red Deer Lineages as a Result of Postglacial Recolonization of the Czech Republic (Central Europe)

Due to a restriction of the distributional range of European red deer (Cervus elaphus L.) during the Quaternary and subsequent recolonization of Europe from different refugia, a clear phylogeographical pattern in genetic structure has been revealed using mitochondrial DNA markers. In Central Europe, 2 distinct, eastern and western, lineages of European red deer are present; however, admixture between them has not yet been studied in detail. We used mitochondrial DNA (control region and cytochrome b gene) sequences and 22 microsatellite loci from 522 individuals to investigate the genetic diversity of red deer in what might be expected to be an intermediate zone. We discovered a high number of unique mtDNA haplotypes belonging to each lineage and high levels of genetic diversity (cyt b H = 0.867, D-loop H = 0.914). The same structuring of red deer populations was also revealed by microsatellite analysis, with results from both analyses thus suggesting a suture zone between the 2 lineages. Despite the fact that postglacial recolonization of Central Europe by red deer occurred more than 10000 years ago, the degree of admixture between the 2 lineages is relatively small, with only 10.8% admixed individuals detected. Direct translocations of animals by humans have slightly blurred the pattern in this region; however, this blurring was more apparent when using maternally inherited markers than nuclear markers.

Extremely low microsatellite diversity but distinct population structure in a long-lived threatened species, the Australian lungfish Neoceratodus forsteri (Dipnoi)

The Australian lungfish is a unique living representative of an ancient dipnoan lineage, listed as ‘vulnerable’ to extinction under Australia’s Environment Protection and Biodiversity Conservation Act 1999. Historical accounts indicate this species occurred naturally in two adjacent river systems in Australia, the Burnett and Mary. Current day populations in other rivers are thought to have arisen by translocation from these source populations. Early genetic work detected very little variation and so had limited power to answer questions relevant for management including how genetic variation is partitioned within and among sub-populations. In this study, we use newly developed microsatellite markers to examine samples from the Burnett and Mary Rivers, as well as from two populations thought to be of translocated origin, Brisbane and North Pine. We test whether there is significant genetic structure among and within river drainages; assign putatively translocated populations to potential source populations; and estimate effective population sizes. Eleven polymorphic microsatellite loci genotyped in 218 individuals gave an average within-population heterozygosity of 0.39 which is low relative to other threatened taxa and for freshwater fishes in general. Based on F_ST values (average over loci = 0.11) and STRUCTURE analyses, we identify three distinct populations in the natural range, one in the Burnett and two distinct populations in the Mary. These analyses also support the hypothesis that the Mary River is the likely source of translocated populations in the Brisbane and North Pine rivers, which agrees with historical published records of a translocation event giving rise to these populations. We were unable to obtain bounded estimates of effective population size, as we have too few genotype combinations, although point estimates were low, ranging from 29 - 129. We recommend that, in order to preserve any local adaptation in the three distinct populations that they be managed separately.

Accurate typing of short tandem repeats from genome-wide sequencing data and its applications

Short tandem repeats (STRs) are implicated in dozens of human genetic diseases and contribute significantly to genome variation and instability. Yet profiling STRs from short-read sequencing data is challenging because of their high sequencing error rates. Here, we developed STR-FM, short tandem repeat profiling using flank-based mapping, a computational pipeline that can detect the full spectrum of STR alleles from short-read data, can adapt to emerging read-mapping algorithms, and can be applied to heterogeneous genetic samples (e.g., tumors, viruses, and genomes of organelles). We used STR-FM to study STR error rates and patterns in publicly available human and in-house generated ultradeep plasmid sequencing data sets. We discovered that STRs sequenced with a PCR-free protocol have up to ninefold fewer errors than those sequenced with a PCR-containing protocol. We constructed an error correction model for genotyping STRs that can distinguish heterozygous alleles containing STRs with consecutive repeat numbers. Applying our model and pipeline to Illumina sequencing data with 100-bp reads, we could confidently genotype several disease-related long trinucleotide STRs. Utilizing this pipeline, for the first time we determined the genome-wide STR germline mutation rate from a deeply sequenced human pedigree. Additionally, we built a tool that recommends minimal sequencing depth for accurate STR genotyping, depending on repeat length and sequencing read length. The required read depth increases with STR length and is lower for a PCR-free protocol. This suite of tools addresses the pressing challenges surrounding STR genotyping, and thus is of wide interest to researchers investigating disease-related STRs and STR evolution.

Characterisation of a novel panel of polymorphic microsatellite loci for the liver fluke, Fasciola hepatica, using a next generation sequencing approach

•

2448 microsatellite loci were identified within 83 Mb of F. hepatica genome sequence data.

•

A panel of 15 polymorphic loci were developed and validated using genomic DNA from 46 parasites.

•

The panel was developed and optimised as a multiplex PCR protocol.

•

All loci could be amplified from several F. hepatica lifecycle stages with the multiplex approach.

The liver fluke, Fasciola hepatica is an economically important pathogen of sheep and cattle and has been described by the WHO as a re-emerging zoonosis. Control is heavily reliant on the use of drugs, particularly triclabendazole and as a result resistance has now emerged. The population structure of F. hepatica is not well known, yet it can impact on host–parasite interactions and parasite control with drugs, particularly regarding the spread of triclabendazole resistance. We have identified 2448 potential microsatellites from 83 Mb of F. hepatica genome sequence using msatfinder. Thirty-five loci were developed and optimised for microsatellite PCR, resulting in a panel of 15 polymorphic loci, with a range of three to 15 alleles. This panel was validated on genomic DNA from 46 adult F. hepatica; 38 liver flukes sourced from a Northwest abattoir, UK and 8 liver flukes from an established isolate (Shrewsbury; Ridgeway Research). Evidence for null alleles was found at four loci (Fh_1, Fh_8, Fh_13 and Fh_14), which showed markedly higher levels of homozygosity than the remaining 11 loci. Of the 38 liver flukes isolated from cattle livers (n = 10) at the abattoir, 37 genotypes were identified. Using a multiplex approach all 15 loci could be amplified from several life cycle stages that typically yield low amounts of DNA, including metacercariae, the infective life cycle stage present on pasture, highlighting the utility of this multiplex microsatellite panel. This study reports the largest panel of microsatellite markers available to date for population studies of F. hepatica and the first multiplex panel of microsatellite markers that can be used for several life cycle stages.

Post-fragmentation population structure in a cooperative breeding Afrotropical cloud forest bird: emergence of a source-sink population network

The impact of demographic parameters on the genetic population structure and viability of organisms is a long-standing issue in the study of fragmented populations. Demographic and genetic tools are now readily available to estimate census and effective population sizes and migration and gene flow rates with increasing precision. Here we analysed the demography and genetic population structure over a recent 15-year time span in five remnant populations of Cabanis's greenbul (Phyllastrephus cabanisi), a cooperative breeding bird in a severely fragmented cloud forest habitat. Contrary to our expectation, genetic admixture and effective population sizes slightly increased, rather than decreased between our two sampling periods. In spite of small effective population sizes in tiny forest remnants, none of the populations showed evidence of a recent population bottleneck. Approximate Bayesian modelling, however, suggested that differentiation of the populations coincided at least partially with an episode of habitat fragmentation. The ratio of meta-Ne to meta-Nc was relatively low for birds, which is expected for cooperative breeding species, while Ne /Nc ratios strongly varied among local populations. While the overall trend of increasing population sizes and genetic admixture may suggest that Cabanis's greenbuls increasingly cope with fragmentation, the time period over which these trends were documented is rather short relative to the average longevity of tropical species. Furthermore, the critically low Nc in the small forest remnants keep the species prone to demographic and environmental stochasticity, and it remains open if, and to what extent, its cooperative breeding behaviour helps to buffer such effects.

Heterozygosity-fitness correlations in a declining seabird population

Loss of genetic diversity is thought to lead to increased risk of extinction in endangered populations due to decreasing fitness of homozygous individuals. Here, we evaluated the presence of inbreeding depression in a long-lived seabird, the European shag (Phalacrocorax aristotelis), after a severe decline in population size by nearly 70%. During three reproductive seasons, 85 breeders were captured and genotyped at seven microsatellite loci. Nest sites were monitored during the breeding season to estimate reproductive success as the number of chicks surviving to full-size-grown per nest. Captured birds were tagged with a ring with an individual code, and resighting data were collected during 7-year period. We found a strong effect of multilocus heterozygosity on female reproductive performance, and a significant, although weaker, effect on breeder survival. However, our matrix population model suggests that this relatively small effect of genetic diversity on breeder survival may have a profound effect on fitness. This highlights the importance of integrating life history consequences in HFC studies. Importantly, heterozygosity was correlated across loci, suggesting that genomewide effects, rather than single loci, are responsible for the observed HFCs. Overall, the HFCs are a worrying symptom of genetic erosion in this declining population. Many long-lived species are prone to extinction, and future studies should evaluate the magnitude of fitness impact of genetic deterioration on key population parameters, such as survival of breeders.

Development and characterization of microsatellite loci for the pseudometallophyte Commelina communis (Commelinaceae)

•

PREMISE OF THE STUDY

Microsatellite primers were developed for the pseudometallophyte Commelina communis (Commelinaceae), an important pioneer plant for phytoremediation of copper-contaminated soil. Two wild populations collected from metalliferous and nonmetalliferous sites were used to assess the polymorphism at each locus. •

METHODS AND RESULTS

Based on the Fast Isolation by AFLP of Sequences COntaining repeats (FIASCO) method, a total of 34 pairs of simple sequence repeat (SSR) markers were designed. When 40 specimens from two populations were screened, 12 microsatellite loci were found to be highly polymorphic. The number of alleles per locus ranged from one to 11 and the observed and expected heterozygosity per locus ranged from 0.000 to 1.000 and from 0.195 to 0.941, respectively. •

CONCLUSIONS

These markers will be useful for examining genetic diversity, population structure, and gene flow in populations of C. communis under different edaphic conditions and guiding sustainable management plans for phytoremediation.

Genetic diversity in aspen and its relation to arthropod abundance

The ecological consequences of biodiversity have become a prominent public issue. Little is known on the effect of genetic diversity on ecosystem services. Here, a diversity experiment was established with European and North American aspen (Populus tremula, P. tremuloides) planted in plots representing either a single deme only or combinations of two, four and eight demes. The goals of this study were to explore the complex inter- and intraspecific genetic diversity of aspen and to then relate three measures for diversity (deme diversity, genetic diversity determined as Shannon index or as expected heterozygosity) to arthropod abundance. Microsatellite and AFLP markers were used to analyze the genetic variation patterns within and between the aspen demes and deme mixtures. Large differences were observed regarding the genetic diversity within demes. An analysis of molecular variance revealed that most of the total genetic diversity was found within demes, but the genetic differentiation among demes was also high. The complex patterns of genetic diversity and differentiation resulted in large differences of the genetic variation within plots. The average diversity increased from plots with only one deme to plots with two, four, and eight demes, respectively and separated plots with and without American aspen. To test whether intra- and interspecific diversity impacts on ecosystem services, arthropod abundance was determined. Increasing genetic diversity of aspen was related to increasing abundance of arthropods. However, the relationship was mainly driven by the presence of American aspen suggesting that species identity overrode the effect of intraspecific variation of European aspen.

Skim-based genotyping by sequencing

Genotyping by sequencing (GBS) is a relatively new method used to determine the differences in the genetic makeup of individuals. Its novelty stems from a combination of two already available methods: genotyping and next-generation sequencing. Depending on the individual study design GBS protocols can take multiple forms, however most share a sequence of core steps that have to be undertaken. These include: sequencing of the DNA from the individuals of interest (usually two parents of a mapping population and their progeny), mapping of the sequencing reads to the reference sequence, SNP calling and filtering, SNP genotyping and imputation, followed by haplotype identification and downstream analysis. GBS has a range of applications from general marker discovery, haplotype identification, and recombination characterization to quantitative trait locus (QTL) analysis, genome-wide association studies (GWAS), and genomic selection (GS). It has already been applied to a range of plant species including: rice, maize, artichoke, and Arabidopsis thaliana. It is a promising approach which is likely to provide new and important insights into plant biology.

Development and validation of SSR markers for the plant-parasitic nematode Subanguina moxae using genome assembly of Illumina pair-end reads

Subanguina moxae, belonging to the subfamily Anguininae, is an obligate parasite of Artemisia plants, which are widely used as cooking herbs and in traditional medicine in East Asia. Because the nematode is distributed throughout East Russia and East Asia, there is concern about the potential for significant damage to commercial farming; however, details about its biology remain unclear. To investigate the genetic diversity of S. moxae, we developed 2243 simple sequence repeat (SSR) markers using Illumina short reads of the genomic DNA. We validated 100 randomly selected markers indicating their robustness and examined polymorphisms among nematode populations sampled from four different locations in Japan. These SSR markers will be a useful tool for understanding the population structure and transmission patterns of this parasitic nematode.

Specific diagnosis of Opisthorchis viverrini using loop-mediated isothermal amplification (LAMP) targeting parasite microsatellites

Opisthorchis viverrini and other food-borne trematode infections are major health problems in Thailand, the Lao People's Democratic Republic, Vietnam and Cambodia. Differential diagnosis of O. viverrini based on the microscopic observation of parasite eggs is difficult in areas where Clonorchis sinensis and minute intestinal flukes coexist. Recently, loop-mediated isothermal amplification (LAMP) has been widely used for detection and identification of trematode for its simple method that is useful in low-resource or field settings. We have reported ITS1-LAMP assay to detect O. viverrini infection from human feces. The sensitivity and specificity of the test was 100% and 61.5%. The sensitivity of the test appeared to be higher than microscopic egg examination; however non-specific amplification from other parasites could not be ruled out. We therefore targeted microsatellites of O. viverrini that is a species specific sequence. By using hydroxyl naphthol blue (HNB)-LAMP, O. viverrini microsatellite 6 (OVMS6) could specifically amplify DNA from O. viverrini genome, but not other parasites such as C. sinensis, Opisthorchis felineus, Centrocestus caninus, Haplorchis taichui, Fasciola gigantica and Haplorchoodes sp. The detection limit of the test is 1 ng genomic DNA, which was 1000 times lower than the ITS1-LAMP, but targeting microstellites showed more specific detection of O. viverrini. In addition, the colorimetric LAMP assay was simple and effective; this makes it potentially applicable for point-of-care diagnosis.

Exome-wide somatic microsatellite variation is altered in cells with DNA repair deficiencies

Microsatellites (MST), tandem repeats of 1–6 nucleotide motifs, are mutational hot-spots with a bias for insertions and deletions (INDELs) rather than single nucleotide polymorphisms (SNPs). The majority of MST instability studies are limited to a small number of loci, the Bethesda markers, which are only informative for a subset of colorectal cancers. In this paper we evaluate non-haplotype alleles present within next-gen sequencing data to evaluate somatic MST variation (SMV) within DNA repair proficient and DNA repair defective cell lines. We confirm that alleles present within next-gen data that do not contribute to the haplotype can be reliably quantified and utilized to evaluate the SMV without requiring comparisons of matched samples. We observed that SMV patterns found in DNA repair proficient cell lines without DNA repair defects, MCF10A, HEK293 and PD20 RV:D2, had consistent patterns among samples. Further, we were able to confirm that changes in SMV patterns in cell lines lacking functional BRCA2, FANCD2 and mismatch repair were consistent with the different pathways perturbed. Using this new exome sequencing analysis approach we show that DNA instability can be identified in a sample and that patterns of instability vary depending on the impaired DNA repair mechanism, and that genes harboring minor alleles are strongly associated with cancer pathways. The MST Minor Allele Caller used for this study is available at https://github.com/zalmanv/MST_minor_allele_caller.

Conservation genetics of a threatened butterfly: comparison of allozymes, RAPDs and microsatellites

Background

Addressing genetic issues in the management of fragmented wild populations of threatened species is one of the most important challenges in conservation biology. Nowadays, a diverse array of molecular methods exists to assess genetic diversity and differentiation of wild populations such as allozymes, dominant markers and co-dominant markers. However it remains worthwhile i) to compare the genetic estimates obtained using those several markers in order to ii) test their relative utility, reliability and relevance and iii) the impact of these results for the design of species-specific conservation measures.

Results

Following the successful isolation of 15 microsatellites loci for the cranberry fritillary butterfly, Boloria aquilonaris, we analyzed the genetic diversity and structure of eight populations located in four different landscapes, at both the regional and the landscape scales. We confront results based on microsatellites to those obtained using allozymes and RAPDs on the same samples. Genetic population analyses using different molecular markers indicate that the B. aquilonaris populations are characterized by a weak genetic variation, likely due to low effective population size and low dispersal at the regional scale. This results in inbreeding in some populations, which may have detrimental consequences on their long term viability. However, gene flow within landscape is limited but not inexistent, with some long range movements resulting in low or no isolation by distance. Spatial structuring was detected among the most isolated populations.

Conclusions

The use of allozymes and RAPD are of very limited value to determine population structuring at small spatial (i.e. landscape) scales, microsatellites giving much higher estimate resolution. The use of RAPD data is also limited for evidencing inbreeding. However, coarse-grain spatial structure (i.e. regional scale), and gene flow estimates based on RAPD and microsatellites data gave congruent results. At a time with increasing development of new molecular methods and markers, dominant markers may still be worthwhile to consider in organisms for which no genomic information is available, and for which limited resources are available.

Genetics reveal the origin and timing of a cryptic insular introduction of muskrats in North America

The muskrat, Ondatra zibethicus, is a semiaquatic rodent native to North America that has become a highly successful invader across Europe, Asia, and South America. It can inflict ecological and economic damage on wetland systems outside of its native range. Anecdotal evidence suggests that, in the early 1900s, a population of muskrats was introduced to the Isles of Shoals archipelago, located within the Gulf of Maine, for the purposes of fur harvest. However, because muskrats are native to the northeastern coast of North America, their presence on the Isles of Shoals could be interpreted as part of the native range of the species, potentially obscuring management planning and biogeographic inferences. To investigate their introduced status and identify a historic source population, muskrats from Appledore Island of the Isles of Shoals, and from the adjacent mainland of Maine and New Hampshire, were compared for mitochondrial cytochrome b sequences and allele frequencies at eight microsatellite loci. Appledore Island muskrats consistently exhibited reduced genetic diversity compared with mainland populations, and displayed signatures of a historic bottleneck. The distribution of mitochondrial haplotypes is suggestive of a New Hampshire source population. The data presented here are consistent with a human-mediated introduction that took place in the early 1900s. This scenario is further supported by the zooarchaeological record and island biogeographic patterns. This is the first genetic study of an introduced muskrat population within US borders and of any island muskrat population, and provides an important contrast with other studies of introduced muskrat populations worldwide.

'True' null allele detection in microsatellite loci: a comparison of methods, assessment of difficulties and survey of possible improvements

Null alleles are alleles that for various reasons fail to amplify in a PCR assay. The presence of null alleles in microsatellite data is known to bias the genetic parameter estimates. Thus, efficient detection of null alleles is crucial, but the methods available for indirect null allele detection return inconsistent results. Here, our aim was to compare different methods for null allele detection, to explain their respective performance and to provide improvements. We applied several approaches to identify the 'true' null alleles based on the predictions made by five different methods, used either individually or in combination. First, we introduced simulated 'true' null alleles into 240 population data sets and applied the methods to measure their success in detecting the simulated null alleles. The single best-performing method was ML-NullFreq_frequency. Furthermore, we applied different noise reduction approaches to improve the results. For instance, by combining the results of several methods, we obtained more reliable results than using a single one. Rule-based classification was applied to identify population properties linked to the false discovery rate. Rules obtained from the classifier described which population genetic estimates and loci characteristics were linked to the success of each method. We have shown that by simulating 'true' null alleles into a population data set, we may define a null allele frequency threshold, related to a desired true or false discovery rate. Moreover, using such simulated data sets, the expected null allele homozygote frequency may be estimated independently of the equilibrium state of the population.

Conserved microsatellites in ants enable population genetic and colony pedigree studies across a wide range of species

Broadly applicable polymorphic genetic markers are essential tools for population genetics, and different types of markers have been developed for this purpose. Microsatellites have been employed as particularly polymorphic markers for over 20 years. However, PCR primers for microsatellite loci are often not useful outside the species for which they were designed. This implies that a new set of loci has to be identified and primers developed for every new study species. To overcome this constraint, we identified 45 conserved microsatellite loci based on the eight currently available ant genomes and designed primers for PCR amplification. Among these loci, we chose 24 for in-depth study in six species covering six different ant subfamilies. On average, 11.16 of these 24 loci were polymorphic and in Hardy-Weinberg equilibrium in any given species. The average number of alleles for these polymorphic loci within single populations of the different species was 4.59. This set of genetic markers will thus be useful for population genetic and colony pedigree studies across a wide range of ant species, supplementing the markers available for previously studied species and greatly facilitating the study of the many ant species lacking genetic markers. Our study shows that it is possible to develop microsatellite loci that are both conserved over a broad range of taxa, yet polymorphic within species. This should encourage researchers to develop similar tools for other large taxonomic groups.

Evolution and origin of sympatric shallow-water morphotypes of Lake Trout, Salvelinus namaycush, in Canada's Great Bear Lake

Range expansion in north-temperate fishes subsequent to the retreat of the Wisconsinan glaciers has resulted in the rapid colonization of previously unexploited, heterogeneous habitats and, in many situations, secondary contact among conspecific lineages that were once previously isolated. Such ecological opportunity coupled with reduced competition likely promoted morphological and genetic differentiation within and among post-glacial fish populations. Discrete morphological forms existing in sympatry, for example, have now been described in many species, yet few studies have directly assessed the association between morphological and genetic variation. Morphotypes of Lake Trout, Salvelinus namaycush, are found in several large-lake systems including Great Bear Lake (GBL), Northwest Territories, Canada, where several shallow-water forms are known. Here, we assess microsatellite and mitochondrial DNA variation among four morphotypes of Lake Trout from the five distinct arms of GBL, and also from locations outside of this system to evaluate several hypotheses concerning the evolution of morphological variation in this species. Our data indicate that morphotypes of Lake Trout from GBL are genetically differentiated from one another, yet the morphotypes are still genetically more similar to one another compared with populations from outside of this system. Furthermore, our data suggest that Lake Trout colonized GBL following dispersal from a single glacial refugium (the Mississippian) and support an intra-lake model of divergence. Overall, our study provides insights into the origins of morphological and genetic variation in post-glacial populations of fishes and provides benchmarks important for monitoring Lake Trout biodiversity in a region thought to be disproportionately susceptible to impacts from climate change.

Isolation and characterization of 16 polymorphic microsatellite loci for Euchiloglanis kishinouyei

A total of 16 novel polymorphic microsatellite marker loci were isolated from a genomic library of Euchiloglanis kishinouyei and further characterized using a sample from a wild population consisting of 40 individuals. The number of alleles among loci ranged from three to 17 and the observed and expected heterozygosities ranged from 0·030 to 0·950 and 0·163 to 0·980, respectively. The average polymorphic information content (PIC) of all loci was 0·450.

Discriminating micropathogen lineages and their reticulate evolution through graph theory-based network analysis: the case of Trypanosoma cruzi, the agent of Chagas disease

Micropathogens (viruses, bacteria, fungi, parasitic protozoa) share a common trait, which is partial clonality, with wide variance in the respective influence of clonality and sexual recombination on the dynamics and evolution of taxa. The discrimination of distinct lineages and the reconstruction of their phylogenetic history are key information to infer their biomedical properties. However, the phylogenetic picture is often clouded by occasional events of recombination across divergent lineages, limiting the relevance of classical phylogenetic analysis and dichotomic trees. We have applied a network analysis based on graph theory to illustrate the relationships among genotypes of Trypanosoma cruzi, the parasitic protozoan responsible for Chagas disease, to identify major lineages and to unravel their past history of divergence and possible recombination events. At the scale of T. cruzi subspecific diversity, graph theory-based networks applied to 22 isoenzyme loci (262 distinct Multi-Locus-Enzyme-Electrophoresis -MLEE) and 19 microsatellite loci (66 Multi-Locus-Genotypes -MLG) fully confirms the high clustering of genotypes into major lineages or "near-clades". The release of the dichotomic constraint associated with phylogenetic reconstruction usually applied to Multilocus data allows identifying putative hybrids and their parental lineages. Reticulate topology suggests a slightly different history for some of the main "near-clades", and a possibly more complex origin for the putative hybrids than hitherto proposed. Finally the sub-network of the near-clade T. cruzi I (28 MLG) shows a clustering subdivision into three differentiated lesser near-clades ("Russian doll pattern"), which confirms the hypothesis recently proposed by other investigators. The present study broadens and clarifies the hypotheses previously obtained from classical markers on the same sets of data, which demonstrates the added value of this approach. This underlines the potential of graph theory-based network analysis for describing the nature and relationships of major pathogens, thereby opening stimulating prospects to unravel the organization, dynamics and history of major micropathogen lineages.

Exploring heterozygosity-survival correlations in a wild songbird population: contrasting effects between juvenile and adult stages

The relationship between genetic diversity and fitness, a major issue in evolutionary and conservation biology, is expected to be stronger in traits affected by many loci and those directly influencing fitness. Here we explore the influence of heterozygosity measured at 15 neutral markers on individual survival, one of the most important parameters determining individual fitness. We followed individual survival up to recruitment and during subsequent adult life of 863 fledgling pied flycatchers born in two consecutive breeding seasons. Mark-recapture analyses showed that individual heterozygosity did not influence juvenile or adult survival. In contrast, the genetic relatedness of parents was negatively associated with the offspring’s survival during the adult life, but this effect was not apparent in the juvenile (from fledgling to recruitment) stage. Stochastic factors experienced during the first year of life in this long-distance migratory species may have swamped a relationship between heterozygosity and survival up to recruitment.

Clonal expansion of the Pseudogymnoascus destructans genotype in North America is accompanied by significant variation in phenotypic expression

Pseudogymnoascus destructans is the causative agent of an emerging infectious disease that threatens populations of several North American bat species. The fungal disease was first observed in 2006 and has since caused the death of nearly six million bats. The disease, commonly known as white-nose syndrome, is characterized by a cutaneous infection with P. destructans causing erosions and ulcers in the skin of nose, ears and/or wings of bats. Previous studies based on sequences from eight loci have found that isolates of P. destructans from bats in the US all belong to one multilocus genotype. Using the same multilocus sequence typing method, we found that isolates from eastern and central Canada also had the same genotype as those from the US, consistent with the clonal expansion of P. destructans into Canada. However, our PCR fingerprinting revealed that among the 112 North American isolates we analyzed, three, all from Canada, showed minor genetic variation. Furthermore, we found significant variations among isolates in mycelial growth rate; the production of mycelial exudates; and pigment production and diffusion into agar media. These phenotypic differences were influenced by culture medium and incubation temperature, indicating significant variation in environmental condition - dependent phenotypic expression among isolates of the clonal P. destructans genotype in North America.

Molecular markers for genetic diversity, gene flow and genetic population structure of freshwater mussel species

Freshwater mussel species are in global decline. Anthropogenic changes of river channels and the decrease of autochthonous fish population, the natural hosts of mussels larval stages (glochidia), are the main causes. Therefore, the conservation of mussel species depends not only on habitat conservation, but also on the availability of the fish host. In Portugal, information concerning most of the mussel species is remarkably scarce. One of the most known species, Unio pictorum is also in decline however, in the basins of the rivers Tua and Sabor (Northeast of Portugal), there is some indication of relatively large populations. The aforementioned rivers can be extremely important for this species conservation not only in Portugal, but also in the remaining Iberian Peninsula. Thus, it is important to obtain data concerning Unio pictorum bioecology (distribution, habitat requirements, population structure, genetic variability, reproductive cycle and recruitment rates), as well as the genetic variability and structure of the population. Concomitantly, information concerning fish population structure, the importance of the different fish species as “glochidia” hosts and their appropriate density to allow effective mussel recruitment, will also be assessed. The achieved data is crucial to obtain information to develop effective management measures in order to promote the conservation of this bivalve species, the conservation of autochthonous fish populations, and consequently the integrity of the river habitats.

Genetic variation in horizontally transmitted fungal endophytes of pine needles reveals population structure in cryptic species

•

PREMISE OF THE STUDY

Fungal endophytes comprise one of the most ubiquitous groups of plant symbionts, inhabiting healthy leaves and stems of all major lineages of plants. Together, they comprise immense species richness, but little is known about the fundamental processes that generate their diversity. Exploration of their population structure is needed, especially with regard to geographic distributions and host affiliations.•

METHODS

We take a multilocus approach to examine genetic variation within and among populations of Lophodermium australe, an endophytic fungus commonly associated with healthy foliage of pines in the southeastern United States. Sampling focused on two pine species ranging from montane to coastal regions of North Carolina and Virginia.•

KEY RESULTS

Our sampling revealed two genetically distinct groups within Lophodermium australe. Our analysis detected less than one migrant per generation between them, indicating that they are distinct species. The species comprising the majority of isolates (major species) demonstrated a panmictic structure, whereas the species comprising the minority of isolates (cryptic species) demonstrated isolation by distance. Distantly related pine species hosted the same Lophodermium species, and host species did not influence genetic structure.•

CONCLUSIONS

We present the first evidence for isolation by distance in a foliar fungal endophyte that is horizontally transmitted. Cryptic species may be common among microbial symbionts and are important to delimit when exploring their genetic structure and microevolutionary processes. The hyperdiversity of endophytic fungi may be explained in part by cryptic species without apparent ecological and morphological differences as well as genetic diversification within rare fungal species across large spatial scales.

Single strand conformation polymorphism of genomic and EST-SSRs marker and its utility in genetic evaluation of sugarcane

Sugarcane is an important crop producing around 75 % of sugar in world and used as first generation biofuel. In present study, the genomic and gene based microsatellite markers were analyzed by low cost Single Strand Confirmation Polymorphism technique for genetic evaluation of 22 selected sugarcane genotypes. Total 16 genomic and 12 Expression Sequence Tag derived markers were able to amplify the selected sugarcane genotypes. Total 138 alleles were amplified of which 99 alleles (72 %) found polymorphic with an average of 4.9 alleles per locus. Microsatellite marker, VCSSR7 and VCSSR 12 showed monomorphic alleles with frequency 7.1 % over the average of 3.5 obtained for polymorphic locus. The level of Polymorphic Information Content (PIC) varied from 0.09 in VCSSR 6 to 0.88 in VCSSR 11 marker respectively with a mean of 0.49. Genomic SSRs showed more polymorphism than EST-SSRs markers on selected sugarcane genotypes whereas, the genetic similarity indices calculated by Jaccard's similarity coefficient varied from 0.55 to 0.81 indicate a high level of genetic similarity among the genotypes that was mainly attributed to intra specific diversity. Hence, the SSR-SSCP technique helped to identify the genetically diverse clones which could be used in crossing program for introgression of sugar and stress related traits in hybrid sugarcane.

Gene flow within and between catchments in the threatened riparian plant Myricaria germanica

One of the major distinctions of riparian habitats is their linearity. In linear habitats, gene flow is predicted to follow a one-dimensional stepping stone model, characterized by bidirectional gene flow between neighboring populations. Here, we studied the genetic structure of Myricaria germanica, a threatened riparian shrub which is capable of both wind and water dispersal. Our data led us to reject the ‘one catchment – one gene pool’ hypothesis as we found support for two gene pools, rather than four as expected in a study area including four catchments. This result also implies that in the history of the studied populations, dispersal across catchments has occurred. Two contemporary catchment-crossing migration events were detected, albeit between spatially proximate catchments. Allelic richness and inbreeding coefficients differed substantially between gene pools. There was significant isolation by distance, and our data confirmed the one-dimensional stepping-stone model of gene flow. Contemporary migration was bidirectional within the studied catchments, implying that dispersal vectors other than water are important for M. germanica.

Isolation of 12 microsatellite markers following a pyrosequencing procedure and cross-priming in two invasive cryptic species, Alexandrium catenella (group IV) and A. tamarense (group III) (Dinophyceae)

Alexandrium catenella (group IV) and Alexandrium tamarense (group III) (Dinophyceae) are two cryptic invasive phytoplankton species belonging to the A. tamarense species complex. Their worldwide spread is favored by the human activities, transportation and climate change. In order to describe their diversity in the Mediterranean Sea and understand their settlements and maintenances in this area, new microsatellite markers were developed based on Thau lagoon (France) samples of A. catenella and A. tamarense strains. In this study twelve new microsatellite markers are proposed. Five of these microsatellite markers show amplifications on A. tamarense and ten on A. catenella. Three of these 12 microsatellite markers allowed amplifications on both cryptic species. Finally, the haplotypic diversity ranged from 0.000 to 0.791 and 0.000 to 0.942 for A. catenella and A. tamarense respectively.

Heterozygosity-fitness correlations in a wild mammal population: accounting for parental and environmental effects

HFCs (heterozygosity–fitness correlations) measure the direct relationship between an individual's genetic diversity and fitness. The effects of parental heterozygosity and the environment on HFCs are currently under-researched. We investigated these in a high-density U.K. population of European badgers (Meles meles), using a multimodel capture–mark–recapture framework and 35 microsatellite loci. We detected interannual variation in first-year, but not adult, survival probability. Adult females had higher annual survival probabilities than adult males. Cubs with more heterozygous fathers had higher first-year survival, but only in wetter summers; there was no relationship with individual or maternal heterozygosity. Moist soil conditions enhance badger food supply (earthworms), improving survival. In dryer years, higher indiscriminate mortality rates appear to mask differential heterozygosity-related survival effects. This paternal interaction was significant in the most supported model; however, the model-averaged estimate had a relative importance of 0.50 and overlapped zero slightly. First-year survival probabilities were not correlated with the inbreeding coefficient (f); however, small sample sizes limited the power to detect inbreeding depression. Correlations between individual heterozygosity and inbreeding were weak, in line with published meta-analyses showing that HFCs tend to be weak. We found support for general rather than local heterozygosity effects on first-year survival probability, and g2 indicated that our markers had power to detect inbreeding. We emphasize the importance of assessing how environmental stressors can influence the magnitude and direction of HFCs and of considering how parental genetic diversity can affect fitness-related traits, which could play an important role in the evolution of mate choice.

Cross-priming of microsatellite loci in subfamily cyprininae (family Cyprinidae): their utility in finding markers for population genetic analysis in three Indian major carps

This study is aimed to identify polymorphic microsatellite markers and establish their potential for population genetics studies in three carp (family cyprinidae; subfamily cyprininae) species, Labeo rohita, Catla catla and Cirrhinus mrigala through use of cyprinid primers. These species have high commercial value and knowledge of genetic variation is important for management of farmed and wild populations. We tested 108 microsatellite primers from 11 species belonging to three different cyprinid subfamilies, Cyprininae, Barbinae and Leuciscinae out of which 63 primers (58.33%) successfully amplified orthologous loci in three focal species. Forty-two loci generated from 29 primers were polymorphic in these three carp species. Sequencing of amplified product confirmed the presence of SSRs in these 42 loci and orthologous nature of the loci. To validate potential of these 42 polymorphic loci in determining the genetic variation, we analyzed 486 samples of three focal species collected from Indus, Ganges and Brahmaputra river systems. Results indicated significant genetic variation, with mean number of alleles per locus ranging from 6.80 to 14.40 and observed heterozygosity ranging from 0.50 to 0.74 in the three focal species. Highly significant (P < 0.00001) allelic homogeneity values revealed that the identified loci can be efficiently used in population genetics analysis of these carp species. Further, thirty-two loci from 19 primers were useful for genotyping in more than one species. The data from the present study was compiled with cross-species amplification data from previous results on eight species of subfamily cyprininae to compare cross-transferability of microsatellite loci. It was revealed that out of 226 heterologous loci amplified, 152 loci that originated from 77 loci exhibited polymorphism and 45 primers were of multispecies utility, common for 2-7 species.

Isolation and characterization of microsatellite markers for Cotinus coggygria Scop. (Anacardiaceae) by 454 pyrosequencing

Cotinus coggygria Scop. (Anacardiaceae) is a deciduous shrub or small tree that is native to a large area covering from southern Europe, east across central Asia, and the Himalayas in northern China. Shotgun 454 pyrosequencing was used to develop microsatellite markers from the genome of C. coggygria. In this study, 349 microsatellite loci were identified from 40,074 individual sequence reads produced by one-sixteenth run, and primer pairs were designed for these loci. To test the primer amplification efficiency, 50 microsatellite primer pairs were tested across 12 individuals from two C. coggygria populations (Wuzhi Mountain: 36°30'N, 113°39'E; Tianlong Mountain: 37°42'N, 112°26'E). Among the 50 tested primer pairs, eight were found to be polymorphic. The average allele number of the microsatellites was 3.5 per locus, with a range from two to five. The inbreeding coefficient ranged from -0.478 to 0.222. The observed and expected heterozygosities varied from 0.167 to 0.750 and from 0.163 to 0.743, respectively. This set of markers is potentially useful for assessing the genetic diversity, as well as for understanding the population structure and phylogeographical and landscape genetic patterns, of C. coggygria.

In silico analysis of SSRs in mitochondrial genomes of fishes

The availability of fish mitochondrial (mt) genomes provides an opportunity to explore the simple sequence repeats. In the present study, mt genomes of 85 fish species reported from Indian subcontinent were downloaded from NCBI and computationally analysed for finding SSRs types, frequency of occurrence, mutation and evolutionary adaptation across species. A total of 92 microsatellites in different nucleotide combinations were detected in 59 species. 26 interspersed SSRs, mostly poly (AT)n were found in the D-loop regions in the species of Cyprinidae. Fifty-six SSRs of 12 bp fixed length were observed in eight genes only. Further, identical repeat motifs were found on the same location in ATP6 and ND4 genes, which were biased towards particular habitat. The comparison of ATP6 and ND4 gene sets to other homologous sequences showed point mutations. This study explores the SSRs discovery and their utility as marker for species and population identification.

Transcriptomic analysis of endangered Chinese salamander: identification of immune, sex and reproduction-related genes and genetic markers

BACKGROUND

The Chinese salamander (Hynobius chinensis), an endangered amphibian species of salamander endemic to China, has attracted much attention because of its value of studying paleontology evolutionary history and decreasing population size. Despite increasing interest in the Hynobius chinensis genome, genomic resources for the species are still very limited. A comprehensive transcriptome of Hynobius chinensis, which will provide a resource for genome annotation, candidate genes identification and molecular marker development should be generated to supplement it.

PRINCIPAL FINDINGS

We performed a de novo assembly of Hynobius chinensis transcriptome by Illumina sequencing. A total of 148,510 nonredundant unigenes with an average length of approximately 580 bp were obtained. In all, 60,388 (40.66%) unigenes showed homologous matches in at least one database and 33,537 (22.58%) unigenes were annotated by all four databases. In total, 41,553 unigenes were categorized into 62 sub-categories by BLAST2GO search, and 19,468 transcripts were assigned to 140 KEGG pathways. A large number of unigenes involved in immune system, local adaptation, reproduction and sex determination were identified, as well as 31,982 simple sequence repeats (SSRs) and 460,923 putative single nucleotide polymorphisms (SNPs).

CONCLUSION

This dataset represents the first transcriptome analysis of the Chinese salamander (Hynobius chinensis), an endangered species, to be also the first time of hynobiidae. The transcriptome will provide valuable resource for further research in discovery of new genes, protection of population, adaptive evolution and survey of various pathways, as well as development of molecule markers in Chinese salamander; and reference information for closely related species.