Comparative morphology of the electrosensory system of the epaulette shark Hemiscyllium ocellatum and brown-banded bamboo shark Chiloscyllium punctatum

We compared the electrosensory system of two benthic elasmobranchs Hemiscyllium ocellatum and Chiloscyllium punctatum. The distribution of the ampullary pores on the head was similar for both species, with a higher density of pores anteriorly and a lower density posteriorly, although C. punctatum generally possessed larger pores. Ampullary canals of the mandibular cluster were quasi-sinusoidal in H. ocellatum, a shape previously found in benthic rays only, whereas ampullary canals in C. punctatum were of a linear morphology as reported for many shark and ray species previously. The ampullae proper were of the lobular type, as occurs in most galean sharks. Chiloscyllium punctatum had six sensory chambers compared with the five per ampulla in H. ocellatum, which were generally smaller than those of C. punctatum. The sensory epithelium comprised flattened receptor cells, compared with the usual pear-shaped receptor cells encountered in other elasmobranchs and their apically nucleated supportive cells did not protrude markedly into the ampullary lumen, unlike those in benthic rays.     

A versatile Rapture (RAD‐Capture) platform for genotyping marine turtles

Advances in high‐throughput sequencing (HTS) technologies coupled with increased interdisciplinary collaboration are rapidly expanding capacity in the scope and scale of wildlife genetic studies. While existing HTS methods can be directly applied to address some evolutionary and ecological questions, certain research goals necessitate tailoring methods to specific study organisms, such as high‐throughput genotyping of the same loci that are comparable over large spatial and temporal scales. These needs are particularly common for studies of highly mobile species of conservation concern like marine turtles, where life history traits, limited financial resources and other constraints require affordable, adaptable methods for HTS genotyping to meet a variety of study goals. Here, we present a versatile marine turtle HTS targeted enrichment platform adapted from the recently developed Rapture (RAD‐Capture) method specifically designed to meet these research needs. Our results demonstrate consistent enrichment of targeted regions throughout the genome and discovery of candidate variants in all species examined for use in various conservation genetics applications. Accurate species identification confirmed the ability of our platform to genotype over 1,000 multiplexed samples and identified areas for future methodological improvement such as optimization for low initial concentration samples. Finally, analyses within green turtles supported the ability of this platform to identify informative SNPs for stock structure, population assignment and other applications over a broad geographic range of interest to management. This platform provides an additional tool for marine turtle genetic studies and broadens capacity for future large‐scale initiatives such as collaborative global marine turtle genetic databases.     

Characterization of 12 polymorphic microsatellite loci for the whitespotted bamboo shark (Chiloscyllium plagiosum Bennett)

Twelve polymorphic microsatellite loci without linkage disequilibrium were isolated and characterized from a microsatellite DNA-enriched DNA library for the whitespotted bamboo shark (Chiloscyllium plagiosum Bennett, 1830). These loci show polymorphism information content ranging from 0.384 to 0.885, allele number ranging from 4 to 12, effective allele number ranging from 1.686 to 9.438, and observed and expected heterozygosities from 0.229 to 1.000 and from 0.407 to 0.894 respectively. Four loci show strong Hardy-Weinberg deviations. We expect that these markers would be useful for population genetic and breeding studies of the whitespotted bamboo shark.     

Degenerate adaptor sequences for detecting PCR duplicates in reduced representation sequencing data improve genotype calling accuracy

RAD‐tag is a powerful tool for high‐throughput genotyping. It relies on PCR amplification of the starting material, following enzymatic digestion and sequencing adaptor ligation. Amplification introduces duplicate reads into the data, which arise from the same template molecule and are statistically nonindependent, potentially introducing errors into genotype calling. In shotgun sequencing, data duplicates are removed by filtering reads starting at the same position in the alignment. However, restriction enzymes target specific locations within the genome, causing reads to start in the same place, and making it difficult to estimate the extent of PCR duplication. Here, we introduce a slight change to the Illumina sequencing adaptor chemistry, appending a unique four‐base tag to the first index read, which allows duplicate discrimination in aligned data. This approach was validated on the Illumina MiSeq platform, using double‐digest libraries of ants (Wasmannia auropunctata) and yeast (Saccharomyces cerevisiae) with known genotypes, producing modest though statistically significant gains in the odds of calling a genotype accurately. More importantly, removing duplicates also corrected for strong sample‐to‐sample variability of genotype calling accuracy seen in the ant samples. For libraries prepared from low‐input degraded museum bird samples (Mixornis gularis), which had low complexity, having been generated from relatively few starting molecules, adaptor tags show that virtually all of the genotypes were called with inflated confidence as a result of PCR duplicates. Quantification of library complexity by adaptor tagging does not significantly increase the difficulty of the overall workflow or its cost, but corrects for differences in quality between samples and permits analysis of low‐input material.     

