MHC-I alleles mediate clearance and antibody response to the zoonotic Lassa virus in Mastomys rodent reservoirs

Abstract

West African Mastomys rodents are the primary reservoir of the zoonotic Lassa virus (LASV). The virus causes haemorrhagic Lassa fever and considerable mortality in humans. To date, the role of Mastomys immunogenetics in resistance to, and persistence of, LASV infections is largely unknown. Here, we investigated the role of Major Histocompatibility Complex class I (MHC-I) on LASV infection status (i.e., active vs. cleared infection, determined via PCR and an immunofluorescence assay on IgG antibodies, respectively) in Mastomys natalensis and M. erythroleucus sampled within southwestern Nigeria. We identified more than 190 and 90 MHC-I alleles by Illumina high throughput-sequencing in M. natalensis and M. erythroleucus, respectively, with different MHC allele compositions and frequencies between LASV endemic and non-endemic sites. In M. natalensis, the MHC allele ManaMHC-I*006 was negatively associated with active infections (PCR-positive) and positively associated with cleared infections (IgG-positive) simultaneously, suggesting efficient immune responses that facilitate LASV clearance in animals carrying this allele. Contrarily, alleles ManaMHC-I*008 and ManaMHC-I*021 in M. natalensis, and MaerMHC-I*008 in M. erythroleucus, were positively associated with active infection, implying susceptibility. Alleles associated with susceptibility shared a glutamic acid at the positively selected codon 57, while ManaMHC-I*006 featured an arginine. There was no link between number of MHC alleles per Mastomys individual and LASV prevalence. Thus, specific alleles, but not MHC diversity per se, seem to mediate antibody responses to viremia. We conclude that co-evolution with LASV likely shaped the MHC-I diversity of the main LASV reservoirs in southwestern Nigeria, and that information on reservoir immunogenetics may hold insights into transmission dynamics and zoonotic spillover risks.

Citation:Olayemi A, Schmid DW, Fleischer R, Wilhelm K, Heni AC, Mueller-Klein N, et al. (2024) MHC-I alleles mediate clearance and antibody response to the zoonotic Lassa virus in Mastomys rodent reservoirs. PLoS Negl Trop Dis 18(2): e0011984. https://doi.org/10.1371/journal.pntd.0011984

Editor:Nam-Hyuk Cho, Seoul National University College of Medicine, REPUBLIC OF KOREA

Received:April 18, 2023; Accepted:February 9, 2024; Published: February 29, 2024

Copyright: © 2024 Olayemi et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability:Raw data and scripts are freely available on https://figshare.com/projects/DATA_MHC-I_alleles_mediate_clearance_and_antibody_response_to_the_zoonotic_Lassa_virus_in_Mastomys_rodent_reservoirs/187320.

Funding:This study was funded by a Georg-Forster Alexander von Humboldt (AvH) fellowship to AO and German Research Foundation grants (DFG, GU 883/4-1 & GU 883/4-2) to SG & AO. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

1. Introduction

Lassa fever, caused by the Lassa virus (LASV), is a lethal zoonotic haemorrhagic disease endemic to West Africa, where it is estimated to cause up to 18,000 human deaths annually [1]. Infection during pregnancy has been reported to result in up to 50% maternal- and 87% foetal mortality [2–4]. About 30% of survivors from this illness suffer sensorineural hearing loss [5]. With a rising number of cases (Nigeria accounting for up to 90–99% of all clinical diagnoses [6]), the possibility of Lassa fever spreading outside the West African region is cause for concern [7,8]. Unsurprisingly, in the recently revised WHO R&D blueprint on emerging infectious diseases, LASV control was marked as high priority [9–11].

The main natural reservoir for LASV are rodents of the family Muridae [12], which do not display significant adverse effects from infection [13]. Humans can contract LASV when they come into contact with excreta of viruric rodents [14]. LASV infections in humans are generally congruent with distribution of the virus across rodent populations. However, this virus distribution tends to be concentrated in particular hotspots [15]. Projecting which rodent individuals and populations are prone to LASV occurrence will help predict the likelihood of emergence and, accordingly, transmission to humans.

The Natal multimammate mouse Mastomys natalensis, the most widespread and abundant rodent across Africa [16], is recognized as the primary LASV reservoir [12]. It is the natural host of a range of additional mammarenaviruses in and beyond West Africa (e.g., Dhati Welel [17], Gairo [18], Luna [19], Mobala-like [20], Mopeia [21] and Morogoro [22] viruses), and also host to other zoonotic pathogens (e.g., Yersinia pestis [23], Leptospira sp. [24], and Trypanosoma sp. [25]). Aside from M. natalensis, the Guinea multimammate mouse M. erythroleucus and other murid rodent species have been demonstrated as independent LASV reservoirs [26]. Though similar in external morphology and ecology, M. natalensis and M. erythroleucus are distinct species possessing different chromosome numbers (2n = 32 and 2n = 38, respectively [27]). Geographically, M. erythroleucus is limited to West- and Central Africa, where it lives in sympatry with the more widely distributed M. natalensis [27]. Yet, there are locations where only one species is present within Nigeria. Even though these multimammate mice function as reservoirs of human LASV, details concerning ecology of the virus in rodents and its connection to host immunogenetics are still unclear.

The Major Histocompatibility Complex (MHC) is a group of genes that mediate the recognition and presentation of self and non-self peptides on cell surfaces in order to trigger adaptive immune responses [28,29]. Glycoproteins coded by MHC-I primarily bind to peptides from intracellular microbes (e.g., viruses), whereas those coded by MHC-II predominantly bind to peptides from extracellular organisms (e.g., bacteria and helminths) [30–32]. MHC-II composition and diversity have been associated with resistance to parasites in a variety of non-model animal hosts [33–35] including wild Muridae (e.g., [36–38]). The allelic diversity at the MHC-II region of Mastomys natalensis was described previously [39]; while, for MHC-I, only orthologous sequences of the Mastomys coucha genome were identified [40]. Generally, descriptions of MHC-I composition and diversity in wild Muridae are scant. Mastomys MHC-I has never been investigated or linked to mammarenaviruses like LASV, even though MHC-I-mediated clearance of particular LASV peptides was already found more than a decade ago [41]. This is surprising given the zoonotic potential of this rodent family that copes well with habitat changes and often lives in close proximity to humans [42,43].

To fill this knowledge gap, we aimed to decipher the allelic MHC-I diversity in two Mastomys species (M. natalensis and M. erythroleucus) captured from two sites within southwestern Nigeria—one where LASV is endemic and the other with no history of infected rodents–in order to identify links between MHC-I genotype and LASV infection. Apart from rodents infected as neonates (which develop a perpetual infection, perhaps for life), the majority of LASV (and other closely-related mammarenavirus) infections in Mastomys are acute and transient, with individuals retaining anti-LASV IgG antibodies after virus clearance [44–47]. Consequently, individuals can be classified as carrying an ongoing LASV infection if PCR-positive or having cleared a past infection if IgG-positive. Regarding active infections, a recent capture-recapture study demonstrated that, even among individuals generally considered to be acutely infected, some carry the virus longer than others [48]. Lingering infections promote rodent-to-rodent transmission, and eventual spillover to humans. But, whether individual MHC-I composition is linked to the capacity to efficiently clear LASV remains an open question.

To investigate the link between immunogenetics and LASV infection, we used a co-occurrence analysis approach, combining Mastomys LASV infection status assessment with high throughput sequencing of MHC-I exon 2 to evaluate whether certain alleles were significantly a) negatively associated with active LASV infections (indicating resistance), b) positively associated with active infections (indicating susceptibility), and c) positively associated with cleared infections (indicating capability to mount an effective antibody response). Since we do not have any information on infection intensity or sero-levels of IgG [49], we define resistance and susceptibility as a rodent’s ability to effectively and timely clear an active LASV infection. We discuss our results in light of site-specific differences in LASV endemism since reservoir immunogenetics may provide insights into transmission dynamics and zoonotic spillover risks.

2. Methods

2.1 Ethics statement

This study analysed tissues obtained from Mastomys rodents trapped as part of our larger project, “Seroprevalence, Incidence and Risk Factors of Lassa Fever, Edo State, Nigeria”, approved by the Health Research Ethics Committee of the Irrua Specialist Teaching Hospital in Edo State, Nigeria (ISTH/HREC/2017/1019/28). No animal experiments were conducted. Rodents were euthanised using the anaesthetic isoflurane before biopsy collection, with strict safety precautions observed during their capture and necropsy [50].

2.2 Study area and sample collection

The geographic distribution of Lassa fever in Africa, the location of the study sites within Nigeria and the schematic representation of the course of an acute infection with LASV and other closely-related mammarenaviruses hosted by natural populations of the rodent Mastomys natalensis is detailed in Fig 1. Mastomys were captured during previous ecological investigations to screen for LASV in southwestern Nigeria [20,51], and comprised 204 specimens of Mastomys natalensis and 69 of M. erythroleucus. Both Mastomys species were captured in Abagboro (7° 32´N, 4° 30´E) during January and September of 2011 and 2012 and in Ebudin (6° 35´N, 6° 10´E, approximately 292 km apart) during July, October, January and April from 2014 to 2016. Ebudin is located in the endemic zone for Lassa fever in Nigeria, while Abagboro is considered to be in the non-epidemic part of the country. Rodent sampling in Ebudin, where LASV is frequently detected in Mastomys rodents [51], allowed us explore the association between the occurrence of MHC-I alleles and LASV prevalence. Additionally, trapping in Abagboro enabled us to compare MHC-I constitution in rodents between endemic and non-endemic localities.

[Figure omitted. See PDF.]

Fig 1. Geographic distribution of Lassa fever, depiction of the study area and the acute nature of most mammarenavirus infections in Mastomys rodents.

A) Countries endemic for Lassa fever across West Africa, and endemic (Ebudin in yellow) and non-endemic (Abagboro in blue) sampling locations for Mastomys rodents within Nigeria. Sample sizes for M. natalensis and M. erythroleucus analysed for MHC-I in Abagboro and Ebudin are indicated below the Nigerian map. B) Schematic representation of the course of an acute infection with LASV and other closely-related mammarenaviruses hosted by natural populations of the rodent Mastomys natalensis. Infection progresses from active infection detected only via PCR (PCR+, t0 –t1), active infection and initiated seroconversion (PCR+ & IgG+, t1- t2), to IgG antibody presence without active infection (IgG+, > t2), indicative of a cleared infection. The table below the graph summarises the results from various mammarenavirus surveys, and how many rodents were identified to be acutely infected without showing an IgG response, infected and showing IgG responses, and only showing an IgG signal. Note: individuals found within the t1 –t2 category are in the middle of seroconversion; or denote chronic infections (not represented in the graph) in rodents infected from when they were newly born. MOPV: Mopeia virus. MOB-like: Mobala-like virus. The figure was created with BioRender.com.

