Statistical properties of new neutrality tests against population growth.

A number of statistical tests for detecting population growth are described. We compared the statistical power of these tests with that of others available in the literature. The tests evaluated fall into three categories: those tests based on the distribution of the mutation frequencies, on the haplotype distribution, and on the mismatch distribution. We found that, for an extensive variety of cases, the most powerful tests for detecting population growth are Fu"s FS test and the newly developed R2 test. The behavior of the R2 test is superior for small sample sizes, whereas FS is better for large sample sizes. We also show that some popular statistics based on the mismatch distribution are very conservative. Key words: population growth, population expansion, coalescent simulations, neutrality tests

Multivariate QST-FST comparisons: a neutrality test for the evolution of the g matrix in structured populations.

Neutrality tests in quantitative genetics provide a statistical framework for the detection of selection on polygenic traits in wild populations. However, the existing method based on comparisons of divergence at neutral markers and quantitative traits (Q(st)-F(st)) suffers from several limitations that hinder a clear interpretation of the results with typical empirical designs. In this article, we propose a multivariate extension of this neutrality test based on empirical estimates of the among-populations (D) and within-populations (G) covariance matrices by MANOVA. A simple pattern is expected under neutrality: D = 2F(st)/(1 - F(st))G, so that neutrality implies both proportionality of the two matrices and a specific value of the proportionality coefficient. This pattern is tested using Flury's framework for matrix comparison [common principal-component (CPC) analysis], a well-known tool in G matrix evolution studies. We show the importance of using a Bartlett adjustment of the test for the small sample sizes typically found in empirical studies. We propose a dual test: (i) that the proportionality coefficient is not different from its neutral expectation [2F(st)/(1 - F(st))] and (ii) that the MANOVA estimates of mean square matrices between and among populations are proportional. These two tests combined provide a more stringent test for neutrality than the classic Q(st)-F(st) comparison and avoid several statistical problems. Extensive simulations of realistic empirical designs suggest that these tests correctly detect the expected pattern under neutrality and have enough power to efficiently detect mild to strong selection (homogeneous, heterogeneous, or mixed) when it is occurring on a set of traits. This method also provides a rigorous and quantitative framework for disentangling the effects of different selection regimes and of drift on the evolution of the G matrix. We discuss practical requirements for the proper application of our test in empirical studies and potential extensions.

Group I introns and associated homing endonuclease genes reveals a clinal structure for Porphyra spiralis var. amplifolia (Bangiales, Rhodophyta) along the Eastern coast of South America

Background: Group I introns are found in the nuclear small subunit ribosomal RNA gene (SSU rDNA) of some species of the genus Porphyra (Bangiales, Rhodophyta). Size polymorphisms in group I introns has been interpreted as the result of the degeneration of homing endonuclease genes (HEG) inserted in peripheral loops of intron paired elements. In this study, intron size polymorphisms were characterized for different Porphyra spiralis var. amplifolia (PSA) populations on the Southern Brazilian coast, and were used to infer genetic relationships and genetic structure of these PSA populations, in addition to cox2-3 and rbcL-S regions. Introns of different sizes were tested qualitatively for in vitro self-splicing. Results: Five intron size polymorphisms within 17 haplotypes were obtained from 80 individuals representing eight localities along the distribution of PSA in the Eastern coast of South America. In order to infer genetic structure and genetic relationships of PSA, these polymorphisms and haplotypes were used as markers for pairwise Fst analyses, Mantel's test and median joining network. The five cox2-3 haplotypes and the unique rbcL-S haplotype were used as markers for summary statistics, neutrality tests Tajima's D and Fu's Fs and for median joining network analyses. An event of demographic expansion from a population with low effective number, followed by a pattern of isolation by distance was obtained for PSA populations with the three analyses. In vitro experiments have shown that introns of different lengths were able to self-splice from pre-RNA transcripts. Conclusion: The findings indicated that degenerated HEGs are reminiscent of the presence of a full-length and functional HEG, once fixed for PSA populations. The cline of HEG degeneration determined the pattern of isolation by distance. Analyses with the other markers indicated an event of demographic expansion from a population with low effective number. The different degrees of degeneration of the HEG do not refrain intron self-splicing. To our knowledge, this was the first study to address intraspecific evolutionary history of a nuclear group I intron; to use nuclear, mitochondrial and chloroplast DNA for population level analyses of Porphyra; and intron size polymorphism as a marker for population genetics.

