The chromosome-scale reference genome for the pinfish (Lagodon rhomboides) provides insights into their evolutionary and demographic history.

The pinfish (Lagodon rhomboides) is an ecologically, economically, and culturally relevant member of family Sparidae, playing crucial roles in marine food webs of the western Atlantic Ocean and Gulf of Mexico. Despite their high abundance and ecological importance, there is a scarcity of genomic resources for this species. We assembled and annotated a chromosome-scale genome for the pinfish, resulting in a highly contiguous 785 Mb assembly of 24 scaffolded chromosomes. The high-quality assembly contains 98.9% complete BUSCOs and shows strong synteny to other chromosome-scale genomes of fish in the family Sparidae, with a limited number of large-scale genomic rearrangements. Leveraging this new genomic resource, we found evidence of significant expansions of dietary gene families over the evolutionary history of the pinfish, which may be associated with an ontogenetic shift from carnivory to herbivory seen in this species. Estimates of historical patterns of population demography using this new reference genome identified several periods of population growth and contraction which were associated with ancient climatic shifts and sea level changes. This genome serves as a valuable reference for future studies of population genomics and differentiation and provides a much-needed genomic resource for this western Atlantic sparid.

Insight into the adaptive role of arachnid genome-wide duplication through chromosome-level genome assembly of the Western black widow spider.

Although spiders are one of the most diverse groups of arthropods, the genetic architecture of their evolutionary adaptations is largely unknown. Specifically, ancient genome-wide duplication occurring during arachnid evolution ~450 mya resulted in a vast assembly of gene families, yet the extent to which selection has shaped this variation is understudied. To aid in comparative genome sequence analyses, we provide a chromosome-level genome of the Western black widow spider (Latrodectus hesperus)-a focus due to its silk properties, venom applications, and as a model for urban adaptation. We used long-read and Hi-C sequencing data, combined with transcriptomes, to assemble 14 chromosomes in a 1.46 Gb genome, with 38,393 genes annotated, and a BUSCO score of 95.3%. Our analyses identified high repetitive gene content and heterozygosity, consistent with other spider genomes, which has led to challenges in genome characterization. Our comparative evolutionary analyses of eight genomes available for species within the Araneoidea group (orb weavers and their descendants) identified 1,827 single-copy orthologs. Of these, 155 exhibit significant positive selection primarily associated with developmental genes, and with traits linked to sensory perception. These results support the hypothesis that several traits unique to spiders emerged from the adaptive evolution of ohnologs-or retained ancestrally duplicated genes-from ancient genome-wide duplication. These comparative spider genome analyses can serve as a model to understand how positive selection continually shapes ancestral duplications in generating novel traits today within and between diverse taxonomic groups.

The genome and population genomics of allopolyploid Coffea arabica reveal the diversification history of modern coffee cultivars.

Coffea arabica, an allotetraploid hybrid of Coffea eugenioides and Coffea canephora, is the source of approximately 60% of coffee products worldwide, and its cultivated accessions have undergone several population bottlenecks. We present chromosome-level assemblies of a di-haploid C. arabica accession and modern representatives of its diploid progenitors, C. eugenioides and C. canephora. The three species exhibit largely conserved genome structures between diploid parents and descendant subgenomes, with no obvious global subgenome dominance. We find evidence for a founding polyploidy event 350,000-610,000 years ago, followed by several pre-domestication bottlenecks, resulting in narrow genetic variation. A split between wild accessions and cultivar progenitors occurred ~30.5 thousand years ago, followed by a period of migration between the two populations. Analysis of modern varieties, including lines historically introgressed with C. canephora, highlights their breeding histories and loci that may contribute to pathogen resistance, laying the groundwork for future genomics-based breeding of C. arabica.

Limited evidence for extensive genetic differentiation between X and Y chromosomes in Hybognathus amarus (Cypriniformes: Leuciscidae).

