Refine
Year of publication
Document Type
- Preprint (2316) (remove)
Has Fulltext
- yes (2316)
Keywords
- Kollisionen schwerer Ionen (33)
- heavy ion collisions (27)
- Deutsch (23)
- Quark-Gluon-Plasma (14)
- equation of state (13)
- QGP (12)
- heavy-ion collisions (11)
- Kongress (10)
- Syntax (10)
- quark-gluon plasma (10)
Institute
- Physik (1424)
- Frankfurt Institute for Advanced Studies (FIAS) (1006)
- Informatik (803)
- Medizin (176)
- Extern (82)
- Biowissenschaften (76)
- Ernst Strüngmann Institut (70)
- Mathematik (48)
- MPI für Hirnforschung (47)
- Psychologie (46)
Endogenous clocks enable organisms to adapt their physiology and behavior to daily variation in environmental conditions. Metabolic processes in cyanobacteria to humans are effected by the circadian clock, and its dysregulation causes metabolic disorders. In mouse and Drosophila were shown that the circadian clock directs translation of factors involved in ribosome biogenesis and synchronizes protein synthesis. However, the role of clocks in Drosophila neurogenesis and the potential impact of clock impairment on neural circuit formation and function is less understood. Here we demonstrate that light stimuli or circadian clock causes a defect in neural stem cell growth and proliferation accompanied by reduced nucleolar size. Further, we define that light and clock independently affect the InR/TOR growth regulatory pathway due to the effect on regulators of protein biosynthesis. Altogether, these data suggest that alterations in growth regulatory pathways induced by light and clock are associated with impaired neural development.
Under natural conditions, the visual system often sees a given input repeatedly. This provides an opportunity to optimize processing of the repeated stimuli. Stimulus repetition has been shown to strongly modulate neuronal-gamma band synchronization, yet crucial questions remained open. Here we used magnetoencephalography in 30 human subjects and find that gamma decreases across ~10 repetitions and then increases across further repetitions, revealing plastic changes of the activated neuronal circuits. Crucially, changes induced by one stimulus did not affect responses to other stimuli, demonstrating stimulus specificity. Changes partially persisted when the inducing stimulus was repeated after 25 minutes of intervening stimuli. They were strongest in early visual cortex and increased interareal feedforward influences. Our results suggest that early visual cortex gamma synchronization enables adaptive neuronal processing of recurring stimuli. These and previously reported changes might be due to an interaction of oscillatory dynamics with established synaptic plasticity mechanisms.
The firing pattern of ventral midbrain dopamine neurons is controlled by afferent and intrinsic activity to generate prediction error signals that are essential for reward-based learning. Given the absence of intracellular in vivo recordings in the last three decades, the subthreshold membrane potential events that cause changes in dopamine neuron firing patterns remain unknown. By establishing stable in vivo whole-cell recordings of >100 spontaneously active midbrain dopamine neurons in anaesthetized mice, we identified the repertoire of subthreshold membrane potential signatures associated with distinct in vivo firing patterns. We demonstrate that dopamine neuron in vivo activity deviates from a single spike pacemaker pattern by eliciting transient increases in firing rate generated by at least two diametrically opposing biophysical mechanisms: a transient depolarization resulting in high frequency plateau bursts associated with a reactive, depolarizing shift in action potential threshold; and a prolonged hyperpolarization preceding slower rebound bursts characterized by a predictive, hyperpolarizing shift in action potential threshold. Our findings therefore illustrate a framework for the biophysical implementation of prediction error and sensory cue coding in dopamine neurons by tuning action potential threshold dynamics.
Genome-wide CRISPR screens are becoming more widespread and allow the simultaneous interrogation of thousands of genomic regions. Although recent progress has been made in the analysis of CRISPR screens, it is still an open problem how to interpret CRISPR mutations in non-coding regions of the genome. Most of the tools concentrate on the interpretation of mutations introduced in gene coding regions. We introduce a computational pipeline that uses epigenomic information about regulatory elements for the interpretation of CRISPR mutations in non-coding regions. We illustrate our approach on the analysis of a genome-wide CRISPR screen in hTERT-RPE-1 cells and reveal novel regulatory elements that mediate chemoresistance against doxorubicin in these cells. We infer links to established and to novel chemoresistance genes. Our approach is general and can be applied on any cell type and with different CRISPR enzymes.
Several studies have probed perceptual performance at different times after a self-paced motor action and found frequency-specific modulations of perceptual performance phase-locked to the action. Such action-related modulation has been reported for various frequencies and modulation strengths. In an attempt to establish a basic effect at the population level, we had a relatively large number of participants (n=50) perform a self-paced button press followed by a detection task at threshold, and we applied both fixed- and random-effects tests. The combined data of all trials and participants surprisingly did not show any significant action-related modulation. However, based on previous studies, we explored the possibility that such modulation depends on the participant’s internal state. Indeed, when we split trials based on performance in neighboring trials, then trials in periods of low performance showed an action-related modulation at ≈17 Hz. When we split trials based on the performance in the preceding trial, we found that trials following a “miss” showed an action-related modulation at ≈17 Hz. Finally, when we split participants based on their false-alarm rate, we found that participants with no false alarms showed an action-related modulation at ≈17 Hz. All these effects were significant in random-effects tests, supporting an inference on the population. Together, these findings indicate that action-related modulations are not always detectable. However, the results suggest that specific internal states such as lower attentional engagement and/or higher decision criterion are characterized by a modulation in the beta-frequency range.
Several recent studies investigated the rhythmic nature of cognitive processes that lead to perception and behavioral report. These studies used different methods, and there has not yet been an agreement on a general standard. Here, we present a way to test and quantitatively compare these methods. We simulated behavioral data from a typical experiment and analyzed these data with several methods. We applied the main methods found in the literature, namely sine-wave fitting, the Discrete Fourier Transform (DFT) and the Least Square Spectrum (LSS). DFT and LSS can be applied both on the averaged accuracy time course and on single trials. LSS is mathematically equivalent to DFT in the case of regular, but not irregular sampling - which is more common. LSS additionally offers the possibility to take into account a weighting factor which affects the strength of the rhythm, such as arousal. Statistical inferences were done either on the investigated sample (fixed-effect) or on the population (random-effect) of simulated participants. Multiple comparisons across frequencies were corrected using False-Discovery-Rate, Bonferroni, or the Max-Based approach. To perform a quantitative comparison, we calculated Sensitivity, Specificity and D-prime of the investigated analysis methods and statistical approaches. Within the investigated parameter range, single-trial methods had higher sensitivity and D-prime than the methods based on the averaged-accuracy-time-course. This effect was further increased for a simulated rhythm of higher frequency. If an additional (observable) factor influenced detection performance, adding this factor as weight in the LSS further improved Sensitivity and D-prime. For multiple comparison correction, the Max-Based approach provided the highest Specificity and D-prime, closely followed by the Bonferroni approach. Given a fixed total amount of trials, the random-effect approach had higher D-prime when trials were distributed over a larger number of participants, even though this gave less trials per participant. Finally, we present the idea of using a dampened sinusoidal oscillator instead of a simple sinusoidal function, to further improve the fit to behavioral rhythmicity observed after a reset event.
Analyzing non-invasive recordings of electroencephalography (EEG) and magnetoencephalography (MEG) directly in sensor space, using the signal from individual sensors, is a convenient and standard way of working with this type of data. However, volume conduction introduces considerable challenges for sensor space analysis. While the general idea of signal mixing due to volume conduction in EEG/MEG is recognized, the implications have not yet been clearly exemplified. Here, we illustrate how different types of activity overlap on the level of individual sensors. We show spatial mixing in the context of alpha rhythms, which are known to have generators in different areas of the brain. Using simulations with a realistic 3D head model and lead field and data analysis of a large resting-state EEG dataset, we show that electrode signals can be differentially affected by spatial mixing by computing a sensor complexity measure. While prominent occipital alpha rhythms result in less heterogeneous spatial mixing on posterior electrodes, central electrodes show a diversity of rhythms present. This makes the individual contributions, such as the sensorimotor mu-rhythm and temporal alpha rhythms, hard to disentangle from the dominant occipital alpha. Additionally, we show how strong occipital rhythms rhythms can contribute the majority of activity to frontal channels, potentially compromising analyses that are solely conducted in sensor space. We also outline specific consequences of signal mixing for frequently used assessment of power, power ratios and connectivity profiles in basic research and for neurofeedback application. With this work, we hope to illustrate the effects of volume conduction in a concrete way, such that the provided practical illustrations may be of use to EEG researchers to in order to evaluate whether sensor space is an appropriate choice for their topic of investigation.
The human brain achieves visual object recognition through multiple stages of nonlinear transformations operating at a millisecond scale. To predict and explain these rapid transformations, computational neuroscientists employ machine learning modeling techniques. However, state-of-the-art models require massive amounts of data to properly train, and to the present day there is a lack of vast brain datasets which extensively sample the temporal dynamics of visual object recognition. Here we collected a large and rich dataset of high temporal resolution EEG responses to images of objects on a natural background. This dataset includes 10 participants, each with 82,160 trials spanning 16,740 image conditions. Through computational modeling we established the quality of this dataset in five ways. First, we trained linearizing encoding models that successfully synthesized the EEG responses to arbitrary images. Second, we correctly identified the recorded EEG data image conditions in a zero-shot fashion, using EEG synthesized responses to hundreds of thousands of candidate image conditions. Third, we show that both the high number of conditions as well as the trial repetitions of the EEG dataset contribute to the trained models’ prediction accuracy. Fourth, we built encoding models whose predictions well generalize to novel participants. Fifth, we demonstrate full end-to-end training of randomly initialized DNNs that output M/EEG responses for arbitrary input images. We release this dataset as a tool to foster research in visual neuroscience and computer vision.
In the course of global climate change, central Europe is experiencing more frequent and prolonged periods of drought. The drought years 2018 and 2019 affected European beeches (Fagus sylvatica L.) differently: even in the same stand, drought damaged trees neighboured healthy trees, suggesting that the genotype rather than the environment was responsible for this conspicuous pattern. We used this natural experiment to study the genomic basis of drought resistance with Pool-GWAS. Contrasting the extreme phenotypes identified 106 significantly associated SNPs throughout the genome. Most annotated genes with associated SNPs (>70%) were previously implicated in the drought reaction of plants. Non-synonymous substitutions led either to a functional amino acid exchange or premature termination. A SNP-assay with 70 loci allowed predicting drought phenotype in 98.6% of a validation sample of 92 trees. Drought resistance in European beech is a moderately polygenic trait that should respond well to natural selection, selective management, and breeding.
The gradual heterogeneity of climatic factors pose varying selection pressures across geographic distances that leave signatures of clinal variation in the genome. Separating signatures of clinal adaptation from signatures of other evolutionary forces, such as demographic processes, genetic drift, and adaptation to non-clinal conditions of the immediate local environment is a major challenge. Here, we examine climate adaptation in five natural populations of the harlequin fly Chironomus riparius sampled along a climatic gradient across Europe. Our study integrates experimental data, individual genome resequencing, Pool-Seq data, and population genetic modelling. Common-garden experiments revealed a positive correlation of population growth rates corresponding to the population origin along the climate gradient, suggesting thermal adaptation on the phenotypic level. Based on a population genomic analysis, we derived empirical estimates of historical demography and migration. We used an FST outlier approach to infer positive selection across the climate gradient, in combination with an environmental association analysis. In total we identified 162 candidate genes as genomic basis of climate adaptation. Enriched functions among these candidate genes involved the apoptotic process and molecular response to heat, as well as functions identified in other studies of climate adaptation in other insects. Our results show that local climate conditions impose strong selection pressures and lead to genomic adaptation despite strong gene flow. Moreover, these results imply that selection to different climatic conditions seems to converge on a functional level, at least between different insect species.
One of the major problems in evolutionary biology is to elucidate the relationships between historical events and the tempo and mode of lineage divergence. The development of relaxed molecular clock models and the increasing availability of DNA sequences resulted in more accurate estimations of taxa divergence times. However, finding the link between competing historical events and divergence is still challenging. Here we investigate assigning constrained-age priors to nodes of interest in a time-calibrated phylogeny as a means of hypothesis comparison. These priors are equivalent to historic scenarios for lineage origin. The hypothesis that best explains the data can be selected by comparing the likelihood values of the competing hypotheses, modelled with different priors. A simulation approach was taken to evaluate the performance of the prior-based method and to compare it with an unconstrained approach. We explored the effect of DNA sequence length and the temporal placement and span of competing hypotheses (i.e. historic scenarios) on selection of the correct hypothesis and the strength of the inference. Competing hypotheses were compared applying a posterior simulation analogue of the Akaike Information Criterion and Bayes factors (obtained after calculation of the marginal likelihood with three estimators: Harmonic Mean, Stepping Stone and Path Sampling). We illustrate the potential application of the prior-based method on an empirical data set to compare competing geological hypotheses explaining the biogeographic patterns in Pleurodeles newts. The correct hypothesis was selected on average 89% times. The best performance was observed with DNA sequence length of 3500-10000 bp. The prior-based method is most reliable when the hypotheses compared are not temporally too close. The strongest inferences were obtained when using the Stepping Stone and Path Sampling estimators. The prior-based approach proved effective in discriminating between competing hypotheses when used on empirical data. The unconstrained analyses performed well but it probably requires additional computational effort. Researchers applying this approach should rely only on inferences with moderate to strong support. The prior-based approach could be applied on biogeographical and phylogeographical studies where robust methods for historical inferences are still lacking.
Active transposable elements (TEs) may result in divergent genomic insertion and abundance patterns among conspecific populations. Upon secondary contact, such divergent genetic backgrounds can theoretically give rise to classical Dobzhansky-Muller incompatibilities (DMI), a way how TEs can contribute to the evolution of endogenous genetic barriers and eventually population divergence. We investigated whether differential TE activity created endogenous selection pressures among conspecific populations of the non-biting midge Chironomus riparius, focussing on a Chironomus-specific TE, the minisatellite-like Cla-element, whose activity is associated with speciation in the genus. Using an improved and annotated draft genome for a genomic study with five natural C. riparius populations, we found highly population-specific TE insertion patterns with many private insertions. A highly significant correlation of pairwise population FST from genome-wide SNPs with the FST estimated from TEs suggests drift as the major force driving TE population differentiation. However, the significantly higher Cla-element FST level due to a high proportion of differentially fixed Cla-element insertions indicates that segregating, i.e. heterozygous insertions are selected against. With reciprocal crossing experiments and fluorescent in-situ hybridisation of Cla-elements to polytene chromosomes, we documented phenotypic effects on female fertility and chromosomal mispairings that might be linked to DMI in hybrids. We propose that the inferred negative selection on heterozygous Cla-element insertions causes endogenous genetic barriers and therefore acts as DMI among C. riparius populations. The intrinsic genomic turnover exerted by TEs, thus, may have a direct impact on population divergence that is operationally different from drift and local adaptation.
Bears are iconic mammals with a complex evolutionary history. Natural bear hybrids and studies of few nuclear genes indicate that gene flow among bears may be more common than expected and not limited to the closely related polar and brown bears. Here we present a genome analysis of the bear family with representatives of all living species. Phylogenomic analyses of 869 mega base pairs divided into 18,621 genome fragments yielded a well-resolved coalescent species tree despite signals for extensive gene flow across species. However, genome analyses using three different statistical methods show that gene flow is not limited to closely related species pairs. Strong ancestral gene flow between the Asiatic black bear and the ancestor to polar, brown and American black bear explains numerous uncertainties in reconstructing the bear phylogeny. Gene flow across the bear clade may be mediated by intermediate species such as the geographically wide-spread brown bears leading to massive amounts of phylogenetic conflict. Genome-scale analyses lead to a more complete understanding of complex evolutionary processes. The increasing evidence for extensive inter-specific gene flow, found also in other animal species, necessitates shifting the attention from speciation processes achieving genome-wide reproductive isolation to the selective processes that maintain species divergence in the face of gene flow.
All giraffe (Giraffa) were previously assigned to a single species (G. Camelopardalis) and nine subspecies. However, multi-locus analyses of all subspecies have shown that there are four genetically distinct clades and suggest four giraffe species. This conclusion might not be fully accepted due to limited data and lack of explicit gene flow analyses. Here we present an extended study based on 21 independent nuclear loci from 137 individuals. Explicit gene flow analyses identify less than one migrant per generation, including between the closely related northern and reticulated giraffe. Thus, gene flow analyses and population genetics of the extended dataset confirm four genetically distinct giraffe clades and support four independent giraffe species. The new findings call for a revision of the IUCN classification of giraffe taxonomy. Three of the four species are threatened with extinction, mostly occurring in politically unstable regions, and as such, require the highest conservation support possible.
Welfare is the largest expenditure category in all advanced democracies. Consequently, much literature has studied partisan effects on total and policy-specific welfare expenditure. Yet, these results cannot be trusted: The methodological standard is to apply time-series cross-section-regressions to annual observation data. But governments hardly change annually. Thus, the number of observations is artificially inflated, leading to incorrect estimates. While this problem has recently been acknowledged, it has not been convincingly resolved. We propose Mixed-Effects Models as a solution, which allow decomposing variance into different levels and permit complex cross-classification data structures. We argue that Mixed-Effects models combine the strengths of existing methodological approaches while alleviating their weaknesses. Empirically, we study partisan effects on total and on disaggregated expenditure in 23 OECD-countries, 1960-2012, using several measures of party preferences.
The debate on effects of globalization on welfare states is extensive. Often couched in terms of a battle between the compensation and the efficiency thesis, the scholarly literature has provided contradictory arguments and findings. This article contributes to the scholarly debate by exploring in greater detail the micro-level foundations of compensation theory. More specifically, we distinguish between individual policy preferences for compensatory social policies (unemployment insurance) and human capital-focused social investment policies (education) and expect globalization to mainly affect demand for educational investment. A multi-level analysis of ISSP survey data provides empirical support for this hypothesis. This finding provides an important revision and extension of the classical analytical perspective of compensation theory, because it shows that citizens value the social investment function of the welfare state above and beyond simple compensation via social transfers. This might be particularly relevant in today's skill-centered knowledge economies.
We cannot imagine a political system without opposition. Despite this crucial position in politics, political science has largely neglected to study oppositions. Attempting to fill this gap, this article analyses the institutional opportunities of parliamentary oppositions. It offers a parsimonious framework by distinguishing two dimensions of opposition influence: Some institutions enable oppositions to control governments, while others offer opportunities to present alternatives. A comparison of oppositions’ opportunities in 21 democracies shows that countries fall into four groups along these dimensions: In majoritarian democracies, weak control mechanisms are countered by excellent opportunities to publicize alternatives. Consociational democracies are characterized by strong control mechanisms, but provide only weak opportunities to present alternatives. Moreover, in Southern Europe, control mechanisms and opportunities to present alternatives are weak, while both are pronounced in Nordic Europe. The results are summarized in three indices that can easily be applied in future research examining oppositions and their power.
Phylogenetic analyses of nuclear and mitochondrial genomes have shown that polar bears captured the mitochondrial genome of brown bears some 160,00 years ago. This hybridization event likely led to an extinction of the original polar bear mitochondrial genome. However, parts of the mitochondrial DNA occasionally integrates into the nuclear genome, forming pseudogenes called numts (nuclear mitochondrial integrations). Screening the polar bear genome for numts, we identified only 13 such integrations. Analyses of whole-genome sequences from additional polar bears, brown and American black bears as well as the giant panda indicates that the discovered numts entered the bear lineage before the initial ursid radiation some 14 million years ago. Our findings suggests a low integration rate of numts in the bear lineage and a complete loss of the original polar bear mitochondrial genome.
Background: Ever decreasing costs along with advances in sequencing and library preparation technologies enable even small research groups to generate chromosome-level assemblies today. Here we report the generation of an improved chromosome-level assembly for the Siamese fighting fish (Betta splendens) that was carried out during a practical university Master’s course. The Siamese fighting fish is a popular aquarium fish and an emerging model species for research on aggressive behaviour. We updated the current genome assembly by generating a new long-read nanopore-based assembly with subsequent scaffolding to chromosome-level using previously published HiC data.
Findings: The use of nanopore-based long-read data sequenced on a MinION platform (Oxford Nanopore Technologies) allowed us to generate a baseline assembly of only 1,276 contigs with a contig N50 of 2.1 Mbp, and a total length of 441 Mbp. Scaffolding using previously published HiC data resulted in 109 scaffolds with a scaffold N50 of 20.7 Mbp. More than 99% of the assembly is comprised in 21 scaffolds. The assembly showed the presence of 95.8% complete BUSCO genes from the Actinopterygii dataset indicating a high quality of the assembly.
Conclusion: We present an improved full chromosome-level assembly of the Siamese fighting fish generated during a university Master’s course. The use of ~35× long-read nanopore data drastically improved the baseline assembly in terms of continuity. We show that relatively in-expensive high-throughput sequencing technologies such as the long-read MinION sequencing platform can be used in educational settings allowing the students to gain practical skills in modern genomics and generate high quality results that benefit downstream research projects.
Borders and edges are salient and behaviourally relevant features for navigating the environment. The brain forms dedicated neural representations of environmental boundaries, which are assumed to serve as a reference for spatial coding. Here we expand this border coding network to include the retrosplenial cortex (RSC) in which we identified neurons that increase their firing near all boundaries of an arena. RSC border cells specifically encode walls, but not objects, and maintain their tuning in the absence of direct sensory detection. Unlike border cells in the medial entorhinal cortex (MEC), RSC border cells are sensitive to the animal’s direction to nearby walls located contralateral to the recorded hemisphere. Pharmacogenetic inactivation of MEC led to a disruption of RSC border coding, but not vice versa, indicating network directionality. Together these data shed light on how information about distance and direction of boundaries is generated in the brain for guiding navigation behaviour.
Improved integration of single cell transcriptome data demonstrated on heart failure in mice and men
(2023)
Biomedical research frequently uses murine models to study disease mechanisms. However, the translation of these findings to human disease remains a significant challenge. In order to improve the comparability of mouse and human data, we present a cross-species integration pipeline for single-cell transcriptomic assays.
The pipeline merges expression matrices and assigns clear orthologous relationships. Starting from Ensembl ortholog assignments, we allocated 82% of mouse genes to unique orthologs by using additional publicly available resources such as Uniprot, and NCBI databases. For genes with multiple matches, we employed the Needleman-Wunsch global alignment based on either amino acid or nucleotide sequence to identify the ortholog with the highest degree of similarity.
The workflow was tested for its functionality and efficiency by integrating scRNA-seq datasets from heart failure patients with the corresponding mouse model. We were able to assign unique human orthologs to up to 80% of the mouse genes, utilizing the known 17,492 orthologous pairs. Curiously, the integration process enabled the identification of both common and unique regulatory pathways between species in heart failure.
In conclusion, our pipeline streamlines the integration process, enhances gene nomenclature alignment and simplifies the translation of mouse models to human disease. We have made the OrthoIntegrate R-package accessible on GitHub (https://github.com/MarianoRuzJurado/OrthoIntegrate), which includes the assignment of ortholog definitions for human and mouse, as well as the pipeline for integrating single cells.
Recently, significant advances have been made by identifying the levels of synchronicity of the underlying dynamics of a given brain state. This research has demonstrated that unconscious dynamics tend to be more synchronous than those found in conscious states, which are more asynchronous. Here we go beyond this dichotomy to demonstrate that the different brain states are always underpinned by spatiotemporal chaos but with dissociable turbulent dynamics. We investigated human neuroimaging data from different brain states (resting state, meditation, deep sleep, and disorders of consciousness after coma) and were able to distinguish between them using complementary model-free and model-based measures of turbulent information transmission. Our model-free approach used recent advances describing a measure of information cascade across spatial scales using tools from turbulence theory. Complementarily, our model-based approach used exhaustive in silico perturbations of whole-brain models fitted to the empirical neuroimaging data, which allowed us to study the information encoding capabilities of the brain states. Overall, the current framework demonstrates that different levels of turbulent dynamics are fundamental for describing and differentiating between brain states.
BOLD signatures of sleep
(2019)
Sleep can be distinguished from wake by changes in brain electrical activity, typically assessed using electroencephalography (EEG). The hallmark of non-rapid-eye-movement sleep are two major EEG events: slow waves and spindles. Here we sought to identify possible signatures of sleep in brain hemodynamic activity, using simultaneous fMRI-EEG. We found that, during the transition from wake to sleep, blood-oxygen-level-dependent (BOLD) activity evolved from a mixed-frequency pattern to one dominated by two distinct oscillations: a low-frequency (~0.05Hz) oscillation prominent in light sleep and a high-frequency (~0.17Hz) oscillation in deep sleep. The two BOLD oscillations correlated with the occurrences of spindles and slow waves, respectively. They were detectable across the whole brain, cortically and subcortically, but had different regional distributions and opposite onset patterns. These spontaneous BOLD oscillations provide fMRI signatures of basic sleep processes, which may be employed to study human sleep at spatial resolution and brain coverage not achievable using EEG.
Summary: Understanding the role of short-interfering RNA (siRNA) in diverse biological processes is of current interest and often approached through small RNA sequencing. However, analysis of these datasets is difficult due to the complexity of biological RNA processing pathways, which differ between species. Several properties like strand specificity, length distribution, and distribution of soft-clipped bases are few parameters known to guide researchers in understanding the role of siRNAs. We present RAPID, a generic eukaryotic siRNA analysis pipeline, which captures information inherent in the datasets and automatically produces numerous visualizations as user-friendly HTML reports, covering multiple categories required for siRNA analysis. RAPID also facilitates an automated comparison of multiple datasets, with one of the normalization techniques dedicated for siRNA knockdown analysis, and integrates differential expression analysis using DESeq2. RAPID is available under MIT license at https://github.com/SchulzLab/RAPID. We recommend using it as a conda environment available from https://anaconda.org/bioconda/rapid.
Several studies suggested that transcription factor (TF) binding to DNA may be impaired or enhanced by DNA methylation. We present MeDeMo, a toolbox for TF motif analysis that combines information about DNA methylation with models capturing intra-motif dependencies. In a large-scale study using ChIP-seq data for 335 TFs, we identify novel TFs that are affected by DNA methylation. Overall, we find that CpG methylation decreases the likelihood of binding for the majority of TFs. For a considerable subset of TFs, we show that intra-motif dependencies are pivotal for accurately modelling the impact of DNA methylation on TF binding.
Background Enhancers play a fundamental role in orchestrating cell state and development. Although several methods have been developed to identify enhancers, linking them to their target genes is still an open problem. Several theories have been proposed on the functional mechanisms of enhancers, which triggered the development of various methods to infer promoter enhancer interactions (PEIs). The advancement of high-throughput techniques describing the three-dimensional organisation of the chromatin, paved the way to pinpoint long-range PEIs. Here we investigated whether including PEIs in computational models for the prediction of gene expression improves performance and interpretability.
Results We have extended our Tepic framework to include DNA contacts deduced from chromatin conformation capture experiments and compared various methods to determine PEIs using predictive modelling of gene expression from chromatin accessibility data and predicted transcription factor (TF) motif data. We found that including long-range PEIs deduced from both HiC and HiChIP data indeed improves model performance. We designed a novel machine learning approach that allows to prioritize TFs in distal loop and promoter regions with respect to their importance for gene expression regulation. Our analysis revealed a set of core TFs that are part of enhancer-promoter loops involving YY1 in different cell lines.
Conclusion: We show that the integration of chromatin conformation data improves gene expression prediction, underlining the importance of enhancer looping for gene expression regulation. Our general approach can be used to prioritize TFs that are involved in distal and promoter-proximal regulation using accessibility, conformation and expression data.
Interest in time-resolved connectivity in fMRI has grown rapidly in recent years. The most widely used technique for studying connectivity changes over time utilizes a sliding windows approach. There has been some debate about the utility of shorter versus longer windows, the use of fixed versus adaptive windows, as well as whether observed resting state dynamics during wakefulness may be predominantly due to changes in sleep state and subject head motion. In this work we use an independent component analysis (ICA)-based pipeline applied to concurrent EEG/fMRI data collected during wakefulness and various sleep stages and show: 1) connectivity states obtained from clustering sliding windowed correlations of resting state functional network time courses well classify the sleep states obtained from EEG data, 2) using shorter sliding windows instead of longer non-overlapping windows improves the ability to capture transition dynamics even at windows as short as 30 seconds, 3) motion appears to be mostly associated with one of the states rather than spread across all of them 4) a fixed tapered sliding window approach outperforms an adaptive dynamic conditional correlation approach, and 5) consistent with prior EEG/fMRI work, we identify evidence of multiple states within the wakeful condition which are able to be classified with high accuracy. Classification of wakeful only states suggest the presence of time-varying changes in connectivity in fMRI data beyond sleep state or motion. Results also inform about advantageous technical choices, and the identification of different clusters within wakefulness that are separable suggest further studies in this direction.
Mitochondrial NADH:ubiquinone oxidoreductase (complex I) is a 1 MDa membrane protein complex with a central role in energy metabolism. Redox-driven proton translocation by complex I contributes substantially to the proton motive force that drives ATP synthase. Several structures of complex I from bacteria and mitochondria have been determined but its catalytic mechanism has remained controversial. We here present the cryo-EM structure of complex I from Yarrowia lipolytica at 2.1 Å resolution, which reveals the positions of more than 1600 protein-bound water molecules, of which ∼100 are located in putative proton translocation pathways. Another structure of the same complex under steady-state activity conditions at 3.4 Å resolution indicates conformational transitions that we associate with proton injection into the central hydrophilic axis. By combining high-resolution structural data with site-directed mutagenesis and large-scale molecular dynamics simulations, we define details of the proton translocation pathways, and offer new insights into the redox-coupled proton pumping mechanism of complex I.
Brookshire (2022) claims that previous analyses of periodicity in detection performance after a reset event suffer from extreme false-positive rates. Here we show that this conclusion is based on an incorrect implemention of a null-hypothesis of aperiodicity, and that a correct implementation confirms low false-positive rates. Furthermore, we clarify that the previously used method of shuffling-in-time, and thereby shuffling-in-phase, cleanly implements the null hypothesis of no temporal structure after the reset, and thereby of no phase locking to the reset. Moving from a corresponding phase-locking spectrum to an inference on the periodicity of the underlying process can be accomplished by parameterizing the spectrum. This can separate periodic from non-periodic components, and quantify the strength of periodicity.
Human behaviour is inextricably linked to the interaction of emotion and cognition. For decades, emotion and cognition were perceived as separable processes, yet with mutual interactions. Recently, this differen-tiation has been challenged by more integrative approaches, but without addressing the exact neurophysiological basis of their interaction. Here, we aimed to uncover neurophysiological mechanisms of emotion-cognition interaction. We used an emotional Flanker task paired with EEG/FEM beamforming in a large cohort (N=121) of healthy human participants, obtaining high temporal and fMRI-equivalent spatial resolution. Spatially, emotion and cognition processing overlapped in the right inferior frontal gyrus (rIFG), specifically in pars triangularis. Temporally, emotion and cognition processing overlapped during the transition from emotional to cognitive processing, with a stronger interaction in β-band power leading to worse behavioral performance. Despite functionally segregated subdivisions in rIFG, frequency-specific information flowed extensively within IFG and top-down to visual areas (V2, Precuneus) – explaining the behavioral interference effect. Thus, for the first time we here show the neural mechanisms of emotion-cognition interaction in space, time, frequency and information transfer with high temporal and spatial resolution, revealing a central role for β-band activity in rIFG. Our results support the idea that rIFG plays a broad role in both inhibitory control and emotional interference inhibition as it is a site of convergence in both processes. Furthermore, our results have potential clinical implications for understanding dysfunctional emotion-cognition interaction and emotional interference inhibition in psychiatric disor-ders, e.g. major depression and substance use disorder, in which patients have difficulties in regulating emotions and executing inhibitory control.
Summary statement When echolocating under demanding conditions e.g. noisy, narrow space, or cluttered environments, frugivorous bats adapt their call pattern by increasing the call rate within biosonar groups.
Abstract For orientation, echolocating bats emit biosonar calls and use echoes arising from call reflections. They often pattern their calls into groups which increases the rate of sensory feedback over time. Insectivorous bats emit call groups at a higher rate when orienting in cluttered compared to uncluttered environments. Frugivorous bats increase the rate of call group emission when they echolocate in noisy environments. Here, calls emitted by conspecifics potentially interfere with the bat’s biosonar signals and complicate the echolocation behavior. To minimize the information loss followed by signal interference, bats may profit from a temporally increased sensory acquisition rate, as it is the case for the call groups. In frugivorous bats, it remains unclear if call group emission represents an exclusive adaptation to avoid interference by signals from other bats or if it represents an adaptation that allows to orient under demanding environmental conditions. Here, we compared the emission pattern of the frugivorous bat Carollia perspicillata when the bats were flying in noisy versus silent, narrow versus wide or cluttered versus non-cluttered corridors. According to our results, the bats emitted larger call groups and they increased the call rate within the call groups when navigating in narrow, cluttered, or noisy environments. Thus, call group emission represents an adaptive behavior when the bats orient in complex environments.
Most mammals rely on the extraction of acoustic information from the environment in order to survive. However, the mechanisms that support sound representation in auditory neural networks involving sensory and association brain areas remain underexplored. In this study, we address the functional connectivity between an auditory region in frontal cortex (the frontal auditory field, FAF) and the auditory cortex (AC) in the bat Carollia perspicillata. The AC is a classic sensory area central for the processing of acoustic information. On the other hand, the FAF belongs to the frontal lobe, a brain region involved in the integration of sensory inputs, modulation of cognitive states, and in the coordination of behavioural outputs. The FAF-AC network was examined in terms of oscillatory coherence (local-field potentials, LFPs), and within an information theoretical framework linking FAF and AC spiking activity. We show that in the absence of acoustic stimulation, simultaneously recorded LFPs from FAF and AC are coherent in low frequencies (1-12 Hz). This “default” coupling was strongest in deep AC layers and was unaltered by acoustic stimulation. However, presenting auditory stimuli did trigger the emergence of coherent auditory-evoked gamma-band activity (>25 Hz) between the FAF and AC. In terms of spiking, our results suggest that FAF and AC engage in distinct coding strategies for representing artificial and natural sounds. Taken together, our findings shed light onto the neuronal coding strategies and functional coupling mechanisms that enable sound representation at the network level in the mammalian brain.
The mammalian frontal and auditory cortices are important for vocal behaviour. Here, using local field potential recordings, we demonstrate for the first time that the timing and spatial pattern of oscillations in the fronto-auditory cortical network of vocalizing bats (Carollia perspicillata) predict the purpose of vocalization: echolocation or communication. Transfer entropy analyses revealed predominantly top-down (frontal-to-auditory cortex) information flow during spontaneous activity and pre-vocal periods. The dynamics of information flow depended on the behavioural role of the vocalization and on the timing relative to vocal onset. Remarkably, we observed the emergence of predominantly bottom-up (auditory-to-frontal cortex) information transfer patterns specific echolocation production, leading to self-directed acoustic feedback. Electrical stimulation of frontal areas selectively enhanced responses to echolocation sounds in auditory cortex. These results reveal unique changes in information flow across sensory and frontal cortices, potentially driven by the purpose of the vocalization in a highly vocal mammalian model.
The brains of black 6 mice (Mus musculus) and Seba’s short-tailed bats (Carollia perspicillata) weigh roughly the same and share the mammalian neocortical laminar architecture. Bats have highly developed sonar calls and social communication and are an excellent neuroethological animal model for auditory research. Mice are olfactory and somatosensory specialists and are used frequently in auditory neuroscience, particularly for their advantage of standardization and genetic tools. Investigating their potentially different general auditory processing principles would advance our understanding of how the ecological needs of a species shape the development and function of the mammalian nervous system. We compared two existing datasets, recorded with linear multichannel electrodes down the depth of the primary auditory cortex (A1) while awake, across both species while presenting repetitive stimulus trains with different frequencies (∼5 and ∼40 Hz). We found that while there are similarities between cortical response profiles in bats and mice, there was a better signal to noise ratio in bats under these conditions, which allowed for a clearer following response to stimuli trains. This was most evident at higher frequency trains, where bats had stronger response amplitude suppression to consecutive stimuli. Phase coherence was far stronger in bats during stimulus response, indicating less phase variability in bats across individual trials. These results show that although both species share cortical laminar organization, there are structural differences in relative depth of layers. Better signal to noise ratio in bats could represent specialization for faster temporal processing shaped by their individual ecological niches.
The ability to vocalize is ubiquitous in vertebrates, but neural networks leading to vocalization production remain poorly understood. Here we performed simultaneous, large scale, neuronal recordings in the frontal cortex and dorsal striatum (caudate nucleus) during the production of echolocation and non-echolocation calls in bats. This approach allows to assess the general aspects underlying vocalization production in mammals and the unique evolutionary adaptations of bat echolocation. Our findings show that distinct intra-areal brain rhythms in the beta (12-30 Hz) and gamma (30-80 Hz) bands of the local field potential can be used to predict the bats’ vocal output and that phase locking between spikes and field potentials occurs prior vocalization production. Moreover, the fronto-striatal network is differentially coupled in the theta-band during the production of echolocation and non-echolocation calls. Overall, our results present evidence for fronto-striatal network oscillations in motor action prediction in mammals.
Consistent individual differences in behaviour (animal personality) are widespread throughout the Animal Kingdom. This includes variation in risk-taking versus risk-averse behavioural tendencies. Variation in several personality dimensions is associated with distinct fitness consequences and thus, may become a target of natural and/or sexual selection. However, the link between animal personality and mate choice—as a major component of sexual selection—remains understudied. We asked (1) whether females and males of the livebearing fish Poecilia mexicana prefer risk-taking mating partners (directional mating preference), (2) or if their preferences are dependent on the choosing individual’s own personality type (assortative mating). We characterized each test subject for its risk-taking behaviour, assessed as the time to emerge from shelter and enter an unknown area. In dichotomous association preference tests, we offered two potential mating partners that differed in risk-taking behaviour but were matched for other phenotypic traits (body size, shape, and colouration). Females, but not males, exhibited a strong directional preference for risk-taking over risk-averse mating partners. At the same time, the strength of females’ preferences correlated positively with their own risk-taking scores. Our study is the first to demonstrate that a strong overall preference for risk-taking mating partners does not preclude effects of choosing individuals’ own personality type on (subtle) individual variation in mating preferences. More generally, two different preferences functions appear to interact to determine the outcome of individual mate choice decisions.
Cognition requires the dynamic modulation of effective connectivity, i.e. the modulation of the postsynaptic neuronal response to a given input. If postsynaptic neurons are rhythmically active, this might entail rhythmic gain modulation, such that inputs synchronized to phases of high gain benefit from enhanced effective connectivity. We show that visually induced gamma-band activity in awake macaque area V4 rhythmically modulates responses to unpredictable stimulus events. This modulation exceeded a simple additive superposition of a constant response onto ongoing gamma-rhythmic firing, demonstrating the modulation of multiplicative gain. Gamma phases leading to strongest neuronal responses also led to shortest behavioral reaction times, suggesting functional relevance of the effect. Furthermore, we find that constant optogenetic stimulation of anesthetized cat area 21a produces gamma-band activity entailing a similar gain modulation. As the gamma rhythm in area 21a did not spread backwards to area 17, this suggests that postsynaptic gamma is sufficient for gain modulation.
Synchronization has been implicated in neuronal communication, but causal evidence remains indirect. We used optogenetics to generate depolarizing currents in pyramidal neurons of cat visual cortex, emulating excitatory synaptic inputs under precise temporal control, while measuring spike output. Cortex transformed constant excitation into strong gamma-band synchronization, revealing the well-known cortical resonance. Increasing excitation with ramps increased the strength and frequency of synchronization. Slow, symmetric excitation profiles revealed hysteresis of power and frequency. Crucially, white-noise input sequences enabled causal analysis of network transmission, establishing that cortical resonance selectively transmits coherent input components. Models composed of recurrently coupled excitatory and inhibitory units uncovered a crucial role of feedback inhibition and suggest that hysteresis can arise through spike-frequency adaptation. The presented approach provides a powerful means to investigate the resonance properties of local circuits and probe how these properties transform input and shape transmission.
The gamma rhythm has been implicated in neuronal communication, but causal evidence remains indirect. We measured spike output of local neuronal networks and emulated their synaptic input through optogenetics. Opsins provide currents through somato-dendritic membranes, similar to synapses, yet under experimental control with high temporal precision. We expressed Channelrhodopsin-2 in excitatory neurons of cat visual cortex and recorded neuronal responses to light with different temporal characteristics. Sine waves of different frequencies entrained neuronal responses with a reliability that peaked for input frequencies in the gamma band. Crucially, we also presented white-noise sequences, because their temporal unpredictability enables analysis of causality. Neuronal spike output was caused specifically by the input’s gamma component. This gamma-specific transfer function is likely an emergent property of in-vivo networks with feedback inhibition. The method described here could reveal the transfer function between the input to any one and the output of any other neuronal group.
The plant family Brassicaceae includes some of the most studied hosts of plant microbiomes, targeting microbial diversity, community assembly rules, and effects on host performance. Compared to bacteria, eukaryotes in the brassicaceous microbiome remain understudied, especially under natural settings. Here, we assessed the impact of host identity and age on the assembly of fungal and oomycete root communities, using DNA metabarcoding of roots and associated soil of three annual co-habiting Brassicaceae collected at two time points. Our results showed that fungal communities are more diverse and structured than those of oomycetes. In both cases, plant identity and sampling time had little influence on community variation, whereas root/soil compartment had a strong effect by exerting control on the entry of soil microorganisms into the roots. The enrichment in roots of specific fungi suggests a specialization towards the asymptomatic colonization of plant tissues, which could be relevant to host’s fitness and health.
Wildfires are relatively rare in subarctic tundra ecosystems, but they can strongly change ecosystem properties. Short-term fire effects on subarctic tundra vegetation are well documented, but long-term vegetation recovery has been studied less. The frequency of tundra fires will increase with climate warming. Understanding the long-term effects of fire is necessary to predict future ecosystem changes.
We used a space-for-time approach to assess vegetation recovery after fire over more than four decades. We studied soil and vegetation patterns on three large fire scars (>44, 28 and 12 years old) in dry, lichen-dominated forest tundra in Western Siberia. On 60 plots, we determined soil temperature and permafrost thaw depth, sampled vegetation and measured plant functional traits. We assessed trends in NDVI to support the field-based results on vegetation recovery.
Soil temperature, permafrost thaw depth and total vegetation cover had recovered to pre-fire levels after >44 years, as well as total vegetation cover. In contrast, after >44 years, functional groups had not recovered to the pre-fire state. Burnt areas had lower lichen and higher bryophyte and shrub cover. The dominating shrub species, Betula nana, exhibited a higher vitality (higher specific leaf area and plant height) on burnt compared with control plots, suggesting a fire legacy effect in shrub growth. Our results confirm patterns of shrub encroachment after fire that were detected before in other parts of the Arctic and Subarctic. In the so far poorly studied Western Siberian forest tundra we demonstrate for the first time, long-term fire-legacies on the functional composition of relatively dry shrub- and lichen-dominated vegetation.
Seed harvesting from wild plant populations is key for ecological restoration, but may threaten the persistence of source populations. Consequently, several countries have set guidelines limiting the proportions of harvestable seeds. Here, we use high-resolution data from 298 plant species to model the demographic consequences of seed harvesting. We find that the current guidelines only protect some species, but are insufficient or overly restrictive for others. We show that the maximum possible fraction of seed harvesting is strongly associated with harvesting frequency and generation time of the target species, ranging from 100% in long-lived species to <1% in the most annuals. Our results provide quantitative basis to guide seed harvesting legislation based on species’ generation time and harvesting regime.
Animals sense ambient temperature so that they can adjust their behavior to the environment; they avoid noxious heat and coldness and stay within a survivable temperature range. C. elegans can sense temperature, memorize past cultivation temperature and navigate towards preferable temperature, for which a thermosensory neuron, AFD, is essential. AFD responds to temperature increase from the past cultivation temperature by increasing intracellular Ca2+ level. We aimed to reveal how AFD encodes and memorizes the information of temperature. Although cGMP synthesis is crucial for thermosensation by AFD, whether and how cGMP level temporally fluctuates in AFD remained elusive. We therefore monitored cGMP level in AFD and found that cGMP dynamically responded to temperature change in a manner dependent on past cultivation temperature. Given that cGMP dynamics is supposed to be upstream of Ca2+ dynamics, our results suggest that AFD’s memory is formed by simpler molecular mechanisms than previously expected from the Ca2+ dynamics. Moreover, we analyzed how guanylyl cyclases and phosphodiesterases, which synthesize and degrade cGMP, respectively, contributed to cGMP and Ca2+ dynamics and thermotaxis behavior.
Signal transfer of visual stimuli to V4 occurs in gamma-rhythmic, pulsed information packages
(2020)
Summary Selective visual attention allows the brain to focus on behaviorally relevant information while ignoring irrelevant signals. As a possible mechanism, routing by synchronization was proposed: neural populations sending attended signals align their gamma-rhythmic activities with receiving populations, such that spikes from the senders arrive at excitability peaks of the receivers, enhancing signal transfer. Conversely, the non-attended signals arrive unaligned to the receiver’s oscillation, reducing signal transfer. Therefore, visual signals should be transferred through periodically pulsed information packages, resulting in a modulation of the stimulus content within the receiver’s activity by its gamma phase and amplitude. To test this prediction, we quantified gamma phase-specific stimulus content within neural activity from area V4 of macaques performing a visual attention task. For the attended stimulus we find enhanced stimulus content reaching its maximum near excitability peaks, with effect magnitude increasing with oscillation amplitude, establishing a functional link between selective processing and gamma activity.
No disease modifying therapy is currently available for Parkinson’s disease (PD), the second most common neurodegenerative disease. The long non-motor prodromal phase of PD is a window of opportunity for early detection and intervention. However, we lack the pathophysiological understanding to develop selective biomarkers and interventions. By developing a mutant α-synuclein selective-overexpression mouse model of prodromal PD, we identified a cell-autonomous selective Kv4 channelopathy in dorsal motor nucleus of the vagus (DMV) neurons. This functional remodeling of intact DMV neurons leads to impaired pacemaker function in vitro and in vivo, which in turn reduces gastrointestinal motility which is a common, very early symptom of prodromal PD. We show for the first time a causal chain of events from α-synuclein via a biophysical dysfunction of specific neuronal populations to a clinically relevant prodromal symptom. These findings can facilitate the rational design of clinical biomarkers to identify people at risk for PD.
Substantia nigra dopamine (SN DA) neurons are progressively lost in Parkinson disease (PD). While the molecular and cellular mechanisms of their differential vulnerability and degeneration have been extensively studied, we still know very little about potential functional adaptations of those SN DA neurons that – at least for some time – manage to survive during earlier stages of PD. We utilized a partial lesion 6-OHDA mouse model to characterize initial electrophysiological impairments and chronic adaptations of surviving identified SN DA neurons, both in vivo and in vitro. Early after lesion (3 weeks), we detected a selective loss of in vivo burst firing in surviving SN DA neurons, which was accompanied by in vitro pacemaker instability. In contrast, late after lesion (>2 months), in vivo firing properties of surviving SN DA neurons had recovered in the presence of 2-fold accelerated pacemaking in vitro. Finally, we show that this chronic cell-autonomous adaptation in surviving SN DA neurons was mediated by Kv4.3 channel downregulation. Our study demonstrates substantial homeostatic plasticity of surviving SN DA neurons after a single-hit non-progressive lesion, which might contribute to the phenotype of initially surviving SN DA neurons in PD.
Parkinson disease (PD), one of the most common neurodegenerative disorder, is believed to be driven by toxic α-synuclein aggregates eventually resulting in selective loss of vulnerable neuron populations, prominent among them, nigrostriatal dopamine (DA) neurons in the lateral substantia nigra (l-SN). How α-synuclein aggregates initiate a pathophysiological cascade selectively in vulnerable neurons is still unclear. Here, we show that the exposure to low nanomolar concentrations of α-synuclein aggregates (i.e. fibrils) but not its monomeric forms acutely and selectively disrupted the electrical pacemaker function of the DA subpopulation most vulnerable in PD. This implies that only dorsolateral striatum projecting l-SN DA neurons were electrically silenced by α-synuclein aggregates, while the activity of neither neighboring DA neurons in medial SN projecting to dorsomedial striatum nor mesolimbic DA neurons in the ventral tegmental area (VTA) were affected. Moreover, we demonstrate functional K-ATP channels comprised of Kir6.2 subunit in DA neurons to be necessary to mediate this acute pacemaker disruption by α-synuclein aggregates. Our study thus identifies a molecularly defined target that quickly translates the presence of α-synuclein aggregates into an immediate impairment of essential neuronal function. This constitutes a novel candidate process how a protein-aggregation-driven sequence in PD is initiated that might eventually lead to selective neurodegeneration.
Mechanisms by which specific histone modifications regulate distinct gene regulatory networks remain little understood. We investigated how H3K79me2, a modification catalyzed by DOT1L and previously considered a general transcriptional activation mark, regulates gene expression in mammalian cardiogenesis. Early embryonic cardiomyocyte ablation of Dot1l revealed that H3K79me2 does not act as a general transcriptional activator, but rather regulates highly specific gene regulatory networks at two critical cardiogenic junctures: left ventricle patterning and postnatal cardiomyocyte cell cycle withdrawal. Mechanistic analyses revealed that H3K79me2 in two distinct domains, gene bodies and regulatory elements, synergized to promote expression of genes activated by DOT1L. Surprisingly, these analyses also revealed that H3K79me2 in specific regulatory elements contributed to silencing genes usually not expressed in cardiomyocytes. As DOT1L mutants had increased numbers of postnatal mononuclear cardiomyocytes and prolonged cardiomyocyte cell cycle activity, controlled inhibition of DOT1L might be a strategy to promote cardiac regeneration post-injury.
Background Transposable elements (TEs) are an important source of genome plasticity across the tree of life. Accumulating evidence suggests that TEs may not be randomly distributed in the genome. Drift and natural selection are important forces shaping TE distribution and accumulation, acting directly on the TE element or indirectly on the host species. Fungi, with their multifaceted phenotypic diversity and relatively small genome size, are ideal models to study the role of TEs in genome evolution and their impact on the host’s ecological and life history traits. Here we present an account of all TEs found in a high-quality reference genome of the lichen-forming fungus Umbilicaria pustulata, a macrolichen species comprising two climatic ecotypes: Mediterranean and cold-temperate. We trace the occurrence of the newly identified TEs in populations along three replicated elevation gradients using a Pool-Seq approach, to identify TE insertions of potential adaptive significance.
Results We found that TEs cover 21.26 % of the 32.9 Mbp genome, with LTR Gypsy and Copia clades being the most common TEs. Out of a total of 182 TE copies we identified 28 insertions displaying consistent insertion frequency differences between the two host ecotypes across the elevation gradients. Most of the highly differentiated insertions were located near genes, indicating a putative function.
Conclusions This pioneering study into the content and climate niche-specific distribution of TEs in a lichen-forming fungus contributes to understanding the roles of TEs in fungal evolution. Particularly, it may serve as a foundation for assessing the impact of TE dynamics on fungal adaptation to the abiotic environment, and the impact of TE activity on the evolution and maintenance of a symbiotic lifestyle.