$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Mycobacterium bovis A cepa BCG (bacilos de Calmette et Guérin) 1173P2 em crescimento exponencial foi submetida a uma série temporal (0, 4, 10 e 20 dias) de falta de nutrientes, seguida por uma ressuscitação de 6 dias em meio rico em nutrientes, conforme apresentado anteriormente em Hu et al.7. Pequenos RNAs foram isolados da cultura bacteriana, com três repetições biológicas, em cada um dos cinco pontos de tempo designados. As bibliotecas Illumina foram construídas usando o fluxo de trabalho de preparação da biblioteca AQRNA-seq descrito acima (Figura 1), seguido de sequenciamento em um sequenciador no BioMicro Center do Massachusetts Institute of Technology. Os dados de sequenciamento foram então processados usando o pipeline de análise de dados AQRNA-seq (Figura 2) personalizado para quantificação da abundância de tRNA.
Após a amplificação por PCR da biblioteca de cDNA com os primers de sequenciamento, observou-se a presença de produtos de PCR com tamanho de 175 pares de bases (pb) em todas as amostras (Figura 3A), sugerindo a formação de dímeros de primer. Para mitigar o transporte de dímeros de primer, produtos de PCR com tamanho superior a 195 pb foram extirpados do gel e purificados ( Figura 3B ).
As leituras de sequência filtradas e aparadas de qualidade foram mapeadas para uma biblioteca de sequência de referência personalizada, incluindo os isoaceitadores de 45 tRNA, o padrão interno e as sequências de controle (ou seja, 23S rRNA, 16S rRNA, 5S rRNA, rnpB e ssr). Os isoaceitadores de tRNA representaram 10,5% a 40,2% do total de leituras mapeadas de uma determinada amostra e mostraram uma abundância muito maior do que as sequências de controle (Figura 4). É importante ressaltar que as proporções de leitura relativamente baixas dos isoaceitadores de tRNA podem ser atribuídas à maior abundância relativa dos padrões internos. Portanto, as proporções de leitura dos isoaceitadores de tRNA em relação aos padrões internos (Figura 4, blocos de cor rosa vs verde) podem ser controladas pelo operador, por meio do ajuste fino da quantidade de padrão interno cravado na reação.
Os dados brutos de abundância de tRNA foram normalizados usando o método da mediana das razões implementado com a versão do pacote DESeq2 (doravante denominada v) 1.36.016 no R Statistical Programming Environment (doravante denominado R) v 4.2.117. Após a normalização, é alcançado um cenário quantitativo de isoaceitadores de tRNA em Mycobacterium bovis BCG durante um curso de tempo de falta de nutrientes e ressuscitação (Figura 5).
Para revelar grupos distintos de amostras com diferentes fenótipos com base em padrões na abundância de isoaceitadores de tRNA, a Análise de Componentes Principais (PCA) foi realizada nos dados de abundância de tRNA normalizados usando o pacote de estatísticas v 4.2.117 em R (Figura 6). A análise distinguiu amostras de dia de fome 0 e dia de ressuscitação 6 de amostras de dias de fome 4, 10 e 20, sugerindo uma diferença considerável na paisagem de tRNA de Mycobacterium bovis BCG cultivado em meio privado de nutrientes e meio rico em nutrientes.
Para traçar o perfil da dinâmica da abundância de cada isoaceitador de tRNA nos cinco pontos de tempo designados, a análise de expressão diferencial foi realizada nos dados de abundância de tRNA normalizados usando o pacote DESeq2 v 1.36.0 em R (Figura 7). A análise revelou que 17 das 20 famílias de isoaceitadores continham isoaceitadores que foram diferencialmente expressos (ou seja, significativamente regulados para cima ou para baixo) em pelo menos um dos pontos de tempo, sugerindo um papel potencial da regulação do pool de tRNA no estado persistente de Mycobacterium bovis BCG durante a tuberculose.

Figura 1: Esquema do fluxo de trabalho de preparação da biblioteca AQRNA-seq. As principais etapas descritas no fluxo de trabalho são listadas no centro do esquema e conectadas às suas respectivas ilustrações gráficas por linhas pontilhadas. A descrição detalhada de cada etapa pode ser encontrada na seção Protocolo. Clique aqui para ver uma versão maior desta figura.

