A harvest of weeds yields insight into a case of contemporary evolution

When Charles Darwin was exploring the idea of evolution via natural selection, he looked to domesticated species, with the opening chapter of The Origin of Species titled ‘Variation Under Domestication’ (Darwin 1859 ). Domesticated species such as crops are a great example of artificial selection, which Darwin realized was analogous to natural selection. But growing among those carefully selected crop varieties are the unwelcome and unwanted plants we call weeds. Despite the importance of weeds and long‐standing interest in their evolution (Baker 1974 ), we still know little about how agricultural weeds evolve, and we often fail to take evolution into account when attempting to manage them (Neve et al. 2009 ). Agricultural weeds are subjected to the unique conditions of farm fields, such as frequent soil disturbance and the addition of water and nutrients. They are also confronted with aggressive attempts at their removal via herbicides and mechanical means. As such, they are under intense demographic and selective pressure and can potentially rapidly evolve in response. In this issue of Molecular Ecology, Kuester and co‐authors make a rare attempt to understand contemporary evolution in an agricultural weed (Kuester et al. 2016 ). They do so using the powerful resurrection approach of comparing ancestors and descendants under common conditions (Franks et al. 2008 ). They sampled multiple populations of the weedy plant Ipomoea purpurea at two points in time. A comparison of these greenhouse‐grown ancestor and descendent populations showed that, over time, populations had lost significant levels of neutral genetic diversity, consistent with genetic bottlenecks. The authors also found a slight increase, on average, of resistance to the herbicide glyphosate, which is the active ingredient in Roundup^®. This work is one of a growing number of studies demonstrating rapid evolution in natural populations (Thompson 2013 ) and also reveals evidence of both selection and drift in populations of an agricultural weed.     

Comparative genomic analysis of clinical and environmental strains provides insight into the pathogenicity and evolution of Vibrio parahaemolyticus

Background

Vibrio parahaemolyticus is a Gram-negative halophilic bacterium. Infections with the bacterium could become systemic and can be life-threatening to immunocompromised individuals. Genome sequences of a few clinical isolates of V. parahaemolyticus are currently available, but the genome dynamics across the species and virulence potential of environmental strains on a genome-scale have not been described before.

Results

Here we present genome sequences of four V. parahaemolyticus clinical strains from stool samples of patients and five environmental strains in Hong Kong. Phylogenomics analysis based on single nucleotide polymorphisms revealed a clear distinction between the clinical and environmental isolates. A new gene cluster belonging to the biofilm associated proteins of V. parahaemolyticus was found in clincial strains. In addition, a novel small genomic island frequently found among clinical isolates was reported. A few environmental strains were found harboring virulence genes and prophage elements, indicating their virulence potential. A unique biphenyl degradation pathway was also reported. A database for V. parahaemolyticus (http://kwanlab.bio.cuhk.edu.hk/vp) was constructed here as a platform to access and analyze genome sequences and annotations of the bacterium.

Conclusions

We have performed a comparative genomics analysis of clinical and environmental strains of V. parahaemolyticus. Our analyses could facilitate understanding of the phylogenetic diversity and niche adaptation of this bacterium.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1135) contains supplementary material, which is available to authorized users.     

Simulated tempering yields insight into the low-resolution Rosetta scoring functions

Rosetta is a structure prediction package that has been employed successfully in numerous protein design and other applications.1 Previous reports have attributed the current limitations of the Rosetta de novo structure prediction algorithm to inadequate sampling, particularly during the low-resolution phase.2-5 Here, we implement the Simulated Tempering (ST) sampling algorithm67 in Rosetta to address this issue. ST is intended to yield canonical sampling by inducing a random walk in temperatures space such that broad sampling is achieved at high temperatures and detailed exploration of local free energy minima is achieved at low temperatures. ST should therefore visit basins in accordance with their free energies rather than their energies and achieve more global sampling than the localized scheme currently implemented in Rosetta. However, we find that ST does not improve structure prediction with Rosetta. To understand why, we carried out a detailed analysis of the low-resolution scoring functions and find that they do not provide a strong bias towards the native state. In addition, we find that both ST and standard Rosetta runs started from the native state are biased away from the native state. Although the low-resolution scoring functions could be improved, we propose that working entirely at full-atom resolution is now possible and may be a better option due to superior native-state discrimination at full-atom resolution. Such an approach will require more attention to the kinetics of convergence, however, as functions capable of native state discrimination are not necessarily capable of rapidly guiding non-native conformations to the native state.     

