Identifying adaptive genetic divergence among populations from genome scans

The identification of signatures of natural selection in genomic surveys has become an area of intense research, stimulated by the increasing ease with which genetic markers can be typed. Loci identified as subject to selection may be functionally important, and hence (weak) candidates for involvement in disease causation. They can also be useful in determining the adaptive differentiation of populations, and exploring hypotheses about speciation. Adaptive differentiation has traditionally been identified from differences in allele frequencies among different populations, summarised by an estimate of F-ST. Low outliers relative to an appropriate neutral population-genetics model indicate loci subject to balancing selection, whereas high outliers suggest adaptive (directional) selection. However, the problem of identifying statistically significant departures from neutrality is complicated by confounding effects on the distribution of F-ST estimates, and current methods have not yet been tested in large-scale simulation experiments. Here, we simulate data from a structured population at many unlinked, diallelic loci that are predominantly neutral but with some loci subject to adaptive or balancing selection. We develop a hierarchical-Bayesian method, implemented via Markov chain Monte Carlo (MCMC), and assess its performance in distinguishing the loci simulated under selection from the neutral loci. We also compare this performance with that of a frequentist method, based on moment-based estimates of F-ST. We find that both methods can identify loci subject to adaptive selection when the selection coefficient is at least five times the migration rate. Neither method could reliably distinguish loci under balancing selection in our simulations, even when the selection coefficient is twenty times the migration rate.

Optical genetic mapping defines regions of chromosomal variation in serovars of S. enterica subsp enterica of concern for human and animal health

Infections involving Salmonella enterica subsp. enterica serovars have serious animal and human health implications; causing gastroenteritis in humans and clinical symptoms, such as diarrhoea and abortion, in livestock. In this study an optical genetic mapping technique was used to screen 20 field isolate strains from four serovars implicated in disease outbreaks. The technique was able to distinguish between the serovars and the available sequenced strains and group them in agreement with similar data from microarrays and PFGE. The optical maps revealed variation in genome maps associated with antimicrobial resistance and prophage content in S. Typhimurium, and separated the S. Newport strains into two clear geographical lineages defined by the presence of prophage sequences. The technique was also able to detect novel insertions that may have had effects on the central metabolism of some strains. Overall optical mapping allowed a greater level of differentiation of genomic content and spatial information than more traditional typing methods.

Quantitative genetic divergence and standing genetic (co)variance in thermal reaction norms along latitude

Although the potential to adapt to warmer climate is constrained by genetic trade-offs, our understanding of how selection and mutation shape genetic (co)variances in thermal reaction norms is poor. Using 71 isofemale lines of the fly Sepsis punctum, originating from northern, central, and southern European climates, we tested for divergence in juvenile development rate across latitude at five experimental temperatures. To investigate effects of evolutionary history in different climates on standing genetic variation in reaction norms, we further compared genetic (co)variances between regions. Flies were reared on either high or low food resources to explore the role of energy acquisition in determining genetic trade-offs between different temperatures. Although the latter had only weak effects on the strength and sign of genetic correlations, genetic architecture differed significantly between climatic regions, implying that evolution of reaction norms proceeds via different trajectories at high latitude versus low latitude in this system. Accordingly, regional genetic architecture was correlated to region-specific differentiation. Moreover, hot development temperatures were associated with low genetic variance and stronger genetic correlations compared to cooler temperatures. We discuss the evolutionary potential of thermal reaction norms in light of their underlying genetic architectures, evolutionary histories, and the materialization of trade-offs in natural environments.

Comparative Genetic Diversity of Wild and Captive Populations of the Bare-Faced Curassow (Crax fasciolata) Based on Cross-Species Microsatellite Markers: Implications for Conservation and Management