A MinION™‐based pipeline for fast and cost‐effective DNA barcoding

DNA barcodes are useful for species discovery and species identification, but obtaining barcodes currently requires a well‐equipped molecular laboratory and is time‐consuming, and/or expensive. We here address these issues by developing a barcoding pipeline for Oxford Nanopore MinION? and demonstrating that one flow cell can generate barcodes for ~500 specimens despite the high basecall error rates of MinION? reads. The pipeline overcomes these errors by first summarizing all reads for the same tagged amplicon as a consensus barcode. Consensus barcodes are overall mismatch‐free but retain indel errors that are concentrated in homopolymeric regions. They are addressed with an optional error correction pipeline that is based on conserved amino acid motifs from publicly available barcodes. The effectiveness of this pipeline is documented by analysing reads from three MinION? runs that represent three different stages of MinION? development. They generated data for (i) 511 specimens of a mixed Diptera sample, (ii) 575 specimens of ants and (iii) 50 specimens of Chironomidae. The run based on the latest chemistry yielded MinION? barcodes for 490 of the 511 specimens which were assessed against reference Sanger barcodes (N = 471). Overall, the MinION? barcodes have an accuracy of 99.3%–100% with the number of ambiguous bases after correction ranging from <0.01% to 1.5% depending on which correction pipeline is used. We demonstrate that it requires ~2 hr of sequencing to gather all information needed for obtaining reliable barcodes for most specimens (>90%). We estimate that up to 1,000 barcodes can be generated in one flow cell and that the cost per barcode can be 相似文献   

High‐throughput SNP discovery in the rabbit (Oryctolagus cuniculus) genome by next‐generation semiconductor‐based sequencing

The European rabbit (Oryctolagus cuniculus) is a domesticated species with one of the broadest ranges of economic and scientific applications and fields of investigation. Rabbit genome information and assembly are available (oryCun2.0), but so far few studies have investigated its variability, and massive discovery of polymorphisms has not been published yet for this species. Here, we sequenced two reduced representation libraries (RRLs) to identify single nucleotide polymorphisms (SNPs) in the rabbit genome. Genomic DNA of 10 rabbits belonging to different breeds was pooled and digested with two restriction enzymes (HaeIII and RsaI) to create two RRLs which were sequenced using the Ion Torrent Personal Genome Machine. The two RRLs produced 2 917 879 and 4 046 871 reads, for a total of 280.51 Mb (248.49 Mb with quality >20) and 417.28 Mb (360.89 Mb with quality >20) respectively of sequenced DNA. About 90% and 91% respectively of the obtained reads were mapped on the rabbit genome, covering a total of 15.82% of the oryCun2.0 genome version. The mapping and ad hoc filtering procedures allowed to reliably call 62 491 SNPs. SNPs in a few genomic regions were validated by Sanger sequencing. The Variant Effect Predictor Web tool was used to map SNPs on the current version of the rabbit genome. The obtained results will be useful for many applied and basic research programs for this species and will contribute to the development of cost‐effective solutions for high‐throughput SNP genotyping in the rabbit.     

The (in)complete organelle genome: exploring the use and nonuse of available technologies for characterizing mitochondrial and plastid chromosomes

