Comprehensive molecular cytogenetic analysis of sorghum genome architecture: distribution of euchromatin, heterochromatin, genes and recombination in comparison to rice

Cyto-molecular characterization of rDNA and chromatin composition in the NOR-associated satellite in Chestnut (Castanea spp.)

The American chestnut (Castanea dentata, 2n = 2x = 24), once known as the "King of the Appalachian Forest", was decimated by chestnut blight during the first half of the twentieth century by an invasive fungus (Cryphonectria parasitica). The Chinese chestnut (C. mollissima, 2n = 2x = 24), in contrast to American chestnut, is resistant to this blight. Efforts are being made to transfer this resistance to American chestnut through backcross breeding and genetic engineering. Both chestnut genomes have been genetically mapped and recently sequenced to facilitate gene discovery efforts aimed at assisting molecular breeding and genetic engineering. To complement and extend this genomic work, we analyzed the distribution and organization of their ribosomal DNAs (35S and 5S rDNA), and the chromatin composition of the nucleolus organizing region (NOR)-associated satellites. Using fluorescent in situ hybridization (FISH), we have identified two 35S (one major and one minor) and one 5S rDNA sites. The major 35S rDNA sites are terminal and sub-terminal in American and Chinese chestnuts, respectively, originating at the end of the short arm of the chromosome, extending through the secondary constriction and into the satellites. An additional 5S locus was identified in certain Chinese chestnut accessions, and it was linked distally to the major 35S site. The NOR-associated satellite in Chinese chestnut was found to comprise a proximal region packed with 35S rDNA and a distinct distal heterochromatic region. In contrast, the American chestnut satellite was relatively small and devoid of the distal heterochromatic region.

Coffea cytogenetics: from the first karyotypes to the meeting with genomics

Coffea karyotype organization and evolution has been uncovered by classical cytogenetics and cytogenomics. We revisit these discoveries and present new karyotype data. Coffea possesses ~ 124 species, including C. arabica and C. canephora responsible for commercial coffee production. We reviewed the Coffea cytogenetics, from the first chromosome counting, encompassing the karyotype characterization, chromosome DNA content, and mapping of chromosome portions and DNA sequences, until the integration with genomics. We also showed new data about Coffea karyotype. The 2n chromosome number evidenced the diploidy of almost all Coffea, and the C. arabica tetraploidy, as well as the polyploidy of other hybrids. Since then, other genomic similarities and divergences among the Coffea have been shown by karyotype morphology, nuclear and chromosomal C-value, AT and GC rich chromosome portions, and repetitive sequence and gene mapping. These cytogenomic data allowed us to know and understand the phylogenetic relations in Coffea, as well as their ploidy level and genomic origin, highlighting the relatively recent allopolyploidy. In addition to the euploidy, the role of the mobile elements in Coffea diversification is increasingly more evident, and the comparative analysis of their structure and distribution on the genome of different species is in the spotlight for future research. An integrative look at all these data is fundamental for a deeper understanding of Coffea karyotype evolution, including the key role of polyploidy in C. arabica origin. The 'Híbrido de Timor', a recent natural allotriploid, is also in the spotlight for its potential as a source of resistance genes and model for plant polyploidy research. Considering this, we also present some unprecedented results about the exciting evolutionary history of these polyploid Coffea.

Chromosome and Genome Diversity in the Genus Trifolium (Fabaceae)

Trifolium L. is an economically important genus that is characterized by variable karyotypes relating to its ploidy level and basic chromosome numbers. The advent of genomic resources combined with molecular cytogenetics provides an opportunity to develop our understanding of plant genomes in general. Here, we summarize the current state of knowledge on Trifolium genomes and chromosomes and review methodologies using molecular markers that have contributed to Trifolium research. We discuss possible future applications of cytogenetic methods in research on the Trifolium genome and chromosomes.

Genetic characterization of a Sorghum bicolor multiparent mapping population emphasizing carbon-partitioning dynamics

Sorghum bicolor, a photosynthetically efficient C₄ grass, represents an important source of grain, forage, fermentable sugars, and cellulosic fibers that can be utilized in myriad applications ranging from bioenergy to bioindustrial feedstocks. Sorghum’s efficient fixation of carbon per unit time per unit area per unit input has led to its classification as a preferred biomass crop highlighted by its designation as an advanced biofuel by the U.S. Department of Energy. Due to its extensive genetic diversity and worldwide colonization, sorghum has considerable diversity for a range of phenotypes influencing productivity, composition, and sink/source dynamics. To dissect the genetic basis of these key traits, we present a sorghum carbon-partitioning nested association mapping (NAM) population generated by crossing 11 diverse founder lines with Grassl as the single recurrent female. By exploiting existing variation among cellulosic, forage, sweet, and grain sorghum carbon partitioning regimes, the sorghum carbon-partitioning NAM population will allow the identification of important biomass-associated traits, elucidate the genetic architecture underlying carbon partitioning and improve our understanding of the genetic determinants affecting unique phenotypes within Poaceae. We contrast this NAM population with an existing grain population generated using Tx430 as the recurrent female. Genotypic data are assessed for quality by examining variant density, nucleotide diversity, linkage decay, and are validated using pericarp and testa phenotypes to map known genes affecting these phenotypes. We release the 11-family NAM population along with corresponding genomic data for use in genetic, genomic, and agronomic studies with a focus on carbon-partitioning regimes.

The Evolutionary Dynamics of Repetitive DNA and Its Impact on the Genome Diversification in the Genus Sorghum

Polyploidization is an evolutionary event leading to structural changes of the genome(s), particularly allopolyploidization, which combines different genomes of distinct species. The tetraploid species, Sorghum halepense, is assumed an allopolyploid species formed by hybridization between diploid S. bicolor and S. propinquum. The repeat profiles of S. bicolor, S. halepense, and their relatives were compared to elucidate the repeats' role in shaping their genomes. The repeat frequencies and profiles of the three diploid accessions (S. bicolor, S. bicolor ssp. verticilliflorum, and S. bicolor var. technicum) and two tetraploid accessions (S. halepense) are similar. However, the polymorphic distribution of the subtelomeric satellites preferentially enriched in the tetraploid S. halepense indicates drastic genome rearrangements after the allopolyploidization event. Verified by CENH3 chromatin immunoprecipitation (ChIP)-sequencing and fluorescence in situ hybridization (FISH) analysis the centromeres of S. bicolor are mainly composed of the abundant satellite SorSat137 (CEN38) and diverse CRMs, Athila of Ty3_gypsy and Ty1_copia-SIRE long terminal repeat (LTR) retroelements. A similar centromere composition was found in S. halepense. The potential contribution of S. bicolor in the formation of tetraploid S. halepense is discussed.

Amand P, Le HQ, Su Z, Bai G. Multiplex restriction amplicon sequencing: a novel next-generation sequencing-based marker platform for high-throughput genotyping

To enable rapid selection of traits in marker-assisted breeding, markers must be technically simple, low-cost, high-throughput and randomly distributed in a genome. We developed such a technology, designated as Multiplex Restriction Amplicon Sequencing (MRASeq), which reduces genome complexity by polymerase chain reaction (PCR) amplification of amplicons flanked by restriction sites. The first PCR primers contain restriction site sequences at 3'-ends, preceded by 6-10 bases of specific or degenerate nucleotide sequences and then by a unique M13-tail sequence which serves as a binding site for a second PCR that adds sequencing primers and barcodes to allow sample multiplexing for sequencing. The sequences of restriction sites and adjacent nucleotides can be altered to suit different species. Physical mapping of MRASeq SNPs from a biparental population of allohexaploid wheat (Triticum aestivum L.) showed a random distribution of SNPs across the genome. MRASeq generated thousands of SNPs from a wheat biparental population and natural populations of wheat and barley (Hordeum vulgare L.). This novel, next-generation sequencing-based genotyping platform can be used for linkage mapping to screen quantitative trait loci (QTL), background selection in breeding and many other genetics and breeding applications of various species.

High Density Genetic Maps of Seashore Paspalum Using Genotyping-By-Sequencing and Their Relationship to The Sorghum Bicolor Genome

