A guided tour of large genome size in animals: what we know and where we are heading

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.

Shark genome size evolution and its relationship with cellular, life-history, ecological, and diversity traits

Among vertebrates, sharks exhibit both large and heterogeneous genome sizes ranging from 2.86 to 17.05 pg. Aiming for a better understanding of the patterns and causalities of shark genome size evolution, we applied phylogenetic comparative methods to published genome-size estimates for 71 species representing the main phylogenetic lineages, life-histories and ecological traits. The sixfold range of genome size variation was strongly traceable throughout the phylogeny, with a major expansion preceding shark diversification during the late Paleozoic and an ancestral state (6.33 pg) close to the present-day average (6.72 pg). Subsequent deviations from this average occurred at higher rates in squalomorph than in galeomorph sharks and were unconnected to evolutionary changes in the karyotype architecture, which were dominated by descending disploidy events. Genome size was positively correlated with cell and nucleus sizes and negatively with metabolic rate. The metabolic constraints on increasing genome size also manifested at higher phenotypic scales, with large genomes associated with slow lifestyles and purely marine waters. Moreover, large genome sizes were also linked to non-placental reproductive modes, which may entail metabolically less demanding embryological developments. Contrary to ray-finned fishes, large genome size was associated neither with the taxonomic diversity of affected clades nor with low genetic diversity.

Effects of Temperature Adaptation on the Metabolism and Physiological Properties of Sturgeon Fish Larvae Cell Line

This study investigated how Atlantic sturgeon cells respond to elevated temperatures, shedding light on the potential impacts of climate change on fish. Atlantic sturgeon (Acipenser oxyrinchus), an IUCN (International Union for Conservation of Nature) Red List species and evolutionarily related to paleonisiform species, may have considerable physiological adaptability, suggesting that this species may be able to cope with changing climatic conditions and higher temperatures. To test this hypothesis, the AOXlar7y cell line was examined at 20 °C (control) and at elevated temperatures of 25 °C and 28 °C. Parameters including proliferation, vitality, morphology, and gene expressions related to proliferation, stemness, and stress were evaluated. Additionally, to achieve a comprehensive understanding of cellular changes, mitochondrial and metabolic activities were assessed using Seahorse XF96. AOXlar7y cells adapted to 28 °C exhibited enhanced mitochondrial adaptability, plasticity, heightened cell proliferation, and increased hsp70 expression. Increased baseline respiration indicated elevated ATP demand, which is potentially linked to higher cell proliferation and heat stress defense. Cells at 28 °C also displayed elevated reserve respiration capacity, suggesting adaptation to energy demands. At 25 °C, AOXlar7y cells showed no changes in basal respiration or mitochondrial capacity, suggesting unchanged ATP demand compared to cells cultivated at 20 °C. Proliferation and glycolytic response to energy requirements were diminished, implying a connection between glycolysis inhibition and proliferation suppression. These research results indicate sturgeon cells are capable of withstanding and adapting to an 8 °C temperature increase. This cellular analysis lays a foundation for future studies aimed at a deeper understanding of fish cell physiological adaptations, which will contribute to a better knowledge of environmental threats facing Atlantic sturgeon and fish populations amid climate change.

Ecological factors and parity mode correlate with genome size variation in squamate reptiles

BACKGROUND

Evidence of correlation between genome size, the nuclear haploid DNA content of a cell, environmental factors and life-history traits have been reported in many animal species. Genome size, however, spans over three orders of magnitude across taxa and such a correlation does not seem to follow a universal pattern. In squamate reptiles, the second most species-rich order of vertebrates, there are currently no studies investigating drivers of genome size variability. We run a series of phylogenetic generalized least-squares models on 227 species of squamates to test for possible relationships between genome size and ecological factors including latitudinal distribution, bioclimatic variables and microhabitat use. We also tested whether genome size variation can be associated with parity mode, a highly variable life history trait in this order of reptiles.

RESULTS

The best-fitting model showed that the interaction between microhabitat use and parity mode mainly accounted for genome size variation. Larger genome sizes were found in live-bearing species that live in rock/sand ecosystems and in egg-laying arboreal taxa. On the other hand, smaller genomes were found in fossorial live-bearing species.

CONCLUSIONS

Environmental factors and species parity mode appear to be among the main parameters explaining genome size variation in squamates. Our results suggest that genome size may favour adaptation of some species to certain environments or could otherwise result from the interaction between environmental factors and parity mode. Integration of genome size and genome sequencing data could help understand the role of differential genome content in the evolutionary process of genome size variation in squamates.

Next-generation development and application of codon model in evolution

To date, numerous nucleotide, amino acid, and codon substitution models have been developed to estimate the evolutionary history of any sequence/organism in a more comprehensive way. Out of these three, the codon substitution model is the most powerful. These models have been utilized extensively to detect selective pressure on a protein, codon usage bias, ancestral reconstruction and phylogenetic reconstruction. However, due to more computational demanding, in comparison to nucleotide and amino acid substitution models, only a few studies have employed the codon substitution model to understand the heterogeneity of the evolutionary process in a genome-scale analysis. Hence, there is always a question of how to develop more robust but less computationally demanding codon substitution models to get more accurate results. In this review article, the authors attempted to understand the basis of the development of different types of codon-substitution models and how this information can be utilized to develop more robust but less computationally demanding codon substitution models. The codon substitution model enables to detect selection regime under which any gene or gene region is evolving, codon usage bias in any organism or tissue-specific region and phylogenetic relationship between different lineages more accurately than nucleotide and amino acid substitution models. Thus, in the near future, these codon models can be utilized in the field of conservation, breeding and medicine.

Contrasting Phylogeographic Patterns of Mitochondrial and Genome-Wide Variation in the Groundwater Amphipod Crangonyx islandicus That Survived the Ice Age in Iceland

The analysis of phylogeographic patterns has often been based on mitochondrial DNA variation, but recent analyses dealing with nuclear DNA have in some instances revealed mito-nuclear discordances and complex evolutionary histories. These enigmatic scenarios, which may involve stochastic lineage sorting, ancestral hybridization, past dispersal and secondary contacts, are increasingly scrutinized with a new generation of genomic tools such as RADseq, which also poses additional analytical challenges. Here, we revisited the previously inconclusive phylogeographic history, showing the mito-nuclear discordance of an endemic groundwater amphipod from Iceland, Crangonyx islandicus, which is the only metazoan known to have survived the Pleistocene beneath the glaciers. Previous studies based on three DNA markers documented a mitochondrial scenario with the main divergence occurring between populations in northern Iceland and an ITS scenario with the main divergence between the south and north. We used double digest restriction-site-associated DNA sequencing (ddRADseq) to clarify this mito-nuclear discordance by applying several statistical methods while estimating the sensitivity to different analytical approaches (data-type, differentiation indices and base call uncertainty). A majority of nuclear markers and methods support the ITS divergence. Nevertheless, a more complex scenario emerges, possibly involving introgression led by male-biased dispersal among northern locations or mitochondrial capture, which may have been further strengthened by natural selection.

