A DNA barcode reference library of Croatian mosquitoes (Diptera: Culicidae): implications for identification and delimitation of species, with notes on the distribution of potential vector species.

BACKGROUND: Mosquitoes pose a risk to human health worldwide, and correct species identification and detection of cryptic species are the most important keys for surveillance and control of mosquito vectors. In addition to traditional identification based on morphology, DNA barcoding has recently been widely used as a complementary tool for reliable identification of mosquito species. The main objective of this study was to create a reference DNA barcode library for the Croatian mosquito fauna, which should contribute to more accurate and faster identification of species, including cryptic species, and recognition of relevant vector species. METHODS: Sampling was carried out in three biogeographical regions of Croatia over six years (2017-2022). The mosquitoes were morphologically identified; molecular identification was based on the standard barcoding region of the mitochondrial COI gene and the nuclear ITS2 region, the latter to identify species within the Anopheles maculipennis complex. The BIN-RESL algorithm assigned the COI sequences to the corresponding BINs (Barcode Index Number clusters) in BOLD, i.e. to putative MOTUs (Molecular Operational Taxonomic Units). The bPTP and ASAP species delimitation methods were applied to the genus datasets in order to verify/confirm the assignment of specimens to specific MOTUs. RESULTS: A total of 405 mosquito specimens belonging to six genera and 30 morphospecies were collected and processed. Species delimitation methods assigned the samples to 31 (BIN-RESL), 30 (bPTP) and 28 (ASAP) MOTUs, with most delimited MOTUs matching the morphological identification. Some species of the genera Culex, Aedes and Anopheles were assigned to the same MOTUs, especially species that are difficult to distinguish morphologically and/or represent species complexes. In total, COI barcode sequences for 34 mosquito species and ITS2 sequences for three species of the genus Anopheles were added to the mosquito sequence database for Croatia, including one individual from the Intrudens Group, which represents a new record for the Croatian mosquito fauna. CONCLUSION: We present the results of the first comprehensive study combining morphological and molecular identification of most mosquito species present in Croatia, including several invasive and vector species. With the exception of some closely related species, this study confirmed that DNA barcoding based on COI provides a reliable basis for the identification of mosquito species in Croatia.

Malassezia pachydermatis from brown bear: A comprehensive analysis reveals novel genotypes and distribution of all detected variants in domestic and wild animals.

Malassezia pachydermatis (phylum Basidiomycota, class Malasseziomycetes) is a zoophilic opportunistic pathogen with recognized potential for invasive infections in humans. Although this pathogenic yeast is widespread in nature, it has been primarily studied in domestic animals, so available data on its genotypes in the wild are limited. In this study, 80 yeast isolates recovered from 42 brown bears (Ursus arctos) were identified as M. pachydermatis by a culture-based approach. MALDI-TOF mass spectrometry (MS) was used to endorse conventional identification. The majority of samples exhibited a high score fluctuation, with 42.5% of isolates generating the best scores in the range confident only for genus identification. However, the use of young biomass significantly improved the identification of M. pachydermatis at the species confidence level (98.8%). Importantly, the same MALDI-TOF MS efficiency would be achieved regardless of colony age if the cut-off value was lowered to ≥1.7. Genotyping of LSU, ITS1, CHS2, and ß-tubulin markers identified four distinct genotypes in M. pachydermatis isolates. The most prevalent among them was the genotype previously found in dogs, indicating its transmission potential and adaptation to distantly related hosts. The other three genotypes are described for the first time in this study. However, only one of the genotypes consisted of all four loci with bear-specific sequences, indicating the formation of a strain specifically adapted to brown bears. Finally, we evaluated the specificity of the spectral profiles of the detected genotypes. MALDI-TOF MS exhibited great potential to detect subtle differences between all M. pachydermatis isolates and revealed distinct spectral profiles of bear-specific genotypes.

DNA barcoding of hard ticks (Ixodidae), notes on distribution of vector species and new faunal record for Croatia.