As a step towards trait mapping in the halophyte seashore paspalum (Paspalum vaginatum Sw.), we developed an F₁ mapping population from a cross between two genetically diverse and heterozygous accessions, 509022 and HI33. Progeny were genotyped using a genotyping-by-sequencing (GBS) approach and sequence reads were analyzed for single nucleotide polymorphisms (SNPs) using the UGbS-Flex pipeline. More markers were identified that segregated in the maternal parent (HA maps) compared to the paternal parent (AH maps), suggesting that 509022 had overall higher levels of heterozygosity than HI33. We also generated maps that consisted of markers that were heterozygous in both parents (HH maps). The AH, HA and HH maps each comprised more than 1000 markers. Markers formed 10 linkage groups, corresponding to the ten seashore paspalum chromosomes. Comparative analyses showed that each seashore paspalum chromosome was syntenic to and highly colinear with a single sorghum chromosome. Four inversions were identified, two of which were sorghum-specific while the other two were likely specific to seashore paspalum. These high-density maps are the first available genetic maps for seashore paspalum. The maps will provide a valuable tool for plant breeders and others in the Paspalum community to identify traits of interest, including salt tolerance.

An Integrated Genotyping-by-Sequencing Polymorphism Map for Over 10,000 Sorghum Genotypes

Mining crop genomic variation can facilitate the genetic research of complex traits and molecular breeding. In sorghum [ L. (Moench)], several large-scale single nucleotide polymorphism (SNP) datasets have been generated using genotyping-by-sequencing of KI reduced representation libraries. However, data reuse has been impeded by differences in reference genome coordinates among datasets. To facilitate reuse of these data, we constructed and characterized an integrated 459,304-SNP dataset for 10,323 sorghum genotypes on the version 3.1 reference genome. The SNP distribution showed high enrichment in subtelomeric chromosome arms and in genic regions (48% of SNPs) and was highly correlated ( = 0.82) to the distribution of KI restriction sites. The genetic structure reflected population differences by botanical race, as well as familial structure among recombinant inbred lines (RILs). Faster linkage disequilibrium decay was observed in the diversity panel than in the RILs, as expected, given the greater opportunity for recombination in diverse populations. To validate the quality and utility of the integrated SNP dataset, we used genome-wide association studies (GWAS) of genebank phenotype data, precisely mapping several known genes (e.g and ) and identifying novel associations for other traits. We further validated the dataset with GWAS of new and published plant height and flowering time data in a nested association mapping population, precisely mapping known genes and identifying epistatic interactions underlying both traits. These findings validate this integrated SNP dataset as a useful genomics resource for sorghum genetics and breeding.

The Agropyron cristatum karyotype, chromosome structure and cross-genome homoeology as revealed by fluorescence in situ hybridization with tandem repeats and wheat single-gene probes

Fluorescence in situ hybridization with probes for 45 cDNAs and five tandem repeats revealed homoeologous relationships of Agropyron cristatum with wheat. The results will contribute to alien gene introgression in wheat improvement. Crested wheatgrass (Agropyron cristatum L. Gaertn.) is a wild relative of wheat and a promising source of novel genes for wheat improvement. To date, identification of A. cristatum chromosomes has not been possible, and its molecular karyotype has not been available. Furthermore, homoeologous relationship between the genomes of A. cristatum and wheat has not been determined. To develop chromosome-specific landmarks, A. cristatum genomic DNA was sequenced, and new tandem repeats were discovered. Their distribution on mitotic chromosomes was studied by fluorescence in situ hybridization (FISH), which revealed specific patterns for five repeats in addition to 5S and 45S ribosomal DNA and rye subtelomeric repeats pSc119.2 and pSc200. FISH with one tandem repeat together with 45S rDNA enabled identification of all A. cristatum chromosomes. To analyze the structure and cross-species homoeology of A. cristatum chromosomes with wheat, probes for 45 mapped wheat cDNAs covering all seven chromosome groups were localized by FISH. Thirty-four cDNAs hybridized to homoeologous chromosomes of A. cristatum, nine hybridized to homoeologous and non-homoeologous chromosomes, and two hybridized to unique positions on non-homoeologous chromosomes. FISH using single-gene probes revealed that the wheat-A. cristatum collinearity was distorted, and important structural rearrangements were observed for chromosomes 2P, 4P, 5P, 6P and 7P. Chromosomal inversions were found for pericentric region of 4P and whole chromosome arm 6PL. Furthermore, reciprocal translocations between 2PS and 4PL were detected. These results provide new insights into the genome evolution within Triticeae and will facilitate the use of crested wheatgrass in alien gene introgression into wheat.

B chromosomes are associated with redistribution of genetic recombination towards lower recombination chromosomal regions in perennial ryegrass

Supernumerary 'B' chromosomes are non-essential components of the genome present in a range of plant and animal species-including many grasses. Within diploid and polyploid ryegrass and fescue species, including the forage grass perennial ryegrass (Lolium perenne L.), the presence of B chromosomes has been reported as influencing both chromosome pairing and chiasma frequencies. In this study, the effects of the presence/absence of B chromosomes on genetic recombination has been investigated through generating DArT (Diversity Arrays Technology) marker genetic maps for six perennial ryegrass diploid populations, the pollen parents of which contained either two B or zero B chromosomes. Through genetic and cytological analyses of these progeny and their parents, we have identified that, while overall cytological estimates of chiasma frequencies were significantly lower in pollen mother cells with two B chromosomes as compared with zero B chromosomes, the recombination frequencies within some marker intervals were actually increased, particularly for marker intervals in lower recombination regions of chromosomes, namely pericentromeric regions. Thus, in perennial ryegrass, the presence of two B chromosomes redistributed patterns of meiotic recombination in pollen mother cells in ways which could increase the range of allelic variation available to plant breeders.

Meiotic Crossing Over in Maize Knob Heterochromatin

There is ample evidence that crossing over is suppressed in heterochromatin associated with centromeres and nucleolus organizers (NORs). This characteristic has been attributed to all heterochromatin, but the generalization may not be justified. To investigate the relationship of crossing over to heterochromatin that is not associated with centromeres or NORs, we used a combination of fluorescence in situ hybridization of the maize 180-bp knob repeat to show the locations of knob heterochromatin and fluorescent immunolocalization of MLH1 protein and AFD1 protein to show the locations of MLH1 foci on maize synaptonemal complexes (SCs, pachytene chromosomes). MLH1 foci correspond to the location of recombination nodules (RNs) that mark sites of crossing over. We found that MLH1 foci occur at similar frequencies per unit length of SC in interstitial knobs and in the 1 µm segments of SC in euchromatin immediately to either side of interstitial knobs. These results indicate not only that crossing over occurs within knob heterochromatin, but also that crossing over is not suppressed in the context of SC length in maize knobs. However, because there is more DNA per unit length of SC in knobs compared to euchromatin, crossing over is suppressed (but not eliminated) in knobs in the context of DNA length compared to adjacent euchromatin.

Homologues of potato chromosome 5 show variable collinearity in the euchromatin, but dramatic absence of sequence similarity in the pericentromeric heterochromatin

BACKGROUND

In flowering plants it has been shown that de novo genome assemblies of different species and genera show a significant drop in the proportion of alignable sequence. Within a plant species, however, it is assumed that different haplotypes of the same chromosome align well. In this paper we have compared three de novo assemblies of potato chromosome 5 and report on the sequence variation and the proportion of sequence that can be aligned.

RESULTS

For the diploid potato clone RH89-039-16 (RH) we produced two linkage phase controlled and haplotype-specific assemblies of chromosome 5 based on BAC-by-BAC sequencing, which were aligned to each other and compared to the 52 Mb chromosome 5 reference sequence of the doubled monoploid clone DM 1-3 516 R44 (DM). We identified 17.0 Mb of non-redundant sequence scaffolds derived from euchromatic regions of RH and 38.4 Mb from the pericentromeric heterochromatin. For 32.7 Mb of the RH sequences the correct position and order on chromosome 5 was determined, using genetic markers, fluorescence in situ hybridisation and alignment to the DM reference genome. This ordered fraction of the RH sequences is situated in the euchromatic arms and in the heterochromatin borders. In the euchromatic regions, the sequence collinearity between the three chromosomal homologs is good, but interruption of collinearity occurs at nine gene clusters. Towards and into the heterochromatin borders, absence of collinearity due to structural variation was more extensive and was caused by hemizygous and poorly aligning regions of up to 450 kb in length. In the most central heterochromatin, a total of 22.7 Mb sequence from both RH haplotypes remained unordered. These RH sequences have very few syntenic regions and represent a non-alignable region between the RH and DM heterochromatin haplotypes of chromosome 5.