A new insight into the pathogenesis of filarial disease

Filariasis is a major public health problem throughout many regions of the tropics. The disease is caused by several species of filarial nematode including Wuchereria bancrofti and Brugia malayi, the agents of lymphatic filariasis, and Onchocerca volvulus, the cause of 'riverblindness'. Disease caused by these worms varies depending on the tissue location of the parasite, and is associated with episodes of acute and chronic inflammation. These pathologies, including elephantiasis and blindness, rank among the most disabling in the world. Studies aimed at characterizing the molecular nature of the inflammatory stimuli derived from filarial nematodes uncovered a long forgotten secret, their symbiont Wolbachia. LPS-like molecules from these intracellular bacteria are responsible for potent inflammatory responses from macrophages and in animal models of filarial disease. Wolbachia has also been associated with severe inflammatory reactions to filarial chemotherapy, being released into the blood following the death of the parasite. Recent studies in animal models even implicate Wolbachia in the onset of lymphodema and blindness. Taken together these studies suggest a major role for Wolbachia in the pathogenesis of filarial disease. It may be possible, through the use of antibiotic therapy, to clear worms of their bacteria, in the hope that this will prevent the onset and development of filarial pathology.     

Conservation analysis of the CydX protein yields insights into small protein identification and evolution

Background

The reliable identification of proteins containing 50 or fewer amino acids is difficult due to the limited information content in short sequences. The 37 amino acid CydX protein in Escherichia coli is a member of the cytochrome bd oxidase complex, an enzyme found throughout Eubacteria. To investigate the extent of CydX conservation and prevalence and evaluate different methods of small protein homologue identification, we surveyed 1095 Eubacteria species for the presence of the small protein.

Results

Over 300 homologues were identified, including 80 unannotated genes. The ability of both closely-related and divergent homologues to complement the E. coli ΔcydX mutant supports our identification techniques, and suggests that CydX homologues retain similar function among divergent species. However, sequence analysis of these proteins shows a great degree of variability, with only a few highly-conserved residues. An analysis of the co-variation between CydX homologues and their corresponding cydA and cydB genes shows a close synteny of the small protein with the CydA long Q-loop. Phylogenetic analysis suggests that the cydABX operon has undergone horizontal gene transfer, although the cydX gene likely evolved in a progenitor of the Alpha, Beta, and Gammaproteobacteria. Further investigation of cydAB operons identified two additional conserved hypothetical small proteins: CydY encoded in CydA_Qlong operons that lack cydX, and CydZ encoded in more than 150 CydA_Qshort operons.

Conclusions

This study provides a systematic analysis of bioinformatics techniques required for the unique challenges present in small protein identification and phylogenetic analyses. These results elucidate the prevalence of CydX throughout the Proteobacteria, provide insight into the selection pressure and sequence requirements for CydX function, and suggest a potential functional interaction between the small protein and the CydA Q-loop, an enigmatic domain of the cytochrome bd oxidase complex. Finally, these results identify other conserved small proteins encoded in cytochrome bd oxidase operons, suggesting that small protein subunits may be a more common component of these enzymes than previously thought.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-946) contains supplementary material, which is available to authorized users.     

Telocytes: new insight into the pathogenesis of gallstone disease

