Combined analysis of metagenomic data revealed consistent changes of gut microbiome structure and function in inflammatory bowel disease

AIMS

To reveal the consistency and discrepancy in the gut microbial structure and function in inflammatory bowel disease (IBD) patients from different regions.

METHODS AND RESULTS

Gut microbes, antibiotic resistance genes (ARGs) and virulence factors genes (VFGs) were analysed using metagenome data from three cohorts. The abundance of Escherichia coli extensively increased in IBD patients, whereas Subdoligranulum unclassified decreased dramatically in IBD patients from three countries. Escherichia coli showed a positive correlation with multiple ARGs and VFGs in cohorts from China and the United States, including multidrug-related resistance genes and Capsule and LOS-related virulence factors genes. Escherichia coli biofilm synthesis pathways significantly enriched in IBD patients from three different regions. Notably, Subdoligranulum unclassified and Eubacterium hallii were negatively related to ARGs and VFGs.

CONCLUSIONS

Consistent changes of microbiome structure and function were observed in IBD patients from three different regions. As pathogenic bacteria, E. coli may accelerate IBD progression through encapsulation in biofilms by upregulating antibiotic resistance in Crohn's disease patients. Subdoligranulum unclassified and E. hallii may be beneficial for IBD patients and could serve as potential probiotics for IBD treatment.

SIGNIFICANCE AND IMPACT OF THE STUDY

This work dispels worries about the regional differences in gut microbial changes in IBD patients and provides useful guidance for more rational microbiome-based therapies.

5' UTR CGG repeat expansion in GIPC1 is associated with oculopharyngodistal myopathy

Oculopharyngodistal myopathy is a late-onset degenerative muscle disorder characterized by ptosis and weakness of the facial, pharyngeal, and distal limb muscles. A recent report suggested a non-coding trinucleotide repeat expansion in LRP12 to be associated with the disease. Here we report a genetic study in a Chinese cohort of 41 patients with the clinical diagnosis of oculopharyngodistal myopathy (21 cases from seven families and 20 sporadic cases). In a large family with 12 affected individuals, combined haplotype and linkage analysis revealed a maximum two-point logarithm of the odds (LOD) score of 3.3 in chromosomal region chr19p13.11-p13.2 and narrowed the candidate region to an interval of 4.5 Mb. Using a comprehensive strategy combining whole-exome sequencing, long-read sequencing, repeat-primed polymerase chain reaction and GC-rich polymerase chain reaction, we identified an abnormal CGG repeat expansion in the 5' UTR of the GIPC1 gene that co-segregated with disease. Overall, the repeat expansion in GIPC1 was identified in 51.9% independent pedigrees (4/7 families and 10/20 sporadic cases), while the repeat expansion in LRP12 was only identified in one sporadic case (3.7%) in our cohort. The number of CGG repeats was <30 in controls but >60 in affected individuals. There was a slight correlation between repeat size and the age at onset. Both repeat expansion and retraction were observed during transmission but somatic instability was not evident. These results further support that non-coding CGG repeat expansion plays an essential role in the pathogenesis of oculopharyngodistal myopathy.

Endocrine and genetic factors affecting egg laying performance in chickens: a review

1. Egg-laying performance reflects the overall reproductive performance of breeding hens. The genetic traits for egg-laying performance have low or medium heritability, and, depending on the period involved, usually ranges from 0.16 to 0.64. Egg-laying in chickens is regulated by a combination of environmental, endocrine and genetic factors. 2. The main endocrine factors that regulate egg-laying are gonadotropin-releasing hormone (GnRH), prolactin (PRL), follicle-stimulating hormone (FSH) and luteinising hormone (LH). 3. In the last three decades, many studies have explored this aspect at a molecular genetic level. Recent studies identified 31 reproductive hormone-based candidate genes that were significantly associated with egg-laying performance. With the development of genome-sequencing technology, 64 new candidate genes and 108 single nucleotide polymorphisms (SNPs) related to egg-laying performance have been found using genome-wide association studies (GWAS), providing novel insights into the molecular genetic mechanisms governing egg production. At the same time, microRNAs that regulate genes responsible for egg-laying in chickens were reviewed. 4. Research on endocrinological and genetic factors affecting egg-laying performance will greatly improve the reproductive performance of chickens and promote the protection, development, and utilisation of poultry. This review summarises studies on the endocrine and genetic factors of egg-laying performance in chickens from 1972 to 2019.

