Insights into the U1 small nuclear ribonucleoprotein complex superfamily

Examining the capacity of human U1 snRNA variants to facilitate pre-mRNA splicing

The human U1 snRNA is encoded by a multigene family consisting of transcribed variants and defective pseudogenes. Many variant U1 (vU1) snRNAs have been demonstrated to not only be transcribed but also processed by the addition of a trimethylated guanosine cap, packaged into snRNPs, and assembled into spliceosomes; however, their capacity to facilitate pre-mRNA splicing has, so far, not been tested. A recent systematic analysis of the human snRNA genes identified 178 U1 snRNA genes that are present in the genome as either tandem arrays or single genes on multiple chromosomes. Of these, 15 were found to be expressed in human tissues and cell lines, although at significantly low levels from their endogenous loci, <0.001% of the canonical U1 snRNA. In this study, we found that placing the variants in the context of the regulatory elements of the RNU1-1 gene improves the expression of many variants to levels comparable to the canonical U1 snRNA. Application of a previously established HeLa cell-based minigene reporter assay to examine the capacity of the vU1 snRNAs to support pre-mRNA splicing revealed that even though the exogenously expressed variant snRNAs were enriched in the nucleus, only a few had a measurable effect on splicing.

Development of Engineered-U1 snRNA Therapies: Current Status

Splicing of pre-mRNA is a crucial regulatory stage in the pathway of gene expression. The majority of human genes that encode proteins undergo alternative pre-mRNA splicing and mutations that affect splicing are more prevalent than previously thought. Targeting aberrant RNA(s) may thus provide an opportunity to correct faulty splicing and potentially treat numerous genetic disorders. To that purpose, the use of engineered U1 snRNA (either modified U1 snRNAs or exon-specific U1s-ExSpeU1s) has been applied as a potentially therapeutic strategy to correct splicing mutations, particularly those affecting the 5' splice-site (5'ss). Here we review and summarize a vast panoply of studies that used either modified U1 snRNAs or ExSpeU1s to mediate gene therapeutic correction of splicing defects underlying a considerable number of genetic diseases. We also focus on the pre-clinical validation of these therapeutic approaches both in vitro and in vivo, and summarize the main obstacles that need to be overcome to allow for their successful translation to clinic practice in the future.

Splicing factor SNRPA associated with microvascular invasion promotes hepatocellular carcinoma metastasis through activating NOTCH1/Snail pathway and is mediated by circSEC62/miR-625-5p axis

Microvascular invasion (MVI) is a crucial risk factor related to the metastasis of hepatocellular carcinoma (HCC), but the underlying mechanisms remain to be revealed. Characterizing the inherent mechanisms of MVI may aid in the development of effective treatment strategies to improve the prognosis of HCC patients with metastasis. Through the Gene Expression Omnibus (GEO) database, we identified that small nuclear ribonucleoprotein polypeptide A (SNRPA) was related to MVI in HCC. SNRPA was overexpressed in MVI-HCC and correlated with poor patient survival. Mechanistically, SNRPA promoted the epithelial-mesenchymal transition (EMT)-like process for HCC cells to accelerate metastasis by activating the NOTCH1/Snail pathway in vitro and in vivo. Importantly, circSEC62 upregulated SNRPA expression in HCC cells via miR-625-5p sponging. Taking these results together, our study identified a novel regulatory mechanism among SNRPA, miR-625-5p, circSEC62 and the NOTCH1/Snail pathway in HCC, which promoted metastasis of HCC and may provide effective suggestions for improving the prognosis of HCC patients with metastasis.

U1 snRNP proteins promote proximal alternative polyadenylation sites by directly interacting with 3' end processing core factors

In eukaryotic cells, both alternative splicing and alternative polyadenylation (APA) play essential roles in the gene regulation network. U1 small ribonucleoprotein particle (U1 snRNP) is a major component of spliceosome, and U1 snRNP complex can suppress proximal APA sites through crosstalking with 3' end processing factors. However, here we show that both knockdown and overexpression of SNRPA, SNRPC, SNRNP70, and SNRPD2, the U1 snRNP proteins, promote the usage of proximal APA sites at the transcriptome level. SNRNP70 can drive the phase transition of PABPN1 from droplet to aggregate, which may reduce the repressive effects of PABPN1 on the proximal APA sites. Additionally, SNRNP70 can also promote the proximal APA sites by recruiting CPSF6, suggesting that the function of CPSF6 on APA is related with other RNA-binding proteins and cell context-dependent. Consequently, these results reveal that, on the contrary to U1 snRNP complex, the free proteins of U1 snRNP complex can promote proximal APA sites through the interaction with 3' end processing machinery.

