Sequence evaluation of four pooled-tissue normalized bovine cDNA libraries and construction of a gene index for cattle

Dissecting the genomic regions of selection on the X chromosome in different cattle breeds

Mammalian X and Y chromosomes independently evolved from various autosomes approximately 300 million years ago (MYA). To fully understand the relationship between genomic composition and phenotypic diversity arising due to the course of evolution, we have scanned regions of selection signatures on the X chromosome in different cattle breeds. In this study, we have prepared the datasets of 184 individuals of different cattle breeds and explored the complete X chromosome by utilizing four within-population and two between-population methods. There were 23, 25, 30, 17, 17, and 12 outlier regions identified in Tajima's D, CLR, iHS, ROH, FST, and XP-EHH. Bioinformatics analysis showed that these regions harbor important candidate genes like AKAP4 for reproduction in Brown Swiss, MBTS2 for production traits in Brown Swiss and Guernsey, CXCR3 and CITED1 for health traits in Jersey and Nelore, and BMX and CD40LG for regulation of X chromosome inactivation in Nelore and Gir. We identified genes shared among multiple methods, such as TRNAC-GCA and IL1RAPL1, which appeared in Tajima's D, ROH, and iHS analyses. The gene TRNAW-CCA was found in ROH, CLR and iHS analyses. The X chromosome exhibits a distinctive interaction between demographic factors and genetic variations, and these findings may provide new insight into the X-linked selection in different cattle breeds.

Association analysis of production traits of Japanese quail (Coturnix japonica) using restriction-site associated DNA sequencing

This study was designed to perform an association analysis and identify SNP markers associated with production traits of Japanese quail using restriction-site-associated DNA sequencing. Weekly body weight data from 805 quail were collected from hatching to 16 weeks of age. A total number of 3990 eggs obtained from 399 female quail were used to assess egg quality traits. Egg-related traits were measured at the beginning of egg production (first stage) and at 12 weeks of age (second stage). Five eggs were analyzed at each stage. Traits, such as egg weight, egg length and short axes, eggshell strength and weight, egg equator thickness, yolk weight, diameter, and colour, albumen weight, age of first egg, total number of laid eggs, and egg production rate, were assessed. A total of 383 SNPs and 1151 associations as well as 734 SNPs and 1442 associations were identified in relation to quail production traits using general linear model (GLM) and mixed linear model (MLM) approaches, respectively. The GLM-identified SNPs were located on chromosomes 1-13, 15, 17-20, 24, 26-28, and Z, underlying phenotypic traits, except for egg and albumen weight at the first stage and yolk yellowness at the second stage. The MLM-identified SNPs were positioned on defined chromosomes associated with phenotypic traits except for the egg long axis at the second stage of egg production. Finally, 35 speculated genes were identified as candidate genes for the targeted traits based on their nearest positions. Our findings provide a deeper understanding and allow a more precise genetic improvement of production traits of Galliformes, particularly in Japanese quail.

At the dawn of the transcriptomic medicine

Progress in genomic analytical technologies has improved our possibilities to obtain information regarding DNA, RNA, and their dynamic changes that occur over time or in response to specific challenges. This information describes the blueprint for cells, tissues, and organisms and has fundamental importance for all living organisms. This review focuses on the technological challenges to analyze the transcriptome and what is the impact of transcriptomics on precision medicine. The transcriptome is a term that covers all RNA present in cells and a substantial part of it will never be translated into protein but is nevertheless functional in determining cell phenotype. Recent developments in transcriptomics have challenged the fundamentals of the central dogma of biology by providing evidence of pervasive transcription of the genome. Such massive transcriptional activity is challenging the definition of a gene and especially the term "pseudogene" that has now been demonstrated in many examples to be both transcribed and translated. We also review the common sources of biomaterials for transcriptomics and justify the suitability of whole blood RNA as the current optimal analyte for clinical transcriptomics. At the end of the review, a brief overview of the clinical implications of transcriptomics in clinical trial design and clinical diagnosis is given. Finally, we introduce the transcriptome as a target for modern drug development as a tool for extending our capacity for precision medicine in multiple diseases.

Symposium review: Advances in sequencing technology herald a new frontier in cattle genomics and genome-enabled selection

The cattle reference genome assembly has underpinned major innovations in beef and dairy genetics through genome-enabled selection, including removal of deleterious recessive variants and selection for favorable alleles affecting quantitative production traits. The initial reference assemblies, up to and including UMD3.1 and Btau4.1, were based on a combination of clone-by-clone sequencing of bacterial artificial chromosome clones generated from blood DNA of a Hereford bull and whole-genome shotgun sequencing of blood DNA from his inbred daughter/granddaughter named L1 Dominette 01449 (Dominette). The approach introduced assembly gaps, misassemblies, and errors, and it limited the ability to assemble regions that undergo rearrangement in blood cells, such as immune gene clusters. Nonetheless, the reference supported the creation of genotyping tools and provided a basis for many studies of gene expression. Recently, long-read sequencing technologies have emerged that facilitated a re-assembly of the reference genome, using lung tissue from Dominette to resolve many of the problems and providing a bridge to place historical studies in common context. The new reference, ARS-UCD1.2, successfully assembled germline immune gene clusters and improved overall continuity (i.e., reduction of gaps and inversions) by over 250-fold. This reference properly places nearly all of the legacy genetic markers used for over a decade in the industry. In this review, we discuss the improvements made to the cattle reference; remaining issues present in the assembly; tools developed to support genome-based studies in beef and dairy cattle; and the emergence of newer genome assembly methods that are producing even higher-quality assemblies for other breeds of cattle at a fraction of the cost. The new frontier for cattle genomics research will likely include a transition from the individual Hereford reference genome, to a "pan-genome" reference, representing all the DNA segments existing in commonly used cattle breeds, bringing the cattle reference into line with the current direction of human genome research.

Copy number variation detection in Chinese indigenous cattle by whole genome sequencing

Copy number variation (CNV) refers to a kind of structural variation, having functional and evolutionary effects on phenotypes. Thus far, further elucidation of the CNVs in different Chinese indigenous cattle breeds by whole genome sequencing have yet not been done. In this study, a comprehensive genomic analysis was performed on 75 cattle individuals including six Chinese indigenous cattle breeds and two non-native specialized beef cattle breeds. Based on the 11,486 CNVRs discovered, population analysis was performed, showed that all the cattle breeds clustered in to three clades, consistent with their lineages Bos taurus, Bos taurus × Bos indicus and Bos indicus. Importantly, a set of CNVRs related genes were found to be associated with the traits of interest, which include meat production or quality (CAST, ACTC1, etc.), adaption (BLA-DQB, EGLN2, etc.) and coat color (KIT, MITF, etc.). These results provide valuable full genome variation resources for Chinese bovine genome research and would be helpful for cattle breeding and selection programs in the future.

De novo assembly and characterization of the Hucho taimen transcriptome

Taimen (Hucho taimen) is an important ecological and economic species that is classified as vulnerable by the IUCN Red List of Threatened Species; however, limited genomic information is available on this species. RNA‐Seq is a useful tool for obtaining genetic information and developing genetic markers for nonmodel species in addition to its application in gene expression profiling. In this study, we performed a comprehensive RNA‐Seq analysis of taimen. We obtained 157 M clean reads (14.7 Gb) and used them to de novo assemble a high‐quality transcriptome with a N50 size of 1,060 bp. In the assembly, 82% of the transcripts were annotated using several databases, and 14,666 of the transcripts contained a full open reading frame. The assembly covered 75% of the transcripts of Atlantic salmon and 57.3% of the protein‐coding genes of rainbow trout. To learn about the genome evolution, we performed a systematic comparative analysis across 11 teleosts including eight salmonids and found 313 unique gene families in taimen. Using Atlantic salmon and rainbow trout transcriptomes as the background, we identified 250 positive selection transcripts. The pathway enrichment analysis revealed a unique characteristic of taimen: It possesses more immune‐related genes than Atlantic salmon and rainbow trout; moreover, some genes have undergone strong positive selection. We also developed a pipeline for identifying microsatellite marker genotypes in samples and successfully identified 24 polymorphic microsatellite markers for taimen. These data and tools are useful for studying conservation genetics, phylogenetics, evolution among salmonids, and selective breeding for threatened taimen.

