Frequencies of Ty1- copia and Ty3- gypsy retroelements within the Triticeae EST databases

Chimeras Linked to Tandem Repeats and Transposable Elements in Tetraploid Hybrid Fish

The formation of the allotetraploid hybrid lineage (4nAT) encompasses both distant hybridization and polyploidization processes. The allotetraploid offspring have two sets of sub-genomes inherited from both parental species, and therefore, it is important to explore its genetic structure. Herein, we construct a bacterial artificial chromosome library of allotetraploids, and then sequence and analyze the full-length sequences of 19 bacterial artificial chromosomes. Sixty-eight DNA chimeras are identified, which are divided into four models according to the distribution of the genomic DNA derived from the parents. Among the 68 genetic chimeras, 44 (64.71%) are linked to tandem repeats (TRs) and 23 (33.82%) are linked to transposable elements (TEs). The chimeras linked to TRs are related to slipped-strand mispairing and double-strand break repair while the chimeras linked to TEs benefit from the intervention of recombinases. In addition, TRs and TEs can also result in insertions/deletions of DNA segments. We conclude that DNA chimeras accompanied by TRs and TEs coordinate a balance between the sub-genomes derived from the parents. It is the first report on the relationship between formation of the DNA chimeras and TRs and TEs in the polyploid animals.

Apomixis frequency under stress conditions in weeping lovegrass (Eragrostis curvula)

To overcome environmental stress, plants develop physiological responses that are triggered by genetic or epigenetic changes, some of which involve DNA methylation. It has been proposed that apomixis, the formation of asexual seeds without meiosis, occurs through the temporal or spatial deregulation of the sexual process mediated by genetic and epigenetic factors influenced by the environment. Here, we explored whether there was a link between the occurrence of apomixis and various factors that generate stress, including drought stress, in vitro culture, and intraspecific hybridization. For this purpose, we monitored the embryo sacs of different weeping lovegrass (Eragrostis curvula [Schrad.] Nees) genotypes after the plants were subjected to these stress conditions. Progeny tests based on molecular markers and genome methylation status were analyzed following the stress treatment. When grown in the greenhouse, the cultivar Tanganyika INTA generated less than 2% of its progeny by sexual reproduction. Plants of this cultivar subjected to different stresses showed an increase of sexual embryo sacs, demonstrating an increased expression of sexuality compared to control plants. Plants of the cv. Tanganyika USDA did not demonstrate the ability to generate sexual embryo sacs under any conditions and is therefore classified as a fully apomictic cultivar. We found that this change in the prevalence of sexuality was correlated with genetic and epigenetic changes analyzed by MSAP and AFLPs profiles. Our results demonstrate that different stress conditions can alter the expression of sexual reproduction in facultative tetraploid apomictic cultivars and when the stress stops the reproductive mode shift back to the apomixis original level. These data together with previous observations allow us to generate a hypothetical model of the regulation of apomixis in weeping lovegrass in which the genetic/s region/s that condition apomixis, is/are affected by ploidy, and is/are subjected to epigenetic control.

RNA-seq studies using wheat PHYTOCHROME B and PHYTOCHROME C mutants reveal shared and specific functions in the regulation of flowering and shade-avoidance pathways

BACKGROUND

In cereal crops such as wheat, an optimal timing of developmental transitions is required to maximize grain yield. Many of these developmental changes are precisely regulated by changes in the duration, intensity or quality of light. Phytochromes are dimeric photoreceptors that absorb light maximally in the red and far-red wavelengths and induce large-scale transcriptional changes in response to variation in light quality. In wheat, PHYC is required for early flowering under long days. However, it is currently unknown whether this function requires the presence of PHYB. In this study, we characterized the role of PHYB in wheat development and used RNA-seq to analyze and compare the transcriptomes of phyB-null and phyC-null TILLING mutants.

RESULTS

Under long-day photoperiods, phyB-null plants exhibit a severe delay in flowering comparable to the delay observed in phyC-null plants. These results demonstrate that both genes are required for the induction of wheat flowering under long days. Using replicated RNA-seq studies we identified 82 genes that are significantly up or down regulated in both the phyB-null and phyC-null mutant relative to their respective wild-type controls. Among these genes are several well-characterized positive regulators of flowering, including PPD1, FT1 and VRN1. Eight-fold more genes were differentially regulated only in the phyB-null mutant (2202) than only in the phyC-null mutant (261). The PHYB-regulated genes were enriched in components of the auxin, gibberellin and brassinosteroid biosynthesis and signaling pathways, and in transcription factors with putative roles in regulating vegetative development and shade-avoidance responses. Several genes involved in abiotic stress tolerance pathways were also found to be regulated by PHYB.

