RNA-Seq analysis reveals expression regulatory divergence of W-linked genes between two contrasting chicken breeds

The regulation of gene expression is a complex process involving organism function and phenotypic diversity, and is caused by cis- and trans- regulation. While prior studies identified the regulatory pattern of the autosome rewiring in hybrids, the role of gene regulation in W sex chromosomes is not clear due to their degradation and sex-limit expression. Here, we developed reciprocal crosses of two chicken breeds, White Leghorn and Cornish Game, which exhibited broad differences of gender-related traits, and assessed the expression of the genes on W chromosome to disentangle the contribution of cis- and trans-factors to expression divergence. We found that there was not appear to be an association between female fecundity and W chromosome gene expression, that 44% of expressed genes had divergent expression between breeds in both tissues, with only 17% of them showing greater expression in White Leghorn. We observed that the proportion of trans-acting genes in W chromosome was higher than cis-regulatory divergence. There were most parental divergence expression genes in muscle, also more heterosis compared with other two tissues. A strong dominant impact of Cornish alleles in brain, while obvious crosses-specific regulatory patterns appeared in liver. Taken together, this work describes the regulatory divergence of W-linked genes between two contrasting breeds and indicates sex chromosomes have a unique regulation and expression mechanism.

Effects of homoeologous exchange on gene expression and alternative splicing in a newly formed allotetraploid wheat

Homoeologous exchange (HE) is a major mechanism generating post-polyploidization genetic variation with important evolutionary consequences. However, the direct impacts of HE without entangling with additional evolutionary forces on gene expression remains to be fully understood. Here, we analyzed high-throughput RNA-seq data of young leaves from individuals of a synthetic allotetraploid wheat (AADD), which contain variable numbers of HEs. We aimed to investigate if and to which extent HE directly impacts gene expression and alternative splicing (AS). We found that HE impacts expression of genes located within HE regions primarily via cis-acting dosage effect, which led to significant changes in the total expression of homoeolog pairs, especially for homoeologs whose original expression was biased. In parallel, HE influences expression of a large amount of genes residing in non-HE regions by trans-regulation leading to convergent expression of homoeolog pairs. Intriguingly, when taking into account of the original relative homoeolog expression states, homoeolog pairs under trans-effect are more prone to showing convergent response to HE whereas those under cis-effect trended to show subgenome-specific expression. Moreover, HE induced quantitative, largely individual-specific, changes of alternative splicing (AS) events. Like homoeologs expression, homoeo-AS events which related to trans effect were more responsive to HE. HE therefore exerts multifaceted immediate effects on gene expression and, to a less extent, also transcript diversity in nascent allopolyploidy.

Comprehensive Analysis of Differential Expression Profiles of Long Non-coding RNAs with Associated Co-expression and Competing Endogenous RNA Networks in the Hippocampus of Patients with Alzheimer's Disease

Background: Alzheimer's disease (AD) is still one of the major threats to human health. Although a satisfactory treatment for AD has not yet been discovered, it is necessary to continue to search for novel approaches to deal with this insidious and debilitating disease. Although numerous studies have shown that long non-coding RNA (lncRNA) occupy a significant role in a variety of diseases, their roles in AD remain unclear. Objectives: Using data analysis to explore the role of lncRNA in the course of AD, to further our understanding of AD, and to look forward to finding a new breakthrough for the treatment of AD. Methods: We downloaded and screened expression data of the hippocampal regions of patients with AD from the Gene Expression Omnibus database. We generated lncRNA-miRNA-mRNA networks based on the competing endogenous RNA (ceRNA) hypothesis, and according to gene expression level, we constructed a coding-noncoding co-expression (CNC) network and then executed cis- and trans-regulation analyses. Results: Through comprehensive and systematic analyses, we found that lncRNAs MALAT1, OIP5-AS1, LINC00657, and lnc-NUMB-1 regulated the expression of the key AD pathogenic genes APP, PSEN1, BACE1; and that these lncRNAs may promote the distribution of β-amyloid (Aβ protein) in the brain through exosomes. In addition, lncRNAs were found to adjust viral transcriptional expression, thereby further supporting viral pathogenesis for AD. Conclusions: The lncRNAs MALAT1, OIP5-AS1, LINC00657, and lnc-NUMB-1 that are present in the hippocampus of AD patients exert an important influence on the development of this disease.

seRNA PAM-1 regulates skeletal muscle satellite cell activation and aging through trans regulation of Timp2 expression synergistically with Ddx5

