Liquid drop of DNA libraries reveals total genome information.

Conventional "bulk" PCR often yields inefficient and nonuniform amplification of complex templates in DNA libraries, introducing unwanted biases. Amplification of single DNA molecules encapsulated in a myriad of emulsion droplets (emulsion PCR, ePCR) allows the mitigation of this problem. Different ePCR regimes were experimentally analyzed to identify the most robust techniques for enhanced amplification of DNA libraries. A phenomenological mathematical model that forms an essential basis for optimal use of ePCR for library amplification was developed. A detailed description by high-throughput sequencing of amplified DNA-encoded libraries highlights the principal advantages of ePCR over bulk PCR. ePCR outperforms PCR, reduces gross DNA errors, and provides a more uniform distribution of the amplified sequences. The quasi single-molecule amplification achieved via ePCR represents the fundamental requirement in case of complex DNA templates being prone to diversity degeneration and provides a way to preserve the quality of DNA libraries.

Mutation at the Site of Hydroxylation in the Ribosomal Protein uL15 (RPL27a) Causes Specific Changes in the Repertoire of mRNAs Translated in Mammalian Cells.

Ribosomal protein uL15 (RPL27a) carries a specific modification, hydroxylation, at the His39 residue, which neighbors the CCA terminus of the E-site-bound tRNA at the mammalian ribosome. Under hypoxia, the level of hydroxylation of this protein decreases. We transiently transfected HEK293T cells with constructs expressing wild-type uL15 or mutated uL15 (His39Ala) incapable of hydroxylation, and demonstrated that ribosomes containing both proteins are competent in translation. By applying RNA-seq to the total cellular and polysome-associated mRNAs, we identified differentially expressed genes (DEGs) in cells containing exogenous uL15 or its mutant form. Analyzing mRNA features of up- and down-regulated DEGs, we found an increase in the level of more abundant mRNAs and shorter CDSs in cells with uL15 mutant for both translated and total cellular mRNAs. The level of longer and rarer mRNAs, on the contrary, decreased. Our data show how ribosome heterogeneity can change the composition of the translatome and transcriptome, depending on the properties of the translated mRNAs.

Transcriptomic Analysis of CRISPR/Cas9-Mediated PARP1-Knockout Cells under the Influence of Topotecan and TDP1 Inhibitor.

Topoisomerase 1 (TOP1) is an enzyme that regulates DNA topology and is essential for replication, recombination, and other processes. The normal TOP1 catalytic cycle involves the formation of a short-lived covalent complex with the 3' end of DNA (TOP1 cleavage complex, TOP1cc), which can be stabilized, resulting in cell death. This fact substantiates the effectiveness of anticancer drugs-TOP1 poisons, such as topotecan, that block the relegation of DNA and fix TOP1cc. Tyrosyl-DNA phosphodiesterase 1 (TDP1) is able to eliminate TOP1cc. Thus, TDP1 interferes with the action of topotecan. Poly(ADP-ribose) polymerase 1 (PARP1) is a key regulator of many processes in the cell, such as maintaining the integrity of the genome, regulation of the cell cycle, cell death, and others. PARP1 also controls the repair of TOP1cc. We performed a transcriptomic analysis of wild type and PARP1 knockout HEK293A cells treated with topotecan and TDP1 inhibitor OL9-119 alone and in combination. The largest number of differentially expressed genes (DEGs, about 4000 both up- and down-regulated genes) was found in knockout cells. Topotecan and OL9-119 treatment elicited significantly fewer DEGs in WT cells and negligible DEGs in PARP1-KO cells. A significant part of the changes caused by PARP1-KO affected the synthesis and processing of proteins. Differences under the action of treatment with TOP1 or TDP1 inhibitors alone were found in the signaling pathways for the development of cancer, DNA repair, and the proteasome. The drug combination resulted in DEGs in the ribosome, proteasome, spliceosome, and oxidative phosphorylation pathways.

mRNA regions where 80S ribosomes pause during translation elongation in vivo interact with protein uS19, a component of the decoding site.

