Search

Climate changes and associated glacial and interglacial cycles during the Quaternary strongly influenced the evolutionary history of countless number of species. Subtropical highland grasslands (SHG) in southern South America constitute a distinct vegetation type with high diversity and endemism rates. The most recent common ancestor of the genus Petunia (Solanaceae) originated in lowland grasslands, and some lineages secondarily migrated to SHG. Here, we describe the evolutionary history of an SHG species, P. altiplana, distributed throughout a wide area in a river-fragmented landscape. We used plastid and nuclear markers to evaluate the role of the Pelotas River and Quaternary climate cycles over the genetic structure and historical demography of this species based on a phylogeographical approach. We found moderate population expansions during the last 25 kyr, with a more recent (c. 1.6 kya) divergence between two groups of populations from opposite river margins, possibly caused by the expansion of the Araucaria Forest along the river valley, effectively isolating distinct grassland patches.

A headwater or river capture is a phenomenon commonly invoked to explain the absence of reciprocal monophyly of genetic lineages among isolated hydrographic basins in freshwater fish. Under the assumption of river capture, a secondary contact between populations previously isolated in different basins explains the observed genetic pattern. However, the absence of reciprocal monophyly could also arise under population isolation through the retention of ancestral of polymorphisms. Here, we applied an approximate Bayesian computation (ABC) framework for estimating the relative probability of scenarios with and without secondary contact. We used Cnesterodon decemmaculatus as a study model because of the multiple possible cases of river capture and the demographic parameters estimated in a previous mitochondrial DNA study that are useful for simulating scenarios to test both hypotheses using the ABC framework. Our results showed that, in general, mitochondrial DNA is useful for distinguishing between these alternative demographic scenarios with reasonable confidence, but in extreme cases (e.g. recent divergence or large population size) there is no power to discriminate between scenarios. Testing hypotheses of drainage rearrangement under a statistically rigorous framework is fundamental for understanding the evolution of freshwater fish fauna as a complement to, or in the absence of, geological evidence.

We report the occurrence of an invasive alien species, palometa Serrasalmus maculatus, in the Patos Lagoon drainage. Primary occurrence data were based on three specimens captured and preserved as vouchers in scientific collections. Additionally, we searched for secondary records from unpublished scientific sources, public agencies reports and media news to find additional reports. We discussed the possible pathways of invasion, suggesting as the vector of introduction transpositions from the Uruguay River basin. Ecological implications for ichthyofauna, environmental impacts and risk of other events of invasion in the adjoining basins are discussed.

The Atlantic spotted dolphin (Stenella frontalis) is endemic to tropical, subtropical, and warm temperate waters of the Atlantic Ocean. Throughout its distribution, both geographic distance and environmental variation may contribute to population structure of the species. In this study, we follow a seascape genetics approach to investigate population differentiation of Atlantic spotted dolphins based on a large worldwide dataset and the relationship with marine environmental variables. The results revealed that the Atlantic spotted dolphin exhibits population genetic structure across its distribution based on mitochondrial DNA control region (mtDNA-CR) data. Analyses based on the contemporary landscape suggested, at both the individual and population level, that the population genetic structure is consistent with the isolation-by-distance model. However, because geography and environmental matrices were correlated, and because in some, but not all analyses, we found a significant effect for the environment, we cannot rule out the addition contribution of environmental factors in structuring genetic variation. Future analyses based on nuclear data are needed to evaluate whether local processes, such as social structure and some level of philopatry within populations, may be contributing to the associations among genetic structure, geographic, and environmental distance.

Brazil is the fifth largest country in the world in area and the fifth most populous. The Brazilian voluntary Bone Marrow Donor Registry is the third largest in terms of number of donors in the world, being a valuable source of HLA genetics to characterize the donor population of Brazil as well. The genetic background of the Brazilian population is quite heterogeneous, resulting from 5 centuries of admixture among Native Americans, Europeans and Africans, making the Brazilian population unique in terms of genetic ancestry. The unique characteristics of populations in different Brazilian regions make them an exciting focus for genetic diversity studies. Studies on HLA genetic diversity of Brazilian populations have been conducted since the late 1980s and, in this review, we highlight the main findings from studies carried out in Brazil based on classical HLA. In addition, we calculated the genetic distance from the molecular data of the studies included in this review in order to have a broader view of the HLA diversity in Brazilian populations. We emphasize that characterization of HLA diversity is not only important for transplantation programs, but can shed a light on ancestry, history and other demographic patterns with or without association with autoimmune disease.

Objectives The aim of this study was to create a comprehensive summary of available mtDNA and Y-chromosome data for Native Americans from North, Central, and South America, including both modern and ancient DNA. To illustrate the usefulness of this dataset we present a broad picture of the genetic variation for both markers across the Americas. Methods We searched PubMed, ResearchGate, Google Scholar for studies about mtDNA or Y-chromosome variation in Native American populations, including geographic, linguistic, ecological (ecoregion), archeological and chronological information. We used AMOVA to estimate the genetic structure associated with language and ecoregion grouping and Mantel tests to evaluate the correlation between genetic and geographic distances. Results Genetic data were obtained from 321 primary sources, including 22,569 individuals from 298 contemporary populations, and 3628 individuals from 202 archeological populations. MtDNA lineages of probable non-Amerindian origin were rare, in contrast with Y-chromosome lineages. Mantel tests showed a statistically significant correlation for the whole continent considering mtDNA but not the Y-chromosome. Genetic structure between groups was always stronger for mtDNA than for the Y-chromosome. Conclusions This study summarizes decades of research conducted in Native American populations for both mtDNA and the Y-chromosome. Continental or sub-continental patterns of variation reveal that most of the genetic variation occurs within populations rather than among linguistic or ecoregional groups, and that isolation by distance is barely detectable in most population sets. The genetic structure among groups was always larger for mtDNA than for the Y-chromosome, suggesting between-sex differences in gene flow.

