CRISPR-based targeted haplotype-resolved assembly of a megabase region.

Constructing high-quality haplotype-resolved genome assemblies has substantially improved the ability to detect and characterize genetic variants. A targeted approach providing readily access to the rich information from haplotype-resolved genome assemblies will be appealing to groups of basic researchers and medical scientists focused on specific genomic regions. Here, using the 4.5 megabase, notoriously difficult-to-assemble major histocompatibility complex (MHC) region as an example, we demonstrated an approach to construct haplotype-resolved assembly of the targeted genomic region with the CRISPR-based enrichment. Compared to the results from haplotype-resolved genome assembly, our targeted approach achieved comparable completeness and accuracy with reduced computing complexity, sequencing cost, as well as the amount of starting materials. Moreover, using the targeted assembled personal MHC haplotypes as the reference both improves the quantification accuracy for sequencing data and enables allele-specific functional genomics analyses of the MHC region. Given its highly efficient use of resources, our approach can greatly facilitate population genetic studies of targeted regions, and may pave a new way to elucidate the molecular mechanisms in disease etiology.

Altered DNA methylation pattern reveals epigenetic regulation of Hox genes in thoracic aortic dissection and serves as a biomarker in disease diagnosis.

BACKGROUND: Thoracic aortic dissection (TAD) is a severe disease with limited understandings in its pathogenesis. Altered DNA methylation has been revealed to be involved in many diseases etiology. Few studies have examined the role of DNA methylation in the development of TAD. This study explored alterations of the DNA methylation landscape in TAD and examined the potential role of cell-free DNA (cfDNA) methylation as a biomarker in TAD diagnosis. RESULTS: Ascending aortic tissues from TAD patients (Stanford type A; n = 6) and healthy controls (n = 6) were first examined via whole-genome bisulfite sequencing (WGBS). While no obvious global methylation shift was observed, numerous differentially methylated regions (DMRs) were identified, with associated genes enriched in the areas of vasculature and heart development. We further confirmed the methylation and expression changes in homeobox (Hox) clusters with 10 independent samples using bisulfite pyrosequencing and quantitative real-time PCR (qPCR). Among these, HOXA5, HOXB6 and HOXC6 were significantly down-regulated in TAD samples relative to controls. To evaluate cfDNA methylation pattern as a biomarker in TAD diagnosis, cfDNA from TAD patients (Stanford type A; n = 7) and healthy controls (n = 4) were examined by WGBS. A prediction model was built using DMRs identified previously from aortic tissues on methylation data from cfDNA. Both high sensitivity (86%) and specificity (75%) were achieved in patient classification (AUC = 0.96). CONCLUSIONS: These findings showed an altered epigenetic regulation in TAD patients. This altered epigenetic regulation and subsequent altered expression of genes associated with vasculature and heart development, such as Hox family genes, may contribute to the loss of aortic integrity and TAD pathogenesis. Additionally, the cfDNA methylation in TAD was highly disease specific, which can be used as a non-invasive biomarker for disease prediction.

A Hu sheep genome with the first ovine Y chromosome reveal introgression history after sheep domestication.

The Y chromosome plays key roles in male fertility and reflects the evolutionary history of paternal lineages. Here, we present a de novo genome assembly of the Hu sheep with the first draft assembly of ovine Y chromosome (oMSY), using nanopore sequencing and Hi-C technologies. The oMSY that we generated spans 10.6 Mb from which 775 Y-SNPs were identified by applying a large panel of whole genome sequences from worldwide sheep and wild Iranian mouflons. Three major paternal lineages (HY1a, HY1b and HY2) were defined across domestic sheep, of which HY2 was newly detected. Surprisingly, HY2 forms a monophyletic clade with the Iranian mouflons and is highly divergent from both HY1a and HY1b. Demographic analysis of Y chromosomes, mitochondrial and nuclear genomes confirmed that HY2 and the maternal counterpart of lineage C represented a distinct wild mouflon population in Iran that diverge from the direct ancestor of domestic sheep, the wild mouflons in Southeastern Anatolia. Our results suggest that wild Iranian mouflons had introgressed into domestic sheep and thereby introduced this Iranian mouflon specific lineage carrying HY2 to both East Asian and Africa sheep populations.

The Nuclear Matrix Protein SAFB Cooperates with Major Satellite RNAs to Stabilize Heterochromatin Architecture Partially through Phase Separation.

