Genomic DNA sequence comparison between two inbred soybean cyst nematode biotypes facilitated by massively parallel 454 micro-bead sequencing

Comparative Transcriptomic Analysis of Soybean Cyst Nematode Inbred Populations Non-adapted or Adapted on Soybean rhg1-a/Rhg4-Mediated Resistance

Soybean cyst nematode (SCN, Heterodera glycines) is most effectively managed through planting resistant soybean cultivars, but the repeated use of the same resistance sources has led to a widespread emergence of virulent SCN populations that can overcome soybean resistance. Resistance to SCN HG type 0 (Race 3) in soybean cultivar Forrest is mediated by an epistatic interaction between the soybean resistance genes rhg1-a and Rhg4. We previously developed two SCN inbred populations by mass-selecting SCN HG type 0 (Race 3) on susceptible and resistant recombinant inbred lines, derived from a cross between Forrest and the SCN-susceptible cultivar Essex, which differ for Rhg4. To identify SCN genes potentially involved in overcoming rhg1-a/Rhg4-mediated resistance, we conducted RNA sequencing on early parasitic juveniles of these two SCN inbred populations infecting their respective hosts, only to discover a handful of differentially expressed genes (DEGs). However, in a comparison with early parasitic juveniles of an avirulent SCN inbred population infecting a resistant host, we discovered 59 and 171 DEGs uniquely up- or downregulated in virulent parasitic juveniles adapted on the resistant host. Interestingly, the proteins coded by these 59 DEGs included vitamin B-associated proteins (reduced folate carrier, biotin synthase, and thiamine transporter) and nematode effectors known to play roles in plant defense suppression, suggesting that virulent SCN may exert a heightened transcriptional response to cope with enhanced plant defenses and an altered nutritional status of a resistant soybean host. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Genome scans for selection signatures identify candidate virulence genes for adaptation of the soybean cyst nematode to host resistance

Plant pathogens are constantly under selection pressure for host resistance adaptation. Soybean cyst nematode (SCN, Heterodera glycines) is a major pest of soybean primarily managed through resistant cultivars; however, SCN populations have evolved virulence in response to selection pressures driven by repeated monoculture of the same genetic resistance. Resistance to SCN is mediated by multiple epistatic interactions between Rhg (for resistance to H. glycines) genes. However, the identity of SCN virulence genes that confer the ability to overcome resistance remains unknown. To identify candidate genomic regions showing signatures of selection for increased virulence, we conducted whole genome resequencing of pooled individuals (Pool-Seq) from two pairs of SCN populations adapted on soybeans with Peking-type (rhg1-a, rhg2, and Rhg4) resistance. Population differentiation and principal component analysis-based approaches identified approximately 0.72-0.79 million SNPs, the frequency of which showed potential selection signatures across multiple genomic regions. Chromosomes 3 and 6 between population pairs showed the greatest density of outlier SNPs with high population differentiation. Conducting multiple outlier detection tests to identify overlapping SNPs resulted in a total of 966 significantly differentiated SNPs, of which 285 exon SNPs were mapped to 97 genes. Of these, six genes encoded members of known stylet-secreted effector protein families potentially involved in host defence modulation including venom-allergen-like, annexin, glutathione synthetase, SPRYSEC, chitinase, and CLE effector proteins. Further functional analysis of identified candidate genes will provide new insights into the genetic mechanisms by which SCN overcomes soybean resistance and inform the development of molecular markers for rapidly screening the virulence profile of an SCN-infested field.

Comparative genomics of two inbred lines of the potato cyst nematode Globodera rostochiensis reveals disparate effector family-specific diversification patterns

BACKGROUND

Potato cyst nematodes belong to the most harmful pathogens in potato, and durable management of these parasites largely depends on host-plant resistances. These resistances are pathotype specific. The current Globodera rostochiensis pathotype scheme that defines five pathotypes (Ro1 - Ro5) is both fundamentally and practically of limited value. Hence, resistant potato varieties are used worldwide in a poorly informed manner.

RESULTS