Engineered U1 snRNAs to modulate alternatively spliced exons

Alternative splicing accounts for a considerable portion of transcriptomic diversity, as most protein-coding genes are spliced into multiple mRNA isoforms. However, errors in splicing patterns can give rise to mis-splicing with pathological consequences, such as the congenital diseases familial dysautonomia, Duchenne muscular dystrophy, and spinal muscular atrophy. Small nuclear RNA (snRNA) components of the U snRNP family have been proposed as a therapeutic modality for the treatment of mis-splicing. U1 snRNAs offer great promise, with prior studies demonstrating in vivo efficacy, suggesting additional preclinical development is merited. Improvements in enabling technologies, including screening methodologies, gene delivery vectors, and relevant considerations from gene editing approaches justify further advancement of U1 snRNA as a therapeutic and research tool. With the goal of providing a user-friendly protocol, we compile and demonstrate a methodological toolkit for sequence-specific targeted perturbation of alternatively spliced pre-mRNA with engineered U1 snRNAs. We observe robust modulation of endogenous pre-mRNA transcripts targeted in two contrasting splicing contexts, SMN2 exon 7 and FAS exon 6, exhibiting the utility and applicability of engineered U1 snRNA to both inclusion and exclusion of targeted exons. We anticipate that these demonstrations will contribute to the usability of U1 snRNA in investigating splicing modulation in eukaryotic cells, increasing accessibility to the broader research community.

Exon Skipping Through Chimeric Antisense U1 snRNAs to Correct Retinitis Pigmentosa GTPase-Regulator (RPGR) Splice Defect

Inherited retinal dystrophies are caused by mutations in more than 250 genes, each of them carrying several types of mutations that can lead to different clinical phenotypes. Mutations in Retinitis Pigmentosa GTPase-Regulator (RPGR) cause X-linked Retinitis pigmentosa (RP). A nucleotide substitution in intron 9 of RPGR causes the increase of an alternatively spliced isoform of the mature mRNA, bearing exon 9a (E9a). This introduces a stop codon, leading to truncation of the protein. Aiming at restoring impaired gene expression, we developed an antisense RNA-based therapeutic approach for the skipping of RPGR E9a. We designed a set of specific U1 antisense snRNAs (U1_asRNAs) and tested their efficacy in vitro, upon transient cotransfection with RPGR minigene reporter systems in HEK-293T, 661W, and PC-12 cell lines.

The upstream 5' splice site remains associated to the transcription machinery during intron synthesis

In the earliest step of spliceosome assembly, the two splice sites flanking an intron are brought into proximity by U1 snRNP and U2AF along with other proteins. The mechanism that facilitates this intron looping is poorly understood. Using a CRISPR interference-based approach to halt RNA polymerase II transcription in the middle of introns in human cells, we discovered that the nascent 5′ splice site base pairs with a U1 snRNA that is tethered to RNA polymerase II during intron synthesis. This association functionally corresponds with splicing outcome, involves bona fide 5′ splice sites and cryptic intronic sites, and occurs transcriptome-wide. Overall, our findings reveal that the upstream 5′ splice sites remain attached to the transcriptional machinery during intron synthesis and are thus brought into proximity of the 3′ splice sites; potentially mediating the rapid splicing of long introns.

We know that most splicing reactions take place co-transcriptionally, but how the transcription machinery facilitate splicing of introns is unknown. Here the authors show that the 5′ splice site remains associated with the transcription machinery during intron synthesis through U1 snRNP, providing a basis for the rapid splicing reaction of introns.

Presenilin 1 mutation likely contributes to U1 small nuclear RNA dysregulation and Alzheimer's disease-like symptoms

