Comparative genomic analysis of the human genome and six bat genomes using unsupervised machine learning: Mb-level CpG and TFBS islands.

BACKGROUND: Emerging infectious disease-causing RNA viruses, such as the SARS-CoV-2 and Ebola viruses, are thought to rely on bats as natural reservoir hosts. Since these zoonotic viruses pose a great threat to humans, it is important to characterize the bat genome from multiple perspectives. Unsupervised machine learning methods for extracting novel information from big sequence data without prior knowledge or particular models are highly desirable for obtaining unexpected insights. We previously established a batch-learning self-organizing map (BLSOM) of the oligonucleotide composition that reveals novel genome characteristics from big sequence data. RESULTS: In this study, using the oligonucleotide BLSOM, we conducted a comparative genomic study of humans and six bat species. BLSOM is an explainable-type machine learning algorithm that reveals the diagnostic oligonucleotides contributing to sequence clustering (self-organization). When unsupervised machine learning reveals unexpected and/or characteristic features, these features can be studied in more detail via the much simpler and more direct standard distribution map method. Based on this combined strategy, we identified the Mb-level enrichment of CG dinucleotide (Mb-level CpG islands) around the termini of bat long-scaffold sequences. In addition, a class of CG-containing oligonucleotides were enriched in the centromeric and pericentromeric regions of human chromosomes. Oligonucleotides longer than tetranucleotides often represent binding motifs for a wide variety of proteins (e.g., transcription factor binding sequences (TFBSs)). By analyzing the penta- and hexanucleotide composition, we observed the evident enrichment of a wide range of hexanucleotide TFBSs in centromeric and pericentromeric heterochromatin regions on all human chromosomes. CONCLUSION: Function of transcription factors (TFs) beyond their known regulation of gene expression (e.g., TF-mediated looping interactions between two different genomic regions) has received wide attention. The Mb-level TFBS and CpG islands are thought to be involved in the large-scale nuclear organization, such as centromere and telomere clustering. TFBSs, which are enriched in centromeric and pericentromeric heterochromatin regions, are thought to play an important role in the formation of nuclear 3D structures. Our machine learning-based analysis will help us to understand the differential features of nuclear 3D structures in the human and bat genomes.

Comparison of structural variants in the whole genome sequences of two Medicago truncatula ecotypes: Jemalong A17 and R108.

BACKGROUND: Structural variants (SVs) constitute a large proportion of the genomic variation that results in phenotypic variation in plants. However, they are still a largely unexplored feature in most plant genomes. Here, we present the whole-genome landscape of SVs between two model legume Medicago truncatula ecotypes-Jemalong A17 and R108- that have been extensively used in various legume biology studies. RESULTS: To catalogue SVs, we first resolved the previously published R108 genome assembly (R108 v1.0) to chromosome-scale using 124 × Hi-C data, resulting in a high-quality genome assembly. The inter-chromosomal reciprocal translocations between chromosomes 4 and 8 were confirmed by performing syntenic analysis between the two genomes. Combined with the Hi-C data, it appears that these translocation events had a significant effect on chromatin organization. Using both whole-genome and short-read alignments, we identified the genomic landscape of SVs between the two genomes, some of which may account for several phenotypic differences, including their differential responses to aluminum toxicity and iron deficiency, and the development of different anthocyanin leaf markings. We also found extensive SVs within the nodule-specific cysteine-rich gene family which encodes antimicrobial peptides essential for terminal bacteroid differentiation during nitrogen-fixing symbiosis. CONCLUSIONS: Our results provide a near-complete R108 genome assembly and the first genomic landscape of SVs obtained by comparing two M. truncatula ecotypes. This may provide valuable genomic resources for the functional and molecular research of legume biology in the future.

Three classes of epigenomic regulators converge to hyperactivate the essential maternal gene deadhead within a heterochromatin mini-domain.