In eukaryotic ribosomes, the conserved protein uS19, formerly known as S15, extends with its C-terminal tail to the decoding site. The cross-linking of uS19 to the A site codon has been detected using synthetic mRNAs bearing 4-thiouridine (s4U) residues. Here, we showed that the A-site tRNA prevents this cross-linking and that the P site codon does not contact uS19. Next, we focused on determining uS19-mRNA interactions in vivo by applying the photoactivatable-ribonucleoside enhancing cross-linking and immunoprecipitation method to a stable HEK293 cell line producing FLAG-tagged uS19 and grown in a medium containing s4U. We found that when translation was stopped by cycloheximide, uS19 was efficiently cross-linked to mRNA regions with a high frequency of Glu, Lys and, more rarely, Arg codons. The results indicate that the complexes, in which the A site codon is not involved in the formation of the mRNA-tRNA duplex, are present among the cycloheximide-arrested 80S complexes, which implies pausing of elongating ribosomes at the above mRNA regions. Thus, our findings demonstrate that the human ribosomal protein uS19 interacts with mRNAs during translation elongation and highlight the regions of mRNAs where ribosome pausing occurs, bringing new structural and functional insights into eukaryotic translation in vivo.

Changes in the Transcriptome Caused by Mutations in the Ribosomal Protein uS10 Associated with a Predisposition to Colorectal Cancer.

A number of mutations in the RPS20 gene encoding the ribosomal protein uS10 have been found to be associated with a predisposition to hereditary non-polyposis colorectal carcinoma (CRC). We transfected HEK293T cells with constructs carrying the uS10 minigene with mutations identical to those mentioned above and examined the effects of the produced proteins on the cellular transcriptome. We showed that uS10 with mutations p.V50SfsX23 or p.L61EfsX11 cannot be incorporated into 40S ribosomal subunits, while the protein with the missense mutation p.V54L functionally replaces the respective endogenous protein in the 40S subunit assembly and the translation process. The comparison of RNA-seq data obtained from cells producing aberrant forms of uS10 with data for those producing the wild-type protein revealed overlapping sets of upregulated and downregulated differently expressed genes (DEGs) related to several pathways. Among the limited number of upregulated DEGs, there were genes directly associated with the progression of CRC, e.g., PPM1D and PIGN. Our findings indicate that the accumulation of the mutant forms of uS10 triggers a cascade of cellular events, similar to that which is triggered when the cell responds to a large number of erroneous proteins, suggesting that this may increase the risk of cancer.

Reorganization of the Landscape of Translated mRNAs in NSUN2-Deficient Cells and Specific Features of NSUN2 Target mRNAs.

The RNA cytosine C5 methyltransferase NSUN2 has a variety of RNA substrates and plays an important role in mRNA metabolism. NSUN2 binds to specific sequences enriched in exosomal mRNAs, suggesting its possible involvement in the sorting of mRNAs into exosomes. We applied the photoactivatable.4-thiouridine-enhanced cross-linking and immunoprecipitation assay involving high-throughput RNA sequencing (RNA-seq) to HEK293T cells to determine NSUN2 mRNA targets. NSUN2 cross-linking sites were found in more than one hundred relatively abundant mRNAs with a high GC content and a pronounced secondary structure. Then, utilizing RNA-seq for the total and polysome-associated mRNA from HEK293T cells with and without the knockdown of NSUN2, we identified differentially expressed genes, as well as genes with altered translational efficiency (GATEs). It turned out that the up-regulated GATE mRNAs were much shorter on average than the down-regulated ones, and their GC content was higher; moreover, they contained motifs with C residues located in GC-rich environments. Our findings reveal the specific features of mRNAs that make them potential targets for NSUN2 and expand our understanding of the role of NSUN2 in controlling translation and, possibly, in mRNA sorting into exosomes implemented through the methylation of cytosine residues.

The Cytogenetic Map of the Nile Crocodile (Crocodylus niloticus, Crocodylidae, Reptilia) with Fluorescence In Situ Localization of Major Repetitive DNAs.

