Genetic variability of populations of the white-eared opossum, Didelphis albiventris Lund 1840 (Didelphimorphia; Didelphidae) in Brazil / Variabilidade genética de populações de gambá de orelha branca, Didelphis albiventris Lund 1840 (Didelphimorphia; Didelphidae) no Brasil

Abstract Didelphis albiventris are found throughout Northeast and Central Brazil to central-southern Uruguay and it was subject of few studies in a population level. Given this, the present study investigated the genetic variability of the species using the mitochondrial molecular marker cytochrome oxidase c subunit I. We analyzed samples from the different biomes within three Brazilian regions: Northeast (Caatinga , Cerrado, and Atlantic Forest), Southeast (Cerrado , Atlantic Forest, Cerrado/Atlantic Forest, and Cerrado/Caatinga ecotones) and South (Pampa and Atlantic Forest). Software BAPs retrieved five distinct demes: dm 1, dm 2, and dm 5 that occurs in South, Northeast and Southeast regions respectively and the dm 3 and dm 4 are wide distributed in Northeast and Southeast. Population analysis performed with AMOVA, haplotype network and Mantel test estimated the veracity of the demes. The FST shows structuring for the five demes, with dm 1 (South region) isolated from the others, however the other analysis showed the Northeast/Southeast demes (dm 2-5) united, diagnosing gene flow between them, mainly at the transitional zones, in areas as far away as areas with similar latitude interval (Southeast vs South) that was not detected gene flow. In the haplotype network, the mutational steps was conclusive in split dm1 from dm 2-5 with 15 mutational steps and the Mantel test was moderated, which is explained by genetic similarity despite the great geographic distances (Northeast/Southeast). Thus, our analysis recognized two different lineages (South and Northeast/Southeast) and indicate that the biomes were not decisive in their isolation. The sharing of demes at the transitional zones and in areas with high latitudinal intervals reflects a recent ancestral polymorphism for D. albiventris. The plasticity in the occupation of the space by this species contributes in its wide dispersion capability, that is, geographical distribution. Our results revealed important implications for the management of D. albiventris in these transitional zones areas where demes were shared.

Resumo Didelphis albiventris é encontrada em todo o Nordeste e região central do Brasil até o centro-sul do Uruguai e foi alvo de poucos estudos em nível populacional. Dessa forma, o presente estudo, investiga a variabilidade genética da espécie usando o marcador molecular citocromo c oxidase subunidade I. Analisou-se amostras de diferentes biomas de três regiões brasileiras: Nordeste (Caatinga, Cerrado e Floresta Atlântica), Sudeste (Cerrado, Floresta Atlântica, ecótonos Cerrado/Floresta Atlântica e Cerrado/Caatinga) e Sul (Pampa e Floresta Atlântica). O software BAPs recuperou cinco demes distintos: dm 1, dm 2 e dm 5, que ocorrem nas regiões Sul, Nordeste e Sudeste, respectivamente, e os dm 3 e dm 4, que são amplamente distribuído no Nordeste e Sudeste. Análises populacionais realizadas com AMOVA, rede de haplótipo e teste de Mantel estimaram a veracidade das demes. O FST mostrou estruturação para as cinco demes, com dm 1 (região Sul) isolada das demais, entretanto as outras análises mostraram as demes Nordeste/Sudeste (dm 2-5) unidos, diagnosticando fluxo gênico entre elas, principalmente em zonas de transição, em áreas tão distante quanto áreas com similar intervalo de latitude (Sudeste e Sul), onde não foram detectado fluxo gênico. Na rede de haplótipo, os passos mutacionais foram conclusivos em separar dm 1 do dm 2-5 com 15 passos mutacionais, e o teste de Mantel foi moderado, o que é explicado pela similaridade genética apesar da grande distância geográfica (Nordeste/Sudeste). Assim, duas linhagens diferentes (Sul e Sudeste/Nordeste) foram encontradas, indicando que os biomas não foram decisivos em seus isolamentos. Os compartilhamentos das demes, em zonas de transição e em áreas com elevados intervalos de latitude, refletem um polimorfismo ancestral recente para D. albiventris. A plasticidade na ocupação do espaço por esta espécie contribui em sua ampla capacidade de dispersão, ou seja, distribuição geográfica. Nossos resultados revelam importantes implicações para o manejo de D. albiventris nessas áreas de zonas de transição, onde as demes são compartilhadas.

