Bacillus anthracis: Molecular taxonomy, population genetics, phylogeny and patho-evolution

Bacillus anthracis, the etiological agent of anthrax, manifests a particular bimodal lifestyle. This bacterial species alternates between short replication phases of 20-40 generations that strictly require infection of the host, normally causing death, interrupted by relatively long, mostly dormant phases as spores in the environment. Hence, the B. anthracis genome is highly homogeneous. This feature and the fact that strains from nearly all parts of the world have been analysed for canonical single nucleotide polymorphisms (canSNPs) and variable number tandem repeats (VNTRs) has allowed the development of molecular epidemiological and molecular clock models to estimate the age of major diversifications in the evolution of B. anthracis and to trace the global spread of this pathogen, which was mostly promoted by movement of domestic cattle with settlers and by international trade of contaminated animal products. From a taxonomic and phylogenetic point of view, B. anthracis is a member of the Bacillus cereus group. The differentiation of B. anthracis from B. cereus sensu strict, solely based on chromosomal markers, is difficult. However, differences in pathogenicity clearly differentiate B. anthracis from B. cereus and are marked by the strict presence of virulence genes located on the two virulence plasmids pXO1 and pXO2, which both are required by the bacterium to cause anthrax. Conversely, anthrax-like symptoms can also be caused by organisms with chromosomal features that are more closely related to B. cereus, but which carry these virulence genes on two plasmids that largely resemble the B. anthracis virulence plasmids. (C) 2011 Elsevier B.V. All rights reserved.

Single Molecular Conductance of Tolanes: Experimental and Theoretical Study on the Junction Evolution Dependent on the Anchoring Group

Employing a scanning tunneling microscopy based beak junction technique and mechanically controlled break junction experiments, we investigated tolane (diphenylacetylene)-type single molecular junctions having four different anchoring groups (SH, pyridyl (PY), NH2, and CN) at a solid/liquid interface. The combination of current–distance and current–voltage measurements and their quantitative statistical analysis revealed the following sequence for junction formation probability and stability: PY > SH > NH2 > CN. For all single molecular junctions investigated, we observed the evolution through multiple junction configurations, with a particularly well-defined binding geometry for PY. The comparison of density functional theory type model calculations and molecular dynamics simulations with the experimental results revealed structure and mechanistic details of the evolution of the different types of (single) molecular junctions upon stretching quantitatively.

Molecular phylogenetic analyses identify Alpine differentiation and dysploid chromosome number changes as major forces for the evolution of the European endemic Phyteuma (Campanulaceae)

Phyteuma is a chromosomally and ecologically diverse vascular plant genus and constitutes an excellent system for studying both the role of chromosomal change for species diversification and the evolution of high-mountain biota. This kind of research is, however, hampered by the lack of a sound phylogenetic framework exacerbated by the notoriously low predictive power of traditional taxonomy with respect to phylogenetic relationships in Campanulaceae. Based on a comprehensive taxon sampling and analyses of nuclear and plastid sequence and AFLP fingerprint data, Phyteuma is confirmed as a monophyletic group sister to the monotypic Physoplexis, which is in line with their peculiar flower morphologies. Within Phyteuma two clades, largely corresponding to previously recognized sections, are consistently found. The traditional circumscription of taxonomic series is largely rejected. Whereas distinctness of the currently recognized species is mostly corroborated, some interspecific relationships remain ambiguous due to incongruences between nuclear and plastid data. Major forces for diversification and evolution of Phyteuma are descending dysploidy (i.e., a decrease in chromosome base number) as well as allopatric and ecological differentiation within the Alps, the genus' center of species diversity. (C) 2013 Elsevier Inc. All rights reserved.

Adaptive Molecular Evolution in the Opsin Genes of Rapidly Speciating Cichlid Species

