SEARCH · PubMed Health

Results for “Purifying selection”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12Linked to original sources

PubMed2005

Evolution of the relaxin-like peptide family.

BACKGROUND: The relaxin-like peptide family belongs in the insulin superfamily and consists of 7 peptides of high structural but low sequence similarity; relaxin-1, 2 and 3, and the insulin-like (INSL) peptides, INSL3, INSL4, INSL5 and INSL6. The functions of relaxin-3, INSL4, INSL5, INSL6 remain uncharacterised. The evolution of this family has been contentious; high sequence variability is seen between closely related species, while distantly related species show high similarity; an invertebrate relaxin sequence has been reported, while a relaxin gene has not been found in the avian and ruminant lineages. RESULTS: Sequence similarity searches of genomic and EST data identified homologs of relaxin-like peptides in mammals, and non-mammalian vertebrates such as fish. Phylogenetic analysis was used to resolve the evolution of the family. Searches were unable to identify an invertebrate relaxin-like peptide. The published relaxin cDNA sequence in the tunicate, Ciona intestinalis was not present in the completed C. intestinalis genome. The newly discovered relaxin-3 is likely to be the ancestral relaxin. Multiple relaxin-3-like sequences are present in fugu fish (Takifugu rubripes) and zebrafish (Danio rerio), but these appear to be specific to the fish lineage. Possible relaxin-1 and INSL5 homologs were also identified in fish and frog species, placing their emergence prior to mammalia, earlier than previously believed. Furthermore, estimates of synonymous and nonsynonymous substitution rates (dN/dS) suggest that the emergence of relaxin-1, INSL4 and INSL6 during mammalia was driven by positive Darwinian selection, hence these peptides are likely to have novel and in the case of relaxin-1, which is still under positive selection in humans and the great apes, possibly still evolving functions. In contrast, relaxin-3 is constrained by strong purifying selection, demonstrating it must have a highly conserved function, supporting its hypothesized important neuropeptide role. CONCLUSIONS: We present a phylogeny describing the evolutionary history of the relaxin-like peptide family and show that positive selection has driven the evolution of the most recent members of the family.

Amino Acid Sequence↗

PubMed2006

Evidence for overdominant selection maintaining X-linked fitness variation in Drosophila melanogaster.

The role of balancing selection in maintaining genetic variation for fitness is largely unresolved. This reflects the inherent difficulty in distinguishing between models of recurrent mutation versus selection, which produce similar patterns of inbreeding depression, as well as the limitations of testing such hypotheses when fitness variation is averaged across the genome. Signatures of X-linked overdominant selection are less likely to be obscured by mutational variation because X-linked mutations are rapidly eliminated by purifying selection in males. Although models maintaining genetic variation for fitness are not necessarily mutually exclusive, a series of predictions for identifying X-linked overdominant selection can be used to separate its contribution from other underlying processes. We consider the role of overdominant selection in maintaining fitness variation in a sample of 12 X chromosomes from a population of Drosophila melanogaster. Substantial variation was observed for male reproductive success and female fecundity, with heterozygous-X genotypes exhibiting the greatest degree of variance, a finding that agrees well with predictions of the overdominance model. The importance of X-linked overdominant selection is discussed along with models of recurrent mutation and sexually antagonistic selection.

Alleles↗

PubMed2003

A novel sensitive high-performance liquid chromatography/electrochemical procedure for measuring formaldehyde produced from oxidative deamination of methylamine and in biological samples.

Formaldehyde is a well-known environmental toxic hazard. It is also a product of oxidative deamination of methylamine catalyzed by semicarbazide-sensitive amine oxidase (SSAO). Increased SSAO-mediated deamination has been implicated in some pathophysiological conditions, such as diabetic complications. The measurement of formaldehyde in the enzymatic reactions and in vivo production using conventional methods was not straightforward due to limitations of selectivity and sensitivity. A novel high-performance liquid chromatography (HPLC)/electrochemical procedure for the measurement of formaldehyde has been developed. The measurement is based on the formation of adducts between formaldehyde and dopamine. These adducts can be selectively purified and concentrated using a batch method of alumina absorption, separated by HPLC, and electrochemically quantified. The method is highly selective and substantially more sensitive, i.e., detection of picomole levels of formaldehyde, than the conventional methods. The procedure not only facilitates the assessment of SSAO activity in vitro but also is useful for assessing formaldehyde in tissues and biological fluids.

Animals↗

PubMed2005

Interchromosomal segmental duplications explain the unusual structure of PRSS3, the gene for an inhibitor-resistant trypsinogen.

