Gonadal transcriptomics elucidate patterns of adaptive evolution within marine rockfishes (Sebastes)

Synthesis and Scope of the Role of Postmating Prezygotic Isolation in Speciation

How barriers to gene flow arise and are maintained are key questions in evolutionary biology. Speciation research has mainly focused on barriers that occur either before mating or after zygote formation. In comparison, postmating prezygotic (PMPZ) isolation-a barrier that acts after gamete release but before zygote formation-is less frequently investigated but may hold a unique role in generating biodiversity. Here we discuss the distinctive features of PMPZ isolation, including the primary drivers and molecular mechanisms underpinning PMPZ isolation. We then present the first comprehensive survey of PMPZ isolation research, revealing that it is a widespread form of prezygotic isolation across eukaryotes. The survey also exposes obstacles in studying PMPZ isolation, in part attributable to the challenges involved in directly measuring PMPZ isolation and uncovering its causal mechanisms. Finally, we identify outstanding knowledge gaps and provide recommendations for improving future research on PMPZ isolation. This will allow us to better understand the nature of this often-neglected reproductive barrier and its contribution to speciation.

Speciation genomics and the role of depth in the divergence of rockfishes (Sebastes) revealed through Pool-seq analysis of enriched sequences

Speciation in the marine environment is challenged by the wide geographic distribution of many taxa and potential for high rates of gene flow through larval dispersal mechanisms. Depth has recently been proposed as a potential driver of ecological divergence in fishes, and yet it is unclear how adaptation along these gradients' shapes genomic divergence. The genus Sebastes contains numerous species pairs that are depth-segregated and can provide a better understanding of the mode and tempo of genomic diversification. Here, we present exome data on two species pairs of rockfishes that are depth-segregated and have different degrees of divergence: S. chlorostictus-S. rosenblatti and S. crocotulus-S. miniatus. We were able to reliably identify "islands of divergence" in the species pair with more recent divergence (S. chlorostictus-S. rosenblatti) and discovered a number of genes associated with neurosensory function, suggesting a role for this pathway in the early speciation process. We also reconstructed demographic histories of divergence and found the best supported model was isolation followed by asymmetric secondary contact for both species pairs. These results suggest past ecological/geographic isolation followed by asymmetric secondary contact of deep to shallow species. Our results provide another example of using rockfish as a model for studying speciation and support the role of depth as an important mechanism for diversification in the marine environment.

Comparative Transcriptomics Reveals Patterns of Adaptive Evolution Associated with Depth and Age Within Marine Rockfishes (Sebastes)

The genetic underpinnings that contribute to ecological adaptation and speciation are not completely understood, especially within marine ecosystems. These evolutionary processes can be elucidated by studying adaptive radiations, because they provide replicates of divergence within a given environment or time-frame. Marine rockfishes (genus Sebastes) are an adaptive radiation and unique model system for studying adaptive evolution in the marine realm. We investigated molecular evolution associated with ecological (depth) and life history (lifespan) divergence in 2 closely related clades of Sebastes. Brain transcriptomes were sequenced via RNA-Seq from 3 species within the subgenus Pteropodus and a pair of related congeners from the subgenus Sebastosomus in order to identify patterns of adaptive evolution. De novo assemblies from these transcriptomes were used to identify 3867 orthologous clusters, and genes subject to positive selection were identified based on all 5 species, depth, and lifespan. Within all our analyses, we identified hemoglobin subunit α to be under strong positive selection and is associated with the depth of occurrence. In our lifespan analysis we identified immune function genes under positive selection in association with maximum lifespan. This study provides insight on the molecular evolution of rockfishes and these candidate genes may provide a better understanding of how these subgenera radiated within the Northeast Pacific.

Population assignment and local adaptation along an isolation-by-distance gradient in Pacific cod (Gadus macrocephalus)

The discernment of populations as management units is a fundamental prerequisite for sustainable exploitation of species. A lack of clear stock boundaries complicates not only the identification of spatial management units, but also the assessment of mixed fisheries by population assignment and mixed stock analysis. Many marine species, such as Pacific cod, are characterized by isolation by distance, showing significant differentiation but no clear stock boundaries. Here, we used restriction-site-associated DNA (RAD) sequencing to investigate population structure and assess power to genetically assign Pacific cod to putative populations of origin. Samples were collected across the species range in the eastern Pacific Ocean, from the Salish Sea to the Aleutian Islands. A total of 6,425 putative biallelic single nucleotide polymorphisms were identified from 276 individuals. We found a strong isolation-by-distance signal along coastlines that mirrored previous microsatellite results and pronounced genetic differentiation between coastal samples and those from the inland waters of the Salish Sea, with no evidence for hybridization between these two populations. Individual assignment success based on two methods was high overall (≥84%) but decreased from south to north. Assignment to geographic location of origin also was successful, with average distance between capture and assignment location of 220 km. Outlier analyses identified more loci potentially under selection along the coast than between Salish Sea and coastal samples, suggesting more diverse adaptation to latitudinal environmental factors than inshore vs. offshore environments. Our results confirm previous observations of sharp genetic differentiation of the Salish Sea population and isolation by distance along the coast, but also highlight the feasibility of using modern genomic techniques to inform stock boundaries and fisheries management in a low FST marine species.

Complete Mitochondrial Genome Sequences of Five Rockfishes (Perciformes: Sebastes)

Rockfishes of the genus Sebastes rank among the longest-lived vertebrate animals. In order to facilitate comparative genomic research in animal longevity, the complete mitochondrial genome sequences are presented for Sebastes aleutianus, Sebastes minor, Sebastes nigrocinctus, Sebastes rubrivinctus, and Sebastes steindachneri.

Comparative transcriptome analysis and identification of candidate effectors in two related rust species (Gymnosporangium yamadae and Gymnosporangium asiaticum)

BACKGROUND

Rust fungi constitute the largest group of plant fungal pathogens. However, a paucity of data, including genomic sequences, transcriptome sequences, and associated molecular markers, hinders the development of inhibitory compounds and prevents their analysis from an evolutionary perspective. Gymnosporangium yamadae and G. asiaticum are two closely related rust fungal species, which are ecologically and economically important pathogens that cause apple rust and pear rust, respectively, proved to be devastating to orchards. In this study, we investigated the transcriptomes of these two Gymnosporangium species during the telial stage of their lifecycles. The aim of this study was to understand the evolutionary patterns of these two related fungi and to identify genes that developed by selection.

RESULTS

The transcriptomes of G. yamadae and G. asiaticum were generated from a mixture of RNA from three biological replicates of each species. We obtained 49,318 and 54,742 transcripts, with N50 values of 1957 and 1664, for G. yamadae and G. asiaticum, respectively. We also identified a repertoire of candidate effectors and other gene families associated with pathogenicity. A total of 4947 pairs of putative orthologues between the two species were identified. Estimation of the non-synonymous/synonymous substitution rate ratios for these orthologues identified 116 pairs with Ka/Ks values greater than1 that are under positive selection and 170 pairs with Ka/Ks values of 1 that are under neutral selection, whereas the remaining 4661 genes are subjected to purifying selection. We estimate that the divergence time between the two species is approximately 5.2 Mya.

CONCLUSION

This study constitutes a de novo assembly and comparative analysis between the transcriptomes of the two rust species G. yamadae and G. asiaticum. The results identified several orthologous genes, and many expressed genes were identified by annotation. Our analysis of Ka/Ks ratios identified orthologous genes subjected to positive or purifying selection. An evolutionary analysis of these two species provided a relatively precise divergence time. Overall, the information obtained in this study increases the genetic resources available for research on the genetic diversity of the Gymnosporangium genus.