Using global remote camera data of a solitary species complex to evaluate the drivers of group formation.

The social system of animals involves a complex interplay between physiology, natural history, and the environment. Long relied upon discrete categorizations of "social" and "solitary" inhibit our capacity to understand species and their interactions with the world around them. Here, we use a globally distributed camera trapping dataset to test the drivers of aggregating into groups in a species complex (martens and relatives, family Mustelidae, Order Carnivora) assumed to be obligately solitary. We use a simple quantification, the probability of being detected in a group, that was applied across our globally derived camera trap dataset. Using a series of binomial generalized mixed-effects models applied to a dataset of 16,483 independent detections across 17 countries on four continents we test explicit hypotheses about potential drivers of group formation. We observe a wide range of probabilities of being detected in groups within the solitary model system, with the probability of aggregating in groups varying by more than an order of magnitude. We demonstrate that a species' context-dependent proclivity toward aggregating in groups is underpinned by a range of resource-related factors, primarily the distribution of resources, with increasing patchiness of resources facilitating group formation, as well as interactions between environmental conditions (resource constancy/winter severity) and physiology (energy storage capabilities). The wide variation in propensities to aggregate with conspecifics observed here highlights how continued failure to recognize complexities in the social behaviors of apparently solitary species limits our understanding not only of the individual species but also the causes and consequences of group formation.

The genetic basis of conspicuous coloration in the Guadeloupean anole: Evolution by sexual and ecological selection.

Understanding how natural selection acts on the genome and contributes to the process of speciation is a primary aim of the study of evolution. Here we used natural variation in two subspecies of the Guadeloupean anole (Anolis marmoratus ssp.), from the island of Guadeloupe in the Lesser Antilles, to explore the genomic basis of adaptation and speciation in Anolis lizards. These subspecies inhabit distinct ecological environments and display marked differences in adult male color and pattern. We sequenced the complete genomes of 20 anoles, 10 from each subspecies, at 1.4× coverage. We used genome-wide scans of population differentiation, allele frequency spectrum, and linkage disequilibrium to characterize the genomic architecture within and between the subspecies. While most of the genome was undifferentiated, we observed five large divergent regions. Within these regions we identified blocks, 5 kb pairs in length, enriched for fixed single nucleotide polymorphisms. These blocks encompass 97 genes, two of which are candidate pigmentation genes. One is melanophilin (mlph), which helps transport melanosomes within melanocytes. The other is a cluster of differentiation 36 (cd36), which regulates carotenoid pigment sequestration. We used high-pressure liquid chromatography to confirm that carotenoid pigments are significantly more abundant in the conspicuous orange-pigmented skin of male A. m. marmoratus suggesting that cd36 may be regulating pigment deposition in this tissue. We identified for the first time a carotenoid gene that is a potential target of divergent sexual selection and may be contributing to the early stages of speciation in Anolis lizards.

Location and Species Matters: Variable Influence of the Environment on the Gene Flow of Imperiled, Native and Invasive Cottontails.

The environment plays an important role in the movement of individuals and their associated genes among populations, which facilitates gene flow. Gene flow can help maintain the genetic diversity both within and between populations and counter the negative impact of genetic drift, which can decrease the fitness of individuals. Sympatric species can have different habitat preferences, and thus can exhibit different patterns of genetic variability and population structure. The specialist-generalist variation hypothesis (SGVH) predicts that specialists will have lower genetic diversity, lower effective population sizes (Ne), and less gene flow among populations. In this study, we used spatially explicit, individual-based comparative approaches to test SGVH predictions in two sympatric cottontail species and identify environmental variables that influence their gene flow. New England cottontail (Sylvilagus transitionalis) is the only native cottontail in the Northeast US, an early successional habitat specialist, and a species of conservation concern. Eastern cottontail (S. floridanus) is an invasive species in the Northeast US and a habitat generalist. We characterized each species' genomic variation by developing double-digest Restriction-site Associated DNA sequence single nucleotide polymorphism markers, quantified their habitat with Geographic Information System environmental variables, and conducted our analyses at multiple scales. Surprisingly, both species had similar levels of genetic diversity and eastern cottontail's Ne was only higher than New England cottontail in one of three subregions. At a regional level, the population clusters of New England cottontail were more distinct than eastern cottontail, but the subregional levels showed more geographic areas of restricted gene flow for eastern cottontail than New England cottontail. In general, the environmental variables had the predicted effect on each species' gene flow. However, the most important environmental variable varied by subregion and species, which shows that location and species matter. Our results provide partial support for the SGVH and the identification of environmental variables that facilitate or impede gene flow can be used to help inform management decisions to conserve New England cottontail.

