TE-Tracker: systematic identification of transposition events through whole-genome resequencing

Background

Transposable elements (TEs) are DNA sequences that are able to move from their location in the genome by cutting or copying themselves to another locus. As such, they are increasingly recognized as impacting all aspects of genome function. With the dramatic reduction in cost of DNA sequencing, it is now possible to resequence whole genomes in order to systematically characterize novel TE mobilization in a particular individual. However, this task is made difficult by the inherently repetitive nature of TE sequences, which in some eukaryotes compose over half of the genome sequence. Currently, only a few software tools dedicated to the detection of TE mobilization using next-generation-sequencing are described in the literature. They often target specific TEs for which annotation is available, and are only able to identify families of closely related TEs, rather than individual elements.

Results

We present TE-Tracker, a general and accurate computational method for the de-novo detection of germ line TE mobilization from re-sequenced genomes, as well as the identification of both their source and destination sequences. We compare our method with the two classes of existing software: specialized TE-detection tools and generic structural variant (SV) detection tools. We show that TE-Tracker, while working independently of any prior annotation, bridges the gap between these two approaches in terms of detection power. Indeed, its positive predictive value (PPV) is comparable to that of dedicated TE software while its sensitivity is typical of a generic SV detection tool. TE-Tracker demonstrates the benefit of adopting an annotation-independent, de novo approach for the detection of TE mobilization events. We use TE-Tracker to provide a comprehensive view of transposition events induced by loss of DNA methylation in Arabidopsis. TE-Tracker is freely available at http://www.genoscope.cns.fr/TE-Tracker.

Conclusions

We show that TE-Tracker accurately detects both the source and destination of novel transposition events in re-sequenced genomes. Moreover, TE-Tracker is able to detect all potential donor sequences for a given insertion, and can identify the correct one among them. Furthermore, TE-Tracker produces significantly fewer false positives than common SV detection programs, thus greatly facilitating the detection and analysis of TE mobilization events.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0377-z) contains supplementary material, which is available to authorized users.     

Expression analysis of heat shock protein genes during Aeromonas hydrophila infection in rohu,Labeo rohita,with special reference to molecular characterization of Grp78

Heat shock protein (Hsp) genes are stress-related genes that activate the host immune system during infection. Hsp genes expression in fish, studied during bacterial infections, is mostly confined to Hsp70 and Hsp90. The present study is an expression analysis of seven Hsp genes: Apg2, Hsp90, Hsp70, glucose-regulated protein 78 (Grp78), heat shock cognate 70 (Hsc70), Grp75, and Hsp30 during Aeromonas hydrophila infection in rohu, Labeo rohita. Forty-eight rohu juveniles were challenged with 3 × 10⁷ cfu bacteria per 20 g of body weight intraperitoneally. The expression of these genes was studied in infected liver, anterior kidney, and spleen tissues of rohu at different time periods: 0, 1, 3, 6, 12, 24, 48, 72 h, 7, and 15 days post-infection by qPCR. The Hsp gene modulation was greater in liver as compared to spleen and kidney tissues. Induced expression of Apg2, Hsp90, Grp78, Grp75, and Hsc70 was noticed during peak periods (3 to 24 h post-challenge) of bacterial infectivity. Hsp70 was found to be down-regulated during the process of infection. In contrast to the other six genes, Hsp30 showed a variable expression pattern in all three tissues. Grp78 was found to be the most highly (1,587-fold) expressed gene in liver at 12 h post-challenge. Further, molecular characterization of Grp78 revealed it to be a highly conserved Hsp gene, essential not only during infection but also during early developmental stages of rohu, and its expression was noticed in all organs studied except in gill tissues, which indicated its potential immune regulatory role during infection process.     

Mutagenesis of Trichoderma reesei endoglucanase I: impact of expression host on activity and stability at elevated temperatures

Background

Trichoderma reesei is a key cellulase source for economically saccharifying cellulosic biomass for the production of biofuels. Lignocellulose hydrolysis at temperatures above the optimum temperature of T. reesei cellulases (~50°C) could provide many significant advantages, including reduced viscosity at high-solids loadings, lower risk of microbial contamination during saccharification, greater compatibility with high-temperature biomass pretreatment, and faster rates of hydrolysis. These potential advantages motivate efforts to engineer T. reesei cellulases that can hydrolyze lignocellulose at temperatures ranging from 60–70°C.

