High-throughput Chromatin Conformation Capture Research Articles

Transcriptional Hubs Within Cliques in Ensemble Hi-C Chromatin Interaction Networks.

Chromatin conformation capture technologies permit the study of chromatin spatial organization on a genome-wide scale at a variety of resolutions. Despite the increasing precision and resolution of high-throughput chromatin conformation capture (Hi-C) methods, it remains challenging to conclusively link transcriptional activity to spatial organizational phenomena. We have developed a clique-based approach for analyzing Hi-C data that helps identify chromosomal hotspots that feature considerable enrichment of chromatin annotations for transcriptional start sites and, building on previously published work, show that these chromosomal hotspots are not only significantly enriched in RNA polymerase II binding sites as identified by the ENCODE project, but also identify a noticeable increase in FANTOM5 and GTEx transcription within our identified cliques across a variety of tissue types. From the obtained data, we surmise that our cliques are a suitable method for identifying transcription factories in Hi-C data, and outline further extensions to the method that may make it useful for locating regions of increased transcriptional activity in datasets where in-depth expression or polymerase data may not be available.

First telomere-to-telomere gapless assembly of the rice blast fungus Pyricularia oryzae

Rice blast caused by Pyricularia oryzae (syn., Magnaporthe oryzae) was one of the most destructive diseases of rice throughout the world. Genome assembly was fundamental to genetic variation identification and critically impacted the understanding of its ability to overcome host resistance. Here, we report a gapless genome assembly of rice blast fungus P. oryzae strain P131 using PacBio, Illumina and high throughput chromatin conformation capture (Hi-C) sequencing data. This assembly contained seven complete chromosomes (43,237,743 bp) and a circular mitochondrial genome (34,866 bp). Approximately 14.31% of this assembly carried repeat sequences, significantly greater than its previous assembled version. This assembly had a 99.9% complement in BUSCO evaluation. A total of 14,982 genes protein-coding genes were predicted. In summary, we assembled the first telomere-to-telomere gapless genome of P. oryzae, which would be a valuable genome resource for future research on the genome evolution and host adaptation.

Modification of the Hi-C Technology for Molecular Genetic Analysis of Formalin-Fixed Paraffin-Embedded Sections of Tumor Tissues.

Molecular genetic analysis of tumor tissues is the most important step towards understanding the mechanisms of cancer development; it is also necessary for the choice of targeted therapy. The Hi-C (high-throughput chromatin conformation capture) technology can be used to detect various types of genomic variants, including balanced chromosomal rearrangements, such as inversions and translocations. We propose a modification of the Hi-C method for the analysis of chromatin contacts in formalin-fixed paraffin-embedded (FFPE) sections of tumor tissues. The developed protocol allows to generate high-quality Hi-C data and detect all types of chromosomal rearrangements. We have analyzed various databases to compile a comprehensive list of translocations that hold clinical importance for the targeted therapy selection. The practical value of molecular genetic testing is its ability to influence the treatment strategies and to provide prognostic insights. Detecting specific chromosomal rearrangements can guide the choice of the targeted therapies, which is a critical aspect of personalized medicine in oncology.

BL-Hi-C reveals the 3D genome structure of Brassica crops with high sensitivity.

High-throughput Chromatin Conformation Capture (Hi-C) technologies can be used to investigate the three-dimensional genomic structure of plants. However, the practical utility of these technologies is impeded by significant background noise, hindering their capability in detecting fine 3D genomic structures. In this study, we optimized the Bridge Linker Hi-C technology (BL-Hi-C) to comprehensively investigate the 3D chromatin landscape of Brassica rapa and Brassica oleracea. The Bouquet configuration of both B. rapa and B. oleracea was elucidated through the construction of a 3D genome simulation. The optimized BL-Hi-C exhibited lower background noise compared to conventional Hi-C methods. Taking this advantage, we used BL-Hi-C to identify FLC gene loops in Arabidopsis, B. rapa, and B. oleracea. We observed that gene loops of FLC2 exhibited conservation across Arabidopsis, B. rapa, and B. oleracea. While gene loops of syntenic FLCs exhibited conservation across B. rapa and B. oleracea, variations in gene loops were evident among multiple paralogs FLCs within the same species. Collectively, our findings highlight the high sensitivity of optimized BL-Hi-C as a powerful tool for investigating the fine 3D genomic organization.

