Simon R. Ellwood

Construction of a comparative genetic map in faba bean (Vicia faba L.); conservation of genome structure with Lens culinaris

BackgroundThe development of genetic markers is complex and costly in species with little pre-existing genomic information. Faba bean possesses one of the largest and least studied genomes among cultivated crop plants and no gene-based genetic maps exist. Gene-based orthologous markers allow chromosomal regions and levels of synteny to be characterised between species, reveal phylogenetic relationships and chromosomal evolution, and enable targeted identification of markers for crop breeding. In this study orthologous codominant cross-species markers have been deployed to produce the first exclusively gene-based genetic linkage map of faba bean (Vicia faba), using an F6 population developed from a cross between the lines Vf6 (equina type) and Vf27 (paucijuga type).ResultsOf 796 intron-targeted amplified polymorphic (ITAP) markers screened, 151 markers could be used to construct a comparative genetic map. Linkage analysis revealed seven major and five small linkage groups (LGs), one pair and 12 unlinked markers. Each LG was comprised of three to 30 markers and varied in length from 23.6 cM to 324.8 cM. The map spanned a total length of 1685.8 cM. A simple and direct macrosyntenic relationship between faba bean and Medicago truncatula was evident, while faba bean and lentil shared a common rearrangement relative to M. truncatula. One hundred and four of the 127 mapped markers in the 12 LGs, which were previously assigned to M. truncatula genetic and physical maps, were found in regions syntenic between the faba bean and M. truncatula genomes. However chromosomal rearrangements were observed that could explain the difference in chromosome numbers between these three legume species. These rearrangements suggested high conservation of M. truncatula chromosomes 1, 5 and 8; moderate conservation of chromosomes 2, 3, 4 and 7 and no conservation with M. truncatula chromosome 6. Multiple PCR amplicons and comparative mapping were suggestive of small-scale duplication events in faba bean. This study also provides a preliminary indication for finer scale macrosynteny between M. truncatula, lentil and faba bean. Markers originally designed from genes on the same M. truncatula BACs were found to be grouped together in corresponding syntenic areas in lentil and faba bean.ConclusionDespite the large size of the faba bean genome, comparative mapping did not reveal evidence for polyploidisation, segmental duplication, or significant rearrangements compared to M. truncatula, although a bias in the use of single locus markers may have limited the detection of duplications. Non-coding repetitive DNA or transposable element content provides a possible explanation for the difference in genome sizes. Similar patterns of rearrangements in faba bean and lentil compared to M. truncatula support phylogenetic studies dividing these species into the tribes Viceae and Trifoliae. However, substantial macrosynteny was apparent between faba bean and M. truncatula, with the exception of chromosome 6 where no orthologous markers were found, confirming previous investigations suggesting chromosome 6 is atypical. The composite map, anchored with orthologous markers mapped in M. truncatula, provides a central reference map for future use of genomic and genetic information in faba bean genetic analysis and breeding.

Comparison of Erysiphe cichoracearum and E. cruciferarum and a Survey of 360 Arabidopsis thaliana Accessions for Resistance to These Two Powdery Mildew Pathogens

In previous work, UEA1 and UCSC1, two geographically distinct, powdery mildew isolates, were recognized for their ability to infect Arabidopsis thaliana. We have clarified the identity of these isolates by determining their host ranges, reexamining their morphology, and comparing their DNA sequences for the 5.8S ribosomal RNA and two flanking internal transcribed spacer sequences. These experiments confirm that UEA1 is a member of Erysiphe cruciferarum and that UCSC1 belongs to E. cichoracearum. Interactions of the two Erysiphe isolates with 360 A. thaliana accessions were examined to provide a comprehensive profile of naturally occurring powdery mildew resistance in this weedy species. The majority of A. thaliana accessions (213) were susceptible to both isolates. Among the accessions exhibiting some degree of resistance, most (84) responded differentially to UEA1 and UCSC1 and the remainder were resistant to both isolates. Notably, resistance to UCSC1 cosegregated with RPW7, a locus previously demonstrated to confer resistance to UEA1 in Ms-0 x Landsberg (erecta) crosses. With this large collection of resistant accessions, questions about species specificity, genetic diversity and the evolution of resistance to powdery mildews can be addressed.

