Search

Null alleles are alleles that for various reasons fail to amplify in a PCR assay. The presence of null alleles in microsatellite data is known to bias the genetic parameter estimates. Thus, efficient detection of null alleles is crucial, but the methods available for indirect null allele detection return inconsistent results. Here, our aim was to compare different methods for null allele detection, to explain their respective performance and to provide improvements. We applied several approaches to identify the 'true' null alleles based on the predictions made by five different methods, used either individually or in combination. First, we introduced simulated 'true' null alleles into 240 population data sets and applied the methods to measure their success in detecting the simulated null alleles. The single best-performing method was ML- Null Freq_frequency. Furthermore, we applied different noise reduction approaches to improve the results. For instance, by combining the results of several methods, we obtained more reliable results than using a single one. Rule-based classification was applied to identify population properties linked to the false discovery rate. Rules obtained from the classifier described which population genetic estimates and loci characteristics were linked to the success of each method. We have shown that by simulating 'true' null alleles into a population data set, we may define a null allele frequency threshold, related to a desired true or false discovery rate. Moreover, using such simulated data sets, the expected null allele homozygote frequency may be estimated independently of the equilibrium state of the population. [ABSTRACT FROM AUTHOR]

The satiety hormone leptin has a key role in the control of energy homeostasis in mammals. Leptin is produced almost exclusively by the adipose tissue and signals the amount of fat stores primarily to the hypothalamus, but also to peripheral tissues. Recently, we identified the long sought leptin (SEROUSSI et al., 2016), in the "dark side of the chicken genome", meaning extremely GC-rich regions that have been difficult to sequence and assemble. Expression profiling of the avian leptins showed a different pattern compared to mammals, surprisingly with no expression in the adipose tissue. Nevertheless, we found that leptin and leptin receptor (LEPR) genes are expressed in tissues that are implicated in the response to leptin signalling in mammals, and their expression pointed to an autocrine/paracrine instead of endocrine mode of action. Therefore, the role of avian leptins may be relevant to the control of reproduction, stress response and metabolism, but in a different way than in mammals, as it is not the signal of the adipose tissue. Following the success of the identification of chicken leptin and realizing its presence in RNA-seq data rather than genomic assemblies, we asked if other missing genes in chicken could be identified using the RNA-seq approach. Extensive RNA-seq analysis of visceral fat, hypothalamus and pituitary, followed by transcript assembly and bioinformatic analysis led to the identification of 191 novel genes in chickens (BORNELOV et al., 2017), all with very high GC content (~70%). These included the tumour necrosis factor (TNF, also known as TNF-alpha) and nephrin (NPHS1), which were long sought after in the chicken genome. Interestingly, 25 of our novel genes were mapped to the chicken genome in the most recent genome assembly Galgal5. This mapping showed that about 70% of the genes were mapped to chicken microchromosomes, which are known to be GC-rich. The same proportion of genes were mapped in gene blocks demonstrating that the indication that missing genes arranged in gene clusters in other vertebrates cannot be taken as an indication for their evolutionary loss. Altogether, these findings suggest that other missing genes in chickens may be positioned in GC-rich gene clusters in the chicken genome. Among our novel chicken genes, we found five genes belonging to the leptin synteny group in vertebrates: RBM28, SND1, LRRC4, and FLNC. We used this identification to map the chicken leptin gene and its syntenic genes to chicken chromosome 1p (SEROUSSI et al., 2017), thus providing the final proof for the correct identification of the chicken leptin gene. [ABSTRACT FROM AUTHOR]

The serine-/arginine-rich splicing factor 2 (SRSF2) plays pivotal roles in pre-mRNA processing and gene transcription. Recurrent mutations, particularly a proline-to-histidine substitution at position 95 (P95H), are common in neoplastic diseases. Here, we assess SRSF2's diverse functions in squamous cell carcinoma. We show that SRSF2 deletion or homozygous P95H mutation both cause extensive DNA damage leading to cell-cycle arrest. Mechanistically, SRSF2 regulates efficient bi-directional transcription of DNA replication and repair genes, independent from its function in splicing. Further, SRSF2 haploinsufficiency induces DNA damage without halting the cell cycle. Exposing mouse skin to tumor-promoting carcinogens enhances the clonal expansion of heterozygous Srsf2 P95H epidermal cells but unexpectedly inhibits tumor formation. To survive carcinogen treatment, Srsf2 P95H +/- cells undergo substantial transcriptional rewiring and restore bi-directional gene expression. Thus, our study underscores SRSF2's importance in regulating transcription to orchestrate the cell cycle and the DNA damage response., Competing Interests: Declaration of interests The authors declare no competing interests., (Crown Copyright © 2024. Published by Elsevier Inc. All rights reserved.)

