Pesticide exposure affects DNA methylation patterns in natural populations of a mayfly.

Chemical pollutants derived from agricultural activities represent a major threat to freshwater biota. Despite growing evidence involving epigenetic processes, such as DNA methylation, in response to pesticide contamination in agroecosystems, research on wild populations of non-model species remains scarce, particularly for endemic freshwater arthropods. Using the MethylRAD method, this study investigates whether exposure to pesticide contamination in natural populations of the endemic mayfly A. torrens produces genome wide changes in levels of DNA methylation. From a total of 1,377,147 MethylRAD markers produced from 285 specimens collected at 30 different study sites along the Limarí watershed of north-central Chile, six showed significant differential methylation between populations exposed and unexposed to pesticides. In all cases the effect of pesticides was positive, independent and stronger than the effects detected for other spatial and environmental factors. Only one candidate marker appeared correlated significantly with additional variables, nitrate and calcium levels, which also reflects the impact of agrichemicals and could additionally suggest, to a lower extent, antagonistic effects of mineral salts concentration for this specific marker. These results suggest that the effect of pesticide exposure on methylation levels is apparent at these six MethylRAD markers in A. torrens populations. Such data is challenging to obtain in natural populations and is, for the most part, lacking in ecotoxicological studies. Our study shows that DNA methylation processes are involved in the response to pesticide contamination in populations of the mayfly A. torrens in their natural habitat, and provides new evidence regarding the impact of pesticide contamination and agricultural activities on the endemic fauna of lotic ecosystems.

Genome-wide genetic diversity yields insights into genomic responses of candidate climate-selected loci in an Andean wetland plant.

Assessing population evolutionary potential has become a central tenet of conservation biology. Since adaptive responses require allelic variation at functional genes, consensus has grown that genetic variation at genes under selection is a better surrogate for adaptive evolutionary potential than neutral genetic diversity. Although consistent with prevailing theory, this argument lacks empirical support and ignores recent theoretical advances questioning the very concept of neutral genetic diversity. In this study, we quantified genome-wide responses of single nucleotide polymorphism loci linked to climatic factors over a strong latitudinal gradient in natural populations of the high Andean wetland plant, Carex gayana, and then assessed whether genetic variation of candidate climate-selected loci better predicted their genome-wide responses than genetic variation of non-candidate loci. Contrary to this expectation, genomic responses of climate-linked loci only related significantly to environmental variables and genetic diversity of non-candidate loci. The effects of genome-wide genetic diversity detected in this study may be a result of either the combined influence of small effect variants or neutral and demographic factors altering the adaptive evolutionary potential of C. gayana populations. Regardless of the processes involved, our results redeem genome-wide genetic diversity as a potentially useful indicator of population adaptive evolutionary potential.

Developmental instability is associated with estrogenic endocrine disruption in the Chilean native fish species, Trichomycterus areolatus.

Endocrine disrupting compounds (EDCs) are widespread contaminants that alter the normal functioning of the endocrine system. While they cause dysfunctions in essential biological processes, it is unclear whether EDCs also impact developmental stability. In the present study, we investigated the occurrence of estrogenic endocrine disrupting compounds in a small watershed of south-central Chile impacted by anthropogenic activities. Then, we assessed their relationship with internal levels of estrogenic active compounds and fluctuating asymmetry (FA), a proxy of developmental stability in organisms with bilateral symmetry, in a native fish species (Trichomycterus areolatus). Yeast estrogenic screen assays were performed to measure estrogenic activity in river sediments and in male fish tissues collected from 17 sites along the Chillán watershed, and geometric morphometrics used to estimate fluctuating asymmetry based on the shapes of 248 fish skulls. Estrogenic activity was detected both in sediments and male fish tissues at concentrations of up to 1005 ng and 83 ng 17ß-estradiol equivalent/kg dw, respectively. No significant correlation was found between the two. However, fish tissue estrogenicity, water temperature and dissolved oxygen explained >80% of the FA population variation. By showing a significant relationship between estrogenic activity and FA of T. areolatus, our results indicate that developmental stability can be altered by estrogenic endocrine disruption, and that FA can be a useful indicator of sub-lethal stress in T. areolatus populations.

Assessment of Gene Expression Biomarkers in the Chilean Pencil Catfish, Trichomycterus areolatus, from the Choapa River Basin, Coquimbo Chile.

