Tissue-specific epigenetic inheritance after paternal heat exposure in male wild guinea pigs.

External temperature change has been shown to modify epigenetic patterns, such as DNA methylation, which regulates gene expression. DNA methylation is heritable, and as such provides a mechanism to convey environmental information to subsequent generations. Studies on epigenetic response to temperature increase are still scarce in wild mammals, even more so studies that compare tissue-specific epigenetic responses. Here, we aim to address differential epigenetic responses on a gene and gene pathway level in two organs, liver and testis. We chose these organs, because the liver is the main metabolic and thermoregulation organ, and epigenetic modifications in testis are potentially transmitted to the F2 generation. We focused on the transmission of DNA methylation changes to naive male offspring after paternal exposure to an ambient temperature increase of 10 °C, and investigated differential methylated regions of sons sired before and after the paternal exposure using Reduced Representation Bisulfite Sequencing. We detected both a highly tissue-specific epigenetic response, reflected in genes involved in organ-specific metabolic pathways, and a more general regulation of single genes epigenetically modified in both organs. We conclude that genomes are context-specifically differentially epigenetically regulated in response to temperature increase. These findings emphasize the epigenetic relevance in cell differentiation, which is essential for the specific function(s) of complex organs, and is represented in a diverse molecular regulation of genes and gene pathways. The results also emphasize the paternal contribution to adaptive processes.

Phylogeography of red muntjacs reveals three distinct mitochondrial lineages.

BACKGROUND: The members of the genus Muntiacus are of particular interest to evolutionary biologists due to their extreme chromosomal rearrangements and the ongoing discussions about the number of living species. Red muntjacs have the largest distribution of all muntjacs and were formerly considered as one species. Karyotype differences led to the provisional split between the Southern Red Muntjac (Muntiacus muntjak) and the Northern Red Muntjac (M. vaginalis), but uncertainties remain as, so far, no phylogenetic study has been conducted. Here, we analysed whole mitochondrial genomes of 59 archival and 16 contemporaneous samples to resolve uncertainties about their taxonomy and used red muntjacs as model for understanding the evolutionary history of other species in Southeast Asia. RESULTS: We found three distinct matrilineal groups of red muntjacs: Sri Lankan red muntjacs (including the Western Ghats) diverged first from other muntjacs about 1.5 Mya; later northern red muntjacs (including North India and Indochina) and southern red muntjacs (Sundaland) split around 1.12 Mya. The diversification of red muntjacs into these three main lineages was likely promoted by two Pleistocene barriers: one through the Indian subcontinent and one separating the Indochinese and Sundaic red muntjacs. Interestingly, we found a high level of gene flow within the populations of northern and southern red muntjacs, indicating gene flow between populations in Indochina and dispersal of red muntjacs over the exposed Sunda Shelf during the Last Glacial Maximum. CONCLUSIONS: Our results provide new insights into the evolution of species in South and Southeast Asia as we found clear genetic differentiation in a widespread and generalist species, corresponding to two known biogeographical barriers: The Isthmus of Kra and the central Indian dry zone. In addition, our molecular data support either the delineation of three monotypic species or three subspecies, but more importantly these data highlight the conservation importance of the Sri Lankan/South Indian red muntjac.

Genetic Structure and Phylogeography of the Leopard Cat (Prionailurus bengalensis) Inferred from Mitochondrial Genomes.

The Leopard cat Prionailurus bengalensis is a habitat generalist that is widely distributed across Southeast Asia. Based on morphological traits, this species has been subdivided into 12 subspecies. Thus far, there have been few molecular studies investigating intraspecific variation, and those had been limited in geographic scope. For this reason, we aimed to study the genetic structure and evolutionary history of this species across its very large distribution range in Asia. We employed both PCR-based (short mtDNA fragments, 94 samples) and high throughput sequencing based methods (whole mitochondrial genomes, 52 samples) on archival, noninvasively collected and fresh samples to investigate the distribution of intraspecific genetic variation. Our comprehensive sampling coupled with the improved resolution of a mitochondrial genome analyses provided strong support for a deep split between Mainland and Sundaic Leopard cats. Although we identified multiple haplogroups within the species' distribution, we found no matrilineal evidence for the distinction of 12 subspecies. In the context of Leopard cat biogeography, we cautiously recommend a revision of the Prionailurus bengalensis subspecific taxonomy: namely, a reduction to 4 subspecies (2 mainland and 2 Sundaic forms).

Paternal intergenerational epigenetic response to heat exposure in male Wild guinea pigs.

