The P10K database: a data portal for the protist 10 000 genomes project.

Protists, a highly diverse group of microscopic eukaryotic organisms distinct from fungi, animals and plants, exert crucial roles within the earth's biosphere. However, the genomes of only a small fraction of known protist species have been published and made publicly accessible. To address this constraint, the Protist 10 000 Genomes Project (P10K) was initiated, implementing a specialized pipeline for single-cell genome/transcriptome assembly, decontamination and annotation of protists. The resultant P10K database (https://ngdc.cncb.ac.cn/p10k/) serves as a comprehensive platform, collating and disseminating genome sequences and annotations from diverse protist groups. Currently, the P10K database has incorporated 2959 genomes and transcriptomes, including 1101 newly sequenced datasets by P10K and 1858 publicly available datasets. Notably, it covers 45% of the protist orders, with a significant representation (53% coverage) of ciliates, featuring nearly a thousand genomes/transcriptomes. Intriguingly, analysis of the unique codon table usage among ciliates has revealed differences compared to the NCBI taxonomy system, suggesting a need to revise the codon tables used for these species. Collectively, the P10K database serves as a valuable repository of genetic resources for protist research and aims to expand its collection by incorporating more sequenced data and advanced analysis tools to benefit protist studies worldwide.

A vitamin-C-derived DNA modification catalysed by an algal TET homologue.

Methylation of cytosine to 5-methylcytosine (5mC) is a prevalent DNA modification found in many organisms. Sequential oxidation of 5mC by ten-eleven translocation (TET) dioxygenases results in a cascade of additional epigenetic marks and promotes demethylation of DNA in mammals1,2. However, the enzymatic activity and function of TET homologues in other eukaryotes remains largely unexplored. Here we show that the green alga Chlamydomonas reinhardtii contains a 5mC-modifying enzyme (CMD1) that is a TET homologue and catalyses the conjugation of a glyceryl moiety to the methyl group of 5mC through a carbon-carbon bond, resulting in two stereoisomeric nucleobase products. The catalytic activity of CMD1 requires Fe(II) and the integrity of its binding motif His-X-Asp, which is conserved in Fe-dependent dioxygenases3. However, unlike previously described TET enzymes, which use 2-oxoglutarate as a co-substrate4, CMD1 uses L-ascorbic acid (vitamin C) as an essential co-substrate. Vitamin C donates the glyceryl moiety to 5mC with concurrent formation of glyoxylic acid and CO2. The vitamin-C-derived DNA modification is present in the genome of wild-type C. reinhardtii but at a substantially lower level in a CMD1 mutant strain. The fitness of CMD1 mutant cells during exposure to high light levels is reduced. LHCSR3, a gene that is critical for the protection of C. reinhardtii from photo-oxidative damage under high light conditions, is hypermethylated and downregulated in CMD1 mutant cells compared to wild-type cells, causing a reduced capacity for photoprotective non-photochemical quenching. Our study thus identifies a eukaryotic DNA base modification that is catalysed by a divergent TET homologue and unexpectedly derived from vitamin C, and describes its role as a potential epigenetic mark that may counteract DNA methylation in the regulation of photosynthesis.

Ciliogenesis requires sphingolipid-dependent membrane and axoneme interaction.

Cilium formation and regeneration requires new protein synthesis, but the underlying cytosolic translational reprogramming remains largely unknown. Using ribosome footprinting, we performed global translatome profiling during cilia regeneration in Chlamydomonas and uncovered that flagellar genes undergo an early transcriptional activation but late translational repression. This pattern guided our identification of sphingolipid metabolism enzymes, including serine palmitoyltransferase (SPT), as essential regulators for ciliogenesis. Cryo-electron tomography showed that ceramide loss abnormally increased the membrane-axoneme distance and generated bulged cilia. We found that ceramides interact with intraflagellar transport (IFT) particle proteins that IFT motors transport along axoneme microtubules (MTs), suggesting that ceramide-IFT particle-IFT motor-MT interactions connect the ciliary membrane with the axoneme to form rod-shaped cilia. SPT-deficient vertebrate cells were defective in ciliogenesis, and SPT mutations from patients with hereditary sensory neuropathy disrupted cilia, which could be restored by sphingolipid supplementation. These results reveal a conserved role of sphingolipid in cilium formation and link compromised sphingolipid production with ciliopathies.

