Haplotype-resolved 3D chromatin architecture of the hybrid pig.

In diploid mammals, allele-specific three-dimensional (3D) genome architecture may lead to imbalanced gene expression. Through ultradeep in situ Hi-C sequencing of three representative somatic tissues (liver, skeletal muscle, and brain) from hybrid pigs generated by reciprocal crosses of phenotypically and physiologically divergent Berkshire and Tibetan pigs, we uncover extensive chromatin reorganization between homologous chromosomes across multiple scales. Haplotype-based interrogation of multi-omic data revealed the tissue dependence of 3D chromatin conformation, suggesting that parent-of-origin-specific conformation may drive gene imprinting. We quantify the effects of genetic variations and histone modifications on allelic differences of long-range promoter-enhancer contacts, which likely contribute to the phenotypic differences between the parental pig breeds. We also observe the fine structure of somatically paired homologous chromosomes in the pig genome, which has a functional implication genome-wide. This work illustrates how allele-specific chromatin architecture facilitates concomitant shifts in allele-biased gene expression, as well as the possible consequential phenotypic changes in mammals.

Comparative three-dimensional genome architectures of adipose tissues provide insight into human-specific regulation of metabolic homeostasis.

Elucidating the regulatory mechanisms of human adipose tissues (ATs) evolution is essential for understanding human-specific metabolic regulation, but the functional importance and evolutionary dynamics of three-dimensional (3D) genome organizations of ATs are not well defined. Here, we compared the 3D genome architectures of anatomically distinct ATs from humans and six representative mammalian models. We recognized evolutionarily conserved and human-specific chromatin conformation in ATs at multiple scales, including compartmentalization, topologically associating domain (TAD), and promoter-enhancer interactions (PEI), which have not been described previously. We found PEI are much more evolutionarily dynamic with respect to compartmentalization and topologically associating domain. Compared to conserved PEIs, human-specific PEIs are enriched for human-specific sequence, and the binding motifs of their potential mediators (transcription factors) are less conserved. Our data also demonstrated that genes involved in the evolutionary dynamics of chromatin organization have weaker transcriptional conservation than those associated with conserved chromatin organization. Furthermore, the genes involved in energy metabolism and the maintenance of metabolic homeostasis are enriched in human-specific chromatin organization, while housekeeping genes, health-related genes, and genetic variations are enriched in evolutionarily conserved compared to human-specific chromatin organization. Finally, we showed extensively divergent human-specific 3D genome organizations among one subcutaneous and three visceral ATs. Together, these findings provide a global overview of 3D genome architecture dynamics between ATs from human and mammalian models and new insights into understanding the regulatory evolution of human ATs.

Comparative 3D genome architecture in vertebrates.

BACKGROUND: The three-dimensional (3D) architecture of the genome has a highly ordered and hierarchical nature, which influences the regulation of essential nuclear processes at the basis of gene expression, such as gene transcription. While the hierarchical organization of heterochromatin and euchromatin can underlie differences in gene expression that determine evolutionary differences among species, the way 3D genome architecture is affected by evolutionary forces within major lineages remains unclear. Here, we report a comprehensive comparison of 3D genomes, using high resolution Hi-C data in fibroblast cells of fish, chickens, and 10 mammalian species. RESULTS: This analysis shows a correlation between genome size and chromosome length that affects chromosome territory (CT) organization in the upper hierarchy of genome architecture, whereas lower hierarchical features, including local transcriptional availability of DNA, are selected through the evolution of vertebrates. Furthermore, conservation of topologically associating domains (TADs) appears strongly associated with the modularity of expression profiles across species. Additionally, LINE and SINE transposable elements likely contribute to heterochromatin and euchromatin organization, respectively, during the evolution of genome architecture. CONCLUSIONS: Our analysis uncovers organizational features that appear to determine the conservation and transcriptional regulation of functional genes across species. These findings can guide ongoing investigations of genome evolution by extending our understanding of the mechanisms shaping genome architecture.

Automatic Recognition of Giant Panda Attributes from Their Vocalizations Based on Squeeze-and-Excitation Network.

