Francisella tularensis is the etiological agent of tularemia, a zoonotic disease with worldwide prevalence. F. tularensis is a highly pathogenic organism and has been designated as a potential biothreat agent. Currently there are four recognized subspecies of F. tularensis: tularensis (type A), holarctica (type B), mediasiatica, and novicida. In addition, genomic studies have further subdivided type A tularensis into two subclassifications, type A.I and type A.II. These two subclassifications differ in geographic distribution with type A.I appearing mainly in the Eastern United States and type A.II appearing mainly in the Western United States. Because of differences of virulence among the subspecies, it is important to be able to quickly identify each of the subspecies rapidly and accurately. This work describes the development of a multiplex real-time polymerase chain reaction (PCR) assay which was shown to be ~98% successful at identifying the known subspecies of F. tularensis. Furthermore, F. tularensis is thought be a genome in decay (losing genes) because of the relatively large number of pseudogenes present in its genome. We hypothesized that the observed frequency of gene loss/pseudogenes may be an artifact of evolution in response to a changing environment, and that genes involved in virulence should be under strong positive selection. Eleven arbitrarily chosen virulence genes were screened for positive selection along with 10 arbitrarily chosen housekeeping genes. Analyses of selection yielded one housekeeping gene and 7 virulence genes which showed significant evidence of positive selection. Our results suggest that while the loss of functional genes through disuse could be accelerated by negative selection, the genome decay in Francisella could also be the byproduct of adaptive evolution, as evidenced by several of its virulence genes which are undergoing strong, positive selection.

The cytochrome bc1 complex of the mitochondrial electron transport chain accomplishes the enzymatic reaction known as the modified Q-cycle. In the Q-cycle the bc1 complex transports protons from the matrix to the intermembrane space of the mitochondria, creating the proton gradient used to make ATP. The energy to move these protons is obtained by shuttling electrons from the coenzyme ubiquinol (QH2) to coenzyme ubiquinone (Q) and the mobile cytochrome c. This well studied complex is ideal for examining molecular adaptation because it consists of ten different subunits, it functions as a dimer, and it includes at least five different active sites. The program TreeSAAP was used to characterize molecular adaptation in the bc1 complex and identify specific amino acid sites that experienced positive destabilizing (radical) selection. Using this information and three-dimensional structures of the protein complex, selection was characterized in terms of coevolution and coadaptation. Coevolution is described as reciprocal local biochemical shifts based on phylogenetic location and results in overall maintenance. Coadaptation, on the other hand, is more dynamic and is described as coordinated local biochemical shifts based on phylogenetic location which results in overall adaptation. In this study both coevolution and coadaptation were identified in various locations on the protein complex near the active sites. Sites in the pore region of cyt c1 were shown to exhibit coevolution, in other words maintenance, of many biochemical properties, whereas sites on helix H of cyt b, which flanks the active sites Qo and Qi, were shown to exhibit coadaptation, in other words coordinated shifts in the specific properties equilibrium constant and solvent accessible reduction ratio. Also, different domains of the protein exhibited significant shifts in drastically different amino acid properties: the protein imbedded in the membrane demonstrated shifts in mainly functional properties, while the part of the complex in the intermembrane space demonstrated shifts in conformational, structural, and energetic properties.