SEARCH · PubMed Health

Results for “chromosome collinearity”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

18 recordsLinked to original sources

PubMed2025

Genome-Wide Analysis of the AT-Hook Gene Family in Malus sieversii and Functional Characterization of MsAHL13.

AT-hook motif nuclear-localized (AHL) proteins are pivotal in plant growth, development, and stress responses. Nevertheless, there is limited research on AHL proteins in Malus sieversii. Our study identified 25 AHL genes from the M. sieversii genome, named MsAHL1-MsAHL25. The encoded protein sequences had lengths ranging from 195 to 554 amino acids, molecular weights from 19.17 to 58.53 kDa, and isoelectric points from 4.67 to 10.09. Chromosomal mapping revealed that these 25 genes were unevenly distributed across 10 chromosomes. Collinearity analysis of AHL genes in M. sieversii implied that gene loss might have occurred during its evolution. The phylogenetic tree classified the AHL proteins of M. sieversii into two subfamilies, showing a close relationship with multiple proteins of M. domestica. Promoter analysis indicated that the AHL genes in M. sieversii harbored numerous stress- and hormone-responsive elements, suggesting their potential role in various stress responses. qRT-PCR analysis of six representative MsAHLs under biotic and abiotic stresses demonstrated that the expression of MsAHL13, MsAHL15, and MsAHL17 was significantly upregulated under salt, drought, and cold stresses, while MsAHL01 expression was inhibited under low-temperature stress. All six MsAHLs were induced by the pathogen Valsa mali. Subcellular localization analysis of the specifically expressed protein MsAHL13 showed its nuclear location. Furthermore, luciferase and yeast two-hybrid assays confirmed the in vitro physical interaction between the MsAHL13 and MsMYB1 proteins. This research offers an important theoretical basis for further exploration of the functional mechanisms of this gene family in responding to environmental stresses.

Malus sieversii↗

PubMed2026

Chromosome-Scale Genome of Zoonotic Eyeworm Thelazia callipaeda from China.

Thelazia callipaeda is a vector-borne zoonotic eyeworm infecting companion animals, wildlife, and humans, but chromosome-scale genomic resources from Chinese clinical material remain limited. We generated a genome supported by Pacific Biosciences (PacBio) high-fidelity (HiFi) sequencing and high-throughput chromosome conformation capture (Hi-C) from 100 adult worms recovered from naturally infected dogs in Beijing and compared its chromosome-scale organization with Portuguese assembly GCA_965194785.1. The final assembly spans 119.53 megabases (Mb) and comprises 115 top-level sequences, including four pseudomolecules totaling 91.26 Mb (76.34%) and 111 unanchored sequences. Genome-mode Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis recovered 98.5% complete chromadorean orthologues, and the representative 11,788-protein gene set recovered 92.6%. Sequence-level alignment resolved Chinese chromosomes 1-4 (chr1-chr4) to Portuguese chr1, chrX, chr3, and chr2, respectively, with retained alignments covering 95.9-99.2% of each Chinese pseudomolecule and estimated sequence identities of 99.75-99.91%. Strong chromosome-scale collinearity was accompanied by localized reverse-collinear regions, including 0.243 Mb and 0.115 Mb intervals on chr2-chrX and chr3-chr3. The anchored sequences contained 96.7% of predicted genes and were substantially more gene-dense than the unanchored sequences. These results establish a clinically sourced Chinese chromosome-scale reference and provide a validated framework for future individual-worm, population-genomic, structural-variation, and comparative genomic studies of this parasite.

Hi-C↗

PubMed2026

A Chromosome-Level Genome Assembly of the Potato Leafhopper Empoasca fabae (Hemiptera: Cicadellidae).

