Lithium influences whole-organism metabolic rate in Drosophila subobscura.

Lithium is widely used to treat bipolar disorder. However, the efficacy and vulnerability as to its side effects are known to differ. Although the specific biochemical mechanism of action is still elusive, lithium may influence mitochondrial function, and consequently, metabolism. Lithium exposure in this study was conducted on a unique set of mito-nuclear introgression lines of Drosophila subobscura to disentangle the independent effects of mitochondrial DNA (mtDNA) against a common nuclear DNA background. The study addressed three issues: (a) whether lithium has a dose-dependent effect on whole-organism metabolic rate, (b) whether mtDNA haplotypes show divergent metabolic efficiency measured by metabolic rate to lithium exposure and (c) whether lithium influences the whole-organism metabolic rate across sexes. The results confirm that lithium influenced the whole-organism metabolic rate, showing a subtle balance between efficacy and adverse effects within a narrow dose range. In addition, lithium exposure was found to influence metabolism differently based on mtDNA haplotypes and sex. This preliminary research may have a range of biological implications for the role of mitochondrial variability in psychiatric disease and treatment by contributing to the understanding and predicting of the lithium treatment response and risk for toxic side effects.

Effects of thyroid hormone on mitochondria and metabolism of human preimplantation embryos.

Thyroid hormones are regarded as the major controllers of metabolic rate and oxygen consumption in mammals. Although it has been demonstrated that thyroid hormone supplementation improves bovine embryo development in vitro, the cellular mechanisms underlying these effects are so far unknown. In this study, we investigated the role of thyroid hormone in development of human preimplantation embryos. Embryos were cultured in the presence or absence of 10-7 M triiodothyronine (T3) till blastocyst stage. Inner cell mass (ICM) and trophectoderm (TE) were separated mechanically and subjected to RNAseq or quantification of mitochondrial DNA copy number. Analyses were performed using DESeq (v1.16.0 on R v3.1.3), MeV4.9 and MitoMiner 4.0v2018 JUN platforms. We found that the exposure of human preimplantation embryos to T3 had a profound impact on nuclear gene transcription only in the cells of ICM (1178 regulated genes-10.5% of 11 196 expressed genes) and almost no effect on cells of TE (38 regulated genes-0.3% of expressed genes). The analyses suggest that T3 induces in ICM a shift in ribosome and oxidative phosphorylation activity, as the upregulated genes are contributing to the composition and organization of the respiratory chain and associated cofactors involved in mitoribosome assembly and stability. Furthermore, a number of genes affecting the citric acid cycle energy production have reduced expression. Our findings might explain why thyroid disorders in women have been associated with reduced fertility and adverse pregnancy outcome. Our data also raise a possibility that supplementation of culture media with T3 may improve outcomes for women undergoing in vitro fertilization.

Canalization and developmental stability of the yellow-necked mouse (Apodemus flavicollis) mandible and cranium related to age and nematode parasitism.

BACKGROUND: Mammalian mandible and cranium are well-established model systems for studying canalization and developmental stability (DS) as two elements of developmental homeostasis. Nematode infections are usually acquired in early life and increase in intensity with age, while canalization and DS of rodent skulls could vary through late postnatal ontogeny. We aimed to estimate magnitudes and describe patterns of mandibular and cranial canalization and DS related to age and parasite intensity (diversity) in adult yellow-necked mice (Apodemus flavicollis). RESULTS: We found the absence of age-related changes in the levels of canalization for mandibular and cranial size and DS for mandibular size. However, individual measures of mandibular and cranial shape variance increased, while individual measures of mandibular shape fluctuating asymmetry (FA) decreased with age. We detected mandibular and cranial shape changes during postnatal ontogeny, but revealed no age-related dynamics of their covariance structure among and within individuals. Categories regarding parasitism differed in the level of canalization for cranial size and the level of DS for cranial shape. We observed differences in age-related dynamics of the level of canalization between non-parasitized and parasitized animals, as well as between yellow-necked mice parasitized by different number of nematode species. Likewise, individual measures of mandibular and cranial shape FA decreased with age for the mandible in the less parasitized category and increased for the cranium in the most parasitized category. CONCLUSIONS: Our age-related results partly agree with previous findings. However, no rodent study so far has explored age-related changes in the magnitude of FA for mandibular size or mandibular and cranial FA covariance structure. This is the first study dealing with the nematode parasitism-related canalization and DS in rodents. We showed that nematode parasitism does not affect mandibular and cranial shape variation and covariance structure among and within individuals. However, parasite intensity (diversity) is related to ontogenetic dynamics of the levels of canalization and DS. Overall, additional studies on animals from natural populations are required before drawing some general conclusions.