Interphase chromatin is hierarchically organized into higher-order architectures that are essential for gene functions, yet the biomolecules that regulate these 3D architectures remain poorly understood. Here, we show that scaffold attachment factor B (SAFB), a nuclear matrix (NM)-associated protein with RNA-binding functions, modulates chromatin condensation and stabilizes heterochromatin foci in mouse cells. SAFB interacts via its R/G-rich region with heterochromatin-associated repeat transcripts such as major satellite RNAs, which promote the phase separation driven by SAFB. Depletion of SAFB leads to changes in 3D genome organization, including an increase in interchromosomal interactions adjacent to pericentromeric heterochromatin and a decrease in genomic compartmentalization, which could result from the decondensation of pericentromeric heterochromatin. Collectively, we reveal the integrated roles of NM-associated proteins and repeat RNAs in the 3D organization of heterochromatin, which may shed light on the molecular mechanisms of nuclear architecture organization.

Building a sequence map of the pig pan-genome from multiple de novo assemblies and Hi-C data.

Pigs were domesticated independently in the Near East and China, indicating that a single reference genome from one individual is unable to represent the full spectrum of divergent sequences in pigs worldwide. Therefore, 12 de novo pig assemblies from Eurasia were compared in this study to identify the missing sequences from the reference genome. As a result, 72.5 Mb of non-redundant sequences (∼3% of the genome) were found to be absent from the reference genome (Sscrofa11.1) and were defined as pan-sequences. Of the pan-sequences, 9.0 Mb were dominant in Chinese pigs, in contrast with their low frequency in European pigs. One sequence dominant in Chinese pigs contained the complete genic region of the tazarotene-induced gene 3 (TIG3) gene which is involved in fatty acid metabolism. Using flanking sequences and Hi-C based methods, 27.7% of the sequences could be anchored to the reference genome. The supplementation of these sequences could contribute to the accurate interpretation of the 3D chromatin structure. A web-based pan-genome database was further provided to serve as a primary resource for exploration of genetic diversity and promote pig breeding and biomedical research.

Towards the Complete Goat Pan-Genome by Recovering Missing Genomic Segments From the Reference Genome.

It is broadly expected that next generation sequencing will ultimately generate a complete genome as is the latest goat reference genome (ARS1), which is considered to be one of the most continuous assemblies in livestock. However, the rich diversity of worldwide goat breeds indicates that a genome from one individual would be insufficient to represent the whole genomic contents of goats. By comparing nine de novo assemblies from seven sibling species of domestic goat with ARS1 and using resequencing and transcriptome data from goats for verification, we identified a total of 38.3 Mb sequences that were absent in ARS1. The pan-sequences contain genic fractions with considerable expression. Using the pan-genome (ARS1 together with the pan-sequences) as a reference genome, variation calling efficacy can be appreciably improved. A total of 56,657 spurious SNPs per individual were repressed and 24,414 novel SNPs per individual on average were recovered as a result of better reads mapping quality. The transcriptomic mapping rate was also increased by ∼1.15%. Our study demonstrated that comparing de novo assemblies from closely related species is an efficient and reliable strategy for finding missing sequences from the reference genome and could be applicable to other species. Pan-genome can serve as an improved reference genome in animals for a better exploration of the underlying genomic variations and could increase the probability of finding genotype-phenotype associations assessed by a comprehensive variation database containing much more differences between individuals. We have constructed a goat pan-genome web interface for data visualization (http://animal.nwsuaf.edu.cn/panGoat).

Determination of local chromatin interactions using a combined CRISPR and peroxidase APEX2 system.

The architecture and function of chromatin are largely regulated by local interacting molecules, such as transcription factors and noncoding RNAs. However, our understanding of these regulatory molecules at a given locus is limited because of technical difficulties. Here, we describe the use of Clustered Regularly Interspaced Short Palindromic Repeats and an engineered ascorbate peroxidase 2 (APEX2) system to investigate local chromatin interactions (CAPLOCUS). We showed that with specific small-guide RNA targets, CAPLOCUS could efficiently identify both repetitive genomic regions and single-copy genomic locus with high resolution. Genome-wide sequencing revealed known and potential long-range chromatin interactions for a specific single-copy locus. CAPLOCUS also identified telomere-associated RNAs. CAPLOCUS, followed by mass spectrometry, identified both known and novel telomere-associated proteins in their native states. Thus, CAPLOCUS may be a useful approach for studying local interacting molecules at any given chromosomal location.