Ecological modeling and pest population management: a possible and necessary connection in a changing world

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Evidence for temporal population replacement and the signature of ecological adaptation in a major Neotropical malaria vector in Amazonian Peru

The major Neotropical malaria vector, Anopheles darlingi, was reintroduced into the Iquitos, Loreto, Peru area during the early 1990s, where it displaced other anophelines and caused a major malaria epidemic. Since then, case numbers in Loreto have fluctuated, but annual increases have been reported since 2012. The population genetic structure of An. darlingi sampled before and after the introduction of long-lasting insecticidal nets (LLINs) was investigated to test the hypothesis of temporal population change (2006 vs. 2012). Current samples of An. darlingi were used to test the hypothesis of ecological adaptation to human modified (highway) compared with wild (riverine) habitat, linked to forest cover. In total, 693 An. darlingi from nine localities in Loreto, Peru area were genotyped using 13 microsatellite loci. To test the hypothesis of habitat differentiation in An. darlingi biting time patterns, HBR and EIR, four collections of An. darlingi from five localities (two riverine and three highway) were analysed. Analyses of microsatellite loci from seven (2006) and nine settlements (2012-2014) in the Iquitos area detected two distinctive populations with little overlap, although it is unclear whether this population replacement event is associated with LLIN distribution or climate. Within the 2012-2014 population two admixed subpopulations, A and B, were differentiated by habitat, with B significantly overrepresented in highway, and both in near-equal proportions in riverine. Both subpopulations had a signature of expansion and there was moderate genetic differentiation between them. Habitat and forest cover level had significant effects on HBR, such that Plasmodium transmission risk, as measured by EIR, in peridomestic riverine settlements was threefold higher than in peridomestic highway settlements. HBR was directly associated with available host biomass rather than forest cover. A population replacement event occurred between 2006 and 2012-2014, concurrently with LLIN distribution and a moderate El Niño event, and prior to an increase in malaria incidence. The likely drivers of this replacement cannot be determined with current data. The present-day An. darlingi population is composed of two highly admixed subpopulations, which appear to be in an early stage of differentiation, triggered by anthropogenic alterations to local habitat.

Spatial and temporal variation in population genetic structure of wild Nile tilapia (Oreochromis niloticus) across Africa

Background: Reconstructing the evolutionary history of a species is challenging. It often depends not only on the past biogeographic and climatic events but also the contemporary and ecological factors, such as current connectivity and habitat heterogeneity. In fact, these factors might interact with each other and shape the current species distribution. However, to what extent the current population genetic structure reflects the past and the contemporary factors is largely unknown. Here we investigated spatio-temporal genetic structures of Nile tilapia (Oreochromis niloticus) populations, across their natural distribution in Africa. While its large biogeographic distribution can cause genetic differentiation at the paleo-biogeographic scales, its restricted dispersal capacity might induce a strong genetic structure at micro-geographic scales. Results: Using nine microsatellite loci and 350 samples from ten natural populations, we found the highest genetic differentiation among the three ichthyofaunal provinces and regions (Ethiopian, Nilotic and Sudano-Sahelian) (R(ST) = 0.38 - 0.69). This result suggests the predominant effect of paleo-geographic events at macro-geographic scale. In addition, intermediate divergences were found between rivers and lakes within the regions, presumably reflecting relatively recent interruptions of gene flow between hydrographic basins (R(ST) = 0.24 - 0.32). The lowest differentiations were observed among connected populations within a basin (R(ST) = 0.015 in the Volta basin). Comparison of temporal sample series revealed subtle changes in the gene pools in a few generations (F = 0 - 0.053). The estimated effective population sizes were 23 - 143 and the estimated migration rate was moderate (m similar to 0.094 - 0.097) in the Volta populations. Conclusions: This study revealed clear hierarchical patterns of the population genetic structuring of O. niloticus in Africa. The effects of paleo-geographic and climatic events were predominant at macro-geographic scale, and the significant effect of geographic connectivity was detected at micro-geographic scale. The estimated effective population size, the moderate level of dispersal and the rapid temporal change in genetic composition might reflect a potential effect of life history strategy on population dynamics. This hypothesis deserves further investigation. The dynamic pattern revealed at micro-geographic and temporal scales appears important from a genetic resource management as well as from a biodiversity conservation point of view.

