Transcription of a centromere-enriched retroelement and local retention of its RNA are significant features of the CENP-A chromatin landscape

Centromeres depend on chromatin containing the conserved histone H3 variant CENP-A for function and inheritance, while the role of centromeric DNA repeats remains unclear. Retroelements are prevalent at centromeres across taxa and represent a potential mechanism for promoting transcription to aid in CENP-A incorporation or for generating RNA transcripts to maintain centromere integrity. Here, we probe into the transcription and RNA localization of the centromere-enriched retroelement G2/Jockey-3 (hereafter referred to as Jockey-3 ) in Drosophila melanogaster , currently the only in vivo model with assembled centromeres. We find that Jockey-3 is a major component of the centromeric transcriptome and produces RNAs that localize to centromeres in metaphase. Leveraging the polymorphism of Jockey-3 and a de novo centromere system, we show that these RNAs remain associated with their cognate DNA sequences in cis , suggesting they are unlikely to perform a sequence-specific function at all centromeres. We show that Jockey-3 transcription is positively correlated with the presence of CENP-A, and that recent Jockey-3 transposition events have occurred preferentially at CENP-A-containing chromatin. We propose that Jockey-3 contributes to the epigenetic maintenance of centromeres by promoting chromatin transcription, while inserting preferentially within these regions, selfishly ensuring its continued expression and transmission. Given the conservation of retroelements as centromere components through evolution, our findings have broad implications in understanding this association in other species.

Unveiling the genetic blueprint of a desert scorpion: A chromosome-level genome of Hadrurus arizonensis provides the first reference for Parvorder Iurida

Over 400 million years old, scorpions represent an ancient group of arachnids and one of the first animals to adapt to life on land. Presently, the lack of available genomes within scorpions hinders research on their evolution. This study leverages ultra-long nanopore sequencing and Pore-C to generate the first chromosome level assembly and annotation for the desert hairy scorpion, Hadrurus arizonensis. The assembled genome is 2.23 Gb in size with an N50 of 280 Mb. Pore-C scaffolding re-oriented 99.6% of bases into nine chromosomes and BUSCO identified 998 (98.6%) complete arthropod single copy orthologs. Repetitive elements represent 54.69% of the assembled bases, including 872,874 (29.39%) LINE elements. A total of 18,996 protein-coding genes and 75,256 transcripts were predicted, and extracted protein sequences yielded a BUSCO score of 97.2%. This is the first genome assembled and annotated within the family Hadruridae, representing a crucial resource for closing gaps in genomic knowledge of scorpions, resolving arachnid phylogeny, and advancing studies in comparative and functional genomics.

Profiling genome-wide methylation in two maples: Fine-scale approaches to detection with nanopore technology

DNA methylation is critical to the regulation of transposable elements and gene expression and can play an important role in the adaptation of stress response mechanisms in plants. Traditional methods of methylation quantification rely on bisulfite conversion that can compromise accuracy. Recent advances in long-read sequencing technologies allow for methylation detection in real time. The associated algorithms that interpret these modifications have evolved from strictly statistical approaches to Hidden Markov Models and, recently, deep learning approaches. Much of the existing software focuses on methylation in the CG context, but methylation in other contexts is important to quantify, as it is extensively leveraged in plants. Here, we present methylation profiles for two maple species across the full range of 5mC sequence contexts using Oxford Nanopore Technologies (ONT) long-reads. Hybrid and reference-guided assemblies were generated for two new Acer accessions: Acer negundo (box elder; 65x ONT and 111X Illumina) and Acer saccharum (sugar maple; 93x ONT and 148X Illumina). The ONT reads generated for these assemblies were re-basecalled, and methylation detection was conducted in a custom pipeline with the published Acer references (PacBio assemblies) and hybrid assemblies reported herein to generate four epigenomes. Examination of the transposable element landscape revealed the dominance of LTR Copia elements and patterns of methylation associated with different classes of TEs. Methylation distributions were examined at high resolution across gene and repeat density and described within the broader angiosperm context, and more narrowly in the context of gene family dynamics and candidate nutrient stress genes.

Salpa genome and developmental transcriptome analyses reveal molecular flexibility enabling reproductive success in a rapidly changing environment

Ocean warming favors pelagic tunicates, such as salps, that exhibit increasingly frequent and rapid population blooms, impacting trophic dynamics and composition and human marine-dependent activities. Salp blooms are a result of their successful reproductive life history, alternating seasonally between asexual and sexual protogynous (i.e. sequential) hermaphroditic stages. While predicting future salp bloom frequency and intensity relies on an understanding of the transitions during the sexual stage from female through parturition and subsequent sex change to male, these transitions have not been explored at the molecular level. Here we report the development of the first complete genome of S. thompsoni and the North Atlantic sister species S. aspera. Genome and comparative analyses reveal an abundance of repeats and G-quadruplex (G4) motifs, a highly stable secondary structure, distributed throughout both salp genomes, a feature shared with other tunicates that perform alternating sexual-asexual reproductive strategies. Transcriptional analyses across sexual reproductive stages for S. thompsoni revealed genes associated with male sex differentiation and spermatogenesis are expressed as early as birth and before parturition, inconsistent with previous descriptions of sequential sexual differentiation in salps. Our findings suggest salp are poised for reproductive success at birth, increasing the potential for bloom formation as ocean temperatures rise.

Kangaroo endogenous retrovirus (KERV) forms megasatellite DNA with a simple repetition pattern in which the provirus structure is retained

The kangaroo endogenous retrovirus (KERV) was previously reported to have undergone a rapid copy number increase in the red-necked wallaby; however, the mode of amplification was left to be clarified. The present study revealed that the long terminal repeat (LTR) (0.6 kb) and internal region (2.0 kb) of a provirus are repeated alternately, forming megasatellite DNA which we named kervRep. This repetition pattern was the same as that observed for walbRep, megasatellite DNA originating from another endogenous retrovirus. Their formation process can be explained using a simple model: pairing slippage followed by homologous recombination. This model features that the initial step is triggered by the presence of two identical sequences within a short distance; the possession of LTRs by endogenous retroviruses fulfills this condition. The discovery of two cases suggests that formation of this type of satellite DNA is one of non-negligible effects of endogenous retroviruses on their host genomes.

Assembly of 43 human Y chromosomes reveals extensive complexity and variation

The prevalence of highly repetitive sequences within the human Y chromosome has prevented its complete assembly to date1 and led to its systematic omission from genomic analyses. Here we present de novo assemblies of 43 Y chromosomes spanning 182,900 years of human evolution and report considerable diversity in size and structure. Half of the male-specific euchromatic region is subject to large inversions with a greater than twofold higher recurrence rate compared with all other chromosomes2. Ampliconic sequences associated with these inversions show differing mutation rates that are sequence context dependent, and some ampliconic genes exhibit evidence for concerted evolution with the acquisition and purging of lineage-specific pseudogenes. The largest heterochromatic region in the human genome, Yq12, is composed of alternating repeat arrays that show extensive variation in the number, size and distribution, but retain a 1:1 copy-number ratio. Finally, our data suggest that the boundary between the recombining pseudoautosomal region 1 and the non-recombining portions of the X and Y chromosomes lies 500 kb away from the currently established1 boundary. The availability of fully sequence-resolved Y chromosomes from multiple individuals provides a unique opportunity for identifying new associations of traits with specific Y-chromosomal variants and garnering insights into the evolution and function of complex regions of the human genome.