Dispersal ability predicts spatial genetic structure in native mammals persisting across an urbanization gradient.

As the rate of urbanization continues to increase globally, a growing body of research is emerging that investigates how urbanization shapes the movement-and consequent gene flow-of species in cities. Of particular interest are native species that persist in cities, either as small relict populations or as larger populations of synanthropic species that thrive alongside humans in new urban environments. In this study, we used genomic sequence data (SNPs) and spatially explicit individual-based analyses to directly compare the genetic structure and patterns of gene flow in two small mammals with different dispersal abilities that occupy the same urbanized landscape to evaluate how mobility impacts genetic connectivity. We collected 215 white-footed mice (Peromyscus leucopus) and 380 big brown bats (Eptesicus fuscus) across an urban-to-rural gradient within the Providence, Rhode Island (U.S.A.) metropolitan area (population =1,600,000 people). We found that mice and bats exhibit clear differences in their spatial genetic structure that are consistent with their dispersal abilities, with urbanization having a stronger effect on Peromyscus mice. There were sharp breaks in the genetic structure of mice within the Providence urban core, as well as reduced rates of migration and an increase in inbreeding with more urbanization. In contrast, bats showed very weak genetic structuring across the entire study area, suggesting a near-panmictic gene pool likely due to the ability to disperse by flight. Genetic diversity remained stable for both species across the study region. Mice also exhibited a stronger reduction in gene flow between island and mainland populations than bats. This study represents one of the first to directly compare multiple species within the same urban-to-rural landscape gradient, an important gap to fill for urban ecology and evolution. Moreover, here we document the impacts of dispersal capacity on connectivity for native species that have persisted as the urban landscape matrix expands.

Mitochondrial DNA analysis of critically endangered Chinese Pangolins (Manis pentadactyla) from Nepal.

Chinese Pangolins (Manis pentadactyla) are Critically Endangered and one of the most illegally traded mammals globally. We generated first COI sequences from five individuals of this species from Nepal. BLASTn search of our 600 bp sequences at GenBank showed pair-wise identity between 99.17% and 100% to M. pentadactyla. There were three haplotypes and a total of five variable sites among five M. pentadactyla sequences. Neighbor-joining tree revealed that all M. pentadactyla from Nepal clustered into same group further splitting into two sub-groups albeit with low bootstrap value, suggesting potential multiple geographic origins. The K2P distance was 0.3% within group and 0.7% between four sequences from Bhaktapur and Kavrepalanchok districts (Mape2, Mape3, Mape5 and Mape6) and museum sample (Mape10). This study has generated reference samples for M. pentadactyla from Nepal and will be helpful in understanding dynamics of illegal trade of this species and in successful identification of M. pentadactyla from Nepal even in the absence of intact specimens.

Microsatellite marker development from next-generation sequencing in the New England cottontail (Sylvilagus transitionalis) and cross-amplification in the eastern cottontail (S. floridanus).

