Dynamic transcriptome and chromatin architecture in granulosa cells during chicken folliculogenesis.

Folliculogenesis is a complex biological process involving a central oocyte and its surrounding somatic cells. Three-dimensional chromatin architecture is an important transcription regulator; however, little is known about its dynamics and role in transcriptional regulation of granulosa cells during chicken folliculogenesis. We investigate the transcriptomic dynamics of chicken granulosa cells over ten follicular stages and assess the chromatin architecture dynamics and how it influences gene expression in granulosa cells at three key stages: the prehierarchical small white follicles, the first largest preovulatory follicles, and the postovulatory follicles. Our results demonstrate the consistency between the global reprogramming of chromatin architecture and the transcriptomic divergence during folliculogenesis, providing ample evidence for compartmentalization rearrangement, variable organization of topologically associating domains, and rewiring of the long-range interaction between promoter and enhancers. These results provide key insights into avian reproductive biology and provide a foundational dataset for the future in-depth functional characterization of granulosa cells.

Three chromosome-level duck genome assemblies provide insights into genomic variation during domestication.

Domestic ducks are raised for meat, eggs and feather down, and almost all varieties are descended from the Mallard (Anas platyrhynchos). Here, we report chromosome-level high-quality genome assemblies for meat and laying duck breeds, and the Mallard. Our new genomic databases contain annotations for thousands of new protein-coding genes and recover a major percentage of the presumed "missing genes" in birds. We obtain the entire genomic sequences for the C-type lectin (CTL) family members that regulate eggshell biomineralization. Our population and comparative genomics analyses provide more than 36 million sequence variants between duck populations. Furthermore, a mutant cell line allows confirmation of the predicted anti-adipogenic function of NR2F2 in the duck, and uncovered mutations specific to Pekin duck that potentially affect adipose deposition. Our study provides insights into avian evolution and the genetics of oviparity, and will be a rich resource for the future genetic improvement of commercial traits in the duck.

Transcriptome Response of Liver and Muscle in Heat-Stressed Laying Hens.

Exposure to high ambient temperature has detrimental effects on poultry welfare and production. Although changes in gene expression due to heat exposure have been well described for broiler chickens, knowledge of the effects of heat on laying hens is still relatively limited. In this study, we profiled the transcriptome for pectoralis major muscle (n = 24) and liver (n = 24), during a 4-week cyclic heating experiment performed on layers in the early phase of egg production. Both heat-control and time-based contrasts were analyzed to determine differentially expressed genes (DEGs). Heat exposure induced different changes in gene expression for the two tissues, and we also observed changes in gene expression over time in the control animals suggesting that metabolic changes occurred during the transition from onset of lay to peak egg production. A total of 73 DEGs in liver were shared between the 3 h heat-control contrast, and the 4-week versus 3 h time contrast in the control group, suggesting a core set of genes that is responsible for maintenance of metabolic homeostasis regardless of the physiologic stressor (heat or commencing egg production). The identified DEGs improve our understanding of the layer's response to stressors and may serve as targets for genetic selection in the future to improve resilience.

Comparative Transcriptome Profiling of Ovary Tissue between Black Muscovy Duck and White Muscovy Duck with High- and Low-Egg Production.