The prevalence of pathogenic forms of Leptospira in natural populations of small wild mammals in Serbia.

The greatest epidemiological significance of leptospirosis in Europe comes from the fact that it is the most widespread zoonosis in the world. However, epizootiological data, especially information on maintenance hosts such as small wild mammals, are largely missing. To fill this gap in data in Serbia, we used RT-PCR for the detection of pathogenic Leptospira species and analysed 107 animals belonging to six species of small wild mammals (Apodemus agrarius, Apodemus flavicollis, Microtus arvalis, Myodes glareolus, Microtus subterraneus and Sorex araneus) collected from two localities. The animals from the first locality that was situated in a tourist area, were collected for four consecutive years (2014-2017). We found persistent incidence of infection from year to year ranging from 6.67% to 78.57%. The average frequency of infected animals was 33.3% with the highest frequency in 2014, the year characterised by a very high number of flooding days. All animals proved to be infected with pathogenic Leptospira species that were collected from the second locality situated in an agricultural area in a single year, 2014. The findings show a variable but constant presence of pathogenic Leptospira species in populations of small wild mammals in the studied areas, which indicates the need for constant monitoring.

Effect of childhood general traumas on suicide attempt depends on TPH2 and ADARB1 variants in psychiatric patients.

Suicidal behavior has been associated with a deficient serotonin neurotransmission which is likely a consequence of individual genetic architecture, exposure to environmental factors and interactions of those factors. We examined whether the interaction of child abuse, TPH2 (tryptophan hydroxylase 2) variant rs4290270, affecting alternative splicing and editing of TPH2 pre-mRNAs, and ADARB1 (adenosine deaminase acting on RNA B1) variants rs4819035 and rs9983925 may influence the risk for suicide attempt in psychiatric patients. TPH2 rs4290270 was genotyped in 165 suicide attempters and 188 suicide non-attempters diagnosed with major depressive disorder, bipolar disorder and schizophrenia. Genotyping data for ADARB1 variants were taken over from our previous study. Child abuse before the age of 18 years was assessed using the Early Trauma Inventory-Self Report. Generalized linear models and backward selection were applied to identify the main and interacting effects of environmental and genetic factors, including psychiatric diagnoses, patients' gender and age as covariates. Childhood general traumas were independently associated with suicide attempt. Two-way interaction between TPH2 rs4290270 and general traumas revealed that TT homozygotes with a history of general traumas had an increased risk for suicide attempt. Three-way interaction of general traumas, TPH2 rs4290270 and ADARB1 rs4819035 indicated that the highest predisposition to suicide attempt was observed in individuals who experienced general traumas and were TT homozygote for rs4290270 and TT homozygote for rs4819035. Our findings suggest that the risk for suicide attempt in psychiatric patients exposed to an adverse childhood environment may depend on TPH2 and ADARB1 variants.

Taking identity-by-descent analysis into the wild: Estimating realized relatedness in free-ranging macaques.

