High-utility conserved avian microsatellite markers enable parentage and population studies across a wide range of species

Development of genetic markers for reproductive management of toucans

Retention of genetic diversity in successive generations is key to successful ex situ programs and will become increasingly important to restore wild populations of threatened animals. When animal genealogy is partly unknown or gaps exist in studbook records, the application of molecular resources facilitates informed breeding. Here, we apply molecular resources to an ex situ breeding population of toucans (Ramphastidae), a bird family zoos commonly maintain. Toucans face population declines from illegal poaching and habitat degradation. We developed novel microsatellite markers using blood samples from 15 Keel-billed Toucans (Ramphastos sulfuratus Lesson 1830). Parentage of two individuals was known a priori, but possible sibship among 13 putative founders-including the parents-was unknown. We compared available avian heterologous and novel microsatellite markers to recover known relationships and reconstruct sibship. Eight of 61 heterologous markers amplified consistently and were polymorphic, but less so than the 18 novel markers. Known sibship (and three sibling pairs whose relatedness was unknown a priori) and paternity-though not maternity except in one case-were well-recovered using both likelihood and pairwise relatedness methods, when incorporating novel but not heterologous markers. Zoo researchers seeking microsatellite primer sets for their breeding toucan populations will likely benefit from our heterologous markers, which can be leveraged both to assess relatedness and select breeding pairs. We recommend that zoo biologists rely on species-specific primers and not optimize heterologous primers for toucan species without molecular resources. We conclude with a brief discussion of modern genotyping methods of interest to zoo researchers.

Genetic diversity and spatial genetic structure support the specialist-generalist variation hypothesis in two sympatric woodpecker species

Species are often arranged along a continuum from “specialists” to “generalists”. Specialists typically use fewer resources, occur in more patchily distributed habitats and have overall smaller population sizes than generalists. Accordingly, the specialist-generalist variation hypothesis (SGVH) proposes that populations of habitat specialists have lower genetic diversity and are genetically more differentiated due to reduced gene flow compared to populations of generalists. Here, expectations of the SGVH were tested by examining genetic diversity, spatial genetic structure and contemporary gene flow in two sympatric woodpecker species differing in habitat specialization. Compared to the generalist great spotted woodpecker (Dendrocopos major), lower genetic diversity was found in the specialist middle spotted woodpecker (Dendrocoptes medius). Evidence for recent bottlenecks was revealed in some populations of the middle spotted woodpecker, but in none of the great spotted woodpecker. Substantial spatial genetic structure and a significant correlation between genetic and geographic distances were found in the middle spotted woodpecker, but only weak spatial genetic structure and no significant correlation between genetic and geographic distances in the great spotted woodpecker. Finally, estimated levels of contemporary gene flow did not differ between the two species. Results are consistent with all but one expectations of the SGVH. This study adds to the relatively few investigations addressing the SGVH in terrestrial vertebrates.

Development of 16 novel EST-SSR markers for species identification and cross-genus amplification in sambar, sika, and red deer

Deer genera around the globe are threatened by anthropogenic interference. The translocation of alien species and their subsequent genetic introgression into indigenous deer populations is particularly harmful to the species of greatest conservation concern. Products derived from deer, including venison and antler velvet, are also at risk of fraudulent labeling. The current molecular markers used to genetically identify deer species were developed from genome sequences and have limited applicability for cross-species amplification. The absence of efficacious diagnostic techniques for identifying deer species has hampered conservation and wildlife crime investigation efforts. Expressed sequence tag-simple sequence repeat (EST-SSR) markers are reliable tools for individual and species identification, especially in terms of cross-species genotyping. We conducted transcriptome sequencing of sambar (Rusa unicolor) antler velvet and acquired 11,190 EST-SSRs from 65,074 newly assembled unigenes. We identified a total of 55 unambiguous amplicons in sambar (n = 45), which were selected as markers to evaluate cross-species genotyping in sika deer (Cervus nippon, n = 30) and red deer (Cervus elaphus, n = 46), resulting in cross-species amplification rates of 94.5% and 89.1%, respectively. Based on polymorphic information content (>0.25) and genotyping fidelity, we selected 16 of these EST-SSRs for species identification. This marker set revealed significant genetic differentiation based on the fixation index and genetic distance values. Principal coordinate analysis and STRUCTURE analysis revealed distinct clusters of species and clearly identified red-sika hybrids. These markers showed applicability across different genera and proved suitable for identification and phylogenetic analyses across deer species.

Fine-scale genetic structure reflects limited and coordinated dispersal in the colonial monk parakeet, Myiopsitta monachus

The genetic structure of animal populations has considerable behavioural, ecological and evolutionary implications and may arise from various demographic traits. Here, we use observational field data and molecular genetics to determine the genetic structure of an invasive population of monk parakeets, Myiopsitta monachus, at a range of spatial scales, and investigate the demographic processes that generate the observed structure. Monk parakeets construct large nests that can house several pairs occupying separate chambers; these nests are often aggregated within nesting trees. We determined patterns of relatedness within compound nests, within nesting trees and between trees. Spatial autocorrelation analyses of pairwise genetic relatedness revealed fine-scale genetic structure with relatives of both sexes spatially clustered within, but not beyond, nesting trees. In addition, males were more related to males sharing their compound nests than to other males occupying the same nesting tree. By contrast, males and females within compound nests were not significantly more closely related than elsewhere in the same tree, and we found no evidence for inbreeding. Adults showed high breeding site fidelity between years despite considerable disturbance of nest sites. Natal dispersal was female-biased, but dispersal distances were relatively short with some natal philopatry observed in both sexes. Sibling coalitions, typically of males, were observed amongst both philopatric and dispersing birds. Our results show significant clustering of kin within compound nests and nesting trees resulting from limited and coordinated natal dispersal, with subsequent breeding site fidelity. The resulting genetic structure has implications for social behaviour in this unusual parrot species.

Development of EST-Molecular Markers from RNA Sequencing for Genetic Management and Identification of Growth Traits in Potato Grouper (Epinephelus tukula)

Simple Summary

The potato grouper is a novel aquaculture species in Taiwan. Due to the lack of genetic information concerning this species, we have developed molecular markers based on transcriptome sequencing and further characterized their association with gene diversity and growth traits of this species. Ultimately, these markers could be utilized as accurate and efficient tools for genetic management and marker-assisted selection of potato grouper with distinct growth traits.

Abstract

