Biomedical subjects

Heng Li

Publications and source records attributed to Heng Li.

At least 19 recordsLinked to original sources

PubMed2026

Silent carriage of tigecycline- and carbapenem-resistant Klebsiella quasipneumoniae co-harbouring tetX(4) and blaNDM-1 in healthy individuals in China.

OBJECTIVES: To investigate the occurrence and genomic characteristics of tigecycline- and carbapenem-resistant Klebsiella quasipneumoniae isolated from healthy individuals in China. METHODS: During a nationwide screening programme, faecal samples from healthy community individuals were cultured for carbapenem-resistant Enterobacterales. Antimicrobial susceptibility testing, whole-genome sequencing and conjugation assays were performed for two K. quasipneumoniae isolates co-harbouring tet(X4) and blaNDM-1. RESULTS: Two K. quasipneumoniae strains carrying tet(X4) and blaNDM-1 were recovered from healthy individuals without recent hospitalisation or antibiotic exposure. Both isolates showed resistance to tigecycline (>8 μg/mL) and carbapenems (≥4 μg/mL). Genomic analysis demonstrated close clonal relatedness between the isolates (ST6460-1LV, KL151). The blaNDM-1 gene was located on an IncX3 plasmid associated with mobile genetic elements, whereas tet(X4) was carried by an IncX1 plasmid. Conjugation assays confirmed successful transfer of both resistance genes, which frequently co-transferred into recipient strains. CONCLUSIONS: To our knowledge, this is the first report of K. quasipneumoniae co-harbouring tet(X4) and blaNDM-1 in healthy individuals in China. The findings indicate that healthy community populations may serve as a hidden reservoir for last-resort antimicrobial resistance genes and underscore the importance of community-based surveillance within a One Health framework.

Klebsiella quasipneumoniae↗

PubMed2026

Large-scale genomic analysis places Chinese CC398 as a persistent human-associated MSSA lineage apart from the dominant global LA-MRSA clade.

Staphylococcus aureus clonal complex (CC)398 has emerged as a dominant livestock-associated methicillin-resistant S. aureus (LA-MRSA) lineage worldwide; however, its evolutionary trajectory and regional diversification remain incompletely understood. We developed a core-genome multilocus sequence typing (cgMLST) scheme with hierarchical clustering and applied it to over 30,000 S. aureus genomes, revealing frequent cross-border transmission of CC398. Subsequent time-calibrated phylogenetic analysis placed the most recent common ancestor at 1942 (95% CI: 1939-1945), with the human-to-livestock host jump around 1969 (95% CI: 1968-1972). Chinese CC398 exhibits a distinct trajectory: unlike the LA-MRSA lineages dominating Europe and North America, Chinese isolates are predominantly human-associated methicillin-susceptible S. aureus (HA-MSSA), forming unique East Asia-specific phylogroups (SAP1, SAP2, and AP1-AP3), with distinct resistance and virulence profiles. The LA lineage remains limited in China, with multinational mixed clusters emerging only after 2019. Analysis of global transmission networks revealed a significant correlation between LA-CC398 spread and international trade in fresh swine products, while no such correlation was observed for the human-associated lineage. Beyond the established lineage markers tet(M) and scn, our analysis identified additional differentially distributed genes, including cadC-a chromosomal cadmium resistance regulator-as a novel HA-lineage-enriched gene whose functional role in host adaptation remains to be determined. This study reveals that CC398 followed fundamentally different evolutionary paths in China versus Western countries, challenging a one-size-fits-all model of its dissemination.IMPORTANCEThis study illustrates how large-scale microbial genomics can resolve the evolutionary origins and regional diversification of bacterial pathogens. By applying a novel cgMLST scheme to over 30,000 S. aureus genomes, we show that CC398 followed fundamentally different evolutionary paths in China versus Western countries-challenging the prevailing model of uniform global dissemination-and that livestock-associated MRSA expansion is closely linked to international trade in fresh pork products. These findings highlight the need for integrated surveillance across human, animal, and trade interfaces to anticipate the emergence and spread of zoonotic pathogens.

