Bayesian clustering of farm types using the mixtures model

A Bayesian method of classifying observations that are assumed to come from a number of distinct subpopulations is outlined. The method is illustrated with simulated data and applied to the classification of farms according to their level and variability of income. The resultant classification shows a greater diversity of technical charactersitics within farm types than is conventionally the case. The range of mean farm income between groups in the new classification is wider than that of the conventional method and the variability of income within groups is narrower. Results show that the highest income group in 2000 included large specialist dairy farmers and pig and poultry producers, whilst in 2001 it included large and small specialist dairy farms and large mixed dairy and arable farms. In both years the lowest income group is dominated by non-milk producing livestock farms.

Modelling Breaks and Clusters in the Steady States of Macroeconomic Variables

Macroeconomists working with multivariate models typically face uncertainty over which (if any) of their variables have long run steady states which are subject to breaks. Furthermore, the nature of the break process is often unknown. In this paper, we draw on methods from the Bayesian clustering literature to develop an econometric methodology which: i) finds groups of variables which have the same number of breaks; and ii) determines the nature of the break process within each group. We present an application involving a five-variate steady-state VAR.

Repeated unidirectional introgression towards Populus balsamifera in contact zones of exotic and native poplars.

As the evolutionary significance of hybridization is largely dictated by its extent beyond the first generation, we broadly surveyed patterns of introgression across a sympatric zone of two native poplars (Populus balsamifera, Populus deltoides) in Quebec, Canada within which European exotic Populus nigra and its hybrids have been extensively planted since the 1800s. Single nucleotide polymorphisms (SNPs) that appeared fixed within each species were characterized by DNA-sequencing pools of pure individuals. Thirty-five of these diagnostic SNPs were employed in a high-throughput assay that genotyped 635 trees of different age classes, sampled from 15 sites with various degrees of anthropogenic disturbance. The degree of admixture within sampled trees was then assessed through Bayesian clustering of genotypes. Hybrids were present in seven of the populations, with 2.4% of all sampled trees showing spontaneous admixture. Sites with hybrids were significantly more disturbed than pure stands, while hybrids comprised both immature juveniles and trees of reproductive age. All three possible F1s were detected. Advanced-generation hybrids were consistently biased towards P. balsamifera regardless of whether hybridization had occurred with P. deltoides or P. nigra. Gene exchange between P. deltoides and P. nigra was not detected beyond the F1 generation; however, detection of a trihybrid demonstrates that even this apparent reproductive isolation does not necessarily result in an evolutionary dead end. Collectively, results demonstrate the natural fertility of hybrid poplars and suggest that introduced genes could potentially affect the genetic integrity of native trees, similar to that arising from introgression between natives.

Genetic and phenotypic population divergence on a microgeographic scale in brown trout.

Salmonid populations of many rivers are rapidly declining. One possible explanation is that habitat fragmentation increases genetic drift and reduces the populations' potential to adapt to changing environmental conditions. We measured the genetic and eco-morphological diversity of brown trout (Salmo trutta) in a Swiss stream system, using multivariate statistics and Bayesian clustering. We found large genetic and phenotypic variation within only 40 km of stream length. Eighty-eight percent of all pairwise F(ST) comparisons and 50% of the population comparisons in body shape were significant. High success rates of population assignment tests confirmed the distinctiveness of populations in both genotype and phenotype. Spatial analysis revealed that divergence increased with waterway distance, the number of weirs, and stretches of poor habitat between sampling locations, but effects of isolation-by-distance and habitat fragmentation could not be fully disentangled. Stocking intensity varied between streams but did not appear to erode genetic diversity within populations. A lack of association between phenotypic and genetic divergence points to a role of local adaptation or phenotypically plastic responses to habitat heterogeneity. Indeed, body shape could be largely explained by topographic stream slope, and variation in overall phenotype matched the flow regimes of the respective habitats.

Hybridization and Population Structure of the Culex Pipiens Complex in the Islands of Macaronesia

The Culex pipiens complex includes two widespread mosquito vector species, Cx. pipiens and Cx. quinquefasciatus. The distribution of these species varies in latitude, with the former being present in temperate regions and the latter in tropical and subtropical regions. However, their distribution range overlaps in certain areas and interspecific hybridization has been documented. Genetic introgression between these species may have epidemiological repercussions for West Nile virus (WNV) transmission. Bayesian clustering analysis based on multilocus genotypes of 12 microsatellites was used to determine levels of hybridization between these two species in Macaronesian islands, the only contact zone described in West Africa. The distribution of the two species reflects both the islands’ biogeography and historical aspects of human colonization. Madeira Island displayed a homogenous population of Cx. pipiens, whereas Cape Verde showed a more intriguing scenario with extensive hybridization. In the islands of Brava and Santiago, only Cx. quinquefasciatus was found, while in Fogo and Maio high hybrid rates (~40%) between the two species were detected. Within the admixed populations, second-generation hybrids (~50%) were identified suggesting a lack of isolation mechanisms. The observed levels of hybridization may locally potentiate the transmission to humans of zoonotic arboviruses such as WNV.

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.

