Genetic assignment methods for the direct, real-time estimation of migration rate: a simulation-based exploration of accuracy and power

Genetic substructure and admixture as important factors in linkage disequilibrium-based estimation of effective number of breeders in recovering wildlife populations

The number of effective breeders (Nb ) and effective population size (Ne ) are population parameters reflective of evolutionary potential, susceptibility to stochasticity, and viability. We have estimated these parameters using the linkage disequilibrium-based approach with LDNE through the latest phase of population recovery of the brown bears (Ursus arctos) in Finland (1993-2010; N = 621). This phase of the recovery was recently documented to be associated with major changes in genetic composition. In particular, differentiation between the northern and the southern genetic cluster declined rapidly within 1.5 generations. Based on this, we have studied effects of the changing genetic structure on Nb and Ne , by comparing estimates for whole Finland with the estimates for the two genetic clusters. We expected a potentially strong relationship between estimate sizes and genetic differentiation, which should disappear as the population recovers and clusters merge. Consistent with this, our estimates for whole Finland were lower than the sum of the estimates of the two genetic clusters and both approaches produced similar estimates in the end. Notably, we also found that admixed genotypes strongly increased the estimates. In all analyses, our estimates for Ne were larger than Nb and likely reflective for brown bears of the larger region of Finland and northwestern Russia. Conclusively, we find that neglecting genetic substructure may lead to a massive underestimation of Nb and Ne . Our results also suggest the need for further empirical analysis focusing on individuals with admixed genotypes and their potential high influence on Nb and Ne .

Genetic population structure of fishers (Pekania pennanti) in the Great Lakes region: remnants and reintroductions

Reintroduction programs have been pivotal in augmenting populations of fishers (Pekania pennanti (Erxleben, 1777)) and re-establishing them to their former range in North America. The majority of reintroduction efforts in fishers have been considered demographically successful, but reintroductions can alter genetic population structure and success has rarely been evaluated in fishers from a genetic standpoint. We used microsatellite data (n = 169) to examine genetic population structure of fishers in the Great Lakes region and comment on the success of past reintroductions at two different spatial scales. We found significant genetic population structure among source and reintroduced populations within the Great Lakes region and large-scale genetic structure between fisher populations located in two geographically distant regions (Great Lakes and Northeast) in the eastern United States. Reintroductions associated with the Great Lakes produced results that were largely consistent with other studies of fisher reintroductions in the Northeast. However, our data are the first to support a measurable impact on genetic population structure in Pekania pennanti pennanti (Erxleben, 1777) from a reintroduction using geographically distant source and reintroduced populations. When feasible, we strongly recommend that reintroduction programs include an investigation of the underlying genetic structure to better define intended goals and supplement measures of demographic success.

Host heterozygosity and genotype rarity affect viral dynamics in an avian subspecies complex

Genetic diversity at community, population and individual levels is thought to influence the spread of infectious disease. At the individual level, inbreeding and heterozygosity are associated with increased risk of infection and disease severity. Host genotype rarity may also reduce infection risk if pathogens are co-adapted to common or local hosts, but to date, no studies have investigated the relative importance of genotype rarity and heterozygosity for infection in a wild, sexually reproducing vertebrate. With beak and feather disease virus (BFDV) infection in a wild parrot (Platycercus elegans), we show that both heterozygosity and genotype rarity of individual hosts predicted infection, but in contrasting ways. Heterozygosity was negatively associated with probability of infection, but not with infection load. In contrast, increased host genotype rarity was associated with lower viral load in infected individuals, but did not predict infection probability. These effects were largely consistent across subspecies, but were not evident at the population level. Subspecies and age were also strongly associated with infection. Our study provides novel insights into infection dynamics by quantifying rarity and diversity simultaneously. We elucidate roles that host genetic diversity can play in infection dynamics, with implications for understanding population divergence, intraspecific diversity and conservation.

Conservation and genetic characterisation of common bean landraces from Cilento region (southern Italy): high differentiation in spite of low genetic diversity

Since its introduction from Central-South America to Italy almost 500 years ago, the common bean (Phaseolus vulgaris L.) was largely cultivated across the peninsula in hundreds of different landraces. However, globalisation and technological modernisation of agricultural practices in the last decades promoted the cultivation of few varieties at the expense of traditional and local agro-ecotypes, which have been confined to local markets or have completely disappeared. The aim of this study was to evaluate the genetic diversity and differentiation in 12 common bean landraces once largely cultivated in the Cilento region (Campania region, southern Italy), and now the object of a recovery program to save them from extinction. The analysis conducted using 13 nuclear microsatellite loci in 140 individuals revealed a high degree of homozygosity within each landrace and a strong genetic differentiation that was reflected in the success in assigning individuals to the source landrace. On the contrary, internal transcribed spacers 1 and 2, analysed in one individual per landrace, were highly similar among common bean landraces but allowed the identification of a cowpea variety (Vigna unguiculata Walp.), a crop largely cultivated in the Old World before the arrival of common bean from Americas. In conclusion, our study highlighted that conservation of landraces is important not only for the cultural and socio-economic value that they have for local communities, but also because the time and conditions in which they have been selected have led to that genetic distinctiveness that is at the basis of many potential agronomical applications and dietary benefits.

