Transposable Elements, Epigenetics, and Genome Evolution

From the cauldron of conflict: Endogenous gene regulation by piRNA and other modes of adaptation enabled by selfish transposable elements

Transposable elements (TEs) provide a prime example of genetic conflict because they can proliferate in genomes and populations even if they harm the host. However, numerous studies have shown that TEs, though typically harmful, can also provide fuel for adaptation. This is because they code functional sequences that can be useful for the host in which they reside. In this review, I summarize the "how" and "why" of adaptation enabled by the genetic conflict between TEs and hosts. In addition, focusing on mechanisms of TE control by small piwi-interacting RNAs (piRNAs), I highlight an indirect form of adaptation enabled by conflict. In this case, mechanisms of host defense that regulate TEs have been redeployed for endogenous gene regulation. I propose that the genetic conflict released by meiosis in early eukaryotes may have been important because, among other reasons, it spurred evolutionary innovation on multiple interwoven trajectories - on the part of hosts and also embedded genetic parasites. This form of evolution may function as a complexity generating engine that was a critical player in eukaryotic evolution.

Assessing genotoxic effects of plastic leachates in Drosophila melanogaster

Plastic polymers were largely added with chemical substances to be utilized in the items and product manufacturing. The leachability of these substances is a matter of concern given the wide amount of plastic waste, particularly in terrestrial environments, where soil represents a sink for these novel contaminants and a possible pathway of human health risk. In this study, we integrated genetic, molecular, and behavioral approaches to comparatively evaluate toxicological effects of plastic leachates, virgin and oxodegradable polypropylene (PP) and polyethylene (PE), in Drosophila melanogaster, a novel in vivo model organism for environmental monitoring studies and (eco)toxicological research. The results of this study revealed that while conventional toxicological endpoints such as developmental times and longevity remain largely unaffected, exposure to plastic leachates induces chromosomal abnormalities and transposable element (TE) activation in neural tissues. The combined effects of DNA damage and TE mobilization contribute to genome instability and increase the likelihood of LOH events, thus potentiating tumor growth and metastatic behavior ofRasV12 clones. Collectively, these findings indicate that plastic leachates exert genotoxic effects in Drosophila thus highlighting potential risks associated with leachate-related plastic pollution and their implications for ecosystems and human health.

Symmetric and asymmetric DNA N6-adenine methylation regulates different biological responses in Mucorales

DNA N6-adenine methylation (6mA) has recently gained importance as an epigenetic modification in eukaryotes. Its function in lineages with high levels, such as early-diverging fungi (EDF), is of particular interest. Here, we investigated the biological significance and evolutionary implications of 6mA in EDF, which exhibit divergent evolutionary patterns in 6mA usage. The analysis of two Mucorales species displaying extreme 6mA usage reveals that species with high 6mA levels show symmetric methylation enriched in highly expressed genes. In contrast, species with low 6mA levels show mostly asymmetric 6mA. Interestingly, transcriptomic regulation throughout development and in response to environmental cues is associated with changes in the 6mA landscape. Furthermore, we identify an EDF-specific methyltransferase, likely originated from endosymbiotic bacteria, as responsible for asymmetric methylation, while an MTA-70 methylation complex performs symmetric methylation. The distinct phenotypes observed in the corresponding mutants reinforced the critical role of both types of 6mA in EDF.

The regulatory potential of transposable elements in maize

Since their initial discovery in maize, transposable elements (TEs) have emerged as being integral to the evolution of maize, accounting for 80% of its genome. However, the repetitive nature of TEs has hindered our understanding of their regulatory potential. Here, we demonstrate that long-read chromatin fiber sequencing (Fiber-seq) permits the comprehensive annotation of the regulatory potential of maize TEs. We uncover that only 94 LTR retrotransposons contain the functional epigenetic architecture required for mobilization within maize leaves. This epigenetic architecture degenerates with evolutionary age, resulting in solo TE enhancers being preferentially marked by simultaneous hyper-CpG methylation and chromatin accessibility, an architecture markedly divergent from canonical enhancers. We find that TEs shape maize gene regulation by creating novel promoters within the TE itself as well as through TE-mediated gene amplification. Lastly, we uncover a pervasive epigenetic code directing TEs to specific loci, including that locus that sparked McClintock's discovery of TEs.

Evolutionary adaptation of an HP1-protein chromodomain integrates chromatin and DNA sequence signals

Members of the diverse heterochromatin protein 1 (HP1) family play crucial roles in heterochromatin formation and maintenance. Despite the similar affinities of their chromodomains for di- and tri-methylated histone H3 lysine 9 (H3K9me2/3), different HP1 proteins exhibit distinct chromatin-binding patterns, likely due to interactions with various specificity factors. Previously, we showed that the chromatin-binding pattern of the HP1 protein Rhino, a crucial factor of the Drosophila PIWI-interacting RNA (piRNA) pathway, is largely defined by a DNA sequence-specific C2H2 zinc finger protein named Kipferl (Baumgartner et al., 2022). Here, we elucidate the molecular basis of the interaction between Rhino and its guidance factor Kipferl. Through phylogenetic analyses, structure prediction, and in vivo genetics, we identify a single amino acid change within Rhino's chromodomain, G31D, that does not affect H3K9me2/3 binding but disrupts the interaction between Rhino and Kipferl. Flies carrying the rhinoG31D mutation phenocopy kipferl mutant flies, with Rhino redistributing from piRNA clusters to satellite repeats, causing pronounced changes in the ovarian piRNA profile of rhinoG31D flies. Thus, Rhino's chromodomain functions as a dual-specificity module, facilitating interactions with both a histone mark and a DNA-binding protein.

Expansions and contractions of repetitive DNA elements reveal contrasting evolutionary responses to the polyploid genome shock hypothesis in Brachypodium model grasses

Brachypodium grass species have been selected as model plants for functional genomics of grass crops, and to elucidate the origins of allopolyploidy and perenniality in monocots, due to their small genome sizes and feasibility of cultivation. However, genome sizes differ greatly between diploid or polyploid Brachypodium lineages. We have used genome skimming sequencing data to uncover the composition, abundance, and phylogenetic value of repetitive elements in 44 representatives of the major Brachypodium lineages and cytotypes. We also aimed to test the possible mechanisms and consequences of the "polyploid genome shock hypothesis" (PGSH) under three different evolutionary scenarios of variation in repeats and genome sizes of Brachypodium allopolyploids. Our data indicated that the proportion of the genome covered by the repeatome in the Brachypodium species showed a 3.3-fold difference between the highest content of B. mexicanum-4x (67.97%) and the lowest of B. stacei-2x (20.77%), and that changes in the sizes of their genomes were a consequence of gains or losses in their repeat elements. LTR-Retand and Tekay retrotransposons were the most frequent repeat elements in the Brachypodium genomes, while Ogre retrotransposons were found exclusively in B. mexicanum. The repeatome phylogenetic network showed a high topological congruence with plastome and nuclear rDNA and transcriptome trees, differentiating the ancestral outcore lineages from the recently evolved core-perennial lineages. The 5S rDNA graph topologies had a strong match with the ploidy levels and nature of the subgenomes of the Brachypodium polyploids. The core-perennial B. sylvaticum presents a large repeatome and characteristics of a potential post-polyploid diploidized origin. Our study evidenced that expansions and contractions in the repeatome were responsible for the three contrasting responses to the PGSH. The exacerbated genome expansion of the ancestral allotetraploid B. mexicanum was a consequence of chromosome-wide proliferation of TEs and not of WGD, the additive repeatome pattern of young allotetraploid B. hybridum of stabilized post-WGD genome evolution, and the genomecontraction of recent core-perennials polyploids (B. pinnatum, B. phoenicoides) of repeat losses through recombination of these highly hybridizing lineages. Our analyses have contributed to unraveling the evolution of the repeatome and the genome size variation in model Brachypodium grasses.

