Archaeological Support for the Three-Stage Expansion of Modern Humans across Northeastern Eurasia and into the Americas

Background

Understanding the dynamics of the human range expansion across northeastern Eurasia during the late Pleistocene is central to establishing empirical temporal constraints on the colonization of the Americas [1]. Opinions vary widely on how and when the Americas were colonized, with advocates supporting either a pre-[2] or post-[1], [3], [4], [5], [6] last glacial maximum (LGM) colonization, via either a land bridge across Beringia [3], [4], [5], a sea-faring Pacific Rim coastal route [1], [3], a trans-Arctic route [4], or a trans-Atlantic oceanic route [5]. Here we analyze a large sample of radiocarbon dates from the northeast Eurasian Upper Paleolithic to identify the origin of this expansion, and estimate the velocity of colonization wave as it moved across northern Eurasia and into the Americas.

Methodology/Principal Findings

We use diffusion models [6], [7] to quantify these dynamics. Our results show the expansion originated in the Altai region of southern Siberia ∼46kBP , and from there expanded across northern Eurasia at an average velocity of 0.16 km per year. However, the movement of the colonizing wave was not continuous but underwent three distinct phases: 1) an initial expansion from 47-32k calBP; 2) a hiatus from ∼32-16k calBP, and 3) a second expansion after the LGM ∼16k calBP. These results provide archaeological support for the recently proposed three-stage model of the colonization of the Americas [8], [9]. Our results falsify the hypothesis of a pre-LGM terrestrial colonization of the Americas and we discuss the importance of these empirical results in the light of alternative models.

Conclusions/Significance

Our results demonstrate that the radiocarbon record of Upper Paleolithic northeastern Eurasia supports a post-LGM terrestrial colonization of the Americas falsifying the proposed pre-LGM terrestrial colonization of the Americas. We show that this expansion was not a simple process, but proceeded in three phases, consistent with genetic data, largely in response to the variable climatic conditions of late Pleistocene northeast Eurasia. Further, the constraints imposed by the spatiotemporal gradient in the empirical radiocarbon record across this entire region suggests that North America cannot have been colonized much before the existing Clovis radiocarbon record suggests.     

Dopamine-Induced Conformational Changes in Alpha-Synuclein

Background

Oligomerization and aggregation of α-synuclein molecules play a major role in neuronal dysfunction and loss in Parkinson''s disease [1]. However, α-synuclein oligomerization and aggregation have mostly been detected indirectly in cells using detergent extraction methods [2], [3], [4]. A number of in vitro studies showed that dopamine can modulate the aggregation of α-synuclein by inhibiting the formation of or by disaggregating amyloid fibrils [5], [6], [7].

Methodology/Principal Findings

Here, we show that α-synuclein adopts a variety of conformations in primary neuronal cultures using fluorescence lifetime imaging microscopy (FLIM). Importantly, we found that dopamine, but not dopamine agonists, induced conformational changes in α-synuclein which could be prevented by blocking dopamine transport into the cell. Dopamine also induced conformational changes in α-synuclein expressed in neuronal cell lines, and these changes were also associated with alterations in oligomeric/aggregated species.

Conclusion/Significance

Our results show, for the first time, a direct effect of dopamine on the conformation of α-synuclein in neurons, which may help explain the increased vulnerability of dopaminergic neurons in Parkinson''s disease.     

Multi-Locus Assortment (MLA) for Transgene Dispersal and Elimination in Mosquito Populations

Background

Replacement of wild-type mosquito populations with genetically modified versions is being explored as a potential strategy to control vector-borne diseases. Due to lower expected relative fitness of transgenic individuals, transgenes must be driven into populations for these scenarios to be successful. Several gene drive mechanisms exist in a theoretical sense but none are currently workable in mosquitoes. Even if strategies were workable, it would be very difficult to recall released transgenes in the event of unforeseen consequences. What is needed is a way to test transgenes in the field for feasibility, efficacy and safety prior to releasing an active drive mechanism.

Methodology/Principal Findings