Genetic variability of populations of the white-eared opossum, Didelphis albiventris Lund 1840 (Didelphimorphia; Didelphidae) in Brazil.

Didelphis albiventris are found throughout Northeast and Central Brazil to central-southern Uruguay and it was subject of few studies in a population level. Given this, the present study investigated the genetic variability of the species using the mitochondrial molecular marker cytochrome oxidase c subunit I. We analyzed samples from the different biomes within three Brazilian regions: Northeast (Caatinga , Cerrado, and Atlantic Forest), Southeast (Cerrado , Atlantic Forest, Cerrado/Atlantic Forest, and Cerrado/Caatinga ecotones) and South (Pampa and Atlantic Forest). Software BAPs retrieved five distinct demes: dm 1, dm 2, and dm 5 that occurs in South, Northeast and Southeast regions respectively and the dm 3 and dm 4 are wide distributed in Northeast and Southeast. Population analysis performed with AMOVA, haplotype network and Mantel test estimated the veracity of the demes. The FST shows structuring for the five demes, with dm 1 (South region) isolated from the others, however the other analysis showed the Northeast/Southeast demes (dm 2-5) united, diagnosing gene flow between them, mainly at the transitional zones, in areas as far away as areas with similar latitude interval (Southeast vs South) that was not detected gene flow. In the haplotype network, the mutational steps was conclusive in split dm1 from dm 2-5 with 15 mutational steps and the Mantel test was moderated, which is explained by genetic similarity despite the great geographic distances (Northeast/Southeast). Thus, our analysis recognized two different lineages (South and Northeast/Southeast) and indicate that the biomes were not decisive in their isolation. The sharing of demes at the transitional zones and in areas with high latitudinal intervals reflects a recent ancestral polymorphism for D. albiventris. The plasticity in the occupation of the space by this species contributes in its wide dispersion capability, that is, geographical distribution. Our results revealed important implications for the management of D. albiventris in these transitional zones areas where demes were shared.

Transcriptomic analysis of skin pigmentation variation in the Virginia opossum (Didelphis virginiana).

Skin pigmentation and coat pigmentation are two of the best-studied examples of traits under natural selection given their quantifiable fitness interactions with the environment (e.g., camouflage) and signalling with other organisms (e.g., warning coloration). Previous morphological studies have found that skin pigmentation variation in the Virginia opossum (Didelphis virginiana) is associated with variation in precipitation and temperatures across its distribution range following Gloger's rule (lighter pigmentation in temperate environments). To investigate the molecular mechanism associated with skin pigmentation variation, we used RNA-Seq and quantified gene expression of wild opossums from tropical and temperate populations. Using differential expression analysis and a co-expression network approach, we found that expression variation in genes with melanocytic and immune functions is significantly associated with the degree of skin pigmentation variation and may be underlying this phenotypic difference. We also found evidence suggesting that the Wnt/ß-catenin signalling pathway might be regulating the depigmentation observed in temperate populations. Based on our study results, we present several alternative hypotheses that may explain Gloger's rule pattern of skin pigmentation variation in opossum, including changes in pathogen diversity supporting a pathogen-resistant hypothesis, thermal stress associated with temperate environments, and pleiotropic and epistatic interactions between melanocytic and immune genes.

Variations of chromosomal structures in Caluromys philander (Didelphimorphia: Didelphidae) from the Amazon region.

Caluromys is considered to be one of the most ancient genera of extant marsupials and is positioned among the basal taxa of the family Didelphidae. At least two species occur in Brazil, C. philander and C. lanatus, both of which have 2n = 14 chromosomes. For the first time, we present evidence of an intrapopulation polymorphism of the sexual chromosome pair in C. philander females from the Central Amazon region. Detailed cytogenetic results of animals from three localities on the Amazon region were analyzed using classical cytogenetics (NOR, C-Band and G-Band) and molecular techniques (18S rDNA and telomere probes). Similar to other conspecific individuals, the diploid number of these animals is 2n = 14, and their fundamental number is 24, with NOR present on the 6th autosomal pair. The X chromosome presented variation detectable by G banding, suggesting a pericentric inversion.