Harbinger transposon insertion in ethylene signaling gene leads to emergence of new sexual forms in cucurbits

In flowering plants, the predominant sexual morph is hermaphroditism, and the emergence of unisexuality is poorly understood. Using Cucumis melo (melon) as a model system, we explore the mechanisms driving sexual forms. We identify a spontaneous mutant exhibiting a transition from bisexual to unisexual male flower, and identify the causal mutation as a Harbinger transposon impairing the expression of Ethylene Insensitive 2 (CmEIN2) gene. Genetics and transcriptomic analysis reveal a dual role of CmEIN2 in both sex determination and fruit shape formation. Upon expression of CmACS11, EIN2 is recruited to repress the expression of the carpel inhibitor, CmWIP1. Subsequently, EIN2 is recruited to mediate stamina inhibition. Following the sex determination phase, EIN2 promotes fruit shape elongation. Genome-wide analysis reveals that Harbinger transposon mobilization is triggered by environmental cues, and integrates preferentially in active chromatin, particularly within promoter regions. Characterization of a large collection of melon germplasm points to active transpositions in the wild, compared to cultivated accessions. Our study underscores the association between chromatin dynamics and the temporal aspects of mobile genetic element insertions, providing valuable insights into plant adaptation and crop genome evolution.

Mutation protocols share with sexual reproduction the physiological role of producing genetic variation within 'constraints that deconstrain'

Because the universe of possible DNA sequences is inconceivably vast, organisms have evolved mechanisms for exploring DNA sequence space while substantially reducing the hazard that would otherwise accrue to any process of random, accidental mutation. One such mechanism is meiotic recombination. Although sexual reproduction imposes a seemingly paradoxical 50% cost to fitness, sex evidently prevails because this cost is outweighed by the advantage of equipping offspring with genetic variation to accommodate environmental vicissitudes. The potential adaptive utility of additional mechanisms for producing genetic variation has long been obscured by a presumption that the vast majority of mutations are deleterious. Perhaps surprisingly, the probability for adaptive variation can be increased by several mechanisms that generate mutations abundantly. Such mechanisms, here called 'mutation protocols', implement implicit 'constraints that deconstrain'. Like meiotic recombination, they produce genetic variation in forms that minimize potential for harm while providing a reasonably high probability for benefit. One example is replication slippage of simple sequence repeats (SSRs); this process yields abundant, reversible mutations, typically with small quantitative effect on phenotype. This enables SSRs to function as adjustable 'tuning knobs'. There exists a clear pathway for SSRs to be shaped through indirect selection favouring their implicit tuning-knob protocol. Several other molecular mechanisms comprise probable components of additional mutation protocols. Biologists might plausibly regard such mechanisms of mutation not primarily as sources of deleterious genetic mistakes but also as potentially adaptive processes for 'exploring' DNA sequence space.

ACMGA: a reference-free multiple-genome alignment pipeline for plant species

BACKGROUND

The short-read whole-genome sequencing (WGS) approach has been widely applied to investigate the genomic variation in the natural populations of many plant species. With the rapid advancements in long-read sequencing and genome assembly technologies, high-quality genome sequences are available for a group of varieties for many plant species. These genome sequences are expected to help researchers comprehensively investigate any type of genomic variants that are missed by the WGS technology. However, multiple genome alignment (MGA) tools designed by the human genome research community might be unsuitable for plant genomes.

RESULTS

To fill this gap, we developed the AnchorWave-Cactus Multiple Genome Alignment (ACMGA) pipeline, which improved the alignment of repeat elements and could identify long (> 50 bp) deletions or insertions (INDELs). We conducted MGA using ACMGA and Cactus for 8 Arabidopsis (Arabidopsis thaliana) and 26 Maize (Zea mays) de novo assembled genome sequences and compared them with the previously published short-read variant calling results. MGA identified more single nucleotide variants (SNVs) and long INDELs than did previously published WGS variant callings. Additionally, ACMGA detected significantly more SNVs and long INDELs in repetitive regions and the whole genome than did Cactus. Compared with the results of Cactus, the results of ACMGA were more similar to the previously published variants called using short-read. These two MGA pipelines identified numerous multi-allelic variants that were missed by the WGS variant calling pipeline.

CONCLUSIONS

Aligning de novo assembled genome sequences could identify more SNVs and INDELs than mapping short-read. ACMGA combines the advantages of AnchorWave and Cactus and offers a practical solution for plant MGA by integrating global alignment, a 2-piece-affine-gap cost strategy, and the progressive MGA algorithm.

The rapidly evolving X-linked MIR-506 family fine-tunes spermatogenesis to enhance sperm competition

Despite rapid evolution across eutherian mammals, the X-linked MIR-506 family miRNAs are located in a region flanked by two highly conserved protein-coding genes (SLITRK2 and FMR1) on the X chromosome. Intriguingly, these miRNAs are predominantly expressed in the testis, suggesting a potential role in spermatogenesis and male fertility. Here, we report that the X-linked MIR-506 family miRNAs were derived from the MER91C DNA transposons. Selective inactivation of individual miRNAs or clusters caused no discernible defects, but simultaneous ablation of five clusters containing 19 members of the MIR-506 family led to reduced male fertility in mice. Despite normal sperm counts, motility, and morphology, the KO sperm were less competitive than wild-type sperm when subjected to a polyandrous mating scheme. Transcriptomic and bioinformatic analyses revealed that these X-linked MIR-506 family miRNAs, in addition to targeting a set of conserved genes, have more targets that are critical for spermatogenesis and embryonic development during evolution. Our data suggest that the MIR-506 family miRNAs function to enhance sperm competitiveness and reproductive fitness of the male by finetuning gene expression during spermatogenesis.

Adaptive expansion of ERVK solo-LTRs is associated with Passeriformes speciation events

Endogenous retroviruses (ERVs) are ancient retroviral remnants integrated in host genomes, and commonly deleted through unequal homologous recombination, leaving solitary long terminal repeats (solo-LTRs). This study, analysing the genomes of 362 bird species and their reptilian and mammalian outgroups, reveals an unusually higher level of solo-LTRs formation in birds, indicating evolutionary forces might have purged ERVs during evolution. Strikingly in the order Passeriformes, and especially the parvorder Passerida, endogenous retrovirus K (ERVK) solo-LTRs showed bursts of formation and recurrent accumulations coinciding with speciation events over past 22 million years. Moreover, our results indicate that the ongoing expansion of ERVK solo-LTRs in these bird species, marked by high transcriptional activity of ERVK retroviral genes in reproductive organs, caused variation of solo-LTRs between individual zebra finches. We experimentally demonstrated that cis-regulatory activity of recently evolved ERVK solo-LTRs may significantly increase the expression level of ITGA2 in the brain of zebra finches compared to chickens. These findings suggest that ERVK solo-LTRs expansion may introduce novel genomic sequences acting as cis-regulatory elements and contribute to adaptive evolution. Overall, our results underscore that the residual sequences of ancient retroviruses could influence the adaptive diversification of species by regulating host gene expression.

Morc1 reestablishes H3K9me3 heterochromatin on piRNA-targeted transposons in gonocytes

