Transposable Element Expression Profiles in Premalignant Pigment Cell Lesions and Melanoma of Xiphophorus.

Transposable elements (TEs) are characterized by their ability to change their genomic position. Through insertion or recombination leading to deletions and other chromosomal aberrations, they can cause genetic instability. The extent to which they thereby exert regulatory influence on cellular functions is unclear. To better characterize TEs in processes such as carcinogenesis, we used the well-established Xiphophorus melanoma model. By transcriptome sequencing, we show that an increasing total number in transposons correlates with progression of malignancy in melanoma samples from Xiphophorus interspecific hybrids. Further, by comparing the presence of TEs in the parental genomes of Xiphophorus maculatus and Xiphophorus hellerii, we could show that even in closely related species, genomic location and spectrum of TEs are considerably different.

THBS1+ myeloid cells expand in SLD hepatocellular carcinoma and contribute to immunosuppression and unfavorable prognosis through TREM1.

Hepatocellular carcinoma (HCC) is an inflammation-associated cancer arising from viral or non-viral etiologies including steatotic liver diseases (SLDs). Expansion of immunosuppressive myeloid cells is a hallmark of inflammation and cancer, but their heterogeneity in HCC is not fully resolved and might underlie immunotherapy resistance. Here, we present a high-resolution atlas of innate immune cells from patients with HCC that unravels an SLD-associated contexture characterized by influx of inflammatory and immunosuppressive myeloid cells, including a discrete population of THBS1+ regulatory myeloid (Mreg) cells expressing monocyte- and neutrophil-affiliated genes. THBS1+ Mreg cells expand in SLD-associated HCC, populate fibrotic lesions, and are associated with poor prognosis. THBS1+ Mreg cells are CD163+ but distinguished from macrophages by high expression of triggering receptor expressed on myeloid cells 1 (TREM1), which contributes to their immunosuppressive activity and promotes HCC tumor growth in vivo. Our data support myeloid subset-targeted immunotherapies to treat HCC.

Integrated Genomic Analyses From Low-Depth Sequencing Help Resolve Phylogenetic Incongruence in the Bamboos (Poaceae: Bambusoideae).

The bamboos (Bambusoideae, Poaceae) comprise a major grass lineage with a complex evolutionary history involving ancient hybridization and allopolyploidy. About 1700 described species are classified into three tribes, Olyreae (herbaceous bamboos), Bambuseae (tropical woody bamboos), and Arundinarieae (temperate woody bamboos). Nuclear analyses strongly support monophyly of the woody tribes, whereas plastome analyses strongly support paraphyly, with Bambuseae sister to Olyreae. Our objectives were to clarify the origin(s) of the woody bamboo tribes and resolve the nuclear vs. plastid conflict using genomic tools. For the first time, plastid and nuclear genomic information from the same bamboo species were combined in a single study. We sampled 51 species of bamboos representing the three tribes, estimated their genome sizes and generated low-depth sample sequence data, from which plastomes were assembled and nuclear repeats were analyzed. The distribution of repeat families was found to agree with nuclear gene phylogenies, but also provides novel insights into nuclear evolutionary history. We infer two early, independent hybridization events, one between an Olyreae ancestor and a woody ancestor giving rise to the two Bambuseae lineages, and another between two woody ancestors giving rise to the Arundinarieae. Retention of the Olyreae plastome associated with differential dominance of nuclear genomes and subsequent diploidization in some lineages explains the paraphyly observed in plastome phylogenetic estimations. We confirm ancient hybridization and allopolyploidy in the origins of the extant woody bamboo lineages and propose biased fractionation and diploidization as important factors in their evolution.

Hepatocellular Carcinoma Immune Landscape and the Potential of Immunotherapies.

Hepatocellular carcinoma (HCC) is the most common liver tumor and among the deadliest cancers worldwide. Advanced HCC overall survival is meager and has not improved over the last decade despite approval of several tyrosine kinase inhibitors (TKi) for first and second-line treatments. The recent approval of immune checkpoint inhibitors (ICI) has revolutionized HCC palliative care. Unfortunately, the majority of HCC patients fail to respond to these therapies. Here, we elaborate on the immune landscapes of the normal and cirrhotic livers and of the unique HCC tumor microenvironment. We describe the molecular and immunological classifications of HCC, discuss the role of specific immune cell subsets in this cancer, with a focus on myeloid cells and pathways in anti-tumor immunity, tumor promotion and immune evasion. We also describe the challenges and opportunities of immunotherapies in HCC and discuss new avenues based on harnessing the anti-tumor activity of myeloid, NK and γδ T cells, vaccines, chimeric antigen receptors (CAR)-T or -NK cells, oncolytic viruses, and combination therapies.