Molecular characterization of an opossum Didelphis albiventris (Marsupialia: Didelphidae) population in an urban fragment of the Brazilian Atlantic rainforest and support to species barcode identification.

We made a molecular study of 40 opossums, Didelphis albiventris, from an urban fragment of the Atlantic Rainforest in southeastern Brazil, analyzing a 653-bp sequence of cytochrome c oxidase, subunit I. We found three close connected haplotypes, with low nucleotide diversity and a haplotype diversity of 59.1% and confirmed sympatry between D. albiventris and D. aurita in this region. The clear phylogenetic separation shows the appropriateness of DNA barcode identification methodology for effectively discriminating between these opossum species.

DNA barcodes effectively identify the morphologically similar Common Opossum (Didelphis marsupialis) and Virginia Opossum (Didelphis virginiana) from areas of sympatry in Mexico.

Two morphologically similar species of opossum from the genus Didelphis-Didelphis virginiana and Didelphis marsupialis-cooccur sympatrically in Mexico. High intraspecific variation complicates their morphological discrimination, under both field and museum conditions. This study aims to evaluate the utility and reliability of using DNA barcodes (short standardized genome fragments used for DNA-based identification) to distinguish these two species. Sequences of the cytochrome c oxidase subunit I (Cox1) mitochondrial gene were obtained from 12 D. marsupialis and 29 D. virginiana individuals and were compared using the neighbor-joining (NJ) algorithm with Kimura's two-parameter (K2P) model of nucleotide substitution. Average K2P distances were 1.56% within D. virginiana and 1.65% in D. marsupialis. Interspecific distances between D. virginiana and D. marsupialis varied from 7.8 to 9.3% and their barcode sequences formed distinct non-overlapping clusters on NJ trees. All sympatric specimens of both species were effectively discriminated, confirming the utility of Cox1 barcoding as a tool for taxonomic identification of these morphologically similar taxa.

A genetic analysis of the Virginia opossum mating system: evidence of multiple paternity in a highly fragmented landscape.

Using molecular techniques, we examined patterns of paternity in Virginia opossums occupying a highly fragmented agricultural landscape in northern Indiana. During 2008, we collected tissue from 64 females and their pouch young in 34 forest patches distributed over a 1100-km(2) region. Using genotypes from 10 microsatellite loci, we determined the minimum number of fathers contributing to each litter using GERUD 1.0. Genotyped offspring with known mothers were then analyzed using CERVUS 3.0, incorporating genotypes from 317 males sampled from 2007-2008 to identify potential fathers. Our analyses revealed that promiscuity was common among females, with 26 (41%) litters having > or = 2 sires. Despite the fact that we intensively sampled forest patches for potential fathers, we only were able to identify 13 fathers contributing to 14 litters, with an average Euclidean distance of 18.7 km between father-offspring pairs found in disparate patches (N = 6). Our inability to identify most (85%) fathers of sampled litters, coupled with the extensive distances observed between putative father-offspring pairs, suggests that opossums may not maintain explicit home ranges in highly fragmented landscapes.

Conformational plasticity of DM43, a metalloproteinase inhibitor from Didelphis marsupialis: chemical and pressure-induced equilibrium (un)folding studies.