The potato leafhopper, Empoasca fabae (Harris, 1841), is a highly polyphagous, migratory insect pest of eastern North America that feeds on more than 200 herbaceous and woody plant species, causing substantial losses to forage and field crops. Despite its agricultural and ecological importance, no genome has been available for this species. Here, we present the first chromosome-level genome assembly of E. fabae, generated from Oxford Nanopore long reads, Illumina short reads, and Omni-C proximity-ligation data. The final assembly spans 908 Mb across 132 scaffolds, with 99.8% of the assembly captured in ten chromosome-length scaffolds (nine autosomes and an X chromosome) with a scaffold N50 of 96.2 Mb. The assembly is highly complete, recovering 92.9% of conserved hemipteran single-copy orthologs from protein annotations, and is composed of 47.6% repetitive sequence, dominated by long terminal repeat retrotransposons and unclassified elements. Read-depth comparison between male and female individuals supports assignment of a single sex-linked chromosome, consistent with an XO sex determination system. BRAKER3 gene annotation predicted 31,406 protein-coding genes after retaining the longest isoform per locus. Comparative genome analysis of the two closest related Typhlocybinae species with genomes available, Matsumurasca onukii and Hebata decipiens, revealed extensive chromosome-scale collinearity while defining a shared core gene repertoire. This reference genome provides a foundation for comparative and population genomic studies and for investigating genetic traits in this economically important crop pest species.

Animals↗

PubMed2026

Genome-wide characterization of BraABCB transporters reveals their potential roles in hormone responses in Brassica rapa var. parachinensis.

Thirty-six BraABCB genes were identified in Brassica rapa var. parachinensis; expression and interaction analyses suggest BraABCB27 and BraABCB28 as hormone-responsive candidates associated with BRI1-related proteins. ABCB transporters are ATP-dependent membrane proteins that mediate the transport of diverse substrates, including phytohormones, and play important roles in plant development and environmental adaptation. Previous studies in Arabidopsis have shown that several ABCB proteins participate in phytohormone transport, including auxin and brassinosteroid transport, whereas the functions of their homologs in Brassica rapa var. parachinensis remain poorly understood. In this study, 36 BraABCB genes were identified and classified into four phylogenetic groups. Conserved domain analysis showed that BraABCB proteins contain typical nucleotide-binding domains and transmembrane domains. Chromosomal distribution, collinearity, and Ka/Ks analyses suggested that the BraABCB family is evolutionarily conserved and mainly subject to purifying constraints. Promoter analysis and RT-qPCR assays of selected Group IV BraABCB genes revealed diverse expression patterns and responses to drought, high temperature, brassinolide, and indole-3-acetic acid treatments. Subcellular localization assays showed that selected Group IV BraABCB proteins exhibited predominant plasma membrane localization. Notably, BraABCB27 and BraABCB28, two close AtABCB1/AtABCB19-related homologs, showed detectable physical associations with BRI1-related proteins in BiFC and split-ubiquitin yeast two-hybrid assays. Together, these results provide a genome-wide characterization of the BraABCB gene family and identify BraABCB27 and BraABCB28 as candidate genes for future studies of their possible associations with BR-related membrane processes and hormone-regulated growth responses in B. rapa var. parachinensis.

Plant Growth Regulators↗

PubMed2026

Contrasting redox-related physiological responses associated with HaGATA23 and HaGATA36 during Orobanche cumana parasitism in sunflower (Helianthus annuus L.).

Helianthus annuus is an economically important Asteraceae species used for seed oil production and ornamental purposes, but its production is seriously affected by the root-parasitic plant Orobanche cumana. GATA transcription factors are zinc-finger DNA-binding regulators involved in plant development and stress adaptation. However, the molecular characteristics of GATA transcription factors in Helianthus annuus and their contribution to Helianthus annuus -Orobanche cumana interaction remain poorly understood. Here, 36 HaGATA members were retrieved from the Helianthus annuus genome and classified into four phylogenetic clades. Chromosomal placement, collinearity, gene structure, motif composition, and promoter elements varied among the 36 HaGATA members, indicating evolutionary conservation coupled with functional diversification. Expression analysis and RT-qPCR analyses revealed differential expression patterns among HaGATA genes under O. cumana stress, with HaGATA23 markedly downregulated and HaGATA36 strongly upregulated. Overexpression of HaGATA23 was associated with increased malondialdehyde (MDA) accumulation and unfavorable changes in antioxidant enzyme activities, whereas its silencing showed the opposite physiological tendency. In contrast, overexpression of HaGATA36 reduced malondialdehyde accumulation, increased peroxidase (POD), catalase (CAT), and superoxide dismutase (SOD) activities, while its silencing showed the reverse tendency. These results indicate that HaGATA23 and HaGATA36 are candidate genes associated with contrasting redox-related physiological responses during O. cumana stress. This work provides evidence that GATA transcription factors are associated with redox-related physiological responses in sunflower under O. cumana treatment and identifies HaGATA23 and HaGATA36 as functionally divergent candidate genes for further validation.