Homo sapiens possess several trypsinogen or trypsinogen-like genes of which three (PRSS1, PRSS2, and PRSS3) produce functional trypsins in the digestive tract. PRSS1 and PRSS2 are located on chromosome 7q35, while PRSS3 is found on chromosome 9p13. Here, we report a variation of the theme of new gene creation by duplication: the PRSS3 gene was formed by segmental duplications originating from chromosomes 7q35 and 11q24. As a result, PRSS3 transcripts display two variants of exon 1. The PRSS3 transcript whose gene organization most resembles PRSS1 and PRSS2 encodes a functional protein originally named mesotrypsinogen. The other variant is a fusion transcript, called trypsinogen IV. We show that the first exon of trypsinogen IV is derived from the noncoding first exon of LOC120224, a chromosome 11 gene. LOC120224 codes for a widely conserved transmembrane protein of unknown function. Comparative analyses suggest that these interchromosomal duplications occurred after the divergence of Old World monkeys and hominids. PRSS3 transcripts consist of a mixed population of mRNAs, some expressed in the pancreas and encoding an apparently functional trypsinogen and others of unknown function expressed in brain and a variety of other tissues. Analysis of the selection pressures acting on the trypsinogen gene family shows that, while the apparently functional genes are under mild to strong purifying selection overall, a few residues appear under positive selection. These residues could be involved in interactions with inhibitors.

Chromosomes, Human↗

PubMed2003

DNA polymorphism and selection at the melanocortin-1 receptor gene in normally pigmented southern African individuals.

Skin pigmentation is a polygenic multifactorial trait determined by the cumulative effects of multiple genetic variants and environmental factors. Melanocortin-1 receptor (MC1R) is one of the genes involved in pigmentation, and has been implicated in the red hair and pale skin phenotype in human Caucasoid individuals. The present study was undertaken to identify variation at the MC1R locus in normally pigmented individuals in two African populations, sub-Saharan Negroids (22 unrelated individuals) and the San (17 unrelated individuals). The study showed considerable MC1R gene sequence variation with the detection of eight synonymous and three nonsynonymous mutations. This is the first report of nonsynonymous mutations in African individuals in the MC1R gene: L99I was found in a single San individual, S47I was detected in a single Negroid individual, and F196L was detected in five Negroid individuals (5/44; 0.11). The functional significance of these mutations is not known. Three of the eight synonymous mutations found, L106L (CTG --> CTA), F300F (TTC --> TTT), and T314T (ACA --> ACG) (also known as A942G), have been reported previously. T314T was the only variant that showed a significant difference between the Negroid and San populations (0.477 and 0.059, respectively; P = 1.6 x 10(-5)). Its low frequency in the San may be the result of random genetic drift in a population of small size, or selection. Several tests of neutrality of the MC1R coding region in these and other African populations were significant, suggesting that purifying selection (functional constraint) had occurred at this gene locus in Africans. This demonstrates that although some nonsynonymous MC1R mutations are tolerated in individuals with dark skin, this gene has likely played a significant role in the maintenance of dark pigmentation in Africans and normal pigment variation in non-African populations.

Africa↗

PubMed1991

The structure of human mitochondrial DNA variation.

Restriction analysis of mitochondrial DNA (mtDNA) of 3065 humans from 62 geographic samples identified 149 haplotypes and 81 polymorphic sites. These data were used to test several aspects of the evolutionary past of the human species. A dendrogram depicting the genetic relatedness of all haplotypes shows that the native African populations have the greatest diversity and, consistent with evidence from a variety of sources, suggests an African origin for our species. The data also indicate that two individuals drawn at random from the entire sample will differ at approximately 0.4% of their mtDNA nucleotide sites, which is somewhat higher than previous estimates. Human mtDNA also exhibits more interpopulation heterogeneity (GST = 0.351 +/- 0.025) than does nuclear DNA (GST = 0.12). Moreover, the virtual absence of intermediate levels of linkage disequilibrium between pairs of sites is consistent with the absence of genetic recombination and places constraints on the rate of mutation. Tests of the selective neutrality of mtDNA variation, including the Ewens-Watterson and Tajima tests, indicate a departure in the direction consistent with purifying selection, but this departure is more likely due to the rapid growth of the human population and the geographic heterogeneity of the variation. The lack of a good fit to neutrality poses problems for the estimation of times of coalescence from human mtDNA data.

Alleles↗

PubMed2005

Natural selection on mitochondrial DNA in Parus and its relevance for phylogeographic studies.