Making global river ecosystem health assessments objective, quantitative and comparable

Assessing and comparing global river ecosystem health in an objective and quantitative way remains a major challenge. In this study the widely-used semi-quantitative methods Rapid Biological assessment Protocols (RBPs) was used to determine the health of rivers. The findings were then compared to the results derived from our new UAV (Unmanned aerial vehicles) orthophotographic imagery method. This method quantitatively and objectively assesses river ecosystem health. As a comparison, our method was used to quantitatively measure distance and areas of a range of hydrological and biological attributes thus improving the accuracy of distance- and area-related indices, consequently avoiding subjective errors in these estimations that is fraught in methods like the RBPs. To strengthen the objectivity of the assessment the weights of these indices were objectively determined using the entropy weighting method. This new method was then tested using 9551 UAV orthophotographs taken over six field campaigns. It performed satisfactorily, showing that in our study area the health status of mountain rivers was the best with the highest score of 0.94 out of 1.0. Temporally, the health of the river was better in summer (0.65) compared with that in autumn (0.40). Changes in river ecosystem health were driven by variations in biology and water quality. In contrast the outputs of RBPs, especially in relation to distance and area indices, had ~ 20% uncertainty due to visual errors and subjectivity in estimations by observers. The UAV orthophotographic imaging method proposed in this study can improve the ability to compare the health of rivers across different periods and regions throughout the globe.

Identification of gene mutations in patients with primary periodic paralysis using targeted next-generation sequencing

Background

Primary periodic paralysis is characterized by recurrent quadriplegia typically associated with abnormal serum potassium levels. The molecular diagnosis of primary PP previously based on Sanger sequencing of hot spots or exon-by-exon screening of the reported genes.

Methods

We developed a gene panel that includes 10 ion channel-related genes and 245 muscular dystrophy- and myopathy-related genes and used this panel to diagnose 60 patients with primary periodic paralysis and identify the disease-causing or risk-associated gene mutations.

Results

Mutations of 5 genes were discovered in 39 patients (65.0%). SCN4A, KCNJ2 and CACNA1S variants accounted for 92.5% of the patients with a genetic diagnosis.

Conclusions

Targeted next-generation sequencing offers a cost-effective approach to expand the genotypes of primary periodic paralysis. A clearer genetic profile enables the prevention of paralysis attacks, avoidance of triggers and the monitoring of complications.

Electronic supplementary material

The online version of this article (10.1186/s12883-019-1322-6) contains supplementary material, which is available to authorized users.

Transcriptome profiling and pathway analysis of the effects of mono-(2-ethylhexyl) phthalate in mouse Sertoli cells

Phthalates are confirmed to have toxic effects on the reproductive system and are likely to have further damaging actions in humans. The present study explored the molecular mechanisms of the toxic effect of mono-(2-ethylhexyl) phthalate (MEHP) on mouse Sertoli cells. Cell apoptosis and proliferation assays were used to assess the effects of MEHP on the TM4 Sertoli cell line derived from mouse testes. TM4 cells were treated with two doses of MEHP or left untreated as a control group, followed by RNA extraction and analysis using high-throughput transcriptome sequencing technology. The gene expression profile obtained was then subjected to a bioinformatics analysis to explore the molecular mechanisms of reproductive toxicity. The results revealed that 528 and 269 genes were upregulated in the high- and low-dose MEHP groups of cells compared with the control group, while 148 and 173 genes were downregulated. Gene ontology (GO) analysis indicated that the differently expressed genes were associated with the GO term 'extracellular region' of the cellular component domain in the high and low MEHP groups. Compared with the control group, eight common pathway changes were identified in the high- and low-dose MEHP groups, including 'terpenoid backbone biosynthesis'. Reverse transcription-quantitative polymerase chain reaction analysis was used to validation, and hermetic effects were observed for certain genes. These results provide an important basis and experimental data for further research into the mechanisms of phthalate-induced toxicity.

A Comprehensive Analysis of the Transcriptomes of Marssonina brunnea and Infected Poplar Leaves to Capture Vital Events in Host-Pathogen Interactions

Background

