Large-scale isolation of candidate virulence genes of Pseudomonas aeruginosa by in vivo selection.

Pseudomonas aeruginosa, an opportunistic human pathogen, is a major causative agent of mortality and morbidity in immunocompromised patients and those with cystic fibrosis genetic disease. To identify new virulence genes of P. aeruginosa, a selection system was developed based on the in vivo expression technology (IVET) that was first reported in Salmonella system. An adenine-requiring auxotrophic mutant strain of P. aeruginosa was isolated and found avirulent on neutropenic mice. A DNA fragment that can complement the mutant strain, containing purEK operon that is required for de novo biosynthesis of purine, was sequenced and used in the IVET vector construction. By applying the IVET selection system to a neutropenic mouse infection model, genetic loci that are specifically induced in vivo were identified. Twenty-two such loci were partially sequenced and analyzed. One of them was a well-studied virulence factor, pyochelin receptor (FptA), that is involved in iron acquisition. Fifteen showed significant homology to reported sequences in GenBank, while the remaining six did not. One locus, designated np20, encodes an open reading frame that shares amino acid sequence homology to transcriptional regulators, especially to the ferric uptake regulator (Fur) proteins of other bacteria. An insertional np20 null mutant strain of P. aeruginosa did not show a growth defect on laboratory media; however, its virulence on neutropenic mice was significantly reduced compared with that of a wild-type parent strain, demonstrating the importance of the np20 locus in the bacterial virulence. The successful isolation of genetic loci that affect bacterial virulence demonstrates the utility of the IVET system in identification of new virulence genes of P. aeruginosa.

Computer aided selection of candidate vaccine antigens

Immunoinformatics is an emergent branch of informatics science that long ago pullulated from the tree of knowledge that is bioinformatics. It is a discipline which applies informatic techniques to problems of the immune system. To a great extent, immunoinformatics is typified by epitope prediction methods. It has found disappointingly limited use in the design and discovery of new vaccines, which is an area where proper computational support is generally lacking. Most extant vaccines are not based around isolated epitopes but rather correspond to chemically-treated or attenuated whole pathogens or correspond to individual proteins extract from whole pathogens or correspond to complex carbohydrate. In this chapter we attempt to review what progress there has been in an as-yet-underexplored area of immunoinformatics: the computational discovery of whole protein antigens. The effective development of antigen prediction methods would significantly reduce the laboratory resource required to identify pathogenic proteins as candidate subunit vaccines. We begin our review by placing antigen prediction firmly into context, exploring the role of reverse vaccinology in the design and discovery of vaccines. We also highlight several competing yet ultimately complementary methodological approaches: sub-cellular location prediction, identifying antigens using sequence similarity, and the use of sophisticated statistical approaches for predicting the probability of antigen characteristics. We end by exploring how a systems immunomics approach to the prediction of immunogenicity would prove helpful in the prediction of antigens.

Selection of candidate probiotic strains from human microbial niches, for the development of tailored foods designed for specific consumers

My PhD project was intended, throughout the selection of probiotics from human milk and healthy vaginal environment, for the development of tailored fermented foods. According to this aim, several activities were carried out. The first one, concerning the isolation of Lactobacillus and Bifidobacterium strains from human milk to find new probiotic candidates to be included in food products showed promising results. Probiotics have been also proposed to improve female genital health and microbial strains isolated and connected with healthy vaginal ecosystem could be used to prevent or treat vaginal dysbiosis. In this context vaginal lactobacilli previously characterized for their technological features and antagonistic activity against several female uro-genital pathogens were investigated for their metabolic aptitude and additional probiotic features, showing interesting results hypothesizing their inclusion in foods. In addition, in order to preserve vaginal strains viability during food processing/digestion it was also evaluated the potential of microencapsulation by spray-drying. In this framework the results obtained were highly promising from the perspective of using encapsulated powders in food formulations. Another activity connected with the main idea to develop a food strategy for the administration of these vaginal strains was carried out. Lactobacillus crispatus BC4, was supplemented in a Squacquerone cheese, and its digestive fate was evaluated adopting SHIME® system. The results showed that during colonic fermentation, L. crispatus BC4 was metabolically active. Additionally, although probiotic delivery to humans has traditionally been associated with fermented dairy foods, recently the demand for non-dairy-alternatives as potential probiotics carrier is increasing. In this framework, my latest activity was connected with the development of fermented soy milks with encapsulated and non-encapsulated L. crispatus BC4 and L. gasseri BC9. The same fermented soy milks were also investigated for their nutritional qualities and after in vitro digestion for their specific functionality on post-menopausal fecal microbiota and protein bioaccessibility.

