Genomic Scaffolds Research Articles

High-coverage sequencing and annotated assemblies of the budgerigar genome.

BackgroundParrots belong to a group of behaviorally advanced vertebrates and have an advanced ability of vocal learning relative to other vocal-learning birds. They can imitate human speech, synchronize their body movements to a rhythmic beat, and understand complex concepts of referential meaning to sounds. However, little is known about the genetics of these traits. Elucidating the genetic bases would require whole genome sequencing and a robust assembly of a parrot genome.FindingsWe present a genomic resource for the budgerigar, an Australian Parakeet (Melopsittacus undulatus) -- the most widely studied parrot species in neuroscience and behavior. We present genomic sequence data that includes over 300× raw read coverage from multiple sequencing technologies and chromosome optical maps from a single male animal. The reads and optical maps were used to create three hybrid assemblies representing some of the largest genomic scaffolds to date for a bird; two of which were annotated based on similarities to reference sets of non-redundant human, zebra finch and chicken proteins, and budgerigar transcriptome sequence assemblies. The sequence reads for this project were in part generated and used for both the Assemblathon 2 competition and the first de novo assembly of a giga-scale vertebrate genome utilizing PacBio single-molecule sequencing.ConclusionsAcross several quality metrics, these budgerigar assemblies are comparable to or better than the chicken and zebra finch genome assemblies built from traditional Sanger sequencing reads, and are sufficient to analyze regions that are difficult to sequence and assemble, including those not yet assembled in prior bird genomes, and promoter regions of genes differentially regulated in vocal learning brain regions. This work provides valuable data and material for genome technology development and for investigating the genomics of complex behavioral traits.

Genetic Mapping Using Genotyping‐by‐Sequencing in the Clonally Propagated Cassava

ABSTRACTCassava (Manihot esculenta L.) is one of the most important food crops in the tropics, but yields are far below their potential. The gene pool of cassava contains natural genetic diversity relevant to many important breeding goals, but breeding progress has been slow, partly because of insufficient genomic resources. As a first step toward implementing genomewide genetic studies that will facilitate rapid genetic gain through breeding, we genotyped‐by‐sequencing a set of 182 full‐sibs population of cassava that segregated in several traits: resistance to the cassava mosaic disease (CMD) and yield under CMD pressure; increased carotenoid content in storage roots; color of stem exterior and anthocyanin pigmentation in the petioles, inner root skin, and apical leaves. Employing a rare‐cutting restriction enzyme, PstI, in a genotyping‐by‐sequencing (GBS) library preparation, we obtained 2478 segregating single nucleotide polymorphisms (SNPs), of which 1257 passed standard filtering for missing genotypes and deviation from expected genotypic frequencies. We mapped 772 SNPs across 19 linkage groups and anchored 313 unique scaffolds from the version 4.1 of the cassava genome assembly. Most of the studied morphological traits as well as resistance to CMD and root carotenoid content showed qualitative inheritance. As expected, quantitative trait loci analysis for these traits revealed single loci surrounded by small confidence intervals. Yield under CMD was associated with the CMD resistance locus. We show that GBS is a powerful genotyping tool that provides a sufficient number of markers for unraveling the genetic architecture of Mendelian traits in cassava in addition to the development of a robust genetic map that can help anchor unassembled genomic scaffolds.

Genome-wide survey of Alternative Splicing in Sorghum Bicolor.

Sorghum bicolor is a member of grass family which is an attractive model plant for genome study due to interesting genome features like low genome size. In this research, we performed comprehensive investigation of Alternative Splicing and ontology aspects of genes those have undergone these events in sorghum bicolor. We used homology based alignments between gene rich transcripts, represented by tentative consensus (TC) transcript sequences, and genomic scaffolds to deduce the structure of genes and identify alternatively spliced transcripts in sorghum. Using homology mapping of assembled expressed sequence tags with genomics data, we identified 2,137 Alternative Splicing events in S. bicolor. Our study showed that complex events and intron retention are the main types of Alternative Splicing events in S. bicolor and highlights the prevalence of splicing site recognition for definition of introns in this plant. Annotations of the alternatively spliced genes revealed that they represent diverse biological process and molecular functions, suggesting a fundamental role for Alternative Splicing in affecting the development and physiology of S. bicolor.

Sequence and expression analyses of ethylene response factors highly expressed in latex cells from Hevea brasiliensis.

