Ambiguity in Social Network Data for Presence,Sensitive-Attribute,Degree and Relationship Privacy Protection

Maintaining privacy in network data publishing is a major challenge. This is because known characteristics of individuals can be used to extract new information about them. Recently, researchers have developed privacy methods based on k-anonymity and l-diversity to prevent re-identification or sensitive label disclosure through certain structural information. However, most of these studies have considered only structural information and have been developed for undirected networks. Furthermore, most existing approaches rely on generalization and node clustering so may entail significant information loss as all properties of all members of each group are generalized to the same value. In this paper, we introduce a framework for protecting sensitive attribute, degree (the number of connected entities), and relationships, as well as the presence of individuals in directed social network data whose nodes contain attributes. First, we define a privacy model that specifies privacy requirements for the above private information. Then, we introduce the technique of Ambiguity in Social Network data (ASN) based on anatomy, which specifies how to publish social network data. To employ ASN, individuals are partitioned into groups. Then, ASN publishes exact values of properties of individuals of each group with common group ID in several tables. The lossy join of those tables based on group ID injects uncertainty to reconstruct the original network. We also show how to measure different privacy requirements in ASN. Simulation results on real and synthetic datasets demonstrate that our framework, which protects from four types of private information disclosure, preserves data utility in tabular, topological and spectrum aspects of networks at a satisfactory level.     

Comparing sources of mobility for modelling the epidemic spread of Zika virus in Colombia

Timely, accurate, and comparative data on human mobility is of paramount importance for epidemic preparedness and response, but generally not available or easily accessible. Mobile phone metadata, typically in the form of Call Detail Records (CDRs), represents a powerful source of information on human movements at an unprecedented scale. In this work, we investigate the potential benefits of harnessing aggregated CDR-derived mobility to predict the 2015-2016 Zika virus (ZIKV) outbreak in Colombia, when compared to other traditional data sources. To simulate the spread of ZIKV at sub-national level in Colombia, we employ a stochastic metapopulation epidemic model for vector-borne diseases. Our model integrates detailed data on the key drivers of ZIKV spread, including the spatial heterogeneity of the mosquito abundance, and the exposure of the population to the virus due to environmental and socio-economic factors. Given the same modelling settings (i.e. initial conditions and epidemiological parameters), we perform in-silico simulations for each mobility network and assess their ability in reproducing the local outbreak as reported by the official surveillance data. We assess the performance of our epidemic modelling approach in capturing the ZIKV outbreak both nationally and sub-nationally. Our model estimates are strongly correlated with the surveillance data at the country level (Pearson’s r = 0.92 for the CDR-informed network). Moreover, we found strong performance of the model estimates generated by the CDR-informed mobility networks in reproducing the local outbreak observed at the sub-national level. Compared to the CDR-informed networks, the performance of the other mobility networks is either comparatively similar or substantially lower, with no added value in predicting the local epidemic. This suggests that mobile phone data captures a better picture of human mobility patterns. This work contributes to the ongoing discussion on the value of aggregated mobility estimates from CDRs data that, with appropriate data protection and privacy safeguards, can be used for social impact applications and humanitarian action.     

Reduced representation bisulfite sequencing for comparative high-resolution DNA methylation analysis

We describe a large-scale random approach termed reduced representation bisulfite sequencing (RRBS) for analyzing and comparing genomic methylation patterns. BglII restriction fragments were size-selected to 500–600 bp, equipped with adapters, treated with bisulfite, PCR amplified, cloned and sequenced. We constructed RRBS libraries from murine ES cells and from ES cells lacking DNA methyltransferases Dnmt3a and 3b and with knocked-down (kd) levels of Dnmt1 (Dnmt[1^kd,3a^−/−,3b^−/−]). Sequencing of 960 RRBS clones from Dnmt[1^kd,3a^−/−,3b^−/−] cells generated 343 kb of non-redundant bisulfite sequence covering 66212 cytosines in the genome. All but 38 cytosines had been converted to uracil indicating a conversion rate of >99.9%. Of the remaining cytosines 35 were found in CpG and 3 in CpT dinucleotides. Non-CpG methylation was >250-fold reduced compared with wild-type ES cells, consistent with a role for Dnmt3a and/or Dnmt3b in CpA and CpT methylation. Closer inspection revealed neither a consensus sequence around the methylated sites nor evidence for clustering of residual methylation in the genome. Our findings indicate random loss rather than specific maintenance of methylation in Dnmt[1^kd,3a^−/−,3b^−/−] cells. Near-complete bisulfite conversion and largely unbiased representation of RRBS libraries suggest that random shotgun bisulfite sequencing can be scaled to a genome-wide approach.     

