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.     

Identification of the infectious source of an unusual outbreak of histoplasmosis,in a hotel in Acapulco,state of Guerrero,Mexico

Three isolates of Histoplasma capsulatum were identified from mice lung, liver, and spleen inoculated with soil samples of the X hotel’s ornamental potted plants that had been fertilized with organic material known as compost. The presence of H. capsulatum in the original compost was detected using the dot-enzyme-linked immunosorbent assay. Nested-PCR, using a specific protein Hcp100 coding gene sequence, confirmed the fungal identification associated with an unusual histoplasmosis outbreak in Acapulco. Although, diversity between the H. capsulatum isolate from the hotel and some clinical isolates from Guerrero (positive controls) was observed using random amplification of polymorphic DNA based-PCR, sequence analyses of H-anti and ole fragment genes revealed a high homology (92–99%) between them.     

Foraging area fidelity for Kemp's ridleys in the Gulf of Mexico

For many marine species, locations of key foraging areas are not well defined. We used satellite telemetry and switching state‐space modeling (SSM) to identify distinct foraging areas used by Kemp's ridley turtles (Lepidochelys kempii) tagged after nesting during 1998–2011 at Padre Island National Seashore, Texas, USA (PAIS; N = 22), and Rancho Nuevo, Tamaulipas, Mexico (RN; N = 9). Overall, turtles traveled a mean distance of 793.1 km (±347.8 SD) to foraging sites, where 24 of 31 turtles showed foraging area fidelity (FAF) over time (N = 22 in USA, N = 2 in Mexico). Multiple turtles foraged along their migratory route, prior to arrival at their “final” foraging sites. We identified new foraging “hotspots” where adult female Kemp's ridley turtles spent 44% of their time during tracking (i.e., 2641/6009 tracking days in foraging mode). Nearshore Gulf of Mexico waters served as foraging habitat for all turtles tracked in this study; final foraging sites were located in water <68 m deep and a mean distance of 33.2 km (±25.3 SD) from the nearest mainland coast. Distance to release site, distance to mainland shore, annual mean sea surface temperature, bathymetry, and net primary production were significant predictors of sites where turtles spent large numbers of days in foraging mode. Spatial similarity of particular foraging sites selected by different turtles over the 13‐year tracking period indicates that these areas represent critical foraging habitat, particularly in waters off Louisiana. Furthermore, the wide distribution of foraging sites indicates that a foraging corridor exists for Kemp's ridleys in the Gulf. Our results highlight the need for further study of environmental and bathymetric components of foraging sites and prey resources contained therein, as well as international cooperation to protect essential at‐sea foraging habitats for this imperiled species.     

Comparing the Response of Birds and Butterflies to Vegetation-Based Mountain Ecotones Using Boundary Detection Approaches

Mountains provide an opportunity to examine changes in biodiversity across environmental gradients and areas of transition (ecotones). Mountain ecotones separate vegetation belts. Here, we aimed to examine whether transition areas for birds and butterflies spatially correspond with ecotones between three previously described altitudinal vegetation belts on Mt. Hermon, northern Israel. These include the Mediterranean Maquis, xero-montane open forest and Tragacanthic mountain steppe vegetation belts. We sampled the abundance of bird and butterfly species in 34 sampling locations along an elevational gradient between 500 and 2200 m. We applied wombling, a boundary-detection technique, which detects rapid changes in a continuous variable, in order to locate the transition areas for bird and butterfly communities and compare the location of these areas with the location of vegetation belts as described in earlier studies of Mt. Hermon. We found some correspondence between the areas of transition of both bird and butterfly communities and the ecotones between vegetation belts. For birds and butterflies, important transitions occurred at the lower vegetation ecotone between Mediterranean maquis and the xero-montane open forest vegetation belts, and between the xero-montane open forest and the mountain steppe Tragacanthic belts. While patterns of species turnover with elevation were similar for birds and butterflies, the change in species richness and diversity with elevation differed substantially between the two taxa. Birds and butterflies responded quite similarly to the elevational gradient and to the shift between vegetation belts in terms of species turnover rates. While the mechanisms generating these patterns may differ, the resulting areas of peak turnover in species show correspondence among three different taxa (plants, birds and butterflies).     

Failure properties of vena cava tissue due to deep penetration during filter insertion

In this work, we use an in-vitro mechanical test to explore the resistance of biaxially stretched vena cava tissue against deep perforation and a methodology which integrates experimental and numerical modeling to identify constitutive fracture properties of the vena cava. Six sheep vena cava were harvested just after killing, and cyclic uniaxial tension tests in longitudinal and circumferential directions and biaxial deep penetration tests were performed. After that, we use a nonlinear finite element model to simulate in vitro penetration of the cava tissue in order to fit the fracture properties under penetration of the vena cava by defining a cohesive fracture zone. An iterative process was developed in order to fit the fracture properties of the vena cava using the previously obtained experimental results. The proposed solutions were obtained with fracture energy of 0.22 or 0.33 N/mm. In comparison with the experimental data, the simulation using \(\delta _{0}=0.01\,\hbox {mm}\), \(\delta _{r}=0.35\,\hbox {mm}\), and \(K=220\, \hbox {N}/\hbox {mm}^{3}\) parameters (\(F_{\hbox {max}}=0.92\)) is in good agreement with results from penetration experiments of cava tissue. It is noticeable that the parameter estimation process of the fracture behavior is more accurate than the estimation process of the elastic behavior for the toe region of the curve.     

