A gene-enriched BAC library for cloning large allele-specific fragments from maize: isolation of a 240-kb contig of the bronze region

Transcriptome Sequencing and Screening of Anthocyanin-Related Genes in the Leaves of Acer truncatum Bunge

Acer truncatum Bunge is generally used as an ornamental tree because of its autumn leaves, although the viewing period is short-approximately 7-15 days. Color improvement of ornamental trees has consistently been an important research topic because color partially determines the value of the commodity; however, a lack of genomic data have limited the progress of molecular breeding research in this area. The purposes of this study were to obtain a transcriptome database for A. truncatum, screen anthocyanin biosynthesis-related genes, and reveal the mechanisms underlying leaf color transformation to provide a basis for increasing the viewing period or breeding cultivars that display red leaves throughout the growing season via gene regulation. In this study, although the use of an Illumina HiSeq 2000 platform and systematic bioinformatics analysis using both young and mature leaves as experimental materials, 233,912,882 clean reads were generated and 121,287 unique transcripts were retrieved. We selected 16 color-related genes (from the transcriptome results) for qRT-PCR to validate the results, and the expression trends of the selected genes were largely consistent with the transcriptome analysis results, with a consistency of 0.875. According to the results of the transcriptome analysis, the validation, and previous studies, we obtained sequences of genes related to anthocyanins, including CHS, CHI, ANS, UFGT, UGT75c1, DFR, BZ1, F3H, F3'H, LAR, ANR, FLS, and those of several transcription factors, including MYB1, BHLH, and WD40. Verifying specific regulation by one or several of these genes in the control of leaf color requires further research. The acquisition of transcriptomic information, especially information concerning anthocyanin biosynthesis-related genes and their base sequences, can provide a theoretical basis for the study of the molecular mechanisms determining changes in leaf color in Acer and is of great importance to the breeding of new cultivars.

The unique genomic landscape surrounding the EPSPS gene in glyphosate resistant Amaranthus palmeri: a repetitive path to resistance

BACKGROUND

The expanding number and global distributions of herbicide resistant weedy species threaten food, fuel, fiber and bioproduct sustainability and agroecosystem longevity. Amongst the most competitive weeds, Amaranthus palmeri S. Wats has rapidly evolved resistance to glyphosate primarily through massive amplification and insertion of the 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) gene across the genome. Increased EPSPS gene copy numbers results in higher titers of the EPSPS enzyme, the target of glyphosate, and confers resistance to glyphosate treatment. To understand the genomic unit and mechanism of EPSPS gene copy number proliferation, we developed and used a bacterial artificial chromosome (BAC) library from a highly resistant biotype to sequence the local genomic landscape flanking the EPSPS gene.

RESULTS

By sequencing overlapping BACs, a 297 kb sequence was generated, hereafter referred to as the "EPSPS cassette." This region included several putative genes, dense clusters of tandem and inverted repeats, putative helitron and autonomous replication sequences, and regulatory elements. Whole genome shotgun sequencing (WGS) of two biotypes exhibiting high and no resistance to glyphosate was performed to compare genomic representation across the EPSPS cassette. Mapping of sequences for both biotypes to the reference EPSPS cassette revealed significant differences in upstream and downstream sequences relative to EPSPS with regard to both repetitive units and coding content between these biotypes. The differences in sequence may have resulted from a compounded-building mechanism such as repetitive transpositional events. The association of putative helitron sequences with the cassette suggests a possible amplification and distribution mechanism. Flow cytometry revealed that the EPSPS cassette added measurable genomic content.

CONCLUSIONS

The adoption of glyphosate resistant cropping systems in major crops such as corn, soybean, cotton and canola coupled with excessive use of glyphosate herbicide has led to evolved glyphosate resistance in several important weeds. In Amaranthus palmeri, the amplification of the EPSPS cassette, characterized by a complex array of repetitive elements and putative helitron sequences, suggests an adaptive structural genomic mechanism that drives amplification and distribution around the genome. The added genomic content not found in glyphosate sensitive plants may be driving evolution through genome expansion.

Dynamic evolution of bz orthologous regions in the Andropogoneae and other grasses

