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.     

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.     

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.     

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.     

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].     

Culture-Independent Microbiological Analysis of Foley Urinary Catheter Biofilms

Background

Prevention of catheter-associated urinary tract infection (CAUTI), a leading cause of nosocomial disease, is complicated by the propensity of bacteria to form biofilms on indwelling medical devices [1], [2], [3], [4], [5].

Methodology/Principal Findings

To better understand the microbial diversity of these communities, we report the results of a culture-independent bacterial survey of Foley urinary catheters obtained from patients following total prostatectomy. Two patient subsets were analyzed, based on treatment or no treatment with systemic fluoroquinolone antibiotics during convalescence. Results indicate the presence of diverse polymicrobial assemblages that were most commonly observed in patients who did not receive systemic antibiotics. The communities typically contained both Gram-positive and Gram-negative microorganisms that included multiple potential pathogens.

Conclusion/Significance

Prevention and treatment of CAUTI must take into consideration the possible polymicrobial nature of any particular infection.     

A Vertebrate-Specific Chp-PAK-PIX Pathway Maintains E-Cadherin at Adherens Junctions during Zebrafish Epiboly

Background

In early vertebrate development, embryonic tissues modulate cell adhesiveness and acto-myosin contractility to correctly orchestrate the complex processes of gastrulation. E-cadherin (E-cadh) is the earliest expressed cadherin and is needed in the mesendodermal progenitors for efficient migration [1], [2]. Regulatory mechanisms involving directed E-cadh trafficking have been invoked downstream of Wnt11/5 signaling [3]. This non-canonical Wnt pathway regulates RhoA-ROK/DAAM1 to control the acto-myosin network. However, in this context nothing is known of the intracellular signals that participate in the correct localization of E-cadh, other than a need for Rab5c signaling [3].

Methodology/Principal Findings

By studying loss of Chp induced by morpholino-oligonucleotide injection in zebrafish, we find that the vertebrate atypical Rho-GTPase Chp is essential for the proper disposition of cells in the early embryo. The underlying defect is not leading edge F-actin assembly (prominent in the cells of the envelope layer), but rather the failure to localize E-cadh and β-catenin at the adherens junctions. Loss of Chp results in delayed epiboly that can be rescued by mRNA co-injection, and phenocopies zebrafish E-cadh mutants [4], [5]. This new signaling pathway involves activation of an effector kinase PAK, and involvement of the adaptor PAK-interacting exchange factor PIX. Loss of signaling by any of the three components results in similar underlying defects, which is most prominent in the epithelial-like envelope layer.

Conclusions/Significance

Our current study uncovers a developmental pathway involving Chp/PAK/PIX signaling, which helps co-ordinate E-cadh disposition to promote proper cell adhesiveness, and coordinate movements of the three major cell layers in epiboly. Our data shows that without Chp signaling, E-cadh shifts to intracellular vesicles rather than the adhesive contacts needed for directed cell movement. These events may mirror the requirement for PAK2 signaling essential for the proper formation of the blood-brain barrier [6], [7].     

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.     

Should We Abandon the t-Test in the Analysis of Gene Expression Microarray Data: A Comparison of Variance Modeling Strategies

High-throughput post-genomic studies are now routinely and promisingly investigated in biological and biomedical research. The main statistical approach to select genes differentially expressed between two groups is to apply a t-test, which is subject of criticism in the literature. Numerous alternatives have been developed based on different and innovative variance modeling strategies. However, a critical issue is that selecting a different test usually leads to a different gene list. In this context and given the current tendency to apply the t-test, identifying the most efficient approach in practice remains crucial. To provide elements to answer, we conduct a comparison of eight tests representative of variance modeling strategies in gene expression data: Welch''s t-test, ANOVA [1], Wilcoxon''s test, SAM [2], RVM [3], limma [4], VarMixt [5] and SMVar [6]. Our comparison process relies on four steps (gene list analysis, simulations, spike-in data and re-sampling) to formulate comprehensive and robust conclusions about test performance, in terms of statistical power, false-positive rate, execution time and ease of use. Our results raise concerns about the ability of some methods to control the expected number of false positives at a desirable level. Besides, two tests (limma and VarMixt) show significant improvement compared to the t-test, in particular to deal with small sample sizes. In addition limma presents several practical advantages, so we advocate its application to analyze gene expression data.     

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.     

Extinction Risk Escalates in the Tropics

The latitudinal biodiversity gradient remains one of the most widely recognized yet puzzling patterns in nature [1]. Presently, the high level of extinction of tropical species, referred to as the “tropical biodiversity crisis”, has the potential to erode this pattern. While the connection between species richness, extinction, and speciation has long intrigued biologists [2], [3], these interactions have experienced increased poignancy due to their relevancy to where we should concentrate our conservation efforts. Natural extinction is a phenomenon thought to have its own latitudinal gradient, with lower extinction rates in the tropics being reported in beetles, birds, mammals, and bivalves [4]–[7]. Processes that have buffered ecosystems from high extinction rates in the past may also buffer ecosystems against disturbance of anthropogenic origin. While potential parallels between historical and present-day extinction patterns have been acknowledged, they remain only superficially explored and plant extinction patterns have been particularly neglected. Studies on the disappearances of animal species have reached conflicting conclusions, with the rate of extinction appearing either higher [8] or lower [9] in species richness hotspots. Our global study of extinction risk in vascular plants finds disproportionately higher extinction risk in tropical countries, even when indicators of human pressure (GDP, population density, forest cover change) are taken into account. Our results are at odds with the notion that the tropics represent a museum of plant biodiversity (places of historically lowered extinction) and we discuss mechanisms that may reconcile this apparent contradiction.     

