ALBERT

Description: The human ZC3HAV1 gene encodes an antiviral protein. The longest splicing isoform of ZC3HAV1 contains a C-terminal PARP-like domain, which has evolved under positive selection in primates. We analyzed the evolutionary history of this same domain in humans and in Pan troglodytes . We identified two variants that segregate in both humans and chimpanzees; one of them (rs3735007) does not occur at a hypermutable site and accounts for a nonsynonymous substitution (Thr851Ile). The probability that the two trans -specific polymorphisms have occurred independently in the two lineages was estimated to be low ( P = 0.0054), suggesting that at least one of them has arisen before speciation and has been maintained by selection. Population genetic analyses in humans indicated that the region surrounding the shared variants displays strong evidences of long-standing balancing selection. Selection signatures were also observed in a chimpanzee population sample. Inspection of 1000 Genomes data confirmed these findings but indicated that search for selection signatures using low-coverage whole-genome data may need masking of repetitive sequences. A case–control study of more than 1,000 individuals from mainland Italy indicated that the Thr851Ile SNP is significantly associated with susceptibility to multiple sclerosis (MS) (odds ratio [OR] = 1.47, 95% confidence intervals [CI]: 1.08–1.99, P = 0.011). This finding was confirmed in a larger sample of 4,416 Sardinians cases/controls (OR = 1.18, 95% CI: 1.037–1.344, P = 0.011), but not in a population from Belgium. We provide one of the first instances of human/chimpanzee trans -specific coding variant located outside the major histocompatibility complex region. The selective pressure is likely to be virus driven; in modern populations, this variant associates with susceptibility to MS, possibly via the interaction with environmental factors.

Description: SUMMARY The study of glacial isostatic adjustment (GIA) is gaining an increasingly important role within the geophysical community. Understanding the response of the Earth to loading is crucial in various contexts, ranging from the interpretation of modern satellite geodetic measurements (e.g. GRACE and GOCE) to the projections of future sea level trends in response to climate change. Modern modelling approaches to GIA are based on various techniques that range from purely analytical formulations to fully numerical methods. Despite various teams independently investigating GIA, we do not have a suitably large set of agreed numerical results through which the methods may be validated; a community benchmark data set would clearly be valuable. Following the example of the mantle convection community, here we present, for the first time, the results of a benchmark study of codes designed to model GIA. This has taken place within a collaboration facilitated through European Cooperation in Science and Technology (COST) Action ES0701. The approaches benchmarked are based on significantly different codes and different techniques. The test computations are based on models with spherical symmetry and Maxwell rheology and include inputs from different methods and solution techniques: viscoelastic normal modes, spectral-finite elements and finite elements. The tests involve the loading and tidal Love numbers and their relaxation spectra, the deformation and gravity variations driven by surface loads characterized by simple geometry and time history and the rotational fluctuations in response to glacial unloading. In spite of the significant differences in the numerical methods employed, the test computations show a satisfactory agreement between the results provided by the participants.

Description: During megathrust earthquakes, great ruptures are accompanied by large scale mass redistribution inside the solid Earth and by ocean mass redistribution due to bathymetry changes. These large scale mass displacements can be detected using the monthly gravity maps of the GRACE satellite mission. In recent years it has become increasingly common to use the long wavelength changes in the Earth's gravity field observed by GRACE to infer seismic source properties for large megathrust earthquakes. An important advantage of space gravimetry is that it is independent from the availability of land for its measurements. This is relevant for observation of megathrust earthquakes, which occur mostly offshore, such as the $M_{\text{w}} \sim 9$ 2004 Sumatra–Andaman, 2010 Maule (Chile) and 2011 Tohoku-Oki (Japan) events. In Broerse et al. , we examined the effect of the presence of an ocean above the rupture on long wavelength gravity changes and showed it to be of the first order. Here we revisit the implementation of an ocean layer through the sea level equation and compare the results with approximated methods that have been used in the literature. One of the simplifications usually lies in the assumption of a globally uniform ocean layer. We show that especially in the case of the 2010 Maule earthquake, due to the closeness of the South American continent, the uniform ocean assumption is not valid and causes errors up to 57 per cent for modelled peak geoid height changes (expressed at a spherical harmonic truncation degree of 40). In addition, we show that when a large amount of slip occurs close to the trench, horizontal motions of the ocean floor play a mayor role in the ocean contribution to gravity changes. Using a slip model of the 2011 Tohoku-Oki earthquake that places the majority of slip close to the surface, the peak value in geoid height change increases by 50 per cent due to horizontal ocean floor motion. Furthermore, we test the influence of the maximum spherical harmonic degree at which the sea level equation is performed for sea level changes occurring along coastlines, which shows to be important for relative sea level changes occurring along the shore. Finally, we demonstrate that ocean floor loading, self-gravitation of water and conservation of water mass are of second order importance for coseismic gravity changes. When GRACE observations are used to determine earthquake parameters such as seismic moment or source depth, the uniform ocean layer method introduces large biases, depending on the location of the rupture with respect to the continent. The same holds for interpreting shallow slip when horizontal motions are not properly accounted for in the ocean contribution. In both cases the depth at which slip occurs will be underestimated.

