Genus-Wide Characterization of Bumblebee Genomes Provides Insights into Their Evolution and Variation in Ecological and Behavioral Traits.

Bumblebees are a diverse group of globally important pollinators in natural ecosystems and for agricultural food production. With both eusocial and solitary life-cycle phases, and some social parasite species, they are especially interesting models to understand social evolution, behavior, and ecology. Reports of many species in decline point to pathogen transmission, habitat loss, pesticide usage, and global climate change, as interconnected causes. These threats to bumblebee diversity make our reliance on a handful of well-studied species for agricultural pollination particularly precarious. To broadly sample bumblebee genomic and phenotypic diversity, we de novo sequenced and assembled the genomes of 17 species, representing all 15 subgenera, producing the first genus-wide quantification of genetic and genomic variation potentially underlying key ecological and behavioral traits. The species phylogeny resolves subgenera relationships, whereas incomplete lineage sorting likely drives high levels of gene tree discordance. Five chromosome-level assemblies show a stable 18-chromosome karyotype, with major rearrangements creating 25 chromosomes in social parasites. Differential transposable element activity drives changes in genome sizes, with putative domestications of repetitive sequences influencing gene coding and regulatory potential. Dynamically evolving gene families and signatures of positive selection point to genus-wide variation in processes linked to foraging, diet and metabolism, immunity and detoxification, as well as adaptations for life at high altitudes. Our study reveals how bumblebee genes and genomes have evolved across the Bombus phylogeny and identifies variations potentially linked to key ecological and behavioral traits of these important pollinators.

Horizontal transfer of a novel Helentron in insects.

Helentrons represent a novel subtype of Helitrons. However, the evolutionary history of Helentrons in organisms is not clearly understood. In this study, we performed structure and autonomous partner analyses, which revealed that bm_455, a TE obtained from the Bombyx mori TE database, BmTEdb, was a member of Helentrons but not a long-terminal repeat (LTR) retrotransposon. Further analyses showed that bm_455 was also present in a wide range of insects including lepidopterans, coleopterans and hymenopterans using a homology-based search strategy. Several lines of evidence (high sequence identity, discontinuous distribution and lack of intense purifying selection) suggested that these elements could have been transferred into these species in part by horizontal transfers (HTs). Because Helentrons can capture host gene fragments, HTs of Helentrons might have a huge impact on their host genome evolution.

The Antioxidant Inhibition of Clove Effective Fraction on Lipid, Protein and Spectra Variation of LDL.

Oxidation of low density lipoprotein (LDL) has been considered as the critical factor which led to atherosclerosis (AS). Lipid and protein in LDL were oxidized to cause change of spectra during oxidation. Clove has been demonstrated to possess the strongest antioxidant capacity among 87 both medical and edible plants proclaimed by China. However, whether LDL oxidation is inhibited by clove? If so, whether it inhibits the oxidation of lipid and protein in LDL, and how does spectral characteristic of LDL change during oxidation when clove was added into LDL. Currently, these questions are still unclear. Therefore, the inhibition efficiency of the effective fraction from clove (EFC) on oxidation of LDL was studied by UV-Visible spectra and fluorescence spectra. The results indicated that EFC might effectively delay propagation of conjugated diene (CD) during LDL oxidation and postpone reaching maximum of its content. Likewise, it might inhibit cholesterol degradation, tryptophan (Trp) fluorescence quenching, lysine (Lys) residues from being modified, and lipofuscins from being generated in peroxidatic reaction among lipid and protein. Besides, EFC also might affect variation of UV-Visible spectra of LDL during oxidation. This study provides reference for future research and development of clove functional food inhibiting AS.

FishTEDB 2.0: an update fish transposable element (TE) database with new functions to facilitate TE research.

