Friday, June 29, 2018

Weekend reads

Long weekend for all the Canadians which means more time to read and certainly no post on Monday.

Knowledge of community structure within an ecosystem is essential when trying to understand the function and importance of the system and when making related management decisions. Within the larger ecosystem, microhabitats play an important role by providing inhabitants with a subset of available resources. On coral reefs, cryptobenthic fishes encompass many groups and make up an important proportion of the biodiversity. However, these fishes are relatively small, exhibit extreme visual or behavioral camouflage, and, therefore, are often overlooked. We examined the differences in fish community structure between three common reef microhabitats (live hard coral, dead coral rubble, and sand) using ichthyocide stations in the central Red Sea. Using a combination of morphological and genetic (cytochrome oxidase I (COI) barcoding) techniques, we identified 326 individuals representing 73 species spread across 17 families, from fifteen 1 m2 quadrats. Fish assemblages in the three microhabitats were significantly different from each other. Rubble microhabitats yielded the highest levels of fish abundance, richness, and diversity, followed by hard coral, and then sand. The results show that benthic composition, even at a small scale, influences cryptobenthic communities. This study also provides new COI sequence data to public databases, in order to further the research of cryptobenthic fishes in the Red Sea region.

The cytochrome c oxidase subunit I (cox1) gene is the main mitochondrial molecular marker playing a pivotal role in phylogenetic research and is a crucial barcode sequence. Folmer's "universal" primers designed to amplify this gene in metazoan invertebrates allowed quick and easy barcode and phylogenetic analysis. On the other hand, the increase in the number of studies on barcoding leads to more frequent publishing of incorrect sequences, due to amplification of non-target taxa, and insufficient analysis of the obtained sequences. Consequently, some sequences deposited in genetic databases are incorrectly described as obtained from invertebrates, while being in fact bacterial sequences. In our study, in which we used Folmer's primers to amplify COI sequences of the crustacean fairy shrimp Branchipus schaefferi (Fischer 1834), we also obtained COI sequences of microbial contaminants from Aeromonas sp. However, when we searched the GenBank database for sequences closely matching these contaminations we found entries described as representatives of Gastrotricha and Mollusca. When these entries were compared with other sequences bearing the same names in the database, the genetic distance between the incorrect and correct sequences amplified from the same species was c.a. 65%. Although the responsibility for the correct molecular identification of species rests on researchers, the errors found in already published sequences data have not been re-evaluated so far. On the basis of the standard sampling technique we have estimated with 95% probability that the chances of finding incorrectly described metazoan sequences in the GenBank depend on the systematic group, and variety from less than 1% (Mollusca and Arthropoda) up to 6.9% (Gastrotricha). Consequently, the increasing popularity of DNA barcoding and metabarcoding analysis may lead to overestimation of species diversity. Finally, the study also discusses the sources of the problems with amplification of non-target sequences.

DNA metabarcoding is increasingly used in dietary studies to estimate diversity, composition, and frequency of occurrence of prey items. However, few studies have assessed how technical and biological replication affect the accuracy of diet estimates. This study addresses these issues using the European free-tailed bat Tadarida teniotis, involving high-throughput sequencing of a small fragment of the COI gene in 15 separate faecal pellets and a 15-pellet pool per each of 20 bats. We investigated how diet descriptors were affected by variability among (i) individuals, (ii) pellets of each individual, and (iii) PCRs of each pellet. In addition, we investigated the impact of (iv) analysing separate pellets versus pellet pools. We found that diet diversity estimates increased steadily with the number of pellets analysed per individual, with seven pellets required to detect ~80% of prey species. Most variation in diet composition was associated with differences among individual bats, followed by pellets per individual, and PCRs per pellet. The accuracy of frequency of occurrence estimates increased with the number of pellets analysed per bat, with the highest error rates recorded for prey consumed infrequently by many individuals. Pools provided poor estimates of diet diversity and frequency of occurrence, which were comparable to analysing a single pellet per individual, and consistently missed the less common prey items. Overall, our results stress that maximizing biological replication is critical in dietary metabarcoding studies, and emphasize that analysing several samples per individual rather than pooled samples produce more accurate results.