The giant panda (Ailuropoda melanoleuca) has long attracted the attention of conservationists as a flagship and umbrella species. Collecting attribute information on the age structure and sex ratio of the wild giant panda populations can support our understanding of their status and the design of more effective conservation schemes. In view of the shortcomings of traditional methods, which cannot automatically recognize the age and sex of giant pandas, we designed a SENet (Squeeze-and-Excitation Network)-based model to automatically recognize the attributes of giant pandas from their vocalizations. We focused on the recognition of age groups (juvenile and adult) and sex of giant pandas. The reason for using vocalizations is that among the modes of animal communication, sound has the advantages of long transmission distances, strong penetrating power, and rich information. We collected a dataset of calls from 28 captive giant panda individuals, with a total duration of 1298.02 s of recordings. We used MFCC (Mel-frequency Cepstral Coefficients), which is an acoustic feature, as inputs for the SENet. Considering that small datasets are not conducive to convergence in the training process, we increased the size of the training data via SpecAugment. In addition, we used focal loss to reduce the impact of data imbalance. Our results showed that the F1 scores of our method for recognizing age group and sex reached 96.46% ± 5.71% and 85.85% ± 7.99%, respectively, demonstrating that the automatic recognition of giant panda attributes based on their vocalizations is feasible and effective. This more convenient, quick, timesaving, and laborsaving attribute recognition method can be used in the investigation of wild giant pandas in the future.

Population genomics identifies patterns of genetic diversity and selection in chicken.

BACKGROUND: There are hundreds of phenotypically distinguishable domestic chicken breeds or lines with highly specialized traits worldwide, which provide a unique opportunity to illustrate how selection shapes patterns of genetic variation. There are many local chicken breeds in China. RESULTS: Here, we provide a population genome landscape of genetic variations in 86 domestic chickens representing 10 phenotypically diverse breeds. Genome-wide analysis indicated that sex chromosomes have less genetic diversity and are under stronger selection than autosomes during domestication and local adaptation. We found an evidence of admixture between Tibetan chickens and other domestic population. We further identified strong signatures of selection affecting genomic regions that harbor genes underlying economic traits (typically related to feathers, skin color, growth, reproduction and aggressiveness) and local adaptation (to high altitude). By comparing the genomes of the Tibetan and lowland fowls, we identified genes associated with high-altitude adaptation in Tibetan chickens were mainly involved in energy metabolism, body size maintenance and available food sources. CONCLUSIONS: The work provides crucial insights into the distinct evolutionary scenarios occurring under artificial selection for agricultural production and under natural selection for success at high altitudes in chicken. Several genes were identified as candidates for chicken economic traits and other phenotypic traits.

Dynamic changes in genes related to glucose uptake and utilization during pig skeletal and cardiac muscle development.

Skeletal and cardiac muscle have important roles in glucose uptake and utilization. However, changes in expression of protein coding genes and miRNAs that participate in glucose metabolism during development are not fully understood. In this study, we investigated the expression of genes related to glucose metabolism during muscle development. We found an age-dependent increase in gene expression in cardiac muscle, with enrichment in heart development- and energy-related metabolic processes. A subset of genes that were up-regulated until 30 or 180 days postnatally, and then down-regulated in psoas major muscle was significantly enriched in mitochondrial oxidative-related processes, while genes that up-regulated in longissimus doris muscle was significantly enriched in glycolysis-related processes. Meanwhile, expression of energy-related microRNAs decreased with increasing age. In addition, we investigated the correlation between microRNAs and mRNAs in three muscle types across different stages of development and found many potential microRNA-mRNA pairs involved in regulating glucose metabolism.

Exploring the dynamic three-dimensional chromatin architecture and transcriptional landscape in goose liver tissues underlying metabolic adaptations induced by a high-fat diet.