PIWI-interacting RNAs (piRNAs) are a class of small non-coding RNAs that are essential for transposon control in animal gonads. In Drosophila ovarian somatic cells, piRNAs are transcribed from large genomic regions called piRNA clusters, which are enriched for transposon fragments and act as a memory of past invasions. Despite being widely present across Drosophila species, somatic piRNA clusters are difficult to identify and study due to their lack of sequence conservation and limited synteny. Current identification methods rely on either extensive manual curation or availability of high-throughput small RNA sequencing data, limiting large-scale comparative studies. We now present FlaHMM, a hidden Markov model developed to automate genomic annotation of flamenco -like unistrand piRNA clusters in Drosophila species, requiring only a genome assembly and transposon annotations. FlaHMM uses transposable element content across 5- or 10-kb bins, which can be calculated from genome sequence alone, and is thus able to detect candidate piRNA clusters without the need to obtain flies and experimentally perform small RNA sequencing. We show that FlaHMM performs on par with piRNA-guided or manual methods, and thus provides a scalable and efficient approach to piRNA cluster annotation in new genome assemblies. FlaHMM is freely available at https://github.com/Hannon-lab/FlaHMM under an MIT licence., (© The Author(s) 2024. Published by Oxford University Press on behalf of NAR Genomics and Bioinformatics.)

The PIWI-interacting RNA (piRNA) pathway prevents endogenous genomic parasites, i.e. transposable elements, from damaging the genetic material of animal gonadal cells. Specific regions in the genome, called piRNA clusters, are thought to define each species' piRNA repertoire and therefore its capacity to recognize and silence specific transposon families. The unistrand cluster flamenco (flam) is essential in the somatic compartment of the Drosophila ovary to restrict Gypsy-family transposons from infecting the neighbouring germ cells. Disruption of flam results in transposon de-repression and sterility, yet it remains unknown whether this silencing mechanism is present more widely. Here, we systematically characterise 119 Drosophila species and identify five additional flam-like clusters separated by up to 45 million years of evolution. Small RNA-sequencing validated these as bona-fide unistrand piRNA clusters expressed in somatic cells of the ovary, where they selectively target transposons of the Gypsy family. Together, our study provides compelling evidence of a widely conserved transposon silencing mechanism that co-evolved with virus-like Gypsy-family transposons., (© 2023. The Author(s).)

Ductal carcinoma in situ (DCIS) is considered a non-invasive precursor to breast cancer, and although associated with an increased risk of developing invasive disease, many women with DCIS will never progress beyond their in situ diagnosis. The path from normal duct to invasive ductal carcinoma (IDC) is not well understood, and efforts to do so are hampered by the substantial heterogeneity that exists between patients, and even within patients. Here we show gene expression analysis from > 2,000 individually micro-dissected ductal lesions representing 145 patients. Combining all samples into one continuous trajectory we show there is a progressive loss in basal layer integrity heading towards IDC, coupled with two epithelial to mesenchymal transitions, one early and a second coinciding with the convergence of DCIS and IDC expression profiles. We identify early processes and potential biomarkers, including CAMK2N1, MNX1, ADCY5, HOXC11 and ANKRD22, whose reduced expression is associated with the progression of DCIS to invasive breast cancer., (© 2022. The Author(s).)

PIWI-interacting RNAs (piRNAs) are small RNAs required to recognize and silence transposable elements. The 5' ends of mature piRNAs are defined through cleavage of long precursor transcripts, primarily by Zucchini (Zuc). Zuc-dependent cleavage typically occurs immediately upstream of a uridine. However, Zuc lacks sequence preference in vitro, pointing towards additional unknown specificity factors. Here, we examine murine piRNAs and reveal a strong and specific enrichment of three sequences (UAA, UAG, UGA)-corresponding to stop codons-at piRNA 5' ends. Stop codon sequences are also enriched immediately after piRNA processing intermediates, reflecting their Zuc-dependent tail-to-head arrangement. Further analyses reveal that a Zuc in vivo cleavage preference at four sequences (UAA, UAG, UGA, UAC) promotes 5' end stop codons. This observation is conserved across mammals and possibly further. Our work provides new insights into Zuc-dependent cleavage and may point to a previously unrecognized connection between piRNA biogenesis and the translational machinery., (© 2022. The Author(s).)

Platelet deficiency, known as thrombocytopenia, can cause hemorrhage and is treated with platelet transfusions. We developed a system for the production of platelet precursor cells, megakaryocytes, from pluripotent stem cells. These cultures can be maintained for >100 days, implying culture renewal by megakaryocyte progenitors (MKPs). However, it is unclear whether the MKP state in vitro mirrors the state in vivo, and MKPs cannot be purified using conventional surface markers. We performed single-cell RNA sequencing throughout in vitro differentiation and mapped each state to its equivalent in vivo. This enabled the identification of five surface markers that reproducibly purify MKPs, allowing us insight into their transcriptional and epigenetic profiles. Last, we performed culture optimization, increasing MKP production. Together, this study has mapped parallels between the MKP states in vivo and in vitro and allowed the purification of MKPs, accelerating the progress of in vitro-derived transfusion products toward the clinic.

