Selection and recombination drive the evolution of MHC class II DRB diversity in ungulates

Major histocompatibility complex (MHC) antigen-presenting genes are the most variable loci in vertebrate genomes. Host-parasite co-evolution is assumed to maintain the excessive polymorphism in the MHC loci. However, the molecular mechanisms underlying the striking diversity in the MHC remain contentious. The extent to which recombination contributes to the diversity at MHC loci in natural populations is still controversial, and there have been only few comparative studies that make quantitative estimates of recombination rates. In this study, we performed a comparative analysis for 15 different ungulates species to estimate the population recombination rate, and to quantify levels of selection. As expected for all species, we observed signatures of strong positive selection, and identified individual residues experiencing selection that were congruent with those constituting the peptide-binding region of the human DRB gene. However, in addition for each species, we also observed recombination rates that were significantly different from zero on the basis of likelihood-permutation tests, and in other non-quantitative analyses. Patterns of synonymous and non-synonymous sequence diversity were consistent with differing demographic histories between species, but recent simulation studies by other authors suggest inference of selection and recombination is likely to be robust to such deviations from standard models. If high rates of recombination are common in MHC genes of other taxa, re-evaluation of many inference-based phylogenetic analyses of MHC loci, such as estimates of the divergence time of alleles and trans-specific polymorphism, may be required.

What a unicellular parasite can tell us about mitochondrial biogenesis?

Most of what we know about mitochondrial biogenesis stems from work in yeast and mammals, which are quite closely related. To understand the conserved features of mitochondria and the evolutionary forces that shaped it, it is important to study a more diverse group of eukaryotes. The parasitic protozoan Trypanosoma brucei and its relatives are excellent systems to do so, since they appear to have diverged from other eukaryotes very early in evolution. This is reflected in a number of unique and extreme features in their mitochondrial biology, including a single continuous mitochondrion that contains a one unit mitochondrial genome that is physically connected across the two membranes with the basal body of the flagellum. Moreover, many mitochondrial transcripts have to be extensively edited in order to become functional mRNAs and organellar translation requires extensive import of cytosolic tRNAs. In my talk I will focus on the discovery and characterization of the elusive mitochondrial protein import system of the mitochondrial outer membrane of trypanosomes. In addition I will present data on a central outer membrane component of the mitochondrial genome inheritance system of T. brucei and compare it to the better characterized system of yeast. - I hope that I can convince you in my talk, that a better understanding of the mitochondrial biology in T. brucei will provide insights into both fundamentally conserved and fundamentally diverged aspects of mitochondrial biogenesis and thus of the evolutionary hstory of mitochondria in general.

What can a unicellular parasite tell us about mitochondrial biogenesis?

Most of what we know about mitochondrial biogenesis stems from work in yeast and mammals, which are quite closely related. To understand the conserved features of mitochondria and the evolutionary forces that shaped it, it is important to study a more diverse group of eukaryotes. The parasitic protozoan Trypanosoma brucei and its relatives are excellent systems to do so, since they appear to have diverged from other eukaryotes very early in evolution. This is reflected in a number of unique and extreme features in their mitochondrial biology, including a single continuous mitochondrion that contains a one unit mitochondrial genome that is physically connected across the two membranes with the basal body of the flagellum. Moreover, many mitochondrial transcripts have to be extensively edited in order to become functional mRNAs and organellar translation requires extensive import of cytosolic tRNAs. In my talk I will focus on the discovery and characterization of the elusive mitochondrial protein import system of the mitochondrial outer membrane of trypanosomes. In addition I will present data on a central outer membrane component of the mitochondrial genome inheritance system of T. brucei and compare it to the better characterized system of yeast. - I hope that I can convince you in my talk, that a better understanding of the mitochondrial biology in T. brucei will provide insights into both fundamentally conserved and fundamentally diverged aspects of mitochondrial biogenesis and thus of the evolutionary history of mitochondria in general.

Multiple avirulence paralogues in cereal powdery mildew fungi may contribute to parasite fitness and defeat of plant resistance