The egg-laying rate is an important indicator for evaluating fertility of poultry. In order to better understand the laying mechanism of Muscovy ducks, gene expression profiles and pathways of ovarian tissues in high- and low-laying black (BH and BL) and white Muscovy ducks (WH and WL) during the peak production period were performed by using RNA-seq. The total number of reads produced for each ovarian sample ranged from 44,344,070 to 47,963,328. A total of 113, 619 and 87 differentially expressed genes (DEGs) were identified in BH-vs-WH, BL-vs-BH and BL-vs-WL, respectively. Among them, 54, 356 and 49 genes were up regulated and 59, 263 and 38 genes were down regulated. In addition, there were only 10 up-regulated genes in WL-vs-WH. In the comparison of DEGs in black and white Muscovy ducks, two co-expressed DEG genes were detected between BH-vs-WH and BL-vs-WL and seven DEGs were co-expressed between BL-vs-BH and WL-vs-WH. The RNA-Seq data were confirmed to be reliable by qPCR. Numerous DEGs known to be involved in ovarian development were identified, including TGFß2, NGFR, CEBPD, CPEB2, POSTN, SMOC1, FGF18, EFNA5 and SDC4. Gene Ontology (GO) annotations indicated that DEGs related to ovarian development were mainly enriched in biological processes of "circadian sleep/wake cycle process," "negative regulation of transforming growth factor-ß secretion," "positive regulation of calcium ion transport" in BH-vs-WH and "cell surface receptor signaling pathway," "Notch signaling pathway" and "calcium ion transport" in BL-vs-BH. Besides, "steroid biosynthetic process," "granulosa cell development" and "egg coat formation" were mainly enriched in BL-vs-WL and "reproduction," "MAPK cascade" and "mitotic cell cycle" were mainly enriched in WL-vs-WH. KEGG pathway analysis showed that the PI3K-Akt signaling pathway and ovarian steroidogenesis were the most enriched in Muscovy duck ovary transcriptome data. This work highlights potential genes and pathways that may affect ovarian development in Muscovy duck.

Transcriptional insights into key genes and pathways controlling muscle lipid metabolism in broiler chickens.

BACKGROUND: Intramuscular fat (IMF) is one of the most important factors positively associated with meat quality. Triglycerides (TGs), as the main component of IMF, play an essential role in muscle lipid metabolism. This transcriptome analysis of pectoralis muscle tissue aimed to identify functional genes and biological pathways likely contributing to the extreme differences in the TG content of broiler chickens. RESULTS: The study included Jingxing-Huang broilers that were significantly different in TG content (5.81 mg/g and 2.26 mg/g, p < 0.01) and deposition of cholesterol also showed the same trend. This RNA sequencing analysis was performed on pectoralis muscle samples from the higher TG content group (HTG) and the lower TG content group (LTG) chickens. A total of 1200 differentially expressed genes (DEGs) were identified between two groups, of which 59 DEGs were related to TG and steroid metabolism. The HTG chickens overexpressed numerous genes related to adipogenesis and lipogenesis in pectoralis muscle tissue, including the key genes ADIPOQ, CD36, FABP4, FABP5, LPL, SCD, PLIN1, CIDEC and PPARG, as well as genes related to steroid biosynthesis (DHCR24, LSS, MSMO1, NSDHL and CH25H). Additionally, key pathways related to lipid storage and metabolism (the steroid biosynthesis and peroxisome proliferator activated receptor (PPAR) signaling pathway) may be the key pathways regulating differential lipid deposition between HTG group and LTG group. CONCLUSIONS: This study showed that increased TG deposition accompanying an increase in steroid synthesis in pectoralis muscle tissue. Our findings of changes in gene expression of steroid biosynthesis and PPAR signaling pathway in HTG and LTG chickens provide insight into genetic mechanisms involved in different lipid deposition patterns in pectoralis muscle tissue.

Molecular cloning, expression and characterization of Pekin duck programmed death-1.

Programmed death-1 (PD-1) has a pivotal role in the attenuation of adaptive immune responses and peripheral tolerance. Here we describe the identification of the Pekin duck programmed death-1 orthologue (duPD-1). The duPD-1 cDNA encodes a 283-amino acid polypeptide that has an amino acid identity of 70%, 32% and 31% with chicken, murine and human PD-1, respectively. The duck PD-1 gene shares five conserved exons with chicken, murine and human PD-1 genes. A cluster of putative regulatory elements within the conserved region B (CR-B) of the basal promotor is conserved. Homology modeling was most compatible with the two ß-sheet IgV domain structure of murine PD-1. Contact residues, shown to be critical for binding of the respective human and murine PD-1 ligands are mostly conserved between avian and mammalian species, whereas residues that define the cytoplasmic immunoreceptor tyrosine-based inhibitory motif (ITIM) and immunoreceptor tyrosine-based switch motif (ITSM) are highly conserved across higher vertebrates and frog. Constitutive expression of duPD-1 transcripts was predominantly found in lymphocyte-rich tissues, and mitogen-stimulation of duck peripheral blood mononuclear cells transiently increased duPD-1 mRNA expression. A soluble duPD-1 protein was expressed and shown to engage the identified duck PD-1 ligands. Our observations show considerable evolutionary conservation between mammalian and avian PD-1 orthologues. This work will facilitate further investigation of the role of PD-1 signaling in adaptive immunity in the Pekin duck, a non-mammalian vertebrate and pathogen host with relevance for human and animal health.