We generated two novel reference genomes of G. rostochiensis inbred lines derived from a Ro1 and a Ro5 population. These genome sequences comprise 173 and 189 scaffolds respectively, marking a ≈ 24-fold reduction in fragmentation as compared to the current reference genome. We provide copy number variations for 19 effector families. Four dorsal gland effector families were investigated in more detail. SPRYSECs, known to be implicated in plant defence suppression, constitute by far the most diversified family studied herein with 60 and 99 variants in Ro1 and Ro5 distributed over 18 and 26 scaffolds. In contrast, CLEs, effectors involved in feeding site induction, show strong physical clustering. The 10 and 16 variants cluster on respectively 2 and 1 scaffolds. Given that pathotypes are defined by their effectoromes, we pinpoint the disparate nature of the contributing effector families in terms of sequence diversification and loss and gain of variants.

CONCLUSIONS

Two novel reference genomes allow for nearly complete inventories of effector diversification and physical organisation within and between pathotypes. Combined with insights we provide on effector family-specific diversification patterns, this constitutes a basis for an effectorome-based virulence scheme for this notorious pathogen.

Recent Advances in Population Genomics of Plant-Parasitic Nematodes

Plant-parasitic nematodes are a costly burden of crop production. Ubiquitous in nature, phytoparasitic nematodes are associated with nearly every important agricultural crop and represent a significant constraint on global food security. Population genetics is a key discipline in plant nematology to understand aspects of the life strategies of these parasites, in particular their modes of reproduction, geographic origins, evolutionary histories, and dispersion abilities. Advances in high-throughput sequencing technologies have enabled a recent but active effort in genomic analyses of plant-parasitic nematodes. Such genomic approaches applied to multiple populations are providing new insights into the molecular and evolutionary processes that underpin the establishment of these nematodes and into a better understanding of the genetic and mechanistic basis of their pathogenicity and adaptation to their host plants. In this review, we attempt to update information about genome resources and genotyping techniques useful for nematologists who are thinking about initiating population genomics or genome sequencing projects. This review is intended also to foster the development of population genomics in plant-parasitic nematodes through highlighting recent publications that illustrate the potential for this approach to identify novel molecular markers or genes of interest and improve our knowledge of the genome variability, pathogenicity, and evolutionary potential of plant-parasitic nematodes.

Soybean Resistance to the Soybean Cyst Nematode Heterodera glycines: An Update

The soybean cyst nematode (SCN), Heterodera glycines, remains a serious threat to soybean production throughout the world. A lack of genetic diversity in resistant soybean cultivars has led to a widespread shift toward virulence in SCN populations, leaving farmers with few proven options other than nonhost rotation to manage this nematode. Recent advances in our understanding of the genes controlling resistance to the nematode have led to improved molecular markers, which are, in turn, increasing the efficiency and precision of the breeding pipeline. A better understanding of the molecular and biochemical basis of SCN resistance and nematode virulence will provide information useful for the development of a long-term strategic plan for diversification and the deployment of cultivars that protect current sources of natural resistance while identifying new targets for engineering novel resistance.

Deciphering the molecular mechanisms behind cellulase production in Trichoderma reesei, the hyper-cellulolytic filamentous fungus

The filamentous fungus Trichoderma reesei is a potent cellulase producer and the best-studied cellulolytic fungus. A lot of investigations not only on glycoside hydrolases produced by T. reesei, but also on the machinery controlling gene expression of these enzyme have made this fungus a model organism for cellulolytic fungi. We have investigated the T. reesei strain including mutants developed in Japan in detail to understand the molecular mechanisms that control the cellulase gene expression, the biochemical and morphological aspects that could favor this phenotype, and have attempted to generate novel strains that may be appropriate for industrial use. Subsequently, we developed recombinant strains by combination of these insights and the heterologous-efficient saccharifing enzymes. Resulting enzyme preparations were highly effective for saccharification of various biomass. In this review, we present some of the salient findings from the recent biochemical, morphological, and molecular analyses of this remarkable cellulase hyper-producing fungus.

Characterization of the reniform nematode genome by shotgun sequencing