The formation of a diploid zygote is a highly complex cellular process that is entirely controlled by maternal gene products stored in the egg cytoplasm. This highly specialized transcriptional program is tightly controlled at the chromatin level in the female germline. As an extreme case in point, the massive and specific ovarian expression of the essential thioredoxin Deadhead (DHD) is critically regulated in Drosophila by the histone demethylase Lid and its partner, the histone deacetylase complex Sin3A/Rpd3, via yet unknown mechanisms. Here, we identified Snr1 and Mod(mdg4) as essential for dhd expression and investigated how these epigenomic effectors act with Lid and Sin3A to hyperactivate dhd. Using Cut&Run chromatin profiling with a dedicated data analysis procedure, we found that dhd is intriguingly embedded in an H3K27me3/H3K9me3-enriched mini-domain flanked by DNA regulatory elements, including a dhd promoter-proximal element essential for its expression. Surprisingly, Lid, Sin3a, Snr1 and Mod(mdg4) impact H3K27me3 and this regulatory element in distinct manners. However, we show that these effectors activate dhd independently of H3K27me3/H3K9me3, and that dhd remains silent in the absence of these marks. Together, our study demonstrates an atypical and critical role for chromatin regulators Lid, Sin3A, Snr1 and Mod(mdg4) to trigger tissue-specific hyperactivation within a unique heterochromatin mini-domain.

Heterochromatin is a quantitative trait associated with spontaneous epiallele formation.

Epialleles are meiotically heritable variations in expression states that are independent from changes in DNA sequence. Although they are common in plant genomes, their molecular origins are unknown. Here we show, using mutant and experimental populations, that epialleles in Arabidopsis thaliana that result from ectopic hypermethylation are due to feedback regulation of pathways that primarily function to maintain DNA methylation at heterochromatin. Perturbations to maintenance of heterochromatin methylation leads to feedback regulation of DNA methylation in genes. Using single base resolution methylomes from epigenetic recombinant inbred lines (epiRIL), we show that epiallelic variation is abundant in euchromatin, yet, associates with QTL primarily in heterochromatin regions. Mapping three-dimensional chromatin contacts shows that genes that are hotspots for ectopic hypermethylation have increases in contact frequencies with regions possessing H3K9me2. Altogether, these data show that feedback regulation of pathways that have evolved to maintain heterochromatin silencing leads to the origins of spontaneous hypermethylated epialleles.

Three-dimensional virtual histology of the human hippocampus based on phase-contrast computed tomography.

We have studied the three-dimensional (3D) cytoarchitecture of the human hippocampus in neuropathologically healthy and Alzheimer's disease (AD) individuals, based on phase-contrast X-ray computed tomography of postmortem human tissue punch biopsies. In view of recent findings suggesting a nuclear origin of AD, we target in particular the nuclear structure of the dentate gyrus (DG) granule cells. Tissue samples of 20 individuals were scanned and evaluated using a highly automated approach of measurement and analysis, combining multiscale recordings, optimized phase retrieval, segmentation by machine learning, representation of structural properties in a feature space, and classification based on the theory of optimal transport. Accordingly, we find that the prototypical transformation between a structure representing healthy granule cells and the pathological state involves a decrease in the volume of granule cell nuclei, as well as an increase in the electron density and its spatial heterogeneity. The latter can be explained by a higher ratio of heterochromatin to euchromatin. Similarly, many other structural properties can be derived from the data, reflecting both the natural polydispersity of the hippocampal cytoarchitecture between different individuals in the physiological context and the structural effects associated with AD pathology.

DNA polymerase zeta contributes to heterochromatin replication to prevent genome instability.

The DNA polymerase zeta (Polζ) plays a critical role in bypassing DNA damage. REV3L, the catalytic subunit of Polζ, is also essential in mouse embryonic development and cell proliferation for reasons that remain incompletely understood. In this study, we reveal that REV3L protein interacts with heterochromatin components including repressive histone marks and localizes in pericentromeric regions through direct interaction with HP1 dimer. We demonstrate that Polζ/REV3L ensures progression of replication forks through difficult-to-replicate pericentromeric heterochromatin, thereby preventing spontaneous chromosome break formation. We also find that Rev3l-deficient cells are compromised in the repair of heterochromatin-associated double-stranded breaks, eliciting deletions in late-replicating regions. Lack of REV3L leads to further consequences that may be ascribed to heterochromatin replication and repair-associated functions of Polζ, with a disruption of the temporal replication program at specific loci. This is correlated with changes in epigenetic landscape and transcriptional control of developmentally regulated genes. These results reveal a new function of Polζ in preventing chromosome instability during replication of heterochromatic regions.

SUV39 SET domains mediate crosstalk of heterochromatic histone marks.