Population assignment in autopolyploids

Understanding the patterns of contemporary gene dispersal within and among populations is of critical importance to population genetics and in managing populations for conservation. In contrast to diploids, there are few studies of gene dispersal in autopolyploids, in part due to complex polysomic inheritance and genotype ambiguity. Here we develop a novel approach for population assignment for codominant markers for autotetraploids and autohexaploids. This method accounts for polysomic inheritance, unreduced gametes and unknown allele dosage. It can also utilise information regarding the origin and genotype of one parent for population assignment of maternal or paternal parents. Using simulations, we demonstrate that our approach achieves high levels of accuracy for assignment even when population divergence is low (F_ST~0.06) and with only 12 microsatellite loci. We also show that substantially higher accuracy is achieved when known maternal information is utilised, regardless of whether allele dosage is known. Although this novel method exhibited near identical levels of accuracy to Structure when population divergence was high, it performed substantially better for most parameters at moderate (F_ST=0.06) to low levels of divergence (F_ST=0.03). These methods fill an important gap in the toolset for autopolyploids and pave the way for investigating contemporary gene dispersal in a widespread group of organisms.

Classification and conservation priority of five Deccani sheep ecotypes of Maharashtra, India

Characterization of Indian livestock breeds has mostly been limited to single breed/population focused on either physical description of traditionally recognized breeds/populations or to their genetic description. Usually, morphological and genetic characterization has taken place in isolation. A parallel morphological characterization of genetically identified breeds or genetic characterization of morphologically described breeds is mostly missing, and their conservation priorities have largely been based on solely considering degree of endangerment. This study uses parallel approach based on morphometric and genetic differentiation for classification of five sheep ecotypes of Maharashtra state, and sets their conservation priority using threat parameters, current utilities/merits and contribution to genetic diversity. A total of 1101 animals were described for 7 body measurements for morphometric characterization. From this sample set, 456 animals were genotyped for 25 microsatellite markers for genetic characterization. Conservation priorities were assessed combining genetic and non-genetic factors. All studied traits varied significantly among ecotypes (p<0.05). All morphometric traits exhibited substantial sexual dimorphism except ear length. Males were 42% heavier than females. Madgyal sheep were the largest amongst the five ecotypes. In the stepwise discriminant analysis, all measured traits were significant and were found to have potential discriminatory power. Tail length was the most discriminatory trait. The Mahalanobis distance of the morphological traits between Kolhapuri and Madgyal was maximum (12.07) while the least differentiation was observed between Madgyal and Solapuri (1.50). Discriminant analysis showed that 68.12% sheep were classified into their source population. The Sangamneri sheep showed least assignment error (22%) whilst Solapuri exhibited maximum error level (41%). A total of 407 alleles were observed, with an average of 16.28 alleles per locus. Sufficient levels of genetic diversity were observed in all the ecotypes with observed heterozygosity values exceeding 0.47 and gene diversity values exceeding 0.76. About 6% of the total genetic variation was explained by population differences (F_ST = 0.059). Pairwise F_ST values indicated least differentiation between Solapuri and Madgyal (0.025). In terms of genetic distances, Kolhapuri and Lonand were most closely related (Ds = 0.177). The most probable structure clustering of the five studied populations was at K = 5. The study showed a fair congruence between the dendrogram constructed on the basis of Mahalanobis distances and Nei’s as well as Reynolds genetic distances. The findings gave highest conservation priority to Lonand and least to Solapuri ecotype.

Genetic diversity of bottlenose dolphin (Tursiops sp.) populations in the western North Pacific and the conservation implications

The evolutionary processes that shape patterns of diversity in highly mobile marine species are poorly understood, but important towards transferable inference on their effective conservation. In this study, bottlenose dolphins (Tursiops sp.) are studied to address this broader question. They exhibit remarkable geographical variation for morphology, life history, and genetic diversity, and this high level of variation has made the taxonomy of the genus controversial. A significant population structure has been reported for the most widely distributed species, the common bottlenose dolphin (T. truncatus), in almost all ocean basins, though no data have been available for the western North Pacific Ocean (WNP). The genetic diversity of bottlenose dolphins in the WNP was investigated based on 20 microsatellite and one mitochondrial DNA markers for samples collected from Taiwanese, Japanese, and Philippine waters (9°-39°N, 120°-140°E) during 1986-2012. The results indicated that there are at least four genetically differentiated populations of common bottlenose dolphins in the western and central North Pacific Ocean. The pattern of differentiation appears to correspond to habitat types, resembling results seen in other populations of the same species. Our analyses also showed that there was no evident gene flow between the two "sister species", the common bottlenose dolphins, and the Indo-Pacific bottlenose dolphins (T. aduncus) occurring sympatrically in our study region.

Population genetics, community of parasites, and resistance to rodenticides in an urban brown rat (Rattus norvegicus) population

