Proliferation of Ty3/gypsy-like retrotransposons in hybrid sunflower taxa inferred from phylogenetic data

<p>Abstract</p> <p>Background</p> <p>Long terminal repeat (LTR) retrotransposons are a class of mobile genetic element capable of autonomous transposition via an RNA intermediate. Their large size and proliferative ability make them important contributors to genome size evolution, especially in plants, where they can reach exceptionally high copy numbers and contribute substantially to variation in genome size even among closely related taxa. Using a phylogenetic approach, we characterize dynamics of proliferation events of <it>Ty3/gypsy</it>-like LTR retrotransposons that led to massive genomic expansion in three <it>Helianthus </it>(sunflower) species of ancient hybrid origin. The three hybrid species are independently derived from the same two parental species, offering a unique opportunity to explore patterns of retrotransposon proliferation in light of reticulate evolutionary events in this species group.</p> <p>Results</p> <p>We demonstrate that <it>Ty3/gypsy</it>-like retrotransposons exist as multiple well supported sublineages in both the parental and hybrid derivative species and that the same element sublineage served as the source lineage of proliferation in each hybrid species' genome. This inference is based on patterns of species-specific element numerical abundance within different phylogenetic sublineages as well as through signals of proliferation events present in the distributions of element divergence values. Employing methods to date paralogous sequences within a genome, proliferation events in the hybrid species' genomes are estimated to have occurred approximately 0.5 to 1 million years ago.</p> <p>Conclusion</p> <p>Proliferation of the same retrotransposon major sublineage in each hybrid species indicates that similar dynamics of element derepression and amplification likely occurred in each hybrid taxon during their formation. Temporal estimates of these proliferation events suggest an earlier origin for these hybrid species than previously supposed.</p

C. elegans Germ Cells Show Temperature and Age-Dependent Expression of Cer1, a Gypsy/Ty3-Related Retrotransposon

Virus-like particles (VLPs) have not been observed in Caenorhabditis germ cells, although nematode genomes contain low numbers of retrotransposon and retroviral sequences. We used electron microscopy to search for VLPs in various wild strains of Caenorhabditis, and observed very rare candidate VLPs in some strains, including the standard laboratory strain of C. elegans, N2. We identified the N2 VLPs as capsids produced by Cer1, a retrotransposon in the Gypsy/Ty3 family of retroviruses/retrotransposons. Cer1 expression is age and temperature dependent, with abundant expression at 15°C and no detectable expression at 25°C, explaining how VLPs escaped detection in previous studies. Similar age and temperature-dependent expression of Cer1 retrotransposons was observed for several other wild strains, indicating that these properties are common, if not integral, features of this retroelement. Retrotransposons, in contrast to DNA transposons, have a cytoplasmic stage in replication, and those that infect non-dividing cells must pass their genomic material through nuclear pores. In most C. elegans germ cells, nuclear pores are largely covered by germline-specific organelles called P granules. Our results suggest that Cer1 capsids target meiotic germ cells exiting pachytene, when free nuclear pores are added to the nuclear envelope and existing P granules begin to be removed. In pachytene germ cells, Cer1 capsids concentrate away from nuclei on a subset of microtubules that are exceptionally resistant to microtubule inhibitors; the capsids can aggregate these stable microtubules in older adults, which exhibit a temperature-dependent decrease in egg viability. When germ cells exit pachytene, the stable microtubules disappear and capsids redistribute close to nuclei that have P granule-free nuclear pores. This redistribution is microtubule dependent, suggesting that capsids that are released from stable microtubules transfer onto new, dynamic microtubules to track toward nuclei. These studies introduce C. elegans as a model to study the interplay between retroelements and germ cell biology

Mosquitoes LTR Retrotransposons: A Deeper View into the Genomic Sequence of Culex quinquefasciatus

A set of 67 novel LTR-retrotransposon has been identified by in silico analyses of the Culex quinquefasciatus genome using the LTR_STRUC program. The phylogenetic analysis shows that 29 novel and putatively functional LTR-retrotransposons detected belong to the Ty3/gypsy group. Our results demonstrate that, by considering only families containing potentially autonomous LTR-retrotransposons, they account for about 1% of the genome of C. quinquefasciatus. In previous studies it has been estimated that 29% of the genome of C. quinquefasciatus is occupied by mobile genetic elements

Contrasting Patterns of Transposable Element Insertions in Drosophila Heat-Shock Promoters

