The Variance of Identity-by-Descent Sharing in the Wright–Fisher Model

Widespread sharing of long, identical-by-descent (IBD) genetic segments is a hallmark of populations that have experienced recent genetic drift. Detection of these IBD segments has recently become feasible, enabling a wide range of applications from phasing and imputation to demographic inference. Here, we study the distribution of IBD sharing in the Wright–Fisher model. Specifically, using coalescent theory, we calculate the variance of the total sharing between random pairs of individuals. We then investigate the cohort-averaged sharing: the average total sharing between one individual and the rest of the cohort. We find that for large cohorts, the cohort-averaged sharing is distributed approximately normally. Surprisingly, the variance of this distribution does not vanish even for large cohorts, implying the existence of “hypersharing” individuals. The presence of such individuals has consequences for the design of sequencing studies, since, if they are selected for whole-genome sequencing, a larger fraction of the cohort can be subsequently imputed. We calculate the expected gain in power of imputation by IBD and subsequently in power to detect an association, when individuals are either randomly selected or specifically chosen to be the hypersharing individuals. Using our framework, we also compute the variance of an estimator of the population size that is based on the mean IBD sharing and the variance in the sharing between inbred siblings. Finally, we study IBD sharing in an admixture pulse model and show that in the Ashkenazi Jewish population the admixture fraction is correlated with the cohort-averaged sharing.IN isolated populations, even purported unrelated individuals often share genetic material that is identical-by-descent (IBD). Traditionally, the term IBD sharing referred to coancestry at a single site (or autozygosity, in the case of a diploid individual) and was widely investigated as a measure of the degree of inbreeding in a population (Hartl and Clark 2006). Recent years have brought dramatic increases in the quantity and density of available genetic data and, together with new computational tools, these data have enabled the detection of IBD sharing of entire genomic segments (see, e.g., Purcell et al. 2007; Kong et al. 2008; Albrechtsen et al. 2009; Gusev et al. 2009; Browning and Browning 2011; Carr et al. 2011; Brown et al. 2012). The availability of IBD detection tools that are efficient enough to detect shared segments in large cohorts has resulted in numerous applications, from demographic inference (Davison et al. 2009; Palamara et al. 2012) and characterization of populations (Gusev et al. 2012a) to selection detection (Albrechtsen et al. 2010), relatedness detection and pedigree reconstruction (Huff et al. 2011; Kirkpatrick et al. 2011; Stevens et al. 2011; Henn et al. 2012), prioritization of individuals for sequencing (Gusev et al. 2012b), inference of HLA type (Setty et al. 2011), detection of haplotypes associated with a disease or a trait (Akula et al. 2011; Gusev et al. 2011; Browning and Thompson 2012), imputation (Uricchio et al. 2012), and phasing (Palin et al. 2011).Recently, some of us used coalescent theory to calculate several theoretical quantities of IBD sharing under a number of demographic histories. Then, shared segments were detected in real populations, and their demographic histories were inferred (Palamara et al. 2012). Here, we expand upon Palamara et al. (2012) to investigate additional aspects of the stochastic variation in IBD sharing. Specifically, we provide a precise calculation for the variance of the total sharing in the Wright–Fisher model, either between a random pair of individuals or between one individual and all others in the cohort.Understanding the variation in IBD sharing is an important theoretical characterization of the Wright–Fisher model, and additionally, it has several practical applications. For example, it can be used to calculate the variance of an estimator of the population size that is based on the sharing between random pairs. In a different domain, the variance in IBD sharing is needed to accurately assess strategies for sequencing study design, specifically, in prioritization of individuals to be sequenced. This is because imputation strategies use IBD sharing between sequenced individuals and genotyped, not-sequenced individuals to increase the number of effective sequences analyzed in the association study (Palin et al. 2011; Gusev et al. 2012b; Uricchio et al. 2012).In the remainder of this article, we first review the derivation of the mean fraction of the genome shared between two individuals (Palamara et al. 2012). We then calculate the variance of this quantity, using coalescent theory with recombination. We provide a number of approximations, one of which results in a surprisingly simple expression, which is then generalized to a variable population size and to the sharing of segments in a length range. We also numerically investigate the pairwise sharing distribution and provide an approximate fit. We then turn to the average total sharing between each individual and the entire cohort. We show that this quantity, which we term the cohort-averaged sharing, is approximately normally distributed, but is much wider than naively expected, implying the existence of hypersharing individuals. We consider several applications: the number of individuals needed to be sequenced to achieve a certain imputation power and the implications to disease mapping, inference of the population size based on the total sharing, and the variance of the sharing between siblings. We finally calculate the mean and the variance of the sharing in an admixture pulse model and show numerically that admixture results in a broader than expected cohort-averaged sharing. Therefore, large variance of the cohort-averaged sharing can indicate admixture. In the Ashkenazi Jewish population, we show that the cohort-averaged sharing is strongly anticorrelated with the fraction of European ancestry.     

