Uncovering the overlapping community structure of complex networks in nature and society

Many complex systems in nature and society can be described in terms of networks capturing the intricate web of connections among the units they are made of. A key question is how to interpret the global organization of such networks as the coexistence of their structural subunits (communities) associated with more highly interconnected parts. Identifying these a priori unknown building blocks (such as functionally related proteins, industrial sectors and groups of people) is crucial to the understanding of the structural and functional properties of networks. The existing deterministic methods used for large networks find separated communities, whereas most of the actual networks are made of highly overlapping cohesive groups of nodes. Here we introduce an approach to analysing the main statistical features of the interwoven sets of overlapping communities that makes a step towards uncovering the modular structure of complex systems. After defining a set of new characteristic quantities for the statistics of communities, we apply an efficient technique for exploring overlapping communities on a large scale. We find that overlaps are significant, and the distributions we introduce reveal universal features of networks. Our studies of collaboration, word-association and protein interaction graphs show that the web of communities has non-trivial correlations and specific scaling properties.Comment: The free academic research software, CFinder, used for the publication is available at the website of the publication: http://angel.elte.hu/clusterin

Separation of Anti-Proliferation and Anti-Apoptotic Functions of Retinoblastoma Protein through Targeted Mutations of Its A/B Domain

BACKGROUND: The human retinoblastoma susceptibility gene encodes a nuclear phosphoprotein RB, which is a negative regulator of cell proliferation. The growth suppression function of RB requires an evolutionarily conserved A/B domain that contains two distinct peptide-binding pockets. At the A/B interface is a binding site for the C-terminal trans-activation domain of E2F. Within the B-domain is a binding site for proteins containing the LxCxE peptide motif. METHODOLOGY/PRINCIPLE FINDINGS: Based on the crystal structure of the A/B domain, we have constructed an RB-K530A/N757F (KN) mutant to disrupt the E2F- and LxCxE-binding pockets. The RB-K530A (K) mutant is sufficient to inactivate the E2F-binding pocket, whereas the RB-N757F (N) mutant is sufficient to inactivate the LxCxE-binding pocket. Each single mutant inhibits cell proliferation, but the RB-KN double mutant is defective in growth suppression. Nevertheless, the RB-KN mutant is capable of reducing etoposide-induced apoptosis. CONCLUSION/SIGNIFICANCE: Previous studies have established that RB-dependent G1-arrest can confer resistance to DNA damage-induced apoptosis. Results from this study demonstrate that RB can also inhibit apoptosis independent of growth suppression

RNA secondary structure prediction from multi-aligned sequences