Staphylococcus aureus↗

PubMed2026

Improving long-read somatic structural variant calling with pangenome and de novo personal genome assembly.

Accurate detection of mosaic and somatic structural variants (SVs) provides early diagnostic and therapeutic evidence for cancers. While long-read whole-genome sequencing leads to more accurate SV detection than short read sequencing, existing long-read SV callers only look at alignment against a single reference genome and are susceptible to systematic false discovery caused by germline differences between the individual genome and the reference genome. Here we develop a new SV filtering method that jointly considers the alignment against a pangenome and the de novo assembly of the germline genome. It dramatically reduces false positive mosaic and somatic SVs in cancer cell lines with little loss in sensitivity for existing long read SV callers. Our study highlights the essential need for pangenome or personal genome assembly to integrate SV calls for both SV discoveries and clinical diagnostics.

Journal Article↗

PubMed2026

FuFiHLA: a tool for full-field HLA typing from long-read data.

MOTIVATION: Allele typing for Human Leukocyte Antigen (HLA) genes has many important clinical applications. Popular short-read typing can only accurately distinguish alleles at the coding sequence level, which potentially limit our understanding of the effect of variants in non-coding region. Long read data has been proved to be useful in typing HLA alleles in full resolution, but only a few tools are publicly available and with significant limitations in practical application. RESULTS: We developed FuFiHLA, a lightweight open-source software, to type HLA alleles. Currently it supports typing alleles of six HLA genes (HLA-A, HLA-B, HLA-C, HLA-DRB1, HLA-DQA1, and HLA-DQB1) from long reads. Evaluation using 233 PacBio HiFi WGS samples from HPRC shows that FuFiHLA achieves 99.6% accuracy in the full field allele typing and QV as 51.8 for consensus allele sequence construction. Additional testing on four Nanopore R10 reads demonstrates slightly reduced accuracy in the fourth field. AVAILABILITY: FuFiHLA is available at https://github.com/jingqing-hu/FuFiHLA under MIT License.

Humans↗

PubMed2025

Improving spliced alignment by modeling splice sites with deep learning.

MOTIVATION: Spliced alignment refers to the alignment of messenger RNA (mRNA) or protein sequences to eukaryotic genomes. It plays a critical role in gene annotation and the study of gene functions. Accurate spliced alignment demands sophisticated modeling of splice sites, but current aligners use simple models, which may affect their accuracy given dissimilar sequences. RESULTS: We implemented minisplice to learn splice signals with a one-dimensional convolutional neural network (1D-CNN) and trained a model with 7,026 parameters for vertebrate and insect genomes. It captures conserved splice signals across phyla and reveals GC-rich introns specific to mammals and birds. We used this model to estimate the empirical splicing probability for every GT and AG in genomes, and modified minimap2 and miniprot to leverage pre-computed splicing probability during alignment. Evaluation on human long-read RNA-seq data and cross-species protein datasets showed our method greatly improves the junction accuracy especially for noisy long RNA-seq reads and proteins of distant homology. AVAILABILITY AND IMPLEMENTATION: https://github.com/lh3/minisplice.

Journal Article↗

PubMed2025

High-resolution metagenome assembly for modern long reads with myloasm.

Long-read metagenome assembly promises complete genomic recovery from microbiomes. However, the complexity of metagenomes poses challenges. We present myloasm, a metagenome assembler for PacBio HiFi and Oxford Nanopore Technologies (ONT) R10.4 long reads. Myloasm uses polymorphic k-mers to construct a high-resolution string graph and then leverages differential abundance for graph simplification. On real-world ONT metagenomes, myloasm assembled three times more complete circular contigs than the next-best assembler. Myloasm can make ONT and HiFi comparable for assembly: for a jointly sequenced gut metagenome, myloasm with ONT assembled more complete circular genomes than any assembler with HiFi. Myloasm recovers previously inaccessible within-species diversity; we recovered six complete Prevotella copri single-contig genomes from a gut metagenome and eight complete TM7 (Saccharibacteria) contigs with > 93% similarity from an oral metagenome. With this improved resolution, we resolved two 98% similar ermF antibiotic resistance genes spreading through distinct strain-specific mobile genetic elements in a human gut.