To maintain fertility, male mice re-repress transposable elements (TEs) that were de-silenced in the early gonocytes before their differentiation into spermatogonia. However, the mechanism of TE silencing re-establishment remains unknown. Here, we found that the DNA-binding protein Morc1, in cooperation with the methyltransferase SetDB1, deposits the repressive histone mark H3K9me3 on a large fraction of activated TEs, leading to heterochromatin. Morc1 also triggers DNA methylation, but TEs targeted by Morc1-driven DNA methylation only slightly overlapped with those repressed by Morc1/SetDB1-dependent heterochromatin formation, suggesting that Morc1 silences TEs in two different manners. In contrast, TEs regulated by Morc1 and Miwi2, the nuclear PIWI-family protein, almost overlapped. Miwi2 binds to PIWI-interacting RNAs (piRNAs) that base-pair with TE mRNAs via sequence complementarity, while Morc1 DNA binding is not sequence specific, suggesting that Miwi2 selects its targets, and then, Morc1 acts to repress them with cofactors. A high-ordered mechanism of TE repression in gonocytes has been identified.

New whole-genome alignment tools are needed for tapping into plant diversity

Genome alignment is one of the most foundational methods for genome sequence studies. With rapid advances in sequencing and assembly technologies, these newly assembled genomes present challenges for alignment tools to meet the increased complexity and scale. Plant genome alignment is technologically challenging because of frequent whole-genome duplications (WGDs) as well as chromosome rearrangements and fractionation, high nucleotide diversity, widespread structural variation, and high transposable element (TE) activity causing large proportions of repeat elements. We summarize classical pairwise and multiple genome alignment (MGA) methods, and highlight techniques that are widely used or are being developed by the plant research community. We also outline the remaining challenges for precise genome alignment and the interpretation of alignment results in plants.

The Rapidly Evolving X-linked miR-506 Family Finetunes Spermatogenesis to Enhance Sperm Competition

Despite rapid evolution across eutherian mammals, the X-linked miR-506 family miRNAs are located in a region flanked by two highly conserved protein-coding genes (Slitrk2 and Fmr1) on the X chromosome. Intriguingly, these miRNAs are predominantly expressed in the testis, suggesting a potential role in spermatogenesis and male fertility. Here, we report that the X-linked miR-506 family miRNAs were derived from the MER91C DNA transposons. Selective inactivation of individual miRNAs or clusters caused no discernable defects, but simultaneous ablation of five clusters containing nineteen members of the miR-506 family led to reduced male fertility in mice. Despite normal sperm counts, motility and morphology, the KO sperm were less competitive than wild-type sperm when subjected to a polyandrous mating scheme. Transcriptomic and bioinformatic analyses revealed that these X-linked miR-506 family miRNAs, in addition to targeting a set of conserved genes, have more targets that are critical for spermatogenesis and embryonic development during evolution. Our data suggest that the miR-506 family miRNAs function to enhance sperm competitiveness and reproductive fitness of the male by finetuning gene expression during spermatogenesis.

Transposable elements: multifunctional players in the plant genome

Transposable elements (TEs) are indispensable components of eukaryotic genomes that play diverse roles in gene regulation, recombination, and environmental adaptation. Their ability to mobilize within the genome leads to gene expression and DNA structure changes. TEs serve as valuable markers for genetic and evolutionary studies and facilitate genetic mapping and phylogenetic analysis. They also provide insight into how organisms adapt to a changing environment by promoting gene rearrangements that lead to new gene combinations. These repetitive sequences significantly impact genome structure, function and evolution. This review takes a comprehensive look at TEs and their applications in biotechnology, particularly in the context of plant biology, where they are now considered "genomic gold" due to their extensive functionalities. The article addresses various aspects of TEs in plant development, including their structure, epigenetic regulation, evolutionary patterns, and their use in gene editing and plant molecular markers. The goal is to systematically understand TEs and shed light on their diverse roles in plant biology.

Deciphering the β-carotene hyperaccumulation in Dunaliella by the comprehensive analysis of Dunaliella salina and Dunaliella tertiolecta under high light conditions

The green microalga Dunaliella salina hyperaccumulates β-carotene in the chloroplast, which turns its cells orange. This does not occur in the sister species Dunaliella tertiolecta. However, the molecular mechanisms of β-carotene hyperaccumulation were still unclear. Here, we discovered the reasons for β-carotene hyperaccumulation by comparing the morphology, physiology, genome, and transcriptome between the carotenogenic D. salina and the noncarotenogenic D. tertiolecta after transfer to high light. The differences in photosynthetic capacity, cell growth, and the concentration of stored carbon suggest that these species regulate the supply and utilization of carbon differently. The number of β-carotene-containing plastid lipid globules increased in both species, but much faster and to a greater extent in D. salina than in D. tertiolecta. Consistent with the accumulation of plastid lipid globules, the expression of the methyl-erythritol-phosphate and carotenoid biosynthetic pathways increased only in D. salina, which explains the de novo synthesis of β-carotene. In D. salina, the concomitantly upregulated expression of the carotene globule proteins suggests that hyperaccumulation of β-carotene also requires a simultaneous increase in its sink capacity. Based on genomic analysis, we propose that D. salina has genetic advantages for routing carbon from growth to carotenoid metabolism.

A review of plant epigenetics through the lens of almond

While genomes were originally seen as static entities that stably held and organized genetic information, recent advances in sequencing have uncovered the dynamic nature of the genome. New conceptualizations of the genome include complex relationships between the environment and gene expression that must be maintained, regulated, and sometimes even transmitted over generations. The discovery of epigenetic mechanisms has allowed researchers to understand how traits like phenology, plasticity, and fitness can be altered without changing the underlying deoxyribonucleic acid sequence. While many discoveries were first made in animal systems, plants provide a particularly complex set of epigenetic mechanisms due to unique aspects of their biology and interactions with human selective breeding and cultivation. In the plant kingdom, annual plants have received the most attention; however, perennial plants endure and respond to their environment and human management in distinct ways. Perennials include crops such as almond, for which epigenetic effects have long been linked to phenomena and even considered relevant for breeding. Recent discoveries have elucidated epigenetic phenomena that influence traits such as dormancy and self-compatibility, as well as disorders like noninfectious bud failure, which are known to be triggered by the environment and influenced by inherent aspects of the plant. Thus, epigenetics represents fertile ground to further understand almond biology and production and optimize its breeding. Here, we provide our current understanding of epigenetic regulation in plants and use almond as an example of how advances in epigenetics research can be used to understand biological fitness and agricultural performance in crop plants.

Enhancers dysfunction in the 3D genome of cancer cells

Eukaryotic genomes are spatially organized inside the cell nucleus, forming a threedimensional (3D) architecture that allows for spatial separation of nuclear processes and for controlled expression of genes required for cell identity specification and tissue homeostasis. Hence, it is of no surprise that mis-regulation of genome architecture through rearrangements of the linear genome sequence or epigenetic perturbations are often linked to aberrant gene expression programs in tumor cells. Increasing research efforts have shed light into the causes and consequences of alterations of 3D genome organization. In this review, we summarize the current knowledge on how 3D genome architecture is dysregulated in cancer, with a focus on enhancer highjacking events and their contribution to tumorigenesis. Studying the functional effects of genome architecture perturbations on gene expression in cancer offers a unique opportunity for a deeper understanding of tumor biology and sets the basis for the discovery of novel therapeutic targets.

Chromosomal Radiation: A model to explain karyotypic diversity in cryptic species