Land planarians (Platyhelminthes) as a model organism for fine-scale phylogeographic studies: understanding patterns of biodiversity in the Brazilian Atlantic Forest hotspot

The Brazilian Atlantic Forest is one of the richest biodiversity hotspots of the world. Paleoclimatic models have predicted two large stability regions in its northern and central parts, whereas southern regions might have suffered strong instability during Pleistocene glaciations. Molecular phylogeographic and endemism studies show, nevertheless, contradictory results: although some results validate these predictions, other data suggest that paleoclimatic models fail to predict stable rainforest areas in the south. Most studies, however, have surveyed species with relatively high dispersal rates whereas taxa with lower dispersion capabilities should be better predictors of habitat stability. Here, we have used two land planarian species as model organisms to analyse the patterns and levels of nucleotide diversity on a locality within the Southern Atlantic Forest. We find that both species harbour high levels of genetic variability without exhibiting the molecular footprint of recent colonization or population expansions, suggesting a long-term stability scenario. The results reflect, therefore, that paleoclimatic models may fail to detect refugia in the Southern Atlantic Forest, and that model organisms with low dispersal capability can improve the resolution of these models.

African signatures of recent positive selection in human FOXI1

Background: The human FOXI1 gene codes for a transcription factor involved in the physiology of the inner ear, testis, and kidney. Using three interspecies comparisons, it has been suggested that this may be a gene underhuman-specific selection. We sought to confirm this finding by using an extended set of orthologous sequences.Additionally, we explored for signals of natural selection within humans by sequencing the gene in 20 Europeans,20 East Asians and 20 Yorubas and by analysing SNP variation in a 2 Mb region centered on FOXI1 in 39worldwide human populations from the HGDP-CEPH diversity panel.Results: The genome sequences recently available from other primate and non-primate species showed that FOXI1divergence patterns are compatible with neutral evolution. Sequence-based neutrality tests were not significant inEuropeans, East Asians or Yorubas. However, the Long Range Haplotype (LRH) test, as well as the iHS and XP-Rsbstatistics revealed significantly extended tracks of homozygosity around FOXI1 in Africa, suggesting a recentepisode of positive selection acting on this gene. A functionally relevant SNP, as well as several SNPs either on theputatively selected core haplotypes or with significant iHS or XP-Rsb values, displayed allele frequencies stronglycorrelated with the absolute geographical latitude of the populations sampled.Conclusions: We present evidence for recent positive selection in the FOXI1 gene region in Africa. Climate mightbe related to this recent adaptive event in humans. Of the multiple functions of FOXI1, its role in kidney-mediatedwater-electrolyte homeostasis is the most obvious candidate for explaining a climate-related adaptation.

Schmidtea mediterranea phylogeography: an old species surviving on a few mediterranean islands?

Schmidtea mediterranea (Platyhelminthes, Tricladida, Continenticola) is found in scattered localities on a few islands and in coastal areas of the western Mediterranean. Although S. mediterranea is the object of many regeneration studies, little is known about its evolutionary history. Its present distribution has been proposed to stem from the fragmentation and migration of the Corsica-Sardinia microplate during the formation of the western Mediterranean basin, which implies an ancient origin for the species. To test this hypothesis, we obtained a large number of samples from across its distribution area. Using known and new molecular markers and, for the first time in planarians, a molecular clock, we analysed the genetic variability and demographic parameters within the species and between its sexual and asexual populations to estimate when they diverged. Results: A total of 2 kb from three markers (COI, CYB and a nuclear intron N13) was amplified from ~200 specimens. Molecular data clustered the studied populations into three groups that correspond to the west, central and southeastern geographical locations of the current distribution of S. mediterranea. Mitochondrial genes show low haplotype and nucleotide diversity within populations but demonstrate higher values when all individuals are considered. The nuclear marker shows higher values of genetic diversity than the mitochondrial genes at the population level, but asexual populations present lower variability than the sexual ones. Neutrality tests are significant for some populations. Phylogenetic and dating analyses show the three groups to be monophyletic, with the west group being the basal group. The time when the diversification of the species occurred is between ~20 and ~4 mya, although the asexual nature of the western populations could have affected the dating analyses. Conclusions: S. mediterranea is an old species that is sparsely distributed in a harsh habitat, which is probably the consequence of the migration of the Corsica-Sardinia block. This species probably adapted to temperate climates in the middle of a changing Mediterranean climate that eventually became dry and hot. These data also suggest that in the mainland localities of Europe and Africa, sexual individuals of S. mediterranea are being replaced by asexual individuals that are either conspecific or are from other species that are better adapted to the Mediterranean climate.