Cichlid fish inhabit a diverse range of environments that vary in the spectral content of light available for vision. These differences should result in adaptive selective pressure on the genes involved in visual sensitivity, the opsin genes. This study examines the evidence for differential adaptive molecular evolution in East African cichlid opsin genes due to gross differences in environmental light conditions. First, we characterize the selective regime experienced by cichlid opsin genes using a likelihood ratio test format, comparing likelihood models with different constraints on the relative rates of amino acid substitution, across sites. Second, we compare turbid and clear lineages to determine if there is evidence of differences in relative rates of substitution. Third, we present evidence of functional diversification and its relationship to the photic environment among cichlid opsin genes. We report statistical evidence of positive selection in all cichlid opsin genes, except short wavelength–sensitive 1 and short wavelength–sensitive 2b. In all genes predicted to be under positive selection, except short wavelength–sensitive 2a, we find differences in selective pressure between turbid and clear lineages. Potential spectral tuning sites are variable among all cichlid opsin genes; however, patterns of substitution consistent with photic environment–driven evolution of opsin genes are observed only for short wavelength–sensitive 1 opsin genes. This study identifies a number of promising candidate-tuning sites for future study by site-directed mutagenesis. This work also begins to demonstrate the molecular evolutionary dynamics of cichlid visual sensitivity and its relationship to the photic environment.

The evolution of antennal courtship in diplazontine parasitoid wasps (Hymenoptera, Ichneumonidae, Diplazontinae)

Background As predicted by theory, traits associated with reproduction often evolve at a comparatively high speed. This is especially the case for courtship behaviour which plays a central role in reproductive isolation. On the other hand, courtship behavioural traits often involve morphological and behavioural adaptations in both sexes; this suggests that their evolution might be under severe constraints, for instance irreversibility of character loss. Here, we use a recently proposed method to retrieve data on a peculiar courtship behavioural trait, i.e. antennal coiling, for 56 species of diplazontine parasitoid wasps. On the basis of a well-resolved phylogeny, we reconstruct the evolutionary history of antennal coiling and associated morphological modifications to study the mode of evolution of this complex character system. Results Our study reveals a large variation in shape, location and ultra-structure of male-specific modifications on the antennae. As for antennal coiling, we find either single-coiling, double-coiling or the absence of coiling; each state is present in multiple genera. Using a model comparison approach, we show that the possession of antennal modifications is highly correlated with antennal coiling behaviour. Ancestral state reconstruction shows that both antennal modifications and antennal coiling are highly congruent with the molecular phylogeny, implying low levels of homoplasy and a comparatively low speed of evolution. Antennal coiling is lost on two independent occasions, and never reacquired. A zero rate of regaining antennal coiling is supported by maximum parsimony, maximum likelihood and Bayesian approaches. Conclusions Our study provides the first comparative evidence for a tight correlation between male-specific antennal modifications and the use of the antennae during courtship. Antennal coiling in Diplazontinae evolved at a comparatively low rate, and was never reacquired in any of the studied taxa. This suggests that the loss of antennal coiling is irreversible on the timescale examined here, and therefore that evolutionary constraints have greatly influenced the evolution of antennal courtship in this group of parasitoid wasps. Further studies are needed to ascertain whether the loss of antennal coiling is irreversible on larger timescales, and whether evolutionary constraints have influenced courtship behavioural traits in a similar way in other groups.

"Pseudo-Beijing": evidence for convergent evolution in the direct repeat region of Mycobacterium tuberculosis