I examined mitochondrial ND2 sequences in species of tits (genus Parus) to determine whether substitutions were neutral or under selection. Haplotype trees in these species are shallow and geographically unstructured, which could be a signature of recent demographic expansion or purifying selection. McDonald-Kreitman tests revealed a pattern of excess replacement polymorphisms for closely related taxa. Replacement substitutions tended to be less common in the basal parts of haplotype trees, suggesting that they are selected against over evolutionary time. Thus, the pattern of selection is consistent with a model of mildly deleterious haplotypes. To investigate the significance of this selection, the ND2 gene was subdivided into surface and transmembrane portions, and synonymous-non-synonymous sites, and F(ST)-values (the amount of variance distributed among populations) computed for each. Among sampling localities in P. major, both partitions gave similar estimates of F(ST), which is consistent with demographic expansion as the cause of the shallow trees. In P. montanus, two localities contained high percentages of individuals that had a single replacement substitution in the transmembrane portion of the gene, which inflated F(ST)-values relative to the surface sites. Thus, the interpretation of the genetic differentiation of these two sites could be either geographical isolation or selection. Given that the sole substitution causing the high F(ST)-values was a replacement one, selection is implicated. However, this substitution occurred in only 2 of 12 localities, revealing that overall the phylogeographic pattern was not biased by selection, once the cause of the two high pairwise F(ST)-comparisons was recognized. Investigators should consider the source of significant genetic heterogeneity within species, to determine if the variation is likely caused by demographic isolation or selection. However, it seems unlikely that most phylogeographic analyses are invalidated by the effects of natural selection.

Animals↗

PubMed2001

A strategy for antitumor vascular therapy by targeting the vascular endothelial growth factor: receptor complex.

Vascular endothelial growth factor (VEGF) is produced by cancer cells in response to hypoxia and is the primary stimulant of vascularization in solid tumors. Endothelial cells lining the blood vessels of these tumors have a high concentration of receptor-bound VEGF on their surface, providing a target for antibody- directed cancer therapy. To obtain a cloned antibody to this target when bound to its receptor on tumor endothelium, we used phage display technology to create a single-chain Fv (sFv) antibody library from mice immunized with the 165-amino acid isoform of human VEGF-A. We selected, purified, and characterized LL4, an anti-VEGF sFv that was shown to react with receptor-bound VEGF. LL4 bound selectively to blood vessel endothelium, as shown by immunohistochemistry on tissue sections of human tumors. Furthermore, using autoradiography and grain counting of histological sections, systemically administered LL4 was shown to localize selectively to the endothelial lining of tumor blood vessels in human colorectal carcinoma xenografts in vivo. This study demonstrates the feasibility of targeting tumor vasculature using recombinant antibodies to the VEGF:receptor complex.

Animals↗

PubMed2004

The effect of selection on genealogies.

The coalescent process can describe the effects of selection at linked loci only if selection is so strong that genotype frequencies evolve deterministically. Here, we develop methods proposed by Kaplan, Darden, and Hudson to find the effects of weak selection. We show that the overall effect is given by an extension to Price's equation: the change in properties such as moments of coalescence times is equal to the covariance between those properties and the fitness of the sample of genes. The distribution of coalescence times differs substantially between allelic classes, even in the absence of selection. However, the average coalescence time between randomly chosen genes is insensitive to the current allele frequency and is affected significantly by purifying selection only if deleterious mutations are common and selection is strong (i.e., the product of population size and selection coefficient, Ns>3). Balancing selection increases mean coalescence times, but the effect becomes large only when mutation rates between allelic classes are low and when selection is extremely strong. Our analysis supports previous simulations that show that selection has surprisingly little effect on genealogies. Moreover, small fluctuations in allele frequency due to random drift can greatly reduce any such effects. This will make it difficult to detect the action of selection from neutral variation alone.

Biological Evolution↗

PubMed2005

Evolutionary dynamics of hepatitis C virus envelope genes during chronic infection.