The complete sequence of a human Y chromosome

The human Y chromosome has been notoriously difficult to sequence and assemble because of its complex repeat structure that includes long palindromes, tandem repeats and segmental duplications1-3. As a result, more than half of the Y chromosome is missing from the GRCh38 reference sequence and it remains the last human chromosome to be finished4,5. Here, the Telomere-to-Telomere (T2T) consortium presents the complete 62,460,029-base-pair sequence of a human Y chromosome from the HG002 genome (T2T-Y) that corrects multiple errors in GRCh38-Y and adds over 30 million base pairs of sequence to the reference, showing the complete ampliconic structures of gene families TSPY, DAZ and RBMY; 41 additional protein-coding genes, mostly from the TSPY family; and an alternating pattern of human satellite 1 and 3 blocks in the heterochromatic Yq12 region. We have combined T2T-Y with a previous assembly of the CHM13 genome4 and mapped available population variation, clinical variants and functional genomics data to produce a complete and comprehensive reference sequence for all 24 human chromosomes.

Deep statistical modelling of nanopore sequencing translocation times reveals latent non-B DNA structures

MOTIVATION

Non-canonical (or non-B) DNA are genomic regions whose three-dimensional conformation deviates from the canonical double helix. Non-B DNA play an important role in basic cellular processes and are associated with genomic instability, gene regulation, and oncogenesis. Experimental methods are low-throughput and can detect only a limited set of non-B DNA structures, while computational methods rely on non-B DNA base motifs, which are necessary but not sufficient indicators of non-B structures. Oxford Nanopore sequencing is an efficient and low-cost platform, but it is currently unknown whether nanopore reads can be used for identifying non-B structures.

RESULTS

We build the first computational pipeline to predict non-B DNA structures from nanopore sequencing. We formalize non-B detection as a novelty detection problem and develop the GoFAE-DND, an autoencoder that uses goodness-of-fit (GoF) tests as a regularizer. A discriminative loss encourages non-B DNA to be poorly reconstructed and optimizing Gaussian GoF tests allows for the computation of P-values that indicate non-B structures. Based on whole genome nanopore sequencing of NA12878, we show that there exist significant differences between the timing of DNA translocation for non-B DNA bases compared with B-DNA. We demonstrate the efficacy of our approach through comparisons with novelty detection methods using experimental data and data synthesized from a new translocation time simulator. Experimental validations suggest that reliable detection of non-B DNA from nanopore sequencing is achievable.

AVAILABILITY AND IMPLEMENTATION

Source code is available at https://github.com/bayesomicslab/ONT-nonb-GoFAE-DND.

Derivation of splice junction-specific antibodies using a unique hapten targeting strategy and directed evolution

Alternative splicing of RNA occurs frequently in eukaryotic cells and can result in multiple protein isoforms that are nearly identical in amino acid sequence, but have unique biological roles. Moreover, the relative abundance of these unique isoforms can be correlative with diseased states and potentially used as biomarkers or therapeutic targets. However, due to high sequence similarities among isoforms, current proteomic methods are incapable of differentiating native protein isoforms derived from most alternative splicing events. Herein, a strategy employing a nonsynonymous, non-native amino acid (nnAA) pseudo-hapten (i.e. an amino acid or amino acid derivative that is different from the native amino acid at a particular position) as a targeting epitope in splice junction-spanning peptides was successful in directed antibody derivation. After isolating nnAA-specific antibodies, directed evolution reduced the antibody's binding dependence on the nnAA pseudo-hapten and improved binding to the native splice junction epitope. The resulting antibodies demonstrated codependent binding affinity to each exon of the splice junction and thus are splice junction- and isoform-specific. Furthermore, epitope scanning demonstrated that positioning of the nnAA pseudo-hapten within a peptide antigen can be exploited to predetermine the isolated antibody's specificity at, or near, amino acid resolution. Thus, this nnAA targeting strategy has the potential to robustly derive splice junction- and site-specific antibodies that can be used in a wide variety of research endeavors to unambiguously differentiate native protein isoforms.

Epigenetic patterns in a complete human genome

The completion of a telomere-to-telomere human reference genome, T2T-CHM13, has resolved complex regions of the genome, including repetitive and homologous regions. Here, we present a high-resolution epigenetic study of previously unresolved sequences, representing entire acrocentric chromosome short arms, gene family expansions, and a diverse collection of repeat classes. This resource precisely maps CpG methylation (32.28 million CpGs), DNA accessibility, and short-read datasets (166,058 previously unresolved chromatin immunoprecipitation sequencing peaks) to provide evidence of activity across previously unidentified or corrected genes and reveals clinically relevant paralog-specific regulation. Probing CpG methylation across human centromeres from six diverse individuals generated an estimate of variability in kinetochore localization. This analysis provides a framework with which to investigate the most elusive regions of the human genome, granting insights into epigenetic regulation.

Complete genomic and epigenetic maps of human centromeres

Existing human genome assemblies have almost entirely excluded repetitive sequences within and near centromeres, limiting our understanding of their organization, evolution, and functions, which include facilitating proper chromosome segregation. Now, a complete, telomere-to-telomere human genome assembly (T2T-CHM13) has enabled us to comprehensively characterize pericentromeric and centromeric repeats, which constitute 6.2% of the genome (189.9 megabases). Detailed maps of these regions revealed multimegabase structural rearrangements, including in active centromeric repeat arrays. Analysis of centromere-associated sequences uncovered a strong relationship between the position of the centromere and the evolution of the surrounding DNA through layered repeat expansions. Furthermore, comparisons of chromosome X centromeres across a diverse panel of individuals illuminated high degrees of structural, epigenetic, and sequence variation in these complex and rapidly evolving regions.

Comparative analyses of gibbon centromeres reveal dynamic genus specific shifts in repeat composition

Centromeres are functionally conserved chromosomal loci essential for proper chromosome segregation during cell division, yet they show high sequence diversity across species. Despite their variation, a near universal feature of centromeres is the presence of repetitive sequences, such as DNA satellites and transposable elements (TEs). Because of their rapidly evolving karyotypes, gibbons represent a compelling model to investigate divergence of functional centromere sequences across short evolutionary timescales. In this study, we use ChIP-seq, RNA-seq, and fluorescence in situ hybridization to comprehensively investigate the centromeric repeat content of the four extant gibbon genera (Hoolock, Hylobates, Nomascus, and Siamang). In all gibbon genera, we find that CENP-A nucleosomes and the DNA-proteins that interface with the inner kinetochore preferentially bind retroelements of broad classes rather than satellite DNA. A previously identified gibbon-specific composite retrotransposon, LAVA, known to be expanded within the centromere regions of one gibbon genus (Hoolock), displays centromere- and species-specific sequence differences, potentially as a result of its co-option to a centromeric function. When dissecting centromere satellite composition, we discovered the presence of the retroelement-derived macrosatellite SST1 in multiple centromeres of Hoolock, whereas alpha-satellites represent the predominate satellite in the other genera, further suggesting an independent evolutionary trajectory for Hoolock centromeres. Finally, using de novo assembly of centromere sequences, we determined that transcripts originating from gibbon centromeres recapitulate the species-specific TE composition. Combined, our data reveal dynamic shifts in the repeat content that define gibbon centromeres and coincide with the extensive karyotypic diversity within this lineage.

