The cis-trans binding strength defined by motif frequencies facilitates statistical inference of transcriptional regulation

Research Note: Identification of core promoter region of the polyunsaturated fatty acid synthesis-related gene family in chicken

Chicken is considered an ideal model species to study the synthesis of polyunsaturated fatty acids (PUFAs) due to its appropriate proportions of fatty acids and abundant content of PUFAs, suitable for human consumption. However, the molecular mechanisms regulating poultry PUFA synthesis remain unclear. Here, we systematically explored the transcriptional regulation activity of the gene family related to PUFA synthesis in chicken by carrying out the Dual-Luciferase Reporter Assay. We identified the core promoter regions of members of the chicken PUFA synthesis-related gene family, including ELOVL1, ELOVL2, ELOVL3, ELOVL4, ELOVL5, ELOVL6, ELOVL7, FADS1, FADS2, FADS6, SCD, and SCD5. Additionally, changes in relative fluorescence values of different truncated segments in the upstream regulatory region of these genes indicate the existence of regulatory regions. Furthermore, we predicted the transcription factors that bind to the identified core promoter regions of multiple genes, including Sp1, NF-1, C/EBPalpha, etc. These findings provide a basis for the molecular mechanisms regulating poultry PUFA synthesis and offer new scientific insight into the potential improvement of poultry meat quality in the future.

MUREN: a robust and multi-reference approach of RNA-seq transcript normalization

Background

Normalization of RNA-seq data aims at identifying biological expression differentiation between samples by removing the effects of unwanted confounding factors. Explicitly or implicitly, the justification of normalization requires a set of housekeeping genes. However, the existence of housekeeping genes common for a very large collection of samples, especially under a wide range of conditions, is questionable.

Results

We propose to carry out pairwise normalization with respect to multiple references, selected from representative samples. Then the pairwise intermediates are integrated based on a linear model that adjusts the reference effects. Motivated by the notion of housekeeping genes and their statistical counterparts, we adopt the robust least trimmed squares regression in pairwise normalization. The proposed method (MUREN) is compared with other existing tools on some standard data sets. The goodness of normalization emphasizes on preserving possible asymmetric differentiation, whose biological significance is exemplified by a single cell data of cell cycle. MUREN is implemented as an R package. The code under license GPL-3 is available on the github platform: github.com/hippo-yf/MUREN and on the conda platform: anaconda.org/hippo-yf/r-muren.

Conclusions

MUREN performs the RNA-seq normalization using a two-step statistical regression induced from a general principle. We propose that the densities of pairwise differentiations are used to evaluate the goodness of normalization. MUREN adjusts the mode of differentiation toward zero while preserving the skewness due to biological asymmetric differentiation. Moreover, by robustly integrating pre-normalized counts with respect to multiple references, MUREN is immune to individual outlier samples.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04288-0.

Dual Eigen-modules of Cis-Element Regulation Profiles and Selection of Cognition-Language Eigen-direction along Evolution in Hominidae

To understand the genomic basis accounting for the phenotypic differences between human and apes, we compare the matrices consisting of the cis-element frequencies in the proximal regulatory regions of their genomes. One such frequency matrix is represented by a robust singular value decomposition. For each singular value, the negative and positive ends of the sorted motif eigenvector correspond to the dual ends of the sorted gene eigenvector, respectively, comprising a dual eigen-module defined by cis-regulatory element frequencies (CREF). The CREF eigen-modules at levels 1, 2, 3, and 6 are highly conserved across humans, chimpanzees, and orangutans. The key biological processes embedded in the top three CREF eigen-modules are reproduction versus embryogenesis, fetal maturation versus immune system, and stress responses versus mitosis. Although the divergence at the nucleotide level between the chimpanzee and human genome was small, their cis-element frequency matrices crossed a singularity point, at which the fourth and fifth singular values were identical. The CREF eigen-modules corresponding to the fourth and fifth singular values were reorganized along the evolution from apes to human. Interestingly, the fourth sorted gene eigenvector encodes the phenotypes unique to human such as long-term memory, language development, and social behavior. The number of motifs present on Alu elements increases substantially at the fourth level. The motif analysis together with the cases of human-specific Alu insertions suggests that mutations related to Alu elements play a critical role in the evolution of the human-phenotypic gene eigenvector.