Results

A B-factor guided approach for improving thermostability was used to engineer variants of endoglucanase I (Cel7B) from T. reesei (TrEGI) that are able to hydrolyze cellulosic substrates more rapidly than the recombinant wild-type TrEGI at temperatures ranging from 50–70°C. When expressed in T. reesei, TrEGI variant G230A/D113S/D115T (G230A/D113S/D115T Tr_TrEGI) had a higher apparent melting temperature (3°C increase in T_m) and improved half-life at 60°C (t_1/2 = 161 hr) than the recombinant (T. reesei host) wild-type TrEGI (t_1/2 = 74 hr at 60°C, Tr_TrEGI). Furthermore, G230A/D113S/D115T Tr_TrEGI showed 2-fold improved activity compared to Tr_TrEGI at 65°C on solid cellulosic substrates, and was as efficient in hydrolyzing cellulose at 60°C as Tr_TrEGI was at 50°C. The activities and stabilities of the recombinant TrEGI enzymes followed similar trends but differed significantly in magnitude depending on the expression host (Escherichia coli cell-free, Saccharomyces cerevisiae, Neurospora crassa, or T. reesei). Compared to N.crassa-expressed TrEGI, S. cerevisiae-expressed TrEGI showed inferior activity and stability, which was attributed to the lack of cyclization of the N-terminal glutamine in Sc_TrEGI and not to differences in glycosylation. N-terminal pyroglutamate formation in TrEGI expressed in S. cerevisiae was found to be essential in elevating its activity and stability to levels similar to the T. reesei or N. crassa-expressed enzyme, highlighting the importance of this ubiquitous modification in GH7 enzymes.

Conclusion

Structure-guided evolution of T. reesei EGI was used to engineer enzymes with increased thermal stability and activity on solid cellulosic substrates. Production of TrEGI enzymes in four hosts highlighted the impact of the expression host and the role of N-terminal pyroglutamate formation on the activity and stability of TrEGI enzymes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12896-015-0118-z) contains supplementary material, which is available to authorized users.     

Genome-wide characterization of long intergenic non-coding RNAs (lincRNAs) provides new insight into viral diseases in honey bees Apis cerana and Apis mellifera

Background

Long non-coding RNAs (lncRNAs) are a class of RNAs that do not encode proteins. Recently, lncRNAs have gained special attention for their roles in various biological process and diseases.

Results

In an attempt to identify long intergenic non-coding RNAs (lincRNAs) and their possible involvement in honey bee development and diseases, we analyzed RNA-seq datasets generated from Asian honey bee (Apis cerana) and western honey bee (Apis mellifera). We identified 2470 lincRNAs with an average length of 1011 bp from A. cerana and 1514 lincRNAs with an average length of 790 bp in A. mellifera. Comparative analysis revealed that 5 % of the total lincRNAs derived from both species are unique in each species. Our comparative digital gene expression analysis revealed a high degree of tissue-specific expression among the seven major tissues of honey bee, different from mRNA expression patterns. A total of 863 (57 %) and 464 (18 %) lincRNAs showed tissue-dependent expression in A. mellifera and A. cerana, respectively, most preferentially in ovary and fat body tissues. Importantly, we identified 11 lincRNAs that are specifically regulated upon viral infection in honey bees, and 10 of them appear to play roles during infection with various viruses.

Conclusions

This study provides the first comprehensive set of lincRNAs for honey bees and opens the door to discover lincRNAs associated with biological and hormone signaling pathways as well as various diseases of honey bee.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1868-7) contains supplementary material, which is available to authorized users.     

Genomic landscape and evolutionary dynamics of mariner transposable elements within the Drosophila genus

Background

The mariner family of transposable elements is one of the most widespread in the Metazoa. It is subdivided into several subfamilies that do not mirror the phylogeny of these species, suggesting an ancient diversification. Previous hybridization and PCR studies allowed a partial survey of mariner diversity in the Metazoa. In this work, we used a comparative genomics approach to access the genus-wide diversity and evolution of mariner transposable elements in twenty Drosophila sequenced genomes.

Results