Muscle satellite cells (SCs) are responsible for muscle homeostasis and regeneration; and lncRNAs play important roles in regulating SC activities. Here in this study, we identify PAM-1 (Pax7 Associated Muscle lncRNA) that is induced in activated SCs to promote SC activation into myoblast cells upon injury. PAM-1 is generated from a myoblast specific super-enhancer (SE); as a seRNA it binds with a number of target genomic loci predominantly in trans. Further studies demonstrate that it interacts with Ddx5 to tether PAM-1 SE to it inter-chromosomal targets Timp2 and Vim to activate the gene expression. Lastly, we show that PAM-1 expression is increased in aging SCs, which leads to enhanced inter-chromosomal interaction and target genes up-regulation. Altogether, our findings identify PAM-1 as a previously unknown lncRNA that regulates both SC activation and aging through its trans gene regulatory activity.

Genetic studies of human-chimpanzee divergence using stem cell fusions

Complete genome sequencing has identified millions of DNA changes that differ between humans and chimpanzees. Although a subset of these changes likely underlies important phenotypic differences between humans and chimpanzees, it is currently difficult to distinguish causal from incidental changes and to map specific phenotypes to particular genome locations. To facilitate further genetic study of human-chimpanzee divergence, we have generated human and chimpanzee auto-tetraploids and allo-tetraploids by fusing induced pluripotent stem cells (iPSCs) of each species. The resulting tetraploid iPSCs can be stably maintained and retain the ability to differentiate along ectoderm, mesoderm, and endoderm lineages. RNA sequencing identifies thousands of genes whose expression differs between humans and chimpanzees when assessed in single-species diploid or auto-tetraploid iPSCs. Analysis of gene expression patterns in inter-specific allo-tetraploid iPSCs shows that human-chimpanzee expression differences arise from substantial contributions of both cis-acting changes linked to the genes themselves, and trans-acting changes elsewhere in the genome. To enable further genetic mapping of species differences, we tested chemical treatments for stimulating genome-wide mitotic recombination between human and chimpanzee chromosomes, and CRISPR methods for inducing species-specific changes on particular chromosomes in allo-tetraploid cells. We successfully generated derivative cells with nested deletions or inter-specific recombination on the X chromosome. These studies identify a long distance cis-regulatory domain of the Fragile X-associated gene (FMR1), confirm an important role for the X chromosome in trans-regulation of other expression differences, and illustrate the potential of this system for more detailed mapping of the molecular basis of human and chimpanzee evolution.

Alternative splicing of viral transcripts: the dark side of the HBV

Regulation of alternative splicing is one of the most efficient mechanisms to enlarge the proteomic diversity in eukaryotic organisms. Many viruses hijack the splicing machinery following infection to accomplish their replication cycle. Regarding the HBV, numerous reports have described alternative splicing events of the long viral transcript (pregenomic RNA), which also acts as a template for viral genome replication. Alternative splicing of HBV pregenomic RNAs allows the synthesis of at least 20 spliced variants. In addition, almost all these spliced forms give rise to defective particles, detected in the blood of infected patients. HBV-spliced RNAs have long been unconsidered, probably due to their uneasy detection in comparison to unspliced forms as well as for their dispensable role during viral replication. However, recent data highlighted the relevance of these HBV-spliced variants through (1) the trans-regulation of the alternative splicing of viral transcripts along the course of liver disease; (2) the ability to generate defective particle formation, putative biomarker of the liver disease progression; (3) modulation of viral replication; and (4) their intrinsic propensity to encode for novel viral proteins involved in liver pathogenesis and immune response. Altogether, tricky regulation of HBV alternative splicing may contribute to modulate multiple viral and cellular processes all along the course of HBV-related liver disease.

Parallel and Population-specific Gene Regulatory Evolution in Cold-Adapted Fly Populations

Changes in gene regulation at multiple levels may comprise an important share of the molecular changes underlying adaptive evolution in nature. However, few studies have assayed within- and between-population variation in gene regulatory traits at a transcriptomic scale, and therefore inferences about the characteristics of adaptive regulatory changes have been elusive. Here, we assess quantitative trait differentiation in gene expression levels and alternative splicing (intron usage) between three closely-related pairs of natural populations of Drosophila melanogaster from contrasting thermal environments that reflect three separate instances of cold tolerance evolution. The cold-adapted populations were known to show population genetic evidence for parallel evolution at the SNP level, and here we find evidence for parallel expression evolution between them, with stronger parallelism at larval and adult stages than for pupae. We also implement a flexible method to estimate cis- versus trans-encoded contributions to expression or splicing differences at the adult stage. The apparent contributions of cis- versus trans-regulation to adaptive evolution vary substantially among population pairs. While two of three population pairs show a greater enrichment of cis-regulatory differences among adaptation candidates, trans-regulatory differences are more likely to be implicated in parallel expression changes between population pairs. Genes with significant cis-effects are enriched for signals of elevated genetic differentiation between cold- and warm-adapted populations, suggesting that they are potential targets of local adaptation. These findings expand our knowledge of adaptive gene regulatory evolution and our ability to make inferences about this important and widespread process.

