The genomic substrate for adaptive radiation in African cichlid fish

Cichlid fishes are famous for large, diverse and replicated adaptive radiations in the Great Lakes of East Africa. To understand the molecular mechanisms underlying cichlid phenotypic diversity, we sequenced the genomes and transcriptomes of five lineages of African cichlids: the Nile tilapia (Oreochromis niloticus), an ancestral lineage with low diversity; and four members of the East African lineage: Neolamprologus brichardi/pulcher (older radiation, Lake Tanganyika), Metriaclima zebra (recent radiation, Lake Malawi), Pundamilia nyererei (very recent radiation, Lake Victoria), and Astatotilapia burtoni (riverine species around Lake Tanganyika). We found an excess of gene duplications in the East African lineage compared to tilapia and other teleosts, an abundance of non-coding element divergence, accelerated coding sequence evolution, expression divergence associated with transposable element insertions, and regulation by novel microRNAs. In addition, we analysed sequence data from sixty individuals representing six closely related species from Lake Victoria, and show genome-wide diversifying selection on coding and regulatory variants, some of which were recruited from ancient polymorphisms. We conclude that a number of molecular mechanisms shaped East African cichlid genomes, and that amassing of standing variation during periods of relaxed purifying selection may have been important in facilitating subsequent evolutionary diversification.

Association of folate metabolic pathway genes with human spina bifida meningomyelocele

Neural tube defects including spina bifida meningomyelocele (SBMM) are common malformations of the brain and spinal cord, and include all abnormalities resulting from lack of closure of the developing neural tube during embryological development.^ The specific aims of this study were to determine if single nucleotide polymorphic variants (SNPs) in the folate/homocysteine metabolic pathway genes confer a risk for NTD susceptibility within this SBMM population.^ In completion of the first specific aim, two novel SNPs were identified in the FOLR1 gene in Chromosome 11of patients including one in non-coding exon 1 with a C → T transition at nucleotide position 71578317 and another in non-coding exon 3 with a T → G transversion at nucleotide position 71579123. It will be important to determine if these variants are present in the respective parents of these individuals. If they are in fact de novo variants, then these SNPs may be more likely to contribute to the birth defect.^ The second project aim was to analyze genotypes associated with SBMM risk by transmission disequilibrium tests (TDT) and association was detected on several SNPs across the folate metabolic pathway genes in this population. SNPs with significant RC-TDT values were found within the DHFR gene (rs1650723), the MTRR gene (rs327592), the FOLR2 gene (rs13908), four tightly linked variants in the FOLR3 gene (rs7925545, rs7926875, rs7926987, rs7926360) and a variant in the SLC19A1 gene (rs1888530). The product of each of these genes performs a vital function in the folate metabolic pathway. It is conceivable, therefore, that if the individual SNP or SNPs can be proven to perturb the function in some way that they may be involved in the disruption of folate metabolism and in the resulting birth defect. Validating the results of this study in other independent populations will further strengthen the evidence that dysfunction of folate enzymes and receptors may confer SBMM risk in humans. ^

Arabidopsis thaliana bZIP44: a transcription factor affecting seed germination and expression of the mannanase-encoding gene AtMAN7.

Endo-β-mannanases (MAN; EC. 3.2.1.78) catalyze the cleavage of β1[RIGHTWARDS ARROW]4 bonds in mannan polymers and have been associated with the process of weakening the tissues surrounding the embryo during seed germination. In germinating Arabidopsis thaliana seeds, the most highly expressed MAN gene is AtMAN7 and its transcripts are restricted to the micropylar endosperm and to the radicle tip just before radicle emergence. Mutants with a T-DNA insertion in AtMAN7 have a slower germination than the wild type. To gain insight into the transcriptional regulation of the AtMAN7 gene, a bioinformatic search for conserved non-coding cis-elements (phylogenetic shadowing) within the Brassicaceae MAN7 gene promoters has been done, and these conserved motifs have been used as bait to look for their interacting transcription factors (TFs), using as a prey an arrayed yeast library from A. thaliana. The basic-leucine zipper TF AtbZIP44, but not the closely related AtbZIP11, has thus been identified and its transcriptional activation upon AtMAN7 has been validated at the molecular level. In the knock-out lines of AtbZIP44, not only is the expression of the AtMAN7 gene drastically reduced, but these mutants have a significantly slower germination than the wild type, being affected in the two phases of the germination process, both in the rupture of the seed coat and in the breakage of the micropylar endosperm cell walls. In the over-expression lines the opposite phenotype is observed.

