RESUMEN
Microproteins encoded by small open reading frames (smORFs) comprise the "dark matter" of proteomes. Although functional microproteins were identified in diverse organisms from all three domains of life, bacterial smORFs remain poorly characterized. In this comprehensive study of intergenic smORFs (ismORFs, 15-70 codons) in 5,668 bacterial genomes of the family Enterobacteriaceae, we identified 67,297 clusters of ismORFs subject to purifying selection. The ismORFs mainly code for hydrophobic, potentially transmembrane, unstructured, or minimally structured microproteins. Using AlphaFold Multimer, we predicted interactions of some of the predicted microproteins encoded by transcribed ismORFs with proteins encoded by neighboring genes, revealing the potential of microproteins to regulate the activity of various proteins, particularly, under stress. We compiled a catalog of predicted microprotein families with different levels of evidence from synteny analysis, structure prediction, and transcription and translation data. This study offers a resource for investigation of biological functions of microproteins.
RESUMEN
Ribosomal DNA (rDNA) repeat units are organized into tandem clusters in eukaryotic cells. In mice, these clusters are located on at least eight chromosomes and show extensive variation in the number of repeats between mouse genomes. To analyze intra- and inter-genomic variation of mouse rDNA repeats, we selectively isolated 25 individual rDNA units using Transformation-Associated Recombination (TAR) cloning. Long-read sequencing and subsequent comparative sequence analysis revealed that each full-length unit comprises an intergenic spacer (IGS) and a â¼13.4 kb long transcribed region encoding the three rRNAs, but with substantial variability in rDNA unit size, ranging from â¼35 to â¼46 kb. Within the transcribed regions of rDNA units, we found 209 variants, 70 of which are in external transcribed spacers (ETSs); but the rDNA size differences are driven primarily by IGS size heterogeneity, due to indels containing repetitive elements and some functional signals such as enhancers. Further evolutionary analysis categorized rDNA units into distinct clusters with characteristic IGS lengths; numbers of enhancers; and presence/absence of two common SNPs in promoter regions, one of which is located within promoter (p)RNA and may influence pRNA folding stability. These characteristic features of IGSs also correlated significantly with 5'ETS variant patterns described previously and associated with differential expression of rDNA units. Our results suggest that variant rDNA units are differentially regulated and open a route to investigate the role of rDNA variation on nucleolar formation and possible associations with pathology.
RESUMEN
Objective: To investigate the relationship between vaccination rates and excess mortality during distinct waves of SARS-CoV-2 variant-specific infections, while considering a state's GDP per capita. Methods: We ranked U.S. states by vaccination rates and GDP and employed the CDC's excess mortality model for regression and odds ratio analysis. Results: Regression analysis reveals that both vaccination and GDP are significant factors related to mortality when considering the entire U.S. population. Notably, in wealthier states (with GDP above $65,000), excess mortality is primarily driven by slow vaccination rates, while in less affluent states, low GDP plays a major role. Odds ratio analysis demonstrates an almost twofold increase in mortality linked to the Delta and Omicron BA.1 virus variants in states with the slowest vaccination rates compared to those with the fastest (OR 1.8, 95% CI 1.7-1.9, p < 0.01). However, this gap disappeared in the post-Omicron BA.1 period. Conclusion: The interplay between slow vaccination and low GDP per capita drives high mortality.
RESUMEN
Aim: To evaluate the effect of vaccination/booster administration dynamics on the reduction of excess mortality during COVID-19 infection waves in European countries. Methods: We selected twenty-nine countries from the OurWorldInData project database according to their population size of more than one million and the availability of information on dominant SARS-CoV-2 variants during COVID-19 infection waves. After selection, we categorized countries according to their "faster" or "slower" vaccination rates. The first category included countries that reached 60% of vaccinated residents by October 2021 and 70% by January 2022. The second or "slower" category included all other countries. In the first or "faster" category, two groups, "boosters faster'' and "boosters slower" were created. Pearson correlation analysis, linear regression, and chi-square test for categorical data were used to identify the association between vaccination rate and excess mortality. We chose time intervals corresponding to the dominance of viral variants: Wuhan, Alpha, Delta, and Omicron BA.1/2. Results and discussion: The "faster" countries, as opposed to the "slower" ones, did better in protecting their residents from mortality during all periods of the SARS-CoV-2 pandemic and even before vaccination. Perhaps higher GDP per capita contributed to their better performance throughout the pandemic. During mass vaccination, when the Delta variant prevailed, the contrast in mortality rates between the "faster" and "slower" categories was strongest. The average excess mortality in the "slower" countries was nearly 5 times higher than in the "faster" countries, and the odds ratio (OR) was 4.9 (95% CI 4.4 to 5.4). Slower booster rates were associated with significantly higher mortality during periods dominated by Omicron BA.1 and BA.2, with an OR of 2.6 (CI 95%. 2.1 to 3.3). Among the European countries we analyzed, Denmark, Norway, and Ireland did best, with a pandemic mortality rate of 0.1% of the population or less. By comparison, Bulgaria, Serbia, and Russia had a much higher mortality rate of up to 1% of the population. Conclusion: Thus, slow vaccination and booster administration was a major factor contributing to an order of magnitude higher excess mortality in "slower" European countries compared to more rapidly immunized countries.
Asunto(s)
COVID-19 , Humanos , SARS-CoV-2 , Pandemias , VacunaciónRESUMEN
Although they are organelles without a limiting membrane, nucleoli have an exclusive structure, built upon the rDNA-rich acrocentric short arms of five human chromosomes (nucleolar organizer regions or NORs). This has raised the question: what are the structural features of a chromosome required for its inclusion in a nucleolus? Previous work has suggested that sequences adjacent to the tandemly repeated rDNA repeat units (DJ, distal junction sequence) may be involved, and we have extended such studies by addressing several issues related to the requirements for the association of NORs with nucleoli. We exploited both a set of somatic cell hybrids containing individual human acrocentric chromosomes and a set of Human Artificial Chromosomes (HACs) carrying different parts of a NOR, including an rDNA unit or DJ or PJ (proximal junction) sequence. Association of NORs with nucleoli was increased when constituent rDNA was transcribed and may be also affected by the status of heterochromatin blocks formed next to the rDNA arrays. Furthermore, our data suggest that a relatively small size DJ region, highly conserved in evolution, is also involved, along with the rDNA repeats, in the localization of p-arms of acrocentric chromosomes in nucleoli. Thus, we infer a cooperative action of rDNA sequence-stimulated by its activity-and sequences distal to rDNA contributing to incorporation into nucleoli. Analysis of NOR sequences also identified LncRNA_038958 in the DJ, a candidate transcript with the region of the suggested promoter that is located close to the DJ/rDNA boundary and contains CTCF binding sites. This LncRNA may affect RNA Polymerase I and/or nucleolar activity. Our findings provide the basis for future studies to determine which RNAs and proteins interact critically with NOR sequences to organize the higher-order structure of nucleoli and their function in normal cells and pathological states.
Asunto(s)
Región Organizadora del Nucléolo , ARN Largo no Codificante , Humanos , Región Organizadora del Nucléolo/genética , Región Organizadora del Nucléolo/metabolismo , ADN Ribosómico/genética , ARN Largo no Codificante/metabolismo , Nucléolo Celular/genética , Nucléolo Celular/metabolismo , Cromosomas Humanos/metabolismoRESUMEN
Aim To evaluate the effect of vaccination/booster administration dynamics on the reduction of excess mortality during COVID-19 infection waves in European countries. Methods We selected twenty-nine countries from the OurWorldInData project database according to their population size of more than one million and the availability of information on dominant SARS-CoV-2 variants during COVID-19 infection waves. After selection, we categorized countries according to their â³fasterâ³ or â³slowerâ³ vaccination rates. The first category included countries that reached 60% of vaccinated residents by October 2021 and 70% by January 2022. The second or â³slowerâ³ category included all other countries. In the first or â³fasterâ³ category, two groups, â³boosters faster'' and â³boosters slowerâ³ were created. Pearson correlation analysis, linear regression, and chi-square test for categorical data were used to identify the association between vaccination rate and excess mortality. We chose time intervals corresponding to the dominance of viral variants: Wuhan, Alpha, Delta, and Omicron BA.1/2. Results The â³fasterâ³ countries, as opposed to the â³slowerâ³ ones, did better in protecting their residents from mortality during all periods of the SARS-CoV-2 pandemic and even before vaccination. Perhaps higher GDP per capita contributed to their better performance throughout the pandemic. During mass vaccination, when the Delta variant prevailed, the contrast in mortality rates between the â³fasterâ³ and â³slowerâ³ categories was strongest. The average excess mortality in the â³slowerâ³ countries was nearly 5 times higher than in the â³fasterâ³ countries, and the odds ratio (OR) was 4.9 (95% CI 4.4 to 5.4). Slower booster rates were associated with significantly higher mortality during periods dominated by Omicron BA.1 and BA.2, with an OR of 2.6 (CI 95%. 2.1 to 3.3). Among the European countries we analyzed, Denmark, Norway, and Ireland did best, with a pandemic mortality rate of 0.1% of the population or less. By comparison, Bulgaria, Serbia, and Russia had a much higher mortality rate of up to 1% of the population. Thus, slow vaccination and booster administration was a major factor contributing to an order of magnitude higher excess mortality in â³slowerâ³ European countries compared to more rapidly immunized countries.
RESUMEN
Increasing numbers of small, regulatory RNAs (sRNAs) corresponding to 3' untranslated regions (UTR) are being discovered in bacteria. One such sRNA, denoted GlnZ, corresponds to the 3' UTR of the Escherichia coli glnA mRNA encoding glutamine synthetase. Several forms of GlnZ, processed from the glnA mRNA, are detected in cells growing with limiting ammonium. GlnZ levels are regulated transcriptionally by the NtrC transcription factor and post-transcriptionally by RNase III. Consistent with the expression, E. coli cells lacking glnZ show delayed outgrowth from nitrogen starvation compared to wild type cells. Transcriptome-wide RNA-RNA interactome datasets indicated that GlnZ binds to multiple target RNAs. Immunoblots and assays of fusions confirmed GlnZ-mediated repression of glnP and sucA, encoding proteins that contribute to glutamine transport and the citric acid cycle, respectively. Although the overall sequences of GlnZ from E. coli K-12, Enterohemorrhagic E. coli and Salmonella enterica have significant differences due to various sequence insertions, all forms of the sRNA were able to regulate the two targets characterized. Together our data show that GlnZ impacts growth of E. coli under low nitrogen conditions by modulating genes that affect carbon and nitrogen flux.
Asunto(s)
Compuestos de Amonio , Microbioma Gastrointestinal , ARN Pequeño no Traducido , Regiones no Traducidas 3' , Carbono/metabolismo , Escherichia coli/genética , Escherichia coli/metabolismo , Regulación Bacteriana de la Expresión Génica , Glutamato-Amoníaco Ligasa/genética , Glutamina/genética , Glutamina/metabolismo , Nitrógeno/metabolismo , ARN Mensajero/genética , ARN Mensajero/metabolismo , ARN Pequeño no Traducido/genética , ARN Pequeño no Traducido/metabolismo , Ribonucleasa III/metabolismo , Factores de Transcripción/metabolismoRESUMEN
The choice of guide RNA (gRNA) for CRISPR-based gene targeting is an essential step in gene editing applications, but the prediction of gRNA specificity remains challenging. Lack of transparency and focus on point estimates of efficiency disregarding the information on possible error sources in the model limit the power of existing Deep Learning-based methods. To overcome these problems, we present a new approach, a hybrid of Capsule Networks and Gaussian Processes. Our method predicts the cleavage efficiency of a gRNA with a corresponding confidence interval, which allows the user to incorporate information regarding possible model errors into the experimental design. We provide the first utilization of uncertainty estimation in computational gRNA design, which is a critical step toward accurate decision-making for future CRISPR applications. The proposed solution demonstrates acceptable confidence intervals for most test sets and shows regression quality similar to existing models. We introduce a set of criteria for gRNA selection based on off-target cleavage efficiency and its variance and present a collection of pre-computed gRNAs for human chromosome 22. Using Neural Network Interpretation methods, we show that our model rediscovers an established biological factor underlying cleavage efficiency, the importance of the seed region in gRNA.
Asunto(s)
Sistemas CRISPR-Cas , Aprendizaje Profundo , Edición Génica , Marcación de Gen , ARN Guía de Kinetoplastida/genética , Algoritmos , Edición Génica/métodos , Marcación de Gen/métodos , Genómica/métodos , Humanos , Redes Neurales de la Computación , Reproducibilidad de los ResultadosRESUMEN
Pervasive transcription of eukaryotic genomes results in expression of long non-coding RNAs (lncRNAs) most of which are poorly conserved in evolution and appear to be non-functional. However, some lncRNAs have been shown to perform specific functions, in particular, transcription regulation. Thousands of small open reading frames (smORFs, <100 codons) located on lncRNAs potentially might be translated into peptides or microproteins. We report a comprehensive analysis of the conservation and evolutionary trajectories of lncRNAs-smORFs from the moss Physcomitrium patens across transcriptomes of 479 plant species. Although thousands of smORFs are subject to substantial purifying selection, the majority of the smORFs appear to be evolutionary young and could represent a major pool for functional innovation. Using nanopore RNA sequencing, we show that, on average, the transcriptional level of conserved smORFs is higher than that of non-conserved smORFs. Proteomic analysis confirmed translation of 82 novel species-specific smORFs. Numerous conserved smORFs containing low complexity regions (LCRs) or transmembrane domains were identified, the biological functions of a selected LCR-smORF were demonstrated experimentally. Thus, microproteins encoded by smORFs are a major, functionally diverse component of the plant proteome.
Asunto(s)
Bryopsida/genética , Sistemas de Lectura Abierta , Proteoma , ARN Largo no Codificante , TranscriptomaRESUMEN
The rDNA clusters and flanking sequences on human chromosomes 13, 14, 15, 21 and 22 represent large gaps in the current genomic assembly. The organization and the degree of divergence of the human rDNA units within an individual nucleolar organizer region (NOR) are only partially known. To address this lacuna, we previously applied transformation-associated recombination (TAR) cloning to isolate individual rDNA units from chromosome 21. That approach revealed an unexpectedly high level of heterogeneity in human rDNA, raising the possibility of corresponding variations in ribosome dynamics. We have now applied the same strategy to analyze an entire rDNA array end-to-end from a copy of chromosome 22. Sequencing of TAR isolates provided the entire NOR sequence, including proximal and distal junctions that may be involved in nucleolar function. Comparison of the newly sequenced rDNAs to reference sequence for chromosomes 22 and 21 revealed variants that are shared in human rDNA in individuals from different ethnic groups, many of them at high frequency. Analysis infers comparable intra- and inter-individual divergence of rDNA units on the same and different chromosomes, supporting the concerted evolution of rDNA units. The results provide a route to investigate further the role of rDNA variation in nucleolar formation and in the empirical associations of nucleoli with pathology.
Asunto(s)
Cromosomas Humanos Par 22/genética , ADN Ribosómico/genética , Genoma Humano/genética , Región Organizadora del Nucléolo/genética , Nucléolo Celular/genética , Clonación Molecular , Heterogeneidad Genética , Genómica , Humanos , Anotación de Secuencia Molecular , Ribosomas/genéticaRESUMEN
The opioid receptor (OPR) family comprises the mu-, delta-, and kappa-opioid, and nociceptin receptors that belong to the superfamily of 7-transmembrane spanning G protein-coupled receptors (GPCRs). The mu-opioid receptor is the main target for clinically used opioid analgesics, and its biology has been extensively studied. The N-terminally truncated 6TM receptors isoform produced through alternative splicing of the OPRM1 gene displays unique signaling and analgesic properties, but it is unclear if other OPRs have the same ability. In this study, we have built a comprehensive map of alternative splicing events that produce 6TM receptor variants in all the OPRs and demonstrated their evolutionary conservation. We then obtained evidence for their translation through ribosomal footprint analysis. We discovered that N-terminally truncated 6TM GPCRs are rare in the human genome and OPRs are overrepresented in this group. Finally, we also observed a significant enrichment of 6TM GPCR genes among genes associated with pain, psychiatric disorders, and addiction. Understanding the biology of 6TM receptors and leveraging this knowledge for drug development should pave the way for novel therapies.
Asunto(s)
Empalme Alternativo/genética , Secuencia Conservada/genética , Receptores Opioides delta/genética , Receptores Opioides kappa/genética , Receptores Opioides mu/genética , Receptores Opioides/genética , Animales , Línea Celular Tumoral , Bases de Datos Genéticas , Variación Genética/genética , Humanos , Macaca , Ratones , Especificidad de la Especie , Receptor de NociceptinaRESUMEN
The effectiveness of oncolytic virotherapy in cancer treatment depends on several factors, including successful virus delivery to the tumor, ability of the virus to enter the target malignant cell, virus replication, and the release of progeny virions from infected cells. The multi-stage process is influenced by the efficiency with which the virus enters host cells via specific receptors. This review describes natural and artificial receptors for two oncolytic paramyxoviruses, nonpathogenic measles, and Sendai viruses. Cell entry receptors are proteins for measles virus (MV) and sialylated glycans (sialylated glycoproteins or glycolipids/gangliosides) for Sendai virus (SeV). Accumulated published data reviewed here show different levels of expression of cell surface receptors for both viruses in different malignancies. Patients whose tumor cells have low or no expression of receptors for a specific oncolytic virus cannot be successfully treated with the virus. Recent published studies have revealed that an expression signature for immune genes is another important factor that determines the vulnerability of tumor cells to viral infection. In the future, a combination of expression signatures of immune and receptor genes could be used to find a set of oncolytic viruses that are more effective for specific malignancies.
RESUMEN
Long noncoding RNAs (lncRNAs) are often associated with polysomes, indicating coding potential. However, only a handful of endogenous proteins encoded by putative lncRNAs have been identified and assigned a function. Here, we report the discovery of a putative gastrointestinal-tract-specific lncRNA (LINC00675) that is regulated by the pioneer transcription factor FOXA1 and encodes a conserved small protein of 79 amino acids which we termed FORCP (FOXA1-Regulated Conserved Small Protein). FORCP transcript is undetectable in most cell types but is abundant in well-differentiated colorectal cancer (CRC) cells where it functions to inhibit proliferation, clonogenicity, and tumorigenesis. The epitope-tagged and endogenous FORCP protein predominantly localizes to the endoplasmic reticulum (ER). In response to ER stress, FORCP depletion results in decreased apoptosis. Our findings on the initial characterization of FORCP demonstrate that FORCP is a novel, conserved small protein encoded by a mis-annotated lncRNA that regulates apoptosis and tumorigenicity in well-differentiated CRC cells.
Asunto(s)
Apoptosis/genética , Carcinogénesis/genética , Neoplasias Colorrectales/genética , Factor Nuclear 3-alfa del Hepatocito/metabolismo , ARN Largo no Codificante/metabolismo , Línea Celular Tumoral , Proliferación Celular , Neoplasias Colorrectales/patología , Retículo Endoplásmico/metabolismo , Estrés del Retículo Endoplásmico , Tracto Gastrointestinal/metabolismo , Genes Reporteros , Factor Nuclear 3-alfa del Hepatocito/genética , Humanos , Anotación de Secuencia Molecular , Especificidad de Órganos , ARN Largo no Codificante/genéticaAsunto(s)
Edición Génica/métodos , Hibridación de Ácido Nucleico/métodos , Artefactos , Disparidad de Par Base/genética , Disparidad de Par Base/fisiología , Sistemas CRISPR-Cas/genética , Sistemas CRISPR-Cas/fisiología , Repeticiones Palindrómicas Cortas Agrupadas y Regularmente Espaciadas/genética , Edición Génica/tendencias , Humanos , Cinética , TermodinámicaRESUMEN
The differences between high risk and low risk human papillomaviruses (HR-HPV and LR-HPV, respectively) that contribute to the tumorigenic potential of HR-HPV are not well understood but can be expected to involve the HPV oncoproteins, E6 and E7. We combine genome comparison and machine learning techniques to identify a previously unnoticed insert near the 3'-end of the E6 oncoprotein gene that is unique to HR-HPV. Analysis of the insert sequence suggests that it exerts a dual effect, by creating a PDZ domain-binding motif at the C-terminus of E6, as well as eliminating the overlap between the E6 and E7 coding regions in HR-HPV. We show that, as a result, the insert might enable coupled termination-reinitiation of the E6 and E7 genes, supported by motifs complementary to the human 18S rRNA. We hypothesize that the added functionality of E6 and positive regulation of E7 expression jointly account for the tumorigenic potential of HR-HPV.
Asunto(s)
Proteínas Oncogénicas Virales , Papillomaviridae , Carcinogénesis , Humanos , Papillomaviridae/genética , Papillomaviridae/patogenicidad , Medición de RiesgoRESUMEN
In response to low levels of magnesium (Mg2+ ), the PhoQP two component system induces the transcription of two convergent genes, one encoding a 31-amino acid protein denoted MgtS and the second encoding a small, regulatory RNA (sRNA) denoted MgrR. Previous studies showed that the MgtS protein interacts with and stabilizes the MgtA Mg2+ importer to increase intracellular Mg2+ levels, while the MgrR sRNA base pairs with the eptB mRNA thus affecting lipopolysaccharide modification. Surprisingly, we found overexpression of the MgtS protein also leads to induction of the PhoRB regulon. Studies to understand this activation showed that MgtS forms a complex with a second protein, PitA, a cation-phosphate symporter. Given that the additive effect of ∆mgtA and ∆mgtS mutations on intracellular Mg2+ concentrations seen previously is lost in the ∆pitA mutant, we suggest that MgtS binds to and prevents Mg2+ leakage through PitA under Mg2+ -limiting conditions. Consistent with a detrimental role of PitA in low Mg2+ , we also observe MgrR sRNA repression of PitA synthesis. Thus, PhoQP induces the expression of two convergent small genes in response to Mg2+ limitation whose products act to modulate PitA at different levels to increase intracellular Mg2+ .
Asunto(s)
Proteínas de Escherichia coli/metabolismo , Escherichia coli/metabolismo , Regulación Bacteriana de la Expresión Génica , Magnesio/metabolismo , Proteínas de la Membrana/metabolismo , Proteínas de Transporte de Fosfato/metabolismo , ARN Pequeño no Traducido/metabolismo , Escherichia coli/genética , Redes Reguladoras de GenesRESUMEN
The delta-opioid receptor (DOPr) participates in mediating the effects of opioid analgesics. However, no selective agonists have entered clinical care despite potential to ameliorate many neurological and psychiatric disorders. In an effort to address the drug development challenges, the functional contribution of receptor isoforms created by alternative splicing of the three-exonic coding gene, OPRD1, has been overlooked. We report that the gene is transcriptionally more diverse than previously demonstrated, producing novel protein isoforms in humans and mice. We provide support for the functional relevance of splice variants through context-dependent expression profiling (tissues, disease model) and conservation of the transcriptional landscape in closely related vertebrates. The conserved alternative transcriptional events have two distinct patterns. First, cassette exon inclusions between exons 1 and 2 interrupt the reading frame, producing truncated receptor fragments comprising only the first transmembrane (TM) domain, despite the lack of exact exon orthologues between distant species. Second, a novel promoter and transcriptional start site upstream of exon 2 produces a transcript of an N-terminally truncated 6TM isoform. However, a fundamental difference in the exonic landscaping as well as translation and translation products poses limits for modelling the human DOPr receptor system in mice.
Asunto(s)
Empalme Alternativo/genética , Receptores Opioides delta/genética , Animales , Línea Celular Tumoral , Dolor Crónico/genética , Dolor Crónico/patología , Secuencia Conservada , Modelos Animales de Enfermedad , Evolución Molecular , Feto/metabolismo , Humanos , Masculino , Ratones Endogámicos C57BL , Biosíntesis de Proteínas , Isoformas de Proteínas/genética , Isoformas de Proteínas/metabolismo , ARN Mensajero/genética , ARN Mensajero/metabolismo , Receptores Opioides delta/metabolismo , Ribosomas/metabolismoRESUMEN
Off-target oligoprobe's interaction with partially complementary nucleotide sequences represents a problem for many bio-techniques. The goal of the study was to identify oligoprobe sequence characteristics that control the ratio between on-target and off-target hybridization. To understand the complex interplay between specific and genome-wide off-target (cross-hybridization) signals, we analyzed a database derived from genomic comparison hybridization experiments performed with an Affymetrix tiling array. The database included two types of probes with signals derived from (i) a combination of specific signal and cross-hybridization and (ii) genomic cross-hybridization only. All probes from the database were grouped into bins according to their sequence characteristics, where both hybridization signals were averaged separately. For selection of specific probes, we analyzed the following sequence characteristics: vulnerability to self-folding, nucleotide composition bias, numbers of G nucleotides and GGG-blocks, and occurrence of probe's k-mers in the human genome. Increases in bin ranges for these characteristics are simultaneously accompanied by a decrease in hybridization specificity-the ratio between specific and cross-hybridization signals. However, both averaged hybridization signals exhibit growing trends along with an increase of probes' binding energy, where the hybridization specific signal increases significantly faster in comparison to the cross-hybridization. The same trend is evident for the S function, which serves as a combined evaluation of probe binding energy and occurrence of probe's k-mers in the genome. Application of S allows extracting a larger number of specific probes, as compared to using only binding energy. Thus, we showed that high values of specific and cross-hybridization signals are not mutually exclusive for probes with high values of binding energy and S. In this study, the application of a new set of sequence characteristics allows detection of probes that are highly specific to their targets for array design and other bio-techniques that require selection of specific probes.
Asunto(s)
Hibridación de Ácido Nucleico , Análisis de Secuencia por Matrices de Oligonucleótidos/métodos , Secuencia de Bases , Bases de Datos Genéticas , Genoma , HumanosRESUMEN
Despite the key role of the human ribosome in protein biosynthesis, little is known about the extent of sequence variation in ribosomal DNA (rDNA) or its pre-rRNA and rRNA products. We recovered ribosomal DNA segments from a single human chromosome 21 using transformation-associated recombination (TAR) cloning in yeast. Accurate long-read sequencing of 13 isolates covering â¼0.82 Mb of the chromosome 21 rDNA complement revealed substantial variation among tandem repeat rDNA copies, several palindromic structures and potential errors in the previous reference sequence. These clones revealed 101 variant positions in the 45S transcription unit and 235 in the intergenic spacer sequence. Approximately 60% of the 45S variants were confirmed in independent whole-genome or RNA-seq data, with 47 of these further observed in mature 18S/28S rRNA sequences. TAR cloning and long-read sequencing enabled the accurate reconstruction of multiple rDNA units and a new, high-quality 44 838 bp rDNA reference sequence, which we have annotated with variants detected from chromosome 21 of a single individual. The large number of variants observed reveal heterogeneity in human rDNA, opening up the possibility of corresponding variations in ribosome dynamics.