We present a concept that explains the pattern of occurrence of widely distributed organisms with large chromosomal diversity, large or small molecular divergence, and the insufficiency or absence of morphological identity. Our model is based on cytogenetic studies associated with molecular and biological data and can be applied to any lineage of sister species, chronospecies, or cryptic species. Through the evaluation of the karyotypic macrostructure, as the physical location of genes e satellites DNAs, in addition to phylogenetic reconstructions from mitochondrial and nuclear genes, per example, we have observed morphologically indistinguishable individuals presenting different locally fixed karyomorphs with phylogeographic discontinuity. The biological process behind this pattern is seen in many groups of cryptic species, in which variation lies mainly in the organization of their genomes but not necessarily in the ecosystems they inhabit or in their external morphology. It's similar to the processes behind other events observed in the distribution of lineages. In this work, we explore the hypothesis of a process analogous to ecological-evolutionary radiation, which we called Chromosomal Radiation. Chromosomal Radiation can be adaptive or non-adaptive and applied to different groups of organisms.

A Kuhnian revolution in molecular biology: Most genes in complex organisms express regulatory RNAs

Thomas Kuhn described the progress of science as comprising occasional paradigm shifts separated by interludes of 'normal science'. The paradigm that has held sway since the inception of molecular biology is that genes (mainly) encode proteins. In parallel, theoreticians posited that mutation is random, inferred that most of the genome in complex organisms is non-functional, and asserted that somatic information is not communicated to the germline. However, many anomalies appeared, particularly in plants and animals: the strange genetic phenomena of paramutation and transvection; introns; repetitive sequences; a complex epigenome; lack of scaling of (protein-coding) genes and increase in 'noncoding' sequences with developmental complexity; genetic loci termed 'enhancers' that control spatiotemporal gene expression patterns during development; and a plethora of 'intergenic', overlapping, antisense and intronic transcripts. These observations suggest that the original conception of genetic information was deficient and that most genes in complex organisms specify regulatory RNAs, some of which convey intergenerational information. Also see the video abstract here: https://youtu.be/qxeGwahBANw.

Young SINEs in pig genomes impact gene regulation, genetic diversity, and complex traits

Transposable elements (TEs) are a major source of genetic polymorphisms and play a role in chromatin architecture, gene regulatory networks, and genomic evolution. However, their functional role in pigs and contributions to complex traits are largely unknown. We created a catalog of TEs (n = 3,087,929) in pigs and found that young SINEs were predominantly silenced by histone modifications, DNA methylation, and decreased accessibility. However, some transcripts from active young SINEs showed high tissue-specificity, as confirmed by analyzing 3570 RNA-seq samples. We also detected 211,067 dimorphic SINEs in 374 individuals, including 340 population-specific ones associated with local adaptation. Mapping these dimorphic SINEs to genome-wide associations of 97 complex traits in pigs, we found 54 candidate genes (e.g., ANK2 and VRTN) that might be mediated by TEs. Our findings highlight the important roles of young SINEs and provide a supplement for genotype-to-phenotype associations and modern breeding in pigs.

Holocentric repeat landscapes: From micro-evolutionary patterns to macro-evolutionary associations with karyotype evolution

Repetitive elements can cause large-scale chromosomal rearrangements, for example through ectopic recombination, potentially promoting reproductive isolation and speciation. Species with holocentric chromosomes, that lack a localized centromere, might be more likely to retain chromosomal rearrangements that lead to karyotype changes such as fusions and fissions. This is because chromosome segregation during cell division should be less affected than in organisms with a localized centromere. The relationships between repetitive elements and chromosomal rearrangements and how they may translate to patterns of speciation in holocentric organisms are though poorly understood. Here, we use a reference-free approach based on low-coverage short-read sequencing data to characterize the repeat landscape of two independently evolved holocentric groups: Erebia butterflies and Carex sedges. We consider both micro- and macro-evolutionary scales to investigate the repeat landscape differentiation between Erebia populations and the association between repeats and karyotype changes in a phylogenetic framework for both Erebia and Carex. At a micro-evolutionary scale, we found population differentiation in repeat landscape that increases with overall intraspecific genetic differentiation among four Erebia species. At a macro-evolutionary scale, we found indications for an association between repetitive elements and karyotype changes along both Erebia and Carex phylogenies. Altogether, our results suggest that repetitive elements are associated with the level of population differentiation and chromosomal rearrangements in holocentric clades and therefore likely play a role in adaptation and potentially species diversification.

Epigenetic regulation of nuclear processes in fungal plant pathogens

Through the association of protein complexes to DNA, the eukaryotic nuclear genome is broadly organized into open euchromatin that is accessible for enzymes acting on DNA and condensed heterochromatin that is inaccessible. Chemical and physical alterations to chromatin may impact its organization and functionality and are therefore important regulators of nuclear processes. Studies in various fungal plant pathogens have uncovered an association between chromatin organization and expression of in planta-induced genes that are important for pathogenicity. This review discusses chromatin-based regulation mechanisms as determined in the fungal plant pathogen Verticillium dahliae and relates the importance of epigenetic transcriptional regulation and other nuclear processes more broadly in fungal plant pathogens.

Selective binding of retrotransposons by ZFP352 facilitates the timely dissolution of totipotency network

Acquisition of new stem cell fates relies on the dissolution of the prior regulatory network sustaining the existing cell fates. Currently, extensive insights have been revealed for the totipotency regulatory network around the zygotic genome activation (ZGA) period. However, how the dissolution of the totipotency network is triggered to ensure the timely embryonic development following ZGA is largely unknown. In this study, we identify the unexpected role of a highly expressed 2-cell (2C) embryo specific transcription factor, ZFP352, in facilitating the dissolution of the totipotency network. We find that ZFP352 has selective binding towards two different retrotransposon sub-families. ZFP352 coordinates with DUX to bind the 2C specific MT2_Mm sub-family. On the other hand, without DUX, ZFP352 switches affinity to bind extensively onto SINE_B1/Alu sub-family. This leads to the activation of later developmental programs like ubiquitination pathways, to facilitate the dissolution of the 2C state. Correspondingly, depleting ZFP352 in mouse embryos delays the 2C to morula transition process. Thus, through a shift of binding from MT2_Mm to SINE_B1/Alu, ZFP352 can trigger spontaneous dissolution of the totipotency network. Our study highlights the importance of different retrotransposons sub-families in facilitating the timely and programmed cell fates transition during early embryogenesis.

The Landscape of the DNA Transposons in the Genome of the Horezu_LaPeri Strain of Drosophila melanogaster

Natural transposons (NTs) represent mobile DNA sequences found in both prokaryotic and eukaryotic genomes. Drosophila melanogaster (the fruit fly) is a eukaryotic model organism with NTs standing for about 20% of its genome and has contributed significantly to the understanding of various aspects of transposon biology. Our study describes an accurate approach designed to map class II transposons (DNA transposons) in the genome of the Horezu_LaPeri fruit fly strain, consecutive to Oxford Nanopore Technology sequencing. A whole genome bioinformatics analysis was conducted using Genome ARTIST_v2, LoRTE and RepeatMasker tools to identify DNA transposons insertions. Then, a gene ontology enrichment analysis was performed in order to evaluate the potential adaptive role of some DNA transposons insertions. Herein, we describe DNA transposon insertions specific for the Horezu_LaPeri genome and a predictive functional analysis of some insertional alleles. The PCR validation of P-element insertions specific for this fruit fly strain, along with a putative consensus sequence for the KP element, is also reported. Overall, the genome of the Horezu_LaPeri strain contains several insertions of DNA transposons associated with genes known to be involved in adaptive processes. For some of these genes, insertional alleles obtained via mobilization of the artificial transposons were previously reported. This is a very alluring aspect, as it suggests that insertional mutagenesis experiments conducting adaptive predictions for laboratory strains may be confirmed by mirroring insertions which are expected to be found at least in some natural fruit fly strains.