Transcription factors of the bzip factors of the BZIP family affecting arabidopsis thaliana germination sensu stricto

During seed germination, the endosperm cell walls (CWs) suffer an important weakening process mainly driven by hydrolytic enzymes, such are endo-?- mannanases (MAN; EC. 3.2.1.78) that catalyze the cleavage of ?1?4 bonds in the mannan-polymers. In Arabidopsis thaliana seeds, endo-?-mannanase activity increases during seed imbibition, decreasing after radicle emergence1. AtMAN7 is the most highly expressed MAN gene in seeds upon germination and their transcripts are restricted to the micropylar endosperm and to the radicle tip just before radicle emergence. Mutants with a T-DNA insertion in this gene (K.O. MAN7) have a slower germination rate than the wild type (t50=34 h versus t50=25 h). To gain insight into the transcriptional regulation of the AtMAN7 gene, a bioinformatic search for conserved non-coding cis-elements (phylogenetic shadowing) within the Brassicaceae orthologous MAN7 gene promoters has been done and these conserved motives have been used as baits to look for their interacting transcription factors (TFs), using as a prey an arrayed yeast library of circa 1,200 TFs from A. thaliana. The basic leucine zipper AtbZIP44, but not its closely related ortholog AtbZIP11, has been thus identified and its regulatory function upon AtMAN7 during seed germination validated by different molecular and physiological techniques, such are RT-qPCR analyses, mRNA Fluorescence in situ Hybridization (FISH) experiments, and by the establishment of the germination kinetics of both over-expression (oex) lines and TDNA insertion mutants in AtbZIP44. The transcriptional combinatorial network through which AtbZIP44 regulates AtMAN7 gene expression during seed germination has been further explored through protein-protein interactions between AtbZIP44 and other bZIP members. In such a way, AtbZIP9 has been identified by yeast two-hybrid experiments and its physiological implication in the control of AtMAN7 expression similarly established.

Regulation of c-maf gene expression by Pax6 in cultured cells

c-Maf is a bZip transcription factor expressed in developmental and cellular differentiation processes. Recently, a c-maf knockout mouse model, showing abnormal lens development, has been reported. In order to study the regulation mechanisms of c-maf gene expression during the differentiation process we have cloned and functionally characterized the rat c-maf (maf-2) gene. The rat c-maf gene is an intronless gene, covering a length of 3.5 kb. Transient transfection analysis of the 5′-flanking region of the c-maf gene using luciferase as the reporter gene shows that Pax6, a master transcription factor for lens development, strongly activates the c-maf promoter construct. Endogenous c-maf is also activated by the Pax6 expression vector. Electrophoresis mobility shift assay and DNase I footprinting analysis show that at least three Pax6-binding sites are located in the 5′-flanking and 5′-non-coding regions of the rat c-maf gene. The c-maf gene was also markedly activated by its own product, c-Maf, through the MARE (Maf recognition element), suggesting that a positive autoregulatory mechanism controls this gene. In situ hybridization histochemical detection of Pax6 and c-Maf in the E14 lens showed that both mRNAs are expressed in the lens equator where lens epithelial cells are differentiating to lens fiber cells. These results suggest that a Pax6/c-Maf transcription factor cascade is working in lens development.

Ruolo di Clusterina nella progressione del tumore prostatico