We launched the initial version of FishTEDB in 2018, which aimed to establish an open-source, user-friendly, data-rich transposable element (TE) database. Over the past 5 years, FishTEDB 1.0 has gained approximately 10 000 users, accumulating more than 450 000 interactions. With the unveiling of extensive fish genome data and the increasing emphasis on TE research, FishTEDB needs to extend the richness of data and functions. To achieve the above goals, we introduced 33 new fish species to FishTEDB 2.0, encompassing a wide array of fish belonging to 48 orders. To make the updated database more functional, we added a genome browser to visualize the positional relationship between TEs and genes and the estimated TE insertion time in different species. In conclusion, we released a new version of the fish TE database, FishTEDB 2.0, designed to assist researchers in the future study of TE functions and promote the progress of biological theories related to TEs. Database URL: https://www.fishtedb.com/.

A telomere-to-telomere reference genome of ficus (Ficus hispida) provides new insights into sex determination.

A high-quality reference genome is indispensable for resolving biologically essential traits. Ficus hispida is a dioecious plant. A complete Ficus reference genome will be crucial for understanding their sex evolution and important biological characteristics, such as aerial roots, mutualistic symbiosis with ficus-wasps, and fruiting from old stems. Here, we generated a telomere-to-telomere (T2T) genome for F. hispida using PacBio HiFi and Oxford Nanopore Ultra-long sequencing technologies. The genome contiguity and completeness has shown improvement compared with the previously released genome, with the annotation of six centromeres and 28 telomeres. We have refined our previously reported 2-Mb male-specific region into a 7.2-Mb genomic region containing 51 newly predicted genes and candidate sex-determination genes AG2 and AG3. Many of these genes showed extremely low expression, likely attributed to hypermethylation in the gene body and promoter regions. Gene regulatory networks (GRNs) revealed that AG2 and AG3 are related to the regulation of stamen development in male flowers, while the AG1 gene is responsible for regulating female flowers' defense responses and secondary metabolite processes. Comparative analysis of GRNs showed that the NAC, WRKY, and MYB transcription factor families dominate the female GRN, whereas the MADS and MYB transcription factor families are prevalent in the male GRN.

Simultaneous extraction of crocin and geniposide from gardenia fruits (Gardenia jasminoides Ellis) by probe-type ultrasound-assisted natural deep eutectic solvents and their inhibition effects on low density lipoprotein oxidation.

The simultaneous extraction of crocin and geniposide from gardenia fruits (Gardenia jasminoides Ellis) was performed by integrating natural deep eutectic solvents (NADES) and ultrasound-assisted extraction (UAE). Among the eight kinds of NADES screened, choline chloride-1,2-propylene glycol was the most suitable extractant. The probe-type ultrasound-assisted NADES extraction system (pr-UAE-NADES) demonstrated higher extraction efficiency compared with plate-type ultrasound-assisted NADES extraction system (pl-UAE-NADES). Orthogonal experimental design and a modified multi-index synthetic weighted scoring method were adopted to optimize pr-UAE-NADES extraction process. The optimal extraction conditions that had a maximum synthetic weighted score of 29.46 were determined to be 25 °C for extraction temperature, 600 W for ultrasonic power, 20 min for extraction time, and 25% (w/w) for water content in NADES, leading to the maximum yields (7.39 ± 0.20 mg/g and 57.99 ± 0.91 mg/g, respectively) of crocin and geniposide. Thirty-three compounds including iridoids, carotenoids, phenolic acids, flavonoids, and triterpenes in the NADES extract were identified by LC-Q-TOF-MS2 coupled with a feature-based molecular networking workflow. The kinetics evaluation of the conjugated dienes generation on Cu2+-induced low density lipoprotein (LDL) oxidation via the four-parameter logistic regression model showed that crocin increased the lag time of LDL oxidation in a concentration-dependent manner (15 µg/mL, 30 µg/mL, 45 µg/mL) by 12.66%, 35.44%, and 73.42%, respectively. The quantitative determination for fluorescence properties alteration of the apolipoprotein B-100 exhibited that crocin effectively inhibited the fluorescence quenching of tryptophan residues and the modification of lysine residues caused by reactive aldehydes and malondialdehydes. The pr-UAE-NADES showed significant efficiency toward the simultaneous extraction of crocin and geniposide from gardenia fruits. And this study demonstrates the potential utility of gardenia fruits in developing anti-atherogenic functional food.

