The Genomic Signature of Population Reconnection Following Isolation: From Theory to HIV.

Ease of worldwide travel provides increased opportunities for organisms not only to colonize new environments but also to encounter related but diverged populations. Such events of reconnection and secondary contact of previously isolated populations are widely observed at different time scales. For example, during the quaternary glaciation, sea water level fluctuations caused temporal isolation of populations, often to be followed by secondary contact. At shorter time scales, population isolation and reconnection of viruses are commonly observed, and such events are often associated with epidemics and pandemics. Here, using coalescent theory and simulations, we describe the temporal impact of population reconnection after isolation on nucleotide differences and the site frequency spectrum, as well as common summary statistics of DNA variation. We identify robust genomic signatures of population reconnection after isolation. We utilize our development to infer the recent evolutionary history of human immunodeficiency virus 1 (HIV-1) in Asia and South America, successfully retrieving the successive HIV subtype colonization events in these regions. Our analysis reveals that divergent HIV-1 subtype populations are currently admixing in these regions, suggesting that HIV-1 may be undergoing a process of homogenization, contrary to popular belief.

Tracking Invasion Histories in the Sea: Facing Complex Scenarios Using Multilocus Data

In recent years, new analytical tools have allowed researchers to extract historical information contained in molecular data, which has fundamentally transformed our understanding of processes ruling biological invasions. However, the use of these new analytical tools has been largely restricted to studies of terrestrial organisms despite the growing recognition that the sea contains ecosystems that are amongst the most heavily affected by biological invasions, and that marine invasion histories are often remarkably complex. Here, we studied the routes of invasion and colonisation histories of an invasive marine invertebrate Microcosmus squamiger (Ascidiacea) using microsatellite loci, mitochondrial DNA sequence data and 11 worldwide populations. Discriminant analysis of principal components, clustering methods and approximate Bayesian computation (ABC) methods showed that the most likely source of the introduced populations was a single admixture event that involved populations from two genetically differentiated ancestral regions - the western and eastern coasts of Australia. The ABC analyses revealed that colonisation of the introduced range of M. squamiger consisted of a series of non-independent introductions along the coastlines of Africa, North America and Europe. Furthermore, we inferred that the sequence of colonisation across continents was in line with historical taxonomic records - first the Mediterranean Sea and South Africa from an unsampled ancestral population, followed by sequential introductions in California and, more recently, the NE Atlantic Ocean. We revealed the most likely invasion history for world populations of M. squamiger, which is broadly characterized by the presence of multiple ancestral sources and non-independent introductions within the introduced range. The results presented here illustrate the complexity of marine invasion routes and identify a cause-effect relationship between human-mediated transport and the success of widespread marine non-indigenous species, which benefit from stepping-stone invasions and admixture processes involving different sources for the spread and expansion of their range.

Divergence with gene flow and fine-scale phylogeographical structure in the wedge-billed woodcreeper, Glyphorynchus spirurus, a Neotropical rainforest bird.

Determining the relative roles of vicariance and selection in restricting gene flow between populations is of central importance to the evolutionary process of population divergence and speciation. Here we use molecular and morphological data to contrast the effect of isolation (by mountains and geographical distance) with that of ecological factors (altitudinal gradients) in promoting differentiation in the wedge-billed woodcreeper, Glyphorynchus spirurus, a tropical forest bird, in Ecuador. Tarsus length and beak size increased relative to body size with altitude on both sides of the Andes, and were correlated with the amount of moss on tree trunks, suggesting the role of selection in driving adaptive divergence. In contrast, molecular data revealed a considerable degree of admixture along these altitudinal gradients, suggesting that adaptive divergence in morphological traits has occurred in the presence of gene flow. As suggested by mitochondrial DNA sequence data, the Andes act as a barrier to gene flow between ancient subspecific lineages. Genome-wide amplified fragment length polymorphism markers reflected more recent patterns of gene flow and revealed fine-scale patterns of population differentiation that were not detectable with mitochondrial DNA, including the differentiation of isolated coastal populations west of the Andes. Our results support the predominant role of geographical isolation in driving genetic differentiation in G. spirurus, yet suggest the role of selection in driving parallel morphological divergence along ecological gradients.