DNA metabarcoding is a technique used to survey biodiversity in many ecological settings, but there are doubts about whether it can provide quantitative results, i.e. the proportions of each species in the mixture as opposed to a species list. While there are several experimental studies that report quantitative metabarcoding results, there are a similar number that fail to do so. Here we provide the rationale to understand under what circumstances the technique can be quantitative. Basically, we simulate a mixture of DNA of S species with a defined initial abundance distribution. In the simulated PCR, each species increases its concentration following a certain amplification efficiency. The final DNA concentration will reflect the initial one when the efficiency is similar for all species; otherwise, the initial and final DNA concentrations would be poorly related. Although there are many known factors that modulate amplification efficiency, we focused on the number of primer-template mismatches, arguably the most important one. We used 15 common primers pairs targeting the mitochondrial COI region and the mitogenomes of ca. 1200 insect species. The results showed that some primers pairs produced quantitative results under most circumstances, whereas some other primers failed to do so. Many species, and a high diversity within the mixture, helped the metabarcoding to be quantitative. In conclusion, depending on the primer pair used in the PCR amplification and on the characteristics of the mixture analysed (i.e., high species richness, low evenness), DNA metabarcoding can provide a quantitative estimate of the relative abundances of different species.

Marine meiofauna comprises up to 22 phyla. Its morphological identification requires time and taxonomists' expertise, and molecular tools can make this task faster. We aim to disentangle meiofaunal diversity patterns at Araçá Bay by applying a model selection approach and estimating the effectiveness of metabarcoding (18S rDNA) and morphological methods for estimating the response of meiofauna diversity in small-scale interactions with environmental variables. A rarefaction curve indicated that ten samples were sufficient for estimating the total number of meiofauna OTUs in a tidal flat. In both approaches, richness was predicted by mean sand percentage, sediment sorting, and bacteria concentration. Nematode genera composition differed significantly between approaches, the result of taxonomic mismatch in the genetic database. The similarity between the model selected for diversity descriptors, the richness of nematode genera and meiofauna composition emphasized the utility of predictive models for metabarcoding estimates to detect small-scale interactions of these organisms.

Background: In light of the current biodiversity crisis, DNA barcoding is developing into an essential tool to quantify state shifts in global ecosystems. Current barcoding protocols often rely on short amplicon sequences, which yield accurate identification of biological entities in a community, but provide limited phylogenetic resolution across broad taxonomic scales. However, the phylogenetic structure of communities is an essential component of biodiversity. Consequently, a barcoding approach is required that unites robust taxonomic assignment power and high phylogenetic utility. A possible solution is offered by sequencing long ribosomal DNA (rDNA) amplicons on the MinION platform (Oxford Nanopore Technologies). Results: Using a dataset of various animal and plant species, with a focus on arthropods, we assemble a pipeline for long rDNA barcode analysis and introduce a new software (MiniBar) to demultiplex dual indexed nanopore reads. We find excellent phylogenetic and taxonomic resolution offered by long rDNA sequences across broad taxonomic scales. We highlight the simplicity of our approach by field barcoding with a miniaturized, mobile laboratory in a remote rainforest. We also test the utility of long rDNA amplicons for analysis of community diversity through metabarcoding and find that they recover highly skewed diversity estimates. Conclusions: Sequencing dual indexed, long rDNA amplicons on the MinION platform is a straightforward, cost effective, portable and universal approach for eukaryote DNA barcoding. Long rDNA amplicons scale up DNA barcoding by enabling the accurate recovery of taxonomic and phylogenetic diversity. However, bulk community analyses using long-read approaches may introduce biases and will require further exploration.

