Genetic Variation among Major Human Geographic Groups Supports a Peculiar Evolutionary Trend in PAX9

A total of 172 persons from nine South Amerindian, three African and one Eskimo populations were studied in relation to the Paired box gene 9 (PAX9) exon 3 (138 base pairs) as well as its 5' and 3' flanking intronic segments (232 bp and 220 bp, respectively) and integrated with the information available for the same genetic region from individuals of different geographical origins. Nine mutations were scored in exon 3 and six in its flanking regions; four of them are new South American tribe-specific singletons. Exon3 nucleotide diversity is several orders of magnitude higher than its intronic regions. Additionally, a set of variants in the PAX9 and 101 other genes related with dentition can define at least some dental morphological differences between Sub-Saharan Africans and non-Africans, probably associated with adaptations after the modern human exodus from Africa. Exon 3 of PAX9 could be a good molecular example of how evolvability works.

Testing Evolutionary and Dispersion Scenarios for the Settlement of the New World

Background: Discussion surrounding the settlement of the New World has recently gained momentum with advances in molecular biology, archaeology and bioanthropology. Recent evidence from these diverse fields is found to support different colonization scenarios. The currently available genetic evidence suggests a ""single migration'' model, in which both early and later Native American groups derive from one expansion event into the continent. In contrast, the pronounced anatomical differences between early and late Native American populations have led others to propose more complex scenarios, involving separate colonization events of the New World and a distinct origin for these groups. Methodology/Principal Findings: Using large samples of Early American crania, we: 1) calculated the rate of morphological differentiation between Early and Late American samples under three different time divergence assumptions, and compared our findings to the predicted morphological differentiation under neutral conditions in each case; and 2) further tested three dispersal scenarios for the colonization of the New World by comparing the morphological distances among early and late Amerindians, East Asians, Australo-Melanesians and early modern humans from Asia to geographical distances associated with each dispersion model. Results indicate that the assumption of a last shared common ancestor outside the continent better explains the observed morphological differences between early and late American groups. This result is corroborated by our finding that a model comprising two Asian waves of migration coming through Bering into the Americas fits the cranial anatomical evidence best, especially when the effects of diversifying selection to climate are taken into account. Conclusions: We conclude that the morphological diversity documented through time in the New World is best accounted for by a model postulating two waves of human expansion into the continent originating in East Asia and entering through Beringia.

Genetic variability and natural selection at the ligand domain of the Duffy binding protein in brazilian Plasmodium vivax populations

Background: Plasmodium vivax malaria is a major public health challenge in Latin America, Asia and Oceania, with 130-435 million clinical cases per year worldwide. Invasion of host blood cells by P. vivax mainly depends on a type I membrane protein called Duffy binding protein (PvDBP). The erythrocyte-binding motif of PvDBP is a 170 amino-acid stretch located in its cysteine-rich region II (PvDBP(II)), which is the most variable segment of the protein. Methods: To test whether diversifying natural selection has shaped the nucleotide diversity of PvDBP(II) in Brazilian populations, this region was sequenced in 122 isolates from six different geographic areas. A Bayesian method was applied to test for the action of natural selection under a population genetic model that incorporates recombination. The analysis was integrated with a structural model of PvDBP(II), and T-and B-cell epitopes were localized on the 3-D structure. Results: The results suggest that: (i) recombination plays an important role in determining the haplotype structure of PvDBP(II), and (ii) PvDBP(II) appears to contain neutrally evolving codons as well as codons evolving under natural selection. Diversifying selection preferentially acts on sites identified as epitopes, particularly on amino acid residues 417, 419, and 424, which show strong linkage disequilibrium. Conclusions: This study shows that some polymorphisms of PvDBP(II) are present near the erythrocyte-binding domain and might serve to elude antibodies that inhibit cell invasion. Therefore, these polymorphisms should be taken into account when designing vaccines aimed at eliciting antibodies to inhibit erythrocyte invasion.

Particle competition for complex network community detection

In many real situations, randomness is considered to be uncertainty or even confusion which impedes human beings from making a correct decision. Here we study the combined role of randomness and determinism in particle dynamics for complex network community detection. In the proposed model, particles walk in the network and compete with each other in such a way that each of them tries to possess as many nodes as possible. Moreover, we introduce a rule to adjust the level of randomness of particle walking in the network, and we have found that a portion of randomness can largely improve the community detection rate. Computer simulations show that the model has good community detection performance and at the same time presents low computational complexity. (C) 2008 American Institute of Physics.

Generalized connectivity between any two nodes in a complex network

