Gene duplication and adaptive evolution of Toll-like receptor genes in birds

Expression profile of Toll-like receptors and cytokines in the cecal tonsil of chickens challenged with Eimeria tenella

Chicken coccidiosis, caused by Eimeria spp., seriously affects the development of the poultry breeding industry. Currently, extensive studies of chicken coccidiosis are mostly focused on acquired immune responses, while information about the innate immune response of chicken coccidiosis is lacking. Toll-like receptor (TLR), the key molecule of the innate immune response, connects innate and adaptive immune responses and induces an immune response against various pathogen infections. Therefore, the quantitative real-time PCR was used to characterize the expression profile of chicken TLRs (chTLRs) and associated cytokines in the cecal tonsil of chickens infected with Eimeria tenella. The results showed that the expression of chTLR1a, chTLR2a, and chTLR5 was significantly upregulated at 3 h post-infection, while chTLR1b, chTLR2b, chTLR3, chTLR7, chTLR15 and chTLR21 was significantly downregulated (p < 0.05). In addition, chTLR1a expression rapidly reached the peaked expression at 3 h post-infection, while chTLR2b and chTLR15 peaked at 168 h post-infection, and chTLR2a expression was highest among chTLRs, peaking at 48 h post-infection (p < 0.05). For cytokines, interleukin (IL)-6 and tumor necrosis factor (TNF)-α peaked at 96 h post-infection, IL-4 and IL-12 peaked at 144 h post-infection, and interferon-γ expression was highest among cytokines at 120 h post-infection. In addition, IL-12 and IL-17 were markedly upregulated at 6 h post-infection (p < 0.05). These results provide insight into innate immune molecules during E. tenella infection in chickens and suggest that innate immune responses may mediate resistance to chicken coccidiosis.

Adaptive evolution of pancreatic ribonuclease gene (RNase1) in Cetartiodactyla

Pancreatic ribonuclease (RNase1), a digestive enzyme produced by the pancreas, is associated with the functional adaptation of dietary habits and is regarded as an attractive model system for studies of molecular evolution. In this study, we identified 218 functional genes and 48 pseudogenes from 114 species that span all four Cetartiodactyla lineages: two herbivorous lineages (Ruminantia and Tylopoda) and two non-herbivorous lineages (Cetancodonta and Suoidea). Multiple RNase1 genes were detected in all species of the two herbivorous lineages, and phylogenetic and genomic location analyses demonstrated that independent gene duplication events occurred in Ruminantia and Tylopoda. In Ruminantia, the gene duplication events occurred in the ancestral branches of the lineage in the Middle Eocene, a time of increasing climatic seasonality during which Ruminantia rapidly radiated. In contrast, only a single RNase1 gene was observed in the species of the two non-herbivorous lineages (Cetancodonta and Suoidea), suggesting that the previous Cetacea-specific loss hypothesis should be rejected. Moreover, the duplicated genes of RNase1 in the two herbivorous lineages (Ruminantia and Tylopoda) may have undergone functional divergence. In combination with the temporal coincidence between gene replication and the enhanced climatic seasonality during the Middle Eocene, this functional divergence suggests that RNase1 gene duplication was beneficial for Ruminantia to use the limited quantities of sparse fibrous vegetation and adapt to seasonal changes in climate. In summary, the findings indicate a complex and intriguing evolutionary pattern of RNase1 in Cetartiodactyla and demonstrate the molecular mechanisms by which organisms adapt to the environment.

A nuclear genome assembly of an extinct flightless bird, the little bush moa

We present a draft genome of the little bush moa (Anomalopteryx didiformis)-one of approximately nine species of extinct flightless birds from Aotearoa, New Zealand-using ancient DNA recovered from a fossil bone from the South Island. We recover a complete mitochondrial genome at 249.9× depth of coverage and almost 900 megabases of a male moa nuclear genome at ~4 to 5× coverage, with sequence contiguity sufficient to identify more than 85% of avian universal single-copy orthologs. We describe a diverse landscape of transposable elements and satellite repeats, estimate a long-term effective population size of ~240,000, identify a diverse suite of olfactory receptor genes and an opsin repertoire with sensitivity in the ultraviolet range, show that the wingless moa phenotype is likely not attributable to gene loss or pseudogenization, and identify potential function-altering coding sequence variants in moa that could be synthesized for future functional assays. This genomic resource should support further studies of avian evolution and morphological divergence.

