genotype and genodive: two programs for the analysis of genetic diversity of asexual organisms

Limited intraspecific variation in drought resistance along a pronounced tropical rainfall gradient

Assessing within-species variation in response to drought is crucial for predicting species' responses to climate change and informing restoration and conservation efforts, yet experimental data are lacking for the vast majority of tropical tree species. We assessed intraspecific variation in response to water availability across a strong rainfall gradient for 16 tropical tree species using reciprocal transplant and common garden field experiments, along with measurements of gene flow and key functional traits linked to drought resistance. Although drought resistance varies widely among species in these forests, we found little evidence for within-species variation in drought resistance. For the majority of functional traits measured, we detected no significant intraspecific variation. The few traits that did vary significantly between drier and wetter origins of the same species all showed relationships opposite to expectations based on drought stress. Furthermore, seedlings of the same species originating from drier and wetter sites performed equally well under drought conditions in the common garden experiment and at the driest transplant site. However, contrary to expectation, wetter-origin seedlings survived better than drier-origin seedlings under wetter conditions in both the reciprocal transplant and common garden experiment, potentially due to lower insect herbivory. Our study provides the most comprehensive picture to date of intraspecific variation in tropical tree species' responses to water availability. Our findings suggest that while drought plays an important role in shaping species composition across moist tropical forests, its influence on within-species variation is limited.

Outbreeding reduces survival during metamorphosis in a headwater stream salamander

Assessing direct fitness effects of individual genetic diversity is challenging due to the intensive and long-term data needed to quantify survival and reproduction in the wild. But resolving these effects is necessary to determine how inbreeding and outbreeding influence eco-evolutionary processes. We used 8 years of capture-recapture data and single nucleotide polymorphism genotypes for 1906 individuals to test for effects of individual heterozygosity on stage-specific survival probabilities in the salamander Gyrinophilus porphyriticus. The life cycle of G. porphyriticus includes an aquatic larval stage followed by metamorphosis into a semi-aquatic adult stage. In our study populations, the larval stage lasts 6-10 years, metamorphosis takes several months, and lifespan can reach 20 years. Previous studies showed that metamorphosis is a sensitive life stage, leading us to predict that fitness effects of individual heterozygosity would occur during metamorphosis. Consistent with this prediction, monthly probability of survival during metamorphosis declined with multi-locus heterozygosity (MLH), from 0.38 at the lowest MLH (0.10) to 0.06 at the highest MLH (0.38), a reduction of 84%. Body condition of larvae also declined significantly with increasing MLH. These relationships were consistent in the three study streams. With evidence of localised inbreeding within streams, these results suggest that outbreeding disrupts adaptations in pre-metamorphic and metamorphic individuals to environmental gradients along streams, adding to evidence that headwater streams are hotspots of microgeographic adaptation. Our results also underscore the importance of incorporating life history in analyses of the fitness effects of individual genetic diversity and suggest that metamorphosis and similar discrete life stage transitions may be critical periods of viability selection.

Multiple paternity, fertilization success, and male quality: Mating system variation in the eelgrass, Zostera marina

Genetic diversity can modulate a population's response to a changing environment and plays a critical role in its ecological function. While multiple processes act to maintain genetic diversity, sexual reproduction remains the primary driving force. Eelgrass (Zostera marina) is an important habitat-forming species found in temperate coastal ecosystems across the globe. Recent increases in sea surface temperatures have resulted in shifts to a mixed-annual life-history strategy (i.e., displaying characteristics of both annual and perennial meadows) at its southern edge-of-range. Given that mating systems are intimately linked to standing levels of genetic variation, understanding the scope of sexual reproduction can illuminate the processes that shape genetic diversity. To characterize edge-of-range eelgrass mating systems, developing seeds on flowering Z. marina shoots were genotyped from three meadows in Topsail, North Carolina. In all meadows, levels of multiple mating were high, with shoots pollinated by an average of eight sires (range: 3-16). The number of fertilized seeds (i.e., reproductive success) varied significantly across sires (range: 1-25) and was positively correlated with both individual heterozygosity and self-fertilization. Outcrossing rates were high (approx. 70%) and varied across spathes. No clones were detected, and kinship among sampled flowering shoots was low, supporting observed patterns of reproductive output. Given the role that genetic diversity plays in enhancing resistance to and resilience from ecological disturbance, disentangling the links between life history, sexual reproduction, and genetic variation will aid in informing the management and conservation of this key foundation species.

Distinct geographic parthenogenesis in spite of niche conservatism and a single ploidy level: A case of Rubus ser. Glandulosi (Rosaceae)

Asexual organisms often differ in their geographic distributions from their sexual relatives. This phenomenon, termed geographic parthenogenesis, has long been known, but the underlying factors behind its diverse patterns have been under dispute. Particularly problematic is an association between asexuality and polyploidy in most taxa. Here, we present a new system of geographic parthenogenesis on the tetraploid level, promising new insights into this complex topic. We used flow cytometric seed screen and microsatellite genotyping to characterise the patterns of distribution of sexuals and apomicts and genotypic distributions in Rubus ser. Glandulosi across its range. Ecological modelling and local-scale vegetation and soil analyses were used to test for niche differentiation between the reproductive groups. Apomicts were detected only in North-western Europe, sexuals in the rest of the range in Europe and West Asia, with a sharp borderline stretched across Central Europe. Despite that, we found no significant differences in ecological niches. Genotypic richness distributions suggested independence of the reproductive groups and a secondary contact. We argue that unless a niche differentiation (resulting from polyploidy and/or hybridity) evolves, the main factors behind the patterns of geographic parthenogenesis in plants are phylogeographic history and neutral microevolutionary processes, such as clonal turnover.

Phylogenomics of American pika (Ochotona princeps) lineage diversification

Quaternary climate oscillations have profoundly influenced current species distributions. For many montane species, these fluctuations were a prominent driver in species range shifts, often resulting in intraspecific diversification, as has been the case for American pikas (Ochotona princeps). Range shifts and population declines in this thermally-sensitive lagomorph have been linked to historical and contemporary environmental changes across its western North American range, with previous research reconstructing five mitochondrial DNA lineages. Here, we paired genome-wide data (25,244 SNPs) with range-wide sampling to re-examine the number and distribution of intra-specific lineages, and investigate patterns of within- and among-lineage divergence and diversity. Our results provide genomic evidence of O. princeps monophyly, reconstructing six distinct lineages that underwent multiple rounds of divergence (0.809-2.81 mya), including a new Central Rocky Mountain lineage. We further found evidence for population differentiation across multiple spatial scales, and reconstructed levels of standing variation comparable to those found in other small mammals. Overall, our findings demonstrate the influence of past glacial cycles on O. princeps lineage diversification, suggest that current subspecific taxonomy may need to be revisited, and provide an important framework for investigations of American pika adaptive potential in the face of anthropogenic climate change.

Population genetics of Sida fallax Walp. (Malvaceae) in the Hawaiian Islands

Introduction

Sida fallax (Malvaceae) is the most widespread and variable taxon of Malvaceae in the Hawaiian Islands, growing with a diversity of morphological forms in different habitats including Midway Atoll, Nihoa, and all the main islands. Morphological variation exists within and among populations. The study aimed to investigate the genetic variation within and among populations from various habitats and geographic locations throughout the Hawaiian range of S. fallax.

Methods