CONCLUSIONS

Our results show that among homologous potato chromosomes large regions are present with dramatic loss of sequence collinearity. This stresses the need for more de novo reference assemblies in order to capture genome diversity in this crop. The discovery of three highly diverged pericentric heterochromatin haplotypes within one species is a novelty in plant genome analysis. The possible origin and cytogenetic implication of this heterochromatin haplotype diversity are discussed.

PAV markers in Sorghum bicolour: genome pattern, affected genes and pathways, and genetic linkage map construction

KEY MESSAGE

5,511 genic small-size PAVs in sorghum were identified and examined, including the pattern and the function enrichment of PAV genes. 325 PAV markers were developed to construct a genetic map. Presence/absence variants (PAVs) correlate closely to the phenotypic variation, by impacting plant genome sizes and the adaption to the environment. To shed more light on their genome-wide patterns, functions and the possibility of using them as molecular markers, we generated next generation genome sequencing data for four sorghum inbred lines and used associated bioinformatic pipelines to identify small-size PAVs (40-10 kb). Five thousand five hundreds and eleven genic PAVs (40-10 kb) were identified and found to affect 3,238 genes. These PAVs were mainly distributed on the sub-telomeric regions, but the highest proportions occurred in the vicinity of the centromeric regions. One of the prominent features of the PAVs is the high occurrence of long terminal repeats retrotransposons and DNA transposons. PAVs caused various alterations to gene structure, primarily including the coding sequence variants, intron variants, transcript ablation, and initiator codon changes. The genes affected by PAVs were significantly enriched in those involved in stress responses and protein modification. We used 325 PAVs polymorphic between two sorghum inbred lines Ji2731 and E-Tian, together with 49 SSR markers, and constructed a genetic map, which consisted of 10 linkage groups corresponding to the 10 chromosomes of sorghum and spanned 1,430.3 cM in length covering 97% of the physical genome. The resources reported here should be useful for genetic study and breeding of sorghum and related species.

Construction of cytogenetic map of Gossypium herbaceum chromosome 1 and its integration with genetic maps

Background

Cytogenetic map can provide not only information of the genome structure, but also can build a solid foundation for genetic research. With the developments of molecular and cytogenetic studies in cotton (Gossypium), the construction of cytogenetic map is becoming more and more imperative.

Results

A cytogenetic map of chromosome 1 (A₁01) of Gossypium herbaceum (A₁) which includes 10 bacterial artificial chromosome (BAC) clones was constructed by using fluorescent in situ hybridization (FISH). Meanwhile, comparison and analysis were made for the cytogenetic map of chromosome 1 (A₁01) of G. herbaceum with four genetic linkage maps of chromosome 1 (A_h01) of G. hirsutum ((AD)₁) and one genetic linkage map of chromosome 1 of (A₁01) G. arboreum (A₂). The 10 BAC clones were also used to be localized on G. raimondii (D₅) chromosome 1 (D₅01), and 2 of them showed clear unique hybridized signals. Furthermore, these 2 BAC clones were also shown localized on chromosome 1 of both A sub-genome and D sub-genome of G. hirsutum.

Conclusion

The comparison of the cytogenetic map with genetic linkage maps showed that most of the identified marker-tagged BAC clones appearing same orders in different maps except three markers showing different positions, which might indicate chromosomal segmental rearrangements. The positions of the 2 BAC clones which were localized on A_h01 and D_h01 chromosomes were almost the same as that on A₁01 and D₅01 chromosomes. The corresponding anchored SSR markers of these 2 BAC clones were firstly found to be localized on chromosome D₅01 (D_h01) as they were not seen mapped like this in any genetic map reported.

The low-recombining pericentromeric region of barley restricts gene diversity and evolution but not gene expression

The low-recombining pericentromeric region of the barley genome contains roughly a quarter of the genes of the species, embedded in low-recombining DNA that is rich in repeats and repressive chromatin signatures. We have investigated the effects of pericentromeric region residency upon the expression, diversity and evolution of these genes. We observe no significant difference in average transcript level or developmental RNA specificity between the barley pericentromeric region and the rest of the genome. In contrast, all of the evolutionary parameters studied here show evidence of compromised gene evolution in this region. First, genes within the pericentromeric region of wild barley show reduced diversity and significantly weakened purifying selection compared with the rest of the genome. Second, gene duplicates (ohnolog pairs) derived from the cereal whole-genome duplication event ca. 60MYa have been completely eliminated from the barley pericentromeric region. Third, local gene duplication in the pericentromeric region is reduced by 29% relative to the rest of the genome. Thus, the pericentromeric region of barley is a permissive environment for gene expression but has restricted gene evolution in a sizeable fraction of barley's genes.

Energy sorghum--a genetic model for the design of C4 grass bioenergy crops

Sorghum is emerging as an excellent genetic model for the design of C4 grass bioenergy crops. Annual energy Sorghum hybrids also serve as a source of biomass for bioenergy production. Elucidation of Sorghum's flowering time gene regulatory network, and identification of complementary alleles for photoperiod sensitivity, enabled large-scale generation of energy Sorghum hybrids for testing and commercial use. Energy Sorghum hybrids with long vegetative growth phases were found to accumulate more than twice as much biomass as grain Sorghum, owing to extended growing seasons, greater light interception, and higher radiation use efficiency. High biomass yield, efficient nitrogen recycling, and preferential accumulation of stem biomass with low nitrogen content contributed to energy Sorghum's elevated nitrogen use efficiency. Sorghum's integrated genetics-genomics-breeding platform, diverse germplasm, and the opportunity for annual testing of new genetic designs in controlled environments and in multiple field locations is aiding fundamental discovery, and accelerating the improvement of biomass yield and optimization of composition for biofuels production. Recent advances in wide hybridization between Sorghum and other C4 grasses could allow the deployment of improved genetic designs of annual energy Sorghums in the form of wide-hybrid perennial crops. The current trajectory of energy Sorghum genetic improvement indicates that it will be possible to sustainably produce biofuels from C4 grass bioenergy crops that are cost competitive with petroleum-based transportation fuels.

Variation block-based genomics method for crop plants

BACKGROUND

In contrast with wild species, cultivated crop genomes consist of reshuffled recombination blocks, which occurred by crossing and selection processes. Accordingly, recombination block-based genomics analysis can be an effective approach for the screening of target loci for agricultural traits.

RESULTS

We propose the variation block method, which is a three-step process for recombination block detection and comparison. The first step is to detect variations by comparing the short-read DNA sequences of the cultivar to the reference genome of the target crop. Next, sequence blocks with variation patterns are examined and defined. The boundaries between the variation-containing sequence blocks are regarded as recombination sites. All the assumed recombination sites in the cultivar set are used to split the genomes, and the resulting sequence regions are termed variation blocks. Finally, the genomes are compared using the variation blocks. The variation block method identified recurring recombination blocks accurately and successfully represented block-level diversities in the publicly available genomes of 31 soybean and 23 rice accessions. The practicality of this approach was demonstrated by the identification of a putative locus determining soybean hilum color.

CONCLUSIONS

We suggest that the variation block method is an efficient genomics method for the recombination block-level comparison of crop genomes. We expect that this method will facilitate the development of crop genomics by bringing genomics technologies to the field of crop breeding.

A high-resolution genetic map of yellow monkeyflower identifies chemical defense QTLs and recombination rate variation

Genotyping-by-sequencing methods have vastly improved the resolution and accuracy of genetic linkage maps by increasing both the number of marker loci as well as the number of individuals genotyped at these loci. Using restriction-associated DNA sequencing, we construct a dense linkage map for a panel of recombinant inbred lines derived from a cross between divergent ecotypes of Mimulus guttatus. We used this map to estimate recombination rate across the genome and to identify quantitative trait loci for the production of several secondary compounds (PPGs) of the phenylpropanoid pathway implicated in defense against herbivores. Levels of different PPGs are correlated across recombinant inbred lines suggesting joint regulation of the phenylpropanoid pathway. However, the three quantitative trait loci identified in this study each act on a distinct PPG. Finally, we map three putative genomic inversions differentiating the two parental populations, including a previously characterized inversion that contributes to life-history differences between the annual/perennial ecotypes.