Reduced genetic diversity, increased isolation and multiple introductions of invasive giant hogweed in the western Swiss Alps.

Abstract The giant hogweed (Heracleum mantegazzianum) has successfully invaded 19 European countries as well as parts of North America. It has become a problematic species due to its ability to displace native flora and to cause public health hazards. Applying population genetics to species invasion can help reconstruct invasion history and may promote more efficient management practice. We thus analysed levels of genetic variation and population genetic structure of H. mantegazzianum in an invaded area of the western Swiss Alps as well as in its native range (the Caucasus), using eight nuclear microsatellite loci together with plastid DNA markers and sequences. On both nuclear and plastid genomes, native populations exhibited significantly higher levels of genetic diversity compared to invasive populations, confirming an important founder event during the invasion process. Invasive populations were also significantly more differentiated than native populations. Bayesian clustering analysis identified five clusters in the native range that corresponded to geographically and ecologically separated groups. In the invaded range, 10 clusters occurred. Unlike native populations, invasive clusters were characterized by a mosaic pattern in the landscape, possibly caused by anthropogenic dispersal of the species via roads and direct collection for ornamental purposes. Lastly, our analyses revealed four main divergent groups in the western Swiss Alps, likely as a consequence of multiple independent establishments of H. mantegazzianum.

Insights into the biogeographical history of the Lower Guinea Forest Domain: evidence for the role of refugia in the intraspecific differentiation of Aucoumea klaineana.

Determining the biogeographical histories of rainforests is central to our understanding of the present distribution of tropical biodiversity. Ice age fragmentation of central African rainforests strongly influenced species distributions. Elevated areas characterized by higher species richness and endemism have been postulated to be Pleistocene forest refugia. However, it is often difficult to separate the effects of history and of present-day ecological conditions on diversity patterns at the interspecific level. Intraspecific genetic variation could yield new insights into history, because refugia hypotheses predict patterns not expected on the basis of contemporary environmental dynamics. Here, we test geographically explicit hypotheses of vicariance associated with the presence of putative refugia and provide clues about their location. We intensively sampled populations of Aucoumea klaineana, a forest tree sensitive to forest fragmentation, throughout its geographical range. Characterizing variation at 10 nuclear microsatellite loci, we were able to obtain phylogeographic data of unprecedented detail for this region. Using Bayesian clustering approaches, we demonstrated the presence of four differentiated genetic units. Their distribution matched that of forest refugia postulated from patterns of species richness and endemism. Our data also show differences in diversity dynamics at leading and trailing edges of the species' shifting distribution. Our results confirm predictions based on refugia hypotheses and cannot be explained on the basis of present-day ecological conditions.

Distribution and population genetic variation of cryptic species of the Alpine mayfly Baetis alpinus (Ephemeroptera: Baetidae) in the Central Alps.

BACKGROUND: Many species contain evolutionarily distinct groups that are genetically highly differentiated but morphologically difficult to distinguish (i.e., cryptic species). The presence of cryptic species poses significant challenges for the accurate assessment of biodiversity and, if unrecognized, may lead to erroneous inferences in many fields of biological research and conservation. RESULTS: We tested for cryptic genetic variation within the broadly distributed alpine mayfly Baetis alpinus across several major European drainages in the central Alps. Bayesian clustering and multivariate analyses of nuclear microsatellite loci, combined with phylogenetic analyses of mitochondrial DNA, were used to assess population genetic structure and diversity. We identified two genetically highly differentiated lineages (A and B) that had no obvious differences in regional distribution patterns, and occurred in local sympatry. Furthermore, the two lineages differed in relative abundance, overall levels of genetic diversity as well as patterns of population structure: lineage A was abundant, widely distributed and had a higher level of genetic variation, whereas lineage B was less abundant, more prevalent in spring-fed tributaries than glacier-fed streams and restricted to high elevations. Subsequent morphological analyses revealed that traits previously acknowledged as intraspecific variation of B. alpinus in fact segregated these two lineages. CONCLUSIONS: Taken together, our findings indicate that even common and apparently ecologically well-studied species may consist of reproductively isolated units, with distinct evolutionary histories and likely different ecology and evolutionary potential. These findings emphasize the need to investigate hidden diversity even in well-known species to allow for appropriate assessment of biological diversity and conservation measures.

Analyse bayésienne et classification pour modèles continus modifiés à zéro

