Contrasting genetic patterns and population histories in three threatened redfin species (Cyprinidae) from the Olifants River system, western South Africa

Differentiation occurred in different areas of the same river system in Pseudobarbus phlegethon and two sister species of the serrated redfin lineage ( Barbus calidus and Barbus erubescens ) of the Western Cape region of South Africa. In an analysis of 27 allozyme loci a deep divergence was found within P. phlegethon between the Olifants and Doring catchments (seven fixed allelic differences; 0·338 <  D  < 0·366; F _ST = 0·925). In contrast, speciation in the serrated redfin lineage occurred within the Doring catchment between B. calidus and B. erubescens (one fixed allelic difference; 0·008 <  D  < 0·052; F _ST = 0·760). Different niche preferences, behaviour and morphology probably played an important role in the formation of these contrasting patterns of genetic structure and suggest that B. calidus may have been better able to disperse between the Olifants and Doring catchments (no fixed allelic differences; 0 <  D  < 0·040; F _ST = 0·281) than P. phlegethon .     

Contrasting patterns of population connectivity between regions in a commercially important mollusc Haliotis rubra: integrating population genetics,genomics and marine LiDAR data

Estimating contemporary genetic structure and population connectivity in marine species is challenging, often compromised by genetic markers that lack adequate sensitivity, and unstructured sampling regimes. We show how these limitations can be overcome via the integration of modern genotyping methods and sampling designs guided by LiDAR and SONAR data sets. Here we explore patterns of gene flow and local genetic structure in a commercially harvested abalone species (Haliotis rubra) from southeastern Australia, where the viability of fishing stocks is believed to be dictated by recruitment from local sources. Using a panel of microsatellite and genomewide SNP markers, we compare allele frequencies across a replicated hierarchical sampling area guided by bathymetric LiDAR imagery. Results indicate high levels of gene flow and no significant genetic structure within or between benthic reef habitats across 1400 km of coastline. These findings differ to those reported for other regions of the fishery indicating that larval supply is likely to be spatially variable, with implications for management and long‐term recovery from stock depletion. The study highlights the utility of suitably designed genetic markers and spatially informed sampling strategies for gaining insights into recruitment patterns in benthic marine species, assisting in conservation planning and sustainable management of fisheries.     

Contrasting patterns of population subdivision and historical demography in three western Mediterranean lizard species inferred from mitochondrial DNA variation

Pleistocene climatic oscillations were a major force shaping genetic variability in many taxa. We analyse the relative effects of the ice ages across a latitudinal gradient in the Western Mediterranean region, testing two main predictions: (i) species with historical distributions in northern latitudes should have experienced greater loss of suitable habitat, resulting in higher extinction of historical lineages than species distributed in southern latitudes, where the effects of the ice ages were not as drastic. This would be reflected in the observation of lower diversity and number of differentiated lineages in northern areas. (ii) a signature of demographic expansion following the climate amelioration should be obvious in northern species, whereas in the south evidence of long-term effective population size stability should be observed. We used as models three species of wall lizards (Podarcis bocagei, Podarcis carbonelli and Podarcis vaucheri) that replace each other along the study area. We investigated the patterns of mitochondrial DNA diversity and subdivision and obtained demographic parameter estimates for each species. Our results suggest that P. bocagei, the northernmost species, bears low genetic diversity, a shallow coalescent history and marks of a demographic expansion. In contrast, P. vaucheri, the species with a southernmost distribution, shows deeper coalescence events, complex geographical substructure and no evidence for population growth. The species with an intermediate distribution, P. carbonelli, shows average levels of diversity, substructure and population growth. Taken together, these results conform to our main predictions and are explained by a differential influence of the ice ages on distinct latitudes.     

Contrasting microsatellite variation between subalpine and western larch, two closely related species with different distribution patterns