Helianthus↗

PubMed2025

Polyploidy-mediated variations in glutamate receptor proteins linked to Fusarium wilt resistance in upland cotton.

Cotton production in the US faces a serious threat from Fusarium oxysporum f. sp. vasinfectum race 4 (FOV4), a soil-borne fungus causing Fusarium wilt by infecting the roots and vascular system of susceptible cotton, leading to rapid wilting and death. Here, we investigate genetic mechanisms of resistance to FOV4 in the highly resistant upland cotton genotype "U1" using an early-generation segregating biparental population ("U1" × "CSX8308") with comprehensive genomic resources. Reference-grade genomic assemblies of the parents revealed minor structural variations between "U1" haplotypes, a high degree of collinearity at chromosome synteny and micro-synteny levels, and significant divergence from "CSX8308" with 8.9 million SNPs. QTL analysis identified significant markers on chromosomes D03 and A02 linked to reduced Fusarium wilt severity. Within these regions, two glutamate-receptor-like (GLR) genes showed structural variation and overlapped between translocated segments on A02 and D03, suggesting a rare but important reinforcing effect of parallel evolution between susceptible and resistant genotypes. Transcriptome profiles of "U1" under FOV4 infection reveal activation of calcium-binding proteins and transcription factors regulating plant hormones (ethylene, abscisic acid, jasmonic acid, and salicylic acid), along with enzymes involved in cell wall remodeling and phytoalexin production. Advancing cotton improvement depends on incorporating durable genetic disease resistance into high-yielding, high-quality cultivars.

Fusarium↗

PubMed2025

Identification of BoRR gene family in cauliflower: roles in curd development and salt tolerance.

BACKGROUND: Cauliflower, as an important vegetable crop, the research on its curd formation mechanism and stress-responsive gene networks is of great significance for improving its quality, yield and abiotic stress tolerance. The response regulator (RR) gene family plays a crucial role in the regulation of various life processes of many organisms. In this research, a comprehensive analysis of the BoRR gene family in cauliflower was carried out. RESULTS: A Total of 57 BoRR genes were identified in cauliflower and classified into seven subtypes (type A/B-I/B-II/B-IV/C/B-PRR/Clock PRR) based on sequence homology. Chromosomal mapping showed even distribution across genomes, while physicochemical analysis revealed diverse protein properties (134-915 amino acids, pI 4.51-9.19) with predominant nuclear localization. Structural analyses found all BoRR proteins contain REC-type domains, with subtype-specific features: type A has REC_typeA_ARR, type B harbors REC_typeB_ARR domains, and Clock PRR shows circadian-related psREC_RR domains. Exon numbers range from 2 to 10, with type A BoRR genes having shorter CDS lengths. Collinearity analysis identified 28 pairs of gene duplicates (26 inter-chromosomal). Comparative analysis showed 133 collinear pairs with Brassica napus, 96 with Brassica. rapa, and only 1 with monocots specie (rice and maize). Promoter analysis identified hormone-responsive motifs (ABRE, TGACG), development-related elements (ARE), and stress-responsive sequences (e.g., MBS for drought tolerance) in the promoters of BoRR genes. GO enrichment linked BoRR genes to phosphorelay signaling, cytokinin/ethylene response, and developmental processes like meristem maintenance. Expression profiling during curd development showed type A genes (BoRR23/27/34/38/45) up-regulated in vegetative-reproductive transition, BoRR3/6/12/32/54 in curd enlargement, and several genes like BoRR49 in flower bud differentiation. Salt stress (1.5% NaCl) induced transient expression in 8 of 9 selected BoRR genes at day 1 after treatment. qRT-PCR validated their roles in developmental regulation and salt tolerance. CONCLUSION: This study provides valuable insights into the BoRR gene family in cauliflower, laying a foundation for further understanding its genetic mechanisms and potentially guiding efforts to enhance curd quality and salt tolerance in cauliflower.