Simplifying Hansen Solubility Parameters for Complex Edible Fats and Oils

Hansen solubility parameters (HSPs), often used to predict the miscibility between two compounds, are an alternative tool in evaluating the ability of the solvent to interact via dispersion, dipole-dipole, and hydrogen bonding interactions. The aim of this paper is to find a simple way to calculate HSPs for complex mixtures of triglycerides (TAGs). HSPs were calculated using two approaches: the first assumes that the contributions to the dispersion, dipole-dipole, and hydrogen bonding interactions may be subdivided into larger functional moieties (i.e., fatty acids and fatty acid methyl esters) that are additive, while the second approach assumes that vegetable oils are comprised of mixtures of simple TAGs in the same mass fractions as the fatty acids. The HSPs obtained using the two approaches are compared to reference values determined using the “Hansen Solubility Parameters in Practice” software (HSPiP) considering the complex TAG profile for each vegetable oil.HSPs for vegetable oils, obtained with the HSPiP software, did not correspond well to the HSPs obtained from the group contribution approach, when using fatty acids, fatty acids + glycerol or fatty acid methyl esters. In contrast, the HSPs calculated for vegetable oils, assuming that all TAGs are simple and in the same mass fractions as the fatty acids, provide similar values to the HSPs obtained from the HSPiP software. Therefore, it is possible to calculate the HSPs for complex oils by simply knowing the fatty acid composition. Knowledge of HSPs may be used to rationalize the ability of certain low molecular weight molecules to develop organogels in vegetable oils as well as the crystallization of triglycerides.     

BAG3 regulates total MAP1LC3B protein levels through a translational but not transcriptional mechanism

Autophagy is mainly regulated by post-translational and lipid modifications of ATG proteins. In some scenarios, the induction of autophagy is accompanied by increased levels of certain ATG mRNAs such as MAP1LC3B/LC3B, ATG5 or ATG12. However, little is known about the regulation of ATG protein synthesis at the translational level. The cochaperone of the HSP70 system BAG3 (BCL2-associated athanogene 3) has been associated to LC3B lipidation through an unknown mechanism. In the present work, we studied how BAG3 controls autophagy in HeLa and HEK293 cells. Our results showed that BAG3 regulates the basal amount of total cellular LC3B protein by controlling its mRNA translation. This effect was apparently specific to LC3B because other ATG protein levels were not affected. BAG3 knockdown did not affect LC3B lipidation induced by nutrient deprivation or proteasome inhibition. We concluded that BAG3 maintains the basal amount of LC3B protein by controlling the translation of its mRNA in HeLa and HEK293 cells.     

Hawksbill turtle terra incognita: conservation genetics of?eastern Pacific rookeries

Prior to 2008 and the discovery of several important hawksbill turtle (Eretmochelys imbricata) nesting colonies in the EP (Eastern Pacific), the species was considered virtually absent from the region. Research since that time has yielded new insights into EP hawksbills, salient among them being the use of mangrove estuaries for nesting. These recent revelations have raised interest in the genetic characterization of hawksbills in the EP, studies of which have remained lacking to date. Between 2008 and 2014, we collected tissue samples from 269 nesting hawksbills at nine rookeries across the EP and used mitochondrial DNA sequences (766 bp) to generate the first genetic characterization of rookeries in the region. Our results inform genetic diversity, population differentiation, and phylogeography of the species. Hawksbills in the EP demonstrate low genetic diversity: We identified a total of only seven haplotypes across the region, including five new and two previously identified nesting haplotypes (pooled frequencies of 58.4% and 41.6%, respectively), the former only evident in Central American rookeries. Despite low genetic diversity, we found strong stock structure between the four principal rookeries, suggesting the existence of multiple populations and warranting their recognition as distinct management units. Furthermore, haplotypes EiIP106 and EiIP108 are unique to hawksbills that nest in mangrove estuaries, a behavior found only in hawksbills along Pacific Central America. The detected genetic differentiation supports the existence of a novel mangrove estuary “reproductive ecotype” that may warrant additional conservation attention. From a phylogeographic perspective, our research indicates hawksbills colonized the EP via the Indo‐Pacific, and do not represent relict populations isolated from the Atlantic by the rising of the Panama Isthmus. Low overall genetic diversity in the EP is likely the combined result of few rookeries, extremely small reproductive populations and evolutionarily recent colonization events. Additional research with larger sample sizes and variable markers will help further genetic understanding of hawksbill turtles in the EP.