Genome structure exhibits remarkable plasticity within Zea mays. To examine how haplotype structure has evolved within the Andropogoneae tribe, we have analyzed the bz gene-rich region of maize (Zea mays), the Zea teosintes mays ssp. mexicana, luxurians and diploperennis, Tripsacum dactyloides, Coix lacryma-jobi and Sorghum propinquum. We sequenced and annotated BAC clones from these species and re-annotated the orthologous Sorghum bicolor region. Gene colinearity in the region is well conserved within the genus Zea. However, the orthologous regions of Coix and Sorghum exhibited several micro-rearrangements relative to Zea, including addition, truncation and deletion of genes. The stc1 gene, involved in the production of a terpenoid insect defense signal, is evolving particularly fast, and its progressive disappearance from some species is occurring by microhomology-mediated recombination. LTR retrotransposons are the main contributors to the dynamic evolution of the bz region. Common transposon insertion sites occur among haplotypes from different Zea mays sub-species, but not outside the species. As in Zea, different patterns of interspersion between genes and retrotransposons are observed in Sorghum. We estimate that the mean divergence times between maize and Tripsacum, Coix and Sorghum are 8.5, 12.1 and 12.4 million years ago, respectively, and that between Coix and Sorghum is 9.3 million years ago. A comparison of the bz orthologous regions of Zea, Sorghum and Coix with those of Brachypodium, Setaria and Oryza allows us to infer how the region has evolved by addition and deletion of genes in the approximately 50 million years since these genera diverged from a common progenitor.

Isolation of a 97-kb minimal essential MHC B locus from a new reverse-4D BAC library of the golden pheasant

The bacterial artificial chromosome (BAC) system is widely used in isolation of large genomic fragments of interest. Construction of a routine BAC library requires several months for picking clones and arraying BACs into superpools in order to employ 4D-PCR to screen positive BACs, which might be time-consuming and laborious. The major histocompatibility complex (MHC) is a cluster of genes involved in the vertebrate immune system, and the classical avian MHC-B locus is a minimal essential one, occupying a 100-kb genomic region. In this study, we constructed a more effective reverse-4D BAC library for the golden pheasant, which first creates sub-libraries and then only picks clones of positive sub-libraries, and identified several MHC clones within thirty days. The full sequencing of a 97-kb reverse-4D BAC demonstrated that the golden pheasant MHC-B locus contained 20 genes and showed good synteny with that of the chicken. The notable differences between these two species were the numbers of class II B loci and NK genes and the inversions of the TAPBP gene and the TAP1-TAP2 region. Furthermore, the inverse TAP2-TAP1 was unique in the golden pheasant in comparison with that of chicken, turkey, and quail. The newly defined genomic structure of the golden pheasant MHC will give an insight into the evolutionary history of the avian MHC.

Construction and characterization of a bacterial artificial chromosome library for the A-genome of cotton (G. arboreum L.)

A bacterial artificial chromosome (BAC) library for the A-genome of cotton has been constructed from the leaves of G. arboreum L cv. Jianglinzhongmian. It is used as elite A-genome germplasm resources in the present cotton breeding program and has been used to build a genetic reference map of cotton. The BAC library consists of 123,648 clones stored in 322 384-well plates. Statistical analysis of a set of 103 randomly selected BAC clones indicated that each clone has an average insert length of 100.2 kb per plasmid, with a range of 30 to 190 kb. Theoretically, this represents 7.2 haploid genome equivalents based on an A-genome size of 1697 Mb. The BAC library has been arranged in column pools and superpools allowing screening with various PCR-based markers. In the future, the A-genome cotton BAC library will serve as both a giant gene resource and a valuable tool for map-based gene isolation, physical mapping and comparative genome analysis.

Haplotype structure strongly affects recombination in a maize genetic interval polymorphic for Helitron and retrotransposon insertions