Since the 1960s, East African athletes, mainly from Kenya and Ethiopia, have dominated long-distance running events in both the male and female categories. Further demographic studies have shown that two ethnic groups are overrepresented among elite endurance runners in each of these countries: the Kalenjin, from Kenya, and the Oromo, from Ethiopia, raising the possibility that this dominance results from genetic or/and cultural factors. However, looking at the life history of these athletes or at loci previously associated with endurance athletic performance, no compelling explanation has emerged. Here, we used a population approach to identify peaks of genetic differentiation for these two ethnicities and compared the list of genes close to these regions with a list, manually curated by us, of genes that have been associated with traits possibly relevant to endurance running in GWAS studies, and found a significant enrichment in both populations (Kalenjin, P = 0.048, and Oromo, P = 1.6x10-5). Those traits are mainly related to anthropometry, circulatory and respiratory systems, energy metabolism, and calcium homeostasis. Our results reinforce the notion that endurance running is a systemic activity with a complex genetic architecture, and indicate new candidate genes for future studies. Finally, we argue that a deterministic relationship between genetics and sports must be avoided, as it is both scientifically incorrect and prone to reinforcing population (racial) stereotyping.

Illegal capture and trade of wild birds has long been a threat to biodiversity. Translocation—the release of individuals from one location into another—is a useful conservation tool in the management of species. However, both health (such as different pathogens) and adaptive (such as local adaptation), differences among populations must be taken into account, as both can impact the recipient population negatively. Here, we provide health and genetic information to support release planning for two of the most trafficked Brazilian wild bird species (Paroaria coronataandSaltator similis). We focused on two fundamental questions: Are there significant differences in pathogen load between wild and captive populations? Is there significant genetic structure among populations? In total, 223 free-living birds were captured, sampled, and released at the same site. Devices and live decoys characteristics were top factors influencing captures. We tested blood, feces, and oropharyngeal swabs from free-ranging (n=101) and confiscated (n=92) birds for Newcastle disease virus,Salmonellaspp., andMycoplasma gallisepticum. Genetic structure among populations was investigated using mtDNA in a subsample of these birds. We found no evidence for Newcastle disease virus andSalmonellaspp. in seized and free-living birds from both species. However, seizedP. coronataandS. similismay be potential sources ofM. gallisepticum. We found significant but low genetic structure among populations occurring in different Biomes (ΦCT=0.26 forP. coronata;ΦCT=0.13 forS. similis) and no significant structure among populations occurring in the Pampa Biome. These results suggest that while it may be important to screen seized birds for avian pathogens, genetic structure among populations seems to be of lesser concern when planning the release of seized songbirds in the wild.

Domestic turkeys present several color phenotypes controlled by at least five genetic loci, but only one of these has been identified precisely: the bronze locus, which turned out to be the melanocortin-1 receptor (MC1R) gene. MC1R variation is important for breeders interested in maintaining or developing different color varieties. In this study, we sequenced most of the MC1R gene from 16 White Holland (the main commercial turkey variety) and 19 pigmented turkeys from southern Brazil with two purposes. The first was to describe the MC1R diversity in White Holland turkeys, which may serve as reservoirs of genetic diversity at this locus. The second was to test whether the traditional color classification used by Brazilian breeders is related to previously known MC1R alleles. White Holland turkeys had four different haplotypes corresponding to the bronze (b +) and black-winged bronze (b 1) alleles. Pigmented turkeys also had four haplotypes corresponding to the b + and b 1 alleles, but different haplotypes represent the most common b + allele in these two groups. The black (B) allele was absent from our samples. Overall, our results suggest that white and pigmented individuals form two different populations, and that the traditional color classification used by Brazilian breeders cannot accurately predict the genotypes at the bronze locus.

Human leukocyte antigen (HLA) genes are very informative in population genetics studies and their variability has been widely used to reconstruct the history of geographic and/or demographic expansions of human populations. The characterization of HLA diversity at the population level is also fundamental in clinical studies, particularly for bone marrow transplantation programs. In this study, we investigated the HLA molecular variation in Rio Grande do Sul, South Brazil, in order to identify possible regional differences across this state. More than 97,000 bone marrow donors were typed at the HLA- A, -B and -DRB1 loci and analyzed by considering two kinds of subdivisions based on both self-identified ethnicity and place of residence: (a) the official geographic subdivision defined by the Brazilian Institute of Geography and Statistics and (b) known information about the colonization history of the state. HLA allele and haplotype frequencies were estimated and compared among the defined subgroups. The results indicate a lack of correlation between genetic variation and geography and thus no clear HLA genetic structure based on geographic criteria. On the other hand, major differences were observed regarding ethnicity. In addition, local populations from Rio Grande do Sul were found to be genetically similar to their corresponding parental European populations from Germany, Italy and Portugal, as documented by historical data. Overall, this study provides a thorough characterization of the HLA genetic variation in Rio Grande do Sul and a better understanding of its demographic history, being most useful for the development of more efficient strategies in bone marrow donors' recruitment.

Characterizing molecular evolution patterns of the Hepatitis B Virus (HBV) is important for a better understanding of the natural history of this infection. However, several molecular evolution estimates are conditioned on tree topology. There is no consensus about the phylogenetic relationships of HBV genotypes, and different studies often find alternative topologies. While most studies consider HBV genotypes F and H as sister to all other human genotypes, a recent study suggested an alternative HBV phylogeny that indicates an accelerated substitution rate for HBV-F/H partially driven by positive selection. In this study, we evaluate the impact of alternative HBV topologies on inferences of HBV phylogeny, rate acceleration, and positive selection on the HBV-F/H branch. Our results indicate that under certain methodological approaches alternative HBV topologies are equally likely. Considering phylogenetic uncertainty, there is no evidence that HBV-F/H had an accelerated substitution rate, even though inferences of positive selection are robust to alternative background topologies. Our results further suggest that, under reasonable assumptions, HBV-F/H most likely represents the sister lineage to all other human/ape HBV genotypes. Understanding the full range of likely topologies will be crucial for elaborating, testing, and refining hypothesis about the evolutionary HBV origins in our species.