Maternal dominance contributes to subgenome differentiation in allopolyploid fishes.

Teleost fishes, which are the largest and most diverse group of living vertebrates, have a rich history of ancient and recent polyploidy. Previous studies of allotetraploid common carp and goldfish (cyprinids) reported a dominant subgenome, which is more expressed and exhibits biased gene retention. However, the underlying mechanisms contributing to observed 'subgenome dominance' remains poorly understood. Here we report high-quality genomes of twenty-one cyprinids to investigate the origin and subsequent subgenome evolution patterns following three independent allopolyploidy events. We identify the closest extant relatives of the diploid progenitor species, investigate genetic and epigenetic differences among subgenomes, and conclude that observed subgenome dominance patterns are likely due to a combination of maternal dominance and transposable element densities in each polyploid. These findings provide an important foundation to understanding subgenome dominance patterns observed in teleost fishes, and ultimately the role of polyploidy in contributing to evolutionary innovations.

[Computational approaches for identification and classification of transposable elements in eukaryotic genomes].

Repetitive sequences (repeats) represent a significant fraction of the eukaryotic genomes and can be divided into tandem repeats, segmental duplications, and interspersed repeats on the basis of their sequence characteristics and how they are formed. Most interspersed repeats are derived from transposable elements (TEs). Eukaryotic TEs have been subdivided into two major classes according to the intermediate they use to move. The transposition and amplification of TEs have a great impact on the evolution of genes and the stability of genomes. However, identification and classification of TEs are complex and difficult due to the fact that their structure and classification are complex and diverse compared with those of other types of repeats. Here, we briefly introduced the function and classification of TEs, and summarized three different steps for identification, classification and annotation of TEs in eukaryotic genomes: (1) assembly of a repeat library, (2) repeat correction and classification, and (3) genome annotation. The existing computational approaches for each step were summarized and the advantages and disadvantages of the approaches were also highlighted in this review. To accurately identify, classify, and annotate the TEs in eukaryotic genomes requires combined methods. This review provides useful information for biologists who are not familiar with these approaches to find their way through the forest of programs.

Horizontal transfer and evolution of transposable elements in vertebrates.

Horizontal transfer of transposable elements (HTT) is an important process shaping eukaryote genomes, yet very few studies have quantified this phenomenon on a large scale or have evaluated the selective constraints acting on transposable elements (TEs) during vertical and horizontal transmission. Here we screen 307 vertebrate genomes and infer a minimum of 975 independent HTT events between lineages that diverged more than 120 million years ago. HTT distribution greatly differs from null expectations, with 93.7% of these transfers involving ray-finned fishes and less than 3% involving mammals and birds. HTT incurs purifying selection (conserved protein evolution) on all TEs, confirming that producing functional transposition proteins is required for a TE to invade new genomes. In the absence of HTT, DNA transposons appear to evolve neutrally within genomes, unlike most retrotransposons, which evolve under purifying selection. This selection regime indicates that proteins of most retrotransposon families tend to process their own encoding RNA (cis-preference), which helps retrotransposons to persist within host lineages over long time periods.

High-quality genome assembly and transcriptome of Ancherythroculter nigrocauda, an endemic Chinese cyprinid species.

Ancherythroculter nigrocauda is a cyprinid fish endemic of the upper reaches of the Yangtze River in China, where it is an important aquaculture and commercial species. It is also a threatened species as a result of overfishing, dam construction and water pollution. In this study, a chromosome-level genome assembly of A. nigrocauda is reported and built using PacBio sequencing and the Hi-C technology. The 1.04-Gb sequenced genome of A. nigrocauda contained 2,403 contigs, with an N50 length of 3.12 Mb. Then, 1,297 contigs, which represented 54.0% of all contigs and 97.2% of the whole content of the genome nucleotide base, were assembled into 24 chromosomes. Combined with transcriptome data from 10 tissues, 27,042 (78.5%) genes were functionally annotated out of 34,414 predicted protein-coding genes. Interestingly, high expression of many positively selected genes and expanded gene families in the brain suggested that these genes might play important roles in brain development in A. nigrocauda. Finally, we found tissue-specific expression of 10,732 genes. Functional analyses showed that they were mainly composed of genes related to (a) environmental information processing, (b) the circulatory system, and (c) development, suggesting they might be important for adaptation to different environments and for development of A. nigrocauda. The high-quality genome obtained in this study not only provides a valuable genomic resource for future studies of A. nigrocauda populations and conservation, but is also an important resource for further functional genomics studies of fishes.