The accuracy and efficiency of marker-assisted selection (MAS) has been proven for economically critical aquaculture species. The potato grouper (Epinephelus tukula), a novel cultured grouper species in Taiwan, shows large potential in aquaculture because of its fast growth rate among other groupers. Because of the lack of genetic information for the potato grouper, the first transcriptome and expressed sequence tag (EST)-derived simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers were developed. Initially, the transcriptome was obtained from seven cDNA libraries by using the Illumina platform. De novo transcriptome of the potato grouper yielded 51.34 Gb and 111,490 unigenes. The EST-derived SSR and SNP markers were applied in genetic management, in parentage analysis, and to discover the functional markers of economic traits. The F₁ juveniles were identified as siblings from one pair of parents (80 broodstocks). Fast- and slow-growth individuals were analyzed using functional molecular markers and through their association with growth performance. The results revealed that two SNPs were correlated with growth traits. The transcriptome database obtained in this study and its derived SSR and SNP markers may be applied not only for MAS but also to maintain functional gene diversity in the novel cultured grouper.

Habitat fragmentation differentially shapes neutral and immune gene variation in a tropical bird species

Habitat fragmentation is a major cause of biodiversity loss, responsible for an alteration of intraspecific patterns of neutral genetic diversity and structure. Although neutral genetic variation can be informative for demographic inferences, it may be a poor predictor of adaptive genetic diversity and thus of the consequences of habitat fragmentation on selective evolutionary processes. In this context, we contrasted patterns of genetic diversity and structure of neutral loci (microsatellites) and immune genes (i.e., toll-like receptors) in an understorey bird species, the wedge-billed woodcreeper Glyphorynchus spirurus. The objectives were (1) to investigate forest fragmentation effects on population genetic diversity, (2) to disentangle the relative role of demography (genetic drift and migration) and selection, and (3) to assess whether immunogenetic patterns could be associated with variation of ectoparasite (i.e., ticks) pressures. Our results revealed an erosion of neutral genetic diversity and a substantial genetic differentiation among fragmented populations, resulting from a decrease in landscape connectivity and leading to the divergence of distinct genetic pools at a small spatial scale. Patterns of genetic diversity observed for TLR4 and TLR5 were concordant with neutral genetic patterns, whereas those observed for TLR3 and TLR21 were discordant. This result underlines that the dominant evolutionary force shaping immunogenetic diversity (genetic drift vs. selection) may be different depending on loci considered. Finally, tick prevalence was higher in fragmented environments. We discussed the hypothesis that pathogen selective pressures may contribute to maintain adaptive genetic diversity despite the negative demographic effect of habitat fragmentation on neutral genetic diversity.

Hemocytes transcriptomes reveal metabolism changes and detoxification mechanisms in response to ammonia stress in Octopus minor

Ammonia is one of the major aquatic environmental pollutants that can bring detrimental effects to the growth and survival of aquatic organisms. However, the molecular mechanisms of ammonia toxicity and ammonia excretion in marine invertebrates especially mollusks are still poorly understood. Cephalopods are exclusively ammonotelic with high protein metabolism and ammonia excretion rate, making this taxonomic group an ideal specimen to explore the ammonia detoxification mechanism. In this study, comparative transcriptomes were employed to investigate the transcriptional changes of O. minor in responses to acute ammonia exposure. A total of 63,237 unigenes with an average length of 811 bp were discovered and 25,708 unigenes were successfully annotated. The transcription of 1845 genes were significantly changed after ammonia stress, including 315 up-regulated genes and 1530 down-regulated genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis based on differentially expressed genes (DEGs) revealed that 44 GO terms and 55 KEGG pathways were over-represented. Notably, a large number of genes involved in immune defense, citric acid (TCA) cycle, oxidative phosphorylation and amino acid metabolisms were significantly down-regulated, indicating the decelerated energy production and amino acid rate in response to acute ammonia stress. These results provide new insights into the potential molecular mechanism of ammonia detoxification on transcriptomic level and will facilitate further mechanism studies on mollusks.

Development of 49 novel microsatellite markers from Next-generation sequencing data and a robust method for parentage tests in the emu (Dromaius novaehollandiae)

The emu is a useful and new breed of poultry, but their genetic improvement has not advanced yet due to their very recent domestication. Pedigree information is difficult to record because of their complex reproduction system (polyandry). To identify parent-offspring relationships in the emu, parentage test based on polymorphic DNA markers have to be developed. In this study, we isolated more than 25,000 microsatellite (simple sequence repeat, SSR) regions from Next-generation sequencing data via the QDD pipeline and developed 49 SSR markers with polymorphism in the Japanese farmed emu. The dinucleotide motifs, (AC)_n, (AT)_n and (AG)_n, were the most frequently detected and were found on 10,167 (38.55%), 8,114 (30.76%) and 4,796 (18.18%) contigs, respectively. Forty-nine novel SSR markers were characterized in 20 individuals and showed N_A ranged from 2 to 12, with an average of 4.2. H_E/H_O ranged from 0.389/0.071 to 0.702/1.000 with an average of 0.601/0.515. PIC value ranged from 0.059 to 0.886 with an average of 0.528, and 17 of 49 markers showed a higher polymorphism than 0.500. Thirty-four individuals were genotyped using 12 markers, and CERVUS simulations based on genotype showed that parents of all offspring were identified with 0.9995-1.0 probability. Thus, 49 novel SSR markers and a robust method for parentage test for the Japanese emu were developed.

Genome-based development of 15 microsatellite markers in fluorescent multiplexes for parentage testing in captive tigers

As one of the most endangered species, tiger (Panthera tigris) inbreeding has become an urgent issue to address. Using a microsatellite (short tandem repeat, STR) identification system, paternity testing may be helpful to avoid inbreeding in captive breeding programs. In this study, we developed a genome-based identification system named tiger pedigree identification multiplex system (TPI-plex). By analyzing the entire tiger genome, 139,967 STR loci were identified and 12.76% of these displayed three to six alleles among three re-sequenced individual tiger genomes. A total of 204 candidate STRs were identified and screened with a reference population containing 31 unrelated captive tigers. Of these, 15 loci were chosen for inclusion in the multiplex panel. The mean allele number and mean expected heterozygosity (He) were 7.3333 and 0.7789, respectively. The cumulative probability of exclusion (CPE) and total probability of discrimination power (TDP) reached 0.999999472 and 0.999999999999995, respectively. The results showed that the TPI-plex system can be applied in routine pedigree identification for captive tigers. We also added a sex identification marker named TAMEL into the TPI-plex for sex determination.