The reniform nematode (RN), a major agricultural pest particularly on cotton in the United States, is among the major plant-parasitic nematodes for which limited genomic information exists. In this study, over 380 Mb of sequence data were generated from pooled DNA of four adult female RNs and assembled into 67,317 contigs, including 25,904 (38.5%) predicted coding contigs and 41,413 (61.5%) noncoding contigs. Most of the characterized repeats were of low complexity (88.9%), and 0.9% of the contigs matched with 53.2% of GenBank ESTs. The most frequent Gene Ontology (GO) terms for molecular function and biological process were protein binding (32%) and embryonic development (20%). Further analysis showed that 741 (1.1%), 94 (0.1%), and 169 (0.25%) RN genomic contigs matched with 1328 (13.9%), 1480 (5.4%), and 1330 (7.4%) supercontigs of Meloidogyne incognita, Brugia malayi, and Pristionchus pacificus, respectively. Chromosome 5 of Caenorhabditis elegans had the highest number of hits to the RN contigs. Seven putative detoxification genes and three carbohydrate-active enzymes (CAZymes) involved in cell wall degradation were studied in more detail. Additionally, kinases, G protein-coupled receptors, and neuropeptides functioning in physiological, developmental, and regulatory processes were identified in the RN genome.

Comparative transcriptome analysis of two races of Heterodera glycines at different developmental stages

The soybean cyst nematode, Heterodera glycines, is an important pest of soybeans. Although resistance is available against this nematode, selection for virulent races can occur, allowing the nematode to overcome the resistance of cultivars. There are abundant field populations, however, little is known about their genetic diversity. In order to elucidate the differences between races, we investigated the transcriptional diversity within race 3 and race 4 inbred lines during their compatible interactions with the soybean host Zhonghuang 13. Six different race-enriched cDNA libraries were constructed with limited nematode samples collected from the three sedentary stages, parasitic J2, J3 and J4 female, respectively. Among 689 putative race-enriched genes isolated from the six libraries with functional annotations, 92 were validated by quantitative RT-PCR (qRT-PCR), including eight putative effector encoding genes. Further race-enriched genes were validated within race 3 and race 4 during development in soybean roots. Gene Ontology (GO) analysis of all the race-enriched genes at J3 and J4 female stages showed that most of them functioned in metabolic processes. Relative transcript level analysis of 13 selected race-enriched genes at four developmental stages showed that the differences in their expression abundance took place at either one or more developmental stages. This is the first investigation into the transcript diversity of H. glycines races throughout their sedentary stages, increasing the understanding of the genetic diversity of H. glycines.

The expression of a naturally occurring, truncated allele of an α-SNAP gene suppresses plant parasitic nematode infection

Transcriptional mapping experiments of the major soybean cyst nematode resistance locus, rhg1, identified expression of the vesicular transport machinery component, α soluble NSF attachment protein (α-SNAP), occurring during defense. Sequencing the α-SNAP coding regions from the resistant genotypes G. max ([Peking/PI 548402]) and G. max ([PI 437654]) revealed they are identical, but differ from the susceptible G. max ([Williams 82/PI 518671]) by the presence of several single nucleotide polymorphisms. Using G. max ([Williams 82/PI 518671]) as a reference, a G → T(2,822) transversion in the genomic DNA sequence at a functional splice site of the α-SNAP([Peking/PI 548402]) allele produced an additional 17 nucleotides of mRNA sequence that contains an in-frame stop codon caused by a downstream G → A(2,832) transition. The G. max ([Peking/PI 548402]) genotype has cell wall appositions (CWAs), structures identified as forming as part of a defense response by the activity of the vesicular transport machinery. In contrast, the 17 nt α-SNAP([Peking/PI 548402]) mRNA motif is not found in G. max ([PI 88788]) that exhibits defense to H. glycines, but lack CWAs. The α-SNAP([PI 88788]) promoter contains sequence elements that are nearly identical to the α-SNAP([Peking/PI 548402]) allele, but differs from the G. max ([Williams 82/PI 518671]) ortholog. Overexpressing the α-SNAP([Peking/PI 548402]) allele in the susceptible G. max ([Williams 82/PI 518671]) genotype suppressed H. glycines infection. The experiments indicate a role for the vesicular transport machinery during infection of soybean by the soybean cyst nematode. However, increased GmEREBP1, PR1, PR2, PR5 gene activity but suppressed PR3 expression accompanied the overexpression of the α-SNAP([Peking/PI 548402]) allele prior to infection.