The AP2/ERF superfamily encodes transcription factors that play a key role in plant development and responses to abiotic and biotic stress. In Hevea brasiliensis, ERF genes have been identified by RNA sequencing. This study set out to validate the number of HbERF genes, and identify ERF genes involved in the regulation of latex cell metabolism. A comprehensive Hevea transcriptome was improved using additional RNA reads from reproductive tissues. Newly assembled contigs were annotated in the Gene Ontology database and were assigned to 3 main categories. The AP2/ERF superfamily is the third most represented compared with other transcription factor families. A comparison with genomic scaffolds led to an estimation of 114 AP2/ERF genes and 1 soloist in Hevea brasiliensis. Based on a phylogenetic analysis, functions were predicted for 26 HbERF genes. A relative transcript abundance analysis was performed by real-time RT-PCR in various tissues. Transcripts of ERFs from group I and VIII were very abundant in all tissues while those of group VII were highly accumulated in latex cells. Seven of the thirty-five ERF expression marker genes were highly expressed in latex. Subcellular localization and transactivation analyses suggested that HbERF-VII candidate genes encoded functional transcription factors.

Fine physical and genetic mapping of powdery mildew resistance gene MlIW172 originating from wild emmer (Triticum dicoccoides).

Powdery mildew, caused by Blumeria graminis f. sp. tritici, is one of the most important wheat diseases in the world. In this study, a single dominant powdery mildew resistance gene MlIW172 was identified in the IW172 wild emmer accession and mapped to the distal region of chromosome arm 7AL (bin7AL-16-0.86-0.90) via molecular marker analysis. MlIW172 was closely linked with the RFLP probe Xpsr680-derived STS marker Xmag2185 and the EST markers BE405531 and BE637476. This suggested that MlIW172 might be allelic to the Pm1 locus or a new locus closely linked to Pm1. By screening genomic BAC library of durum wheat cv. Langdon and 7AL-specific BAC library of hexaploid wheat cv. Chinese Spring, and after analyzing genome scaffolds of Triticum urartu containing the marker sequences, additional markers were developed to construct a fine genetic linkage map on the MlIW172 locus region and to delineate the resistance gene within a 0.48 cM interval. Comparative genetics analyses using ESTs and RFLP probe sequences flanking the MlIW172 region against other grass species revealed a general co-linearity in this region with the orthologous genomic regions of rice chromosome 6, Brachypodium chromosome 1, and sorghum chromosome 10. However, orthologous resistance gene-like RGA sequences were only present in wheat and Brachypodium. The BAC contigs and sequence scaffolds that we have developed provide a framework for the physical mapping and map-based cloning of MlIW172.

SNP markers‐based map construction and genome‐wide linkage analysis in Brassica napus

An Illumina Infinium array comprising 5306 single nucleotide polymorphism (SNP) markers was used to genotype 175 individuals of a doubled haploid population derived from a cross between Skipton and Ag-Spectrum, two Australian cultivars of rapeseed (Brassica napus L.). A genetic linkage map based on 613 SNP and 228 non-SNP (DArT, SSR, SRAP and candidate gene markers) covering 2514.8 cM was constructed and further utilized to identify loci associated with flowering time and resistance to blackleg, a disease caused by the fungus Leptosphaeria maculans. Comparison between genetic map positions of SNP markers and the sequenced Brassica rapa (A) and Brassica oleracea (C) genome scaffolds showed several genomic rearrangements in the B.napus genome. A major locus controlling resistance to L.maculans was identified at both seedling and adult plant stages on chromosome A07. QTL analyses revealed that up to 40.2% of genetic variation for flowering time was accounted for by loci having quantitative effects. Comparative mapping showed Arabidopsis and Brassica flowering genes such as Phytochrome A/D, Flowering Locus C and agamous-Like MADS box gene AGL1 map within marker intervals associated with flowering time in a DH population from Skipton/Ag-Spectrum. Genomic regions associated with flowering time and resistance to L.maculans had several SNP markers mapped within 10cM. Our results suggest that SNP markers will be suitable for various applications such as trait introgression, comparative mapping and high-resolution mapping of loci in B.napus.

OMIGA: Optimized Maker-Based Insect Genome Annotation

