Proteomic identification of putative microRNA394 target genes in Arabidopsis thaliana identifies MLP family members critical for normal development

Expression of the F-Box protein Leaf Curling Responsiveness (LCR) is regulated by microRNA, miR394, and alterations to this interplay in Arabidopsis thaliana produce defects in leaf polarity and shoot apical meristem (SAM) organisation. Although the miR394-LCR node has been documented in Arabidopsis, the identification of proteins targeted by LCR F-box itself has proven problematic. Here, a proteomic analysis of shoot apices from plants with altered LCR levels identified a member of the Major Latex Protein (MLP) family gene as a potential LCR F-box target. Bioinformatic and molecular analyses also suggested that other MLP family members are likely to be targets for this post-translational regulation. Direct interaction between LCR F-Box and MLP423 was validated. Additional MLP members had reduction in protein accumulation, in varying degrees, mediated by LCR F-Box. Transgenic Arabidopsis lines, in which MLP28 expression was reduced through an artificial miRNA technology, displayed severe developmental defects, including changes in leaf patterning and morphology, shoot apex defects, and eventual premature death. These phenotypic characteristics resemble those of Arabidopsis plants modified to over-express LCR. Taken together, the results demonstrate that MLPs are driven to degradation by LCR, and indicate that MLP gene family is target of miR394-LCR regulatory node, representing potential targets for directly post-translational regulation mediated by LCR F-Box. In addition, MLP28 family member is associated with the LCR regulation that is critical for normal Arabidopsis development.

Genome-wide expression profiling identifies deregulated miRNAs in malignant astrocytoma

Malignant astrocytoma includes anaplastic astrocytoma (grade III) and glioblastoma (grade IV). Among them, glioblastoma is the most common primary brain tumor with dismal responses to all therapeutic modalities. We performed a large-scale, genome-wide microRNA (miRNA) (n=756) expression profiling of 26 glioblastoma, 13 anaplastic astrocytoma and 7 normal brain samples with an aim to find deregulated miRNA in malignant astrocytoma. We identified several differentially regulated miRNAs between these groups, which could differentiate glioma grades and normal brain as recognized by PCA. More importantly, we identified a most discriminatory 23-miRNA expression signature, by using PAM, which precisely distinguished glioblastoma from anaplastic astrocytoma with an accuracy of 95%. The differential expression pattern of nine miRNAs was further validated by real-time RT-PCR on an independent set of malignant astrocytomas (n-72) and normal samples (n=7). Inhibition of two glioblastoma-upregulated miRNAs (miR-21 and miR-23a) and exogenous overexpression of two glioblastoma-downregulated miRNAs (miR-218 and miR-219-5p) resulted in reduced soft agar colony formation but showed varying effects on cell proliferation and chemosensitivity. Thus we have identified the miRNA expression signature for malignant astrocytoma, in particular glioblastoma, and showed the miRNA involvement and their importance in astrocytoma development. Modern Pathology (2010) 23, 1404-1417; doi:10.1038/modpathol.2010.135; published online 13 August 2010

Genomic alterations in Myeloproliferative Neoplasms and Myelodysplastic Syndromes