Transcriptomics of a giant freshwater prawn (Macrobrachium rosenbergii): de novo assembly, annotation and marker discovery

BACKGROUND

Giant freshwater prawn (Macrobrachium rosenbergii or GFP), is the most economically important freshwater crustacean species. However, as little is known about its genome, 454 pyrosequencing of cDNA was undertaken to characterise its transcriptome and identify genes important for growth.

METHODOLOGY AND PRINCIPAL FINDINGS

A collection of 787,731 sequence reads (244.37 Mb) obtained from 454 pyrosequencing analysis of cDNA prepared from muscle, ovary and testis tissues taken from 18 adult prawns was assembled into 123,534 expressed sequence tags (ESTs). Of these, 46% of the 8,411 contigs and 19% of 115,123 singletons possessed high similarity to sequences in the GenBank non-redundant database, with most significant (E value < 1e(-5)) contig (80%) and singleton (84%) matches occurring with crustacean and insect sequences. KEGG analysis of the contig open reading frames identified putative members of several biological pathways potentially important for growth. The top InterProScan domains detected included RNA recognition motifs, serine/threonine-protein kinase-like domains, actin-like families, and zinc finger domains. Transcripts derived from genes such as actin, myosin heavy and light chain, tropomyosin and troponin with fundamental roles in muscle development and construction were abundant. Amongst the contigs, 834 single nucleotide polymorphisms, 1198 indels and 658 simple sequence repeats motifs were also identified.

CONCLUSIONS

The M. rosenbergii transcriptome data reported here should provide an invaluable resource for improving our understanding of this species' genome structure and biology. The data will also instruct future functional studies to manipulate or select for genes influencing growth that should find practical applications in aquaculture breeding programs.

Identification of ovule transcripts from the Apospory-Specific Genomic Region (ASGR)-carrier chromosome

Background

Apomixis, asexual seed production in plants, holds great potential for agriculture as a means to fix hybrid vigor. Apospory is a form of apomixis where the embryo develops from an unreduced egg that is derived from a somatic nucellar cell, the aposporous initial, via mitosis. Understanding the molecular mechanism regulating aposporous initial specification will be a critical step toward elucidation of apomixis and also provide insight into developmental regulation and downstream signaling that results in apomixis. To discover candidate transcripts for regulating aposporous initial specification in P. squamulatum, we compared two transcriptomes derived from microdissected ovules at the stage of aposporous initial formation between the apomictic donor parent, P. squamulatum (accession PS26), and an apomictic derived backcross 8 (BC₈) line containing only the Apospory-Specific Genomic Region (ASGR)-carrier chromosome from P. squamulatum. Toward this end, two transcriptomes derived from ovules of an apomictic donor parent and its apomictic backcross derivative at the stage of apospory initiation, were sequenced using 454-FLX technology.

Results

Using 454-FLX technology, we generated 332,567 reads with an average read length of 147 base pairs (bp) for the PS26 ovule transcriptome library and 363,637 reads with an average read length of 142 bp for the BC₈ ovule transcriptome library. A total of 33,977 contigs from the PS26 ovule transcriptome library and 26,576 contigs from the BC₈ ovule transcriptome library were assembled using the Multifunctional Inertial Reference Assembly program. Using stringent in silico parameters, 61 transcripts were predicted to map to the ASGR-carrier chromosome, of which 49 transcripts were verified as ASGR-carrier chromosome specific. One of the alien expressed genes could be assigned as tightly linked to the ASGR by screening of apomictic and sexual F₁s. Only one transcript, which did not map to the ASGR, showed expression primarily in reproductive tissue.

Conclusions