Insects are one of the largest classes of animals on Earth and constitute more than half of all living species. The i5k initiative has begun sequencing of more than 5,000 insect genomes, which should greatly help in exploring insect resource and pest control. Insect genome annotation remains challenging because many insects have high levels of heterozygosity. To improve the quality of insect genome annotation, we developed a pipeline, named Optimized Maker-Based Insect Genome Annotation (OMIGA), to predict protein-coding genes from insect genomes. We first mapped RNA-Seq reads to genomic scaffolds to determine transcribed regions using Bowtie, and the putative transcripts were assembled using Cufflink. We then selected highly reliable transcripts with intact coding sequences to train de novo gene prediction software, including Augustus. The re-trained software was used to predict genes from insect genomes. Exonerate was used to refine gene structure and to determine near exact exon/intron boundary in the genome. Finally, we used the software Maker to integrate data from RNA-Seq, de novo gene prediction, and protein alignment to produce an official gene set. The OMIGA pipeline was used to annotate the draft genome of an important insect pest, Chilo suppressalis, yielding 12,548 genes. Different strategies were compared, which demonstrated that OMIGA had the best performance. In summary, we present a comprehensive pipeline for identifying genes in insect genomes that can be widely used to improve the annotation quality in insects. OMIGA is provided at http://ento.njau.edu.cn/omiga.html .

Isolation of a Variant Strain of Pleurotus eryngii and the Development of Specific DNA Markers to Identify the Variant Strain

A degenerated strain of Pleurotus eryngii KNR2312 was isolated from a commercial farm. Random amplified polymorphic DNA analysis performed on the genomic DNA of the normal and degenerated strains of this species revealed differences in the DNA banding pattern. A unique DNA fragment (1.7 kbp), which appeared only in the degenerated strain, was isolated and sequenced. Comparing this sequence with the KNR2312 genomic sequence showed that the sequence of the degenerated strain comprised three DNA regions that originated from nine distinct scaffolds of the genomic sequence, suggesting that chromosome-level changes had occurred in the degenerated strain. Using the unique sequence, three sets of PCR primers were designed that targeted the full length, the 5' half, and the 3' half of the DNA. The primer sets P2-1 and P2-2 yielded 1.76 and 0.97 kbp PCR products, respectively, only in the case of the degenerated strain, whereas P2-3 generated a 0.8 kbp product in both the normal and the degenerated strains because its target region was intact in the normal strain as well. In the case of the P2-1 and P2-2 sets, the priming regions of the forward and reverse primers were located at distinct genomic scaffolds in the normal strain. These two primer sets specifically detected the degenerate strain of KNR2312 isolated from various mushrooms including 10 different strains of P. eryngii, four strains of P. ostreatus, and 11 other wild mushrooms.

Single Cell Genomics of Uncultured, Health-Associated Tannerella BU063 (Oral Taxon 286) and Comparison to the Closely Related Pathogen Tannerella forsythia

The uncultivated bacterium Tannerella BU063 (oral taxon 286) is the closest relative to the periodontal pathogen Tannerella forsythia, but is not disease-associated itself. Using a single cell genomics approach, we isolated 12 individual BU063 cells by flow cytometry, and we amplified and sequenced their genomes. Comparative analyses of the assembled genomic scaffolds and their gene contents allowed us to study the diversity of this taxon within the oral community of a single human donor that provided the sample. Eight different BU063 genotypes were represented, all about 5% divergent at the nucleotide level. There were 2 pairs of cells and one group of three that were more highly identical, and may represent clonal populations. We did pooled assemblies on the nearly identical genomes to increase the assembled genomic coverage. The presence of a set of 66 “core” housekeeping genes showed that two of the single cell assemblies and the assembly derived from the three putatively identical cells were essentially complete. As expected, the genome of BU063 is more similar to Tannerella forsythia than any other known genome, although there are significant differences, including a 44% difference in gene content, changes in metabolic pathways, loss of synteny, and an 8–9% difference in GC content. Several identified virulence genes of T. forsythia are not found in BU063 including karilysin, prtH, and bspA. The absence of these genes may explain the lack of periodontal pathogenesis by this species and provides a new foundation to further understand the genome evolution and mechanisms of bacterial-host interaction in closely related oral microbes with different pathogenicity potential.

De novo transcriptome assembly reveals sex-specific selection acting on evolving neo-sex chromosomes in Drosophila miranda