Background Mycobacterium tuberculosis has a global population structure consisting of six main phylogenetic lineages associated with specific geographic regions and human populations. One particular M. tuberculosis genotype known as “Beijing” has repeatedly been associated with drug resistance and has been emerging in some parts of the world. “Beijing” strains are traditionally defined based on a characteristic spoligotyping pattern. We used three alternative genotyping techniques to revisit the phylogenetic classification of M. tuberculosis complex (MTBC) strains exhibiting the typical “Beijing” spoligotyping pattern. Methods and Findings MTBC strains were obtained from an ongoing molecular epidemiological study in Switzerland and Nepal. MTBC genotyping was performed based on SNPs, genomic deletions, and 24-loci MIRU-VNTR. We identified three MTBC strains from patients originating from Tibet, Portugal and Nepal which exhibited a spoligotyping patterns identical to the classical Beijing signature. However, based on three alternative molecular markers, these strains were assigned to Lineage 3 (also known as Delhi/CAS) rather than to Lineage 2 (also known as East-Asian lineage). Sequencing of the RD207 in one of these strains showed that the deletion responsible for this “Pseudo-Beijing” spoligotype was about 1,000 base pairs smaller than the usual deletion of RD207 in classical “Beijing” strains, which is consistent with an evolutionarily independent deletion event in the direct repeat (DR) region of MTBC. Conclusions We provide an example of convergent evolution in the DR locus of MTBC, and highlight the limitation of using spoligotypes for strain classification. Our results indicate that a proportion of “Beijing” strains may have been misclassified in the past. Markers that are more phylogenetically robust should be used when exploring strain-specific differences in experimental or clinical phenotypes.

Sarcomatous evolution of oligodendroglioma ("oligosarcoma"): confirmatory report of an uncommon pattern of malignant progression in oligodendroglial tumors

By analogy to gliosarcoma, the neologism "oligosarcoma" is to describe an uncommon form of biphasic central nervous system tumor composed of contiguous neuroepithelial and mesenchymal elements, each of which individually meet the criteria of oligodendroglioma and sarcoma, respectively. By virtue of its distinctive genotype (codeletion 1p/19q), oligodendroglioma is a particularly inviting paradigm to test the assumption that such mixed tumors are clonally derived from a glial primary. We observed this constellation in a 41-year-old male who underwent two resection procedures for a recurring right frontal tumor at five years' interval. On imaging, both lesions were contrast-enhancing, and measured 7 cm × 7 cm × 6.8 cm and 7 cm × 6.5 cm × 4cm, respectively. Following the first operation, temozolomide monotherapy was administered. Whereas initial histology showed conventional anaplastic oligodendroglioma, the recurrence consisted mostly of a fibrosarcoma-like, fascicular neoplasm that was immunoreactive for vimentin, smooth muscle actin, S100 protein, and focally epithelial membrane antigen. In between, a subset of otherwise indistinguishable spindle cells expressed GFAP, and focally merged with residues of oligodendroglioma. Molecular testing for loss of heterozygosity confirmed codeletion of 1p/19q in both the primary tumor and the sarcomatous recurrence. Similarly, generalized immunoreactivity for the mutant R132H form of isocitrate dehydrogenase in both lesions indicated an identical mutation of the IDH1 gene. By the above standards, biologically consistent "oligosarcomas" are felt to be exceedingly rare, and possibly participate of a nosologically heterogeneous group of combined glial/mesenchymal lesions that may also include iatrogenically induced second malignancies as well as true collision tumors.

Evolution of bombesin conjugates for targeted PET imaging of tumors

Bombesin receptors are under intense investigation as molecular targets since they are overexpressed in several prevalent solid tumors. We rationally designed and synthesized a series of modified bombesin (BN) peptide analogs to study the influence of charge and spacers at the N-terminus, as well as amino acid substitutions, on both receptor binding affinity and pharmacokinetics. This enabled development of a novel (64/67)Cu-labeled BN peptide for PET imaging and targeted radiotherapy of BN receptor-positive tumors. Our results show that N-terminally positively charged peptide ligands had significantly higher affinity to human gastrin releasing peptide receptor (GRPr) than negatively charged or uncharged ligands (IC(50): 3.2±0.5 vs 26.3±3.5 vs 41.5±2.5 nM). The replacement of Nle(14) by Met, and deletion of D-Tyr(6), further resulted in 8-fold higher affinity. Contrary to significant changes to human GRPr binding, modifications at the N-terminal and at the 6(th), 11(th), and 14(th) position of BN induced only slight influences on affinity to mouse GRPr. [Cu(II)]-CPTA-[βAla(11)] BN(7-14) ([Cu(II)]-BZH7) showed the highest internalization rate into PC-3 cells with relatively slow efflux because of its subnanomolar affinity to GRPr. Interestingly, [(64/67)Cu]-BZH7 also displayed similar affinities to the other 2 human BN receptor subtypes. In vivo studies showed that [(64/67)Cu]-BZH7 had a high accumulation in PC-3 xenografts and allowed for clear-cut visualization of the tumor in PET imaging. In addition, a CPTA-glycine derivative, forming a hippurane-type spacer, enhanced kidney clearance of the radiotracer. These data indicate that the species variation of BN receptor plays an important role in screening radiolabeled BN. As well, the positive charge from the metallated complex at the N-terminal significantly increases affinity to human GRPr. Application of these observations enabled the novel ligand [(64/67)Cu]-BZH7 to clearly visualize PC-3 tumors in vivo. This study provides a strong starting point for optimizing radiopeptides for targeting carcinomas that express any of the BN receptor subtypes.

