Worldwide genetic variation at the 3 ' untranslated region of the HLA-G gene: balancing selection influencing genetic diversity

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Papel da hidrofobicidade na evolução de proteínas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Genotyping of Human parvovirus B19 among Brazilian patients with hemoglobinopathies

Human parvovirus B19 (B19V) infection can be a life-threatening condition among patients with hereditary (chronic) hemolytic anemias. Our objective was to characterize the infection molecularly among patients with sickle cell disease and thalassemia. Forty-seven patients (37 with sickle cell disease, and 10 with beta-thalassemia major) as well as 47 healthy blood donors were examined for B19V infection by anti-B19V IgG enzyme immunoassay, quantitative PCR, which detects all B19V genotypes, and DNA sequencing. B19V viremia was documented in nine patients (19.1%) as two displayed acute infection and the rest had a low titre viremia (mean 3.4 x 10(4) copies/mL). All donors were negative for B19V DNA. Anti-B19V IgG was detected in 55.3% of the patients and 57.4% among the donors. Based on partial NS1 fragments, all patient isolates were classified as genotype 1 and subgenotype 1A. The evolutionary events of the examined partial NS1 gene sequence were associated with a lack of positive selection. The quantification of all B19V genotypes by a single hydrolytic probe is a technically useful method, but it is difficult to establish relationships between B19V sequence characteristics and infection outcome.

Epidemiological and functional characterization of Streptococcus pneumoniae pili

Streptococcus pneumoniae is an important life threatening human pathogen causing agent of invasive diseases such as otitis media, pneumonia, sepsis and meningitis, but is also a common inhabitant of the respiratory tract of children and healthy adults. Likewise most streptococci, S. pneumoniae decorates its surface with adhesive pili, composed of covalently linked subunits and involved in the attachment to epithelial cells and virulence. The pneumococcal pili are encoded by two genomic regions, pilus islet 1 (PI-1), and pilus islet-2 (PI-2), which are present in about 30% and 16% of the pneumococcal strains, respectively. PI-1 exists in three clonally related variants, whereas PI-2 is highly conserved. The presence of the islets does not correlate with the serotype of the strains, but with the genotype (as determined by Multi Locus Sequence Typing). The prevalence of PI-1 and PI-2 positive strains is similar in isolates from invasive disease and carriage. To better dissect a possible association between PIs presence and disease we evaluated the distribution of the two PIs in a panel of 113 acute otitis media (AOM) clinical isolates from Israel. PI-1 was present in 30.1% (N=34) of the isolates tested, and PI-2 in 7% (N=8). We found that 50% of the PI-1 positive isolates belonged to the international clones Spain9V-3 (ST156) and Taiwan19F-14 (ST236), and that PI-2 was not present in the absence of Pl-1. In conclusion, there was no correlation between PIs presence and AOM, and, in general, the observed differences in PIs prevalence are strictly dependent upon regional differences in the distribution of the clones. Finally, in the AOM collection the prevalence of PI-1 was higher among antibiotic resistant isolates, confirming previous indications obtained by the in silico analysis of the MLST database collection. Since the pilus-1 subunits were shown to confer protection in mouse models of infection both in active and passive immunization studies, and were regarded as potential candidates for a new generation of protein-based vaccines, the functional characterization was mainly focused on S. pneumoniae pilus -1 components. The pneumococcal pilus-1 is composed of three subunits, RrgA, RrgB and RrgC, each stabilized by intra-molecular isopeptide bonds and covalently polymerized by means of inter-molecular isopeptide bonds to form an extended fibre. The pilus shaft is a multimeric structure mainly composed by the RrgB backbone subunit. The minor ancillary proteins are located at the tip and at the base of the pilus, where they have been proposed to act as the major adhesin (RrgA) and as the pilus anchor (RrgC), respectively. RrgA is protective in in vivo mouse models, and exists in two variants (clades I and II). Mapping of the sequence variability onto the RrgA structure predicted from X-ray data showed that the diversity was restricted to the “head” of the protein, which contains the putative binding domains, whereas the elongated “stalk” was mostly conserved. To investigate whether this variability could influence the adhesive capacity of RrgA and to map the regions important for binding, two full-length protein variants and three recombinant RrgA portions were tested for adhesion to lung epithelial cells and to purified extracellular matrix (ECM) components. The two RrgA variants displayed similar binding abilities, whereas none of the recombinant fragments adhered at levels comparable to those of the full-length protein, suggesting that proper folding and structural arrangement are crucial to retain protein functionality. Furthermore, the two RrgA variants were shown to be cross-reactive in vitro and cross-protective in vivo in a murine model of passive immunization. Taken together, these data indicate that the region implicated in adhesion and the functional epitopes responsible for the protective ability of RrgA may be conserved and that the considerable level of variation found within the “head” domain of RrgA may have been generated by immunologic pressure without impairing the functional integrity of the pilus.