A total of 124 samples, with up to five samples per population where possible, were collected from 26 populations across six of the main Hawaiian Islands (Kaua'i, O'ahu, Maui, Moloka'i, Lāna'i, and Hawai'i) and Nihoa in the Northwestern Hawaiian Islands. The sampling strategy encompassed collecting populations from different habitats and geographic locations, including coastal and mountain ecotypes, with many intermediate morphological forms. Multiplexed ISSR genotyping by sequencing (MIG-seq) was used to detect single nucleotide polymorphisms (SNP) and genetic differences among individuals and populations were evaluated using PCO analyses.

Results

The relationship of FST with the geographical distance between the populations was assessed using the Mantel test. The results showed that populations on a single island were more closely related to each other and to populations on islands within their respective groups than they were to populations on other islands.

Discussion

The overall genetic relationships among islands were, to a large extent, predictive based on island position within the chain and, to a lesser extent, within island topography.

Comparative analysis of molecular and morphological diversity in two diploid Paspalum species (Poaceae) with contrasting mating systems

KEY MESSAGE

Interspecific comparison of two Paspalum species has demonstrated that mating systems (selfing and outcrossing) contribute to variation (genetically and morphologically) within species through similar but mutually exclusive processes. Mating systems play a key role in the genetic dynamics of populations. Studies show that populations of selfing plants have less genetic diversity than outcrossing plants. Yet, many such studies have ignored morphological diversity. Here, we compared the morphological and molecular diversity patterns in populations of two phylogenetically-related sexual diploids that differ in their mating system: self-sterile Paspalum indecorum and self-fertile P. pumilum. We assessed the morphological variation using 16 morpho-phenological characters and the molecular diversity using three combinations of AFLPs. We compared the morphological and molecular diversity within and among populations in each mating system. Contrary to expectations, selfers showed higher morphological variation within populations, mainly in vegetative and phenological traits, compared to outcrossers. The high morphological variation within populations of selfers led to a low differentiation among populations. At molecular level, selfing populations showed lower levels of genotypic and genetic diversity than outcrossing populations. As expected, selfers showed higher population structure than outcrossers (PhiST = 0.301 and PhiST = 0.108, respectively). Increased homozygous combinations for the same trait/locus enhance morphological variation and reduce molecular variation within populations in selfing P. pumilum. Thus, selfing outcomes are opposite when comparing morphological and molecular variation in P. pumilum. Meanwhile, pollen flow in obligate outcrossing populations of P. indecorum increases within-population molecular variation, but tends to homogenize phenotypes within-population. Pollen flow in obligate outcrossers tends to merge geographically closer populations; but isolation by distance can lead to a weak differentiation among distant populations of P. indecorum.

Population genomic analyses reveal that salinity and geographic isolation drive diversification in a free-living protist

Protists make up the vast diversity of eukaryotic life and play a critical role in biogeochemical cycling and in food webs. Because of their small size, cryptic life cycles, and large population sizes, our understanding of speciation in these organisms is very limited. We performed population genomic analyses on 153 strains isolated from eight populations of the recently radiated dinoflagellate genus Apocalathium, to explore the drivers and mechanisms of speciation processes. Species of this genus inhabit both freshwater and saline habitats, lakes and seas, and are found in cold temperate environments across the world. RAD sequencing analyses revealed that the populations were overall highly differentiated, but morphological similarity was not congruent with genetic similarity. While geographic isolation was to some extent coupled to genetic distance, this pattern was not consistent. Instead, we found evidence that the environment, specifically salinity, is a major factor in driving ecological speciation in Apocalathium. While saline populations were unique in loci coupled to genes involved in osmoregulation, freshwater populations appear to lack these. Our study highlights that adaptation to freshwater through loss of osmoregulatory genes may be an important speciation mechanism in free-living aquatic protists.

Molecular genetics, seed morphology and fatty acids diversity in castor (Ricinus communis L., Euphorbiaceae) Iranian populations

BACKGROUND

Castor (Ricinus communis L.) seeds contain a large amount of oil that has several biological activities. In the current research, phytogeographic distribution, seed morphological characteristics, molecular genetic diversity and structure, and fatty acid composition were investigated in nine Iranian castor populations.

METHODS AND RESULTS

The cetyltrimethylammonium bromide (CTAB) protocol was used to extract the nuclear genomes. These were later amplified using 13 SCoT molecular primers. The phytogeographic distribution was determined based on the Zohary mapping, GC apparatus determined the fatty acid composition of the seeds. GenAlex, STRUCTURE, GenoDive, PopGene, and PopART software were used for the statistical analyzes. On phytogeographic mapping, the harvested populations belonged to different districts of the Euro-Siberian and Irano-Turanian regions (Holarctic kingdom). Most of the quantitative morphological traits of the seeds differed significantly (P ≤ 0.05) between the populations. The AMOVA test demonstrated a large proportion of significant genetic diversity assigned among populations, which were approved by some estimated parameters of genetic diversity such as Nm, Ht, Hs, and Gst. Nei's genetic distance and structure analysis confirmed the existence of two main genotype groups and some intermediates. However, there was no isolation by distance between the genotypes. Unsaturated fatty acids were detected as the main component of seed oil with linoleic and ricinoleic acids. Significant correlations were detected between the main fatty acids of seed oil with seed morphological traits, geographic distance and the geographic parameters of habitats. According to the composition of the seed fatty acids, four chemotypes groups were detected.

CONCLUSIONS

The classification patterns of the populations based on molecular genetic data, fatty acid composition, and phytogeographic mapping were not identical. These findings indicated that Iranian castor populations had unusual seed fatty acid composition which strongly depended on habitat geographic factors and seed morphological traits. However, the identified chemotypes and genotypes can be used in future breeding programs.

Habitat and climate influence hybridization among three genetically distinct Canada jay (Perisoreus canadensis) morphotypes in an avian hybrid zone complex

Examining the frequency and distribution of hybrids across contact zones provide insights into the factors mediating hybridization. In this study, we examined the effect of habitat and climate on hybridization patterns for three phenotypically, genetically, and ecologically distinct groups of the Canada jay (Perisoreus canadensis) in a secondary contact zone in western North America. Additionally, we tested whether the frequency of hybridization involving the three groups (referred to as Boreal, Pacific and Rocky Mountain morphotypes) is similar across the hybrid zones or whether some pairs have hybridized more frequently than others. We reanalyzed microsatellite, mtDNA and plumage data, and new microsatellite and plumage data for 526 individuals to identify putative genetic and phenotypic hybrids. The genetically and phenotypically distinct groups are associated with different habitats and occupy distinct climate niches across the contact zone. Most putative genetic hybrids (86%) had Rocky Mountain ancestry. Hybrids were observed most commonly in intermediate climate niches and in habitats where Engelmann spruce (Picea engelmannii) overlaps broadly with boreal and subalpine tree species. Our finding that hybrids occupy intermediate climate niches relative to parental morphotypes matches patterns for other plant and animal species found in this region. This study demonstrates how habitat and climate influence hybridization patterns in areas of secondary contact and adds to the growing body of research on tri-species hybrid zones.

Development of genetic markers for reproductive management of toucans

Retention of genetic diversity in successive generations is key to successful ex situ programs and will become increasingly important to restore wild populations of threatened animals. When animal genealogy is partly unknown or gaps exist in studbook records, the application of molecular resources facilitates informed breeding. Here, we apply molecular resources to an ex situ breeding population of toucans (Ramphastidae), a bird family zoos commonly maintain. Toucans face population declines from illegal poaching and habitat degradation. We developed novel microsatellite markers using blood samples from 15 Keel-billed Toucans (Ramphastos sulfuratus Lesson 1830). Parentage of two individuals was known a priori, but possible sibship among 13 putative founders-including the parents-was unknown. We compared available avian heterologous and novel microsatellite markers to recover known relationships and reconstruct sibship. Eight of 61 heterologous markers amplified consistently and were polymorphic, but less so than the 18 novel markers. Known sibship (and three sibling pairs whose relatedness was unknown a priori) and paternity-though not maternity except in one case-were well-recovered using both likelihood and pairwise relatedness methods, when incorporating novel but not heterologous markers. Zoo researchers seeking microsatellite primer sets for their breeding toucan populations will likely benefit from our heterologous markers, which can be leveraged both to assess relatedness and select breeding pairs. We recommend that zoo biologists rely on species-specific primers and not optimize heterologous primers for toucan species without molecular resources. We conclude with a brief discussion of modern genotyping methods of interest to zoo researchers.

Molecular phylogeography and historical demography of a widespread herbaceous species from eastern North America, Podophyllum peltatum

PREMISE

Glacial/interglacial cycles and topographic complexity are both considered to have shaped today's diverse phylogeographic patterns of taxa from unglaciated eastern North America (ENA). However, few studies have focused on the phylogeography and population dynamics of wide-ranging ENA herbaceous species occurring in forest understory habitat. We examined the phylogeographic pattern and evolutionary history of Podophyllum peltatum L., a widely distributed herb inhabiting deciduous forests of ENA.

METHODS

Using chloroplast DNA (cpDNA) sequences and nuclear microsatellite loci, we investigated the population structure and genetic diversity of the species. Molecular dating, demographic history analyses, and ecological niche modeling were also performed to illustrate the phylogeographic patterns.

RESULTS

Our cpDNA results identified three main groups that are largely congruent with boundaries along the Appalachian Mountains and the Mississippi River, two major geographic barriers in ENA. Populations located to the east of the Appalachians and along the central Appalachians exhibited relatively higher levels of genetic diversity. Extant lineages may have diverged during the late Miocene, and range expansions of different groups may have happened during the Pleistocene glacial/interglacial cycles.

CONCLUSIONS

Our findings indicate that geographic barriers may have started to facilitate the population divergence in P. peltatum before the Pleistocene. Persistence in multiple refugia, including areas around the central Appalachians during the Quaternary glacial period, and subsequent expansions under hospitable climatic condition, especially westward expansion, are likely responsible for the species' contemporary genetic structure and phylogeographic pattern.

Genetic Diversity of Aspergillus flavus on Maize in Guatemala

Aflatoxin contamination of maize is a leading threat to health in Guatemala. This contamination is the result of infection from Aspergillus flavus and has been effectively reduced in other countries through application of nonaflatoxigenic, indigenous strains of A. flavus. We collected 82 maize samples from throughout Guatemala in two years and isolated 272 A. flavus from these samples, including 126 unique genotypes. We provide here a phenotypic and simple sequence repeat (SSR)-based genotypic description of these isolates, as well as an analysis of the diversity of this population. High levels of genetic diversity were observed with the nonaflatoxigenic isolates in this study, but this information contributes to the development of indigenous aflatoxin biocontrol products.

Discrimination of Camellia cultivars using iD-NA analysis

Recently, many new cultivars have been taken abroad illegally, which is now considered an international issue. Botanical evidence found at a crime scene provides valuable information about the origin of the sample. However, botanical resources for forensic evidence remain underutilized because molecular markers, such as microsatellites, are not available without a limited set of species. Multiplexed intersimple sequence repeat (ISSR) genotyping by sequencing (MIG-seq) and its analysis method, identification of not applicable (iD-NA), have been used to determine several genome-wide genetic markers, making them applicable to all plant species, including those with limited available genetic information. Camellia cultivars are popular worldwide and are often planted in many gardens and bred to make new cultivars. In this study, we aimed to analyze Camellia cultivars/species through MIG-seq. MIG-seq could discriminate similar samples, such as bud mutants and closely related samples that could not be distinguished based on morphological features. This discrimination was consistent with that of a previous study that classified cultivars based on short tandem repeat (STR) markers, indicating that MIG-seq has the same or higher discrimination ability as STR markers. Furthermore, we observed unknown phylogenetic relationships. Because MIG-seq can be applied to unlimited species and low-quality DNA, it may be useful in various scientific fields.

Sex-Associated SNP Confirmation of Sex-Reversed Male Farmed Japanese Flounder Paralichthys olivaceus

Female Japanese flounder Paralichthys olivaceus grow more rapidly than the male. The goal of all-female commercial production requires an efficient method of genetic sex identification. We conducted genome-wide association analysis of female and male farmed Japanese flounder (n = 24 per phenotypic sex) and found all regions of chromosome 24 to be significantly associated with phenotypic sex, suggesting it as the sex chromosome. Genetic sex was identified based on single nucleotide polymorphisms (SNP) on chromosome 24 (n = 3568) using multidimensional scaling analysis, and individuals were clearly separated according to sex by the first dimension. The 61 SNPs most highly associated with sex were selected, and an amplicon-based SNP panel was developed. This was used to determine genetic sex of 39 females and 40 males. Eleven phenotypic males were assigned as female with XX genotype, suggesting sex reversal. Genetic sex was also assessed based on the indel of the amh gene promoter, which is the major candidate sex gene of Japanese flounder. We found four SNPs perfectly associated with genotypic sex in the sex-associated SNP panel, one of which was located in exon 2 of the amh gene. Along with the indel of the amh gene promoter, the sex-associated SNP panel will be of value in identifying genetic sex of farmed Japanese flounder. Molecular sexing will facilitate all-female production by breeding sex-reversed males.

Double reduction estimation and equilibrium tests in natural autopolyploid populations

Many bioinformatics pipelines include tests for equilibrium. Tests for diploids are well studied and widely available, but extending these approaches to autopolyploids is hampered by the presence of double reduction, the comigration of sister chromatid segments into the same gamete during meiosis. Though a hindrance for equilibrium tests, double reduction rates are quantities of interest in their own right, as they provide insights about the meiotic behavior of autopolyploid organisms. Here, we develop procedures to (i) test for equilibrium while accounting for double reduction, and (ii) estimate the double reduction rate given equilibrium. To do so, we take two approaches: a likelihood approach, and a novel U-statistic minimization approach that we show generalizes the classical equilibrium χ2 test in diploids. For small sample sizes and uncertain genotypes, we further develop a bootstrap procedure based on our U-statistic to test for equilibrium. We validate our methods on both simulated and real data.

Clones on the run: The genomics of a recently expanded partially clonal species

Why species that in their core areas mainly reproduce sexually become enriched with clones in marginal populations ("geographic parthenogenesis") remains unclear. Earlier hypotheses have emphasized that selection might promote clonality because it protects locally adapted genotypes. On the other hand, it also hampers recombination and adaptation to changing conditions. The aim of the present study was to investigate the early stages of range expansion in a partially clonal species and what drives an increase in cloning during such expansion. We used genome-wide sequencing to investigate the origin and evolution of large clones formed in a macroalgal species (Fucus vesiculosus) during a recent expansion into the postglacial Baltic Sea. We found low but persistent clonality in core populations, while at range margins, large dominant clonal lineages had evolved repeatedly from different sexual populations. A range expansion model showed that even when asexual recruitment is less favourable than sexual recruitment in core populations, repeated bottlenecks at the expansion front can establish a genetically eroded clonal wave that spreads ahead of a sexual wave into the new area. Genetic variation decreases by drift following repeated bottlenecks at the expansion front. This results in the emerging clones having low expected heterozygosity, which corroborated our empirical observations. We conclude that Baker's Law (clones being favoured by uniparental reproductive assurance in new areas) can play an important role during range expansion in partially clonal species, resulting in a complex spatiotemporal mosaic of clonal and sexual lineages that might persist during thousands of generations.

Genome-wide data suggest a revision in management of the Caspian Sea trout Salmo caspius

Salmo caspius Kessler, 1877 is one of the most commercially important species of Salmonidae in the southern basin of the Caspian Sea. The occurrence of its wild populations has undergone sever reduction during the last years. In spite of the yearly restocking activity, still no progress on the recovery of its wild population has been observed. Hence, the present study was done in order to assess the efficiency of the current restocking activity in the southern Caspian basin in term of genetic diversity. DNA extracts of 32 S. caspius from four different groups were screened using 62621 genome-wide single nucleotide polymorphisms (SNP). The overal genetic diversity and Fst values were 0.18 and 0.08, respectively. Considering the observed admixture pattern and the positive values for inbreeding coeficient it seems that S. caspius suffers from its small effective population size. In order to obtain the maximum performance, alonside with expanding the size of brood stocks, rehabilitation of the habitats and spawning rivers of this nationally endangered species is of great importance.

Genomic evidence indicates small island-resident populations and sex-biased behaviors of Hawaiian reef Manta Rays

BACKGROUND

Reef manta rays (Mobula alfredi) are globally distributed in tropical and subtropical seas. Their life history traits (slow growth, late maturity, low reproductive output) make them vulnerable to perturbations and therefore require informed management strategies. Previous studies have reported wide-spread genetic connectivity along continental shelves suggesting high gene flow along continuous habitats spanning hundreds of kilometers. However, in the Hawaiian Islands, tagging and photo-identification evidence suggest island populations are isolated despite proximity, a hypothesis that has not yet been evaluated with genetic data.

RESULTS

This island-resident hypothesis was tested by analyzing whole mitogenome haplotypes and 2048 nuclear single nucleotide polymorphisms (SNPs) between M. alfredi (n = 38) on Hawai'i Island and Maui Nui (the 4-island complex of Maui, Moloka'i, Lāna'i and Kaho'olawe). Strong divergence in the mitogenome (ΦST = 0.488) relative to nuclear genome-wide SNPs (neutral FST = 0.003; outlier FST = 0.186), and clustering of mitochondrial haplotypes among islands provides robust evidence that female reef manta rays are strongly philopatric and do not migrate between these two island groups. Combined with restricted male-mediated migration, equivalent to a single male moving between islands every 2.2 generations (~ 64 years), we provide evidence these populations are significantly demographically isolated. Estimates of contemporary effective population size (Ne) are 104 (95% CI: 99-110) in Hawai'i Island and 129 (95% CI: 122-136) in Maui Nui.

