Physiol. Genomics Genetics/Genomics of Vascular Disease Workshop
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Physiol. Genomics 25: 405-413, 2006. First published March 22, 2006; doi:10.1152/physiolgenomics.00307.2005
1094-8341/06 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental Tables and Figures
Right arrow All Versions of this Article:
25/3/405    most recent
00307.2005v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (12)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Larson, J. H.
Right arrow Articles by Lewin, H. A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Larson, J. H.
Right arrow Articles by Lewin, H. A.
Received 15 December 2005; accepted in final form 23 February 2006.
Physiological Genomics 25:405-413 (2006)
1094-8341/06 $8.00 © 2006 American Physiological Society

Discovery of eight novel divergent homologs expressed in cattle placenta

Joshua H. Larson, Charu G. Kumar, Robin E. Everts, Cheryl A. Green, Annelie Everts-van der Wind, Mark R. Band and Harris A. Lewin

Department of Animal Sciences, University of Illinois at Urbana-Champaign, Urbana, Illinois


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
Ten divergent homologs were identified using a subtractive bioinformatic analysis of 12,614 cattle placenta expressed sequence tags followed by comparative, evolutionary, and gene expression studies. Among the 10 divergent homologs, 8 have not been identified previously. These were named as follows: cattle cerebrum and skeletal muscle-specific transcript 1 (CSSMST1), cattle intestine-specific transcript 1 (CIST1), hepatitis A virus cellular receptor 1 amino-terminal domain-containing protein (HAVCRNDP), prolactin-related proteins 8, 9, and 11 (PRP8, PRP9, and PRP11, respectively) and secreted and transmembrane protein 1A and 1B (SECTM1A and SECTM1B, respectively). In addition, two previously known divergent genes were identified, trophoblast Kunitz domain protein 1 (TKDP1) and a new splice variant of TKDP4. Nucleotide substitution analysis provided evidence for positive selection in members of the PRP gene family, SECTM1A and SECTM1B. Gene expression profiles, motif predictions, and annotations of homologous sequences indicate immunological and reproductive functions of the divergent homologs. The genes identified in this study are thus of evolutionary and physiological importance and may have a role in placental adaptations.

molecular evolution; bovine genome; comparative genomics


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
THE PLACENTA DEMONSTRATES remarkable variation among eutherian mammals, and it is believed that the vast diversity of placental architectures seen in different mammalian lineages is the result of adaptive evolution (8, 47). Two major selective forces likely driving adaptation of the placenta are fetal-maternal conflict over nutrient resource allocation and maintenance of immunologic tolerance of the fetal semi-allograft (45, 66). Cattle possess a synepitheliochorial placenta arranged as placentomes consisting of a fetal cotyledon interdigitating with a maternal caruncle. The cattle embryonal trophoblast contains granulated binucleate cells at the microvillar junction with cytoplasmic projections that penetrate the maternal uterine endometrium (70, 72). Ruminant-specific genes with expression localized to the trophoblast have been identified, including the interferon-{tau} genes (54), placental lactogens (PLs; Ref. 6), pregnancy-associated glycoproteins (PAGs; Ref. 71), prolactin-related proteins (PRPs; Refs. 46, 67, 72), and trophoblast Kunitz domain proteins (TKDPs; Ref. 42). The expression of these genes at sites of direct contact between fetal and maternal cells suggests that they represent molecular adaptations mediating maternal-fetal interaction during pregnancy (8, 55).

Comparisons between the whole genome sequences of human, mouse, and rat revealed extensive lineage-specific expansions in gene families associated with reproduction, immunity, chemosensation, detoxification, and proteolysis (48, 52). The identification of orthologous relationships for ~80% of genes in the human and mouse genomes indicates that lineage-specific gene family expansion likely represents a major route of mammalian evolution. These lineage-specific paralogs are believed to represent adaptive innovations on the basis of evidence that positive selection has operated on a subset of proteins functioning in reproduction and immunity (48). It is likely that the Cetartiodactyla experienced similar numbers of lineage-specific gene family expansions, given that cetartiodactyls diverged from the common ancestor of the Rodentia and primates ~94 million years ago (MYA), before the rodent-primate divergence ~87 MYA (62). It is important to determine which genes or gene families have been subjected to diversification and what roles these genes played in generating adaptive cetartiodactyl phenotypes.