Our results suggest that a strategy of comparative sequencing of transcriptomes between donor parent and backcross lines containing an alien chromosome of interest can be an efficient method of identifying transcripts derived from an alien chromosome in a chromosome addition line.

Differences in gene expression amplitude overlie a conserved transcriptomic program occurring between the rapid and potent localized resistant reaction at the syncytium of the Glycine max genotype Peking (PI 548402) as compared to the prolonged and potent resistant reaction of PI 88788

Glycine max L. Merr. (soybean) resistance to Heterodera glycines Ichinohe occurs at the site of infection, a nurse cell known as the syncytium. Resistance is classified into two cytologically-defined responses, the G. max ([Peking])- and G. max ([PI 88788])-types. Each type represents a cohort of G. max genotypes. Resistance in G. max ([Peking]) occurs by a potent and rapid localized response, affecting parasitic second stage juveniles (p-J2). In contrast, resistance occurs by a potent but more prolonged reaction in the genotype G. max ([PI 88788]) that affects nematode development at the J3 and J4 stages. Microarray analyses comparing these cytologically and developmentally distinct resistant reactions reveal differences in gene expression in pericycle and surrounding cells even before infection. The differences include higher relative levels of the differentially expressed in response to arachidonic acid 1 gene (DEA1 [Gm-DEA1]) (+224.19-fold) and a protease inhibitor (+68.28-fold) in G. max ([Peking/PI 548402]) as compared to G. max ([PI 88788]). Gene pathway analyses compare the two genotypes (1) before, (2) at various times during, (3) constitutively throughout the resistant reaction and (4) at all time points prior to and during the resistant reaction. The amplified levels of transcriptional activity of defense genes may explain the rapid and potent reaction in G. max ([Peking/PI 548402]) as compared to G. max ([PI 88788]). In contrast, the shared differential expression levels of genes in G. max ([Peking/PI 548402]) and G. max ([PI 88788]) may indicate a conserved genomic program underlying the G. max resistance on which the genotype-specific gene expression programs are built off.

Genomic diversity and evolution of Mycobacterium ulcerans revealed by next-generation sequencing

Mycobacterium ulcerans is the causative agent of Buruli ulcer, the third most common mycobacterial disease after tuberculosis and leprosy. It is an emerging infectious disease that afflicts mainly children and youths in West Africa. Little is known about the evolution and transmission mode of M. ulcerans, partially due to the lack of known genetic polymorphisms among isolates, limiting the application of genetic epidemiology. To systematically profile single nucleotide polymorphisms (SNPs), we sequenced the genomes of three M. ulcerans strains using 454 and Solexa technologies. Comparison with the reference genome of the Ghanaian classical lineage isolate Agy99 revealed 26,564 SNPs in a Japanese strain representing the ancestral lineage. Only 173 SNPs were found when comparing Agy99 with two other Ghanaian isolates, which belong to the two other types previously distinguished in Ghana by variable number tandem repeat typing. We further analyzed a collection of Ghanaian strains using the SNPs discovered. With 68 SNP loci, we were able to differentiate 54 strains into 13 distinct SNP haplotypes. The average SNP nucleotide diversity was low (average 0.06–0.09 across 68 SNP loci), and 96% of the SNP locus pairs were in complete linkage disequilibrium. We estimated that the divergence of the M. ulcerans Ghanaian clade from the Japanese strain occurred 394 to 529 thousand years ago. The Ghanaian subtypes diverged about 1000 to 3000 years ago, or even much more recently, because we found evidence that they evolved significantly faster than average. Our results offer significant insight into the evolution of M. ulcerans and provide a comprehensive report on genetic diversity within a highly clonal M. ulcerans population from a Buruli ulcer endemic region, which can facilitate further epidemiological studies of this pathogen through the development of high-resolution tools.