Description: AIM2-like receptors (ALRs) are a family of nucleic acid sensors essential for innate immune responses against viruses and bacteria. We performed an evolutionary analysis of ALR genes ( MNDA, PYHIN1, IFI16 , and AIM2 ) by analyzing inter- and intraspecies diversity. Maximum-likelihood analyses indicated that IFI16 and AIM2 evolved adaptively in primates, with branch-specific selection at the catarrhini lineage for IFI16 . Application of a population genetics–phylogenetics approach also allowed identification of positive selection events in the human lineage. Positive selection in primates targeted sites located at the DNA-binding interface in both IFI16 and AIM2. In IFI16, several sites positively selected in primates and in the human lineage were located in the PYD domain, which is involved in protein–protein interaction and is bound by a human cytomegalovirus immune evasion protein. Finally, positive selection was found to target nuclear localization signals in IFI16 and the spacer region separating the two HIN domains. Population genetic analysis in humans revealed that an IFI16 genic region has been a target of long-standing balancing selection, possibly acting on two nonsynonymous polymorphisms located in the spacer region. Data herein indicate that ALRs have been repeatedly targeted by natural selection. The balancing selection region in IFI16 carries a variant with opposite risk effect for distinct autoimmune diseases, suggesting antagonistic pleiotropy. We propose that the underlying scenario is the result of an ancestral and still ongoing host–pathogen arms race and that the maintenance of susceptibility alleles for autoimmune diseases at IFI16 represents an evolutionary trade-off.

Description: The protein product of the myxovirus resistance 2 ( MX2 ) gene restricts HIV-1 and simian retroviruses. We demonstrate that MX2 evolved adaptively in mammals with distinct sites representing selection targets in distinct branches; selection mainly involved residues in loop 4, previously shown to carry antiviral determinants. Modeling data indicated that positively selected sites form a continuous surface on loop 4, which folds into two antiparallel α-helices protruding from the stalk domain. A population genetics–phylogenetics approach indicated that the coding region of MX2 mainly evolved under negative selection in the human lineage. Nonetheless, population genetic analyses demonstrated that natural selection operated on MX2 during the recent history of human populations: distinct selective events drove the frequency increase of two haplotypes in the populations of Asian and European ancestry. The Asian haplotype carries a susceptibility allele for melanoma; the European haplotype is tagged by rs2074560, an intronic variant. Analyses performed on three independent European cohorts of HIV-1-exposed seronegative individuals with different geographic origin and distinct exposure route showed that the ancestral (G) allele of rs2074560 protects from HIV-1 infection with a recessive effect (combined P = 1.55 x 10 –4 ). The same allele is associated with lower in vitro HIV-1 replication and increases MX2 expression levels in response to IFN-α. Data herein exploit evolutionary information to identify a novel host determinant of HIV-1 infection susceptibility.

Description: Little is known regarding the post-transcriptional networks that control gene expression in eukaryotes. Additionally, we still need to understand how these networks evolve, and the relative role played in them by their sequence-dependent regulatory factors, non-coding RNAs (ncRNAs) and RNA-binding proteins (RBPs). Here, we used an approach that relied on both phylogenetic sequence sharing and conservation in the whole mapped 3'-untranslated regions (3'-UTRs) of vertebrate species to gain knowledge on core post-transcriptional networks. The identified human hyper conserved elements (HCEs) were predicted to be preferred binding sites for RBPs and not for ncRNAs, namely microRNAs and long ncRNAs. We found that the HCE map identified a well-known network that post-transcriptionally regulates histone mRNAs. We were then able to discover and experimentally confirm a translational network composed of RNA Recognition Motif (RRM)-type RBP mRNAs that are positively controlled by HuR, another RRM-type RBP. HuR shows a preference for these RBP mRNAs bound in stem–loop motifs, confirming its role as a ‘regulator of regulators’. Analysis of the transcriptome-wide HCE distribution revealed a profile of prevalently small clusters separated by unconserved intercluster RNA stretches, which predicts the formation of discrete small ribonucleoprotein complexes in the 3'-UTRs.

