Whole genome sequencing and methylome analysis of the wild guinea pig

Bank vole genomics links determinate and indeterminate growth of teeth

BACKGROUND

Continuously growing teeth are an important innovation in mammalian evolution, yet genetic regulation of continuous growth by stem cells remains incompletely understood. Dental stem cells responsible for tooth crown growth are lost at the onset of tooth root formation. Genetic signaling that initiates this loss is difficult to study with the ever-growing incisor and rooted molars of mice, the most common mammalian dental model species, because signals for root formation overlap with signals that pattern tooth size and shape (i.e., cusp patterns). Bank and prairie voles (Cricetidae, Rodentia, Glires) have evolved rooted and unrooted molars while retaining similar size and shape, providing alternative models for studying roots.

RESULTS

We assembled a de novo genome of Myodes glareolus, a vole with high-crowned, rooted molars, and performed genomic and transcriptomic analyses in a broad phylogenetic context of Glires (rodents and lagomorphs) to assess differential selection and evolution in tooth forming genes. Bulk transcriptomics comparisons of embryonic molar development between bank voles and mice demonstrated overall conservation of gene expression levels, with species-specific differences corresponding to the accelerated and more extensive patterning of the vole molar. We leverage convergent evolution of unrooted molars across the clade to examine changes that may underlie the evolution of unrooted molars. We identified 15 dental genes with changing synteny relationships and six dental genes undergoing positive selection across Glires, two of which were undergoing positive selection in species with unrooted molars, Dspp and Aqp1. Decreased expression of both genes in prairie voles with unrooted molars compared to bank voles supports the presence of positive selection and may underlie differences in root formation.

CONCLUSIONS

Our results support ongoing evolution of dental genes across Glires and identify candidate genes for mechanistic studies of root formation. Comparative research using the bank vole as a model species can reveal the complex evolutionary background of convergent evolution for ever-growing molars.

Vole genomics links determinate and indeterminate growth of teeth

Continuously growing teeth are an important innovation in mammalian evolution, yet genetic regulation of continuous growth by stem cells remains incompletely understood. Dental stem cells responsible for tooth crown growth are lost at the onset of tooth root formation. Genetic signaling that initiates this loss is difficult to study with the ever-growing incisor and rooted molars of mice, the most common mammalian dental model species, because signals for root formation overlap with signals that pattern tooth size and shape (i.e., cusp patterns). Different species of voles (Cricetidae, Rodentia, Glires) have evolved rooted and unrooted molars that have similar size and shape, providing alternative models for studying roots. We assembled a de novo genome of Myodes glareolus, a vole with high-crowned, rooted molars, and performed genomic and transcriptomic analyses in a broad phylogenetic context of Glires (rodents and lagomorphs) to assess differential selection and evolution in tooth forming genes. We identified 15 dental genes with changing synteny relationships and six dental genes undergoing positive selection across Glires, two of which were undergoing positive selection in species with unrooted molars, Dspp and Aqp1. Decreased expression of both genes in prairie voles with unrooted molars compared to bank voles supports the presence of positive selection and may underlie differences in root formation. Bulk transcriptomics analyses of embryonic molar development in bank voles also demonstrated conserved patterns of dental gene expression compared to mice, with species-specific variation likely related to developmental timing and morphological differences between mouse and vole molars. Our results support ongoing evolution of dental genes across Glires, revealing the complex evolutionary background of convergent evolution for ever-growing molars.

Integration of fluorescence in situ hybridization and chromosome-length genome assemblies revealed synteny map for guinea pig, naked mole-rat, and human

Descriptions of karyotypes of many animal species are currently available. In addition, there has been a significant increase in the number of sequenced genomes and an ever-improving quality of genome assembly. To close the gap between genomic and cytogenetic data we applied fluorescent in situ hybridization (FISH) and Hi-C technology to make the first full chromosome-level genome comparison of the guinea pig (Cavia porcellus), naked mole-rat (Heterocephalus glaber), and human. Comparative chromosome maps obtained by FISH with chromosome-specific probes link genomic scaffolds to individual chromosomes and orient them relative to centromeres and heterochromatic blocks. Hi-C assembly made it possible to close all gaps on the comparative maps and to reveal additional rearrangements that distinguish the karyotypes of the three species. As a result, we integrated the bioinformatic and cytogenetic data and adjusted the previous comparative maps and genome assemblies of the guinea pig, naked mole-rat, and human. Syntenic associations in the two hystricomorphs indicate features of their putative ancestral karyotype. We postulate that the two approaches applied in this study complement one another and provide complete information about the organization of these genomes at the chromosome level.