The genetic basis of color-related local adaptation in a ring-like colonization around the Mediterranean.

Uncovering the genetic basis of phenotypic variation and the population history under which it established is key to understand the trajectories along which local adaptation evolves. Here, we investigated the genetic basis and evolutionary history of a clinal plumage color polymorphism in European barn owls (Tyto alba). Our results suggest that barn owls colonized the Western Palearctic in a ring-like manner around the Mediterranean and meet in secondary contact in Greece. Rufous coloration appears to be linked to a recently evolved nonsynonymous-derived variant of the melanocortin 1 receptor (MC1R) gene, which according to quantitative genetic analyses evolved under local adaptation during or following the colonization of Central Europe. Admixture patterns and linkage disequilibrium between the neutral genetic background and color found exclusively within the secondary contact zone suggest limited introgression at secondary contact. These results from a system reminiscent of ring species provide a striking example of how local adaptation can evolve from derived genetic variation.

Genomics of clinal variation in Drosophila: disentangling the interactions of selection and demography.

Clines in phenotypes and genotype frequencies across environmental gradients are commonly taken as evidence for spatially varying selection. Classical examples include the latitudinal clines in various species of Drosophila, which often occur in parallel fashion on multiple continents. Today, genomewide analysis of such clinal systems provides a fantastic opportunity for unravelling the genetics of adaptation, yet major challenges remain. A well-known but often neglected problem is that demographic processes can also generate clinality, independent of or coincident with selection. A closely related issue is how to identify true genic targets of clinal selection. In this issue of Molecular Ecology, three studies illustrate these challenges and how they might be met. Bergland et al. report evidence suggesting that the well-known parallel latitudinal clines in North American and Australian D. melanogaster are confounded by admixture from Africa and Europe, highlighting the importance of distinguishing demographic from adaptive clines. In a companion study, Machado et al. provide the first genomic comparison of latitudinal differentiation in D. melanogaster and its sister species D. simulans. While D. simulans is less clinal than D. melanogaster, a significant fraction of clinal genes is shared between both species, suggesting the existence of convergent adaptation to clinaly varying selection pressures. Finally, by drawing on several independent sources of evidence, Bo?ičević et al. identify a functional network of eight clinal genes that are likely involved in cold adaptation. Together, these studies remind us that clinality does not necessarily imply selection and that separating adaptive signal from demographic noise requires great effort and care.

Complex patterns of global spread in invasive insects: eco-evolutionary and management consequences

The advent of simple and affordable tools for molecular identification of novel insect invaders and assessment of population diversity has changed the face of invasion biology in recent years. The widespread application of these tools has brought with it an emerging understanding that patterns in biogeography, introduction history and subsequent movement and spread of many invasive alien insects are far more complex than previously thought. We reviewed the literature and found that for a number of invasive insects, there is strong and growing evidence that multiple introductions, complex global movement, and population admixture in the invaded range are commonplace. Additionally, historical paradigms related to species and strain identities and origins of common invaders are in many cases being challenged. This has major consequences for our understanding of basic biology and ecology of invasive insects and impacts quarantine, management and biocontrol programs. In addition, we found that founder effects rarely limit fitness in invasive insects and may benefit populations (by purging harmful alleles or increasing additive genetic variance). Also, while phenotypic plasticity appears important post-establishment, genetic diversity in invasive insects is often higher than expected and increases over time via multiple introductions. Further, connectivity among disjunct regions of global invasive ranges is generally far higher than expected and is often asymmetric, with some populations contributing disproportionately to global spread. We argue that the role of connectivity in driving the ecology and evolution of introduced species with multiple invasive ranges has been historically underestimated and that such species are often best understood in a global context.

Genomic evidence for adaptive inversion clines in Drosophila melanogaster.