The bare-faced curassow (Crax fasciolata) is a large Neotropical bird that suffers anthropogenic pressure across much of its range. A captive population is maintained for conservation management, although there has been no genetic screening of stocks. Based on the six microsatellite markers developed for Crax globulosa, the genetic variability of C. fasciolata and possible differences between a wild and a captive population were investigated. Only three loci were polymorphic, with a total of 27 alleles. More than half of these alleles were private to the wild (n = 8) or captive (n = 7) populations. Significant deviations from Hardy-Weinberg equilibrium were restricted to the captive population. Despite the number of private alleles, genetic drift has probably promoted differentiation between populations. Our results indicate that wild C. fasciolata populations are genetically impoverished and structured, but species-specific microsatellite markers will be necessary for a more reliable assessment of the species` genetic diversity.

Evidence of a genetic bottleneck in an El Nino affected population of South American fur seals, Arctocephalus australis

The South American fur seal, Arctocephalus australis, was one of the earliest otariid seals to be exploited by humans: at least 6000 years ago on the Atlantic coast and 4000 on the Pacific coast of South America. More than 750,000 fur seals were killed in Uruguay until 1991. However, a climatological phenomenon-the severe 1997-1998 El Nino Southern Oscillation (ENSO)-was responsible for the decline of 72% Of the Peruvian fur seal population due to starvation as a consequence of warming of sea-surface temperatures and primary productivity reduction. Currently, there is no precise information on global population size or on the species` conservation status. The present study includes the first bottleneck test for the Pacific and Atlantic populations of A. australis based on the analysis of seven microsatellite loci. Genetic bottleneck compromises the evolutionary potential of a population to respond to environmental changes. The perspective becomes even more alarming due to current global warming models that predict stronger and more frequent ENSO events in the future. Our analysis found moderate support for deviation from neutrality-equilibrium for the Pacific population of fur seals and none for the Atlantic population. This difference among population reflects different demographic histories, and is consistent with a greater reduction in population size in the Pacific. Such an event could be a result of the synergic effects of recurrent ENSO events and the anthropogenic impact (sealing and prey overfishing) on this population.

East-West cranial differentiation in pre-Columbian populations from Central and North America

In a recent study we found that crania from South Amerindian populations on each side of the Andes differ significantly in terms of craniofacial shape. Western populations formed one morphological group, distributed continuously over 14,000 km from the Fuegian archipelago (southern Chile) to the Zulia region (northwestern Venezuela). Easterners formed another group, distributed from the Atlantic Coast up to the eastern foothills of the Andes. This differentiation is further supported by several genetic studies, and indirectly by ecological and archaeological studies. Some authors suggest that this dual biological pattern is consistent with differential rates of gene flow and genetic drift operating on both sides of the Cordillera due to historical reasons. Here we show that such East-West patterning is also observable in North America. We suggest that the ""ecological zones model"" proposed by Dixon, explaining the spread of the early Americans along a Pacific dispersal corridor, combined with the evolution of different population dynamics in both regions, is the most parsimonious mechanism to explain the observed patterns of within- and between-group craniofacial variability. (c) 2007 Elsevier Ltd. All rights reserved.

Craniometric Similarities Within and Between Human Populations in Comparison with Neutral Genetic Data

The statement that pairs of individuals from different populations are often more genetically similar than pairs from the same population is a widespread idea inside and outside the scientific community. Witherspoon et al. [""Genetic similarities within and between human populations,"" Genetics 176:351-359 (2007)] proposed an index called the dissimilarity fraction (omega) to access in a quantitative way the validity of this statement for genetic systems. Witherspoon demonstrated that, as the number of loci increases, omega decreases to a point where, when enough sampling is available, the statement is false. In this study, we applied the dissimilarity fraction to Howells`s craniometric database to establish whether or not similar results are obtained for cranial morphological traits. Although in genetic studies thousands of loci are available, Howells`s database provides no more than 55 metric traits, making the contribution of each variable important. To cope with this limitation, we developed a routine that takes this effect into consideration when calculating. omega Contrary to what was observed for the genetic data, our results show that cranial morphology asymptotically approaches a mean omega of 0.3 and therefore supports the initial statement-that is, that individuals from the same geographic region do not form clear and discrete clusters-further questioning the idea of the existence of discrete biological clusters in the human species. Finally, by assuming that cranial morphology is under an additive polygenetic model, we can say that the population history signal of human craniometric traits presents the same resolution as a neutral genetic system dependent on no more than 20 loci.