An ecologist's guide for studying DNA methylation variation in wild vertebrates

The field of molecular biology is advancing fast with new powerful technologies, sequencing methods and analysis software being developed constantly. Commonly used tools originally developed for research on humans and model species are now regularly used in ecological and evolutionary research. There is also a growing interest in the causes and consequences of epigenetic variation in natural populations. Studying ecological epigenetics is currently challenging, especially for vertebrate systems, because of the required technical expertise, complications with analyses and interpretation, and limitations in acquiring sufficiently high sample sizes. Importantly, neglecting the limitations of the experimental setup, technology and analyses may affect the reliability and reproducibility, and the extent to which unbiased conclusions can be drawn from these studies. Here, we provide a practical guide for researchers aiming to study DNA methylation variation in wild vertebrates. We review the technical aspects of epigenetic research, concentrating on DNA methylation using bisulfite sequencing, discuss the limitations and possible pitfalls, and how to overcome them through rigid and reproducible data analysis. This review provides a solid foundation for the proper design of epigenetic studies, a clear roadmap on the best practices for correct data analysis and a realistic view on the limitations for studying ecological epigenetics in vertebrates. This review will help researchers studying the ecological and evolutionary implications of epigenetic variation in wild populations.

Genome evolution of blind subterranean mole rats: Adaptive peripatric versus sympatric speciation

Speciation mechanisms remain controversial. Two speciation models occur in Israeli subterranean mole rats, genus Spalax: a regional speciation cline southward of four peripatric climatic chromosomal species and a local, geologic-edaphic, genic, and sympatric speciation. Here we highlight their genome evolution. The five species were separated into five genetic clusters by single nucleotide polymorphisms, copy number variations (CNVs), repeatome, and methylome in sympatry. The regional interspecific divergence correspond to Pleistocene climatic cycles. Climate warmings caused chromosomal speciation. Triple effective population size, N _e , declines match glacial cold cycles. Adaptive genes evolved under positive selection to underground stresses and to divergent climates, involving interspecies reproductive isolation. Genomic islands evolved mainly due to adaptive evolution involving ancient polymorphisms. Repeatome, including both CNV and LINE1 repetitive elements, separated the five species. Methylation in sympatry identified geologically chalk-basalt species that differentially affect thermoregulation, hypoxia, DNA repair, P53, and other pathways. Genome adaptive evolution highlights climatic and geologic-edaphic stress evolution and the two speciation models, peripatric and sympatric.

The mitochondrial genome of Cavia aperea

Cavia aperea is a wild guinea pig found throughout South America. The previously published mitochondrial sequence for C. aperea was highly divergent from the C. porcellus sequence and contained stop codons within open reading frames. Here we resequenced the mitochondrial genomes of C. aperea and C. porcellus. Both sequences reflect gene organization typical for mammalian mitochondrial DNA. Our C. aperea mtDNA sequence shows that all of the open reading frames are intact, but confirms the strikingly low level of sequence identity (92.7%) with the closely related C. porcellus mtDNA.

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.

DNA Methylation Patterns in the Round Goby Hypothalamus Support an On-The-Spot Decision Scenario for Territorial Behavior

The question as to how early life experiences are stored on a molecular level and affect traits later in life is highly topical in ecology, medicine, and epigenetics. In this study, we use a fish model to investigate whether DNA methylation mediates early life experiences and predetermines a territorial male reproductive phenotype. In fish, adult reproductive phenotypes frequently depend on previous life experiences and are often associated with distinct morphological traits. DNA methylation is an epigenetic mechanism which is both sensitive to environmental conditions and stably inherited across cell divisions. We therefore investigate early life predisposition in the round goby Neogobius melanostomus by growth back-calculations and then study DNA methylation by MBD-Seq in the brain region controlling vertebrate reproductive behavior, the hypothalamus. We find a link between the territorial reproductive phenotype and high growth rates in the first year of life. However, hypothalamic DNA methylation patterns reflect the current behavioral status independently of early life experiences. Together, our data suggest a non-predetermination scenario in the round goby, in which indeterminate males progress to a non-territorial status in the spawning season, and in which some males then assume a specialized territorial phenotype if current conditions are favorable.

