Short interspersed nuclear elements (SINEs) are abundant in Solanaceae and have a family‐specific impact on gene structure and genome organization

Short interspersed nuclear elements (SINEs) are highly abundant non‐autonomous retrotransposons that are widespread in plants. They are short in size, non‐coding, show high sequence diversity, and are therefore mostly not or not correctly annotated in plant genome sequences. Hence, comparative studies on genomic SINE populations are rare. To explore the structural organization and impact of SINEs, we comparatively investigated the genome sequences of the Solanaceae species potato (Solanum tuberosum), tomato (Solanum lycopersicum), wild tomato (Solanum pennellii), and two pepper cultivars (Capsicum annuum). Based on 8.5 Gbp sequence data, we annotated 82 983 SINE copies belonging to 10 families and subfamilies on a base pair level. Solanaceae SINEs are dispersed over all chromosomes with enrichments in distal regions. Depending on the genome assemblies and gene predictions, 30% of all SINE copies are associated with genes, particularly frequent in introns and untranslated regions (UTRs). The close association with genes is family specific. More than 10% of all genes annotated in the Solanaceae species investigated contain at least one SINE insertion, and we found genes harbouring up to 16 SINE copies. We demonstrate the involvement of SINEs in gene and genome evolution including the donation of splice sites, start and stop codons and exons to genes, enlargement of introns and UTRs, generation of tandem‐like duplications and transduction of adjacent sequence regions.     

Evolutionary modes of emergence of short interspersed nuclear element (SINE) families in grasses

Short interspersed nuclear elements (SINEs) are non‐autonomous transposable elements which are propagated by retrotransposition and constitute an inherent part of the genome of most eukaryotic species. Knowledge of heterogeneous and highly abundant SINEs is crucial for de novo (or improvement of) annotation of whole genome sequences. We scanned Poaceae genome sequences of six important cereals (Oryza sativa, Triticum aestivum, Hordeum vulgare, Panicum virgatum, Sorghum bicolor, Zea mays) and Brachypodium distachyon to examine the diversity and evolution of SINE populations. We comparatively analyzed the structural features, distribution, evolutionary relation and abundance of 32 SINE families and subfamilies within grasses, comprising 11 052 individual copies. The investigation of activity profiles within the Poaceae provides insights into their species‐specific diversification and amplification. We found that Poaceae SINEs (PoaS) fall into two length categories: simple SINEs of up to 180 bp and dimeric SINEs larger than 240 bp. Detailed analysis at the nucleotide level revealed that multimerization of related and unrelated SINE copies is an important evolutionary mechanism of SINE formation. We conclude that PoaS families diversify by massive reshuffling between SINE families, likely caused by insertion of truncated copies, and provide a model for this evolutionary scenario. Twenty‐eight of 32 PoaS families and subfamilies show significant conservation, in particular either in the 5′ or 3′ regions, across Poaceae species and share large sequence stretches with one or more other PoaS families.     

The conserved 3′ Angio‐domain defines a superfamily of short interspersed nuclear elements (SINEs) in higher plants

Repetitive sequences are ubiquitous components of eukaryotic genomes affecting genome size and evolution as well as gene regulation. Among them, short interspersed nuclear elements (SINEs) are non‐coding retrotransposons usually shorter than 1000 bp. They contain only few short conserved structural motifs, in particular an internal promoter derived from cellular RNAs and a mostly AT‐rich 3′ tail, whereas the remaining regions are highly variable. SINEs emerge and vanish during evolution, and often diversify into numerous families and subfamilies that are usually specific for only a limited number of species. In contrast, at the 3′ end of multiple plant SINEs we detected the highly conserved ‘Angio‐domain’. This 37 bp segment defines the Angio‐SINE superfamily, which encompasses 24 plant SINE families widely distributed across 13 orders within the plant kingdom. We retrieved 28 433 full‐length Angio‐SINE copies from genome assemblies of 46 plant species, frequently located in genes. Compensatory mutations in and adjacent to the Angio‐domain imply selective restraints maintaining its RNA structure. Angio‐SINE families share segmental sequence similarities, indicating a modular evolution with strong Angio‐domain preservation. We suggest that the conserved domain contributes to the evolutionary success of Angio‐SINEs through either structural interactions between SINE RNA and proteins increasing their transpositional efficiency, or by enhancing their accumulation in genes.     