Brown rats are one of the most widespread urban species worldwide. Despite the nuisances they induce and their potential role as a zoonotic reservoir, knowledge on urban rat populations remains scarce. The main purpose of this study was to characterize an urban brown rat population from Chanteraines park (Hauts-de-Seine, France), with regards to haematology, population genetics, immunogenic diversity, resistance to anticoagulant rodenticides, and community of parasites. Haematological parameters were measured. Population genetics was investigated using 13 unlinked microsatellite loci. Immunogenic diversity was assessed for Mhc-Drb. Frequency of the Y139F mutation (conferring resistance to rodenticides) and two linked microsatellites were studied, concurrently with the presence of anticoagulant residues in the liver. Combination of microscopy and molecular methods were used to investigate the occurrence of 25 parasites. Statistical approaches were used to explore multiple parasite relationships and model parasite occurrence. Eighty-six rats were caught. The first haematological data for a wild urban R. norvegicus population was reported. Genetic results suggested high genetic diversity and connectivity between Chanteraines rats and surrounding population(s). We found a high prevalence (55.8%) of the mutation Y139F and presence of rodenticide residues in 47.7% of the sampled individuals. The parasite species richness was high (16). Seven potential zoonotic pathogens were identified, together with a surprisingly high diversity of Leptospira species (4). Chanteraines rat population is not closed, allowing gene flow and making eradication programs challenging, particularly because rodenticide resistance is highly prevalent. Parasitological results showed that co-infection is more a rule than an exception. Furthermore, the presence of several potential zoonotic pathogens, of which four Leptospira species, in this urban rat population raised its role in the maintenance and spread of these pathogens. Our findings should stimulate future discussions about the development of a long-term rat-control management program in Chanteraines urban park.

City-Dwellers and Country Folks: Lack of Population Differentiation Along an Urban-Rural Gradient in the Mosquito Culex pipiens (Diptera: Culicidae)

Mosquitoes (Diptera, Culicidae) occur in natural, urban, and peri-urban areas throughout the globe. Although the characteristics of urban and peri-urban habitats differ from those of natural habitats in many ways (e.g., fragmentation, pollution, noise, and light), few studies have examined the population connectivity of mosquitoes in urban areas. To obtain an overview of the species composition, we sampled mosquitoes from 23 sites in and around the city of Berlin, Germany. Of 23 species, five occurred in urban, 10 in peri-urban, and 20 in rural areas. Culex pipiens Linnaeus (Diptera: Culicidae) was the most common species collected (75% of all individuals) and occurred in all habitats. Hence this species was selected to be analysed at 10 microsatellite markers. There was no significant differentiation (FST = 0.016, P = 0.9) or isolation by distance (P = 0.06) among Cx. pipiens populations along an urban-rural gradient. The only significant differences detected were between Cx. pipiens and a laboratory population of Cx. pipiens f. molestus (pairwise FST = 0.114-0.148, P ≤ 0.001 in all comparisons), suggesting that the markers chosen were suitable for the identification of population differentiation. Our results indicate that Cx. pipiens gene flow is widespread within and among urban, peri-urban, and rural areas and that urban habitat does not necessarily impede or enhance gene flow among these populations.

Using genetic tools to estimate the prevalence of non-native red deer (Cervus elaphus) in a Western European population

Game species like the red deer have been subjected to anthropogenic impacts for centuries. Translocations are often carried out—sometimes illegally—not only for sporting purposes, but also to increase trophy quality, reduce inbreeding, or mitigate bottlenecks after excessive persecution. Apart from the blurring of large‐scale genetic structure, translocations without adequate quarantine measure risk introducing pathogens into potentially immunologically naïve populations. It is therefore important to understand the frequency of clandestine translocations. Identification of non‐autochthonous animals and their potential origin is often difficult and, in red deer, has been hampered by the lack of large‐scale genotypic datasets for comparison. In the present study, we make use of a recently published European‐wide microsatellite dataset to detect and quantify the presence of non‐autochthonous red deer in a large population sample (n = 1,780) from Central Europe (Belgium). Using factorial correspondence analysis, assignment tests and Bayesian clustering algorithms we arrive at an estimate of 3.7% non‐autochthonous animals (or their descendants). Some of these animals were assigned to a nearby French population and may have immigrated into Belgium naturally, but the large majority must have been introduced by humans. Our analysis pointed to the British Isles and Germany/Poland as the potential origin of many introduced deer, regions known to have been source populations for translocations in Europe and beyond. We found evidence for recreational hunters using carcasses from farmed deer to fulfill mandatory hunting quotas. Our study is the first to quantify the extent of human‐mediated introductions in a European game species at such a large scale with large and representative sample sizes.

Allometric trajectories of body and head morphology in three sympatric Arctic charr (Salvelinus alpinus (L.)) morphs

A study of body and head development in three sympatric reproductively isolated Arctic charr (Salvelinus alpinus (L.)) morphs from a subarctic lake (Skogsfjordvatn, northern Norway) revealed allometric trajectories that resulted in morphological differences. The three morphs were ecologically assigned to a littoral omnivore, a profundal benthivore and a profundal piscivore, and this was confirmed by genetic analyses (microsatellites). Principal component analysis was used to identify the variables responsible for most of the morphological variation of the body and head shape. The littoral omnivore and the profundal piscivore morph had convergent allometric trajectories for the most important head shape variables, developing bigger mouths and relatively smaller eyes with increasing head size. The two profundal morphs shared common trajectories for the variables explaining most of the body and head shape variation, namely head size relative to body size, placement of the dorsal and pelvic fins, eye size and mouth size. In contrast, the littoral omnivore and the profundal benthivore morphs were not on common allometric trajectories for any of the examined variables. The findings suggest that different selective pressures could have been working on traits related to their trophic niche such as habitat and diet utilization of the three morphs, with the two profundal morphs experiencing almost identical environmental conditions.

