Rapid in situ identification of biological specimens via DNA amplicon sequencing using miniaturized laboratory equipment.

In many parts of the world, human-mediated environmental change is depleting biodiversity faster than it can be characterized, while invasive species cause agricultural damage, threaten human health and disrupt native habitats. Consequently, the application of effective approaches for rapid surveillance and identification of biological specimens is increasingly important to inform conservation and biosurveillance efforts. Taxonomic assignments have been greatly advanced using sequence-based applications, such as DNA barcoding, a diagnostic technique that utilizes PCR and DNA sequence analysis of standardized genetic regions. However, in many biodiversity hotspots, endeavors are often hindered by a lack of laboratory infrastructure, funding for biodiversity research and restrictions on the transport of biological samples. A promising development is the advent of low-cost, miniaturized scientific equipment. Such tools can be assembled into functional laboratories to carry out genetic analyses in situ, at local institutions, field stations or classrooms. Here, we outline the steps required to perform amplicon sequencing applications, from DNA isolation to nanopore sequencing and downstream data analysis, all of which can be conducted outside of a conventional laboratory environment using miniaturized scientific equipment, without reliance on Internet connectivity. Depending on sample type, the protocol (from DNA extraction to full bioinformatic analyses) can be completed within 10 h, and with appropriate quality controls can be used for diagnostic identification of samples independent of core genomic facilities that are required for alternative methods.

Developmental validation of Oxford Nanopore Technology MinION sequence data and the NGSpeciesID bioinformatic pipeline for forensic genetic species identification.

Species identification of non-human biological evidence through DNA nucleotide sequencing is routinely used for forensic genetic analysis to support law enforcement. The gold standard for forensic genetics is conventional Sanger sequencing; however, this is gradually being replaced by high-throughput sequencing (HTS) approaches which can generate millions of individual reads in a single experiment. HTS sequencing, which now dominates molecular biology research, has already been demonstrated for use in a number of forensic genetic analysis applications, including species identification. However, the generation of HTS data to date requires expensive equipment and is cost-effective only when large numbers of samples are analysed simultaneously. The Oxford Nanopore Technologies (ONT) MinION™ is an affordable and small footprint DNA sequencing device with the potential to quickly deliver reliable and cost effective data. However, there has been no formal validation of forensic species identification using high-throughput (deep read) sequence data from the MinION making it currently impractical for many wildlife forensic end-users. Here, we present a MinION deep read sequence data validation study for species identification. First, we tested whether the clustering-based bioinformatics pipeline NGSpeciesID can be used to generate an accurate consensus sequence for species identification. Second, we systematically evaluated the read variation distribution around the generated consensus sequences to understand what confidence we have in the accuracy of the resulting consensus sequence and to determine how to interpret individual sample results. Finally, we investigated the impact of differences between the MinION consensus and Sanger control sequences on correct species identification to understand the ability and accuracy of the MinION consensus sequence to differentiate the true species from the next most similar species. This validation study establishes that ONT MinION sequence data used in conjunction with the NGSpeciesID pipeline can produce consensus DNA sequences of sufficient accuracy for forensic genetic species identification.

BURMESE ROOFED TURTLE (BATAGUR TRIVITTATA) DISEASE SCREENING IN MYANMAR.

The Burmese roofed turtle (Batagur trivittata), a critically endangered freshwater turtle, is endemic to Myanmar. Once thought to be extinct, remnant wild populations were discovered in 2001 and limited captive individuals identified in pagoda ponds or confiscated from fishers in Myanmar. These and their offspring are maintained in five facilities in Myanmar and form the basis of a conservation program (habitat protection, captive breeding, nest protection, egg collection, head-starting, and release). Prerelease health screenings were performed in 2014 and 2018 at Yadanabon Zoological Gardens, a head-starting facility in Limpha Village, and Lawkanandar Wildlife Park. One hundred forty-three turtles were assessed (37 male, 50 female, 56 juveniles [too young to determine sex]; two females were assessed in both years), age range of 1 to 12 y (one unknown age adult founder), and body mass range of 0.111 to 32.72 kg. Health evaluations both years included physical examination and combined choanal/cloacal swab samples for polymerase chain reaction testing of the potential chelonian pathogens intranuclear coccidia, Mycoplasma, Herpesvirus, Ranavirus, and Adenovirus (not all tests performed each year). In 2018, cloacal swabs from 30 and 20 turtles at the Yadanabon Zoological Gardens and Lawkanandar Wildlife Park, respectively, were cultured for Salmonella. All turtles were assessed as healthy based on normal physical examination findings, and all had negative test results. Prerelease health screening, such as performed in this study, is an important component of release, reintroduction, and translocation projects to prevent introduction of novel pathogens into naïve wild populations.

MinION-Based DNA Barcoding of Preserved and Non-Invasively Collected Wildlife Samples.

