Natural selection maintains the transcribed LTR retrotransposons in Nosema bombycis

Eight intact LTR retrotransposons (Nbrl-Nbr8) have been previously characterized from the genome of Nosema bombycis, a eu-karyotic parasite with a compact and reduced genome. Here we describe six novel transcribed Nbr elements (Nbr9-Nbrl4) identified through either cDNA library or RT-PCR. Like previously determined ones, all of them belong to the Ty3/Gypsy superfamily. Retro-transposon diversity and incomplete domains with insertions {Nbrll), deletions (Nbrll) and in-frame stop codons in coding regions (Nbr9) were detected, suggesting that both defective and loss events of LTR retrotransposon have happened in N. bombycis genome.Analysis of selection showed that strong purifying selection acts on all elements except Nbrll. This implies that selective pressure keeps both these Nbrs and their functions in genome. Interestingly, Nbrll is under positive selection and some positively selected codons were identified, indicating that new functionality might have evolved in the Nbrll retrotransposon. Unlike other transposable elements, Nbrll has integrated into a conserved syntenic block and probably resulted in the inversion of both flanking regions. This demonstrates that transposable element is an important factor for the reshuffling and evolution of their host genomes, and may be maintained under natural selection.     

Natural selection maintains the transcribed LTR retrotransposons in Nosema bombycis

Eight intact LTR retrotransposons(Nbr1?Nbr8)have been previously characterized from the genome of Nosema bombycis,a eu-karyotic parasite with a compact and reduced genome.Here we describe six novel transcribed Nbr elements(Nbr9?Nbr14)identified through either cDNA library or RT-PCR.Like previously determined ones,all of them belong to the Ty3/Gypsy superfamily.Retrotransposon diversity and incomplete domains with insertions(Nbr12),deletions(Nbr11)and in-frame stop codons in coding regions(Nbr9)were detected...     

Codon-based tests of positive selection,branch lengths,and the evolution of mammalian immune system genes

Using basic probability theory, we show that there is a substantial likelihood that even in the presence of strong purifying selection, there will be a number of codons in which the number of synonymous nucleotide substitutions per site (d (S)) exceeds the number of non-synonymous nucleotide substitutions per site (d (N)). In an empirical study, we examined the numbers of synonymous (b (S)) and non-synonymous substitutions (b (N)) along branches of the phylogenies of 69 single-copy orthologous genes from seven species of mammals. A pattern of b (N) > b (S) was most commonly seen in the shortest branches of the tree and was associated with a high coefficient of variation in both b (N) and b (S), suggesting that high stochastic error in b (N) and b (S) on short branches, rather than positive Darwinian selection, is the explanation of most cases where b (N) is greater than b (S) on a given branch. The branch-site method of Zhang et al. (Zhang, Nielsen, Yang, Mol Biol Evol, 22:2472-2479, 2005) identified 117 codons on 35 branches as "positively selected," but a majority of these codons lacked synonymous substitutions, while in the others, synonymous and non-synonymous differences per site occurred in approximately equal frequencies. Thus, it was impossible to rule out the hypothesis that chance variation in the pattern of mutation across sites, rather than positive selection, accounted for the observed pattern. Our results showed that b (N)/b (S) was consistently elevated in immune system genes, but neither the search for branches with b (N) > b (S) nor the branch-site method revealed this trend.     

Two families of non-LTR retrotransposons, Syrinx and Daphne, from the Darwinulid ostracod, Darwinula stevensoni

Two novel families of non-LTR retrotransposons, named Syrinx and Daphne, were cloned and characterized in a putative ancient asexual ostracod Darwinula stevensoni. Phylogenetic analysis reveals that Daphne is the founding member of a novel clade of non-LTR retroelements, which also contains retrotransposon families from the sea urchin and the silkworm and forms a sister clade to L2-like elements. The Syrinx family of non-LTR retrotransposons exhibits evidence of relatively recent activity, manifested in high levels of sequence similarity between individual copies and a three- to ten-fold excess of synonymous substitutions, which is indicative of purifying selection. The Daphne family may have very few copies with intact open reading frames, and exhibits neutral within-family ratio of non-synonymous to synonymous substitutions. It can additionally be characterized by formation of inverted truncated head-to-head structures. All of these features make recent activity less likely than in the Syrinx family. Our results are discussed in light of the evolutionary consequences of long-term asexuality in general and in D. stevensoni in particular.     

Wake up of transposable elements following Drosophila simulans worldwide colonization.

