A CitationRank algorithm inheriting Google technology designed to highlight genes responsible for serious adverse drug reaction

DeepSide: A Deep Learning Approach for Drug Side Effect Prediction

Drug failures due to unforeseen adverse effects at clinical trials pose health risks for the participants and lead to substantial financial losses. Side effect prediction algorithms have the potential to guide the drug design process. LINCS L1000 dataset provides a vast resource of cell line gene expression data perturbed by different drugs and creates a knowledge base for context specific features. The state-of-the-art approach that aims at using context specific information relies on only the high-quality experiments in LINCS L1000 and discards a large portion of the experiments. In this study, our goal is to boost the prediction performance by utilizing this data to its full extent. We experiment with 5 deep learning architectures. We find that a multi-modal architecture produces the best predictive performance among multi-layer perceptron-based architectures when drug chemical structure (CS), and the full set of drug perturbed gene expression profiles (GEX) are used as modalities. Overall, we observe that the CS is more informative than the GEX. A convolutional neural network-based model that uses only SMILES string representation of the drugs achieves the best results and provides 13.0% macro-AUC and 3.1% micro-AUC improvements over the state-of-the-art. We also show that the model is able to predict side effect-drug pairs that are reported in the literature but was missing in the ground truth side effect dataset. DeepSide is available at http://github.com/OnurUner/DeepSide.

Aqueous extract from Madhuca indica bark protects cells from oxidative stress caused by electron beam radiation: in vitro, in vivo and in silico approach

In an endeavor to find the novel natural radioprotector to secure normal cells surrounding cancerous cell during radiation exposure, Madhuca indica (M. indica) aqueous stem bark extract was evaluated for radioprotective activity using in vitro, in vivo, and in silico models. M. indica extract exhibited concentration dependent protective effect on electron beam radiation (EBR) induced damage to pBR322 DNA; the highest protection was achieved at 150 μg concentrations. Similarly, M. indica extract (400 mg/kg) administrated to mice prior to irradiation protected DNA from the radiation damage, which was confirmed by inhibiting comet parameters. The study showed a significant increase in the levels of glutathione and superoxide dismutase levels. The study also revealed that administration of M. Indica at the different dose to mice significantly reduced EBR induced MDA, sialic acid and nitric acid levels. Further extract prevented histophatological changes of skin and liver. In contrast, protein-protein interaction studies were performed to find the hub protein, involved in radiation-induced DNA damage. Among 437 proteins that are found expressed during radiation, p53 was found to be a master protein regulating the whole pathway. Molecular interaction between p53 and M. indica extract was predicted by quantitative structure-activity relationship and ADMET properties. Biomolecules such as quercetin, myricetin, and 7-hydroxyflavone were found to be promising inhibitors of p53 protein and may help in the protection of EBR induced DNA damage during cancer treatment.

Identifying the common genetic networks of ADR (adverse drug reaction) clusters and developing an ADR classification model

Adverse drug reactions (ADRs) are one of the major concerns threatening public health and have resulted in failures in drug development. Thus, predicting ADRs and discovering the mechanisms underlying ADRs have become important tasks in pharmacovigilance. Identification of potential ADRs by computational approaches in the early stages would be advantageous in drug development. Here we propose a computational method that elucidates the action mechanisms of ADRs and predicts potential ADRs by utilizing ADR genes, drug features, and protein-protein interaction (PPI) networks. If some ADRs share similar features, there is a high possibility that they may appear together in a drug and share analogous mechanisms. Proceeding from this assumption, we clustered ADRs according to interactions of ADR genes in the PPI networks and the frequency of co-occurrence of ADRs in drugs. ADR clusters were verified based on a side effect database and literature data regarding whether ADRs have relevance to other ADRs in the same cluster. Gene networks shared by ADRs in each cluster were constructed by cumulating the shortest paths between drug target genes and ADR genes in the PPI network. We developed a classification model to predict potential ADRs using these gene networks shared by ADRs and calculated cross-validation AUC (area under the curve) values for each ADR cluster. In addition, in order to demonstrate correlations between gene networks shared by ADRs and ADRs in a cluster, we applied the Wilcoxon rank sum statistical test to the literature data and results of a Google query search. We attained statistically meaningful p-values (<0.05) for every ADR cluster. The results suggest that our approach provides insights into discovering the action mechanisms of ADRs and is a novel attempt to predict ADRs in a biological aspect.