Selection And Validation Of Reference Genes For Functional Studies In The Calliphoridae Family.

The genera Cochliomyia and Chrysomya contain both obligate and saprophagous flies, which allows the comparison of different feeding habits between closely related species. Among the different strategies for comparing these habits is the use of qPCR to investigate the expression levels of candidate genes involved in feeding behavior. To ensure an accurate measure of the levels of gene expression, it is necessary to normalize the amount of the target gene with the amount of a reference gene having a stable expression across the compared species. Since there is no universal gene that can be used as a reference in functional studies, candidate genes for qPCR data normalization were selected and validated in three Calliphoridae (Diptera) species, Cochliomyia hominivorax Coquerel, Cochliomyia macellaria Fabricius, and Chrysomya albiceps Wiedemann . The expression stability of six genes ( Actin, Gapdh, Rp49, Rps17, α -tubulin, and GstD1) was evaluated among species within the same life stage and between life stages within each species. The expression levels of Actin, Gapdh, and Rp49 were the most stable among the selected genes. These genes can be used as reliable reference genes for functional studies in Calliphoridae using similar experimental settings.

Selection of suitable housekeeping genes for expression analysis in glioblastoma using quantitative RT-PCR

Background: Considering the broad variation in the expression of housekeeping genes among tissues and experimental situations, studies using quantitative RT-PCR require strict definition of adequate endogenous controls. For glioblastoma, the most common type of tumor in the central nervous system, there was no previous report regarding this issue. Results: Here we show that amongst seven frequently used housekeeping genes TBP and HPRT1 are adequate references for glioblastoma gene expression analysis. Evaluation of the expression levels of 12 target genes utilizing different endogenous controls revealed that the normalization method applied might introduce errors in the estimation of relative quantities. Genes presenting expression levels which do not significantly differ between tumor and normal tissues can be considered either increased or decreased if unsuitable reference genes are applied. Most importantly, genes showing significant differences in expression levels between tumor and normal tissues can be missed. We also demonstrated that the Holliday Junction Recognizing Protein, a novel DNA repair protein over expressed in lung cancer, is extremely over-expressed in glioblastoma, with a median change of about 134 fold. Conclusion: Altogether, our data show the relevance of previous validation of candidate control genes for each experimental model and indicate TBP plus HPRT1 as suitable references for studies on glioblastoma gene expression.

Ecological criteria for evaluating candidate sites for marine reserves

Several schemes have been developed to help select the locations of marine reserves. All of them combine social, economic, and biological criteria, and few offer any guidance as to how to prioritize among the criteria identified. This can imply that the relative weights given to different criteria are unimportant. Where two sites are of equal value ecologically; then socioeconomic criteria should dominate the choice of which should be protected. However, in many cases, socioeconomic criteria are given equal or greater weight than ecological considerations in the choice of sites. This can lead to selection of reserves with little biological value that fail to meet many of the desired objectives. To avoid such a possibility, we develop a series of criteria that allow preliminary evaluation of candidate sites according to their relative biological values in advance of the application of socioeconomic criteria. We include criteria that,. while not strictly biological, have a strong influence on the species present or ecological processes. Out scheme enables sites to be assessed according to their biodiversity, the processes which underpin that diversity, and the processes that support fisheries and provide a spectrum of other services important to people. Criteria that capture biodiversity values include biogeographic representation, habitat representation and heterogeneity, and presence of species or populations of special interest (e.g., threatened species). Criteria that capture sustainability of biodiversity and fishery values include the size of reserves necessary to protect viable habitats, presence of exploitable species, vulnerable life stages, connectivity among reserves, links among ecosystems, and provision of ecosystem services to people. Criteria measuring human and natural threats enable candidate sites to be eliminated from consideration if risks are too great, but also help prioritize among sites where threats can be mitigated by protection. While our criteria can be applied to the design of reserve networks, they also enable choice of single reserves to be made in the context of the attributes of existing protected areas. The overall goal of our scheme is to promote the development of reserve networks that will maintain biodiversity and ecosystem functioning at large scales. The values of eco-system goods and services for people ultimately depend on meeting this objective.