Journal Article↗

PubMed2025

Defining and cataloging variants in pangenome graphs.

Structural variation causes some human haplotypes to align poorly with the linear reference genome, leading to 'reference bias'. A pangenome reference graph could ameliorate this bias by relating a sample to multiple reference assemblies. However, this approach requires a new definition of a 'genetic variant.' We introduce a definition of pangenome variants and a method, pantree, to identify them. Our approach involves a pangenome reference tree which includes all nodes (sequences) of the pangenome graph, but only a subset of its edges; non-reference edges are variant edges. Our variants are biallelic and have well-defined positions. Analyzing the Minigraph-Cactus draft human pangenome reference graph, we identified 29.6 million genetic variants. Most variants (99.2%) are small, and most small variants (73.9%) are SNPs. 3.5 million variants (11.7%) have a reference allele which is not on GRCh38; these variants are difficult to detect without a pangenome reference, or with existing pangenome-based approaches. They tend to be embedded within tangled, multiallelic regions. We analyze two medically relevant regions, around the HLA-A and RHD genes, identifying thousands of small variants embedded within several large insertions, deletions, and inversions. We release an open-source software tool together with a VCF variant catalogue.

Journal Article↗

PubMed2025

Pitfalls of bacterial pan-genome analysis approaches: a case study of Mycobacterium tuberculosis and two less clonal bacterial species.

SUMMARY: Pan-genome analysis is a fundamental tool for studying bacterial genome evolution; however, the variety in methods used to define and measure the pan-genome poses challenges to the interpretation and reliability of results. Using Mycobacterium tuberculosis, a clonally evolving bacterium with a small accessory genome, as a model system, we systematically evaluated sources of variability in pan-genome estimates. Our analysis revealed that differences in assembly type (short-read versus hybrid), annotation pipeline, and pan-genome software, significantly impact predictions of core and accessory genome size. Extending our analysis to two additional bacterial species, Escherichia coli and Staphylococcus aureus, we observed consistent tool-dependent biases but species-specific patterns in pan-genome variability. Our findings highlight the importance of integrating nucleotide- and protein-level analyses to improve the reliability and reproducibility of pan-genome studies across diverse bacterial populations. AVAILABILITY AND IMPLEMENTATION: Panqc is freely available under an MIT license at https://github.com/maxgmarin/panqc.

Genome, Bacterial↗

PubMed2024

Neotelomeres and telomere-spanning chromosomal arm fusions in cancer genomes revealed by long-read sequencing.

Alterations in the structure and location of telomeres are pivotal in cancer genome evolution. Here, we applied both long-read and short-read genome sequencing to assess telomere repeat-containing structures in cancers and cancer cell lines. Using long-read genome sequences that span telomeric repeats, we defined four types of telomere repeat variations in cancer cells: neotelomeres where telomere addition heals chromosome breaks, chromosomal arm fusions spanning telomere repeats, fusions of neotelomeres, and peri-centromeric fusions with adjoined telomere and centromere repeats. These results provide a framework for the systematic study of telomeric repeats in cancer genomes, which could serve as a model for understanding the somatic evolution of other repetitive genomic elements.

Humans↗

PubMed2023

Neotelomeres and Telomere-Spanning Chromosomal Arm Fusions in Cancer Genomes Revealed by Long-Read Sequencing.