Fern Spores-"Ready-to-Use" Standards for Plant Genome Size Estimation Using a Flow Cytometric Approach

Spores and pollen of plants were used as flow cytometric materials to efficiently infer genome sizes. Given this advantage, they hold great potential for various flow cytometric applications, particularly as plant genome size standards. To develop such novel standards, we investigated conditions of pretreatment (bead vortex), buffer, and reliable genome sizes of three fern spore collections-Cibotium taiwanense "Kuo4395", Sphaeropteris lepifera "Tang0001", and Alsophila metteniana "Lee s.n.". Additionally, up to 30 year-old spore collections were obtained from herbarium specimens and from samples stored at 4 °C; their spore nuclei were extracted, and the quality and quantity of these nucleus extractions through storage ages were examined. Nuclear extractions with a longer bead vortex duration or lower spore/bead ratio generally resulted in a higher recovered quantity but a lower quality or purity. For each spore standard, the protocol optimization was determined by their performance in bead vortex conditions, and a 1C genome size was further inferred by linear regression (C. taiwanense "Kuo4395" = 5.058 pg; S. lepifera "Tang0001" = 7.117 pg; and A. metteniana "Lee s.n." = 19.379 pg). Spore nucleus quality and quantity are significantly negatively correlated with storage ages. Nuclear extractions of 10-year-old refrigerated spores remained qualified as a genome size standard; however, none of the herbarium spore collections fit such criteria. Our study is the first to develop and apply dried and refrigerated spores for genome size standards. These standards are ready to use, easy to manipulate, and feature long-term storage in comparison with traditionally used standards of fresh leaves.

Transposable element expansion and low-level piRNA silencing in grasshoppers may cause genome gigantism

Background

Transposable elements (TEs) have been likened to parasites in the genome that reproduce and move ceaselessly in the host, continuously enlarging the host genome. However, the Piwi-interacting RNA (piRNA) pathway defends animal genomes against the harmful consequences of TE invasion by imposing small-RNA-mediated silencing. Here we compare the TE activity of two grasshopper species with different genome sizes in Acrididae (Locusta migratoria manilensis♀1C = 6.60 pg, Angaracris rhodopa♀1C = 16.36 pg) to ascertain the influence of piRNAs.

Results

We discovered that repetitive sequences accounted for 74.56% of the genome in A. rhodopa, more than 56.83% in L. migratoria, and the large-genome grasshopper contained a higher TEs proportions. The comparative analysis revealed that 41 TEs (copy number > 500) were shared in both species. The two species exhibited distinct “landscapes” of TE divergence. The TEs outbreaks in the small-genome grasshopper occurred at more ancient times, while the large-genome grasshopper maintains active transposition events in the recent past. Evolutionary history studies on TEs suggest that TEs may be subject to different dynamics and resistances in these two species. We found that TE transcript abundance was higher in the large-genome grasshopper and the TE-derived piRNAs abundance was lower than in the small-genome grasshopper. In addition, we found that the piRNA methylase HENMT, which is underexpressed in the large-genome grasshopper, impedes the piRNA silencing to a lower level.

Conclusions

Our study revealed that the abundance of piRNAs is lower in the gigantic genome grasshopper than in the small genome grasshopper. In addition, the key gene HENMT in the piRNA biogenesis pathway (Ping-Pong cycle) in the gigantic genome grasshopper is underexpressed. We hypothesize that low-level piRNA silencing unbalances the original positive correlation between TEs and piRNAs, and triggers TEs to proliferate out of control, which may be one of the reasons for the gigantism of grasshopper genomes.

Body mass and cell size shape the tolerance of fishes to low oxygen in a temperature-dependent manner

Aerobic metabolism generates 15-20 times more energy (ATP) than anaerobic metabolism, which is crucial in maintaining energy budgets in animals, fueling metabolism, activity, growth and reproduction. For ectothermic water-breathers such as fishes, low dissolved oxygen may limit oxygen uptake and hence aerobic metabolism. Here, we assess, within a phylogenetic context, how abiotic and biotic drivers explain the variation in hypoxia tolerance observed in fishes. To do so, we assembled a database of hypoxia tolerance, measured as critical oxygen tensions (P_crit ) for 195 fish species. Overall, we found that hypoxia tolerance has a clear phylogenetic signal and is further modulated by temperature, body mass, cell size, salinity and metabolic rate. Marine fishes were more susceptible to hypoxia than freshwater fishes. This pattern is consistent with greater fluctuations in oxygen and temperature in freshwater habitats. Fishes with higher oxygen requirements (e.g. a high metabolic rate relative to body mass) also were more susceptible to hypoxia. We also found evidence that hypoxia and warming can act synergistically, as hypoxia tolerance was generally lower in warmer waters. However, we found significant interactions between temperature and the body and cell size of a fish. Constraints in oxygen uptake related to cellular surface area to volume ratios and effects of viscosity on the thickness of the boundary layers enveloping the gills could explain these thermal dependencies. The lower hypoxia tolerance in warmer waters was particularly pronounced for fishes with larger bodies and larger cell sizes. Previous studies have found a wide diversity in the direction and strength of relationships between P_crit and body mass. By including interactions with temperature, our study may help resolve these divergent findings, explaining the size dependency of hypoxia tolerance in fish.

Development of RNAi Methods for the Mormon Cricket, Anabrus simplex (Orthoptera: Tettigoniidae)

Mormon crickets are a major rangeland pest in the western United States and are currently managed by targeted applications of non-specific chemical insecticides, which can potentially have negative effects on the environment. In this study, we took the first steps toward developing RNAi methods for Mormon crickets as a potential alternative to traditional broad-spectrum insecticides. To design an effective RNAi-based insecticide, we first generated a de novo transcriptome for the Mormon cricket and developed dsRNAs that could silence the expression of seven housekeeping genes. We then characterized the RNAi efficiencies and time-course of knockdown using these dsRNAs, and assessed their ability to induce mortality. We have demonstrated that it is possible to elicit RNAi responses in the Mormon cricket by injection, but knockdown efficiencies and the time course of RNAi response varied according to target genes and tissue types. We also show that one of the reasons for the poor knockdown efficiencies could be the presence of dsRNA-degrading enzymes in the hemolymph. RNAi silencing is possible in Mormon cricket, but more work needs to be done before it can be effectively used as a population management method.

Ecological, physiological and life-history traits correlate with genome sizes in decapod crustaceans

Crustaceans are characterized by some of the most variable genome sizes among animals. Significant relationships between genome size and specific eco-physiological and morphological features have been described in many crustacean taxa, such as Amphipoda, Ostracoda, Cladocera, and Copepoda. A consistent pattern of genome size variation is yet to be found, however, in one of the most specious order of crustaceans, the Decapoda. We investigated how genome size is related to selected eco-physiological and life-history traits in species belonging to this order. We analyzed the respiratory, excretory, and developmental adaptations of 179 decapod species inhabiting shallow marine, deep marine, freshwater, intertidal, and terrestrial habitats. Our results show that the genome size is significantly larger in decapods that have a direct development. Moreover, in Anomura, Astacidea, and Brachyura we found larger genome sizes in species (i) living in freshwater habitats, (ii) using gills as a strictly water-breathing organs, and (iii) presenting a direct development. While species that (i) live in shallow waters, and intertidal and terrestrial habitats, (ii) have some degrees of air-breathing adaptations, namely gills that can exchange oxygen with air if wet or branchiostegal lungs, and (iii) possess an indirect or extended development are characterized by significantly smaller genomes. Our analyses show that developmental complexity affects genome size in decapods, and that multiple eco-physiological and life-history traits correlate with genome size in Anomura, Astacidea, and Brachyura.