The genetic basis of morphological change in convergent evolution of natural populations: identifying candidate genes behind convergent evolution in blind cavefish astyanax mexicanus

Dissertation presented to obtain the Ph.D. degree in Biology at the Instituto de Tecnologia Química e Biológica, Universidade Nova de Lisboa.

Molecular variation at a candidate gene implicated in the regulation of fire ant social behavior

The fire ant Solenopsis invicta and its close relatives display an important social polymorphism involving differences in colony queen number. Colonies are headed by either a single reproductive queen (monogyne form) or multiple queens (polygyne form). This variation in social organization is associated with variation at the gene Gp-9, with monogyne colonies harboring only B-like allelic variants and polygyne colonies always containing b-like variants as well. We describe naturally occurring variation at Gp-9 in fire ants based on 185 full-length sequences, 136 of which were obtained from S. invicta collected over much of its native range. While there is little overall differentiation between most of the numerous alleles observed, a surprising amount is found in the coding regions of the gene, with such substitutions usually causing amino acid replacements. This elevated coding-region variation may result from a lack of negative selection acting to constrain amino acid replacements over much of the protein, different mutation rates or biases in coding and non-coding sequences, negative selection acting with greater strength on non-coding than coding regions, and/or positive selection acting on the protein. Formal selection analyses provide evidence that the latter force played an important role in the basal b-like lineages coincident with the emergence of polygyny. While our data set reveals considerable paraphyly and polyphyly of S. invicta sequences with respect to those of other fire ant species, the b-like alleles of the socially polymorphic species are monophyletic. An expanded analysis of colonies containing alleles of this clade confirmed the invariant link between their presence and expression of polygyny. Finally, our discovery of several unique alleles bearing various combinations of b-like and B-like codons allows us to conclude that no single b-like residue is completely predictive of polygyne behavior and, thus, potentially causally involved in its expression. Rather, all three typical b-like residues appear to be necessary.

Identification of Schistosoma mansoni candidate antigens for diagnosis of schistosomiasis

The development of a more sensitive diagnostic test for schistosomiasis is needed to overcome the limitations of the use of stool examination in low endemic areas. Using parasite antigens in enzyme linked immunosorbent assay is a promising strategy, however a more rational selection of parasite antigens is necessary. In this study we performed in silico analysis of the Schistosoma mansoni genome, using SchistoDB database and bioinformatic tools for screening immunogenic antigens. Based on evidence of expression in all parasite life stage within the definitive host, extracellular or plasmatic membrane localization, low similarity to human and other helminthic proteins and presence of predicted B cell epitopes, six candidates were selected: a glycosylphosphatidylinositol-anchored 200 kDa protein, two putative cytochrome oxidase subunits, two expressed proteins and one hypothetical protein. The recognition in unidimensional and bidimensional Western blot of protein with similar molecular weight and isoelectric point to the selected antigens by sera from S. mansoni infected mice indicate a good correlation between these two approaches in selecting immunogenic proteins.

Selection of Reference Genes for Quantitative Real Time PCR (qPCR) Assays in Tissue from Human Ascending Aorta.

Dilatation of the ascending aorta (AAD) is a prevalent aortopathy that occurs frequently associated with bicuspid aortic valve (BAV), the most common human congenital cardiac malformation. The molecular mechanisms leading to AAD associated with BAV are still poorly understood. The search for differentially expressed genes in diseased tissue by quantitative real-time PCR (qPCR) is an invaluable tool to fill this gap. However, studies dedicated to identify reference genes necessary for normalization of mRNA expression in aortic tissue are scarce. In this report, we evaluate the qPCR expression of six candidate reference genes in tissue from the ascending aorta of 52 patients with a variety of clinical and demographic characteristics, normal and dilated aortas, and different morphologies of the aortic valve (normal aorta and normal valve n = 30; dilated aorta and normal valve n = 10; normal aorta and BAV n = 4; dilated aorta and BAV n = 8). The expression stability of the candidate reference genes was determined with three statistical algorithms, GeNorm, NormFinder and Bestkeeper. The expression analyses showed that the most stable genes for the three algorithms employed were CDKN1β, POLR2A and CASC3, independently of the structure of the aorta and the valve morphology. In conclusion, we propose the use of these three genes as reference genes for mRNA expression analysis in human ascending aorta. However, we suggest searching for specific reference genes when conducting qPCR experiments with new cohort of samples.