A comprehensive genetic map of sugarcane that provides enhanced map coverage and integrates high-throughput Diversity Array Technology (DArT) markers

BACKGROUND

Sugarcane genetic mapping has lagged behind other crops due to its complex autopolyploid genome structure. Modern sugarcane cultivars have from 110-120 chromosomes and are in general interspecific hybrids between two species with different basic chromosome numbers: Saccharum officinarum (2n = 80) with a basic chromosome number of 10 and S. spontaneum (2n = 40-128) with a basic chromosome number of 8. The first maps that were constructed utilised the single dose (SD) markers generated using RFLP, more recent maps generated using AFLP and SSRs provided at most 60% genome coverage. Diversity Array Technology (DArT) markers are high throughput allowing greater numbers of markers to be generated.

RESULTS

Progeny from a cross between a sugarcane variety Q165 and a S. officinarum accession IJ76-514 were used to generate 2467 SD markers. A genetic map of Q165 was generated containing 2267 markers, These markers formed 160 linkage groups (LGs) of which 147 could be placed using allelic information into the eight basic homology groups (HGs) of sugarcane. The HGs contained from 13 to 23 LGs and from 204 to 475 markers with a total map length of 9774.4 cM and an average density of one marker every 4.3 cM. Each homology group contained on average 280 markers of which 43% were DArT markers 31% AFLP, 16% SSRs and 6% SNP markers. The multi-allelic SSR and SNP markers were used to place the LGs into HGs.

CONCLUSIONS

The DArT array has allowed us to generate and map a larger number of markers than ever before and consequently to map a larger portion of the sugarcane genome. This larger number of markers has enabled 92% of the LGs to be placed into the 8 HGs that represent the basic chromosome number of the ancestral species, S. spontaneum. There were two HGs (HG2 and 8) that contained larger numbers of LGs verifying the alignment of two sets of S. officinarum chromosomes with one set of S. spontaneum chromosomes and explaining the difference in basic chromosome number between the two ancestral species. There was also evidence of more complex structural differences between the two ancestral species.

A genome-wide association study of seed protein and oil content in soybean

BACKGROUND

Association analysis is an alternative to conventional family-based methods to detect the location of gene(s) or quantitative trait loci (QTL) and provides relatively high resolution in terms of defining the genome position of a gene or QTL. Seed protein and oil concentration are quantitative traits which are determined by the interaction among many genes with small to moderate genetic effects and their interaction with the environment. In this study, a genome-wide association study (GWAS) was performed to identify quantitative trait loci (QTL) controlling seed protein and oil concentration in 298 soybean germplasm accessions exhibiting a wide range of seed protein and oil content.

RESULTS

A total of 55,159 single nucleotide polymorphisms (SNPs) were genotyped using various methods including Illumina Infinium and GoldenGate assays and 31,954 markers with minor allele frequency >0.10 were used to estimate linkage disequilibrium (LD) in heterochromatic and euchromatic regions. In euchromatic regions, the mean LD (r2) rapidly declined to 0.2 within 360 Kbp, whereas the mean LD declined to 0.2 at 9,600 Kbp in heterochromatic regions. The GWAS results identified 40 SNPs in 17 different genomic regions significantly associated with seed protein. Of these, the five SNPs with the highest associations and seven adjacent SNPs were located in the 27.6-30.0 Mbp region of Gm20. A major seed protein QTL has been previously mapped to the same location and potential candidate genes have recently been identified in this region. The GWAS results also detected 25 SNPs in 13 different genomic regions associated with seed oil. Of these markers, seven SNPs had a significant association with both protein and oil.

CONCLUSIONS

This research indicated that GWAS not only identified most of the previously reported QTL controlling seed protein and oil, but also resulted in narrower genomic regions than the regions reported as containing these QTL. The narrower GWAS-defined genome regions will allow more precise marker-assisted allele selection and will expedite positional cloning of the causal gene(s).

High-throughput genomics in sorghum: from whole-genome resequencing to a SNP screening array

With its small, diploid and completely sequenced genome, sorghum (Sorghum bicolor L. Moench) is highly amenable to genomics-based breeding approaches. Here, we describe the development and testing of a robust single-nucleotide polymorphism (SNP) array platform that enables polymorphism screening for genome-wide and trait-linked polymorphisms in genetically diverse S. bicolor populations. Whole-genome sequences with 6× to 12× coverage from five genetically diverse S. bicolor genotypes, including three sweet sorghums and two grain sorghums, were aligned to the sorghum reference genome. From over 1 million high-quality SNPs, we selected 2124 Infinium Type II SNPs that were informative in all six source genomes, gave an optimal Assay Design Tool (ADT) score, had allele frequencies of 50% in the six genotypes and were evenly spaced throughout the S. bicolor genome. Furthermore, by phenotype-based pool sequencing, we selected an additional 876 SNPs with a phenotypic association to early-stage chilling tolerance, a key trait for European sorghum breeding. The 3000 attempted bead types were used to populate half of a dual-species Illumina iSelect SNP array. The array was tested using 564 Sorghum spp. genotypes, including offspring from four unrelated recombinant inbred line (RIL) and F2 populations and a genetic diversity collection. A high call rate of over 80% enabled validation of 2620 robust and polymorphic sorghum SNPs, underlining the efficiency of the array development scheme for whole-genome SNP selection and screening, with diverse applications including genetic mapping, genome-wide association studies and genomic selection.

Extensive variation in the density and distribution of DNA polymorphism in sorghum genomes

Sorghum genotypes currently used for grain production in the United States were developed from African landraces that were imported starting in the mid-to-late 19^th century. Farmers and plant breeders selected genotypes for grain production with reduced plant height, early flowering, increased grain yield, adaptation to drought, and improved resistance to lodging, diseases and pests. DNA polymorphisms that distinguish three historically important grain sorghum genotypes, BTx623, BTx642 and Tx7000, were characterized by genome sequencing, genotyping by sequencing, genetic mapping, and pedigree-based haplotype analysis. The distribution and density of DNA polymorphisms in the sequenced genomes varied widely, in part because the lines were derived through breeding and selection from diverse Kafir, Durra, and Caudatum race accessions. Genomic DNA spanning dw1 (SBI-09) and dw3 (SBI-07) had identical haplotypes due to selection for reduced height. Lower SNP density in genes located in pericentromeric regions compared with genes located in euchromatic regions is consistent with background selection in these regions of low recombination. SNP density was higher in euchromatic DNA and varied >100-fold in contiguous intervals that spanned up to 300 Kbp. The localized variation in DNA polymorphism density occurred throughout euchromatic regions where recombination is elevated, however, polymorphism density was not correlated with gene density or DNA methylation. Overall, sorghum chromosomes contain distal euchromatic regions characterized by extensive, localized variation in DNA polymorphism density, and large pericentromeric regions of low gene density, diversity, and recombination.

A high-resolution cucumber cytogenetic map integrated with the genome assembly

Background

High-resolution cytogenetic map can provide not only important biological information on genome organization but also solid foundation for genetic and genomic research. The progress in the molecular and cytogenetic studies has created the basis for developing the cytogenetic map in cucumber (Cucumis sativus L.).

Results

Here, the cytogenetic maps of four cucumber chromosomes (chromosomes 1, 3–5) were constructed by fluorescence in situ hybridization (FISH) analysis on cucumber pachytene chromosomes. Together with our previously constructed cytogenetic maps of three cucumber chromosomes (chromosomes 2, 6–7), cucumber has a complete cytogenetic map with 76 anchoring points between the genetic, the cytogenetic and the draft genome assembly maps. To compare our pachytene FISH map directly to the genetic linkage and draft genome assembly maps, we used a standardized map unit—relative map position (RMP) to produce the comparative map alignments. The alignments allowed a global view of the relationship of genetic and physical distances along each cucumber chromosome, and accuracy and coverage of the draft genome assembly map.

Conclusions

We demonstrated a good correlation between positions of the markers in the linkage and physical maps, and essentially complete coverage of chromosome arms by the draft genome assembly. Our study not only provides essential information for the improvement of sequence assembly but also offers molecular tools for cucumber genomics research, comparative genomics and evolutionary study.