Chromosome-level genome assembly of the Stoliczka’s Asian trident bat (Aselliscus stoliczkanus)

Stoliczka’s Asian trident bat (Aselliscus stoliczkanus) is a small-bodied species and very sensitive to climate change. Here, we presented a chromosome-level genome assembly of A. stoliczkanus by combining Illumina sequencing, Nanopore sequencing and high-throughput chromatin conformation capture (Hi-C) sequencing technology. The genome assembly was 2.18 Gb in size with 98.26% of the genome sequences anchored onto 14 autosomes and two sex chromosomes (X and Y). The quality of the genome assembly is very high with a contig and scaffold N50 of 72.98 and 162 Mb, respectively, Benchmarking Universal Single-Copy Orthologs (BUSCO) score of 96.6%, and the consensus quality value (QV) of 47.44. A total of 20,567 genes were predicted and 98.8% of these genes were functionally annotated. Syntenic blocks between A. stoliczkanus and Homo sapiens, together with previous comparative cytogenetic studies, provide valuable foundations for further comparative genomic and cytogenetic studies in mammals. The reference-quality genome of A. stoliczkanus contributes an important resource for conservative genomics and landscape genomics in predicting adaptation and vulnerability to climate change.

Hi-C sequencing unravels dynamic three-dimensional chromatin interactions in muntjac lineage: insights from chromosome fusions in Fea's muntjac genome.

Eukaryotes have varying numbers and structures of characteristic chromosomes across lineages or species. The evolutionary trajectory of species may have been affected by spontaneous genome rearrangements. Chromosome fusion drastically alters karyotypes. However, the mechanisms and consequences of chromosome fusions, particularly in muntjac species, are poorly understood. Recent research-based advancements in three-dimensional (3D) genomics, particularly high-throughput chromatin conformation capture (Hi-C) sequencing, have allowed for the identification of chromosome fusions and provided mechanistic insights into three muntjac species: Muntiacus muntjak, M. reevesi, and M. crinifrons. This study aimed to uncover potential genome rearrangement patterns in the threatened species Fea's muntjac (Muntiacus feae), which have not been previously examined for such characteristics. Deep Hi-C sequencing (31.42 × coverage) was performed to reveal the 3D chromatin architecture of the Fea's muntjac genome. Patterns of repeated chromosome fusions that were potentially mediated by high-abundance transposable elements were identified. Comparative Hi-C maps demonstrated linkage homology between the sex chromosomes in Fea's muntjac and autosomes in M. reevesi, indicating that fusions may have played a crucial role in the evolution of the sex chromosomes of the lineage. The species-level dynamics of topologically associated domains (TADs) suggest that TAD organization could be altered by differential chromosome interactions owing to repeated chromosome fusions. However, research on the effect of TADs on muntjac genome evolution is insufficient. This study generated Hi-C data for the Fea's muntjac, providing a genomic resource for future investigations of the evolutionary patterns of chromatin conformation at the chromosomal level.

Accessible gene borders establish a core structural unit for chromatin architecture in Arabidopsis.

Three-dimensional (3D) chromatin structure is linked to transcriptional regulation in multicellular eukaryotes including plants. Taking advantage of high-resolution Hi-C (high-throughput chromatin conformation capture), we detected a small structural unit with 3D chromatin architecture in the Arabidopsis genome, which lacks topologically associating domains, and also in the genomes of tomato, maize, and Marchantia polymorpha. The 3D folding domain unit was usually established around an individual gene and was dependent on chromatin accessibility at the transcription start site (TSS) and transcription end site (TES). We also observed larger contact domains containing two or more neighboring genes, which were dependent on accessible border regions. Binding of transcription factors to accessible TSS/TES regions formed these gene domains. We successfully simulated these Hi-C contact maps via computational modeling using chromatin accessibility as input. Our results demonstrate that gene domains establish basic 3D chromatin architecture units that likely contribute to higher-order 3D genome folding in plants.

SnapHiC-D: a computational pipeline to identify differential chromatin contacts from single-cell Hi-C data.