Effects of let-7 microRNA on Cell Growth and Differentiation of Papillary Thyroid Cancer

Papillary thyroid carcinoma (PTC) is the most common endocrine malignancy and RET/PTC rearrangements represent key genetic events frequently associated to this cancer, enhancing proliferation and dedifferentiation by activation of the RET/PTC-RAS-BRAF-mitogen-activated protein kinase (MAPK) pathway. Recently, let-7 microRNA was found to reduce RAS levels in lung cancer, acting as a tumor suppressor gene. Here, we report that RET/PTC3 oncogenic activation in PCCL3 rat thyroid cells markedly reduces let-7f expression. Moreover, stable transfection of let-7 microRNA in TPC-1 cells, which harbor RET/PTC1 rearrangement, inhibits MAPK activation. As a result, let-7f was capable of reducing TPC-1 cell growth, and this might be explained, at least in part, by decreased messenger RNA (mRNA) expression of cell cycle stimulators such as MYC and CCND1 (cyclin D1) and increased P21 cell cycle inhibitor mRNA. In addition, let-7 enhanced transcriptional expression of molecular markers of thyroid differentiation such as TITF1 and TG. Thus, reduced expression of let-7f might be an essential molecular event in RET/PTC malignant transformation. Moreover, let-7f effects on thyroid growth and differentiation might attenuate neoplastic process of RET/PTC papillary thyroid oncogenesis through impairment of MAPK signaling pathway activation. This is the first functional demonstration of an association of let-7 with thyroid cancer cell growth and differentiation.

Population genetic analysis of large sequence polymorphisms in Plasmodium falciparum blood-stage antigens

Plasmodium falciparum, the causative agent of human malaria, invades host erythrocytes using several proteins on the surface of the invasive merozoite, which have been proposed as potential vaccine candidates. Members of the multi-gene PfRh family are surface antigens that have been shown to play a central role in directing merozoites to alternative erythrocyte receptors for invasion. Recently, we identified a large structural polymorphism, a 0.58 Kb deletion, in the C-terminal region of the PfRh2b gene, present at a high frequency in parasite populations from Senegal. We hypothesize that this region is a target of humoral immunity. Here, by analyzing 371 P. falciparum isolates we show that this major allele is present at varying frequencies in different populations within Senegal, Africa, and throughout the world. For allelic dimorphisms in the asexual stage antigens, Msp-2 and EBA-175, we find minimal geographic differentiation among parasite populations from Senegal and other African localities, suggesting extensive gene flow among these populations and/or immune-mediated frequency-dependent balancing selection. In contrast, we observe a higher level of inter-population divergence (as measured by F(st)) for the PfRh2b deletion, similar to that observed for SNPs from the sexual stage Pfs45/48 loci, which is postulated to be under directional selection. We confirm that the region containing the PfRh2b polymorphism is a target of humoral immune responses by demonstrating antibody reactivity of endemic sera. Our analysis of inter-population divergence suggests that in contrast to the large allelic dimorphisms in EBA-175 and Msp-2, the presence or absence of the large PfRh2b deletion may not elicit frequency-dependent immune selection, but may be under positive immune selection, having important implications for the development of these proteins as vaccine candidates. (C) 2009 Elsevier B.V. All rights reserved.

Symmetry in biology: from genetic code to stochastic gene regulation

