Positive diversifying selection in avian Mx genes

Myxovirus resistance (Mx) Gene Diversity in Avian Influenza Virus Infections

Avian influenza viruses (AIVs) pose threats to animal and human health. Outbreaks from the highly pathogenic avian influenza virus (HPAIV) in indigenous chickens in Bangladesh are infrequent. This could be attributed to the Myxovirus resistance (Mx) gene. To determine the impact of Mx gene diversity on AIV infections in chicken, we assessed the Mx genes, AIVs, and anti-AIV antibodies. DNA from blood cells, serum, and cloacal swab samples was isolated from non-vaccinated indigenous chickens and vaccinated commercial chickens. Possible relationships were assessed using the general linear model (GLM) procedure. Three genotypes of the Mx gene were detected (the resistant AA type, the sensitive GG type, and the heterozygous AG type). The AA genotype (0.48) was more prevalent than the GG (0.19) and the AG (0.33) genotypes. The AA genotype was more prevalent in indigenous than in commercial chickens. A total of 17 hemagglutinating viruses were isolated from the 512 swab samples. AIVs were detected in two samples (2/512; 0.39%) and subtyped as H1N1, whereas Newcastle disease virus (NDV) was detected in the remaining samples. The viral infections did not lead to apparent symptoms. Anti-AIV antibodies were detected in 44.92% of the samples with levels ranging from 27.37% to 67.65% in indigenous chickens and from 26% to 87.5% in commercial chickens. The anti-AIV antibody was detected in 40.16%, 65.98%, and 39.77% of chickens with resistant, sensitive, and heterozygous genotypes, respectively. The genotypes showed significant association (p < 0.001) with the anti-AIV antibodies. The low AIV isolation rates and high antibody prevalence rates could indicate seroconversion resulting from exposure to the virus as it circulates. Results indicate that the resistant genotype of the Mx gene might not offer anti-AIV protection for chickens.

Four Mx Genes Identified in Andrias davidianus and Characterization of Their Response to Chinese Giant Salamander Iridovirus Infection

Amphibians, including Andrias davidianus, are declining worldwide partly due to infectious diseases. The Myxovirus resistance (Mx) gene is a typical interferon (IFN)-stimulated gene (ISG) involved in the antiviral immunity. Therefore, knowledge regarding the antiviral immunity of A. davidianus can be used for improved reproduction in captivity and protection in the wild. In this study, we amplified and characterized four different A. davidianus Mx genes (adMx) and generated temporal mRNA expression profiles in healthy and Chinese giant salamander iridovirus (GSIV) infected A. davidianus by qualitative real-time PCR (qPCR). The four adMx genes ranged in length from 2008 to 2840 bp. The sequences revealed conserved protein domains including the dynamin superfamily signature motif and the tripartite guanosine-5-triphosphate (GTP)-binding motif. Gene and deduced amino acid sequence alignment revealed relatively high sequence identity with the Mx genes and proteins of other vertebrates. In phylogenetic analysis, the adMx genes clustered together, but also clustered closely with those of fish species. The four adMx genes were broadly expressed in healthy A. davidianus, but were differentially expressed in the spleen during the GSIV infection. Our results show that the adMx genes share major structural features with their homologs, suggesting similar functions to those in other species.

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.

Cloning, characterization and expression of GTPase effecter domain of chicken Mx1 gene

.

Exosomal miRNA profiling from H5N1 avian influenza virus-infected chickens