Evolution of gene expression changes in newborn rats after mechanical ventilation with reversible intubation

Mechanical ventilation (MV) is life-saving but potentially harmful for lungs of premature infants. So far, animal models dealt with the acute impact of MV on immature lungs, but less with its delayed effects. We used a newborn rodent model including non-surgical and therefore reversible intubation with moderate ventilation and hypothesized that there might be distinct gene expression patterns after a ventilation-free recovery period compared to acute effects directly after MV. Newborn rat pups were subjected to 8 hr of MV with 60% oxygen (O(2)), 24 hr after injection of lipopolysaccharide (LPS), intended to create a low inflammatory background as often recognized in preterm infants. Animals were separated in controls (CTRL), LPS injection (LPS), or full intervention with LPS and MV with 60% O(2) (LPS + MV + O(2)). Lungs were recovered either directly following (T:0 hr) or 48 hr after MV (T:48 hr). Histologically, signs of ventilator-induced lung injury (VILI) were observed in LPS + MV + O(2) lungs at T:0 hr, while changes appeared similar to those known from patients with chronic lung disease (CLD) with fewer albeit larger gas exchange units, at T:48 hr. At T:0 hr, LPS + MV + O(2) increased gene expression of pro-inflammatory MIP-2. In parallel anti-inflammatory IL-1Ra gene expression was increased in LPS and LPS + MV + O(2) groups. At T:48 hr, pro- and anti-inflammatory genes had returned to their basal expression. MMP-2 gene expression was decreased in LPS and LPS + MV + O(2) groups at T:0 hr, but no longer at T:48 hr. MMP-9 gene expression levels were unchanged directly after MV. However, at T:48 hr, gene and protein expression increased in LPS + MV + O(2) group. In conclusion, this study demonstrates the feasibility of delayed outcome measurements after a ventilation-free period in newborn rats and may help to further understand the time-course of molecular changes following MV. The differences obtained from the two time points could be interpreted as an initial transitory increase of inflammation and a delayed impact of the intervention on structure-related genes.

Selection and recombination drive the evolution of MHC class II DRB diversity in ungulates

Major histocompatibility complex (MHC) antigen-presenting genes are the most variable loci in vertebrate genomes. Host-parasite co-evolution is assumed to maintain the excessive polymorphism in the MHC loci. However, the molecular mechanisms underlying the striking diversity in the MHC remain contentious. The extent to which recombination contributes to the diversity at MHC loci in natural populations is still controversial, and there have been only few comparative studies that make quantitative estimates of recombination rates. In this study, we performed a comparative analysis for 15 different ungulates species to estimate the population recombination rate, and to quantify levels of selection. As expected for all species, we observed signatures of strong positive selection, and identified individual residues experiencing selection that were congruent with those constituting the peptide-binding region of the human DRB gene. However, in addition for each species, we also observed recombination rates that were significantly different from zero on the basis of likelihood-permutation tests, and in other non-quantitative analyses. Patterns of synonymous and non-synonymous sequence diversity were consistent with differing demographic histories between species, but recent simulation studies by other authors suggest inference of selection and recombination is likely to be robust to such deviations from standard models. If high rates of recombination are common in MHC genes of other taxa, re-evaluation of many inference-based phylogenetic analyses of MHC loci, such as estimates of the divergence time of alleles and trans-specific polymorphism, may be required.