Myeloproliferative neoplasms (MPN) and myelodysplastic syndromes (MDS) are a heterogeneous group of clonal hematopoietic disorders whose etiology and molecular pathogenesis are poorly understood. During the past decade, enormous developments in microarray technology and bioinformatics methods have made it possible to mine novel molecular alterations in a large number of malignancies, including MPN and MDS, which has facilitated the detection of new prognostic, predictive and therapeutic biomarkers for disease stratification. By applying novel microarray techniques, we profiled copy number alterations and microRNA (miRNA) expression changes in bone marrow aspirate and blood samples. In addition, we set up and validated an miRNA expression test for bone marrow core biopsies in order to utilize the large archive material available in many laboratories. We also tested JAK2 mutation status and compare it with the in vitro growth pattern of hematologic progenitors cells. In the study focusing on 100 MPN cases, we detected a Janus kinase 2 (JAK2) mutation in 71 cases. We observed spontaneous erythroid colony growth in all mutation-positive cases in addition to nine mutation negative cases. Interestingly, seven JAK2V167F negative ET cases showed spontaneous megakaryocyte colony formation, one case of which also harbored a myeloproliferative leukemia virus oncogene (MPL) mutation. We studied copy number alterations in 35 MPN and 37 MDS cases by using oligonucleotide-based array comparative hybridization (array CGH). Only one essential thrombocythemia (ET) case presented copy number alterations in chromosomes 1q and 13q. In contrast, MDS cases were characterized by numerous novel cryptic chromosomal aberrations with the most common copy number losses at 5q21.3q33.1 and 7q22.1q33, while the most common copy number gain was trisomy 8. As for the study of the bone marrow core biopsy samples, we showed that even though these samples were embedded in paraffin and underwent decalcification, they were reliable sources of miRNA and suitable for array expression analysis. Further, when studying the miRNA expression profiles of the 19 MDS cases, we found that, compared to controls, two miRNAs (one human Epstein-Barr virus (miR-BART13) miRNA and one human (has-miR-671-5p) miRNA) were downregulated, whereas two other miRNAs (hsa-miR-720 and hsa-miR-21) were upregulated. However, we could find no correlation between copy number alterations and microRNA expression when integrating these two data. This thesis brings to light new information about genomic changes implicated in the development of MPN and MDS, and also underlines the power of applying genome-wide array screening techniques in neoplasias. Rapid advances in molecular techniques and the integration of different genomic data will enable the discovery of the biological contexts of many complex disorders, including myeloid neoplasias.

Circulating miRNA profile in HCV infected serum: novel insight into pathogenesis

Changes in circulating miRNA profiles have been associated with different diseases. Here we demonstrate the circulating miRNA profile in serum of HCV infected individuals using a microRNA array that profiles the expression of 940 miRNAs. Serum samples from two HCV genotype -1 and two HCV genotype -3 infected individuals were compared with healthy controls. Expression levels of miR-134, miR-198, miR-320c and miR-483-5p that were commonly upregulated in case of both genotypes were validated in 36 individual patient serum samples. Serum miR-134, miR-320c and miR-483-5p were significantly upregulated during HCV infection. miR-320c and miR-483-5p were also upregulated in HCV-JFH1 infected cells and cell culture supernatant. Pathway analysis of putative target genes of these miRNAs indicated involvement of PI3K-Akt, NFKB and MAPK signaling pathways. Results revealed novel insights on the role of circulating miRNAs in mediating pathogenesis in HCV-infected cells.

Ac2PIM-responsive miR-150 and miR-143 Target Receptor-interacting Protein Kinase 2 and Transforming Growth Factor Beta-activated Kinase 1 to Suppress NOD2-induced Immunomodulators

Specific and coordinated regulation of innate immune receptor-driven signaling networks often determines the net outcome of the immune responses. Here, we investigated the cross-regulation of toll-like receptor (TLR)2 and nucleotide-binding oligomerization domain (NOD)2 pathways mediated by Ac2PIM, a tetra-acylated form of mycobacterial cell wall component and muramyl dipeptide (MDP), a peptidoglycan derivative respectively. While Ac2PIM treatment of macrophages compromised their ability to induce NOD2-dependent immunomodulators like cyclooxygenase (COX)-2, suppressor of cytokine signaling (SOCS)-3, and matrix metalloproteinase (MMP)-9, no change in the NOD2-responsive NO, TNF-alpha, VEGF-A, and IL-12 levels was observed. Further, genome-wide microRNA expression profiling identified Ac2PIM-responsive miR-150 and miR-143 to target NOD2 signaling adaptors, RIP2 and TAK1, respectively. Interestingly, Ac2PIM was found to activate the SRC-FAK-PYK2-CREB cascade via TLR2 to recruit CBP/P300 at the promoters of miR-150 and miR-143 and epigenetically induce their expression. Loss-of-function studies utilizing specific miRNA inhibitors establish that Ac2PIM, via the miRNAs, abrogate NOD2-induced PI3K-PKC delta-MAPK pathway to suppress beta-catenin-mediated expression of COX-2, SOCS-3, and MMP-9. Our investigation has thus underscored the negative regulatory role of Ac2PIM-TLR2 signaling on NOD2 pathway which could broaden our understanding on vaccine potential or adjuvant utilities of Ac2PIM and/or MDP.