Successful Invasions of Short Internally Deleted Elements (SIDEs) and Its Partner CR1 in Lepidoptera Insects.

Although DNA transposons often generated internal deleted derivatives such as miniature inverted-repeat transposable elements, short internally deleted elements (SIDEs) derived from nonlong terminal-repeat retrotransposons are rare. Here, we found a novel SIDE, named Persaeus, that originated from the chicken repeat 1 (CR1) retrotransposon Zenon and it has been found widespread in Lepidoptera insects. Our findings suggested that Persaeus and the partner Zenon have experienced a transposition burst in their host genomes and the copy number of Persaeus and Zenon in assayed genomes are significantly correlated. Accordingly, the activity though age analysis indicated that the replication wave of Persaeus coincided with that of Zenon. Phylogenetic analyses suggested that Persaeus may have evolved at least four times independently, and that it has been vertically transferred into its host genomes. Together, our results provide new insights into the evolution dynamics of SIDEs and its partner non-LTRs.

Evolutionary dynamics of transposable elements during silkworm domestication.

Although there are some documented examples on population dynamics of transposable elements (TEs) in model organisms, the evolutionary dynamics of TEs in domesticated species has not been systematically investigated. The objective of this study is to understand population dynamics of TEs during silkworm domestication. In this work, using transposon-display we examined the polymorphism of seven TE families [they represent about 59% of silkworm (Bombyx mori) total TE content] in four domesticated silkworm populations and one wild silkworm population. Maximum likelihood (ML) was used to estimate selection pressure. Population differentiation and structure were performed by using AMOVA analysis and program DISTRUCT, respectively. The results of transposon-display showed that significant differentiation occurred between the domesticated silkworm and wild silkworm. These TEs have experienced expansions and fixation in the domesticated silkworm but not in wild silkworm. Furthermore, the ML results indicated that purifying selection of TEs in the domesticated silkworm were significantly weaker than that in the wild silkworm. Interestingly, an adaptation insertion induced by BmMITE-2 was found, and this insertion can reduce the polymorphism of the flanking regions of its neighboring COQ7 gene. Our results suggested that TEs expanded and were fixed in the domesticated silkworm might result from demographic effects and artificial selection during domestication. We concluded that the data presented in this study have general implication in animal and crop improvements as well as in domestication of new species.

Genome-wide identification and evolution of TC1/Mariner in the silkworm (Bombyx mori) genome.

TC1/Mariner transposons belong to class II transposable elements (TEs) that use DNA-mediated "cut and paste" mechanism to transpose, and they have been identified in almost all organisms. Although silkworm (Bombyx mori) has a large amount of TC1/Mariner elements, the genome wide information of this superfamily in the silkworm is unknown. In this study, we have identified 2670 TC1/Mariner (Bmmar) elements in the silkworm genome. All the TEs were classified into 22 families by means of fgclust, a tool of repetitive sequence classification, seven of which was first reported in this study. Phylogenetic and structure analyses based on the catalytic domain (DDxD/E) of transposase sequences indicated that all members of TC1/Mariner were grouped into five subgroups: Mariner, Tc1, maT, DD40D and DD41D/E. Of these five subgroups, maT rather than Mariner possessed most members of TC1/Mariner (51.23%) in the silkworm genome. In particular, phylogenetic analysis and structure analysis revealed that Bmmar15 (DD40D) formed a new basal subgroup of TC1/Mariner element in insects, which was referred to as bmori. Furthermore, we concluded that DD40D appeared to intermediate between mariner and Tc1. Finally, we estimated the insertion time for each copy of TC1/Mariner in the silkworm and found that most of members were dramatically amplified during a period from 0 to 1 mya. Moreover, the detailed functional data analysis showed that Bmmar1, Bmmar6 and Bmmar9 had EST evidence and intact transposases. These implied that TC1/Mariner might have potential transpositional activity. In conclusion, this study provides some new insights into the landscape, origin and evolution of TC1/Mariner in the insect genomes.