Alterations in the structure and location of telomeres are key events in cancer genome evolution. However, previous genomic approaches, unable to span long telomeric repeat arrays, could not characterize the nature of these alterations. Here, we applied both long-read and short-read genome sequencing to assess telomere repeat-containing structures in cancers and cancer cell lines. Using long-read genome sequences that span telomeric repeat arrays, we defined four types of telomere repeat variations in cancer cells: neotelomeres where telomere addition heals chromosome breaks, chromosomal arm fusions spanning telomere repeats, fusions of neotelomeres, and peri-centromeric fusions with adjoined telomere and centromere repeats. Analysis of lung adenocarcinoma genome sequences identified somatic neotelomere and telomere-spanning fusion alterations. These results provide a framework for systematic study of telomeric repeat arrays in cancer genomes, that could serve as a model for understanding the somatic evolution of other repetitive genomic elements.

Telomere↗

PubMed2006

Snap: an integrated SNP annotation platform.

Snap (Single Nucleotide Polymorphism Annotation Platform) is a server designed to comprehensively analyze single genes and relationships between genes basing on SNPs in the human genome. The aim of the platform is to facilitate the study of SNP finding and analysis within the framework of medical research. Using a user-friendly web interface, genes can be searched by name, description, position, SNP ID or clone name. Several public databases are integrated, including gene information from Ensembl, protein features from Uniprot/SWISS-PROT, Pfam and DAS-CBS. Gene relationships are fetched from BIND, MINT, KEGG and are integrated with ortholog data from TreeFam to extend the current interaction networks. Integrated tools for primer-design and mis-splicing analysis have been developed to facilitate experimental analysis of individual genes with focus on their variation. Snap is available at http://snap.humgen.au.dk/ and at http://snap.genomics.org.cn/.

Databases, Nucleic Acid↗

PubMed2006

Demonstration of phycobilisome mobility by the time- and space-correlated fluorescence imaging of a cyanobacterial cell.

The cell-wide mobility of PBSs was confirmed by synchronously monitoring the fluorescence recovery after photobleaching (FRAP) and the fluorescence loss in photobleaching (FLIP). On the other hand, a fluorescence recovery was still observed even if PBSs were immobile (PBSs fixed on the membranes by betaine and isolated PBSs fixed on the agar plate) or PBS mobility was unobservable (cell wholly bleached). Furthermore, it was proved that some artificial factors were involved not only in FRAP but also in FLIP, including renaturation of the reversibly denatured proteins, laser scanning-induced fluorescence loss and photo-damage to the cell. With consideration of the fast renaturation component in fluorescence recovery, the diffusion coefficient was estimated to be tenfold smaller than that without the component. Moreover, it was observed that the fluorescence intensity on the bleached area was always lower than that on the non-bleached area, even after 20 min, while it should be equal if PBSs were mobile freely. Based on the increasing proportion of the PBSs anti-washed to Triton X-100 (1%) with prolonged laser irradiation to the cells locked in light state 1 by PBQ, it was concluded that some PBSs became immobile due to photo-linking to PSII.

Cyanobacteria↗

PubMed2006

PigGIS: Pig Genomic Informatics System.

Pig Genomic Information System (PigGIS) is a web-based depository of pig (Sus scrofa) genomic learning mainly engineered for biomedical research to locate pig genes from their human homologs and position single nucleotide polymorphisms (SNPs) in different pig populations. It utilizes a variety of sequence data, including whole genome shotgun (WGS) reads and expressed sequence tags (ESTs), and achieves a successful mapping solution to the low-coverage genome problem. With the data presently available, we have identified a total of 15 700 pig consensus sequences covering 18.5 Mb of the homologous human exons. We have also recovered 18 700 SNPs and 20 800 unique 60mer oligonucleotide probes for future pig genome analyses. PigGIS can be freely accessed via the web at http://www.piggis.org/ and http://pig.genomics.org.cn/.

Animals↗

PubMed2006

The state transition mechanism - simply depending on light-on and -off in Spirulina platensis.

