Estimates of conserved microsynteny among the genomes of Glycine max, Medicago truncatula and Arabidopsis thaliana

A methodological scoping review of the integration of fMRI to guide dMRI tractography. What has been done and what can be improved: A 20-year perspective

Combining MRI modalities is a growing trend in neurosciences. It provides opportunities to investigate the brain architecture supporting cognitive functions. Integrating fMRI activation to guide dMRI tractography offers potential advantages over standard tractography methods. A quick glimpse of the literature on this topic reveals that this technique is challenging, and no consensus or "best practices" currently exist, at least not within a single document. We present the first attempt to systematically analyze and summarize the literature of 80 studies that integrated task-based fMRI results to guide tractography, over the last two decades. We report 19 findings that cover challenges related to sample size, microstructure modelling, seeding methods, multimodal space registration, false negatives/positives, specificity/validity, gray/white matter interface and more. These findings will help the scientific community (1) understand the strengths and limitations of the approaches, (2) design studies using this integrative framework, and (3) motivate researchers to fill the gaps identified. We provide references toward best practices, in order to improve the overall result's replicability, sensitivity, specificity, and validity.

Structural complexity and functional diversity of plant NADPH oxidases

Plant NADPH oxidases also known as respiratory burst oxidase homologs (Rbohs) are a family of membrane-bound enzymes that play diverse roles in the defense response and morphogenetic processes via regulated generation of reactive oxygen species. Rbohs are associated with a variety of functions, although the reason for this is not clear. To evaluate using bioinformatics, the possible mechanisms for the observed functional diversity within the plant kingdom, 127 Rboh protein sequences representing 26 plant species were analyzed. Multiple clusters were identified with gene duplications that were both dicot as well as monocot-specific. The N-terminal sequences were observed to be highly variable. The conserved cysteine (equivalent of Cys890) in C-terminal of AtRbohD suggested that the redox-based modification like S-nitrosylation may regulate the activity of other Rbohs. Three-dimensional models corresponding to the N-terminal domain for Rbohs from Arabidopsis thaliana and Oryza sativa were constructed and molecular dynamics studies were carried out to study the role of Ca2+ in the folding of Rboh proteins. Certain mutations indicated possibly affect the structure and function of the plant NADPH oxidases, thereby providing the rationale for further experimental validation.

Fighting Asian Soybean Rust

Phakopsora pachyrhizi is a biotrophic fungus provoking SBR disease. SBR poses a major threat to global soybean production. Though several R genes provided soybean immunity to certain P. pachyrhizi races, the pathogen swiftly overcame this resistance. Therefore, fungicides are the only current means to control SBR. However, insensitivity to fungicides is soaring in P. pachyrhizi and, therefore, alternative measures are needed for SBR control. In this article, we discuss the different approaches for fighting SBR and their potential, disadvantages, and advantages over other measures. These encompass conventional breeding for SBR resistance, transgenic approaches, exploitation of transcription factors, secondary metabolites, and antimicrobial peptides, RNAi/HIGS, and biocontrol strategies. It seems that an integrating approach exploiting different measures is likely to provide the best possible means for the effective control of SBR.

Medicago truncatula and Glycine max: Different Drought Tolerance and Similar Local Response of the Root Nodule Proteome

Legume crops present important agronomical and environmental advantages mainly due to their capacity to reduce atmospheric N2 to ammonium via symbiotic nitrogen fixation (SNF). This process is very sensitive to abiotic stresses such as drought, but the mechanism underlying this response is not fully understood. The goal of the current work is to compare the drought response of two legumes with high economic impact and research importance, Medicago truncatula and Glycine max, by characterizing their root nodule proteomes. Our results show that, although M. truncatula exhibits lower water potential values under drought conditions compared to G. max, SNF declined analogously in the two legumes. Both of their nodule proteomes are very similar, and comparable down-regulation responses in the diverse protein functional groups were identified (mainly proteins related to the metabolism of carbon, nitrogen, and sulfur). We suggest lipoxygenases and protein turnover as newly recognized players in SNF regulation. Partial drought conditions applied to a split-root system resulted in the local down-regulation of the entire proteome of drought-stressed nodules in both legumes. The high degree of similarity between both legume proteomes suggests that the vast amount of research conducted on M. truncatula could be applied to economically important legume crops, such as soybean.

Dynamic rearrangements determine genome organization and useful traits in soybean

Soybean (Glycine max) is a paleopolyploid whose genome has gone through at least two rounds of polyploidy and subsequent diploidization events. Several studies have investigated the changes in genome structure produced by the relatively recent polyploidy event, but little is known about the ancient polyploidy due to the high frequency of gene loss after duplication. Our previous study, regarding a region responsible for bacterial leaf pustule, reported two homeologous Rxp regions produced by the recent whole-genome duplication event. In this study, we identified the full set of four homeologous Rxp regions (ranging from 1.96 to 4.60 Mb) derived from both the recent and ancient polyploidy events, and this supports the quadruplicated structure of the soybean genome. Among the predicted genes on chromosome 17 (linkage group D2), 71% of them were conserved in a recently duplicated region, while 21% and 24% of duplicated genes were retained in two homeologous regions formed by the ancient polyploidy. Furthermore, comparative analysis showed a 2:1 relationship between soybean and Medicago truncatula, since M. truncatula did not undergo the recent polyploidy event that soybean did. Unlike soybean, M. truncatula homeologous regions were highly fractionated and their synteny did not exist, revealing different rates of diploidization process between the two species. Our data show that extensive synteny remained in the four homeologous regions in soybean, even though the soybean genome experienced dynamic genome rearrangements following paleopolyploidy events. Moreover, multiple Rxp quantitative trait loci on different soybean chromosomes actually comprise homeologous regions produced by two rounds of polyploidy events.