GENETIC VARIATION, LOCAL ADAPTATION AND POPULATION STRUCTURE IN NORTH AMERICAN RED OAK SPECIES, QUERCUS RUBRA L. AND Q. ELLIPSOIDALIS E. J. HILL

Forest trees, like oaks, rely on high levels of genetic variation to adapt to varying environmental conditions. Thus, genetic variation and its distribution are important for the long-term survival and adaptability of oak populations. Climate change is projected to lead to increased drought and fire events as well as a northward migration of tree species, including oaks. Additionally, decline in oak regeneration has become increasingly concerning since it may lead to decreased gene flow and increased inbreeding levels. This will in turn lead to lowered levels of genetic diversity, negatively affecting the growth and survival of populations. At the same time, populations at the species’ distribution edge, like those in this study, could possess important stores of genetic diversity and adaptive potential, while also being vulnerable to climatic or anthropogenic changes. A survey of the level and distribution of genetic variation and identification of potentially adaptive genes is needed since adaptive genetic variation is essential for their long-term survival. Oaks possess a remarkable characteristic in that they maintain their species identity and specific environmental adaptations despite their propensity to hybridize. Thus, in the face of interspecific gene flow, some areas of the genome remain differentiated due to selection. This characteristic allows the study of local environmental adaptation through genetic variation analyses. Furthermore, using genic markers with known putative functions makes it possible to link those differentiated markers to potential adaptive traits (e.g., flowering time, drought stress tolerance). Demographic processes like gene flow and genetic drift also play an important role in how genes (including adaptive genes) are maintained or spread. These processes are influenced by disturbances, both natural and anthropogenic. An examination of how genetic variation is geographically distributed can display how these genetic processes and geographical disturbances influence genetic variation patterns. For example, the spatial clustering of closely related trees could promote inbreeding with associated negative effects (inbreeding depression), if gene flow is limited. In turn this can have negative consequences for a species’ ability to adapt to changing environmental conditions. In contrast, interspecific hybridization may also allow the transfer of genes between species that increase their adaptive potential in a changing environment. I have studied the ecologically divergent, interfertile red oaks, Quercus rubra and Q. ellipsoidalis, to identify genes with potential roles in adaptation to abiotic stress through traits such as drought tolerance and flowering time, and to assess the level and distribution of genetic variation. I found evidence for moderate gene flow between the two species and low interspecific genetic differences at most genetic markers (Lind and Gailing 2013). However, the screening of genic markers with potential roles in phenology and drought tolerance led to the identification of a CONSTANS-like (COL) gene, a candidate gene for flowering time and growth. This marker, located in the coding region of the gene, was highly differentiated between the two species in multiple geographical areas, despite interspecific gene flow, and may play a role in reproductive isolation and adaptive divergence between the two species (Lind-Riehl et al. 2014). Since climate change could result in a northward migration of trees species like oaks, this gene could be important in maintaining species identity despite increased contact zones between species (e.g., increased gene flow). Finally I examined differences in spatial genetic structure (SGS) and genetic variation between species and populations subjected to different management strategies and natural disturbances. Diverse management activities combined with various natural disturbances as well as species specific life history traits influenced SGS patterns and inbreeding levels (Lind-Riehl and Gailing submitted).

Genetic Breeding And Diversity Of The Genus Passiflora: Progress And Perspectives In Molecular And Genetic Studies.