Not long ago, scientists paid dearly in time, money and skill for every nucleotide that they sequenced. Today, DNA sequencing technologies epitomize the slogan ‘faster, easier, cheaper and more’, and in many ways, sequencing an entire genome has become routine, even for the smallest laboratory groups. This is especially true for mitochondrial and plastid genomes. Given their relatively small sizes and high copy numbers per cell, organelle DNAs are currently among the most highly sequenced kind of chromosome. But accurately characterizing an organelle genome and the information it encodes can require much more than DNA sequencing and bioinformatics analyses. Organelle genomes can be surprisingly complex and can exhibit convoluted and unconventional modes of gene expression. Unravelling this complexity can demand a wide assortment of experiments, from pulsed‐field gel electrophoresis to Southern and Northern blots to RNA analyses. Here, we show that it is exactly these types of ‘complementary’ analyses that are often lacking from contemporary organelle genome papers, particularly short ‘genome announcement’ articles. Consequently, crucial and interesting features of organelle chromosomes are going undescribed, which could ultimately lead to a poor understanding and even a misrepresentation of these genomes and the genes they express. High‐throughput sequencing and bioinformatics have made it easy to sequence and assemble entire chromosomes, but they should not be used as a substitute for or at the expense of other types of genomic characterization methods.     

SNP discovery and genotyping using restriction‐site‐associated DNA sequencing in chickens

Single nucleotide polymorphisms (SNPs) are essential to the understanding of population genetic variation and diversity. Here, we performed restriction‐site‐associated DNA sequencing (RAD‐seq) on 72 individuals from 13 Chinese indigenous and three introduced chicken breeds. A total of 620 million reads were obtained using an Illumina Hiseq2000 sequencer. An average of 75 587 SNPs were identified from each individual. Further filtering strictly validated 28 895 SNPs candidates for all populations. When compared with the NCBI dbSNP (chicken_9031), 15 404 SNPs were new discoveries. In this study, RAD‐seq was performed for the first time on chickens, implicating the remarkable effectiveness and potential applications on genetic analysis and breeding technique for whole‐genome selection in chicken and other agricultural animals.     

Mapping migration in a songbird using high‐resolution genetic markers

Neotropic migratory birds are declining across the Western Hemisphere, but conservation efforts have been hampered by the inability to assess where migrants are most limited—the breeding grounds, migratory stopover sites or wintering areas. A major challenge has been the lack of an efficient, reliable and broadly applicable method for measuring the strength of migratory connections between populations across the annual cycle. Here, we show how high‐resolution genetic markers can be used to identify genetically distinct groups of a migratory bird, the Wilson's warbler (Cardellina pusilla), at fine enough spatial scales to facilitate assessing regional drivers of demographic trends. By screening 1626 samples using 96 highly divergent single nucleotide polymorphisms selected from a large pool of candidates (~450 000), we identify novel region‐specific migratory routes and timetables of migration along the Pacific Flyway. Our results illustrate that high‐resolution genetic markers are more reliable, precise and amenable to high throughput screening than previously described intrinsic marking techniques, making them broadly applicable to large‐scale monitoring and conservation of migratory organisms.     

Genetic sex assignment in wild populations using genotyping‐by‐sequencing data: A statistical threshold approach

Establishing the sex of individuals in wild systems can be challenging and often requires genetic testing. Genotyping‐by‐sequencing (GBS) and other reduced‐representation DNA sequencing (RRS) protocols (e.g., RADseq, ddRAD) have enabled the analysis of genetic data on an unprecedented scale. Here, we present a novel approach for the discovery and statistical validation of sex‐specific loci in GBS data sets. We used GBS to genotype 166 New Zealand fur seals (NZFS, Arctocephalus forsteri) of known sex. We retained monomorphic loci as potential sex‐specific markers in the locus discovery phase. We then used (i) a sex‐specific locus threshold (SSLT) to identify significantly male‐specific loci within our data set; and (ii) a significant sex‐assignment threshold (SSAT) to confidently assign sex in silico the presence or absence of significantly male‐specific loci to individuals in our data set treated as unknowns (98.9% accuracy for females; 95.8% for males, estimated via cross‐validation). Furthermore, we assigned sex to 86 individuals of true unknown sex using our SSAT and assessed the effect of SSLT adjustments on these assignments. From 90 verified sex‐specific loci, we developed a panel of three sex‐specific PCR primers that we used to ascertain sex independently of our GBS data, which we show amplify reliably in at least two other pinniped species. Using monomorphic loci normally discarded from large SNP data sets is an effective way to identify robust sex‐linked markers for nonmodel species. Our novel pipeline can be used to identify and statistically validate monomorphic and polymorphic sex‐specific markers across a range of species and RRS data sets.     

