Multifaceted Hi-C benchmarking: what makes a difference in chromosome-scale genome scaffolding?

A high-quality genome assembly and annotation of Thielaviopsis punctulata DSM102798

Black scorch disease (BSD), caused by the fungal pathogen Thielaviopsis punctulata (Tp) DSM102798, poses a significant threat to date palm cultivation in the United Arab Emirates (UAE). In this study, Chicago and Hi-C libraries were prepared as input for the Dovetail HiRise pipeline to scaffold the genome of Tp DSM102798. We generated an assembly with a total length of 28.23 Mb comprising 1,256 scaffolds, and the assembly had a contig N50 of 18.56 kb, L50 of three, and a BUSCO completeness score of 98.6% for 758 orthologous genes. Annotation of this assembly produced 7,169 genes and 3,501 Gene Ontology (GO) terms. Compared to five other Thielaviopsis genomes, Tp DSM102798 exhibited the highest continuity with a cumulative size of 27.598 Mb for the first seven scaffolds, surpassing the assemblies of all examined strains. These findings offer a foundation for targeted strategies that enhance date palm resistance against BSD, and foster more sustainable and resilient agricultural systems.

Chromosomal DNA sequences of the Pacific saury genome: versatile resources for fishery science and comparative biology

Pacific saury (Cololabis saira) is a commercially important small pelagic fish species in Asia. In this study, we conducted the first-ever whole genome sequencing of this species, with single molecule, real-time (SMRT) sequencing technology. The obtained high-fidelity (HiFi) long-read sequence data, which amount to ~30-folds of its haploid genome size that was measured with quantitative PCR (1.17 Gb), were assembled into contigs. Scaffolding with Hi-C reads yielded a whole genome assembly containing 24 chromosome-scale sequences, with a scaffold N50 length of 47.7 Mb. Screening of repetitive elements including telomeric repeats was performed to characterize possible factors that need to be resolved towards 'telomere-to-telomere' sequencing. The larger genome size than in medaka, a close relative in Beloniformes, is at least partly explained by larger repetitive element quantity, which is reflected in more abundant tRNAs, in the Pacific saury genome. Protein-coding regions were predicted using transcriptome data, which resulted in 22,274 components. Retrieval of Pacific saury homologs of aquaporin (AQP) genes known from other teleost fishes validated high completeness and continuity of the genome assembly. These resources are available at https://treethinkers.nig.ac.jp/saira/ and will assist various molecular-level studies in fishery science and comparative biology.

Elasmobranch genome sequencing reveals evolutionary trends of vertebrate karyotype organization

Genomic studies of vertebrate chromosome evolution have long been hindered by the scarcity of chromosome-scale DNA sequences of some key taxa. One of those limiting taxa has been the elasmobranchs (sharks and rays), which harbor species often with numerous chromosomes and enlarged genomes. Here, we report the chromosome-scale genome assembly for the zebra shark Stegostoma tigrinum, an endangered species that has a relatively small genome among sharks (3.71 Gb), as well as for the whale shark Rhincodon typus Our analysis using a male-female comparison identified an X Chromosome, the first genomically characterized shark sex chromosome. The X Chromosome harbors the Hox C cluster whose intact linkage has not been shown for an elasmobranch fish. The sequenced shark genomes show a gradualism of chromosome length with remarkable length-dependent characteristics-shorter chromosomes tend to have higher GC content, gene density, synonymous substitution rate, and simple tandem repeat content as well as smaller gene length and lower interspersed repeat content. We challenge the traditional binary classification of karyotypes as with and without so-called microchromosomes. Even without microchromosomes, the length-dependent characteristics persist widely in nonmammalian vertebrates. Our investigation of elasmobranch karyotypes underpins their unique characteristics and provides clues for understanding how vertebrate karyotypes accommodate intragenomic heterogeneity to realize a complex readout. It also paves the way to dissecting more genomes with variable sizes to be sequenced at high quality.

GreenHill: a de novo chromosome-level scaffolding and phasing tool using Hi-C

