Initial Upper Palaeolithic humans in Europe had recent Neanderthal ancestry.

Modern humans appeared in Europe by at least 45,000 years ago1-5, but the extent of their interactions with Neanderthals, who disappeared by about 40,000 years ago6, and their relationship to the broader expansion of modern humans outside Africa are poorly understood. Here we present genome-wide data from three individuals dated to between 45,930 and 42,580 years ago from Bacho Kiro Cave, Bulgaria1,2. They are the earliest Late Pleistocene modern humans known to have been recovered in Europe so far, and were found in association with an Initial Upper Palaeolithic artefact assemblage. Unlike two previously studied individuals of similar ages from Romania7 and Siberia8 who did not contribute detectably to later populations, these individuals are more closely related to present-day and ancient populations in East Asia and the Americas than to later west Eurasian populations. This indicates that they belonged to a modern human migration into Europe that was not previously known from the genetic record, and provides evidence that there was at least some continuity between the earliest modern humans in Europe and later people in Eurasia. Moreover, we find that all three individuals had Neanderthal ancestors a few generations back in their family history, confirming that the first European modern humans mixed with Neanderthals and suggesting that such mixing could have been common.

Weill-Marchesani syndrome: natural history and genotype-phenotype correlations from 18 news cases and review of literature.

BACKGROUND: Weill-Marchesani syndrome (WMS) belongs to the group of acromelic dysplasias, defined by short stature, brachydactyly and joint limitations. WMS is characterised by specific ophthalmological abnormalities, although cardiovascular defects have also been reported. Monoallelic variations in FBN1 are associated with a dominant form of WMS, while biallelic variations in ADAMTS10, ADAMTS17 and LTBP2 are responsible for a recessive form of WMS. OBJECTIVE: Natural history description of WMS and genotype-phenotype correlation establishment. MATERIALS AND METHODS: Retrospective multicentre study and literature review. INCLUSION CRITERIA: clinical diagnosis of WMS with identified pathogenic variants. RESULTS: 61 patients were included: 18 individuals from our cohort and 43 patients from literature. 21 had variants in ADAMTS17, 19 in FBN1, 19 in ADAMTS10 and 2 in LTBP2. All individuals presented with eye anomalies, mainly spherophakia (42/61) and ectopia lentis (39/61). Short stature was present in 73% (from -2.2 to -5.5 SD), 10/61 individuals had valvulopathy. Regarding FBN1 variants, patients with a variant located in transforming growth factor (TGF)-ß-binding protein-like domain 5 (TB5) domain were significantly smaller than patients with FBN1 variant outside TB5 domain (p=0.0040). CONCLUSION: Apart from the ophthalmological findings, which are mandatory for the diagnosis, the phenotype of WMS seems to be more variable than initially described, partially explained by genotype-phenotype correlation.

Environmental Drivers of the Moonmilk Microbiome Diversity in Some Temperate and Tropical Caves.

Moonmilk is a cave deposit that was used for medical and cosmetic purposes and has lately raised interest for its antimicrobial potential. We studied five moonmilk samples from four caves with different microclimatic conditions, two temperate in north-western and northern Romania (Ferice, Fața Apei, and Izvorul Taușoarelor caves) and one tropical in Minas Gerais, Brazil (Nestor Cave). The physicochemical and mineralogical analyses confirmed the presence of calcite and dolomite as the main phase in the moonmilk. A 16S rRNA gene-based metabarcoding approach showed the most abundant bacteria phyla Proteobacteria, GAL15, Actinobacteriota, and Acidobacteriota. The investigated caves differed in the dominant orders of bacteria, with the highest distance between the Romanian and Nestor Cave samples. Climate and, implicitly, the soil microbiome can be responsible for some differences we found between all the samples. However, other factors can be involved in shaping the moonmilk microbiome, as differences were found between samples in the same cave (Ferice). In our five moonmilk samples, 1 phylum, 70 orders (~ 36%), and 252 genera (~ 47%) were unclassified, which hints at the great potential of cave microorganisms for future uses.