The core promoter controls basal and inducible expression of duck retinoic acid inducible gene-I (RIG-I).

Retinoic acid inducible gene-I (RIG-I) is a cytoplasmic RNA sensor for detecting a variety of RNA viruses including influenza A viruses. Detection ultimately produces Type I interferon (IFN), which stimulates expression of interferon stimulated genes (ISGs), including RIG-I itself in a positive feedback loop. The structure and function of RIG-I is conserved across phylogeny, despite significant protein sequence divergence, however, the promoter sequences do not show the expected phylogenetic relationships and it is not known whether they are similarly regulated. We previously cloned duck RIG-I and showed it is highly induced during influenza A infection consistent with induction by the interferon produced. Here, we identified the Pekin duck RIG-I promoter and constructed promoter reporter vectors, which we transfected into duck embryonic fibroblasts or chicken DF-1 cells and tested in dual luciferase assays. We showed that activation of the Mitochondrial Antiviral Signalling (MAVS) pathway using the constitutively active N-terminal region of RIG-I or polyinosinic-polycytidylic acid (poly I:C) led to stimulation of duck RIG-I promoter activity. Using deletion constructs we showed the core promoter lies in the proximal 250 basepairs, and we identified essential cis-regulatory elements, a GC-box and an interferon-sensitive response element (ISRE), responsible for basal and inducible expression, respectively. Using mCherry-tagged interferon regulatory factors (IRFs) cloned from chickens and ducks, we show overexpression of chIRF7 induced the duck RIG-I promoter, and this required the ISRE site. Finally, we also demonstrated that overexpressed chIRF7 translocated to the nucleus, which was augmented by MAVS activation using RIG-I 2CARD. Our findings demonstrate that RIG-I expression is induced by chIRF7, in a positive regulatory loop. These studies show that the duck RIG-I promoter is appropriately regulated in chicken cells, necessary for the potential generation of transgenic chickens expressing RIG-I.

Identification and functional characterization, including cytokine production modulation, of the novel chicken Interleukin-11.

Interleukin (IL)-11 plays an important role in the immune system. However, IL-11 has not yet been characterized in avian species, including chickens. This study is the first to clone and functionally characterize chicken IL-11 (chIL-11). Multiple alignments and phylogenetic tree comparisons of chIL-11 with IL-11 proteins from other species revealed high levels of conservation and a close relationship between chicken and Japanese quail IL-11. Our results demonstrate that chIL-11 was a functional ligand of IL-11RA and IL-6ST in chicken HD11 and OU2 cell lines, as well as activated and regulated JAK-STAT, NF-κB, PI3K/AKT, and MAPK signaling pathways in chicken cell lines. In addition, chIL-11 inhibited nitric oxide production, affected proliferation of both tested cell lines, inhibited Type 1 and 17 T helper (Th) cytokine and IL-26, IL-12, and IL-17A-induced interferon-γ production, and enhanced Th2 cytokine (IL-4 and IL-10) production. Taken together, functional analysis of chIL-11 revealed it bound to IL-11RA and IL-6ST and activated the JAK-STAT, NF-κB, and MAPK signaling pathways, which resulted in modulation of Th1/Th17 and Th2 cytokine production in chicken HD11 and OU2 cell lines. Overall, this indicates chIL-11 has a role in both the innate and adaptive immune system.

Birds Generally Carry a Small Repertoire of Bitter Taste Receptor Genes.

