Observational study of the heterogeneous global meteotsunami generated after the Hunga Tonga-Hunga Ha'apai Volcano eruption.

The Hunga Tonga-Hunga Ha'apai volcano eruption of January 15th 2022 generated a global atmospheric and oceanic response that was recorded by an unprecedented amount of sensors. The eruption caused an atmospheric perturbation that travelled as a Lamb wave surrounding the Earth at least 3 times, and was recorded by hundreds of barographs worldwide. The atmospheric wave showed complex patterns of amplitude and spectral energy content, although most of the energy was concentrated in the band (2-120 min). Simultaneously to each passage of the atmospheric wave and after, significant Sea Level Oscillations (SLOs) in the tsunami frequency band were recorded by tide gauges located all around the globe, in what it can be referred to as a global meteotsunami. The amplitude and dominant frequency of the recorded SLOs showed a high spatial heterogeneity. Our point is that the geometry of continental shelves and harbours acted as tuners for the surface waves generated by the atmospheric disturbance at open sea, amplifying the signal at the eigenmodes of each shelf and harbour.

Atmospheric-ocean coupling drives prevailing and synchronic dispersal patterns of marine species with long pelagic durations.

Dispersal shapes population connectivity and plays a critical role in marine metacommunities. Prominent species for coastal socioecological systems, such as jellyfish and spiny lobsters, feature long pelagic dispersal phases (LPDPs), which have long been overlooked. Here, we use a cross-scale approach combining field surveys of these species with a high-resolution hydrodynamic model to decipher the underlying mechanisms of LPDP patterns in northwestern Mediterranean shores. We identified basin-scale prevailing dispersal routes and synchronic year-to-year patterns tightly linked to prominent circulation features typical of marginal seas and semienclosed basins, with an outstanding role of a retentive source area replenishing shores and potentially acting as a pelagic nursery area. We show how the atmospheric forcing of the ocean, a marked hydrological driver of the Mediterranean Sea, modulates dispersal routes and sources of LPDP at interannual scales. These findings represent a crucial advance in our understanding of the functioning of metapopulations of species with LPDP in marginal seas and may contribute to the effective management of coastal ecosystem services in the face of climate change.

Environmentally-induced sex reversal in fish with chromosomal vs. polygenic sex determination.

Sex ratio depends on sex determination mechanisms and is a key demographic parameter determining population viability and resilience to natural and anthropogenic stressors. There is increasing evidence that the environment can alter sex ratio even in genetically sex-determined species (GSD), as elevated temperature can cause female-to-male sex reversal (neomales). Alarmingly, neomales are being discovered in natural populations of several fish, amphibian and reptile species worldwide. Understanding the basis of neomale development is important for conservation biology. Among GSD species, it is unknown whether those with chromosomal sex determination (CSD), the most common system, will better resist the influence of high temperature than those with polygenic sex determination (PSD). Here, we compared the effects of elevated temperature in two wild zebrafish strains, Nadia (NA) and Ekkwill (EKW), which have CSD with a ZZ/ZW system, against the AB laboratory strain, which has PSD. First, we uncovered novel sex genotypes and the results showed that, at control temperature, the masculinization rate roughly doubled with the addition of each Z chromosome, while some ZW and WW fish of the wild strains became neomales. Surprisingly, we found that at elevated temperatures WW fish were just as likely as ZW fish to become neomales and that all strains were equally susceptible to masculinization. These results demonstrate that the Z chromosome is not essential for male development and that the dose of W buffers masculinization at the control temperature but not at elevated temperature. Furthermore, at the elevated temperature the testes of neomales, but not of normal males, contained more spermatozoa than at the control temperature. Our results show in an unprecedented way that, in a global warming scenario, CSD species may not necessarily be better protected against the masculinizing effect of elevated temperature than PSD species, and reveal genotype-by-temperature interactions in male sex determination and spermatogenesis.

Coastal flooding and mean sea-level rise allowances in atoll island.

Atoll islands are among the places most vulnerable to climate change due to their low elevation above mean sea level. Even today, some of these islands suffer from severe flooding generated by wind-waves, that will be exacerbated with mean sea-level rise. Wave-induced flooding is a complex physical process that requires computationally-expensive numerical models to be reliably estimated, thus limiting its application to single island case studies. Here we present a new model-based parameterisation for wave setup and a set of numerical simulations for the wave-induced flooding in coral reef islands as a function of their morphology, the Manning friction coefficient, wave characteristics and projected mean sea level that can be used for rapid, broad scale (e.g. entire atoll island nations) flood risk assessments. We apply this new approach to the Maldives to compute the increase in wave hazard due to mean sea-level rise, as well as the change in island elevation or coastal protection required to keep wave-induced flooding constant. While future flooding in the Maldives is projected to increase drastically due to sea-level rise, we show that similar impacts in nearby islands can occur decades apart depending on the exposure to waves and the topobathymetry of each island. Such assessment can be useful to determine on which islands adaptation is most urgently needed.

