Inferring social network structure from bacterial sequence data

Using DNA sequence data from pathogens to infer transmission networks has traditionally been done in the context of epidemics and outbreaks. Sequence data could analogously be applied to cases of ubiquitous commensal bacteria; however, instead of inferring chains of transmission to track the spread of a pathogen, sequence data for bacteria circulating in an endemic equilibrium could be used to infer information about host contact networks. Here, we show--using simulated data--that multilocus DNA sequence data, based on multilocus sequence typing schemes (MLST), from isolates of commensal bacteria can be used to infer both local and global properties of the contact networks of the populations being sampled. Specifically, for MLST data simulated from small-world networks, the small world parameter controlling the degree of structure in the contact network can robustly be estimated. Moreover, we show that pairwise distances in the network--degrees of separation--correlate with genetic distances between isolates, so that how far apart two individuals in the network are can be inferred from MLST analysis of their commensal bacteria. This result has important consequences, and we show an example from epidemiology: how this result could be used to test for infectious origins of diseases of unknown etiology.     

Ancient DNA reveals kinship burial patterns of a pre-Columbian Andean community

ABSTRACT: BACKGROUND: A detailed genetic study of the pre-Columbian population inhabiting the Tompullo 2 archaeological site (department Arequipa, Peru) was undertaken to resolve the kin relationships between individuals buried in six different chullpas. Kin relationships were an important factor shaping the social organization in the pre-Columbian Andean communities, centering on the ayllu, a group of relatives that shared a common land and responsibilities. The aim of this study was to evaluate whether this Andean model of a social organization had an influence on mortuary practices, in particular to determine whether chullpas served as family graves. RESULTS: The remains of forty-one individuals were analyzed with both uniparental (mtDNA, Y-chromosome) and biparental (autosomal microsatellites) markers. Reproducible HVRI sequences, autosomal and Y chromosomal STR profiles were obtained for 24, 16 and 11 individuals, respectively. Mitochondrial DNA diversity was comparable to that of ancient and contemporary Andean populations. The Tompullo 2 population exhibited the closest relationship with the modern population from the same region. A kinship analysis revealed complex pattern of relations within and between the graves. However mean relatedness coefficients regarding the pairs of individuals buried in the same grave were significantly higher than those regarding pairs buried in different graves. The Y chromosome profiles of 11 males suggest that only members of one male line were buried in the same grave. CONCLUSIONS: Genetic investigation of the population that inhabited Tompullo 2 site shows continuity between pre-Columbian and modern Native Amerindian populations inhabiting the Arequipa region. This suggests that no major demographic processes have influenced the mitochondrial DNA diversity of these populations during the past five hundred years. The kinship analysis involving uni- and biparental markers suggests that the community that inhabited the Tompullo 2 site was organized into extended family groups that were buried in different graves. This finding is in congruence with known models of social organization of Andean communities.     

Fluorescence-based sensors to monitor localization and functions of linear and k63-linked ubiquitin chains in cells

Ubiquitin chains modify a major subset of the proteome, but detection of ubiquitin signaling dynamics and localization is limited due to a lack of appropriate tools. Here, we employ ubiquitin-binding domain (UBD)-based fluorescent sensors to monitor linear and K63-linked chains in?vitro and in?vivo. We utilize the UBD in NEMO and ABIN (UBAN) for detection of linear chains, and RAP80 ubiquitin-interacting motif (UIM) and TAB2 Npl4 zinc finger (NZF) domains to detect K63 chains. Linear and K63 sensors decorated the ubiquitin coat surrounding cytosolic Salmonella during bacterial autophagy, whereas K63 sensors selectively monitored Parkin-induced mitophagy and DNA damage responses in fixed and living cells. In addition, linear and K63 sensors could be used to monitor endogenous signaling pathways, as demonstrated by their ability to differentially interfere with TNF- and IL-1-induced NF-κB pathway. We propose that UBD-based biosensors could serve as prototypes to track and trace other chain types and ubiquitin-like signals in?vivo.     

Within-host competition and drug resistance in the human malaria parasite Plasmodium falciparum

Infections with the malaria parasite Plasmodium falciparum typically comprise multiple strains, especially in high-transmission areas where infectious mosquito bites occur frequently. However, little is known about the dynamics of mixed-strain infections, particularly whether strains sharing a host compete or grow independently. Competition between drug-sensitive and drug-resistant strains, if it occurs, could be a crucial determinant of the spread of resistance. We analysed 1341 P. falciparum infections in children from Angola, Ghana and Tanzania and found compelling evidence for competition in mixed-strain infections: overall parasite density did not increase with additional strains, and densities of individual chloroquine-sensitive (CQS) and chloroquine-resistant (CQR) strains were reduced in the presence of competitors. We also found that CQR strains exhibited low densities compared with CQS strains (in the absence of chloroquine), which may underlie observed declines of chloroquine resistance in many countries following retirement of chloroquine as a first-line therapy. Our observations support a key role for within-host competition in the evolution of drug-resistant malaria. Malaria control and resistance-management efforts in high-transmission regions may be significantly aided or hindered by the effects of competition in mixed-strain infections. Consideration of within-host dynamics may spur development of novel strategies to minimize resistance while maximizing the benefits of control measures.     