Chromosome-level haplotype-resolved genome assembly is an important resource in molecular biology. However, current de novo haplotype assemblers require parental data or reference genomes and often fail to provide chromosome-level results. We present GreenHill, a novel scaffolding and phasing tool that considers various assemblers' contigs as input to reconstruct chromosome-level haplotypes using Hi-C without parental or reference data. Its unique functions include new error correction based on Hi-C contacts and the simultaneous use of Hi-C and long reads. Benchmarks reveal that GreenHill outperforms other approaches in contiguity and phasing accuracy, and the majority of chromosome arms are entirely phased.

Epigenetic plasticity safeguards heterochromatin configuration in mammals

Heterochromatin is a key architectural feature of eukaryotic chromosomes critical for cell type-specific gene expression and genome stability. In the mammalian nucleus, heterochromatin segregates from transcriptionally active genomic regions and exists in large, condensed, and inactive nuclear compartments. However, the mechanisms underlying the spatial organization of heterochromatin need to be better understood. Histone H3 lysine 9 trimethylation (H3K9me3) and lysine 27 trimethylation (H3K27me3) are two major epigenetic modifications that enrich constitutive and facultative heterochromatin, respectively. Mammals have at least five H3K9 methyltransferases (SUV39H1, SUV39H2, SETDB1, G9a and GLP) and two H3K27 methyltransferases (EZH1 and EZH2). In this study, we addressed the role of H3K9 and H3K27 methylation in heterochromatin organization using a combination of mutant cells for five H3K9 methyltransferases and an EZH1/2 dual inhibitor, DS3201. We showed that H3K27me3, which is normally segregated from H3K9me3, was redistributed to regions targeted by H3K9me3 after the loss of H3K9 methylation and that the loss of both H3K9 and H3K27 methylation resulted in impaired condensation and spatial organization of heterochromatin. Our data demonstrate that the H3K27me3 pathway safeguards heterochromatin organization after the loss of H3K9 methylation in mammalian cells.

A global high-density chromatin interaction network reveals functional long-range and trans-chromosomal relationships

BACKGROUND

Chromatin contacts are essential for gene-expression regulation; however, obtaining a high-resolution genome-wide chromatin contact map is still prohibitively expensive owing to large genome sizes and the quadratic scale of pairwise data. Chromosome conformation capture (3C)-based methods such as Hi-C have been extensively used to obtain chromatin contacts. However, since the sparsity of these maps increases with an increase in genomic distance between contacts, long-range or trans-chromatin contacts are especially challenging to sample.

RESULTS

Here, we create a high-density reference genome-wide chromatin contact map using a meta-analytic approach. We integrate 3600 human, 6700 mouse, and 500 fly Hi-C experiments to create species-specific meta-Hi-C chromatin contact maps with 304 billion, 193 billion, and 19 billion contacts in respective species. We validate that meta-Hi-C contact maps are uniquely powered to capture functional chromatin contacts in both cis and trans. We find that while individual dataset Hi-C networks are largely unable to predict any long-range coexpression (median 0.54 AUC), meta-Hi-C networks perform comparably in both cis and trans (0.65 AUC vs 0.64 AUC). Similarly, for long-range expression quantitative trait loci (eQTL), meta-Hi-C contacts outperform all individual Hi-C experiments, providing an improvement over the conventionally used linear genomic distance-based association. Assessing between species, we find patterns of chromatin contact conservation in both cis and trans and strong associations with coexpression even in species for which Hi-C data is lacking.

CONCLUSIONS

We have generated an integrated chromatin interaction network which complements a large number of methodological and analytic approaches focused on improved specificity or interpretation. This high-depth "super-experiment" is surprisingly powerful in capturing long-range functional relationships of chromatin interactions, which are now able to predict coexpression, eQTLs, and cross-species relationships. The meta-Hi-C networks are available at https://labshare.cshl.edu/shares/gillislab/resource/HiC/ .

Squalomix: shark and ray genome analysis consortium and its data sharing platform