Mitochondrial and nuclear DNA reveals reticulate evolution in hares (Lepus spp., Lagomorpha, Mammalia) from Ethiopia

For hares (Lepus spp., Leporidae, Lagomorpha, Mammalia) from Ethiopia no conclusive molecular phylogenetic data are available. To provide a first molecular phylogenetic model for the Abyssinian Hare (Lepus habessinicus), the Ethiopian Hare (L. fagani), and the Ethiopian Highland Hare (L. starcki) and their evolutionary relationships to hares from Africa, Eurasia, and North America, we phylogenetically analysed mitochondrial ATPase subunit 6 (ATP6; n = 153 / 416bp) and nuclear transferrin (TF; n = 155 / 434bp) sequences of phenotypically determined individuals. For the hares from Ethiopia, genotype composition at twelve microsatellite loci (n = 107) was used to explore both interspecific gene pool separation and levels of current hybridization, as has been observed in some other Lepus species. For phylogenetic analyses ATP6 and TF sequences of Lepus species from South and North Africa (L. capensis, L. saxatilis), the Anatolian peninsula and Europe (L. europaeus, L. timidus) were also produced and additional TF sequences of 18 Lepus species retrieved from GenBank were included as well. Median joining networks, neighbour joining, maximum likelihood analyses, as well as Bayesian inference resulted in similar models of evolution of the three species from Ethiopia for the ATP6 and TF sequences, respectively. The Ethiopian species are, however, not monophyletic, with signatures of contemporary uni- and bidirectional mitochondrial introgression and/ or shared ancestral polymorphism. Lepus habessinicus carries mtDNA distinct from South African L. capensis and North African L. capensis sensu lato; that finding is not in line with earlier suggestions of its conspecificity with L. capensis. Lepus starcki has mtDNA distinct from L. capensis and L. europaeus, which is not in line with earlier suggestions to include it either in L. capensis or L. europaeus. Lepus fagani shares mitochondrial haplotypes with the other two species from Ethiopia, despite its distinct phenotypic and microsatellite differences; moreover, it is not represented by a species-specific mitochondrial haplogroup, suggesting considerable mitochondrial capture by the other species from Ethiopia or species from other parts of Africa. Both mitochondrial and nuclear sequences indicate close phylogenetic relationships among all three Lepus species from Ethiopia, with L. fagani being surprisingly tightly connected to L. habessinicus. TF sequences suggest close evolutionary relationships between the three Ethiopian species and Cape hares from South and North Africa; they further suggest that hares from Ethiopia hold a position ancestral to many Eurasian and North American species.

Population genetics of the Schistosoma snail host Bulinus truncatus in Egypt

The tropical freshwater snail Bulinus truncatus serves as an important intermediate host of several human and cattle Schistosoma species in many African regions. Despite some ecological and malacological studies, there is no information on the genetic diversity of B. truncatus in Egypt. Here, we sampled 70-100 snails in ten localities in Upper Egypt and the Nile Delta. Per locality, we sequenced 10 snails at a partial fragment of the cytochrome c oxidase subunit 1 gene (cox1) and we genotyped 25-30 snails at six microsatellite markers. A total of nine mitochondrial haplotypes were detected, of which five were unique to the Nile Delta and three were unique to Upper Egypt, indicating that snail populations may have evolved independently in both regions. Bayesian clustering and hierarchical F-statistics using microsatellite markers further revealed strong population genetic structure at the level of locality. Observed heterozygosity was much lower compared to what is expected under random mating, which could be explained by high selfing rates, population size reductions and to a lesser extent by the Wahlund effect. Despite these observations, we found signatures of gene flow and cross-fertilization, even between snails from the Nile Delta and Upper Egypt, indicating that B. truncatus can travel across large distances in Egypt. These observations could have serious consequences for disease epidemiology, as it means that infected snails from one region could rapidly and unexpectedly spark a new epidemic in another distant region. This could be one of the factors explaining the rebound of human Schistosoma infections in the Nile Delta, despite decades of sustained schistosomiasis control.

Spatial connectivity in an adult-sedentary reef fish with extended pelagic larval phase

