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.

Detection and sequence analysis of a nickel-induced mutation in a retroviral model system

We have developed a novel way to assess the mutagenicity of environmentally important metal carcinogens, such as nickel, by creating a positive selection system based upon the conditional expression of a retroviral transforming gene. The target gene is the v-mos gene in MuSVts110, a murine retrovirus possessing a growth temperature dependent defect in expression of the transforming gene due to viral RNA splicing. In normal rat kidney cells infected with MuSVts110 (6m2 cells), splicing of the MuSVts110 RNA to form the mRNA from which the transforming protein, p85$\sp{\rm gag-mos}$, is translated is growth-temperature dependent, occurring at 33 C and below but not at 39 C and above. This splicing "defect" is mediated by cis-acting viral sequences. Nickel chloride treatment of 6m2 cells followed by growth at 39 C, allowed the selection of "revertant" cells which constitutively express p85$\sp{\rm gag-mos}$ due to stable changes in the viral RNA splicing phenotype, suggesting that nickel, a carcinogen whose mutagenicity has not been well established, could induce mutations in mammalian genes. We also show by direct sequencing of PCR-amplified integrated MuSVts110 DNA from a 6m2 nickel-revertant cell line that the nickel-induced mutation affecting the splicing phenotype is a cis-acting 70-base duplication of a region of the viral DNA surrounding the 3$\sp\prime$ splice site. These findings provide the first example of the molecular basis for a nickel-induced DNA lesion and establish the mutagenicity of this potent carcinogen. ^

Adaptive Molecular Evolution in the Opsin Genes of Rapidly Speciating Cichlid Species

Cichlid fish inhabit a diverse range of environments that vary in the spectral content of light available for vision. These differences should result in adaptive selective pressure on the genes involved in visual sensitivity, the opsin genes. This study examines the evidence for differential adaptive molecular evolution in East African cichlid opsin genes due to gross differences in environmental light conditions. First, we characterize the selective regime experienced by cichlid opsin genes using a likelihood ratio test format, comparing likelihood models with different constraints on the relative rates of amino acid substitution, across sites. Second, we compare turbid and clear lineages to determine if there is evidence of differences in relative rates of substitution. Third, we present evidence of functional diversification and its relationship to the photic environment among cichlid opsin genes. We report statistical evidence of positive selection in all cichlid opsin genes, except short wavelength–sensitive 1 and short wavelength–sensitive 2b. In all genes predicted to be under positive selection, except short wavelength–sensitive 2a, we find differences in selective pressure between turbid and clear lineages. Potential spectral tuning sites are variable among all cichlid opsin genes; however, patterns of substitution consistent with photic environment–driven evolution of opsin genes are observed only for short wavelength–sensitive 1 opsin genes. This study identifies a number of promising candidate-tuning sites for future study by site-directed mutagenesis. This work also begins to demonstrate the molecular evolutionary dynamics of cichlid visual sensitivity and its relationship to the photic environment.

Identification of cadmium-induced mutations in a mammalian cell line