This article focuses on the identification of the number of paths with different lengths between pairs of nodes in complex networks and how these paths can be used for characterization of topological properties of theoretical and real-world complex networks. This analysis revealed that the number of paths can provide a better discrimination of network models than traditional network measurements. In addition, the analysis of real-world networks suggests that the long-range connectivity tends to be limited in these networks and may be strongly related to network growth and organization.

Finding Quasi-Optimal Network Topologies for Information Transmission in Active Networks

This work clarifies the relation between network circuit (topology) and behaviour (information transmission and synchronization) in active networks, e.g. neural networks. As an application, we show how one can find network topologies that are able to transmit a large amount of information, possess a large number of communication channels, and are robust under large variations of the network coupling configuration. This theoretical approach is general and does not depend on the particular dynamic of the elements forming the network, since the network topology can be determined by finding a Laplacian matrix (the matrix that describes the connections and the coupling strengths among the elements) whose eigenvalues satisfy some special conditions. To illustrate our ideas and theoretical approaches, we use neural networks of electrically connected chaotic Hindmarsh-Rose neurons.

Network of recurrent events for the Olami-Feder-Christensen model

We numerically study the dynamics of a discrete spring-block model introduced by Olami, Feder, and Christensen (OFC) to mimic earthquakes and investigate to what extent this simple model is able to reproduce the observed spatiotemporal clustering of seismicity. Following a recently proposed method to characterize such clustering by networks of recurrent events [J. Davidsen, P. Grassberger, and M. Paczuski, Geophys. Res. Lett. 33, L11304 (2006)], we find that for synthetic catalogs generated by the OFC model these networks have many nontrivial statistical properties. This includes characteristic degree distributions, very similar to what has been observed for real seismicity. There are, however, also significant differences between the OFC model and earthquake catalogs, indicating that this simple model is insufficient to account for certain aspects of the spatiotemporal clustering of seismicity.

Automatic Network Fingerprinting through Single-Node Motifs

Complex networks have been characterised by their specific connectivity patterns (network motifs), but their building blocks can also be identified and described by node-motifs-a combination of local network features. One technique to identify single node-motifs has been presented by Costa et al. (L. D. F. Costa, F. A. Rodrigues, C. C. Hilgetag, and M. Kaiser, Europhys. Lett., 87, 1, 2009). Here, we first suggest improvements to the method including how its parameters can be determined automatically. Such automatic routines make high-throughput studies of many networks feasible. Second, the new routines are validated in different network-series. Third, we provide an example of how the method can be used to analyse network time-series. In conclusion, we provide a robust method for systematically discovering and classifying characteristic nodes of a network. In contrast to classical motif analysis, our approach can identify individual components (here: nodes) that are specific to a network. Such special nodes, as hubs before, might be found to play critical roles in real-world networks.

Gene Expression Complex Networks: Synthesis, Identification, and Analysis

Thanks to recent advances in molecular biology, allied to an ever increasing amount of experimental data, the functional state of thousands of genes can now be extracted simultaneously by using methods such as cDNA microarrays and RNA-Seq. Particularly important related investigations are the modeling and identification of gene regulatory networks from expression data sets. Such a knowledge is fundamental for many applications, such as disease treatment, therapeutic intervention strategies and drugs design, as well as for planning high-throughput new experiments. Methods have been developed for gene networks modeling and identification from expression profiles. However, an important open problem regards how to validate such approaches and its results. This work presents an objective approach for validation of gene network modeling and identification which comprises the following three main aspects: (1) Artificial Gene Networks (AGNs) model generation through theoretical models of complex networks, which is used to simulate temporal expression data; (2) a computational method for gene network identification from the simulated data, which is founded on a feature selection approach where a target gene is fixed and the expression profile is observed for all other genes in order to identify a relevant subset of predictors; and (3) validation of the identified AGN-based network through comparison with the original network. The proposed framework allows several types of AGNs to be generated and used in order to simulate temporal expression data. The results of the network identification method can then be compared to the original network in order to estimate its properties and accuracy. Some of the most important theoretical models of complex networks have been assessed: the uniformly-random Erdos-Renyi (ER), the small-world Watts-Strogatz (WS), the scale-free Barabasi-Albert (BA), and geographical networks (GG). The experimental results indicate that the inference method was sensitive to average degree k variation, decreasing its network recovery rate with the increase of k. The signal size was important for the inference method to get better accuracy in the network identification rate, presenting very good results with small expression profiles. However, the adopted inference method was not sensible to recognize distinct structures of interaction among genes, presenting a similar behavior when applied to different network topologies. In summary, the proposed framework, though simple, was adequate for the validation of the inferred networks by identifying some properties of the evaluated method, which can be extended to other inference methods.