Unexpected invasion of miniature inverted-repeat transposable elements in viral genomes.

BACKGROUND: Transposable elements (TEs) are common and often present with high copy numbers in cellular genomes. Unlike in cellular organisms, TEs were previously thought to be either rare or absent in viruses. Almost all reported TEs display only one or two copies per viral genome. In addition, the discovery of pandoraviruses with genomes up to 2.5-Mb emphasizes the need for biologists to rethink the fundamental nature of the relationship between viruses and cellular life. RESULTS: Herein, we performed the first comprehensive analysis of miniature inverted-repeat transposable elements (MITEs) in the 5170 viral genomes for which sequences are currently available. Four hundred and fifty one copies of ten miniature inverted-repeat transposable elements (MITEs) were found and each MITE had reached relatively large copy numbers (some up to 90) in viruses. Eight MITEs belonging to two DNA superfamilies (hobo/Activator/Tam3 and Chapaev-Mirage-CACTA) were for the first time identified in viruses, further expanding the organismal range of these two superfamilies. TEs may play important roles in shaping the evolution of pandoravirus genomes, which were here found to be very rich in MITEs. We also show that putative autonomous partners of seven MITEs are present in the genomes of viral hosts, suggesting that viruses may borrow the transpositional machinery of their cellular hosts' autonomous elements to spread MITEs and colonize their own genomes. The presence of seven similar MITEs in viral hosts, suggesting horizontal transfers (HTs) as the major mechanism for MITEs propagation. CONCLUSIONS: Our discovery highlights that TEs contribute to shape genome evolution of pandoraviruses. We concluded that as for cellular organisms, TEs are part of the pandoraviruses' diverse mobilome.

iMITEdb: the genome-wide landscape of miniature inverted-repeat transposable elements in insects.