https://doi.org/10.1371/journal.pntd.0011984.g001

Rodents were captured in houses and surrounding vegetation using Sherman live-traps. Specimens were euthanised with the anaesthetic isoflurane and blood drawn by cardiac puncture. Using whole-blood samples, active LASV infection was determined via PCR [20,46,51] and previous infection by an Immunofluorescence Assay (IFA) to detect IgG antibodies [46,52]. The virus strains from all PCR-positive Mastomys individuals in this study belong to LASV lineage II [51], known to circulate in southern Nigeria. Because Mastomys species can be difficult to distinguish by external morphology, individual rodents were identified by sequencing a segment of mitochondrial Cytochrome b DNA from ethanol-preserved rodent kidneys, extracted using the DNeasy Blood & Tissue kit (Qiagen, Germany) according to manufacturer’s instructions. Genomic DNA extracted from rodent kidneys was also used for MHC-I characterization.

2.3 Primer design and library preparation

To the best of our knowledge, MHC-I has not been sequenced previously for rodents of the genus Mastomys. Aligning the few publicly available MHC-I sequences from the rodent family Muridae (e.g., [40,53]) in Geneious Prime version 2021.2 (Biomatters Inc., USA), we designed primers to amplify MHC-I exon 2 of Mastomys specimens using Primer3 [54]. The primer sequences Ma16F (forward primer, YCCCAGGCWCACACTCRMTGCG) and Ma257R (reverse primer, CTCACCGYCCTCGCTCTGG) are located mainly within the 270 bp of exon 2, with only 6 and 5 bp, respectively, extending into adjacent introns.

In order to generate an Illumina MiSeq DNA Amplicon library, a two-step PCR was carried out following the four-primer amplicon tagging scheme of Fluidigm (Access Array System for Illumina Sequencing Systems, Fluidigm, USA). For the first target specific (TS) PCR the primers were appended with either a CS1 adapter (CS1-4N-Ma16F) or CS2 adapter (CS2-Ma257). Four random bases (4N) were added to the forward primer to facilitate cluster identification on the Illumina flow cell. During the second barcoding PCR a sample-specific barcode (10bp) and Illumina adapter sequences were added to the CS adapters. The first TS-PCR of 10 μl-units included 1 μl template DNA, 5 μl AmpliGold 360 MasterMix (Thermo Fisher Scientific, USA), 1 μl GC buffer, 300nM of each CS primer and 2.4 μl ultrapure water. The temperature profile consisted of an initial denaturation phase at 95°C for 600 s; 30 cycles of denaturation at 95°C for 30 s, annealing at 59°C for 30 s, elongation at 72°C for 45 s; and a final elongation at 72°C for 420 s. Our MHC-I fragment targeted by PCR (240 bp) was examined for strength of amplification and purity on a QIAxcel Advanced System (Qiagen, Germany).

The second 20 μl PCR contained 2 μl TS-PCR product, 10 μl AmpliGold 360 MasterMix, 4 μl unique Fluidigm barcode primers, and 4 μl ultrapure water. Apart from a reduction to 10 amplification cycles and annealing at 60°C, the PCR protocol was equal to the first PCR. Distinct barcodes were also assigned to samples serving as replicates (M. natalensis: 28 out of 204 individuals; M. erythroleucus: 8 out of 69 individuals) (S2 Appendix). The barcoded PCR product was cleaned by using magnetic beads (NucleoMag NGS Clean-up and Size Select, Macherey-Nagel, Germany) and quality checked on the QIAxcel Advanced System. Concentrations for each amplicon were determined with QuantiFlour dsDNA System (Promega, USA) on a microplate reader (Tecan, Switzerland), and equimolar DNA amounts were pooled for the final library. A paired-end sequencing run was performed on an Illumina Miseq using a V2 reagent kit (2x250 cycles).

2.4 MHC allele calling

Our Illumina run produced 8,299,036 raw reads for M. natalensis and 3,082,197 for M. erythroleucus (S1 Appendix). Sequences were analyzed using the “Allele Calling Procedure for Illumina Amplicon Sequencing” pipeline (ACACIA; [55,56]) running both species separately. Essentially, ACACIA designates alleles from raw reads while removing sequences that represent chimeras, reading errors or singletons. In brief, we used the default settings provided by the pipeline for quality control, forward and reverse read merging, primer filtering and discarding of artifacts (singletons and chimeras). The remaining sequences were blasted against a reference database consisting of 754 MHC-I sequences (principally exon 2, with some sequences including exon 3) from 16 murid species extracted from NCBI (www.ncbi.nih.gov), resulting in 3,986,903 reads for M. natalensis and 1,465,859 for M. erythroleucus. For allele calling, we increased the entropy level from 0.2 to 0.35 to deal with extremely high genetic diversity and set the minimum number of reads for retaining a sequence as an MHC allele in an individual to 20. We reached a replicability of 93.83% and 99.12% for the 28 and 8 replicate samples from M. natalensis and M. erythroleucus, respectively.

Based on the ACACIA results, and in keeping with our conservative approach to oligotyping, we excluded nine individuals of M. natalensis and two of M. erythroleucus with less than 15,000 raw reads generated in the Illumina sequencing. We removed an additional five M. natalensis, which we considered to have an abnormally low number of alleles per individual (2–8), in light of the extent of MHC-I loci detected for this species (up to 22). Ultimately, 189 individuals of M. natalensis and 67 M. erythroleucus were retained for subsequent analyses (S1 Appendix). Details such as the nucleotide sequence and frequency of each MHC-I allele are contained in S2 Appendix.

2.5 MHC-I supertyping

Apart from individual alleles, we sought to investigate whether LASV prevalence was linked to the occurrence of MHC-I supertypes (STs), i.e., alleles grouped according to pathogen-recognition properties [57,58]. Identification of positive selection across codon sites in this study was conducted separately for M. erythroleucus and M. natalensis using CodeML in the PamlX graphical user interface [59]. The analysis detected 16 positively selected sites (PSSs), 12 of which are common to both species (S3 Appendix). MHC-I alleles were then clustered into supertypes (STs) based on five z-values [60] that described the physicochemical binding properties of each amino acid at the 12 shared PSSs using DAPC clustering (DAPC, [61]). Hence, MHC supertypes summarize alleles with putatively similar antigen binding characteristics. The optimal number of STs for our z matrix was suggested by the Bayesian Information Criterion (BIC) curve (S4 Appendix). We retained only alleles that were convincingly assigned (i.e., in at least 7 out of 10 DAPC seed runs) within a particular ST. This resulted in 225 M. erythroleucus and M. natalensis alleles allocated to 24 distinct supertypes (S4 Appendix).

2.6 Statistical analyses

MHC allele and supertype composition between Mastomys species (M. natalensis and M. erythroleucus) and between localities (Abagboro and Ebudin) were compared using a PERMANOVA computed with the adonis() function in the ‘vegan’ package [62, 63]. Significant association between the occurrence of individual MHC alleles, supertypes and active or previous LASV infection was assessed with the cooccur() function from the ‘cooccur’ package [64]. Only alleles and STs present in more than 10% of individuals were included in the co-occurrence analyses. Lastly, we confirmed correlative results from the co-occurrence analysis by employing generalised linear mixed effect models (GLMMs) with presence of active or previous LASV infection as response variable and presence of specific MHC alleles/STs and MHC allele/ST diversity as explanatory variables, controlling for eye lens weight in milligrams (a proxy for rodent age; [46]), sex and sampling date as random effect. We used a binomial GLMM with Logit-link function, with no indication of violation of model assumptions. In cases where the co-occurrence analysis identified alleles that were correlated, GLMMs were only computed for the more frequent allele [65]. To sum up, our criteria for recognizing an allele or ST as convincingly associated with LASV were:

1. Significant association with LASV PCR and/or IgG data in the co-occurrence analyses,

2. having a higher frequency compared with correlated alleles/STs also associated with LASV in the co-occurrence analysis,

3. and being significantly linked to LASV occurrence in the complementary GLMM analysis.

All statistical analyses were performed in R using the R Studio computing environment version 2022.07.2 [66].

3. Results

3.1 MHC-I allele and supertype composition between Mastomys species and localities

A total of 193 nucleotide alleles and 21 supertypes (STs) were obtained from 189 M. natalensis rodents, while 96 alleles and 15 STs were obtained from 67 M. erythroleucus (Table 1). Between both species, the average number of alleles per individual (28.27 ± 6.33 vs. 19.69 ± 4.97) and STs per individual (13.16 ± 2.59 vs. 8.97 ± 1.84) was significantly higher for M. natalensis, with the pattern holding true for both sampling locations. Our results indicate a difference in the number of MHC-I loci between the two species, with M. natalensis suggested to have at least 22 MHC-I loci (maximum number of alleles for an individual = 43), whereas M. erythroleucus was suggested to carry at least 15 loci (maximum number of alleles for an individual = 30). Rarefaction analysis demonstrated that the full extent of allelic diversity for M. natalensis was likely not reached, while the lower sampling effort likely captured the natural allelic diversity in M. erythroleucus (S5 Appendix). The majority of alleles were exclusive to either M. natalensis (166) or M. erythroleucus (69), though 27 sequences were shared between both species (Table 1, Fig 2; and S2 Appendix).

[Figure omitted. See PDF.]

Table 1. Summarised distribution of MHC-I alleles and supertypes in two sympatric Mastomys species (M. natalensis and M. erythroleucus).

Rodents were trapped in two distinct sites in Nigeria: Abagboro, a non-endemic site for Lassa fever, and Ebudin, an endemic site for Lassa fever.

https://doi.org/10.1371/journal.pntd.0011984.t001

[Figure omitted. See PDF.]

Fig 2. Indications of Trans-Species Polymorphism (TSP): A gene tree clade illustrating certain nucleotide sequences that are more similar between than within M. natalensis and M. erythroleucus rodents.

E.g., ManaMHC-I*001 & MaerMHC-I*001. M. natalensis alleles are represented by black squares and M. erythroleucus hollow ones.

https://doi.org/10.1371/journal.pntd.0011984.g002