Molecular characterization,genomic distribution and evolutionary dynamics of Short INterspersed Elements in the termite genome

Short INterspersed Elements (SINEs) in invertebrates, and especially in animal inbred genomes such that of termites, are poorly known; in this paper we characterize three new SINE families (Talub, Taluc and Talud) through the analyses of 341 sequences, either isolated from the Reticulitermes lucifugus genome or drawn from EST Genbank collection. We further add new data to the only isopteran element known so far, Talua. These SINEs are tRNA-derived elements, with an average length ranging from 258 to 372 bp. The tails are made up by poly(A) or microsatellite motifs. Their copy number varies from 7.9 × 10³ to 10⁵ copies, well within the range observed for other metazoan genomes. Species distribution, age and target site duplication analysis indicate Talud as the oldest, possibly inactive SINE originated before the onset of Isoptera (~150 Myr ago). Taluc underwent to substantial sequence changes throughout the evolution of termites and data suggest it was silenced and then re-activated in the R. lucifugus lineage. Moreover, Taluc shares a conserved sequence block with other unrelated SINEs, as observed for some vertebrate and cephalopod elements. The study of genomic environment showed that insertions are mainly surrounded by microsatellites and other SINEs, indicating a biased accumulation within non-coding regions. The evolutionary dynamics of Talu~ elements is explained through selective mechanisms acting in an inbred genome; in this respect, the study of termites’ SINEs activity may provide an interesting framework to address the (co)evolution of mobile elements and the host genome.     

Extensive allele‐level remodeling of histone methylation modification in reciprocal F1 hybrids of rice subspecies

Epigenetic mechanisms play a major role in heterosis, partly as a result of the remodeling of epigenetic modifications in F₁ hybrids. Based on chromatin immunoprecipitation‐sequencing (ChIP‐Seq) analyses, we show that at the allele level extensive histone methylation remodeling occurred for a subset of genomic loci in reciprocal F₁ hybrids of Oryza sativa (rice) cultivars Nipponbare and 93‐11, representing the two subspecies japonica and indica. Globally, the allele modification‐altered loci in leaf or root of the reciprocal F₁ hybrids involved ?12–43% or more of the genomic regions carrying either of two typical histone methylation markers, H3K4me3 (>21 000 genomic regions) and H3K27me3 (>11 000 genomic regions). Nevertheless, at the total modification level, the majority (from ?43 to >90%) of the modification‐altered alleles lay within the range of parental additivity in the hybrids because of concerted alteration in opposite directions, consistent with an overall attenuation of allelic differences in the modifications. Importantly, of the genomic regions that did show non‐additivity in total modification level by either marker in the two tissues of hybrids, >80% manifested transgressivity, which involved genes enriched in specific functional categories. Extensive allele‐level alteration of H3K4me3 alone was positively correlated with genome‐wide changes in allele‐level gene expression, whereas at the total level, both H3K4me3 and H3K27me3 remodeling, although affecting just a small number of genes, contributes to the overall non‐additive gene expression to variable extents, depending on tissue/marker combinations. Our results emphasize the importance of allele‐level analysis in hybrids to assess the remodeling of epigenetic modifications and their relation to changes in gene expression.     

Sexually divergent DNA methylation patterns with hippocampal aging

