Dynamics of binding ability prediction between spike protein and human ACE2 reveals the adaptive strategy of SARS-CoV-2 in humans.

SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) is a novel coronavirus causing the COVID-19 pandemic in 2020. High adaptive plasticity on the spike protein of SASR-CoV-2 enables it to transmit across different host species. In the present study, we collected 2092 high-quality genome sequences of SARS-CoV-2 from 160 regions in over 50 countries and reconstructed their phylogeny. We also analyzed the polymorphic interaction between spike protein and human ACE2 (hACE2). Phylogenetic analysis of SARS-CoV-2 suggests that SARS-CoV-2 is probably originated from a recombination event on the spike protein between a bat coronavirus and a pangolin coronavirus that endows it humans infectivity. Compared with other regions in the S gene of SARS-CoV-2, the direct-binding sites of the receptor-binding domain (RBD) is more conserved. We focused on 3,860 amino acid mutations in spike protein RBD (T333-C525) of SARS-CoV-2 and simulated their differential stability and binding affinity to hACE2 (S19-D615). The results indicate no preference for SARS-CoV-2 infectivity on people of different ethnic groups. The variants in the spike protein of SARS-CoV-2 may also be a good indicator demonstrating the transmission route of SARS-CoV-2 from its natural reservoir to human hosts.

Genetic Epidemiology of Medication Safety and Efficacy Related Variants in the Central Han Chinese Population With Whole Genome Sequencing.

Medication safety and efficacy-related pharmacogenomic research play a critical role in precision medicine. This study comprehensively analyzed the pharmacogenomic profiles of the central Han Chinese population in the context of medication safety and efficacy and compared them with other global populations. The ultimate goal is to improve medical treatment guidelines. We performed whole-genome sequencing in 487 Han Chinese individuals and investigated the allele frequencies of pharmacogenetic variants in 1,731 drug response-related genes. We identified 2,139 (81.18%) previously reported variants in our population with annotations in the PharmGKB database. The allele frequencies of these 2,139 clinical-related variants were similar to those in other East Asian populations but different from those in other global populations. We predicted the functional effects of nonsynonymous variants in the 1,731 pharmacogenes and identified 1,281 novel and 4,442 previously reported deleterious variants. Of the 1,281 novel deleterious variants, five are common variants with an allele frequency >5%, and the rest are rare variants with an allele frequency <5%. Of the 4,442 known deleterious variants, the allele frequencies were found to differ from those in other populations, of which 146 are common variants. In addition, we found many variants in non-coding regions, the functions of which require further investigation. This study compiled a large amount of data on pharmacogenomic variants in the central Han Chinese population. At the same time, it provides insight into the role of pharmacogenomic variants in clinical medication safety and efficacy.

Genome-wide characterization of Toll-like receptor gene family in common carp (Cyprinus carpio) and their involvement in host immune response to Aeromonas hydrophila infection.

The Toll-like receptor (TLR) gene family is a class of conserved pattern recognition receptors, which play an essential role in innate immunity providing efficient defense against invading microbial pathogens. Although TLRs have been extensively characterized in both invertebrates and vertebrates, a comprehensive analysis of TLRs in common carp is lacking. In the present study, we have conducted the first genome-wide systematic analysis of common carp (Cyprinus carpio) TLR genes. A set of 27 common carp TLR genes were identified and characterized. Sequence similarity analysis, functional domain prediction and phylogenetic analysis supported their annotation and orthologies. By examining the gene copy number of TLR genes across several vertebrates, gene duplications and losses were observed. The expression patterns of TLR genes were examined during early developmental stages and in various healthy tissues, and the results showed that TLR genes were ubiquitously expressed, indicating a likely role in maintaining homeostasis. Moreover, the differential expression of TLRs was examined after Aeromons hydrophila infection, and showed that most TLR genes were induced, with diverse patterns. TLR1, TLR4-2, TLR4-3, TLR22-2, TLR22-3 were significantly up-regulated at minimum one timepoint, whereas TLR2-1, TLR4-1, TLR7-1 and TLR7-2 were significantly down-regulated. Our results suggested that TLR genes play critical roles in the common carp immune response. Collectively, our findings provide fundamental genomic resources for future studies on fish disease management and disease-resistance selective breeding strategy development.