Tyrosinase and Layer-by-Layer supported tyrosinases in the synthesis of lipophilic catechols with antiinfluenza activity

Catechol derivatives with lipophilic properties have been selectively synthesized by tyrosinase in high yield avoiding long and tedious protection/deprotection steps usually required in traditional procedures. The synthesis was effective also with immobilized tyrosinase able to perform for more runs. The novel catechols were evaluated against influenza A virus, that continue to represent a severe threat worldwide. A significant antiviral activity was observed in derivatives characterized by antioxidant activity and long carbon alkyl side-chains, suggesting the possibility of a new inhibition mechanism based on both redox and lipophilic properties.     

The effects of demographic stochasticity and parameter uncertainty on predicting the establishment of introduced species

Invasive species are a serious threat to biodiversity worldwide and predicting whether an introduced species will first establish and then become invasive can be useful to preserve ecosystem services. Establishment is influenced by multiple factors, such as the interactions between the introduced individuals and the resident community, and demographic and environmental stochasticity. Field observations are often incomplete or biased. This, together with an imperfect knowledge of the ecological traits of the introduced species, makes the prediction of establishment challenging. Methods that consider the combined effects of these factors on our ability to predict the establishment of an introduced species are currently lacking. We develop an inference framework to assess the combined effects of demographic stochasticity and parameter uncertainty on our ability to predict the probability of establishment following the introduction of a small number of individuals. We find that even moderate levels of demographic stochasticity influence both the probability of establishment, and, crucially, our ability to correctly predict that probability. We also find that estimation of the demographic parameters of an introduced species is fundamental to obtain precise estimates of the interaction parameters. For typical values of demographic stochasticity, the drop in our ability to predict an establishment can be 30% when having priors on the demographic parameters compared to having their accurate values. The results from our study illustrate how demographic stochasticity may bias the prediction of the probability of establishment. Our method can be applied to estimate probability of establishment of introduced species in field scenarios, where time series data and prior information on the demographic traits of the introduced species are available.     

The conformation of peptide thymosin α1 in solution and in a membrane-like environment by circular dichroism and NMR spectroscopy. a possible model for its interaction with the lymphocyte membrane

The 28-residue peptide thymosin α1 was studied by circular dichroism and two-dimensional NMR. Circular dichroism indicates that thymosin α1 in water solution does not assume a preferred conformation, while in the presence of small unilamellar vesicles of dimiristoylphosphatidylcholine and dimiristoylphosphatidic acid (10:1) and in sodium dodecyl sulphate, it assumes a partly structured conformation. Presence of zinc ions produces similar effects. In a more hydrophobic environment like a solution of a mixed solvent water-2,2,2 trifluoroethanol, it adopts a structured conformation. NMR spectra indicated that in this mixture as solvent, thymosin α1 has a structure characterized by two regions. A β-turn is present between residue 5 and residue 8, while the region between residues 17 and 24 shows an α helix conformation. These changes of conformation in different environments may be considered structural requirements in the steps of its interaction with the lymphocyte membrane. In fact, these conformational changes may correspond to the first event of the mechanism of lymphocyte activation in the immune response modulation by thymosin α1.     

