Margaret E. Hunter

Validation of eDNA surveillance sensitivity for detection of Asian carps in controlled and field experiments.

In many North American rivers, populations of multiple species of non-native cyprinid fishes are present, including black carp (Mylpharyngodon piceus), grass carp (Ctenopharyngodon idella), bighead carp (Hypophthalmichthys nobilis), silver carp (Hypophthalmichthys molitrix), common carp (Cyprinus carpio), and goldfish (Carassius auratus). All six of these species are found in the Mississippi River basin and tracking their invasion has proven difficult, particularly where abundance is low. Knowledge of the location of the invasion front is valuable to natural resource managers because future ecological and economic damages can be most effectively prevented when populations are low. To test the accuracy of environmental DNA (eDNA) as an early indicator of species occurrence and relative abundance, we applied eDNA technology to the six non-native cyprinid species putatively present in a 2.6 river mile stretch of the Chicago (IL, USA) canal system that was subsequently treated with piscicide. The proportion of water samples yielding positive detections increased with relative abundance of the six species, as indicated by the number of carcasses recovered after poisoning. New markers for black carp, grass carp, and a common carp/goldfish are reported and details of the marker testing to ensure specificity are provided.

Role of a polymorphism in a Hox/Pax-responsive enhancer in the evolution of the vertebrate spine

Patterning of the vertebrate skeleton requires the coordinated activity of Hox genes. In particular, Hox10 proteins are essential to set the transition from thoracic to lumbar vertebrae because of their rib-repressing activity. In snakes, however, the thoracic region extends well into Hox10-expressing areas of the embryo, suggesting that these proteins are unable to block rib formation. Here, we show that this is not a result of the loss of rib-repressing properties by the snake proteins, but rather to a single base pair change in a Hox/Paired box (Pax)-responsive enhancer, which prevents the binding of Hox proteins. This polymorphism is also found in Paenungulata, such as elephants and manatees, which have extended rib cages. In vivo, this modified enhancer failed to respond to Hox10 activity, supporting its role in the extension of rib cages. In contrast, the enhancer could still interact with Hoxb6 and Pax3 to promote rib formation. These results suggest that a polymorphism in the Hox/Pax-responsive enhancer may have played a role in the evolution of the vertebrate spine by differently modulating its response to rib-suppressing and rib-promoting Hox proteins.

Environmental DNA (eDNA) Sampling Improves Occurrence and Detection Estimates of Invasive Burmese Pythons

Environmental DNA (eDNA) methods are used to detect DNA that is shed into the aquatic environment by cryptic or low density species. Applied in eDNA studies, occupancy models can be used to estimate occurrence and detection probabilities and thereby account for imperfect detection. However, occupancy terminology has been applied inconsistently in eDNA studies, and many have calculated occurrence probabilities while not considering the effects of imperfect detection. Low detection of invasive giant constrictors using visual surveys and traps has hampered the estimation of occupancy and detection estimates needed for population management in southern Florida, USA. Giant constrictor snakes pose a threat to native species and the ecological restoration of the Florida Everglades. To assist with detection, we developed species-specific eDNA assays using quantitative PCR (qPCR) for the Burmese python (Python molurus bivittatus), Northern African python (P. sebae), boa constrictor (Boa constrictor), and the green (Eunectes murinus) and yellow anaconda (E. notaeus). Burmese pythons, Northern African pythons, and boa constrictors are established and reproducing, while the green and yellow anaconda have the potential to become established. We validated the python and boa constrictor assays using laboratory trials and tested all species in 21 field locations distributed in eight southern Florida regions. Burmese python eDNA was detected in 37 of 63 field sampling events; however, the other species were not detected. Although eDNA was heterogeneously distributed in the environment, occupancy models were able to provide the first estimates of detection probabilities, which were greater than 91%. Burmese python eDNA was detected along the leading northern edge of the known population boundary. The development of informative detection tools and eDNA occupancy models can improve conservation efforts in southern Florida and support more extensive studies of invasive constrictors. Generic sampling design and terminology are proposed to standardize and clarify interpretations of eDNA-based occupancy models.

Marsh rabbit mortalities tie pythons to the precipitous decline of mammals in the Everglades