Objectives The major aim of this article was to estimate the demographic impact of European arrival and colonization over Native American populations from southern Brazil and Uruguay. We also compared the mitochondrial DNA (mtDNA) genetic diversity, structure, and demography of Native American lineages present in current indigenous (Natives) and nonindigenous admixed (Admixed) populations to estimate the effective population size (Ne ) of contemporary and ancestral (pre-Columbian) Native American populations. Methods We retrieved published mtDNA sequences from Native (n = 396) and Admixed (n = 309) populations from southern Brazil, Uruguay, and surrounding areas. We conducted genetic diversity, structure, and demographic analyses. Finally, we used Approximate Bayesian Computation to estimate the Ne for current Native, Admixed, and pre-Columbian Native American populations. Results We found higher Native American mtDNA genetic diversity in admixed rather than in indigenous populations (131/309 vs 27/396 different haplotypes, respectively). Only Admixed populations maintained ancient signals of the Native American population expansion approximately 14 to 17 kya, which have decayed in Natives. Our Ne estimates suggest that Natives represent only 0.33% (0.18%-1.19%) of the Ne for ancestral pre-Columbian indigenous populations. Conclusions Admixed populations represent an important genetic reservoir of Native American lineages, many of which are extinct in contemporary indigenous populations. In addition, the Native American lineages present in Admixed populations retain part of the past demographic history of Native Americans. The intensity of the reduction is congruent with historical accounts of strong indigenous depopulation during the colonization process.

Background HLA genes are the most polymorphic of the human genome and have distinct allelic frequencies in populations of different geographical regions of the world, serving as genetic markers in ancestry studies. In addition, specific HLA alleles may be associated with various autoimmune and infectious diseases. The bone marrow donor registry in Brazil is the third largest in the world, and it counts with genetic typing of HLA-A, -B, and -DRB1. Since 1991 Brazil has maintained the DATASUS database, a system fed with epidemiological and health data from compulsory registration throughout the country. Methods In this work, we perform spatial analysis and georeferencing of HLA genetic data from more than 86,000 bone marrow donors from Rio Grande do Sul (RS) and data of hospitalization for rheumatoid arthritis, multiple sclerosis and Crohn’s disease in RS, comprising the period from 1995 to 2016 obtained through the DATASUS system. The allele frequencies were georeferenced using Empirical Bayesian Kriging; the diseases prevalence were georeferenced using Inverse Distance Weighted and cluster analysis for both allele and disease were performed using Getis-Ord Gi* method. Spearman’s test was used to test the correlation between each allele and disease. Results The results indicate a HLA genetic structure compatible with the history of RS colonization, where it is possible to observe differentiation between regions that underwent different colonization processes. Spatial analyzes of autoimmune disease hospitalization data were performed revealing clusters for different regions of the state for each disease analyzed. The correlation test between allelic frequency and the occurrence of autoimmune diseases indicated a significant correlation between the HLA-B*08 allele and rheumatoid arthritis. Conclusions Genetic mapping of populations and the spatial analyzes such as those performed in this work have great economic relevance and can be very useful in the formulation of public health campaigns and policies, contributing to the planning and adjustment of clinical actions, as well as informing and educating professionals and the population. Electronic supplementary material The online version of this article (10.1186/s12942-018-0154-8) contains supplementary material, which is available to authorized users.

Elasmobranchs, the group of cartilaginous fishes that include sharks and rays, are especially vulnerable to overfishing due to low fecundity and late sexual maturation. A significant number of elasmobranch species are currently overexploited or threatened by fisheries activities. Additionally, several recent reports have indicated that there has been a reduction in regional elasmobranch population sizes. Brazil is an important player in elasmobranch fisheries and one of the largest importers of shark meat. However, carcasses entering the shark meat market have usually had their fins and head removed, which poses a challenge to reliable species identification based on the morphology of captured individuals. This is further complicated by the fact that the internal Brazilian market trades several different elasmobranch species under a common popular name: "cacao." The use of such imprecise nomenclature, even among governmental agencies, is problematic for both controlling the negative effects of shark consumption and informing the consumer about the origins of the product. In this study, we used DNA barcoding (mtDNA, COI gene) to identify, at the species level, "cacao" samples available in local markets from Southern Brazil. We collected 63 samples traded as "cacao," which we found to correspond to 20 different species. These included two teleost species: Xiphias gladius (n = 1) and Genidens barbus (n = 6), and 18 species from seven elasmobranch orders (Carcharhiniformes, n = 42; Squaliformes, n = 3; Squatiniformes, n = 2; Rhinopristiformes, n = 4; Myliobatiformes, n = 3; Rajiformes, n = 1; and Torpediniformes, n = 1). The most common species in our sample were Prionace glauca (n = 15) and Sphyrna lewini (n = 14), while all other species were represented by four samples or less. Considering IUCN criteria, 47% of the elasmobranch species found are threatened at the global level, while 53% are threatened and 47% are critically endangered in Brazil. These results underline that labeling the meat of any shark species as "cacao" is problematic for monitoring catch allocations from the fishing industry and discourages consumer engagement in conservationist practices through informed decision-making.

