K-r/K-c but not d(N)/d(S) correlates positively with body mass in birds, raising implications for inferring lineage-specific selection

Background: The ratio of the rates of non-synonymous and synonymous substitution (d(N)/d(S)) is commonly used to estimate selection in coding sequences. It is often suggested that, all else being equal, d(N)/d(S) should be lower in populations with large effective size (Ne) due to increased efficacy of purifying selection. As N-e is difficult to measure directly, life history traits such as body mass, which is typically negatively associated with population size, have commonly been used as proxies in empirical tests of this hypothesis. However, evidence of whether the expected positive correlation between body mass and d(N)/d(S) is consistently observed is conflicting. Results: Employing whole genome sequence data from 48 avian species, we assess the relationship between rates of molecular evolution and life history in birds. We find a negative correlation between dN/dS and body mass, contrary to nearly neutral expectation. This raises the question whether the correlation might be a method artefact. We therefore in turn consider non-stationary base composition, divergence time and saturation as possible explanations, but find no clear patterns. However, in striking contrast to d(N)/d(S), the ratio of radical to conservative amino acid substitutions (K-r/K-c) correlates positively with body mass. Conclusions: Our results in principle accord with the notion that non-synonymous substitutions causing radical amino acid changes are more efficiently removed by selection in large populations, consistent with nearly neutral theory. These findings have implications for the use of d(N)/d(S) and suggest that caution is warranted when drawing conclusions about lineage-specific modes of protein evolution using this metric.

Functional characterization of the Pneumocystis jirovecii potential drug targets dhfs and abz2 involved in folate biosynthesis.

Pneumocystis species are fungal parasites colonizing mammal lungs with strict host specificity. Pneumocystis jirovecii is the human-specific species and can turn into an opportunistic pathogen causing severe pneumonia in immunocompromised individuals. This disease is currently the second most frequent life-threatening invasive fungal infection worldwide. The most efficient drug, cotrimoxazole, presents serious side effects, and resistance to this drug is emerging. The search for new targets for the development of new drugs is thus of utmost importance. The recent release of the P. jirovecii genome sequence opens a new era for this task. It can now be carried out on the actual targets to be inhibited instead of on those of the relatively distant model Pneumocystis carinii, the species infecting rats. We focused on the folic acid biosynthesis pathway because (i) it is widely used for efficient therapeutic intervention, and (ii) it involves several enzymes that are essential for the pathogen and have no human counterparts. In this study, we report the identification of two such potential targets within the genome of P. jirovecii, the dihydrofolate synthase (dhfs) and the aminodeoxychorismate lyase (abz2). The function of these enzymes was demonstrated by the rescue of the null allele of the orthologous gene of Saccharomyces cerevisiae.

Genomic and transcriptomic analysis of the streptomycin-dependent Mycobacterium tuberculosis strain 18b.

The ability of Mycobacterium tuberculosis to establish a latent infection (LTBI) in humans confounds the treatment of tuberculosis. Consequently, there is a need to discover new therapeutic agents that can kill M. tuberculosis both during active disease and LTBI. The streptomycin-dependent strain of M. tuberculosis, 18b, provides a useful tool for this purpose since upon removal of streptomycin (STR) it enters a non-replicating state that mimics latency both in vitro and in animal models. The 4.41 Mb genome sequence of M. tuberculosis 18b was determined and this revealed the strain to belong to clade 3 of the ancient ancestral lineage of the Beijing family. STR-dependence was attributable to insertion of a single cytosine in the 530 loop of the 16S rRNA and to a single amino acid insertion in the N-terminal domain of initiation factor 3. RNA-seq was used to understand the genetic programme activated upon STR-withdrawal and hence to gain insight into LTBI. This revealed reconfiguration of gene expression and metabolic pathways showing strong similarities between non-replicating 18b and M. tuberculosis residing within macrophages, and with the core stationary phase and microaerophilic responses. The findings of this investigation confirm the validity of 18b as a model for LTBI, and provide insight into both the evolution of tubercle bacilli and the functioning of the ribosome.