Estimating genetic diversity and population structure of 22 chicken breeds in Asia using microsatellite markers

Objective

Estimating the genetic diversity and structures, both within and among chicken breeds, is critical for the identification and conservation of valuable genetic resources. In chickens, microsatellite (MS) marker polymorphisms have previously been widely used to evaluate these distinctions. Our objective was to analyze the genetic diversity and relationships among 22 chicken breeds in Asia based on allelic frequencies.

Methods

We used 469 genomic DNA samples from 22 chicken breeds from eight Asian countries (South Korea, KNG, KNB, KNR, KNW, KNY, KNO; Laos, LYO, LCH, LBB, LOU; Indonesia, INK, INS, ING; Vietnam, VTN, VNH; Mongolia, MGN; Kyrgyzstan, KGPS; Nepal, NPS; Sri Lanka, SBC) and three imported breeds (RIR, Rhode Island Red; WLG, White Leghorn; CON, Cornish). Their genetic diversity and phylogenetic relationships were analyzed using 20 MS markers.

Results

In total, 193 alleles were observed across all 20 MS markers, and the number of alleles ranged from 3 (MCW0103) to 20 (LEI0192) with a mean of 9.7 overall. The NPS breed had the highest expected heterozygosity (H_exp, 0.718±0.027) and polymorphism information content (PIC, 0.663±0.030). Additionally, the observed heterozygosity (H_obs) was highest in LCH (0.690±0.039), whereas WLG showed the lowest H_exp (0.372±0.055), H_obs (0.384±0.019), and PIC (0.325±0.049). Nei’s DA genetic distance was the closest between VTN and VNH (0.086), and farthest between KNG and MGN (0.503). Principal coordinate analysis showed similar results to the phylogenetic analysis, and three axes explained 56.2% of the variance (axis 1, 19.17%; 2, 18.92%; 3, 18.11%). STRUCTURE analysis revealed that the 22 chicken breeds should be divided into 20 clusters, based on the highest ΔK value (46.92).

Conclusion

This study provides a basis for future genetic variation studies and the development of conservation strategies for 22 chicken breeds in Asia.

Genetic diversity and population structure of indigenous chicken in Rwanda using microsatellite markers

Rwanda has about 4.5 million of indigenous chicken (IC) that are very low in productivity. To initiate any genetic improvement programme, IC needs to be accurately characterized. The key purpose of this study was to ascertain the genetic diversity of IC in Rwanda using microsatellite markers. Blood samples of IC sampled from 5 agro-ecological zones were collected from which DNA was extracted, amplified by PCR and genotyped using 28 microsatellite markers. A total of 325 (313 indigenous and 12 exotic) chickens were genotyped and revealed a total number of 305 alleles varying between 2 and 22 with a mean of 10.89 per locus. One hundred eighty-six (186) distinct alleles and 60 private alleles were also observed. The frequency of private alleles was highest in samples from the Eastern region, whereas those from the North West had the lowest. The influx of genes was lower in the Eastern agro-ecological zone than the North West. The mean observed heterozygosity was 0.6155, whereas the average expected heterozygosity was 0.688. The overall inbreeding coefficient among the population was 0.040. Divergence from the Hardy-Weinberg equilibrium was significant (p<0.05) in 90% of loci in all the populations. The analysis of molecular variance revealed that about 92% of the total variation originated from variation within populations. Additionally, the study demonstrated that IC in Rwanda could be clustered into four gene groups. In conclusion, there was considerable genetic diversity in IC in Rwanda, which represents a crucial genetic resource that can be conserved or optimized through genetic improvement.

Phenotypic and Genetic Variation of an Interspecific Centaurium Hybrid (Gentianaceae) and Its Parental Species

Interspecific hybridization is one of the major actuators of evolutionary changes in plants. As the result of allopolyploid hybridization, offspring may gain different ploidy levels in comparison to parental species, which can provide them instant reproductive isolation. Two tetraploid sister species, Centaurium erythraea and C. littorale, readily cross-fertilize, resulting in hybrids of various ploidy. In northern Serbia, two stable populations of a hexaploid taxon C. pannonicum have been documented. It has been proposed previously that this taxon emerged after an interspecific hybridization event between two tetraploid sister-species: C. erythraea and C. littorale subsp. compressum. The existing populations of the hybridogenic taxon, as well as neighboring populations of the two parental taxa were here characterized by both morphometrics and molecular markers (EST-SSR and trnL-F). Three leaf and two flower characteristics were found to be informative in delimitation of the parental taxa and in their discernment from hybrid individuals, the latter having intermediate values. Eight microsatellite markers were found to have good ability to distinguish studied taxa, placing C. pannonicum in closer relationship with C. erythraea. Conversely, trnL-F plastid marker nominated C. littorale subsp. compressum to be the donor of the C. pannonicum plastid DNA. Reproductive isolation of the hexaploid hybrid individuals from the parental species should be examined as the next logical step in describing the new species.

Transcriptome analysis and discovery of genes involved in immune pathways in Solen strictus (Gould, 1861) under Vibrio anguillarum

The Solen strictus (Gould, 1861) has been recognized as an important marine economic bivalve in Eastern and Southeast Asia. To gain a better understanding of the S. strictus immune system and its related genes in response to bacterial infections, we performed a comparative gene transcription analysis from S. strictus with Vibrio anguillarum through RNA-Seq technology, meanwhile the differentially expressed genes (DEGs) were investigated. After assembly, a total of 195,774 transcripts with an average length of 996 bp were obtained. Total 153,038 unigenes were annotated in the nr, Swiss-Prot, KEGG, COG, KOG, GO and Pfam databases, and 56,597 unigenes (36.98%) were annotated in at least one database. After bacterial challenge, there were 1588 significant differentially expressed genes (DEGs) between the challenged and control groups, including 999 up-regulated and 589 down-regulated genes. All the DEGs were classified into three gene ontology categories, and allocated to 225 KEGG pathways. Immune-related genes were detected from immune system pathways among the top 20 enriched pathways, such as Toll-like receptor signaling, RIG-I-like receptor signaling and NOD-like receptor signaling pathway. In addition, 56,079 potential simple sequence repeats (SSRs) and 1,031,521 candidate single nucleotide polymorphisms (SNPs) were detected and identified in the S. strictus transcriptome. Results of the present study will provide valuable theoretical resources for future genetic and genomic research on S. strictus. The research results will be helpful for improving the efficiency and quality of artificial breeding, establishing genetic linkage map, and enhancing health management for this species.

