Microsatellite loci and the complete mitochondrial DNA sequence characterized through next generation sequencing and de novo genome assembly for the critically endangered orange-bellied parrot, Neophema chrysogaster

Temporal Loss of Genome-Wide and Immunogenetic Diversity in a Near-Extinct Parrot

Loss of genetic diversity threatens a species' adaptive potential and long-term resilience. Predicted to be extinct by 2038, the orange-bellied parrot (Neophema chrysogaster) is a critically endangered migratory bird threatened by numerous viral, bacterial and fungal diseases. The species has undergone multiple population crashes, reaching a low of three wild-born females and 13 males in 2016, and is now represented by only a single wild population and individuals in the captive breeding program. Here we used our high-quality long-read reference genome, and contemporary (N = 19) and historical (N = 16) resequenced genomes from as early as 1829, to track the long-term genomic erosion and immunogenetic diversity decline in this species. 62% of genomic diversity was lost between historical (mean autosomal heterozygosity = 0.00149 ± 0.000699 SD) and contemporary (0.00057 ± 0.000026) parrots. A greater number and length of runs of homozygosity in contemporary samples were also observed. A temporal reduction in the number of alleles at Toll-like receptor genes was found (historical average alleles = 5.78 ± 2.73; contemporary = 3.89 ± 2.10), potentially exacerbating disease susceptibility in the contemporary population. Of particular concern is the new threat of avian influenza strain (HPAI) to Australia. We discuss the conservation implications of our findings and propose that hybridisation and synthetic biology may be required to address the catastrophic loss of genetic diversity that has occurred in this species in order to prevent extinction.

Phylogeography of the white-crowned parrot (Pionus senilis)

Abstract The white-crowned parrot Pionus senilis (von Spix, 1824) is distributed throughout Middle America, inhabiting the Gulf of Mexico coastal area from Tamaulipas (Mexico) to northern Panama. We used mitochondrial data (COI, ND2 and ND4) from 55 specimens to infer phylogenetic relationships, and analyzed the phylogeographic structure, genetic diversity, divergence periods, and historical demography to explore phylogeographic patterns. We found three divergent lineages: two geographically separated by the Isthmus of Tehuantepec, and the third, in Costa Rica by the Nicaragua Depression. The analysis of molecular variance and statistical analyses were consistent with genetically distinct populations. The Central American lineage diverged 1.33 million years ago, whereas the other two lines branched off 1.19 million years ago. This phylogenetic pattern has been reported in other species of Middle American birds.

A Literature Synthesis of Actions to Tackle Illegal Parrot Trade

The order Psittaciformes is one of the most prevalent groups in the illegal wildlife trade. Efforts to understand this threat have focused on describing the elements of the trade itself: actors, extraction rates, and routes. However, the development of policy-oriented interventions also requires an understanding of how research aims and actions are distributed across the trade chain, regions, and species. We used an action-based approach to review documents published on illegal Psittaciformes trade at a global scale to analyze patterns in research aims and actions. Research increased exponentially in recent decades, recording 165 species from 46 genera, with an over representation of American and Australasian genera. Most of the research provided basic knowledge for the intermediary side of the trade chain. Aims such as the identification of network actors, zoonosis control, and aiding physical detection had numerous but scarcely cited documents (low growth rate), while behavior change had the highest growth rate. The Americas had the highest diversity of research aims, contributing with basic knowledge, implementation, and monitoring across the whole trade chain. Better understanding of the supply side dynamics in local markets, actor typology, and actor interactions are needed. Protecting areas, livelihood incentives, and legal substitutes are actions under-explored in parrots, while behavior change is emerging.

Characterization and Comparative Analysis of Complete Mitogenomes of Three Cacatua Parrots (Psittaciformes: Cacatuidae)