Expression of Toll-like receptors and host defence peptides in the cecum of chicken challenged with Eimeria tenella

Chicken coccidiosis, caused by Eimeria protozoa, affects poultry farming. Toll-like receptors (TLRs) and host defence peptides (HDPs) help host innate immune responses to eliminate invading pathogens, but their roles in Eimeria tenella infection remain poorly understood. Herein, 14-day-old chickens were treated orally with 50,000 E. tenella oocysts and the cecum was dissected at different timepoints. mRNA expression of 10 chicken TLRs (chTLRs) and five HDPs was measured by quantitative real-time PCR. chTLR7 and chTLR15 were upregulated significantly at 3 h post-infection while other chTLRs were downregulated (p < .05). chTLR1a, chTLR1b, chTLR2b and chTLR4 peaked at 36 h post-infection, chTLR3, chTLR5 and chTLR15 peaked at 72 h post-infection and chTLR21 expression was highest among chTLRs, peaking at 48 h post-infection (p < 0.05). For HDPs, cathelicidin (CATH) 1 to 3 and B1 peaked at 48 h post-infection, liver-expressed antimicrobial peptide 2 peaked at 96 h post-infection, and CATH 2 expression was highest among HDPs. CATH2 and CATH3 were markedly upregulated at 3 h post-infection (p < .05). The results provide insight into innate immune molecules during E. tenella infection in chicken, and indicate that innate immune responses may mediate resistance to chicken coccidiosis.

Genetic variation at innate and adaptive immune genes - contrasting patterns of differentiation and local adaptation in a wild gull

Immunogenetic variation in natural vertebrate populations is expected to respond to spatial and temporal fluctuations in pathogen assemblages. While spatial heterogeneity in pathogen-driven selection enhances local immunogenetic adaptations and population divergence, different immune genes may yield contrasting responses to the environment. Here, we investigated population differentiation at the key pathogen recognition genes of the innate and adaptive immune system in a colonial bird species, the black-headed gull Chroicocephalus ridibundus. We assessed genetic variation at three toll-like receptor (TLR) genes (innate immunity) and the major histocompatibility complex (MHC) class I and II genes (adaptive immunity) in gulls from seven colonies scattered across Poland. As expected, we found much greater polymorphism at the MHC than TLRs. Population differentiation at the MHC class II, but not MHC-I, was significantly stronger than at neutral microsatellite loci, suggesting local adaptation. This could reflect spatial variation in the composition of extracellular parasite communities (e.g., helminths), possibly driven by sharp differences in habitat structure between colonies. Despite contrasting patterns of population differentiation, both MHC classes showed similar regimes of diversifying selection. Some significant population differentiation was also observed at TLRs, suggesting that innate immune receptors may respond to fine-scale spatial variation in pathogen pressure, although this pattern could have been enhanced by drift. Our results suggested that local adaptation at the pathogen recognition immune genes can be maintained at relatively small or moderate spatial scales in species with high dispersal potential and they highlighted the complexity of immunogenetic responses of animals to heterogeneous environments.

Evolutionary trade-off between innate and acquired immune defences in birds

BACKGROUND

The development, maintenance, and use of immune defences are costly. Therefore, animals face trade-offs in terms of resource allocation within their immune system and between their immune system and other physiological processes. To maximize fitness, evolution may favour investment in one immunological defence or subsystem over another in a way that matches a species broader life history strategy. Here, we used phylogenetically-informed comparative analyses to test for relationships between two immunological components. Natural antibodies and complement were used as proxies for the innate branch; structural complexity of the major histocompatibility complex (MHC) region was used for the acquired branch.

RESULTS

We found a negative association between the levels of natural antibodies (i.e., haemagglutination titre) and the total MHC gene copy number across the avian phylogeny, both at the species and family level. The family-level analysis indicated that this association was apparent for both MHC-I and MHC-II, when copy numbers within these two MHC regions were analysed separately. The association remained significant after controlling for basic life history components and for ecological traits commonly linked to pathogen exposure.

CONCLUSION