DNA methylation is a central regulator of genome function, and altered methylation patterns are indicative of biological aging and mortality. Age‐related cellular, biochemical, and molecular changes in the hippocampus lead to cognitive impairments and greater vulnerability to neurodegenerative disease that varies between the sexes. The role of hippocampal epigenomic changes with aging in these processes is unknown as no genome‐wide analyses of age‐related methylation changes have considered the factor of sex in a controlled animal model. High‐depth, genome‐wide bisulfite sequencing of young (3 month) and old (24 month) male and female mouse hippocampus revealed that while total genomic methylation amounts did not change with aging, specific sites in CG and non‐CG (CH) contexts demonstrated age‐related increases or decreases in methylation that were predominantly sexually divergent. Differential methylation with age for both CG and CH sites was enriched in intergenic and intronic regions and under‐represented in promoters, CG islands, and specific enhancer regions in both sexes, suggesting that certain genomic elements are especially labile with aging, even if the exact genomic loci altered are predominantly sex‐specific. Lifelong sex differences in autosomal methylation at CG and CH sites were also observed. The lack of genome‐wide hypomethylation, sexually divergent aging response, and autosomal sex differences at CG sites was confirmed in human data. These data reveal sex as a previously unappreciated central factor of hippocampal epigenomic changes with aging. In total, these data demonstrate an intricate regulation of DNA methylation with aging by sex, cytosine context, genomic location, and methylation level.     

Epigenetic silencing of DACH1 induces the invasion and metastasis of gastric cancer by activating TGF‐β signalling

Gastric cancer (GC) is the fourth most common malignancy in males and the fifth most common malignancy in females worldwide. DACH1 is frequently methylated in hepatic and colorectal cancer. To further understand the regulation and mechanism of DACH1 in GC, eight GC cell lines, eight cases of normal gastric mucosa, 98 cases of primary GC and 50 cases of adjacent non‐tumour tissues were examined. Methylation‐specific PCR, western blot, transwell assay and xenograft mice were used in this study. Loss of DACH1 expression correlated with promoter region methylation in GC cells, and re‐expression was induced by 5‐Aza‐2′‐deoxyazacytidine. DACH1 is methylated in 63.3% (62/98) of primary GC and 38% (19/50) of adjacent non‐tumour tissues, while no methylation was found in normal gastric mucosa. Methylation of DACH1 correlated with reduced expression of DACH1 (P < 0.01), late tumour stage (stage III/IV) (P < 0.01) and lymph node metastasis (P < 0.05). DACH1 expression inhibited epithelial–mesenchymal transition and metastasis by inhibiting transforming growth factor (TGF)‐β signalling and suppressed GC cell proliferation through inducing G2/M phase arrest. The tumour size is smaller in DACH1‐expressed BGC823 cell xenograft mice than in unexpressed group (P < 0.01). Restoration of DACH1 expression also sensitized GC cells to docetaxel. These studies suggest that DACH1 is frequently methylated in human GC and expression of DACH1 was controlled by promoter region methylation. DACH1 suppresses GC proliferation, invasion and metastasis by inhibiting TGF‐β signalling pathways both in vitro and in vivo. Epigenetic silencing DACH1 may induce GC cells' resistance to docetaxel.     

RNA‐directed DNA methylation efficiency depends on trigger and target sequence identity

RNA‐directed DNA methylation (RdDM) in plants has been extensively studied, but the RNA molecules guiding the RdDM machinery to their targets are still to be characterized. It is unclear whether these molecules require full complementarity with their target. In this study, we have generated Nicotiana tabacum (Nt) plants carrying an infectious tomato apical stunt viroid (TASVd) transgene (Nt‐TASVd) and a non‐infectious potato spindle tuber viroid (PSTVd) transgene (Nt‐SB2). The two viroid sequences exhibit 81% sequence identity. Nt‐TASVd and Nt‐SB2 plants were genetically crossed. In the progeny plants (Nt‐SB2/TASVd), deep sequencing of small RNAs (sRNAs) showed that TASVd infection was associated with the accumulation of abundant small interfering RNAs (siRNAs) that mapped along the entire TASVd but only partially matched the SB2 transgene. TASVd siRNAs efficiently targeted SB2 RNA for degradation, but no transitivity was detectable. Bisulfite sequencing in the Nt‐SB2/TASVd plants revealed that the TASVd transgene was targeted for dense cis‐RdDM along its entire sequence. In the same plants, the SB2 transgene was targeted for trans‐RdDM. The SB2 methylation pattern, however, was weak and heterogeneous, pointing to a positive correlation between trigger–target sequence identity and RdDM efficiency. Importantly, trans‐RdDM on SB2 was also detected at sites where no homologous siRNAs were detected. Our data indicate that RdDM efficiency depends on the trigger–target sequence identity, and is not restricted to siRNA occupancy. These findings support recent data suggesting that RNAs with sizes longer than 24 nt (>24‐nt RNAs) trigger RdDM.     