Pyrenophora teres: profile of an increasingly damaging barley pathogen

UNLABELLED Pyrenophora teres, causal agent of net blotch of barley, exists in two forms, designated P. teres f. teres and P. teres f. maculata, which induce net form net blotch (NFNB) and spot form net blotch (SFNB), respectively. Significantly more work has been performed on the net form than on the spot form although recent activity in spot form research has increased because of epidemics of SFNB in barley-producing regions. Genetic studies have demonstrated that NFNB resistance in barley is present in both dominant and recessive forms, and that resistance/susceptibility to both forms can be conferred by major genes, although minor quantitative trait loci have also been identified. Early work on the virulence of the pathogen showed toxin effector production to be important in disease induction by both forms of pathogen. Since then, several laboratories have investigated effectors of virulence and avirulence, and both forms are complex in their interaction with the host. Here, we assemble recent information from the literature that describes both forms of this important pathogen and includes reports describing the host-pathogen interaction with barley. We also include preliminary findings from a genome sequence survey. TAXONOMY Pyrenophora teres Drechs. Kingdom Fungi; Phylum Ascomycota; Subphylum Pezizomycotina; Class Dothideomycete; Order Pleosporales; Family Pleosporaceae; Genus Pyrenophora, form teres and form maculata. IDENTIFICATION To date, no clear morphological or life cycle differences between the two forms of P. teres have been identified, and therefore they are described collectively. Towards the end of the growing season, the fungus produces dark, globosely shaped pseudothecia, about 1-2mm in diameter, on barley. Ascospores measuring 18-28µm × 43-61µm are light brown and ellipsoidal and often have three to four transverse septa and one or two longitudinal septa in the median cells. Conidiophores usually arise singly or in groups of two or three and are lightly swollen at the base. Conidia measuring 30-174µm × 15-23µm are smoothly cylindrical and straight, round at both ends, subhyaline to yellowish brown, often with four to six pseudosepta. Morphologically, P. teres f. teres and P. teres f. maculata are indistinguishable. HOST RANGE Comprehensive work on the host range of P. teres f. teres has been performed; however, little information on the host range of P. teres f. maculata is available. Hordeum vulgare and H. vulgare ssp. spontaneum are considered to be the primary hosts for P. teres. However, natural infection by P. teres has been observed in other wild Hordeum species and related species from the genera Bromus, Avena and Triticum, including H. marinum, H. murinum, H. brachyantherum, H. distichon, H. hystrix, B. diandrus, A. fatua, A. sativa and T. aestivum (Shipton et al., 1973, Rev. Plant Pathol. 52:269-290). In artificial inoculation experiments under field conditions, P. teres f. teres has been shown to infect a wide range of gramineous species in the genera Agropyron, Brachypodium, Elymus, Cynodon, Deschampsia, Hordelymus and Stipa (Brown et al., 1993, Plant Dis. 77:942-947). Additionally, 43 gramineous species were used in a growth chamber study and at least one of the P. teres f. teres isolates used was able to infect 28 of the 43 species tested. However, of these 28 species, 14 exhibited weak type 1 or 2 reactions on the NFNB 1-10 scale (Tekauz, 1985). These reaction types are small pin-point lesions and could possibly be interpreted as nonhost reactions. In addition, the P. teres f. teres host range was investigated under field conditions by artificially inoculating 95 gramineous species with naturally infected barley straw. Pyrenophora teres f. teres was re-isolated from 65 of the species when infected leaves of adult plants were incubated on nutrient agar plates; however, other than Hordeum species, only two of the 65 host species exhibited moderately susceptible or susceptible field reaction types, with most species showing small dark necrotic lesions indicative of a highly resistant response to P. teres f. teres. Although these wild species have the potential to be alternative hosts, the high level of resistance identified for most of the species makes their role as a source of primary inoculum questionable. DISEASE SYMPTOMS Two types of symptom are caused by P. teres. These are net-type lesions caused by P. teres f. teres and spot-type lesions caused by P. teres f. maculata. The net-like symptom, for which the disease was originally named, has characteristic narrow, dark-brown, longitudinal and transverse striations on infected leaves. The spot form symptom consists of dark-brown, circular to elliptical lesions surrounded by a chlorotic or necrotic halo of varying width.