The SUV39 class of methyltransferase enzymes deposits histone H3 lysine 9 di- and trimethylation (H3K9me2/3), the hallmark of constitutive heterochromatin. How these enzymes are regulated to mark specific genomic regions as heterochromatic is poorly understood. Clr4 is the sole H3K9me2/3 methyltransferase in the fission yeast Schizosaccharomyces pombe, and recent evidence suggests that ubiquitination of lysine 14 on histone H3 (H3K14ub) plays a key role in H3K9 methylation. However, the molecular mechanism of this regulation and its role in heterochromatin formation remain to be determined. Our structure-function approach shows that the H3K14ub substrate binds specifically and tightly to the catalytic domain of Clr4, and thereby stimulates the enzyme by over 250-fold. Mutations that disrupt this mechanism lead to a loss of H3K9me2/3 and abolish heterochromatin silencing similar to clr4 deletion. Comparison with mammalian SET domain proteins suggests that the Clr4 SET domain harbors a conserved sensor for H3K14ub, which mediates licensing of heterochromatin formation.

The exploration of N6-deoxyadenosine methylation in mammalian genomes.

N6-methyladenine (N6-mA, m6dA, or 6mA), a prevalent DNA modification in prokaryotes, has recently been identified in higher eukaryotes, including mammals. Although 6mA has been well-studied in prokaryotes, the function and regulatory mechanism of 6mA in eukaryotes are still poorly understood. Recent studies indicate that 6mA can serve as an epigenetic mark and play critical roles in various biological processes, from transposable-element suppression to environmental stress response. Here, we review the significant advances in methodology for 6mA detection and major progress in understanding the regulation and function of this non-canonical DNA methylation in eukaryotes, predominantly mammals.

Genome-wide profiles of UV lesion susceptibility, repair, and mutagenic potential in melanoma.

Exposure to the ultraviolet (UV) radiation in sunlight creates DNA lesions, which if left unrepaired can induce mutations and contribute to skin cancer. The two most common UV-induced DNA lesions are the cis-syn cyclobutane pyrimidine dimers (CPDs) and pyrimidine (6-4) pyrimidone photoproducts (6-4PPs), both of which can initiate mutations. Interestingly, mutation frequency across the genomes of many cancers is heterogenous with significant increases in heterochromatin. Corresponding increases in UV lesion susceptibility and decreases in repair are observed in heterochromatin versus euchromatin. However, the individual contributions of CPDs and 6-4PPs to mutagenesis have not been systematically examined in specific genomic and epigenomic contexts. In this study, we compared genome-wide maps of 6-4PP and CPD lesion abundances in primary cells and conducted comprehensive analyses to determine the genetic and epigenetic features associated with susceptibility. Overall, we found a high degree of similarity between 6-4PP and CPD formation, with an enrichment of both in heterochromatin regions. However, when examining the relative levels of the two UV lesions, we found that bivalent and Polycomb-repressed chromatin states were uniquely more susceptible to 6-4PPs. Interestingly, when comparing UV susceptibility and repair with melanoma mutation frequency in these regions, disparate patterns were observed in that susceptibility was not always inversely associated with repair and mutation frequency. Functional enrichment analysis hint at mechanisms of negative selection for these regions that are essential for cell viability, immune function and induce cell death when mutated. Ultimately, these results reveal both the similarities and differences between UV-induced lesions that contribute to melanoma.

Telomere-to-telomere assembly of a fish Y chromosome reveals the origin of a young sex chromosome pair.

BACKGROUND: The origin of sex chromosomes requires the establishment of recombination suppression between the proto-sex chromosomes. In many fish species, the sex chromosome pair is homomorphic with a recent origin, providing species for studying how and why recombination suppression evolved in the initial stages of sex chromosome differentiation, but this requires accurate sequence assembly of the X and Y (or Z and W) chromosomes, which may be difficult if they are recently diverged. RESULTS: Here we produce a haplotype-resolved genome assembly of zig-zag eel (Mastacembelus armatus), an aquaculture fish, at the chromosomal scale. The diploid assembly is nearly gap-free, and in most chromosomes, we resolve the centromeric and subtelomeric heterochromatic sequences. In particular, the Y chromosome, including its highly repetitive short arm, has zero gaps. Using resequencing data, we identify a ~7 Mb fully sex-linked region (SLR), spanning the sex chromosome centromere and almost entirely embedded in the pericentromeric heterochromatin. The SLRs on the X and Y chromosomes are almost identical in sequence and gene content, but both are repetitive and heterochromatic, consistent with zero or low recombination. We further identify an HMG-domain containing gene HMGN6 in the SLR as a candidate sex-determining gene that is expressed at the onset of testis development. CONCLUSIONS: Our study supports the idea that preexisting regions of low recombination, such as pericentromeric regions, can give rise to SLR in the absence of structural variations between the proto-sex chromosomes.