Pol V produced RNA facilitates transposable element excision site repair in Arabidopsis

The plant-specific RNA Polymerase V (Pol V) plays a key role in gene silencing, but its role in repair of double stranded DNA breaks is unclear. Excision of the transposable element mPing creates double stranded breaks that are repaired by NHEJ. We measured mPing excision site repair in multiple DNA methylation mutants including pol V using an mPing : GFP reporter. Two independent mutant alleles of pol V showed less GFP expression, indicating that the Pol V protein plays a role in excision site repair. Sequence analysis of the pol V excision sites indicated an elevated rate of large deletions consistent with less efficient repair. These results clarify the role of Pol V, but not other RNA-directed DNA methylation proteins (Pol IV) or maintenance DNA methylation pathways ( MET1 ), in the repair of double-strand DNA breaks.

Retrotransposon-mediated evolutionary rewiring of a pathogen response orchestrates a resistance phenotype in an insect host

Ongoing host-pathogen interactions can trigger a coevolutionary arms race, while genetic diversity within the host can facilitate its adaptation to pathogens. Here, we used the diamondback moth (Plutella xylostella) and its pathogen Bacillus thuringiensis (Bt) as a model for exploring an adaptive evolutionary mechanism. We found that insect host adaptation to the primary Bt virulence factors was tightly associated with a short interspersed nuclear element (SINE - named SE2) insertion into the promoter of the transcriptionally activated MAP4K4 gene. This retrotransposon insertion coopts and potentiates the effect of the transcription factor forkhead box O (FOXO) in inducing a hormone-modulated Mitogen-activated protein kinase (MAPK) signaling cascade, leading to an enhancement of a host defense mechanism against the pathogen. This work demonstrates that reconstructing a cis-trans interaction can escalate a host response mechanism into a more stringent resistance phenotype to resist pathogen infection, providing a new insight into the coevolutionary mechanism of host organisms and their microbial pathogens.

The enormous repetitive Antarctic krill genome reveals environmental adaptations and population insights

Antarctic krill (Euphausia superba) is Earth's most abundant wild animal, and its enormous biomass is vital to the Southern Ocean ecosystem. Here, we report a 48.01-Gb chromosome-level Antarctic krill genome, whose large genome size appears to have resulted from inter-genic transposable element expansions. Our assembly reveals the molecular architecture of the Antarctic krill circadian clock and uncovers expanded gene families associated with molting and energy metabolism, providing insights into adaptations to the cold and highly seasonal Antarctic environment. Population-level genome re-sequencing from four geographical sites around the Antarctic continent reveals no clear population structure but highlights natural selection associated with environmental variables. An apparent drastic reduction in krill population size 10 mya and a subsequent rebound 100 thousand years ago coincides with climate change events. Our findings uncover the genomic basis of Antarctic krill adaptations to the Southern Ocean and provide valuable resources for future Antarctic research.

Evolution of woody plants to the land-sea interface - The atypical genomic features of mangroves with atypical phenotypic adaptation

How plants adapt and diverge in extreme environments is a key question of plant evolution and ecology. Mangrove invasion of intertidal environments is facilitated by adaptive phenotypes such as aerial roots, salt-secreting leaf, and viviparity, and genomic mechanisms including whole genome duplication and transposable element number reduction. However, a number of mangroves lack these typical phenotypes. The question we ask is whether these phenotypically atypical mangroves also have distinct genomic features? The sibling mangrove species Lumnitzera littorea and Lumnitzera racemosa provide a model to study this question. We sequenced and assembled their genomes to chromosome level, together with a closely related species Combretum micranthum. While most mangroves have small genomes, the genomes of both Lumnitzera species are large (1443 and 1317 Mb) and carry a high proportion of repeat sequences (~75%). Moreover, Lumnitzera species have not undergone post-gamma whole-genome duplications. Their genome size increased mainly due to the expansion of repeat sequences in their ancestors. However, Lumnitzera genomes have reduced transposable elements by constraining the proliferation of new LTR-RTs. Meanwhile, the two species have more gene families contracted than expanded, and some gene families with reversed size change may underlie their differentiation in root morphology and local distribution. We identified 86 chromosomal inversions, five of which are measured between 6.5 and 12.8 megabases. A number of genes located in these inversions function in pigment biosynthesis, a process likely involved in flower colour differentiation between the Lumnitzera species. We conclude that the mangroves with atypical phenotypes also have atypical genomic evolution.

A chromosome-level genome assembly of Ostrea denselamellosa provides initial insights into its evolution

The oyster Ostrea denselamellosa is a live-bearing species with a sharp decline in the natural population. Despite recent breakthroughs in long-read sequencing, high quality genomic data are very limited in O. denselamellosa. Here, we carried out the first whole genome sequencing at the chromosome-level in O. denselamellosa. Our studies yielded a 636 Mb assembly with scaffold N50 around 71.80 Mb. 608.3 Mb (95.6% of the assembly) were anchored to 10 chromosomes. A total of 26,412 protein-coding genes were predicted, of which 22,636 (85.7%) were functionally annotated. By comparative genomics, we found that long interspersed nuclear element (LINE) and short interspersed nuclear element (SINE) made up a larger proportion in O. denselamellosa genome than in other oysters'. Moreover, gene family analysis showed some initial insight into its evolution. This high-quality genome of O. denselamellosa provides a valuable genomic resource for studies of evolution, adaption and conservation in oysters.

Disentangling Crocus Series Verni and Its Polyploids

Spring crocuses, the eleven species within Crocus series Verni (Iridaceae), consist of di- and tetraploid cytotypes. Among them is a group of polyploids from southeastern Europe with yet-unclear taxonomic affiliation. Crocuses are generally characterized by complex dysploid chromosome number changes, preventing a clear correlation between these numbers and ploidy levels. To reconstruct the evolutionary history of series Verni and particularly its polyploid lineages associated with C. heuffelianus, we used an approach combining phylogenetic analyses of two chloroplast regions, 14 nuclear single-copy genes plus rDNA spacers, genome-wide genotyping-by-sequencing (GBS) data, and morphometry with ploidy estimations through genome size measurements, analysis of genomic heterozygosity frequencies and co-ancestry, and chromosome number counts. Chromosome numbers varied widely in diploids with 2n = 8, 10, 12, 14, 16, and 28 and tetraploid species or cytotypes with 2n = 16, 18, 20, and 22 chromosomes. Crocus longiflorus, the diploid with the highest chromosome number, possesses the smallest genome (2C = 3.21 pg), while the largest diploid genomes are in a range of 2C = 7-8 pg. Tetraploid genomes have 2C values between 10.88 pg and 12.84 pg. Heterozygosity distribution correlates strongly with genome size classes and allows discernment of di- and tetraploid cytotypes. Our phylogenetic analyses showed that polyploids in the C. heuffelianus group are allotetraploids derived from multiple and partly reciprocal crosses involving different genotypes of diploid C. heuffelianus (2n = 10) and C. vernus (2n = 8). Dysploid karyotype changes after polyploidization resulted in the tetraploid cytotypes with 20 and 22 chromosomes. The multi-data approach we used here for series Verni, combining evidence from nuclear and chloroplast phylogenies, genome sizes, chromosome numbers, and genomic heterozygosity for ploidy estimations, provides a way to disentangle the evolution of plant taxa with complex karyotype changes that can be used for the analysis of other groups within Crocus and beyond. Comparing these results with morphometric analysis results in characters that can discern the different taxa currently subsumed under C. heuffelianus.