Understanding host-pathogen interaction mechanisms helps to elucidate the entire infection process and focus on important events, and it is a promising approach for improvement of disease control and selection of treatment strategy. Time-course host-pathogen transcriptome analyses and network inference have been applied to unravel the direct or indirect relationships of gene expression alterations. However, time series analyses can suffer from absent time points due to technical problems such as RNA degradation, which limits the application of algorithms that require strict sequential sampling. Here, we introduce an efficient method using independence test to infer an independent network that is exclusively concerned with the frequency of gene expression changes.

Results

Highly resistant NL895 poplar leaves and weakly resistant NL214 leaves were infected with highly active and weakly active Marssonina brunnea, respectively, and were harvested at different time points. The independent network inference illustrated the top 1,000 vital fungus-poplar relationships, which contained 768 fungal genes and 54 poplar genes. These genes could be classified into three categories: a fungal gene surrounded by many poplar genes; a poplar gene connected to many fungal genes; and other genes (possessing low degrees of connectivity). Notably, the fungal gene M6_08342 (a metalloprotease) was connected to 10 poplar genes, particularly including two disease-resistance genes. These core genes, which are surrounded by other genes, may be of particular importance in complicated infection processes and worthy of further investigation.

Conclusions

We provide a clear framework of the interaction network and identify a number of candidate key effectors in this process, which might assist in functional tests, resistant clone selection, and disease control in the future.

Comparing Mycobacterium tuberculosis genomes using genome topology networks

Background

Over the last decade, emerging research methods, such as comparative genomic analysis and phylogenetic study, have yielded new insights into genotypes and phenotypes of closely related bacterial strains. Several findings have revealed that genomic structural variations (SVs), including gene gain/loss, gene duplication and genome rearrangement, can lead to different phenotypes among strains, and an investigation of genes affected by SVs may extend our knowledge of the relationships between SVs and phenotypes in microbes, especially in pathogenic bacteria.

Results

In this work, we introduce a ‘Genome Topology Network’ (GTN) method based on gene homology and gene locations to analyze genomic SVs and perform phylogenetic analysis. Furthermore, the concept of ‘unfixed ortholog’ has been proposed, whose members are affected by SVs in genome topology among close species. To improve the precision of 'unfixed ortholog' recognition, a strategy to detect annotation differences and complete gene annotation was applied. To assess the GTN method, a set of thirteen complete M. tuberculosis genomes was analyzed as a case study. GTNs with two different gene homology-assigning methods were built, the Clusters of Orthologous Groups (COG) method and the orthoMCL clustering method, and two phylogenetic trees were constructed accordingly, which may provide additional insights into whole genome-based phylogenetic analysis. We obtained 24 unfixable COG groups, of which most members were related to immunogenicity and drug resistance, such as PPE-repeat proteins (COG5651) and transcriptional regulator TetR gene family members (COG1309).

Conclusions

The GTN method has been implemented in PERL and released on our website. The tool can be downloaded from http://homepage.fudan.edu.cn/zhouyan/gtn/, and allows re-annotating the ‘lost’ genes among closely related genomes, analyzing genes affected by SVs, and performing phylogenetic analysis. With this tool, many immunogenic-related and drug resistance-related genes were found to be affected by SVs in M. tuberculosis genomes. We believe that the GTN method will be suitable for the exploration of genomic SVs in connection with biological features of bacterial strains, and that GTN-based phylogenetic analysis will provide additional insights into whole genome-based phylogenetic analysis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1259-0) contains supplementary material, which is available to authorized users.

Evidence of positive selection at signal peptide region of interferon gamma

Interferon gamma (IFNG) is a major cytokine and plays crucial roles in pathogen clearance. About the course of evolution of IFNG, it has been reported that IFNG is being subjected to adaptive selection, which is proved at the level of gene. Neighbor-joining method was used to reconstruct the phylogenetic tree of all IFNG protein-coding sequences. The pair-wise computation of Ka/Ks between every exon homologs, branch-specific model, and site-specific model of the likelihood method were performed to detect positive selection of IFNG. We reported, for the first time, that the signal peptide region of IFNG is under significant positive selection, evolving faster than other parts. We provide evidence at the level of individual exon and individual amino acid site that IFNG is under adaptive evolution, which establishes the basis for further researches about IFNG.

The evolution of alternative splicing exons in vascular endothelial growth factor A