The proximal promoter regions of heat-shock genes harbor a remarkable number of P transposable element (TE) insertions relative to both positive and negative control proximal promoter regions in natural populations of Drosophila melanogaster. We have screened the sequenced genomes of 12 species of Drosophila to test whether this pattern is unique to these populations. In the 12 species' genomes, transposable element insertions are no more abundant in promoter regions of single-copy heat-shock genes than in promoters with similar or dissimilar architecture. Also, insertions appear randomly distributed across the promoter region, whereas insertions clustered near the transcription start site in promoters of single-copy heat-shock genes in D. melanogaster natural populations. Hsp70 promoters exhibit more TE insertions per promoter than all other genesets in the 12 species, similarly to in natural populations of D. melanogaster. Insertions in the Hsp70 promoter region, however, cluster away from the transcription start site in the 12 species, but near it in natural populations of D. melanogaster. These results suggest that D. melanogaster heat-shock promoters are unique in terms of their interaction with transposable elements, and confirm that Hsp70 promoters are distinctive in TE insertions across Drosophila

Expansion and subfunctionalisation of flavonoid 3',5'-hydroxylases in the grapevine lineage

<p>Abstract</p> <p>Background</p> <p>Flavonoid 3',5'-hydroxylases (F3'5'Hs) and flavonoid 3'-hydroxylases (F3'Hs) competitively control the synthesis of delphinidin and cyanidin, the precursors of blue and red anthocyanins. In most plants, <it>F3'5'H </it>genes are present in low-copy number, but in grapevine they are highly redundant.</p> <p>Results</p> <p>The first increase in <it>F3'5'H </it>copy number occurred in the progenitor of the eudicot clade at the time of the γ triplication. Further proliferation of <it>F3'5'H</it>s has occurred in one of the paleologous loci after the separation of Vitaceae from other eurosids, giving rise to 15 paralogues within 650 kb. Twelve reside in 9 tandem blocks of ~35-55 kb that share 91-99% identity. The second paleologous <it>F3'5'H </it>has been maintained as an orphan gene in grapevines, and lacks orthologues in other plants. Duplicate <it>F3'5'H</it>s have spatially and temporally partitioned expression profiles in grapevine. The orphan <it>F3'5'H </it>copy is highly expressed in vegetative organs. More recent duplicate <it>F3'5'H</it>s are predominately expressed in berry skins. They differ only slightly in the coding region, but are distinguished in the structure of the promoter. Differences in <it>cis</it>-regulatory sequences of promoter regions are paralleled by temporal specialisation of gene transcription during fruit ripening. Variation in anthocyanin profiles consistently reflects changes in the <it>F3'5'H </it>mRNA pool across different cultivars. More <it>F3'5'H </it>copies are expressed at high levels in grapevine varieties with 93-94% of 3'5'-OH anthocyanins. In grapevines depleted in 3'5'-OH anthocyanins (15-45%), fewer <it>F3'5'H </it>copies are transcribed, and at lower levels. Conversely, only two copies of the gene encoding the competing F3'H enzyme are present in the grape genome; one copy is expressed in both vegetative and reproductive organs at comparable levels among cultivars, while the other is transcriptionally silent.</p> <p>Conclusions</p> <p>These results suggest that expansion and subfunctionalisation of <it>F3'5'H</it>s have increased the complexity and diversification of the fruit colour phenotype among red grape varieties.</p

Evolutionary Dynamics of the Ty3/Gypsy LTR Retrotransposons in the Genome of Anopheles gambiae

Ty3/gypsy elements represent one of the most abundant and diverse LTR-retrotransposon (LTRr) groups in the Anopheles gambiae genome, but their evolutionary dynamics have not been explored in detail. Here, we conduct an in silico analysis of the distribution and abundance of the full complement of 1045 copies in the updated AgamP3 assembly. Chromosomal distribution of Ty3/gypsy elements is inversely related to arm length, with densities being greatest on the X, and greater on the short versus long arms of both autosomes. Taking into account the different heterochromatic and euchromatic compartments of the genome, our data suggest that the relative abundance of Ty3/gypsy LTRrs along each chromosome arm is determined mainly by the different proportions of heterochromatin, particularly pericentric heterochromatin, relative to total arm length. Additionally, the breakpoint regions of chromosomal inversion 2La appears to be a haven for LTRrs. These elements are underrepresented more than 7-fold in euchromatin, where 33% of the Ty3/gypsy copies are associated with genes. The euchromatin on chromosome 3R shows a faster turnover rate of Ty3/gypsy elements, characterized by a deficit of proviral sequences and the lowest average sequence divergence of any autosomal region analyzed in this study. This probably reflects a principal role of purifying selection against insertion for the preservation of longer conserved syntenyc blocks with adaptive importance located in 3R. Although some Ty3/gypsy LTRrs show evidence of recent activity, an important fraction are inactive remnants of relatively ancient insertions apparently subject to genetic drift. Consistent with these computational predictions, an analysis of the occupancy rate of putatively older insertions in natural populations suggested that the degenerate copies have been fixed across the species range in this mosquito, and also are shared with the sibling species Anopheles arabiensis

Evolution of Stress-Regulated Gene Expression in Duplicate Genes of Arabidopsis thaliana