Genetic Architecture and Selection of Chinese Cattle Revealed by Whole Genome Resequencing

The bovine genetic resources in China are diverse, but their value and potential are yet to be discovered. To determine the genetic diversity and population structure of Chinese cattle, we analyzed the whole genomes of 46 cattle from six phenotypically and geographically representative Chinese cattle breeds, together with 18 Red Angus cattle genomes, 11 Japanese black cattle genomes and taurine and indicine genomes available from previous studies. Our results showed that Chinese cattle originated from hybridization between Bos taurus and Bos indicus. Moreover, we found that the level of genetic variation in Chinese cattle depends upon the degree of indicine content. We also discovered many potential selective sweep regions associated with domestication related to breed-specific characteristics, with selective sweep regions including genes associated with coat color (ERCC2, MC1R, ZBTB17, and MAP2K1), dairy traits (NCAPG, MAPK7, FST, ITFG1, SETMAR, PAG1, CSN3, and RPL37A), and meat production/quality traits (such as BBS2, R3HDM1, IGFBP2, IGFBP5, MYH9, MYH4, and MC5R). These findings substantially expand the catalogue of genetic variants in cattle and reveal new insights into the evolutionary history and domestication traits of Chinese cattle.

Identification of an Alternative Splicing Product of the Otx2 Gene Expressed in the Neural Retina and Retinal Pigmented Epithelial Cells

To investigate the complexity of alternative splicing in the retina, we sequenced and analyzed a total of 115,706 clones from normalized cDNA libraries from mouse neural retina (66,217) and rat retinal pigmented epithelium (49,489). Based upon clustering the cDNAs and mapping them with their respective genomes, the estimated numbers of genes were 9,134 for the mouse neural retina and 12,050 for the rat retinal pigmented epithelium libraries. This unique collection of retinal of messenger RNAs is maintained and accessible through a web-base server to the whole community of retinal biologists for further functional characterization. The analysis revealed 3,248 and 3,202 alternative splice events for mouse neural retina and rat retinal pigmented epithelium, respectively. We focused on transcription factors involved in vision. Among the six candidates suitable for functional analysis, we selected Otx2S, a novel variant of the Otx2 gene with a deletion within the homeodomain sequence. Otx2S is expressed in both the neural retina and retinal pigmented epithelium, and encodes a protein that is targeted to the nucleus. OTX2S exerts transdominant activity on the tyrosinase promoter when tested in the physiological environment of primary RPE cells. By overexpressing OTX2S in primary RPE cells using an adeno associated viral vector, we identified 10 genes whose expression is positively regulated by OTX2S. We find that OTX2S is able to bind to the chromatin at the promoter of the retinal dehydrogenase 10 (RDH10) gene.

Profiling Gene Expression in Germinating Brassica Roots

Based on previously developed solid-phase gene extraction (SPGE) we examined the mRNA profile in primary roots of Brassica rapa seedlings for highly expressed genes like ACT7 (actin7), TUB (tubulin1), UBQ (ubiquitin), and low expressed GLK (glucokinase) during the first day post-germination. The assessment was based on the mRNA load of the SPGE probe of about 2.1 ng. The number of copies of the investigated genes changed spatially along the length of primary roots. The expression level of all genes differed significantly at each sample position. Among the examined genes ACT7 expression was most even along the root. UBQ was highest at the tip and root-shoot junction (RS). TUB and GLK showed a basipetal gradient. The temporal expression of UBQ was highest in the MZ 9 h after primary root emergence and higher than at any other sample position. Expressions of GLK in EZ and RS increased gradually over time. SPGE extraction is the result of oligo-dT and oligo-dA hybridization and the results illustrate that SPGE can be used for gene expression profiling at high spatial and temporal resolution. SPGE needles can be used within two weeks when stored at 4 °C. Our data indicate that gene expression studies that are based on the entire root miss important differences in gene expression that SPGE is able to resolve for example growth adjustments during gravitropism.

Localization, developmental regulation and stress induction of p8-L1 and p8-L2 in zebrafish

The transcription cofactor, p8, has a basic helix-loop-helix motif and is involved in several metabolic processes. It has previously been reported that p8 plays an important role in stress response in zebrafish. We expanded on these studies to further investigate the roles of two p8-like cDNAs found in zerbrafish, p8-L1 and p8-L2. Zebrafish p8-L1 and p8-L2 cDNAs consist of 742 bp and 813 bp, respectively. Both have the same open reading frame despite different genomic structures and encode a 76 amino acid polypeptide with conserved a Phospho_p8 domain. These two cDNA’s were abundant in different tissues in the zebrafish: p8-L1 was high in intestine and muscle tissue and in low abundance in the backbone, whereas the highest amount of p8-L2 was found in the backbone, similar to p8. During embryogenesis, both p8-L1 and p8-L2 were abundant at the cleavage stage and decreased gradually in abundance from blastula to pharyngula stage. However, p8-L1 abundance increased during hatching as observed in p8. Quantitative real-time PCR assay suggested that p8-L1 and p8-L2 were both up-regulated significantly under osmotic pressure and pH value challenges, suggesting an important role of p8-L1 and p8-L2 genes in stress response.

Transcriptome analysis reveals the time of the fourth round of genome duplication in common carp (Cyprinus carpio)

Background

Common carp (Cyprinus carpio) is thought to have undergone one extra round of genome duplication compared to zebrafish. Transcriptome analysis has been used to study the existence and timing of genome duplication in species for which genome sequences are incomplete. Large-scale transcriptome data for the common carp genome should help reveal the timing of the additional duplication event.

Results

We have sequenced the transcriptome of common carp using 454 pyrosequencing. After assembling the 454 contigs and the published common carp sequences together, we obtained 49,669 contigs and identified genes using homology searches and an ab initio method. We identified 4,651 orthologous pairs between common carp and zebrafish and found 129,984 paralogous pairs within the common carp. An estimation of the synonymous substitution rate in the orthologous pairs indicated that common carp and zebrafish diverged 120 million years ago (MYA). We identified one round of genome duplication in common carp and estimated that it had occurred 5.6 to 11.3 MYA. In zebrafish, no genome duplication event after speciation was observed, suggesting that, compared to zebrafish, common carp had undergone an additional genome duplication event. We annotated the common carp contigs with Gene Ontology terms and KEGG pathways. Compared with zebrafish gene annotations, we found that a set of biological processes and pathways were enriched in common carp.

Conclusions

The assembled contigs helped us to estimate the time of the fourth-round of genome duplication in common carp. The resource that we have built as part of this study will help advance functional genomics and genome annotation studies in the future.

A genome-wide SNP scan reveals novel loci for egg production and quality traits in white leghorn and brown-egg dwarf layers

Availability of the complete genome sequence as well as high-density SNP genotyping platforms allows genome-wide association studies (GWAS) in chickens. A high-density SNP array containing 57,636 markers was employed herein to identify associated variants underlying egg production and quality traits within two lines of chickens, i.e., White Leghorn and brown-egg dwarf layers. For each individual, age at first egg (AFE), first egg weight (FEW), and number of eggs (EN) from 21 to 56 weeks of age were recorded, and egg quality traits including egg weight (EW), eggshell weight (ESW), yolk weight (YW), eggshell thickness (EST), eggshell strength (ESS), albumen height(AH) and Haugh unit(HU) were measured at 40 and 60 weeks of age. A total of 385 White Leghorn females and 361 brown-egg dwarf dams were selected to be genotyped. The genome-wide scan revealed 8 SNPs showing genome-wise significant (P<1.51E-06, Bonferroni correction) association with egg production and quality traits under the Fisher's combined probability method. Some significant SNPs are located in known genes including GRB14 and GALNT1 that can impact development and function of ovary, but more are located in genes with unclear functions in layers, and need to be studied further. Many chromosome-wise significant SNPs were also detected in this study and some of them are located in previously reported QTL regions. Most of loci detected in this study are novel and the follow-up replication studies may be needed to further confirm the functional significance for these newly identified SNPs.

EST analysis on the gonad development related organs and microarray screen for differentially expressed genes in mature ovary and testis of Scylla paramamosain