We outline a method, termed Multi-locus assortment (MLA), to spread transgenes into vector populations by the release of genetically-modified mosquitoes carrying multiple stable transgene inserts. Simulations indicate that [1] insects do not have to carry transgenes at more than 4 loci, [2] transgenes can be maintained at high levels by sequential small releases, the frequency of which depends on the construct fitness cost, and [3] in the case of unforeseen negative non-target effects, transgenes can be eliminated from the population by halting transgenic releases and/or mass releases of wild-type insects. We also discuss potential methods to create MLA mosquito strains in the laboratory.

Conclusions/Significance

While not as efficient as active drive mechanisms, MLA has other advantages: [1] MLA strains can be constructed for some mosquito species with currently-available technology, [2] MLA will allow the ecological components of transgenic mosquito releases to be tested before actual gene drive mechanisms are ready to be deployed, [3] since MLA is not self-propagating, the risk of an accidental premature release into nature is minimized, and [4] in the case that active gene drive mechanisms prove impossible to develop, the MLA approach can be used as a back-up transgene dispersal mechanism for disease control efforts in some systems.     

A Universal Model of Commuting Networks

We show that a recently proposed model generates accurate commuting networks on 80 case studies from different regions of the world (Europe and United-States) at different scales (e.g. municipalities, counties, regions). The model takes as input the number of commuters coming in and out of each geographic unit and generates the matrix of commuting flows between the units. The single parameter of the model follows a universal law that depends only on the scale of the geographic units. We show that our model significantly outperforms two other approaches proposing a universal commuting model [1], [2], particularly when the geographic units are small (e.g. municipalities).     

Inhibition of mTOR by Rapamycin Abolishes Cognitive Deficits and Reduces Amyloid-β Levels in a Mouse Model of Alzheimer's Disease

Background

Reduced TOR signaling has been shown to significantly increase lifespan in a variety of organisms [1], [2], [3], [4]. It was recently demonstrated that long-term treatment with rapamycin, an inhibitor of the mTOR pathway[5], or ablation of the mTOR target p70S6K[6] extends lifespan in mice, possibly by delaying aging. Whether inhibition of the mTOR pathway would delay or prevent age-associated disease such as AD remained to be determined.

Methodology/Principal Findings

We used rapamycin administration and behavioral tools in a mouse model of AD as well as standard biochemical and immunohistochemical measures in brain tissue to provide answers for this question. Here we show that long-term inhibition of mTOR by rapamycin prevented AD-like cognitive deficits and lowered levels of Aβ₄₂, a major toxic species in AD[7], in the PDAPP transgenic mouse model. These data indicate that inhibition of the mTOR pathway can reduce Aβ₄₂ levels in vivo and block or delay AD in mice. As expected from the inhibition of mTOR, autophagy was increased in neurons of rapamycin-treated transgenic, but not in non-transgenic, PDAPP mice, suggesting that the reduction in Aβ and the improvement in cognitive function are due in part to increased autophagy, possibly as a response to high levels of Aβ.

Conclusions/Significance

Our data suggest that inhibition of mTOR by rapamycin, an intervention that extends lifespan in mice, can slow or block AD progression in a transgenic mouse model of the disease. Rapamycin, already used in clinical settings, may be a potentially effective therapeutic agent for the treatment of AD.     

Phenylbutyric Acid Rescues Endoplasmic Reticulum Stress-Induced Suppression of APP Proteolysis and Prevents Apoptosis in Neuronal Cells

Background

The familial and sporadic forms of Alzheimer''s disease (AD) have an identical pathology with a severe disparity in the time of onset [1]. The pathological similarity suggests that epigenetic processes may phenocopy the Familial Alzheimer''s disease (FAD) mutations within sporadic AD. Numerous groups have demonstrated that FAD mutations in presenilin result in ‘loss of function’ of γ-secretase mediated APP cleavage [2], [3], [4], [5]. Accordingly, ER stress is prominent within the pathologically impacted brain regions in AD patients [6] and is reported to inhibit APP trafficking through the secretory pathway [7], [8]. As the maturation of APP and the cleaving secretases requires trafficking through the secretory pathway [9], [10], [11], we hypothesized that ER stress may block trafficking requisite for normal levels of APP cleavage and that the small molecular chaperone 4-phenylbutyrate (PBA) may rescue the proteolytic deficit.

Methodology/Principal Findings