RADSex: A computational workflow to study sex determination using restriction site-associated DNA sequencing data.

The study of sex determination and sex chromosome organization in nonmodel species has long been technically challenging, but new sequencing methodologies now enable precise and high-throughput identification of sex-specific genomic sequences. In particular, restriction site-associated DNA sequencing (RAD-Seq) is being extensively applied to explore sex determination systems in many plant and animal species. However, software specifically designed to search for and visualize sex-biased markers using RAD-Seq data is lacking. Here, we present RADSex, a computational analysis workflow designed to study the genetic basis of sex determination using RAD-Seq data. RADSex is simple to use, requires few computational resources, makes no prior assumptions about the type of sex-determination system or structure of the sex locus, and offers convenient visualization through a dedicated R package. To demonstrate the functionality of RADSex, we re-analysed a published data set of Japanese medaka, Oryzias latipes, where we uncovered a previously unknown Y chromosome polymorphism. We then used RADSex to analyse new RAD-Seq data sets from 15 fish species spanning multiple taxonomic orders. We identified the sex determination system and sex-specific markers in six of these species, five of which had no known sex-markers prior to this study. We show that RADSex greatly facilitates the study of sex determination systems in nonmodel species thanks to its speed of analyses, low resource usage, ease of application and visualization options. Furthermore, our analysis of new data sets from 15 species provides new insights on sex determination in fish.

The Developmental and Genetic Architecture of the Sexually Selected Male Ornament of Swordtails.

Sexual selection results in sex-specific characters like the conspicuously pigmented extension of the ventral tip of the caudal fin-the "sword"-in males of several species of Xiphophorus fishes. To uncover the genetic architecture underlying sword formation and to identify genes that are associated with its development, we characterized the sword transcriptional profile and combined it with genetic mapping approaches. Results showed that the male ornament of swordtails develops from a sexually non-dimorphic prepattern of transcription factors in the caudal fin. Among genes that constitute the exclusive sword transcriptome and are located in the genomic region associated with this trait we identify the potassium channel, Kcnh8, as a sword development gene. In addition to its neural function kcnh8 performs a known role in fin growth. These findings indicate that during evolution of swordtails a brain gene has been co-opted for an additional novel function in establishing a male ornament.

Chromonomer: A Tool Set for Repairing and Enhancing Assembled Genomes Through Integration of Genetic Maps and Conserved Synteny.

The pace of the sequencing and computational assembly of novel reference genomes is accelerating. Though DNA sequencing technologies and assembly software tools continue to improve, biological features of genomes such as repetitive sequence as well as molecular artifacts that often accompany sequencing library preparation can lead to fragmented or chimeric assemblies. If left uncorrected, defects like these trammel progress on understanding genome structure and function, or worse, positively mislead this research. Fortunately, integration of additional, independent streams of information, such as a marker-dense genetic map and conserved orthologous gene order from related taxa, can be used to scaffold together unlinked, disordered fragments and to restructure a reference genome where it is incorrectly joined. We present a tool set for automating these processes, one that additionally tracks any changes to the assembly and to the genetic map, and which allows the user to scrutinize these changes with the help of web-based, graphical visualizations. Chromonomer takes a user-defined reference genome, a map of genetic markers, and, optionally, conserved synteny information to construct an improved reference genome of chromosome models: a "chromonome". We demonstrate Chromonomer's performance on genome assemblies and genetic maps that have disparate characteristics and levels of quality.

3D hotspots of marine litter in the Mediterranean: A modeling study.

The 3D dispersion of marine litter (ML) over the Mediterranean basin has been simulated using the velocity fields from a high resolution circulation model as base to run a 3D Lagrangian model. Three simulations have been performed to mimic the evolution of ML with density lower, similar, or higher than seawater. In all cases a realistic distribution of ML sources was used. Our results show that the accumulation/dispersion areas of the floating and buoyancy neutral particles are practically the same, although the latter are distributed in the water column, 80% of them found in the photic layer (average depth of 35m). Regarding to the densest particles, they rapidly sink and reach the seafloor close to their source. The regions of higher temporal variability mostly coincide with the ML accumulation regions. Weak seasonal variability occurs at a sub-basin scale as a result of the particles redistribution induced by the seasonal variability of the current field.