OBJECTIVE: The New England cottontail (Sylvilagus transitionalis) is a species of high conservation priority in the Northeastern United States, and was a candidate for federal listing under the Endangered Species Act until a recent decision determined that conservation actions were sufficient to preclude listing. The aim of this study was to develop a suite of microsatellite loci to guide future research efforts such as the analysis of population genetic structure, genetic variation, dispersal, and genetic mark-recapture population estimation. RESULTS: Thirty-five microsatellite markers containing tri- and tetranucleotide sequences were developed from shotgun genomic sequencing of tissue from S. transitionalis, S. obscurus, and S. floridanus. These loci were screened in n = 33 wild S. transitionalis sampled from a population in eastern Massachusetts, USA. Thirty-two of the 35 loci were polymorphic with 2-6 alleles, and observed heterozygosities of 0.06-0.82. All loci conformed to Hardy-Weinberg Equilibrium proportions and there was no evidence of linkage disequilibrium or null alleles. Primers for 33 of the 35 loci amplified DNA extracted from n = 6 eastern cottontail (S. floridanus) samples, of which nine revealed putative species-diagnostic alleles. These loci will provide a useful tool for conservation genetics investigations of S. transitionalis and a potential diagnostic species assay for differentiating sympatric eastern and New England cottontails.

Comparative mtDNA analyses of three sympatric macropodids from a conservation area on the Huon Peninsula, Papua New Guinea.

Matschie's tree kangaroo (Dendrolagus matschiei), New Guinea pademelon (Thylogale browni), and small dorcopsis (Dorcopsulus vanheurni) are sympatric macropodid taxa, of conservation concern, that inhabit the Yopno-Urawa-Som (YUS) Conservation Area on the Huon Peninsula, Papua New Guinea. We sequenced three partial mitochondrial DNA (mtDNA) genes from the three taxa to (i) investigate network structure; and (ii) identify conservation units within the YUS Conservation Area. All three taxa displayed a similar pattern in the spatial distribution of their mtDNA haplotypes and the Urawa and Som rivers on the Huon may have acted as a barrier to maternal gene flow. Matschie's tree kangaroo and New Guinea pademelon within the YUS Conservation Area should be managed as single conservation units because mtDNA nucleotides were not fixed for a given geographic area. However, two distinct conservation units were identified for small dorcopsis from the two different mountain ranges within the YUS Conservation Area.

An Analysis of Overstory Tree Canopy Cover in Sites Occupied by Native and Introduced Cottontails in the Northeastern United States with Recommendations for Habitat Management for New England Cottontail.

The New England cottontail (Sylvilagus transitionalis) is a high conservation priority in the Northeastern United States and has been listed as a candidate species under the Endangered Species Act. Loss of early successional habitat is the most common explanation for the decline of the species, which is considered to require habitat with dense low vegetation and limited overstory tree canopy. Federal and state wildlife agencies actively encourage landowners to create this habitat type by clearcutting blocks of forest. However, there are recent indications that the species also occupies sites with moderate overstory tree canopy cover. This is important because many landowners have negative views about clearcutting and are more willing to adopt silvicultural approaches that retain some overstory trees. Furthermore, it is possible that clearcuts with no overstory canopy cover may attract the eastern cottontail (S. floridanus), an introduced species with an expanding range. The objective of our study was to provide guidance for future efforts to create habitat that would be more favorable for New England cottontail than eastern cottontail in areas where the two species are sympatric. We analyzed canopy cover at 336 cottontail locations in five states using maximum entropy modelling and other statistical methods. We found that New England cottontail occupied sites with a mean overstory tree canopy cover of 58% (SE±1.36), and was less likely than eastern cottontail to occupy sites with lower overstory canopy cover and more likely to occupy sites with higher overstory canopy cover. Our findings suggest that silvicultural approaches that retain some overstory canopy cover may be appropriate for creating habitat for New England cottontail. We believe that our results will help inform critical management decisions for the conservation of New England cottontail, and that our methodology can be applied to analyses of habitat use of other critical wildlife species.