To address the ongoing debate over the impact of invasive species on native terrestrial wildlife, we conducted a large-scale experiment to test the hypothesis that invasive Burmese pythons (Python molurus bivittatus) were a cause of the precipitous decline of mammals in Everglades National Park (ENP). Evidence linking pythons to mammal declines has been indirect and there are reasons to question whether pythons, or any predator, could have caused the precipitous declines seen across a range of mammalian functional groups. Experimentally manipulating marsh rabbits, we found that pythons accounted for 77% of rabbit mortalities within 11 months of their translocation to ENP and that python predation appeared to preclude the persistence of rabbit populations in ENP. On control sites, outside of the park, no rabbits were killed by pythons and 71% of attributable marsh rabbit mortalities were classified as mammal predations. Burmese pythons pose a serious threat to the faunal communities and ecological functioning of the Greater Everglades Ecosystem, which will probably spread as python populations expand their range.

Low genetic diversity and minimal population substructure in the endangered Florida manatee: implications for conservation

Abstract Species of management concern that have been affected by human activities typically are characterized by low genetic diversity, which can adversely affect their ability to adapt to environmental changes. We used 18 microsatellite markers to genotype 362 Florida manatees (Trichechus manatus latirostris), and investigated genetic diversity, population structure, and estimated genetically effective population size (Ne). The observed and expected heterozygosity and average number of alleles were 0.455 ± 0.04, 0.479 ± 0.04, and 4.77 ± 0.51, respectively. All measures of Florida manatee genetic diversity were less than averages reported for placental mammals, including fragmented or nonideal populations. Overall estimates of differentiation were low, though significantly greater than zero, and analysis of molecular variance revealed that over 95% of the total variance was among individuals within predefined management units or among individuals along the coastal subpopulations, with only minor portions of variance explained by between group variance. Although genetic issues, as inferred by neutral genetic markers, appear not to be critical at present, the Florida manatee continues to face demographic challenges due to anthropogenic activities and stochastic factors such as red tides, oil spills, and disease outbreaks; these can further reduce genetic diversity of the manatee population.

Detection limits of quantitative and digital PCR assays and their influence in presence‐absence surveys of environmental DNA

A set of universal guidelines is needed to determine the limit of detection (LOD) in PCR‐based analyses of low‐concentration DNA. In particular, environmental DNA (eDNA) studies require sensitive and reliable methods to detect rare and cryptic species through shed genetic material in environmental samples. Current strategies for assessing detection limits of eDNA are either too stringent or subjective, possibly resulting in biased estimates of species’ presence. Here, a conservative LOD analysis grounded in analytical chemistry is proposed to correct for overestimated DNA concentrations predominantly caused by the concentration plateau, a nonlinear relationship between expected and measured DNA concentrations. We have used statistical criteria to establish formal mathematical models for both quantitative and droplet digital PCR. To assess the method, a new Grass Carp (Ctenopharyngodon idella) TaqMan assay was developed and tested on both PCR platforms using eDNA in water samples. The LOD adjustment reduced Grass Carp occupancy and detection estimates while increasing uncertainty—indicating that caution needs to be applied to eDNA data without LOD correction. Compared to quantitative PCR, digital PCR had higher occurrence estimates due to increased sensitivity and dilution of inhibitors at low concentrations. Without accurate LOD correction, species occurrence and detection probabilities based on eDNA estimates are prone to a source of bias that cannot be reduced by an increase in sample size or PCR replicates. Other applications also could benefit from a standardized LOD such as GMO food analysis and forensic and clinical diagnostics.

Rapid microsatellite marker development using next generation pyrosequencing to inform invasive Burmese python -- Python molurus bivittatus -- management

Invasive species represent an increasing threat to native ecosystems, harming indigenous taxa through predation, habitat modification, cross-species hybridization and alteration of ecosystem processes. Additionally, high economic costs are associated with environmental damage, restoration and control measures. The Burmese python, Python molurus bivittatus, is one of the most notable invasive species in the US, due to the threat it poses to imperiled species and the Greater Everglades ecosystem. To address population structure and relatedness, next generation sequencing was used to rapidly produce species-specific microsatellite loci. The Roche 454 GS-FLX Titanium platform provided 6616 di-, tri- and tetra-nucleotide repeats in 117,516 sequences. Using stringent criteria, 24 of 26 selected tri- and tetra-nucleotide loci were polymerase chain reaction (PCR) amplified and 18 were polymorphic. An additional six cross-species loci were amplified, and the resulting 24 loci were incorporated into eight PCR multiplexes. Multi-locus genotypes yielded an average of 61% (39%–77%) heterozygosity and 3.7 (2–6) alleles per locus. Population-level studies using the developed microsatellites will track the invasion front and monitor population-suppression dynamics. Additionally, cross-species amplification was detected in the invasive Ball, P. regius, and Northern African python, P. sebae. These markers can be used to address the hybridization potential of Burmese pythons and the larger, more aggressive P. sebae.