Progress in modeling quality in aquaculture: an application of the Self‐Organizing Map to the study of skeletal anomalies and meristic counts in gilthead seabream (Sparus aurata,L. 1758)

One of the most common drawbacks of artificial life conditions imposed by aquaculture is the quite high presence of skeletal anomalies (SAs) in reared fish, which reduce both functional performances and marketing image/commercial value of the reared lots. Thus, skeletal malformations and their incidence are one of the most important factors affecting fish farmer’s production costs, and several efforts have been due to develop appropriate tools in detecting patterns of co‐variation among rearing parameters and fish quality. In this paper we explore the advantages of using Self‐Organized Maps (SOMs) when dealing with the analysis of correlations between the pattern of SA presence and rearing parameters in gilthead seabream (Sparus aurata L.), that is a largely reared fish of high commercial value. SOM, which is one of the best known neural networks with unsupervised learning rules, were applied to develop a model of the occurrence of SAs, both in terms of type and quantity, in seabream lots from different rearing approaches (extensive, semi‐intensive and intensive). The trained SOMs classified lots according to the variation observed in the different weights of SAs, but also allows the detection of a series of correspondence, namely between: (i) the patter of SAs occurrence and the different rearing approach currently used in seabream aquaculture; and (ii) the total SAs incidence and the variability of meristic counts, represent a completely independent dataset. Mesocosms resulted the best rearing approach to produce wild‐like fish, whereas intensive rearing is characterized by the large presence of SA. Globally, results suggested that this approach is reliable to be used for estimate the distance between aquaculture products and the wild‐like phenotype used as quality reference.     

Physicochemical characterization of casein phosphopeptide-amorphous calcium phosphate nanocomplexes

Milk caseins stabilize calcium and phosphate ions and make them available to the neonate. Tryptic digestion of the caseins yields phosphopeptides from their polar N-terminal regions that contain clusters of phosphorylated seryl residues. These phosphoseryl clusters have been hypothesized to be responsible for the interaction between the caseins and calcium phosphate that lead to the formation of casein micelles. The casein phosphopeptides stabilize calcium and phosphate ions through the formation of complexes. The calcium phosphate in these complexes is biologically available for intestinal absorption and remineralization of subsurface lesions in tooth enamel. We have studied the structure of the complexes formed by the casein phosphopeptides with calcium phosphate using a range of physicochemical techniques including x-ray powder diffraction, scanning electron microscopy, transmission electron microscopy, and equilibrium binding analyses. The amorphous nature of the calcium phosphate phase was confirmed by two independent methods: x-ray powder diffraction and selected area diffraction. In solution, the ion activity product of a basic amorphous calcium phosphate phase was the only ion product that was a function of bound phosphate independent of pH, consistent with basic amorphous calcium phosphate being the phase stabilized by the casein phosphopeptides. Detailed investigations of calcium and calcium phosphate binding using a library of synthetic homologues and analogues of the casein phosphopeptides have revealed that although the fully phosphorylated seryl-cluster motif is pivotal for the interaction with calcium and phosphate, other factors are also important. In particular, calcium binding and calcium phosphate stabilization by the peptides was influenced by peptide net charge, length, and sequence.     

Infectious Agents and Neurodegeneration

A growing body of epidemiologic and experimental data point to chronic bacterial and viral infections as possible risk factors for neurodegenerative diseases, including Alzheimer??s disease, Parkinson??s disease and amyotrophic lateral sclerosis. Infections of the central nervous system, especially those characterized by a chronic progressive course, may produce multiple damage in infected and neighbouring cells. The activation of inflammatory processes and host immune responses cause chronic damage resulting in alterations of neuronal function and viability, but different pathogens can also directly trigger neurotoxic pathways. Indeed, viral and microbial agents have been reported to produce molecular hallmarks of neurodegeneration, such as the production and deposit of misfolded protein aggregates, oxidative stress, deficient autophagic processes, synaptopathies and neuronal death. These effects may act in synergy with other recognized risk factors, such as aging, concomitant metabolic diseases and the host??s specific genetic signature. This review will focus on the contribution given to neurodegeneration by herpes simplex type-1, human immunodeficiency and influenza viruses, and by Chlamydia pneumoniae.     