Mycobacterium ulcerans is the causative agent of Buruli ulcer (BU), a necrotizing skin disease and the third most common mycobacterial disease after tuberculosis and leprosy. It is an emerging infectious disease that afflicts mainly children and youths in West Africa. The disease is also found in tropical and subtropical regions of Asia, the Western Pacific, and Latin America. Limited knowledge of this neglected tropical disease is partially due to the lack of known genetic polymorphisms among isolates, which hinder the study of transmission, epidemiology, and evolution of M. ulcerans. Our aim is to systematically profile genetic diversity among M. ulcerans isolates by sequencing and comparing the genomes of selected strains. We identified single nucleotide polymorphisms (SNPs) within a highly clonal M. ulcerans population from a Buruli ulcer endemic region. Based on the SNPs discovered, we developed SNP typing assays and were able to differentiate a collection of M. ulcerans isolates from this Buruli ulcer endemic region into 13 SNP haplotypes. Our results lay the ground for developing a highly discriminatory and cost-effective tool to study M. ulcerans evolution and epidemiology at a population level.

From parasite genomes to one healthy world: Are we having fun yet?

In 1990, the Human Genome Sequencing Project was established. This laid the ground work for an explosion of sequence data that has since followed. As a result of this effort, the first complete genome of an animal, Caenorhabditis elegans was published in 1998. The sequence of Drosophila melanogaster was made available in March, 2000 and in the following year, working drafts of the human genome were generated with the completed sequence (92%) being released in 2003. Recent advancements and next-generation technologies have made sequencing common place and have infiltrated every aspect of biological research, including parasitology. To date, sequencing of 32 apicomplexa and 24 nematode genomes are either in progress or near completion, and over 600k nematode EST and 200k apicomplexa EST submissions fill the databases. However, the winds have shifted and efforts are now refocusing on how best to store, mine and apply these data to problem solving. Herein we tend not to summarize existing X-omics datasets or present new technological advances that promise future benefits. Rather, the information to follow condenses up-to-date-applications of existing technologies to problem solving as it relates to parasite research. Advancements in non-parasite systems are also presented with the proviso that applications to parasite research are in the making.

Population-specific gene expression in the plant pathogenic nematode Heterodera glycines exists prior to infection and during the onset of a resistant or susceptible reaction in the roots of the Glycine max genotype Peking

BACKGROUND

A single Glycine max (soybean) genotype (Peking) reacts differently to two different populations of Heterodera glycines (soybean cyst nematode) within the first twelve hours of infection during resistant (R) and susceptible (S) reactions. This suggested that H. glycines has population-specific gene expression signatures. A microarray analysis of 7539 probe sets representing 7431 transcripts on the Affymetrix soybean GeneChip were used to identify population-specific gene expression signatures in pre-infective second stage larva (pi-L2) prior to their infection of Peking. Other analyses focused on the infective L2 at 12 hours post infection (i-L2(12h)), and the infective sedentary stages at 3 days post infection (i-L2(3d)) and 8 days post infection (i-L2/L3(8d)).

RESULTS

Differential expression and false discovery rate (FDR) analyses comparing populations of pi-L2 (i.e., incompatible population, NL1-RHg to compatible population, TN8) identified 71 genes that were induced in NL1-RHg as compared to TN8. These genes included putative gland protein G23G12, putative esophageal gland protein Hgg-20 and arginine kinase. The comparative analysis of pi-L2 identified 44 genes that were suppressed in NL1-RHg as compared to TN8. These genes included a different Hgg-20 gene, an EXPB1 protein and a cuticular collagen. By 12 h, there were 7 induced genes and 0 suppressed genes in NL1-RHg. By 3d, there were 9 induced and 10 suppressed genes in NL1-RHg. Substantial changes in gene expression became evident subsequently. At 8d there were 13 induced genes in NL1-RHg. This included putative gland protein G20E03, ubiquitin extension protein, putative gland protein G30C02 and beta-1,4 endoglucanase. However, 1668 genes were found to be suppressed in NL1-RHg. These genes included steroid alpha reductase, serine proteinase and a collagen protein.

CONCLUSION

These analyses identify a genetic expression signature for these two populations both prior to and subsequently as they undergo an R or S reaction. The identification of genes like steroid alpha reductase and serine proteinase that are involved in feeding and nutritional uptake as being highly suppressed during the R response at 8d may indicate genes that the plant is targeting. The analyses also identified numerous putative parasitism genes that are differentially expressed. The 1668 genes that are suppressed in NL1-RHg, and hence induced in TN8 may represent genes that are important during the parasitic stages of H. glycines development. The potential for different arrays of putative parasitism genes to be expressed in different nematode populations may indicate how H. glycines evolve mechanisms to overcome resistance.