The major mechanisms of gallstone formation include biliary cholesterol hypersecretion, supersaturation and crystallization, mucus hypersecretion, gel formation and bile stasis. Gallbladder hypomotility seems to be a key event that triggers the precipitation of cholesterol microcrystals from supersaturated lithogenic bile. Telocytes, a new type of interstitial cells, have been recently identified in many organs, including gallbladder. Considering telocyte functions, it is presumed that these cells might be involved in the signalling processes. The purpose of this study was to correlate the quantity of telocytes in the gallbladder with the lithogenicity of bile. Gallbladder specimens were collected from 24 patients who underwent elective laparoscopic cholecystectomy for symptomatic gallstone disease. The control group consisted of 25 consecutive patients who received elective treatment for pancreatic head tumours. Telocytes were visualized in paraffin sections of gallbladders with double immunofluorescence using primary antibodies against c‐Kit (anti‐CD117) and anti‐mast cell tryptase. Cholesterol, phospholipid and bile acid levels were measured in gallbladder bile. The number of telocytes in the gallbladder wall was significantly lower in the study group than that in the control group (3.03 ± 1.43 versus 6.34 ± 1.66 cell/field of view in the muscularis propria, P < 0.001) and correlated with a significant increase in the cholesterol saturation index. The glycocholic and taurocholic acid levels were significantly elevated in the control subjects compared with the study group. The results suggest that bile composition may play an important role in the reduction in telocytes density in the gallbladder.     

Comparative pathogenesis of Mycobacterium marinum and Mycobacterium tuberculosis

A thorough understanding of Mycobacterium tuberculosis pathogenesis in humans has been elusive in part because of imperfect surrogate laboratory hosts, each with its own idiosyncrasies. Mycobacterium marinum is the closest genetic relative of the M. tuberculosis complex and is a natural pathogen of ectotherms. In this review, we present evidence that the similar genetic programmes of M. marinum and M. tuberculosis and the corresponding host immune responses reveal a conserved skeleton of Mycobacterium host–pathogen interactions. While both species have made niche-specific refinements, an essential framework has persisted. We highlight genetic comparisons of the two organisms and studies of M. marinum in the developing zebrafish. By pairing M. marinum with the simplified immune system of zebrafish embryos, many of the defining mechanisms of mycobacterial pathogenesis can be distilled and investigated in a tractable host/pathogen pair.     

Comparative MD analysis of the stability of transthyretin providing insight into the fibrillation mechanism

Proteins can misfold and aggregate, which is believed to be the cause of a variety of diseases, affecting very diverse organs in the body. Many questions about the nature of aggregation and the proteins that are involved in these events are still left unanswered. One of the proteins that is known to form amyloids is transthyretin (TTR), the secondary transporter of thyroxine, and transporter of retinol-binding protein. Several experimental results have helped to explain this aberrant behavior of TTR; however, structural insights of the amyloidgenic process are still lacking. Therefore, we have used all-atom MD simulation and free energy calculations to study the initial phase of this process. We have calculated the free energy changes of the initial tetramer dissociation under different conditions and in the presence of thyroxine. We show that tetramer formation is indeed only thermodynamically favorable in neutral pH conditions. We find that binding of two thyroxine molecules stabilizes the complex, and that this occurs with negative cooperativity. In addition to the energetic calculations, we have also investigated the dominant motions of the TTR and found that only the dimeric form of the protein could undergo the initial fibril formation.     

Serum proteome analysis of vivax malaria: An insight into the disease pathogenesis and host immune response

