Hi-C Sequencing Research Articles

Chromosome-level genome assemblies of vulnerable male and female elongate loach (Leptobotia elongata)

Endemic to the upper and middle reaches of the Yangtze River in China, elongate loach (Leptobotia elongata) has become a vulnerable species mainly due to overfishing and habitat destruction. Thus far, no genome data of this species are reported. As a result, lacking of such genomic information has restricted practical conservation and utilization of this economic fish. Here, we constructed chromosome-level genome assemblies for both male and female elongate loach by integration of MGI, PacBio HiFi and Hi-C sequencing technologies. Two primary genome assemblies (586-Mb and 589-Mb) were obtained for female and male fishes, respectively. Indeed, 98.22% and 98.61% of the contig sequences were anchored onto 25 chromosomes, with identification of 26.22% and 25.92% repeat contents in both assembled genomes. Meanwhile, a total of 25,215 and 25,253 protein-coding genes were annotated, of which 97.41% and 98.8% could be predicted with functions. Taken together, our genome data presented here provide a valuable genomic resource for in-depth evolutionary and functional research, as well as molecular breeding and conservation of this economic fish species.

Evaluation of Hi-C Sequencing for Detection of Gene Fusions in Hematologic and Solid Tumor Pediatric Cancer Samples.

Hi-C sequencing is a DNA-based next-generation sequencing method that preserves the 3D genome conformation and has shown promise in detecting genomic rearrangements in translational research studies. To evaluate Hi-C as a potential clinical diagnostic platform, analytical concordance with routine laboratory testing was assessed using primary pediatric leukemia and sarcoma specimens. Archived viable and non-viable frozen leukemic cells and formalin-fixed paraffin-embedded (FFPE) tumor specimens were analyzed. Pediatric acute myeloid leukemia (AML) and alveolar rhabdomyosarcoma (A-RMS) specimens with known genomic rearrangements were subjected to Hi-C to assess analytical concordance. Subsequently, a discovery cohort consisting of AML and acute lymphoblastic leukemia (ALL) cases without known genomic rearrangements based on prior clinical diagnostic testing was evaluated to determine whether Hi-C could detect rearrangements. Using a standard sequencing depth of 50 million raw read-pairs per sample, or approximately 5X raw genomic coverage, we observed 100% concordance between Hi-C and previous clinical cytogenetic and molecular testing. In the discovery cohort, a clinically relevant gene fusion was detected in 45% of leukemia cases (5/11). This study provides an institutional proof of principle evaluation of Hi-C sequencing to medical diagnostic testing as it identified several clinically relevant rearrangements, including those that were missed by current clinical testing workflows.

A high-quality chromosome-level genome assembly of the topmouth culter (Culter alburnus Basilewsky, 1855)

Culter alburnus is extensively distributed in various rivers and lakes across China. As a widely adaptive fish species, it has significant economic values and special ecological roles. To meet research demands and provide better genomic resources, in this research, a chromosome-level genome assembly was constructed using HiFi long-reads and Hi-C sequencing data. Compared with the published versions, our genome assembly is of higher quality with only 31 gaps and closer to its true structure and sequence. The genome size was 1.052 Gb, with a contig N50 of 32.92 Mb and a scaffold N50 of 43.09 Mb. 55 contigs were anchored to 24 chromosomes on the basis of Hi-C data. A total of 598.23 Mb of repetitive sequences were annotated and 28,228 protein-coding genes were predicted. Additionally, BUSCO assessment indicated assembly and annotation scores of 98.3% and 99.2%, respectively. This high-quality genome will provide scientific support for excavating the species characteristics of C. alburnus and exploring its molecular mechanisms in response to environmental changes and stress.

Chromosomal-level genome assembly of Hylurgus ligniperda: insights into host adaptation and environmental tolerance