Despite four decades of effort, robust propagation of pluripotent stem cells from livestock animals remains challenging. The requirements for self-renewal are unclear and the relationship of cultured stem cells to pluripotent cells resident in the embryo uncertain. Here, we avoided using feeder cells or serum factors to provide a defined culture microenvironment. We show that the combination of activin A, fibroblast growth factor and the Wnt inhibitor XAV939 (AFX) supports establishment and continuous expansion of pluripotent stem cell lines from porcine, ovine and bovine embryos. Germ layer differentiation was evident in teratomas and readily induced in vitro. Global transcriptome analyses highlighted commonality in transcription factor expression across the three species, while global comparison with porcine embryo stages showed proximity to bilaminar disc epiblast. Clonal genetic manipulation and gene targeting were exemplified in porcine stem cells. We further demonstrated that genetically modified AFX stem cells gave rise to cloned porcine foetuses by nuclear transfer. In summary, for major livestock mammals, pluripotent stem cells related to the formative embryonic disc are reliably established using a common and defined signalling environment. This article has an associated 'The people behind the papers' interview., Competing Interests: Competing interests M.K., R.A. and A.S. may benefit from commercial licensing of EDSCs by the University of Nottingham and the University of Cambridge., (© 2021. Published by The Company of Biologists Ltd.)

Color patterns within individual feathers are common in birds but little is known about the genetic mechanisms causing such patterns. Here, we investigate the genetic basis for autosomal barring in chicken, a horizontal striping pattern on individual feathers. Using an informative backcross, we demonstrate that the MC1R locus is strongly associated with this phenotype. A deletion at SOX10, underlying the dark brown phenotype on its own, affects the manifestation of the barring pattern. The coding variant L133Q in MC1R is the most likely causal mutation for autosomal barring in this pedigree. Furthermore, a genetic screen across six different breeds showing different patterning phenotypes revealed that the most striking shared characteristics among these breeds were that they all carried the MC1R alleles Birchen or brown. Our data suggest that the presence of activating MC1R mutations enhancing pigment synthesis is an important mechanism underlying pigmentation patterns on individual feathers in chicken. We propose that MC1R and its antagonist ASIP play a critical role for determining within-feather pigmentation patterns in birds by acting as activator and inhibitor possibly in a Turing reaction-diffusion model., (© 2021 The Authors. Pigment Cell & Melanoma Research published by John Wiley & Sons Ltd.)

DDX3X is a ubiquitously expressed RNA helicase involved in multiple stages of RNA biogenesis. DDX3X is frequently mutated in Burkitt lymphoma, but the functional basis for this is unknown. Here, we show that loss-of-function DDX3X mutations are also enriched in MYC-translocated diffuse large B cell lymphoma and reveal functional cooperation between mutant DDX3X and MYC. DDX3X promotes the translation of mRNA encoding components of the core translational machinery, thereby driving global protein synthesis. Loss-of-function DDX3X mutations moderate MYC-driven global protein synthesis, thereby buffering MYC-induced proteotoxic stress during early lymphomagenesis. Established lymphoma cells restore full protein synthetic capacity by aberrant expression of DDX3Y, a Y chromosome homolog, the expression of which is normally restricted to the testis. These findings show that DDX3X loss of function can buffer MYC-driven proteotoxic stress and highlight the capacity of male B cell lymphomas to then compensate for this loss by ectopic DDX3Y expression., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

Pausing of RNA polymerase II (Pol II) close to promoters is a common regulatory step in RNA synthesis, and is coordinated by a ribonucleoprotein complex scaffolded by the noncoding RNA RN7SK. The function of RN7SK-regulated gene transcription in adult tissue homoeostasis is currently unknown. Here, we deplete RN7SK during mouse and human epidermal stem cell differentiation. Unexpectedly, loss of this small nuclear RNA specifically reduces transcription of numerous cell cycle regulators leading to cell cycle exit and differentiation. Mechanistically, we show that RN7SK is required for efficient transcription of highly expressed gene pairs with bidirectional promoters, which in the epidermis co-regulated cell cycle and chromosome organization. The reduction in transcription involves impaired splicing and RNA decay, but occurs in the absence of chromatin remodelling at promoters and putative enhancers. Thus, RN7SK is directly required for efficient Pol II transcription of highly transcribed bidirectional gene pairs, and thereby exerts tissue-specific functions, such as maintaining a cycling cell population in the epidermis., (© 2021. The Author(s).)

