Hi-C Sequencing Research Articles

Chromosome-level genome assembly of the northern Pacific seastar Asterias amurensis

Asterias amurensis has attracted widespread concern because of its population outbreaks, which has impacted fisheries and aquaculture, as well as disrupting local ecosystems. A high-quality reference genome is necessary to better investigate mechanisms of outbreak and adaptive changes. Combining PacBio HiFi and Hi-C sequencing data, we generated a chromosome-level A. amurensis genome with a size of 491.53 Mb. The contig N50 and scaffold N50 were 8.05 and 23.75 Mb, respectively. The result of BUSCO analysis revealed a completeness score of 98.85%. A total of 16,531 protein-coding genes were predicted in the genome, of which 94.63% were functionally annotated. The high-quality genome assembly resulting from this study will provide a valuable genetic resource for future research on the mechanism of population outbreaks and invasion ecology.

Chromosome-level genome assembly and annotation of Zicaitai (Brassica rapa var. purpuraria)

Zicaitai is a seasonal vegetable known for its high anthocyanin content in both stalks and leaves, yet its reference genome has not been published to date. Here, we generated the first chromosome-level genome assembly of Zicaitai using a combination of PacBio long-reads, Illumina short-reads, and Hi-C sequencing techniques. The final genome length is 474.12 Mb with a scaffold N50 length of 43.82 Mb, a BUSCO score of 99.30% and the LAI score of 10.14. Repetitive elements accounted for 60.89% (288.72 Mb) of the genome, and Hi-C data enabled the allocation of 430.87 Mb of genome sequences to ten pseudochromosomes. A total of 42,051 protein-coding genes were successfully predicted using multiple methods, of which 99.74% were functionally annotated. Notably, comparing the genome of Zicaitai with seven other species in the Cruciferae family revealed strong conservation in terms of gene numbers and structures. Overall, the high-quality genome assembly provides a critical resource for studying the genetic basis of important agronomic traits in Zicaitai.

Debugging and consolidating multiple synthetic chromosomes reveals combinatorial genetic interactions

The Sc2.0 project is building a eukaryotic synthetic genome from scratch. A major milestone has been achieved with all individual Sc2.0 chromosomes assembled. Here, we describe the consolidation of multiple synthetic chromosomes using advanced endoreduplication intercrossing with tRNA expression cassettes to generate a strain with 6.5 synthetic chromosomes. The 3D chromosome organization and transcript isoform profiles were evaluated using Hi-C and long-read direct RNA sequencing. We developed CRISPR Directed Biallelic URA3-assisted Genome Scan, or "CRISPR D-BUGS," to map phenotypic variants caused by specific designer modifications, known as "bugs." We first fine-mapped a bug in synthetic chromosome II (synII) and then discovered a combinatorial interaction associated with synIII and synX, revealing an unexpected genetic interaction that links transcriptional regulation, inositol metabolism, and tRNASerCGA abundance. Finally, to expedite consolidation, we employed chromosome substitution to incorporate the largest chromosome (synIV), thereby consolidating >50% of the Sc2.0 genome in one strain.

Chromosome-level genome assembly of sea scallop Placopecten magellanicus provides insights into the genetic characteristics and adaptive evolution of large scallops

Placopecten magellanicus (Gmelin, 1791), a deep-sea Atlantic scallop, holds significant commercial value as a benthic marine bivalve along the northwest Atlantic coast. Recognizing its economic importance, the need to reconstruct its genome assembly becomes apparent, fostering insights into natural resources and generic breeding potential. This study reports a high-quality chromosome-level genome of P. magellanicus, achieved through the integration of Illumina short read sequencing, PacBio HiFi sequencing, and Hi-C sequencing techniques. The resulting assembly spans 1778 Mb with a scaffold N50 of 86.71 Mb. An intriguing observation arises – the genome size of P. magellanicus surpasses that of its Pectinidae family peers by 1.80 to 2.46 times. Within this genome, 28,111 protein-coding genes were identified. Comparative genomic analysis involving five scallop species unveils the critical determinant of this expanded genome: the proliferation of repetitive sequences recently inserted, contributing to its enlarged size. The landscape of whole genome collinearity sheds light on the relationships among scallop species, enhancing our broader understanding of their genomic framework. This genome provides genomic resources for future molecular biology research on scallops and serves as a guide for the exploration of longevity-related genes in scallops.