Manipulation of cytosine methylation does not remove latitudinal clines in two invasive goldenrod species in Central Europe

Invasive species frequently differentiate phenotypically in novel environments within a few generations, often even with limited genetic variation. For the invasive plants Solidago canadensis and S. gigantea, we tested whether such differentiation might have occurred through heritable epigenetic changes in cytosine methylation. In a 2‐year common‐garden experiment, we grew plants from seeds collected along a latitudinal gradient in their non‐native Central European range to test for trait differentiation and whether differentiation disappeared when seeds were treated with the demethylation agent zebularine. Microsatellite markers revealed no population structure along the latitudinal gradient in S. canadensis, but three genetic clusters in S. gigantea. Solidago canadensis showed latitudinal clines in flowering phenology and growth. In S. gigantea, the number of clonal offspring decreased with latitude. Although zebularine had a significant effect on early growth, probably through effects on cytosine methylation, latitudinal clines remained (or even got stronger) in plants raised from seeds treated with zebularine. Thus, our experiment provides no evidence that epigenetic mechanisms by selective cytosine methylation contribute to the observed phenotypic differentiation in invasive goldenrods in Central Europe.     

Effects of single and multiple herbivory by host and non‐host caterpillars on the attractiveness of herbivore‐induced volatiles of sugarcane to the generalist parasitoid Cotesia flavipes

It is well known that parasitoids are attracted to volatiles emitted by host‐damaged plants; however, this tritrophic interaction may change if plants are attacked by more than one herbivore species. The larval parasitoid Cotesia flavipesCameron (Hymenoptera: Braconidae) has been used intensively in Brazil to control the sugarcane borer, Diatraea saccharalisFabricius (Lepidoptera: Pyralidae) in sugarcane crops, where Spodoptera frugiperda (JE Smith) (Lepidoptera: Noctuidae), a non‐stemborer lepidopteran, is also a pest. Here, we investigated the ability of C. flavipes to discriminate between an unsuitable host (S. frugiperda) and a suitable host (D. saccharalis) based on herbivore‐induced plant volatiles (HIPVs) emitted by sugarcane, and whether multiple herbivory (D. saccharalis feeding on stalk + S. frugiperda feeding on leaves) in sugarcane affected the attractiveness of HIPVs to C. flavipes. Olfactometer assays indicated that volatiles of host and non‐host‐damaged plants were attractive to C. flavipes. Even though host‐ and non‐host‐damaged plants emitted considerably different volatile blends, neither naïve nor experienced wasps discriminated suitable and unsuitable hosts by means of HIPVs emitted by sugarcane. With regard to multiple herbivory, wasps innately preferred the odor blend emitted by sugarcane upon non‐host + host herbivory over host‐only damaged plants. Multiple herbivory caused a suppression of some volatiles relative to non‐host‐damaged sugarcane that may have resulted from the unaltered levels of jasmonic acid in host‐damaged plants, or from reduced palatability of host‐damaged plants to S. frugiperda. In conclusion, our study showed that C. flavipes responds to a wide range of plant volatile blends, and does not discriminate host from non‐host and non‐stemborer caterpillars based on HIPVs emitted from sugarcane. Moreover, we showed that multiple herbivory by the sugarcane borer and fall armyworm increases the attractiveness of sugarcane plants to the parasitoids.     

Fine‐scale population epigenetic structure in relation to gastrointestinal parasite load in red grouse (Lagopus lagopus scotica)

