Structural variants (SVs, including large-scale insertions, deletions, inversions, and translocations) significantly impact the functions of genes in the microbial genome, and SVs in the microbiome are associated with diverse biological processes and human diseases. With the advancements in sequencing and bioinformatics technologies, increasingly, sequencing data and analysis tools are already being extensively utilized for microbiome SV analyses, leading to a higher demand for more dedicated SV analysis workflows. Moreover, due to the unique detection biases of various sequencing technologies, including short-read sequencing (such as Illumina platforms) and long-read sequencing (e.g., Oxford Nanopore and PacBio), SV discovery based on multiple platforms is necessary to comprehensively identify the wide variety of SVs. Here, we establish an integrated pipeline MetaSVs combining Nanopore long reads and Illumina short reads to analyze SVs in the microbial genomes from gut microbiome and further identify differential SVs that can be reflective of metabolic differences. Our pipeline provides researchers easy access to SVs and relevant metabolites in the microbial genomes without the requirement of specific technical expertise, which is particularly useful to researchers interested in metagenomic SVs but lacking sophisticated bioinformatic knowledge.

Structural variants (SVs) are infrequently observed in Duchenne muscular dystrophy (DMD), a condition mainly marked by deletions and point mutations in the DMD gene. SVs in DMD remain difficult to reliably detect due to the limited SV-detection capacity of conventionally used short-read sequencing technology. Herein, we present a family, a boy and his mother, with clinical signs of muscular dystrophy, elevated creatinine kinase levels, and intellectual disability. A muscle biopsy from the boy showed dystrophin deficiency. Routine molecular techniques failed to detect abnormalities in the DMD gene, however, dystrophin mRNA transcripts analysis revealed an absence of exons 59 to 79. Subsequent long-read whole-genome sequencing identified a rare complex structural variant, a 77 kb novel intragenic inversion, and a balanced translocation t(X;1)(p21.2;p13.3) rearrangement within the DMD gene, expanding the genetic spectrum of dystrophinopathy. Our findings suggested that SVs should be considered in cases where conventional molecular techniques fail to identify pathogenic variants.

Alfalfa (Medicago sativa L.) is one of the most important forage legumes in the world, including autotetraploid (M. sativa ssp. sativa) and diploid alfalfa (M. sativa ssp. caerulea, progenitor of autotetraploid alfalfa). Here, we reported a high-quality genome of ZW0012 (diploid alfalfa, 769 Mb, contig N50 = 5.5 Mb), which was grouped into the Northern group in population structure analysis, suggesting that our genome assembly filled a major gap among the members of M. sativa complex. During polyploidization, large phenotypic differences occurred between diploids and tetraploids, and the genetic information underlying its massive phenotypic variations remains largely unexplored. Extensive structural variations (SVs) were identified between ZW0012 and XinJiangDaYe (an autotetraploid alfalfa with released genome). We identified 71 ZW0012-specific PAV genes and 1296 XinJiangDaYe-specific PAV genes, mainly involved in defence response, cell growth, and photosynthesis. We have verified the positive roles of MsNCR1 (a XinJiangDaYe-specific PAV gene) in nodulation using an Agrobacterium rhizobia-mediated transgenic method. We also demonstrated that MsSKIP23_1 and MsFBL23_1 (two XinJiangDaYe-specific PAV genes) regulated leaf size by transient overexpression and virus-induced gene silencing analysis. Our study provides a high-quality reference genome of an important diploid alfalfa germplasm and a valuable resource of variation landscape between diploid and autotetraploid, which will facilitate the functional gene discovery and molecular-based breeding for the cultivars in the future.

Genomic structural variants (SVs) constitute a significant proportion of genetic variation in the genome. The rapid development of long-reads sequencing has facilitated the detection of long-fragment SVs. There is no published study to detect SVs using long-read data from sheep. We applied a long-read mapping approach to detect SVs and characterized a total of 30,771 insertions, deletions, inversions and translocations. We identified 716, 916, 842 and 303 specific SVs in Southdown sheep, Alpine merino sheep, Qilian White Tibetan sheep and Oula sheep, respectively. We annotated these SVs and found that these SV-related genes were primarily enriched in the well-established pathways involved in the regulation of the immune system, growth and development and environmental adaptability. We detected and annotated SVs based on NGS resequencing data to validate the accuracy based on third-generation detection. Moreover, five candidate SVs were verified using the PCR method in 50 sheep. Our study is the first to use a long-reads sequencing approach to construct a novel structural variation map in sheep. We have completed a preliminary exploration of the potential effects of SVs on sheep.