CONCLUSIONS

Concordant with evidence from photo identification and tagging studies, these genetic results indicate reef manta rays in Hawai'i have small, genetically-isolated resident island populations. We hypothesize that due to the Island Mass Effect, large islands provide sufficient resources to support resident populations, thereby making crossing deep channels separating island groups unnecessary. Small effective population size, low genetic diversity, and k-selected life history traits make these isolated populations vulnerable to region-specific anthropogenic threats, which include entanglement, boat strikes, and habitat degradation. The long-term persistence of reef manta rays in the Hawaiian Islands will require island-specific management strategies.

The Correlation between Genotype Richness of Submerged Macrophytes and Periphyton Biomass: A Mesocosm Study Based on Five Dominant Submerged Macrophytes from Yangtze River

Submerged macrophyte and periphyton are main primary producers which strongly interact with each other in clear water shallow lakes. In this study, the effects of genetic variation of the macrophyte species on periphyton biomass were studied in five submerged species. A two-year mesocosm study was conducted with four levels of genetic diversity (1, 4, 8 and 16 genotypes) for each submerged macrophyte, including 1600 individuals and 320 boxes in 20 mesocosms. Of the five submerged species, only Vallisneria spinulosa showed a positive correlation between its levels of genotype richness and the periphyton biomass. The correlation between genetic distance of genotypes and periphyton biomass was tested, which varied with the difference of seasons and species. In summary, we found that in freshwater mesocosms, the genetic diversity of submerged macrophytes may play a role in regulating the periphyton biomass, but the interaction between genetic diversity of macrophytes and periphyton biomass was not straightforward. This study will provide new insights into the interaction dynamics between the two primary producers in shallow lakes.