It has been well accepted that the RNA secondary structures of most functional non-coding RNAs (ncRNAs) are closely related to their functions and are conserved during evolution. Hence, prediction of conserved secondary structures from evolutionarily related sequences is one important task in RNA bioinformatics; the methods are useful not only to further functional analyses of ncRNAs but also to improve the accuracy of secondary structure predictions and to find novel functional RNAs from the genome. In this review, I focus on common secondary structure prediction from a given aligned RNA sequence, in which one secondary structure whose length is equal to that of the input alignment is predicted. I systematically review and classify existing tools and algorithms for the problem, by utilizing the information employed in the tools and by adopting a unified viewpoint based on maximum expected gain (MEG) estimators. I believe that this classification will allow a deeper understanding of each tool and provide users with useful information for selecting tools for common secondary structure predictions.Comment: A preprint of an invited review manuscript that will be published in a chapter of the book `Methods in Molecular Biology'. Note that this version of the manuscript may differ from the published versio

Genetic and environmental influence on thyroid gland volume and thickness of thyroid isthmus: a twin study.

Objectives Decreased thyroid volume has been related to increased prevalence of thyroid cancer.Subjects and methods One hundred and fourteen Hungarian adult twin pairs (69 monozygotic, 45 dizygotic) with or without known thyroid disorders underwent thyroid ultrasound. Thickness of the thyroid isthmus was measured at the thickest portion of the gland in the midline using electronic calipers at the time of scanning. Volume of the thyroid lobe was computed according to the following formula: thyroid height*width*depth*correction factor (0.63).Results Age-, sex-, body mass index- and smoking-adjusted heritability of the thickness of thyroid isthmus was 50% (95% confidence interval [CI], 35 to 66%). Neither left nor right thyroid volume showed additive genetic effects, but shared environments were 68% (95% CI, 48 to 80%) and 79% (95% CI, 72 to 87%), respectively. Magnitudes of monozygotic and dizygotic co-twin correlations were not substantially impacted by the correction of covariates of body mass index and smoking. Unshared environmental effects showed a moderate influence on dependent parameters (24-50%).Conclusions Our analysis support that familial factors are important for thyroid measures in a general twin population. A larger sample size is needed to show whether this is because of common environmental (e.g. intrauterine effects, regional nutrition habits, iodine supply) or genetic effects

Skeeter Buster: A Stochastic, Spatially Explicit Modeling Tool for Studying Aedes aegypti Population Replacement and Population Suppression Strategies

Dengue is a viral disease that affects approximately 50 million people annually, and is estimated to result in 12,500 fatalities. Dengue viruses are vectored by mosquitoes, predominantly by the species Aedes aegypti. Because there is currently no vaccine or specific treatment, the only available strategy to reduce dengue transmission is to control the populations of these mosquitoes. This can be achieved by traditional approaches such as insecticides, or by recently developed genetic methods that propose the release of mosquitoes genetically engineered to be unable to transmit dengue viruses. The expected outcome of different control strategies can be compared by simulating the population dynamics and genetics of mosquitoes at a given location. Development of optimal control strategies can then be guided by the modeling approach. To that end, we introduce a new modeling tool called Skeeter Buster. This model describes the dynamics and the genetics of Ae. aegypti populations at a very fine scale, simulating the contents of individual houses, and even the individual water-holding containers in which mosquito larvae reside. Skeeter Buster can be used to compare the predicted outcomes of multiple control strategies, traditional or genetic, making it an important tool in the fight against dengue

Variations in the Difference between Mean Sea Level measured either side of Cape Hatteras and Their Relation to the North Atlantic Oscillation

We consider the extent to which the difference in mean sea level (MSL) measured on the North American Atlantic coast either side of Cape Hatteras varies as a consequence of dynamical changes in the ocean caused by fluctuations in the North Atlantic Oscillation (NAO). From analysis of tide gauge data, we know that changes in MSL-difference and NAO index are correlated on decadal to century timescales enabling a scale factor of MSL-difference change per unit change in NAO index to be estimated. Changes in trend in the NAO index have been small during the past few centuries (when measured using windows of order 60–120 years). Therefore, if the same scale factor applies through this period of time, the corresponding changes in trend in MSL-difference for the past few centuries should also have been small. It is suggested thereby that the sea level records for recent centuries obtained from salt marshes (adjusted for long-term vertical land movements) should have essentially the same NAO-driven trends south and north of Cape Hatteras, only differing due to contributions from other processes such as changes in the Meridional Overturning Circulation or ‘geophysical fingerprints’. The salt marsh data evidently support this interpretation within their uncertainties for the past few centuries, and perhaps even for the past millennium. Recommendations are made on how greater insight might be obtained by acquiring more measurements and by improved modelling of the sea level response to wind along the shelf

Novel Cβ–Cγ Bond Cleavages of Tryptophan-Containing Peptide Radical Cations

In this study, we observed unprecedented cleavages of the Cβ–Cγ bonds of tryptophan residue side chains in a series of hydrogen-deficient tryptophan-containing peptide radical cations (M•+) during low-energy collision-induced dissociation (CID). We used CID experiments and theoretical density functional theory (DFT) calculations to study the mechanism of this bond cleavage, which forms [M – 116]+ ions. The formation of an α-carbon radical intermediate at the tryptophan residue for the subsequent Cβ–Cγ bond cleavage is analogous to that occurring at leucine residues, producing the same product ions; this hypothesis was supported by the identical product ion spectra of [LGGGH – 43]+ and [WGGGH – 116]+, obtained from the CID of [LGGGH]•+ and [WGGGH]•+, respectively. Elimination of the neutral 116-Da radical requires inevitable dehydrogenation of the indole nitrogen atom, leaving the radical centered formally on the indole nitrogen atom ([Ind]•-2), in agreement with the CID data for [WGGGH]•+ and [W1-CH3GGGH]•+; replacing the tryptophan residue with a 1-methyltryptophan residue results in a change of the base peak from that arising from a neutral radical loss (116 Da) to that arising from a molecule loss (131 Da), both originating from Cβ–Cγ bond cleavage. Hydrogen atom transfer or proton transfer to the γ-carbon atom of the tryptophan residue weakens the Cβ–Cγ bond and, therefore, decreases the dissociation energy barrier dramatically

Single-nucleotide polymorphisms in the RB1 gene and association with breast cancer in the British population

A substantial proportion of the familial risk of breast cancer may be attributable to genetic variants each contributing a small effect. pRb controls the cell cycle and polymorphisms within it are candidates for such low penetrance susceptibility alleles, since the gene has been implicated in several human tumours, particularly breast cancer. The purpose of this study was to determine whether common variants in the RB1 gene are associated with breast cancer risk. We assessed 15 tagging single-nucleotide polymorphisms (SNPs) using a case–control study design (n⩽4474 cases and n⩽4560 controls). A difference in genotype frequencies was found between cases and controls for rs2854344 in intron 17 (P-trend=0.007) and rs198580 in intron 19 (P-trend=0.018). Carrying the minor allele of these SNPs appears to confer a protective effect on breast cancer risk (odd ratio (OR)=0.86 (0.76–0.96) for rs2854344 and OR=0.80 (0.66–0.96) for rs198580). However, after adjusting for multiple testing these associations were borderline with an adjusted P-trend=0.068 for the most significant SNP (rs2854344). The RB1 gene is not known to contain any coding SNPs with allele frequencies ⩾5% but several intronic variants are in perfect linkage disequilibrium with the associated SNPs. Replication studies are needed to confirm the associations with breast cancer