Biological relatedness is a key consideration in studies of behavior, population structure, and trait evolution. Except for parent-offspring dyads, pedigrees capture relatedness imperfectly. The number and length of DNA segments that are identical-by-descent (IBD) yield the most precise estimates of relatedness. Here, we leverage novel methods for estimating locus-specific IBD from low coverage whole genome resequencing data to demonstrate the feasibility and value of resolving fine-scaled gradients of relatedness in free-living animals. Using primarily 4-6× coverage data from a rhesus macaque (Macaca mulatta) population with available long-term pedigree data, we show that we can call the number and length of IBD segments across the genome with high accuracy even at 0.5× coverage. The resulting estimates demonstrate substantial variation in genetic relatedness within kin classes, leading to overlapping distributions between kin classes. They identify cryptic genetic relatives that are not represented in the pedigree and reveal elevated recombination rates in females relative to males, which allows us to discriminate maternal and paternal kin using genotype data alone. Our findings represent a breakthrough in the ability to understand the predictors and consequences of genetic relatedness in natural populations, contributing to our understanding of a fundamental component of population structure in the wild.

Response prediction in patients with gastric and esophagogastric adenocarcinoma under neoadjuvant chemotherapy using targeted gene expression analysis and next-generation sequencing in pre-therapeutic biopsies.

OBJECTIVES: Perioperative chemo-(radio-) therapy is the accepted standard in European patients with locally advanced adenocarcinoma of the esophagogastric junction or stomach (AEG/AS). However, 30-85% of patients do not respond to this treatment. The aim of our study was the identification of predictive biomarkers in pre-therapeutic endoscopic tumor biopsies from patients with histopathologic response (Becker-1) versus non-response (Becker-2/3) to preoperative chemotherapy. METHODS: Formalin-fixed paraffin-embedded biopsies from 36 Caucasian patients (Becker-1 n = 11, Becker-2 n = 7, Becker-3 n = 18) with AEG/AS, taken prior to neoadjuvant chemotherapy were selected. For RNA expression analysis, we employed the NanoString nCounter System. To identify genomic alterations like single nucleotide variants (SNV), copy number variation (CNV) and fusion events, we used Illumina TST170 gene panel. For HER2 and FGFR2 protein expression, immunostaining was performed. Furthermore, we analyzed the microsatellite instability (MSI) and Epstein-Barr virus (EBV) infection status by EBER in situ hybridization. RESULTS: Heat map and principal component analyses showed no clustering by means of gene expression according to regression grade. Concerning two recently proposed predictive markers, our data showed equal distribution for MSI (Becker-1: 2; Becker-2: 1; Becker-3: 3; out of 29 tested) and EBV infection was rare (1/32). We could not reveal discriminating target genes concerning SNV, but found a higher mutational burden in non-responders versus responders and fusion (in 6/14) and CNV events (in 5/14) exclusively in Becker-3. CONCLUSIONS: Although we could not identify discriminating target genes, our data suggest that molecular alterations are in general more prevalent in patients with AEG/AS belonging to the non-responding Becker group 3.

Positive Selection in Gene Regulatory Factors Suggests Adaptive Pleiotropic Changes During Human Evolution.

Gene regulatory factors (GRFs), such as transcription factors, co-factors and histone-modifying enzymes, play many important roles in modifying gene expression in biological processes. They have also been proposed to underlie speciation and adaptation. To investigate potential contributions of GRFs to primate evolution, we analyzed GRF genes in 27 publicly available primate genomes. Genes coding for zinc finger (ZNF) proteins, especially ZNFs with a Krüppel-associated box (KRAB) domain were the most abundant TFs in all genomes. Gene numbers per TF family differed between all species. To detect signs of positive selection in GRF genes we investigated more than 3,000 human GRFs with their more than 70,000 orthologs in 26 non-human primates. We implemented two independent tests for positive selection, the branch-site-model of the PAML suite and aBSREL of the HyPhy suite, focusing on the human and great ape branch. Our workflow included rigorous procedures to reduce the number of false positives: excluding distantly similar orthologs, manual corrections of alignments, and considering only genes and sites detected by both tests for positive selection. Furthermore, we verified the candidate sites for selection by investigating their variation within human and non-human great ape population data. In order to approximately assign a date to positively selected sites in the human lineage, we analyzed archaic human genomes. Our work revealed with high confidence five GRFs that have been positively selected on the human lineage and one GRF that has been positively selected on the great ape lineage. These GRFs are scattered on different chromosomes and have been previously linked to diverse functions. For some of them a role in speciation and/or adaptation can be proposed based on the expression pattern or association with human diseases, but it seems that they all contributed independently to human evolution. Four of the positively selected GRFs are KRAB-ZNF proteins, that induce changes in target genes co-expression and/or through arms race with transposable elements. Since each positively selected GRF contains several sites with evidence for positive selection, we suggest that these GRFs participated pleiotropically to phenotypic adaptations in humans.