We have asked here how the remarkable variation in maize haplotype structure affects recombination. We compared recombination across a genetic interval of 9S in 2 highly dissimilar heterozygotes that shared 1 parent. The genetic interval in the common haplotype is approximately 100 kb long and contains 6 genes interspersed with gene-fragment-bearing Helitrons and retrotransposons that, together, comprise 70% of its length. In one heterozygote, most intergenic insertions are homozygous, although polymorphic, enabling us to determine whether any recombination junctions fall within them. In the other, most intergenic insertions are hemizygous and, thus, incapable of homologous recombination. Our analysis of the frequency and distribution of recombination in the interval revealed that: (i) Most junctions were circumscribed to the gene space, where they showed a highly nonuniform distribution. In both heterozygotes, more than half of the junctions fell in the stc1 gene, making it a clear recombination hotspot in the region. However, the genetic size of stc1 was 2-fold lower when flanked by a hemizygous 25-kb retrotransposon cluster. (ii) No junctions fell in the hypro1 gene in either heterozygote, making it a genic recombination coldspot. (iii) No recombination occurred within the gene fragments borne on Helitrons nor within retrotransposons, so neither insertion class contributes to the interval's genetic length. (iv) Unexpectedly, several junctions fell in an intergenic region not shared by all 3 haplotypes. (v) In general, the ability of a sequence to recombine correlated inversely with its methylation status. Our results show that haplotypic structural variability strongly affects the frequency and distribution of recombination events in maize.

Maize genome structure variation: interplay between retrotransposon polymorphisms and genic recombination

Although maize (Zea mays) retrotransposons are recombinationally inert, the highly polymorphic structure of maize haplotypes raises questions regarding the local effect of intergenic retrotransposons on recombination. To examine this effect, we compared recombination in the same genetic interval with and without a large retrotransposon cluster. We used three different bz1 locus haplotypes, McC, B73, and W22, in the same genetic background. We analyzed recombination between the bz1 and stc1 markers in heterozygotes that differ by the presence and absence of a 26-kb intergenic retrotransposon cluster. To facilitate the genetic screen, we used Ds and Ac markers that allowed us to identify recombinants by their seed pigmentation. We sequenced 239 recombination junctions and assigned them to a single nucleotide polymorphism-delimited interval in the region. The genetic distance between the markers was twofold smaller in the presence of the retrotransposon cluster. The reduction was seen in bz1 and stc1, but no recombination occurred in the highly polymorphic intergenic region of either heterozygote. Recombination within genes shuffled flanking retrotransposon clusters, creating new chimeric haplotypes and either contracting or expanding the physical distance between markers. Our findings imply that haplotype structure will profoundly affect the correlation between genetic and physical distance for the same interval in maize.

Remarkable variation in maize genome structure inferred from haplotype diversity at the bz locus

Maize is probably the most diverse of all crop species. Unexpectedly large differences among haplotypes were first revealed in a comparison of the bz genomic regions of two different inbred lines, McC and B73. Retrotransposon clusters, which comprise most of the repetitive DNA in maize, varied markedly in makeup, and location relative to the genes in the region and genic sequences, later shown to be carried by two helitron transposons, also differed between the inbreds. Thus, the allelic bz regions of these Corn Belt inbreds shared only a minority of the total sequence. To investigate further the variation caused by retrotransposons, helitrons, and other insertions, we have analyzed the organization of the bz genomic region in five additional cultivars selected because of their geographic and genetic diversity: the inbreds A188, CML258, and I137TN, and the land races Coroico and NalTel. This vertical comparison has revealed the existence of several new helitrons, new retrotransposons, members of every superfamily of DNA transposons, numerous miniature elements, and novel insertions flanked at either end by TA repeats, which we call TAFTs (TA-flanked transposons). The extent of variation in the region is remarkable. In pairwise comparisons of eight bz haplotypes, the percentage of shared sequences ranges from 25% to 84%. Chimeric haplotypes were identified that combine retrotransposon clusters found in different haplotypes. We propose that recombination in the common gene space greatly amplifies the variability produced by the retrotransposition explosion in the maize ancestry, creating the heterogeneity in genome organization found in modern maize.

Creation of a BAC resource to study the structure and evolution of the banana (Musa balbisiana) genome