Genome-wide survey and genetic characteristics of Ophichthus evermanni based on Illumina sequencing platform

Ophichthidae fishes limit to continental shelf of all tropical and subtropical oceans and contain more than 350 species, representing the greatest specialization diversity in the order Anguiliformes. In the present study, we conducted a genome survey sequencing (GSS) analysis of Ophichthus evermanni by Illumina sequencing platform to briefly reveal its genomic characteristics and phylogenetic relationship. The first de novo assembled 1.97 Gb draft genome of O. evermanni was predicted based on K-mer analysis without obvious nucleotide bias. The heterozygosity ratio was 0.70%, and the sequence repeat ratio was calculated to be 43.30%. A total of 9016 putative coding genes were successfully predicted, in which 3587 unigenes were identified by gene ontology (GO) analysis and 4375 unigenes were classified into cluster of orthologous groups for enkaryotic complete genomes (KOG) functional categories. About 2,812,813 microsatellite motifs including mono-, di-, tri-, tetra-, penta- and hexanucleotide motifs were identified, with an occurrence frequency of 23.32%. The most abundant type was dinucleotide repeat motifs, accounting for 49.19% of the total repeat types. The mitochondrial genome, as a byproduct of GSS, was assembled to investigate the evolutionary relationships between O. evermanni and its relatives. Bayesian inference (BI) phylogenetic tree inferring from concatenated 12 protein-coding genes (PCGs) showed complicated relationships among Ophichthidae species, indicating a polyphyletic origin of the family. The results would achieve more thorough genetic information of snake eels and provide a theoretical basis and reference for further genome-wide analysis of O. evermanni.

Genome size evolution in the diverse insect order Trichoptera

Background

Genome size is implicated in the form, function, and ecological success of a species. Two principally different mechanisms are proposed as major drivers of eukaryotic genome evolution and diversity: polyploidy (i.e., whole-genome duplication) or smaller duplication events and bursts in the activity of repetitive elements. Here, we generated de novo genome assemblies of 17 caddisflies covering all major lineages of Trichoptera. Using these and previously sequenced genomes, we use caddisflies as a model for understanding genome size evolution in diverse insect lineages.

Results

We detect a ∼14-fold variation in genome size across the order Trichoptera. We find strong evidence that repetitive element expansions, particularly those of transposable elements (TEs), are important drivers of large caddisfly genome sizes. Using an innovative method to examine TEs associated with universal single-copy orthologs (i.e., BUSCO genes), we find that TE expansions have a major impact on protein-coding gene regions, with TE-gene associations showing a linear relationship with increasing genome size. Intriguingly, we find that expanded genomes preferentially evolved in caddisfly clades with a higher ecological diversity (i.e., various feeding modes, diversification in variable, less stable environments).

Conclusion

Our findings provide a platform to test hypotheses about the potential evolutionary roles of TE activity and TE-gene associations, particularly in groups with high species, ecological, and functional diversities.

Dissecting cricket genomes for the advancement of entomology and entomophagy

Significant advances in biophysical methods such as next-generation sequencing technologies have now opened the way to conduct evolutionary and applied research based on the genomic information of greatly diverse insects. Crickets belonging to Orthoptera (Insecta: Polyneoptera), one of the most flourishing groups of insects, have contributed to the development of multiple scientific fields including developmental biology and neuroscience and have been attractive targets in evolutionary ecology for their diverse ecological niches. In addition, crickets have recently gained recognition as food and feed. However, the genomic information underlying their biological basis and application research toward breeding is currently underrepresented. In this review, we summarize the progress of genomics of crickets. First, we outline the phylogenetic position of crickets in insects and then introduce recent studies on cricket genomics and transcriptomics in a variety of fields. Furthermore, we present findings from our analysis of polyneopteran genomes, with a particular focus on their large genome sizes, chromosome number, and repetitive sequences. Finally, how the cricket genome can be beneficial to the food industry is discussed. This review is expected to enhance greater recognition of how important the cricket genomes are to the multiple biological fields and how basic research based on cricket genome information can contribute to tackling global food security.

Interacting Effects of Cell Size and Temperature on Gene Expression, Growth, Development and Swimming Performance in Larval Zebrafish

Cell size may be important in understanding the thermal biology of ectotherms, as the regulation and consequences of cell size appear to be temperature dependent. Using a recently developed model system of triploid zebrafish (which have around 1.5-fold larger cells than their diploid counterparts) we examine the effects of cell size on gene expression, growth, development and swimming performance in zebrafish larvae at different temperatures. Both temperature and ploidy affected the expression of genes related to metabolic processes (citrate synthase and lactate dehydrogenase), growth and swimming performance. Temperature also increased development rate, but there was no effect of ploidy level. We did find interactive effects between ploidy and temperature for gene expression, body size and swimming performance, confirming that the consequences of cell size are temperature dependent. Triploids with larger cells performed best at cool conditions, while diploids performed better at warmer conditions. These results suggest different selection pressures on ectotherms and their cell size in cold and warm habitats.

A haploid diamondback moth (Plutella xylostella L.) genome assembly resolves 31 chromosomes and identifies a diamide resistance mutation

The diamondback moth, Plutella xylostella (L.), is a highly mobile brassica crop pest with worldwide distribution and can rapidly evolve resistance to insecticides, including group 28 diamides. Reference genomes assembled using Illumina sequencing technology have provided valuable resources to advance our knowledge regarding the biology, origin and movement of diamondback moth, and more recently with its sister species, Plutella australiana. Here we apply a trio binning approach to sequence and annotate a chromosome level reference genome of P. xylostella using PacBio Sequel and Dovetail Hi-C sequencing technology and identify a point mutation that causes resistance to commercial diamides. A P. xylostella population collected from brassica crops in the Lockyer Valley, Australia (LV-R), was reselected for chlorantraniliprole resistance then a single male was crossed to a P. australiana female and a hybrid pupa sequenced. A chromosome level 328 Mb P. xylostella genome was assembled with 98.1% assigned to 30 autosomes and the Z chromosome. The genome was highly complete with 98.4% of BUSCO Insecta genes identified and RNAseq informed protein prediction annotated 19,002 coding genes. The LV-R strain survived recommended field application doses of chlorantraniliprole, flubendiamide and cyclaniliprole. Some hybrids also survived these doses, indicating significant departure from recessivity, which has not been previously documented for diamides. Diamide chemicals modulate insect Ryanodine Receptors (RyR), disrupting calcium homeostasis, and we identified an amino acid substitution (I4790K) recently reported to cause diamide resistance in a strain from Japan. This chromosome level assembly provides a new resource for insect comparative genomics and highlights the emergence of diamide resistance in Australia. Resistance management plans need to account for the fact that resistance is not completely recessive.

