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.

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.

Whole-genome sequence and assembly of the Javan gibbon (Hylobates moloch).

The Javan gibbon, Hylobates moloch, is an endangered gibbon species restricted to the forest remnants of western and central Java, Indonesia, and one of the rarest of the Hylobatidae family. Hylobatids consist of 4 genera (Holoock, Hylobates, Symphalangus, and Nomascus) that are characterized by different numbers of chromosomes, ranging from 38 to 52. The underlying cause of this karyotype plasticity is not entirely understood, at least in part, due to the limited availability of genomic data. Here we present the first scaffold-level assembly for H. moloch using a combination of whole-genome Illumina short reads, 10X Chromium linked reads, PacBio, and Oxford Nanopore long reads and proximity-ligation data. This Hylobates genome represents a valuable new resource for comparative genomics studies in primates.

Co-option of the lineage-specific LAVA retrotransposon in the gibbon genome.

Co-option of transposable elements (TEs) to become part of existing or new enhancers is an important mechanism for evolution of gene regulation. However, contributions of lineage-specific TE insertions to recent regulatory adaptations remain poorly understood. Gibbons present a suitable model to study these contributions as they have evolved a lineage-specific TE called LAVA (LINE-AluSz-VNTR-AluLIKE), which is still active in the gibbon genome. The LAVA retrotransposon is thought to have played a role in the emergence of the highly rearranged structure of the gibbon genome by disrupting transcription of cell cycle genes. In this study, we investigated whether LAVA may have also contributed to the evolution of gene regulation by adopting enhancer function. We characterized fixed and polymorphic LAVA insertions across multiple gibbons and found 96 LAVA elements overlapping enhancer chromatin states. Moreover, LAVA was enriched in multiple transcription factor binding motifs, was bound by an important transcription factor (PU.1), and was associated with higher levels of gene expression in cis We found gibbon-specific signatures of purifying/positive selection at 27 LAVA insertions. Two of these insertions were fixed in the gibbon lineage and overlapped with enhancer chromatin states, representing putative co-opted LAVA enhancers. These putative enhancers were located within genes encoding SETD2 and RAD9A, two proteins that facilitate accurate repair of DNA double-strand breaks and prevent chromosomal rearrangement mutations. Co-option of LAVA in these genes may have influenced regulation of processes that preserve genome integrity. Our findings highlight the importance of considering lineage-specific TEs in studying evolution of gene regulatory elements.

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.

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.

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.

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.

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.

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.

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.

Transcriptomic imprints of adaptation to fresh water: parallel evolution of osmoregulatory gene expression in the Alewife.

Comparative approaches in physiological genomics offer an opportunity to understand the functional importance of genes involved in niche exploitation. We used populations of Alewife (Alosa pseudoharengus) to explore the transcriptional mechanisms that underlie adaptation to fresh water. Ancestrally anadromous Alewives have recently formed multiple, independently derived, landlocked populations, which exhibit reduced tolerance of saltwater and enhanced tolerance of fresh water. Using RNA-seq, we compared transcriptional responses of an anadromous Alewife population to two landlocked populations after acclimation to fresh (0 ppt) and saltwater (35 ppt). Our results suggest that the gill transcriptome has evolved in primarily discordant ways between independent landlocked populations and their anadromous ancestor. By contrast, evolved shifts in the transcription of a small suite of well-characterized osmoregulatory genes exhibited a strong degree of parallelism. In particular, transcription of genes that regulate gill ion exchange has diverged in accordance with functional predictions: freshwater ion-uptake genes (most notably, the 'freshwater paralog' of Na+ /K+ -ATPase α-subunit) were more highly expressed in landlocked forms, whereas genes that regulate saltwater ion secretion (e.g. the 'saltwater paralog' of NKAα) exhibited a blunted response to saltwater. Parallel divergence of ion transport gene expression is associated with shifts in salinity tolerance limits among landlocked forms, suggesting that changes to the gill's transcriptional response to salinity facilitate freshwater adaptation.

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.

Hypermorphic expression of centromeric retroelement-encoded small RNAs impairs CENP-A loading.

The proper functioning of centromeres requires a complex cascade of epigenetic events involving chromatin and kinetochore assembly; however, the precise mechanism by which this cascade proceeds is unknown. The pivotal event during kinetochore formation is the "loading," or deposition, of CENP-A. This histone H3 variant is specific to centromeres and replaces conventional H3 in centromeric chromatin. Failure to load CENP-A into mammalian centromeres in late telophase/early G1 of the cell cycle leads to malsegregation and cell division defects in subsequent cell cycles. Mounting evidence supports the hypothesis that an RNA component is involved, although how RNAs participate in centromere formation in mammals has remained unknown. Using the marsupial model, the tammar wallaby, we show that centromeric retroelements produce small RNAs and that hypermorphic expression of these centromeric small RNAs results in disruption of CENP-A localization. We propose that tight regulation of the processing of this new class of small RNAs, crasiRNAs, is an integral component of the epigenetic framework necessary for centromere establishment.

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.

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.

Conserving a threatened North American walnut: a chromosome-scale reference genome for butternut (Juglans cinerea).

With the advent of affordable and more accurate third-generation sequencing technologies, and the associated bioinformatic tools, it is now possible to sequence, assemble, and annotate more species of conservation concern than ever before. Juglans cinerea, commonly known as butternut or white walnut, is a member of the walnut family, native to the Eastern United States and Southeastern Canada. The species is currently listed as Endangered on the IUCN Red List due to decline from an invasive fungus known as Ophiognomonia clavigignenti-juglandacearum (Oc-j) that causes butternut canker. Oc-j creates visible sores on the trunks of the tree which essentially starves and slowly kills the tree. Natural resistance to this pathogen is rare. Conserving butternut is of utmost priority due to its critical ecosystem role and cultural significance. As part of an integrated undergraduate and graduate student training program in biodiversity and conservation genomics, the first reference genome for Juglans cinerea is described here. This chromosome-scale 539 Mb assembly was generated from over 100 × coverage of Oxford Nanopore long reads and scaffolded with the Juglans mandshurica genome. Scaffolding with a closely related species oriented and ordered the sequences in a manner more representative of the structure of the genome without altering the sequence. Comparisons with sequenced Juglandaceae revealed high levels of synteny and further supported J. cinerea's recent phylogenetic placement. Comparative assessment of gene family evolution revealed a significant number of contracting families, including several associated with biotic stress response.