BackgroundThe Drosophila miranda neo-sex chromosome system is a useful resource for studying recently evolved sex chromosomes. However, the neo-Y genomic assembly is fragmented due to the accumulation of repetitive sequence. Furthermore, the separate assembly of the neo-X and neo-Y chromosomes into genomic scaffolds has proven to be difficult, due to their low level of sequence divergence, which in coding regions is about 1.5%. Here, we de novo assemble the transcriptome of D. miranda using RNA-seq data from several male and female tissues, and develop a bioinformatic pipeline to separately reconstruct neo-X and neo-Y transcripts.ResultsWe obtain 2,141 transcripts from the neo-X and 1,863 from the neo-Y. Neo-Y transcripts are generally shorter than their homologous neo-X transcripts (N50 of 2,048-bp vs. 2,775-bp) and expressed at lower levels. We find that 24% of expressed neo-Y transcripts harbor nonsense mutation within their open reading frames, yet most non-functional neo-Y genes are expressed throughout all of their length. We find evidence of gene loss of male-specific genes on the neo-X chromosome, and transcriptional silencing of testis-specific genes from the neo-X.ConclusionsNonsense mediated decay (NMD) has been implicated to degrade transcripts containing pre-mature termination codons (PTC) in Drosophila, but rampant description of neo-Y genes with pre-mature stop codons suggests that it does not play a major role in down-regulating transcripts from the neo-Y. Loss or transcriptional down-regulation of genes from the neo-X with male-biased function provides evidence for beginning demasculinization of the neo-X. Thus, evolving sex chromosomes can rapidly shift their gene content or patterns of gene expression in response to their sex-biased transmission, supporting the idea that sex-specific or sexually antagonistic selection plays a major role in the evolution of heteromorphic sex chromosomes.

Mapping Genes to Chicken Microchromosome 16 and Discovery of Olfactory and Scavenger Receptor Genes Near the Major Histocompatibility Complex

Trisomy mapping is a powerful method for assigning genes to chicken microchromosome 16 (GGA 16). The single chicken nucleolar organizer region (NOR), the 2 major histocompatibility complex regions (MHC-Y and MHC-B), and CD1 genes were all previously assigned to GGA 16 using trisomy mapping. Here, we combined array comparative genomic hybridization with trisomy mapping to screen unassigned genomic scaffolds (consigned temporarily to chrUn_random) for sequences originating from GGA 16. A number of scaffolds mapped to GGA 16. Among these were scaffolds that contain genes for olfactory (OR) and cysteine-rich domain scavenger (SRCR) receptors, along with a number of genes that encode putative immunoglobulin-like receptors and other molecules. We used high-resolution cytogenomic analyses to confirm assignment of OR and SRCR genes to GGA 16 and to pinpoint members of these gene families to the q-arm in partially overlapping regions between the centromere and the NOR. Southern blots revealed sequence polymorphism within the OR/SRCR region and linkage with the MHC-Y region, thereby providing evidence for conserved linkage between OR genes and the MHC within birds. This work localizes OR genes to the vicinity of the chicken MHC and assigns additional genes, including immune defense genes, to GGA 16.

High-throughput genome scaffolding from in vivo DNA interaction frequency

Despite advances in DNA-sequencing technology, assembly of complex genomes remains a major challenge, particularly for genomes sequenced using short reads, which yield highly fragmented assemblies. Here we show that genome-wide in vivo chromatin interaction frequency data, which are measurable with chromosome conformation capture–based experiments, can be used as genomic distance proxies to accurately position individual contigs without requiring any sequence overlap. We also use these data to construct approximate genome scaffolds de novo. Applying our approach to incomplete regions of the human genome, we predict the positions of 65 previously unplaced contigs, in agreement with alternative methods in 26/31 cases attempted in common. Our approach can theoretically bridge any gap size and should be applicable to any species for which global chromatin interaction data can be generated.

Comparative Genomics Analysis of Trichoderma reesei Strains

Abstract Trichoderma reesei is a key fungus for industrial production of lignocellulolytic enzymes. The genome sequences of the T. reesei hyper-cellulolytic strain RUT-C30 and its parental strain QM6a were compared at the nucleotide level. Approximately 97% of the 87 genomic-sequence scaffolds of T. reesei QM6a (33Mb) were found to have the corresponding nucleotide in the 182 genome-sequence scaffolds of RUT-C30 (32Mb). There are 455 loci within the QM6 sequence not detected in the RUT-C30 sequence. Regions at the termini of QM6a scaffolds as well as 14 small scaffolds do not have corresponding regions in RUT-C30 genomic scaffolds. Seventy-eight protein-encoding genes are included within these regions. Mutated nucleotide(s) in 2,371 positions, including short insertion/deletions (indels), were detected in the aligned regions. The predicted protein-coding regions of 97 gene models contain mutations, 34 of which were not previously described. Twenty-seven out of 34 newly discovered genes were found to have ...