The taxon Elasmobranchii (sharks and rays) contains one of the long-established evolutionary lineages of vertebrates with a tantalizing collection of species occupying critical aquatic habitats. To overcome the current limitation in molecular resources, we launched the Squalomix Consortium in 2020 to promote a genome-wide array of molecular approaches, specifically targeting shark and ray species. Among the various bottlenecks in working with elasmobranchs are their elusiveness and low fecundity as well as the large and highly repetitive genomes. Their peculiar body fluid composition has also hindered the establishment of methods to perform routine cell culturing required for their karyotyping. In the Squalomix consortium, these obstacles are expected to be solved through a combination of in-house cytological techniques including karyotyping of cultured cells, chromatin preparation for Hi-C data acquisition, and high fidelity long-read sequencing. The resources and products obtained in this consortium, including genome and transcriptome sequences, a genome browser powered by JBrowse2 to visualize sequence alignments, and comprehensive matrices of gene expression profiles for selected species are accessible through https://github.com/Squalomix/info.

Fully resolved assembly of Cryptosporidium parvum

Background

Cryptosporidium parvum is an apicomplexan parasite commonly found across many host species with a global infection prevalence in human populations of 7.6%. Understanding its diversity and genomic makeup can help in fighting established infections and prohibiting further transmission. The basis of every genomic study is a high-quality reference genome that has continuity and completeness, thus enabling comprehensive comparative studies.

Findings

Here, we provide a highly accurate and complete reference genome of Cryptosporidium parvum. The assembly is based on Oxford Nanopore reads and was improved using Illumina reads for error correction. We also outline how to evaluate and choose from different assembly methods based on 2 main approaches that can be applied to other Cryptosporidium species. The assembly encompasses 8 chromosomes and includes 13 telomeres that were resolved. Overall, the assembly shows a high completion rate with 98.4% single-copy BUSCO genes.

Conclusions

This high-quality reference genome of a zoonotic IIaA17G2R1 C. parvum subtype isolate provides the basis for subsequent comparative genomic studies across the Cryptosporidium clade. This will enable improved understanding of diversity, functional, and association studies.

Highly rigid H3.1/H3.2-H3K9me3 domains set a barrier for cell fate reprogramming in trophoblast stem cells

Here, Hada et al. comprehensively analyzed epigenomic features of mouse trophoblast stem cells (TSCs). They used genome-wide, high-throughput analyses to show that the TSC genome contains large-scale (>1-Mb) rigid heterochromatin architectures that have a high degree of histone H3.1/3.2–H3K9me3 accumulation, termed TSC-defined highly heterochromatinized domains (THDs), and are uniquely developed in placental lineage cells that serve to protect them from fate reprogramming to stably maintain placental function.

The placenta is a highly evolved, specialized organ in mammals. It differs from other organs in that it functions only for fetal maintenance during gestation. Therefore, there must be intrinsic mechanisms that guarantee its unique functions. To address this question, we comprehensively analyzed epigenomic features of mouse trophoblast stem cells (TSCs). Our genome-wide, high-throughput analyses revealed that the TSC genome contains large-scale (>1-Mb) rigid heterochromatin architectures with a high degree of histone H3.1/3.2–H3K9me3 accumulation, which we termed TSC-defined highly heterochromatinized domains (THDs). Importantly, depletion of THDs by knockdown of CAF1, an H3.1/3.2 chaperone, resulted in down-regulation of TSC markers, such as Cdx2 and Elf5, and up-regulation of the pluripotent marker Oct3/4, indicating that THDs maintain the trophoblastic nature of TSCs. Furthermore, our nuclear transfer technique revealed that THDs are highly resistant to genomic reprogramming. However, when H3K9me3 was removed, the TSC genome was fully reprogrammed, giving rise to the first TSC cloned offspring. Interestingly, THD-like domains are also present in mouse and human placental cells in vivo, but not in other cell types. Thus, THDs are genomic architectures uniquely developed in placental lineage cells, which serve to protect them from fate reprogramming to stably maintain placental function.

The chromosome-scale genome assembly, annotation and evolution of Rhododendron henanense subsp. lingbaoense