Sex determination systems and genetic sex differentiation across fishes are highly diverse but are unknown for most Cypriniformes, including Rio Grande silvery minnow (Hybognathus amarus). In this study, we aimed to detect and validate sex-linked markers to infer sex determination system and to demonstrate the utility of combining several methods for sex-linked marker detection in nonmodel organisms. To identify potential sex-linked markers, Nextera-tagmented reductively amplified DNA (nextRAD) libraries were generated from 66 females, 64 males, and 60 larvae of unknown sex. These data were combined with female and male de novo genomes from Nanopore long-read sequences. We identified five potential unique male nextRAD-tags and one potential unique male contig, suggesting an XY sex determination system. We also identified two single-nucleotide polymorphisms (SNPs) in the same contig with values of FST, allele frequencies, and heterozygosity conforming with expectations of an XY system. Through PCR we validated the marker containing the sex-linked SNPs and a single nextRAD-tag sex-associated marker but it was not male specific. Instead, more copies of this locus in the male genome were suggested by enhanced amplification in males. Results are consistent with an XY system with low differentiation between sex-determining regions. Further research is needed to confirm the level of differentiation between the sex chromosomes. Nonetheless, this study highlighted the power of combining reduced representation and whole-genome sequencing for identifying sex-linked markers, especially when reduced representation sequencing does not include extensive variation between sexes, either because such variation is not present or not captured.

Lacustrine speciation associated with chromosomal inversion in a lineage of riverine fishes.

Geographic isolation is the primary driver of speciation in many vertebrate lineages. This trend is exemplified by North American darters, a clade of freshwater fishes where nearly all sister species pairs are allopatric and separated by millions of years of divergence. One of the only exceptions is the Lake Waccamaw endemic Etheostoma perlongum and its riverine sister species Etheostoma maculaticeps, which have no physical barriers to gene flow. Here we show that lacustrine speciation of E. perlongum is characterized by morphological and ecological divergence likely facilitated by a large chromosomal inversion. While E. perlongum is phylogenetically nested within the geographically widespread E. maculaticeps, there is a sharp genetic and morphological break coinciding with the lake-river boundary in the Waccamaw River system. Despite recent divergence, an active hybrid zone, and ongoing gene flow, analyses using a de novo reference genome reveal a 9 Mb chromosomal inversion with elevated divergence between E. perlongum and E. maculaticeps. This region exhibits striking synteny with known inversion supergenes in two distantly related fish lineages, suggesting deep evolutionary convergence of genomic architecture. Our results illustrate that rapid, ecological speciation with gene flow is possible even in lineages where geographic isolation is the dominant mechanism of speciation.

Navigating the CoGe Online Software Suite for Polyploidy Research.

The CoGe software suite at genomevolution.org hosts a number of tools that facilitate genomic research on plant and animal whole-genome multiplication-polyploidy. SynMap permits analysis and visualization of two-way syntenic dotplot alignments of genomes, includes many options and data/graphics download possibilities, and even permits three-genome synteny maps and interactive views. FractBias is a tool that operates within SynMap that permits calculation and graphic display of genome fragments (such as chromosomes) of one species mapped to another, displaying both blockwise homology depths and the extent of syntenic gene (syntelog) loss following polyploidy events. SynMap macrosynteny results can segue into the microsynteny tool GEvo, which provides genome-browser-like views of homologous genome blocks. CoGe FeatView allows call-up of given gene features already stored in the CoGe resource, and CoGeBlast permits searches for additional features that can be analyzed or downloaded further. Links from these tools can be fed into SynFind, which can find syntenic blocks surrounding a feature across multiple specified genomes while also simultaneously providing overall genome-wide syntenic depth calculations that can be interpreted to reflect polyploidy levels. Here, we describe basic use of these tools on the CoGe software suite.

Concordant patterns of morphological, stable isotope and genetic variation in a recent ecological radiation (Salmonidae: Coregonus spp.).

Groups of sympatric taxa with low interspecific genetic differentiation, but considerable ecological differences, offer great opportunities to study the dynamics of divergence and speciation. This is the case of ciscoes (Coregonus spp.) in the Laurentian Great Lakes, which are characterized by a complex evolutionary history and are commonly described as having undergone an adaptive radiation. In this study, morphometrics, stable isotopes and transcriptome sequencing were used to study the relationships within the Coregonus artedi complex in western Lake Superior. We observed general concordance for morphological, ecological and genomic variation, but the last was more taxonomically informative as it showed less overlap among species in multivariate space. Low levels of genetic differentiation were observed between individuals morphologically identified as Coregonus hoyi and C. zenithicus, which could be evidence of incomplete lineage sorting or recent hybridization between the two groups. Transcriptome-based single nucleotide polymorphisms exhibited significant divergence for genes associated with vision, development, metabolism and immunity among species that occupy different habitats. This study highlights the importance of using an integrative approach when studying groups of taxa with a complex evolutionary history, as individual-level analyses of multiple independent data sets can provide a clearer picture of the patterns and processes associated with the origins of biodiversity.