A total of 5160 high quality ESTs (expressed sequence tags) averaging 357 bp were collected from normalized cDNA libraries created from testis, ovary and mixed organs of mud crab Scylla paramamosain. Clustering and assembly of these ESTs resulted in a total of 3837 unique sequences with 576 overlapping contigs and 3261 singletons. Comparisons with the GenBank non-redundant (Nr) protein database (BLASTx, e-values <10(-5)) revealed putative functions or matched homologs from other organisms for 847 (22%) of the ESTs. Several gonad development related genes such as cathepsin C, thioredoxin peroxidase, vitellogenin receptor precursor, 50S ribosomal protein L24 and ubiquitin-conjugating enzyme E2 isoform 2 were identified from this EST project and demonstrated as gonad differential expression genes by rqRT-PCR. Sixty five different types of SSRs (simple sequence repeats) were identified from the total 411 EST-SSR motifs. A home-made cDNA microarray containing 5664 spots was developed and the hybridization results indicated that 39 unique transcripts were differentially expressed in testis and ovaries (P<0.05). The expression levels of eleven unique transcripts examined by rqRT-PCR were matched with microarray fairly. These results will provide a useful resource for functional genomic studies on the biology of reproduction of mud crab.

A 4 Mb high resolution BAC contig on bovine chromosome 1q12 and comparative analysis with human chromosome 21q22

The bovine RPCI-42 BAC library was screened to construct a sequence-ready ~4 Mb single contig of 92 BAC clones on BTA 1q12. The contig covers the region between the genes KRTAP8P1 and CLIC6. This genomic segment in cattle is of special interest as it contains the dominant gene responsible for the hornless or polled phenotype in cattle. The construction of the BAC contig was initiated by screening the bovine BAC library with heterologous cDNA probes derived from 12 human genes of the syntenic region on HSA 21q22. Contig building was facilitated by BAC end sequencing and chromosome walking. During the construction of the contig, 165 BAC end sequences and 109 single-copy STS markers were generated. For comparative mapping of 25 HSA 21q22 genes, genomic PCR primers were designed from bovine EST sequences and the gene-associated STSs mapped on the contig. Furthermore, bovine BAC end sequence comparisons against the human genome sequence revealed significant matches to HSA 21q22 and allowed the in silico mapping of two new genes in cattle. In total, 31 orthologues of human genes located on HSA 21q22 were directly mapped within the bovine BAC contig, of which 16 genes have been cloned and mapped for the first time in cattle. In contrast to the existing comparative bovine–human RH maps of this region, these results provide a better alignment and reveal a completely conserved gene order in this 4 Mb segment between cattle, human and mouse. The mapping of known polled linked BTA 1q12 microsatellite markers allowed the integration of the physical contig map with existing linkage maps of this region and also determined the exact order of these markers for the first time. Our physical map and transcript map may be useful for positional cloning of the putative polled gene in cattle. The nucleotide sequence data reported in this paper have been submitted to EMBL and have been assigned Accession Numbers AJ698510–AJ698674.

Using the DFCI gene index databases for biological discovery

The DFCI Gene Index Web pages provide access to analyses of ESTs and gene sequences for nearly 114 species, as well as a number of resources derived from these. Each species-specific database is presented using a common format with a home page. A variety of methods exist that allow users to search each species-specific database. Methods implemented currently include nucleotide or protein sequence queries using WU-BLAST, text-based searches using various sequence identifiers, searches by gene, tissue and library name, and searches using functional classes through Gene Ontology assignments. This protocol provides guidance for using the Gene Index Databases to extract information.

Annotation of novel transcripts putatively relevant for bovine fat metabolism

Two bovine transcripts encoded by the interleukin-1 receptor-associated kinase 1 (IRAK1) gene and the locus LOC618944 predicted as similar to human chromosome 6 open reading frame 52 (C6orf52) gene had indicated divergent expression in bovine skeletal muscle containing different amount of intramuscular fat in a pilot screening experiment. However, for both loci any role in the regulation of energy or fat metabolism is not yet described. In this study, we validated and refined gene structure, screened for mRNA splice variants and analyzed the tissue-specific gene expression patterns of both loci as a prerequisite to elucidate their potential physiological function. Based on comparative sequence analysis, a new full-length gene model for the bovine IRAK1 gene was developed and confirmed experimentally. Expression of IRAK1 mRNA was found in a variety of tissues, and a splice variant was identified in skeletal muscle caused by an in-frame deleted segment of 210 bp affecting regions of intrinsic disorder in the respective protein. For the locus LOC618944, our data contributed to a revised gene model and its assignment to BTA23 (bovine chromosome 23) on the current bovine genome assembly supported by comparative similarity analysis between the bovine and human genomes and experimental data. Furthermore, we identified several splice variants in mammary gland, fat and skeletal muscle tissue and detected a highly similar processed pseudogene on BTA26. All transcript variants of LOC618944 detected in the analyzed tissues represent noncoding RNAs. For both loci, our results suggest yet undetected physiological functions in tissues relevant for fat or energy metabolism in cattle.

Molecular mechanisms regulating bovine ovarian follicular selection

Transcription profiling of ovarian follicles. Understanding the mechanisms by which a single follicle is selected for further ovulation is important to control fertility in mammals. However, development of new treatments is limited by our poor understanding of molecular mechanisms regulating follicular selection. Our hypothesis is that genes involved in the control of cell proliferation and apoptosis are differentially regulated during follicular selection. Our objective was to identify these new genes. Bovine follicles were collected and gene expression levels were measured using microarrays. First, follicles were allocated to three groups, according to the time spent from the initiation of follicular wave to surgery (24 H, 36 H, and 48-60 H). Fifty-seven genes are differentially expressed at a false discovery rate of 5%. These genes are involved in the control of lipid metabolism (P-value = 0.0005), cell proliferation (0.007), cell death (0.003), cell morphology (0.003), and immune response (0.003). Follicles were also grouped into four categories, according to the expected time of deviation (early deviation; 8 mm, mid-deviation; 8.5 mm, late deviation; 9 mm, dominant follicles; >or=10 mm). One hundred and twenty eight genes are differentially expressed between these four groups, including genes involved in cell proliferation (0.00002), cell death (0.0006), cell-to-cell signaling (0.003), cell morphology (0.003), lipid metabolism (0.0004), and immune response (0.00007). The expression levels of 10 genes were confirmed using quantitative real time PCR. As expected, we identified new differentially regulated genes involved in the control of cell growth and apoptosis. We also discovered a potential role for immune cells, and in particular macrophages, in follicular selection.

Use of a candidate gene array to delineate gene expression patterns in cattle selected for resistance or susceptibility to intestinal nematodes

In the present study, we use microarray technology to investigate the expression patterns of 381 genes with known association to host immune responses. Hybridization targets were derived from previously characterized bovine cDNAs. A total of 576 reporters (473 sequence-validated cDNAs and 77 controls) were spotted onto glass slides in two sets of four replicates. Two color, comparative hybridizations across both mesenteric lymph node (MLN) and small intestine mucosa (SIM) RNA samples were done between animals with previously demonstrated phenotypic differences based on natural exposure to gastrointestinal (GI) nematodes over a 6-month exposure period. A total of 138 significant hybridization differences were detected by mixed model analysis of variance. A subset of these significant differences was validated by quantitative, real-time RT-PCR to assay transcript levels for 18 genes. These results confirmed that in the SIM, susceptible animals showed significantly higher levels in the genes encoding IGHG1, CD3E, ACTB, IRF1, CCL5 and C3, while in the MLN of resistant animals, higher levels of expression were confirmed for PTPRC, CD1D and ITGA4. Combined, the results indicate that immune responses against GI nematode infections involve multiple response pathways. Higher levels of expression for IgE receptor, integrins, complement, monocyte/macrophage and tissue factors are related to resistance. In contrast, higher levels of expression for immunoglobulin chains and TCRs are related to susceptibility. Identification of these genes provides a framework to better understand the genetic variation underlying parasite resistance.

Genomic resources for a commercial flatfish, the Senegalese sole (Solea senegalensis): EST sequencing, oligo microarray design, and development of the Soleamold bioinformatic platform

Background

The Senegalese sole, Solea senegalensis, is a highly prized flatfish of growing commercial interest for aquaculture in Southern Europe. However, despite the industrial production of Senegalese sole being hampered primarily by lack of information on the physiological mechanisms involved in reproduction, growth and immunity, very limited genomic information is available on this species.

Results