Identifying Fishes through DNA Barcodes and Microarrays

Background

International fish trade reached an import value of 62.8 billion Euro in 2006, of which 44.6% are covered by the European Union. Species identification is a key problem throughout the life cycle of fishes: from eggs and larvae to adults in fisheries research and control, as well as processed fish products in consumer protection.

Methodology/Principal Findings

This study aims to evaluate the applicability of the three mitochondrial genes 16S rRNA (16S), cytochrome b (cyt b), and cytochrome oxidase subunit I (COI) for the identification of 50 European marine fish species by combining techniques of “DNA barcoding” and microarrays. In a DNA barcoding approach, neighbour Joining (NJ) phylogenetic trees of 369 16S, 212 cyt b, and 447 COI sequences indicated that cyt b and COI are suitable for unambiguous identification, whereas 16S failed to discriminate closely related flatfish and gurnard species. In course of probe design for DNA microarray development, each of the markers yielded a high number of potentially species-specific probes in silico, although many of them were rejected based on microarray hybridisation experiments. None of the markers provided probes to discriminate the sibling flatfish and gurnard species. However, since 16S-probes were less negatively influenced by the “position of label” effect and showed the lowest rejection rate and the highest mean signal intensity, 16S is more suitable for DNA microarray probe design than cty b and COI. The large portion of rejected COI-probes after hybridisation experiments (>90%) renders the DNA barcoding marker as rather unsuitable for this high-throughput technology.

Conclusions/Significance

Based on these data, a DNA microarray containing 64 functional oligonucleotide probes for the identification of 30 out of the 50 fish species investigated was developed. It represents the next step towards an automated and easy-to-handle method to identify fish, ichthyoplankton, and fish products.     

Pyrazolo[3,4-d]pyrimidine nucleic acids: adjustment of dA-dT to dG-dC base pair stability

Oligonucleotides incorporating 8-aza-7-deazapurin-2,6-diamine (pyrazolo[3,4-d]pyrimidin-4,6-diamine) nucleoside 2a or its 7-bromo derivative 2b show enhanced duplex stability compared to those containing dA. While incorporation of 2a opposite dT increases the T_m value only slightly, the 7-bromo compound 2b forms a very stable base pair which is as strong as the dG-dC pair. Compound 2b shows a similar base discrimination in duplex DNA as dA. The base-modified nucleosides 2a,b have a significantly more stable N-glycosylic bond than the rather labile purin-2,6-diamine 2′-deoxyribonucleoside 1. Base protection with acyl groups, with which we had difficulties in the case of purine nucleoside 1, was effective with pyrazolo[3,4-d]-pyrimidine nucleosides 2a,b. Oligonucleotides containing 2a,b were obtained by solid phase synthesis employing phosphoramidite chemistry. Compound 2b harmonizes the stability of DNA duplexes. Their stability is no longer dependent on the base pair composition while they still maintain their sequence specificity. Thus, they have the potential to reduce the number of mispairs when hybridized in solution or immobilized on arrays.     

Cloning and Overexpression of Alkaline Phosphatase PhoK from Sphingomonas sp. Strain BSAR-1 for Bioprecipitation of Uranium from Alkaline Solutions

Cells of Sphingomonas sp. strain BSAR-1 constitutively expressed an alkaline phosphatase, which was also secreted in the extracellular medium. A null mutant lacking this alkaline phosphatase activity was isolated by Tn5 random mutagenesis. The corresponding gene, designated phoK, was cloned and overexpressed in Escherichia coli strain BL21(DE3). The resultant E. coli strain EK4 overexpressed cellular activity 55 times higher and secreted extracellular PhoK activity 13 times higher than did BSAR-1. The recombinant strain very rapidly precipitated >90% of input uranium in less than 2 h from alkaline solutions (pH, 9 ± 0.2) containing 0.5 to 5 mM of uranyl carbonate, compared to BSAR-1, which precipitated uranium in >7 h. In both strains BSAR-1 and EK4, precipitated uranium remained cell bound. The EK4 cells exhibited a much higher loading capacity of 3.8 g U/g dry weight in <2 h compared to only 1.5 g U/g dry weight in >7 h in BSAR-1. The data demonstrate the potential utility of genetically engineering PhoK for the bioprecipitation of uranium from alkaline solutions.     