Using a model-based framework for analysing genetic diversity during germination and heterotrophic growth of Medicago truncatula

BACKGROUND AND AIMS

The framework provided by an emergence model was used: (1) for phenotyping germination and heterotrophic growth of Medicago truncatula in relation to two major environmental factors, temperature and water potential; and (2) to evaluate the extent of genetic differences in emergence-model parameters.

METHODS

Eight cultivars and natural accessions of M. trunculata were studied. Germination was recorded from 5 to 30 degrees C and from 0 to -0.75 MPa, and seedling growth from 10 to 20 degrees C.

KEY RESULTS

Thermal time to reach 50 % germination was very short (15 degrees Cd) and almost stable between genotypes, while base temperature (2-3 degrees C) and base water potential for germination (-0.7 to -1.3 MPa) varied between genotypes. Only 35 degrees Cd after germination were required to reach 30 mm hypocotyl length with significant differences among genotypes. Base temperature for elongation varied from 5.5 to 7.5 degrees C. Low temperatures induced a general shortening of the seedling, with some genotypes more responsive than others. No relationship with initial seed mass or seed reserve distribution was observed, which might have explained differences between genotypes and the effects of low temperatures.

CONCLUSIONS

The study provides a set of reference values for M. trunculata users. The use of the ecophysiological model allows comparison of these values between such non-crop species and other crops. It has enabled phenotypic variability in response to environmental conditions related to the emergence process to be identified. The model will allow simulation of emergence differences between genotypes in a range of environments using these parameter values. Genomic tools available for the model species M. trunculata will make it possible to analyse the genetic and molecular determinants of these differences.

The lipoxygenase gene family: a genomic fossil of shared polyploidy between Glycine max and Medicago truncatula

BACKGROUND

Soybean lipoxygenases (Lxs) play important roles in plant resistance and in conferring the distinct bean flavor. Lxs comprise a multi-gene family that includes GmLx1, GmLx2 and GmLx3, and many of these genes have been characterized. We were interested in investigating the relationship between the soybean lipoxygenase isozymes from an evolutionary perspective, since soybean has undergone two rounds of polyploidy. Here we report the tetrad genome structure of soybean Lx regions produced by ancient and recent polyploidy. Also, comparative genomics with Medicago truncatula was performed to estimate Lxs in the common ancestor of soybean and Medicago.

RESULTS

Two Lx regions in Medicago truncatula showing synteny with soybean were analyzed. Differential evolutionary rates between soybean and Medicago were observed and the median Ks values of Mt-Mt, Gm-Mt, and Gm-Gm paralogs were determined to be 0.75, 0.62, and 0.46, respectively. Thus the comparison of Gm-Mt paralogs (Ks = 0.62) and Gm-Mt orthologs (Ks = 0.45) supports the ancient duplication of Lx regions in the common ancestor prior to the Medicago-Glycine split. After speciation, no Lx regions generated by another polyploidy were identified in Medicago. Instead tandem duplication of Lx genes was observed. On the other hand, a lineage-specific duplication occurred in soybean resulting in two pairs of Lx regions. Each pair of soybean regions was co-orthologous to one Lx region in Medicago. A total of 34 Lx genes (15 MtLxs and 19 GmLxs) were divided into two groups by phylogenetic analysis. Our study shows that the Lx gene family evolved from two distinct Lx genes in the most recent common ancestor.

CONCLUSION

This study analyzed two pairs of Lx regions generated by two rounds of polyploidy in soybean. Each pair of soybean homeologous regions is co-orthologous to one region of Medicago, demonstrating the quartet structure of the soybean genome. Differential evolutionary rates between soybean and Medicago were observed; thus optimized rates of Ks per year should be applied for accurate estimation of coalescence times to each case of comparison: soybean-soybean, soybean-Medicago, or Medicago-Medicago. In conclusion, the soybean Lx gene family expanded by ancient polyploidy prior to taxon divergence, followed by a soybean- specific duplication and tandem duplications, respectively.

Sequence level analysis of recently duplicated regions in soybean [Glycine max (L.) Merr.] genome