Epigenetic modifications, of which DNA methylation is the best studied one, can convey environmental information through generations via parental germ lines. Past studies have focused on the maternal transmission of epigenetic information to the offspring of isogenic mice and rats in response to external changes, whereas heterogeneous wild mammals as well as paternal epigenetic effects have been widely neglected. In most wild mammal species, males are the dispersing sex and have to cope with differing habitats and thermal changes. As temperature is a major environmental factor we investigated if genetically heterogeneous Wild guinea pig (Cavia aperea) males can adapt epigenetically to an increase in temperature and if that response will be transmitted to the next generation(s). Five adult male guinea pigs (F0) were exposed to an increased ambient temperature for 2 months, i.e. the duration of spermatogenesis. We studied the liver (as the main thermoregulatory organ) of F0 fathers and F1 sons, and testes of F1 sons for paternal transmission of epigenetic modifications across generation(s). Reduced representation bisulphite sequencing revealed shared differentially methylated regions in annotated areas between F0 livers before and after heat treatment, and their sons' livers and testes, which indicated a general response with ecological relevance. Thus, paternal exposure to a temporally limited increased ambient temperature led to an 'immediate' and 'heritable' epigenetic response that may even be transmitted to the F2 generation. In the context of globally rising temperatures epigenetic mechanisms may become increasingly relevant for the survival of species.

Environmental Change-Dependent Inherited Epigenetic Response.

Epigenetic modifications are a mechanism conveying environmental information to subsequent generations via parental germ lines. Research on epigenetic responses to environmental changes in wild mammals has been widely neglected, as well as studies that compare responses to changes in different environmental factors. Here, we focused on the transmission of DNA methylation changes to naive male offspring after paternal exposure to either diet (~40% less protein) or temperature increase (10 °C increased temperature). Because both experiments focused on the liver as the main metabolic and thermoregulation organ, we were able to decipher if epigenetic changes differed in response to different environmental changes. Reduced representation bisulfite sequencing (RRBS) revealed differentially methylated regions (DMRs) in annotated genomic regions in sons sired before (control) and after the fathers' treatments. We detected both a highly specific epigenetic response dependent on the environmental factor that had changed that was reflected in genes involved in specific metabolic pathways, and a more general response to changes in outer stimuli reflected by epigenetic modifications in a small subset of genes shared between both responses. Our results indicated that fathers prepared their offspring for specific environmental changes by paternally inherited epigenetic modifications, suggesting a strong paternal contribution to adaptive processes.

Human-mediated introduction of introgressed deer across Wallace's line: Historical biogeography of Rusa unicolor and R. timorensis.

In this study we compared the phylogeographic patterns of two Rusa species, Rusa unicolor and Rusa timorensis, in order to understand what drove and maintained differentiation between these two geographically and genetically close species and investigated the route of introduction of individuals to the islands outside of the Sunda Shelf. We analyzed full mitogenomes from 56 archival samples from the distribution areas of the two species and 18 microsatellite loci in a subset of 16 individuals to generate the phylogeographic patterns of both species. Bayesian inference with fossil calibration was used to estimate the age of each species and major divergence events. Our results indicated that the split between the two species took place during the Pleistocene, ~1.8 Mya, possibly driven by adaptations of R. timorensis to the drier climate found on Java compared to the other islands of Sundaland. Although both markers identified two well-differentiated clades, there was a largely discrepant pattern between mitochondrial and nuclear markers. While nDNA separated the individuals into the two species, largely in agreement with their museum label, mtDNA revealed that all R. timorensis sampled to the east of the Sunda shelf carried haplotypes from R. unicolor and one Rusa unicolor from South Sumatra carried a R. timorensis haplotype. Our results show that hybridization occurred between these two sister species in Sundaland during the Late Pleistocene and resulted in human-mediated introduction of hybrid descendants in all islands outside Sundaland.

Raritas: a program for counting high diversity categorical data with highly unequal abundances.

Acquiring data on the occurrences of many types of difficult to identify objects are often still made by human observation, for example, in biodiversity and paleontologic research. Existing computer counting programs used to record such data have various limitations, including inflexibility and cost. We describe a new open-source program for this purpose-Raritas. Raritas is written in Python and can be run as a standalone app for recent versions of either MacOS or Windows, or from the command line as easily customized source code. The program explicitly supports a rare category count mode which makes it easier to collect quantitative data on rare categories, for example, rare species which are important in biodiversity surveys. Lastly, we describe the file format used by Raritas and propose it as a standard for storing geologic biodiversity data. 'Stratigraphic occurrence data' file format combines extensive sample metadata and a flexible structure for recording occurrence data of species or other categories in a series of samples.