Summary: CONCUR is a standalone tool for codon usage analysis in ribosome profiling experiments. CONCUR uses the aligned reads in BAM format to estimate codon counts at the ribosome E-, P- and A-sites and at flanking positions., Availability and Implementation: CONCUR is written in Perl and is freely available at https://github.com/susbo/concur., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2020. Published by Oxford University Press.)

The nuclear pore complex (NPC) is the principal gateway between nucleus and cytoplasm that enables exchange of macromolecular cargo. Composed of multiple copies of ~30 different nucleoporins (Nups), the NPC acts as a selective portal, interacting with factors which individually license passage of specific cargo classes. Here we show that two Nups of the inner channel, Nup54 and Nup58, are essential for transposon silencing via the PIWI-interacting RNA (piRNA) pathway in the Drosophila ovary. In ovarian follicle cells, loss of Nup54 and Nup58 results in compromised piRNA biogenesis exclusively from the flamenco locus, whereas knockdowns of other NPC subunits have widespread consequences. This provides evidence that some Nups can acquire specialised roles in tissue-specific contexts. Our findings consolidate the idea that the NPC has functions beyond simply constituting a barrier to nuclear/cytoplasmic exchange as genomic loci subjected to strong selective pressure can exploit NPC subunits to facilitate their expression., Competing Interests: MM, VL, AP, SM, SB, IH, MS, GH, BC No competing interests declared, (© 2021, Munafò et al.)

Background: Machine learning involves strategies and algorithms that may assist bioinformatics analyses in terms of data mining and knowledge discovery. In several applications, viz. in Life Sciences, it is often more important to understand how a prediction was obtained rather than knowing what prediction was made. To this end so-called interpretable machine learning has been recently advocated. In this study, we implemented an interpretable machine learning package based on the rough set theory. An important aim of our work was provision of statistical properties of the models and their components., Results: We present the R.ROSETTA package, which is an R wrapper of ROSETTA framework. The original ROSETTA functions have been improved and adapted to the R programming environment. The package allows for building and analyzing non-linear interpretable machine learning models. R.ROSETTA gathers combinatorial statistics via rule-based modelling for accessible and transparent results, well-suited for adoption within the greater scientific community. The package also provides statistics and visualization tools that facilitate minimization of analysis bias and noise. The R.ROSETTA package is freely available at https://github.com/komorowskilab/R.ROSETTA . To illustrate the usage of the package, we applied it to a transcriptome dataset from an autism case-control study. Our tool provided hypotheses for potential co-predictive mechanisms among features that discerned phenotype classes. These co-predictors represented neurodevelopmental and autism-related genes., Conclusions: R.ROSETTA provides new insights for interpretable machine learning analyses and knowledge-based systems. We demonstrated that our package facilitated detection of dependencies for autism-related genes. Although the sample application of R.ROSETTA illustrates transcriptome data analysis, the package can be used to analyze any data organized in decision tables.

In animal gonads, the PIWI-interacting RNA (piRNA) pathway guards genome integrity in part through the co-transcriptional gene silencing of transposon insertions. In Drosophila ovaries, piRNA-loaded Piwi detects nascent transposon transcripts and instructs heterochromatin formation through the P anoramix- i nduced c o- t ranscriptional s ilencing (PICTS) complex, containing Panoramix, Nxf2 and Nxt1. Here, we report that the highly conserved dynein light chain LC8/Cut-up (Ctp) is an essential component of the PICTS complex. Loss of Ctp results in transposon de-repression and a reduction in repressive chromatin marks specifically at transposon loci. In turn, Ctp can enforce transcriptional silencing when artificially recruited to RNA and DNA reporters. We show that Ctp drives dimerisation of the PICTS complex through its interaction with conserved motifs within Panoramix. Artificial dimerisation of Panoramix bypasses the necessity for its interaction with Ctp, demonstrating that conscription of a protein from a ubiquitous cellular machinery has fulfilled a fundamental requirement for a transposon silencing complex., Competing Interests: EE, KJ, SB, MM, VF, EK, EB, BC, GH No competing interests declared, (© 2021, Eastwood et al.)