Clines in chromosomal inversion polymorphisms-presumably driven by climatic gradients-are common but there is surprisingly little evidence for selection acting on them. Here we address this long-standing issue in Drosophila melanogaster by using diagnostic single nucleotide polymorphism (SNP) markers to estimate inversion frequencies from 28 whole-genome Pool-seq samples collected from 10 populations along the North American east coast. Inversions In(3L)P, In(3R)Mo, and In(3R)Payne showed clear latitudinal clines, and for In(2L)t, In(2R)NS, and In(3R)Payne the steepness of the clinal slopes changed between summer and fall. Consistent with an effect of seasonality on inversion frequencies, we detected small but stable seasonal fluctuations of In(2R)NS and In(3R)Payne in a temperate Pennsylvanian population over 4 years. In support of spatially varying selection, we observed that the cline in In(3R)Payne has remained stable for >40 years and that the frequencies of In(2L)t and In(3R)Payne are strongly correlated with climatic factors that vary latitudinally, independent of population structure. To test whether these patterns are adaptive, we compared the amount of genetic differentiation of inversions versus neutral SNPs and found that the clines in In(2L)t and In(3R)Payne are maintained nonneutrally and independent of admixture. We also identified numerous clinal inversion-associated SNPs, many of which exhibit parallel differentiation along the Australian cline and reside in genes known to affect fitness-related traits. Together, our results provide strong evidence that inversion clines are maintained by spatially-and perhaps also temporally-varying selection. We interpret our data in light of current hypotheses about how inversions are established and maintained.

Detectings low introgression of invasive alleles in an extensively restocked game bird

Interbreeding of two species in the wild implies introgression of alleles from one species into the other only when admixed individuals survive and successfully backcross with the parental species. Consequently, estimating the proportion of first generation hybrids in a population may not inform about the evolutionary impact of hybridization. Samples obtained over a long time span may offer a more accurate view of the spreading of introgressed alleles in a species" gene pool. Common quail (Coturnix coturnix) populations in Europe have been restocked extensively with farm quails of hybrid origin (crosses with Japanese quails, C. japonica). We genetically monitored a common quail population over 15 years to investigate whether genetic introgression is occurring and used simulations to investigate our power to detect it. Our results revealed that some introgression has occurred, but we did not observe a significant increase over time in the proportion of admixed individuals. However, simulations showed that the degree of admixture may be larger than anticipated due to the limited power of analyses over a short time span, and that observed data was compatible with a low rate of introgression, probably resulting from reduced fitness of admixed individuals. Simulations predicted this could result in extensive admixture in the near future.

A simple model for the estimation of congenital malformation frequency in racially mixed populations

A simple model is proposed, using the method of maximum likelihood to estimate malformation frequencies in racial groups based on data obtained from hospital services. This model uses the proportions of racial admixture, and the observed malformation frequency. It was applied to two defects: postaxial polydactyly and cleft lip, the frequencies of which are recognizedly heterogeneous among racial groups. The frequencies estimated in each racial group were those expected for these malformations, which proves the applicability of the method.

CCR5 genotype and plasma ß-chemokine concentration of Brazilian HIV-infected individuals

The 32-bp deletion in the HIV-1 co-receptor CCR5 confers a high degree of resistance to HIV-1 infection in homozygous individuals for the deleted allele and partial protection against HIV-1 during disease progression in heterozygotes. Natural ligands for CCR5, MIP-1alpha, MIP-1ß and RANTES, have been shown to inhibit HIV replication in CD4+ T cells. In the present study, we examined the CCR5 genotype by PCR and the plasma levels of RANTES and MIP-1alpha by ELISA among blood donors (N = 26) and among HIV-1-infected individuals (N = 129). The control group consisted of healthy adult volunteers and HIV-1-infected subjects were an asymptomatic and heterogeneous group of individuals with regard to immunologic and virologic markers of HIV-1 disease. The frequency of the CCR5 mutant allele (delta32ccr5) in this population was 0.032; however, no delta32ccr5 homozygote was detected. These results could be related to the intense ethnic admixture of the Brazilian population. There was no correlation between circulating ß-chemokines (MIP-1alpha, RANTES) and viral load in HIV-infected individuals. RANTES concentrations in plasma samples from HIV+ patients carrying the homozygous CCR5 allele (CCR5/CCR5) (28.23 ng/ml) were higher than in the control samples (16.07 ng/ml; P<0.05); however, this HIV+ patient group (mean 26.23 pg/ml) had significantly lower concentrations of MIP-1alpha than those observed in control samples (mean 31.20 pg/ml; P<0.05). Both HIV-1-infected and uninfected individuals heterozygous for the delta32ccr5 allele had significantly lower concentrations of circulating RANTES (mean 16.07 and 6.11 ng/ml, respectively) than CCR5/CCR5 individuals (mean 28.23 and 16.07 ng/ml, respectively; P<0.05). These findings suggest that the CCR5 allele and ß-chemokine production may affect the immunopathogenesis of HIV-1.