Correlated evolution of larval development, egg size and genome size across two genera of snapping shrimp

Across plants and animals, genome size is often correlated with life-history traits: large genomes are correlated with larger seeds, slower development, larger body size and slower cell division. Among decapod crustaceans, caridean shrimps are among the most variable both in terms of genome size variation and life-history characteristics such as larval development mode and egg size, but the extent to which these traits are associated in a phylogenetic context is largely unknown. In this study, we examine correlations among egg size, larval development and genome size in two different genera of snapping shrimp, Alpheus and Synalpheus, using phylogenetically informed analyses. In both Alpheus and Synalpheus, egg size is strongly linked to larval development mode: species with abbreviated development had significantly larger eggs than species with extended larval development. We produced the first comprehensive dataset of genome size in Alpheus (n = 37 species) and demonstrated that genome size was strongly and positively correlated with egg size in both Alpheus and Synalpheus. Correlated trait evolution analyses showed that in Alpheus, changes in genome size were clearly dependent on egg size. In Synalpheus, evolutionary path analyses suggest that changes in development mode (from extended to abbreviated) drove increases in egg volume; larger eggs, in turn, resulted in larger genomes. These data suggest that variation in reproductive traits may underpin the high degree of variation in genome size seen in a wide variety of caridean shrimp groups more generally.

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.

Eusociality in snapping shrimps is associated with larger genomes and an accumulation of transposable elements

Despite progress uncovering the genomic underpinnings of sociality, much less is known about how social living affects the genome. In different insect lineages, for example, eusocial species show both positive and negative associations between genome size and structure, highlighting the dynamic nature of the genome. Here, we explore the relationship between sociality and genome architecture in Synalpheus snapping shrimps that exhibit multiple origins of eusociality and extreme interspecific variation in genome size. Our goal is to determine whether eusociality leads to an accumulation of repetitive elements and an increase in genome size, presumably due to reduced effective population sizes resulting from a reproductive division of labor, or whether an initial accumulation of repetitive elements leads to larger genomes and independently promotes the evolution of eusociality through adaptive evolution. Using phylogenetically informed analyses, we find that eusocial species have larger genomes with more transposable elements (TEs) and microsatellite repeats than noneusocial species. Interestingly, different TE subclasses contribute to the accumulation in different species. Phylogenetic path analysis testing alternative causal relationships between sociality and genome architecture is most consistent with the hypothesis that TEs modulate the relationship between sociality and genome architecture. Although eusociality appears to influence TE accumulation, ancestral state reconstruction suggests moderate TE abundances in ancestral species could have fueled the initial transitions to eusociality. Ultimately, we highlight a complex and dynamic relationship between genome and social evolution, demonstrating that sociality can influence the evolution of the genome, likely through changes in demography related to patterns of reproductive skew.

Genome Size Covaries More Positively with Propagule Size than Adult Size: New Insights into an Old Problem

The body size and (or) complexity of organisms is not uniformly related to the amount of genetic material (DNA) contained in each of their cell nuclei ('genome size'). This surprising mismatch between the physical structure of organisms and their underlying genetic information appears to relate to variable accumulation of repetitive DNA sequences, but why this variation has evolved is little understood. Here, I show that genome size correlates more positively with egg size than adult size in crustaceans. I explain this and comparable patterns observed in other kinds of animals and plants as resulting from genome size relating strongly to cell size in most organisms, which should also apply to single-celled eggs and other reproductive propagules with relatively few cells that are pivotal first steps in their lives. However, since body size results from growth in cell size or number or both, it relates to genome size in diverse ways. Relationships between genome size and body size should be especially weak in large organisms whose size relates more to cell multiplication than to cell enlargement, as is generally observed. The ubiquitous single-cell 'bottleneck' of life cycles may affect both genome size and composition, and via both informational (genotypic) and non-informational (nucleotypic) effects, many other properties of multicellular organisms (e.g., rates of growth and metabolism) that have both theoretical and practical significance.

Shrinking body sizes in response to warming: explanations for the temperature-size rule with special emphasis on the role of oxygen

Body size is central to ecology at levels ranging from organismal fecundity to the functioning of communities and ecosystems. Understanding temperature-induced variations in body size is therefore of fundamental and applied interest, yet thermal responses of body size remain poorly understood. Temperature-size (T-S) responses tend to be negative (e.g. smaller body size at maturity when reared under warmer conditions), which has been termed the temperature-size rule (TSR). Explanations emphasize either physiological mechanisms (e.g. limitation of oxygen or other resources and temperature-dependent resource allocation) or the adaptive value of either a large body size (e.g. to increase fecundity) or a short development time (e.g. in response to increased mortality in warm conditions). Oxygen limitation could act as a proximate factor, but we suggest it more likely constitutes a selective pressure to reduce body size in the warm: risks of oxygen limitation will be reduced as a consequence of evolution eliminating genotypes more prone to oxygen limitation. Thus, T-S responses can be explained by the 'Ghost of Oxygen-limitation Past', whereby the resulting (evolved) T-S responses safeguard sufficient oxygen provisioning under warmer conditions, reflecting the balance between oxygen supply and demands experienced by ancestors. T-S responses vary considerably across species, but some of this variation is predictable. Body-size reductions with warming are stronger in aquatic taxa than in terrestrial taxa. We discuss whether larger aquatic taxa may especially face greater risks of oxygen limitation as they grow, which may be manifested at the cellular level, the level of the gills and the whole-organism level. In contrast to aquatic species, terrestrial ectotherms may be less prone to oxygen limitation and prioritize early maturity over large size, likely because overwintering is more challenging, with concomitant stronger end-of season time constraints. Mechanisms related to time constraints and oxygen limitation are not mutually exclusive explanations for the TSR. Rather, these and other mechanisms may operate in tandem. But their relative importance may vary depending on the ecology and physiology of the species in question, explaining not only the general tendency of negative T-S responses but also variation in T-S responses among animals differing in mode of respiration (e.g. water breathers versus air breathers), genome size, voltinism and thermally associated behaviour (e.g. heliotherms).

Are acute and acclimated thermal effects on metabolic rate modulated by cell size? A comparison between diploid and triploid zebrafish larvae

Being composed of small cells may carry energetic costs related to maintaining ionic gradients across cell membranes as well as benefits related to diffusive oxygen uptake. Here, we test the hypothesis that these costs and benefits of cell size in ectotherms are temperature dependent. To study the consequences of cell size for whole-organism metabolic rate, we compared diploid and triploid zebrafish larvae differing in cell size. A fully factorial design was applied combining three different rearing and test temperatures that allowed us to distinguish acute from acclimated thermal effects. Individual oxygen consumption rates of diploid and triploid larvae across declining levels of oxygen availability were measured. We found that both acute and acclimated thermal effects affected the metabolic response. In comparison with triploids, diploids responded more strongly to acute temperatures, especially when reared at the highest temperature. These observations support the hypothesis that animals composed of smaller cells (i.e. diploids) are less vulnerable to oxygen limitation in warm aquatic habitats. Furthermore, we found slightly improved hypoxia tolerance in diploids. By contrast, warm-reared triploids had higher metabolic rates when they were tested at acute cold temperature, suggesting that being composed of larger cells may provide metabolic advantages in the cold. We offer two mechanisms as a potential explanation of this result, related to homeoviscous adaptation of membrane function and the mitigation of developmental noise. Our results suggest that being composed of larger cells provides metabolic advantages in cold water, while being composed of smaller cells provides metabolic advantages in warm water.