The genetic history of Ice Age Europe.

Modern humans arrived in Europe ~45,000 years ago, but little is known about their genetic composition before the start of farming ~8,500 years ago. Here we analyse genome-wide data from 51 Eurasians from ~45,000-7,000 years ago. Over this time, the proportion of Neanderthal DNA decreased from 3-6% to around 2%, consistent with natural selection against Neanderthal variants in modern humans. Whereas there is no evidence of the earliest modern humans in Europe contributing to the genetic composition of present-day Europeans, all individuals between ~37,000 and ~14,000 years ago descended from a single founder population which forms part of the ancestry of present-day Europeans. An ~35,000-year-old individual from northwest Europe represents an early branch of this founder population which was then displaced across a broad region, before reappearing in southwest Europe at the height of the last Ice Age ~19,000 years ago. During the major warming period after ~14,000 years ago, a genetic component related to present-day Near Easterners became widespread in Europe. These results document how population turnover and migration have been recurring themes of European prehistory.

Combining exome/genome sequencing with data repository analysis reveals novel gene-disease associations for a wide range of genetic disorders.

PURPOSE: Within this study, we aimed to discover novel gene-disease associations in patients with no genetic diagnosis after exome/genome sequencing (ES/GS). METHODS: We followed two approaches: (1) a patient-centered approach, which after routine diagnostic analysis systematically interrogates variants in genes not yet associated to human diseases; and (2) a gene variant centered approach. For the latter, we focused on de novo variants in patients that presented with neurodevelopmental delay (NDD) and/or intellectual disability (ID), which are the most common reasons for genetic testing referrals. Gene-disease association was assessed using our data repository that combines ES/GS data and Human Phenotype Ontology terms from over 33,000 patients. RESULTS: We propose six novel gene-disease associations based on 38 patients with variants in the BLOC1S1, IPO8, MMP15, PLK1, RAP1GDS1, and ZNF699 genes. Furthermore, our results support causality of 31 additional candidate genes that had little published evidence and no registered OMIM phenotype (56 patients). The phenotypes included syndromic/nonsyndromic NDD/ID, oral-facial-digital syndrome, cardiomyopathies, malformation syndrome, short stature, skeletal dysplasia, and ciliary dyskinesia. CONCLUSION: Our results demonstrate the value of data repositories which combine clinical and genetic data for discovering and confirming gene-disease associations. Genetic laboratories should be encouraged to pursue such analyses for the benefit of undiagnosed patients and their families.

An early modern human from Romania with a recent Neanderthal ancestor.

Neanderthals are thought to have disappeared in Europe approximately 39,000-41,000 years ago but they have contributed 1-3% of the DNA of present-day people in Eurasia. Here we analyse DNA from a 37,000-42,000-year-old modern human from Pestera cu Oase, Romania. Although the specimen contains small amounts of human DNA, we use an enrichment strategy to isolate sites that are informative about its relationship to Neanderthals and present-day humans. We find that on the order of 6-9% of the genome of the Oase individual is derived from Neanderthals, more than any other modern human sequenced to date. Three chromosomal segments of Neanderthal ancestry are over 50 centimorgans in size, indicating that this individual had a Neanderthal ancestor as recently as four to six generations back. However, the Oase individual does not share more alleles with later Europeans than with East Asians, suggesting that the Oase population did not contribute substantially to later humans in Europe.

Simultaneous Determination of As, Bi, Sb, Se, Te, Hg, Pb and Sn by Small-Sized Electrothermal Vaporization Capacitively Coupled Plasma Microtorch Optical Emission Spectrometry Using Direct Liquid Microsampling.