50. Analysis of gene rearrangements in neoplasms with Hi-C sequencing using fresh-frozen and FFPE specimens

50. Analysis of gene rearrangements in neoplasms with Hi-C sequencing using fresh-frozen and FFPE specimens

Arbuscular mycorrhizal fungi heterokaryons have two nuclear populations with distinct roles in host-plant interactions.

Arbuscular mycorrhizal fungi (AMF) are prominent root symbionts that can carry thousands of nuclei deriving from two parental strains in a large syncytium. These co-existing genomes can also vary in abundance with changing environmental conditions. Here we assemble the nuclear genomes of all four publicly available AMF heterokaryons using PacBio high-fidelity and Hi-C sequencing. We find that the two co-existing genomes of these strains are phylogenetically related but differ in structure, content and epigenetics. We confirm that AMF heterokaryon genomes vary in relative abundance across conditions and show this can lead to nucleus-specific differences in expression during interactions with plants. Population analyses also reveal signatures of genetic exchange indicative of past events of sexual reproduction in these strains. This work uncovers the origin and contribution of two nuclear genomes in AMF heterokaryons and opens avenues for the improvement and environmental application of these strains.

The first chromosome-level genome assembly and transcriptome sequencing provide insights into cantharidin production of the blister beetles.

Blister beetles (Coleoptera: Meloidae) produce a natural defensive toxin cantharidin (CTD), which has been used for various cancer treatments and other diseases. Currently, the lack of chromosome-level reference genomes in Meloidae limits further understanding of the mechanism of CTD biosynthesis and environmental adaptation. In this study, the chromosome-level genome assembly of Mylabris phalerata was generated based on PacBio and Hi-C sequencing. This reference genome was about 136.68 Mb in size with contig N50 of 9.17 Mb and composed of 12 chromosomes. In comparison to six other Coleoptera insects, M. phalerata exhibited multiple expanded gene families enriched in juvenile hormone (JH) biosynthetic process pathway, farnesol dehydrogenase activity, and cytochrome P450, which may be related to CTD biosynthesis. Consistently, the transcriptomic analysis suggested the "terpenoid backbone biosynthesis" pathway and "the juvenile hormone" as putative core pathways of CTD biosynthesis and presented eight up-regulated differential expression genes in male adults as candidate genes. It is possible that the restricted feeding niche and lifestyle of M. phalerata were the cause of the gene family's contraction of odorant binding proteins. The ABC transporters (ABCs) related to exporting bound toxins out of the cell and the resistance to the self-secreted toxins (e.g. CTD) were also contracted, possibly due to other self-protection strategies in M. phalerata. A foundation of understanding CTD biosynthesis and environmental adaptation of blister beetles will be established by our reference genome and discoveries.

Whole-genome and genome-wide association studies improve key agricultural traits of safflower for industrial and medicinal use.

Safflower (Carthamus tinctorius) is widely cultivated around the world for its seeds and flowers. The presence of linoleic acid (LA) in its seeds and hydroxysafflor yellow A (HSYA) in its flowers are the crucial traits that enable safflower to be used for industrial and medicinal purposes. Understanding the genetic control of these traits is essential for optimizing the quality of safflower and its breeding. To further this research, we present a chromosome-scale assembly of the genome of the safflower variety 'Chuanhonghua 1', which was achieved using an integrated strategy combining Illumina, Oxford Nanopore, and Hi-C sequencing. We obtained a 1.17-Gb assembly with a contig N50 of 1.08Mb, and all assembled sequences were assigned to 12 pseudochromosomes. Safflower's evolution involved the core eudicot γ-triplication event and a whole-genome duplication event, which led to large-scale genomic rearrangements. Extensive genomic shuffling has occurred since the divergence of the ancestor of dicotyledons. We conducted metabolite and transcriptome profiles with time- and part-dependent changes and screened candidate genes that significantly contribute to seed lipid biosynthesis. We also analyzed key gene families that participate in LA and HSYA biosynthesis. Additionally, we re-sequenced 220 safflower lines and carried out a genome-wide association study using high-quality SNP data for eight agronomic traits. We identified SNPs related to important traits in safflower. Besides, the candidate gene HH_034464 (CtCGT1) was shown to be involved in the biosynthesis of HSYA. Overall, we provide a high-quality reference genome and elucidate the genetic basis of LA and HSYA biosynthesis in safflower. This vast amount of data will benefit further research for functional gene mining and breeding in safflower.