Variations of plasmid content in Rickettsia felis.

BACKGROUND: Since its first detection, characterization of R. felis has been a matter of debate, mostly due to the contamination of an initial R. felis culture by R. typhi. However, the first stable culture of R. felis allowed its precise phenotypic and genotypic characterization, and demonstrated that this species belonged to the spotted fever group rickettsiae. Later, its genome sequence revealed the presence of two forms of the same plasmid, physically confirmed by biological data. In a recent article, Gillespie et al. (PLoS One. 2007;2(3):e266.) used a bioinformatic approach to refute the presence of the second plasmid form, and proposed the creation of a specific phylogenetic group for R. felis. METHODOLOGY/PRINCIPAL FINDINGS: In the present report, we, and five independent international laboratories confirmed unambiguously by PCR the presence of two plasmid forms in R. felis strain URRWXCal(2) (T), but observed that the plasmid content of this species, from none to 2 plasmid forms, may depend on the culture passage history of the studied strain. We also demonstrated that R. felis does not cultivate in Vero cells at 37 degrees C but generates plaques at 30 degrees C. Finally, using a phylogenetic study based on 667 concatenated core genes, we demonstrated the position of R. felis within the spotted fever group. SIGNIFICANCE: We demonstrated that R. felis, which unambiguously belongs to the spotted fever group rickettsiae, may contain up to two plasmid forms but this plasmid content is unstable.

Criblamydia sequanensis Harbors a Megaplasmid Encoding Arsenite Resistance.

Criblamydia sequanensis is an amoeba-resisting bacterium recently isolated from the Seine River. This Chlamydia-related bacterium harbors a genome of approximately 3 Mbp and a megaplasmid of 89,525 bp. The plasmid encodes several efflux systems and an operon for arsenite resistance. This first genome sequence within the Criblamydiaceae family enlarges our view on the evolution and the ecology of this important bacterial clade largely understudied so far.

Partial characterization of nif genes from the bacterium Azospirillum amazonense

Azospirillum amazonense revealed genomic organization patterns of the nitrogen fixation genes similar to those of the distantly related species A. brasilense. Our work suggests that A. brasilense nifHDK, nifENX, fixABC operons and nifA and glnB genes may be structurally homologous to the counterpart genes of A. amazonense. This is the first analysis revealing homology between A. brasilense nif genes and the A. amazonense genome. Sequence analysis of PCR amplification products revealed similarities between the amino acid sequences of the highly conserved nifD and glnB genes of A. amazonense and related genes of A. brasilense and other bacteria. However, the A. amazonense non-coding regions (the upstream activator sequence region and the region between the nifH and nifD genes) differed from related regions of A. brasilense even in nitrogenase structural genes which are highly conserved among diazotrophic bacteria. The feasibility of the 16S ribosomal RNA gene-based PCR system for specific detection of A. amazonense was shown. Our results indicate that the PCR primers for 16S rDNA defined in this article are highly specific to A. amazonense and can distinguish this species from A. brasilense.

Genomics and X-ray microanalysis indicate that Ca2+ and thiols mediate the aggregation and adhesion of Xylella fastidiosa

The availability of the genome sequence of the bacterial plant pathogen Xylella fastidiosa, the causal agent of citrus variegated chlorosis, is accelerating important investigations concerning its pathogenicity. Plant vessel occlusion is critical for symptom development. The objective of the present study was to search for information that would help to explain the adhesion of X. fastidiosa cells to the xylem. Scanning electron microscopy revealed that adhesion may occur without the fastidium gum, an exopolysaccharide produced by X. fastidiosa, and X-ray microanalysis demonstrated the presence of elemental sulfur both in cells grown in vitro and in cells found inside plant vessels, indicating that the sulfur signal is generated by the pathogen surface. Calcium and magnesium peaks were detected in association with sulfur in occluded vessels. We propose an explanation for the adhesion and aggregation process. Thiol groups, maintained by the enzyme peptide methionine sulfoxide reductase, could be active on the surface of the bacteria and appear to promote cell-cell aggregation by forming disulfide bonds with thiol groups on the surface of adjacent cells. The enzyme methionine sulfoxide reductase has been shown to be an auxiliary component in the adhesiveness of some human pathogens. The negative charge conferred by the ionized thiol group could of itself constitute a mechanism of adhesion by allowing the formation of divalent cation bridges between the negatively charged bacteria and predominantly negatively charged xylem walls.