Relative Genetic Homogeneity within a Phenotypically Diverse group: The Case of Lake Tana Labeobarbus (Cyprinidae) Species Flock, Ethiopia

The Lake Tana Labeobarbus species flock represents one of the world's most famous examples of lacustrine species radiations. Previous studies of this group have resulted in the description of at least 15 species based on their differences in functional morphology and definition of two clades (lacustrine and riverine spawning clades) based on life history traits. A total of 166 fish representing 14 Labeobarbus species were genotyped using 10 lineage-specific hexaploid microsatellite loci. Six of these loci were developed for this study based on DNA sequence contigs derived from a microsatellite-enriched genomic library of Labeobarbus intermedius from Lake Tana; the remaining four loci were obtained from a previous study. The genotypes of the 10 loci were analyzed to examine genetic diversity and population structure within Lake Tana Labeobarbus. Overall mean allelic richness (NA) was 17.6 alleles per locus and observed (Ho) and expected (He) heterozygosities were 0.84 ± 0.14 and 0.73 ± 0.09, respectively, across all Lake Tana Labeobarbus samples examined. Our analyses reveal that there is little genetic differentiation among species (FST = 0.020-0.099; only 10 of 91 species comparisons were significant), but moderate differentiation (FST = 0.11, p < 0.05) between lacustrine and riverine spawning populations. Relative to previous phylogenetic hypotheses, our phenetic analysis employing the R-based Analysis of Phylogenetics and Evolution (APE) program seems to perform marginally better in revealing lineages within Lake Tana Labeobarbus. Herein, our results are compared to a previous microsatellite-based study of the same populations.

Geographic isolation and long-distance gene flow influence the genetic structure of the blue fan palm Brahea armata (Arecaceae)

The formation of the Baja California Peninsula (BCP) has impacted the microevolutionary dynamics of different species in ways that depend on biological traits such as dispersal capacity. Plants with relatively low levels of vagility have exhibited high genetic divergence between the BCP and Continental mainland. Brahea armata (Arecaceae) is a palm species inhabiting the northern part of the BCP and Sonora; its distribution occurs in isolated oases of vegetation. We aimed to evaluate the influence of the formation of the BCP on the genetic structure of B. armata using nuclear microsatellites and chloroplast markers (cpDNA) to compare patterns of genetic diversity and structure with previous published studies. Because gene flow through seeds is usually more limited compared to pollen flow, we expect to find stronger genetic structure at (cpDNA) than at nuclear markers. Moreover, larger genetic structure might also be explained by the smaller effective population size of cpDNA. We analyzed six microsatellite markers and two cpDNA regions. The main results indicated high levels of genetic differentiation among isolated populations located in the BCP, while low genetic differentiation was found between southern populations of the BCP and Sonora, suggesting long distance gene flow. In contrast, chloroplast markers indicated high levels of genetic structure between BCP and Sonora populations, suggesting asymmetrical gene flow between pollen (measured by nuclear microsatellites) and seed (cpDNA markers). This study provides valuable information on genetic diversity of B. armata that can be relevant for conservation and management; and develops microsatellites markers that can be transferred to other Brahea species.

Sponge diversification in marine lakes: Implications for phylogeography and population genomic studies on sponges

The relative influence of geography, currents, and environment on gene flow within sessile marine species remains an open question. Detecting subtle genetic differentiation at small scales is challenging in benthic populations due to large effective population sizes, general lack of resolution in genetic markers, and because barriers to dispersal often remain elusive. Marine lakes can circumvent confounding factors by providing discrete and replicated ecosystems. Using high-resolution double digest restriction-site-associated DNA sequencing (4826 Single Nucleotide Polymorphisms, SNPs), we genotyped populations of the sponge Suberites diversicolor (n = 125) to test the relative importance of spatial scales (1-1400 km), local environmental conditions, and permeability of seascape barriers in shaping population genomic structure. With the SNP dataset, we show strong intralineage population structure, even at scales <10 km (average F ST = 0.63), which was not detected previously using single markers. Most variation was explained by differentiation between populations (AMOVA: 48.8%) with signatures of population size declines and bottlenecks per lake. Although the populations were strongly structured, we did not detect significant effects of geographic distance, local environments, or degree of connection to the sea on population structure, suggesting mechanisms such as founder events with subsequent priority effects may be at play. We show that the inclusion of morphologically cryptic lineages that can be detected with the COI marker can reduce the obtained SNP set by around 90%. Future work on sponge genomics should confirm that only one lineage is included. Our results call for a reassessment of poorly dispersing benthic organisms that were previously assumed to be highly connected based on low-resolution markers.