As they belong to the most species-rich class of tetrapod vertebrates, birds have long been believed to possess an inferior taste system. However, the bitter taste is fundamental in birds to recognize dietary toxins (which are typically bitter) in potential food sources. To characterize the evolution of avian bitter taste receptor genes (Tas2rs) and to test whether dietary toxins have shaped the repertoire size of avian Tas2rs, we examined 48 genomes representing all but 3 avian orders. The total number of Tas2r genes was found to range from 1 in the domestic pigeon to 12 in the bar-tailed trogon, with an average of 4, which suggested that a much smaller Tas2r gene repertoire exists in birds than in other vertebrates. Furthermore, we uncovered a positive correlation between the number of putatively functional Tas2rs and the abundance of potential toxins in avian diets. Because plant products contain more toxins than animal tissues and insects release poisonous defensive secretions, we hypothesized that herbivorous and insectivorous birds may demand more functional Tas2rs than carnivorous birds feeding on noninsect animals. Our analyses appear to support this hypothesis and highlight the critical role of taste perception in birds.

Top-down-assisted bottom-up method for homologous protein sequencing: hemoglobin from 33 bird species.

Ticks are vectors for disease transmission because they are indiscriminant in their feeding on multiple vertebrate hosts, transmitting pathogens between their hosts. Identifying the hosts on which ticks have fed is important for disease prevention and intervention. We have previously shown that hemoglobin (Hb) remnants from a host on which a tick fed can be used to reveal the host's identity. For the present research, blood was collected from 33 bird species that are common in the U.S. as hosts for ticks but that have unknown Hb sequences. A top-down-assisted bottom-up mass spectrometry approach with a customized searching database, based on variability in known bird hemoglobin sequences, has been devised to facilitate fast and complete sequencing of hemoglobin from birds with unknown sequences. These hemoglobin sequences will be added to a hemoglobin database and used for tick host identification. The general approach has the potential to sequence any set of homologous proteins completely in a rapid manner. Graphical Abstract ᅟ.

Molecular characterization and functional analysis of duck TRAF6.

Tumor necrosis factor receptor-associated factor 6 (TRAF6) plays a pivotal role in activating various signaling cascades as an intracellular signal transducer. Although significant progress has been made clarifying TRAF6 function in mammals, the role of TRAF6 in ducks (duTRAF6) remains poorly understood. In the present study, we cloned the full-length duTRAF6 cDNA from duck embryo fibroblasts (DEFs) for the first time. Real-time quantitative reverse transcription-polymerase chain reaction assays showed that duTRAF6 was widely expressed in different tissues. Overexpression of duTRAF6 activated nuclear factor kappa B (NF-κB) and induced interferon-ß expression. Furthermore, a deletion mutant analysis revealed that the duTRAF6 region between aa 115 and 375 was essential for activating NF-κB. In addition, duTRAF6 knockdown by RNA interference significantly reduced poly(I:C)- and Sendai virus-induced NF-κB activation in DEFs. Taken together, our results demonstrate that duTRAF6 plays a crucial immunoregulatory role in the duck innate immune response.

Molecular cloning, tissue distribution, and immune function of goose TLR7.

TLR7 is a transmembrane endosomal protein that plays an essential role in innate antiviral responses via the recognition of conserved viral molecular patterns. Here, we cloned the full-length cDNA of goose TLR7 and carried out a molecular characterization of goose TLR7. The goose TLR7 gene is 3900 bp and encodes a 1045 amino acid protein with high homology to poultry (93% to duck and 83% to chicken). Similar conclusions were made by phylogenetic analysis. The predicted protein secondary structure of goose TLR7 contained a conserved Toll/interleukin-1 receptor domain and characteristic leucine-rich repeat regions, which has also been reported for duck TLR7. Additionally, the tissue distribution of goose TLR7 suggests that immune-associated tissues, especially the cecal tonsil and bursa of Fabricius, have high goose TLR7 expression levels. Goose TLR7 is abundantly expressed in lung tissues, which is distinct from its expression in chickens. Similar to duck TLR7, goose spleen mononuclear cells (MNCs) exposed to the mammalian TLR7 agonists R848 and Imiquimod showed significant induction of the production of proinflammatory cytokines and IFN-α. New type gosling viral enteritis virus (NGVEV) infection resulted in high mRNA expression levels of goose TLR7 in the spleen. By contrast, no direct interaction between NGVEV and goose TLR7 was detected after infecting goose spleen MNCs with NGVEV in vitro. However, triggering of goose TLR7 resulted in the rapid up-regulation of proinflammatory cytokines and anti-viral molecules, suggesting that goose TLR7 plays an important role in anti-viral defense.