Digital genotyping of sorghum - a diverse plant species with a large repeat-rich genome

Background

Rapid acquisition of accurate genotyping information is essential for all genetic marker-based studies. For species with relatively small genomes, complete genome resequencing is a feasible approach for genotyping; however, for species with large and highly repetitive genomes, the acquisition of whole genome sequences for the purpose of genotyping is still relatively inefficient and too expensive to be carried out on a high-throughput basis. Sorghum bicolor is a C₄ grass with a sequenced genome size of ~730 Mb, of which ~80% is highly repetitive. We have developed a restriction enzyme targeted genome resequencing method for genetic analysis, termed Digital Genotyping (DG), to be applied to sorghum and other grass species with large repeat-rich genomes.

Results

DG templates are generated using one of three methylation sensitive restriction enzymes that recognize a nested set of 4, 6 or 8 bp GC-rich sequences, enabling varying depth of analysis and integration of results among assays. Variation in sequencing efficiency among DG markers was correlated with template GC-content and length. The expected DG allele sequence was obtained 97.3% of the time with a ratio of expected to alternative allele sequence acquisition of >20:1. A genetic map aligned to the sorghum genome sequence with an average resolution of 1.47 cM was constructed using 1,772 DG markers from 137 recombinant inbred lines. The DG map enhanced the detection of QTL for variation in plant height and precisely aligned QTL such as Dw3 to underlying genes/alleles. Higher-resolution NgoMIV-based DG haplotypes were used to trace the origin of DNA on SBI-06, spanning Ma1 and Dw2 from progenitors to BTx623 and IS3620C. DG marker analysis identified the correct location of two miss-assembled regions and located seven super contigs in the sorghum reference genome sequence.

Conclusion

DG technology provides a cost-effective approach to rapidly generate accurate genotyping data in sorghum. Currently, data derived from DG are used for many marker-based analyses, including marker-assisted breeding, pedigree and QTL analysis, genetic map construction, map-based gene cloning and association studies. DG in combination with whole genome resequencing is dramatically accelerating all aspects of genetic analysis of sorghum, an important genetic reference for C₄ grass species.

Individual chromosome identification, chromosomal collinearity and genetic-physical integrated map in Gossypium darwinii and four D genome cotton species revealed by BAC-FISH

The study was conducted on the individual chromosome identification in Gossypium darwinii (A(d)D(d)), G. klotzschianum (D(3k)), G. davidsonii (D(3d)), G. armourianum (D(2-1)) and G. aridum (D(4)) using a multi-probe fluorescence of in situ hybridization (FISH) system. Comparative analysis on their genetic maps with that of physical maps was made as well. The FISH probes used contained two sets of bacterial artificial chromosome (BAC) clones [one is specific to 26 individual chromosomes from A and D subgenomes of G. hirsutum (A(h) and D(h)) while the other is a D genome centromere-specific BAC clone 150D24], 45S and 5S rDNA clones. The results showed that all A(d) chromosomes were marked with the 13 A(h) chromosome-specific BAC clones, whilst all D(d), D(3k), D(3d), D(2-1) and D(4) chromosomes and chromosomal arms were identified with the 13 D(h) chromosome-specific BAC clones and the D genome centromere-specific BAC. According to the homology within D subgenomes which are between A (D) genome and A (D) subgenome, the systematic nomenclature for individual chromosome in the five species was established. The chromosomes of A (D) subgenomes in G. darwinii were named as A(d)01-A(d)13 (D(d)01-D(d)13). The chromosomes in D(3k), D(3d), D(2-1) and D(4) were named as D(3k)01-D(3k)13, D(3d)01-D(3d)13, D(2-1)01-D(2-1)13 and D(4)01-D(4)13, respectively. Based on the successful identification for individual chromosomes, 45S and 5S rDNA were located to the special chromosomes and chromosomal arms for all five species. And there appeared chromosomal collinearity from the BAC clones among different species by comparing BACs positions, which suggested that the majority of chromosome segment homology may exist between D genomes and D subgenome. Moreover, as the genetic map and physical map were integrated, the orientations of genetic maps for D(d) and D genomes of diploid cotton were established. The orientations of some of chromosomes in genetic maps (D(d)03, D(d)04, D(d)06, D(d)09, D(d)10 and D(d)12) were found switched. The SSR marker in the middle of linkage group 04 was corrected at nearby the end of chromosome 04 by FISH. The study will be helpful to establish a theoretical basis using the wild gene bank to exploit more genes aiming for cotton breeding and will provide powerful evidences both for the evolution of Gossypium and assembling the sequences of the obtained and as well the on-going whole genome sequencing projects of cotton.

Genetic structure and linkage disequilibrium in a diverse, representative collection of the C4 model plant, Sorghum bicolor

To facilitate the mapping of genes in sorghum [Sorghum bicolor (L.) Moench] underlying economically important traits, we analyzed the genetic structure and linkage disequilibrium in a sorghum mini core collection of 242 landraces with 13,390 single-nucleotide polymorphims. The single-nucleotide polymorphisms were produced using a highly multiplexed genotyping-by-sequencing methodology. Genetic structure was established using principal component, Neighbor-Joining phylogenetic, and Bayesian cluster analyses. These analyses indicated that the mini-core collection was structured along both geographic origin and sorghum race classification. Examples of the former were accessions from Southern Africa, East Asia, and Yemen. Examples of the latter were caudatums with widespread geographical distribution, durras from India, and guineas from West Africa. Race bicolor, the most primitive and the least clearly defined sorghum race, clustered among other races and formed only one clear bicolor-centric cluster. Genome-wide linkage disequilibrium analyses showed linkage disequilibrium decayed, on average, within 10-30 kb, whereas the short arm of SBI-06 contained a linkage disequilibrium block of 20.33 Mb, confirming a previous report of low recombination on this chromosome arm. Four smaller but equally significant linkage disequilibrium blocks of 3.5-35.5 kb were detected on chromosomes 1, 2, 9, and 10. We examined the genes encoded within each block to provide a first look at candidates such as homologs of GS3 and FT that may indicate a selective sweep during sorghum domestication.

Integration of high-resolution physical and genetic map reveals differential recombination frequency between chromosomes and the genome assembling quality in cucumber

Cucumber is an important model crop and the first species sequenced in Cucurbitaceae family. Compared to the fast increasing genetic and genomics resources, the molecular cytogenetic researches in cucumber are still very limited, which results in directly the shortage of relation between plenty of physical sequences or genetic data and chromosome structure. We mapped twenty-three fosmids anchored by SSR markers from LG-3, the longest linkage group, and LG-4, the shortest linkage group on pachytene chromosomes 3 and 4, using uorescence in situ hybridization (FISH). Integrated molecular cytogenetic maps of chromosomes 3 and 4 were constructed. Except for three SSR markers located on heterochromatin region, the cytological order of markers was concordant with those on the linkage maps. Distinct structural differences between chromosomes 3 and 4 were revealed by the high resolution pachytene chromosomes. The extreme difference of genetic length between LG-3 and LG-4 was mainly attributed to the difference of overall recombination frequency. The significant differentiation of heterochromatin contents in chromosomes 3 and 4 might have a direct correlation with recombination frequency. Meanwhile, the uneven distribution of recombination frequency along chromosome 4 was observed, and recombination frequency of the long arm was nearly 3.5 times higher than that of the short arm. The severe suppression of recombination was exhibited in centromeric and heterochromatin domains of chromosome 4. Whereas a close correlation between the gene density and recombination frequency was observed in chromosome 4, no significant correlation was observed between them along chromosome 3. The comparison between cytogenetic and sequence maps revealed a large gap on the pericentromeric heterochromatin region of sequence map of chromosome 4. These results showed that integrated molecular cytogenetic maps can provide important information for the study of genetic and genomics in cucumber.

A sugar beet (Beta vulgaris L.) reference FISH karyotype for chromosome and chromosome-arm identification, integration of genetic linkage groups and analysis of major repeat family distribution