Background. Knowledge on the globally outstanding Amazonian biodiversity and its environmental determinants stems almost exclusively from aboveground organisms, notably plants. In contrast, the environmental factors and habitat preferences that drive diversity patterns for micro-organisms in the ground remain elusive, despite the fact that micro-organisms constitute the overwhelming majority of life forms in any given location, in terms of both diversity and abundance. Here we address how the diversity and community turnover of operational taxonomic units (OTU) of micro-organisms in soil and litter respond to soil physicochemical properties; whether OTU diversities and community composition in soil and litter are correlated with each other; and whether they respond in a similar way to soil properties. Methods. We used recently inferred OTUs from high-throughput metabarcoding of the 16S (prokaryotes) and 18S (eukaryotes) genes to estimate OTU diversity (OTU richness and effective number of OTUs) and community composition for prokaryotes and eukaryotes in soil and litter across four localities in Brazilian Amazonia. All analyses were run separately for prokaryote and eukaryote OTUs, and for each group using both presence-absence and abundance data. Combining these with novel data on soil chemical and physical properties, we identify abiotic correlates of soil and litter micro-organism diversity and community structure using regression, ordination, and variance partitioning analysis. Results. Soil organic carbon content was the strongest factor explaining OTU diversity (negative correlation) and pH was the strongest factor explaining turnover for prokaryotes and eukaryotes in both soil and litter. We found significant effects also for other soil variables, including both chemical and physical properties. The correlation between OTU diversity in litter and in soil was non-significant for eukaryotes and weak for prokaryotes, suggesting that diversity of in one substrate should not be used as a proxy for diversity in the other. The community compositions of both prokaryotes and eukaryotes were more separated for habitat type than for substrate (soil and litter). Discussion. In spite of the limited sampling (four localities, 39 plots), our results provide a broad-scale view of the physical and chemical correlations of soil and litter biodiversity in a longitudinal transect across the world’s largest rainforest. Our methods help to understand links between soil properties, OTU diversity patterns, and community composition and turnover. The lack of strong correlation between OTU diversity in litter and in soil suggests independence of diversity drives of these substrates and highlights the importance of including both measures in biodiversity assessments. Massive sequencing of soil and litter samples holds the potential to complement traditional biological inventories in advancing our understanding of the factors affecting tropical diversity.

Constructing networks has become an indispensable approach in understanding how different taxa interact. However, methodologies vary widely among studies, potentially limiting our ability to meaningfully compare results. In particular, how network architecture is influenced by the extent to which nodes are resolved to either taxa or taxonomic units is poorly understood. To address this, here we collate nine datasets of ecological interactions, from both observations and DNA metabarcoding, and construct networks under a range of commonly-used node resolutions. We demonstrate that small changes in node resolution can cause wide variation in almost all key metric values, including robustness and nestedness. Moreover, relative values of metrics such as robustness were seen to fluctuate continuously with node resolution, thereby potentially confounding comparisons of networks, as well as interpretations concerning their constituent ecological interactions. These findings highlight the need for care when comparing networks, especially where these differ with respect to node resolution.




Tuesday, June 26, 2018

Three opportunities to learn about metabarcoding

In case you are looking for ways to learn about metabarcoding, there are actually three different courses offered this year. All of them differ in approach, content focus, and venue but what they have in common is the fact that they provide participants with a comprehensive package that enables them make informed decisions when it comes to organizing experiments, field work and analytics.

This program will provide an overview of the state of current technology and the various platforms used. The course consists of a series of online lectures and research exercises introducing different aspects of metabarcoding and environmental DNA research. We will also touch on the suite of bioinformatics tools available for sequence analysis and data interpretation.

This course will focus on eDNA metabarcoding, however targeted single species detection and other alternatives will also be explored, as they can sometimes be suitable metabarcoding alternatives.

The lectures will cover different aspects of DNA metabarcoding. The bioinformatics practicals will introduce data analysis from raw sequences to basic ecological conclusions. The molecular ecology practical will present basic techniques for DNA extraction in the field and DNA amplification by PCR.

Monday, June 25, 2018

DNA barcoding for pollen forecasting

PollerGEN is a group of interdisciplinary researchers funded by NERC to understand grass pollen deposition. We aim to revolutionise the way that pollen is measured, model spatial and temporal deposition from different species of grass pollen and identify linkages to human health.

DNA barcoding of pollen is not a new invention. It is not easy either but has been shown to provide extremely valuable information, e.g. for understanding plant-pollinator interactionshoney bee foraging, or the characterization of honeybee pollen pellets. It should come to no surprise that researchers are also working on an application that intends to improve a forecasting system that has become more and more important for a large portion of the human population - pollen forecasting for hay fever and other allergic reactions.

At this point most forecasts are build using data from a network of pollen traps which operate throughout the main pollen seasons. These traps measure how many pollen grains are present on a daily basis and identifications of species are done using morphology-based methods. The latter is extremely challenging when it comes to species with very uniform appearance, e.g. grasses. However, the species identity often makes a big difference. It is fairly rare that somebody is allergic to all grass pollen but we are having difficulties to tell which pollen in the mix is the culprit.