BackgroundHylurgus ligniperda (Coleoptera: Curculionidae) is a worldwide forest quarantine pest. It is widely distributed, has many host tree species, and possesses strong adaptability. To explore its environmental adaptability and the related molecular mechanisms, we conducted chromosome-level genome sequencing and analyzed the transcriptome under different environmental factors, identifying key expressed genes.ResultsWe employed PacBio, Illumina, and Hi-C sequencing techniques to assemble a 520 Mb chromosomal-level genome of H. ligniperda, obtaining an N50 of 39.97 Mb across 138 scaffolds. A total of 10,765 protein-coding genes were annotated after repeat masking. Fourteen chromosomes were identified, among which Hyli14 was determined to be the sex chromosome. Survival statistics were tested over various growth periods under high temperature and low humidity conditions. The maximum survival period of adults reached 292 days at 25 °C, 65% relative humidity. In comparison, the maximum survival period was 14 days under 35 °C, 65% relative humidity, and 106 days under 25°C, 40% relative humidity. This indicated that environmental stress conditions significantly reduced adults' survival period. We further conducted transcriptome analysis to screen for potentially influential differentially expressed genes, such as CYP450 and Histone. Subsequently, we performed gene family analysis to gain insights into their functions and interactions, such as CYP450 and Histone. CYP450 genes affected the detoxification metabolism of enzymes in the Cytochrome P450 pathway to adapt to different environments. Histone genes are involved in insect hormone biosynthesis and longevity-regulating pathways in H. ligniperda to adapt to environmental stress.ConclusionsThe genome at the chromosome level of H. ligniperda was assembled for the first time. The mortality of H. ligniperda increased significantly at 35 ℃, 65% RH, and 25 ℃, 40% RH. CYP450 and Histone genes played an important role in response to environmental stress. This genome offers a substantial genetic resource for investigating the molecular mechanisms behind beetle invasion and spread.

Chromosome-level genome assembly of the invasive pest Pseudococcus jackbeardsleyi (Hemiptera: Pseudococcidae)

Among over 2,000 species of mealybugs (Hemiptera: Pseudococcidae), only 13 genomes have been published so far, seriously limiting the researches on the phylogeny and adaptive evolution of this group. The continuous publication of mealybug genomes will significantly facilitate our exploration of the biological characteristics, detrimental attributes, and control strategies of the Pseudococcidae family. Jack Beardsley mealybug (Pseudococcus jackbeardsleyi) as one of the hazardous invasive pests, it could cause enormous losses to the fruit and vegetable industries worldwide. Herein, we combined Nanopore long-read, short-read Illumina and Hi-C sequencing, generating a high-quality chromosome-level genome assembly of P. jackbeardsleyi. The genome size was determined to be 334.818 Mb, which was assembled into 5 linkage groups with a N50 of 67.233 Mb. The BUSCO analysis demonstrated the completeness of the genome assembly and annotation are 95.7% and 92.8%, respectively. The developed high-quality genome will serve as an asset for delving into the genetic mechanisms underlying the invasiveness of P. jackbeardsleyi, thereby offering a crucial theoretical foundation for the prevention and management of Pseudococcidae pests.

The chromosome-level genome assembly of the red swamp crayfish Procambarus clarkii

Red swamp crayfish, Procambarus clarkii, is the most cultured freshwater crayfish species. It attracts significant research attention due to its considerable economic importance. However, the limited availability of genome information has impeded further genetic studies and breeding programs. By utilizing Illumina, PacBio, and Hi-C sequencing technologies, we present a more comprehensive and continuous chromosome-level assembly for P. clarkii than the published one. The final genome size is 4.03 Gb, consisting of 2,358 scaffolds with a N50 of 42.87 Mb. Notably, 3.68 Gb, corresponding to 91.42% of the genome, was anchored to 94 chromosomes. The assembly comprises 70.64% repetitive sequences, including 5.21% tandem repeats and 65.40% transposable elements. Additionally, a total of 4,456 non-coding RNAs and 28,852 protein-coding genes were predicted in the P. clarkii genome, with 96.26% of the genes were annotated. This high-quality genome assembly not only represents a significant improvement for the genome of P. clarkii and provides insights into the unique genome evolution, but also offers valuable information for developing freshwater aquaculture and accelerating genetic breeding.

Distributions, interactions, and dynamics of prokaryotes and phages in a hybrid biological wastewater treatment system