Miniature inverted-repeat transposable elements (MITEs) have attracted much attention due to their widespread occurrence and high copy numbers in eukaryotic genomes. However, the systematic knowledge about MITEs in insects and other animals is still lacking. In this study, we identified 6012 MITE families from 98 insect species genomes. Comparison of these MITEs with known MITEs in the NCBI non-redundant database and Repbase showed that 5701(∼95%) of 6012 MITE families are novel. The abundance of MITEs varies drastically among different insect species, and significantly correlates with genome size. In general, larger genomes contain more MITEs than small genomes. Furthermore, all identified MITEs were included in a newly constructed database (iMITEdb) (http://gene.cqu.edu.cn/iMITEdb/), which has functions such as browse, search, BLAST and download. Overall, our results not only provide insight on insect MITEs but will also improve assembly and annotation of insect genomes. More importantly, the results presented in this study will promote studies of MITEs function, evolution and application in insects. DATABASE URL: http://gene.cqu.edu.cn/iMITEdb/.

TRT, a Vertebrate and Protozoan Tc1-Like Transposon: Current Activity and Horizontal Transfer.

We report a Danio rerio transposon named DrTRT, for D. rerio Transposon Related to Tc1 The complete sequence of the DrTRT transposon is 1,563 base pairs (bp) in length, and its transposase putatively encodes a 338-amino acid protein that harbors a DD37E motif in its catalytic domain. We present evidence based on searches of publicly available genomes that TRT elements commonly occur in vertebrates and protozoa. Phylogenetic and functional domain comparisons confirm that TRT constitutes a new subfamily within the Tc1 family. Hallmark features of having no premature termination codons within the transposase, the presence of all expected functional domains, and its occurrence in the bony fish transcriptome suggest that TRT might have current or recent activity in these species. Further analysis showed that the activity of TRT elements in these species might have arisen about between 4 and 19 Ma. Interestingly, our results also implied that the widespread distribution of TRT among fishes, frog, and snakes is the result of multiple independent HT events, probably from bony fishes to snakes or frog. Finally, the mechanisms underlying horizontal transfer of TRT elements are discussed.

Repeated horizontal transfers of four DNA transposons in invertebrates and bats.

BACKGROUND: Horizontal transfer (HT) of transposable elements (TEs) into a new genome is considered as an important force to drive genome variation and biological innovation. However, most of the HT of DNA transposons previously described occurred between closely related species or insects. RESULTS: In this study, we carried out a detailed analysis of four DNA transposons, which were found in the first sequenced twisted-wing parasite, Mengenilla moldrzyki. Through the homology-based strategy, these transposons were also identified in other insects, freshwater planarian, hydrozoans, and bats. The phylogenetic distribution of these transposons was discontinuous, and they showed extremely high sequence identities (>87%) over their entire length in spite of their hosts diverging more than 300 million years ago (Mya). Additionally, phylogenies and comparisons of transposons versus orthologous gene identities demonstrated that these transposons have transferred into their hosts by independent HTs. CONCLUSIONS: Here, we provided the first documented example of HT of CACTA transposons, which have been so far extensively studied in plants. Our results demonstrated that bats had continuously acquired new DNA elements via HT. This implies that predation on a large quantity of insects might increase bat exposure to HT. In addition, parasite-host interaction might facilitate exchanging of their genetic materials.

The diversification of PHIS transposon superfamily in eukaryotes.

BACKGROUND: PHIS transposon superfamily belongs to DNA transposons and includes PIF/Harbinger, ISL2EU, and Spy transposon groups. These three groups have similar DDE domain-containing transposases; however, their coding capacity, species distribution, and target site duplications (TSDs) are significantly different. RESULTS: In this study, we systematically identified and analyzed PHIS transposons in 836 sequenced eukaryotic genomes using transposase homology search and structure approach. In total, 380 PHIS families were identified in 112 genomes and 168 of 380 families were firstly reported in this study. Besides previous identified PIF/Harbinger, ISL2EU, and Spy groups, three new types (called Pangu, NuwaI, and NuwaII) of PHIS superfamily were identified; each has its own distinctive characteristics, especially in TSDs. Pangu and NuwaII transposons are characterized by 5'-ANT-3' and 5'-C|TNA|G-3' TSDs, respectively. Both transposons are widely distributed in plants, fungi, and animals; the NuwaI transposons are characterized by 5'-CWG-3' TSDs and mainly distributed in animals. CONCLUSIONS: Here, in total, 380 PHIS families were identified in eukaryotes. Among these 380 families, 168 were firstly reported in this study. Furthermore, three new types of PHIS superfamily were identified. Our results not only enrich the transposon diversity but also have extensive significance for improving genome sequence assembly and annotation of higher organisms.

Recurrent horizontal transfers of Chapaev transposons in diverse invertebrate and vertebrate animals.

Horizontal transfer (HT) of a transposable element (TE) into a new genome is regarded as an important force to drive genome variation and biological innovation. In addition, HT also plays an important role in the persistence of TEs in eukaryotic genomes. Here, we provide the first documented example for the repeated HT of three families of Chapaev transposons in a wide range of animal species, including mammals, reptiles, jawed fishes, lampreys, insects, and in an insect bracovirus. Multiple alignments of the Chapaev transposons identified in these species revealed extremely high levels of nucleotide sequence identity (79-99%), which are inconsistent with vertical evolution given the deep divergence time separating these host species. Rather, the discontinuous distribution amongst species and lack of purifying selection acting on these transposons strongly suggest that they were independently and horizontally transferred into these species lineages. The detection of Chapaev transposons in an insect bracovirus indicated that these viruses might act as a possible vector for the horizontal spread of Chapaev transposons. One of the Chapaev families was also shared by lampreys and some of their common hosts (such as sturgeon and paddlefish), which suggested that parasite-host interaction might facilitate HTs.