Genomic Characterization of Large Heterochromatic Gaps in the Human Genome Assembly

Abstract

The largest gaps in the human genome assembly correspond to multi-megabase heterochromatic regions composed primarily of two related families of tandem repeats, Human Satellites 2 and 3 (HSat2,3). The abundance of repetitive DNA in these regions challenges standard mapping and assembly algorithms, and as a result, the sequence composition and potential biological functions of these regions remain largely unexplored. Furthermore, existing genomic tools designed to predict consensus-based descriptions of repeat families cannot be readily applied to complex satellite repeats such as HSat2,3, which lack a consistent repeat unit reference sequence. Here we present an alignment-free method to characterize complex satellites using whole-genome shotgun read datasets. Utilizing this approach, we classify HSat2,3 sequences into fourteen subfamilies and predict their chromosomal distributions, resulting in a comprehensive satellite reference database to further enable genomic studies of heterochromatic regions. We also identify 1.3 Mb of non-repetitive sequence interspersed with HSat2,3 across 17 unmapped assembly scaffolds, including eight annotated gene predictions. Finally, we apply our satellite reference database to high-throughput sequence data from 396 males to estimate array size variation of the predominant HSat3 array on the Y chromosome, confirming that satellite array sizes can vary between individuals over an order of magnitude (7 to 98 Mb) and further demonstrating that array sizes are distributed differently within distinct Y haplogroups. In summary, we present a novel framework for generating initial reference databases for unassembled genomic regions enriched with complex satellite DNA, and we further demonstrate the utility of these reference databases for studying patterns of sequence variation within human populations.

Multicolor CRISPR labeling of chromosomal loci in human cells

“…There is resurgent interest in the repeated sequence elements of the human genome as their "post-Human Genome Project" characterization continues (21). For many applications, we envision the need to identify sets of repeated DNA sequences that are unique to a single locus on a chromosome or are present at multiple sites on only one chromosome in the complement (i.e., as a bar code).…”

‘Satellite DNA transcripts have diverse biological roles in Drosophila’

“…But the observation that at least a subset of satDNA families are involved in important biological roles makes the study of these highly abundant and fast-evolving components of the eukaryote genome highly relevant for structural, functional and evolutionary genomics. It is also worth mentioning that in the past few years, new and efficient bioinformatic tools are becoming available for the identification of satDNAs from sequenced genomes (Novák et al, 2014), whereas longtemplate sequencing and new computational approaches (Altemose et al, 2014) have been fostering the assembly of satDNA repeats. With all these new findings and tools, genomic studies, including those related to wholesequenced genomes, will certainly benefit by a 'satDNA recall'.…”

Influence of Data Geometry in Random Subset Feature Selection

“…The normal distribution is a real valued distribution whose mean, median and mode are equal [2]. This is symmetrical about mean; and not suitable for the variables which are inherently positive or negatively skewed.…”