Rhododendron henanense subsp. lingbaoense (hereafter referred to as R. henanense) is an endemic species naturally distributed in the Henan province, China, with high horticultural, ornamental and medicinal value. Herein, we report a de novo genome assembly for R. henanense using a combination of PacBio long read and Illumina short read sequencing technologies. In total, we assembled 634.07 Mb with a contig N50 of 2.5 Mb, representing ~96.93% of the estimated genome size. By applying Hi-C data, 13 pseudochromosomes of R. henanense genome were assembled, covering ~98.21% of the genome assembly. The genome was composed of ~65.76% repetitive sequences and 31,098 protein-coding genes, 88.77% of which could be functionally annotated. Rhododendron henanense displayed a high level of synteny with other Rhododendron species from the Hymenanthes subgenus. Our data also suggests that R. henanense genes related to stress responses have undergone expansion, which may underly the unique abiotic and biotic stress resistance of the species. This alpine Rhododendron chromosome-scale genome assembly provides fundamental molecular resources for germplasm conservation, breeding efforts, evolutionary studies, and elucidating the unique biological characteristics of R. henanense.

qc3C: Reference-free quality control for Hi-C sequencing data

Hi-C is a sample preparation method that enables high-throughput sequencing to capture genome-wide spatial interactions between DNA molecules. The technique has been successfully applied to solve challenging problems such as 3D structural analysis of chromatin, scaffolding of large genome assemblies and more recently the accurate resolution of metagenome-assembled genomes (MAGs). Despite continued refinements, however, preparing a Hi-C library remains a complex laboratory protocol. To avoid costly failures and maximise the odds of successful outcomes, diligent quality management is recommended. Current wet-lab methods provide only a crude assay of Hi-C library quality, while key post-sequencing quality indicators used have—thus far—relied upon reference-based read-mapping. When a reference is accessible, this reliance introduces a concern for quality, where an incomplete or inexact reference skews the resulting quality indicators. We propose a new, reference-free approach that infers the total fraction of read-pairs that are a product of proximity ligation. This quantification of Hi-C library quality requires only a modest amount of sequencing data and is independent of other application-specific criteria. The algorithm builds upon the observation that proximity ligation events are likely to create k-mers that would not naturally occur in the sample. Our software tool (qc3C) is to our knowledge the first to implement a reference-free Hi-C QC tool, and also provides reference-based QC, enabling Hi-C to be more easily applied to non-model organisms and environmental samples. We characterise the accuracy of the new algorithm on simulated and real datasets and compare it to reference-based methods.

The Hi-C sequencing technique offers the potential for significant scientific insight about the spatial arrangement of DNA, however achieving such outcomes is highly dependent on the quality of the resulting sequencing library. Unlike conventional next-gen sequencing, only a fraction of a given Hi-C library contains this useful spatial information (the signal) with the remainder being effectively noise. As Hi-C remains a challenging laboratory technique, signal strength of resulting libraries can vary greatly. As a quality metric, the quantification a library’s signal content is an essential asset in any quality mitigation strategy. Quality assessment of Hi-C data has until now relied on access to a (ideally refined) reference sequence, by which indirect indicators of quality are determined. Here we describe qc3C, a software tool capable of the direct, reference-free estimation of the signal content of a Hi-C library. In doing so, not only can researchers make informed decisions on how to progress based on library information content, but eliminating the reference also enables Hi-C quality management for non-model organism and metagenomics researchers.

Chromosome-scale reference genome assembly of a diploid potato clone derived from an elite variety

Potato (Solanum tuberosum L.) is one of the most important crops with a worldwide production of 370 million metric tons. The objectives of this study were (1) to create a high-quality consensus sequence across the two haplotypes of a diploid clone derived from a tetraploid elite variety and assess the sequence divergence from the available potato genome assemblies, as well as among the two haplotypes; (2) to evaluate the new assembly’s usefulness for various genomic methods; and (3) to assess the performance of phasing in diploid and tetraploid clones, using linked-read sequencing technology. We used PacBio long reads coupled with 10x Genomics reads and proximity ligation scaffolding to create the dAg1_v1.0 reference genome sequence. With a final assembly size of 812 Mb, where 750 Mb are anchored to 12 chromosomes, our assembly is larger than other available potato reference sequences and high proportions of properly paired reads were observed for clones unrelated by pedigree to dAg1. Comparisons of the new dAg1_v1.0 sequence to other potato genome sequences point out the high divergence between the different potato varieties and illustrate the potential of using dAg1_v1.0 sequence in breeding applications.