The C-terminus alternative splicing in VEGFA (vascular endothelial growth factor A) is known for its impact on physiological and pathological angiogenesis. Based on our prediction and RT-PCR verification, we identified anti-angiogenic VEGFA165b isoforms in mouse and rabbit for the first time. We also found that the relative expression level of VEGFA165b isoform had been increasing from rodents to human, and exon8b may have experienced a minor-to-major form exon conversion, possibly correlated with its gain-of-function. It is suggested that introduction of alternative splicing exons (esp. exon6 and exon8b) made important contributions to the transcriptional diversity of VEGFA and played a crucial role in the evolution of its regulatory mechanism.

An intronic polymorphism rs2237062 in the CXCL14 gene influences HBV-related HCC progression in Chinese population

CXCL14 (C-X-C motif chemokine ligand 14) is a conserved member of chemokine family and functions as a chemoattractant for multiplicate immunocytes. CXCL14 expression is constitutive in normal tissues, but absent in wide range of epithelial tumors. Many reports have claimed its important role in tumorigenesis and vascularization. An association between rs2237062 polymorphism and hepatocellular carcinoma (HCC) susceptibility was found in patients with chronic HCV infection in Japanese population. Here we analyzed, by using a polymerase chain reaction-ligation detection reaction (PCR-LDR), the polymorphism in 202 non-HCC patients with HBV infection, 361 HBV-related HCC patients and 407 healthy controls. The aim was to detect the possible association of this single-nucleotide polymorphism (SNP) with HBV-related HCC susceptibility and progression. However, no association was found between rs2237062 polymorphism and susceptibility to HBV infection or HBV-related HCC. Intriguingly, our stratification analysis revealed that HBV-related HCC patients in advanced phase (TNM-II-IV stage) had significantly higher C allele frequency at this polymorphism than patients at early stage (TNM-I stage) (33.5% vs. 25.7%), and its odds ratio reached 1.47 (95% CI 1.06-2.04, P = 0.021). These results suggest that the rs2237062 polymorphism in the CXCL14 gene might influence HBV-related HCC progression in Chinese population.

BMI1 promotes the progression of laryngeal squamous cell carcinoma

BMI1 is highly expressed in several malignant tumors, and its expression level is associated with tumor progression, proliferation, and prognosis. However, no published studies have examined the role of BMI1 in laryngeal squamous cell carcinoma (SCC). Expression of BMI1 in primary tumors was analyzed by immunofluorescence staining, real-time PCR, and Western blotting. BMI1 was knocked down, and proliferation, apoptosis, and cell cycle assays were performed. Sensitivity to radiochemotherapy was evaluated, and tumorigenicity assays were performed in vivo. BMI1 was highly expressed in laryngeal SCCs. BMI1 promoted cell proliferation and tumor progression, and inhibited apoptosis due to influences on the cell cycle. More importantly, BMI1 suppressed the sensitization of laryngeal Hep2 cells to radiochemotherapy. BMI1 is essential to maintain the proliferation and progression of laryngeal SCCs. Therefore, depletion of BMI1 may be a potential therapeutic option for cancer management.

Dysregulation of Claudin family genes in colorectal cancer in a Chinese population

Claudins play an important role in tumor metastasis and in invasiveness of colorectal cancer (CRC). We have evaluated the relationship between CRC and expression of the claudin genes in Chinese patients with CRC. We measured CLDN1 and CLDN7 mRNA using quantitative PCR, and protein levels with immunohistochemistry in cancer tissues and adjacent normal tissue. Cancer tissues had significantly higher levels of CLDN1, and significantly lower levels of CLDN3, CLDN4, and CLDN7 than did normal tissue. CLDN3, CLDN4, and CLDN7 expression levels were higher in CRC of the protruded type than in CRC of the infiltrative type. Expression of CLDN7 correlated with lymph node metastasis. Stage N0 cancer tissues had higher levels of CLDN7 than did stages N1 and N2, suggesting that CLDN7 expression was closely related to the extent of lymph node metastasis. CLDN1 protein was upregulated, but CLDN7 protein was downregulated in cancer tissues when compared with expression in adjacent normal tissues. In conclusion, CLDN3, CLDN4, and CLDN7 were significantly downregulated, whereas CLDN1 was significantly upregulated in CRC. The altered expression of claudin genes may play a role in the initiation and development of CRC.

Evolution of peroxisome proliferator-activated receptor gamma alternative splicing