Epigenetic modification of cytosine methylation states can be elicited by environmental stresses and may be a key process affecting phenotypic plasticity and adaptation. Parasites are potent stressors with profound physiological and ecological effects on their host, but there is little understanding in how parasites may influence host methylation states. Here, we estimate epigenetic diversity and differentiation among 21 populations of red grouse (Lagopus lagopus scotica) in north‐east Scotland and test for association of gastrointestinal parasite load (caecal nematode Trichostrongylus tenuis) with hepatic genome‐wide and locus‐specific methylation states. Following methylation‐sensitive AFLP (MSAP), 129 bands, representing 73 methylation‐susceptible and 56 nonmethylated epiloci, were scored across 234 individuals. The populations differed significantly in genome‐wide methylation levels and were also significantly epigenetically (F_SC = 0.0227; P < 0.001) and genetically (F_SC = 0.0058; P < 0.001) differentiated. Parasite load was not associated with either genome‐wide methylation levels or epigenetic differentiation. Instead, we found eight disproportionately differentiated epilocus‐specific methylation states (F_ST outliers) using bayescan software and significant positive and negative association of 35 methylation states with parasite load from bespoke generalized estimating equations (GEE), simple logistic regression (sam ) and Bayesian environmental analysis (bayenv 2). Following Sanger sequencing, genome mapping and geneontology (go ) annotation, some of these epiloci were linked to genes involved in regulation of cell cycle, signalling, metabolism, immune system and notably rRNA methylation, histone acetylation and small RNAs. These findings demonstrate an epigenetic signature of parasite load in populations of a wild bird and suggest intriguing physiological effects of parasite‐associated cytosine methylation.     

DNA methylation of retrotransposons,DNA transposons and genes in sugar beet (Beta vulgaris L.)

The methylation of cytosines shapes the epigenetic landscape of plant genomes, coordinates transgenerational epigenetic inheritance, represses the activity of transposable elements (TEs), affects gene expression and, hence, can influence the phenotype. Sugar beet (Beta vulgaris ssp. vulgaris), an important crop that accounts for 30% of worldwide sugar needs, has a relatively small genome size (758 Mbp) consisting of approximately 485 Mbp repetitive DNA (64%), in particular satellite DNA, retrotransposons and DNA transposons. Genome‐wide cytosine methylation in the sugar beet genome was studied in leaves and leaf‐derived callus with a focus on repetitive sequences, including retrotransposons and DNA transposons, the major groups of repetitive DNA sequences, and compared with gene methylation. Genes showed a specific methylation pattern for CG, CHG (H = A, C, and T) and CHH sites, whereas the TE pattern differed, depending on the TE class (class 1, retrotransposons and class 2, DNA transposons). Along genes and TEs, CG and CHG methylation was higher than that of adjacent genomic regions. In contrast to the relatively low CHH methylation in retrotransposons and genes, the level of CHH methylation in DNA transposons was strongly increased, pointing to a functional role of asymmetric methylation in DNA transposon silencing. Comparison of genome‐wide DNA methylation between sugar beet leaves and callus revealed a differential methylation upon tissue culture. Potential epialleles were hypomethylated (lower methylation) at CG and CHG sites in retrotransposons and genes and hypermethylated (higher methylation) at CHH sites in DNA transposons of callus when compared with leaves.     

High‐resolution melt analysis does not reveal mutagenic risk in sexed sperm and in vitro‐derived bovine embryos

The objectives of the present work were to verify whether simultaneous exposure to Hoechst 33342 and UV irradiation during sorting by flow cytometry may induce gene point mutations in bovine sperm and to assess whether the dye incorporated in the sperm may imply a mutagenic effect during the embryonic development. To this aim, high‐resolution melt analysis (HRMA) was used to discriminate variations of single nucleotides in sexed vs. non‐sexed control samples. Three batches of sorted and non‐sorted commercial semen of seven bulls (42 samples) were subjected to HRMA. A set of 139 genes located on all the chromosomes was selected, and 407 regions of the genome covering a total of 83 907 bases were analyzed. Thereafter, sperm of one sexed and one non‐sexed batch of each bull was used in in vitro fertilization, and the derived embryos were analyzed (n = 560). One hundred and thirty‐three regions of the bovine genome, located in 40 genes, were screened for a total coverage of 23 397 bases. The comparison between the frequencies of variations, with respect to the sequences deposited, observed in the sexed and non‐sexed sperm (843 vs. 770) and embryos (246 vs. 212) showed no significant differences (P > 0.05), as measured by chi‐square tests. It can be concluded that staining with Hoechst 33342 and exposure to UV during sorting does not lead to significant changes in the frequencies of variants in the commercial sexed semen and in embryos produced in vitro with the same treated sperm.     