Les modèles à sur-représentation de zéros discrets et continus ont une large gamme d'applications et leurs propriétés sont bien connues. Bien qu'il existe des travaux portant sur les modèles discrets à sous-représentation de zéro et modifiés à zéro, la formulation usuelle des modèles continus à sur-représentation -- un mélange entre une densité continue et une masse de Dirac -- empêche de les généraliser afin de couvrir le cas de la sous-représentation de zéros. Une formulation alternative des modèles continus à sur-représentation de zéros, pouvant aisément être généralisée au cas de la sous-représentation, est présentée ici. L'estimation est d'abord abordée sous le paradigme classique, et plusieurs méthodes d'obtention des estimateurs du maximum de vraisemblance sont proposées. Le problème de l'estimation ponctuelle est également considéré du point de vue bayésien. Des tests d'hypothèses classiques et bayésiens visant à déterminer si des données sont à sur- ou sous-représentation de zéros sont présentées. Les méthodes d'estimation et de tests sont aussi évaluées au moyen d'études de simulation et appliquées à des données de précipitation agrégées. Les diverses méthodes s'accordent sur la sous-représentation de zéros des données, démontrant la pertinence du modèle proposé. Nous considérons ensuite la classification d'échantillons de données à sous-représentation de zéros. De telles données étant fortement non normales, il est possible de croire que les méthodes courantes de détermination du nombre de grappes s'avèrent peu performantes. Nous affirmons que la classification bayésienne, basée sur la distribution marginale des observations, tiendrait compte des particularités du modèle, ce qui se traduirait par une meilleure performance. Plusieurs méthodes de classification sont comparées au moyen d'une étude de simulation, et la méthode proposée est appliquée à des données de précipitation agrégées provenant de 28 stations de mesure en Colombie-Britannique.

Characterisation of microsatellite markers for fig-pollinating wasps in the Pleistodontes imperialis species complex

We characterised a set of nine polymorphic microsatellite loci for Pleistodontes imperialis sp. 1, the pollinator wasp of Port Jackson fig (Ficus rubiginosa) in south-eastern Australia. Characterisation was performed on 30 female individuals collected from a population in Sydney, Australia. The average number of alleles per locus was 7.33, and eight loci were not in Hardy–Weinberg equilibrium. This was expected as fig wasps are known to be highly inbred. A test of genetic differentiation between two natural populations of P. imperialis sp. 1 (Sydney and Newcastle, Australia – some 120 km apart) yielded a very low FST value of 0.012, suggesting considerable gene flow. Bayesian clustering analysis using TESS 2.3.1, which does not assume Hardy–Weinberg equilibrium, however, indicated potential spatial substructuring between the Sydney and Newcastle populations, as well as within the Sydney population. The described loci were also characterised for two other species in the P. imperialis complex: P. imperialis sp. 2 (Townsville, Australia) and P. imperialis sp. 4 (Brisbane, Australia). Seven and six of the nine loci were polymorphic for P. imperialis sp. 2 and P.imperialis sp. 4, respectively.

Analysis of the genetic diversity and structure across a wide range of germplasm reveals prominent gene flow in apple at the European level

Abstract Background: The amount and structure of genetic diversity in dessert apple germplasm conserved at a European level is mostly unknown, since all diversity studies conducted in Europe until now have been performed on regional or national collections. Here, we applied a common set of 16 SSR markers to genotype more than 2,400 accessions across 14 collections representing three broad European geographic regions (North+East, West and South) with the aim to analyze the extent, distribution and structure of variation in the apple genetic resources in Europe. Results: A Bayesian model-based clustering approach showed that diversity was organized in three groups, although these were only moderately differentiated (FST=0.031). A nested Bayesian clustering approach allowed identification of subgroups which revealed internal patterns of substructure within the groups, allowing a finer delineation of the variation into eight subgroups (FST=0.044). The first level of stratification revealed an asymmetric division of the germplasm among the three groups, and a clear association was found with the geographical regions of origin of the cultivars. The substructure revealed clear partitioning of genetic groups among countries, but also interesting associations between subgroups and breeding purposes of recent cultivars or particular usage such as cider production. Additional parentage analyses allowed us to identify both putative parents of more than 40 old and/or local cultivars giving interesting insights in the pedigree of some emblematic cultivars. Conclusions: The variation found at group and sub-group levels may reflect a combination of historical processes of migration/selection and adaptive factors to diverse agricultural environments that, together with genetic drift, have resulted in extensive genetic variation but limited population structure. The European dessert apple germplasm represents an important source of genetic diversity with a strong historical and patrimonial value. The present work thus constitutes a decisive step in the field of conservation genetics. Moreover, the obtained data can be used for defining a European apple core collection useful for further identification of genomic regions associated with commercially important horticultural traits in apple through genome-wide association studies.

Cryptic Genetic Diversity Is Paramount in Small-Bodied Amphibians of the Genus Euparkerella (Anura: Craugastoridae) Endemic to the Brazilian Atlantic Forest

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Genetic diversity and population structure of the New World screwworm fly from the Amazon region of Brazil

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

Prospecção de marcadores microssatélites, análise de viriabilidade e estrutura genética populacional de arara-azul-grande Anodorhynchus hyacinthinus (Latham, 1790) com ênfase na região do mosaico dos Carajás-PA

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)