Transcriptome signatures in common carp spleen in response to Aeromonas hydrophila infection.

The common carp is an important aquaculture species that is worldwide distributed. Nowadays, intensive rearing in aquaculture increases the susceptibility of fish to various pathogens such as Aeromonas hydrophila, which has caused severe damage to carp production. However, systematic analysis on the host response of common carp against A. hydrophila is less studied. In order to better understand the common carp immune response process against bacteria at the global gene expression level, we examined transcriptional profiles of the common carp spleen at three timepoints following experimental infection with A. hydrophila. A total of 545 million 125-bp paired end reads were generated, and all trimmed clean reads were mapped onto the common carp whole genome sequence. Comparison of the transcriptomes between the treatment and control group fish revealed 2900 unigenes with significantly differential expression, including 732, 936, 928 genes up-regulated, and 248, 475, 700 genes down-regulated at 4 h, 12 h, 24 h post infection respectively. The captured significantly differentially expressed genes are mainly involved in the pathways including junction/adhesion, pathogen recognition, cell surface receptor signaling, and immune system process/defense response. Our study will provide fundamental information on molecular mechanism underlying the immune response of teleost against bacterial infection and might suggest strategies for selection of resistant strains of common carp in aquaculture.

Comparative transcriptome analysis between aquatic and aerial breathing organs of Channa argus to reveal the genetic basis underlying bimodal respiration.

Aerial breathing in fish was an important adaption for successful survival in hypoxic water. All aerial breathing fish are bimodal breathers. It is intriguing that they can obtain oxygen from both air and water. However, the genetic basis underlying bimodal breathing has not been extensively studied. In this study, we performed next-generation sequencing on a bimodal breathing fish, the Northern snakehead, Channa argus, and generated a transcriptome profiling of C. argus. A total of 53,591 microsatellites and 26,378 SNPs were identified and classified. A Ka/Ks analysis of the unigenes indicated that 63 genes were under strong positive selection. Furthermore, the transcriptomes from the aquatic breathing organ (gill) and the aerial breathing organ (suprabranchial chamber) were sequenced and compared, and the results showed 1,966 genes up-regulated in the gill and 2,727 genes up-regulated in the suprabranchial chamber. A gene pathway analysis concluded that four functional categories were significant, of which angiogenesis and elastic fibre formation were up-regulated in the suprabranchial chamber, indicating that the aerial breathing organ may be more efficient for gas exchange due to its highly vascularized and elastic structure. In contrast, ion uptake and transport and acid-base balance were up-regulated in the gill, indicating that the aquatic breathing organ functions in ion homeostasis and acid-base balance, in addition to breathing. Understanding the genetic mechanism underlying bimodal breathing will shed light on the initiation and importance of aerial breathing in the evolution of vertebrates.

Genome wide identification of scavenger receptors class A in common carp (Cyprinus carpio) and their expression following Aeromonas hydrophila infection.