CONCLUSIONS

PHYB and PHYC are both required for the photoperiodic induction of wheat flowering, whereas PHYB alone regulates a large number of genes involved in hormone biosynthesis and signaling, shade-avoidance response, and abiotic stress tolerance. Our analysis provides a comprehensive overview of the PHYB- and PHYC-mediated transcriptional changes during light signaling, and an initial step towards the dissection of this regulatory gene network in wheat. This further dissection will be required to explore the individual phytochrome-mediated developmental responses and to evaluate their potential to improve wheat adaptation to changing environments.

Transcriptionally active LTR retrotransposons in Eucalyptus genus are differentially expressed and insertionally polymorphic

BACKGROUND

In Eucalyptus genus, studies on genome composition and transposable elements (TEs) are particularly scarce. Nearly half of the recently released Eucalyptus grandis genome is composed by retrotransposons and this data provides an important opportunity to understand TE dynamics in Eucalyptus genome and transcriptome.

RESULTS

We characterized nine families of transcriptionally active LTR retrotransposons from Copia and Gypsy superfamilies in Eucalyptus grandis genome and we depicted genomic distribution and copy number in two Eucalyptus species. We also evaluated genomic polymorphism and transcriptional profile in three organs of five Eucalyptus species. We observed contrasting genomic and transcriptional behavior in the same family among different species. RLC_egMax_1 was the most prevalent family and RLC_egAngela_1 was the family with the lowest copy number. Most families of both superfamilies have their insertions occurring <3 million years, except one Copia family, RLC_egBianca_1. Protein theoretical models suggest different properties between Copia and Gypsy domains. IRAP and REMAP markers suggested genomic polymorphisms among Eucalyptus species. Using EST analysis and qRT-PCRs, we observed transcriptional activity in several tissues and in all evaluated species. In some families, osmotic stress increases transcript values.

CONCLUSION

Our strategy was successful in isolating transcriptionally active retrotransposons in Eucalyptus, and each family has a particular genomic and transcriptional pattern. Overall, our results show that retrotransposon activity have differentially affected genome and transcriptome among Eucalyptus species.

Analysis of plant LTR-retrotransposons at the fine-scale family level reveals individual molecular patterns

BACKGROUND

Sugarcane is an important crop worldwide for sugar production and increasingly, as a renewable energy source. Modern cultivars have polyploid, large complex genomes, with highly unequal contributions from ancestral genomes. Long Terminal Repeat retrotransposons (LTR-RTs) are the single largest components of most plant genomes and can substantially impact the genome in many ways. It is therefore crucial to understand their contribution to the genome and transcriptome, however a detailed study of LTR-RTs in sugarcane has not been previously carried out.

RESULTS

Sixty complete LTR-RT elements were classified into 35 families within four Copia and three Gypsy lineages. Structurally, within lineages elements were similar, between lineages there were large size differences. FISH analysis resulted in the expected pattern of Gypsy/heterochromatin, Copia/euchromatin, but in two lineages there was localized clustering on some chromosomes. Analysis of related ESTs and RT-PCR showed transcriptional variation between tissues and families. Four distinct patterns were observed in sRNA mapping, the most unusual of which was that of Ale1, with very large numbers of 24nt sRNAs in the coding region. The results presented support the conclusion that distinct small RNA-regulated pathways in sugarcane target the lineages of LTR-RT elements.

CONCLUSIONS

Individual LTR-RT sugarcane families have distinct structures, and transcriptional and regulatory signatures. Our results indicate that in sugarcane individual LTR-RT families have distinct behaviors and can potentially impact the genome in diverse ways. For instance, these transposable elements may affect nearby genes by generating a diverse set of small RNA's that trigger gene silencing mechanisms. There is also some evidence that ancestral genomes contribute significantly different element numbers from particular LTR-RT lineages to the modern sugarcane cultivar genome.

Effect of the down-regulation of the high Grain Protein Content (GPC) genes on the wheat transcriptome during monocarpic senescence

BACKGROUND