Selection of a reference gene for studies on lipid‐related aquatic adaptations of toothed whales (Grampus griseus)

Toothed whales are one group of marine mammals that has developed special adaptations, such as echolocation for predation, to successfully live in a dynamic aquatic environment. Their fat metabolism may differ from that of other mammals because toothed whales have acoustic fats. Gene expression in the metabolic pathways of animals can change with respect to their evolution and environment. A real‐time quantitative polymerase chain reaction (RT‐qPCR) is a reliable technique for studying the relative expressions of genes. However, since the accuracy of RT‐qPCR data is totally dependent on the reference gene, the selection of the reference gene is an essential step. In this study, 10 candidate reference genes (ZC3H10, FTL, LGALS1, RPL27, GAPDH, FTH1, DCN, TCTP, NDUS5, and UBIM) were initially tested for amplification efficiency using RT‐qPCR. After excluding DCN, the remaining nine genes, which are nearly 100% efficient, were selected for the gene stability analysis. Stable reference genes across eight different fat tissue, liver, and muscle samples from Grampus griseus were identified by four algorithms, which were provided in Genorm, NormFinder, BestKeeper, and Delta CT. Finally, a RefFinder comprehensive ranking was performed based on the stability values, and the nine genes were ranked as follows: LGALS1 > FTL > GAPDH > ZC3H10 > FTH1 > NDUS5 > TCTP > RPL27 > UBIM. The LGALS1 and FTL genes were identified as the most stable novel reference genes. The third‐ranked gene, GAPDH, is a well‐known housekeeping gene for mammals. Ultimately, we suggest the use of LGALS1 as a reliable novel reference gene for genomics studies on the lipid‐related aquatic adaptations of toothed whales.     

Estimating Finite Rate of Population Increase for Sharks Based on Vital Parameters

The vital parameter data for 62 stocks, covering 38 species, collected from the literature, including parameters of age, growth, and reproduction, were log-transformed and analyzed using multivariate analyses. Three groups were identified and empirical equations were developed for each to describe the relationships between the predicted finite rates of population increase (λ’) and the vital parameters, maximum age (T_max), age at maturity (T_m), annual fecundity (f/R_c)), size at birth (L_b), size at maturity (L_m), and asymptotic length (L_∞). Group (1) included species with slow growth rates (0.034 yr^-1 < k < 0.103 yr^-1) and extended longevity (26 yr < T_max < 81 yr), e.g., shortfin mako Isurus oxyrinchus, dusky shark Carcharhinus obscurus, etc.; Group (2) included species with fast growth rates (0.103 yr^-1 < k < 0.358 yr^-1) and short longevity (9 yr < T_max < 26 yr), e.g., starspotted smoothhound Mustelus manazo, gray smoothhound M. californicus, etc.; Group (3) included late maturing species (L_m/L_∞ ≧ 0.75) with moderate longevity (T_max < 29 yr), e.g., pelagic thresher Alopias pelagicus, sevengill shark Notorynchus cepedianus. The empirical equation for all data pooled was also developed. The λ’ values estimated by these empirical equations showed good agreement with those calculated using conventional demographic analysis. The predictability was further validated by an independent data set of three species. The empirical equations developed in this study not only reduce the uncertainties in estimation but also account for the difference in life history among groups. This method therefore provides an efficient and effective approach to the implementation of precautionary shark management measures.     

Quantitative Assessment of the Association between rs2046210 at 6q25.1 and Breast Cancer Risk