Targeted resequencing of coding DNA sequences for SNP discovery in nonmodel species.

Targeted capture coupled with high-throughput sequencing can be used to gain information about nuclear sequence variation at hundreds to thousands of loci. Divergent reference capture makes use of molecular data of one species to enrich target loci in other (related) species. This is particularly valuable for nonmodel organisms, for which often no a priori knowledge exists regarding these loci. Here, we have used targeted capture to obtain data for 809 nuclear coding DNA sequences (CDS) in a nonmodel organism, the Eurasian lynx Lynx lynx, using baits designed with the help of the published genome of a related model organism (the domestic cat Felis catus). Using this approach, we were able to survey intraspecific variation at hundreds of nuclear loci in L. lynx across the species' European range. A large set of biallelic candidate SNPs was then evaluated using a high-throughput SNP genotyping platform (Fluidigm), which we then reduced to a final 96 SNP-panel based on assay performance and reliability; validation was carried out with 100 additional Eurasian lynx samples not included in the SNP discovery phase. The 96 SNP-panel developed from CDS performed very successfully in the identification of individuals and in population genetic structure inference (including the assignment of individuals to their source population). In keeping with recent studies, our results show that genic SNPs can be valuable for genetic monitoring of wildlife species.

An expanded mammal mitogenome dataset from Southeast Asia.

Southeast (SE) Asia is 1 of the most biodiverse regions in the world, and it holds approximately 20% of all mammal species. Despite this, the majority of SE Asia's genetic diversity is still poorly characterized. The growing interest in using environmental DNA to assess and monitor SE Asian species, in particular threatened mammals-has created the urgent need to expand the available reference database of mitochondrial barcode and complete mitogenome sequences. We have partially addressed this need by generating 72 new mitogenome sequences reconstructed from DNA isolated from a range of historical and modern tissue samples. Approximately 55 gigabases of raw sequence were generated. From this data, we assembled 72 complete mitogenome sequences, with an average depth of coverage of ×102.9 and ×55.2 for modern samples and historical samples, respectively. This dataset represents 52 species, of which 30 species had no previous mitogenome data available. The mitogenomes were geotagged to their sampling location, where known, to display a detailed geographical distribution of the species. Our new database of 52 taxa will strongly enhance the utility of environmental DNA approaches for monitoring mammals in SE Asia as it greatly increases the likelihoods that identification of metabarcoding sequencing reads can be assigned to reference sequences. This magnifies the confidence in species detections and thus allows more robust surveys and monitoring programmes of SE Asia's threatened mammal biodiversity. The extensive collections of historical samples from SE Asia in western and SE Asian museums should serve as additional valuable material to further enrich this reference database.

Evolutionary Relationships among Extinct and Extant Sloths: The Evidence of Mitogenomes and Retroviruses.

Macroevolutionary trends exhibited by retroviruses are complex and not entirely understood. The sloth endogenized foamy-like retrovirus (SloEFV), which demonstrates incongruence in virus-host evolution among extant sloths (Order Folivora), has not been investigated heretofore in any extinct sloth lineages and its premodern history within folivorans is therefore unknown. Determining retroviral coevolutionary trends requires a robust phylogeny of the viral host, but the highly reduced modern sloth fauna (6 species in 2 genera) does not adequately represent what was once a highly diversified clade (∼100 genera) of placental mammals. At present, the amount of molecular data available for extinct sloth taxa is limited, and analytical results based on these data tend to conflict with phylogenetic inferences made on the basis of morphological studies. To augment the molecular data set, we applied hybridization capture and next-generation Illumina sequencing to two extinct and three extant sloth species to retrieve full mitochondrial genomes (mitogenomes) from the hosts and the polymerase gene of SloEFV. The results produced a fully resolved and well-supported phylogeny that supports dividing crown families into two major clades: 1) The three-toed sloth, Bradypus, and Nothrotheriidae and 2) Megalonychidae, including the two-toed sloth, Choloepus, and Mylodontidae. Our calibrated time tree indicates that the Miocene epoch (23.5 Ma), particularly its earlier part, was an important interval for folivoran diversification. Both extant and extinct sloths demonstrate multiple complex invasions of SloEFV into the ancestral sloth germline followed by subsequent introgressions across different sloth lineages. Thus, sloth mitogenome and SloEFV evolution occurred separately and in parallel among sloths.

Two species of Southeast Asian cats in the genus Catopuma with diverging histories: an island endemic forest specialist and a widespread habitat generalist.