Increasing the nutrient concentration of wheat grains is important to ameliorate nutritional deficiencies in many parts of the world. Proteins and nutrients in the wheat grain are largely derived from the remobilization of degraded leaf molecules during monocarpic senescence. The down-regulation of the NAC transcription factor Grain Protein Content (GPC) in transgenic wheat plants delays senescence (>3 weeks) and reduces the concentration of protein, Zn and Fe in the grain (>30%), linking senescence and nutrient remobilization.Based on the early and rapid up-regulation of GPC in wheat flag leaves after anthesis, we hypothesized that this transcription factor is an early regulator of monocarpic senescence. To test this hypothesis, we used high-throughput mRNA-seq technologies to characterize the effect of the GPC down-regulation on the wheat flag-leaf transcriptome 12 days after anthesis. At this early stage of senescence GPC transcript levels are significantly lower in transgenic GPC-RNAi plants than in the wild type, but there are still no visible phenotypic differences between genotypes.

RESULTS

We generated 1.4 million 454 reads from early senescing flag leaves (average ~350 nt) and assembled 1.2 million into 30,497 contigs that were used as a reference to map 145 million Illumina reads from three wild type and four GPC-RNAi plants. Following normalization and statistical testing, we identified a set of 691 genes differentially regulated by GPC (431 ≥ 2-fold change). Transcript level ratios between transgenic and wild type plants showed a high correlation (R = 0.83) between qRT-PCR and Illumina results, providing independent validation of the mRNA-seq approach. A set of differentially expressed genes were analyzed across an early senescence time-course.

CONCLUSIONS

Monocarpic senescence is an active process characterized by large-scale changes in gene expression which begins considerably before the appearance of visual symptoms of senescence. The mRNA-seq approach used here was able to detect small differences in transcript levels during the early stages of senescence. This resulted in an extensive list of GPC-regulated genes, which includes transporters, hormone regulated genes, and transcription factors. These GPC-regulated genes, particularly those up-regulated during senescence, provide valuable entry points to dissect the early stages of monocarpic senescence and nutrient remobilization in wheat.

Use of methylation filtration and C(0)t fractionation for analysis of genome composition and comparative genomics in bread wheat

We investigated the compositional and structural differences in sequences derived from different fractions of wheat genomic DNA obtained using methylation filtration and C(0)t fractionation. Comparative analysis of these sequences revealed large compositional and structural variations in terms of GC content, different structural elements including repeat sequences (e.g., transposable elements and simple sequence repeats), protein coding genes, and non-coding RNA genes. A correlation between methylation status [determined on the basis of selective inclusion/exclusion in methylation-filtered (MF) library] of different repeat elements and expression level was observed. The expression levels were determined by comparing MF sequences with expressed sequence tags (ESTs) available in the public domain. Only a limited overlap among MF, high C(0)t (HC), and ESTs was observed, suggesting that these sequences may largely either represent the low-copy non-transcribed sequences or include genes with low expression levels. Thus, these results indicated a need to study MF and HC sequences along with ESTs to fully appreciate complexity of wheat gene space.

Small RNAs, DNA methylation and transposable elements in wheat

Background

More than 80% of the wheat genome is composed of transposable elements (TEs). Since active TEs can move to different locations and potentially impose a significant mutational load, their expression is suppressed in the genome via small non-coding RNAs (sRNAs). sRNAs guide silencing of TEs at the transcriptional (mainly 24-nt sRNAs) and post-transcriptional (mainly 21-nt sRNAs) levels. In this study, we report the distribution of these two types of sRNAs among the different classes of wheat TEs, the regions targeted within the TEs, and their impact on the methylation patterns of the targeted regions.

Results

We constructed an sRNA library from hexaploid wheat and developed a database that included our library and three other publicly available sRNA libraries from wheat. For five completely-sequenced wheat BAC contigs, most perfectly matching sRNAs represented TE sequences, suggesting that a large fraction of the wheat sRNAs originated from TEs. An analysis of all wheat TEs present in the Triticeae Repeat Sequence database showed that sRNA abundance was correlated with the estimated number of TEs within each class. Most of the sRNAs perfectly matching miniature inverted repeat transposable elements (MITEs) belonged to the 21-nt class and were mainly targeted to the terminal inverted repeats (TIRs). In contrast, most of the sRNAs matching class I and class II TEs belonged to the 24-nt class and were mainly targeted to the long terminal repeats (LTRs) in the class I TEs and to the terminal repeats in CACTA transposons. An analysis of the mutation frequency in potentially methylated sites revealed a three-fold increase in TE mutation frequency relative to intron and untranslated genic regions. This increase is consistent with wheat TEs being preferentially methylated, likely by sRNA targeting.