Cacatua alba, Cacatua galerita, and Cacatua goffiniana are parrots of the family Cacatuidae. Wild populations of these species are declining with C. alba listed by the International Union for the Conservation of Nature and Natural Resources (IUCN) as Endangered. In this study, complete mitogenomes were sequenced for a comparative analysis among the Cacatua species, and a detailed analysis of the control region. Mitogenome lengths of C. alba,C. galerita, and C. goffiniana were 18,894, 18,900, and 19,084 bp, respectively. They included 13 protein coding genes, two ribosomal RNA genes, 24 transfer RNA genes, three degenerated genes, and two control regions. Ten conserved motifs were found in three domains within each of the two control regions. For an evolution of duplicated control regions of Cacatua, domain I and the 3′ end of domain III experienced an independent evolution, while domain II and most of the regions of domain III was subjected to a concerted evolution. Based on a phylogenetic analysis of 37 mitochondrial genes, the genus Cacatua formed a well-supported, monophyletic, crown group within the Cacatuidae. Molecular dating results showed that Cacatua diverged from other genera of Cacatuinae in the middle of Miocene.

Microsatellite loci development for three catfish species from northwestern South America

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

Next generation sequencing: an application in forensic sciences?

CONTEXT

Over the last few decades, advances in sequencing have improved greatly. One of the most important achievements of Next Generation Sequencing (NGS) is to produce millions of sequence reads in a short period of time, and to produce large sequences of DNA in fragments of any size. Libraries can be generated from whole genomes or any DNA or RNA region of interest without the need to know its sequence beforehand. This allows for looking for variations and facilitating genetic identification.

OBJECTIVES

A deep analysis of current NGS technologies and their application, especially in forensics, including a discussion about the pros and cons of these technologies in genetic identification.

METHODS

A systematic literature search in PubMed, Science Direct and Scopus electronic databases was performed for the period of December 2012 to June 2015.

RESULTS

In the forensic field, one of the main problems is the limited amount of sample available, as well as its degraded state. If the amount of DNA input required for preparing NGS libraries continues to decrease, nearly any sample could be sequenced; therefore, the maximum information from any biological remains could be obtained. Additionally, microbiome typification could be an interesting application to study for crime scene characterisation.

CONCLUSIONS

NGS technologies are going to be crucial for DNA human typing in cases like mass disasters or other events where forensic specimens and samples are compromised and degraded. With the use of NGS it will be possible to achieve the simultaneous analysis of the standard autosomal DNA (STRs and SNPs), mitochondrial DNA, and X and Y chromosomal markers.

Null alleles are ubiquitous at microsatellite loci in the Wedge Clam (Donax trunculus)

Recent studies have reported an unusually high frequency of nonamplifying alleles at microsatellite loci in bivalves. Null alleles have been associated with heterozygous deficits in many studies. While several studies have tested for its presence using different analytical tools, few have empirically tested for its consequences in estimating population structure and differentiation. We characterised 16 newly developed microsatellite loci and show that null alleles are ubiquitous in the wedge clam, Donax trunculus. We carried out several tests to demonstrate that the large heterozygous deficits observed in the newly characterised loci were most likely due to null alleles. We tested the robustness of microsatellite genotyping for population assignment by showing that well-recognised biogeographic regions of the south Atlantic and south Mediterranean coast of Spain harbour genetically different populations.

Transcriptome Analysis of the Tadpole Shrimp (Triops longicaudatus) by Illumina Paired-End Sequencing: Assembly, Annotation, and Marker Discovery

The tadpole shrimp (Triops longicaudatus) is an aquatic crustacean that helps control pest populations. It inhabits freshwater ponds and pools and has been described as a living fossil. T. longicaudatus was officially declared an endangered species South Korea in 2005; however, through subsequent protection and conservation management, it was removed from the endangered species list in 2012. The limited number of available genetic resources on T. longicaudatus makes it difficult to obtain valuable genetic information for marker-aided selection programs. In this study, whole-transcriptome sequencing of T. longicaudatus generated 39.74 GB of clean data and a total of 269,822 contigs using the Illumina HiSeq 2500 platform. After clustering, a total of 208,813 unigenes with an N₅₀ length of 1089 bp were generated. A total of 95,105 unigenes were successfully annotated against Protostome (PANM), Unigene, Eukaryotic Orthologous Groups (KOG), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases using BLASTX with a cut-off of 1E−5. A total of 57,731 unigenes were assigned to GO terms, and 7247 unigenes were mapped to 129 KEGG pathways. Furthermore, 1595 simple sequence repeats (SSRs) were detected from the unigenes with 1387 potential SSR markers. This is the first report of high-throughput transcriptome analysis of T. longicaudatus, and it provides valuable insights for genetic research and molecular-assisted breeding of this important species.