The objective of this study was to describe changes in the gene expression in the Chilean catfish, Trichomycterus areolatus, based on their geographic location within the Choapa River. Genes of choice included those that are biomarkers of exposure to metals, oxidative stress, and endocrine disruption. Male and female T. areolatus were sampled from four sites in January 2015 differently impacted by human activities. In males, but not females, hepatic gene expression of heat shock protein (HSP70) and cytochrome P450 1A (CYP1A) were significantly elevated at the site adjacent to the small city of Salamanca, relative to the other sites. In females, hepatic HSP70, the aryl hydrocarbon receptor (AHR), and the estrogen responsive genes, vitellogenin (VTG) and estrogen receptor alpha (ERα), were significantly lower at the site located furthest downstream. A similar downstream pattern of lower expression levels also was found in ovarian tissue for the genes, HSP70 and ERα. Gill gene expression showed a unique pattern in females as levels of metallothionein were elevated at the site furthest downstream. While analytical chemistry of water samples provided limited evidence of agrichemical contamination, the gene expression data are consistent with an exposure to agrichemicals and metals. T. areolatus may be a valuable sentinel organism and its use as a bioindicator species in some rivers within Chile can provide considerable insight, particularly in situations analytical chemistry is limited by environmental constraints.

Pesticide contamination drives adaptive genetic variation in the endemic mayfly Andesiops torrens within a semi-arid agricultural watershed of Chile.

Agrichemical contamination can provoke evolutionary responses in freshwater populations. It is a particularly relevant issue in semi-arid regions due to the sensitivity of endemic species to pollutants and to interactions with temperature stress. This paper investigates the presence of pesticides in rivers within a semi-arid agricultural watershed of Chile, testing for their effects on population genetic characteristics of the endemic mayfly Andesiops torrens (Insecta, Ephemeroptera). Pesticides were detected in sediment samples in ten out of the 30 sites analyzed throughout the upper part of the Limarí watershed. To study the evolutionary impact of such contamination on A. torrens, we used a genome-wide approach and analyzed 2056 single nucleotide polymorphisms (SNPs) loci in 551 individuals from all sites. Genetic differentiation was weak between populations, suggesting high gene flow across the study area. While we did not find evidence of pesticide effects on genetic diversity nor on population differentiation, the allele frequency of three outlier SNP loci correlated significantly with pesticide occurrence. Interrogation of genomic resources indicates that two of these SNPs are located within functional genes that encode for the low-density lipoprotein receptor-related protein 2 and Dumpy, both potentially involved in insect cuticle resistance processes. Such genomic signatures of local adaptation are indicative of past adverse effects of pesticide exposure on the locally adapted populations. Our results reveal that A. torrens is sensitive to pesticide exposure, but that a high gene flow may confer resilience to contamination. This research supports the contention that A. torrens is an ideal model organism to study evolutionary responses induced by pesticides on non-target, endemic species.

Unexpected population fragmentation in an endangered seabird: the case of the Peruvian diving-petrel.

In less than one century, the once-abundant Peruvian diving petrel has become the first endangered seabird of the Humboldt Current System (HCS). This small endemic petrel of the South American Pacific coast is now an important indicator of ongoing habitat loss and of the success of local conservation policies in the HCS - an ecoregion designated as a priority for the conservation of global biodiversity. Yet so far, poorly understood life history traits such as philopatry or dispersal ability may strongly influence the species' response to ecosystem changes, but also our capacity to assess and interpret this response. To address this question, we explore the range-wide population structure of the Peruvian diving petrel, and show that this small seabird exhibits extreme philopatric behavior at the island level. Mitochondrial DNA sequences and genome-wide SNP data reveal significant isolation and low migration at very short distances, and provide strong evidence for questioning the alleged recovery in the Peruvian and Chilean populations of this species. Importantly, the full demographic independence between colonies makes local population rescue through migration unlikely. As a consequence, the Peruvian diving petrel appears to be particularly vulnerable to ongoing anthropogenic pressure. By excluding immigration as a major factor of demographic recovery, our results highlight the unambiguously positive impact of local conservation measures on breeding populations; yet at the same time they also cast doubt on alleged range-wide positive population trends. Overall, the protection of independent breeding colonies, and not only of the species as a whole, remains a major element in the conservation strategy for endemic seabirds. Finally, we underline the importance of considering the philopatric behavior and demographic independence of breeding populations, even at very fine spatial scales, in spatial planning for marine coastal areas.

Isolation and cross-amplification of the first set of polymorphic microsatellite markers of two high-Andean cushion plants.