The highly abundant N6-methyladenosine (m6A) RNA modification affects most aspects of mRNA function, yet the precise function of the rarer 5-methylcytidine (m5C) remains largely unknown. Here, we map m5C in the human transcriptome using methylation-dependent individual-nucleotide resolution cross-linking and immunoprecipitation (miCLIP) combined with RNA bisulfite sequencing. We identify NSUN6 as a methyltransferase with strong substrate specificity towards mRNA. NSUN6 primarily targeted three prime untranslated regions (3'UTR) at the consensus sequence motif CTCCA, located in loops of hairpin structures. Knockout and rescue experiments revealed enhanced mRNA and translation levels when NSUN6-targeted mRNAs were methylated. Ribosome profiling further demonstrated that NSUN6-specific methylation correlated with translation termination. While NSUN6 was dispensable for mouse embryonic development, it was down-regulated in human tumours and high expression of NSUN6 indicated better patient outcome of certain cancer types. In summary, our study identifies NSUN6 as a methyltransferase targeting mRNA, potentially as part of a quality control mechanism involved in translation termination fidelity., (© The Author(s) 2020. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Posttranscriptional modifications in transfer RNA (tRNA) are often critical for normal development because they adapt protein synthesis rates to a dynamically changing microenvironment. However, the precise cellular mechanisms linking the extrinsic stimulus to the intrinsic RNA modification pathways remain largely unclear. Here, we identified the cytosine-5 RNA methyltransferase NSUN2 as a sensor for external stress stimuli. Exposure to oxidative stress efficiently repressed NSUN2, causing a reduction of methylation at specific tRNA sites. Using metabolic profiling, we showed that loss of tRNA methylation captured cells in a distinct catabolic state. Mechanistically, loss of NSUN2 altered the biogenesis of tRNA-derived noncoding fragments (tRFs) in response to stress, leading to impaired regulation of protein synthesis. The intracellular accumulation of a specific subset of tRFs correlated with the dynamic repression of global protein synthesis. Finally, NSUN2-driven RNA methylation was functionally required to adapt cell cycle progression to the early stress response. In summary, we revealed that changes in tRNA methylation profiles were sufficient to specify cellular metabolic states and efficiently adapt protein synthesis rates to cell stress., Competing Interests: I have read the journal's policy and the authors of this manuscript have the following competing interests: MF is consulting for Storm Therapeutics.

The presence and absence of RNA modifications regulates RNA metabolism by modulating the binding of writer, reader, and eraser proteins. For 5-methylcytosine (m 5 C) however, it is largely unknown how it recruits or repels RNA-binding proteins. Here, we decipher the consequences of m 5 C deposition into the abundant non-coding vault RNA VTRNA1.1. Methylation of cytosine 69 in VTRNA1.1 occurs frequently in human cells, is exclusively mediated by NSUN2, and determines the processing of VTRNA1.1 into small-vault RNAs (svRNAs). We identify the serine/arginine rich splicing factor 2 (SRSF2) as a novel VTRNA1.1-binding protein that counteracts VTRNA1.1 processing by binding the non-methylated form with higher affinity. Both NSUN2 and SRSF2 orchestrate the production of distinct svRNAs. Finally, we discover a functional role of svRNAs in regulating the epidermal differentiation programme. Thus, our data reveal a direct role for m 5 C in the processing of VTRNA1.1 that involves SRSF2 and is crucial for efficient cellular differentiation.

Background: The uneven use of synonymous codons in the transcriptome regulates the efficiency and fidelity of protein translation rates. Yet, the importance of this codon bias in regulating cell state-specific expression programmes is currently debated. Here, we ask whether different codon usage controls gene expression programmes in self-renewing and differentiating embryonic stem cells., Results: Using ribosome and transcriptome profiling, we identify distinct codon signatures during human embryonic stem cell differentiation. We find that cell state-specific codon bias is determined by the guanine-cytosine (GC) content of differentially expressed genes. By measuring the codon frequencies at the ribosome active sites interacting with transfer RNAs (tRNA), we further discover that self-renewing cells optimize translation of codons that depend on the inosine tRNA modification in the anticodon wobble position. Accordingly, inosine levels are highest in human pluripotent embryonic stem cells. This effect is conserved in mice and is independent of the differentiation stimulus., Conclusions: We show that GC content influences cell state-specific mRNA levels, and we reveal how translational mechanisms based on tRNA modifications change codon usage in embryonic stem cells.

The gene regulatory network (GRN) of naive mouse embryonic stem cells (ESCs) must be reconfigured to enable lineage commitment. TCF3 sanctions rewiring by suppressing components of the ESC transcription factor circuitry. However, TCF3 depletion only delays and does not prevent transition to formative pluripotency. Here, we delineate additional contributions of the ETS-family transcription factor ETV5 and the repressor RBPJ. In response to ERK signaling, ETV5 switches activity from supporting self-renewal and undergoes genome relocation linked to commissioning of enhancers activated in formative epiblast. Independent upregulation of RBPJ prevents re-expression of potent naive factors, TBX3 and NANOG, to secure exit from the naive state. Triple deletion of Etv5, Rbpj, and Tcf3 disables ESCs, such that they remain largely undifferentiated and locked in self-renewal, even in the presence of differentiation stimuli. Thus, genetic elimination of three complementary drivers of network transition stalls developmental progression, emulating environmental insulation by small-molecule inhibitors., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)