Understanding the spatial scale of demographic connectivity in marine reef fishes dispersing pelagic larvae is a challenging task because of the technical difficulties associated with tagging and monitoring the movements of progeny at early life stages. Several studies highlighted a strong importance of local retention with levels of dispersal of ecological significance restricted to short distances. To date little information is available in species where pelagic dispersal lasts for long periods of time. In this work, population structure and connectivity were studied in the grey triggerfish, Balistes capriscus. Grey triggerfish larvae and juveniles remain associated with floating Sargassum sp. beds for an estimated period of 4-7 months before settling on benthic habitats where they remain sedentary as adults. Analysis of genetic variation among populations along the continental shelf of the northern Gulf of Mexico and U.S. east coast, encompassing over 3,100 km of coastline, revealed homogeneous allele frequencies and a weak isolation-by-distance pattern. Moment and maximum-likelihood estimates of dispersal parameters both indicated occurrence of large neighbourhoods with estimates of the dispersal distribution parameter σ of 914 and 780 km, respectively. Simulated distributions of dispersal distances using several distribution functions all featured substantial fractions of long-distance dispersal events with the 90% percentiles of travel distance prior to settlement averaging 1,809 km. These results suggest a high dependency of local recruitment on the output of nonlocal spawning stocks located hundreds of kilometres away and a reduced role of local retention in this species.

Tracking the return of Aedes aegypti to Brazil, the major vector of the dengue, chikungunya and Zika viruses

BACKGROUND

Aedes aegypti, commonly known as "the yellow fever mosquito", is of great medical concern today primarily as the major vector of dengue, chikungunya and Zika viruses, although yellow fever remains a serious health concern in some regions. The history of Ae. aegypti in Brazil is of particular interest because the country was subjected to a well-documented eradication program during 1940s-1950s. After cessation of the campaign, the mosquito quickly re-established in the early 1970s with several dengue outbreaks reported during the last 30 years. Brazil can be considered the country suffering the most from the yellow fever mosquito, given the high number of dengue, chikungunya and Zika cases reported in the country, after having once been declared "free of Ae. aegypti".

METHODOLOGY/PRINCIPAL FINDINGS

We used 12 microsatellite markers to infer the genetic structure of Brazilian Ae. aegypti populations, genetic variability, genetic affinities with neighboring geographic areas, and the timing of their arrival and spread. This enabled us to reconstruct their recent history and evaluate whether the reappearance in Brazil was the result of re-invasion from neighboring non-eradicated areas or re-emergence from local refugia surviving the eradication program. Our results indicate a genetic break separating the northern and southern Brazilian Ae. aegypti populations, with further genetic differentiation within each cluster, especially in southern Brazil.

CONCLUSIONS/SIGNIFICANCE

Based on our results, re-invasions from non-eradicated regions are the most likely scenario for the reappearance of Ae. aegypti in Brazil. While populations in the northern cluster are likely to have descended from Venezuela populations as early as the 1970s, southern populations seem to have derived more recently from northern Brazilian areas. Possible entry points are also revealed within both southern and northern clusters that could inform strategies to control and monitor this important arbovirus vector.

Decline of heterozygosity in a large but isolated population: a 45-year examination of moose genetic diversity on Isle Royale

Wildlife conservation and management approaches typically focus on demographic measurements to assess population viability over both short and long periods. However, genetic diversity is an important predictor of long term population vitality. We investigated the pattern of change in genetic diversity in a large and likely isolated moose (Alces alces) population on Isle Royale (Lake Superior) from 1960–2005. We characterized samples, partitioned into five different 5-year periods, using nine microsatellite loci and a portion of the mtDNA control region. We also simulated the moose population to generate a theoretical backdrop of genetic diversity change. In the empirical data, we found that the number of alleles was consistently low and that observed heterozygosity notably declined from 1960 to 2005 (p = 0.08, R² = 0.70). Furthermore, inbreeding coefficients approximately doubled from 0.08 in 1960–65 to 0.16 in 2000–05. Finally, we found that the empirical rate of observed heterozygosity decline was faster than the rate of observed heterozygosity loss in our simulations. Combined, these data suggest that genetic drift and inbreeding occurred in the Isle Royale moose populations over the study period, leading to significant losses in heterozygosity. Although inbreeding can be mitigated by migration, we found no evidence to support the occurrence of recent migrants into the population using analysis of our mtDNA haplotypes nor microsatellite data. Therefore, the Isle Royale moose population illustrates that even large populations are subjected to inbreeding in the absence of migration.

Investigating gene flow between the blind cavefish Garra barreimiae and its conspecific surface populations

Cave-dwelling taxa often share the same phenotypic modifications like absence of eyes and pigmentation. These “troglomorphic characters” are expressed in the populations of Garra barreimiae from the Al Hoota Cave and nearby Hoti Pit in Northern Oman. Surface morphotypes of this cyprinid species are common throughout the distribution area. Very rarely individuals with intermediate phenotypes can be found. In the present study, potential gene flow between cave and surface populations was tested and population structure within five sampling sites was assessed. Overall, 213 individuals were genotyped at 18 microsatellite loci. We found that the cave populations have lower genetic diversity and are clearly isolated from the surface populations, which seem to be sporadically in contact with each other. The results indicate a recent genetic bottleneck in the cave populations. Thus, it can be assumed that during climatic changes the connection between cave and surface water bodies was disjoined, leaving a subpopulation trapped inside. Nevertheless, occasional gene flow between the morphotypes is detectable, but hybridisation seems only possible in cave habitat with permanent connection to surface water. Individuals from surface sites bearing intermediate phenotypes but cave genotypes imply that phenotypic plasticity might play a role in the development of the phenotype.

Hybridization between two cryptic filamentous brown seaweeds along the shore: analysing pre- and postzygotic barriers in populations of individuals with varying ploidy levels