Description: We present relative astrometric measurements of visual binaries, made in 2011 with the Pupil Interferometry Speckle camera and COronagraph (PISCO) at the 102-cm Zeiss telescope of the Brera Astronomical Observatory in Merate (Italy). Our observing list contains orbital couples as well as binaries whose motion is still uncertain. We obtained new measurements of 469 objects, with angular separations in the range 0.14–8.1 arcsec, and an average accuracy of 0.02 arcsec. The mean error on the position angles is 0 $_{.}^{\circ}$ 7. Most of the position angles were determined without the usual 180° ambiguity with the application of triple-correlation techniques and/or by inspection of the long integration files. Thanks to a new low-magnification option included in PISCO, we have been able to observe fainter stars than previously. The limiting magnitude of our instrumentation on the Zeiss telescope is now close to m V  = 10–12, which permits the observation of some red dwarfs. Finally, we present new revised orbits for ADS 8739, 9182 Aa,Ab, 9626 Ba,Bb, 12880 and 14412, partly derived from those observations. The corresponding estimated values for the masses of those systems are compatible with the spectral types.

Description: Activation of the contact system leads to the cleavage of kininogen by plasma kallikrein resulting in kinin release and in the initiation of the intrinsic pathway of coagulation. Proteolysis of kininogen also generates antimicrobial peptides (AMPs) and can be induced by diverse pathogens. Thus, the contact system is regarded as a branch of innate immunity. We performed an evolutionary analysis of contact system genes by analyzing both inter- and intraspecies diversity. Results indicated that mammalian kininogen genes evolved adaptively. Positively selected sites are located in all protein domains with the exclusion of the bradykinin region and also involve AMP sequences (including the highly effective NAT26 peptide); positively selected sites also occur at alternative cleavage sites for neutrophil-released kinins. Population genetic analysis in humans indicated that a region of the kininogen gene ( KNG1 ) has been a target of long-standing multiallelic balancing selection and that the coalescence time of the haplotype phylogeny dates back to the split between the humans and chimpanzees. No selection signature was detected in the Pan troglodytes KNG1 gene or in human genes encoding other components of the contact system. The selection targets in human KNG1 might be accounted for by variants with transcriptional regulatory activity. Results herein indicate a continuum in selective pressure acting on different timescales and targeting KNG1 . This is in line with evidences suggesting a central role for kininogen in modulating of immune response and with its being a target of an extremely diverse array of pathogen species.

Description: Toxic metals are part of our environment, and undue exposure to them leads to a variety of pathologies. In response, most organisms adapt their metabolism and have evolved systems to limit this toxicity and to acquire tolerance. Ribosome biosynthesis being central for protein synthesis, we analyzed in yeast the effects of a moderate concentration of cadmium (Cd 2+ ) on Pol I transcription that represents 〉60% of the transcriptional activity of the cells. We show that Cd 2+ rapidly and drastically shuts down the expression of the 35S rRNA. Repression does not result from a poisoning of any of the components of the class I transcriptional machinery by Cd 2+ , but rather involves a protein phosphatase 2A (PP2A)-dependent cellular signaling pathway that targets the formation/dissociation of the Pol I–Rrn3 complex. We also show that Pol I transcription is repressed by other toxic metals, such as Ag + and Hg 2+ , which likewise perturb the Pol I–Rrn3 complex, but through PP2A-independent mechanisms. Taken together, our results point to a central role for the Pol I–Rrn3 complex as molecular switch for regulating Pol I transcription in response to toxic metals.

Description: Previous studies indicated that a few risk variants for autoimmune diseases are subject to pathogen-driven selection. Nonetheless, the proportion of risk loci that has been targeted by pathogens and the type of infectious agent(s) that exerted the strongest pressure remain to be evaluated. We assessed whether different pathogens exerted a pressure on known Crohn's disease (CD) risk variants and demonstrate that these single-nucleotide polymorphisms (SNPs) are preferential targets of protozoa-driven selection ( P = 0.008). In particular, 19% of SNPs associated with CD have been subject to protozoa-driven selective pressure. Analysis of P values from genome-wide association studies (GWASs) and meta-analyses indicated that protozoan-selected SNPs display significantly stronger association with CD compared with nonselected variants. This same behavior was not observed for GWASs of other autoimmune diseases. Thus, we integrated selection signatures and meta-analysis results to prioritize five genic SNPs for replication in an Italian cohort. Three SNPs were significantly associated with CD risk, and combination with meta-analysis results yielded P values 〈 4 x 10 –6 . The bona fide risk alleles are located in ARHGEF2 , an interactor of NOD2, NSF , a gene involved in autophagy, and HEBP1 , encoding a possible mediator of inflammation. Pathway analysis indicated that ARHGEF2 and NSF participate in a molecular network, which also contains VAMP3 (previously associated to CD) and is centered around miR-31 (known to be disregulated in CD). Thus, we show that protozoa-driven selective pressure had a major role in shaping predisposition to CD. We next used this information for the identification of three bona fide novel susceptibility loci.