Characterization of the Z lineage Major histocompatability complex class I genes in zebrafish

Zebrafish (Danio rerio) are a valuable model for studying immunity, infection, and hematopoietic disease and have recently been employed for transplantation assays. However, the lack of syngeneic zebrafish creates challenges with identifying immune-matched individuals. The MHC class I genes, which mediate allogeneic recognition in mammals, have been grouped into three broad lineages in zebrafish: the classical U genes on chromosome 19, the Z genes which have been reported to map to chromosome 1, and the L genes that map to multiple loci. Transplantations between individual zebrafish that are matched at the U locus fail to consistently engraft suggesting that additional loci contribute to allogeneic recognition. Although two full-length zebrafish Z transcripts have been described, the genomic organization and diversity of these genes have not been reported. Herein we define ten Z genes on chromosomes 1 and 3 and on an unplaced genomic scaffold. We report that neither of the Z transcripts previously described match the current genome assembly and classify these transcripts as additional gene loci. We characterize full-length transcripts for 9 of these 12 genes. We demonstrate a high level of expression variation of the Z genes between individual zebrafish suggestive of haplotypic variation. We report low level sequence variation for individual Z genes between individual zebrafish reflecting a possible nonclassical function, although these molecules may still contribute to allogeneic recognition. Finally, we present a gene nomenclature system for the Z genes consistent with MHC nomenclature in other species and with the zebrafish gene nomenclature guidelines.

Generation of Physical Map Contig-Specific Sequences Useful for Whole Genome Sequence Scaffolding

Along with the rapid advances of the nextgen sequencing technologies, more and more species are added to the list of organisms whose whole genomes are sequenced. However, the assembled draft genome of many organisms consists of numerous small contigs, due to the short length of the reads generated by nextgen sequencing platforms. In order to improve the assembly and bring the genome contigs together, more genome resources are needed. In this study, we developed a strategy to generate a valuable genome resource, physical map contig-specific sequences, which are randomly distributed genome sequences in each physical contig. Two-dimensional tagging method was used to create specific tags for 1,824 physical contigs, in which the cost was dramatically reduced. A total of 94,111,841 100-bp reads and 315,277 assembled contigs are identified containing physical map contig-specific tags. The physical map contig-specific sequences along with the currently available BAC end sequences were then used to anchor the catfish draft genome contigs. A total of 156,457 genome contigs (~79% of whole genome sequencing assembly) were anchored and grouped into 1,824 pools, in which 16,680 unique genes were annotated. The physical map contig-specific sequences are valuable resources to link physical map, genetic linkage map and draft whole genome sequences, consequently have the capability to improve the whole genome sequences assembly and scaffolding, and improve the genome-wide comparative analysis as well. The strategy developed in this study could also be adopted in other species whose whole genome assembly is still facing a challenge.

L_RNA_scaffolder: scaffolding genomes with transcripts

BackgroundGeneration of large mate-pair libraries is necessary for de novo genome assembly but the procedure is complex and time-consuming. Furthermore, in some complex genomes, it is hard to increase the N50 length even with large mate-pair libraries, which leads to low transcript coverage. Thus, it is necessary to develop other simple scaffolding approaches, to at least solve the elongation of transcribed fragments.ResultsWe describe L_RNA_scaffolder, a novel genome scaffolding method that uses long transcriptome reads to order, orient and combine genomic fragments into larger sequences. To demonstrate the accuracy of the method, the zebrafish genome was scaffolded. With expanded human transcriptome data, the N50 of human genome was doubled and L_RNA_scaffolder out-performed most scaffolding results by existing scaffolders which employ mate-pair libraries. In these two examples, the transcript coverage was almost complete, especially for long transcripts. We applied L_RNA_scaffolder to the highly polymorphic pearl oyster draft genome and the gene model length significantly increased.ConclusionsThe simplicity and high-throughput of RNA-seq data makes this approach suitable for genome scaffolding. L_RNA_scaffolder is available at http://www.fishbrowser.org/software/L_RNA_scaffolder.

The cephalochordate Branchiostoma genome contains 26 intermediate filament (IF) genes: Implications for evolution of chordate IF proteins