Molecular methods are increasingly being utilized for accurate identification of ticks (Acari: Ixodidae), especially in cases of morphologically highly similar species. In this study, we performed molecular research of the tick fauna in Croatia using DNA barcoding method. Ticks were sampled in three biogeographical regions and thirteen species were recorded: Dermacentor marginatus, Dermacentor reticulatus, Haemaphysalis concinna, Haemaphysalis inermis, Haemaphysalis punctata, Hyalomma marginatum, Ixodes frontalis, Ixodes hexagonus, Ixodes kaiseri, Ixodes ricinus, Rhipicephalus bursa, Rhipicephalus sanguineus s.l. and Rhipicephalus turanicus. Ixodes kaiseri is for the first time recorded in the fauna of Croatia. Of the thirteen hard tick species analyzed in this study, pathogens from different groups (bacteria, protozoa and viruses) have been detected in eight species in Croatia so far. For the important vector species R. sanguineus s.s., new distributional data for Croatia are given. The standard COI barcoding region was amplified, and the sequences were analyzed by species delimitation methods together with the sequences of conspecific and congeneric species from the public BOLD database. Our specimens of H. punctata represent a new, genetically distinct MOTU. A brief overview of the available public DNA barcoding data for Ixodidae is presented, highlighting the need for an integrative approach for the clarification of the taxonomic status of problematic Ixodid taxa. The results provide a basis for the establishment of a molecular data platform for the Ixodidae of the Croatian fauna.

Satellite DNA-Mediated Gene Expression Regulation: Physiological and Evolutionary Implication.

Satellite DNAs are tandemly repeated sequences organized in large clusters within (peri)centromeric and/or subtelomeric heterochromatin. However, in many species, satellite DNAs are not restricted to heterochromatin but are also dispersed as short arrays within euchromatin. Such genomic organization together with transcriptional activity seems to be a prerequisite for the gene-modulatory effect of satellite DNAs which was first demonstrated in the beetle Tribolium castaneum upon heat stress. Namely, enrichment of a silent histone mark at euchromatic repeats of a major beetle satellite DNA results in epigenetic silencing of neighboring genes. In addition, human satellite III transcripts induced by heat shock contribute to genome-wide gene silencing, providing protection against stress-induced cell death. Gene silencing mediated by satellite RNA was also shown to be fundamental for the early embryonic development of the mosquito Aedes aegypti. Apart from a physiological role during embryogenesis and heat stress response, activation of satellite DNAs in terms of transcription and proliferation can have an evolutionary impact. Spreading of satellite repeats throughout euchromatin promotes the variation of epigenetic landscapes and gene expression diversity, contributing to the evolution of gene regulatory networks and to genome adaptation in fluctuating environmental conditions.

Diversity of mosquito fauna (Diptera, Culicidae) in higher-altitude regions of Croatia.

Global climate change and the accompanying rise in temperature could affect the biology and ecology of a number of vectors, including mosquitoes. High altitude areas that were previously unsuitable for the spread of mosquito vector populations could become suitable. The aim of this research was to study the distribution of mosquito species in higher altitude regions of Croatia. Samples were collected in three areas: Slavonian Mountains, Gorski Kotar, and Middle Velebit. Specimens were morphologically determined and confirmed by DNA barcoding and other genetic markers and showed the presence of 16 species belonging to six genera. The most abundant species were the Culex pipiens complex with 50% of the collected specimens. Both pipiens (Linnaeus, 1758) and molestus (Forskal, 1775) biotypes and their hybrids were identified within the complex, followed by Culex torrentium (Martini, 1925) (20.2%), Culiseta longiareolata (Macquart, 1838) (8.5%), and the invasive species Aedes japonicus (Theobald, 1901) (7.8% of the total number of collected specimens). The remaining 12 species made up 14.7% of the collected specimens. Intraspecific COI p-distances were within the standard barcoding threshold for OTUs, while interspecific genetic distances were much higher, confirming the existence of barcoding gaps. Mosquito fauna of Croatian mountains showed a moderate variety and made 30.8% of the total number of recorded mosquito species in Croatia thus far.