Subalpine larch (Larix lyallii Parl.) and western larch (Larix occidentalis Nutt.) represent two closely related species with contrasting abundance and distribution patterns in Western North America. Genetic diversity at seven informative microsatellite loci was determined for 19 populations of subalpine larch and nine populations of western larch. Contrasting genetic diversity and patterns of population differentiation were observed between the two species. The overall within-population genetic diversity parameters were lower in subalpine larch (A = 3.2; A(P) = 3.6; H(E) = 0.418) than in western larch (A(P) = 5.51; H(E) = 0.580), a pattern that is likely related to historical or demographic factors. No evidence of interspecific hybridization was observed. Significantly more population differentiation (theta = 0.15; R(ST) = 0.07), consistent with more restricted gene flow, was observed for subalpine larch as compared to western larch (theta = 0.05; R(ST) = 0.04). Under the assumption of an infinite allele mutation model, 12 of the 19 subalpine larch populations showed signs of deviation from the mutation-drift equilibrium, which suggests Holocene population bottlenecks and fluctuations in effective population size for this species. None of the western larch populations deviated significantly from the mutation-drift equilibrium. For both species, Mantel's test revealed a significant positive relationship between geographical and genetic distances indicative of isolation by distance. A similar geographical structure was detected in both species, suggesting at least two genetically distinct glacial populations in each species. The various implications for gene conservation are discussed.     

Contrasting patterns of population structure and demographic history in cryptic species of Bostrychia intricata (Rhodomelaceae,Rhodophyta) from New Zealand

Spatial patterns of genetic diversity provide insight into the demography and history of species. Morphologically similar but genetically distinct “cryptic” species are increasingly being recognized in marine organisms through molecular analyses. Such species are, on closer inspection, often discovered to display contrasting life histories or occasionally minor morphological differences; molecular tools can thus be useful indicators of diversity. Bostrychia intricata, a marine red alga, is widely distributed throughout the Southern Hemisphere and comprises many cryptic species. We used mitochondrial cytochrome c oxidase I gene sequences to assess the genetic variation, population genetic structure, and demographic history of B. intricata in New Zealand. Our results supported the existence of three cryptic species of B. intricata (N2, N4, and N5) in New Zealand. Cryptic species N4, which was found throughout New Zealand, showed a higher genetic diversity and wider distribution than the other two species, which were only found in the North Island and northern South Island. Our analyses showed low to moderate genetic differentiation among eastern North Island populations for cryptic species N2, but high differentiation among North and South Island populations for N4, suggesting different population structure between these cryptic species. Data also indicated that N2 has recently undergone population expansion, probably since the Last Glacial Maximum (LGM), while the higher genetic diversity in N4 populations suggests persistence in situ through the LGM. The contrasting population structures and inferred demographic histories of these species highlight that life history can vary greatly even among morphologically indistinguishable taxa.     

Estimates of gene flow, genetic substructure and population heterogeneity in bracken (Pteridium aquilinum)

Bracken [ Pteridium aquilinum (L.) Kuhn] is a cosmopolitan species and is a noxious weed in many areas. Because of its abundance, particularly in Britain, bracken affords an ideal system for investigating various aspects of population genetics and evolution. High mobility of dispersal units (spores) suggests that rates of gene flow among distant populations should be high. Gene flow is a major evolutionary force that influences the genetic structure of populations. To examine the effects of gene flow on population heterogeneity and population substructuring in bracken, starch gel electrophoresis of enzymes was used to provide the necessary genetic database. Allele frequency data at 21 loci were obtained for seven British populations, one Majorcan and one from the eastern United States. A model was employed to estimate the amount of gene flow ( Nm ) at several levels. Gene flow among British populations was extremely high ( Nm = 36.51), one of the highest estimates reported for plants. Among eight European populations gene flow was lower (but still considered high) at Nm = 2.47. Trans-Atlantic gene flow was low ( Nm = 0.0926).
F -statistics were used to assess population heterogeneity and substructuring. The data indicate that, compared with other species, there is very little genetic differentiation among British populations of bracken. Indeed, it appears that the whole island is behaving as a single randommating population. This result is consistent with high levels of gene flow. Only one population (on the Isle of Arran) showed statistically significant genetic substructuring. Habitat heterogeneity on the island and age structure are hypothesized as possible causes of this result.
The data reported here support previous studies demonstrating that bracken is genetically polymorphic and is an outcrossing species.     

Contrasting population structure and demographic history of cereal aphids in different environmental and agricultural landscapes