Inhibition of p300 lysine acetyltransferase activity by luteolin reduces tumor growth in head and neck squamous cell carcinoma (HNSCC) xenograft mouse model

Chromatin acetylation is attributed with distinct functional relevance with respect to gene expression in normal and diseased conditions thereby leading to a topical interest in the concept of epigenetic modulators and therapy. We report here the identification and characterization of the acetylation inhibitory potential of an important dietary flavonoid, luteolin. Luteolin was found to inhibit p300 acetyltransferase with competitive binding to the acetyl CoA binding site. Luteolin treatment in a xenografted tumor model of head and neck squamous cell carcinoma (HNSCC), led to a dramatic reduction in tumor growth within 4 weeks corresponding to a decrease in histone acetylation. Cells treated with luteolin exhibit cell cycle arrest and decreased cell migration. Luteolin treatment led to an alteration in gene expression and miRNA profile including up-regulation of p53 induced miR-195/215, let7C; potentially translating into a tumor suppressor function. It also led to down regulation of oncomiRNAs such as miR-135a, thereby reflecting global changes in the microRNA network. Furthermore, a direct correlation between the inhibition of histone acetylation and gene expression was established using chromatin immunoprecipitation on promoters of differentially expressed genes. A network of dysregulated genes and miRNAs was mapped along with the gene ontology categories, and the effects of luteolin were observed to be potentially at multiple levels: at the level of gene expression, miRNA expression and miRNA processing.

High-Throughput Sequencing of microRNAs in Peripheral Blood Mononuclear Cells: Identification of Potential Weight Loss Biomarkers

INTRODUCTION: MicroRNAs (miRNAs) are being increasingly studied in relation to energy metabolism and body composition homeostasis. Indeed, the quantitative analysis of miRNAs expression in different adiposity conditions may contribute to understand the intimate mechanisms participating in body weight control and to find new biomarkers with diagnostic or prognostic value in obesity management. OBJECTIVE: The aim of this study was the search for miRNAs in blood cells whose expression could be used as prognostic biomarkers of weight loss. METHODS: Ten Caucasian obese women were selected among the participants in a weight-loss trial that consisted in following an energy-restricted treatment. Weight loss was considered unsuccessful when <5% of initial body weight (non-responders) and successful when >5% (responders). At baseline, total miRNA isolated from peripheral blood mononuclear cells (PBMC) was sequenced with SOLiD v4. The miRNA sequencing data were validated by RT-PCR. RESULTS: Differential baseline expression of several miRNAs was found between responders and non-responders. Two miRNAs were up-regulated in the non-responder group (mir-935 and mir-4772) and three others were down-regulated (mir-223, mir-224 and mir-376b). Both mir-935 and mir-4772 showed relevant associations with the magnitude of weight loss, although the expression of other transcripts (mir-874, mir-199b, mir-766, mir-589 and mir-148b) also correlated with weight loss. CONCLUSIONS: This research addresses the use of high-throughput sequencing technologies in the search for miRNA expression biomarkers in obesity, by determining the miRNA transcriptome of PBMC. Basal expression of different miRNAs, particularly mir-935 and mir-4772, could be prognostic biomarkers and may forecast the response to a hypocaloric diet.

Revealing, illuminating, and modifying proteins in human diseases using noncanonical amino acids

To better understand human diseases, much recent work has focused on proteins to either identify disease targets through proteomics or produce therapeutics via protein engineering. Noncanonical amino acids (ncAAs) are tools for altering the chemical and physical properties of proteins, providing a facile strategy not only to label proteins but also to engineer proteins with novel properties. My thesis research has focused on the development and applications of noncanonical amino acids in identifying, imaging, and engineering proteins for studying human diseases. Chapter 1 introduces the concept of ncAAs and reveals insights to how I chose my thesis projects.

ncAAs have been incorporated to tag and enrich newly synthesized proteins for mass spectrometry through a method termed BONCAT, or bioorthogonal noncanonical amino acid tagging. Chapter 2 describes the investigation of the proteomic response of human breast cancer cells to induced expression of tumor suppressor microRNA miR-126 by combining BONCAT with another proteomic method, SILAC or stable isotope labeling by amino acids in cell culture. This proteomic analysis led to the discovery of a direct target of miR-126, shedding new light on its role in suppressing cancer metastasis.