The simultaneous determination of chemical vapor-generating elements involving derivatization is difficult even by inductively coupled plasma optical emission spectrometry or mass spectrometry. This study proposes a new direct liquid microsampling method for the simultaneous determination of As, Bi, Se, Te, Hg, Pb, and Sn, using a fully miniaturized set-up based on electrothermal vaporization capacitively coupled plasma microtorch optical emission spectrometry. The method is cost-effective, free from non-spectral interference, and easy to run by avoiding derivatization. The method involves the vaporization of analytes from the 10 µL sample and recording of episodic spectra generated in low-power (15 W) and low-Ar consumption (150 mL min-1) plasma microtorch interfaced with low-resolution microspectrometers. Selective vaporization at 1300 °C ensured the avoidance of non-spectral effects and allowed the use of external calibration. Several spectral lines for each element even in the range 180-210 nm could be selected. Generally, this spectral range is examined with large-scale instrumentation. Even in the absence of derivatization, the obtained detection limits were low (0.02-0.75 mg kg-1) and allowed analysis of environmental samples, such as cave and river sediments. The recovery was in the range of 86-116%, and the accuracy was better than 10%. The method is of general interest and could be implemented on any miniaturized or classical laboratory spectrometric instrumentation.

Application of Inductively Coupled Plasma Spectrometric Techniques and Multivariate Statistical Analysis in the Hydrogeochemical Profiling of Caves-Case Study Cloșani, Romania.

The aim of the study was to develop the hydrogeochemical profiling of caves based on the elemental composition of water and silty soil samples and a multivariate statistical analysis. Major and trace elements, including rare earths, were determined in the water and soil samples. The general characteristics of water, anions content, inorganic and organic carbon fractions and nitrogen species (NO3- and NH4+) were also considered. The ANOVA-principal component analysis (PCA) and two-way joining analysis were applied on samples collected from Cloșani Cave, Romania. The ANOVA-PCA revealed that the hydrogeochemical characteristics of Ca2+-HCO3- water facies were described by five factors, the strongest being associated with water-carbonate rock interactions and the occurrence of Ca, Mg and HCO3- (43.4%). Although organic carbon fractions have a lower influence (20.1%) than inorganic ones on water characteristics, they are involved in the chemical processes of nitrogen and of the elements involved in redox processes (Fe, Mn, Cr and Sn). The seasonal variability of water characteristics, especially during the spring, was observed. The variability of silty soil samples was described by four principal components, the strongest influence being attributed to rare earth elements (52.2%). The ANOVA-PCA provided deeper information compared to Gibbs and Piper diagrams and the correlation analysis.

Local- versus broad-scale environmental drivers of continental ß-diversity patterns in subterranean spider communities across Europe.

Macroecologists seek to identify drivers of community turnover (ß-diversity) through broad spatial scales. However, the influence of local habitat features in driving broad-scale ß-diversity patterns remains largely untested, owing to the objective challenges of associating local-scale variables to continental-framed datasets. We examined the relative contribution of local- versus broad-scale drivers of continental ß-diversity patterns, using a uniquely suited dataset of cave-dwelling spider communities across Europe (35-70° latitude). Generalized dissimilarity modelling showed that geographical distance, mean annual temperature and size of the karst area in which caves occurred drove most of ß-diversity, with differential contributions of each factor according to the level of subterranean specialization. Highly specialized communities were mostly influenced by geographical distance, while less specialized communities were mostly driven by mean annual temperature. Conversely, local-scale habitat features turned out to be meaningless predictors of community change, which emphasizes the idea of caves as the human accessible fraction of the extended network of fissures that more properly represents the elective habitat of the subterranean fauna. To the extent that the effect of local features turned to be inconspicuous, caves emerge as experimental model systems in which to study broad biological patterns without the confounding effect of local habitat features.

Macrodactyly in tuberous sclerosis complex: Case report and review of the literature.

Macrodactyly in the context of tuberous sclerosis complex (TSC) is a known but rare manifestation. We report the case of a boy diagnosed with TSC at 2 years and 4 months of age, presenting with bilateral macrodactyly of the first three fingers of both hands, with underlying radiographic changes, in whom molecular analysis identified a frameshift mutation on the TSC1 gene (encoding hamartin), leading to a premature stop codon. We reviewed the literature for reported cases of TSC patients with the same manifestation. In four of 14 patients, including ours, macrodactyly caused some type of joint limitation or flexion deformity, thus contradicting the established idea that this is a finding without clinical significance. Our patient is, to our knowledge, the first reported to have clear bilateral involvement. We briefly discuss the underlying mechanism for this phenomenon, which has yet to be fully elucidated, although somatic mosaicism for loss of heterozygosity at TSC loci is a plausible explanation. © 2016 Wiley Periodicals, Inc.