PollerGen, a project run out of Bangor University wants to change this by using a DNA-based approach. 

The colleagues are now working on a way to detect airborne pollen from different species of allergenic grass. We’re also developing new pollen source maps, and modelling how pollen grains likely move across landscapes, as well as identifying which species are linked with the exacerbation of asthma and hay fever.

We’re going to be using a new UK plant DNA barcode library, as well as environmental genomic technologies to identify complex mixtures of tree and grass pollens from a molecular genetic perspective. By combining this information with detailed source maps and aerobiological modelling, we hope to redefine how pollen forecasts are measured and reported in the future.

We have just started the third year of pollen collection and hope to road test the combined forecasting methods over the next year. In the long run, our vision is to be able to provide specific pollen forecasts for grass, and unravel which species of grass pollen are most likely causing allergic responses. More broadly, we also want to provide information to healthcare professionals and charities, who can translate this information to help pollen allergy sufferers live healthier and more productive lives.

Pretty cool.

Friday, June 22, 2018

Weekend reads

Here we go again, another week has passed quickly. Light on posting, mainly because I had some days off and no chance to do digging for blog posts. Nevertheless, here your weekly dose of interesting papers. Really good stuff.

Genetic taxonomic assignment can be more sensitive than morphological taxonomic assignment, particularly for small, cryptic or rare species. Sequence processing is essential to taxonomic assignment, but can also produce errors because optimal parameters are not known a priori. Here, we explored how sequence processing parameters influence taxonomic assignment of 18S sequences from bulk zooplankton samples produced by 454 pyrosequencing. We optimized a sequence processing pipeline for two common research goals, estimation of species richness and early detection of aquatic invasive species (AIS), and then tested most optimal models' performances through simulations. We tested 1,050 parameter sets on 18S sequences from 20 AIS to determine optimal parameters for each research goal. We tested optimized pipelines' performances (detectability and sensitivity) by computationally inoculating sequences of 20 AIS into ten bulk zooplankton samples from ports across Canada. We found that optimal parameter selection generally depends on the research goal. However, regardless of research goal, we found that metazoan 18S sequences produced by 454 pyrosequencing should be trimmed to 375-400 bp and sequence quality filtering should be relaxed (1.5 ≤ maximum expected error ≤ 3.0, Phred score = 10). Clustering and denoising were only viable for estimating species richness, because these processing steps made some species undetectable at low sequence abundances which would not be useful for early detection of AIS. With parameter sets optimized for early detection of AIS, 90% of AIS were detected with fewer than 11 target sequences, regardless of whether clustering or denoising was used. Despite developments in next-generation sequencing, sequence processing remains an important issue owing to difficulties in balancing false-positive and false-negative errors in metabarcoding data.

DNA metabarcoding has been introduced as a revolutionary way to identify organisms and monitor ecosystems. However, the potential of this approach for biomonitoring remains partially unfulfilled because a significant part of the sampled DNA cannot be affiliated to species due to incomplete reference libraries. Thus, biotic indices which are based on the estimated abundances of species in a community and their ecological profiles can be inaccurate. We propose to compute biotic indices using phylogenetic imputation of OTUs' ecological profiles (OTU-PITI approach). Firstly, OTUs sequences are inserted within a reference phylogeny. Secondly, OTUs' ecological profiles are estimated on the basis of their phylogenetic relationships with reference species whose ecology is known. Based on these ecological profiles, biotic indices can be computed using all available OTUs. Using freshwater diatoms as a case study, we show that short DNA barcodes can be placed accurately within a phylogeny and their ecological preferences estimated with a satisfactory level of precision. In light of these results, we tested the approach with a dataset of 139 environmental samples of benthic river diatoms for which the same biotic index (IPS) was calculated using (i) traditional microscopy, (ii) OTUs with taxonomic assignment approach, (iii) OTUs with phylogenetic estimation of ecological profiles (OTU-PITI), and (iv) OTU with taxonomic assignment completed by the phylogenetic approach (OTU-PITI) for unclassified OTUs. Using traditional microscopy as a reference, we found that the combination of the OTUs' taxonomic assignment completed by the phylogenetic method performed satisfactorily and substantially better than the other methods tested.