Genetic diversity of populations has important ecological and evolutionary consequences, whose understanding is fundamental to improve the sustainability of agricultural production. Studies of how differences in agricultural management and environment influence the population structure of insect pests are central to predict outbreaks and optimize control programs. Here, we have studied the population genetic diversity and evolution of Sitobion avenae and Sitobion miscanthi, which are among the most relevant aphid pests of cereals across Europe and Asia, respectively. We have used genotyping by sequencing (GBS) to identify genome‐wide single nucleotide polymorphisms (SNPs) to infer the geographic structure and migration patterns. In the present study, we show that the population structure in present‐day populations is different from that described in previous studies, which suggest that they have evolved recently possibly as a response to human‐induced changes in agriculture. This study shows that S. avenae in England is predominantly parthenogenetic and there has been a demographic and spatial expansion of a single genetic cluster, which could correspond with the insecticide resistance superclone identified in previous studies. Conversely, in China, S. miscanthi populations are mostly cyclical parthenogenetic, with one sexual stage in autumn to produce overwintering eggs, and there are six genetically differentiated subpopulations and high genetic differentiation between geographic locations, which suggests that further taxonomical research is needed. Unlike S. avenae in England, there is no evidence for insecticide resistance and there is no predominance of a single lineage in S. miscanthi in China.     

Patterns of morphological diversification of mainland Anolis lizards from northwestern South America

Anolis lizards are one of the most diverse vertebrate genera and are the classic example of adaptive radiation and convergent evolution. Anoles exhibit great morphological diversity produced by the ecological opportunity to exploit several arboreal niches. Anole radiation in the Caribbean islands is well studied, but the mainland radiation is less understood. We used a large morphological data set and a molecular phylogeny to describe the morphological diversification of anoles from northwestern South America, a region with the highest anole diversity on a mainland. We describe morphological diversity as summarized by ten morphotypes, defined mainly by body size, limb proportions, and subdigital lamellae. We show that some morphotypes are limited to forested lowlands and others to Andean highlands; by contrast, Anolis assemblages from tropical rainforests are comprised of the same four morphotypes. We demonstrate that morphological diversification followed a pattern of adaptive radiation across a landscape of adaptive peaks. Our results are consistent with the most recent hypothesis of convergence stated for Caribbean radiation, and demonstrate convergence between mainland morphotypes and Caribbean ecomorphs, which suggests that common processes are driving both radiations. © 2016 The Linnean Society of London     

Rapid evolution and gene-specific patterns of selection for three genes of spermatogenesis in Drosophila

Hybrid males resulting from crosses between closely related species of Drosophila are sterile. The F1 hybrid sterility phenotype is mainly due to defects occurring during late stages of development that relate to sperm individualization, and so genes controlling sperm development may have been subjected to selective diversification between species. It is also possible that genes of spermatogenesis experience selective constraints given their role in a developmental pathway. We analyzed the molecular evolution of three genes playing a role during the sperm developmental pathway in Drosophila at an early (bam), a mid (aly), and a late (dj) stage. The complete coding region of these genes was sequenced in different strains of Drosophila melanogaster and Drosophila simulans. All three genes showed rapid divergence between species, with larger numbers of nonsynonymous to synonymous differences between species than polymorphisms. Although this could be interpreted as evidence for positive selection at all three genes, formal tests of selection do not support such a conclusion. Departures from neutrality were detected only for dj and bam but not aly. The role played by selection is unique and determined by gene-specific characteristics rather than site of expression. In dj, the departure was due to a high proportion of neutral synonymous polymorphisms in D. simulans, and there was evidence of purifying selection maintaining a high lysine amino acid protein content that is characteristic of other DNA-binding proteins. The earliest spermatogenesis gene surveyed, which plays a role in both male and female gametogenesis, was bam, and its significant departure from neutrality was due to an excess of nonsynonymous substitutions between species. Bam is degraded at the end of mitosis, and rapid evolutionary changes among species might be a characteristic shared with other degradable transient proteins. However, the large number of nonsynonymous changes between D. melanogaster and D. simulans and a phylogenetic comparative analysis among species confirms evidence of positive selection driving the evolution of Bam and suggests an yet unknown germ cell line developmental adaptive change between these two species.     

Phylogeographical patterns shed light on evolutionary process in South America

