Insights on the HLA-G Evolutionary History Provided by a Nearby Alu Insertion

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

TNFRSF13B Genetic variability an anthropological - evolutionary approach to Biomedical Research

In the recent years TNFRSF13B coding variants have been implicated by clinical genetics studies in Common Variable Immunodeficiency (CVID), the most common clinically relevant primary immunodeficiency in individuals of European ancestry, but their functional effects in relation to the development of the disease have not been entirely established. To examine the potential contribution of such variants to CVID, the more comprehensive perspective of an evolutionary approach was applied in this study, underling the belief that evolutionary genetics methods can play a role in dissecting the origin, causes and diffusion of human diseases, representing a powerful tool also in human health research. For this purpose, TNFRSF13B coding region was sequenced in 451 healthy individuals belonging to 26 worldwide populations, in addition to 96 control, 77 CVID and 38 Selective IgA Deficiency (IgAD) individuals from Italy, leading to the first achievement of a global picture of TNFRSF13B nucleotide diversity and haplotype structure and making suggestion of its evolutionary history possible. A slow rate of evolution, within our species and when compared to the chimpanzee, low levels of genetic diversity geographical structure and the absence of recent population specific selective pressures were observed for the examined genomic region, suggesting that geographical distribution of its variability is more plausibly related to its involvement also in innate immunity rather than in adaptive immunity only. This, together with the extremely subtle disease/healthy samples differences observed, suggests that CVID might be more likely related to still unknown environmental and genetic factors, rather than to the nature of TNFRSF13B variants only.

Transposable elements dynamics in taxa with different reproductive strategies or speciation rate

In recent years the advances in genomics allowed to understand the importance of Transposable Elements (TE) in the evolution of eukaryotic genomes. In this thesis I face two aspects of the TE impact on the in the animal kingdom. The first part is a comparison of the dynamics of the TE dynamics in three species of stick-insects of the Genus Bacillus. I produced three random genomic libraries of 200 Kbps for the three parental species of the taxon: a gonochoric population of Bacillus rossius (facultative parthenogenetic), Bacillus grandii (gonochoric) and Bacillus atticus (obligate parthenogenetic). The unisexual taxon Bacillus atticus does not shows dramatic differences in TE total content and activity with respect to Bacillus grandii and Bacillus rossius. This datum does not confirm the trend observed in other animal models in which unisexual taxa tend to repress the activity of TE to escape the extinction by accumulation of harmful mutations. In the second part I tried to add a contribute to the debate initiated in recent years about the possibility that a high TE content is linked to a high rate of speciation. I designed an evolutionary framework to establish the different rate of speciation among two or more taxa, then I compared TE dynamics considering the different rates of speciation. The species dataset comprises: 29 mammals, four birds, two fish and two insects. On the whole the majority of comparisons confirms the expected trend. In particular the amount of species analyzed in Mammalia allowed me to get a statistical support (p<0,05) of the fact that the TE activity of recently mobilized elements is positively related with the rate of speciation.

Evolutionary forces shaping genomic islands of population differentiation in humans

Background Levels of differentiation among populations depend both on demographic and selective factors: genetic drift and local adaptation increase population differentiation, which is eroded by gene flow and balancing selection. We describe here the genomic distribution and the properties of genomic regions with unusually high and low levels of population differentiation in humans to assess the influence of selective and neutral processes on human genetic structure. Methods Individual SNPs of the Human Genome Diversity Panel (HGDP) showing significantly high or low levels of population differentiation were detected under a hierarchical-island model (HIM). A Hidden Markov Model allowed us to detect genomic regions or islands of high or low population differentiation. Results Under the HIM, only 1.5% of all SNPs are significant at the 1% level, but their genomic spatial distribution is significantly non-random. We find evidence that local adaptation shaped high-differentiation islands, as they are enriched for non-synonymous SNPs and overlap with previously identified candidate regions for positive selection. Moreover there is a negative relationship between the size of islands and recombination rate, which is stronger for islands overlapping with genes. Gene ontology analysis supports the role of diet as a major selective pressure in those highly differentiated islands. Low-differentiation islands are also enriched for non-synonymous SNPs, and contain an overly high proportion of genes belonging to the 'Oncogenesis' biological process. Conclusions Even though selection seems to be acting in shaping islands of high population differentiation, neutral demographic processes might have promoted the appearance of some genomic islands since i) as much as 20% of islands are in non-genic regions ii) these non-genic islands are on average two times shorter than genic islands, suggesting a more rapid erosion by recombination, and iii) most loci are strongly differentiated between Africans and non-Africans, a result consistent with known human demographic history.