Transcriptome Characterization for Non-Model Endangered Lycaenids, Protantigius superans and Spindasis takanosis, Using Illumina HiSeq 2500 Sequencing

The Lycaenidae butterflies, Protantigius superans and Spindasis takanosis, are endangered insects in Korea known for their symbiotic association with ants. However, necessary genomic and transcriptomics data are lacking in these species, limiting conservation efforts. In this study, the P. superans and S. takanosis transcriptomes were deciphered using Illumina HiSeq 2500 sequencing. The P. superans and S. takanosis transcriptome data included a total of 254,340,693 and 245,110,582 clean reads assembled into 159,074 and 170,449 contigs and 107,950 and 121,140 unigenes, respectively. BLASTX hits (E-value of 1.0 × 10⁻⁵) against the known protein databases annotated a total of 46,754 and 51,908 transcripts for P. superans and S. takanosis. Approximately 41.25% and 38.68% of the unigenes for P. superans and S. takanosis found homologous sequences in Protostome DB (PANM-DB). BLAST2GO analysis confirmed 18,611 unigenes representing Gene Ontology (GO) terms and a total of 5259 unigenes assigned to 116 pathways for P. superans. For S. takanosis, a total of 6697 unigenes were assigned to 119 pathways using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database. Additionally, 382,164 and 390,516 Simple Sequence Repeats (SSRs) were compiled from the unigenes of P. superans and S. takanosis, respectively. This is the first report to record new genes and their utilization for conservation of lycaenid species population and as a reference information for closely related species.

Characterization of the complete mitochondrial genome and a set of polymorphic microsatellite markers through next-generation sequencing for the brown brocket deer Mazama gouazoubira

The complete mitochondrial genome of the brown brocket deer Mazama gouazoubira and a set of polymorphic microsatellite markers were identified by 454-pyrosequencing. De novo genome assembly recovered 98% of the mitochondrial genome with a mean coverage of 9-fold. The mitogenome consisted of 16,356 base pairs that included 13 protein-coding genes, two ribosomal subunit genes, 22 transfer RNAs and the control region, as found in other deer. The genetic divergence between the mitogenome described here and a previously published report was ∼0.5%, with the control region and ND5 gene showing the highest intraspecific variation. Seven polymorphic loci were characterized using 15 unrelated individuals; there was moderate genetic variation across most loci (mean of 5.6 alleles/locus, mean expected heterozygosity = 0.70), with only one locus deviating significantly from Hardy-Weinberg equilibrium, probably because of null alleles. Marker independence was confirmed with tests for linkage disequilibrium. The genetic variation of the mitogenome and characterization of microsatellite markers will provide useful tools for assessing the phylogeography and population genetic patterns in M. gouazoubira, particularly in the context of habitat fragmentation in South America.

Population genetic diversity and hybrid detection in captive zebras

Zebras are members of the horse family. There are three species of zebras: the plains zebra Equus quagga, the Grevy’s zebra E. grevyi and the mountain zebra E. zebra. The Grevy’s zebra and the mountain zebra are endangered, and hybridization between the Grevy’s zebra and the plains zebra has been documented, leading to a requirement for conservation genetic management within and between the species. We characterized 28 microsatellite markers in Grevy’s zebra and assessed cross-amplification in plains zebra and two of its subspecies, as well as mountain zebra. A range of standard indices were employed to examine population genetic diversity and hybrid populations between Grevy’s and plains zebra were simulated to investigate subspecies and hybrid detection. Microsatellite marker polymorphism was conserved across species with sufficient variation to enable individual identification in all populations. Comparative diversity estimates indicated greater genetic variation in plains zebra and its subspecies than Grevy’s zebra, despite potential ascertainment bias. Species and subspecies differentiation were clearly demonstrated and F1 and F2 hybrids were correctly identified. These findings provide insights into captive population genetic diversity in zebras and support the use of these markers for identifying hybrids, including the known hybrid issue in the endangered Grevy’s zebra.