Our results provide the first phylogenetically robust evidence for an evolutionary trade-off within the avian immune system, with a more developed acquired immune system (i.e., more complex MHC architecture) in more derived bird lineages (e.g., passerines) being accompanied by an apparent downregulation of the innate immune system.

Divergent evolution drives high diversity of toll-like receptors (TLRs) in passerine birds: Buntings and finches

Toll-like receptors (TLRs) form a key component of animal innate immunity, being responsible for recognition of conserved microbial structures. As such, TLRs may be subject to diversifying and balancing selection, which maintains allelic variation both within and between populations. However, most research on TLRs in non-model avian species is focused on bottlenecked populations with depleted genetic variation. Here, we assessed variation at the extracellular domains of three TLR genes (TLR1LA, TLR3, TLR4) across eleven species from two passerine families of buntings (Emberizidae) and finches (Fringillidae), all having large breeding population sizes (millions of individuals). We found extraordinary TLR polymorphism in our study taxa, with >100 alleles detected at TLR1LA and TLR4 across species and high haplotype diversity (>0.75) in several species. Despite recent species divergence, no nucleotide allelic variants were shared between species, suggesting rapid TLR evolution. Higher variation at TLR1LA and TLR4 than TLR3 was associated with a stronger signal of diversifying selection, as measured with nucleotide substitutions rates and the number of positively selected sites (PSS). Structural protein modelling of TLRs showed that some PSS detected within TLR1LA and TLR4 were previously recognized as functionally important sites or were located in their proximity, possibly affecting ligand recognition. Furthermore, we identified PSS responsible for major surface electrostatic charge clustering, which may indicate their adaptive importance. Our study provides compelling evidence for the divergent evolution of TLR genes in buntings and finches and indicates that high TLR variation may be adaptively maintained via diversifying selection acting on functional ligand binding sites.

Regulation of mRNA and miRNA in the response to Salmonella enterica serovar Enteritidis infection in chicken cecum

BACKGROUND

Salmonella enterica, serovar Enteritidis (SE) is a food-borne pathogen, which can cause great threat to human health through consumption of the contaminated poultry products. Chicken is the main host of SE. The mRNA and microRNA (miRNA) expression profiles were analyzed on cecum of Shouguang chicken via next-generation sequencing and bioinformatics approaches. The treated group was inoculated SE, and the control group was inoculated with phosphate buffer saline (PBS).

RESULTS

There were 1760 differentially expressed mRNAs in the SE-infected group, of which 1046 were up-regulated mRNA, and 714 were down-regulated mRNA. In addition, a total of 821 miRNAs were identified, and 174 miRNAs were differentially expressed, of which 100 were up-regulated and 74 were down-regulated. Functional enrichment of differentially expressed mRNAs was similar to miRNA target genes. The functional analysis results of differentially expressed mRNAs and miRNAs were performed. Immune-related processes and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways were enriched by up-regulated mRNA. The down-regulated mRNAs were enriched in tissue development and metabolic-related KEGG pathways. The functional analysis of up-regulated miRNA target genes was similar to the down-regulated mRNAs. The down-regulated miRNA target genes were enriched in metabolic-related GO (Gene Ontology) -BP (Biological process) terms and KEGG pathways. The overlap of the up-regulated mRNA and the up-regulated miRNA target genes (class I) was 325, and the overlap of the down-regulated miRNA target genes (class II) was 169. The class I enriched in the immune-related GO-BP terms and KEGG pathways. The class II mainly enriched in metabolic-related GO-BP terms and KEGG pathways. Then we detected the expression of mRNA and miRNA through qRT-PCR. The results shown that the expression of HHIP, PGM1, HTR2B, ITGB5, RELN, SFRP1, TCF7L2, SCNN1A, NEK7, miR-20b-5p, miR-1662, miR-15a, miR-16-1-3p was significantly different between two groups. Dual-luciferase reporter assay was used to detect the relationship between miR-20b-5p and SCNN1A. The result indicated that miR-20b-5p regulate immune or metabolic responses after SE infection in Shouguang chickens by directly targeting SCNN1A.

CONCLUSIONS

The findings here contribute to the further analysis of the mechanism of mRNA and miRNA defense against SE infection, and provide a theoretical foundation for the molecular disease-resistant breeding of chickens.

