Comparative analysis of methylthioalkylmalate synthase (MAM) gene family and flanking DNA sequences in Brassica oleracea and Arabidopsis thaliana

Gene BoGSL-PRO is associated with presence of 3-carbon side-chain glucosinolates (GSL). This gene is a member of the methylthioalkylmalate synthase (MAM) gene family. A BAC clone of Brassica oleracea, B21F5, containing this gene, was sequenced, annotated and compared to its corresponding region in Arabidopsis thaliana. Twelve protein-coding genes and 10 transposable elements were found in this clone. The corresponding region in A. thaliana chromosome I has 14 genes and no transposable elements. Analysis of MAM gene family in both species, which also include genes controlling 4-carbon side-chain GSL, separated the genes in two groups based on exon numbers and function. Phylogenetic analysis of the amino acid sequences encoded by these genes suggest that these two groups were produced by a duplication that must have occurred before the divergence of the Rosid and Asterid lineages of angiosperms. Comparison with putative orthologs from several prokaryotes further suggest that the members of the gene family with 10 exons, which encode proteins involved in 4-carbon side-chain GSL biosynthesis, were derived via truncation of the 3′ end from ancestral genes more similar in length to those with 12 exons, which encode proteins involved in 3-carbon side-chain GSL biosynthesis. Lower gene density in B. oleracea compared to A. thaliana is due in part to presence of transposable elements (TE) mostly in inter-genic regions.     

Genome ecosystem and transposable elements species

Transposable elements are known to be “selfish DNA” sequences able to spread and be maintained in all genomes analyzed so far. Their evolution depends on the interaction they have with the other components of the genome, including genes and other transposable elements. These relationships are complex and have often been compared to those of species living and competing in an ecosystem. The aim of this current work is a proposition to fill the conceptual gap existing between genome biology and ecology, assuming that genomic components, such as transposable elements families, can be compared to species interacting in an ecosystem. Using this framework, some of the main models defined in the population genetics of transposable elements can then been reformulated, and some new kinds of realistic relationships, such as symbiosis between different genomic components, can then be modelled and explored.     

The splicing of transposable elements and its role in intron evolution

Recent studies have demonstrated that transposable elements in maize and Drosophila are spliced from pre-mRNA. These transposable element introns represent the first examples of recent addition of introns into nuclear genes. The eight reported examples of transposable element splicing include members of the maize Ac/Ds and Spm/dSpm and the Drosophila P and 412 element families. The details of the splicing of these transposable elements and their relevance to models of intron origin are discussed.     

Diverse origins of waxy foxtail millet crops in East and Southeast Asia mediated by multiple transposable element insertions

The naturally occurring waxy and low-amylose variants of foxtail millet and other cereals, like rice and barley, originated in East and Southeast Asia under human selection for sticky foods. Mutations in the GBSS1 gene for granule-bound starch synthase 1 are known to be associated with these traits. We have analyzed the gene in foxtail millet, and found that, in this species, these traits were originated by multiple independent insertions of transposable elements and by subsequent secondary insertions into these elements or deletion of parts of the elements. The structural analysis of transposable elements inserted in the GBSS1 gene revealed that the non-waxy was converted to the low-amylose phenotype once, while shifts from non-waxy to waxy occurred three times, from low amylose to waxy once and from waxy to low amylose once. The present results, and the geographical distribution of different waxy molecular types, strongly suggest that these types originated independently and were dispersed into their current distribution areas. The patterns of GBSS1 variation revealed here suggest that foxtail millet may serve as a key to solving the mystery of the origin of waxy-type cereals in Asia. The GBSS1 gene in foxtail millet provides a new example of the evolution of a gene involved in the processes of domestication and its post-domestication fate under the influence of human selection. Electronic Supplementary Material Supplementary material is available for this article at     

Rapid evolution and complex structural organization in genomic regions harboring multiple prolamin genes in the polyploid wheat genome

Genes encoding wheat prolamins belong to complicated multi-gene families in the wheat genome. To understand the structural complexity of storage protein loci, we sequenced and analyzed orthologous regions containing both gliadin and LMW-glutenin genes from the A and B genomes of a tetraploid wheat species, Triticum turgidum ssp. durum. Despite their physical proximity to one another, the gliadin genes and LMW-glutenin genes are organized quite differently. The gliadin genes are found to be more clustered than the LMW-glutenin genes which are separated from each other by much larger distances. The separation of the LMW-glutenin genes is the result of both the insertion of large blocks of repetitive DNA owing to the rapid amplification of retrotransposons and the presence of genetic loci interspersed between them. Sequence comparisons of the orthologous regions reveal that gene movement could be one of the major factors contributing to the violation of microcolinearity between the homoeologous A and B genomes in wheat. The rapid sequence rearrangements and differential insertion of repetitive DNA has caused the gene islands to be not conserved in compared regions. In addition, we demonstrated that the i-type LMW-glutenin originated from a deletion of 33-bps in the 5′ coding region of the m-type gene. Our results show that multiple rounds of segmental duplication of prolamin genes have driven the amplification of the ω-gliadin genes in the region; such segmental duplication could greatly increase the repetitive DNA content in the genome depending on the amount of repetitive DNA present in the original duplicate region. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Shuangcheng Gao and Yong Qiang Gu contributed equally to the work.     

