BACKGROUND
bioinformatics tools currently allow research to focus on the integration of large-scale data generated by sequencing, differential expression analysis, gene interaction studies and others.
several initiatives exist to organize this knowledge in secondary databases, thus allowing easier access and visualization.
databases containing interaction information are a good source for novel research.
ihop  <cit>  allows users to tag gene names of interest and browse through the related pubmed literature with highlighted keywords.
another interaction database is string  <cit> , which contains physical interactions and functional associations between proteins and integrates data retrieved from literature , genomic context, large scale experiments and conserved co-expression.
text-mining, therefore, has a fundamental role in these tools and allows access to interactions spread throughout the literature.
the extraction of biological events from literature through text-mining tools is essential to not only update the interaction databases but also for the creation and annotation of pathways.
metabolic and regulatory pathways are an example of organized knowledge that allow a better visualization of a complex system and can be found in databases such as ipath  <cit> , biocyc  <cit>  or kegg pathways  <cit> .
when orthology information is added to pathways, the same process can be represented in different organisms.
orthology is also an important tool for sequence annotation.
current orthologue databases such as cog and kog  <cit> , eggnog  <cit> , orthomcl  <cit>  and kegg orthology  <cit>  all provide a good source for manually curated clusters of orthologues defined for organisms with complete genomes.
we developed a procedure to enrich the cog database with uniref <dig> clusters from the uniprot database  <cit> , creating the uecog database  <cit> .
recently, a similar procedure was applied to the kegg orthology database creating the enriched ueko database .
the available tools described raise the possibility of integrating current information and generating complex regulatory pathways.
previous publications individually reported the regulatory interactions that control preimplantation embryo development  <cit> .
however, a complete preimplantation development regulatory pathway has never been built.
in humans, the preimplantation phase of embryonic development is a period of approximately six days after fertilization prior to attachment of the embryo to the uterine wall.
implantation can occur before or in the seventh embryonic day , a time during which the uterus is receptive  <cit> .
mammalian embryonic development has been thoroughly studied in mice and the blastomeres remain totipotent, able to generate any other cell, up to the eight-cell stage, unlike other animals  <cit> .
after fertilization, successive cleavages take place during the first two days of development, resulting in the eight-cell embryo.
the next stage of development is called the morula stage.
an increase in cell-cell contact results in formation of a compacted morula.
the subsequent divisions increase the complexity of the embryo and cells may be located on the inside, surrounded by other cells, or on the outside, in contact with the environment.
the identification of the initial cells for each lineage has shown that the trophectoderm  is derived mostly from the outer cells, whereas the inner cells give rise to the inner cell mass .
later, the icm divides into the primitive endoderm  and the epiblast .
during the differentiation of the te from the icm, the blastocoel is formed through a process of cavitation.
the embryo is called a blastocyst when all three structures are present .
twenty-four hours after blastocyst formation occurs, the last stage of preimplantation development takes place when the pe differentiates from the icm.
the three lineages thus formed in preimplantation development present different fates during subsequent embryonic development.
while the epiblast, which forms from the icm following implantation, is still undifferentiated and will give rise to the fetus itself, the trophectoderm will become the fetal portion of the placenta and the primitive endoderm  will form the yolk sac  <cit> .
complex regulatory processes such as animal development are a result of the interaction of many different gene products and elements that control the expression of these genes.
traditional experiments that determine the function of one or a few genes are essential, but do not result in a comprehensive view of complex systems.
a complex regulatory network should be able to portray specific and general aspects of development, such as the embryonic fate of certain cells  <cit> .
in this work, we noticed the absence in databases of a pathway describing the preimplantation phase of embryo development and sought to develop the given pathway using text-mining tools, complementing it with orthology information.
the resulting pathway comprises  <dig> genes and the interactions between them.
clusters of orthologous groups were generated for each gene represented in the pathway and provided the necessary information to determine the last common ancestor.
this determination revealed that the preimplantation development pathway is an ancient chordata pathway with addition of modern elements throughout evolution.
RESULTS
text-mining
initially, we used the pubmed platform to search for articles related to the embryo preimplantation development  and obtained  <dig> entries as a result.
to obtain a more efficient set of articles with relevancy to our work, the result entries were submitted to medlineranker  <cit> .
this software computes discriminating words by comparing a set of user selected abstracts indicated as highly relevant to a background set and then scores any abstracts in terms of their content of those discriminating words.
after the classification, we selected the top  <dig> abstracts for further analysis, which presented p-value lower than  <dig>  and by manual inspection provide large amount of information when uploaded in pescador.
since human and mouse embryo development are highly similar, it was plausible to use abstracts from work on both organisms as source of information for the preimplantation pathway construction, paying attention to any possible conflict.
using these  <dig> highly informative abstracts as our input, pescador  an online platform for friendly operation of the laitor software  <cit> ) was used for tagging of gene names and biointeractions extraction from each abstract.
as a result,  <dig> gene names were tagged and  <dig> type  <dig> biointeractions were highlighted as well as other informative biointeractions.
biointeractions are classified by laitor  <cit>  as type  <dig> when in the same sentence the software encounters a gene name, a biointeraction word and another gene name, in that order .
from these tagged abstracts we manually curated the information and constructed the pathway for the preimplantation embryo development describing  <dig> genes and numerous interactions between them during the early developmental stages, trophectoderm differentiation from the inner cell mass and posterior extraembrionary endoderm differentiation.
a sample abstract tagged by pescador and the manual extraction of the information it contains is exemplified in figure  <dig>  the pathway shown in figure  <dig> was constructed according to kgml .
the large decrease in the initial number of genes tagged in the abstracts is mainly due to redundancy between abstracts  and also to genes tagged in type  <dig> and  <dig> biointeractions, which not always result in pathway building information.
preimplantation pathway
the pathway obtained after the analysis of all the abstracts from the pescador output is represented in figure  <dig> and the regulations are reviewed below.
first embryonic cleavages
the oncogene c-myc is an important transcriptional regulator and its expression is observed in the initial stages of development, where it is present in embryonic cells until the morula stage and repressed thereafter  <cit> .
two additional genes recently associated with these early developmental stages are boris and ecsa.
boris is involved in early development following fertilization and soon afterwards repressed, and ecsa, expression begins in the blastocyst exclusively in the cells of the inner cell mass .
the presence of these genes was compared to the expression pattern of the oct <dig> transcription factor, which is present in the early cleavages, repressed after this initial stage, and then its expression is afterwards stimulated again in the blastocyst  <cit> .
the expression of the gametogenesis associated gene gse was also recently identified in cells of the early embryo; later this protein is found only in the icm, suggesting a role in the specification of cell lineage  <cit> .
methylation patterns and correct preimplantation development
genomic methylation patterns in mammalian cells depend on dnmt <dig> .
in the mouse, an embryo-specific variant called dnmt1o is expressed in the early stages of development.
in the 8-cell stage this protein relocates to the cell nucleus where maintains essential methylation patterns, allowing embryos to complete early developmental events  <cit> .
it was recently shown  <cit>  that the inability of dnmt1o to properly relocate not only results in a developmental arrest at the 5- <dig> cell stage, but is also responsible for the downregulation of five genes involved in the formation of gap and tight junctions .
these junctions are crucial for early processes such as compaction of the 8-cell embryo and cavitation of the blastocoel.
te versus icm dichotomy: key role of lats controlling tead <dig> co-activator yap
cells destined to become part of the icm are marked by repression of two genes   <cit>  and by upregulation of sox <dig>  <cit> .
in these cells, the major pluripotency transcription factors, including nanog and oct <dig>  remain active due to the expression of an important player and member of the hippo signaling pathway: lats.
this serine/threonine protein kinase is responsible for phosphorylating yap, leading to its cytoplasmic localization and thus preventing its association with the transcription factor tead <dig>
triggering te differentiation: tead4/yap target cdx <dig> to repress nanog and oct4
conversely, in the outer cells that will differentiate and form the trophectoderm, yap is unphosphorylated, remains in the nucleus and associates with tead <dig>  leading to the activation of cdx <dig>  a key repressor of nanog and oct <dig>  <cit> .
repression of oct <dig> and nanog transcription by cdx <dig> then releases the inhibition that these two key factors were exerting on many different genes, in turn activating these targets  <cit> .
activation of cdx <dig> requires release from basal repression; nanog  <cit>  and oct <dig>  <cit>  repress basal levels of cdx <dig> and induction of higher levels of cdx <dig> by tead4/yap overcomes this repression, allowing cdx <dig> to play its role  <cit> .
tead <dig> was also recently determined to activate another trophectoderm differentiation factor, gata <dig>  <cit> , which acts alongside cdx <dig> and affects transcription of a number of genes independent of cdx <dig>  the tead4-dependent activation of gata <dig> seems to be independent of yap, suggesting tead <dig> interacts with another partner as well as yap.
also required for high level expression of cdx <dig> in trophectoderm cells is the cell motility protein arp3; experiments with complete knockdown of this protein show trophoblast cells unable to develop properly, possibly undergoing apoptosis as a result of loss of cdx <dig>  <cit> .
the tgfbeta pathway is another important pathway for trophectoderm differentiation; tgfbeta signaling is stimulated by bmp <dig>  which leads to the activation of smad proteins.
these proteins can also stimulate transcription of cdx <dig>  <cit> , and bmp <dig> is known to inhibit id <dig>  an inhibitor of differentiation  <cit> , and to activate hand <dig>  which is involved in trophoblast cell differentiation  <cit> .
in the absence of oct <dig> and nanog
the downregulation of oct <dig> in the outer cells of the embryo leads to the activation of a positive regulator of te cell fate, eomes   <cit> , which is also a possible cdx <dig> target  <cit> .
the subsequent differentiation of these cells into trophectoderm is accompanied by the expression of several genes, such as the glycoprotein psg <dig>  <cit>  and the marker krt <dig>  psg <dig> and krt <dig> expression are among the first signs that a blastomere has lost its totipotent competence, prior to any visible differentiation  <cit> .
removal of oct4-dependent repression also results in activation of genes such as etif2b and rps <dig>  <cit> , allowing these cells to engage in an intense translation routine.
knockdown studies targeting oct <dig> also show that it represses the expression of gcm <dig>  which is normally placenta specific  <cit> , and of the hcg hormone’s beta chain  <cit> .
concurrently, nanog downregulation allows the expression of a number of genes associated with both trophectoderm  and extraembryonic endoderm   <cit> .
these latter genes will in turn initiate the formation of tissues such as the primitive endoderm, a component of the yolk sac.
from the early blastocyst stage on, desmosomes are assembled in the trophectoderm in response to desmocollin , which is also not expressed in the icm  <cit> .
thus, tead4/yap activation of cdx <dig>  accompanied by the subsequent repression of nanog and oct <dig>  describes a scenario for the te differentiation.
underneath the maintained activation of oct <dig> and nanog
back in the icm, the main pluripotency genes remain active and form a complex regulation pathway.
recently it was discovered that transcription of nanog is further stimulated by the presence of compounds such as retinol  <cit> .
klf <dig>  klf <dig> and klf <dig> exert a redundant role in the activation of nanog.
these krüppel-like factors were described as essential for the maintenance of pluripotency.
indeed, klf <dig> was already known for this role and is commonly used in reprogramming of differentiated cells into induced pluripotent stem cells.
however, only the simultaneous depletion of klf <dig>   <dig> and  <dig> results in the differentiation of stem cells, indicating functional redundancy  <cit> .
other proteins known to activate nanog include the two other main pluripotency regulators, oct <dig>  <cit>  and sox <dig>  <cit> .
the estrogen receptor esrrb is also reported to be involved in the activation of nanog by oct <dig> and sox <dig>  <cit> .
conversely, nanog can activate oct <dig>  <cit> , and esrrb is necessary to maintain oct <dig> promoter activity  <cit> .
each of the three key factors, oct <dig>  sox <dig> and nanog, also act as self-activators, e.g.
the partners oct <dig> and sox <dig> bind and activate oct <dig> transcription  <cit> .
another key transcription factor involved in the maintenance of cell pluripotency is sall <dig>  <cit> .
sall <dig> binds to the conserved regulatory region in the pou5f <dig>  distal enhancer and activates its transcription  <cit> .
studies with microrna interference of sall <dig> show that the loss of this factor leads to reduction of oct <dig> mrna levels and significant expression of cdx <dig> in the icm  <cit> .
b-myb, a gene expressed in proliferating cells, is also a positive regulator of oct <dig> and studies report early differentiation of icm in the absence of b-myb  <cit> .
the notch signaling pathway is a conserved pathway that is involved in cellular communication processes and correct cell fate decisions that also has a role in icm development  <cit> .
nle protein, a direct regulator of this pathway, is essential for survival of the icm  <cit> .
another protein associated with development and survival of the icm is tbn , whose absence promotes cell apoptosis in the icm  <cit> .
expression of the platelet and endothelial cell adhesion molecule  was detected by immunofluorescence confocal microscopy in the blastocyst and restricted to the icm cells.
subsequently, pecam <dig> remains only in the pluripotent epiblast cells, disappearing the moment these cells undergo differentiation  <cit> , and indicating a new role for this molecule during embryo development.
activation, but with moderation
other control pathways maintain expression of these genes at a steady-state concentration and balance these many mechanisms for activation and upregulation of transcription.
a complex regulation feedback loop consists of foxd <dig>  nanog and oct <dig>  <cit> .
to keep oct <dig> and nanog expression within steady-state levels, these three genes interact so that  expression of nanog activates foxd <dig> and oct <dig> but not above steady-state levels due to oct <dig> exerted repression; and  foxd <dig> and nanog activate oct <dig> expression but not above steady-state levels due to oct <dig> self-repression.
dax <dig> is an orphan nuclear hormone receptor recently identified as a repressor of oct <dig> transcription  <cit> .
dax <dig> expression was also capable of reducing nanog and rex <dig> expression.
assays show that dax <dig> binds to oct <dig> and abolishes its dna binding activity, thus decreasing the transcription of nanog and rex <dig>  targets of oct <dig> activation.
another repressor in the icm is tcf <dig>  a wnt signaling pathway effector.
tle <dig>  and ctbp  are key partners of tcf <dig> in mediating this repressive effect.
tcf <dig> binds to and represses the oct <dig> promoter, and this repressive effect requires both the groucho and ctbp interacting domains of tcf <dig>  <cit> .
tcf <dig> also limits the steady-state levels of nanog mrna, protein, and promoter activity in self-renewing embryonic stem cells ; the tcf <dig> groucho domain is involved in this repression  <cit> .
thus, tcf <dig> is critical for maintaining the appropriate levels of both oct <dig> and nanog in escs.
experiments show that loss of tcf <dig> by rna interference  knockdown blocks the ability of escs to differentiate  <cit> , emphasizing the importance of this interaction.
downstream of oct <dig> and sox2
oct <dig> activates embryonic stem cell-specific gene  <dig> , which encodes an rna binding protein present in the icm that is responsible for regulating several specific target transcripts  <cit> .
oct <dig> and sox <dig> are also responsible for the regulation of the fibroblast growth factor  <dig>   <cit> .
expression of fgf <dig>  therefore, requires the combined activity of these two transcription factors that bind to adjacent sites on the fgf <dig> enhancer dna region  <cit> .
once expressed, the fgf <dig> protein can interact with its receptor fgfr <dig> and activate icm and adjacent te cell proliferation, activating extraembryonic endoderm cells as well in later stages.
several other genes with important functions in embryonic development are also targets of oct4-dependent activation.
these include growth factor tdgf <dig>  growth inhibitor sap <dig>  regulator of nonsense transcripts rent <dig>  two proteins involved in stem cell self-renewal dppa <dig> and dppa <dig> , anterior visceral endoderm  markers lefty <dig> and lefty <dig>  surface antigen thy <dig>  and other genes encoding proteins involved in specialized cellular processes  and hypothetical proteins   <cit> .
the master regulation exerted by sox <dig> and oct <dig> during mammalian embryogenesis is believed to operate through their cooperative binding to dna regulatory regions composed of adjacent hmg and pou motifs   <cit> .
exemplifying this arrangement, dppa <dig> is one such gene with the presence of an hmg/pou cassette in its promoter region  <cit> .
downstream of nanog and stat3
activation of jak/stat pathway also has an important contribution to pluripotency.
in mice, the lif/stat <dig> pathway  <cit>  for maintenance of cell pluripotency comprises lif and lif receptor, which deliver intracellular signaling through stat <dig>  stat <dig>  a signal transducer and activator of transcription is activated by the jak <dig> kinase and binds to several promoters inducing transcription of pluripotency related genes  <cit> .
nanog and stat <dig> were found to bind to and synergistically activate stat3-dependent promoters  <cit> .
nanog also functions as a transcriptional inhibitor to nfκb, a factor known to have pro-differentiation activity  <cit> .
nanog is also responsible for smad <dig> repression, thereby preventing bmp4-induced differentiation through the tgfbeta signaling pathway, for which smad <dig> is a key signal transducer  <cit> .
extraembryonic endoderm differentiation from icm cells
prior to embryo implantation one more differentiation takes place.
certain cells from the icm give rise to the primitive endoderm, the first morphologically distinct cell type of the extraembryonic endoderm.
the extraembryonic endoderm comprises the primitive, parietal and visceral endoderm components and will become the yolk sac during posterior development stages.
wnt <dig> was recently identified as an inducer of primitive endoderm and this induction is accompanied by translocation of beta-catenin  and snail <dig> to the nucleus  <cit> .
this study also showed that up-regulation of protein kinase a  induces markers of parietal endoderm.
another wnt family member, wnt9a, is expressed only in icm cells that surround the bastocoel  <cit>  and induces repositioning of the cells expressing gata <dig>  which is necessary for formation of primitive endoderm  <cit> .
sox <dig> plays a major role in parietal endoderm differentiation.
through studies with short interfering rna molecules, it was established that sox <dig> is responsible for transcription induction of gata <dig> and gata <dig>  <cit> .
individual or combined silencing of sox <dig>  gata <dig> and gata <dig> result in suppression of cell shape changes and production of laminin- <dig> , characteristic changes present in parietal endoderm differentiation  <cit> .
gata <dig> was previously identified as a transcription factor responsible for the activation of fgf <dig>  <cit> .
sox <dig> also activates the fgf <dig> promoter.
conversely, sox <dig> can negatively modulate the gata4-dependent activation of fgf <dig>  which is supported by the role of this factor in icm pluripotency  <cit> .
another sox family member, sox <dig>  is responsible for the differentiation of the extraembryonic endoderm in the final steps of preimplantation development  <cit> .
the runx <dig> factor is associated with the expression of sox <dig> and is also specific for the extraembryonic endoderm  <cit> .
hnf <dig> is a transcription factor specific of the extraembryonic endoderm with subsequent roles in post-implantation development and organogenesis  <cit> .
its expression may result from bmp4-induced differentiation  <cit> .
finally, the dab <dig> protein is indispensable for the development of visceral endoderm; though its exact role is still not established, it is perhaps related to correct cell positioning  <cit> .
the expression of cer <dig>  a marker of the anterior visceral endoderm , commences before embryo implantation in the subset of cells that comprise the primitive endoderm.
this ancestral population includes both cells expressing cer <dig> together with cells in which cer <dig> expression begins after implantation and formation of the ave  <cit> .
search for homologues
to establish an ortholog database and provide sequence information to the genes contained in the preimplantation pathway, aminoacid sequences corresponding to the human and mouse gene products were used as seed for the software seedserver .
in fact, only the uniprot identifier for these proteins is necessary to execute seedserver - gene symbols were verified in the ncbi gene database and converted to the corresponding geneid, and the desired identifiers were obtained afterwards from the uniprot database.
for each gene a cluster of homologues was generated comprising from  <dig> to  <dig> sequences .
the recruited sequences contained in each cluster can be swiss-prot annotated or unrevised trembl sequences.
in total, 25% of the cluster sequences are swiss-prot, the great majority of clusters being comprised of trembl sequences .
the search for homologues through seedserver provides therefore a large amount of candidates for manual curation in swiss-prot.
furthermore, seedserver can recruit sequences from organisms without a complete genome due to its use of ueko  and bidirectional best hit  searches conducted by seedlinkage  <cit> , and in fact only 27% of the sequences present in all clusters are from organisms with a complete genome.
the ortholog clustering by seedserver was only performed for genes that had a corresponding swissprot annotated gene product to be used as seed, therefore hrscp and dppa <dig>  which are described in the pathway, did not go through this analysis.
pathway ancestry
we then focused on the putative origin of these genes, determining which clade in the human lineage  shares each gene.
the generation of ortholog clusters allowed for the determination of the last common ancestor  for each of the genes in the pathway.
figure  <dig> shows the genes according to their origin.
genes were arbitrarily considered ancient for this analysis if their last common ancestor originated before the divergence of the clade euteleostomi and are coloured grey.
genes with a lca belonging to the clade euteleostomi or originated after divergence of euteleostomi are considered recent genes and are coloured blue.
ancient origin genes with an ortholog in drosophila melanogaster are marked with a red asterisk.
this arbitrary classification was meant to attract attention to the two key pluripotency controlling genes, nanog  and oct <dig> .
the graph shown in figure  <dig> represents the distribution of all the genes in the pathway according to their origin respect to clades of the human lineage.
it may be observed that a large quantity of genes originates in certain periods as seen in eumetazoa, coelomata, euteleostomi and eutheria.
the reasons for this wavelike origin need to be further analysed.
on the other hand, the apparent origin of complex structures, that characterize all descendents from a certain moment of evolution, might have occurred simultaneously to the specialization of gene groups.
the coverage of genomic sequences in the database is far from homogeneous and can influence the shape of this graph  <cit> .
in any case, the pattern observed agrees with the expansion of protein families related to stem cell markers observed in the ray-finned fish, that is, after divergence of the euteleostomi  <cit> .
furthermore, we searched for functional information related to the d. melanogaster orthologues in order to determine if these functions are somehow similar or related to the functions of the corresponding pathway genes.
this was done through a second text mining approach similar to the first and from the information recovered a secondary pathway was generated simply to illustrate the ortholog genes and their relative functional roles .
the regulatory pathways in which these genes are involved show us that these genes are all related to some part of drosophila embryo development, some of them with highly conserved functions still observed in the preimplantation pathway described.
an example is the hippo signalling pathway, which is extremely conserved, showing wts  phosphorylating yki ; this modification prevents yki interaction with sd .
the correlation between the human gene names and corresponding d. melanogaster ortholog names can be found in additional file  <dig> and also the pmid reference for the gene function in drosophila development.
discussion
the use of text-mining tools for the generation of regulatory pathways is an effective approach and it is important for the current interest of gathering data related to an organism or biological process.
the search for information related to a specific concept such as “preimplantation development” resulted in the selection of data related to this process only.
when other tools such as ihop  <cit>  and string  <cit>  are used for the search of biointeractions, it is necessary to know the names for the genes you are interested on and the information is then retrieved.
moreover in the case of ihop, the information retrieved consists of a large list of papers related to the gene of interest, which need to be manually analysed to extract the information related to the specific process.
in the case of string, the result of a query is a network of direct associations to other genes, which can be activations, repressions, or unknown, but for which it is not possible to perform a search restricting the query to a specific process for which you seek to determine the involvement of a given gene.
the approach described in this work  summarized in figure  <dig>  allows the researcher to initiate the study of a pathway without knowing exactly the genes involved, simply by selecting the published information related to the process of interest.
the manual curation required to create a pathway through this approach is significantly smaller.
however, the verification of all the interactions highlighted by the tool is essential.
text-mining is not able to eliminate the selection of false interaction pairs; in the case of laitor , the type  <dig> and  <dig> interactions can present genes with no association specified in the text  <cit> .
the text-mining data contribute the complete description of the pathway in the form of a literature review, a necessary step for the validation of the regulations represented, and for the inclusion of the pathway in a specific database, such as kegg pathway.
the establishment of this procedure for pathway generation allows future work to enlarge the knowledge on subjects still not approached, such as regulatory pathways for several types of cancer, mechanisms of pathogen resistance in plants and response to abiotic stresses in plants, among other themes of interest.
the inclusion of the preimplantation pathway in databases such as the kegg database will allow automatic annotation for several other organisms, as it is usually done in this database.
concurrently, a laboratory with a specific interest can promptly build a similar pathway for its local use.
from the  <dig> genes present in the pathway,  <dig> do not possess entries in kegg orthology and would constitute important additions.
considering that the contribution of kegg for the sequence recruitment in the seedserver clusters is only 25% of the total number of sequences, some organisms evolutionarily divergent from the ones represented in kegg begin to play a more relevant role for a more efficient annotation of new sequences.
it is relevant to stress that only the seedlinkage and ueko components of seedserver are capable of clustering sequences proceeding from organisms without a complete genome project.
moreover, linkage of recruited to seed sequences are verified with psi-blast.
another important contribution from the ortholog clustering by seedserver is the identification of candidates for swiss-prot annotation.
swiss-prot annotation depends on the correct association of sequences to gene families and proteins with known function, using the available literature as a reference.
the annotation is facilitated since each of the genes is associated with pubmed identifiers  stored in the pescador tool, which are important references for the related orthologs.
the search for functional information for the d. melanogaster orthologues revealed the involvement of the genes in processes related to the embryonic development and was also a good validation for the clustering by seedserver, since all sequences from d. melanogaster that clustered to the initial human and mouse genes present an embryo development related function.
generation of correct clusters is essential for the correct determination of gene ancestry, but it is not the sole limiting factor.
sequencing of key organisms from taxonomic outgroups relative to the ones with complete genome sequences available will be a crucial source of sequences that will allow a revaluation of gene ancestry.
meanwhile, additional sequences clustered by software  and database enrichment  improve the inspection of ancestry.
determination of the ancestry for the genes in the preimplantation pathway was nonetheless a central analysis, given the expectancy that this pathway would be mainly formed by more contemporary components.
our data suggest that an ancient fraction of the pathway including nanog and sox <dig> originated before chordata, whereas a modern fraction including oct <dig> and lif has appeared near the origin of eutheria, the placentary organisms.
thus, an important transcriptional pathway comprising ancient and modern members has been characterized with text mining, and homologues search with seedserver promptly allowed lca determination.
CONCLUSIONS
generation of regulatory pathways through text-mining tools allows integration of data generated by previous studies for a more complete view of a biological process.
if the genes present in this pathway are associated with clusters of orthologues this information is added to the pathway making the visualization of the same process available for different organisms.
the analysis of orthology also permits determination of the ancestry of the genes involved in the process leading to a better understanding of the evolution of such process.

 bioinformatics tools currently allow research to focus on the integration of large-scale data generated by sequencing, differential expression analysis, gene interaction studies and others. several initiatives exist to organize this knowledge in secondary databases, thus allowing easier access and visualization. databases containing interaction information are a good source for novel research. ihop allows users to tag gene names of interest and browse through the related pubmed literature with highlighted keywords. another interaction database is string, which contains physical interactions and functional associations between proteins and integrates data retrieved from literature, genomic context, large scale experiments and conserved co-expression. text-mining, therefore, has a fundamental role in these tools and allows access to interactions spread throughout the literature. the extraction of biological events from literature through text-mining tools is essential to not only update the interaction databases but also for the creation and annotation of pathways. metabolic and regulatory pathways are an example of organized knowledge that allow a better visualization of a complex system and can be found in databases such as ipath, biocyc or kegg pathways. when orthology information is added to pathways, the same process can be represented in different organisms. orthology is also an important tool for sequence annotation. current orthologue databases such as cog and kog, eggnog, orthomcl and kegg orthology all provide a good source for manually curated clusters of orthologues defined for organisms with complete genomes. we developed a procedure to enrich the cog database with uniref clusters from the uniprot database, creating the uecog database. recently, a similar procedure was applied to the kegg orthology database creating the enriched ueko database. the available tools described raise the possibility of integrating current information and generating complex regulatory pathways. previous publications individually reported the regulatory interactions that control preimplantation embryo development. however, a complete preimplantation development regulatory pathway has never been built. in humans, the preimplantation phase of embryonic development is a period of approximately six days after fertilization prior to attachment of the embryo to the uterine wall. implantation can occur before or in the seventh embryonic day, a time during which the uterus is receptive. mammalian embryonic development has been thoroughly studied in mice and the blastomeres remain totipotent, able to generate any other cell, up to the eight-cell stage, unlike other animals. after fertilization, successive cleavages take place during the first two days of development, resulting in the eight-cell embryo. the next stage of development is called the morula stage. an increase in cell-cell contact results in formation of a compacted morula. the subsequent divisions increase the complexity of the embryo and cells may be located on the inside, surrounded by other cells, or on the outside, in contact with the environment. the identification of the initial cells for each lineage has shown that the trophectoderm is derived mostly from the outer cells, whereas the inner cells give rise to the inner cell mass. later, the icm divides into the primitive endoderm and the epiblast. during the differentiation of the te from the icm, the blastocoel is formed through a process of cavitation. the embryo is called a blastocyst when all three structures are present. twenty-four hours after blastocyst formation occurs, the last stage of preimplantation development takes place when the pe differentiates from the icm. the three lineages thus formed in preimplantation development present different fates during subsequent embryonic development. while the epiblast, which forms from the icm following implantation, is still undifferentiated and will give rise to the fetus itself, the trophectoderm will become the fetal portion of the placenta and the primitive endoderm will form the yolk sac. complex regulatory processes such as animal development are a result of the interaction of many different gene products and elements that control the expression of these genes. traditional experiments that determine the function of one or a few genes are essential, but do not result in a comprehensive view of complex systems. a complex regulatory network should be able to portray specific and general aspects of development, such as the embryonic fate of certain cells. in this work, we noticed the absence in databases of a pathway describing the preimplantation phase of embryo development and sought to develop the given pathway using text-mining tools, complementing it with orthology information. the resulting pathway comprises genes and the interactions between them. clusters of orthologous groups were generated for each gene represented in the pathway and provided the necessary information to determine the last common ancestor. this determination revealed that the preimplantation development pathway is an ancient chordata pathway with addition of modern elements throughout evolution.  text-mining initially, we used the pubmed platform to search for articles related to the embryo preimplantation development and obtained entries as a result. to obtain a more efficient set of articles with relevancy to our work, the result entries were submitted to medlineranker. this software computes discriminating words by comparing a set of user selected abstracts indicated as highly relevant to a background set and then scores any abstracts in terms of their content of those discriminating words. after the classification, we selected the top abstracts for further analysis, which presented p-value lower than and by manual inspection provide large amount of information when uploaded in pescador. since human and mouse embryo development are highly similar, it was plausible to use abstracts from work on both organisms as source of information for the preimplantation pathway construction, paying attention to any possible conflict. using these highly informative abstracts as our input, pescador an online platform for friendly operation of the laitor software ) was used for tagging of gene names and biointeractions extraction from each abstract. as a result, gene names were tagged and type biointeractions were highlighted as well as other informative biointeractions. biointeractions are classified by laitor as type when in the same sentence the software encounters a gene name, a biointeraction word and another gene name, in that order. from these tagged abstracts we manually curated the information and constructed the pathway for the preimplantation embryo development describing genes and numerous interactions between them during the early developmental stages, trophectoderm differentiation from the inner cell mass and posterior extraembrionary endoderm differentiation. a sample abstract tagged by pescador and the manual extraction of the information it contains is exemplified in figure the pathway shown in figure was constructed according to kgml. the large decrease in the initial number of genes tagged in the abstracts is mainly due to redundancy between abstracts and also to genes tagged in type and biointeractions, which not always result in pathway building information. preimplantation pathway the pathway obtained after the analysis of all the abstracts from the pescador output is represented in figure and the regulations are reviewed below. first embryonic cleavages the oncogene c-myc is an important transcriptional regulator and its expression is observed in the initial stages of development, where it is present in embryonic cells until the morula stage and repressed thereafter. two additional genes recently associated with these early developmental stages are boris and ecsa. boris is involved in early development following fertilization and soon afterwards repressed, and ecsa, expression begins in the blastocyst exclusively in the cells of the inner cell mass. the presence of these genes was compared to the expression pattern of the oct transcription factor, which is present in the early cleavages, repressed after this initial stage, and then its expression is afterwards stimulated again in the blastocyst. the expression of the gametogenesis associated gene gse was also recently identified in cells of the early embryo; later this protein is found only in the icm, suggesting a role in the specification of cell lineage. methylation patterns and correct preimplantation development genomic methylation patterns in mammalian cells depend on dnmt. in the mouse, an embryo-specific variant called dnmt1o is expressed in the early stages of development. in the 8-cell stage this protein relocates to the cell nucleus where maintains essential methylation patterns, allowing embryos to complete early developmental events. it was recently shown that the inability of dnmt1o to properly relocate not only results in a developmental arrest at the 5- cell stage, but is also responsible for the downregulation of five genes involved in the formation of gap and tight junctions. these junctions are crucial for early processes such as compaction of the 8-cell embryo and cavitation of the blastocoel. te versus icm dichotomy: key role of lats controlling tead co-activator yap cells destined to become part of the icm are marked by repression of two genes and by upregulation of sox. in these cells, the major pluripotency transcription factors, including nanog and oct remain active due to the expression of an important player and member of the hippo signaling pathway: lats. this serine/threonine protein kinase is responsible for phosphorylating yap, leading to its cytoplasmic localization and thus preventing its association with the transcription factor tead triggering te differentiation: tead4/yap target cdx to repress nanog and oct4 conversely, in the outer cells that will differentiate and form the trophectoderm, yap is unphosphorylated, remains in the nucleus and associates with tead leading to the activation of cdx a key repressor of nanog and oct. repression of oct and nanog transcription by cdx then releases the inhibition that these two key factors were exerting on many different genes, in turn activating these targets. activation of cdx requires release from basal repression; nanog and oct repress basal levels of cdx and induction of higher levels of cdx by tead4/yap overcomes this repression, allowing cdx to play its role. tead was also recently determined to activate another trophectoderm differentiation factor, gata, which acts alongside cdx and affects transcription of a number of genes independent of cdx the tead4-dependent activation of gata seems to be independent of yap, suggesting tead interacts with another partner as well as yap. also required for high level expression of cdx in trophectoderm cells is the cell motility protein arp3; experiments with complete knockdown of this protein show trophoblast cells unable to develop properly, possibly undergoing apoptosis as a result of loss of cdx. the tgfbeta pathway is another important pathway for trophectoderm differentiation; tgfbeta signaling is stimulated by bmp which leads to the activation of smad proteins. these proteins can also stimulate transcription of cdx, and bmp is known to inhibit id an inhibitor of differentiation, and to activate hand which is involved in trophoblast cell differentiation. in the absence of oct and nanog the downregulation of oct in the outer cells of the embryo leads to the activation of a positive regulator of te cell fate, eomes, which is also a possible cdx target. the subsequent differentiation of these cells into trophectoderm is accompanied by the expression of several genes, such as the glycoprotein psg and the marker krt psg and krt expression are among the first signs that a blastomere has lost its totipotent competence, prior to any visible differentiation. removal of oct4-dependent repression also results in activation of genes such as etif2b and rps, allowing these cells to engage in an intense translation routine. knockdown studies targeting oct also show that it represses the expression of gcm which is normally placenta specific, and of the hcg hormone’s beta chain. concurrently, nanog downregulation allows the expression of a number of genes associated with both trophectoderm and extraembryonic endoderm. these latter genes will in turn initiate the formation of tissues such as the primitive endoderm, a component of the yolk sac. from the early blastocyst stage on, desmosomes are assembled in the trophectoderm in response to desmocollin, which is also not expressed in the icm. thus, tead4/yap activation of cdx accompanied by the subsequent repression of nanog and oct describes a scenario for the te differentiation. underneath the maintained activation of oct and nanog back in the icm, the main pluripotency genes remain active and form a complex regulation pathway. recently it was discovered that transcription of nanog is further stimulated by the presence of compounds such as retinol. klf klf and klf exert a redundant role in the activation of nanog. these krüppel-like factors were described as essential for the maintenance of pluripotency. indeed, klf was already known for this role and is commonly used in reprogramming of differentiated cells into induced pluripotent stem cells. however, only the simultaneous depletion of klf and results in the differentiation of stem cells, indicating functional redundancy. other proteins known to activate nanog include the two other main pluripotency regulators, oct and sox. the estrogen receptor esrrb is also reported to be involved in the activation of nanog by oct and sox. conversely, nanog can activate oct, and esrrb is necessary to maintain oct promoter activity. each of the three key factors, oct sox and nanog, also act as self-activators, e.g. the partners oct and sox bind and activate oct transcription. another key transcription factor involved in the maintenance of cell pluripotency is sall. sall binds to the conserved regulatory region in the pou5f distal enhancer and activates its transcription. studies with microrna interference of sall show that the loss of this factor leads to reduction of oct mrna levels and significant expression of cdx in the icm. b-myb, a gene expressed in proliferating cells, is also a positive regulator of oct and studies report early differentiation of icm in the absence of b-myb. the notch signaling pathway is a conserved pathway that is involved in cellular communication processes and correct cell fate decisions that also has a role in icm development. nle protein, a direct regulator of this pathway, is essential for survival of the icm. another protein associated with development and survival of the icm is tbn, whose absence promotes cell apoptosis in the icm. expression of the platelet and endothelial cell adhesion molecule was detected by immunofluorescence confocal microscopy in the blastocyst and restricted to the icm cells. subsequently, pecam remains only in the pluripotent epiblast cells, disappearing the moment these cells undergo differentiation, and indicating a new role for this molecule during embryo development. activation, but with moderation other control pathways maintain expression of these genes at a steady-state concentration and balance these many mechanisms for activation and upregulation of transcription. a complex regulation feedback loop consists of foxd nanog and oct. to keep oct and nanog expression within steady-state levels, these three genes interact so that expression of nanog activates foxd and oct but not above steady-state levels due to oct exerted repression; and foxd and nanog activate oct expression but not above steady-state levels due to oct self-repression. dax is an orphan nuclear hormone receptor recently identified as a repressor of oct transcription. dax expression was also capable of reducing nanog and rex expression. assays show that dax binds to oct and abolishes its dna binding activity, thus decreasing the transcription of nanog and rex targets of oct activation. another repressor in the icm is tcf a wnt signaling pathway effector. tle and ctbp are key partners of tcf in mediating this repressive effect. tcf binds to and represses the oct promoter, and this repressive effect requires both the groucho and ctbp interacting domains of tcf. tcf also limits the steady-state levels of nanog mrna, protein, and promoter activity in self-renewing embryonic stem cells; the tcf groucho domain is involved in this repression. thus, tcf is critical for maintaining the appropriate levels of both oct and nanog in escs. experiments show that loss of tcf by rna interference knockdown blocks the ability of escs to differentiate, emphasizing the importance of this interaction. downstream of oct and sox2 oct activates embryonic stem cell-specific gene, which encodes an rna binding protein present in the icm that is responsible for regulating several specific target transcripts. oct and sox are also responsible for the regulation of the fibroblast growth factor. expression of fgf therefore, requires the combined activity of these two transcription factors that bind to adjacent sites on the fgf enhancer dna region. once expressed, the fgf protein can interact with its receptor fgfr and activate icm and adjacent te cell proliferation, activating extraembryonic endoderm cells as well in later stages. several other genes with important functions in embryonic development are also targets of oct4-dependent activation. these include growth factor tdgf growth inhibitor sap regulator of nonsense transcripts rent two proteins involved in stem cell self-renewal dppa and dppa, anterior visceral endoderm markers lefty and lefty surface antigen thy and other genes encoding proteins involved in specialized cellular processes and hypothetical proteins. the master regulation exerted by sox and oct during mammalian embryogenesis is believed to operate through their cooperative binding to dna regulatory regions composed of adjacent hmg and pou motifs. exemplifying this arrangement, dppa is one such gene with the presence of an hmg/pou cassette in its promoter region. downstream of nanog and stat3 activation of jak/stat pathway also has an important contribution to pluripotency. in mice, the lif/stat pathway for maintenance of cell pluripotency comprises lif and lif receptor, which deliver intracellular signaling through stat stat a signal transducer and activator of transcription is activated by the jak kinase and binds to several promoters inducing transcription of pluripotency related genes. nanog and stat were found to bind to and synergistically activate stat3-dependent promoters. nanog also functions as a transcriptional inhibitor to nfκb, a factor known to have pro-differentiation activity. nanog is also responsible for smad repression, thereby preventing bmp4-induced differentiation through the tgfbeta signaling pathway, for which smad is a key signal transducer. extraembryonic endoderm differentiation from icm cells prior to embryo implantation one more differentiation takes place. certain cells from the icm give rise to the primitive endoderm, the first morphologically distinct cell type of the extraembryonic endoderm. the extraembryonic endoderm comprises the primitive, parietal and visceral endoderm components and will become the yolk sac during posterior development stages. wnt was recently identified as an inducer of primitive endoderm and this induction is accompanied by translocation of beta-catenin and snail to the nucleus. this study also showed that up-regulation of protein kinase a induces markers of parietal endoderm. another wnt family member, wnt9a, is expressed only in icm cells that surround the bastocoel and induces repositioning of the cells expressing gata which is necessary for formation of primitive endoderm. sox plays a major role in parietal endoderm differentiation. through studies with short interfering rna molecules, it was established that sox is responsible for transcription induction of gata and gata. individual or combined silencing of sox gata and gata result in suppression of cell shape changes and production of laminin-, characteristic changes present in parietal endoderm differentiation. gata was previously identified as a transcription factor responsible for the activation of fgf. sox also activates the fgf promoter. conversely, sox can negatively modulate the gata4-dependent activation of fgf which is supported by the role of this factor in icm pluripotency. another sox family member, sox is responsible for the differentiation of the extraembryonic endoderm in the final steps of preimplantation development. the runx factor is associated with the expression of sox and is also specific for the extraembryonic endoderm. hnf is a transcription factor specific of the extraembryonic endoderm with subsequent roles in post-implantation development and organogenesis. its expression may result from bmp4-induced differentiation. finally, the dab protein is indispensable for the development of visceral endoderm; though its exact role is still not established, it is perhaps related to correct cell positioning. the expression of cer a marker of the anterior visceral endoderm, commences before embryo implantation in the subset of cells that comprise the primitive endoderm. this ancestral population includes both cells expressing cer together with cells in which cer expression begins after implantation and formation of the ave. search for homologues to establish an ortholog database and provide sequence information to the genes contained in the preimplantation pathway, aminoacid sequences corresponding to the human and mouse gene products were used as seed for the software seedserver. in fact, only the uniprot identifier for these proteins is necessary to execute seedserver - gene symbols were verified in the ncbi gene database and converted to the corresponding geneid, and the desired identifiers were obtained afterwards from the uniprot database. for each gene a cluster of homologues was generated comprising from to sequences. the recruited sequences contained in each cluster can be swiss-prot annotated or unrevised trembl sequences. in total, 25% of the cluster sequences are swiss-prot, the great majority of clusters being comprised of trembl sequences. the search for homologues through seedserver provides therefore a large amount of candidates for manual curation in swiss-prot. furthermore, seedserver can recruit sequences from organisms without a complete genome due to its use of ueko and bidirectional best hit searches conducted by seedlinkage, and in fact only 27% of the sequences present in all clusters are from organisms with a complete genome. the ortholog clustering by seedserver was only performed for genes that had a corresponding swissprot annotated gene product to be used as seed, therefore hrscp and dppa which are described in the pathway, did not go through this analysis. pathway ancestry we then focused on the putative origin of these genes, determining which clade in the human lineage shares each gene. the generation of ortholog clusters allowed for the determination of the last common ancestor for each of the genes in the pathway. figure shows the genes according to their origin. genes were arbitrarily considered ancient for this analysis if their last common ancestor originated before the divergence of the clade euteleostomi and are coloured grey. genes with a lca belonging to the clade euteleostomi or originated after divergence of euteleostomi are considered recent genes and are coloured blue. ancient origin genes with an ortholog in drosophila melanogaster are marked with a red asterisk. this arbitrary classification was meant to attract attention to the two key pluripotency controlling genes, nanog and oct. the graph shown in figure represents the distribution of all the genes in the pathway according to their origin respect to clades of the human lineage. it may be observed that a large quantity of genes originates in certain periods as seen in eumetazoa, coelomata, euteleostomi and eutheria. the reasons for this wavelike origin need to be further analysed. on the other hand, the apparent origin of complex structures, that characterize all descendents from a certain moment of evolution, might have occurred simultaneously to the specialization of gene groups. the coverage of genomic sequences in the database is far from homogeneous and can influence the shape of this graph. in any case, the pattern observed agrees with the expansion of protein families related to stem cell markers observed in the ray-finned fish, that is, after divergence of the euteleostomi. furthermore, we searched for functional information related to the d. melanogaster orthologues in order to determine if these functions are somehow similar or related to the functions of the corresponding pathway genes. this was done through a second text mining approach similar to the first and from the information recovered a secondary pathway was generated simply to illustrate the ortholog genes and their relative functional roles. the regulatory pathways in which these genes are involved show us that these genes are all related to some part of drosophila embryo development, some of them with highly conserved functions still observed in the preimplantation pathway described. an example is the hippo signalling pathway, which is extremely conserved, showing wts phosphorylating yki; this modification prevents yki interaction with sd. the correlation between the human gene names and corresponding d. melanogaster ortholog names can be found in additional file and also the pmid reference for the gene function in drosophila development. discussion the use of text-mining tools for the generation of regulatory pathways is an effective approach and it is important for the current interest of gathering data related to an organism or biological process. the search for information related to a specific concept such as “preimplantation development” resulted in the selection of data related to this process only. when other tools such as ihop and string are used for the search of biointeractions, it is necessary to know the names for the genes you are interested on and the information is then retrieved. moreover in the case of ihop, the information retrieved consists of a large list of papers related to the gene of interest, which need to be manually analysed to extract the information related to the specific process. in the case of string, the result of a query is a network of direct associations to other genes, which can be activations, repressions, or unknown, but for which it is not possible to perform a search restricting the query to a specific process for which you seek to determine the involvement of a given gene. the approach described in this work summarized in figure allows the researcher to initiate the study of a pathway without knowing exactly the genes involved, simply by selecting the published information related to the process of interest. the manual curation required to create a pathway through this approach is significantly smaller. however, the verification of all the interactions highlighted by the tool is essential. text-mining is not able to eliminate the selection of false interaction pairs; in the case of laitor, the type and interactions can present genes with no association specified in the text. the text-mining data contribute the complete description of the pathway in the form of a literature review, a necessary step for the validation of the regulations represented, and for the inclusion of the pathway in a specific database, such as kegg pathway. the establishment of this procedure for pathway generation allows future work to enlarge the knowledge on subjects still not approached, such as regulatory pathways for several types of cancer, mechanisms of pathogen resistance in plants and response to abiotic stresses in plants, among other themes of interest. the inclusion of the preimplantation pathway in databases such as the kegg database will allow automatic annotation for several other organisms, as it is usually done in this database. concurrently, a laboratory with a specific interest can promptly build a similar pathway for its local use. from the genes present in the pathway, do not possess entries in kegg orthology and would constitute important additions. considering that the contribution of kegg for the sequence recruitment in the seedserver clusters is only 25% of the total number of sequences, some organisms evolutionarily divergent from the ones represented in kegg begin to play a more relevant role for a more efficient annotation of new sequences. it is relevant to stress that only the seedlinkage and ueko components of seedserver are capable of clustering sequences proceeding from organisms without a complete genome project. moreover, linkage of recruited to seed sequences are verified with psi-blast. another important contribution from the ortholog clustering by seedserver is the identification of candidates for swiss-prot annotation. swiss-prot annotation depends on the correct association of sequences to gene families and proteins with known function, using the available literature as a reference. the annotation is facilitated since each of the genes is associated with pubmed identifiers stored in the pescador tool, which are important references for the related orthologs. the search for functional information for the d. melanogaster orthologues revealed the involvement of the genes in processes related to the embryonic development and was also a good validation for the clustering by seedserver, since all sequences from d. melanogaster that clustered to the initial human and mouse genes present an embryo development related function. generation of correct clusters is essential for the correct determination of gene ancestry, but it is not the sole limiting factor. sequencing of key organisms from taxonomic outgroups relative to the ones with complete genome sequences available will be a crucial source of sequences that will allow a revaluation of gene ancestry. meanwhile, additional sequences clustered by software and database enrichment improve the inspection of ancestry. determination of the ancestry for the genes in the preimplantation pathway was nonetheless a central analysis, given the expectancy that this pathway would be mainly formed by more contemporary components. our data suggest that an ancient fraction of the pathway including nanog and sox originated before chordata, whereas a modern fraction including oct and lif has appeared near the origin of eutheria, the placentary organisms. thus, an important transcriptional pathway comprising ancient and modern members has been characterized with text mining, and homologues search with seedserver promptly allowed lca determination.  generation of regulatory pathways through text-mining tools allows integration of data generated by previous studies for a more complete view of a biological process. if the genes present in this pathway are associated with clusters of orthologues this information is added to the pathway making the visualization of the same process available for different organisms. the analysis of orthology also permits determination of the ancestry of the genes involved in the process leading to a better understanding of the evolution of such process. 