African signatures of recent positive selection in human FOXI1

Background: The human FOXI1 gene codes for a transcription factor involved in the physiology of the inner ear, testis, and kidney. Using three interspecies comparisons, it has been suggested that this may be a gene underhuman-specific selection. We sought to confirm this finding by using an extended set of orthologous sequences.Additionally, we explored for signals of natural selection within humans by sequencing the gene in 20 Europeans,20 East Asians and 20 Yorubas and by analysing SNP variation in a 2 Mb region centered on FOXI1 in 39worldwide human populations from the HGDP-CEPH diversity panel.Results: The genome sequences recently available from other primate and non-primate species showed that FOXI1divergence patterns are compatible with neutral evolution. Sequence-based neutrality tests were not significant inEuropeans, East Asians or Yorubas. However, the Long Range Haplotype (LRH) test, as well as the iHS and XP-Rsbstatistics revealed significantly extended tracks of homozygosity around FOXI1 in Africa, suggesting a recentepisode of positive selection acting on this gene. A functionally relevant SNP, as well as several SNPs either on theputatively selected core haplotypes or with significant iHS or XP-Rsb values, displayed allele frequencies stronglycorrelated with the absolute geographical latitude of the populations sampled.Conclusions: We present evidence for recent positive selection in the FOXI1 gene region in Africa. Climate mightbe related to this recent adaptive event in humans. Of the multiple functions of FOXI1, its role in kidney-mediatedwater-electrolyte homeostasis is the most obvious candidate for explaining a climate-related adaptation.

Adaptive model selection using empirical complexities

Given $n$ independent replicates of a jointly distributed pair $(X,Y)\in {\cal R}^d \times {\cal R}$, we wish to select from a fixed sequence of model classes ${\cal F}_1, {\cal F}_2, \ldots$ a deterministic prediction rule $f: {\cal R}^d \to {\cal R}$ whose risk is small. We investigate the possibility of empirically assessingthe {\em complexity} of each model class, that is, the actual difficulty of the estimation problem within each class. The estimated complexities are in turn used to define an adaptive model selection procedure, which is based on complexity penalized empirical risk.The available data are divided into two parts. The first is used to form an empirical cover of each model class, and the second is used to select a candidate rule from each cover based on empirical risk. The covering radii are determined empirically to optimize a tight upper bound on the estimation error. An estimate is chosen from the list of candidates in order to minimize the sum of class complexity and empirical risk. A distinguishing feature of the approach is that the complexity of each model class is assessed empirically, based on the size of its empirical cover.Finite sample performance bounds are established for the estimates, and these bounds are applied to several non-parametric estimation problems. The estimates are shown to achieve a favorable tradeoff between approximation and estimation error, and to perform as well as if the distribution-dependent complexities of the model classes were known beforehand. In addition, it is shown that the estimate can be consistent,and even possess near optimal rates of convergence, when each model class has an infinite VC or pseudo dimension.For regression estimation with squared loss we modify our estimate to achieve a faster rate of convergence.

The stationary distribution of a continuously varying strategy in a class-structured population under mutation-selection-drift balance.