We have investigated the folding of DM43, a homodimeric metalloproteinase inhibitor isolated from the serum of the South American opossum Didelphis marsupialis. Denaturation of the protein induced by GdnHCl (guanidine hydrochloride) was monitored by extrinsic and intrinsic fluorescence spectroscopy. While the equilibrium (un)folding of DM43 followed by tryptophan fluorescence was well described by a cooperative two-state transition, bis-ANS (4,4'-dianilino-1,1'-binaphthyl-5,5'-disulfonic acid) fluorescence measurements revealed an intensity maximum at the midpoint of the unfolding transition (2 M GdnHCl), indicating a partially folded intermediate state. We further investigated the DM43 intermediate stabilized at 2 M GdnHCl using size exclusion chromatography. This analysis revealed that the folding intermediate can be best described as partially folded DM43 monomers. Thermodynamic analysis of the GdnHCl-induced denaturation of DM43 revealed Gibbs free-energy changes of 13.57 kcal/mol for dimer dissociation and 1.86 kcal/mol for monomer unfolding, pointing to a critical role of dimerization as a determinant of the structure and stability of this protein. In addition, by using hydrostatic pressure (up to 3.5 kbar) we were able to stabilize partially folded states different from those stabilized in the presence of GdnHCl. Taken together, these results indicate that the conformational plasticity of DM43 could provide this protein with the ability to adapt its conformation to a variety of different environments and biological partners during its biological lifetime.

Genes flanking Xist in mouse and human are separated on the X chromosome in American marsupials.

X inactivation, the transcriptional silencing of one of the two X chromosomes in female mammals, achieves dosage compensation of X-linked genes relative to XY males. In eutherian mammals X inactivation is regulated by the X-inactive specific transcript (Xist), a cis-acting non-coding RNA that triggers silencing of the chromosome from which it is transcribed. Marsupial mammals also undergo X inactivation but the mechanism is relatively poorly understood. We set out to analyse the X chromosome in Monodelphis domestica and Didelphis virginiana, focusing on characterizing the interval defined by the Chic1 and Slc16a2 genes that in eutherians flank the Xist locus. The synteny of this region is retained on chicken chromosome 4 where other loci belonging to the evolutionarily ancient stratum of the human X chromosome, the so-called X conserved region (XCR), are also located. We show that in both M. domestica and D. virginiana an evolutionary breakpoint has separated the Chic1 and Slc16a2 loci. Detailed analysis of opossum genomic sequences revealed linkage of Chic1 with the Lnx3 gene, recently proposed to be the evolutionary precursor of Xist, and Fip1, the evolutionary precursor of Tsx, a gene located immediately downstream of Xist in eutherians. We discuss these findings in relation to the evolution of Xist and X inactivation in mammals.

Marsupials from space: fluctuating asymmetry, geographical information systems and animal conservation.

We report the development of a new quantitative method of assessing the effects of anthropogenic impacts on living beings; this method allows us to assess actual impacts and to travel backwards in time to assess impacts. In this method, we have crossed data on fluctuating asymmetry (FA, a measure of environmental or genetic stress), using Didelphis albiventris as a model, with geographical information systems data relating to environmental composition. Our results show that more impacted environments resulted in statistically higher levels of FA. Our method appears to be a useful and flexible conservation tool for assessing anthropogenic impacts.

Different patterns of selection on the nuclear genes IRBP and DMP-1 affect the efficiency but not the outcome of phylogeny estimation for didelphid marsupials.

Selection at the protein-level can influence nucleotide substitution patterns for protein-coding genes, which in turn can affect their performance as phylogenetic characters. In this study, we compare two protein-coding nuclear genes that appear to have evolved under markedly different selective constraints and evaluate how selection has shaped their phylogenetic signal. We sequenced 1,100+ bp of exon 6 of the gene encoding dentin matrix protein 1 (DMP1) from most of the currently recognized genera of New World opossums (family: Didelphidae) and compared these data to an existing matrix of sequences from the interphotoreceptor retinoid-binding protein gene (IRBP) and morphological characters. In comparison to IRBP, DMP1 has far fewer sites under strong purifying selection and exhibits a number of sites under positive directional selection. Furthermore, selection on the DMP1 protein appears to conserve short, acidic, serine-rich domains rather than primary amino acid sequence; as a result, DMP1 has significantly different nucleotide substitution patterns from IRBP. Using Bayesian methods, we determined that DMP1 evolves almost 30% faster than IRBP, has 2.5 times more variable sites, has less among-site rate heterogeneity, is skewed toward A and away from CT (IRBP has relatively even base frequencies), and has a significantly lower rate of change between adenine and any other nucleotide. Despite these different nucleotide substitution patterns, estimates of didelphid relationships based on separate phylogenetic analyses of these genes are remarkably congruent whether patterns of nucleotide substitution are explicitly modeled or not. Nonetheless, DMP1 contains more phylogenetically informative characters per unit sequence and resolves more nodes with higher support than does IRBP. Thus, for these two genes, relaxed functional constraints and positive selection appear to improve the efficiency of phylogenetic estimation without compromising its accuracy.

Fibroblast growth factor 9: cloning and immunolocalisation during tooth development in Didelphis albiventris.

There are no reports in literature about functional roles of fibroblast growth factor 9 (FGF-9) in tooth development in animals with complete tooth pattern. The classical model for studying tooth development is the mouse, which has small number of teeth and distinctive incisor and molar patterns. The opossum Didelphis albiventris with five upper and four lower incisors, one canine, three premolars, and four molars, on each side of the jaw, seems to be a convenient model to test results obtained in the mouse. Molecular expression studies indicate that FGF-9 participates in murine tooth initiation and regulation of morphogenesis. Searching for similarities and differences in FGF-9 expression between the opossum and the mouse, amino acid sequence and expression pattern of FGF-9 in the developing first molars of D. albiventris were characterised. FGF-9 cDNA sequence was obtained using RT-PCR and expressed in bacterial system for recombinant protein production and analysis of immunoreactivity. FGF-9 expression during tooth development was investigated by immunoperoxidase method. FGF-9 protein consists in a 209-residue polypeptide with a predicted molecular mass of 23.5 kDa. FGF-9 amino acid sequence has 98% of sequence identity to human and 97% to rodents. During tooth development, epithelial FGF-9 expression was seen at the dental lamina stage. Mesenchymal expression was seen at the bud stage and at the cap stage. No significant expression was found in the enamel knot. While in rodents FGF-9 is involved in initiation and regulation of tooth shape, it is suggested that it is only involved in tooth initiation in D. albiventris.

The program of sex chromosome pairing in meiosis is highly conserved across marsupial species: implications for sex chromosome evolution.

Marsupials present a series of genetic and chromosomal features that are highly conserved in very distant species. One of these features is the absence of a homologous region between X and Y chromosomes. According to this genetic differentiation, sex chromosomes do not synapse during the first meiotic prophase in males, and a special structure, the dense plate, maintains sex chromosome association. In this report we present results on the process of meiotic sex chromosome pairing obtained from three different species, Thylamys elegans, Dromiciops gliroides, and Rhyncholestes raphanurus, representing the three orders of American marsupials. We have investigated the relationships between the axial structures organized along sex chromosomes and the formation of the dense plate. We found that in the three species the dense plate arises as a modification of sex chromosomal axial elements, but without the involvement of other meiotic axial structures, such as the cohesin axes. Considering the phylogenetic relationships among the marsupials studied here, our data reinforce the idea that the dense plate emerged early in marsupial evolution as an efficient mechanism to ensure the association of the nonhomologous sex chromosomes. This situation could have influenced the further evolution of sex chromosomes in marsupials.

Molecular characterization and evolution of X and Y-borne ATRX homologues in American marsupials.

In eutherians, the sex-reversing ATRX gene on the X has no homologue on the Y chromosome. However, testis-specific and ubiquitously expressed X-borne genes have been identified in Australian marsupials. We studied nucleotide sequence and chromosomal location of ATRX homologues in two American marsupials, the opossums Didelphis virginiana and Monodelphis domestica. A PCR fragment of M. domestica ATRX was used to probe Southern blots and to screen male genomic libraries. Southern analysis demonstrated ATRX homologues on both X and Y in D. virginiana, and two clones were isolated which hybridized to a single position on the Y chromosome in male-derived cells but to multiple sites of the X in female cells. In M. domestica, there was a single clone that mapped to the X but not to the Y, suggesting that it represents the M. domestica ATRX. However a male-specific band was detected in Southern blots probed with the D. virginiana ATRY and with a mouse ATRX clone, which implies that the Y copy in M. domestica has diverged further from other ATRX homologues. Thus there appears to be a Y-borne copy of ATRY in American, as well as Australian marsupials, although it has diverged in sequence, as have other Y genes that are testis-specific in both eutherian and marsupial lineages.