Evolutionary History of Alpha Satellite DNA Repeats Dispersed within Human Genome Euchromatin.

Major human alpha satellite DNA repeats are preferentially assembled within (peri)centromeric regions but are also dispersed within euchromatin in the form of clustered or short single repeat arrays. To study the evolutionary history of single euchromatic human alpha satellite repeats (ARs), we analyzed their orthologous loci across the primate genomes. The continuous insertion of euchromatic ARs throughout the evolutionary history of primates starting with the ancestors of Simiformes (45-60 Ma) and continuing up to the ancestors of Homo is revealed. Once inserted, the euchromatic ARs were stably transmitted to the descendant species, some exhibiting copy number variation, whereas their sequence divergence followed the species phylogeny. Many euchromatic ARs have sequence characteristics of (peri)centromeric alpha repeats suggesting heterochromatin as a source of dispersed euchromatic ARs. The majority of euchromatic ARs are inserted in the vicinity of other repetitive elements such as L1, Alu, and ERV or are embedded within them. Irrespective of the insertion context, each AR insertion seems to be unique and once inserted, ARs do not seem to be subsequently spread to new genomic locations. In spite of association with (retro)transposable elements, there is no indication that such elements play a role in ARs proliferation. The presence of short duplications at most of ARs insertion sites suggests site-directed recombination between homologous motifs in ARs and in the target genomic sequence, probably mediated by extrachromosomal circular DNA, as a mechanism of spreading within euchromatin.

Stenelmis koreana Satô, 1978 (Coleoptera: Elmidae): confirmed as a wide-spread species by DNA-sequencing.

Stenelmis koreana Satô, 1978 (Coleoptera: Elmidae) is here recorded for the first time from Kyrgyzstan and Western Siberia. It was hitherto thought to be confined to Korea and the Russian Far East. The identification of a specimen from Kyrgyzstan was confirmed by DNA-sequencing after comparison with two sequences of S. koreana from Korea. The COI haplotype of the Kyrgyzstan specimen has very low sequence divergence (0.53 % or 0.0053 uncorrected p-distance) with respect to the sequences of the Korean specimens, which is within the standard intraspecific sequence divergence for COI in beetles.

Dispersion Profiles and Gene Associations of Repetitive DNAs in the Euchromatin of the Beetle Tribolium castaneum.

Satellite DNAs are tandemly repeated sequences clustered within heterochromatin. However, in some cases, such as the major TCAST1 satellite DNA from the beetle Tribolium castaneum, they are found partially dispersed within euchromatin. Such organization together with transcriptional activity enables TCAST1 to modulate the activity of neighboring genes. In order to explore if other T. castaneum repetitive families have features that could provide them with a possible gene-modulatory role, we compare here the structure, organization, dispersion profiles, and transcription activity of 10 distinct TCAST repetitive families including TCAST1. The genome organization of TCAST families exhibit either satellite-like or transposon-like characteristics. In addition to heterochromatin localization, bioinformatic searches of the assembled genome have revealed dispersion of all families within euchromatin, preferentially in the form of single repeats. Dispersed TCAST repeats are mutually correlated in distribution and are grouped in distinct regions of euchromatin. The repeats are associated with genes, are enriched in introns relative to intergenic regions, and very rarely overlap exons. In spite of the different mechanisms of repeat proliferation, such as transposition and homologous recombination, all TCAST families share a similar frequency of spreading as well as dispersion and gene association profiles. Additionally, TCAST families are transcribed and their transcription is significantly activated by heat stress. A possibility that such common features of TCAST families might be related to their potential gene-modulatory role is discussed.

Satellite DNA-like elements associated with genes within euchromatin of the beetle Tribolium castaneum.