Use of high-density tiling microarrays to identify mutations globally and elucidate mechanisms of drug resistance in Plasmodium falciparum

Using tiling microarrays, a mechanism by which Plasmodium falciparum parasites acquire resistance to the antimalarial fosmidomycin has been elucidated

Background

The identification of genetic changes that confer drug resistance or other phenotypic changes in pathogens can help optimize treatment strategies, support the development of new therapeutic agents, and provide information about the likely function of genes. Elucidating mechanisms of phenotypic drug resistance can also assist in identifying the mode of action of uncharacterized but potent antimalarial compounds identified in high-throughput chemical screening campaigns against Plasmodium falciparum.

Results

Here we show that tiling microarrays can detect de novo a large proportion of the genetic changes that differentiate one genome from another. We show that we detect most single nucleotide polymorphisms or small insertion deletion events and all known copy number variations that distinguish three laboratory isolates using readily accessible methods. We used the approach to discover mutations that occur during the selection process after transfection. We also elucidated a mechanism by which parasites acquire resistance to the antimalarial fosmidomycin, which targets the parasite isoprenoid synthesis pathway. Our microarray-based approach allowed us to attribute in vitro derived fosmidomycin resistance to a copy number variation event in the pfdxr gene, which enables the parasite to overcome fosmidomycin-mediated inhibition of isoprenoid biosynthesis.

Conclusions

We show that newly emerged single nucleotide polymorphisms can readily be detected and that malaria parasites can rapidly acquire gene amplifications in response to in vitro drug pressure. The ability to define comprehensively genetic variability in P. falciparum with a single overnight hybridization creates new opportunities to study parasite evolution and improve the treatment and control of malaria.

RNAi and functional genomics in plant parasitic nematodes

Plant nematology is currently undergoing a revolution with the availability of the first genome sequences as well as comprehensive expressed sequence tag (EST) libraries from a range of nematode species. Several strategies are being used to exploit this wealth of information. Comparative genomics is being used to explore the acquisition of novel genes associated with parasitic lifestyles. Functional analyses of nematode genes are moving toward larger scale studies including global transcriptome profiling. RNA interference (RNAi) has been shown to reduce expression of a range of plant parasitic nematode genes and is a powerful tool for functional analysis of nematode genes. RNAi-mediated suppression of genes essential for nematode development, survival, or parasitism is revealing new targets for nematode control. Plant nematology in the genomics era is now facing the challenge to develop RNAi screens adequate for high-throughput functional analyses.

High-resolution melt analysis to identify and map sequence-tagged site anchor points onto linkage maps: a white lupin (Lupinus albus) map as an exemplar

* The provision of sequence-tagged site (STS) anchor points allows meaningful comparisons between mapping studies but can be a time-consuming process for nonmodel species or orphan crops. * Here, the first use of high-resolution melt analysis (HRM) to generate STS markers for use in linkage mapping is described. This strategy is rapid and low-cost, and circumvents the need for labelled primers or amplicon fractionation. * Using white lupin (Lupinus albus, x = 25) as a case study, HRM analysis was applied to identify 91 polymorphic markers from expressed sequence tag (EST)-derived and genomic libraries. Of these, 77 generated STS anchor points in the first fully resolved linkage map of the species. The map also included 230 amplified fragment length polymorphisms (AFLP) loci, spanned 1916 cM (84.2% coverage) and divided into the expected 25 linkage groups. * Quantitative trait loci (QTL) analyses performed on the population revealed genomic regions associated with several traits, including the agronomically important time to flowering (tf), alkaloid synthesis and stem height (Ph). Use of HRM-STS markers also allowed us to make direct comparisons between our map and that of the related crop, Lupinus angustifolius, based on the conversion of RFLP, microsatellite and single nucleotide polymorphism (SNP) markers into HRM markers.