Chromosomal-level genome assembly of the high-quality Xian/Indica rice (Oryza sativa L.) Xiangyaxiangzhan

The indica rice variety XYXZ carries elite traits including appearance and eating quality. Here, we report the de novo assembly of XYXZ using Illumine paired-end whole-genome shotgun sequencing and Nanopore sequencing. We annotated 39,722 protein-coding genes in the 395.04 Mb assembly. In comparison to other cultivars, XYXZ showed a larger gene size including the transcripts and introns, and more exons per gene. And hundreds of ultra-long genes were also detected. A total of 4362 complete LTRs were annotated, and among them, many were located next to or in protein-coding genes including several genes related to rice quality. We observed the different distributions of LTRs in these genes among XYXZ, Nipponbare, and R498, implying these LTRs might potentially affect expressions of the proximal genes and rice quality. Overall, This chromosome-length genome assembly of XYXZ provides a valuable resource for gene discovery, genetic variation and evolution, and the breeding of high-quality rice.

Down, then up: non-parallel genome size changes and a descending chromosome series in a recent radiation of the Australian allotetraploid plant species, Nicotiana section Suaveolentes (Solanaceae)

BACKGROUND AND AIMS

The extent to which genome size and chromosome numbers evolve in concert is little understood, particularly after polyploidy (whole-genome duplication), when a genome returns to a diploid-like condition (diploidization). We study this phenomenon in 46 species of allotetraploid Nicotiana section Suaveolentes (Solanaceae), which formed <6 million years ago and radiated in the arid centre of Australia.

METHODS

We analysed newly assessed genome sizes and chromosome numbers within the context of a restriction site-associated nuclear DNA (RADseq) phylogenetic framework.

KEY RESULTS

RADseq generated a well-supported phylogenetic tree, in which multiple accessions from each species formed unique genetic clusters. Chromosome numbers and genome sizes vary from n = 2x = 15 to 24 and 2.7 to 5.8 pg/1C nucleus, respectively. Decreases in both genome size and chromosome number occur, although neither consistently nor in parallel. Species with the lowest chromosome numbers (n = 15-18) do not possess the smallest genome sizes and, although N. heterantha has retained the ancestral chromosome complement, n = 2x = 24, it nonetheless has the smallest genome size, even smaller than that of the modern representatives of ancestral diploids.

CONCLUSIONS

The results indicate that decreases in genome size and chromosome number occur in parallel down to a chromosome number threshold, n = 20, below which genome size increases, a phenomenon potentially explained by decreasing rates of recombination over fewer chromosomes. We hypothesize that, more generally in plants, major decreases in genome size post-polyploidization take place while chromosome numbers are still high because in these stages elimination of retrotransposons and other repetitive elements is more efficient. Once such major genome size change has been accomplished, then dysploid chromosome reductions take place to reorganize these smaller genomes, producing species with small genomes and low chromosome numbers such as those observed in many annual angiosperms, including Arabidopsis.

piRNA pathway evolution beyond gonad context: Perspectives from apicomplexa and trypanosomatids

piRNAs function as genome defense mechanisms against transposable elements insertions within germ line cells. Recent studies have unraveled that piRNA pathways are not limited to germ cells as initially reckoned, but are instead also found in non-gonadal somatic contexts. Moreover, these pathways have also been reported in bacteria, mollusks and arthropods, associated with safeguard of genomes against transposable elements, regulation of gene expression and with direct consequences in axon regeneration and memory formation. In this Perspective we draw attention to early branching parasitic protozoa, whose genome preservation is an essential function as in late eukaryotes. However, little is known about the defense mechanisms of these genomes. We and others have described the presence of putative PIWI-related machinery members in protozoan parasites. We have described the presence of a PIWI-like protein in Trypanosoma cruzi, bound to small non-coding RNAs (sRNAs) as cargo of secreted extracellular vesicles relevant in intercellular communication and host infection. Herein, we put forward the presence of members related to Argonaute pathways in both Trypanosoma cruzi and Toxoplasma gondii. The presence of PIWI-like machinery in Trypansomatids and Apicomplexa, respectively, could be evidence of an ancestral piRNA machinery that evolved to become more sophisticated and complex in multicellular eukaryotes. We propose a model in which ancient PIWI proteins were expressed broadly and had functions independent of germline maintenance. A better understanding of current and ancestral PIWI/piRNAs will be relevant to better understand key mechanisms of genome integrity conservation during cell cycle progression and modulation of host defense mechanisms by protozoan parasites.

Genome-wide identification and characterization of exapted transposable elements in the large genome of sunflower (Helianthus annuus L.)

Transposable elements (TEs) are an important source of genome variability, playing many roles in the evolution of eukaryotic species. Besides well-known phenomena, TEs may undergo the exaptation process and generate the so-called exapted transposable element genes (ETEs). Here we present a genome-wide survey of ETEs in the large genome of sunflower (Helianthus annuus L.), in which the massive amount of TEs, provides a significant source for exaptation. A library of sunflower TEs was used to build TE-specific Hidden Markov Model profiles, to search for all available sunflower gene products. In doing so, 20 016 putative ETEs were identified and further investigated for the characteristics that distinguish TEs from genes, leading to the validation of 3530 ETEs. The analysis of ETEs transcription patterns under different stress conditions showed a differential regulation triggered by treatments mimicking biotic and abiotic stress; furthermore, the distribution of functional domains of differentially regulated ETEs revealed a relevant presence of domains involved in many aspects of cellular functions. A comparative genomic investigation was performed including species representative of Asterids and appropriate outgroups: the bulk of ETEs that resulted were specific to the sunflower, while few ETEs presented orthologues in the genome of all analyzed species, making the hypothesis of a conserved function. This study highlights the crucial role played by exaptation, actively contributing to species evolution.

THERMOSENSITIVE BARREN PANICLE (TAP) is required for rice panicle and spikelet development at high ambient temperature

In cereal plants, the size of the panicle (inflorescence) is a critical factor for yield. Panicle size is determined by a complex interplay of genetic and environmental factors, but the mechanisms underlying adaptations to temperature stress during panicle development remain largely unknown. We identify the rice THERMOSENSITIVE BARREN PANICLE (TAP) gene, which encodes a transposase-derived FAR1-RELATED SEQUENCE (FRS) protein and is responsible for regulating panicle and spikelet development at high ambient temperature. The tap mutants display high temperature-dependent reproductive abnormalities, including compromised secondary branch and spikelet initiation and pleiotropic floral organ defects. Consistent with its thermosensitive phenotype, TAP expression is induced by high temperature. TAP directly promotes the expression of OsYABBY3 (OsYAB3), OsYAB4, and OsYAB5, which encode key transcriptional regulators in panicle and spikelet development. In addition, TAP physically interacts with OsYAB4 and OsYAB5 proteins; phenotypic analysis of osyab4 tap-1 and osyab5 tap-1 double mutants indicates that TAP-OsYAB4/OsYAB5 complexes act to maintain normal panicle and spikelet development. Taken together, our study reveals the novel role of a TE-derived transcription factor in controlling rice panicle development under high ambient temperatures, shedding light on the molecular mechanism underlying the adaptation of cereal crops to increasing environmental temperatures.