Analysis of the immunological and functional features of the Neisserial Heparin Binding Antigen (NHBA)

Neisserial Heparin Binding Antigen (NHBA) is a surface-exposed lipoprotein ubiquitously expressed by genetically diverse Neisseria meningitidis strains and is an antigen of the multicomponent protein-based 4CMenB vaccine, able to induce bactericidal antibodies in humans and to bind heparin-like molecules. The aim of this study is to characterize the immunological and functional properties of NHBA. To evaluate immunogenicity and the contribution of aminoacid sequence variability to vaccine coverage, we constructed recombinant isogenic strains that are susceptible to bactericidal killing only by anti-NHBA antibodies and engineered them to express equal levels of selected NHBA peptides. In these recombinant strains, we observed different titres associated with the different peptide variants. These recombinant strains were then further engineered to express NHBA chimeric proteins to investigate the regions important for immunogenicity. In natural strains, anti-NHBA antibodies were found to be cross-protective against strains expressing different peptides. To investigate the functional properties of this antigen, the recombinant purified NHBA protein was tested in in vitro binding studies and was found to be able to bind epithelial cells. The binding was abolished when cells were treated specifically with heparinase III, suggesting that the interaction with the cells is mediated by heparan sulfate proteoglycans (HSPG). Mutation of the Arg-rich tract of NHBA abrogated the binding, confirming the importance of this region in mediating the binding to heparin-like molecules. In a panel of N. meningitidis strains, the deletion of nhba resulted in a reduction of adhesion with respect to each isogenic wild type strain. Furthermore, the adhesion of the wild-type strain was prevented by using anti-NHBA polyclonal sera, demonstrating the specificity of the interaction. These results suggest that NHBA could be a novel meningococcal adhesin contributing to host-cell interaction. Moreover, we analysed NHBA NalP-mediated cleavage in different NHBA peptides and showed that not all NHBA peptides are cleaved.

Evidence for balancing selection operating at the het-c heterokaryon incompatibility locus in a group of filamentous fungi

In filamentous fungi, het loci (for heterokaryon incompatibility) are believed to regulate self/nonself-recognition during vegetative growth. As filamentous fungi grow, hyphal fusion occurs within an individual colony to form a network. Hyphal fusion can occur also between different individuals to form a heterokaryon, in which genetically distinct nuclei occupy a common cytoplasm. However, heterokaryotic cells are viable only if the individuals involved have identical alleles at all het loci. One het locus, het-c, has been characterized at the molecular level in Neurospora crassa and encodes a glycine-rich protein. In an effort to understand the role of this locus in filamentous fungi, we chose to study its evolution by analyzing het-c sequence variability in species within Neurospora and related genera. We determined that the het-c locus was polymorphic in a field population of N. crassa with close to equal frequency of each of the three allelic types. Different species and even genera within the Sordariaceae shared het-c polymorphisms, indicating that these polymorphisms originated in an ancestral species. Finally, an analysis of the het-c specificity region shows a high occurrence of nonsynonymous substitution. The persistence of allelic lineages, the nearly equal allelic distribution within populations, and the high frequency of nonsynonymous substitutions in the het-c specificity region suggest that balancing selection has operated to maintain allelic diversity at het-c. Het-c shares this particular evolutionary characteristic of departing from neutrality with other self/nonself-recognition systems such as major histocompatibility complex loci in mammals and the S (self-incompatibility) locus in angiosperms.