We aimed to study the importance of hybridization between two cryptic species of the genus Ectocarpus, a group of filamentous algae with haploid-diploid life cycles that include the principal genetic model organism for the brown algae. In haploid-diploid species, the genetic structure of the two phases of the life cycle can be analysed separately in natural populations. Such life cycles provide a unique opportunity to estimate the frequency of hybrid genotypes in diploid sporophytes and meiotic recombinant genotypes in haploid gametophytes allowing the effects of reproductive barriers preventing fertilization or preventing meiosis to be untangle. The level of hybridization between E. siliculosus and E. crouaniorum was quantified along the European coast. Clonal cultures (568 diploid, 336 haploid) isolated from field samples were genotyped using cytoplasmic and nuclear markers to estimate the frequency of hybrid genotypes in diploids and recombinant haploids. We identified admixed individuals using microsatellite loci, classical assignment methods and a newly developed Bayesian method (XPloidAssignment), which allows the analysis of populations that exhibit variations in ploidy level. Over all populations, the level of hybridization was estimated at 8.7%. Hybrids were exclusively observed in sympatric populations. More than 98% of hybrids were diploids (40% of which showed signs of aneuploidy) with a high frequency of rare alleles. The near absence of haploid recombinant hybrids demonstrates that the reproductive barriers are mostly postzygotic and suggests that abnormal chromosome segregation during meiosis following hybridization of species with different genome sizes could be a major cause of interspecific incompatibility in this system.

The eastern part of the Fertile Crescent concealed an unexpected route of olive (Olea europaea L.) differentiation

BACKGROUND AND AIMS

Olive is considered a native plant of the eastern side of the Mediterranean basin, from where it should have spread westward along the Mediterranean shores, while little is known about its diffusion in the eastern direction.

METHODS

Genetic diversity levels and population genetic structure of a wide set of olive ecotypes and varieties collected from several provinces of Iran, representing a high percentage of the entire olive resources present in the area, was screened with 49 chloroplast and ten nuclear simple sequence repeat markers, and coupled with archaeo-botanical and historical data on Mediterranean olive varieties. Approximate Bayesian Computation was applied to define the demographic history of olives including Iranian germplasm, and species distribution modelling was performed to understand the impact of the Late Quaternary on olive distribution.

KEY RESULTS

The results of the present study demonstrated that: (1) the climatic conditions of the last glacial maximum had an important role on the actual olive distribution, (2) all Iranian olive samples had the same maternal inheritance as Mediterranean cultivars, and (3) the nuclear gene flow from the Mediterranean basin to the Iranian plateau was almost absent, as well as the contribution of subspecies cuspidata to the diversity of Iranian olives.

CONCLUSIONS

Based on this evidence, a new scenario for the origin and distribution of this important fruit crop has been traced. The evaluation of olive trees growing in the eastern part of the Levant highlighted a new perspective on the spread and distribution of olive, suggesting two routes of olive differentiation, one westward, spreading along the Mediterranean basin, and another moving towards the east and reaching the Iranian plateau before its domestication.

Genetic diversity and population structure of the tsetse fly Glossina fuscipes fuscipes (Diptera: Glossinidae) in Northern Uganda: Implications for vector control

Uganda is the only country where the chronic and acute forms of human African Trypanosomiasis (HAT) or sleeping sickness both occur and are separated by < 100 km in areas north of Lake Kyoga. In Uganda, Glossina fuscipes fuscipes is the main vector of the Trypanosoma parasites responsible for these diseases as well for the animal African Trypanosomiasis (AAT), or Nagana. We used highly polymorphic microsatellite loci and a mitochondrial DNA (mtDNA) marker to provide fine scale spatial resolution of genetic structure of G. f. fuscipes from 42 sampling sites from the northern region of Uganda where a merger of the two disease belts is feared. Based on microsatellite analyses, we found that G. f. fuscipes in northern Uganda are structured into three distinct genetic clusters with varying degrees of interconnectivity among them. Based on genetic assignment and spatial location, we grouped the sampling sites into four genetic units corresponding to northwestern Uganda in the Albert Nile drainage, northeastern Uganda in the Lake Kyoga drainage, western Uganda in the Victoria Nile drainage, and a transition zone between the two northern genetic clusters characterized by high level of genetic admixture. An analysis using HYBRIDLAB supported a hybrid swarm model as most consistent with tsetse genotypes in these admixed samples. Results of mtDNA analyses revealed the presence of 30 haplotypes representing three main haplogroups, whose location broadly overlaps with the microsatellite defined clusters. Migration analyses based on microsatellites point to moderate migration among the northern units located in the Albert Nile, Achwa River, Okole River, and Lake Kyoga drainages, but not between the northern units and the Victoria Nile drainage in the west. Effective population size estimates were variable with low to moderate sizes in most populations and with evidence of recent population bottlenecks, especially in the northeast unit of the Lake Kyoga drainage. Our microsatellite and mtDNA based analyses indicate that G. f. fuscipes movement along the Achwa and Okole rivers may facilitate northwest expansion of the Rhodesiense disease belt in Uganda. We identified tsetse migration corridors and recommend a rolling carpet approach from south of Lake Kyoga northward to minimize disease dispersal and prevent vector re-colonization. Additionally, our findings highlight the need for continuing tsetse monitoring efforts during and after control.