Exosomes are membrane vesicles containing proteins, lipids, DNA, mRNA, and micro RNA (miRNA). Exosomal miRNA from donor cells can regulate the gene expression of recipient cells. Here, Ri chickens were divided into resistant (Mx/A; BF2/B21) and susceptible (Mx/G; BF2/B13) trait by genotyping of Mx and BF2 genes. Then, Ri chickens were infected with H5N1, a highly pathogenic avian influenza virus (HPAIV). Exosomes were purified from blood serum of resistant chickens for small RNA sequencing. Sequencing data were analysed using FastQCv0.11.7, Cutadapt 1.16, miRBase v21, non-coding RNA database, RNAcentral 10.0, and miRDeep2. Differentially expressed miRNAs were determined using statistical methods, including fold-change, exactTest using edgeR, and hierarchical clustering. Target genes were predicted using miRDB. Gene ontology analysis was performed using gProfiler. Twenty miRNAs showed significantly different expression patterns between resistant control and infected chickens. Nine miRNAs were up-regulated and 11 miRNAs were down-regulated in the infected chickens compared with that in the control chickens. In target gene analysis, various immune-related genes, such as cytokines, chemokines, and signalling molecules, were detected. In particular, mitogen-activated protein kinase (MAPK) pathway molecules were highly controlled by differentially expressed miRNAs. The result of qRT-PCR for miRNAs was identical with sequencing data and miRNA expression level was higher in resistant than susceptible chickens. This study will help to better understand the host immune response, particularly exosomal miRNA expression against HPAIV H5N1 and could help to determine biomarkers for disease resistance.

Characterization of Myxovirus resistance protein in birds showing different susceptibilities to highly pathogenic influenza virus

We compared the Mx expression and anti-viral function and the 3D structure of Mx protein in four species: chicken (Gallus gallus), whooper swan (Cygnus cygnus), jungle crow (Corvus macrorhynchos), and rock dove (Columba livia). We observed different mortalities associated with highly pathogenic avian influenza virus (HPAIV) infection to understand the relationship between Mx function as an immune response factor and HPAIV proliferation in bird cells. Different levels of Mx were observed among the different bird species after virus infection. Strong Mx expression was confirmed in the rock dove and whooper swan 6 hr after viral infection. The lowest virus copy numbers were observed in rock dove. The virus infectivity was significantly reduced in the BALB/3T3 cells expressing rock dove and jungle crow Mx. These results suggested that high Mx expression and significant Mx-induced anti-viral effects might result in the rock dove primary cells having the lowest virus copy number. Comparison of the expected 3D structure of Mx protein in all four bird species demonstrated that the structure of loop L4 varied among the investigated species. It was reported that differences in amino acid sequence in loop L4 affect antiviral activity in human and mouse cells, and a significant anti-viral effect was observed in the rock dove Mx. Thus, the amino acid sequence of loop L4 in rock dove might represent relatively high anti-viral activity.

Genetic variation of major histocompatibility complex (MHC) in wild Red Junglefowl (Gallus gallus)

The major histocompatibility complex (MHC) is a multi-family gene cluster that encodes proteins with immuno-responsive function. While studies of MHC in domesticated poultry are relatively common, very little is known about this highly polymorphic locus in wild Red Junglefowl (Gallus gallus), the natural progenitor of domestic chickens. We investigated the diversity of MHC within and among four wild Red Junglefowl populations across diversified natural habitats in South Central Vietnam. Based on a SNP panel of 84 sites spanning 210 Kb of the MHC-B locus, we identified 310 unique haplotypes in 398 chromosomes. None of these haplotypes have been described before and we did not observe any of the wild Red Junglefowl haplotypes in domesticated chickens. Analysis of molecular variance (AMOVA) revealed that 94.51% of observed haplotype variation was accounted for at the within individual level. Little genetic variance was apportioned within and among populations, the latter accounting only for 0.83%. We also found evidence of increased recombination, including numerous hotspots, and limited linkage disequilibrium among the 84 SNP sites. Compared to an average haplotype diversity of 3.55% among seventeen lines of domestic chickens, our results suggest extraordinarily high haplotype diversity remains in wild Red Junglefowl and is consistent with a pattern of balancing selection. Wild Red Junglefowl in Vietnam, therefore, represent a rich resource of natural genomic variation independent from artificial selection.