Clusterina (CLU) è una proteina ubiquitaria, presente nella maggior parte dei fluidi corporei e implicata in svariati processi fisiologici. Dalla sua scoperta fino ad oggi, CLU è risultata essere una proteina enigmatica, la cui funzione non è ancora stata compresa appieno. Il gene codifica per 3 varianti trascrizionali identificate nel database NCBI con i codici: NM_001831 (CLU 1 in questo lavoro di tesi), NR_038335 (CLU 2 in questo lavoro di tesi) e NR_045494 (CLU 3 in questo lavoro di tesi). Tutte le varianti sono trascritte come pre-mRNA contenenti 9 esoni e 8 introni e si differenziano per l’esone 1, la cui sequenza è unica e caratteristica di ogni variante. Sebbene in NCBI sia annotato che le varianti CLU 2 e CLU 3 non sono codificanti, tramite analisi bioinformatica è stato predetto che da tutti e tre i trascritti possono generarsi proteine di differente lunghezza e localizzazione cellulare. Tra tutte le forme proteiche ipotizzate, l’unica a essere stata isolata e sequenziata è quella tradotta dall’AUG presente sull’esone 2 che dà origine a una proteina di 449 aminoacidi. Il processo di maturazione prevede la formazione di un precursore citoplasmatico (psCLU) che subisce modificazioni post-traduzionali tra cui formazione di ponti disolfuro, glicosilazioni, taglio in due catene denominate β e α prima di essere secreta come eterodimero βα (sCLU) nell’ambiente extracellulare, dove esercita la sua funzione di chaperone ATP-indipendente. Oltre alla forma extracellulare, è possibile osservare una forma intracellulare con localizzazione citosolica la cui funzione non è stata ancora completamente chiarita. Questo lavoro di tesi si è prefissato lo scopo di incrementare le conoscenze in merito ai trascritti CLU 1 e CLU 2 e alla loro regolazione, oltre ad approfondire il ruolo della forma citosolica della proteina in relazione al signaling di NF-kB che svolge un ruolo importante nel processo di sviluppo e metastatizzazione del tumore. Nella prima parte, uno screening di differenti linee cellulari, quali cellule epiteliali di prostata e di mammella, sia normali sia tumorali, fibroblasti di origine polmonare e linfociti di tumore non-Hodgkin, ha permesso di caratterizzare i trascritti CLU 1 e CLU 2. Dall’analisi è emerso che la sequenza di CLU 1 è più corta al 5’ rispetto a quella depositata in NCBI con l’identificativo NM_001831 e il primo AUG disponibile per l’inizio della traduzione è localizzato sull’esone 2. È stato dimostrato che CLU 2, al contrario di quanto riportato in NCBI, è tradotto in proteina a partire dall’AUG presente sull’esone 2, allo stesso modo in cui viene tradotto CLU 1. Inoltre, è stato osservato che i livelli d’espressione dei trascritti variano notevolmente tra le diverse linee cellulari e nelle cellule epiteliali CLU 2 è espressa sempre a bassi livelli. In queste cellule, l’espressione di CLU 2 è silenziata per via epigenetica e la somministrazione di farmaci capaci di rendere la cromatina più accessibile, quali tricostatina A e 5-aza-2’-deossicitidina, è in grado di incrementarne l’espressione. Nella seconda parte, un’analisi bioinformatica seguita da saggi di attività in vitro in cellule epiteliali prostatiche trattate con farmaci epigenetici, hanno permesso di identificare, per la prima volta in uomo, una seconda regione regolatrice denominata P2, capace di controllare l’espressione di CLU 2. Rispetto a P1, il classico promotore di CLU già ampiamente studiato da altri gruppi di ricerca, P2 è un promotore debole, privo di TATA box, che nelle cellule epiteliali prostatiche è silente in condizioni basali e la cui attività incrementa in seguito alla somministrazione di farmaci epigenetici capaci di alterare le modificazioni post-traduzionali delle code istoniche nell’intorno di P2. Ne consegue un rilassamento della cromatina e un successivo aumento di trascrizione di CLU 2. La presenza di un’isola CpG differentemente metilata nell’intorno di P1 spiegherebbe, almeno in parte, i differenti livelli di espressione di CLU che si osservano tra le diverse linee cellulari. Nella terza parte, l’analisi del pathway di NF-kB in un modello sperimentale di tumore prostatico in cui CLU è stata silenziata o sovraespressa, ha permesso di capire come la forma citosolica di CLU abbia un ruolo inibitorio nei confronti dell’attività del fattore trascrizionale NF-kB. CLU inibisce la fosforilazione e l’attivazione di p65, il membro più rappresentativo della famiglia NF-kB, con conseguente riduzione della trascrizione di alcuni geni da esso regolati e coinvolti nel rimodellamento della matrice extracellulare, quali l’urochinasi attivatrice del plasminogeno, la catepsina B e la metallo proteinasi 9. È stato dimostrato che tale inibizione non è dovuta a un’interazione fisica diretta tra CLU e p65, per cui si suppone che CLU interagisca con uno dei componenti più a monte della via di segnalazione responsabile della fosforilazione ed attivazione di p65.

Análise transcriptômica de miRNAs em pacientes sob dupla terapia de antiagregação