iDVEIP: A computer-aided approach for the prediction of viral entry inhibitory peptides.

With the notable surge in therapeutic peptide development, various peptides have emerged as potential agents against virus-induced diseases. Viral entry inhibitory peptides (VEIPs), a subset of antiviral peptides (AVPs), offer a promising avenue as entry inhibitors (EIs) with distinct advantages over chemical counterparts. Despite this, a comprehensive analytical platform for characterizing these peptides and their effectiveness in blocking viral entry remains lacking. In this study, we introduce a groundbreaking in silico approach that leverages bioinformatics analysis and machine learning to characterize and identify novel VEIPs. Cross-validation results demonstrate the efficacy of a model combining sequence-based features in predicting VEIPs with high accuracy, validated through independent testing. Additionally, an EI type model has been developed to distinguish peptides specifically acting as Eis from AVPs with alternative activities. Notably, we present iDVEIP, a web-based tool accessible at http://mer.hc.mmh.org.tw/iDVEIP/, designed for automatic analysis and prediction of VEIPs. Emphasizing its capabilities, the tool facilitates comprehensive analyses of peptide characteristics, providing detailed amino acid composition data for each prediction. Furthermore, we showcase the tool's utility in identifying EIs against severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2).

Efficient methods for multiple types of precise gene-editing in Chlamydomonas.

Precise gene-editing using CRISPR/Cas9 technology remains a long-standing challenge, especially for genes with low expression and no selectable phenotypes in Chlamydomonas reinhardtii, a classic model for photosynthesis and cilia research. Here, we developed a multi-type and precise genetic manipulation method in which a DNA break was generated by Cas9 nuclease and the repair was mediated using a homologous DNA template. The efficacy of this method was demonstrated for several types of gene editing, including inactivation of two low-expression genes (CrTET1 and CrKU80), the introduction of a FLAG-HA epitope tag into VIPP1, IFT46, CrTET1 and CrKU80 genes, and placing a YFP tag into VIPP1 and IFT46 for live-cell imaging. We also successfully performed a single amino acid substitution for the FLA3, FLA10 and FTSY genes, and documented the attainment of the anticipated phenotypes. Lastly, we demonstrated that precise fragment deletion from the 3'-UTR of MAA7 and VIPP1 resulted in a stable knock-down effect. Overall, our study has established efficient methods for multiple types of precise gene editing in Chlamydomonas, enabling substitution, insertion and deletion at the base resolution, thus improving the potential of this alga in both basic research and industrial applications.

iDVIP: identification and characterization of viral integrase inhibitory peptides.

Antiretroviral peptides are a kind of bioactive peptides that present inhibitory activity against retroviruses through various mechanisms. Among them, viral integrase inhibitory peptides (VINIPs) are a class of antiretroviral peptides that have the ability to block the action of integrase proteins, which is essential for retroviral replication. As the number of experimentally verified bioactive peptides has increased significantly, the lack of in silico machine learning approaches can effectively predict the peptides with the integrase inhibitory activity. Here, we have developed the first prediction model for identifying the novel VINIPs using the sequence characteristics, and the hybrid feature set was considered to improve the predictive ability. The performance was evaluated by 5-fold cross-validation based on the training dataset, and the result indicates the proposed model is capable of predicting the VINIPs, with a sensitivity of 85.82%, a specificity of 88.81%, an accuracy of 88.37%, a balanced accuracy of 87.32% and a Matthews correlation coefficient value of 0.64. Most importantly, the model also consistently provides effective performance in independent testing. To sum up, we propose the first computational approach for identifying and characterizing the VINIPs, which can be considered novel antiretroviral therapy agents. Ultimately, to facilitate further research and development, iDVIP, an automatic computational tool that predicts the VINIPs has been developed, which is now freely available at http://mer.hc.mmh.org.tw/iDVIP/.