Transposable elements (TEs) make up around 10%-15% of the Drosophila melanogaster genome, but its sibling species Drosophila simulans carries only one third as many such repeat sequences. We do not, however, have an overall view of copy numbers of the various classes of TEs (long terminal repeat [LTR] retrotransposons, non-LTR retrotransposons, and transposons) in genomes of natural populations of both species. We analyzed 34 elements in individuals from various natural populations of these species. We show that D. melanogaster has higher average chromosomal insertion site numbers per genome than D. simulans for all TEs except five. The LTR retrotransposons gypsy, ZAM, and 1731 and the transposon bari-1 present similar low copy numbers in both species. The transposon hobo has a large number of insertion sites, with significantly more sites in D. simulans. High variation between populations in number of insertion sites of some elements of D. simulans suggests that these elements can invade the genome of the entire species starting from a local population. We propose that TEs in the D. simulans genome are being awakened and amplified as they had been a long time ago in D. melanogaster.     

A computational study of the dynamics of LTR retrotransposons in the Populus trichocarpa genome

Retrotransposons are an ubiquitous component of plant genomes, especially abundant in species with large genomes. Populus trichocarpa has a relatively small genome, which was entirely sequenced; however, studies focused on poplar retrotransposons dynamics are rare. With the aim to study the retrotransposon component of the poplar genome, we have scanned the complete genome sequence searching full-length long-terminal repeat (LTR) retrotransposons, i.e., characterised by two long terminal repeats at the 5′ and 3′ ends. A computational approach based on detection of conserved structural features, on building multiple alignments, and on similarity searches was used to identify 1,479 putative full-length LTR retrotransposons. Ty1-copia elements were more numerous than Ty3-gypsy. However, many LTR retroelements were not assigned to any superfamily because lacking of diagnostic features and non-autonomous. LTR retrotransposon remnants were by far more numerous than full-length elements, indicating that during the evolution of poplar, large amplification of these elements was followed by DNA loss. Within superfamilies, Ty3-gypsy families are made of more members than Ty1-copia ones. Retrotransposition occurred with increasing frequency following the separation of Populus sections, with different waves of retrotransposition activity between Ty3-gypsy and Ty1-copia elements. Recently inserted elements appear more frequently expressed than older ones. Finally, different levels of activity of retrotransposons were observed according to their position and their density in the linkage groups. On the whole, the results support the view of retrotransposons as a community of different organisms in the genome, whose activity (both retrotransposition and DNA loss) has heavily impacted and probably continues to impact poplar genome structure and size.     

Genetic diversity of near genome-wide hepatitis C virus sequences during chronic infection: evidence for protein structural conservation over time

Infection with hepatitis C virus (HCV) is one of the leading causes of chronic hepatitis, liver cirrhosis and end-stage liver disease worldwide. The genetics of HCV infection in humans and the disease course of chronic hepatitis C are both remarkably variable. Although the response to interferon treatment is largely dependent on HCV genotypes, whether or not a relationship exists between HCV genome variability and clinical course of hepatitis C disease still remains unknown. To more thoroughly understand HCV genome evolution over time in association with disease course, near genome-wide HCV genomes present in 9 chronically infected participants over 83 total study years were sequenced. Overall, within HCV genomes, the number of synonymous substitutions per synonymous site (d(S)) significantly exceeded the number of non-synonymous substitutions per site (d(N)). Although both d(S) and d(N) significantly increased with duration of chronic infection, there was a highly significant decrease in d(N)/d(S) ratio in HCV genomes over time. These results indicate that purifying selection acted to conserve viral protein structure despite persistence of high level of nucleotide mutagenesis inherent to HCV replication. Based on liver biopsy fibrosis scores, HCV genomes from participants with advanced fibrosis had significantly greater d(S) values and lower d(N)/d(S) ratios compared to participants with mild liver disease. Over time, viral genomes from participants with mild disease had significantly greater annual changes in d(N), along with higher d(N)/d(S) ratios, compared to participants with advanced fibrosis. Yearly amino acid variations in the HCV p7, NS2, NS3 and NS5B genes were all significantly lower in participants with severe versus mild disease, suggesting possible pathogenic importance of protein structural conservation for these viral gene products.     

Structural and evolutionary analyses of the Ty3/gypsy group of LTR retrotransposons in the genome of Anopheles gambiae