Rearrangement and evolution of mitochondrial genomes in parrots

Mitochondrial genome rearrangements that result in control region duplication have been described for a variety of birds, but the mechanisms leading to their appearance and maintenance remain unclear, and their effect on sequence evolution has not been explored. A recent survey of mitochondrial genomes in the Psittaciformes (parrots) found that control region duplications have arisen independently at least six times across the order. We analyzed complete mitochondrial genome sequences from 20 parrot species, including representatives of each lineage with control region duplications, to document the gene order changes and to examine effects of genome rearrangements on patterns of sequence evolution. The gene order previously reported for Amazona parrots was found for four of the six independently derived genome rearrangements, and a previously undescribed gene order was found in Prioniturus luconensis, representing a fifth clade with rearranged genomes; the gene order resulting from the remaining rearrangement event could not be confirmed. In all rearranged genomes, two copies of the control region are present and are very similar at the sequence level, while duplicates of the other genes involved in the rearrangement show signs of degeneration or have been lost altogether. We compared rates of sequence evolution in genomes with and without control region duplications and did not find a consistent acceleration or deceleration associated with the duplications. This could be due to the fact that most of the genome rearrangement events in parrots are ancient, and additionally, to an effect of body size on evolutionary rate that we found for mitochondrial but not nuclear sequences. Base composition analyses found that relative to other birds, parrots have unusually strong compositional asymmetry (AT- and GC-skew) in their coding sequences, especially at fourfold degenerate sites. Furthermore, we found higher AT skew in species with control region duplications. One potential cause for this compositional asymmetry is that parrots have unusually slow mtDNA replication. If this is the case, then any replicative advantage provided by having a second control region could result in selection for maintenance of both control regions once duplicated.

pSTR Finder: a rapid method to discover polymorphic short tandem repeat markers from whole-genome sequences

Background

Whole-genome sequencing is performed routinely as a means to identify polymorphic genetic loci such as short tandem repeat loci. We have developed a simple tool, called pSTR Finder, which is freely available as a means of identifying putative polymorphic short tandem repeat (STR) loci from data generated from genome-wide sequences. The program performs cross comparisons on the STR sequences generated using the Tandem Repeats Finder based on multiple-genome samples in a FASTA format. These comparisons generate reports listing identical, polymorphic, and different STR loci when comparing two samples.

Methods

The web site http://forensic.mc.ntu.edu.tw:9000/PSTRWeb/Default has been developed as a means to identify polymorphic STR loci within complex mass genome sequences. The program was developed to generate a series of user-friendly reports.

Results

As proof of concept for the program, four FASTA genome sequence samples of human chromosome X (AC_000155.1, CM000685.1, NC_018934.2, and CM000274.1) were obtained from GenBank and were analyzed for the presence of putative STR regions. The sequences within AC-000155.1 were used as an initial reference sequence from which there were 5443 identical and 4305 polymorphic STR loci identified using a repeat unit of 1–6 and 10 bp as the flanking sequence either side of the putative STR loci. A reliability test was used to compare five FASTA samples, which had sections of DNA sequence removed to mimic partial or fragmented DNA sequences, to determine whether pSTR Finder can efficiently and consistently find identical, polymorphic, and different STR loci.

Conclusions

From the mass of DNA sequence data, the project was found to reproducibly identify polymorphic STR loci and generate user-friendly reports detailing the number and location of these potential polymorphic loci. This freely available program was found to be a useful tool to find polymorphic STR within whole-genome sequence data in forensic genetic studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s13323-015-0027-x) contains supplementary material, which is available to authorized users.