BACKGROUND: High throughput DNA sequencing of bulk invertebrate samples or metabarcoding is becoming increasingly used to provide profiles of biological communities for environmental monitoring. As metabarcoding becomes more widely applied, new reference DNA barcodes linked to individual specimens identified by taxonomists are needed. This can be achieved through using DNA extraction methods that are not only suitable for metabarcoding but also for building reference DNA barcode libraries.
METHODS: In this study, we test the suitability of a rapid non-destructive DNA extraction method for metabarcoding of freshwater invertebrate samples.
RESULTS: This method resulted in detection of taxa from many taxonomic groups, comparable to results obtained with two other tissue-based extraction methods. Most taxa could also be successfully used for subsequent individual-based DNA barcoding and taxonomic identification. The method was successfully applied to field-collected invertebrate samples stored for taxonomic studies in 70% ethanol at room temperature, a commonly used storage method for freshwater samples.
DISCUSSION: With further refinement and testing, non-destructive extraction has the potential to rapidly characterise species biodiversity in invertebrate samples, while preserving specimens for taxonomic investigation.

Marine plankton populate 70% of Earth's surface, providing the energy that fuels ocean food webs and contributing to global biogeochemical cycles. Plankton communities are extremely diverse and geographically variable, and are overwhelmingly composed of low-abundance species. The role of this rare biosphere and its ecological underpinnings are however still unclear. Here, we analyse the extensive dataset generated by the Tara Oceans expedition for marine microbial eukaryotes (protists) and use an adaptive algorithm to explore how metabarcoding-based abundance distributions vary across plankton communities in the global ocean. We show that the decay in abundance of non-dominant operational taxonomic units, which comprise over 99% of local richness, is commonly governed by a power-law. Despite the high spatial turnover in species composition, the power-law exponent varies by less than 10% across locations and shows no biogeographical signature, but is weakly modulated by cell size. Such striking regularity suggests that the assembly of plankton communities in the dynamic and highly variable ocean environment is governed by large-scale ubiquitous processes. Understanding their origin and impact on plankton ecology will be important for evaluating the resilience of marine biodiversity in a changing ocean.

MOTIVATION: Correct taxonomic identification of DNA sequences is central to studies of biodiversity using both shotgun metagenomic and metabarcoding approaches. However, no genetic marker gives sufficient performance across all the biological kingdoms, hampering studies of taxonomic diversity in many groups of organisms. This has led to the adoption of a range of genetic markers for DNA metabarcoding. While many taxonomic classification software tools can be re-trained on these genetic markers, they are often designed with assumptions that impair their utility on genes other than the SSU and LSU rRNA. Here, we present an update to Metaxa2 that enables the use of any genetic marker for taxonomic classification of metagenome and amplicon sequence data.
RESULTS: We evaluated the Metaxa2 Database Builder on eleven commonly used barcoding regions and found that while there are wide differences in performance between different genetic markers, our software performs satisfactorily provided that the input taxonomy and sequence data are of high quality.
AVAILABILITY: Freely available on the web as part of the Metaxa2 package at http://microbiology.se/software/metaxa2/.

BACKGROUND: The world's herbaria contain millions of specimens, collected and named by thousands of researchers, over hundreds of years. However, this treasure has remained largely inaccessible to genetic studies, because of both generally limited success of DNA extraction and the challenges associated with PCR-amplifying highly degraded DNA. In today's next-generation sequencing world, opportunities and prospects for historical DNA have changed dramatically, as most NGS methods are actually designed for taking short fragmented DNA molecules as templates.
RESULTS: As a practical test of routine recovery of rDNA and plastid genome sequences from herbarium specimens, we sequenced 25 herbarium specimens up to 80 years old from 16 different Angiosperm families. Paired-end reads were generated, yielding successful plastid genome assemblies for 23 species and nuclear rDNAs for 24 species, respectively. These data showed that genome skimming can be used to generate genomic information from herbarium specimens as old as 80 years and using as little as 500 pg of degraded starting DNA.
CONCLUSIONS: The routine plastome sequencing from herbarium specimens is feasible and cost-effective (compare with Sanger sequencing or plastome-enrichment approaches), and can be performed with limited sample destruction.

