An efficient pipeline for ancient DNA mapping and recovery of endogenous ancient DNA from whole‐genome sequencing data

Ancient DNA research has developed rapidly over the past few decades due to improvements in PCR and next‐generation sequencing (NGS) technologies, but challenges still exist. One major challenge in relation to ancient DNA research is to recover genuine endogenous ancient DNA sequences from raw sequencing data. This is often difficult due to degradation of ancient DNA and high levels of contamination, especially homologous contamination that has extremely similar genetic background with that of the real ancient DNA. In this study, we collected whole‐genome sequencing (WGS) data from 6 ancient samples to compare different mapping algorithms. To further explore more effective methods to separate endogenous DNA from homologous contaminations, we attempted to recover reads based on ancient DNA specific characteristics of deamination, depurination, and DNA fragmentation with different parameters. We propose a quick and improved pipeline for separating endogenous ancient DNA while simultaneously decreasing homologous contaminations to very low proportions. Our goal in this research was to develop useful recommendations for ancient DNA mapping and for separation of endogenous DNA to facilitate future studies of ancient DNA.     

Study protocol to examine the effects of spaceflight and a spaceflight analog on neurocognitive performance: extent,longevity, and neural bases

Background

Long duration spaceflight (i.e., 22 days or longer) has been associated with changes in sensorimotor systems, resulting in difficulties that astronauts experience with posture control, locomotion, and manual control. The microgravity environment is an important causal factor for spaceflight induced sensorimotor changes. Whether spaceflight also affects other central nervous system functions such as cognition is yet largely unknown, but of importance in consideration of the health and performance of crewmembers both in- and post-flight. We are therefore conducting a controlled prospective longitudinal study to investigate the effects of spaceflight on the extent, longevity and neural bases of sensorimotor and cognitive performance changes. Here we present the protocol of our study.

Methods/design

This study includes three groups (astronauts, bed rest subjects, ground-based control subjects) for which each the design is single group with repeated measures. The effects of spaceflight on the brain will be investigated in astronauts who will be assessed at two time points pre-, at three time points during-, and at four time points following a spaceflight mission of six months. To parse out the effect of microgravity from the overall effects of spaceflight, we investigate the effects of seventy days head-down tilted bed rest. Bed rest subjects will be assessed at two time points before-, two time points during-, and three time points post-bed rest. A third group of ground based controls will be measured at four time points to assess reliability of our measures over time. For all participants and at all time points, except in flight, measures of neurocognitive performance, fine motor control, gait, balance, structural MRI (T1, DTI), task fMRI, and functional connectivity MRI will be obtained. In flight, astronauts will complete some of the tasks that they complete pre- and post flight, including tasks measuring spatial working memory, sensorimotor adaptation, and fine motor performance. Potential changes over time and associations between cognition, motor-behavior, and brain structure and function will be analyzed.

Discussion

This study explores how spaceflight induced brain changes impact functional performance. This understanding could aid in the design of targeted countermeasures to mitigate the negative effects of long-duration spaceflight.

     

Climate change and the optimal flowering time of annual plants in seasonal environments

Long‐term phenology monitoring has documented numerous examples of changing flowering dates during the last century. A pivotal question is whether these phenological responses are adaptive or not under directionally changing climatic conditions. We use a classic dynamic growth model for annual plants, based on optimal control theory, to find the fitness‐maximizing flowering time, defined as the switching time from vegetative to reproductive growth. In a typical scenario of global warming, with advanced growing season and increased productivity, optimal flowering time advances less than the start of the growing season. Interestingly, increased temporal spread in production over the season may either advance or delay the optimal flowering time depending on overall productivity or season length. We identify situations where large phenological changes are necessary for flowering time to remain optimal. Such changes also indicate changed selection pressures. In other situations, the model predicts advanced phenology on a calendar scale, but no selection for early flowering in relation to the start of the season. We also show that the optimum is more sensitive to increased productivity when productivity is low than when productivity is high. All our results are derived using a general, graphical method to calculate the optimal flowering time applicable for a large range of shapes of the seasonal production curve. The model can thus explain apparent maladaptation in phenological responses in a multitude of scenarios of climate change. We conclude that taking energy allocation trade‐offs and appropriate time scales into account is critical when interpreting phenological patterns.     

Mycobacterium tuberculosis Bacteremia in a Cohort of HIV-Infected Patients Hospitalized with Severe Sepsis in Uganda–High Frequency,Low Clinical Sand Derivation of a Clinical Prediction Score

