Liming Ren

The duck genome and transcriptome provide insight into an avian influenza virus reservoir species

The duck (Anas platyrhynchos) is one of the principal natural hosts of influenza A viruses. We present the duck genome sequence and perform deep transcriptome analyses to investigate immune-related genes. Our data indicate that the duck possesses a contractive immune gene repertoire, as in chicken and zebra finch, and this repertoire has been shaped through lineage-specific duplications. We identify genes that are responsive to influenza A viruses using the lung transcriptomes of control ducks and ones that were infected with either a highly pathogenic (A/duck/Hubei/49/05) or a weakly pathogenic (A/goose/Hubei/65/05) H5N1 virus. Further, we show how the ducks defense mechanisms against influenza infection have been optimized through the diversification of its β-defensin and butyrophilin-like repertoires. These analyses, in combination with the genomic and transcriptomic data, provide a resource for characterizing the interaction between host and influenza viruses.

Expression of IgM, IgD, and IgY in a Reptile, Anolis carolinensis

The reptiles are the last major group of jawed vertebrates in which the organization of the IGH locus and its encoded Ig H chain isotypes have not been well characterized. In this study, we show that the green anole lizard (Anolis carolinensis) expresses three Ig H chain isotypes (IgM, IgD, and IgY) but no IgA. The presence of the δ gene in the lizard demonstrates an evolutionary continuity of IgD from fishes to mammals. Although the germline δ gene contains 11 CH exons, only the first 4 are used in the expressed IgD membrane-bound form. The μ chain lacks the cysteine in CH1 that forms a disulfide bond between H and L chains, suggesting that (as in IgM of some amphibians) the H and L polypeptide chains are not covalently associated. Although conventional IgM transcripts (four CH domains) encoding both secreted and membrane-bound forms were detected, alternatively spliced transcripts encoding a short membrane-bound form were also observed and shown to lack the first two CH domains (VDJ-CH3-CH4-transmembrane region). Similar to duck IgY, lizard IgY H chain (υ) transcripts encoding both full-length and truncated (IgYΔFc) forms (with two CH domains) were observed. The absence of an IgA-encoding gene in the lizard IGH locus suggests a complex evolutionary history for IgA in the saurian lineage leading to modern birds, lizards, and their relatives.

Ornithorhynchus anatinus (Platypus) Links the Evolution of Immunoglobulin Genes in Eutherian Mammals and Nonmammalian Tetrapods

The evolutionary origins of mammalian immunoglobulin H chain isotypes (IgM, IgD, IgG, IgE, and IgA) are still incompletely understood as these isotypes differ considerably in structure and number from their counterparts in nonmammalian tetrapods. We report in this study that the platypus (Ornithorhynchus anatinus) Ig H chain constant region gene locus contains eight Ig encoding genes, which are arranged in an μ-δ-ο-γ2-γ1-α1-ε-α2 order, spanning a total of ∼200 kb DNA, encoding six distinct isotypes. The ο (ο for Ornithorhynchus) gene encodes a novel Ig H chain isotype that consists of four constant region domains and a hinge, and is structurally different from any of the five known mammalian Ig classes. This gene is phylogenetically related to υ (ε) and γ, and thus appears to be a structural intermediate between these two genes. The platypus δ gene encodes ten heavy chain constant region domains, lacks a hinge region and is similar to IgD in amphibians and fish, but strikingly different from that in eutherian mammals. The platypus Ig H chain isotype repertoire thus shows a unique combination of genes that share similarity both to those of nonmammalian tetrapods and eutherian animals and demonstrates how phylogenetically informative species can be used to reconstruct the evolutionary history of functionally important genes.

Molecular evolution of the vertebrate TLR1 gene family - a complex history of gene duplication, gene conversion, positive selection and co-evolution