Hepatitis C virus (HCV) envelope glycoproteins E1 and E2 are important targets for the host immune response. The genes encoding these proteins exhibit a high degree of variability that gives rise to differing phenotypic traits, including alterations in receptor-binding affinity and immune recognition and escape. In order to elucidate patterns of adaptive evolution during chronic infection, a panel of full-length E1E2 clones was generated from sequential serum samples obtained from four chronically infected individuals. By using likelihood-based methods for phylogenetic inference, the evolutionary dynamics of circulating HCV quasispecies populations were assessed and a site-by-site analysis of the d(N)/d(S) ratio was performed, to identify specific codons undergoing diversifying positive selection. HCV phylogenies, coupled with the number and distribution of selected sites, differed markedly between patients, highlighting that HCV evolution during chronic infection is a patient-specific phenomenon. This analysis shows that purifying selection is the major force acting on HCV populations in chronic infection. Whilst no significant evidence for positive selection was observed in E1, a number of sites under positive selection were identified within the ectodomain of the E2 protein. All of these sites were located in regions hypothesized to be exposed to the selective environment of the host, including a number of functionally defined domains that have been reported to be involved in immune evasion and receptor binding. Dated-tip methods for estimation of underlying HCV mutation rates were also applied to the data, enabling prediction of the most recent common ancestor for each patient's quasispecies.

Amino Acid Sequence↗

PubMed2003

Protein evolution and codon usage bias on the neo-sex chromosomes of Drosophila miranda.

The neo-sex chromosomes of Drosophila miranda constitute an ideal system to study the effects of recombination on patterns of genome evolution. Due to a fusion of an autosome with the Y chromosome, one homolog is transmitted clonally. Here, I compare patterns of molecular evolution of 18 protein-coding genes located on the recombining neo-X and their homologs on the nonrecombining neo-Y chromosome. The rate of protein evolution has significantly increased on the neo-Y lineage since its formation. Amino acid substitutions are accumulating uniformly among neo-Y-linked genes, as expected if all loci on the neo-Y chromosome suffer from a reduced effectiveness of natural selection. In contrast, there is significant heterogeneity in the rate of protein evolution among neo-X-linked genes, with most loci being under strong purifying selection and two genes showing evidence for adaptive evolution. This observation agrees with theory predicting that linkage limits adaptive protein evolution. Both the neo-X and the neo-Y chromosome show an excess of unpreferred codon substitutions over preferred ones and no difference in this pattern was observed between the chromosomes. This suggests that there has been little or no selection maintaining codon bias in the D. miranda lineage. A change in mutational bias toward AT substitutions also contributes to the decline in codon bias. The contrast in patterns of molecular evolution between amino acid mutations and synonymous mutations on the neo-sex-linked genes can be understood in terms of chromosome-specific differences in effective population size and the distribution of selective effects of mutations.

Animals↗

PubMed2003

Mutation accumulation in growing asexual lineages.

The stochastic loss of entire classes of individuals bearing the fewest number of mutations-a process known as Muller's ratchet-is studied in asexual populations growing unconstrained from a single founder. In the neutral regime, where mutations have zero effect on fitness, we derive a recursion equation for the probability distribution of the minimum number of mutations carried by individuals in the least-loaded class, and obtain an explicit condition for the halting of the ratchet. Next, we consider the case of deleterious mutations, and show that weak selection can actually accelerate the ratchet beyond that achieved for the neutral regime. This effect is transitory, however, as our results suggest that even weak purifying selection will eventually lead to the complete cessation of the ratchet. These results may have important implications for problems in biology and the medical sciences.

Models, Genetic↗

PubMed2002

On the adsorption behaviour of saliva and purified salivary proteins at solid/liquid interfaces.

Salivary proteinaceous substances are known to play important roles in the formation of the salivary pellicle. The aim of this study was to investigate some aspects of the interfacial behaviour of selected purified salivary proteins, as well as human saliva secretions, using time-resolved in situ ellipsometry. Hydrophobic methylated silica and hydrophilic pure silica were used as test substrates. Experiments were performed in vitro, preferentially in the low concentration range, with samples of fresh human whole resting saliva, parotid resting saliva and submandibular/sublingual resting saliva. The protein fractions investigated were human MUC5B, PRP-1, PRP-3 and statherin, as well as bovine submaxillary mucin (BSM). The results indicate that the amount of material adsorbed was strongly related to the protein concentration in the range investigated for both pure proteins and secretions. Generally, a larger amount of material was adsorbed onto hydrophobic surfaces than onto hydrophilic ones. However, pure PRP-1 adsorbed in similar amounts to both hydrophilic and hydrophobic surfaces in the concentration range investigated and BSM adsorbed in larger amounts at high concentrations on hydrophilic surfaces. Comparison of the observed adsorption rates for salivary secretions and calculated diffusion rates for individual proteins suggested initial adsorption of low-molecular-weight proteins/peptides. On hydrophilic surfaces the data indicate adsorption of proteins with diffusion rates corresponding to those of statherin, PRP-3 and PRP-1. MUC5B adsorbs at a later stage from both HWS and the individual secretions, which can be explained by a "Vroman effect"-like phenomenon. On hydrophobic surfaces, adsorption rates were found to be faster than those calculated for any of the proteins, and thus smaller proteins/peptides appear to be involved. The similar adsorption behaviour of PRP-1 and parotid saliva (HPS) on hydrophilic surfaces may suggest that long aPRPs account for a substantial portion of the film-forming capacity of HPS. Effects of added electrolyte could be explained by general screening effects and specific Ca2+ binding to serine phosphates in aqueous solutions, but were complex in phosphate buffer. Inter-individual differences in amounts adsorbed from HWS, HPS and HSMSLS, respectively, were not found to be statistically significant.