Hybrid identification for Glycine max and Glycine soja with SSR markers and analysis of salt tolerance

Glycine max cultivars Lee68, Nannong 1138-2, and Nannong 8831 were used as the female parents, and hybrid lines (F5) 4,111, 4,076 (N23674 × BB52), 3,060 (Lee68 × N23227), and 185 (Jackson × BB52) that selected for salt tolerance generation by generation from the cross combination of G. max and G. soja were used as the male parents, 11 (A-K) backcrosses or three-way crosses were designed and 213 single hybrids were harvested. The optimized soybean simple sequence repeat (SSR)-polymerase chain reaction (PCR) system was used to analyze the SSR polymorphism of above parental lines and get the parental co-dominant SSR markers for hybrid identification, and in which 30 true hybrids were gained. The true hybrids (G1, G3, G9, G12, G13, G16) of G cross combination were chosen as the representative for the salt tolerance test, and the results showed that, as exposed to salt stress, the seedlings of G9 line displayed higher salt tolerant coefficient, relative growth rate, and dry matter accumulation, when compared with their female parent Nannong 1138-2, and even performed equally strong salt tolerance as the male parent 3,060. It provides a feasible method of the combination of molecular SSR markers and simple physiological parameters to identify the true hybrids of G. max and G. soja, and to innovate the salt-tolerant soybean germplasms.

Liver transcriptome analysis of the Sparus macrocephlus in response to Vibrio parahaemolyticus infection

The black seabream (Sparus macrocephlus) is an economically pivotal aquaculture species cultured in China and Southeast Asian countries. To understand the molecular immune mechanisms underlying the response to Vibrio parahaemolyticus, a comparative gene transcription analysis were performed with utilized fresh livers of V. parahaemolyticus-immunized Sparus macrocephlus with a control group through RNA-Seq technology. A total of 256663 contigs were obtained after excluded the low-quality sequences and assembly. The average length of contigs collected from this research is 1066.93 bp. Furthermore, blast analysis indicates 30747 contigs were annotated based on homology with matches in the NT, NR, gene, and string databases. A gene ontology analysis was employed to classify 21598 genes according to three major functional categories: molecular function, cellular component, and biological process. A total of 14470 genes were discovered in 303 KEGG pathways. RSEM and EdgeR were introduced to estimate 3841 genes significantly different expressed (False Discovery Rate<0.001) which includes 4072 up-regulated genes and 3771 down-regulated genes. A significant enrichment analysis of these differentially expressed genes and isogenes were conducted to reveal the major immune-related pathways which refer to the toll-like receptor, complement, coagulation cascades, and chemokine signaling pathways. In addition, 92175 potential simple sequence repeats (SSRs) and 121912 candidate single nucleotide polymorphisms (SNPs) were detected and identified sequencely in the Sparus macrocephlus liver transcriptome. This research characterized a gene expression pattern for normal and the V. parahaemolyticus -immunized Sparus macrocephlus for the first time and not only sheds new light on the molecular mechanisms underlying the host-V. parahaemolyticus interaction but contribute to facilitate future studies on Sparus macrocephlus gene expression and functional genomics.

Population diversity and relatedness in Sugarbirds (Promeropidae: Promerops spp.)

Sugarbirds are a family of two socially-monogamous passerine species endemic to southern Africa. Cape and Gurney’s Sugarbird (Promerops cafer and P. gurneyi) differ in abundance, dispersion across their range and in the degree of sexual dimorphism in tail length, factors that affect breeding systems and potentially genetic diversity. According to recent data, P. gurneyi are in decline and revision of the species’ IUCN conservation status to a threatened category may be warranted. It is therefore necessary to understand genetic diversity and risk of inbreeding in this species. We used six polymorphic microsatellite markers and one mitochondrial gene (ND2) to compare genetic diversity in P. cafer from Helderberg Nature Reserve and P. gurneyi from Golden Gate Highlands National Park, sites at the core of each species distribution. We describe novel universal avian primers which amplify the entire ND2 coding sequence across a broad range of bird orders. We observed high mitochondrial and microsatellite diversity in both sugarbird populations, with no detectable inbreeding and large effective population sizes.

Transcriptome Analysis of Epithelioma Papulosum Cyprini Cells Infected by Reovirus Isolated from Allogynogenetic Silver Crucian Carp

The present study aimed to identify differentially expressed genes (DEGs) and major signal transduction pathways that were related to the immune response of epithelioma papulosum cyprinid (EPC) cells to reoviruses isolated from allogynogenetic silver crucian carp. The study also lays a theoretical foundation for the pathogenesis and immunity of the reovirus, which is helpful to the breeding of cyprinids fish. Reovirus infected and uninfected EPC cells were analyzed by using a new-generation high-throughput sequencing technology. DEGs were identified, annotated, and classified, and the signal pathways involved in the response to reovirus infection were identified by using bioinformatics tool. The data were assembled into 92,101 contigs with an average length of 835.24 bp and an N50 value of 1432 nt. Differential expression analysis of all the genes identified 3316 DEGs at a false discovery rate (FDR) of <0.01 and a fold-change of ≥3, of which 1691 were upregulated genes, 1625 were downregulated, and about 305 were immune-related genes. Gene Ontology (GO) enrichment analysis resulted in the annotation of 3941 GO terms, including 2719 biological processes (37,810 unigenes), 376 cell components (7943 unigenes), and 846 molecular functions (11,750 unigenes). KEGG metabolic pathway analysis matched the DEGs from pre-and post-infection EPC cells to 193 pathways, of which 35 were immune-related, including the Toll-like receptor, cytokine-cytokine receptor interaction, and the JAK-STAT signaling pathways.

Phenotypic divergence despite low genetic differentiation in house sparrow populations

Studying patterns of phenotypic variation among populations can shed light on the drivers of evolutionary processes. The house sparrow (Passer domesticus) is one of the world's most ubiquitous bird species, as well as a successful invader. We investigated phenotypic variation in house sparrow populations across a climatic gradient and in relation to a possible scenario of an invasion. We measured variation in morphological, coloration, and behavioral traits (exploratory behavior and neophobia) and compared it to the neutral genetic variation. We found that sparrows were larger and darker in northern latitudes, in accordance with Bergmann's and Gloger's biogeographic rules. Morphology and behavior mostly differed between the southernmost populations and the other regions, supporting the possibility of an invasion. Genetic differentiation was low and diversity levels were similar across populations, indicating high gene flow. Nevertheless, the southernmost and northern populations differed genetically to some extent. Furthermore, genetic differentiation (F _ST) was lower in comparison to phenotypic variation (P _ST), indicating that the phenotypic variation is shaped by directional selection or by phenotypic plasticity. This study expands our knowledge on evolutionary mechanisms and biological invasions.