Ancient human footprints in Ciur-Izbuc Cave, Romania.

In 1965, Ciur-Izbuc Cave in the Carpathian Mountains of Romania was discovered to contain about 400 ancient human footprints. At that time, researchers interpreted the footprints to be those of a man, woman and child who entered the cave by an opening which is now blocked but which was usable in antiquity. The age of the prints (≈10-15 ka BP) was based partly on their association with cave bear (Ursus spelaeus) footprints and bones, and the belief that cave bears became extinct near the end of the last ice age. Since their discovery, the human and bear evidence and the cave itself have attracted spelunkers and other tourists, with the result that the ancient footprints are in danger of destruction by modern humans. In an effort to conserve the footprints and information about them and to reanalyze them with modern techiques, Ciur-Izbuc Cave was restudied in summer of 2012. Modern results are based on fewer than 25% of the originally described human footprints, the rest having been destroyed. It is impossible to confirm some of the original conclusions. The footprints do not cluster about three different sizes, and the number of individuals is estimated to be six or seven. Two cases of bears apparently overprinting humans help establish antiquity, and C-14 dates suggest a much greater age than originally thought. Unfortunately, insufficient footprints remain to measure movement variables such as stride length. However, detailed three-dimensional mapping of the footprints does allow a more precise description of human movements within the cave.

A familial case of diffuse periventricular nodular heterotopia identified prenatally: Filamin A defect as the probable cause.