Chromatin remodeling complexes play essential roles in metazoan development through widespread control of gene expression, but the precise molecular mechanisms by which they do this in vivo remain ill defined. Using an inducible system with fine temporal resolution, we show that the nucleosome remodeling and deacetylation (NuRD) complex controls chromatin architecture and the protein binding repertoire at regulatory regions during cell state transitions. This is primarily exerted through its nucleosome remodeling activity while deacetylation at H3K27 follows changes in gene expression. Additionally, NuRD activity influences association of RNA polymerase II at transcription start sites and subsequent nascent transcript production, thereby guiding the establishment of lineage-appropriate transcriptional programs. These findings provide a detailed molecular picture of genome-wide modulation of lineage-specific transcription by an essential chromatin remodeling complex as well as insight into the orchestration of molecular events involved in transcriptional transitions in vivo. VIDEO ABSTRACT., (Copyright © 2018 The Author(s). Published by Elsevier Inc. All rights reserved.)

Background: The mammalian adipose tissue plays a central role in energy-balance control, whereas the avian visceral fat hardly expresses leptin, the key adipokine in mammals. Therefore, to assess the endocrine role of adipose tissue in birds, we compared the transcriptome and proteome between two metabolically different types of chickens, broilers and layers, bred towards efficient meat and egg production, respectively., Results: Broilers and layer hens, grown up to sexual maturation under free-feeding conditions, differed 4.0-fold in weight and 1.6-fold in ovarian-follicle counts, yet the relative accumulation of visceral fat was comparable. RNA-seq and mass-spectrometry (MS) analyses of visceral fat revealed differentially expressed genes between broilers and layers, 1106 at the mRNA level (FDR ≤ 0.05), and 203 at the protein level (P ≤ 0.05). In broilers, Ingenuity Pathway Analysis revealed activation of the PTEN-pathway, and in layers increased response to external signals. The expression pattern of genes encoding fat-secreted proteins in broilers and layers was characterized in the RNA-seq and MS data, as well as by qPCR on visceral fat under free feeding and 24 h-feed deprivation. This characterization was expanded using available RNA-seq data of tissues from red junglefowl, and of visceral fat from broilers of different types. These comparisons revealed expression of new adipokines and secreted proteins (LCAT, LECT2, SERPINE2, SFTP1, ZP1, ZP3, APOV1, VTG1 and VTG2) at the mRNA and/or protein levels, with dynamic gene expression patterns in the selected chicken lines (except for ZP1; FDR/P ≤ 0.05) and feed deprivation (NAMPT, SFTPA1 and ZP3) (P ≤ 0.05). In contrast, some of the most prominent adipokines in mammals, leptin, TNF, IFNG, and IL6 were expressed at a low level (FPKM/RPKM< 1) and did not show differential mRNA expression neither between broiler and layer lines nor between fed vs. feed-deprived chickens., Conclusions: Our study revealed that RNA and protein expression in visceral fat changes with selective breeding, suggesting endocrine roles of visceral fat in the selected phenotypes. In comparison to gene expression in visceral fat of mammals, our findings points to a more direct cross talk of the chicken visceral fat with the reproductive system and lower involvement in the regulation of appetite, inflammation and insulin resistance.

Correction: After the publication of this work [1] an error was noticed in one of the author surnames. The author name Leif Anderson should be spelt as Leif Andersson.

Background: Misidentification of the chicken leptin gene has hampered research of leptin signaling in this species for almost two decades. Recently, the genuine leptin gene with a GC-rich (~70%) repetitive-sequence content was identified in the chicken genome but without indicating its genomic position. This suggests that such GC-rich sequences are difficult to sequence and therefore substantial regions are missing from the current chicken genome assembly., Results: A radiation hybrid panel of chicken-hamster Wg3hCl2 cells was used to map the genome location of the chicken leptin gene. Contrary to our expectations, based on comparative genome mapping and sequence characteristics, the chicken leptin was not located on a microchromosome, which are known to contain GC-rich and repetitive regions, but at the distal tip of the largest chromosome (1p). Following conserved synteny with other vertebrates, we also mapped five additional genes to this genomic region (ARF5, SND1, LRRC4, RBM28, and FLNC), bridging the genomic gap in the current Galgal5 build for this chromosome region. All of the short scaffolds containing these genes were found to consist of GC-rich (54 to 65%) sequences comparing to the average GC-content of 40% on chromosome 1. In this syntenic group, the RNA-binding protein 28 (RBM28) was in closest proximity to leptin. We deduced the full-length of the RBM28 cDNA sequence and profiled its expression patterns detecting a negative correlation (R = - 0.7) between the expression of leptin and of RBM28 across tissues that expressed at least one of the genes above the average level. This observation suggested a local regulatory interaction between these genes. In adipose tissues, we observed a significant increase in RBM28 mRNA expression in breeds with lean phenotypes., Conclusion: Mapping chicken leptin together with a cluster of five syntenic genes provided the final proof for its identification as the true chicken ortholog. The high GC-content observed for the chicken leptin syntenic group suggests that other similar clusters of genes in GC-rich genomic regions are missing from the current genome assembly (Galgal5), which should be resolved in future assemblies of the chicken genome.

