Bacterial communities of the oviduct of turkeys | Scientific Reports - Nature.com
Abstract
Bacterial communities in the reproductive tract of avian species play an important role in keeping birds healthy and encouraging growth. Infection can occur during egg formation with pathogens that can be transmitted to the embryo. In this study, we investigated the bacterial composition in the turkey reproductive tract using a taxa identification based on the amplicon sequence of the V3–V4 region of the 16S rRNA gene. The microbial composition and relative abundance of bacteria differed between individual birds. Among the 19 phyla detected in turkey oviduct were unique taxa like Planctomycetes or Petescibacteria. Differences in composition of bacterial diversity were found at the family and genus level. Oviducts contained also several genus with well-recognized avian pathogens like Escherichia-Shigella, Enterococcus, Staphylococcus, and Ornithobacterium. Some of the bacteria described in this study have not been so far identified in turkeys. The objective of this study was to identify bacterial communities in the turkey oviduct and compared the composition of the oviduct with that in chickens broadening the knowledge of the microbial composition in the reproductive tract of poultry.
Introduction
The microbial community found throughout the organism and in systems important for physical development play an important role in the health and correct growth of birds1,2,3. The role of the avian oviduct is to secrete the components surrounding the yolk, provide the anatomical structure where fertilization takes place, and transport the egg through the reproductive tract4. In the avian reproductive tract, the bacterial community can be associated with the reproductive organs and also be sexually transmitted or opportunistic5. These bacteria can affect the host's defense, immunity, and future colonization with other microbes. The reproductive tract may also contain environmental contaminants that have spread via copulation6,7.
Little is known about the bacterial composition of the female turkey reproductive tract. Many pathogens associated with it can lead to the development of clinical signs which will result in significant economic losses8,9. Infections with some bacteria cause poorer weight gain, internal and external egg quality, and hatchability rate as well as greater mortality and susceptibility to opportunistic bacterial infections10,11. The occurrence of diseases in the reproductive tract is also associated with a drop in egg production12,13. In some cases, opportunistic bacteria have the potential to become pathogenic when other factors act upon the immune system of the host8,14. Additionally the potential for vertical transmission of pathogens like Salmonella, Campylobacter species, or Escherichia coli in the oviduct by transferring these pathogens to the embryo through the egg white makes the maintenance of a healthy oviduct microbiome very important8,14.
Considering the possibility of pathogen transmission both through the yolk and on the surface of the eggshell, it is important to know the composition of the bacterial community in the oviduct of birds. The limited knowledge so far is based on studies describing the reproductive system of chickens14,15,16,17. Addressing the lack of research on the bacterial composition of the reproductive tract in poultry species other than chickens, the main objective of the present study was to determine and characterization the microbial community in the laying turkey reproductive system. The bacterial composition of turkey oviducts was also compared with that of chickens. In addition, we investigated the influence of different farming conditions on the composition of the microbiome by evaluating samples from two different farms from birds located in the same geographical region kept under similar environmental conditions.
Knowledge of the bacterial composition in the turkey oviduct can be useful in controlling, diagnosing, and treating commercial turkey flocks. The microbial community in the reproductive tract can influence the physiology of the oviduct, laying, and egg quality. To the best of our knowledge, this is the first study describing the bacterial composition in the laying turkey oviduct.
Results
In this study, we characterized the bacterial composition of fifteen turkey oviducts. Samples T-OVI-1 to T-OVI-5 were from turkeys at 49 weeks of age from farm A and samples T-OVI-6 to T-OVI-15 were from turkeys at 52 weeks of age from farm B. Among the 19 phyla detected in oviducts, the dominant phyla was Proteobacteria (35.07% ± 30.48%). In descending order Firmicutes were the next phyla (33.05% ± 26.72%) followed by Actinobacteria and Bacteroidetes with an abundance of 14.82% ± 11.56% and 12.24% ± 23.72% respectively (Fig. 1a). Phyla Acidobacteriota, Planctomycetes, Verrucomicrobia, Fusobacteria, Desulfobactera, Campylobacteriota, Gemmatimonadetes, Myxococcotawere and also six phyla with abundance lower than 0.1 (Fig. 1a and Supplementary Table 1).
At the order level the turkey oviduct is more diverse and is dominated by Lactobacillales 17.50% ± 22.08%, Enterobacteriales 14.6% ± 31.1%, Pseudomonadales 11.71% ± 21.02%, Flavobacteriales 6.4% ± 23.52%, Bacteroidales 5.3% ± 7.89%, Burkholderiales 4.15% ± 5.88% and Micrococcales 3.97% ± 4.53% (Fig. 1b).
At the family level, the most common taxa was Enterobacteriaceae (14.58% ± 31.11%). Next there were the Moraxeceae (11.2% ± 21.17), Lactobacillaceae (7.3%% ± 15.54%), Weeksellaceae (6.16% ± 23.56%), Enterococcaceae (6.13% ± 14.1%), Lachnospiraceae (3.66% ± 4.93%), and Streptococcaceae (3.4% ± 5.53%) families. Their relative abundances were highly variable between individual birds. The Shannon index, showed large individual variance of oviduct samples. To estimate the richness of bacterial species that were identified in the samples the Chao index was used (Fig. 2a, b). In the PCoA plot, the most oviduct samples from turkey from farm B were separated from the oviduct from farm A. No significant differences in these two farms were noted (Kruskal–Wallis test; p < 0.05) (Fig. 2c). Venn diagrams were used to compare the microbial composition of turkey oviducts at different levels on two different farms. At the family level, farm A and B have 113 shared taxa out of the 117 found in birds on farm A and 181 on farm B. At genus level from 236 taxa in turkeys oviducts from farm A and 761 taxa in oviduct from farm B 192 taxa were shared (Fig. 3).
The Escherichia-Shigella (14.38% ± 31.2%), Enhydrobacter (9.97% ± 21.54%), Lactobacillus (7.3% ± 15.54%), Enterococcus (6.13% ± 14.1%), Chryseobacrium (6.12% ± 23.57%), and Streptococcus (3.32% ± 5.55%) were the most abundant genus in the turkey oviducts. Many genus ranged below 1% relative abundance (Fig. 4b and Supplementary Table 3).
Discussion
Sequencing of the 16S rRNA gene is widely used to study the genomic diversity of bacterial communities and the microbiome in both healthy and various diseased states15,18. The bacterial composition of the gut and respiratory tract of poultry has been relatively well described, whereas little is known about the microbiome of the turkey oviduct2,19,20,21,22. Pathogenic bacteria in this organ can impair the turkey's defense mechanism against microorganisms, weaken its immunity, and facilitate future colonization. Pathogens transmitted vertically can inflict significant economic losses on the poultry industry because they decrease egg production and the hatchability rate whilst increasing mortality8,14,23,24. Clinical signs in the oviduct caused by bacterial infections that are present in the host's body for a long time may be troublesome to treat because the bacteria are resistant to some antibiotics. Therefore, knowledge of the bacterial composition in the healthy turkey oviduct would enable better treatment and a reduction in the use of antibiotics in turkey flocks. The 16S high-throughput sequencing of the turkey oviduct provides an opportunity to discover the bacterial communities present in the reproductive tract of turkey hens.
In this study, we characterized the bacterial composition of fifteen turkey oviductsfrom two commercial farms. Taxa which were detected in greater than 16% of oviducts were considered to analysis. Differences in the microbial composition of individual birds were noted. Among the 19 phyla detected in oviducts, the dominant phyla were Proteobacteria, Firmicutes, Actinobacteria and Bacteroidetes (Fig. 1a and Supplementary Table 1). Slightly different results were reported in egg-laying hens in Korea and meat-type hens in Israel where the predominant type was Firmicutes14,15. Comparing the bacterial composition of turkey oviducts with chicken oviducts a small amount of some phyla were detected in both (Planctomycetes and Verrucomicrobia)15. Parts of these phyla were described in the avian respiratory tract21,25. Some of these bacteria were also found in water and soil samples (Nitrospira, Desulfobactera, and Campylobacteriota)26,27,28. There were also bacteria from superphyla such as those from the Planctomycetes or Patescibacteria29,30. A few of these were recently established bacteria such as Gemmatimonadetes and Myxococcota31,32,33,34. As far as we could ascertain, the described phyla have not yet been identified in turkeys.
To study the effect of different farming conditions on the composition of the microbiome, we examined samples from two different farms from birds located in the same geographical region kept under similar environmental conditions. The differences between the microbial composition in the oviducts of turkeys from farm A and farm B are most evident in the number of shared taxa at family and genus level (Fig. 3). Part of the differences in taxa abundance may be due to various numbers of samples. Other factors that can affect abundance are other differences in poultry farming such as the feed for turkey or feed additives used during laying. The diversity of bacteria composition was varied in the turkey oviduct. Comparing the composition of the oviduct microbiota in chickens and turkeys, it seems that although they represent different species and locations, both have a similar composition in the main groups of bacteria. Some differences were more shown at the order level. The chickens oviduct is dominated by Lactobacillales, Bacteroidales, Clostridiales, Burkholderiales, Pseudomonadales, and Bacillales. In turkeys, Lactobacillales, Pseudomonadales, and Bacillalesare at a similar level but the others are of much lower value (Table 1). At the order level the turkey oviduct is more diverse and is dominated by Lactobacillales, Enterobacteriales, Pseudomonadales, Flavobacteriales, Bacteroidales, Burkholderiales and Micrococcales.
At the family level, the most common taxa was Enterobacteriaceae. Next there were the Moraxeceae, Lactobacillaceae, Weeksellaceae, Enterococcaceae, Lachnospiraceae and Streptococcaceae families. Their relative abundances were highly variable between individual birds. Most turkey hens had a unique family pool in the oviduct (Fig. 4a and Supplementary Table 2).
The most abundant genus in the turkey oviducts were Escherichia-Shigella, Enhydrobacter, Lactobacillus, Enterococcus, Chryseobacrium and Streptococcus (Fig. 4b and Supplementary Table 3). Lactobacillus species were not seen to dominate the reproductive tract of turkeys, which is a similar finding to that of a study of the chicken oviduct15. For turkey samples, taxa representing a few species of Lactobacillus (L. aviarius, L. salivarius, unclassified Lactobacillus species, L. panis, L. plantrum, and L. ingluviei) were found. The opposite results were observed in the chicken reproductive tract, where L. salivarius was the dominant lactobacilli and L. aviarius was absent15.
The likelihood of transfer of some bacteria from the gut to the oviduct, and even the full length of the oviduct up to the infundibulum, have been investigated in chickens14. Our results indicate that this is also possible in turkeys. Low abundance of Alistipes, Bacteroides, Blautia, Butyricicoccus, Helicobacter, Phascolarctobacterium, Pseudoflavonifractor, Roseburia, Ruminococcus, Slackia, Subdoligranulum, and unclassified bacteria presence in turkey caecum were also in the oviduct (Supplementary Table 3)35. These results may confirm the conclusions that despite the environment of the oviduct containing lysozymes and other antimicrobials, it is exposed to bacteria from the cloaca or from the external environment16,36.
During formation, the egg could be contaminated with bacteria which can be vertically transferred to the embryos13,37. Some of these bacteria, as a result of existence in the oviduct tissue, may result in eggshell anomalies which reduce the quality of the eggs. In the present study, Escherichia-Shigella, and Ornithobacterium are at least two genus with pathogenic potential to be vertically transmitted24. Ornithobacterium is an emergent poultry breeding respiratory pathogen, acting opportunistically in most cases. Clinical symptoms may include, but are not limited to, reduced egg production. The presence of Ornithobacterium rhinotracheale was confirmed via the 16S rRNA gene identification by real-time PCR (data not showed). In this work, we also detected the presence of Enterococcus spp. mostly unclassified below genus level, but a small volume of these were Enterococcus faecium, E. faecalis, E. duran, and E. cecorum which most frequently are associated with clinical avian diseases38. They are also an emerging pathogen in the poultry industry. Enterococcus species was identified by matrix-assisted laser desorption/ionization–time-of-flight mass spectrometry (MALDI TOF MS) (data not showed). The bacterial species in the oviduct that may also affect avian health are Avibacterium paragallinarum, and Escherichia coli which was also identified by MALDI TOF MS (data not showed). Some of these microbial species are commonly found in both healthy and infected poultry respiratory tracts in turkeys and chickens.
Of course, not all bacteria are pathogenic. A recent study observed that oviduct microbiota in the chicken is involved in the synthesis and deposition of eggshell pigments. Vaginal Staphylococcus and Ralsonia might affect the content of pigments of the shell cuticle16.This study also observed the presence of these bacteria at very low levels.
Our preliminary study characterizes the microbial communities of the turkey oviduct, hence broadening the knowledge of the reproductive tract in poultry. Further studies based on larger sample sizes are necessary to provide an investigation on the role of microorganisms in egg formation and their good quality and effect on the host13,16. The present results emphasize that the bacterial composition of the turkey oviduct varies between individual birds. The differences between birds from farms A and B can be seen at lower bacteria taxonomic levels. The results also suggest that the bacterial composition of turkey oviducts differs from the community found in chicken oviducts. The results of this preliminary study may provide new insights for further identification of novel pathogens and interrelation with the reproduction process and laying performance in turkeys as well as provide information on the differences in microbial composition and bacterial diversity in poultry oviducts. To the best of our knowledge, this is the first study to examine the bacterial composition of the reproductive tract of turkeys.
Methods
Sample collection
At the end of the production cycle (49–52 weeks of age), fifteen turkey hens BUT-6 from two turkey breeding farms (farm A n = 5 and farm B n = 10) were delivered to the National Veterinary Research Institute in Poland as part of diagnostic tests. After euthanizing the birds by decapitation, necropsy was conducted under aseptic conditions during which samples of the oviduct were collected. Reproductive tract samples comprising the infundibulum, magnum, and uterus (middle section of each part) were aseptically collected for analysis of the microbiota. All birds were kept under similar environmental conditions on both farms geographically located in the same province with access to water ad libitum and fed regally. The turkey females were artificially inseminated using a sterile injector. At the time of sampling, the birds did not display lesions signifying any potential clinical disease.
DNA extraction and16S rRNA gene sequencing
All oviduct samples were suspended in Tris buffer (10 mM, pH 8, 5; Eurx, GdaÅ„sk, Poland) homogenized using a manual LabGEN 125 device (Cole-Parmer, Vernon Hills, IL, USA). Every part of the homogenized tissue (3 parts of oviduct: infundibulum, magnum, and uterus) in the birds were pooled (n = 15) and stored for one day at − 20 °C until DNA extraction. Before starting extraction, 50 µl lysozyme (10 mg/ml, Novazym), 6 µl mutanolysin (5KU/ml, Sigma-Aldrich), and 8 µl lysostaphin (5 g/ml, Sigma-Aldrich) were added to the samples followed by incubation for 45 min at 37 °C. DNA was extracted using a Maxwell RSC Tissue Pathogen Kit (Promega, Madison, WI, USA). DNA extraction from the Tris buffer used for sample preparation was conducted as a negative control.
After extracting DNA the V3–V4 hypervariable regions of the 16S rRNA gene were amplified using 341F and 785R primers and the library was prepared39. Sequencing of samples was performed externally as a commercial service (Genomed, Warsaw, Poland) using MiSeq paired-end 2 × 300 bp technology in a v3 kit (Illumina, San Diego, CA, USA).
Bacterial composition of turkey oviducts analysis
Raw reads (obtained after sequencing) were subjected to quality control in Cutadapt software where adapter, primer sequences, and low quality bases were removed40. Sequences, that were too short, were filtered out. Sequences were processed and taxonomy assigned using QIIME241 amplicon sequence variants (ASVs) were determined with DADA2 using the denoise-paired method. The SILVA 138 release was used as the reference database for the taxonomic assignment42,43. Further analysis was carried out using the R program, phyloseq and vegan packages44,45,46. We focused on taxa that were present in at least two oviducts in birds on farm A or B. We did not consider taxa occurring only in individual birds. From a total of 106 taxa at the order level, 185 taxa at the family level and 805 taxa at the genus level, the 20 most abundant taxa at each level were extracted and are shown in Figs. 1 and 4. The relative taxa abundance of the flocks is presented as a mean % value. Alpha diversity was measured using the Shannon and Chao index. The Kruskal–Wallis test was performed on samples from turkeys from farm A and B for statistically significant differences in alpha diversity values. The value of p < 0.05 was considered statistically significant. Beta diversity principal coordinate analysis (PCoA) plots were generated based on the Bray–Curtis method47. Venn diagrams were constructed showing the number of taxa at the phylum, order, family, and genus levels (bioinformatics.psb.ugent.be/webtools/Venn/).
Ethics declarations
The samples were collected from commercial turkeys farm by authorized veterinarians during clinical studies following standard procedures for diagnostic examination. All methods used in this study were carried out in accordance with relevant guidelines and regulations. The study was conducted in accordance with the guidelines of Animal Research: Reporting of In Vivo Experiments. According to the local law and Directive 2010/63/EU on the protection of animals used for scientific purposes (Chapter I, article 1, p. 5 b, d; article 3 p.1) the formal ethical approval is not required for this kind of study.
Data availability
The sequences were deposited in the NCBI Sequence Read Archive under BioProject accession Number: PRJNA810919. The datasets generated during the current study are available from the corresponding author on reasonable request.
References
Grond, K., Sandercock, B. K., Jumpponen, A. & Zeglin, L. H. The avian gut microbiota: Community, physiology and function in wild birds. J. Avian Biol. 49, 1–19 (2018).
Stanley, D., Hughes, R. J. & Moore, R. J. Microbiota of the chicken gastrointestinal tract: Influence on health, productivity and disease. Appl. Microbiol. Biotechnol. 98, 4301–4310 (2014).
Wei, S., Morrison, M. & Yu, Z. Bacterial census of poultry intestinal microbiome. Poult. Sci. 92, 671–683 (2013).
Bakst, M. R. & Akuffo, V. Morphology of the turkey vagina with and without an egg mass in the uterus. Poult. Sci. 88, 631–635 (2009).
Hunton, P. Research on eggshell structure and quality: An historical overview. Rev. Bras. Ciência AvÃcola 7, 67–71 (2005).
Bellinvia, S., Johnston, P. R., Reinhardt, K. & Otti, O. Bacterial communities of the reproductive organs of virgin and mated common bedbugs, cimex lectularius. Ecol. Entomol. 45, 142–154 (2020).
Kers, J. G. et al. Host and environmental factors affecting the intestinal microbiota in chickens. Front. Microbiol. 9, 1–14 (2018).
Gast, R. K., Guraya, R., Guard-Bouldin, J., Holt, P. S. & Moore, R. W. Colonization of specific regions of the reproductive tract and deposition at different locations inside eggs laid by hens infected with salmonella enteritidis or salmonella heidelberg. Avian Dis. 51, 40–44 (2007).
Srinivasan, P., Balasubramaniam, G. A., Gopala Krishna Murthy, T. R. & Balachandran, P. Prevalence and pathology of oviduct impaction in commercial white leghorn layer chicken in Namakkal region of India. Vet. World 7, 553–558 (2014).
Turan, N. & Ak, S. Investigation of the presence of ornithobacterium rhinotracheale in chickens in turkey and determination of the seroprevalance of the infection using the enzyme-linked immunosorbent assay. Avian Dis. 46, 442–446 (2002).
Landman, W. J. M. Is Mycoplasma synoviae outrunning mycoplasma gallisepticum? A viewpoint from the Netherlands. Avian Pathol. 43, 2–8 (2014).
Roque, L. & Soares, M. C. Effects of eggshell quality and broiler breeder age on hatchability. Poult. Sci. 73, 1838–1845 (1994).
Kursa, O., PakuÅ‚a, A., Tomczyk, G., PaÅ›ko, S. & Sawicka, A. Eggshell apex abnormalities caused by two different mycoplasma synoviae genotypes and evaluation of eggshell anomalies by full-field optical coherence tomography. BMC Vet. Res. 15, 1–8 (2019).
Shterzer, N. et al. Large overlap between the intestinal and reproductive tract microbiomes of chickens. Front. Microbiol. 11, 1–9 (2020).
Lee, S. et al. Characterization of microbial communities in the chicken oviduct and the origin of chicken embryo gut microbiota. Sci. Rep. 9, 1–11 (2019).
Wen, C. et al. Microbiota continuum along the chicken oviduct and its association with host genetics and egg formation. Poult. Sci. 100, 101104 (2021).
Lee, S. J. et al. Comparison of microbiota in the cloaca, colon, and magnum of layer chicken. PLoS ONE 15, 1–15 (2020).
Moon, C. D., Cookson, A. L., Young, W., Maclean, P. H. & Bermingham, E. N. Metagenomic insights into the roles of proteobacteria in the gastrointestinal microbiomes of healthy dogs and cats. Microbiol. Open https://doi.org/10.1002/mbo3.677 (2018).
Wei, S., Lilburn, M. & Yu, Z. The Bacteriomes of Ileal mucosa and cecal content of broiler chickens and turkeys as revealed by metagenomic analysis. Int. J. Microbiol. 2016, 1–12 (2016).
Johnson, T. J. et al. A consistent and predictable commercial broiler chicken bacterial microbiota in antibiotic-free production displays strong correlations with performance. Appl. Environ. Microbiol. 84, 1–18 (2018).
Shabbir, M. Z. et al. Microbial communities present in the lower respiratory tract of clinically healthy birds in Pakistan. Poult. Sci. 94, 612–620 (2014).
Konieczka, P. et al. Increased dietary inclusion levels of lysine are more effective than arginine in supporting the functional status of the gut in growing turkeys. Animal 11, 1–15 (2021).
Kursa, O., Paku, A., Tomczyk, G. & Sawicka, A. Eggshell apex abnormalities caused by two different mycoplasma synoviae genotypes and evaluation of eggshell anomalies by full-field optical coherence tomography. BMC Vet. Res. 15(1), 1–8 (2019).
Barbosa, E. V. et al. Ornithobacterium rhinotracheale: An update review about an emerging poultry pathogen. Vet. Sci. 7(1), 3 (2020).
Kursa, O., Tomczyk, G., Sawicka-Durkalec, A., Giza, A. & SÅ‚omiany-Szwarc, M. Bacterial communities of the upper respiratory tract of turkeys. Sci. Rep. 11(1), 1–11 (2021).
Daims, H. & Wagner, M. Nitrospira. Trends Microbiol. https://doi.org/10.1016/j.tim.2018.02.001 (2018).
Wang, C. et al. Study on the microbial characteristics of seawater pond based on enteromorpha bloom. Haiyang Xuebao 39, 107–116 (2017).
Kalashnikova, V. A., Dzhikidze, E. K., Stasilevich, Z. K., Krylova, R. I. & Kebu, T. I. Campilobacteria in the etiology of acute intestinal infections in primates. Vestn. Ross. Akad. Meditsinskikh Nauk 1, 6–10 (2006).
Tian, R. et al. Small and mighty: Adaptation of superphylum patescibacteria to groundwater environment drives their genome simplicity. Microbiome 8, 1–15 (2020).
Ward, N. et al. The order planctomycetales, Including the genera planctomyces, pirellula, gemmata and isosphaera and the candidatus genera brocadia, kuenenia and scalindua. The Prokaryotes https://doi.org/10.1007/0-387-30747-8_31 (2006).
Uzun, M., Alekseeva, L., Krutkina, M., Koziaeva, V. & Grouzdev, D. Unravelling the diversity of magnetotactic bacteria through analysis of open genomic databases. Sci. Data 7, 1–13 (2020).
Zeng, Y. et al. Gemmatimonas groenlandica sp. nov. Is an aerobic anoxygenic phototroph in the phylum gemmatimonadetes. Front. Microbiol. 11, 1–18 (2021).
Viljakainen, V. R. & Hug, L. A. The phylogenetic and global distribution of bacterial polyhydroxyalkanoate bioplastic-degrading genes. Environ. Microbiol. 23, 1717–1731 (2021).
Sun, R. et al. Long-term N fertilization decreased diversity and altered the composition of soil bacterial and archaeal communities. Agronomy https://doi.org/10.3390/agronomy9100574 (2019).
Wilkinson, T. J. et al. Characterization of the microbiome along the gastrointestinal tract of growing turkeys. Front. Microbiol. 8, 1–11 (2017).
Panheleux, M. et al. Organic matrix composition and ultrastructure of eggshell: A comparative study. Br. Poult. Sci. 40, 240–252 (1999).
Wang, C., Pors, S. E., Olsen, R. H. & Bojesen, A. M. Transmission and pathogenicity of gallibacterium anatis and escherichia coli in embryonated eggs. Vet. Microbiol. 217, 76–81 (2018).
Landman, W. J. M., Feberwee, A., Veldman, K. T. & Mevius, D. J. Study on the vertical transmission of arthropathic and amyloidogenic enterococcus faecalis in a flock of brown layer chickens. Vet. Q. 23, 88–91 (2001).
Klindworth, A. et al. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 41, 1–11 (2013).
Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. https://doi.org/10.14806/ej.17.1.200 (2011).