Divergence in coloration and ecological speciation in the Anolis marmoratus species complex.

Adaptive divergence in coloration is expected to produce reproductive isolation in species that use colourful signals in mate choice and species recognition. Indeed, many adaptive radiations are characterized by differentiation in colourful signals, suggesting that divergent selection acting on coloration may be an important component of speciation. Populations in the Anolis marmoratus species complex from the Caribbean island of Guadeloupe display striking divergence in the colour and pattern of adult males that occurs over small geographic distances, suggesting strong divergent selection. Here we test the hypothesis that divergence in coloration results in reduced gene flow among populations. We quantify variation in adult male coloration across a habitat gradient between mesic and xeric habitats, use a multilocus coalescent approach to infer historical demographic parameters of divergence, and examine gene flow and population structure using microsatellite variation. We find that colour variation evolved without geographic isolation and in the face of gene flow, consistent with strong divergent selection and that both ecological and sexual selection are implicated. However, we find no significant differentiation at microsatellite loci across populations, suggesting little reproductive isolation and high levels of contemporary gene exchange. Strong divergent selection on loci affecting coloration probably maintains clinal phenotypic variation despite high gene flow at neutral loci, supporting the notion of a porous genome in which adaptive portions of the genome remain fixed whereas neutral portions are homogenized by gene flow and recombination. We discuss the impact of these findings for studies of colour evolution and ecological speciation.

Genetic evaluation of the Association of Zoos and Aquariums Matschie's tree kangaroo (Dendrolagus matschiei) captive breeding program.

Matschie's tree kangaroo (Dendrolagus matschiei) is an endangered species that has been bred in captivity since the 1970s. In 1992, the Tree Kangaroo Species Survival Plan(®) (TKSSP) was established to coordinate the captive management of Association of Zoos and Aquariums (AZA) D. matschiei. The TKSSP makes annual breeding recommendations primarily based on the mean kinship (MK) strategy. Captive breeding programs often use the MK strategy to preserve genetic diversity in small populations-to avoid the negative consequences of inbreeding and retain their adaptive potential. The ability of a captive breeding program to retain the population's genetic diversity over time can be evaluated by comparing the genetic diversity of the captive population to wild populations. We analyzed DNA extracted from blood and fecal samples from AZA (n = 71), captive (n = 28), and wild (n = 22) D. matschiei using eight microsatellite markers and sequenced the partial mitochondrial DNA control region gene. AZA D. matschiei had a similar expected heterozygosity (H(e) = 0.595 ± 0.184) compared with wild D. matschiei (H(e) = 0.628 ± 0.143), but they had different allelic frequencies (F(ST) = 0.126; P < 0.001). AZA D. matschiei haplotype diversity was almost two times lower than wild D. matschiei H = 0.740 ± 0.063. These data will assist management of AZA D. matschiei and serve as a baseline for AZA and wild D. matschiei genetic diversity values that could be used to monitor future changes in their genetic diversity.

Population genetic diversity and fitness in multiple environments.