We identified 36 different mariner lineages belonging to six distinct subfamilies, including a subfamily not described previously. Wide variation in lineage abundance and copy number were observed among species and among mariner lineages, suggesting continuous turn-over. Most mariner lineages are inactive and contain a high proportion of damaged copies. We showed that, in addition to substitutions that rapidly inactivate copies, internal deletion is a major mechanism contributing to element decay and the generation of non-autonomous sublineages. Hence, 23% of copies correspond to several Miniature Inverted-repeat Transposable Elements (MITE) sublineages, the first ever described in Drosophila for mariner. In the most successful MITEs, internal deletion is often associated with internal rearrangement, which sheds light on the process of MITE origin. The estimation of the transposition rates over time revealed that all lineages followed a similar progression consisting of a rapid amplification burst followed by a rapid decrease in transposition. We detected some instances of multiple or ongoing transposition bursts. Different amplification times were observed for mariner lineages shared by different species, a finding best explained by either horizontal transmission or a reactivation process. Different lineages within one species have also amplified at different times, corresponding to successive invasions. Finally, we detected a preference for insertion into short TA-rich regions, which appears to be specific to some subfamilies.

Conclusions

This analysis is the first comprehensive survey of this family of transposable elements at a genus scale. It provides precise measures of the different evolutionary processes that were hypothesized previously for this family based on PCR data analysis. mariner lineages were observed at almost all “life cycle” stages: recent amplification, subsequent decay and potential (re)-invasion or invasion of genomes.

Electronic supplementary material

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

Endogenous siRNAs and piRNAs derived from transposable elements and genes in the malaria vector mosquito Anopheles gambiae

Background

The siRNA and piRNA pathways have been shown in insects to be essential for regulation of gene expression and defence against exogenous and endogenous genetic elements (viruses and transposable elements). The vast majority of endogenous small RNAs produced by the siRNA and piRNA pathways originate from repetitive or transposable elements (TE). In D. melanogaster, TE-derived endogenous siRNAs and piRNAs are involved in genome surveillance and maintenance of genome integrity. In the medically relevant malaria mosquito Anopheles gambiae TEs constitute 12-16% of the genome size. Genetic variations induced by TE activities are known to shape the genome landscape and to alter the fitness in An. gambiae.

Results

Here, using bioinformatics approaches we analyzed the small RNA data sets from 6 libraries formally reported in a previous study and examined the expression of the mixed germline/somatic siRNAs and piRNAs produced in adult An. gambiae females. We characterized a large population of TE-derived endogenous siRNAs and piRNAs, which constitutes 56-60% of the total siRNA and piRNA reads in the analysed libraries. Moreover, we identified a number of protein coding genes producing gene-specific siRNAs and piRNAs that were generally expressed at much lower levels than the TE-associated small RNAs. Detailed sequence analysis revealed that An. gambiae piRNAs were produced by both “ping-pong” dependent (TE-associated piRNAs) and independent mechanisms (genic piRNAs). Similarly to D. melanogaster, more than 90% of the detected piRNAs were produced from TE-associated clusters in An. gambiae. We also found that biotic stress as blood feeding and infection with Plasmodium parasite, the etiological agent of malaria, modulated the expression levels of the endogenous siRNAs and piRNAs in An. gambiae.

Conclusions

We identified a large and diverse set of the endogenously derived siRNAs and piRNAs that share common and distinct aspects of small RNA expression across insect species, and inferred their impact on TE and gene activity in An. gambiae. The TE-specific small RNAs produced by both the siRNA and piRNA pathways represent an important aspect of genome stability and genetic variation, which might have a strong impact on the evolution of the genome and vector competence in the malaria mosquitoes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1436-1) contains supplementary material, which is available to authorized users.     

Dietary polyunsaturated fatty acids and the Pro12Ala polymorphisms of PPARG regulate serum lipids through divergent pathways: a randomized crossover clinical trial