In spite of many genetic studies that contributed for a deep knowledge about the peopling of the Americas, no consensus has emerged about important parameters such as the effective size of the Native Americans founder population. Previous estimates based on genomic datasets may have been biased by the use of admixed individuals from Latino populations, while other recent studies using samples from Native American individuals relied on approximated analytical approaches. In this study we use resequencing data for nine independent regions in a set of Native American and Siberian individuals and a full-likelihood approach based on isolation-with-migration scenarios accounting for recent flow between Asian and Native American populations. Our results suggest that, in agreement with previous studies, the effective size of the Native American population was small, most likely in the order of a few hundred individuals, with point estimates close to 250 individuals, even though credible intervals include a number as large as ~4,000 individuals. Recognizing the size of the genetic bottleneck during the peopling of the Americas is important for determining the extent of genetic markers needed to characterize Native American populations in genome-wide studies and to evaluate the adaptive potential of genetic variants in this population.

The Pampas is a Neotropical biome formed primarily by low altitude grasslands and encompasses the southernmost portion of Brazil, Uruguay, and part of Argentina. Despite the high level of endemism, and its significant environmental heterogeneity, Pampean species are underrepresented in phylogeographic studies, especially aquatic organisms. The Pampean hydrological system resulted from a long history of tectonism, climate, and sea level changes since the Neogene. In this study, we examined the population genetic structure of Cnesterodon decemmaculatus, a freshwater fish species that occurs throughout most of the Pampa biome. We characterized mitochondrial and autosomal genetic lineages in populations sampled from Southern Brazil and Uruguay to investigate (1) the correspondence between current drainage systems and evolutionary lineages, (2) the demographic history for each genetic lineage, and (3) the temporal depth of these lineages. Overall, we found that the major evolutionary lineages in this species are strongly related to the main Pampean drainage systems, even though stream capture events may have affected the distribution of genetic lineages among drainages. There was evidence for recent population growth in the lineages occupying drainages closest to the shore, which may indicate the effect of quaternary sea-level changes. In general, divergence time estimates among evolutionary lineages were shallow, ranging from 20,000 to 800,000 years before present, indicating a geologically recent history for this group, as previously reported in other Pampean species. A Bayesian phylogeographical reconstruction suggested that an ancestral lineage probably colonized the Uruguay River Basin, and then expanded throughout the Pampas. This evolutionary scenario may represent useful starting models for other freshwater species having a similar distribution.

The use of traps is extremely important in several types of ecological studies, and may assist in the capture of individuals in areas that are difficult to access. In the present study, we compared the effectiveness of wooden (Schramm) versus "Tomahawk" traps to capture Salvator merianae (Duméril & Bibron, 1839) lizards. The study was conducted in Eldorado do Sul, Rio Grande do Sul, Brazil. Field data were collected from August 2013 to March 2015, during the reproductive period of the species. The study involved two types of baited traps: i) "Tomahawk", made of galvanized steel; and ii) Schramm, a wooden trap. The capture rate of the Schramm wooden traps was 1.63 individuals/day, and of the "Tomahawk" was 0.36 individuals/day. These results are important for researchers working with large lizards and may help to increase sampling efficiency for these organisms.

Freshwater fishes often display a marked phylogeographic structure strongly associated with historical and ecological changes in the aquatic environment. Different ecological conditions in the same river drainage may act as permeable barriers to dispersion and gene flow. Previous studies recognized two discrete spatial components for the ichthyofauna in the freshwater coastal drainages of southern Brazil: the lowland fish fauna in the lagoons and the fish fauna of the rivers flowing in the valleys. In order to test if the coastal lagoons may limit the dispersion of a riverine species, we describe the phylogeographic structure among populations of Cyanocharax itaimbe, a species endemic to this region. We analysed 55 specimens characterized for two mitochondrial and one nuclear genes. Sequences were analysed using gene trees and species tree approaches, together with standard population genetics methods. Molecular analyses indicated three evolutionary groups which diverged from each other between an estimated 1,600,000 and 450,000 years before the present. However, two currently isolated river systems share the same evolutionary clade, whereas a single drainage contains two different lineages. Our results indicate strong genetic structure among populations along with generally conserved morphology. The strong genetic structure among populations living in the same drainage system may be explained by ecological differences between lagoons and rivers (or palaeochannels) that act as barriers to dispersion.

The FOXP subfamily is probably the most extensively characterized subfamily of the forkhead superfamily, playing important roles in development and homeostasis in vertebrates. Intrinsically disorder protein regions (IDRs) are protein segments that exhibit multiple physical interactions and play critical roles in various biological processes, including regulation and signaling. IDRs in proteins may play an important role in the evolvability of genetic systems. In this study, we analyzed 77 orthologous FOXP genes/proteins from Tetrapoda, regarding protein disorder content and evolutionary rate. We also predicted the number and type of short linear motifs (SLIMs) in the IDRs. Similar levels of protein disorder (approximately 70%) were found for FOXP1, FOXP2, and FOXP4. However, for FOXP3, which is shorter in length and has a more specific function, the disordered content was lower (30%). Mammals showed higher protein disorders for FOXP1 and FOXP4 than non-mammals. Specific analyses related to linear motifs in the four genes showed also a clear differentiation between FOXPs in mammals and non-mammals. We predicted for the first time the role of IDRs and SLIMs in the FOXP gene family associated with possible adaptive novelties within Tetrapoda. For instance, we found gain and loss of important phosphorylation sites in the Homo sapiens FOXP2 IDR regions, with possible implication for the evolution of human speech.