Salt Tolerance↗

PubMed2026

The Key Trichoderma-Induced Gene Encoding a DUF568 Domain-Containing Protein Mediates Defense Responses in Wheat.

Genes encoding DUF568 domain-containing proteins participate in plant stress adaptation. To elucidate the functional role of DUF568 domain-containing genes in Trichoderma-induced wheat defense responses against wheat Fusarium crown rot, we performed a genome-wide identification and characterization of the TaDUF568 gene family in hexaploid wheat (Triticum aestivum L.). In this study, a total of 33 TaDUF568 family genes were systematically identified and characterized at the genome-wide level, exhibiting uneven chromosomal distribution and diverse physicochemical properties. Phylogenetic, structural, and collinearity analyses revealed conserved family characteristics among monocot species. Segmental duplication was verified as the primary driver of gene family expansion. Expression profiling revealed divergent tissue-specific expression patterns among TaDUF568 family members, among which TaDUF568.18 was strongly induced by Trichoderma M2. Subcellular localization assays confirmed that TaDUF568.18 is a plasma membrane-localized protein. Functional validation via stable transgenes demonstrated that overexpression of TaDUF568.18 restricted lesion expansion, improved agronomic traits, and enhanced disease resistance. This study is the first to characterize the wheat DUF568 family and confirm that TaDUF568.18 (annotated as TaAIR12) acts as a positive regulator of Trichoderma-mediated wheat defense, providing a valuable gene resource for wheat disease-resistance breeding.

DUF568↗

PubMed2026

Comparative analysis of olfactory receptor repertoires reveals evolutionary dynamics and high-altitude adaptation in Schizopygopsis younghusbandi based on the chromosome-level genomes.

The olfactory receptor (OR) gene represent a significant multigene family in vertebrates, forming the core molecular basis of olfactory perception and playing a crucial role in the environmental adaptation of species. High-altitude ecosystems represent extreme habitats characterized by specific abiotic stresses, including low oxygen levels, low temperatures, and intense ultraviolet radiation. These environments also exhibit low aquatic biodiversity and a limited variety of odor molecules, factors that have influenced the adaptive evolution of the sensory systems in endemic species. However, the genetic mechanisms underlying olfactory adaptation in high-altitude freshwater fish remained inadequately understood. In this study, we performed comparative genomics analyses to reveal the evolutionary processes underlying the adaptive and functional evolution of OR genes in S. younghusbandi, a cyprinid fish endemic to the Qinghai-Xizang Plateau. The results indicated that, compared to their low-altitude relatives, S. younghusbandi possessed a significantly smaller number of OR genes, with only 98 genes, which revealed the contraction of the gene family. Phylogenetic analysis revealed that the OR genes of cyprinid fish could be categorized into two major lineages: type I and type II. The η and δ families, which perceive water-soluble odors, in S. younghusbandi underwent significant and specific expansion, while the ε family was completely absent. This pattern reflected adaptive changes in olfactory recognition to accommodate the simplified odor spectrum of high-altitude water bodies. Chromosomal localization analysis demonstrated that OR genes were clustered, and collinearity analysis confirmed the presence of conserved genomic fragments among species. Selection pressure analysis revealed that the Ka/Ks values of all homologous gene pairs were less than 1, indicating that the OR genes of S. younghusbandi underwent strong purifying selection as a group to preserve core olfactory function. A few genes exhibited relaxed selection characteristics, which may have facilitated the fine-tuning of adaptability to high-altitude environments. In conclusion, this study elucidated the evolutionary dynamics and adaptive characteristics of the OR gene in S. younghusbandi, offering a new perspective on the molecular mechanisms underlying olfactory adaptation at high altitudes and enriching the research on sensory evolution in vertebrates.

Schizopygopsis younghusbandi↗

PubMed2026