Population structure of the migratory fish Prochilodus lineatus (Characiformes) from rio Grande basin (Brazil), an area fragmented by dams

1. Prochilodus lineatus (Prochilodontidae, Characiformes) is a migratory species of great economic importance both in fisheries and aquaculture that is found throughout the Jacui, Paraiba do Sul, Parana, Paraguay and Uruguay river basins in South America. Earlier population studies of P. lineatus in the rio Grande basin (Parana basin) indicated the existence of a single population; however, the range of this species has been fragmented by the construction of several dams. Such dams modified the environmental conditions and could have constrained the reproductive migration of P. lineatus, possibly leading to changes in the population genetic structure. 2. In order to evaluate how genetic diversity is allocated in the rio Grande basin, 141 specimens of P. lineatus from eight collection sites were analysed using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) with 15 restriction enzymes. 3. Forty-six haplotypes were detected, and 70% of them are restricted. The mean genetic variability indexes (h = 0.7721 and pi = 1.6%) were similar to those found in natural populations with a large effective size. Fst and Exact Test values indicated a lack of structuring among the samples, and the model of isolation by distance was tested and rejected. 4. The haplotype network indicated that this population of P. lineatus has been maintained as a single variable stock with some differences in the genetic composition (haplotypes) between samples. Indications of population expansion were detected, and this finding was supported by neutrality tests and mismatch distribution analyses. 5. The present study focused on regions between dams to serve as a parameter for further evaluations of genetic variability and the putative impact of dams and repopulation programmes in natural populations of P. lineatus. Copyright (C) 2011 John Wiley & Sons, Ltd.

Balanced polymorphism in bottlenecked populations: The case of the CCR5 5 ` cis-regulatory region in Amazonian Amerindians

The 5` cis-regulatory region of the CCR5 gene exhibits a strong signature of balancing selection in several human populations. Here we analyze the polymorphism of this region in Amerindians from Amazonia, who have a complex demographic history, including recent bottlenecks that are known to reduce genetic variability. Amerindians show high nucleotide diversity (pi = 0.27%) and significantly positive Tajima`s D, and carry haplotypes associated with weak and strong gene expression. To evaluate whether these signatures of balancing selection could be explained by demography, we perform neutrality tests based on empiric and simulated data. The observed Tajima`s D was higher than that of other world populations: higher than that found for 18 noncoding regions of South Amerindians, and higher than 99.6% of simulated genealogies, which assume nonequilibrium conditions. Moreover, comparing Amerindians and Asians, the Fst for CCR5 cis-regulatory region was unusually low, in relation to neutral markers. These findings indicate that, despite their complex demographic history, South Amerindians carry a detectable signature of selection on the CCR5 cis-regulatory region. (C) 2010 American Society for Histocompatibility and Immunogenetics. Published by Elsevier Inc. All rights reserved.

Limited global diversity of the Plasmodium vivax merozoite surface protein 4 gene

Merozoite surface proteins (MSPs) of the malaria parasites are major candidates for vaccine development targeting asexual blood stages. However, the diverse antigenic repertoire of these antigens that induce strain-specific protective immunity in human is a major challenge for vaccine design and often determines the efficacy of a vaccine. Here we further assessed the genetic diversity of Plasmodium vivax MSP4 (PvMSP4) protein using 195 parasite samples collected mostly from Thailand, Indonesia and Brazil. Overall, PvMSP4 is highly conserved with only eight amino acid substitutions. The majority of the haplotype diversity was restricted to the two short tetrapeptide repeat arrays in exon 1 and 2, respectively. Selection and neutrality tests indicated that exon 1 and the entire coding region of PvMSP4 were under purifying selection. Despite the limited nucleotide polymorphism of PvMSP4, significant genetic differentiation among the three major parasite populations was detected. Moreover, microgeographical heterogeneity was also evident in the parasite populations from different endemic areas of Thailand. (C) 2009 Elsevier B.V. All rights reserved.

Genetic diversity of the HLA-G coding region in Amerindian populations from the Brazilian Amazon: a possible role of natural selection