Structural and functional regulation of Chlamydomonas lysosome-related organelles during environmental changes.

Lysosome-related organelles (LROs) are a class of heterogeneous organelles conserved in eukaryotes that primarily play a role in storage and secretion. An important function of LROs is to mediate metal homeostasis. Chlamydomonas reinhardtii is a model organism for studying metal ion metabolism; however, structural and functional analyses of LROs in C. reinhardtii are insufficient. Here, we optimized a method for purifying these organelles from 2 populations of cells: stationary phase or overloaded with iron. The morphology, elemental content, and lysosomal activities differed between the 2 preparations, even though both have phosphorus and metal ion storage functions. LROs in stationary phase cells had multiple non-membrane-bound polyphosphate granules to store phosphorus. Those in iron-overloaded cells were similar to acidocalcisomes (ACs), which have a boundary membrane and contain 1 or 2 large polyphosphate granules to store more phosphorus. We established a method for quantifying the capacity of LROs to sequester individual trace metals. Based on a comparative proteomic analysis of these 2 types of LROs, we present a comprehensive AC proteome and identified 113 putative AC proteins. The methods and protein inventories provide a framework for studying the biogenesis and modification of LROs and the mechanisms by which they participate in regulating metal ion metabolism.

Heme-binding protein CYB5D1 is a radial spoke component required for coordinated ciliary beating.

Coordinated beating is crucial for the function of multiple cilia. However, the molecular mechanism is poorly understood. Here, we characterize a conserved ciliary protein CYB5D1 with a heme-binding domain and a cordon-bleu ubiquitin-like domain. Mutation or knockdown of Cyb5d1 in zebrafish impaired coordinated ciliary beating in the otic vesicle and olfactory epithelium. Similarly, the two flagella of an insertional mutant of the CYB5D1 ortholog in Chlamydomonas (Crcyb5d1) showed an uncoordinated pattern due to a defect in the cis-flagellum. Biochemical analyses revealed that CrCYB5D1 is a radial spoke stalk protein that binds heme only under oxidizing conditions. Lack of CrCYB5D1 resulted in a reductive shift in flagellar redox state and slowing down of the phototactic response. Treatment of Crcyb5d1 with oxidants restored coordinated flagellar beating. Taken together, these data suggest that CrCYB5D1 may integrate environmental and intraciliary signals and regulate the redox state of cilia, which is crucial for the coordinated beating of multiple cilia.

Integrating In Silico and In Vitro Approaches to Identify Natural Peptides with Selective Cytotoxicity against Cancer Cells.

Anticancer peptides (ACPs) are bioactive compounds known for their selective cytotoxicity against tumor cells via various mechanisms. Recent studies have demonstrated that in silico machine learning methods are effective in predicting peptides with anticancer activity. In this study, we collected and analyzed over a thousand experimentally verified ACPs, specifically targeting peptides derived from natural sources. We developed a precise prediction model based on their sequence and structural features, and the model's evaluation results suggest its strong predictive ability for anticancer activity. To enhance reliability, we integrated the results of this model with those from other available methods. In total, we identified 176 potential ACPs, some of which were synthesized and further evaluated using the MTT colorimetric assay. All of these putative ACPs exhibited significant anticancer effects and selective cytotoxicity against specific tumor cells. In summary, we present a strategy for identifying and characterizing natural peptides with selective cytotoxicity against cancer cells, which could serve as novel therapeutic agents. Our prediction model can effectively screen new molecules for potential anticancer activity, and the results from in vitro experiments provide compelling evidence of the candidates' anticancer effects and selective cytotoxicity.