The South American continent is composed of several biogeographical regions harbouring the highest biodiversity on the globe, encompassing five of the world's biodiversity ‘hot spots’. Nonetheless, the patterns and processes responsible for shaping its astonishing species diversity are largely unknown. Here, we present a review of current South American phylogeographical knowledge based on published articles on this topic. An appraisal of the literature reveals emerging phylogeographical patterns in the biota of South America. The striking phylogeographical divergence observed among organism lineages in South American studies is suggestive of high levels of undocumented species diversity. The interplay between Pleistocene climatic oscillations and Pliocene/Miocene orogenic events has contributed to shaping the current diversity and distribution of modern lineages in both the tropical and temperate regions of South America. Although older divergence times were observed for a range of species, most herpetofauna underwent an intraspecific lineage split much earlier than other organisms. The geographical ranges of species associated with forest habitats were reduced mainly during glacial cycles, whereas species associated with open vegetation domains have shown variable responses to climatic oscillations. The results suggest a highly complex mosaic of phylogeographical patterns in South America. We suggest future research directions to promote a better understanding of the origin and maintenance of the South American biota.     

Spatial patterns of caddisflies from Austral South America

Trichoptera are relevant members of the freshwater benthos that can exploit a broad variety of habitats. We explored their distribution patterns in Austral South America using a network quantitative approach applied on point occurrence data. The primary goal was to recognize groups of species connected by strong links of sympatry and secondarily to evaluate the structure of the sympatry network (patterns of connectivity) at higher taxonomic levels. We compiled 2,522 geographic points associated with 446 species. The strength of sympatry links were inferred directly from dot maps. Using the taxonomic resolution of species, we identified two major groups of co-occurring species that mimic in a geographic context the classic zoogeographic division of the study area: the Andean–Patagonian complex (of “cold-adapted” organisms) versus the Extra-Andean domain (of “warm-adapted” organisms) spanning the subtropical belt of the region. Within these major divisions, groups of highly co-distributed species were also found which in turn show a variety of overlapping spatial configurations. Overlap was centred at the following pivoting areas: Yungas of NW Argentina, Paranaense forest of NE Argentina and Valdivian temperate forest of Patagonia. Three very interesting findings emerged from the analyses: (i) Sierras Centrales of Córdoba correspond to the southern portion of many subtropical elements occurring in NW Argentina, (ii) some elements fit the disjunction between NW and NE Argentina and (iii) Uruguay is more closely related to NE Argentina than to Buenos Aires, suggesting that the Pampas region (Uruguay + Buenos Aires) could be an ill-defined biogeographic entity. Patterns at the species level become progressively blurred in going upward through the hierarchical classification.     

Contrasting abundance and contribution of clonal proliferation to the population structure of the corkscrew sea anemone Bartholomea annulata in the tropical Western Atlantic

Clonal propagation is an important life history trait for many sessile organisms, and often leads to the formation of monoclonal aggregations. In the marine environment, sea anemones have been model species for testing theory regarding the evolution of sex and understanding the contribution of sexual versus asexual reproduction to the population structure in facultatively clonal animals. However, little attention has been paid to tropical actiniarians. The corkscrew anemone Bartholomea annulata is common in tropical marine habitats in the western Atlantic and Caribbean; it forms small aggregations (2–4 anemones) on coral reefs and larger aggregations (>10 anemones) in mangrove habitats. We used field surveys and molecular analyses to investigate patterns of distribution, abundance, and genetic structure of aggregations formed by B. annulata on a reef in the US Virgin Islands and in a unique mangrove habitat in the Florida Keys. Abundance was greatest at the abandoned rock quarry mangrove habitat in the Florida Keys, where anemones formed continuously distributed aggregations carpeting the exposed limestone walls. Genetic diversity assessed via intersimple sequence repeats (ISSRs) and six microsatellite loci revealed that asexual reproduction plays only a minor role in the formation of both small and large anemone aggregations. Specifically, ISSR analyses showed that only ~10% of anemone aggregations were clonal in the US Virgin Islands, while microsatellite genotyping identified clonality in only 1 of 35 aggregations. In the Florida Keys, only four clonal genotypes were recovered within aggregations, but eight clones, representing 33% of the total surveyed population, had individuals in multiple aggregations. Thus, population structure of B. annulata appears to rely primarily on sexual reproduction, although asexual reproduction may play a nontrivial role in some environments. Mechanisms that drive the formation of genotypically diverse aggregations remain unknown, but may include attraction toward conspecifics, shared use of preferred habitats, or the local retention of larvae in partially enclosed habitats.     