The recent availability of the genome of Anopheles gambiae offers an extraordinary opportunity for comparative studies of the diversity of transposable elements (TEs) and their evolutionary dynamics between two related species, taking advantage of the existing information from Drosophila melanogaster. To this goal, we screened the genome of A. gambiae for elements belonging to the Ty3/gypsy group of long-terminal repeat (LTR) retrotransposons. The A. gambiae genome displays a rich diversity of LTR retrotransposons, clearly greater than D. melanogaster. We have characterized in detail 63 families, belonging to five of the nine main lineages of the Ty3/gypsy group. The Mag lineage is the most diverse and abundant, with more than 30 families. In sharp contrast with this finding, a single family belonging to this lineage has been found in D. melanogaster, here reported for the first time in the literature, most probably consisting of old inactive elements. The CsRn1 lineage is also abundant in A. gambiae but almost absent from D. melanogaster. Conversely, the Osvaldo lineage has been detected in Drosophila but not in Anopheles. Comparison of structural characteristics of different families led to the identification of several lineage-specific features such as the primer-binding site (PBS), the gag-pol translational recoding signal (TRS), which is extraordinarily diverse within the Ty3/gypsy retrotransposons of A. gambiae, or the presence/absence of specific amino acid motifs. Interestingly, some of these characteristics, although in general well conserved within lineages, may have evolved independently in particular branches of the phylogenetic tree. We also show evidence of recent activity for around 75% of the families. Nevertheless, almost all families contain a high proportion of degenerate members and solitary LTRs (solo LTRs), indicative of a lower turnover rate of retrotransposons belonging to the Ty3/gypsy group in A. gambiae than in D. melanogaster. Finally, we have detected significant overrepresentations of insertions on the X chromosome versus autosomes and of putatively active insertions on euchromatin versus heterochromatin.     

The ingi and RIME non-LTR retrotransposons are not randomly distributed in the genome of Trypanosoma brucei

The ingi (long and autonomous) and RIME (short and nonautonomous) non--long-terminal repeat retrotransposons are the most abundant mobile elements characterized to date in the genome of the African trypanosome Trypanosoma brucei. These retrotransposons were thought to be randomly distributed, but a detailed and comprehensive analysis of their genomic distribution had not been performed until now. To address this question, we analyzed the ingi/RIME sequences and flanking sequences from the ongoing T. brucei genome sequencing project (TREU927/4 strain). Among the 81 ingi/RIME elements analyzed, 60% are complete, and 7% of the ingi elements (approximately 15 copies per haploid genome) appear to encode for their own transposition. The size of the direct repeat flanking the ingi/RIME retrotransposons is conserved (i.e., 12-bp), and a strong 11-bp consensus pattern precedes the 5'-direct repeat. The presence of a consensus pattern upstream of the retroelements was confirmed by the analysis of the base occurrence in 294 GSS containing 5'-adjacent ingi/RIME sequences. The conserved sequence is present upstream of ingis and RIMEs, suggesting that ingi-encoded enzymatic activities are used for retrotransposition of RIMEs, which are short nonautonomous retroelements. In conclusion, the ingi and RIME retroelements are not randomly distributed in the genome of T. brucei and are preceded by a conserved sequence, which may be the recognition site of the ingi-encoded endonuclease.     

Whole-Genome Sequencing of Theileria parva Strains Provides Insight into Parasite Migration and Diversification in the African Continent

The disease caused by the apicomplexan protozoan parasite Theileria parva, known as East Coast fever or Corridor disease, is one of the most serious cattle diseases in Eastern, Central, and Southern Africa. We performed whole-genome sequencing of nine T. parva strains, including one of the vaccine strains (Kiambu 5), field isolates from Zambia, Uganda, Tanzania, or Rwanda, and two buffalo-derived strains. Comparison with the reference Muguga genome sequence revealed 34 814–121 545 single nucleotide polymorphisms (SNPs) that were more abundant in buffalo-derived strains. High-resolution phylogenetic trees were constructed with selected informative SNPs that allowed the investigation of possible complex recombination events among ancestors of the extant strains. We further analysed the dN/dS ratio (non-synonymous substitutions per non-synonymous site divided by synonymous substitutions per synonymous site) for 4011 coding genes to estimate potential selective pressure. Genes under possible positive selection were identified that may, in turn, assist in the identification of immunogenic proteins or vaccine candidates. This study elucidated the phylogeny of T. parva strains based on genome-wide SNPs analysis with prediction of possible past recombination events, providing insight into the migration, diversification, and evolution of this parasite species in the African continent.     

Whole-genome and chromosome evolution associated with host adaptation and speciation of the wheat pathogen Mycosphaerella graminicola