BackgroundUnderstanding the interactions and dynamics of microbiotas within biological wastewater treatment systems is essential for ensuring their stability and long-term sustainability. In this study, we developed a systematic framework employing multi-omics and Hi-C sequencing to extensively investigate prokaryotic and phage communities within a hybrid biofilm and activated sludge system.ResultsWe uncovered distinct distribution patterns, metabolic capabilities, and activities of functional prokaryotes through the analysis of 454 reconstructed prokaryotic genomes. Additionally, we reconstructed a phage catalog comprising 18,645 viral operational taxonomic units (vOTUs) with high length and contiguity using hybrid assembly, and a distinct distribution of phages was depicted between activated sludge (AS) and biofilm. Importantly, 1340 host-phage pairs were established using Hi-C and conventional in silico methods, unveiling the host-determined phage prevalence. The majority of predicted hosts were found to be involved in various crucial metabolic processes, highlighting the potential vital roles of phages in influencing substance metabolism within this system. Moreover, auxiliary metabolic genes (AMGs) related to various categories (e.g., carbohydrate degradation, sulfur metabolism, transporter) were predicted. Subsequent activity analysis emphasized their potential ability to mediate host metabolism during infection. We also profiled the temporal dynamics of phages and their associated hosts using 13-month time-series metagenomic data, further demonstrating their tight interactions. Notably, we observed lineage-specific infection patterns, such as potentially host abundance- or phage/host ratio-driven phage population changes.ConclusionsThe insights gained from this research contribute to the growing body of knowledge surrounding interactions and dynamics of host-phage and pave the way for further exploration and potential applications in the field of microbial ecology.7mTZubWt51_LqpKWPutADaVideo

LncRNA LOC730101 Promotes Darolutamide Resistance in Prostate Cancer by Suppressing miR-1-3p.

Antiandrogen is part of the standard-of-care treatment option for metastatic prostate cancer. However, prostate cancers frequently relapse, and the underlying resistance mechanism remains incompletely understood. This study seeks to investigate whether long non-coding RNAs (lncRNAs) contribute to the resistance against the latest antiandrogen drug, darolutamide. Our RNA sequencing analysis revealed significant overexpression of LOC730101 in darolutamide-resistant cancer cells compared to the parental cells. Elevated LOC730101 levels were also observed in clinical samples of metastatic castration-resistant prostate cancer (CRPC) compared to primary prostate cancer samples. Silencing LOC730101 with siRNA significantly impaired the growth of darolutamide-resistant cells. Additional RNA sequencing analysis identified a set of genes regulated by LOC730101, including key players in the cell cycle regulatory pathway. We further demonstrated that LOC730101 promotes darolutamide resistance by competitively inhibiting microRNA miR-1-3p. Moreover, by Hi-C sequencing, we found that LOC730101 is located in a topologically associating domain (TAD) that undergoes specific gene induction in darolutamide-resistant cells. Collectively, our study demonstrates the crucial role of the lncRNA LOC730101 in darolutamide resistance and its potential as a target for overcoming antiandrogen resistance in CRPC.

Chromosome-level genome assembly of Plagiognathops microlepis based on PacBio HiFi and Hi-C sequencing

Plagiognathops microlepis is an economic freshwater fish in the subfamily Xenocyprinae of Cyprinidae. It is widely distributed in the freshwater ecosystem of China, with moderate economic value and broad development prospects. However, the lack of genomic resources has limited our understanding on the genetic basis, phylogenetic status and adaptive evolution strategies of this fish. Here, we assembled a chromosome-level reference genome of P. microlepis by integrating Pacbio HiFi long-reads, Illumina short-reads and Hi-C sequencing data. The size of this genome is 1004.34 Mb with a contig N50 of 38.80 Mb. Using Hi-C sequencing data, 99.59% of the assembled sequences were further anchored to 24 chromosomes. A total of 578.91 Mb repeat sequences and 28,337 protein-coding genes were predicted in the current genome, of which, 26,929 genes were functionally annotated. This genome provides valuable information for investigating the phylogeny and evolutionary history of cyprinid fishes, as well as the genetic basis of adaptive strategies and special traits in P. microlepis.

Chromosome-level genome assembly of the two-spotted spider mite Tetranychus urticae