Multiple independent chromosomal fusions accompanied the radiation of the Antarctic teleost genus Trematomus (Notothenioidei:Nototheniidae).

BACKGROUND: Chromosomal rearrangements are thought to be an important driving force underlying lineage diversification, but their link to speciation continues to be debated. Antarctic teleost fish of the family Nototheniidae (Notothenioidei) diversified in a changing environmental context, which led to ecological, morphological, and genetic differentiation among populations. In addition, extensive chromosomal repatterning accompanied species divergence in several clades. The most striking karyotypic changes involved the recent species radiation (about 10 My) of the genus Trematomus, with chromosomal pair numbers ranging between 29 and 12. These dramatic reductions in chromosome number resulted mostly from large-scale chromosome fusions. Multiple centric and/or tandem fusions have been hypothesized in at least seven of the twelve recognized Trematomus species. To reconstruct their evolutionary history, we employed comparative cytogenomics (BAC-FISH and chromosome painting) to reveal patterns of interspecific chromosomal orthologies across several notothenioid clades. RESULTS: We defined orthologous chromosomal segments of reference, termed Structural Units (SUs). SUs were identified in a total of 18 notothenioid species. We demonstrated for the first time that SUs were strongly conserved across every specimen examined, with chromosomal syntenies highlighting a paucity of intrachromosomal macro-rearrangements. Multiple independent fusions of these SUs were inferred in the Trematomus species, in contrast to the shared SU fusions in species of the sister lineage Notothenia. CONCLUSIONS: The SU segments were defined units of chromosomal rearrangement in the entire family Nototheiidae, which diverged from the other notothenioid families 20 My ago. Some of the identified chromosomal syntenies within the SUs were even conserved in their closest relatives, the family Eleginopsidae. Comparing the timing of acquisition of the fusions in the closely related genera Notothenia and Trematomus of the nototheniid species family, we conclude that they exhibit distinct chromosomal evolutionary histories, which may be relevant to different speciation scenarios.

Tracking a mass mortality outbreak of pen shell Pinna nobilis populations: A collaborative effort of scientists and citizens.

A mass mortality event is devastating the populations of the endemic bivalve Pinna nobilis in the Mediterranean Sea from early autumn 2016. A newly described Haplosporidian endoparasite (Haplosporidium pinnae) is the most probable cause of this ecological catastrophe placing one of the largest bivalves of the world on the brink of extinction. As a pivotal step towards Pinna nobilis conservation, this contribution combines scientists and citizens' data to address the fast- and vast-dispersion and prevalence outbreaks of the pathogen. Therefore, the potential role of currents on parasite expansion was addressed by means of drift simulations of virtual particles in a high-resolution regional currents model. A generalized additive model was implemented to test if environmental factors could modulate the infection of Pinna nobilis populations. The results strongly suggest that the parasite has probably dispersed regionally by surface currents, and that the disease expression seems to be closely related to temperatures above 13.5 °C and to a salinity range between 36.5-39.7 psu. The most likely spread of the disease along the Mediterranean basin associated with scattered survival spots and very few survivors (potentially resistant individuals), point to a challenging scenario for conservation of the emblematic Pinna nobilis, which will require fast and strategic management measures and should make use of the essential role citizen science projects can play.

Epigenetic factors Dnmt1 and Uhrf1 coordinate intestinal development.

Intestinal tract development is a coordinated process involving signaling among the progenitors and developing cells from all three germ layers. Development of endoderm-derived intestinal epithelium has been shown to depend on epigenetic modifications, but whether that is also the case for intestinal tract cell types from other germ layers remains unclear. We found that functional loss of a DNA methylation machinery component, ubiquitin-like protein containing PHD and RING finger domains 1 (uhrf1), leads to reduced numbers of ectoderm-derived enteric neurons and severe disruption of mesoderm-derived intestinal smooth muscle. Genetic chimeras revealed that Uhrf1 functions both cell-autonomously in enteric neuron precursors and cell-non-autonomously in surrounding intestinal cells, consistent with what is known about signaling interactions between these cell types that promote one another's development. Uhrf1 recruits the DNA methyltransferase Dnmt1 to unmethylated DNA during replication. Dnmt1 is also expressed in enteric neurons and smooth muscle progenitors. dnmt1 mutants have fewer enteric neurons and disrupted intestinal smooth muscle compared to wildtypes. Because dnmt1;uhrf1 double mutants have a similar phenotype to dnmt1 and uhrf1 single mutants, Dnmt1 and Uhrf1 must function together during enteric neuron and intestinal muscle development. This work shows that genes controlling epigenetic modifications are important to coordinate intestinal tract development, provides the first demonstration that these genes influence development of the ENS, and advances uhrf1 and dnmt1 as potential new Hirschsprung disease candidates.