Genetic and Molecular Evidence of a Tetrapolar Mating System in the Edible Mushroom Grifola frondosa.

Grifola frondosa is a valuable edible fungus with high nutritional and medicinal values. The mating systems of fungi not only offer practical strategies for breeding, but also have far-reaching effects on genetic variability. Grifola frondosa has been considered as a sexual species with a tetrapolar mating system based on little experimental data. In the present study, one group of test crosses and six groups of three-round mating experiments from two parental strains were conducted to determine the mating system in G. frondosa. A chi-squared test of the results of the test-cross mating experiments indicated that they satisfied Mendelian segregation, while a series of three-round mating experiments showed that Mendelian segregation was not satisfied, implying a segregation distortion phenomenon in G. frondosa. A genomic map of the G. frondosa strain, y59, grown from an LMCZ basidiospore, with 40.54 Mb and 12 chromosomes, was generated using genome, transcriptome and Hi-C sequencing technology. Based on the genomic annotation of G. frondosa, the mating-type loci A and B were located on chromosomes 1 and 11, respectively. The mating-type locus A coded for the β-fg protein, HD1, HD2 and MIP, in that order. The mating-type locus B consisted of six pheromone receptors (PRs) and five pheromone precursors (PPs) in a crossed order. Moreover, both HD and PR loci may have only one sublocus that determines the mating type in G. frondosa. The nonsynonymous SNP and indel mutations between the A1B1 and A2B2 mating-type strains and the reference genome of y59 only occurred on genes HD2 and PR1/2, preliminarily confirming that the mating type of the y59 strain was A1B2 and not A1B1. Based on the genetic evidence and the more reliable molecular evidence, the results reveal that the mating system of G. frondosa is tetrapolar. This study has important implications for the genetics and hybrid breeding of G. frondosa.

Integrating chromatin conformation information in a self-supervised learning model improves metagenome binning.

Metagenome binning is a key step, downstream of metagenome assembly, to group scaffolds by their genome of origin. Although accurate binning has been achieved on datasets containing multiple samples from the same community, the completeness of binning is often low in datasets with a small number of samples due to a lack of robust species co-abundance information. In this study, we exploited the chromatin conformation information obtained from Hi-C sequencing and developed a new reference-independent algorithm, Metagenome Binning with Abundance and Tetra-nucleotide frequencies-Long Range (metaBAT-LR), to improve the binning completeness of these datasets. This self-supervised algorithm builds a model from a set of high-quality genome bins to predict scaffold pairs that are likely to be derived from the same genome. Then, it applies these predictions to merge incomplete genome bins, as well as recruit unbinned scaffolds. We validated metaBAT-LR's ability to bin-merge and recruit scaffolds on both synthetic and real-world metagenome datasets of varying complexity. Benchmarking against similar software tools suggests that metaBAT-LR uncovers unique bins that were missed by all other methods. MetaBAT-LR is open-source and is available at https://bitbucket.org/project-metabat/metabat-lr.

Topological reorganization and functional alteration of distinct genomic components in gallbladder cancer.

