The heptad repeat region is a major selection target in MERS-CoV and related coronaviruses

Middle East respiratory syndrome coronavirus (MERS-CoV) originated in bats and spread to humans via zoonotic transmission from camels. We analyzed the evolution of the spike (S) gene in betacoronaviruses (betaCoVs) isolated from different mammals, in bat coronavirus populations, as well as in MERS-CoV strains from the current outbreak. Results indicated several positively selected sites located in the region comprising the two heptad repeats (HR1 and HR2) and their linker. Two sites (R652 and V1060) were positively selected in the betaCoVs phylogeny and correspond to mutations associated with expanded host range in other coronaviruses. During the most recent evolution of MERS-CoV, adaptive mutations in the HR1 (Q/R/H1020) arose in camels or in a previous host and spread to humans. We determined that different residues at position 1020 establish distinct inter-and intra-helical interactions and affect the stability of the six-helix bundle formed by the HRs. A similar effect on stability was observed for a nearby mutation (T1015N) that increases MERS-CoV infection efficiency in vitro. Data herein indicate that the heptad repeat region was a major target of adaptive evolution in MERS-CoV-related viruses; these results are relevant for the design of fusion inhibitor peptides with antiviral function

Recipient Criteria Predictive of Graft Failure in Kidney Transplantation

Several classifications systems have been developed to predict outcomes of kidney transplantation based on donor variables. This study aims to identify kidney transplant recipient variables that would predict graft outcome irrespective of donor characteristics. All U.S. kidney transplant recipients between October 25,1999 and January 1, 2007 were reviewed. Cox proportional hazards regression was used to model time until graft failure. Death-censored and nondeath-censored graft survival models were generated for recipients of live and deceased donor organs. Recipient age, gender, body mass index (BMI), presence of cardiac risk factors, peripheral vascular disease, pulmonary disease, diabetes, cerebrovascular disease, history of malignancy, hepatitis B core antibody, hepatitis C infection, dialysis status, panel-reactive antibodies (PRA), geographic region, educational level, and prior kidney transplant were evaluated in all kidney transplant recipients. Among the 88,284 adult transplant recipients the following groups had increased risk of graft failure: younger and older recipients, increasing PRA (hazard ratio [HR],1.03-1.06], increasing BMI (HR, 1.04-1.62), previous kidney transplant (HR, 1.17-1.26), dialysis at the time of transplantation (HR, 1.39-1.51), hepatitis C infection (HR, 1.41-1.63), and educational level (HR, 1.05-1.42). Predictive criteria based on recipient characteristics could guide organ allocation, risk stratification, and patient expectations in planning kidney transplantation

Signatures of Environmental Genetic Adaptation Pinpoint Pathogens as the Main Selective Pressure through Human Evolution

Previous genome-wide scans of positive natural selection in humans have identified a number of non-neutrally evolving genes that play important roles in skin pigmentation, metabolism, or immune function. Recent studies have also shown that a genome-wide pattern of local adaptation can be detected by identifying correlations between patterns of allele frequencies and environmental variables. Despite these observations, the degree to which natural selection is primarily driven by adaptation to local environments, and the role of pathogens or other ecological factors as selective agents, is still under debate. To address this issue, we correlated the spatial allele frequency distribution of a large sample of SNPs from 55 distinct human populations to a set of environmental factors that describe local geographical features such as climate, diet regimes, and pathogen loads. In concordance with previous studies, we detected a significant enrichment of genic SNPs, and particularly non-synonymous SNPs associated with local adaptation. Furthermore, we show that the diversity of the local pathogenic environment is the predominant driver of local adaptation, and that climate, at least as measured here, only plays a relatively minor role. While background demography by far makes the strongest contribution in explaining the genetic variance among populations, we detected about 100 genes which show an unexpectedly strong correlation between allele frequencies and pathogenic environment, after correcting for demography. Conversely, for diet regimes and climatic conditions, no genes show a similar correlation between the environmental factor and allele frequencies. This result is validated using low-coverage sequencing data for multiple populations. Among the loci targeted by pathogen-driven selection, we found an enrichment of genes associated to autoimmune diseases, such as celiac disease, type 1 diabetes, and multiples sclerosis, which lends credence to the hypothesis that some susceptibility alleles for autoimmune diseases may be maintained in human population due to past selective processes