The fungus Mycosphaerella graminicola has been a pathogen of wheat since host domestication 10,000-12,000 years ago in the Fertile Crescent. The wheat-infecting lineage emerged from closely related Mycosphaerella pathogens infecting wild grasses. We use a comparative genomics approach to assess how the process of host specialization affected the genome structure of M. graminicola since divergence from the closest known progenitor species named M. graminicola S1. The genome of S1 was obtained by Illumina sequencing resulting in a 35 Mb draft genome sequence of 32X. Assembled contigs were aligned to the previously sequenced M. graminicola genome. The alignment covered >90% of the non-repetitive portion of the M. graminicola genome with an average divergence of 7%. The sequenced M. graminicola strain is known to harbor thirteen essential chromosomes plus eight dispensable chromosomes. We found evidence that structural rearrangements significantly affected the dispensable chromosomes while the essential chromosomes were syntenic. At the nucleotide level, the essential and dispensable chromosomes have evolved differently. The average synonymous substitution rate in dispensable chromosomes is considerably lower than in essential chromosomes, whereas the average non-synonymous substitution rate is three times higher. Differences in molecular evolution can be related to different transmission and recombination patterns, as well as to differences in effective population sizes of essential and dispensable chromosomes. In order to identify genes potentially involved in host specialization or speciation, we calculated ratios of synonymous and non-synonymous substitution rates in the >9,500 aligned protein coding genes. The genes are generally under strong purifying selection. We identified 43 candidate genes showing evidence of positive selection, one encoding a potential pathogen effector protein. We conclude that divergence of these pathogens was accompanied by structural rearrangements in the small dispensable chromosomes, while footprints of positive selection were present in only a small number of protein coding genes.     

TEnest: automated chronological annotation and visualization of nested plant transposable elements