The two-spotted spider mite, Tetranychus urticae Koch (Acari: Tetranychidae), is a notorious pest in agriculture that has developed resistance to almost all chemical types used for its control. Here, we assembled a chromosome-level genome for the TSSM using Illumina, Nanopore, and Hi-C sequencing technologies. The assembled contigs had a total length of 103.94 Mb with an N50 of 3.46 Mb, with 87.7 Mb of 34 contigs anchored to three chromosomes. The chromosome-level genome assembly had a BUSCO completeness of 94.8%. We identified 15,604 protein-coding genes, with 11,435 genes that could be functionally annotated. The high-quality genome provides invaluable resources for the genetic and evolutionary study of TSSM.

High-quality reference genome of cowpea beetle Callosobruchus maculatus

Callosobruchus maculatus is one of the most competitive stored grain pests, which causes a great loss to agricultural economy. However, due to an inadequacy of high-quality reference genome, the molecular mechanisms for olfactory and hypoxic adaptations to stored environments are unknown and require to be revealed urgently, which will contribute to the detection and prevention of the invasive pests C. maculatus. Here, we presented a high-quality chromosome-level genome of C. maculatus based on Illumina, Nanopore and Hi-C sequencing data. The total size was 1.2 Gb, and 65.17% (797.47 Mb) of it was identified to be repeat sequences. Among assembled chromosomes, chromosome 10 was considered the X chromosome according to the evidence of reads coverage and homologous genes among species. The current version of high-quality genome provides preferable data resources for the adaptive evolution research of C. maculatus.

Integrated analyses highlight interactions between the three-dimensional genome and DNA, RNA and epigenomic alterations in metastatic prostate cancer

The impact of variations in the three-dimensional structure of the genome has been recognized, but solid cancer tissue studies are limited. Here, we performed integrated deep Hi-C sequencing with matched whole-genome sequencing, whole-genome bisulfite sequencing, 5-hydroxymethylcytosine (5hmC) sequencing and RNA sequencing across a cohort of 80 biopsy samples from patients with metastatic castration-resistant prostate cancer. Dramatic differences were present in gene expression, 5-methylcytosine/5hmC methylation and in structural variation versus mutation rate between A and B (open and closed) chromatin compartments. A subset of tumors exhibited depleted regional chromatin contacts at the AR locus, linked to extrachromosomal circular DNA (ecDNA) and worse response to AR signaling inhibitors. We also identified topological subtypes associated with stark differences in methylation structure, gene expression and prognosis. Our data suggested that DNA interactions may predispose to structural variant formation, exemplified by the recurrent TMPRSS2–ERG fusion. This comprehensive integrated sequencing effort represents a unique clinical tumor resource.

A near complete genome assembly of the East Friesian sheep genome

Advancements in sequencing have enabled the assembly of numerous sheep genomes, significantly advancing our understanding of the link between genetic variation and phenotypic traits. However, the genome of East Friesian sheep (Ostfriesisches Milchschaf), a key high-yield milk breed, remains to be fully assembled. Here, we constructed a near-complete and gap-free East Friesian genome assembly using PacBio HiFi, ultra-long ONT and Hi-C sequencing. The resulting genome assembly spans approximately 2.96 Gb, with a contig N50 length of 104.1 Mb and only 164 unplaced sequences. Remarkably, our assembly has captured 41 telomeres and 24 centromeres. The assembled sequence is of high quality on completeness (BUSCO score: 97.1%) and correctness (QV: 69.1). In addition, a total of 24,580 protein-coding genes were predicted, of which 97.2% (23,891) carried at least one conserved functional domain. Collectively, this assembly provides not only a near T2T gap-free genome, but also provides a valuable genetic resource for comparative genome studies of sheep and will serve as an important tool for the sheep research community.

Chromosome-level genome assembly of an oligophagous leaf beetle Ophraella communa (Coleoptera: Chrysomelidae)