Background

When manifested as Mycobacterium tuberculosis (MTB) bacteremia, disseminated MTB infection clinically mimics other serious blood stream infections often hindering early diagnosis and initiation of potentially life-saving anti-tuberculosis therapy. In a cohort of hospitalized HIV-infected Ugandan patients with severe sepsis, we report the frequency, management and outcomes of patients with MTB bacteremia and propose a risk score based on clinical predictors of MTB bacteremia.

Methods

We prospectively enrolled adult patients with severe sepsis at two Ugandan hospitals and obtained blood cultures for MTB identification. Multivariable logistic regression modeling was used to determine predictors of MTB bacteremia and to inform the stratification of patients into MTB bacteremia risk categories based on relevant patient characteristics.

Results

Among 368 HIV-infected patients with a syndrome of severe sepsis, eighty-six (23%) had MTB bacteremia. Patients with MTB bacteremia had a significantly lower median CD4 count (17 vs 64 lymphocytes/mm³, p<0.001) and a higher 30-day mortality (53% vs 32%, p = 0.001) than patients without MTB bacteremia. A minority of patients with MTB bacteremia underwent standard MTB diagnostic testing (24%) or received empiric anti-tuberculosis therapy (15%). Independent factors associated with MTB bacteremia included male sex, increased heart rate, low CD4 count, absence of highly active anti-retroviral therapy, chief complaint of fever, low serum sodium and low hemoglobin. A risk score derived from a model containing these independent predictors had good predictive accuracy [area under the curve = 0.85, 95% CI 0.80–0.89].

Conclusions

Nearly 1 in 4 adult HIV-infected patients hospitalized with severe sepsis in 2 Ugandan hospitals had MTB bacteremia. Among patients in whom MTB was suspected, standard tests for diagnosing pulmonary MTB were inaccurate for correctly classifying patients with or without bloodstream MTB infection. A MTB bacteremia risk score can improve early diagnosis of MTB bacteremia particularly in settings with increased HIV and MTB co-infection.     

The Impact of Motor Axon Misdirection and Attrition on Behavioral Deficit Following Experimental Nerve Injuries

Peripheral nerve transection and neuroma-in-continuity injuries are associated with permanent functional deficits, often despite successful end-organ reinnervation. Axonal misdirection with non-specific reinnervation, frustrated regeneration and axonal attrition are believed to be among the anatomical substrates that underlie the poor functional recovery associated with these devastating injuries. Yet, functional deficits associated with axonal misdirection in experimental neuroma-in-continuity injuries have not yet been studied. We hypothesized that experimental neuroma-in-continuity injuries would result in motor axon misdirection and attrition with proportional persistent functional deficits. The femoral nerve misdirection model was exploited to assess major motor pathway misdirection and axonal attrition over a spectrum of experimental nerve injuries, with neuroma-in-continuity injuries simulated by the combination of compression and traction forces in 42 male rats. Sciatic nerve injuries were employed in an additional 42 rats, to evaluate the contribution of axonal misdirection to locomotor deficits by a ladder rung task up to 12 weeks. Retrograde motor neuron labeling techniques were utilized to determine the degree of axonal misdirection and attrition. Characteristic histological neuroma-in-continuity features were demonstrated in the neuroma-in-continuity groups and poor functional recovery was seen despite successful nerve regeneration and muscle reinnervation. Good positive and negative correlations were observed respectively between axonal misdirection (p<.0001, r²=.67), motor neuron counts (attrition) (p<.0001, r²=.69) and final functional deficits. We demonstrate prominent motor axon misdirection and attrition in neuroma-in-continuity and transection injuries of mixed motor nerves that contribute to the long-term functional deficits. Although widely accepted in theory, to our knowledge, this is the first experimental evidence to convincingly demonstrate these correlations with data inclusive of the neuroma-in-continuity spectrum. This work emphasizes the need to focus on strategies that promote both robust and accurate nerve regeneration to optimize functional recovery. It also demonstrates that clinically relevant neuroma-in-continuity injuries can now also be subjected to experimental investigation.     

Validation of a Remote Sensing Model to Identify Simulium damnosum s.l. Breeding Sites in Sub-Saharan Africa

Background

Recently, most onchocerciasis control programs have begun to focus on elimination. Developing an effective elimination strategy relies upon accurately mapping the extent of endemic foci. In areas of Africa that suffer from a lack of infrastructure and/or political instability, developing such accurate maps has been difficult. Onchocerciasis foci are localized near breeding sites for the black fly vectors of the infection. The goal of this study was to conduct ground validation studies to evaluate the sensitivity and specificity of a remote sensing model developed to predict S. damnosum s.l. breeding sites.