B2 SINE RNA as a novel regulator of glucocorticoid receptor transcriptional activity

Glucocorticoids are a key component to the cellular response to stress. Glucocorticoids act via glucocorticoid receptors found ubiquitously in the brain and body. Glucocorticoid receptors can bind to response elements throughout the genome to elicit changes in transcription, an adaptation observed at the cellular level. Yet, the transcriptional changes as a consequence of glucocorticoid receptor activation are variable across brain regions, stress conditions and recurrent bouts of glucocorticoid exposure. Here we describe a non-coding RNA, B2 SINE, which is regulated by glucocorticoids and can in turn regulate glucocorticoid receptor transcriptional activity. We show that activated glucocorticoid receptors interact directly with B2 SINE RNA via a decoy response element contained within the transcript sequence and alter receptor binding to response elements in the genome and, subsequently, changes in loci expression.

WiFi Related Radiofrequency Electromagnetic Fields Promote Transposable Element Dysregulation and Genomic Instability in Drosophila melanogaster

Exposure to artificial radio frequency electromagnetic fields (RF-EMFs) has greatly increased in recent years, thus promoting a growing scientific and social interest in deepening the biological impact of EMFs on living organisms. The current legislation governing the exposure to RF-EMFs is based exclusively on their thermal effects, without considering the possible non-thermal adverse health effects from long term exposure to EMFs. In this study we investigated the biological non-thermal effects of low-level indoor exposure to RF-EMFs produced by WiFi wireless technologies, using Drosophila melanogaster as the model system. Flies were exposed to 2.4 GHz radiofrequency in a Transverse Electromagnetic (TEM) cell device to ensure homogenous controlled fields. Signals were continuously monitored during the experiments and regulated at non thermal levels. The results of this study demonstrate that WiFi electromagnetic radiation causes extensive heterochromatin decondensation and thus a general loss of transposable elements epigenetic silencing in both germinal and neural tissues. Moreover, our findings provide evidence that WiFi related radiofrequency electromagnetic fields can induce reactive oxygen species (ROS) accumulation, genomic instability, and behavioural abnormalities. Finally, we demonstrate that WiFi radiation can synergize with Ras^V12 to drive tumor progression and invasion. All together, these data indicate that radiofrequency radiation emitted from WiFi devices could exert genotoxic effects in Drosophila and set the stage to further explore the biological effects of WiFi electromagnetic radiation on living organisms.

The nematode Oscheius tipulae as a genetic model for programmed DNA elimination

Programmed DNA elimination (PDE) is a notable exception to the paradigm of genome integrity. In metazoa, PDE often occurs coincident with germline to somatic cell differentiation. During PDE, portions of genomic DNA are lost, resulting in reduced somatic genomes. Prior studies have described the sequences lost, as well as chromosome behavior, during metazoan PDE. However, a system for studying the mechanisms and consequences of PDE in metazoa is lacking. Here, we present a functional and genetic model for PDE in the free-living Rhabditidae nematode Oscheius tipulae, a family that also includes Caenorhabditis elegans. O. tipulae was recently suggested to eliminate DNA. Using staged embryos and DNA FISH, we showed that O. tipulae PDE occurs during embryogenesis at the 8-16 cell stages. We identified a conserved motif, named Sequence For Elimination (SFE), for all 12 break sites on the six chromosomes at the junctions of retained and eliminated DNA. SFE mutants exhibited a "fail-to-eliminate" phenotype only at the modified sites. END-seq revealed that breaks can occur at multiple positions within the SFE, with extensive end resection followed by telomere addition to both retained and eliminated ends. We identified many functional SFEs at the chromosome ends through END-seq in the wild-type embryos, genome sequencing of SFE mutants, and comparative genomics of 23 wild isolates. We suggest that these alternative SFEs provide flexibility in the sequences eliminated and a fail-safe mechanism for PDE. These studies establish O. tipulae as a new, attractive model for studying the mechanisms and consequences of PDE in a metazoan.

Legume Pangenome: Status and Scope for Crop Improvement

In the last decade, legume genomics research has seen a paradigm shift due to advances in genome sequencing technologies, assembly algorithms, and computational genomics that enabled the construction of high-quality reference genome assemblies of major legume crops. These advances have certainly facilitated the identification of novel genetic variants underlying the traits of agronomic importance in many legume crops. Furthermore, these robust sequencing technologies have allowed us to study structural variations across the whole genome in multiple individuals and at the species level using 'pangenome analysis.' This review updates the progress of constructing pangenome assemblies for various legume crops and discusses the prospects for these pangenomes and how to harness the information to improve various traits of economic importance through molecular breeding to increase genetic gain in legumes and tackle the increasing global food crisis.

Alu Deletions in LAMA2 and CDH4 Genes Are Key Components of Polygenic Predictors of Longevity

Longevity is a unique human phenomenon and a highly stable trait, characterized by polygenicity. The longevity phenotype occurs due to the ability to successfully withstand the age-related genomic instability triggered by Alu elements. The purpose of our cross-sectional study was to evaluate the combined contribution of ACE*Ya5ACE, CDH4*Yb8NBC516, COL13A1*Ya5ac1986, HECW1*Ya5NBC182, LAMA2*Ya5-MLS19, PLAT*TPA25, PKHD1L1*Yb8AC702, SEMA6A*Yb8NBC597, STK38L*Ya5ac2145 and TEAD1*Ya5ac2013 Alu elements to longevity. The study group included 2054 unrelated individuals aged from 18 to 113 years who are ethnic Tatars from Russia. We analyzed the dynamics of the allele and genotype frequencies of the studied Alu polymorphic loci in the age groups of young (18-44 years old), middle-aged (45-59 years old), elderly (60-74 years old), old seniors (75-89 years old) and long-livers (90-113 years old). Most significant changes in allele and genotype frequencies were observed between the long-livers and other groups. The search for polygenic predictors of longevity was performed using the APSampler program. Attaining longevity was associated with the combinations LAMA2*ID + CDH4*D (OR = 2.23, PBonf = 1.90 × 10-2) and CDH4*DD + LAMA2*ID + HECW1*D (OR = 4.58, PBonf = 9.00 × 10-3) among persons aged between 18 and 89 years, LAMA2*ID + CDH4*D + SEMA6A*I for individuals below 75 years of age (OR = 3.13, PBonf = 2.00 × 10-2), LAMA2*ID + HECW1*I for elderly people aged 60 and older (OR = 3.13, PBonf = 2.00 × 10-2) and CDH4*DD + LAMA2*D + HECW1*D (OR = 4.21, PBonf = 2.60 × 10-2) and CDH4*DD + LAMA2*D + ACE*I (OR = 3.68, PBonf = 1.90 × 10-2) among old seniors (75-89 years old). The key elements of combinations associated with longevity were the deletion alleles of CDH4 and LAMA2 genes. Our results point to the significance for human longevity of the Alu polymorphic loci in CDH4, LAMA2, HECW1, SEMA6A and ACE genes, involved in the integration systems.

Genomic signatures associated with maintenance of genome stability and venom turnover in two parasitoid wasps