Genome Size of 17 Species From Caelifera (Orthoptera) and Determination of Internal Standards With Very Large Genome Size in Insecta

Comparative studies of insect genome size show that Orthoptera is a unique group of Insecta with a significantly enlarged genome. To determine a suitable internal standard for a large genome and to compare the effects of different internal standards on estimates of genome size, we used four internal standards to estimate nuclear DNA content in nine insect species with large genomes. The results showed that the combination of two internal standards, Locusta migratoria (♂1C = 6.20 pg, ♀1C = 6.60 pg) and Periplaneta americana♂ (1C = 3.41 pg), was suitable for estimating large genome of Caelifera by flow cytometry. Using these two internal standards, we estimated the genome sizes of 17 species of Caelifera (12 genera in Acrididae, 2 genera in Pamphagidae, 1 genus in Pyrgomorphidae) using flow cytometry. Genomes ranged from 6.57 pg (Shirakiacris shirakii) to 18.64 pg (Bryodemella holdereri), the largest described in insects to date. These species showed significant genomic dimorphism based on sex: females had a 0.56 pg larger genome than males on average, which might be due to the sex chromosome determinism mechanism of X0(♂)/XX(♀). To test the results obtained by flow cytometry, we used k-mers of Illumina sequencing data to gauge the C-value of Calliptamus abbreviatus and Haplotropis brunneriana. The results of the two methods are slightly different. Genomes were estimated to be about 0.28 and 0.26 pg smaller, respectively, than the flow cytometry values. Furthermore, we also reconstructed the evolutionary relationships of these taxa and discuss the genome size evolution in a phylogenetic framework.

Using asexual vertebrates to study genome evolution and animal physiology: Banded (Fundulus diaphanus) x Common Killifish (F. heteroclitus) hybrid lineages as a model system

Wild, asexual, vertebrate hybrids have many characteristics that make them good model systems for studying how genomes evolve and epigenetic modifications influence animal physiology. In particular, the formation of asexual hybrid lineages is a form of reproductive incompatibility, but we know little about the genetic and genomic mechanisms by which this mode of reproductive isolation proceeds in animals. Asexual lineages also provide researchers with the ability to produce genetically identical individuals, enabling the study of autonomous epigenetic modifications without the confounds of genetic variation. Here, we briefly review the cellular and molecular mechanisms leading to asexual reproduction in vertebrates and the known genetic and epigenetic consequences of the loss of sex. We then specifically discuss what is known about asexual lineages of Fundulus diaphanus x F. heteroclitus to highlight gaps in our knowledge of the biology of these clones. Our preliminary studies of F. diaphanus and F. heteroclitus karyotypes from Porter's Lake (Nova Scotia, Canada) agree with data from other populations, suggesting a conserved interspecific chromosomal arrangement. In addition, genetic analyses suggest that: (a) the same major clonal lineage (Clone A) of F. diaphanus x F. heteroclitus has remained dominant over the past decade, (b) some minor clones have also persisted, (c) new clones may have recently formed, and iv) wild clones still mainly descend from F. diaphanus ♀ x F. heteroclitus ♂ crosses (96% in 2017-2018). These data suggest that clone formation may be a relatively rare, but continuous process, and there are persistent environmental or genetic factors causing a bias in cross direction. We end by describing our current research on the genomic causes and consequences of a transition to asexuality and the potential physiological consequences of epigenetic variation.

Triploidy in zebrafish larvae: Effects on gene expression, cell size and cell number, growth, development and swimming performance

There is renewed interest in the regulation and consequences of cell size adaptations in studies on understanding the ecophysiology of ectotherms. Here we test if induction of triploidy, which increases cell size in zebrafish (Danio rerio), makes for a good model system to study consequences of cell size. Ideally, diploid and triploid zebrafish should differ in cell size, but should otherwise be comparable in order to be suitable as a model. We induced triploidy by cold shock and compared diploid and triploid zebrafish larvae under standard rearing conditions for differences in genome size, cell size and cell number, development, growth and swimming performance and expression of housekeeping genes and hsp70.1. Triploid zebrafish have larger but fewer cells, and the increase in cell size matched the increase in genome size (+ 50%). Under standard conditions, patterns in gene expression, ontogenetic development and larval growth were near identical between triploids and diploids. However, under demanding conditions (i.e. the maximum swimming velocity during an escape response), triploid larvae performed poorer than their diploid counterparts, especially after repeated stimuli to induce swimming. This result is consistent with the idea that larger cells have less capacity to generate energy, which becomes manifest during repeated physical exertion resulting in increased fatigue. Triploidy induction in zebrafish appears a valid method to increase specifically cell size and this provides a model system to test for consequences of cell size adaptation for the energy budget and swimming performance of this ectothermic vertebrate.

Within-Population Genome Size Variation is Mediated by Multiple Genomic Elements That Segregate Independently during Meiosis

Within-species variation in genome size has been documented in many animals and plants. Despite its importance for understanding eukaryotic genome diversity, there is only sparse knowledge about how individual-level processes mediate genome size variation in populations. Here, we study a natural population of the rotifer Brachionus asplanchnoidis whose members differ up to 1.9-fold in diploid genome size, but were still able to interbreed and produce viable offspring. We show that genome size is highly heritable and can be artificially selected up or down, but not below a certain basal diploid genome size for this species. Analyses of segregation patterns in haploid males reveal that large genomic elements (several megabases in size) provide the substrate of genome size variation. These elements, and their segregation patterns, explain the generation of new genome size variants, the short-term evolutionary potential of genome size change in populations, and some seemingly paradoxical patterns, like an increase in genome size variation among highly inbred lines. Our study suggests that a conceptual model involving only two variables, 1) a basal genome size of the population, and 2) a vector containing information on additional elements that may increase genome size in this population (size, number, and meiotic segregation behavior), can effectively address most scenarios of short-term evolutionary change of genome size in a population.

Sex- and tissue-specific transcriptome analyses and expression profiling of olfactory-related genes in Ceracris nigricornis Walker (Orthoptera: Acrididae)

BACKGROUND

The sophisticated insect olfactory system plays an important role in recognizing external odors and enabling insects to adapt to environment. Foraging, host seeking, mating, ovipositing and other forms of chemical communication are based on olfaction, which requires the participation of multiple olfactory genes. The exclusive evolutionary trend of the olfactory system in Orthoptera insects is an excellent model for studying olfactory evolution, but limited olfaction research is available for these species. The olfactory-related genes of Ceracris nigricornis Walker (Orthoptera: Acrididae), a severe pest of bamboos, have not yet been reported.