Tandemly arranged and dispersed repetitive DNA sequences are important structural and functional elements that make up a significant portion of vertebrate genomes. Using high throughput, low coverage whole genome sequencing followed by bioinformatics analysis, we have identified seven major tandem repetitive DNAs and two fragments of LTR retrotransposons in the genome of the Nile crocodile (Crocodylus niloticus, 2n = 32). The repeats showed great variability in structure, genomic organization, and chromosomal distribution as revealed by fluorescence in situ hybridization (FISH). We found that centromeric and pericentromeric heterochromatin of C. niloticus is composed of previously described in Crocodylus siamensis CSI-HindIII and CSI-DraI repetitive sequence families, a satellite revealed in Crocodylus porosus, and additionally contains at least three previously unannotated tandem repeats. Both LTR sequences identified here belong to the ERV1 family of endogenous retroviruses. Each pericentromeric region was characterized by a diverse set of repeats, with the exception of chromosome pair 4, in which we found only one type of satellite. Only a few repeats showed non-centromeric signals in addition to their centromeric localization. Mapping of 18S-28S ribosomal RNA genes and telomeric sequences (TTAGGG)n did not demonstrate any co-localization of these sequences with revealed centromeric and pericentromeric heterochromatic blocks.

Flow-Seq Evaluation of Translation Driven by a Set of Natural Escherichia coli 5'-UTR of Variable Length.

Flow-seq is a method that combines fluorescently activated cell sorting and next-generation sequencing to deduce a large amount of data about translation efficiency from a single experiment. Here, we constructed a library of fluorescent protein-based reporters preceded by a set of 648 natural 5'-untranslated regions (5'-UTRs) of Escherichia coli genes. Usually, Flow-seq libraries are constructed using uniform-length sequence elements, in contrast to natural situations, where functional elements are of heterogenous lengths. Here, we demonstrated that a 5'-UTR library of variable length could be created and analyzed with Flow-seq. In line with previous Flow-seq experiments with randomized 5'-UTRs, we observed the influence of an RNA secondary structure and Shine-Dalgarno sequences on translation efficiency; however, the variability of these parameters for natural 5'-UTRs in our library was smaller in comparison with randomized libraries. In line with this, we only observed a 30-fold difference in translation efficiency between the best and worst bins sorted with this factor. The results correlated with those obtained with ribosome profiling.

Deficiency of the Ribosomal Protein uL5 Leads to Significant Rearrangements of the Transcriptional and Translational Landscapes in Mammalian Cells.

Protein uL5 (formerly called L11) is an integral component of the large (60S) subunit of the human ribosome, and its deficiency in cells leads to the impaired biogenesis of 60S subunits. Using RNA interference, we reduced the level of uL5 in HEK293T cells by three times, which caused an almost proportional decrease in the content of the fraction corresponding to 80S ribosomes, without a noticeable diminution in the level of polysomes. By RNA sequencing of uL5-deficient and control cell samples, which were those of total mRNA and mRNA from the polysome fraction, we identified hundreds of differentially expressed genes (DEGs) at the transcriptome and translatome levels and revealed dozens of genes with altered translational efficiency (GATEs). Transcriptionally up-regulated DEGs were mainly associated with rRNA processing, pre-mRNA splicing, translation and DNA repair, while down-regulated DEGs were genes of membrane proteins; the type of regulation depended on the GC content in the 3' untranslated regions of DEG mRNAs. The belonging of GATEs to up-regulated and down-regulated ones was determined by the coding sequence length of their mRNAs. Our findings suggest that the effects observed in uL5-deficient cells result from an insufficiency of translationally active ribosomes caused by a deficiency of 60S subunits.

Knockdown of the Ribosomal Protein eL38 in HEK293 Cells Changes the Translational Efficiency of Specific Genes.

The protein eL38 is one of the smallest proteins of the mammalian ribosome, which is a component of its large (60S) subunit. The haploinsufficiency of eL38 in mice leads to the Tail-short mutant phenotype characterized by defects in the development of the axial skeleton caused by the poor translation of mRNA subsets of Hox genes. Using the ribosome profiling assay applied to HEK293 cells knocked down of eL38, we examined the effects of the lack of eL38 in 60S subunits on gene expression at the level of translation. A four-fold decrease in the cell content of eL38 was shown to result in significant changes in the translational efficiencies of 150 genes. Among the genes, whose expression at the level of translation was enhanced, there were mainly those associated with basic metabolic processes; namely, translation, protein folding, chromosome organization, splicing, and others. The set of genes with reduced translation efficiencies contained those that are mostly involved in the processes related to the regulation of transcription, including the activation of Hox genes. Thus, we demonstrated that eL38 insufficiency significantly affects the expression of certain genes at the translational level. Our findings facilitate understanding the possible causes of some anomalies in eL38-deficient animals.

Human ribosomal protein eS1 is engaged in cellular events related to processing and functioning of U11 snRNA.