Technical considerations in Hi-C scaffolding and evaluation of chromosome-scale genome assemblies

The recent development of ecological studies has been fueled by the introduction of massive information based on chromosome‐scale genome sequences, even for species for which genetic linkage is not accessible. This was enabled mainly by the application of Hi‐C, a method for genome‐wide chromosome conformation capture that was originally developed for investigating the long‐range interaction of chromatins. Performing genomic scaffolding using Hi‐C data is highly resource‐demanding and employs elaborate laboratory steps for sample preparation. It starts with building a primary genome sequence assembly as an input, which is followed by computation for genome scaffolding using Hi‐C data, requiring careful validation. This article presents technical considerations for obtaining optimal Hi‐C scaffolding results and provides a test case of its application to a reptile species, the Madagascar ground gecko (Paroedura picta). Among the metrics that are frequently used for evaluating scaffolding results, we investigate the validity of the completeness assessment of chromosome‐scale genome assemblies using single‐copy reference orthologues.

Homeotic Genes: Clustering, Modularity, and Diversity

Hox genes code for transcription factors and are evolutionarily conserved. They regulate a plethora of downstream targets to define the anterior-posterior (AP) body axis of a developing bilaterian embryo. Early work suggested a possible role of clustering and ordering of Hox to regulate their expression in a spatially restricted manner along the AP axis. However, the recent availability of many genome assemblies for different organisms uncovered several examples that defy this constraint. With recent advancements in genomics, the current review discusses the arrangement of Hox in various organisms. Further, we revisit their discovery and regulation in Drosophila melanogaster. We also review their regulation in different arthropods and vertebrates, with a significant focus on Hox expression in the crustacean Parahyale hawaiensis. It is noteworthy that subtle changes in the levels of Hox gene expression can contribute to the development of novel features in an organism. We, therefore, delve into the distinct regulation of these genes during primary axis formation, segment identity, and extra-embryonic roles such as in the formation of hair follicles or misregulation leading to cancer. Toward the end of each section, we emphasize the possibilities of several experiments involving various organisms, owing to the advancements in the field of genomics and CRISPR-based genome engineering. Overall, we present a holistic view of the functioning of Hox in the animal world.

De novo genome assembly of the potent medicinal plant Rehmannia glutinosa using nanopore technology

Rehmannia glutinosa is a potent medicinal plant with a significant importance in traditional Chinese medicine. Its root is enriched with various bioactive molecules mainly iridoids, possessing important pharmaceutical properties. However, the molecular biology and evolution of R. glutinosa have been largely unexplored. Here, we report a reference genome of R. glutinosa using Nanopore technology, Illumina and Hi-C sequencing. The assembly genome is 2.49 Gb long with a scaffold N50 length of 70 Mb and high heterozygosity (2%). Since R. glutinosa is an autotetraploid (4n = 56), the difference between each set of chromosomes is very small, and it is difficult to distinguish the two sets of chromosomes using Hi-C. Hence, only one set of the genome size was mounted to the chromosome level. Scaffolds covering 52.61% of the assembled genome were anchored on 14 pseudochromosomes. Over 67% of the genome consists of repetitive sequences dominated by Copia long terminal repeats and 48,475 protein-coding genes were predicted. Phylogenetic analysis corroborates the placement of R. glutinosa in the Orobanchaceae family. Our results indicated an independent and very recent whole genome duplication event that occurred 3.64 million year ago in the R. glutinosa lineage. Comparative genomics analysis demonstrated expansion of the UDP-dependent glycosyltransferases and terpene synthase gene families, known to be involved in terpenoid biosynthesis and diversification. Furthermore, the molecular biosynthetic pathway of iridoids has been clarified in this work. Collectively, the generated reference genome of R. glutinosa will facilitate discovery and development of important pharmacological compounds.