BACKGROUND: Goose, descendants of migratory ancestors, have undergone extensive selective breeding, resulting in their remarkable ability to accumulate fat in the liver and exhibit a high tolerance for significant energy intake. As a result, goose offers an excellent model for studying obesity, metabolic disorders, and liver diseases in mammals. Although the impact of the three-dimensional arrangement of chromatin within the cell nucleus on gene expression and transcriptional regulation is widely acknowledged, the precise functions of chromatin architecture reorganization during fat deposition in goose liver tissues still need to be fully comprehended. RESULTS: In this study, geese exhibited more pronounced changes in the liver index and triglyceride (TG) content following the consumption of the high-fat diet (HFD) than mice without significant signs of inflammation. Additionally, we performed comprehensive analyses on 10 goose liver tissues (5 HFD, 5 normal), including generating high-resolution maps of chromatin architecture, conducting whole-genome gene expression profiling, and identifying H3K27ac peaks in the livers of geese and mice subjected to the HFD. Our results unveiled a multiscale restructuring of chromatin architecture, encompassing Compartment A/B, topologically associated domains, and interactions between promoters and enhancers. The dynamism of the three-dimensional genome architecture, prompted by the HFD, assumed a pivotal role in the transcriptional regulation of crucial genes. Furthermore, we identified genes that regulate chromatin conformation changes, contributing to the metabolic adaptation process of lipid deposition and hepatic fat changes in geese in response to excessive energy intake. Moreover, we conducted a cross-species analysis comparing geese and mice exposed to the HFD, revealing unique characteristics specific to the goose liver compared to a mouse. These chromatin conformation changes help elucidate the observed characteristics of fat deposition and hepatic fat regulation in geese under conditions of excessive energy intake. CONCLUSIONS: We examined the dynamic modifications in three-dimensional chromatin architecture and gene expression induced by an HFD in goose liver tissues. We conducted a cross-species analysis comparing that of mice. Our results contribute significant insights into the chromatin architecture of goose liver tissues, offering a novel perspective for investigating mammal liver diseases.

Building Haplotype-Resolved 3D Genome Maps of Chicken Skeletal Muscle.

Haplotype-resolved 3D chromatin architecture related to allelic differences in avian skeletal muscle development has not been addressed so far, although chicken husbandry for meat consumption has been prevalent feature of cultures on every continent for more than thousands of years. Here, high-resolution Hi-C diploid maps (1.2-kb maximum resolution) are generated for skeletal muscle tissues in chicken across three developmental stages (embryonic day 15 to day 30 post-hatching). The sequence features governing spatial arrangement of chromosomes and characterize homolog pairing in the nucleus, are identified. Multi-scale characterization of chromatin reorganization between stages from myogenesis in the fetus to myofiber hypertrophy after hatching show concordant changes in transcriptional regulation by relevant signaling pathways. Further interrogation of parent-of-origin-specific chromatin conformation supported that genomic imprinting is absent in birds. This study also reveals promoter-enhancer interaction (PEI) differences between broiler and layer haplotypes in skeletal muscle development-related genes are related to genetic variation between breeds, however, only a minority of breed-specific variations likely contribute to phenotypic divergence in skeletal muscle potentially via allelic PEI rewiring. Beyond defining the haplotype-specific 3D chromatin architecture in chicken, this study provides a rich resource for investigating allelic regulatory divergence among chicken breeds.

Transcriptomic and lipidomic profiling of subcutaneous and visceral adipose tissues in 15 vertebrates.

The storage of lipids as energy in adipose tissue (AT) has been conserved over the course of evolution. However, substantial differences in ATs physiological activities were reported among species. Hence, establishing the mechanisms shaping evolutionarily divergence in ATs transcriptomes could provide a deeper understanding of AT regulation and its roles in obesity-related diseases. While previous studies performed anatomical, physiological and morphological comparisons between ATs across different species, little is currently understood at the molecular phenotypic levels. Here, we characterized transcriptional and lipidomic profiles of available subcutaneous and visceral ATs samples across 15 vertebrate species, spanning more than 300 million years of evolution, including placental mammals, birds and reptiles. We provide detailed descriptions of the datasets produced in this study and report gene expression and lipid profiles across samples. We demonstrate these data are robust and reveal the AT transcriptome and lipidome vary greater among species than within the same species. These datasets may serve as a resource for future studies on the functional differences among ATs in vertebrate species.

Reorganization of 3D genome architecture across wild boar and Bama pig adipose tissues.