Ribosomal proteins are involved in many cellular processes through interactions with various RNAs. Here, applying the photoactivatable-ribonucleoside-enhanced cross-linking and immunoprecipitation approach to HEK293 cells overproducing ribosomal protein (rp) eS1, we determined the products of RNU5A-1 and RNU11 genes encoding U5 and U11 snRNAs as the RNA partners of ribosome-unbound rp eS1. U11 pre-snRNA-associated rp eS1 was revealed in the cytoplasm and nucleus where rp eS1-bound U11/U12 di-snRNP was also found. Utilizing recombinant rp eS1 and 4-thiouridine-containing U11 snRNA transcript, we identified an N-terminal peptide contacting the U-rich sequence in the Sm site-containing RNA region. We also showed that the rp eS1 binding site on U11 snRNA is located in the cleft between stem-loops I and III and that its structure mimics the respective site on the 18S rRNA. It was found that cell depletion of rp eS1 leads to a decrease in the splicing efficiency of minor introns and to an increase in the level of U11 pre-snRNA with the unprocessed 3' terminus. Our findings demonstrate the engagement of human rp eS1 in events related to the U11 snRNA processing and to minor-class splicing. Contacts of rp eS1 with U5 snRNA in the minor pre-catalytic spliceosome are discussed.

Recognition but no repair of abasic site in single-stranded DNA by human ribosomal uS3 protein residing within intact 40S subunit.

Isolated human ribosomal protein uS3 has extra-ribosomal functions including those related to base excision DNA repair, e.g. AP lyase activity that nicks double-stranded (ds) DNA 3΄ to the abasic (AP) site. However, the ability of uS3 residing within ribosome to recognize and cleave damaged DNA has never been addressed. Here, we compare interactions of single-stranded (ss) DNA and dsDNA bearing AP site with human ribosome-bound uS3 and with the isolated protein, whose interactions with ssDNA were not yet studied. The AP lyase activity of free uS3 was much higher with ssDNA than with dsDNA, whereas ribosome-bound uS3 was completely deprived of this activity. Nevertheless, an exposed peptide of ribosome-bound uS3 located far away from the putative catalytic center previously suggested for isolated uS3 cross-linked to full-length uncleaved ssDNA, but not to dsDNA. In contrast, free uS3 cross-linked mainly to the 5΄-part of the damaged DNA strand after its cleavage at the AP site. ChIP-seq analysis showed preferential uS3 binding to nucleolus-associated chromatin domains. We conclude that free and ribosome-bound uS3 proteins interact with AP sites differently, exhibiting their non-translational functions in DNA repair in and around the nucleolus and in regulation of DNA damage response in looped DNA structures, respectively.

Cytosolic YB-1 and NSUN2 are the only proteins recognizing specific motifs present in mRNAs enriched in exosomes.

Exosomes, membranous vesicles secreted by various cells, are involved in intercellular communication and carry vast repertoires of RNAs and proteins. Processes mediating RNA sorting into exosomes are currently poorly understood. Using bioinformatics approaches, three structural motifs ACCAGCCU, CAGUGAGC and UAAUCCCA have been discovered as enriched in exosomal mRNAs and long noncoding RNAs. Here, utilizing short RNA hairpins, each containing one of the motifs, in a pull-down assay of cytosolic extract of human embryonic kidney 293 (HEK293) cells, we prove that multifunctional RNA-binding protein YB-1 specifically interacts with all three motifs, whereas methyltransferase NSUN2 recognizes only the motif CAGUGAGC. RNA hairpins other than those mentioned above pull out neither YB-1 nor NSUN2. Both these proteins are found in exosomes secreted by HEK293 cells. YB-1 for all that is detected as a form having a slightly higher electrophoretic mobility than that of YB-1 associated with the above RNA hairpins, assuming changes in posttranslational modifications of the protein during its transfer from cytoplasm into exosomes. Next generation sequencing of total exosomal RNA (eRNA) reveals a large representative set of RNA species, including mRNAs containing the above-mentioned motifs. The degree of enrichment in exosomes with this kind of mRNAs strongly depends on the locations of eRNA-specific motifs within the mRNA sequences. Altogether, our findings point to YB-1 and NSUN2 as possible mediators of the process of transfer of specific mRNAs into exosomes, allowing us to speculate on an involvement of these proteins in the mRNA sorting via the recognition of the above motifs.