Previous studies showed that U1 small nuclear RNA (snRNA) was selectively enriched in the brain of individuals with familial Alzheimer's disease (AD), resulting in widespread changes in RNA splicing. Our study further reported that presenilin-1 (PSEN1) induced an increase in U1 snRNA expression, accompanied by changed amyloid precursor protein expression, β-amyloid level, and cell death in SH-SY5Y cells. However, the effect of U1 snRNA overexpression on learning and memory is still unclear. In the present study, we found that neuronal U1 snRNA overexpression could generate U1 snRNA aggregates in the nuclear, accompanied by the widespread alteration of RNA splicing, resulting in the impairments of synaptic plasticity and spatial memory. In addition, more U1 snRNAs is bound to the intron binding sites accompanied by an increased intracellular U1 snRNA level. This suggests that U1 snRNA overexpression regulates RNA splicing and gene expression in neurons by manipulating the recruitment of the U1 snRNA to the nascent transcripts. Using in situ hybridization staining of human central nervous system-type neurons, we identified nuclear aggregates of U1 snRNA in neurons by upregulating the U1 snRNA level. Quantitative polymerase chain reaction analysis showed U1 snRNA accumulation in the insoluble fraction of neurons with PSEN1 mutation neurons rather than other types of U snRNAs. These results show an independent function of U1 snRNA in regulating RNA splicing, suggesting that aberrant RNA processing may mediate neurodegeneration induced by PSEN1 mutation.

U1 snRNP telescripting: molecular mechanisms and beyond

U1 snRNP is one of the most abundant ribonucleoprotein (RNP) complexes in eukaryotic cells and is estimated to be approximately 1 million copies per cell. Apart from its canonical role in mRNA splicing, this complex has emerged as a key regulator of eukaryotic mRNA length via inhibition of mRNA 3'-end processing at numerous intronic polyadenylation sites, in a process that is also termed 'U1 snRNP telescripting'. Several reviews have extensively described the concept of U1 telescripting and subsequently highlighted its potential impacts in mRNA metabolism. Here, we review what is currently known regarding the underlying mechanisms of this important phenomenon and discuss open questions and future challenges.

Noncoding RNAs Set the Stage for RNA Polymerase II Transcription

Effective synthesis of mammalian messenger (m)RNAs depends on many factors that together direct RNA polymerase II (pol II) through the different stages of the transcription cycle and ensure efficient cotranscriptional processing of mRNAs. In addition to the many proteins involved in transcription initiation, elongation, and termination, several noncoding (nc)RNAs also function as global transcriptional regulators. Understanding the mode of action of these non-protein regulators has been an intense area of research in recent years. Here, we describe how these ncRNAs influence key regulatory steps of the transcription process, to affect large numbers of genes. Through direct association with pol II or by modulating the activity of transcription or RNA processing factors, these regulatory RNAs perform critical roles in gene expression.

5' and 3' splicing signals evolution in vertebrates: Analysis in a conserved gene family

The mitochondrial solute carrier genes (SLC25) are highly conserved during vertebrate evolution. In most SLC25 genes of zebrafish, chicken, mouse, and human, the introns are located at exactly superimposable positions. In these topographically corresponding introns we studied the composition of the initial and terminal hexanucleotides (5'ss and 3'ss) which are instrumental in splicing signaling, focusing on the evolutionary conservation/mutation dynamics of these genetically related sequences. At each position, the per cent conservation of zebrafish individual nucleotides in chicken, mouse and human is proportional to their percent frequency in zebrafish; furthermore, nucleotide mutations are biased in favor of the more represented nucleotides, thus compensating for those highly represented zebrafish nucleotides which have not been conserved. As a result of these evolutionary dynamics, the general nucleotide composition at each position has remained relatively conserved throughout vertebrates. At 5'ss, following the canonical GT, A and G are largely prevailing at position +3, A at +4 and G at +5 (GT[A/G]AGx). At 3'ss, T and C are largely prevailing at positions -6, -5 and -3, preceding the canonical intron terminal AG ([C/T] [C/T]x[C/T]AG). However, the actual composition of the tetranucleotides at 5' and 3' often does not conform to the above scheme. At 5'ss the more canonical sequence is completely expressed in 63% of cases and partially (2 or 1 matches) in 37 % of cases. At 3'ss the more canonical sequence is completely expressed in 71 % of cases and partially (2 or 1 matches) in 29 % of cases. The nucleotide conservation loss (nucleotide mutation) is higher in the evolution from fish to the last common ancestor of birds and mammals (58 %), then diminishes in the successive evolution steps up to the mammalian common ancestor (10 %), and becomes still lower at the divergence of rodents and primates (5 %).

RNA Therapeutics: How Far Have We Gone?