A Global Analysis of the Effectiveness of Marine Protected Areas in Preventing Coral Loss

Background

A variety of human activities have led to the recent global decline of reef-building corals [1], [2]. The ecological, social, and economic value of coral reefs has made them an international conservation priority [2], [3]. The success of Marine Protected Areas (MPAs) in restoring fish populations [4] has led to optimism that they could also benefit corals by indirectly reducing threats like overfishing, which cause coral degradation and mortality [2], [5]. However, the general efficacy of MPAs in increasing coral reef resilience has never been tested.

Methodology/Principal Findings

We compiled a global database of 8534 live coral cover surveys from 1969–2006 to compare annual changes in coral cover inside 310 MPAs to unprotected areas. We found that on average, coral cover within MPAs remained constant, while coral cover on unprotected reefs declined. Although the short-term differences between unprotected and protected reefs are modest, they could be significant over the long-term if the effects are temporally consistent. Our results also suggest that older MPAs were generally more effective in preventing coral loss. Initially, coral cover continued to decrease after MPA establishment. Several years later, however, rates of coral cover decline slowed and then stabilized so that further losses stopped.

Conclusions/Significance

These findings suggest that MPAs can be a useful tool not only for fisheries management, but also for maintaining coral cover. Furthermore, the benefits of MPAs appear to increase with the number of years since MPA establishment. Given the time needed to maximize MPA benefits, there should be increased emphasis on implementing new MPAs and strengthening the enforcement of existing MPAs.     

Genetic Diversity of Circulating Rotavirus Strains in Tanzania Prior to the Introduction of Vaccination

Background

Tanzania currently rolls out vaccination against rotavirus-diarrhea, a major cause of child illness and death. As the vaccine covers a limited number of rotavirus variants, this study describes the molecular epidemiology of rotavirus among children under two years in Dar es Salaam, Tanzania, prior to implementation of vaccination.

Methods

Stool specimens, demographic and clinical information, were collected from 690 children admitted to hospital due to diarrhea (cases) and 545 children without diarrhea (controls) during one year. Controls were inpatient or children attending child health clinics. Rotavirus antigen was detected using ELISA and positive samples were typed by multiplex semi-nested PCR and sequencing.

Results

The prevalence of rotavirus was higher in cases (32.5%) than in controls (7.7%, P<0.001). The most common G genotypes were G1 followed by G8, G12, and G4 in cases and G1, G12 and G8 in controls. The Tanzanian G1 variants displayed 94% similarity with the Rotarix vaccine G1 variant. The commonest P genotypes were P[8], P[4] and P[6], and the commonest G/P combination G1 P[8] (n = 123), G8 P[4] and G12 P[6]. Overall, rotavirus prevalence was higher in cool (23.9%) than hot months (17.1%) of the year (P = 0.012). We also observed significant seasonal variation of G genotypes. Rotavirus was most frequently found in the age group of four to six months. The prevalence of rotavirus in cases was lower in stunted children (28.9%) than in non-stunted children (40.1%, P = 0.003) and lower in HIV-infected (15.4%, 4/26) than in HIV-uninfected children (55.3%, 42/76, P<0.001).

Conclusion

This pre-vaccination study shows predominance of genotype G1 in Tanzania, which is phylogenetically distantly related to the vaccine strains. We confirm the emergence of genotype G8 and G12. Rotavirus infection and circulating genotypes showed seasonal variation. This study also suggests that rotavirus may not be an opportunistic pathogen in children infected with HIV.     

Cue-Recruitment for Extrinsic Signals after Training with Low Information Stimuli

Cue-recruitment occurs when a previously ineffective signal comes to affect the perceptual appearance of a target object, in a manner similar to the trusted cues with which the signal was put into correlation during training [1], [2]. Jain, Fuller and Backus [3] reported that extrinsic signals, those not carried by the target object itself, were not recruited even after extensive training. However, recent studies have shown that training using weakened trusted cues can facilitate recruitment of intrinsic signals [4]–[7]. The current study was designed to examine whether extrinsic signals can be recruited by putting them in correlation with weakened trusted cues. Specifically, we tested whether an extrinsic visual signal, the rotary motion direction of an annulus of random dots, and an extrinsic auditory signal, direction of an auditory pitch glide, can be recruited as cues for the rotation direction of a Necker cube. We found learning, albeit weak, for visual but not for auditory signals. These results extend the generality of the cue-recruitment phenomenon to an extrinsic signal and provide further evidence that the visual system learns to use new signals most quickly when other, long-trusted cues are unavailable or unreliable.     

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.