Aflatoxin contamination of maize and groundnut in Burundi: Distribution of contamination, identification of causal agents and potential biocontrol genotypes of Aspergillus flavus

Aflatoxin contamination of the staples maize and groundnut is a concern for health and economic impacts across sub-Saharan Africa. The current study (i) determined aflatoxin levels in maize and groundnut collected at harvest in Burundi, (ii) characterized populations of Aspergillus section Flavi associated with the two crops, and (iii) assessed aflatoxin-producing potentials among the recovered fungi. A total of 120 groundnut and 380 maize samples were collected at harvest from eight and 16 provinces, respectively. Most of the groundnut (93%) and maize (87%) contained aflatoxin below the European Union threshold, 4 μg/kg. Morphological characterization of the recovered Aspergillus section Flavi fungi revealed that the L-morphotype of A. flavus was the predominant species. Aflatoxin production potentials of the L-morphotype isolates were evaluated in maize fermentations. Some isolates produced over 137,000 μg/kg aflatoxin B1. Thus, despite the relatively low aflatoxin levels at harvest, the association of both crops with highly toxigenic fungi poses significant risk of post-harvest aflatoxin contamination and suggests measures to mitigate aflatoxin contamination in Burundi should be developed. Over 55% of the L-morphotype A. flavus did not produce aflatoxins. These atoxigenic L-morphotype fungi were characterized using molecular markers. Several atoxigenic genotypes were detected across the country and could be used as biocontrol agents. The results from the current study hold promise for developing aflatoxin management strategies centered on biocontrol for use in Burundi to reduce aflatoxin contamination throughout the value chain.

Reproductive and Agronomic Characterization of Novel Apomictic Hybrids of Paspalum (Poaceae)

The tetraploid germplasm of Paspalum contains a large diversity that can be used to generate better forages. The objective was to evaluate a group of Paspalum notatum and Paspalum simplex apomictic hybrids for a set of agronomic traits and apomixis expressivity. Forage yield, cold tolerance, winter regrowth, and seed yield were evaluated. The expressivity of apomixis was evaluated in P. simplex hybrids by flow cytometry. Progeny testing with molecular markers was used to determine the genotypic variability in the progeny. Differences within P. notatum and P. simplex hybrids were observed for all traits, and some of them were superior in comparison with the controls. The accumulated forage yield during three years was 988 g m−2 in the P. notatum hybrids, whereas, in P. simplex, the average forage yield per harvest (40 days of regrowth) was 180 g m−2. In P. simplex, the apomixis expressivity varied between 0 and 100%, and 65% of the hybrids showed high apomixis expressivity (superior to 70%). The genotypic mean homogeneity in the progeny was 76% and 85% in P. notatum and P. simplex, respectively. The generation of hybrids with high apomixis expressivity that combine good agronomic performance and homogeneity in the offspring is possible in tetraploid P. notatum and P. simplex.

Genetic differentiation and genetic structure of mixed-ploidy Camellia hainanica populations

Camellia hainanica, which is common in China's Hainan Province, is an important woody olive tree species. Due to many years of geographic isolation, C. hainanica has not received the attention it deserves, which limits the exploitation of germplasm resources. Therefore, it is necessary to study population genetic characteristics for further utilization and conservation of C. hainanica. In this study, 96 individuals in six wild Camellia hainanica populations were used for ploidy analysis of the chromosome number, and the genetic diversity and population structure were investigated using 12 pairs of SSR primers. The results show complex ploidy differentiation in C. hainanica species. The ploidy of wild C. hainanica includes tetraploid, pentaploid, hexaploid, heptaploid, octoploid and decaploid species. Genetic analysis shows that genetic diversity and genetic differentiation among populations are low. Populations can be divided into two clusters based on their genetic structure, which matches their geographic location. Finally, to further maintain the genetic diversity of C. hainanica, ex-situ cultivation and in-situ management measures should be considered to protect it in the future.

A taxonomic revision of ten whitefish species from the lakes Lucerne, Sarnen, Sempach and Zug, Switzerland, with descriptions of seven new species (Teleostei, Coregonidae)

The taxonomy of the endemic whitefish of the lakes of the Reuss River system (Lucerne, Sarnen, Zug) and Lake Sempach, Switzerland, is reviewed and revised. Lake Lucerne harbours five species. Coregonus intermundiasp. nov. and C. suspensussp. nov., are described. Coregonus nobilis Haack, 1882, C. suidteri Fatio, 1885, and C. zugensis Nüsslin, 1882, are redescribed. Genetic studies have shown that C. suidteri and C. zugensis are composed of several distinct species endemic to different lakes. The names C. suidteri and C. zugensis are restricted to the species of lakes Sempach and Zug, respectively. The whitefish populations previously referred to as C. suidteri and C. zugensis from Lake Lucerne are described as C. litoralissp. nov. and C. muellerisp. nov., respectively. Furthermore, the whitefish from Lake Zug that were previously referred to as C. suidteri are described as C. supersumsp. nov. A holotype is designated for C. supersum that was previously one of two syntypes of C. zugensis. The other syntype is retained for C. zugensis. Coregonus obliterussp. nov. is described from Lake Zug, and C. obliterus and C. zugensis from Lake Zug are extinct. Finally, we describe C. sarnensissp. nov. from lakes Sarnen and Alpnach. Coregonus suidteri from Lake Sempach shows strong signals of introgression from deliberately translocated non-native whitefish species, which questions if the extant population still carries a genetic legacy from the original species and thus may need to be considered extinct. Coregonus suspensus is genetically partially of allochthonous origin, closely related to the radiation of Lake Constance. It is therefore compared to all known and described species of Lake Constance: C. wartmanni Bloch, 1784, C. macrophthalmus Nüsslin, 1882, C. arenicolus Kottelat,1997, and C. gutturosus Gmelin, 1818.

Wild snapdragon plant pedigree sheds light on limited connectivity enhanced by higher migrant reproductive success in a fragmented landscape

Background: In contrast with historical knowledge, a recent view posits that a non-negligible proportion of populations thrive in a fragmented landscape. One underlying mechanism is the maintenance of functional connectivity, i.e., the net flow of individuals or their genes moving among suitable habitat patches. Alternatively, functional connectivity might be typically limited but enhanced by a higher reproductive success of migrants. Methods: We tested for this hypothesis in wild snapdragon plants inhabiting six patches separated by seawater in a fragmented Mediterranean scrubland landscape. We reconstructed their pedigree by using a parentage assignment method based on microsatellite genetic markers. We then estimated functional connectivity and the reproductive success of plants resulting from between-patch dispersal events. Results: We found that wild snapdragon plants thrived in this fragmented landscape, although functional connectivity between habitat patches was low (i.e. 2.9%). The progeny resulting from between-patch dispersal events had a higher reproductive success than residents. Conclusion: Our findings imply that low functional connectivity in a fragmented landscapes may have been enhanced by higher reproductive success after migration. This original mechanisms might be partly compensating the negative impact of fragmentation.

Pollinator type strongly impacts gene flow within and among plant populations for six Neotropical species

Animal pollinators directly affect plant gene flow by transferring pollen grains between individuals. Pollinators with restricted mobility are predicted to limit gene flow within and among populations, whereas pollinators that fly longer distances are likely to promote genetic cohesion. These predictions, however, remain poorly tested. We examined population genetic structure and fine-scale spatial genetic structure (FSGS) in six perennial understory angiosperms in Andean cloud forests of northwestern Ecuador. Species belong to three families (Gesneriaceae, Melastomataceae, and Rubiaceae), and within each family we paired one insect-pollinated with one hummingbird-pollinated species, predicting that insect-pollinated species have greater population differentiation (as quantified with the F_ST statistic) and stronger FSGS (as quantified with the S_P statistic) than hummingbird-pollinated species. We confirmed putative pollinators through a literature review and fieldwork, and inferred population genetic parameters with a genome-wide genotyping approach. In two of the three species pairs, insect-pollinated species had much greater (>2-fold) population-level genetic differentiation and correspondingly steeper declines in fine-scale genetic relatedness. In the Gesneriaceae pair, however, F_ST and S_P values were similar between species and to those of the other hummingbird-pollinated plants. In this pair, the insect pollinators are euglossine bees (as opposed to small bees and flies in the other pairs), which are thought to forage over large areas, and therefore may provide similar levels of gene flow as hummingbirds. Overall, our results shed light on how different animal pollination modes influence the spatial scale of plant gene flow, suggesting that small insects strongly decrease genetic cohesion.