Despite the ecological and economic importance of passion fruit (Passiflora spp.), molecular markers have only recently been utilized in genetic studies of this genus. In addition, both basic genetic researches related to population studies and pre-breeding programs of passion fruit remain scarce for most Passiflora species. Considering the number of Passiflora species and the increasing use of these species as a resource for ornamental, medicinal, and food purposes, the aims of this review are the following: (i) to present the current condition of the passion fruit crop; (ii) to quantify the applications and effects of using molecular markers in studies of Passiflora; (iii) to present the contributions of genetic engineering for passion fruit culture; and (iv) to discuss the progress and perspectives of this research. Thus, the present review aims to summarize and discuss the relationship between historical and current progress on the culture, breeding, and molecular genetics of passion fruit.

TCF7L2 variant genotypes and type 2 diabetes risk in Brazil: significant association, but not a significant tool for risk stratification in the general population

Background: Genetic polymorphisms of the TCF7L2 gene are strongly associated with large increments in type 2 diabetes risk in different populations worldwide. In this study, we aimed to confirm the effect of the TCF7L2 polymorphism rs7903146 on diabetes risk in a Brazilian population and to assess the use of this genetic marker in improving diabetes risk prediction in the general population. Methods: We genotyped the single nucleotide polymorphisms (SNP) rs7903146 of the TCF7L2 gene in 560 patients with known coronary disease enrolled in the MASS II (Medicine, Angioplasty, or Surgery Study) Trial and in 1,449 residents of Vitoria, in Southeast Brazil. The associations of this gene variant to diabetes risk and metabolic characteristics in these two different populations were analyzed. To access the potential benefit of using this marker for diabetes risk prediction in the general population we analyzed the impact of this genetic variant on a validated diabetes risk prediction tool based on clinical characteristics developed for the Brazilian general population. Results: SNP rs7903146 of the TCF7L2 gene was significantly associated with type 2 diabetes in the MASS-II population (OR = 1.57 per T allele, p = 0.0032), confirming, in the Brazilian population, previous reports of the literature. Addition of this polymorphism to an established clinical risk prediction score did not increased model accuracy (both area under ROC curve equal to 0.776). Conclusion: TCF7L2 rs7903146 T allele is associated with a 1.57 increased risk for type 2 diabetes in a Brazilian cohort of patients with known coronary heart disease. However, the inclusion of this polymorphism in a risk prediction tool developed for the general population resulted in no improvement of performance. This is the first study, to our knowledge, that has confirmed this recent association in a South American population and adds to the great consistency of this finding in studies around the world. Finally, confirming the biological association of a genetic marker does not guarantee improvement on already established screening tools based solely on demographic variables.

Successful maintenance of a stingless bee population despite a severe genetic bottleneck

Stingless bees play an important ecological role as pollinators of many wild plant species in the tropics and have significant potential for the pollination of agricultural crops. Nevertheless, conservation efforts as well as commercial breeding programmes require better guidelines on the amount of genetic variation that is needed to maintain viable populations. In this context, we carried out a long-term genetic study on the stingless bee Melipona scutellaris to evaluate the population viability consequences of prolonged breeding from a small number of founder colonies. In particular, it was artificially imposed a genetic bottleneck by setting up a population starting from only two founder colonies, and continued breeding from it for a period of over 10 years in a location outside its natural area of occurrence. We show that despite a great reduction in the number of alleles present at both neutral microsatellite loci and the sex-determining locus relative to its natural source population, and an increased frequency in the production of sterile diploid males, the genetically impoverished population could be successfully bred and maintained for at least 10 years. This shows that in stingless bees, breeding from a small stock of colonies may have less severe consequences than previously suspected. In addition, we provide a simulation model to determine the number of colonies that are needed to maintain a certain number of sex alleles in a population, thereby providing useful guidelines for stingless bee breeding and conservation efforts.

Social structures, genetic structures and dispersal strategies in Australian rabbit (Oryctolagus cuniculus) populations