Drug-induced adverse events prediction with the LINCS L1000 data

MOTIVATION

Adverse drug reactions (ADRs) are a central consideration during drug development. Here we present a machine learning classifier to prioritize ADRs for approved drugs and pre-clinical small-molecule compounds by combining chemical structure (CS) and gene expression (GE) features. The GE data is from the Library of Integrated Network-based Cellular Signatures (LINCS) L1000 dataset that measured changes in GE before and after treatment of human cells with over 20 000 small-molecule compounds including most of the FDA-approved drugs. Using various benchmarking methods, we show that the integration of GE data with the CS of the drugs can significantly improve the predictability of ADRs. Moreover, transforming GE features to enrichment vectors of biological terms further improves the predictive capability of the classifiers. The most predictive biological-term features can assist in understanding the drug mechanisms of action. Finally, we applied the classifier to all  >20 000 small-molecules profiled, and developed a web portal for browsing and searching predictive small-molecule/ADR connections.

AVAILABILITY AND IMPLEMENTATION

The interface for the adverse event predictions for the  >20 000 LINCS compounds is available at http://maayanlab.net/SEP-L1000/ CONTACT: avi.maayan@mssm.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Targets of drugs are generally, and targets of drugs having side effects are specifically good spreaders of human interactome perturbations

Network-based methods are playing an increasingly important role in drug design. Our main question in this paper was whether the efficiency of drug target proteins to spread perturbations in the human interactome is larger if the binding drugs have side effects, as compared to those which have no reported side effects. Our results showed that in general, drug targets were better spreaders of perturbations than non-target proteins, and in particular, targets of drugs with side effects were also better spreaders of perturbations than targets of drugs having no reported side effects in human protein-protein interaction networks. Colorectal cancer-related proteins were good spreaders and had a high centrality, while type 2 diabetes-related proteins showed an average spreading efficiency and had an average centrality in the human interactome. Moreover, the interactome-distance between drug targets and disease-related proteins was higher in diabetes than in colorectal cancer. Our results may help a better understanding of the network position and dynamics of drug targets and disease-related proteins, and may contribute to develop additional, network-based tests to increase the potential safety of drug candidates.

Deep sequencing reveals important roles of microRNAs in response to drought and salinity stress in cotton

Drought and salinity are two major environmental factors adversely affecting plant growth and productivity. However, the regulatory mechanism is unknown. In this study, the potential roles of small regulatory microRNAs (miRNAs) in cotton response to those stresses were investigated. Using next-generation deep sequencing, a total of 337 miRNAs with precursors were identified, comprising 289 known miRNAs and 48 novel miRNAs. Of these miRNAs, 155 miRNAs were expressed differentially. Target prediction, Gene Ontology (GO)-based functional classification, and Kyoto Encyclopedia of Genes and Genomes (KEGG)-based functional enrichment show that these miRNAs might play roles in response to salinity and drought stresses through targeting a series of stress-related genes. Degradome sequencing analysis showed that at least 55 predicted target genes were further validated to be regulated by 60 miRNAs. CitationRank-based literature mining was employed to determinhe the importance of genes related to drought and salinity stress. The NAC, MYB, and MAPK families were ranked top under the context of drought and salinity, indicating their important roles for the plant to combat drought and salinity stress. According to target prediction, a series of cotton miRNAs are associated with these top-ranked genes, including miR164, miR172, miR396, miR1520, miR6158, ghr-n24, ghr-n56, and ghr-n59. Interestingly, 163 cotton miRNAs were also identified to target 210 genes that are important in fibre development. These results will contribute to cotton stress-resistant breeding as well as understanding fibre development.

Prediction of drug-target interactions for drug repositioning only based on genomic expression similarity