HLA-G has an important role in the modulation of the maternal immune system during pregnancy, and evidence that balancing selection acts in the promoter and 3′UTR regions has been previously reported. To determine whether selection acts on the HLA-G coding region in the Amazon Rainforest, exons 2, 3 and 4 were analyzed in a sample of 142 Amerindians from nine villages of five isolated tribes that inhabit the Central Amazon. Six previously described single-nucleotide polymorphisms (SNPs) were identified and the Expectation-Maximization (EM) and PHASE algorithms were used to computationally reconstruct SNP haplotypes (HLA-G alleles). A new HLA-G allele, which originated in Amerindian populations by a crossing-over event between two widespread HLA-G alleles, was identified in 18 individuals. Neutrality tests evidenced that natural selection has a complex part in the HLA-G coding region. Although balancing selection is the type of selection that shapes variability at a local level (Native American populations), we have also shown that purifying selection may occur on a worldwide scale. Moreover, the balancing selection does not seem to act on the coding region as strongly as it acts on the flanking regulatory regions, and such coding signature may actually reflect a hitchhiking effect.Genes and Immunity advance online publication, 3 October 2013; doi:10.1038/gene.2013.47.

Genetic variability and demographic history of Heliothis virescens (Lepidoptera: Noctuidae) populations from Brazil inferred by mtDNA sequences

Intra-and inter-population genetic variability and the demographic history of Heliothis virescens (F.) populations were evaluated by using mtDNA markers (coxI, coxII and nad6) with samples from the major cotton-and soybean-producing regions in Brazil in the growing seasons 2007/08, 2008/09 and 2009/10. AMOVA indicated low and non-significant genetic structure, regardless of geographical scale, growing season or crop, with most of genetic variation occurring within populations. Clustering analyzes also indicated low genetic differentiation. The haplotype network obtained with combined datasets resulted in 35 haplotypes, with 28 exclusive occurrences, four of them sampled only from soybean fields. The minimum spanning network showed star-shaped structures typical of populations that underwent a recent demographic expansion. The recent expansion was supported by other demographic analyzes, such as the Bayesian skyline plot, the unimodal distribution of paired differences among mitochondrial sequences, and negative and significant values of neutrality tests for the Tajima's D and Fu's F-S parameters. In addition, high values of haplotype diversity ((H) over cap) and low values of nucleotide diversity (pi), combined with a high number of low frequency haplotypes and values of theta(pi)<theta(W), suggested a recent demographic expansion of H. virescens populations in Brazil. This demographic event could be responsible for the low genetic structure currently found; however, haplotypes present uniquely at the same geographic regions and from one specific host plant suggest an initial differentiation among H. virescens populations within Brazil.

Non-classical HLA-E gene variability in Brazilians: a nearly invariable locus surrounded by the most variable genes in the human genome

The non-classical human leukocyte antigen (HLA) class I genes present a very low rate of variation. So far, only 10 HLA-E alleles encoding three proteins have been described, but only two are frequently found in worldwide populations. Because of its historical background, Brazilians are very suitable for population genetic studies. Therefore, 104 bone marrow donors from Brazil were evaluated for HLA-E exons 14. Seven variation sites were found, including two known single nucleotide polymorphisms (SNPs) at positions +424 and +756 and five new SNPs at positions +170 (intron 1), +1294 (intron 3), +1625, +1645 and +1857 (exon 4). Haplotyping analysis did show eight haplotypes, three of them known as E*01:01:01, E*01:03:01 and E*01:03:02:01 and five HLA-E new alleles that carry the new variation sites. The HLA-E*01:01:01 allele was the predominant haplotype (62.50%), followed by E*01:03:02:01 (24.52%). Selective neutrality tests have disclosed an interesting pattern of selective pressures in which balancing selection is probably shaping allele frequency distributions at an SNP at exon 3 (codon 107), sequence diversity at exon 4 and the non-coding regions is facing significant purifying pressure. Even in an admixed population such as the Brazilian one, the HLA-E locus is very conserved, presenting few polymorphic SNPs in the coding region.

Phylogeography of an Atlantic forest passerine reveals demographic stability through the last glacial maximum

