Experimental viral evolution to specific host MHC genotypes reveals fitness and virulence trade-offs in alternative MHC types

Genetic Diversity and Sequence Conservation of Peptide-Binding Regions of MHC Class I Genes in Pig, Cattle, Chimpanzee, and Human

Comparative analyses of MHC gene diversity and evolution across different species could offer valuable insights into the evolution of MHC genes. Intra- and inter-species sequence diversity and conservation of 12 classical major histocompatibility complex (MHC) class I genes from cattle, chimpanzees, pigs, and humans was analyzed using 20 representative allelic groups for each gene. The combined analysis of paralogous loci for each species revealed that intra-locus amino-acid sequence variations in the peptide-binding region (PBR) of MHC I genes did not differ significantly between species, ranging from 8.44% for SLA to 10.75% for BoLA class I genes. In contrast, intraspecies differences in the non-PBRs of these paralogous genes were more pronounced, varying from 4.59% for SLA to 16.89% for HLA. Interestingly, the Shannon diversity index and rate of nonsynonymous substitutions for PBR were significantly higher in SLA and BoLA than those in Patr and HLA. Analysis of peptide-binding pockets across all analyzed MHC class I genes of the four species indicated that pockets A and E showed the lowest and highest diversity, respectively. The estimated divergence times suggest that primate and artiodactyl MHC class I genes diverged 60.41 Mya, and BoLA and SLA genes diverged 35.34 Mya. These results offer new insights into the conservation and diversity of MHC class I genes in various mammalian species.

Cross-population genetic analysis revealed genetic variation and selection in the Ovar-DRB1 gene of Indian sheep breeds