Since only four M. erythroleucus were captured in Abagboro, variation in MHC-I composition between sites could only be statistically evaluated for M. natalensis. Though 85/193 alleles and 21/22 STs were shared between sites (Table 2), dissemination of these two components of M. natalensis MHC-I composition differed significantly between Abagboro and Ebudin (Fig 3; PERMANOVA, alleles, R = 0.503, P = 0.001; STs, R = 0.2055, P = 0.001; S5 Appendix). Also, Ebudin possessed a higher proportion of private alleles (49%, 83/168) than Abagboro (22.73%, 25/110). Individual MHC allelic diversity, however, was similar between locations (Table 2 and S5 Appendix).

[Figure omitted. See PDF.]

Table 2. Mastomys natalensis MHC-I composition between two localities in Nigeria.

Abagboro is a non-endemic site for Lassa fever, and Ebudin is an endemic site for Lassa fever.

https://doi.org/10.1371/journal.pntd.0011984.t002

[Figure omitted. See PDF.]

Fig 3. MHC-I allele and supertype distribution between Mastomys species and localities.

Alleles *001 to *080 are shown here. Rare alleles (labelled > *080) are displayed in S2 Appendix. Colours represent the distinct localities, Abagboro (in blue) is considered non-endemic for LASV, Ebudin (in yellow) is endemic for LASV. The dotted lines represent the minimum frequency threshold of alleles and supertypes to be included in the co-occurrence analyses.

https://doi.org/10.1371/journal.pntd.0011984.g003

3.2 MHC-I genotype and LASV infection

M. natalensis and M. erythroleucus in Ebudin, the endemic zone for Lassa fever, were detected to carry active or previous LASV infections (Table 3). Conversely, both rodent species in Abagboro, the non-endemic area, were LASV-negative, except for one M. natalensis individual, which was IgG-positive. Because LASV infections were only recorded in Ebudin, our assessment of the relationship between MHC-I constitution and LASV focused on this location. Several alleles were found to be significantly associated with active and previous LASV infection in the co-occurrence analyses (Fig 4 and S6 Appendix). When MHC alleles were correlated with each other, only the allele occurring at higher frequency was investigated further in the GLMMs. Specifically, alleles ManaMHC-I*012, *039 & *104 were correlated with *006, but were less frequent (S2 and S6 Appendices). The pattern was the same for alleles ManaMHC-I*009, *056 & *137 with *008; and allele MaerMHC-I*024 with *MaerMHC-I*008. Furthermore, we only regarded LASV-allele associations from the co-occurrence analyses as reliable when they were supported by GLMM results (Fig 4 and S7 Appendix). Ultimately, alleles that satisfied these stipulations were ManaMHC-I*006, *008, *021 and MaerMHC-I*008 (Fig 4 and S6 and S7 Appendices).

[Figure omitted. See PDF.]

Table 3. LASV infection profile of two Mastomys rodents including the number of individuals that were tested negative for LASV and those tested positive for an active LASV infection in PCRs and/or positive for LASV specific IgG antibodies.

https://doi.org/10.1371/journal.pntd.0011984.t003

[Figure omitted. See PDF.]

Fig 4. Residual co-occurrence values and sequence patterns of Mastomys alleles found to be significantly associated with LASV infection.

A) M. natalensis alleles significantly associated with LASV PCR. B) M. natalensis alleles significantly associated with LASV IgG. C) M. erythroleucus alleles significantly associated with LASV PCR. D) M. erythroleucus alleles significantly associated with LASV IgG. Asterisks denote significant LASV-allele relationships from both the co-occurrence analyses and General Linear Mixed Effect Models (GLMMs). Mana*009, *056 and *137 were correlated with *008; Mana*012, *039 and *104 were correlated with *006; and Maer*024 was correlated with Maer*008; In this case, GLMMs were only computed for the more frequent alleles (see Methods). E) Amino acid sequences of alleles significantly associated with LASV occurrence in both the co-occurrence analyses and GLMMs, showing positively selected sites common to both M. natalensis (n) and M. erythroleucus (e).

https://doi.org/10.1371/journal.pntd.0011984.g004

Allele ManaMHC-I*006 was negatively associated with active (PCR-positive) LASV infection in M. natalensis, suggesting resistance (Fig 4A). The same allele was positively associated with previous infection (indicated by an IgG-positive result), possibly suggesting competent clearance (Fig 4B). The alleles ManaMHC-I*008 & *021 were positively associated with active infection, suggesting susceptibility (Fig 4A). In the case of M. erythroleucus, MaerMHC-I*008 was positively associated with active LASV infection (Fig 4C). Associations between STs and LASV infection status detected by co-occurrence analyses were not supported by GLMMs (S6 and S7 Appendices). Apart from the alleles above, no fixed effects or control variables (except eye lens weight, which was positively linked to previous LASV infection) were significantly associated with active or previous LASV infections (S7 Appendix).

3.3 Amino acid differences associated with LASV infection

Among the PSSs of the alleles linked to LASV infections (Fig 4), the three susceptibility alleles shared a Glutamic acid at position 57, whereas the resistance allele ManaMHC-I*006 presented an Arginine at the same position. A similar motif was observed at three sequence positions outside the PSS: ManaMHC-I*006 possessed residue Glutamine at positions 26 and 41, while the susceptibility alleles coded Glutamic acid and Proline, respectively. ManaMHC-I*006 also showed Lysine and the susceptibility alleles Arginine at position 74.

Bearing in mind the trans-species polymorphism hinted at in this study (Fig 2 and S2 Appendix), it is important to note that ManaMHC-I*006 & ManaMHC-I*021 are not private alleles, but have counterparts in M. erythroleucus (MaerMHC-I*020 & MaerMHC-I*062, respectively). However, the M. erythroleucus alleles were not significantly associated with LASV.

Even though we cannot directly test whether differences in allele and/or ST composition may contribute to differences in LASV prevalence between sites, it is noteworthy that ManaMHC-I*006, ManaMHC-I*008 & ManaMHC-I*021 were relatively abundant in Ebudin, where they are associated with LASV. These alleles also occurred with a high frequency in Abagboro, where LASV was not detected (Fig 3).

4. Discussion

Considering the importance of African multimammate mice as preeminent reservoirs of several actual and potential zoonotic pathogens that include mammarenaviruses, very little is known about the link to the immunogenetics of these animals [39]. This is key to gauge the likelihood for any such virus, such as LASV, to be transmitted between individual rodents, to emerge across populations and to spillover to humans. This study is the first to describe MHC-I sequence characteristics, immunogenetic diversity, and distribution of MHC-I alleles and supertypes between Mastomys natalensis and M. erythroleucus species, and between populations from the endemic and non-endemic zone for Lassa fever. Importantly, we demonstrated that specific MHC-I alleles in Mastomys rodents within Nigeria are linked to resistance, competent clearance, and susceptibility to LASV.

4.1 MHC-I constitution in Mastomys rodents

Information about diversity of the MHC class I in wild murid rodents is still rather limited. Wild-caught Mus musculus across various sites in Europe and Iran yielded 29 and 183 MHC-I alleles from exon 2 of the H2-D and H2-K loci, respectively [53], compared to the 30–40 genes found in conventional laboratory Mus musculus. Wild African Pygmy mice Mus (Nannomys) setulosus were even found to harbour several thousand MHC-I genes [67]. Our study offers a first glimpse into similar diversity of the MHC-I in Mastomys rodents: 193 and 96 alleles were detected in M. natalensis and M. erythroleucus, respectively. Moreover, at least 22 MHC-I loci were estimated in M. natalensis, and 15 in M. erythroleucus. Such exceptional diversity is emblematic of the MHC and commonly coupled to gene duplication [68–70]. The allelic MHC-I variability observed between M. natalensis and M. erythroleucus is not just a result of larger sample size, but more likely due to M. natalensis’ wide geographic distribution, its ecological adaptability and well-documented interaction with a spectrum of mammarenaviruses and pathogens beyond that of any other African rodent. A similarly expansive allelic MHC-I repertoire has been reported for birds and bats with a particularly wide range [71,72]. The fact that M. natalensis and M. erythroleucus share 27 alleles implies trans-species polymorphism possibly maintained by balancing selection acting from a million years ago, when both taxa are estimated to have split from their common ancestor [73]. Furthermore, 12 of 16 positively selected sites are shared by M. natalensis and M. erythroleucus, which supports the idea that both species experience similar pathogen-mediated selection.

4.2 Immunogenetic correlates of acute and cleared LASV infection

The majority of LASV and similar mammarenavirus infections in Mastomys are acute [21,44,46] and last between 2–4 weeks [45,47]. IgG antibodies commonly appear in serum 1–2 weeks after infection and are retained possibly for life [45,47]. Field data for various mammarenaviruses show that M. natalensis which simultaneously test positive by PCR and IgG assays (i.e., those that carry a chronic infection from shortly after birth, or are seroconverting) represent a minor proportion of infected individuals [21,44,46]. Our findings line up with these earlier results, as the great majority of individuals were either PCR or IgG positive, with only very few exceptions. It begs the question whether immunogenetic differences are at the root of differences in the speed of progression of virus infection as well as the swiftness of IgG responses. In M. natalensis, allele ManaMHC-I*006 co-occurred negatively with the LASV PCR- and positively with the IgG-signal. We propose that individuals bearing this allele, when challenged with LASV, clear the virus efficiently. These rodents are therefore less likely to be encountered with viremia, but more likely with mounted IgG responses. By contrast, positive associations between LASV-PCR and ManaMHC-I*008, ManaMHC-I*021 as well as MaerMHC-I*008 implied individuals were more likely to be captured with viremia, possibly due to inability to promptly clear the virus and develop protective IgG antibodies.

A mark-recapture study involving the dynamics of Morogoro virus (MORV, considered a LASV surrogate) in a natural M. natalensis population demonstrated that 15 out of 57 individuals, which tested PCR-positive with an active infection, had not turned IgG-positive upon recapture by (or, in some cases, beyond) the second week [48]. The authors hypothesize that IgG antibody development might have been delayed or was expressed in titres below the detection threshold. Our work offers a plausible explanation: the identity of MHC alleles might determine the speed at which an infection is cleared while concurrently affecting the IgG antibody response. For instance, White Leghorn chickens (Gallus gallus domesticus) that possessed the MHC haplotype B13 showed relatively high antibody titres in a lineage normally known to display low antibody response to Brucella abortus [74]. Likewise, minimum and maximum antibody titres to Salmonella enteritidis, sheep red blood cell and Brucella abortus are connected to MHC class I & II SNPs in chickens [75]. The association of the allele ManaMHC-I*006 with both acute infection likelihood and antibody detection points towards a mechanistic link between IgG production and MHC constitution. Future research should explore whether a link exists between MHC genotype, the timing of IgG antibody development post-infection and the amplitude of titre response.