Overexpression of the sulfate transporter-encoding SULTR2 increases chromium accumulation in Chlamydomonas reinhardtii.

Hexavalent chromium [Cr(Ⅵ)] is a highly toxic contaminant in aquatic systems, and microalgae represent promising bioremediators of metal-containing wastewater. However, the metal-binding capacity of algal cells is limited. Therefore, we improved the cellular Cr(Ⅵ) biosorption capacity of Chlamydomonas reinhardtii by overexpressing the sulfate transporter gene SULTR2. SULTR2 was predominantly located in the cytoplasm of the cell, and few proteins mobilized to the cell membrane as a Cr transporter under Cr stress conditions. Intracellular Cr accumulation was almost doubled in SULTR2-overexpressing transgenic strains after exposure to 30 µM K2 Cr2 O7 for 4 d. Alginate-based immobilization increased the rate of Cr removal from 43.81% to 88.15% for SULTR2-overexpressing transgenic strains after exposure to 10 µM K2 Cr2 O7 for 6 d. The immobilized cells also displayed a significant increase in nutrient removal efficiency compared to that of free-swimming cells. Therefore, SULTR2 overexpression in algae has a great potential for the bioremediation of Cr(Ⅵ)-containing wastewater.

Transdermal nanolipoplex simultaneously inhibits subcutaneous melanoma growth and suppresses systemically metastatic melanoma by activating host immunity.

Benefit for clinical melanoma treatments, the transdermal neoadjuvant therapy could reduce surgery region and increase immunotherapy efficacy. Using lipoplex (Lipo-PEG-PEI-complex, LPPC) encapsulated doxorubicin (DOX) and carrying CpG oligodeoxynucleotide; the transdermally administered nano-liposomal drug complex (LPPC-DOX-CpG) would have high cytotoxicity and immunostimulatory activity to suppress systemic metastasis of melanoma. LPPC-DOX-CpG dramatically suppressed subcutaneous melanoma growth by inducing tumor cell apoptosis and recruiting immune cells into the tumor area. Animal studies further showed that the colonization and growth of spontaneously metastatic melanoma cells in the liver and lung were suppressed by transdermal LPPC-DOX-CpG. Furthermore, NGS analysis revealed IFN-γ and NF-κB pathways were triggered to recruit and activate the antigen-presenting-cells and effecter cells, which could activate the anti-tumor responses as the major mechanism responsible for the therapeutic effect of LPPC-DOX-CpG. Finally, we have successfully proved transdermal LPPC-DOX-CpG as a promising penetrative carrier to activate systemic anti-tumor immunity against subcutaneous and metastatic tumor.

Identification and analysis of smORFs in Chlamydomonas reinhardtii.

Small open reading frames (smORFs) have been acknowledged as an important partner in organism functions ranging from bacteria to higher eukaryotes. However, there is a lack of investigation of smORFs in green algae, despite their importance in ecology and evolution. We applied bioinformatic analysis, ribosome profiling, and small peptide proteomics to provide a genome-wide and high-confident smORF database in the model green alga Chlamydomonas reinhardtii. The whole genome was screened first to mine potential coding smORFs. Then conservative analysis, ribosome profiling, and proteomics data were processed to identify conserved smORFs and generate translation evidence. The combination of procedures resulted in 2014 smORFs that might exist in the C. reinhardtii genome. The expression of smORFs in Cd treatment suggested that two smORFs might participate in redox reaction, three in inorganic phosphate transport, and one in DNA repair under stress. Our study built a genome-widely database in C. reinhardtii, providing target smORFs for further research.

The spectrin-based membrane skeleton is asymmetric and remodels during neural development in C. elegans.