We analyzed the draft genome of the cephalochordate Branchiostoma floridae (B. floridae) for genes encoding intermediate filament (IF) proteins. From 26 identified IF genes 13 were not reported before. Four of the new IF genes belong to the previously established Branchiostoma IF group A, four to the Branchiostoma IF group B, one is homologous to the type II keratin E2 while the remaining four new IF sequences N1 to N4 could not be readily classified in any of the previously established Branchiostoma IF groups. All eleven identified A and B2-type IF genes are located on the same genomic scaffold and arose due to multiple cephalochordate-specific duplications. Another IF gene cluster, identified in the B. floridae genome, contains three keratins (E1, Y1, D1), two keratin-like IF genes (C2, X1), one new IF gene (N1) and one IF unrelated gene, but does not show any similarities to the well defined vertebrate type I or type II keratin gene clusters. In addition, some type III sequence features were documented in the new IF protein N2, which, however, seems to share a common ancestry with the Branchiostoma keratins D1 and two keratin-related genes C. Thus, a few type I and type II keratin genes existed in a common ancestor of cephalochordates and vertebrates, which after separation of these two lineages gave rise to the known complexities of the vertebrate cytoplasmic type I–IV IF proteins, as well as to the multiple keratin and related IF genes in cephalochordates, due to multiple gene duplications, deletions and sequence divergences.

Choosing a Benchtop Sequencing Machine to Characterise Helicobacter pylori Genomes

The fully annotated genome sequence of the European strain, 26695 was first published in 1997 and, in 1999, it was directly compared to the USA isolate J99, promoting two standard laboratory isolates for Helicobacter pylori (H. pylori) research. With the genomic scaffolds available from these important genomes and the advent of benchtop high-throughput sequencing technology, a bacterial genome can now be sequenced within a few days. We sequenced and analysed strains J99 and 26695 using the benchtop-sequencing machines Ion Torrent PGM and the Illumina MiSeq Nextera and Nextera XT methodologies. Using publically available algorithms, we analysed the raw data and interrogated both genomes by mapping the data and by de novo assembly. We compared the accuracy of the coding sequence assemblies to the originally published sequences. With the Ion Torrent PGM, we found an inherently high-error rate in the raw sequence data. Using the Illumina MiSeq, we found significantly more non-covered nucleotides when using the less expensive Illumina Nextera XT compared with the Illumina Nextera library creation method. We found the most accurate de novo assemblies using the Nextera technology, however, extracting an accurate multi-locus sequence type was inconsistent compared to the Ion Torrent PGM. We found the cagPAI failed to assemble onto a single contig in all technologies but was more accurate using the Nextera. Our results indicate the Illumina MiSeq Nextera method is the most accurate for de novo whole genome sequencing of H. pylori.

Genomer — A Swiss Army Knife for Genome Scaffolding

The increasing accessibility and reduced costs of sequencing has made genome analysis accessible to more and more researchers. Yet there remains a steep learning curve in the subsequent computational steps required to process raw reads into a database-deposited genome sequence. Here we describe “Genomer,” a tool to simplify the manual tasks of finishing and uploading a genome sequence to a database. Genomer can format a genome scaffold into the common files required for submission to GenBank. This software also simplifies updating a genome scaffold by allowing a human-readable YAML format file to be edited instead of large sequence files. Genomer is written as a command line tool and is an effort to make the manual process of genome scaffolding more robust and reproducible. Extensive documentation and video tutorials are available at http://next.gs.

Genomic Analysis and mRNA Expression of Equine Type I Interferon Genes

This study aimed at identifying all of the type I interferon (IFN) genes of the horse and at monitoring their expression in equine cells on in vitro induction. We identified 32 putative type I IFN loci on horse chromosome 23 and an unplaced genomic scaffold. A phylogentic analysis characterized these into 8 different type I IFN classes, that is, putative functional genes for 6 IFN-α, 4 IFN-β, 8 IFN-ω (plus 4 pseudogenes), 3 IFN-δ (plus 1 pseudogene), 1 IFN-κ and 1 IFN-ε, plus 1 IFN-ν pseudogene, and 3 loci belonging to what has previously been called IFN-αω. Our analyses indicate that the IFN-αω genes are quite distinct from both IFN-α and IFN-ω, and we refer to this type I IFN as IFN-μ. Results from cell cultures showed that leukocytes readily expressed IFN-α, IFN-β, IFN-δ, IFN-μ, and IFN-ω mRNA on induction with, for example, live virus; while fibroblasts only expressed IFN-β mRNA on stimulation. IFN-κ or IFN-ε expression was not consistently induced in these cell cultures. Thus, the equine type I IFN family comprised 8 classes, 7 of which had putative functional genes, and mRNA expression of 5 was induced in vitro. Moreover, a relatively low number of IFN-α subtypes was found in the horse compared with other eutherian mammals.