Genomic and evolutionary basis of parthenogenesis in a disease-vector tick species.

Haemaphysalis longicornis is an important tick species and pathogen vector characterized by the co-circulation of triploid parthenogenetic and diploid bisexual strains. However, the evolutionary basis of parthenogenesis in this species is unclear. Here we report reference-quality, haplotype-resolved genome assemblies of the parthenogenetic strain and two reference-quality genomes of the bisexual strains. Comparative genomic analysis revealed high collinearity between the parthenogenetic and bisexual genomes, with a stable chromosomal architecture maintained among the three haplotypes of the parthenogenetic strain. The parthenogenetic H. longicornis genome exhibited a major expansion in cell cycle-related gene families, including the inhibitor of apoptosis protein (IAP) family, but was characterized by a contraction in other gene families. Population resequencing of 179 individuals revealed two distinct subpopulations, with chromosome 7 harbouring high genetic differentiation and several candidate genes probably associated with parthenogenesis. Functional experiments showed that knockdown of the BIRC5 gene, a member of the IAP family, suppressed oviposition in both strains, with the parthenogenetic strain exhibiting milder adverse effects probably due to a stronger transcriptional response. Overall, our results reveal the genomic and evolutionary features associated with polyploid parthenogenesis in H. longicornis.

Animals↗

PubMed2026

Genome-wide characterization of the bZIP gene family in Rattus norvegicus and expression profiling analysis during brain development.

BACKGROUND: The brown rat (Rattus norvegicus) serves as a cornerstone model organism in biomedical research, particularly for understanding physiological homeostasis and stress responses. The basic leucine zipper (bZIP) transcription factor family is a pivotal regulatory network involved in growth, organogenesis, and neurodevelopment. Despite its importance, a systematic characterization of the bZIP gene family in rats has remained elusive. RESULTS: In this study, we performed a genome-wide identification of 61 RnbZIP genes, which were categorized into 10 distinct subfamilies based on phylogenetic relationships and chromosomal localization. Structural analysis revealed conserved motif arrangements within subfamilies, while collinearity analysis identified significant gene duplication events-predominantly tandem and segmental duplications-that have driven the evolutionary expansion of the RnbZIP family. Quantitative analysis showed that members within the same subfamily shared 45%-92% sequence similarity (calculated using the BLOSUM62 scoring matrix), and all duplicated gene pairs underwent strong purifying selection (Ka/Ks < 1). Comparative genomics across seven rodent species further underscored the evolutionary conservation and divergence of these factors. Expression profiling across diverse organs and brain developmental stages indicated that RnbZIP genes exhibit high tissue specificity. Notably, 10 candidate genes, including RnbZIP01, RnbZIP02, and RnbZIP08, demonstrated dynamic expression patterns during brain maturation, suggesting their essential roles in neurodevelopmental processes. CONCLUSIONS: Our findings provide a comprehensive structural and evolutionary framework for the RnbZIP gene family, highlighting their potential regulatory functions in rat organogenesis and brain development. This study establishes a valuable resource for further functional characterization of specific bZIP members in mammalian neurological systems.

Animals↗

PubMed2025

Genome-Wide Identification of the TIFY Family in Cannabis sativa L. and Its Potential Functional Analysis in Response to Alkaline Stress and in Cannabinoid Metabolism.