In the red flour beetle Tribolium castaneum the major TCAST satellite DNA accounts for 35% of the genome and encompasses the pericentromeric regions of all chromosomes. Because of the presence of transcriptional regulatory elements and transcriptional activity in these sequences, TCAST satellite DNAs also have been proposed to be modulators of gene expression within euchromatin. Here, we analyze the distribution of TCAST homologous repeats in T. castaneum euchromatin and study their association with genes as well as their potential gene regulatory role. We identified 68 arrays composed of TCAST-like elements distributed on all chromosomes. Based on sequence characteristics the arrays were composed of two types of TCAST-like elements. The first type consists of TCAST satellite-like elements in the form of partial monomers or tandemly arranged monomers, up to tetramers, whereas the second type consists of TCAST-like elements embedded with a complex unit that resembles a DNA transposon. TCAST-like elements were also found in the 5' untranslated region (UTR) of the CR1-3_TCa retrotransposon, and therefore retrotransposition may have contributed to their dispersion throughout the genome. No significant difference in the homogenization of dispersed TCAST-like elements was found either at the level of local arrays or chromosomes nor among different chromosomes. Of 68 TCAST-like elements, 29 were located within introns, with the remaining elements flanked by genes within a 262 to 404,270 nt range. TCAST-like elements are statistically overrepresented near genes with immunoglobulin-like domains attesting to their nonrandom distribution and a possible gene regulatory role.

Taxonomic position of Eunapius subterraneus (Porifera, Spongillidae) inferred from molecular data--a revised classification needed?

The freshwater sponge Eunapius subterraneus was described in 1984 on the basis of its morphology and unique ecological features. It inhabits caves in the Ogulin karst area as the only known stygobitic sponge, and an endangered karst species. We used three genetic markers with different evolutionary rates in phylogenetic analyses of E. subterraneus. All of the markers exclude this sponge from the genus Eunapius. Based on our results, we emphasize the need for revision of the taxonomic classification of E. subterraneus as well as the need for a thorough re-evaluation of freshwater sponge systematics.

Wide distribution of related satellite DNA families within the genus Pimelia (Tenebrionidae).

Major satellites of species in the genus Pimelia comprise large portions of their genomes and belong to seven major satellite families which all originate from a common ancestral sequence. Here we present the results of comprehensive screening of 26 Pimelia species belonging to three distinct geographic groups (Ibero-Balearic, African and Canary Islands) for the presence of different Pimelia satellite families in their genomes. Dot-blot hybridization experiments suggest that together with one dominant, highly abundant satellite family, other families are also present in genomes of the majority of examined Pimelia species, but as low-copy number repeats. The estimated abundance of these underrepresented repeats is about 4,000 copies per haploid genome. Signals of highly abundant satellite family from P. scabrosa (PSCA) in examined congeneric species, obtained after PCR amplification and Southern hybridization under high stringency conditions, corroborate sequence preservation of low-copy representatives of satellite families. PRINS localized low-copy repeats within the pericentromeric regions of all chromosomes. These results point to the existence of an extensive library of repetitive DNAs that was already present in the genome of the common ancestor of extant Pimelia taxa, and shifts the period of diversification of Pimelia satellites far in the history of this genus.

Phylogeny of pholcid spiders (Araneae: Pholcidae): combined analysis using morphology and molecules.

The spider family Pholcidae comprises a large number of mainly tropical, web-weaving spiders, and is among the most diverse and dominant spider groups in the world. The phylogeny of this family has so far been investigated exclusively using morphological data. Here, we present the first molecular data for the family analyzed in a phylogenetic context. Four different gene regions (12S rRNA, 16S rRNA, cytochrome c oxidase subunit I, 28S rRNA) and 45 morphological characters were scored for 31 pholcid and three outgroup taxa. The data were analyzed both for individual genes, combined molecular data, and molecular plus morphological data, using parsimony, maximum likelihood, and Bayesian methods. Some of the phylogenetic hypotheses obtained previously using morphology alone were also supported by our results, like the monophyly of pholcines and of the New World clade. On the other hand, some of the previous hypotheses could be discarded with some confidence (monophyly of holocnemines, the position of Priscula), and still others need further investigation (the position of holocnemines, ninetines, and Metagonia). The data obtained provide an excellent basis for future investigations of phylogenetic patterns both within the family and among spider families.