In this study, the pattern of movement of young male and female rabbits and the genetic structures present in adult male and female populations in four habitats was examined. The level of philopatry in young animals was found to vary between 18-90% for males and 32-95% for females in different populations. It was skewed, with more males dispersing than females in some but not all populations. Analysis of allozyme data using spatial autocorrelation showed that adult females from the same social group, unlike males, were significantly related in four of the five populations studied. Changes in genetic structure and rate of dispersal were measured before and during the recovery of a population that was artificially reduced in size. There were changes in the rate and distance of dispersal with density and sex. Subadults of both sexes moved further in the first year post crash (low density) than in the following years. While the level of dispersal for females was lower than that of the males for the first 3 years, thereafter (high density) both sexes showed similar, low levels of dispersal (20%). The density at which young animals switch behaviour between dispersal and philopatry differed for males and females. The level of genetic structuring in adult females was high in the precrash population, reduced in the first year post crash and undetectable in the second year. Dispersal behaviour of rabbits both affects the genetic structure of the population and changes with conditions. Over a wide range of levels of philopatry, genetic structuring is present in the adult female, but not the male population. Consequently, though genetic structuring is present, it does not lead to inbreeding. More long-distance movements are found in low-density populations, even though vacant warrens are available near birth warrens. The distances moved decreased as density increased. Calculation of the effective population size (N-e) shows that changes in dispersal distance offset changes in density, so that N-e remains constant.

Population structure in high shore littorinids: contrast between riprap and rocky shores

World Congress of Malacology, Ponta Delgada, July 22-28, 2013.

Homicide and impunity: an ecological analysis at state level in Brazil

OBJECTIVE: To assess a new impunity index and variables that have been found to predict variation in homicide rates in other geographical levels as predictive of state-level homicide rates in Brazil. METHODS: This was a cross-sectional ecological study. Data from the mortality information system relating to the 27 Brazilian states for the years 1996 to 2005 were analyzed. The outcome variables were taken to be homicide victim rates in 2005, for the entire population and for men aged 20-29 years. Measurements of economic and social development, economic inequality, demographic structure and life expectancy were analyzed as predictors. An "impunity index", calculated as the total number of homicides between 1996 and 2005 divided by the number of individuals in prison in 2007, was constructed. The data were analyzed by means of simple linear regression and negative binomial regression. RESULTS: In 2005, state-level crude total homicide rates ranged from 11 to 51 per 100,000; for young men, they ranged from 39 to 241. The impunity index ranged from 0.4 to 3.5 and was the most important predictor of this variability. From negative binomial regression, it was estimated that the homicide victim rate among young males increased by 50% for every increase of one point in this ratio. CONCLUSIONS: Classic predictive factors were not associated with homicides in this analysis of state-level variation in Brazil. However, the impunity index indicated that the greater the impunity, the higher the homicide rate.

Ecological-economic modelling of interactions between wild and commercial bees and pesticide use

The decline in extent of wild pollinators in recent years has been partly associated with changing farm practices and in particular with increase of pesticide use. In this paper we combine ecological modelling with economic analysis of a single farm output under the assumption that both pollination and pest control are essential inputs. We show that the drive to increase farm output can lead to a local decline in the wild bee population. Commercial bees are often considered an alternative to wild pollinators, but we show that their introduction can lead to further decline and finally local extinction of wild bees. The transitions between different outcomes are characterised by threshold behaviour and are potentially difficult to predict and detect in advance. Small changes in economic (input prices) and ecological (wild bees carrying capacity and effect of pesticides on bees) can move the economic-ecological system beyond the extinction threshold. We also show that increasing the pesticide price or decreasing the commercial bee price might lead to reestablishment of wild bees following their local extinction. Thus, we demonstrate the importance of combining ecological modelling with economics to study the provision of ecosystem services and to inform sustainable management of ecosystem service providers.

How life history and demography promote or inhibit the evolution of helping behaviours.