Bryconamericus lethostigmus is the type-species of the monotypic genus Odontostoechus, diagnosed in part based on the presence of a unique tooth series in the premaxilla. Recently a new proposal of classification of the Stevardiinae placed Odontostoechus as a junior synonym of a monophyletic genus Bryconamericus sensu stricto, a genus characterized by the presence of two tooth series. Bryconamericus lethostigmus is redescribed herein and the single tooth series in the premaxilla is demonstrated to originate from merging of the external tooth row with the inner row during ontogeny refuting primary hypothesis of homology between the mouth morphology of B. lethostigmus and the genera Bryconacidnus, Ceratobranchia, Monotocheirodon, Othonocheirodus, Rhinopetitia and Rhinobrycon. A phylogeographic analysis indicated that the pattern described for the sympatric species Diapoma itaimbe is not mirrored by B. lethostigmus. The results also do not support the hypothesis of a new species in the rio Araranguá drainage. RESUMO Bryconamericus lethostigmus é a espécie tipo do gênero monotípico Odontostoechus, diagnosticado em parte pela presença de uma única série de dentes na pré-maxila. Recentemente uma nova proposta de classificação de Stevardiinae considerou Odontostoechus como sinônimo júnior do gênero monofilético Bryconamericus sensu stricto caracterizado pela presença de duas séries de dentes. Bryconamericus lethostigmus é redescrito e demonstra-se que a única série de dentes no premaxilar se origina pela junção da série externa de dentes com a série interna durante a ontogenia, refutando a hipótese de homologia primária entre a morfologia da boca de B. lethostigmus e os gêneros Bryconacidnus, Ceratobranchia, Monotocheirodon, Othonocheirodus, Rhinopetitia e Rhinobrycon. A análise filogeográfica indica que o padrão descrito para a espécie simpátrica Diapoma itaimbe não se repete em B. lethostigmus. Os resultados também não suportam a hipótese de uma nova espécie para a bacia do rio Araranguá.

Understanding the spatiotemporal distribution of genetic variation and the ways in which this distribution is connected to the ecological context of natural populations is fundamental for understanding the nature and mode of intraspecific and, ultimately,interspecific differentiation. The Petunia axillaris complex is endemic to the grasslands of southern South America and includes three subspecies: P. a. axillaris, P. a. parodii and P. a. subandina. These subspecies are traditionally delimited based on both geography and floral morphology, although the latter is highly variable. Here, we determined the patterns of genetic (nuclear and cpDNA), morphological and ecological (bioclimatic) variation of a large number of P. axillaris populations and found that they are mostly coincident with subspecies delimitation. The nuclear data suggest that the subspecies are likely independent evolutionary units, and their morphological differences may be associated with local adaptations to diverse climatic and/or edaphic conditions and population isolation. The demographic dynamics over time estimated by skyline plot analyses showed different patterns for each subspecies in the last 100 000 years, which is compatible with a divergence time between 35 000 and 107 000 years ago between P. a. axillaris and P. a. parodii, as estimated with the IMa program. Coalescent simulation tests using Approximate Bayesian Computation do not support previous suggestions of extensive gene flow between P. a. axillaris and P. a. parodii in their contact zone. Fil: Turchetto, Caroline. Universidade Federal Do Rio Grande Do Sul; Brasil; Fil: Fagundes, Nelson J. R.. Universidade Federal Do Rio Grande Do Sul; Brasil; Fil: Segatto, Ana L. A.. Universidade Federal Do Rio Grande Do Sul; Brasil; Fil: Kuhlemeier, Cris. Institute of Plant Science; Suiza; Fil: Solis Neffa, Viviana Griselda. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Nordeste. Instituto de Botánica del Nordeste (i); Argentina Fil: Speranza, Pablo R. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Nordeste. Instituto de Botánica del Nordeste (i); Argentina Fil: Bonatto, Sandro L.. Pontificia Universidade Catolica Do Rio Grande Do Sul; Brasil; Fil: Freitas, Loreta B.. Universidade Federal Do Rio Grande Do Sul; Brasil

In the last decades researchers have been able to determine the molecular basis of some phenotypes, to test for evidence of natural selection upon them, and to demonstrate that the same genes or genetic pathways can be associated with convergent traits. Colour traits are often subject to natural selection because even small changes in these traits can have a large effect on fitness via camouflage, sexual selection or other mechanisms. The melanocortin-1 receptor locus (MC1R) is frequently associated with intraspecific coat colour variation in vertebrates, but it has been far harder to demonstrate that this locus is involved in adaptive interspecific colour differences. Here, we investigate the contribution of the MC1R gene to the colour diversity found in toucans (Ramphastidae). We found divergent selection on MC1R in the clade represented by the genus Ramphastos and that this coincided with the evolution of darker plumage in members of this genus. Using phylogenetically-corrected correlations we show significant and specific relationships between the rate of non-synonymous change in MC1R (dN) and plumage darkness across Ramphastidae, and also between the rate of functionally significant amino acid changes in MC1R and plumage darkness. Furthermore, three of the seven amino acid changes in MC1R that occurred in the ancestral Ramphastos branch are associated with melanism in other birds. Taken together, our results suggest that the dark colour of Ramphastos toucans was related to non-synonymous substitutions in MC1R that may have been subject to positive selection or to a relaxation of selective pressure. These results also demonstrate a quantitative relationship between gene and phenotype evolution, representing an example of how MC1R molecular evolution may affect macroevolution of plumage phenotypes. This article is protected by copyright. All rights reserved.

Native American populations present the highest prevalence of Hepatitis B Virus (HBV) infection in the Americas, which may be associated to severe disease outcomes. Ten HBV genotypes (A–J) have been described, displaying a remarkable geographic structure, which most likely reflects historic patterns of human migrations. In this study, we characterize the HBV strains circulating in a historical sample of Native South Americans to characterize the historical viral dynamics in this population. The sample consisted of 1070 individuals belonging to 38 populations collected between 1965 and 1997. Presence of HBV DNA was checked by quantitative real-time PCR, and determination of HBV genotypes and subgenotypes was performed through sequencing and phylogenetic analysis of a fragment including part of HBsAg and Pol coding regions (S/Pol). A Bayesian Skyline Plot analysis was performed to compare the viral population dynamics of HBV/A1 strains found in Native Americans and in the general Brazilian population. A total of 109 individuals were positive for HBV DNA (~ 10%), and 70 samples were successfully sequenced and genotyped. Subgenotype A1 (HBV/A1), related to African populations and the African slave trade, was the most prevalent (66–94%). The Skyline Plot analysis showed a marked population expansion of HBV/A1 in Native Americans occurring more recently (1945–1965) than in the general Brazilian population. Our results suggest that historic processes that contributed to formation of HBV/A1 circulating in Native American are related with more recent migratory waves towards the Amazon basin, which generated a different viral dynamics in this region.