AdmixKJump: identifying population structure in recently diverged groups

Motivation

Correctly modeling population structure is important for understanding recent evolution and for association studies in humans. While pre-existing knowledge of population history can be used to specify expected levels of subdivision, objective metrics to detect population structure are important and may even be preferable for identifying groups in some situations. One such metric for genomic scale data is implemented in the cross-validation procedure of the program ADMIXTURE, but it has not been evaluated on recently diverged and potentially cryptic levels of population structure. Here, I develop a new method, AdmixKJump, and test both metrics under this scenario.

Findings

I show that AdmixKJump is more sensitive to recent population divisions compared to the cross-validation metric using both realistic simulations, as well as 1000 Genomes Project European genomic data. With two populations of 50 individuals each, AdmixKJump is able to detect two populations with 100% accuracy that split at least 10KYA, whereas cross-validation obtains this 100% level at 14KYA. I also show that AdmixKJump is more accurate with fewer samples per population. Furthermore, in contrast to the cross-validation approach, AdmixKJump is able to detect the population split between the Finnish and Tuscan populations of the 1000 Genomes Project.

Conclusion

AdmixKJump has more power to detect the number of populations in a cohort of samples with smaller sample sizes and shorter divergence times.

Availability

A java implementation can be found at https://sites.google.com/site/igsevolgenomicslab/home/downloads     

Superoxide dismutase, catalase, glutathione peroxidase and gluthatione S-transferases M1 and T1 gene polymorphisms in three Brazilian population groups

Antioxidants such as superoxide dismutase (SOD), catalase (CAT) and glutathione peroxidase (GPX1) reduce the oxidation rates in the organism. Gluthatione S-transferases (GSTs) play a vital role in phase 2 of biotransformation of many substances. Variation in the expression of these enzymes suggests individual differences for the degree of antioxidant protection and geographical differences in the distribution of these variants. We described the distribution frequency of CAT (21A/T), SOD2 (Ala9Val), GPX1 (Pro198Leu), GSTM1 and GSTT1 polymorphisms in three Brazilian population groups: Kayabi Amerindians (n = 60), Kalunga Afro-descendants (n = 72), and an urban mixed population from Federal District (n = 162). Frequencies of the variants observed in Kalunga (18% to 58%) and Federal District (33% to 63%) were similar to those observed in Euro and Afro-descendants, while in Kayabi (3% to 68%), depending on the marker, frequencies were similar to the ones found in different ethnic groups. Except for SOD2 in all population groups studied here, and for GPX1 in Kalunga, the genotypic distributions were in accordance with Hardy-Weinberg Equilibrium. These data can clarify the contribution of different ethnicities in the formation of mixed populations, such as that of Brazil. Moreover, outcomes will be valuable resources for future functional studies and for genetic studies in specific populations. If these studies are designed to comprehensively explore the role of these genetic polymorphisms in the etiology of human diseases they may help to prevent inconsistent genotype-phenotype associations in pharmacogenetic studies.     

Application of maximum likelihood statistics for the population and family studies of inheritance of the chimpanzee R-C-E-F and V-A-B-D blood groups

Maximum likelihood statistics were applied to the analysis of serological data to confirm the originally proposed genetic models of the chimpanzee R-C-E-F and V-A-B-D systems. Five hundred ninety-nine chimpanzees, including 81 parents of 114 offspring, were tested for R-C-E-F, and 60 parents of 80 offspring were tested for V-A-B-D blood groups. An estimation-maximization procedure was used to obtain maximum likelihood estimates and support intervals of the haplotype frequencies. For each haplotype, the null hypothesis of nonexistence was evaluated. The frequencies obtained by this method do not differ significantly from those calculated by the square root formula, but put these estimates on a statistically more rigorous footing.     

Identification and population genetic comparison of three ascidian species based on mtDNA sequences