SSR analysis of the Medicago truncatula SARDI core collection reveals substantial diversity and unusual genotype dispersal throughout the Mediterranean basin

The world’s oldest and largest Medicago truncatula collection is housed at the South Australian Research and Development Institute (SARDI). We used six simple sequence repeat (SSR) loci to analyse the genetic diversity and relationships between randomly selected individuals from 192 accessions in the core collection. M. truncatula is composed of three subspecies (ssp.): ssp. truncatula, ssp. longeaculeata, and ssp. tricycla. Analysis at the level of six SSR loci supports the concept of ssp. tricycla, all the samples of which showed unique alleles at two loci. Contingency Chi-squared tests were significant between ssp. tricycla and ssp. truncatula at four loci, suggesting a barrier to gene flow between these subspecies. In accessions defined as ssp. longeaculeata, no unique allelic distribution or diagnostic sizes were observed, suggesting this apparent ssp. is a morphological variant of ssp. truncatula. The data also suggest M. truncatula that exhibits unusually wide genotype dispersal throughout its native Mediterranean region, possibly due to animal and trade-related movements. Our results showed the collection to be highly diverse, exhibiting an average of 25 SSR alleles per locus, with over 90% of individuals showing discrete genotypes. The rich diversity of the SARDI collection provides an invaluable resource for studying natural allelic variation of M. truncatula. To efficiently exploit the variation in the SARDI collection, we have defined a subset of accessions (n=61) that maximises the diversity.

The first genetic and comparative map of white lupin (Lupinus albus l.): Identification of QTLs for anthracnose resistance and flowering time, and a locus for alkaloid content

Abstract We report the first genetic linkage map of white lupin (Lupinus albus L.). An F8 recombinant inbred line population developed from Kiev mutant × P27174 was mapped with 220 amplified fragment length polymorphism and 105 gene-based markers. The genetic map consists of 28 main linkage groups (LGs) that varied in length from 22.7 cM to 246.5 cM and spanned a total length of 2951 cM. There were seven additional pairs and 15 unlinked markers, and 12.8% of markers showed segregation distortion at P < 0.05. Syntenic relationships between Medicago truncatula and L. albus were complex. Forty-five orthologous markers that mapped between M. truncatula and L. albus identified 17 small syntenic blocks, and each M. truncatula chromosome aligned to between one and six syntenic blocks in L. albus. Genetic mapping of three important traits: anthracnose resistance, flowering time, and alkaloid content allowed loci governing these traits to be defined. Two quantitative trait loci (QTLs) with significant effects were identified for anthracnose resistance on LG4 and LG17, and two QTLs were detected for flowering time on the top of LG1 and LG3. Alkaloid content was mapped as a Mendelian trait to LG11.

A first genome assembly of the barley fungal pathogen Pyrenophora teres f. teres

BackgroundPyrenophora teres f. teres is a necrotrophic fungal pathogen and the cause of one of barleys most important diseases, net form of net blotch. Here we report the first genome assembly for this species based solely on short Solexa sequencing reads of isolate 0-1. The assembly was validated by comparison to BAC sequences, ESTs, orthologous genes and by PCR, and complemented by cytogenetic karyotyping and the first genome-wide genetic map for P. teres f. teres.ResultsThe total assembly was 41.95 Mbp and contains 11,799 gene models of 50 amino acids or more. Comparison against two sequenced BACs showed that complex regions with a high GC content assembled effectively. Electrophoretic karyotyping showed distinct chromosomal polymorphisms between isolates 0-1 and 15A, and cytological karyotyping confirmed the presence of at least nine chromosomes. The genetic map spans 2477.7 cM and is composed of 243 markers in 25 linkage groups, and incorporates simple sequence repeat markers developed from the assembly. Among predicted genes, non-ribosomal peptide synthetases and efflux pumps in particular appear to have undergone a P. teres f. teres-specific expansion of non-orthologous gene families.ConclusionsThis study demonstrates that paired-end Solexa sequencing can successfully capture coding regions of a filamentous fungal genome. The assembly contains a plethora of predicted genes that have been implicated in a necrotrophic lifestyle and pathogenicity and presents a significant resource for examining the bases for P. teres f. teres pathogenicity.

