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 genomes of all lungfish inform on genome expansion and tetrapod evolution.

The genomes of living lungfishes can inform on the molecular-developmental basis of the Devonian sarcopterygian fish-tetrapod transition. We de novo sequenced the genomes of the African (Protopterus annectens) and South American lungfishes (Lepidosiren paradoxa). The Lepidosiren genome (about 91 Gb, roughly 30 times the human genome) is the largest animal genome sequenced so far and more than twice the size of the Australian (Neoceratodus forsteri)1 and African2 lungfishes owing to enlarged intergenic regions and introns with high repeat content (about 90%). All lungfish genomes continue to expand as some transposable elements (TEs) are still active today. In particular, Lepidosiren's genome grew extremely fast during the past 100 million years (Myr), adding the equivalent of one human genome every 10 Myr. This massive genome expansion seems to be related to a reduction of PIWI-interacting RNAs and C2H2 zinc-finger and Krüppel-associated box (KRAB)-domain protein genes that suppress TE expansions. Although TE abundance facilitates chromosomal rearrangements, lungfish chromosomes still conservatively reflect the ur-tetrapod karyotype. Neoceratodus' limb-like fins still resemble those of their extinct relatives and remained phenotypically static for about 100 Myr. We show that the secondary loss of limb-like appendages in the Lepidosiren-Protopterus ancestor was probably due to loss of sonic hedgehog limb-specific enhancers.

The African coelacanth genome provides insights into tetrapod evolution.

The discovery of a living coelacanth specimen in 1938 was remarkable, as this lineage of lobe-finned fish was thought to have become extinct 70 million years ago. The modern coelacanth looks remarkably similar to many of its ancient relatives, and its evolutionary proximity to our own fish ancestors provides a glimpse of the fish that first walked on land. Here we report the genome sequence of the African coelacanth, Latimeria chalumnae. Through a phylogenomic analysis, we conclude that the lungfish, and not the coelacanth, is the closest living relative of tetrapods. Coelacanth protein-coding genes are significantly more slowly evolving than those of tetrapods, unlike other genomic features. Analyses of changes in genes and regulatory elements during the vertebrate adaptation to land highlight genes involved in immunity, nitrogen excretion and the development of fins, tail, ear, eye, brain and olfaction. Functional assays of enhancers involved in the fin-to-limb transition and in the emergence of extra-embryonic tissues show the importance of the coelacanth genome as a blueprint for understanding tetrapod evolution.

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.

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.

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.

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.

Evolutionary active transposable elements in the genome of the coelacanth.

The apparent morphological stasis in the lineage of the coelacanth, which has been called a "living fossil" by many, has been suggested to be causally related to a slow evolution of its genome, with strongly reduced activity of transposable elements (TEs). Analysis of the African coelacanth showed that at least 25% of its genome is constituted of transposable elements including retrotransposons, endogenous retroviruses and DNA transposons, with a strong predominance of non-Long Terminal Repeat (non-LTR) retrotransposons. The coelacanth genome has been shaped by four major general bursts of transposition during evolution, with major contributions of LINE1, LINE2, CR1, and Deu non-LTR retrotransposons. Many transposable elements are expressed in different tissues and might be active. The number of TE families in coelacanth, but also in lungfish, is lower than in teleost fish, but is higher than in chicken and human. This observation is in agreement with the hypothesis of a sequential elimination of many TE families in the sarcopterygian lineage during evolution. Taken together, our analysis indicates that the coelacanth contains more TE families than birds and mammals, and that these elements have been active during the evolution of the coelacanth lineage. Hence, at the level of transposable element activity, the coelacanth genome does not appear to evolve particularly slowly.

Transcriptional activity of transposable elements in coelacanth.

The morphological stasis of coelacanths has long suggested a slow evolutionary rate. General genomic stasis might also imply a decrease of transposable elements activity. To evaluate the potential activity of transposable elements (TEs) in "living fossil" species, transcriptomic data of Latimeria chalumnae and its Indonesian congener Latimeria menadoensis were compared through the RNA-sequencing mapping procedures in three different organs (liver, testis, and muscle). The analysis of coelacanth transcriptomes highlights a significant percentage of transcribed TEs in both species. Major contributors are LINE retrotransposons, especially from the CR1 family. Furthermore, some particular elements such as a LF-SINE and a LINE2 sequences seem to be more expressed than other elements. The amount of TEs expressed in testis suggests possible transposition burst in incoming generations. Moreover, significant amount of TEs in liver and muscle transcriptomes were also observed. Analyses of elements displaying marked organ-specific expression gave us the opportunity to highlight exaptation cases, that is, the recruitment of TEs as new cellular genes, but also to identify a new Latimeria-specific family of Short Interspersed Nuclear Elements called CoeG-SINEs. Overall, transcriptome results do not seem to be in line with a slow-evolving genome with poor TE activity.

A multicopy Y-chromosomal SGNH hydrolase gene expressed in the testis of the platyfish has been captured and mobilized by a Helitron transposon.

BACKGROUND: Teleost fish present a high diversity of sex determination systems, with possible frequent evolutionary turnover of sex chromosomes and sex-determining genes. In order to identify genes involved in male sex determination and differentiation in the platyfish Xiphophorus maculatus, bacterial artificial chromosome contigs from the sex-determining region differentiating the Y from the X chromosome have been assembled and analyzed. RESULTS: A novel three-copy gene called teximY (for testis-expressed in Xiphophorus maculatus on the Y) was identified on the Y but not on the X chromosome. A highly related sequence called texim1, probably at the origin of the Y-linked genes, as well as three more divergent texim genes were detected in (pseudo)autosomal regions of the platyfish genome. Texim genes, for which no functional data are available so far in any organism, encode predicted esterases/lipases with a SGNH hydrolase domain. Texim proteins are related to proteins from very different origins, including proteins encoded by animal CR1 retrotransposons, animal platelet-activating factor acetylhydrolases (PAFah) and bacterial hydrolases. Texim gene distribution is patchy in animals. Texim sequences were detected in several fish species including killifish, medaka, pufferfish, sea bass, cod and gar, but not in zebrafish. Texim-like genes are also present in Oikopleura (urochordate), Amphioxus (cephalochordate) and sea urchin (echinoderm) but absent from mammals and other tetrapods. Interestingly, texim genes are associated with a Helitron transposon in different fish species but not in urochordates, cephalochordates and echinoderms, suggesting capture and mobilization of an ancestral texim gene in the bony fish lineage. RT-qPCR analyses showed that Y-linked teximY genes are preferentially expressed in testis, with expression at late stages of spermatogenesis (late spermatids and spermatozeugmata). CONCLUSIONS: These observations suggest either that TeximY proteins play a role in Helitron transposition in the male germ line in fish, or that texim genes are spermatogenesis genes mobilized and spread by transposable elements in fish genomes.

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.

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.

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.

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.

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.