Sequencing of a S. senegalensis multi-tissue normalized cDNA library, from adult tissues (brain, stomach, intestine, liver, ovary, and testis), larval stages (pre-metamorphosis, metamorphosis), juvenile stages (post-metamorphosis, abnormal fish), and undifferentiated gonads, generated 10,185 expressed sequence tags (ESTs). Clones were sequenced from the 3'-end to identify isoform specific sequences. Assembly of the entire EST collection into contigs gave 5,208 unique sequences of which 1,769 (34%) had matches in GenBank, thus showing a low level of redundancy. The sequence of the 5,208 unigenes was used to design and validate an oligonucleotide microarray representing 5,087 unique Senegalese sole transcripts. Finally, a novel interactive bioinformatic platform, Soleamold, was developed for the Senegalese sole EST collection as well as microarray and ISH data.

Conclusion

New genomic resources have been developed for S. senegalensis, an economically important fish in aquaculture, which include a collection of expressed genes, an oligonucleotide microarray, and a publicly available bioinformatic platform that can be used to study gene expression in this species. These resources will help elucidate transcriptional regulation in wild and captive Senegalese sole for optimization of its production under intensive culture conditions.

Using the TIGR gene index databases for biological discovery

The TIGR Gene Index web pages provide access to analyses of ESTs and gene sequences for nearly 60 species, as well as a number of resources derived from these. Each species-specific database is presented using a common format with a homepage. A variety of methods exist that allow users to search each species-specific database. Methods implemented currently include nucleotide or protein sequence queries using WU-BLAST, text-based searches using various sequence identifiers, searches by gene, tissue and library name, and searches using functional classes through Gene Ontology assignments. This protocol provides guidance for using the Gene Index Databases to extract information.

Large-scale identification of human genes implicated in epidermal barrier function

Identification of genes expressed in epidermal granular keratinocytes by ORESTES, including a number that are highly specific for these cells.

Background

During epidermal differentiation, keratinocytes progressing through the suprabasal layers undergo complex and tightly regulated biochemical modifications leading to cornification and desquamation. The last living cells, the granular keratinocytes (GKs), produce almost all of the proteins and lipids required for the protective barrier function before their programmed cell death gives rise to corneocytes. We present here the first analysis of the transcriptome of human GKs, purified from healthy epidermis by an original approach.

Results

Using the ORESTES method, 22,585 expressed sequence tags (ESTs) were produced that matched 3,387 genes. Despite normalization provided by this method (mean 4.6 ORESTES per gene), some highly transcribed genes, including that encoding dermokine, were overrepresented. About 330 expressed genes displayed less than 100 ESTs in UniGene clusters and are most likely to be specific for GKs and potentially involved in barrier function. This hypothesis was tested by comparing the relative expression of 73 genes in the basal and granular layers of epidermis by quantitative RT-PCR. Among these, 33 were identified as new, highly specific markers of GKs, including those encoding a protease, protease inhibitors and proteins involved in lipid metabolism and transport. We identified filaggrin 2 (also called ifapsoriasin), a poorly characterized member of the epidermal differentiation complex, as well as three new lipase genes clustered with paralogous genes on chromosome 10q23.31. A new gene of unknown function, C1orf81, is specifically disrupted in the human genome by a frameshift mutation.

Conclusion

These data increase the present knowledge of genes responsible for the formation of the skin barrier and suggest new candidates for genodermatoses of unknown origin.

Identification of genetic markers for fat deposition and meat tenderness on bovine chromosome 5: development of a low-density single nucleotide polymorphism map

As genetic markers, SNP are well suited for the development of genetic tests for production traits in livestock. They are stable through many generations and can provide direct assessment of individual animal's genetic merit if they are in linkage disequilibrium and phase with functional genetic variation. Bovine chromosome 5 has been shown to harbor genetic variation affecting production traits in multiple cattle populations; thus, this chromosome was targeted for SNP-based marker development and subsequent association analysis with carcass and growth phenotypes. Discovery of SNP was performed in a panel of 16 sires representing two sires from each of seven beef breeds and two Holstein sires by PCR amplification and sequencing using primers designed from genomic sequence obtained by low-coverage sequencing of bacterial artificial chromosome (BAC) clones. From 550 SNP, 296 (54%) were tentatively identified as having a minor allele frequency >10%. Forty-five SNP derived from 15 BAC were chosen based on minor allele frequency and were genotyped in 564 steers and their sires. Production and carcass data were collected on the steers as a part of the Germplasm Evaluation (GPE), Cycle VII Project at the U.S. Meat Animal Research Center (Clay Center, NE), which involves of the evaluation of sires from seven of the most popular U.S. breeds. Haplotypes based on seven SNP derived from a BAC containing the bovine genes HEM1 and PDE1B were associated with traits related to carcass fat. Steers homozygous for the major haplotype had 0.15 +/- 0.04 cm less subcutaneous fat, 0.57 +/- 0.18 kg less rib fat, 0.18 +/- 0.07 lower yield grade, 1.11 +/- 0.35% less predicted fat yield, and 0.79 +/- 0.3% greater predicted retail product yield than heterozygotes. The frequency of the major haplotype was 0.70 in the steers, and it ranged from 0.44 (Limousin) to 0.98 (Simmental and Gelbvieh) in a panel consisting of an average of 20 purebred sires from each of the seven breeds. A second set of haplotypes based on four SNP derived from a BAC containing the genes NOL1 and CHD4 was associated with Warner-Bratzler shear force. Steers homozygous for the major haplotype had 0.27 +/- 0.11 kg greater shear force than those heterozygous for the major haplotype and one of two minor haplotypes. The frequency of the major haplotype was 0.59 in the steers and ranged from 0.27 (Hereford) to approximately 0.95 (Angus and Red Angus) in the panel of purebred sires. These results demonstrate the feasibility of targeting QTL regions for SNP-based marker development and that a low level of coverage can identify markers associated with phenotypic traits.

Resistance and susceptibility to a protozoan parasite of cattle—Gene expression differences in macrophages from different breeds of cattle

Cattle infected with the tick-borne protozoan, Theileria annulata, usually undergo severe morbidity, and mortality ensues in a high proportion of animals. However, we have shown that a Bos indicus breed, the Sahiwal, which originates in a T. annulata endemic area, is more resistant to the parasite. Although Sahiwals become infected, the breed exhibits fewer clinical signs and recovers from a dose of parasite which is fatal in the Holstein B. taurus breed. The Sahiwals have a significantly lower fever response, and lower levels of parasite than the Holsteins. One unusual feature of this disease is the production of acute phase proteins (APP), indicating that the parasite induces high systemic levels of pro-inflammatory cytokines. In the Holsteins there is prolonged production of the APP, alpha1-glycoprotein, which, in contrast, is only slightly elevated in the Sahiwals. As the parasite infects macrophages (mphi), our hypothesis is that the Sahiwals can control the excessive production of pro-inflammatory cytokines in response to infection, and that this control is expressed at the level of the mphi. We thus reasoned that the genes underlying the observed difference in resistance to tropical theileriosis, might be identified by investigating gene expression differences in mphi from both breeds. It is possible that relevant polymorphisms might in themselves result in gene expression differences, so initially we targeted likely candidates. However, we detected no differences in expression of the pro-inflammatory cytokines, tumour necrosis factor-alpha (TNFalpha), interleukin-1beta (IL-1beta) or IL-6, in infected mphi. As it is more likely that polymorphisms in candidate genes influence the expression of other genes involved in interrelated pathways, we undertook a more global approach. We designed a bovine mphi specific cDNA microarray, which contains representatives of 5000 different genes expressed in mphi, and investigated the transcriptional responses of mphi from both breeds in response to a variety of stimuli, including infection with T. annulata. Our results indicate that there are fundamental differences in gene expression in mphi from both breeds in the way they respond to infection, and even in their pre-infection resting state.

Gene expression profiling of muscle tissue in Brahman steers during nutritional restriction