Simulium reptans (Linnaeus, 1758) and Simulium reptantoides Carlsson, 1962 from the Balkan Peninsula.

Simulium reptans (Linnaeus, 1758) and Simulium reptantoides Carlsson, 1962 are two species of the Simulium reptans group whose distribution is unclear because of their confusing taxonomy and systematics. Their genetic variability is well known for populations in northern and central Europe and shows that both species have two forms; however, the genetic variability of these species in southern and eastern Europe is unknown. To identify the status of these two species in southeast Europe, mtDNA was extracted from 19 individuals from 12 localities across the Balkan Peninsula. Phylogenetic analysis confirmed the existence of two species with 7.38-7.94% divergence. Each species was comprised of two clades, with 2.31% and 1.43% interclade divergence for S. reptans and S. reptantoides, respectively. This study revealed the presence of both species across the Balkans and that S. reptans occurs in this area in only one form (S. reptans B), while S. reptantoides is found in two genetic forms (A and B).

Population genetic structure of the Mediterranean horseshoe bat Rhinolophus euryale in the central Balkans.

Migratory behaviour, sociality and roost selection have a great impact on the population structure of one species. Many bat species live in groups, and movements between summer and hibernation sites are common in temperate bats. The Mediterranean horseshoe bat Rhinolophus euryale is a cave-dwelling species that exhibits roost philopatry and undertakes seasonal movements which are usually shorter than 50 km. Its distribution in Serbia is restricted to karstic areas in western and eastern parts of the country, with a lack of known roosts between them. In this study, microsatellite markers were used to evaluate genetic variation in this species in the Central Balkans. Specifically, spatial genetic structuring between geographic regions and relatedness within different colony types were assessed. All analysed loci were polymorphic, and there was no significant inbreeding coefficient recorded. A moderate degree of genetic differentiation among the sampled colonies was found, and significant isolation by distance was recorded. Our results revealed that populations show a tendency to segregate into three clusters. Unexpectedly, populations from Montenegro and Eastern Serbia tended to group into one cluster, while populations from Western Serbia and Slovenia represented second and third cluster, respectively. The majority of variance was partitioned within colonies, and only a small but significant portion among clusters. Average relatedness within colony members was close to zero, did not differ significantly between the different colony types, and kinship is unlikely to be a major grouping mechanism in this species.

Phylogeography of Simulium Subgenus Wilhelmia (Diptera: Simuliidae)-Insights From Balkan Populations.

Many morphologically similar species of the simuliid (Diptera: Simuliidae) subgenus Wilhelmia, Enderlein are difficult to distinguish. Thus, the revision of the subgenus using various morphological, cytogenetic, and genetic analyses has been attempted. Neglected until now, the Balkan Peninsula, a crossroad between Europe and Anatolia, provides insight which could resolve problematic interrelationships of the taxa within this subgenus. To uncover the status and relations within the subgenus Wilhelmia, mtDNA was extracted from 47 individuals of six morphospecies: Simulium balcanicum (Enderlein, 1924), Simulium turgaicum Rubtsov, 1940, Simulium lineatum (Meigen, 1804), Simulium pseudequinum Séguy, 1921, Simulium equinum (Linnaeus, 1758), and Simulium paraequinum Puri, 1933 from 21 sites throughout the Balkan Peninsula. Phylogenetic analysis of the Wilhelmia species using mitochondrial DNA barcoding (COI) gene showed two major branches, the lineatum branch, which includes the lineages sergenti, paraequinum, and lineatum, and the equinum branch. In the equinum branch, the mtDNA sequences formed six clades, with high genetic distances, suggesting the existence of different species. Historically, the clades of the equinum branch appeared at numerous islands, perhaps as a result of allopatric speciation. The paraequinum lineage (lineatum branch) is composed of two species. However, six clades of the lineatum lineage overlapped with intra- and interspecific genetic distances. Our results revealed that the species S. balcanicum, S. pseudequinum B, and S. equinum were omnipresent in the Balkans. The results point to not only the fair diversity of Wilhelmia species in the Balkans, but also indicate that most Wilhelmia species live in sympatry.