Pathometagenomics reveals susceptibility to intestinal infection by Morganella to be mediated by the blood group-related B4galnt2 gene in wild mice

Infectious disease is widely considered to be a major driver of evolution. A preponderance of signatures of balancing selection at blood group-related genes is thought to be driven by inherent trade-offs in susceptibility to disease. B4galnt2 is subject to long-term balancing selection in house mice, where two divergent allele classes direct alternative tissue-specific expression of a glycosyltransferase in the intestine versus blood vessels. The blood vessel allele class leads to prolonged bleeding times similar to von Willebrand disease in humans, yet has been maintained for millions of years. Based on in vivo functional studies in inbred lab strains, it is hypothesized that the cost of prolonged bleeding times may be offset by an evolutionary trade-off involving susceptibility to a yet unknown pathogen(s). To identify candidate pathogens for which resistance could be mediated by B4galnt2 genotype, we here employed a novel "pathometagenomic" approach in a wild mouse population, which combines bacterial 16S rRNA gene-based community profiling with histopathology of gut tissue. Through subsequent isolation, genome sequencing and controlled experiments in lab mice, we show that the presence of the blood vessel allele is associated with resistance to a newly identified subspecies of Morganella morganii, a clinically important opportunistic pathogen. Given the increasing importance of zoonotic events, the approach outlined here may find useful application in the detection of emerging diseases in wild animal populations.

Population Genomics Study and Implications for the Conservation of Zabelia tyaihyonii Based on Genotyping-By-Sequencing

Zabelia tyaihyonii (Nakai) Hisauti and H. Hara is a perennial shrub endemic to Republic of Korea that grows naturally in only a very limited region of the dolomite areas of Gangwon-do and Chungcheongbuk-do Provinces in the Republic of Korea. Given its geographical characteristics, it is more vulnerable than more widely distributed species. Despite the need for comprehensive information to support conservation, population genetic information for this species is very scarce. In this study, we analyzed the genetic diversity and population structure of 94 individuals from six populations of Z. tyaihyonii using a genotyping-by-sequencing (GBS) approach to provide important information for proper conservation and management. Our results, based on 3088 single nucleotide polymorphisms (SNPs), showed a mean expected heterozygosity (He) of 0.233, no sign of within-population inbreeding (G_IS that was close to or even below zero in all populations), and a high level of genetic differentiation (F_ST = 0.170). Analysis of molecular variance (AMOVA) indicated that the principal molecular variance existed within populations (84.5%) rather than among populations (17.0%). We suggested that six management units were proposed for conservation considering Bayesian structure analysis and phylogenetic analysis, and given the various current situations faced by Z. tyaihyonii, it is believed that not only the in situ conservation but also the ex situ conservation should be considered.

Application of genetic and Spatially Explicit Capture-Recapture analyses to design adaptive feral cat control in a large inhabited island

Faunas of oceanic islands have a high proportion of endemic species which contribute to the uniqueness of island communities. Island species are particularly naïve and vulnerable to alien predators, such as cats (Felis catus). On large, inhabited islands, where the complete eradication of feral cat populations is not considered feasible, control represents the best management option to lower their detrimental effects on native fauna. The first objective of our study was to investigate population genetics of feral cats of Réunion Island. The second objective was to understand the space use of feral cats established near the breeding colonies of the two endemic and endangered seabirds of Réunion Island, the Barau’s Petrel (Pterodroma baraui) and the Mascarene Petrel (Pseudobulweria aterrima). We evaluated genetic diversity, population structure and gene flow amongst six groups of feral cats located at a maximum of 10 km from known petrel colonies. We also analysed the behaviour and space use of one of these feral cat groups using camera-trap data and Spatially Explicit Capture-Recapture (SECR) models. Genetic analyses revealed that feral cats were structured in three genetic clusters explained mostly by the island topography. Two clusters were observed at five sampled sites, suggesting high connectivity amongst these sites. The last cluster was found in only one site, suggesting high isolation. This site was a remote mountain area located in the vicinity of one of the main Barau’s Petrel colonies. The behavioural study was conducted on this isolated feral cat population. Mark recapture analysis suggested that feral cats were present at low density and had large home ranges, which is probably explained by reduced food availability. Finally, we make several recommendations for refining feral cat management programmes on inhabited islands.

Possible northern persistence of Siebold's beech, Fagus crenata, at its northernmost distribution limit on an island in Japan Sea: Okushiri Island, Hokkaido

Siebold's beech, Fagus crenata, is widely distributed across the Japanese Archipelago and islands in Japan Sea. Similar to the northern limit of the geographical distribution of F. crenata on the mainland of Hokkaido, the northern limit of the distribution of F. crenata on islands in the Japan Sea is observed on Okushiri Island (ca 42°N). To understand the genetic relationships of F. crenata on Okushiri Island, we examined chloroplast (cp) DNA haplotypes and 11 nuclear microsatellite (SSR) loci among 1,838 individuals from 44 populations from Okushiri Island, mainland Hokkaido, and the northern part of the Tohoku region on Honshu Island. We identified 2 cpDNA haplotypes, which represent not only populations on the Japan Sea coast but also those on the Pacific coast and this suggested the Okushiri Island populations might not be formed by single colonization. Genetic diversity of the Okushiri Island populations of nuclear SSR was not lower than the mainland and the STRUCTURE analysis revealed the Okushiri Island individuals were admixed between Hokkaido and Tohoku clusters. Approximate Bayesian computation inferred that divergence between Tohoku and Hokkaido, and admixture between two populations which generated Okushiri populations occurred before the last glacial maximum (LGM), that is, 7,890 (95% hyper probability density (HPD): 3,420 - 9,910) and 3,870 (95% HPD: 431- 8,540) generations ago, respectively. These inferences were well supported by a geological history which suggested an isolation of Okushiri Island from Hokkaido started prior to the Middle Pleistocene. We discuss the possible persistence of F. crenata during the last glacial maximum on northern islands in the Japan Sea such as Okushiri Island.

vcfpop: Performing population genetics analyses for autopolyploids and aneuploids based on next-generation sequencing data sets

Polyploids are cells or organisms with a genome consisting of more than two sets of homologous chromosomes. Polyploid plants have important traits that facilitate speciation and are thus often model systems for evolutionary, molecular ecology and agricultural studies. However, due to their unusual mode of inheritance and double-reduction, diploid models of population genetic analysis cannot properly be applied to autopolyploids. To overcome this problem, we developed a software package entitled vcfpop to perform a variety of population genetic analyses for autopolyploids, such as parentage analysis, analysis of molecular variance, principal coordinates analysis, hierarchical clustering analysis and Bayesian clustering. We used three data sets to evaluate the capability of vcfpop to analyse large data sets on a desktop computer. This software is freely available at http://github.com/huangkang1987/vcfpop.

Seascape genomics of common dolphins (Delphinus delphis) reveals adaptive diversity linked to regional and local oceanography

Background

High levels of standing genomic variation in wide-ranging marine species may enhance prospects for their long-term persistence. Patterns of connectivity and adaptation in such species are often thought to be influenced by spatial factors, environmental heterogeneity, and oceanographic and geomorphological features. Population-level studies that analytically integrate genome-wide data with environmental information (i.e., seascape genomics) have the potential to inform the spatial distribution of adaptive diversity in wide-ranging marine species, such as many marine mammals. We assessed genotype-environment associations (GEAs) in 214 common dolphins (Delphinus delphis) along > 3000 km of the southern coast of Australia.