BACKGROUND: A growing body of evidence has revealed that the mammalian genome is organized into hierarchical layers that are closely correlated with and may even be causally linked with variations in gene expression. Recent studies have characterized chromatin organization in various porcine tissues and cell types and compared them among species and during the early development of pigs. However, how chromatin organization differs among pig breeds is poorly understood. RESULTS: In this study, we investigated the 3D genome organization and performed transcriptome characterization of two adipose depots (upper layer of backfat [ULB] and greater omentum [GOM]) in wild boars and Bama pigs; the latter is a typical indigenous pig in China. We found that over 95% of the A/B compartments and topologically associating domains (TADs) are stable between wild boars and Bama pigs. In contrast, more than 70% of promoter-enhancer interactions (PEIs) are dynamic and widespread, involving over a thousand genes. Alterations in chromatin structure are associated with changes in the expression of genes that are involved in widespread biological functions such as basic cellular functions, endocrine function, energy metabolism and the immune response. Approximately 95% and 97% of the genes associated with reorganized A/B compartments and PEIs in the two pig breeds differed between GOM and ULB, respectively. CONCLUSIONS: We reported 3D genome organization in adipose depots from different pig breeds. In a comparison of Bama pigs and wild boar, large-scale compartments and TADs were mostly conserved, while fine-scale PEIs were extensively reorganized. The chromatin architecture in these two pig breeds was reorganized in an adipose depot-specific manner. These results contribute to determining the regulatory mechanism of phenotypic differences between Bama pigs and wild boar.

Giant panda age recognition based on a facial image deep learning system.

The conservation of the giant panda (Ailuropoda melanoleuca), as an iconic vulnerable species, has received great attention in the past few decades. As an important part of the giant panda population survey, the age distribution of giant pandas can not only provide useful instruction but also verify the effectiveness of conservation measures. The current methods for determining the age groups of giant pandas are mainly based on the size and length of giant panda feces and the bite value of intact bamboo in the feces, or in the case of a skeleton, through the wear of molars and the growth line of teeth. These methods have certain flaws that limit their applications. In this study, we developed a deep learning method to study age group classification based on facial images of captive giant pandas and achieved an accuracy of 85.99% on EfficientNet. The experimental results show that the faces of giant pandas contain some age information, which mainly concentrated between the eyes of giant pandas. In addition, the results also indicate that it is feasible to identify the age groups of giant pandas through the analysis of facial images.

Dynamic transcriptome and chromatin architecture in granulosa cells during chicken folliculogenesis.

Folliculogenesis is a complex biological process involving a central oocyte and its surrounding somatic cells. Three-dimensional chromatin architecture is an important transcription regulator; however, little is known about its dynamics and role in transcriptional regulation of granulosa cells during chicken folliculogenesis. We investigate the transcriptomic dynamics of chicken granulosa cells over ten follicular stages and assess the chromatin architecture dynamics and how it influences gene expression in granulosa cells at three key stages: the prehierarchical small white follicles, the first largest preovulatory follicles, and the postovulatory follicles. Our results demonstrate the consistency between the global reprogramming of chromatin architecture and the transcriptomic divergence during folliculogenesis, providing ample evidence for compartmentalization rearrangement, variable organization of topologically associating domains, and rewiring of the long-range interaction between promoter and enhancers. These results provide key insights into avian reproductive biology and provide a foundational dataset for the future in-depth functional characterization of granulosa cells.

Generation and characterization of stable pig pregastrulation epiblast stem cell lines.

Pig epiblast-derived pluripotent stem cells are considered to have great potential and broad prospects for human therapeutic model development and livestock breeding. Despite ongoing attempts since the 1990s, no stably defined pig epiblast-derived stem cell line has been established. Here, guided by insights from a large-scale single-cell transcriptome analysis of pig embryos from embryonic day (E) 0 to E14, specifically, the tracing of pluripotency changes during epiblast development, we developed an in vitro culture medium for establishing and maintaining stable pluripotent stem cell lines from pig E10 pregastrulation epiblasts (pgEpiSCs). Enabled by chemical inhibition of WNT-related signaling in combination with growth factors in the FGF/ERK, JAK/STAT3, and Activin/Nodal pathways, pgEpiSCs maintain their pluripotency transcriptome features, similar to those of E10 epiblast cells, and normal karyotypes after more than 240 passages and have the potential to differentiate into three germ layers. Strikingly, ultradeep in situ Hi-C analysis revealed functional impacts of chromatin 3D-spatial associations on the transcriptional regulation of pluripotency marker genes in pgEpiSCs. In practice, we confirmed that pgEpiSCs readily tolerate at least three rounds of successive gene editing and generated cloned gene-edited live piglets. Our findings deliver on the long-anticipated promise of pig pluripotent stem cells and open new avenues for biological research, animal husbandry, and regenerative biomedicine.

Association of female reproductive tract microbiota with egg production in layer chickens.