Individual knock out of glycine receptor alpha subunits identifies a specific requirement of glra1 for motor function in zebrafish.

Glycine receptors (GlyRs) are ligand-gated chloride channels mediating inhibitory neurotransmission in the brain stem and spinal cord. They function as pentamers composed of alpha and beta subunits for which 5 genes have been identified in human (GLRA1, GLRA2, GLRA3, GLRA4, GLRB). Several in vitro studies showed that the pentameric subtype composition as well as its stoichiometry influence the distribution and the molecular function of the receptor. Moreover, mutations in some of these genes are involved in different human conditions ranging from tinnitus to epilepsy and hyperekplexia, suggesting distinct functions of the different subunits. Although the beta subunit is essential for synaptic clustering of the receptor, the specific role of each alpha subtype is still puzzling in vivo. The zebrafish genome encodes for five glycine receptor alpha subunits (glra1, glra2, glra3, glra4a, glra4b) thus offering a model of choice to investigate the respective role of each subtype on general motor behaviour. After establishing a phylogeny of GlyR subunit evolution between human and zebrafish, we checked the temporal expression pattern of these transcripts during embryo development. Interestingly, we found that glra1 is the only maternally transmitted alpha subunit. We also showed that the expression of the different GlyR subunits starts at different time points during development. Lastly, in order to decipher the role of each alpha subunit on the general motor behaviour of the fish, we knocked out individually each alpha subunit by CRISPR/Cas9-targeted mutagenesis. Surprisingly, we found that knocking out any of the alpha2, 3, a4a or a4b subunit did not lead to any obvious developmental or motor phenotype. However, glra1-/- (hitch) embryos depicted a strong motor dysfunction from 3 days, making them incapable to swim and thus leading to their premature death. Our results infer a strong functional redundancy between alpha subunits and confirm the central role played by glra1 for proper inhibitory neurotransmission controlling locomotion. The genetic tools we developed here will be of general interest for further studies aiming at dissecting the role of GlyRs in glycinergic transmission in vivo and the hitch mutant (hic) is of specific relevance as a new model of hyperekplexia.

A High-Quality Reference Genome for the Invasive Mosquitofish Gambusia affinis Using a Chicago Library.

The western mosquitofish, Gambusia affinis, is a freshwater poecilid fish native to the southeastern United States but with a global distribution due to widespread human introduction. Gambusia affinis has been used as a model species for a broad range of evolutionary and ecological studies. We sequenced the genome of a male G. affinis to facilitate genetic studies in diverse fields including invasion biology and comparative genetics. We generated Illumina short read data from paired-end libraries and in vitro proximity-ligation libraries. We obtained 54.9× coverage, N50 contig length of 17.6 kb, and N50 scaffold length of 6.65 Mb. Compared to two other species in the Poeciliidae family, G. affinis has slightly fewer genes that have shorter total, exon, and intron length on average. Using a set of universal single-copy orthologs in fish genomes, we found 95.5% of these genes were complete in the G. affinis assembly. The number of transposable elements in the G. affinis assembly is similar to those of closely related species. The high-quality genome sequence and annotations we report will be valuable resources for scientists to map the genetic architecture of traits of interest in this species.

Clonal polymorphism and high heterozygosity in the celibate genome of the Amazon molly.

The extreme rarity of asexual vertebrates in nature is generally explained by genomic decay due to absence of meiotic recombination, thus leading to extinction of such lineages. We explore features of a vertebrate asexual genome, the Amazon molly, Poecilia formosa, and find few signs of genetic degeneration but unique genetic variability and ongoing evolution. We uncovered a substantial clonal polymorphism and, as a conserved feature from its interspecific hybrid origin, a 10-fold higher heterozygosity than in the sexual parental species. These characteristics seem to be a principal reason for the unpredicted fitness of this asexual vertebrate. Our data suggest that asexual vertebrate lineages are scarce not because they are at a disadvantage, but because the genomic combinations required to bypass meiosis and to make up a functioning hybrid genome are rarely met in nature.

