UNIVERSIDADE ESTADUAL PAULISTA “JULIO DE MESQUITA FILHO” FACULDADE DE CIÊNCIAS AGRÁRIAS E VETERINÁRIAS CÂMPUS DE JABOTICABAL MAPEAMENTO DE GENES DE RESISTÊNCIA DA SOJA À FERRUGEM ASIÁTICA E ANÁLISE TRANSCRICIONAL NA INTERAÇÃO PATÓGENO-HOSPEDEIRO Danielle Cristina Gregorio da Silva Orientador: Prof. Dr. Antônio Orlando Di Mauro Co-orientador: Dr. Ricardo Vilela Abdelnoor Tese apresentada à Faculdade de Ciências Agrárias e Veterinárias - Unesp, Campus de Jaboticabal, como parte das exigências para a obtenção do título de Doutor em Agronomia (Genética e Melhoramento de Plantas). JABOTICABAL – SP Fevereiro de 2008 ii Silva, Danielle Cristina Gregorio da S586m Mapeamento de genes de resistência da soja à ferrugem asiática e análise transcricional na interação patógeno-hospedeiro / Danielle Cristina Gregorio da Silva. – – Jaboticabal, 2008 xiii, 153 f. : il. ; 28 cm Tese (doutorado) - Universidade Estadual Paulista, Faculdade de Ciências Agrárias e Veterinárias, 2008 Orientador: Antonio Orlando Di Mauro Banca examinadora: Manoel Victor Franco Lemos, José Baldin Pinheiro, Francismar Correa Marcelino, Ricardo César de Paula Bibliografia 1. Glycine max. 2. Phakopsora pachyrhizi. 3. Genes de resistência. I. Título. II. Jaboticabal-Faculdade de Ciências Agrárias e Veterinárias. CDU 633.34:631.52 Ficha catalográfica elaborada pela Seção Técnica de Aquisição e Tratamento da Informação – Serviço Técnico de Biblioteca e Documentação - UNESP, Câmpus de Jaboticabal. iii DADOS CURRICULARES DO AUTOR DANIELLE CRISTINA GREGORIO DA SILVA - nascida em Toledo-PR, em 02 de outubro de 1978. Possui graduação em Ciências Biológicas pela Universidade Estadual de Londrina, em Londrina-PR (conclusão em 1999) e mestrado em Genética e Melhoramento pela mesma instituição (conclusão em 2002). É docente nas Faculdades Luís Meneghel, em Bandeirantes-PR (desde 2002), onde ministra as disciplinas ‘Genética Humana e Evolução’, para o curso de Enfermagem, e ‘Genética Molecular’ e ‘Evolução’ para o Curso de Ciências Biológicas. Foi membro da comissão executiva do colegiado do curso de Ciências Biológicas nas Faculdades Luís Meneghel (2004) e atualmente é membro da Comissão Coordenadora de Concurso Público na mesma instituição. iv AGRADECIMENTOS Agradeço aos Drs. Ricardo Vilela Abdelnoor, Naoki Yamanaka, Alexandre Nepomuceno e Antônio Orlando Di Mauro, pelo incentivo e orientação. Aos membros da banca examinadora Drs. Francismar C. Marcelino, José Baldin Pinheiro, Rinaldo C. de Paula e Manoel Victor F. Lemos, por dedicarem seu precioso tempo na melhoria deste trabalho. Aos pesquisadores da Embrapa Drs. Álvaro M. R. Almeida, Eliseu Bineck, Carlos A. A. Arias e Rodrigo L. Brogin, pelo apoio técnico-científico e incentivo. Aos co-executores deste trabalho: André, João, Lívia, Luana, Adriana, Selma, Noelle, Renata e Rodrigo. Sem seus esforços a conclusão deste não teria sido possível. A Lizandra L. Catelli e Noelle G. L. Torres, pela revisão da tese e pelo imenso apoio. A todo o pessoal do laboratório de biotecnologia da Embrapa Soja e aos amigos da FALM, pelo auxílio técnico, apoio emocional e amizade. À diretoria da FALM, chefes de departamento Sandremir e Luci, e coordenadores de curso Fabiano e Ana Cecília, pelo apoio durante minha ausência, suprindo minhas faltas na instituição durante a realização do doutorado. Aos colegas de Jaboticabal, em especial a Daniela, Thaiza, Morcelli, Fernanda, Débora e Franco, pela carinhosa recepção, amizade e auxílio. Ao pessoal da Universidade Estadual de Iowa, em especial a Martijn Van De Mortel, Steven A. Witham e Thomas J. Baum, pela recepção, orientação e apoio. Aos amigos brasileiros de Iowa, em especial a Isabel e Emerson. Vocês são parte de minha família. A Eliana Lemos e à equipe do Laboratório de Bioquímica de Microrganismos e Plantas da UNESP-Jaboticabal, pela parceria na execução do experimento de microarranjos. A EMBRAPA, UNESP, JIRCAS, CAPES, CNPQ e FINEP, pelo apoio financeiro, concessão de bolsas e disponibilização de estrutura física para a realização deste trabalho. v Ao meu noivo, Fernando. Sem sua compreensão, apoio, incentivo e amor, o caminho teria sido muito mais difícil. Enfim, a todos que de alguma forma contribuíram para a realização deste trabalho, meus sinceros agradecimentos. vi SUMÁRIO Pág. CAPÍTULO 1 - CONSIDERAÇÕES GERAIS...................................................... 01 1 INTRODUÇÃO............................................................................................. 01 2 REVISÃO DE LITERATURA....................................................................... 03 2.1 O Agronegócio da Soja no Brasil......................................................... 03 2.2 A Ferrugem Asiática da Soja............................................................... 05 2.3 Disseminação e Controle da Ferrugem Asiática da Soja............. ....... 09 2.4 Respostas Moleculares de Defesa Mediadas por Genes de Resistência......................................................................................... ....... 11 2.5 Mapeamento Genético de Genes de Resistência........................ ....... 18 2.6 Métodos de Análise da Expressão Gênica Diferencial em Situação de Doença................................................................................... 23 2.7 Estudos Moleculares da Resistência da Soja à Ferrugem Asiática....................................................................................................... 28 3 OBJETIVOS................................................................................................. 31 3.1 Objetivo Geral...................................................................................... 31 3.2 Objetivos Específicos........................................................................... 31 4 REFERÊNCIAS BIBLIOGRÁFICAS............................................................. 32 CAPÍTULO 2 - MAPEAMENTO MOLECULAR DE DOIS LOCOS QUE CONFEREM RESISTÊNCIA À FERRUGEM ASIÁTICA DA SOJA............ ....... 48 Resumo.................................................................................................... ....... 48 Introdução........................................................................................................ 48 Material e Métodos................................................................................... ....... 51 Material Vegetal......................................................................................... 51 Inoculação e Avaliação.............................................................................. 51 Extração de DNA e Estratégia de Mapeamento........................................ 52 Análise dos Dados..................................................................................... 53 Resultados....................................................................................................... 53 vii Avaliação Fenotípica.......................................................................... ....... 53 Amplificação de SSR e Análise de Ligação............................................... 54 Discussão........................................................................................................ 57 Referências..................................................................................................... 60 CAPÍTULO 3 - CARACTERIZAÇÃO DE GENES DA SOJA ATIVADOS DURANTE INTERAÇÃO COM A FERRUGEM ASIÁTICA................................. 67 Resumo........................................................................................................... 67 Introdução....................................................................................................... 67 Material e Métodos......................................................................................... 70 Inoculação e Delineamento Experimental.......................................... ....... 70 Construção das Bibliotecas Subtrativas de cDNA..................................... 72 Sequenciamento e Análise de ESTs.................................................. ....... 73 Construção dos Microarranjos de cDNA.................................................... 75 Marcação Fluorescente dos cDNAs-Alvo.................................................. 75 Desnaturação e Hibridização..................................................................... 76 Obtenção das imagens e Análise dos Dados............................................ 77 Resultados e Discussão................................................................................... 78 Caracterização das Bibliotecas Subtrativas............................................... 80 Obtenção e Caracterização do Conjunto Completo de Sequências Únicas........................................................................................................ 83 Análise de Expressão Gênica.................................................................... 90 Rotas Metabólicas Potencialmente Ativadas em Resposta à Ferrugem Asiática...................................................................................... 91 Proteínas com Estrutura Análoga a Genes de Resistência....................... 96 Sinalização das Respostas de Defesa Mediada por Proteínas................. 96 Sinalização Hormonal................................................................................ 97 Mudanças na Expressão Gênica Controladas ao Nível Transcricional e Traducional...................................................................... 98 Produção de Espécies Reativas de Oxigênio............................................ 99 Indução de Morte Celular........................................................................... 100 viii Degradação Protéica Associada à Via Ubiquitina/Proteossomo............... 100 Chaperonas................................................................................................ 102 Produção de Fitoalexinas........................................................................... 102 Proteínas Antimicrobianas......................................................................... 103 Reforço na Parede Celular......................................................................... 104 Alterações no Metabolismo de Aminoácidos............................................. 104 Alterações no Transporte de Substâncias................................................. 105 Senescência............................................................................................... 107 Outras Rotas Metabólicas.......................................................................... 107 Estratégia de Defesa Mediada pelo Rpp2................................................. 108 Referências Bibliográficas.............................................................................. 111 CAPÍTULO 4 - IMPLICAÇÕES........................................................................... 125 APÊNDICE.......................................................................................................... 129 ix LISTA DE FIGURAS Pág. CAPÍTULO 1- CONSIDERAÇÕES GERAIS....................................................... 01 Figura 1: Desenvolvimento de fungos biotróficos, como oídios e ferrugens............................................................................................................. 06 Figura 2: Sintomas de plantas suscetíveis e resistentes à ferrugem asiática da soja.................................................................................................. 08 Figura 3: Visão geral das vias de transdução de sinal que ativam e coordenam as respostas de defesa das plantas......................................... ....... 18 Figura 4: Diagrama esquemático da SSH................................................... ...... 27 CAPÍTULO 2 - MAPEAMENTO MOLECULAR DE DOIS LOCOS QUE CONFEREM RESISTÊNCIA À FERRUGEM ASIÁTICA DA SOJA.................... 48 Figura 1: Mapas de ligação das regiões genômicas do Rpp2 (A) e Rpp4 (B).............................................................................................................. 56 CAPÍTULO 3 - CARACTERIZAÇÃO DE GENES DA SOJA ATIVADOS DURANTE INTERAÇÃO COM A FERRUGEM ASIÁTICA................................. 67 Figura 1: Estratégia de subtração utilizada para a montagem das bibliotecas de cDNA............................................................................................ 73 Figura 2: Lesões de ferrugem asiática da soja observadas através de microscópio estereoscópico na face abaxial de folha de soja aos 12 dias após a inoculação........................................................................... ....... 80 Figura 3: Espécies mais representadas entre os similares encontrados no GenBank para as sequências das quatro bibliotecas e os respectivos números de similares encontrados..................................................................... 82 Figura 4: Padrão de co-expressão das sequências entre as quatro bibliotecas E024, E192, P024 e P192......................................................... ....... 84 Figura 5: Classificação das supostas funções gênicas das ESTs das bibliotecas E024, E192, P024 e P192......................................................... ....... 89 Figura 6: Potencial estratégia de defesa mediada pelo gene Rpp2................... 110 x LISTA DE TABELAS Pág. CAPÍTULO 2 - MAPEAMENTO MOLECULAR DE DOIS LOCOS QUE CONFEREM RESISTÊNCIA À FERRUGEM ASIÁTICA DA SOJA............ ....... 48 Tabela 1: Análise de qui-quadrado da segregação fenotípica quanto à resistência à ferrugem asiática da soja (geração F2:3) e marcadores SSR (geração F2) nas populações de mapeamento do Rpp2 e do Rpp4................... 55 CAPÍTULO 3 - CARACTERIZAÇÃO DE GENES DA SOJA ATIVADOS DURANTE INTERAÇÃO COM A FERRUGEM ASIÁTICA................................. 67 Tabela 1: Descrição de cada uma das quatro bibliotecas subtrativas e do conjunto completo de sequências quanto ao número, qualidade, agrupamento e redundância............................................................................... 82 Tabela 2: Funções moleculares expressas exclusivamente em cada uma das quatro bibliotecas e nas combinações de bibliotecas.................................. 86 Tabela 3: ESTs da soja inéditas identificadas neste trabalho............................ 92 Tabela 4: Seqüências diferencialmente expressas segundo análsie de microarranjos, para as quais foi encontrado um similar significativo (valor-e máximo 1E-8).......................................................................................... 93 xi LISTA DE ABREVIATURAS AFLP - “Amplified Fragment Length Polymorphism” ou Polimorfismo de Comprimento de Fragmento Amplificado Avr - avirulência AVRDC - “Asian Vegetable Research and Development Center” ou Centro Asiático de Pesquisa e Desenvolvimento Vegetal BAC - “Bacterial Artificial Chromosome” ou Cromossomo Artificial Bacteriano BSR - “Brown Stem Rot” ou Podridão Marrom do Caule BSA - “Bulked Segregant Analysis” ou Análise de Agrupamentos Segregantes CC - “Coiled-Coil” ou Bobina Enrolada CONAB - Companhia Nacional de Abastecimento E024 - amostra de folha do genótipo Embrapa 48 às 24 horas após inoculação E192 - amostra de folha do genótipo Embrapa 48 às 192 horas após inoculação EMBRAPA - Empresa Brasileira de Pesquisa Agropecuária EST - “Expressed Sequence Tags” ou Etiquetas de Sequências Expressas HR - “Hipersensitive Response” ou Resposta de Hipersensibilidade ITS - “Internal Transcribed Sequences” ou Seqüências Internas Transcritas JA - “Jasmonic Acid” ou Ácido Jasmônico LOX - LipOXigenases LRR - “Leucine Rich Repeats” ou Repetições Ricas em Leucina MAS - “Marker Assisted Selection” ou Seleção Assistida por Marcadores MPSS - “Massively Pararel Signature Sequencing” ou Sequenciamento de Assinaturas Massivamente Paralelo NB - “Nucleotide Binding” ou Ligação a Nucletídio NCBI - “National Center for Biotechnology Information” ou Centro Nacional (Americano) de Informação Biotecnológica NO - “Nitric Oxide” ou Óxido Nítrico ORF - “Open Reading Frame” ou Matriz Aberta de Leitura xii P024 - amostra de folha do genótipo PI230970 às 24 horas após inoculação P192 - amostra de folha do genótipo PI230970 às 192 horas após inoculação PAL - “Phenilalanine Amonialyase” ou Fenilalanina Amonialiase PCR - “Polymerase Chain Reaction” ou Reação em Cadeia da Polimerase PR - “Pathogenesis Related” ou Relacionada à Patogênese QTL - “Quantitative Trait Locus” ou Loco de Característica Quantitativa R - “Resistance” ou Resistência RAPD - “Random Amplified Polymorphic DNA” ou DNA Polimórfico Amplificado ao Acaso RB - “Reddish-Brown” ou Marrom-Avermelhado RFLP - “Restriction Fragment Length Polymorphism” ou Polimorfismo de Comprimento de Fragmento de Restrição RGAs - “Resistance Gene Analogs” ou Análogos de Genes de Resistência RILs - “Recombinant Imbred Lines” ou Linhagens Endogâmicas Recombinantes ROS - “Reactive Oxygen Species” ou Espécies Reativas de Oxigênio RPA - “Ribonuclease Protection Assay” ou Ensaio de Proteção de Ribonuclease SA - “Salicilic Acid” ou Ácido Salicílico SAGE - “Serial Analysis of Gene Expression” ou Análise Serial da Expressão Gênica SAR - “Systemic Aquired Resistance” ou Resistência Sistêmica Adquirida SCN - “Soybean Cist Nematode” ou Nematóide do Cisto da Soja SDS - Sódio Duodecil Sulfato Ser/Thr - “Serine/Threonine” ou Serina/Treonina SOD - Superóxido Dismutase SSC - Cloreto de Sódio/Sulfato de Sódio SSH - “Supressive Subtractive Hibridization” ou Hibridização Subtrativa Supressiva SSR - “Simple Sequence Repeat” ou Sequência Simples Repetida TIR- “Toll-Interleukin Receptor” ou Receptor Toll-Interleucina USDA - “United States Department of Agriculture” ou Departamento de Agricultura dos Estados Unidos xiii MAPEAMENTO DE GENES DE RESISTÊNCIA À FERRUGEM ASIÁTICA DA SOJA E ANÁLISE TRANSCRICIONAL NA INTERAÇÃO PATÓGENO-HOSPEDEIRO RESUMO - A ferrugem asiática da soja é causada pelo fungo Phakopsora pachyrhizi. Quatro genes de resistência à ferrugem da soja distintos, Rpp1 a 4, foram descritos. O objetivo deste trabalho foi mapear geneticamente os genes de resistência Rpp2 e Rpp4 e identificar genes induzidos pela resposta de defesa mediada pelo gene Rpp2. Duas populações F2:3 derivadas dos cruzamentos entre as linhagens resistentes PI 230970 (Rpp2) e PI 459025 (Rpp4) com a cultivar suscetível BRS 184, e marcadores SSR, foram utilizados neste estudo. Os locos Rpp2 e Rpp4 foram mapeados nos grupos de ligação J e G da soja, respectivamente, e os marcadores associados terão grande valor no processo de seleção assistida por marcadores para este caráter. Quatro bibliotecas de cDNA, derivadas de folhas do genótipo resistente PI 230970 (Rpp2) e do genótipo suscetível Embrapa 48, obtidas às 24 e 192 horas após a inoculação com esporos de P. pachyrhizi, foram geradas utilizando a metodologia SSH e tiveram seu padrão de expressão avaliado por análise de microarranjos de cDNA. Foram produzidas 3.807 sequências viáveis, das quais 670 foram sequências únicas, sendo 149 destas (22,2%) novas sequências da soja. A categoria funcional mais representativa para as quatro bibliotecas foi proteção celular, defesa e virulência. Apenas 65 trnascritos foram diferencialmente expressos, estando eles envolvidos na geração de espécies reativas de oxigênio, fitoalexinas e proteínas antimicrobianas, morte celular e senescência, modificação, estabilização e degradação protéica, controle da expressão gênica e reforço de parede celular. Este trabalho permitiu a identificação de novas seqüências da soja e contribuiu para a elucidação da estratégia de defesa mediada pelo gene Rpp2. Palavras-Chave: Glycine max, Phakopsora pachyrhizi, genes de resistência, transcritoma. xiv MAPPING OF GENES FOR RESISTANCE TO SOYBEAN RUST AND TRANSCRIPTIONAL ANALYSIS ON THE PLANT-PATHOGEN INTERACTION SUMMARY - Asian soybean rust is caused by the fungus Phakopsora pachyrhizi. Four resistance genes to this disease, Rpp1 to 4, were previously identified. The aim of this work was to identify SSR markers linked to Rpp2 and Rpp4, and to characterize genes up-regulated under soybean-Asian rust interaction. Two F2:3 populations derived from the crosses PI 230970 (Rpp2) x BRS 184 and PI 459025 (Rpp4) x BRS 184, and SSR markers were used for mapping. Rpp2 and Rpp4 loci were mapped on the soybean linkage groups J and G, respectively, and the associated markers will be highly valuable to assist the introduction of these genes into soybean cultivars. Four cDNA libraries derived from leaf samples of the resistant genotype PI 230970 and the susceptible genotype Embrapa 48, obtained at 24 and 192 hours after inoculation, were generated using SSH and evaluated by cDNA microarray analysis. A total of 3,807 reliable ESTs were obtained, 670 from them were unique. The most representative functional category in all libraries was “cell rescue, defense and virulence”. Only 65 transcripts were differentially expressed among treatments, and were involved in the generation of reactive oxygen species, phytoalexins and antimicrobial proteins; cell death and senescence; protein fate; gene expression control; and reinforcement of cell wall. This work produced novel soybean sequences and contributed to the comprehension of the soybean mechanism of resistance to Asian rust mediated by Rpp2. Keywords: Glycine max, Phakopsora pachyrhizi, resistance genes, transcriptome. CAPÍTULO 1 - CONSIDERAÇÕES GERAIS 1 INTRODUÇÃO O Brasil é o segundo maior produtor mundial de soja [Glycine max (L.) Merrill], produzindo 58,4 milhões de toneladas do grão na safra de 2006/07. A soja é a cultura que recebe maior investimento em área plantada no país e é um dos principais produtos nacionais de exportação (CONAB, 2007). No entanto, fatores bióticos e abióticos que limitam a produtividade dessa cultura colocam em risco toda essa riqueza. Atualmente, a principal doença que afeta a produção de soja no Brasil é a ferrugem asiática, causada pelo fungo Phakopsora pachyrhizi Sydow & Sydow. Para a ferrugem asiática da soja, foram descritos na literatura quatro genes de resistência dominantes e independentes, denominados Rpp1 a 4 (CHENG & CHAN, 1968; HIDAYAT & SOMAATMADJA, 1977; BROMFIELD & HARTWIG, 1980; HARTWIG, 1986). Quando a doença foi primeiramente detectada no Brasil, todos estes genes eram efetivos. Entretanto, em 2003, uma nova raça fisiológica de P. pachyrhizi quebrou a resistência conferida pelos genes Rpp1 e Rpp3. Somente o Rpp2 e o Rpp4 permanecem resistentes à raça fisiológica do fungo presente atualmente no Brasil. Os métodos tradicionais utilizados para monitorar a transferência de genes a cultivares são demorados e trabalhosos. Nos últimos anos, a seleção facilitada por marcadores moleculares tem recebido muita atenção como método rápido e viável para o melhoramento da resistência a patógenos nas plantas cultivadas (FRANCIA et al., 2005). A utilização de perfis moleculares que indiquem a presença de um dado genótipo em um indivíduo eliminam o efeito ambiental no processo de seleção, aumentando a eficiência do programa de melhoramento. Ainda, quando diferentes locos codificam para o mesmo fenótipo, a utilização de marcadores moleculares possibilita a fácil identificação dos indivíduos portadores de mais do que um loco para o caráter, facilitando assim a piramidação destes em cultivares. Além de facilitar a piramidação 2 dos genes de resistência, a utilização de marcadores moleculares pode ainda permitir a localização genômica dos locos de resistência, o que cria a possibilidade de clonagem dos mesmos através da estratégia de clonagem posicional. Recentemente, foram identificados marcadores SSR (“Simple Sequence Repeats” ou Sequências Simples Repetidas) ligados a locos de resistência à ferrugem asiática da soja, como exemplo o Rpp1, presente na PI200492, e o Rpp presente na cultivar japonesa Hyuuga, para o qual não foi determinada a relação com os genes Rpp1 a 4 (HYTEN et al., 2007; MONTEROS et al., 2007). Entretanto, os locos ainda efetivos no Brasil Rpp2 e Rpp4 não foram mapeados até o momento. Ainda, pouco se sabe quanto ao mecanismo molecular de interação entre a soja e a ferrugem asiática e as complexas vias metabólicas que são ativadas em genótipos suscetíveis e resistentes em função desta interação, culminando na efetiva infecção da planta pelo patógeno ou no impedimento da proliferação do mesmo. Tendo em vista o grande potencial de perdas na produção de soja que a ferrugem asiática representa para o país e considerando a necessidade de desenvolvimento de cultivares resistentes e da falta de conhecimentos sobre a interação entre soja e ferrugem, tornam-se prioritários o estudo e a pesquisa sobre a resistência a esta doença. A caracterização molecular dos genes de resistência e dos genes que com eles interagem poderá auxiliar o desenvolvimento de métodos eficientes de controle do fungo, além de promover um avanço na compreensão da reação da soja a outros patógenos. 3 2 REVISÃO DE LITERATURA 2.1 O Agronegócio da Soja no Brasil Apesar de conhecida e explorada no Oriente há mais de 5.000 anos, a soja [Glycine max (L.) Merrill] teve seu cultivo ignorado pelo Ocidente até a segunda década do século XX, quando os Estados Unidos da América (EUA) iniciaram sua exploração comercial (primeiro como forrageira e, posteriormente, como grão). Em 1940, no auge do seu cultivo como forrageira, foram cultivados, naquele país, cerca de 2 milhões de ha com tal propósito. A partir de 1941, a área cultivada para grãos superou a cultivada para forragem, cujo cultivo declinou rapidamente, até desaparecer em meados dos anos 60, enquanto a área cultivada para a produção de grãos crescia de forma exponencial, não apenas nos EUA, como também no Brasil e na Argentina, principalmente (BERTRAND et al., 1987; MORAIS & SILVA, 1996; MULLER, 1981). A soja é um grão que dá origem a diversos produtos e subprodutos muito utilizados pela agroindústria, indústria química e de alimentos. Na alimentação humana, seu uso mais conhecido é como óleo refinado, mas também entra na composição de vários produtos embutidos, como salsichas, chocolates, temperos para saladas, maioneses, entre outros produtos. Sua proteína é a base de ingredientes como massas, produtos de carne, cereais, misturas preparadas, bebidas, alimentação para bebês e alimentos dietéticos. Também é muito usada pela indústria de adesivos e nutrientes, alimentação animal, adubos, formulador de espumas, fabricação de fibra, revestimento, papel e emulsão de água para tintas. Recentemente, sua utilização vem crescendo também como fonte alternativa de combustível. O biodiesel de soja já vem sendo testado por instituições de pesquisa brasileiras (EMBRAPA, 2007c). A produtividade da soja brasileira, ao longo dos anos, evoluiu muito devido ao avanço das tecnologias de produção que foram adaptadas às condições tropicais e subtropicais do Brasil (HASSE & BUENO, 1996). Quando comparada com as demais 4 culturas brasileiras, a soja foi a que mais cresceu em área cultivada. Segundo a Companhia Nacional de Abastecimento - CONAB, houve um investimento em área plantada com a soja no Brasil de 20,7 milhões ha (44,8% da área cultivada total) na safra de 2006/07. A produção de soja foi de 58,4 milhões de toneladas (44,4% do total de grãos) na safra de 2006/07. A região Centro-Oeste foi a mais expressiva em área plantada com soja, representando 44,0% da área total plantada com soja no Brasil, em produtividade, atingindo patamares de 2.910 kg/ha, e produção, com 45,4% da produção brasileira dessa cultura. Entretanto, os mais altos valores de produtividade por Estado da Federação em 2006/07 foram alcançados por Roraima (o único Estado da Federação livre da ferrugem asiática), com 3.070 kg/ha (CONAB, 2007). A soja sustenta o agronegócio brasileiro, sendo o maior produto nacional de exportação, com 31,9 milhões de toneladas de janeiro a setembro de 2007, correspondendo a 8,97 bilhões de dólares ou 21,8% do total de exportações de produtos agropecuários, segundo a CONAB. Neste mesmo período, o agronegócio trouxe um saldo positivo à balança comercial brasileira de 35,12 bilhões de dólares, sendo que o saldo positivo da balança comercial total foi de 30,93 bilhões de doláres. Em outras palavras, o saldo positivo gerado pelo agronegócio compensou o saldo negativo gerado pelos outros setores da economia; a participação do agronegócio na balança comercial total foi de 113,56%. A China tem sido nos últimos anos o maior importador da soja brasileira em grão, seguida dos países baixos. O Brasil também exporta o farelo e o óleo de soja, porém em quantidades menores que o grão (CONAB, 2007). Segundo estimativas da safra 2007/08, feitas pelo Departamento de Agricultura dos Estados Unidos em novembro de 2007 (USDA, 2007), o Brasil continua sendo o segundo maior produtor mundial de soja com 62 milhões de toneladas, atrás apenas dos EUA, com 70,6 milhões de toneladas. Na terceira posição mundial encontra-se a Argentina, com 47 milhões de toneladas e em quarto a China, com 14,3 milhões de toneladas. Na safra anterior, de 2006/07, os EUA produziram 86,77 milhões, o Brasil 59 milhões, a Argentina 47,2 milhões e a China 16,2 milhões de toneladas de soja. 5 2.2 A Ferrugem Asiática da Soja A soja pode ser infectada por duas espécies de fungo que causam ferrugem: o Phakopsora meibomiae (Arthur), nativo do continente americano, e o Phakopsora pachyrhizi Sydow e Sydow, nativo do oriente, também denominado ferrugem asiática. A ferrugem americana ocorre em condições de temperaturas médias abaixo de 25°C e umidade relativa elevada, raramente causa perdas significativas e ocorre nas regiões dos Cerrados e Sul do País. Por outro lado, o fungo da ferrugem asiática é capaz de se desenvolver em temperaturas variando de 15 °C a 30 °C, pode se adaptar a qualquer região brasileira e pode causar severas perdas na cultura da soja (EMBRAPA, 2004). O fungo P. pachyrhizi, um basidiomiceto da ordem Uredinales, foi descrito pela primeira vez no Japão, em 1902, e em 1914 já havia se espalhado por diversos países do sudeste da Ásia (HENNING & GODOY, 2006). Nas Américas, a doença chegou em 2001, mais especificamente no Brasil, Paraguai e Bolívia. No hemisfério norte, foi constatada somente em 2004, na Colômbia e nos Estados Unidos (SCHNEIDER et al., 2005). As ferrugens são parasitas obrigatórios (biotróficos), isto é, não matam seus hospedeiros, mas desenvolvem estratégias efetivas para explorar células vivas como fontes de alimento. O ciclo de vida de P. pachyrhizi inicia-se com uredósporos, oriundos de urédias, que atingem as folhas de soja na face superior ou inferior e germinam, se a temperatura for favorável e houver pelo menos seis horas de molhamento foliar. Dos uredósporos germinados forma-se um tubo germinativo, capaz de penetrar os tecidos foliares. Em uma infecção por P. pachyrhizi, a penetração ocorre diretamente pela epiderme intacta, mas pode também ocorrer pelos estômatos (ZAMBOLIN, 2006). 6 Figura 1: Desenvolvimento de fungos biotróficos, como oídios e ferrugens. A e C: Os uredósporos germinam formando um tubo germinativo que desenvolve um apressório; B e C: a hifa de infecção entra em contato com células do parênquima, se diferenciando em célula mãe do haustório, a qual produz o haustório dentro da célula vegetal. A partir deste ponto o patógeno começa a absorver nutrientes da planta. Adaptado de MENTABERRY (2007). As fotos são de Peronospora parasitica desenvolvendo-se em folha de Arabidopsis thaliana. Nas ferrugens em geral, o tubo germinativo forma um apressório, uma estrutura que facilita a penetração da hifa. Dentro da folha, a hifa de infecção entra em contato 7 com uma célula do mesofilo e é formada uma célula mãe do haustório, a qual realiza a penetração da célula hospedeira e a diferenciação do haustório (Figura 1). Este é uma estrutura altamente ramificada, delimitada pela membrana celular da célula hospedeira, denominada de membrana extra-haustorial, uma matriz extra-haustorial e a parede celular do fungo. O haustório proporciona uma ampla superfície de contato com a célula hospedeira e através dele são adquiridas quantidades massivas de açúcares e aminoácidos por mecanismos de simporte a favor de gradiente de prótons. Propõe-se que a formação desta estrutura dependa de compostos voláteis, como o nonanal, decanal e acetato de hexenil, emitidos pelo hospedeiro (PANSTRUGA, 2003; MENDGEN et al., 2006). As lesões produzidas por P. pachyrhizi são angulares, podendo apresentar coloração bronze (“tan”), indicando suscetibilidade, ou castanho-avermelhada (“reddish- brown” ou RB), caracterizando resistência (Figura 2) (BROMFIELD, 1984). Os primeiros sintomas são minúsculos pontos de coloração clara, e com o avanço da doença são formadas as urédias com minúsculos poros de onde são liberados os esporos. Nos genótipos suscetíveis, ocorre a senescência precoce, resultando em adiantamento do ciclo e prejuízo no enchimento dos grãos. O grau de esporulação nas lesões de genótipos suscetíveis é elevado, e em genótipos resistentes, baixo, quando existente. O período de latência (dias entre a inoculação e o aparecimento das lesões) e a viabilidade de esporos variam entre os genótipos (SOARES, 2007). As lesões da ferrugem asiática podem surgir em 5 a 7 dias e os esporos em 9 a 12 dias. A infecção ocorre em umidade relativa superior a 60%. Chuvas constantes de baixa intensidade e temperaturas entre 18 e 25 ºC por vários dias favorecem a epidemia da doença (ZAMBOLIN, 2006). As chuvas sempre possibilitam condições ideais para a infecção com ferrugem, seja por liberar os uredósporos, seja por fornecer água para a germinação, ou seja por manter a temperatura para a infecção próxima da ideal (BERGAMIN FILHO, 2006). O fungo da ferrugem asiática da soja apresenta uma alta diversidade genética. Vinte raças de ferrugem asiática foram identificadas em amostras coletadas em plantas de soja e hospedeiros selvagens no Japão (YAMAOKA et al., 2002), nove raças em 42 8 isolados em Taiwan (AVRDC, 1985) e 59 raças em 69 isolados na Tailândia (POONPOLGUL, 2004). No Brasil, embora não haja comprovação científica, devido à inexistência de hospedeiros diferenciadores de raças até o momento, admite-se que haja também inúmeras raças fisiológicas do patógeno em campos de soja (ZAMBOLIN 2006). A raça predominante atualmente nos campos de soja do norte do Paraná é muito similar à raça proveniente do Zimbábue, conforme atestam resultados de seqüenciamento de regiões ITS (“Internal Transcribed Sequences”) do genoma do fungo (Ricardo Vilela Abdelnoor, comunicação pessoal). Figura 2: Sintomas de plantas suscetíveis e resistentes à ferrugem asiática da soja. A: face inferior de folha de planta suscetível. As lesões são de cor bronzeada (“tan”), ocorre senescência precocemente e há alta concentração de lesões com urédias esporulando. B: face inferior de folha de planta resistente. As lesões são de cor marrom-avermelhada (RB), a senescência é reduzida e as lesões possuem poucas ou nenhuma urédia com baixa esporulação. Foto: Danielle C. G. Silva, 2006. 9 A ferrugem asiática possui uma ampla gama de hospedeiros, que inclui além de Glycine max, outras espécies do gênero Glycine e outros gêneros de leguminosas. Uma coleção de 95 espécies de plantas hospedeiras, distribuídas em 42 gêneros da família Fabaceae (CALDWELL & MACLAREN, 2004) são hospedeiras da ferrugem. Assim, apesar do fungo ser biotrófico, o inóculo pode ser mantido nas épocas de entressafra, parasitando, além de soja voluntária, outras espécies de planta presentes no mesmo ambiente, como por exemplo, o feijão de lima (Phaseolus vulgaris var. macrocarpus), o desmódio (Desmodium sp.), o kudzu (Pueraria lobata) e a soja perene (Neonotonia wightii). Ainda, pode ser disseminado pelo vento a partir de cultivos de soja a quilômetros de distância. 2.3 Disseminação e Controle da Ferrugem Asiática da Soja Até o ano de 1999, a ferrugem asiática era limitada aos países orientais produtores de soja. A doença foi primeiramente detectada nas Américas em 2001, no sul do Brasil, Paraguai e Bolívia, e desde então tem se espalhado por todos os estados brasileiros produtores de soja, com exceção do estado de Roraima (EMBRAPA, 2007a). A ferrugem da soja é uma doença extremamente importante para as regiões tropicais e subtropicais, fato que a coloca, no quesito danos à produção, próxima aos modelos mais destrutivos de doença tropical (BERGAMIN FILHO, 2006). A situação da ferrugem está mais preocupante no Brasil devido ao fato da soja ser plantada como uma monocultura, as áreas de cultivo serem extensas e contínuas e o clima ser favorável ao desenvolvimento da doença, principalmente no Brasil central (HENNING & GODOY, 2006). A cada ano, a severidade da doença tem aumentado nas regiões de clima mais favorável e, consequentemente, também tem aumentado o custo do controle (YORINORI, 2006). As perdas em grãos provocadas pela ferrugem asiática da soja totalizaram aproximadamente 4,5% da produção brasileira de soja em 2006/07, o que equivale a 10 2,67 milhões de toneladas de grãos ou US$ 615,7 milhões, considerando o preço médio de US$ 230,6 por tonelada de soja. Somando-se o custo da operação de controle, cuja média nacional ficou em 2,3 aplicações por ha, representando US$1,58 bilhão, o custo total da ferrugem asiática no Brasil na safra 2006/07 é de US$ 2,19 bilhões (EMBRAPA, 2007ab). Desde 2001 a 2006, a doença causou prejuízos acumulados (incluindo perdas de grãos, custos com o controle químico e perdas de arrecadação) estimados em US$ 5,1 bilhões (YORINORI, 2006). Assim, o custo total com a doença de 2001 a 2007, foi de cerca de US$ 7,3 bilhões. Desde a detecção da doença no Brasil, o principal método de controle utilizado é a aplicação de fungicidas, os quais oneram substancialmente a produção, além de serem nocivos ao meio ambiente. Na safra de 2005/06, além da estiagem que afetou extensas áreas de soja, o controle da ferrugem nas regiões onde ocorreram chuvas foi dificultada por diversos fatores: ocorrência da doença desde o início do desenvolvimento da soja; dificuldade na identificação inicial; conceito arraigado da aplicação preventiva ou protetora nos estádios R1 a R3 (FEHR & CAVINNES, 1977), resultando em atraso no controle inicial; ineficiência na cobertura das folhas com fungicidas; falhas na tecnologia de aplicação e, finalmente, a dificuldade financeira dos produtores por falta de crédito para compra dos insumos. Nos cerrados brasileiros, diversas propriedades ficaram sem cultivo de soja, mesmo com arrendamentos gratuitos; alguns produtores trocaram herbicidas e inseticidas por fungicidas e, outros, semearam a soja sem condições de pulverizar, o que aumentou o potencial de inóculo na região (YORINORI, 2006). Um dos métodos mais promissores de controle é o uso de resistência e/ou tolerância genética. Existem quatro genes R dominantes e independentes que conferem resistência à ferrugem asiática descritos na literatura: Rpp1, encontrado na PI 200492 (CHENG & CHAN, 1968), Rpp2, encontrado na PI 230970 (HIDAYAT & SOMAATMADJA, 1977), Rpp3, encontrado na PI 462312 (BROMFIELD & HARTWIG, 1980) e Rpp4, encontrado na PI 459025 (HARTWIG, 1986). Além destes genes, em vários outros genótipos foram encontrados locos de resistência para os quais não se sabe a relação com os genes previamente descritos. Uma avaliação de mais de 1000 11 acessos do banco de germoplasma da Embrapa Soja, Londrina-PR, quanto à resistência à ferrugem asiática, revelou 64 linhagens avançadas e cultivares brasileiras com níveis variáveis de resistência/tolerância (YORINORI & NUNES JUNIOR, 2007). Em outro trabalho, uma avaliação de 263 acessos japoneses e chineses, resultou num total de nove variedades resistentes, tolerantes ou parcialmente tolerantes (PASSIANOTTO et al., 2007). Em outra análise, mais ampla, 16.595 acessos do banco de germoplasma do USDA, na Universidade de Illinois, foram avaliados quanto à resistência à ferrugem (MILES et al., 2006). Este trabalho resultou em 805 genótipos promissores como fontes de resistência à ferrugem asiática. No momento, não existem cultivares resistentes/tolerantes à ferrugem asiática disponíveis para a produção de grãos no Brasil. Entretanto, a Fundação MT e a Embrapa divulgaram a disponibilização de variedades resistentes para as safras 2009/2010 e 2010/2011, respectivamente (NÓRCIO, 2007; AGÊNCIA DA NOTÍCIA, 2008). 2.4 Respostas Moleculares de Defesa Mediadas por Genes de Resistência A resistência de um hospedeiro, dentro do contexto da fisiologia do parasitismo, pode ser definida como a capacidade da planta em atrasar ou evitar a entrada e/ou a subseqüente atividade de um patógeno em seus tecidos. Embora as plantas, na natureza, estejam normalmente expostas a um número incalculável de patógenos, a resistência das mesmas mostra-se como regra, enquanto que a suscetibilidade mostra- se como exceção (PASCHOLATI & LEITE, 1995). O patógeno falha em se estabelecer em uma planta devido a uma entre três possíveis razões: a planta é incapaz de suportar os requerimentos fisiológicos de um patógeno potencial e, portanto, não pode ser considerada um hospedeiro; a planta possui barreiras estruturais pré-formadas ou compostos tóxicos que impedem o desenvolvimento do patógeno, novamente sendo considerada uma planta não- hospedeira; e a planta reconhece o ataque do patógeno, quando são ativados 12 mecanismos de defesa que impedem a proliferação do mesmo. Estes três tipos de interação são ditas incompatíveis, mas somente a última depende da indução de resposta na planta (HAMMOND-KOSACK & JONES, 1996). A resposta de defesa induzida por genes de resistência raça-específicos, particularmente aquelas envolvendo parasitas biotróficos, é frequentemente explicada pela teoria gene-a-gene. O reconhecimento de um patógeno resulta da interação direta do produto de um gene R (de resistência), geralmente dominante, com o produto do gene Avr (de avirulência) correspondente do patógeno (FLOR, 1971). O reconhecimento da proteína Avr leva a uma cascata de eventos de sinalização que irá ativar genes da célula vegetal hospedeira, resultando em uma rápida resposta de defesa. A maioria dos genes R codificam para cinco classes principais de proteínas de resistência: proteínas com domínio transmembrana e Repetições Ricas em Leucina (“Leucine Rich Repeats” ou LRR) localizadas no lado extracelular; proteínas com sinal transmembrana ligado a um domínio citoplasmático ‘bobina enrolada’ (“coiled-coil” ou CC); proteínas com domínio citoplasmático serina/treonina quinase (“Serine/Threonine” ou Ser/Thr) provavelmente associadas a membrana; proteínas com domínio citoplasmático Ser/Thr quinase, domínio transmembrana e LRR; e proteínas citoplasmáticas com domínio de ligação a nucleotídio (“Nucleotide Binding” ou NB) associado a LRR (NB-LRR) (DANGL & JONES, 2001). A classe NB-LRR pode ser subdividida com base na composição da região amino-terminal dos genes. Muitos têm um domínio com homologia aos domínios intracelulares de sinalização “Toll” de drosófila e receptores de interleucina “(IL)-1” de mamíferos, sendo chamados de TIR (“Toll-Interleukin Receptor”). Os genes com esta característica amino-terminal são denominados TIR-NB-LRR. Outros genes contêm domínios CC na região amino- terminal, sendo chamados de CC-NB-LRR (DANGL e JONES 2001). Apesar destas serem as classes mais importantes, foram encontrados genes de resistência apresentando outras estruturas, como o Hm1 de milho, cujo produto é uma toxina HC redutase (JOHAL & BRIGGS, 1992), e o RGA1 de arroz, cujo produto possui domínio de proteína G (ASSMANN, 2005), e outros. 13 O reconhecimento das proteínas Avr pelas proteínas R em interações entre plantas e patógenos biotróficos possivelmente ocorre no citoplasma. As proteínas AvrL567 da ferrugem do linho (Melampsora lini), as quais são expressas em haustórios e contêm peptídios sinais, sugerindo sua secreção na matriz extra-haustorial, são reconhecidas por proteínas R citoplasmáticas do tipo NB-LRR. Este fato sugere a existência de um mecanismo específico de translocação que envie estas proteínas através da membrana extra-haustorial para dentro do citoplasma da célula hospedeira (CATANZARITI et al., 2006). Além deste dado, uma proteína de função desconhecida secretada pelo haustório da ferrugem da fava Uromyces fabae foi imunolocalizada no citoplasma e no núcleo da célula hospedeira, comprovando a possibilidade de que proteínas Avr expressas nos haustórios sejam transportadas para dentro das células vegetais hospedeiras (KEMEN et al., 2005). Apesar destas e de outras evidências, o modo pelo qual esse processo ocorre ainda é desconhecido (ELLIS et al., 2006). Entre as principais respostas de defesa desencadeadas por genes R estão: liberação de radicais livres, resposta de hipersensibilidade, reforço da parede celular por lignificação e produção de calose, produção de compostos antimicrobianos (fitoalexinas e enzimas hidrolíticas) e indução de proteínas relacionadas à patogênese (“Pathogenesis Related” ou PR) (STASKAWICZ et al., 1995; KNOGGE, 1996). O principal componente da resposta de defesa ativada pela interação entre o produto do gene R e o produto do gene Avr correspondente é a Resposta de Hipersensibilidade (“Hipersensitive Response” ou HR). A HR é definida com a morte das células do hospedeiro dentro de poucas horas após o contato com o patógeno. Em interações com patógenos biotróficos, os quais desenvolvem íntima associação com as células do hospedeiro através dos haustórios, a morte celular priva o patógeno do acesso a nutrientes, controlando seu desenvolvimento. A HR é uma resposta de defesa local, ativada somente nos pontos de contato entre a planta e o patógeno. As células mortas neste processo contém altas concentrações de moléculas antimicrobianas e, por isso, não são subsequentemente atacadas por patógenos necrotróficos. A morte celular na HR parece resultar de dois processos principais: a produção de espécies reativas de 14 oxigênio, as quais são tóxicas para a célula e levam à morte celular, e a ativação de um processo de morte celular programada (HAMMOND-KOSACK & JONES, 1996). As espécies reativas de oxigênio (“Reactive Oxygen Species” ou ROS) correspondem a uma das primeiras respostas ativadas em muitas interações incompatíveis. Existem inúmeras fontes potenciais de ROS: as fontes protoplásticas, representadas pelas mitocondrias, cloroplastos e peroxissomos; as NADPH oxidases e as superóxido sintases NADH-dependentes da membrana celular; e fontes apoplásticas de geração de peróxido de hidrogênio (H2O2), incluindo peroxidases, amino oxidases e oxalato sintases. Em células de feijão, o peróxido de hidrogênio parece ser originado de uma peroxidase ligada a parede celular em um sistema de três componentes, requerendo fluxo de íons e levando à alcalinização extracelular e à liberação de moléculas redutoras, como os tióis (BOWEL et al., 2002). Os ânions superóxido (O2 -), gerados na membrana ou nos protoplastos e secretados ao ambiente extracelular, seriam rapidamente dismutados de um modo não enzimático ou através da ação da superóxido dismutase (SOD), produzindo peróxido de hidrogênio. Este, por não ser eletricamente carregado, pode atravessar facilmente a membrana plasmática, atingindo o ambiente intracelular. Os ânios superóxido e o peróxido de hidrogênio não são muito reativos, mas podem ser convertidos aos radicais altamente destrutivos hidroxiperoxil (HO2) e hidroxila (OH-), respectivamente, sendo este último formado na presença de íons Fe2+. Todas estas ROS podem resultar em dano considerável ao hospedeiro e ao patógeno, ao reagirem com proteínas, reduzindo a atividade de enzimas; com lipídeos, desestabilizando as membranas; e com o DNA, podendo causar mutações (HAMMOND-KOSACK & JONES, 1996). As ROS, além de desencadearem a morte celular vinculada à HR e controlarem o desenvolvimento do patógeno, ainda: reforçam a parede celular, uma vez que peróxido de hidrogênio é essencial para a formação de precursores de lignina, e radicais de oxigênio podem oxidar as glicoproteínas da parede celular, as quais são então ligadas à mesma através de ligações cruzadas; ativam uma cascata de sinalização mediada pelo ácido salicílico (“Salicilic Acid” ou SA), uma vez que o peróxido de hidrogênio aumenta a atividade da enzima BA2-H envolvida na biossíntese 15 do SA; ativa mecanimos de proteção celular contra sua própria toxicidade; e levam à alteração do potencial redox na célula, o qual é conhecido por ativar fatores de transcrição, como o NF-kB e o AP-1, assim desempenhando uma atividade de sinalização (HAMMOND-KOSACK & JONES, 1996). Como um resultado do acúmulo de ROS, são ativados mecanismos de proteção celular contra os seus efeitos destrutivos. Entre estes mecanismos estão a produção de SODs; o ciclo do ascorbato/glutationa, que remove as ROS dos protoplastos; a produção de catalases que convertem o H2O2 em H2O e O2 nos peroxissomos e glioxissomos, entre outros (RESENDE et al., 2003). O óxido nítrico (“Nitric Oxide” ou NO) é outra importante molécula que se acumula durante interações incompatíveis. O NO pode auxiliar a formação da HR por potencializar a indução da morte celular mediada por ROS. É possível que as ROS e o NO reajam formando o altamente tóxico peroxinitrito (ONOO-), mas esta idéia ainda deve ser esclarecida (GRANT & LOAKE, 2000). Um tipo de reforço da parede celular que ocorre rapidamente em resposta à invasão fúngica é a formação de papilas. Elas são formadas imediatamente abaixo do sítio de penetração da hifa de infecção, bloqueando a penetração fúngica das células hospedeiras (BAYLES et al., 1990). É considerada a possibilidade das papilas, ao invés de desempenharem uma função antagônica ao patógeno fúngico, serem induzidas por ele para auxiliar a penetração. A deposição de calose, as anteriomente citadas ligações cruzadas de glicoproteínas à parede celular e lignificação adicional, entre outros, são modificações da parede celular que foram relatadas como parte das respostas iniciais ou tardias ativadas em interações incompatíveis (HAMMOND-KOSACK & JONES, 1996). A produção de ácido salicílico (“Salicilic Acid” ou SA), ácido benzóico (“Benzoic Acid” ou BA), etileno e de ácido jasmônico (“Jasmonic Acid” ou JA), está relacionada à resposta a fitopatógenos e insetos. Em locais em que ocorre a incompatibilidade mediada pelo reconhecimento R-Avr, níveis elevados de SA garantem que respostas de defesa, como a indução de genes PR, sejam ativadas, inibindo o desenvolvimento do patógeno. Uma proteína envolvida na produção de SA em plantas é a fenilalanina 16 amonialiase (“phenilalanine amonialyase” ou PAL), a qual tem sua produção aumentada após ataque com fungos (TAIZ & ZEIGER, 2004). Entrentanto, a forma mais observada de produção de SA envolve a liberação de BA, que induz a expressão de um citocromo P450 monoxigenase (BA2-H), o qual converte BA em SA. A atividade de BA2-H é fortemente aumentada após HR (LÉON et al., 1995). O etileno e o JA são moléculas sinalizadoras envolvidas na resposta de defesa a estresses bióticos e abióticos. Eles induzem à expressão de compostos antimicrobianos não induzidos por SA. Acredita-se que estas moléculas atuem em uma cascata antagonista à cascata envolvendo SA, embora possa haver cruzamento entre estas duas rotas de sinalização (DONG, 1998). Alguns exemplos do papel do SA, JA e etileno são dados a seguir. Os transcritos de inúmeros genes PR acumulam dentro de minutos a horas do ataque pelo patógeno. Geralmente os mesmos genes são induzidos em interações compatíveis e incompatíveis, entretanto mais fraca e vagarosamente nas interações compatíveis. Algumas proteínas PR são quitinases e glucanases, enzimas que degradam polissacarídeos estruturais das paredes das células fúngicas e provavelmente reduzem o crescimento fúngico. O etileno e o SA parecem agir sinergicamente levando ao aumento das proteínas PR. Outras proteínas PR, como as lipoxigenases (LOX), geram moléculas sinalizadoras secundárias tais como JA e peróxidos lipídicos, e produzem metabólitos secundários tóxicos voláteis e não voláteis. As defensinas são também proteínas com atividade antimicrobiana. A via de sinalização que controla o acúmulo de defensinas é mediada por etileno e JA (HAMMOND- KOSACK e JONES, 1996). As fitoalexinas são um grupo de metabólitos secundários quimicamente diverso que se acumula em torno do local de infecção e apresenta atividade antimicrobiana. Os isoflavonóides constituem um grupo de fitoalexinas comuns em leguminosas. Embora as fitoalexinas acumulem-se em concentrações tóxicas aos patógenos em bioensaios, o significado destes compostos para a defesa da planta ainda não é bem compreendido (TAIZ & ZEIGER, 2004). 17 Poucos minutos após o reconhecimento R-Avr, rotas complexas de sinalização são ativadas. Um elemento inicial comum dessas cascatas é a mudança transitória da permeabilidade iônica da membrana plasmática vegetal. A ativação dos genes R estimula a entrada dos íons Ca2+ e H+ na célula e a saída de K+ e Cl- (NÜMBERGER & SCHEEL, 2001). O rápido acúmulo de Ca2+ citosólico regula o complexo NADPH oxidase diretamente ou indiretamente, através de aumentada produção de NADPH pela proteína de ligação a Ca2+, calmodulina. Assim, o acúmulo de Ca2+ estaria envolvido na produção de ROS (GRANT & LOAKE, 2000). Presume-se que o Ca2+ também possa sinalizar outras reações de defesa. Quando uma planta sobrevive ao ataque de um patógeno, muitas vezes mostra um aumento na sua resistência a ataques subsequentes em qualquer outra região não atacada e apresenta proteção contra uma variedade de espécies de patógenos. O fenômeno chamado de resistência sistêmica adquirida (“Systemic Aquired Resistance” ou SAR) desenvolve-se após um período de vários dias após a infecção inicial (RYALS et al., 1996). A SAR parece resultar dos níveis aumentados de certos compostos de defesa já mencionados, incluindo quitinases e outras enzimas hidrolíticas. Apesar dos mecanismos de indução de SAR serem ainda desconhecidos, um dos sinais endógenos provavelmente é o SA. Embora o SA não seja um sinal móvel, o seu éster metil, o metil salicilato (no inglês MeSA), age como um sinal volátil que é transmitido para as partes distantes da planta e até para plantas adjacentes (SHULAEV et al., 1997). Este processo parece funcionar contra ataques de vírus, fungos, bactérias e nematóides. Um resumo das principais vias de transdução de sinal necessárias para ativar e coordenar a indução localizada de respostas de defesa é apresentado na Figura 3. Apesar de uma série de mudanças metabólicas e estruturais em diversos patossistemas serem conhecidas, os mecanismos envolvidos na resposta de defesa da soja à ferrugem asiática ainda não foram elucidados. Estas informações serão muito úteis para a elaboração de estratégias efetivas de melhoramento de cultivares para a resistência durável. 18 Figura 3: Visão geral das vias de transdução de sinal que ativam e coordenam as respostas de defesa das plantas. Adaptado de HAMMOND-KOSACK & JONES (2000). 2.5 Mapeamento Genético de Genes de Resistência O mapa genético de uma espécie é um modelo abstrato do arranjo linear de um grupo de genes ou marcadores (SCHUSTER & CRUZ, 2004). Marcadores são fatores morfológicos, fisiológicos, bioquímicos ou genéticos passíveis de serem identificados, herdáveis e que permitem o estudo comparativo de genótipos e suas progênies. Dentre os tipos de marcadores existentes, os marcadores de DNA são os mais interessantes, devido à ausência de efeito ambiental, abundância e alto polimorfismo (FERREIRA & GRATTAPAGLIA, 1998). O desenvolvimento e uso de marcadores moleculares para a detecção e exploração de polimorfismos de DNA em sistemas animais e vegetais é um dos avanços mais significativos no campo da biologia molecular e biotecnologia. Ele 19 impulsionou a pesquisa genômica de plantas durante o último quarto de século, e tornou o uso de marcadores moleculares uma área de pesquisa em ascensão na genética vegetal (GUPTA et al., 2005). Os marcadores moleculares começaram a ser empregados na década de 80 e foram extensivamente utilizados para acessar a diversidade genética de populações, caracterizar germoplasmas, mapear locos de característica qualitativas e quantitativas, e auxiliar na introgressão de genes em cultivares, processo denominado Seleção Assistida por Marcadores (“Marker Assisted Selection” ou MAS). Os marcadores de DNA mais utilizados até o momento para estes propósitos podem ser classificados em dois tipos principais: marcadores moleculares loco não específicos dominantes, como o RAPD (“Random Amplified Polymorphic DNA” - WILLIAMS et al., 1990) e o AFLP (“Amplified Fragment Length Polymorphism” - VOS et al., 1995); e marcadores moleculares loco-específicos codominantes, como o RFLP (“Restriction Fragment Length Polymorphism”- GRODZIKER et al., 1974) e os SSR (“Simple Sequence Repeats” - LITT & LUTTY, 1989). Apesar destes terem sido os mais utilizados, atualmente existem outros tipos de marcadores, como os SNPS (“Single Nucleotide Polymorphism”) e os ESTs-SSR (EST: “Expressed Sequence Tags” ou Etiquetas de Seqüências Expressas), os quais têm sido aplicados com sucesso no mapeamento de características importantes. Os SSR ou microssatélites baseiam-se na amplificação via PCR (“Polymerase Chain Reaction” ou Reação em Cadeia da Polimerase) de seqüências simples repetidas utilizando iniciadores específicos complementares às regiões que flanqueiam essas seqüências. No genoma dos organismos eucariontes existem muitas seqüências repetitivas de DNA que se localizam em regiões de microssatélites. As seqüências de DNA que flanqueiam essas regiões contendo repetições curtas são conservadas dentro de uma espécie, permitindo a seleção de iniciadores de PCR que possam ser utilizados para amplificar os SSRs. Após a amplificação dessas regiões específicas, o polimorfismo no tamanho dos fragmentos obtidos é devido ao número diferente de repetições das seqüências simples (FERREIRA & GRATTAPAGLIA, 1998). Em um 20 estudo de SSR utilizando 34 espécies de plantas, o elemento repetido mais comum foi o dinucleotídio AT (MORGANTE & OLIVIERI, 1993) Os marcadores SSR são altamente polimórficos, abundantes e distribuídos por todo o genoma, e sua análise é passível de automação. Com o desenvolvimento e a publicação dos iniciadores para amplificação de SSRs da soja, marcadores para todos os grupos de ligação puderam ser mapeados (CREGAN et al., 1999; SONG et al., 2004). Em soja, os marcadores SSR têm sido muito utilizados para o mapeamento de genes específicos que determinam características agronomicamente importantes, como resistência a doenças e, também, para a identificação de locos de características quantitativas (“Quantitative Trait Loci” ou QTLs) de importância econômica, envolvidos na produtividade de grãos e na resistência genética a pragas e doenças (FRANCIA et al., 2005). A grande disponibilidade de marcadores neutros, altamente polimórficos, cuja herança não sofre influência ambiental, aliada a procedimentos estatísticos complexos, tem permitido a construção de mapas de ligação para a maioria das espécies vegetais de interesse agronômico, até mesmo para aquelas de longo ciclo de vida, como as florestais e as frutíferas. A construção de mapas de ligação tem como base a análise de segregação de centenas de marcadores e, por isso, é computadorizada. Programas como “Mapmaker” (LANDER et al., 1987), “Linkage 1” (SUITER et al., 1983), “Gmendel” (LIU & KNAPP, 1990) e GQMOL 9.1 (CRUZ & SCHUSTER, 2006), foram desenvolvidos para auxiliar na análise genética dos dados visando à construção de mapas genéticos (CARNEIRO & VIEIRA, 2002). A seleção da população para mapeamento envolve a escolha de genitores, os quais devem ser suficientemente contrastantes para garantir o máximo de polimorfismo, e a determinação do tipo de população, sendo considerada uma etapa crítica para o sucesso da construção do mapa (STAUB et al., 1996). Diversos tipos de populações podem ser utilizadas para o mapeamento. Para espécies vegetais passíveis de autofecundação, como a soja, podem ser utilizadas as populações F2, F3 ou Fn, os duplo haplóides, as linhagens endogâmicas recombinantes (“Recombinant Imbred Lines” ou RILs) e populações derivadas de retrocruzamento. A utilização de população 21 F2 tendo cada indivíduo uma progênie F3 individual (F2:3) tem as vantagens de facilidade e rapidez de obtenção e disponibilidade de informação sobre os três possíveis genótipos para um loco (AA, Aa, aa), o que aumenta a precisão dos componentes genéticos (SCHUSTER & CRUZ, 2004). Somente é possível identificar os indivíduos heterozigotos (Aa) para caracteres que apresentam herança de dominância completa (como a resistência à ferrugem asiática) quando a progênie de cada indivíduo F2 é avaliada. Se na progênie de um dado indivíduo F2 for observado somente o fenótipo dominante, este indivíduo é homozigoto dominante; se for detectado somente o fenótipo recessivo, o indivíduo F2 é homozigoto recessivo; mas, se houver segregação para os fenótipos dominante e recessivo, o indivíduo F2 é heterozigoto. Para garantir a confiabilidade deste teste de progênie, elas devem ser formadas por um número adequado de indivíduos. Para caracteres de herança monogênica completamente dominante, 15 indivíduos na progênie gera uma confiança de aproximadamente 98,7% na avaliação. Para se obter esta estimativa, utilisa-se a fórmula: n=[ln(1-p)]/[ln(1-P)], onde “n” é o tamanho da população, “p” é a probabilidade de que o genótipo almejado ocorra (que, neste caso, é o genótipo recessivo, com probabilidade igual a 25%) e “P” é a probabilidade desejada de ocorrência de pelo menos um indivíduo com o genótipo almejado (RAMALHO et al., 2004). Populações F2:3 têm sido frequentemente utilizadas no mapeamento de uma série de genes em diversas culturas. A escolha do tipo de população mais apropriada ao mapeamento está relacionada com o tipo de marcador empregado (TANKSLEY et al., 1988). O máximo de informação genética é alcançado quando se usa uma população F2 (ou F2:3) e marcadores codominantes (STAUB et al., 1996). O tamanho da população também é algo determinante para o sucesso do trabalho. O tamanho da população de mapeamento pode ser definido considerando-se a precisão desejada para a estimativa de frequência de recombinação. Quanto maior o tamanho da população, maior será a resolução do mapa e, consequentemente, menores frequências de recombinação poderão ser estimadas. Populações F2 com marcadores codominantes são aquelas que demandam o menor tamanho de população. Para populações F2 de cerca de 150 indivíduos avaliadas com marcadores codominantes, pode se obter frequências de 22 recombinação de 5%, num intervalo de confiança de 5% e num nível de confiança de 95% (SCHUSTER & CRUZ, 2004). A busca por marcadores ligados ao caráter que se pretende mapear é uma etapa decisiva no trabalho de mapeamento. Uma estratégia interessante é a análise de segregantes agrupados (“Bulked Segregant Analysis” ou BSA - MICHELMORE et al., 1991). Este método consiste na avaliação de dois grupos de indivíduos obtidos de pontos extremos da curva de distribuição do caráter em estudo. Espera-se que os alelos dos marcadores que estiverem associados aos alelos do caráter que está sendo mapeado estejam distribuídos de forma desigual entre os dois grupos, possibilitando sua detecção. A BSA tornou a identificação de marcadores moleculares ligados a caracteres de interesse um processo mais eficiente. Este método tem sido utilizado com sucesso na detecção de marcadores moleculares associados a genes de resistência a doenças (ZHAO et al., 2005). Os caracteres mapeados devem distribuir-se de acordo com um padrão de segregação mendeliana. Este padrão pode ser verificado através do teste do Qui- quadrado (�2). É recomendado, preferencialmente, o descarte dos locos que apresentam distorções da segregação mendeliana para não comprometer a qualidade do mapa (BEARZOTI, 2000). A falta de aditividade das frequências de recombinação levou ao desenvolvimento das funções de mapeamento. Elas são utilizadas para converter frequências de recombinação, expressas em unidades de mapa (centiMorgan ou cM) em medidas de distância com propriedades mais interessantes para o ordenamento dos locos. As mais conhecidas são as de Haldane (HALDANE, 1919) e de Kosambi (KOSAMBI, 1944). A função de Haldane admite a independência das permutas nos intervalos adjacentes, enquanto a função de Kosambi considera a interferência. Uma vez estimadas as frequências de recombinação e estabelecida a função de mapeamento, a etapa seguinte é formar os grupos de ligação contendo os marcadores ligados. Dois marcadores são considerados ligados quando a frequência de recombinação entre eles for inferior a um limite predefinido e/ou o LOD escore (“Logarithm of ODds” ou Logarítmo dos Nós) for inferior a um limite também predefinido. 23 A distância máxima de ligação utilizada como padrão em vários programas de mapeamento é 37,2 cM e o LOD escore igual a 3,0. Um LOD escore igual a 3 indica que a ocorrência de ligação é mil vezes mais provável que a de segregação independente (BEARZOTI, 2000; SCHUSTER & CRUZ, 2004). 2.6 Métodos de Análise da Expressão Gênica Diferencial em Situação de Doença A expressão de genes eucarióticos que produzem proteínas pode ser dividida em três etapas: 1) a transcrição de um filamento de pré-mRNA utilizando a seqüência de um gene como molde, 2) o processamento do pré-mRNA anterior à sua translocação para o citoplasma, 3) a tradução do filamento de mRNA em uma cadeia polipeptídica, e 4) as modificações pós-traducionais, incluindo o dobramento do polipeptídeo na conformação correta para sua atividade biológica e outras modificações, como a adição específica de grupos de carboidratos. A primeira etapa da expressão gênica tem sido a mais estudada com o objetivo de qualificar e quantificar um gene ou conjuntos de genes diferencialmente expressos por uma célula, tecido ou organismo numa dada condição. A coleção completa de mRNAs transcritos em uma célula, tecido ou órgão em um momento particular é denominada de transcritoma e o estudo dos mRNAs transcritos é denominado transcritômica (DALE & SCHANTZ, 2002). O estudo do transcritoma de uma célula ou tecido representa o primeiro passo na caracterização funcional de genes de interesse e identificação de locos corregulados. As metodologias utilizadas para estabelecer perfis transcricionais geralmente baseiam-se na geração de populações de cDNA a partir da população de mRNA expressa. O primeiro grupo de técnicas de análise de transcritomas baseia-se nas várias formas de acessar a quantidade de um transcrito específico em uma dada amostra em um determinado momento. Dentre as técnicas que compõem este grupo destacam-se os “Northern Blots” (ALWINE et al., 1977), Ensaio de Proteção de Ribonuclease (“Ribonuclease Protection Assay” ou RPA - MELTON et al., 1984; WINTER et al., 24 1985), Transcrição Reversa (“Reverse Transcription” ou RT) seguida de PCR em Tempo Real (“Real Time RT-PCR” - TaqMan - GIBSON et al., 1996) e Ensaio de Extensão de Iniciador (“Primer Extension Assay” - BOORSTEIN & CRAIG, 1989). O segundo grupo de técnicas se preocupa com a análise global de transcritos. Dentre as técnicas compreendidas neste grupo, é possível citar a Análise Serial da Expressão Gênica (“Serial Analysis of Gene Expression” ou SAGE - VELCULESCU et al., 1995), Seqüenciamento de Assinaturas Massivamente Paralelo (“Massively Pararel Signature Sequencing” ou MPSS - BRENNER et al., 2000) e Piroseqüenciamento (WEBER et al., 2007). Um terceiro grupo de técnicas abrange a comparação de diferentes transcritomas e a ciência envolvendo estas técnicas pode ser chamada de transcritômica comparativa. Neste grupo, destacam-se as técnicas “Differential Display” (LIANG & PARDEE, 1992), “RNA fingerprinting” (WELSH et al., 1992), Polimorfismo de Comprimento de Fragmento Amplificado baseado em cDNA (cDNA-AFLP - BACHEM et al., 1996), “GeneFishing” (KIM et al., 2004), Hibridização Subtrativa (BAUTZ & REILLY, 1966; DIATCHENKO et al., 1996) e Microarranjos de DNA (SCHENA et al., 1995). Tais técnicas permitem isolar transcritos que sejam expressos em um dado tecido, fase do desenvolvimento ou condição, como, por exemplo, uma situação de doença sistêmica na planta, utilizando uma amostra controle como referência para comparação. Dentre os métodos apresentados acima, a hibridização subtrativa é o único capaz de enriquecer uma população de cDNAs com transcritos diferencialmente expressos. O método foi primeiramente utilizado por BAUTZ & REILLY (1966) para purificar o mRNA do fago T4. A amostra que possui os transcritos de interesse é chamada testadora (“tester”) e aquela utilizada como controle é chamada referência (“driver”). Embora haja variações no método, a teoria básica por trás da subtração é a mesma: ambas as amostras de mRNA são convertidas a cDNA, os quais são postos para hibridizar; as seqüências híbridas são removidas e as restantes, não-hibridizadas, representam os genes que foram expressos na amostra testadora, mas não foram na referência. Logo, estes são os transcritos diferencialmente expressos entre as amostras (DALE & SCHANTZ, 2002; JI & CAI, 2004). 25 As metodologias tradicionais de hibridização subtrativa são de uso limitado, porque elas requerem grandes quantidades de mRNA inicial e produzem uma ínfima quantidade de cDNA ao final do processo, dificultando sua clonagem. A Hibridização Subtrativa Supressiva (“Supressive Subtractive Hibridization” ou SSH - DIATCHENKO et al., 1996), um método baseado em PCR, contornou as limitações das técnicas anteriores, resultando em um enriquecimento de 1000 vezes na quantidade de transcritos diferencialmente expressos, sendo assim adequado à identificação de transcritos raros (JI & CAI, 2004) Na SSH, o cDNA é sintetizado a partir de mRNA obtido das duas amostras que se pretende comparar: a testadora e a referência. A amostra testadora é subdividida em duas partes, sendo cada uma delas ligada a um adaptador diferente (a seqüência deste adaptador é complementar aos iniciadores que serão utilizados posteriormente para amplificação por PCR). Duas hibridizações são então realizadas. Na primeira, um excesso de cDNA referência é adicionado a cada amostra do testador. As amostras são então desnaturadas por calor e deixadas se complementar, gerando as moléculas tipo ‘a’, ‘b’, ‘c’ e ‘d’ em cada amostra (Figura 4). As moléculas tipo ‘a’ são significativamente enriquecidas com seqüências diferencialmente expressas, enquanto os cDNAs que não são diferencialmente expressos formam as moléculas tipo ‘c’ mediante complementação com o cDNA referência (DIATCHENKO e al., 1996). Na segunda hibridização, as duas amostras da primeira hibridização são misturadas sem serem desnaturadas. Agora, somente os cDNA testadores unifilamentares subtraídos e equalizados podem se reassociar, formando um novo tipo de híbridos: as moléculas tipo ‘e’ (Figura 4). Estes novos híbridos serão amplificados por PCR, gerando a população de cDNAs enriquecida com seqüências diferencialmente expressas. As moléculas tipo ‘a’ e ‘d’ não podem ser amplificadas, por não possuírem o sítio de complementação dos iniciadores (Figura 4). Devido ao efeito de supressão da PCR, a maioria das moléculas tipo ‘b’ formam uma estrutura em grampo, que previne sua amplificação. As moléculas tipo ‘c’ possuem somente um sítio de complementação de primer, e são amplificadas linearmente. As moléculas tipo ‘e’ amplificadas podem ser 26 diretamente inseridas em um vetor T/A, o qual é adaptado para produtos de PCR, para a clonagem e obtenção de uma biblioteca subtrativa (DIATCHENKO e al., 1996). A hibridização subtrativa tem sido muito empregada na identificação de transcritos ativados em situação de doença (DALE & SCHANTZ, 2002). Vários trabalhos têm sido publicados utilizando a SSH para identificar com sucesso transcritos diferencialmente expressos (AKOPIAN & WOOD, 2000; FERNÁNDEZ et al., 2002; HEIN et al., 2004; DEGENHARDT et al., 2005; WANG et al., 2007; ZOUARI et al., 2007). Em alguns destes trabalhos, os cDNA clonados são impressos em micro ou macroarranjos e hibridizados com cDNAs-alvo adequados para a confirmação da expressão diferencial (YANG et al., 1999; ZENGH et al., 2004; RISHI et al., 2004; LI et al., 2007). Na técnica de microarranjos de cDNA (SCHENA et al., 1995), seqüências de cDNA de milhares de genes são organizadas sobre lâminas de vidro. As duas populações de mRNA que se quer comparar são transcritas a cDNA, cada uma das quais na presença de nucleotídeos marcados com fluorocromos diferentes. A seguir, os cDNAs marcados são desnaturados, misturados na mesma proporção e hibridizados às sondas fixadas à lâmina. As lâminas são expostas à luz de comprimentos de onda apropriados e os sinais são registrados em imagens. A comparação entre as duas imagens geradas permite a identificação dos genes diferencialmente expressos. O advento da impressão robótica de sondas em alta densidade e da análise eletrônica de imagens tornou tecnicamente simples produzir tais arranjos e praticamente ilimitada a quantidade de seqüências que podem ser analisadas (DALE & SCHANTZ, 2002). 27 Figura 4: Diagrama esquemático da SSH. Linhas sólidas representam os cDNAs digeridos com Rsa I. Caixas preenchidas representam a parte externa do adaptador 1 e 2R, que são idênticas e complementares ao iniciador da primeira PCR. Caixas vazias representam a parte interna do adaptador 1 e são complementares ao iniciador interno 1 da segunda PCR. Caixas sombreadas representam a parte interna do adaptador 2R e são complementares ao iniciador interno 2R da segunda PCR (BD BIOSCIENCES, 2004). 28 2.7 Estudos Moleculares da Resistência da Soja à Ferrugem Asiática Trabalhos de mapeamento dos locos de resistência à ferrugem asiática têm sido conduzidos. O locus Rpp1 foi mapeado entre os marcadores SSR Sct_187 e Sat_164, à distância de 0,4 cM de ambos os marcadores, no grupo de ligação G da soja (HYTEN et al., 2007). Além do Rpp1, alelos identificados em outras variedades também foram mapeados. O loco identificado na cultivar japonesa Hyuuga foi mapeado entre os marcadores SSR Satt460 e Satt307, no grupo de ligação C2, a distâncias de 0,8 cM e 2,4 cM, respectivamente (MONTEROS et al., 2007). Nesta mesma região genômica, foi mapeado o loco presente na cultivar FT-2, localizando-se entre os marcadores SSR Satt460 e Satt079, a distâncias de 20,3 cM e 7,8 cM, respectivamente (BROGIN et al. 2004; POLIZEL et al., 2006). Segundo trabalho desenvolvido por GILLI et al. (2006), com base em análises de similaridade genética, o alelo presente na FT-2 pode ser o Rpp3. Esforços anteriores na Tailândia também localizaram um loco de resistência no grupo de ligação G, em uma associação estreita com os marcadores Satt012 and Satt472 (NUNTAPUNT et al., 2004), em genótipos de soja previamente submetidos a agentes mutagênicos. Não foi identificada a relação entre os locos presentes na Hyuuga, FT-2 e nas linhas mutagênicas tailandesas, e os genes Rpp descritos anteriormente. E até o momento, não foi publicado trabalho relatando o mapeamento dos genes Rpp2 e Rpp4, os quais são os únicos, dentre os previamente descritos na literatura, a conferir resistência contra a raça fisiológica da ferrugem presente no Brasil. Os genes de resistência podem estar distribuídos aleatoriamente pelo genoma ou em agrupamentos (“clusters”). Vários dos genes de resistência a fitopatógenos clonados até o momento residem em agrupamentos, indicando que um mecanismo genético comum esteja envolvido em sua evolução. Em geral, duplicações gênicas causadas por pareamento desigual seguido de recombinação e subsequente acúmulo de mutações de ponto, estão envolvidos na criação de novos genes de resistência e em seu agrupamento numa dada região genômica. A ligação genética de genes de resistência foi relatada em milho para o agrupamento de genes Rp1, em cevada para o 29 agrupamento de genes Mla, em alface para os locos Dm, em aveia para os genes Pc e em linho para os genes L (KANAZIN et al. 1996). Com base no conhecimento das seqüências mais comumente encontradas nos genes de resistência, foi proposta e confirmada a existência de análogos de genes de resistência (“Resistance Gene Analogs” ou RGAs). Os RGAs são clonados por PCR utilizando iniciadores cuja seqüência é complementar a regiões conservadas nos genes de resistência. Em trabalho desenvolvido por KANAZIN et al. (1996) em soja, 19 RGAS foram clonados e mapeados nos grupos de ligação C2, D1, H, J, L, M, N1 e P. Alguns trabalhos de genômica funcional da soja sob infecção com ferrugem asiática foram realizados recentemente. A análise de expressão gênica de uma cultivar americana de soja através de microarranjo de DNA revelou que a maioria dos genes que mudaram sua expressão em situação de infecção precoce foram genes de resposta de defesa (PANTHEE et al., 2007). Este trabalho, entretanto, foi conduzido utilizando um genótipo para o qual não se conhece a resposta à ferrugem, uma vez que as plantas inoculadas não apresentaram sintomas de infecção. A infecção foi dada como positiva através da amplificação de um fragmento específico do fungo nas amostras de cDNA das folhas inoculadas. Deste modo, não é possível relacionar os genes expressos diferencialmente com um mecanismo de resistência, mas sim com um mecanismo de resposta a infecção. Em outro trabalho de análise de expressão gênica utilizando o mesmo microarranjo comercial (MORTEL et al., 2007), foram analisados dois genótipos de soja, a Embrapa 48 (suscetível) e a PI230970 (resistente, possuidora do gene Rpp2), no decorrer dos primeiros sete dias após a inoculação em intervalos variáveis de 6 h a 1 dia. Foram detectadas mudanças bifásicas distintas na expressão dos mRNAs, com um pico de expressão de genes relacionados à defesa dentro das primeiras 12 h após a infecção, em ambos os genótipos, e outro próximo de 96 h após infecção no genótipo suscetível e próximo de 72 h após infecção no genótipo resistente, indicando uma expressão gênica relativamente precoce nos mecanismos mediados pelo gene Rpp2. Em ambos os trabalhos, a análise de microarranjos foi conduzida utilizando microarranjos de oligonucleotídeos comerciais (“Soybean Genome GeneChip Array”, 30 Affymetrix), os quais representam mais de 37.500 transcritos da soja. As sondas de oligonucleotídeos presentes no microarranjo comercial são restritas a certos germoplasmas e condições de cultivo, o que representa uma limitação da técnica. Apesar de terem sido incluídas na fabricação destes microarranjos oligonucleotídeos representando genes de resposta ao nematóide do cisto (Heterodera glycines) e à podridão da raiz (Phytophthora sojae), não foram representados genes de defesa a fungos biotróficos, para os quais as plantas apresentam mecanismos específicos de defesa. Assim, é provavel que estes microarranjos não contemplem algumas ou várias das seqüências ativadas em interações compatíveis e incompatíveis entre soja e ferrugem asiática. Até o momento, não foram publicados trabalhos relatando a clonagem de genes envolvidos na reação da soja à ferrugem asiática. Trabalhos de genômica funcional do fungo também foram realizados. Uma biblioteca de cDNA unidirecional foi construída utilizando mRNA de uredósporos germinados de P. pachyrhizi. Foram identificadas 488 ESTs únicas, 189 das quais apresentaram similaridade significativa com seqüências depositadas na base de dados de proteínas não-redundantes do NCBI (“National Center for Biotechnology Information” - http://www.ncbi.nlm.nih.gov/ - POSADA-BUITRAGO & FREDERICK, 2005). A maior parte dos transcritos identificados correspondeu, nesta ordem, a proteínas do metabolismo, expressão gênica, estrutura celular e crescimento, divisão celular, sinalização celular e defesa. 31 3 OBJETIVOS No presente trabalho é descrita a utilização de ferramentas da biologia molecular para auxiliar a compreensão dos mecanismos envolvidos na resistência à ferrugem asiática e para facilitar o melhoramento de cultivares de soja para a resistência durável a este patógeno. Assim, os objetivos deste trabalho foram: 3.1 Objetivo Geral - Acessar informações moleculares da soja ao nível do genoma e transcritoma, a fim de acelerar a obtenção de variedades resistentes à ferrugem asiática. 3.2 Objetivos Específicos - Identificar marcadores moleculares ligados aos genes de resistência à ferrugem asiática Rpp2 e Rpp4, possibilitanto seu mapeamento nos grupos de ligação da soja. - Caracterizar as mudanças transcricionais ativadas durante interações compatível e incompatível da planta de soja com o fungo da ferrugem asiática, bem como isolar novas seqüências da soja que respondam ao ataque por este fungo. 32 4 REFERÊNCIAS BIBLIOGRÁFICAS AGÊNCIA DA NOTÍCIA. Soja resistente à ferrugem asiática será apresentada a produtores de MT. Disponível em . Acesso em 13 de março de 2008. AKOPIAN, A. N.; WOOD, J. Construction and screening of a subtractive cDNA library. In: HUNT, S.; LIVESEY, F. (Ed.). Functional genomics: a practical approach. New York: Oxford, 2000. p. 272. ALWINE, J. C.; KEMP, D. J.; STARK, G. R. Method for detection of specific RNAs in agarose gels by transfer to diabenzyloxymethyl-paper and hybridization with DNA probes. Proceedings of the National Academy of Sciences of the United States of America, Washington, v. 74, n. 12, p. 5350-5354, Dec. 1977. ASSMANN, S. M. G protein regulation of disease resistance during infection of rice with rice blast fungus. Science, Washington, v. 2005, n. 310, p. 13-14, Nov. 2005. AVRDC. Physiological races of soybean rust. Progress Report 1983. Shanhua: ROC, 1985. p. 253-257. BACHEM, C.; VAN DER HOEVEN, R. S.; DE BRUIJN, S. M.; VREUQDENHIL, D.; ZABEAU, M.; VISSER, R. G. Visualization of differential gene expression using a novel method of RNA finger-printing based on AFLP: analysis of gene expression during potato tuber development. Plant Journal, Oxford, v. 9, n. 5, p. 745-753, May 1996. BAUTZ, E. K.; REILLY, E. Gene-specific messenger RNA: isolation by the deletion method. Science, Washington, v. 151, n. 21, p. 328, Jan. 1966. 33 BAYLES, C. J.; GHEMAWAT, M. S.; AIST, J. R. lnhibition by 2-deoxy-D-glucose of callose formation, papillae deposition, and resistance to powdery mildew in an mlo barley mutant. Physiological and MoIecular Plant Pathology, London, v. 36, n. 1, p. 63-72, Jan. 1990. BD BIOSCIENCES. BD PCR-Select™ cDNA Subtraction Kit User Manual: Cat. No. 637401. Palo Alto, 2004. Disponível em . Acesso em 29 nov. 2005. BEARZOTI, E. Mapeamento de QTL. In: PINHEIRO, J.B.; CARNEIRO, I.F (Ed.). Análise de QTL no melhoramento de plantas. Goiânia: FUNAPE, 2000. p.63-223. BERGAMIN FILHO, A. Epidemiologia comparativa: ferrugem da soja e outras doenças. In: ZAMBOLIN, L. (Ed.), Ferrugem asiática da soja. Viçosa: Suprema, 2006. 139 p. BERTRAND, J-P.; LAURENT, C.; LECLERCQ, V. O mundo da soja. São Paulo: HUCITEC, 1987. 139 p. BOWEL G. P.;BINDSCHEDLER, L. V.; BLEE, K. A.; BUTT, V. S.; DAVIES, D. R.; GARDNER, S. L.; GERRISH, C.; MINIBAYEVA, F. The apoplastic oxidative burst in response to biotic stress in plants: a three-component system. Journal of Experimental Botany, Oxford, v. 53, n. 372, p. 1367-1376, May 2002. BOORSTEIN, W. R.; CRAIG, E.A. Primer extension analysis of RNA. Methods of Enzymology, New York, v. 180, p. 347–69, 1989. BRASIL. Ministério da Agricultura Pecuária e Abastecimento. Serviço Nacional de Proteção de Cultivares: certificados de proteção concedidos no período de 01/01/1998 a 07/11/2007. Brasília: MAPA, 2007. Disponível em: . Acesso em: 1º mar. 2007. 34 BRENNER, S.; JOHNSON, M.; BRIDQHAM, J. GOLDA, G.; LLOYD, D. H.; JOHNSON, D.; LUO, S.; MCCURDY, S.; FOY, M.; EWAN, M.; ROTH, R.; GEORGE, D.; ELETR, S.; ALBRECHT, G.; VERMAAS, E.; WILLIAMS, S. R.; MOON, K.; BURCHAM, T.; PALLAS, M.; DUBRIDGE, R. B.; KIRCHNER, J.; FEARON, K.; MAO, J.; CORCORAN, K. Gene expression analysis by massively parallel signature sequencing (MPSS) on microbead arrays. Nature Biotechnology, New York, v. 18, n. 6, p. 630-634, June 2000. BROGIN, R. L.; ARIAS, C. A. A.; VELLO, N. A.; TOLEDO, J. F. F.; PIPOLO, A. E.; CATELLI, L. L.; MARIN, S. R. R. Molecular mapping of a gene conferring resistance to soybean rust In: World Soybean Research Conference, 7., 2004, Iguasu Falls. Abstracts… Iguasu Falls: Embrapa, 2004, p. 318. BROMFIELD, K.R. Soybean rust. Saint Paul: American Phytopathological Society, 1984. 65 p. BROMFIELD, K. R.; HARTWIG, E. E. Resistance to soybean rust and mode of inheritance. Crop Science, Madison, v. 20, p. 254-255, Mar.-Apr. 1980. CALDWELL, P.M. & MACLAREN. Soybean rust research in South Africa. In: World Soybean Research Conference, 7., 2004, Iguasu Falls. Proceedings… Iguasu Falls: Embrapa, 2004. p. 354-60. CARNEIRO, M. S.; VIEIRA, M. L. C. Mapas genéticos em plantas. Bragantia, Campinas, v. 61, n. 2, p. 89-100, 2002. CATANZARITI, A.-M.; DODDS, P. N.; LAWRENCE, G. J.; AYLIFFE, M. A.; ELLISA, J. G. Haustorially Expressed Secreted Proteins from Flax Rust Are Highly Enriched for Avirulence Elicitors. Plant Cell, Rockville, v. 18, n.1, p. 243–256, Jan. 2006. CHENG, Y. W.; CHAN, K. L. The breeding of ‘Tainung 3’soybean. Journal of Taiwan Agricultural Research, Taiwan, v. 17, p. 30-35, 1968. 35 CONAB. Indicadores da agropecuária: outubro de 2007. Brasília: CONAB/ DIGEM/SUINF/GEINT, Ano XVI, n. 10., 2007. 64 p. CREGAN, P. B.; JARVIK, T.; BUSH, A. L.; SHOEMAKER, R. C.; LARK, K. G.; KAHLER, A. L.; KAYA, N.; VANTOAI, T. T.; LOHNES, D. G.; CHUNG, L.; SPECHT, J. E. An integrated genetic linkage map of the soybean genome. Crop Science, Madison, v. 39, n. 5, p. 1464-1490, Set.-Oct. 1999. CRUZ, C.D.; SCHUSTER, I. Estatística genômica: aplicada a populações derivadas de cruzamentos controlados. Viçosa: UFV, 2004. p. 568. CRUZ, C. D.; SCHUSTER, I. GQMOL: aplicativo computacional para análise de dados moleculares e de suas associações com caracteres quantitativos, Versão 9.1. 2006 Disponível em Acesso 01 dez. 2007. DALE, J. W.; VON SCHANTZ, M. From genes to genomes: concepts and applications of DNA technology. New York: Wiley, 2002. p. 360. DANGL, J. L.; JONES, J. D. G. Plant pathogens and integrated defence responses to infection. Nature, London, v. 411, n. 14, p. 826-833, June 2001. DEGENHARDT, J.; AL-MASRI, A. N.; OGLU, S. K. R.; SZANKOWSKI, I.; GAU, A. E. Characterization by suppression subtractive hybridization of transcripts that are differentially expressed in leaves of apple scab-resistant and susceptible cultivars of Malus domestica. Molecular Genetics and Genomics, Berlin, v. 173, p. 326-335, 2005. DIATCHENKO, L.; LAU, Y. F.; CAMPBELL, A. P.; CHENCHIK, A.; MOQADAM, F.; HUANG, B.; LUKYANOV, S.; LUKYANOV, K.; GURSKAYA, N.; SVERDLOV, E. D.; SIEBERT, P. D. Suppression subtractive hybridization: a method for generating differentially regulated or tissue-specific cDNA probes and libraries. Proceedings of the 36 National Academy of Sciences of the United States of America, Washington, v. 93, n. 11, p. 6025-6030, June 1996. DONG, X. SA, JA, ethylene, and resistance in plants. Current Opinion on Plant Biology, v.1, n. 4, p. 316-323, Aug. 1998. ELLIS, J.; CATANZARITI, A.-M.; DODDS, P. The problem of how fungal and oomycete avirulence proteins enter plant cells. Trends in Plant Science, London, v. 11, n. 2, p. 61-63, Feb. 2006. EMBRAPA. Situação da Ferrugem Asiática da Soja no Brasil e na América do Sul (Documentos). Londrina: Embrapa Soja, 2004. 30 p. EMBRAPA. Tabela de custo ferrugem. Disponível em . Acesso 5 set. 2007a. EMBRAPA. Prejuízo com ferrugem da soja se mantém estável nesta safra. Disponível em Acessado em 05 de setembro de 2007b. EMBRAPA. Diferentes usos do grão. Disponível em Acesso 1 dez. 2007c. FEHR, W. R.; CAVINESS, C. E. Stages of soybean development: Special Report 80. Ames: ISU, 1977. FERNÁNDEZ, C.; GREGORY, W. F.; LOKE, P.; MAIZEL, S. R. M. Full-length-enriched cDNA libraries from Echinococcus granulosus contain separate populations of oligo- capped and trans-spliced transcripts and a high level of predicted signal peptide 37 sequences. Molecular and Biochemical Parasitology, Amsterdam, v. 122, n. 2, p. 171-180, Jul. 2002. FERREIRA, M.E; GRATTAPAGLIA, D. 1998. Introdução ao uso de marcadores moleculares em análise genética. 3. ed., Brasília: Embrapa-Cenargen. FLOR, H. H. Current status of gene-a-gene concept. Annual Review of Phytopathology, Palo Alto, v. 9, p. 275-296, Sep. 1971. FRANCIA, E.; TACCONI, G.; CROSATTI, C.; BARABASCHI, D.; BULGARELLI, D.; DALL’AGLIO, E.; VALE, G. Marker assisted selection in crop plants. Plant Cell, Tissue and Organ Culture, Dordrecht, v. 82, n. 3, p. 317–342, Sep 2005. GIBSON, U. E.; HEID, C. A.; WILLIAMS, P. M. A novel method for real time quantitative RT-PCR. Genome Research, v. 6, p. 995-1001, Oct. 1996. GILLI, J. R.; SEQUIN, L.; SOLDINI, D. O.; MASIERO, B. Estudios de dos regiones del genoma de soja relacionadas con la resistencia genética a la roya asiática. In: Congreso de Soja Del Mercosur, 3., 2006, Rosário. Resumenes… Rosário: ACSOJA, 2006. GRANT, J. J.; LOAKE, G. J. Role of reactive oxygen intermediates and cognate redox signaling in disease resistance. Plant Physiology, Washington, v. 124, n. 1, p. 21-29, Sep 2000. GRODZICKER, T.; WILLIAMS, J.; SHARP, P.; SAMBROOK, J. Physical mapping of temperature-sensitive mutations of adenoviruses. Cold Spring Harbor Symposium of Quantitative Biology, v. 39, p. 439-446, 1974. 38 GUPTA, P. K.; RUSTGI, S.; KULWAL, P. L. Linkage disequilibrium and association studies in higher plants: present status and future prospects. Plant Molecular Biology, Dordrecht, v. 57, n. 4, 461-485, Mar. 2005. HALDANE, J.B.S. The combination of linkage values, and the calculation of distance between the loci of linked factors. Journal of Genetics, Bangarole, v.8, p.299-309, 1919. HAMMOND-KOSACK, K. E.; JONES, J. D. G. Resistance gene-dependent plant defense response. Plant Cell, Rockville, v. 8, n. 10, p.1773-1791, Oct. 1996. HAMMOND-KOSACK, K. E.; JONES, J. D. G. Responses to plant pathogens. In: BUCHANAN, B. B.; GRUISSEM, W.; JONES, R. L. (Ed.). Biochemistry and Molecular Biology of Plants. Rockville: American Society of Plant Physiologists, 2000.1367p. HASSE, G.; BUENO, F. O Brasil da soja: abrindo fronteiras, semeando cidades. Porto Alegre: Ceval Alimentos, 1996. 256 p. HARTWIG, E. E. Identification of a fourth major gene conferring resistance to soybean rust. Crop Science, Madison, v. 26, n. 6, p. 1135-1136, 1986. HEIN, I.; CAMPBELL, E. I.; WOODHEAD, M.; HEDLEY, P. E.; YOUNG, V.; MORRIS, W. L.; RAMSAY, L.; STOCKHAUS, J.; LYON, G.D.; NEWTON, A. C.; BIRCH, P. R. J. Characterization of early transcriptional changes involving multiple signaling pathways in the Mla13 barley interaction with powdery mildew (Blumeria graminis f. sp. hordei). Planta, Berlin, v. 218, n. 5, p. 803–813, Mar. 2004. HENNING, A.A.; GODOY, C. V. Situação da ferrugem da soja no Brasil e no mundo. In: ZAMBOLIN, L. (Ed.). Ferrugem asiática da soja. Viçosa MG: Suprema, 2006. 139 p. 39 HIDAYAT, O. O.; SOMAATMADJA, S. Screening of soybean breeding lines for resistance to soybean rust (Phakopsora pachyrhizi Sydow). Soybean Rust Newsletter, v. 1, p. 9-22, 1977. HYTEN, D. L.; HARTMAN, G. L.; NELSON, R. L.; FREDERICK, R. D.; CONCIBIDO, V. C.; NARVEL, J. M.; CREGAN, P. B. Map location of the Rpp1 locus that confers resistance to soybean rust in soybean. Crop Science, Madison, 47, p. 837-838, Mar.- Apr. 2007. JI, W.; CAI, L. Suppression Subtractive Hybridization PCR. In: WEISSENSTEINER, T.; GRIFFIN, H. G.; GRIFFIN, A. (Ed.). PCR Technology: Current Innovations, 2. ed., 2004. p. 287-298. JOHAL, G. S.; BRIGGS, S. P. Reductase activity encoded by the HM1 disease resistance gene in maize. Science, Washington, v. 258, n. 5084, p. 985-987, Nov. 1992. KATO, M.; YORINORI, J. T.; PAIVA, W. M.; YAMAOKA, Y.; ALMEIDA, A. M. R. Epidemiological studies of soybean rust in South America: host suscetibility and pathogenic races. Jircas Working Report, Tsukuba, v. 51, p. 35-46, p. 2007. KANAZIN, V.; MAREK, L. F.; SHOEMAKER, R.C. Resistance Gene Analogs are conserved and clustered in soybean. Proceedings of the National Academy of Sciences of the United States of America, Washington, v. 93, n. 21, p. 11746-11750, Oct. 1996. KEMEN, E.; KEMEN, A. C.; RAFIQI, M.; HEMPEL, U.; MENDGEN, K.; HAHN, M.; VOEGELE, R. Identification of a protein from rust fungi transferred from haustoria into infected plant cells. Molecular Plant-Microbe Interaction, Saint Paul, v. 18, n. 11, 1130-1130, Nov. 2005. 40 KIM, Y.-J.; KWAK, C. -I.; GU, Y.-Y.; HWANG, I.-T.; CHUN, J.-Y. Annealing control primer system for identification of differentially expressed genes on agarose gels. Biotechniques, Natick, v. 36, n. 3, p. 424-434, Mar. 2004. KNOGGE, W. Fungal infection of plants. Plant Cell, Rockville, v. 8, n. 10, 1711-1722, Oct. 1996. KOSAMBI, D.D. The estimation of map distance from recombination values. Annuaire of Eugenetics, v. 12, 172-175, 1944. LANDER, E. S.; GREEN, P.; ABRAHAMSON, J.; BARLOW, A.; DALY, J. M.; LINCOLN, S. E.; NEWBERG, L. Mapmaker: an interactive computer package for constructing primary genetic linkage maps of experimental and natural populations. Genomics, San Diego, v.1, n. 2, p.174-181, Oct. 1987. LEON, J.; LAWTON, M.A.; RASKIN, I. Hydrogen peroxide stimulates salicylic acid biosynthesis in tobacco. Plant Physiology, Washington, 108, n. 4, p. 1673-1678, 1995. LI, H.-Y.; WANG, T.-Y.; SHI, Y.-S.; FU, J.-J.; SONG, Y.-C.; WANG, G.-Y.; LI, Y. Isolation and characterization of induced genes under drought stress at the flowering stage in maize (Zea mays). DNA Sequence, v. 18, n. 6, p. 445-460, 2007. LIANG, P.; PARDEE, A. B., Differential display of eukaryotic messenger RNA by means of the polymerase chain reaction. Science, Washington, v. 257, n. 14, p. 967-971, Aug. 1992. LITT, M.; LUTY, J. A. A hipervariable microsatellite revealed by in vitro amplification of a dinucleotide repeat within the cardiac muscle actin gene. American Journal of Human Genetic, Chicago, v. 44, p. 398-401, 1989. 41 LIU, B.H.; KNAPP, S.J. GMENDEL: a program for Mendelian segregation and linkage analysis of individual or multiple progeny populations using log-likelihood ratios. Journal of Heredity, Cary, v. 81, n. 5, p. 407-418, Sep. 1990. MELTON, D.A.; KRIEG, P. A.; REBAGLIATI, M. R.; MANIATIS, T.; ZINN, K.; GREEN, M. R. Efficient In Vitro Synthesis of Biologically Active RNA and RNA Hybridization Probes From Plasmids Containing a Bacteriophage SP6 Promoter. Nucleic Acids Research, Oxford, v. 12, n. 18, p. 7035-7056, Sep. 1984. MENDGEN, K.; KRIEG, P. A.; REBAGLIATI, M. R.; MANIATIS, T.; ZINN, K.; GREEN, M. R. Volatiles modulate the development of plant pathogenic rust fungi. Planta, Berlin, v. 224, n. 6, p. 1353-1361, Nov. 2006. MENTABERRY, A. Resistência a hongos fitopatógenos mediante ingeniería genética. 2007. Disponível em . Acesso 01 Dez. 2007. MILES, M. R.; FREDERICK, R. D.; HARTMAN, G. L. Evaluation of soybean germplasm for resistance to Phakopsora pachyrhizi. Plant Health Progress On line, Jan. 2006. doi:10.1094/PHP-2006-0104-01-RS. Disponível em: . Acesso 20 Dez. 2007. MICHELMORE, R.; PARAN, I.; KESSELI, R. V. Identification of markers linked to disease-resistance genes by bulked segregant analysis - a rapid method to detect markers in specific genomic regions by using segregating populations. Proceedings of the National Academy of Sciences of the United States of America, Washington, v. 88, n. 21, p. 9828-9832, Nov. 1991. 42 MONTEROS, M. J.; MISSAOUI, A. M.; PHILLIPS, D. V.; WALKER, D. R.; BOERMA, H. R. Mapping and Confirmation of the 'Hyuuga' Red-Brown Lesion Resistance Gene for Asian Soybean Rust. Crop Science, Madison, 47, p. 829-834, Mar.-Apr. 2007. MORAIS, A.A.C.; SILVA, A.L.. Soja: suas aplicações. São Paulo: MEDSI, 1996. 259 p. MORGANTE, M.; OLIVIERI, A.M. PCR-amplified microsatellites as markers in plant genetics. Plant Journal, Oxford, v. 3, n. 1, p.175-182, Jan. 1993. MORTEL, M. V. D.; RECKNOR, J. C.; GRAHAM, M. A.; NETTLETON, D.; DITTMAN, J. D.; NELSON, R. T.; GODOY, C. V.; ABDELNOOR, R. V.; ALMEIDA, A. M. R.; BAUM, T. J.; WHITHAM, S. A. Distinct biphasic mRNA changes in response to Asian Soybean Rust Infection. Molecular Plant-Microbe Interactions, Sain Paul, v. 20, n. 8, p. 887– 899, Aug. 2007. MULLER, L. Taxonomia e morfologia. In: MIYASAKA, S; MEDINA, J.C. (Ed.). A soja no Brasil. Campinas: ITAL, 1981. p. 65-108. NÓRCIO, L. Grãos de soja resistente à ferrugem asiática devem ser comercializados para a safra de 2010. Disponível em . Acesso em 13 de março de 2008. NÜRNBERGER, T., SCHEEL, D. Signal transmission in the plant imunne response. Trends in Plant Science, London, v. 6, n. 8, p. 372-379, Aug. 2001. NUNTAPUNT, M.; SRISOMBUN, S.; CHUNWONGSE, J. Soybean breeding for rust resistance in Thailand and extent of rust resistant cultivars used. In: World Soybean Research Conference, 7., 2004, Iguasu Falls. Abstracts… Iguasu Falls: Embrapa, 2004, p. 423-430. 43 PASSIANOTTO, A. L. L.; SILVA, D. C. G.; NOGUEIRA, L. M.; SANTOS, J. V. M.; NEPOMUCENO, A. L.; ABDELNOOR, R. V.; YAMANAKA, N. Caracterização de dois genes de resistência à ferrugem asiàtica da soja, Rpp2 e Rpp4, utilizando marcadores moleculares. In: Congresso Brasileiro de Fitopatologia, 40., 2007, Maringá. Anais... Maringá: Sociedade Brasileira de Fitopatologia, 2007, CD ROM. PANSTRUGA, R. Establishing compatibility between plants and obligate biotrophic pathogens. Current Opinion in Plant Biology, v. 6, n. 4, p. 320-326, Aug. 2003. PANTHEE, D. R.; YUAN, J. S.; WRIGHT, D. L.; MAROIS, J. J.; MAILHOT, D.; STEWART, J. R. C.N. Gene expression analysis in soybean in response to the causal agent of Asian soybean rust (Phakopsora pachyrhizi Sydow) in an early growth stage. Functional and Integrative Genomics, Berlin, v. 7, n. 4, p. 291-301, Oct. 2007. PASCHOLATI, S. F.; LEITE, B. Hospedeiro: mecanismos de resistência. In: BERGAMIM FILHO, A.; K