the generation of regulatory pathways through text-mining tools allows the integration of data generated by several studies for a more complete visualization of complex biological processes. the clustering of homologous genes together with determination of the ancestry leads to a better understanding of the evolution of such process. we noticed the absence of a chart or pathway describing the well-studied preimplantation development stages; furthermore, not all genes involved in the process have entries in kegg orthology, important information for knowledge application with relation to other organisms. these homologues allowed the determination of the last common ancestor for each gene and revealed that the preimplantation development pathway consists of a conserved ancient core of genes with the addition of modern elements. the integration of sequencing and gene interaction data and subsequent generation of pathways and networks contained in databases such as kegg pathway is essential for the comprehension of complex biological processes.
 the integration of sequencing and gene interaction data and subsequent generation of pathways and networks contained in databases such as kegg pathway is essential for the comprehension of complex biological processes. we noticed the absence of a chart or pathway describing the well-studied preimplantation development stages; furthermore, not all genes involved in the process have entries in kegg orthology, important information for knowledge application with relation to other organisms.  in this work we sought to develop the regulatory pathway for the preimplantation development stage using text-mining tools such as medline ranker and pescador to reveal biointeractions among the genes involved in this process. the genes present in the resulting pathway were also used as seeds for software developed by our group called seedserver to create clusters of homologous genes. these homologues allowed the determination of the last common ancestor for each gene and revealed that the preimplantation development pathway consists of a conserved ancient core of genes with the addition of modern elements.  the generation of regulatory pathways through text-mining tools allows the integration of data generated by several studies for a more complete visualization of complex biological processes. using the genes in this pathway as “seeds” for the generation of clusters of homologues, the pathway can be visualized for other organisms. the clustering of homologous genes together with determination of the ancestry leads to a better understanding of the evolution of such process. 15- november 6th international conference of the brazilian association for bioinformatics and computational biology ouro preto, brazil 