4.3 Does Mastomys MHC-I explain LASV endemicity?

Can we explain differences in LASV prevalence between the endemic zone of Ebudin and the non-endemic area of Abagboro based on immunogenetic disparities between sites? First, Mastomys natalensis did not differ significantly in average number of alleles or supertypes per individual between Ebudin and Abagboro, nor was there a link between LASV and MHC diversity. Yet, there was a higher proportion of private alleles in Ebudin than Abagboro. In this sense, the endemic site contains a larger pool of rare alleles originating from diversifying selection, likely maintained by a cyclical “arms race” [30] between host immunogenetic defenses and evolving LASV lineage II variants (and/or many other pathogens) [51]. The LASV-associated alleles (ManaMHC-I*006, *008 & 021) are among the most abundant in M. natalensis (S2 Appendix). Intriguingly, however, these alleles generally registered lower percentages in Ebudin (43/113, 38%; 39/113, 35%; 37/113, 33%, respectively) than they did in Abagboro (39/76, 51%; 34/76, 45%; 26/76, 34%, respectively). Although our data cannot more than hint at possible frequency-dependent selection mechanisms in this population, changes in MHC allele frequency across space and time can be rapid in short-lived mammals [76,77] such as Mastomys rodents, which live to about 339 days [78]. In concert with more extensive sampling in non-endemic sites or transitional localities, extending our surveys in these populations by several more Mastomys generations could improve insight into how the frequency of MHC alleles fluctuates with LASV prevalence. Also, the study of T-cell lymphocytes, specially TCD8+, would be very desirable for a more complete picture of the role of immunogenetics in LASV infections.

5. Conclusion

We were able to identify MHC-I allele associations with LASV that conferred resistance and competent clearance and, contrariwise, susceptibility in the most important rodent reservoirs, M. natalensis and M. erythroleucus. The prevalence of mammarenaviruses in Mastomys populations are immediate public health concerns. Hence, expanding our surveys spatially and temporally will provide a broader view of the varied MHC-I alleles related to LASV occurrence, and their frequency distribution across Mastomys populations within and outside Nigeria. A better understanding of the ecology and evolution of the host, particularly with respect to interactions with humans in a populous country such as Nigeria, might help to identify spillover hotspots and predict emergence, as was done in bat-borne viruses for example [79]. Lastly, our results open the door to future research on MHC-I allele-LASV peptide combinations that advance the development of a LASV vaccine for rodents in order to regulate LASV prevalence in its most important reservoir [80,81].

Supporting information

S1 Appendix. Data obtained by illumina high-thoughput sequencing.

https://doi.org/10.1371/journal.pntd.0011984.s001

(DOCX)

S2 Appendix. MHC diversity in Mastomys.

https://doi.org/10.1371/journal.pntd.0011984.s002

(XLSX)

S3 Appendix. Mastomys MHC-I exon 2 sequence attributes.

https://doi.org/10.1371/journal.pntd.0011984.s003

(DOCX)

S4 Appendix. Allocation of MHC I alleles into supertypes.

https://doi.org/10.1371/journal.pntd.0011984.s004

(XLSX)

S5 Appendix. Statistical tests comparing MHC-I composition between Mastomys species and between sampling localities.

https://doi.org/10.1371/journal.pntd.0011984.s005

(DOCX)

S6 Appendix. Co-occurrence analyses of LASV and specific MHC alleles/supertypes in Mastomys.

https://doi.org/10.1371/journal.pntd.0011984.s006

(DOCX)

S7 Appendix. Generalized linear mixed effect models of LASV and specific MHC alleles/supertypes asssociations in Mastomys.

https://doi.org/10.1371/journal.pntd.0011984.s007

(DOCX)

Acknowledgments

We are grateful to U. Stehle for laboratory support.

References

1. 1. Basinski AJ, Fichet-Calvet E, Sjodin AR, Varrelman TJ, Remien CH, Layman NC, et al. Bridging the gap: Using reservoir ecology and human serosurveys to estimate Lassa virus spillover in West Africa. PLoS Comput Biol. 2021;17(3):e1008811. Epub 2021/03/04. pmid:33657095

* View Article

* PubMed/NCBI

* Google Scholar

2. 2. Keane E, Gilles HM. Lassa fever in Panguma Hospital, Sierra Leone, 1973–6. Br Med J. 1977;1(6073):1399–402. Epub 1977/05/28. pmid:861652

* View Article

* PubMed/NCBI

* Google Scholar

3. 3. Price ME, Fisher-Hoch SP, Craven RB, McCormick JB. A prospective study of maternal and fetal outcome in acute Lassa fever infection during pregnancy. BMJ. 1988;297(6648):584–7. Epub 1988/09/03. pmid:3139220

* View Article

* PubMed/NCBI

* Google Scholar

4. 4. Shaffer JG, Grant DS, Schieffelin JS, Boisen ML, Goba A, Hartnett JN, et al. Lassa fever in post-conflict sierra leone. PLoS Negl Trop Dis. 2014;8(3):e2748. Epub 2014/03/22. pmid:24651047 Hemorrhagic Fever Consortium (vhfc.org) is a partnership of academic and industry scientists who are developing diagnostics, therapeutics and vaccines for LF and other severe diseases. Tulane University and various industry partners have filed United States and foreign patent applications on behalf of the consortium for several of these technologies. If commercial products are developed, consortium members may receive royalties or profits. This does not alter our adherence to all PLOS policies on sharing data and materials.

* View Article

* PubMed/NCBI

* Google Scholar

5. 5. Mateer EJ, Huang C, Shehu NY, Paessler S. Lassa fever-induced sensorineural hearing loss: A neglected public health and social burden. PLoS Negl Trop Dis. 2018;12(2):e0006187. Epub 2018/02/23. pmid:29470486

* View Article

* PubMed/NCBI

* Google Scholar

6. 6. Arruda LB, Haider N, Olayemi A, Simons D, Ehichioya D, Yinka-Ogunleye A, et al. The niche of One Health approaches in Lassa fever surveillance and control. Ann Clin Microbiol Antimicrob. 2021;20(1):29. Epub 2021/04/26. pmid:33894784

* View Article

* PubMed/NCBI

* Google Scholar

7. 7. Grange ZL, Goldstein T, Johnson CK, Anthony S, Gilardi K, Daszak P, et al. Ranking the risk of animal-to-human spillover for newly discovered viruses. Proc Natl Acad Sci U S A. 2021;118(15). Epub 2021/04/07. pmid:33822740

* View Article

* PubMed/NCBI

* Google Scholar

8. 8. Mylne AQ, Pigott DM, Longbottom J, Shearer F, Duda KA, Messina JP, et al. Mapping the zoonotic niche of Lassa fever in Africa. Trans R Soc Trop Med Hyg. 2015;109(8):483–92. Epub 2015/06/19. pmid:26085474

* View Article

* PubMed/NCBI

* Google Scholar

9. 9. WHO. Annual Review of Diseases Prioritized under the Research and Development Blueprint Informal Consultation 24–25 January 2017. Geneva, Switzerland.. 2017;http://www.who.int/blueprint/what/research-development/2017-Prioritization-Long-Report.pdf?ua=1.

10. 10. Mehand MS, Al-Shorbaji F, Millett P, Murgue B. The WHO R&D Blueprint: 2018 review of emerging infectious diseases requiring urgent research and development efforts. Antiviral Res. 2018;159:63–7. Epub 2018/09/28. pmid:30261226

* View Article

* PubMed/NCBI

* Google Scholar

11. 11. Garry RF. Lassa fever—the road ahead. Nat Rev Microbiol. 2022:1–10. Epub 2022/09/14. pmid:36097163 countermeasures for Lassa virus (LASV) and other emerging viruses.

* View Article

* PubMed/NCBI

* Google Scholar

12. 12. Monath TP, Newhouse VF, Kemp GE, Setzer HW, Cacciapuoti A. Lassa virus isolation from Mastomys natalensis rodents during an epidemic in Sierra Leone. Science. 1974;185(4147):263–5. Epub 1974/07/19. pmid:4833828.

* View Article

* PubMed/NCBI

* Google Scholar

13. 13. Marien J, Borremans B, Gryseels S, Soropogui B, De Bruyn L, Bongo GN, et al. No measurable adverse effects of Lassa, Morogoro and Gairo arenaviruses on their rodent reservoir host in natural conditions. Parasit Vectors. 2017;10(1):210. Epub 2017/04/30. pmid:28449693

* View Article

* PubMed/NCBI

* Google Scholar

14. 14. Wood R, Bangura U, Marien J, Douno M, Fichet-Calvet E. Detection of Lassa virus in wild rodent feces: Implications for Lassa fever burden within households in the endemic region of Faranah, Guinea. One Health. 2021;13:100317. Epub 2021/09/16. pmid:34522759

* View Article

* PubMed/NCBI

* Google Scholar

15. 15. Marien J, Lo Iacono G, Rieger T, Magassouba N, Gunther S, Fichet-Calvet E. Households as hotspots of Lassa fever? Assessing the spatial distribution of Lassa virus-infected rodents in rural villages of Guinea. Emerg Microbes Infect. 2020;9(1):1055–64. Epub 2020/05/28. pmid:32459576

* View Article

* PubMed/NCBI

* Google Scholar

16. 16. Happold DCD. Mammals of Africa Volume III—Rodents, Hares and Rabbits. Bloomsbury, London. 2013:784 pp.

17. 17. Gouy de Bellocq JB, A; Martynov A. A.; Lavrenchenko L. A. Dhati Welel virus, the missing mammarenavirus of the widespread Mastomys natalensis. Journal of Vertebrate Biology. 2020;69(20018):1–11.

* View Article

* Google Scholar

18. 18. Gryseels S, Rieger T, Oestereich L, Cuypers B, Borremans B, Makundi R, et al. Gairo virus, a novel arenavirus of the widespread Mastomys natalensis: Genetically divergent, but ecologically similar to Lassa and Morogoro viruses. Virology. 2015;476:249–56. Epub 2015/01/07. pmid:25559385.

* View Article

* PubMed/NCBI

* Google Scholar

19. 19. Ishii A, Thomas Y, Moonga L, Nakamura I, Ohnuma A, Hang’ombe B, et al. Novel arenavirus, Zambia. Emerg Infect Dis. 2011;17(10):1921–4. Epub 2011/10/18. pmid:22000372

* View Article

* PubMed/NCBI

* Google Scholar