Parasitoid wasps are rapidly developing as a model for evolutionary biology. Here we present chromosomal genomes of two Anastatus wasps, A. japonicus and A. fulloi, and leverage these genomes to study two fundamental questions-genome size evolution and venom evolution. Anastatus shows a much larger genome than is known among other wasps, with unexpectedly recent bursts of LTR retrotransposons. Importantly, several genomic innovations, including Piwi gene family expansion, ubiquitous Piwi expression profiles, as well as transposable element-piRNA coevolution, have likely emerged for transposable element silencing to maintain genomic stability. Additionally, we show that the co-option evolution arose by expression shifts in the venom gland plays a dominant role in venom turnover. We also highlight the potential importance of non-venom genes that are coexpressed with venom genes during venom evolution. Our findings greatly advance the current understanding of genome size evolution and venom evolution, and these genomic resources will facilitate comparative genomics studies of insects in the future.

Genome analysis of Phrixothrix hirtus (Phengodidae) railroad worm shows the expansion of odorant-binding gene families and positive selection on morphogenesis and sex determination genes

Among bioluminescent beetles of the Elateroidea superfamily, Phengodidae is the third largest family, with 244 bioluminescent species distributed only in the Americas, but is still the least studied from the phylogenetic and evolutionary points of view. The railroad worm Phrixothrix hirtus is an essential biological model and symbolic species due to its bicolor bioluminescence, being the only organism that produces true red light among bioluminescent terrestrial species. Here, we performed partial genome assembly of P. hirtus, combining short and long reads generated with Illumina sequencing, providing the first source of genomic information and a framework for comparative analyses of the bioluminescent system in Elateroidea. This is the largest genome described in the Elateroidea superfamily, with an estimated size of ∼3.4 Gb, displaying 32 % GC content, and 67 % transposable elements. Comparative genomic analyses showed a positive selection of genes and gene family expansion events of growths and morphogenesis gene products, which could be associated with the atypical anatomical development and morphogenesis found in paedomorphic females and underdeveloped males. We also observed gene family expansion among distinct odorant-binding receptors, which could be associated with the pheromone communication system typical of these beetles, and retrotransposable elements. Common genes putatively regulating bioluminescence production and control, including two luciferase genes corresponding to lateral lanterns green-emitting and head lanterns red-emitting luciferases with 7 exons and 6 introns, and genes potentially involved in luciferin biosynthesis were found, indicating that there are no clear differences about the presence or absence of gene families associated with bioluminescence in Elateroidea.

Krüppel-like factor 5 rewires NANOG regulatory network to activate human naive pluripotency specific LTR7Ys and promote naive pluripotency

Endogenous retroviruses (ERVs) have been reported to participate in pre-implantation development of mammalian embryos. In early human embryogenesis, different ERV sub-families are activated in a highly stage-specific manner. How the specificity of ERV activation is achieved remains largely unknown. Here, we demonstrate the mechanism of how LTR7Ys, the human morula-blastocyst-specific HERVH long terminal repeats, are activated by the naive pluripotency transcription network. We find that KLF5 interacts with and rewires NANOG to bind and regulate LTR7Ys; in contrast, the primed-specific LTR7s are preferentially bound by NANOG in the absence of KLF5. The specific activation of LTR7Ys by KLF5 and NANOG in pluripotent stem cells contributes to human-specific naive pluripotency regulation. KLF5-LTR7Y axis also promotes the expression of trophectoderm genes and contributes to the expanded cell potential toward extra-embryonic lineage. Our study suggests that HERVs are activated by cell-state-specific transcription machinery and promote stage-specific transcription network and cell potency.

Evolution of complex genome architecture in gymnosperms

Gymnosperms represent an ancient lineage that diverged from early spermatophytes during the Devonian. The long fossil records and low diversity in living species prove their complex evolutionary history, which included ancient radiations and massive extinctions. Due to their ultra-large genome size, the whole-genome assembly of gymnosperms has only generated in the past 10 years and is now being further expanded into more taxonomic representations. Here, we provide an overview of the publicly available gymnosperm genome resources and discuss their assembly quality and recent findings in large genome architectures. In particular, we describe the genomic features most related to changes affecting the whole genome. We also highlight new realizations relative to repetitive sequence dynamics, paleopolyploidy, and long introns. Based on the results of relevant genomic studies of gymnosperms, we suggest additional efforts should be made toward exploring the genomes of medium-sized (5–15 gigabases) species. Lastly, more comparative analyses among high-quality assemblies are needed to understand the genomic shifts and the early species diversification of seed plants.

Environmental stress and transposons in plants

Transposons were once thought to be junk repetitive DNA in the genome. However, their importance gradually became apparent as it became clear that they regulate gene expression, which is essential for organisms to survive, and that they are important factors in the driving force of evolution. Since there are multiple transposons in the genomes of all organisms, transposons have likely been activated and increased in copy number throughout their long history. This review focuses on environmental stress as a factor in transposon activation, paying particular attention to transposons in plants that are activated by environmental stresses. It is now known that plants respond to environmental stress in various ways, and correspondingly, many transposons respond to stress. The relationship between environmental stress and transposons is reviewed, including the mechanisms of their activation and the effects of transposon activation on host plants.

Complexification of eukaryote phenotype: Adaptive immuno-cognitive systems as unique Gödelian block chain distributed ledger

The digitization of inheritable information in the genome has been called the 'algorithmic take-over of biology'. The McClintock discovery that viral software based transposable elements that conduct cut-paste (transposon) and copy-paste (retrotransposon) operations are needed for genomic evolvability underscores the truism that only software can change software and also that viral hacking by internal and external bio-malware is the Achilles heel of genomic digital systems. There was a paradigm shift in genomic information processing with the Adaptive Immune System (AIS) 500 mya followed by the Mirror Neuron System (MNS), latterly mostly in primate brains, which reaches its apogee in human social cognition. The AIS and MNS involve distinctive Gödelian features of self-reference (Self-Ref) and offline virtual self-representation (Self-Rep) for complex self-other interaction with prodigious open-ended capacity for anticipative malware detection and novelty production within a unique blockchain distributed ledger (BCDL). The role of self-referential information processing, often considered to be central to the sentient self with origins in the immune system 'Thymic self', is shown to be part of the Gödel logic behind a generator-selector framework at a molecular level, which exerts stringent selection criteria to maintain genomic BCDL. The latter manifests digital and decentralized record keeping where no internal or external bio-malware can compromise the immutability of the life's building blocks and no novel blocks can be added that is not consistent with extant blocks. This is demonstrated with regard to somatic hypermutation with novel anti-body production in the face of external non-self antigen attacks.

Evolution and origin of bread wheat

Bread wheat (Triticum aestivum, genome BBAADD) is a young hexaploid species formed only 8,500-9,000 years ago through hybridization between a domesticated free-threshing tetraploid progenitor, genome BBAA, and Aegilops tauschii, the diploid donor of the D subgenome. Very soon after its formation, it spread globally from its cradle in the fertile crescent into new habitats and climates, to become a staple food of humanity. This extraordinary global expansion was probably enabled by allopolyploidy that accelerated genetic novelty through the acquisition of new traits, new intergenomic interactions, and buffering of mutations, and by the attractiveness of bread wheat's large, tasty, and nutritious grain with high baking quality. New genome sequences suggest that the elusive donor of the B subgenome is a distinct (unknown or extinct) species rather than a mosaic genome. We discuss the origin of the diploid and tetraploid progenitors of bread wheat and the conflicting genetic and archaeological evidence on where it was formed and which species was its free-threshing tetraploid progenitor. Wheat experienced many environmental changes throughout its evolution, therefore, while it might adapt to current climatic changes, efforts are needed to better use and conserve the vast gene pool of wheat biodiversity on which our food security depends.