Analysis of the spotted gar genome suggests absence of causative link between ancestral genome duplication and transposable element diversification in teleost fish.

Teleost fish have been shown to contain many superfamilies of transposable elements (TEs) that are absent from most tetrapod genomes. Since theories predict an increase in TE activity following polyploidization, such diversity might be linked to the 3R whole-genome duplication that occurred approximately 300 million years ago before the teleost radiation. To test this hypothesis, we have analyzed the genome of the spotted gar Lepisosteus oculatus, which diverged from the teleost lineage before the 3R duplication. Our results indicate that TE diversity and copy numbers are similar in gar and teleost genomes, suggesting that TE diversity was ancestral and not linked to the 3R whole-genome duplication. We propose that about 25 distinct superfamilies of TEs were present in the last ancestor of gars and teleost fish about 300 million years ago in the ray-finned fish lineage.

Germ cell and tumor associated piRNAs in the medaka and Xiphophorus melanoma models.

BACKGROUND: A growing number of studies report an abnormal expression of Piwi-interacting RNAs (piRNAs) and the piRNA processing enzyme Piwi in many cancers. Whether this finding is an epiphenomenon of the chaotic molecular biology of the fast dividing, neoplastically transformed cells or is functionally relevant to tumorigenesisis is difficult to discern at present. To better understand the role of piRNAs in cancer development small laboratory fish models can make a valuable contribution. However, little is known about piRNAs in somatic and neoplastic tissues of fish. RESULTS: To identify piRNA clusters that might be involved in melanoma pathogenesis, we use several transgenic lines of medaka, and platyfish/swordtail hybrids, which develop various types of melanoma. In these tumors Piwi, is expressed at different levels, depending on tumor type. To quantify piRNA levels, whole piRNA populations of testes and melanomas of different histotypes were sequenced. Because no reference piRNA cluster set for medaka or Xiphophorus was yet available we developed a software pipeline to detect piRNA clusters in our samples and clusters were selected that were enriched in one or more samples. We found several loci to be overexpressed or down-regulated in different melanoma subtypes as compared to hyperpigmented skin. Furthermore, cluster analysis revealed a clear distinction between testes, low-grade and high-grade malignant melanoma in medaka. CONCLUSIONS: Our data imply that dysregulation of piRNA expression may be associated with development of melanoma. Our results also reinforce the importance of fish as a suitable model system to study the role of piRNAs in tumorigenesis.

The spotted gar genome illuminates vertebrate evolution and facilitates human-teleost comparisons.

To connect human biology to fish biomedical models, we sequenced the genome of spotted gar (Lepisosteus oculatus), whose lineage diverged from teleosts before teleost genome duplication (TGD). The slowly evolving gar genome has conserved in content and size many entire chromosomes from bony vertebrate ancestors. Gar bridges teleosts to tetrapods by illuminating the evolution of immunity, mineralization and development (mediated, for example, by Hox, ParaHox and microRNA genes). Numerous conserved noncoding elements (CNEs; often cis regulatory) undetectable in direct human-teleost comparisons become apparent using gar: functional studies uncovered conserved roles for such cryptic CNEs, facilitating annotation of sequences identified in human genome-wide association studies. Transcriptomic analyses showed that the sums of expression domains and expression levels for duplicated teleost genes often approximate the patterns and levels of expression for gar genes, consistent with subfunctionalization. The gar genome provides a resource for understanding evolution after genome duplication, the origin of vertebrate genomes and the function of human regulatory sequences.

X. couchianus and X. hellerii genome models provide genomic variation insight among Xiphophorus species.