Transcriptome analysis of yellow catfish (Pelteobagrus fulvidraco) liver challenged with polyriboinosinic polyribocytidylic acid (poly I:C)

Yellow catfish (Pelteobagrus fulvidraco) is one of the most important economic freshwater species in China. However, infection by bacterial pathogenic diseases has caused high mortality and great economic loss in aquaculture. It is necessary for disease control to know more about the P. fulvidraco immune system and its related genes in response to bacterial or viral infections. In this study, the transcriptomic profiles of liver from P. fulvidraco stimulated by polyriboinosinic polyribocytidylic acid (poly I:C) was analyzed using high-throughput sequencing method. After assembly and annotation, total 67,447 unigenes were acquired, with an average length of 1091 bp. Under the infection of poly I:C, 522 differentially expressed genes (DEGs) were identified, including 307 up-regulated genes and 215 down-regulated genes. To further investigate the immune-related DEGs, Gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment were performed. The result of GO enrichment indicated gene response to external stimulus, regulation of response to stimulus, cellular response to stimulus, immune response and immune system progress. Significant KEGG enrichment analysis identified major immune related pathways. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) analysis revealed that 13 immune response genes were identified to be up-regulated after 12 h of poly I:C stimulation compared to controls. Taken together, the results of our study are beneficial for better understanding of the immune system and defense mechanisms of yellow catfish in response to poly I:C infection.

Transcriptome analysis of Ruditapes philippinarum hepatopancreas provides insights into immune signaling pathways under Vibrio anguillarum infection

The Manila clam, Ruditapes philippinarum, is one of the most economically important aquatic clams that are harvested on a large scale by the mariculture industry in China. However, increasing reports of bacterial pathogenic diseases have had a negative effect on the aquaculture industry of R. philippinarum. In the present study, the two transcriptome libraries of untreated (termed H) and challenged Vibrio anguillarum (termed HV) hepatopancreas were constructed and sequenced from Manila clam using an Illumina-based paired-end sequencing platform. In total, 75,302,886 and 66,578,976 high-quality clean reads were assembled from 101,080,746 and 99,673,538 raw data points from the two transcriptome libraries described above, respectively. Furthermore, 156,116 unigenes were generated from 210,685 transcripts, with an N50 length of 1125 bp, and from the annotated SwissProt, NR, NT, KO, GO, KOG and KEGG databases. Moreover, a total of 4071 differentially expressed unigenes (HV vs H) were detected, including 903 up-regulated and 3168 down-regulated genes. Among these differentially expressed unigenes, 226 unigenes were annotated using KEGG annotation in 16 immune-related signaling pathways, including Toll-like receptor, NF-kappa B, MAPK, NOD-like receptor, RIG-I-like receptor, and the TNF and chemokine signaling pathways. Finally, 20,341 simple sequence repeats (SSRs) and 214,430 potential single nucleotide polymorphisms (SNPs) were detected from the H and HV transcriptome libraries. In conclusion, these studies identified many candidate immune-related genes and signaling pathways and conducted a comparative analysis of the differentially expressed unigenes from Manila clam hepatopancreas in response to V. anguillarum stimulation. These data laid the foundation for studying the innate immune systems and defense mechanisms in R. philippinarum.

Assessment of genetic diversity and phylogenetic relationships of Korean native chicken breeds using microsatellite markers

Objective

This study was conducted to investigate the basic information on genetic structure and characteristics of Korean Native chickens (NC) and foreign breeds through the analysis of the pure chicken populations and commercial chicken lines of the Hanhyup Company which are popular in the NC market, using the 20 microsatellite markers.

Methods

In this study, the genetic diversity and phylogenetic relationships of 445 NC from five different breeds (NC, Leghorn [LH], Cornish [CS], Rhode Island Red [RIR], and Hanhyup [HH] commercial line) were investigated by performing genotyping using 20 microsatellite markers.

Results

The highest genetic distance was observed between RIR and LH (18.9%), whereas the lowest genetic distance was observed between HH and NC (2.7%). In the principal coordinates analysis (PCoA) illustrated by the first component, LH was clearly separated from the other groups. The correspondence analysis showed close relationship among individuals belonging to the NC, CS, and HH lines. From the STRUCTURE program, the presence of 5 clusters was detected and it was found that the proportion of membership in the different clusters was almost comparable among the breeds with the exception of one breed (HH), although it was highest in LH (0.987) and lowest in CS (0.578). For the cluster 1 it was high in HH (0.582) and in CS (0.368), while for the cluster 4 it was relatively higher in HH (0.392) than other breeds.

Conclusion

Our study showed useful genetic diversity and phylogenetic relationship data that can be utilized for NC breeding and development by the commercial chicken industry to meet consumer demands.

The use of microsatellite markers in Neotropical studies of wild birds: a literature review

Despite extensive habitat fragmentation, the Neotropical region possesses 30% of the world´s bird species. Microsatellites have remained one of the most popular genetic markers and have been used in ecological and conservation studies since the 1990's. We conducted a literature review comparing the number of papers published from January 1990 to July 2015 that used microsatellite markers for studies of wild birds in the Neotropical region, USA and some European countries. We assigned the articles to three categories of studies: population genetics, animal behavior/kinship analysis and the development of species-specific bird microsatellite markers. We also compared the studies in the Neotropics that used heterologous versus species-specific markers and provide a list of heterologous markers of utility in multiple birds. Despite the rich bird fauna in the Neotropics, the number of articles published represents only 5.6% of that published by the USA and selected European countries. Within the Neotropical region, Brazil possessed 60.5% of the total papers published, with the remaining 39.5% shared between five countries. We conclude that the lack of specialized laboratories and resources still represents a limit to microsatellite-based genetic studies of birds within the Neotropical region. To overcome these limitations, we suggest the use of heterologous microsatellite markers as a cost-effective and time-effective tool to assist ecological studies of wild birds.

Level of genetic differentiation affects relative performances of expressed sequence tag and genomic SSRs