Coping with genetic diversity: the contribution of pathogen and human genomics to modern vaccinology

Vaccine development faces major difficulties partly because of genetic variation in both infectious organisms and humans. This causes antigenic variation in infectious agents and a high interindividual variability in the human response to the vaccine. The exponential growth of genome sequence information has induced a shift from conventional culture-based to genome-based vaccinology, and allows the tackling of challenges in vaccine development due to pathogen genetic variability. Additionally, recent advances in immunogenetics and genomics should help in the understanding of the influence of genetic factors on the interindividual and interpopulation variations in immune responses to vaccines, and could be useful for developing new vaccine strategies. Accumulating results provide evidence for the existence of a number of genes involved in protective immune responses that are induced either by natural infections or vaccines. Variation in immune responses could be viewed as the result of a perturbation of gene networks; this should help in understanding how a particular polymorphism or a combination thereof could affect protective immune responses. Here we will present: i) the first genome-based vaccines that served as proof of concept, and that provided new critical insights into vaccine development strategies; ii) an overview of genetic predisposition in infectious diseases and genetic control in responses to vaccines; iii) population genetic differences that are a rationale behind group-targeted vaccines; iv) an outlook for genetic control in infectious diseases, with special emphasis on the concept of molecular networks that will provide a structure to the huge amount of genomic data.

The patentability of human genetic material in China : a comparative analysis

"Mémoire Présenté à la Faculté des Études Supérieures en vue de l'obtention du Grade de Maîtrise En Droit Option Recherche"

The evolution of inter-genomic variation in arbuscular mycorrhizal fungi

Contexte: Les champignons mycorhiziens à arbuscules (AMF) établissent des relations symbiotiques avec la plupart des plantes grâce à leurs réseaux d’hyphes qui s’associent avec les racines de leurs hôtes. De précédentes études ont révélé des niveaux de variation génétique extrêmes pour des loci spécifiques permettant de supposer que les AMF peuvent contenir des milliers de noyaux génétiquement divergents dans un même cytoplasme. Si aucun processus de reproduction sexuée n’a jusqu’ici été observé chez ces mycorhizes, on constate cependant que des niveaux élevés de variation génétique peuvent être maintenus à la fois par l’échange de noyaux entre hyphes et par des processus fréquents de recombinaison entre noyaux. Les AMF se propagent par l’intermédiaire de spores qui contiennent chacune un échantillon d’une population initiale de noyaux hétérogènes, directement hérités du mycélium parent. À notre connaissance les AMF sont les seuls organismes qui ne passent jamais par un stade mononucléaire, ce qui permet aux noyaux de diverger génétiquement dans un même cytoplasme. Ces aspects singuliers de la biologie des AMF rendent l’estimation de leur diversité génétique problématique. Ceci constitue un défi majeur pour les écologistes sur le terrain mais également pour les biologistes moléculaires dans leur laboratoire. Au-delà même des problématiques de diversité spécifique, l’amplitude du polymorphisme entre noyaux mycorhiziens est mal connue. Le travail proposé dans ce manuscrit de thèse explore donc les différents aspects de l’architecture génomique singulière des AMF. Résultats L’ampleur du polymorphisme intra-isolat a été déjà observée pour la grande sous-unité d’ARN ribosomal de l’isolat Glomus irregulare DAOM-197198 (précédemment identifié comme G. intraradices) et pour le gène de la polymerase1-like (PLS) de Glomus etunicatum isolat NPI. Dans un premier temps, nous avons pu confirmer ces résultats et nous avons également pu constater que ces variations étaient transcrites. Nous avons ensuite pu mettre en évidence la présence d’un goulot d’étranglement génétique au moment de la sporulation pour le locus PLS chez l’espèce G. etunicatum illustrant les importants effets d’échantillonnage qui se produisaient entre chaque génération de spore. Enfin, nous avons estimé la différentiation génétique des AMF en utilisant à la fois les réseaux de gènes appliqués aux données de séquençage haut-débit ainsi que cinq nouveaux marqueurs génomiques en copie unique. Ces analyses révèlent que la différenciation génomique est présente de manière systématique dans deux espèces (G. irregulare et G. diaphanum). Conclusions Les résultats de cette thèse fournissent des preuves supplémentaires en faveur du scénario d’une différenciation génomique entre noyaux au sein du même isolat mycorhizien. Ainsi, au moins trois membres du genre Glomus, G. irregulare, G. diaphanum and G. etunicatum, apparaissent comme des organismes dont l’organisation des génomes ne peut pas être décrit d’après un modèle Mendélien strict, ce qui corrobore l’hypothèse que les noyaux mycorhiziens génétiquement différenciés forment un pangenome.