Expression profiling using microarrays allows for the detailed characterization of the gene networks that regulate an animal's response to environmental stresses. During nutritional restriction, processes such as protein turnover, connective tissue remodeling, and muscle atrophy take place in the skeletal muscle of the animal. These processes and their regulation are of interest in the context of managing livestock for optimal production efficiency and product quality. Here we expand on recent research applying complementary DNA (cDNA) microarray technology to the study of the effect of nutritional restriction on bovine skeletal muscle. Using a custom cDNA microarray of 9,274 probes from cattle muscle and s.c. fat libraries, we examined the differential gene expression profile of the LM from 10 Brahman steers under three different dietary treatments. The statistical approach was based on mixed-model ANOVA and model-based clustering of the BLUP solutions for the gene x diet interaction effect. From the results, we defined a transcript profile of 156 differentially expressed array elements between the weight loss and weight gain diet substrates. After sequence and annotation analyses, the 57 upregulated elements represented 29 unique genes, and the 99 downregulated elements represented 28 unique genes. Most of these co-regulated genes cluster into groups with distinct biological function related to protein turnover and cytoskeletal metabolism and contribute to our mechanistic understanding of the processes associated with remodeling of muscle tissue in response to nutritional stress.

A microarray analysis for genes regulated by interferon-tau in ovine luminal epithelial cells

Interferon-tau (IFNT) is released by preimplantation conceptuses of ruminant species and prepares the mother for pregnancy. Although one important function is to protect the corpus luteum from the luteolytic activity of prostaglandin-F 2alpha, IFNT most likely regulates a range of other physiological processes in endometrium. Here, an immortalized cell line from ovine uterine luminal epithelial cells was treated with IFNT for either 8 or 24 h. RNA was subjected to cDNA microarray analysis, with RNA from untreated cells as the reference standard. Of 15 634 genes, 1274 (8%) were IFNT responsive at P<0.01 and 585 at P<0.001 to at least one treatment. Of the latter, 356 were up-regulated and 229 down-regulated. Increasing IFNT concentrations from 10 ng/ml to 10 microg/ml had minor effects, and most genes up- or down-regulated at 8 h were regulated similarly at 24 h. Although IFNT influences many genes implicated in antiviral activity and apoptosis, its action also likely regulates prostaglandin metabolism, growth factors and their receptors, apoptosis and the nuclear factor (NF)-kappaB cascade, extracellular matrix accretion, angiogenesis, blood coagulation, and inflammation. In particular, it increased mRNA concentrations of genes related to the vascular endothelial growth factor R2 pathway of angiogenesis and down-regulated ones associated with hypoxia. Two genes implicated in the antiluteolytic actions of IFNT (encoding cyclooxygenase-2 and the oxytocin receptor respectively) were down-regulated in response to all treatments. IFNT targets a complex range of physiological processes during the establishment of pregnancy.

Using genomic approaches to unravel livestock (host)-tick-pathogen interactions

Ticks and tick-borne diseases are a major constraint on livestock farming in many developing countries, which has a huge impact on their economies. Genomic information is becoming more abundant for many of the species involved, which if exploited successfully could be used to develop new control strategies. Here, we review the genomic resources that are now available and discuss how this information is currently being harnessed or can be used in the future to explore the complex interplay that occurs between livestock hosts, tick vectors and tick-borne pathogens.

Genetic variation and responses to vaccines

Disease is a major source of economic loss to the livestock industry. Understanding the role of genetic factors in immune responsiveness and disease resistance should provide new approaches to the control of disease through development of safe synthetic subunit vaccines and breeding for disease resistance. The major histocompatibility complex (MHC) has been an important candidate locus for immune responsiveness studies. However, it is clear that other loci play an important role. Identifying these and quantifying the relative importance of MHC and non-MHC genes should result in new insights into host–pathogen interactions, and information that can be exploited by vaccine designers. The rapidly increasing information available about the bovine genome and the identification of polymorphisms in immune-related genes will offer potential candidates that control immune responses to vaccines. The bovine MHC,BoLA, encodes two distinct isotypes of class II molecules, DR and DQ, and in about half the common haplotypes theDQgenes are duplicated and expressed. DQ molecules are composed of two polymorphic chains whereas DR consists of one polymorphic and one non-polymorphic chain. Although, it is clear that MHC polymorphism is related to immune responsiveness, it is less clear how different allelic and locus products influence the outcome of an immune response in terms of generating protective immunity in outbred animals. A peptide derived from foot-and-mouth disease virus (FMDV) was used as a probe for BoLA class II function. Both DR and DQ are involved in antigen presentation. In an analysis of T-cell clones specific for the peptide, distinct biases to particular restriction elements were observed. In addition inter-haplotype pairings of DQA and DQB molecules produced functional molecules, which greatly increases the numbers of possible restriction elements, compared with the number of genes, particularly in cattle with duplicatedDQgenes. In a vaccine trial with several peptides derived from FMDV,BoLAclass IIDRB3polymorphisms were correlated with both protection and non-protection. Although variation in immune responsiveness to the FMDV peptide between different individuals is partly explainable byBoLAclass II alleles, other genetic factors play an important role. In a quantitative trait locus project, employing a second-generation cross between Charolais and Holstein cattle, significant sire and breed effects were also observed in T-cell, cytokine and antibody responses to the FMDV peptide. These results suggest that both MHC and non-MHC genes play a role in regulating bovine immune traits of relevance to vaccine design. Identifying these genes and quantifying their relative contributions is the subject of further studies.

Analysis of the immune-inducible genes of Plutella xylostella using expressed sequence tags and cDNA microarray

In the present study, the complex gene expression responses of Plutella xylostella to microbial challenges and injury were surveyed using a newly constructed expressed sequence tag (EST) clone collection and cDNA microarray analysis. A total of 1132 P. xylostella ESTs were cloned, annotated and categorized by their putative functions; these included proteases, protease inhibitors, recognition molecules and anti-microbial peptides. GeneOntology revealed that 4% of the P. xylostella ESTs corresponded to immunity-related genes potentially involved in innate immunity. We then used microarray analysis to identify 44 genes that were differentially expressed with at least a two-fold expression difference in P. xylostella before and after pathogen challenge. Together, our EST categorization and microarray profiling analyses allowed us to identify 70 genes that should be considered candidate immune response genes, providing important new insights into the molecular events that occur during the innate immune response in P. xylostella.

Evolutionarily conserved and non-conserved retrovirus restriction activities of artiodactyl APOBEC3F proteins

The APOBEC3 proteins are unique to mammals. Many inhibit retrovirus infection through a cDNA cytosine deamination mechanism. HIV-1 neutralizes this host defense through Vif, which triggers APOBEC3 ubiquitination and degradation. Here, we report an APOBEC3F-like, double deaminase domain protein from three artiodactyls: cattle, pigs and sheep. Like their human counterparts, APOBEC3F and APOBEC3G, the artiodactyl APOBEC3F proteins are DNA cytosine deaminases that locate predominantly to the cytosol and can inhibit the replication of HIV-1 and MLV. Retrovirus restriction is attributable to deaminase-dependent and -independent mechanisms, as deaminase-defective mutants retain significant anti-retroviral activity. However, unlike human APOBEC3F and APOBEC3G, the artiodactyl APOBEC3F proteins have an active N-terminal DNA cytosine deaminase domain, which elicits a broader dinucleotide deamination preference, and they are resistant to HIV-1 Vif. These data indicate that DNA cytosine deamination; sub-cellular localization and retrovirus restriction activities are conserved in mammals, whereas active site location, local mutational preferences and Vif susceptibility are not. Together, these studies indicate that some properties of the mammal-specific, APOBEC3-dependent retroelement restriction system are necessary and conserved, but others are simultaneously modular and highly adaptable.

A gene coexpression network for bovine skeletal muscle inferred from microarray data

We present the application of large-scale multivariate mixed-model equations to the joint analysis of nine gene expression experiments in beef cattle muscle and fat tissues with a total of 147 hybridizations, and we explore 47 experimental conditions or treatments. Using a correlation-based method, we constructed a gene network for 822 genes. Modules of muscle structural proteins and enzymes, extracellular matrix, fat metabolism, and protein synthesis were clearly evident. Detailed analysis of the network identified groupings of proteins on the basis of physical association. For example, expression of three components of the z-disk, MYOZ1, TCAP, and PDLIM3, was significantly correlated. In contrast, expression of these z-disk proteins was not highly correlated with the expression of a cluster of thick (myosins) and thin (actin and tropomyosins) filament proteins or of titin, the third major filament system. However, expression of titin was itself not significantly correlated with the cluster of thick and thin filament proteins and enzymes. Correlation in expression of many fast-twitch muscle structural proteins and enzymes was observed, but slow-twitch-specific proteins were not correlated with the fast-twitch proteins or with each other. In addition, a number of significant associations between genes and transcription factors were also identified. Our results not only recapitulate the known biology of muscle but have also started to reveal some of the underlying associations between and within the structural components of skeletal muscle.