Conclusions

Our study examines the wheat epigenome in relation to known TEs. sRNA-directed transcriptional and post-transcriptional silencing plays important roles in the short-term suppression of TEs in the wheat genome, whereas DNA methylation and increased mutation rates may provide a long-term mechanism to inactivate TEs.

Control of flowering time and spike development in cereals: the earliness per se Eps-1 region in wheat, rice, and Brachypodium

The earliness per se gene Eps-A^m1 from diploid wheat Triticum monococcum affects heading time, spike development, and spikelet number. In this study, the Eps1 orthologous regions from rice, Aegilops tauschii, and Brachypodium distachyon were compared as part of current efforts to clone this gene. A single Brachypodium BAC clone spanned the Eps-A^m1 region, but a gap was detected in the A. tauschii physical map. Sequencing of the Brachypodium and A. tauschii BAC clones revealed three genes shared by the three species, which showed higher identity between wheat and Brachypodium than between them and rice. However, most of the structural changes were detected in the wheat lineage. These included an inversion encompassing the wg241-VatpC region and the presence of six unique genes. In contrast, only one unique gene (and one pseudogene) was found in Brachypodium and none in rice. Three genes were present in both Brachypodium and wheat but were absent in rice. Two of these genes, Mot1 and FtsH4, were completely linked to the earliness per se phenotype in the T. monococcum high-density genetic map and are candidates for Eps-A^m1. Both genes were expressed in apices and developing spikes, as expected for Eps-A^m1 candidates. The predicted MOT1 protein showed amino acid differences between the parental T. monococcum lines, but its effect is difficult to predict. Future steps to clone the Eps-A^m1 gene include the generation of mot1 and ftsh4 mutants and the completion of the T. monococcum physical map to test for the presence of additional candidate genes.

Isolation and characterization of genomic retrotransposon sequences from octoploid strawberry (Fragaria x ananassa Duch.)

Strawberry (Fragaria spp.) is a kind of herbaceous perennial plant that propagates vegetatively. The conserved domains of reverse transcriptase (RT) genes of Ty1-copia and Ty3-gypsy groups of LTR retrotransposons were amplified from the cultivated strawberry (Fragaria x ananassa Duch.). Sequence analysis of clones demonstrated that 5 of 19 Ty1-copia group unique sequences and 2 of 10 Ty3-gypsy unique sequences in F. x ananassa genome possessed either stop codon or frameshift. Ty1-copia group sequences are highly heterogeneous (divergence ranged from 1 to 69.8%), but the Ty3-gypsy group sequences are less (divergence ranged from 1 to 10%). Southern dot blot hybridization result suggested that both of the LTR retrotransposons are present in the genome of cultivated strawberry with high copy number (Ty1-copia group 2,875 Ty3-gypsy group 348). RT-PCR amplification from total RNA, which was extracted from leaves of micropropagated strawberry plants, did not yield either of the RT fragments. This is the first report on the presence of RT sequences of Ty1-copia and Ty3-gypsy group retrotransposons in F. x ananassa genome.

Reme1, a Copia retrotransposon in melon, is transcriptionally induced by UV light

For the first time, numerous sequences of Copia and Gypsy retrotransposons from the Cucumis melo genome have been obtained and analyzed. Phylogenetic analyses of sequences of both types of long terminal repeat (LTR) retrotransposons were carried out. The melon genome contains approximately 20,000 Gypsy and 6,800 Copia elements, comprising about 26% of its total size. Starting from a retrotransposon fragment, we have cloned and characterized an entire melon retrotransposon, named Reme1, which is 5,149 bp long. Reme1 belongs to the Superfamily Copia retrotransposons by its protein domain order and sequence similarity to other Copia elements of dicotyledons. The haploid genome of melon (var. "Piel de Sapo") contains about 120 copies of Reme1. Several copies of Reme1 are transcriptionally active, although at low levels, in melon leaves as analyzed by reverse-transcription PCR (RT-PCR) and sequencing. However, the transcript pool is considerably increased when melon leaves are treated with UV light, as has been seen for various retroelements in many organisms. The cDNAs of Reme1 transcripts showed less diversity than do Reme1 genomic sequences, suggesting that a subfamily of these elements is differentially responsive to UV.

Studies on genome relationship and species-specific PCR marker for Dasypyrum breviaristatum in Triticeae