Single-cell high-throughput chromatin conformation capture technologies (scHi-C) has been used to map chromatin spatial organization in complex tissues. However, computational tools to detect differential chromatin contacts (DCCs) from scHi-C datasets in development and through disease pathogenesis are still lacking. Here, we present SnapHiC-D, a computational pipeline to identify DCCs between two scHi-C datasets. Compared to methods designed for bulk Hi-C data, SnapHiC-D detects DCCs with high sensitivity and accuracy. We used SnapHiC-D to identify cell-type-specific chromatin contacts at 10 Kb resolution in mouse hippocampal and human prefrontal cortical tissues, demonstrating that DCCs detected in the hippocampal and cortical cell types are generally associated with cell-type-specific gene expression patterns and epigenomic features. SnapHiC-D is freely available at https://github.com/HuMingLab/SnapHiC-D.

HiBrowser: an interactive and dynamic browser for synchronous Hi-C data visualization.

With the development of chromosome conformation capture technology, the genome-wide investigation of higher-order chromatin structure by using high-throughput chromatin conformation capture (Hi-C) technology is emerging as an important component for understanding the mechanism of gene regulation. Considering genetic and epigenetic differences are typically used to explore the pathological reasons on the chromosome and gene level, visualizing multi-omics data and performing an intuitive analysis by using an interactive browser become a powerful and welcomed way. In this paper, we develop an effective sequence and chromatin interaction data display browser called HiBrowser for visualizing and analyzing Hi-C data and their associated genetic and epigenetic annotations. The advantages of HiBrowser are flexible multi-omics navigation, novel multidimensional synchronization comparisons and dynamic interaction system. In particular, HiBrowser first provides an out of the box web service and allows flexible and dynamic reconstruction of custom annotation tracks on demand during running. In order to conveniently and intuitively analyze the similarities and differences among multiple samples, such as visual comparisons of normal and tumor tissue samples, and pan genomes of multiple (consanguineous) species, HiBrowser develops a clone mode to synchronously display the genome coordinate positions or the same regions of multiple samples on the same page of visualization. HiBrowser also supports a pluralistic and precise search on correlation data of distal cis-regulatory elements and navigation to any region on Hi-C heatmap of interest according to the searching results. HiBrowser is a no-build tool, and could be easily deployed in local server. The source code is available at https://github.com/lyotvincent/HiBrowser.

3D genome structural variations play important roles in regulating seed oil content of Brassica napus

Dissecting the complex regulatory mechanism of seed oil content (SOC) is one of the main research goals in Brassica napus. Increasing evidence suggests that genome architecture is linked to multiple biological functions. However, the effect of genome architecture on SOC regulation remains unclear. Here, we used high-throughput chromatin conformation capture to characterize differences in the three-dimensional (3D) landscape of genome architecture of seeds from two B. napus lines, N53-2 (with high SOC) and Ken-C8 (with low SOC). Bioinformatics analysis demonstrated that differentially accessible regions and differentially expressed genes between N53-2 and Ken-C8 were preferentially enriched in regions with quantitative trait loci (QTLs)/associated genomic regions (AGRs) for SOC. A multi-omics analysis demonstrated that expression of SOC-related genes was tightly correlated with genome structural variations in QTLs/AGRs of B. napus. The candidate gene BnaA09g48250D, which showed structural variation in a QTL/AGR on chrA09, was identified by fine-mapping of a KN double-haploid population derived from hybridization of N53-2 and Ken-C8. Overexpression and knockout of BnaA09g48250D led to significant increases and decreases in SOC, respectively, in the transgenic lines. Taken together, our results reveal the 3D genome architecture of B. napus seeds and the roles of genome structural variations in SOC regulation, enriching our understanding of the molecular mechanisms of SOC regulation from the perspective of spatial chromatin structure.

Spatio-temporal transcriptome dynamics coordinate rapid transition of core crop functions in 'lactating' pigeon.

Pigeons (Columba livia) are among a select few avian species that have developed a specialized reproductive mode wherein the parents produce a 'milk' in their crop to feed newborn squabs. Nonetheless, the transcriptomic dynamics and role in the rapid transition of core crop functions during 'lactation' remain largely unexplored. Here, we generated a de novo pigeon genome assembly to construct a high resolution spatio-temporal transcriptomic landscape of the crop epithelium across the entire breeding stage. This multi-omics analysis identified a set of 'lactation'-related genes involved in lipid and protein metabolism, which contribute to the rapid functional transitions in the crop. Analysis of in situ high-throughput chromatin conformation capture (Hi-C) sequencing revealed extensive reorganization of promoter-enhancer interactions linked to the dynamic expression of these 'lactation'-related genes between stages. Moreover, their expression is spatially localized in specific epithelial layers, and can be correlated with phenotypic changes in the crop. These results illustrate the preferential de novo synthesis of 'milk' lipids and proteins in the crop, and provides candidate enhancer loci for further investigation of the regulatory elements controlling pigeon 'lactation'.