In the southern Andes mountains (27-39◦S) Azorella madreporica and Laretia acaulis, two Apiaceae cushion plant species commonly known as yaretas, conform a well-established altitudinal vegetation belt along the lower Andean zone. These species have been considered as fundamental components of several ecological dynamics within their communities; however, high mountain ecosystems are increasingly threatened worldwide by natural and anthropogenic pressures and the southern Andes are not the exception. Recognizing that genetic information is crucial for the success of any conservation or restoration initiative inwild populations, we developed and cross-amplified 28 specifically designed microsatellite markers (14 in A. madreporica and 14 in L. acaulis), and also tested the cross amplification of 25 markers from the related species Azorella selago. In a region which is particularly vulnerable to global change trends, this new polymorphic microsatellite loci will be useful in the study of the genetic diversity of these high-mountain cushion plants, which are pivotal in the structuring of their native ecosystems.

Landscape connectivity among remnant populations of guanaco (Lama guanicoe Müller, 1776) in an arid region of Chile impacted by global change.

Connectivity between populations plays a key role in the long-term persistence of species in fragmented habitats. This is of particular concern for biodiversity preservation in drylands, since water limited landscapes are typically characterized by little suitable habitat cover, high habitat fragmentation, harsh matrices, and are being rapidly degraded at a global scale. In this study, we modelled landscape connectivity between 11 guanaco Lama guanicoe populations in Chile's arid Norte Chico, a region that supports the last remnant coastal populations of this emblematic herbivore indigenous to South America. We produced a habitat suitability model to derive a regional surface resistance map, and used circuit theory to map functional connectivity, investigate the relative isolation between populations, and identify those that contribute most to the patch connectivity network. Predicted suitable habitat for L. guanicoe represented about 25% of the study region (i.e., 29,173 km2) and was heterogeneously distributed along a continuous stretch along the Andes, and discontinuous patches along the coast. As a result, we found that high connectivity current flows in the mid and high Andes formed a wide, continuous connectivity corridor, enabling connectivity between all high Andean populations. Coastal populations, in contrast, were more isolated. These groups demonstrate no inter-population connectivity between themselves, only with higher altitude populations, and for two of them, animal movement was linked to the effectiveness of wildlife crossings along the Pan-American highway. Our results indicate that functional connectivity is an issue of concern for L. guanicoe in Chile's Norte Chico, implying that future conservation and management plans should emphasize strategies aimed at conserving functional connectivity between coastal and Andean populations, as well as the protection of habitat patches likely to act as stepping stones within the connectivity network.

Partitioning genetic and species diversity refines our understanding of species-genetic diversity relationships.

Disentangling the origin of species-genetic diversity correlations (SGDCs) is a challenging task that provides insight into the way that neutral and adaptive processes influence diversity at multiple levels. Genetic and species diversity are comprised by components that respond differently to the same ecological processes. Thus, it can be useful to partition species and genetic diversity into their different components to infer the mechanisms behind SGDCs. In this study, we applied such an approach using a high-elevation Andean wetland system, where previous evidence identified neutral processes as major determinants of the strong and positive covariation between plant species richness and AFLP genetic diversity of the common sedge Carex gayana. To tease apart putative neutral and non-neutral genetic variation of C. gayana, we identified loci putatively under selection from a dataset of 1,709 SNPs produced using restriction site-associated DNA sequencing (RAD-seq). Significant and positive relationships between local estimates of genetic and species diversities (α-SGDCs) were only found with the putatively neutral loci datasets and with species richness, confirming that neutral processes were primarily driving the correlations and that the involved processes differentially influenced local species diversity components (i.e., richness and evenness). In contrast, SGDCs based on genetic and community dissimilarities (ß-SGDCs) were only significant with the putative non-neutral datasets. This suggests that selective processes influencing C. gayana genetic diversity were involved in the detected correlations. Together, our results demonstrate that analyzing distinct components of genetic and species diversity simultaneously is useful to determine the mechanisms behind species-genetic diversity relationships.

Cynomolgus macaques naturally infected with Trypanosoma cruzi-I exhibit an overall mixed pro-inflammatory/modulated cytokine signature characteristic of human Chagas disease.