Development and validation of a bovine macrophage specific cDNA microarray

Background

The response of macrophages to danger signals is an important early stage in the immune response. Our understanding of this complex event has been furthered by microarray analysis, which allows the simultaneous investigation of the expression of large numbers of genes. However, the microarray resources available to study these events in livestock animals are limited.

Results

Here we report the development of a bovine macrophage specific (BoMP) cDNA microarray. The BoMP microarray contains 5026 sequence elements (printed in duplicate) and numerous controls. The majority of the clones incorporated on the microarray were derived from the BoMP cDNA library generated from bovine myeloid cells subjected to various stimuli, including over 900 sequences unique to the library. Additional clones representing immunologically important genes have been included on the BoMP microarray. The microarray was validated by investigating the response of bovine monocytes to stimulation with interferon-γ and lipopolysaccharide using amplified RNA. At 2 and 16 hours post stimulation 695 genes exhibited statistically significant differential expression, including; 26 sequences unique to the BoMP library, interleukin 6, prion protein and toll-like receptor 4.

Conclusion

A 5 K cDNA microarray has been successfully developed to investigate gene expression in bovine myeloid cells. The BoMP microarray is available from the ARK-Genomics Centre for Functional Genomics in Farm Animals, UK.

Generation of a large scale repertoire of Expressed Sequence Tags (ESTs) from normalised rainbow trout cDNA libraries

Background

Within the framework of a genomics project on livestock species (AGENAE), we initiated a high-throughput DNA sequencing program of Expressed Sequence Tags (ESTs) in rainbow trout, Oncorhynchus mykiss.

Results

We constructed three cDNA libraries including one highly complex pooled-tissue library. These libraries were normalized and subtracted to reduce clone redundancy. ESTs sequences were produced, and 96 472 ESTs corresponding to high quality sequence reads were released on the international database, currently representing 42.5% of the overall sequence knowledge in this species. All these EST sequences and other publicly available ESTs in rainbow trout have been included on a publicly available Website (SIGENAE) and have been clustered into a total of 52 930 clusters of putative transcripts groups, including 24 616 singletons. 57.1% of these 52 930 clusters are represented by at least one Agenae EST and 14 343 clusters (27.1%) are only composed by Agenae ESTs. Sequence analysis also reveals that normalization and especially subtraction were effective in decreasing redundancy, and that the pooled-tissue library was representative of the initial tissue complexity.

Conclusion

Due to present work on the construction of rainbow trout normalized cDNA libraries and their extensive sequencing, along with other large scale sequencing programs, rainbow trout is now one of the major fish models in term of EST sequences available in a public database, just after Zebrafish, Danio rerio. This information is now used for the selection of a non redundant set of clones for producing DNA micro-arrays in order to examine global gene expression.

Molecular and functional characterization of bovine beta-defensin-1

We report the biochemical and functional properties of a novel bovine beta-defensin (bBD-1). Cloned from bovine mammary papillary duct epithelia, the bBD-1 cDNA predicts a 69 amino acid propeptide that is much more similar to human beta-defensin-1 (hBD-1) than to other bovine defensins. The bBD-1 gene contains two exons and one 8.5 kb intron. Using RT-PCR, we detected the bBD-1 transcript in the teat mucosa, kidney, vagina, ovary, oviduct, and colon. A synthetic bBD-1 peptide demonstrates potent antibacterial activity against Escherichia coli. The widespread expression of bBD-1 mRNA indicates that bBD-1 may play an important role in the bovine host defense against infections.

The genomic structure and a novel alternatively spliced form of porcine pTalpha chain

A complete genomic nucleotide sequence for porcine pTalpha gene was obtained from a BAC clone, which revealed a novel exon 2 missing in human and murine counterparts. Cattle and dog genomic sequences showed the counterparts corresponding to porcine exon 2. Using thymocyte RNA and RT-PCR, three types of porcine pTalpha-chain cDNA sequences, pTalpha1, pTalpha2 and pTalpha3, were obtained. These three different cDNA sequences were alternatively spliced products with pTalpha1 consisting of exons 1, 2, 3, 4, and 5, pTalpha2 consisting of exons 1, 2, 4, and 5, and pTalpha3 consisting of exons 1, 2, 3 and the intron down stream of exon 3. pTalpha1 and pTalpha2 correspond to previously reported pTalphaa, and pTalphab, respectively, and pTalpha3 is reported for the first time. Using RT-PCR, pTalpha3 appeared expressed predominantly in the thymocyte RNA. The chromosome location of pTalpha was investigated using Radiation Hybrid Map and FISH, both of which revealed the location at SSC7q11-q12.

Chicken genomics resource: sequencing and annotation of 35,407 ESTs from single and multiple tissue cDNA libraries and CAP3 assembly of a chicken gene index