Many traits and/or strategies expressed by organisms are quantitative phenotypes. Because populations are of finite size and genomes are subject to mutations, these continuously varying phenotypes are under the joint pressure of mutation, natural selection and random genetic drift. This article derives the stationary distribution for such a phenotype under a mutation-selection-drift balance in a class-structured population allowing for demographically varying class sizes and/or changing environmental conditions. The salient feature of the stationary distribution is that it can be entirely characterized in terms of the average size of the gene pool and Hamilton's inclusive fitness effect. The exploration of the phenotypic space varies exponentially with the cumulative inclusive fitness effect over state space, which determines an adaptive landscape. The peaks of the landscapes are those phenotypes that are candidate evolutionary stable strategies and can be determined by standard phenotypic selection gradient methods (e.g. evolutionary game theory, kin selection theory, adaptive dynamics). The curvature of the stationary distribution provides a measure of the stability by convergence of candidate evolutionary stable strategies, and it is evaluated explicitly for two biological scenarios: first, a coordination game, which illustrates that, for a multipeaked adaptive landscape, stochastically stable strategies can be singled out by letting the size of the gene pool grow large; second, a sex-allocation game for diploids and haplo-diploids, which suggests that the equilibrium sex ratio follows a Beta distribution with parameters depending on the features of the genetic system.

Postglacial recolonization at a snail's pace (Trochulus villosus): confronting competing refugia hypotheses using model selection.

The localization of Last Glacial Maximum (LGM) refugia is crucial information to understand a species' history and predict its reaction to future climate changes. However, many phylogeographical studies often lack sampling designs intensive enough to precisely localize these refugia. The hairy land snail Trochulus villosus has a small range centred on Switzerland, which could be intensively covered by sampling 455 individuals from 52 populations. Based on mitochondrial DNA sequences (COI and 16S), we identified two divergent lineages with distinct geographical distributions. Bayesian skyline plots suggested that both lineages expanded at the end of the LGM. To find where the origin populations were located, we applied the principles of ancestral character reconstruction and identified a candidate refugium for each mtDNA lineage: the French Jura and Central Switzerland, both ice-free during the LGM. Additional refugia, however, could not be excluded, as suggested by the microsatellite analysis of a population subset. Modelling the LGM niche of T. villosus, we showed that suitable climatic conditions were expected in the inferred refugia, but potentially also in the nunataks of the alpine ice shield. In a model selection approach, we compared several alternative recolonization scenarios by estimating the Akaike information criterion for their respective maximum-likelihood migration rates. The 'two refugia' scenario received by far the best support given the distribution of genetic diversity in T. villosus populations. Provided that fine-scale sampling designs and various analytical approaches are combined, it is possible to refine our necessary understanding of species responses to environmental changes.

RTS,S/AS02A malaria vaccine does not induce parasite CSP T cell epitope selection and reduces multiplicity of infection

Objective: The candidate malaria vaccine RTS,S/AS02A is a recombinant protein containing part of the circumsporozoite protein (CSP) sequence of Plasmodium falciparum, linked to the hepatitis B surface antigen and formulated in the proprietary adjuvant system AS02A. In a recent trial conducted in children younger than age five in southern Mozambique, the vaccinedemonstrated significant and sustained efficacy against both infection and clinical disease. In a follow-up study to the main trial, breakthrough infections identified in the trial were examined to determine whether the distribution of csp sequences was affected by the vaccine and to measure the multiplicity of infecting parasite genotypes. Design: P. falciparum DNA from isolates collected during the trial was used for genotype studies. Setting: The main trial was carried out in the Manhiça district, Maputo province, Mozambique, between April 2003 and May 2004. Participants: Children from the two cohorts of the main trial provided parasite isolates as follows: children from Cohort 1 who were admitted to hospital with clinical malaria; children from Cohort 1 who were parasite-positive in a cross-sectional survey at study month 8.5; children from Cohort 2 identified as parasite-positive during follow-up by active detection of infection. Outcome: Divergence of DNA sequence encoding the CSP T cell-epitope region sequence from that of the vaccine sequence was measured in 521 isolates. The number of distinct P. falciparum genotypes was also determined. Results: We found no evidence that parasite genotypes from children in the RTS,S/AS02A arm were more divergent than those receiving control vaccines. For Cohort 1 (survey at studymonth 8.5) and Cohort 2, infections in the vaccine group contained significantly fewer genotypes than those in the control group, (p 1/4 0.035, p 1/4 0.006), respectively, for the two cohorts. This was not the case for children in Cohort 1 who were admitted to hospital (p 1/4 0.478). Conclusions: RTS,S/AS02A did not select for genotypes encoding divergent T cell epitopes in the C-terminal region of CSP in this trial. In both cohorts, there was a modest reduction in the mean number of parasite genotypes harboured by vaccinated children compared with controls, but only among those with asymptomatic infections.