20. 20. Olayemi A, Obadare A, Oyeyiola A, Igbokwe J, Fasogbon A, Igbahenah F, et al. Arenavirus Diversity and Phylogeography of Mastomys natalensis Rodents, Nigeria. Emerg Infect Dis. 2016;22(4):694–7. Epub 2016/03/18. pmid:26982388

* View Article

* PubMed/NCBI

* Google Scholar

21. 21. Borremans B, Leirs H, Gryseels S, Gunther S, Makundi R, de Bellocq JG. Presence of Mopeia virus, an African arenavirus, related to biotope and individual rodent host characteristics: implications for virus transmission. Vector Borne Zoonotic Dis. 2011;11(8):1125–31. Epub 2010/12/15. pmid:21142956.

* View Article

* PubMed/NCBI

* Google Scholar

22. 22. Gunther S, Hoofd G, Charrel R, Roser C, Becker-Ziaja B, Lloyd G, et al. Mopeia virus-related arenavirus in natal multimammate mice, Morogoro, Tanzania. Emerg Infect Dis. 2009;15(12):2008–12. Epub 2009/12/08. pmid:19961688

* View Article

* PubMed/NCBI

* Google Scholar

23. 23. Haikukutu L, Lyaku JR, Lyimo C, Kasanga CJ, Kandusi SE, Rahelinirina S, et al. Plague in Tanzania: first report of sylvatic plague in Morogoro region, persistence in Mbulu focus, and ongoing quiescence in Lushoto and Iringa foci. IJID Reg. 2022;4:105–10. Epub 2022/07/27. pmid:35880003

* View Article

* PubMed/NCBI

* Google Scholar

24. 24. Thomas CMV.; V-Houtte N.; Mhamphi G.; Mariën J.; Sabuni C.; Makundi I.; Nzalawahe J.; Machang’u R.; Leirs H. Molecular detection of Leptospira and Bartonella in Mastomys natalensis and its ectoparasites in Morogoro, Tanzania. Mammalia. 2023;8(5):539–44.

* View Article

* Google Scholar

25. 25. Ortiz PA, Garcia HA, Lima L, da Silva FM, Campaner M, Pereira CL, et al. Diagnosis and genetic analysis of the worldwide distributed Rattus-borne Trypanosoma (Herpetosoma) lewisi and its allied species in blood and fleas of rodents. Infect Genet Evol. 2018;63:380–90. Epub 2017/09/09. pmid:28882517.

* View Article

* PubMed/NCBI

* Google Scholar

26. 26. Olayemi A, Cadar D, Magassouba N, Obadare A, Kourouma F, Oyeyiola A, et al. New Hosts of The Lassa Virus. Sci Rep. 2016;6:25280. Epub 2016/05/04. pmid:27140942

* View Article

* PubMed/NCBI

* Google Scholar

27. 27. Olayemi A, Fichet-Calvet E. Systematics, Ecology, and Host Switching: Attributes Affecting Emergence of the Lassa Virus in Rodents across Western Africa. Viruses. 2020;12(3). Epub 2020/03/19. pmid:32183319

* View Article

* PubMed/NCBI

* Google Scholar

28. 28. Apanius V, Penn D, Slev PR, Ruff LR, Potts WK. The nature of selection on the major histocompatibility complex. Crit Rev Immunol. 1997;17(2):179–224. Epub 1997/01/01. pmid:9094452.

* View Article

* PubMed/NCBI

* Google Scholar

29. 29. Janeway CA, Travers P., Walport M., & Sclomchik M. J. Immunobiology: the immune system in health and disease. New York: Garland Science Publishing. 2005.

30. 30. Kaufman J. Generalists and Specialists: A New View of How MHC Class I Molecules Fight Infectious Pathogens. Trends Immunol. 2018;39(5):367–79. Epub 2018/02/06. pmid:29396014

* View Article

* PubMed/NCBI

* Google Scholar

31. 31. Radwan J, Babik W, Kaufman J, Lenz TL, Winternitz J. Advances in the Evolutionary Understanding of MHC Polymorphism. Trends Genet. 2020;36(4):298–311. Epub 2020/02/12. pmid:32044115.

* View Article

* PubMed/NCBI

* Google Scholar

32. 32. Sommer S. The importance of immune gene variability (MHC) in evolutionary ecology and conservation. Front Zool. 2005;2:16. Epub 2005/10/26. pmid:16242022

* View Article

* PubMed/NCBI

* Google Scholar

33. 33. Hawley DM, Fleischer RC. Contrasting epidemic histories reveal pathogen-mediated balancing selection on class II MHC diversity in a wild songbird. PLoS One. 2012;7(1):e30222. Epub 2012/02/01. pmid:22291920

* View Article

* PubMed/NCBI

* Google Scholar

34. 34. Savage AE, Zamudio KR. MHC genotypes associate with resistance to a frog-killing fungus. Proc Natl Acad Sci U S A. 2011;108(40):16705–10. Epub 2011/09/29. pmid:21949385

* View Article

* PubMed/NCBI

* Google Scholar

35. 35. Wegner KM, Kalbe M, Kurtz J, Reusch TB, Milinski M. Parasite selection for immunogenetic optimality. Science. 2003;301(5638):1343. Epub 2003/09/06. pmid:12958352.

* View Article

* PubMed/NCBI

* Google Scholar

36. 36. Froeschke G, Sommer S. MHC class II DRB variability and parasite load in the striped mouse (Rhabdomys pumilio) in the Southern Kalahari. Mol Biol Evol. 2005;22(5):1254–9. Epub 2005/02/11. pmid:15703235.

* View Article

* PubMed/NCBI

* Google Scholar

37. 37. Lenz TL, Wells K, Pfeiffer M, Sommer S. Diverse MHC IIB allele repertoire increases parasite resistance and body condition in the Long-tailed giant rat (Leopoldamys sabanus). BMC Evol Biol. 2009;9:269. Epub 2009/11/26. pmid:19930637

* View Article

* PubMed/NCBI

* Google Scholar

38. 38. Stewart PW, Chapes SK. Role of major histocompatibility complex class II in resistance of mice to naturally acquired infection with Syphacia obvelata. Comp Med. 2003;53(1):70–4. Epub 2003/03/11. pmid:12625509.

* View Article

* PubMed/NCBI

* Google Scholar

39. 39. Gouy de Bellocq J, Leirs H. Polymorphism and signatures of selection in the multimammate rat DQB gene. Immunogenetics. 2010;62(1):59–64. Epub 2009/12/03. pmid:19953242.

* View Article

* PubMed/NCBI

* Google Scholar

40. 40. Hardin A, Nevonen KA, Eckalbar WL, Carbone L, Ahituv N. Comparative Genomic Characterization of the Multimammate Mouse Mastomys coucha. Mol Biol Evol. 2019;36(12):2805–12. Epub 2019/08/20. pmid:31424545

* View Article

* PubMed/NCBI

* Google Scholar

41. 41. Botten J, Alexander J, Pasquetto V, Sidney J, Barrowman P, Ting J, et al. Identification of protective Lassa virus epitopes that are restricted by HLA-A2. J Virol. 2006;80(17):8351–61. Epub 2006/08/17. pmid:16912286

* View Article

* PubMed/NCBI

* Google Scholar

42. 42. Gibb R, Redding DW, Chin KQ, Donnelly CA, Blackburn TM, Newbold T, et al. Zoonotic host diversity increases in human-dominated ecosystems. Nature. 2020;584(7821):398–402. Epub 2020/08/08. pmid:32759999.

* View Article

* PubMed/NCBI

* Google Scholar

43. 43. Johnson CK, Hitchens PL, Pandit PS, Rushmore J, Evans TS, Young CCW, et al. Global shifts in mammalian population trends reveal key predictors of virus spillover risk. Proc Biol Sci. 2020;287(1924):20192736. Epub 2020/04/08. pmid:32259475

* View Article

* PubMed/NCBI

* Google Scholar

44. 44. Fichet-Calvet E, Becker-Ziaja B, Koivogui L, Gunther S. Lassa serology in natural populations of rodents and horizontal transmission. Vector Borne Zoonotic Dis. 2014;14(9):665–74. Epub 2014/09/18. pmid:25229705

* View Article

* PubMed/NCBI

* Google Scholar

45. 45. Borremans B, Vossen R, Becker-Ziaja B, Gryseels S, Hughes N, Van Gestel M, et al. Shedding dynamics of Morogoro virus, an African arenavirus closely related to Lassa virus, in its natural reservoir host Mastomys natalensis. Sci Rep. 2015;5:10445. Epub 2015/05/30. pmid:26022445

* View Article

* PubMed/NCBI

* Google Scholar

46. 46. Olayemi A, Oyeyiola A, Obadare A, Igbokwe J, Adesina AS, Onwe F, et al. Widespread arenavirus occurrence and seroprevalence in small mammals, Nigeria. Parasit Vectors. 2018;11(1):416. Epub 2018/07/15. pmid:30005641

* View Article

* PubMed/NCBI

* Google Scholar

47. 47. Hoffmann C, Wurr S, Pallasch E, Bockholt S, Rieger T, Gunther S, et al. Experimental Morogoro Virus Infection in Its Natural Host, Mastomys natalensis. Viruses. 2021;13(5). Epub 2021/06/03. pmid:34067011

* View Article

* PubMed/NCBI

* Google Scholar

48. 48. Marien J, Borremans B, Gryseels S, Broecke BV, Becker-Ziaja B, Makundi R, et al. Arenavirus Dynamics in Experimentally and Naturally Infected Rodents. Ecohealth. 2017;14(3):463–73. Epub 2017/06/16. pmid:28616660.

* View Article

* PubMed/NCBI

* Google Scholar

49. 49. Westerdahl H, Asghar M, Hasselquist D, Bensch S. Quantitative disease resistance: to better understand parasite-mediated selection on major histocompatibility complex. Proc Biol Sci. 2012;279(1728):577–84. Epub 2011/07/08. pmid:21733902

* View Article

* PubMed/NCBI

* Google Scholar

50. 50. Mills JNY T. L.; Childs J.; Parmenter R. R.; Ksiazek T. G.; Rollin P. E; et al. Guidelines for working with rodents potentially infected with hantavirus. J Mammal. 1995;76:716–22.

* View Article

* Google Scholar

51. 51. Adesina AS, Oyeyiola A, Obadare A, Igbokwe J, Abejegah C, Akhilomen P, et al. Circulation of Lassa virus across the endemic Edo-Ondo axis, Nigeria, with cross-species transmission between multimammate mice. Emerg Microbes Infect. 2023;12(1):2219350. Epub 2023/06/08. pmid:37288752

* View Article