Genetic variation within the Mx gene of commercially selected chicken lines reveals multiple haplotypes, recombination and a protein under selection pressure

The Mx protein is one of the best-characterized interferon-stimulated antiviral mediators. Mx homologs have been identified in most vertebrates examined; however, their location within the cell, their level of activity, and the viruses they inhibit vary widely. Recent studies have demonstrated multiple Mx alleles in chickens and some reports have suggested a specific variant (S631N) within exon 14 confers antiviral activity. In the current study, the complete genome of nine elite egg-layer type lines were sequenced and multiple variants of the Mx gene identified. Within the coding region and upstream putative promoter region 36 SNP variants were identified, producing a total of 12 unique haplotypes. Each elite line contained from one to four haplotypes, with many of these haplotypes being found in only one line. Observation of changes in haplotype frequency over generations, as well as recombination, suggested some unknown selection pressure on the Mx gene. Trait association analysis with either individual SNP or haplotypes showed a significant effect of Mx haplotype on several egg production related traits, and on mortality following Marek's disease virus challenge in some lines. Examination of the location of the various SNP within the protein suggests synonymous SNP tend to be found within structural or enzymatic regions of the protein, while non-synonymous SNP are located in less well defined regions. The putative resistance variant N631 was found in five of the 12 haplotypes with an overall frequency of 47% across the nine lines. Two Mx recombinants were identified within the elite populations, indicating that novel variation can arise and be maintained within intensively selected lines. Collectively, these results suggest the conflicting reports in the literature describing the impact of the different SNP on chicken Mx function may be due to the varying context of haplotypes present in the populations studied.

Pathogen-driven adaptive evolution of myxovirus resistance (Mx) genes in fishes

Myxovirus resistance (Mx) proteins, which belong to the dynamin super-family, are known to inhibit RNA viral replication in a wide range of taxonomic groups, including fishes. Given their crucial role in host immune defense, the key amino acid residues in the GTP effector domain (GED) near the C-terminus are expected to evolve adaptively in order to protect the host against invading viral pathogens. The present study reveals the role of recombination and positive selection in the evolution of Mx proteins in fishes. While the GTP-binding domain in the N-terminal domain has experienced purifying selection, several amino acid residues in GED have evolved under positive selection, thus indicating adaptive evolution. Given the antiviral activity of GED, the adaptive evolutionary changes that were observed in this region are therefore predicted to be pathogen-driven.

Evolution-guided identification of antiviral specificity determinants in the broadly acting interferon-induced innate immunity factor MxA

MxA is an interferon-induced dynamin-like GTPase with wide-ranging antiviral activity, which hinges upon detection of unique viral structures that differ across virus families. Despite elucidation of its structure, the basis of MxA antiviral specificity remains enigmatic. We used an evolution-guided approach to identify the loop L4 of MxA as a hotspot for recurrent positive selection in primates. Further, we show that single amino acid changes in L4 are necessary and sufficient to explain dramatic differences in species-specific antiviral activity of primate MxA proteins against the orthomyxoviruses Thogoto virus and influenza A virus. Taken together, our findings identify a genetic determinant of MxA target recognition and suggest a model by which MxA achieves antiviral breadth without compromising viral specificity.

Insight into alternative approaches for control of avian influenza in poultry, with emphasis on highly pathogenic H5N1

Highly pathogenic avian influenza virus (HPAIV) of subtype H5N1 causes a devastating disease in poultry but when it accidentally infects humans it can cause death. Therefore, decrease the incidence of H5N1 in humans needs to focus on prevention and control of poultry infections. Conventional control strategies in poultry based on surveillance, stamping out, movement restriction and enforcement of biosecurity measures did not prevent the virus spreading, particularly in developing countries. Several challenges limit efficiency of the vaccines to prevent outbreaks of HPAIV H5N1 in endemic countries. Alternative and complementary approaches to reduce the current burden of H5N1 epidemics in poultry should be encouraged. The use of antiviral chemotherapy and natural compounds, avian-cytokines, RNA interference, genetic breeding and/or development of transgenic poultry warrant further evaluation as integrated intervention strategies for control of HPAIV H5N1 in poultry.