In recent years, the RNA molecule became one of the most promising targets for therapeutic intervention. Currently, a large number of RNA-based therapeutics are being investigated both at the basic research level and in late-stage clinical trials. Some of them are even already approved for treatment. RNA-based approaches can act at pre-mRNA level (by splicing modulation/correction using antisense oligonucleotides or U1snRNA vectors), at mRNA level (inhibiting gene expression by siRNAs and antisense oligonucleotides) or at DNA level (by editing mutated sequences through the use of CRISPR/Cas). Other RNA approaches include the delivery of in vitro transcribed (IVT) mRNA or the use of oligonucleotides aptamers. Here we review these approaches and their translation into clinics trying to give a brief overview also on the difficulties to its application as well as the research that is being done to overcome them.

A novel role of U1 snRNP: Splice site selection from a distance

Removal of introns by pre-mRNA splicing is fundamental to gene function in eukaryotes. However, understanding the mechanism by which exon-intron boundaries are defined remains a challenging endeavor. Published reports support that the recruitment of U1 snRNP at the 5'ss marked by GU dinucleotides defines the 5'ss as well as facilitates 3'ss recognition through cross-exon interactions. However, exceptions to this rule exist as U1 snRNP recruited away from the 5'ss retains the capability to define the splice site, where the cleavage takes place. Independent reports employing exon 7 of Survival Motor Neuron (SMN) genes suggest a long-distance effect of U1 snRNP on splice site selection upon U1 snRNP recruitment at target sequences with or without GU dinucleotides. These findings underscore that sequences distinct from the 5'ss may also impact exon definition if U1 snRNP is recruited to them through partial complementarity with the U1 snRNA. In this review we discuss the expanded role of U1 snRNP in splice-site selection due to U1 ability to be recruited at more sites than predicted solely based on GU dinucleotides.

Regulation of expression of human RNA polymerase II-transcribed snRNA genes

In addition to protein-coding genes, RNA polymerase II (pol II) transcribes numerous genes for non-coding RNAs, including the small-nuclear (sn)RNA genes. snRNAs are an important class of non-coding RNAs, several of which are involved in pre-mRNA splicing. The molecular mechanisms underlying expression of human pol II-transcribed snRNA genes are less well characterized than for protein-coding genes and there are important differences in expression of these two gene types. Here, we review the DNA features and proteins required for efficient transcription of snRNA genes and co-transcriptional 3′ end formation of the transcripts.

QM/MM Calculations on Protein-RNA Complexes: Understanding Limitations of Classical MD Simulations and Search for Reliable Cost-Effective QM Methods

Although atomistic explicit-solvent Molecular Dynamics (MD) is a popular tool to study protein-RNA recognition, satisfactory MD description of protein-RNA complexes is not always achieved. Unfortunately, it is often difficult to separate MD simulation instabilities primarily caused by the simple point-charge molecular mechanics (MM) force fields from problems related to the notorious uncertainties in the starting structures. Herein, we report a series of large-scale QM/MM calculations on the U1A protein-RNA complex. This experimentally well-characterized system has an intricate protein-RNA interface, which is very unstable in MD simulations. The QM/MM calculations identify several H-bonds poorly described by the MM method and thus indicate the sources of instabilities of the U1A interface in MD simulations. The results suggest that advanced QM/MM computations could be used to indirectly rationalize problems seen in MM-based MD simulations of protein-RNA complexes. As the most accurate QM method, we employ the computationally demanding meta-GGA density functional TPSS-D3(BJ)/def2-TZVP level of theory. Because considerably faster methods would be needed to extend sampling and to study even larger protein-RNA interfaces, a set of low-cost QM/MM methods is compared to the TPSS-D3(BJ)/def2-TZVP data. The PBEh-3c and B97-3c density functional composite methods appear to be suitable for protein-RNA interfaces. In contrast, HF-3c and the tight-binding Hamiltonians DFTB3-D3 and GFN-xTB perform unsatisfactorily and do not provide any advantage over the MM description. These conclusions are supported also by similar analysis of a simple HutP protein-RNA interface, which is well-described by MD with the exception of just one H-bond. Some other methodological aspects of QM/MM calculations on protein-RNA interfaces are discussed.

SNRPA enhances tumour cell growth in gastric cancer through modulating NGF expression

OBJECTIVES

SNRPA is a protein component of U1 small nuclear ribonucleoprotein (U1 snRNP) complex, which takes part in the splicing of pre-mRNAs. Its expression and function in tumour remain unknown. Herein, we elucidated the functional contribution of SNRPA to the progression of gastric cancer (GC).

