Widespread germline genetic heterogeneity of human ribosomal RNA genes

Abstract

Polymorphism drives survival under stress and provides adaptability. Genetic polymorphism of ribosomal RNA (rRNA) genes derives from internal repeat variation of this multicopy gene, and from interindividual variation. A considerable amount of rRNA sequence heterogeneity has been proposed but has been challenging to estimate given the scarcity of accurate reference sequences. We identified four rDNA copies on chromosome 21 (GRCh38) with 99% similarity to recently introduced reference sequence KY962518.1. Pairwise alignment of the rRNA coding sequences of these copies showed differences in sequence and length. We customized a GATK bioinformatics pipeline using the four rDNA loci, spanning a total 145 kb, for variant calling. We employed whole genome sequencing (WGS) data from the 1000 Genomes Project phase 3 and analyzed variants in 2,504 individuals from 26 populations. Using the pipeline, we identified a total of 3,790 variant positions. The variants positioned non-randomly on the rRNA gene. Invariant regions included the promoter, early 5 ETS, 5.8S, ITS1 and certain regions of the 28S rRNA, and large areas of the intragenic spacer. 18S rRNA coding region had very few variants, while a total of 470 variant positions were observed on 28S rRNA. The majority of the 28S rRNA variants located on highly flexible human-expanded rRNA helical folds ES7L and ES27L, suggesting that these represent positions of diversity and are potentially under continuous evolution. These findings provide a genetic view for rRNA heterogeneity and raise the need to functional assess how the 28S rRNA variants affect ribosome functions.