A terapia antiagregante é comumente indicada na prevenção e tratamento de doenças cardiovasculares. A dupla antiagregação com clopidrogrel e ácido acetilsalicílico (AAS) tem sido frequentemente adotada em pacientes com Doença Arterial Coronariana (DAC), mas apresenta ineficácia em uma parcela significativa da população com genótipo de respondedores. Essa falha terapêutica nos leva a questionar se outros mecanismos moleculares podem estar influenciando na resposta a esses fármacos. Recentes estudos sugerem que pequenas sequências de RNA não codificantes denominadas microRNAs (miRNAs) podem estar fortemente relacionadas com resposta ao tratamento fármaco-terapêutico, controlando as proteínas envolvidas na farmacocinética e farmacodinâmica. Entretanto, os principais miRNAs que atuam na dinâmica da resposta medicamentosa ainda não foram bem definidos. O objetivo deste estudo foi avaliar o perfil de miRNAs no sangue total periférico, procurando melhor esclarecer os mecanismos envolvidos na resposta aos antiagregantes plaquetários AAS e clopidogrel. Para isso, selecionou-se pacientes com DAC, os quais apresentavam diferentes respostas à dupla terapia de antiagregação determinadas pelo teste de agregação plaquetária. Baseados nos fenótipos, os perfis de expressão de miRNAs foram comparados entre os valores da taxa de agregação categorizados em tercis (T) de resposta. O grupo T1 foi constituído de pacientes respondedores, o T2 de respondedores intermediários e o T3 de não respondedores. Os perfis de miRNAs foram obtidos após sequenciamento de última geração e os dados obtidos foram analisados pelo pacote Deseq2. Os resultados mostraram 18 miRNAs diferentemente expressos entre os dois tercis extremos. Dentre esses miRNAs, 10 deles apresentaram importantes alvos relacionados com vias de ativação e agregação plaquetária quando analisados pelo software Ingenuity®. Dos 10 miRNAs, 4 deles, os quais apresentaram-se menos expressos no sequenciamento, demonstraram os mesmos perfis de expressão quando analisados pela reação em cadeia pela polimerase quantitativa (qPCR): hsa-miR-423-3p, hsa-miR-744-5p, hsa-miR- 30a-5p e hsa-let-7g-5p. A partir das análises de predição de alvos, pôde-se observar que os quatro miRNAs, quando menos expressos simultaneamente, predizem ativação da agregação plaquetária. Além disso, os miRNAs hsa-miR- 423-5p, hsa-miR-744-5p e hsa-let-7g-5p mostraram correlação com o perfil lipídico dos pacientes que, por sua vez, apresentou influência nos valores de agregação compreendidos no T3 de resposta a ambos os medicamentos. Sendo assim, conclui-se que maiores taxas de agregação plaquetária podem estar indiretamente relacionadas com os padrões de expressão de hsa-miR- 423-3p, hsa-miR-744-5p e hsa-let-7g-5p. Sugere-se que a avaliação do perfil de expressão destes 3 miRNAs no sangue periférico de pacientes com DAC possa predizer resposta terapêutica inadequada ao AAS e ao clopidogrel

The prion protein gene: Identifying regulatory signals using marsupial sequence

The function of the prion protein gene (PRNP) and its normal product PrPC is elusive. We used comparative genomics as a strategy to understand the normal function of PRNP. As the reliability of comparisons increases with the number of species and increased evolutionary distance, we isolated and sequenced a 66.5 kb BAC containing the PRNP gene from a distantly related mammal, the model Australian marsupial Macropus eugenii (tammar wallaby). Marsupials are separated from eutherians such as human and mouse by roughly 180 million years of independent evolution. We found that tammar PRNP, like human PRNP, has two exons. Prion proteins encoded by the tammar wallaby and a distantly related marsupial, Monodelphis domestica (Brazilian opossum) PRNP contain proximal PrP repeats with a distinct, marsupial-specific composition and a variable number. Comparisons of tammar wallaby PRNP with PRNPs from human, mouse, bovine and ovine allowed us to identify non-coding gene regions conserved across the marsupial-eutherian evolutionary distance, which are candidates for regulatory regions. In the PRNP 3' UTR we found a conserved signal for nuclear-specific polyadenylation and the putative cytoplasmic polyadenylation element (CPE), indicating that post-transcriptional control of PRNP mRNA activity is important. Phylogenetic footprinting revealed conserved potential binding sites for the MZF-1 transcription factor in both upstream promoter and intron/intron 1, and for the MEF2, MyTI, Oct-1 and NFAT transcription factors in the intron(s). The presence of a conserved NFAT-binding site and CPE indicates involvement of PrPC in signal transduction and synaptic plasticity. (c) 2004 Elsevier B.V. All rights reserved.