The discovery and characterization of informative intraspecific genetic markers is fundamental for evolutionary and conservation genetics studies. Here, we used nuclear ribosomal ITS sequences to access intraspecific genetic diversity in 23 species of the genus Passiflora L. Some degree of variation was detected in 21 of these. The Passiflora and Decaloba (DC.) Rchb. subgenera showed significant differences in the sizes of the two ITS regions and in GC content, which can be related to reproductive characteristics of species in these subgenera. Furthermore, clear geographical patterns in the spatial distribution of sequence types were identified in six species. The results indicate that ITS may be a useful tool for the evaluation of intraspecific genetic variation in Passiflora.

To assess the extent of genetic variability of rDNA intergenic spacer (IGS) in Metarhizium sp., 34 strains (27 isolated in Brazil) were sequenced and analyzed together with an additional 20 Metarhizium anisopliae var. anisopliae sequences retrieved from GenBank. Overall, the global nucleotide diversity for the region under study was of 0.090, while for the Brazilian isolates it was only 0.016. Phylogenetic analyses showed four well-supported groups (A, B, C, and D), one of which (D) has not been previously identified. All but one of the Brazilian strains cluster in this novel D phylogroup, suggesting that the genetic variation found in Brazil is a subset of the worldwide M. anisopiliae var. anisopliae variation.

In the Southwestern Atlantic Ocean, humpback whales migrate every winter to the Brazilian coast for breeding and calving in the Abrolhos Bank. This breeding stock represents the remnants of a larger population heavily exploited during the beginning of the 20th century. Despite its relevance to conservation efforts, the degree of current genetic variation and the migratory relationship with Antarctic feeding areas for this population are still largely unknown. To examine these questions, we sequenced ∼400 bp of the mitochondrial DNA control region from samples taken off the Brazilian coast (n = 171) and near the Antarctic Peninsula (n = 77). The genetic variability of the Brazilian humpback whale breeding population was high and similar to that found in other Southern Hemisphere breeding grounds. Phylogenetic analysis suggested the existence of a new mitochondrial clade that exists at low frequency among Southern Hemisphere populations. Direct comparison between the Brazilian and the Colombia breeding populations and the Antarctic Peninsula feeding population showed no genetic differentiation between this feeding region and the Colombian breeding area or between feeding Areas I and II near the Antarctic Peninsula. In contrast, these populations were genetically distinct from the Brazilian population. Two humpback whales sampled off South Georgia Islands, in the Scotia Sea, shared identical haplotypes to whales from Brazil. Our results, supported by photo-identification and satellite telemetry data, suggest that the main feeding area of the Southern Hemisphere humpback whale population is likely to be located near the South Georgia/South Sandwich Islands area and not in the Antarctic Peninsula.

The Pampas is a biologically rich South American biome, but is poorly represented in phylogeographic studies. While the Pleistocene glacial cycles may have affected the evolutionary history of species distributed in forested biomes, little is known about their effects on the habitats that remained stable through glacial cycles. The South American Pampas have been covered by grasslands during both glacial and interglacial periods and therefore represent an interesting system to test whether the genetic structure in such environments is less pronounced. In this study, we sampled Pampean populations of Homonota uruguayensis from Southern Brazil and Uruguay to assess the tempo and mode of population divergence, using both morphological measurements and molecular markers. Our results indicate that, in spite of its narrow geographic distribution, populations of H. uruguayensis show high levels of genetic structure. We found four major well-supported mtDNA clades with strong geographic associations. Estimates of their divergence times fell between 3.16 and 1.82 million years before the present. Populations from the central portion of the species distribution, on the border between Uruguay and Brazil, have high genetic diversity and may have undergone a population expansion approximately 250,000 years before the present. The high degree of genetic structure is reflected in the analyses of morphological characters, and most individuals could be correctly assigned to their parental population based on morphology alone. Finally, we discuss the biogeographic and conservation implications of these findings.

Summary The low density lipoprotein receptor gene (LDLR) contains many Alu insertions, and is especially Alu-rich at its 3 � -untranslated region (3 � -UTR). Previous studies suggested that the LDLR 3 � -UTR could regulate gene expression by the stabilization of its mRNA. Given the faster Alu evolutionary rate, and wondering about its consequences in a possibly regulatory locus, we have studied ∼ 800 bp of 222 chromosomes from individuals of African, Asian, Caucasian and Amerind ancestry, to better understand the evolution of the worldwide genetic diversity at this locus. Twenty-one polymorphic sites, distributed in 15 haplotypes, were found. High genetic diversity was observed, concentrated in one Alu insertion (Alu U), which also shows a fast evolutionary rate. Genetic diversity is similar in all populations except Amerinds, suggesting a bottleneck during the peopling of the American continent. Three haplotype clusters (A, B, C) are distinguished, cluster A being the most recently formed (∼ 500,000 years ago). No clear geographic structure emerges from the haplotype network, the global Fst (0.079) being lower than the average for the human genome. When ancestral population growth is taken into account, neutrality statistics are higher than expected, possibly suggesting the action of balancing selection worldwide.