BACKGROUND: Xiphophorus fishes are represented by 26 live-bearing species of tropical fish that express many attributes (e.g., viviparity, genetic and phenotypic variation, ecological adaptation, varied sexual developmental mechanisms, ability to produce fertile interspecies hybrids) that have made attractive research models for over 85 years. Use of various interspecies hybrids to investigate the genetics underlying spontaneous and induced tumorigenesis has resulted in the development and maintenance of pedigreed Xiphophorus lines specifically bred for research. The recent availability of the X. maculatus reference genome assembly now provides unprecedented opportunities for novel and exciting comparative research studies among Xiphophorus species. RESULTS: We present sequencing, assembly and annotation of two new genomes representing Xiphophorus couchianus and Xiphophorus hellerii. The final X. couchianus and X. hellerii assemblies have total sizes of 708 Mb and 734 Mb and correspond to 98 % and 102 % of the X. maculatus Jp 163 A genome size, respectively. The rates of single nucleotide change range from 1 per 52 bp to 1 per 69 bp among the three genomes and the impact of putatively damaging variants are presented. In addition, a survey of transposable elements allowed us to deduce an ancestral TE landscape, uncovered potential active TEs and document a recent burst of TEs during evolution of this genus. CONCLUSIONS: Two new Xiphophorus genomes and their corresponding transcriptomes were efficiently assembled, the former using a novel guided assembly approach. Three assembled genome sequences within this single vertebrate order of new world live-bearing fishes will accelerate our understanding of relationship between environmental adaptation and genome evolution. In addition, these genome resources provide capability to determine allele specific gene regulation among interspecies hybrids produced by crossing any of the three species that are known to produce progeny predisposed to tumor development.

Insights into Sex Chromosome Evolution and Aging from the Genome of a Short-Lived Fish.