* PubMed/NCBI

* Google Scholar

52. 52. Wulff H, Lange JV. Indirect immunofluorescence for the diagnosis of Lassa fever infection. Bull World Health Organ. 1975;52(4–6):429–36. Epub 1975/01/01. pmid:821624

* View Article

* PubMed/NCBI

* Google Scholar

53. 53. Linnenbrink M, Teschke M, Montero I, Vallier M, Tautz D. Meta-populational demes constitute a reservoir for large MHC allele diversity in wild house mice (Mus musculus). Front Zool. 2018;15:15. Epub 2018/05/04. pmid:29721030

* View Article

* PubMed/NCBI

* Google Scholar

54. 54. Rozen S, Skaletsky H. Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000;132:365–86. Epub 1999/11/05. pmid:10547847.

* View Article

* PubMed/NCBI

* Google Scholar

55. 55. Gillingham MAF, Montero BK, Wihelm K, Grudzus K, Sommer S, Santos PSC. A novel workflow to improve genotyping of multigene families in wildlife species: An experimental set-up with a known model system. Mol Ecol Resour. 2021;21(3):982–98. Epub 2020/10/29. pmid:33113273.

* View Article

* PubMed/NCBI

* Google Scholar

56. 56. Fleischer R, Schmid DW, Wasimuddin , Brandel SD, Rasche A, Corman VM, et al. Interaction between MHC diversity and constitution, gut microbiota and Astrovirus infections in a neotropical bat. Mol Ecol. 2022;31(12):3342–59. Epub 2022/05/06. pmid:35510794.

* View Article

* PubMed/NCBI

* Google Scholar

57. 57. Sepil I, Lachish S, Hinks AE, Sheldon BC. Mhc supertypes confer both qualitative and quantitative resistance to avian malaria infections in a wild bird population. Proc Biol Sci. 2013;280(1759):20130134. Epub 2013/03/22. pmid:23516242

* View Article

* PubMed/NCBI

* Google Scholar

58. 58. Sidney J, Grey HM, Kubo RT, Sette A. Practical, biochemical and evolutionary implications of the discovery of HLA class I supermotifs. Immunol Today. 1996;17(6):261–6. Epub 1996/06/01. pmid:8962628.

* View Article

* PubMed/NCBI

* Google Scholar

59. 59. Xu B, Yang Z. PAMLX: a graphical user interface for PAML. Mol Biol Evol. 2013;30(12):2723–4. Epub 2013/10/10. pmid:24105918.

* View Article

* PubMed/NCBI

* Google Scholar

60. 60. Sandberg M, Eriksson L, Jonsson J, Sjostrom M, Wold S. New chemical descriptors relevant for the design of biologically active peptides. A multivariate characterization of 87 amino acids. J Med Chem. 1998;41(14):2481–91. Epub 1998/07/03. pmid:9651153.

* View Article

* PubMed/NCBI

* Google Scholar

61. 61. Jombart T, Devillard S, Balloux F. Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet. 2010;11:94. Epub 2010/10/19. pmid:20950446

* View Article

* PubMed/NCBI

* Google Scholar

62. 62. Anderson MJ. Permutational Multivariate Analysis of Variance. WileyStatsRef: Statistics Reference Online. 2017.

63. 63. Oksanen J, Simpson, G. L., Blanchet, G. F., Roeland Kindt, Pierre, Legendre PRM, R.B. O’Hara, Peter Solymos, M. Henry H. Stevens, Eduard Szoecs HW, Matt Barbour, Michael Bedward, Ben Bolker, Daniel Borcard GC, Michael Chirico, Miquel De Caceres, Sebastien Durand HBAE, Rich FitzJohn, Michael Friendly BF, Geoffrey Hannigan, Mark O. Hill, Leo, et al. vegan: Community Ecology Package. R package version 2.6–2. 2022.

64. 64. Griffith DMM, C. J.;. cooccur: Probabilistic Species Co.Occurrence Analysis in R. Journal of Statistical Software. 2016;131664617.

65. 65. Raberg L, Clough D, Hagstrom A, Scherman K, Andersson M, Drews A, et al. MHC class II genotype-by-pathogen genotype interaction for infection prevalence in a natural rodent-Borrelia system. Evolution. 2022;76(9):2067–75. Epub 2022/08/02. pmid:35909235

* View Article

* PubMed/NCBI

* Google Scholar

66. 66. R-Core-Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2022.

67. 67. Delarbre C, Jaulin C, Kourilsky P, Gachelin G. Evolution of the major histocompatibility complex: a hundred-fold amplification of MHC class I genes in the African pigmy mouse Nannomys setulosus. Immunogenetics. 1992;37(1):29–38. Epub 1992/01/01. pmid:1428009.

* View Article

* PubMed/NCBI

* Google Scholar

68. 68. Axtner J, Sommer S. Gene duplication, allelic diversity, selection processes and adaptive value of MHC class II DRB genes of the bank vole, Clethrionomys glareolus. Immunogenetics. 2007;59(5):417–26. Epub 2007/03/14. pmid:17351770.

* View Article

* PubMed/NCBI

* Google Scholar

69. 69. Miller HC, Lambert DM. Gene duplication and gene conversion in class II MHC genes of New Zealand robins (Petroicidae). Immunogenetics. 2004;56(3):178–91. Epub 2004/05/13. pmid:15138734.

* View Article

* PubMed/NCBI

* Google Scholar

70. 70. Reusch TB, Schaschl H, Wegner KM. Recent duplication and inter-locus gene conversion in major histocompatibility class II genes in a teleost, the three-spined stickleback. Immunogenetics. 2004;56(6):427–37. Epub 2004/08/24. pmid:15322775.

* View Article

* PubMed/NCBI

* Google Scholar

71. 71. Biedrzycka A, O’Connor E, Sebastian A, Migalska M, Radwan J, Zajac T, et al. Extreme MHC class I diversity in the sedge warbler (Acrocephalus schoenobaenus); selection patterns and allelic divergence suggest that different genes have different functions. BMC Evol Biol. 2017;17(1):159. Epub 2017/07/07. pmid:28679358

* View Article

* PubMed/NCBI

* Google Scholar

72. 72. Qurkhuli T, Schwensow N, Brandel SD, Tschapka M, Sommer S. Can extreme MHC class I diversity be a feature of a wide geographic range? The example of Seba’s short-tailed bat (Carollia perspicillata). Immunogenetics. 2019;71(8–9):575–87. Epub 2019/09/15. pmid:31520134

* View Article

* PubMed/NCBI

* Google Scholar

73. 73. Hánová AK A.; Mikula O.; Bryjová A., Šumbera R.; Bryja J. Diversity, distribution, and evolutionary history of the most studied African rodents, multimammate mice of the genus Mastomys: An overview after a quarter of century of using DNA sequencing. Journal of Systematics and Evolutionary Research. 2021;59(8):2500–18.

* View Article

* Google Scholar

74. 74. Dunnington EA, Larsen CT, Gross WB, Siegel PB. Antibody responses to combinations of antigens in white Leghorn chickens of different background genomes and major histocompatibility complex genotypes. Poult Sci. 1992;71(11):1801–6. Epub 1992/11/01. pmid:1437966.

* View Article

* PubMed/NCBI

* Google Scholar

75. 75. Zhou H, Lamont SJ. Chicken MHC class I and II gene effects on antibody response kinetics in adult chickens. Immunogenetics. 2003;55(3):133–40. Epub 2003/05/14. pmid:12743657.

* View Article

* PubMed/NCBI

* Google Scholar

76. 76. Migalska M, Przesmycka K, Alsarraf M, Bajer A, Behnke-Borowczyk J, Grzybek M, et al. Long term patterns of association between MHC and helminth burdens in the bank vole support Red Queen dynamics. Mol Ecol. 2022;31(12):3400–15. Epub 2022/05/06. pmid:35510766

* View Article

* PubMed/NCBI

* Google Scholar

77. 77. Oliver MK, Lambin X, Cornulier T, Piertney SB. Spatio-temporal variation in the strength and mode of selection acting on major histocompatibility complex diversity in water vole (Arvicola terrestris) metapopulations. Mol Ecol. 2009;18(1):80–92. Epub 2009/01/15. pmid:19140966.

* View Article

* PubMed/NCBI

* Google Scholar

78. 78. Coetzee CG. The biology, behaviour, and ecology of Mastomys natalensis in southern Africa. Bull World Health Organ. 1975;52(4–6):637–44. Epub 1975/01/01. pmid:1085219

* View Article

* PubMed/NCBI

* Google Scholar

79. 79. Eby P, Peel AJ, Hoegh A, Madden W, Giles JR, Hudson PJ, et al. Pathogen spillover driven by rapid changes in bat ecology. Nature. 2023;613(7943):340–4. Epub 2022/11/18. pmid:36384167

* View Article

* PubMed/NCBI

* Google Scholar

80. 80. Mehra NK, Kaur G. MHC-based vaccination approaches: progress and perspectives. Expert Rev Mol Med. 2003;5(7):1–17. Epub 2004/02/28. pmid:14987410.

* View Article

* PubMed/NCBI

* Google Scholar

81. 81. Bettencourt P, Muller J, Nicastri A, Cantillon D, Madhavan M, Charles PD, et al. Identification of antigens presented by MHC for vaccines against tuberculosis. NPJ Vaccines. 2020;5(1):2. Epub 2020/01/08. pmid:31908851

* View Article

* PubMed/NCBI

* Google Scholar

Citation: Olayemi A, Schmid DW, Fleischer R, Wilhelm K, Heni AC, Mueller-Klein N, et al. (2024) MHC-I alleles mediate clearance and antibody response to the zoonotic Lassa virus in Mastomys rodent reservoirs. PLoS Negl Trop Dis 18(2): e0011984. https://doi.org/10.1371/journal.pntd.0011984

References

1. Basinski AJ, Fichet-Calvet E, Sjodin AR, Varrelman TJ, Remien CH, Layman NC, et al. Bridging the gap: Using reservoir ecology and human serosurveys to estimate Lassa virus spillover in West Africa. PLoS Comput Biol. 2021;17(3):e1008811. Epub 2021/03/04. pmid:33657095

2. Keane E, Gilles HM. Lassa fever in Panguma Hospital, Sierra Leone, 1973–6. Br Med J. 1977;1(6073):1399–402. Epub 1977/05/28. pmid:861652

3. Price ME, Fisher-Hoch SP, Craven RB, McCormick JB. A prospective study of maternal and fetal outcome in acute Lassa fever infection during pregnancy. BMJ. 1988;297(6648):584–7. Epub 1988/09/03. pmid:3139220