Vivax malaria is the most widely distributed human malaria resulting in 80-300 million clinical cases every year. It causes severe infection and mortality but is generally regarded as a benign disease and has not been investigated in detail. The present study aimed to perform human serum proteome analysis in a malaria endemic area in India to identify potential serum biomarkers for vivax malaria and understand host response. The proteomic analysis was performed on 16 age and gender matched subjects (vivax patients and control) in duplicate. Protein extraction protocols were optimized for large coverage of the serum proteome and to obtain high-resolution data. Identification of 67 differentially expressed and statistically significant (Student's t-test; p<0.05) protein spots was established by MALDI-TOF/TOF mass spectrometry. Many of the identified proteins such as apolipoprotein A and E, serum amyloid A and P, haptoglobin, ceruloplasmin, and hemopexin are interesting from a diagnostic point of view and could further be studied as potential serum biomarkers. The differentially expressed serum proteins in vivax malaria identified in this study were subjected to functional pathway analysis using multiple software, including Ingenuity Pathway Analysis (IPA), Protein ANalysis THrough Evolutionary Relationships (PANTHER) and Database for Annotation, Visualization and Integrated Discovery (DAVID) functional annotation tool for better understanding of the biological context of the identified proteins, their involvement in various physiological pathways and association with disease pathogenesis. Functional pathway analysis of the differentially expressed proteins suggested the modulation of multiple vital physiological pathways, including acute phase response signaling, complement and coagulation cascades, hemostasis and vitamin D metabolism pathway due to this parasitic infection. This article is part of a Special Issue entitled: Proteomics: The clinical link.     

Molecular modelling of the GIR1 branching ribozyme gives new insight into evolution of structurally related ribozymes

Twin-ribozyme introns contain a branching ribozyme (GIR1) followed by a homing endonuclease (HE) encoding sequence embedded in a peripheral domain of a group I splicing ribozyme (GIR2). GIR1 catalyses the formation of a lariat with 3 nt in the loop, which caps the HE mRNA. GIR1 is structurally related to group I ribozymes raising the question about how two closely related ribozymes can carry out very different reactions. Modelling of GIR1 based on new biochemical and mutational data shows an extended substrate domain containing a GoU pair distinct from the nucleophilic residue that dock onto a catalytic core showing a different topology from that of group I ribozymes. The differences include a core J8/7 region that has been reduced and is complemented by residues from the pre-lariat fold. These findings provide the basis for an evolutionary mechanism that accounts for the change from group I splicing ribozyme to the branching GIR1 architecture. Such an evolutionary mechanism can be applied to other large RNAs such as the ribonuclease P.     

Comparative mapping of human Chromosome 19 with the chicken shows conserved synteny and gives an insight into chromosomal evolution

Human Chromosome 19 (HSA19) is virtually completely sequenced. A complete physical contig map made up of BACs and cosmids is also available for this chromosome. It is, therefore, a rich source of information that we have used as the basis for a comparative mapping study with the chicken. Various orthologs of genes known to map to HSA19 have been mapped in the chicken. Five chicken microchromosomes (two of which were previously undefined) are seen to show conserved synteny with this chromosome, along with individual gene homologs on Chr 1 and another tiny microchromosome. Compared with the mouse, which has 12 chromosomal regions homologous to HSA19, the chicken genotype displays fewer evolutionary rearrangements. The ancestral nature of the chicken karyotype is demonstrated and may prove to be an excellent tool for studying genome evolution.     

Molecular phylogenetic analysis of Commiphora (Burseraceae) yields insight on the evolution and historical biogeography of an "impossible" genus

Expansion of the arid zone of sub-Saharan tropical Africa during the Miocene is posited as a significant contributing factor in the evolution of contemporary African flora. Nevertheless, few molecular phylogenetic studies have tested this hypothesis using reconstructed historical biogeographies of plants within this zone. Here, we present a molecular phylogeny of Commiphora, a predominantly tropical African, arid-adapted tree genus, in order to test the monophyly of its taxonomic sections and identify clades that will help direct future study of this species-rich and geographically widespread taxon. We then use multiple fossil calibrations of Commiphora phylogeny to determine the timing of well-supported diversification events within the genus and interpret these age estimates to determine the relative contribution of vicariance and dispersal in the expansion of Commiphora's geographic range. We find that Commiphora is sister to Vietnamese Bursera tonkinensis and that its crown group radiation corresponds with the onset of the Miocene.     

New insight into the role of dendritic cells in malaria immune pathogenesis

The mechanism by which the host develops protective immunity to malaria remains poorly understood. Dendritic cells (DCs) are central to the initiation and regulation of the adaptive immune response. Modulation of DC function might enable Plasmodium to evade the immune system. Millington et al. propose one mechanism by which malaria inhibits DC-T-cell interactions without interfering directly with T-cell receptor engagement. The consequence is a decrease in the co-stimulation required to develop an effective immune response.     

Comparative phylogenetic analysis of aquaporins provides insight into the gene family expansion and evolution in plants and their role in drought tolerant and susceptible chickpea cultivars

Aquaporins (AQPs) are water channel proteins that play a significant role in drought stress. Although the AQPs identified in multiple plant species, there is no detailed evolutionary and comparative study of AQPs regarding chickpea plant. The current study involved evolutionary analyses coupled with promoter and expression analyses of chickpea AQPs (CaAQPs). A total of 924 non-redundant AQPs were studied in 24 plant species including algae, mosses, lycophytes, monocots and dicots. Phylogenetic analysis demonstrated a clear divergence of eight AQP subfamilies (LIPs, SIPs, GIPs, NIPs, XIPs, PIPs, HIPs and TIPs). The comparative phylogenetic trees of AQP subfamilies among Arabidopsis, soybean, common bean, maize and chickpea demonstrated that the AQPs were highly species-specific. Interestingly, the dual NPA motif was conserved in all species. However, the ar/R selectivity filter signatures [W/T/S/N/G/A]-[V/S/L/I/A]-[S/G/A]-R (in NIPs), F-H-T-R (in PIPs), [H/N/Q/S]-[A/I/L/S/V]-[A/G]-[A/C/L/M/R/V] (in TIPs) and [V/I/L/M]-[V/I/A/F/M]-[A/S/F/C]-[N/F/L/I/A/S (in SIPs) were found in five species. Moreover, the Froger's positions (P1-P5) were found as [F/L/Y]-[S/T]-A-Y-[L/I/M/V/F] (in NIPs), [Q/E/M]-S-A-F-W (in PIPs), [A/L/S/T/V]-[A/C/N/S/T/V]-[P/R/S]-[Y/N/F]-[W/Q] (in TIPs) and [I/M/F]-[A/V]-[A/V]-Y-W (in SIPs). The MEME motif analyses showed that most of the motifs were specific to subfamily and subgroups. Tissue-specific expression profiling of CaAQPs revealed that CaTIPs and CaPIPs are highly expressed in most of the tissues, while CaNIPs and CaSIPs have low expression. In promoter analysis of CaAQPs, multiple stress-related cis-acting elements e.g. MYB, MYC, ABRE, etc. were found. Semi-quantitative RT-PCR analysis showed that CaPIP2;3 and CaNIP3;1 are positive regulator, while CaSIP1;1 and CaPIP2;1 have a negative role in drought tolerance. The findings and implications of this study are discussed in detail.     

The amphioxus genome provides unique insight into the evolution of immunity

Immune systems evolve as essential strategies to maintain homeostasis with the environment, prevent microbial assault and recycle damaged host tissues. The immune system is composed of two components, innate and adaptive immunity. The former is common to all animals while the latter consists of a vertebrate-specific system that relies on somatically derived lymphocytes and is associated with near limitless genetic diversity as well as long-term memory. Deuterostome invertebrates provide a view of immune repertoires in phyla that immediately predate the origins of vertebrates. Genomic studies in amphioxus, a cephalochordate, have revealed homologs of genes encoding most innate immune receptors found in vertebrates; however, many of the gene families have undergone dramatic expansions, greatly increasing the innate immune repertoire. In addition, domain-swapping accounts for the innovation of new predicted pathways of receptor function. In both amphioxus and Ciona, a urochordate, the VCBPs (variable region containing chitin-binding proteins), which consist of immunoglobulin V (variable) and chitin binding domains, mediate recognition through the V domains. The V domains of VCBPs in amphioxus exhibit high levels of allelic complexity that presumably relate to functional specificity. Various features of the amphioxus immune repertoire reflect novel selective pressures, which likely have resulted in innovative strategies. Functional genomic studies underscore the value of amphioxus as a model for studying innate immunity and may help reveal how unique relationships between innate immune receptors and both pathogens and symbionts factored in the evolution of adaptive immune systems.