Genome-wide association studies (GWAS) have identified several genetic susceptibility loci for breast cancer (BC). One of them, conducted among Chinese women, found an association of rs2046210 at 6q25.1 with the risk of BC recently. Since then, numerous association studies have been carried out to investigate the relationship between this polymorphism and BC risk in various populations. However, these have yielded contradictory results. We therefore performed a meta-analysis to clarify this inconsistency. Overall, a total of 235003 subjects based on 13 studies were included in our study. Significantly increased BC risk was detected in the pooled analysis [allele contrast: OR = 1.13, 95%CI = 1.10–1.17, P(Z) <10⁻⁵, P(Q) <10⁻⁴; dominant model: OR = 1.21, 95%CI = 1.14–1.27, P(Z) <10⁻⁵, P(Q) <10⁻⁴; recessive model: OR = 1.18, 95%CI = 1.12–1.24, P(Z) <10⁻⁵, P(Q) = 0.04]. In addition, our data revealed that rs2046210 conferred greater risk in estrogen receptor (ER)-negative tumors [OR = 1.27, 95%CI = 1.15–1.40, P(Z) <10⁻⁵, P(Q) <10⁻⁴] than in ER-positive ones [OR = 1.18, 95%CI = 1.09–1.28, P(Z) <10⁻⁴, P(Q) = 0.0003]. When stratified by ethnicity, significant associations were found in Caucasian and Asian populations, but not detected among Africans. There was evidence of heterogeneity (P<0.05), however, the heterogeneity largely disappeared after stratification by ethnicity. The present meta-analysis demonstrated that the rs2046210 polymorphism may be associated with increased BC susceptibility, but this association varies in different ethnicities.     

Behavioral Response of Nothanguina phyllobia to Selected Plant Species

The silver-leaf nightshade nenmtode, Nothanguina phyllobia, is a promising biological control agent for its only reported host, Solanum elaeagnifolium Cav. When infective larvae of N. phyllobia and stem tissue of 39 econmnically important plant species were suspended in 0.5% water agar, nematodes aggregated about S. elaeagnifolium, Solanum carolinense L., Solanum melongena L., Solanum tuberosum L., and Prunus caroliniana (Mill.) Ait. Nematodes responded to Solanum spp. via positive chemotaxis and/or klinokinesis, but aggregated near tissue of P. caroliniana as a result of orthokinetic effects. Nematodes aggregated away from tissue of Hibiscus esculentus L., Triticum aestivum L., Santolina sp., Rosa sp., and Kochia scoparia (L.) Schrad. in the absence of orthokinetic effects. Experiments that excluded light and maintained relative humidity at 100% showed N. phyllobia to ascend the stems of 35 plant species to a height of > 9 cm within 12 h. Differences in stem ascension were not attributable to stem surface characteristics.     

EpiCollect: Linking Smartphones to Web Applications for Epidemiology,Ecology and Community Data Collection

Background

Epidemiologists and ecologists often collect data in the field and, on returning to their laboratory, enter their data into a database for further analysis. The recent introduction of mobile phones that utilise the open source Android operating system, and which include (among other features) both GPS and Google Maps, provide new opportunities for developing mobile phone applications, which in conjunction with web applications, allow two-way communication between field workers and their project databases.

Methodology

Here we describe a generic framework, consisting of mobile phone software, EpiCollect, and a web application located within www.spatialepidemiology.net. Data collected by multiple field workers can be submitted by phone, together with GPS data, to a common web database and can be displayed and analysed, along with previously collected data, using Google Maps (or Google Earth). Similarly, data from the web database can be requested and displayed on the mobile phone, again using Google Maps. Data filtering options allow the display of data submitted by the individual field workers or, for example, those data within certain values of a measured variable or a time period.

Conclusions

Data collection frameworks utilising mobile phones with data submission to and from central databases are widely applicable and can give a field worker similar display and analysis tools on their mobile phone that they would have if viewing the data in their laboratory via the web. We demonstrate their utility for epidemiological data collection and display, and briefly discuss their application in ecological and community data collection. Furthermore, such frameworks offer great potential for recruiting ‘citizen scientists’ to contribute data easily to central databases through their mobile phone.     

High Variety of Known and New RNA and DNA Viruses of Diverse Origins in Untreated Sewage