Results

We identified 747 candidate adaptive SNPs out of a filtered panel of 17,327 SNPs, and five putatively locally-adapted populations with high levels of standing genomic variation were disclosed along environmentally heterogeneous coasts. Current velocity, sea surface temperature, salinity, and primary productivity were the key environmental variables associated with genomic variation. These environmental variables are in turn related to three main oceanographic phenomena that are likely affecting the dispersal of common dolphins: (1) regional oceanographic circulation, (2) localised and seasonal upwellings, and (3) seasonal on-shelf circulation in protected coastal habitats. Signals of selection at exonic gene regions suggest that adaptive divergence is related to important metabolic traits.

Conclusion

To the best of our knowledge, this represents the first seascape genomics study for common dolphins (genus Delphinus). Information from the associations between populations and their environment can assist population management in forecasting the adaptive capacity of common dolphins to climate change and other anthropogenic impacts.

Genetic diversity, asexual reproduction and conservation of the edible fruit tree Spondias purpurea L. (Anacardiaceae) in the Costa Rican tropical dry forest

The term circa situm has been used to describe different conservation strategies within agricultural landscapes. Circa situm conserves planted or remnant species in farmlands, where natural vegetation has been modified through anthropogenic intervention. It has been proposed that trees planted or retained under circa situm conditions may contribute to maintaining genetic diversity, however information on the role of this strategy in preserving genetic diversity is scarce. The aim of this study was to determine the levels of genetic diversity and structure, and mating patterns in planted and unmanaged stands of the tropical fruit tree Spondias purpurea L. in north western Costa Rica. In three localities, we used seven polymorphic microsatellite loci and genotyped 201 adults and 648 seeds from planted and wild stands. We found no differences in genetic diversity among planted and wild stands. Genetic structure analysis revealed that gene flow occurs among planted and wild stands within localities. Clones were present and their diversity and evenness were both high and similar between planted and wild stands. The number of pollen donors per progeny array was low (Nep = 1.01) which resulted in high levels of correlated paternity (rp = 0.9). Asexual seeds were found in 4.6% of the progeny arrays, which had multilocus genotypes that were identical to the maternal trees. Our results show that although planted stands under circa situm conditions can maintain similar levels of genetic diversity than wild stands, the low number of sires and asexual seed formation could threaten the long term persistence of populations.

Origins of high latitude introductions of Aedes aegypti to Nebraska and Utah during 2019

Aedes aegypti (L.), the yellow fever mosquito, is also an important vector of dengue and Zika viruses, and an invasive species in North America. Aedes aegypti inhabits tropical and sub-tropical areas of the world and in North America is primarily distributed throughout the southern US states and Mexico. The northern range of Ae. aegypti is limited by cold winter months and establishment in these areas has been mostly unsuccessful. However, frequent introductions of Ae. aegypti to temperate, non-endemic areas during the warmer months can lead to seasonal activity and disease outbreaks. Two Ae. aegypti incursions were reported in the late summer of 2019 into York, Nebraska and Moab, Utah. These states had no history of established populations of this mosquito and no evidence of previous seasonal activity. We genotyped a subset of individuals from each location at 12 microsatellite loci and ~ 14,000 single nucleotide polymorphic markers to determine their genetic affinities to other populations worldwide and investigate their potential source of introduction. Our results support a single origin for each of the introductions from different sources. Aedes aegypti from Utah likely derived from Tucson, Arizona, or a nearby location. Nebraska specimen results were not as conclusive, but point to an origin from southcentral or southeastern US. In addition to an effective, efficient, and sustainable control of invasive mosquitoes, such as Ae. aegypti, identifying the potential routes of introduction will be key to prevent future incursions and assess their potential health threat based on the ability of the source population to transmit a particular virus and its insecticide resistance profile, which may complicate vector control.

Genetic structure of wild rice Zizania latifolia in an expansive heterogeneous landscape along a latitudinal gradient

Global aquatic habitats are undergoing rapid degradation and fragmentation as a result of climate change and changes in land use. Understanding the genetic variability and adaptive potential of aquatic plant species is thus important for conservation purposes. In this study, we investigated the genetic diversity and structure of the extant natural populations of Zizania latifolia from five river basins in China based on 46 microsatellite markers. We tested isolation by environment (IBE), isolation by resistance (IBR), and isolation by distance (IBD) patterns using a reciprocal causal model (RCM). Furthermore, we elucidated the impact of the environment on Z. latifolia genetic diversity using generalized linear models (GLMs) and spatially explicit mixed models. Low genetic diversity (HE = 0.125–0.433) and high genetic differentiation (FST = 0.641, Øpt = 0.654) were found. Higher historical gene flow (MH = 0.212–2.354) than contemporary gene flow (MC = 0.0112–0.0247) and significant bottlenecks in almost all populations were identified, highlighting the negative impact of wetland fragmentation. The IBE model was exclusively supported for all populations and in three river basins. The IBD and IBR models were supported in one river basin each. The maximum temperature of the warmest month and precipitation seasonality were the plausible environmental parameters responsible for the observed pattern of genetic diversity. Local adaptation signatures were found, with nine loci identified as outliers, four of which were gene-linked and associated with environmental variables. Based on these findings, IBE is more important than IBD and IBR in shaping the genetic structure of Z. latifolia.

Contrasting patterns of genetic diversity and differentiation across the continental disjunct range of a sexually polymorphic aquatic plant

BACKGROUND AND AIMS

Reproductive systems enabling opportunities for self-fertilization influence population genetic structure and play a key role in colonization and genetic differentiation during range expansion. Because of their well-developed powers of dispersal, aquatic plants often have widespread disjunct geographical distributions, providing opportunities to investigate the role of reproductive systems in structuring genetic variation between parts of the range that differ in migration history and ecology.

METHODS

We compared reproductive systems and spatial genetic structure of the freshwater aquatic macrophyte Sagittaria latifolia between disjunct western and eastern ranges of North America (NA). Populations of this species are most commonly either monoecious or dioecious. We examined chloroplast DNA haplotype variation to test the hypothesis that the western range of this species represents a secondary colonization from the east, and evaluated the roles of reproductive system differences and geography in structuring contemporary patterns of genetic variation at 11 polymorphic SSR (simple sequence repeat) loci.

KEY RESULTS

Chloroplast haplotyping revealed a single haplotype in western NA compared to numerous haplotypes in eastern NA, consistent with a genetic bottleneck during westward migration. Estimates of genetic diversity in eastern NA populations differed significantly between reproductive systems, but this pattern was not evident in the western range. Eastern populations could be reliably assigned to genetic clusters based on their reproductive systems, whereas western populations clustered primarily by geographical location.

CONCLUSIONS

The sparser distribution of aquatic habitats in the drier western range of S. latifolia, combined with secondary colonization of this region, probably cause the lower genetic diversity and increased differentiation among populations, possibly overriding the effects of reproductive system evident in the eastern portion of the range. Our findings demonstrate that the complex interplay between migratory history, reproductive systems and habitat availability plays an important role in structuring spatial patterns of genetic variation in disjunct plant populations.

Genetic population structure and reproductive system of two invasive Asian earthworms, Amynthas tokioensis and Amynthas agrestis