Powdery mildews, obligate biotrophic fungal parasites on a wide range of important crops, can be controlled by plant resistance (R) genes, but these are rapidly overcome by parasite mutants evading recognition. It is unknown how this rapid evolution occurs without apparent loss of parasite fitness. R proteins recognize avirulence (AVR) molecules from parasites in a gene-for-gene manner and trigger defense responses. We identify AVRa10 and AVRk1 of barley powdery mildew fungus, Blumeria graminis f sp hordei (Bgh), and show that they induce both cell death and naccessibility when transiently expressed in Mla10 and Mlk1 barley (Hordeum vulgare) varieties, respectively. In contrast with other reported fungal AVR genes, AVRa10 and AVRk1 encode proteins that lack secretion signal peptides and enhance infection success on susceptible host plant cells. AVRa10 and AVRk1 belong to a large family with mayor que30 paralogues in the genome of Bgh, and homologous sequences are present in other formae speciales of the fungus infecting other grasses. Our findings imply that the mildew fungus has a repertoire of AVR genes, which may function as effectors and contribute to parasite virulence. Multiple copies of related but distinct AVR effector paralogues might enable populations of Bgh to rapidly overcome host R genes while maintaining virulence.

Population structure and recent evolution of Plasmodium falciparum

Plasmodium falciparum is the agent of malignant malaria, one of mankind's most severe maladies. The parasite exhibits antigenic polymorphisms that have been postulated to be ancient. We have proposed that the extant world populations of P. falciparum have derived from one single parasite, a cenancestor, within the last 5,000–50,000 years. This inference derives from the virtual or complete absence of synonymous nucleotide polymorphisms at genes not involved in immune or drug responses. Seeking to conciliate this claim with extensive antigenic polymorphism, we first note that allele substitutions or polymorphisms can arise very rapidly, even in a single generation, in large populations subject to strong natural selection. Second, new alleles can arise not only by single-nucleotide mutations, but also by duplication/deletion of short simple-repeat DNA sequences, a process several orders of magnitude faster than single-nucleotide mutation. We analyze three antigenic genes known to be extremely polymorphic: Csp, Msp-1, and Msp-2. We identify regions consisting of tandem or proximally repetitive short DNA sequences, including some previously unnoticed. We conclude that the antigenic polymorphisms are consistent with the recent origin of the world populations of P. falciparum inferred from the analysis of nonantigenic genes.

Evolution of resistance to sulfadoxine-pyrimethamine in Plasmodium falciparum

The development of resistance to sulfadoxine-pyrimethamine by Plasmodium parasites is a major problem for the effective treatment of malaria, especially P. falciparum malaria. Although the molecular basis for parasite resistance is known, the factors promoting the development and transmission of these resistant parasites are less clear. This paper reports the results of a quantitative comparison of factors previously hypothesized as important for the development of drug resistance, drug dosage, time of treatment, and drug elimination half-life, with an in-host dynamics model of P. falciparum malaria in a malaria-naive host. The results indicate that the development of drug resistance can be categorized into three stages. The first is the selection of existing parasites with genetic mutations in the dihydrofolate reductase or dihydropteroate synthetase gene. This selection is driven by the long half-life of the sulfadoxine-pyrimethamine combination. The second stage involves the selection of parasites with allelic types of higher resistance within the host during an infection. The timing of treatment relative to initiation of a specific anti-P. falciparum EMP1 immune response is an important factor during this stage, as is the treatment dosage. During the third stage, clinical treatment failure becomes prevalent as the parasites develop sufficient resistance mutations to survive therapeutic doses of the drug combination. Therefore, the model output reaffirms the importance of correct treatment of confirmed malaria cases in slowing the development of parasite resistance to sulfadoxine-pyrimethamine.

Species accumulation curves and their applications in parasite ecology

Species accumulation curves (SACs) chart the increase in recovery of new species as a function of some measure of sampling effort. Studies of parasite diversity can benefit from the application of SACs, both as empirical tools to guide sampling efforts and predict richness, and because their properties are informative about community patterns and the structure of parasite diversity. SACs can be used to infer interactivity in parasite infra-communities, to partition species richness into contributions from different spatial scales and different levels of the host hierarchy (individuals, populations and communities) or to identify modes of community assembly (niche versus dispersal). A historical tendency to treat individual hosts as statistically equivalent replicates (quadrats) seemingly satisfies the sample-based subgroup of SACs but care is required in this because of the inequality of hosts as sampling units. Knowledge of the true distribution of parasite richness over multiple host-derived and spatial scales is far from complete but SACs can improve the understanding of diversity patterns in parasite assemblages.