Small drug molecules usually bind to multiple protein targets or even unintended off-targets. Such drug promiscuity has often led to unwanted or unexplained drug reactions, resulting in side effects or drug repositioning opportunities. So it is always an important issue in pharmacology to identify potential drug-target interactions (DTI). However, DTI discovery by experiment remains a challenging task, due to high expense of time and resources. Many computational methods are therefore developed to predict DTI with high throughput biological and clinical data. Here, we initiatively demonstrate that the on-target and off-target effects could be characterized by drug-induced in vitro genomic expression changes, e.g. the data in Connectivity Map (CMap). Thus, unknown ligands of a certain target can be found from the compounds showing high gene-expression similarity to the known ligands. Then to clarify the detailed practice of CMap based DTI prediction, we objectively evaluate how well each target is characterized by CMap. The results suggest that (1) some targets are better characterized than others, so the prediction models specific to these well characterized targets would be more accurate and reliable; (2) in some cases, a family of ligands for the same target tend to interact with common off-targets, which may help increase the efficiency of DTI discovery and explain the mechanisms of complicated drug actions. In the present study, CMap expression similarity is proposed as a novel indicator of drug-target interactions. The detailed strategies of improving data quality by decreasing the batch effect and building prediction models are also effectively established. We believe the success in CMap can be further translated into other public and commercial data of genomic expression, thus increasing research productivity towards valid drug repositioning and minimal side effects.

Small drug molecules usually bind to unintended off-targets, leading to unexpected drug responses such as side effects or drug repositioning opportunities. Thus, identifying unintended drug-target interactions (DTI) is particularly required for understanding complicated drug actions. It remains expensive nowadays to experimentally determine DTI, so various computational methods are developed. In this study, we initiatively demonstrated that target binding is directly correlated with drug induced genomic expression profiles in Connectivity Map (CMap). By improving data quality of CMap, we illustrated three important facts: (1) Drugs binding to common targets show higher gene-expression similarity than random compounds, indicating that upstream ligand binding could be characterized by downstream gene-expression change. (2) It is found that some targets are better characterized by CMap than others. To guarantee efficiency of DTI discovery, prediction models should be specifically built for those well characterized targets. (3) It is broadly observed in the predicted DTI that ligands for the same target may collectively interact with common off-target. This observation is consistent with published experimental evidence and can help illustrate the mechanisms of unexplained drug reactions. Based on CMap, our work established an efficient pipeline of identifying potential DTI. By extending the success in CMap to other genomic data sources, we believe more DTI would be discovered.

Optimal scaling of digital transcriptomes

Deep sequencing of transcriptomes has become an indispensable tool for biology, enabling expression levels for thousands of genes to be compared across multiple samples. Since transcript counts scale with sequencing depth, counts from different samples must be normalized to a common scale prior to comparison. We analyzed fifteen existing and novel algorithms for normalizing transcript counts, and evaluated the effectiveness of the resulting normalizations. For this purpose we defined two novel and mutually independent metrics: (1) the number of “uniform” genes (genes whose normalized expression levels have a sufficiently low coefficient of variation), and (2) low Spearman correlation between normalized expression profiles of gene pairs. We also define four novel algorithms, one of which explicitly maximizes the number of uniform genes, and compared the performance of all fifteen algorithms. The two most commonly used methods (scaling to a fixed total value, or equalizing the expression of certain ‘housekeeping’ genes) yielded particularly poor results, surpassed even by normalization based on randomly selected gene sets. Conversely, seven of the algorithms approached what appears to be optimal normalization. Three of these algorithms rely on the identification of “ubiquitous” genes: genes expressed in all the samples studied, but never at very high or very low levels. We demonstrate that these include a “core” of genes expressed in many tissues in a mutually consistent pattern, which is suitable for use as an internal normalization guide. The new methods yield robustly normalized expression values, which is a prerequisite for the identification of differentially expressed and tissue-specific genes as potential biomarkers.

Structure and dynamics of molecular networks: a novel paradigm of drug discovery: a comprehensive review

Despite considerable progress in genome- and proteome-based high-throughput screening methods and in rational drug design, the increase in approved drugs in the past decade did not match the increase of drug development costs. Network description and analysis not only give a systems-level understanding of drug action and disease complexity, but can also help to improve the efficiency of drug design. We give a comprehensive assessment of the analytical tools of network topology and dynamics. The state-of-the-art use of chemical similarity, protein structure, protein-protein interaction, signaling, genetic interaction and metabolic networks in the discovery of drug targets is summarized. We propose that network targeting follows two basic strategies. The "central hit strategy" selectively targets central nodes/edges of the flexible networks of infectious agents or cancer cells to kill them. The "network influence strategy" works against other diseases, where an efficient reconfiguration of rigid networks needs to be achieved by targeting the neighbors of central nodes/edges. It is shown how network techniques can help in the identification of single-target, edgetic, multi-target and allo-network drug target candidates. We review the recent boom in network methods helping hit identification, lead selection optimizing drug efficacy, as well as minimizing side-effects and drug toxicity. Successful network-based drug development strategies are shown through the examples of infections, cancer, metabolic diseases, neurodegenerative diseases and aging. Summarizing >1200 references we suggest an optimized protocol of network-aided drug development, and provide a list of systems-level hallmarks of drug quality. Finally, we highlight network-related drug development trends helping to achieve these hallmarks by a cohesive, global approach.