Comparison of three standard anatomical reference frames for the tibia-fibula complex

Definition of anatomical reference frames is necessary both for in vitro biomechanical testing, and for in vivo human movement analyses. Different reference frames have been proposed in the literature for the lower limb, and in particular for the tibia–fibula complex. The scope of this work was to compare the three most commonly referred proposals (proposed by [Ruff, C.B., Hayes, W.C., 1983. Cross-sectional geometry at Pecos Pueblo femora and tibiae —A biomechanical investigation: I. method and general patterns of variation. American Journal of Physical Anthropology 60, pp. 359–381.], by [Cappozzo, A., Catani, F., Della Croce, U., Leardini, A., 1995. Position and orientation in space of bones during movement: anatomical frame definition and determination. Clinical Biomechanics (Bristol, Avon) 10, pp. 171–178.], and by the Standardization and Terminology Committee of the International Society of Biomechanics, [Wu, G., Siegler, S., Allard, P., Kirtley, C., Leardini, A., Rosenbaum, D., Whittle, M., D'Lima, D.D., Cristofolini, L., Witte, H., Schmid, O., Stokes, I., 2002. ISB recommendation on definitions of joint coordinate system of various joints for reporting of human joint motion—part I: ankle, hip and spine. International Society of Biomechanics. Journal of Biomechanics 35, pp. 543–548.]). These three frames were identified on six cadaveric tibia–fibula specimens based on the relevant anatomical landmarks, using a high-precision digitizer. The intra-operator (ten repetitions) and inter-operator (three operators) repeatability were investigated in terms of reference frame orientation. The three frames had similar intra-operator repeatability. The reference frame proposed by Ruff et al. had a better inter-operator repeatability (this must be put in relation with the original context of interest, i.e. in vitro measurements on dissected bones). The reference frames proposed by Ruff et al. and by ISB had a similar alignment; the frame proposed by Cappozzo et al. was considerably externally rotated and flexed with respect to the other two. Thus, the reference frame proposed by Ruff et al. is preferable when the full bone surface is accessible (typically during in vitro tests). Conversely, no advantage in terms of repeatability seems to exist between the reference frames proposed by Cappozzo et al. and ISB.     

Is the hydrophobic core a universal structural element in proteins?

The hydrophobic core, when subjected to analysis based on the fuzzy oil drop model, appears to be a universal structural component of proteins irrespective of their secondary, supersecondary, and tertiary conformations. A study has been performed on a set of nonhomologous proteins representing a variety of CATH categories. The presence of a well-ordered hydrophobic core has been confirmed in each case, regardless of the protein’s biological function, chain length or source organism. In light of fuzzy oil drop (FOD) analysis, various supersecondary forms seem to share a common structural factor in the form of a hydrophobic core, emerging either as part of the whole protein or a specific domain. The variable status of individual folds with respect to the FOD model reflects their propensity for conformational changes, frequently associated with biological function. Such flexibility is expressed as variable stability of the hydrophobic core, along with specific encoding of potential conformational changes which depend on the properties of helices and β-folds.     

Joint non-uniform sampling of all incremented time delays for quicker acquisition in protein relaxation studies

NMR relaxometry plays crucial role in studies of protein dynamics. The measurement of longitudinal and transverse relaxation rates of \(^{15}\)N is the main source of information on backbone motions. However, even the most basic approach exploiting a series of \(^{15}\)N HSQC spectra can require several hours of measurement time. Standard non-uniform sampling (NUS), i.e. random under-sampling of indirect time domain, typically cannot reduce this by more than 2–4\(\times\) due to relatively low “compressibility” of these spectra. In this paper we propose an extension of NUS to relaxation delays. The two-dimensional space of \(t_1\)/\(t_{relax}\) is sampled in a way similar to NUS of \(t_1\)/\(t_2\) domain in 3D spectra. The signal is also processed in a way similar to that known from 3D NUS spectra i.e. using one of the most popular compressed sensing algorithms, iterative soft thresholding. The 2D Fourier transform matrix is replaced with mixed inverse Laplace-Fourier transform matrix. The peak positions in resulting 3D spectrum are characterized by two frequency coordinates and relaxation rate and thus no additional fitting of exponential curves is required. The method is tested on three globular proteins, providing satisfactory results in a time corresponding to acquisition of two conventional \(^{15}\)N HSQC spectra.