Elevated retrotransposon activity and genomic instability in primed pluripotent stem cells.

BACKGROUND: Naïve and primed pluripotent stem cells (PSCs) represent two different pluripotent states. Primed PSCs following in vitro culture exhibit lower developmental potency as evidenced by failure in germline chimera assays, unlike mouse naïve PSCs. However, the molecular mechanisms underlying the lower developmental competency of primed PSCs remain elusive. RESULTS: We examine the regulation of telomere maintenance, retrotransposon activity, and genomic stability of primed PSCs and compare them with naïve PSCs. Surprisingly, primed PSCs only minimally maintain telomeres and show fragile telomeres, associated with declined DNA recombination and repair activity, in contrast to naïve PSCs that robustly elongate telomeres. Also, we identify LINE1 family integrant L1Md_T as naïve-specific retrotransposon and ERVK family integrant IAPEz to define primed PSCs, and their transcription is differentially regulated by heterochromatic histones and Dnmt3b. Notably, genomic instability of primed PSCs is increased, in association with aberrant retrotransposon activity. CONCLUSIONS: Our data suggest that fragile telomere, retrotransposon-associated genomic instability, and declined DNA recombination repair, together with reduced function of cell cycle and mitochondria, increased apoptosis, and differentiation properties may link to compromised developmental potency of primed PSCs, noticeably distinguishable from naïve PSCs.

Quadruplex Folding Promotes the Condensation of Linker Histones and DNAs via Liquid-Liquid Phase Separation.

Liquid-liquid phase separation (LLPS) of proteins and DNA has recently emerged as a possible mechanism underlying the dynamic organization of chromatin. We herein report the role of DNA quadruplex folding in liquid droplet formation via LLPS induced by interactions between DNA and linker histone H1 (H1), a key regulator of chromatin organization. Fluidity measurements inside the droplets, binding assays using G-quadruplex-selective probes, and structural analyses based on circular dichroism demonstrated that quadruplex DNA structures, such as the G-quadruplex and i-motif, promote droplet formation with H1 and decrease molecular motility within droplets. The dissolution of the droplets in the presence of additives and the LLPS of the DNA structural units indicated that, in addition to electrostatic interactions between the DNA and the intrinsically disordered region of H1, π-π stacking between quadruplex DNAs could potentially drive droplet formation, unlike in the electrostatically driven LLPS of duplex DNA and H1. According to phase diagrams of anionic molecules with various conformations, the high LLPS ability associated with quadruplex folding arises from the formation of interfaces consisting of organized planes of guanine bases and the side surfaces with a high charge density. Given that DNA quadruplex structures are well-documented in heterochromatin regions, it is imperative to understand the role of DNA quadruplex folding in the context of intranuclear LLPS.

An enhancer screen identifies new suppressors of small-RNA-mediated epigenetic gene silencing.

Small non-protein coding RNAs are involved in pathways that control the genome at the level of chromatin. In Schizosaccharomyces pombe, small interfering RNAs (siRNAs) are required for the faithful propagation of heterochromatin that is found at peri-centromeric repeats. In contrast to repetitive DNA, protein-coding genes are refractory to siRNA-mediated heterochromatin formation, unless siRNAs are expressed in mutant cells. Here we report the identification of 20 novel mutant alleles that enable de novo formation of heterochromatin at a euchromatic protein-coding gene by using trans-acting siRNAs as triggers. For example, a single amino acid substitution in the pre-mRNA cleavage factor Yth1 enables siRNAs to trigger silent chromatin formation with unparalleled efficiency. Our results are consistent with a kinetic nascent transcript processing model for the inhibition of small-RNA-directed de novo formation of heterochromatin and lay a foundation for further mechanistic dissection of cellular activities that counteract epigenetic gene silencing.