Chromosome-Level Genome Assembly of Chinese Sucker (Myxocyprinus asiaticus) Reveals Strongly Conserved Synteny Following a Catostomid-Specific Whole-Genome Duplication.

Fishes of the family Catostomidae ("suckers"; Teleostei: Cypriniformes) are hypothesized to have undergone an allopolyploidy event approximately 60 Ma. However, genomic evidence has previously been unavailable to assess this hypothesis. We sequenced and assembled the first chromosome-level catostomid genome, Chinese sucker (Myxocyprinus asiaticus), and present clear evidence of a catostomid-specific whole-genome duplication (WGD) event ("Cat-4R"). Our results reveal remarkably strong, conserved synteny since this duplication event, as well as between Myxocyprinus and an unduplicated outgroup, zebrafish (Danio rerio). Gene content and repetitive elements are also approximately evenly distributed across homeologous chromosomes, suggesting that both subgenomes retain some function, with no obvious bias in gene fractionation or subgenome dominance. The Cat-4R duplication provides another independent example of genome evolution following WGD in animals, in this case at the extreme end of conserved genome architecture over at least 25.2 Myr since the duplication. The M. asiaticus genome is a useful resource for researchers interested in understanding genome evolution following WGD in animals.

Nanopore Amplicon Sequencing Reveals Molecular Convergence and Local Adaptation of Rhodopsin in Great Lakes Salmonids.

Local adaptation can drive diversification of closely related species across environmental gradients and promote convergence of distantly related taxa that experience similar conditions. We examined a potential case of adaptation to novel visual environments in a species flock (Great Lakes salmonids, genus Coregonus) using a new amplicon genotyping protocol on the Oxford Nanopore Flongle and MinION. We sequenced five visual opsin genes for individuals of Coregonus artedi, Coregonus hoyi, Coregonus kiyi, and Coregonus zenithicus. Comparisons revealed species-specific differences in a key spectral tuning amino acid in rhodopsin (Tyr261Phe substitution), suggesting local adaptation of C. kiyi to the blue-shifted depths of Lake Superior. Ancestral state reconstruction demonstrates that parallel evolution and "toggling" at this amino acid residue has occurred several times across the fish tree of life, resulting in identical changes to the visual systems of distantly related taxa across replicated environmental gradients. Our results suggest that ecological differences and local adaptation to distinct visual environments are strong drivers of both evolutionary parallelism and diversification.

FiCli, the Fish and Climate Change Database, informs climate adaptation and management for freshwater fishes.

Inland fishes provide important ecosystem services to communities worldwide and are especially vulnerable to the impacts of climate change. Fish respond to climate change in diverse and nuanced ways, which creates challenges for practitioners of fish conservation, climate change adaptation, and management. Although climate change is known to affect fish globally, a comprehensive online, public database of how climate change has impacted inland fishes worldwide and adaptation or management practices that may address these impacts does not exist. We conducted an extensive, systematic primary literature review to identify peer-reviewed journal publications describing projected and documented examples of climate change impacts on inland fishes. From this standardized Fish and Climate Change database, FiCli (pronounced fick-lee), researchers and managers can query fish families, species, response types, or geographic locations to obtain summary information on inland fish responses to climate change and recommended management actions. The FiCli database is updatable and provides access to comprehensive published information to inform inland fish conservation and adaptation planning in a changing climate.

Comparative transcriptomics of cyprinid minnows and carp in a common wild setting: a resource for ecological genomics in freshwater communities.

Comparative transcriptomics can now be conducted on organisms in natural settings, which has greatly enhanced understanding of genome-environment interactions. Here, we demonstrate the utility and potential pitfalls of comparative transcriptomics of wild organisms, with an example from three cyprinid fish species (Teleostei:Cypriniformes). We present extensively filtered and annotated transcriptome assemblies that provide a valuable resource for studies of genome evolution (e.g. polyploidy), ecological and morphological diversification, speciation, and shared and unique responses to environmental variation in cyprinid fishes. Our results and analyses address the following points: (i) 'essential developmental genes' are shown to be ubiquitously expressed in a diverse suite of tissues across later ontogenetic stages (i.e. juveniles and adults), making these genes are useful for assessing the quality of transcriptome assemblies, (ii) the influence of microbiomes and other exogenous DNA, (iii) potentially novel, species-specific genes, and (iv) genomic rearrangements (e.g. whole genome duplication). The data we present provide a resource for future comparative work in cypriniform fishes and other taxa across a variety of sub-disciplines, including stress response, morphological diversification, community ecology, ecotoxicology, and climate change.

