Abstract:Abstract Background Human X and Y chromosomes share an evolutionary origin and, as a consequence, sequence similarity. We investigated whether the sequence homology between the X and Y chromosomes affects the alignment of RNA-Seq reads and estimates of differential expression. We tested the effects of using reference genomes and reference transcriptomes informed by the sex chromosome complement of the sample’s genome on the measurements of RNA-Seq abundance and sex differences in expression. Results The default genome includes the entire human reference genome (GRCh38), including the entire sequence of the X and Y chromosomes. We created two sex chromosome complement informed reference genomes. One sex chromosome complement informed reference genome was used for samples that lacked a Y chromosome; for this reference genome version, we hard-masked the entire Y chromosome. For the other sex chromosome complement informed reference genome, to be used for samples with a Y chromosome, we hard-masked only the pseudoautosomal regions of the Y chromosome, because these regions are duplicated identically in the reference genome on the X chromosome. We analyzed the transcript abundance in the whole blood, brain cortex, breast, liver, and thyroid tissues from 20 genetic female (46, XX) and 20 genetic male (46, XY) samples. Each sample was aligned twice: once to the default reference genome and then independently aligned to a reference genome informed by the sex chromosome complement of the sample, repeated using two different read aligners, HISAT and STAR. We then quantified sex differences in gene expression using featureCounts to get the raw count estimates followed by Limma/Voom for normalization and differential expression. We additionally created sex chromosome complement informed transcriptome references for use in pseudo-alignment using Salmon. Transcript abundance was quantified twice for each sample: once to the default target transcripts and then independently to target transcripts informed by the sex chromosome complement of the sample. Conclusions We show that regardless of the choice of the read aligner, using an alignment protocol informed by the sex chromosome complement of the sample results in higher expression estimates on the pseudoautosomal regions of the X chromosome in both genetic male and genetic female samples, as well as an increased number of unique genes being called as differentially expressed between the sexes. We additionally show that using a pseudo-alignment approach informed on the sex chromosome complement of the sample eliminates Y-linked expression in female XX samples.

Computational protocol for the identification of X-linked genes contributing to X chromosome upregulation from RNA-sequencing datasets

Integrated analysis of Xist upregulation and X-chromosome inactivation with single-cell and single-allele resolution

A self-enhanced transport mechanism through long noncoding RNAs for X chromosome inactivation

Using Population-Scale Transcriptomic and Genomic Data to Map 3′ UTR Alternative Polyadenylation Quantitative Trait Loci

Integrated analysis of Xist upregulation and gene silencing at the onset of random X-chromosome inactivation at high temporal and allelic resolution

Protocol for Quantitative Analysis of RNA 3′-End Processing Induced by Disassociated Subunits Using Chromatin-Associated RNA-seq Data

A Small Proportion of X-linked Genes Contribute to X Chromosome Upregulation in Early Embryos Via BRD4-mediated Transcriptional Activation

Reference genome and transcriptome informed by the sex chromosome complement of the sample increase ability to detect sex differences in gene expression from RNA-Seq data

Inferring genes that escape X-Chromosome inactivation reveals important contribution of variable escape genes to sex-biased diseases

Characterization of X Chromosome Inactivation Using Integrated Analysis of Whole-Exome and mRNA Sequencing

Who's afraid of the X? Incorporating the X and Y chromosomes into the analysis of DNA methylation array data

Identifying down and up-regulated chromosome regions using RNA-Seq data

Elastic dosage compensation by X-chromosome upregulation

Random X-chromosome inactivation dynamics in vivo by single-cell RNA sequencing

RSQ: a Statistical Method for Quantification of Isoform-Specific Structurome Using Transcriptome-Wide Structural Profiling Data

Contribution Value: A Indicator for Measuring the Contribution of Ncrnas to Transcriptome

RNA sequencing shows no dosage compensation of the active X-chromosome

RNA Binding Density on X-chromosome Differing from that on 22 Autosomes in Human

Stable and robust Xi and Y transcriptomes drive cell-type-specific autosomal and Xa responses in vivo and in vitro in four human cell types

Stable and robust Xi and Y transcriptomes drive cell-type-specific autosomal and Xa responses and in four human cell types

Eight quick tips for including chromosome X in genome-wide association studies