In natural populations, dispersal tends to be limited so that individuals are in local competition with their neighbours. As a consequence, most behaviours tend to have a social component, e.g. they can be selfish, spiteful, cooperative or altruistic as usually considered in social evolutionary theory. How social behaviours translate into fitness costs and benefits depends considerably on life-history features, as well as on local demographic and ecological conditions. Over the last four decades, evolutionists have been able to explore many of the consequences of these factors for the evolution of social behaviours. In this paper, we first recall the main theoretical concepts required to understand social evolution. We then discuss how life history, demography and ecology promote or inhibit the evolution of helping behaviours, but the arguments developed for helping can be extended to essentially any social trait. The analysis suggests that, on a theoretical level, it is possible to contrast three critical benefit-to-cost ratios beyond which costly helping is selected for (three quantitative rules for the evolution of altruism). But comparison between theoretical results and empirical data has always been difficult in the literature, partly because of the perennial question of the scale at which relatedness should be measured under localized dispersal. We then provide three answers to this question.

Distinct population structure for co-occurring Anopheles goeldii and Anopheles triannulatus in Amazonian Brazil

To evaluate whether environmental heterogeneity contributes to the genetic heterogeneity in Anopheles triannulatus, larval habitat characteristics across the Brazilian states of Roraima and Pará and genetic sequences were examined. A comparison with Anopheles goeldii was utilised to determine whether high genetic diversity was unique to An. triannulatus. Student t test and analysis of variance found no differences in habitat characteristics between the species. Analysis of population structure of An. triannulatus and An. goeldii revealed distinct demographic histories in a largely overlapping geographic range. Cytochrome oxidase I sequence parsimony networks found geographic clustering for both species; however nuclear marker networks depicted An. triannulatus with a more complex history of fragmentation, secondary contact and recent divergence. Evidence of Pleistocene expansions suggests both species are more likely to be genetically structured by geographic and ecological barriers than demography. We hypothesise that niche partitioning is a driving force for diversity, particularly in An. triannulatus.

From the Apennines to the Alps: colonization genetics of the naturally expanding Italian wolf (Canis lupus) population

Wolves in Italy strongly declined in the past and were confined south of the Alps since the turn of the last century, reduced in the 1970s to approximately 100 individuals surviving in two fragmented subpopulations in the central-southern Apennines. The Italian wolves are presently expanding in the Apennines, and started to recolonize the western Alps in Italy, France and Switzerland about 16 years ago. In this study, we used a population genetic approach to elucidate some aspects of the wolf recolonization process. DNA extracted from 3068 tissue and scat samples collected in the Apennines (the source populations) and in the Alps (the colony), were genotyped at 12 microsatellite loci aiming to assess (i) the strength of the bottleneck and founder effects during the onset of colonization; (ii) the rates of gene flow between source and colony; and (iii) the minimum number of colonizers that are needed to explain the genetic variability observed in the colony. We identified a total of 435 distinct wolf genotypes, which showed that wolves in the Alps: (i) have significantly lower genetic diversity (heterozygosity, allelic richness, number of private alleles) than wolves in the Apennines; (ii) are genetically distinct using pairwise F(ST) values, population assignment test and Bayesian clustering; (iii) are not in genetic equilibrium (significant bottleneck test). Spatial autocorrelations are significant among samples separated up to c. 230 km, roughly correspondent to the apparent gap in permanent wolf presence between the Alps and north Apennines. The estimated number of first-generation migrants indicates that migration has been unidirectional and male-biased, from the Apennines to the Alps, and that wolves in southern Italy did not contribute to the Alpine population. These results suggest that: (i) the Alps were colonized by a few long-range migrating wolves originating in the north Apennine subpopulation; (ii) during the colonization process there has been a moderate bottleneck; and (iii) gene flow between sources and colonies was moderate (corresponding to 1.25-2.50 wolves per generation), despite high potential for dispersal. Bottleneck simulations showed that a total of c. 8-16 effective founders are needed to explain the genetic diversity observed in the Alps. Levels of genetic diversity in the expanding Alpine wolf population, and the permanence of genetic structuring, will depend on the future rates of gene flow among distinct wolf subpopulation fragments.