A single recessive gene, rxp, on linkage group (LG) D2 controls bacterial leaf-pustule resistance in soybean. We identified two homoeologous contigs (GmA and GmA') composed of five bacterial artificial chromosomes (BACs) during the selection of BAC clones around Rxp region. With the recombinant inbred line population from the cross of Pureunkong and Jinpumkong 2, single-nucleotide polymorphism and simple sequence repeat marker genotyping were able to locate GmA' on LG A1. On the basis of information in the Soybean Breeders Toolbox and our results, parts of LG A1 and LG D2 share duplicated regions. Alignment and annotation revealed that many homoeologous regions contained kinases and proteins related to signal transduction pathway. Interestingly, inserted sequences from GmA and GmA' had homology with transposase and integrase. Estimation of evolutionary events revealed that speciation of soybean from Medicago and the recent divergence of two soybean homoeologous regions occurred at 60 and 12 million years ago, respectively. Distribution of synonymous substitution patterns, K(s), yielded a first secondary peak (mode K(s) = 0.10-0.15) followed by two smaller bulges were displayed between soybean homologous regions. Thus, diploidized paleopolyploidy of soybean genome was again supported by our study.

Sequencing and analysis of the gene-rich space of cowpea

Background

Cowpea, Vigna unguiculata (L.) Walp., is one of the most important food and forage legumes in the semi-arid tropics because of its drought tolerance and ability to grow on poor quality soils. Approximately 80% of cowpea production takes place in the dry savannahs of tropical West and Central Africa, mostly by poor subsistence farmers. Despite its economic and social importance in the developing world, cowpea remains to a large extent an underexploited crop. Among the major goals of cowpea breeding and improvement programs is the stacking of desirable agronomic traits, such as disease and pest resistance and response to abiotic stresses. Implementation of marker-assisted selection and breeding programs is severely limited by a paucity of trait-linked markers and a general lack of information on gene structure and organization. With a nuclear genome size estimated at ~620 Mb, the cowpea genome is an ideal target for reduced representation sequencing.

Results

We report here the sequencing and analysis of the gene-rich, hypomethylated portion of the cowpea genome selectively cloned by methylation filtration (MF) technology. Over 250,000 gene-space sequence reads (GSRs) with an average length of 610 bp were generated, yielding ~160 Mb of sequence information. The GSRs were assembled, annotated by BLAST homology searches of four public protein annotation databases and four plant proteomes (A. thaliana, M. truncatula, O. sativa, and P. trichocarpa), and analyzed using various domain and gene modeling tools. A total of 41,260 GSR assemblies and singletons were annotated, of which 19,786 have unique GenBank accession numbers. Within the GSR dataset, 29% of the sequences were annotated using the Arabidopsis Gene Ontology (GO) with the largest categories of assigned function being catalytic activity and metabolic processes, groups that include the majority of cellular enzymes and components of amino acid, carbohydrate and lipid metabolism. A total of 5,888 GSRs had homology to genes encoding transcription factors (TFs) and transcription associated factors (TAFs) representing about 5% of the total annotated sequences in the dataset. Sixty-two (62) of the 64 well-characterized plant transcription factor (TF) gene families are represented in the cowpea GSRs, and these families are of similar size and phylogenetic organization to those characterized in other plants. The cowpea GSRs also provides a rich source of genes involved in photoperiodic control, symbiosis, and defense-related responses. Comparisons to available databases revealed that about 74% of cowpea ESTs and 70% of all legume ESTs were represented in the GSR dataset. As approximately 12% of all GSRs contain an identifiable simple-sequence repeat, the dataset is a powerful resource for the design of microsatellite markers.

Conclusion

The availability of extensive publicly available genomic data for cowpea, a non-model legume with significant importance in the developing world, represents a significant step forward in legume research. Not only does the gene space sequence enable the detailed analysis of gene structure, gene family organization and phylogenetic relationships within cowpea, but it also facilitates the characterization of syntenic relationships with other cultivated and model legumes, and will contribute to determining patterns of chromosomal evolution in the Leguminosae. The micro and macrosyntenic relationships detected between cowpea and other cultivated and model legumes should simplify the identification of informative markers for marker-assisted trait selection and map-based gene isolation necessary for cowpea improvement.

Gene duplication and paleopolyploidy in soybean and the implications for whole genome sequencing

BACKGROUND

Soybean, Glycine max (L.) Merr., is a well documented paleopolyploid. What remains relatively under characterized is the level of sequence identity in retained homeologous regions of the genome. Recently, the Department of Energy Joint Genome Institute and United States Department of Agriculture jointly announced the sequencing of the soybean genome. One of the initial concerns is to what extent sequence identity in homeologous regions would have on whole genome shotgun sequence assembly.

RESULTS

Seventeen BACs representing approximately 2.03 Mb were sequenced as representative potential homeologous regions from the soybean genome. Genetic mapping of each BAC shows that 11 of the 20 chromosomes are represented. Sequence comparisons between homeologous BACs shows that the soybean genome is a mosaic of retained paleopolyploid regions. Some regions appear to be highly conserved while other regions have diverged significantly. Large-scale "batch" reassembly of all 17 BACs combined showed that even the most homeologous BACs with upwards of 95% sequence identity resolve into their respective homeologous sequences. Potential assembly errors were generated by tandemly duplicated pentatricopeptide repeat containing genes and long simple sequence repeats. Analysis of a whole-genome shotgun assembly of 80,000 randomly chosen JGI-DOE sequence traces reveals some new soybean-specific repeat sequences.

CONCLUSION

This analysis investigated both the structure of the paleopolyploid soybean genome and the potential effects retained homeology will have on assembling the whole genome shotgun sequence. Based upon these results, homeologous regions similar to those characterized here will not cause major assembly issues.

Gene duplication and paleopolyploidy in soybean and the implications for whole genome sequencing

BACKGROUND

Soybean, Glycine max (L.) Merr., is a well documented paleopolyploid. What remains relatively under characterized is the level of sequence identity in retained homeologous regions of the genome. Recently, the Department of Energy Joint Genome Institute and United States Department of Agriculture jointly announced the sequencing of the soybean genome. One of the initial concerns is to what extent sequence identity in homeologous regions would have on whole genome shotgun sequence assembly.

RESULTS

Seventeen BACs representing approximately 2.03 Mb were sequenced as representative potential homeologous regions from the soybean genome. Genetic mapping of each BAC shows that 11 of the 20 chromosomes are represented. Sequence comparisons between homeologous BACs shows that the soybean genome is a mosaic of retained paleopolyploid regions. Some regions appear to be highly conserved while other regions have diverged significantly. Large-scale "batch" reassembly of all 17 BACs combined showed that even the most homeologous BACs with upwards of 95% sequence identity resolve into their respective homeologous sequences. Potential assembly errors were generated by tandemly duplicated pentatricopeptide repeat containing genes and long simple sequence repeats. Analysis of a whole-genome shotgun assembly of 80,000 randomly chosen JGI-DOE sequence traces reveals some new soybean-specific repeat sequences.

CONCLUSION

This analysis investigated both the structure of the paleopolyploid soybean genome and the potential effects retained homeology will have on assembling the whole genome shotgun sequence. Based upon these results, homeologous regions similar to those characterized here will not cause major assembly issues.

Comparison between a coffee single copy chromosomal region and Arabidopsis duplicated counterparts evidenced high level synteny between the coffee genome and the ancestral Arabidopsis genome

The Arabidopsis thaliana genome sequence provides a catalogue of reference genes that can be used for comparative analysis of other species thereby facilitating map-based cloning in economically important crops. We made use of a coffee bacterial artificial chromosome (BAC) contig linked to the S(H)3 leaf rust resistance gene to assess microsynteny between coffee (Coffea arabica L.) and Arabidopsis. Microsynteny was revealed and the matching counterparts to C. arabica contigs were seen to be scattered throughout four different syntenic segments of Arabidopsis on chromosomes (Ath) I, III, IV and V. Coffee BAC filter hybridizations were performed using coffee putative conserved orthologous sequences to Arabidopsis predicted genes located on the different Arabidopsis syntenic regions. The coffee BAC contig related to the S(H)3 region was successfully consolidated and later on validated by fingerprinting. Furthermore, the anchoring markers appeared in same order on the coffee BAC contigs and in all Arabidopsis segments with the exception of a single inversion on AtIII and AtIV Arabidopsis segments. However, the S(H)3 coffee region appears to be closer to the ancestral genome segment (before the divergence of Arabidopsis and coffee) than any of the duplicated counterparts in the present-day Arabidopsis genome. The genome duplication events at the origin of its Arabidopsis counterparts occurred most probably after the separation (i.e. 94 million years ago) of Euasterid (Coffee) and Eurosid (Arabidopsis).

Structure of two melon regions reveals high microsynteny with sequenced plant species

In this study, two melon bacterial artificial chromosome (BAC) clones have been sequenced and annotated. BAC 1-21-10 spans 92 kb and contains the nsv locus conferring resistance to the Melon Necrotic Spot Virus (MNSV) in melon linkage group 11. BAC 13J4 spans 98 kb and belongs to a BAC contig containing resistance gene homologues, extending a previous sequenced region of 117 kb in linkage group 4. Both regions have microsyntenic relationships to the model plant species Arabidopsis thaliana, and to Medicago truncatula and Populus trichocarpa. The network of synteny found between melon and each of the sequenced genomes reflects the polyploid structure of Arabidopsis, Populus, and Medicago genomes due to whole genome duplications (WGD). A detailed analysis revealed that both melon regions have a lower relative syntenic quality with Arabidopsis (eurosid II) than when compared to Populus and Medicago (eurosid I). Although phylogenetically Cucurbitales seem to be closer to Fabales than to Malphigiales, synteny was higher between both melon regions and Populus. Presented data imply that the recently completed Populus genome sequence could preferentially be used to obtain positional information in melon, based on microsynteny.

Sequence conservation of homeologous bacterial artificial chromosomes and transcription of homeologous genes in soybean (Glycine max L. Merr.)

The paleopolyploid soybean genome was investigated by sequencing homeologous BAC clones anchored by duplicate N-hydroxycinnamoyl/benzoyltransferase (HCBT) genes. The homeologous BACs were genetically mapped to linkage groups C1 and C2. Annotation of the 173,747- and 98,760-bp BACs showed that gene conservation in both order and orientation is high between homeologous regions with only a single gene insertion/deletion and local tandem duplications differing between the regions. The nucleotide sequence conservation extends into intergenic regions as well, probably due to conserved regulatory sequences. Most of the homeologs appear to have a role in either transcription/DNA binding or cellular signaling, suggesting a potential preference for retention of duplicate genes with these functions. Reverse transcriptase-PCR analysis of homeologs showed that in the tissues sampled, most homeologs have not diverged greatly in their transcription profiles. However, four cases of changes in transcription were identified, primarily in the HCBT gene cluster. Because a mapped locus corresponds to a soybean cyst nematode (SCN) QTL, the potential role of HCBT genes in response to SCN is discussed. These results are the first sequenced-based analysis of homeologous BACs in soybean, a diploidized paleopolyploid.

LegumeDB1 bioinformatics resource: comparative genomic analysis and novel cross-genera marker identification in lupin and pasture legume species

The identification of markers in legume pasture crops, which can be associated with traits such as protein and lipid production, disease resistance, and reduced pod shattering, is generally accepted as an important strategy for improving the agronomic performance of these crops. It has been demonstrated that many quantitative trait loci (QTLs) identified in one species can be found in other plant species. Detailed legume comparative genomic analyses can characterize the genome organization between model legume species (e.g., Medicago truncatula, Lotus japonicus) and economically important crops such as soybean (Glycine max), pea (Pisum sativum), chickpea (Cicer arietinum), and lupin (Lupinus angustifolius), thereby identifying candidate gene markers that can be used to track QTLs in lupin and pasture legume breeding. LegumeDB is a Web-based bioinformatics resource for legume researchers. LegumeDB analysis of Medicago truncatula expressed sequence tags (ESTs) has identified novel simple sequence repeat (SSR) markers (16 tested), some of which have been putatively linked to symbiosome membrane proteins in root nodules and cell-wall proteins important in plant-pathogen defence mechanisms. These novel markers by preliminary PCR assays have been detected in Medicago truncatula and detected in at least one other legume species, Lotus japonicus, Glycine max, Cicer arietinum, and (or) Lupinus angustifolius (15/16 tested). Ongoing research has validated some of these markers to map them in a range of legume species that can then be used to compile composite genetic and physical maps. In this paper, we outline the features and capabilities of LegumeDB as an interactive application that provides legume genetic and physical comparative maps, and the efficient feature identification and annotation of the vast tracks of model legume sequences for convenient data integration and visualization. LegumeDB has been used to identify potential novel cross-genera polymorphic legume markers that map to agronomic traits, supporting the accelerated identification of molecular genetic factors underpinning important agronomic attributes in lupin.

Paleopolyploidy and gene duplication in soybean and other legumes

Two of the most important observations from whole-genome sequences have been the high rate of gene birth and death and the prevalence of large-scale duplication events, including polyploidy. There is also a growing appreciation that polyploidy is more than the sum of the gene duplications it creates, in part because polyploidy duplicates the members of entire regulatory networks. Thus, it may be important to distinguish paralogs that are produced by individual gene duplications from the homoeologous sequences produced by (allo)polyploidy. This is not a simple task, for several reasons, including the chromosomally cryptic nature of many duplications and the variable rates of gene evolution. Recent progress has been made in understanding patterns of gene and genome duplication in the legume family, specifically in soybean.

The Soybean Genome Database (SoyGD): a browser for display of duplicated, polyploid, regions and sequence tagged sites on the integrated physical and genetic maps of Glycine max

Genomes that have been highly conserved following increases in ploidy (by duplication or hybridization) like Glycine max (soybean) present challenges during genome analysis. At http://soybeangenome.siu.edu the Soybean Genome Database (SoyGD) genome browser has, since 2002, integrated and served the publicly available soybean physical map, bacterial artificial chromosome (BAC) fingerprint database and genetic map associated genomic data. The browser shows both build 3 and build 4 contiguous sets of clones (contigs) of the soybean physical map. Build 4 consisted of 2854 contigs that encompassed 1.05 Gb and 404 high-quality DNA markers that anchored 742 contigs. Many DNA markers anchored sets of 2-8 different contigs. Each contig in the set represented a homologous region of related sequences. GBrowse was adapted to show sets of homologous contigs at all potential anchor points, spread laterally and prevented from overlapping. About 8064 minimum tiling path (MTP2) clones provided 13,473 BAC end sequences (BES) to decorate the physical map. Analyses of BES placed 2111 gene models, 40 marker anchors and 1053 new microsatellite markers on the map. Estimated sequence tag probes from 201 low-copy gene families located 613 paralogs. The genome browser portal showed each data type as a separate track. Tetraploid, octoploid, diploid and homologous regions are shown clearly in relation to an integrated genetic and physical map.

Significant microsynteny with new evolutionary highlights is detected between Arabidopsis and legume model plants despite the lack of macrosynteny

The increased amount of data produced by large genome sequencing projects allows scientists to carry out important syntenic studies to a great extent. Detailed genetic maps and entirely or partially sequenced genomes are compared, and macro- and microsyntenic relations can be determined for different species. In our study, the syntenic relationships between key legume plants and two model plants, Arabidopsis thaliana and Populus trichocarpa were investigated. The comparison of the map position of 172 gene-based Medicago sativa markers to the organization of homologous A. thaliana genes could not identify any sign of macrosynteny between the two genomes. A 276 kb long section of chromosome 5 of the model legume Medicago truncatula was used to investigate potential microsynteny with the other legume Lotus japonicus, as well as with Arabidopsis and Populus. Besides the overall correlation found between the legume plants, the comparison revealed several microsyntenic regions in the two more distant plants with significant resemblance. Despite the large phylogenetic distance, clear microsyntenic regions between Medicago and Arabidopsis or Populus were detected unraveling new intragenomic evolutionary relations in Arabidopsis.

Highly syntenic regions in the genomes of soybean, Medicago truncatula, and Arabidopsis thaliana

BACKGROUND

Recent genome sequencing enables mega-base scale comparisons between related genomes. Comparisons between animals, plants, fungi, and bacteria demonstrate extensive synteny tempered by rearrangements. Within the legume plant family, glimpses of synteny have also been observed. Characterizing syntenic relationships in legumes is important in transferring knowledge from model legumes to crops that are important sources of protein, fixed nitrogen, and health-promoting compounds.

RESULTS

We have uncovered two large soybean regions exhibiting synteny with M. truncatula and with a network of segmentally duplicated regions in Arabidopsis. In all, syntenic regions comprise over 500 predicted genes spanning 3 Mb. Up to 75% of soybean genes are colinear with M. truncatula, including one region in which 33 of 35 soybean predicted genes with database support are colinear to M. truncatula. In some regions, 60% of soybean genes share colinearity with a network of A. thaliana duplications. One region is especially interesting because this 500 kbp segment of soybean is syntenic to two paralogous regions in M. truncatula on different chromosomes. Phylogenetic analysis of individual genes within these regions demonstrates that one is orthologous to the soybean region, with which it also shows substantially denser synteny and significantly lower levels of synonymous nucleotide substitutions. The other M. truncatula region is inferred to be paralogous, presumably resulting from a duplication event preceding speciation.

CONCLUSION

The presence of well-defined M. truncatula segments showing orthologous and paralogous relationships with soybean allows us to explore the evolution of contiguous genomic regions in the context of ancient genome duplication and speciation events.

Exploratory integration of peanut genetic and physical maps and possible contributions from Arabidopsis

Arachis hypogaea is a widely cultivated crop both as an oilseed and protein source. The genomic analysis of Arachis species hitherto has been limited to the construction of genetic maps; the most comprehensive one contains 370 loci over 2,210 cM in length. However, no attempt has been made to analyze the physical structure of the peanut genome. To investigate the practicality of physical mapping in peanut, we applied a total of 117 oligonucleotide-based probes ("overgos") derived from genetically mapped RFLP probes onto peanut BAC filters containing 182,784 peanut large-insert DNA clones in a multiplex experimental design; 91.5% of the overgos identified at least one BAC clone. In order to gain insights into the potential value of Arabidopsis genome sequence for studies in divergent species with complex genomes such as peanut, we employed 576 Arabidopsis-derived overgos selected on the basis of maximum homology to orthologous sequences in other plant taxa to screen the peanut BAC library. A total of 353 (61.3%) overgos detected at least one peanut BAC clone. This experiment represents the first steps toward the creation of a physical map in peanut and illustrates the potential value of leveraging information from distantly related species such as Arabidopsis for both practical applications such as comparative map-based cloning and shedding light on evolutionary relationships. We also evaluated the possible correlation between functional categories of Arabidopsis overgos and their success rates in hybridization to the peanut BAC library.

Placing Paleopolyploidy in Relation to Taxon Divergence: A Phylogenetic Analysis in Legumes Using 39 Gene Families

Young polyploid events are easily diagnosed by various methods, but older polyploid events become increasingly difficult to identify as chromosomal rearrangements, tandem gene or partial chromosome duplications, changes in substitution rates among duplicated genes, pseudogenization or locus loss, and interlocus interactions complicate the means of inferring past genetic events. Genomic data have provided valuable information about the polyploid history of numerous species, but on their own fail to show whether related species, each with a polyploid past, share a particular polyploid event. A phylogenetic approach provides a powerful method to determine this but many processes may mislead investigators. These processes can affect individual gene trees, but most likely will not affect all genes, and almost certainly will not affect all genes in the same way. Thus, a multigene approach, which combines the large-scale aspect of genomics with the resolution of phylogenetics, has the power to overcome these difficulties and allow us to infer genomic events further into the past than would otherwise be possible. Previous work using synonymous distances among gene pairs within species has shown evidence for large-scale duplications in the legumes Glycine max and Medicago truncatula. We present a case study using 39 gene families, each with three or four members in G. max and the putative orthologues in M. truncatula, rooted using Arabidopsis thaliana. We tested whether the gene duplications in these legumes occurred separately in each lineage after their divergence (Hypothesis 1), or whether they share a round of gene duplications (Hypothesis 2). Many more gene family topologies supported Hypothesis 2 over Hypothesis 1 (11 and 2, respectively), even after synonymous distance analysis revealed that some topologies were providing misleading results. Only ca. 33% of genes examined support either hypothesis, which strongly suggests that single gene family approaches may be insufficient when studying ancient events with nuclear DNA. Our results suggest that G. max and M. truncatula, along with approximately 7000 other legume species from the same clade, share an ancient round of gene duplications, either due to polyploidy or to some other process.

Cross-species amplification of Medicago truncatula microsatellites across three major pulse crops

Model plants are facilitating the genetic characterization and comparative mapping of a number of traditional crops. Medicago truncatula has been widely accepted as a model plant to this end as it provides the essential tools for multiple aspects of legume genetics and genomics. A large set of markers from highly conserved M. truncatula gene regions is being created and used to establish a worldwide framework for comparative genomic studies in legumes. We have investigated the potential for cross-species amplification of 209 expressed sequence tag (EST)-based and 33 bacterial artificial chromosome (BAC)-based microsatellites from M. truncatula in the three most important European legume pulses-pea, faba bean and chickpea-that might facilitate future comparative mapping. Our results revealed significant transferability of M. truncatula microsatellites to the three pulses (40% in faba bean, 36.3% in chickpea and 37.6% in pea). The percentage of M. truncatula EST-SSRs (simple sequence repeats) amplified in the three crops (39-43%) was twofold higher than that of the genomic SSRs (21-24%). Sequence analysis determined that the level of conservation in the microsatellite motif was very low, while the flanking regions were generally well conserved. The variations in the sequences were mainly due to changes in the number of repeat motifs in the microsatellite region combined with indel and base substitutions. None of the functional microsatellites showed direct polymorphism among the parental genotypes tested, consequently preventing their immediate use for mapping purposes.

Genomic analysis of a region encompassing QRfs1 and QRfs2: genes that underlie soybean resistance to sudden death syndrome

Candidate genes were identified for two loci, QRfs2 providing resistance to the leaf scorch called soybean (Glycine max (L.) Merr.) sudden death syndrome (SDS) and QRfs1 providing resistance to root infection by the causal pathogen Fusarium solani f.sp. glycines. The 7.5 +/- 0.5 cM region of chromosome 18 (linkage group G) was shown to encompass a cluster of resistance loci using recombination events from 4 near-isogenic line populations and 9 DNA markers. The DNA markers anchored 9 physical map contigs (7 are shown on the soybean Gbrowse, 2 are unpublished), 45 BAC end sequences (41 in Gbrowse), and contiguous DNA sequences of 315, 127, and 110 kbp. Gene density was high at 1 gene per 7 kbp only around the already sequenced regions. Three to 4 gene-rich islands were inferred to be distributed across the entire 7.5 cM or 3.5 Mbp showing that genes are clustered in the soybean genome. Candidate resistance genes were identified and a molecular basis for interactions among the disease resistance genes in the cluster inferred.

Comparative genome analysis of monocots and dicots, toward characterization of angiosperm diversity

The importance of angiosperms to sustaining humanity by providing a wide range of 'ecosystem services' warrants increased exploration of their genomic diversity. The nearly completed sequences for two species representing the major angiosperm subclasses, specifically the dicot Arabidopsis thaliana and the monocot Oryza sativa, provide a foundation for comparative analysis across the angiosperms. The angiosperms also exemplify some challenges to be faced as genomics makes new inroads into describing biotic diversity, in particular polyploidy (genome-wide chromatin duplication), and much larger genome sizes than have been studied to date.

Estimating genome conservation between crop and model legume species

Legumes are simultaneously one of the largest families of crop plants and a cornerstone in the biological nitrogen cycle. We combined molecular and phylogenetic analyses to evaluate genome conservation both within and between the two major clades of crop legumes. Genetic mapping of orthologous genes identifies broad conservation of genome macrostructure, especially within the galegoid legumes, while also highlighting inferred chromosomal rearrangements that may underlie the variation in chromosome number between these species. As a complement to comparative genetic mapping, we compared sequenced regions of the model legume Medicago truncatula with those of the diploid Lotus japonicus and the polyploid Glycine max. High conservation was observed between the genomes of M. truncatula and L. japonicus, whereas lower levels of conservation were evident between M. truncatula and G. max. In all cases, conserved genome microstructure was punctuated by significant structural divergence, including frequent insertion/deletion of individual genes or groups of genes and lineage-specific expansion/contraction of gene families. These results suggest that comparative mapping may have considerable utility for basic and applied research in the legumes, although its predictive value is likely to be tempered by phylogenetic distance and genome duplication.

Comparative mapping in the Pinaceae

A comparative genetic map was constructed between two important genera of the family Pinaceae. Ten homologous linkage groups in loblolly pine (Pinus taeda L.) and Douglas fir (Pseudotsuga menziesii [Mirb.] Franco) were identified using orthologous expressed sequence tag polymorphism (ESTP) and restriction fragment length polymorphism (RFLP) markers. The comparative mapping revealed extensive synteny and colinearity between genomes of the Pinaceae, consistent with the hypothesis of conservative chromosomal evolution in this important plant family. This study reports the first comparative map in forest trees at the family taxonomic level and establishes a framework for comparative genomics in Pinaceae.

The Iccare web server: an attempt to merge sequence and mapping information for plant and animal species

The Iccare web server, http://genopole.toulouse.inra.fr/bioinfo/Iccare, provides a simple yet efficient tool for crude EST (expressed sequence tag) annotation specifically dedicated to comparative mapping approaches. Iccare uses all the EST and mRNA sequences from public databases for an organism of interest (query species) and compares them to all the transcripts of one reference organism (Homo sapiens or Arabidopsis thaliana). The results are displayed according to the location of the genes on the chromosomes of the reference organism. Gene structure information and sequence similarities are combined in a graphical representation in order to pinpoint the nature of the transcript query sequence. The user can subsequently design primers or probes for the purpose of physical or genetic mapping. In addition to the query organisms already available in Iccare, users can perform a tailor-made search with their own sequences against the animal or plant reference organism genes.

Widespread paleopolyploidy in model plant species inferred from age distributions of duplicate genes

It is often anticipated that many of today's diploid plant species are in fact paleopolyploids. Given that an ancient large-scale duplication will result in an excess of relatively old duplicated genes with similar ages, we analyzed the timing of duplication of pairs of paralogous genes in 14 model plant species. Using EST contigs (unigenes), we identified pairs of paralogous genes in each species and used the level of synonymous nucleotide substitution to estimate the relative ages of gene duplication. For nine of the investigated species (wheat [Triticum aestivum], maize [Zea mays], tetraploid cotton [Gossypium hirsutum], diploid cotton [G. arboretum], tomato [Lycopersicon esculentum], potato [Solanum tuberosum], soybean [Glycine max], barrel medic [Medicago truncatula], and Arabidopsis thaliana), the age distributions of duplicated genes contain peaks corresponding to short evolutionary periods during which large numbers of duplicated genes were accumulated. Large-scale duplications (polyploidy or aneuploidy) are strongly suspected to be the cause of these temporal peaks of gene duplication. However, the unusual age profile of tandem gene duplications in Arabidopsis indicates that other scenarios, such as variation in the rate at which duplicated genes are deleted, must also be considered.

National Science Foundation-sponsored workshop report. Draft plan for soybean genomics

Recent efforts to coordinate and define a research strategy for soybean (Glycine max) genomics began with the establishment of a Soybean Genetics Executive Committee, which will serve as a communication focal point between the soybean research community and granting agencies. Secondly, a workshop was held to define a strategy to incorporate existing tools into a framework for advancing soybean genomics research. This workshop identified and ranked research priorities essential to making more informed decisions as to how to proceed with large scale sequencing and other genomics efforts. Most critical among these was the need to finalize a physical map and to obtain a better understanding of genome microstructure. Addressing these research needs will require pilot work on new technologies to demonstrate an ability to discriminate between recently duplicated regions in the soybean genome and pilot projects to analyze an adequate amount of random genome sequence to identify and catalog common repeats. The development of additional markers, reverse genetics tools, and bioinformatics is also necessary. Successful implementation of these goals will require close coordination among various working groups.

Comparative physical mapping reveals features of microsynteny between Glycine max, Medicago truncatula, and Arabidopsis thaliana

To gain insight into genomic relationships between soybean (Glycine max) and Medicago truncatula, eight groups of bacterial artificial chromosome (BAC) contigs, together spanning 2.60 million base pairs (Mb) in G. max and 1.56 Mb in M. truncatula, were compared through high-resolution physical mapping combined with sequence and hybridization analysis of low-copy BAC ends. Cross-hybridization among G. max and M. truncatula contigs uncovered microsynteny in six of the contig groups and extensive microsynteny in three. Between G. max homoeologous (within genome duplicate) contigs, 85% of coding and 75% of noncoding sequences were conserved at the level of cross-hybridization. By contrast, only 29% of sequences were conserved between G. max and M. truncatula, and some kilobase-scale rearrangements were also observed. Detailed restriction maps were constructed for 11 contigs from the three highly microsyntenic groups, and these maps suggested that sequence order was highly conserved between G. max duplicates and generally conserved between G. max and M. truncatula. One instance of homoeologous BAC contigs in M. truncatula was also observed and examined in detail. A sequence similarity search against the Arabidopsis thaliana genome sequence identified up to three microsyntenic regions in A. thaliana for each of two of the legume BAC contig groups. Together, these results confirm previous predictions of one recent genome-wide duplication in G. max and suggest that M. truncatula also experienced ancient large-scale genome duplications.

Legume genomes: more than peas in a pod

A growing array of sequence-based tools is helping to reveal the organization, evolution and syntenic relationships of legume genomes. The results indicate that legumes form a coherent taxonomic group with frequent and widespread macro- and microsynteny. This is good news for two model legume systems, Medicago truncatula and Lotus japonicus. Indeed, both models have recently been used to clone and characterize genes for nodulation-related receptors that were originally described in legumes with more complex genomes. Studies of legume genomes have also provided insight into genome size, gene clustering, genome duplications and repetitive elements. To understand legume genomes better, it will be necessary to develop tools for studying under-represented taxa beyond the relatively small group of economically important species that have been examined so far.