Drift and selection influence geographic variation at immune loci of prairie-chickens

Previous studies of immunity in wild populations have focused primarily on genes of the major histocompatibility complex (MHC); however, studies of model species have identified additional immune-related genes that also affect fitness. In this study, we sequenced five non-MHC immune genes in six greater prairie-chicken (Tympanuchus cupido) populations that have experienced varying degrees of genetic drift as a consequence of population bottlenecks and fragmentation. We compared patterns of geographic variation at the immune genes with six neutral microsatellite markers to investigate the relative effects of selection and genetic drift. Global F(ST) outlier tests identified positive selection on just one of five immune genes (IAP-1) in one population. In contrast, at other immune genes, standardized G'(ST) values were lower than those at microsatellites for a majority of pairwise population comparisons, consistent with balancing selection or with species-wide positive or purifying selection resulting in similar haplotype frequencies across populations. The effects of genetic drift were also evident as summary statistics (e.g., Tajima's D) did not differ from neutrality for the majority of cases, and immune gene diversity (number of haplotypes per gene) was correlated positively with population size. In summary, we found that both genetic drift and selection shaped variation at the five immune genes, and the strength and type of selection varied among genes. Our results caution that neutral forces, such as drift, can make it difficult to detect current selection on genes.

Association of Mx1 Asn631 variant alleles with reductions in morbidity, early mortality, viral shedding, and cytokine responses in chickens infected with a highly pathogenic avian influenza virus

Myxovirus-resistance (Mx) proteins are produced by host cells in response to type I interferons, and some members of the Mx gene family in mammals have been shown to limit replication of influenza and other viruses. According to an early report, chicken Mx1 variants encoding Asn at position 631 have antiviral activity, whereas variants with Ser at 631 lack activity in experiments evaluating Mx1 complementary DNA (cDNA) expressed ectopically in a cell line. We evaluated whether the Mx1 631 dimorphism influenced pathogenesis of highly pathogenic avian influenza virus (HPAIV) infection in chickens of two commercial broiler lines, each segregating for Asn631 and Ser631 variants. Following intranasal infection with HPAIV strain A/Chicken/Queretaro/14588-19/1995 H5N2, chickens homozygous for Asn631 allele were significantly more resistant to disease based on early mortality, morbidity, or virus shedding than Ser631 homozygotes. Higher amounts of splenic cytokine transcripts were observed in the Ser631 birds after infection, consistent with higher viral loads seen in this group and perhaps contributing to their higher morbidity. Nucleotide sequence determination of Mx1 cDNAs demonstrated that the Asn631 variants in the two chicken lines differed at several amino acid positions outside 631. In vitro experiments with a different influenza strain (low pathogenicity) failed to demonstrate an effect of Mx1 Asn631 on viral replication suggesting that in vivo responses may differ markedly from in vitro, or that choice of virus strain may be critical in demonstrating effects of chicken Mx1. Overall, these studies provide the first evidence that Mx1 has antiviral effects in chickens infected with influenza virus.

Evolutionary analysis and expression profiling of zebra finch immune genes

Genes of the immune system are generally considered to evolve rapidly due to host-parasite coevolution. They are therefore of great interest in evolutionary biology and molecular ecology. In this study, we manually annotated 144 avian immune genes from the zebra finch (Taeniopygia guttata) genome and conducted evolutionary analyses of these by comparing them with their orthologs in the chicken (Gallus gallus). Genes classified as immune receptors showed elevated d(N)/d(S) ratios compared with other classes of immune genes. Immune genes in general also appear to be evolving more rapidly than other genes, as inferred from a higher d(N)/d(S) ratio compared with the rest of the genome. Furthermore, ten genes (of 27) for which sequence data were available from at least three bird species showed evidence of positive selection acting on specific codons. From transcriptome data of eight different tissues, we found evidence for expression of 106 of the studied immune genes, with primary expression of most of these in bursa, blood, and spleen. These immune-related genes showed a more tissue-specific expression pattern than other genes in the zebra finch genome. Several of the avian immune genes investigated here provide strong candidates for in-depth studies of molecular adaptation in birds.