The killifish Nothobranchius furzeri is the shortest-lived vertebrate that can be bred in the laboratory. Its rapid growth, early sexual maturation, fast aging, and arrested embryonic development (diapause) make it an attractive model organism in biomedical research. Here, we report a draft sequence of its genome that allowed us to uncover an intra-species Y chromosome polymorphism representing-in real time-different stages of sex chromosome formation that display features of early mammalian XY evolution "in action." Our data suggest that gdf6Y, encoding a TGF-ß family growth factor, is the master sex-determining gene in N. furzeri. Moreover, we observed genomic clustering of aging-related genes, identified genes under positive selection, and revealed significant similarities of gene expression profiles between diapause and aging, particularly for genes controlling cell cycle and translation. The annotated genome sequence is provided as an online resource (http://www.nothobranchius.info/NFINgb).

The coelacanth: Can a "living fossil" have active transposable elements in its genome?

The coelacanth has long been regarded as a "living fossil," with extant specimens looking very similar to fossils dating back to the Cretaceous period. The hypothesis of a slowly or even not evolving genome has been proposed to account for this apparent morphological stasis. While this assumption seems to be sustained by different evolutionary analyses on protein-coding genes, recent studies on transposable elements have provided more conflicting results. Indeed, the coelacanth genome contains many transposable elements and has been shaped by several major bursts of transposition during evolution. In addition, comparison of orthologous genomic regions from the genomes of the 2 extant coelacanth species L. chalumnae and L. menadoensis revealed multiple species-specific insertions, indicating transposable element recent activity and contribution to post-speciation genome divergence. These observations, which do not support the genome stasis hypothesis, challenge either the impact of transposable elements on organismal evolution or the status of the coelacanth as a "living fossil." Closer inspection of fossil and molecular data indicate that, even if coelacanths might evolve more slowly than some other lineages due to demographic and/or ecological factors, this variation is still in the range of a "non-fossil" vertebrate species.

Transposable elements and early evolution of sex chromosomes in fish.

In many organisms, the sex chromosome pair can be recognized due to heteromorphy; the Y and W chromosomes have often lost many genes due to the absence of recombination during meiosis and are frequently heterochromatic. Repetitive sequences are found at a high proportion on such heterochromatic sex chromosomes and the evolution and emergence of sex chromosomes has been connected to the dynamics of repeats and transposable elements. With an amazing plasticity of sex determination mechanisms and numerous instances of independent emergence of novel sex chromosomes, fish represent an excellent lineage to investigate the early stages of sex chromosome differentiation, where sex chromosomes often are homomorphic and not heterochromatic. We have analyzed the composition, distribution, and relative age of TEs from available sex chromosome sequences of seven teleost fish. We observed recent bursts of TEs and simple repeat accumulations around young sex determination loci. More strikingly, we detected transposable element (TE) amplifications not only on the sex determination regions of the Y and W sex chromosomes, but also on the corresponding regions of the X and Z chromosomes. In one species, we also clearly demonstrated that the observed TE-rich sex determination locus originated from a TE-poor genomic region, strengthening the link between TE accumulation and emergence of the sex determination locus. Altogether, our results highlight the role of TEs in the initial steps of differentiation and evolution of sex chromosomes.

Evolutionary impact of transposable elements on genomic diversity and lineage-specific innovation in vertebrates.

Since their discovery, a growing body of evidence has emerged demonstrating that transposable elements are important drivers of species diversity. These mobile elements exhibit a great variety in structure, size and mechanisms of transposition, making them important putative actors in organism evolution. The vertebrates represent a highly diverse and successful lineage that has adapted to a wide range of different environments. These animals also possess a rich repertoire of transposable elements, with highly diverse content between lineages and even between species. Here, we review how transposable elements are driving genomic diversity and lineage-specific innovation within vertebrates. We discuss the large differences in TE content between different vertebrate groups and then go on to look at how they affect organisms at a variety of levels: from the structure of chromosomes to their involvement in the regulation of gene expression, as well as in the formation and evolution of non-coding RNAs and protein-coding genes. In the process of doing this, we highlight how transposable elements have been involved in the evolution of some of the key innovations observed within the vertebrate lineage, driving the group's diversity and success.

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.

Guidelines for the nomenclature of genetic elements in tunicate genomes.

Tunicates are invertebrate members of the chordate phylum, and are considered to be the sister group of vertebrates. Tunicates are composed of ascidians, thaliaceans, and appendicularians. With the advent of inexpensive high-throughput sequencing, the number of sequenced tunicate genomes is expected to rise sharply within the coming years. To facilitate comparative genomics within the tunicates, and between tunicates and vertebrates, standardized rules for the nomenclature of tunicate genetic elements need to be established. Here we propose a set of nomenclature rules, consensual within the community, for predicted genes, pseudogenes, transcripts, operons, transcriptional cis-regulatory regions, transposable elements, and transgenic constructs. In addition, the document proposes guidelines for naming transgenic and mutant lines.

Interspecies insertion polymorphism analysis reveals recent activity of transposable elements in extant coelacanths.

Coelacanths are lobe-finned fish represented by two extant species, Latimeria chalumnae in South Africa and Comoros and L. menadoensis in Indonesia. Due to their intermediate phylogenetic position between ray-finned fish and tetrapods in the vertebrate lineage, they are of great interest from an evolutionary point of view. In addition, extant specimens look similar to 300 million-year-old fossils; because of their apparent slowly evolving morphology, coelacanths have been often described as « living fossils ¼. As an underlying cause of such a morphological stasis, several authors have proposed a slow evolution of the coelacanth genome. Accordingly, sequencing of the L. chalumnae genome has revealed a globally low substitution rate for protein-coding regions compared to other vertebrates. However, genome and gene evolution can also be influenced by transposable elements, which form a major and dynamic part of vertebrate genomes through their ability to move, duplicate and recombine. In this work, we have searched for evidence of transposition activity in coelacanth genomes through the comparative analysis of orthologous genomic regions from both Latimeria species. Comparison of 5.7 Mb (0.2%) of the L. chalumnae genome with orthologous Bacterial Artificial Chromosome clones from L. menadoensis allowed the identification of 27 species-specific transposable element insertions, with a strong relative contribution of CR1 non-LTR retrotransposons. Species-specific homologous recombination between the long terminal repeats of a new coelacanth endogenous retrovirus was also detected. Our analysis suggests that transposon activity is responsible for at least 0.6% of genome divergence between both Latimeria species. Taken together, this study demonstrates that coelacanth genomes are not evolutionary inert: they contain recently active transposable elements, which have significantly contributed to post-speciation genome divergence in Latimeria.

The cavefish genome reveals candidate genes for eye loss.

Natural populations subjected to strong environmental selection pressures offer a window into the genetic underpinnings of evolutionary change. Cavefish populations, Astyanax mexicanus (Teleostei: Characiphysi), exhibit repeated, independent evolution for a variety of traits including eye degeneration, pigment loss, increased size and number of taste buds and mechanosensory organs, and shifts in many behavioural traits. Surface and cave forms are interfertile making this system amenable to genetic interrogation; however, lack of a reference genome has hampered efforts to identify genes responsible for changes in cave forms of A. mexicanus. Here we present the first de novo genome assembly for Astyanax mexicanus cavefish, contrast repeat elements to other teleost genomes, identify candidate genes underlying quantitative trait loci (QTL), and assay these candidate genes for potential functional and expression differences. We expect the cavefish genome to advance understanding of the evolutionary process, as well as, analogous human disease including retinal dysfunction.