We developed a reference karyotype for B. vulgaris which is applicable to all beet cultivars and provides a consistent numbering of chromosomes and genetic linkage groups. Linkage groups of sugar beet were assigned to physical chromosome arms by FISH (fluorescent in situ hybridization) using a set of 18 genetically anchored BAC (bacterial artificial chromosome) markers. Genetic maps of sugar beet were correlated to chromosome arms, and North-South orientation of linkage groups was established. The FISH karyotype provides a technical platform for genome studies and can be applied for numbering and identification of chromosomes in related wild beet species. The discrimination of all nine chromosomes by BAC probes enabled the study of chromosome-specific distribution of the major repetitive components of sugar beet genome comprising pericentromeric, intercalary and subtelomeric satellites and 18S-5.8S-25S and 5S rRNA gene arrays. We developed a multicolor FISH procedure allowing the identification of all nine sugar beet chromosome pairs in a single hybridization using a pool of satellite DNA probes. Fiber-FISH was applied to analyse five chromosome arms in which the furthermost genetic marker of the linkage group was mapped adjacently to terminal repetitive sequences on pachytene chromosomes. Only on two arms telomere arrays and the markers are physically linked, hence these linkage groups can be considered as terminally closed making the further identification of distal informative markers difficult. The results support genetic mapping by marker localization, the anchoring of contigs and scaffolds for the annotation of the sugar beet genome sequence and the analysis of the chromosomal distribution patterns of major families of repetitive DNA.

An integrated cytogenetic and physical map reveals unevenly distributed recombination spots along the papaya sex chromosomes

Papaya is a model system for the study of sex chromosome evolution in plants. However, the cytological structures of the papaya chromosomes remain largely unknown and chromosomal features have not been linked with any genetic or genomic data. We constructed a cytogenetic map of the papaya sex chromosome (chromosome 1) by hybridizing 16 microsatellite markers and 2 cytological feature-associated markers on pachytene chromosomes using fluorescence in situ hybridization (FISH). Except for three markers, the order of the markers was concordant to that of marker loci along the linkage map. This discrepancy was likely caused by skewed segregation in the highly heterochromatic or centromeric regions. The papaya sex chromosome is largely euchromatic, its heterochromatin spans about 15 % of the Y chromosome and is mostly restricted to the centromeric and pericentromeric regions. Analysis of the recombination frequency along the papaya sex chromosome revealed a complete suppression of recombination in the centromere and pericentromere region and 60 % higher recombination rate in the long arm than in the short arm. The uneven distribution of recombination events might be caused by differences in sequence composition. Sequence analysis of 18 scaffolds in total length of 15 Mb revealed higher gene density towards the telomeres and lower gene density towards the centromere, and a relatively higher gene density in the long arm than in the short arm. In an opposite trend, the centromeric and pericentromeric region contained the highest repetitive sequences and the long arm showed the lowest repetitive sequences. This cytogenetic map provides essential information for evolutionary study of sex chromosomes in Caricaceae and will facilitate the analysis of papaya sex chromosomes.

Preparations of meiotic pachytene chromosomes and extended DNA fibers from cotton suitable for fluorescence in situ hybridization

Fluorescence in situ hybridization (FISH) has become one of the most important techniques applied in plant molecular cytogenetics. However, the application of this technique in cotton has lagged behind because of difficulties in chromosome preparation. The focus of this article was FISH performed not only on cotton pachytene chromosomes, but also on cotton extended DNA fibers. The cotton pollen mother cells (PMCs) instead of buds or anthers were directly digested in enzyme to completely breakdown the cell wall. Before the routine acetic acid treatment, PMCs were incubated in acetic acid and enzyme mixture to remove the cytoplasm and clear the background. The method of ice-cold Carnoy's solution spreading chromosome was adopted instead of nitrogen removed method to avoid chromosomes losing and fully stretch chromosome. With the above-improved steps, the high-quality well-differentiated pachytene chromosomes with clear background were obtained. FISH results demonstrated that a mature protocol of cotton pachytene chromosomes preparation was presented. Intact and no debris cotton nuclei were obtained by chopping from etiolation cotyledons instead of the conventional liquid nitrogen grinding method. After incubating the nuclei with nucleus lysis buffer on slide, the parallel and clear background DNA fibers were acquired along the slide. This method overcomes the twist, accumulation and fracture of DNA fibers compared with other methods. The entire process of DNA fibers preparation requires only 30 min, in contrast, it takes 3 h with routine nitrogen grinding method. The poisonous mercaptoethanol in nucleus lysis buffer is replaced by nonpoisonous dithiothreitol. PVP40 in nucleus isolation buffer is used to prevent oxidation. The probability of success in isolating nuclei for DNA fiber preparation is almost 100% tested with this method in cotton. So a rapid, safe, and efficient method for the preparation of cotton extended DNA fibers suitable for FISH was established.

Integrating sorghum whole genome sequence information with a compendium of sorghum QTL studies reveals uneven distribution of QTL and of gene-rich regions with significant implications for crop improvement

A comprehensive analysis was conducted using 48 sorghum QTL studies published from 1995 to 2010 to make information from historical sorghum QTL experiments available in a form that could be more readily used by sorghum researchers and plant breeders. In total, 771 QTL relating to 161 unique traits from 44 studies were projected onto a sorghum consensus map. Confidence intervals (CI) of QTL were estimated so that valid comparisons could be made between studies. The method accounted for the number of lines used and the phenotypic variation explained by individual QTL from each study. In addition, estimated centimorgan (cM) locations were calculated for the predicted sorghum gene models identified in Phytozome (JGI GeneModels SBI v1.4) and compared with QTL distribution genome-wide, both on genetic linkage (cM) and physical (base-pair/bp) map scales. QTL and genes were distributed unevenly across the genome. Heterochromatic enrichment for QTL was observed, with approximately 22% of QTL either entirely or partially located in the heterochromatic regions. Heterochromatic gene enrichment was also observed based on their predicted cM locations on the sorghum consensus map, due to suppressed recombination in heterochromatic regions, in contrast to the euchromatic gene enrichment observed on the physical, sequence-based map. The finding of high gene density in recombination-poor regions, coupled with the association with increased QTL density, has implications for the development of more efficient breeding systems in sorghum to better exploit heterosis. The projected QTL information described, combined with the physical locations of sorghum sequence-based markers and predicted gene models, provides sorghum researchers with a useful resource for more detailed analysis of traits and development of efficient marker-assisted breeding strategies.

An integrated molecular cytogenetic map of Cucumis sativus L. chromosome 2

Background

Integration of molecular, genetic and cytological maps is still a challenge for most plant species. Recent progress in molecular and cytogenetic studies created a basis for developing integrated maps in cucumber (Cucumis sativus L.).

Results

In this study, eleven fosmid clones and three plasmids containing 45S rDNA, the centromeric satellite repeat Type III and the pericentriomeric repeat CsRP1 sequences respectively were hybridized to cucumber metaphase chromosomes to assign their cytological location on chromosome 2. Moreover, an integrated molecular cytogenetic map of cucumber chromosomes 2 was constructed by fluorescence in situ hybridization (FISH) mapping of 11 fosmid clones together with the cucumber centromere-specific Type III sequence on meiotic pachytene chromosomes. The cytogenetic map was fully integrated with genetic linkage map since each fosmid clone was anchored by a genetically mapped simple sequence repeat marker (SSR). The relationship between the genetic and physical distances along chromosome was analyzed.

Conclusions

Recombination was not evenly distributed along the physical length of chromosome 2. Suppression of recombination was found in centromeric and pericentromeric regions. Our results also indicated that the molecular markers composing the linkage map for chromosome 2 provided excellent coverage of the chromosome.

Variation in crossover rates across a 3-Mb contig of bread wheat (Triticum aestivum) reveals the presence of a meiotic recombination hotspot

In bread wheat (Triticum aestivum L.), initial studies using deletion lines indicated that crossover (CO) events occur mainly in the telomeric regions of the chromosomes with a possible correlation with the presence of genes. However, little is known about the distribution of COs at the sequence level. To investigate this, we studied in detail the pattern of COs along a contig of 3.110 Mb using two F2 segregating populations (Chinese Spring × Renan (F2-CsRe) and Chinese Spring × Courtot (F2-CsCt)) each containing ~2,000 individuals. The availability of the sequence of the contig from Cs enabled the development of 318 markers among which 23 co-dominant polymorphic markers (11 SSRs and 12 SNPs) were selected for CO distribution analyses. The distribution of CO events was not homogeneous throughout the contig, ranging from 0.05 to 2.77 cM/Mb, but was conserved between the two populations despite very different contig recombination rate averages (0.82 cM/Mb in F2-CsRe vs 0.35 cM/Mb in F2-CsCt). The CO frequency was correlated with the percentage of coding sequence in Cs and with the polymorphism rate between Cs and Re or Ct in both populations, indicating an impact of these two factors on CO distribution. At a finer scale, COs were found in a region covering 2.38 kb, spanning a gene coding for a glycosyl transferase (Hga3), suggesting the presence of a CO hotspot. A non-crossover event covering at least 453 bp was also identified in the same interval. From these results, we can conclude that gene content could be one of the factors driving recombination in bread wheat.