The ability to sequence a variety of wildlife samples with portable, field-friendly equipment will have significant impacts on wildlife conservation and health applications. However, the only currently available field-friendly DNA sequencer, the MinION by Oxford Nanopore Technologies, has a high error rate compared to standard laboratory-based sequencing platforms and has not been systematically validated for DNA barcoding accuracy for preserved and non-invasively collected tissue samples. We tested whether various wildlife sample types, field-friendly methods, and our clustering-based bioinformatics pipeline, SAIGA, can be used to generate consistent and accurate consensus sequences for species identification. Here, we systematically evaluate variation in cytochrome b sequences amplified from scat, hair, feather, fresh frozen liver, and formalin-fixed paraffin-embedded (FFPE) liver. Each sample was processed by three DNA extraction protocols. For all sample types tested, the MinION consensus sequences matched the Sanger references with 99.29%-100% sequence similarity, even for samples that were difficult to amplify, such as scat and FFPE tissue extracted with Chelex resin. Sequencing errors occurred primarily in homopolymer regions, as identified in previous MinION studies. We demonstrate that it is possible to generate accurate DNA barcode sequences from preserved and non-invasively collected wildlife samples using portable MinION sequencing, creating more opportunities to apply portable sequencing technology for species identification.

Last exit before the brink: Conservation genomics of the Cambodian population of the critically endangered southern river terrapin.

The southern river terrapin, Batagur affinis is one of the world's 25 most endangered freshwater turtle species. The major portion of the global population is currently found in peninsular Malaysia, with the only remnant Indochinese population in southern Cambodia. For more than a decade, wild nests in this remnant Cambodian population have been fenced and hatchlings reared in captivity. Here we amplified 10 microsatellite markers from all 136 captive individuals, obtained 2,658 presumably unlinked and neutral single nucleotide polymorphisms from 72 samples with ddRAD-seq, and amplified 784 bp of mtDNA from 50 samples. Our results reveal that the last Indochinese population comprised only four kinship groups as of 2012, with all offspring sired from <10 individuals in the wild. We demonstrate an obvious decrease in genetic contributions of breeders in the wild from 2006-2012 and identify high-value breeders instrumental for ex-situ management of the contemporary genetic stock of the species.

Re-programming of C. elegans male epidermal precursor fates by Wnt, Hox, and LIN-12/Notch activities.

In Caenorhabditiselegans males, different subsets of ventral epidermal precursor (Pn.p) cells adopt distinct fates in a position-specific manner: three posterior cells, P(9-11).p, comprise the hook sensillum competence group (HCG) with three potential fates (1 degrees , 2 degrees , or 3 degrees ), while eight anterior cells, P(1-8).p, fuse with the hyp7 epidermal syncytium. Here we show that activation of the canonical BAR-1 beta-catenin pathway of Wnt signaling alters the competence of P(3-8).p and specifies ectopic HCG-like fates. This fate transformation requires the Hox gene mab-5. In addition, misexpression of mab-5 in P(1-8).p is sufficient to establish HCG competence among these cells, as well as to generate ectopic HCG fates in combination with LIN-12 or EGF signaling. While increased Wnt signaling induces predominantly 1 degrees HCG fates, increased LIN-12 or EGF signaling in combination with MAB-5 overexpression promotes 2 degrees HCG fates in anterior Pn.p cells, suggesting distinctive functions of Wnt, LIN-12, and EGF signaling in specification of HCG fates. Lastly, wild-type mab-5 function is necessary for normal P(9-11).p fate specification, indicating that regulation of ectopic HCG fate formation revealed in anterior Pn.p cells reflect mechanisms of pattern formation during normal hook development.

Wnt and EGF pathways act together to induce C. elegans male hook development.

Comparative studies of vulva development between Caenorhabditis elegans and other nematode species have provided some insight into the evolution of patterning networks. However, molecular genetic details are available only in C. elegans and Pristionchus pacificus. To extend our knowledge on the evolution of patterning networks, we studied the C. elegans male hook competence group (HCG), an equivalence group that has similar developmental origins to the vulval precursor cells (VPCs), which generate the vulva in the hermaphrodite. Similar to VPC fate specification, each HCG cell adopts one of three fates (1 degree, 2 degrees, 3 degrees), and 2 degrees HCG fate specification is mediated by LIN-12/Notch. We show that 2 degrees HCG specification depends on the presence of a cell with the 1 degree fate. We also provide evidence that Wnt signaling via the Frizzled-like Wnt receptor LIN-17 acts to specify the 1 degree and 2 degrees HCG fate. A requirement for EGF signaling during 1 degree fate specification is seen only when LIN-17 activity is compromised. In addition, activation of the EGF pathway decreases dependence on LIN-17 and causes ectopic hook development. Our results suggest that WNT plays a more significant role than EGF signaling in specifying HCG fates, whereas in VPC specification EGF signaling is the major inductive signal. Nonetheless, the overall logic is similar in the VPCs and the HCG: EGF and/or WNT induce a 1 degree lineage, and LIN-12/NOTCH induces a 2 degrees lineage. Wnt signaling is also required for execution of the 1 degree and 2 degrees HCG lineages. lin-17 and bar-1/beta-catenin are preferentially expressed in the presumptive 1 degree cell P11.p. The dynamic subcellular localization of BAR-1-GFP in P11.p is concordant with the timing of HCG fate determination.