Chromosomal-scale de novo genome assemblies of Cynomolgus Macaque and Common Marmoset

Cynomolgus macaque (Macaca fascicularis) and common marmoset (Callithrix jacchus) have been widely used in human biomedical research. Long-standing primate genome assemblies used the human genome as a reference for ordering and orienting the assembled fragments into chromosomes. Here we performed de novo genome assembly of these two species without any human genome-based bias observed in the genome assemblies released earlier. We assembled PacBio long reads, and the resultant contigs were scaffolded with Hi-C data, which were further refined based on Hi-C contact maps and alternate de novo assemblies. The assemblies achieved scaffold N50 lengths of 149 Mb and 137 Mb for cynomolgus macaque and common marmoset, respectively. The high fidelity of our assembly is also ascertained by BAC-end concordance in common marmoset. Our assembly of cynomolgus macaque outperformed all the available assemblies of this species in terms of contiguity. The chromosome-scale genome assemblies produced in this study are valuable resources for non-human primate models and provide an important baseline in human biomedical research.

Regulation of mammalian 3D genome organization and histone H3K9 dimethylation by H3K9 methyltransferases

Histone H3 lysine 9 dimethylation (H3K9me2) is a highly conserved silencing epigenetic mark. Chromatin marked with H3K9me2 forms large domains in mammalian cells and overlaps well with lamina-associated domains and the B compartment defined by Hi-C. However, the role of H3K9me2 in 3-dimensional (3D) genome organization remains unclear. Here, we investigated genome-wide H3K9me2 distribution, transcriptome, and 3D genome organization in mouse embryonic stem cells following the inhibition or depletion of H3K9 methyltransferases (MTases): G9a, GLP, SETDB1, SUV39H1, and SUV39H2. We show that H3K9me2 is regulated by all five MTases; however, H3K9me2 and transcription in the A and B compartments are regulated by different MTases. H3K9me2 in the A compartments is primarily regulated by G9a/GLP and SETDB1, while H3K9me2 in the B compartments is regulated by all five MTases. Furthermore, decreased H3K9me2 correlates with changes to more active compartmental state that accompanied transcriptional activation. Thus, H3K9me2 contributes to inactive compartment setting.

Ten Years of Collaborative Progress in the Quest for Orthologs

Accurate determination of the evolutionary relationships between genes is a foundational challenge in biology. Homology-evolutionary relatedness-is in many cases readily determined based on sequence similarity analysis. By contrast, whether or not two genes directly descended from a common ancestor by a speciation event (orthologs) or duplication event (paralogs) is more challenging, yet provides critical information on the history of a gene. Since 2009, this task has been the focus of the Quest for Orthologs (QFO) Consortium. The sixth QFO meeting took place in Okazaki, Japan in conjunction with the 67th National Institute for Basic Biology conference. Here, we report recent advances, applications, and oncoming challenges that were discussed during the conference. Steady progress has been made toward standardization and scalability of new and existing tools. A feature of the conference was the presentation of a panel of accessible tools for phylogenetic profiling and several developments to bring orthology beyond the gene unit-from domains to networks. This meeting brought into light several challenges to come: leveraging orthology computations to get the most of the incoming avalanche of genomic data, integrating orthology from domain to biological network levels, building better gene models, and adapting orthology approaches to the broad evolutionary and genomic diversity recognized in different forms of life and viruses.

The chromosome-level genome of dragon fruit reveals whole-genome duplication and chromosomal co-localization of betacyanin biosynthetic genes