Scavenger receptors class A (SCARAs) is a subgroup of diverse families of pattern recognition receptors that bind a range of ligands, and play important roles in innate immune processes through pathogens detection, adhesion, endocytosis, and phagocytosis. However, most studies of SCARAs have focused on mammals, and much less is known of SCARAs in fish species. In this study, we identified 7 SCARAs across the common carp genome, which were classified into four subclasses according to comparative genomic analysis including sequence similarities analysis, gene structure and functional domain prediction. Further phylogenetic and syntenic analysis supported their annotation and orthologies. Through examining gene copy number of SCARA genes across several vertebrates, SCARA2, SCARA3 and SCARA4 were found have undergone gene duplication. The expression patterns of SCARAs in common carp were examined during early developmental stages, in healthy tissues, and after Aeromonas hydrophila infection. Most SCARA genes were ubiquitously expressed during common carp early developmental stages, and presented diverse patterns in various healthy tissues, with relatively high expression levels in spleen, liver, intestine, gill and brain, indicating their critical roles likely in maintaining homeostasis and host immune response activities. After A. hydrophila infection, most SCARA genes were up-regulated at 4 h post infection in mucosal tissue intestine, while generally up-regulated at 12 h post infection in spleen, suggesting a tissue-specific pattern of regulation. Taken together, all these results suggested that SCARA genes played important roles in host immune response to A. hydrophila infection in common carp, and provided important genomic resources for future studies on fish disease management.

Genome-Wide Identification, Characterization and Phylogenetic Analysis of ATP-Binding Cassette (ABC) Transporter Genes in Common Carp (Cyprinus carpio).

The ATP-binding cassette (ABC) gene family is considered to be one of the largest gene families in all forms of prokaryotic and eukaryotic life. Although the ABC transporter genes have been annotated in some species, detailed information about the ABC superfamily and the evolutionary characterization of ABC genes in common carp (Cyprinus carpio) are still unclear. In this research, we identified 61 ABC transporter genes in the common carp genome. Phylogenetic analysis revealed that they could be classified into seven subfamilies, namely 11 ABCAs, six ABCBs, 19 ABCCs, eight ABCDs, two ABCEs, four ABCFs, and 11 ABCGs. Comparative analysis of the ABC genes in seven vertebrate species including common carp, showed that at least 10 common carp genes were retained from the third round of whole genome duplication, while 12 duplicated ABC genes may have come from the fourth round of whole genome duplication. Gene losses were also observed for 14 ABC genes. Expression profiles of the 61 ABC genes in six common carp tissues (brain, heart, spleen, kidney, intestine, and gill) revealed extensive functional divergence among the ABC genes. Different copies of some genes had tissue-specific expression patterns, which may indicate some gene function specialization. This study provides essential genomic resources for future studies in common carp.

Genome wide identification, phylogeny and expression of zinc transporter genes in common carp.

BACKGROUND: Zinc is an essential trace element in organisms, which serves as a cofactor for hundreds of enzymes that are involved in many pivotal biological processes including growth, development, reproduction and immunity. Therefore, the homeostasis of zinc in the cell is fundamental. The zinc transporter gene family is a large gene family that encodes proteins which regulate the movement of zinc across cellular and intracellular membranes. However, studies on teleost zinc transporters are mainly limited to model species. METHODOLOGY/PRINCIPAL FINDINGS: We identified a set of 37 zinc transporters in common carp genome, including 17 from SLC30 family (ZnT), and 20 from SLC39 family (ZIP). Phylogenetic and syntenic analysis revealed that most of the zinc transporters are highly conserved, though recent gene duplication and gene losses do exist. Through examining the copy number of zinc transporter genes across several vertebrate genomes, thirteen zinc transporters in common carp are found to have undergone the gene duplications, including SLC30A1, SLC30A2, SLC30A5, SLC30A7, SLC30A9, SLC30A10, SLC39A1, SLC39A3, SLC39A4, SLC39A5, SLC39A6, SLC39A7 and SLC39A9. The expression patterns of all zinc transporters were established in various tissues, including blood, brain, gill, heart, intestine, liver, muscle, skin, spleen and kidney, and showed that most of the zinc transporters were ubiquitously expressed, indicating the critical role of zinc transporters in common carp. CONCLUSIONS: To some extent, examination of gene families with detailed phylogenetic or orthology analysis could verify the authenticity and accuracy of assembly and annotation of the recently published common carp whole genome sequences. The gene families are also considered as a unique source for evolutionary studies. Moreover, the whole set of common carp zinc transporters provides an important genomic resource for future biochemical, toxicological and physiological studies of zinc in teleost.