Rapid evolution of the DNA-binding site in LAGLIDADG homing endonucleases

Sequence analysis of chloroplast and mitochondrial large subunit rRNA genes from over 75 green algae disclosed 28 new group I intron-encoded proteins carrying a single LAGLIDADG motif. These putative homing endonucleases form four subfamilies of homologous enzymes, with the members of each subfamily being encoded by introns sharing the same insertion site. We showed that four divergent endonucleases from the I-CreI subfamily cleave the same DNA substrates. Mapping of the 66 amino acids that are conserved among the members of this subfamily on the 3-dimensional structure of I-CreI bound to its recognition sequence revealed that these residues participate in protein folding, homodimerization, DNA recognition and catalysis. Surprisingly, only seven of the 21 I-CreI amino acids interacting with DNA are conserved, suggesting that I-CreI and its homologs use different subsets of residues to recognize the same DNA sequence. Our sequence comparison of all 45 single-LAGLIDADG proteins identified so far suggests that these proteins share related structures and that there is a weak pressure in each subfamily to maintain identical protein–DNA contacts. The high sequence variability we observed in the DNA-binding site of homologous LAGLIDADG endonucleases provides insight into how these proteins evolve new DNA specificity.

Human brain contains high levels of heteroplasmy in the noncoding regions of mitochondrial DNA.

We have analyzed the level of intraindividual sequence variability (heteroplasmy) of mtDNA in human brain by denaturing gradient gel electrophoresis and sequencing. Single base substitutions, as well as insertions or deletions of single bases, were numerous in the noncoding control region (D-loop), and 35-45% of the molecules from a single tissue showed sequence differences. By contrast, heteroplasmy in coding regions was not detected. The lower level of heteroplasmy in the coding regions is indicative of selection against deleterious mutations. Similar levels of heteroplasmy were found in two brain regions from the same individual, while no heteroplasmy was detected in blood. Thus, heteroplasmy seems to be more frequent in nonmitotic tissues. We observed a 7.7-fold increase in the frequency of deletions/insertions and a 2.2-fold increase in the overall frequency of heteroplasmic mutations in two individuals aged 96 and 99, relative to an individual aged 28. Our results show that intraindividual sequence variability occurs at a high frequency in the noncoding regions of normal human brain and indicate that small insertions and deletions might accumulate with age at a lower rate than large rearrangements.

Myc and Max: molecular evolution of a family of proto-oncogene products and their dimerization partner.

The myc gene family encodes a group of transcription factors that regulate cell proliferation and differentiation. These genes are widely studied because of their importance as proto-oncogenes. Phylogenetic analyses are described here for 45 Myc protein sequences representing c-, N-, L-, S-, and B-myc genes. A gene duplication early in vertebrate evolution produced the c-myc lineage and another lineage that later gave rise to the N- and L-myc lineages by another gene duplication. Evolutionary divergence in the myc gene family corresponds closely to the known branching order of the major vertebrate groups. The patterns of sequence evolution are described for five separate highly conserved regions, and these analyses show that differential rates of sequence divergence (= mosaic evolution) have occurred among conserved motifs. Further, the closely related dimerization partner protein Max exhibits significantly less sequence variability than Myc. It is suggested that the reduced variability in max stems from natural selection acting to preserve dimerization capability with products of myc and related genes.

Genetic variation of the scleractinian coral Stylophora pistillata, from western Pacific reefs