In addition to mass spectrometry, ncAAs can also be utilized to fluorescently label proteins. Chapter 3 details the synthesis of a set of cell-permeant cyclooctyne probes and demonstration of selective labeling of newly synthesized proteins in live mammalian cells using azidohomoalanine. Similar to live cell imaging, the ability to selectively label a particular cell type within a mixed cell population is important to interrogating many biological systems, such as tumor microenvironments. By taking advantage of the metabolic differences between cancer and normal cells, Chapter 5 discusses efforts to develop selective labeling of cancer cells using a glutamine analogue.

Furthermore, Chapter 4 describes the first demonstration of global replacement at polar amino acid positions and its application in developing an alternative PEGylation strategy for therapeutic proteins. Polar amino acids typically occupy solvent-exposed positions on the protein surface, and incorporation of noncanonical amino acids at these positions should allow easier modification and cause less perturbation compared to replacements at the interior positions of proteins.

Genetic Variants in MiRNA Processing Genes and Pre-MiRNAs Are Associated with the Risk of Chronic Lymphocytic Leukemia

Genome wide association studies (GWAS) have identified several low-penetrance susceptibility alleles in chronic lymphocytic leukemia (CLL). Nevertheless, these studies scarcely study regions that are implicated in non-coding molecules such as microRNAs (miRNAs). Abnormalities in miRNAs, as altered expression patterns and mutations, have been described in CLL, suggesting their implication in the development of the disease. Genetic variations in miRNAs can affect levels of miRNA expression if present in pre-miRNAs and in miRNA biogenesis genes or alter miRNA function if present in both target mRNA and miRNA sequences. Therefore, the present study aimed to evaluate whether polymorphisms in pre-miRNAs, and/or miRNA processing genes contribute to predisposition for CLL. A total of 91 SNPs in 107 CLL patients and 350 cancer-free controls were successfully analyzed using TaqMan Open Array technology. We found nine statistically significant associations with CLL risk after FDR correction, seven in miRNA processing genes (rs3805500 and rs6877842 in DROSHA, rs1057035 in DICER1, rs17676986 in SND1, rs9611280 in TNRC6B, rs784567 in TRBP and rs11866002 in CNOT1) and two in pre-miRNAs (rs11614913 in miR196a2 and rs2114358 in miR1206). These findings suggest that polymorphisms in genes involved in miRNAs biogenesis pathway as well as in pre-miRNAs contribute to the risk of CLL. Large-scale studies are needed to validate the current findings.

Investigação de alterações em genes de microRNAs expressos no cérebro como causa de deficiência intelectual ligada ao cromossomo X