Early-generation germplasm introgression from Sorghum macrospermum into sorghum (S. bicolor)

Sorghum has been improved by public and private breeding programs utilizing germplasm mostly from within the species Sorghum bicolor . Recently, hybridization with an Australian species, S. macrospermum (AAB1B1YYZZ), has been demonstrated and the genomic relationship to S. bicolor (AAB1B1) shown to be partially compatible. For this species to be potentially useful to sorghum improvement programs, there must be documented introgression into an S. bicolor background. Fifteen BC1F1 progeny were recovered using the interspecific hybrid as a female and embryo rescue. In these progeny, chromosome numbers ranged from 35 to 70 and all were essentially male-sterile. Repeated backcrossing with S. bicolor pollen produced BC2F1 seed on 3 of the 15 BC1F1 plants. BC2F1 progeny had varying levels of male fertility; selfed seed set ranged from 0% to 95%, with only 2 individuals being completely male-sterile. Using AFLP and SSR markers, genomic introgression of S. macrospermum ranged from 0% to 18.6%. Cytogenetic analysis of 19 individuals revealed that chromosome numbers were 20, except for a single backcross that had 21 chromosomes. Molecular cytogenetic analysis confirmed the presence of recombinant introgression chromosomes as well as alien addition and alien substitution chromosomes within the BC2F1s.

Cytogenetic map of common bean (Phaseolus vulgaris L.)

A cytogenetic map of common bean was built by in situ hybridization of 35 bacterial artificial chromosomes (BACs) selected with markers mapping to eight linkage groups, plus two plasmids for 5S and 45S ribosomal DNA and one bacteriophage. Together with three previously mapped chromosomes (chromosomes 3, 4, and 7), 43 anchoring points between the genetic map and the cytogenetic map of the species are now available. Furthermore, a subset of four BAC clones was proposed to identify the 11 chromosome pairs of the standard cultivar BAT93. Three of these BACs labelled more than a single chromosome pair, indicating the presence of repetitive DNA in their inserts. A repetitive distribution pattern was observed for most of the BACs; for 38% of them, highly repetitive pericentromeric or subtelomeric signals were observed. These distribution patterns corresponded to pericentromeric and subtelomeric heterochromatin blocks observed with other staining methods. Altogether, the results indicate that around half of the common bean genome is heterochromatic and that genes and repetitive sequences are intermingled in the euchromatin and heterochromatin of the species.

Insights from the comparison of plant genome sequences

The next decade will see essentially completed sequences for multiple branches of virtually all angiosperm clades that include major crops and/or botanical models. These sequences will provide a powerful framework for relating genome-level events to aspects of morphological and physiological variation that have contributed to the colonization of much of the planet by angiosperms. Clarification of the fundamental angiosperm gene set, its arrangement, lineage-specific variations in gene repertoire and arrangement, and the fates of duplicated gene pairs will advance knowledge of functional and regulatory diversity and perhaps shed light on adaptation by lineages to whole-genome duplication, which is a distinguishing feature of angiosperm evolution. Better understanding of the relationships among angiosperm genomes promises to provide a firm foundation upon which to base translational genomics: the leveraging of hard-won structural and functional genomic information from crown botanical models to dissect novel and, in some cases, economically important features in many additional organisms.

Exploiting rice-sorghum synteny for targeted development of EST-SSRs to enrich the sorghum genetic linkage map

The sequencing and detailed comparative functional analysis of genomes of a number of select botanical models open new doors into comparative genomics among the angiosperms, with potential benefits for improvement of many orphan crops that feed large populations. In this study, a set of simple sequence repeat (SSR) markers was developed by mining the expressed sequence tag (EST) database of sorghum. Among the SSR-containing sequences, only those sharing considerable homology with rice genomic sequences across the lengths of the 12 rice chromosomes were selected. Thus, 600 SSR-containing sorghum EST sequences (50 homologous sequences on each of the 12 rice chromosomes) were selected, with the intention of providing coverage for corresponding homologous regions of the sorghum genome. Primer pairs were designed and polymorphism detection ability was assessed using parental pairs of two existing sorghum mapping populations. About 28% of these new markers detected polymorphism in this 4-entry panel. A subset of 55 polymorphic EST-derived SSR markers were mapped onto the existing skeleton map of a recombinant inbred population derived from cross N13 x E 36-1, which is segregating for Striga resistance and the stay-green component of terminal drought tolerance. These new EST-derived SSR markers mapped across all 10 sorghum linkage groups, mostly to regions expected based on prior knowledge of rice-sorghum synteny. The ESTs from which these markers were derived were then mapped in silico onto the aligned sorghum genome sequence, and 88% of the best hits corresponded to linkage-based positions. This study demonstrates the utility of comparative genomic information in targeted development of markers to fill gaps in linkage maps of related crop species for which sufficient genomic tools are not available.

Assignment of genetic linkage maps to diploid Solanum tuberosum pachytene chromosomes by BAC-FISH technology

A cytogenetic map has been developed for diploid potato (Solanum tuberosum), in which the arms of the 12 potato bivalents can be identified in pachytene complements using multicolor fluorescence in situ hybridization (FISH) with a set of 60 genetically anchored bacterial artificial chromosome (BAC) clones from the RHPOTKEY BAC library. This diagnostic set of selected BACs (five per chromosome) hybridizes to euchromatic regions and corresponds to well-defined loci in the ultradense genetic map, and with these probes a new detailed and reliable pachytene karyotype could be established. Chromosome size has been estimated both from microscopic length measurements and from 4′,6-diamidino-2-phenylindole fluorescence-based DNA content measurements. In both approaches, chromosome 1 is the largest (100–115 Mb) and chromosome 11 the smallest (49–53 Mb). Detailed measurements of mega-base-pair to micrometer ratios have been obtained for chromosome 5, with average values of 1.07 Mb/μm for euchromatin and 3.67 Mb/μm for heterochromatin. In addition, our FISH results helped to solve two discrepancies in the potato genetic map related to chromosomes 8 and 12. Finally, we discuss the significance of the potato cytogenetic map for extended FISH studies in potato and related Solanaceae, which will be especially beneficial for the potato genome-sequencing project.

Cytogenetic mapping of common bean chromosomes reveals a less compartmentalized small-genome plant species

Cytogenetic maps of common bean chromosomes 3, 4 and 7 were constructed by fluorescence in-situ hybridization (FISH) of BAC and a few other genomic clones. Although all clones were selected with genetically mapped markers, mostly with single-copy RFLPs, a large subset of BACs, from 13 different genomic regions, contained repetitive sequences, as concluded from the regional distribution patterns of multiple FISH signals on chromosomes: pericentromeric, subtelomeric and dispersed. Pericentromeric repeats were present in all 11 chromosome pairs with different intensities, whereas subtelomeric repeats were present in several chromosome ends, but with different signal intensities depending on the BAC, suggesting that the terminal heterochromatin fraction of this species may be composed of different repeats. The correlation of genetic and physical distances along the three studied chromosomes was obtained for 23 clones. This correlation suggests suppression of recombination around extended pericentromeric regions in a similar way to that previously reported for plant species with larger genomes. These results indicate that a relatively small plant genome may also possess a large proportion of repeats interspersed with single-copy sequences in regions other than the pericentromeric heterochromatin and, nevertheless, exhibit lower recombination around the pericentromeric fraction of the genome.

Microcollinearity in an ethylene receptor coding gene region of the Coffea canephora genome is extensively conserved with Vitis vinifera and other distant dicotyledonous sequenced genomes

BACKGROUND