Evolution of gonadotropin-inhibitory hormone receptor and its ligand.

Gonadotropin-inhibitory hormone (GnIH) is a neuropeptide inhibitor of gonadotropin secretion, which was first identified in the Japanese quail hypothalamus. GnIH peptides share a C-terminal LPXRFamide (X=L or Q) motif in most vertebrates. The receptor for GnIH (GnIHR) is the seven-transmembrane G protein-coupled receptor 147 (GPR147) that inhibits cAMP production. GPR147 is also named neuropeptide FF (NPFF) receptor 1 (NPFFR1), because it also binds NPFF that has a C-terminal PQRFamide motif. To understand the evolutionary history of the GnIH system in the animal kingdom, we searched for receptors structurally similar to GnIHR in the genome of six mammals (human, mouse, rat, cattle, cat, and rabbit), five birds (pigeon, chicken, turkey, budgerigar, and zebra finch), one reptile (green anole), one amphibian (Western clawed flog), six fishes (zebrafish, Nile tilapia, Fugu, coelacanth, spotted gar, and lamprey), one hemichordate (acorn worm), one echinoderm (purple sea urchin), one mollusk (California sea hare), seven insects (pea aphid, African malaria mosquito, honey bee, buff-tailed bumblebee, fruit fly, jewel wasp, and red flour beetle), one cnidarian (hydra), and constructed phylogenetic trees by neighbor joining (NJ) and maximum likelihood (ML) methods. A multiple sequence alignment of the receptors showed highly conserved seven-transmembrane domains as well as disulfide bridge sites between the first and second extracellular loops, including the receptor of hydra. Both NJ and ML analyses grouped the receptors of vertebrates into NPFFR1 and NPFFR2 (GPR74), and the receptors of insects into the receptor for SIFamide peptides that share a C-terminal YRKPPFNGSIFamide motif. Although human, quail and zebrafish GnIHR (NPFFR1) were most structurally similar to SIFamide receptor of fruit fly in the Famide peptide (FMRFamide, neuropeptide F, short neuropeptide F, drosulfakinin, myosuppressin, SIFamide) receptor families, the amino acid sequences and the peptide coding regions of GnIH precursors were most similar to FMRFamide precursor of fruit fly in the precursors of Famide peptide families. Chromosome synteny analysis of the precursor genes of human, quail and zebrafish GnIH and fruit fly Famide peptides further identified conserved synteny in vertebrate GnIH and fruit fly FMRFa precursor genes as well as other Famide peptide precursor genes. These results suggest that GnIH and its receptor pair and SIFamide and its receptor pair may have diverged and co-evolved independently in vertebrates and insects, respectively, from their ancestral Famide peptide and its receptor pair, during diversification and evolution of deuterostomian and protostomian species.

Cloning, expression and antiviral bioactivity of red-crowned crane interferon-α.