Figura 2: Esquema do pipeline de análise de dados AQRNA-seq. As principais etapas descritas no pipeline são listadas no centro do esquema e conectadas às suas respectivas ilustrações gráficas por linhas pontilhadas. A descrição detalhada de cada etapa está disponível no GitHub (https://github.com/Chenrx9293/AQRNA-seq-JoVE.git). Clique aqui para ver uma versão maior desta figura.

Figura 3: Eletroforese em gel de agarose dos fragmentos de cDNA após amplificação por PCR com primers de sequenciamento. (A) Imagem do gel antes da extração e purificação do gel. As pistas 7 e 14 do lado esquerdo contêm 5 μL da escada de DNA de 50 pb, enquanto as outras pistas contêm 20 μL de cada uma das 15 amostras. A localização do tamanho dos produtos de PCR indica sua concentração mais alta dentro da faixa de 175 pb (dímeros de primer) a 300 pb (dois primers + 120 pb 5S rRNA). (B) Imagem do gel após extração e purificação do gel. Para cada amostra, o bloco de gel entre 200 pb e 400 pb foi excisado para minimizar a contaminação dos dímeros de primer na biblioteca de sequenciamento. Clique aqui para ver uma versão maior desta figura.

Figura 4: Número de leituras de sequência mapeadas com êxito para a biblioteca de sequências de referência. O eixo x mostra os nomes das amostras (por exemplo, D18-69XX) agrupadas por ponto de tempo (por exemplo, Dia de Fome 0). Para cada amostra, a contagem de leitura associada a várias categorias de assunto de destino é representada usando blocos de cores empilhados uns sobre os outros. Os números localizados no centro dos blocos de cores representam as proporções de leituras correspondentes aos respectivos assuntos de destino em uma determinada amostra. Clique aqui para ver uma versão maior desta figura.

Figura 5: Cenário quantitativo dos isoaceitadores de tRNA de Mycobacterium bovis BCG em vários pontos de tempo ao longo do curso do tempo de fome e ressuscitação. Os dados brutos de abundância de tRNA foram normalizados usando o método da mediana das razões. Aqui, cada linha descreve as abundâncias normalizadas de tRNA (eixo y) como média ± erro padrão para 3 réplicas biológicas em cada ponto de tempo. No eixo x, isoaceitadores da mesma família foram agrupados e marcados com o aminoácido correspondente. Clique aqui para ver uma versão maior desta figura.

Figura 6: Gráfico de cosseno quadrado de amostras derivadas da análise de componentes principais (PCA). A PCA foi realizada com base na abundância normalizada de tRNA. O cosseno quadrado indica a importância dos componentes principais para as amostras, e as amostras foram plotadas em relação ao cosseno quadrado dos dois primeiros componentes principais. As amostras foram rotuladas usando IDs de amostra e codificadas por cores por ponto de tempo. Clique aqui para ver uma versão maior desta figura.

Figura 7: Expressão diferencial de isoaceitadores de tRNA em diferentes pontos de tempo. As abundâncias de tRNA normalizadas foram resumidas como médias (nós de linha) ± erro padrão (barras de erro) em 3 réplicas biológicas. Devido à limitação de espaço, as condições foram abreviadas da seguinte forma: S0-S20 = dias de fome 0-20; R6 = 6º dia de ressuscitação. A análise de expressão diferencial foi realizada para cada isoaceitador de tRNA, comparando vários pontos no tempo de forma pareada usando o teste da razão de verossimilhança e o teste de Wald. Letras compactas foram empregadas para representar significância estatística, onde as abundâncias de um determinado isoaceitador de tRNA em pontos de tempo compartilhando pelo menos uma letra comum não eram significativamente diferentes umas das outras. Por exemplo, a abundância de tRNA-Lys-CTT-1-1 (no painel Lysine) foi significativamente regulada de S0 para S4 e de S4 para S10, mas não de S10 para S20. Em seguida, foi significativamente regulado de S20 para R6. Clique aqui para ver uma versão maior desta figura.
Tabela 1: Oligonucleotídeos envolvidos no fluxo de trabalho de preparação da biblioteca AQRNA-seq. O padrão interno é o RNA, enquanto todos os outros oligonucleotídeos são DNA. Os primers de PCR e os primers de sequenciamento personalizados listados são específicos para as plataformas de sequenciamento. Primers de PCR adicionais podem ser projetados com novas sequências de índice. Clique aqui para baixar esta tabela.