Compensating for literature annotation bias when predicting novel drug-disease relationships through Medical Subject Heading Over-representation Profile (MeSHOP) similarity

Background

Using annotations to the articles in MEDLINE^®/PubMed^®, over six thousand chemical compounds with pharmacological actions have been tracked since 1996. Medical Subject Heading Over-representation Profiles (MeSHOPs) quantitatively leverage the literature associated with biological entities such as diseases or drugs, providing the opportunity to reposition known compounds towards novel disease applications.

Methods

A MeSHOP is constructed by counting the number of times each medical subject term is assigned to an entity-related research publication in the MEDLINE database and calculating the significance of the count by comparing against the count of the term in a background set of publications. Based on the expectation that drugs suitable for treatment of a disease (or disease symptom) will have similar annotation properties to the disease, we successfully predict drug-disease associations by comparing MeSHOPs of diseases and drugs.

Results

The MeSHOP comparison approach delivers an 11% improvement over bibliometric baselines. However, novel drug-disease associations are observed to be biased towards drugs and diseases with more publications. To account for the annotation biases, a correction procedure is introduced and evaluated.

Conclusions

By explicitly accounting for the annotation bias, unexpectedly similar drug-disease pairs are highlighted as candidates for drug repositioning research. MeSHOPs are shown to provide a literature-supported perspective for discovery of new links between drugs and diseases based on pre-existing knowledge.

Systematic drug repositioning based on clinical side-effects

Drug repositioning helps fully explore indications for marketed drugs and clinical candidates. Here we show that the clinical side-effects (SEs) provide a human phenotypic profile for the drug, and this profile can suggest additional disease indications. We extracted 3,175 SE-disease relationships by combining the SE-drug relationships from drug labels and the drug-disease relationships from PharmGKB. Many relationships provide explicit repositioning hypotheses, such as drugs causing hypoglycemia are potential candidates for diabetes. We built Naïve Bayes models to predict indications for 145 diseases using the SEs as features. The AUC was above 0.8 in 92% of these models. The method was extended to predict indications for clinical compounds, 36% of the models achieved AUC above 0.7. This suggests that closer attention should be paid to the SEs observed in trials not just to evaluate the harmful effects, but also to rationally explore the repositioning potential based on this “clinical phenotypic assay”.

Identification of hub genes related to the recovery phase of irradiation injury by microarray and integrated gene network analysis

BACKGROUND

Irradiation commonly causes long-term bone marrow injury charactertized by defective HSC self-renewal and a decrease in HSC reserve. However, the effect of high-dose IR on global gene expression during bone marrow recovery remains unknown.

METHODOLOGY

Microarray analysis was used to identify differentially expressed genes that are likely to be critical for bone marrow recovery. Multiple bioinformatics analyses were conducted to identify key hub genes, pathways and biological processes.

PRINCIPAL FINDINGS

1) We identified 1302 differentially expressed genes in murine bone marrow at 3, 7, 11 and 21 days after irradiation. Eleven of these genes are known to be HSC self-renewal associated genes, including Adipoq, Ccl3, Ccnd1, Ccnd2, Cdkn1a, Cxcl12, Junb, Pten, Tal1, Thy1 and Tnf; 2) These 1302 differentially expressed genes function in multiple biological processes of immunity, including hematopoiesis and response to stimuli, and cellular processes including cell proliferation, differentiation, adhesion and signaling; 3) Dynamic Gene Network analysis identified a subgroup of 25 core genes that participate in immune response, regulation of transcription and nucleosome assembly; 4) A comparison of our data with known irradiation-related genes extracted from literature showed 42 genes that matched the results of our microarray analysis, thus demonstrated consistency between studies; 5) Protein-protein interaction network and pathway analyses indicated several essential protein-protein interactions and signaling pathways, including focal adhesion and several immune-related signaling pathways.