BACKGROUND: The microbiota of the female reproductive tract is increasingly recognized as playing fundamental roles in animal reproduction. To explore the relative contribution of reproductive tract microbiomes to egg production in chickens, we investigated the microbiota in multiple reproductive and digestive tract sites from 128 female layer (egg-producing) chickens in comparable environments. RESULTS: We identified substantial differences between the diversity, composition, and predicted function of site-associated microbiota. Differences in reproductive tract microbiota were more strongly associated with egg production than those in the digestive tract. We identified 4 reproductive tract microbial species, Bacteroides fragilis, Bacteroides salanitronis, Bacteroides barnesiae, and Clostridium leptum, that were related to immune function and potentially contribute to enhanced egg production. CONCLUSIONS: These findings provide insights into the diverse microbiota characteristics of reproductive and digestive tracts and may help in designing strategies for controlling and manipulating chicken reproductive tract microbiota to improve egg production.

A pig BodyMap transcriptome reveals diverse tissue physiologies and evolutionary dynamics of transcription.

A comprehensive transcriptomic survey of pigs can provide a mechanistic understanding of tissue specialization processes underlying economically valuable traits and accelerate their use as a biomedical model. Here we characterize four transcript types (lncRNAs, TUCPs, miRNAs, and circRNAs) and protein-coding genes in 31 adult pig tissues and two cell lines. We uncover the transcriptomic variability among 47 skeletal muscles, and six adipose depots linked to their different origins, metabolism, cell composition, physical activity, and mitochondrial pathways. We perform comparative analysis of the transcriptomes of seven tissues from pigs and nine other vertebrates to reveal that evolutionary divergence in transcription potentially contributes to lineage-specific biology. Long-range promoter-enhancer interaction analysis in subcutaneous adipose tissues across species suggests evolutionarily stable transcription patterns likely attributable to redundant enhancers buffering gene expression patterns against perturbations, thereby conferring robustness during speciation. This study can facilitate adoption of the pig as a biomedical model for human biology and disease and uncovers the molecular bases of valuable traits.

Genomic analyses provide insights into breed-of-origin effects from purebreds on three-way crossbred pigs.

Crossbreeding is widely used aimed at improving crossbred performance for poultry and livestock. Alleles that are specific to different purebreds will yield a large number of heterozygous single-nucleotide polymorphisms (SNPs) in crossbred individuals, which are supposed to have the power to alter gene function or regulate gene expression. For pork production, a classic three-way crossbreeding system of Duroc × (Landrace × Yorkshire) (DLY) is generally used to produce terminal crossbred pigs with stable and prominent performance. Nonetheless, little is known about the breed-of-origin effects from purebreds on DLY pigs. In this study, we first estimated the distribution of heterozygous SNPs in three kinds of three-way crossbred pigs via whole genome sequencing data originated from three purebreds. The result suggested that DLY is a more effective strategy for three-way crossbreeding as it could yield more stably inherited heterozygous SNPs. We then sequenced a DLY pig family and identified 95, 79, 132 and 42 allele-specific expression (ASE) genes in adipose, heart, liver and skeletal muscle, respectively. Principal component analysis and unrestricted clustering analyses revealed the tissue-specific pattern of ASE genes, indicating the potential roles of ASE genes for development of DLY pigs. In summary, our findings provided a lot of candidate SNP markers and ASE genes for DLY three-way crossbreeding system, which may be valuable for pig breeding and production in the future.

Genome-Wide Chromatin Structure Changes During Adipogenesis and Myogenesis.

The recently developed high-throughput chromatin conformation capture (Hi-C) technology enables us to explore the spatial architecture of genomes, which is increasingly considered an important regulator of gene expression. To investigate the changes in three-dimensional (3D) chromatin structure and its mediated gene expression during adipogenesis and myogenesis, we comprehensively mapped 3D chromatin organization for four cell types (3T3-L1 pre-adipocytes, 3T3-L1-D adipocytes, C2C12 myoblasts, and C2C12-D myotubes). We demonstrate that the dynamic spatial genome architecture affected gene expression during cell differentiation. A considerable proportion (~22%) of the mouse genome underwent compartment A/B rearrangement during adipogenic and myogenic differentiation, and most (~80%) upregulated marker genes exhibited an active chromatin state with B to A switch or stable A compartment. More than half (65.4%-73.2%) of the topologically associating domains (TADs) are dynamic. The newly formed TAD and intensified local interactions in the Fabp gene cluster indicated more precise structural regulation of the expression of pro-differentiation genes during adipogenesis. About half (32.39%-59.04%) of the differential chromatin interactions (DCIs) during differentiation are promoter interactions, although these DCIs only account for a small proportion of genome-wide interactions (~9.67% in adipogenesis and ~4.24% in myogenesis). These differential promoter interactions were enriched with promoter-enhancer interactions (PEIs), which were mediated by typical adipogenic and myogenic transcription factors. Differential promoter interactions also included more differentially expressed genes than nonpromoter interactions. Our results provide a global view of dynamic chromatin interactions during adipogenesis and myogenesis and are a resource for studying long-range chromatin interactions mediating the expression of pro-differentiation genes.