4. Shaffer JG, Grant DS, Schieffelin JS, Boisen ML, Goba A, Hartnett JN, et al. Lassa fever in post-conflict sierra leone. PLoS Negl Trop Dis. 2014;8(3):e2748. Epub 2014/03/22. pmid:24651047 Hemorrhagic Fever Consortium (vhfc.org) is a partnership of academic and industry scientists who are developing diagnostics, therapeutics and vaccines for LF and other severe diseases. Tulane University and various industry partners have filed United States and foreign patent applications on behalf of the consortium for several of these technologies. If commercial products are developed, consortium members may receive royalties or profits. This does not alter our adherence to all PLOS policies on sharing data and materials.

5. Mateer EJ, Huang C, Shehu NY, Paessler S. Lassa fever-induced sensorineural hearing loss: A neglected public health and social burden. PLoS Negl Trop Dis. 2018;12(2):e0006187. Epub 2018/02/23. pmid:29470486

6. Arruda LB, Haider N, Olayemi A, Simons D, Ehichioya D, Yinka-Ogunleye A, et al. The niche of One Health approaches in Lassa fever surveillance and control. Ann Clin Microbiol Antimicrob. 2021;20(1):29. Epub 2021/04/26. pmid:33894784

7. Grange ZL, Goldstein T, Johnson CK, Anthony S, Gilardi K, Daszak P, et al. Ranking the risk of animal-to-human spillover for newly discovered viruses. Proc Natl Acad Sci U S A. 2021;118(15). Epub 2021/04/07. pmid:33822740

8. Mylne AQ, Pigott DM, Longbottom J, Shearer F, Duda KA, Messina JP, et al. Mapping the zoonotic niche of Lassa fever in Africa. Trans R Soc Trop Med Hyg. 2015;109(8):483–92. Epub 2015/06/19. pmid:26085474

9. WHO. Annual Review of Diseases Prioritized under the Research and Development Blueprint Informal Consultation 24–25 January 2017. Geneva, Switzerland.. 2017;http://www.who.int/blueprint/what/research-development/2017-Prioritization-Long-Report.pdf?ua=1.

10. Mehand MS, Al-Shorbaji F, Millett P, Murgue B. The WHO R&D Blueprint: 2018 review of emerging infectious diseases requiring urgent research and development efforts. Antiviral Res. 2018;159:63–7. Epub 2018/09/28. pmid:30261226

11. Garry RF. Lassa fever—the road ahead. Nat Rev Microbiol. 2022:1–10. Epub 2022/09/14. pmid:36097163 countermeasures for Lassa virus (LASV) and other emerging viruses.

12. Monath TP, Newhouse VF, Kemp GE, Setzer HW, Cacciapuoti A. Lassa virus isolation from Mastomys natalensis rodents during an epidemic in Sierra Leone. Science. 1974;185(4147):263–5. Epub 1974/07/19. pmid:4833828.

13. Marien J, Borremans B, Gryseels S, Soropogui B, De Bruyn L, Bongo GN, et al. No measurable adverse effects of Lassa, Morogoro and Gairo arenaviruses on their rodent reservoir host in natural conditions. Parasit Vectors. 2017;10(1):210. Epub 2017/04/30. pmid:28449693

14. Wood R, Bangura U, Marien J, Douno M, Fichet-Calvet E. Detection of Lassa virus in wild rodent feces: Implications for Lassa fever burden within households in the endemic region of Faranah, Guinea. One Health. 2021;13:100317. Epub 2021/09/16. pmid:34522759

15. Marien J, Lo Iacono G, Rieger T, Magassouba N, Gunther S, Fichet-Calvet E. Households as hotspots of Lassa fever? Assessing the spatial distribution of Lassa virus-infected rodents in rural villages of Guinea. Emerg Microbes Infect. 2020;9(1):1055–64. Epub 2020/05/28. pmid:32459576

16. Happold DCD. Mammals of Africa Volume III—Rodents, Hares and Rabbits. Bloomsbury, London. 2013:784 pp.

17. Gouy de Bellocq JB, A; Martynov A. A.; Lavrenchenko L. A. Dhati Welel virus, the missing mammarenavirus of the widespread Mastomys natalensis. Journal of Vertebrate Biology. 2020;69(20018):1–11.

18. Gryseels S, Rieger T, Oestereich L, Cuypers B, Borremans B, Makundi R, et al. Gairo virus, a novel arenavirus of the widespread Mastomys natalensis: Genetically divergent, but ecologically similar to Lassa and Morogoro viruses. Virology. 2015;476:249–56. Epub 2015/01/07. pmid:25559385.

19. Ishii A, Thomas Y, Moonga L, Nakamura I, Ohnuma A, Hang’ombe B, et al. Novel arenavirus, Zambia. Emerg Infect Dis. 2011;17(10):1921–4. Epub 2011/10/18. pmid:22000372

20. Olayemi A, Obadare A, Oyeyiola A, Igbokwe J, Fasogbon A, Igbahenah F, et al. Arenavirus Diversity and Phylogeography of Mastomys natalensis Rodents, Nigeria. Emerg Infect Dis. 2016;22(4):694–7. Epub 2016/03/18. pmid:26982388

21. Borremans B, Leirs H, Gryseels S, Gunther S, Makundi R, de Bellocq JG. Presence of Mopeia virus, an African arenavirus, related to biotope and individual rodent host characteristics: implications for virus transmission. Vector Borne Zoonotic Dis. 2011;11(8):1125–31. Epub 2010/12/15. pmid:21142956.

22. Gunther S, Hoofd G, Charrel R, Roser C, Becker-Ziaja B, Lloyd G, et al. Mopeia virus-related arenavirus in natal multimammate mice, Morogoro, Tanzania. Emerg Infect Dis. 2009;15(12):2008–12. Epub 2009/12/08. pmid:19961688

23. Haikukutu L, Lyaku JR, Lyimo C, Kasanga CJ, Kandusi SE, Rahelinirina S, et al. Plague in Tanzania: first report of sylvatic plague in Morogoro region, persistence in Mbulu focus, and ongoing quiescence in Lushoto and Iringa foci. IJID Reg. 2022;4:105–10. Epub 2022/07/27. pmid:35880003

24. Thomas CMV.; V-Houtte N.; Mhamphi G.; Mariën J.; Sabuni C.; Makundi I.; Nzalawahe J.; Machang’u R.; Leirs H. Molecular detection of Leptospira and Bartonella in Mastomys natalensis and its ectoparasites in Morogoro, Tanzania. Mammalia. 2023;8(5):539–44.

25. Ortiz PA, Garcia HA, Lima L, da Silva FM, Campaner M, Pereira CL, et al. Diagnosis and genetic analysis of the worldwide distributed Rattus-borne Trypanosoma (Herpetosoma) lewisi and its allied species in blood and fleas of rodents. Infect Genet Evol. 2018;63:380–90. Epub 2017/09/09. pmid:28882517.

26. Olayemi A, Cadar D, Magassouba N, Obadare A, Kourouma F, Oyeyiola A, et al. New Hosts of The Lassa Virus. Sci Rep. 2016;6:25280. Epub 2016/05/04. pmid:27140942

27. Olayemi A, Fichet-Calvet E. Systematics, Ecology, and Host Switching: Attributes Affecting Emergence of the Lassa Virus in Rodents across Western Africa. Viruses. 2020;12(3). Epub 2020/03/19. pmid:32183319

28. Apanius V, Penn D, Slev PR, Ruff LR, Potts WK. The nature of selection on the major histocompatibility complex. Crit Rev Immunol. 1997;17(2):179–224. Epub 1997/01/01. pmid:9094452.

29. Janeway CA, Travers P., Walport M., & Sclomchik M. J. Immunobiology: the immune system in health and disease. New York: Garland Science Publishing. 2005.

30. Kaufman J. Generalists and Specialists: A New View of How MHC Class I Molecules Fight Infectious Pathogens. Trends Immunol. 2018;39(5):367–79. Epub 2018/02/06. pmid:29396014

31. Radwan J, Babik W, Kaufman J, Lenz TL, Winternitz J. Advances in the Evolutionary Understanding of MHC Polymorphism. Trends Genet. 2020;36(4):298–311. Epub 2020/02/12. pmid:32044115.

32. Sommer S. The importance of immune gene variability (MHC) in evolutionary ecology and conservation. Front Zool. 2005;2:16. Epub 2005/10/26. pmid:16242022

33. Hawley DM, Fleischer RC. Contrasting epidemic histories reveal pathogen-mediated balancing selection on class II MHC diversity in a wild songbird. PLoS One. 2012;7(1):e30222. Epub 2012/02/01. pmid:22291920

34. Savage AE, Zamudio KR. MHC genotypes associate with resistance to a frog-killing fungus. Proc Natl Acad Sci U S A. 2011;108(40):16705–10. Epub 2011/09/29. pmid:21949385

35. Wegner KM, Kalbe M, Kurtz J, Reusch TB, Milinski M. Parasite selection for immunogenetic optimality. Science. 2003;301(5638):1343. Epub 2003/09/06. pmid:12958352.

36. Froeschke G, Sommer S. MHC class II DRB variability and parasite load in the striped mouse (Rhabdomys pumilio) in the Southern Kalahari. Mol Biol Evol. 2005;22(5):1254–9. Epub 2005/02/11. pmid:15703235.

37. Lenz TL, Wells K, Pfeiffer M, Sommer S. Diverse MHC IIB allele repertoire increases parasite resistance and body condition in the Long-tailed giant rat (Leopoldamys sabanus). BMC Evol Biol. 2009;9:269. Epub 2009/11/26. pmid:19930637

38. Stewart PW, Chapes SK. Role of major histocompatibility complex class II in resistance of mice to naturally acquired infection with Syphacia obvelata. Comp Med. 2003;53(1):70–4. Epub 2003/03/11. pmid:12625509.

39. Gouy de Bellocq J, Leirs H. Polymorphism and signatures of selection in the multimammate rat DQB gene. Immunogenetics. 2010;62(1):59–64. Epub 2009/12/03. pmid:19953242.

40. Hardin A, Nevonen KA, Eckalbar WL, Carbone L, Ahituv N. Comparative Genomic Characterization of the Multimammate Mouse Mastomys coucha. Mol Biol Evol. 2019;36(12):2805–12. Epub 2019/08/20. pmid:31424545

41. Botten J, Alexander J, Pasquetto V, Sidney J, Barrowman P, Ting J, et al. Identification of protective Lassa virus epitopes that are restricted by HLA-A2. J Virol. 2006;80(17):8351–61. Epub 2006/08/17. pmid:16912286