BACKGROUND: When a large number of alleles are lost from a population, increases in individual homozygosity may reduce individual fitness through inbreeding depression. Modest losses of allelic diversity may also negatively impact long-term population viability by reducing the capacity of populations to adapt to altered environments. However, it is not clear how much genetic diversity within populations may be lost before populations are put at significant risk. Development of tools to evaluate this relationship would be a valuable contribution to conservation biology. To address these issues, we have created an experimental system that uses laboratory populations of an estuarine crustacean, Americamysis bahia with experimentally manipulated levels of genetic diversity. We created replicate cultures with five distinct levels of genetic diversity and monitored them for 16 weeks in both permissive (ambient seawater) and stressful conditions (diluted seawater). The relationship between molecular genetic diversity at presumptive neutral loci and population vulnerability was assessed by AFLP analysis. RESULTS: Populations with very low genetic diversity demonstrated reduced fitness relative to high diversity populations even under permissive conditions. Population performance decreased in the stressful environment for all levels of genetic diversity relative to performance in the permissive environment. Twenty percent of the lowest diversity populations went extinct before the end of the study in permissive conditions, whereas 73% of the low diversity lines went extinct in the stressful environment. All high genetic diversity populations persisted for the duration of the study, although population sizes and reproduction were reduced under stressful environmental conditions. Levels of fitness varied more among replicate low diversity populations than among replicate populations with high genetic diversity. There was a significant correlation between AFLP diversity and population fitness overall; however, AFLP markers performed poorly at detecting modest but consequential losses of genetic diversity. High diversity lines in the stressful environment showed some evidence of relative improvement as the experiment progressed while the low diversity lines did not. CONCLUSIONS: The combined effects of reduced average fitness and increased variability contributed to increased extinction rates for very low diversity populations. More modest losses of genetic diversity resulted in measurable decreases in population fitness; AFLP markers did not always detect these losses. However when AFLP markers indicated lost genetic diversity, these losses were associated with reduced population fitness.

Microsatellite marker development and Mendelian analysis in the Matschie's tree kangaroo (Dendrolagus matschiei).

Matschie's tree kangaroo (Dendrolagus matschiei) is an endangered arboreal macropodid endemic to the Huon Peninsula, Papua New Guinea (PNG). We developed 5 microsatellite markers for D. matschiei, which are the first markers developed for Dendrolagus. We screened 17 additional markers that were developed for other marsupial taxa and identified 3 that were polymorphic in D. matschiei. We estimated allelic and genetic diversity with the set of 8 markers by analyzing 22 D. matschiei from Wasaunon on the Huon Peninsula, PNG. The number of alleles ranged from 2 to 9 and expected heterozygosity ranged from 0.440 to 0.794. We tested for null alleles and Mendelian inheritance by analyzing 19 pairs of D. matschiei parents and offspring from Association of Zoos and Aquariums institutions. Null alleles were not detected and Mendelian inheritance was followed for all 8 markers. We also evaluated the reliability of using the markers to amplify DNA extracted from D. matschiei fecal samples and the ability of the markers to amplify DNA samples from Goodfellow's tree kangaroo (Dendrolagus goodfellowi ssp.), Doria's tree kangaroo (Dendrolagus dorianus ssp.), and Grizzled tree kangaroo (Dendrolagus inustus ssp.). Microsatellite markers can be used to inform management decisions to conserve D. matschiei in captivity and the wild.

Genetic diversity in captive and wild Matschie's tree kangaroo (Dendrolagus matschiei) from Huon Peninsula, Papua New Guinea, based on mtDNA control region sequences.

The Association of Zoos and Aquariums (AZA) Matschie's tree kangaroo (Dendrolagus matschiei) population is at a critical point for assessing long-term viability. This population, established from 19 genetically uncharacterized D. matschiei, has endured a founder effect because only four individuals contributed the majority of offspring. The highly variable mitochondrial DNA (mtDNA) control region was sequenced for five of the female-founders by examining extant representatives of their maternal lineage and compared with wild (n = 13) and captive (n = 18) D. matschiei from Papua New Guinea (PNG). AZA female-founder D. matschiei control region haplotype diversity was low, compared with captive D. matschiei held in PNG. AZA D. matschiei have only two control region haplotypes because four out of five AZA female-founder D. matschiei had an identical sequence. Both AZA haplotypes were identified among the 17 wild and captive D. matschiei haplotypes from PNG. Genomic DNA extracted from wild D. matschiei fecal samples was a reliable source of mtDNA that could be used for a larger scale study. We recommend a nuclear DNA genetic analysis to more fully characterize AZA D. matschiei genetic diversity and to assist their Species Survival Plan((R)). An improved understanding of D. matschiei genetics will contribute substantially to the conservation of these unique animals both in captivity and the wild.