Applied evolution: An integrated approach to studying life history traits in response to drug selection

The use of chemical control measures to reduce the impact of parasite and pest species has frequently resulted in the development of resistance. Thus, resistance management has become a key concern in human and veterinary medicine, and in agricultural production. Although it is known that factors such as gene flow between susceptible and resistant populations, drug type, application methods, and costs of resistance can affect the rate of resistance evolution, less is known about the impacts of density-dependent eco-evolutionary processes that could be altered by drug-induced mortality. The overall aim of this thesis was to take an experimental evolution approach to assess how life history traits respond to drug selection, using a free-living dioecious worm (Caenorhabditis remanei) as a model. In Chapter 2, I defined the relationship between C. remanei survival and Ivermectin dose over a range of concentrations, in order to control the intensity of selection used in the selection experiment described in Chapter 4. The dose-response data were also used to appraise curve-fitting methods, using Akaike Information Criterion (AIC) model selection to compare a series of nonlinear models. The type of model fitted to the dose response data had a significant effect on the estimates of LD50 and LD99, suggesting that failure to fit an appropriate model could give misleading estimates of resistance status. In addition, simulated data were used to establish that a potential cost of resistance could be predicted by comparing survival at the upper asymptote of dose-response curves for resistant and susceptible populations, even when differences were as low as 4%. This approach to dose-response modeling ensures that the maximum amount of useful information relating to resistance is gathered in one study. In Chapter 3, I asked how simulations could be used to inform important design choices used in selection experiments. Specifically, I focused on the effects of both within- and between-line variation on estimated power, when detecting small, medium and large effect sizes. Using mixed-effect models on simulated data, I demonstrated that commonly used designs with realistic levels of variation could be underpowered for substantial effect sizes. Thus, use of simulation-based power analysis provides an effective way to avoid under or overpowering a study designs incorporating variation due to random effects. In Chapter 4, I 3 investigated how Ivermectin dosage and changes in population density affect the rate of resistance evolution. I exposed replicate lines of C. remanei to two doses of Ivermectin (high and low) to assess relative survival of lines selected in drug-treated environments compared to untreated controls over 10 generations. Additionally, I maintained lines where mortality was imposed randomly to control for differences in density between drug treatments and to distinguish between the evolutionary consequences of drug treatment versus ecological processes affected by changes in density-dependent feedback. Intriguingly, both drug-selected and random-mortality lines showed an increase in survivorship when challenged with Ivermectin; the magnitude of this increase varied with the intensity of selection and life-history stage. The results suggest that interactions between density-dependent processes and life history may mediate evolved changes in susceptibility to control measures, which could result in misleading conclusions about the evolution of heritable resistance following drug treatment. In Chapter 5, I investigated whether the apparent changes in drug susceptibility found in Chapter 4 were related to evolved changes in life-history of C. remanei populations after selection in drug-treated and random-mortality environments. Rapid passage of lines in the drug-free environment had no effect on the measured life-history traits. In the drug-free environment, adult size and fecundity of drug-selected lines increased compared to the controls but drug selection did not affect lifespan. In the treated environment, drug-selected lines showed increased lifespan and fecundity relative to controls. Adult size of randomly culled lines responded in a similar way to drug-selected lines in the drug-free environment, but no change in fecundity or lifespan was observed in either environment. The results suggest that life histories of nematodes can respond to selection as a result of the application of control measures. Failure to take these responses into account when applying control measures could result in adverse outcomes, such as larger and more fecund parasites, as well as over-estimation of the development of genetically controlled resistance. In conclusion, my thesis shows that there may be a complex relationship between drug selection, density-dependent regulatory processes and life history of populations challenged with control measures. This relationship could have implications for how resistance is monitored and managed if life histories of parasitic species show such eco-evolutionary responses to drug application.

Prognostic Value Of Dna And Mrna E6/e7 Of Human Papillomavirus In The Evolution Of Cervical Intraepithelial Neoplasia Grade 2.