Dragon fruits are tropical fruits economically important for agricultural industries. As members of the family of Cactaceae, they have evolved to adapt to the arid environment. Here we report the draft genome of Hylocereus undatus, commercially known as the white-fleshed dragon fruit. The chromosomal level genome assembly contains 11 longest scaffolds corresponding to the 11 chromosomes of H. undatus. Genome annotation of H. undatus found ~29,000 protein-coding genes, similar to Carnegiea gigantea (saguaro). Whole-genome duplication (WGD) analysis revealed a WGD event in the last common ancestor of Cactaceae followed by extensive genome rearrangements. The divergence time between H. undatus and C. gigantea was estimated to be 9.18 MYA. Functional enrichment analysis of orthologous gene clusters (OGCs) in six Cactaceae plants found significantly enriched OGCs in drought resistance. Fruit flavor-related functions were overrepresented in OGCs that are significantly expanded in H. undatus. The H. undatus draft genome also enabled the discovery of carbohydrate and plant cell wall-related functional enrichment in dragon fruits treated with trypsin for a longer storage time. Lastly, genes of the betacyanin (a red-violet pigment and antioxidant with a very high concentration in dragon fruits) biosynthetic pathway were found to be co-localized on a 12 Mb region of one chromosome. The consequence may be a higher efficiency of betacyanin biosynthesis, which will need experimental validation in the future. The H. undatus draft genome will be a great resource to study various cactus plants.

Erythrocytes 3D genome organization in vertebrates

Generation of mature red blood cells, consisting mainly of hemoglobin, is a remarkable example of coordinated action of various signaling networks. Chromatin condensation is an essential step for terminal erythroid differentiation and subsequent nuclear expulsion in mammals. Here, we profiled 3D genome organization in the blood cells from ten species belonging to different vertebrate classes. Our analysis of contact maps revealed a striking absence of such 3D interaction patterns as loops or TADs in blood cells of all analyzed representatives. We also detect large-scale chromatin rearrangements in blood cells from mammals, birds, reptiles and amphibians: their contact maps display strong second diagonal pattern, representing an increased frequency of long-range contacts, unrelated to TADs or compartments. This pattern is completely atypical for interphase chromosome structure. We confirm that these principles of genome organization are conservative in vertebrate erythroid cells.

Three-dimensional chromatin in infectious disease-A role for gene regulation and pathogenicity?

The recent Coronavirus Disease 2019 pandemic has once again reminded us the importance of understanding infectious diseases. One important but understudied area in infectious disease research is the role of nuclear architecture or the physical arrangement of the genome in the nucleus in controlling gene regulation and pathogenicity. Recent advances in research methods, such as Genome-wide chromosome conformation capture using high-throughput sequencing (Hi-C), have allowed for easier analysis of nuclear architecture and chromosomal reorganization in both the infectious disease agents themselves as well as in their host cells. This review will discuss broadly on what is known about nuclear architecture in infectious disease, with an emphasis on chromosomal reorganization, and briefly discuss what steps are required next in the field.

In this review, we examine the current state of nuclear architecture in infectious diseases with an emphasis on chromosomal reorganization. Nuclear architecture plays an important role in regulation of transcription for several pathogens, as well as inflammatory responses in their host. Recent advances in technologies such as Hi-C have allowed in-depth studies of chromosomal reorganization during infectious disease development and provided insights into transcription mechanisms and pathogenicity. In addition, it has been demonstrated that pathogens can also affect/utilize the hosts nuclear architecture. These areas are heavily understudied in pathogens, and we hope this review will provide a comprehensive review on the current state of nuclear architecture in infectious diseases and provide an additional avenue for eradication efforts.

Cell culture-based karyotyping of orectolobiform sharks for chromosome-scale genome analysis

Karyotyping, traditionally performed using cytogenetic techniques, is indispensable for validating genome assemblies whose sequence lengths can be scaled up to chromosome sizes using modern methods. Karyotype reports of chondrichthyans are scarce because of the difficulty in cell culture. Here, we focused on carpet shark species and the culture conditions for fibroblasts and lymphocytes. The utility of the cultured cells enabled the high-fidelity characterization of their karyotypes, namely 2n = 102 for the whale shark (Rhincodon typus) and zebra shark (Stegostoma fasciatum), and 2n = 106 for the brownbanded bamboo shark (Chiloscyllium punctatum) and whitespotted bamboo shark (C. plagiosum). We identified heteromorphic XX/XY sex chromosomes for the two latter species and demonstrated the first-ever fluorescence in situ hybridization of shark chromosomes prepared from cultured cells. Our protocols are applicable to diverse chondrichthyan species and will deepen the understanding of early vertebrate evolution at the molecular level.