CONCLUSIONS

Comparisons to other gene array datasets indicate that global gene expression profiles of irradiation damaged bone marrow show significant differences between injury and recovery phases. Our data suggest that immune response (including hematopoiesis) can be considered as a critical biological process in bone marrow recovery. Several critical hub genes that are key members of significant pathways or gene networks were identified by our comprehensive analysis.

An algorithmic framework for predicting side effects of drugs

One of the critical stages in drug development is the identification of potential side effects for promising drug leads. Large-scale clinical experiments aimed at discovering such side effects are very costly and may miss subtle or rare side effects. Previous attempts to systematically predict side effects are sparse and consider each side effect independently. In this work, we report on a novel approach to predict the side effects of a given drug, taking into consideration information on other drugs and their side effects. Starting from a query drug, a combination of canonical correlation analysis and network-based diffusion is applied to predict its side effects. We evaluate our method by measuring its performance in a cross validation setting using a comprehensive data set of 692 drugs and their known side effects derived from package inserts. For 34% of the drugs, the top scoring side effect matches a known side effect of the drug. Remarkably, even on unseen data, our method is able to infer side effects that highly match existing knowledge. In addition, we show that our method outperforms a prediction scheme that considers each side effect separately. Our method thus represents a promising step toward shortcutting the process and reducing the cost of side effect elucidation.

DRAR-CPI: a server for identifying drug repositioning potential and adverse drug reactions via the chemical-protein interactome

Identifying new indications for existing drugs (drug repositioning) is an efficient way of maximizing their potential. Adverse drug reaction (ADR) is one of the leading causes of death among hospitalized patients. As both new indications and ADRs are caused by unexpected chemical-protein interactions on off-targets, it is reasonable to predict these interactions by mining the chemical-protein interactome (CPI). Making such predictions has recently been facilitated by a web server named DRAR-CPI. This server has a representative collection of drug molecules and targetable human proteins built up from our work in drug repositioning and ADR. When a user submits a molecule, the server will give the positive or negative association scores between the user's molecule and our library drugs based on their interaction profiles towards the targets. Users can thus predict the indications or ADRs of their molecule based on the association scores towards our library drugs. We have matched our predictions of drug-drug associations with those predicted via gene-expression profiles, achieving a matching rate as high as 74%. We have also successfully predicted the connections between anti-psychotics and anti-infectives, indicating the underlying relevance of anti-psychotics in the potential treatment of infections, vice versa. This server is freely available at http://cpi.bio-x.cn/drar/.

Exploring off-targets and off-systems for adverse drug reactions via chemical-protein interactome--clozapine-induced agranulocytosis as a case study

In the era of personalized medical practice, understanding the genetic basis of patient-specific adverse drug reaction (ADR) is a major challenge. Clozapine provides effective treatments for schizophrenia but its usage is limited because of life-threatening agranulocytosis. A recent high impact study showed the necessity of moving clozapine to a first line drug, thus identifying the biomarkers for drug-induced agranulocytosis has become important. Here we report a methodology termed as antithesis chemical-protein interactome (CPI), which utilizes the docking method to mimic the differences in the drug-protein interactions across a panel of human proteins. Using this method, we identified HSPA1A, a known susceptibility gene for CIA, to be the off-target of clozapine. Furthermore, the mRNA expression of HSPA1A-related genes (off-target associated systems) was also found to be differentially expressed in clozapine treated leukemia cell line. Apart from identifying the CIA causal genes we identified several novel candidate genes which could be responsible for agranulocytosis. Proteins related to reactive oxygen clearance system, such as oxidoreductases and glutathione metabolite enzymes, were significantly enriched in the antithesis CPI. This methodology conducted a multi-dimensional analysis of drugs' perturbation to the biological system, investigating both the off-targets and the associated off-systems to explore the molecular basis of an adverse event or the new uses for old drugs.