This study aimed at evaluating whether human papillomavirus (HPV) groups and E6/E7 mRNA of HPV 16, 18, 31, 33, and 45 are prognostic of cervical intraepithelial neoplasia (CIN) 2 outcome in women with a cervical smear showing a low-grade squamous intraepithelial lesion (LSIL). This cohort study included women with biopsy-confirmed CIN 2 who were followed up for 12 months, with cervical smear and colposcopy performed every three months. Women with a negative or low-risk HPV status showed 100% CIN 2 regression. The CIN 2 regression rates at the 12-month follow-up were 69.4% for women with alpha-9 HPV versus 91.7% for other HPV species or HPV-negative status (P < 0.05). For women with HPV 16, the CIN 2 regression rate at the 12-month follow-up was 61.4% versus 89.5% for other HPV types or HPV-negative status (P < 0.05). The CIN 2 regression rate was 68.3% for women who tested positive for HPV E6/E7 mRNA versus 82.0% for the negative results, but this difference was not statistically significant. The expectant management for women with biopsy-confirmed CIN 2 and previous cytological tests showing LSIL exhibited a very high rate of spontaneous regression. HPV 16 is associated with a higher CIN 2 progression rate than other HPV infections. HPV E6/E7 mRNA is not a prognostic marker of the CIN 2 clinical outcome, although this analysis cannot be considered conclusive. Given the small sample size, this study could be considered a pilot for future larger studies on the role of predictive markers of CIN 2 evolution.

Effects Of Therapeutic Approach On The Neonatal Evolution Of Very Low Birth Weight Infants With Patent Ductus Arteriosus.

To analyze the effects of treatment approach on the outcomes of newborns (birth weight [BW] < 1,000 g) with patent ductus arteriosus (PDA), from the Brazilian Neonatal Research Network (BNRN) on: death, bronchopulmonary dysplasia (BPD), severe intraventricular hemorrhage (IVH III/IV), retinopathy of prematurity requiring surgical (ROPsur), necrotizing enterocolitis requiring surgery (NECsur), and death/BPD. This was a multicentric, cohort study, retrospective data collection, including newborns (BW < 1000 g) with gestational age (GA) < 33 weeks and echocardiographic diagnosis of PDA, from 16 neonatal units of the BNRN from January 1, 2010 to Dec 31, 2011. Newborns who died or were transferred until the third day of life, and those with presence of congenital malformation or infection were excluded. Groups: G1 - conservative approach (without treatment), G2 - pharmacologic (indomethacin or ibuprofen), G3 - surgical ligation (independent of previous treatment). Factors analyzed: antenatal corticosteroid, cesarean section, BW, GA, 5 min. Apgar score < 4, male gender, Score for Neonatal Acute Physiology Perinatal Extension (SNAPPE II), respiratory distress syndrome (RDS), late sepsis (LS), mechanical ventilation (MV), surfactant (< 2 h of life), and time of MV. death, O2 dependence at 36 weeks (BPD36wks), IVH III/IV, ROPsur, NECsur, and death/BPD36wks. Student's t-test, chi-squared test, or Fisher's exact test; Odds ratio (95% CI); logistic binary regression and backward stepwise multiple regression. Software: MedCalc (Medical Calculator) software, version 12.1.4.0. p-values < 0.05 were considered statistically significant. 1,097 newborns were selected and 494 newborns were included: G1 - 187 (37.8%), G2 - 205 (41.5%), and G3 - 102 (20.6%). The highest mortality was observed in G1 (51.3%) and the lowest in G3 (14.7%). The highest frequencies of BPD36wks (70.6%) and ROPsur were observed in G3 (23.5%). The lowest occurrence of death/BPD36wks occurred in G2 (58.0%). Pharmacological (OR 0.29; 95% CI: 0.14-0.62) and conservative (OR 0.34; 95% CI: 0.14-0.79) treatments were protective for the outcome death/BPD36wks. The conservative approach of PDA was associated to high mortality, the surgical approach to the occurrence of BPD36wks and ROPsur, and the pharmacological treatment was protective for the outcome death/BPD36wks.

Evolution Of The Deaths Registry System In Brazil: Associations With Changes In The Mortality Profile, Under-registration Of Death Counts, And Ill-defined Causes Of Death.