A Deficiência Intelectual (DI) é uma condição definida como um funcionamento intelectual significativamente prejudicado, expresso juntamente com limitações em pelo menos duas áreas do comportamento adaptativo que se manifestam antes dos 18 anos de idade. A prevalência estimada da DI na população em geral é de 2-3% e um número expressivo de casos permanece sem um diagnóstico definitivo. Há um consenso geral de que a DI é mais comum em indivíduos do sexo masculino em relação aos do sexo feminino. Entre as explicações para este excesso está a concentração de genes específicos para a habilidade cognitiva no cromossomo X. MicroRNAs (miRNAs) são pequenas moléculas de RNA não codificador que modulam a expressão gênica pós-transcricional de RNAs mensageiros alvo. Recentemente, estudos têm demonstrado a importância essencial dos miRNAs para o desenvolvimento e funcionamento cerebrais e sabe-se que o cromossomo X tem uma alta densidade de genes de miRNAs. Neste contexto, os miRNAs são candidatos potenciais como fatores genéticos envolvidos na Deficiência Intelectual Ligada ao X (DILX). Neste estudo, foram analisadas as regiões genômicas de 17 genes de miRNAs expressos no cérebro localizados no cromossomo X, com o objetivo de investigar o possível envolvimento de variantes na sequência destes miRNAs na DILX. Para este fim, selecionamos amostras de DNA genômico (sangue periférico) de 135 indivíduos do sexo masculino portadores de DI sugestiva de DILX de um grupo de mais de 1.100 pacientes com DI encaminhados ao Serviço de Genética Humana da UERJ. O critério de inclusão para este estudo era de que os probandos apresentassem um ou mais parentes do sexo masculino afetados pela DI que fossem interligados por via materna. As amostras de DNA dos pacientes foram amplificadas utilizando a técnica de reação em cadeia da polimerase, seguida por purificação e sequenciamento direto pelo método de Sanger dos fragmentos amplificados. Para avaliar a conservação dos 17 miRNAs foi realizada uma análise filogenética in silico incluindo sequências dos miRNAs selecionados de humanos e de outras 8 espécies de primatas estreitamente relacionadas. Não foram encontradas alterações nas sequências nos genes de 17 miRNAs analisados, mesmo diante do padrão genético altamente heterogêneo da população brasileira. Adicionalmente, a análise filogenética destes miRNAs revelou uma alta conservação entre as espécies comparadas. Considerando o papel dos miRNAs como reguladores da expressão gênica, a ausência de alterações e a alta conservação entre primatas sugerem uma forte pressão seletiva sobre estas moléculas, reforçando a sua importância funcional para o organismo em geral. Apesar de não termos encontrado variantes de sequência nos miRNAs estudados, o envolvimento de miRNAs na DI não pode ser completamente descartado. Alterações fora da molécula de miRNA precursor, nos fatores de processamento, nos sítios alvo e variações no número de cópias de genes de miRNAs podem implicar em alteração na expressão dos miRNAs e, consequentemente, na funcionalidade do miRNA maduro. Sendo assim, uma análise sistemática da expressão de miRNAs em pacientes com DILX é urgentemente necessária, a fim de desvendar novos genes/mecanismos moleculares relacionados a esta condição.

MicroRNAs in breast cancer progression and DNA damage response

Os tumores de mama são caracterizados pela sua alta heterogeneidade. O câncer de mama é uma doença complexa, que possui o seu desenvolvimento fortemente influenciado por fatores ambientais, combinada a uma progressiva acumulação de mutações genéticas e desregulação epigenética de vias críticas. Alterações nos padrões de expressão gênica podem ser resultado de uma desregulação no controle de eventos epigenéticos, assim como, na regulação pós-transcricional pelo mecanismo de RNA de interferência endógeno via microRNA (miRNA). Estes eventos são capazes de levar à iniciação, à promoção e à manutenção da carcinogênese, como também ter implicações no desenvolvimento da resistência à terapia Os miRNAs formam uma classe de RNAs não codificantes, que durante os últimos anos surgiram como um dos principais reguladores da expressão gênica, através da sua capacidade de regular negativamente a atividade de RNAs mensageiros (RNAms) portadores de uma seqüencia parcialmente complementar. A importância da regulação mediada por miRNAs foi observada pela capacidade destas moléculas em regular uma vasta gama de processos biológicos incluindo a proliferação celular, diferenciação e a apoptose. Para avaliar a expressão de miRNAs durante a progressão tumoral, utilizamos como modelo experimental a série 21T que compreende 5 linhagens celulares originárias da mesma paciente diagnosticada com um tumor primário de mama do tipo ErbB2 e uma posterior metástase pulmonar. Essa série é composta pela linhagem obtida a partir do tecido normal 16N, pelas linhagens correspondentes ao carcinoma primário 21PT e 21NT e pelas linhagens obtidas um ano após o diagnóstico inicial, a partir da efusão pleural no sítio metastatico 21MT1 e 21MT2. O miRNAoma da série 21T revelou uma redução significativa nos níveis de miR-205 e nos níveis da proteina e-caderina e um enriquecimento do fator pró-metastático ZEB-1 nas células 21MT. Considerando a importância dos miRNAs na regulação da apoptose, e que a irradiação em diferentes espectros é comumente usada em procedimentos de diagnóstico como mamografia e na radioterapia, avaliamos a expressão de miRNAs após irradiação de alta e baixa energia e do tratamento doxorrubicina. Para os ensaios foram utilizados as linhagens não tumorais MCF-10A e HB-2 e as linhagens de carcinoma da mama MCF-7 e T-47D. Observou-se que raios-X de baixa energia são capazes de promover quebras na molécula do DNA e apoptose assim como, alterar sensivelmente miRNAs envolvidos nessas vias como o let-7a, miR-34a e miR-29b. No que diz respeito à resposta a danos genotóxicos, uma regulação positiva sobre a expressão de miR-29b, o qual em condições normais é regulado negativamente foi observada uma regulação positiva sobre miR-29b expressão após todos os tratamentos em células tumorais. Nossos resultados indicam que miR-29b é um possível biomarcador de estresse genotóxico e que miR-205 pode participar no potencial metastático das células 21T.