RESULTS

We sequenced and analyzed the transcriptomes from different tissues of C. nigricornis and obtained 223.76 Gb clean data that were assembled into 43,603 unigenes with an N50 length of 2235 bp. Among the transcripts, 66.79% of unigenes were annotated. Based on annotation and tBLASTn results, 112 candidate olfactory-related genes were identified for the first time, including 20 odorant-binding proteins (OBPs), 10 chemosensory-binding proteins (CSPs), 71 odorant receptors (ORs), eight ionotropic receptors (IRs) and three sensory neuron membrane proteins (SNMPs). The fragments per kilobase per million mapped fragments (FPKM) values showed that most olfactory-related differentially expressed genes (DEGs) were enriched in the antennae, and these results were confirmed by detecting the expression of olfactory-related genes with quantitative real-time PCR (qRT-PCR). Among these antennae-enriched genes, some were sex-biased, indicating their different roles in the olfactory system of C. nigricornis.

CONCLUSIONS

This study provides the first comprehensive list and expression profiles of olfactory-related genes in C. nigricornis and a foundation for functional studies of these olfactory-related genes at the molecular level.

TERAD: Extraction of transposable element composition from RADseq data

Transposable elements (TEs) - selfish DNA sequences that can move within the genome - comprise a large proportion of the genomes of many organisms. Although low-coverage whole-genome sequencing can be used to survey TE composition, it is noneconomical for species with large quantities of DNA. Here, we utilize restriction-site associated DNA sequencing (RADSeq) as an alternative method to survey TE composition. First, we demonstrate in silico that double digest restriction-site associated DNA sequencing (ddRADseq) markers contain the same TE compositions as whole genome assemblies across arthropods. Next, we show empirically using eight Synalpheus snapping shrimp species with large genomes that TE compositions from ddRADseq and low-coverage whole-genome sequencing are comparable within and across species. Finally, we develop a new bioinformatic pipeline, TERAD, to extract TE compositions from RADseq data. Our study expands the utility of RADseq to study the repeatome, making comparative studies of genome structure for species with large genomes more tractable and affordable.

Ecological constraints associated with genome size across salamander lineages

Salamanders have some of the largest, and most variable, genome sizes among the vertebrates. Larger genomes have been associated with larger cell sizes, lower metabolic rates, and longer embryonic and larval durations in many different taxonomic groups. These life-history traits are often important for dictating fitness under different environmental conditions, suggesting that a species' genome size may have the potential to constrain its ecological distribution. We test how genome size varies with the ephemerality of larval habitat across the salamanders, predicting that species with larger genomes will be constrained to more permanent habitats that permit slower development, while species with smaller genomes will be more broadly distributed across the gradient of habitat ephemerality. We found that salamanders with larger genomes are almost exclusively associated with permanent aquatic habitats. In addition, the evolutionary transition rate between permanent and ephemeral larval habitats is much higher in salamander lineages with smaller genome sizes. These patterns suggest that genome size may act as an evolutionary constraint on the ecological habitats of salamanders, restricting those species with large genomes and slower development to habitats with permanent sources of water.

Draft Aphaenogaster genomes expand our view of ant genome size variation across climate gradients

Given the abundance, broad distribution, and diversity of roles that ants play in many ecosystems, they are an ideal group to serve as ecosystem indicators of climatic change. At present, only a few whole-genome sequences of ants are available (19 of >16,000 species), mostly from tropical and sub-tropical species. To address this limited sampling, we sequenced genomes of temperate-latitude species from the genus Aphaenogaster, a genus with important seed dispersers. In total, we sampled seven colonies of six species: Aphaenogaster ashmeadi, Aphaenogaster floridana, Aphaenogaster fulva, Aphaenogaster miamiana, Aphaenogaster picea, and Aphaenogaster rudis. The geographic ranges of these species collectively span eastern North America from southern Florida to southern Canada, which encompasses a latitudinal gradient in which many climatic variables are changing rapidly. For the six genomes, we assembled an average of 271,039 contigs into 47,337 scaffolds. The Aphaenogaster genomes displayed high levels of completeness with 96.1% to 97.6% of Hymenoptera BUSCOs completely represented, relative to currently sequenced ant genomes which ranged from 88.2% to 98.5%. Additionally, the mean genome size was 370.5 Mb, ranging from 310.3 to 429.7, which is comparable to that of other sequenced ant genomes (212.8-396.0 Mb) and flow cytometry estimates (210.7-690.4 Mb). In an analysis of currently sequenced ant genomes and the new Aphaenogaster sequences, we found that after controlling for both spatial autocorrelation and phylogenetics ant genome size was marginally correlated with sample site climate similarity. Of all examined climate variables, minimum temperature, and annual precipitation had the strongest correlations with genome size, with ants from locations with colder minimum temperatures and higher levels of precipitation having larger genomes. These results suggest that climate extremes could be a selective force acting on ant genomes and point to the need for more extensive sequencing of ant genomes.

Comparative analysis of mitochondrial genomes of the superfamily Grylloidea (Insecta, Orthoptera) reveals phylogenetic distribution of gene rearrangements

To further characterize mitochondrial genome (mitogenome) features of the superfamily Grylloidea (Insecta, Orthoptera), mitogenomes of Cacoplistes rogenhoferi and Meloimorpha japonica representing the family Mogoplistidae and three Ornebius species of Phalangopsidae were sequenced. A repeat-containing control region (CR) and 37 genes were present in these mitogenomes. Unusual start codons (TCG, CCG, and CTG) of cox1 and, in Ornebius, a partial stop codon (T) of nad1 followed by a 15-17-bp intergenic spacer were proposed based on transcript information and sequence alignments. The mitogenome-based phylogenetic trees suggest strongly the familial relationships as (((Phalangopsidae + Gryllidae) + Trigonidiidae) + Mogoplistidae). The exclusive occurrence of the trnE-trnS1-trnN rearrangement in Phalangopsidae, Gryllidae, and Trigonidiidae is suggestive of its appearance in the common ancestor of these families after the separation of Mogoplistidae. The trnV transposition in O. bimaculatus and formerly sequenced Trigonidium sjostedti (Trigonidiidae) indicates a potential consequence of parallel evolution. This study offers novel insights into mitogenome evolution, especially gene rearrangements, of Grylloidea.

Genome size estimation of brackishwater fishes and penaeid shrimps by flow cytometry

