Comparative Analysis of Genomic Repeat Content in Gomphocerine Grasshoppers Reveals Expansion of Satellite DNA and Helitrons in Species with Unusually Large Genomes

Transcriptional dynamics and tissue-specific expression patterns of transposable elements in orthopteran insects

Transposable elements (TEs) are prevalent in the genomes of orthopteran insects, contributing significantly to their genome evolution and diversity. In light of the existing gap in our understanding of TEs transcript dynamics in orthopteran insects, we recognize the critical need to undertake comprehensive analyses in this area. Therefore, we have decided to delve into the characterization of TE transcripts, their abundance profiles, and the formation of chimeric transcripts using RNA-seq data and genome assemblies. The transcript analysis of TEs across various species revealed significant differences in TE abundance and expression patterns. In particular, Schistocerca americana exhibited twice the number of transcripts within the genus Schistocerca compared to the average of other species, while Gryllus bimaculatus displayed the lowest number of transcripts. Despite this, all Schistocerca species shared similar fractions of TE transcripts at the clade level, with DNA transposons (45%) being the most abundant, followed by LINE (19%) and LTR elements (18%). Interestingly, Acrida cinerea displayed different TE abundance patterns compared to Schistocerca species, particularly with an increased proportion of LTR transcripts, accounting for 31% of the total transcripts. Further analysis revealed tissue-specific transcriptional activity of TE clades, with notable differences between male and female specimens. In Gryllus bimaculatus, TEs were highly transcribed across ovaries and gut tissues in females compared to male testes and gut. Conversely, Gastrimargus marmoratus displayed higher TE transcription in male tissues compared to females, indicating species-specific expression patterns. A similar pattern has been observed in Acrida cinerea, except in female gonads, where 4618 TEs were transcribed compared to 3757 in male gonads. Despite these variations, no correlation was found between genome size and TE transcript abundance. Additionally, highly conserved TEs were involved in the formation of chimeric transcripts, indicating potential regulatory roles in gene expression. The expression quantification analysis of chimeric TEs and genes revealed tissue-specific expression patterns, and TEs do not control the overall expression of all genes except some, suggesting regulatory roles of TEs in gene expression. Overall, our study underscores tissue-specific variations in TE expression and transcript abundance among different species. Additionally, our findings suggest the involvement of highly conserved TEs in the formation of chimeric transcripts across different species.

Repetitive DNAs: the 'invisible' regulators of insect adaptation and speciation

Like other eukaryotes, insect genomes contain a large portion of repetitive sequences, particularly transposable elements and satellite DNAs. This review highlights key studies on repetitive DNAs and examines their structural, functional, and evolutionary impact on insect genomes. Repetitive sequences promote genetic diversification through mutations and large-scale rearrangements, playing a crucial role in shaping genomic architecture, aiding organismal adaptation, and driving speciation. We also explore the influence of repeats in genome size variation and species incompatibilities, along with their contribution to adaptive phenotypes and gene regulation. Studying repetitive DNA in insects not only provides insights into basic genomic features but also offers valuable information for conservation strategies, pest control, and advancements in genetics, ecology, and evolutionary biology.

Expansion of three types of transposon superfamilies within 25 Mya lead to large genome size of a rice insect pest

The brown planthoppers (BPH, Nilaparvata lugens), white backed planthopper (WBPH, Sogatella furcifera) and small brown planthopper (SBPH, Laodelphax striatellus) are widely distributed rice insect pests, causing huge annual yield loss of rice production. Though these three planthoppers belong to the same family, Delphacidae of Hemiptera, their genome sizes (GS) are very different, ranging from 541 to 1088 Mb. To uncover the main factors driving GS changes of three planthoppers, we first estimated the GS of their ancestor Fulgoroidea, to be 794.33 Mb, indicating GS expansion in BPH but contraction in SBPH and WBPH. Next, we identified repetitive sequences and compared the TE landscapes, showed that three types of transposon superfamilies, hAT, Tc1-Mariner and Gypsy, expanded within 25 Mya in BPH. In addition, BPH kept ancient TEs of Fulgoroidea dated back to 175 Mya, while SBPH and WBPH have lost most of these ancient TEs. Here, we present evidence that the gain of recently expanded TEs driving the GS expansion and loss of ancient TEs leading to the GS contraction, providing new insights into the mechanism of GS variation.

Evolutionary dynamics of repetitive elements and their relationship with genome size in Acrididae

It is widely accepted that repetitive elements (REs) represent the primary mechanism driving genome size variation across eukaryotes. The observed genome sizes and REs of 59 species within the Acrididae were obtained and characterized. The genome sizes observed ranged from 6.60 pg to 19.35 pg, while the proportion of REs varied from 57.92 % to 83.58 %. The primary contributors were identified as LTR (2.34 % ∼ 20.98 %) and LINEs (6.70 % ∼ 16.33 %). The results of ancestral reconstruction indicated that the proportion of REs in ancestral nodes was 69.53 %, which suggests that they have undergone extensive genome expansion or contraction. A significant positive correlation was identified between the proportion of REs and genome size. Transposable elements were found to account for approximately 41 % of the observed variation in genome size. Moreover, the LTR was identified as the most significant RE type in relation to genome size expansion within the Acrididae.