The leaf beetle Ophraella communa LeSage (Coleoptera: Chrysomelidae) is an effective biological control agent of the common ragweed. Here, we assembled a chromosome-level genome of the O. communa by combining Illumina, Nanopore, and Hi-C sequencing technologies. The genome size of the final genome assembly is 733.1 Mb, encompassing 17 chromosomes, with an improved contig N50 of 7.05 Mb compared to the original version. Genome annotation reveals 25,873 protein-coding genes, with functional annotations available for 22,084 genes (85.35%). Non-coding sequence annotation identified 204 rRNAs, 626 tRNAs, and 1791 small RNAs. Repetitive elements occupy 414.41 Mb, constituting 57.76% of the genome. This high-quality genome is fundamental for advancing biological control strategies employing O. communa.

Genome Assembly and Structural Variation Analysis of Luffa acutangula Provide Insights on Flowering Time and Ridge Development.

Luffa spp. is an important worldwide cultivated vegetable and medicinal plant from the Cucurbitaceae family. In this study, we report a high-quality chromosome-level genome of the high-generation inbred line SG261 of Luffa acutangula. The genomic sequence was determined by PacBio long reads, Hi-C sequencing reads, and 10× Genomics sequencing, with an assembly size of 739.82 Mb, contig N50 of 18.38 Mb, and scaffold N50 of 56.08 Mb. The genome of L. acutangula SG261 was predicted to contain 27,312 protein-coding genes and 72.56% repetitive sequences, of which long terminal repeats (LTRs) were an important form of repetitive sequences, accounting for 67.84% of the genome. Phylogenetic analysis reveals that L. acutangula evolved later than Luffa cylindrica, and Luffa is closely related to Momodica charantia. Comparing the genome of L. acutangula SG261 and L. cylindrica with PacBio data, 67,128 high-quality structural variations (SVs) and 55,978 presence-absence variations (PAVs) were identified in SG261, resulting in 2424 and 1094 genes with variation in the CDS region, respectively, and there are 287 identical genes affected by two different structural variation analyses. In addition, we found that the transcription factor FY (FLOWERING LOCUS Y) families had a large expansion in L. acutangula SG261 (flowering in the morning) compared to L. cylindrica (flowering in the afternoon), which may result in the early flowering time in L. acutangula SG261. This study provides valuable reference for the breeding of and pan-genome research into Luffa species.

Chromosome-Scale Genome Assembly Provides Insights into Fresh Pine Wood Decay Strategies of the Wolfiporia hoelen

The sclerotia of Wolfiporia hoelen (Fr.) Y.C. Dai & V. Papp is an important traditional Chinese medicine with diverse pharmacological properties. This study utilized a combination of PacBio Long-Read Sequencing, Illumina Short-Read Sequencing, and Hi-C Sequencing to generate a high-quality chromosome-level genome assembly of a W. hoelen strain Minling A5. There were 112 contigs in the genome, with 62.95 Mb in total length and 4.21 Mb in length for the contig N50. The average GC content was 51.89%. Based on Hi-C data, we corrected the CCS data and scaffolded them into 14 pseudo-chromosomes. The genome contained 44.37% repetitive sequences and 12,670 protein-coding genes, 86.53% (10,963) of which could be functionally annotated in at least one of the KOG, GO, Pfam, Swissprot, TrEMBL, NR, and KEGG databases. In addition, 240 transfer RNAs, 97 ribosomal RNAs, and 103 other non-coding RNAs were identified in the W. hoelen genome. A total of 755 pseudogenes were also identified, with an average length of 2665.51 bp. Further, there were 398, 100, 2837, 519, and 2068 genes annotated by CAZymes, TCDB, PHI, P450, and DFVF databases, respectively. One notable attribute of W. hoelen is its capacity to thrive in a substrate of fresh pine sawdust. Through an analysis of the growth on various pure wood sawdust culture media, we found that the growth of W. hoelen and Sparassis latifolia on pine sawdust was similar to that on broad-leaved wood sawdust, while the growth of Pleurotus ostreatus, P. eryngii, and Cyclocybe aegerita was slower than that on broad-leaved wood sawdust. By the functional annotation analysis of orthogroups in these five mushroom-forming fungi, it was determined that 645 orthogroups were specifically common in W. hoelen and S. latifolia. The genes in these specific orthogroups were significantly enriched in 12 pathways, including steroid biosynthesis, biosynthesis of antibiotics, and tyrosine metabolism. The high-quality genome and comparative genome analysis results significantly contribute to advancing our foundational knowledge of W. hoelen biology, while also offering valuable insights for the development of innovative biotechnological approaches aimed at enhancing the efficient and sustainable utilization of Pinus.