Statistical evaluation of alternative models of human evolution

An appropriate model of recent human evolution is not only important to understand our own history, but it is necessary to disentangle the effects of demography and selection on genome diversity. Although most genetic data support the view that our species originated recently in Africa, it is still unclear if it completely replaced former members of the Homo genus, or if some interbreeding occurred during its range expansion. Several scenarios of modern human evolution have been proposed on the basis of molecular and paleontological data, but their likelihood has never been statistically assessed. Using DNA data from 50 nuclear loci sequenced in African, Asian and Native American samples, we show here by extensive simulations that a simple African replacement model with exponential growth has a higher probability (78%) as compared with alternative multiregional evolution or assimilation scenarios. A Bayesian analysis of the data under this best supported model points to an origin of our species approximately 141 thousand years ago (Kya), an exit out-of-Africa approximately 51 Kya, and a recent colonization of the Americas approximately 10.5 Kya. We also find that the African replacement model explains not only the shallow ancestry of mtDNA or Y-chromosomes but also the occurrence of deep lineages at some autosomal loci, which has been formerly interpreted as a sign of interbreeding with Homo erectus.

Parallel evolution of facial stripe patterns in the Neolamprologus brichardilpulcher species complex endemic to Lake Tanganyika

Colour pattern diversity can be due to random processes or to natural or sexual selection. Consequently, similarities in colour patterns are not always correlated with common ancestry, but may result from convergent evolution under shared selection pressures or drift. Neolamprologus brichardi and Neolamprologus pulcher have been described as two distinct species based on differences in the arrangement of two dark bars on the operculum. Our study uses DNA sequences of the mitochondrial control region to show that relatedness of haplotypes disagrees with species assignment based on head colour pattern. This suggests repeated parallel evolution of particular stripe patterns. The complete lack of shared haplotypes between populations of the same or different phenotypes reflects strong philopatric behaviour, possibly induced by the cooperative breeding mode in which offspring remain in their natal territory and serve as helpers until they disperse to nearby territories or take over a breeding position. Concordant phylogeographic patterns between N. brichardi/N. pulcher populations and other rock-dwelling cichlids suggest that the same colonization routes have been taken by sympatric species and that these routes were affected by lake level fluctuations in the past. (C) 2007 Elsevier Inc. All rights reserved.

Single-Molecule Conductance of Functionalized Oligoynes: Length Dependence and Junction Evolution

We report a combined experimental and theoretical investigation of the length dependence and anchor group dependence of the electrical conductance of a series of oligoyne molecular wires in single-molecule junctions with gold contacts. Experimentally, we focus on the synthesis and properties of diaryloligoynes with n = 1, 2, and 4 triple bonds and the anchor dihydrobenzo[b]thiophene (BT). For comparison, we also explored the aurophilic anchor group cyano (CN), amino (NH2), thiol (SH), and 4-pyridyl (PY). Scanning tunneling microscopy break junction (STM-BJ) and mechanically controllable break junction (MCBJ) techniques are employed to investigate single-molecule conductance characteristics. The BT moiety is superior as compared to traditional anchoring groups investigated so far. BT-terminated oligoynes display a 100% probability of junction formation and possess conductance values which are the highest of the oligoynes studied and, moreover, are higher than other conjugated molecular wires of similar length. Density functional theory (DFT)-based calculations are reported for oligoynes with n = 1−4 triple bonds. Complete conductance traces and conductance distributions are computed for each family of molecules. The sliding of the anchor groups leads to oscillations in both the electrical conductance and the binding energies of the studied molecular wires. In agreement with experimental results, BT-terminated oligoynes are predicted to have a high electrical conductance. The experimental attenuation constants βH range between 1.7 nm−1 (CN) and 3.2 nm−1 (SH) and show the following trend: βH(CN) < βH(NH2) < βH(BT) < βH(PY) ≈ βH(SH). DFT-based calculations yield lower values, which range between 0.4 nm−1 (CN) and 2.2 nm−1 (PY).