Peroxisome proliferator-activated receptor gamma (PPAR gamma) plays an important role in the control of energy balance and lipid and glucose homeostasis. Different transcript variants were investigated not only in human but also in other vertebrates. To look into the evolutionary changes of these variants, we analyzed the genomic sequences of PPAR gamma genes from several vertebrate species, as well as their mRNA and EST data. Several potential alternative splicing exons at the 5'-end of the PPAR gamma gene were identified. The 5'-end of the PPAR gamma gene is discovered to be evolutionarily active and recruits new exons via different strategies. Moreover, it is shown that the only coding alternative exon (exon B) processes much higher Ka/Ks compared with its constitutive counterparts. In addition, its Ka/Ks is greater than 1 in the rat, mouse, and rabbit, indicating adaptive evolution and possible energy storage related gain-of-function for the exon.

A polymorphism of microRNA196a genome region was associated with decreased risk of glioma in Chinese population

MicroRNAs (miRNAs) are small non-coding RNAs that play important roles in regulation of eukaryotic gene expression. Aberrant expression and structural alternation of miRNAs are considered to participate in tumorigenesis and cancer development. Recently, different genotypes of miR-196a polymorphisms (SNP, rs11614913) were found to be associated with the survival of patients with lung cancer and increased risk of breast cancer. To further investigate whether this polymorphism may influence glioma risk or not, we examined the SNP allele frequency in Chinese population. Our data shows the genotype CC of miR-196a (rs11614913) polymorphism is associated with decreased risk of glioma in the Chinese population (OR = 0.74, 95% CI:0.56-0.98). Furthermore, a significant association was observed between this genotype and glioma risk in the subgroups of adult glioma (OR = 0.73, 95% CI:0.55-0.98), male glioma (OR = 0.69, 95% CI:0.48-0.99) and patients with glioblastoma (OR = 0.58, 95% CI:0.37-0.91). This was the first study investigating the association between the miR-196a rs11614913 and glioma risk. Compared with the results from previous studies in lung cancer and breast cancer, our data suggest a different genotype association in glioma. This may be related to the diversity on the tissue origin, tumor type, tumorigenesis, and developing process.

Genome-wide prediction of G protein-coupled receptors in Verticillium spp

G protein-coupled receptors (GPCRs) are critical factors in regulating morphogenesis, mating, infection and virulence in fungi. In this study, various computational strategies were applied to identify GPCR-like proteins from the genomes of both Verticillium dahliae and Verticillium albo-atrum. The putative GPCRs were distributed over 13 classes, and significantly, three of those represented novel classes of GPCR-like proteins in fungi. The three novel GPCRs had high levels of identity to their counterparts in higher eukaryotes, including Homo sapiens. The numbers of GPCR-like proteins in the two Verticillium spp. were similar to those seen in other filamentous fungi, such as Magnaporthe grisea, Neurospora crassa and Fusarium graminearum. Additionally, the carbon/amino acid receptors were divided into three different subclasses, indicating that differences among the GPCRs existed not only among different classes but also within classes. In conclusion, the identification and classification of GPCRs and their homology to some well-studied fungi will be an important starting point for future research in Verticillium spp.

An intronic polymorphism in the corticotropin-releasing hormone receptor 2 gene increases susceptibility to HBV-related hepatocellular carcinoma in Chinese population

Corticotropin-releasing hormone receptor 2 (CRHR2) plays a role in both the central nervous system (CNS) and the peripheral nervous system. CRHR2 together with its ligands, urocortins (Ucns) and corticotropin-releasing hormone (CRH), functions as a mediator of inflammatory response and inhibitor of angiogenesis. Recently, it has been reported to be expressed in many human cancers. An association between rs2267716 polymorphism in the CRHR2 gene and susceptibility to hepatocellular carcinoma (HCC) was found in patients with chronic hepatitis C virus (HCV) infection. In the present study we analyzed, using a polymerase chain reaction-ligation detection reaction (PCR-LDR), the rs2267716 polymorphism in 364 hepatitis B virus (HBV)-related HCC patients, 196 non-HCC patients with HBV infection, and 404 healthy controls. The aim was to detect the possible association of this single-nucleotide polymorphism (SNP) with susceptibility to HBV-related HCC. Significant differences of rs2267716 allele were detected between HBV-related HCC patients and healthy controls (OR = 1.55, 95% CI 1.13-2.15, P = 0.007) or non-HCC patients with HBV infection (OR = 1.61, 95% CI 1.13-2.31, P = 0.009). These results suggest that the rs2267716 polymorphism in the CRHR2 gene might influence the risk of developing HCC in patients with HBV infection in Chinese population.