Seq'ing identity and function in a repeat-derived noncoding RNA world

Innovations in high-throughout sequencing approaches are being marshaled to both reveal the composition of the abundant and heterogeneous noncoding RNAs that populate cell nuclei and lend insight to the mechanisms by which noncoding RNAs influence chromosome biology and gene expression. This review focuses on some of the recent technological developments that have enabled the isolation of nascent transcripts and chromatin-associated and DNA-interacting RNAs. Coupled with emerging genome assembly and analytical approaches, the field is poised to achieve a comprehensive catalog of nuclear noncoding RNAs, including those derived from repetitive regions within eukaryotic genomes. Herein, particular attention is paid to the challenges and advances in the sequence analyses of repeat and transposable element-derived noncoding RNAs and in ascribing specific function(s) to such RNAs.

Evolution of Marsupial Genomes

Marsupial genomes, which are packaged into large chromosomes, provide a powerful resource for studying the mechanisms of genome evolution. The extensive and valuable body of work on marsupial cytogenetics, combined more recently with genome sequence data, has enabled prediction of the 2n = 14 karyotype ancestral to all marsupial families. The application of both chromosome biology and genome sequencing, or chromosomics, has been a necessary approach for various aspects of mammalian genome evolution, such as understanding sex chromosome evolution and the origin and evolution of transmissible tumors in Tasmanian devils. The next phase of marsupial genome evolution research will employ chromosomics approaches to begin addressing fundamental questions in marsupial genome evolution and chromosome evolution more generally. The answers to these complex questions will impact our understanding across a broad range of fields, including the genetics of speciation, genome adaptation to environmental stressors, and species management.

A centromere satellite concomitant with extensive karyotypic diversity across the Peromyscus genus defies predictions of molecular drive

A common feature of eukaryotic centromeres is the presence of large tracts of tandemly arranged repeats, known as satellite DNA. However, these centromeric repeats appear to experience rapid evolution under forces such as molecular drive and centromere drive, seemingly without consequence to the integrity of the centromere. Moreover, blocks of heterochromatin within the karyotype, including the centromere, are hotspots for chromosome rearrangements that may drive speciation events by contributing to reproductive isolation. However, the relationship between the evolution of heterochromatic sequences and the karyotypic dynamics of these regions remains largely unknown. Here, we show that a single conserved satellite DNA sequence in the order Rodentia of the genus Peromyscus localizes to recurrent sites of chromosome rearrangements and heterochromatic amplifications. Peromyscine species display several unique features of chromosome evolution compared to other Rodentia, including stable maintenance of a strict chromosome number of 48 among all known species in the absence of any detectable interchromosomal rearrangements. Rather, the diverse karyotypes of Peromyscine species are due to intrachromosomal variation in blocks of repeated DNA content. Despite wide variation in the copy number and location of repeat blocks among different species, we find that a single satellite monomer maintains a conserved sequence and homogenized tandem repeat structure, defying predictions of molecular drive. The conservation of this satellite monomer results in common, abundant, and large blocks of chromatin that are homologous among chromosomes within one species and among diverged species. Thus, such a conserved repeat may have facilitated the retention of polymorphic chromosome variants within individuals and intrachromosomal rearrangements between species-both factors that have previously been hypothesized to contribute towards the extremely wide range of ecological adaptations that this genus exhibits.

FURIN variant associations with postexercise hypotension are intensity and race dependent