TIFY transcription factors play crucial regulatory roles in secondary metabolism and stress response. However, the expression patterns of the Cannabis sativa L. TIFY gene family under alkali stress, their involvement in cannabinoid metabolism, and their underlying genetic evolutionary mechanisms remain largely unexplored. In this study, we used bioinformatics approaches to conduct genome-wide identification and functional characterization of the C. sativa TIFY gene family. Fourteen TIFY genes were identified and mapped onto seven chromosomes. These genes were classified into four subfamilies: TIFY, JAZ, ZML, and PPD, with the JAZ subfamily further subdivided into five distinct branches. Collinearity analysis suggested that gene duplication events contributed to the expansion of the TIFY gene family in C. sativa. Weighted gene coexpression network analysis (WGCNA) revealed that CsJAZ2, CsJAZ3, and CsJAZ6 participated in the cannabinoid regulatory network. Cis-element analysis indicated that the promoter regions of TIFY genes were enriched in hormone- and stress-responsive elements. Furthermore, transcriptome and RT-qPCR analyses were conducted to examine gene expression patterns under alkaline stress (the RNA employed in RT-qPCR was extracted from the apical leaves of samples subjected to short-duration alkaline stress treatment). The results showed that CsJAZ5 and CsJAZ6 were downregulated, whereas CsPPD1, CsTIFY1, and CsZML1 were upregulated in response to alkali stress. In summary, CsJAZ5, CsPPD1, and CsTIFY1 may serve as candidate genes for the development of alkali-tolerant cultivars, while CsJAZ2 and CsJAZ3 may be valuable targets for enhancing cannabinoid production. This study provides important molecular insights and a theoretical basis for future research on the evolutionary dynamics and functional roles of TIFY transcription factors, particularly in stress adaptation and cannabinoid metabolism.

Cannabis↗

PubMed2026

Whole-genome characterization of seven multidrug-resistant Neisseria gonorrhoeae isolates from a single tertiary center in Beijing.

BACKGROUND: To characterize the whole-genome features of Neisseria gonorrhoeae clinical isolates collected from a tertiary medical institution in Beijing, with a focus on the genomic basis of ceftriaxone non-susceptibility and multidrug resistance. METHODS: Clinical isolates were collected from April 2023 to November 2024. Of 14 collected isolates, seven were successfully subcultured after revival and included in subsequent analyses. Minimum inhibitory concentrations (MICs) were determined by the Etest method. Whole-genome data were obtained using a combination of second- and third-generation sequencing technologies. The isolates were combined with global and Chinese reference datasets to construct a core-genome single-nucleotide polymorphism (core-SNP) phylogenetic tree. Chromosomal resistance-associated mutations and plasmid characteristics were subsequently analyzed. RESULTS: The seven isolates displayed genomic diversity at the whole-genome level. Four isolates (8087, 8423, 8461, and 8801) carried penA 60.001 and belonged to distinct sequence types, including ST7365, ST8123, and ST7367. One additional isolate (8726) carried penA 273.001; both alleles encode PBP2 proteins sharing the core substitutions A311V, I312M, V316T, and T483S. All five isolates were non-susceptible to ceftriaxone (MIC 0.25-0.5 mg/L). Ceftriaxone non-susceptibility was associated with the co-occurrence of mutations at core penA positions and additional mutations in porB and ponA, with an mtrR mutation present in one isolate. Plasmid collinearity analysis revealed that several multidrug-resistant isolates simultaneously harbored an intact conjugative plasmid and an African-type resistance plasmid carrying bla TEM-1. CONCLUSION: The multidrug-resistant phenotype of Neisseria gonorrhoeae results from the co-existence of chromosomal multi-locus mutations and resistance plasmids. The penA 60.001 isolates in this study did not originate from a single source. This allele appeared in multiple local clonal lineages. This pattern is consistent with horizontal gene transfer of this resistance determinant into multiple endemic lineages.

Neisseria gonorrhoeae↗

PubMed2025

Chromosome-level Genome Assembly of the Halophytic Turfgrass Zoysia macrostachya.

Zoysia macrostachya Franch. & Sav. is a halophytic perennial turfgrass in the Poaceae family, commonly found in the coastal regions of Korea, Japan, and East Asia. Z. macrostachya thrives in high-salinity environments, making it an excellent model for studying abiotic stress resilience. In this study, we present a chromosome-level genome assembly of Z. macrostachya, constructed using Oxford Nanopore long reads, Illumina short reads, and Omni-C sequencing data. The assembly spans 329.78 Mb across 20 chromosomes, with a scaffold N50 of 19.24 Mb, and includes complete telomeric sequences at both ends. The assembly showed 97.8% complete BUSCOs, indicating high genome completeness. Repeat element and gene annotation identified 44.03% of the genome as repetitive elements and 33,474 protein-coding genes. The gene annotation showed 97.1% complete BUSCOs and 86.92% functionally characterized genes. Macrosynteny analysis highlighted highly collinear relationships with related species, providing a foundational understanding of the Z. macrostachya genomic structure. This high-quality genome serves as a valuable resource for advancing salinity tolerance research and improving the genetic diversity of Zoysia species.