Antarctic blackfin icefish genome reveals adaptations to extreme environments.

Icefishes (suborder Notothenioidei; family Channichthyidae) are the only vertebrates that lack functional haemoglobin genes and red blood cells. Here, we report a high-quality genome assembly and linkage map for the Antarctic blackfin icefish Chaenocephalus aceratus, highlighting evolved genomic features for its unique physiology. Phylogenomic analysis revealed that Antarctic fish of the teleost suborder Notothenioidei, including icefishes, diverged from the stickleback lineage about 77 million years ago and subsequently evolved cold-adapted phenotypes as the Southern Ocean cooled to sub-zero temperatures. Our results show that genes involved in protection from ice damage, including genes encoding antifreeze glycoprotein and zona pellucida proteins, are highly expanded in the icefish genome. Furthermore, genes that encode enzymes that help to control cellular redox state, including members of the sod3 and nqo1 gene families, are expanded, probably as evolutionary adaptations to the relatively high concentration of oxygen dissolved in cold Antarctic waters. In contrast, some crucial regulators of circadian homeostasis (cry and per genes) are absent from the icefish genome, suggesting compromised control of biological rhythms in the polar light environment. The availability of the icefish genome sequence will accelerate our understanding of adaptation to extreme Antarctic environments.

Comparison of Xiphophorus and human melanoma transcriptomes reveals conserved pathway interactions.

Comparative analysis of human and animal model melanomas can uncover conserved pathways and genetic changes that are relevant for the biology of cancer cells. Spontaneous melanoma in Xiphophorus interspecies backcross hybrid progeny may be informative in identifying genes and functional pathways that are similarly related to melanoma development in all vertebrates, including humans. To assess functional pathways involved in the Xiphophorus melanoma, we performed gene expression profiling of the melanomas produced in interspecies BC1 and successive backcross generations (i.e., BC5 ) of the cross: X. hellerii × [X. maculatus Jp 163 A × X. hellerii]. Using RNA-Seq, we identified genes that are transcriptionally co-expressed with the driver oncogene, xmrk. We determined functional pathways in the fish melanoma that are also present in human melanoma cohorts that may be related to dedifferentiation based on the expression levels of pigmentation genes. Shared pathways between human and Xiphophorus melanomas are related to inflammation, cell migration, cell proliferation, pigmentation, cancer development, and metastasis. Our results suggest xmrk co-expressed genes are associated with dedifferentiation and highlight these signaling pathways as playing important roles in melanomagenesis.

Cold Fusion: Massive Karyotype Evolution in the Antarctic Bullhead Notothen Notothenia coriiceps.

Half of all vertebrate species share a series of chromosome fusions that preceded the teleost genome duplication (TGD), but we do not understand the causative evolutionary mechanisms. The "Robertsonian-translocation hypothesis" suggests a regular fusion of each ancestral acro- or telocentric chromosome to just one other by centromere fusions, thus halving the karyotype. An alternative "genome-stirring hypothesis" posits haphazard and repeated fusions, inversions, and reciprocal and nonreciprocal translocations. To study large-scale karyotype reduction, we investigated the decrease of chromosome numbers in Antarctic notothenioid fish. Most notothenioids have 24 haploid chromosomes, but bullhead notothen (Notothenia coriiceps) has 11. To understand mechanisms, we made a RAD-tag meiotic map with ∼10,000 polymorphic markers. Comparative genomics aligned about a thousand orthologs of platyfish and stickleback genes along bullhead chromosomes. Results revealed that 9 of 11 bullhead chromosomes arose by fusion of just two ancestral chromosomes and two others by fusion of three ancestral chromosomes. All markers from each ancestral chromosome remained contiguous, implying no inversions across fusion borders. Karyotype comparisons support a history of: (1) Robertsonian fusions of 22 ancestral chromosomes in pairs to yield 11 fused plus two small unfused chromosomes, like N. angustata; (2) fusion of one of the remaining two ancestral chromosomes to a preexisting fused pair, giving 12 chromosomes like N. rossii; and (3) fusion of the remaining ancestral chromosome to another fused pair, giving 11 chromosomes in N. coriiceps These results raise the question of what selective forces promoted the systematic fusion of chromosomes in pairs and the suppression of pericentric inversions in this lineage, and provide a model for chromosome fusions in stem teleosts.

Molecular genetic analysis of the melanoma regulatory locus in Xiphophorus interspecies hybrids.