Abstract A054: Accurate prediction of cohesin-mediated 3D genome organization using 2D chromatin features

Abstract The three-dimensional (3D) genome organization influences diverse nuclear processes. Chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) and High-throughput chromatin conformation capture (Hi-C) are powerful methods used to study the 3D genome organization. However, these two experiments are costly, time-consuming, require tens to hundreds of millions of cells, and challenging to optimize and analyze. Predicting ChIA-PET/Hi-C data using cheaper ChIP-Seq data and other easily obtainable features is a useful alternative. Also, it is well-established that the cohesin protein complex is a key determinant of 3D genome organization. Considering this, we present Chromatin Interaction Predictor (ChIPr), a suite of regression models based on deep neural networks (DNN), random forest, and gradient boosting. ChIPr predicts cohesin-mediated chromatin interaction strength between any two given loci in the genome. It was trained mainly using two top-tiered ENCODE cell lines GM12878 and K562 and tested on other cell lines (H1 and HepG2). Comprehensive tests showed that ChIPr predictions correlated well with the original ChIA-PET data at the peak-level resolution and at bin sizes 25kb and 5kb. In addition, ChIPr accurately captured most of the cell-type-dependent loops identified by both, ChIA-PET and Hi-C data, respectively. Extensive feature testing highlighted genomic distance and RAD21 (a cohesin component) ChIP-Seq signals as the most important inputs contributing to ChIPr, in determining the chromatin interaction strength. A standard ChIPr model requires three experimental inputs: ChIP-Seq signals for RAD21, H3K27ac (enhancer/active chromatin mark), and H3K27me3 (inactive chromatin mark). While a reduced model requires a single experiment input: ChIP-Seq signals for RAD21 and performs equally well. Additionally, the use of ChIPr was further extended to predict the contact maps for several prostate cancer (PCa) cell lines whose ChIA-PET data were not available on the ENCODE portal, for instance, RWPE1 and VCaP. We generated the necessary ChIP-seq data for these cell lines to use as input to the model to further examine the changes in the 3D interactions of PCa driver genes and also their isolated neighbourhoods and decipher useful biological insights about the gene regulatory networks involved. Integrative analysis revealed novel insights into the role of CTCF motif, its orientation, and CTCF binding on the prevalence and strength of cohesin-mediated chromatin interactions (~50% - 80% of the RAD21 interactions were enriched with CTCF, < 15% of the interactions had no CTCF ChIP-seq binding in both of the two peaks). Lastly, we also noticed a subset of RAD21 interactions with CTCF binding in only one or none of the two anchor peaks, which were significantly enriched with enhancers. Taken together, this study outlines the general features of genome folding and opens new avenues to analyze spatial genome organization in specimens with limited cell numbers. Citation Format: Khyati Chandratre. Accurate prediction of cohesin-mediated 3D genome organization using 2D chromatin features [abstract]. In: Proceedings of the AACR Special Conference: Advances in Prostate Cancer Research; 2023 Mar 15-18; Denver, Colorado. Philadelphia (PA): AACR; Cancer Res 2023;83(11 Suppl):Abstract nr A054.

Conserved noncoding sequences correlate with distant gene contacts in Arabidopsis and Brassica.

Physical contact between genes distant on chromosomes is a potentially important way for genes to coordinate their expressions. To investigate the potential importance of distant contacts, we performed high-throughput chromatin conformation capture (Hi-C) experiments on leaf nuclei isolated from Brassica rapa and Brassica oleracea. We then combined our results with published Hi-C data from Arabidopsis thaliana. We found that distant genes come into physical contact and do so preferentially between the proximal promoter of one gene and the downstream region of another gene. Genes with higher numbers of conserved noncoding sequences (CNSs) nearby were more likely to have contact with distant genes. With more CNSs came higher numbers of transcription factor binding sites and more histone modifications associated with the activity. In addition, for the genes we studied, distant contacting genes with CNSs were more likely to be transcriptionally coordinated. These observations suggest that CNSs may enrich active histone modifications and recruit transcription factors, correlating with distant contacts to ensure coordinated expression. This study advances our knowledge of gene contacts and provides insights into the relationship between CNSs and distant gene contacts in plants.