Most scleractinian coral species are widely distributed across the tropical and subtropical Indo-Pacific. However, the genetic connectivity between populations of corals separated by large distances (thousands of kilometers) is not well known. We analyzed variability in the nucleotide sequence of the internal transcribed spacer-1 (ITS-1) of the nuclear ribosomal gene unit in the ubiquitous coral Stylophora pistillata, across the western Pacific Ocean. Eight populations from Japan, Malaysia, and the northern and southern Great Barrier Reef (GBR) were studied. Phylogenetic analyses and analysis of molecular variance (AMOVA) clearly revealed that there is panmixia among these coral populations. AMOVA showed that ITS-1 sequence variability was greater within populations (78.37%) than among populations (12.06%). These patterns strongly suggest high levels of connectivity across the species' latitudinal distribution range in the western Pacific, as is seen in many marine invertebrates.

Similarity between class A and class B G-protein-coupled receptors exemplified through calcitonin gene-related peptide receptor modelling and mutagenesis studies

Modelling class B G-protein-coupled receptors (GPCRs) using class A GPCR structural templates is difficult due to lack of homology. The plant GPCR, GCR1, has homology to both class A and class B GPCRs. We have used this to generate a class A-class B alignment, and by incorporating maximum lagged correlation of entropy and hydrophobicity into a consensus score, we have been able to align receptor transmembrane regions. We have applied this analysis to generate active and inactive homology models of the class B calcitonin gene-related peptide (CGRP) receptor, and have supported it with site-directed mutagenesis data using 122 CGRP receptor residues and 144 published mutagenesis results on other class B GPCRs. The variation of sequence variability with structure, the analysis of polarity violations, the alignment of group-conserved residues and the mutagenesis results at 27 key positions were particularly informative in distinguishing between the proposed and plausible alternative alignments. Furthermore, we have been able to associate the key molecular features of the class B GPCR signalling machinery with their class A counterparts for the first time. These include the [K/R]KLH motif in intracellular loop 1, [I/L]xxxL and KxxK at the intracellular end of TM5 and TM6, the NPXXY/VAVLY motif on TM7 and small group-conserved residues in TM1, TM2, TM3 and TM7. The equivalent of the class A DRY motif is proposed to involve Arg(2.39), His(2.43) and Glu(3.46), which makes a polar lock with T(6.37). These alignments and models provide useful tools for understanding class B GPCR function.

Detection of foveation windows and analysis of foveation sequences in congenital nystagmus

Congenital nystagmus (CN) is an ocular-motor disorder characterised by involuntary, conjugated ocular oscillations, that can arise since the first months of life. Pathogenesis of congenital nystagmus is still under investigation. In general, CN patients show a considerable decrease of their visual acuity: image fixation on the retina is disturbed by nystagmus continuous oscillations, mainly horizontal. However, image stabilisation is still achieved during the short periods in which eye velocity slows down while the target image is placed onto the fovea (called foveation intervals). To quantify the extent of nystagmus, eye movement recording are routinely employed, allowing physicians to extract and analyse nystagmus main features such as shape, amplitude and frequency. Using eye movement recording, it is also possible to compute estimated visual acuity predictors: analytical functions which estimates expected visual acuity using signal features such as foveation time and foveation position variability. Use of those functions add information to typical visual acuity measurement (e.g. Landolt C test) and could be a support for therapy planning or monitoring. This study focus on robust detection of CN patients' foveations. Specifically, it proposes a method to recognize the exact signal tracts in which a subject foveates, This paper also analyses foveation sequences. About 50 eyemovement recordings, either infrared-oculographic or electrooculographic, from different CN subjects were acquired. Results suggest that an exponential interpolation for the slow phases of nystagmus could improve foveation time computing and reduce influence of breaking saccades and data noise. Moreover a concise description of foveation sequence variability can be achieved using non-fitting splines. © 2009 Springer Berlin Heidelberg.

Selection of conserved epitopes from hepatitis C virus for pan-populational stimulation of T-cell responses