Performance Assessment in Fingerprinting and Multi Component Quantitative NMR Analyses

An interlaboratory comparison (ILC) was organized with the aim to set up quality control indicators suitable for multicomponent quantitative analysis by nuclear magnetic resonance (NMR) spectroscopy. A total of 36 NMR data sets (corresponding to 1260 NMR spectra) were produced by 30 participants using 34 NMR spectrometers. The calibration line method was chosen for the quantification of a five-component model mixture. Results show that quantitative NMR is a robust quantification tool and that 26 out of 36 data sets resulted in statistically equivalent calibration lines for all considered NMR signals. The performance of each laboratory was assessed by means of a new performance index (named Qp-score) which is related to the difference between the experimental and the consensus values of the slope of the calibration lines. Laboratories endowed with a Qp-score falling within the suitable acceptability range are qualified to produce NMR spectra that can be considered statistically equivalent in terms of relative intensities of the signals. In addition, the specific response of nuclei to the experimental excitation/relaxation conditions was addressed by means of the parameter named NR. NR is related to the difference between the theoretical and the consensus slopes of the calibration lines and is specific for each signal produced by a well-defined set of acquisition parameters

A Functional Variant in ERAP1 Predisposes to Multiple Sclerosis

The ERAP1 gene encodes an aminopeptidase involved in antigen processing. A functional polymorphism in the gene (rs30187, Arg528Lys) associates with susceptibility to ankylosying spondylitis (AS), whereas a SNP in the interacting ERAP2 gene increases susceptibility to another inflammatory autoimmune disorder, Crohn's disease (CD). We analysed rs30187 in 572 Italian patients with CD and in 517 subjects suffering from multiple sclerosis (MS); for each cohort, an independent sex- and age-matched control group was genotyped. The frequency of the 528Arg allele was significantly higher in both disease cohorts compared to the respective control population (for CD, OR = 1.20 95%CI: 1.01–1.43, p = 0.036; for RRMS, OR = 1.26; 95%CI: 1.04–1.51, p = 0.01). Meta-analysis with the Wellcome Trust Cases Control Consortium GWAS data confirmed the association with MS (pmeta = 0.005), but not with CD. In AS, the rs30187 variant has a predisposing effect only in an HLA-B27 allelic background. It remains to be evaluated whether interaction between ERAP1 and distinct HLA class I alleles also affects the predisposition to MS, and explains the failure to provide definitive evidence for a role of rs30187 in CD. Results herein support the emerging concept that a subset of master-regulatory genes underlay the pathogenesis of autoimmunity

A regulatory polymorphism in HAVCR2 modulates susceptibility to HIV-1 infection