DNA double-strand breaks (DSBs) are harmful to mammalian cells and a few of them can cause cell death. Accumulating DSBs in these cells to analyze their genomic distribution and their potential impact on chromatin structure is difficult. In this study, we used CRISPR to generate Ku80-/- human cells and arrested the cells in G1 phase to accumulate DSBs before conducting END-seq and Nanopore analysis. Our analysis revealed that DNA with high methylation level accumulates DSB hotspots in Ku80-/- human cells. Furthermore, we identified chromosome structural variants (SVs) using Nanopore sequencing and observed a higher number of SVs in Ku80-/- human cells. Based on our findings, we suggest that the high efficiency of Ku80 knockout in human HCT116 cells makes it a promising model for characterizing SVs in the context of 3D chromatin structure and studying the alternative-end joining (Alt-EJ) DSB repair pathway.

Xishuangbanna (XIS) cucumber (Cucumis sativus L. var. xishuangbannanesis) is a semiwild variety originating from low latitude tropic areas, and therefore shows extreme cold sensitivity and heat tolerance. Here, we mapped the quantitative trait loci (QTLs) that control the cold sensitivity and heat tolerance of XIS cucumber seedlings. Using bulked segregant analysis (BSA), we identified three QTLs (HTT1.1, HTT3.1, and HTT3.2, with a total length of 11.98 Mb) for heat tolerance and two QTLs (LTT6.1 and LTT6.2, with a total length of 8.74 Mb) for cold sensitivity. The QTL LTT6.1 was then narrowed down to a length of 641 kb by using kompetitive allele-specific PCR (KASP) markers. Based on structural variants (SVs) and single-nucleotide polymorphisms (SNPs), we found the LTT6.1 is covered by a high divergent region including a 50 kb deletion in the XIS49 genome, which affects the gene structure of lipase abhydrolase domain containing 6 (ABHD6, Csa_6G032560). Accordingly, there is a very big difference in lipid composition, but not in other osmoprotectants like free amino acids and fatty acids, between XIS49 and cultivated cucumber CL. Moreover, we calculated the composite likelihood ratio (CLR) and identified selective sweeps from 115 resequencing data, and found that lipid- and fatty-acid-related processes are major aspects in the domestication of the XIS group cucumber. LTT6.1 is a particularly special region positioned nearby lipid-related selective sweeps. These studies above suggested that the lipid-related domestication of XIS cucumbers should account for their extreme cold sensitivity.

The domestication of Brassica oleracea has resulted in diverse morphological types with distinct patterns of organ development. Here we report a graph-based pan-genome of B. oleracea constructed from high-quality genome assemblies of different morphotypes. The pan-genome harbors over 200 structural variant hotspot regions enriched in auxin- and flowering-related genes. Population genomic analyses revealed that early domestication of B. oleracea focused on leaf or stem development. Gene flows resulting from agricultural practices and variety improvement were detected among different morphotypes. Selective-sweep and pan-genome analyses identified an auxin-responsive small auxin up-regulated RNA gene and a CLAVATA3/ESR-RELATED family gene as crucial players in leaf-stem differentiation during the early stage of B. oleracea domestication and the BoKAN1 gene as instrumental in shaping the leafy heads of cabbage and Brussels sprouts. Our pan-genome and functional analyses further revealed that variations in the BoFLC2 gene play key roles in the divergence of vernalization and flowering characteristics among different morphotypes, and variations in the first intron of BoFLC3 are involved in fine-tuning the flowering process in cauliflower. This study provides a comprehensive understanding of the pan-genome of B. oleracea and sheds light on the domestication and differential organ development of this globally important crop species.