The chromosome-scale genome assembly of Jasminum sambac var. unifoliatum provides insights into the formation of floral fragrance

Jasmine [Jasminum sambac (L.) Ait.], a tropical and subtropical plant emits a sweet, heady fragrance during flower opening. However, the molecular mechanisms underlying this phenomenon remain largely unknown. In the present study, integrated Illumina sequencing, Pacbio sequencing, and high-throughput chromatin conformation capture (Hi-C) scaffolding was used to generate a 495.60 Mb genome assembly of J. sambac var. unifoliatum cultivar ‘Fuzhou Single-petal’ (JSU-FSP), with contig N50 of 16.88 Mb; 96.23% of the assembly was assigned to 13 pseudochromosomes. The genome harbors 30 989 protein-coding genes, and 49.47% of the assembled sequences are repetitive sequences. The analysis of duplication modes showed that 51% of genes were duplicated through dispersed duplication, and expanded gene families are mainly involved in photosynthesis, which may be responsible for the light-loving characteristic specific to jasmine. Transcriptome analysis revealed that at least 35 structural genes involved in the biosynthesis of volatile terpenes (VTs), volatile phenylpropanoid/benzenoids (VPBs), fatty acid-derived volatiles (FADVs), and indole were highly expressed in the flower-opening stage, both preharvest and postharvest, and are proposed to be important in endowing flower aroma. Additionally, at least 28 heat shock protein (HSP) and 11 β-glucosidase (BGLU) genes may be involved in the formation of floral fragrance. These findings provide insights into the formation of the floral fragrance of jasmine and will promote germplasm utilization for breeding improved jasmine varieties.

A gap-free and haplotype-resolved lemon genome provides insights into flavor synthesis and huanglongbing (HLB) tolerance.

The lemon (Citrus limon; family Rutaceae) is one of the most important and popular fruits worldwide. Lemon also tolerates huanglongbing (HLB) disease, which is a devastating citrus disease. Here we produced a gap-free and haplotype-resolved chromosome-scale genome assembly of the lemon by combining Pacific Biosciences circular consensus sequencing, Oxford Nanopore 50-kb ultra-long, and high-throughput chromatin conformation capture technologies. The assembly contained nine-pair chromosomes with a contig N50 of 35.6Mb and zero gaps, while a total of 633.0Mb genomic sequences were generated. The origination analysis identified 338.5Mb genomic sequences originating from citron (53.5%), 147.4Mb from mandarin (23.3%), and 147.1Mb from pummelo (23.2%). The genome included 30 528 protein-coding genes, and most of the assembled sequences were found to be repetitive sequences. Several significantly expanded gene families were associated with plant-pathogen interactions, plant hormone signal transduction, and the biosynthesis of major active components, such as terpenoids and flavor compounds. Most HLB-tolerant genes were expanded in the lemon genome, such as 2-oxoglutarate (2OG)/Fe(II)-dependent oxygenase and constitutive disease resistance 1, cell wall-related genes, and lignin synthesis genes. Comparative transcriptomic analysis showed that phloem regeneration and lower levels of phloem plugging are the elements that contribute to HLB tolerance in lemon. Our results provide insight into lemon genome evolution, active component biosynthesis, and genes associated with HLB tolerance.

HPTAD: A computational method to identify topologically associating domains from HiChIP and PLAC-seq datasets

High-throughput chromatin conformation capture technologies, such as Hi-C and Micro-C, have enabled genome-wide view of chromatin spatial organization. Most recently, Hi-C-derived enrichment-based technologies, including HiChIP and PLAC-seq, offer attractive alternatives due to their high signal-to-noise ratio and low cost. While a series of computational tools have been developed for Hi-C data, methods tailored for HiChIP and PLAC-seq data are still under development. Here we present HPTAD, a computational method to identify topologically associating domains (TADs) from HiChIP and PLAC-seq data. We performed comprehensive benchmark analysis to demonstrate its superior performance over existing TAD callers designed for Hi-C data. HPTAD is freely available at https://github.com/yunliUNC/HPTAD.