Development of spontaneous melanoma in Xiphophorus interspecies backcross hybrid progeny, (X. hellerii × [X. maculatus Jp 163 A × X. hellerii]) is due to Mendelian segregation of a oncogene (xmrk) and a molecularly uncharacterized locus, called R(Diff), on LG5. R(Diff) is thought to suppresses the activity of xmrk in healthy X. maculatus Jp 163 A parental species that rarely develop melanoma. To better understand the molecular genetics of R(Diff), we utilized RNA-Seq to study allele-specific gene expression of spontaneous melanoma tumors and corresponding normal skin samples derived from 15 first generation backcross (BC1 ) hybrids and 13 fifth generation (BC5 ) hybrids. Allele-specific expression was determined for all genes and assigned to parental allele inheritance for each backcross hybrid individual. Results showed that genes residing in a 5.81 Mbp region on LG5 were exclusively expressed from the X. hellerii alleles in tumor-bearing BC1 hybrids. This observation indicates this region is consistently homozygous for X. hellerii alleles in tumor bearing animals, and therefore defines this region to be the R(Diff) locus. The R(Diff) locus harbors 164 gene models and includes the previously characterized R(Diff) candidate, cdkn2x. Twenty-one genes in the R(Diff) region show differential expression in the tumor samples compared to normal skin tissue. These results further characterize the R(Diff) locus and suggest tumor suppression may require a multigenic region rather than a single gene variant. Differences in gene expression between tumor and normal skin tissue in this region may indicate interactions among several genes are required for backcross hybrid melanoma development.

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.

A Genetic Map for the Only Self-Fertilizing Vertebrate.

The mangrove killifish Kryptolebias marmoratus, and its close relative Kryptolebias hermaphroditus, are the only vertebrate species known to reproduce by self-fertilization due to functional ovotestis development. To improve our understanding of their genomes, we constructed a genetic map. First, a single F1 fish was made by artificial fertilization between K. marmoratus and K. hermaphroditus strains. F2 progeny were then obtained by self-fertilization of the F1 fish. We used RAD-seq to query genomic DNAs from the two parental strains, the F1 individual and 49 F2 progeny. Results identified 9904 polymorphic RAD-tags (DNA markers) that mapped to 24 linkage groups, corresponding to the haploid chromosome number of these species. The total length of the map was 1248 cM, indicating that about one recombination occurred for each of the 24 homologous chromosome pairs in each meiosis. Markers were not evenly distributed along the chromosomes: in all chromosomes, many markers (> 8% of the total markers for each chromosome) mapped to chromosome tips. Centromeres suppress recombination, and this uneven distribution is probably due to the species' acrocentric chromosomes. Mapped marker sequences were compared to genomic sequences of medaka and platyfish, the next most closely related species with sequenced genomes that are anchored to genetic maps. Results showed that each mangrove killifish chromosome corresponds to a single chromosome of both platyfish and medaka, suggesting strong conservation of chromosomes over 100 million years of evolution. Our genetic map provides a framework for the K. marmoratus/K. hermaphroditus genome sequence and an important resource for understanding the biology of hermaphroditism.

Wild sex in zebrafish: loss of the natural sex determinant in domesticated strains.

Sex determination can be robustly genetic, strongly environmental, or genetic subject to environmental perturbation. The genetic basis of sex determination is unknown for zebrafish (Danio rerio), a model for development and human health. We used RAD-tag population genomics to identify sex-linked polymorphisms. After verifying this "RAD-sex" method on medaka (Oryzias latipes), we studied two domesticated zebrafish strains (AB and TU), two natural laboratory strains (WIK and EKW), and two recent isolates from nature (NA and CB). All four natural strains had a single sex-linked region at the right tip of chromosome 4, enabling sex genotyping by PCR. Genotypes for the single nucleotide polymorphism (SNP) with the strongest statistical association to sex suggested that wild zebrafish have WZ/ZZ sex chromosomes. In natural strains, "male genotypes" became males and some "female genotypes" also became males, suggesting that the environment or genetic background can cause female-to-male sex reversal. Surprisingly, TU and AB lacked detectable sex-linked loci. Phylogenomics rooted on D. nigrofasciatus verified that all strains are monophyletic. Because AB and TU branched as a monophyletic clade, we could not rule out shared loss of the wild sex locus in a common ancestor despite their independent domestication. Mitochondrial DNA sequences showed that investigated strains represent only one of the three identified zebrafish haplogroups. Results suggest that zebrafish in nature possess a WZ/ZZ sex-determination mechanism with a major determinant lying near the right telomere of chromosome 4 that was modified during domestication. Strains providing the zebrafish reference genome lack key components of the natural sex-determination system but may have evolved variant sex-determining mechanisms during two decades in laboratory culture.