Microsatellites, also called simple sequence repeats (SSRs), are markers of choice to estimate relevant parameters for conservation genetics, such as migration rates, effective population size and kinship. Cross-amplification of SSRs is the simplest way to obtain sets of markers, and highly conserved SSRs have recently been developed from expressed sequence tags (EST) to improve SSR cross-species utility. As EST-SSRs are located in coding regions, the higher stability of their flanking regions reduces the frequency of null alleles and improves cross-species amplification. However, EST-SSRs have generally less allelic variability than genomic SSRs, potentially leading to differences in estimates of population genetic parameters such as genetic differentiation. To assess the potential of EST-SSRs in studies of within-species genetic diversity, we compared the relative performance of EST- and genomic SSRs following a multispecies approach on passerine birds. We tested whether patterns and levels of genetic diversity within and between populations assessed from EST- and from genomic SSRs are congruent, and we investigated how the relative efficiency of EST- and genomic SSRs is influenced by levels of differentiation. EST- and genomic SSRs ensured comparable inferences of population genetic structure in cases of strong genetic differentiation, and genomic SSRs performed slightly better than EST-SSRs when differentiation is moderate. However and interestingly, EST-SSRs had a higher power to detect weak genetic structure compared to genomic SSRs. Our study attests that EST-SSRs may be valuable molecular markers for conservation genetic studies in taxa such as birds, where the development of genomic SSRs is impeded by their low frequency.

The broad-snouted caiman population recovery in Argentina. A case of genetics conservation

Caiman latirostriswild populations have suffered a drastic reduction in the past, and for that reason, a management and monitoring plan was applied since 1990 in Santa Fe, Argentina in order to achieve population recovery. Although ranching system has a noteworthy success in terms of population size recovering, there is no information about the estimation of population genetic parameters. In particular, the consequence of the bottleneck underwent by these populations has not been assessed. We evaluated variability and genetic structure ofC. latirostrispopulations from Santa Fe through time, using microsatellites and mitochondrial DNA. Population genetic parameters were compared among four sites and three different periods to assess the impact of management activities, and effective population size was estimated in order to detect bottleneck events. We observed an increase in microsatellite variability and low genetic variability in mitochondrial lineages through time. Variability estimates are similar among sites in each sampling period; and there is scarce differentiation among them. The genetic background of each sampling site has changed through time; we assume this fact may be due to entry of individuals of different origin, through management and repopulation activities. Moreover, taking into account the expected heterozygosity and effective population size values, it can be assumed that bottleneck events indeed have occurred in the recent past. Our results suggest that, in addition to increasing population size, genetic variability of the species has been maintained. However, the information is still incomplete, and regular monitoring should continue in order to arrive to solid conclusions.

Habitat specialization predicts genetic response to fragmentation in tropical birds

Habitat fragmentation is one of the most severe threats to biodiversity as it may lead to changes in population genetic structure, with ultimate modifications of species evolutionary potential and local extinctions. Nonetheless, fragmentation does not equally affect all species and identifying which ecological traits are related to species sensitivity to habitat fragmentation could help prioritization of conservation efforts. Despite the theoretical link between species ecology and extinction proneness, comparative studies explicitly testing the hypothesis that particular ecological traits underlies species-specific population structure are rare. Here, we used a comparative approach on eight bird species, co-occurring across the same fragmented landscape. For each species, we quantified relative levels of forest specialization and genetic differentiation among populations. To test the link between forest specialization and susceptibility to forest fragmentation, we assessed species responses to fragmentation by comparing levels of genetic differentiation between continuous and fragmented forest landscapes. Our results revealed a significant and substantial population structure at a very small spatial scale for mobile organisms such as birds. More importantly, we found that specialist species are more affected by forest fragmentation than generalist ones. Finally, our results suggest that even a simple habitat specialization index can be a satisfying predictor of genetic and demographic consequences of habitat fragmentation, providing a reliable practical and quantitative tool for conservation biology.

A Single Transcriptome of a Green Toad (Bufo viridis) Yields Candidate Genes for Sex Determination and -Differentiation and Non-Anonymous Population Genetic Markers

Large genome size, including immense repetitive and non-coding fractions, still present challenges for capacity, bioinformatics and thus affordability of whole genome sequencing in most amphibians. Here, we test the performance of a single transcriptome to understand whether it can provide a cost-efficient resource for species with large unknown genomes. Using RNA from six different tissues from a single Palearctic green toad (Bufo viridis) specimen and Hiseq2000, we obtained 22,5 Mio reads and publish >100,000 unigene sequences. To evaluate efficacy and quality, we first use this data to identify green toad specific candidate genes, known from other vertebrates for their role in sex determination and differentiation. Of a list of 37 genes, the transcriptome yielded 32 (87%), many of which providing the first such data for this non-model anuran species. However, for many of these genes, only fragments could be retrieved. In order to allow also applications to population genetics, we further used the transcriptome for the targeted development of 21 non-anonymous microsatellites and tested them in genetic families and backcrosses. Eleven markers were specifically developed to be located on the B. viridis sex chromosomes; for eight markers we can indeed demonstrate sex-specific transmission in genetic families. Depending on phylogenetic distance, several markers, which are sex-linked in green toads, show high cross-amplification success across the anuran phylogeny, involving nine systematic anuran families. Our data support the view that single transcriptome sequencing (based on multiple tissues) provides a reliable genomic resource and cost-efficient method for non-model amphibian species with large genome size and, despite limitations, should be considered as long as genome sequencing remains unaffordable for most species.

Contrasted patterns of genetic differentiation across eight bird species in the Lesser Antilles

Archipelagoes are considered as "natural laboratories" for studying processes that shape the distribution of diversity. The Lesser Antilles provide a favorable geographical context for divergence to occur. However, although morphological subspecies have been described across this archipelago in numerous avian species, the potential for the Lesser Antilles in driving intra-specific genetic divergence in highly mobile organisms such as birds remains understudied. Here, we assessed level of intra-specific genetic diversity and differentiation between three islands of the Lesser Antilles (Guadeloupe, Dominica and Martinique) using a multi-species approach on eight bird species. For each species, we built a set of microsatellite markers from cross-species amplifications. Significant patterns of inter-island and/or within-island genetic differentiation were detected in all species. However, levels of intra-specific genetic differentiation among the eight bird species were not always consistent with the boundaries of subspecies previously described in the sampled islands. These results suggest different histories of colonization/expansion and/or different species-specific ecological traits affecting gene flow, advocating for multi-species studies of historical and contemporary factors shaping the distribution of diversity on islands.

Transcriptome analysis of the Larimichthys crocea liver in response to Cryptocaryon irritans