The first bacterial artificial chromosome (BAC) library of the banana species Musa balbisiana 'Pisang Klutuk Wulung' (PKW BAC library) was constructed and characterized. One improved and one novel protocol for nuclei isolation were employed to overcome problems caused by high levels of polyphenols and polysaccharides present in leaf tissues. The use of flow cytometry to purify cell nuclei eliminated contamination with secondary metabolites and plastid DNA. Furthermore, the usefulness of the inducible pCC1BAC vector to obtain a higher amount of BAC DNA was demonstrated. The PKW BAC library represents nine haploid genome equivalents of M. balbisiana and its mean insert size is 135 kb. It consists of two sublibraries, of which the first one (SN sublibrary with 24,960 clones) was prepared according to an improved standard nuclei isolation protocol, whereas the second (FN sublibrary with 11,904 clones) was obtained from flow-sorted nuclei. Screening with 12 RFLP probes, which were genetically anchored to 8 genetic linkage groups of the banana species Musa acuminata, revealed an average of 11 BAC clones per probe, thus confirming the genome coverage estimated based on the insert size, as well as a high level of conservation between the two species of Musa. Localization of selected BAC clones to mitotic chromosomes using FISH indicated that the BAC library represented a useful resource for cytogenetic mapping. As the first step in map-based cloning of a genetic factor that is involved in the activation of integrated pararetroviral sequences of Banana streak virus (BSV), the BSV expressed locus (BEL) was physically delimited. The PKW BAC library represents a publicly available tool, and is currently used to reveal the integration and activation mechanisms of BSV sequences and to study banana genome structure and evolution.

Direct targeting and rapid isolation of BAC clones spanning a defined chromosome region

To isolate genes of interest in plants, it is essential to construct bacterial artificial chromosome (BAC) libraries from specific genotypes. Construction and organisation of BAC libraries is laborious and costly, especially from organisms with large and complex genomes. In the present study, we developed the pooled BAC library strategy that allows rapid and low cost generation and screening of genomic libraries from any genotype of interest. The BAC library is constructed, directly organised into a few pools and screened for BAC clones of interest using PCR and hybridisation steps, without requiring organization into individual clones. As a proof of concept, a pooled BAC library of approximately 177,000 recombinant clones has been constructed from the barley cultivar Cebada Capa that carries the Rph7 leaf rust resistance gene. The library has an average insert size of 140 kb, a coverage of six barley genome equivalents and is organised in 138 pools of about 1,300 clones each. We rapidly established a single contig of six BAC clones spanning 230 kb at the Rph7 locus on chromosome 3HS. The described low-cost cloning strategy is fast and will greatly facilitate direct targeting of genes and large-scale intra- and inter-species comparative genome analysis.

A new single-locus cytogenetic mapping system for maize (Zea mays L.): overcoming FISH detection limits with marker-selected sorghum (S. propinquum L.) BAC clones

The development of a cytogenetic map for maize (Zea mays L.) is shown to be feasible by means of a combination of resources from sorghum and oat that overcome limitations of single-copy gene detection. A maize chromosome-addition line of oat, OMAd9.2, provided clear images of optically isolated pachytene chromosomes through a chromosome spread and painting technique. A direct labeled oligonucleotide fluorescence in situ hybridization (FISH) probe MCCY specifically stained the centromere. The arm ratio (long/short) for maize chromosome 9 in the addition line was 1.7, comparable to the range of 1.6-2.1 previously reported for maize chromosome 9. A sorghum (Sorghum propinquum L.) BAC library was screened by hybridization with each of three maize core-bin-marker (CBM) probes: umc109 (CBM9.01), umc192/bz1 (CBM9.02), and csu54b (CBM9.08). A single BAC clone for each marker was chosen; designated sCBM9.1, sCBM9.2, or sCBM9.8; and used as a FISH probe on pachytene spreads from OMAd9.2. In each case, discrete FISH signals were observed, and their cytogenetic positions were determined to be 9S.79 (at position 79% of the length of chromosome 9 short arm) for sCBM9.1, 9S.65 for sCBM9.2, and approximately 9L.95 for sCBM9.8. These map positions were co-linear with linkage-map positions for these and other loci common to the linkage and cytogenetic maps. This work represents a major breakthrough for cytogenetic mapping of the maize genome, and also provides a general strategy that can be applied to cytogenetic mapping of other plant species with relatively large and complex genomes.

Intraspecific violation of genetic colinearity and its implications in maize