Interferon-α (IFN-α) genes have been cloned from a variety of animals, but information regarding crane IFN-α has not been reported to date. In this study, we cloned a full-length Red-crowned Crane interferon-α (crIFN-α) gene sequence consisting of a 486bp partial 5' UTR, 741bp complete ORF and 559bp partial 3' UTR. This gene encodes a protein of 246 amino acids and shares 60 to 80% identity with avian IFN-α and less than 45% identity with mammalian IFN-α. The expression of crIFN-α with an N-terminal His-tag was investigated in Escherichia coli, and the protein was purified on a nickel column. To obtain activated proteins, crIFN-α inclusion bodies were renatured by dialysis. In vitro cytopathic inhibition assays indicated that the recombinant crIFN-α could inhibit the replication of vesicular stomatitis virus in chicken fibroblasts. These antiviral activities were abrogated by rabbit anti-crIFN-α antibodies in vitro. In addition, an immunofluorescence assay indicated that crIFN-α could be expressed in chicken fibroblasts and was primarily located in the cytoplasm. Taken together, our results suggest that the crIFN-α gene may play an important role in inhibiting the replication of viruses.

Ontogeny of the avian intestinal immunoglobulin repertoire: modification in CDR3 length and conserved VH-pseudogene usage.

Immunoglobulins play an important role in maintenance of mucosal homeostasis in the gut. The antigen binding specificity of these immunoglobulins depends for a large part on the hypervariable CDR3 region. To gain knowledge about isotype-specific development of the CDR3 repertoire we examined CDR3 spectratypes at multiple time points between 4 and 70 days post hatch. In order to identify clonal expansions deviation from the normal distribution (SS) and the average CDR3 length was calculated. IgA-CDR3 regions were studied in more detail by DNA sequence analysis at day 7 and 70 and preferential VH pseudogene usage was estimated. The SS of CDR3 repertoires of the IgM, IgG and IgA isotypes successively increased, but for each isotype this increase was transiently. The length of the CDR3 regions decreased with age for IgM becoming similar to the CDR3 length of IgA at day 70. The IgA- and IgG-CDR3 lengths did not change with age. On average, the CDR3 length of IgA was the shortest. IgA CDR3 sequences were similar between animals aged 7 and 70 days. A limited number of pseudogenes was used, and no differences in pseudogene usage were observed between animals aged 7 and 70 days. Of the identified VH pseudogenes, half of the sequences used VH15, whilst a number of the pseudogenes were not used at all. We conclude that CDR3 spectratype profiles change during aging, whilst at the CDR3-sequence level, variation in VH pseudogene usage for ileal IgA is limited suggesting conservation during ontogeny.

Expression patterns of the opsin 5-related genes in the developing chicken retina.

The opsin gene family encodes G protein-coupled seven-transmembrane proteins that bind to a retinaldehyde chromophore for photoreception. It has been reported that opsin 5 is expressed in mammalian neural tissue, but its function has been elusive. As a first step to understand the function for opsin 5 in the developing eye, we searched for chicken opsin 5-related genes in the genome by a bioinformatic approach and isolated opsin 5 cDNA fragments from the embryonic retina by RT-PCR. We found that there are three opsin 5-related genes, designated cOpn5m (chicken opsin 5, mammalian type), cOpn5L1 (chicken opsin 5-like 1), and cOpn5L2 (chicken opsin 5-like 2), in the chicken genome. Quantitative PCR analysis has revealed that cOpn5m is the most abundant in the developing and early posthatching neural retina. In situ hybridization analysis has shown that cOpn5m is specifically expressed in subsets of differentiating ganglion cells and amacrine cells. These results suggest that the mammalian type opsin 5 may contribute to the development of these retinal cells in the chicken.

Natural selection in avian protein-coding genes expressed in brain.

The evolution of birds from theropod dinosaurs took place approximately 150 million years ago, and was associated with a number of specific adaptations that are still evident among extant birds, including feathers, song and extravagant secondary sexual characteristics. Knowledge about the molecular evolutionary background to such adaptations is lacking. Here, we analyse the evolution of > 5000 protein-coding gene sequences expressed in zebra finch brain by comparison to orthologous sequences in chicken. Mean d(N)/d(S) is 0.085 and genes with their maximal expression in the eye and central nervous system have the lowest mean d(N)/d(S) value, while those expressed in digestive and reproductive tissues exhibit the highest. We find that fast-evolving genes (those which have higher than expected rate of nonsynonymous substitution, indicative of adaptive evolution) are enriched for biological functions such as fertilization, muscle contraction, defence response, response to stress, wounding and endogenous stimulus, and cell death. After alignment to mammalian orthologues, we identify a catalogue of 228 genes that show a significantly higher rate of protein evolution in the two bird lineages than in mammals. These accelerated bird genes, representing candidates for avian-specific adaptations, include genes implicated in vocal learning and other cognitive processes. Moreover, colouration genes evolve faster in birds than in mammals, which may have been driven by sexual selection for extravagant plumage characteristics.