Ethnicity and glutathione S-transferase (GSTM1/GSTT1) polymorphisms in a Brazilian population

The distribution of polymorphisms related to glutathione S-transferases (GST) has been described in different populations, mainly for white individuals. We evaluated the distribution of GST mu (GSTM1) and theta (GSTT1) genotypes in 594 individuals, by multiplex PCR-based methods, using amplification of the exon 7 of CYP1A1 gene as an internal control. In São Paulo, 233 whites, 87 mulattos, and 137 blacks, all healthy blood-donor volunteers, were tested. In Bahia, where black and mulatto populations are more numerous, 137 subjects were evaluated. The frequency of the GSTM1 null genotype was significantly higher among whites (55.4%) than among mulattos (41.4%; P = 0.03) and blacks (32.8%; P < 0.0001) from São Paulo, or Bahian subjects in general (35.7%; P = 0.0003). There was no statistically different distribution among any non-white groups. The distribution of GSTT1 null genotype among groups did not differ significantly. The agreement between self-reported and interviewer classification of skin color in the Bahian group was low. The interviewer classification indicated a gradient of distribution of the GSTM1 null genotype from whites (55.6%) to light mulattos (40.4%), dark mulattos (32.0%) and blacks (28.6%). However, any information about race or ethnicity should be considered with caution regarding the bias introduced by different data collection techniques, specially in countries where racial admixture is intense, and ethnic definition boundaries are loose. Because homozygous deletions of GST gene might be associated with cancer risk, a better understanding of chemical metabolizing gene distribution can contribute to risk assessment of humans exposed to environmental carcinogens.

Prevalence of deltaF508, G551D, G542X, and R553X mutations among cystic fibrosis patients in the North of Brazil

Cystic fibrosis (CF) is the most common genetic disease among Caucasians and is rare among sub-Saharan Africans. The Brazilian population is not ethnically homogeneous but it is the result of three-way ethnic admixture of Europeans, Africans and Amerindians in varying proportions, depending on the region. In the present study, we investigated 33 patients who had been diagnosed and are currently under treatment for CF at the University Hospital João de Barros Barreto, Belém, Pará State. The molecular analysis for G542X, G551D and R553X mutations was performed by PCR followed by RFLP using BstNI, HincII and MboI, respectively, in polyacrylamide gel eletrophoresis and stained with AgNO3. ThedeltaF508 mutation (a deletion of 3 bp) was only analyzed by polyacrylamide gel electrophoresis and stained with AgNO3. Each sample was analyzed for regions of interest in the CFTR gene using amplified by PCR and specific primers. The deltaF508 and G551D mutations presented frequencies of 22.7 and 3%, respectively. In 74.3% of the remaining patients, none of the mutations investigated was found. The present study characterized in a sample of patients with an established clinical diagnosis of CF (asthma, repeated bronchopneumonia, disorders of nutritional status, etc.) the most frequent mutation ( deltaF508) in the North region of Brazil and is also the first report of the G551D mutation. In spite of the wide spectrum of CF mutations and the heterogeneous ethnic origin of the Amazon population, the molecular diagnosis is a helpful additional tool for the diagnosis and treatment of CF patients.

Association of polymorphisms at the ADIPOR1 regulatory region with type 2 diabetes and body mass index in a Brazilian population with European or African ancestry