Plasmodium vivax lineages: geographical distribution, tandem repeat polymorphism, and phylogenetic relationship

Background: Multi-drug resistance and severe/ complicated cases are the emerging phenotypes of vivax malaria, which may deteriorate current anti-malarial control measures. The emergence of these phenotypes could be associated with either of the two Plasmodium vivax lineages. The two lineages had been categorized as Old World and New World, based on geographical sub-division and genetic and phenotypical markers. This study revisited the lineage hypothesis of P. vivax by typing the distribution of lineages among global isolates and evaluated their genetic relatedness using a panel of new mini-satellite markers. Methods: 18S SSU rRNA S-type gene was amplified from 420 Plasmodium vivax field isolates collected from different geographical regions of India, Thailand and Colombia as well as four strains each of P. vivax originating from Nicaragua, Panama, Thailand (Pak Chang), and Vietnam (ONG). A mini-satellite marker panel was then developed to understand the population genetic parameters and tested on a sample subset of both lineages. Results: 18S SSU rRNA S-type gene typing revealed the distribution of both lineages (Old World and New World) in all geographical regions. However, distribution of Plasmodium vivax lineages was highly variable in every geographical region. The lack of geographical sub-division between lineages suggests that both lineages are globally distributed. Ten mini-satellites were scanned from the P. vivax genome sequence; these tandem repeats were located in eight of the chromosomes. Mini-satellites revealed substantial allelic diversity (7-21, AE = 14.6 +/- 2.0) and heterozygosity (He = 0.697-0.924, AE = 0.857 +/- 0.033) per locus. Mini-satellite comparison between the two lineages revealed high but similar pattern of genetic diversity, allele frequency, and high degree of allele sharing. A Neighbour-Joining phylogenetic tree derived from genetic distance data obtained from ten mini-satellites also placed both lineages together in every cluster. Conclusions: The global lineage distribution, lack of genetic distance, similar pattern of genetic diversity, and allele sharing strongly suggested that both lineages are a single species and thus new emerging phenotypes associated with vivax malaria could not be clearly classified as belonging to a particular lineage on basis of their geographical origin.

A genomewide survey of developmentally relevant genes in Ciona intestinalis - III. Genes for Fox, ETS, nuclear receptors and NF kappa B

