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.     

Differentially Expressed miRNAs Influence Metabolic Processes in Pituitary Oncocytoma

Neurochemical Research - Spindle cell oncocytomas (SCO) of the pituitary are rare tumors accounting for 0.1–0.4% of all sellar tumors. Due to their rarity, little information is available...     

The effects of clonidine on the kidney function in the anesthetized dog

Studies were performed to determine the mechanism by which the antihypertensive agent clonidine increased urine flow. The response of the kidney has been examined in four combinations. The parameters of renal function have been compared during volume expansion by 1.5-2.0% body weight Ringer solution. In the control animals, volume expansion by 2% body weight, resulted in a slight increase in sodium excretion and urine flow. In 10 anesthetized dogs 1.0 microgram/kg/min of clonidine infused i.v. during 30 minutes (the total amount of clonidine infused was 30 micrograms/kg) decreased the arterial blood pressure from 136 +/- 13 mmHg to 127 +/- 12 mmHg and elevated urine flow from 2.95 +/- 1.65 ml/min to 4.34 +/- 1.77 ml/min while the urine osmolality diminished from 399 +/- 107 mosm/l to 265 +/- 90 mosm/l and the glomerular filtration remained constant. In 5 animals 0.1 microgram/kg/min of clonidine was infused into the left renal artery (this dose is corresponding to the renal fraction of the cardiac output) without any effects in the left kidney. 1.0 microgram/kg/min of clonidine infused directly into the left renal artery produced vasoconstriction in the ipsilateral kidney, decreased the glomerular filtration rate and the urine flow. By contrast in the right kidney the urine flow rose without hemodynamic changes, and the urine osmolality became hypoosmotic compared to the plasma. In ten dogs 1.0 microgram/kg/min of clonidine and 1 mU/kg/min of arginine-vasopressin were infused intravenously. The vasopressin infusion superimposed on the clonidine could not inhibit the increase of the urine excretion, and the fall of the urine osmolality. The results suggest that the clonidine increases the renal medullary blood flow possibly via a direct mechanism, decreases the sympathetic outflow to the kidney and via an indirect pathway, mediated by the renin-angiotensin system. The renal medullary flow increase produces a washout of the medullary osmotic gradient, and the water reabsorption diminishes.     

An efficient haplotyping method with DNA pools

Determination of haplotype frequencies (the joint distribution of genetic markers) in large population samples is a powerful tool for association studies. This is due to their greater extent of polymorphism since any two bi-allelic single nucleotide polymorphisms (SNPs) generate a potential four-allele genetic marker. Therefore, a haplotype may capture a given functional polymorphism with higher statistical power than its SNP components. The statistical estimation of haplotype frequencies, usually employed in linkage disequilibrium studies, requires individual genotyping for each SNP in the haplotype, thus making it an expensive process. In this study, we describe a new method for direct measurement of haplotype frequencies in DNA pools by allele-specific, long-range haplotype amplification. The proposed method allows the efficient determination of haplotypes composed of two SNPs in close vicinity (up to 20 kb).     

Replication delay along FRA7H, a common fragile site on human chromosome 7, leads to chromosomal instability

Common fragile sites are specific chromosomal loci that show gaps, breaks, or rearrangements in metaphase chromosomes under conditions that interfere with DNA replication. The mechanism underlying the chromosomal instability at fragile sites was hypothesized to associate with late replication time. Here, we aimed to investigate the replication pattern of the common fragile site FRA7H, encompassing 160 kb on the long arm of human chromosome 7. Using in situ hybridization on interphase nuclei, we revealed that the replication of this region is initiated relatively early, before 30% of S phase is completed. However, a high fraction ( approximately 35%) of S-phase nuclei showed allelic asynchrony, indicating that the replication of FRA7H is accomplished at different times in S phase. This allelic asynchrony is not the result of a specific replication time of each FRA7H allele. Analysis of the replication pattern of adjacent clones along FRA7H by using cell population and two-color fluorescent in situ hybridization analyses showed significant differences in the replication of adjacent clones, under normal growth condition and upon aphidicolin treatment. This pattern significantly differed from that of two nonfragile regions which showed a coordinated replication under both conditions. These results indicate that aphidicolin is enhancing an already existing difference in the replication time along the FRA7H region. Based on our replication analysis of FRA7H and on previous analysis of the common fragile site FRA3B, we suggest that delayed replication is underlying the fragility at aphidicolin-induced common fragile sites.     

