Exploration of tissue-specific gene expression patterns underlying timing of breeding in contrasting temperature environments in a song bird

Avian migration clocks in a changing world

Avian long-distance migration requires refined programming to orchestrate the birds' movements on annual temporal and continental spatial scales. Programming is particularly important as long-distance movements typically anticipate future environmental conditions. Hence, migration has long been of particular interest in chronobiology. Captivity studies using a proxy, the shift to nocturnality during migration seasons (i.e., migratory restlessness), have revealed circannual and circadian regulation, as well as an innate sense of direction. Thanks to rapid development of tracking technology, detailed information from free-flying birds, including annual-cycle data and actograms, now allows relating this mechanistic background to behaviour in the wild. Likewise, genomic approaches begin to unravel the many physiological pathways that contribute to migration. Despite these advances, it is still unclear how migration programmes are integrated with specific environmental conditions experienced during the journey. Such knowledge is imminently important as temporal environments undergo rapid anthropogenic modification. Migratory birds as a group are not dealing well with the changes, yet some species show remarkable adjustments at behavioural and genetic levels. Integrated research programmes and interdisciplinary collaborations are needed to understand the range of responses of migratory birds to environmental change, and more broadly, the functioning of timing programmes under natural conditions.

A cross-taxonomic perspective on the integration of temperature cues in vertebrate seasonal neuroendocrine pathways

The regulation of seasonality has been an area of interest for decades, yet global climate change has created extra urgency in the quest to understand how sensory circuits and neuroendocrine control systems interact to generate flexibility in biological timekeeping. The capacity of temperature to alter endogenous or photoperiod-regulated neuroendocrine mechanisms driving seasonality, either as a direct cue or through temperature-dependent effects on energy and metabolism, is at the heart of this phenological flexibility. However, until relatively recently, little research had been done on the integration of temperature information in canonical seasonal neuroendocrine pathways, particularly in vertebrates. We review recent advances from research in vertebrates that deepens our understanding of how temperature cues are perceived and integrated into seasonal hypothalamic thyroid hormone (TH) signaling, which is a critical regulator of downstream seasonal phenotypic changes such as those regulated by the BPG (brain-pituitary-gonadal) axis. Temperature perception occurs through cutaneous transient receptor potential (TRP) neurons, though sensitivity of these neurons varies markedly across taxa. Although photoperiod is the dominant cue used to trigger seasonal physiology or entrain circannual clocks, across birds, mammals, fish, reptiles and amphibians, seasonality appears to be temperature sensitive and in at least some cases this appears to be related to phylogenetically conserved TH signaling in the hypothalamus. Nevertheless, the exact mechanisms through which temperature modulates seasonal neuroendocrine pathways remains poorly understood.

The human "contaminome": bacterial, viral, and computational contamination in whole genome sequences from 1000 families

The unmapped readspace of whole genome sequencing data tends to be large but is often ignored. We posit that it contains valuable signals of both human infection and contamination. Using unmapped and poorly aligned reads from whole genome sequences (WGS) of over 1000 families and nearly 5000 individuals, we present insights into common viral, bacterial, and computational contamination that plague whole genome sequencing studies. We present several notable results: (1) In addition to known contaminants such as Epstein-Barr virus and phiX, sequences from whole blood and lymphocyte cell lines contain many other contaminants, likely originating from storage, prep, and sequencing pipelines. (2) Sequencing plate and biological sample source of a sample strongly influence contamination profile. And, (3) Y-chromosome fragments not on the human reference genome commonly mismap to bacterial reference genomes. Both experiment-derived and computational contamination is prominent in next-generation sequencing data. Such contamination can compromise results from WGS as well as metagenomics studies, and standard protocols for identifying and removing contamination should be developed to ensure the fidelity of sequencing-based studies.

The time is ripe for functional genomics: Can epigenetic changes mediate reproductive timing?