The large yellow croaker (Larimichthys crocea) is an economically important marine fish cultured in China and East Asian countries and is facing a serious threat from Cryptocaryon irritans, which is a protozoan ectoparasite that infects most reared marine fish species. To understand the molecular immune mechanisms underlying the response to C. irritans, we first performed a comparative gene transcription analysis using livers from C. irritans-immunized L. croceas and from a control group through RNA-Seq technology. After the removal of low-quality sequences and assembly, 51360 contigs were obtained, with an average length of 1066.93 bp. Further, a blast analysis indicates that 30747 contigs can be annotated based on homology with matches in the NT, NR, gene, and string databases. A gene ontology analysis was used to classify 21598 genes according to three major functional categories: molecular function, cellular component, and biological process. Moreover, 14470 genes were found in 303 KEGG pathways. We used RSEM and EdgeR to determine that 3841 genes were significantly differentially expressed (FDR < 0.001), including 2129 up-regulated genes and 1712 down-regulated genes. A significant enrichment analysis of these differentially expressed genes and isogenes revealed major immune-related pathways, including the toll-like receptor, complement and coagulation cascades, and chemokine signaling pathways. In addition, 28748 potential simple sequence repeats (SSRs) were detected from 12776 transcripts, and 62992 candidate single nucleotide polymorphisms (SNPs) were identified in the L. croceas liver transcriptome. This study characterized a gene expression pattern for normal and C. irritans-immunized L. croceas for the first time and not only sheds new light on the molecular mechanisms underlying the host-C. irritans interaction but also facilitates future studies on L. croceas gene expression and functional genomics.

Antagonistic effect of helpers on breeding male and female survival in a cooperatively breeding bird

1. Cooperatively breeding species are typically long lived and hence, according to theory, are expected to maximize their lifetime reproductive success through maximizing survival. Under these circumstances, the presence of helpers could be used to lighten the effort of current reproduction for parents to achieve higher survival. 2. In addition, individuals of different sexes and ages may follow different strategies, but whether male and female breeders and individuals of different ages benefit differently from the presence of helpers has often been overlooked. Moreover, only one study that investigated the relationship between parental survival and the presence of helpers used capture-mark-recapture analyses (CMR). These methods are important since they allow us to account for the non-detection of individuals that are alive in the population but not detected, and thus, the effects on survival and recapture probability to be disentangled. 3. Here, we used multi-event CMR methods to investigate whether the number of helpers was associated with an increase in survival probability for male and female breeders of different ages in the sociable weaver Philetairus socius. In this species, both sexes reduce their feeding rate in the presence of helpers. We therefore predicted that the presence of helpers should increase the breeders' survival in both sexes, especially early in life when individuals potentially have more future breeding opportunities. In addition, sociable weaver females reduce their investment in eggs in the presence of helpers, so we predicted a stronger effect of helpers on female than male survival. 4. As expected we found that females had a higher survival probability when breeding with more helpers. Unexpectedly, however, male survival probability decreased with increasing number of helpers. This antagonistic effect diminished as the breeders grew older. 5. These results illustrate the complexity of fitness costs and benefits underlying cooperative behaviours and how these may vary with the individuals' sex and age. They also highlight the need for further studies on the sex-specific effects of helpers on survival.

Isolation and characterization of new microsatellite markers for application in population genetic studies of Caiman latirostris and related species

Wild populations ofCaiman latirostrisare subject to sustainable use programs in Argentina, becoming a species with important impact in the regional economy, based in their skin and meat. Genetic studies are fundamental to acquire information on important parameters for conservation and management, which may be obtained from analysis of molecular markers. Some microsatellites have been previously isolated in this species, but due to some difficulties in using them, we obtained new ones using Next Generation Sequencing approach. This study reports eight new microsatellites forC. latirostrisand tests their utility in a related species,Caiman yacare, with successful application in population genetics and mating systems studies. In addition, we shared data about a novel and fast bioinformatics tool to find microsatellites and to design their corresponding primers.

Recent insertion/deletion (reINDEL) mutations: increasing awareness to boost molecular-based research in ecology and evolution

Today, the comparative analysis of DNA molecules mainly uses information inferred from nucleotide substitutions. Insertion/deletion (INDEL) mutations, in contrast, are largely considered uninformative and discarded, due to our lacking knowledge on their evolution. However, including rather than discarding INDELs would be relevant to any research area in ecology and evolution that uses molecular data. As a practical approach to better understanding INDEL evolution in general, we propose the study of recent INDEL (reINDEL) mutations – mutations where both ancestral and derived state are seen in the sample. The precondition for reINDEL identification is knowledge about the pedigree of the individuals sampled. Sound reINDEL knowledge will allow the improved modeling needed for including INDELs in the downstream analysis of molecular data. Both microsatellites, currently still the predominant marker system in the analysis of populations, and sequences generated by next-generation sequencing, a promising and rapidly developing range of technologies, offer the opportunity for reINDEL identification. However, a 2013 sample of animal microsatellite studies contained unexpectedly few reINDELs identified. As most likely explanation, we hypothesize that reINDELs are underreported rather than absent and that this underreporting stems from common reINDEL unawareness. If our hypothesis applies, increased reINDEL awareness should allow gathering data rapidly. We recommend the routine reporting of either the absence or presence of reINDELs together with standardized key information on the nature of mutations when they are detected and the use of the keyword “reINDEL” to increase visibility in both instances of successful and unsuccessful search.

Continental-scale patterns of pathogen prevalence: a case study on the corncrake

Pathogen infections can represent a substantial threat to wild populations, especially those already limited in size. To determine how much variation in the pathogens observed among fragmented populations is caused by ecological factors, one needs to examine systems where host genetic diversity is consistent among the populations, thus controlling for any potentially confounding genetic effects. Here, we report geographic variation in haemosporidian infection among European populations of corncrake. This species now occurs in fragmented populations, but there is little genetic structure and equally high levels of genetic diversity among these populations. We observed a longitudinal gradient of prevalence from western to Eastern Europe negatively correlated with national agricultural yield, but positively correlated with corncrake census population sizes when only the most widespread lineage is considered. This likely reveals a possible impact of local agriculture intensity, which reduced host population densities in Western Europe and, potentially, insect vector abundance, thus reducing the transmission of pathogens. We conclude that in the corncrake system, where metapopulation dynamics resulted in variations in local census population sizes, but not in the genetic impoverishment of these populations, anthropogenic activity has led to a reduction in host populations and pathogen prevalence.

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.

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.