Deep sequencing of untreated sewage provides an opportunity to monitor enteric infections in large populations and for high-throughput viral discovery. A metagenomics analysis of purified viral particles in untreated sewage from the United States (San Francisco, CA), Nigeria (Maiduguri), Thailand (Bangkok), and Nepal (Kathmandu) revealed sequences related to 29 eukaryotic viral families infecting vertebrates, invertebrates, and plants (BLASTx E score, <10⁻⁴), including known pathogens (>90% protein identities) in numerous viral families infecting humans (Adenoviridae, Astroviridae, Caliciviridae, Hepeviridae, Parvoviridae, Picornaviridae, Picobirnaviridae, and Reoviridae), plants (Alphaflexiviridae, Betaflexiviridae, Partitiviridae, Sobemovirus, Secoviridae, Tombusviridae, Tymoviridae, Virgaviridae), and insects (Dicistroviridae, Nodaviridae, and Parvoviridae). The full and partial genomes of a novel kobuvirus, salivirus, and sapovirus are described. A novel astrovirus (casa astrovirus) basal to those infecting mammals and birds, potentially representing a third astrovirus genus, was partially characterized. Potential new genera and families of viruses distantly related to members of the single-stranded RNA picorna-like virus superfamily were genetically characterized and named Picalivirus, Secalivirus, Hepelivirus, Nedicistrovirus, Cadicistrovirus, and Niflavirus. Phylogenetic analysis placed these highly divergent genomes near the root of the picorna-like virus superfamily, with possible vertebrate, plant, or arthropod hosts inferred from nucleotide composition analysis. Circular DNA genomes distantly related to the plant-infecting Geminiviridae family were named Baminivirus, Nimivirus, and Niminivirus. These results highlight the utility of analyzing sewage to monitor shedding of viral pathogens and the high viral diversity found in this common pollutant and provide genetic information to facilitate future studies of these newly characterized viruses.     

Spectrophotometric Studies on Intact Muscle : II. Recovery from contractile activity

The kinetics of the mitochondrial respiratory chain of intact muscle and the concomitant changes of the intercellular pH were investigated. Addition of lactate and pyruvate under resting conditions produces reductions of DPN and cytochrome b, and, occasionally, of cytochrome c and flavoprotein. Succinate gives similar but smaller changes. In recently excised muscles moderate contractile activity produces a reduction of cytochrome c and oxidations of DPNH, cytochrome b, and sometimes of the flavoproteins. Tetanic contractions and larger numbers of twitches produce reductions of DPN and of cytochromes b and c. In sartorii of the tropical toad, stored for approximately 2 days at 0–3°C, contractile activity always gives rise to long lasting oxidations of DPNH and cytochrome b. Addition of pyruvate or lactate shortens these oxidation cycles with a concomitant reduction of cytochrome c. These responses to contractions agree with those of mitochondria isolated from leg muscles of the toad upon the addition of ADP. Apparently the mitochondria in resting, excised muscles are not supplied with an excess of substrate. Measurements on the intercellular pH showed that even limited activity ( < 5 twitches) initiates glycolysis. The primary control of respiration resides, nevertheless, in the ADP concentration, rather than in the levels of substrate or inorganic phosphate. The results are quantitatively consistent with the view that ATP is the primary energy donor for muscular contraction.     

Quantitative Trait Loci for Mercury Accumulation in Maize (Zea mays L.) Identified Using a RIL Population

To investigate the genetic mechanism of mercury accumulation in maize (Zea mays L.), a population of 194 recombinant inbred lines derived from an elite hybrid Yuyu 22, was used to identify quantitative trait loci (QTLs) for mercury accumulation at two locations. The results showed that the average Hg concentration in the different tissues of maize followed the order: leaves > bracts > stems > axis > kernels. Twenty-three QTLs for mercury accumulation in five tissues were detected on chromosomes 1, 4, 7, 8, 9 and 10, which explained 6.44% to 26.60% of the phenotype variance. The QTLs included five QTLs for Hg concentration in kernels, three QTLs for Hg concentration in the axis, six QTLs for Hg concentration in stems, four QTLs for Hg concentration in bracts and five QTLs for Hg concentration in leaves. Interestingly, three QTLs, qKHC9a, qKHC9b, and qBHC9 were in linkage with two QTLs for drought tolerance. In addition, qLHC1 was in linkage with two QTLs for arsenic accumulation. The study demonstrated the concentration of Hg in Hg-contaminated paddy soil could be reduced, and maize production maintained simultaneously by selecting and breeding maize Hg pollution-safe cultivars (PSCs).     

The Effect of Double – Blind Carbohydrate Ingestion during 60 km of Self-Paced Exercise in Warm Ambient Conditions