Perturbation of spectrin-based membrane mechanics causes hereditary elliptocytosis and spinocerebellar ataxia, but the underlying cellular basis of pathogenesis remains unclear. Here, we introduced conserved disease-associated spectrin mutations into the Caenorhabditis elegans genome and studied the contribution of spectrin to neuronal migration and dendrite formation in developing larvae. The loss of spectrin resulted in ectopic actin polymerization outside of the existing front and secondary membrane protrusions, leading to defective neuronal positioning and dendrite morphology in adult animals. Spectrin accumulated in the lateral region and rear of migrating neuroblasts and redistributes from the soma into the newly formed dendrites, indicating that the spectrin-based membrane skeleton is asymmetric and remodels to regulate actin assembly and cell shape during development. We affinity-purified spectrin from C. elegans and showed that its binding partner ankyrin functions with spectrin. Asymmetry and remodeling of the membrane skeleton might enable spatiotemporal modulation of membrane mechanics for distinct developmental events.

Anticancer peptide Q7 suppresses the growth and migration of human endometrial cancer by inhibiting DHCR24 expression and modulating the AKT-mediated pathway.

Endometrial cancer is one of the most common malignancy affecting women in developed countries. Resection uterus or lesion area is usually the first option for a simple and efficient therapy. Therefore, it is necessary to find a new therapeutic drug to reduce surgery areas to preserve fertility. Anticancer peptides (ACP) are bioactive amino acids with lower toxicity and higher specificity than chemical drugs. This study is to address an ACP, herein named Q7, which could downregulate 24-Dehydrocholesterol Reductase (DHCR24) to disrupt lipid rafts formation, and sequentially affect the AKT signal pathway of HEC-1-A cells to suppress their tumorigenicity such as proliferation and migration. Moreover, lipo-PEI-PEG-complex (LPPC) was used to enhance Q7 anticancer activity in vitro and efficiently show its effects on HEC-1-A cells. Furthermore, LPPC-Q7 exhibited a synergistic effect in combination with doxorubicin or paclitaxel. To summarize, Q7 was firstly proved to exhibit an anticancer effect on endometrial cancer cells and combined with LPPC efficiently improved the cytotoxicity of Q7.

miRTarBase 2020: updates to the experimentally validated microRNA-target interaction database.

MicroRNAs (miRNAs) are small non-coding RNAs (typically consisting of 18-25 nucleotides) that negatively control expression of target genes at the post-transcriptional level. Owing to the biological significance of miRNAs, miRTarBase was developed to provide comprehensive information on experimentally validated miRNA-target interactions (MTIs). To date, the database has accumulated >13,404 validated MTIs from 11,021 articles from manual curations. In this update, a text-mining system was incorporated to enhance the recognition of MTI-related articles by adopting a scoring system. In addition, a variety of biological databases were integrated to provide information on the regulatory network of miRNAs and its expression in blood. Not only targets of miRNAs but also regulators of miRNAs are provided to users for investigating the up- and downstream regulations of miRNAs. Moreover, the number of MTIs with high-throughput experimental evidence increased remarkably (validated by CLIP-seq technology). In conclusion, these improvements promote the miRTarBase as one of the most comprehensively annotated and experimentally validated miRNA-target interaction databases. The updated version of miRTarBase is now available at http://miRTarBase.cuhk.edu.cn/.

Polyubiquitylation of α-tubulin at K304 is required for flagellar disassembly in Chlamydomonas.

Cilia/flagella are structurally conserved and dynamic organelles; their assembly and disassembly are coordinated with the cell cycle and cell differentiation. Several post-translational modifications, including acetylation, methylation, phosphorylation and ubiquitylation, participate in ciliary disassembly. However, the detailed mechanism and the role of ubiquitylation in ciliary disassembly are unclear. This study identified 20 proteins that were ubiquitylated in shortening flagella of Chlamydomonas α-Tubulin was the most abundant ubiquitylated protein and it was labeled with K63 polyubiquitin chains primarily at K304. Expression of an α-tubulin mutant (K304R), which could not be ubiquitylated, decreased the rate of flagellar disassembly and resulted in an enrichment of the mutant form in the axoneme, suggesting that ubiquitylation of α-tubulin is required for the normal kinetics of axonemal disassembly. Immunoprecipitation and glutathione-S-transferase pulldown assays demonstrated that the retrograde intraflagellar transport (IFT) protein, IFT139, interacted with a variety of ubiquitylated proteins, including α-tubulin, suggesting that IFT-A was responsible for transporting ubiquitylated proteins out of the flagella. Our data suggest an important role for ubiquitylation and retrograde IFT in ciliary disassembly.This article has an associated First Person interview with the first author of the paper.