Genome-wide computational identification of microRNAs and their targets in the deep-branching eukaryote Giardia lamblia

Using a combined computational program. we identified 50 potential microRNAs (miRNAs) in Giardia lamblia. one of the most primitive unicellular eukaryotes. These miRNAs are unique to G. lamblia and no homologues have been found in other organisms; miRNAs.

Identifying cancer subtypes in glioblastoma by combining genomic, transcriptomic and epigenomic data

We present a nonparametric Bayesian method for disease subtype discovery in multi-dimensional cancer data. Our method can simultaneously analyse a wide range of data types, allowing for both agreement and disagreement between their underlying clustering structure. It includes feature selection and infers the most likely number of disease subtypes, given the data. We apply the method to 277 glioblastoma samples from The Cancer Genome Atlas, for which there are gene expression, copy number variation, methylation and microRNA data. We identify 8 distinct consensus subtypes and study their prognostic value for death, new tumour events, progression and recurrence. The consensus subtypes are prognostic of tumour recurrence (log-rank p-value of $3.6 \times 10^{-4}$ after correction for multiple hypothesis tests). This is driven principally by the methylation data (log-rank p-value of $2.0 \times 10^{-3}$) but the effect is strengthened by the other 3 data types, demonstrating the value of integrating multiple data types. Of particular note is a subtype of 47 patients characterised by very low levels of methylation. This subtype has very low rates of tumour recurrence and no new events in 10 years of follow up. We also identify a small gene expression subtype of 6 patients that shows particularly poor survival outcomes. Additionally, we note a consensus subtype that showly a highly distinctive data signature and suggest that it is therefore a biologically distinct subtype of glioblastoma. The code is available from https://sites.google.com/site/multipledatafusion/

大脑表达 miRNA 在大脑进化中的作用及其遗传机制初步研究

人类区别于其它动物的最本质特征是其拥有其他动物所无法比拟的大脑容量及高 级的认知能力。即使与其近亲－非人灵长类相比，人也拥有比非人灵长类大好几倍的脑 容量和更为发达的认知能力。现在一般认为，人类大脑的形成是适应性选择(达尔文正 选择)的结果。但是到目前为止，对人类起源过程中大脑容量增大及认知能力提高的遗 传学机制却知之甚少。以前的研究表明，几个与大脑发育相关的蛋白质编码基因在人类 起源中受到了正向选择。同时，也有证据表明人类大脑的进化也可能是基因表达调控变 化的结果。因此寻找人与非人灵长类大脑表达基因的调控差异或许能够进一步为人与非 人灵长类为何有如此巨大的差异提供分子生物学水平上的解释。 MicroRNA(miRNA)是一类在转录后水平调控基因表达的不编码蛋白质的小 RNA(长度20－24个碱基)。通过调控靶基因的表达，miRNA参与了众多的生理过程。而 很多大脑表达的miRNA的表达量在大脑发育过程中呈显著的变化，这表明miRNA参与 了大脑的发育调控。由此可知，miRNA对其靶基因调控效率的改变很可能引起大脑发育 调控的改变。 本文通过寻找大脑表达的保守microRNA 及对其靶基因的预测，同时用比较基因组 学的方法发现：有很多microRNA 的靶基因3’UTR 的靶位点中存在人类特异突变位点。 我们推测这些人类特异突变位点可能改变microRNA 对其靶基因的表达调控效应，而调 控效应的改变可能在进化过程中对认知能力的提高发挥重要作用。利用体外报告基因系 统，我们发现有几个预测的靶基因是对应miRNA 的真实靶基因。其中，miR-127 的靶 基因（SEMA3F）3’UTR 的靶位点中所含的一个人类特有的突变增强了miR-127 对 SEMA3F 的调控效率。将该位点突变回复为黑猩猩的位点使得miR-127 对SEMA3F 的 调控效率降低到黑猩猩的水平。这表明miR-127 对SEMA3F 的调控效率的改变确实是 由该人类特异突变位点引起的。我们提供了人类特异突变位点能够引起miR-127 对SEMA3F 的调控效率的改变的体外证据，但是体内的调控模式是否如此尚需进一步的工 作。 总之，本文通过体外试验表明，miRNA靶基因3’UTR的序列变异具有功能效应，它 有可能是人类中枢神经系统在起源和演化中发挥关键作用的重要遗传机制之一。