Diet changes alter paternally inherited epigenetic pattern in male Wild guinea pigs

Epigenetic modifications, of which DNA methylation is the most stable, are a mechanism conveying environmental information to subsequent generations via parental germ lines. The paternal contribution to adaptive processes in the offspring might be crucial, but has been widely neglected in comparison to the maternal one. To address the paternal impact on the offspring's adaptability to changes in diet composition, we investigated if low protein diet (LPD) in F0 males caused epigenetic alterations in their subsequently sired sons. We therefore fed F0 male Wild guinea pigs with a diet lowered in protein content (LPD) and investigated DNA methylation in sons sired before and after their father's LPD treatment in both, liver and testis tissues. Our results point to a 'heritable epigenetic response' of the sons to the fathers' dietary change. Because we detected methylation changes also in the testis tissue, they are likely to be transmitted to the F2 generation. Gene-network analyses of differentially methylated genes in liver identified main metabolic pathways indicating a metabolic reprogramming ('metabolic shift'). Epigenetic mechanisms, allowing an immediate and inherited adaptation may thus be important for the survival of species in the context of a persistently changing environment, such as climate change.

Perkinsus infection is associated with alterations in the level of global DNA methylation of gills and gastrointestinal tract of the oyster Crassostrea gasar

Bivalves are filter feeders that obtain food from seawater that may contain infectious agents, such as the protozoan parasites Perkinsus marinus and P. olseni that are associated with massive mortalities responsible for losses in the aquaculture industry. Despite all physical and chemical barriers, microorganisms cross epithelia and infect host tissues to cause pathologies. Epigenetics mechanisms play important roles in a variety of human processes, from embryonic development to cell differentiation and growth. It is currently emerging as crucial mechanism involved in modulation of host-parasite interactions and pathogenesis, promoting discovery of targets for drug treatment. In bivalves, little is known about epigenetic mechanism in host parasite interactions. The objective of the present study was to evaluate the effect of Perkinsus sp. infections on DNA methylation levels in tissues of Crassostrea gasar oysters. Samples were collected in 2015 and 2016 in the Mamanguape River estuary (PB). Oyster gills were removed and used for Perkinsus sp.

DIAGNOSIS

Gills (G) and gastrointestinal tract (GT), as well as cultured P. marinus trophozoites were preserved in 95% ethanol for DNA extractions. DNA methylation levels were estimated from G and GT tissues of uninfected (n=60) and infected oysters (n=60), and from P. marinus trophozoites, by ELISA assays. Results showed that the mean prevalence of Perkinsus sp. infections was high (87.3%) in 2015 and moderate (59.6%) in 2016. DNA methylation levels of G and GT tissues were significantly lower in infected oyster than in uninfected oysters, suggesting that infections are associated with hypomethylation. Methylation level was significantly higher in G than in GT tissues, indicating a likely tissue-specific mechanism. P. marinus trophozoites showed 33% methylation. This was the first study that confirms alterations of DNA methylation in two tissues of C. gasar oysters in association with Perkinsus sp. infections.

Expression Profile of the Integrin Receptor Subunits in the Guinea Pig Sclera

PURPOSE

The ocular dimensional changes in myopia reflect increased scleral remodeling, and in high myopia, loss of scleral integrity leads to biomechanical weakening and continued scleral creep. As integrins, a type of cell surface receptors, have been linked to scleral remodeling, they represent potential targets for myopia therapies. As a first step, this study aimed to characterize the integrin subunits at the messenger RNA level in the sclera of the guinea pig, a more recently added but increasingly used animal model for myopia research.

METHODS

Primers for α and β integrin subunits were designed using NCBI/UCSC Genome Browser and Primer3 software tools. Total RNA was extracted from normal scleral tissue and isolated cultured scleral fibroblasts, as well as liver and lung, as reference tissues, all from guinea pig. cDNA was produced by reverse transcription, PCR was used to amplify products of predetermined sizes, and products were sequenced using standard methods.

RESULTS