This paper examines the spatial pattern of ill-defined causes of death across Brazilian regions, and its relationship with the evolution of completeness of the deaths registry and changes in the mortality age profile. We make use of the Brazilian Health Informatics Department mortality database and population censuses from 1980 to 2010. We applied demographic methods to evaluate the quality of mortality data for 137 small areas and correct for under-registration of death counts when necessary. The second part of the analysis uses linear regression models to investigate the relationship between, on the one hand, changes in death counts coverage and age profile of mortality, and on the other, changes in the reporting of ill-defined causes of death. The completeness of death counts coverage increases from about 80% in 1980-1991 to over 95% in 2000-2010 at the same time the percentage of ill-defined causes of deaths reduced about 53% in the country. The analysis suggests that the government's efforts to improve data quality are proving successful, and they will allow for a better understanding of the dynamics of health and the mortality transition.

Digestive Peptidase Evolution In Holometabolous Insects Led To A Divergent Group Of Enzymes In Lepidoptera.

Trypsins and chymotrypsins are well-studied serine peptidases that cleave peptide bonds at the carboxyl side of basic and hydrophobic l-amino acids, respectively. These enzymes are largely responsible for the digestion of proteins. Three primary processes regulate the activity of these peptidases: secretion, precursor (zymogen) activation and substrate-binding site recognition. Here, we present a detailed phylogenetic analysis of trypsins and chymotrypsins in three orders of holometabolous insects and reveal divergent characteristics of Lepidoptera enzymes in comparison with those of Coleoptera and Diptera. In particular, trypsin subsite S1 was more hydrophilic in Lepidoptera than in Coleoptera and Diptera, whereas subsites S2-S4 were more hydrophobic, suggesting different substrate preferences. Furthermore, Lepidoptera displayed a lineage-specific trypsin group belonging only to the Noctuidae family. Evidence for facilitated trypsin auto-activation events were also observed in all the insect orders studied, with the characteristic zymogen activation motif complementary to the trypsin active site. In contrast, insect chymotrypsins did not seem to have a peculiar evolutionary history with respect to their mammal counterparts. Overall, our findings suggest that the need for fast digestion allowed holometabolous insects to evolve divergent groups of peptidases with high auto-activation rates, and highlight that the evolution of trypsins led to a most diverse group of enzymes in Lepidoptera.

Time evolution of the South Atlantic Magnetic Anomaly

The South Atlantic Magnetic Anomaly (SAMA) is one of the most outstanding anomalies of the geomagnetic field. The SAMA secular variation was obtained and compared to the evolution of other anomalies using spherical harmonic field models for the 1590-2005 period. An analysis of data from four South American observatories shows how this large scale anomaly affected their measurements. Since SAMA is a low total field anomaly, the field was separated into its nondipolar, quadrupolar and octupolar parts. The time evolution of the non-dipole/total, quadrupolar/total and octupolar/total field ratios yielded increasingly high values for the South Atlantic since 1750. The SAMA evolution is compared to the evolution of other large scale surface geomagnetic features like the North and the South Pole and the Siberia High, and this comparison shows the intensity equilibrium between these anomalies in both hemispheres. The analysis of non-dipole fields in historical period suggests that SAMA is governed by (i) quadrupolar field for drift, and (ii) quadrupolar and octupolar fields for intensity and area of influence. Furthermore, our study reinforces the possibility that SAMA may be related to reverse fluxes in the outer core under the South Atlantic region.

Leishmania spp. parasite isolation through inoculation of patient biopsy macerates in interferon gamma knockout mice

Isolation of Leishmania parasite and species identification are important for confirmation and to help define the epidemiology of the leishmaniasis. Mice are often used to isolate pathogens, but the most common mouse strains are resistant to infection with parasites from the Leishmania (Viannia) subgenus. In this study we tested the inoculation of interferon gamma knockout (IFNγ KO) mice with biopsy macerates from Leishmania-infected patients to increase the possibility of isolating parasites. Biopsies from twenty five patients with clinical signs of leishmaniasis were taken and tested for the presence of parasites. Immunohistochemical assay (IHC) and conventional histopathology detected the parasite in 88% and 83% of the patients, respectively. Leishmania sp. were isolated in biopsy macerates from 52% of the patients by culture in Grace's insect medium, but 13% of isolates were lost due to contamination. Inoculation of macerates in IFNγ KO mice provides isolation of parasites in 31.8% of the biopsies. Most isolates belong to L. (Viannia) subgenus, as confirmed by PCR, except one that belongs to L. (Leishmania) subgenus. Our preliminary results support the use of IFNγ KO mice to improve the possibility to isolate New World Leishmania species.