Flow cytometry was used for estimating the genome size of five brackishwater finfish and four shrimp species. The genome size for Lutjanus argentimaculatus was 0.95 ± 0.10 and 0.79 ± 0.01 pg for Scatophagus argus. The genome sizes for Chanos chanos (0.72 ± 0.01 pg), Etroplus suratensis (1.71 ± 0.16 pg) and Liza macrolepis (0.87 ± 0.02 pg) which are important aquaculture species are reported for the first time in this study. The phylogenetic tree constructed using sixty-seven sequence accessions of cytochrome c oxidase subunit 1 (COI) gene of Lates calcarifer revealed two separate clades. The Indian Lates calcarifer species with estimated genome size of 0.44 ± 0.02 pg belonged to a clade different than that of South East Asia and Australia reported to have larger genome size. The genome size for the four major species of genus Penaeus (Penaeus monodon, Penaeus indicus, Penaeus vannamei and Penaeus japonicus) were found in similar range. The genome size of female shrimps ranged from 2.91 ± 0.03 pg (P. monodon) to 2.14 ± 0.02 pg (P. japonicus). In male shrimps, the genome size ranged from 2.86 ± 0.06 pg (P. monodon) to 2.19 ± 0.02 pg (P. indicus). Significant difference was observed in the genome size between male and female shrimp of all species except in P. monodon. The highest relative difference of 12.78% was observed in the genome size between the either sex in P. indicus. The interspecific relative difference of 30.59% in genome size was highest between the male shrimps of P. monodon and P. indicus and 35.98% between the female shrimps of P. monodon and P. japonicus. The stored gills and pleopod tissues could be successfully used up to 3 weeks to estimate the genome size in shrimps.

Genome size variation in deep-sea amphipods

Genome size varies considerably across taxa, and extensive research effort has gone into understanding whether variation can be explained by differences in key ecological and life-history traits among species. The extreme environmental conditions that characterize the deep sea have been hypothesized to promote large genome sizes in eukaryotes. Here we test this supposition by examining genome sizes among 13 species of deep-sea amphipods from the Mariana, Kermadec and New Hebrides trenches. Genome sizes were estimated using flow cytometry and found to vary nine-fold, ranging from 4.06 pg (4.04 Gb) in Paralicella caperesca to 34.79 pg (34.02 Gb) in Alicella gigantea. Phylogenetic independent contrast analysis identified a relationship between genome size and maximum body size, though this was largely driven by those species that display size gigantism. There was a distinct shift in the genome size trait diversification rate in the supergiant amphipod A. gigantea relative to the rest of the group. The variation in genome size observed is striking and argues against genome size being driven by a common evolutionary history, ecological niche and life-history strategy in deep-sea amphipods.

Genome size in arthropods; different roles of phylogeny, habitat and life history in insects and crustaceans

Despite the major role of genome size for physiology, ecology, and evolution, there is still mixed evidence with regard to proximate and ultimate drivers. The main causes of large genome size are proliferation of noncoding elements and/or duplication events. The relative role and interplay between these proximate causes and the evolutionary patterns shaped by phylogeny, life history traits or environment are largely unknown for the arthropods. Genome size shows a tremendous variability in this group, and it has a major impact on a range of fitness-related parameters such as growth, metabolism, life history traits, and for many species also body size. In this study, we compared genome size in two major arthropod groups, insects and crustaceans, and related this to phylogenetic patterns and parameters affecting ambient temperature (latitude, depth, or altitude), insect developmental mode, as well as crustacean body size and habitat, for species where data were available. For the insects, the genome size is clearly phylogeny-dependent, reflecting primarily their life history and mode of development, while for crustaceans there was a weaker association between genome size and phylogeny, suggesting life cycle strategies and habitat as more important determinants. Maximum observed latitude and depth, and their combined effect, showed positive, and possibly phylogenetic independent, correlations with genome size for crustaceans. This study illustrate the striking difference in genome sizes both between and within these two major groups of arthropods, and that while living in the cold with low developmental rates may promote large genomes in marine crustaceans, there is a multitude of proximate and ultimate drivers of genome size.

Nuclear DNA content correlates with depth, body size, and diversification rate in amphipod crustaceans from ancient Lake Baikal, Russia

Lake Baikal in Russia is a large, ancient lake that has been the site of a major radiation of amphipod crustaceans. Nearly 400 named species are known in this single lake, and it is thought that many more await description. The size and depth of Lake Baikal, in particular, may have contributed to the radiation of endemic amphipods by providing a large number of microhabitats for species to invade and subsequently experience reproductive isolation. Here we investigate the possibility that large-scale genomic changes have also accompanied diversification in these crustaceans. Specifically, we report genome size estimates for 36 species of Baikal amphipods, and examine the relationship between genome size, body size, and the maximum depths at which the amphipods are found in the lake. Genome sizes ranged nearly 8-fold in this sample of amphipod species, from 2.15 to 16.63 pg, and there were significant, positive, phylogenetically corrected relationships between genome size, body size, maximum depth, and diversification rate among these species. Our results suggest that major genomic changes, including transposable element proliferation, have accompanied speciation that was driven by selection for differences in body size and habitat preference in Lake Baikal amphipods.

Evolutionary characterization of Ty3/gypsy-like LTR retrotransposons in the parasitic cestode Echinococcus granulosus

Cyclophyllidean cestodes including Echinococcus granulosus have a smaller genome and show characteristics such as loss of the gut, a segmented body plan, and accelerated growth rate in hosts compared with other tissue-invading helminths. In an effort to address the molecular mechanism relevant to genome shrinkage, the evolutionary status of long-terminal-repeat (LTR) retrotransposons, which are known as the most potent genomic modulators, was investigated in the E. granulosus draft genome. A majority of the E. granulosus LTR retrotransposons were classified into a novel characteristic clade, named Saci-2, of the Ty3/gypsy family, while the remaining elements belonged to the CsRn1 clade of identical family. Their nucleotide sequences were heavily corrupted by frequent base substitutions and segmental losses. The ceased mobile activity of the major retrotransposons and the following intrinsic DNA loss in their inactive progenies might have contributed to decrease in genome size. Apart from the degenerate copies, a gag gene originating from a CsRn1-like element exhibited substantial evidences suggesting its domestication including a preserved coding profile and transcriptional activity, the presence of syntenic orthologues in cestodes, and selective pressure acting on the gene. To my knowledge, the endogenized gag gene is reported for the first time in invertebrates, though its biological function remains elusive.

Para-allopatry in hybridizing fire-bellied toads (Bombina bombina and B. variegata): Inference from transcriptome-wide coalescence analyses

Ancient origins, profound ecological divergence, and extensive hybridization make the fire‐bellied toads Bombina bombina and B. variegata (Anura: Bombinatoridae) an intriguing test case of ecological speciation. Previous modeling has proposed that the narrow Bombina hybrid zones represent strong barriers to neutral introgression. We test this prediction by inferring the rate of gene exchange between pure populations on either side of the intensively studied Kraków transect. We developed a method to extract high confidence sets of orthologous genes from de novo transcriptome assemblies, fitted a range of divergence models to these data and assessed their relative support with analytic likelihood calculations. There was clear evidence for postdivergence gene flow, but, as expected, no perceptible signal of recent introgression via the nearby hybrid zone. The analysis of two additional Bombina taxa (B. v. scabra and B. orientalis) validated our parameter estimates against a larger set of prior expectations. Despite substantial cumulative introgression over millions of years, adaptive divergence of the hybridizing taxa is essentially unaffected by their lack of reproductive isolation. Extended distribution ranges also buffer them against small‐scale environmental perturbations that have been shown to reverse the speciation process in other, more recent ecotypes.

Transposons, Genome Size, and Evolutionary Insights in Animals