Mathematical models, as instruments for understanding the workings of nature, are a traditional tool of physics, but they also play an ever increasing role in biology - in the description of fundamental processes as well as that of complex systems. In this review, the authors discuss two examples of the application of group theoretical methods, which constitute the mathematical discipline for a quantitative description of the idea of symmetry, to genetics. The first one appears, in the form of a pseudo-orthogonal (Lorentz like) symmetry, in the stochastic modelling of what may be regarded as the simplest possible example of a genetic network and, hopefully, a building block for more complicated ones: a single self-interacting or externally regulated gene with only two possible states: ` on` and ` off`. The second is the algebraic approach to the evolution of the genetic code, according to which the current code results from a dynamical symmetry breaking process, starting out from an initial state of complete symmetry and ending in the presently observed final state of low symmetry. In both cases, symmetry plays a decisive role: in the first, it is a characteristic feature of the dynamics of the gene switch and its decay to equilibrium, whereas in the second, it provides the guidelines for the evolution of the coding rules.

Genetic diversity of two Brazilian populations of the Pampas deer (Ozotoceros bezoarticus, Linnaeus 1758)

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Evidence of ecotypic differentiation between populations of the tree species Parapiptadenia rigida due to flooding

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

Strain differentiation of Trichophyton rubrum by randomly amplified polymorphic DNA and analysis of rDNA nontranscribed spacer

Trichophyton rubrum is the most common pathogen causing dermatophytosis. Molecular strain-typing methods have recently been developed to tackle epidemiological questions and the problem of relapse following treatment. A total of 67 strains of T rubrum were screened for genetic variation by randomly amplified polymorphic DNA (RAPD) analysis, with two primers, 5'-d[GGTGCGGGAA]-3' and 5'-d[CCCGTCAGCA]-3', as well as by subrepeat element analysis of the nontranscribed spacer of rDNA, using the repetitive subelements TRS-1 and TRS-2. A total of 12 individual patterns were recognized with the first primer and 11 with the second. Phylogenetic analysis of the RAPID products showed a high degree of similarity (> 90 %) among the epidemiologically related clinical isolates, while the other strains possessed 60% similarity. Specific amplification of TRS-1 produced three strain-characteristic banding patterns (PCR types); simple patterns representing one copy of TRS-1 and two copies of TRS-2 accounted for around 85 % of all isolates. It is concluded that molecular analysis has important implications for epidemiological studies, and RAPID analysis is especially suitable for molecular typing in T. rubrum.

Genetic Structure of Populations of the Rice-Infecting Pathogen Rhizoctonia solani AG-1 IA from China

Sheath blight disease (SBD) on rice, caused by Rhizoctonia solani AG-1 IA, is one of the most devastating rice diseases on a global basis, including China (in Eastern Asia), the world's largest rice-growing country. We analyzed the population genetics of nine rice-infecting populations from China using nine microsatellite loci. One allopatric population from India (Southern Asia) was included in the analyses. In total, 300 different multilocus genotypes were found among 572 fungal isolates. Clonal fractions within rice fields were 16 to 95%, suggesting that sclerotia were a major source of primary inoculum in some fields. Global Phi(ST) statistics (Phi(ST) = 42.49; P <= 0.001) were consistent with a relatively high level of differentiation among populations overall; however, pairwise comparisons gave nonsignificant R(ST) values, consistent with contemporary gene flow among five of the populations. Four of these populations were located along the Yangtze River tributary network. Gene flow followed an isolation-by-distance model consistent with restricted long-distance migration. Historical migration rates were reconstructed and yielded values that explained the current levels of population subdivision. Except for one population which appeared to be strictly clonal, all populations showed evidence of a mixed reproductive mode, including both asexual and sexual reproduction. One population had a strictly recombining structure (all loci were in Hardy-Weinberg equilibrium) but the remaining populations from China and the one from India exhibited varying degrees of sexual reproduction. Six populations showed significant F(IS) values consistent with inbreeding.

Allozyme differentiation of two populations of the genus Neoplecostomus Eigenmann & Eigenmann, 1888 (Teleostei, Loricariidae) from the upper Paraná River basin, Brazil

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