Assembly and variation analyses of Clarias batrachus mitogenome retrieved from WGS data and its phylogenetic relationship with other catfishes

Whole genome sequencing (WGS) using next generation sequencing technologies paves the way to sequence the mitochondrial genomes with greater ease and lesser time. Here, we used the WGS data of Clarias batrachus, generated from Roche 454 and Ion Torrent sequencing platforms, to assemble the complete mitogenome using both de novo and reference based approaches. Both the methods yielded almost similar results and the best assembled mitogenome was of 16,510 bp size (GenBank Acc. No. KM259918). The mitogenome annotation resulted in 13 coding genes, 22 tRNA genes, 2 rRNA genes and one control region, and the gene order was found to be identical with other catfishes. Variation analyses between assembled and the reference (GenBank Acc. No. NC_023923) mitogenome revealed 51 variations. The phylogenetic analysis of coding DNA sequences and tRNA supports the monophyly of catfishes. Two SSRs were identified in C. batrachus mitogenome, out of which one was unique to this species. Based on the relative rate of gene evolution, protein coding mitochondrial genes were found to evolve at a much faster pace than the d-loop, which in turn are followed by the rRNAs; the tRNAs showed wide variability in the rate of sequence evolution, and on average evolve the slowest. Among the coding genes, ND2 evolves most rapidly. The variations present in the coding regions of the mitogenome and their comparative analyses with other catfish species may be useful in species conservation and management programs.

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

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

An alternative suite of universal primers for genotyping in multiplex PCR

The universal primer three-primer approach can dramatically reduce the cost when genotyping the microsatellites. One former research reported four universal primers that can be used in singleplex and multiplex genotyping. In this study, we proposed an alternative suite of universal primers with four dyes for genotyping 8–12 loci in one single run. This multiplex method was tested on Tetranychus truncatus. Published microsatellite loci of T. kanzawai, Frankliniella occidentalis and Nilaparvata lugens were modified as needed and also tested. The robustness of the method was confirmed by comparing with singleplex using multiple fluorophores and genotyping two populations of T. truncatus. This method showed lower signal strength than the singleplex three-primer system, but it was still sufficient to determine the fragment length. The cost of such a project can be reduced dramatically when many loci of different species are involved. In this way, laboratories performing population genetic analyses or studying several different species may benefit from the use of this cost-effective protocol.

Integrated shotgun sequencing and bioinformatics pipeline allows ultra-fast mitogenome recovery and confirms substantial gene rearrangements in Australian freshwater crayfishes

BACKGROUND

Although it is possible to recover the complete mitogenome directly from shotgun sequencing data, currently reported methods and pipelines are still relatively time consuming and costly. Using a sample of the Australian freshwater crayfish Engaeus lengana, we demonstrate that it is possible to achieve three-day turnaround time (four hours hands-on time) from tissue sample to NCBI-ready submission file through the integration of MiSeq sequencing platform, Nextera sample preparation protocol, MITObim assembly algorithm and MITOS annotation pipeline.

RESULTS

The complete mitochondrial genome of the parastacid freshwater crayfish, Engaeus lengana, was recovered by modest shotgun sequencing (1.2 giga bases) using the Illumina MiSeq benchtop sequencing platform. Genome assembly using the MITObim mitogenome assembler recovered the mitochondrial genome as a single contig with a 97-fold mean coverage (min. = 17; max. = 138). The mitogenome consists of 15,934 base pairs and contains the typical 37 mitochondrial genes and a non-coding AT-rich region. The genome arrangement is similar to the only other published parastacid mitogenome from the Australian genus Cherax.

CONCLUSIONS

We infer that the gene order arrangement found in Cherax destructor is common to Australian crayfish and may be a derived feature of the southern hemisphere family Parastacidae. Further, we report to our knowledge, the simplest and fastest protocol for the recovery and assembly of complete mitochondrial genomes using the MiSeq benchtop sequencer.