The invasive Asian earthworms, Amynthas tokioensis and A. agrestis, have been successful in entering North American forests in recent decades, with significant damage to both soils and above-ground environments. This success could be driven in part by a polyploid genetic system and parthenogenetic reproduction, often suggested as benefits for invasive species. Therefore, we assessed the genetic population structure, genetic diversity, and reproductive system of both species using morphological traits and panels of microsatellite markers. A total of 216 A. tokioensis and 196 A. agrestis from six sites in Vermont USA were analyzed. Although all worms were morphologically hermaphroditic, all the A. agrestis lacked the male pore (the structure allowing pass of sperm between individuals), and only 19% of the A. tokioensis possessed the male pore. All A. tokioensis earthworms were triploid (scored for three alleles for at least 1 locus, and usually several), and A. agrestis was a mix of triploid and diploid individuals. Notable was the high proportion (80%) of A. agrestis earthworms that were diploid at one site. There was clearly clonal reproduction, with identical seven- locus genotypes observed for earthworms from each site, with as many as 45 individuals with the identical genotype at one site. However, the earthworms were also genetically diverse, with 14 genotypes observed for A. tokioensis and 54 for A. agrestis, and with many singleton genotypes (a single individual). Most genotypes (71% for A. tokioensis and 92% for A. agrestis) were found at a single site. The greatest number of genotypes was found at a commercial nursery where fully 23/26 A. agrestis earthworms were singleton genotypes. As expected for the pattern of private clone alleles at sites, several measures of geographic genetic differentiation were positive, and as expected for triploid systems, an AMOVA analysis showed high within-individual genetic diversity. The paradox of clear clonal reproduction, but with a great number of genotypes for each species, and the mix of triploid and diploid individuals could be explained if the worms have been sexually reproductive, with the switch to the uniparental system only recently (or even if sexual reproduction is episodic). Last, a large number of microsatellite loci were recovered for each species and there sequence and suggested PCR primers are provided for free use by other researchers.

Elevated inbreeding in Heliconia tortuosa is determined by tropical forest stand age, isolation, and loss of hummingbird functional diversity

Forest conversion and habitat loss are major threats to biological diversity. Forest regeneration can mitigate the negative effects of old growth forest loss on species diversity, but less is known about the extent to which forest loss reduces genetic diversity in remnant populations and whether secondary forests play a role in the maintenance of genetic diversity. We quantified genetic diversity in a tropical hummingbird-pollinated understory herb, Heliconia tortuosa, across a landscape mosaic of primary and secondary forest regrowth. Using microsatellite genotypes from >850 adult and juvenile plants within 33 forest patches and extensive bird surveys, we examined the effect of contemporary and historical landscape features including forest age (primary vs. secondary forest), stand isolation, and pollinator assemblages on genetic diversity and levels of inbreeding in H. tortuosa. We found that inbreeding was up to 3x higher in secondary forest, and this effect was amplified with reductions in primary forest in the surrounding landscape through reduced observed heterozygosity in isolated fragments. Inbreeding in forest patches was negatively correlated with the local frequency of specialist long-distance foraging traplining hummingbirds. Traplining hummingbirds therefore appear to facilitate mating among unrelated plants - an inference we tested using empirically parameterized simulations. Higher levels of inbreeding in H. tortuosa are therefore associated with reduced functional diversity of hummingbirds in secondary forests and forest patches isolated from primary forests. Our findings suggest a cryptic consequence of primary forest loss and secondary forest regeneration through the disruption of mutualistic interactions resulting in the erosion of genetic diversity in a common understory plant.

Unique genetic structure of the human tapeworm Dibothriocephalus latus from the Alpine lakes region - a successful adaptation?

Dibothriocephalus latus is the most frequent causative agent of fish-borne zoonosis (diphyllobothriosis) in Europe, where it is currently circulating mainly in the Alpine lakes region (ALR) and Russia. Three mitochondrial genes (cox1, cob and nad3) and 6 microsatellite loci were analysed to determine how is the recently detected triploidy/parthenogenesis in tapeworms from ALR displayed at the DNA level. A geographically distant population from the Krasnoyarsk Reservoir in Russia (RU-KR) was analysed as a comparative population. One or 2 alleles of each microsatellite locus was detected in plerocercoids from RU-KR, corresponding to the microsatellite pattern of a diploid organism. In contrast, 1–3 alleles were observed in tapeworms from ALR, in accordance with their triploidy. The high diversity of mitochondrial haplotypes in D. latus from RU-KR implied an original and relatively stable population, but the identical structure of mitochondrial genes of tapeworms from ALR was probably a consequence of a bottleneck typical of introduced populations. These results indicated that the diploid/sexually reproducing population from RU-KR was ancestral, located within the centre of the distribution of the species, and the triploid/parthenogenetically reproducing subalpine population was at the margin of the distribution. The current study revealed the allelic structure of the microsatellite loci in the triploid tapeworm for the first time.

Recent speciation and phenotypic plasticity within a parthenogenetic lineage of Levantine whip spiders (Chelicerata: Amblypygi: Charinidae)

Caves constitute ideal study systems for investigating adaptation and speciation, as the abiotic conditions shared by aphotic habitats exert a set of environmental filters on their communities. Arachnids constitute an important component of many cave ecosystems worldwide. We investigated the population genomics of two whip spider species: Sarax ioanniticus, a widely distributed parthenogenetic species found across the eastern Mediterranean; and S. israelensis, a recently described troglomorphic species that is endemic to caves in Israel. Here, we show that S. israelensis is completely genetically distinct from S. ioanniticus and most likely also constitutes a parthenogen. Counterintuitively, despite the lack of genetic variability within S. ioanniticus and S. israelensis, we discovered considerable variation in the degree of median eye reduction, particularly in the latter species. Natural history data from captive-bred specimens of S. israelensis validated the interpretation of parthenogenesis. Our results are most consistent with a scenario of a sexual ancestral species that underwent speciation, followed by independent transitions to apomictic parthenogenesis in each of the two daughter species. Moreover, the lack of genetic variability suggests that variation in eye morphology in S. israelensis is driven exclusively by epigenetic mechanisms.

Genetic divergence and demography of pudu deer (Pudu puda) in five provinces of southern Chile, analyzed through latitudinal and longitudinal ranges

Pudu deer (Pudu puda) is endemic to the temperate rainforests of Chile. Genetic studies at different geographic scales for this species are required to better determine the genetic divergence within and among populations and their demography across the distribution range. These data can provide unique insights into the species or population status for conservation plans and decision-makers. We analyzed the mtDNA control region (CR) and cytochrome b (Cyt b) sequences of pudu deer in five provinces of southern Chile located at different latitudinal locations (Cautín, Valdivia, Osorno, Llanquihue and Chiloé Island) and three geographic areas within the studied provinces, representative of different longitudinal sites (Andes range, Central Valley and Coastal Range), to understand their genetic divergence and demography. The haplotype (H) and nucleotide (Π) diversities of CR and Cyt b ranged from 0.64286 to 0.98333 and from 0.00575 to 0.01022, respectively. CR diversity was significantly different among provinces, with Valdivia showing higher values than Llanquihue and Chiloé Island (H = 0.98333 vs. 0.64286–0.92727, P &lt; 0.05). Cyt b variation also showed significant differences among provinces, particularly, among Cautín and Llanquihue (H = 1.000 vs. 0.222, P &lt; 0.05). Genetic structuring among provinces was relatively high, as indicated by the FST index (FST = 0.41905). Clustering analysis indicated the presence of a distinctive cluster for Chiloé Island individuals. Fu’s FS and Tajima’s D based on CR revealed significant, negative deviations from equilibrium for Chiloé Island (D = -1.65898), Valdivia (Fs = -7.75335) and Llanquihue (Fs = -3.93267), suggesting population expansion in these provinces. Analysis at the longitudinal range showed significant differences among areas based on Π (P &lt; 0.05), with the Andes range and Central Valley showing higher diversity than the Coastal Range. Neither population structuring (FST = 0.01360, P &gt; 0.05) nor distinctive clusters in the longitudinal range were observed. Fu’s Fs and Tajima’s D were negative and significant for the Coastal Range based on CR (Fs = -6.64752, P &lt; 0.001) and Cyt b (D = -1.74110, P &lt; 0.05), suggesting the existence of population expansion. Our results suggest that pudu deer in the analyzed provinces is a genetically structured species, which could be associated with reduced panmixia among populations. The genetic divergence pattern and the population expansion recorded are likely to be associated with past processes of recolonization after Pleistocene glaciation events.