Epidemiological studies have shown cadmium to induce cancer in humans, while experimental studies have proven this metal to be a potent tumor inducer in animals. However, cadmium appears nonmutagenic in most prokaryotic and eukaryotic mutagenesis assays. In this study, we present the identification of mutations in normal rat kidney cells infected with the mutant MuSVts110 retrovirus (6m2 cells) as a result of treatment with cadmium chloride. The detection of these mutations was facilitated by the use of a novel mutagenesis assay established in this laboratory. The 6m2 reversion assay is a positive selection system based on the conditional expression of the MuSVts110 v-mos gene. In MuSVts110 the gag and mos genes are fused out of frame, thus the translation of the v-mos sequence requires a frameshift in the genomic RNA. In 6m2 cells this frameshift is accomplished by the temperature-dependent splicing of the primary MuSVts110 transcript. Splicing of MuSVts110, which is mediated by cis-acting sequences, occurs when 6m2 cells are grown at 33$\sp\circ$C and below, but not at 39$\sp\circ$C. Therefore, 6m2 cells appear transformed at low growth temperatures, but take on a morphologically normal appearance when grown at high temperatures. The treatment of 6m2 cells with cadmium chloride resulted in the outgrowth of a number of cells that reverted to the transformed state at high growth temperatures. Analysis of the viral proteins expressed in these cadmium-induced 6m2 revertants suggested that they contained mutations in their MuSVts110 DNA. Sequencing of the viral DNA from three revertants that constitutively expressed the P85$\sp{gag{-}mos}$ transforming protein revealed five different mutations. The Cd-B2 revertant contained three of those mutations: an A-to-G transition 48 bases downstream of the MuSVts110 3$\sp\prime$ splice site, plus a G-to-T and an A-to-T transversion 84 and 100 bases downstream of the 5$\sp\prime$ splice site, respectively. The Cd-15-5 revertant also contained a point mutation, a T-to-C transition 46 bases downstream of the 5$\sp\prime$ splice site, while Cd-10-5 contained a three base deletion of MuSVts110 11 bases upstream of the 3$\sp\prime$ splice site. A fourth revertant, Cd-10, expressed a P100$\sp{gag{-}mos}$ transforming protein, and was found to have a two base deletion. This deletion accomplished the frameshift necessary for v-mos expression, but did not alter MuSVts110 RNA splicing and the expression of p85$\sp{gag{-}mos}.$ Lastly, sequencing of the MuSVts110 DNA from three spontaneous revertants revealed the same G to T transversion in each one. This was the same mutation that was found in the Cd-B2 revertant. These findings provide the first example of mutations resulting from exposure to cadmium and suggest, by the difference in each mutation, the complexity of the mechanism utilized by cadmium to induce DNA damage. ^

Variation in life-history traits and their plasticities to elevational transplantation among seed families suggests potential for adaptative evolution of 15 tropical plant species to climate change

• Premise of the study: Because not all plant species will be able to move in response to global warming, adaptive evolution matters largely for plant persistence. As prerequisites for adaptive evolution, genetic variation in and selection on phenotypic traits are needed, but these aspects have not been studied in tropical species. We studied how plants respond to transplantation to different elevations on Mt. Kilimanjaro, Tanzania, and whether there is quantitative genetic (among-seed family) variation in and selection on life-history traits and their phenotypic plasticity to the different environments. • Methods: We reciprocally transplanted seed families of 15 common tropical, herbaceous species of the montane and savanna vegetation zone at Mt. Kilimanjaro to a watered experimental garden in the montane (1450 m) and in the savanna (880 m) zone at the mountain’s slope and measured performance, reproductive, and phenological traits. • Results: Plants generally performed worse in the savanna garden, indicating that the savanna climate was more stressful and thus that plants may suffer from future climate warming. We found significant quantitative genetic variation in all measured performance and reproductive traits in both gardens and for several measures of phenotypic plasticity in response to elevational transplantation. Moreover, we found positive selection on traits at low and intermediate trait values levelling to neutral or negative selection at high values. • Conclusions: We conclude that common plants at Mt. Kilimanjaro express quantitative genetic variation in fitness-relevant traits and in their plasticities, suggesting potential to adapt evolutionarily to future climate warming and increased temperature variability.

Prehistoric genomes reveal the genetic foundation and cost of horse domestication

The domestication of the horse revolutionized warfare, trade, and the exchange of people and ideas. This at least 5,500-y-long process, which ultimately transformed wild horses into the hundreds of breeds living today, is difficult to reconstruct from archeological data and modern genetics alone. We therefore sequenced two complete horse genomes, predating domestication by thousands of years, to characterize the genetic footprint of domestication. These ancient genomes reveal predomestic population structure and a significant fraction of genetic variation shared with the domestic breeds but absent from Przewalski’s horses. We find positive selection on genes involved in various aspects of locomotion, physiology, and cognition. Finally, we show that modern horse genomes contain an excess of deleterious mutations, likely representing the genetic cost of domestication.

Evolutionary history and global spread of the Mycobacterium tuberculosis Beijing lineage.