In sheep, MHC variability is studied widely to explore disease association. The aim of the current study was to explore the genetic diversity of Ovar-DRB diversity across sheep breeds of India. Here, Ovar-DRB1 locus was studied across 20 sheep breeds. DRB1 was amplified (301 bp) and sequenced using a PCR-sequence-based typing approach. Results revealed a high degree of heterozygosity across breeds (mean: 73.99%). Overall mean distance for DRB1 was highest in Sangamneri (0.18) and lowest in Madgyal sheep (0.10). There was a higher rate of transition, across breeds. Further, 39 alleles were isolated in different breeds, out of which 10 were new. To allow easy access and use of the immune-polymorphic database, an online database management system was launched (http://www.mhcdbms.in/). Nucleotide content across breeds for the DRB1 region revealed the richness of GC content (59.26%). Wu-Kabat index revealed vast genetic variation across peptide binding sites (PBS) of DRB1. Residues 6, 66, 69, 52, and 81, were polymorphic showing utility for antigen presentation. All breeds were under positive selection for DRB1 locus (dN > dS). Study revealed the importance of DRB locus diversity for beta chain specifically at PBS across sheep breeds of the Indian subcontinent and presented evidence of positive selection for DRB owing to its evolutionary significance.

The evolution of host specialization in an insect pathogen

Niche breadth coevolution between biotic partners underpins theories of diversity and co-existence and influences patterns of disease emergence and transmission in host-parasite systems. Despite these broad implications, we still do not fully understand how the breadth of parasites' infectivity evolves, the nature of any associated costs, or the genetic basis of specialization. Here, we serially passage a granulosis virus on multiple inbred populations of its Plodia interpunctella host to explore the dynamics and outcomes of specialization. In particular, we collect time series of phenotypic and genetic data to explore the dynamics of host genotype specialization throughout the course of experimental evolution and examine two fitness components. We find that the Plodia interpunctella granulosis virus consistently evolves and increases in overall specialization, but that our two fitness components evolve independently such that lines can specialize in productivity or infectivity. Furthermore, we find that specialization in our experiment is a highly polygenic trait best explained by a combination of evolutionary mechanisms. These results are important for understanding the evolution of specialization in host-parasite interactions and its broader implications for co-existence, diversification, and infectious disease management.

MHC class II genotype-by-pathogen genotype interaction for infection prevalence in a natural rodent-Borrelia system

MHC genes are extraordinarily polymorphic in most taxa. Host-pathogen coevolution driven by negative frequency-dependent selection (NFDS) is one of the main hypotheses for the maintenance of such immunogenetic variation. Here, we test a critical but rarely tested assumption of this hypothesis-that MHC alleles affect resistance/susceptibility to a pathogen in a strain-specific way, that is, there is a host genotype-by-pathogen genotype interaction. In a field study of bank voles naturally infected with the tick-transmitted bacterium Borrelia afzelii, we tested for MHC class II (DQB) genotype-by-B. afzelii strain interactions for infection prevalence between 10 DQB alleles and seven strains. One allele (DQB*37) showed an interaction, such that voles carrying DQB*37 had higher prevalence of two strains and lower prevalence of one strain than individuals without the allele. These findings were corroborated by analyses of strain composition of infections, which revealed an effect of DQB*37 in the form of lower β diversity among infections in voles carrying the allele. Taken together, these results provide rare support at the molecular genetic level for a key assumption of models of antagonistic coevolution through NFDS.

Long term patterns of association between MHC and helminth burdens in the bank vole support Red Queen dynamics

Major histocompatibility complex (MHC) genes encode proteins crucial for adaptive immunity of vertebrates. Negative frequency-dependent selection (NFDS), resulting from adaptation of parasites to common MHC types, has been hypothesized to maintain high, functionally relevant polymorphism of MHC, but demonstration of this relationship has remained elusive. In particular, differentiation of NFDS from fluctuating selection, resulting from changes in parasite communities in time and space (FS), has proved difficult in short-term studies. Here, we used temporal data, accumulated through long-term monitoring of helminths infecting bank voles (Myodes glareolus), to test specific predictions of NFDS on MHC class II. Data were collected in three, moderately genetically differentiated subpopulations in Poland, which were characterized by some stable spatiotemporal helminth communities but also events indicating introduction of new species and loss of others. We found a complex association between individual MHC diversity and species richness, where intermediate numbers of DRB supertypes correlated with lowest species richness, but the opposite was true for DQB supertypes - arguing against universal selection for immunogenetic optimality. We also showed that particular MHC supertypes explain a portion of the variance in prevalence and abundance of helminths, but this effect was subpopulation-specific, which is consistent with both NFDS and FS. Finally, in line with NFDS, we found that certain helminths that have recently colonized or spread in a given subpopulation, more frequently or intensely infected voles with MHC supertypes that have been common in the recent past. Overall, our results highlight complex spatial and temporal patterns of MHC-parasite associations, the latter being consistent with Red Queen coevolutionary dynamics.

Friend virus severity is associated with male mouse social status and environmental temperature

Pathogen virulence is highly variable within populations, and although many factors contributing to virulence differences are known, there is still much variation left unexplained. Identifying and characterizing environmental conditions associated with different virulence levels is therefore an important undertaking in infectious disease research. One factor considered to be a major determinant of overall health and susceptibility to disease in social animals is social status. Health differences associated with social status are thought to be caused by different levels of chronic stress in higher- versus lower-status individuals. There is considerable evidence that these effects extend to the standing immune profile and that social status directly influences susceptibility to pathogens. Here we examined the association between dominance status in male wild-derived house mice, Mus musculus, and susceptibility to Friend virus complex in the context of seminatural populations with intense male-male competition and no predation. Due to an interruption in our facility's heating system, we were unexpectedly presented with the opportunity to assess how reduced ambient temperature influences the association of host social status and pathogen virulence. Environmental temperature has been implicated as a contributor to pathogen virulence, giving us a unique chance to examine its role in a previously unexamined pathogen system, while the added context of social status can expand our understanding of how the interaction of different environmental conditions affects virulence. We found that pathogen virulence and replication were lower in socially dominant hosts compared to nondominant hosts. When temperature was reduced, cool enclosure-housed dominant males were more susceptible to infection than their warm enclosure-housed counterparts. The mechanistic underpinnings that link infectious disease and social status remain difficult to disentangle from their associated factors, but this study opens the door for future experiments using a novel approach in the most well-studied mammalian model available.

Coinfection with a virus constrains within-host infection load but increases transmission potential of a highly virulent fungal plant pathogen

The trade‐off between within‐host infection rate and transmission to new hosts is predicted to constrain pathogen evolution, and to maintain polymorphism in pathogen populations. Pathogen life‐history stages and their correlations that underpin infection development may change under coinfection with other parasites as they compete for the same limited host resources. Cross‐kingdom interactions are common among pathogens in both natural and cultivated systems, yet their impacts on disease ecology and evolution are rarely studied. The host plant Plantago lanceolata is naturally infected by both Phomopsis subordinaria, a seed killing fungus, as well as Plantago lanceolata latent virus (PlLV) in the Åland Islands, SW Finland. We performed an inoculation assay to test whether coinfection with PlLV affects performance of two P. subordinaria strains, and the correlation between within‐host infection rate and transmission potential. The strains differed in the measured life‐history traits and their correlations. Moreover, we found that under virus coinfection, within‐host infection rate of P. subordinaria was smaller but transmission potential was higher compared to strains under single infection. The negative correlation between within‐host infection rate and transmission potential detected under single infection became positive under coinfection with PlLV. To understand whether within‐host and between‐host dynamics are correlated in wild populations, we surveyed 260 natural populations of P. lanceolata for P. subordinaria infection occurrence. When infections were found, we estimated between‐hosts dynamics by determining pathogen population size as the proportion of infected individuals, and within‐host dynamics by counting the proportion of infected flower stalks in 10 infected plants. In wild populations, the proportion of infected flower stalks was positively associated with pathogen population size. Jointly, our results suggest that the trade‐off between within‐host infection load and transmission may be strain specific, and that the pathogen life‐history that underpin epidemics may change depending on the diversity of infection, generating variation in disease dynamics.

Optimal immune specificity at the intersection of host life history and parasite epidemiology

Hosts diverge widely in how, and how well, they defend themselves against infection and immunopathology. Why are hosts so heterogeneous? Both epidemiology and life history are commonly hypothesized to influence host immune strategy, but the relationship between immune strategy and each factor has commonly been investigated in isolation. Here, we show that interactions between life history and epidemiology are crucial for determining optimal immune specificity and sensitivity. We propose a demographically-structured population dynamics model, in which we explore sensitivity and specificity of immune responses when epidemiological risks vary with age. We find that variation in life history traits associated with both reproduction and longevity alters optimal immune strategies-but the magnitude and sometimes even direction of these effects depends on how epidemiological risks vary across life. An especially compelling example that explains previously-puzzling empirical observations is that depending on whether infection risk declines or rises at reproductive maturity, later reproductive maturity can select for either greater or lower immune specificity, potentially illustrating why studies of lifespan and immune variation across taxa have been inconclusive. Thus, the sign of selection on the life history-immune specificity relationship can be reversed in different epidemiological contexts. Drawing on published life history data from a variety of chordate taxa, we generate testable predictions for this facet of the optimal immune strategy. Our results shed light on the causes of the heterogeneity found in immune defenses both within and among species and the ultimate variability of the relationship between life history and immune specificity.

Comparative Analysis of SLA-1, SLA-2, and DQB1 Genetic Diversity in Locally-Adapted Kenyan Pigs and Their Wild Relatives, Warthogs

Swine leukocyte antigen (SLA) plays a central role in controlling the immune response by discriminating self and foreign antigens and initiating an immune response. Studies on SLA polymorphism have demonstrated associations between SLA allelic variants, immune response, and disease resistance. The SLA polymorphism is due to host-pathogen co-evolution resulting in improved adaptation to diverse environments making SLA a crucial genomic region for comparative diversity studies. Although locally-adapted African pigs have small body sizes, they possess increased resilience under harsh environmental conditions and robust immune systems with reported tolerance to some diseases, including African swine fever. However, data on the SLA diversity in these pigs are not available. We characterized the SLA of unrelated locally-adapted domestic pigs from Homa Bay, Kenya, alongside exotic pigs and warthogs. We undertook SLA comparative diversity of the functionally expressed SLA class I (SLA-1, SLA-2) and II (DQB1) repertoires in these three suids using the reverse transcription polymerase chain reaction (RT-PCR) sequence-based typing (SBT) method. Our data revealed higher genetic diversity in the locally-adapted pigs and warthogs compared to the exotic pigs. The nucleotide substitution rates were higher in the peptide-binding regions of the SLA-1, SLA-2, and DQB1 loci, indicative of adaptive evolution. We obtained high allele frequencies in the three SLA loci, including some breed-specific private alleles, which could guide breeders to increase their frequency through selection if confirmed to be associated with enhanced resilience. Our study contributes to the growing body of knowledge on genetic diversity in free-ranging animal populations in their natural environment, availing the first DQB1 gene data from locally-adapted Kenyan pigs.

Deep Sequencing of MHC-Adapted Viral Lines Reveals Complex Recombinational Exchanges With Endogenous Retroviruses Leading to High-Frequency Variants

Experimental evolution (serial passage) of Friend virus complex (FVC) in mice demonstrates phenotypic adaptation to specific host major histocompatibility complex (MHC) genotypes. These evolved viral lines show increased fitness and virulence in their host-genotype-of-passage, but display fitness and virulence tradeoffs when infecting unfamiliar host MHC genotypes. Here, we deep sequence these viral lines in an attempt to discover the genetic basis of FVC adaptation. The principal prediction for genotype-specific adaptation is that unique mutations would rise to high frequency in viral lines adapted to each host MHC genotype. This prediction was not supported by our sequencing data as most observed high-frequency variants were present in each of our independently evolved viral lines. However, using a multi-variate approach to measure divergence between viral populations, we show that populations of replicate evolved viral lines from the same MHC congenic mouse strain were more similar to one another than to lines derived from different MHC congenic mouse strains, suggesting that MHC genotype does predictably act on viral evolution in our model. Sequence analysis also revealed rampant recombination with endogenous murine leukemia virus sequences (EnMuLVs) that are encoded within the BALB/c mouse genome. The highest frequency variants in all six lines contained a 12 bp insertion from a recombinant EnMuLV source, suggesting such recombinants were either being favored by selection or were contained in a recombinational hotspot. Interestingly, they did not reach fixation, as if they are low fitness. The amount of background mutations linked to FVC/EnMuLV variable sites indicated that FVC/EnMuLV recombinants had not reached mutation selection equilibrium and thus, that EnMuLV sequences are likely continuously introgressing into the replicating viral population. These discoveries raise the question: is the expression of EnMuLV sequences in mouse splenocytes that permit recombination with exogenous FVC a pathogen or host adaptation?

Host genotype and genetic diversity shape the evolution of a novel bacterial infection

Pathogens continue to emerge from increased contact with novel host species. Whilst these hosts can represent distinct environments for pathogens, the impacts of host genetic background on how a pathogen evolves post-emergence are unclear. In a novel interaction, we experimentally evolved a pathogen (Staphylococcus aureus) in populations of wild nematodes (Caenorhabditis elegans) to test whether host genotype and genetic diversity affect pathogen evolution. After ten rounds of selection, we found that pathogen virulence evolved to vary across host genotypes, with differences in host metal ion acquisition detected as a possible driver of increased host exploitation. Diverse host populations selected for the highest levels of pathogen virulence, but infectivity was constrained, unlike in host monocultures. We hypothesise that population heterogeneity might pool together individuals that contribute disproportionately to the spread of infection or to enhanced virulence. The genomes of evolved populations were sequenced, and it was revealed that pathogens selected in distantly-related host genotypes diverged more than those in closely-related host genotypes. S. aureus nevertheless maintained a broad host range. Our study provides unique empirical insight into the evolutionary dynamics that could occur in other novel infections of wildlife and humans.

Pathogen-mediated selection favours the maintenance of innate immunity gene polymorphism in a widespread wild ungulate

Toll-like receptors (TLR) play a central role in recognition and host frontline defence against a wide range of pathogens. A number of recent studies have shown that TLR genes (Tlrs) often exhibit large polymorphism in natural populations. Yet, there is little knowledge on how this polymorphism is maintained and how it influences disease susceptibility in the wild. In previous work, we showed that some Tlrs exhibit similarly high levels of genetic diversity as genes of the Major Histocompatibility Complex (MHC), and signatures of contemporary balancing selection in roe deer (Capreolus capreolus), the most abundant cervid species in Europe. Here, we investigated the evolutionary mechanisms by which pathogen-mediated selection could shape this innate immunity genetic diversity by examining the relationships between Tlr (Tlr2, Tlr4 and Tlr5) genotypes (heterozygosity status and presence of specific alleles) and infections with Toxoplasma and Chlamydia, two widespread intracellular pathogens known to cause reproductive failure in ungulates. We showed that Toxoplasma and Chlamydia exposures vary significantly across years and landscape features with few co-infection events detected and that the two pathogens exert antagonistic selection on Tlr2 polymorphism. By contrast, we found limited support for Tlr heterozygote advantage. Our study confirmed the importance of looking beyond Mhc genes in wildlife immunogenetic studies. It also emphasized the necessity to consider multiple pathogen challenges and their spatiotemporal variation to improve our understanding of vertebrate defence evolution against pathogens.

Experimental virus evolution in cancer cell monolayers, spheroids, and tissue explants

Viral laboratory evolution has been used for different applications, such as modeling viral emergence, drug-resistance prediction, and therapeutic virus optimization. However, these studies have been mainly performed in cell monolayers, a highly simplified environment, raising concerns about their applicability and relevance. To address this, we compared the evolution of a model virus in monolayers, spheroids, and tissue explants. We performed this analysis in the context of cancer virotherapy by performing serial transfers of an oncolytic vesicular stomatitis virus (VSV-Δ51) in 4T1 mouse mammary tumor cells. We found that VSV-Δ51 gained fitness in each of these three culture systems, and that adaptation to the more complex environments (spheroids or explants) correlated with increased fitness in monolayers. Most evolved lines improved their ability to suppress β-interferon secretion compared to the VSV-Δ51 founder, suggesting that the selective pressure exerted by antiviral innate immunity was important in the three systems. However, system-specific patterns were also found. First, viruses evolved in monolayers remained more oncoselective that those evolved in spheroids, since the latter showed concomitant adaptation to non-tumoral mouse cells. Second, deep sequencing indicated that viral populations evolved in monolayers or explants tended to be more genetically diverse than those evolved in spheroids. Finally, we found highly variable outcomes among independent evolutionary lines propagated in explants. We conclude that experimental evolution in monolayers tends to be more reproducible than in spheroids or explants, and better preserves oncoselectivity. Our results also suggest that monolayers capture at least some relevant selective pressures present in more complex systems.

What evolutionary processes maintain MHC IIꞵ diversity within and among populations of stickleback?

Major Histocompatibility Complex (MHC) genes code for proteins that recognize foreign protein antigens to initiate T-cell-mediated adaptive immune responses. They are often the most polymorphic genes in vertebrate genomes. How evolution maintains this diversity remains of debate. Three main hypotheses seek to explain the maintenance of MHC diversity by invoking pathogen-mediated selection: heterozygote advantage, frequency-dependent selection, and fluctuating selection across landscapes or through time. Here, we use a large-scale field parasite survey in a stickleback metapopulation to test predictions derived from each of these hypotheses. We identify over 1000 MHC IIβ variants (alleles spanning paralogous genes) and find that many of them covary positively or negatively with parasite load, suggesting that these genes contribute to resistance or susceptibility. However, despite our large sample-size, we find no evidence for the widely cited stabilizing selection on MHC heterozygosity, in which individuals with an intermediate number of MHC variants have the lowest parasite burden. Nor do we observe a rare-variant advantage, or widespread fluctuating selection across populations. In contrast, we find that MHC diversity is best predicted by neutral genome-wide heterozygosity and between-population genomic divergence, suggesting neutral processes are important in shaping the pattern of metapopulation MHC diversity. Thus, although MHC IIβ is highly diverse and relevant to the type and intensity of macroparasite infection in these populations of stickleback, the main models of MHC evolution still provide little explanatory power in this system.

Allele-Specific Quantification of HLA-DRB1 Transcripts Reveals Imbalanced Allelic Expression That Modifies the Amino Acid Effects in HLA-DRβ1

OBJECTIVE

HLA association fine-mapping studies have shown the effects of missense variants in HLA-DRB1 on rheumatoid arthritis (RA) susceptibility, prognosis, and autoantibody production. However, the phenotypic effects of expression changes in HLA-DRB1 remain poorly understood. Therefore, we investigated the allele-specific expression of HLA-DRB1 and its effect on an HLA-DRβ1 structure-associated trait in RA.

METHODS

We quantified the allele-specific expression of each HLA-DRB1 3-field classic allele in 48 Korean RA patients with anti-citrullinated protein antibodies (ACPAs) and 319 healthy European subjects by using both RNA sequencing and HLA-DRB1 genotype data to calculate the relative expression strength of multiple HLA-DRB1 alleles (n = 14 in Koreans and n = 25 in Europeans) in each population. The known association between ACPA level and alanine at position 74 of HLA-DRβ1 in ACPA-positive RA was revisited to understand the phenotypic effect of allele-specific expression of HLA-DRB1 by modeling multivariate logistic regression with the genomic dosage or relative expression dosage of Ala-74 in 2 independent sets of 1,723 Korean RA patients with ACPA.

RESULTS

The relative expression strength was highly allele-specific, causing imbalanced allelic expression in HLA-DRB1 heterozygotes. The association between HLA-DRβ1 Ala-74 and ACPA level in RA was better explained by relative expression dosage of Ala-74 than by the genomic dosage (change in Akaike's information criterion = -6.98). Moreover, the expression variance of Ala-74 in Ala-74 heterozygotes with no genomic variance of Ala-74 was significantly associated with ACPA level (P = 2.26 × 10^-3 ).

CONCLUSION

Our findings illustrate the advantage of integrating quantitative and qualitative changes in HLA-DRB1 into a single model for understanding HLA-DRB1 associations.

Hidden paths to endless forms most wonderful: Complexity of bacterial motility shapes diversification of latent phenotypes

BACKGROUND

Evolution in one selective environment often latently generates phenotypic change that is manifested only later in different environments, but the complexity of behavior important to fitness in the original environment might influence the character of such latent-phenotype evolution. Using Myxococcus xanthus, a bacterium possessing two motility systems differing in effectiveness on hard vs. soft surfaces, we test (i) whether and how evolution while swarming on one surface-the selective surface-latently alters motility on the alternative surface type and (ii) whether patterns of such latent-phenotype evolution depend on the complexity of ancestral motility, specific ancestral motility genotypes and/or the selective surface of evolution. We analysze an experiment in which populations established from three ancestral genotypes-one with both motility systems intact and two others with one system debilitated-evolved while swarming across either hard or soft agar in six evolutionary treatments. We then compare motility-phenotype patterns across selective vs. alternative surface types.

RESULTS

Latent motility evolution was pervasive but varied in character as a function of the presence of one or two functional motility systems and, for some individual-treatment comparisons, the specific ancestral genotype and/or selective surface. Swarming rates on alternative vs. selective surfaces were positively correlated generally among populations with one functional motility system but not among those with two. This suggests that opportunities for pleiotropy and epistasis generated by increased genetic complexity underlying behavior can alter the character of latent-phenotype evolution. No tradeoff between motility performance across surface types was detected in the dual-system treatments, even after adaptation on a surface on which one motility system dominates strongly over the other in driving movement, but latent-phenotype evolution was instead idiosyncratic in these treatments. We further find that the magnitude of stochastic diversification at alternative-surface swarming among replicate populations greatly exceeded diversification of selective-surface swarming within some treatments and varied across treatments.

CONCLUSION

Collectively, our results suggest that increases in the genetic and mechanistic complexity of behavior can increase the complexity of latent-phenotype evolution outcomes and illustrate that diversification manifested during evolution in one environment can be augmented greatly by diversification of latent phenotypes manifested later.

Genetic variation and selection in the major histocompatibility complex Class II gene in the Guizhou pony

The Guizhou pony (GZP) is an indigenous species of equid found in the mountains of the Guizhou province in southwest China. We selected four regions of the equine leukocyte antigen (ELA), including DQA, DRA, DQB, and DRB, and used them to assess the diversity of the major histocompatibility complex (MHC) class II gene using direct sequencing technology. DRA had the lowest d_N/d_S ratio (0.560) compared with the other three loci, indicating that DRA was conserved and could be conserved after undergoing selective processes. Nine DQA, five DQB, nine DRA, and seven DRB codons were under significant positive selection at the antigen binding sites (ABS), suggesting that the selected residues in ABS may play a significant role in the innate immune system of the GZP. Two GZP alleles were shared with Przewalski’s horse, and six older GZP haplotypes had a better relationship with other horse species by one or two mutational steps, indicating that the GZP may be a natural ancient variety of equid. The specific diversity of ABS and the numbers of unique haplotypes in the evolutionary process affords this species a better genetic fitness and ability to adapt to the native environment.

Effects of global warming on sex ratios in fishes

In fishes, sex is determined by genetics, the environment or an interaction of both. Temperature is among the most important environmental factors that can affect sex determination. As a consequence, changes in temperature at critical developmental stages can induce biases in primary sex ratios in some species. However, early sex ratios can also be biased by sex-specific tolerances to environmental stresses that may, in some cases, be amplified by changes in water temperature. Sex-specific reactions to environmental stress have been observed at early larval stages before gonad formation starts. It is therefore necessary to distinguish between temperature effects on sex determination, generally acting through the stress axis or epigenetic mechanisms, and temperature effects on sex-specific mortality. Both are likely to affect sex ratios and hence population dynamics. Moreover, in cases where temperature effects on sex determination lead to genotype-phenotype mismatches, long-term effects on population dynamics are possible, for example temperature-induced masculinization potentially leading to the loss of Y chromosomes or feminization to male-biased operational sex ratios in future generations. To date, most studies under controlled conditions conclude that if temperature affects sex ratios, elevated temperatures mostly lead to a male bias. The few studies that have been performed on wild populations seem to confirm this general trend. Recent findings suggest that transgenerational plasticity could mitigate the effects of warming on sex ratios in some populations.

Exome-wide association study reveals largely distinct gene sets underlying specific resistance to dengue virus types 1 and 3 in Aedes aegypti

Although specific interactions between host and pathogen genotypes have been well documented in invertebrates, the identification of host genes involved in discriminating pathogen genotypes remains a challenge. In the mosquito Aedes aegypti, the main dengue virus (DENV) vector worldwide, statistical associations between host genetic markers and DENV types or strains were previously detected, but the host genes underlying this genetic specificity have not been identified. In particular, it is unknown whether DENV type- or strain-specific resistance relies on allelic variants of the same genes or on distinct gene sets. Here, we investigated the genetic architecture of DENV resistance in a population of Ae. aegypti from Bakoumba, Gabon, which displays a stronger resistance phenotype to DENV type 1 (DENV-1) than to DENV type 3 (DENV-3) infection. Following experimental exposure to either DENV-1 or DENV-3, we sequenced the exomes of large phenotypic pools of mosquitoes that are either resistant or susceptible to each DENV type. Using variation in single-nucleotide polymorphism (SNP) frequencies among the pools, we computed empirical p values based on average gene scores adjusted for the differences in SNP counts, to identify genes associated with infection in a DENV type-specific manner. Among the top 5% most significant genes, 263 genes were significantly associated with resistance to both DENV-1 and DENV-3, 287 genes were only associated with DENV-1 resistance and 290 were only associated with DENV-3 resistance. The shared significant genes were enriched in genes with ATP binding activity and sulfur compound transmembrane transporter activity, whereas the genes uniquely associated with DENV-3 resistance were enriched in genes with zinc ion binding activity. Together, these results indicate that specific resistance to different DENV types relies on largely non-overlapping sets of genes in this Ae. aegypti population and pave the way for further mechanistic studies.

Compatibility between hosts and pathogens is often genetically specific in invertebrates but host genes underlying this genetic specificity have not been elucidated. We investigated the genetic architecture of dengue virus type-specific resistance in the mosquito vector Aedes aegypti. We used a natural population of Ae. aegypti from Bakoumba, Gabon, which is differentially resistant to dengue virus type 1 and dengue virus type 3. We surveyed genetic variation in protein-coding regions of the mosquito genome and compared the frequency of genetic polymorphisms between groups of mosquitoes that are either resistant or susceptible to each dengue virus type. We found that the Ae. aegypti genes associated with resistance to dengue virus type 1 or dengue virus type 3 were largely non-overlapping. This finding indicates that different sets of host genes, rather than different variants of the same genes, confer pathogen-specific resistance in this population. This study is an important step towards identification of mechanisms underlying the genetic specificity of invertebrate host-pathogen interactions.

Major histocompatibility complex class I diversity limits the repertoire of T cell receptors

Major histocompatibility complex (MHC) genes encode proteins that initiate adaptive immune responses through the presentation of foreign antigens to T cells. The high polymorphism found at these genes, thought to be promoted and maintained by pathogen-mediated selection, contrasts with the limited number of MHC loci found in most vertebrates. Although expressing many diverse MHC genes should broaden the range of detectable pathogens, it has been hypothesized to also cause deletion of larger fractions of self-reactive T cells, leading to a detrimental reduction of the T cell receptor (TCR) repertoire. However, a key prediction of this TCR depletion hypothesis, that the TCR repertoire should be inversely related to the individual MHC diversity, has never been tested. Here, using high-throughput sequencing and advanced sequencing error correction, we provide evidence of such an association in a rodent species with high interindividual variation in the number of expressed MHC molecules, the bank vole (Myodes glareolus). Higher individual diversity of MHC class I, but not class II, was associated with smaller TCR repertoires. Our results thus provide partial support for the TCR depletion model, while also highlighting the complex, potentially MHC class-specific mechanisms by which autoreactivity may trade off against evolutionary expansion of the MHC gene family.

Demographic, environmental and genetic determinants of mating success in captive koalas (Phascolarctos cinereus)

Many factors have been shown to affect mating behavior. For instance, genes of the major histocompatibility complex (MHC) are known to influence mate choice in a wide variety of vertebrate species. The genetic management of captive populations can be confounded if intrinsic mate choice reduces or eliminates reproductive success between carefully chosen breeding pairs. For example, the San Diego Zoo koala colony only has a 45% copulation rate for matched individuals. Herein, we investigated determinants of koala mating success using breeding records (1984-2010) and genotypes for 52 individuals at four MHC markers. We quantified MHC diversity according to functional amino acids, heterozygosity, and the probability of producing a heterozygous offspring. We then used categorical analysis and logistic regression to investigate both copulation and parturition success. In addition, we also examined age, day length, and average pairwise kinship. Our post-hoc power analysis indicates that at a power level of 1-β = 0.8, we should have been able to detect strong MHC preferences. However, we did not find a significant MHC effect on either copulation or parturition success with one exception: pairs with lower or no production of a joey had significantly lower MHC functional amino acid diversity in the categorical analysis. In contrast, day length and dam age (or age difference of the pair) consistently had an effect on mating success. These findings may be leveraged to improve the success of attempted pairs, conserve resources, and facilitate genetic management.

Specific MHC class I supertype associated with parasite infection and color morph in a wild lizard population

The major histocompatibility complex (MHC) is a large gene family that plays a central role in the immune system of all jawed vertebrates. Nonavian reptiles are underrepresented within the MHC literature and little is understood regarding the mechanisms maintaining MHC diversity in this vertebrate group. Here, we examined the relative roles of parasite-mediated selection and sexual selection in maintaining MHC class I diversity of a color polymorphic lizard. We discovered evidence for parasite-mediated selection acting via rare-allele advantage or fluctuating selection as ectoparasite load was significantly lower in the presence of a specific MHC supertype (functional clustering of alleles): supertype four. Based on comparisons between ectoparasite prevalence and load, and assessment of the impact of ectoparasite load on host fitness, we suggest that supertype four confers quantitative resistance to ticks or an intracellular tickborne parasite. We found no evidence for MHC-associated mating in terms of pair genetic distance, number of alleles, or specific supertypes. An association was uncovered between supertype four and male throat color morph. However, it is unlikely that male throat coloration acts as a signal of MHC genotype to conspecifics because we found no evidence to suggest that male throat coloration predicts male mating status. Overall, our results suggest that parasite-mediated selection plays a role in maintaining MHC diversity in this population via rare-allele advantage and/or fluctuating selection. Further work is required to determine whether sexual selection also plays a role in maintaining MHC diversity in agamid lizards.

Analysis of peptide-SLA binding by establishing immortalized porcine alveolar macrophage cells with different SLA class II haplotypes

Primary porcine alveolar macrophages (PAM) are useful for studying viral infections and immune response in pigs; however, long-term use of these cells is limited by the cells’ short lifespan. We immortalized primary PAMs by transfecting them with both hTERT and SV40LT and established two immortalized cell lines (iPAMs) actively proliferating even after 35 passages. These cells possessed the characteristics of primary PAMs, including strong expression of swine leukocyte antigen (SLA) class II genes and the inability to grow anchorage-independently. We characterized their SLA genes and subsequently performed peptide-SLA binding assays using a peptide from porcine circovirus type 2 open reading frame 2 to experimentally measure the binding affinity of the peptide to SLA class II. The number of peptides bound to cells measured by fluorescence was very low for PK15 cells (7.0% ± 1.5), which are not antigen-presenting cells, unlike iPAM61 (33.7% ± 3.4; SLA-DQA*0201/0303, DQB1*0201/0901, DRB1*0201/1301) and iPAM303 (73.3% ± 5.4; SLA DQA*0106/0201, DQB1*0202/0701, DRB1*0402/0602). The difference in peptide binding between the two iPAMs was likely due to the allelic differences between the SLA class II molecules that were expressed. The development of an immortal PAM cell panel harboring diverse SLA haplotypes and the use of an established method in this study can become a valuable tool for evaluating the interaction between antigenic peptides and SLA molecules and is important for many applications in veterinary medicine including vaccine development.

Rapid Antagonistic Coevolution in an Emerging Pathogen and Its Vertebrate Host

Host-pathogen coevolution is assumed to play a key role in eco-evolutionary processes, including epidemiological dynamics and the evolution of sexual reproduction [1-4]. Despite this, direct evidence for host-pathogen coevolution is exceptional [5-7], particularly in vertebrate hosts. Indeed, although vertebrate hosts have been shown to evolve in response to pathogens or vice versa [8-12], there is little evidence for the necessary reciprocal changes in the success of both antagonists over time [13]. Here, we generate a time-shift experiment to demonstrate adaptive, reciprocal changes in North American house finches (Haemorhous mexicanus) and their emerging bacterial pathogen, Mycoplasma gallisepticum [14-16]. Our experimental design is made possible by the existence of disease-exposed and unexposed finch populations, which were known to exhibit equivalent responses to experimental inoculation until the recent spread of genetic resistance in the former [14, 17]. Whereas inoculations with pathogen isolates from epidemic outbreak caused comparable sub-lethal eye swelling in hosts from exposed (hereafter adapted) and unexposed (hereafter ancestral) populations, inoculations with isolates sampled after the spread of resistance were threefold more likely to cause lethal symptoms in hosts from ancestral populations. Similarly, the probability that pathogens successfully established an infection in the primary host and, before inducing death, transmitted to an uninfected sentinel was highest when recent isolates were inoculated in hosts from ancestral populations and lowest when early isolates were inoculated in hosts from adapted populations. Our results demonstrate antagonistic host-pathogen coevolution, with hosts and pathogens displaying increased resistance and virulence in response to each other over time.

Immunogenetic novelty confers a selective advantage in host-pathogen coevolution

The major histocompatibility complex (MHC) is crucial to the adaptive immune response of vertebrates and is among the most polymorphic gene families known. Its high diversity is usually attributed to selection imposed by fast-evolving pathogens. Pathogens are thought to evolve to escape recognition by common immune alleles, and, hence, novel MHC alleles, introduced through mutation, recombination, or gene flow, are predicted to give hosts superior resistance. Although this theoretical prediction underpins host-pathogen "Red Queen" coevolution, it has not been demonstrated in the context of natural MHC diversity. Here, we experimentally tested whether novel MHC variants (both alleles and functional "supertypes") increased resistance of guppies (Poecilia reticulata) to a common ectoparasite (Gyrodactylus turnbulli). We used exposure-controlled infection trials with wild-sourced parasites, and Gyrodactylus-naïve host fish that were F2 descendants of crossed wild populations. Hosts carrying MHC variants (alleles or supertypes) that were new to a given parasite population experienced a 35-37% reduction in infection intensity, but the number of MHC variants carried by an individual, analogous to heterozygosity in single-locus systems, was not a significant predictor. Our results provide direct evidence of novel MHC variant advantage, confirming a fundamental mechanism underpinning the exceptional polymorphism of this gene family and highlighting the role of immunogenetic novelty in host-pathogen coevolution.

Experimental manipulation of population-level MHC diversity controls pathogen virulence evolution in Mus musculus

The virulence levels attained by serial passage of pathogens through similar host genotypes are much higher than observed in natural systems; however, it is unknown what keeps natural virulence levels below these empirically demonstrated maximum levels. One hypothesis suggests that host diversity impedes pathogen virulence, because adaptation to one host genotype carries trade-offs in the ability to replicate and cause disease in other host genotypes. To test this hypothesis, with the simplest level of population diversity within the loci of the major histocompatibility complex (MHC), we serially passaged Friend virus complex (FVC) through two rounds, in hosts with either the same MHC genotypes (pure passage) or hosts with different MHC genotypes (alternated passage). Alternated passages showed a significant overall reduction in viral titre (31%) and virulence (54%) when compared to pure passages. Furthermore, a resistant host genotype initially dominated any effects due to MHC diversity; however, when FVC was allowed to adapt to the resistant host genotype, predicted MHC effects emerged; that is, alternated lines show reduced virulence. These data indicate serial exposure to diverse MHC genotypes is an impediment to pathogen adaptation, suggesting genetic variation at MHC loci is important for limiting virulence in a rapidly evolving pathogen and supports negative frequency-dependent selection as a force maintaining MHC diversity in host populations.

Interpreting T-Cell Cross-reactivity through Structure: Implications for TCR-Based Cancer Immunotherapy

Immunotherapy has become one of the most promising avenues for cancer treatment, making use of the patient’s own immune system to eliminate cancer cells. Clinical trials with T-cell-based immunotherapies have shown dramatic tumor regressions, being effective in multiple cancer types and for many different patients. Unfortunately, this progress was tempered by reports of serious (even fatal) side effects. Such therapies rely on the use of cytotoxic T-cell lymphocytes, an essential part of the adaptive immune system. Cytotoxic T-cells are regularly involved in surveillance and are capable of both eliminating diseased cells and generating protective immunological memory. The specificity of a given T-cell is determined through the structural interaction between the T-cell receptor (TCR) and a peptide-loaded major histocompatibility complex (MHC); i.e., an intracellular peptide–ligand displayed at the cell surface by an MHC molecule. However, a given TCR can recognize different peptide–MHC (pMHC) complexes, which can sometimes trigger an unwanted response that is referred to as T-cell cross-reactivity. This has become a major safety issue in TCR-based immunotherapies, following reports of melanoma-specific T-cells causing cytotoxic damage to healthy tissues (e.g., heart and nervous system). T-cell cross-reactivity has been extensively studied in the context of viral immunology and tissue transplantation. Growing evidence suggests that it is largely driven by structural similarities of seemingly unrelated pMHC complexes. Here, we review recent reports about the existence of pMHC “hot-spots” for cross-reactivity and propose the existence of a TCR interaction profile (i.e., a refinement of a more general TCR footprint in which some amino acid residues are more important than others in triggering T-cell cross-reactivity). We also make use of available structural data and pMHC models to interpret previously reported cross-reactivity patterns among virus-derived peptides. Our study provides further evidence that structural analyses of pMHC complexes can be used to assess the intrinsic likelihood of cross-reactivity among peptide-targets. Furthermore, we hypothesize that some apparent inconsistencies in reported cross-reactivities, such as a preferential directionality, might also be driven by particular structural features of the targeted pMHC complex. Finally, we explain why TCR-based immunotherapy provides a special context in which meaningful T-cell cross-reactivity predictions can be made.

Evolutionary and functional implications of hypervariable loci within the skin virome

Localized genomic variability is crucial for the ongoing conflicts between infectious microbes and their hosts. An understanding of evolutionary and adaptive patterns associated with genomic variability will help guide development of vaccines and antimicrobial agents. While most analyses of the human microbiome have focused on taxonomic classification and gene annotation, we investigated genomic variation of skin-associated viral communities. We evaluated patterns of viral genomic variation across 16 healthy human volunteers. Human papillomavirus (HPV) and Staphylococcus phages contained 106 and 465 regions of diversification, or hypervariable loci, respectively. Propionibacterium phage genomes were minimally divergent and contained no hypervariable loci. Genes containing hypervariable loci were involved in functions including host tropism and immune evasion. HPV and Staphylococcus phage hypervariable loci were associated with purifying selection. Amino acid substitution patterns were virus dependent, as were predictions of their phenotypic effects. We identified diversity generating retroelements as one likely mechanism driving hypervariability. We validated these findings in an independently collected skin metagenomic sequence dataset, suggesting that these features of skin virome genomic variability are widespread. Our results highlight the genomic variation landscape of the skin virome and provide a foundation for better understanding community viral evolution and the functional implications of genomic diversification of skin viruses.

Functionally Diverse NK-Like T Cells Are Effectors and Predictors of Successful Aging

The fundamental challenge of aging and long-term survivorship is maintenance of functional independence and compression of morbidity despite a life history of disease. Inasmuch as immunity is a determinant of individual health and fitness, unraveling novel mechanisms of immune homeostasis in late life is of paramount interest. Comparative studies of young and old persons have documented age-related atrophy of the thymus, the contraction of diversity of the T cell receptor (TCR) repertoire, and the intrinsic inefficiency of classical TCR signaling in aged T cells. However, the elderly have highly heterogeneous health phenotypes. Studies of defined populations of persons aged 75 and older have led to the recognition of successful aging, a distinct physiologic construct characterized by high physical and cognitive functioning without measurable disability. Significantly, successful agers have a unique T cell repertoire; namely, the dominance of highly oligoclonal αβT cells expressing a diverse array of receptors normally expressed by NK cells. Despite their properties of cell senescence, these unusual NK-like T cells are functionally active effectors that do not require engagement of their clonotypic TCR. Thus, NK-like T cells represent a beneficial remodeling of the immune repertoire with advancing age, consistent with the concept of immune plasticity. Significantly, certain subsets are predictors of physical/cognitive performance among older adults. Further understanding of the roles of these NK-like T cells to host defense, and how they integrate with other physiologic domains of function are new frontiers for investigation in Aging Biology. Such pursuits will require a research paradigm shift from the usual young-versus-old comparison to the analysis of defined elderly populations. These endeavors may also pave way to age-appropriate, group-targeted immune interventions for the growing elderly population.

MHC-correlated preferences in diestrous female horses (Equus caballus)

Genes of the major histocompatibility complex (MHC) have been shown to influence communication in many vertebrates, possibly with context-specific MHC-correlated reactions. Here we test for MHC-linked female preferences in the polygynous horse (Equus caballus) by repeatedly exposing 19 mares to a group of seven sexually experienced stallions. Each mare was tested four times during two consecutive reproductive cycles, twice during estrus and twice during diestrus. Male plasma testosterone concentrations were determined from weekly blood samples, and equine leukocyte antigen (ELA) class I and II alleles were determined serologically at the end of the experiments. Perception of male attractiveness was strongly dependent on estrous cycle: mean preference scores did not correlate for mares in diestrus and estrus and varied more during estrus than during diestrus. We found elevated female interests for MHC-dissimilar stallions, but only during diestrus, not during estrus. Female preferences were not significantly predicted by mean male testosterone plasma concentrations. However, testosterone concentrations changed during the 11 weeks of the experiment. By the end of the experiment, average testosterone concentration was significantly correlated to the average number of MHC alleles the stallions shared with the mares. We conclude that the MHC affects female preferences for stallions, but non-MHC linked male characteristics can overshadow effects of the MHC during estrus.

Extensive Allelic Diversity of MHC Class I in Wild Mallard Ducks

MHC class I is critically involved in defense against viruses, and diversity from polygeny and polymorphism contributes to the breadth of the immune response and health of the population. In this article, we examine MHC class I diversity in wild mallard ducks, the natural host and reservoir of influenza A viruses. We previously showed domestic ducks predominantly use UAA, one of five MHC class I genes, but whether biased expression is also true for wild mallards is unknown. Using RT-PCR from blood, we examined expressed MHC class I alleles from 38 wild mallards (Anas platyrhynchos) and identified 61 unique alleles, typically 1 or 2 expressed alleles in each individual. To determine whether expressed alleles correspond to UAA adjacent to TAP2 as in domestic ducks, we cloned and sequenced genomic UAA-TAP2 fragments from all mallards, which matched transcripts recovered and allowed us to assign most alleles as UAA Allelic differences are primarily located in α1 and α2 domains in the residues known to interact with peptide in mammalian MHC class I, suggesting the diversity is functional. Most UAA alleles have unique residues in the cleft predicting distinct specificity; however, six alleles have an unusual conserved cleft with two cysteine residues. Residues that influence peptide-loading properties and tapasin involvement in chicken are fixed in duck alleles and suggest tapasin independence. Biased expression of one MHC class I gene may make viral escape within an individual easy, but high diversity in the population places continual pressure on the virus in the reservoir species.

Structures of Preferred Human IgV Genes-Based Protective Antibodies Identify How Conserved Residues Contact Diverse Antigens and Assign Source of Specificity to CDR3 Loop Variation

The human Ab response to certain pathogens is oligoclonal, with preferred IgV genes being used more frequently than others. A pair of such preferred genes, IGVK3-11 and IGVH3-30, contributes to the generation of protective Abs directed against the 23F serotype of the pneumonococcal capsular polysaccharide of Streptococcus pneumoniae and against the AD-2S1 peptide of the gB membrane protein of human CMV. Structural analyses of Fab fragments of mAbs 023.102 and pn132p2C05 in complex with portions of the 23F polysaccharide revealed five germline-encoded residues in contact with the key component, l-rhamnose. In the case of the AD-2S1 peptide, the KE5 Fab fragment complex identified nine germline-encoded contact residues. Two of these germline-encoded residues, Arg91L and Trp94L, contact both the l-rhamnose and the AD-2S1 peptide. Comparison of the respective paratopes that bind to carbohydrate and protein reveals that stochastic diversity in both CDR3 loops alone almost exclusively accounts for their divergent specificity. Combined evolutionary pressure by human CMV and the 23F serotype of S. pneumoniae acted on the IGVK3-11 and IGVH3-30 genes as demonstrated by the multiple germline-encoded amino acids that contact both l-rhamnose and AD-2S1 peptide.

Experimental evolution in silico: a custom-designed mathematical model for virulence evolution of Bacillus thuringiensis

Most mathematical models on the evolution of virulence are based on epidemiological models that assume parasite transmission follows the mass action principle. In experimental evolution, however, mass action is often violated due to controlled infection protocols. This "theory-experiment mismatch" raises the question whether there is a need for new mathematical models to accommodate the particular characteristics of experimental evolution. Here, we explore the experimental evolution model system of Bacillus thuringiensis as a parasite and Caenorhabditis elegans as a host. Recent experimental studies with strict control of parasite transmission revealed that one-sided adaptation of B. thuringiensis with non-evolving hosts selects for intermediate or no virulence, sometimes coupled with parasite extinction. In contrast, host-parasite coevolution selects for high virulence and for hosts with strong resistance against B. thuringiensis. In order to explain the empirical results, we propose a new mathematical model that mimics the basic experimental set-up. The key assumptions are: (i) controlled parasite transmission (no mass action), (ii) discrete host generations, and (iii) context-dependent cost of toxin production. Our model analysis revealed the same basic trends as found in the experiments. Especially, we could show that resistant hosts select for highly virulent bacterial strains. Moreover, we found (i) that the evolved level of virulence is independent of the initial level of virulence, and (ii) that the average amount of bacteria ingested significantly affects the evolution of virulence with fewer bacteria ingested selecting for highly virulent strains. These predictions can be tested in future experiments. This study highlights the usefulness of custom-designed mathematical models in the analysis and interpretation of empirical results from experimental evolution.

INTEGRATING PARASITES AND PATHOGENS INTO THE STUDY OF GEOGRAPHIC RANGE LIMITS

The geographic distributions of all species are limited, and the determining factors that set these limits are of fundamental importance to the fields of ecology and evolutionary biology. Plant and animal ranges have been of primary concern, while those of parasites, which represent much of the Earth's biodiversity, have been neglected. Here, we review the determinants of the geographic ranges of parasites and pathogens, and explore how parasites provide novel systems with which to investigate the ecological and evolutionary processes governing host/parasite spatial distributions. Although there is significant overlap in the causative factors that determine range borders of parasites and free-living species, parasite distributions are additionally constrained by the geographic range and ecology of the host species' population, as well as by evolutionary factors that promote host-parasite coevolution. Recently, parasites have been used to infer population demographic and ecological information about their host organisms and we conclude that this strategy can be further exploited to understand geographic range limitations of both host and parasite populations.

Genetics of Autoimmune Thyroiditis in Type 1 Diabetes Reveals a Novel Association With DPB1*0201: Data From the Type 1 Diabetes Genetics Consortium

BACKGROUND

Autoimmune thyroiditis occurs in 10-25% of patients with type 1 diabetes (T1D). Most of these patients are also positive for thyroid peroxidase (TPO) antibodies. Thyroid dysfunction complicates T1D metabolic control and is a component of the autoimmune polyglandular syndrome (APS, type 2 or 3). Previous studies of isolated T1D and of T1D combined with other autoimmune disorders showed genetic susceptibility for alleles in HLA-DQB1 and -DRB1 and also CTLA4 and PTPN22.

RESEARCH DESIGN AND METHODS

We analyzed the Type 1 Diabetes Genetics Consortium Autoantibody Workshop data by differentiating those T1D probands with and without TPO antibodies or thyroid disease with respect to polymorphisms in HLA, CTLA4, INS, PTPN22, and VDR, taking into account the ethnic origin. Genotype and clinical/immunogenic phenotype data were analyzed by gene counting methods and logistic regression analysis.

RESULTS

The presence of TPO antibodies (25.2%) and thyroid disease (8.4%) was associated with older age, female sex, and presence of other autoantibodies (GAD65, ATPase, 21-OH) (all P<0.001). The highest prevalence was in patients of Hispanic ancestry (31%) and the lowest in those of African ancestry (8%). In T1D non-Hispanic whites, HLA-DRB1*0101 is significantly (P<0.0001) less frequent in TPO-positive than in TPO-negative individuals, whereas HLA-DRB1*0404, -DQB1*0301, and -DPB1*0201 are significantly (P<0.0001) more frequent. Subjects with a high titer of TPO autoantibodies and with thyroid disease were associated with female sex and older age and negatively associated with DRB1*0401-DQB1*0302 (P<0.0001). No significant differences were observed for an association of TPO positivity or thyroid disease with single nucleotide polymorphisms in the INS, CTLA4, or VDR loci, with nominal significance (P=0.01) for PTPN22 R620W variant.

CONCLUSIONS

Thyroid autoimmunity is highly prevalent in T1D patients of non-Hispanic white, Asian, or Hispanic origin. The strongest disease risk is conferred by female sex and older age. This risk is modulated by HLA-DRB1 and HLA-DPB1 loci. The immunogenetic profile for T1D with thyroid autoimmunity may identify distinct pathways regulating polyglandular autoimmunity and disease.

Serial infection of diverse host (Mus) genotypes rapidly impedes pathogen fitness and virulence

Reduced genetic variation among hosts may favour the emergence of virulent infectious diseases by enhancing pathogen replication and its associated virulence due to adaptation to a limited set of host genotypes. Here, we test this hypothesis using experimental evolution of a mouse-specific retroviral pathogen, Friend virus (FV) complex. We demonstrate rapid fitness (i.e. viral titre) and virulence increases when FV complex serially infects a series of inbred mice representing the same genotype, but not when infecting a diverse array of inbred mouse strains modelling the diversity in natural host populations. Additionally, a single infection of a different host genotype was sufficient to constrain the emergence of a high fitness/high virulence FV complex phenotype in these experiments. The potent inhibition of viral fitness and virulence was associated with an observed loss of the defective retroviral genome (spleen focus-forming virus), whose presence exacerbates infection and drives disease in susceptible mice. Results from our experiments provide an important first step in understanding how genetic variation among vertebrate hosts influences pathogen evolution and suggests that serial exposure to different genotypes within a single host species may act as a constraint on pathogen adaptation that prohibits the emergence of more virulent infections. From a practical perspective, these results have implications for low-diversity host populations such as endangered species and domestic animals.

No more non-model species: the promise of next generation sequencing for comparative immunology

Next generation sequencing (NGS) allows for the rapid, comprehensive and cost effective analysis of entire genomes and transcriptomes. NGS provides approaches for immune response gene discovery, profiling gene expression over the course of parasitosis, studying mechanisms of diversification of immune receptors and investigating the role of epigenetic mechanisms in regulating immune gene expression and/or diversification. NGS will allow meaningful comparisons to be made between organisms from different taxa in an effort to understand the selection of diverse strategies for host defence under different environmental pathogen pressures. At the same time, it will reveal the shared and unique components of the immunological toolkit and basic functional aspects that are essential for immune defence throughout the living world. In this review, we argue that NGS will revolutionize our understanding of immune responses throughout the animal kingdom because the depth of information it provides will circumvent the need to concentrate on a few "model" species.

Evolutionary constraints in variable environments, from proteins to networks

Environmental changes can not only trigger a regulatory response, but also impose evolutionary pressures that can modify the underlying regulatory network. Here, we review recent approaches that are beginning to disentangle this complex interplay between regulatory and evolutionary responses. Systematic genetic reconstructions have shown how evolutionary constraints arise from epistatic interactions between mutations in fixed environments. This approach is now being extended to more complex environments and systems. The first results suggest that epistasis is affected dramatically by environmental changes and, hence, can profoundly affect the course of evolution. Thus, external environments not only define the selection of favored phenotypes, but also affect the internal constraints that can limit the evolution of these phenotypes. These findings also raise new questions relating to the conditions for evolutionary transitions and the evolutionary potential of regulatory networks.

Host adaptation to viruses relies on few genes with different cross-resistance properties

Host adaptation to one parasite may affect its response to others. However, the genetics of these direct and correlated responses remains poorly studied. The overlap between these responses is instrumental for the understanding of host evolution in multiparasite environments. We determined the genetic and phenotypic changes underlying adaptation of Drosophila melanogaster to Drosophila C virus (DCV). Within 20 generations, flies selected with DCV showed increased survival after DCV infection, but also after cricket paralysis virus (CrPV) and flock house virus (FHV) infection. Whole-genome sequencing identified two regions of significant differentiation among treatments, from which candidate genes were functionally tested with RNAi. Three genes were validated--pastrel, a known DCV-response gene, and two other loci, Ubc-E2H and CG8492. Knockdown of Ubc-E2H and pastrel also led to increased sensitivity to CrPV, whereas knockdown of CG8492 increased susceptibility to FHV infection. Therefore, Drosophila adaptation to DCV relies on few major genes, each with different cross-resistance properties, conferring host resistance to several parasites.