Genomic organization of human arylamine N-acetyltransferase Type I reveals alternative promoters that generate different 5'-UTR splice variants with altered translational activities

In humans, a polymorphic gene encodes the drug-metabolizing enzyme NATI (arylamine N-acetyltransferase Type 1), which is widely expressed throughout the body. While the protein-coding region of NATI is contained within a single exon, examination of the human EST (expressed sequence tag) database at the NCBI revealed the presence of nine separate exons, eight of which were located in the 5'non-coding region of NATI. Differential splicing produced at least eight unique mRNA isoforms that could be grouped according to the location of the first exon, which suggested that NATI expression occurs from three alternative promoters. Using RT (reverse transcriptase)-PCR, we identified one major transcript in various epithelial cells derived from different tissues. In contrast, multiple transcripts were observed in blood-derived cell lines (CEM, THP-1 and Jurkat), with a novel variant, not identified in the EST database, found in CEM cells only. The major splice variant increased gene expression 9-11-fold in a luciferase reporter assay, while the other isoforrns were similar or slightly greater than the control. We examined the upstream region of the most active splice variant in a promoter-reporter assay, and isolated a 257 bp sequence that produced maximal promoter activity. This sequence lacked a TATA box, but contained a consensus Sp1 site and a CAAT box, as well as several other putative transcription-factor-binding sites. Cell-specific expression of the different NATI transcripts may contribute to the variation in NATI activity in vivo.

RNAdb - a comprehensive mammalian noncoding RNA database

In recent years, there have been increasing numbers of transcripts identified that do not encode proteins, many of which are developmentally regulated and appear to have regulatory functions. Here, we describe the construction of a comprehensive mammalian noncoding RNA database (RNAdb) which contains over 800 unique experimentally studied noncoding RNAs (ncRNAs), including many associated with diseases and/or developmental processes. The database is available at http://research.imb.uq. edu.au/RNAdb and is searchable by many criteria. It includes microRNAs and snoRNAs, but not infrastructural RNAs, such as rRNAs and tRNAs, which are catalogued elsewhere. The database also includes over 1100 putative antisense ncRNAs and almost 20000 putative ncRNAs identified in high-quality murine and human cDNA libraries, with more to be added in the near future. Many of these RNAs are large, and many are spliced, some alternatively. The database will be useful as a foundation for the emerging field of RNomics and the characterization of the roles of ncRNAs in mammalian gene expression and regulation.

Teaching principles of genetics through SNP analysis

An understanding of inheritance requires comprehension of genetic processes at all levels, from molecules to populations. Frequently genetics courses are separated into molecular and organismal genetics and students may fail to see the relationships between them. This is particularly true with human genetics, because of the difficulties in designing experimental approaches which are consistent with ethical restrictions, student abilities and background knowledge, and available time and materials. During 2005 we used analysis of single nucleotide polymorphisms (SNPs) in two genetic regions to enhance student learning and provide a practical experience in human genetics. Students scanned databases to discover SNPs in a gene of interest, used software to design PCR primers and a restriction enzyme based assay for the alleles, and carried out an analysis of the SNP on anonymous individual and family DNAs. The project occupied eight to ten hours per week for one semester, with some time spent in the laboratory and some spent in database searching, reading and writing the report. In completing their projects, students acquired a knowledge of Mendel’s first law (through looking at inheritance patterns), Mendel’s second law and the exceptions (the concepts of linkage and linkage disequilibrium), DNA structure (primer design and restriction enzyme analysis) and function (SNPs in coding and non-coding regions), population genetics and the statistical analysis of allele frequencies, genomics, bioinformatics and the ethical issues associated with the use of human samples. They also developed skills in presentation of results by publication and conference participation. Deficiencies in their understanding (for example of inheritance patterns, gene structure, statistical approaches and report writing) were detected and guidance given during the project. SNP analysis was found to be a powerful approach to enhance and integrate student understanding of genetic concepts.

Novel expression of microRNAs in serum samples of Iraqi breast cancer women