Altered three-dimensional architecture of chromatin influences various genomic regulators and subsequent gene expression in human cancer. However, knowledge of the topological rearrangement of genomic hierarchical layers in cancer is largely limited. Here, by taking advantage of in situ Hi-C, RNA-sequencing, and chromatin immunoprecipitation sequencing (ChIP-seq), we investigated structural reorganization and functional changes in chromosomal compartments, topologically associated domains (TADs), and CCCTC binding factor (CTCF)-mediated loops in gallbladder cancer (GBC) tissues and cell lines. We observed that the chromosomal compartment A/B switch was correlated with CTCF binding levels and gene expression changes. Increased inter-TAD interactions with weaker TAD boundaries were identified in cancer cell lines relative to normal controls. Furthermore, the chromatin short loops and cancer unique loops associated with chromatin remodeling and epithelial-mesenchymal transition activation were enriched in cancer compared with their control counterparts. Cancer-specific enhancer-promoter loops, which contain multiple transcription factor binding motifs, acted as a central element to regulate aberrant gene expression. Depletion of individual enhancers in each loop anchor that connects with promoters led to the inhibition of their corresponding gene expressions. Collectively, our data offer the landscape of hierarchical layers of cancer genome and functional alterations that contribute to the development of GBC.

Chromosome-level Dinobdella ferox genome provided a molecular model for its specific parasitism

BackgroundDinobdella ferox is the most frequently reported leech species parasitizing the mammalian nasal cavity. However, the molecular mechanism of this special parasitic behavior has remained largely unknown.MethodsPacBio long-read sequencing, next-generation sequencing (NGS), and Hi-C sequencing were employed in this study to generate a novel genome of D. ferox, which was annotated with strong certainty using bioinformatics methods. The phylogenetic and genomic alterations of D. ferox were then studied extensively alongside the genomes of other closely related species. The obligatory parasitism mechanism of D. ferox was investigated using RNA-seq and proteomics data.ResultsPacBio long-read sequencing and NGS yielded an assembly of 228 Mb and contig N50 of 2.16 Mb. Along Hi-C sequencing, 96% of the sequences were anchored to nine linkage groups and a high-quality chromosome-level genome was generated. The completed genome included 19,242 protein-coding genes. For elucidating the molecular mechanism of nasal parasitism, transcriptome data were acquired from the digestive tract and front/rear ends of D. ferox. Examining secretory proteins in D. ferox saliva helped to identify intimate connections between these proteins and membrane proteins in nasal epithelial cells. These interacting proteins played important roles in extracellular matrix (ECM)–receptor interaction, tight junction, focal adhesion, and adherens junction. The interaction between D. ferox and mammalian nasal epithelial cells included three major steps of pattern recognition, mucin connection and breakdown, and repair of ECM. The remodeling of ECM between epithelial cells of the nasal mucosa and epithelial cells of D. ferox may produce a stable adhesion environment for parasitism.ConclusionsOur study represents the first-ever attempt to propose a molecular model for specific parasitism. This molecular model may serve as a practical reference for parasitism models of other species and a theoretical foundation for a molecular process of parasitism.Graphical abstract

High-quality wild barley genome assemblies and annotation with Nanopore long reads and Hi-C sequencing data

Wild barley, from “Evolution Canyon (EC)” in Mount Carmel, Israel, are ideal models for cereal chromosome evolution studies. Here, the wild barley EC_S1 is from the south slope with higher daily temperatures and drought, while EC_N1 is from the north slope with a cooler climate and higher relative humidity, which results in a differentiated selection due to contrasting environments. We assembled a 5.03 Gb genome with contig N50 of 3.53 Mb for wild barley EC_S1 and a 5.05 Gb genome with contig N50 of 3.45 Mb for EC_N1 using 145 Gb and 160.0 Gb Illumina sequencing data, 295.6 Gb and 285.35 Gb Nanopore sequencing data and 555.1 Gb and 514.5 Gb Hi-C sequencing data, respectively. BUSCOs and CEGMA evaluation suggested highly complete assemblies. Using full-length transcriptome data, we predicted 39,179 and 38,373 high-confidence genes in EC_S1 and EC_N1, in which 93.6% and 95.2% were functionally annotated, respectively. We annotated repetitive elements and non-coding RNAs. These two wild barley genome assemblies will provide a rich gene pool for domesticated barley.

The completed genome sequence of Pestalotiopsis versicolor, a pathogenic ascomycete fungus with implications for bayberry production