Orthopteran Neo-Sex Chromosomes Reveal Dynamics of Recombination Suppression and Evolution of Supergenes

The early evolution of sex chromosomes has remained obscure for more than a century. The Vandiemenella viatica species group of morabine grasshoppers is highly suited for studying the early stages of sex chromosome divergence and degeneration of the Y chromosome. This stems from the fact that neo-XY sex chromosomes have independently evolved multiple times by X-autosome fusions with different autosomes. Here, we generated new chromosome-level assemblies for two chromosomal races representing karyotypes with and without neo-sex chromosomes (P24XY and P24X0), and sequence data of a third chromosomal race with a different neo-XY chromosome system (P25XY). Interestingly, these two neo-XY chromosomal races are formed by different X-autosome fusions (involving chr1 and chrB, respectively), and we found that both neo-Y chromosomes have partly ceased to recombine with their neo-X counterpart. We show that the neo-XY chromosomes have diverged through accumulation of SNPs and structural mutations, and that many neo-Y-linked genes have degenerated since recombination ceased. However, the non-recombining regions of neo-Y chromosomes host non-degenerated genes crucial for sex determination, such as sex-lethal and transformer, alongside genes associated with spermatogenesis, fertility, and reproduction, illustrating their integrative role as a masculinizing supergene. Contrary to expectations, the neo-Y chromosomes showed (slightly) lower density of transposable elements (TEs) compared to other genomic regions. The study reveals the unique dynamics of young sex chromosomes, with evolution of recombination suppression and pronounced decay of (some) neo-sex chromosome genes, and provides a compelling case illustrating how chromosomal fusions and post-fusion mutational processes contribute to the evolution of supergenes.

The evolutionary dynamics of genome sizes and repetitive elements in Ensifera (Insecta: Orthoptera)

BACKGROUND

In evolutionary biology, identifying and quantifying inter-lineage genome size variation and elucidating the underlying causes of that variation have long been goals. Repetitive elements (REs) have been proposed and confirmed as being among the most important contributors to genome size variation. However, the evolutionary implications of genome size variation and RE dynamics are not well understood.

RESULTS

A total of 35 Ensifera insects were collected from different areas in China, including nine species of crickets and 26 species of katydids. The genome sizes of seven species were then determined using flow cytometry. The RepeatExplorer2 pipeline was employed to retrieve the repeated sequences for each species, based on low-coverage (0.1 X) high-throughput Illumina unassembled short reads. The genome sizes of the 35 Ensifera insects exhibited a considerable degree of variation, ranging from 1.00 to 18.34 pg. This variation was more than 18-fold. Similarly, the RE abundances exhibited considerable variation, ranging from 13.66 to 61.16%. In addition, the Tettigonioidea had larger genomes and contained significantly more REs than did the Grylloidea genomes. Analysis of the correlation between RE abundance and the genome size of 35 Ensifera insects revealed that the abundance of REs, transposable elements (TEs), long terminal repeats (LTRs), and long interspersed nuclear elements (LINEs) are significantly correlated with genome size. Notably, there is an inflection point in this correlation, where species with increasingly large genomes (e.g., > 5-10 pg) have repeats that contribute less to genome expansion than expected. Furthermore, this study revealed contrasting evolutionary directions between the Tettigonioidea and Grylloidea clades in terms of the expansion of REs. Tettigonioidea species exhibit a gradual increase in ancestral genome size and RE abundance as they diverge, while Grylloidea species experience sustained genome contraction.

CONCLUSIONS

This study reveals extensive variation in genome size and RE abundance in Ensifera insects, with distinct evolutionary patterns across two major groups, Tettigonioidea and Grylloidea. This provides valuable insights into the variation in genome size and RE abundance in Ensifera insects, offering a comprehensive understanding of their evolutionary history.

Repetitive element expansions contribute to genome size gigantism in Pamphagidae: A comparative study (Orthoptera, Acridoidea)

Pamphagidae is a family of Acridoidea that inhabits the desert steppes of Eurasia and Africa. This study employed flow cytometry to estimate the genome size of eight species in the Pamphagidae. The results indicate that the genome size of the eight species ranged from 13.88 pg to 14.66 pg, with an average of 14.26 pg. This is the largest average genome size recorded for the Orthoptera families, as well as for the entire Insecta. Furthermore, the study explored the role of repetitive sequences in the genome, including their evolutionary dynamics and activity, using low-coverage next-generation sequencing data. The genome is composed of 14 different types of repetitive sequences, which collectively make up between 59.9% and 68.17% of the total genome. The Pamphagidae family displays high levels of transposable element (TE) activity, with the number of TEs increasing and accumulating since the family's emergence. The study found that the types of repetitive sequences contributing to the TE outburst events are similar across species. Additionally, the study identified unique repetitive elements for each species. The differences in repetitive sequences among the eight Pamphagidae species correspond to their phylogenetic relationships. The study sheds new light on genome gigantism in the Pamphagidae and provides insight into the correlation between genome size and repetitive sequences within the family.