MATERIALS AND METHODS

SNRPA expression was investigated in a GC tissue microarray by immunohistochemical staining. Cell proliferation was evaluated by CCK-8, colony formation and EdU incorporation assays. A mouse xenograft model was used to detect the tumourigenicity. Gene expression profiling was performed and then the potential target genes were verified by quantitative real-time PCR and western blot analyses. The functional relevance between SNRPA and its target gene was examined by cell growth assays.

RESULTS

SNRPA expression was higher in tumour tissues than in matched normal gastric mucosa tissues, and it was positively correlated with the tumour size and progression. High SNRPA expression indicated poor prognosis of GC patients. Silencing SNRPA in GC cells markedly inhibited cell proliferation in vitro and tumour growth in a xenograft model, while overexpressing SNRPA exhibited opposite results. Moreover, we identified NGF (Nerve growth factor) as a downstream effector of SNRPA and further proved that NGF was crucial for SNRPA-mediated GC cell growth.

CONCLUSIONS

These findings suggested that SNRPA may contribute to GC progression via NGF and could be a prognostic biomarker for GC.

Gene expression profiling reveals U1 snRNA regulates cancer gene expression

U1 small nuclear RNA (U1 snRNA), as one of the most abundant ncRNAs in human cells, plays an important role in splicing of pre-mRNAs. Compared to previous studies which have focused on the primary function of U1 snRNA and the neurodegenerative diseases caused by abnormalities of U1 snRNA, this study is to investigate how U1 snRNA over-expression affects the expression of mammal genes on a genome-wide scale. By comparing the gene expression profiles of U1 snRNA over-expressed cells with those of their controls using microarray experiments, 916 genes or loci were identified significantly Differentially Expressed (DE). These 595 up-regulated DE genes and 321 down-regulated DE genes were analyzed using annotations from GO categories and pathways from the KEGG database. As a result, three of 12 enriched pathways were well-known cancer pathways, while the other nine pathways were associated to cancers in previous studies. The further analysis of 73 genes involved in 12 pathways suggested that U1 snRNA could regulate cancer gene expression. The microarray data under the GEO Series accession number GSE84304 is available in the NCBI GEO database.

Overexpression of U1 snRNA induces decrease of U1 spliceosome function associated with Alzheimer's disease

We recently reported that presenilin-1 (PS1) induced an increase of U1 snRNA expression accompanied with the change of amyloid precursor protein expression, β-amyloid level and cell death. In the present study, our data showed that both overexpression and knockdown of U1 snRNA could cause the loss in the function of U1 snRNA and resulted in PCPA as well as the same downstream phenomena including the expression changes of genes specific to AD, tau hyperphosphorylation on the site of Thr212, the decrease of acetylated α-tubulin, the reduction of cell viability and upregulation of RIPK1, RIPK3 and caspase8. These findings not only helped researchers better understand the functions of U1 snRNA, but also paved the way to reveal the mechanisms of AD from a different point of view and may find a new therapeutic target for the disease.

How introns enhance gene expression

In many eukaryotes, including mammals, plants, yeast, and insects, introns can increase gene expression without functioning as a binding site for transcription factors. This phenomenon was termed 'intron-mediated enhancement'. Introns can increase transcript levels by affecting the rate of transcription, nuclear export, and transcript stability. Moreover, introns can also increase the efficiency of mRNA translation. This review discusses the current knowledge about these processes. The role of splicing in IME and the significance of intron position relative to the sites of transcription and translation initiation are elaborated. Particular emphasis is placed on the question why different introns, present at the same location of the same genes and spliced at a similar high efficiency, can have very different impacts on expression - from almost no effect to considerable stimulation. This situation can be at least partly accounted for by the identification of splicing-unrelated intronic elements with a special ability to enhance mRNA accumulation or translational efficiency. The many factors that could lead to the large variation observed between the impact of introns in different genes and experimental systems are highlighted. It is suggested that there is no sole, definite answer to the question "how do introns enhance gene expression". Rather, each intron-gene combination might undergo its own unique mixture of processes that lead to the perceptible outcome.

Variant U1 snRNAs are implicated in human pluripotent stem cell maintenance and neuromuscular disease