Cloning and Characterization of a novel splice variant of humanRab18gene (RAB18)

Rab GTPase proteins are a kind of small GTP-binding proteins, which functions mainly focus on regulating interacellular trafficking pathways during vesicular transport. To date, 60 distinct human RAB proteins have been identified. RAB18 gene is discovered from endothelial cells. Its function is considered as endosomes and plasma membrane recycling. Research indicates RAB18 may relate to inflammation and some kinds of tumor. Here we report a splice variant of RAB18, which is 2571 bp in length and has an open reading frame coding a predicted 235 amino-acids protein. RT-PCR shows that the cDNA has different expression pattern with RAB18 and is highly expressed in testis.

Cloning and characterization of a novel splice variant of human U2AF1L3 gene

Pre-mRNA splicing allows individual genes to produce multiple protein isoforms with diverse functions. Recognition of functional splice sites in pre-mRNAs is very important in this splicing process and requires some protein auxiliary factors such as U2 small nuclear ribonucleoprotein auxiliary factor small subunit (U2AF35, encoded by U2AF1). By its RNA binding domains, U2AF35 interacts with U2AF65 to bind 3' splice site of pre-mRNA and initiates splicing. Another protein, which is named as U2AF1-like3 (U2AF1L3), shows high similarity with U2AF35 and may have related function in pre-mRNA splicing. Here, we report a splice variant of U2AF1L3, which is 767 bp in length and has an open reading frame (ORF) coding a predicted 181 amino acids protein. Reverse transcription-PCR (RT-PCR) shows that this isoform has different expression pattern with U2AF1L3 and is highly expressed in heart, brain and lung.

Isolation and characterization of ubiquitin-activating enzyme E1-domain containing 1, UBE1DC1

Ubiquitin and other ubiquitin-like proteins play important roles in post-translational modification. They are phylogenetically well-conserved in eukaryotes. Activated by other proteins, ubiquitin and ubiquitin-like proteins can covalently modify target proteins. The enzymes responsible for the activation of this modification have been known to include UBA1, SAE2, UBA3, SAE1 and ULA1. Here we report a new ubiquitin activating enzyme like cDNA, named ubiquitin activating enzyme E1-domain containing 1 (UBE1DC1), whose cDNA is 2654 base pairs in length and contains an open reading frame encoding 404 amino acids. The UBE1DC1 gene consists of 12 exons and is located at human chromosome 3q22. The result of RT-PCR showed that UBE1DC1 is expressed in most of human tissues.

Co-evolutionary analysis of insulin/insulin like growth factor 1 signal pathway in vertebrate species

Insulin/insulin like growth factor 1 signaling pathway plays an essential role on the regulation of variant, interrelated and fundamental processes, such as metabolism, growth, reproduction and aging. This pathway is highly conserved during evolution. We analyzed eight vertebrate species' insulin/insulin like growth factor 1 pathway and examined the co-evolutionary relationship between different protein components by quantifying similarity of phylogenetic trees. The collected species include Takifugu rubripes, Danio rerio, Xenopus laevis, Gallus gallus, Mus musculus, Rattous norvegicus, Pan troglodytes and Homo sapiens. Results show that interacting proteins in this pathway share highly co-evolutionary relationship while contrarily the proteins without interaction have low co-evolutionary relationship. We also predict some receptor and ligand partners enjoy highly binding potential in corresponding species' pathway.

Co-evolution analysis on endocrine research: a methodological approach

The rapid growth of different kinds of biological information allows a good opportunity to analyze the co-evolutionary characteristics in endocrine regulatory pathways. Data ranging from kinds of species' genome, gene sequence, protein structure, and expression profile of different organisms can reveal the inner co-evolutionary relationship of ligands, receptors, and other related molecules. In return, these co-evolutionary characteristics can help us determine uncharacterized ligands and receptors, annotate gene functions, highlight amino acid residues with biochemical significance, and identify regulated genes in the endocrine process. Encouraging examples in this field, although at their starting stage, have emerged. Here we focus on recent progress in endocrine-related co-evolution research from a methodological approach.