Structural variants of high-nuclearity clusters are extremely important for their modular assembly study and functional expansion, yet the synthesis of such giant structural variants remains a great challenge. Herein, we prepared a lantern-type giant polymolybdate cluster (L-Mo132 ) containing equal metal nuclearity with the famous Keplerate type Mo132 (K-Mo132 ). The skeleton of L-Mo132 features a rare truncated rhombic triacontrahedron, which is totally different with the truncated icosahedral K-Mo132 . To the best of our knowledge, this is the first time to observe such structural variants in high-nuclearity cluster built up of more than 100 metal atoms. Scanning transmission electron microscopy reveals that L-Mo132 has good stability. More importantly, because the pentagonal [Mo6 O27 ]n- building blocks in L-Mo132 are concave instead of convex in the outer face, it contains multiple terminal coordinated water molecules on its outer surface, which make it expose more active metal sites to display superior phenol oxidation performance, which is more higher than that of K-Mo132 coordinated in M=O bonds on the outer surface.

Understanding how genetic variants impact molecular phenotypes is a key goal of functional genomics, currently hindered by reliance on a single haploid reference genome. Here, we present the EN-TEx resource of 1,635 open-access datasets from four donors (∼30 tissues × ∼15 assays). The datasets are mapped to matched, diploid genomes with long-read phasing and structural variants, instantiating a catalog of >1 million allele-specific loci. These loci exhibit coordinated activity along haplotypes and are less conserved than corresponding, non-allele-specific ones. Surprisingly, a deep-learning transformer model can predict the allele-specific activity based only on local nucleotide-sequence context, highlighting the importance of transcription-factor-binding motifs particularly sensitive to variants. Furthermore, combining EN-TEx with existing genome annotations reveals strong associations between allele-specific and GWAS loci. It also enables models for transferring known eQTLs to difficult-to-profile tissues (e.g., from skin to heart). Overall, EN-TEx provides rich data and generalizable models for more accurate personal functional genomics.

Due to the difficulty in accurately identifying structural variants (SVs) across genomes, their impact on cis-regulatory divergence of closely related species, especially fish, remains to be explored. Recently identified broad H3K4me3 domains are essential for the regulation of genes involved in several biological processes. However, the role of broad H3K4me3 domains in phenotypic divergence remains poorly understood. Siniperca chuatsi and S. scherzeri are closely related but divergent in several phenotypic traits, making them an ideal model to study cis-regulatory evolution in sister species. Here, we generated chromosome-level genomes of S. chuatsi and S. scherzeri, with assembled genome sizes of 716.35 and 740.54 Mb, respectively. The evolutionary histories of S. chuatsi and S. scherzeri were studied by inferring dynamic changes in ancestral population sizes. To explore the genetic basis of adaptation in S. chuatsi and S. scherzeri, we performed gene family expansion and contraction analysis and identified positively selected genes (PSGs). To investigate the role of SVs in cis-regulatory divergence of closely related fish species, we identified high-quality SVs as well as divergent H3K27ac and H3K4me3 domains in the genomes of S. chuatsi and S. scherzeri. Integrated analysis revealed that cis-regulatory divergence caused by SVs played an essential role in phenotypic divergence between S. chuatsi and S. scherzeri. Additionally, divergent broad H3K4me3 domains were mostly associated with cancer-related genes in S. chuatsi and S. scherzeri and contributed to their phenotypic divergence.
由于难以在基因组中全面地鉴定结构变异，我们仍然缺乏对结构变异在近缘物种，特别是鱼类近缘物种顺式调控趋异中作用的认知。近期的研究发现，宽H3K4me3峰在多个生物学过程的基因调控中起关键作用。然而，宽H3K4me3峰在物种表型趋异中的作用仍然未知。翘嘴鳜（ Siniperca chuatsi）和斑鳜（ Siniperca scherzeri）是一对姐妹种，但是它们在多个性状上有很大差异。这使得它们成为研究顺式调控演化的良好模型。我们构建了染色体级别的翘嘴鳜和斑鳜基因组，组装大小分别为716.35 Mb和740.54 Mb。通过推断历史种群大小，我们分析了翘嘴鳜和斑鳜的演化历史。通过基因家族扩张和收缩分析，以及鉴定正选择基因，我们解析了翘嘴鳜和斑鳜适应环境的遗传学基础。为了阐明结构变异在近缘鱼类物种顺式调控趋异中的作用，我们在翘嘴鳜和斑鳜基因组中鉴定了高质量的结构变异，以及差异H3K27ac和H3K4me3峰。联合分析发现，由结构变异引起的顺式调控趋异在翘嘴鳜和斑鳜表型趋异中起关键作用。此外，在翘嘴鳜和斑鳜中，差异宽H3K4me3峰主要与癌症相关基因相关联，并导致这两个物种的表型趋异。.