Genome, Plant↗

PubMed2026

Genome-Wide Characterization of the Apple HD-Zip IV Gene Family and Functional Validation of MdHDZIV3 Under PEG-Induced Osmotic Stress.

The homeodomain-leucine zipper IV (HD-Zip IV) transcription factor subfamily plays essential roles in epidermal development, cuticle formation, lipid metabolism, and environmental adaptation in plants. Despite its biological importance, the HD-Zip IV family has not been systematically characterized in apple (Malus domestica). Here, we identified 17 apple HD-Zip IV genes and named them MdHDZIV1-MdHDZIV17 based on their locations on the chromosomes. The 17 genes showed a nonuniform distribution on eight chromosomes, while the occurrence of both tandem and segmental duplications indicated that family expansion involved more than one duplication mechanism. All MdHDZIV proteins contained the conserved HD, LZ, START, and SAD domains but lacked the MEKHLA domain, consistent with typical HD-Zip IV structural features. Phylogenetic analysis classified MdHDZIV proteins into five groups together with HD-Zip IV members from Arabidopsis thaliana and rice, indicating evolutionary conservation of this subfamily. Collinearity and Ka/Ks analyses revealed that duplicated MdHDZIV gene pairs were mainly subjected to purifying selection. Promoter scanning revealed diverse cis-regulatory motifs associated with hormonal signaling, environmental stress, light response, and epidermal regulation, including ABRE, ARE, W-box, MYC, G-box, and L1-box motifs. Integration of transcriptomic profiling with qRT-PCR validation revealed pronounced tissue-dependent differences in the expression of MdHDZIV genes in leaf, fruit skin, and branch bark. Under PEG6000-induced osmotic stress and NaCl-induced salt stress, 10 candidate MdHDZIV genes displayed gene-specific and stress type-specific expression patterns, with MdHDZIV3 showing strong induction under PEG6000 treatment. Functional validation in apple calli showed that MdHDZIV3 overexpression enhanced PEG tolerance, increased fresh weight, elevated SOD and POD activities, and reduced MDA accumulation under osmotic stress. These findings provide a genome-wide framework for understanding the apple HD-Zip IV gene family.

abiotic stress↗

PubMed2026

Reference-Guided Chromosome-Scale Genome Assembly With Insights on Population Genomics of the Atlantic Goliath Grouper (Epinephelus itajara), Islas del Rosario, Colombia.

Epinephelus itajara, commonly known as the Atlantic Goliath grouper, is the largest species among the western North Atlantic groupers and is critically endangered. This species plays a crucial ecological, cultural, and economic role and has been the focus of captive breeding efforts at the Oceanario of the Rosario Islands, Colombia. However, despite its ecological and conservation importance, genomic resources and population genomic data for E. itajara remain scarce, particularly in the Colombian Caribbean. This study presents a reference-guided chromosome-scale genome assembly and an analysis of the population genomic structure of E. itajara using PacBio HiFi sequencing and Illumina technologies. The assembled genome has a total size of 1.12 Gb, with a contig N50 of 42.69 Mb and a scaffold N50 of 46.30 Mb. A total of 22,692 protein-coding genes were identified after masking 46% of the genome, which consists of repetitive elements. Comparative genomic analyses revealed a high degree of collinearity with closely related Epinephelus species and identified E. lanceolatus as the closest relative, supporting recent divergence and conserved genome architecture within the genus. Additionally, a population genomics analysis was conducted using 7706 high-quality SNPs to assess the genomic structure of captive populations. The results revealed four distinct genomic lineages, with moderate genetic differentiation among the sampled individuals. In the Colombian Caribbean, two unique lineages were identified, associated with the localities of Bahía Cispatá and Bahía Barbacoas, suggesting possible geographic isolation. These genomic resources provide valuable tools and new opportunities to better understand the genomic diversity, evolutionary history, and reproductive mechanisms of E. itajara. Moreover, they serve as a foundation for conservation strategies, including selective breeding programs aimed at increasing genomic diversity in captive populations and guiding restoration efforts in its natural habitat.