Mycobacterium tuberculosis strains of the Beijing lineage are globally distributed and are associated with the massive spread of multidrug-resistant (MDR) tuberculosis in Eurasia. Here we reconstructed the biogeographical structure and evolutionary history of this lineage by genetic analysis of 4,987 isolates from 99 countries and whole-genome sequencing of 110 representative isolates. We show that this lineage initially originated in the Far East, from where it radiated worldwide in several waves. We detected successive increases in population size for this pathogen over the last 200 years, practically coinciding with the Industrial Revolution, the First World War and HIV epidemics. Two MDR clones of this lineage started to spread throughout central Asia and Russia concomitantly with the collapse of the public health system in the former Soviet Union. Mutations identified in genes putatively under positive selection and associated with virulence might have favored the expansion of the most successful branches of the lineage.

Implications of Longterm Diameter-Limit Harvesting: Effects on Radial Growth of Red Spruce (Picea rubens) and Genetic Diversity of White Pine (Pinus strobus)

For over 3 centuries, diameter-limit harvesting has been a predominant logging method in the northeastern United States. Silvicultural theory asserts that such intensively selective harvesting can lead to genetic degradation. A decrease in softwood productivity has recently been reported in Maine - has a long history of dysgenic selection degraded the genetic resources of Maine softwoods, contributing to a decrease in growth and productivity? This study examines two aspects of potential implications of diameter-limit harvesting: effects on residual phenotypes of red spruce and impacts on genetic diversity of white pine. Radial growth of residual red spruce trees in stands experiencing 50 years of fixed diameter-limit harvesting was measured using annual increment rings and compared with residual red spruce trees in positive selection stands. Trees remaiaing after several rounds of diameter-limit harvesting exhibited sigdicantl y smaller radial sizes throughout their lives, and displayed significantly slower growth rates for the first 80 years of measured growth. These results strongly suggest that the largest and fastest-growing genotypes and their respective gene complexes determining good radial growth have been removed from the diameter-limit stand. Dysgenic selection can be observed in fixed diarneter-limit stands, resulting in a diminished genetic resource and decreased residual stand value. To examine more direct genetic implications of long-term diameter-limit harvesting, microsatellite DNA markers were implemented to study genetic diversity of eastern white pine in Maine. Three age groups of trees were studied: mature trees older than 200 years, juvenile trees 5-30 years old, and embryos. Trees were genotyped at 10 microsatellite loci. Overall genetic diversity levels of eastern white pine in Maine were extremely high, with an average observed heterozygosity of 0.762. Genetic differentiation was minimal among and between all three age groups, although an excess of heterozygotes was shown in the mature and juvenile groups that was not reflected in the embryo group, which actually had a slight heterozygote deficiency. Allele frequencies did not differ significantly between age groups, but did reveal more rare and low frequency alleles in the embryo groups than in the mature group. Overall, low frequency alleles comprise the largest portion of alleles in the sample population, with no common alleles evident overall. These results suggest that significant genetic degradation has either not occurred for white pine, or that the results of dysgenic selection have not yet emerged. It is clear, however, that selective harvesting could result in a loss of low frequency alleles, which are a primary reserve of evolutionary potential in a species. Implications of these studies affect industrial forestry, regional economics, and ecological concerns for the northeast. Long-term diameter-limit harvesting can lead to a degradation of residual phenotypes, and an overall decrease in stand quality. Potentially, a loss of low frequency, locally adapted alleles could result in a decrease of allelic richness and degradation of the regidnal genetic resource. Decreased genetic variation can lead to seriously limited evolutionary potential of species and ecosystems, particularly in rapidly changing environments. Based on these findings, I recommend a reassessment of any harvesting prescription that includes fixed diameter-limit removals, particularly for species that have low natural genetic diversity levels or a limited natural range, such as red spruce. Maintenance of a healthy genetic reserve can avoid effects of dysgenic harvesting.

Characterization of anti -sulfatide autoantibodies in demyelinating diseases