Organisms with a high density of transposable elements (TEs) exhibit nesting, with subsequent repeats found inside previously inserted elements. Nesting splits the sequence structure of TEs and makes annotation of repetitive areas challenging. We present TEnest, a repeat identification and display tool made specifically for highly repetitive genomes. TEnest identifies repetitive sequences and reconstructs separated sections to provide full-length repeats and, for long-terminal repeat (LTR) retrotransposons, calculates age since insertion based on LTR divergence. TEnest provides a chronological insertion display to give an accurate visual representation of TE integration history showing timeline, location, and families of each TE identified, thus creating a framework from which evolutionary comparisons can be made among various regions of the genome. A database of repeats has been developed for maize (Zea mays), rice (Oryza sativa), wheat (Triticum aestivum), and barley (Hordeum vulgare) to illustrate the potential of TEnest software. All currently finished maize bacterial artificial chromosomes totaling 29.3 Mb were analyzed with TEnest to provide a characterization of the repeat insertions. Sixty-seven percent of the maize genome was found to be made up of TEs; of these, 95% are LTR retrotransposons. The rate of solo LTR formation is shown to be dissimilar across retrotransposon families. Phylogenetic analysis of TE families reveals specific events of extreme TE proliferation, which may explain the high quantities of certain TE families found throughout the maize genome. The TEnest software package is available for use on PlantGDB under the tools section (http://www.plantgdb.org/prj/TE_nest/TE_nest.html); the source code is available from (http://wiselab.org).     

Characterization and phylogenetic relationships among microsporidia infecting silkworm, Bombyx mori, using inter simple sequence repeat (ISSR) and small subunit rRNA (SSU-rRNA) sequence analysis.

This study is the first report on the genetic characterization and relationships among different microsporidia infecting the silkworm, Bombyx mori, using inter simple sequence repeat PCR (ISSR-PCR) analysis. Six different microsporidians were distinguished through molecular DNA typing using ISSR-PCR. Thus, ISSR-PCR analysis can be a powerful tool to detect polymorphisms and identify microsporidians, which are difficult to study with microscopy because of their extremely small size. Of the 100 ISSR primers tested, only 28 primers had reproducibility and high polymorphism (93%). A total of 24 ISSR primers produced 55 unique genetic markers, which could be used to differentiate the microsporidians from each other. Among the 28 SSRs tested, the most abundant were (CA)n, (GA)n, and (GT)n repeats. The degree of band sharing was used to evaluate genetic similarity between different microsporidian isolates and to construct a phylogenetic tree using Jaccard's similarity coefficient. The results indicate that the DNA profiles based on ISSR markers can be used as diagnostic tools to identify different microsporidia with considerable accuracy. In addition, the small subunit ribosomal RNA (SSU-rRNA) sequence gene was amplified, cloned, and sequenced from each of the 6 microsporidian isolates. These sequences were compared with 20 other microsporidian SSU-rRNA sequences to develop a phylogenetic tree for the microsporidia isolated from the silkworms. This method was found to be useful in establishing the phylogenetic relationships among the different microsporidians isolated from silkworms. Of the 6 microsporidian isolates, NIK-1s revealed an SSU-rRNA gene sequence similar to Nosema bombycis, indicating that NIK-1s is similar to N. bombycis; the remaining 5 isolates, which differed from each other and from N. bombycis, were considered to be different variants belonging to the species N. bombycis.     

Phylogenetic relationships of three new microsporidian isolates from the silkworm, Bombyx mori

The pathogenicity, mode of transmission, tissue specificity of infection and the small subunit rRNA (SSU-rRNA) gene sequences of the three new microsporidian isolates from the silkworm Bombyx mori were studied. Out of the three, NIK-2r revealed life cycle features and SSU-rRNA gene sequence similar to Nosema bombycis, suggesting that it is N. bombycis. The other two, NIK-4m and NIK-3h, differed from each other as well as from N. bombycis. NIK-4m was highly pathogenic and did not show any vertical transmission, in accordance with the apparent lack of gonadal infection, whereas NIK-3h was less pathogenic and vertical transmission was not detected but could not be excluded. Phylogenetic analysis based on SSU-rRNA gene sequence placed NIK-3h and NIK-4m in a distinct clade that included almost all the Vairimorpha species and Nosema species that infect lepidopteran and non-lepidopteran hosts, while NIK-2r was included in a clade containing almost all the Nosema isolates that infect only lepidopteran hosts. Thus, we have presented molecular evidence that one of the three isolates is in fact the type species N. bombycis, while the other two isolates are Vairimorpha spp. There was distinct separation of microsporidian isolates infecting only lepidopteran hosts and those infecting lepidopteran and non-lepidopteran hosts, reflecting possible co-evolution of hosts and microsporidian isolates.     

家蚕微孢子虫中小RNAs的全基因组鉴定与分析

【目的】家蚕Bombyx mori微粒子病是蚕业生产上的毁灭性病害,家蚕微孢子虫Nosema bombycis是该病的病原,可经卵垂直传播和经口水平传播。为了探索家蚕微孢子虫中对重复元件的抵御以及对基因转录调控的潜在方式,本研究拟在基因组水平上对该物种的小RNAs进行全面系统的分析,鉴定与转座子相关的小RNAs和潜在的miRNAs。【方法】从感染家蚕微孢子虫的家蚕中肠中提取总RNA,分离小片段RNA并反转录后,进行Solexa高通量测序。通过生物信息学方法对小RNAs进行分类及功能注释,鉴定起源于家蚕微孢子虫不同类型转座子的小RNAs,并对潜在的miRNA进行预测分析。【结果】家蚕微孢子虫小RNAs的长度主要是24和25 nt,其中大部分序列表现出5′末端的尿嘧啶偏好性。家蚕微孢子虫中存在丰富的与转座子相关联的小RNAs,并且与转座子标准序列匹配的反义小RNAs明显多于正义小RNAs。同时,鉴定获得了31个候选miRNAs,部分为Nosema属的其他孢子虫中所共有,暗示其在微孢子虫基因组进化上具有保守性。【结论】首次鉴定到家蚕微孢子虫的转座子相关性小RNAs,暗示小RNAs在家蚕微孢子虫基因组对转座子防御过程中起到作用,31个潜在的miRNAs为家蚕微孢子虫miRNAs的功能验证提供了后续靶标。     

SearchDOGS Bacteria,Software That Provides Automated Identification of Potentially Missed Genes in Annotated Bacterial Genomes

We report the development of SearchDOGS Bacteria, software to automatically detect missing genes in annotated bacterial genomes by combining BLAST searches with comparative genomics. Having successfully applied the approach to yeast genomes, we redeveloped SearchDOGS to function as a standalone, downloadable package, requiring only a set of GenBank annotation files as input. The software automatically generates a homology structure using reciprocal BLAST and a synteny-based method; this is followed by a scan of the entire genome of each species for unannotated genes. Results are provided in a HTML interface, providing coordinates, BLAST results, syntenic location, omega values (Ka/Ks, where Ks is the number of synonymous substitutions per synonymous site and Ka is the number of nonsynonymous substitutions per nonsynonymous site) for protein conservation estimates, and other information for each candidate gene. Using SearchDOGS Bacteria, we identified 155 gene candidates in the Shigella boydii sb227 genome, including 56 candidates of length < 60 codons. SearchDOGS Bacteria has two major advantages over currently available annotation software. First, it outperforms current methods in terms of sensitivity and is highly effective at identifying small or highly diverged genes. Second, as a freely downloadable package, it can be used with unpublished or confidential data.