Due to the selection pressure imposed by highly variable environmental conditions, stress sensing and regulatory response mechanisms in plants are expected to evolve rapidly. One potential source of innovation in plant stress response mechanisms is gene duplication. In this study, we examined the evolution of stress-regulated gene expression among duplicated genes in the model plant Arabidopsis thaliana. Key to this analysis was reconstructing the putative ancestral stress regulation pattern. By comparing the expression patterns of duplicated genes with the patterns of their ancestors, duplicated genes likely lost and gained stress responses at a rapid rate initially, but the rate is close to zero when the synonymous substitution rate (a proxy for time) is >∼0.8. When considering duplicated gene pairs, we found that partitioning of putative ancestral stress responses occurred more frequently compared to cases of parallel retention and loss. Furthermore, the pattern of stress response partitioning was extremely asymmetric. An analysis of putative cis-acting DNA regulatory elements in the promoters of the duplicated stress-regulated genes indicated that the asymmetric partitioning of ancestral stress responses are likely due, at least in part, to differential loss of DNA regulatory elements; the duplicated genes losing most of their stress responses were those that had lost more of the putative cis-acting elements. Finally, duplicate genes that lost most or all of the ancestral responses are more likely to have gained responses to other stresses. Therefore, the retention of duplicates that inherit few or no functions seems to be coupled to neofunctionalization. Taken together, our findings provide new insight into the patterns of evolutionary changes in gene stress responses after duplication and lay the foundation for testing the adaptive significance of stress regulatory changes under highly variable biotic and abiotic environments

Modes of Gene Duplication Contribute Differently to Genetic Novelty and Redundancy, but Show Parallels across Divergent Angiosperms

BACKGROUND: Both single gene and whole genome duplications (WGD) have recurred in angiosperm evolution. However, the evolutionary effects of different modes of gene duplication, especially regarding their contributions to genetic novelty or redundancy, have been inadequately explored. RESULTS: In Arabidopsis thaliana and Oryza sativa (rice), species that deeply sample botanical diversity and for which expression data are available from a wide range of tissues and physiological conditions, we have compared expression divergence between genes duplicated by six different mechanisms (WGD, tandem, proximal, DNA based transposed, retrotransposed and dispersed), and between positional orthologs. Both neo-functionalization and genetic redundancy appear to contribute to retention of duplicate genes. Genes resulting from WGD and tandem duplications diverge slowest in both coding sequences and gene expression, and contribute most to genetic redundancy, while other duplication modes contribute more to evolutionary novelty. WGD duplicates may more frequently be retained due to dosage amplification, while inferred transposon mediated gene duplications tend to reduce gene expression levels. The extent of expression divergence between duplicates is discernibly related to duplication modes, different WGD events, amino acid divergence, and putatively neutral divergence (time), but the contribution of each factor is heterogeneous among duplication modes. Gene loss may retard inter-species expression divergence. Members of different gene families may have non-random patterns of origin that are similar in Arabidopsis and rice, suggesting the action of pan-taxon principles of molecular evolution. CONCLUSION: Gene duplication modes differ in contribution to genetic novelty and redundancy, but show some parallels in taxa separated by hundreds of millions of years of evolution

BEL/Pao retrotransposons in metazoan genomes

Background Long terminal repeat (LTR) retrotransposons are a widespread kind of transposable element present in eukaryotic genomes. They are a major factor in genome evolution due to their ability to create large scale mutations and genome rearrangements. Compared to other transposable elements, little attention has been paid to elements belonging to the metazoan BEL/Pao subclass of LTR retrotransposons. No comprehensive characterization of these elements is available so far. The aim of this study was to describe all BEL/Pao elements in a set of 62 sequenced metazoan genomes, and to analyze their phylogenetic relationship. Results We identified a total of 7,861 BEL/Pao elements in 53 of our 62 study genomes. We identified BEL/Pao elements in 20 genomes where such elements had not been found so far. Our analysis shows that BEL/Pao elements are the second-most abundant class of LTR retrotransposons in the genomes we study, more abundant than Ty1/Copia elements, and second only to Ty3/Gypsy elements. They occur in multiple phyla, including basal metazoan phyla, suggesting that BEL/Pao elements arose early in animal evolution. We confirm findings from previous studies that BEL/Pao elements do not occur in mammals. The elements we found can be grouped into more than 1725 families, 1623 of which are new, previously unknown families. These families fall into seven superfamilies, only five of which have been characterized so far. One new superfamily is a major subdivision of the Pao superfamily which we propose to call Dan, because it is restricted to the genome of the zebrafish Danio rerio. The other new superfamily comprises 83 elements and is restricted to lower aquatic eumetazoans. We propose to call this superfamily Flow. BEL/Pao elements do not show any signs of recent horizontal gene transfer between distantly related species. Conclusions In sum, our analysis identifies thousands of new BEL/Pao elements and provides new insights into their distribution, abundance, and evolution