Comparative genomics of the waterfowl innate immune system

Animal species differ considerably in their ability to fight off infections. Finding the genetic basis of these differences is not easy, as the immune response is comprised of a complex network of proteins that interact with one another to defend the body against infection. Here, we used population- and comparative genomics to study the evolutionary forces acting on the innate immune system in natural hosts of the avian influenza virus (AIV). For this purpose, we used a combination of hybrid capture, next- generation sequencing and published genomes to examine genetic diversity, divergence, and signatures of selection in 127 innate immune genes at a micro- and macroevolutionary time scale in 26 species of waterfowl. We show across multiple immune pathways (AIV-, toll-like-, and RIG-I -like receptors signalling pathways) that genes involved genes in pathogen detection (i.e., toll-like receptors) and direct pathogen inhibition (i.e., antimicrobial peptides and interferon-stimulated genes), as well as host proteins targeted by viral antagonist proteins (i.e., mitochondrial antiviral-signaling protein, [MAVS]) are more likely to be polymorphic, genetically divergent, and under positive selection than other innate immune genes. Our results demonstrate that selective forces vary across innate immune signaling signalling pathways in waterfowl, and we present candidate genes that may contribute to differences in susceptibility and resistance to infectious diseases in wild birds, and that may be manipulated by viruses. Our findings improve our understanding of the interplay between host genetics and pathogens, and offer the opportunity for new insights into pathogenesis and potential drug targets.

Evolution of RNA sensing receptors in birds

Birds are important hosts for many RNA viruses, including influenza A virus, Newcastle disease virus, West Nile virus and coronaviruses. Innate defense against RNA viruses in birds involves detection of viral RNA by pattern recognition receptors. Several receptors of different classes are involved, such as endosomal toll-like receptors and cytoplasmic retinoic acid-inducible gene I-like receptors, and their downstream adaptor proteins. The function of these receptors and their antagonism by viruses is well established in mammals; however, this has received less attention in birds. These receptors have been characterized in a few bird species, and the completion of avian genomes will permit study of their evolution. For each receptor, functional work has established ligand specificity and activation by viral infection. Engagement of adaptors, regulation by modulators and the supramolecular organization of proteins required for activation are incompletely understood in both mammals and birds. These receptors bind conserved nucleic acid agonists such as single- or double-stranded RNA and generally show purifying selection, particularly the ligand binding regions. However, in birds, these receptors and adaptors differ between species, and between individuals, suggesting that they are under selection for diversification over time. Avian receptors and signalling pathways, like their mammalian counterparts, are targets for antagonism by a variety of viruses, intent on escape from innate immune responses.

Evolution of innate and adaptive immune genes in a non-model waterbird, the common tern

Toll-like receptors (TLRs) and the Major Histocompatibility Complex (MHC) are the key pathogen-recognition genes of vertebrate immune system and they have a crucial role in the initiation of innate and adaptive immune response, respectively. Recent advancements in sequencing technology sparked research on highly duplicated MHC genes in non-model species, but TLR variation in natural vertebrate populations has remained little studied and comparisons of polymorphism across both TLRs and MHC are scarce. Here, we aimed to compare variation across innate (four TLR loci) and adaptive (MHC class I and class II) immune genes in a non-model avian species, the common tern Sterna hirundo. We detected relatively high allelic richness at TLR genes (9-48 alleles per locus), which was similar to or even higher than the estimated per locus allelic richness at the MHC (24-30 alleles at class I and 13-16 alleles at class II under uniform sample sizes). Despite this, the total number of MHC alleles across all duplicated loci (four class I and three class II) was much higher and MHC alleles showed greater sequence divergence than TLRs. Positive selection targeted relatively more sites at the MHC than TLRs, but the strength of selection (dN/dS ratios) at TLRs was higher when compared to MHC class I. There were also differences in the signature of positive selection and recombination (gene conversion) between MHC class I and II (stronger signature at class II), suggesting that mechanisms maintaining variation at the MHC may vary between both classes. Our study indicates that allelic richness of both innate and adaptive immune receptors may be maintained at relatively high levels in viable avian populations and we recommend a transition from the traditional gene-specific to multi-gene approach in studying molecular evolution of vertebrate immune system.