Although a lot of hard work against cancer to reduces its spread but it still continues to kill with abandon. The need for a biomarker for cancer early detection becomes the most mind concentrated scientists. MicroRNAs the tiny non coding RNA molecules opened new path for the scientists to determine the cancer in its early stages. Expression of microRNAs profiles has been investigated to be involved in cancer development. Here we determined the expression of microRNAs in serum of Iraqi healthy volunteers and other women diagnosed with breast cancer. MicroRNAs expression has been determined by using real time qPCR and delta method has been used. Four of thirteen microRNAs were shown to be expressed in serum of Iraqi breast cancer women. Let-7a and miR-21 were shown to be significantly over expressed in serum of breast cancer compared with healthy serum volunteers (P= 0.022 and 0.026) respectively. While miR-26b and miR-429 found to be significantly down expressed in serum of breast cancer women (P= 0.0034 and 0.031) respectively. The result concluded that these expressed microRNAs in serum of breast cancer women could be used as a first indicator of breast cancer occurrence.

The handedness-associated PCSK6 locus spans an intronic promoter regulating novel transcripts

We recently reported the association of the PCSK6 gene with handedness through a quantitative genome-wide association study (GWAS; P < 0.5 × 10(-8)) for a relative hand skill measure in individuals with dyslexia. PCSK6 activates Nodal, a morphogen involved in regulating left-right body axis determination. Therefore, the GWAS data suggest that the biology underlying the patterning of structural asymmetries may also contribute to behavioural laterality, e.g. handedness. The association is further supported by an independent study reporting a variable number tandem repeat (VNTR) within the same PCSK6 locus to be associated with degree of handedness in a general population cohort. Here, we have conducted a functional analysis of the PCSK6 locus combining further genetic analysis, in silico predictions and molecular assays. We have shown that the previous GWAS signal was not tagging a VNTR effect, suggesting that the two markers have independent effects. We demonstrated experimentally that one of the top GWAS-associated markers, rs11855145, directly alters the binding site for a nuclear factor. Furthermore, we have shown that the predicted regulatory region adjacent to rs11855415 acts as a bidirectional promoter controlling the expression of novel RNA transcripts. These include both an antisense long non-coding RNA (lncRNA) and a short PCSK6 isoform predicted to be coding. This is the first molecular characterization of a handedness-associated locus that supports the role of common variants in non-coding sequences in influencing complex phenotypes through gene expression regulation.

The isolation, characterization, and application of WRKY genes as useful molecular markers in tropical trees

The improvement of tropical tree crops using conventional breeding methods faces challenges due to the length of time involved. Thus, like most crops, there is an effort to utilize molecular genetic markers in breeding programs to select for desirable agronomic traits. Known as marker assisted breeding or marker assisted selection, genetic markers associated with a phenotype of interest are used to screen and select material reducing the time necessary to evaluate candidates. As the focus of this research was improving disease resistance in tropical trees, the usefulness of the WRKY gene superfamily was investigated as candidates for generating useful molecular genetic markers. WRKY genes encode plant-specific transcriptional factors associated with regulating plants' responses to both biotic and abiotic stress. ^ One pair of degenerate primers amplified 48 WRKY gene fragments from three taxonomically distinct, economically important, tropical tree crop species: 18 from Theobroma cacao L., 21 from Cocos nucifera L. and 9 from Persea americana Mill. Several loci from each species were polymorphic because of single nucleotide substitutions present within a putative non-coding region of the loci. Capillary array electrophoresis-single strand conformational polymorphism (CAE-SSCP) mapped four WRKY loci onto a genetic linkage map of a T. cacao F2 population segregating for resistance to witches' broom disease. Additionally, PCR primers specific for four T. cacao loci successfully amplified WRKY loci from 15 members of the Byttneriae tribe. A method was devised to allow the reliable discrimination of alleles by CAE-SSCP using only the mobility assigned to the sample peaks. Once this method was validated, the diversity of three WRKY loci was evaluated in a germplasm collection of T. cacao . One locus displayed high diversity in the collection, with at least 18 alleles detected from mobility differences of the product peaks. The number of WRKY loci available within the genome, ease of isolation by degenerate PCR, codominant segregation demonstrated in the F2 population, and usefulness for screening germplasm collections and closely related wild species demonstrates that the WRKY superfamily of genes are excellent candidates for developing a number of genetic molecular markers for breeding purposes in tropical trees. ^