The state transition in cyanobacteria is a long-discussed topic of how the photosynthetic machine regulates the excitation energy distribution in balance between the two photosystems. In the current work, whether the state transition is realized by "mobile phycobilisome (PBS)" or "energy spillover" has been clearly answered by monitoring the spectral responses of the intact cells of the cyanobacterium Spirulina platensis. Firstly, light-induced state transition depends completely on a movement of PBSs toward PSI or PSII while the redox-induced one on not only the "mobile PBS" but also an "energy spillover". Secondly, the "energy spillover" is triggered by dissociation of PSI trimers into the monomers which specially occurs under a case from light to dark, while the PSI monomers will re-aggregate into the trimers under a case from dark to light, i.e., the PSI oligomerization is reversibly regulated by light switch on and off. Thirdly, PSI oligomerization is regulated by the local H(+) concentration on the cytosol side of the thylakoid membranes, which in turn is regulated by light switch on and off. Fourthly, PSI oligomerization change is the only mechanism for the "energy spillover". Thus, it can be concluded that the "mobile PBS" is a common rule for light-induced state transition while the "energy spillover" is only a special case when dark condition is involved.

Biochemistry↗

PubMed2006

Phylogeny and classification of Paris (Melanthiaceae) inferred from DNA sequence data.

BACKGROUND AND AIMS: Paris (Melanthiaceae) is a temperate genus of about 24 perennial herbaceous species distributed from Europe to eastern Asia. The delimitation of the genus and its subdivisions are unresolved questions in the taxonomy of Paris. The objective of this study is to test the generic and infrageneric circumscription of Paris with DNA sequence data. METHODS: Phylogenetic analysis of 21 species of Paris based on nuclear ITS and plastid psbA-trnH and trnL-trnF DNA sequence data, alone and in combination, was employed to assess previous classifications. KEY RESULTS: Paris is monophyletic in all analyses. Neither of the two traditionally recognized subgenera (Paris and Daiswa) are monophyletic. Sections Axiparis, Kinugasa, Paris and Thibeticae are monophyletic in only some of the analyses. Species of sections Dunnianae, Fargesianae and Marmoratae are consistently intercalated among species of section Euthyra in all analyses. Strong discordance between nuclear and plastid lineages is detected. CONCLUSIONS: The data support the classification of Paris as a single genus rather than as three genera (Daiswa, Kinugasa and Paris sensu stricto). They provide justification for the transfer of section Axiparis from subgenus Paris to subgenus Daiswa and for the combination of sections Dunnianae, Fargesianae and Marmoratae into section Euthyra. The nuclear-plastid discordance is interpreted as the result of interspecific hybridization among sympatric species.

Biological Evolution↗

PubMed2006

Intercontinental biogeography of subfamily Orontioideae (Symplocarpus, Lysichiton, and Orontium) of Araceae in Eastern Asia and North America.

Symplocarpus, Lysichiton, and Orontium (Orontioideae) are three of the few north temperate genera of the primarily tropical Araceae. Symplocarpus is disjunctly distributed in eastern Asia (3 spp.) and eastern North America (1 sp.); Lysichiton has an intercontinental discontinuous distribution in eastern Asia (1 sp.) and northwestern North America (1 sp.); and the monotypic Orontium is restricted to eastern North America. Phylogenetic analysis of the trnL-F and ndhF sequences supports (1) the monophyly of both Symplocarpus and Lysichiton, (2) the sister-group relationship of Symplocarpus and Lysichiton, and (3) the clade of Orontium, Symplocarpus, and Lysichiton. Although Symplocarpus shows a much wider disjunction than Lysichiton, the estimated divergence time of the former [4.49+/-1.69 or 6.88+/-4.18 million years ago (mya)] was similar to that of the latter (4.02+/-1.60 or 7.18+/-4.33 mya) based on the penalized likelihood and the Bayesian dating methods, respectively. Eastern Asia was suggested to be the ancestral area of the Symplocarpus-Lysichiton clade based on the dispersal-vicariance analysis. Our biogeographic results support independent migrations of Symplocarpus and Lysichiton across the Bering land bridge in the late Tertiary (Pliocene/late Miocene). Fossil evidence suggests Orontioideae dated back to the late Cretaceous in the temperate Northern Hemisphere (72 mya). The relative rate test shows similar substitution rates of the trnL-F sequences between the proto and the true aroids, although the latter has substantially higher species diversity. The proto Araceae perhaps suffered from a higher rate of extinction in the temperate zone associated with periods of climatic cooling in the Tertiary.