Dasypyrum breviaristatum and nine related species in Triticeae were analyzed using the random amplified polymorphic DNA (RAPD) technique, in order to understand the genetic relationship and to develop species specific markers. The genome relationship dendrogram shows that D. breviaristatum and D. villosum could not be grouped together, indicating that D. breviaristatum was unlikely to be directly derived from D. villosum, while D. breviaristatum was closest to Thinopyrum intermedium, which implied that they might have similar breeding behaviors when introducing their chromatins into wheat. A D. breviaristatum genome specific RAPD product of 1182bp, was cloned and designated as pDb12H. Sequence analysis revealed that pDb12H was strongly homologuos to a long terminal repeat (LTR) Sabrina retrotransposon newly reported in Hordeum. The pDb12H was converted into a PCR based marker, which allows effectively monitoring the D. breviaristatum chromatin introgression into wheat. Fluorescence in situ hybridization (FISH) suggested that pDb12H was specifically hybridized throughout all D. breviaristatum chromosomes arms except for the terminal and centromeric regions, which can be used to characterize wheat -D. breviaristatum chromosome translocation. The genomes repetitive element will also be useful to study gene interactions between the wheat and alien genomes after the polyploidization.

The centromeric retrotransposons of rice are transcribed and differentially processed by RNA interference

Retrotransposons consist of significant portions of many complex eukaryotic genomes and are often enriched in heterochromatin. The centromeric retrotransposon (CR) family in grass species is colonized in the centromeres and highly conserved among species that have been diverged for >50 MY. These unique characteristics have inspired scientists to speculate about the roles of CR elements in organization and function of centromeric chromatin. Here we report that the CRR (CR of rice) elements in rice are highly enriched in chromatin associated with H3K9me2, a hallmark for heterochromatin. CRR elements were transcribed in root, leaf, and panicle tissues, suggesting a constitutive transcription of this retrotransposon family. However, the overall transcription level was low and the CRR transcripts appeared to be derived from relatively few loci. The majority of the CRR transcripts had chimerical structures and contained only partial CRR sequences. We detected small RNAs (smRNAs) cognate to nonautonomous CRR1 (noaCRR1) and CRR1, but not CRR2 elements. This result was also confirmed by in silico analysis of rice smRNA sequences. These results suggest that different CRR subfamilies may play different roles in the RNAi-mediated pathway for formation and maintenance of centromeric heterochromatin.

Analysis of the wheat endosperm transcriptome

Among the cereals, wheat is the most widely grown geographically and is part of the staple diet in much of the world. Understanding how the cereal endosperm develops and functions will help generate better tools to manipulate grain qualities important to end-users. We used a genomics approach to identify and characterize genes that are expressed in the wheat endosperm. We analyzed the 17,949 publicly available wheat endosperm EST sequences to identify genes involved in the biological processes that occur within this tissue. Clustering and assembly of the ESTs resulted in the identification of 6,187 tentative unique genes, 2,358 of which formed contigs and 3,829 remained as singletons. A BLAST similarity search against the NCBI non-redundant sequence database revealed abundant messages for storage proteins, putative defense proteins, and proteins involved in starch and sucrose metabolism. The level of abundance of the putatively identified genes reflects the physiology of the developing endosperm. Half of the identified genes have unknown functions. Approximately 61% of the endosperm ESTs has been tentatively mapped in the hexaploid wheat genome. Using microarrays for global RNA profiling, we identified endosperm genes that are specifically up regulated in the developing grain.

Identification of stress-responsive genes in an indica rice (Oryza sativa L.) using ESTs generated from drought-stressed seedlings

The impacts of drought on plant growth and development limit cereal crop production worldwide. Rice (Oryza sativa) productivity and production is severely affected due to recurrent droughts in almost all agroecological zones. With the advent of molecular and genomic technologies, emphasis is now placed on understanding the mechanisms of genetic control of the drought-stress response. In order to identify genes associated with water-stress response in rice, ESTs generated from a normalized cDNA library, constructed from drought-stressed leaf tissue of an indica cultivar, Nagina 22 were used. Analysis of 7794 cDNA sequences led to the identification of 5815 rice ESTs. Of these, 334 exhibited no significant sequence homology with any rice ESTs or full-length cDNAs in public databases, indicating that these transcripts are enriched during drought stress. Analysis of these 5815 ESTs led to the identification of 1677 unique sequences. To characterize this drought transcriptome further and to identify candidate genes associated with the drought-stress response, the rice data were compared with those for abiotic stress-induced sequences obtained from expression profiling studies in Arabidopsis, barley, maize, and rice. This comparative analysis identified 589 putative stress-responsive genes (SRGs) that are shared by these diverse plant species. Further, the identified leaf SRGs were compared to expression profiles for a drought-stressed rice panicle library to identify common sequences. Significantly, 125 genes were found to be expressed under drought stress in both tissues. The functional classification of these 125 genes showed that a majority of them are associated with cellular metabolism, signal transduction, and transcriptional regulation.