Ascidians are sessile marine chordate invertebrates found along seashores worldwide and are typically regarded as invasive organisms. Knowledge concerning their global genetic structure and subsequent invasive potential is limited. Here, we identified three ascidians—Ciona robusta, Ciona savignyi, and Styela clava from the northeast region of China using morphological characteristics and mitochondrial cytochrome c oxidase subunit I (cox1) as genetic marker. We additionally used phylogenetics to aid in the identification of these three species. The results of a population genetic analysis showed that among the three species, the level of haplotype diversity was particularly high within C. savignyi, and nucleotide diversity varied moderately. We divided the three species separately into native and invasive populations using 170 cox1 sequences from global resources to explore population genetic structure and invasive potential. Although in the network analysis Ciona spp. formed haplogroups of native and invasive populations, some haplotypes were still shared. We found that the haplotypes did not cluster within the network of S. clava. Our AMOVA results also showed that Ciona spp. had a weak genetic structure, and less genetic differentiation was present in S. clava. These data suggest that there are extensive incursions of these three ascidians into different geographical regions. Global comparisons of ascidian populations will help in the understanding of their population genetic structure and invasive potential, hence providing important insights regarding conservation as well as management.     

Numerical variation of the presacral vertebral column in three population groups in North America

A total of 1,239 skeletons from among Mongoloid, Caucasoid and Negroid population groups in North America was examined for variations in the number of presacral vertebrae. The overall incidence of variation was 11%; 6% with 23 and 5% with 25 presacral vertebrae. Differences in total variation among the three groups were not significant, but differences in the incidence of 23 and 25 presacral vertebrae among the groups were highly significant. Numerical vertebral variation occurred in 11% of both sexes, but with males having a higher frequency of 25 presacral vertebrae and females having a higher frequency of 23 presacral vertebrae. The incidence of the specific variation of 23 or 25 presacral vertebrae was not significantly different among the males of the three groups studied, whereas 23 presacral vertebrae were found significantly more often in the Negroid females. Numerical variation of vertebrae was not associated with age. The data strongly support the conclusion that the total frequency of variation in the number of presacral vertebrae is a specific characteristic of any particular population group and that there is a tendency in all population groups toward an increase in number in males and a decrease in number in females.     

Contrasting patterns of mitochondrial and microsatellite population structure in fragmented populations of greater prairie-chickens

Greater prairie-chickens (Tympanuchus cupido pinnatus) were once found throughout the tallgrass prairie of midwestern North America but over the last century these prairies have been lost or fragmented by human land use. As a consequence, many current populations of prairie-chickens have become isolated and small. This fragmentation of populations is expected to lead to reductions in genetic variation as a result of random genetic drift and a decrease in gene flow. As expected, we found that genetic variation at both microsatellite DNA and mitochondrial DNA (mtDNA) markers was reduced in smaller populations, particularly in Wisconsin. There was relatively little range-wide geographical structure (FST) when we examined mtDNA haplotypes but there was a significant positive relationship between genetic (FST) and geographical distance (isolation by distance). In contrast, microsatellite DNA loci revealed significant geographical structure (FST) and a weak effect of isolation by distance throughout the range. These patterns were much stronger when populations with reduced levels of genetic variability (Wisconsin) were removed from the analyses. This suggests that the effects of genetic drift were stronger than gene flow at microsatellite loci, whereas these forces were in range-wide equilibrium at mtDNA markers. These differences between the two molecular markers may be explained by a larger effective population size (Ne) for mtDNA, which is expected in species such as prairie-chickens that have female-biased dispersal and high levels of polygyny. Our results suggest that historic populations of prairie-chickens were once interconnected by gene flow but current populations are now isolated. Thus, maintaining gene flow may be important for the long-term persistence of prairie-chicken populations.     

Application of mitochondrial control region in population genetic studies of the shrimp Penaeus

This study reports a primer set for amplifying a partial fragment of about 610 bp in the fast mutating mitochondrial control region in shrimps of the genus Penaeus (Decapoda: Penaeidae). The utility of this amplified fragment for studying population differentiation and structuring, compared with more conservative mitochondrial genes (16S rRNA and COI), was explored in P. merguiensis populations over a vast geographical range based on sequence and RFLP analyses. The results indicate that the mitochondrial control region provides more informative sites and reveals more haplotypes, making it most useful for evaluating genetic variations within and between populations of Penaeus species.