Through RNA-Seq analyses, we identified 137 genes that are missing in chicken, including the long-sought-after nephrin and tumor necrosis factor genes. These genes tended to cluster in GC-rich regions that have poor coverage in genome sequence databases. Hence, the occurrence of syntenic groups of vertebrate genes that have not been observed in Aves does not prove the evolutionary loss of such genes.Please see related Research article: http://dx.doi.org/10.1186/s13059-014-0565-1 and Please see response from Lovell et al: https://www.dx.doi.org/10.1186/s13059-017-1234-y.

Defensins constitute an evolutionary conserved family of cationic antimicrobial peptides that play a key role in host innate immune responses to infection. Defensin genes generally reside in complex genomic regions that are prone to structural variation, and defensin genes exhibit extensive copy number variation in humans and in other species. Copy number variation of defensin genes was examined in inbred lines of Leghorn and Fayoumi chickens, and a duplication of defensin7 was discovered in the Fayoumi breed. Analysis of junction sequences confirmed the occurrence of a simple tandem duplication of defensin7 with sequence identity at the junction, suggesting nonallelic homologous recombination between defensin7 and defensin6 The duplication event generated two chimeric promoters that are best explained by gene conversion followed by homologous recombination. Expression of defensin7 was not elevated in animals with two genes despite both genes being transcribed in the tissues examined. Computational prediction of promoter regions revealed the presence of several putative transcription factor binding sites generated by the duplication event. These data provide insight into the evolution and possible function of large gene families and specifically, the defensins., Competing Interests: The authors declare no conflict of interest.

More than 20 years after characterization of the key regulator of mammalian energy balance, leptin, we identified the leptin (LEP) genes of chicken (Gallus gallus) and duck (Anas platyrhynchos). The extreme guanine-cytosine content (∼70%), the location in a genomic region with low-complexity repetitive and palindromic sequence elements, the relatively low sequence conservation, and low level of expression have hampered the identification of these genes until now. In vitro-expressed chicken and duck leptins specifically activated signaling through the chicken leptin receptor in cell culture. In situ hybridization demonstrated expression of LEP mRNA in granular and Purkinje cells of the cerebellum, anterior pituitary, and in embryonic limb buds, somites, and branchial arches, suggesting roles in adult brain control of energy balance and during embryonic development. The expression patterns of LEP and the leptin receptor (LEPR) were explored in chicken, duck, and quail (Coturnix japonica) using RNA-sequencing experiments available in the Short Read Archive and by quantitative RT-PCR. In adipose tissue, LEP and LEPR were scarcely transcribed, and the expression level was not correlated to adiposity. Our identification of the leptin genes in chicken and duck genomes resolves a long lasting controversy regarding the existence of leptin genes in these species. This identification was confirmed by sequence and structural similarity, conserved exon-intron boundaries, detection in numerous genomic, and transcriptomic datasets and characterization by PCR, quantitative RT-PCR, in situ hybridization, and bioassays. Our results point to an autocrine/paracrine mode of action for bird leptin instead of being a circulating hormone as in mammals.

Background: Several post-translational histone modifications are mainly found in gene promoters and are associated with the promoter activity. It has been hypothesized that histone modifications regulate the transcription, as opposed to the traditional view with transcription factors as the key regulators. Promoters of most active genes do not only initiate transcription of the coding sequence, but also a substantial amount of transcription of the antisense strand upstream of the transcription start site (TSS). This promoter feature has generally not been considered in previous studies of histone modifications and transcription factor binding., Results: We annotated protein-coding genes as bi- or unidirectional depending on their mode of transcription and compared histone modifications and transcription factor occurrences between them. We found that H3K4me3, H3K9ac, and H3K27ac were significantly more enriched upstream of the TSS in bidirectional genes compared with the unidirectional ones. In contrast, the downstream histone modification signals were similar, suggesting that the upstream histone modifications might be a consequence of transcription rather than a cause. Notably, we found well-positioned CTCF and RAD21 peaks approximately 60-80 bp upstream of the TSS in the unidirectional genes. The peak heights were related to the amount of antisense transcription and we hypothesized that CTCF and cohesin act as a barrier against antisense transcription., Conclusions: Our results provide insights into the distribution of histone modifications at promoters and suggest a novel role of CTCF and cohesin as regulators of transcriptional direction.