Although allelic sequences can vary extensively, it is generally assumed that each gene in one individual will have an allelic counterpart in another individual of the same species. We report here that this assumption does not hold true in maize. We have sequenced over 100 kb from the bz genomic region of two different maize lines and have found dramatic differences between them. First, the retrotransposon clusters, which comprise most of the repetitive DNA in maize, differ markedly in make-up and location relative to the genes in the bz region. Second, and more importantly, the genes themselves differ between the two lines, demonstrating that genetic microcolinearity can be violated within the same species. Our finding has bearing on the underlying genetic basis of hybrid vigor in maize, and possibly other organisms, and on the measurement of genetic distances.

Recombination rates between adjacent genic and retrotransposon regions in maize vary by 2 orders of magnitude

Genetic map length and gene number in eukaryotes vary considerably less than genome size, giving rise to the hypothesis that recombination is restricted to genes. The complex genome of maize contains a large fraction of repetitive DNA, composed principally of retrotransposons arranged in clusters. Here, we assess directly the contribution of retrotransposon clusters and genes to genetic length. We first measured recombination across adjacent homozygous genetic intervals on either side of the bronze (bz) locus. We then isolated and characterized two bacterial artificial chromosome clones containing those intervals. Recombination was almost 2 orders of magnitude higher in the distal side, which is gene-dense and lacks retrotransposons, than in the proximal side, which is gene-poor and contains a large cluster of methylated retrotransposons. We conclude that the repetitive retrotransposon DNA in maize, which constitutes the bulk of the genome, most likely contributes little if any to genetic length.

The highly recombinogenic bz locus lies in an unusually gene-rich region of the maize genome

The bronze (bz) locus exhibits the highest rate of recombination of any gene in higher plants. To investigate the possible basis of this high rate of recombination, we have analyzed the physical organization of the region around the bz locus. Two adjacent bacterial artificial chromosome clones, comprising a 240-kb contig centered around the Bz-McC allele, were isolated, and 60 kb of contiguous DNA spanning the two bacterial artificial chromosome clones was sequenced. We find that the bz locus lies in an unusually gene-rich region of the maize genome. Ten genes, at least eight of which are shown to be transcribed, are contained in a 32-kb stretch of DNA that is uninterrupted by retrotransposons. We have isolated nearly full length cDNAs corresponding to the five proximal genes in the cluster. The average intertranscript distance between them is just 1 kb, revealing a surprisingly compact packaging of adjacent genes in this part of the genome. At least 11 small insertions, including several previously described miniature inverted repeat transposable elements, were detected in the introns and 3' untranslated regions of genes and between genes. The gene-rich region is flanked at the proximal and distal ends by retrotransposon blocks. Thus, the maize genome appears to have scattered regions of high gene density similar to those found in other plants. The unusually high rate of intragenic recombination seen in bz may be related to the very high gene density of the region.

A maize sesquiterpene cyclase gene induced by insect herbivory and volicitin: characterization of wild-type and mutant alleles

Plants can defend themselves from herbivorous insects by emitting volatile chemical signals that attract natural enemies of the herbivore. For example, maize seedlings attacked by beet armyworm larvae (Spodoptera exigua) produce a mixture of terpenoid and indole volatiles that serve to attract parasitic wasps. A key step in terpenoid biosynthesis is the conversion of acyclic prenyl diphosphates to terpenoid compounds by specific terpenoid synthases (cyclases). We have cloned a maize sesquiterpene cyclase gene, stc1, by transposon tagging and have identified two deletion mutations of the gene. The stc1 gene is located on chromosome 9S and does not seem to have a closely related ortholog in the maize genome. It is induced 15- to 30-fold in maize leaves by beet armyworm larvae feeding or by application of purified volicitin, the insect-derived elicitor, at a mechanically wounded site. stc1 induction is systemic, because undamaged leaves of the same plant show a similar increase in stc1 transcription. Analysis of volatiles from volicitin-treated seedlings revealed that a major naphthalene-based sesquiterpene was present in wild-type seedlings but absent in the Ac-insertion and x-ray-deletion mutants. Therefore, we have identified a maize gene that responds to caterpillar herbivory by producing a chemical defense signal that most likely serves to attract natural enemies of the herbivore.