The APP-Gal4VP16/Gal4-reporter screen was stably incorporated into neuroblastoma cells in order to assay γ-secretase mediated APP proteolysis under normal and pharmacologically induced ER stress conditions. Three unrelated pharmacological agents (tunicamycin, thapsigargin and brefeldin A) all repressed APP proteolysis in parallel with activation of unfolded protein response (UPR) signaling—a biochemical marker of ER stress. Co-treatment of the γ-secretase reporter cells with PBA blocked the repressive effects of tunicamycin and thapsigargin upon APP proteolysis, UPR activation, and apoptosis. In unstressed cells, PBA stimulated γ-secretase mediated cleavage of APP by 8–10 fold, in the absence of any significant effects upon amyloid production, by promoting APP trafficking through the secretory pathway and the stimulation of the non-pathogenic α/γ-cleavage.

Conclusions/Significance

ER stress represses γ-secretase mediated APP proteolysis, which replicates some of the proteolytic deficits associated with the FAD mutations. The small molecular chaperone PBA can reverse ER stress induced effects upon APP proteolysis, trafficking and cellular viability. Pharmaceutical agents, such as PBA, that stimulate α/γ-cleavage of APP by modifying intracellular trafficking should be explored as AD therapeutics.     

Functional Analysis of Prognostic Gene Expression Network Genes in Metastatic Breast Cancer Models

Identification of conserved co-expression networks is a useful tool for clustering groups of genes enriched for common molecular or cellular functions [1]. The relative importance of genes within networks can frequently be inferred by the degree of connectivity, with those displaying high connectivity being significantly more likely to be associated with specific molecular functions [2]. Previously we utilized cross-species network analysis to identify two network modules that were significantly associated with distant metastasis free survival in breast cancer. Here, we validate one of the highly connected genes as a metastasis associated gene. Tpx2, the most highly connected gene within a proliferation network specifically prognostic for estrogen receptor positive (ER+) breast cancers, enhances metastatic disease, but in a tumor autonomous, proliferation-independent manner. Histologic analysis suggests instead that variation of TPX2 levels within disseminated tumor cells may influence the transition between dormant to actively proliferating cells in the secondary site. These results support the co-expression network approach for identification of new metastasis-associated genes to provide new information regarding the etiology of breast cancer progression and metastatic disease.     

Common minor histocompatibility antigen discovery based upon patient clinical outcomes and genomic data

Background

Minor histocompatibility antigens (mHA) mediate much of the graft vs. leukemia (GvL) effect and graft vs. host disease (GvHD) in patients who undergo allogeneic stem cell transplantation (SCT) [1], [2], [3], [4]. Therapeutic decision making and treatments [5] based upon mHAs will require the evaluation of multiple candidate mHAs and the selection of those with the potential to have the greatest impact on clinical outcomes. We hypothesized that common, immunodominant mHAs, which are presented by HLA-A, B, and C molecules, can mediate clinically significant GvL and/or GvHD, and that these mHAs can be identified through association of genomic data with clinical outcomes.

Methodology/Principal Findings

Because most mHAs result from donor/recipient cSNP disparities, we genotyped 57 myeloid leukemia patients and their donors at 13,917 cSNPs [6]. We correlated the frequency of genetically predicted mHA disparities with clinical evidence of an immune response and then computationally screened all peptides mapping to the highly associated cSNPs for their ability to bind to HLA molecules. As proof-of-concept, we analyzed one predicted antigen, T4A, whose mHA mismatch trended towards improved overall and disease free survival in our cohort. T4A mHA mismatches occurred at the maximum theoretical frequency for any given SCT. T4A-specific CD8+ T lymphocytes (CTLs) were detected in 3 of 4 evaluable post-transplant patients predicted to have a T4A mismatch.

Conclusions/Significance

Our method is the first to combine clinical outcomes data with genomics and bioinformatics methods to predict and confirm a mHA. Refinement of this method should enable the discovery of clinically relevant mHAs in the majority of transplant patients and possibly lead to novel immunotherapeutics [5].     

Reference-Free Comparative Genomics of 174 Chloroplasts