The roles of EGF and Wnt signaling during patterning of the C. elegans Bgamma/delta Equivalence Group.

BACKGROUND: During development, different signaling pathways interact to specify cell fate by regulating transcription factors necessary for fate specification and morphogenesis. In Caenorhabditis elegans, the EGF-Ras and Wnt signaling pathways have been shown to interact to specify cell fate in three equivalence groups: the vulval precursor cells (VPCs), the hook competence group (HCG) and P11/12. In the VPCs, HCG and P11/12 pair, EGF and Wnt signaling positively regulate different Hox genes, each of which also functions during fate specification. In the male, EGF-Ras signaling is required to specify the Bgamma fate within the Bgamma/delta equivalence pair, while Notch signaling is required for Bdelta fate specification. In addition, TGF-beta signaling by dbl-1/dpp controls ceh-13/labial/Hox1 expression in Bgamma. RESULTS: We show that EGF-Ras signaling is required for Bgamma expression of ceh-13/labial/Hox1. The transcription factors lin-1/ETS and lin-31/Forkhead, functioning downstream of the EGF pathway, as well as sur-2/MED23 (a component of the Mediator complex) also control ceh-13 expression in Bgamma. In addition, our results indicate that lin-44/Wnt, mom-2/Wnt and lin-17/Fz are necessary to maintain the division of Bgamma along a longitudinal axis. We also show that dbl-1/dpp acts either in parallel or downstream of EGF pathway to regulate ceh-13/Hox1 expression in Bgamma. Lastly, we find that a dbl-1/dpp null mutation did not cause any vulval or P12 defects and did not enhance vulval and P12 defects of reduction-of-function mutations of components of the EGF pathway. CONCLUSIONS: ceh-13/labial/Hox1 expression in Bgamma is regulated by the EGF pathway and downstream factors lin-1/ETS lin-31/Forkhead and sur-2/MED23. Wnt signaling is required for proper Bgamma division, perhaps to orient the Bgamma mitotic spindle. Our results suggest that dbl-1/dpp is not required for VPC and P12 specification, highlighting another difference among these EGF-dependent equivalence groups.

Systems level circuit model of C. elegans undulatory locomotion: mathematical modeling and molecular genetics.

To establish the relationship between locomotory behavior and dynamics of neural circuits in the nematode C. elegans we combined molecular and theoretical approaches. In particular, we quantitatively analyzed the motion of C. elegans with defective synaptic GABA and acetylcholine transmission, defective muscle calcium signaling, and defective muscles and cuticle structures, and compared the data with our systems level circuit model. The major experimental findings are: (1) anterior-to-posterior gradients of body bending flex for almost all strains both for forward and backward motion, and for neuronal mutants, also analogous weak gradients of undulatory frequency, (2) existence of some form of neuromuscular (stretch receptor) feedback, (3) invariance of neuromuscular wavelength, (4) biphasic dependence of frequency on synaptic signaling, and (5) decrease of frequency with increase of the muscle time constant. Based on (1) we hypothesize that the Central Pattern Generator (CPG) is located in the head both for forward and backward motion. Points (1) and (2) are the starting assumptions for our theoretical model, whose dynamical patterns are qualitatively insensitive to the details of the CPG design if stretch receptor feedback is sufficiently strong and slow. The model reveals that stretch receptor coupling in the body wall is critical for generation of the neuromuscular wave. Our model agrees with our behavioral data (3), (4), and (5), and with other pertinent published data, e.g., that frequency is an increasing function of muscle gap-junction coupling.

Conservation rules, their breakdown, and optimality in Caenorhabditis sinusoidal locomotion.

Undulatory locomotion is common to nematodes as well as to limbless vertebrates, but its control is not understood in spite of the identification of hundred of genes involved in Caenorhabditis elegans locomotion. To reveal the mechanisms of nematode undulatory locomotion, we quantitatively analysed the movement of C. elegans with genetic perturbations to neurons, muscles, and skeleton (cuticle). We also compared locomotion of different Caenorhabditis species. We constructed a theoretical model that combines mechanics and biophysics, and that is constrained by the observations of propulsion and muscular velocities, as well as wavelength and amplitude of undulations. We find that normalized wavelength is a conserved quantity among wild-type C. elegans individuals, across mutants, and across different species. The velocity of forward propulsion scales linearly with the velocity of the muscular wave and the corresponding slope is also a conserved quantity and almost optimal; the exceptions are in some mutants affecting cuticle structure. In theoretical terms, the optimality of the slope is equivalent to the exact balance between muscular and visco-elastic body reaction bending moments. We find that the amplitude and frequency of undulations are inversely correlated and provide a theoretical explanation for this fact. These experimental results are valid both for young adults and for all larval stages of wild-type C. elegans. In particular, during development, the amplitude scales linearly with the wavelength, consistent with our theory. We also investigated the influence of substrate firmness on motion parameters, and found that it does not affect the above invariants. In general, our biomechanical model can explain the observed robustness of the mechanisms controlling nematode undulatory locomotion.