Chromosome-level Assembly, Dosage Compensation and Sex-biased Gene Expression in the Small Brown Planthopper, Laodelphax striatellus.

In insects, sex chromosome differentiation often results in unequal gene dosages between sexes. Dosage compensation mechanisms evolve to balance gene expression, but the degree and mechanism of regulation often vary by insect species. In hemipteran species, the small brown planthopper (SBPH), Laodelphax striatellus, is an injurious crop pest, with a sex chromosome type XX in females and XO in males. This species offers the opportunity to study dosage compensation and sex-biased gene expression. In this study, we generated a chromosome-level genome of SBPH using Oxford Nanopore Technologies and high-throughput chromatin conformation capture (Hi-C) technology. We also sequenced RNA-seq data from 16 tissue samples to annotate the genome and analyze gene dosage compensation. We finally obtained a 510.2 megabases (Mb) genome with 99.12% of the scaffolds anchored on 15 chromosomes (14 autosomes and 1 X chromosome) and annotated 16,160 protein-coding genes based on full-length cDNA sequencing data. Furthermore, we found complete dosage compensation in all L. striatellus somatic tissues, but lack of dosage compensation in gonad tissue testis. We also found that female-biased genes were significantly enriched on the X chromosome in all tissues, whereas male-biased genes in gonad tissues were enriched on autosomes. This study not only provides a high-quality genome assembly but also lays a foundation for a better understanding of the sexual regulatory network in hemipteran insects.

A chromosome-level genome assembly of Styphnolobium japonicum combined with comparative genomic analyses offers insights on the evolution of flavonoid and lignin biosynthesis

Styphnolobium japonicum is an important industrial crop that is widely used as a source of flavonoids. In this study, the first chromosome-scale reference genome of S. japonicum was provided by combining Illumina and nanopore sequencing and high-throughput chromatin conformation capture technology. The genome of S. japonicum is 520 Mb and has an N50 contig that is 4.77 Mb. The contig could be clustered into 14 pseudochromosomes. A total of 32,551 genes that encode proteins were annotated in the genome. Comparative and phylogenetic analyses of the genome revealed the phylogenetic placement of this plant, which indicates that S. japonicum experienced two whole-genome duplication (WGD) events. The genes retained after the WGD were shared by members of the papilionoideae. A comparative genomic analysis revealed that the expansion and contraction of the flavonoid pathway gene families was the reason that the production of flavonoids was enhanced. This reference genome enabled the prediction of 17 genes that are involved in the flavonoid and lignin biosynthetic gene families. Collectively, S. japonicum could have undergone adaptive evolution to promote its growth by enhancing the production of flavonoids and altering the composition of lignin.

Normalization and de-noising of single-cell Hi-C data with BandNorm and scVI-3D

Single-cell high-throughput chromatin conformation capture methodologies (scHi-C) enable profiling of long-range genomic interactions. However, data from these technologies are prone to technical noise and biases that hinder downstream analysis. We develop a normalization approach, BandNorm, and a deep generative modeling framework, scVI-3D, to account for scHi-C specific biases. In benchmarking experiments, BandNorm yields leading performances in a time and memory efficient manner for cell-type separation, identification of interacting loci, and recovery of cell-type relationships, while scVI-3D exhibits advantages for rare cell types and under high sparsity scenarios. Application of BandNorm coupled with gene-associating domain analysis reveals scRNA-seq validated sub-cell type identification.

A chromosome-level genome assembly of the potato grouper (Epinephelus tukula)

The potato grouper, Epinephelus tukula, is one of the largest coral reef teleost, and it is an important germplasm resource for selection and cross breeding. Here we report a potato grouper genome assembly generated using PacBio long-read sequencing, Illumina sequencing and high-throughput chromatin conformation capture (Hi-C) technology. The genome size was 1.13 Gb, with a total of 508 contigs anchored into 24 chromosomes. The scaffold N50 was 42.65 Mb. For the genome models, our assembled genome contained 98.11% complete BUSCO with the vertebrata_odb9 database. One more copies of Gh and Hsp90b1 were identified in the E. tukula genome, which might contribute to its fast growth and high resistance to stress. In addition, 435 putative antimicrobial peptide (AMP) genes were identified in the potato grouper. This study provides a good reference for whole genome selective breeding of the potato grouper and for future development of novel marine drugs.