Feature selection environment for genomic applications

Background: Feature selection is a pattern recognition approach to choose important variables according to some criteria in order to distinguish or explain certain phenomena (i.e., for dimensionality reduction). There are many genomic and proteomic applications that rely on feature selection to answer questions such as selecting signature genes which are informative about some biological state, e. g., normal tissues and several types of cancer; or inferring a prediction network among elements such as genes, proteins and external stimuli. In these applications, a recurrent problem is the lack of samples to perform an adequate estimate of the joint probabilities between element states. A myriad of feature selection algorithms and criterion functions have been proposed, although it is difficult to point the best solution for each application. Results: The intent of this work is to provide an open-source multiplataform graphical environment for bioinformatics problems, which supports many feature selection algorithms, criterion functions and graphic visualization tools such as scatterplots, parallel coordinates and graphs. A feature selection approach for growing genetic networks from seed genes ( targets or predictors) is also implemented in the system. Conclusion: The proposed feature selection environment allows data analysis using several algorithms, criterion functions and graphic visualization tools. Our experiments have shown the software effectiveness in two distinct types of biological problems. Besides, the environment can be used in different pattern recognition applications, although the main concern regards bioinformatics tasks.

Identification of Novel Immunoregulatory Molecules in Human Thymic Regulatory CD4(+)CD25(+) T Cells by Phage Display

Thymic CD4(+)CD25(+) cells play an important role in immune regulation and are continuously developed in the thymus as an independent lineage. How these cells are generated, what are their multiple pathways of suppressive activity and which are their specific markers are questions that remain unanswered. To identify molecules involved in the function and development of human CD4(+)CD25(+) T regulatory cells we targeted thymic CD4(+)CD25(+) cells by peptide phage display. A phage library containing random peptides was screened ex vivo for binding to human thymic CD4(+)CD25(+) T cells. After four rounds of selection on CD4(+)CD25(+) enriched populations of thymocytes, we sequenced several phage displayed peptides and selected one with identity to the Vitamin D Receptor (VDR). We confirmed the binding of the VDR phage to active Vitamin D in vitro, as well as the higher expression of VDR in CD4(+)CD25(+) cells. We suggest that differential expression of VDR on natural Tregs may be related to the relevance of Vitamin D in function and ontogeny of these cells.

Successful maintenance of a stingless bee population despite a severe genetic bottleneck

Stingless bees play an important ecological role as pollinators of many wild plant species in the tropics and have significant potential for the pollination of agricultural crops. Nevertheless, conservation efforts as well as commercial breeding programmes require better guidelines on the amount of genetic variation that is needed to maintain viable populations. In this context, we carried out a long-term genetic study on the stingless bee Melipona scutellaris to evaluate the population viability consequences of prolonged breeding from a small number of founder colonies. In particular, it was artificially imposed a genetic bottleneck by setting up a population starting from only two founder colonies, and continued breeding from it for a period of over 10 years in a location outside its natural area of occurrence. We show that despite a great reduction in the number of alleles present at both neutral microsatellite loci and the sex-determining locus relative to its natural source population, and an increased frequency in the production of sterile diploid males, the genetically impoverished population could be successfully bred and maintained for at least 10 years. This shows that in stingless bees, breeding from a small stock of colonies may have less severe consequences than previously suspected. In addition, we provide a simulation model to determine the number of colonies that are needed to maintain a certain number of sex alleles in a population, thereby providing useful guidelines for stingless bee breeding and conservation efforts.

Angiotensin receptor blockade improves the net balance of cardiac Ca(2+) handling-related proteins in sympathetic hyperactivity-induced heart failure