Epinephelus itajara↗

PubMed2026

Prophage landscapes in clinical MRSA: safety profiling and discovery of Lys81, a broad-spectrum bacteriolytic enzyme.

INTRODUCTION: Methicillin-resistant Staphylococcus aureus (MRSA) poses a significant threat to global healthcare, requiring novel therapeutic strategies. Prophages, latent phage genomes integrated into bacterial chromosomes, are important resources for antimicrobial development due to their genomic stability and genetic engineering potential. METHODS: In this study, we performed genomewide sequencing on 329 MRSA isolates to predict prophage sequences, followed by analyses of these prophages-including examinations of virulence genes, antibiotic resistance genes, homologous proteins of pathogenic MRSA phages, and functional predictions of these homologous proteins-to evaluate their safety and value as genetic engineering scaffolds and to screen for novel broadspectrum bacteriolytic enzymes. RESULTS: Our data indicate that 85.7% (282/329) of strains carried complete prophage sequences; 64 strains lacked virulence factors or genes, meeting the core criteria for safe vectors. Resistance screening found only 6 prophages carried msrA, confirming the biosafety of the remaining strains. A significant correlation existed between prophage virulence gene capacity and genomic structure (R2 = 0.99986684, p = 3.64e-69). High-virulence clusters (>10 factors) showed high structural similarity; 10 characteristic sequences linked to S. aureus phages and their prevalence patterns were identified via conserved motif analysis. Collinearity analysis with reference to virulent MRSA phages and 3D structural predictions of orthologous proteins identified two lysozymes and a host-recognition device. Notably, Lys81, an N-acetylmuramoyl-L-alanine amidase ortholog, was prioritized and characterized as a broad-spectrum lytic enzyme. Our data show Lys81 has key properties: (1) Broad-spectrum antibacterial activity, lysing 52.3% (23/44) of clinical S. aureus strains and cross-acting against Gram-positive bacteria such as Pseudomonas aeruginosa and Listeria; (2) Excellent environmental adaptability, maintaining activity at pH 5.0 and 0°C, with 25 mM Na+ and Ca2 + enhancing function; (3) Potent biofilm clearance, achieving 83% MRSA biofilm reduction at 50 μg/mL; and (4) Favorable in vivo safety/efficacy, eradicating MRSA infections in lung organoid models with minimal cytotoxicity. DISCUSSION: This study establishes a theoretical foundation for the clinical translation of MRSA prophages, positioning Lys81 as a novel candidate for treating drug-resistant bacterial infections.

Lys81↗

PubMed2026

Pangenomic analyses in the cultivated grapevine confirm high genomic collinearity and extensive dispensable gene content likely involved in adaptation.

Pangenomes have now been developed for several horticultural crops, yet the extent to which genome diversity in sequence and organization contribute to plant adaptation and major agronomic traits remains poorly understood. Here, we assembled the genomes of 9 cultivated grapevine varieties and compared the genomes of 15 cultivated grapevine varieties for variation in gene and TE content. We found that genomic collinearity is highly conserved among varieties. We still observed substantial variation across genomes. Notably, we identified across varieties 55,662 orthologous genes, of which 55.3% appears to be dispensable. Dispensable genes are enriched for functions related to adaptation to biotic and abiotic constraints, suggesting that they may play a role in adaptation. Comparing our results with a recently published study, we found substantial differences with ∼12.6% of the genes we classified as core genes being classified as dispensable genes in this other study. We then constructed a pangenome graph and used it to performed genome-wide association studies for 3 important traits in grapevine production, which allowed us to include large structural variants as markers in the analyses. We identified 32 loci that we did not detect when we used the PN40024 genome as a reference, 20 of which are newly reported associations. Overall, our results indicates that despite recent advances in characterizing plant pangenomes, current gene classification into core and dispensable gene categories should be taken with caution. They also highlight the value of incorporating structural variants into GWAS, to better characterize the genetic architecture of agronomic traits.

Vitis↗