Multiple sclerosis (MS) is the most common autoimmune disease of the central nerve system and Guillain Barré Syndrome (GBS) is an inflammatory neuropathy involving the peripheral nerves. Anti-myelin immunoglobins may play a role in the demyelination processes of the both diseases. Sulfatide is an abundant glycolipid on myelin and is a candidate target antigen for disease related autoantibodies. The objective of this study was to characterize anti-sulfatide antibodies and compare antibodies from GBS and MS patients with fetal antibodies. Our hypothesis is that some B cells producing disease-associated autoantibodies are derived from or related to B cells of the fetal repertoire. Here we report that reactivity of plasma IgM against sulfatide was elevated in twelve MS patients compared with twelve normal subjects. This result implies that anti-sulfatide antibodies are disease-related. A total of sixteen human B lymphocyte clones producing anti-sulfatide autoantibodies were isolated from MS patients, GBS patients and a human fetus. Seven of the clones were from three MS patients, four of the clones were from three GBS patients and five were from the spleen of a twenty-week human fetus. Sequences have been obtained for the heavy and light chain variable regions (VDJ and VJ regions) of all of the anti-sulfatide immunoglobulins. Seven of the sixteen antibodies used VH3 for the variable region gene of the heavy chain consistent with the rate of VH3 usage in randomly selected B cells. Somatic mutations were significantly more frequent in the patient antibodies than in the fetus and somatic mutations in CDR's (Complementarity Determining Region) were significantly more frequent than in framework regions. No significant difference was found between patients and fetus in length of VH CDRIII. However, it is reported that antibodies from randomly selected normal adult B cells have longer CDRIII lengths than those of the fetus (Sanz I, 1991 Journal of Immunology Sep 1;147(5):1720-9). Our results are consistent with derivation of the precursors of B cells producing these autoantibodies from B cells related to those of the fetal repertoire. These findings are consistent with a model in which quiescent B cells from clones produced early in development undergo proliferation in dysregulated disease states, accumulating somatic mutations and increasing in reactivity toward self-antigens. ^ Epitope mapping and molecular modeling were done to elucidate the relationships between antibody structure and binding characteristics. The autoantibodies were tested for binding activity to three different antigens: sulfatide, galactoceramide and ceramide. Molecular modeling suggests that antibodies with positive charge surrounded by or adjacent to hydrophobic groups in the binding pocket bind to the head of sulfatide via the sulfate group through electrostatic interactions. However, the antibodies with hydrophobic groups separated from positive charges appear to bind to the hydrophobic tail of sulfatide. This observation was supported by a study of the effect of NaCl concentration on antigen binding. The result suggested that electrostatic interactions played a major role in sulfate group binding and that hydrophobic interactions were of greater importance for binding to the ceramide group. Our three-dimensional structure data indicated that epitope specificity of these antibodies is more predictable at the level of tertiary than primary structure and suggested positive selection based on structure occurred in the. formation of those autoantibodies. ^

DNA alignment and Bayesian trees of partial sequences of COI, COIII and the haemocyanin functional unit f-g of 28 octopod species

Background: Octopods have successfully colonised the world's oceans from the tropics to the poles. Yet, successful persistence in these habitats has required adaptations of their advanced physiological apparatus to compensate impaired oxygen supply. Their oxygen transporter haemocyanin plays a major role in cold tolerance and accordingly has undergone functional modifications to sustain oxygen release at sub-zero temperatures. However, it remains unknown how molecular properties evolved to explain the observed functional adaptations. We thus aimed to assess whether natural selection affected molecular and structural properties of haemocyanin that explains temperature adaptation in octopods. Results: Analysis of 239 partial sequences of the haemocyanin functional units (FU) f and g of 28 octopod species of polar, temperate, subtropical and tropical origin revealed natural selection was acting primarily on charge properties of surface residues. Polar octopods contained haemocyanins with higher net surface charge due to decreased glutamic acid content and higher numbers of basic amino acids. Within the analysed partial sequences, positive selection was present at site 2545, positioned between the active copper binding centre and the FU g surface. At this site, methionine was the dominant amino acid in polar octopods and leucine was dominant in tropical octopods. Sites directly involved in oxygen binding or quaternary interactions were highly conserved within the analysed sequence. Conclusions: This study has provided the first insight into molecular and structural mechanisms that have enabled octopods to sustain oxygen supply from polar to tropical conditions. Our findings imply modulation of oxygen binding via charge-charge interaction at the protein surface, which stabilize quaternary interactions among functional units to reduce detrimental effects of high pH on venous oxygen release. Of the observed partial haemocyanin sequence, residue 2545 formed a close link between the FU g surface and the active centre, suggesting a role as allosteric binding site. The prevalence of methionine at this site in polar octopods, implies regulation of oxygen affinity via increased sensitivity to allosteric metal binding. High sequence conservation of sites directly involved in oxygen binding indicates that functional modifications of octopod haemocyanin rather occur via more subtle mechanisms, as observed in this study.