42. Gibb R, Redding DW, Chin KQ, Donnelly CA, Blackburn TM, Newbold T, et al. Zoonotic host diversity increases in human-dominated ecosystems. Nature. 2020;584(7821):398–402. Epub 2020/08/08. pmid:32759999.

43. Johnson CK, Hitchens PL, Pandit PS, Rushmore J, Evans TS, Young CCW, et al. Global shifts in mammalian population trends reveal key predictors of virus spillover risk. Proc Biol Sci. 2020;287(1924):20192736. Epub 2020/04/08. pmid:32259475

44. Fichet-Calvet E, Becker-Ziaja B, Koivogui L, Gunther S. Lassa serology in natural populations of rodents and horizontal transmission. Vector Borne Zoonotic Dis. 2014;14(9):665–74. Epub 2014/09/18. pmid:25229705

45. Borremans B, Vossen R, Becker-Ziaja B, Gryseels S, Hughes N, Van Gestel M, et al. Shedding dynamics of Morogoro virus, an African arenavirus closely related to Lassa virus, in its natural reservoir host Mastomys natalensis. Sci Rep. 2015;5:10445. Epub 2015/05/30. pmid:26022445

46. Olayemi A, Oyeyiola A, Obadare A, Igbokwe J, Adesina AS, Onwe F, et al. Widespread arenavirus occurrence and seroprevalence in small mammals, Nigeria. Parasit Vectors. 2018;11(1):416. Epub 2018/07/15. pmid:30005641

47. Hoffmann C, Wurr S, Pallasch E, Bockholt S, Rieger T, Gunther S, et al. Experimental Morogoro Virus Infection in Its Natural Host, Mastomys natalensis. Viruses. 2021;13(5). Epub 2021/06/03. pmid:34067011

48. Marien J, Borremans B, Gryseels S, Broecke BV, Becker-Ziaja B, Makundi R, et al. Arenavirus Dynamics in Experimentally and Naturally Infected Rodents. Ecohealth. 2017;14(3):463–73. Epub 2017/06/16. pmid:28616660.

49. Westerdahl H, Asghar M, Hasselquist D, Bensch S. Quantitative disease resistance: to better understand parasite-mediated selection on major histocompatibility complex. Proc Biol Sci. 2012;279(1728):577–84. Epub 2011/07/08. pmid:21733902

50. Mills JNY T. L.; Childs J.; Parmenter R. R.; Ksiazek T. G.; Rollin P. E; et al. Guidelines for working with rodents potentially infected with hantavirus. J Mammal. 1995;76:716–22.

51. Adesina AS, Oyeyiola A, Obadare A, Igbokwe J, Abejegah C, Akhilomen P, et al. Circulation of Lassa virus across the endemic Edo-Ondo axis, Nigeria, with cross-species transmission between multimammate mice. Emerg Microbes Infect. 2023;12(1):2219350. Epub 2023/06/08. pmid:37288752

52. Wulff H, Lange JV. Indirect immunofluorescence for the diagnosis of Lassa fever infection. Bull World Health Organ. 1975;52(4–6):429–36. Epub 1975/01/01. pmid:821624

53. Linnenbrink M, Teschke M, Montero I, Vallier M, Tautz D. Meta-populational demes constitute a reservoir for large MHC allele diversity in wild house mice (Mus musculus). Front Zool. 2018;15:15. Epub 2018/05/04. pmid:29721030

54. Rozen S, Skaletsky H. Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000;132:365–86. Epub 1999/11/05. pmid:10547847.

55. Gillingham MAF, Montero BK, Wihelm K, Grudzus K, Sommer S, Santos PSC. A novel workflow to improve genotyping of multigene families in wildlife species: An experimental set-up with a known model system. Mol Ecol Resour. 2021;21(3):982–98. Epub 2020/10/29. pmid:33113273.

56. Fleischer R, Schmid DW, Wasimuddin , Brandel SD, Rasche A, Corman VM, et al. Interaction between MHC diversity and constitution, gut microbiota and Astrovirus infections in a neotropical bat. Mol Ecol. 2022;31(12):3342–59. Epub 2022/05/06. pmid:35510794.

57. Sepil I, Lachish S, Hinks AE, Sheldon BC. Mhc supertypes confer both qualitative and quantitative resistance to avian malaria infections in a wild bird population. Proc Biol Sci. 2013;280(1759):20130134. Epub 2013/03/22. pmid:23516242

58. Sidney J, Grey HM, Kubo RT, Sette A. Practical, biochemical and evolutionary implications of the discovery of HLA class I supermotifs. Immunol Today. 1996;17(6):261–6. Epub 1996/06/01. pmid:8962628.

59. Xu B, Yang Z. PAMLX: a graphical user interface for PAML. Mol Biol Evol. 2013;30(12):2723–4. Epub 2013/10/10. pmid:24105918.

60. Sandberg M, Eriksson L, Jonsson J, Sjostrom M, Wold S. New chemical descriptors relevant for the design of biologically active peptides. A multivariate characterization of 87 amino acids. J Med Chem. 1998;41(14):2481–91. Epub 1998/07/03. pmid:9651153.

61. Jombart T, Devillard S, Balloux F. Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet. 2010;11:94. Epub 2010/10/19. pmid:20950446

62. Anderson MJ. Permutational Multivariate Analysis of Variance. WileyStatsRef: Statistics Reference Online. 2017.

63. Oksanen J, Simpson, G. L., Blanchet, G. F., Roeland Kindt, Pierre, Legendre PRM, R.B. O’Hara, Peter Solymos, M. Henry H. Stevens, Eduard Szoecs HW, Matt Barbour, Michael Bedward, Ben Bolker, Daniel Borcard GC, Michael Chirico, Miquel De Caceres, Sebastien Durand HBAE, Rich FitzJohn, Michael Friendly BF, Geoffrey Hannigan, Mark O. Hill, Leo, et al. vegan: Community Ecology Package. R package version 2.6–2. 2022.

64. Griffith DMM, C. J.;. cooccur: Probabilistic Species Co.Occurrence Analysis in R. Journal of Statistical Software. 2016;131664617.

65. Raberg L, Clough D, Hagstrom A, Scherman K, Andersson M, Drews A, et al. MHC class II genotype-by-pathogen genotype interaction for infection prevalence in a natural rodent-Borrelia system. Evolution. 2022;76(9):2067–75. Epub 2022/08/02. pmid:35909235

66. R-Core-Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2022.

67. Delarbre C, Jaulin C, Kourilsky P, Gachelin G. Evolution of the major histocompatibility complex: a hundred-fold amplification of MHC class I genes in the African pigmy mouse Nannomys setulosus. Immunogenetics. 1992;37(1):29–38. Epub 1992/01/01. pmid:1428009.

68. Axtner J, Sommer S. Gene duplication, allelic diversity, selection processes and adaptive value of MHC class II DRB genes of the bank vole, Clethrionomys glareolus. Immunogenetics. 2007;59(5):417–26. Epub 2007/03/14. pmid:17351770.

69. Miller HC, Lambert DM. Gene duplication and gene conversion in class II MHC genes of New Zealand robins (Petroicidae). Immunogenetics. 2004;56(3):178–91. Epub 2004/05/13. pmid:15138734.

70. Reusch TB, Schaschl H, Wegner KM. Recent duplication and inter-locus gene conversion in major histocompatibility class II genes in a teleost, the three-spined stickleback. Immunogenetics. 2004;56(6):427–37. Epub 2004/08/24. pmid:15322775.

71. Biedrzycka A, O’Connor E, Sebastian A, Migalska M, Radwan J, Zajac T, et al. Extreme MHC class I diversity in the sedge warbler (Acrocephalus schoenobaenus); selection patterns and allelic divergence suggest that different genes have different functions. BMC Evol Biol. 2017;17(1):159. Epub 2017/07/07. pmid:28679358

72. Qurkhuli T, Schwensow N, Brandel SD, Tschapka M, Sommer S. Can extreme MHC class I diversity be a feature of a wide geographic range? The example of Seba’s short-tailed bat (Carollia perspicillata). Immunogenetics. 2019;71(8–9):575–87. Epub 2019/09/15. pmid:31520134

73. Hánová AK A.; Mikula O.; Bryjová A., Šumbera R.; Bryja J. Diversity, distribution, and evolutionary history of the most studied African rodents, multimammate mice of the genus Mastomys: An overview after a quarter of century of using DNA sequencing. Journal of Systematics and Evolutionary Research. 2021;59(8):2500–18.

74. Dunnington EA, Larsen CT, Gross WB, Siegel PB. Antibody responses to combinations of antigens in white Leghorn chickens of different background genomes and major histocompatibility complex genotypes. Poult Sci. 1992;71(11):1801–6. Epub 1992/11/01. pmid:1437966.

75. Zhou H, Lamont SJ. Chicken MHC class I and II gene effects on antibody response kinetics in adult chickens. Immunogenetics. 2003;55(3):133–40. Epub 2003/05/14. pmid:12743657.

76. Migalska M, Przesmycka K, Alsarraf M, Bajer A, Behnke-Borowczyk J, Grzybek M, et al. Long term patterns of association between MHC and helminth burdens in the bank vole support Red Queen dynamics. Mol Ecol. 2022;31(12):3400–15. Epub 2022/05/06. pmid:35510766

77. Oliver MK, Lambin X, Cornulier T, Piertney SB. Spatio-temporal variation in the strength and mode of selection acting on major histocompatibility complex diversity in water vole (Arvicola terrestris) metapopulations. Mol Ecol. 2009;18(1):80–92. Epub 2009/01/15. pmid:19140966.

78. Coetzee CG. The biology, behaviour, and ecology of Mastomys natalensis in southern Africa. Bull World Health Organ. 1975;52(4–6):637–44. Epub 1975/01/01. pmid:1085219

79. Eby P, Peel AJ, Hoegh A, Madden W, Giles JR, Hudson PJ, et al. Pathogen spillover driven by rapid changes in bat ecology. Nature. 2023;613(7943):340–4. Epub 2022/11/18. pmid:36384167

80. Mehra NK, Kaur G. MHC-based vaccination approaches: progress and perspectives. Expert Rev Mol Med. 2003;5(7):1–17. Epub 2004/02/28. pmid:14987410.

81. Bettencourt P, Muller J, Nicastri A, Cantillon D, Madhavan M, Charles PD, et al. Identification of antigens presented by MHC for vaccines against tuberculosis. NPJ Vaccines. 2020;5(1):2. Epub 2020/01/08. pmid:31908851