BACKGROUND: Non-human primates have been shown to be useful models for Chagas disease. We previously reported that natural T. cruzi infection of cynomolgus macaques triggers clinical features and immunophenotypic changes of peripheral blood leukocytes resembling those observed in human Chagas disease. In the present study, we further characterize the cytokine-mediated microenvironment to provide supportive evidence of the utility of cynomolgus macaques as a model for drug development for human Chagas disease. METHODS AND FINDINGS: In this cross-sectional study design, flow cytometry and systems biology approaches were used to characterize the ex vivo and in vitro T. cruzi-specific functional cytokine signature of circulating leukocytes from TcI-T. cruzi naturally infected cynomolgus macaques (CH). Results showed that CH presented an overall CD4+-derived IFN-γ pattern regulated by IL-10-derived from CD4+ T-cells and B-cells, contrasting with the baseline profile observed in non-infected hosts (NI). Homologous TcI-T. cruzi-antigen recall in vitro induced a broad pro-inflammatory cytokine response in CH, mediated by TNF from innate/adaptive cells, counterbalanced by monocyte/B-cell-derived IL-10. TcIV-antigen triggered a more selective cytokine signature mediated by NK and T-cell-derived IFN-γ with modest regulation by IL-10 from T-cells. While NI presented a cytokine network comprised of small number of neighborhood connections, CH displayed a complex cross-talk amongst network elements. Noteworthy, was the ability of TcI-antigen to drive a complex global pro-inflammatory network mediated by TNF and IFN-γ from NK-cells, CD4+ and CD8+ T-cells, regulated by IL-10+CD8+ T-cells, in contrast to the TcIV-antigens that trigger a modest network, with moderate connecting edges. CONCLUSIONS: Altogether, our findings demonstrated that CH present a pro-inflammatory/regulatory cytokine signature similar to that observed in human Chagas disease. These data bring additional insights that further validate these non-human primates as experimental models for Chagas disease.

Genetic variation of loci potentially under selection confounds species-genetic diversity correlations in a fragmented habitat.

Positive species-genetic diversity correlations (SGDCs) are often thought to result from the parallel influence of neutral processes on genetic and species diversity. Yet, confounding effects of non-neutral mechanisms have not been explored. Here, we investigate the impact of non-neutral genetic diversity on SGDCs in high Andean wetlands. We compare correlations between plant species diversity and genetic diversity (GD) calculated with and without loci potentially under selection (outlier loci). The study system includes 2188 specimens from five species (three common aquatic macroinvertebrate and two dominant plant species) that were genotyped for 396 amplified fragment length polymorphism loci. We also appraise the importance of neutral processes on SGDCs by investigating the influence of habitat fragmentation features. Significant positive SGDCs were detected for all five species (mean SGDC = 0.52 ± 0.05). While only a few outlier loci were detected in each species, they resulted in significant decreases in GD and in SGDCs. This supports the hypothesis that neutral processes drive species-genetic diversity relationships in high Andean wetlands. Unexpectedly, the effects on genetic diversity GD of the habitat fragmentation characteristics in this study increased with the presence of outlier loci in two species. Overall, our results reveal pitfalls in using habitat features to infer processes driving SGDCs and show that a few loci potentially under selection are enough to cause a significant downward bias in SGDC. Investigating confounding effects of outlier loci thus represents a useful approach to evidence the contribution of neutral processes on species-genetic diversity relationships.

De novo Assembly and Analysis of the Chilean Pencil Catfish Trichomycterus areolatus Transcriptome.

Trichomycterus areolatus is an endemic species of pencil catfish that inhabits the riffles and rapids of many freshwater ecosystems of Chile. Despite its unique adaptation to Chile's high gradient watersheds and therefore potential application in the investigation of ecosystem integrity and environmental contamination, relatively little is known regarding the molecular biology of this environmental sentinel. Here, we detail the assembly of the Trichomycterus areolatus transcriptome, a molecular resource for the study of this organism and its molecular response to the environment. RNA-Seq reads were obtained by next-generation sequencing with an Illumina® platform and processed using PRINSEQ. The transcriptome assembly was performed using TRINITY assembler. Transcriptome validation was performed by functional characterization with KOG, KEGG, and GO analyses. Additionally, differential expression analysis highlights sex-specific expression patterns, and a list of endocrine and oxidative stress related transcripts are included.

Lizards on ice: evidence for multiple refugia in Liolaemus pictus (Liolaemidae) during the last glacial maximum in the Southern Andean beech forests.