Chromatin-interacting transposon RNAs linking to the core trans-inhibition circuitry for embryonic stem cell identity

Transposable DNA sequences constitute more than half of the human and mouse genomes. A large number of non-coding RNAs, named as transposon RNAs, are derived from these transposable elements. The cis-regulatory function of transposable DNA elements, such as LINE-1 and Alu has been largely explored. But the biological roles of transposon RNAs aren't well understood. Here, investigations of RNA-chromatin interactions provide us with comprehensive evidence that specific families of transposon RNAs play roles in trans-regulation linking to the core inhibition circuitry for embryonic stem cell identity. Alternative modes of the RNA-DNA hybrid duplex and protein-recruited RNA scaffold are required for the regulatory activities of transposon RNAs. LINE-1 RNAs co-locating with KAP1 form a negative feedback loop stabilizing the transcription of LINE-1 DNA elements via RNA-DNA hybrids. In another way LINE-1 RNAs, together with the reprogramming three factors and Polycomb repressive complexes, participate in the inhibition on dozens of differentiation-relative genes.

Sequence and phylogenetic analysis revealed structurally conserved domains and motifs in lincRNA-p21

Long Intergenic Non-coding RNAs (lincRNAs) are the largest class of long non-coding RNAs in the eukaryotes, which originate from the intergenic regions of the genome. A ~4kb long lincRNA-p21 is derived from a transcription unit next to the p21/Cdkn1a gene locus. LincRNA-p21 plays key regulatory roles in p53 dependent transcriptional repression and translational repression through its physical association with proteins such as hnRNP-K and HuR. It is also involved in the aberrant gene expression in different cancers. However, detailed information on its structure, recognition, and trans-regulation by proteins is not well known. In this study, we have carried out a complete gene analysis and annotation of lincRNA-p21. This analysis showed that lincRNA-p21 is highly conserved in primates, and its conservation drops significantly in lower organisms. Furthermore, our analysis has revealed two structurally conserved domains in the 5’ and 3’ terminal regions of lincRNA-p21. Phylogenetic analysis has revealed discrete evolutionary dynamics in these conserved domains for orthologous sequences of lincRNA-p21, which have evolved slowly across primates compared to other mammals. Using Infernal based covariance analysis, we have computed the secondary structures of these domains. The secondary structures were further validated by energy minimization criteria for individual orthologous sequences as well as the full-length human lincRNA-p21. In summary, this analysis has led to the identification of sequence and structural motifs in the conserved fragments, indicating the functional importance for these regions.

Transcriptome-wide association study of risk of recurrence in Black and White breast cancer patients

Background: Continuous risk of recurrence scores (CRS) based on PAM50 gene expression are vital prognostic tools for breast cancer (BC). Studies have shown that Black women (BW) have higher CRS than White women (WW). Although systemic injustices contribute substantially to BC disparities, evidence for biological and germline contributions is also emerging. We investigated germline genetic associations with CRS and CRS disparity through a Transcriptome-Wide Association Study (TWAS). Methods: In the Carolina Breast Cancer Study, using race-specific predictive models of tumor expression from germline genetics, we performed race-stratified (N=1,043 WW, 1083 BW) linear regressions of three CRS (ROR-S: PAM50 subtype score; Proliferation Score; ROR-P: ROR-S plus Proliferation Score) on imputed Genetically-Regulated tumor eXpression (GReX). Using Bayesian multivariate regression and adaptive shrinkage, we tested TWAS-significant genes for associations with PAM50 tumor expression and subtype to elucidate patterns of germline regulation underlying TWAS-gene and CRS associations. Results: At FDR-adjusted P < 0.10, we detected 7 TWAS-genes among WW and 1 TWAS-gene among BW. Among WW, CRS showed positive associations with MCM10, FAM64A, CCNB2, and MMP1 GReX and negative associations with VAV3, PCSK6, and GNG11 GReX. Among BW, higher MMP1 GReX predicted lower Proliferation score and ROR-P. TWAS-gene and PAM50 tumor expression associations highlighted potential mechanisms for TWAS-gene to CRS associations. Conclusions: Among BC patients, we find differential germline associations with three CRS by race, underscoring the need for larger, more diverse datasets in molecular studies of BC. Our findings also suggest possible germline trans-regulation of PAM50 tumor expression, with potential implications for interpreting CRS in clinical settings.