A survey against the draft genome sequence and the cDNA/EST database of Ciona intestinalis identified a number of genes encoding transcription factors regulating a variety of processes including development. In the present study, we describe almost complete sets of genes for Fox, ETS-domain transcription factors, nuclear receptors, and NFkappaB as well as other factors regulating NFkappaB activity, with their phylogenetic nature. Vertebrate Fox transcription factors are currently delineated into 17 subfamilies: FoxA to FoxQ. The present survey yielded 29 genes of this family in the Ciona genome, 24 of which were Ciona orthologues of known Fox genes. In addition, we found 15 ETS aenes, 17 nuclear receptor genes, and several NFkappaB signaling pathway genes in the Ciona genome. The number of Ciona genes in each family is much smaller than that of vertebrates, which represents a simplified feature of the ascidian genome. For example, humans have two NFkappaB genes, three Rel genes, and five NFAT genes, while Ciona has one gene for each family. The Ciona genome also contains smaller numbers of genes for the NFkappaB regulatory system, i.e. after the split of ascidians/vertebrates, vertebrates evolved a more complex NFkappaB system. The present results therefore provide molecular information for the investigation of complex developmental processes, and an insight into chordate evolution.

Physiological changes in Campylobacter jejuni on entry into stationary phase

Campylobacter jejuni NCTC 11168 does not exhibit the general increase in cellular stress resistance on entry into stationary phase that is seen in most other bacteria. This is consistent with the lack of global stationary phase regulatory elements in this organism. deduced from an analysis of its genome sequence. We now show that C. jejuni NCTC 11168 does undergo certain changes in stationary phase, of a pattern not previously described. As cells entered stationary phase there was a change in membrane fatty acid composition, principally a decrease in the proportion of unsaturated fatty acids and an increase in the content of cyclopropane and short-chain fatty acids. These changes in membrane composition were accompanied by an increase in the resilience of the cell membrane towards loss of integrity caused by pressure and an increase in cellular pressure resistance. By contrast. there were no major changes in resistance to acid or heat treatment. A similar pattern of changes in stress resistance on entry, into stationary phase was seen in C. jejuni NCTC 11351, the type strain. These changes appear to represent a restricted Physiological response to the conditions existing in stationary phase cultures, in an organism having limited capacity for genetic regulation and adaptation to environment. © 2004 Elsevier B.V. All rights reserved.

Type III secretion homologs are present in Brucella melitensis, B-ovis, and B-suis biovars 1, 2, and 3

Protein sequences from characterized type III secretion (TTS) systems were used as probes in silico to identify several TTS gene homologs in the genome sequence of Brucella suis biovar 1 strain 1330. Four of the genes, named flhB, fliP, fliR, and fliF on the basis of greatest homologies to known flagellar apparatus proteins, were targeted in PCR and hybridization assays to determine their distribution among other Brucella nomen species and biovars. The results indicated that flhB, fliP, fliR and fliF are present in Brucella melitensis, Brucella ovis, and Brucella suis biovars 1, 2 and 3. Similar homologos have been reported previously in Brucella abortus. Using RT-PCR assays, we were unable to detect any expression of these genes. It is not yet known whether the genes are the cryptic remnants of a flagellar system or are actively involved in a process contributing to pathogenicity or previously undetected motility, but they are distributed widely in Brucella and merit further study to determine their role.

Regulatory hotspots are associated with plant gene expression under varying soil phosphorus supply in Brassica rapa

Gene expression is a quantitative trait that can be mapped genetically in structured populations to identify expression quantitative trait loci (eQTL). Genes and regulatory networks underlying complex traits can subsequently be inferred. Using a recently released genome sequence, we have defined cis- and trans-eQTL and their environmental response to low phosphorus (P) availability within a complex plant genome and found hotspots of trans-eQTL within the genome. Interval mapping, using P supply as a covariate, revealed 18,876 eQTL. trans-eQTL hotspots occurred on chromosomes A06 and A01 within Brassica rapa; these were enriched with P metabolism-related Gene Ontology terms (A06) as well as chloroplast-and photosynthesis-related terms (A01). We have also attributed heritability components to measures of gene expression across environments, allowing the identification of novel gene expression markers and gene expression changes associated with low P availability. Informative gene expression markers were used to map eQTL and P use efficiency-related QTL. Genes responsive to P supply had large environmental and heritable variance components. Regulatory loci and genes associated with P use efficiency identified through eQTL analysis are potential targets for further characterization and may have potential for crop improvement.