The pathogenic fungus Pestalotiopsis versicolor is a major etiological agent of fungal twig blight disease affecting bayberry trees. However, the lack of complete genome sequence information for this crucial pathogenic fungus hinders the molecular and genetic investigation of its pathogenic mechanism. To address this knowledge gap, we have generated the complete genome sequence of P. versicolor strain XJ27, employing a combination of Illumina, PacBio, and Hi-C sequencing technologies. This comprehensive genome sequence, comprising 7 chromosomes with an N50 contig size of 7,275,017 bp, a GC content ratio of 50.16%, and a total size of 50.80 Mb, encompasses 13,971 predicted coding genes. By performing comparative genomic analysis between P. versicolor and the genomes of eleven plant-pathogenic fungi, as well as three closely related fungi within the same group, we have gained initial insights into its evolutionary trajectory, particularly through gene family analysis. These findings shed light on the distinctive characteristics and evolutionary history of P. versicolor. Importantly, the availability of this high-quality genetic resource will serve as a foundational tool for investigating the biology, molecular pathogenesis, and virulence of P. versicolor. Furthermore, it will facilitate the development of more potent antifungal medications by uncovering potential vulnerabilities in its genetic makeup.

A chromosome-level genome assembly for the astaxanthin-producing microalga Haematococcus pluvialis

The green microalga Haematococcus pluvialis can synthesize high amounts of astaxanthin, which is a valuable antioxidant that has been utilized in human health, cosmetics, and aquaculture. To illustrate detailed molecular clues to astaxanthin yield, we performed PacBio HIFI along with Hi-C sequencing to construct an improved chromosome-level haplotypic genome assembly with 32 chromosomes and a genome size of 316.0 Mb. Its scaffold N50 (942.6 kb) and contig N50 (304.8 kb) have been upgraded remarkably from our previous genome draft, and a total of 32,416 protein-coding genes were predicted. We also established a high-evidence phylogenetic tree from seven representative algae species, with the main aim to calculate their divergence times and identify expanded/contracted gene families. We also characterized genome-wide localizations on chromosomes of some important genes such as five BKTs (encoding beta-carotene ketolases) that are putatively involved in astaxanthin production. In summary, we reported the first chromosome-scale map of H. pluvialis, which provides a valuable genetic resource for in-depth biomedical investigations on this momentous green alga and commercial astaxanthin bioproduction.

Genome assembly of two diploid and one auto-tetraploid Cyclocarya paliurus genomes

Cyclocarya paliurus, an endemic species in the genus Juglandaceae with the character of heterodichogamy, is one of triterpene-rich medicinal plants in China. To uncover the genetic mechanisms behind the special characteristics, we sequenced the genomes of two diploid (protandry, PA-dip and protogyny, PG-dip) and one auto-tetraploid (PA-tetra) C. paliurus genomes. Based on 134.9 (~225x), 75.5 (~125x) and 271.8 Gb (~226x) subreads of PacBio platform sequencing data, we assembled 586.62 Mb (contig N50 = 1.9 Mb), 583.45 Mb (contig N50 = 1.4 Mb), and 2.38 Gb (contig N50 = 430.9 kb) for PA-dip, PG-dip and PA-tetra genome, respectively. Furthermore, 543.53, 553.87, and 2168.65 Mb in PA-dip, PG-dip, and PA-tetra, were respectively anchored to 16, 16, and 64 pseudo-chromosomes using over 65.4 Gb (~109x), 68 Gb (~113x), and 264 (~220x) Hi-C sequencing data. Annotation of PA-dip, PG-dip, and PA-tetra genome assembly identified 34,699, 35,221, and 34,633 protein-coding genes (90,752 gene models) or allele-defined genes, respectively. In addition, 45 accessions from nine locations were re-sequenced, and more than 10 × coverage reads were generated.

A chromosome-level genome assembly of a deep-sea starfish (Zoroaster cf. ophiactis)