All males are not created equal: Fertility differences depend on gamete recognition polymorphisms in sea urchins

Behaviors, morphologies, and genetic loci directly involved in reproduction have been increasingly shown to be polymorphic within populations. Explaining how such variants are maintained by selection is crucial to understanding the genetic basis of fertility differences, but direct tests of how alleles at reproductive loci affect fertility are rare. In the sea urchin genus Echinometra, the protein bindin mediates sperm attachment to eggs, evolves quickly, and is polymorphic within species. Eggs exposed to experimental sperm mixtures show strong discrimination on the basis of the males’ bindin genotype. Different females produce eggs that nonrandomly select sperm from different males, showing that variable egg–sperm interactions determine fertility. Eggs select sperm with a bindin genotype similar to their own, suggesting strong linkage between female choice and male trait loci. These experiments demonstrate that alleles at a single locus can have a strong effect on fertilization and that reproductive loci may retain functional polymorphisms through epistatic interactions between male and female traits. They also suggest that positive selection at gamete recognition loci like bindin involves strong selection within species on mate choice interactions.

Self-antigen does not accelerate immature B cell apoptosis, but stimulates receptor editing as a consequence of developmental arrest

In pre-B lymphocytes, productive rearrangement of Ig light chain genes allows assembly of the B cell receptor (BCR), which selectively promotes further developmental maturation through poorly defined transmembrane signaling events. Using a novel in vitro system to study immune tolerance during development, we find that BCR reactivity to auto-antigen blocks this positive selection, preventing down-regulation of light chain gene recombination and promoting secondary light chain gene rearrangements that often alter BCR specificity, a process called receptor editing. Under these experimental conditions, self-antigen induces secondary light chain gene rearrangements in at least two-thirds of autoreactive immature B cells, but fails to accelerate cell death at this stage. These data suggest that in these cells the mechanism of immune tolerance is receptor selection rather than clonal selection.

Efficient allelic exchange and transposon mutagenesis in Mycobacterium tuberculosis

A better understanding of Mycobacterium tuberculosis virulence mechanisms is highly dependent on the design of efficient mutagenesis systems. A system enabling the positive selection of insertional mutants having lost the delivery vector was developed. It uses ts-sacB vectors, which combine the counterselective properties of the sacB gene and a mycobacterial thermosensitive origin of replication and can therefore be efficiently counterselected on sucrose at 39°C. This methodology allowed the construction of M. tuberculosis transposition mutant libraries. Greater than 106 mutants were obtained, far exceeding the number theoretically required to obtain at least one insertion in every nonessential gene. This system is also efficient for gene exchange mutagenesis as demonstrated with the purC gene: 100% of the selected clones were allelic exchange mutants. Therefore, a single, simple methodology has enabled us to develop powerful mutagenesis systems, the lack of which was a major obstacle to the genetic characterization of M. tuberculosis.

Ecological factors rather than temporal factors dominate the evolution of vesicular stomatitis virus

Vesicular stomatitis New Jersey virus (VSV-NJ) is a rhabdovirus that causes economically important disease in cattle and other domestic animals in endemic areas from southeastern United States to northern South America. Its negatively stranded RNA genome is capable of undergoing rapid evolution, which allows phylogenetic analysis and molecular epidemiology studies to be performed. Previous epidemiological studies in Costa Rica showed the existence of at least two distinct ecological zones of high VSV-NJ activity, one located in the highlands (premontane tropical moist forest) and the other in the lowlands (tropical dry forest). We wanted to test the hypothesis that the viruses circulating in these ecological zones were genetically distinct. For this purpose, we sequenced the hypervariable region of the phosphoprotein gene for 50 VSV-NJ isolates from these areas. Phylogenetic analysis showed that viruses from each ecological zone had distinct genotypes. These genotypes were maintained in each area for periods of up to 8 years. This evolutionary pattern of VSV-NJ suggests an adaptation to ecological factors that could exert selective pressure on the virus. As previous data indicated an absence of virus adaptation to factors related to the bovine host (including immunological pressure), it appears that VSV genetic divergence represents positive selection to adapt to specific vectors and/or reservoirs at each ecological zone.

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.