dbAMP: an integrated resource for exploring antimicrobial peptides with functional activities and physicochemical properties on transcriptome and proteome data.

Antimicrobial peptides (AMPs), naturally encoded from genes and generally contained 10-100 amino acids, are crucial components of the innate immune system and can protect the host from various pathogenic bacteria, as well as viruses. In recent years, the widespread use of antibiotics has inspired the rapid growth of antibiotic-resistant microorganisms that usually induce critical infection and pathogenesis. An increasing interest therefore was motivated to explore natural AMPs that enable the development of new antibiotics. With the potential of AMPs being as new drugs for multidrug-resistant pathogens, we were thus motivated to develop a database (dbAMP, http://csb.cse.yzu.edu.tw/dbAMP/) by accumulating comprehensive AMPs from public domain and manually curating literature. Currently in dbAMP there are 12 389 unique entries, including 4271 experimentally verified AMPs and 8118 putative AMPs along with their functional activities, supported by 1924 research articles. The advent of high-throughput biotechnologies, such as mass spectrometry and next-generation sequencing, has led us to further expand dbAMP as a database-assisted platform for providing comprehensively functional and physicochemical analyses for AMPs based on the large-scale transcriptome and proteome data. Significant improvements available in dbAMP include the information of AMP-protein interactions, antimicrobial potency analysis for 'cryptic' region detection, annotations of AMP target species, as well as AMP detection on transcriptome and proteome datasets. Additionally, a Docker container has been developed as a downloadable package for discovering known and novel AMPs on high-throughput omics data. The user-friendly visualization interfaces have been created to facilitate peptide searching, browsing, and sequence alignment against dbAMP entries. All the facilities integrated into dbAMP can promote the functional analyses of AMPs and the discovery of new antimicrobial drugs.

dbPTM in 2019: exploring disease association and cross-talk of post-translational modifications.

The dbPTM (http://dbPTM.mbc.nctu.edu.tw/) has been maintained for over 10 years with the aim to provide functional and structural analyses for post-translational modifications (PTMs). In this update, dbPTM not only integrates more experimentally validated PTMs from available databases and through manual curation of literature but also provides PTM-disease associations based on non-synonymous single nucleotide polymorphisms (nsSNPs). The high-throughput deep sequencing technology has led to a surge in the data generated through analysis of association between SNPs and diseases, both in terms of growth amount and scope. This update thus integrated disease-associated nsSNPs from dbSNP based on genome-wide association studies. The PTM substrate sites located at a specified distance in terms of the amino acids encoded from nsSNPs were deemed to have an association with the involved diseases. In recent years, increasing evidence for crosstalk between PTMs has been reported. Although mass spectrometry-based proteomics has substantially improved our knowledge about substrate site specificity of single PTMs, the fact that the crosstalk of combinatorial PTMs may act in concert with the regulation of protein function and activity is neglected. Because of the relatively limited information about concurrent frequency and functional relevance of PTM crosstalk, in this update, the PTM sites neighboring other PTM sites in a specified window length were subjected to motif discovery and functional enrichment analysis. This update highlights the current challenges in PTM crosstalk investigation and breaks the bottleneck of how proteomics may contribute to understanding PTM codes, revealing the next level of data complexity and proteomic limitation in prospective PTM research.

iDPGK: characterization and identification of lysine phosphoglycerylation sites based on sequence-based features.

BACKGROUND: Protein phosphoglycerylation, the addition of a 1,3-bisphosphoglyceric acid (1,3-BPG) to a lysine residue of a protein and thus to form a 3-phosphoglyceryl-lysine, is a reversible and non-enzymatic post-translational modification (PTM) and plays a regulatory role in glucose metabolism and glycolytic process. As the number of experimentally verified phosphoglycerylated sites has increased significantly, statistical or machine learning methods are imperative for investigating the characteristics of phosphoglycerylation sites. Currently, research into phosphoglycerylation is very limited, and only a few resources are available for the computational identification of phosphoglycerylation sites. RESULT: We present a bioinformatics investigation of phosphoglycerylation sites based on sequence-based features. The TwoSampleLogo analysis reveals that the regions surrounding the phosphoglycerylation sites contain a high relatively of positively charged amino acids, especially in the upstream flanking region. Additionally, the non-polar and aliphatic amino acids are more abundant surrounding phosphoglycerylated lysine following the results of PTM-Logo, which may play a functional role in discriminating between phosphoglycerylation and non-phosphoglycerylation sites. Many types of features were adopted to build the prediction model on the training dataset, including amino acid composition, amino acid pair composition, positional weighted matrix and position-specific scoring matrix. Further, to improve the predictive power, numerous top features ranked by F-score were considered as the final combination for classification, and thus the predictive models were trained using DT, RF and SVM classifiers. Evaluation by five-fold cross-validation showed that the selected features was most effective in discriminating between phosphoglycerylated and non-phosphoglycerylated sites. CONCLUSION: The SVM model trained with the selected sequence-based features performed well, with a sensitivity of 77.5%, a specificity of 73.6%, an accuracy of 74.9%, and a Matthews Correlation Coefficient value of 0.49. Furthermore, the model also consistently provides the effective performance in independent testing set, yielding sensitivity of 75.7% and specificity of 64.9%. Finally, the model has been implemented as a web-based system, namely iDPGK, which is now freely available at http://mer.hc.mmh.org.tw/iDPGK/ .

Genome-wide discovery of viral microRNAs based on phylogenetic analysis and structural evolution of various human papillomavirus subtypes.

In mammals, microRNAs (miRNAs) play key roles in controlling posttranscriptional regulation through binding to the mRNAs of target genes. Recently, it was discovered that viral miRNAs may be involved in human cancers and diseases. It is likely that viral miRNAs help viruses enter the latent phase of their life cycle and become undetected by the host's immune system, while increasing the host's risk for cancer development. Cervical cancer is typically related to the infection of human papillomavirus (HPV) through sexual transmission. To further understand the molecular mechanisms underlying the associations of HPV infection with genital diseases, we developed a systematic method for viral miRNA identification and viral miRNA-mediated regulatory network construction based on genome-wide sequence analysis. The complete genomes of certain high-risk HPV subtypes were used to predict putative viral pre-miRNAs by bioinformatics approaches. In addition, small RNA libraries in human cervical lesions from existing publications were collected to validate the predicted HPV pre-miRNAs. For the construction of virally encoded miRNA-mediated regulatory network of HPV infection, cervical squamous epithelial carcinoma gene expression data were extracted from the RNA sequencing platform in The Cancer Genome Atlas; the differentially expressed genes were used to identify the putative targets of viral miRNAs. Predicted cellular target genes of HPV-encoded miRNAs provide an overview of these viral miRNA's putative functions. Finally, a large-scale genome analysis was carried out to examine the phylogenetic relationship and structural evolution among genital HPV types that have the potential to cause genital cancer. In this study, we discovered putative HPV-encoded miRNAs, which were validated against the small RNA libraries in human cervical lesions. Furthermore, as indicated by their biological functions, host genes targeted by HPV-encoded miRNAs may play significant roles in virus infection and carcinogenesis. These viral miRNAs pose as promising candidates for the development of antiviral drugs. More importantly, the identified subtype-specific miRNAs have the potential to be used as biomarkers for HPV subtype determination.