Divergence of 3′ ends as a driver of short interspersed nuclear element (SINE) evolution in the Salicaceae

Short interspersed nuclear elements (SINEs) are small, non‐autonomous and heterogeneous retrotransposons that are widespread in plants. To explore the amplification dynamics and evolutionary history of SINE populations in representative deciduous tree species, we analyzed the genomes of the six following Salicaceae species: Populus deltoides, Populus euphratica, Populus tremula, Populus tremuloides, Populus trichocarpa, and Salix purpurea. We identified 11 Salicaceae SINE families (SaliS‐I to SaliS‐XI), comprising 27 077 full‐length copies. Most of these families harbor segmental similarities, providing evidence for SINE emergence by reshuffling or heterodimerization. We observed two SINE groups, differing in phylogenetic distribution pattern, similarity and 3′ end structure. These groups probably emerged during the ‘salicoid duplication’ (~65 million years ago) in the Salix–Populus progenitor and during the separation of the genus Salix (45–65 million years ago), respectively. In contrast to conserved 5′ start motifs across species and SINE families, the 3′ ends are highly variable in sequence and length. This extraordinary 3′‐end variability results from mutations in the poly(A) tail, which were fixed by subsequent amplificational bursts. We show that the dissemination of newly evolved 3′ ends is accomplished by a displacement of older motifs, leading to various 3′‐end subpopulations within the SaliS families.     

Effects of quantitative and qualitative differences in volatiles from host‐ and non‐host‐infested maize on the attraction of the larval parasitoid Cotesia kariyai

Cotesia kariyai Watanabe (Hymenoptera: Braconidae) is a specialist larval parasitoid of Mythimna separata Walker (Lepidoptera: Noctuidae). Cotesia kariyai wasps use herbivore‐induced plant volatiles (HIPVs) to locate hosts. However, complex natural habitats are full of volatiles released by both herbivorous host‐ and non‐host‐infested plants at various levels of intensity. Therefore, the presence of non‐hosts may affect parasitoid decisions while foraging. Here, the host‐finding efficiency of naive C. kariyai from HIPVs influenced by host‐ and non‐host‐infested maize [Zea mays L. (Poaceae)] plants was investigated with a four‐arm olfactometer. Ostrinia furnacalis Guenée (Lepidoptera: Crambidae) was selected as a non‐host species. One unit (1 U) of host‐ or non‐host‐infested plant was prepared by infesting a potted plant with five host or seven non‐host larvae. In two‐choice bioassays, host‐infested plants fed upon by different numbers of larvae, and various units of host‐ and non‐host‐infested plants (infestation units; 1 U, 2 U, and 3 U) were arranged to examine the effects of differences in volatile quantity and quality on the olfactory responses of C. kariyai with the assumption that volatile quantity and quality changes with differences in numbers of insects and plants. Cotesia kariyai was found to perceive quantitative differences in volatiles from host‐infested plants, preferring larger quantities of volatiles from larger numbers of larvae or plants. Also, the parasitoids discriminated between healthy plants, host‐infested plants, and non‐host‐infested plants by recognising volatiles released from those plants. Cotesia kariyai showed a reduced preference for host‐induced volatiles, when larger numbers of non‐host‐infested plants were present. Therefore, quantitative and qualitative differences in volatiles from host‐ and non‐host‐infested plants appear to affect the decision of C. kariyai during host‐habitat searching in multiple tritrophic systems.