Limited Cross-Linking of 4-1BB by 4-1BB Ligand and the Agonist Monoclonal Antibody Utomilumab.

Monoclonal antibodies (mAbs) targeting the co-stimulatory molecule 4-1BB are of interest for tumor immunotherapy. We determined the complex structures of human 4-1BB with 4-1BB ligand (4-1BBL) or utomilumab to elucidate the structural basis of 4-1BB activation. The 4-1BB/4-1BBL complex displays a typical TNF/TNFR family binding mode. The structure of utomilumab/4-1BB complex shows that utomilumab binds to dimeric 4-1BB with a distinct but partially overlapping binding area with 4-1BBL. Competitive binding analysis demonstrates that utomilumab blocks the 4-1BB/4-1BBL interaction, indicating the interruption of ligand-mediated signaling. The binding profiles of 4-1BBL and utomilumab to monomeric or dimeric 4-1BB indicate limited cross-linking of 4-1BB molecules. These findings provide mechanistic insight into the binding of 4-1BB with its ligand and its agonist mAb, which may facilitate the future development of anti-4-1BB biologics for tumor immunotherapy.

mRNA N6-methyladenosine methylation of postnatal liver development in pig.

N6-methyladenosine (m6A) is a ubiquitous reversible epigenetic RNA modification that plays an important role in the regulation of post-transcriptional protein coding gene expression. Liver is a vital organ and plays a major role in metabolism with numerous functions. Information concerning the dynamic patterns of mRNA m6A methylation during postnatal development of liver has been long overdue and elucidation of this information will benefit for further deciphering a multitude of functional outcomes of mRNA m6A methylation. Here, we profile transcriptome-wide m6A in porcine liver at three developmental stages: newborn (0 day), suckling (21 days) and adult (2 years). About 33% of transcribed genes were modified by m6A, with 1.33 to 1.42 m6A peaks per modified gene. m6A was distributed predominantly around stop codons. The consensus motif sequence RRm6ACH was observed in 78.90% of m6A peaks. A negative correlation (average Pearson's r = -0.45, P < 10-16) was found between levels of m6A methylation and gene expression. Functional enrichment analysis of genes consistently modified by m6A methylation at all three stages showed genes relevant to important functions, including regulation of growth and development, regulation of metabolic processes and protein catabolic processes. Genes with higher m6A methylation and lower expression levels at any particular stage were associated with the biological processes required for or unique to that stage. We suggest that differential m6A methylation may be important for the regulation of nutrient metabolism in porcine liver.

Remarkably similar CTLA-4 binding properties of therapeutic ipilimumab and tremelimumab antibodies.

Monoclonal antibody based immune checkpoint blockade therapies have achieved clinical successes in management of malignant tumors. As the first monoclonal antibody targeting immune checkpoint molecules entered into clinics, the molecular basis of ipilimumab-based anti-CTLA-4 blockade has not yet been fully understood. In the present study, we report the complex structure of ipilimumab and CTLA-4. The complex structure showed similar contributions from VH and VL of ipilimumab in binding to CTLA-4 front ß-sheet strands. The blockade mechanism of ipilimumab is that the strands of CTLA-4 contributing to the binding to B7-1 or B7-2 were occupied by ipilimumab and thereafter prevents the binding of B7-1 or B7-2 to CTLA-4. Though ipilimumab binds to the same epitope with tremelimumab on CTLA-4 with similar binding affinity, the higher dissociation rate of ipilimumab may indicate the dynamic binding to CTLA-4, which may affect its pharmacokinetics. The molecular basis of ipilimumab-based anti-CTLA-4 blockade and comparative study of the binding characteristics of ipilimumab and tremelimumab would shed light for the discovery of small molecular inhibitors and structure-based monoclonal antibody optimization or new biologics.