FURIN is a proprotein convertase subtilisin/kexin enzyme important in pro-renin receptor processing, and FURIN (furin, paired basic amino acid cleaving enzyme) variants are involved in multiple aspects of blood pressure (BP) regulation. Therefore, we examined associations among FURIN variants and the immediate blood pressure (BP) response to bouts of aerobic exercise, termed postexercise hypotension (PEH). Obese (30.9 ± 3.6 kg  m-2 ) Black- (n = 14) and White- (n = 9) adults 42.0 ± 9.8 year with hypertension (139.8 ± 10.4/84.6 ± 6.2 mmHg) performed three random experiments: bouts of vigorous (VIGOROUS) and moderate (MODERATE) intensity cycling and control. Subjects were then attached to an ambulatory BP monitor for 19 h. We performed deep-targeted exon sequencing with the Illumina TruSeq Custom Amplicon kit. FURIN genotypes were coded as the number of minor alleles (#MA) and selected for additional statistical analysis based upon Bonferonni or Benjamini-Yekutieli multiple testing corrected P-values under time-adjusted linear models for 19 hourly BP measurements. After VIGOROUS over 19 h, as FURIN #MA increased in rs12917264 (P = 2.4E-04) and rs75493298 (P = 6.4E-04), systolic BP (SBP) decreased 30.4-33.7 mmHg; and in rs12917264 (P = 1.6E-03) and rs75493298 (P = 9.7E-05), diastolic BP (DBP) decreased 17.6-20.3 mmHg among Blacks only. In addition, after MODERATE over 19 h in FURIN rs74037507 (P = 8.0E-04), as #MA increased, SBP increased 20.8 mmHg among Blacks only. Whereas, after MODERATE over the awake hours in FURIN rs1573644 (P = 6.2E-04), as #MA increased, DBP decreased 12.5 mmHg among Whites only. FURIN appears to exhibit intensity and race-dependent associations with PEH that merit further exploration among a larger, ethnically diverse sample of adults with hypertension.

Continuous-Trait Probabilistic Model for Comparing Multi-species Functional Genomic Data

A large amount of multi-species functional genomic data from high-throughput assays are becoming available to help understand the molecular mechanisms for phenotypic diversity across species. However, continuous-trait probabilistic models, which are key to such comparative analysis, remain under-explored. Here we develop a new model, called phylogenetic hidden Markov Gaussian processes (Phylo-HMGP), to simultaneously infer heterogeneous evolutionary states of functional genomic features in a genome-wide manner. Both simulation studies and real data application demonstrate the effectiveness of Phylo-HMGP. Importantly, we applied Phylo-HMGP to analyze a new cross-species DNA replication timing (RT) dataset from the same cell type in five primate species (human, chimpanzee, orangutan, gibbon, and green monkey). We demonstrate that our Phylo-HMGP model enables discovery of genomic regions with distinct evolutionary patterns of RT. Our method provides a generic framework for comparative analysis of multi-species continuous functional genomic signals to help reveal regions with conserved or lineage-specific regulatory roles.

Adaptation and conservation insights from the koala genome

The koala, the only extant species of the marsupial family Phascolarctidae, is classified as 'vulnerable' due to habitat loss and widespread disease. We sequenced the koala genome, producing a complete and contiguous marsupial reference genome, including centromeres. We reveal that the koala's ability to detoxify eucalypt foliage may be due to expansions within a cytochrome P450 gene family, and its ability to smell, taste and moderate ingestion of plant secondary metabolites may be due to expansions in the vomeronasal and taste receptors. We characterized novel lactation proteins that protect young in the pouch and annotated immune genes important for response to chlamydial disease. Historical demography showed a substantial population crash coincident with the decline of Australian megafauna, while contemporary populations had biogeographic boundaries and increased inbreeding in populations affected by historic translocations. We identified genetically diverse populations that require habitat corridors and instituting of translocation programs to aid the koala's survival in the wild.

Tissue-Specific Transcriptome for Poeciliopsis prolifica Reveals Evidence for Genetic Adaptation Related to the Evolution of a Placental Fish

The evolution of the placenta is an excellent model to examine the evolutionary processes underlying adaptive complexity due to the recent, independent derivation of placentation in divergent animal lineages. In fishes, the family Poeciliidae offers the opportunity to study placental evolution with respect to variation in degree of post-fertilization maternal provisioning among closely related sister species. In this study, we present a detailed examination of a new reference transcriptome sequence for the live-bearing, matrotrophic fish, Poeciliopsis prolifica, from multiple-tissue RNA-seq data. We describe the genetic components active in liver, brain, late-stage embryo, and the maternal placental/ovarian complex, as well as associated patterns of positive selection in a suite of orthologous genes found in fishes. Results indicate the expression of many signaling transcripts, "non-coding" sequences and repetitive elements in the maternal placental/ovarian complex. Moreover, patterns of positive selection in protein sequence evolution were found associated with live-bearing fishes, generally, and the placental P. prolifica, specifically, that appear independent of the general live-bearer lifestyle. Much of the observed patterns of gene expression and positive selection are congruent with the evolution of placentation in fish functionally converging with mammalian placental evolution and with the patterns of rapid evolution facilitated by the teleost-specific whole genome duplication event.

Epigenetic maintenance of topological domains in the highly rearranged gibbon genome

The relationship between evolutionary genome remodeling and the three-dimensional structure of the genome remain largely unexplored. Here, we use the heavily rearranged gibbon genome to examine how evolutionary chromosomal rearrangements impact genome-wide chromatin interactions, topologically associating domains (TADs), and their epigenetic landscape. We use high-resolution maps of gibbon–human breaks of synteny (BOS), apply Hi-C in gibbon, measure an array of epigenetic features, and perform cross-species comparisons. We find that gibbon rearrangements occur at TAD boundaries, independent of the parameters used to identify TADs. This overlap is supported by a remarkable genetic and epigenetic similarity between BOS and TAD boundaries, namely presence of CpG islands and SINE elements, and enrichment in CTCF and H3K4me3 binding. Cross-species comparisons reveal that regions orthologous to BOS also correspond with boundaries of large (400–600 kb) TADs in human and other mammalian species. The colocalization of rearrangement breakpoints and TAD boundaries may be due to higher chromatin fragility at these locations and/or increased selective pressure against rearrangements that disrupt TAD integrity. We also examine the small portion of BOS that did not overlap with TAD boundaries and gave rise to novel TADs in the gibbon genome. We postulate that these new TADs generally lack deleterious consequences. Last, we show that limited epigenetic homogenization occurs across breakpoints, irrespective of their time of occurrence in the gibbon lineage. Overall, our findings demonstrate remarkable conservation of chromatin interactions and epigenetic landscape in gibbons, in spite of extensive genomic shuffling.

An updated genetic map of Peromyscus with chromosomal assignment of linkage groups

Species across the rodent genus Peromyscus have become prominent models for studying diverse mechanistic and evolutionary processes, including chromosome evolution, infectious disease transmission and human health, ecological adaptation, coat color variation, and parental care. Supporting such diverse research programs has been the development of genetic and genomic resources for species within this genus, including genome data, interspecific chromosome homologies, and a recently developed genetic map. Based on interspecific hybrids between the deer mouse (Peromyscus maniculatus bairdii) and the old-field, or beach, mouse (Peromyscus polionotus) and backcross progeny to Peromyscus maniculatus, a linkage map was developed based on 190 genes and 141 microsatellite loci. However, resolution of several linkage groups with respect to chromosome assignment was lacking and four chromosomes (8, 16, 20, and 21) were not clearly delineated with linkage data alone. The recent development of a high-density map for Peromyscus proved ineffective in resolving chromosome linkage for these four chromosomes. Herein we present an updated linkage map for Peromyscus maniculatus, including linkage group-chromosome assignments, using fluorescence in situ hybridization mapping of BACs and whole chromosome paints. We resolve the previously conflicting chromosome assignment of linkage groups to Chromosomes 8, 16, 20, and 21, and confirm the assignment of linkage groups to Chromosomes 18 and 22. This updated linkage map with validated chromosome assignment provides a solid foundation for chromosome nomenclature for this species.

Sex chromosome repeats tip the balance towards speciation

Because sex chromosomes, by definition, carry genes that determine sex, mutations that alter their structural and functional stability can have immediate consequences for the individual by reducing fertility, but also for a species by altering the sex ratio. Moreover, the sex-specific segregation patterns of heteromorphic sex chromosomes make them havens for selfish genetic elements that not only create suboptimal sex ratios but can also foster sexual antagonism. Compensatory mutations to mitigate antagonism or return sex ratios to a Fisherian optimum can create hybrid incompatibility and establish reproductive barriers leading to species divergence. The destabilizing influence of these selfish elements is often manifest within populations as copy number variants (CNVs) in satellite repeats and transposable elements (TE) or as CNVs involving sex-determining genes, or genes essential to fertility and sex chromosome dosage compensation. This review catalogs several examples of well-studied sex chromosome CNVs in Drosophilids and mammals that underlie instances of meiotic drive, hybrid incompatibility and disruptions to sex differentiation and sex chromosome dosage compensation. While it is difficult to pinpoint a direct cause/effect relationship between these sex chromosome CNVs and speciation, it is easy to see how their effects in creating imbalances between the sexes, and the compensatory mutations to restore balance, can lead to lineage splitting and species formation.

Replication timing kept in LINE

O'Neill and O'Neill discuss Platt et al.’s findings that LINE1 elements are key to control of replication timing by ASAR long noncoding RNAs.

Accurate and synchronous replication timing between chromosome homologues is essential for maintaining chromosome stability, yet how this is achieved has remained a mystery. In this issue, Platt et al. (2018. J. Cell Biol.https://doi.org/10.1083/jcb.201707082) identify antisense LINE (L1) transcripts within long noncoding RNAs as the critical factor in maintaining synchronous chromosome-wide replication timing.

Transposable elements: genome innovation, chromosome diversity, and centromere conflict

Although it was nearly 70 years ago when transposable elements (TEs) were first discovered “jumping” from one genomic location to another, TEs are now recognized as contributors to genomic innovations as well as genome instability across a wide variety of species. In this review, we illustrate the ways in which active TEs, specifically retroelements, can create novel chromosome rearrangements and impact gene expression, leading to disease in some cases and species-specific diversity in others. We explore the ways in which eukaryotic genomes have evolved defense mechanisms to temper TE activity and the ways in which TEs continue to influence genome structure despite being rendered transpositionally inactive. Finally, we focus on the role of TEs in the establishment, maintenance, and stabilization of critical, yet rapidly evolving, chromosome features: eukaryotic centromeres. Across centromeres, specific types of TEs participate in genomic conflict, a balancing act wherein they are actively inserting into centromeric domains yet are harnessed for the recruitment of centromeric histones and potentially new centromere formation.

Genome of the Tasmanian tiger provides insights into the evolution and demography of an extinct marsupial carnivore

The Tasmanian tiger or thylacine (Thylacinus cynocephalus) was the largest carnivorous Australian marsupial to survive into the modern era. Despite last sharing a common ancestor with the eutherian canids ~160 million years ago, their phenotypic resemblance is considered the most striking example of convergent evolution in mammals. The last known thylacine died in captivity in 1936 and many aspects of the evolutionary history of this unique marsupial apex predator remain unknown. Here we have sequenced the genome of a preserved thylacine pouch young specimen to clarify the phylogenetic position of the thylacine within the carnivorous marsupials, reconstruct its historical demography and examine the genetic basis of its convergence with canids. Retroposon insertion patterns placed the thylacine as the basal lineage in Dasyuromorphia and suggest incomplete lineage sorting in early dasyuromorphs. Demographic analysis indicated a long-term decline in genetic diversity starting well before the arrival of humans in Australia. In spite of their extraordinary phenotypic convergence, comparative genomic analyses demonstrated that amino acid homoplasies between the thylacine and canids are largely consistent with neutral evolution. Furthermore, the genes and pathways targeted by positive selection differ markedly between these species. Together, these findings support models of adaptive convergence driven primarily by cis-regulatory evolution.

Deep-targeted sequencing of endothelial nitric oxide synthase gene exons uncovers exercise intensity and ethnicity-dependent associations with post-exercise hypotension

In previous studies, we found an endothelial nitric oxide synthase gene (NOS3) variant rs2070744 associated with the ambulatory blood pressure (BP) response following bouts of moderate and vigorous intensity acute exercise, termed post-exercise hypotension (PEH). In a validation cohort, we sequenced NOS3 exons for associations with PEH Obese (30.9 ± 3.6 kg.m-2) African American (n = 14) [AF] and Caucasian (n = 9) adults 42.0 ± 9.8 years with hypertension (139.8 ± 10.4/84.6 ± 6.2 mmHg) performed three random experiments: bouts of vigorous and moderate intensity cycling and control. Subjects were attached to an ambulatory BP monitor for 19 h. We performed deep-targeted exon sequencing with the Illumina TruSeq Custom Amplicon kit. Variant genotypes were coded as number of minor alleles (#MA) and selected for additional statistical analysis based upon Bonferonni or Benjamini-Yekutieli multiple testing-corrected P-values under time-adjusted linear models for 19 hourly BP measurements for each subject. After vigorous intensity over 19 h, among NOS3 variants passing multiple testing thresholds, as the #MA increased in rs891512 (P = 6.4E-04), rs867225 (P = 6.5E-04), rs743507 (P = 2.6E-06), and rs41483644 (P = 2.4E-04), systolic (SBP) decreased from 17.5 to 33.7 mmHg; and in rs891512 (P = 9.7E-05), rs867225 (P = 2.6E-05), rs41483644 (P = 1.6E-03), rs3730009 (P = 2.6E-04), and rs77325852 (P = 5.6E-04), diastolic BP decreased from 11.1 mmHg to 20.3 mmHg among AF only. In contrast, after moderate intensity over 19 h in NOS3 rs3918164, as the #MA increased, SBP increased by 16.6 mmHg (P = 2.4E-04) among AF only. NOS3 variants exhibited associations with PEH after vigorous, but not moderate intensity exercise among AF only. NOS3 should be studied further for its effects on PEH in a large, ethnically diverse sample of adults with hypertension to confirm our findings.

RNA-dependent stabilization of SUV39H1 at constitutive heterochromatin

Heterochromatin formed by the SUV39 histone methyltransferases represses transcription from repetitive DNA sequences and ensures genomic stability. How SUV39 enzymes localize to their target genomic loci remains unclear. Here, we demonstrate that chromatin-associated RNA contributes to the stable association of SUV39H1 with constitutive heterochromatin in human cells. We find that RNA associated with mitotic chromosomes is concentrated at pericentric heterochromatin, and is encoded, in part, by repetitive α-satellite sequences, which are retained in cis at their transcription sites. Purified SUV39H1 directly binds nucleic acids through its chromodomain; and in cells, SUV39H1 associates with α-satellite RNA transcripts. Furthermore, nucleic acid binding mutants destabilize the association of SUV39H1 with chromatin in mitotic and interphase cells – effects that can be recapitulated by RNase treatment or RNA polymerase inhibition – and cause defects in heterochromatin function. Collectively, our findings uncover a previously unrealized function for chromatin-associated RNA in regulating constitutive heterochromatin in human cells.

DOI:http://dx.doi.org/10.7554/eLife.25299.001

Each cell in a human body contains the same DNA sequence, which serves as a set of instructions for how the body should develop and operate. However, only certain sections of DNA are “active” at any particular time and in any given type of cell. When a section of DNA is active, cells make many copies of it using a molecule called RNA. When a section of DNA in inactive, very little RNA is made. Some sections of DNA must always be kept inactive to avoid damaging the cell.

DNA is packaged around proteins called histones, and enzymes that modify histones control which sections of DNA are switched on or off. One such modifying enzyme, called SUV39H1, is important for inactivating sections of DNA that could cause harm to the cell if they are active. Previous studies showed that the loss of SUV39H1 and related proteins cause abnormalities and cancer in mice. However, it is not clear how this enzyme identifies and inactivates the DNA it needs to target.

Johnson, Yewdell et al. studied SUV39H1 in human cells. The experiments show that RNA binds to the SUV39H1 enzyme and controls how it interacts with DNA. Specifically, Johnson, Yewdell et al. found that sections of DNA that are inactive can still make a small amount of RNA, and that this RNA tethers SUV39H1 to the DNA to keep the DNA switched off. Mutant forms of SUV39H1 that are unable to interact with RNA fall off the DNA, which allows DNA sequences that are normally switched off to become active.

The findings of Johnson, Yewdell et al. reveal a new role for RNAs in regulating whether DNA is switched on or off. The next step is to determine whether other enzymes that can also modify histones use the same mechanism to activate or inactivate DNA. Differences in how the activity of DNA is regulated between individuals plays a crucial role in generating the diversity we see in nature. Therefore, this work helps us to understand our basic biology and may provide new opportunities for treating disease.

DOI:http://dx.doi.org/10.7554/eLife.25299.002

Rapid Evolutionary Rates and Unique Genomic Signatures Discovered in the First Reference Genome for the Southern Ocean Salp, Salpa thompsoni (Urochordata, Thaliacea)

A preliminary genome sequence has been assembled for the Southern Ocean salp, Salpa thompsoni (Urochordata, Thaliacea). Despite the ecological importance of this species in Antarctic pelagic food webs and its potential role as an indicator of changing Southern Ocean ecosystems in response to climate change, no genomic resources are available for S. thompsoni or any closely related urochordate species. Using a multiple-platform, multiple-individual approach, we have produced a 318,767,936-bp genome sequence, covering >50% of the estimated 602 Mb (±173 Mb) genome size for S. thompsoni. Using a nonredundant set of predicted proteins, >50% (16,823) of sequences showed significant homology to known proteins and ∼38% (12,151) of the total protein predictions were associated with Gene Ontology functional information. We have generated 109,958 SNP variant and 9,782 indel predictions for this species, serving as a resource for future phylogenomic and population genetic studies. Comparing the salp genome to available assemblies for four other urochordates, Botryllus schlosseri, Ciona intestinalis, Ciona savignyi and Oikopleura dioica, we found that S. thompsoni shares the previously estimated rapid rates of evolution for these species. High mutation rates are thus independent of genome size, suggesting that rates of evolution >1.5 times that observed for vertebrates are a broad taxonomic characteristic of urochordates. Tests for positive selection implemented in PAML revealed a small number of genes with sites undergoing rapid evolution, including genes involved in ribosome biogenesis and metabolic and immune process that may be reflective of both adaptation to polar, planktonic environments as well as the complex life history of the salps. Finally, we performed an initial survey of small RNAs, revealing the presence of known, conserved miRNAs, as well as novel miRNA genes; unique piRNAs; and mature miRNA signatures for varying developmental stages. Collectively, these resources provide a genomic foundation supporting S. thompsoni as a model species for further examination of the exceptional rates and patterns of genomic evolution shown by urochordates. Additionally, genomic data will allow for the development of molecular indicators of key life history events and processes and afford new understandings and predictions of impacts of climate change on this key species of Antarctic pelagic ecosystems.

The repetitive DNA element BncDNA, enriched in the B chromosome of the cichlid fish Astatotilapia latifasciata, transcribes a potentially noncoding RNA

Supernumerary chromosomes have been studied in many species of eukaryotes, including the cichlid fish, Astatotilapia latifasciata. However, there are many unanswered questions about the maintenance, inheritance, and functional aspects of supernumerary chromosomes. The cichlid family has been highlighted as a model for evolutionary studies, including those that focus on mechanisms of chromosome evolution. Individuals of A. latifasciata are known to carry up to two B heterochromatic isochromosomes that are enriched in repetitive DNA and contain few intact gene sequences. We isolated and characterized a transcriptionally active repeated DNA, called B chromosome noncoding DNA (BncDNA), highly represented across all B chromosomes of A. latifasciata. BncDNA transcripts are differentially processed among six different tissues, including the production of smaller transcripts, indicating transcriptional variation may be linked to B chromosome presence and sexual phenotype. The transcript lengths and lack of similarity with known protein/gene sequences indicate BncRNA might represent a novel long noncoding RNA family (lncRNA). The potential for interaction between BncRNA and known miRNAs were computationally predicted, resulting in the identification of possible binding of this sequence in upregulated miRNAs related to the presence of B chromosomes. In conclusion, Bnc is a transcriptionally active repetitive DNA enriched in B chromosomes with potential action over B chromosome maintenance in somatic cells and meiotic drive in gametic cells.

Establishment of Centromeric Chromatin by the CENP-A Assembly Factor CAL1 Requires FACT-Mediated Transcription

Centromeres are essential chromosomal structures that mediate accurate chromosome segregation during cell division. Centromeres are specified epigenetically by the heritable incorporation of the centromeric histone H3 variant CENP-A. While many of the primary factors that mediate centromeric deposition of CENP-A are known, the chromatin and DNA requirements of this process have remained elusive. Here, we uncover a role for transcription in Drosophila CENP-A deposition. Using an inducible ectopic centromere system that uncouples CENP-A deposition from endogenous centromere function and cell-cycle progression, we demonstrate that CENP-A assembly by its loading factor, CAL1, requires RNAPII-mediated transcription of the underlying DNA. This transcription depends on the CAL1 binding partner FACT, but not on CENP-A incorporation. Our work establishes RNAPII passage as a key step in chaperone-mediated CENP-A chromatin establishment and propagation.

Identification of a recently active mammalian SINE derived from ribosomal RNA

Complex eukaryotic genomes are riddled with repeated sequences whose derivation does not coincide with phylogenetic history and thus is often unknown. Among such sequences, the capacity for transcriptional activity coupled with the adaptive use of reverse transcription can lead to a diverse group of genomic elements across taxa, otherwise known as selfish elements or mobile elements. Short interspersed nuclear elements (SINEs) are nonautonomous mobile elements found in eukaryotic genomes, typically derived from cellular RNAs such as tRNAs, 7SL or 5S rRNA. Here, we identify and characterize a previously unknown SINE derived from the 3'-end of the large ribosomal subunit (LSU or 28S rDNA) and transcribed via RNA polymerase III. This new element, SINE28, is represented in low-copy numbers in the human reference genome assembly, wherein we have identified 27 discrete loci. Phylogenetic analysis indicates these elements have been transpositionally active within primate lineages as recently as 6 MYA while modern humans still carry transcriptionally active copies. Moreover, we have identified SINE28s in all currently available assembled mammalian genome sequences. Phylogenetic comparisons indicate that these elements are frequently rederived from the highly conserved LSU rRNA sequences in a lineage-specific manner. We propose that this element has not been previously recognized as a SINE given its high identity to the canonical LSU, and that SINE28 likely represents one of possibly many unidentified, active transposable elements within mammalian genomes.

Maternal methyl supplemented diets and effects on offspring health

Women seeking to become pregnant and pregnant women are currently advised to consume high amounts of folic acid and other methyl donors to prevent neural tube defects in their offspring. These diets can alter methylation patterns of several biomolecules, including nucleic acids, and histone proteins. Limited animal model data suggests that developmental exposure to these maternal methyl supplemented (MS) diets leads to beneficial epimutations. However, other rodent and humans studies have yielded opposing findings with such diets leading to promiscuous epimutations that are likely associated with negative health outcomes. Conflict exists to whether these maternal diets are preventative or exacerbate the risk for Autism Spectrum Disorders (ASD) in children. This review will discuss the findings to date on the potential beneficial and aversive effects of maternal MS diets. We will also consider how other factors might influence the effects of MS diets. Current data suggest that there is cause for concern as maternal MS diets may lead to epimutations that underpin various diseases, including neurobehavioral disorders. Further studies are needed to explore the comprehensive effects maternal MS diets have on the offspring epigenome and subsequent overall health.

Transcriptome resources for the white-footed mouse (Peromyscus leucopus): new genomic tools for investigating ecologically divergent urban and rural populations

Genomic resources are important and attainable for examining evolutionary change in divergent natural populations of nonmodel species. We utilized two next-generation sequencing (NGS) platforms, 454 and SOLiD 5500XL, to assemble low-coverage transcriptomes of the white-footed mouse (Peromyscus leucopus), a widespread and abundant native rodent in eastern North America. We sequenced liver mRNA transcripts from multiple individuals collected from urban populations in New York City and rural populations in undisturbed protected areas nearby and assembled a reference transcriptome using 1 080 065 954 SOLiD 5500XL (75 bp) reads and 3 052 640 454 GS FLX + reads. The reference contained 40 908 contigs with a N50 = 1044 bp and a total content of 30.06 Megabases (Mb). Contigs were annotated from Mus musculus (39.96% annotated) Uniprot databases. We identified 104 655 high-quality single nucleotide polymorphisms (SNPs) and 65 single sequence repeats (SSRs) with flanking primers. We also used normalized read counts to identify putative gene expression differences in 10 genes between populations. There were 19 contigs significantly differentially expressed in urban populations compared to rural populations, with gene function annotations generally related to the translation and modification of proteins and those involved in immune responses. The individual transcriptomes generated in this study will be used to investigate evolutionary responses to urbanization. The reference transcriptome provides a valuable resource for the scientific community using North American Peromyscus species as emerging model systems for ecological genetics and adaptation.

Relaxed selection causes microevolution of seawater osmoregulation and gene expression in landlocked Alewives

Ecological transitions from marine to freshwater environments have been important in the creation of diversity among fishes. Evolutionary changes associated with these transitions likely involve modifications of osmoregulatory function. In particular, relaxed selection on hypo-osmoregulation should strongly affect animals that transition into novel freshwater environments. We used populations of the Alewife (Alosa pseudoharengus) to study evolutionary shifts in hypo-osmoregulatory capacity and ion regulation associated with freshwater transitions. Alewives are ancestrally anadromous, but multiple populations in Connecticut have been independently restricted to freshwater lakes; these landlocked populations complete their entire life cycle in freshwater. Juvenile landlocked and anadromous Alewives were exposed to three salinities (1, 20 and 30 ppt) in small enclosures within the lake. We detected strong differentiation between life history forms: landlocked Alewives exhibited reduced seawater tolerance and hypo-osmoregulatory performance compared to anadromous Alewives. Furthermore, gill Na(+)/K(+)-ATPase activity and transcription of genes for seawater osmoregulation (NKCC-Na(+)/K(+)/2Cl(-) cotransporter and CFTR-cystic fibrosis transmembrane conductance regulator) exhibited reduced responsiveness to seawater challenge. Our study demonstrates that adaptations of marine-derived species to completely freshwater life cycles involve partial loss of seawater osmoregulatory performance mediated through changes to ion regulation in the gill.

A genetic map of Peromyscus with chromosomal assignment of linkage groups (a Peromyscus genetic map)

The rodent genus Peromyscus is the most numerous and species-rich mammalian group in North America. The naturally occurring diversity within this genus allows opportunities to investigate the genetic basis of adaptation, monogamy, behavioral and physiological phenotypes, growth control, genomic imprinting, and disease processes. Increased genomic resources including a high quality genetic map are needed to capitalize on these opportunities. We produced interspecific hybrids between the prairie deer mouse (P. maniculatus bairdii) and the oldfield mouse (P. polionotus) and scored meiotic recombination events in backcross progeny. A genetic map was constructed by genotyping of backcross progeny at 185 gene-based and 155 microsatellite markers representing all autosomes and the X-chromosome. Comparison of the constructed genetic map with the molecular maps of Mus and Rattus and consideration of previous results from interspecific reciprocal whole chromosome painting allowed most linkage groups to be unambiguously assigned to specific Peromyscus chromosomes. Based on genomic comparisons, this Peromyscus genetic map covers ~83% of the Rattus genome and 79% of the Mus genome. This map supports previous results that the Peromyscus genome is more similar to Rattus than Mus. For example, coverage of the 20 Rattus autosomes and the X-chromosome is accomplished with only 28 segments of the Peromyscus map, but coverage of the 19 Mus autosomes and the X-chromosome requires 40 chromosomal segments of the Peromyscus map. Furthermore, a single Peromyscus linkage group corresponds to about 91% of the rat and only 76% of the mouse X-chromosomes.

Retroelement demethylation associated with abnormal placentation in Mus musculus x Mus caroli hybrids

The proper functioning of the placenta requires specific patterns of methylation and the appropriate regulation of retroelements, some of which have been co-opted by the genome for placental-specific gene expression. Our inquiry was initiated to determine the causes of the placental defects observed in crosses between two species of mouse, Mus musculus and Mus caroli. M. musculus × M. caroli fetuses are rarely carried to term, possibly as a result of genomic incompatibility in the placenta. Taking into account that placental dysplasia is observed in Peromyscus and other Mus hybrids, and that endogenous retroviruses are expressed in the placental transcriptome, we hypothesized that these placental defects could result, in part, from failure of the genome defense mechanism, DNA methylation, to regulate the expression of retroelements. Hybrid M. musculus × M. caroli embryos were produced by artificial insemination, and dysplastic placentas were subjected to microarray and methylation screens. Aberrant overexpression of an X-linked Mus retroelement in these hybrid placentas is consistent with local demethylation of this retroelement, concomitant with genome instability, disruption of gene regulatory pathways, and dysgenesis. We propose that the placenta is a specific site of control that is disrupted by demethylation and retroelement activation in interspecific hybridization that occur as a result of species incompatibility of methylation machinery. To our knowledge, the present data provide the first report of retroelement activation linked to decreased methylation in a eutherian hybrid system.

Making a long story short: noncoding RNAs and chromosome change

As important as the events that influence selection for specific chromosome types in the derivation of novel karyotypes, are the events that initiate the changes in chromosome number and structure between species, and likewise polymorphisms, variants and disease states within species. Although once thought of as transcriptional 'noise', noncoding RNAs (ncRNAs) are now recognized as important mediators of epigenetic regulation and chromosome stability. Here we highlight recent work that illustrates the influence short and long ncRNAs have as participants in the function and stability of chromosome regions such as centromeres, telomeres, evolutionary breakpoints and fragile sites. We summarize recent evidence that ncRNAs can facilitate chromosome change and present mechanisms by which ncRNAs create DNA breaks. Finally, we present hypotheses on how they may create novel karyotypes and thus affect chromosome evolution.

Genome sequence of an Australian kangaroo, Macropus eugenii, provides insight into the evolution of mammalian reproduction and development

BACKGROUND

We present the genome sequence of the tammar wallaby, Macropus eugenii, which is a member of the kangaroo family and the first representative of the iconic hopping mammals that symbolize Australia to be sequenced. The tammar has many unusual biological characteristics, including the longest period of embryonic diapause of any mammal, extremely synchronized seasonal breeding and prolonged and sophisticated lactation within a well-defined pouch. Like other marsupials, it gives birth to highly altricial young, and has a small number of very large chromosomes, making it a valuable model for genomics, reproduction and development.

RESULTS

The genome has been sequenced to 2 × coverage using Sanger sequencing, enhanced with additional next generation sequencing and the integration of extensive physical and linkage maps to build the genome assembly. We also sequenced the tammar transcriptome across many tissues and developmental time points. Our analyses of these data shed light on mammalian reproduction, development and genome evolution: there is innovation in reproductive and lactational genes, rapid evolution of germ cell genes, and incomplete, locus-specific X inactivation. We also observe novel retrotransposons and a highly rearranged major histocompatibility complex, with many class I genes located outside the complex. Novel microRNAs in the tammar HOX clusters uncover new potential mammalian HOX regulatory elements.

CONCLUSIONS

Analyses of these resources enhance our understanding of marsupial gene evolution, identify marsupial-specific conserved non-coding elements and critical genes across a range of biological systems, including reproduction, development and immunity, and provide new insight into marsupial and mammalian biology and genome evolution.

Centromere conversion and retention in somatic cell hybrids

The generation of somatic cell hybridization-derived cell lines between highly divergent species affords the opportunity to examine the concept of 'genome dominance' in the context of genetic and epigenetic changes. While whole-scale genome dominance has been well documented in natural hybrids among closely related species, an examination of centromere position and sequence retention in 2 marsupial-eutherian hybrids has revealed a mechanism for 'centromere dominance' as a driving force in the generation of stable somatic cell hybrids following an initial period of genomic instability. While one somatic cell hybrid cell line appeared to retain marsupial centromere sequences which remained competent to recruit the centromere-specific histone variant CENP-A in a Chinese hamster background, fusion events between marsupial and mouse-derived chromosomes in another hybrid line led to a centromere sequence conversion from one species to the other. We postulate that the necessity to maintain an epigenetically defined centromere following genome hybridization may be responsible for retention of specific chromosomes and may result in rapid sequence turnover to facilitate the recruitment of CENP-A containing histones.

Enhancing genome assemblies by integrating non-sequence based data

Introduction

Many genome projects were underway before the advent of high-throughput sequencing and have thus been supported by a wealth of genome information from other technologies. Such information frequently takes the form of linkage and physical maps, both of which can provide a substantial amount of data useful in de novo sequencing projects. Furthermore, the recent abundance of genome resources enables the use of conserved synteny maps identified in related species to further enhance genome assemblies.

Methods

The tammar wallaby (Macropus eugenii) is a model marsupial mammal with a low coverage genome. However, we have access to extensive comparative maps containing over 14,000 markers constructed through the physical mapping of conserved loci, chromosome painting and comprehensive linkage maps. Using a custom Bioperl pipeline, information from the maps was aligned to assembled tammar wallaby contigs using BLAT. This data was used to construct pseudo paired-end libraries with intervals ranging from 5-10 MB. We then used Bambus (a program designed to scaffold eukaryotic genomes by ordering and orienting contigs through the use of paired-end data) to scaffold our libraries. To determine how map data compares to sequence based approaches to enhance assemblies, we repeated the experiment using a 0.5× coverage of unique reads from 4 KB and 8 KB Illumina paired-end libraries. Finally, we combined both the sequence and non-sequence-based data to determine how a combined approach could further enhance the quality of the low coverage de novo reconstruction of the tammar wallaby genome.

Results

Using the map data alone, we were able order 2.2% of the initial contigs into scaffolds, and increase the N50 scaffold size to 39 KB (36 KB in the original assembly). Using only the 0.5× paired-end sequence based data, 53% of the initial contigs were assigned to scaffolds. Combining both data sets resulted in a further 2% increase in the number of initial contigs integrated into a scaffold (55% total) but a 35% increase in N50 scaffold size over the use of sequence-based data alone.

Conclusions

We provide a relatively simple pipeline utilizing existing bioinformatics tools to integrate map data into a genome assembly which is available at http://www.mcb.uconn.edu/fac.php?name=paska. While the map data only contributed minimally to assigning the initial contigs to scaffolds in the new assembly, it greatly increased the N50 size. This process added structure to our low coverage assembly, greatly increasing its utility in further analyses.

Abundant human DNA contamination identified in non-primate genome databases

During routine screens of the NCBI databases using human repetitive elements we discovered an unlikely level of nucleotide identity across a broad range of phyla. To ascertain whether databases containing DNA sequences, genome assemblies and trace archive reads were contaminated with human sequences, we performed an in depth search for sequences of human origin in non-human species. Using a primate specific SINE, AluY, we screened 2,749 non-primate public databases from NCBI, Ensembl, JGI, and UCSC and have found 492 to be contaminated with human sequence. These represent species ranging from bacteria (B. cereus) to plants (Z. mays) to fish (D. rerio) with examples found from most phyla. The identification of such extensive contamination of human sequence across databases and sequence types warrants caution among the sequencing community in future sequencing efforts, such as human re-sequencing. We discuss issues this may raise as well as present data that gives insight as to how this may be occurring.

Chromosomes, conflict, and epigenetics: chromosomal speciation revisited

Since Darwin first noted that the process of speciation was indeed the "mystery of mysteries," scientists have tried to develop testable models for the development of reproductive incompatibilities-the first step in the formation of a new species. Early theorists proposed that chromosome rearrangements were implicated in the process of reproductive isolation; however, the chromosomal speciation model has recently been questioned. In addition, recent data from hybrid model systems indicates that simple epistatic interactions, the Dobzhansky-Muller incompatibilities, are more complex. In fact, incompatibilities are quite broad, including interactions among heterochromatin, small RNAs, and distinct, epigenetically defined genomic regions such as the centromere. In this review, we will examine both classical and current models of chromosomal speciation and describe the "evolving" theory of genetic conflict, epigenetics, and chromosomal speciation.

A unique late-replicating XY to autosome translocation in Peromyscus melanophrys

We report on the characterization of the Peromyscus melanophrys karyotype and sex chromosome system. Classic studies reported the sex chromosome system of this species may be as complex as an X(1)X(1)X(2)X(2)/X(1)X(2)Y(1)Y(2) and provided conflicting identification of the X chromosome. Using Peromyscus maniculatus chromosome paints, we have positively identified the sex chromosomes and clarified the sex determining system that once perplexed Peromyscus researchers. The sex chromosomes are characterized by a unique autosomal translocation of DNA shared between both the X and Y chromosomes. The translocated material is late replicating and heterochromatic yet retains the active chromatin conformation. Thus, autosomal regions derived from translocations involving repeat-rich material may retain some epigenetic marks specific to the sex chromosomes despite loss of epigenetic silencing activity.

The role of ncRNA in centromeres: a lesson from marsupials

Though centromeres have been thought to be comprised of repetitive, transcriptionally inactive DNA, new evidence suggests that eukaryotic centromeres produce a variety of transcripts and that RNA is essential for centromere competence. It has been proposed that centromere satellite transcripts play an essential role in centromere function through demarcation of the kinetochore-binding domain. However, the regional limits and regulation of transcription within the mammalian centromere are unknown. Analysis of transcriptional domains within the centromere in mammalian models is impeded by the unbridgeable expanse of satellite monomers throughout the pericentromere. The comparatively small size of the wallaby centromere and the evolutionary role of the centromere in marsupial speciation events position the wallaby centromere as a tractable and valuable mammalian centromere model. We highlight the current understanding of the wallaby centromere and the role of transcription in centromere function.