Populations are under strong selection to match reproductive timing with favourable environmental conditions. This becomes particularly important and challenging with increasing interannual environmental variability. Adjusting reproductive timing requires the ability to sense and interpret relevant environmental cues, while responding flexibly to their interannual variation. For instance, in seasonal species, reproductive timing is often dependent on photoperiod and temperature. Although many genes influencing the timing of reproduction have been identified, far less attention has been paid to the gene-regulatory cascades orchestrating these complex gene-environment interactions. In a From the Cover article in this issue of Molecular Ecology, Lindner, Laine, et al. (2021) addressed this knowledge gap by investigating the role of DNA methylation in mediating reproductive timing in the seasonally breeding great tit (Parus major). Using a clever blood sampling design, they investigated genome-wide DNA methylation changes following individual female birds across multiple reproductive stages. This approach revealed 10 candidate genes with a strong correlation between promoter methylation and reproductive status. Some of these genes are known to be involved in reproductive timing (e.g., MYLK-like or NR5A1), yet for others this function was previously unknown (Figure 1). Interestingly, NR5A1 is a key transcription factor, which may affect other genes that are part of the same regulatory network. The findings of Lindner, Laine, et al. (2021) provide a strong case for studying DNA methylation to uncover how gene-environment interactions influence important life-history traits, such as reproductive timing.

Epigenetic effects of parasites and pesticides on captive and wild nestling birds

Anthropogenic changes to the environment challenge animal populations to adapt to new conditions and unique threats. While the study of adaptation has focused on genetic variation, epigenetic mechanisms may also be important. DNA methylation is sensitive to environmental stressors, such as parasites and pesticides, which may affect gene expression and phenotype. We studied the effects of an invasive ectoparasite, Philornis downsi, on DNA methylation of Galápagos mockingbirds (Mimus parvulus). We used the insecticide permethrin to manipulate P. downsi presence in nests of free-living mockingbirds and tested for effects of parasitism on nestling mockingbirds using epiGBS, a reduced-representation bisulfite sequencing (RRBS) approach. To distinguish the confounding effects of insecticide exposure, we conducted a matching experiment exposing captive nestling zebra finches (Taeniopygia guttata) to permethrin. We used zebra finches because they were the closest model organism to mockingbirds that we could breed in controlled conditions. We identified a limited number of differentially methylated cytosines (DMCs) in parasitized versus nonparasitized mockingbirds, but the number was not more than expected by chance. In contrast, we saw clear effects of permethrin on methylation in captive zebra finches. DMCs in zebra finches paralleled documented effects of permethrin exposure on vertebrate cellular signaling and endocrine function. Our results from captive birds indicate a role for epigenetic processes in mediating sublethal nontarget effects of pyrethroid exposure in vertebrates. Environmental conditions in the field were more variable than the laboratory, which may have made effects of both parasitism and permethrin harder to detect in mockingbirds. RRBS approaches such as epiGBS may be a cost-effective way to characterize genome-wide methylation profiles. However, our results indicate that ecological epigenetic studies in natural populations should consider the number of cytosines interrogated and the depth of sequencing in order to have adequate power to detect small and variable effects.

Effect of Post-Hatch Heat-Treatment in Heat-Stressed Transylvanian Naked Neck Chicken

Simple Summary

Heat stress due to high environmental temperature negatively influences animal productivity. Extensive studies have been carried out to evaluate the mechanisms of heat stress in chickens. It was shown that the expression level of heat-shock factors (HSFs) and heat-shock proteins (HSPs) were affected. Tissue-specific responses to the thermal challenge were also found in the heart, liver and muscle. Our study examined the changes in primary production parameters and four heat-shock factor and two heat-shock protein expression profiles in chicken gonads. In the first experiment, 24 h after hatching, 80 Transylvanian Naked Neck chickens were heat-treated at 38.5 °C ambient temperature with 60% humidity for 12 h. In this experiment, we studied the primary productivity parameters of matured chickens after the performed heat stress. In the second experiment, the heat treatment was the same, and we examined the expression pattern of heat-shock factors and heat-shock proteins in the control and treated gonads. We collected the samples immediately after the heat-treatment in case of half of the treated and control group. We found a significant difference in egg production, and increased expression level of HSP90 and HSF4 in heat-treated female gonads.

Abstract