Background. The bay cat Catopuma badia is endemic to Borneo, whereas its sister species the Asian golden cat Catopuma temminckii is distributed from the Himalayas and southern China through Indochina, Peninsular Malaysia and Sumatra. Based on morphological data, up to five subspecies of the Asian golden cat have been recognized, but a taxonomic assessment, including molecular data and morphological characters, is still lacking. Results. We combined molecular data (whole mitochondrial genomes), morphological data (pelage) and species distribution projections (up to the Late Pleistocene) to infer how environmental changes may have influenced the distribution of these sister species over the past 120 000 years. The molecular analysis was based on sequenced mitogenomes of 3 bay cats and 40 Asian golden cats derived mainly from archival samples. Our molecular data suggested a time of split between the two species approximately 3.16 Ma and revealed very low nucleotide diversity within the Asian golden cat population, which supports recent expansion of the population. Discussion. The low nucleotide diversity suggested a population bottleneck in the Asian golden cat, possibly caused by the eruption of the Toba volcano in Northern Sumatra (approx. 74 kya), followed by a continuous population expansion in the Late Pleistocene/Early Holocene. Species distribution projections, the reconstruction of the demographic history, a genetic isolation-by-distance pattern and a gradual variation of pelage pattern support the hypothesis of a post-Toba population expansion of the Asian golden cat from south China/Indochina to Peninsular Malaysia and Sumatra. Our findings reject the current classification of five subspecies for the Asian golden cat, but instead support either a monotypic species or one comprising two subspecies: (i) the Sunda golden cat, distributed south of the Isthmus of Kra: C. t. temminckii and (ii) Indochinese, Indian, Himalayan and Chinese golden cats, occurring north of the Isthmus: C. t. moormensis.

Comprehensive profiling of retroviral integration sites using target enrichment methods from historical koala samples without an assembled reference genome.

Background. Retroviral integration into the host germline results in permanent viral colonization of vertebrate genomes. The koala retrovirus (KoRV) is currently invading the germline of the koala (Phascolarctos cinereus) and provides a unique opportunity for studying retroviral endogenization. Previous analysis of KoRV integration patterns in modern koalas demonstrate that they share integration sites primarily if they are related, indicating that the process is currently driven by vertical transmission rather than infection. However, due to methodological challenges, KoRV integrations have not been comprehensively characterized. Results. To overcome these challenges, we applied and compared three target enrichment techniques coupled with next generation sequencing (NGS) and a newly customized sequence-clustering based computational pipeline to determine the integration sites for 10 museum Queensland and New South Wales (NSW) koala samples collected between the 1870s and late 1980s. A secondary aim of this study sought to identify common integration sites across modern and historical specimens by comparing our dataset to previously published studies. Several million sequences were processed, and the KoRV integration sites in each koala were characterized. Conclusions. Although the three enrichment methods each exhibited bias in integration site retrieval, a combination of two methods, Primer Extension Capture and hybridization capture is recommended for future studies on historical samples. Moreover, identification of integration sites shows that the proportion of integration sites shared between any two koalas is quite small.

Comparative analysis of miRNAs and their targets across four plant species.

BACKGROUND: MicroRNA (miRNA) mediated regulation of gene expression has been recognized as a major posttranscriptional regulatory mechanism also in plants. We performed a comparative analysis of miRNAs and their respective gene targets across four plant species: Arabidopsis thaliana (Ath), Medicago truncatula(Mtr), Brassica napus (Bna), and Chlamydomonas reinhardtii (Cre). RESULTS: miRNAs were obtained from mirBase with 218 miRNAs for Ath, 375 for Mtr, 46 for Bna, and 73 for Cre, annotated for each species respectively. miRNA targets were obtained from available database annotations, bioinformatic predictions using RNAhybrid as well as predicted from an analysis of mRNA degradation products (degradome sequencing) aimed at identifying miRNA cleavage products. On average, and considering both experimental and bioinformatic predictions together, every miRNA was associated with about 46 unique gene transcripts with considerably variation across species. We observed a positive and linear correlation between the number miRNAs and the total number of transcripts across different plant species suggesting that the repertoire of miRNAs correlates with the size of the transcriptome of an organism. Conserved miRNA-target pairs were found to be associated with developmental processes and transcriptional regulation, while species-specific (in particular, Ath) pairs are involved in signal transduction and response to stress processes. Conserved miRNAs have more targets and higher expression values than non-conserved miRNAs. We found evidence for a conservation of not only the sequence of miRNAs, but their expression levels as well. CONCLUSIONS: Our results support the notion of a high birth and death rate of miRNAs and that miRNAs serve many species specific functions, while conserved miRNA are related mainly to developmental processes and transcriptional regulation with conservation operating at both the sequence and expression level.