Leveraging Distant Relatedness to Quantify Human Mutation and Gene-Conversion Rates

The rate at which human genomes mutate is a central biological parameter that has many implications for our ability to understand demographic and evolutionary phenomena. We present a method for inferring mutation and gene-conversion rates by using the number of sequence differences observed in identical-by-descent (IBD) segments together with a reconstructed model of recent population-size history. This approach is robust to, and can quantify, the presence of substantial genotyping error, as validated in coalescent simulations. We applied the method to 498 trio-phased sequenced Dutch individuals and inferred a point mutation rate of 1.66 × 10⁻⁸ per base per generation and a rate of 1.26 × 10⁻⁹ for <20 bp indels. By quantifying how estimates varied as a function of allele frequency, we inferred the probability that a site is involved in non-crossover gene conversion as 5.99 × 10⁻⁶. We found that recombination does not have observable mutagenic effects after gene conversion is accounted for and that local gene-conversion rates reflect recombination rates. We detected a strong enrichment of recent deleterious variation among mismatching variants found within IBD regions and observed summary statistics of local sharing of IBD segments to closely match previously proposed metrics of background selection; however, we found no significant effects of selection on our mutation-rate estimates. We detected no evidence of strong variation of mutation rates in a number of genomic annotations obtained from several recent studies. Our analysis suggests that a mutation-rate estimate higher than that reported by recent pedigree-based studies should be adopted in the context of DNA-based demographic reconstruction.     

Synthetic prions with novel strain-specified properties

Prions are infectious proteins that possess multiple self-propagating structures. The information for strains and structural specific barriers appears to be contained exclusively in the folding of the pathological isoform, PrP^Sc. Many recent studies determined that de novo prion strains could be generated in vitro from the structural conversion of recombinant (rec) prion protein (PrP) into amyloidal structures. Our aim was to elucidate the conformational diversity of pathological recPrP amyloids and their biological activities, as well as to gain novel insights in characterizing molecular events involved in mammalian prion conversion and propagation. To this end we generated infectious materials that possess different conformational structures. Our methodology for the prion conversion of recPrP required only purified rec full-length mouse (Mo) PrP and common chemicals. Neither infected brain extracts nor amplified PrP^Sc were used. Following two different in vitro protocols recMoPrP converted to amyloid fibrils without any seeding factor. Mouse hypothalamic GT1 and neuroblastoma N2a cell lines were infected with these amyloid preparations as fast screening methodology to characterize the infectious materials. Remarkably, a large number of amyloid preparations were able to induce the conformational change of endogenous PrP^C to harbor several distinctive proteinase-resistant PrP forms. One such preparation was characterized in vivo habouring a synthetic prion with novel strain specified neuropathological and biochemical properties.     

"Shock and kill" effects of class I-selective histone deacetylase inhibitors in combination with the glutathione synthesis inhibitor buthionine sulfoximine in cell line models for HIV-1 quiescence

Background

Bovine Leukemia virus (BLV) is a deltaretrovirus that induces lymphoproliferation and leukemia in ruminants. In ex vivo cultures of B lymphocytes isolated from BLV-infected sheep show that spontaneous apoptosis is reduced. Here, we investigated the involvement of reactive oxygen species (ROS) in this process.

Results

We demonstrate that (i) the levels of ROS and a major product of oxidative stress (8-OHdG) are reduced, while the thioredoxin antioxidant protein is highly expressed in BLV-infected B lymphocytes, (ii) induction of ROS by valproate (VPA) is pro-apoptotic, (iii) inversely, the scavenging of ROS with N-acetylcysteine inhibits apoptosis, and finally (iv) the levels of ROS inversely correlate with the proviral loads.

Conclusion

Together, these observations underline the importance of ROS in the mechanisms of inhibition of apoptosis linked to BLV infection.