Understanding of adaptation and evolution of organisms in the deep sea requires more genomic resources. Zoroaster cf. ophiactis is a sea star in the family Zoroasteridae occurring exclusively in the deep sea. In this study, a chromosome-level genome assembly for Z. cf. ophiactis was generated by combining Nanopore long-read, Illumina short-read, and Hi-C sequencing data. The final assembly was 1,002.0 Mb in length, with a contig N50 of 376 Kb and a scaffold N50 of 40.4 Mb, and included 22 pseudo-chromosomes, covering 92.3% of the assembly. Completeness analysis evaluated with BUSCO revealed that 95.91% of the metazoan conserved genes were complete. Additionally, 39,426 protein-coding genes were annotated for this assembly. This chromosome-level genome assembly represents the first high-quality genome for the deep-sea Asteroidea, and will provide a valuable resource for future studies on evolution and adaptation of deep-sea echinoderms.

High-Quality Genome Assembly and Genome-Wide Association Study of Male Sterility Provide Resources for Flax Improvement.

Flax is an economic crop with a long history. It is grown worldwide and is mainly used for edible oil, industry, and textiles. Here, we reported a high-quality genome assembly for "Neiya No. 9", a popular variety widely grown in China. Combining PacBio long reads, Hi-C sequencing, and a genetic map reported previously, a genome assembly of 473.55 Mb was constructed, which covers ~94.7% of the flax genome. These sequences were anchored onto 15 chromosomes. The N50 lengths of the contig and scaffold were 0.91 Mb and 31.72 Mb, respectively. A total of 32,786 protein-coding genes were annotated, and 95.9% of complete BUSCOs were found. Through morphological and cytological observation, the male sterility of flax was considered dominant nuclear sterility. Through GWAS analysis, the gene LUSG00017705 (cysteine synthase gene) was found to be closest to the most significant SNP, and the expression level of this gene was significantly lower in male sterile plants than in fertile plants. Among the significant SNPs identified in the GWAS analysis, only two were located in the coding region, and these two SNPs caused changes in the protein encoded by LUSG00017565 (cysteine protease gene). It was speculated that these two genes may be related to male sterility in flax. This is the first time the molecular mechanism of male sterility in flax has been reported. The high-quality genome assembly and the male sterility genes revealed, provided a solid foundation for flax breeding.

A genome assembly for the southern Pacific rattlesnake, Crotalus oreganus helleri, in the western rattlesnake species complex.

Rattlesnakes play important roles in their ecosystems by regulating prey populations, are involved in complex coevolutionary dynamics with their prey, and exhibit a variety of unusual adaptations, including maternal care, heat-sensing pit organs, hinged fangs, and medically-significant venoms. The western rattlesnake (Crotalus oreganus) is one of the widest ranging rattlesnake species, with a distribution from British Columbia, where it is listed as threatened, to Baja California and east across the Great Basin to western Wyoming, Colorado and New Mexico. Here, we report a new reference genome assembly for one of six currently recognized subspecies, C. oreganus helleri, as part of the California Conservation Genomics Project (CCGP). Consistent with the reference genomic sequencing strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology to produce a de novo assembled genome. The assembly comprises a total of 698 scaffolds spanning 1,564,812,557 base pairs, has a contig N50 of 64.7 Mb, a scaffold N50 of 110.8 Mb, and BUSCO complete score of 90.5%. This reference genome will be valuable for studies on the genomic basis of venom evolution and variation within Crotalus, in resolving the taxonomy of C. oreganus and its relatives, and for the conservation and management of rattlesnakes in general.

Chromosome-level genome assembly of watershield (Brasenia schreberi)

Watershield (Brasenia schreberi) is an aquatic plant that belongs to the basal angiosperm family Cabombaceae. This species has been cultivated as an aquatic vegetable for more than 3000 years in East Asia, but the natural populations have greatly declined in recent decades and have become endangered in several countries of East Asia. In this study, by using PacBio long reads, Illumina short reads, and Hi-C sequencing data, we assembled the genome of B. schreberi, which was approximately 1170.4 Mb in size with a contig N50 of 7.1 Mb. Of the total assembled sequences, 93.6% were anchored to 36 pseudochromosomes with a scaffold N50 of 28.9 Mb. A total of 74,699 protein-coding genes were predicted in the B. schreberi genome, and 558 Mb of repetitive elements occupying 47.69% of the genome were identified. BUSCO analysis yielded a completeness score of 95.8%. The assembled high-quality genome of B. schreberi will be a valuable reference for the study of conservation, evolution and molecular breeding in this species.