Idiosyncratic drug reactions (IDR) generally cannot be identified until after a drug is taken by a large population, but usually result in restricted use or withdrawal. Clozapine provides the most effective treatment for schizophrenia but its use is limited because of a life-threatening IDR, i.e., the agranulocytosis. A high impact clinical study demonstrated the necessity of moving clozapine from 3^rd line to 1^st line drug; therefore, intensive research has aimed at identifying genes responsible for clozapine-induced agranulocytosis (CIA). Olanzapine, an analog of clozapine, has much lower incidence of agranulocytosis. Based on this phenomenon, we proposed an in silico methodology termed as antithesis chemical-protein interactome (CPI), which mimics the differences in the drug-protein interactions of the two drugs across a panel of human proteins. e.g., HSPA1A was identified to be targeted by clozapine not olanzapine. Furthermore, the gene expression of the HSPA1A-related gene system was also found up-regulated after clozapine treatment. This approach can examine the system's perturbation in terms of both the off-target and the off-system's interaction with the drug, providing theoretical basis for decoding the adverse drug reactions or the new uses for old drugs.

The role of translational bioinformatics in drug discovery

The application of translational approaches (e.g. from bed to bench and back) is gaining momentum in the pharmaceutical industry. By utilizing the rapidly increasing volume of data at all phases of drug discovery, translational bioinformatics is poised to address some of the key challenges faced by the industry. Indeed, computational analysis of clinical data and patient records has informed decision-making in multiple aspects of drug discovery and development. Here, we review key examples of translational bioinformatics approaches to emphasize its potential to enhance the quality of drug discovery pipelines, reduce attrition rates and, ultimately, lead to more effective treatments.

It's the machine that matters: Predicting gene function and phenotype from protein networks

Increasing knowledge about the organization of proteins into complexes, systems, and pathways has led to a flowering of theoretical approaches for exploiting this knowledge in order to better learn the functions of proteins and their roles underlying phenotypic traits and diseases. Much of this body of theory has been developed and tested in model organisms, relying on their relative simplicity and genetic and biochemical tractability to accelerate the research. In this review, we discuss several of the major approaches for computationally integrating proteomics and genomics observations into integrated protein networks, then applying guilt-by-association in these networks in order to identify genes underlying traits. Recent trends in this field include a rising appreciation of the modular network organization of proteins underlying traits or mutational phenotypes, and how to exploit such protein modularity using computational approaches related to the internet search algorithm PageRank. Many protein network-based predictions have recently been experimentally confirmed in yeast, worms, plants, and mice, and several successful approaches in model organisms have been directly translated to analyze human disease, with notable recent applications to glioma and breast cancer prognosis.

Role of systems pharmacology in understanding drug adverse events

Systems pharmacology involves the application of systems biology approaches, combining large-scale experimental studies with computational analyses, to the study of drugs, drug targets, and drug effects. Many of these initial studies have focused on identifying new drug targets, new uses of known drugs, and systems-level properties of existing drugs. This review focuses on systems pharmacology studies that aim to better understand drug side effects and adverse events. By studying the drugs in the context of cellular networks, these studies provide insights into adverse events caused by off-targets of drugs as well as adverse events-mediated complex network responses. This allows rapid identification of biomarkers for side effect susceptibility. In this way, systems pharmacology will lead to not only newer and more effective therapies, but safer medications with fewer side effects.

ReCGiP, a database of reproduction candidate genes in pigs based on bibliomics

BACKGROUND

Reproduction in pigs is one of the most economically important traits. To improve the reproductive performances, numerous studies have focused on the identification of candidate genes. However, it is hard for one to read all literatures thoroughly to get information. So we have developed a database providing candidate genes for reproductive researches in pig by mining and processing existing biological literatures in human and pigs, named as ReCGiP.

DESCRIPTION

Based on text-mining and comparative genomics, ReCGiP presents diverse information of reproduction-relevant genes in human and pig. The genes were sorted by the degree of relevance with the reproduction topics and were visualized in a gene's co-occurrence network where two genes were connected if they were co-cited in a PubMed abstract. The 'hub' genes which had more 'neighbors' were thought to be have more important functions and could be identified by the user in their web browser. In addition, ReCGiP provided integrated GO annotation, OMIM and biological pathway information collected from the Internet. Both pig and human gene information can be found in the database, which is now available.

CONCLUSIONS

ReCGiP is a unique database providing information on reproduction related genes for pig. It can be used in the area of the molecular genetics, the genetic linkage map, and the breeding of the pig and other livestock. Moreover, it can be used as a reference for human reproduction research.