Although numerous studies reported the effects of heat stress in chickens, it was not investigated in the Transylvanian Naked Neck breed. In our research, Transylvanian Naked Neck chickens, 24 h after hatching, were heat-treated at 38.5 °C for 12 h. We compared the control and heat-treated adult chickens’ productivity parameters following 12 weeks of heat-stress at 30 °C. We found that the heat-treated layers had significantly higher egg production in heat stress, but in cockerels, the sperm quality did not differ significantly between the two groups. To detect the effect of heat-treatment on a molecular level, the expression of two heat-shock proteins and four heat-shock factors were analysed in the gonads of control and heat-treated chickens. We found that the expression level of HSP90 and HSF4 increased significantly in heat-treated female chicken gonads. Still, in adult females, the expression of HSF2 and HSF3 were substantially lower compared to the control. In adult heat-treated males, the HSP70, HSF1 and HSF3 expression levels showed a significant increase in both gonads compared to the control. We think that the presented significant differences in egg production might be related to the increased expression level of HSP90 and HSF4 in heat-treated female gonads.

Rapid changes in DNA methylation associated with the initiation of reproduction in a small songbird

Species with a circannual life cycle need to match the timing of their life history events to the environment to maximize fitness. However, our understanding of how circannual traits such as timing of reproduction are regulated on a molecular level remains limited. Recent studies have implicated that epigenetic mechanisms can be an important part in the processes that regulate circannual traits. Here, we explore the role of DNA methylation in mediating reproductive timing in a seasonally breeding bird species, the great tit (Parus major), using genome‐wide DNA methylation data from individual females that were blood sampled repeatedly throughout the breeding season. We demonstrate rapid and directional changes in DNA methylation within the promoter region of several genes, including a key transcription factor (NR5A1) known from earlier studies to be involved in the initiation of timing of reproduction. Interestingly, the observed changes in DNA methylation at NR5A1 identified here are in line with earlier gene expression studies of reproduction in chicken, indicating that the observed shifts in DNA methylation at this gene can have a regulatory role. Our findings provide an important step towards elucidating the genomic mechanism that mediates seasonal timing of a key life history traits and provide support for the idea that epigenetic mechanisms may play an important role in circannual traits.

see also the Perspective by Melanie J. Heckwolf and Britta S. Meyer

Temporal changes in DNA methylation and RNA expression in a small song bird: within- and between-tissue comparisons

Background

DNA methylation is likely a key mechanism regulating changes in gene transcription in traits that show temporal fluctuations in response to environmental conditions. To understand the transcriptional role of DNA methylation we need simultaneous within-individual assessment of methylation changes and gene expression changes over time. Within-individual repeated sampling of tissues, which are essential for trait expression is, however, unfeasible (e.g. specific brain regions, liver and ovary for reproductive timing). Here, we explore to what extend between-individual changes in DNA methylation in a tissue accessible for repeated sampling (red blood cells (RBCs)) reflect such patterns in a tissue unavailable for repeated sampling (liver) and how these DNA methylation patterns are associated with gene expression in such inaccessible tissues (hypothalamus, ovary and liver). For this, 18 great tit (Parus major) females were sacrificed at three time points (n = 6 per time point) throughout the pre-laying and egg-laying period and their blood, hypothalamus, ovary and liver were sampled.

Results

We simultaneously assessed DNA methylation changes (via reduced representation bisulfite sequencing) and changes in gene expression (via RNA-seq and qPCR) over time. In general, we found a positive correlation between changes in CpG site methylation in RBCs and liver across timepoints. For CpG sites in close proximity to the transcription start site, an increase in RBC methylation over time was associated with a decrease in the expression of the associated gene in the ovary. In contrast, no such association with gene expression was found for CpG site methylation within the gene body or the 10 kb up- and downstream regions adjacent to the gene body.

Conclusion

Temporal changes in DNA methylation are largely tissue-general, indicating that changes in RBC methylation can reflect changes in DNA methylation in other, often less accessible, tissues such as the liver in our case. However, associations between temporal changes in DNA methylation with changes in gene expression are mostly tissue- and genomic location-dependent. The observation that temporal changes in DNA methylation within RBCs can relate to changes in gene expression in less accessible tissues is important for a better understanding of how environmental conditions shape traits that temporally change in expression in wild populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07329-9.