TARGeT: a web-based pipeline for retrieving and characterizing gene and transposable element families from genomic sequences

Gene families compose a large proportion of eukaryotic genomes. The rapidly expanding genomic sequence database provides a good opportunity to study gene family evolution and function. However, most gene family identification programs are restricted to searching protein databases where data are often lagging behind the genomic sequence data. Here, we report a user-friendly web-based pipeline, named TARGeT (Tree Analysis of Related Genes and Transposons), which uses either a DNA or amino acid ‘seed’ query to: (i) automatically identify and retrieve gene family homologs from a genomic database, (ii) characterize gene structure and (iii) perform phylogenetic analysis. Due to its high speed, TARGeT is also able to characterize very large gene families, including transposable elements (TEs). We evaluated TARGeT using well-annotated datasets, including the ascorbate peroxidase gene family of rice, maize and sorghum and several TE families in rice. In all cases, TARGeT rapidly recapitulated the known homologs and predicted new ones. We also demonstrated that TARGeT outperforms similar pipelines and has functionality that is not offered elsewhere.     

Characterization and potential evolutionary impact of transposable elements in the genome of Cochliobolus heterostrophus

Background

Cochliobolus heterostrophus is a dothideomycete that causes Southern Corn Leaf Blight disease. There are two races, race O and race T that differ by the absence (race O) and presence (race T) of ~ 1.2-Mb of DNA encoding genes responsible for the production of T-toxin, which makes race T much more virulent than race O. The presence of repetitive elements in fungal genomes is considered to be an important source of genetic variability between different species.

Results

A detailed analysis of class I and II TEs identified in the near complete genome sequence of race O was performed. In total in race O, 12 new families of transposons were identified. In silico evidence of recent activity was found for many of the transposons and analyses of expressed sequence tags (ESTs) demonstrated that these elements were actively transcribed. Various potentially active TEs were found near coding regions and may modify the expression and structure of these genes by acting as ectopic recombination sites. Transposons were found on scaffolds carrying polyketide synthase encoding genes, responsible for production of T-toxin in race T. Strong evidence of ectopic recombination was found, demonstrating that TEs can play an important role in the modulation of genome architecture of this species. The Repeat Induced Point mutation (RIP) silencing mechanism was shown to have high specificity in C. heterostrophus, acting only on transposons near coding regions.

Conclusions

New families of transposons were identified. In C. heterostrophus, the RIP silencing mechanism is efficient and selective. The co-localization of effector genes and TEs, therefore, exposes those genes to high rates of point mutations. This may accelerate the rate of evolution of these genes, providing a potential advantage for the host. Additionally, it was shown that ectopic recombination promoted by TEs appears to be the major event in the genome reorganization of this species and that a large number of elements are still potentially active. So, this study provides information about the potential impact of TEs on the evolution of C. heterostrophus.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-536) contains supplementary material, which is available to authorized users.     

Can transposable element copy number distribution parameters be estimated from natural populations of Drosophila melanogaster?

The copy frequency distribution of a transposable element family in a Drosophila melanogaster natural population is generally characterised by the values of the Charlesworths' model parameters α and β (Charlesworth & Charlesworth, 1983). The estimation of these parameters is made using the observed distribution of the occupied sites in a population sample. Several results have been interpreted as due either to the influence of stochastic factors or to deterministic factors (transposition, excision, selection…). The accuracy of this method was tested by estimations performed on samples from simulated populations. The results show that with the sample size usually used for natural population studies, the confidence intervals are too large to reasonably deduce either the element copy number distribution or the values of transposition and excision rate and selective coefficients.     

Evidence for interspecific transfer of the transposable element mariner betweenDrosophila andZaprionus