The cytoplasmic location of chicken mx is not the determining factor for its lack of antiviral activity

Background

Chicken Mx belongs to the Mx family of interferon-induced dynamin-like GTPases, which in some species possess potent antiviral properties. Conflicting data exist for the antiviral capability of chicken Mx. Reports of anti-influenza activity of alleles encoding an Asn631 polymorphism have not been supported by subsequent studies. The normal cytoplasmic localisation of chicken Mx may influence its antiviral capacity. Here we report further studies to determine the antiviral potential of chicken Mx against Newcastle disease virus (NDV), an economically important cytoplasmic RNA virus of chickens, and Thogoto virus, an orthomyxovirus known to be exquisitely sensitive to the cytoplasmic MxA protein from humans. We also report the consequences of re-locating chicken Mx to the nucleus.

Methodology/Principal Findings

Chicken Mx was tested in virus infection assays using NDV. Neither the Asn631 nor Ser631 Mx alleles (when transfected into 293T cells) showed inhibition of virus-directed gene expression when the cells were subsequently infected with NDV. Human MxA however did show significant inhibition of NDV-directed gene expression. Chicken Mx failed to inhibit a Thogoto virus (THOV) minireplicon system in which the cytoplasmic human MxA protein showed potent and specific inhibition. Relocalisation of chicken Mx to the nucleus was achieved by inserting the Simian Virus 40 large T antigen nuclear localisation sequence (SV40 NLS) at the N-terminus of chicken Mx. Nuclear re-localised chicken Mx did not inhibit influenza (A/PR/8/34) gene expression during virus infection in cell culture or influenza polymerase activity in A/PR/8/34 or A/Turkey/50-92/91 minireplicon systems.

Conclusions/Significance

The chicken Mx protein (Asn631) lacks inhibitory effects against THOV and NDV, and is unable to suppress influenza replication when artificially re-localised to the cell nucleus. Thus, the natural cytoplasmic localisation of the chicken Mx protein does not account for its lack of antiviral activity.

Mx gene diversity and influenza association among five wild dabbling duck species (Anas spp.) in Alaska

Mx (myxovirus-resistant) proteins are induced by interferon and inhibit viral replication as part of the innate immune response to viral infection in many vertebrates. Influenza A virus appears to be especially susceptible to Mx antiviral effects. We characterized exon 13 and the 3' UTR of the Mx gene in wild ducks, the natural reservoir of influenza virus and explored its potential relevance to influenza infection. We observed a wide range of intra- and interspecies variations. Total nucleotide diversity per site was 0.0014, 0.0027, 0.0044, 0.0051, and 0.0061 in mallards, northern shovelers, northern pintails, American wigeon, and American green-winged teals, respectively. There were 61 haplotypes present across all five species and four were shared among species. Additionally, we observed a significant association between Mx haplotype and influenza infection status in northern shovelers. However, we found no evidence of balancing or diversifying selection in this region of the Mx gene. Characterization of the duck Mx gene is an important step in understanding how the gene may affect disease resistance or susceptibility in wild populations. Furthermore, given that waterfowl act as a natural reservoir for influenza virus, the Mx gene could be an important determinant in the ecology of the virus.

PCR-RFLP genotyping protocol for chicken Mx gene G/A polymorphism associated with the S631N mutation