Evolutionary dynamics of Anolis sex chromosomes revealed by sequencing of flow sorting-derived microchromosome-specific DNA.

Squamate reptiles show a striking diversity in modes of sex determination, including both genetic (XY or ZW) and temperature-dependent sex determination systems. The genomes of only a handful of species have been sequenced, analyzed and assembled including the genome of Anolis carolinensis. Despite a high genome coverage, only macrochromosomes of A. carolinensis were assembled whereas the content of most microchromosomes remained unclear. Most of the Anolis species have homomorphic XY sex chromosome system. However, some species have large heteromorphic XY chromosomes (e.g., A. sagrei) and even multiple sex chromosomes systems (e.g. A. pogus), that were shown to be derived from fusions of the ancestral XY with microautosomes. We applied next generation sequencing of flow sorting-derived chromosome-specific DNA pools to characterize the content and composition of microchromosomes in A. carolinensis and A. sagrei. Comparative analysis of sequenced chromosome-specific DNA pools revealed that the A. sagrei XY sex chromosomes contain regions homologous to several microautosomes of A. carolinensis. We suggest that the sex chromosomes of A. sagrei are derived by fusions of the ancestral sex chromosome with three microautosomes and subsequent loss of some genetic content on the Y chromosome.

Deficiency of the ribosomal protein uS10 (RPS20) reorganizes human cells translatome according to the abundance, CDS length and GC content of mRNAs.

Ribosomal protein uS10, a product of the RPS20 gene, is an essential constituent of the small (40S) subunit of the human ribosome. Disruptive mutations in its gene are associated with a predisposition to hereditary colorectal carcinoma. Here, using HEK293T cells, we show that a deficiency of this protein leads to a decrease in the level of ribosomes (ribosomal shortage). RNA sequencing of the total and polysome-associated mRNA samples reveals hundreds of genes differentially expressed in the transcriptome (t)DEGs and translatome (p)DEGs under conditions of uS10 deficiency. We demonstrate that the (t)DEG and (p)DEG sets partially overlap, determine genes with altered translational efficiency (TE) and identify cellular processes affected by uS10 deficiency-induced ribosomal shortage. We reveal that translated mRNAs of upregulated (p)DEGs and genes with altered TE in uS10-deficient cells are generally more abundant and that their GC contents are significantly lower than those of the respective downregulated sets. We also observed that upregulated (p)DEGs have longer coding sequences. Based on our findings, we propose a combinatorial model describing the process of reorganization of mRNA translation under conditions of ribosomal shortage. Our results reveal rules according to which ribosomal shortage reorganizes the transcriptome and translatome repertoires of actively proliferating cells.

Bioluminescent aptamer-based microassay for detection of melanoma inhibitory activity protein (MIA).

Melanoma inhibitory activity protein (MIA) does obviously offer the potential to reveal clinical manifestations of melanoma. Despite a pressing need for effective diagnosis of this highly fatal disease, there are no clinically approved MIA detection ELISA kits available. A recommended MIA threshold has not yet been defined, mostly by reason of variability in immunoglobulins' affinity and stability, the difference in sample preparation and assay conditions. Here we present a pair of high-affinity DNA aptamers developed as an alternative recognition and binding element for MIA detection. Their stability and reproducible synthesis are expected to ensure this analysis under standard conditions. The devised aptamer-based solid-phase microassay of model standard and control human sera involves luciferase NLuc as a highly sensitive reporter. Bioluminescence dependence on MIA concentration ranges in a linear manner from 2.5 to 250 ng mL-1, providing a MIA detection limit of 1.67 ± 0.57 ng mL-1.

Genetic Diversity of the Human Adenovirus C Isolated from Hospitalized Children in Russia (2019-2022).