Northern Uganda is an epidemiologically important region affected by human African trypanosomiasis (HAT) because it harbors both forms of the HAT disease (T. b. gambiense and T. b. rhodesiense). The geographic location of this region creates the risk that these distinct foci could merge, which would complicate diagnosis and treatment, and may result in recombination between the two parasite strains with as yet unknown consequences. Both strains require a tsetse fly vector for transmission, and in Uganda, G. f. fuscipes is the major vector of HAT. Controlling the vector remains one of the most effective strategies for controlling trypanosome parasites. However, vector control efforts may not be sustainable in terms of long term reduction in G. f. fuscipes populations due to population rebounds. Population genetics data can allow us to determine the likely source of population rebounds and to establish a more robust control strategy. In this study, we build on a previous broad spatial survey of G. f. fuscipes genetic structure in Uganda by adding more than 30 novel sampling sites that are strategically spaced across a region of northern Uganda that, for historical and political reasons, was severely understudied and faces particularly high disease risk. We identify natural population breaks, migration corridors and a hybrid zone with evidence of free interbreeding of G. f. fuscipes across the geographic region that spans the two HAT disease foci. We also find evidence of low effective population sizes and population bottlenecks in some areas that have been subjects of past control but remain regions of high tsetse density, which stresses the risk of population rebounds if monitoring is not explicitly incorporated into the control strategy. We use these results to make suggestions that will enhance the design and implementation of control activities in northern Uganda.

Population genetic and field-ecological analyses return similar estimates of dispersal over space and time in an endangered amphibian

The explosive growth of empirical population genetics has seen a proliferation of analytical methods leading to a steady increase in our ability to accurately measure key population parameters, including genetic isolation, effective population size, and gene flow, in natural systems. Assuming they yield similar results, population genetic methods offer an attractive complement to, or replacement of, traditional field‐ecological studies. However, empirical assessments of the concordance between direct field‐ecological and indirect population genetic studies of the same populations are uncommon in the literature. In this study, we investigate genetic isolation, rates of dispersal, and population sizes for the endangered California tiger salamander, Ambystoma californiense, across multiple breeding seasons in an intact vernal pool network. We then compare our molecular results to a previously published study based on multiyear, mark–recapture data from the same breeding sites. We found that field and genetic estimates of population size were only weakly correlated, but dispersal rates were remarkably congruent across studies and methods. In fact, dispersal probability functions derived from genetic data and traditional field‐ecological data were a significant match, suggesting that either method can be used effectively to assess population connectivity. These results provide one of the first explicit tests of the correspondence between landscape genetic and field‐ecological approaches to measuring functional population connectivity and suggest that even single‐year genetic samples can return biologically meaningful estimates of natural dispersal and gene flow.

Fine-scale population genetic structure of the Bengal tiger (Panthera tigris tigris) in a human-dominated western Terai Arc Landscape, India

Despite massive global conservation strategies, tiger populations continued to decline until recently, mainly due to habitat loss, human-animal conflicts, and poaching. These factors are known to affect the genetic characteristics of tiger populations and decrease local effective population sizes. The Terai Arc Landscape (TAL) at the foothills of the Himalaya is one of the 42 source sites of tigers around the globe. Therefore, information on how landscape features and anthropogenic factors affect the fine-scale spatial genetic structure and variation of tigers in TAL is needed to develop proper management strategies for achieving long-term conservation goals. We document, for the first time, the genetic characteristics of this tiger population by genotyping 71 tiger samples using 13 microsatellite markers from the western region of TAL (WTAL) of 1800 km2. Specifically, we aimed to estimate the genetic variability, population structure, and gene flow. The microsatellite markers indicated that the levels of allelic diversity (MNA = 6.6) and genetic variation (Ho = 0.50, HE = 0.64) were slightly lower than those reported previously in other Bengal tiger populations. We observed moderate gene flow and significant genetic differentiation (FST= 0.060) and identified the presence of cryptic genetic structure using Bayesian and non-Bayesian approaches. There was low and significantly asymmetric migration between the two main subpopulations of the Rajaji Tiger Reserve and the Corbett Tiger Reserve in WTAL. Sibship relationships indicated that the functionality of the corridor between these subpopulations may be retained if the quality of the habitat does not deteriorate. However, we found that gene flow is not adequate in view of changing land use matrices. We discuss the need to maintain connectivity by implementing the measures that have been suggested previously to minimize the level of human disturbance, including relocation of villages and industries, prevention of encroachment, and banning sand and boulder mining in the corridors.

Rethinking the influence of hydroelectric development on gene flow in a long-lived fish, the Lake Sturgeon Acipenser fulvescens