A total of 91 Ayoreo individuals previously studied for blood groups and protein markers, living in two Bolivian and one Paraguayan communities, were extensively investigated in relation to the hypervariable 1 segment of the control and eight coding regions of their mitochondrial DNA (mtDNA). They show an extremely reduced mtDNA variability, the observed haplotypes being classifiable in just two haplogroups (C and D). They were also variously studied in relation to six Alu insertions; in this case, however, the prevalences found did not depart markedly from those obtained in other populations of this ethnic group. To assess the Ayoreo position in relation to these populations, 11 other groups that had also been studied for these systems and for blood groups and proteins were selected. The dendrograms obtained with two of the three sets of markers showed distinct patterns, but the Ayoreo were placed in a central position in both the blood group + proteins and Alu insertions trees. Therefore, they are clearly distinct in relation to their mtDNA only, suggesting a strong founder effect and/or random loss of variability in this system. Am. J. Hum. Biol. 16:479–488, 2004.© 2004 Wiley-Liss, Inc.

The Brazil Ministry of Health maintains a Registry of Bone Marrow Donors that corresponds to approximately 12% of the Bone Marrow Donors Worldwide registry. This registry contains information on ethnicity (by self-assessment of color) and HLA-A, -B, and -DRB1 type. The self-assessment of color tool has been extensively used for admixed population characterization. In this context, Brazil represents a highly admixed population, resulting from 5 centuries of colonization and interbreeding, mainly, but not exclusively, among Native Americans, Europeans, and Africans. Here we evaluated self-assessed skin color and HLA genetic information from 71,291 bone marrow donors of southern Brazil to verify how likely is the HLA profiling correspondence within and between self-assessed color groups. We found that HLA itself was a better ancestry indicator than was self-assessed color. Therefore, self-assessment of color in highly admixed populations, such as that of Brazil, is not indicative of higher correspondence in the HLA profiles within skin color groups.

Background Archaeology reports millenary cultural contacts between Peruvian Coast-Andes and the Amazon Yunga, a rainforest transitional region between Andes and Lower Amazonia. To clarify the relationships between cultural and biological evolution of these populations, in particular between Amazon Yungas and Andeans, we used DNA-sequence data, a model-based Bayesian approach and several statistical validations to infer a set of demographic parameters. Results We found that the genetic diversity of the Shimaa (an Amazon Yunga population) is a subset of that of Quechuas from Central-Andes. Using the Isolation-with-Migration population genetics model, we inferred that the Shimaa ancestors were a small subgroup that split less than 5300 years ago (after the development of complex societies) from an ancestral Andean population. After the split, the most plausible scenario compatible with our results is that the ancestors of Shimaas moved toward the Peruvian Amazon Yunga and incorporated the culture and language of some of their neighbors, but not a substantial amount of their genes. We validated our results using Approximate Bayesian Computations, posterior predictive tests and the analysis of pseudo-observed datasets. Conclusions We presented a case study in which model-based Bayesian approaches, combined with necessary statistical validations, shed light into the prehistoric demographic relationship between Andeans and a population from the Amazon Yunga. Our results offer a testable model for the peopling of this large transitional environmental region between the Andes and the Lower Amazonia. However, studies on larger samples and involving more populations of these regions are necessary to confirm if the predominant Andean biological origin of the Shimaas is the rule, and not the exception. Electronic supplementary material The online version of this article (doi:10.1186/s12862-014-0174-3) contains supplementary material, which is available to authorized users.

The Magnificent Frigatebird Fregata magnificens has a pantropical distribution, nesting on islands along the Atlantic and Pacific coasts. In the Caribbean, there is little genetic structure among colonies; however, the genetic structure among the colonies off Brazil and its relationship with those in the Caribbean are unknown. In this study, we used mtDNA and microsatellite markers to infer population structure and evolutionary history in a sample of F. magnificens individuals collected in Brazil, Grand Connétable (French Guyana), and Barbuda. Virtually all Brazilian individuals had the same mtDNA haplotype. There was no haplotype sharing between Brazil and the Caribbean, though Grand Connétable shared haplotypes with both regions. A Bayesian clustering analysis using microsatellite data found two genetic clusters: one associated with Barbuda and the other with the Brazilian populations. Grand Connétable was more similar to Barbuda but had ancestry from both clusters, corroborating its "intermediate" position. The Caribbean and Grand Connétable populations showed higher genetic diversity and effective population size compared to the Brazilian population. Overall, our results are in good agreement with an effect of marine winds in isolating the Brazilian meta-population.

Cândido Godói (CG) is a small town in South Brazil, which has the highest prevalence of twin births in Brazil. Recently, a number of studies have shown that p53 plays an important role in reproduction through blastocyst implantation and intra utero embryo survival. Thus, gene polymorphisms in the p53 pathway were investigated in this population.Single nucleotide polymorphisms from five genes in the p53 pathway were investigated, as well as background characteristics of 42 mothers of twins (cases) and 101 mothers of singletons (controls), all residents from CG.Mothers of twins have higher number of pregnancies and higher frequencies of P72 allele at TP53 and T allele at MDM4 genes compared with controls. Logistic regression shows that both TP53 and number of pregnancies maintained their association with twinning (P =0.004 and P =0.002, respectively), with TP53 having a higher odds ratio than number of pregnancies (2.73 versus 1.70, respectively). No interactive effect between TP53 and MDM4 (P =0.966) is observed. As expected, mothers of twins have three times more cases of cancer in their first-degree relatives than control mothers (P =0.011).Our results suggest that the P72 allele of TP53 is a strong risk factor for twinning in CG, while the number of pregnancies and the T allele at MDM4 may represent weaker risk factors. These two alleles are associated with infertility, but the anti-apoptotic effect of low levels of p53 in general, and of the P72 allele in particular, may play a role after implantation, enhancing the chance for a double pregnancy to succeed to term.