In this study we analyzed the phylogeographic pattern and historical demography of an endemic Atlantic forest (AF) bird, Basileuterus leucoblepharus, and test the influence of the last glacial maximum (LGM) on its population effective size using coalescent simulations. We address two main questions: (i) Does B. leucoblepharus present population genetic structure congruent with the patterns observed for other AF organisms? (ii) How did the LGM affect the effective population size of B. leucoblepharus? We sequenced 914 bp of the mitochondrial gene cytochrome b and 512 bp of the nuclear intron 5 of beta-fibrinogen of 62 individuals from 15 localities along the AF. Both molecular markers revealed no genetic structure in B. leucoblepharus. Neutrality tests based on both loci showed significant demographic expansion. The extended Bayesian skyline plot showed that the species seems to have experienced demographic expansion starting around 300,000 years ago, during the late Pleistocene. This date does not coincide with the LGM and the dynamics of population size showed stability during the LGM. To further test the effect of the LGM on this species, we simulated seven demographic scenarios to explore whether populations suffered specific bottlenecks. The scenarios most congruent with our data were population stability during the LGM with bottlenecks older than this period. This is the first example of an AF organism that does not show phylogeographic breaks caused by vicariant events associated to climate change and geotectonic activities in the Quaternary. Differential ecological, environmental tolerances and habitat requirements are possibly influencing the different evolutionary histories of these organisms. Our results show that the history of organism diversification in this megadiverse Neotropical forest is complex. Crown Copyright (c) 2012 Published by Elsevier Inc. All rights reserved.

PHYLOGENETIC ANALYSES, GENE DOSAGE EFFECT AND SELECTION IN MEMBERS OF THE LATERAL BOUNDARY DOMAIN (LBD) GENES TRANSCRIPTION FACTOR FAMILY

Wood formation is an economically and environmentally important process and has played a significant role in the evolution of terrestrial plants. Despite its significance, the molecular underpinnings of the process are still poorly understood. We have previously shown that four Lateral Boundary Domain (LBD) transcription factors have important roles in the regulation of wood formation with two (LBD1 and LBD4) involved in secondary phloem and ray cell development and two (LBD15 and LBD18) in secondary xylem formation. Here, we used comparative phylogenetic analyses to test potential roles of the four LBD genes in the evolution of woodiness. We studied the copy number and variation in DNA and amino acid sequences of the four LBDs in a wide range of woody and herbaceous plant taxa with fully sequenced and annotated genomes. LBD1 showed the highest gene copy number across the studied species, and LBD1 gene copy number was strongly and significantly correlated with the level of ray seriation. The lianas, cucumber and grape, with multiseriate ray cells showed the highest gene copy number (12 and 11, respectively). Because lianas’ growth habit requires significant twisting and bending, the less lignified ray parenchyma cells likely facilitate stem flexibility and maintenance of xylem conductivity. We further demonstrate conservation of amino acids in the LBD18 protein sequences that are specific to woody taxa. Neutrality tests showed evidence for strong purifying selection on these gene regions across various orders, indicating adaptive convergent evolution of LBD18. Structural modeling demonstrates that the conserved amino acids have a significant impact on the tertiary protein structure and thus are likely of significant functional importance.

Mitochondrial demographic history of the Egyptian mongoose (Herpestes ichneumon), an expanding carnivore in the Iberian Peninsula

Describing the genetic patterns and the demographic history of expanding species is essential for providing insights into the processes linked with range dynamics. We analysed the mitochondrial diversity of the Egyptian mongoose (Herpestes ichneumon) across the Iberian Peninsula, where the species is currently expanding northwest. A total of 242 individuals were analysed, together with nine representatives from the North African dispersal source. Haplotype segregation and strong differentiation between Iberian and North African populations confirmed the longterm presence of the species in the Iberian Peninsula. The distribution of mitochondrial diversity fitted the pattern of a historically diversified population in southern Iberia, from which the recent dispersals into northern areas may have occurred. Higher levels of haplotype and nucleotide diversities in the northern areas, together with the heterogeneous distribution of pairwise population differentiations and the weak signal for isolationbydistance suggest the existence of longdispersal migrants across the Iberian Peninsula. Sudden and spatial expansion scenarios of H. ichneumon in the Iberian Peninsula were supported by mismatch analysis and marginally supported by neutrality tests. However, the precise time of occurrence of the detected expansion remains unclear. Future studies should incorporate additional markers in order to further clarify the population dynamics of the Egyptian mongoose in its Iberian range.