Direct analysis of unassembled genomic data could greatly increase the power of short read DNA sequencing technologies and allow comparative genomics of organisms without a completed reference available. Here, we compare 174 chloroplasts by analyzing the taxanomic distribution of short kmers across genomes [1]. We then assemble de novo contigs centered on informative variation. The localized de novo contigs can be separated into two major classes: tip = unique to a single genome and group = shared by a subset of genomes. Prior to assembly, we found that ∼18% of the chloroplast was duplicated in the inverted repeat (IR) region across a four-fold difference in genome sizes, from a highly reduced parasitic orchid [2] to a massive algal chloroplast [3], including gnetophytes [4] and cycads [5]. The conservation of this ratio between single copy and duplicated sequence was basal among green plants, independent of photosynthesis and mechanism of genome size change, and different in gymnosperms and lower plants. Major lineages in the angiosperm clade differed in the pattern of shared kmers and de novo contigs. For example, parasitic plants demonstrated an expected accelerated overall rate of evolution, while the hemi-parasitic genomes contained a great deal more novel sequence than holo-parasitic plants, suggesting different mechanisms at different stages of genomic contraction. Additionally, the legumes are diverging more quickly and in different ways than other major families. Small duplicated fragments of the rrn23 genes were deeply conserved among seed plants, including among several species without the IR regions, indicating a crucial functional role of this duplication. Localized de novo assembly of informative kmers greatly reduces the complexity of large comparative analyses by confining the analysis to a small partition of data and genomes relevant to the specific question, allowing direct analysis of next-gen sequence data from previously unstudied genomes and rapid discovery of informative candidate regions.     

Land-Applied Goat Manure as a Source of Human Q-Fever in the Netherlands, 2006–2010

Studies have shown a link between Q-fever positive farms (QFPFs) and community cases of human Q-fever. Our study is the first to investigate the potential role of contaminated land-applied manure in human Q-fever, based on a large set of nationwide notification and farm management data. Time between manure application and disease onset in geographically linked notified human cases coincided with the incubation period of Q-fever. Proximity of contaminated land parcels predicted human cases better than proximity of QFPFs (80% vs. 58%, 0–5 km in 2009). Incidence around QFPFs and contaminated land parcels decreased with distance, but not around non-contaminated land parcels. Incidence was higher around contaminated land parcels than non-contaminated land parcels (RR = [10],95%CI = [7], [1]–[14,2]). Our findings deliver evidence that, apart from QFPFs, land-applied contaminated manure may be another source of human Q-fever.     

Epidemiology and Genetic Diversity of Rotavirus Strains in Children with Acute Gastroenteritis in Lahore,Pakistan

Pakistan harbors high disease burden of gastro-enteric infections with majority of these caused by rotavirus. Unfortunately, lack of proper surveillance programs and laboratory facilities have resulted in scarcity of available data on rotavirus associated disease burden and epidemiological information in the country. We investigated 1306 stool samples collected over two years (2008–2009) from hospitalized children under 5 years of age for the presence of rotavirus strains and its genotypic diversity in Lahore. The prevalence rate during 2008 and 2009 was found to be 34% (n = 447 out of 1306). No significant difference was found between different age groups positive for rotavirus (p>0.05). A subset of EIA positive samples was further screened for rotavirus RNA through RT-PCR and 44 (49.43%) samples, out of total 89 EIA positive samples, were found positive. G and P type prevalence was found as follows: G1P [4] = 3(6.81%); G1P [6] = 9(20.45%); G1P [8] = 1(2.27%); G2P [4] = 21(47.72%); G2P [8] = 1(2.27%); G9P [4] = 1(2.27%); G9P [6] = 1(2.27%) and G9P [8] = 7(15.90%). Phylogenetic analysis revealed that the VP7 and VP4 sequences clustered closely with the previously detected strains in the country as well as Belgian rotaviruses. Antigenic characterization was performed by analyzing major epitopes in the immunodominant VP7 and VP4 gene segments. Although the neutralization conferring motifs were found variable between the Pakistani strains and the two recommended vaccines strains (Rotarix™ and RotaTeq™), we validate the use of rotavirus vaccine in Pakistan based on the proven and recognized vaccine efficacy across the globe. Our findings constitute the first report on rotavirus’ genotype diversity, their phylogenetic relatedness and epidemiology during the pre-vaccination era in Lahore, Pakistan and support the immediate introduction of rotavirus vaccine in the routine immunization program of the country.     

Time-Delayed Mutual Information of the Phase as a Measure of Functional Connectivity