A movable feast: diverse retrotransposons and their contribution to barley genome dynamics

Cellular genes comprise at most 5% of the barley genome; the rest is occupied primarily by retrotransposons. Retrotransposons move intracellularly by a replicative mechanism similar to that of retroviruses. We describe the major classes of retrotransposons in barley, including the two nonautonomous groups that were recently identified, and detail the evidence supporting our current understanding of their life cycle. Data from analyses of long contiguous segments of the barley genome, as well as surveys of the prevalence of full-length retrotransposons and their solo LTR derivatives in the genus Hordeum, indicate that integration and recombinational loss of retrotransposons are major factors shaping the genome. The sequence conservation and integrative capacity of barley retrotransposons have made them excellent sources for development of molecular marker systems.

Phenotypic impacts of repetitive DNA in flowering plants

The discovery that nuclear DNA content varies widely among species, and even within species, was unexpected because it was thought that the number of genes required for an organism should be common across taxa. We now know that the bulk of nuclear DNA content variation is caused by repetitive DNA sequences characterized according to the nature of repeat (tandem vs dispersed) or chromosomal location/mechanism of replication (pericentromeric, telomeric or subtelomeric, microsatellites, minisatellites, satellites, transposable elements, retroelements). Variation in repetitive DNA, manifested as variation in nuclear DNA content, has been shown to have broad ecological and life-history consequences. For example, large genome size appears to limit fitness in extreme environmental conditions. Within species, variation in DNA content has been coupled to growth and development, such as maturation time in crop species. In Silene latifolia, DNA content is negatively correlated with flower size, a character that, in turn, has well documented ecological significance. These intraspecific studies suggest a connection between repetitive DNA and quantitative genetic determination of continuous characters. Novel insights into mechanisms by which repetitive DNA influences phenotype will lead to models of evolutionary change that extend well beyond the conventional view of evolution by allelic substitution.

Tempos of gene locus deletions and duplications and their relationship to recombination rate during diploid and polyploid evolution in the Aegilops-Triticum alliance

The origin of tetraploid wheat and the divergence of diploid ancestors of wheat A and D genomes were estimated to have occurred 0.36 and 2.7 million years ago, respectively. These estimates and the evolutionary history of 3159 gene loci were used to estimate the rates with which gene loci have been deleted and duplicated during the evolution of wheat diploid ancestors and during the evolution of polyploid wheat. During diploid evolution, the deletion rate was 2.1 x 10(-3) locus(-1) MY(-1) for single-copy loci and 1.0 x 10(-2) locus(-1) MY(-1) for loci in paralogous sets. Loci were duplicated with a rate of 2.9 x 10(-3) locus(-1) MY(-1) during diploid evolution. During polyploid evolution, locus deletion and locus duplication rates were 1.8 x 10(-2) and 1.8 x 10(-3) locus(-1) MY(-1), respectively. Locus deletion and duplication rates correlated positively with the distance of the locus from the centromere and the recombination rate during diploid evolution. The functions of deleted and duplicated loci were inferred to gain insight into the surprisingly high rate of deletions of loci present apparently only once in a genome. The significance of these findings for genome evolution at the diploid and polyploid level is discussed.

A 2500-locus bin map of wheat homoeologous group 5 provides insights on gene distribution and colinearity with rice

We constructed high-density deletion bin maps of wheat chromosomes 5A, 5B, and 5D, including 2338 loci mapped with 1052 EST probes and 217 previously mapped loci (total 2555 loci). This information was combined to construct a consensus chromosome bin map of group 5 including 24 bins. A relatively higher number of loci were mapped on chromosome 5B (38%) compared to 5A (34%) and 5D (28%). Differences in the levels of polymorphism among the three chromosomes were partially responsible for these differences. A higher number of duplicated loci was found on chromosome 5B (42%). Three times more loci were mapped on the long arms than on the short arms, and a significantly higher number of probes, loci, and duplicated loci were mapped on the distal halves than on the proximal halves of the chromosome arms. Good overall colinearity was observed among the three homoeologous group 5 chromosomes, except for the previously known 5AL/4AL translocation and a putative small pericentric inversion in chromosome 5A. Statistically significant colinearity was observed between low-copy-number ESTs from wheat homoeologous group 5 and rice chromosomes 12 (88 ESTs), 9 (72 ESTs), and 3 (84 ESTs).