Stable isotope evidence for trophic overlap of sympatric Mexican Lake Chapala silversides (Teleostei: Atherinopsidae: Chirostoma spp)

We explore the trophic role that a diverse sympatric group of fishes in the genus Chirostoma play in a large, shallow lake in central Mexico, Lake Chapala. We use δ13C and δ15N stable isotope - based food web analyses to explore how they relate to other components of the Lake Chapala ecosystem. We find five Chirostoma species in top trophic levels of the Chapala food web compared to other fishes, relying on a combination of zooplankton, fish and benthic resources as energy sources. Food web metric analyses showed generally overlapping trophic niches for members of Chirostoma, especially in terms of δ13C. However, C. jordani had lower mean δ15N isotopic values than C. promelas. As a group, "pescados blancos" (C. sphyraena and C. promelas) also had higher δ15N signatures than "charales" (C. consocium, C. jordani and C. labarcae) reflecting greater piscivory, but these differences were not strong for all food web metrics used. Trophic overlap among species of Chirostoma in Lake Chapala raises questions about the forces that might have led to a morphologically diverse but functionally similar and monophyletic group of species.

Exploramos el papel trófico de un diverso grupo de peces (género Chirostoma) que habita en simpatría en el Lago Chapala, México central. Utilizamos isótopos estables de δ13C y δ15N para explorar la relación que guardan éstos peces con otros componentes del ecosistema. Encontramos a Chirostoma en niveles altos de la red trófica de Chapala, dependiendo energéticamente de zooplankton, peces y recursos del bentos. Análisis de parámetros de la red trófica demostraron traslape de nicho trófico para cinco miembros de Chirostoma, especialmente con relación a δ13C. Sin embargo, C. jordani tuvo un menor valor promedio de δ15N que C. promelas. Al ser analizados como grupo, los "pescados blancos" (C. sphyraena, C. promelas) tuvieron niveles promedio de δ15N mayores que los "charales" (C. consocium, C. jordani, C. labarcae), indicando mayor ictiofagia, pero las diferencias no fueron significativas para todos los parámetros de red trófica utilizados. El traslape trófico entre las especies de Chirostoma en Chapala da origen a cuestionamientos acerca de las fuerzas que pueden haber intervenido en el surgimiento de un morfológicamente diverso pero funcionalmente similar grupo monofilético de peces.

Clock gene evolution: seasonal timing, phylogenetic signal, or functional constraint?

Genetic determinants of seasonal reproduction are not fully understood but may be important predictors of organism responses to climate change. We used a comparative approach to study the evolution of seasonal timing within a fish community in a natural common garden setting. We tested the hypothesis that allelic length variation in the PolyQ domain of a circadian rhythm gene, Clock1a, corresponded to interspecific differences in seasonal reproductive timing across 5 native and 1 introduced cyprinid fishes (n = 425 individuals) that co-occur in the Rio Grande, NM, USA. Most common allele lengths were longer in native species that initiated reproduction earlier (Spearman's r = -0.70, P = 0.23). Clock1a allele length exhibited strong phylogenetic signal and earlier spawners were evolutionarily derived. Aside from length variation in Clock1a, all other amino acids were identical across native species, suggesting functional constraint over evolutionary time. Interestingly, the endangered Rio Grande silvery minnow (Hybognathus amarus) exhibited less allelic variation in Clock1a and observed heterozygosity was 2- to 6-fold lower than the 5 other (nonimperiled) species. Reduced genetic variation in this functionally important gene may impede this species' capacity to respond to ongoing environmental change.

A general hypothesis-testing framework for stable isotope ratios in ecological studies.

We propose a framework for hypothesis-testing of stable isotope ratios in ecological studies. Statistical procedures are based on analysis of nested linear models and a residual permutation procedure (RPP) that is employed to evaluate probabilities associated with test statistics. We used simulated examples and a real data set to illustrate the utility and generality of the method. First, we developed a test for differences in centroid location and dispersion of delta13C and delta15N values within and among groups of isotopic data. Second, we evaluated magnitude and direction of change in centroid position (termed "path") of a pair of isotopic samples separated in space/time relative to paths of other paired sample sets. Third, we compared attributes of path trajectories (size, direction, and shape) over sample sets containing more than two samples to provide a quantitative description of how patterns of isotopic ratios change in response to spatial and temporal gradients. Examples are limited to the bivariate case (delta13C-delta15N biplots), but the statistical method can readily be applied to univariate and multivariate cases.