We propose a time-delayed mutual information of the phase for detecting nonlinear synchronization in electrophysiological data such as MEG. Palus already introduced the mutual information as a measure of synchronization [1]. To obtain estimates on small data-sets as reliably as possible, we adopt the numerical implementation as proposed by Kraskov and colleagues [2]. An embedding with a parametric time-delay allows a reconstruction of arbitrary nonstationary connective structures – so-called connectivity patterns – in a wide class of systems such as coupled oscillatory or even purely stochastic driven processes [3]. By using this method we do not need to make any assumptions about coupling directions, delay times, temporal dynamics, nonlinearities or underlying mechanisms. For verifying and refining the methods we generate synthetic data-sets by a mutual amplitude coupled network of Rössler oscillators with an a-priori known connective structure. This network is modified in such a way, that the power-spectrum forms a power law, which is also observed in electrophysiological recordings. The functional connectivity measure is tested on robustness to additive uncorrelated noise and in discrimination of linear mixed input data. For the latter issue a suitable de-correlation technique is applied. Furthermore, the compatibility to inverse methods for a source reconstruction in MEG such as beamforming techniques is controlled by dedicated dipole simulations. Finally, the method is applied on an experimental MEG recording.     

Characterization of Antibacterial and Hemolytic Activity of Synthetic Pandinin 2 Variants and Their Inhibition against Mycobacterium tuberculosis

The contention and treatment of Mycobacterium tuberculosis and other bacteria that cause infectious diseases require the use of new type of antibiotics. Pandinin 2 (Pin2) is a scorpion venom antimicrobial peptide highly hemolytic that has a central proline residue. This residue forms a structural “kink” linked to its pore-forming activity towards human erythrocytes. In this work, the residue Pro14 of Pin2 was both substituted and flanked using glycine residues (P14G and P14GPG) based on the low hemolytic activities of antimicrobial peptides with structural motifs Gly and GlyProGly such as magainin 2 and ponericin G1, respectively. The two Pin2 variants showed antimicrobial activity against E. coli, S. aureus, and M. tuberculosis. However, Pin2 [GPG] was less hemolytic (30%) than that of Pin2 [G] variant. In addition, based on the primary structure of Pin2 [G] and Pin2 [GPG], two short peptide variants were designed and chemically synthesized keeping attention to their physicochemical properties such as hydrophobicity and propensity to adopt alpha-helical conformations. The aim to design these two short antimicrobial peptides was to avoid the drawback cost associated to the synthesis of peptides with large sequences. The short Pin2 variants named Pin2 [14] and Pin2 [17] showed antibiotic activity against E. coli and M. tuberculosis. Besides, Pin2 [14] presented only 25% of hemolysis toward human erythrocytes at concentrations as high as 100 µM, while the peptide Pin2 [17] did not show any hemolytic effect at the same concentration. Furthermore, these short antimicrobial peptides had better activity at molar concentrations against multidrug resistance M. tuberculosis than that of the conventional antibiotics ethambutol, isoniazid and rifampicin. Therefore, Pin2 [14] and Pin2 [17] have the potential to be used as an alternative antibiotics and anti-tuberculosis agents with reduced hemolytic effects.     

Fear Inhibition in High Trait Anxiety

Trait anxiety is recognized as an individual risk factor for the development of anxiety disorders but the neurobiological mechanisms remain unknown. Here we test whether trait anxiety is associated with impaired fear inhibition utilizing the AX+/BX- conditional discrimination procedure that allows for the independent evaluation of startle fear potentiation and inhibition of fear [1]. Sixty undergraduate students participated in the study - High Trait Anxious: n = 28 and Low Trait Anxious: n = 32. We replicated earlier findings that a transfer of conditioned inhibition for startle responses requires contingency awareness. However, contrary to the fear inhibition hypothesis, our data suggest that high trait anxious individuals show a normal fear inhibition of conditioned startle responding. Only at the cognitive level the high trait anxious individuals showed evidence for impaired inhibitory learning of the threat cue. Together with other findings where impaired fear inhibition was only observed in those PTSD patients who were either high on hyperarousal symptoms [2] or with current anxiety symptoms [3], we question whether impaired fear inhibition is a biomarker for the development of anxiety disorders.