Sequence composition, organization, and evolution of the core Triticeae genome

We investigated the composition and the basis of genome expansion in the core Triticeae genome using Aegilops tauschii, the D-genome donor of bread wheat. We sequenced an unfiltered genomic shotgun (trs) and a methylation-filtration (tmf) library of A. tauschii, and analyzed wheat expressed sequence tags (ESTs) to estimate the expression of genes and transposable elements (TEs). The sampled D-genome sequences consisted of 91.6% repetitive elements, 2.5% known genes, and 5.9% low-copy sequences of unknown function. TEs constituted 68.2% of the D-genome compared with 50% in maize and 14% in rice. The DNA transposons constituted 13% of the D-genome compared with 2% in maize. TEs were methylated unevenly within and among elements and families, and most were transcribed which contributed to genome expansion in the core Triticeae genome. The copy number of a majority of repeat families increased gradually following polyploidization. Certain TE families occupied discrete chromosome territories. Nested insertions and illegitimate recombination occurred extensively between the TE families, and a majority of the TEs contained internal deletions. The GC content varied significantly among the three sequence sets examined ranging from 42% in tmf to 46% in trs and 52% in the EST. Based on enrichment of genic sequences, methylation-filtration offers one option, although not as efficient as in maize, for isolating gene-rich regions from the large genome of wheat.

Genomic organization of the complex alpha-gliadin gene loci in wheat

To better understand the molecular evolution of the large alpha-gliadin gene family, a half-million bacterial artificial chromosome (BAC) library clones from tetraploid durum wheat, Triticum turgidum ssp. durum (2n = 4x = 28, genome AB), were screened for large genomic segments carrying the alpha-gliadin genes of the Gli-2 loci on the group 6 homoeologous chromosomes. The resulting 220 positive BAC clones--each containing between one and four copies of alpha-gliadin sequences--were fingerprinted for contig assembly to produce contiguous chromosomal regions covering the Gli-2 loci. While contigs consisting of as many as 21 BAC clones and containing up to 17 alpha-gliadin genes were formed, many BAC clones remained as singletons. The accuracy of the order of BAC clones in the contigs was verified by Southern hybridization analysis of the BAC fingerprints using an alpha-gliadin probe. These results indicate that alpha-gliadin genes are not evenly dispersed in the Gli-2 locus regions. Hybridization of these BACs with probes for long terminal repeat retrotransposons was used to determine the abundance and distribution of repetitive DNA in this region. Sequencing of BAC ends indicated that 70% of the sequences were significantly similar to different classes of retrotransposons, suggesting that these elements are abundant in this region. Several mechanisms underlying the dynamic evolution of the Gli-2 loci are discussed.

Sequencing of the Triticum monococcum Hardness locus reveals good microcolinearity with rice

The Hardness ( Ha) locus on chromosome 5D is the main determinant of grain texture in hexaploid wheat. The related genes Puroindoline-a and -b ( Pina-D1, Pinb-D1) and Grain Softness Protein ( Gsp-D1) are tightly linked at this locus. Mutations in the Pina-D1 and Pinb-D1 genes are associated with increased grain hardness. We report here the complete sequence of a 101-kb BAC clone from Triticum monococcum (A(m ) genome) which includes these three genes, and its comparison with the orthologous region in rice. The genes Gsp-A(m) 1, Pina-A(m) 1 and Pinb-A(m) 1 are separated by 37 kb and 32 kb, respectively, and are organized in the same transcriptional orientation. Four additional genes, including a pair of duplicated genes, were identified upstream of Gsp-A(m) 1 within a high-density gene island. These additional genes were found in the same order and orientation, and the same relative distances apart as similar genes previously annotated on rice chromosome 12. An interesting discovery was a small unannotated putative rice gene that was similar to the Gsp-A(m) 1 gene of T. monococcum (65% similarity at the protein level), and that was disposed in the same orientation, and located in the same position relative to the other orthologous genes. The high gene density observed in this BAC (1 gene per 14 kb) was expected for a distal chromosome region, but the level of microcolinearity with rice was higher than that reported in similar distal regions of other wheat chromosomes. Most of the BAC sequence (40%) was represented by repetitive elements, mainly concentrated in regions adjacent to the genes Pina-A(m) 1 and Pinb-A(m) 1. Rearrangements among these repetitive elements might provide an explanation for the frequent deletions observed at this locus in the genomes of the polyploid wheat species.