Statistical Models for the Analysis and Design of Digital Polymerase Chain Reaction (dPCR) Experiments.

Statistical methods for the analysis and design of experiments using digital PCR (dPCR) have received only limited attention and have been misused in many instances. To address this issue and to provide a more general approach to the analysis of dPCR data, we describe a class of statistical models for the analysis and design of experiments that require quantification of nucleic acids. These models are mathematically equivalent to generalized linear models of binomial responses that include a complementary, log-log link function and an offset that is dependent on the dPCR partition volume. These models are both versatile and easy to fit using conventional statistical software. Covariates can be used to specify different sources of variation in nucleic acid concentration, and a models parameters can be used to quantify the effects of these covariates. For purposes of illustration, we analyzed dPCR data from different types of experiments, including serial dilution, evaluation of copy number variation, and quantification of gene expression. We also showed how these models can be used to help design dPCR experiments, as in selection of sample sizes needed to achieve desired levels of precision in estimates of nucleic acid concentration or to detect differences in concentration among treatments with prescribed levels of statistical power.

The Florida manatee (Trichechus manatus latirostris) immunoglobulin heavy chain suggests the importance of clan III variable segments in repertoire diversity

Abstract Manatees are a vulnerable, charismatic sentinel species from the evolutionarily divergent Afrotheria. Manatee health and resistance to infectious disease is of great concern to conservation groups, but little is known about their immune system. To develop manatee‐specific tools for monitoring health, we first must have a general knowledge of how the immunoglobulin heavy (IgH) chain locus is organized and transcriptionally expressed. Using the genomic scaffolds of the Florida manatee (Trichechus manatus latirostris), we characterized the potential IgH segmental diversity and constant region isotypic diversity and performed the first Afrotherian repertoire analysis. The Florida manatee has low V(D)J combinatorial diversity (3744 potential combinations) and few constant region isotypes. They also lack clan III V segments, which may have caused reduced VH segment numbers. However, we found productive somatic hypermutation concentrated in the complementarity determining regions. In conclusion, manatees have limited IGHV clan and combinatorial diversity. This suggests that clan III V segments are essential for maintaining IgH locus diversity. HighlightsT. m. latirostris has only four functional IgH constant region isotypes.There are few functional VH, no clan III VH, and low combinatorial diversity.The absence of clan III VH correlates with low VH numbers in three other species.This is the first IgH repertoire analysis within Afrotheria.Nonsynonymous somatic hypermutation is higher in the CDR regions.

DEVELOPMENT AND VALIDATION OF QUANTITATIVE PCR ASSAYS TO MEASURE CYTOKINE TRANSCRIPT LEVELS IN THE FLORIDA MANATEE (TRICHECHUS MANATUS LATIROSTRIS)

Abstract:  Cytokines have important roles in the mammalian response to viral and bacterial infections, trauma, and wound healing. Because of early cytokine production after physiologic stresses, the regulation of messenger RNA (mRNA) transcripts can be used to assess immunologic responses before changes in protein production. To detect and assess early immune changes in endangered Florida manatees (Trichechus manatus latirostris), we developed and validated a panel of quantitative PCR assays to measure mRNA transcription levels for the cytokines interferon (IFN)-γ; interleukin (IL)-2, -6, and -10; tumor necrosis factor-α; and the housekeeping genes glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and β-actin (reference genes). Assays were successfully validated using blood samples from free-ranging, apparently healthy manatees from the east and west coasts of central Florida, US. No cytokine or housekeeping gene transcription levels were significantly different among age classes or sexes. However, the transcription levels for GAPDH, IL-2, IL-6, and IFN-γ were significantly higher (P<0.05) in manatees from the east coast of Florida than they were from those from the west coast. We found IL-10 and β-actin to be consistent between sites and identified β-actin as a good candidate for use as a reference gene in future studies. Our assays can aid in the investigation of manatee immune response to physical trauma and novel or ongoing environmental stressors.