Animals↗

PubMed2006

Inhibitory effects of L- and T-type calcium antagonists on contractions of human detrusor muscle.

The inhibitory and relaxant effects of the L-type calcium antagonists nifedipine, nimodipine, verapamil and diltiazem, and of the T-type calcium antagonist mibefradil, on contractions of isolated human detrusor muscle were investigated. The tissue was obtained from 10 patients undergoing cystectomy due to bladder cancer. Effects of the calcium antagonists at different concentrations on the concentration-response curves for carbachol were investigated. Furthermore, concentration-relaxation curves were performed using potassium-precontracted muscle strips. All L-type calcium antagonists suppressed the mean concentration-response curve of carbachol significantly at a concentration of 10(-6) M. Mibefradil up to 10(-5) M did not significantly suppress it. Nifedipine significantly reduced the carbachol-induced maximum contraction to 75% and 44%, verapamil to 75% and 67% of the appropriate control value at concentrations of 10(-7) and 10(-6) M, respectively. Diltiazem reduced it insignificantly to 96% and 71% at the above-mentioned concentrations. The concentration-relaxation experiments revealed following pD2-values and maximum relaxations of nifedipine, nimodipine, verapamil and diltiazem, respectively: 6.23, 6.37, 5.66, 5.81 and 85%, 83%, 82%, 90%. Maximum relaxations and pD2-values were not significantly different from each other. The lowest concentration, for which a significant effect compared to control in Student;s t-test was found, amounted to 10(-10) M, 10(-9) M, 10(-7) M, 10(-6.5) M and 10(-4) M for nimodipine, nifedipine, diltiazem, verapamil and mibefradil, respectively. L-type calcium antagonists are very potent relaxant agents of the human detrusor muscle in vitro.

Adult↗

PubMed2006

Inhibitory effects of various L-type and T-type calcium antagonists on electrically generated, potassium-induced and carbachol-induced contractions of porcine detrusor muscle.

The inhibitory effects of different calcium antagonists on contractions of isolated porcine detrusor muscle were investigated. Suppression of the maximum potassium-induced contraction and electrically generated contractions by nifedipine, verapamil and diltiazem were investigated. Furthermore, concentration-response curves of carbachol after pretreatment with the L-type antagonists nifedipine, verapamil, diltiazem, nimodipine and the T-type antagonist mibefradil at different concentrations were performed. Nifedipine significantly reduced the potassium-induced maximum contraction to 89, 60, 21, 8 and 4% (10(-9)-10(-5) M). Verapamil and diltiazem significantly reduced it to 64, 30 and 5% (10(-7)-10(-5) M) or 79, 27, 7 and 1% (10(-7)-10(-4) M), respectively. Nifedipine, verapamil and diltiazem significantly reduced the electrically generated contraction to 55, 36, 34 and 25% (10(-7)-10(-4) M), 71, 32 and 2% (10(-6)-10(-4) M), 96, 78, 38 and 5% (10(-7)-10(-4) M), respectively. pD2 values of nifedipine, verapamil and diltiazem amounted to 7.07, 5.56 and 5.40 and differed significantly. After pretreatment with nifedipine at 10(-6) M, the concentration-response curve of carbachol was nearly suppressed. The effects of nimodipine, verapamil and diltiazem were smaller. Mibefradil caused only at 10(-5) M a significant reduction. All investigated L-type calcium antagonists were strong inhibitors of the examined contractions. Nifedipine showed the biggest inhibitory effect.

Animals↗