The human adenovirus (HAdV) is a common pathogen in children that can cause acute respiratory virus infection (ARVI). However, the molecular epidemiological and clinical information relating to HAdV among hospitalized children with ARVI is rarely reported in Russia. A 4-year longitudinal (2019-2022) study among hospitalized children (0-17 years old) with ARVI in Novosibirsk, Russia, was conducted to evaluate the epidemiological and molecular characteristics of HAdV. Statistically significant differences in the detection rates of epidemiological and virological data of all positive viral detections of HAdV were analyzed using a two-tailed Chi-square test. The incidence of HAdV and other respiratory viruses such as human influenza A and B viruses, respiratory syncytial virus, coronavirus, parainfluenza virus, metapneumovirus, rhinovirus, bocavirus, and SARS-CoV-2 was investigated among 3190 hospitalized children using real-time polymerase chain reaction. At least one of these respiratory viruses was detected in 74.4% of hospitalized cases, among which HAdV accounted for 4%. A total of 1.3% co-infections with HAdV were also registered. We obtained full-genome sequences of 12 HAdVs, which were isolated in cell cultures. Genetic analysis revealed the circulation of adenovirus of genotypes C1, C2, C5, C89, and 108 among hospitalized children in the period from 2019-2022.

Deconvolution of B cell receptor repertoire in multiple sclerosis patients revealed a delay in tBreg maturation.

Background: B lymphocytes play a pivotal regulatory role in the development of the immune response. It was previously shown that deficiency in B regulatory cells (Bregs) or a decrease in their anti-inflammatory activity can lead to immunological dysfunctions. However, the exact mechanisms of Bregs development and functioning are only partially resolved. For instance, only a little is known about the structure of their B cell receptor (BCR) repertoires in autoimmune disorders, including multiple sclerosis (MS), a severe neuroinflammatory disease with a yet unknown etiology. Here, we elucidate specific properties of B regulatory cells in MS. Methods: We performed a prospective study of the transitional Breg (tBreg) subpopulations with the CD19+CD24highCD38high phenotype from MS patients and healthy donors by (i) measuring their content during two diverging courses of relapsing-remitting MS: benign multiple sclerosis (BMS) and highly active multiple sclerosis (HAMS); (ii) analyzing BCR repertoires of circulating B cells by high-throughput sequencing; and (iii) measuring the percentage of CD27+ cells in tBregs. Results: The tBregs from HAMS patients carry the heavy chain with a lower amount of hypermutations than tBregs from healthy donors. The percentage of transitional CD24highCD38high B cells is elevated, whereas the frequency of differentiated CD27+ cells in this transitional B cell subset was decreased in the MS patients as compared with healthy donors. Conclusions: Impaired maturation of regulatory B cells is associated with MS progression.

Extended Evaluation of Viral Diversity in Lake Baikal through Metagenomics.

Lake Baikal is a unique oligotrophic freshwater lake with unusually cold conditions and amazing biological diversity. Studies of the lake's viral communities have begun recently, and their full diversity is not elucidated yet. Here, we performed DNA viral metagenomic analysis on integral samples from four different deep-water and shallow stations of the southern and central basins of the lake. There was a strict distinction of viral communities in areas with different environmental conditions. Comparative analysis with other freshwater lakes revealed the highest similarity of Baikal viromes with those of the Asian lakes Soyang and Biwa. Analysis of new data, together with previously published data allowed us to get a deeper insight into the diversity and functional potential of Baikal viruses; however, the true diversity of Baikal viruses in the lake ecosystem remains still unknown. The new metaviromic data will be useful for future studies of viral composition, distribution, and the dynamics associated with global climatic and anthropogenic impacts on this ecosystem.

Gene Analysis, Cloning, and Heterologous Expression of Protease from a Micromycete Aspergillus ochraceus Capable of Activating Protein C of Blood Plasma.

Micromycetes are known to secrete numerous enzymes of biotechnological and medical potential. Fibrinolytic protease-activator of protein C (PAPC) of blood plasma from micromycete Aspergillus ochraceus VKM-F4104D was obtained in recombinant form utilising the bacterial expression system. This enzyme, which belongs to the proteinase-K-like proteases, is similar to the proteases encoded in the genomes of Aspergillus fumigatus ATCC MYA-4609, A. oryzae ATCC 42149 and A. flavus 28. Mature PAPC-4104 is 282 amino acids long, preceded by the 101-amino acid propeptide necessary for proper folding and maturation. The recombinant protease was identical to the native enzyme from micromycete in terms of its biological properties, including an ability to hydrolyse substrates of activated protein C (pGlu-Pro-Arg-pNA) and factor Xa (Z-D-Arg-Gly-Arg-pNA) in conjugant reactions with human blood plasma. Therefore, recombinant PAPC-4104 can potentially be used in medicine, veterinary science, diagnostics, and other applications.