Coffea canephora, also called Robusta, belongs to the Rubiaceae, the fourth largest angiosperm family. This diploid species (2x = 2n = 22) has a fairly small genome size of approximately 690 Mb and despite its extreme economic importance, particularly for developing countries, knowledge on the genome composition, structure and evolution remain very limited. Here, we report the 160 kb of the first C. canephora Bacterial Artificial Chromosome (BAC) clone ever sequenced and its fine analysis.

RESULTS

This clone contains the CcEIN4 gene, encoding an ethylene receptor, and twenty other predicted genes showing a high gene density of one gene per 7.8 kb. Most of them display perfect matches with C. canephora expressed sequence tags or show transcriptional activities through PCR amplifications on cDNA libraries. Twenty-three transposable elements, mainly Class II transposon derivatives, were identified at this locus. Most of these Class II elements are Miniature Inverted-repeat Transposable Elements (MITE) known to be closely associated with plant genes. This BAC composition gives a pattern similar to those found in gene rich regions of Solanum lycopersicum and Medicago truncatula genomes indicating that the CcEIN4 regions may belong to a gene rich region in the C. canephora genome. Comparative sequence analysis indicated an extensive conservation between C. canephora and most of the reference dicotyledonous genomes studied in this work, such as tomato (S. lycopersicum), grapevine (V. vinifera), barrel medic M. truncatula, black cottonwood (Populus trichocarpa) and Arabidopsis thaliana. The higher degree of microcollinearity was found between C. canephora and V. vinifera, which belong respectively to the Asterids and Rosids, two clades that diverged more than 114 million years ago.

CONCLUSION

This study provides a first glimpse of C. canephora genome composition and evolution. Our data revealed a remarkable conservation of the microcollinearity between C. canephora and V. vinifera and a high conservation with other distant dicotyledonous reference genomes. Altogether, these results provide valuable information to identify candidate genes in C. canephora genome and serve as a foundation to establish strategies for whole genome sequencing. Future large-scale sequence comparison between C. canephora and reference sequenced genomes will help in understanding the evolutionary history of dicotyledonous plants.

A consensus genetic map of sorghum that integrates multiple component maps and high-throughput Diversity Array Technology (DArT) markers

BACKGROUND

Sorghum genome mapping based on DNA markers began in the early 1990s and numerous genetic linkage maps of sorghum have been published in the last decade, based initially on RFLP markers with more recent maps including AFLPs and SSRs and very recently, Diversity Array Technology (DArT) markers. It is essential to integrate the rapidly growing body of genetic linkage data produced through DArT with the multiple genetic linkage maps for sorghum generated through other marker technologies. Here, we report on the colinearity of six independent sorghum component maps and on the integration of these component maps into a single reference resource that contains commonly utilized SSRs, AFLPs, and high-throughput DArT markers.

RESULTS

The six component maps were constructed using the MultiPoint software. The lengths of the resulting maps varied between 910 and 1528 cM. The order of the 498 markers that segregated in more than one population was highly consistent between the six individual mapping data sets. The framework consensus map was constructed using a "Neighbours" approach and contained 251 integrated bridge markers on the 10 sorghum chromosomes spanning 1355.4 cM with an average density of one marker every 5.4 cM, and were used for the projection of the remaining markers. In total, the sorghum consensus map consisted of a total of 1997 markers mapped to 2029 unique loci (1190 DArT loci and 839 other loci) spanning 1603.5 cM and with an average marker density of 1 marker/0.79 cM. In addition, 35 multicopy markers were identified. On average, each chromosome on the consensus map contained 203 markers of which 58.6% were DArT markers. Non-random patterns of DNA marker distribution were observed, with some clear marker-dense regions and some marker-rare regions.

CONCLUSION

The final consensus map has allowed us to map a larger number of markers than possible in any individual map, to obtain a more complete coverage of the sorghum genome and to fill a number of gaps on individual maps. In addition to overall general consistency of marker order across individual component maps, good agreement in overall distances between common marker pairs across the component maps used in this study was determined, using a difference ratio calculation. The obtained consensus map can be used as a reference resource for genetic studies in different genetic backgrounds, in addition to providing a framework for transferring genetic information between different marker technologies and for integrating DArT markers with other genomic resources. DArT markers represent an affordable, high throughput marker system with great utility in molecular breeding programs, especially in crops such as sorghum where SNP arrays are not publicly available.

Comparative sequence analysis of MONOCULM1-orthologous regions in 14 Oryza genomes

Comparative genomics is a powerful tool to decipher gene and genome evolution. Placing multiple genome comparisons in a phylogenetic context improves the sensitivity of evolutionary inferences. In the genus Oryza, this comparative approach can be used to investigate gene function, genome evolution, domestication, polyploidy, and ecological adaptation. A large genomic region surrounding the MONOCULM1 (MOC1) locus was chosen for study in 14 Oryza species, including 10 diploids and 4 allotetraploids. Sequencing and annotation of 18 bacterial artificial chromosome clones for these species revealed highly conserved gene colinearity and structure in the MOC1 region. Since the Oryza radiation about 14 Mya, differences in transposon amplification appear to be responsible for the different current sizes of the Oryza genomes. In the MOC1 region, transposons were only conserved between genomes of the same type (e.g., AA or BB). In addition to the conserved gene content, several apparent genes have been generated de novo or uniquely retained in the AA lineage. Two different 3-gene segments have been inserted into the MOC1 region of O. coarctata (KK) or O. sativa by unknown mechanism(s). Large and apparently noncoding sequences flanking the MOC1 gene were observed to be under strong purifying selection. The allotetraploids Oryza alta and Oryza minuta were found to be products of recent polyploidization, less than 1.6 and 0.4 Mya, respectively. In allotetraploids, pseudogenization of duplicated genes was common, caused by large deletions, small frame-shifting insertions/deletions, or nonsense mutations.

Chromatin structure and physical mapping of chromosome 6 of potato and comparative analyses with tomato

Potato (Solanum tuberosum) has the densest genetic linkage map and one of the earliest established cytogenetic maps among all plant species. However, there has been limited effort to integrate these maps. Here, we report fluorescence in situ hybridization (FISH) mapping of 30 genetic marker-anchored bacterial artificial chromosome (BAC) clones on the pachytene chromosome 6 of potato. The FISH mapping results allowed us to define the genetic positions of the centromere and the pericentromeric heterochromatin and to relate chromatin structure to the distribution of recombination along the chromosome. A drastic reduction of recombination was associated with the pericentromeric heterochromatin that accounts for approximately 28% of the physical length of the pachytene chromosome. The pachytene chromosomes 6 of potato and tomato (S. lycopersicum) share a similar morphology. However, distinct differences of heterochromatin distribution were observed between the two chromosomes. FISH mapping of several potato BACs on tomato pachytene chromosome 6 revealed an overall colinearity between the two chromosomes. A chromosome inversion was observed in the euchromatic region of the short arms. These results show that the potato and tomato genomes contain more chromosomal rearrangements than those reported previously on the basis of comparative genetic linkage mapping.

Analysis of repetitive DNA distribution patterns in the Tribolium castaneum genome

BACKGROUND

Insect genomes vary widely in size, a large fraction of which is often devoted to repetitive DNA. Re-association kinetics indicate that up to 42% of the genome of the red flour beetle, Tribolium castaneum, is repetitive. Analysis of the abundance and distribution of repetitive DNA in the recently sequenced genome of T. castaneum is important for understanding the structure and function of its genome.

RESULTS

Using TRF, TEpipe and RepeatScout we found that approximately 30% of the T. castaneum assembled genome is composed of repetitive DNA. Of this, 17% is found in tandem arrays and the remaining 83% is dispersed, including transposable elements, which in themselves constitute 5-6% of the genome. RepeatScout identified 31 highly repetitive DNA elements with repeat units longer than 100 bp, which constitute 7% of the genome; 65% of these highly repetitive elements and 74% of transposable elements accumulate in regions representing 40% of the assembled genome that is anchored to chromosomes. These regions tend to occur near one end of each chromosome, similar to previously described blocks of pericentric heterochromatin. They contain fewer genes with longer introns, and often correspond with regions of low recombination in the genetic map.

CONCLUSION

Our study found that transposable elements and other repetitive DNA accumulate in certain regions in the assembled T. castaneum genome. Several lines of evidence suggest these regions are derived from the large blocks of pericentric heterochromatin in T. castaneum chromosomes.