Background: The use of classification algorithms is becoming increasingly important for the field of computational biology. However, not only the quality of the classification, but also its biological interpretation is important. This interpretation may be eased if interacting elements can be identified and visualized, something that requires appropriate tools and methods., Results: We developed a new approach to detecting interactions in complex systems based on classification. Using rule-based classifiers, we previously proposed a rule network visualization strategy that may be applied as a heuristic for finding interactions. We now complement this work with Ciruvis, a web-based tool for the construction of rule networks from classifiers made of IF-THEN rules. Simulated and biological data served as an illustration of how the tool may be used to visualize and interpret classifiers. Furthermore, we used the rule networks to identify feature interactions, compared them to alternative methods, and computationally validated the findings., Conclusions: Rule networks enable a fast method for model visualization and provide an exploratory heuristic to interaction detection. The tool is made freely available on the web and may thus be used to aid and improve rule-based classification.

Both genetic and environmental factors are important for the development of allergic diseases. However, a detailed understanding of how such factors act together is lacking. To elucidate the interplay between genetic and environmental factors in allergic diseases, we used a novel bioinformatics approach that combines feature selection and machine learning. In two materials, PARSIFAL (a European cross-sectional study of 3113 children) and BAMSE (a Swedish birth-cohort including 2033 children), genetic variants as well as environmental and lifestyle factors were evaluated for their contribution to allergic phenotypes. Monte Carlo feature selection and rule based models were used to identify and rank rules describing how combinations of genetic and environmental factors affect the risk of allergic diseases. Novel interactions between genes were suggested and replicated, such as between ORMDL3 and RORA, where certain genotype combinations gave odds ratios for current asthma of 2.1 (95% CI 1.2-3.6) and 3.2 (95% CI 2.0-5.0) in the BAMSE and PARSIFAL children, respectively. Several combinations of environmental factors appeared to be important for the development of allergic disease in children. For example, use of baby formula and antibiotics early in life was associated with an odds ratio of 7.4 (95% CI 4.5-12.0) of developing asthma. Furthermore, genetic variants together with environmental factors seemed to play a role for allergic diseases, such as the use of antibiotics early in life and COL29A1 variants for asthma, and farm living and NPSR1 variants for allergic eczema. Overall, combinations of environmental and life style factors appeared more frequently in the models than combinations solely involving genes. In conclusion, a new bioinformatics approach is described for analyzing complex data, including extensive genetic and environmental information. Interactions identified with this approach could provide useful hints for further in-depth studies of etiological mechanisms and may also strengthen the basis for risk assessment and prevention.

Background: Steatohepatitis occurs in alcoholic liver disease and may progress to liver cirrhosis and hepatocellular carcinoma. Its molecular pathogenesis is to a large degree unknown. Histone modifications play a key role in transcriptional regulations as marks for silencing and activation of gene expression and as marks for functional elements. Many transcription factors (TFs) are crucial for the control of the genes involved in metabolism, and abnormality in their function may lead to disease., Methods: We performed ChIP-seq of the histone modifications H3K4me1, H3K4me3 and H3K27ac and a candidate transcription factor (USF1) in liver tissue from patients with steatohepatitis and normal livers and correlated results to mRNA-expression and genotypes., Results: We found several regions that are differentially enriched for histone modifications between disease and normal tissue, and qRT-PCR results indicated that the expression of the tested genes strongly correlated with differential enrichment of histone modifications but is independent of USF1 enrichment. By gene ontology analysis of differentially modified genes we found many disease associated genes, some of which had previously been implicated in the etiology of steatohepatitis. Importantly, the genes associated to the strongest histone peaks in the patient were over-represented in cancer specific pathways suggesting that the tissue was on a path to develop to cancer, a common complication to the disease. We also found several novel SNPs and GWAS catalogue SNPs that are candidates to be functional and therefore needs further study., Conclusion: In summary we find that analysis of chromatin features in tissue samples provides insight into disease mechanisms.

Splicing is a complex process regulated by sequence at the classical splice sites and other motifs in exons and introns with an enhancing or silencing effect. In addition, specific histone modifications on nucleosomes positioned over the exons have been shown to correlate both positively and negatively with exon expression. Here, we trained a model of "IF … THEN …" rules to predict exon inclusion levels in a transcript from histone modification patterns. Furthermore, we showed that combinations of histone modifications, in particular those residing on nucleosomes preceding or succeeding the exon, are better predictors of exon inclusion levels than single modifications. The resulting model was evaluated with cross validation and had an average accuracy of 72% for 27% of the exons, which demonstrates that epigenetic signals substantially mark alternative splicing.