The Mx (myxovirus resistance) gene codes for a protein with antiviral activity. Non-synonymous G/A polymorphism at position 2032 of chicken Mx cDNA results in a change at amino acid 631 of the Mx protein. This mutation has been shown to affect the antiviral activity of the Mx molecule, although recent studies have not confirmed this effect in response to some influenza strains. Nevertheless, the G/A polymorphism could be important for the chicken's response to other viruses. A robust PCR-RFLP protocol for genotyping chicken Mx gene polymorphism associated with the S631N mutation was developed. The F primer anneals to the last intron of the Mx gene, and the R primer anneals to the last exon of the gene, with an expected PCR product of 299 bp. PCR products were digested with Hpy8I. This enzyme cuts the sequence 5'-GTN|NAC-3', 2 bp downstream of the Mx polymorphism for the G allele, whereas the fragment containing the A allele is not cleaved. One hundred and twenty-seven chickens (commercial broilers, White Leghorn and New Hampshire) were genotyped using this protocol, and genotyping data were validated by sequencing. Full identity of results between the two genotyping methods was observed for all 127 samples, proving the reliability and robustness of this PCR-RFLP protocol.

The differential evolutionary dynamics of avian cytokine and TLR gene classes

The potential for investigating immune gene diversity has been greatly enhanced by recent advances in sequencing power. In this study, variation at two categories of avian immune genes with differing functional roles, pathogen detection and mediation of immune mechanisms, was examined using high-throughput sequencing. TLRs identify and alert the immune system by detecting molecular motifs that are conserved among pathogenic microorganisms, whereas cytokines act as mediators of resulting inflammation and immunity. Nine genes from each class were resequenced in a panel of domestic chickens and wild jungle fowl (JF). Tests on population-wide genetic variation between the gene classes indicated that allele frequency spectra at each group were distinctive. TLRs showed evidence pointing toward directional selection, whereas cytokines had signals more suggestive of frequency-dependent selection. This difference persisted between the distributions considering only coding sites, suggesting functional relevance. The unique patterns of variation at each gene class may be constrained by their different functional roles in the immune response. TLRs identify a relatively limited number of exogeneous pathogenic-related patterns and would be required to adapt quickly in response to evolving novel microbes encountered in new environmental niches. In contrast, cytokines interact with many molecules in mediating the power of immune mechanisms, and accordingly respond to the selective stimuli of many infectious diseases. Analyses also indicated that a general pattern of high variability has been enhanced by widespread genetic exchange between chicken and red JF, and possibly between chicken and gray JF at TLR1LA and TLR2A.

Evidence of balanced diversity at the chicken interleukin 4 receptor alpha chain locus

BACKGROUND

The comparative analysis of genome sequences emerging for several avian species with the fully sequenced chicken genome enables the genome-wide investigation of selective processes in functionally important chicken genes. In particular, because of pathogenic challenges it is expected that genes involved in the chicken immune system are subject to particularly strong adaptive pressure. Signatures of selection detected by inter-species comparison may then be investigated at the population level in global chicken populations to highlight potentially relevant functional polymorphisms.

RESULTS

Comparative evolutionary analysis of chicken (Gallus gallus) and zebra finch (Taeniopygia guttata) genes identified interleukin 4 receptor alpha-chain (IL-4Ralpha), a key cytokine receptor as a candidate with a significant excess of substitutions at nonsynonymous sites, suggestive of adaptive evolution. Resequencing and detailed population genetic analysis of this gene in diverse village chickens from Asia and Africa, commercial broilers, and in outgroup species red jungle fowl (JF), grey JF, Ceylon JF, green JF, grey francolin and bamboo partridge, suggested elevated and balanced diversity across all populations at this gene, acting to preserve different high-frequency alleles at two nonsynonymous sites.

CONCLUSION

Haplotype networks indicate that red JF is the primary contributor of diversity at chicken IL-4Ralpha: the signature of variation observed here may be due to the effects of domestication, admixture and introgression, which produce high diversity. However, this gene is a key cytokine-binding receptor in the immune system, so balancing selection related to the host response to pathogens cannot be excluded.