BackgroundThe Toll-like receptors represent a large superfamily of type I transmembrane glycoproteins, some common to a wide range of species and others are more restricted in their distribution. Most members of the Toll-like receptor superfamily have few paralogues; the exception is the TLR1 gene family with four closely related genes in mammals TLR1, TLR2, TLR6 and TLR10, and four in birds TLR1A, TLR1B, TLR2A and TLR2B. These genes were previously thought to have arisen by a series of independent gene duplications. To understand the evolutionary pattern of the TLR1 gene family in vertebrates further, we cloned the sequences of TLR1A, TLR1B, TLR2A and TLR2B in duck and turkey, constructed phylogenetic trees, predicted codons under positive selection and identified co-evolutionary amino acid pairs within the TLR1 gene family using sequences from 4 birds, 28 mammals, an amphibian and a fish.ResultsThis detailed phylogenetic analysis not only clarifies the gene gains and losses within the TLR1 gene family of birds and mammals, but also defines orthologues between these vertebrates. In mammals, we predict amino acid sites under positive selection in TLR1, TLR2 and TLR6 but not TLR10. We detect co-evolution between amino acid residues in TLR2 and the other members of this gene family predicted to maintain their ability to form functional heterodimers. In birds, we predict positive selection in the TLR2A and TLR2B genes at functionally significant amino acid residues. We demonstrate that the TLR1 gene family has mostly been subject to purifying selection but has also responded to directional selection at a few sites, possibly in response to pathogen challenge.ConclusionsOur phylogenetic and structural analyses of the vertebrate TLR1 family have clarified their evolutionary origins and predict amino acid residues likely to be important in the hosts defense against invading pathogens.

Extensive diversification of IgH subclass-encoding genes and IgM subclass switching in crocodilians

Crocodilians are a group of reptiles that are closely related to birds and are thought to possess a strong immune system. Here we report that the IgH locus in the Siamese crocodile and the Chinese alligator contains multiple μ genes, in contrast to other tetrapods. Both the μ2 and μ3 genes are expressed through class-switch recombination involving the switch region and germline transcription. Both IgM1 and IgM2 are present in the serum as polymers, which implies that IgM class switching may have significant roles in humoural immunity. The crocodilian α genes are the first IgA-encoding genes identified in reptiles, and these genes show an inverted transcriptional orientation similar to that of birds. The identification of both α and δ genes in crocodilians suggests that the IgH loci of modern living mammals, reptiles and birds share a common ancestral organization.

Difference in microRNA expression and editing profile of lung tissues from different pig breeds related to immune responses to HP-PRRSV.

Porcine reproductive and respiratory syndrome (PRRS) is one of the most devastating diseases for the pig industry. Our goal was to identify microRNAs involved in the host immune response to PRRS. We generated microRNA expression profiles of lung tissues from Tongcheng or Landrace pigs infected with a highly pathogenic PRRS virus (PRRSV) at 3, 5, 7 dpi (day post infection) and control individuals from these two breeds. Our data showed that 278 known and 294 novel microRNAs were expressed in these combined microRNA transcriptomes. Compared with control individuals, almost half of the known microRNAs (116 in Tongcheng and 153 in Landrace) showed significantly differential expression (DEmiRNAs) at least once. The numbers of down-regulated DEmiRNAs were larger than the corresponding number of up-regulated DEmiRNAs in both breeds. Interestingly, miR-2320-5p, which was predicted to bind to conserved sequences of the PRRSV genome, was down-regulated significantly at 3 dpi after PRRSV infection in both breeds. In addition, PRRSV infection induced a significant increase of microRNA editing level in both breeds. Our results provide novel insight into the role of microRNA in response to PRRSV infection in vivo, which will aid the research for developing novel therapies against PRRSV.

Analysis of the reptile CD1 genes: evolutionary implications

CD1, as the third family of antigen-presenting molecules, is previously only found in mammals and chickens, which suggests that the chicken and mammalian CD1 shared a common ancestral gene emerging at least 310 million years ago. Here, we describe CD1 genes in the green anole lizard and Crocodylia, demonstrating that CD1 is ubiquitous in mammals, birds, and reptiles. Although the reptilian CD1 protein structures are predicted to be similar to human CD1d and chicken CD1.1, CD1 isotypes are not found to be orthologous between mammals, birds, and reptiles according to phylogenetic analyses, suggesting an independent diversification of CD1 isotypes during the speciation of mammals, birds, and reptiles. In the green anole lizard, although the single CD1 locus and MHC I gene are located on the same chromosome, there is an approximately 10-Mb-long sequence in between, and interestingly, several genes flanking the CD1 locus belong to the MHC paralogous region on human chromosome 19. The CD1 genes in Crocodylia are located in two loci, respectively linked to the MHC region and MHC paralogous region (corresponding to the MHC paralogous region on chromosome 19). These results provide new insights for studying the origin and evolution of CD1.