Many hydroelectric dams have been in place for 50 - >100 years, which for most fish species means that enough generations have passed for fragmentation induced divergence to have accumulated. However, for long-lived species such as Lake Sturgeon, Acipenser fulvescens, it should be possible to discriminate between historical population structuring and contemporary gene flow and improve the broader understanding of anthropogenic influence. On the Winnipeg River, Manitoba, two hypotheses were tested: 1) Measureable quantities of former reservoir dwelling Lake Sturgeon now reside downstream of the Slave Falls Generating Station, and 2) genetically differentiated populations of Lake Sturgeon occur upstream and downstream, a result of historical structuring. Genetic methods based on ten microsatellite markers were employed, and simulations were conducted to provide context. With regards to contemporary upstream to downstream contributions, the inclusion of length-at-age data proved informative. Both pairwise relatedness and Bayesian clustering analysis substantiated that fast-growing outliers, apparently entrained after residing in the upstream reservoir for several years, accounted for ~15% of the Lake Sturgeon 525–750 mm fork length captured downstream. With regards to historical structuring, upstream and downstream populations were found to be differentiated (F_ST = 0.011, and 0.013–0.014 when fast-growing outliers were excluded), and heterozygosity metrics were higher for downstream versus upstream juveniles. Historical asymmetric (downstream) gene flow in the vicinity of the generating station was the most logical explanation for the observed genetic structuring. In this section of the Winnipeg River, construction of a major dam does not appear to have fragmented a previously panmictic Lake Sturgeon population, but alterations to habitat may be influencing upstream to downstream contributions in unexpected ways.

Testing of microsatellite multiplexes for individual identification of Cape Parrots (Poicephalus robustus): paternity testing and monitoring trade

Background

Illegal trade in rare wildlife species is a major threat to many parrot species around the world. Wildlife forensics plays an important role in the preservation of endangered or threatened wildlife species. Identification of illegally harvested or traded animals through DNA techniques is one of the many methods used during forensic investigations. Natural populations of the South African endemic Cape Parrot (Poicephalus robustus) are negatively affected by the removal of eggs and chicks for the pet trade.

Methods

In this study, 16 microsatellite markers specifically designed for the South African endemic Cape Parrot (P. robustus) are assessed for their utility in forensic casework. Using these 16 loci, the genetic diversity of a subset of the captive Cape Parrot population was also assessed and compared to three wild Cape Parrot populations.

Results

It was determined that the full 16 locus panel has sufficient discriminatory power to be used in parentage analyses and can be used to determine if a bird has been bred in captivity and so can be legally traded or if it has been illegally removed from the wild. In cases where birds have been removed from the wild, this study suggests that a reduced 12 locus microsatellite panel has sufficient power to assign confiscated birds to geographic population of origin.

Discussion

The level of genetic diversity observed within the captive Cape Parrot population was similar to that observed in the wild populations, which suggests that the captive population is not suffering from decreased levels of genetic diversity. The captive Cape Parrots did however have double the number of private alleles compared to that observed in the most genetically diverse wild population. This is probably due to the presence of rare alleles present in the founder population, which has not been lost due to genetic drift, as many of the individuals tested in this study are F1–F3 wild descendants. The results from this study provide a suit of markers that can be used to aid conservation and law enforcement authorities to better control legal and illegal trade of this South African endemic.

Visualizations for genetic assignment analyses using the saddlepoint approximation method

We propose a method for visualizing genetic assignment data by characterizing the distribution of genetic profiles for each candidate source population. This method enhances the assignment method of Rannala and Mountain (1997) by calculating appropriate graph positions for individuals for which some genetic data are missing. An individual with missing data is positioned in the distributions of genetic profiles for a population according to its estimated quantile based on its available data. The quantiles of the genetic profile distribution for each population are calculated by approximating the cumulative distribution function (CDF) using the saddlepoint method, and then inverting the CDF to get the quantile function. The saddlepoint method also provides a way to visualize assignment results calculated using the leave-one-out procedure. This new method offers an advance upon assignment software such as geneclass2, which provides no visualization method, and is biologically more interpretable than the bar charts provided by the software structure. We show results from simulated data and apply the methods to microsatellite genotype data from ship rats (Rattus rattus) captured on the Great Barrier Island archipelago, New Zealand. The visualization method makes it straightforward to detect features of population structure and to judge the discriminative power of the genetic data for assigning individuals to source populations.

The use of microsatellite markers in Neotropical studies of wild birds: a literature review

Despite extensive habitat fragmentation, the Neotropical region possesses 30% of the world´s bird species. Microsatellites have remained one of the most popular genetic markers and have been used in ecological and conservation studies since the 1990's. We conducted a literature review comparing the number of papers published from January 1990 to July 2015 that used microsatellite markers for studies of wild birds in the Neotropical region, USA and some European countries. We assigned the articles to three categories of studies: population genetics, animal behavior/kinship analysis and the development of species-specific bird microsatellite markers. We also compared the studies in the Neotropics that used heterologous versus species-specific markers and provide a list of heterologous markers of utility in multiple birds. Despite the rich bird fauna in the Neotropics, the number of articles published represents only 5.6% of that published by the USA and selected European countries. Within the Neotropical region, Brazil possessed 60.5% of the total papers published, with the remaining 39.5% shared between five countries. We conclude that the lack of specialized laboratories and resources still represents a limit to microsatellite-based genetic studies of birds within the Neotropical region. To overcome these limitations, we suggest the use of heterologous microsatellite markers as a cost-effective and time-effective tool to assist ecological studies of wild birds.