The increasing popularity of cytochrome c oxidase subunit 1 (COI) DNA metabarcoding warrants a careful look at the underlying reference databases used to make high-throughput taxonomic assignments. The objectives of this study are to document trends and assess the future usability of COI records for metabarcode identification. Over 2.5 million COI sequences were found in GenBank, half of which were fully identified to the species rank. From 2003 to 2017, the number of COI Eukaryote records deposited has grown by two orders of magnitude representing a nearly 42-fold increase in unique species. For fully identified records, 92% are at least 500 bp in length, 74% have a country annotation, and 51% have latitude-longitude annotations. To ensure the future usability of COI records in GenBank we suggest: 1) Improving the geographic representation of COI records 2) Improving the cross-referencing of COI records in the Barcode of Life Data System and GenBank to facilitate consolidation and incorporation into existing bioinformatic pipelines, 3) Adherence to the minimum information about a marker gene sequence guidelines, and 4) Integrating metabarcodes from eDNA and mixed community studies with existing sequences. COI metabarcoders are normally considered consumers of taxonomic data. Here we discuss the potential for taxonomists to reverse this pattern and instead mine metabarcode data to guide species discovery. The growth of COI reference records over the past 15 years has been substantial and is likely to be a resource across many fields for years to come.

Thursday, June 14, 2018

#BadStockPhotosOfMyJob

Ever seen anything in relation to the hashtag #BadStockPhotosOfMyJob? If not you should check out Twitter or search for it on Google because it really shows some ridiculously funny photos that exhibit some of the worst stereotypes people have when thinking about other's jobs. Especially the perception of what scientists do is almost tragic. I thought its a good idea to show a few examples including ironic comments by the real scientists. It's funny indeed but sometimes also just sad to see what others think we scientists do for a living.


 


  I have no words for those four.

Wednesday, June 13, 2018

Interview with a vampire

In this study, we show for the first time that it is possible to use DNA meta-barcoding to generate data on both diet and the predator's population structure. And we more or less get this additional information for free because the vampire bat's DNA is found in the DNA that we extract from blood meal and faecal samples

When the sun sets in South and Central America, the vampire bats wake up and fly out in search of prey. The vampire bat's diet consists of blood. It prefers to feed on domestic animals such as cows and pigs, but when it does so, there is a risk of transmitting pathogens such as rabies. In order to control rabies transmitted by vampire bats, it is crucial to have a method that allows large-scale assessment of vampire bat prey. A study published back in April led by researchers from Denmark and the UK, shows that metabarcoding can do just that.

The colleagues analysed vampire bat blood meal and faecal samples collected in Peru, along the coast, in the Andes and in the Amazon. In diet studies, the metabarcoding is normally only used to assess diet, but in this study, the researchers went one step further and gathered information on the vampire bat's population structure. The latter is an approach very similar to work my group has been doing in collaboration with researchers in Germany. This 'free of charge' data can help researchers understand how the landscape influences the connectivity of vampire bat populations, which could influence the spread of pathogens. 

We are slowly beginning to understand that all the metabarcoding data we generate to better understand community composition of a given environment contains several layers of information. It is perhaps much richer than an OTU table. That being said it is an entire different story on how to release let alone disentangle all that information.

It is great to gain insight into both predator and prey from DNA in droppings and blood meals. Apart from feeding on domestic animals, vampire bats occasionally took blood from wild tapirs, so the method may be useful for determining the distribution of elusive mammal prey. It is also of note that we found no evidence of vampire bats feeding on humans from the DNA left over from their dinners.





Tuesday, June 12, 2018

Citizen science vs giant slugs

Citizen science is a powerful tool to combat the challenges created by invasive species. Our study emphasizes the importance of collaborations between researchers, government administration, and citizen volunteers. 

The giant slug Limax maximus is an invasive species which made its way from northern Europe all the way to Japan and other regions of the world. It is a notorious pest of horticultural and agricultural crops. 

Recently a Japanese research team found that a certain set of weather conditions could be a reliable short-term indicator of how often giant slugs would appear on a set mountain path. The findings showed that the slugs were more likely to appear on days with higher humidity, lower windspeed and lower precipitation than the 20-year average. These observations can be used to predict future  outbreaks of the pest. 

This study was actually made possible by citizen science. In order to survey the number of slugs present on the mountain path chosen for the study (Mt. Maruyama route, in Sapporo, Japan), a volunteer naturalist hiked the path at 5:00 AM nearly every day for two years. The colleagues collected weather data obtained from a nearby meteorological station and combined them with observational data to calculate correlations between slug appearances and complex weather conditions.