Summary The transposable element mariner occurs widely in themelanogaster species group ofDrosophila. However, in drosophilids outside of themelanogaster species group, sequences showing strong DNA hybridization with mariner are found only in the genusZaprionus. the mariner sequence obtained fromZaprionus tuberculatus is 97% identical with that fromDrosophila mauritiana, a member of themelanogaster species subgroup, whereas a mariner sequence isolated fromDrosophila tsacasi is only 92% identical with that fromD. mauritiana. BecauseD. tsacasi is much more closely related toD. mauritiana than isZaprionus, the presence of mariner inZaprionus may result from horizontal transfer. In order to confirm lack of a close phylogenetic relationship between the genusZaprionus and themelanogaster species group, we compared the alcohol dehydrogenase (Adh) sequences among these species. The results show that the coding region of Adh is only 82% identical betweenZ. tuberculatus andD. mauritiana, as compared with 90% identical betweenD. tsacasi andD. mauritiana. Furthermore, the mariner gene phylogeny obtained by maximum likelihood and maximum parsimony analyses is discordant with the species phylogeny estimated by using the Adh genes. The only inconsistency in the mariner gene phylogeny is in the placement of theZaprionus mariner sequence, which clusters with mariner fromDrosophila teissieri andDrosophila yakuba in themelanogaster species subgroup. These results strongly suggest horizontal transfer.     

Next-generation sequencing of flow-sorted wheat chromosome 5D reveals lineage-specific translocations and widespread gene duplications

Background

The ~17 Gb hexaploid bread wheat genome is a high priority and a major technical challenge for genomic studies. In particular, the D sub-genome is relatively lacking in genetic diversity, making it both difficult to map genetically, and a target for introgression of agriculturally useful traits. Elucidating its sequence and structure will therefore facilitate wheat breeding and crop improvement.

Results

We generated shotgun sequences from each arm of flow-sorted Triticum aestivum chromosome 5D using 454 FLX Titanium technology, giving 1.34× and 1.61× coverage of the short (5DS) and long (5DL) arms of the chromosome respectively. By a combination of sequence similarity and assembly-based methods, ~74% of the sequence reads were classified as repetitive elements, and coding sequence models of 1314 (5DS) and 2975 (5DL) genes were generated. The order of conserved genes in syntenic regions of previously sequenced grass genomes were integrated with physical and genetic map positions of 518 wheat markers to establish a virtual gene order for chromosome 5D.

Conclusions

The virtual gene order revealed a large-scale chromosomal rearrangement in the peri-centromeric region of 5DL, and a concentration of non-syntenic genes in the telomeric region of 5DS. Although our data support the large-scale conservation of Triticeae chromosome structure, they also suggest that some regions are evolving rapidly through frequent gene duplications and translocations.

Sequence accessions

EBI European Nucleotide Archive, Study no. ERP002330

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1080) contains supplementary material, which is available to authorized users.     

<Emphasis Type="Italic">Helitron</Emphasis> mediated amplification of cytochrome P450 monooxygenase gene in maize

The mass movement of gene sequences by Helitrons has significantly contributed to the lack of gene collinearity reported between different maize inbred lines. However, Helitron captured-genes reported to date represent truncated versions of their progenitor genes. In this report, we provide evidence that maize CYP72A27-Zm gene represents a cytochrome P450 monooxygenase (P450) gene recently captured by a Helitron and transposed into an Opie-2 retroposon. The four exons of the CYP72A27 gene contained within the element contain a putative open reading frame (ORF) for 428 amino acid residues. We provide evidence that Helitron captured CYP72A27-Zm is transcribed. To identify the progenitor gene and the evolutionary time of capture, we searched the plant genome database and discovered other closely related CYP72A27-Zm genes in maize and grasses. Our analysis indicates that CYP72A27-Zm represents an almost complete copy of maize CYP72A26-Zm gene captured by a Helitron about 3.1 million years ago (mya). The Helitron-captured gene then duplicated twice, approximately 1.5-1.6 mya giving rise to CYP72A36-Zm and CYP72A37-Zm. These data provide evidence that Helitrons can capture and mobilize intact genes that are transcribed and potentially encode biologically relevant proteins.     

Somatic mutations caused by excision of the transposable element, Tpn1, from the DFR gene for pigmentation in sub-epidermal layer of periclinally chimeric flowers of Japanese morning glory and their germinal transmission to their progeny