Periventricular nodular heterotopia (PNH) is a neuronal migration defect characterized by the presence of ectopic grey matter nodules adjacent to the walls of the lateral ventricles. The main genetic etiology of PNH are variants in the Filamin A gene (FLNA, MIM #300049), located in the X chromosome. It affects mostly females (embryonic lethality in males), with about 50% of cases inherited from healthy mothers or with a mild phenotype. It is associated with epilepsy (75%-90%), cardiovascular (65%) and pulmonary pathologies (25%). A 28-year-old primigravida was referred for prenatal care in obstetrics department because of personal history of obliterative bronchiolitis. She has a family history of asthma (mother and sister) and adulthood-onset epilepsy (father). The pregnancy was uneventful up to 20 weeks and 3 days when bilateral periventricular irregularities and mega cisterna magna were identified on ultrasound in a female fetus. Neurosonography was performed, which led to the hypothesis of diffuse PNH, supported by MRI. The hypothesis of PNH associated to the FLNA gene was made. Brain MRI on the pregnant woman was requested, which confirmed a similar pattern of PNH. The arrayCGH (PerkinElmer, Prenatal filter 37K) was normal, and whole exome sequencing identified the likely pathogenic c.1554del p.(Val519fs*) variant in the FLNA gene. We present a case of X-linked hereditary PNH that highlights the value of fetal neurosonography in making a putative diagnosis. The diagnosis was supported by MRI in both fetus and mother. The investigation was supplemented by genetic studies, which confirmed the diagnosis.

A collaborative backbone resource for comparative studies of subterranean evolution: The World Asellidae database.

Transition to novel environments, such as groundwater colonization by surface organisms, provides an excellent research ground to study phenotypic evolution. However, interspecific comparative studies on evolution to groundwater life are few because of the challenge in assembling large ecological and molecular resources for species-rich taxa comprised of surface and subterranean species. Here, we make available to the scientific community an operational set of working tools and resources for the Asellidae, a family of freshwater isopods containing hundreds of surface and subterranean species. First, we release the World Asellidae database (WAD) and its web application, a sustainable and FAIR solution to producing and sharing data and biological material. WAD provides access to thousands of species occurrences, specimens, DNA extracts and DNA sequences with rich metadata ensuring full scientific traceability. Second, we perform a large-scale dated phylogenetic reconstruction of Asellidae to support phylogenetic comparative analyses. Of 424 terminal branches, we identify 34 pairs of surface and subterranean species representing independent replicates of the transition from surface water to groundwater. Third, we exemplify the usefulness of WAD for documenting phenotypic shifts associated with colonization of subterranean habitats. We provide the first phylogenetically controlled evidence that body size of males decreases relative to that of females upon groundwater colonization, suggesting competition for rare receptive females selects for smaller, more agile males in groundwater. By making these tools and resources widely accessible, we open up new opportunities for exploring how phenotypic traits evolve in response to changes in selective pressures and trade-offs during groundwater colonization.

Role of TBX20 Truncating Variants in Dilated Cardiomyopathy and Left Ventricular Noncompaction.

BACKGROUND: Less than 40% of patients with dilated cardiomyopathy (DCM) have a pathogenic/likely pathogenic genetic variant identified. TBX20 has been linked to congenital heart defects; although an association with left ventricular noncompaction (LVNC) and DCM has been proposed, it is still considered a gene with limited evidence for these phenotypes. This study sought to investigate the association between the TBX20 truncating variant (TBX20tv) and DCM/LVNC. METHODS: TBX20 was sequenced by next-generation sequencing in 7463 unrelated probands with a diagnosis of DCM or LVNC, 22 773 probands of an internal comparison group (hypertrophic cardiomyopathy, channelopathies, or aortic diseases), and 124 098 external controls (individuals from the gnomAD database). Enrichment of TBX20tv in DCM/LVNC was calculated, cosegregation was determined in selected families, and clinical characteristics and outcomes were analyzed in carriers. RESULTS: TBX20tv was enriched in DCM/LVNC (24/7463; 0.32%) compared with internal (1/22 773; 0.004%) and external comparison groups (4/124 098; 0.003%), with odds ratios of 73.23 (95% CI, 9.90-541.45; P<0.0001) and 99.76 (95% CI, 34.60-287.62; P<0.0001), respectively. TBX20tv was cosegregated with DCM/LVNC phenotype in 21 families for a combined logarythm of the odds score of 4.53 (strong linkage). Among 57 individuals with TBX20tv (49.1% men; mean age, 35.9±20.8 years), 41 (71.9%) exhibited DCM/LVNC, of whom 14 (34.1%) had also congenital heart defects. After a median follow-up of 6.9 (95% CI, 25-75:3.6-14.5) years, 9.7% of patients with DCM/LVNC had end-stage heart failure events and 4.8% experienced malignant ventricular arrhythmias. CONCLUSIONS: TBX20tv is associated with DCM/LVNC; congenital heart defect is also present in around one-third of cases. TBX20tv-associated DCM/LVNC is characterized by a nonaggressive phenotype, with a low incidence of major cardiovascular events. TBX20 should be considered a definitive gene for DCM and LVNC and routinely included in genetic testing panels for these phenotypes.

Oribatid mite fossils from pre-Quaternary sediments in Slovenian caves II. Amiracarus pliocennatus n.gen., n.sp. (Microzetidae) from Pliocene, with comments on the other species of the genus.

Amiracarus pliocennatus n. gen., n. sp. is described based on fossils found in Pliocene and Pleistocene clastic sediments in caves of the Slovenian Classical Karst (Kras). Diagnosis and relationships of the new genus are given and relationships within the extant species of the genus are discussed as well as variability and reliability of some characters, used for species differentiation. Amiracarus senensis (Bernini 1975) n. comb. is proposed as a type species of the new genus, and other four extant species are newly placed in Amiracarus n. gen.: A. abeloosi (Lions 1978) n. comb., A. discrepans (Mahunka 1966) n. comb., A. grootaerti (Wauthy &Ducarme 2011) n. comb. and A. similis (Subìas & Iturrondobeitia 1978) n. comb. Finding of a fossil individual of M. senensis Bernini 1975 in Ponicova Cave of South-West Carpathians in Romania is reported and a description of this individual is given. Relationships with the genus Miracarus Kunst 1959 are discussed and a new diagnosis of this genus is also given, along with the redescription of its type species, M. hurkai Kunst 1959.

Karst spring microbiome: Diversity, core taxa, and community response to pathogens and antibiotic resistance gene contamination.

Karst aquifers are important water resources for drinking water supplies worldwide. Although they are susceptible to anthropogenic contamination due to their high permeability, there is a lack of detailed knowledge on the stable core microbiome and how contamination may affect these communities. In this study, eight karst springs (distributed across three different regions in Romania) were sampled seasonally for one year. The core microbiota was analysed by 16S rRNA gene amplicon sequencing. To identify bacteria carrying antibiotic resistance genes and mobile genetic elements, an innovative method was applied, consisting of high-throughput antibiotic resistance gene quantification performed on potential pathogen colonies cultivated on Compact Dry™ plates. A taxonomically stable bacterial community consisting of Pseudomonadota, Bacteroidota, and Actinomycetota was revealed. Core analysis reaffirmed these results and revealed primarily freshwater-dwelling, psychrophilic/psychrotolerant species affiliated to Rhodoferax, Flavobacterium, and Pseudomonas genera. Both sequencing and cultivation methods indicated that more than half of the springs were contaminated with faecal bacteria and pathogens. These samples contained high levels of sulfonamide, macrolide, lincosamide and streptogramins B, and trimethoprim resistance genes spread primarily by transposase and insertion sequences. Differential abundance analysis found Synergistota, Mycoplasmatota, and Chlamydiota as suitable candidates for pollution monitoring in karst springs. This is the first study highlighting the applicability of a combined approach based on high-throughput SmartChip™ antibiotic resistance gene quantification and Compact Dry™ pathogen cultivation for estimating microbial contaminants in karst springs and other challenging low biomass environments.

Subterranean transport of microplastics as evidenced in karst springs and their characterization using Raman spectroscopy.

The increasing use of plastic materials has led to accumulation of large amounts of plastic waste in environment and a global challenge to be tackled with. The natural process of macro-plastics aging generates a multitude of secondary microplastic fragments accumulating in all areas of the planet. The pollution with microplastics of large water bodies, such as rivers, seas and oceans was already proven, but the presence of microplastics even in karst spring water was not reported yet. In this study, Raman micro-spectroscopy was used to confirm the presence of microplastics in the spring water samples collected from two rural karst springs in the Apuseni Mountains (Țarina and Josani), North-Western Romania. Two sets of water samples of 1000 L collected in spring time 2021 and one in autumn 2021 were filtered and analyzed. Using the Python programming language and combining two separate Raman databases, one for plastics and the other for pigments, we established a customized database to unambiguously identify the type of plastic and pigment present in the discovered micro-fragments. The generated reference pigment-plastic spectra were compared to those of potential microplastics found on filters and Pearson's coefficient was used to measure the level of similarity. The presence of microplastics in karst spring waters was confirmed and a quantitative estimation expressed as number of fragments or fibers per liter was 0.034 in Josani and 0.06 in Țarina karst spring. Five months later sampling (autumn 2021) revealed 0.05 microplastics per liter. The spectral results revealed that most microplastics found were dominated by polyethylene terephthalate (PET), followed by polypropylene and interestingly, abundant blue micro-fragments were identified according to their copper phthalocyanine pigments (pigment Blue 15) or indigo carmine (pigment Blue 63) characteristic spectral fingerprints, which surpassed the inherent spectral background level characteristic for the Raman spectra of naturally contaminated waste micro-samples. Their origin in mountain karst spring waters and potential decrease in time is discussed.

Characterization of three adults and an adolescent with Osteogenesis Imperfecta type VI and a novel founder SERPINF1 variant.

Osteogenesis imperfecta (OI) type VI is an extremely rare form of OI caused by biallelic variants in the SERPINF1 gene, which codes for the pigment-epithelium derived factor (PEDF). We report on four patients (three adults and one adolescent) with a severe deforming form of OI. All patients presented no abnormalities at birth, frequent long bone and vertebrae fractures (mainly during childhood), marked short stature, severe bone deformities, chronic mild to moderate pain, and severe limitation of mobility, with three being completely wheelchair bound. Blue sclera and dentinogenesis imperfecta were absent, although some patients presented tooth, ophthalmological, and/or cardiac features. Radiographic findings included, among others, thin diaphysis and popcorn calcifications, both of which are non-specific to this type of OI. The novel homozygous variants c.816_819del (p.Met272Ilefs*8) and c.283+2T > G in SERPINF1 were identified in three and one patient, respectively. The three patients carrying the frameshift variant were born in nearby regions suggesting a founder effect. Describing the long-term outcomes of four patients with OI type VI, this cohort adds relevant data on the clinical features and prognosis of this type of OI.

Diversity, distribution and organic substrates preferences of microbial communities of a low anthropic activity cave in North-Western Romania.

Introduction: Karst caves are characterized by relatively constant temperature, lack of light, high humidity, and low nutrients availability. The diversity and functionality of the microorganisms dwelling in caves micro-habitats are yet underexplored. Therefore, in-depth investigations of these ecosystems aid in enlarging our understanding of the microbial interactions and microbially driven biogeochemical cycles. Here, we aimed at evaluating the diversity, abundance, distribution, and organic substrate preferences of microbial communities from Peștera cu Apa din Valea Leșului (Leșu Cave) located in the Apuseni Mountains (North-Western Romania). Materials and Methods: To achieve this goal, we employed 16S rRNA gene amplicon sequencing and community-level physiological profiling (CLPP) paralleled by the assessment of environmental parameters of cave sediments and water. Results and Discussion: Pseudomonadota (synonym Proteobacteria) was the most prevalent phylum detected across all samples whereas the abundance detected at order level varied among sites and between water and sediment samples. Despite the general similarity at the phylum-level in Leșu Cave across the sampled area, the results obtained in this study suggest that specific sites drive bacterial community at the order-level, perhaps sustaining the enrichment of unique bacterial populations due to microenvironmental conditions. For most of the dominant orders the distribution pattern showed a positive correlation with C-sources such as putrescine, γ-amino butyric acid, and D-malic acid, while particular cases were positively correlated with polymers (Tween 40, Tween 80 and α-cyclodextrin), carbohydrates (α-D-lactose, i-erythritol, D-mannitol) and most of the carboxylic and ketonic acids. Physicochemical analysis reveals that sediments are geochemically distinct, with increased concentration of Ca, Fe, Al, Mg, Na and K, whereas water showed low nitrate concentration. Our PCA indicated the clustering of different dominant orders with Mg, As, P, Fe, and Cr. This information serves as a starting point for further studies in elucidating the links between the taxonomic and functional diversity of subterranean microbial communities.

The gut microbiome mediates adaptation to scarce food in Coleoptera.

Beetles are ubiquitous cave invertebrates worldwide that adapted to scarce subterranean resources when they colonized caves. Here, we investigated the potential role of gut microbiota in the adaptation of beetles to caves from different climatic regions of the Carpathians. The beetles' microbiota was host-specific, reflecting phylogenetic and nutritional adaptation. The microbial community structure further resolved conspecific beetles by caves suggesting microbiota-host coevolution and influences by local environmental factors. The detritivore species hosted a variety of bacteria known to decompose and ferment organic matter, suggesting turnover and host cooperative digestion of the sedimentary microbiota and allochthonous-derived nutrients. The cave Carabidae, with strong mandibula, adapted to predation and scavenging of animal and plant remains, had distinct microbiota dominated by symbiotic lineages Spiroplasma or Wolbachia. All beetles had relatively high levels of fermentative Carnobacterium and Vagococcus involved in lipid accumulation and a reduction of metabolic activity, and both features characterize adaptation to caves.