Evolutionary Dynamics of Satellite DNA Repeats across the Tettigoniidae Family: Insights from Genomic Analysis

Satellite DNA repeats are repetitive DNA sequences found in eukaryotic genomes, typically consisting of short DNA motifs repeated in tandem arrays. Despite the vast body of literature on satellite DNA repeats in other taxa, investigations specifically targeting Tettigoniidae remain conspicuously absent. Our study aims to fill a critical gap in our understanding of satellitome evolutionary processes shaping Tettigoniidae genomes. Repeatome analysis revealed that the Meconema thalassinum genome comprises 92%, and Phryganogryllacris superangulata had the lowest value of 34%, with an average of 67% in other Tettigoniidae species. The analysis reveals significant variation in the number of satellite DNA repeats across species of the Tettigoniidae family, with M. thalassinum exhibiting the highest count, 246, reported in insects to date and the lowest count, 10, in Pholidoptera griseoptera. Ruspolia dubia and Ruspolia yunnana, which are congeneric species, showcase distinct counts of 104 and 84 families, respectively. Satellite DNA repeats in R. dubia exhibit the highest abundance, constituting 17.2% of the total genome, while the lowest abundance was reported in P. griseoptera, at 5.65%. The genome size correlates weakly with the satellite DNA family count (rs = 0.42, p = 0.29), but a strong correlation exists between satellite abundance and family number (rs = 0.73, p = 0.03). Moreover, the analysis of satellite DNA gain and loss patterns provides insights into the amplification and homogenization of satellite DNA families within the genome, with species-specific repeats exhibiting a positive trend toward amplification. The chromosomal distribution in M. thalassinum displayed that the highest accumulation was observed on Chr12, Chr01, and Chr04, constituting 17.79%, 17.4%, and 17.22% of the total chromosome size, respectively. The chromosome-specific propagation of satellite DNA families was evident, with MthSat01 solely on chromosome 1 and MthSat170 on chromosome 2, sharing 1.64% and 2.33%. The observed conservation and variations in satellite DNA number and abundances, along with distinct patterns of gain and loss, indicate the influence of potentially diverse evolutionary processes shaping the genomic landscape of these insects, which requires further investigation. Furthermore, the differential accumulation of satellite DNA on specific chromosomes implies that potential chromosome-specific functions or structural features influence the retention and proliferation of satellite sequences.

Comprehensive analysis of the Xya riparia genome uncovers the dominance of DNA transposons, LTR/Gypsy elements, and their evolutionary dynamics

Transposable elements (TEs) are DNA sequences that can move or replicate within a genome, and their study has become increasingly important in understanding genome evolution and function. The Tridactylidae family, including Xya riparia (pygmy mole cricket), harbors a variety of transposable elements (TEs) that have been insufficiently investigated. Further research is required to fully understand their diversity and evolutionary characteristics. Hence, we conducted a comprehensive repeatome analysis of X. riparia species using the chromosome-level assembled genome. The study aimed to comprehensively analyze the abundance, distribution, and age of transposable elements (TEs) in the genome. The results indicated that the genome was 1.67 Gb, with 731.63 Mb of repetitive sequences, comprising 27% of Class II (443.25 Mb), 16% of Class I (268.45 Mb), and 1% of unknown TEs (19.92 Mb). The study found that DNA transposons dominate the genome, accounting for approximately 60% of the total repeat size, with retrotransposons and unknown elements accounting for 37% and 3% of the genome, respectively. The members of the Gypsy superfamily were the most abundant amongst retrotransposons, accounting for 63% of them. The transposable superfamilies (LTR/Gypsy, DNA/nMITE, DNA/hAT, and DNA/Helitron) collectively constituted almost 70% of the total repeat size of all six chromosomes. The study further unveiled a significant linear correlation (Pearson correlation: r = 0.99, p-value = 0.00003) between the size of the chromosomes and the repetitive sequences. The average age of DNA transposon and retrotransposon insertions ranges from 25 My (million years) to 5 My. The satellitome analysis discovered 13 satellite DNA families that comprise about 0.15% of the entire genome. In addition, the transcriptional analysis of TEs found that DNA transposons were more transcriptionally active than retrotransposons. Overall, the study suggests that the genome of X. riparia is complex, characterized by a substantial portion of repetitive elements. These findings not only enhance our understanding of TE evolution within the Tridactylidae family but also provide a foundation for future investigations into the genomic intricacies of related species.

Genome sizes and repeatome evolution in zoantharians (Cnidaria: Hexacorallia: Zoantharia)