Mouse inbred strain sequence information and yin-yang crosses for quantitative trait locus fine mapping

The shared ancestry of mouse inbred strains, together with the availability of sequence and phenotype information, is a resource that can be used to map quantitative trait loci (QTL). The difficulty in using only sequence information lies in the fact that in most instances the allelic state of the QTL cannot be unambiguously determined in a given strain. To overcome this difficulty, the performance of multiple crosses between various inbred strains has been proposed. Here we suggest and evaluate a general approach, which consists of crossing the two strains used initially to map the QTL and any new strain. We have termed these crosses "yin-yang," because they are complementary in nature as shown by the fact that the QTL will necessarily segregate in only one of the crosses. We used the publicly available SNP database of chromosome 16 to evaluate the mapping resolution achievable through this approach. Although on average the improvement of mapping resolution using only four inbred strains was relatively small (i.e., reduction of the QTL-containing interval by half at most), we found a great degree of variability among different regions of chromosome 16 with regard to mapping resolution. This suggests that with a large number of strains in hand, selecting a small number of strains may provide a significant contribution to the fine mapping of QTL.     

A genome-wide scan of Ashkenazi Jewish Crohn's disease suggests novel susceptibility loci

Crohn''s disease (CD) is a complex disorder resulting from the interaction of intestinal microbiota with the host immune system in genetically susceptible individuals. The largest meta-analysis of genome-wide association to date identified 71 CD–susceptibility loci in individuals of European ancestry. An important epidemiological feature of CD is that it is 2–4 times more prevalent among individuals of Ashkenazi Jewish (AJ) descent compared to non-Jewish Europeans (NJ). To explore genetic variation associated with CD in AJs, we conducted a genome-wide association study (GWAS) by combining raw genotype data across 10 AJ cohorts consisting of 907 cases and 2,345 controls in the discovery stage, followed up by a replication study in 971 cases and 2,124 controls. We confirmed genome-wide significant associations of 9 known CD loci in AJs and replicated 3 additional loci with strong signal (p<5×10⁻⁶). Novel signals detected among AJs were mapped to chromosomes 5q21.1 (rs7705924, combined p = 2×10⁻⁸; combined odds ratio OR = 1.48), 2p15 (rs6545946, p = 7×10⁻⁹; OR = 1.16), 8q21.11 (rs12677663, p = 2×10⁻⁸; OR = 1.15), 10q26.3 (rs10734105, p = 3×10⁻⁸; OR = 1.27), and 11q12.1 (rs11229030, p = 8×10⁻⁹; OR = 1.15), implicating biologically plausible candidate genes, including RPL7, CPAMD8, PRG2, and PRG3. In all, the 16 replicated and newly discovered loci, in addition to the three coding NOD2 variants, accounted for 11.2% of the total genetic variance for CD risk in the AJ population. This study demonstrates the complementary value of genetic studies in the Ashkenazim.     

Determination of zinc in vegetal tissue microsamples by platinum-wire loop in flame atomization atomic absorption spectrometry

The zinc content of 3 μL of vegetal samples (tree leaves, lichens and grape sap) atomized from a Pt-wire in the methane–air flame has been determined by atomic absorption spectrometry. The effect of gas flow rates and the atomization height in the flame on the absorption of zinc was evaluated at 213.9 nm. The best results were obtained at a height of 5 mm and gas flow rates of 200 L/h air and 26 L/h methane, respectively. The effect of Na, K, Ca, Mg, SO₄²⁻, and PO₄³⁻ on the absorption of zinc was studied too. The detection limit of 0.40 ± 0.21 ng was obtained at a significance level of 0.05, using the two-step Neyman–Pearson criterion. The zinc content of the samples has been determined with continuous nebulization and by atomization from the Pt-wire, using both the standard calibration curve and the standard addition method. The results of the two procedures agree within the determination errors.