Methodology/Principal Findings

Remote sensing images from Togo were analyzed to identify areas containing signature characteristics of S. damnosum s.l. breeding habitat. All 30 sites with the spectral signature were found to contain S. damnosum larvae, while 0/52 other sites judged as likely to contain larvae were found to contain larvae. The model was then used to predict breeding sites in Northern Uganda. This area is hyper-endemic for onchocerciasis, but political instability had precluded mass distribution of ivermectin until 2009. Ground validation revealed that 23/25 sites with the signature contained S. damnosum larvae, while 8/10 sites examined lacking the signature were larvae free. Sites predicted to have larvae contained significantly more larvae than those that lacked the signature.

Conclusions/Significance

This study suggests that a signature extracted from remote sensing images may be used to predict the location of S. damnosum s.l. breeding sites with a high degree of accuracy. This method should be of assistance in predicting communities at risk for onchocerciasis in areas of Africa where ground-based epidemiological surveys are difficult to implement.     

Reconstructing the Population Genetic History of the Caribbean

The Caribbean basin is home to some of the most complex interactions in recent history among previously diverged human populations. Here, we investigate the population genetic history of this region by characterizing patterns of genome-wide variation among 330 individuals from three of the Greater Antilles (Cuba, Puerto Rico, Hispaniola), two mainland (Honduras, Colombia), and three Native South American (Yukpa, Bari, and Warao) populations. We combine these data with a unique database of genomic variation in over 3,000 individuals from diverse European, African, and Native American populations. We use local ancestry inference and tract length distributions to test different demographic scenarios for the pre- and post-colonial history of the region. We develop a novel ancestry-specific PCA (ASPCA) method to reconstruct the sub-continental origin of Native American, European, and African haplotypes from admixed genomes. We find that the most likely source of the indigenous ancestry in Caribbean islanders is a Native South American component shared among inland Amazonian tribes, Central America, and the Yucatan peninsula, suggesting extensive gene flow across the Caribbean in pre-Columbian times. We find evidence of two pulses of African migration. The first pulse—which today is reflected by shorter, older ancestry tracts—consists of a genetic component more similar to coastal West African regions involved in early stages of the trans-Atlantic slave trade. The second pulse—reflected by longer, younger tracts—is more similar to present-day West-Central African populations, supporting historical records of later transatlantic deportation. Surprisingly, we also identify a Latino-specific European component that has significantly diverged from its parental Iberian source populations, presumably as a result of small European founder population size. We demonstrate that the ancestral components in admixed genomes can be traced back to distinct sub-continental source populations with far greater resolution than previously thought, even when limited pre-Columbian Caribbean haplotypes have survived.     

Sh3glb1/Bif-1 and mitophagy

Evasion of apoptosis, which enables cells to survive and proliferate under metabolic stress, is one of the hallmarks of cancer. We have recently reported that SH3GLB1/Bif-1 functions as a haploinsufficient tumor suppressor to prevent the acquisition of apoptosis resistance and malignant transformation during Myc-driven lymphomagenesis. SH3GLB1 is a membrane curvature-inducing protein that interacts with BECN1 though UVRAG and regulates the post-Golgi trafficking of membrane-integrated ATG9A for autophagy. At the premalignant stage, allelic loss of Sh3glb1 enhances Myc-induced chromosomal instability and results in the upregulation of anti-apoptotic proteins, including MCL1 and BCL2L1. Notably, we found that Sh3glb1 haploinsufficiency increases mitochondrial mass in overproliferated prelymphomatous Eμ-Myc cells. Moreover, loss of Sh3glb1 suppresses autophagy-dependent mitochondrial clearance (mitophagy) in PARK2/Parkin-expressing mouse embryonic fibroblasts (MEFs) treated with the mitochondrial uncoupler CCCP. Interestingly, PARK2-expressing Sh3glb1-deficient cells accumulate ER-associated immature autophagosome-like structures after treatment with CCCP. Taken together, we propose a model of mitophagy in which SH3GLB1 together with the class III phosphatidylinositol 3-kinase complex II (PIK3C3CII) (PIK3R4-PIK3C3-BECN1-UVRAG) regulates the trafficking of ATG9A-containing Golgi-derived membranes (A9⁺GDMs) to damaged mitochondria for autophagosome formation to counteract oncogene-driven tumorigenesis.