Across eukaryotes, large variations of genome sizes have been observed even between closely related species. Transposable elements as part of the repeated DNA have been proposed and confirmed as one of the most important contributors to genome size variation. However, the evolutionary implications of genome size variation and transposable element dynamics are not well understood. Together with phenotypic traits, they are commonly referred to as the "C-value enigma". The order Zoantharia are benthic cnidarians found from intertidal zones to the deep sea, and some species are particularly abundant in coral reefs. Despite their high ecological relevance, zoantharians have yet to be largely studied from the genomic point of view. This study aims at investigating the role of the repeatome (total content of repeated elements) in genome size variations across the order Zoantharia. To this end, whole-genomes of 32 zoantharian species representing five families were sequenced. Genome sizes were estimated and the abundances of different repeat classes were assessed. In addition, the repeat overlap between species was assessed by a sequence clustering method. The genome sizes in the dataset varied up to 2.4 fold magnitude. Significant correlations between genome size, repeated DNA content and transposable elements, respectively (Pearson's correlation test R2 = 0.47, p = 0.0016; R2 = 0.22, p = 0.05) were found, suggesting their involvement in the dynamics of genome expansion and reduction. In all species, long interspersed nuclear elements and DNA transposons were the most abundant identified elements. These transposable elements also appeared to have had a recent expansion event. This was in contrast to the comparative clustering analysis which revealed species-specific patterns of satellite elements' amplification. In summary, the genome sizes of zoantharians likely result from the complex dynamics of repeated elements. Finally, the majority of repeated elements (up to 70%) could not be annotated to a known repeat class, highlighting the need to further investigate non-model cnidarian genomes. More research is needed to understand how repeated DNA dynamics relate to zoantharian evolution and their biology.

Genome sequencing provides insights into the evolution of gene families encoding plant cell wall-degrading enzymes in longhorned beetles

With more than 36,000 species, the longhorned beetles (family Cerambycidae) are a mega-diverse lineage of mostly xylophagous insects, all of which are represented by the sole sequenced genome of the Asian longhorned beetle (Anoplophora glabripennis; Lamiinae). Their successful radiation has been linked to their ability to degrade plant cell wall components using a range of so-called plant cell wall-degrading enzymes (PCWDEs). Our previous analysis of larval gut transcriptomes demonstrated that cerambycid beetles horizontally acquired genes encoding PCWDEs from various microbial donors; these genes evolved through multiple duplication events to form gene families. To gain further insights into the evolution of these gene families during the Cerambycidae radiation, we assembled draft genomes for four beetle species belonging to three subfamilies using long-read nanopore sequencing. All the PCWDE-encoding genes we annotated from the corresponding larval gut transcriptomes were present in these draft genomes. We confirmed that the newly discovered horizontally acquired glycoside hydrolase family 7 (GH7), subfamily 26 of GH43 (GH43_26), and GH53 (all of which are absent from the A. glabripennis genome) were indeed encoded by these beetles' genome. Most of the PCWDE-encoding genes of bacterial origin gained introns after their transfer into the beetle genome. Altogether, we show that draft genome assemblies generated from nanopore long-reads offer meaningful information to the study of the evolution of gene families in insects. We anticipate that our data will support studies aiming to better understand the biology of the Cerambycidae and other beetles in general.

The draft genome of the microscopic Nemertoderma westbladi sheds light on the evolution of Acoelomorpha genomes

Background: Xenacoelomorpha is a marine clade of microscopic worms that is an important model system for understanding the evolution of key bilaterian novelties, such as the excretory system. Nevertheless, Xenacoelomorpha genomics has been restricted to a few species that either can be cultured in the lab or are centimetres long. Thus far, no genomes are available for Nemertodermatida, one of the group's main clades and whose origin has been dated more than 400 million years ago. Methods: DNA was extracted from a single specimen and sequenced with HiFi following the PacBio Ultra-Low DNA Input protocol. After genome assembly, decontamination, and annotation, the genome quality was benchmarked using two acoel genomes and one Illumina genome as reference. The gene content of three cnidarians, three acoelomorphs, four deuterostomes, and eight protostomes was clustered in orthogroups to make inferences of gene content evolution. Finally, we focused on the genes related to the ultrafiltration excretory system to compare patterns of presence/absence and gene architecture among these clades. Results: We present the first nemertodermatid genome sequenced from a single specimen of Nemertoderma westbladi. Although genome contiguity remains challenging (N50: 60 kb), it is very complete (BUSCO: 80.2%, Metazoa; 88.6%, Eukaryota) and the quality of the annotation allows fine-detail analyses of genome evolution. Acoelomorph genomes seem to be relatively conserved in terms of the percentage of repeats, number of genes, number of exons per gene and intron size. In addition, a high fraction of genes present in both protostomes and deuterostomes are absent in Acoelomorpha. Interestingly, we show that all genes related to the excretory system are present in Xenacoelomorpha except Osr, a key element in the development of these organs and whose acquisition seems to be interconnected with the origin of the specialised excretory system. Conclusion: Overall, these analyses highlight the potential of the Ultra-Low Input DNA protocol and HiFi to generate high-quality genomes from single animals, even for relatively large genomes, making it a feasible option for sequencing challenging taxa, which will be an exciting resource for comparative genomics analyses.

High-quality haploid genomes corroborate 29 chromosomes and highly conserved synteny of genes in Hyles hawkmoths (Lepidoptera: Sphingidae)

BACKGROUND

Morphological and traditional genetic studies of the young Pliocene genus Hyles have led to the understanding that despite its importance for taxonomy, phenotypic similarity of wing patterns does not correlate with phylogenetic relationship. To gain insights into various aspects of speciation in the Spurge Hawkmoth (Hyles euphorbiae), we assembled a chromosome-level genome and investigated some of its characteristics.

RESULTS

The genome of a male H. euphorbiae was sequenced using PacBio and Hi-C data, yielding a 504 Mb assembly (scaffold N50 of 18.2 Mb) with 99.9% of data represented by the 29 largest scaffolds forming the haploid chromosome set. Consistent with this, FISH analysis of the karyotype revealed n = 29 chromosomes and a WZ/ZZ (female/male) sex chromosome system. Estimates of chromosome length based on the karyotype image provided an additional quality metric of assembled chromosome size. Rescaffolding the published male H. vespertilio genome resulted in a high-quality assembly (651 Mb, scaffold N50 of 22 Mb) with 98% of sequence data in the 29 chromosomes. The larger genome size of H. vespertilio (average 1C DNA value of 562 Mb) was accompanied by a proportional increase in repeats from 45% in H. euphorbiae (measured as 472 Mb) to almost 55% in H. vespertilio. Several wing pattern genes were found on the same chromosomes in the two species, with varying amounts and positions of repetitive elements and inversions possibly corrupting their function.

CONCLUSIONS

Our two-fold comparative genomics approach revealed high gene synteny of the Hyles genomes to other Sphingidae and high correspondence to intact Merian elements, the ancestral linkage groups of Lepidoptera, with the exception of three simple fusion events. We propose a standardized approach for genome taxonomy using nucleotide homology via scaffold chaining as the primary tool combined with Oxford plots based on Merian elements to infer and visualize directionality of chromosomal rearrangements. The identification of wing pattern genes promises future understanding of the evolution of forewing patterns in the genus Hyles, although further sequencing data from more individuals are needed. The genomic data obtained provide additional reliable references for further comparative studies in hawkmoths (Sphingidae).

Satellite DNAs-From Localized to Highly Dispersed Genome Components

According to the established classical view, satellite DNAs are defined as abundant non-coding DNA sequences repeated in tandem that build long arrays located in heterochromatin. Advances in sequencing methodologies and development of specialized bioinformatics tools enabled defining a collection of all repetitive DNAs and satellite DNAs in a genome, the repeatome and the satellitome, respectively, as well as their reliable annotation on sequenced genomes. Supported by various non-model species included in recent studies, the patterns of satellite DNAs and satellitomes as a whole showed much more diversity and complexity than initially thought. Differences are not only in number and abundance of satellite DNAs but also in their distribution across the genome, array length, interspersion patterns, association with transposable elements, localization in heterochromatin and/or in euchromatin. In this review, we compare characteristic organizational features of satellite DNAs and satellitomes across different animal and plant species in order to summarize organizational forms and evolutionary processes that may lead to satellitomes' diversity and revisit some basic notions regarding repetitive DNA landscapes in genomes.

New estimates of genome size in Orthoptera and their evolutionary implications

Animal genomes vary widely in size, and much of their architecture and content remains poorly understood. Even among related groups, such as orders of insects, genomes may vary in size by orders of magnitude-for reasons unknown. The largest known insect genomes were repeatedly found in Orthoptera, e.g., Podisma pedestris (1C = 16.93 pg), Stethophyma grossum (1C = 18.48 pg) and Bryodemella holdereri (1C = 18.64 pg). While all these species belong to the suborder of Caelifera, the ensiferan Deracantha onos (1C = 19.60 pg) was recently found to have the largest genome. Here, we present new genome size estimates of 50 further species of Ensifera (superfamilies Gryllidea, Tettigoniidea) and Caelifera (Acrididae, Tetrigidae) based on flow cytometric measurements. We found that Bryodemella tuberculata (Caelifera: Acrididae) has the so far largest measured genome of all insects with 1C = 21.96 pg (21.48 gBp). Species of Orthoptera with 2n = 16 and 2n = 22 chromosomes have significantly larger genomes than species with other chromosome counts. Gryllidea genomes vary between 1C = 0.95 and 2.88 pg, and Tetrigidae between 1C = 2.18 and 2.41, while the genomes of all other studied Orthoptera range in size from 1C = 1.37 to 21.96 pg. Reconstructing ancestral genome sizes based on a phylogenetic tree of mitochondrial genomic data, we found genome size values of >15.84 pg only for the nodes of Bryodemella holdereri / B. tuberculata and Chrysochraon dispar / Euthystira brachyptera. The predicted values of ancestral genome sizes are 6.19 pg for Orthoptera, 5.37 pg for Ensifera, and 7.28 pg for Caelifera. The reasons for the large genomes in Orthoptera remain largely unknown, but a duplication or polyploidization seems unlikely as chromosome numbers do not differ much. Sequence-based genomic studies may shed light on the underlying evolutionary mechanisms.

Combined CRISPR toolkits reveal the domestication landscape and function of the ultra-long and highly repetitive silk genes

Highly repetitive sequences play a major structural and function role in the genome. In the present study, we developed Cas9-assisted cloning and SMRT sequencing of long repetitive sequences (CACS) to sequence and manipulate highly repetitive genes from eukaryotic genomes. CACS combined Cas9-mediated cleavage of a target segment from an intact genome, Gibson assembly cloning, and PacBio SMRT sequencing. Applying CACS, we directly cloned and sequenced the complete sequences of fibroin heavy chain (FibH) genes from 17 domesticated (Bombyx mori) and 7 wild (Bombyx mandarina) silkworms. Our analysis revealed the unique fine structure organization, genetic variations, and domestication dynamics of FibH. We also demonstrated that the length of the repetitive regions determined the mechanical properties of silk fiber, which was further confirmed by Cas9 editing of FibH. CACS is a simple, robust, and efficient approach, providing affordable accessibility to highly repetitive regions of a genome. STATEMENT OF SIGNIFICANCE: Silkworm silk is the earliest and most widely used animal fiber, and its excellent performance mainly depends on the fibroin heavy chain (FibH) protein. The FibH gene is the main breakthrough in understanding the formation mechanism and improvement of silk fiber. In the study, we developed a CACS method for characterizing the fine structure and domestication landscape of 24 silkworm FibH genes. We used CRISPR/Cas9 to edit the repetitive sequence of FibH genes, revealing the relationship between FibH genes and mechanical properties of silkworm silk. Our study is helpful in modifying silk genes to manipulate other valuable highly repetitive sequences, and provides insight for silkworm breeding.

Assembly-Free Detection and Quantification of Transposable Elements with dnaPipeTE

The detection and quantification of transposable elements (TE) are notoriously challenging despite their relevance in evolutionary genomics and molecular ecology. The main hurdle is caused by the dependence of numerous tools on genome assemblies, whose level of completion directly affects the comparability of the results across species or populations. dnaPipeTE, whose use is demonstrated here, tackles this issue by directly performing TE detection, classification, and quantification from unassembled short reads. This chapter details all the required steps to perform a comparative analysis of the TE content between two related species, starting from the installation of a recently containerized version of the program to the post-processing of the outputs.

Transposons and non-coding regions drive the intrafamily differences of genome size in insects

Genome size (GS) can vary considerably between phylogenetically close species, but the landscape of GS changes in insects remain largely unclear. To better understand the specific evolutionary factors that determine GS in insects, we examined flow cytometry-based published GS data from 1,326 insect species, spanning 700 genera, 155 families, and 21 orders. Model fitting showed that GS generally followed an Ornstein-Uhlenbeck adaptive evolutionary model in Insecta overall. Ancestral reconstruction indicated a likely GS of 1,069 Mb, suggesting that most insect clades appeared to undergo massive genome expansions or contractions. Quantification of genomic components in 56 species from nine families in four insect orders revealed that the proliferation of transposable elements contributed to high variation in GS between close species, such as within Coleoptera. This study sheds lights on the pattern of GS variation in insects and provides a better understanding of insect GS evolution.

Myriapod genomes reveal ancestral horizontal gene transfer and hormonal gene loss in millipedes

Animals display a fascinating diversity of body plans. Correspondingly, genomic analyses have revealed dynamic evolution of gene gains and losses among animal lineages. Here we sequence six new myriapod genomes (three millipedes, three centipedes) at key phylogenetic positions within this major but understudied arthropod lineage. We combine these with existing genomic resources to conduct a comparative analysis across all available myriapod genomes. We find that millipedes generally have considerably smaller genomes than centipedes, with the repeatome being a major contributor to genome size, driven by independent large gains of transposons in three centipede species. In contrast to millipedes, centipedes gained a large number of gene families after the subphyla diverged, with gains contributing to sensory and locomotory adaptations that facilitated their ecological shift to predation. We identify distinct horizontal gene transfer (HGT) events from bacteria to millipedes and centipedes, with no identifiable HGTs shared among all myriapods. Loss of juvenile hormone O-methyltransferase, a key enzyme in catalysing sesquiterpenoid hormone production in arthropods, was also revealed in all millipede lineages. Our findings suggest that the rapid evolution of distinct genomic pathways in centipede and millipede lineages following their divergence from the myriapod ancestor, was shaped by differing ecological pressures.

Myriapods play an important ecological role in soil and forest ecosystems. Here the authors analyse nine myriapod genomes, showing rapid evolution of distinct genomic pathways in centipede and millipede lineages, shaped by differing ecological pressures.

Satellitome Analysis and Transposable Elements Comparison in Geographically Distant Populations of Spodoptera frugiperda

Spodoptera frugiperda (fall armyworm) is a member of the superfamily Noctuoidea that accounts for more than a third of all Lepidoptera and includes a considerable number of agricultural and forest pest species. Spodoptera frugiperda is a polyphagous species that is a significant agricultural pest worldwide, emphasizing its economic importance. Spodoptera frugiperda's genome size, assembly, phylogenetic classification, and transcriptome analysis have all been previously described. However, the different studies reported different compositions of repeated DNA sequences that occupied the whole assembled genome, and the Spodoptera frugiperda genome also lacks the comprehensive study of dynamic satellite DNA. We conducted a comparative analysis of repetitive DNA across geographically distant populations of Spodoptera frugiperda, particularly satellite DNA, using publicly accessible raw genome data from eight different geographical regions. Our results showed that most transposable elements (TEs) were commonly shared across all geographically distant samples, except for the Maverick and PIF/Harbinger elements, which have divergent repeat copies. The TEs age analysis revealed that most TEs families consist of young copies 1-15 million years old; however, PIF/Harbinger has some older/degenerated copies of 30-35 million years old. A total of seven satellite DNA families were discovered, accounting for approximately 0.65% of the entire genome of the Spodoptera frugiperda fall armyworm. The repeat profiling analysis of satellite DNA families revealed differential read depth coverage or copy numbers. The satellite DNA families range in size from the lowest 108 bp SfrSat06-108 families to the largest (1824 bp) SfrSat07-1824 family. We did not observe a statistically significant correlation between monomer length and K2P divergence, copy number, or abundance of each satellite family. Our findings suggest that the satellite DNA families identified in Spodoptera frugiperda account for a considerable proportion of the genome's repetitive fraction. The satellite DNA families' repeat profiling revealed a point mutation along the reference sequences. Limited TEs differentiation exists among geographically distant populations of Spodoptera frugiperda.

Comparative Analysis of Transposable Elements in Genus Calliptamus Grasshoppers Revealed That Satellite DNA Contributes to Genome Size Variation

Simple Summary

Calliptamus is a genus of grasshoppers belonging to the family Acrididae. The genus Calliptamus includes approximately 17 recognized species. Calliptamus abbreviatus, Calliptamus italicus, and Calliptamus barbarus are three species that are widely found in northern China. These species are polyphagous, feeding on a variety of wild plants as well as crops, particularly legumes. The genome sizes, phylogenetic position, and transcriptome analysis of the genus Calliptamus were already known previous to this research. The repeatome analysis of these species was missing, which is directly linked to the larger genome sizes of the grasshoppers. Here, we classified repetitive DNA sequences at the level of superfamilies and sub-families, and found that LINE, TcMar-Tc1 and Ty3-gypsy LTR retrotransposons dominated the repeatomes of all genomes, accounting for 16–34% of the total genomes of these species. Satellite DNA dynamic evolutionary changes in all three genomes played a role in genome size evolution. This study would be a valuable source for future genome assemblies.

Abstract

Transposable elements (TEs) play a significant role in both eukaryotes and prokaryotes genome size evolution, structural changes, duplication, and functional variabilities. However, the large number of different repetitive DNA has hindered the process of assembling reference genomes, and the genus level TEs diversification of the grasshopper massive genomes is still under investigation. The genus Calliptamus diverged from Peripolus around 17 mya and its species divergence dated back about 8.5 mya, but their genome size shows rather large differences. Here, we used low-coverage Illumina unassembled short reads to investigate the effects of evolutionary dynamics of satDNAs and TEs on genome size variations. The Repeatexplorer2 analysis with 0.5X data resulted in 52%, 56%, and 55% as repetitive elements in the genomes of Calliptamus barbarus, Calliptamus italicus, and Calliptamus abbreviatus, respectively. The LINE and Ty3-gypsy LTR retrotransposons and TcMar-Tc1 dominated the repeatomes of all genomes, accounting for 16–35% of the total genomes of these species. Comparative analysis unveiled that most of the transposable elements (TEs) except satDNAs were highly conserved across three genomes in the genus Calliptamus grasshoppers. Out of a total of 20 satDNA families, 17 satDNA families were commonly shared with minor variations in abundance and divergence between three genomes, and 3 were Calliptamus barbarus specific. Our findings suggest that there is a significant amplification or contraction of satDNAs at genus phylogeny which is the main cause that made genome size different.

Comparative analysis reveals within-population genome size variation in a rotifer is driven by large genomic elements with highly abundant satellite DNA repeat elements

BACKGROUND

Eukaryotic genomes are known to display an enormous variation in size, but the evolutionary causes of this phenomenon are still poorly understood. To obtain mechanistic insights into such variation, previous studies have often employed comparative genomics approaches involving closely related species or geographically isolated populations within a species. Genome comparisons among individuals of the same population remained so far understudied-despite their great potential in providing a microevolutionary perspective to genome size evolution. The rotifer Brachionus asplanchnoidis represents one of the most extreme cases of within-population genome size variation among eukaryotes, displaying almost twofold variation within a geographic population.

RESULTS

Here, we used a whole-genome sequencing approach to identify the underlying DNA sequence differences by assembling a high-quality reference genome draft for one individual of the population and aligning short reads of 15 individuals from the same geographic population including the reference individual. We identified several large, contiguous copy number variable regions (CNVs), up to megabases in size, which exhibited striking coverage differences among individuals, and whose coverage overall scaled with genome size. CNVs were of remarkably low complexity, being mainly composed of tandemly repeated satellite DNA with only a few interspersed genes or other sequences, and were characterized by a significantly elevated GC-content. CNV patterns in offspring of two parents with divergent genome size and CNV patterns in several individuals from an inbred line differing in genome size demonstrated inheritance and accumulation of CNVs across generations.

CONCLUSIONS

By identifying the exact genomic elements that cause within-population genome size variation, our study paves the way for studying genome size evolution in contemporary populations rather than inferring patterns and processes a posteriori from species comparisons.

The salmon louse genome: Copepod features and parasitic adaptations

Copepods encompass numerous ecological roles including parasites, detrivores and phytoplankton grazers. Nonetheless, copepod genome assemblies remain scarce. Lepeophtheirus salmonis is an economically and ecologically important ectoparasitic copepod found on salmonid fish. We present the 695.4 Mbp L. salmonis genome assembly containing ≈60% repetitive regions and 13,081 annotated protein-coding genes. The genome comprises 14 autosomes and a ZZ-ZW sex chromosome system. Assembly assessment identified 92.4% of the expected arthropod genes. Transcriptomics supported annotation and indicated a marked shift in gene expression after host attachment, including apparent downregulation of genes related to circadian rhythm coinciding with abandoning diurnal migration. The genome shows evolutionary signatures including loss of genes needed for peroxisome biogenesis, presence of numerous FNII domains, and an incomplete heme homeostasis pathway suggesting heme proteins to be obtained from the host. Despite repeated development of resistance against chemical treatments L. salmonis exhibits low numbers of many genes involved in detoxification.

Impact of Repetitive DNA Elements on Snake Genome Biology and Evolution

The distinctive biology and unique evolutionary features of snakes make them fascinating model systems to elucidate how genomes evolve and how variation at the genomic level is interlinked with phenotypic-level evolution. Similar to other eukaryotic genomes, large proportions of snake genomes contain repetitive DNA, including transposable elements (TEs) and satellite repeats. The importance of repetitive DNA and its structural and functional role in the snake genome, remain unclear. This review highlights the major types of repeats and their proportions in snake genomes, reflecting the high diversity and composition of snake repeats. We present snakes as an emerging and important model system for the study of repetitive DNA under the impact of sex and microchromosome evolution. We assemble evidence to show that certain repetitive elements in snakes are transcriptionally active and demonstrate highly dynamic lineage-specific patterns as repeat sequences. We hypothesize that particular TEs can trigger different genomic mechanisms that might contribute to driving adaptive evolution in snakes. Finally, we review emerging approaches that may be used to study the expression of repetitive elements in complex genomes, such as snakes. The specific aspects presented here will stimulate further discussion on the role of genomic repeats in shaping snake evolution.

Simple inheritance of color and pattern polymorphism in the steppe grasshopper Chorthippus dorsatus

The green-brown polymorphism of grasshoppers and bush-crickets represents one of the most penetrant polymorphisms in any group of organisms. This poses the question of why the polymorphism is shared across species and how it is maintained. There is mixed evidence for whether and in which species it is environmentally or genetically determined in Orthoptera. We report breeding experiments with the steppe grasshopper Chorthippus dorsatus, a polymorphic species for the presence and distribution of green body parts. Morph ratios did not differ between sexes, and we find no evidence that the rearing environment (crowding and habitat complexity) affected the polymorphism. However, we find strong evidence for genetic determination for the presence/absence of green and its distribution. Results are most parsimoniously explained by three autosomal loci with two alleles each and simple dominance effects: one locus influencing the ability to show green color, with a dominant allele for green; a locus with a recessive allele suppressing green on the dorsal side; and a locus with a recessive allele suppressing green on the lateral side. Our results contribute to the emerging contrast between the simple genetic inheritance of green-brown polymorphisms in the subfamily Gomphocerinae and environmental determination in other subfamilies of grasshoppers. In three out of four species of Gomphocerinae studied so far, the results suggest one or a few loci with a dominance of alleles allowing the occurrence of green. This supports the idea that brown individuals differ from green individuals by homozygosity for loss-of-function alleles preventing green pigment production or deposition.

The Evolutionary Patterns of Genome Size in Ensifera (Insecta: Orthoptera)

Genomic size variation has long been a focus for biologists. However, due to the lack of genome size data, the mechanisms behind this variation and the biological significance of insect genome size are rarely studied systematically. The detailed taxonomy and phylogeny of the Ensifera, as well as the extensive documentation concerning their morphological, ecological, behavioral, and distributional characteristics, make them a strong model for studying the important scientific problem of genome size variation. However, data on the genome size of Ensifera are rather sparse. In our study, we used flow cytometry to determine the genome size of 32 species of Ensifera, the smallest one being only 1C = 0.952 pg with the largest species up to 1C = 19.135 pg, representing a 20-fold range. This provides a broader blueprint for the genome size variation of Orthoptera than was previously available. We also completed the assembly of nine mitochondrial genomes and combined mitochondrial genome data from public databases to construct phylogenetic trees containing 32 species of Ensifera and three outgroups. Based on these inferred phylogenetic trees, we detected the phylogenetic signal of genome size variation in Ensifera and found that it was strong in both males and females. Phylogenetic comparative analyses revealed that there were no correlations between genome size and body size or flight ability in Tettigoniidae. Reconstruction of ancestral genome size revealed that the genome size of Ensifera evolved in a complex pattern, in which the genome size of the grylloid clade tended to decrease while that of the non-grylloid clade expanded significantly albeit with fluctuations. However, the evolutionary mechanisms underlying variation of genome size in Ensifera are still unknown.