The relationship between genome size and the percentage of transposons in 161 animal species evidenced that variations in genome size are linked to the amplification or the contraction of transposable elements. The activity of transposable elements could represent a response to environmental stressors. Indeed, although with different trends in protostomes and deuterostomes, comprehensive changes in genome size were recorded in concomitance with particular periods of evolutionary history or adaptations to specific environments. During evolution, genome size and the presence of transposable elements have influenced structural and functional parameters of genomes and cells. Changes of these parameters have had an impact on morphological and functional characteristics of the organism on which natural selection directly acts. Therefore, the current situation represents a balance between insertion and amplification of transposons and the mechanisms responsible for their deletion or for decreasing their activity. Among the latter, methylation and the silencing action of small RNAs likely represent the most frequent mechanisms.

Information in the Biosphere: Biological and Digital Worlds

Evolution has transformed life through key innovations in information storage and replication, including RNA, DNA, multicellularity, and culture and language. We argue that the carbon-based biosphere has generated a cognitive system (humans) capable of creating technology that will result in a comparable evolutionary transition. Digital information has reached a similar magnitude to information in the biosphere. It increases exponentially, exhibits high-fidelity replication, evolves through differential fitness, is expressed through artificial intelligence (AI), and has facility for virtually limitless recombination. Like previous evolutionary transitions, the potential symbiosis between biological and digital information will reach a critical point where these codes could compete via natural selection. Alternatively, this fusion could create a higher-level superorganism employing a low-conflict division of labor in performing informational tasks.

What's in a genome? The C-value enigma and the evolution of eukaryotic genome content

Some notable exceptions aside, eukaryotic genomes are distinguished from those of Bacteria and Archaea in a number of ways, including chromosome structure and number, repetitive DNA content, and the presence of introns in protein-coding regions. One of the most notable differences between eukaryotic and prokaryotic genomes is in size. Unlike their prokaryotic counterparts, eukaryotes exhibit enormous (more than 60,000-fold) variability in genome size which is not explained by differences in gene number. Genome size is known to correlate with cell size and division rate, and by extension with numerous organism-level traits such as metabolism, developmental rate or body size. Less well described are the relationships between genome size and other properties of the genome, such as gene content, transposable element content, base pair composition and related features. The rapid expansion of 'complete' genome sequencing projects has, for the first time, made it possible to examine these relationships across a wide range of eukaryotes in order to shed new light on the causes and correlates of genome size diversity. This study presents the results of phylogenetically informed comparisons of genome data for more than 500 species of eukaryotes. Several relationships are described between genome size and other genomic parameters, and some recommendations are presented for how these insights can be extended even more broadly in the future.

A population genomic scan in Chorthippus grasshoppers unveils previously unknown phenotypic divergence

Understanding the genetics of speciation and the processes that drive it is a central goal of evolutionary biology. Grasshoppers of the Chorthippus species group differ strongly in calling song (and corresponding female preferences) but are exceedingly similar in other characteristics such as morphology. Here, we performed a population genomic scan on three Chorthippus species (Chorthippus biguttulus, C. mollis and C. brunneus) to gain insight into the genes and processes involved in divergence and speciation in this group. Using an RNA-seq approach, we examined functional variation between the species by calling SNPs for each of the three species pairs and using FST -based approaches to identify outliers. We found approximately 1% of SNPs in each comparison to be outliers. Between 37% and 40% of these outliers were nonsynonymous SNPs (as opposed to a global level of 17%) indicating that we recovered loci under selection. Among the outliers were several genes that may be involved in song production and hearing as well as genes involved in other traits such as food preferences and metabolism. Differences in food preferences between species were confirmed with a behavioural experiment. This indicates that multiple phenotypic differences implicating multiple evolutionary processes (sexual selection and natural selection) are present between the species.

Comparative analysis of transposable elements highlights mobilome diversity and evolution in vertebrates

Transposable elements (TEs) are major components of vertebrate genomes, with major roles in genome architecture and evolution. In order to characterize both common patterns and lineage-specific differences in TE content and TE evolution, we have compared the mobilomes of 23 vertebrate genomes, including 10 actinopterygian fish, 11 sarcopterygians, and 2 nonbony vertebrates. We found important variations in TE content (from 6% in the pufferfish tetraodon to 55% in zebrafish), with a more important relative contribution of TEs to genome size in fish than in mammals. Some TE superfamilies were found to be widespread in vertebrates, but most elements showed a more patchy distribution, indicative of multiple events of loss or gain. Interestingly, loss of major TE families was observed during the evolution of the sarcopterygian lineage, with a particularly strong reduction in TE diversity in birds and mammals. Phylogenetic trends in TE composition and activity were detected: Teleost fish genomes are dominated by DNA transposons and contain few ancient TE copies, while mammalian genomes have been predominantly shaped by nonlong terminal repeat retrotransposons, along with the persistence of older sequences. Differences were also found within lineages: The medaka fish genome underwent more recent TE amplification than the related platyfish, as observed for LINE retrotransposons in the mouse compared with the human genome. This study allows the identification of putative cases of horizontal transfer of TEs, and to tentatively infer the composition of the ancestral vertebrate mobilome. Taken together, the results obtained highlight the importance of TEs in the structure and evolution of vertebrate genomes, and demonstrate their major impact on genome diversity both between and within lineages.

The Genome 10K Project: a way forward

The Genome 10K Project was established in 2009 by a consortium of biologists and genome scientists determined to facilitate the sequencing and analysis of the complete genomes of 10,000 vertebrate species. Since then the number of selected and initiated species has risen from ∼26 to 277 sequenced or ongoing with funding, an approximately tenfold increase in five years. Here we summarize the advances and commitments that have occurred by mid-2014 and outline the achievements and present challenges of reaching the 10,000-species goal. We summarize the status of known vertebrate genome projects, recommend standards for pronouncing a genome as sequenced or completed, and provide our present and future vision of the landscape of Genome 10K. The endeavor is ambitious, bold, expensive, and uncertain, but together the Genome 10K Consortium of Scientists and the worldwide genomics community are moving toward their goal of delivering to the coming generation the gift of genome empowerment for many vertebrate species.

A generalized mechanistic codon model

Models of codon evolution have attracted particular interest because of their unique capabilities to detect selection forces and their high fit when applied to sequence evolution. We described here a novel approach for modeling codon evolution, which is based on Kronecker product of matrices. The 61 × 61 codon substitution rate matrix is created using Kronecker product of three 4 × 4 nucleotide substitution matrices, the equilibrium frequency of codons, and the selection rate parameter. The entities of the nucleotide substitution matrices and selection rate are considered as parameters of the model, which are optimized by maximum likelihood. Our fully mechanistic model allows the instantaneous substitution matrix between codons to be fully estimated with only 19 parameters instead of 3,721, by using the biological interdependence existing between positions within codons. We illustrate the properties of our models using computer simulations and assessed its relevance by comparing the AICc measures of our model and other models of codon evolution on simulations and a large range of empirical data sets. We show that our model fits most biological data better compared with the current codon models. Furthermore, the parameters in our model can be interpreted in a similar way as the exchangeability rates found in empirical codon models.