Genetic and phenotypic divergence in an island bird: isolation by distance, by colonization or by adaptation?

Discerning the relative roles of adaptive and nonadaptive processes in generating differences among populations and species, as well as how these processes interact, is a fundamental aim in biology. Both genetic and phenotypic divergence across populations can be the product of limited dispersal and gradual genetic drift across populations (isolation by distance), of colonization history and founder effects (isolation by colonization) or of adaptation to different environments preventing migration between populations (isolation by adaptation). Here, we attempt to differentiate between these processes using island populations of Berthelot's pipit (Anthus berthelotii), a passerine bird endemic to three Atlantic archipelagos. Using microsatellite markers and approximate Bayesian computation, we reveal that the northward colonization of this species ca. 8500 years ago resulted in genetic bottlenecks in the colonized archipelagos. We then show that high levels of genetic structure exist across archipelagos and that these are consistent with a pattern of isolation by colonization, but not with isolation by distance or adaptation. Finally, we show that substantial morphological divergence also exists and that this is strongly concordant with patterns of genetic structure and bottleneck history, but not with environmental differences or geographic distance. Overall, our data suggest that founder effects are responsible for both genetic and phenotypic changes across archipelagos. Our findings provide a rare example of how founder effects can persist over evolutionary timescales and suggest that they may play an important role in the early stages of speciation.

Characterization of the house sparrow (Passer domesticus) transcriptome: a resource for molecular ecology and immunogenetics

The house sparrow (Passer domesticus) is an important model species in ecology and evolution. However, until recently, genomic resources for molecular ecological projects have been lacking in this species. Here, we present transcriptome sequencing data (RNA-Seq) from three different house sparrow tissues (spleen, blood and bursa). These tissues were specifically chosen to obtain a diverse representation of expressed genes and to maximize the yield of immune-related gene functions. After de novo assembly, 15,250 contigs were identified, representing sequence data from a total of 8756 known avian genes (as inferred from the closely related zebra finch). The transcriptome assembly contain sequence data from nine manually annotated MHC genes, including an almost complete MHC class I coding sequence. There were 407, 303 and 68 genes overexpressed in spleen, blood and bursa, respectively. Gene ontology terms related to ribosomal function were associated with overexpression in spleen and oxygen transport functions with overexpression in blood. In addition to the transcript sequences, we provide 327 gene-linked microsatellites (SSRs) with sufficient flanking sequences for primer design, and 3177 single-nucleotide polymorphisms (SNPs) within genes, that can be used in follow-up molecular ecology studies of this ecological well-studied species.

Stronger transferability but lower variability in transcriptomic- than in anonymous microsatellites: evidence from Hylid frogs

A simple way to quickly optimize microsatellites in nonmodel organisms is to reuse loci available in closely related taxa; however, this approach can be limited by the stochastic and low cross-amplification success experienced in some groups (e.g. amphibians). An efficient alternative is to develop loci from transcriptome sequences. Transcriptomic microsatellites have been found to vary in their levels of cross-species amplification and variability, but this has to date never been tested in amphibians. Here, we compare the patterns of cross-amplification and levels of polymorphism of 18 published anonymous microsatellites isolated from genomic DNA vs. 17 loci derived from a transcriptome, across nine species of tree frogs (Hyla arborea and Hyla cinerea group). We established a clear negative relationship between divergence time and amplification success, which was much steeper for anonymous than transcriptomic markers, with half-lives (time at which 50% of the markers still amplify) of 1.1 and 37 My, respectively. Transcriptomic markers are significantly less polymorphic than anonymous loci, but remain variable across diverged taxa. We conclude that the exploitation of amphibian transcriptomes for developing microsatellites seems an optimal approach for multispecies surveys (e.g. analyses of hybrid zones, comparative linkage mapping), whereas anonymous microsatellites may be more informative for fine-scale analyses of intraspecific variation. Moreover, our results confirm the pattern that microsatellite cross-amplification is greatly variable among amphibians and should be assessed independently within target lineages. Finally, we provide a bank of microsatellites for Palaearctic tree frogs (so far only available for H. arborea), which will be useful for conservation and evolutionary studies in this radiation.

Advancing avian behavioral neuroendocrinology through genomics

Genome technologies are transforming all areas of biology, including the study of hormones, brain and behavior. Annotated reference genome assemblies are rapidly being produced for many avian species. Here we briefly review the basic concepts and tools used in genomics. We then consider how these are informing the study of avian behavioral neuroendocrinology, focusing in particular on lessons from the study of songbirds. We discuss the impact of having a complete "parts list" for an organism; the transformational potential of studying large sets of genes at once instead one gene at a time; the growing recognition that environmental and behavioral signals trigger massive shifts in gene expression in the brain; and the prospects for using comparative genomics to uncover the genetic roots of behavioral variation. Throughout, we identify promising new directions for bolstering the application of genomic information to further advance the study of avian brain and behavior.

Characterization and high cross-species transferability of microsatellite markers from the floral transcriptome of Aspidistra saxicola (Asparagaceae)

Recent studies utilizing transcriptome sequences generated by next-generation sequencing (NGS) technologies have demonstrated the ability to rapidly detect and characterize thousands of gene-based microsatellites from different plants. However, these simple sequence repeats (SSRs) were seldom used directly to test interspecific transferability in the populations of closely related species. Aspidistra Ker-Gawl. is a monocot genus with high species richness and diversity in flower structure, but its fresh floral materials are not easy to obtain. Until now, little is known about genetic background in the species of Aspidistra, quite apart from the fearful reduction of their natural habitats. In this study, the floral transcriptome of Aspidistra saxicola was obtained using NGS. Based on these data, a total of 5527 SSRs were identified in the unigenes. Among these SSRs, the proportions of di- and tri-nucleotide repeats were quite close (49.6% verse 46.8%), and the most tri-nucleotide repeats were AGG/CCT followed by AAG/CTT and AGC/GCT in A. saxicola, showing distinct differences with other angiosperm species. To assess genetic diversity in the species of Aspidistra, 48 SSR loci were tested in four available populations of A. elatior. The results revealed that more than a third of the loci were polymorphic. The majority of these primers could be amplified in 24 species representing the main clades of Aspidistra. The primer subsets from transcriptome data proved highly useful for detecting polymorphisms in the related species, supporting the finding that NGS is an efficient approach to molecular marker development at both intra- and interspecies levels, especially in endangered nonmodel species.