人类基因组中保守二级结构的纯净化选择及其在转录调控网络中的作用

保守序列是一种跨物种保守的基因组序列，而且绝大多数为非蛋白编码序 列。保守序列在人类遗传疾病中发挥着重要作用。其中，一部分保守序列能够 折叠形成二级结构。已鉴定的一些保守二级结构编码一些RNA 分子，如 microRNA、RNA 编辑序列和组蛋白mRNA 3’端非翻译区茎环结构等。但是，对 于绝大部分的保守二级结构，它们的生物学功能以及作用于它们上面的进化作 用力依然是未知的。 群体的SNP 数据在分析序列上的进化作用力时非常有效。SNP 在群体中的 频率会因为受到不同的进化作用力而表现出差异，而与其是否位于基因组中的 突变热点无关。对于受纯净化选择作用的SNP，它们的频率一般会比中性SNP 具有低的新生型等位基因频率（DAF）。我们运用生物信息学的方法，在人类基 因组保守二级结构中找到746 个SNP。这746 个SNP 与基因组其它区段的SNP 在突变模式上并不存在显著差异，在保守二级结构内同样存在突变热点。通过 与侧翼序列SNP 的分布比较发现，保守二级结构上SNP 密度约为其侧翼序列的 2/3。相比于侧翼序列SNP，有更高比例的保守二级结构SNP 具有低的DAF 值。 这些结果提示，有很多保守二级结构上的SNP 因为受到纯净化选择作用而在现 代人群中被剔除了。保守二级结构与侧翼序列在SNP 密度和DAF 上的差异要高 于保守序列与非保守序列之间的差异，提示保守二级结构是受到纯净化选择作 用最为严格的一类保守序列。我们发现，在保守二级结构内部，纯净化选择作 用的强度也有差异。茎区比环区具有更低的SNP 密度，而且有更高比例的茎区 SNP 具有低的DAF 值。这个结果提示，保守二级结构上的纯净化选择力主要作 用于茎区上的位点。我们推测，这可能是茎区上的突变往往比环区的突变对二级结构的造成更大的影响导致的。 我们通过寻找保守二级结构与转录因子SOX2、OCT4、NANOG、SUZ12 和C-MYC 结合位点之间的重叠，还分析了保守二级结构在转录调控网络中的作用。结果 显示，很多保守二级结构是作为转录因子的结合位点调控了许多与发育相关的 转录因子编码基因的表达。转录因子与保守二级结构之间的结合模式非常复杂， 可以有多个转录因子结合到同一个保守二级结构上，也可以是一个转录因子结 合到自身编码基因相关的保守二级结构上。不同的转录因子和保守二级结构结 合可以主导靶基因的特异模式，当绝大多数相关的保守二级结构与SUZ12 结合 时，基因表达受到抑制，而当绝大多数相关的保守二级结构不与SUZ12 结合时， 基因表达受到激活。在转录调控网络中，约有30%的保守二级结构是作为启动 子来调控基因的表达。因为转录因子SOX2、OCT4、NANOG、SUZ12 和C-MYC 仅仅 只结合到很小一部分保守二级结构上，提示可能还有更多的转录因子会结合到 保守二级结构上。因此，保守二级结构介导的转录调控网络要比目前已知的复 杂得多。