Aims: The clinical benefits of angiotensin II type 1 (AT1) receptor blockers (ARB) in heart failure (HF) include cardiac anti-remodeling and improved ventricular function. However, the cellular mechanisms underlying the benefits of ARB on ventricular function need to be better clarified. In the present manuscript, we evaluated the effects of AT1 receptor blockade on the net balance of Ca(2+) handling proteins in hearts of mice lacking alpha(2A) and alpha(2C) adrenoceptors (alpha(2A)/alpha(2C)ARKO), which develop sympathetic hyperactivity (SH) induced-HF. Main methods: A cohort of male wild-type (WT) and congenic alpha(2A)/alpha(2C)ARKO mice in a C57BL6/J genetic background (5-7 mo of age) was randomly assigned to receive either placebo or ARB (Losartan, 10 mg/kg for 8wks). Ventricular function (VF) was assessed by echocardiography, and cardiac myocyte width and ventricular fibrosis by a computer-assisted morphometric system. Sarcoplasmic reticulum Ca(2+) ATPase (SERCA2), phospholamban (PLN), phospho-Ser(16)-PLN, phospho-Thr(17)-PLN, phosphatase 1 (PP1), Na(+)-Ca(2+) exchanger (NCX), Ca(2+)/calmodulin-dependent protein kinase 11 (CaMKII) and phospho-Thr(286)-CaMKII were analyzed by Western blot. Key findings: alpha(2A)/alpha(2C)ARKO mice displayed ventricular dysfunction, cardiomyocyte hypertrophy and cardiac fibrosis paralleled by decreased SERCA2 and increased phospho-Thr(17)-PLN, CaMKII, phospho-Thr(286)-CaMKII and NCX levels. ARB induced anti-cardiac remodeling effect and improved VF in alpha(2A)/alpha(2C)ARKO associated with increased SERCA2 and phospho-Ser(16)-PLN levels, and SERCA2:NCX ratio. Additionally, ARB decreased phospho-Thr(17)-PLN levels as well as reestablished NCX, CaMKII and phospho-Thr(286)-CaMKII toward WT levels. Significance: Altogether, these data provide new insights on intracellular Ca(2+) regulatory mechanisms underlying improved ventricular function by ARB therapy in HF. (c) 2011 Elsevier Inc. All rights reserved.

Exercise training reduces cardiac angiotensin II levels and prevents cardiac dysfunction in a genetic model of sympathetic hyperactivity-induced heart failure in mice

The role of exercise training (ET) on cardiac renin-angiotensin system (RAS) was investigated in 3-5 month-old mice lacking alpha(2A-) and alpha(2C-)adrenoceptors (alpha(2A)/alpha(2C)ARKO) that present heart failure (HF) and wild type control (WT). ET consisted of 8-week running sessions of 60 min, 5 days/week. In addition, exercise tolerance, cardiac structural and function analysis were made. At 3 months, fractional shortening and exercise tolerance were similar between groups. At 5 months, alpha(2A)/alpha(2C)ARKO mice displayed ventricular dysfunction and fibrosis associated with increased cardiac angiotensin (Ang) II levels (2.9-fold) and increased local angiotensin-converting enzyme activity (ACE 18%). ET decreased alpha(2A)/alpha(2C)ARKO cardiac Ang II levels and ACE activity to age-matched untrained WT mice levels while increased ACE2 expression and prevented exercise intolerance and ventricular dysfunction with little impact on cardiac remodeling. Altogether, these data provide evidence that reduced cardiac RAS explains, at least in part, the beneficial effects of ET on cardiac function in a genetic model of HF.

Intracellular mechanisms of specific beta-adrenoceptor antagonists involved in improved cardiac function and survival in a genetic model of heart failure

beta-blockers, as class, improve cardiac function and survival in heart failure (HF). However, the molecular mechanisms underlying these beneficial effects remain elusive. In the present study, metoprolol and carvedilol were used in doses that display comparable heart rate reduction to assess their beneficial effects in a genetic model of sympathetic hyperactivity-induced HF (alpha(2A)/alpha(2C)-ARKO mice). Five month-old HF mice were randomly assigned to receive either saline, metoprolol or carvedilol for 8 weeks and age-matched wild-type mice (WT) were used as controls. HF mice displayed baseline tachycardia, systolic dysfunction evaluated by echocardiography, 50% mortality rate, increased cardiac myocyte width (50%) and ventricular fibrosis (3-fold) compared with WT. All these responses were significantly improved by both treatments. Cardiomyocytes from HF mice showed reduced peak [Ca(2+)](i) transient (13%) using confocal microscopy imaging. Interestingly, while metoprolol improved [Ca(2+)](i) transient, carvedilol had no effect on peak [Ca(2+)](i) transient but also increased [Ca(2+)] transient decay dynamics. We then examined the influence of carvedilol in cardiac oxidative stress as an alternative target to explain its beneficial effects. Indeed, HF mice showed 10-fold decrease in cardiac reduced/oxidized glutathione ratio compared with WT, which was significantly improved only by carvedilol treatment. Taken together, we provide direct evidence that the beneficial effects of metoprolol were mainly associated with improved cardiac Ca(2+) transients and the net balance of cardiac Ca(2+) handling proteins while carvedilol preferentially improved cardiac redox state. (C) 2008 Elsevier Inc. All rights reserved.