A unique set of 11,008 onion expressed sequence tags reveals expressed sequence and genomic differences between the monocot orders Asparagales and Poales

Enormous genomic resources have been developed for plants in the monocot order Poales; however, it is not clear how representative the Poales are for the monocots as a whole. The Asparagales are a monophyletic order sister to the lineage carrying the Poales and possess economically important plants such as asparagus, garlic, and onion. To assess the genomic differences between the Asparagales and Poales, we generated 11,008 unique ESTs from a normalized cDNA library of onion. Sequence analyses of these ESTs revealed microsatellite markers, single nucleotide polymorphisms, and homologs of transposable elements. Mean nucleotide similarity between rice and the Asparagales was 78% across coding regions. Expressed sequence and genomic comparisons revealed strong differences between the Asparagales and Poales for codon usage and mean GC content, GC distribution, and relative GC content at each codon position, indicating that genomic characteristics are not uniform across the monocots. The Asparagales were more similar to eudicots than to the Poales for these genomic characteristics.

Identification and genomic distribution of gypsy like retrotransposons in Citrus and Poncirus

Transposable elements might be importantly involved in citrus genetic instability and genome evolution. The presence of gypsy like retrotransposons, their heterogeneity and genomic distribution in Citrus and Poncirus, have been investigated. Eight clones containing part of the POL coding region of gypsy like retrotransposons have been isolated from a commercial variety of Citrus clementina, one of the few sexual species in Citrus. Four of the eight clones might correspond to active elements given that they present all the conserved motifs described in the literature as essential for activity, no in-frame stop codon and no frame-shift mutation. High homology has been found between some of these citrus elements and retroelements within a resistance-gene cluster from potato, another from Poncirus trifoliata and two putative resistance polyproteins from rice. Nested copies of gypsy like elements are scattered along the Citrus and Poncirus genomes. The results on genomic distribution show that these elements were introduced before the divergence of both genera and evolved separately thereafter. IRAPs based on gypsy and copia types of retrotransposons seem to distribute differently, therefore gypsy based IRAPs prove a new, complementary set of molecular markers in Citrus to study and map genetic variability, especially for disease resistance. Similarly to copia-derived IRAPs, the number of copies and heterozygosity values found for gypsy derived IRAPs are lower in Poncirus than in Citrus aurantium, which is less apomictic and the most usual rootstock for clementines until 1970.

Annotation and BAC/PAC localization of nonredundant ESTs from drought-stressed seedlings of an indica rice

To decipher the genes associated with drought stress response and to identify novel genes in rice, we utilized 1540 high-quality expressed sequence tags (ESTs) for functional annotation and mapping to rice genomic sequences. These ESTs were generated earlier by 3'-end single-pass sequencing of 2000 cDNA clones from normalized cDNA libraries constructed form drought-stressed seedlings of an indica rice. A rice UniGene set of 1025 transcripts was constructed from this collection through the BLASTN algorithm. Putative functions of 559 nonredundant ESTs were identified by BLAST similarity search against public databases. Putative functions were assigned at a stringency E value of 10(-6) in BLASTN and BLASTX algorithms. To understand the gene structure and function further, we have utilized the publicly available finished and unfinished rice BAC/PAC (BAC, bacterial artificial chromosome; PAC, P1 artificial chromosome) sequences for similarity search using the BLASTN algorithm. Further, 603 nonredundant ESTs have been mapped to BAC/PAC clones. BAC clones were assigned by a homology of above 95% identity along 90% of EST sequence length in the aligned region. In all, 700 ESTs showed rice EST hits in GenBank. Of the 325 novel ESTs, 128 were localized to BAC clones. In addition, 127 ESTs with identified putative functions but with no homology in IRGSP (International Rice Genome Sequencing Program) BAC/PAC sequences were mapped to the Chinese WGS (whole genome shotgun contigs) draft sequence of the rice genome. Functional annotation uncovered about a hundred candidate ESTs associated with abiotic stress in rice and Arabidopsis that were previously reported based on microarray analysis and other studies. This study is a major effort in identifying genes associated with drought stress response and will serve as a resource to rice geneticists and molecular biologists.