Induction of wheat straw delignification by Trametes species.

Wheat straw is the major crop residue in European countries which makes it the most promising material for bioconversion into biofuels. However, cellulose and hemicellulose are protected with lignin, so delignification is an inevitable phase in lignocellulose processing. The organisms predominantly responsible for its degradation are white-rot fungi and among them Trametes species represent promising degraders due to a well-developed ligninolytic enzyme system. Although numerous studies have confirmed that low molecular weight compounds can induce the production and activity of ligninolytic enzymes it is not clear how this reflects on the extent of delignification. The aim of the study was to assess the capacity of p-anisidine and veratryl alcohol to induce the production and activity of Mn-oxidizing peroxidases and laccases, and wheat straw delignification by six Trametes species. Significant inter- and intraspecific variations in activity and features of these enzymes were found, as well as differences in the potential of lignocellulose degradation in the presence or absence of inducers. Differences in the catalytic properties of synthesized enzyme isoforms strongly affected lignin degradation. Apart from enhanced lignin degradation, the addition of p-anisidine could significantly improve the selectivity of wheat straw ligninolysis, which was especially evident for T. hirsuta strains.

Tryptophan Hydroxylase 1 Variant rs1800532 is Associated with Suicide Attempt in Serbian Psychiatric Patients but does not Moderate the Effect of Recent Stressful Life Events.

Tryptophan hydroxylase 1 (TPH1) gene, coding for serotonin synthesizing enzyme, and recent stressful life events (SLEs) have been commonly associated with suicidal behavior. TPH1 has been also hypothesized to be involved in stress-response mechanisms. The aim of this study was to assess TPH1 variant rs1800532 and its possible interaction with recent SLEs as risk factors for suicide attempt (SA) in Serbian psychiatric patients, including 165 suicide attempters and 188 suicide nonattempters. rs1800532 and recent SLEs were independently associated with SA, while rs1800532 did not moderate the effect of recent SLEs on SA vulnerability among Serbian psychiatric patients.

Possible influence of B chromosomes on genes included in immune response and parasite burden in Apodemus flavicollis.

Genetic background underlying wild populations immune response to different parasites is still not well understood. We studied immune response to multiple infections and to competition between different parasite species at different developmental stages in population of yellow-necked mouse, Apodemus flavicollis. Quantitative real-time PCR was used to investigate associations of MHC II-DRB, IL-10 and Tgf-ß genes expressions with presence of intestinal parasites at different developmental stages. Furthermore, we were interested whether the host related characteristics (sex, age, body condition, presence of B chromosomes or expression of other genes) or characteristics of present parasites (number of adult parasites of each identified species, egg count of each parasite genus, total number of nematode individuals) affect differential expression of the studied genes. A significant invert association between the expression of MHC II-DRB and Tgf-ß gene was found, which together with absence of IL-10 association confirmed modified Th2 as the main type of immune response to nematode infections. Effect of recorded parasites and parasite life-cycle stage on expression levels of MHC II-DRB gene was detected only through interactions with host-related characteristics such as sex, age, and the presence of B chromosomes. The presence of B chromosomes is associated with lower expression level of Tgf-ß gene. Although the influence of host genetic background on parasite infection has already been well documented, this is the first study in mammals that gave presence of B chromosomes on immune response full consideration.