Identification of key genes and expression profiles in osteoarthritis by co-expressed network analysis

BACKGROUND

The underlying molecular characteristics of osteoarthritis (OA), a common age-related joint disease, remains elusive. Here, we aimed to identify potential early diagnostic biomarkers and elucidate underlying mechanisms of OA using weighted gene co-expression network analysis (WGCNA).

MATERIAL AND METHODS

We obtained the gene expression profile dataset GSE55235, GSE55457, and GSE55584, from the Gene Expression Omnibus. WGCNA was used to investigate the changes in co-expressed genes between normal and OA synovial membrane samples. Modules that were highly correlated to OA were subjected to functional enrichment analysis using the R clusterProfiler package. Differentially expressed genes (DEGs) between the two samples were screened using the "limma" package in R. A Venn diagram was constructed to intersect the genes in significant modules and DEGs. RT -PCR was used to further verify the hub gene expression levels between normal and OA samples.

RESULTS

The preserved significant module was found to be highly associated with OA development and progression (P < 1e-200, correlation = 0.92). Functional enrichment analysis suggested that the antiquewhite4 module was highly correlated to FoxO signaling pathway, and the metabolism of fatty acids and 2-oxocarboxylic acid. A total of 13 hub genes were identified based on significant module network topology and DEG analysis, and RT-PCR confirmed that these genes were significantly increased in OA samples compared with that in normal samples.

CONCLUSIONS

We identified 13 hub genes correlated to the development and progression of OA, which may provide new biomarkers and drug targets for OA.

Influence of heat stress on reference genes stability in heart and liver of two chickens genotypes

Introduction

Real-time polymerase chain reaction (RT-qPCR) is an important tool for analyzing gene expression. However, before analyzing the expression of target genes, it is crucial to normalize the reference genes, in order to find the most stable gene to be used as an endogenous control. A gene that remains stable in all samples under different treatments is considered a suitable normalizer. In this sense, we aimed to identify stable reference genes for normalization of target genes in the heart and liver tissues from two genetically divergent groups of chickens (Cobb 500^® commercial line and Peloco backyard chickens) under comfort and acute heat stress environmental conditions. Eight reference genes (ACTB, HPRT1, RPL5, EEF1, MRPS27, MRPS30, TFRC and LDHA) were analyzed for expression stability. The samples were obtained from 24 chickens, 12 from the backyard Peloco and 12 from the Cobb 500® line, exposed to two environmental conditions (comfort and heat stress). Comfort temperature was 23°C and heat stress temperature was 39.5°C for one hour. Subsequently, the animals were euthanized, and heart and liver tissue fragments were collected for RNA extraction and amplification. To determine the stability rate of gene expression, three different statistical algorithms were applied: BestKeeper, geNorm and NormFinder, and to obtain an aggregated stability list, the RankAgregg package of R software was used.

Results

The most stable genes using BestKeeper tool, including the two factors (genetic group and environmental condition), were LDHA, RPL5 and MRPS27 for heart tissue, and TFRC, RPL5 and EEF1 for liver tissue. Applying geNorm algorithm, the best reference genes were RPL5, EEF1 and MRPS30 for heart tissue and LDHA, EEF1 and RPL5 for liver. Using the NormFinder algorithm, the best normalizer genes were EEF1, RPL5 and LDHA in heart, and EEF1, RPL5 and ACTB in liver tissue. In the overall ranking obtained by RankAggreg package, considering the three algorithms, the RPL5, EEF1 and LDHA genes were the most stable for heart tissue, whereas RPL5, EEF1 and ACTB were the most stable for liver tissue.

Conclusion

According to the RankAggreg tool classification based on the three different algorithms (BestKeeper, geNorm and NormFinder), the most stable genes were RPL5, EEF1 and LDHA for heart tissue and RPL5, EEF1 and ACTB for liver tissue of chickens subjected to comfort and acute heat stress environmental conditions. However, the best reference genes may vary depending on the experimental conditions of each study, such as different breeds, environmental stressors, and tissues analyzed. Therefore, the need to perform priori studies to assay the best reference genes at the outset of each study is emphasized.