Human and animal studies suggest an interaction between the Pro12Ala polymorphism of PPARG and dietary fat. In this randomized crossover clinical trial, we investigated whether subjects with the Pro12Pro and Ala12Ala genotypes of PPARG respond differently to a diet supplemented with high saturated (SAFA) or polyunsaturated fatty acid (PUFA).We recruited non-diabetic men from a population-based METSIM study (including 10,197 men) to obtain men with the Ala12Ala and the Pro12Pro genotypes matched for age and body mass index. Seventeen men with the Pro12Pro genotype and 14 with the Ala12Ala genotype were randomized to both a PUFA diet and a SAFA diet for 8 weeks in a crossover setting. Serum lipids and adipose tissue mRNA expression were measured during the diet intervention. At baseline, subjects with the Ala12Ala genotype had higher levels of HDL cholesterol and lower levels of LDL cholesterol, total triglycerides, and apolipoprotein B compared to those subjects with the Pro12Pro genotype (P < 0.05, FDR < 0.1). The Ala12Ala genotype also associated with higher mRNA expression of PPARG2, LPIN1, and SREBP-1c compared to participants with the Pro12Pro genotype (FDR < 0.001). On the other hand, PUFA diet resulted in lower levels of fasting glucose, total cholesterol, total triglycerides, and apolipoprotein B (P < 0.05, FDR < 0.1) but did not affect PPARG2 mRNA expression in adipose tissue. We conclude that individuals with the Pro12Pro genotype, with higher triglyceride levels at baseline, are more likely to benefit from the PUFA diet. However, the beneficial effects of dietary PUFA and the Ala12Ala genotype of PPARG on serum lipids are mediated through divergent mechanisms.

Electronic supplementary material

The online version of this article (doi:10.1007/s12263-015-0493-z) contains supplementary material, which is available to authorized users.     

Characterization of Arabidopsis sterol glycosyltransferase TTG15/UGT80B1 role during freeze and heat stress

Sterol glycosyltransferases regulate the properties of sterols by catalyzing the transfer of carbohydrate molecules to the sterol moiety for the synthesis of steryl glycosides and acyl steryl glycosides. We have analyzed the functional role of TTG15/UGT80B1 gene of Arabidopsis thaliana in freeze/thaw and heat shock stress using T-DNA insertional sgt knockout mutants. Quantitative study of spatial as well as temporal gene expression showed tissue-specific and dynamic expression patterns throughout the growth stages. Comparative responses of Col-0, TTG15/UGT80B1 knockout mutant and p35S:TTG15/UGT80B1 restored lines were analyzed under heat and freeze stress conditions. Heat tolerance was determined by survival of plants at 42°C for 3 h, MDA analysis and chlorophyll fluorescence image (CFI) analysis. Freezing tolerance was determined by survival of the plants at -1°C temperature in non-acclimatized (NA) and cold acclimatized (CA) conditions and also by CFI analysis, which revealed that, p35S:TTG15/UGT80B1 restored plants were more adapted to freeze stress than TTG15/UGT80B1 knockout mutant under CA condition. HPLC analysis of the plants showed reduced sterol glycoside in mutant seedlings as compared to other genotypes. Following CA condition, both β-sitosterol and sitosterol glycoside quantity was more in Col-0 and p35S:TTG15/UGT80B1 restored lines, whereas it was significantly less in TTG15/UGT80B1 knockout mutants. From these results, it may be concluded that due to low content of free sterols and sterol glycosides, the physiology of mutant plants was more affected during both, the chilling and heat stress.     

Rates of transposition in Escherichia coli

The evolutionary role of transposable elements (TEs) is still highly controversial. Two key parameters, the transposition rate (u and w, for replicative and non-replicative transposition) and the excision rate (e) are fundamental to understanding their evolution and maintenance in populations. We have estimated u, w and e for six families of TEs (including eight members: IS1, IS2, IS3, IS4, IS5, IS30, IS150 and IS186) in Escherichia coli, using a mutation accumulation (MA) experiment. In this experiment, mutations accumulate essentially at the rate at which they appear, during a period of 80 500 (1610 generations × 50 lines) generations, and spontaneous transposition events can be detected. This differs from other experiments in which insertions accumulated under strong selective pressure or over a limited genomic target. We therefore provide new estimates for the spontaneous rates of transposition and excision in E. coli. We observed 25 transposition and three excision events in 50 MA lines, leading to overall rate estimates of u ∼ 1.15 × 10^–5, w ∼ 4 × 10⁻⁸ and e ∼ 1.08 × 10⁻⁶ (per element, per generation). Furthermore, extensive variation between elements was found, consistent with previous knowledge of the mechanisms and regulation of transposition for the different elements.