Aligning a New Reference Genetic Map of Lupinus angustifolius with the Genome Sequence of the Model Legume, Lotus japonicus

We have developed a dense reference genetic map of Lupinus angustifolius (2n = 40) based on a set of 106 publicly available recombinant inbred lines derived from a cross between domesticated and wild parental lines. The map comprised 1090 loci in 20 linkage groups and three small clusters, drawing together data from several previous mapping publications plus almost 200 new markers, of which 63 were gene-based markers. A total of 171 mainly gene-based, sequence-tagged site loci served as bridging points for comparing the Lu. angustifolius genome with the genome sequence of the model legume, Lotus japonicus via BLASTn homology searching. Comparative analysis indicated that the genomes of Lu. angustifolius and Lo. japonicus are highly diverged structurally but with significant regions of conserved synteny including the region of the Lu. angustifolius genome containing the pod-shatter resistance gene, lentus. We discuss the potential of synteny analysis for identifying candidate genes for domestication traits in Lu. angustifolius and in improving our understanding of Fabaceae genome evolution.

CodingQuarry: highly accurate hidden Markov model gene prediction in fungal genomes using RNA-seq transcripts

BackgroundThe impact of gene annotation quality on functional and comparative genomics makes gene prediction an important process, particularly in non-model species, including many fungi. Sets of homologous protein sequences are rarely complete with respect to the fungal species of interest and are often small or unreliable, especially when closely related species have not been sequenced or annotated in detail. In these cases, protein homology-based evidence fails to correctly annotate many genes, or significantly improve ab initio predictions. Generalised hidden Markov models (GHMM) have proven to be invaluable tools in gene annotation and, recently, RNA-seq has emerged as a cost-effective means to significantly improve the quality of automated gene annotation. As these methods do not require sets of homologous proteins, improving gene prediction from these resources is of benefit to fungal researchers. While many pipelines now incorporate RNA-seq data in training GHMMs, there has been relatively little investigation into additionally combining RNA-seq data at the point of prediction, and room for improvement in this area motivates this study.ResultsCodingQuarry is a highly accurate, self-training GHMM fungal gene predictor designed to work with assembled, aligned RNA-seq transcripts. RNA-seq data informs annotations both during gene-model training and in prediction. Our approach capitalises on the high quality of fungal transcript assemblies by incorporating predictions made directly from transcript sequences. Correct predictions are made despite transcript assembly problems, including those caused by overlap between the transcripts of adjacent gene loci.Stringent benchmarking against high-confidence annotation subsets showed CodingQuarry predicted 91.3% of Schizosaccharomyces pombe genes and 90.4% of Saccharomyces cerevisiae genes perfectly. These results are 4-5% better than those of AUGUSTUS, the next best performing RNA-seq driven gene predictor tested. Comparisons against whole genome Sc. pombe and S. cerevisiae annotations further substantiate a 4-5% improvement in the number of correctly predicted genes.ConclusionsWe demonstrate the success of a novel method of incorporating RNA-seq data into GHMM fungal gene prediction. This shows that a high quality annotation can be achieved without relying on protein homology or a training set of genes. CodingQuarry is freely available (https://sourceforge.net/projects/codingquarry/), and suitable for incorporation into genome annotation pipelines.

Differences in syntenic complexity between Medicago truncatula with Lens culinaris and Lupinus albus

Orthologous markers transferable between distantly related legume species allow for the rapid generation of genetic maps in species where there is little pre-existing genomic or EST information. We are using the model legume Medicago truncatula Gaertn. to develop such markers in legumes of importance to Australian agriculture. This will enable the construction of comparative genetic maps, help to determine patterns of chromosomal evolution in the legume family, and characterise syntenic relationships between M. truncatula and cultivated legumes. This information can then be used to identify markers that are tightly linked to the genes of interest, candidate gene(s) for a trait, and expedite the isolation of such genes. Among the Papilionoideae, we compared ESTs from the phylogenetically distant species, M. truncatula, Lupinus albus and Glycine max, to produce 500 intron-targeted amplified polymorphic markers (ITAPs). In addition to 126 M. truncatula cross-species markers from Department of Plant Pathology, University of California (USA), these markers were used to generate comparative genetic maps of lentil (Lens culinaris Medik.) and white lupin (Lupinus albus Linn.). Our results showed that 90% of the ITAPs markers amplified genomic DNA in M. truncatula, 80% in Lupinus albus, and 70% in Lens culinaris. The comparative map of Lens culinaris was constructed based on 79 ITAP markers. The Lupinus albus comparative map was developed from 105 gene-based markers together with 223 AFLP markers. Although a direct and simple syntenic relationship was observed between M. truncatula and Lens culinaris genomes, there is evidence of moderate chromosomal rearrangement. This may account for the different chromosome numbers in the two species. A more complicated pattern among homologous blocks was apparent between the Lupinus albus and M. truncatula genomes.

Two alternative recessive quantitative trait loci influence resistance to spring black stem and leaf spot in Medicago truncatula

BackgroundKnowledge of the genetic basis of plant resistance to necrotrophic pathogens is incomplete and has been characterised in relatively few pathosystems. In this study, the cytology and genetics of resistance to spring black stem and leaf spot caused by Phoma medicaginis, an economically important necrotrophic pathogen of Medicago spp., was examined in the model legume M. truncatula.ResultsMacroscopically, the resistant response of accession SA27063 was characterised by small, hypersensitive-like spots following inoculation while the susceptible interaction with accessions A17 and SA3054 showed necrotic lesions and spreading chlorosis. No unique cytological differences were observed during early infection (<48 h) between the resistant and susceptible genotypes, except pathogen growth was restricted to one or a few host cells in SA27063. In both interactions reactive oxygen intermediates and phenolic compounds were produced, and cell death occurred. Two F2 populations segregating for resistance to spring black stem and leaf spot were established between SA27063 and the two susceptible accessions, A17 and SA3054. The cross between SA27063 and A17 represented a wider cross than between SA27063 and SA3054, as evidenced by higher genetic polymorphism, reduced fertility and aberrant phenotypes of F2 progeny. In the SA27063 × A17 F2 population a highly significant quantitative trait locus (QTL, LOD = 7.37; P < 0.00001) named resistance to the necrotroph Phomamedicaginis one (rnpm1) genetically mapped to the top arm of linkage group 4 (LG4). rnpm1 explained 33.6% of the phenotypic variance in the populations response to infection depicted on a 1–5 scale and was tightly linked to marker AW256637. A second highly significant QTL (LOD = 6.77; P < 0.00001), rnpm2, was located on the lower arm of LG8 in the SA27063 × SA3054 map. rnpm2 explained 29.6% of the phenotypic variance and was fine mapped to a 0.8 cM interval between markers h2_16a6a and h2_21h11d. rnpm1 is tightly linked to a cluster of Toll/Interleukin1 receptor-nucleotide binding site-leucine-rich repeat (TIR-NBS-LRR) genes and disease resistance protein-like genes, while no resistance gene analogues (RGAs) are apparent in the genomic sequence of the reference accession A17 at the rnpm2 locus.ConclusionThe induction of defence responses and cell death in the susceptible interaction following infection by P. medicaginis suggested this pathogen is not negatively affected by these responses and may promote them. A QTL for resistance was revealed in each of two populations derived from crosses between a resistant accession and two different susceptible accessions. Both loci are recessive in nature, and the simplest explanation for the existence of two separate QTLs is the occurrence of host genotype-specific susceptibility loci that may interact with undetermined P. medicaginis virulence factors.