In the present study, eight new divergent homologs expressed in cattle placenta were identified using a subtractive bioinformatic approach previously developed for the identification of the cattle major histocompatibility complex class I-like gene family A genes (38). The combined approaches of sequence alignment and comparative mapping were used to establish homologous relationships. Biological functions were inferred from the comparison of sequence motif predictions and microarray expression profiles to existing annotations for homologous proteins. Lastly, nucleotide sequence substitution analysis provided evidence for positive selection on two groups of divergent homologs.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
Sequences and bioinformatics.
More than 12,600 cattle expressed sequence tags (ESTs) from term placenta (dbEST library no. 3743) were interrogated by basic local alignment search tool (BLAST)N (1) to identify similarity to sequences in the GenBank nonredundant (February 15, 2005 release; cetartiodactyl sequences removed), dbEST (February 25, 2005 release; cetartiodactyl sequences removed), and human genome draft sequence (NCBI build 35.1) databases. Cattle ESTs with significant alignments with these sources were removed from the working set to enrich for highly divergent sequences. An E-value significance threshold of <e–5 was used for GenBank nonredundant and dbEST databases; E-value <e–10 was the significance threshold for human genome draft. Repetitive elements were masked using RepeatMasker (http://www.repeatmasker.org/).

The divergent 5'-cattle placenta ESTs were aligned by BLASTN to all other cattle GenBank sequences (excluding cattle placenta library no. 3743). Database sequences that aligned significantly were clustered with the working set of divergent 5'-ESTs using an in-house program that required 95% identity and 40-bp overlap. Consensus sequences were assembled using CAP3 (24). The BLASTN pipeline (Supplemental Table S1; available at the Physiological Genomics web site)1 was used to discard sequences with significant alignments.

Sequencing the 3'-ends of each cDNA clone with a divergent 5'-EST generated full-clone sequences. Bacterial stocks containing cDNA clones with divergent 5'-ESTs were cultured and rearrayed using a Biorobot 8000 (Qiagen, Valencia, CA), and template plasmid DNA was isolated according to the manufacturer's protocol. Cycle sequencing and electrophoresis were performed using an established protocol (38), substituting an oligo(d)T18V sequencing primer. PHRED (15) was used to remove low-quality sequence (PHRED score <20).

The BLASTN pipeline was implemented to discard clones with 3'-ESTs displaying significant similarity to sequences in the databases, using the significance thresholds described above. The remaining divergent 3'-ESTs were clustered with the divergent 5'-ESTs and all cattle ESTs in dbEST (excluding cattle placenta library no. 3743). The clusters were assembled, and clones with nonoverlapping 5'- and 3'-regions were sequenced by primer walking (38). The BLASTN pipeline was then used to remove full-clone sequences with similarity to database sequences.

The full-length divergent cDNA sequences were interrogated by BLASTX (2) against the GenBank coding sequence database (including cetartiodactyl sequences; January 9, 2005) to identify distant homologies. Cattle sequences that aligned to database sequences with E-value <e–10 were designated as divergent homologs. Full-clone sequences were conceptually translated using SIXFRAME (http://workbench.sdsc.edu/). Usage of all homology-related terms is according to Koonin (36).

Multiple alignments were constructed using CLUSTALX (27). Signal peptides, transmembrane (TM) domains, and conserved domains were predicted using PSORTII (49), TMPred (23), TMHMM 2.0 (http://www.cbs.dtu.dk), and the conserved domain database (44). Glycosylation, glycophosphatidylinositol (GPI) anchors, and additional predictions were carried out using NetNGlyc and NetOGlyc (http://www.cbs.dtu.dk), big-Pi predictor (10), and PredictProtein (56). Homology modeling was performed using Swiss-Model and Swiss-PdbViewer (17). Sequence data have been deposited with the DDBJ/EMBL/GenBank data libraries under the following accession numbers: AY563705AY563905, AY741372AY741374, BK005429BK005441, and CK393994CK394219.

Radiation hybrid mapping.
To confirm predicted homologous relationships, we performed radiation hybrid (RH) mapping and compared the results with positions in the human genome on the basis of comparative mapping information (14). The divergent homologs were mapped using a 5,000-rad cattle-hamster RH panel (69). Primer design, PCR, and mapping were carried out as previously described (Ref. 14 and Supplemental Table S2). Homologous positions in the human genome (NCBI build 35) were identified using the UCSC Genome Browser (31).

Microarray analysis.
To elucidate gene expression patterns and to gain insight into biological function, expression profiles were determined for 7 of the 10 transcripts using a 7,653-element cattle cDNA microarray and RNAs collected from placenta and 17 additional cattle tissues (Fig. 1 and M. R. Band, R. E. Everts, Z. L. Liu, D. E. Morin, J. U. Peled, S. L. Rodriguez-Zas, and H. A. Lewin, unpublished data; NCBI GEO: GSE3029). Total RNA was isolated from a Hereford cow term placentome using an established protocol (38). Twenty micrograms of total RNA and titrated RNA spiking controls (5) were used as templates for the synthesis of aminoallyl-labeled cDNA (aa-cDNA) using an indirect labeling protocol (20). First-strand aa-cDNA was purified and fluorescently labeled as previously described (68). The aa-cDNA synthesis and labeling procedures were also carried out on a standard reference RNA sample developed in our laboratory for microarray studies (13). Labeled cDNA samples were purified and hybridized to a microarray consisting of 7,653 cattle spleen- and placenta-derived cDNAs as described previously (13, 68). Slides were scanned using an Axon 4000B scanner, images were analyzed with GenePix 4.0 (Molecular Devices, Union City, CA), and data were normalized with GeneSpring 5.0 (Silicon Genetics, Redwood City, CA) using the approach of Band et al. (5). Microarray data have been deposited with the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) under the accession no. GSE3029.


Figure 1
View larger version (21K):
[in this window]
[in a new window]
 
Fig. 1. Microarray expression data for 7 divergent homologs (for details, see RESULTS), using a diverse panel of cattle tissues. Expression is shown as the log10 of the normalized ratio of signal intensities in each cattle tissue compared with a reference standard. The tissues are indicated as follows: CR, cerebrum; HE, heart; KI, kidney; LI, large intestine; LV, liver; MU, skeletal muscle; PL, placentome; RU, rumen; SI, small intestine; SK, skin. Data were excluded when transcript expression level did not exceed background by 3 standard deviations in a particular tissue. Significance of transcript expression level difference between tissue and reference standard was determined by t-test; associated P values are indicated: +P < 0.05, ++P < 0.01, +++P < 0.005, and ++++P < 0.001. Only tissues with significant effects are shown (data not shown for adrenal gland, cerebellum, lung, lymph node, spleen, testis, thalamus, and thymus).

 
Substitution analysis.
Ratios of nonsynonymous to synonymous substitutions (dN/dS, or {omega}) were determined for the divergent homologs using the phylogenetic analysis by maximum likelihood (PAML) software package (74) to identify evidence of positive Darwinian selection as a contributing factor in their sequence divergence. Sequences used for the substitution analyses included the following: cattle CCSMST1 (BK005429), chimpanzee LOC453815 (XM_510728), mouse BC003965 (XM_128528), cattle HAVCRNDP (BK005433), cattle predicted HAVCR1 (BF043390), cattle SECTM1A (BK005436), cattle SECTM1B (BK005437), cattle PRP1 (NM_174159), cattle PRP2 (M27239), cattle PRP3 (NM_174160), cattle PRP4 (M33269), cattle PRP5 (NM_212448), cattle PRP6 (X59504), cattle PRP7 (AB187564), cattle PRP8 (BK005439), cattle PRP9 (BK005440), and cattle PRP11 (BK005438). Signal peptide and TM regions were omitted. The YN00 program in PAML was used to estimate global dN/dS ({omega}t) for each group of aligned sequences using the method of Yang and Nielsen (75). Values of {omega}t >1.0 indicate that positive selection acted on the sequences analyzed (50). For the PRP genes, the CODEML program in PAML was used to identify variation in selection intensity. The data were modeled using maximum likelihood methods (76), and the results were compared to obtain a test statistic. Three comparisons were performed. Model M1, a neutrality model that constrained {omega} to be either 0 or 1, was compared with both M2, a selection model that added an additional {omega} ratio class estimated from the data, and M3, a selection model that used an unconstrained discrete distribution to model classes of {omega} ratios. Our analysis used three discrete classes for M3. We also compared M7, a continuous distribution neutrality model that estimates {omega} using a beta function limited to the interval from 0 to 1, with M8, a continuous distribution selection model that adds an additional class of sites with {omega} estimated from the data and not constrained to the interval between 0 and 1. A test statistic of twice the negative value of the difference between the log likelihood values generated under each model was compared with a {chi}2 distribution with degrees of freedom (df) calculated from the difference in the number of model parameters (M2 vs. M1, df = 2; M3 vs. M1, df = 4; M8 vs. M7, df = 2). Posterior probabilities for PRP sites under positive selection were generated under M8.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
A starting set of 12,614 term cattle placenta 5'-ESTs was reduced to 2,379 divergent 5'-ESTs after BLASTN subtraction and further reduced to 1,595 divergent sequences after BLASTN subtraction was implemented for the assembled consensus sequences (Supplemental Table S1). Sequencing of 3'-ESTs was successful for 1,180 clones with 5'-ESTs in the working set, and BLASTN subtraction of the 3'-ESTs reduced the number of putative divergent clones to 776. BLASTN subtraction of full-clone sequences reduced the working set to 368 high-probability divergent transcripts with no significant nucleotide similarity to sequences derived from organisms outside the Cetartiodactyla. Distant peptide homologies were identified for 10 transcripts using BLASTX (Supplemental Table S3). Conceptual translations of the 10 divergent cattle homologs were characterized with respect to predicted signal peptides, conserved domains, and TM regions (Table 1). Orthologous and/or paralogous relationships are proposed below for each cattle sequence on the basis of their alignment and comparative mapping data.


View this table:
[in this window]
[in a new window]
 
Table 1. Bioinformatic characterization of the divergent homolog protein sequences

 
CCSMST1.
A translation of the cattle transcript BK005429 shares 64 and 59% sequence identity with predicted proteins XP_510728 of chimpanzee and XP_128528 of mouse, respectively (Supplemental Fig. S1). The genes associated with these protein sequences were thus classified as putative orthologs, and orthology was subsequently confirmed by RH mapping (Supplemental Table S4). The cattle gene was named cattle cerebrum and skeletal muscle-specific transcript 1 (CCSMST1) on the basis of its transcript expression profile (Fig. 1). Alignment of the CCSMST1 transcript to the current draft of the cattle genome (build 2.0) showed identity to genomic contig NW_930368 assigned to BTA25, supporting the RH mapping data. The first 52 bases of the transcript make up the first exon, which contains the 5'-untranslated region (UTR) and the majority of the signal peptide. A 78-bp intron separates the first exon from a second exon containing the remainder of the coding sequence and 3'-UTR. Substitution analysis revealed {omega}t = 0.556 for the CCSMST1 orthologs, suggesting that nonsynonymous substitutions may be selected against in this gene.

CIST1.
A translation of the cattle transcript BK005441 shares 49 and 30% sequence identity with the predicted human CIST1 ortholog and chimpanzee protein XP_530016, respectively (Supplemental Fig. S2). The cattle ortholog was named cattle intestine-specific transcript 1 (CIST1) on the basis of its transcript expression profile (Fig. 1). Query of the cattle draft genomic sequence with the CIST1 transcript sequence revealed identity to genomic contig NW_931809 assigned to BTA7, supporting the RH mapping data (Supplemental Table S4). The genomic structure consists of a 175-bp exon consisting of the 5'-UTR and signal peptide, a 2,603-bp intron, a 390-bp exon, a 1,375-bp intron, a 152-bp exon, a 118-bp intron, and a 114-bp exon containing the remainder of the coding sequence and 3'-UTR, respectively. Cattle CIST1 contains an extracellular serine/threonine/proline (S/T/P)-rich mucin-like domain; these three residue types constitute 56% of this region. Nearly all of the serine and threonine residues in this region represent O-linked glycosylation sites with significant probability. CIST1 and its orthologs were not evaluated by substitution analysis methods because of extensive nonconserved regions in the coding sequences.

HAVCRNDP.
A translation of the cattle transcript BK005433 shares 48, 47, and 58% sequence identity with human hepatitis A virus cellular receptor 1 (HAVCR1), mouse T cell immunoglobulin and mucin domain-containing molecule 1 (TIM-1), and a conceptual translation of cattle EST BF043390, (the putative cattle HAVCR1 ortholog), respectively (Supplemental Fig. S3). Query of dbEST with the BK005433 transcript revealed numerous ESTs aligning with 100% identity and no evidence of alternate splicing. Query of the cattle draft genomic sequence with the BK005433 transcript showed identity to the unassigned cattle genomic contig NW_973071. However, RH mapping demonstrated HAVCRNDP to be located on BTA7, in proximity to IL12B (Supplemental Table S4). The genomic structure consists of a 210-bp exon containing 5'-UTR, a 5,146-bp intron, a 68-bp exon containing 5'-UTR and the majority of the signal peptide, a 1,954-bp intron, a 330-bp exon, a 2,779-bp intron, and a 239-bp exon consisting of the remainder of the coding sequence and 3'-UTR, respectively. The BK005433 transcript does not align with the cattle genomic sequence containing cattle HAVCR1, indicating that these loci represent distant paralogs. The cattle gene was thus given the name hepatitis A virus cellular receptor 1 amino-terminal domain-containing protein (HAVCRNDP). The HAVCRNDP molecule terminates at approximately the start of the mucin-like region found in the HAVCR1 proteins (16, 25). HAVCRNDP was significantly overexpressed relative to the reference standard in liver, small intestine, and skin tissues (Fig. 1). Substitution analysis of the cattle HAVCRNDP and HAVCR1 paralogs was performed on the conserved immunoglobulin (Ig) domain region only. The HAVCRNDP paralogs had {omega}t = 0.770, indicating that positive selection does not act on this gene or that nonsynonymous substitutions are weakly selected against.

PRP8, PRP9, and PRP11.
The PRP molecules described here were designated PRP8, PRP9, and PRP11 to conform with the revised PRP nomenclature (Refs. 9, 58; Fig. 2). Cattle PRP8, encoded by transcript BK005439, is most closely related to PRP1 (98% identity), whereas PRP9, encoded by transcript BK005440, is most closely related to PRP5 (97% identity), and PRP11, encoded by transcript BK005438, is most closely related to PRP3 (85% identity). Each of the newly described PRP transcripts have been uniquely localized to genomic sequence with 100% identity (PRP8 to AAFC01015348, PRP9 to AAFC01043996, PRP11 to AAFC01039561), indicating that these sequences represent paralogs. However, definitive proof of paralogy or allelism for all the PRPs awaits complete sequencing of this region in the cattle genome. PRP genes have not been identified in the human genome; however, the cattle PRP genes on BTA23 (Supplemental Table S4) share a conserved chromosome position with the human prolactin (PRL) locus on HSA6p22.3 (9). PRP8 was significantly overexpressed relative to the reference standard in heart, skeletal muscle, and placentome tissues, and PRP11 was significantly overexpressed relative to the reference standard in placentome tissue (Fig. 1). There was no sequence element on the array corresponding to PRP9, so its expression pattern could not be evaluated.


Figure 2
View larger version (32K):
[in this window]
[in a new window]
 
Fig. 2. Protein alignment of the cattle prolactin-related protein (PRP) family including PRP1 (NP_776584), PRP2 (AAA30728), PRP3 (NP_776585), PRP4 (AAA30740), PRP5 (NP_997613), PRP6 (CAA42092), PRP7 (BAD74059), PRP8 (DAA05611), PRP9 (DAA05612), and PRP11 (DAA05610). Alignment gaps are represented by periods, and aligned PRP sites sharing identity with cattle PRP1 are represented by dashes. The signal peptide region (sp) is represented by a solid underscore appearing above the alignment. Universally conserved cysteine and tryptophan residues are annotated with Cs and Ws, respectively, beneath the alignment. Predicted NH2-glycosylation motifs are represented by underlined text. Positively selected sites are indicated by asterisks above the alignment, and posterior probabilities associated with each positively selected site (see Table 2) are represented within the alignment by normal text (probability >0.90), bold text (probability >0.95), and italic and bold text (probability >0.99).

 
Substitution analysis showed the PRP paralogs to have {omega}t = 1.457, indicating that these genes have been subjected to positive selection (Ref. 50 and Supplemental Fig. S4). We found significant heterogeneity in selection intensity within the PRP molecules and identified 21 molecular sites appearing to be under a high degree of positive selection (Table 2, Fig. 2). Twenty of the twenty-one selected sites were mapped onto the tertiary structure of human PRL [Molecular Modeling Database (MMDB): 1N9D; Fig. 3 \.]. Fifteen selected residues are outwardly directed, indicating that the positive selection observed at these residue positions is likely exerted by an interaction with another molecule. There are no outwardly directed positively selected sites between PRP1 residues 117 and 180, a region with little evidence for positive selection and corresponding to helices 2 and 3 of the predicted PRP molecule on the basis of homology modeling (11, 30). Two relative "hotspots" of positively selected residues flank this region. The first hotspot, from residues 108 to 116, corresponds to the extended loop preceding helix 2. The second hotspot, from residues 181 to 198, is contained within the extended loop distal to helix 3. Thirty-two prolactin receptor-binding residues identified in the ovine PL molecule (11) were mapped onto the PRL molecule, and only four PL receptor-binding residues were found to overlap with the PRP sites under positive selection.


View this table:
[in this window]
[in a new window]
 
Table 2. Likelihood ratio tests of {omega} variation and identification of molecular sites under positive selection in the cattle prolactin-related proteins*

 

Figure 3
View larger version (48K):
[in this window]
[in a new window]
 
Fig. 3. Tertiary structure of human prolactin (PRL) [Molecular Modeling Database (MMDB): 1N9D] showing the spatial arrangement of homologous PRP residues under positive selection with >90% probability. The PRL backbone appears in blue; 20 of 21 PRP sites under positive selection are mapped onto the structure and appear as red space-filling amino acids.

 
SECTM1A and SECTM1B.
Cattle secreted and transmembrane protein 1A (SECTM1A), encoded by transcript BK005436, shares 32 and 27% sequence identity with the human and mouse SECTM1 proteins, respectively (Supplemental Fig. S5). Cattle SECTM1B, encoded by transcript BK005437, shares 30 and 26% sequence identity with the human and mouse molecules, respectively. SECTM1A and SECTM1B are 72% identical. The loci were named SECTM1A and SECTM1B to differentiate the two cattle paralogs, because only a single SECTM1 gene has been described in other species. SECTM1A contains a hydrophobic region from residues 167 to 184 preceded by a GPI anchor site at residue 163. Only partial sequences of SECTM1A and SECTM1B are present within the cattle genome draft sequence contigs (GenBank: NW_929537 and NW_929532, respectively), so determination of the genomic organization of the SECTM1 paralogs awaits completion of the cattle genome sequence. RH mapping demonstrated both SECTM1 genes to be located on BTA19, in proximity to tubulin-specific chaperone D (TBCD) (Supplemental Table S4). There were no sequence elements on the array corresponding to SECTM1A or SECTM1B, so expression patterns could not be evaluated. Substitution analysis revealed that the SECTM1 paralogs had {omega}t = 1.137, suggesting that these genes have been subjected to positive selection (50).

TKDP1 and TKDP4.
Cattle trophoblast Kunitz domain protein 1 (TKDP1) and TKDP4 were recently characterized (37, 42). The TKDP4 molecule identified in our study, encoded by transcript BK005435, is a splice variant lacking residues 32–59 of the previously described TKDP4 protein (42). TKDP genes have not been identified in the human genome; however, the TKDP genes on BTA13 (Supplemental Table S4) share a conserved chromosomal position with SPINLW1 on HSA20q13.12, which encodes a protein with a COOH-terminal Kunitz domain (53) and may share a common ancestry with the cattle TKDP genes. TKDP4 was significantly overexpressed relative to the reference standard in cerebrum, skeletal muscle, and placentome tissues, but TKDP1 did not show any preference in tissue expression pattern (Fig. 1). The TKDP paralogs were excluded from substitution analysis due to heterogeneous domain organization.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
The identification of genes that are highly divergent among closely related species is an important step in understanding the phenotypic changes that accompany adaptive evolution. The advent of high-throughput DNA sequencing and associated genome databases permit a more global search for those genes that diversify through copy number and primary sequence. In the absence of complete genome sequence, EST collections offer a means to search for divergent homologs as well as transcripts and genes that are unique to a given species (38). In the present study, cattle placenta was chosen as a source of ESTs, because the placenta is among the most phenotypically divergent mammalian tissues and has been shown to express novel paralogs such as the PAG and TKDP genes (42, 71). Our subtractive BLAST approach yielded eight previously unknown divergent homologs, the previously identified TKDP1 gene, and a novel splice variant of TKDP4. On the basis of cattle-human comparative mapping data (14), each divergent cattle gene shows conserved chromosomal position with respect to its predicted human ortholog (Supplemental Table S4). Gene expression analyses using microarrays proved to be highly effective in determining their tissue distributions and for primary functional annotations. For example, expression patterns were used to name CIST1 and CCSMST1, two genes with unannotated orthologs in other species. Our results support these computationally predicted orthologous gene models, thus demonstrating that our approach is also useful for annotating chimpanzee, human, mouse, and other genomes.

CCSMST1 appears to be expressed on the cell surface on the basis of its signal peptide and TM region (Table 1, Supplemental Fig. S1). Although there is insufficient evidence for positive selection over the entire molecule, the {omega}t is a stringent statistic and not likely to detect positive selection in subregions (12). Detection of heterogeneous selection intensity in the CCSMST1 genes is not possible at this time because of a lack of sufficient orthologs for a meaningful analysis (3).

The CIST1 signal peptide and TM domain indicate that the protein is expressed on the cell surface (Table 1, Supplemental Fig. S2). The high expression of CIST1 in large intestine, small intestine, rumen, and kidney tissues suggests that the protein has a gastrointestinal and genitourinary function (Fig. 1). Numerous O-glycosylation sites in the S/T/P-rich mucin-like domain indicate that CIST1 is highly glycosylated. Low peptide sequence conservation in the extracellular domain among CIST1 orthologs suggests that CIST1 function is mediated by its O-linked carbohydrate moieties.

HAVCR1 of humans and African green monkeys is an integral membrane protein bound by the hepatitis A virus (HAV) (16, 29). Cattle HAVCRNDP is the first HAVCR1 homolog described consisting of only a single secreted NH2-terminal HAVCR1-like Ig domain (Supplemental Fig. S3). Because of the absence of the mucin-like, TM, and intracellular HAVCR1 domains, the function of HAVCRNDP is likely to be specialized and dependent on the Ig domain. The Ig domain of African green monkey HAVCR1 is necessary and sufficient for HAV binding and neutralization, and the HAVCR1 mucin-like domain is necessary for viral genome uncoating (59, 60, 65). A secreted, soluble, HAVCR1-type Ig domain molecule could potentially bind a HAV-like virus, prevent viral uncoating, and limit viral infectivity. Thus HAVCRNDP may function as an antiviral molecule in cattle. Certain Picornaviruses (of which HAV is a member) infect cattle, including the apthoviruses that cause foot and mouth disease and are currently regarded as among the most economically important veterinary pathogens worldwide (35). The high expression of HAVCRNDP in liver and small intestine (Fig. 1) may neutralize such pathogens at sites of infectivity and clear them from the animal. Substitution analysis showed a {omega}t <1, but this does not eliminate the possibility that molecular subregions could be subject to positive selection. Given current data, we hypothesize that HAVCRNDP arose by duplication and divergence of the HAVCR1 locus during ruminant evolution as an adaptive response to an enteric pathogen(s).

The cattle PRP family contains seven previously described genes expressed in cattle placenta throughout gestation (26, 33, 57, 63, 64, 67, 73). PRP proteins are believed to function in implantation, placental development, and fetomaternal interaction throughout gestation, yet no known receptor has been identified (19, 3234, 46, 72). PRP1 interacts with {alpha}2-macroglobulin, which may facilitate presentation of PRP molecules to their molecular targets or allow PRP molecules to regulate the activity of other {alpha}2-macroglobulin-bound molecules (28, 34). Heterogeneity in glycosylation is the major structural difference among the PRPs and may ultimately contribute to the functional characteristics of each paralog (34, 77). High expression of PRP8 and PRP11 in term placenta is consistent with the expression profiles observed for other PRPs (Fig. 1). Surprisingly, PRP8 was highly expressed in heart and skeletal muscle.

Substitution analysis demonstrated that the PRPs have been subject to positive selection (Supplemental Fig. S4). Twenty-one high-probability (>90%) selected sites were identified, revealing two hotspots of positive selection flanking PRP helices 2 and 3, helices that may be under purifying selection (Figs. 2 and 3, Table 2). This structural association suggests that purifying selection resulted from a constraining influence such as that imposed by a receptor or ligand on a docking surface necessary for molecular interaction. The {alpha}2-macroglobulin molecule may exert this influence, as it performs multiple intrinsic and carrier molecule functions (28, 34) that likely constrain peptide sequence mutability. The positions of the positive selection hotspots next to this putative docking site suggest that the hotspot residues may introduce significant diversity at a functional site(s) in the complex. There appears to be no association between PRP sites under positive selection and receptor-binding sites identified in ovine PL.

Cattle SECTM1A and SECTM1B represent the first discovery of paralogous SECTM1 loci in any species. Tight linkage and high sequence similarity compared with human and mouse SECTM1 suggest that SECTM1A and SECTM1B originated by duplication of an ancestral SECTM1 locus in the cetartiodactyl lineage. In humans and mice, SECTM1 is a predicted TM molecule with an NH2-terminal Ig-like extracellular domain that interacts with CD7 (40, 61). The signal peptide and NH2-glycosylation sites suggest that cattle SECTM1A and SECTM1B are extracellular glycoproteins (Supplemental Fig. S5), although the GPI anchor may allow SECTM1A to function as a cell surface receptor or soluble ligand depending on lipase activity and accessibility (18). Given the apparent role of CD7 in the activation of {gamma}{delta} T cells (4, 7), the expression of SECTM1 molecules by term placenta may induce proliferation of {gamma}{delta} T cells to the elevated levels observed in the ruminant neonate (21, 41). The SECTM1 molecules may also modulate maternal immunity (22, 39).

Substitution analysis indicated that the cattle SECTM1 paralogs have been subjected to positive selection. The existence of two cattle SECTM1 proteins may represent the need for both secreted and membrane-associated CD7 ligands in cattle. Alternately, one of the cattle SECTM1 molecules may function as the CD7 ligand, while the other acquired novel function. The latter possibility is more consistent with the sequence divergence of the cattle molecules and the evidence for positive selection.

The TKDP family in cattle is a group of five molecules containing COOH-terminal Kunitz-type domains and novel NH2-terminal domains (37, 42). In the present study, we identified an alternately spliced TKDP4 transcript, suggesting that other TKDP transcripts may be alternatively spliced, as proposed for TKDP1 (42). The TKDP expression data support previous analyses showing that TKDP1 is no longer expressed in term placenta or in any other tissue and that TKDP4 is highly expressed in near-term placenta (37, 42). The expression of TKDP4 in cerebrum and skeletal muscle is of unknown functional significance. The TKDPs are hypothesized to regulate cellular differentiation by inhibiting growth factor proteinases or to protect the conceptus by inhibition of maternal T cell granzymes (37, 42). Cattle TKDP4 inhibits trypsin and plasmin in vitro, but other TKDPs showed no proteinase inhibition, suggesting either a high degree of target proteinase specificity or an alternate function (43). An additional hypothesis, formulated by MacLean et al. (43), is that TKDP molecules regulate the activity of voltage-gated ion channels similar to their Kunitz domain-containing homologs in snake venom. The expression of TKDP4 in cerebrum and skeletal muscle, two tissues reliant on the conduction of action potentials, adds credence to this hypothesis and suggests that the study of extraplacental expression of TKDP genes could provide clues to additional functions.

One consistent theme in our results is that the divergent homologs are predicted to be expressed extracellularly without exception, despite clear differences in their ancestry. The divergent paralogs show a greater degree of overt positive selection compared with the divergent orthologs, likely because single-copy orthologs are under a greater degree of functional constraint against nonsynonymous sequence substitution than are duplicated paralogous genes. The latter may perform redundant functions and create a more permissive environment for nonsynonymous substitution (51). More sophisticated analyses may be required to detect positive selection within small orthologous sequence groupings and probably require structural and functional data as a guide. Substitution analysis of the CIST1 and TKDP genes could not be performed because of the extensive diversity of domain sequences and organization, but their discovery as divergent homologs suggests that gross divergence in domain structure contributes unique species- or lineage-specific molecules and may represent an adaptive modality in mammals.

Although the functions discussed for the divergent genes are based on bioinformatic predictions and derived functional annotations of homologs, the two functional attributes prevalent in our discoveries are immunity and reproduction. These two biological functions represent major foci of genomic diversity identified in other mammalian species (48, 52). Functional characterization of the divergent proteins identified in this study is required to evaluate practical applications of our findings and is expected to provide insight into the molecular events that guided ruminant evolution.


    GRANTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
This work was made possible by a grant from the United States Department of Agriculture, National Research Initiative, Project No. AG2002-35205-11625. J. Larson was supported in part by a Cargill Research Excellence Fellowship in Animal Genomics.

While this paper was in review, the discovery of two novel prolactin-related protein paralogs, designated PRP8 and PRP9, was reported (Ushizawa K, Takahashi T, Hosoe M, Kaneyama K, and Hashizume K. Cloning and expression of two new prolactin-related proteins, prolactin-related protein-VIII and -IX, in bovine placenta. Reprod Biol Endocrinol 3: 68, 2005). To avoid confusion in nomenclature, the PRP8 gene reported in this study (GenBank accession no. BK005439) will be renamed PRP12. In addition, the PRP9 gene reported in this study (GenBank accession no. BK005440) will be renamed PRP13.


    ACKNOWLEDGMENTS
 
We thank John Morse, Chris Wright, and Laura Guest (W. M. Keck Center for Comparative and Functional Genomics, University of Illinois at Urbana-Champaign) for clone rearraying and sequencing services.


    FOOTNOTES
 
Article published online before print. See web site for date of publication (http://physiolgenomics.physiology.org).

Address for reprint requests and other correspondence: H. A. Lewin, 206A Edward R. Madigan Laboratory, 1201 W. Gregory, Urbana, IL, 61801 (e-mail: h-lewin{at}uiuc.edu).

1 The Supplemental Material for this article (Supplemental Tables and Figures) is available online at http://physiolgenomics.physiology.org/cgi/content/full/00307.2005/DC1. Back


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 

  1. Altschul SF, Gish W, Miller W, Myers EW, and Lipman DJ. Basic local alignment search tool. J Mol Biol 215: 403–410, 1990.[CrossRef][Web of Science][Medline]
  2. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, and Lipman DJ. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25: 3389–3402, 1997.[Abstract/Free Full Text]
  3. Anisimova M, Bielawski JP, and Yang Z. Accuracy and power of the likelihood ratio test in detecting adaptive molecular evolution. Mol Biol Evol 18: 1585–1592, 2001.[Abstract/Free Full Text]
  4. Aruffo A and Seed B. Molecular cloning of two CD7 (T-cell leukemia antigen) cDNAs by a COS cell expression system. EMBO J 6: 3313–3316, 1987.[Web of Science][Medline]
  5. Band MR, Olmstead C, Everts RE, Liu ZL, and Lewin HA. A 3800 gene microarray for cattle functional genomics: comparison of gene expression in spleen, placenta, and brain. Anim Biotechnol 13: 163–172, 2002.[CrossRef][Web of Science][Medline]
  6. Byatt JC, Shimomura K, Duello TM, and Bremel RD. Isolation and characterization of multiple forms of bovine placental lactogen from secretory granules of the fetal cotyledon. Endocrinology 119: 1343–1350, 1986.[Abstract/Free Full Text]
  7. Carrel S, Salvi S, Rafti F, Favrot M, Rapin C, and Sekaly RP. Direct involvement of CD7 (gp40) in activation of TcR {gamma}/{delta}+ T cells. Eur J Immunol 21: 1195–1200, 1991.[Medline]
  8. Crespi B and Semeniuk C. Parent-offspring conflict in the evolution of vertebrate reproductive mode. Am Nat 163: 635–653, 2004.[CrossRef][Web of Science][Medline]
  9. Dietz AB, Georges M, Threadgill DW, Womack JE, and Schuler LA. Somatic cell mapping, polymorphism, and linkage analysis of bovine prolactin-related proteins and placental lactogen. Genomics 14: 137–143, 1992.[CrossRef][Web of Science][Medline]
  10. Eisenhaber B, Bork P, and Eisenhaber F. Prediction of potential GPI-modification sites in proprotein sequences. J Mol Biol 292: 741–758, 1999.[CrossRef][Web of Science][Medline]
  11. Elkins PA, Christinger HW, Sandowski Y, Sakal E, Gertler A, de Vos AM, and Kossiakoff AA. Ternary complex between placental lactogen and the extracellular domain of the prolactin receptor. Nat Struct Biol 7: 808–815, 2000.[CrossRef][Web of Science][Medline]
  12. Endo T, Ikeo K, and Gojobori T. Large-scale search for genes on which positive selection may operate. Mol Biol Evol 13: 685–690, 1996.[Abstract]
  13. Everts RE, Band MR, Liu ZL, Kumar CG, Liu L, Loor JJ, Oliveira R, and Lewin HA. A 7872 cDNA microarray and its use in bovine functional genomics. Vet Immunol Immunopathol 105: 235–245, 2005.[CrossRef][Web of Science][Medline]
  14. Everts-van der Wind A, Larkin DM, Green CA, Elliott JS, Olmstead CA, Chiu R, Schein JE, Marra MA, Womack JE, and Lewin HA. A high-resolution whole-genome cattle-human comparative map reveals details of mammalian chromosome evolution. Proc Natl Acad Sci USA 102: 18526–18531, 2005.[Abstract/Free Full Text]
  15. Ewing B, Hillier L, Wendl MC, and Green P. Base-calling of automated sequencer traces using PHRED. I. Accuracy assessment. Genome Res 8: 175–185, 1998.[Abstract/Free Full Text]
  16. Feigelstock D, Thompson P, Mattoo P, Zhang Y, and Kaplan GG. The human homolog of HAVcr-1 codes for a hepatitis A virus cellular receptor. J Virol 72: 6621–6628, 1998.[Abstract/Free Full Text]
  17. Guex N and Peitsch MC. SWISS-MODEL and the Swiss-PdbViewer: an environment for comparative protein modeling. Electrophoresis 18: 2714–2723, 1997.[CrossRef][Web of Science][Medline]
  18. Hart GW. Glycophospholipid anchors. In: Essentials of Glycobiology, edited by Varki A, Cummings R, Esko J, Freeze H, Hart G, and Marth J. Cold Spring Harbor, NY: Cold Spring Harbor Laboratory, 1999, p. 131–144.
  19. Hashizume K, Ishiwata H, Kizaki K, Yamada O, Takahashi T, Imai K, Patel OV, Akagi S, Shimizu M, Takahashi S, Katsuma S, Shiojima S, Hirasawa A, Tsujimoto G, Todoroki J, and Izaike Y. Implantation and placental development in somatic cell clone recipient cows. Cloning Stem Cells 4: 197–209, 2002.[CrossRef][Medline]
  20. Hedge P, Qi R, Abernathy K, Gay C, Dharap S, Gaspard R, Hughes JE, Snesrud E, Lee N, and Quackenbush J. A concise guide to cDNA microarray analysis. Biotechniques 29: 548–562, 2000.[Web of Science][Medline]
  21. Hein WR and Mackay CR. Prominence of {gamma}{delta} T cells in the ruminant immune system. Immunol Today 12: 30–34, 1991.[CrossRef][Web of Science][Medline]
  22. Heyborne K, Fu YX, Nelson A, Farr A, O'Brien R, and Born W. Recognition of trophoblasts by {gamma}{delta} T cells. J Immunol 153: 2918–2926, 1994.[Abstract]
  23. Hofmann K and Stoffel W. Tmbase—a database of membrane spanning protein segments (Abstract). Biol Chem Hoppe-Seyler 374: 166, 1993.
  24. Huang X and Madan A. CAP3: A DNA sequence assembly program. Genome Res 9: 868–877, 1999.[Abstract/Free Full Text]
  25. Ichimura T, Bonventre JV, Bailly V, Wei H, Hession CA, Cate RL, and Sanicola M. Kidney injury molecule-1 (KIM-1), a putative epithelial cell adhesion molecule containing a novel immunoglobulin domain, is up-regulated in renal cells after injury. J Biol Chem 273: 4135–4142, 1998.[Abstract/Free Full Text]
  26. Ishiwata H, Katsuma S, Kizaki K, Patel OV, Nakano H, Takahashi T, Imai K, Hirasawa A, Shiojima S, Ikawa H, Suzuki Y, Tsujimoto G, Iziake Y, Todoroki J, and Hashizume K. Characterization of gene expression profiles in early bovine pregnancy using a custom cDNA microarray. Mol Reprod Dev 65: 9–18, 2003.[CrossRef][Web of Science][Medline]
  27. Jeanmougin F, Thompson JD, Gouy M, Higgins DG, and Gibson TJ. Multiple sequence alignment with Clustal X. Trends Biochem Sci 23: 403–405, 1998.[CrossRef][Web of Science][Medline]
  28. Kan CC, Solomon E, Belt KT, Chain AC, Hiorns LR, and Fey G. Nucleotide sequence of cDNA encoding human {alpha}2-macroglobulin and assignment of the chromosomal locus. Proc Natl Acad Sci USA 82: 2282–2286, 1985.[Abstract/Free Full Text]
  29. Kaplan G, Totsuka A, Thompson P, Akatsuka T, Moritsugu Y, and Feinstone SM. Identification of a surface glycoprotein on African green monkey kidney cells as a receptor for hepatitis A virus. EMBO J 15: 4282–4296, 1996.[Web of Science][Medline]
  30. Keeler C, Dannies PS, and Hodsdon ME. The tertiary structure and backbone dynamics of human prolactin. J Mol Biol 328: 1105–1121, 2003.[CrossRef][Web of Science][Medline]
  31. Kent WJ. The human genome browser at UCSC. Genome Res 12: 996–1006, 2002.[Abstract/Free Full Text]
  32. Kessler MA, Duello TM, and Schuler LA. Expression of prolactin-related hormones in the early bovine conceptus, and potential for paracrine effect on the endometrium. Endocrinology 129: 1885–1895, 1991.[Abstract/Free Full Text]
  33. Kessler MA, Milosavljevic M, Zieler CG, and Schuler LA. A subfamily of bovine prolactin-related transcripts distinct from placental lactogen in the fetal placenta. Biochemistry 28: 5154–5161, 1989.[CrossRef][Medline]
  34. Kessler MA and Schuler LA. Purification and properties of placental prolactin-related protein-I. Placenta 18: 29–36, 1997.[CrossRef][Web of Science][Medline]
  35. Knowles NJ and Samuel AR. Molecular epidemiology of foot-and-mouth disease virus. Virus Res 91: 65–80, 2003.[CrossRef][Medline]
  36. Koonin EV. Orthologs, paralogs, and evolutionary genomics. Annu Rev Genet 39: 309–338, 2005.[CrossRef][Web of Science][Medline]
  37. Kramer KK, Duffy JY, Klemann SW, Bixby JA, Low BG, Pope WF, and Roberts RM. Selective cloning of cDNA for secretory proteins of early embryos. Identification of a transiently expressed Kunitz domain protein from preimplantation sheep trophoblast. J Biol Chem 269: 7255–7261, 1994.[Abstract/Free Full Text]
  38. Larson JH, Rebeiz MJ, Stiening CM, Windish RL, Beever JE, and Lewin HA. MHC class I-like genes in cattle, MHCLA, with similarity to genes encoding NK cell stimulatory ligands. Immunogenetics 55: 16–22, 2003.[Medline]
  39. Liu WJ, Gottshall SL, and Hansen PJ. Increased expression of cell surface markers on endometrial {gamma}{delta} T-cell receptor+ intraepithelial lymphocytes induced by the local presence of the sheep conceptus. Am J Reprod Immunol 37: 199–205, 1997.[Medline]
  40. Lyman SD, Escobar S, Rousseau AM, Armstrong A, and Fanslow WC. Identification of CD7 as a cognate of the human K12 (SECTM1) protein. J Biol Chem 275: 3431–3437, 2000.[Abstract/Free Full Text]
  41. Mackay CR and Hein WR. A large proportion of bovine T cells express the {gamma}{delta} T cell receptor and show a distinct tissue distribution and surface phenotype. Int Immunol 1: 540–545, 1989.[Abstract/Free Full Text]
  42. MacLean JA 2nd, Chakrabarty A, Xie S, Bixby JA, Roberts RM, and Green JA. Family of Kunitz proteins from trophoblast: expression of the trophoblast Kunitz domain proteins (TKDP) in cattle and sheep. Mol Reprod Dev 65: 30–40, 2003.[CrossRef][Web of Science][Medline]
  43. MacLean JA 2nd, Roberts RM, and Green JA. Atypical Kunitz-type serine proteinase inhibitors produced by the ruminant placenta. Biol Reprod 71: 455–463, 2004.[Abstract/Free Full Text]
  44. Marchler-Bauer A, Panchenko AR, Shoemaker BA, Thiessen PA, Geer LY, and Bryant SH. CDD: a database of conserved domain alignments with links to domain three-dimensional structure. Nucleic Acids Res 30: 281–283, 2002.[Abstract/Free Full Text]
  45. Medawar PD. Some immunological and endocrinological problems raised by the evolution of viviparity in vertebrates. Symp Soc Exp Biol 7: 320–338, 1953.[Web of Science]
  46. Milosavljevic M, Duello TM, and Schuler LA. In situ localization of two prolactin-related messenger ribonucleic acids to binucleate cells of bovine placentomes. Endocrinology 125: 883–889, 1989.[Abstract/Free Full Text]
  47. Mossman HW. Vertebrate Fetal Membranes: Comparative Ontogeny and Morphology, Evolution, Phylogenetic Significance, Basic Functions, Research Opportunities. New Brunswick, NJ: Rutgers Univ. Press, 1987.
  48. Mouse Genome Sequencing Consortium. Initial sequencing and comparative analysis of the mouse genome. Nature 420: 520–562, 2002.[CrossRef][Medline]
  49. Nakai K and Kanehisa M. A knowledge base for predicting protein localization sites in eukaryotic cells. Genomics 14: 897–911, 1992.[CrossRef][Web of Science][Medline]
  50. Nei M and Gojobori T. Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol Biol Evol 3: 418–426, 1986.[Abstract]
  51. Ohno S. Evolution by Gene Duplication. New York: Springer-Verlag, 1970.
  52. Rat Genome Sequencing Project Consortium. Genome sequence of the Brown Norway rat yields insights into mammalian evolution. Nature 428: 493–521, 2004.[CrossRef][Medline]
  53. Richardson RT, Sivashanmugam P, Hall SH, Hamil KG, Moore PA, Ruben SM, French FS, and O'Rand M. Cloning and sequencing of human Eppin: a novel family of protease inhibitors expressed in the epididymis and testis. Gene 270: 93–102, 2001.[CrossRef][Web of Science][Medline]
  54. Roberts RM, Ealy AD, Alexenko AP, Han CS, and Ezashi T. Trophoblast interferons. Placenta 20: 259–264, 1999.[Web of Science][Medline]
  55. Roberts RM, Xie S, and Mathialagan N. Maternal recognition of pregnancy. Biol Reprod 54: 294–302, 1996.[Abstract]
  56. Rost B and Liu J. The PredictProtein server. Nucleic Acids Res 31: 3300–3304, 2003.[Abstract/Free Full Text]
  57. Schuler LA and Hurley WL. Molecular cloning of a prolactin-related mRNA expressed in bovine placenta. Proc Natl Acad Sci USA 84: 5650–5654, 1987.[Abstract/Free Full Text]
  58. Schuler LA, Kessler MA, Tanaka M, and Nakashima K. Nomenclature clarification for the bovine placental prolactin-related hormones. Endocrinology 129: 2057, 1991.[Abstract/Free Full Text]
  59. Silberstein E, Dveksler G, and Kaplan GG. Neutralization of hepatitis A virus (HAV) by an immunoadhesin containing the cysteine-rich region of HAV cellular receptor-1. J Virol 75: 717–725, 2001.[Abstract/Free Full Text]
  60. Silberstein E, Xing L, van de Beek W, Lu J, Cheng H, and Kaplan GG. Alteration of hepatitis A virus (HAV) particles by a soluble form of HAV cellular receptor 1 containing the immunoglobin- and mucin-like regions. J Virol 77: 8765–8774, 2003.[Abstract/Free Full Text]
  61. Slentz-Kesler KA, Hale LP, and Kaufman RE. Identification and characterization of K12 (SECTM1), a novel human gene that encodes a Golgi-associated protein with transmembrane and secreted isoforms. Genomics 47: 327–340, 1998.[CrossRef][Medline]
  62. Springer MS, Murphy WJ, Eizirik E, and O'Brien SJ. Placental mammal diversification and the Cretaceous-Tertiary boundary. Proc Natl Acad Sci USA 100: 1056–1061, 2003.[Abstract/Free Full Text]
  63. Tanaka M, Minoura H, Ushiro H, and Nakashima K. A novel cDNA clone encoding a prolactin-like protein that lacks the two C-terminal cysteine residues isolated from bovine placenta. Biochim Biophys Acta 1088: 385–389, 1991.[Medline]
  64. Tanaka M, Yamakawa M, Watahiki M, Yamamoto M, and Nakashima K. Isolation of a novel prolactin-like cDNA clone from bovine placenta: occurrence of new family members. Biochim Biophys Acta 1008: 193–197, 1989.[Medline]
  65. Thompson P, Lu J, and Kaplan GG. The Cys-rich region of hepatitis A virus cellular receptor 1 is required for binding of hepatitis A virus and protective monoclonal antibody 190/4. J Virol 72: 3751–3761, 1998.[Abstract/Free Full Text]
  66. Trivers RL. Parent-offspring conflict. Am Zool 14: 249–264, 1974.[Web of Science]
  67. Ushizawa K, Kaneyama K, Takahashi T, Tokunaga T, Tsunoda Y, and Hashizume K. Cloning and expression of a new member of prolactin-related protein in bovine placenta: bovine prolactin-related protein-VII. Biochem Biophys Res Commun 326: 435–441, 2005.[CrossRef][Web of Science][Medline]
  68. Wang HY, Malek RL, Kwitek AE, Greene AS, Luu TV, Behbahani B, Frank B, Quackenbush J, and Lee NH. Assessing unmodified 70-mer oligonucleotide probe performance on glass-slide microarrays. Genome Biol 4: R5.1–R5.13, 2003.[CrossRef][Medline]
  69. Womack JE, Johnson JS, Owens EK, Rexroad CE 3rd, Schlapfer J, and Yang YP. A whole-genome radiation hybrid panel for bovine gene mapping. Mamm Genome 8: 854–856, 1997.[CrossRef][Web of Science][Medline]
  70. Wooding FB and Wathes DC. Binucleate cell migration in the bovine placentome. J Reprod Fertil 59: 425–430, 1980.[Abstract/Free Full Text]
  71. Xie S, Green J, Bixby JB, Szafranska B, DeMartini JC, Hecht S, and Roberts RM. The diversity and evolutionary relationships of the pregnancy-associated glycoproteins, an aspartic proteinase subfamily consisting of many trophoblast-expressed genes. Proc Natl Acad Sci USA 94: 12809–12816, 1997.[Abstract/Free Full Text]
  72. Yamada O, Todoroki J, Kizaki K, Takahashi T, Imai K, Patel OV, Schuler LA, and Hashizume K. Expression of prolactin-related protein I at the fetomaternal interface during the implantation period in cows. Reproduction 124: 427–437, 2002.[Abstract]
  73. Yamakawa M, Tanaka M, Koyama M, Kagesato Y, Watahiki M, Yamamoto M, and Nakashima K. Expression of new members of the prolactin growth hormone gene family in bovine placenta. Isolation and characterization of two prolactin-like cDNA clones. J Biol Chem 265: 8915–8920, 1990.[Abstract/Free Full Text]
  74. Yang Z. PAML: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci 13: 555–556, 1997.[Free Full Text]
  75. Yang Z and Nielsen R. Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models. Mol Biol Evol 17: 32–43, 2000.[Abstract/Free Full Text]
  76. Yang Z, Nielsen R, Goldman N, and Pedersen AM. Codon-substitution models for heterogeneous selection pressure at amino acid sites. Genetics 155: 431–449, 2000.[Abstract/Free Full Text]
  77. Zieler CG, Kessler MA, and Schuler LA. Characterization of a novel prolactin-related protein from bovine fetal placenta. Endocrinology 126: 2377–2382, 1990.[Abstract/Free Full Text]



This article has been cited by other articles:


Home page
J EndocrinolHome page
J. F Trott, K. C Horigan, J. M Gloviczki, K. M Costa, B. A Freking, C. Farmer, K. Hayashi, T. Spencer, J. E Morabito, and R. C Hovey
Tissue-specific regulation of porcine prolactin receptor expression by estrogen, progesterone, and prolactin
J. Endocrinol., July 1, 2009; 202(1): 153 - 166.
[Abstract] [Full Text] [PDF]


Home page
ScienceHome page
The Bovine Genome Sequencing and Analysis Consorti, C. G. Elsik, R. L. Tellam, and K. C. Worley
The Genome Sequence of Taurine Cattle: A Window to Ruminant Biology and Evolution
Science, April 24, 2009; 324(5926): 522 - 528.
[Abstract] [Full Text] [PDF]


Home page
J EndocrinolHome page
T. Takahashi, O. Yamada, M. J Soares, and K. Hashizume
Bovine prolactin-related protein-I is anchored to the extracellular matrix through interactions with type IV collagen
J. Endocrinol., February 1, 2008; 196(2): 225 - 234.
[Abstract] [Full Text] [PDF]


Home page
ReproductionHome page
S L Rodriguez-Zas, K Schellander, and H A Lewin
Biological interpretations of transcriptomic profiles in mammalian oocytes and embryos
Reproduction, February 1, 2008; 135(2): 129 - 139.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental Tables and Figures
Right arrow All Versions of this Article:
25/3/405    most recent
00307.2005v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (12)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Larson, J. H.
Right arrow Articles by Lewin, H. A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Larson, J. H.
Right arrow Articles by Lewin, H. A.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online
Copyright © 2006 by the American Physiological Society.