Its accessibility, unique evolutionary position, and recently assembled genome sequence have advanced the chicken to the forefront of comparative genomics and developmental biology research as a model organism. Several chicken expressed sequence tag (EST) projects have placed the chicken in 10th place for accrued ESTs among all organisms in GenBank. We have completed the single-pass 5′-end sequencing of 37,557 chicken cDNA clones from several single and multiple tissue cDNA libraries and have entered 35,407 EST sequences into GenBank. Our chicken EST sequences and those found in public databases (on July 1, 2004) provided a total of 517,727 public chicken ESTs and mRNAs. These sequences were used in the CAP3 assembly of a chicken gene index composed of 40,850 contigs and 79,192 unassembled singlets. The CAP3 contigs show a 96.7% match to the chicken genome sequence. The University of Delaware (UD) EST collection (43,928 clones) was assembled into 19,237 nonredundant sequences (13,495 contigs and 5,742 unassembled singlets). The UD collection contains 6,223 unique sequences that are not found in other public EST collections but show a 76% match to the chicken genome sequence. Our chicken contig and singlet sequences were annotated according to the highest BlastX and/or BlastN hits. The UD CAP3 contig assemblies and singlets are searchable by nucleotide sequence or key word ( http://cogburn.dbi.udel.edu ), and the cDNA clones are readily available for distribution from the chick EST website and clone repository ( http://www.chickest.udel.edu ). The present paper describes the construction and normalization of single and multiple tissue chicken cDNA libraries, high-throughput EST sequencing from these libraries, the CAP3 assembly of a chicken gene index from all public ESTs, and the identification of several nonredundant chicken gene sets for production of custom DNA microarrays.

Novel Urothelium Specific Gene Expression Identified by Differential Display Reverse Transcriptase-Polymerase Chain Reaction

PURPOSE

Understanding the molecular basis of differential gene expression among different tissues at various developmental stages and in neoplastic transformation is an important biological goal. The potential clinical applications of this improved understanding are more precise diagnosis of disease, prediction of prognosis, novel targeted therapies and prediction of response to therapy.

MATERIALS AND METHODS

Differential display reverse transcriptase-polymerase chain reaction was used to compare gene expression in bovine urothelium to that in autologous lung, esophagus, liver and spleen. Products that appeared to have urothelial specific expression were sequenced and assessed for homology with known sequences. Ribonuclease protection assays were used to further confirm the expression pattern.

RESULTS

A total of 32 discrete cDNAs were identified, including 3 products from genes known to be urothelium specific in their expression, 16 with significant homology to bovine, human or mouse expressed sequence tags and 5 with no sequence homology to any currently available sequence. Urothelium specific mRNA expression was confirmed for 3 genes by ribonuclease protection assays and one (Udd06) was further characterized as a urea transporter.

CONCLUSIONS

The use of differential display reverse transcriptase-polymerase chain reaction and other complementary techniques for parallel gene expression analysis will permit the complete characterization of the urothelial transcriptome and help identify potential molecular targets for rationally targeted therapy.

Characterization of 954 bovine full-CDS cDNA sequences

Background

Genome assemblies rely on the existence of transcript sequence to stitch together contigs, verify assembly of whole genome shotgun reads, and annotate genes. Functional genomics studies also rely on transcript sequence to create expression microarrays or interpret digital tag data produced by methods such as Serial Analysis of Gene Expression (SAGE). Transcript sequence can be predicted based on reconstruction from overlapping expressed sequence tags (EST) that are obtained by single-pass sequencing of random cDNA clones, but these reconstructions are prone to errors caused by alternative splice forms, transcripts from gene families with related sequences, and expressed pseudogenes. These errors confound genome assembly and annotation. The most useful transcript sequences are derived by complete insert sequencing of clones containing the entire length, or at least the full protein coding sequence (CDS) portion, of the source mRNA. While the bovine genome sequencing initiative is nearing completion, there is currently a paucity of bovine full-CDS mRNA and protein sequence data to support bovine genome assembly and functional genomics studies. Consequently, the production of high-quality bovine full-CDS cDNA sequences will enhance the bovine genome assembly and functional studies of bovine genes and gene products. The goal of this investigation was to identify and characterize the full-CDS sequences of bovine transcripts from clones identified in non-full-length enriched cDNA libraries. In contrast to several recent full-length cDNA investigations, these full-CDS cDNAs were selected, sequenced, and annotated without the benefit of the target organism's genomic sequence, by using comparison of bovine EST sequence to existing human mRNA to identify likely full-CDS clones for full-length insert cDNA (FLIC) sequencing.

Results

The predicted bovine protein lengths, 5' UTR lengths, and Kozak consensus sequences from 954 bovine FLIC sequences (bFLICs; average length 1713 nt, representing 762 distinct loci) are all consistent with previously sequenced mammalian full-length transcripts.

Conclusion

In most cases, the bFLICs span the entire CDS of the genes, providing the basis for creating predicted bovine protein sequences to support proteomics and comparative evolutionary research as well as functional genomics and genome annotation. The results demonstrate the utility of the comparative approach in obtaining predicted protein sequences in other species.

Construction and validation of a Bovine Innate Immune Microarray

BACKGROUND

Microarray transcript profiling has the potential to illuminate the molecular processes that are involved in the responses of cattle to disease challenges. This knowledge may allow the development of strategies that exploit these genes to enhance resistance to disease in an individual or animal population.

RESULTS

The Bovine Innate Immune Microarray developed in this study consists of 1480 characterised genes identified by literature searches, 31 positive and negative control elements and 5376 cDNAs derived from subtracted and normalised libraries. The cDNA libraries were produced from 'challenged' bovine epithelial and leukocyte cells. The microarray was found to have a limit of detection of 1 pg/microg of total RNA and a mean slide-to-slide correlation co-efficient of 0.88. The profiles of differentially expressed genes from Concanavalin A (ConA) stimulated bovine peripheral blood lymphocytes were determined. Three distinct profiles highlighted 19 genes that were rapidly up-regulated within 30 minutes and returned to basal levels by 24 h; 76 genes that were up-regulated between 2-8 hours and sustained high levels of expression until 24 h and 10 genes that were down-regulated. Quantitative real-time RT-PCR on selected genes was used to confirm the results from the microarray analysis. The results indicate that there is a dynamic process involving gene activation and regulatory mechanisms re-establishing homeostasis in the ConA activated lymphocytes. The Bovine Innate Immune Microarray was also used to determine the cross-species hybridisation capabilities of an ovine PBL sample.

CONCLUSION

The Bovine Innate Immune Microarray has been developed which contains a set of well-characterised genes and anonymous cDNAs from a number of different bovine cell types. The microarray can be used to determine the gene expression profiles underlying innate immune responses in cattle and sheep.

Comparative gene-based in silico analysis of transcriptomes in different bovine tissues and (or) organs

A gene-based approach was used to annotate 322,168 cattle expressed sequence tags (ESTs) based on human genes in order to census the transcriptomes, analyze their expression similarities, and identify genes preferentially expressed in different bovine tissues and (or) organs. Of the 34,157 human coding genes used in a standalone BLAST search, 14,928 could be matched with provisional orthologous sequences in a total of 230,135 bovine ESTs. The remaining 92,033 bovine ESTs were estimated to represent an additional 5970 genes in cattle. On average, approximately 8600 genes were estimated to be expressed in a single tissue and (or) organ and 13,000 in a pooled tissue library. On the basis of the estimated numbers of genes, no more than 3% of genes would be missed when approximately 34,000 ESTs were sequenced from a single tissue and (or) organ library and approximately 40,000 ESTs from a pooled source, respectively. Cluster analyses of the gene expression patterns among 12 single tissues and (or) organs in cattle revealed that their expression similarities would depend on physiological functions. In addition, a total of 1502 genes were identified as preferentially expressed genes in these 12 single tissues and (or) organs with LOD (logarithm of the odds, base 10) > or = 3.0. Therefore, our study provides some insights for further investigating the developmental and functional relations of various tissues and organs in mammals.

Characterization of a normalized cDNA library from bovine intestinal muscle and epithelial tissues

Tissue-specific cDNA library sequences (expressed sequence tags, or EST) yield a detailed snapshot of gene expression and are useful in developing second-generation molecular resources (i.e., microarrays) for gene expression profiling. The objective of this study was to develop and characterize an intestine-specific cDNA library to examine the transcriptome of the bovine gut and identify expressed genes that influence ruminant nutrition and health. We describe BARC-8BOV, a normalized cDNA library developed from mRNA isolated from four distinct intestinal locations (duodenal, jejunal and ileal small intestine, colon) of Holstein dairy cattle resulting in 19,110 5'-EST deposited into the NCBI GenBank EST database. Assembly and clustering of these 19,110 clone sequences yielded 11,208 unique elements (3,419 contigs and 7,789 singletons) with an average length of 695 base pairs. Analysis strongly suggests normalization and tissue pooling were effective at increasing the discovery rate of new bovine sequence. A total of 1,123 sequence elements not previously identified in cattle, but with similarity to known genes in other animal species, were identified and shown to be involved in numerous critical biological processes. An additional 745 transcripts were not previously represented as EST in nucleotide or protein databases, and further analysis of these could lead to the identification of gut-specific transcript variants of known genes or potentially the discovery of novel bovine genes. Of the 11,208 assembled sequences, 11,034, or 98.4%, match sequences present in the bovine DNA trace archive at NCBI, and add to a bovine EST database previously lacking significant gut tissue representation. Ultimately, these data will also contribute in efforts to annotate the bovine genome.

The expanding role of microarrays in the investigation of macrophage responses to pathogens

In the last few years, microarray technology has emerged as the method of choice for large-scale gene expression studies. It provides an efficient and rapid method to investigate the entire transcriptome of a cell. No research field has benefited more from microarray technology than the study of the exquisite interplay between pathogens and hosts. Numerous microarray studies have now been published in this field, which have provided insights into the mechanisms of host defence and the tactics employed by pathogens to circumvent these protection strategies. These studies have led to a more comprehensive understanding of the host immune response and identified new avenues of research for potential control strategies against pathogens. In the past, research has concentrated on human and mouse microarrays to investigate host-pathogen interactions, regardless of the host species. This trend is changing with the ever-expanding sequence resources now available for many pathogen and host species, including livestock animals. The use of species-specific microarrays has furthered our understanding of host-pathogen interactions for particular organisms and aided in the annotation of unknown genes. Macrophages play a central role in the host's innate and adaptive immune responses to pathogens. These cells are in the first line of defence and interact with a wide range of pathogens; many of which have evolved strategies to circumvent the macrophage defence mechanisms and survive within these cells. In this report, we review the wealth of studies using microarray technology to investigate the response of macrophages to pathogens. These studies illustrate how microarray technology has expanded our understanding of the dialogue between macrophage and pathogen and provide examples of the benefits and pitfalls of using this technique. Furthermore, we discuss the resources available to use microarray analysis to study the immune response of a non-human, non-rodent species, the cow.

Generation of a bovine oocyte cDNA library and microarray: resources for identification of genes important for follicular development and early embryogenesis

The oocyte is a key regulator of ovarian folliculogenesis and early embryonic development. However, the composition of the oocyte transcriptome and identities and functions of key oocyte-specific genes involved in the above processes are relatively unknown. Using a PCR-based cDNA amplification method (SMART technology), we constructed a bovine oocyte cDNA library. Analysis of 230 expressed sequence tags (ESTs) from this library identified 102 unique sequences. Although some correspond to housekeeping genes (e.g., ribosomal protein L15) and some represent genes previously known to be expressed in oocytes and other tissues, most encode for genes whose expression in mammalian oocytes has not been reported previously (e.g., cocaine- and amphetamine-regulated transcript) or genes of unknown function. Sixteen did not show significant sequence similarity to any entries in the GenBank database and were classified as novel. Using over 2,000 unsequenced, randomly selected cDNA clones from the library, we constructed an oocyte microarray and performed experiments to identify genes preferentially expressed in fetal ovary (an enriched source of oocytes) relative to somatic tissues. Eleven clones were identified by microarray analysis with consistently higher expression in fetal ovaries (collected from animals at days 210-260 of gestation) compared with spleen and liver. DNA sequence analysis of these clones revealed that two correspond to JY-1, a novel bovine oocyte-specific gene. The remaining nine clones represent five identified genes and one additional completely novel gene. Increased abundance of mRNA in fetal ovary for five of the six genes identified was confirmed by real-time PCR. Results demonstrate the potential utility of these unique resources for identification of oocyte-expressed genes potentially important for reproductive function.

Correlation of amino acid preference and mammalian viral genome type

MOTIVATION

In the event of an outbreak of a disease caused by an initially unknown pathogen, the ability to characterize anonymous sequences prior to isolation and culturing of the pathogen will be helpful. We show that it is possible to classify viral sequences by genome type (dsDNA, ssDNA, ssRNA positive strand, ssRNA negative strand, retroid) using amino acid distribution.

RESULTS

In this paper we describe the results of analysis of amino acid preference in mammalian viruses. The study was carried out at the genome level as well as two shorter sequence levels: short (300 amino acids) and medium length (660 amino acids). The analysis indicates a correlation between the viral genome types dsDNA, ssDNA, ssRNA positive strand, ssRNA negative strand and retroid and amino acid preference. We investigated three different models of amino acid preference. The simplest amino acid preference model, 1-AAP, is a normalized description of the frequency of amino acids in genomes of a viral genome type. A slightly more complex model is the ordered pair amino acid preference model (2-AAP), which characterizes genomes of different viral genome types by the frequency of ordered pairs of amino acids. The most complex and accurate model is the ordered triple amino acid preference model (3-AAP), which is based on ordered triples of amino acids. The results demonstrate that mammalian viral genome types differ in their amino acid preference.

AVAILABILITY

The tools used to format and analyze data and supplementary material are available at http://www.cse.sc.edu/~rose/aminoPreference/index.html

CONTACT

rose@cse.sc.edu.

Interspecific chromosome-wide transcription profiles reveal the existence of mammalian-specific and species-specific chromosome domains

A long-range exploration of expression levels through wide chromosome territories was carried out in three species (pig, cattle, and chicken) by aligning EST counts against the human genome. This strategy made it possible to produce expression profiles that were very similar between pig and cattle and that were significantly correlated with chicken levels of expression. In parallel with these alignments, we developed a statistical approach enabling us to screen genomic regions for both underexpression and overexpression at the chromosome level within a given species, as well as interspecifically. The observed correlations are indicative of the existence of interspecifically conserved domains of gene expression, not only for housekeeping genes (which are highly expressed), but also for regions where genes are significantly underexpressed. Furthermore, our strategy made it possible to point out regions that are differentially regulated between species. These expression data were crossed with available comparative mapping information for pigs and cattle, suggesting that coregulated regions are syntenic in various mammals.

A collection of bovine cDNA probes for gene expression profiling in muscle

Array technology has been increasingly used to monitor global gene expression patterns in various tissues and cell types. However, applications to muscle development and pathology as well as meat production in livestock species have been hampered by the lack of appropriate cDNA collections. To overcome this problem, a directed cDNA library was constructed starting from 23 muscles of meat-producing bovines to derive a collection of 3573 clones. A preliminary sequence characterization of this collection indicated that the most abundant transcripts correspond to genes encoding proteins involved in energy metabolism (COX and NADH dehydrogenase subunits) and belonging to the contractile apparatus (myosin chains and troponin isoforms). From this cDNA library, we selected a set of 435 clones representing 340 unique genes, of which 24 were novel. This collection was subsequently completed with 75 specific cDNA probes for genes of interest already studied in our laboratory. The bovine 'muscle' cDNA repertoire thus designed was spotted onto a nylon membrane (macroarray) in order to test its utility to further investigate the transcriptome of bovine muscles in relation to meat quality traits. It is also anticipated that this type of collection might be useful for the study of chronic myologic diseases in other mammalian species, including humans.

cDNA microarray analysis of bovine embryo gene expression profiles during the pre-implantation period

Background

After fertilization, embryo development involves differentiation, as well as development of the fetal body and extra-embryonic tissues until the moment of implantation. During this period various cellular and molecular changes take place with a genetic origin, e.g. the elongation of embryonic tissues, cell-cell contact between the mother and the embryo and placentation. To identify genetic profiles and search for new candidate molecules involved during this period, embryonic gene expression was analyzed with a custom designed utero-placental complementary DNA (cDNA) microarray.

Methods

Bovine embryos on days 7, 14 and 21, extra-embryonic membranes on day 28 and fetuses on days 28 were collected to represent early embryo, elongating embryo, pre-implantation embryo, post-implantation extra-embryonic membrane and fetus, respectively. Gene expression at these different time points was analyzed using our cDNA microarray. Two clustering algorithms such as k-means and hierarchical clustering methods identified the expression patterns of differentially expressed genes across pre-implantation period. Novel candidate genes were confirmed by real-time RT-PCR.

Results

In total, 1,773 individual genes were analyzed by complete k-means clustering. Comparison of day 7 and day 14 revealed most genes increased during this period, and a small number of genes exhibiting altered expression decreased as gestation progressed. Clustering analysis demonstrated that trophoblast-cell-specific molecules such as placental lactogens (PLs), prolactin-related proteins (PRPs), interferon-tau, and adhesion molecules apparently all play pivotal roles in the preparation needed for implantation, since their expression was remarkably enhanced during the pre-implantation period. The hierarchical clustering analysis and RT-PCR data revealed new functional roles for certain known genes (dickkopf-1, NPM, etc) as well as novel candidate genes (AW464053, AW465434, AW462349, AW485575) related to already established trophoblast-specific genes such as PLs and PRPs.

Conclusions

A large number of genes in extra-embryonic membrane increased up to implantation and these profiles provide information fundamental to an understanding of extra-embryonic membrane differentiation and development. Genes in significant expression suggest novel molecules in trophoblast differentiation.

Integrating linkage and radiation hybrid mapping data for bovine chromosome 15

BACKGROUND

Bovine chromosome (BTA) 15 contains a quantitative trait loci (QTL) for meat tenderness, as well as several breaks in synteny with human chromosome (HSA) 11. Both linkage and radiation hybrid (RH) maps of BTA 15 are available, but the linkage map lacks gene-specific markers needed to identify genes underlying the QTL, and the gene-rich RH map lacks associations with marker genotypes needed to define the QTL. Integrating the maps will provide information to further explore the QTL as well as refine the comparative map between BTA 15 and HSA 11. A recently developed approach to integrating linkage and RH maps uses both linkage and RH data to resolve a consensus marker order, rather than aligning independently constructed maps. Automated map construction procedures employing this maximum-likelihood approach were developed to integrate BTA RH and linkage data, and establish comparative positions of BTA 15 markers with HSA 11 homologs.

RESULTS

The integrated BTA 15 map represents 145 markers; 42 shared by both data sets, 36 unique to the linkage data and 67 unique to RH data. Sequence alignment yielded comparative positions for 77 bovine markers with homologs on HSA 11. The map covers approximately 32% of HSA 11 sequence in five segments of conserved synteny, another 15% of HSA 11 is shared with BTA 29. Bovine and human order are consistent in portions of the syntenic segments, but some rearrangement is apparent. Comparative positions of gene markers near the meat tenderness QTL indicate the region includes separate segments of HSA 11. The two microsatellite markers flanking the QTL peak are between defined syntenic segments.

CONCLUSIONS

Combining data to construct an integrated map not only consolidates information from different sources onto a single map, but information contributed from each data set increases the accuracy of the map. Comparison of bovine maps with well annotated human sequence can provide useful information about genes near mapped bovine markers, but bovine gene order may be different than human. Procedures to connect genetic and physical mapping data, build integrated maps for livestock species, and connect those maps to more fully annotated sequence can be automated, facilitating the maintenance of up-to-date maps, and providing a valuable tool to further explore genetic variation in livestock.