Wolbachia are endosymbiotic bacteria that commonly infect arthropods, inducing certain phenotypes in their hosts. So far, no endemic South American species of terrestrial isopods have been investigated for Wolbachia infection. In this work, populations from two species of Balloniscus (B. sellowii and B. glaber) were studied through a diagnostic PCR assay. Fifteen new Wolbachia 16S rDNA sequences were detected. Wolbachia found in both species were generally specific to one population, and five populations hosted two different Wolbachia 16S rDNA sequences. Prevalence was higher in B. glaber than in B. sellowii, but uninfected populations could be found in both species. Wolbachia strains from B. sellowii had a higher genetic variation than those isolated from B. glaber. AMOVA analyses showed that most of the genetic variance was distributed among populations of each species rather than between species, and the phylogenetic analysis suggested that Wolbachia strains from Balloniscus cluster within Supergroup B, but do not form a single monophyletic clade, suggesting multiple infections for this group. Our results highlight the importance of studying Wolbachia prevalence and genetic diversity in Neotropical species and suggest that South American arthropods may harbor a great number of diverse strains, providing an interesting model to investigate the evolution of Wolbachia and its hosts.

Cândido Godói (CG) is a small municipality in South Brazil with approximately 6,000 inhabitants. It is known as the "Twins' Town" due to its high rate of twin births. Recently it was claimed that such high frequency of twinning would be connected to experiments performed by the German Nazi doctor Joseph Mengele. It is known, however, that this town was founded by a small number of families and therefore a genetic founder effect may represent an alternatively explanation for the high twinning prevalence in CG. In this study, we tested specific predictions of the "Nazi's experiment" and of the "founder effect" hypotheses. We surveyed a total of 6,262 baptism records from 1959-2008 in CG catholic churches, and identified 91 twin pairs and one triplet. Contrary to the "Nazi's experiment hypothesis", there is no spurt in twinning between the years (1964-1968) when Mengele allegedly was in CG (P = 0.482). Moreover, there is no temporal trend for a declining rate of twinning since the 1960s (P = 0.351), and no difference in twinning among CG districts considering two different periods: 1927-1958 and 1959-2008 (P = 0.638). On the other hand, the "founder effect hypothesis" is supported by an isonymy analysis that shows that women who gave birth to twins have a higher inbreeding coefficient when compared to women who never had twins (0.0148, 0.0081, respectively, P = 0.019). In summary, our results show no evidence for the "Nazi's experiment hypothesis" and strongly suggest that the "founder effect hypothesis" is a much more likely alternative for explaining the high prevalence of twinning in CG. If this hypothesis is correct, then this community represents a valuable population where genetic factors linked to twinning may be identified.

e-mail: Laurent.excoffier@iee.unibe.ch; International audience; Recent papers have promoted the view that model-based methods in general, and those based on Approximate Bayesian Computation (ABC) in particular, are flawed in a number of ways, and are therefore inappropriate for the analysis of phylogeographic data. These papers further argue that Nested Clade Phylogeographic Analysis (NCPA) offers the best approach in statistical phylogeography. In order to remove the confusion and misconceptions introduced by these papers, we justify and explain the reasoning behind model-based inference. We argue that ABC is a statistically valid approach, alongside other computational statistical techniques that have been successfully used to infer parameters and compare models in population genetics. We also examine the NCPA method and highlight numerous deficiencies, either when used with single or multiple loci. We further show that the ages of clades are carelessly used to infer ages of demographic events, that these ages are estimated under a simple model of panmixia and population stationarity but are then used under different and unspecified models to test hypotheses, a usage the invalidates these testing procedures. We conclude by encouraging researchers to study and use model-based inference in population genetics.

Although there is agreement in that the Bering Strait was the entry point for the initial colonization of the American continent, there is considerable uncertainty regarding the timing and pattern of human migration from Asia to America. In order to perform a statistical assessment of the relative probability of alternative migration scenarios and to estimate key demographic parameters associated with them, we used an approximate Bayesian computation framework to analyze a data set of 401 autosomal microsatellite loci typed in 29 native American populations. A major finding is that a single, discrete, wave of colonization is highly inconsistent with observed levels of genetic diversity. A scenario with two discrete migration waves is also not supported by the data. The current genetic diversity of Amerindian populations is best explained by a third model involving recurrent gene flow between Asia and America, after initial colonization. We estimate that this colonization involved about 100 individuals and occurred some 13,000 years ago, in agreement with well-established archeological data.

Background and objective Uncoupling protein 2 (UCP2) plays a role in controlling reactive oxygen species (ROS) production by mitochondria. As ROS overproduction is related to diabetic retinopathy (DR), UCP2 gene polymorphisms might be involved in the development of this complication. We investigated whether the -866G/A (rs659366), Ala55Val (rs660339) and 45 bp insertion/deletion (Ins/Del) polymorphisms in the UCP2 gene might be associated with proliferative DR (PDR). Design and methods In this case-control study, we analysed 501 type 2 diabetic patients (242 patients with PDR and 259 subjects without any degree of DR) and 196 type 1 diabetic patients (85 cases with PDR and 111 without DR). Haplotypes constructed from the combination of the three UCP2 polymorphisms were inferred using a Bayesian statistical method. Results In the type 2 diabetic group, multivariate analyses confirmed that the haplotype [A Val Ins] was an independent risk factor for PDR when present in one [adjusted odds ratio (aOR) = 2.12; P = 0.006], at least one (aOR = 2.75; P = 0.00001), or two copies (aOR = 5.30; P = 0.00001), suggesting an additive model of inheritance. Nevertheless, in type 1 diabetic patients, the association of this haplotype with PDR was confirmed only when it was present in at least one (aOR = 2.68; P = 0.014) or two copies (aOR = 6.02; P = 0.005). Conclusions The haplotype [A Val Ins] seems to be an important risk factor associated with PDR in both type 2 and 1 diabetic groups.