A complete, telomere-to-telomere human genome sequence presents new opportunities for evolutionary genomics

MangroveDB: A Comprehensive Online Database for Mangroves Based on Multi-Omics Data

Mangroves are dominant flora of intertidal zones along tropical and subtropical coastline around the world that offer important ecological and economic value. Recently, the genomes of mangroves have been decoded, and massive omics data were generated and deposited in the public databases. Reanalysis of multi-omics data can provide new biological insights excluded in the original studies. However, the requirements for computational resource and lack of bioinformatics skill for experimental researchers limit the effective use of the original data. To fill this gap, we uniformly processed 942 transcriptome data, 386 whole-genome sequencing data, and provided 13 reference genomes and 40 reference transcriptomes for 53 mangroves. Finally, we built an interactive web-based database platform MangroveDB (https://github.com/Jasonxu0109/MangroveDB), which was designed to provide comprehensive gene expression datasets to facilitate their exploration and equipped with several online analysis tools, including principal components analysis, differential gene expression analysis, tissue-specific gene expression analysis, GO and KEGG enrichment analysis. MangroveDB not only provides query functions about genes annotation, but also supports some useful visualization functions for analysis results, such as volcano plot, heatmap, dotplot, PCA plot, bubble plot, population structure, and so on. In conclusion, MangroveDB is a valuable resource for the mangroves research community to efficiently use the massive public omics datasets.

Integrated analysis of the complete sequence of a macaque genome

The crab-eating macaques (Macaca fascicularis) and rhesus macaques (Macaca mulatta) are pivotal in biomedical and evolutionary research1-3. However, their genomic complexity and interspecies genetic differences remain unclear4. Here, we present a complete genome assembly of a crab-eating macaque, revealing 46% fewer segmental duplications and 3.83 times longer centromeres than those of humans5,6. We also characterize 93 large-scale genomic differences between macaques and humans at a single-base-pair resolution, highlighting their impact on gene regulation in primate evolution. Using ten long-read macaque genomes, hundreds of short-read macaque genomes and full-length transcriptome data, we identified roughly 2 Mbp of fixed-genetic variants, roughly 240 Mbp of complex loci, 16.76 Mbp genetic differentiation regions and 110 alternative splice events, potentially associated with various phenotypic differences between the two macaque species. In summary, the integrated genetic analysis enhances understanding of lineage-specific phenotypes, adaptation and primate evolution, thereby improving their biomedical applications in human disease research.

Advancing evolutionary medicine with complete primate genomes and advanced biotechnologies

Evolutionary medicine, which integrates evolutionary biology and medicine, significantly enhances our understanding of human traits and disease susceptibility. However, previous studies in this field have often focused on single-nucleotide variants due to technological limitations in characterizing complex genomic regions, hindering the comprehensive analyses of their evolutionary origins and clinical significance. In this review, we summarize recent advancements in complete telomere-to-telomere (T2T), primate genomes and other primate resources, and illustrate how these resources facilitate the research of complex regions. We focus on several biomedically relevant regions to examine the relationship between primate genome evolution and human diseases. We also highlight the potentials of high-throughput functional genomic technologies for assessing candidate loci. Finally, we discuss future directions for primate research within the context of evolutionary medicine.

A refined analysis of Neanderthal-introgressed sequences in modern humans with a complete reference genome

BACKGROUND

Leveraging long-read sequencing technologies, the first complete human reference genome, T2T-CHM13, corrects assembly errors in previous references and resolves the remaining 8% of the genome. While studies on archaic admixture in modern humans have so far relied on the GRCh37 reference due to the availability of archaic genome data, the impact of T2T-CHM13 in this field remains unexplored.

RESULTS

We remap the sequencing reads of the high-quality Altai Neanderthal and Denisovan genomes onto GRCh38 and T2T-CHM13. Compared to GRCh37, we find that T2T-CHM13 significantly improves read mapping quality in archaic samples. We then apply IBDmix to identify Neanderthal-introgressed sequences in 2504 individuals from 26 geographically diverse populations using different reference genomes. We observe that commonly used pre-phasing filtering strategies in public datasets substantially influence archaic ancestry determination, underscoring the need for careful filter selection. Our analysis identifies approximately 51 Mb of Neanderthal sequences unique to T2T-CHM13, predominantly in genomic regions where GRCh38 and T2T-CHM13 assemblies diverge. Additionally, we uncover novel instances of population-specific archaic introgression in diverse populations, spanning genes involved in metabolism, olfaction, and ion-channel function. Finally, to facilitate the exploration of archaic alleles and adaptive signals in human genomics and evolutionary research, we integrate these introgressed sequences and adaptive signals across all reference genomes into a visualization database, ASH ( www.arcseqhub.com ).

CONCLUSIONS

Our study enhances the detection of archaic variations in modern humans, highlights the importance of utilizing the T2T-CHM13 reference, and provides novel insights into the functional consequences of archaic hominin admixture.

CUT&Tag Identifies Repetitive Genomic Loci that are Excluded from ChIP Assays

Determining the genomic localization of chromatin features is an essential aspect of investigating gene expression control, and ChIP-Seq has long been the gold standard technique for interrogating chromatin landscapes. Recently, the development of alternative methods, such as CUT&Tag, have provided researchers with alternative strategies that eliminate the need for chromatin purification, and allow for in situ investigation of histone modifications and chromatin bound factors. Mindful of technical differences, we set out to investigate whether distinct chromatin modifications were equally compatible with these different chromatin interrogation techniques. We found that ChIP-Seq and CUT&Tag performed similarly for modifications known to reside at gene regulatory regions, such as promoters and enhancers, but major differences were observed when we assessed enrichment over heterochromatin-associated loci. Unlike ChIP-Seq, CUT&Tag detects robust levels of H3K9me3 at a substantial number of repetitive elements, with especially high sensitivity over evolutionarily young retrotransposons. IAPEz-int elements for example, exhibited underrepresentation in mouse ChIP-Seq datasets but strong enrichment using CUT&Tag. Additionally, we identified several euchromatin-associated proteins that co-purify with repetitive loci and are similarly depleted when applying ChIP-based methods. This study reveals that our current knowledge of chromatin states across the heterochromatin portions of the mammalian genome is extensively incomplete, largely due to limitations of ChIP-Seq. We also demonstrate that newer in situ chromatin fragmentation-based techniques, such as CUT&Tag and CUT&RUN, are more suitable for studying chromatin modifications over repetitive elements and retrotransposons.

EZH2 inhibition enhances T cell immunotherapies by inducing lymphoma immunogenicity and improving T cell function

T cell-based immunotherapies have demonstrated effectiveness in treating diffuse large B cell lymphoma (DLBCL) and follicular lymphoma (FL) but predicting response and understanding resistance remains a challenge. To address this, we developed syngeneic models reflecting the genetics, epigenetics, and immunology of human FL and DLBCL. We show that EZH2 inhibitors reprogram these models to re-express T cell engagement genes and render them highly immunogenic. EZH2 inhibitors do not harm tumor-controlling T cells or CAR-T cells. Instead, they reduce regulatory T cells, promote memory chimeric antigen receptor (CAR) CD8 phenotypes, and reduce exhaustion, resulting in a decreased tumor burden. Intravital 2-photon imaging shows increased CAR-T recruitment and interaction within the tumor microenvironment, improving lymphoma cell killing. Therefore, EZH2 inhibition enhances CAR-T cell efficacy through direct effects on CAR-T cells, in addition to rendering lymphoma B cells immunogenic. This approach is currently being evaluated in two clinical trials, NCT05934838 and NCT05994235, to improve immunotherapy outcomes in B cell lymphoma patients.

Analyses of Human Genetic Data to Identify Clinically Relevant Domains of Neuroligins

Background/Objectives: Neuroligins (NLGNs) are postsynaptic adhesion molecules critical for neuronal development that are highly associated with autism spectrum disorder (ASD). Here, we provide an overview of the literature on NLGN rare variants. In addition, we introduce a new approach to analyze human variation within NLGN genes to identify sensitive regions that have an increased frequency of ASD-associated variants to better understand NLGN function. Methods: To identify critical protein subdomains within the NLGN gene family, we developed an algorithm that assesses tolerance to missense mutations in human genetic variation by comparing clinical variants from ClinVar to reference variants from gnomAD. This approach provides tolerance values to subdomains within the protein. Results: Our algorithm identified several critical regions that were conserved across multiple NLGN isoforms. Importantly, this approach also identified a previously reported cluster of pathogenic variants in NLGN4X (also conserved in NLGN1 and NLGN3) as well as a region around the highly characterized NLGN3 R451C ASD-associated mutation. Additionally, we highlighted other, as of yet, uncharacterized regions enriched with mutations. Conclusions: The systematic analysis of NLGN ASD-associated variants compared to variants identified in the unaffected population (gnomAD) reveals conserved domains in NLGN isoforms that are tolerant to variation or are enriched in clinically relevant variants. Examination of databases also allows for predictions of the presumed tolerance to loss of an allele. The use of the algorithm we developed effectively allowed the evaluation of subdomains of NLGNs and can be used to examine other ASD-associated genes.

Comprehensive evaluation and guidance of structural variation detection tools in chicken whole genome sequence data

BACKGROUND

Structural variations (SVs) are widespread across genome and have a great impact on evolution, disease, and phenotypic diversity. Despite the development of numerous bioinformatic tools, commonly referred to as SV callers, tailored for detecting SVs using whole genome sequence (WGS) data and employing diverse algorithms, their performance necessitates rigorous evaluation with real data and validated SVs. Moreover, a considerable proportion of these tools have been primarily designed and optimized using human genome data. Consequently, their applicability and performance in Avian species, characterized by smaller genomes and distinct genomic architectures, remain inadequately assessed.

RESULTS

We performed a comprehensive assessment of the performance of ten widely used SV callers using population-level real genomic data with the validated five common types of SVs. The performance of SV callers varies with the types and sizes of SVs. As compared with other tools, GRIDSS, Lumpy, Wham, and Manta present better detection accuracy. Pindel can detect more small SVs than others. CNVnator and CNVkit can detect more medium and large copy number variations. Given the poor consistency among different SV callers, the combination calling strategy is not recommended. All tools show poor ability in the detection of insertions (especially with size > 150 bp). At least 50× read depth is required to detect more than 80% of the SVs for most tools.

CONCLUSIONS

This study highlights the importance and necessity of using real sequencing data, rather than simulated data only, with validated SVs for SV caller evaluation. Some practical guidance and suggestions are provided for SV detection in future researches.

Proteogenomic characterization of skull-base chordoma

Skull-base chordoma is a rare, aggressive bone cancer with a high recurrence rate. Despite advances in genomic studies, its molecular characteristics and effective therapies remain unknown. Here, we conduct integrative genomics, transcriptomics, proteomics, and phosphoproteomics analyses of 187 skull-base chordoma tumors. In our study, chromosome instability is identified as a prognostic predictor and potential therapeutic target. Multi-omics data reveals downstream effects of chromosome instability, with RPRD1B as a putative target for radiotherapy-resistant patients. Chromosome 1q gain, associated with chromosome instability and upregulated mitochondrial functions, lead to poorer clinical outcomes. Immune subtyping identify an immune cold subtype linked to chromosome 9p/10q loss and immune evasion. Proteomics-based classification reveals subtypes (P-II and P-III) with high chromosome instability and immune cold features, with P-II tumors showing increased invasiveness. These findings, confirmed in 17 paired samples, provide insights into the biology and treatment of skull-base chordoma.

Gene expression analysis identifies hub genes and pathways distinguishing fatal from survivor outcomes of Ebola virus disease

The Ebola virus poses a severe public health threat, yet understanding factors influencing disease outcomes remains incomplete. Our study aimed to identify critical pathways and hub genes associated with fatal and survivor Ebola disease outcomes. We analyzed differentially expressed hub genes (DEGs) between groups with fatal and survival outcomes, as well as a healthy control group. We conducted additional analysis to determine the functions and pathways associated with these DEGs. We found 13,198 DEGs in the fatal and 12,039 DEGs in the survival group compared to healthy controls, and 1873 DEGs in the acute fatal and survivor groups comparison. Upregulated DEGs in the comparison between the acute fatal and survivor groups were linked to ECM receptor interaction, complement and coagulation cascades, and PI3K-Akt signaling. Upregulated hub genes identified from the acute fatal and survivor comparison (FGB, C1QA, SERPINF2, PLAT, C9, SERPINE1, F3, VWF) were enriched in complement and coagulation cascades; the downregulated hub genes (IL1B, 1L17RE, XCL1, CXCL6, CCL4, CD8A, CD8B, CD3D) were associated with immune cell processes. Hub genes CCL2 and F2 were unique to fatal outcomes, while CXCL1, HIST1H4F, and IL1A were upregulated hub genes unique to survival outcomes compared to healthy controls. Our results demonstrate for the first time the association of EVD outcomes to specific hub genes and their associated pathways and biological processes. The identified hub genes and pathways could help better elucidate Ebola disease pathogenesis and contribute to the development of targeted interventions and personalized treatment for distinct EVD outcomes.

Low-input PacBio sequencing generates high-quality individual fly genomes and characterizes mutational processes

Long-read sequencing, exemplified by PacBio, revolutionizes genomics, overcoming challenges like repetitive sequences. However, the high DNA requirement ( > 1 µg) is prohibitive for small organisms. We develop a low-input (100 ng), low-cost, and amplification-free library-generation method for PacBio sequencing (LILAP) using Tn5-based tagmentation and DNA circularization within one tube. We test LILAP with two Drosophila melanogaster individuals, and generate near-complete genomes, surpassing preexisting single-fly genomes. By analyzing variations in these two genomes, we characterize mutational processes: complex transpositions (transposon insertions together with extra duplications and/or deletions) prefer regions characterized by non-B DNA structures, and gene conversion of transposons occurs on both DNA and RNA levels. Concurrently, we generate two complete assemblies for the endosymbiotic bacterium Wolbachia in these flies and similarly detect transposon conversion. Thus, LILAP promises a broad PacBio sequencing adoption for not only mutational studies of flies and their symbionts but also explorations of other small organisms or precious samples.

Genomic analyses provide insights into sex differentiation of tetraploid strawberry (Fragaria moupinensis)

The strawberry genus, Fragaria, exhibits a wide range of sexual systems and natural ploidy variation. Nearly, all polyploid strawberry species exhibit separate sexes (dioecy). Research has identified the sex-determining sequences as roughly conserved but with repeatedly changed genomic locations across octoploid strawberries. However, it remains unclear whether tetraploid wild strawberries evolved dioecy independently or shared a common origin with octoploid strawberries. In this study, we investigated the sex determinants of F. moupinensis, a dioecious plant with heterogametic females (ZW). Utilizing a combination of haplotype-resolved genome sequencing of the female F. moupinensis, k-mer-based and coverage-based genome-wide association studies (GWAS), and transcriptomic analysis, we discovered a non-recombining, approximately 33.6 kb W-specific region on chromosome 2a. Within this region, only one candidate sex-determining gene (FmoAFT) was identified. Furthermore, an extensive resequencing of the entire Fragaria genus indicated that the W-specific region displays conservative female specificity across all tetraploid species. This observation suggests that dioecy evolved independently in tetraploid and octoploid strawberries. Moreover, employing virus-induced gene silencing (VIGS), we knocked down the expression of the FmoAFT homologue transcript in cultivated strawberries, revealing its potential role in promoting female functions during early carpel development. We also applied DNA affinity purification sequencing (DAP-seq) and yeast one-hybrid assays to identify potential direct targets of FmoAFT. These insights shed new light on the genetic basis and evolutionary history of sex determination in strawberries, thereby facilitating the formulation of strategies to manipulate sex determination in breeding programs.

NextPolish2: A Repeat-aware Polishing Tool for Genomes Assembled Using HiFi Long Reads

The high-fidelity (HiFi) long-read sequencing technology developed by PacBio has greatly improved the base-level accuracy of genome assemblies. However, these assemblies still contain base-level errors, particularly within the error-prone regions of HiFi long reads. Existing genome polishing tools usually introduce overcorrections and haplotype switch errors when correcting errors in genomes assembled from HiFi long reads. Here, we describe an upgraded genome polishing tool - NextPolish2, which can fix base errors remaining in those "highly accurate" genomes assembled from HiFi long reads without introducing excessive overcorrections and haplotype switch errors. We believe that NextPolish2 has a great significance to further improve the accuracy of telomere-to-telomere (T2T) genomes. NextPolish2 is freely available at https://github.com/Nextomics/NextPolish2.

Heterobifunctional small molecules to modulate RNA function

RNA has diverse cellular functionality, including regulating gene expression, protein translation, and cellular response to stimuli, due to its intricate structures. Over the past decade, small molecules have been discovered that target functional structures within cellular RNAs and modulate their function. Simple binding, however, is often insufficient, resulting in low or even no biological activity. To overcome this challenge, heterobifunctional compounds have been developed that can covalently bind to the RNA target, alter RNA sequence, or induce its cleavage. Herein, we review the recent progress in the field of RNA-targeted heterobifunctional compounds using representative case studies. We identify critical gaps and limitations and propose a strategic pathway for future developments of RNA-targeted molecules with augmented functionalities.

Comparative genomics of macaques and integrated insights into genetic variation and population history

The crab-eating macaques ( Macaca fascicularis ) and rhesus macaques ( M. mulatta ) are widely studied nonhuman primates in biomedical and evolutionary research. Despite their significance, the current understanding of the complex genomic structure in macaques and the differences between species requires substantial improvement. Here, we present a complete genome assembly of a crab-eating macaque and 20 haplotype-resolved macaque assemblies to investigate the complex regions and major genomic differences between species. Segmental duplication in macaques is ∼42% lower, while centromeres are ∼3.7 times longer than those in humans. The characterization of ∼2 Mbp fixed genetic variants and ∼240 Mbp complex loci highlights potential associations with metabolic differences between the two macaque species (e.g., CYP2C76 and EHBP1L1 ). Additionally, hundreds of alternative splicing differences show post-transcriptional regulation divergence between these two species (e.g., PNPO ). We also characterize 91 large-scale genomic differences between macaques and humans at a single-base-pair resolution and highlight their impact on gene regulation in primate evolution (e.g., FOLH1 and PIEZO2 ). Finally, population genetics recapitulates macaque speciation and selective sweeps, highlighting potential genetic basis of reproduction and tail phenotype differences (e.g., STAB1 , SEMA3F , and HOXD13 ). In summary, the integrated analysis of genetic variation and population genetics in macaques greatly enhances our comprehension of lineage-specific phenotypes, adaptation, and primate evolution, thereby improving their biomedical applications in human diseases.

A chromosome-scale fishing cat reference genome for the evaluation of potential germline risk variants

The fishing cat, Prionailurus viverrinus, faces a population decline, increasing the importance of maintaining healthy zoo populations. Unfortunately, zoo-managed individuals currently face a high prevalence of transitional cell carcinoma (TCC), a form of bladder cancer. To investigate the genetics of inherited diseases among captive fishing cats, we present a chromosome-scale assembly, generate the pedigree of the zoo-managed population, reaffirm the close genetic relationship with the Asian leopard cat (Prionailurus bengalensis), and identify 7.4 million single nucleotide variants (SNVs) and 23,432 structural variants (SVs) from whole genome sequencing (WGS) data of healthy and TCC cats. Only BRCA2 was found to have a high recurrent number of missense mutations in fishing cats diagnosed with TCC when compared to inherited human cancer risk variants. These new fishing cat genomic resources will aid conservation efforts to improve their genetic fitness and enhance the comparative study of feline genomes.

Structurally divergent and recurrently mutated regions of primate genomes

We sequenced and assembled using multiple long-read sequencing technologies the genomes of chimpanzee, bonobo, gorilla, orangutan, gibbon, macaque, owl monkey, and marmoset. We identified 1,338,997 lineage-specific fixed structural variants (SVs) disrupting 1,561 protein-coding genes and 136,932 regulatory elements, including the most complete set of human-specific fixed differences. We estimate that 819.47 Mbp or ∼27% of the genome has been affected by SVs across primate evolution. We identify 1,607 structurally divergent regions wherein recurrent structural variation contributes to creating SV hotspots where genes are recurrently lost (e.g., CARD, C4, and OLAH gene families) and additional lineage-specific genes are generated (e.g., CKAP2, VPS36, ACBD7, and NEK5 paralogs), becoming targets of rapid chromosomal diversification and positive selection (e.g., RGPD gene family). High-fidelity long-read sequencing has made these dynamic regions of the genome accessible for sequence-level analyses within and between primate species.

Long-term monitoring of ultratrace nucleic acids using tetrahedral nanostructure-based NgAgo on wearable microneedles

Real-time and continuous monitoring of nucleic acid biomarkers with wearable devices holds potential for personal health management, especially in the context of pandemic surveillance or intensive care unit disease. However, achieving high sensitivity and long-term stability remains challenging. Here, we report a tetrahedral nanostructure-based Natronobacterium gregoryi Argonaute (NgAgo) for long-term stable monitoring of ultratrace unamplified nucleic acids (cell-free DNAs and RNAs) in vivo for sepsis on wearable device. This integrated wireless wearable consists of a flexible circuit board, a microneedle biosensor, and a stretchable epidermis patch with enrichment capability. We comprehensively investigate the recognition mechanism of nucleic acids by NgAgo/guide DNA and signal transformation within the Debye distance. In vivo experiments demonstrate the suitability for real-time monitoring of cell-free DNA and RNA with a sensitivity of 0.3 fM up to 14 days. These results provide a strategy for highly sensitive molecular recognition in vivo and for on-body detection of nucleic acid.

Pangenomics: A new era in the field of neurodegenerative diseases

A pangenome is composed of all the genetic variability of a group of individuals, and its application to the study of neurodegenerative diseases may provide valuable insights into the underlying aspects of genetic heterogenetiy for these complex ailments, including gene expression, epigenetics, and translation mechanisms. Furthermore, a reference pangenome allows for the identification of previously undetected structural commonalities and differences among individuals, which may help in the diagnosis of a disease, support the prediction of what will happen over time (prognosis) and aid in developing novel treatments in the perspective of personalized medicine. Therefore, in the present review, the application of the pangenome concept to the study of neurodegenerative diseases will be discussed and analyzed for its potential to enable an improvement in diagnosis and prognosis for these illnesses, leading to the development of tailored treatments for individual patients from the knowledge of the genomic composition of a whole population.

Telomere-to-telomere assembly of cassava genome reveals the evolution of cassava and divergence of allelic expression

Cassava is a crucial crop that makes a significant contribution to ensuring human food security. However, high-quality telomere-to-telomere cassava genomes have not been available up to now, which has restricted the progress of haploid molecular breeding for cassava. In this study, we constructed two nearly complete haploid resolved genomes and an integrated, telomere-to-telomere gap-free reference genome of an excellent cassava variety, 'Xinxuan 048', thereby providing a new high-quality genomic resource. Furthermore, the evolutionary history of several species within the Euphorbiaceae family was revealed. Through comparative analysis of haploid genomes, it was found that two haploid genomes had extensive differences in linear structure, transcriptome features, and epigenetic characteristics. Genes located within the highly divergent regions and differentially expressed alleles are enriched in the functions of auxin response and the starch synthesis pathway. The high heterozygosity of cassava 'Xinxuan 048' leads to rapid trait segregation in the first selfed generation. This study provides a theoretical basis and genomic resource for molecular breeding of cassava haploids.

Chromosome-Level Genome Assembly of Protosalanx chinensis and Response to Air Exposure Stress

Protosalanx chinensis is a suitable particular species for genetic studies on nearly scaleless skin, transparency and high sensitivity to hypoxia stress. Here, we generated a high-quality chromosome-level de novo assembly of P. chinensis. The final de novo assembly yielded 379.47 Mb with 28 pseudo-chromosomes and a scaffold N50 length of 14.52 Mb. In total, 21,074 protein-coding genes were predicted. P. chinensis, Esox lucius and Hypomesus transpacificus had formed a clade, which diverged about 115.5 million years ago. In the air exposure stress experiment, we found that some genes play an essential role during P. chinensis hypoxia, such as bhlh, Cry1, Clock, Arntl and Rorb in the circadian rhythm pathway. These genomic data offer a crucial foundation for P. chinensis ecology and adaptation studies, as well as a deeper understanding of the response to air exposure stress.

Characterization of large-scale genomic differences in the first complete human genome

BACKGROUND

The first telomere-to-telomere (T2T) human genome assembly (T2T-CHM13) release is a milestone in human genomics. The T2T-CHM13 genome assembly extends our understanding of telomeres, centromeres, segmental duplication, and other complex regions. The current human genome reference (GRCh38) has been widely used in various human genomic studies. However, the large-scale genomic differences between these two important genome assemblies are not characterized in detail yet.

RESULTS

Here, in addition to the previously reported "non-syntenic" regions, we find 67 additional large-scale discrepant regions and precisely categorize them into four structural types with a newly developed website tool called SynPlotter. The discrepant regions (~ 21.6 Mbp) excluding telomeric and centromeric regions are highly structurally polymorphic in humans, where the deletions or duplications are likely associated with various human diseases, such as immune and neurodevelopmental disorders. The analyses of a newly identified discrepant region-the KLRC gene cluster-show that the depletion of KLRC2 by a single-deletion event is associated with natural killer cell differentiation in ~ 20% of humans. Meanwhile, the rapid amino acid replacements observed within KLRC3 are probably a result of natural selection in primate evolution.

CONCLUSION

Our study provides a foundation for understanding the large-scale structural genomic differences between the two crucial human reference genomes, and is thereby important for future human genomics studies.

Human Pangenomics: Promises and Challenges of a Distributed Genomic Reference

A pangenome is a collection of the common and unique genomes that are present in a given species. It combines the genetic information of all the genomes sampled, resulting in a large and diverse range of genetic material. Pangenomic analysis offers several advantages compared to traditional genomic research. For example, a pangenome is not bound by the physical constraints of a single genome, so it can capture more genetic variability. Thanks to the introduction of the concept of pangenome, it is possible to use exceedingly detailed sequence data to study the evolutionary history of two different species, or how populations within a species differ genetically. In the wake of the Human Pangenome Project, this review aims at discussing the advantages of the pangenome around human genetic variation, which are then framed around how pangenomic data can inform population genetics, phylogenetics, and public health policy by providing insights into the genetic basis of diseases or determining personalized treatments, targeting the specific genetic profile of an individual. Moreover, technical limitations, ethical concerns, and legal considerations are discussed.

A review of the pangenome: how it affects our understanding of genomic variation, selection and breeding in domestic animals?

As large-scale genomic studies have progressed, it has been revealed that a single reference genome pattern cannot represent genetic diversity at the species level. While domestic animals tend to have complex routes of origin and migration, suggesting a possible omission of some population-specific sequences in the current reference genome. Conversely, the pangenome is a collection of all DNA sequences of a species that contains sequences shared by all individuals (core genome) and is also able to display sequence information unique to each individual (variable genome). The progress of pangenome research in humans, plants and domestic animals has proved that the missing genetic components and the identification of large structural variants (SVs) can be explored through pangenomic studies. Many individual specific sequences have been shown to be related to biological adaptability, phenotype and important economic traits. The maturity of technologies and methods such as third-generation sequencing, Telomere-to-telomere genomes, graphic genomes, and reference-free assembly will further promote the development of pangenome. In the future, pangenome combined with long-read data and multi-omics will help to resolve large SVs and their relationship with the main economic traits of interest in domesticated animals, providing better insights into animal domestication, evolution and breeding. In this review, we mainly discuss how pangenome analysis reveals genetic variations in domestic animals (sheep, cattle, pigs, chickens) and their impacts on phenotypes and how this can contribute to the understanding of species diversity. Additionally, we also go through potential issues and the future perspectives of pangenome research in livestock and poultry.

Structurally divergent and recurrently mutated regions of primate genomes

To better understand the pattern of primate genome structural variation, we sequenced and assembled using multiple long-read sequencing technologies the genomes of eight nonhuman primate species, including New World monkeys (owl monkey and marmoset), Old World monkey (macaque), Asian apes (orangutan and gibbon), and African ape lineages (gorilla, bonobo, and chimpanzee). Compared to the human genome, we identified 1,338,997 lineage-specific fixed structural variants (SVs) disrupting 1,561 protein-coding genes and 136,932 regulatory elements, including the most complete set of human-specific fixed differences. Across 50 million years of primate evolution, we estimate that 819.47 Mbp or ~27% of the genome has been affected by SVs based on analysis of these primate lineages. We identify 1,607 structurally divergent regions (SDRs) wherein recurrent structural variation contributes to creating SV hotspots where genes are recurrently lost (CARDs, ABCD7, OLAH) and new lineage-specific genes are generated (e.g., CKAP2, NEK5) and have become targets of rapid chromosomal diversification and positive selection (e.g., RGPDs). High-fidelity long-read sequencing has made these dynamic regions of the genome accessible for sequence-level analyses within and between primate species for the first time.

USAT: a bioinformatic toolkit to facilitate interpretation and comparative visualization of tandem repeat sequences

Background

Tandem repeats (TR), highly variable genomic variants, are widely used in individual identification, disease diagnostics, and evolutionary studies. The recent advances in sequencing technologies and bioinformatic tools facilitate calling TR haplotypes genome widely. Both length-based and sequence-based TR alleles are used in different applications. However, sequence-based TR alleles could provide the highest precision in characterizing TR haplotypes. The need to identify the differences at the single nucleotide level between or among TR haplotypes with an easy-use bioinformatic tool is essential.

Results

In this study, we developed a Universal STR Allele Toolkit (USAT) for TR haplotype analysis, which takes TR haplotype output from existing tools to perform allele size conversion, sequence comparison of haplotypes, figure plotting, comparison for allele distribution, and interactive visualization. An exemplary application of USAT for analysis of the CODIS core STR loci for DNA forensics with benchmarking human individuals demonstrated the capabilities of USAT. USAT has user-friendly graphic interfaces and runs fast in major computing operating systems with parallel computing enabled.

Conclusion

USAT is a user-friendly bioinformatics software for interpretation, visualization, and comparisons of TRs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-05021-1.