Adsorption↗

PubMed2004

Evidence of selection on the domesticated ERVWE1 env retroviral element involved in placentation.

The human endogenous retrovirus HERV-W multicopy family includes a unique proviral locus, termed ERVWE1, which contains gag and pol pseudogenes and has retained a full-length envelope open reading frame (ORF). This Env protein (syncytin) is a highly fusogenic membrane glycoprotein and has been proposed to be involved in hominoid placental physiology. To track the hallmarks of natural selection acting on the ERVWE1 env gene, the pattern of substitutions and indels was analyzed within all human HERV-W elements and along the ERVWE1 orthologous loci in chimpanzee, gorilla, orangutan, and gibbon. The comparison of ERVWE1 and paralogous HERV-W copies revealed an ERVWE1-specific signature consisting of a four amino acid deletion in the intracytoplasmic tail of the glycoprotein. We show that this deletion is crucial for the envelope fusogenic activity. The comparison of the human ERVWE1 locus with its orthologs demonstrates the existence of a selective pressure to maintain the env reading frame open. Notably, the 3' part of the env gene, encoding regions required for the fusion process, is under purifying selection. The identification of selective constraints on env ERVWE1 confirms that this retroviral locus has been recruited in the hominoid lineage to become a bona fide gene.

Amino Acid Sequence↗

PubMed2004

Different regulatory mechanisms underlie similar transposable element profiles in pufferfish and fruitflies.

Comparative analysis of recently sequenced eukaryotic genomes has uncovered extensive variation in transposable element (TE) abundance, diversity, and distribution. The TE profile in the sequenced pufferfish genomes is more similar to that of Drosophila melanogaster than to human or mouse, in that pufferfish TEs exhibit low overall abundance, high family diversity, and localization in the heterochromatin. It has been suggested that selection against the deleterious effects of ectopic recombination between TEs has structured the TE profile in Drosophila and pufferfish but not in humans. We test this hypothesis by measuring the sample frequency of 48 euchromatic TE insertions in the genome of the green spotted pufferfish (Tetraodon nigroviridis). We estimate the strength of selection acting on recent insertions by analyzing the site frequency spectrum using a maximum-likelihood approach. We show that in contrast to Drosophila, euchromatic TE insertions in Tetraodon are selectively neutral and that the low copy number and compartmentalized distribution of TEs in the Tetraodon genome must be caused by regulation by means other than purifying selection acting on recent insertions. Inference of regulatory processes governing TE profiles should take into account factors such as effective population size, incidence of inbreeding/outcrossing, and other species-specific traits.

Animals↗

PubMed2006

Inferring the distribution of mutational effects on fitness in Drosophila.

The properties of the distribution of deleterious mutational effects on fitness (DDME) are of fundamental importance for evolutionary genetics. Since it is extremely difficult to determine the nature of this distribution, several methods using various assumptions about the DDME have been developed, for the purpose of parameter estimation. We apply a newly developed method to DNA sequence polymorphism data from two Drosophila species and compare estimates of the parameters of the distribution of the heterozygous fitness effects of amino acid mutations for several different distribution functions. The results exclude normal and gamma distributions, since these predict too few effectively lethal mutations and power-law distributions as a result of predicting too many lethals. Only the lognormal distribution appears to fit both the diversity data and the frequency of lethals. This DDME arises naturally in complex systems when independent factors contribute multiplicatively to an increase in fitness-reducing damage. Several important parameters, such as the fraction of effectively neutral non-synonymous mutations and the harmonic mean of non-neutral selection coefficients, are robust to the form of the DDME. Our results suggest that the majority of non-synonymous mutations in Drosophila are under effective purifying selection.

Animals↗

PubMed2026