Establishment of H3K9me3-dependent heterochromatin during embryogenesis in Drosophila miranda.

Heterochromatin is a key architectural feature of eukaryotic genomes crucial for silencing of repetitive elements. During Drosophila embryonic cellularization, heterochromatin rapidly appears over repetitive sequences, but the molecular details of how heterochromatin is established are poorly understood. Here, we map the genome-wide distribution of H3K9me3-dependent heterochromatin in individual embryos of Drosophila miranda at precisely staged developmental time points. We find that canonical H3K9me3 enrichment is established prior to cellularization and matures into stable and broad heterochromatin domains through development. Intriguingly, initial nucleation sites of H3K9me3 enrichment appear as early as embryonic stage 3 over transposable elements (TEs) and progressively broaden, consistent with spreading to neighboring nucleosomes. The earliest nucleation sites are limited to specific regions of a small number of recently active retrotransposon families and often appear over promoter and 5' regions of LTR retrotransposons, while late nucleation sites develop broadly across the entirety of most TEs. Interestingly, early nucleating TEs are strongly associated with abundant maternal piRNAs and show early zygotic transcription. These results support a model of piRNA-associated co-transcriptional silencing while also suggesting additional mechanisms for site-restricted H3K9me3 nucleation at TEs in pre-cellular Drosophila embryos.

Tet1 regulates epigenetic remodeling of the pericentromeric heterochromatin and chromocenter organization in DNA hypomethylated cells.

Pericentromeric heterochromatin (PCH), the constitutive heterochromatin of pericentromeric regions, plays crucial roles in various cellular events, such as cell division and DNA replication. PCH forms chromocenters in the interphase nucleus, and chromocenters cluster at the prophase of meiosis. Chromocenter clustering has been reported to be critical for the appropriate progression of meiosis. However, the molecular mechanisms underlying chromocenter clustering remain elusive. In this study, we found that global DNA hypomethylation, 5hmC enrichment in PCH, and chromocenter clustering of Dnmt1-KO ESCs were similar to those of the female meiotic germ cells. Tet1 is essential for the deposition of 5hmC and facultative histone marks of H3K27me3 and H2AK119ub at PCH, as well as chromocenter clustering. RING1B, one of the core components of PRC1, is recruited to PCH by TET1, and PRC1 plays a critical role in chromocenter clustering. In addition, the rearrangement of the chromocenter under DNA hypomethylated condition was mediated by liquid-liquid phase separation. Thus, we demonstrated a novel role of Tet1 in chromocenter rearrangement in DNA hypomethylated cells.

Crosstalk between H2A variant-specific modifications impacts vital cell functions.

Selection of C-terminal motifs participated in evolution of distinct histone H2A variants. Hybrid types of variants combining motifs from distinct H2A classes are extremely rare. This suggests that the proximity between the motif cases interferes with their function. We studied this question in flowering plants that evolved sporadically a hybrid H2A variant combining the SQ motif of H2A.X that participates in the DNA damage response with the KSPK motif of H2A.W that stabilizes heterochromatin. Our inventory of PTMs of H2A.W variants showed that in vivo the cell cycle-dependent kinase CDKA phosphorylates the KSPK motif of H2A.W but only in absence of an SQ motif. Phosphomimicry of KSPK prevented DNA damage response by the SQ motif of the hybrid H2A.W/X variant. In a synthetic yeast expressing the hybrid H2A.W/X variant, phosphorylation of KSPK prevented binding of the BRCT-domain protein Mdb1 to phosphorylated SQ and impaired response to DNA damage. Our findings illustrate that PTMs mediate interference between the function of H2A variant specific C-terminal motifs. Such interference could explain the mutual exclusion of motifs that led to evolution of H2A variants.

Single-nucleus transcriptomic landscape of primate hippocampal aging.