Guinea pig scleral tissue expressed all known integrin alpha subunits except αD and αE. The latter integrin subunits were also not expressed by cultured guinea pig scleral fibroblasts; however, their expression was confirmed in guinea pig liver. In addition, isolated cultured fibroblasts did not express integrin subunits αL, αM, and αX. This difference between results for cultured cells and intact sclera presumably reflects the presence in the latter of additional cell types. Both guinea pig scleral tissue and isolated scleral fibroblasts expressed all known integrin beta subunits. All results were verified through sequencing.

CONCLUSION

The possible contributions of integrins to scleral remodeling make them plausible targets for myopia prevention. Data from this study will help guide future ex vivo and in vitro studies directed at understanding the relationship between scleral integrins and ocular growth regulation in the guinea pig model for myopia.

Efficiency of methylated DNA immunoprecipitation bisulphite sequencing for whole-genome DNA methylation analysis

AIMS

We compared four common methods for measuring DNA methylation levels and recommended the most efficient method in terms of cost and coverage.

MATERIALS & METHODS

The DNA methylation status of liver and stomach tissues was profiled using four different methods, whole-genome bisulphite sequencing (WG-BS), targeted bisulphite sequencing (Targeted-BS), methylated DNA immunoprecipitation sequencing (MeDIP-seq) and methylated DNA immunoprecipitation bisulphite sequencing (MeDIP-BS). We calculated DNA methylation levels using each method and compared the results.

RESULTS

MeDIP-BS yielded the most similar DNA methylation profile to WG-BS, with 20 times less data, suggesting remarkable cost savings and coverage efficiency compared with the other methods.

CONCLUSION

MeDIP-BS is a practical cost-effective method for analyzing whole-genome DNA methylation that is highly accurate at base-pair resolution.

Paternal heat exposure causes DNA methylation and gene expression changes of Stat3 in Wild guinea pig sons

Epigenetic mechanisms convey environmental information through generations and can regulate gene expression. Epigenetic studies in wild mammals are rare, but enable understanding adaptation processes as they may occur in nature. In most wild mammal species, males are the dispersing sex and thus often have to cope with differing habitats and thermal changes more rapidly than the often philopatric females. As temperature is a major environmental selection factor, we investigated whether genetically heterogeneous Wild guinea pig (Cavia aperea) males adapt epigenetically to an increase in temperature, whether that response will be transmitted to the next generation(s), and whether it regulates mRNA expression. Five (F0) adult male guinea pigs were exposed to an increased ambient temperature for 2 months, corresponding to the duration of the species' spermatogenesis. To study the effect of heat, we focused on the main thermoregulatory organ, the liver. We analyzed CpG‐methylation changes of male offspring (F1) sired before and after the fathers' heat treatment (as has recently been described in Weyrich et al. [Mol. Ecol., 2015]). Transcription analysis was performed for the three genes with the highest number of differentially methylated changes detected: the thermoregulation gene Signal Transducer and Activator of Transcription 3 (Stat3), the proteolytic peptidase gene Cathepsin Z (Ctsz), and Sirtuin 6 (Sirt6) with function in epigenetic regulation. Stat3 gene expression was significantly reduced (P < 0.05), which indicated a close link between CpG‐methylation and expression levels for this gene. The two other genes did not show gene expression changes. Our results indicate the presence of a paternal transgenerational epigenetic effect. Quick adaptation to climatic changes may become increasingly relevant for the survival of wildlife species as global temperatures are rising.

Differential DNA Methylation Analysis without a Reference Genome

Genome-wide DNA methylation mapping uncovers epigenetic changes associated with animal development, environmental adaptation, and species evolution. To address the lack of high-throughput methods for DNA methylation analysis in non-model organisms, we developed an integrated approach for studying DNA methylation differences independent of a reference genome. Experimentally, our method relies on an optimized 96-well protocol for reduced representation bisulfite sequencing (RRBS), which we have validated in nine species (human, mouse, rat, cow, dog, chicken, carp, sea bass, and zebrafish). Bioinformatically, we developed the RefFreeDMA software to deduce ad hoc genomes directly from RRBS reads and to pinpoint differentially methylated regions between samples or groups of individuals (http://RefFreeDMA.computational-epigenetics.org). The identified regions are interpreted using motif enrichment analysis and/or cross-mapping to annotated genomes. We validated our method by reference-free analysis of cell-type-specific DNA methylation in the blood of human, cow, and carp. In summary, we present a cost-effective method for epigenome analysis in ecology and evolution, which enables epigenome-wide association studies in natural populations and species without a reference genome.