Historical climate changes and orogenesis are two important factors that have shaped intraspecific biodiversity patterns worldwide. Although southern South America has experienced such complex events, there is a paucity of studies examining the effects on intraspecific diversification in this part of the world. Liolaemus pictus is the southernmost distributed lizard in the Chilean temperate forest, whose genetic structure has likely been influenced by Pleistocene glaciations. We conducted a phylogeographic study of L. pictus in Chile and Argentina based on one mitochondrial and two nuclear genes recovering two strongly divergent groups, Northern and Southern clades. The first group is distributed from the northernmost limit of the species to the Araucanía region while the second group is distributed throughout the Andes and the Chiloé archipelago in Southern Chile. Our results suggest that L. pictus originated 751 Kya, with divergence between the two clades occurring in the late Pleistocene. Demographic reconstructions for the Northern and Southern clades indicate a decrease in effective population sizes likely associated with Pleistocene glaciations. Surprisingly, patterns of genetic variation, clades age and historical gene flow in populations distributed within the limits of the Last Glacial Maximum (LGM) are not explained by recent colonization. We propose an "intra-Andean multiple refuge" hypothesis, along with the classical refuge hypothesis previously proposed for the biota of the Chilean Coastal range and Eastern Andean Cordillera. Our hypothesis is supported by niche modelling analysis suggesting the persistence of fragments of suitable habitat for the species within the limits of the LGM ice shield. This type of refuge hypothesis is proposed for the first time for an ectothermic species.

The role of spatial processes and environmental determinants in microgeographic shell variation of the freshwater snail Chilina dombeyana (Bruguière, 1789).

Wildlife data often show spatial organization, demonstrating positive correlations either as a result of processes occurring over the landscape or due to the influence of spatially structured environmental variables. It is, thus, essential to consider non-random spatial structure when evaluating the underlying causes of biological variation. In this study, we analyzed the population structure of Chilina dombeyana shell morphology of 14 populations that are close geographically and belong to the same hydrographic basin. We utilized a variation partitioning approach to evaluate the importance of spatial processes, such as migration, acting over the landscape, and environmental characteristics, including habitat and hydrologic characteristics, and the occurrence of aquatic predators in promoting between population variation. Our results demonstrate spatially structured variation in C. dombeyana shell morphology, with populations living near each other having more similar shell sizes than populations living farther apart. The shell size variation partition indicated that both spatially structured environmental factors and genetic relationships resulting from migration or shared common ancestry may explain this pattern. Shell shape variation, in contrast, was found to be essentially under the influence of non-spatially structured environmental factors, with habitat and water characteristics accounting for about half of the total variation among populations. The large proportion of the variation in shell size that is spatially structured demonstrates that spatial structure on morphological traits might be strong and highlights the need to consider such phenomenon in intraspecific studies of phenotypic evolution.

Permanent genetic resources added to Molecular Ecology Resources Database 1 February 2011-31 March 2011.

This article documents the addition of 111 microsatellite marker loci to the Molecular Ecology Resources Database. Loci were developed for the following species: Acipenser oxyrinchus desotoi, Anopheles nuneztovari sensu lato, Asellus aquaticus, Calopteryx splendens, Calopteryx virgo, Centaurea aspera, Centaurea seridis, Chilina dombeyana, Proctoeces cf. lintoni and Pyrenophora teres f. teres.

Localization of genes for V+LDL plasma cholesterol levels on two diets in the opossum Monodelphis domestica.

Plasma cholesterol levels among individuals vary considerably in response to diet. However, the genes that influence this response are largely unknown. Non-HDL (V+LDL) cholesterol levels vary dramatically among gray, short-tailed opossums fed an atherogenic diet, and we previously reported that two quantitative trait loci (QTLs) influenced V+LDL cholesterol on two diets. We used hypothesis-free, genome-wide linkage analyses on data from 325 pedigreed opossums and located one QTL for V+LDL cholesterol on the basal diet on opossum chromosome 1q [logarithm of the odds (LOD) = 3.11, genomic P = 0.019] and another QTL for V+LDL on the atherogenic diet (i.e., high levels of cholesterol and fat) on chromosome 8 (LOD = 9.88, genomic P = 5 x 10(-9)). We then employed a novel strategy involving combined analyses of genomic resources, expression analysis, sequencing, and genotyping to identify candidate genes for the chromosome 8 QTL. A polymorphism in ABCB4 was strongly associated (P = 9 x 10(-14)) with the plasma V+LDL cholesterol concentrations on the high-cholesterol, high-fat diet. The results of this study indicate that genetic variation in ABCB4, or closely linked genes, is responsible for the dramatic differences among opossums in their V+LDL cholesterol response to an atherogenic diet.

A microsatellite-based, physically anchored linkage map for the gray, short-tailed opossum (Monodelphis domestica).