Long‐term persistence of wildlife populations in a pastoral area

Facilitating coexistence between people and wildlife is a major conservation challenge in East Africa. Some conservation models aim to balance the needs of people and wildlife, but the effectiveness of these models is rarely assessed. Using a case‐study approach, we assessed the ecological performance of a pastoral area in northern Tanzania (Manyara Ranch) and established a long‐term wildlife population monitoring program (carried out intermittently from 2003 to 2008 and regularly from 2011 to 2019) embedded in a distance sampling framework. By comparing density estimates of the road transect‐based long‐term monitoring to estimates derived from systematically distributed transects, we found that the bias associated with nonrandom placement of transects was nonsignificant. Overall, cattle and sheep and goat reached the greatest densities and several wildlife species occurred at densities similar (zebra, wildebeest, waterbuck, Kirk's dik‐dik) or possibly even greater (giraffe, eland, lesser kudu, Grant's gazelle, Thomson's gazelle) than in adjacent national parks in the same ecosystem. Generalized linear mixed models suggested that most wildlife species (8 out of 14) reached greatest densities during the dry season, that wildlife population densities either remained constant or increased over the 17‐year period, and that herbivorous livestock species remained constant, while domestic dog population decreased over time. Cross‐species correlations did not provide evidence for interference competition between grazing or mixed livestock species and wildlife species but indicate possible negative relationships between domestic dog and warthog populations. Overall, wildlife and livestock populations in Manyara Ranch appear to coexist over the 17‐year span. Most likely, this is facilitated by existing connectivity to adjacent protected areas, effective anti‐poaching efforts, spatio‐temporal grazing restrictions, favorable environmental conditions of the ranch, and spatial heterogeneity of surface water and habitats. This long‐term case study illustrates the potential of rangelands to simultaneously support wildlife conservation and human livelihood goals if livestock grazing is restricted in space, time, and numbers.     

DNA barcodes from century‐old type specimens using next‐generation sequencing

Type specimens have high scientific importance because they provide the only certain connection between the application of a Linnean name and a physical specimen. Many other individuals may have been identified as a particular species, but their linkage to the taxon concept is inferential. Because type specimens are often more than a century old and have experienced conditions unfavourable for DNA preservation, success in sequence recovery has been uncertain. This study addresses this challenge by employing next‐generation sequencing (NGS) to recover sequences for the barcode region of the cytochrome c oxidase 1 gene from small amounts of template DNA. DNA quality was first screened in more than 1800 century‐old type specimens of Lepidoptera by attempting to recover 164‐bp and 94‐bp reads via Sanger sequencing. This analysis permitted the assignment of each specimen to one of three DNA quality categories – high (164‐bp sequence), medium (94‐bp sequence) or low (no sequence). Ten specimens from each category were subsequently analysed via a PCR‐based NGS protocol requiring very little template DNA. It recovered sequence information from all specimens with average read lengths ranging from 458 bp to 610 bp for the three DNA categories. By sequencing ten specimens in each NGS run, costs were similar to Sanger analysis. Future increases in the number of specimens processed in each run promise substantial reductions in cost, making it possible to anticipate a future where barcode sequences are available from most type specimens.     

Age and growth of the tropical oviparous shark,Chiloscyllium punctatum in Indonesian waters

The brown-banded bamboo shark, Chiloscyllium punctatum, is the most common shark caught in coastal commercial fisheries throughout Southeast Asia, yet there is a lack of the life-history information necessary for reliable stock assessments. The authors estimated growth rates and age at maturity using analysis of growth bands in vertebral centra. They trialled four different techniques to enhance the visibility and improve identification of the putative annual growth bands necessary for age estimation. The authors found that the burn method on whole vertebral centra provided the most readable and consistent results for age analysis. The logistic model was chosen as the best-fit growth model for age estimation of 330 individual C. punctatum from Indonesia. Several age verification methods, including marginal increment ratio and length-frequency analysis, were performed with the support of age validation through the use of calcein-labelled vertebrae from two sharks maintained in captivity. This study found that C. punctatum from Indonesian waters is a fast-growing species that can grow up to 18 cm year⁻¹, reach an estimated maximum total length of 1 m, mature at c. 6.5 years and live for up to 14 years.     

Recent progress and challenges in population genetics of polyploid organisms: an overview of current state‐of‐the‐art molecular and statistical tools

Despite the importance of polyploidy and the increasing availability of new genomic data, there remain important gaps in our knowledge of polyploid population genetics. These gaps arise from the complex nature of polyploid data (e.g. multiple alleles and loci, mixed inheritance patterns, association between ploidy and mating system variation). Furthermore, many of the standard tools for population genetics that have been developed for diploids are often not feasible for polyploids. This review aims to provide an overview of the state‐of‐the‐art in polyploid population genetics and to identify the main areas where further development of molecular techniques and statistical theory is required. We review commonly used molecular tools (amplified fragment length polymorphism, microsatellites, Sanger sequencing, next‐generation sequencing and derived technologies) and their challenges associated with their use in polyploid populations: that is, allele dosage determination, null alleles, difficulty of distinguishing orthologues from paralogues and copy number variation. In addition, we review the approaches that have been used for population genetic analysis in polyploids and their specific problems. These problems are in most cases directly associated with dosage uncertainty and the problem of inferring allele frequencies and assumptions regarding inheritance. This leads us to conclude that for advancing the field of polyploid population genetics, most priority should be given to development of new molecular approaches that allow efficient dosage determination, and to further development of analytical approaches to circumvent dosage uncertainty and to accommodate ‘flexible’ modes of inheritance. In addition, there is a need for more simulation‐based studies that test what kinds of biases could result from both existing and novel approaches.     

How quantitative is metabarcoding: A meta‐analytical approach

Metabarcoding has been used in a range of ecological applications such as taxonomic assignment, dietary analysis and the analysis of environmental DNA. However, after a decade of use in these applications there is little consensus on the extent to which proportions of reads generated corresponds to the original proportions of species in a community. To quantify our current understanding, we conducted a structured review and meta‐analysis. The analysis suggests that a weak quantitative relationship may exist between the biomass and sequences produced (slope = 0.52 ± 0.34, p < 0.01), albeit with a large degree of uncertainty. None of the tested moderators, sequencing platform type, the number of species used in a trial or the source of DNA, were able to explain the variance. Our current understanding of the factors affecting the quantitative performance of metabarcoding is still limited: additional research is required before metabarcoding can be confidently utilized for quantitative applications. Until then, we advocate the inclusion of mock communities when metabarcoding as this facilitates direct assessment of the quantitative ability of any given study.     

Methyl CpG Binding Domain Ultra‐Sequencing: a novel method for identifying inter‐individual and cell‐type‐specific variation in DNA methylation

Experience‐dependent changes in DNA methylation can exert profound effects on neuronal function and behaviour. A single learning event can induce a variety of DNA modifications within the neuronal genome, some of which may be common to all individuals experiencing the event, whereas others may occur in a subset of individuals. Variations in experience‐induced DNA methylation may subsequently confer increased vulnerability or resilience to the development of neuropsychiatric disorders. However, the detection of experience‐dependent changes in DNA methylation in the brain has been hindered by the interrogation of heterogeneous cell populations, regional differences in epigenetic states and the use of pooled tissue obtained from multiple individuals. Methyl CpG Binding Domain Ultra‐Sequencing (MBD Ultra‐Seq) overcomes current limitations on genome‐wide epigenetic profiling by incorporating fluorescence‐activated cell sorting and sample‐specific barcoding to examine cell‐type‐specific CpG methylation in discrete brain regions of individuals. We demonstrate the value of this method by characterizing differences in 5‐methylcytosine (5mC) in neurons and non‐neurons of the ventromedial prefrontal cortex of individual adult C57BL/6 mice, using as little as 50 ng of genomic DNA per sample. We find that the neuronal methylome is characterized by greater CpG methylation as well as the enrichment of 5mC within intergenic loci. In conclusion, MBD Ultra‐Seq is a robust method for detecting DNA methylation in neurons derived from discrete brain regions of individual animals. This protocol will facilitate the detection of experience‐dependent changes in DNA methylation in a variety of behavioural paradigms and help identify aberrant experience‐induced DNA methylation that may underlie risk and resiliency to neuropsychiatric disease.