Pigmentation in flowers of Japanese morning glory is intense in the epidermal layer, lighter in the subepidermis, and much lighter in the internal tissues; by contrast coloration in stems occurs only in the sub-epidermal layer. The a-3 ^f mutant of Japanese morning glory bears white flowers with normal-colored flecks and sectors, and its variegation also occurs in leaves and stems. The mutable line can produce chimeric flowers pigmented uniformly in the sub-epidermal tissue and variegated in the epidermal layer, and stems of these flowers are also pigmented. Since they give selfed progeny that segregate to give a ratio of three germinal revertants bearing fully colored flowers to one flecked mutant, it has been [OR Imai (1934) has] postulated that somatic mutations in the sub-epidermal layer can be transmitted to the next generation and that the germ cells in the reproductive organs must form from the cells of the sub-epidermal layer. Recently, we found that the 6.4-kb En/Spm-related transposable element, Tpn1, resides within the DFR-B gene for anthocyanin biosynthesis in the mutable a-3 ^f line. To test whether somatic mutations caused by Tpn1 excision from the DFR-B gene in the subepidermis of periclinally chimeric flowers are transmissible to their progeny, we have examined the structure of the DFR-B region in the germinal revertants derived from the chimeric flowers and compared the sequences generated by the somatic excision of Tpn1 in periclinally chimeric flowers with those in their germinal revertants. Our results confirm that somatic mutations caused by Tpn1 excision from the DFR-B gene in the sub-epidermal tissue of chimeric flowers can be transmitted to their progeny, which results in the generation of germinal revertants.     

Merlin, a new superfamily of DNA transposons identified in diverse animal genomes and related to bacterial IS1016 insertion sequences

Several new families of DNA transposons were identified by computer-assisted searches in a wide range of animal species that includes nematodes, flat worms, mosquitoes, sea squirt, zebrafish, and humans. Many of these elements have coding capacity for transposases, which are related to each other and to those encoded by the IS1016 group of bacterial insertion sequences. Although these transposases display a motif similar to the DDE motif found in many transposases and integrases, they cannot be directly allied to any of the previously described eukaryotic transposases. Other common features of the new eukaryotic and bacterial transposons include similarities in their terminal inverted repeats and 8-bp or 9-bp target-site duplications. Together, these data indicate that these elements belong to a new superfamily of DNA transposons, called Merlin/IS1016, which is common in many eubacterial and animal genomes. We also present evidence that these transposons have been recently active in several animal species. This evidence is particularly strong in the parasitic blood fluke Schistosoma mansoni, in which Merlin is also the first described DNA transposon family.     

Influence of the transposable element neighborhood on human gene expression in normal and tumor tissues

Transposable elements (TEs) are genomic sequences able to replicate themselves, and to move from one chromosomal position to another within the genome. Many TEs contain their own regulatory regions, which means that they may influence the expression of neighboring genes. TEs may also be activated and transcribed in various cancers. We therefore tested whether gene expression in normal and tumor tissues is influenced by the neighboring TEs. To do this, we associated all human genes to the nearest TEs. We analyzed the expression of these genes in normal and tumor tissues using SAGE and EST data, and related this to the presence and type of TEs in their vicinity. We confirmed that TEs tend to be located in antisense orientation relative to their hosting genes. We found that the average number of tissues where a gene is expressed varies depending on the type of TEs located near the gene, and that the difference in expression level between normal and tumor tissues is greatest for genes that host SINE elements. This deregulation increases with the number of SINE copies in the gene vicinity. This suggests that SINE elements might contribute to the cascade of gene deregulation in cancer cells.     

Reconstruction of ancient genome and gene order from complete microbial genome sequences

Microbial genome sequences provide us with the fossil records for inferring their origination and evolution. Assuming that current microbial genomes are the evolutionary results of ancient genomes or fragments and the neighboring genes in ancient genomes are more likely neighbors in current genomes, in this paper we proposed a paleontological algorithm and assembled the orthologous gene groups from 66 complete and current microbial genome sequences into a pseudo-ancient genome, which consists of continuous fragments of various sizes. We performed bootstrap resampling and correlation analyses and the results showed that the assembled ancient genome and fragments are statistically significant and the genes of the same fragment are inherently related and likely derived from common ancestors. This method provides a new computational tool for studying microbial genome structure and evolution.     

Molecular and functional characterization of Slide, anAc-like autonomous transposable element from tobacco

A new transposable element of tobacco, Slide, was isolated from thetl mutant line, which shows somatic instability, after its transposition into a locus encoding nitrate reductase (NR). The Slide-124 element is 3733 bp long and its coding sequences show similarities with conserved domains of the transposases ofAc, Tam3 andhobo. Excision from the NR locus is detectable in somatic leaf tissues and Slide mobility is triggered by in vitro tissue culture. Slide excision events create footprints similar to those left byAc and Tam3. Tobacco lines derived from thetl mutant line seem characterized by unmethylated copies of a few members of the highly repetitive Slide family. Slide mobility was monitored in transient expression assays. In wild-type tobacco protoplasts, the complete Slide element, as well as a defective copy, is able to excise. The complete Slide element, but not the defective version, is able to excise in protoplasts of the heterologous species lettuce (Lactuca sativa). These results show that Slide carries the functions required for its own mobility, and represents the first autonomousAc-like element characterized inSolanaceae species.