The genome of the gray, short-tailed opossum, Monodelphis domestica, will be the first of any marsupial to be fully sequenced. The utility of this sequence will be greatly enhanced by construction and integration of detailed genetic and physical maps. Therefore, it is important to verify the unusual recombinational characteristics that were suggested by the 'first-generation' M. domestica linkage map; specifically, very low levels of recombination and severely reduced female recombination, both of which are contrary to patterns in other vertebrates. We constructed a new linkage map based on a different genetic cross, using a new and much larger set of map markers, and physically anchored and oriented the linkage groups onto chromosomes via fluorescence in-situ hybridization mapping. This map includes 150 loci in eight autosomal linkage groups corresponding to the eight autosome pairs, and spans 86-89% of the autosomal genome. The sex-averaged autosomal map covers 715 cM, with a full-length estimate of 866 cM; the shortest full-length linkage map reported for any vertebrate. The sex-specific maps confirmed severely reduced female recombination in all linkage groups, and an overall F/M map ratio = 0.54. These results greatly extend earlier findings, and provide an improved microsatellite-based linkage map for this species.

Modo-UG, a marsupial nonclassical MHC class I locus.

Modo-UG is a class I gene located in the MHC of the marsupial Monodelphis domestica, the gray, short-tailed opossum. Modo-UG is expressed as three alternatively spliced mRNA forms, all of which encode a transmembrane form with a short cytoplasmic tail that lacks phosphorylation sites typically found in classical class I molecules. The three alternative mRNAs would encode a full-length form, an isoform lacking the alpha2 domain, and one lacking both alpha2 and alpha3 domains. Genotyping both captive-bred and wild M. domestica from different geographic regions revealed no variation in the residues that make up Modo-UG's peptide-binding groove. Modo-UG's low polymorphism is contrasting to that of a nearby class I locus, Modo-UA1, which has a highly polymorphic peptide-binding region. Absence of functional polymorphism in Modo-UG is therefore not a general feature of opossum class I genes but the result of negative selection. Modo-UG is the first MHC linked marsupial class I to be described that appears to clearly have nonclassical features.

Reconstructing an ancestral mammalian immune supercomplex from a marsupial major histocompatibility complex.

The first sequenced marsupial genome promises to reveal unparalleled insights into mammalian evolution. We have used the Monodelphis domestica (gray short-tailed opossum) sequence to construct the first map of a marsupial major histocompatibility complex (MHC). The MHC is the most gene-dense region of the mammalian genome and is critical to immunity and reproductive success. The marsupial MHC bridges the phylogenetic gap between the complex MHC of eutherian mammals and the minimal essential MHC of birds. Here we show that the opossum MHC is gene dense and complex, as in humans, but shares more organizational features with non-mammals. The Class I genes have amplified within the Class II region, resulting in a unique Class I/II region. We present a model of the organization of the MHC in ancestral mammals and its elaboration during mammalian evolution. The opossum genome, together with other extant genomes, reveals the existence of an ancestral "immune supercomplex" that contained genes of both types of natural killer receptors together with antigen processing genes and MHC genes.

Development, inheritance, and linkage-group assignment of 60 novel microsatellite markers for the gray, short-tailed opossum Monodelphis domestica.

Short-tandem-repeat (SSR) or microsatellite polymorphisms are some of the most extensively employed genetic markers in contemporary linkage mapping studies. To date, only a limited number of microsatellites have been isolated in the gray, short-tailed opossum Monodelphis domestica, a South American marsupial widely used for comparative biological and biomedical research. To increase the number of potentially useful mapping markers, we screened 2 microsatellite-enriched genomic libraries containing alternatively (CA)n or (GA)n repeats. A total of 184 clones were sequenced, from which 60 polymorphic microsatellite markers were successfully optimized. The efficiency of this enrichment protocol for M. domestica microsatellite isolation is discussed, and suggestions to improve the outcome are made. All 60 loci showed high allelic diversity, with allele numbers ranging from 2 to 10 in a subset of 33 unrelated animals. Normal Mendelian inheritance was confirmed for all loci by analyzing allelic segregation in 5 two-generation families. One microsatellite appeared to be X linked, and null alleles were found in 5 others. Two-point linkage analyses were implemented using the data on the 5 families, leading to the assignment of 59 of these loci to the existing linkage groups. The 60 novel microsatellites developed in this study will contribute significantly to the M. domestica linkage map, and further QTL mapping studies.