The HAVCR2 gene encodes TIM-3, an immunoglobulin superfamily member expressed by exhausted CD8+ T cells during chronic viral infection. We investigated whether genetic variation at HAVCR2 modulates the susceptibility to HIV-1 acquisition; specifically we focused on a 3\u2032 UTR variant (rs4704846, A/G) that represents a natural selection target. We genotyped rs4704846 in three independent cohorts of HIV-1 exposed seronegative (HESN) individuals with different geographic origin (Italy and Spain) and distinct route of exposure to HIV-1 (sexual and injection drug use). Matched HIV-1 positive subjects and healthy controls were also analyzed. In all case-control cohorts the minor G allele at rs4704846 was more common in HIV-1 infected individuals than in HESN, with healthy controls showing intermediate frequency. Results from the three association analyses were combined through a random effect meta-analysis, which revealed no heterogeneity among samples (Cochrane's Q, p value = 0.89, I2 = 0) and yielded a p value of 6.8 710 124. The minor G allele at rs4704846 was found to increase HAVCR2 expression after in vitro HIV-1 infection. Thus, a positively selected polymorphism in the 3\u2032 UTR, which modulates HAVCR2 expression, is associated with the susceptibility to HIV-1 infection. These data warrant further investigation into the role of TIM-3 in the prevention and treatment of HIV-1/AIDS

An evolutionary analysis of antigen processing and presentation across different timescales reveals pervasive selection

The antigenic repertoire presented by MHC molecules is generated by the antigen processing and presentation (APP) pathway. We analyzed the evolutionary history of 45 genes involved in APP at the inter- and intra-species level. Results showed that 11 genes evolved adaptively in mammals. Several positively selected sites involve positions of fundamental importance to the protein function (e.g. the TAP1 peptide-binding domains, the sugar binding interface of langerin, and the CD1D trafficking signal region). In CYBB, all selected sites cluster in two loops protruding into the endosomal lumen; analysis of missense mutations responsible for chronic granulomatous disease (CGD) showed the action of different selective forces on the very same gene region, as most CGD substitutions involve aminoacid positions that are conserved in all mammals. As for ERAP2, different computational methods indicated that positive selection has driven the recurrent appearance of protein-destabilizing variants during mammalian evolution. Application of a population-genetics phylogenetics approach showed that purifying selection represented a major force acting on some APP components (e.g. immunoproteasome subunits and chaperones) and allowed identification of positive selection events in the human lineage. We also investigated the evolutionary history of APP genes in human populations by developing a new approach that uses several different tests to identify the selection target, and that integrates low-coverage whole-genome sequencing data with Sanger sequencing. This analysis revealed that 9 APP genes underwent local adaptation in human populations. Most positive selection targets are located within noncoding regions with regulatory function in myeloid cells or act as expression quantitative trait loci. Conversely, balancing selection targeted nonsynonymous variants in TAP1 and CD207 (langerin). Finally, we suggest that selected variants in PSMB10 and CD207 contribute to human phenotypes. Thus, we used evolutionary information to generate experimentally-testable hypotheses and to provide a list of sites to prioritize in follow-up analyses

Genome-Wide Identification of Susceptibility Alleles for Viral Infections through a Population Genetics Approach

Viruses have exerted a constant and potent selective pressure on human genes throughout evolution. We utilized the marks left by selection on allele frequency to identify viral infection-associated allelic variants. Virus diversity (the number of different viruses in a geographic region) was used to measure virus-driven selective pressure. Results showed an excess of variants correlated with virus diversity in genes involved in immune response and in the biosynthesis of glycan structures functioning as viral receptors; a significantly higher than expected number of variants was also seen in genes encoding proteins that directly interact with viral components. Genome-wide analyses identified 441 variants significantly associated with virus-diversity; these are more frequently located within gene regions than expected, and they map to 139 human genes. Analysis of functional relationships among genes subjected to virus-driven selective pressure identified a complex network enriched in viral products-interacting proteins. The novel approach to the study of infectious disease epidemiology presented herein may represent an alternative to classic genome-wide association studies and provides a large set of candidate susceptibility variants for viral infections

Complex Interplay of Evolutionary Forces in the ladybird Homeobox Genes of Drosophila melanogaster

Tandemly arranged paralogous genes lbe and lbl are members of the Drosophila NK homeobox family. We analyzed population samples of Drosophila melanogaster from Africa, Europe, North and South America, and single strains of D. sechellia, D. simulans, and D. yakuba within two linked regions encompassing partial sequences of lbe and lbl. The evolution of lbe and lbl is highly constrained due to their important regulatory functions. Despite this, a variety of forces have shaped the patterns of variation in lb genes: recombination, intragenic gene conversion and natural selection strongly influence background variation created by linkage disequilibrium and dimorphic haplotype structure. The two genes exhibited similar levels of nucleotide diversity and positive selection was detected in the noncoding regions of both genes. However, synonymous variability was significantly higher for lbe: no nonsynonymous changes were observed in this gene. We argue that balancing selection impacts some synonymous sites of the lbe gene. Stability of mRNA secondary structure was significantly different between the lbe (but not lbl) haplotype groups and may represent a driving force of balancing selection in epistatically interacting synonymous sites. Balancing selection on synonymous sites may be the first, or one of a few such observations, in Drosophila. In contrast, recurrent positive selection on lbl at the protein level influenced evolution at three codon sites. Transcription factor binding-site profiles were different for lbe and lbl, suggesting that their developmental functions are not redundant. Combined with our previous results on nucleotide variation in esterase and other homeobox genes, these results suggest that interplay of balancing and directional selection may be a general feature of molecular evolution in Drosophila and other eukaryote genomes