Multifaceted Hi-C benchmarking: what makes a difference in chromosome-scale genome scaffolding?

Background

Hi-C is derived from chromosome conformation capture (3C) and targets chromatin contacts on a genomic scale. This method has also been used frequently in scaffolding nucleotide sequences obtained by de novo genome sequencing and assembly, in which the number of resultant sequences rarely converges to the chromosome number. Despite its prevalent use, the sample preparation methods for Hi-C have not been intensively discussed, especially from the standpoint of genome scaffolding.

Results

To gain insight into the best practice of Hi-C scaffolding, we performed a multifaceted methodological comparison using vertebrate samples and optimized various factors during sample preparation, sequencing, and computation. As a result, we identified several key factors that helped improve Hi-C scaffolding, including the choice and preparation of tissues, library preparation conditions, the choice of restriction enzyme(s), and the choice of scaffolding program and its usage.

Conclusions

This study provides the first comparison of multiple sample preparation kits/protocols and computational programs for Hi-C scaffolding by an academic third party. We introduce a customized protocol designated “inexpensive and controllable Hi-C (iconHi-C) protocol,” which incorporates the optimal conditions identified in this study, and demonstrate this technique on chromosome-scale genome sequences of the Chinese softshell turtle Pelodiscus sinensis.

Early vertebrate origin of CTCFL, a CTCF paralog, revealed by proximity-guided shark genome scaffolding

The nuclear protein CCCTC-binding factor (CTCF) contributes as an insulator to chromatin organization in diverse animals. The gene encoding this protein has a paralog which was first identified to be expressed exclusively in the testis in mammals and designated as CTCFL (also called BORIS). CTCFL orthologs were reported only among amniotes, and thus CTCFL was once thought to have arisen in the amniote lineage. In this study, we identified elasmobranch CTCFL orthologs, and investigated its origin with the aid of a shark genome assembly improved by proximity-guided scaffolding. Our analysis employing evolutionary interpretation of syntenic gene location suggested an earlier timing of the gene duplication between CTCF and CTCFL than previously thought, that is, around the common ancestor of extant vertebrates. Also, our transcriptomic sequencing revealed a biased expression of the catshark CTCFL in the testis, suggesting the origin of the tissue-specific localization in mammals more than 400 million years ago. To understand the historical process of the functional consolidation of the long-standing chromatin regulator CTCF, its additional paralogs remaining in some of the descendant lineages for spatially restricted transcript distribution should be taken into consideration.

Apomixis Technology: Separating the Wheat from the Chaff

Projections indicate that current plant breeding approaches will be unable to incorporate the global crop yields needed to deliver global food security. Apomixis is a disruptive innovation by which a plant produces clonal seeds capturing heterosis and gene combinations of elite phenotypes. Introducing apomixis into hybrid cultivars is a game-changing development in the current plant breeding paradigm that will accelerate the generation of high-yield cultivars. However, apomixis is a developmentally complex and genetically multifaceted trait. The central problem behind current constraints to apomixis breeding is that the genomic configuration and molecular mechanism that initiate apomixis and guide the formation of a clonal seed are still unknown. Today, not a single explanation about the origin of apomixis offer full empirical coverage, and synthesizing apomixis by manipulating individual genes has failed or produced little success. Overall evidence suggests apomixis arise from a still unknown single event molecular mechanism with multigenic effects. Disentangling the genomic basis and complex genetics behind the emergence of apomixis in plants will require the use of novel experimental approaches benefiting from Next Generation Sequencing technologies and targeting not only reproductive genes, but also the epigenetic and genomic configurations associated with reproductive phenotypes in homoploid sexual and apomictic carriers. A comprehensive picture of most regulatory changes guiding apomixis emergence will be central for successfully installing apomixis into the target species by exploiting genetic modification techniques.