The hippocampus plays a crucial role in learning and memory, and its progressive deterioration with age is functionally linked to a variety of human neurodegenerative diseases. Yet a systematic profiling of the aging effects on various hippocampal cell types in primates is still missing. Here, we reported a variety of new aging-associated phenotypic changes of the primate hippocampus. These include, in particular, increased DNA damage and heterochromatin erosion with time, alongside loss of proteostasis and elevated inflammation. To understand their cellular and molecular causes, we established the first single-nucleus transcriptomic atlas of primate hippocampal aging. Among the 12 identified cell types, neural transiently amplifying progenitor cell (TAPC) and microglia were most affected by aging. In-depth dissection of gene-expression dynamics revealed impaired TAPC division and compromised neuronal function along the neurogenesis trajectory; additionally elevated pro-inflammatory responses in the aged microglia and oligodendrocyte, as well as dysregulated coagulation pathways in the aged endothelial cells may contribute to a hostile microenvironment for neurogenesis. This rich resource for understanding primate hippocampal aging may provide potential diagnostic biomarkers and therapeutic interventions against age-related neurodegenerative diseases.

HP1 drives de novo 3D genome reorganization in early Drosophila embryos.

Fundamental features of 3D genome organization are established de novo in the early embryo, including clustering of pericentromeric regions, the folding of chromosome arms and the segregation of chromosomes into active (A-) and inactive (B-) compartments. However, the molecular mechanisms that drive de novo organization remain unknown1,2. Here, by combining chromosome conformation capture (Hi-C), chromatin immunoprecipitation with high-throughput sequencing (ChIP-seq), 3D DNA fluorescence in situ hybridization (3D DNA FISH) and polymer simulations, we show that heterochromatin protein 1a (HP1a) is essential for de novo 3D genome organization during Drosophila early development. The binding of HP1a at pericentromeric heterochromatin is required to establish clustering of pericentromeric regions. Moreover, HP1a binding within chromosome arms is responsible for overall chromosome folding and has an important role in the formation of B-compartment regions. However, depletion of HP1a does not affect the A-compartment, which suggests that a different molecular mechanism segregates active chromosome regions. Our work identifies HP1a as an epigenetic regulator that is involved in establishing the global structure of the genome in the early embryo.

The Dynamic Influence of Linker Histone Saturation within the Poly-Nucleosome Array.

Linker histones bind to nucleosomes and modify chromatin structure and dynamics as a means of epigenetic regulation. Biophysical studies have shown that chromatin fibers can adopt a plethora of conformations with varying levels of compaction. Linker histone condensation, and its specific binding disposition, has been associated with directly tuning this ensemble of states. However, the atomistic dynamics and quantification of this mechanism remains poorly understood. Here, we present molecular dynamics simulations of octa-nucleosome arrays, based on a cryo-EM structure of the 30-nm chromatin fiber, with and without the globular domains of the H1 linker histone to determine how they influence fiber structures and dynamics. Results show that when bound, linker histones inhibit DNA flexibility and stabilize repeating tetra-nucleosomal units, giving rise to increased chromatin compaction. Furthermore, upon the removal of H1, there is a significant destabilization of this compact structure as the fiber adopts less strained and untwisted states. Interestingly, linker DNA sampling in the octa-nucleosome is exaggerated compared to its mono-nucleosome counterparts, suggesting that chromatin architecture plays a significant role in DNA strain even in the absence of linker histones. Moreover, H1-bound states are shown to have increased stiffness within tetra-nucleosomes, but not between them. This increased stiffness leads to stronger long-range correlations within the fiber, which may result in the propagation of epigenetic signals over longer spatial ranges. These simulations highlight the effects of linker histone binding on the internal dynamics and global structure of poly-nucleosome arrays, while providing physical insight into a mechanism of chromatin compaction.

SINE-B1 Distribution and Chromosome Rearrangements in the South American Proechimys gr. goeldii (Echimyidae, Rodentia).

Proechimys species are remarkable for their extensive chromosome rearrangements, representing a good model to understand genome evolution. Herein, we cytogenetically analyzed 3 different cytotypes of Proechimys gr. goeldii to assess their evolutionary relationship. We also mapped the transposable element SINE-B1 on the chromosomes of P. gr. goeldii in order to investigate its distribution among individuals and evaluate its possible contribution to karyotype remodeling in this species. SINE-B1 showed a dispersed distribution along chromosome arms and was also detected at the pericentromeric regions of some chromosomes, including pair 1 and the sex chromosomes, which are involved in chromosome rearrangements. In addition, we describe a new cytotype for P. gr. goeldii, reinforcing the significant role of gross chromosomal rearrangements during the evolution of the genus. The results of FISH with SINE-B1 suggest that this issue should be more deeply investigated for a better understanding of its role in the mechanisms involved in the wide variety of Proechimys karyotypes.