The hepatitis C virus (HCV) is able to persist as a chronic infection, which can lead to cirrhosis and liver cancer. There is evidence that clearance of HCV is linked to strong responses by CD8 cytotoxic T lymphocytes (CTLs), suggesting that eliciting CTL responses against HCV through an epitope-based vaccine could prove an effective means of immunization. However, HCV genomic plasticity as well as the polymorphisms of HLA I molecules restricting CD8 T-cell responses challenges the selection of epitopes for a widely protective vaccine. Here, we devised an approach to overcome these limitations. From available databases, we first collected a set of 245 HCV-specific CD8 T-cell epitopes, all known to be targeted in the course of a natural infection in humans. After a sequence variability analysis, we next identified 17 highly invariant epitopes. Subsequently, we predicted the epitope HLA I binding profiles that determine their potential presentation and recognition. Finally, using the relevant HLA I-genetic frequencies, we identified various epitope subsets encompassing 6 conserved HCV-specific CTL epitopes each predicted to elicit an effective T-cell response in any individual regardless of their HLA I background. We implemented this epitope selection approach for free public use at the EPISOPT web server. © 2013 Magdalena Molero-Abraham et al.

High natural gene expression variation in the reef-building coral Acropora millepora: potential for acclimative and adaptive plasticity

Background: Ecosystems worldwide are suffering the consequences of anthropogenic impact. The diverse ecosystem of coral reefs, for example, are globally threatened by increases in sea surface temperatures due to global warming. Studies to date have focused on determining genetic diversity, the sequence variability of genes in a species, as a proxy to estimate and predict the potential adaptive response of coral populations to environmental changes linked to climate changes. However, the examination of natural gene expression variation has received less attention. This variation has been implicated as an important factor in evolutionary processes, upon which natural selection can act. Results: We acclimatized coral nubbins from six colonies of the reef-building coral Acropora millepora to a common garden in Heron Island (Great Barrier Reef, GBR) for a period of four weeks to remove any site-specific environmental effects on the physiology of the coral nubbins. By using a cDNA microarray platform, we detected a high level of gene expression variation, with 17% (488) of the unigenes differentially expressed across coral nubbins of the six colonies (jsFDR-corrected, p < 0.01). Among the main categories of biological processes found differentially expressed were transport, translation, response to stimulus, oxidation-reduction processes, and apoptosis. We found that the transcriptional profiles did not correspond to the genotype of the colony characterized using either an intron of the carbonic anhydrase gene or microsatellite loci markers. Conclusion: Our results provide evidence of the high inter-colony variation in A. millepora at the transcriptomic level grown under a common garden and without a correspondence with genotypic identity. This finding brings to our attention the importance of taking into account natural variation between reef corals when assessing experimental gene expression differences. The high transcriptional variation detected in this study is interpreted and discussed within the context of adaptive potential and phenotypic plasticity of reef corals. Whether this variation will allow coral reefs to survive to current challenges remains unknown.

Correlation between pre-treatment quasispecies complexity and treatment outcome in chronic HCV genotype 3a

Pre-treatment HCV quasispecies complexity and diversity may predict response to interferon based anti-viral therapy. The objective of this study was to retrospectively (1) examine temporal changes in quasispecies prior to the start of therapy and (2) investigate extensively quasispecies evolution in a group of 10 chronically infected patients with genotype 3a, treated with pegylated alpha 2a-Interferon and ribavirin. The degree of sequence heterogeneity within the hypervariable region 1 was assessed by analyzing 20-30 individual clones in serial serum samples. Genetic parameters, including amino acid Shannon entropy, Hamming distance and genetic distance were calculated for each sample. Treatment outcome was divided into (1) sustained virological responders (SVR) and (2) treatment failure (TF).Our results indicate, (1) quasispecies complexity and diversity are lower in the SVR group, (2) quasispecies vary temporally and (3) genetic heterogeneity at baseline can be used to predict treatment outcome. We discuss the results from the perspective of replicative homeostasis. We discuss the results from the perspective of replicative homeostasis.