Expression of the chick Sizzled gene in progenitors of the cardiac outflow tract.

The secreted Frizzled receptor related proteins (Sfrp's) belong to a protein family that comprises antagonists and modifiers of Wnt and BMP signalling. Here we report the isolation and expression pattern of the Sfrp gene "Sizzled" in the chick. Sizzled genes, as well as the closely related crescent genes, exist in the genomes of fishes, frogs and chicks, but not of mammals. The chicken Sizzled gene (Szl) is initially expressed in the anterior endoderm of gastrulating and early head fold embryos. An additional, separate expression domain develops at the posterior end of the embryo from the head process stage onwards. Szl transcripts are then detected in precardial mesodermal cells, are transiently transcribed in the straight heart tube, and later prominently in the splanchnic mesoderm surrounding the arterial pole of the developing heart, the anterior heart field. These cells are subsequently recruited to form the cardiac outflow tract. cSzl expression is downregulated when the septation of the outflow tract by neural crest derived cells begins.

Evolutionary patterns of MHC class II B in owls and their implications for the understanding of avian MHC evolution.

Owing to its special mode of evolution and central role in the adaptive immune system, the major histocompatibility complex (MHC) has become the focus of diverse disciplines such as immunology, evolutionary ecology, and molecular evolution. MHC evolution has been studied extensively in diverse vertebrate lineages over the last few decades, and it has been suggested that birds differ from the established mammalian norm. Mammalian MHC genes evolve independently, and duplication history (i.e., orthology) can usually be traced back within lineages. In birds, this has been observed in only 3 pairs of closely related species. Here we report strong evidence for the persistence of orthology of MHC genes throughout an entire avian order. Phylogenetic reconstructions of MHC class II B genes in 14 species of owls trace back orthology over tens of thousands of years in exon 3. Moreover, exon 2 sequences from several species show closer relationships than sequences within species, resembling transspecies evolution typically observed in mammals. Thus, although previous studies suggested that long-term evolutionary dynamics of the avian MHC was characterized by high rates of concerted evolution, resulting in rapid masking of orthology, our results question the generality of this conclusion. The owl MHC thus opens new perspectives for a more comprehensive understanding of avian MHC evolution.

Identification and gene expression profiling of the Pum1 and Pum2 members of the Pumilio family in the chicken.

Members of the Pumilio (Pum) family of RNA-binding proteins act as translational repressors and are required for germ cell development and asymmetric division. We identified the chicken Pum1 and Pum2 genes and analyzed their expression patterns in various tissues. Comparative sequence analysis of the Pum1 and Pum2 proteins from the drosophila, chicken, mouse, and human revealed a high degree of evolutionary conservation in terms of the levels of homology of the peptide sequences and the structure of Pumilio homology domain (PUM-HD), C-terminal RNA-binding domain, with similar spacing between the adjacent Pum eight tandem repeats. In addition, phylogenetic patterns of pumilio family showed that Pum 1 and 2 of chicken are more closely related to those of mouse and human than other species and Pum1 is more conserved than Pum2. Using real-time RT-PCR, the expression levels of the Pum1 and Pum2 genes were found to be highest in hatched female gonads, and high-level expression of Pum2 was detected in 12-day and hatched gonads among the various chicken embryonic tissues tested. In adult tissues, the expression levels of Pum1 and Pum2 were expressed at higher levels in the testis and muscle than in any other tissue. The characteristics of the tissue-specific expression of Pum genes suggest that Pum1 and Pum2 have effects crucially in particular stage during development of chicken gonads depending on sexual maturation.