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.     

Long-lasting inhibition of EGFR autophosphorylation in A549 tumor cells by intracellular accumulation of non-covalent inhibitors

In the present study, a small set of reversible or irreversible 4-anilinoquinazoline EGFR inhibitors was tested in A549 cells at early (1 h) and late (8 h) time points after inhibitor removal from culture medium. A combination of assays was employed to explain the observed long-lasting inhibition of EGFR autophosphorylation. We found that EGFR inhibition at 8 h can be due, besides to the covalent interaction of the inhibitor with Cys797, as for PD168393 (2) and its prodrug 4, to the intracellular accumulation of non-covalent inhibitors by means of an active cell uptake, as for 5 and 6. Compounds 5–6 showed similar potency and duration of inhibition of EGFR autophosphorylation as the covalent inhibitor 2, while being devoid of reactive groups forming covalent bonds with protein thiols.     

Linkage of infantile Bartter syndrome with sensorineural deafness to chromosome 1p.

Bartter syndrome (BS) is a family of disorders manifested by hypokalemic hypochloremic metabolic alkalosis with normotensive hyperreninemic hyperaldosteronism. We evaluated a unique, inbred Bedouin kindred in which sensorineural deafness (SND) cosegregates with an infantile variant of the BS phenotype. Using a DNA-pooling strategy, we screened the human genome and successfully demonstrated linkage of this unique syndrome to chromosome 1p31. The genes for two kidney-specific chloride channels and a sodium/hydrogen antiporter, located near this region, were excluded as candidate genes. Although the search for the disease-causing gene in this family continues, this linkage further demonstrates the genetic heterogeneity of BS. In addition, the cosegregation of these phenotypes allows us to postulate that a single genetic alteration may be responsible for the SND and the BS phenotype. The identification and characterization of this gene would lead to a better understanding of the normal physiology of the kidney and the inner ear.     

Capping structures of simian virus 40 19S and 16S mRNAs.

In vivo [methyl 3H]-labeled SV40 19S and 16S mRNA species were purified and their internal methylation as well as their capping structures analyzed. SV40 viral mRNA sedimenting in the 19S region contains approximately equal proportions of m7GpppAm and m7Gppm6Am, while the 16S mRNA contains mainly m7Gpppm6Am. N6 methyl adenosine is located internally within the RNA chains of both the 19S and 16S species.     

Interplay between structure and dynamics in chitosan films investigated with solid-state NMR, dynamic mechanical analysis, and X-ray diffraction

Modern solid-state NMR techniques, combined with X-ray diffraction, revealed the molecular origin of the difference in mechanical properties of self-associated chitosan films. Films cast from acidic aqueous solutions were compared before and after neutralization, and the role of the counterion (acetate vs Cl(-)) was investigated. There is a competition between local structure and long-range order. Hydrogen bonding gives good mechanical strength to neutralized films, which lack long-range organization. The long-range structure is better defined in films cast from acidic solutions in which strong electrostatic interactions cause rotational distortion around the chitosan chains. Plasticization by acetate counterions enhances long-range molecular organization and film flexibility. In contrast, Cl(-) counterions act as a defect and impair the long-range organization by immobilizing hydration water. Molecular motion and proton exchange are restricted, resulting in brittle films despite the high moisture content.     

Site-specific excisional recombination strategies for elimination of undesirable transgenes from crop plants

A major limitation of crop biotechnology and breeding is the lack of efficient molecular technologies for precise engineering of target genomic loci. While transformation procedures have become routine for a growing number of plant species, the random introduction of complex transgenenic DNA into the plant genome by current methods generates unpredictable effects on both transgene and homologous native gene expression. The risk of transgene transfer into related plant species and consumers is another concern associated with the conventional transformation technologies. Various approaches to avoid or eliminate undesirable transgenes, most notably selectable marker genes used in plant transformation, have recently been developed. These approaches include cotransformation with two independent T-DNAs or plasmid DNAs followed by their subsequent segregation, transposon-mediated DNA elimination, and most recently, attempts to replace bacterial T-DNA borders and selectable marker genes with functional equivalents of plant origin. The use of site-specific recombination to remove undesired DNA from the plant genome and concomitantly, via excision-mediated DNA rearrangement, switch-activate by choice transgenes of agronomical, food or feed quality traits provides a versatile “transgene maintenance and control” strategy that can significantly contribute to the transfer of transgenic laboratory developments into farming practice. This review focuses on recent reports demonstrating the elimination of undesirable transgenes (essentially selectable marker and recombinase genes) from the plant genome and concomitant activation of a silent transgene (e.g., a reporter gene) mediated by different site-specific recombinases driven by constitutive or chemically, environmentally or developmentally regulated promoters. These reports indicate major progress in excision strategies which extends application of the technology from annual, sexually propagated plants towards perennial, woody and vegetatively propagated plants. Current trends and future prospects for optimization of excision-activation machinery and its practical implementation for the generation of transgenic plants and plant products free of undesired genes are discussed.     

Cloning and functional expression of alkaline alpha-galactosidase from melon fruit: similarity to plant SIP proteins uncovers a novel family of plant glycosyl hydrolases

Raffinose and stachyose are ubiquitous galactosyl-sucrose oligosaccharides in the plant kingdom which play major roles, second only to sucrose, in photoassimilate translocation and seed carbohydrate storage. These sugars are initially metabolised by alpha-galactosidases (alpha-gal). We report the cloning and functional expression of the first genes, CmAGA1 and CmAGA2, encoding for plant alpha-gals with alkaline pH optima from melon fruit (Cucumis melo L.), a raffinose and stachyose translocating species. The alkaline alpha-gal genes show very high sequence homology with a family of undefined 'seed imbibition proteins' (SIPs) which are present in a wide range of plant families. In order to confirm the function of SIP proteins, a representative SIP gene, from tomato, was expressed and shown to have alkaline alpha-gal activity. Phylogenetic analysis based on amino acid sequences shows that the family of alkaline alpha-gals shares little homology with the known prokaryotic and eukaryotic alpha-gals of glycosyl hydrolase families 27 and 36, with the exception of two cross-family conserved sequences containing aspartates which probably function in the catalytic step. This previously uncharacterised, plant-specific alpha-gal family of glycosyl hydrolases, with optimal activity at neutral-alkaline pH likely functions in key processes of galactosyl-oligosaccharide metabolism, such as during seed germination and translocation of RFO photosynthate.     

Modeling valley fever (coccidioidomycosis) incidence on the basis of climate conditions

Valley fever (coccidioidomycosis) is a disease endemic to arid regions within the Western Hemisphere, and is caused by a soil-dwelling fungus, Coccidioides immitis. Incidence data for Pima County, reported to the Arizona Department of Health Services as new cases of valley fever, were used to conduct exploratory analyses and develop monthly multivariate models of relationships between valley fever incidence and climate conditions and variability in Pima County, Arizona, USA. Bivariate and compositing analyses conducted during the exploratory portion of the study revealed that antecedent temperature and precipitation in different seasons are important predictors of incidence. These results were used in the selection of candidate variables for multivariate predictive modeling, which was designed to predict deviation from mean incidence on the basis of past, current, and forecast climate conditions. The models were specified using a backward stepwise procedure, and were most sensitive to key predictor variables in the winter season and variables that were time-lagged 1 year or more prior to the month being predicted. Model accuracy was generally moderate ( r(2) values for the monthly models, tested on independent data, ranged from 0.15 to 0.50), and months with high incidence can be predicted more accurately than months with low incidence.