Association studies between ADIPOR1 genetic variants and predisposition to type 2 diabetes (DM2) have provided contradictory results. We determined if two single nucleotide polymorphisms (SNP c.-8503G>A and SNP c.10225C>G) in regulatory regions of ADIPOR1 in 567 Brazilian individuals of European (EA; N = 443) or African (AfA; N = 124) ancestry from rural (quilombo remnants; N = 439) and urban (N = 567) areas. We detected a significant effect of ethnicity on the distribution of the allelic frequencies of both SNPs in these populations (EA: -8503A = 0.27; AfA: -8503A = 0.16; P = 0.001 and EA: 10225G = 0.35; AfA: 10225G = 0.51; P < 0.001). Neither of the polymorphisms were associated with DM2 in the case-control study in EA (SNP c.-8503G>A: DM2 group -8503A = 0.26; control group -8503A = 0.30; P = 0.14/SNP 10225C>G: DM2 group 10225G = 0.37; control group 10225G = 0.32; P = 0.40) and AfA populations (SNP c.-8503G>A: DM2 group -8503A = 0.16; control group -8503A = 0.15; P = 0.34/SNP 10225C>G: DM2 group 10225G = 0.51; control group 10225G = 0.52; P = 0.50). Similarly, none of the polymorphisms were associated with metabolic/anthropometric risk factors for DM2 in any of the three populations, except for HDL cholesterol, which was significantly higher in AfA heterozygotes (GC = 53.75 ± 17.26 mg/dL) than in homozygotes. We conclude that ADIPOR1 polymorphisms are unlikely to be major risk factors for DM2 or for metabolic/anthropometric measurements that represent risk factors for DM2 in populations of European and African ancestries.

DNA tests probe the genomic ancestry of Brazilians

We review studies from our laboratories using different molecular tools to characterize the ancestry of Brazilians in reference to their Amerindian, European and African roots. Initially we used uniparental DNA markers to investigate the contribution of distinct Y chromosome and mitochondrial DNA lineages to present-day populations. High levels of genetic admixture and strong directional mating between European males and Amerindian and African females were unraveled. We next analyzed different types of biparental autosomal polymorphisms. Especially useful was a set of 40 insertion-deletion polymorphisms (indels) that when studied worldwide proved exquisitely sensitive in discriminating between Amerindians, Europeans and Sub-Saharan Africans. When applied to the study of Brazilians these markers confirmed extensive genomic admixture, but also demonstrated a strong imprint of the massive European immigration wave in the 19th and 20th centuries. The high individual ancestral variability observed suggests that each Brazilian has a singular proportion of Amerindian, European and African ancestries in his mosaic genome. In Brazil, one cannot predict the color of persons from their genomic ancestry nor the opposite. Brazilians should be assessed on a personal basis, as 190 million human beings, and not as members of color groups.

Prevalence of the STK15 F31I polymorphism and its relationship with mammographic density

Several studies have identified the single nucleotide polymorphism STK15 F31I as a low-penetrance risk allele for breast cancer, but its prevalence and risk association in the Brazilian population have not been determined. The goal of this study was to identify the frequency of this polymorphism in the Brazilian setting. Considering the high degree of admixture of our population, it is of fundamental importance to validate the results already reported in the literature and also to verify the relationship between this variant and breast cancer risk. A total of 750 women without breast cancer were genotyped using the TaqMan PCR assay for STK15 F31I polymorphism. Clinical information was obtained from review of the medical records and mammographic density from the images obtained using the BI-RADS System. The estimated risk of developing cancer was calculated according to the Gail model. The genotypic frequencies observed in this study were 4.5, 38.7, and 56.6%, respectively, for the STK15 F31I AA, AT and TT genotypes. The AT and AA genotypes were encountered significantly more often in premenopausal women with moderately dense, dense and heterogeneously dense breast tissue (P = 0.023). In addition, the presence of the TT genotype was significantly associated with age at menarche ≥12 years (P = 0.023). High mammographic density, associated with increased breast cancer risk, was encountered more frequently in premenopausal women with the risk genotypes STK15 F31I AA and AT. The genotypic frequencies observed in our Brazilian sample were similar to those described in other predominantly European populations.