ENT3C: an entropy-based similarity measure for Hi-C and micro-C derived contact matrices.

Hi-C and micro-C sequencing have shed light on the profound importance of 3D genome organization in cellular function by probing 3D contact frequencies across the linear genome. The resulting contact matrices are extremely sparse and susceptible to technical- and sequence-based biases, making their comparison challenging. The development of reliable, robust and efficient methods for quantifying similarity between contact matrices is crucial for investigating variations in the 3D genome organization in different cell types or under different conditions, as well as evaluating experimental reproducibility. We present a novel method, ENT3C, which measures the change in pattern complexity in the vicinity of contact matrix diagonals to quantify their similarity. ENT3C provides a robust, user-friendly Hi-C or micro-C contact matrix similarity metric and a characteristic entropy signal that can be used to gain detailed biological insights into 3D genome organization.

Chromosome-level genome assembly of Guide Black-Fur sheep (Ovis aries)

Guide Black-Fur sheep (GD) is a breed of Tibetan sheep (Ovis aries) that lives in the Qinghai–Tibetan plateau region at an altitude of over 4,000 m. However, a lack of genomic information has made it difficult to understand the high-altitude adaptation of these sheep. We sequenced and assembled the GD reference genome using PacBio, Hi-C, and Illumina sequencing technologies. The final assembled genome size was 2.73 Gb, with a contig N50 of 20.30 Mb and a scaffold N50 of 107.63 Mb. The genome is predicted to contain 20,759 protein-coding genes, of which 98.42 have functional annotations. Repeat elements account for approximately 52.2% of the genomic landscape. The completeness of the GD genome assembly is highlighted by a BUSCO score of 93.1%. This high-quality genome assembly provides a critical resource for future molecular breeding and genetic improvement of Tibetan sheep.

Chromosome genome assembly of the Camphora longepaniculata (Gamble) with PacBio and Hi-C sequencing data.

Camphora longepaniculata, a crucial commercial crop and a fundamental component of traditional Chinese medicine, is renowned for its abundant production of volatile terpenoids. However, the lack of available genomic information has hindered pertinent research efforts in the past. To bridge this gap, the present study aimed to use PacBio HiFi, short-read, and highthroughput chromosome conformation capture sequencing to construct a chromosome-level assembly of the C. longepaniculata genome. With twelve chromosomes accounting for 99.82% (766.69 Mb) of the final genome assembly, which covered 768.10 Mb, it was very complete. Remarkably, the assembly's contig and scaffold N50 values are exceptional as well-41.12 and 63.78 Mb, respectively-highlighting its excellent quality and intact structure. Furthermore, a total of 39,173 protein-coding genes were predicted, with 38,766 (98.96%) of them being functionally annotated. The completeness of the genome was confirmed by the Benchmarking Universal Single-Copy Ortholog evaluation, which revealed 99.01% of highly conserved plant genes. As the first comprehensive assembly of the C. longepaniculata genome, it provides a crucial starting point for deciphering the complex pathways involved in terpenoid production. Furthermore, this excellent genome serves as a vital resource for upcoming research on the breeding and genetics of C. longepaniculata.

Haplotype-resolved assembly of a pig genome using single-sperm sequencing

Single gamete cell sequencing together with long-read sequencing can reliably produce chromosome-level phased genomes. In this study, we employed PacBio HiFi and Hi-C sequencing on a male Landrace pig, coupled with single-sperm sequencing of its 102 sperm cells. A haplotype assembly method was developed based on long-read sequencing and sperm-phased markers. The chromosome-level phased assembly showed higher phasing accuracy than methods that rely only on HiFi reads. The use of single-sperm sequencing data enabled the construction of a genetic map, successfully mapping the sperm motility trait to a specific region on chromosome 1 (105.40–110.70 Mb). Furthermore, with the assistance of Y chromosome-bearing sperm data, 26.16 Mb Y chromosome sequences were assembled. We report a reliable approach for assembling chromosome-level phased genomes and reveal the potential of sperm population in basic biology research and sperm phenotype research.