Comparative analysis of olfactory receptor repertoires reveals evolutionary dynamics and high-altitude adaptation in Schizopygopsis younghusbandi based on the chromosome-level genomes.

The olfactory receptor (OR) gene represent a significant multigene family in vertebrates, forming the core molecular basis of olfactory perception and playing a crucial role in the environmental adaptation of species. High-altitude ecosystems represent extreme habitats characterized by specific abiotic stresses, including low oxygen levels, low temperatures, and intense ultraviolet radiation. These environments also exhibit low aquatic biodiversity and a limited variety of odor molecules, factors that have influenced the adaptive evolution of the sensory systems in endemic species. However, the genetic mechanisms underlying olfactory adaptation in high-altitude freshwater fish remained inadequately understood. In this study, we performed comparative genomics analyses to reveal the evolutionary processes underlying the adaptive and functional evolution of OR genes in S. younghusbandi, a cyprinid fish endemic to the Qinghai-Xizang Plateau. The results indicated that, compared to their low-altitude relatives, S. younghusbandi possessed a significantly smaller number of OR genes, with only 98 genes, which revealed the contraction of the gene family. Phylogenetic analysis revealed that the OR genes of cyprinid fish could be categorized into two major lineages: type I and type II. The η and δ families, which perceive water-soluble odors, in S. younghusbandi underwent significant and specific expansion, while the ε family was completely absent. This pattern reflected adaptive changes in olfactory recognition to accommodate the simplified odor spectrum of high-altitude water bodies. Chromosomal localization analysis demonstrated that OR genes were clustered, and collinearity analysis confirmed the presence of conserved genomic fragments among species. Selection pressure analysis revealed that the Ka/Ks values of all homologous gene pairs were less than 1, indicating that the OR genes of S. younghusbandi underwent strong purifying selection as a group to preserve core olfactory function. A few genes exhibited relaxed selection characteristics, which may have facilitated the fine-tuning of adaptability to high-altitude environments. In conclusion, this study elucidated the evolutionary dynamics and adaptive characteristics of the OR gene in S. younghusbandi, offering a new perspective on the molecular mechanisms underlying olfactory adaptation at high altitudes and enriching the research on sensory evolution in vertebrates.

Schizopygopsis younghusbandi↗

PubMed2006

Sca1, a previously undescribed paralog from autotransporter protein-encoding genes in Rickettsia species.

BACKGROUND: Among the 17 genes encoding autotransporter proteins of the "surface cell antigen" (sca) family in the currently sequenced Rickettsia genomes, ompA, sca5 (ompB) and sca4 (gene D), have been extensively used for identification and phylogenetic purposes for Rickettsia species. However, none of these genes is present in all 20 currently validated Rickettsia species. Of the remaining 14 sca genes, sca1 is the only gene to be present in all nine sequenced Rickettsia genomes. To estimate whether the sca1 gene is present in all Rickettsia species and its usefulness as an identification and phylogenetic tool, we searched for sca1genes in the four published Rickettsia genomes and amplified and sequenced this gene in the remaining 16 validated Rickettsia species. RESULTS: Sca1 is the only one of the 17 rickettsial sca genes present in all 20 Rickettsia species. R. prowazekii and R. canadensis exhibit a split sca1 gene whereas the remaining species have a complete gene. Within the sca1 gene, we identified a 488-bp variable sequence fragment that can be amplified using a pair of conserved primers. Sequences of this fragment are specific for each Rickettsia species. The phylogenetic organization of Rickettsia species inferred from the comparison of sca1 sequences strengthens the classification based on the housekeeping gene gltA and is similar to those obtained from the analyses of ompA, sca5 and sca4, thus suggesting similar evolutionary constraints. We also observed that Sca1 protein sequences have evolved under a dual selection pressure: with the exception of typhus group rickettsiae, the amino-terminal part of the protein that encompasses the predicted passenger domain, has evolved under positive selection in rickettsiae. This suggests that the Sca1 protein interacts with the host. In contrast, the C-terminal portion containing the autotransporter domain has evolved under purifying selection. In addition, sca1 is transcribed in R. conorii, and might therefore be functional in this species. CONCLUSION: The sca1 gene, encoding an autotransporter protein that evolves under dual evolution pressure, is the only sca-family gene to be conserved by all Rickettsia species. As such, it is a valuable identification target for these bacteria, especially because rickettsial isolates can be identified by amplification and sequencing of a discriminatory gene fragment using a single primer pair. It may also be used as a phylogenetic tool. However, its current functional status remains to be determined although it was found expressed in R. conorii.

Antigens, Surface↗