The U1 small nuclear (sn)RNA (U1) is a multifunctional ncRNA, known for its pivotal role in pre-mRNA splicing and regulation of RNA 3' end processing events. We recently demonstrated that a new class of human U1-like snRNAs, the variant (v)U1 snRNAs (vU1s), also participate in pre-mRNA processing events. In this study, we show that several human vU1 genes are specifically upregulated in stem cells and participate in the regulation of cell fate decisions. Significantly, ectopic expression of vU1 genes in human skin fibroblasts leads to increases in levels of key pluripotent stem cell mRNA markers, including NANOG and SOX2. These results reveal an important role for vU1s in the control of key regulatory networks orchestrating the transitions between stem cell maintenance and differentiation. Moreover, vU1 expression varies inversely with U1 expression during differentiation and cell re-programming and this pattern of expression is specifically de-regulated in iPSC-derived motor neurons from Spinal Muscular Atrophy (SMA) type 1 patient's. Accordingly, we suggest that an imbalance in the vU1/U1 ratio, rather than an overall reduction in Uridyl-rich (U)-snRNAs, may contribute to the specific neuromuscular disease phenotype associated with SMA.

Double-target Antisense U1snRNAs Correct Mis-splicing Due to c.639+861C>T and c.639+919G>A GLA Deep Intronic Mutations

Fabry disease is a rare X-linked lysosomal storage disorder caused by deficiency of the α-galactosidase A (α-Gal A) enzyme, which is encoded by the GLA gene. GLA transcription in humans produces a major mRNA encoding α-Gal A and a minor mRNA of unknown function, which retains a 57-nucleotide-long cryptic exon between exons 4 and 5, bearing a premature termination codon. NM_000169.2:c.639+861C>T and NM_000169.2:c.639+919G>A GLA deep intronic mutations have been described to cause Fabry disease by inducing overexpression of the alternatively spliced mRNA, along with a dramatic decrease in the major one. Here, we built a wild-type GLA minigene and two minigenes that carry mutations c.639+861C>T and c.639+919G>A. Once transfected into cells, the minigenes recapitulate the molecular patterns observed in patients, at the mRNA, protein, and enzymatic level. We constructed a set of specific double-target U1asRNAs to correct c.639+861C>T and c.639+919G>A GLA mutations. Efficacy of U1asRNAs in inducing the skipping of the cryptic exon was evaluated upon their transient co-transfection with the minigenes in COS-1 cells, by real-time polymerase chain reaction (PCR), western blot analysis, and α-Gal A enzyme assay. We identified a set of U1asRNAs that efficiently restored α-Gal A enzyme activity and the correct splicing pathways in reporter minigenes. We also identified a unique U1asRNA correcting both mutations as efficently as the mutation-specific U1asRNAs. Our study proves that an exon skipping-based approach recovering α-Gal A activity in the c.639+861C>T and c.639+919G>A GLA mutations is active.

Large Scale Gene Expression Meta-Analysis Reveals Tissue-Specific, Sex-Biased Gene Expression in Humans

The severity and prevalence of many diseases are known to differ between the sexes. Organ specific sex-biased gene expression may underpin these and other sexually dimorphic traits. To further our understanding of sex differences in transcriptional regulation, we performed meta-analyses of sex biased gene expression in multiple human tissues. We analyzed 22 publicly available human gene expression microarray data sets including over 2500 samples from 15 different tissues and 9 different organs. Briefly, by using an inverse-variance method we determined the effect size difference of gene expression between males and females. We found the greatest sex differences in gene expression in the brain, specifically in the anterior cingulate cortex, (1818 genes), followed by the heart (375 genes), kidney (224 genes), colon (218 genes), and thyroid (163 genes). More interestingly, we found different parts of the brain with varying numbers and identity of sex-biased genes, indicating that specific cortical regions may influence sexually dimorphic traits. The majority of sex-biased genes in other tissues such as the bladder, liver, lungs, and pancreas were on the sex chromosomes or involved in sex hormone production. On average in each tissue, 32% of autosomal genes that were expressed in a sex-biased fashion contained androgen or estrogen hormone response elements. Interestingly, across all tissues, we found approximately two-thirds of autosomal genes that were sex-biased were not under direct influence of sex hormones. To our knowledge this is the largest analysis of sex-biased gene expression in human tissues to date. We identified many sex-biased genes that were not under the direct influence of sex chromosome genes or sex hormones. These may provide targets for future development of sex-specific treatments for diseases.