This study evaluated double blind ingestions of placebo (PLA) versus 6% carbohydrate (CHO) either as capsules (c) or beverage (b) during 60 km self-paced cycling in the heat (32°C and 50% relative humidity). Ten well-trained males (mean ± SD: 26±3 years; 64.5±7.7 kg and 70.7±8.8 ml.kg⁻¹.min⁻¹ maximal oxygen consumption) completed four separate 60 km time trials (TT) punctuated by 1 km sprints (14, 29, 44, 59 km) whilst ingesting either PLA_b or PLA_c or CHO_b or CHO_c. The TT was not different among treatments (PLA_b 130.2±11.2 min, CHO_b 140.5±18.1 min, PLA_c 143.1±29.2 min, CHO_c 137.3±20.1 min; P>0.05). Effect size (Cohen’s d) for time was only moderate when comparing CHO_b – PLA_b (d = 0.68) and PLA_b – PLA_c (d = 0.57) whereas all other ES were ‘trivial’ to ‘small’. Mean speed throughout the trial was significantly higher for PLA_b only (P<0.05). Power output was only different (P<0.05) between the sprints and low intensity efforts within and across conditions. Core and mean skin temperatures were similar among trials. We conclude that CHO ingestion is of little or no benefit as a beverage compared with placebo during 60 km TT in the heat.     

Whole Exome Sequencing Reveals Homozygous Mutations in RAI1, OTOF,and SLC26A4 Genes Associated with Nonsyndromic Hearing Loss in Altaian Families (South Siberia)

Hearing loss (HL) is one of the most common sensorineural disorders and several dozen genes contribute to its pathogenesis. Establishing a genetic diagnosis of HL is of great importance for clinical evaluation of deaf patients and for estimating recurrence risks for their families. Efforts to identify genes responsible for HL have been challenged by high genetic heterogeneity and different ethnic-specific prevalence of inherited deafness. Here we present the utility of whole exome sequencing (WES) for identifying candidate causal variants for previously unexplained nonsyndromic HL of seven patients from four unrelated Altaian families (the Altai Republic, South Siberia). The WES analysis revealed homozygous missense mutations in three genes associated with HL. Mutation c.2168A>G (SLC26A4) was found in one family, a novel mutation c.1111G>C (OTOF) was revealed in another family, and mutation c.5254G>A (RAI1) was found in two families. Sanger sequencing was applied for screening of identified variants in an ethnically diverse cohort of other patients with HL (n = 116) and in Altaian controls (n = 120). Identified variants were found only in patients of Altaian ethnicity (n = 93). Several lines of evidences support the association of homozygosity for discovered variants c.5254G>A (RAI1), c.1111C>G (OTOF), and c.2168A>G (SLC26A4) with HL in Altaian patients. Local prevalence of identified variants implies possible founder effect in significant number of HL cases in indigenous population of the Altai region. Notably, this is the first reported instance of patients with RAI1 missense mutation whose HL is not accompanied by specific traits typical for Smith-Magenis syndrome. Presumed association of RAI1 gene variant c.5254G>A with isolated HL needs to be proved by further experimental studies.     

Functional Expression of Drosophila para Sodium Channels : Modulation by the Membrane Protein TipE and Toxin Pharmacology

The Drosophila para sodium channel α subunit was expressed in Xenopus oocytes alone and in combination with tipE, a putative Drosophila sodium channel accessory subunit. Coexpression of tipE with para results in elevated levels of sodium currents and accelerated current decay. Para/TipE sodium channels have biophysical and pharmacological properties similar to those of native channels. However, the pharmacology of these channels differs from that of vertebrate sodium channels: (a) toxin II from Anemonia sulcata, which slows inactivation, binds to Para and some mammalian sodium channels with similar affinity (K _d ≅ 10 nM), but this toxin causes a 100-fold greater decrease in the rate of inactivation of Para/TipE than of mammalian channels; (b) Para sodium channels are >10-fold more sensitive to block by tetrodotoxin; and (c) modification by the pyrethroid insecticide permethrin is >100-fold more potent for Para than for rat brain type IIA sodium channels. Our results suggest that the selective toxicity of pyrethroid insecticides is due at least in part to the greater affinity of pyrethroids for insect sodium channels than for mammalian sodium channels.