Perinatal adaptation in mammals: the impact of metabolic rate

Mammalian birth is accompanied by profound changes in metabolic rate that can be described in terms of body size relationship (Kleiber's rule). Whereas the fetus, probably as an adaptation to the low intrauterine pO2, exhibits an "inappropriately" low, adult-like specific metabolic rate, the term neonate undergoes a rapid metabolic increase up to the level to be expected from body size. A similar, albeit slowed, "switching-on" of metabolic size allometry is found in human preterm neonates whereas animals that are normally born in a very immature state are able to retard or even suppress the postnatal metabolic increase in favor of weight gain and O2 supply. Moreover, small immature mammalian neonates exhibit a temporary oxyconforming behavior which enhances their hypoxia tolerance, yet is lost to the extent by which the size-adjusted metabolic rate is "locked" by increasing mitochondrial density. Beyond the perinatal period, there are no other deviations from metabolic size allometry among mammals except in hibernation where the temporary "switching-off" of Kleiber's rule is accompanied by a deep reduction in tissue pO2. This gives support to the hypothesis that the postnatal metabolic increase represents an "escape from oxygen" similar to the evolutionary roots of mitochondrial respiration, and that the overall increase in specific metabolic rate with decreasing size might contribute to prevent tissues from O2 toxicity.

Molluscan memory of injury: evolutionary insights into chronic pain and neurological disorders.

Molluscan preparations have yielded seminal discoveries in neuroscience, but the experimental advantages of this group have not, until now, been complemented by adequate molecular or genomic information for comparisons to genetically defined model organisms in other phyla. The recent sequencing of the transcriptome and genome of Aplysia californica, however, will enable extensive comparative studies at the molecular level. Among other benefits, this will bring the power of individually identifiable and manipulable neurons to bear upon questions of cellular function for evolutionarily conserved genes associated with clinically important neural dysfunction. Because of the slower rate of gene evolution in this molluscan lineage, more homologs of genes associated with human disease are present in Aplysia than in leading model organisms from Arthropoda (Drosophila) or Nematoda (Caenorhabditis elegans). Research has hardly begun in molluscs on the cellular functions of gene products that in humans are associated with neurological diseases. On the other hand, much is known about molecular and cellular mechanisms of long-term neuronal plasticity. Persistent nociceptive sensitization of nociceptors in Aplysia displays many functional similarities to alterations in mammalian nociceptors associated with the clinical problem of chronic pain. Moreover, in Aplysia and mammals the same cell signaling pathways trigger persistent enhancement of excitability and synaptic transmission following noxious stimulation, and these highly conserved pathways are also used to induce memory traces in neural circuits of diverse species. This functional and molecular overlap in distantly related lineages and neuronal types supports the proposal that fundamental plasticity mechanisms important for memory, chronic pain, and other lasting alterations evolved from adaptive responses to peripheral injury in the earliest neurons. Molluscan preparations should become increasingly useful for comparative studies across phyla that can provide insight into cellular functions of clinically important genes.

Eco-evolutionary Dynamics of Individual-Based Food Webs

The past decade has seen the rise of high resolution datasets. One of the main surprises of analysing such data has been the discovery of a large genetic, phenotypic and behavioural variation and heterogeneous metabolic rates among individuals within natural populations. A parallel discovery from theory and experiments has shown a strong temporal convergence between evolutionary and ecological dynamics, but a general framework to analyse from individual-level processes the convergence between ecological and evolutionary dynamics and its implications for patterns of biodiversity in food webs has been particularly lacking. Here, as a first approximation to take into account intraspecific variability and the convergence between the ecological and evolutionary dynamics in large food webs, we develop a model from population genomics and microevolutionary processes that uses sexual reproduction, genetic-distance-based speciation and trophic interactions. We confront the model with the prey consumption per individual predator, species-level connectance and prey–predator diversity in several environmental situations using a large food web with approximately 25,000 sampled prey and predator individuals. We show higher than expected diversity of abundant species in heterogeneous environmental conditions and strong deviations from the observed distribution of individual prey consumption (i.e. individual connectivity per predator) in all the environmental conditions. The observed large variance in individual prey consumption regardless of the environmental variability collapsed species-level connectance after small increases in sampling effort. These results suggest (1) intraspecific variance in prey–predator interactions has a strong effect on the macroscopic properties of food webs and (2) intraspecific variance is a potential driver regulating the speed of the convergence between ecological and evolutionary dynamics in species-rich food webs. These results also suggest that genetic–ecological drift driven by sexual reproduction, equal feeding rate among predator individuals, mutations and genetic-distance-based speciation can be used as a neutral food web dynamics test to detect the ecological and microevolutionary processes underlying the observed patterns of individual and species-based food webs at local and macroecological scales.

Experimental Evidence of an Eco-evolutionary Feedback during Adaptive Divergence

Differences in how organisms modify their environment can evolve rapidly and might influence adaptive population divergence [1, 2]. In a common garden experiment in aquatic mesocosms, we found that adult stickleback from a recently diverged pair of lake and stream populations had contrasting effects on ecosystem metrics. These modifications were caused by both genetic and plastic differences between populations and were sometimes comparable in magnitude to those caused by the presence/ absence of stickleback. Lake and streamfish differentially affected the biomass of zooplankton and phytoplankton, the concentration of phosphorus, and the abundance of several prey (e.g., copepods) and non-prey (e.g., cyanobacteria) species. The adult mediated effects on mesocosm ecosystems influenced the survival and growth of a subsequent generation of juvenile stickleback reared in the same mesocosms. The prior presence of adults decreased the overall growth rate of juveniles, and the prior presence of stream adults lowered overall juvenile survival. Among the survivors, lake juveniles grew faster than co-occurring stream juveniles, except in mesocosm ecosystems previously modified by adult lake fish that were reared on plankton. Overall, our results provide evidence for reciprocal interactions between ecosystem dynamics and evolutionary change (i.e., eco-evolutionary feedbacks) in the early stages of adaptive population divergence.

STUDIES ON THE PATTERNS OF NUCLEOTIDE AND AMINO ACID SUBSTITUTION (ELECTROPHORESIS, MUTATION RATE, SELECTION, BASE, COMPOSITION)

Theoretical and empirical studies were conducted on the pattern of nucleotide and amino acid substitution in evolution, taking into account the effects of mutation at the nucleotide level and purifying selection at the amino acid level. A theoretical model for predicting the evolutionary change in electrophoretic mobility of a protein was also developed by using information on the pattern of amino acid substitution. The specific problems studied and the main results obtained are as follows: (1) Estimation of the pattern of nucleotide substitution in DNA nuclear genomes. The pattern of point mutations and nucleotide substitutions among the four different nucleotides are inferred from the evolutionary changes of pseudogenes and functional genes, respectively. Both patterns are non-random, the rate of change varying considerably with nucleotide pair, and that in both cases transitions occur somewhat more frequently than transversions. In protein evolution, substitution occurs more often between amino acids with similar physico-chemical properties than between dissimilar amino acids. (2) Estimation of the pattern of nucleotide substitution in RNA genomes. The majority of mutations in retroviruses accumulate at the reverse transcription stage. Selection at the amino acid level is very weak, and almost non-existent between synonymous codons. The pattern of mutation is very different from that in DNA genomes. Nevertheless, the pattern of purifying selection at the amino acid level is similar to that in DNA genomes, although selection intensity is much weaker. (3) Evaluation of the determinants of molecular evolutionary rates in protein-coding genes. Based on rates of nucleotide substitution for mammalian genes, the rate of amino acid substitution of a protein is determined by its amino acid composition. The content of glycine is shown to correlate strongly and negatively with the rate of substitution. Empirical formulae, called indices of mutability, are developed in order to predict the rate of molecular evolution of a protein from data on its amino acid sequence. (4) Studies on the evolutionary patterns of electrophoretic mobility of proteins. A theoretical model was constructed that predicts the electric charge of a protein at any given pH and its isoelectric point from data on its primary and quaternary structures. Using this model, the evolutionary change in electrophoretic mobilities of different proteins and the expected amount of electrophoretically hidden genetic variation were studied. In the absence of selection for the pI value, proteins will on the average evolve toward a mildly basic pI. (Abstract shortened with permission of author.) ^

Two hemoglobin genes in Arabidopsis thaliana: The evolutionary origins of leghemoglobins

We cloned two hemoglobin genes from Arabidopsis thaliana. One gene, AHB1, is related in sequence to the family of nonsymbiotic hemoglobin genes previously identified in a number of plant species (class 1). The second hemoglobin gene, AHB2, represents a class of nonsymbiotic hemoglobin (class 2) related in sequence to the symbiotic hemoglobin genes of legumes and Casuarina. The properties of these two hemoglobins suggest that the two families of nonsymbiotic hemoglobins may differ in function from each other and from the symbiotic hemoglobins. AHB1 is induced, in both roots and rosette leaves, by low oxygen levels. Recombinant AHB1 has an oxygen affinity so high as to make it unlikely to function as an oxygen transporter. AHB2 is expressed at a low level in rosette leaves and is low temperature-inducible. AHB2 protein has a lower affinity for oxygen than AHB1 but is similar to AHB1 in having an unusually low, pH-sensitive oxygen off-rate.

Social status regulates growth rate: Consequences for life-history strategies

The life-history strategies of organisms are sculpted over evolutionary time by the relative prospects of present and future reproductive success. As a consequence, animals of many species show flexible behavioral responses to environmental and social change. Here we show that disruption of the habitat of a colony of African cichlid fish, Haplochromis burtoni (Günther) caused males to switch social status more frequently than animals kept in a stable environment. H. burtoni males can be either reproductively active, guarding a territory, or reproductively inactive (nonterritorial). Although on average 25–50% of the males are territorial in both the stable and unstable environments, during the 20-week study, nearly two-thirds of the animals became territorial for at least 1 week. Moreover, many fish changed social status several times. Surprisingly, the induced changes in social status caused changes in somatic growth. Nonterritorial males and animals ascending in social rank showed an increased growth rate whereas territorial males and animals descending in social rank slowed their growth rate or even shrank. Similar behavioral and physiological changes are caused by social change in animals kept in stable environmental conditions, although at a lower rate. This suggests that differential growth, in interaction with environmental conditions, is a central mechanism underlying the changes in social status. Such reversible phenotypic plasticity in a crucial life-history trait may have evolved to enable animals to shift resources from reproduction to growth or vice versa, depending on present and future reproductive prospects.

The diversity and evolutionary relationships of the pregnancy-associated glycoproteins, an aspartic proteinase subfamily consisting of many trophoblast-expressed genes

The pregnancy-associated glycoproteins (PAGs) are structurally related to the pepsins, thought to be restricted to the hooved (ungulate) mammals and characterized by being expressed specifically in the outer epithelial cell layer (chorion/trophectoderm) of the placenta. At least some PAGs are catalytically inactive as proteinases, although each appears to possess a cleft capable of binding peptides. By cloning expressed genes from ovine and bovine placental cDNA libraries, by Southern genomic blotting, by screening genomic libraries, and by using PCR to amplify portions of PAG genes from genomic DNA, we estimate that cattle, sheep, and most probably all ruminant Artiodactyla possess many, possibly 100 or more, PAG genes, many of which are placentally expressed. The PAGs are highly diverse in sequence, with regions of hypervariability confined largely to surface-exposed loops. Nonsynonymous (replacement) mutations in the regions of the genes coding for these hypervariable loop segments have accumulated at a higher rate than synonymous (silent) mutations. Construction of distance phylograms, based on comparisons of PAG and related aspartic proteinase amino acid sequences, suggests that much diversification of the PAG genes occurred after the divergence of the Artiodactyla and Perissodactyla, but that at least one gene is represented outside the hooved species. The results also suggest that positive selection of duplicated genes has acted to provide considerable functional diversity among the PAGs, whose presence at the interface between the placenta and endometrium and in the maternal circulation indicates involvement in fetal–maternal interactions.

The heme-copper oxidases of Thermus thermophilus catalyze the reduction of nitric oxide: Evolutionary implications

We show that the heme-copper terminal oxidases of Thermus thermophilus (called ba3 and caa3) are able to catalyze the reduction of nitric oxide (NO) to nitrous oxide (N2O) under reducing anaerobic conditions. The rate of NO consumption and N2O production were found to be linearly dependent on enzyme concentration, and activity was abolished by enzyme denaturation. Thus, contrary to the eukaryotic enzyme, both T. thermophilus oxidases display a NO reductase activity (3.0 ± 0.7 mol NO/mol ba3 × min and 32 ± 8 mol NO/mol caa3 × min at [NO] ≈ 50 μM and 20°C) that, though considerably lower than that of bona fide NO reductases (300–4,500 mol NO/mol enzyme × min), is definitely significant. We also show that for ba3 oxidase, NO reduction is associated to oxidation of cytochrome b at a rate compatible with turnover, suggesting a mechanism consistent with the stoichiometry of the overall reaction. We propose that the NO reductase activity of T. thermophilus oxidases may depend on a peculiar CuB+ coordination, which may be revealed by the forthcoming three-dimensional structure. These findings support the hypothesis of a common phylogeny of aerobic respiration and bacterial denitrification, which was proposed on the basis of structural similarities between the Pseudomonas stutzeri NO reductase and the cbb3 terminal oxidases. Our findings represent functional evidence in support of this hypothesis.

Evolutionary instability of the major histocompatibility complex class I loci in New World primates

Homologues of the human major histocompatibility complex (MHC) HLA-A, -B, -E, -F, and -G loci are present in all the Catarrhini (Old World primates, apes, and humans), and some of their allelic lineages have survived several speciation events. Analysis of 26 MHC class I cDNAs from seven different genera of New World primates revealed that the Callitrichinae (tamarins and marmosets) are an exception to these rules of MHC stability. In gene trees of primate MHC class I genes, sequences from the Callitrichinae cluster in a genus-specific fashion, whereas in the other genera of New World primates, as in the Catarrhini, they cluster in a transgeneric way. The genus-specific clustering of the Callitrichinae cDNAs indicates that there is no orthology between MHC class I loci in genera of this phyletic group. Additionally, the Callitrichinae genera exhibit limited variability of their MHC class I genes, in contrast to the high variability displayed by all other primates. Each Callitrichinae genus, therefore, expresses its own set of MHC class I genes, suggesting that an unusually high rate of turnover of loci occurs in this subfamily. The limited variability of MHC class I genes in the Callitrichinae is likely the result of the recent origin of these loci.

Mutation rate varies among alleles at a microsatellite locus: Phylogenetic evidence

The understanding of the mutational mechanism that generates high levels of variation at microsatellite loci lags far behind the application of these genetic markers. A phylogenetic approach was developed to study the pattern and rate of mutations at a dinucleotide microsatellite locus tightly linked to HLA-DQB1 (DQCAR). A random Japanese population (n = 129) and a collection of multiethnic samples (n = 941) were typed at the DQB1 and DQCAR loci. The phylogeny of DQB1 alleles was then reconstructed and DQCAR alleles were superimposed onto the phylogeny. This approach allowed us to group DQCAR alleles that share a common ancestor. The results indicated that the DQCAR mutation rate varies drastically among alleles within this single microsatellite locus. Some DQCAR alleles never mutated during a long period of evolutionary time. Sequencing of representative DQCAR alleles showed that these alleles lost their ability to mutate because of nucleotide substitutions that shorten the length of uninterrupted CA repeat arrays; in contrast, all mutating alleles had relatively longer perfect CA repeat sequences.