Internal Duplications of DH, JH, and C Region Genes Create an Unusual IgH Gene Locus in Cattle

It has been suspected for many years that cattle possess two functional IgH gene loci, located on Bos taurus autosome (BTA) 21 and BTA11, respectively. In this study, based on fluorescence in situ hybridization and additional experiments, we showed that all functional bovine IgH genes were located on BTA21, and only a truncated μCH2 exon was present on BTA11. By sequencing of seven bacterial artificial chromosome clones screened from a Hostein cow bacterial artificial chromosome library, we generated a 678-kb continuous genomic sequence covering the bovine IGHV, IGHD, IGHJ, and IGHC genes, which are organized as IGHVn-IGHDn-IGHJn-IGHM1-(IGHDP-IGHV3-IGHDn)3-IGHJn-IGHM2-IGHD-IGHG3-IGHG1-IGHG2-IGHE-IGHA. Although both of two functional IGHM genes, IGHM1 and IGHM2, can be expressed via independent VDJ recombinations, the IGHM2 can also be expressed through class switch recombination. Likely because more IGHD segments can be involved in the expression of IGHM2, the IGHM2 gene was shown to be dominantly expressed in most tissues throughout different developmental stages. Based on the length and identity of the coding sequence, the 23 IGHD segments identified in the locus could be divided into nine subgroups (termed IGHD1 to IGHD9). Except two members of IGHD9 (14 nt in size), all other functional IGHD segments are longer than 30 nt, with the IGHD8 gene (149 bp) to be the longest. These remarkably long germline IGHD segments play a pivotal role in generating the exceptionally great H chain CDR 3 length variability in cattle.

Characterization of the MHC class II α-chain gene in ducks

In humans, classical MHC class II molecules include DQ, DR, and DP, which are similar in structure but consist of distinct α- and β-chains. The genes encoding these molecules are all located in the MHC class II gene region. In non-mammalian vertebrates such as chickens, only a single class II α-chain gene corresponding to the human DRA has been identified. Here, we report a characterization of the duck MHC class II α-chain (Anpl-DRA) encoding gene, which contains four exons encoding a typical signal peptide, a peptide-binding α1 domain, an immunoglobulin-like α2 domain, and Tm/Cyt, respectively. This gene is present in the duck genome as a single copy and is highly expressed in the spleen. Sequencing of cDNA and genomic DNA of the Anpl-DRA of different duck individuals/strains revealed low levels of genetic polymorphism, especially in the same strain, although most duck individuals have two different alleles. Otherwise, we found that the duck gene is located next to class II β genes, which is the same as in humans but different from the situation in chickens.

A comprehensive analysis of the germline and expressed TCR repertoire in White Peking duck

Recently, many immune-related genes have been extensively studied in ducks, but relatively little is known about their TCR genes. Here, we determined the germline and expressed repertoire of TCR genes in White Peking duck. The genomic organization of the duck TCRα/δ, TCRγ and unconventional TCRδ2 loci are highly conserved with their counterparts in mammals or chickens. By contrast, the duck TCRβ locus is organized in an unusual pattern, (Vβ)n-Dβ-(Jβ)2-Cβ1-(Jβ)4-Cβ2, which differs from the tandem-aligned clusters in mammals or the translocon organization in some teleosts. Excluding the first exon encoding the immunoglobulin domain, the subsequent exons of the two Cβ show significant diversity in nucleotide sequence and exon structure. Based on the nucleotide sequence identity, 49 Vα, 30 Vδ, 13 Vβ and 15 Vγ unique gene segments are classified into 3 Vα, 5 Vδ, 4 Vβ and 6 Vγ subgroups, respectively. Phylogenetic analyses revealed that most duck V subgroups, excluding Vβ1, Vγ5 and Vγ6, have closely related orthologues in chicken. The coding joints of all cDNA clones demonstrate conserved mechanisms that are used to increase junctional diversity. Collectively, these data provide insight into the evolution of TCRs in vertebrates and improve our understanding of the avian immune system.